Source author record

Pengjie Zhang

Pengjie Zhang 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

48works
9topics
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

48 published item(s)

preprint2025arXiv

Introduction to the Chinese Space Station Survey Telescope (CSST)

The Chinese Space Station Survey Telescope (CSST) is an upcoming Stage-IV sky survey telescope, distinguished by its large field of view (FoV), high image quality, and multi-band observation capabilities. It can simultaneously conduct precise measurements of the Universe by performing multi-color photometric imaging and slitless spectroscopic surveys. The CSST is equipped with five scientific instruments, i.e. Multi-band Imaging and Slitless Spectroscopy Survey Camera (SC), Multi-Channel Imager (MCI), Integral Field Spectrograph (IFS), Cool Planet Imaging Coronagraph (CPI-C), and THz Spectrometer (TS). Using these instruments, CSST is expected to make significant contributions and discoveries across various astronomical fields, including cosmology, galaxies and active galactic nuclei (AGN), the Milky Way and nearby galaxies, stars, exoplanets, Solar System objects, astrometry, and transients and variable sources. This review aims to provide a comprehensive overview of the CSST instruments, observational capabilities, data products, and scientific potential.

preprint2022arXiv

About One-point Statistics of the Ratio of Two Fourier-transformed Cosmic Fields and an Application

The Fourier transformation is an effective and efficient operation of Gaussianization at the one-point level. Using a set of N-body simulation data, we verified that the one-point distribution functions of the dark matter momentum divergence and density fields closely follow complex Gaussian distributions. The one-point distribution function of the quotient of two complex Gaussian variables is introduced and studied. Statistical theories are then applied to model one-point statistics about the growth of individual Fourier mode of the dark matter density field, which can be obtained by the ratio of two Fourier transformed cosmic fields. Our simulation results proved that the models based on the Gaussian approximation are impressively accurate, and our analysis revealed many interesting aspects about the growth of dark matter's density fluctuation in Fourier space.

preprint2022arXiv

Cross-correlation of Planck CMB lensing with DESI galaxy groups

We measure the cross-correlation between galaxy groups constructed from DESI Legacy Imaging Survey DR8 and \emph{Planck} CMB lensing, over overlapping sky area of 16876 $\rm deg^2$. The detections are significant and consistent with the expected signal of the large-scale structure of the universe, over group samples of various redshift, mass, richness $N_{\rm g}$ and over various scale cuts. The overall S/N is 40 for a conservative sample with $N_{\rm g}\geq 5$, and increases to $50$ for the sample with $N_{\rm g}\geq 2$. Adopting the \emph{Planck} 2018 cosmology, we constrain the density bias of groups with $N_{\rm g}\geq 5$ as $b_{\rm g}=1.31\pm 0.10$, $2.22\pm 0.10$, $3.52\pm 0.20$ at $0.1<z\leq 0.33$, $0.33<z\leq 0.67$, $0.67<z\leq1$ respectively. The group catalog provides the estimation of group halo mass and therefore allows us to detect the dependence of bias on group mass with high significance. It also allows us to compare the measured bias with the theoretically predicted one using the estimated group mass. We find excellent agreement for the two high redshift bins. However, it is lower than the theory by $\sim 3σ$ for the lowest redshift bin. Another interesting finding is the significant impact of the thermal Sunyaev Zel'dovich (tSZ). It contaminates the galaxy group-CMB lensing cross-correlation at $\sim 30\%$ level, and must be deprojected first in CMB lensing reconstruction.

preprint2022arXiv

Forecasts on CMB lensing observations with AliCPT-1

AliCPT-1 is the first Chinese CMB experiment aiming for high precision measurement of Cosmic Microwave Background B-mode polarization. The telescope, currently under deployment in Tibet, will observe in two frequency bands centered at 90 and 150 GHz. We forecast the CMB lensing reconstruction, lensing-galaxy as well as lensing-CIB (Cosmic Infrared Background) cross correlation signal-to-noise ratio (SNR) for AliCPT-1. We consider two stages with different integrated observation time, namely "4 module*yr" (first stage) and "48 module*yr" (final stage). For lensing reconstruction, we use three different quadratic estimators, namely temperature-only, polarization-only and minimum-variance estimators, using curved sky geometry. We take into account the impact of inhomogeneous hit counts as well as of the mean-field bias due to incomplete sky coverage. In the first stage, our results show that the 150 GHz channel is able to measure the lensing signal at $15σ$ significance with the minimum-variance estimator. In the final stage, the measurement significance will increase to $31σ$. We also combine the two frequency data in the harmonic domain to optimize the SNR. Our result show that the coadding procedure can significantly reduce the reconstruction bias in the multiple range l>800. Thanks to the high quality of the polarization data in the final stage of AliCPT-1, the EB estimator will dominate the lensing reconstruction in this stage. We also estimate the SNR of cross-correlations between AliCPT-1 CMB lensing and other tracers of the large scale structure of the universe. For its cross-correlation with DESI galaxies/quasars, we report the cross-correlation SNR = 10-20 for the 4 redshift bins at 0.05<z<2.1. In the first stage, the total SNR is about $32$. In the final stage, the lensing-galaxy cross-correlation can reach SNR=52.

preprint2021arXiv

Numerical investigation of non-Gaussianities in the phase and modulus of density Fourier modes

We numerically investigate non-Gaussianities in the late-time cosmological density field in Fourier space. We explore various statistics, including the two-point and three-point probability distribution function (PDF) of phase and modulus, and two \& three-point correlation function of of phase and modulus. We detect significant non-Gaussianity for certain configurations. We compare the simulation results with the theoretical expansion series of \citet{2007ApJS..170....1M}. We find that the $\mathcal{O}(V^{-1/2})$ order term alone is sufficiently accurate to describe all the measured non-Gaussianities in not only the PDFs, but also the correlations. We also numerically find that the phase-modulus cross-correlation contributes $\sim 50\%$ to the bispectrum, further verifying the accuracy of the $\mathcal{O}(V^{-1/2})$ order prediction. This work demonstrates that non-Gaussianity of the cosmic density field is simpler in Fourier space, and may facilitate the data analysis in the era of precision cosmology.

preprint2021arXiv

The design of the Ali CMB Polarization Telescope receiver

Ali CMB Polarization Telescope (AliCPT-1) is the first CMB degree-scale polarimeter to be deployed on the Tibetan plateau at 5,250m above sea level. AliCPT-1 is a 90/150 GHz 72 cm aperture, two-lens refracting telescope cooled down to 4 K. Alumina lenses, 800mm in diameter, image the CMB in a 33.4° field of view on a 636mm wide focal plane. The modularized focal plane consists of dichroic polarization-sensitive Transition-Edge Sensors (TESes). Each module includes 1,704 optically active TESes fabricated on a 150mm diameter silicon wafer. Each TES array is read out with a microwave multiplexing readout system capable of a multiplexing factor up to 2,048. Such a large multiplexing factor has allowed the practical deployment of tens of thousands of detectors, enabling the design of a receiver that can operate up to 19 TES arrays for a total of 32,376 TESes. AliCPT-1 leverages the technological advancements in the detector design from multiple generations of previously successful feedhorn-coupled polarimeters, and in the instrument design from BICEP-3, but applied on a larger scale. The cryostat receiver is currently under integration and testing. During the first deployment year, the focal plane will be populated with up to 4 TES arrays. Further TES arrays will be deployed in the following years, fully populating the focal plane with 19 arrays on the fourth deployment year. Here we present the AliCPT-1 receiver design, and how the design has been optimized to meet the experimental requirements.

preprint2020arXiv

First detection of the GI-type of intrinsic alignments of galaxies using the self-calibration method in a photometric galaxy survey

Weak gravitational lensing is one of the most promising cosmological probes to constrain dark matter, dark energy, and the nature of gravity at cosmic scales. Intrinsic alignments (IAs) of galaxies have been recognized as one of the most serious systematic effects facing gravitational lensing. Such alignments must be isolated and removed to obtain a pure lensing signal. Furthermore, the alignments are related to the processes of galaxy formation, so their extracted signal can help in understanding such formation processes and improving their theoretical modeling. We report in this Letter the first detection of the gravitational shear--intrinsic shape (GI) correlation and the intrinsic shape--galaxy density (Ig) correlation using the self-calibration method in a photometric redshift survey. These direct measurements are made from the KiDS-450 photometric galaxy survey with a significance of 3.65$σ$ in the third bin for the Ig correlation, and 3.51$σ$ for the GI cross-correlation between the third and fourth bins. The self-calibration method uses the information available from photometric surveys without needing to specify an IA model and will play an important role in validating IA models and IA mitigation in future surveys such as the Rubin Observatory Legacy Survey of Space and Time, Euclid, and WFIRST.

preprint2020arXiv

Forecast for FAST: from Galaxies Survey to Intensity Mapping

The Five-Hundred-Meter Aperture Spherical Radio Telescope(FAST) is the largest single-dish radio telescope in the world. In this paper, we make forecast on the FAST HI large scale structure survey by mock observations. We consider a drift scan survey with the L-band 19 beam receiver, which may be commensal with the pulsar search and Galactic HI survey. We also consider surveys at lower frequency, either using the current single feed wide band receiver, or a future multi-beam phased array feed (PAF) in the UHF band. We estimate the number density of detected HI galaxies and the measurement error in positions, the precision of the surveys are evaluated using both Fisher matrix and simulated observations. The measurement error in the HI galaxy power spectrum is estimated, and we find that the error is relatively large even at moderate redshifts, as the number of positively detected galaxies drops drastically with increasing redshift. However, good cosmological measurement could be obtained with the intensity mapping technique where the large scale HI distribution is measured without resolving individual galaxies. The figure of merit (FoM) for the dark energy equation of state with different observation times are estimated, we find that with the existing L-band multi-beam receiver, a good measurement of low redshift large scale structure can be obtained, which complements the existing optical surveys. With a PAF in the UHF band, the constraint can be much stronger, reaching the level of a dark energy task force (DETF) stage IV experiment.

preprint2020arXiv

Neutrino effects on the morphology of cosmic large-scale structure

In this work, we propose a powerful probe of neutrino effects on the large-scale structure (LSS) of the Universe, i.e., Minkowski functionals (MFs). The morphology of LSS can be fully described by four MFs. This tool, with strong statistical power, is robust to various systematics and can comprehensively probe all orders of N-point statistics. By using a pair of high-resolution N-body simulations, for the first time, we comprehensively studied the subtle neutrino effects on the morphology of LSS. For an ideal LSS survey of volume $\sim1.73$ Gpc$^3$/$h^3$, neutrino signals are mainly detected from void regions with a significant level up to $\thicksim 10σ$ and $\thicksim 300σ$ for CDM and total matter density fields, respectively. This demonstrates its enormous potential for much improving the neutrino mass constraint in the data analysis of up-coming ambitious LSS surveys.

preprint2020arXiv

Separating the Intrinsic Alignment Signal and the Lensing Signal using Self-Calibration in Photo-z Surveys with KiDS450 and KV450 Data

To reach the full potential for the next generation of weak lensing surveys, it is necessary to mitigate the contamination of intrinsic alignments (IA) of galaxies in the observed cosmic shear signal. The self calibration (SC) of intrinsic alignments provides an independent method to measure the IA signal from the survey data and the photometric redshift information. It operates differently from the marginalization method based on IA modeling. In this work, we present the first application of SC to the KiDS450 data and the KV450 data, to split directly the intrinsic shape - galaxy density (Ig) correlation signal and the gravitational shear - galaxy density (Gg) correlation signal, using the information from photometric redshift (photo-z). We achieved a clear separation of the two signals and performed several validation tests. Our measured signals are found to be in general agreement with the KiDS450 cosmic shear best-fit cosmology, for both lensing and IA measurements. For KV450, we use partial (high-z) data, and our lensing measurements are also in good agreement with KV450 cosmic shear best-fit, however, our IA signal suggests a larger IA amplitude. We discussed the impact of photo-z quality on IA detection and several other potential systematic biases. Finally, we discuss the potential application of the information extracted for both the lensing signal and the IA signal in future surveys.

preprint2020arXiv

Strong lensing as a giant telescope to localize the host galaxy of gravitational wave event

Standard siren cosmology of gravitational wave (GW) merger events relies on the identification of host galaxies and their redshifts. But this can be highly challenging due to numerous candidates of galaxies in the GW localization area. We point out that the number of candidates can be reduced by orders of magnitude for strongly lensed GW events, due to extra observational constraints. For the next-generation GW detectors like Einstein Telescope (ET), we estimate that this number is usually significantly less than one, as long as the GW localization uncertainty is better than $\sim 10\, \rm deg^2$. This implies that the unique identification of the host galaxy of lensed GW event detected by ET and Cosmic Explorer (CE) is possible. This provides us a promising opportunity to measure the redshift of the GW event and facilitate the standard siren cosmology. We also discuss its potential applications in understanding the evolution process and environment of the GW event.

preprint2020arXiv

The copula of the cosmological matter density field is non-Gaussian

Non-Gaussianity of the cosmological matter density field can be largely reduced by a local Gaussianization transformation (and its approximations such as the logrithmic transformation). Such behavior can be recasted as the Gaussian copula hypothesis, and has been verified to very high accuracy at two-point level. On the other hand, statistically significant non-Gaussianities in the Gaussianized field have been detected in simulations. We point out that, this apparent inconsistency is caused by the very limited degrees of freedom in the copula function, which make it misleading as a diagnosis of residual non-Gaussianity in the Gaussianized field. Using the copula density, we highlight the departure from Gaussianity. We further quantify its impact in the predicted n-point correlation functions. We explore a remedy of the Gaussian copula hypothesis, which alleviates but not completely solves the above problems.

preprint2016arXiv

Fast generation of weak lensing maps by the inverse-Gaussianization method

To take full advantage of the unprecedented power of upcoming weak lensing surveys, understanding the noise, such as cosmic variance and geometry/mask effects, is as important as understanding the signal itself. Accurately quantifying the noise requires a large number of statistically independent mocks for a variety of cosmologies. This is impractical for weak lensing simulations, which are costly for simultaneous requirements of large box size (to cover a significant fraction of the past light cone) and high resolution (to robustly probe the small scale where most lensing signal resides). Therefore fast mock generation methods are desired and are under intensive investigation. We propose a new fast weak lensing map generation method, named the inverse-Gaussianization method, based on the finding that a lensing convergence field can be Gaussianized to excellent accuracy by a local transformation [Yu et al, Phys. Rev. D 84, 023523 (2011)]. Given a simulation, it enables us to produce as many as infinite statistically independent lensing maps as fast as producing the simulation initial conditions. The proposed method is tested against simulations for each tomography bin centered at lens redshift $z \sim 0.5$, 1, and 2, with various statistics. We find that the lensing maps generated by our method have reasonably accurate power spectra, bispectra, and power spectrum covariance matrix. Therefore, it will be useful for weak lensing surveys to generate realistic mocks. As an example of application, we measure the probability distribution function of the lensing power spectrum, from 16384 lensing maps produced by the inverse-Gaussianization method.

preprint2015arXiv

Determination of the large scale volume weighted halo velocity bias in simulations

A profound assumption in peculiar velocity cosmology is $b_v=1$ at sufficiently large scales, where $b_v$ is the volume weighted halo(galaxy) velocity bias with respect to the matter velocity field. However, this fundamental assumption has not been robustly verified in numerical simulations. Furthermore, it is challenged by structure formation theory (BBKS, 1986, ApJ; Desjacques and Sheth, 2010, PRD), which predicts the existence of velocity bias (at least for proto-halos) due to the fact that halos reside in special regions (local density peaks). The major obstacle to measure the volume weighted velocity from N-body simulations is an unphysical sampling artifact. It is entangled in the measured velocity statistics and becomes significant for sparse populations. With recently improved understanding of the sampling artifact (Zhang, Zheng and Jing, 2015, PRD; Zheng, Zhang and Jing, 2015, PRD), for the first time we are able to {\it appropriately correct this sampling artifact and then robustly measure the volume weighted halo velocity bias}. (1) We verify $b_v=1$ within $2\%$ model uncertainty at $k\lesssim 0.1h/$Mpc and $z=0$-$2$ for halos of mass $\sim 10^{12}$-$10^{13} h^{-1} M_\odot$, and, therefore, consolidates a foundation of the peculiar velocity cosmology. (2) We also find statistically significant signs of $b_v\neq 1$ at $k\gtrsim 0.1h/$Mpc. Unfortunately, whether this is real or caused by residual sampling artifact requires further investigation. Nevertheless, cosmology based on $k\gtrsim 0.1h/$Mpc velocity data shall be careful this potential velocity bias.

preprint2015arXiv

Diagnosing multiplicative error by lensing magnification of type Ia supernovae

Weak lensing causes spatially coherent fluctuations in flux of type Ia supernovae (SNe Ia). This lensing magnification allows for weak lensing measurement independent of cosmic shear. It is free of shape measurement errors associated with cosmic shear and can therefore be used to diagnose and calibrate multiplicative error. Although this lensing magnification is difficult to measure accurately in auto correlation, its cross correlation with cosmic shear and galaxy distribution in overlapping area can be measured to significantly higher accuracy. Therefore these cross correlations can put useful constraint on multiplicative error, and the obtained constraint is free of cosmic variance in weak lensing field. We present two methods implementing this idea and estimate their performances. We find that, with $\sim 1$ million SNe Ia that can be achieved by the proposed D2k survey with the LSST telescope (Zhan et al. 2008), multiplicative error of $\sim 0.5\%$ for source galaxies at $z_s\sim 1$ can be detected and larger multiplicative error can be corrected to the level of $0.5\%$. It is therefore a promising approach to control the multiplicative to the sub-percent level required for stage IV projects. The combination of the two methods even has the potential to diagnose and calibrate galaxy intrinsic alignment, which is another major systematic error in cosmic shear cosmology.

preprint2015arXiv

Kriging Interpolating Cosmic Velocity Field

[abridged] Volume-weighted statistics of large scale peculiar velocity is preferred by peculiar velocity cosmology, since it is free of uncertainties of galaxy density bias entangled in mass-weighted statistics. However, measuring the volume-weighted velocity statistics from galaxy (halo/simulation particle) velocity data is challenging. For the first time, we apply the Kriging interpolation to obtain the volume-weighted velocity field. Kriging is a minimum variance estimator. It predicts the most likely velocity for each place based on the velocity at other places. We test the performance of Kriging quantified by the E-mode velocity power spectrum from simulations. Dependences on the variogram prior used in Kriging, the number $n_k$ of the nearby particles to interpolate and the density $n_P$ of the observed sample are investigated. First, we find that Kriging induces $1\%$ and $3\%$ systematics at $k\sim 0.1h{\rm Mpc}^{-1}$ when $n_P\sim 6\times 10^{-2} ({\rm Mpc}/h)^{-3}$ and $n_P\sim 6\times 10^{-3} ({\rm Mpc}/h)^{-3}$, respectively. The deviation increases for decreasing $n_P$ and increasing $k$. When $n_P\lesssim 6\times 10^{-4} ({\rm Mpc}/h)^{-3}$, a smoothing effect dominates small scales, causing significant underestimation of the velocity power spectrum. Second, increasing $n_k$ helps to recover small scale power. However, for $n_P\lesssim 6\times 10^{-4} ({\rm Mpc}/h)^{-3}$ cases, the recovery is limited. Finally, Kriging is more sensitive to the variogram prior for lower sample density. The most straightforward application of Kriging on the cosmic velocity field does not show obvious advantages over the nearest-particle method (Zheng et al. 2013) and could not be directly applied to cosmology so far. However, whether potential improvements may be achieved by more delicate versions of Kriging is worth further investigation.

preprint2015arXiv

Sampling Artifact in Volume Weighted Velocity Measurement.--- I. Theoretical Modelling

Cosmology based on large scale peculiar velocity preferes volume weighted velocity statistics. However, measuring the volume weighted velocity statistics from inhomogeneously distributed galaxies (simulation particles/halos) suffer from an inevitable and significant sampling artifact. We study this sampling artifact in the velocity power spectrum measured by the nearest-particle (NP) velocity assignment method(Zheng et al. 2013, PRD). We derive the analytical expression of leading and higher order terms. We find that the sampling artifact suppresses the $z=0$ E-mode velocity power spectrum by $\sim 10\%$ at $k=0.1h/$Mpc , for samples with number density $10^{-3}({\rm Mpc}/h)^{-3}$. This suppression becomes larger for larger $k$ and for sparser samples. We argue that, this source of systematic errors in peculiar velocity cosmology, albeit severe, can be self-calibrated in the framework of our theoretical modelling. We also work out the sampling artifact in the density-velocity cross power spectrum measurement. More robust evaluation of related statistics through simulations will be presented in a companion paper (Zheng, Zhang & Jing, 2015, PRD). We also argue that similar sampling artifact exists in other velocity assignment methods and hence must be carefully corrected to avoid systematic bias in peculiar velocity cosmology.

preprint2015arXiv

Sampling Artifact in Volume Weighted Velocity Measurement.--- II. Detection in simulations and comparison with theoretical modelling

Measuring the volume weighted velocity power spectrum suffers from a severe systematic error, due to imperfect sampling of the velocity field from inhomogeneous distribution of dark matter particles/halos in simulations or galaxies with velocity measurement. This "sampling artifact" depends on both the mean particle number density $\bar{n}_P$ and the intrinsic large scale structure (LSS) fluctuation in the particle distribution. (1) We report robust detection of this sampling artifact in N-body simulations. It causes $\sim 12$% underestimation of the velocity power spectrum at $k=0.1$h/Mpc for samples with $\bar{n}_P=6\times10^{-3}$ (Mpc/h)$^{-3}$. This systematic underestimation increases with decreasing $\bar{n}_P$ and increasing $k$. Its dependence on the intrinsic LSS fluctuations is also robustly detected. (2) All these findings are expected by our theoretical modelling in paper I \cite{Zhang14}. In particular, the leading order theoretical approximation agrees quantitatively well with simulation result for $\bar{n}_P\gtrsim6\times 10^{-4}$(Mpc/h)$^{-3}$. Furthermore, we provide an ansatz to take high order terms into account. It improves the model accuracy to $\lesssim1$% at $k\lesssim0.1$h/Mpc over 3 orders of magnitude in $\bar{n}_P$ and over typical LSS clustering from $z=0$ to $z=2$. (3) The sampling artifact is determined by the deflection ${\bf D}$ field, which is straightforwardly available in both simulations and data of galaxy velocity. Hence the sampling artifact in the velocity power spectrum measurement can be self-calibrated within our framework. By applying such self-calibration in simulations, it becomes promising to determine the {\it real} large scale velocity bias of $10^{13}M_\odot$ halos with $\sim 1$% accuracy, and that of lower mass halos by better accuracy. ...[abridged]

preprint2015arXiv

Testing eternal inflation with the kinetic Sunyaev Zel'dovich effect

Perhaps the most controversial idea in modern cosmology is that our observable universe is contained within one bubble among many, all inhabiting the eternally inflating multiverse. One of the few way to test this idea is to look for evidence of the relic inhomogeneities left by the collisions between other bubbles and our own. Such relic inhomogeneities induces a coherent bulk flow over gigaparsec scales. Therefore, bubble collisions leave unique imprints in the cosmic microwave background (CMB) through the kinetic Sunyaev Zel'dovich (kSZ) effect, temperature anisotropies induced by the scattering of photons from coherently moving free electrons in the diffuse intergalactic medium. The kSZ signature produced by bubble collisions has a unique directional dependence and is tightly correlated with the galaxy distribution; it can therefore be distinguished from other contributions to the CMB anisotropies. An important advantage of the kSZ signature is that it peaks on arcminute angular scales, where the limiting factors in making a detection are instrumental noise and foreground subtraction. This is in contrast to the collision signature in the primary CMB, which peaks on angular scales much larger than one degree, and whose detection is therefore limited by cosmic variance. In this paper, we examine the prospects for probing the inhomogeneities left by bubble collisions using the kSZ effect. We provide a forecast for detection using cross-correlations between CMB and galaxy surveys, finding that the detectability using the kSZ effect can be competitive with constraints from CMB temperature and polarization data.

preprint2015arXiv

Testing the tensor-vector-scalar Theory with the latest cosmological observations

The tensor-vector-scalar (TeVeS) model is considered a viable theory of gravity. It produces the Milgrom's modified Newtonian dynamics in the nonrelativistic weak field limit and is free from ghosts. This model has been tested against various cosmological observations. Here we investigate whether new observations such as the galaxy velocity power spectrum measured by 6dF and the kinetic Sunyaev Zel'dovich effect power spectrum measured by ACT/SPT can put further constraints on the TeVeS model. Furthermore, we perform the test of TeVeS cosmology with a sterile neutrino by confronting to Planck data, and find that it is ruled out by cosmic microwave background measurements from the Planck mission.

preprint2015arXiv

The source-lens clustering effect in the context of lensing tomography and its self-calibration

Cosmic shear can only be measured where there are galaxies. This source-lens clustering (SLC) effect has two sources, intrinsic source clustering and cosmic magnification (magnification/size bias). Lensing tomography can suppress the former. However, this reduction is limited by the existence of photo-z error and nonzero redshift bin width. Furthermore, SLC induced by cosmic magnification cannot be reduced by lensing tomography. Through N-body simulations, we quantify the impact of SLC on the lensing power spectrum in the context of lensing tomography. We consider both the standard estimator and the pixel-based estimator. We find that none of them can satisfactorily handle both sources of SLC. (1) For the standard estimator, SLC induced by both sources can bias the lensing power spectrum by O(1)-O(10)%. Intrinsic source clustering also increases statistical uncertainties in the measured lensing power spectrum. However, the standard estimator suppresses intrinsic source clustering in the cross-spectrum. (2) In contrast, the pixel-based estimator suppresses SLC through cosmic magnification. However, it fails to suppress SLC through intrinsic source clustering and the measured lensing power spectrum can be biased low by O(1)-O(10)%. In short, for typical photo-z errors (sigma_z/(1+z)=0.05) and photo-z bin sizes (Delta_z^P=0.2), SLC alters the lensing E-mode power spectrum by 1-10%, with ell~10^3$ and z_s~1 being of particular interest to weak lensing cosmology. Therefore the SLC is a severe systematic for cosmology in Stage-IV lensing surveys. We present useful scaling relations to self-calibrate the SLC effect.

preprint2013arXiv

Nonlinearities in modified gravity cosmology. II. Impacts of modified gravity on the halo properties

The statistics of dark matter halos is an essential component of understanding the nonlinear evolution in modified gravity cosmology. Based on a series of modified gravity N-body simulations, we investigate the halo mass function, concentration and bias. We model the impact of modified gravity by a single parameter ζ, which determines the enhancement of particle acceleration with respect to GR, given the identical mass distribution (ζ=1 in GR). We select snapshot redshifts such that the linear matter power spectra of different gravity models are identical, in order to isolate the impact of gravity beyond modifying the linear growth rate. At the baseline redshift corresponding to z_S=1.2 in the standard ΛCDM, for a 10% deviation from GR(|ζ-1|=0.1), the measured halo mass function can differ by about 5-10%, the halo concentration by about 10-20%, while the halo bias differs significantly less. These results demonstrate that the halo mass function and/or the halo concentration are sensitive to the nature of gravity and may be used to make interesting constraints along this line.

preprint2013arXiv

Peculiar Velocity Decomposition, Redshift Space Distortion and Velocity Reconstruction in Redshift Surveys - I. The Methodology

Massive spectroscopic surveys will measure the redshift space distortion (RSD) induced by galaxy peculiar velocity to unprecedented accuracy. We develop a new method to improve the RSD modeling and to carry out robust reconstruction of the 3D large scale peculiar velocity through galaxy redshift surveys, in light of RSD. (1) We propose a mathematically unique and physically motivated decomposition of peculiar velocity into three eigen-components of different origins, different scale dependences and different impacts on RSD. (2) This decomposition has the potential to simplify and improve the RSD modeling. We derive a new formula for the redshift space power spectrum. Under the velocity decomposition scheme, all high order Gaussian corrections and non-Gaussian correction of order $δ^3$ can be taken into account without introducing extra model uncertainties. We find that the leading order contribution $\propto (1+f\tilde{W}(k)u^2)^2$ differs from the Kaiser formula by a function $\tilde{W}(k)\leq 1$. We hence identify a significant systematic error causing underestimation of the structure growth parameter $f$ by as much as $O(10%)$ even at relatively large scale $k=0.1h/$Mpc. The velocity decomposition also reveals the three origins of the finger of God (FOG) effect. (3) The velocity decomposition clarifies issues in peculiar velocity reconstruction through 3D galaxy distribution. We discuss two possible ways to carry out the velocity reconstruction. Both use the otherwise troublesome RSD in velocity reconstruction as a valuable source of information. Both have the advantage to render the reconstruction of a stochastic 3D field into the reconstruction of a deterministic window function $W^s(k,u)$ of limited degrees of freedom...[abridged]

preprint2013arXiv

The effect of Dark Matter and Dark Energy interactions on the peculiar velocity field and the kinetic Sunyaev-Zel'dovich effect

The interaction between Dark Matter and Dark Energy has been proposed as a mechanism to alleviate the coincidence problem. We analyze the effect of the interaction on the evolution of the gravitational field and propose two new observables based on its effect on the matter peculiar velocity field. We find that for different model parameters the matter peculiar velocity could be factor 2 times larger or 5 times smaller than the amplitude of velocity perturbations in the concordance LCDM cosmological model at the same scales. We compare the effect on the peculiar velocities with those on the Integrated Sachs-Wolfe effect and we show that velocities can potentially provide constraints on the strength of the interaction stronger than those currently available. We show that the current upper limits on the amplitude of the kinetic Sunyaev-Zel'dovich power spectrum provide constraints on the coupling within the dark sectors that are consistent with those obtained previously from the CMB and galaxy clusters. In particular, we show that current upper limits from the Atacama Cosmology Telescope and the South Pole Telescope favor Dark Energy decaying into Dark Matter, as required to solve the coincidence problem.

preprint2013arXiv

Weak lensing reconstruction through cosmic magnification. II. Improved power spectrum determination and map-making

The existence of galaxy intrinsic clustering severely hampers the weak lensing reconstruction from cosmic magnification. In paper I \citep{Yang2011}, we proposed a minimal variance estimator to overcome this problem. By utilizing the different dependences of cosmic magnification and galaxy intrinsic clustering on galaxy flux, we demonstrated that the otherwise overwhelming galaxy intrinsic clustering can be significantly suppressed such that lensing maps can be reconstructed with promising accuracy. This procedure relies heavily on the accuracy of determining the galaxy bias from the same data. Paper I adopts an iterative approach, which degrades toward high redshift. The current paper presents an alternative method, improving over paper I. We prove that the measured galaxy clustering between flux bins allows for simultaneous determination of the lensing power spectrum and the flux dependence of galaxy bias, at this redshift bin. Comparing to paper I, the new approach is not only more straightforward, but also more robust. It identifies an ambiguity in determining the galaxy bias and further discovers a mathematically robust way to suppress this ambiguity to non-negligible level ($\sim 0.1%$). The accurately determined galaxy bias can then be applied to the minimal variance estimator proposed in paper I to improve the lensing map-making. The gain at high redshift is significant. These maps can be used to measure other statistics, such as cluster finding and peak statistics. Furthermore, by including galaxy clustering measurement between different redshift bins, we can also determine the lensing cross power spectrum between these bins, up to a small and correctable multiplicative factor.

preprint2012arXiv

Bulk flow of halos in ΛCDM simulation

Analysis of the Pangu N-body simulation validates that the bulk flow of halos follows a Maxwellian distribution which variance is consistent with the prediction of the linear theory of structure formation. We propose that the consistency between the observed bulk velocity and theories should be examined at the effective scale of the radius of a spherical top-hat window function yielding the same smoothed velocity variance in linear theory as the sample window function does. We compared some recently estimated bulk flows from observational samples with the prediction of the ΛCDM model we used; some results deviate from expectation at a level of ~ 3σbut the discrepancy is not as severe as previously claimed. We show that bulk flow is only weakly correlated with the dipole of the internal mass distribution, the alignment angle between the mass dipole and the bulk flow has a broad distribution peaked at ~ 30-50 deg., and also that the bulk flow shows little dependence on the mass of the halos used in the estimation. In a simulation of box size 1Gpc/h, for a cell of radius 100 Mpc/h the maximal bulk velocity is >500 km/s, dipoles of the environmental mass outside the cell are not tightly aligned with the bulk flow, but are rather located randomly around it with separation angles ~ 20-40 deg. In the fastest cell there is a slightly smaller number of low-mass halos; however halos inside are clustered more strongly at scales > ~ 20 Mpc/h, which might be a significant feature since the correlation between bulk flow and halo clustering actually increases in significance beyond such scales.

preprint2012arXiv

Gaussianizing the non-Gaussian lensing convergence field II: the applicability to noisy data

In paper I (Yu et al. [1]), we show through N-body simulation that a local monotonic Gaussian transformation can significantly reduce non-Gaussianity in a noise-free lensing convergence field. This makes the Gaussianization a promising theoretical tool to understand high-order lensing statistics. Here we present a study of its applicability in lensing data analysis, in particular when shape measurement noise is presented in lensing convergence maps. (i) We find that shape measurement noise significantly degrades the Gaussianization performance and the degradation increases for shallower surveys. (ii) The Wiener filter is efficient in reducing the impact of shape measurement noise. The Gaussianization of the Wiener-filtered lensing maps is able to suppress skewness, kurtosis, and the 5th- and 6th-order cumulants by a factor of 10 or more. It also works efficiently to reduce the bispectrum to zero.

preprint2012arXiv

Observational Consequences of Dark Energy Decay

We consider the generic scenario of dark energy which arises through the latent heat of a hidden sector first order cosmological phase transition. This field could account for the extra radiation degree of freedom suggested by the CMB. We present the bubble nucleation solution for the viscous limit. The decay rate of the field is constrained by published KSZ data, and may be an explanation of current excess ISW correlations. Cross correlation of current and future surveys can further constrain or test the parameter space. The decay model is plausibly in the observable range, and avoids anthropic problems. This class of models is not well constrained by the popular dark energy figure of merit.

preprint2012arXiv

Weak lensing reconstruction through cosmic magnification I: a minimal variance map reconstruction

We present a concept study on weak lensing map reconstruction through the cosmic magnification effect in galaxy number density distribution. We propose a minimal variance linear estimator to minimize both the dominant systematical and statistical errors in the map reconstruction. It utilizes the distinctively different flux dependences to separate the cosmic magnification signal from the overwhelming galaxy intrinsic clustering noise. It also minimizes the shot noise error by an optimal weighting scheme on the galaxy number density in each flux bin. Our method is in principle applicable to all galaxy surveys with reasonable redshift information. We demonstrate its applicability against the planned Square Kilometer Array survey, under simplified conditions. Weak lensing maps reconstructed through our method are complementary to that from cosmic shear and CMB and 21cm lensing. They are useful for cross checking over systematical errors in weak lensing reconstruction and for improving cosmological constraints.

preprint2011arXiv

Confirmation of the Copernican principle at Gpc radial scale and above from the kinetic Sunyaev Zel'dovich effect power spectrum

The Copernican principle, a cornerstone of modern cosmology, remains largely unproven at Gpc radial scale and above. Here we will show that, violations of this type will inevitably cause a first order anisotropic kinetic Sunyaev Zel'dovich (kSZ) effect. If large scale radial inhomogeneities have amplitude large enough to explain the "dark energy" phenomena, the induced kSZ power spectrum will be much larger than the ACT/SPT upper limit. This single test confirms the Copernican principle and rules out the adiabatic void model as a viable alternative to dark energy.

preprint2011arXiv

Gaussianizing the non-Gaussian lensing convergence field I: the performance of the Gaussianization

Motivated by recent works of Neyrinck et al. 2009 and Scherrer et al. 2010, we proposed a Gaussianization transform to Gaussianize the non-Gaussian lensing convergence field $κ$. It performs a local monotonic transformation $κ\rightarrow y$ pixel by pixel to make the unsmoothed one-point probability distribution function of the new variable $y$ Gaussian. We tested whether the whole $y$ field is Gaussian against N-body simulations. (1) We found that the proposed Gaussianization suppresses the non-Gaussianity by orders of magnitude, in measures of the skewness, the kurtosis, the 5th- and 6th-order cumulants of the $y$ field smoothed over various angular scales relative to that of the corresponding smoothed $κ$ field. The residual non-Gaussianities are often consistent with zero within the statistical errors. (2) The Gaussianization significantly suppresses the bispectrum. Furthermore, the residual scatters around zero, depending on the configuration in the Fourier space. (3) The Gaussianization works with even better performance for the 2D fields of the matter density projected over $\sim 300 \mpch$ distance interval centered at $z\in(0,2)$, which can be reconstructed from the weak lensing tomography. (4) We identified imperfectness and complexities of the proposed Gaussianization. We noticed weak residual non-Gaussianity in the $y$ field. We verified the widely used logarithmic transformation as a good approximation to the Gaussianization transformation. However, we also found noticeable deviations.

preprint2011arXiv

Spherical Collapse in f(R) Gravity

We use 1-dimensional numerical simulations to study spherical collapse in the f(R) gravity models. We include the nonlinear coupling of the gravitational potential to the scalar field in the theory and use a relaxation scheme to follow the collapse. We find an unusual enhancement in density near the virial radius which may provide observable tests of gravity. We also use the estimated collapse time to calculate the critical overdensity $δ_c$ used in calculating the mass function and bias of halos. We find that analytical approximations previously used in the literature do not capture the complexity of nonlinear spherical collapse.

preprint2011arXiv

The exact analytical solution of the linear structure growth rate in LCDM cosmology and its cosmological applications

We derive the exact analytical solution of the linear structure growth rate in LCDM cosmology with flat or curved geometry, under the Newtonian gauge. Unlike the well known solution under the Newtonian limit (Heath 1977), our solution takes all general relativistic corrections into account and is hence valid at both the sub- and super-horizon scales. With this exact solution, we evaluate cosmological impacts induced by these relativistic corrections. (1) General relativistic corrections alter the density growth from z=100 to z=0 by 10% at k=0.01h/Mpc and the impact becomes stronger toward larger scales. We caution the readers that the overdensity is not gauge invariant and the above statement is restrained to the Newtonian gauge. (2) Relativistic corrections introduce a k^{-2} scale dependence in the density fluctuation. It mimics a primordial non-Gaussianity of the local type with $f^{\rm local}_{\rm NL}\sim 1$. This systematical error may become non-negligible for future all sky deep galaxy surveys. (3) Cosmological simulations with box size greater than 1Gpc are also affected by these relativistic corrections. We provide a post-processing recipe to correct for these effects. (4) These relativistic corrections affect the redshift distortion. However, at redshifts and scales relevant to redshift distortion measurements, such effect is negligible.

preprint2011arXiv

The kinetic SZ tomography with spectroscopic redshift surveys

The kinetic Sunyaev Zel'dovich effect (kSZ) effect is a potentially powerful probe to the missing baryons. However, the kSZ signal is overwhelmed by various contaminations and the cosmological application is hampered by loss of redshift information due to the projection effect. We propose a kSZ tomography method to alleviate these problems, with the aid of galaxy spectroscopic redshift surveys. We propose to estimate the large scale peculiar velocity through the 3D galaxy distribution, weigh it by the 3D galaxy density and adopt the product projected along the line of sight with a proper weighting as an estimator of the true kSZ temperature fluctuation $Θ$. We thus propose to measure the kSZ signal through the $\HatΘ$-$Θ$ cross correlation. This approach has a number of advantages (see details in the abstract of the paper). We test the proposed kSZ tomography against non-adiabatic and adiabatic hydrodynamical simulations. We confirm that $\hatΘ$ is indeed tightly correlated with $Θ$ at $k\la 1h/$Mpc, although nonlinearities in the density and velocity fields and nonlinear redshift distortion do weaken the tightness of the $\hatΘ$-$Θ$ correlation. We further quantify the reconstruction noise in $\HatΘ$ from galaxy distribution shot noise. Based on these results, we quantify the applicability of the proposed kSZ tomography for future surveys. We find that, in combination with the BigBOSS-N spectroscopic redshift survey, the PLANCK CMB experiment will be able to detect the kSZ with an overall significance of $\sim 50σ$ and further measure its redshift distribution at many redshift bins over $0<z<2$.

preprint2011arXiv

The optimal weighting function for cosmic magnification measurement through foreground galaxy-background galaxy (quasar) cross correlation

Cosmic magnification has been detected through cross correlation between foreground and background populations (galaxies or quasars). It has been shown that weighing each background object by its $α-1$ can significantly improve the cosmic magnification measurement \citep{Menard02,Scranton05}. Here, $α$ is the logarithmic slope of the luminosity function of background populations. However, we find that this weighting function is optimal only for sparse background populations in which intrinsic clustering is negligible with respect to shot noise. We derive the optimal weighting function for general case including scale independent and scale dependent weights. The optimal weighting function improves the S/N (signal to noise ratio) by $\sim 20%$ for a BigBOSS-like survey and the improvement can reach a factor of $\sim 2$ for surveys with much denser background populations.

preprint2011arXiv

The thermal SZ tomography

The thermal Sunyaev-Zel'dovich (tSZ) effect directly measures the thermal pressure of free electrons integrated along the line of sight and thus contains valuable information on the thermal history of the universe. However, the redshift information is entangled in the projection along the line of sight. This projection effect severely degrades the power of the tSZ effect to reconstruct the thermal history. We investigate the tSZ tomography technique to recover this otherwise lost redshift information by cross correlating the tSZ effect with galaxies of known redshifts, or alternatively with matter distribution reconstructed from weak lensing tomography. We investigate in detail the 3D distribution of the gas thermal pressure and its relation with the matter distribution, through our adiabatic hydrodynamic simulation and the one with additional gastrophysics including radiative cooling, star formation and supernova feedback. (1) We find a strong correlation between the gas pressure and matter distribution, with a typical cross correlation coefficient r ~ 0.7 at k . 3h/Mpc and z < 2. This tight correlation will enable robust cross correlation measurement between SZ surveys such as Planck, ACT and SPT and lensing surveys such as DES and LSST, at ~20-100σ level. (2) We propose a tomography technique to convert the measured cross correlation into the contribution from gas in each redshift bin to the tSZ power spectrum. Uncertainties in gastrophysics may affect the reconstruction at ~ 2% level, due to the ~ 1% impact of gastrophysics on r, found in our simulations. However, we find that the same gastrophysics affects the tSZ power spectrum at ~ 40% level, so it is robust to infer the gastrophysics from the reconstructed redshift resolved contribution.

preprint2010arXiv

A proposal on the galaxy intrinsic alignment self-calibration in weak lensing surveys

The galaxy intrinsic alignment causes the galaxy ellipticity-ellipticity power spectrum between two photometric redshifts to decrease faster with respect to the redshift separation $Δz^P$, for fixed mean redshift. This offers a valuable diagnosis on the intrinsic alignment. We show that the distinctive dependences of the GG, II and GI correlations on $Δz^P$ over the range $|Δz^P|\la 0.2$ can be understood robustly without strong assumptions on the intrinsic alignment. This allows us to measure the intrinsic alignment within each conventional photo-z bin of typical size $\ga 0.2$, through lensing tomography of photo-z bin size $\sim 0.01$. Both the statistical and systematical errors in the lensing cosmology can be reduced by this self-calibration technique.

preprint2010arXiv

Magnification as a Tool in Weak Lensing

Weak lensing surveys exploit measurements of galaxy ellipticities. These measurements are subject to errors which degrade the cosmological information that can be extracted from the surveys. Here we propose a way of using the galaxy data themselves to calibrate the measurement errors. In particular, the cosmic shear field, which causes the galaxies to appear elliptical, also changes their sizes and fluxes. Information about the sizes and fluxes of the galaxies can be added to the shape information to obtain more robust information about the cosmic shear field. The net result will be tighter constraints on cosmological parameters such as those which describe dark energy.

preprint2010arXiv

Nonlinearities in modified gravity cosmology I: signatures of modified gravity in the nonlinear matter power spectrum

A large fraction of cosmological information on dark energy and gravity is encoded in the nonlinear regime. Precision cosmology thus requires precision modeling of nonlinearities in general dark energy and modified gravity models. We modify the Gadget-2 code and run a series of N-body simulations on modified gravity cosmology to study the nonlinearities. The modified gravity model that we investigate in the present paper is characterized by a single parameter ζ, which determines the enhancement of particle acceleration with respect to general relativity (GR), given the identical mass distribution (ζ= 1 in GR). The first nonlinear statistics we investigate is the nonlinear matter power spectrum at k < 3h/Mpc, which is the relevant range for robust weak lensing power spectrum modeling at l < 2000. In this study, we focus on the relative difference in the nonlinear power spectra at corresponding redshifts where different gravity models have the same linear power spectra. This particular statistics highlights the imprint of modified gravity in the nonlinear regime and the importance to include the nonlinear regime in testing GR. By design, it is less susceptible to the sample variance and numerical artifacts. We adopt a mass assignment method based on wavelet to improve the power spectrum measurement. We run a series of tests to determine the suitable simulation specifications (particle number, box size and initial redshift). We find that, the nonlinear power spectra can differ by ~30% for 10% deviation from GR (|ζ-1| = 0.1) where the rms density fluctuations reach 10. This large difference, on one hand, shows the richness of information on gravity in the corresponding scales, and on the other hand, invalidates simple extrapolations of some existing fitting formulae to modified gravity cosmology.

preprint2010arXiv

Self calibration of gravitational shear-galaxy intrinsic ellipticity correlation in weak lensing surveys

The galaxy intrinsic alignment is a severe challenge to precision cosmic shear measurement. We propose to self-calibrate the induced gravitational shear-galaxy intrinsic ellipticity correlation (the GI correlation, \citealt{Hirata04b}) in weak lensing surveys with photometric redshift measurement. (1) We propose a method to extract the intrinsic ellipticity-galaxy density cross correlation (I-g) from the galaxy ellipticity-density measurement in the same redshift bin. (2) We also find a generic scaling relation to convert the extracted I-g correlation to the demanded GI correlation. We perform concept study under simplified conditions and demonstrate its capability to significantly reduce the GI contamination. We discuss the impact of various complexities on the two key ingredients of the self-calibration technique, namely the method to extract the I-g correlation and the scaling relation between the I-g and the GI correlation. We expect none of them is likely able to completely invalidate the proposed self-calibration technique.

preprint2010arXiv

Self calibration of photometric redshift scatter in weak lensing surveys

Photo-z errors, especially catastrophic errors, are a major uncertainty for precision weak lensing cosmology. We find that the shear-(galaxy number) density and density-density cross correlation measurements between photo-z bins, available from the same lensing surveys, contain valuable information for self-calibration of the scattering probabilities between the true-z and photo-z bins. The self-calibration technique we propose does not rely on cosmological priors nor parameterization of the photo-z probability distribution function, and preserves all of the cosmological information available from shear-shear measurement. We estimate the calibration accuracy through the Fisher matrix formalism. We find that, for advanced lensing surveys such as the planned stage IV surveys, the rate of photo-z outliers can be determined with statistical uncertainties of 0.01-1% for $z<2$ galaxies. Among the several sources of calibration error that we identify and investigate, the {\it galaxy distribution bias} is likely the most dominant systematic error, whereby photo-z outliers have different redshift distributions and/or bias than non-outliers from the same bin. This bias affects all photo-z calibration techniques based on correlation measurements. Galaxy bias variations of $O(0.1)$ produce biases in photo-z outlier rates similar to the statistical errors of our method, so this galaxy distribution bias may bias the reconstructed scatters at several-$σ$ level, but is unlikely to completely invalidate the self-calibration technique.

preprint2010arXiv

The dark flow induced small scale kinetic Sunyaev Zel'dovich effect

Recently Kashlinsky et al. 2008, 2010 reported a discovery of a $\sim 10^3$ km/$s$ bulk flow of the universe out to $z\simeq 0.3$, through the dark flow induced CMB dipole in directions of clusters. We point out that, if this dark flow exists, it will also induce observable CMB temperature fluctuations at multipole $\ell\sim 10^3$, through modulation of the inhomogeneous electron distribution on the uniform dark flow. The induced small scale kinetic Sunyaev Zel'dovich (SZ) effect will reach $\sim 1\muk^2$ at multipole $10^3\la \ell\la 10^4$, only a factor of $\sim 2$ smaller than the conventional kinetic SZ effect. Furthermore, it will be correlated with the large scale structure (LSS) and its correlation with 2MASS galaxy distribution reaches $0.3 μ$K at $\ell=10^3$, under a directional dependent optimal weighting scheme. We estimate that, WMAP plus 2MASS should already be able to detect this dark flow induced small scale kinetic SZ effect with $\sim 6σ$ confidence. Deeper galaxy surveys such as SDSS can further improve the measurement. Planck plus existing galaxy surveys can reach $\ga 14σ$ detection. Existing CMB-LSS cross correlation measurements shall be reanalyzed to test the existence of the dark flow and, if it exists, shall be used to eliminate possible bias on the integrated Sachs-Wolfe effect measurement through the CMB-LSS cross correlation.

preprint2010arXiv

The ergodicity bias in the observed galaxy distribution

The spatial distribution of galaxies we observed is subject to the given condition that we, human beings are sitting right in a galaxy -- the Milky Way. Thus the ergodicity assumption is questionable in interpretation of the observed galaxy distribution. The resultant difference between observed statistics (volume average) and the true cosmic value (ensemble average) is termed as the ergodicity bias. We perform explicit numerical investigation of the effect for a set of galaxy survey depths and near-end distance cuts. It is found that the ergodicity bias in observed two- and three-point correlation functions in most cases is insignificant for modern analysis of samples from galaxy surveys and thus close a loophole in precision cosmology. However, it may become non-negligible in certain circumstances, such as those applications involving three-point correlation function at large scales of local galaxy samples. Thus one is reminded to take extra care in galaxy sample construction and interpretation of the statistics of the sample, especially when the characteristic redshift is low.

preprint2008arXiv

A new method of measuring the cluster peculiar velocity power spectrum

We propose to use spatial correlations of the kinetic Sunyaev-Zeldovich (KSZ) flux as an estimator of the peculiar velocity power spectrum. In contrast with conventional techniques, our new method does not require measurements of the thermal SZ signal or the X-ray temperature. Moreover, this method has the special advantage that the expected systematic errors are always sub-dominant to statistical errors on all scales and redshifts of interest. We show that future large sky coverage KSZ surveys may allow a peculiar velocity power spectrum estimates of an accuracy reaching ~10%.

preprint2007arXiv

Detecting First Star Lyman-$α$ Spheres through Gravitational Telescopes

Lyman-$α$ spheres, i.e. regions around the first stars which are illuminated by Lyman-$α$ photons and show 21cm absorption feature against the CMB, are smoking guns at the dawn of the reionization epoch. Though overwhelming radio foreground makes their detections extremely difficult, we pointed out that, strong gravitational lensing can significantly improve their observational feasibility. Since Lyman-$α$ spheres have ~10" sizes, comparable to the caustic size of galaxy clusters, individual images of each strongly lensed Lyman-$α$ sphere often merge together and form single structures in the 21cm sky with irregular shapes. Using high-resolution N-body LCDM simulations, we found that the lensing probability to have magnification bigger than 10 is ~10^{-5}. This results in $\ga 10^6$ strongly lensed Lyman-$α$ spheres across the sky, which should be the primary targets for first detections of Lyman-$α$ spheres. Although the required total radio array collecting area for their detection is large (~100 km^2), the design of long fixed cylindrical reflectors can significantly reduce the total cost of such array to the level of the square kilometer array (SKA) and makes the detection of these very first objects feasible.

preprint2007arXiv

Impact of astrophysical processes on the gamma-ray background from dark matter annihilations

We study the impact of astrophysical processes on the gamma-ray background produced by the annihilation of dark matter particles in cosmological halos, with particular attention to the consequences of the formation of supermassive black holes. In scenarios where these objects form adiabatically from the accretion of matter on small seeds, dark matter is first compressed into very dense ``spikes'', then its density progressively decreases due to annihilations and scattering off of stellar cusps. With respect to previous analyses, based on non-evolving halos, the predicted annihilation signal is higher and significantly distorted at low energies, reflecting the large contribution to the total flux from unevolved spikes at high redshifts. The peculiar spectral feature arising from the specific redshift distribution of the signal, would discriminate the proposed scenario from more conventional astrophysical explanations. We discuss how this affects the prospects for detection and demonstrate that the gamma-ray background from DM annihilations might be detectable even in absence of a signal from the Galactic center.

preprint2005arXiv

The Concordance Cosmic Star Formation Rate: Implications from and for the Supernova Neutrino and Gamma Ray Backgrounds

We constrain the Cosmic Star Formation Rate (CSFR) by requiring that massive stars produce the observed UV, optical, and IR light while at the same time not overproduce the Diffuse Supernova Neutrino Background as bounded by Super-Kamiokande. With the massive star component so constrained we then show that a reasonable choice of stellar Initial Mass Function and other parameters results in SNIa rates and iron yields in good agreement with data. In this way we define a `concordance' CSFR that predicts the optical SNII rate and the SNIa contribution to the MeV Cosmic Gamma-Ray Background. The CSFR constrained to reproduce these and other proxies of intermediate and massive star formation is more clearly delineated than if it were measured by any one technique and has the following testable consequences: (1) SNIa contribute only a small fraction of the MeV Cosmic Gamma-Ray Background, (2) massive star core-collapse is nearly always accompanied by a successful optical SNII, and (3) the Diffuse Supernova Neutrino Background is tantalizingly close to detectability.

preprint2002arXiv

Hydrodynamical Simulations of the IGM at High Mach Numbers

We present a new approach to doing Eulerian computational fluid dynamics that is designed to work at high Mach numbers encountered in hydrodynamical simulations of the IGM. In conventional Eulerian CFD, the thermal energy is poorly tracked in supersonic bulk flows where local fluid variables cannot be accurately separated from the much larger bulk flow components. We described a method in which local fluid quantities can be directly tracked and the Eulerian fluid equations solved in a local frame moving with the flow. The new algorithm has been used to run large hydrodynamical simulations on a 1024^3 grid to study the kinetic SZ effect. The KSZ power spectrum is broadly peaked at l~10^4 with temperature fluctuations on micro Kelvin levels.