We investigate void properties in $f(R)$ models using N-body simulations, focusing on their differences from General Relativity (GR) and their detectability. In the Hu-Sawicki $f(R)$ modified gravity (MG) models, the halo number density profiles of voids are not distinguishable from GR. In contrast, the same $f(R)$ voids are more empty of dark matter, and their profiles are steeper. This can in principle be observed by weak gravitational lensing of voids, for which the combination of a spectroscopic redshift and a lensing photometric redshift survey over the same sky is required. Neglecting the lensing shape noise, the $f(R)$ model parameter amplitudes $|f_{R0}|=10^{-5}$ and $10^{-4}$ may be distinguished from GR using the lensing tangential shear signal around voids by 4 and 8$sigma$ for a volume of 1~(Gpc/$h$)$^3$. The line-of-sight projection of large-scale structure is the main systematics that limits the significance of this signal for the near future wide angle and deep lensing surveys. For this reason, it is challenging to distinguish $|f_{R0}|=10^{-6}$ from GR. We expect that this can be overcome with larger volume. The halo void abundance being smaller and the steepening of dark matter void profiles in $f(R)$ models are unique features that can be combined to break the degeneracy between $|f_{R0}|$ and $sigma_8$.
We explore voids in dark matter and halo fields from simulations of $Lambda$CDM and Hu-Sawicki $f(R)$ models. In $f(R)$ gravity, dark matter void abundances are greater than that of general relativity (GR). However, when using haloes to identify voids, the differences of void abundances become much smaller, but can still be told apart, in principle, at the 2, 6 and 14 $sigma$ level for the $f(R)$ model parameter amplitudes of $|f_{R0}|=10^{-6}$, $10^{-5}$ and $10^{-4}$. In contrast, the abundance of large voids found using haloes in $f(R)$ gravity is lower than in GR. The more efficient halo formation in underdense regions makes $f(R)$ voids less empty of haloes. This counter intuitive result suggests that voids are not necessarily emptier in $f(R)$ if one looks at galaxies in voids. Indeed, the halo number density profiles of voids are not distinguishable from GR. However, the same $f(R)$ voids are more empty of dark matter. This can in principle be observed by weak gravitational lensing of voids, for which the combination of a spec-$z$ and a photo-$z$ survey over the same sky is necessary. For a volume of 1~(Gpc/$h$)$^3$, neglecting the lensing shape noise, $|f_{R0}|=10^{-5}$ and $10^{-4}$ may be distinguished from GR using the lensing tangential shear signal around voids by 4 and 8$sigma$. The line-of-sight projection of large-scale structure is the main systematics that limits the significance of this signal, limiting the constraining power for $|f_{R0}|=10^{-6}$. The halo void abundance being smaller and the steepening of dark matter void profiles in $f(R)$ models are unique features that can be combined to break the degeneracy between $|f_{R0}|$ and $sigma_8$. The outflow of mass from void centers and velocity dispersions are greater in $f(R)$. Model differences in velocity profiles imply potential powerful constraints of the model in phase space and in redshift space.
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.
We measure the average temperature decrement on the cosmic microwave background (CMB) produced by voids selected in the SDSS DR7 spectroscopic redshift galaxy catalog, spanning redshifts $0<z<0.44$. We find an imprint of amplitude between 2.6 and 2.9$mu K$ as viewed through a compensated top-hat filter scaled to the radius of each void; we assess the statistical significance of the imprint at ~2$sigma$. We make crucial use of $N$-body simulations to calibrate our analysis. As expected, we find that large voids produce cold spots on the CMB through the Integrated Sachs-Wolfe (ISW) effect. However, we also find that small voids in the halo density field produce hot spots, because they reside in contracting, larger-scale overdense regions. This is an important effect to consider when stacking CMB imprints from voids of different radius. We have found that the same filter radius that gives the largest ISW signal in simulations also yields close to the largest detected signal in the observations. However, although it is low in significance, our measured signal is much higher-amplitude than expected from ISW in the concordance $Lambda$CDM universe. The discrepancy is also at the ~2$sigma$ level. We have demonstrated that our result is robust against the varying of thresholds over a wide range.
