A number of codes for general-relativistic simulations of cosmological structure formation have been developed in recent years. Here we demonstrate that a sample of these codes produce consistent results beyond the Newtonian regime. We simulate solutions to Einstein's equations dominated by gravitomagnetism -- a vector-type gravitational field that doesn't exist in Newtonian gravity and produces frame-dragging, the leading-order post-Newtonian effect. We calculate the coordinate-invariant effect on intersecting null geodesics by performing ray tracing in each independent code. With this observable quantity, we assess and compare each code's ability to compute relativistic effects.

We address the generation of initial conditions (ICs) for GRAMSES, a code for nonlinear general relativistic (GR) $N$-body cosmological simulations recently introduced in Ref. [1]. GRAMSES adopts a constant mean curvature slicing with a minimal distortion gauge, where the linear growth rate is scale-dependent, and the standard method for realising initial particle data is not straightforwardly applicable. A new method is introduced, in which the initial positions of particles are generated from the displacement field realised for a matter power spectrum as usual, but the velocity is calculated by finite-differencing the displacement fields around the initial redshift. In this way, all the information required for setting up the initial conditions is drawn from three consecutive input matter power spectra, and additional assumptions such as scale-independence of the linear growth factor and growth rate are not needed. We implement this method in a modified 2LPTic code, and demonstrate that in a Newtonian setting it can reproduce the velocity field given by the default 2LPTic code with subpercent accuracy. We also show that the matter and velocity power spectra of the initial particle data generated for GRAMSES simulations using this method agree very well with the linear-theory predictions in the particular gauge used by GRAMSES. Finally, we discuss corrections to the finite difference calculation of the velocity when radiation is present, as well as additional corrections implemented in GRAMSES to ensure consistency. This method can be applied in ICs generation for GR simulations in generic gauges, and simulations of cosmological models with scale-dependent linear growth rate.

We study the imprints of an effective dark energy fluid in the large scale structure of the universe through the observed angular power spectrum of galaxies in the relativistic regime. We adopt the phenomenological approach that introduces two parameters $\{Q,\eta\}$ at the level of linear perturbations and allow to take into account the modified clustering (or effective gravitational constant) and anisotropic stress appearing in models beyond $\Lambda$CDM. We characterize the effective dark energy fluid by an equation of state parameter $w=-0.95$ and various sound speed cases in the range $10^{-6}\leq c^2_s\leq 1$, thus covering K-essence and quintessence cosmologies. We calculate the angular power spectra of standard and relativistic effects for these scenarios under the $\{Q,\eta\}$ parametrization, and we compare these relative to a fiducial $\Lambda$CDM cosmology. We find that, overall, deviations relative to $\Lambda$CDM are stronger at low redshift since the behavior of the dark energy fluid can mimic the cosmological constant during matter domination era but departs during dark energy domination. In particular, at $z=0.1$ the matter density fluctuations are suppressed by up to $\sim3\%$ for the quintessence-like case, while redshift-space distortions and Doppler effect can be enhanced by $\sim15\%$ at large scales for the lowest sound speed scenario. On the other hand, at $z=2$ we find deviations of up to $\sim5\%$ in gravitational lensing, whereas the Integrated Sachs-Wolfe effect can deviate up to $\sim17\%$. Furthermore, when considering an imperfect dark energy fluid scenario, we find that all effects are insensitive to the presence of anisotropic stress at low redshift, and only the Integrated Sachs-Wolfe effect can detect this feature at $z=2$ and very large scales.

We present GRAMSES, a new pipeline for nonlinear cosmological $N$-body simulations in General Relativity (GR). This code adopts the Arnowitt-Deser-Misner (ADM) formalism of GR, with constant mean curvature and minimum distortion gauge fixings, which provides a fully nonlinear and background independent framework for relativistic cosmology. Employing a fully constrained formulation, the Einstein equations are reduced to a set of ten elliptical equations which are solved using multigrid relaxation with adaptive mesh refinements (AMR), and three hyperbolic equations for the evolution of tensor degrees of freedom. The current version of GRAMSES neglects the latter by using the conformal flatness approximation, which allows it to compute the two scalar and two vector degrees of freedom of the metric. In this paper we describe the methodology, implementation, code tests and first results for cosmological simulations in a $\Lambda$CDM universe, while the generation of initial conditions and physical results will be discussed elsewhere. Inheriting the efficient AMR and massive parallelisation infrastructure from the publicly-available $N$-body and hydrodynamic simulation code RAMSES, GRAMSES is ideal for studying the detailed behaviour of spacetime inside virialised cosmic structures and hence accurately quantifying the impact of backreaction effects on the cosmic expansion, as well as for investigating GR effects on cosmological observables using cosmic-volume simulations.