We present the galaxy power spectrum in general relativity. Using a novel approach to computing the power spectrum, we derive the galaxy power spectrum taking into account all the relativistic effects in observations. We first demonstrate that the presence of an infrared divergence in previous calculations of the galaxy power spectrum is inconsistent with the equivalence principle, and we show that this unphysical behavior on large scales vanishes if we consider all the relativistic effects consistently. As a consequence, the relativistic effects in the galaxy power spectrum yield no corrections to the measurement of the primordial non-Gaussianity. We present our numerical computation of the full galaxy power spectrum, with particular emphasis on the deviations from the standard redshift-space power spectrum on large scales. We conclude that, as relativistic effects significantly alter the galaxy power spectrum at $k\lesssim k_{eq}$, they need to be taken into account in the analysis of large-scale data.

We perform theoretical and numerical studies of the full relativistic two-point galaxy correlation function, considering the linear-order scalar and tensor perturbation contributions and the wide-angle effects. Using the gauge-invariant relativistic description of galaxy clustering and accounting for the contributions at the observer position, we demonstrate that the complete theoretical expression is devoid of any long-mode contributions from scalar or tensor perturbations and it lacks the infrared divergences in agreement with the equivalence principle. By showing that the gravitational potential contribution to the correlation function converges in the infrared, our study justifies an IR cut-off $(k_{\text{IR}} \leq H_0)$ in computing the gravitational potential contribution. Using the full gauge-invariant expression, we numerically compute the galaxy two-point correlation function and study the individual contributions in the conformal Newtonian gauge. We find that the terms at the observer position such as the coordinate lapses and the observer velocity (missing in the standard formalism) dominate over the other relativistic contributions in the conformal Newtonian gauge such as the source velocity, the gravitational potential, the integrated Sachs-Wolf effect, the Shapiro time-delay and the lensing convergence. Compared to the standard Newtonian theoretical predictions that consider only the density fluctuation and redshift-space distortions, the relativistic effects in galaxy clustering result in a few percent-level systematic errors beyond the scale of the baryonic acoustic oscillation. Our theoretical and numerical study provides a comprehensive understanding of the relativistic effects in the galaxy two-point correlation function, as it proves the validity of the theoretical prediction and accounts for effects that are often neglected in its numerical evaluation.

We present the second-order expression for the observed redshift, accounting for all the relativistic effects from the light propagation and from the frame change at the observer and the source positions. We derive the generic gauge-transformation law that any observable quantities should satisfy, and we verify our second-order expression for the observed redshift by explicitly checking its gauge transformation property. This is the first time an explicit verification is made for the second-order calculations of observable quantities. We present our results in popular gauge choices for easy use and discuss the origin of disagreements in previous calculations.

We present the correlation function of the luminosity distances in a flat $\Lambda$CDM universe. Decomposing the luminosity distance fluctuation into the velocity, the gravitational potential, and the lensing contributions, we study their individual contributions to the correlation function. The lensing contribution is important at large redshift ($z\gtrsim 0.5$) but only for small angular separation ($\theta \lesssim 3^\circ$), while the velocity contribution dominates over the other contributions at low redshift or at larger separation. However, the gravitational potential contribution is always subdominant at all scale, if the correct gauge-invariant expression is used. The correlation function of the luminosity distances depends significantly on the matter content, especially for the lensing contribution, thus providing a novel tool of estimating cosmological parameters.

Measurements of the luminosity distance have played a key role in discovering the late-time cosmic acceleration. However, when accounting for inhomogeneities in the Universe, its interpretation has been plagued with infrared divergences in its theoretical predictions, which are in some cases used to explain the cosmic acceleration without dark energy. The divergences in most calculations are artificially removed by imposing an infrared cut-off scale. For the first time, we show that a gauge-invariant calculation of the luminosity distance is devoid of such divergences and consistent with the equivalence principle, eliminating the need to impose a cut-off scale. We present proper numerical calculations of the luminosity distance using the gauge-invariant expression and demonstrate that the numerical results with an ad hoc cut-off scale in previous calculations have negligible systematic errors as long as the cut-off scale is larger than the horizon scale. We discuss the origin of infrared divergences and their cancellation in the luminosity distance.

We investigate the effects of homogeneous general dark energy on the non-linear matter perturbation in fully general relativistic context. The equation for the density contrast contains even at linear order new contributions which are non-zero for general dark energy. Taking into account the next-leading-order corrections, we derive the total power spectrum in real and redshift spaces. We find that the observable galaxy power spectrum deviates from the LambdaCDM spectrum, which is nearly identical to that in the Einstein-de Sitter universe, and the relative difference is about 10% on a scale of the baryon acoustic oscillations.

We study the non-linear corrections to the matter and velocity power spectra in the synchronous gauge (SG). We consider the perturbations up to third order in a zero-pressure fluid in flat cosmological background, which is relevant for the non-linear growth of cosmic structure. As a result, we point out that the SG is an inappropriate coordinate choice when handling the non-linear growth of the large-scale structure. Although the equations in the SG happen to coincide with those in the comoving gauge (CG) to linear order, they differ from second order. In particular, the second order hydrodynamic equations in the the SG are apparently in the Lagrangian form, whereas those in the CG are in the Eulerian form. Thus, the non-linear power spectra naively presented in the original SG show strange behavior quite different from the result of the Newtonian theory even on sub-horizon scales. The power spectra in the SG show regularized behaviors only after we introduce convective terms in the second order so that the equations in two gauges coincide to the second order.

Under the standard perturbation theory (SPT), we obtain the fully consistent third-order density fluctuation and kernels for the general dark energy models without using the Einstein-de Sitter (EdS) universe assumption for the first time. We also show that even though the temporal and spatial components of the SPT solutions can not be separable, one can find the exact solutions to any order in general dark energy models. With these exact solutions, we obtain the less than \% error correction of one-loop matter power spectrum compared to that obtained from the EdS assumption for $k = 0.1 {\rm h\, Mpc}^{-1}$ mode at $z = 0$ (1, 1.5). Thus, the EdS assumption works very well at this scale. However, if one considers the correction for $P_{13}$, the error is about 6 (9, 11) \% for the same mode at $z = 0$ (1, 1.5). One absorbs $P_{13}$ into the linear power spectrum in the renormalized perturbation theory (RPT) and thus one should use the exact solution instead of the approximation one. The error on the resummed propagator $N$ of RPT is about 14 (8, 6) \% at $z =0$ (1, 1.5) for $k = 0.4 {\rm h\, Mpc}^{-1}$. For $k = 1 {\rm h\, Mpc}^{-1}$, the error correction of the total matter power spectrum is about 3.6 (4.6, 4.5) \% at $z = 0$ (1, 1.5). Upcoming observation is required to archive the sub-percent accuracy to provide the strong constraint on the dark energy and this consistent solution is prerequisite for the model comparison.

We show that the relativistic effects are negligibly small in the non-linear density and velocity bispectra. Although the non-linearities of Einstein equation introduce additional non-linear terms to the Newtonian fluid equations, the corrections to the bispectrum only show up on super-horizon scales. We show this with the next-to-leading order non-linear bispectrum for a pressureless fluid in a flat Friedmann-Robertson-Walker background, by calculating the density and velocity fields up to fourth order. We work in the comoving gauge, where the dynamics is identical to the Newtonian up to second order. We also discuss the leading order matter bispectrum in various gauges, and show yet another relativistic effect near horizon scales that the matter bispectrum strongly depends on the gauge choice.