Source author record

Yan-Chuan Cai

Yan-Chuan Cai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

Cosmological forecasts with the clustering of weak lensing peaks

Maximising the information that can be extracted from weak lensing measurements is a key goal for upcoming surveys such as LSST and Euclid. This is typically achieved through statistics that are complementary to the cosmic shear two-point correlation function, the most well established of which is the weak lensing peak abundance. In this work, we study the clustering of weak lensing peaks, and present parameter constraint forecasts for an LSST-like survey. We use the cosmoslics $w$CDM simulations to measure the peak two-point correlation function for a range of cosmological parameters, and use the simulation data to train a Gaussian process regression emulator which is applied to generate likelihood contours and provide parameter constraint forecasts from mock observations. We investigate the dependence of the peak two-point correlation function on the peak height, and find that the clustering of low amplitude peaks is complementary to that of high amplitude peaks. Consequently, their combination gives significantly tighter constraints than the clustering of high peaks alone. The peak two-point correlation function is significantly more sensitive to the cosmological parameters $h$ and $w_0$ than the peak abundance, and when the probes are combined, constraints on $Ω_{\rm m}$, $S_8$, $h$ and $w_0$ improve by at least a factor of two, relative to the peak abundance alone. Finally, we compare the forecasts for weak lensing peaks and weak lensing voids, and show that the two are also complementary; both probes can offer better constraints on $S_8$ and $w_0$ than the shear correlation function by roughly a factor of two.

preprint2022arXiv

Looking for a twist: probing the cosmological gravitomagnetic effect via weak lensing-kSZ cross correlations

General relativity predicts that the rotational momentum flux of matter twists the spacetime via a vector gravitomagnetic (frame-dragging) field, which remains undetected in cosmology. This vector field induces an additional gravitational lensing effect; at the same time, the momentum field sources the kinetic Sunyaev-Zel'dovich (kSZ) effect. The common origin of these two effects allows us to probe the gravitomagnetic signal via their cross-correlations. In this paper, we explore the possibility of detecting the gravitomagnetic field in $Λ$CDM by cross-correlating the weak-lensing convergence field with the CMB temperature map, which is imprinted with the kSZ signal. This approach allows us to extract the gravitomagnetic effect because the cross correlation between the standard Newtonian contribution to the weak-lensing convergence field, $κ_Φ$, and the kSZ effect is expected to vanish. We study the cross correlations with a suite of large-volume Newtonian $N$-body simulations and a small-volume, high-resolution, general-relativistic counterpart. We show that insufficient simulation resolution can introduce significant spurious correlations between $κ_Φ$ and kSZ. From the high-resolution simulation, we find that the cumulative signal-to-noise ratio (SNR) of the kSZ-gravitomagnetic convergence field can reach almost 15 (30) at $\ell\simeq5000$ ($10^4$) for the lensing source redshift $z_s=0.83$, if only cosmic variance is considered. We make forecast for next-generation lensing surveys such as Euclid and LSST, and CMB experiments such as Simons Observatory and CMB-S4, and find that, for $z_s=1.4$, the cumulative SNR can exceed 5 (9) at $\ell\simeq5000$ ($10^4$), indicating that the cosmological gravitomagnetic effect can be detected, if several foreground contaminations can be removed.

preprint2022arXiv

Understanding the relation between thermal Sunyaev-Zeldovich decrement and halo mass using the SIMBA and TNG simulations

The relation between the integrated thermal Sunyaev-Zeldovich (tSZ) $y$-decrement versus halo mass ($Y$--$M$) can potentially constrain galaxy formation models, if theoretical and observational systematics can be properly assessed. We investigate the $Y$--$M$ relation in the SIMBA and IllustrisTNG-100 cosmological hydrodynamic simulations, quantifying the effects of feedback, line-of-sight projection, and beam convolution. We find that SIMBA's AGN jet feedback generates strong deviations from self-similar expectations for the $Y$--$M$ relation, especially at $M_{500}<10^{13}M_{\odot}$. In SIMBA, this is driven by suppressed in-halo $y$ contributions owing to lowered halo baryon fractions. IllustrisTNG results more closely resemble SIMBA without jets. Projections of line-of-sight structures weaken these model differences slightly, but they remain significant -- particularly at group and lower halo masses. In contrast, beam smearing at $\textit{Planck}$ resolution makes the models indistinguishable, and both models appear to agree well with $\textit{Planck}$ data down to the lowest masses probed. We show that the arcminute resolution expected from forthcoming facilities would retain the differences between model predictions, and thereby provide strong constraints on AGN feedback.

preprint2021arXiv

FORGE -- the f(R) gravity cosmic emulator project I: Introduction and matter power spectrum emulator

We present a large suite of cosmological simulations, the FORGE (F-of-R Gravity Emulator) simulation suite, which is designed to build accurate emulators for cosmological observables in galaxy clustering, weak gravitational lensing and galaxy clusters, for the $f(R)$ gravity model. The total of 200 simulations explore the cosmological parameter space around the Planck(2018) cosmology with a Latin hypercube, for 50 combinations of $\bar{f}_{R0}$, $Ω_m$, $σ_8$ and $h$ with all other parameters fixed. For each parameter combination, or node, we ran four independent simulations, one pair using $1024^3$ particles in $500 h^{-1} Mpc$ simulation boxes to cover small scales, and another pair using $512^3$ simulation particles in $1500 h^{-1} Mpc$ boxes for larger scales. Each pair of initial conditions are selected such that sample variance on large scales is minimised on average. In this work we present an accurate emulator for the matter power spectrum in $f(R)$ gravity trained on FORGE. We have verified, using the cross-validation technique, that the emulator accuracy is better than $2.5\%$ for the majority of nodes, particularly around the center of the explored parameter space, up to scales of $k = 10 h Mpc^{-1}$. We have also checked the power spectrum emulator against simulations which are not part of our training set and found excellent agreement. Due to its high accuracy on small scales, the FORGE matter power spectrum emulator is well suited for weak lensing analysis and can play a key tool in constraining $f(R)$ gravity using current and future observational data.

preprint2020arXiv

Baryon acoustic oscillations reconstruction using convolutional neural networks

We propose a new scheme to reconstruct the baryon acoustic oscillations (BAO) signal, which contains key cosmological information, based on deep convolutional neural networks (CNN). Trained with almost no fine-tuning, the network can recover large-scale modes accurately in the test set: the correlation coefficient between the true and reconstructed initial conditions reaches $90\%$ at $k\leq 0.2 h\mathrm{Mpc}^{-1}$, which can lead to significant improvements of the BAO signal-to-noise ratio down to $k\simeq0.4h\mathrm{Mpc}^{-1}$. Since this new scheme is based on the configuration-space density field in sub-boxes, it is local and less affected by survey boundaries than the standard reconstruction method, as our tests confirm. We find that the network trained in one cosmology is able to reconstruct BAO peaks in the others, i.e. recovering information lost to non-linearity independent of cosmology. The accuracy of recovered BAO peak positions is far less than that caused by the difference in the cosmology models for training and testing, suggesting that different models can be distinguished efficiently in our scheme. It is very promising that Our scheme provides a different new way to extract the cosmological information from the ongoing and future large galaxy surveys.

preprint2019arXiv

A gravitational lensing detection of filamentary structures connecting luminous red galaxies

We present a weak lensing detection of filamentary structures in the cosmic web, combining data from the Kilo-Degree Survey, the Red Cluster Sequence Lensing Survey and the Canada-France-Hawaii Telescope Lensing Survey. The line connecting luminous red galaxies with a separation of $3 - 5\, h^{-1}\text{Mpc}$ is chosen as a proxy for the location of filaments. We measure the average weak lensing shear around $\sim$11,000 candidate filaments selected in this way from the Sloan Digital Sky Survey. After nulling the shear induced by the dark matter haloes around each galaxy, we report a $3.4\,σ$ detection of an anisotropic shear signal from the matter that connects them. Adopting a filament density profile, motivated from $N$-body simulations, the average density at the centre of these filamentary structures is found to be $15 \pm 4$ times the critical density.

preprint2016arXiv

Redshift-space distortions around voids

We have derived estimators for the linear growth rate of density fluctuations using the cross-correlation function of voids and haloes in redshift space, both directly and in Fourier form. In linear theory, this cross-correlation contains only monopole and quadrupole terms. At scales greater than the void radius, linear theory is a good match to voids traced out by haloes in N-body simulations; small-scale random velocities are unimportant at these radii, only tending to cause small and often negligible elongation of the redshift-space cross-correlation function near its origin. By extracting the monopole and quadrupole from the cross-correlation function, we measure the linear growth rate without prior knowledge of the void profile or velocity dispersion. We recover the linear growth parameter $β$ to 9% precision from an effective volume of 3(Gpc/h)^3 using voids with radius greater than 25Mpc/h. Smaller voids are predominantly sub-voids, which may be more sensitive to the random velocity dispersion; they introduce noise and do not help to improve the measurement. Adding velocity dispersion as a free parameter allows us to use information at radii as small as half of the void radius. The precision on $β$ is reduced to approximately 5%. Contrary to the simple redshift-space distortion pattern in overdensities, voids show diverse shapes in redshift space, and can appear either elongated or flattened along the line of sight. This can be explained by the competing amplitudes of the local density contrast, plus the radial velocity profile and its gradient, with the latter two factors being determined by the cumulative density profile of voids. The distortion pattern is therefore determined solely by the void profile and is different for void-in-cloud and void-in-void. This diversity of redshift-space void morphology complicates measurements of the Alcock-Paczynski effect using voids.

preprint2016arXiv

Testing the spherical evolution of cosmic voids

We study the spherical evolution model for voids in $Λ$CDM, where the evolution of voids is governed by dark energy at an earlier time than that for the whole universe or in overdensities. We show that the presence of dark energy suppresses the growth of peculiar velocities, causing void shell-crossing to occur at progressively later epochs as $Ω_Λ$ increases. We apply the spherical model to evolve the initial conditions of N-body simulated voids and compare the resulting final void profiles. We find that the model is successful in tracking the evolution of voids with radii greater than $30 h^{-1} \rm Mpc$, implying that void profiles could be used to constrain dark energy. We find that the initial peculiar velocities of voids play a significant role in shaping their evolution. Excluding the peculiar velocity in the evolution model delays the time of shell crossing.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2016arXiv

The lensing and temperature imprints of voids on the Cosmic Microwave Background

We have searched for the signature of cosmic voids in the CMB, in both the Planck temperature and lensing-convergence maps; voids should give decrements in both. We use zobov voids from the DR12 SDSS CMASS galaxy sample. We base our analysis on N-body simulations, to avoid a posteriori bias. For the first time, we detect the signature of voids in CMB lensing: the significance is $3.2σ$, close to $Λ$CDM in both amplitude and projected density-profile shape. A temperature dip is also seen, at modest significance ($2.3σ$), with amplitude about 6 times the prediction. This temperature signal is induced mostly by voids with radius between 100 and 150 Mpc/h, while the lensing signal is mostly contributed by smaller voids -- as expected; lensing relates directly to density, while ISW depends on gravitational potential. The void abundance in observations and simulations agree, as well. We also repeated the analysis excluding lower-significance voids: no lensing signal is detected, with an upper limit of about twice the $Λ$CDM prediction. But the mean temperature decrement now becomes non-zero at the $3.7σ$ level (similar to that found by Granett et al.), with amplitude about 20 times the prediction. However, the observed dependence of temperature on void size is in poor agreement with simulations, whereas the lensing results are consistent with $Λ$CDM theory. Thus, the overall tension between theory and observations does not favour non-standard theories of gravity, despite the hints of an enhanced amplitude for the ISW effect from voids.

preprint2016arXiv

The view from the boundary: a new void stacking method

We introduce a new method for stacking voids and deriving their profile that greatly increases the potential of voids as a tool for precision cosmology. Given that voids are highly non-spherical and have most of their mass at their edge, voids are better described relative to their boundary rather than relative to their centre, as in the conventional spherical stacking approach. The boundary profile is obtained by computing the distance of each volume element from the void boundary. Voids can then be stacked and their profiles computed as a function of this boundary distance. This approach enhances the weak lensing signal of voids, both shear and convergence, by a factor of two when compared to the spherical stacking method. It also results in steeper void density profiles that are characterised by a very slow rise inside the void and a pronounced density ridge at the void boundary. The resulting boundary density profile is self-similar when rescaled by the thickness of the density ridge, implying that the average rescaled profile is independent of void size. The boundary velocity profile is characterized by outflows in the inner regions whose amplitude scales with void size, and by a strong inflow into the filaments and walls delimiting the void. This new picture enables a straightforward discrimination between collapsing and expanding voids both for individual objects as well as for stacked samples.

preprint2014arXiv

A Possible Cold Imprint of Voids on the Microwave Background Radiation

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$μ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$σ$. 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 $Λ$CDM universe. The discrepancy is also at the ~2$σ$ level. We have demonstrated that our result is robust against the varying of thresholds over a wide range.

preprint2014arXiv

Testing Gravity using Cosmic Voids

We explore voids in dark matter and halo fields from simulations of $Λ$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 $σ$ 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$σ$. 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 $σ_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.

preprint2014arXiv

Testing Gravity using Void Profiles

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$σ$ 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 $σ_8$.

preprint2014arXiv

The Integrated Sachs-Wolfe effect in $f(R)$ gravity

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 $\dotΦ$ (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 $\dotΦ$ 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.

preprint2014arXiv

Void Dynamics

Cosmic voids are becoming key players in testing the physics of our Universe. Here we concentrate on the abundances and the dynamics of voids as these are among the best candidates to provide information on cosmological parameters. Cai, Padilla \& Li (2014) use the abundance of voids to tell apart Hu \& Sawicki $f(R)$ models from General Relativity. An interesting result is that even though, as expected, voids in the dark matter field are emptier in $f(R)$ gravity due to the fifth force expelling away from the void centres, this result is reversed when haloes are used to find voids. The abundance of voids in this case becomes even lower in $f(R)$ compared to GR for large voids. Still, the differences are significant and this provides a way to tell apart these models. The velocity field differences between $f(R)$ and GR, on the other hand, are the same for halo voids and for dark matter voids. Paz et al. (2013), concentrate on the velocity profiles around voids. First they show the necessity of four parameters to describe the density profiles around voids given two distinct void populations, voids-in-voids and voids-in-clouds. This profile is used to predict peculiar velocities around voids, and the combination of the latter with void density profiles allows the construction of model void-galaxy cross-correlation functions with redshift space distortions. When these models are tuned to fit the measured correlation functions for voids and galaxies in the Sloan Digital Sky Survey, small voids are found to be of the void-in-cloud type, whereas larger ones are consistent with being void-in-void. This is a novel result that is obtained directly from redshift space data around voids. These profiles can be used to remove systematics on void-galaxy Alcock-Pacinsky tests coming from redshift-space distortions.

preprint2014arXiv

Voids in Modified Gravity Reloaded: Eulerian Void Assignment

We revisit the excursion set approach to calculate void abundances in chameleon-type modified gravity theories, which was previously studied by Clampitt, Cai and Li (2013). We focus on properly accounting for the void-in-cloud effect, i.e., the growth of those voids sitting in over-dense regions may be restricted by the evolution of their surroundings. This effect may change the distribution function of voids hence affect predictions on the differences between modified gravity and GR. We show that the thin-shell approximation usually used to calculate the fifth force is qualitatively good but quantitatively inaccurate. Therefore, it is necessary to numerically solve the fifth force in both over-dense and under-dense regions. We then generalise the Eulerian void assignment method of Paranjape, Lam and Sheth (2012) to our modified gravity model. We implement this method in our Monte Carlo simulations and compare its results with the original Lagrangian methods. We find that the abundances of small voids are significantly reduced in both modified gravity and GR due to the restriction of environments. However, the change in void abundances for the range of void radii of interest for both models is similar. Therefore, the difference between models remains similar to the results from the Lagrangian method, especially if correlated steps of the random walks are used. As Clampitt, Cai and Li (2013), we find that the void abundance is much more sensitive to modified gravity than halo abundances. Our method can then be a faster alternative to N-body simulations for studying the qualitative behaviour of a broad class of theories. We also discuss the limitations and other practical issues associated with its applications.

preprint2013arXiv

Voids in Modified Gravity: Excursion Set Predictions

We investigate the behavior of the fifth force in voids in chameleon models using the spherical collapse method. Contrary to Newtonian gravity, we find the fifth force is repulsive in voids. The strength of the fifth force depends on the density inside and outside the void region as well as its radius. It can be many times larger than the Newtonian force and their ratio is in principle unbound. This is very different from the case in halos, where the fifth force is no more than 1/3 of gravity. The evolution of voids is governed by the Newtonian gravity, the effective dark energy force and the fifth force. While the first two forces are common in both LCDM and chameleon universes, the fifth force is unique to the latter. Driven by the outward-pointing fifth force, individual voids in chameleon models expand faster and grow larger than in a LCDM universe. The expansion velocity of the void shell can be 20% to 30% larger for voids of a few Mpc/h in radius, while their sizes can be larger by ~10%. These differences are smaller for larger voids of the same density. We compare void statistics using excursion set theory; for voids of the same size, their number density is larger in chameleon models. The fractional difference increases with void size. The chance of having voids of radius ~25 Mpc/h can be 2.5 times larger. This difference is about 10 times larger than that in the halo mass function. We find strong environmental dependence of void properties in chameleon models. The differences in size and expansion velocity with GR are both larger for small voids in high density regions. In general, the difference between chameleon models and LCDM in void properties (size, expansion velocity and distribution function) are larger than the corresponding quantities for halos. This suggests that voids might be better candidates than halos for testing gravity.

preprint2012arXiv

Combining weak lensing tomography and spectroscopic redshift surveys

Redshift space distortion (RSD) is a powerful way of measuring the growth of structure and testing General Relativity, but it is limited by cosmic variance and the degeneracy between galaxy bias b and the growth rate factor f. The cross-correlation of lensing shear with the galaxy density field can in principle measure b in a manner free from cosmic variance limits, breaking the f-b degeneracy and allowing inference of the matter power spectrum from the galaxy survey. We analyze the growth constraints from a realistic tomographic weak lensing photo-z survey combined with a spectroscopic galaxy redshift survey over the same sky area. For sky coverage f_sky=0.5, analysis of the transverse modes measures b to 2-3% accuracy per Δz=0.1 bin at z<1 when \sim10 galaxies per square arcmin are measured in the lensing survey and all halos with M>M_min=10^13 h^{-1}M_solar have spectra. For the gravitational growth parameter parameter γ(f=Ω_m^γ), combining the lensing information with RSD analysis of non-transverse modes yields accuracy σ(γ)\sim 0.01. Adding lensing information to the RSD survey improves σ(γ) by an amount equivalent to a 3 times (10 times) increase in RSD survey area when the spectroscopic survey extends down to halo mass 10^13.5 (10^14) h^{-1} M_solar. We also find that the σ(γ) of overlapping surveys is equivalent to that of surveys 1.5-2 times larger if they are separated on the sky. This gain is greatest when the spectroscopic mass threshold is 10^13 - 10^14 h^{-1} M_solar, similar to LRG surveys. The gain of overlapping surveys is reduced for very deep or very shallow spectroscopic surveys, but any practical surveys are more powerful when overlapped than when separated. The gain of overlapped surveys is larger in the case when the primordial power spectrum normalization is uncertain by >0.5%.

preprint2011arXiv

Cosmology without cosmic variance

We examine the improvements in constraints on the linear growth factor G and its derivative f=d ln G / dln a that are available from the combination of a large-scale galaxy redshift survey with a weak gravitational lensing survey of background sources. In the linear perturbation theory limit, the bias-modulation method of McDonald & Seljak allows one to distinguish the real-space galaxy clustering from the peculiar velocity signal in each Fourier mode. The ratio of lensing signal to galaxy clustering in transverse modes yields the bias factor b of each galaxy subset (as per Pen 2004), hence calibrating the conversion from galaxy real-space density to matter density in every mode. In combination these techniques permit measure of the growth rate f in each Fourier mode. This yields a measure of the growth rate free of sample variance, i.e. the uncertainty in f can be reduced without bound by increasing the number of redshifts within a finite volume. In practice, the gain from the absence of sample variance is bounded by the limited range of bias modulation among dark-matter halos. Nonetheless, the addition of background weak lensing data to a redshift survey increases information on G and f by an amount equivalent to a 10-fold increase in the volume of a standard redshift-space distortion measurement---if the lensing signal can be measured to sub-percent accuracy. This argues that a combined lensing and redshift survey over a common low-redshift volume is a more powerful test of general relativity than an isolated redshift survey over larger volume at high redshift. An example case is that a survey of ~10^6 redshifts over half the sky in the redshift range $z=0.5\pm 0.05$ can determine the growth exponent γfor the model $f=Ω_m^γ$ to an accuracy of $\pm 0.015$, using only modes with k<0.1h/Mpc, but only if a weak lensing survey is conducted in concert. [Abridged]

preprint2010arXiv

Full-sky map of the ISW and Rees-Sciama effect from Gpc simulations

We present a new method for constructing maps of the secondary temperature fluctuations imprinted on the cosmic microwave background (CMB) radiation by photons propagating through the evolving cosmic gravitational potential. Large cosmological N-body simulations are used to calculate the complete non-linear evolution of the peculiar gravitational potential. Tracing light rays back through the past lightcone of a chosen observer accurately captures the temperature perturbations generated by linear (the integrated Sachs-Wolfe or ISW effect) and non-linear (the Rees-Sciama or RS effect) evolution. These effects give rise to three kinds of non-linear features in the temperature maps. (a) In overdense regions, converging flows of matter induce cold spots of order 100 Mpc in extent which can dominate over the ISW effect at high redshift, and are surrounded by hot rings. (b) In underdense regions, the RS effect enhances ISW cold spots which can be surrounded by weak hot rings. (c) Transverse motions of large lumps of matter produce characteristic dipole features, consisting of adjacent hot and cold spots separated by a few tens of Megaparsecs. These non-linear features are not easily detectable; they modulate the ISW sky maps at about the 10 percent level. The RS effect causes the angular power spectrum to deviate from linear theory at l~50 and generates non-Gaussianity, skewing the one-point distribution function to negative values. Cold spots of similar angular size, but much smaller amplitude than the CMB cold spot reported by Cruz et al. are produced. Joint analysis of our maps and the corresponding galaxy distribution may enable techniques to be developed to detect these non-linear, non-Gaussian features. Our maps are available at http://astro.dur.ac.uk/~cai/ISW

preprint2009arXiv

Towards accurate modelling of the ISW effect, the non-linear contribution

In a universe with a cosmological constant, the large-scale gravitational potential varies in time and this is, in principle, observable. Using an N-body simulation of a $Λ$CDM universe, we show that linear theory is not sufficiently accurate to predict the power spectrum of the time derivative, $\dotΦ$, needed to compute the imprint of large-scale structure on the cosmic microwave background (CMB). The linear part of the $\dotΦ$ power spectrum (the integrated Sachs-Wolfe effect or ISW) drops quickly as the relative importance of $Ω_Λ$ diminishes at high redshift, while the non-linear part (the Rees-Sciama effect or RS) evolves more slowly with redshift. Therefore, the deviation of the total power spectrum from linear theory occurs at larger scales at higher redshifts. The deviation occurs at $k\sim 0.1 $ $h$ Mpc$^{-1}$ at $z=0$. The cross-correlation power spectrum of the density $δ$ with $\dotΦ$ behaves differently to the power spectrum of $\dotΦ$. Firstly, the deviation from linear theory occurs at smaller scales ($k\sim 1 $ $h$ Mpc$^{-1}$ at $z=0$). Secondly, the correlation becomes negative when the non-linear effect dominates. For the cross-correlation power spectrum of galaxy samples with the CMB, the non-linear effect becomes significant at $l\sim 500$ and rapidly makes the cross power spectrum negative. For high redshift samples, the cross-correlation is expected to be suppressed by $5-10%$ on arcminute scales. The RS effect makes a negligible contribution to the large-scale ISW cross-correlation measurement. However, on arc-minute scales it will contaminate the expected cross-correlation signal induced by the Sunyaev-Zel'dovich effect.

preprint2006arXiv

Note on Redshift Distortion in Fourier Space

We explore features of redshift distortion in Fourier analysis of N-body simulations. The phases of the Fourier modes of the dark matter density fluctuation are generally shifted by the peculiar motion along the line of sight, the induced phase shift is stochastic and has probability distribution function (PDF) symmetric to the peak at zero shift while the exact shape depends on the wave vector, except on very large scales where phases are invariant by linear perturbation theory. Analysis of the phase shifts motivates our phenomenological models for the bispectrum in redshift space. Comparison with simulations shows that our toy models are very successful in modeling bispectrum of equilateral and isosceles triangles at large scales. In the second part we compare the monopole of the power spectrum and bispectrum in the radial and plane-parallel distortion to test the plane-parallel approximation. We confirm the results of Scoccimarro (2000) that difference of power spectrum is at the level of 10%, in the reduced bispectrum such difference is as small as a few percents. However, on the plane perpendicular to the line of sight of k_z=0, the difference in power spectrum between the radial and plane-parallel approximation can be more than 10%, and even worse on very small scales. Such difference is prominent for bispectrum, especially for those configurations of tilted triangles. The non-Gaussian signals under radial distortion on small scales are systematically biased downside than that in plane-parallel approximation, while amplitudes of differences depend on the opening angle of the sample to the observer. The observation gives warning to the practice of using the power spectrum and bispectrum measured on the k_z=0 plane as estimation of the real space statistics.