Comprehensive Two-Point Analyses of Weak Gravitational Lensing Surveys

We present a framework for analyzing weak gravitational lensing survey data, including lensing and source-density observables, plus spectroscopic redshift calibration data. All two-point observables are predicted in terms of parameters of a perturbed Robertson-Walker metric, making the framework independent of the models for gravity, dark energy, or galaxy properties. For Gaussian fluctuations the 2-point model determines the survey likelihood function and allows Fisher-matrix forecasting. The framework includes nuisance terms for the major systematic errors: shear measurement errors, magnification bias and redshift calibration errors, intrinsic galaxy alignments, and inaccurate theoretical predictions. We propose flexible parameterizations of the many nuisance parameters related to galaxy bias and intrinsic alignment. For the first time we can integrate many different observables and systematic errors into a single analysis. As a first application of this framework, we demonstrate that: uncertainties in power-spectrum theory cause very minor degradation to cosmological information content; nearly all useful information (excepting baryon oscillations) is extracted with ~3 bins per decade of angular scale; and the rate at which galaxy bias varies with redshift substantially influences the strength of cosmological inference. The framework will permit careful study of the interplay between numerous observables, systematic errors, and spectroscopic calibration data for large weak-lensing surveys.Comment: submitted to Ap

The Local Nanohertz Gravitational-Wave Landscape From Supermassive Black Hole Binaries

Supermassive black hole binaries (SMBHBs) in the 10 million to 10 billion $M_\odot$ range form in galaxy mergers, and live in galactic nuclei with large and poorly constrained concentrations of gas and stars. There are currently no observations of merging SMBHBs--- it is in fact possible that they stall at their final parsec of separation and never merge. While LIGO has detected high frequency GWs, SMBHBs emit GWs in the nanohertz to millihertz band. This is inaccessible to ground-based interferometers, but possible with Pulsar Timing Arrays (PTAs). Using data from local galaxies in the 2 Micron All-Sky Survey, together with galaxy merger rates from Illustris, we find that there are on average $91\pm7$ sources emitting GWs in the PTA band, and $7\pm2$ binaries which will never merge. Local unresolved SMBHBs can contribute to GW background anisotropy at a level of $\sim20\%$, and if the GW background can be successfully isolated, GWs from at least one local SMBHB can be detected in 10 years.Comment: submitted to Nature Astronomy (reformatted for arXiv

Large scale anisotropies on halo infall

We perform a statistical analysis of the peculiar velocity field around dark matter haloes in numerical simulations. We examine different properties of the infall of material onto haloes and its relation to central halo shapes and the shape of the large scale surrounding regions (LSSR). We find that the amplitude of the infall velocity field along the halo shape minor axis is larger than that along the major axis. This is consistent for general triaxial haloes, and for both prolate and oblate systems. We also report a strong anisotropy of the velocity field along the principal axes of the LSSR. The infall velocity field around dark matter haloes reaches a maximum value along the direction of the minor axis of the LSSR, whereas along the direction of its major axis, it exhibits the smallest velocities. We also analyse the dependence of the matter velocity field on the local environment. The amplitude of the infall velocity at high local density regions is larger than at low local density regions. The velocity field tends to be more laminar along the direction towards the minor axis of the LSSR, where the mean ratio between flow velocity and velocity dispersion is of order unity and nearly constant up to scales of 15 Mpc/h. We also detect anisotropies in the outflowing component of the velocity field, showing a maximum amplitude along the surrounding LSSR major axis.Comment: 12 pages, 9 figures, accepted for publication in Monthly Notices of the Royal Astronomical Societ

Catastrophic photometric redshift errors: weak lensing survey requirements

We study the sensitivity of weak lensing surveys to the effects of catastrophic redshift errors - cases where the true redshift is misestimated by a significant amount. To compute the biases in cosmological parameters, we adopt an efficient linearized analysis where the redshift errors are directly related to shifts in the weak lensing convergence power spectra. We estimate the number Nspec of unbiased spectroscopic redshifts needed to determine the catastrophic error rate well enough that biases in cosmological parameters are below statistical errors of weak lensing tomography. While the straightforward estimate of Nspec is ~10^6 we find that using only the photometric redshifts with z<=2.5 leads to a drastic reduction in Nspec to ~30,000 while negligibly increasing statistical errors in dark energy parameters. Therefore, the size of spectroscopic survey needed to control catastrophic errors is similar to that previously deemed necessary to constrain the core of the z_s-z_p distribution. We also study the efficacy of the recent proposal to measure redshift errors by cross-correlation between the photo-z and spectroscopic samples. We find that this method requires ~10% a priori knowledge of the bias and stochasticity of the outlier population, and is also easily confounded by lensing magnification bias. The cross-correlation method is therefore unlikely to supplant the need for a complete spectroscopic redshift survey of the source population.Comment: 14 pages, 3 figure

Neutrino mass bounds from cosmology

Cosmology is at present one of the most powerful probes of neutrino properties. The advent of precision data from the cosmic microwave background and large scale structure has allowed for a very strong bound on the neutrino mass. Here, I review the status of cosmological bounds on neutrino properties with emphasis on mass bounds on light neutrinos.Comment: 6 pages, 3 figures, talk given at NOW2004 workshop, Conca Specchiulla, Italy, September 11-17, 200

Average Heating Rate of Hot Atmospheres in Distant Clusters by Radio AGN: Evidence for Continuous AGN Heating

We examine atmospheric heating by radio active galactic nuclei (AGN) in distant X-ray clusters by cross correlating clusters selected from the 400 Square Degree (400SD) X-ray Cluster survey with radio sources in the NRAO VLA Sky Survey. Roughly 30% of the clusters show radio emission above a flux threshold of 3 mJy within a projected radius of 250 kpc. The radio emission is presumably associated with the brightest cluster galaxy. The mechanical jet power for each radio source was determined using scaling relations between radio power and cavity (mechanical) power determined for nearby clusters, groups, and galaxies with hot atmospheres containing X-ray cavities. The average jet power of the central radio AGN is approximately $2\times 10^{44}$\ergs. We find no significant correlation between radio power, hence mechanical jet power, and the X-ray luminosities of clusters in the redshift range 0.1 -- 0.6. This implies that the mechanical heating rate per particle is higher in lower mass, lower X-ray luminosity clusters. The jet power averaged over the sample corresponds to an atmospheric heating of approximately 0.2 keV per particle within R$_{500}$. Assuming the current AGN heating rate does not evolve but remains constant to redshifts of 2, the heating rate per particle would rise by a factor of two. We find that the energy injected from radio AGN contribute substantially to the excess entropy in hot atmospheres needed to break self-similarity in cluster scaling relations. The detection frequency of radio AGN is inconsistent with the presence of strong cooling flows in 400SD clusters, but does not exclude weak cooling flows. It is unclear whether central AGN in 400SD clusters are maintained by feedback at the base of a cooling flow. Atmospheric heating by radio AGN may retard the development of strong cooling flows at early epochs.Comment: ApJ in pres

The Lyman-alpha forest at redshifts 0.1 -- 1.6: good agreement between a large hydrodynamic simulation and HST spectra

We give a comprehensive statistical description of the Lyman-alpha absorption from the intergalactic medium in a hydrodynamic simulation at redshifts 0.1-1.6, the range of redshifts covered by HST spectra of QSOs. We use the ENZO code to make a 76 comoving Mpc cube simulation using 75 kpc cells, for a Hubble constant of 71 km/s/Mpc. The best prior work, by \citet{dave99},used an SPH simulation in a 15.6 Mpc box with an effective resolution of 245 kpc and slightly different cosmological parameters. At redshifts z=2 this simulation is different from data. \citet{tytler07b} found that the simulated spectra at z=2 have too little power on large scales, Lyman-alpha lines are too wide, there is a lack high column density lines, and there is a lack of pixels with low flux. Here we present statistics at z<1.6, including the flux distribution, the mean flux, the effective opacity, and the power and correlation of the flux. We also give statistics of the lyman alpha lines including the line width distribution, the column density distribution, the number of lines per unit equivalent width and redshift, and the correlation between the line width and column density. We find that the mean amount of absorption in the simulated spectra changes smoothly with redshift with DA(z)=0.01(1+z)^{2.25}. Both the trend and absolute values are close to measurements of HST spectra by \citet{kirkman07a}. The column density and line width distributions are also close to those measured from HST spectra by \citet{janknecht06a}, except for the mode of the line width distribution which is smaller in the HST spectra. Although some differences that we saw at z=2 are too subtle to be seen in existing HST spectra, overall, the simulation gives an good description of HST spectra at 0.1<z<1.6

Measuring neutrino masses with a future galaxy survey

We perform a detailed forecast on how well a Euclid-like photometric galaxy and cosmic shear survey will be able to constrain the absolute neutrino mass scale. Adopting conservative assumptions about the survey specifications and assuming complete ignorance of the galaxy bias, we estimate that the minimum mass sum of sum m_nu ~ 0.06 eV in the normal hierarchy can be detected at 1.5 sigma to 2.5 sigma significance, depending on the model complexity, using a combination of galaxy and cosmic shear power spectrum measurements in conjunction with CMB temperature and polarisation observations from Planck. With better knowledge of the galaxy bias, the significance of the detection could potentially reach 5.4 sigma. Interestingly, neither Planck+shear nor Planck+galaxy alone can achieve this level of sensitivity; it is the combined effect of galaxy and cosmic shear power spectrum measurements that breaks the persistent degeneracies between the neutrino mass, the physical matter density, and the Hubble parameter. Notwithstanding this remarkable sensitivity to sum m_nu, Euclid-like shear and galaxy data will not be sensitive to the exact mass spectrum of the neutrino sector; no significant bias (< 1 sigma) in the parameter estimation is induced by fitting inaccurate models of the neutrino mass splittings to the mock data, nor does the goodness-of-fit of these models suffer any significant degradation relative to the true one (Delta chi_eff ^2< 1).Comment: v1: 29 pages, 10 figures. v2: 33 pages, 12 figures; added sections on shape evolution and constraints in more complex models, accepted for publication in JCA

The removal of shear-ellipticity correlations from the cosmic shear signal: Influence of photometric redshift errors on the nulling technique

Cosmic shear is regarded one of the most powerful probes to reveal the properties of dark matter and dark energy. To fully utilize its potential, one has to be able to control systematic effects down to below the level of the statistical parameter errors. Particularly worrisome in this respect is intrinsic alignment, causing considerable parameter biases via correlations between the intrinsic ellipticities of galaxies and the gravitational shear, which mimic lensing. In an earlier work we have proposed a nulling technique that downweights this systematic, only making use of its well-known redshift dependence. We assess the practicability of nulling, given realistic conditions on photometric redshift information. For several simplified intrinsic alignment models and a wide range of photometric redshift characteristics we calculate an average bias before and after nulling. Modifications of the technique are introduced to optimize the bias removal and minimize the information loss by nulling. We demonstrate that one of the presented versions is close to optimal in terms of bias removal, given high quality of photometric redshifts. For excellent photometric redshift information, i.e. at least 10 bins with a small dispersion, a negligible fraction of catastrophic outliers, and precise knowledge about the redshift distributions, one version of nulling is capable of reducing the shear-intrinsic ellipticity contamination by at least a factor of 100. Alternatively, we describe a robust nulling variant which suppresses the systematic signal by about 10 for a very broad range of photometric redshift configurations. Irrespective of the photometric redshift quality, a loss of statistical power is inherent to nulling, which amounts to a decrease of the order 50% in terms of our figure of merit.Comment: 26 pages, including 16 figures; minor changes to match accepted version; published in Astronomy and Astrophysic

Testing Cold Dark Matter Models At Moderate to High Redshift

The COBE microwave background temperature fluctuations and the abundance of local rich clusters of galaxies provide the two most powerful constraints on cosmological models. When all variants of the standard cold dark matter (CDM) model are subject to the combined constraint, the power spectrum of any model is fixed to $\sim 10$ accuracy in both the shape and overall amplitude. These constrained models are not expected to differ dramatically in their local large-scale structure properties. However, their evolutionary histories differ, resulting in dramatic differences towards high redshift. We examine in detail six standardized, COBE and cluster normalized CDM models with respect to a large set of independent observations. The observations include correlation function of rich clusters of galaxies, galaxy power spectrum, evolution of rich cluster abundance, gravitational lensing by moderate -to-high redshift clusters, \lya forest, damped \lya systems, high redshift galaxies, reionization of the universe and future CMB experiments. It seems that each of the independent observations examined is or potentially is capable of distinguishing between at least some of the models. The combined power of several or all of these observations is tremendous. Thus, we appear to be on the verge of being able to make dramatic tests of all models in the near future using a rapidly growing set of observations, mostly at moderate to high redshift. Consistency or inconsistency between different observed phenomena on different scales and/or at different epochs with respect to the models will have profound implications for theory of growth of cosmic structure.Comment: ApJ in press (1998), 26 emulateapj page