We present the merging of the Particle-Mesh (PM) relativistic GEVOLUTION code with the TreePM GADGET-4 code, with the aim of studying general relativity effects in cosmology. Our code, called GRGADGET, is able to track the evolution of metric perturbations in the weak field limit by using GEVOLUTION's implementation of a relativistic PM in the Poisson gauge. To achieve this, starting from GEVOLUTION, we have written a C++ library called LIBGEVOLUTION, which allows a code to access and use the same abstractions and resources that GEVOLUTION uses for its PM-only N-body simulations. The code works under the assumption that particle interactions at short distances can be approximated as Newtonian, so that we can combine the forces computed with a Newtonian Tree with those computed with a relativistic PM. The result is a TreePM simulation code that represents metric perturbations at the scales where they are relevant while resolving non-linear structures. We validate our code by closely matching GADGET-4 forces, computed with the Tree switched off, with those computed with LIBGEVOLUTION in the Newtonian limit. With GRGADGET, we obtain a matter power spectrum that is compatible with Newtonian GADGET-4 at small scales and contains GR features at large scales that are consistent with results obtained with GEVOLUTION. We demonstrate that, due to the better resolution of the highly non-linear regime, the representation of the relativistic fields sampled on the mesh improves with respect to the PM-only simulations.

GrGadget: an N-body TreePM relativistic code for cosmological simulations

Eduardo Quintana-Miranda
;
Pierluigi Monaco;Luca Tornatore
2023-01-01

Abstract

We present the merging of the Particle-Mesh (PM) relativistic GEVOLUTION code with the TreePM GADGET-4 code, with the aim of studying general relativity effects in cosmology. Our code, called GRGADGET, is able to track the evolution of metric perturbations in the weak field limit by using GEVOLUTION's implementation of a relativistic PM in the Poisson gauge. To achieve this, starting from GEVOLUTION, we have written a C++ library called LIBGEVOLUTION, which allows a code to access and use the same abstractions and resources that GEVOLUTION uses for its PM-only N-body simulations. The code works under the assumption that particle interactions at short distances can be approximated as Newtonian, so that we can combine the forces computed with a Newtonian Tree with those computed with a relativistic PM. The result is a TreePM simulation code that represents metric perturbations at the scales where they are relevant while resolving non-linear structures. We validate our code by closely matching GADGET-4 forces, computed with the Tree switched off, with those computed with LIBGEVOLUTION in the Newtonian limit. With GRGADGET, we obtain a matter power spectrum that is compatible with Newtonian GADGET-4 at small scales and contains GR features at large scales that are consistent with results obtained with GEVOLUTION. We demonstrate that, due to the better resolution of the highly non-linear regime, the representation of the relativistic fields sampled on the mesh improves with respect to the PM-only simulations.
2023
21-apr-2023
Pubblicato
https://academic.oup.com/mnras/article/522/4/5238/7136172
File in questo prodotto:
File Dimensione Formato  
Quintana2023.pdf

accesso aperto

Tipologia: Documento in Versione Editoriale
Licenza: Copyright dell'editore
Dimensione 2.29 MB
Formato Adobe PDF
2.29 MB Adobe PDF Visualizza/Apri
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11368/3048178
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 1
  • ???jsp.display-item.citation.isi??? 1
social impact