Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOn measuring the covariance matrix of the nonlinear power spectrum from simulations
77
0
0.0
(
0
)
Authors
A. J. S. Hamilton
Ask ChatGPT about the research
No Arabic abstract
We show how to estimate the covariance of the power spectrum of a statistically homogeneous and isotropic density field from a single periodic simulation, by applying a set of weightings to the density field, and by measuring the scatter in power spectra between different weightings. We recommend a specific set of 52 weightings containing only combinations of fundamental modes, constructed to yield a minimum variance estimate of the covariance of power. Numerical tests reveal that at nonlinear scales the variance of power estimated by the weightings method substantially exceeds that estimated from a simple ensemble method. We argue that the discrepancy is caused by beat-coupling, in which products of closely spaced Fourier modes couple by nonlinear gravitational growth to the beat mode between them. Beat-coupling appears whenever nonlinear power is measured from Fourier modes with a finite spread of wavevector, and is therefore present in the weightings method but not the ensemble method. Beat-coupling inevitably affects real galaxy surveys, whose Fourier modes have finite width. Surprisingly, the beat-coupling contribution dominates the covariance of power at nonlinear scales, so that, counter-intuitively, it is expected that the covariance of nonlinear power in galaxy surveys is dominated not by small scale structure, but rather by beat-coupling to the largest scales of the survey.
rate research
Read More
We seek to improve estimates of the power spectrum covariance matrix from a limited number of simulations by employing a novel statistical technique known as shrinkage estimation. The shrinkage technique optimally combines an empirical estimate of the covariance with a model (the target) to minimize the total mean squared error compared to the true underlying covariance. We test this technique on N-body simulations and evaluate its performance by estimating cosmological parameters. Using a simple diagonal target, we show that the shrinkage estimator significantly outperforms both the empirical covariance and the target individually when using a small number of simulations. We find that reducing noise in the covariance estimate is essential for properly estimating the values of cosmological parameters as well as their confidence intervals. We extend our method to the jackknife covariance estimator and again find significant improvement, though simulations give better results. Even for thousands of simulations we still find evidence that our method improves estimation of the covariance matrix. Because our method is simple, requires negligible additional numerical effort, and produces superior results, we always advocate shrinkage estimation for the covariance of the power spectrum and other large-scale structure measurements when purely theoretical modeling of the covariance is insufficient.
We use 5000 cosmological N-body simulations of 1(Gpc/h)^3 box for the concordance LCDM model in order to study the sampling variances of nonlinear matter power spectrum. We show that the non-Gaussian errors can be important even on large length scales relevant for baryon acoustic oscillations (BAO). Our findings are (1) the non-Gaussian errors degrade the cumulative signal-to-noise ratios (S/N) for the power spectrum amplitude by up to a factor of 2 and 4 for redshifts z=1 and 0, respectively. (2) There is little information on the power spectrum amplitudes in the quasi-nonlinear regime, confirming the previous results. (3) The distribution of power spectrum estimators at BAO scales, among the realizations, is well approximated by a Gaussian distribution with variance that is given by the diagonal covariance component. (4) For the redshift-space power spectrum, the degradation in S/N by non-Gaussian errors is mitigated due to nonlinear redshift distortions. (5) For an actual galaxy survey, the additional shot noise contamination compromises the cosmological information inherent in the galaxy power spectrum, but also mitigates the impact of non-Gaussian errors. The S/N is degraded by up to 30% for a WFMOS-type survey. (6) The finite survey volume causes additional non-Gaussian errors via the correlations of long-wavelength fluctuations with the fluctuations we want to measure, further degrading the S/N values by about 30% even at high redshift z=3.
We compute the effects induced by the use of small CMB maps on the measurement of the $cl{l}$ coefficients of the angular power spectrum and show that small systematic effects have to be taken into account. We also compute numerically the cosmic variance and covariance of the $cl{l}$ spectrum for various spherical cap like maps. Comparisons with simulations are presented. The calculations are done using the standard method based on the spherical harmonic transform or using the temperature angular correlation spectrum.
The large-scale structure of the Universe should soon be measured at high redshift during the Epoch of Reionization (EoR) through line-intensity mapping. A number of ongoing and planned surveys are using the 21 cm line to trace neutral hydrogen fluctuations in the intergalactic medium (IGM) during the EoR. These may be fruitfully combined with separate efforts to measure large-scale emission fluctuations from galactic lines such as [CII], CO, H-$alpha$, and Ly-$alpha$ during the same epoch. The large scale power spectrum of each line encodes important information about reionization, with the 21 cm power spectrum providing a relatively direct tracer of the ionization history. Here we show that the large scale 21 cm power spectrum can be extracted using only cross-power spectra between the 21 cm fluctuations and each of two separate line-intensity mapping data cubes. This technique is more robust to residual foregrounds than the usual 21 cm auto-power spectrum measurements and so can help in verifying auto-spectrum detections. We characterize the accuracy of this method using numerical simulations and find that the large-scale 21 cm power spectrum can be inferred to an accuracy of within 5% for most of the EoR, reaching 0.6% accuracy on a scale of $ksim0.1,text{Mpc}^{-1}$ at $left< x_i right> = 0.36$ ($z = 8.34$ in our model). An extension from two to $N$ additional lines would provide $N(N-1)/2$ cross-checks on the large-scale 21 cm power spectrum. This work strongly motivates redundant line-intensity mapping surveys probing the same cosmological volumes.
Primordial gravitational waves generated during inflation lead to the B-mode polarization in the cosmic microwave background and a stochastic gravitational wave background in the Universe. We will explore the current constraint on the tilt of primordial gravitational-wave spectrum, and forecast how the future observations can improve the current constraint.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
