My research focus on cosmology and large-scale structure evolution and formation.
My research focus on cosmology and large-scale structure evolution and formation.
To obtain precise cosmological parameter constraints, we need to know the accurate prediction for the statistics in every cosmological model, and then run the MCMC process on the measured statistics from observations. For previous surveys, we use the theoretical preditions such as halo model. While for DESI, CSST and after, we need more accurate predicitons. Running simulations is the ultimate way. But due to the limited computational resources, we can only sample the whole (6D~8D) cosmological parameter space using O(100) simulations. We complete this sampling by Halton sequence, and ran 129 simulations, each with particle number about 30 billion. The emulator for various statistics including those non-Gaussian ones containing information beyond two point, is under design.
Chen Zhao, Yu Yu +, in preparation.
Given a nonlinear density map, one can find a distorted mesh on which the density is constant. This mesh provides an estimation of the displacement field.
One can obtain the density field in Lagrangian space, i.e. the divergence of the displacement, which is referred as the reconstructed density field. This reconstructed density field has much better cross-correlation with the initial condition and thus sharpened baryon acoustic oscillation peaks.
Nonlinear reconstruction [ArXiV], [ADS]
Halo nonlinear reconstruction [ArXiV], [ADS]
Nonlinear reconstruction of redshift space distortions [ArXiV], [ADS]
Biased Tracer Reconstruction with Halo Mass Information [ArXiV], [ADS]
Weak lensing signal can be measured by looking at a large number of galaxy shapes, since the galaxy shapes are slightly changed by the weak lensing field, so-called cosmic shear. Weak lensing also changes the galaxy flux/size, and (de)magnifies the sky portions, leading to a change of the galaxy number density, so-called cosmic magnification.
Extracting the lensing signal from the galaxy distribution is very challenge, due to the intrinsic clustering signal which is several orders of magnitude larger. In future observation, the dependence of the cosmic magnification on the galaxy luminosity function can be measured precisely. Taking use of this dependence, one can separate cosmic magnification signal from the intrinsic clustering.
Weak lensing reconstruction through cosmic magnification - II. Improved power spectrum determination and map-making [ArXiV], [ADS]
Weak-lensing Power Spectrum Reconstruction by Counting Galaxies. I. The ABS Method [ArXiV], [ADS]
Weak-lensing Power Spectrum Reconstruction by Counting Galaxies. II. Improving the ABS Method with the Shift Parameter [ArXiV], [ADS]
Weak Lensing Magnification Reconstruction with the Modified Internal Linear Combination Method [ArXiV], [ADS]
After a logarithmic transform, the matter density field is Gaussianized. This process changes the information hierarchy structure, and the amount of information in power spectrum is increased.
Weak lensing field shares the similar property. In the first work we proved that the weak lensing field could be effectively Gaussianized by a local transform. However, the severe shape noise (the intrinsic shape of galaxies) hampers the application of the local transform. In the second work, wiener filtering is proposed to overcome this problem and recover the ability of the local transform.
Based on the finding that weak lensing field can be effectively Gaussianized by a local transform, we proposed the inverse-Gaussianization method to fast generate large amount of lensing mocks. These mocks have correct power spectrum and reasonably good power spectrum covariance. We also find that the peak height counts and peak steepness counts are well produced. In the covariance of peak statistics, super sample variance is a key ingredient and we also extend the inverse-Gaussianization method to include the super sample variance.
Gaussianizing the non-Gaussian lensing convergence field: The performance of the Gaussianization [ArXiV], [ADS]
Gaussianizing the non-Gaussian lensing convergence field II. The applicability to noisy data [ArXiV], [ADS]
Fast generation of weak lensing maps by the inverse-Gaussianization method [ArXiV], [ADS]
Weak-lensing Peak Statistics in Mocks by the inverse-Gaussianization method [ArXiV], [ADS]
Due to the late-time nonlinear evolution of the Universe, the cosmic fields are non-Gaussian. The information is transferred into high-order statistics, which is complicated and hard to model and measure. However, the local transform is found to be effective to Gaussianize the fields and pull the information back to power spectrum. But the local transform is imperfect. We can find some residual non-Gaussianity in the small scale. The perfectness of the local transform methods rely on the Gaussian Copula Hypothesis. We revisit this Hypothesis 10 years after it was proposed first time. Using high resolution simulation, we find obvious deviation from Gaussian copula in its differential form.
The Copula of the Cosmological Matter Density Field is Non-Gaussian [ArXiV], [ADS]
Numerical investigation of non-Gaussianities in the phase and modulus of density Fourier modes [ArXiV]
Volume-weighted statistics of large-scale peculiar velocity is preferred by peculiar velocity cosmology, since it is free of the uncertainties of galaxy density bias entangled in observed number density-weighted statistics. However, measuring the volume-weighted velocity statistics from galaxy (halo/simulation particle) velocity data is challenging.
In this series of works, we attempt to obtain volume-weighted velocity power spectrum by isotropic Kriging, an interpolation method widely used in geoscience. Kriging provides the most likely velocity field based on the data around and the prior correlation function.
Kriging interpolating cosmic velocity field [ArXiV], [ADS]
Kriging interpolating cosmic velocity field - II: improvement from more delicate Kriging [ArXiV], [ADS]