Correlations in the Spatial Power Spectrum Inferred from Angular Clustering: Methods and Application to APM

We reconsider the inference of spatial power spectra from angular clustering data and show how to include correlations in both the angular correlation function and the spatial power spectrum. Inclusion of the full covariance matrices loosens the constraints on large-scale structure inferred from the APM survey by over a factor of two. We present a new inversion technique based on singular value decomposition that allows one to propagate the covariance matrix on the angular correlation function through to that of the spatial power spectrum and to reconstruct smooth power spectra without underestimating the errors. Within a parameter space of the CDM shape Gamma and the amplitude sigma_8, we find that the angular correlations in the APM survey constrain Gamma to be 0.19-0.37 at 68% confidence when fit to scales larger than k=0.2h Mpc^-1. A downturn in power at k<0.04h Mpc^-1 is significant at only 1-sigma. These results are optimistic as we include only Gaussian statistical errors and neglect any boundary effects.Comment: 37 pages, LaTex, 9 figures. Submitted to Ap

Faint counts as a function of morphological type in a hierarchical merger model

The unprecedented resolution of the refurbished Wide Field and Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST) has led to major advances in our understanding of galaxy formation. The high image quality in the Medium Deep Survey and Hubble Deep Field has made it possible, for the first time, to classify faint distant galaxies according to morphological type. These observations have revealed a large population of galaxies classed as irregulars or which show signs of recent merger activity. Their abundance rises steeply with apparent magnitude, providing a likely explanation for the large number of blue galaxies seen at faint magnitudes. We demonstrate that such a population arises naturally in a model in which structure forms hierarchically and which is dynamically dominated by cold dark matter. The number counts of irregular, spiral and elliptical galaxies as a function of magnitude seen in the HST data are well reproduced in this model.We present detailed predictions for the outcome of spectroscopic follow-up observations of the HST surveys. By measuring the redshift distributions of faint galaxies of different morphological types, these programmes will provide a test of the hierarchical galaxy formation paradigm and might distinguish between models with different cosmological parameters.Comment: 5 pages, 3 postscript figures included. To be published as a Letter in Monthly Notices of the RAS. Postscript version available at http://star-www.dur.ac.uk/~cmb/counts.htm

Estimation of Primordial Spectrum with post-WMAP 3 year data

In this paper we implement an improved (error sensitive) Richardson-Lucy deconvolution algorithm on the measured angular power spectrum from the WMAP 3 year data to determine the primordial power spectrum assuming different points in the cosmological parameter space for a flat LCDM cosmological model. We also present the preliminary results of the cosmological parameter estimation by assuming a free form of the primordial spectrum, for a reasonably large volume of the parameter space. The recovered spectrum for a considerably large number of the points in the cosmological parameter space has a likelihood far better than a `best fit' power law spectrum up to \Delta \chi^2_{eff} \approx -30. We use Discrete Wavelet Transform (DWT) for smoothing the raw recovered spectrum from the binned data. The results obtained here reconfirm and sharpen the conclusion drawn from our previous analysis of the WMAP 1st year data. A sharp cut off around the horizon scale and a bump after the horizon scale seem to be a common feature for all of these reconstructed primordial spectra. We have shown that although the WMAP 3 year data prefers a lower value of matter density for a power law form of the primordial spectrum, for a free form of the spectrum, we can get a very good likelihood to the data for higher values of matter density. We have also shown that even a flat CDM model, allowing a free form of the primordial spectrum, can give a very high likelihood fit to the data. Theoretical interpretation of the results is open to the cosmology community. However, this work provides strong evidence that the data retains discriminatory power in the cosmological parameter space even when there is full freedom in choosing the primordial spectrum.Comment: 13 pages, 4 figures, uses Revtex4, new analysis and results, references added, matches version accepted to Phys. Rev.

How BAO measurements can fail to detect quintessence

We model the nonlinear growth of cosmic structure in different dark energy models, using large volume N-body simulations. We consider a range of quintessence models which feature both rapidly and slowly varying dark energy equations of state, and compare the growth of structure to that in a universe with a cosmological constant. The adoption of a quintessence model changes the expansion history of the universe, the form of the linear theory power spectrum and can alter key observables, such as the horizon scale and the distance to last scattering. The difference in structure formation can be explained to first order by the difference in growth factor at a given epoch; this scaling also accounts for the nonlinear growth at the 15% level. We find that quintessence models which feature late $(z<2)$, rapid transitions towards $w=-1$ in the equation of state, can have identical baryonic acoustic oscillation (BAO) peak positions to those in $\Lambda$CDM, despite being very different from $\Lambda$CDM both today and at high redshifts $(z \sim 1000)$. We find that a second class of models which feature non-negligible amounts of dark energy at early times cannot be distinguished from $\Lambda$CDM using measurements of the mass function or the BAO. These results highlight the need to accurately model quintessence dark energy in N-body simulations when testing cosmological probes of dynamical dark energy.Comment: 10 pages, 7 figures, to appear in the Invisible Univers International Conference AIP proceedings serie

Cosmic Variance In the Transparency of the Intergalactic Medium After Reionization

Following the completion of cosmic reionization, the mean-free-path of ionizing photons was set by a population of Ly-limit absorbers. As the mean-free-path steadily grew, the intensity of the ionizing background also grew, thus lowering the residual neutral fraction of hydrogen in ionization equilibrium throughout the diffuse intergalactic medium (IGM). Ly-alpha photons provide a sensitive probe for tracing the distribution of this residual hydrogen at the end of reionization. Here we calculate the cosmic variance among different lines-of-sight in the distribution of the mean Ly-alpha optical depths. We find fractional variations in the effective post-reionization optical depth that are of order unity on a scale of ~100 co-moving Mpc, in agreement with observations towards high-redshift quasars. Significant contributions to these variations are provided by the cosmic variance in the density contrast on the scale of the mean-free-path for ionizing photons, and by fluctuations in the ionizing background induced by delayed or enhanced structure formation. Cosmic variance results in a highly asymmetric distribution of transmission through the IGM, with fractional fluctuations in Ly-alpha transmission that ar larger than in Ly-beta transmission.Comment: 7 pages 3 figures. Replaced with version accepted for publication in Ap

Biased-estimations of the Variance and Skewness

Nonlinear combinations of direct observables are often used to estimate quantities of theoretical interest. Without sufficient caution, this could lead to biased estimations. An example of great interest is the skewness $S_3$ of the galaxy distribution, defined as the ratio of the third moment \xibar_3 and the variance squared \xibar_2^2. Suppose one is given unbiased estimators for \xibar_3 and \xibar_2^2 respectively, taking a ratio of the two does not necessarily result in an unbiased estimator of $S_3$. Exactly such an estimation-bias affects most existing measurements of $S_3$. Furthermore, common estimators for \xibar_3 and \xibar_2 suffer also from this kind of estimation-bias themselves: for \xibar_2, it is equivalent to what is commonly known as the integral constraint. We present a unifying treatment allowing all these estimation-biases to be calculated analytically. They are in general negative, and decrease in significance as the survey volume increases, for a given smoothing scale. We present a re-analysis of some existing measurements of the variance and skewness and show that most of the well-known systematic discrepancies between surveys with similar selection criteria, but different sizes, can be attributed to the volume-dependent estimation-biases. This affects the inference of the galaxy-bias(es) from these surveys. Our methodology can be adapted to measurements of analogous quantities in quasar spectra and weak-lensing maps. We suggest methods to reduce the above estimation-biases, and point out other examples in LSS studies which might suffer from the same type of a nonlinear-estimation-bias.Comment: 28 pages of text, 9 ps figures, submitted to Ap

Clustering of the Diffuse Infrared Light from the COBE DIRBE maps. III. Power spectrum analysis and excess isotropic component of fluctuations

The cosmic infrared background (CIB) radiation is the cosmic repository for energy release throughout the history of the universe. Using the all-sky data from the COBE DIRBE instrument at wavelengths 1.25 - 100 mic we attempt to measure the CIB fluctuations. In the near-IR, foreground emission is dominated by small scale structure due to stars in the Galaxy. There we find a strong correlation between the amplitude of the fluctuations and Galactic latitude after removing bright foreground stars. Using data outside the Galactic plane ($|b| > 20\deg$) and away from the center ($90\deg< l <270\deg$) we extrapolate the amplitude of the fluctuations to cosec$|b|=0$. We find a positive intercept of $\delta F_{\rm rms} = 15.5^{+3.7}_{-7.0},5.9^{+1.6}_{-3.7}, 2.4^{+0.5}_{-0.9}, 2.0^{+0.25}_{-0.5}$ nW/m2/sr at 1.25, 2.2,3.5 and 4.9 mic respectively, where the errors are the range of 92% confidence limits. For color subtracted maps between band 1 and 2 we find the isotropic part of the fluctuations at $7.6^{+1.2}_{-2.4}$ nW/m2/sr. Based on detailed numerical and analytic models, this residual is not likely to originate from the Galaxy, our clipping algorithm, or instrumental noise. We demonstrate that the residuals from the fit used in the extrapolation are distributed isotropically and suggest that this extra variance may result from structure in the CIB. For 2\deg< \theta < 15^\deg, a power-spectrum analysis yields firm upper limits of (\theta/5^\deg) \times\delta F_{\rm rms} (\theta) < 6, 2.5, 0.8, 0.5 nW/m2/sr at 1.25, 2.2, 3.5 and 4.9 mic respectively. From 10-100 mic, the upper limits <1 nW/m2/sr.Comment: Ap.J., in press. 69 pages including 24 fig

Extending the halo mass resolution of NN-body simulations

We present a scheme to extend the halo mass resolution of N-body simulations of the hierarchical clustering of dark matter. The method uses the density field of the simulation to predict the number of sub-resolution dark matter haloes expected in different regions. The technique requires as input the abundance of haloes of a given mass and their average clustering, as expressed through the linear and higher order bias factors. These quantities can be computed analytically or, more accurately, derived from a higher resolution simulation as done here. Our method can recover the abundance and clustering in real- and redshift-space of haloes with mass below $\sim 7.5 \times 10^{13}h^{-1}M_{\odot}$ at $z=0$ to better than 10%. We demonstrate the technique by applying it to an ensemble of 50 low resolution, large-volume $N$-body simulations to compute the correlation function and covariance matrix of luminous red galaxies (LRGs). The limited resolution of the original simulations results in them resolving just two thirds of the LRG population. We extend the resolution of the simulations by a factor of 30 in halo mass in order to recover all LRGs. With existing simulations it is possible to generate a halo catalogue equivalent to that which would be obtained from a $N$-body simulation using more than 20 trillion particles; a direct simulation of this size is likely to remain unachievable for many years. Using our method it is now feasible to build the large numbers of high-resolution large volume mock galaxy catalogues required to compute the covariance matrices necessary to analyse upcoming galaxy surveys designed to probe dark energy.Comment: 11 pages, 7 Figure

Gravity's smoking gun?

We present a new constraint on the biased galaxy formation picture. Gravitational instability theory predicts that the two-point mass density correlation function, \xi(r), has an inflection point at the separation r=r_0, corresponding to the boundary between the linear and nonlinear regime of clustering, \xi = 1. We show how this feature can be used to constrain the square of the biasing parameter, b^2 = \xi_g / \xi on scales r = r_0, where \xi_g is the galaxy-galaxy correlation function, allowed to differ from \xi. We apply our method to real data: the \xi_g(r), estimated from the APM galaxy survey. Our results suggest that the APM galaxies trace the mass at separations r > 5 Mpc/h, where h is the Hubble constant in units of 100 km/s Mpc. The present results agree with earlier studies, based on comparing higher order correlations in the APM with weakly non-linear perturbation theory. Both approaches constrain the "b" factor to be within 20% of unity. If the existence of the feature we identified in the APM \xi_g(r) -- the inflection point near \xi_g = 1 -- is confirmed by more accurate surveys, we may have discovered gravity's smoking gun: the long awaited ``shoulder'' in \xi, predicted by Gott and Rees 25 years ago.Comment: 4 pages, 2 figures, minor changes and references added, matches version published in ApJ letter