No Arabic abstract
We reexamine the screening mechanism in $f(R)$ gravity using N-body simulations. By explicitly examining the relation between the extra scalar field $delta f_R$ and the gravitational potential $phi$ in the perturbed Universe, we find that the relation between these two fields plays an important role in understanding the screening mechanism. We show that the screening mechanism in $f(R)$ gravity depends mainly on the depth of the potential well, and find a useful condition for identifying unscreened halos in simulations. We also discuss the potential application of our results to real galaxy surveys.
We use large volume N-body simulations to predict the clustering of dark matter in redshift space in f(R) modified gravity cosmologies. This is the first time that the nonlinear matter and velocity fields have been resolved to such a high level of accuracy over a broad range of scales in this class of models. We find significant deviations from the clustering signal in standard gravity, with an enhanced boost in power on large scales and stronger damping on small scales in the f(R) models compared to GR at redshifts z<1. We measure the velocity divergence (P_theta theta) and matter (P_delta delta) power spectra and find a large deviation in the ratios sqrt{P_theta theta/P_delta delta} and P_delta theta/P_deltadelta, between the f(R) models and GR for 0.03<k/(h/Mpc)<0.5. In linear theory these ratios equal the growth rate of structure on large scales. Our results show that the simulated ratios agree with the growth rate for each cosmology (which is scale dependent in the case of modified gravity) only for extremely large scales, k<0.06h/Mpc at z=0. The velocity power spectrum is substantially different in the f(R) models compared to GR, suggesting that this observable is a sensitive probe of modified gravity. We demonstrate how to extract the matter and velocity power spectra from the 2D redshift space power spectrum, P(k,mu), and can recover the nonlinear matter power spectrum to within a few percent for k<0.1h/Mpc. However, the model fails to describe the shape of the 2D power spectrum demonstrating that an improved model is necessary in order to reconstruct the velocity power spectrum accurately. The same model can match the monopole moment to within 3% for GR and 10% for the f(R) cosmology at k<0.2 h/Mpc at z=1. Our results suggest that the extraction of the velocity power spectrum from future galaxy surveys is a promising method to constrain deviations from GR.
We study the late-time Integrated Sachs-Wolfe (ISW) effect in $f(R)$ gravity using N-body simulations. In the $f(R)$ model under study, the linear growth rate is larger than that in general relativity (GR). This slows down the decay of the cosmic potential and induces a smaller ISW effect on large scales. Therefore, the $dotPhi$ (time derivative of the potential) power spectrum at $k<0.1h$/Mpc is suppressed relative to that in GR. In the non-linear regime, relatively rapid structure formation in $f(R)$ gravity boosts the non-linear ISW effect relative to GR, and the $dotPhi$ power spectrum at $k>0.1h$/Mpc is increased (100$%$ greater on small scales at $z=0$). We explore the detectability of the ISW signal via stacking supercluster and supervoids. The differences in the corresponding ISW cold or hot spots are $sim 20%$ for structures of $sim 100$Mpc/$h$. Such differences are greater for smaller structures, but the amplitude of the signal is lower. The high amplitude of ISW signal detected by Granett et al. can not explained in the $f(R)$ model. We find relatively big differences between $f(R)$ and GR in the transverse bulk motion of matter, and discuss its detectability via the relative frequency shifts of photons from multiple lensed images.
Based on thermodynamics, we discuss the galactic clustering of expanding Universe by assuming the gravitational interaction through the modified Newtons potential given by $f(R)$ gravity. We compute the corrected $N$-particle partition function analytically. The corrected partition function leads to more exact equations of states of the system. By assuming that system follows quasi-equilibrium, we derive the exact distribution function which exhibits the $f(R)$ correction. Moreover, we evaluate the critical temperature and discuss the stability of the system. We observe the effects of correction of $f(R)$ gravity on the power law behavior of particle-particle correlation function also. In order to check feasibility of an $f(R)$ gravity approach to the clustering of galaxies, we compare our results with an observational galaxy cluster catalog.
Modifications of the equations of general relativity at large distances offer one possibility to explain the observed properties of our Universe without invoking a cosmological constant. Numerous proposals for such modified gravity cosmologies exist, but often their consequences for structure formation in the non-linear sector are not yet accurately known. In this work, we employ high-resolution numerical simulations of f(R)-gravity models coupled with a semi-analytic model (SAM) for galaxy formation to obtain detailed predictions for the evolution of galaxy properties. The f(R)-gravity models imply the existence of a `fifth-force, which is however locally suppressed, preserving the successes of general relativity on solar system scales. We show that dark matter haloes in f(R)-gravity models are characterized by a modified virial scaling with respect to the LCDM scenario, reflecting a higher dark matter velocity dispersion at a given mass. This effect is taken into account in the SAM by an appropriate modification of the mass--temperature relation. We find that the statistical properties predicted for galaxies (such as the stellar mass function and the cosmic star formation rate) in f(R)-gravity show generally only very small differences relative to LCDM, smaller than the dispersion between the results of different SAM models, which can be viewed as a measure of their systematic uncertainty. We also demonstrate that galaxy bias is not able to disentangle between f(R)-gravity and the standard cosmological scenario. However, f(R)-gravity imprints modifications in the linear growth rate of cosmic structures at large scale, which can be recovered from the statistical properties of large galaxy samples.
Using N-body simulations, we measure the power spectrum of the effective dark matter density field, which is defined through the modified Poisson equation in $f(R)$ cosmologies. We find that when compared to the conventional dark matter power spectrum, the effective power spectrum deviates more significantly from the $Lambda$CDM model. For models with $f_{R0}=-10^{-4}$, the deviation can exceed 150% while the deviation of the conventional matter power spectrum is less than 50%. Even for models with $f_{R0}=-10^{-6}$, for which the conventional matter power spectrum is very close to the $Lambda$CDM prediction, the effective power spectrum shows sizeable deviations. Our results indicate that traditional analyses based on the dark matter density field may seriously underestimate the impact of $f(R)$ gravity on galaxy clustering. We therefore suggest the use of the effective density field in such studies. In addition, based on our findings, we also discuss several possible methods of making use of the differences between the conventional and effective dark matter power spectra in $f(R)$ gravity to discriminate the theory from the $Lambda$CDM model.