No Arabic abstract
Future large-scale structure surveys will measure three-point statistics with high statistical significance. This will offer significant improvements on our understanding of gravity, provided we can model these statistics accurately. We assess the performance of several schemes for theoretical modelling of the matter bispectrum, including halo-model based approaches and fitting formulae. We compare the model predictions against N-body simulations, considering scales up to $k_{rm max} = 4 h/{rm Mpc}$, well into non-linear regime of structure formation. Focusing on the equilateral configuration, we conduct this analysis for three theories of gravity: general relativity, $f(R)$ gravity, and the DGP braneworld model. Additionally, we compute the lensing convergence bispectrum for these models. We find that all current modelling prescriptions in modified gravity, in particular for theories with scale-dependent linear growth, fail to attain the accuracy required by the precision of the Stage IV surveys such as emph{Euclid}. Among these models, we find that a halo-model corrected fitting formula achieves the best overall performance.
We study perturbation theory for large-scale structure in the most general scalar-tensor theories propagating a single scalar degree of freedom, which include Horndeski theories and beyond. We model the parameter space using the effective field theory of dark energy. For Horndeski theories, the gravitational field and fluid equations are invariant under a combination of time-dependent transformations of the coordinates and fields. This symmetry allows one to construct a physical adiabatic mode which fixes the perturbation-theory kernels in the squeezed limit and ensures that the well-known consistency relations for large-scale structure, originally derived in general relativity, hold in modified gravity as well. For theories beyond Horndeski, instead, one generally cannot construct such an adiabatic mode. Because of this, the perturbation-theory kernels are modified in the squeezed limit and the consistency relations for large-scale structure do not hold. We show, however, that the modification of the squeezed limit depends only on the linear theory. We investigate the observational consequences of this violation by computing the matter bispectrum. In the squeezed limit, the largest effect is expected when considering the cross-correlation between different tracers. Moreover, the individual contributions to the 1-loop matter power spectrum do not cancel in the infrared limit of the momentum integral, modifying the power spectrum on non-linear scales.
Using a perturbative approach we solve stellar structure equations for low-density (solar-type) stars whose interior is described with a polytropic equation of state in scenarios involving a subset of modified gravity theories. Rather than focusing on particular theories, we consider a model-independent approach in which deviations from General Relativity are effectively described by a single parameter $xi$. We find that for length scales below those set by stellar General Relativistic radii the modifications introduced by modified gravity can affect the computed values of masses and radii. As a consequence, the stellar luminosity is also affected. We discuss possible further implications for higher density stars and observability of the effects before described.
We demonstrate that baryonification algorithms, which displace particles in gravity-only simulations according to physically-motivated prescriptions, can simultaneously capture the impact of baryonic physics on the 2 and 3-point statistics of matter. Specifically, we show that our implementation of a baryonification algorithm jointly fits the changes induced by baryons on the power spectrum and equilateral bispectrum on scales up to k < 5 h/Mpc and redshifts z<2, as measured in six different cosmological hydrodynamical simulations. The accuracy of our fits are typically 1% for the power spectrum, and for the equilateral and squeezed bispectra, which somewhat degrades to 3% for simulations with extreme feedback prescriptions. Our results support the physical assumptions underlying baryonification approaches, and encourage their use in interpreting weak gravitational lensing and other cosmological observables.
We study the properties of dark matter haloes in a wide range of modified gravity models, namely, $f(R)$, DGP, and interacting dark energy models. We study the effects of modified gravity and dark energy on the internal properties of haloes, such as the spin and the structural parameters. We find that $f(R)$ gravity enhance the median value of the Bullock spin parameter, but could not detect such effects for DGP and coupled dark energy. $f(R)$ also yields a lower median sphericity and oblateness, while coupled dark energy has the opposite effect. However, these effects are very small. We then study the interaction rate of haloes in different gravity, and find that only strongly coupled dark energy models enhance the interaction rate. We then quantify the enhancement of the alignment of the spins of interacting halo pairs by modified gravity. Finally, we study the alignment of the major axes of haloes with the large-scale structures. The alignment of the spins of interacting pairs of haloes in DGP and coupled dark energy models show no discrepancy with GR, while $f(R)$ shows a weaker alignment. Strongly coupled dark energy shows a stronger alignment of the halo shape with the large-scale structures.
Model-independent constraints on modified gravity models hitherto exist mainly on linear scales. A recently developed formalism presented a consistent parameterisation that is valid on all scales. Using this approach, we perform model-independent modified gravity $N$-body simulations on all cosmological scales with a time-dependent $mu$. We present convergence tests of our simulations, and we examine how well existing fitting functions reproduce the non-linear matter power spectrum of the simulations. We find that although there is a significant variation in the accuracy of all of the fitting functions over the parameter space of our simulations, the ReACT framework delivers the most consistent performance for the matter power spectrum. We comment on how this might be improved to the level required for future surveys such as Euclid and the Vera Rubin Telescope (LSST). We also show how to compute weak-lensing observables consistently from the simulated matter power spectra in our approach, and show that ReACT also performs best when fitting the weak-lensing observables. This paves the way for a full model-independent test of modified gravity using all of the data from such upcoming surveys.