Source author record

Bhuvnesh Jain

Bhuvnesh Jain 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

59works
10topics
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

59 published item(s)

preprint2026arXiv

Tidal alignment and tidal torquing modeling for the cosmic shear three-point correlation function and mass aperture skewness

We present a model for the intrinsic alignment contamination of the shear three-point correlation function and skewness of the mass aperture statistic using the tidal alignment and tidal torquing (TATT) formalism. We compute the intrinsic alignment bispectra components in terms of the TATT model parameters. We consider two effective field theory approaches in the literature, relate them to the TATT model parameters and an extension to TATT that includes the velocity-shear (VS) parameter. We compare the impact of changing between NLA, TATT, and TATT+VS on the theoretical computation of the 3PCF using the best fit parameters and tomographic redshift distributions from Dark Energy Survey Year 3. We find that the TATT model significantly impacts the skewed triangle configurations of the 3PCF. Additionally, including the higher-order effects from TATT can introduce opposite effects on the two-point function and on the mass aperture skewness, damping the signal of the former while boosting the signal of the latter. We argue that a joint 2PCF+3PCF analysis with the TATT model can help break the degeneracy between its model parameters and provide more robust constraints on both cosmology and intrinsic alignment amplitude parameters. We show that typical values of order unity for the intrinsic alignment parameters introduce differences of around $10\%$ between NLA and TATT predictions.

preprint2022arXiv

Cosmology with the Roman Space Telescope -- Synergies with CMB lensing

We explore synergies between the Nancy Grace Roman Space Telescope and CMB lensing data to constrain dark energy and modified gravity scenarios. A simulated likelihood analysis of the galaxy clustering and weak lensing data from the Roman Space Telescope High Latitude Survey combined with CMB lensing data from the Simons Observatory is undertaken, marginalizing over important astrophysical effects and calibration uncertainties. Included in the modeling are the effects of baryons on small-scale clustering, scale-dependent growth suppression by neutrinos, as well as uncertainties in the galaxy clustering biases, in the intrinsic alignment contributions to the lensing signal, in the redshift distributions, and in the galaxy shape calibration. The addition of CMB lensing roughly doubles the dark energy figure-of-merit from Roman photometric survey data alone, varying from a factor of 1.7 to 2.4 improvement depending on the particular Roman survey configuration. Alternatively, the inclusion of CMB lensing information can compensate for uncertainties in the Roman galaxy shape calibration if it falls below the design goals. Furthermore, we report the first forecast of Roman constraints on a model-independent structure growth, parameterized by $σ_8 (z)$, and on the Hu-Sawicki f(R) gravity as well as an improved forecast of the phenomenological $(Σ_0,μ_0)$ model. We find that CMB lensing plays a crucial role in constraining $σ_8(z)$ at z>2, with percent-level constraints forecasted out to z=4. CMB lensing information does not improve constraints on the f(R) models substantially. It does, however, increase the $(Σ_0,μ_0)$ figure-of-merit by a factor of about 1.5.

preprint2022arXiv

Perturbation theory models for LSST-era galaxy clustering: tests with sub-percent mock catalog measurements in Fourier and configuration space

We analyze the clustering of galaxies using the z=1.006 snapshot of the CosmoDC2 simulation, a high-fidelity synthetic galaxy catalog designed to validate analysis methods for the Large Synoptic Survey Telescope (LSST). We present sub-percent measurements of the galaxy auto-correlation and galaxy-dark matter cross correlations in Fourier space and configuration space for a magnitude-limited galaxy sample. At these levels of precision, the statistical errors of the measurement are comparable to the systematic effects present in the simulation and measurement procedure; nevertheless, using a hybrid-PT model, we are able to model non-linear galaxy bias with 0.5% precision up to scales of $k_{\rm max}= 0.5 \ h/{\rm Mpc}$ and $r_{\rm min}= 4 \ {\rm Mpc}/h$. While the linear bias parameter is measured with 0.01% precision, other bias parameters are determined with considerably weaker constraints and sometimes bimodal posterior distributions. We compare our fiducial model with lower dimensional models where the higher-order bias parameters are fixed at their co-evolution values and find that leaving these parameters free provides significant improvements in our ability to model small scale information. We also compare bias parameters for galaxy samples defined using different magnitude bands and find agreement between samples containing equal numbers of galaxies. Finally, we compare bias parameters between Fourier space and configuration space and find moderate tension between the two approaches. Although our model is often unable to fit the CosmoDC2 galaxy samples within the 0.1% precision of our measurements, our results suggest that the hybrid-PT model used in this analysis is capable of modeling non-linear galaxy bias within the percent level precision needed for upcoming galaxy surveys.

preprint2022arXiv

Snowmass2021 CMB-HD White Paper

CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10 h Mpc^(-1)), which probes dark matter particle properties. It will also allow 2.) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and velocity, another probe of cosmic structure. In addition, CMB-HD would allow us to cross critical thresholds: 3.) ruling out or detecting any new, light (< 0.1 eV) particles that were in thermal equilibrium with known particles in the early Universe, 4.) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe, and 5.) ruling out or detecting inflationary magnetic fields. CMB-HD would also provide world-leading constraints on 6.) axion-like particles, 7.) cosmic birefringence, 8.) the sum of the neutrino masses, and 9.) the dark energy equation of state. The CMB-HD survey would be delivered in 7.5 years of observing 20,000 square degrees of sky, using two new 30-meter-class off-axis crossed Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

preprint2021arXiv

A 5% measurement of the gravitational constant in the Large Magellanic Cloud

We perform a novel test of General Relativity by measuring the gravitational constant in the Large Magellanic Cloud (LMC). The LMC contains six well-studied Cepheid variable stars in detached eclipsing binaries. Radial velocity and photometric observations enable a complete orbital solution, and precise measurements of the Cepheids' periods permit detailed stellar modelling. Both are sensitive to the strength of gravity, the former via Kepler's third law and the latter through the gravitational free-fall time. We jointly fit the observables for stellar parameters and the gravitational constant. Performing a full Markov Chain Monte Carlo analysis of the parameter space including all relevant nuisance parameters, we constrain the gravitational constant in the Large Magellanic Cloud relative to the Solar System to be $G_\text{LMC}/G_\text{SS} = 0.93^{+0.05}_{-0.04}$. We discuss the implications of this 5% measurement of Newton's constant in another galaxy for dark energy and modified gravity theories. This result excludes one Cepheid, CEP-1812, which is an outlier and needs further study: it is either a highly unusual system to which our model does not apply, or it prefers $G_\text{LMC}<G_\text{SS}$ at $2.6σ$. We also obtain new bounds on critical parameters that appear in semi-analytic descriptions of stellar processes. In particular, we measure the mixing length parameter to be $α=0.90^{+0.36}_{-0.26}$ (when assumed to be constant across our sample), and obtain constraints on the parameters describing turbulent dissipation and convective flux.

preprint2021arXiv

Chameleon Early Dark Energy and the Hubble Tension

Early dark energy (EDE) offers a particularly interesting theoretical approach to the Hubble tension, albeit one that introduces its own set of challenges, including a new `why then' problem related to the EDE injection time at matter-radiation equality, and a mild worsening of the large-scale structure (LSS) tension. Both these challenges center on the properties of dark matter, which becomes the dominant component of the Universe at EDE injection and is also responsible for seeding LSS. Motivated by this, we present the potential of couplings between EDE and dark matter to address these challenges, focusing on a mechanism similar to chameleon dark energy theories, deeming this chameleon early dark energy (CEDE). We present relevant background and perturbation equations and study the dynamics of the case of a quartic scalar potential and an exponential coupling.

preprint2021arXiv

Connecting galaxy evolution in clusters with their radial profiles and phase space distribution: results from the IllustrisTNG hydrodynamical simulations

We study the population of galaxies around galaxy clusters in the hydrodynamic simulation suite IllustrisTNG 300-1 to study the signatures of their evolutionary history on observable properties. We measure the radial number density profile, phase space distribution, and splashback radius for galaxies of different masses and colors over the redshift range $z=0-1$. The three primary physical effects which shape the galaxy distribution within clusters are the galaxy quenching, angular momentum distribution and dynamical friction. We find three distinct populations of galaxies by applying a Gaussian mixture model to their distribution in color and mass. They have distinct evolutionary histories and leave distinct signatures on their distribution around cluster halos. We find that low-mass red galaxies show the most concentrated distribution in clusters and the largest splashback radius, while high-mass red galaxies show a less concentrated distribution and a smaller splashback radius. Blue galaxies, which mostly quench into the low-mass red population, have the shallowest distribution within the clusters, with those on radial orbits quenched rapidly before reaching pericenter. Comparison with the distribution of galaxies from the Dark Energy Survey (DES) survey around Sunyaev-Zeldovich (SZ) clusters from the Atacama Cosmology Telescope (ACT) and South Pole Telescope (SPT) surveys shows evidence for differences in galaxy evolution between simulations and data.

preprint2021arXiv

The Novel Probes Project -- Tests of Gravity on Astrophysical Scales

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments.

preprint2021arXiv

What does a cosmological experiment really measure? Covariant posterior decomposition with normalizing flows

We present methods to rigorously extract parameter combinations that are constrained by data from posterior distributions. The standard approach uses linear methods that apply to Gaussian distributions. We show the limitations of the linear methods for current surveys, and develop non-linear methods that can be used with non-Gaussian distributions, and are independent of the parameter basis. These are made possible by the use of machine-learning models, normalizing flows, to learn posterior distributions from their samples. These models allow us to obtain the local covariance of the posterior at all positions in parameter space and use its inverse, the Fisher matrix, as a local metric over parameter space. The posterior distribution can then be non-linearly decomposed into the leading constrained parameter combinations via parallel transport in the metric space. We test our methods on two non-Gaussian, benchmark examples, and then apply them to the parameter posteriors of the Dark Energy Survey and Planck CMB lensing. We illustrate how our method automatically learns the survey-specific, best constrained effective amplitude parameter $S_8$ for cosmic shear alone, cosmic shear and galaxy clustering, and CMB lensing. We also identify constrained parameter combinations in the full parameter space, and as an application we estimate the Hubble constant, $H_0$, from large-structure data alone.

preprint2020arXiv

A local resolution of the Hubble tension: The impact of screened fifth forces on the cosmic distance ladder

The discrepancy between the values of the Hubble constant $H_0$ derived from the local distance ladder and the cosmic microwave background provides a tantalising hint of new physics. We explore a potential resolution involving screened fifth forces in the local Universe, which alter the Cepheid calibration of supernova distances. In particular, if the Cepheids with direct distance measurements from parallax or water masers are screened but a significant fraction of those in other galaxies are not, neglecting the difference between their underlying period--luminosity relations biases the local $H_0$ measurement high. This difference derives from a reduction in the Cepheid pulsation period and possible increase in luminosity under a fifth force. We quantify the internal and environmental gravitational properties of the Riess et al. distance ladder galaxies to assess their degrees of screening under a range of phenomenological models, and propagate this information into the $H_0$ posterior as a function of fifth force strength. We consider well-studied screening models in scalar--tensor gravity theories such as chameleon, K-mouflage and Vainshtein, along with a recently-proposed mechanism based on baryon--dark matter interactions in which screening is governed by local dark matter density. We find that a fifth force strength $\sim5-30\%$ that of gravity can alleviate (though not resolve) the $H_0$ tension in some scenarios, around the sensitivity level at which tests based on other distance ladder data can constrain this strength. Although our analysis is exploratory and based on screening models not necessarily realised in full theories, our results demonstrate that new physics-based local resolutions of the $H_0$ tension are possible, supplementing those already known in the pre-recombination era.

preprint2020arXiv

CMB-HD: Astro2020 RFI Response

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

preprint2020arXiv

Screened Fifth Forces Mediated by Dark Matter-Baryon Interactions: Theory and Astrophysical Probes

We derive the details of a new screening mechanism where the interactions of baryons and dark matter can be screened according to the local dark matter density. In this mechanism, the value of Newton's constant is dark matter density-dependent, allowing for the possibility that astrophysical phenomena are very different in galaxies less dense than the Milky Way. The parameterized post Newtonian parameter $γ$, which quantifies the difference between kinematical and lensing probes, also depends on dark matter density. We calculate the effects of varying $G$ on various stages of stellar evolution, focusing on observables that impact cosmology: the Cepheid period--luminosity relation and the supernova Ia magnitude--redshift relation. Other potential tests of the model are also investigated including main-sequence, post-main sequence, and low mass dwarf stars. Finally, we discuss how extragalactic tests of $γ$ could provide complementary constraints.

preprint2020arXiv

Statistics of the Chemical Composition of Solar Analog Stars and Links to Planet Formation

The Sun has been found to be depleted in refractory (rock-forming) elements relative to nearby solar analogs, suggesting a potential indicator of planet formation. Given the small amplitude of the depletion, previous analyses have primarily relied on high signal-to-noise stellar spectra and a strictly differential approach to determine elemental abundances. We present an alternative, likelihood-based approach that can be applied to much larger samples of stars with lower precision abundance determinations. We utilize measurements of about 1700 solar analogs from the Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) and the stellar parameter and chemical abundance pipeline (ASPCAP DR16). By developing a hierarchical mixture model for the data, we place constraints on the statistical properties of the elemental abundances, including correlations with condensation temperature and the fraction of stars with refractory element depletions. We find evidence for two distinct populations: a depleted population of stars that makes up the majority of solar analogs including the Sun, and a not-depleted population that makes up between 10-30% of our sample. We find correlations with condensation temperature generally in agreement with higher precision surveys of a smaller sample of stars. Such trends, if robustly linked to the formation of planetary systems, provide a means to connect stellar chemical abundance patterns to planetary systems over large samples of Milky Way stars.

preprint2020arXiv

The Statistics of Extended Debris Disks Measured with Gaia and Planck

Thermal emission from debris disks around stars has been measured using targeted and resolved observations. We present an alternative, likelihood-based approach in which temperature maps from the Planck CMB survey at 857 and 545 GHz are analyzed in conjunction with stellar positions from Gaia to estimate the fraction of stars hosting disks and the thermal emission from the disks. The debris disks are not resolved (or even necessarily detected individually) but their statistical properties and the correlations with stellar properties are measured for several thousand stars. We compare our findings with higher sensitivity surveys of smaller samples of stars. For dimmer stars, in particular K and M-dwarfs, we find about 10 percent of stars within 80 pc have emission consistent with debris disks. We also report on 80 candidate disks, the majority of which are not previously identified. We have previously constrained the properties of Exo-Oort clouds using Planck data -- with future CMB surveys both components can be measured for different stellar types, providing a new avenue to study the outer parts of planetary systems.

preprint2019arXiv

A model independent comparison of supernova and strong lensing cosmography: implications for the Hubble constant tension

We use supernovae measurements, calibrated by the local determination of the Hubble constant $H_0$ by SH0ES, to interpolate the distance-redshift relation using Gaussian process regression. We then predict, independent of the cosmological model, the distances that are measured with strong lensing time delays. We find excellent agreement between these predictions and the measurements. The agreement holds when we consider only the redshift dependence of the distance-redshift relation, independent of the value of $H_0$. Our results disfavor the possibility that lens mass modeling contributes a 10\% bias or uncertainty in the strong lensing analysis, as suggested recently in the literature. In general our analysis strengthens the case that residual systematic errors in both measurements are below the level of the current discrepancy with the CMB determination of $H_0$, and supports the possibility of new physical phenomena on cosmological scales. With additional data our methodology can provide more stringent tests of unaccounted for systematics in the determinations of the distance-redshift relation in the late universe.

preprint2019arXiv

CMB-HD: An Ultra-Deep, High-Resolution Millimeter-Wave Survey Over Half the Sky

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

preprint2018arXiv

Probing Oort clouds around Milky Way stars with CMB surveys

Long-period comets observed in our solar system are believed to originate from the Oort cloud, which is estimated to extend from roughly a few thousand to $10^5$ AU from the Sun. Despite many theoretical arguments for its existence, no direct observations of the cloud have been reported. Here, we explore the possibility of measuring Oort clouds around other stars through their emission at submillimeter wavelengths. Observations with the 545 and 857 GHz bands of the Planck satellite are well matched to the expected temperatures of Oort cloud bodies (on the order of 10 K). By correlating the Planck maps with catalogs of stars observed by the Gaia mission, we are able to constrain interesting regions of the exo-Oort cloud parameter space, placing limits on the total mass and the minimum size of grains in the cloud. We compare our measurements with known debris disk systems -- in the case of Vega and Fomalhaut we find a significant excess that is in agreement with measurements from Herschel. We use the measurements around Fomalhaut to constrain a possible exo-Oort cloud of that system. We explore an observed excess around the brightest and nearest stars in our sample as arising from possible exo-Oort clouds or other extended sources of thermal emission. We argue that future CMB surveys and targeted observations with far-infrared and millimeter wavelength telescopes have the potential to detect exo-Oort clouds or other extended sources of thermal emission beyond $\sim 1000$ AU from the parent stars.

preprint2016arXiv

Clustering and Bias Measurements of SDSS Voids

Using a void catalog from the SDSS survey, we present the first measurements of void clustering and the corresponding void bias. Over the range 30-200 Mpc/h the void auto-correlation is detected at 5-sigma significance for voids of radius 15-20 Mpc/h. We also measure the void-galaxy cross-correlation at higher signal-to-noise and compare the inferred void bias with the autocorrelation results. Void bias is constant with scale for voids of a given size, but its value falls from 5.6 +/- 1.0 to below zero as the void radius increases from 15 to 30 Mpc/h. The comparison of our measurements with carefully matched galaxy mock catalogs, with no free parameters related to the voids, shows that model predictions can be reliably made for void correlations. We study the dependence of void bias on tracer density and void size with a view to future applications. In combination with our previous lensing measurements of void mass profiles, these clustering measurements provide another step towards using voids as cosmological tracers.

preprint2016arXiv

Constraining the Mass-Richness Relationship of redMaPPer Clusters with Angular Clustering

The potential of using cluster clustering for calibrating the mass-observable relation of galaxy clusters has been recognized theoretically for over a decade. Here, we demonstrate the feasibility of this technique to achieve high precision mass calibration using redMaPPer clusters in the Sloan Digital Sky Survey North Galactic Cap. By including cross-correlations between several richness bins in our analysis we significantly improve the statistical precision of our mass constraints. The amplitude of the mass-richness relation is constrained to 7% statistical precision. However, the error budget is systematics dominated, reaching an 18% total error that is dominated by theoretical uncertainty in the bias-mass relation for dark matter halos. We perform a detailed treatment of the effects of assembly bias on our analysis, finding that the contribution of such effects to our parameter uncertainties is somewhat greater than that of measurement noise. We confirm the results from Miyatake et al. (2015) that the clustering amplitude of redMaPPer clusters depends on galaxy concentration, and provide additional evidence in support of this effect being due to some form of assembly bias. The results presented here demonstrate the power of cluster clustering for mass calibration and cosmology provided the current theoretical systematics can be ameliorated.

preprint2016arXiv

Cosmology from large scale galaxy clustering and galaxy-galaxy lensing with Dark Energy Survey Science Verification data

We present cosmological constraints from the Dark Energy Survey (DES) using a combined analysis of angular clustering of red galaxies and their cross-correlation with weak gravitational lensing of background galaxies. We use a 139 square degree contiguous patch of DES data from the Science Verification (SV) period of observations. Using large scale measurements, we constrain the matter density of the Universe as Omega_m = 0.31 +/- 0.09 and the clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/- 0.13 after marginalizing over seven nuisance parameters and three additional cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16} = 0.74 +/- 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 = 0.78 +/- 0.09 using two bins over the range 0.2 <z< 0.5. We study the robustness of the results under changes in the data vectors, modelling and systematics treatment, including photometric redshift and shear calibration uncertainties, and find consistency in the derived cosmological parameters. We show that our results are consistent with previous cosmological analyses from DES and other data sets and conclude with a joint analysis of DES angular clustering and galaxy-galaxy lensing with Planck CMB data, Baryon Accoustic Oscillations and Supernova type Ia measurements.

preprint2016arXiv

Detection of Stacked Filament Lensing Between SDSS Luminous Red Galaxies

We search for the lensing signal of massive filaments between 135,000 pairs of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey. We develop a new estimator that cleanly removes the much larger shear signal of the neighboring LRG halos, relying only on the assumption of spherical symmetry. We consider two models: a thick filament model constructed from ray-tracing simulations for $Λ$CDM model, and a thin filament model which models the filament by a string of halos along the line connecting the two LRGs. We show that the filament lensing signal is in nice agreement with the thick simulation filament, while strongly disfavoring the thin model. The magnitude of the lensing shear due to the filament is below $10^{-4}$. Employing the likelihood ratio test, we find a 4.5$σ$ significance for the detection of the filament lensing signal, corresponding to a null hypothesis fluctuation probability of $3 \times 10^{-6}$. We also carried out several null tests to verify that the residual shear signal from neighboring LRGs and other shear systematics are minimized.

preprint2016arXiv

Lensing Measurements of the Ellipticity of LRG Dark Matter Halos

Lensing measurements of the shapes of dark matter halos can provide tests of gravity theories and possible dark matter interactions. We measure the quadrupole weak lensing signal from the elliptical halos of 70,000 SDSS Luminous Red Galaxies. We use a new estimator that nulls the spherical halo lensing signal, isolating the shear due to anisotropy in the dark matter distribution. One of the two Cartesian components of our estimator is insensitive to the primary systematic, a spurious alignment of lens and source ellipticities, allowing us to make robust measurements of halo ellipticity. Our best-fit value for the ellipticity of the surface mass density is $0.24 \pm 0.06$, which translates to an axis ratio of 0.78. We rule out the hypothesis of no ellipticity at the $4σ$ confidence level, and ellipticity < 0.12 (axis ratio > 0.89) at the $2σ$ level. We discuss how our measurements of halo ellipticity are revised to higher values using estimates of the misalignment of mass and light from simulations. Finally, we apply the same techniques to a smaller sample of redMaPPer galaxy clusters and obtain a $3σ$ measurement of cluster ellipticity. We discuss how the improved signal to noise properties of our estimator can enable studies of halo shapes for different galaxy populations with upcoming surveys.

preprint2016arXiv

Tidal stripping as a test of satellite quenching in redMaPPer clusters

When dark matter halos are accreted by massive host clusters, strong gravitational tidal forces begin stripping mass from the accreted subhalos. This stripping eventually removes all mass beyond a subhalo's tidal radius, but the unbound mass remains in the vicinity of the satellite for at least a dynamical time t_dyn. The N-body subhalo study of Chamberlain et al. verified this picture and pointed out a useful observational consequence: measurements of subhalo correlations beyond the tidal radius are sensitive to the infall time, t_infall, of the subhalo onto its host. We perform this cross-correlation measurement using ~ 160,000 red satellite galaxies in SDSS redMaPPer clusters and find evidence that subhalo correlations do persist well beyond the tidal radius, suggesting that many of the observed satellites fell into their current host less than a dynamical time ago, t_infall < t_dyn. Combined with estimated dynamical times t_dyn ~ 3-5 Gyr and SED fitting results for the time at which satellites stopped forming stars, t_quench ~ 6 Gyr, we infer that for a significant fraction of the satellites, star formation quenched before those satellites entered their current hosts. The result holds for red satellites over a large range of cluster-centric distances 0.1 - 0.6 Mpc/h. We discuss the implications of this result for models of galaxy formation.

preprint2015arXiv

Cosmic Discordance: Are Planck CMB and CFHTLenS weak lensing measurements out of tune?

We examine the level of agreement between low redshift weak lensing data and the CMB using measurements from the CFHTLenS and Planck+WMAP polarization. We perform an independent analysis of the CFHTLenS six bin tomography results of Heymans et al. (2013). We extend their systematics treatment and find the cosmological constraints to be relatively robust to the choice of non-linear modeling, extension to the intrinsic alignment model and inclusion of baryons. We find that the 90% confidence contours of CFHTLenS and Planck+WP do not overlap even in the full 6-dimensional parameter space of $Λ$CDM, so the two datasets are discrepant. Allowing a massive active neutrino or tensor modes does not significantly resolve the disagreement in the full n-dimensional parameter space. Our results differ from some in the literature because we use the full tomographic information in the weak lensing data and marginalize over systematics. We note that adding a sterile neutrino to $Λ$CDM does bring the 8-dimensional 64% contours to overlap, mainly due to the extra effective number of neutrino species, which we find to be 0.84 $\pm$ 0.35 (68%) greater than standard on combining the datasets. We discuss why this is not a completely satisfactory resolution, leaving open the possibility of other new physics or observational systematics as contributing factors. We provide updated cosmology fitting functions for the CFHTLenS constraints and discuss the differences from ones used in the literature.

preprint2014arXiv

Beyond the Cosmological Standard Model

After a decade and a half of research motivated by the accelerating universe, theory and experiment have a reached a certain level of maturity. The development of theoretical models beyond Λ, or smooth dark energy, often called modified gravity, has led to broader insights into a path forward, and a host of observational and experimental tests have been developed. In this review we present the current state of the field and describe a framework for anticipating developments in the next decade. We identify the guiding principles for rigorous and consistent modifications of the standard model, and discuss the prospects for empirical tests. We begin by reviewing attempts to consistently modify Einstein gravity in the infrared, focusing on the notion that additional degrees of freedom introduced by the modification must screen themselves from local tests of gravity. We categorize screening mechanisms into three broad classes: mechanisms which become active in regions of high Newtonian potential, those in which first derivatives become important, and those for which second derivatives are important. Examples of the first class, such as f(R) gravity, employ the familiar chameleon or symmetron mechanisms, whereas examples of the last class are galileon and massive gravity theories, employing the Vainshtein mechanism. In each case, we describe the theories as effective theories. We describe experimental tests, summarizing laboratory and solar system tests and describing in some detail astrophysical and cosmological tests. We discuss future tests which will be sensitive to different signatures of new physics in the gravitational sector. Parts that are more relevant to theorists vs. observers/experimentalists are clearly indicated, in the hope that this will serve as a useful reference for both audiences, as well as helping those interested in bridging the gap between them.

preprint2014arXiv

Testing Gravity Theories Using Stars

Modified theories of gravity have received a renewed interest due to their ability to account for the cosmic acceleration. In order to satisfy the solar system tests of gravity, these theories need to include a screening mechanism that hides the modifications on small scales. One popular and well-studied theory is chameleon gravity. Our own galaxy is necessarily screened, but less dense dwarf galaxies may be unscreened and their constituent stars can exhibit novel features. In particular, unscreened stars are brighter, hotter and more ephemeral than screened stars in our own galaxy. They also pulsate with a shorter period. In this essay, we exploit these new features to constrain chameleon gravity to levels three orders of magnitude lower the previous measurements. These constraints are currently the strongest in the literature.

preprint2013arXiv

Astrophysical Tests of Modified Gravity: the Morphology and Kinematics of Dwarf Galaxies

This paper is the third in a series on tests of gravity using observations of stars and nearby dwarf galaxies. We carry out four distinct tests using published data on the kinematics and morphology of dwarf galaxies, motivated by the theoretical work of Hui et al. (2009) and Jain and Vanderplas (2011). In a wide class of gravity theories a scalar field couples to matter and provides an attractive fifth force. Due to their different self-gravity, stars and gas may respond differently to the scalar force leading to several observable deviations from standard gravity. HI gas, red giant stars and main sequence stars can be displaced relative to each other, and the stellar disk can display warps or asymmetric rotation curves aligned with external potential gradients. To distinguish the effects of modified gravity from standard astrophysical phenomena, we use a control sample of galaxies that are expected to be screened from the fifth force. In all cases we find no significant deviation from the null hypothesis of general relativity. The limits obtained from dwarf galaxies are not yet competitive with the limits from cepheids obtained in our first paper, but can be improved to probe regions of parameter space that are inaccessible using other tests. We discuss how our methodology can be applied to new radio and optical observations of nearby galaxies.

preprint2013arXiv

Growth of Cosmic Structure: Probing Dark Energy Beyond Expansion

The quantity and quality of cosmic structure observations have greatly accelerated in recent years. Further leaps forward will be facilitated by imminent projects, which will enable us to map the evolution of dark and baryonic matter density fluctuations over cosmic history. The way that these fluctuations vary over space and time is sensitive to the nature of dark matter and dark energy. Dark energy and gravity both affect how rapidly structure grows; the greater the acceleration, the more suppressed the growth of structure, while the greater the gravity, the more enhanced the growth. While distance measurements also constrain dark energy, the comparison of growth and distance data tests whether General Relativity describes the laws of physics accurately on large scales. Modified gravity models are able to reproduce the distance measurements but at the cost of altering the growth of structure (these signatures are described in more detail in the accompanying paper on Novel Probes of Gravity and Dark Energy). Upcoming surveys will exploit these differences to determine whether the acceleration of the Universe is due to dark energy or to modified gravity. To realize this potential, both wide field imaging and spectroscopic redshift surveys play crucial roles. Projects including DES, eBOSS, DESI, PFS, LSST, Euclid, and WFIRST are in line to map more than a 1000 cubic-billion-light-year volume of the Universe. These will map the cosmic structure growth rate to 1% in the redshift range 0<z<2, over the last 3/4 of the age of the Universe.

preprint2013arXiv

Information content of weak lensing power spectrum and bispectrum: including the non-Gaussian error covariance matrix

We address the amount of information in the non-Gaussian regime of weak lensing surveys by modelling all relevant covariances of the power spectra and bispectra, using 1000 ray-tracing simulation realizations for a Lambda-CDM model and an analytical halo model. We develop a formalism to describe the covariance matrices of power spectra and bispectra of all triangle configurations. In addition to the known contributions which extend up to six-point correlation functions, we propose a new contribution `the halo sample variance (HSV)' arising from the coupling of the lensing Fourier modes with large-scale mass fluctuations on scales comparable with the survey region via halo bias theory. We show that the model predictions are in good agreement with the simulation once we take the HSV into account. The HSV gives a dominant contribution to the covariance matrices at multipoles l > 10^3, which arises from massive haloes with a mass of > 10^14 solar mass and at relatively low redshifts z < 0.4. Since such haloes are easily identified from a multi-colour imaging survey, the effect can be estimated from the data. By adding the bispectrum to the power spectrum, the total information content or the cumulative signal-to-noise ratio up to a certain maximum multipole of a few 10^3 is improved by 20--50 per cent, which is equivalent to a factor of 1.4--2.3 larger survey area for the power spectrum measurement alone. However, it is still smaller than the case of a Gaussian field by a factor of 3 mostly due to the HSV. Thus bispectrum measurements are useful for cosmology, but using information from upcoming surveys requires that non-Gaussian covariances are carefully estimated.

preprint2013arXiv

Novel Probes of Gravity and Dark Energy

The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure.

preprint2012arXiv

Astrophysical Tests of Modified Gravity: A Screening Map of the Nearby Universe

Astrophysical tests of modified modified gravity theories in the nearby universe have been emphasized recently by Hui, Nicolis and Stubbs (2009) and Jain and VanderPlas (2011). A key element of such tests is the screening mechanism whereby general relativity is restored in massive halos or high density environments like the Milky Way. In chameleon theories of gravity, including all f(R) models, field dwarf galaxies may be unscreened and therefore feel an extra force, as opposed to screened galaxies. The first step to study differences between screened and unscreened galaxies is to create a 3D screening map. We use N-body simulations to test and calibrate simple approximations to determine the level of screening in galaxy catalogs. Sources of systematic errors in the screening map due to observational inaccuracies are modeled and their contamination is estimated. We then apply our methods to create a map out to 200 Mpc in the Sloan Digital Sky Survey footprint using data from the Sloan survey and other sources. In two companion papers this map will be used to carry out new tests of gravity using distance indicators and the disks of dwarf galaxies. We also make our screening map publicly available.

preprint2012arXiv

Astrophysical Tests of Modified Gravity: Constraints from Distance Indicators in the Nearby Universe

We use distance measurements in the nearby universe to carry out new tests of gravity, surpassing other astrophysical tests by over two orders of magnitude for chameleon theories. The three nearby distance indicators -- cepheids, tip of the red giant branch (TRGB) stars, and water masers -- operate in gravitational fields of widely different strengths. This enables tests of scalar-tensor gravity theories because they are screened from enhanced forces to different extents. Inferred distances from cepheids and TRGB stars are altered (in opposite directions) over a range of chameleon gravity theory parameters well below the sensitivity of cosmological probes. Using published data we have compared cepheid and TRGB distances in a sample of unscreened dwarf galaxies within 10 Mpc. As a control sample we use a comparable set of screened galaxies. We find no evidence for the order unity force enhancements expected in these theories. Using a two-parameter description of the models (the coupling strength and background field value) we obtain constraints on chameleon and symmetron screening scenarios. In particular we show that f(R) models with background field values fR0 above 5e-7 are ruled out at the 95% confidence level. We also compare TRGB and maser distances to the galaxy NGC 4258 as a second test for larger field values. While there are several approximations and caveats in our study, our analysis demonstrates the power of gravity tests in the local universe. We discuss the prospects for additional, improved tests with future observations.

preprint2012arXiv

The impact of camera optical alignments on weak lensing measures for the Dark Energy Survey

Telescope Point Spread Function (PSF) quality is critical for realising the potential of cosmic weak lensing observations to constrain dark energy and test General Relativity. In this paper we use quantitative weak gravitational lensing measures to inform the precision of lens optical alignment, with specific reference to the Dark Energy Survey (DES). We compute optics spot diagrams and calculate the shear and flexion of the PSF as a function of position on the focal plane. For perfect optical alignment we verify the high quality of the DES optical design, finding a maximum PSF contribution to the weak lensing shear of 0.04 near the edge of the focal plane. However this can be increased by a factor of approximately three if the lenses are only just aligned within their maximum specified tolerances. We calculate the E and B-mode shear and flexion variance as a function of de-centre or tilt of each lens in turn. We find tilt accuracy to be a few times more important than de-centre, depending on the lens considered. Finally we consider the compound effect of de-centre and tilt of multiple lenses simultaneously, by sampling from a plausible range of values of each parameter. We find that the compound effect can be around twice as detrimental as when considering any one lens alone. Furthermore, this combined effect changes the conclusions about which lens is most important to align accurately. For DES, the tilt of the first two lenses is the most important.

preprint2011arXiv

Designing Surveys for Tests of Gravity

Modified gravity theories may provide an alternative to dark energy to explain cosmic acceleration. We argue that the observational program developed to test dark energy needs to be augmented to capture new tests of gravity on astrophysical scales. Several distinct signatures of gravity theories exist outside the linear regime, especially owing to the screening mechanism that operates inside halos like the Milky Way to ensure that gravity tests in the solar system are satisfied. This opens up several decades in length scale and new classes of galaxies at low-redshift that can be exploited by surveys. While theoretical work on models of gravity is in the early stages, we can already identify new regimes which cosmological surveys could target to test gravity. These include: 1. A small scale component that focuses on the interior and vicinity of galaxy and cluster halos. 2. Spectroscopy of low redshift galaxies, especially galaxies smaller than the Milky Way, in environments that range from voids to clusters. 3. A program of combining lensing and dynamical information, from imaging and spectroscopic surveys respectively, on the same (or statistically identical) sample of galaxies.

preprint2011arXiv

Halo Scale Predictions of Symmetron Modified Gravity

We offer predictions of symmetron modified gravity in the neighborhood of realistic dark matter halos. The predictions for the fifth force are obtained by solving the nonlinear symmetron equation of motion in the spherical NFW approximation. In addition, we compare the three major known screening mechanisms: Vainshtein, Chameleon, and Symmetron around such dark matter sources, emphasizing the significant differences between them and highlighting observational tests which exploit these differences. Finally, we demonstrate the host halo environmental screening effect ("blanket screening") on smaller satellite halos by solving for the modified forces around a density profile which is the sum of satellite and approximate host components.

preprint2011arXiv

Interpolating Masked Weak Lensing Signal with Karhunen-Loeve Analysis

We explore the utility of Karhunen Loeve (KL) analysis in solving practical problems in the analysis of gravitational shear surveys. Shear catalogs from large-field weak lensing surveys will be subject to many systematic limitations, notably incomplete coverage and pixel-level masking due to foreground sources. We develop a method to use two dimensional KL eigenmodes of shear to interpolate noisy shear measurements across masked regions. We explore the results of this method with simulated shear catalogs, using statistics of high-convergence regions in the resulting map. We find that the KL procedure not only minimizes the bias due to masked regions in the field, it also reduces spurious peak counts from shape noise by a factor of ~ 3 in the cosmologically sensitive regime. This indicates that KL reconstructions of masked shear are not only useful for creating robust convergence maps from masked shear catalogs, but also offer promise of improved parameter constraints within studies of shear peak statistics.

preprint2011arXiv

Microlensing of Kepler Stars as a Method of Detecting Primordial Black Hole Dark Matter

If the Dark Matter consists of primordial black holes (PBHs), we show that gravitational lensing of stars being monitored by NASA's Kepler search for extra-solar planets can cause significant numbers of detectable microlensing events. A search through the roughly 150,000 lightcurves would result in large numbers of detectable events for PBHs in the mass range $5 \ten{-10}\msun$ to $\aten{-4}\msun$. Non-detection of these events would close almost two orders of magnitude of the mass window for PBH dark matter. The microlensing rate is higher than previously noticed due to a combination of the exceptional photometric precision of the Kepler mission and the increase in cross section due to the large angular sizes of the relatively nearby Kepler field stars. We also present a new formalism for calculating optical depth and microlensing rates in the presence of large finite-source effects.

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

Tests of Modified Gravity with Dwarf Galaxies

In modified gravity theories that seek to explain cosmic acceleration, dwarf galaxies in low density environments can be subject to enhanced forces. The class of scalar-tensor theories, which includes f(R) gravity, predict such a force enhancement (massive galaxies like the Milky Way can evade it through a screening mechanism that protects the interior of the galaxy from this "fifth" force). We study observable deviations from GR in the disks of late-type dwarf galaxies moving under gravity. The fifth-force acts on the dark matter and HI gas disk, but not on the stellar disk owing to the self-screening of main sequence stars. We find four distinct observable effects in such disk galaxies: 1. A displacement of the stellar disk from the HI disk. 2. Warping of the stellar disk along the direction of the external force. 3. Enhancement of the rotation curve measured from the HI gas compared to that of the stellar disk. 4. Asymmetry in the rotation curve of the stellar disk. We estimate that the spatial effects can be up to 1 kpc and the rotation velocity effects about 10 km/s in infalling dwarf galaxies. Such deviations are measurable: we expect that with a careful analysis of a sample of nearby dwarf galaxies one can improve astrophysical constraints on gravity theories by over three orders of magnitude, and even solar system constraints by one order of magnitude. Thus effective tests of gravity along the lines suggested by Hui et al (2009) and Jain (2011) can be carried out with low-redshift galaxies, though care must be exercised in understanding possible complications from astrophysical effects.

preprint2010arXiv

3D Reconstruction of the Density Field: An SVD Approach to Weak Lensing Tomography

We present a new method for constructing three-dimensional mass maps from gravitational lensing shear data. We solve the lensing inversion problem using truncation of singular values (within the context of generalized least squares estimation) without a priori assumptions about the statistical nature of the signal. This singular value framework allows a quantitative comparison between different filtering methods: we evaluate our method beside the previously explored Wiener filter approaches. Our method yields near-optimal angular resolution of the lensing reconstruction and allows cluster sized halos to be de-blended robustly. It allows for mass reconstructions which are 2-3 orders-of-magnitude faster than the Wiener filter approach; in particular, we estimate that an all-sky reconstruction with arcminute resolution could be performed on a time-scale of hours. We find however that linear, non-parametric reconstructions have a fundamental limitation in the resolution achieved in the redshift direction.

preprint2010arXiv

Cosmological Tests of Gravity

Modifications of general relativity provide an alternative explanation to dark energy for the observed acceleration of the universe. We review recent developments in modified gravity theories, focusing on higher dimensional approaches and chameleon/f(R) theories. We classify these models in terms of the screening mechanisms that enable such theories to approach general relativity on small scales (and thus satisfy solar system constraints). We describe general features of the modified Friedman equation in such theories. The second half of this review describes experimental tests of gravity in light of the new theoretical approaches. We summarize the high precision tests of gravity on laboratory and solar system scales. We describe in some detail tests on astrophysical scales ranging from ~kpc (galaxy scales) to ~Gpc (large-scale structure). These tests rely on the growth and inter-relationship of perturbations in the metric potentials, density and velocity fields which can be measured using gravitational lensing, galaxy cluster abundances, galaxy clustering and the Integrated Sachs-Wolfe effect. A robust way to interpret observations is by constraining effective parameters, such as the ratio of the two metric potentials. Currently tests of gravity on astrophysical scales are in the early stages --- we summarize these tests and discuss the interesting prospects for new tests in the coming decade.

preprint2010arXiv

Lensing Magnification: Implications for Counts of Submillimeter Galaxies and SZ Clusters

We study lensing magnification of source galaxies by intervening galaxy groups and clusters using a halo model. Halos are modeled with truncated NFW profiles with ellipticity added to their lensing potential and propagated to observable lensing statistics. We present the formalism to calculate observable effects due to a distribution of halos of different masses at different redshifts along the l ine of sight. We calculate the effects of magnification on the number counts of high-redshift galaxies. Using BLAST survey data for submillimeter galaxies (SMGs), we find that magnification affects the steep, high flux par t of the counts by about 60%. The effect becomes much stronger if the intrinsic distribution is signi ficantly steeper than observed. We also consider the effect of this high-redshift galaxy population on contaminating the Sunyaev-Zel'dovich (SZ) signal of massive clusters using the halo model approach. We find that for the majority of clusters expected to be detected with ongoing SZ surveys, there is significant contamination from the Poisson noise due to background SMGs. This contr ibution can be comparable to the SZ increment for typical clusters and can also contaminate the SZ decrement of low mass clusters. Thus SZ observations, especially for the increment part of the SZ spectrum, need to include careful modeling of this irreducible contamination for mass estimation. Lensing further enhances the contamination, especially close to the cores of massive clusters and for very disturbed clusters with large magnification cross-section.

preprint2010arXiv

Magnification Effects on Source Counts and Fluxes

We consider the effect of lensing magnification on high redshift sources in the case that magnification varies on the sky, as expected in wide fields of view or within observed galaxy clusters. We give expressions for number counts, flux and flux variance as integrals over the probability distribution of the magnification. We obtain these through a simple mapping between averages over the observed sky and over the magnification probability distribution in the source plane. Our results clarify conflicting expressions in the literature and can be used to calculate a variety of magnification effects. We highlight two applications: 1. Lensing of high-z galaxies by galaxy clusters can provide the dominant source of scatter in SZ observations at frequencies larger than the SZ null. 2. The number counts of high-z galaxies with a Schechter-like luminosity function will be changed at high luminosities to a power law, with significant enhancement of the observed counts at L > 10 L*.

preprint2010arXiv

Re-capturing cosmic information

Gravitational lensing of distant galaxies can be exploited to infer the convergence field as a function of angular position on the sky. The statistics of this field, much like that of the cosmic microwave background (CMB), can be studied to extract information about fundamental parameters in cosmology, most notably the dark energy in the Universe. Unlike the CMB, the distribution of matter in the Universe which determines the convergence field is highly non-Gaussian, reflecting the nonlinear processes which accompanied structure formation. Much of the cosmic information contained in the initial field is therefore unavailable to the standard power spectrum measurements. Here we propose a method for re-capturing cosmic information by using the power spectrum of a simple function of the observed (nonlinear) convergence field. We adapt the approach of Neyrinck et al. (2009) to lensing by using a modified logarithmic transform of the convergence field. The Fourier transform of the log-transformed field has modes that are nearly uncorrelated, which allows for additional cosmological information to be extracted from small-scale modes.

preprint2010arXiv

Submillimeter Galaxy Number Counts and Magnification by Galaxy Clusters

We present an analytical model which reproduces measured galaxy number counts from surveys in the wavelength range of 500 micron to 2 mm. The model involves a single high-redshift galaxy population with a Schechter luminosity function which has been gravitationally lensed by galaxy clusters in the mass range 10^13 to 10^15 Msun. This simple model reproduces both the low flux and the high flux end of the number counts reported by the BLAST, SCUBA, AzTEC and the SPT surveys. In particular, our model accounts for the most luminous galaxies detected by SPT as the result of high magnifications by galaxy clusters (magnification factors of 10-30). This interpretation implies that submillimeter and millimeter surveys of this population may prove to be a useful addition to ongoing cluster detection surveys. The model also implies that the bulk of submillimeter galaxies detected at wavelengths larger than 500 micron lie at redshifts greater than 2.

preprint2009arXiv

Galaxy-CMB and galaxy-galaxy lensing on large scales: sensitivity to primordial non-Gaussianity

A convincing detection of primordial non-Gaussianity in the local form of the bispectrum, whose amplitude is given by the fNL parameter, offers a powerful test of inflation. In this paper we calculate the modification of two-point cross-correlation statistics of weak lensing - galaxy-galaxy lensing and galaxy-Cosmic Microwave Background (CMB) cross-correlation - due to fNL. We derive and calculate the covariance matrix of galaxy-galaxy lensing including cosmic variance terms. We focus on large scales (l<100) for which the shape noise of the shear measurement becomes irrelevant and cosmic variance dominates the error budget. For a modest degree of non-Gaussianity, fNL=+/-50, modifications of the galaxy-galaxy lensing signal at the 10% level are seen on scales R~300 Mpc, and grow rapidly toward larger scales as \propto R^2. We also see a clear signature of the baryonic acoustic oscillation feature in the matter power spectrum at ~150 Mpc, which can be measured by next-generation lensing experiments. In addition we can probe the local-form primordial non-Gaussianity in the galaxy-CMB lensing signal by correlating the lensing potential reconstructed from CMB with high-z galaxies. For example, for fNL=+/-50, we find that the galaxy-CMB lensing cross power spectrum is modified by ~10% at l~40, and by a factor of two at l~10, for a population of galaxies at z=2 with a bias of 2. The effect is greater for more highly biased populations at larger z; thus, high-z galaxy surveys cross-correlated with CMB offer a yet another probe of primordial non-Gaussianity.

preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

preprint2009arXiv

Tests of Gravity from Imaging and Spectroscopic Surveys

Tests of gravity on large-scales in the universe can be made using both imaging and spectroscopic surveys. The former allow for measurements of weak lensing, galaxy clustering and cross-correlations such as the ISW effect. The latter probe galaxy dynamics through redshift space distortions. We use a set of basic observables, namely lensing power spectra, galaxy-lensing and galaxy-velocity cross-spectra in multiple redshift bins (including their covariances), to estimate the ability of upcoming surveys to test gravity theories. We use a two-parameter description of gravity that allows for the Poisson equation and the ratio of metric potentials to depart from general relativity. We find that the combination of imaging and spectroscopic observables is essential in making robust tests of gravity theories. The range of scales and redshifts best probed by upcoming surveys is discussed. We also compare our parametrization to others used in the literature, in particular the gamma parameter modification of the growth factor.

preprint2009arXiv

The Impact of Non-Gaussian Errors on Weak Lensing Surveys

The weak lensing power spectrum carries cosmological information via its dependence on the growth of structure and on geometric factors. Since much of the cosmological information comes from scales affected by nonlinear clustering, measurements of the lensing power spectrum can be degraded by non-Gaussian covariances. Recently there have been conflicting studies about the level of this degradation. We use the halo model to estimate it and include new contributions related to the finite size of lensing surveys, following Rimes and Hamilton's study of 3D simulations. We find that non-Gaussian correlations between different multipoles can degrade the cumulative signal-to-noise for the power spectrum amplitude by up to a factor of 2 (or 5 for a worst-case model that exceeds current N-body simulation predictions). However, using an eight-parameter Fisher analysis we find that the marginalized errors on individual parameters are degraded by less than 10% (or 20% for the worst-case model). The smaller degradation in parameter accuracy is primarily because: individual parameters in a high-dimensional parameter space are degraded much less than the volume of the full Fisher ellipsoid; lensing involves projections along the line of sight, which reduce the non-Gaussian effect; some of the cosmological information comes from geometric factors which are not degraded at all. We contrast our findings with those of Lee & Pen (2008) who suggested a much larger degradation in information content. Finally, our results give a useful guide for exploring survey design by giving the cosmological information returns for varying survey area, depth and the level of some systematic errors.

preprint2008arXiv

Weak Gravitational Lensing and its Cosmological Applications

Weak gravitational lensing is a unique probe of the dark side of the universe: it provides a direct way to map the distribution of dark matter around galaxies, clusters of galaxies and on cosmological scales. Furthermore, the measurement of lensing induced distortions of the shapes of distant galaxies is a powerful probe of dark energy. In this review we describe how lensing measurements are made and interpreted. We discuss various systematic effects that can hamper progress and how they may be overcome. We review some of the recent results in weak lensing by galaxies, galaxy clusters and cosmic shear and discuss the prospects for dark energy measurements from planned surveys.

preprint2007arXiv

Color Tomography

Lensing tomography with multi-color imaging surveys can probe dark energy and the cosmological power spectrum. However accurate photometric redshifts for tomography out to high redshift require imaging in five or more bands, which is expensive to carry out over thousands of square degrees. Since lensing makes coarse, statistical use of redshift information, we explore the prospects for tomography using limited color information from two or three band imaging. With an appropriate calibration sample, we find that it is feasible to create up to four redshift bins using imaging data in just the g, r and i bands. We construct such redshift sub-samples from mock catalogs by clustering galaxies in color space and discarding regions with poorly-defined redshift distributions. The loss of galaxy number density decreases the accuracy of lensing measurements, but even losing half or more of the galaxies is not a severe loss for large area surveys. We estimate the errors on lensing power spectra and dark energy parameters with color tomography and discuss trade-offs in survey area and filter choice. We discuss the systematic errors that may change our conclusions, especially the information needed to tackle intrinsic alignments.

preprint2005arXiv

Cosmological parameters from lensing power spectrum and bispectrum tomography

We examine how lensing tomography with the bispectrum and power spectrum can constrain cosmological parameters and the equation of state of dark energy. Our analysis uses the full information at the two- and three-point level from angular scales of a few degrees to 5 arcminutes (50 < l < 3000), which will be probed by lensing surveys. We use all triangle configurations, cross-power spectra and bispectra constructed from up to three redshift bins with photometric redshifts, and relevant covariances in our analysis. We find that the parameter constraints from bispectrum tomography are comparable to those from power spectrum tomography. Combining the two improves parameter accuracies by a factor of three due to their complementarity. For the dark energy parameterization w(a) = w0 + wa(1-a), the marginalized errors from lensing alone are sigma(w0) = 0.03 fsky^{-1/2} and sigma(wa) = 0.1 fsky^{-1/2}. We show that these constraints can be further improved when combined with measurements of the cosmic microwave background or Type Ia supernovae. The amplitude and shape of the mass power spectrum are also shown to be precisely constrained. We use hyper-extended perturbation theory to compute the nonlinear lensing bispectrum for dark energy models. Accurate model predictions of the bispectrum in the moderately nonlinear regime, calibrated with numerical simulations, will be needed to realize the parameter accuracy we have estimated. Finally, we estimate how well the lensing bispectrum can constrain a model with primordial non-Gaussianity.

preprint2005arXiv

Systematic Errors in Future Weak Lensing Surveys: Requirements and Prospects for Self-Calibration

We study the impact of systematic errors on planned weak lensing surveys and compute the requirements on their contributions so that they are not a dominant source of the cosmological parameter error budget. The generic types of error we consider are multiplicative and additive errors in measurements of shear, as well as photometric redshift errors. In general, more powerful surveys have stronger systematic requirements. For example, for a SNAP-type survey the multiplicative error in shear needs to be smaller than 1%(fsky/0.025)^{-1/2} of the mean shear in any given redshift bin, while the centroids of photometric redshift bins need to be known to better than 0.003(fsky/0.025)^{-1/2}. With about a factor of two degradation in cosmological parameter errors, future surveys can enter a self-calibration regime, where the mean systematic biases are self-consistently determined from the survey and only higher-order moments of the systematics contribute. Interestingly, once the power spectrum measurements are combined with the bispectrum, the self-calibration regime in the variation of the equation of state of dark energy w_a is attained with only a 20-30% error degradation.

preprint2004arXiv

Constraints on Dark Energy Models from Galaxy Clusters with Multiple Arcs

We make an exploratory study of how well dark energy models can be constrained using lensed arcs at different redshifts behind cluster lenses. Arcs trace the critical curves of clusters, and the growth of critical curves with source redshift is sensitive to the dark-energy equation of state. Using analytical models and numerically simulated clusters, we explore the key factors involved in using cluster arcs as a probe of dark energy. We quantify the sensitivity to lens mass, concentration and ellipticity with analytical models that include the effects of dark energy on halo structure. We show with simple examples how degeneracies between mass models and cosmography may be broken using arcs at multiple redshifts or additional constraints on the lens density profile. However we conclude that the requirements on the data are so stringent that it is very unlikely that robust constraints can be obtained from individual clusters. We argue that surveys of clusters, analyzed in conjunction with numerical simulations, are a more promising prospect for arc-cosmography. We use such numerically simulated clusters to estimate how large a sample of clusters/arcs could provide interesting constraints on dark energy models. We focus on the scatter produced by differences in the mass distribution of individual clusters. We find from our sample of simulated clusters that at least 1000 pairs of arcs are needed to obtain constraints if the mass distribution of individual clusters is taken to be undetermined. We discuss several unsolved problems that need study to fully develop this method for precision studies with future surveys.

preprint2004arXiv

Three-point Correlation Functions of SDSS Galaxies in Redshift Space: Morphology, Color, and Luminosity Dependence

We present measurements of the redshift--space three-point correlation function of galaxies in the Sloan Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy morphology, color and luminosity. In order to control systematics due to selection effects, we used $r$--band, volume-limited samples of galaxies, constructed from the magnitude-limited SDSS data ($14.5<r<17.5$), and further divided the samples into two morphological types (early and late) or two color populations (red and blue). The three-point correlation function of SDSS galaxies follow the hierarchical relation well and the reduced three-point amplitudes in redshift--space are almost scale-independent ($Q_z=0.5\sim1.0$). In addition, their dependence on the morphology, color and luminosity is not statistically significant. Given the robust morphological, color and luminosity dependences of the two-point correlation function, this implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that simple deterministic linear relation with the underlying mass could not explain our measurements on these scales.

preprint2002arXiv

The Kurtosis of the Cosmic Shear Field

We study the fourth-order moment of the cosmic shear field using the dark matter halo approach to describe the nonlinear gravitational evolution of structure in the universe. Since the third-order moment of the shear field vanishes because of symmetry, non-Gaussian signatures in its one-point statistics emerge at the fourth-order level. We argue that the shear kurtosis parameter S_4 = <gamma^4>/<gamma^2>^3 may be more directly applicable to realistic data than the well-studied higher-order statistics of the convergence field, since obtaining the convergence requires a non-local reconstruction from the measured shear field. We compare our halo model predictions for the variance, skewness and kurtosis of lensing fields with ray-tracing simulations of cold dark matter models and find good agreement. The shear kurtosis calculation is made tractable by developing approximations for fast and accurate evaluations of the 8-dimensional integrals needed to obtain the kurtosis. We show that on small scales it is dominated by correlations within halos more massive than 10^14 solar masses. The shear kurtosis is sensitive to the mass density parameter of the universe, Omega, and has relatively weak dependences on other parameters. The approximations we develop for the third- and fourth-order moments allow for accurate halo model predictions for the 3-dimensional mass distribution as well. We demonstrate their accuracy in the small scale regime, below 2 Mpc, where analytical approaches used in the literature so far cease to be accurate.

preprint2000arXiv

Cosmic shear and halo abundances: Analytical versus numerical results

The aperture mass has been shown in a series of recent publications to be a useful quantitative tool for weak lensing studies, ranging from cosmic shear to the detection of a mass-selected sample of dark matter haloes. Quantitative analytical predictions for the aperture mass have been based on a number of simplifying assumptions. In this paper, we test the reliability of these assumptions and the quality of the analytic approximations, using ray-tracing simulations through a cosmological density field generated by very large N-body simulations. We find that those analytic predictions which take into account the non-linear evolution of the matter distribution, such as the dispersion of the aperture mass and the halo abundance, are surprisingly accurately reproduced with our numerical results, whereas the predictions for the skewness, based on quasi-linear theory, are rather imprecise. In particular, we verify numerically that the probability distribution of the aperture mass decreases exponentially for values much larger than the rms. Given the good overall agreement, comparisons between the observed distribution of the aperture mass and the theoretical values provide a powerful tool for testing cosmological models.

preprint1996arXiv

Does Gravitational Clustering Stabilize On Small Scales?

The stable clustering hypothesis is a key analytical anchor on the nonlinear dynamics of gravitational clustering in cosmology. It states that on sufficiently small scales the mean pair velocity approaches zero, or equivalently, that the mean number of neighbours of a particle remains constant in time at a given physical separation. In this paper we use N-body simulations of scale free spectra P(k) \propto k^n with -2 \leq n \leq 0 and of the CDM spectrum to test for stable clustering using the time evolution and shape of the correlation function ξ(x,t), and the mean pair velocity on small scales. For all spectra the results are consistent with the stable clustering predictions on the smallest scales probed, x < 0.07 x_{nl}(t), where x_{nl}(t) is the correlation length. The measured stable clustering regime corresponds to a typical range of 200 \lsim ξ\lsim 2000, though spectra with more small scale power approach the stable clustering asymptote at larger values of ξ. We test the amplitude of ξpredicted by the analytical model of Sheth \& Jain (1996), and find agreement to within 20\% in the stable clustering regime for nearly all spectra. For the CDM spectrum the nonlinear ξis accurately approximated by this model with n \simeq -2 on physical scales \lsim 100-300 h^{-1} kpc for σ_8 = 0.5-1, and on smaller scales at earlier times. The growth of ξfor CDM-like models is discussed in the context of a power law parameterization often used to describe galaxy clustering at high redshifts. The growth parameter εis computed as a function of time and length scale, and found to be larger than 1 in the moderately nonlinear regime -- thus the growth of ξis much faster on scales of interest than is commonly assumed.

preprint1996arXiv

The non-linear correlation function and the shapes of virialized halos

The correlation function xi(r) of matter in the non-linear regime is assumed to be determined by the density profiles rho(r) and the mass distribution n(M) of virialized halos. The Press--Schechter approach is used to compute n(M), and the stable clustering hypothesis is used to determine the density profiles of these Press--Schechter halos. Thus, the shape and amplitude of xi(r) on small scales is related to the initial power spectrum of density fluctuations. The case of clustering from scale-free initial conditions is treated in detail. If n is the slope of the initial power spectrum of density fluctuations, then stable clustering requires that xi(r)\propto r^{-gamma}, where gamma is a known function of n. If halo--halo correlations can be neglected, then rho(r)\propto r^{-epsilon}, where epsilon = (gamma+3)/2 = 3(4+n)/(5+n). For all values of n of current interest, this slope is steeper than the value 3(3+n)/(4+n) that was obtained by Hoffman & Shaham in their treatment of the shapes of the outer regions of collapsed halos. Our main result is a prediction for the amplitude of the non-linear correlation function. The predicted amplitude and its dependence on n are in good quantitative agreement with N-body simulations of self-similar clustering. If stable clustering is a good approximation only inside the half-mass radii of Press--Schechter halos, then the density contrast required for the onset of stable clustering can be estimated. This density contrast is in the range ~300-600 and increases with the initial slope n, in agreement with estimates from N-body simulations.