Source author record

Wayne Hu

Wayne Hu 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

93works
7topics
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

93 published item(s)

preprint2022arXiv

Accurate effective fluid approximation for ultralight axions

Ultralight axions are theoretically interesting and phenomenologically rich dark sector candidates, but they are difficult to track across cosmological timescales because of their fast oscillations. We resolve this problem by developing a novel method to evolve them efficiently and accurately. We first construct an exact effective fluid which at late times matches the axion but which evolves in a simple way. We then approximate this evolution with a carefully chosen equation of state and sound speed. With our scheme we find that we can obtain subpercent accuracy for the linear theory suppression of axion density fluctuations relative to that of cold dark matter without tracking even a single complete oscillation of the axion field. We use our technique to test other approximation schemes and to provide a fitting formula for the transfer function for the matter power spectrum in linear theory in axion models. Implementing our approach in existing cosmological axion codes is straightforward and will help unleash the potential of high-precision next-generation experiments.

preprint2022arXiv

Amplification of Primordial Perturbations from the Rise or Fall of the Inflaton

The next generation of cosmic microwave background, gravitational wave, and large scale structure, experiments will provide an unprecedented opportunity to probe the primordial power spectrum on small scales. An exciting possibility for what lurks on small scales is a sharp rise in the primordial power spectrum: This can lead to the formation of primordial black holes, providing a dark matter candidate or the black holes observed by the LIGO-Virgo collaboration. In this work we develop a mechanism for the amplification of the small-scale primordial power spectrum, in the context of single-field inflation with a step-like feature in the inflaton potential. Specifically, we consider both the upward and the downward step in the potential. We also discuss the possibility of the strong coupling between perturbations because the rapid changes of the potential derivatives with the time-dependent field value, caused by the step-like feature, could make the coupling stronger. As a result, we find that the perturbations can remain weakly coupled yet sufficiently enhanced if the step realizes the rapid changes of the potential derivatives in some fraction of an e-fold, $\mathcal O(\mathcal P_{\mathcal R}^{1/2}) \lesssim ΔN < 1$, where $\mathcal P_\mathcal R$ is the power spectrum of the curvature perturbation at that time. We also discuss the PBH formation rate from the inflaton trapping at the local minimum, which can occur in the potential with an upward step.

preprint2022arXiv

Modified gravitational wave propagation with higher modes and its degeneracies with lensing

Low-energy alternatives to General Relativity (GR) generically modify the phase of gravitational waves (GWs) during their propagation. As detector sensitivities increase, it becomes key to understand how these modifications affect the GW higher modes and to disentangle possible degeneracies with astrophysical phenomena. We apply a general formalism -- the WKB approach -- for solving analytically wave propagation in the spatial domain with a modified dispersion relation (MDR). We compare this WKB approach to applying a stationary phase approximation (SPA) in the temporal domain with time delays associated to the group or particle velocity. To this end, we extend the SPA to generic signals with higher modes, keeping careful track of reference phases and arrival times. We find that the WKB approach coincides with the SPA using the group velocity, in agreement with the principles of wave propagation. We then explore the degeneracies between a GW propagation with an MDR and a strongly-lensed GW in GR, since the latter can introduce a frequency-independent phase shift which is not degenerate with source parameters in the presence of higher modes. We find that for a particular MDR there is an exact degeneracy for wave propagation, unlike with the SPA for particle propagation. For the other cases, we search for the values of the MDR parameters that minimize the $χ^2$ and conclude that strongly-lensed GR GWs could be misinterpreted as GWs in modified gravity. Future MDR constraints with higher mode GWs should include the possibility of frequency-independent phase shifts, allowing for the identification of modified gravity and strong lensing distortions at the same time.

preprint2022arXiv

The Early Dark Sector, the Hubble Tension, and the Swampland

We consider the interplay of the Early Dark Energy (EDE) model, the Swampland Distance Conjecture (SDC), and cosmological parameter tensions. EDE is a proposed resolution of the Hubble tension relying upon a near-Planckian scalar field excursion, while the SDC predicts an exponential sensitivity of masses of other fields to such an excursion, $m\propto e^{-c|Δϕ|/M_{\rm pl}}$ with $c\sim{\cal O}(1)$. Meanwhile, EDE is in tension with large-scale structure (LSS) data, due to shifts in the standard $Λ$CDM parameters necessary to fit the cosmic microwave background (CMB). One might hope that a proper treatment of the model, e.g., accounting for the SDC, may ameliorate the tension with LSS. Motivated by these considerations, we introduce the Early Dark Sector (EDS) model, wherein the mass of dark matter is exponentially sensitive to super-Planckian field excursions of the EDE scalar. The EDS model exhibits new phenomenology in both the early and late universe, the latter due to an EDE-mediated dark matter self-interaction. This dark matter-philic "fifth force", while constrained to be small, remains active in the late universe and is not screened in virialized halos. We find that the new interaction with dark matter partially resolves the LSS tension. However, the marginalized posteriors are nonetheless consistent with $f_{\rm EDE}=0$ at 95$\%$ CL once the Dark Energy Survey Year 3 measurement of $S_8$ is included. We study constraints on the model from Atacama Cosmology Telescope data, and find a factor of two improvement on the error bar on the SDC parameter $c$, along with an increased preference for the EDE component. We discuss the implications of these constraints for the SDC, and find the tightest observational constraints to date on a swampland parameter, suggesting that an EDE description of cosmological data is in tension with the SDC.

preprint2021arXiv

Hubble constant difference between CMB lensing and BAO measurements

We apply a tension metric $Q_\textrm{UDM}$, the update difference in mean parameters, to understand the source of the difference in the measured Hubble constant $H_0$ inferred with cosmic microwave background lensing measurements from the Planck satellite ($H_0=67.9^{+1.1}_{-1.3}\, \mathrm{km/s/Mpc}$) and from the South Pole Telescope ($H_0=72.0^{+2.1}_{-2.5}\, \mathrm{km/s/Mpc}$) when both are combined with baryon acoustic oscillation (BAO) measurements with priors on the baryon density (BBN). $Q_\textrm{UDM}$ isolates the relevant parameter directions for tension or concordance where the two data sets are both informative, and aids in the identification of subsets of data that source the observed tension. With $Q_\textrm{UDM}$, we uncover that the difference in $H_0$ is driven by the tension between Planck lensing and BAO+BBN, at probability-to-exceed of 6.6%. Most of this mild tension comes from the galaxy BAO measurements parallel to the line of sight. The redshift dependence of the parallel BAOs pulls both the matter density $Ω_m$ and $H_0$ high in $Λ$CDM, but these parameter anomalies are usually hidden when the BAO measurements are combined with other cosmological data sets with much stronger $Ω_m$ constraints.

preprint2021arXiv

Primordial Black Holes Arise When The Inflaton Falls

Primordial Black Holes (PBHs) have entered the forefront of theoretical cosmology, due their potential role in phenomena ranging from gravitational waves, to dark matter, to galaxy formation. While producing PBHs from inflationary fluctuations naively would seem to require a large deceleration of the inflaton from its velocity at the horizon exit of CMB scales, in this work we demonstrate that an acceleration from a relatively small downward step in the potential that is transited in much less than an e-fold amplifies fluctuations as well. Depending on the location of the step, such PBHs could explain dark matter or the black holes detected by the gravitational wave interferometers. The perturbation enhancement has a natural interpretation as particle production due to the non-adiabatic transition associated with the step.

preprint2021arXiv

Spectral Distortion Anisotropy from Inflation for Primordial Black Holes

Single field inflationary models that seek to greatly enhance small scale power in order to form primordial black holes predict both a squeezed bispectrum that is enhanced by this small scale power and a potentially detectable enhancement of CMB spectral distortions. Despite this combination, spectral distortion anisotropy on CMB scales remains small since the squeezed bispectrum represents an unobservable modulation of the scale rather than local amplitude for the short wavelength acoustic power that dissipates and forms the $μ$ spectral distortion. The leading order amplitude effect comes from the local modulation of acoustic dissipation at the beginning of the $μ$ epoch at the end of thermalization by a long wavelength mode that is correlated with CMB anisotropy itself. Compensating factors from the suppression by the square of the ratio the comoving horizon at thermalization to the smallest detectable primary CMB scales ($\sim 0.0005$) and maximal allowed enhancement of $μ$ ($\sim 5000$) leaves a signal in the $μT$ cross spectrum that is still well beyond the capabilities of PIXIE or LiteBIRD due to sensitivity and resolution while remaining much larger than in single field slow roll inflation and potentially observable.

preprint2020arXiv

Apparent Superluminality of Lensed Gravitational Waves

We consider gravitational wave (GW) sources with an associated electromagnetic (EM) counterpart, and analyze the time delay between both signals in the presence of lensing. If GWs have wavelengths comparable to the Schwarzschild radius of astrophysical lenses, they must be treated with wave optics, whereas EM waves are typically well within the approximation of geometric optics. With concrete examples, we confirm that the GW signal never arrives before its EM counterpart, if both are emitted at the same time. However, during the inspiral of a binary, peaks of the GW waveform can arrive before their EM counterpart. We stress this is only an apparent superluminality since the GW waveform is both distorted and further delayed with respect to light. In any case, measuring the multi-messenger time delay and correctly interpreting it has important implications for unveiling the distribution of lenses, testing the nature of gravity, and probing the cosmological expansion history.

preprint2020arXiv

Can Late Dark Energy Transitions Raise the Hubble constant?

Late times dark energy transitions at redshifts $z \ll 0.1$ can raise the predicted value of the Hubble constant to the SH0ES value, $74.03\pm 1.42$ (km s$^{-1}$ Mpc$^{-1})$ or more, while providing an equally good fit as $Λ$CDM at $67.73 \pm 0.41$ to higher redshift data, in particular from the cosmic microwave background and baryon acoustic oscillations. These models however do not fully resolve the true source of tension between the distance ladder and high redshift observations: the local calibration of supernovae luminosities well out into the Hubble flow. When tested in this manner by transferring the SH0ES calibration to the Pantheon supernovae dataset, the ability of such transitions to raise the Hubble constant is reduced to $69.17 \pm 1.09$. Such an analysis should also be used when testing any dynamical dark energy model which can produce similarly fine features in redshift or local void models.

preprint2020arXiv

Effective field theory of degenerate higher-order inflation

We extend the effective field theory of inflation to a general Lagrangian constructed from Arnowitt-Deser-Misner variables that encompasses the most general interactions with up to second derivatives of the scalar field whose background breaks temporal diffeomorphism invariance. Degeneracy conditions, corresponding to 8 distinct types -- only one of which corresponds to known degenerate higher-order scalar-tensor models -- provide necessary conditions for eliminating the Ostrogradsky ghost in a covariant theory at the level of the quadratic action in unitary gauge. Novel implications of the degenerate higher-order system for the Cauchy problem are illustrated with the phase space portrait of an explicit inflationary example: not all field configurations lead to physical solutions for the metric even for positive potentials; solutions are unique for a given configuration only up to a branch choice; solutions on one branch can apparently end at nonsingular points of the metric and their continuation on alternate branches lead to nonsingular bouncing solutions; unitary gauge perturbations can go unstable even when degenerate terms in the Lagrangian are infinitesimal. The attractor solution leads to an inflationary scenario where slow-roll parameters vary and running of the tilt can be large even with no explicit features in the potential far from the end of inflation, requiring the optimized slow-roll approach for predicting observables.

preprint2020arXiv

Lensing-like tensions in the Planck legacy release

We analyze the final release of the Planck satellite data to constrain the gravitational lensing potential in a model-independent manner. The amount of lensing determined from the smoothing of the acoustic peaks in the temperature and polarization power spectra is 2$σ$ too high when compared with the measurements using the lensing reconstruction and 2.8$σ$ too high when compared with $Λ$CDM expectation based on the "unlensed" portion of the temperature and polarization power spectra. The largest change from the previous data release is the $Λ$CDM expectation, driven by improved constraints to the optical depth to reionization. The anomaly still is inconsistent with actual gravitational lensing, given that the lensing reconstruction constraints are discrepant independent of the model. Within the context of $Λ$CDM, improvements in its parameter constraints from lensing reconstruction bring this tension to 2.1$σ$ and from further adding baryon acoustic oscillation and Pantheon supernova data to a marginally higher 2.2$σ$. Once these other measurements are included, marginalizing this lensing-like anomaly cannot substantially resolve tensions with low-redshift measurements of $H_0$ and $S_8$ in $Λ$CDM, $Λ$CDM+$N_\mathrm{eff}$ or $Λ$CDM+$\sum m_ν$; furthermore the artificial strengthening of constraints on $\sum m_ν$ is less than 20%.

preprint2020arXiv

Primordial Black Holes as Dark Matter through Higgs Criticality

We study the dynamics of a spectator Higgs field which stochastically evolves during inflation onto near-critical trajectories on the edge of a runaway instability. We show that its fluctuations do not produce primordial black holes (PBHs) in sufficient abundance to be the dark matter, nor do they produce significant second-order gravitational waves. First we show that the Higgs produces larger fluctuations on CMB scales than on PBH scales, itself a no-go for a viable PBH scenario. Then we track the superhorizon perturbations nonlinearly through reheating using the delta N formalism to show that they are not converted to large curvature fluctuations. Our conclusions hold regardless of any fine-tuning of the Higgs field for both the Standard Model Higgs and for Higgs potentials modified to prevent unbounded runaway.

preprint2019arXiv

Quantifying concordance of correlated cosmological data sets

We develop estimators of agreement and disagreement between correlated cosmological data sets. These account for data correlations when computing the significance of both tensions and excess confirmation while remaining statistically optimal. We discuss and thoroughly characterize different approaches commenting on the ones that have the best behavior in practical applications. We complement the calculation of their statistical distribution within the Gaussian model with one estimator that takes non-Gaussianities fully into account. To illustrate the use of our techniques, we apply these estimators to supernovae measurements of the distance-redshift relation, absolutely calibrated by the local distance ladder. The suite of best estimators that we discuss finds results that are in excellent agreement between estimators and find no indications of significant internal inconsistencies in this data set above the $1\%$ probability threshold. This shows the robustness of local determinations of the Hubble constant to features in the distance-redshift relation.

preprint2016arXiv

CMB-S4 Science Book, First Edition

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

preprint2016arXiv

Double screening

Attempts to modify gravity in the infrared typically require a screening mechanism to ensure consistency with local tests of gravity. These screening mechanisms fit into three broad classes; we investigate theories which are capable of exhibiting more than one type of screening. Specifically, we focus on a simple model which exhibits both Vainshtein and kinetic screening. We point out that due to the two characteristic length scales in the problem, the type of screening that dominates depends on the mass of the sourcing object, allowing for different phenomenology at different scales. We consider embedding this double screening phenomenology in a broader cosmological scenario and show that the simplest examples that exhibit double screening are radiatively stable.

preprint2016arXiv

Lensing Bias to CMB Measurements of Compensated Isocurvature Perturbations

Compensated isocurvature perturbations (CIPs) are modes in which the baryon and dark matter density fluctuations cancel. They arise in the curvaton scenario as well as some models of baryogenesis. While they leave no observable effects on the cosmic microwave background (CMB) at linear order, they do spatially modulate two-point CMB statistics and can be reconstructed in a manner similar to gravitational lensing. Due to the similarity between the effects of CMB lensing and CIPs, lensing contributes nearly Gaussian random noise to the CIP estimator that approximately doubles the reconstruction noise power. Additionally, the cross correlation between lensing and the integrated Sachs-Wolfe (ISW) effect generates a correlation between the CIP estimator and the temperature field even in the absence of a correlated CIP signal. For cosmic-variance limited temperature measurements out to multipoles $l \leq 2500$, subtracting a fixed lensing bias degrades the detection threshold for CIPs by a factor of $1.3$, whether or not they are correlated with the adiabatic mode.

preprint2016arXiv

Quintessential Scale Dependence from Separate Universe Simulations

By absorbing fluctuations into a local background, separate universe simulations provide a powerful technique to characterize the response of small-scale observables to the long-wavelength density fluctuations, for example those of the power spectrum and halo mass function which lead to the squeezed-limit $n$-point function and halo bias, respectively. Using quintessence dark energy as the paradigmatic example, we extend these simulation techniques to cases where non-gravitational forces in other sectors establish a Jeans scale across which the growth of density fluctuations becomes scale dependent. By characterizing the separate universes with matching background expansion histories, we show that the power spectrum and mass function responses depend on whether the long-wavelength mode is above or below the Jeans scale. Correspondingly, the squeezed bispectrum and halo bias also become scale dependent. Models of bias that are effectively local in the density field at a single epoch, initial or observed, cannot describe this effect which highlights the importance of temporal nonlocality in structure formation. Validated by these quintessence tests, our techniques are applicable to a wide range of models where the complex dynamics of additional fields affect the clustering of matter in the linear regime and it would otherwise be difficult to simulate their impact in the nonlinear regime.

preprint2016arXiv

Self-accelerating Massive Gravity: Hidden Constraints and Characteristics

Self-accelerating backgrounds in massive gravity provide an arena to explore the Cauchy problem for derivatively coupled fields that obey complex constraints which reduce the phase space degrees of freedom. We present here an algorithm based on the Kronecker form of a matrix pencil that finds all hidden constraints, for example those associated with derivatives of the equations of motion, and characteristic curves for any 1+1 dimensional system of linear partial differential equations. With the Regge-Wheeler-Zerilli decomposition of metric perturbations into angular momentum and parity states, this technique applies to fully 3+1 dimensional perturbations of massive gravity around any isotropic self-accelerating background. Five spin modes of the massive graviton propagate once the constraints are imposed: two spin-2 modes with luminal characteristics present in the massless theory as well as two spin-1 modes and one spin-0 mode. Although the new modes all possess the same - typically spacelike - characteristic curves, the spin-1 modes are parabolic while the spin-0 modes are hyperbolic. The joint system, which remains coupled by non-derivative terms, cannot be solved as a simple Cauchy problem from a single non-characteristic surface. We also illustrate the generality of the algorithm with other cases where derivative constraints reduce the number of propagating degrees of freedom or order of the equations.

preprint2016arXiv

Separate Universe Consistency Relation and Calibration of Halo Bias

Linear halo bias is the response of dark matter halo number density to a long wavelength fluctuation in the dark matter density. Using abundance matching between separate universe simulations which absorb the latter into a change in the background, we test the consistency relation between the change in a one point function, the halo mass function, and a two point function, the halo-matter cross correlation in the long wavelength limit. We find excellent agreement between the two at the $1-2\%$ level for average halo biases between $1 \lesssim \bar b_1 \lesssim 4$ and no statistically significant deviations at the $4-5\%$ level out to $\bar b_1 \approx 8$. Halo bias inferred assuming instead a universal mass function is significantly different and inaccurate at the 10\% level or more. The separate universe technique provides a way of calibrating linear halo bias efficiently for even highly biased rare halos in the $Λ$CDM model. Observational violation of the consistency relation would indicate new physics, e.g.~in the dark matter, dark energy or primordial non-Gaussianity sectors.

preprint2016arXiv

Separating the Universe into the Real and Fake

The separate universe technique provides a means of establishing consistency relations between short wavelength observables and the long wavelength matter density fluctuations within which they evolve by absorbing the latter into the cosmological background. We extend it to cases where non-gravitational forces introduce a Jeans scale in other species like dynamical dark energy or massive neutrinos. The technique matches the synchronous gauge matter density fluctuations to the local expansion using the acceleration equation and accounts for the temporal nonlocality and scale dependence of the long wavelength response of small scale matter observables, e.g. the nonlinear power spectrum, halo abundance and the implied halo bias, and $N$-point correlation functions. Above the Jeans scale, the local Friedmann equation relates the expansion to real energy densities and a curvature that is constant in comoving coordinates. Below the Jeans scale, the curvature evolves and acts like a fake density component. In all cases, the matter evolution on small scales is correctly modeled as we illustrate using scalar field dark energy with adiabatic or isocurvature initial conditions across the Jeans scale set by its finite sound speed.

preprint2015arXiv

Compensated isocurvature perturbations in the curvaton model

Primordial fluctuations in the relative number densities of particles, or isocurvature perturbations, are generally well constrained by cosmic microwave background (CMB) data. A less probed mode is the compensated isocurvature perturbation (CIP), a fluctuation in the relative number densities of cold dark matter and baryons. In the curvaton model, a subdominant field during inflation later sets the primordial curvature fluctuation $ζ$. In some curvaton-decay scenarios, the baryon and cold dark matter isocurvature fluctuations nearly cancel, leaving a large CIP correlated with $ζ$. This correlation can be used to probe these CIPs more sensitively than the uncorrelated CIPs considered in past work, essentially by measuring the squeezed bispectrum of the CMB for triangles whose shortest side is limited by the sound horizon. Here, the sensitivity of existing and future CMB experiments to correlated CIPs is assessed, with an eye towards testing specific curvaton-decay scenarios. The planned CMB Stage 4 experiment could detect the largest CIPs attainable in curvaton scenarios with more than 3$σ$ significance. The significance could improve if small-scale CMB polarization foregrounds can be effectively subtracted. As a result, future CMB observations could discriminate between some curvaton-decay scenarios in which baryon number and dark matter are produced during different epochs relative to curvaton decay. Independent of the specific motivation for the origin of a correlated CIP perturbation, cross-correlation of CIP reconstructions with the primary CMB can improve the signal-to-noise ratio of a CIP detection. For fully correlated CIPs the improvement is a factor of $\sim$2$-$3.

preprint2015arXiv

Running from Features: Optimized Evaluation of Inflationary Power Spectra

In models like axion monodromy, temporal features during inflation which are not associated with its ending can produce scalar, and to a lesser extent, tensor power spectra where deviations from scale-free power law spectra can be as large as the deviations from scale invariance itself. Here the standard slow-roll approach breaks down since its parameters evolve on an efolding scale $ΔN$ much smaller than the efolds to the end of inflation. Using the generalized slow-roll approach, we show that the expansion of observables in a hierarchy of potential or Hubble evolution parameters comes from a Taylor expansion of the features around an evaluation point that can be optimized. Optimization of the leading-order expression provides a sufficiently accurate approximation for current data as long as the power spectrum can be described over the well-observed few efolds by the local tilt and running. Standard second-order approaches, often used in the literature, ironically are worse than leading-order approaches due to inconsistent evaluation of observables. We develop a new optimized next-order approach which predicts observables to $10^{-3}$ even for $ΔN\sim 1$ where all parameters in the infinite hierarchy are of comparable magnitude. For models with $ΔN \ll 1$, the generalized slow-roll approach provides integral expressions that are accurate to second order in the deviation from scale invariance. Their evaluation in the monodromy model provides highly accurate explicit relations between the running oscillation amplitude, frequency and phase in the curvature spectrum and parameters of the potential.

preprint2015arXiv

Self-accelerating Massive Gravity: Superluminality, Cauchy Surfaces and Strong Coupling

Self-accelerating solutions in massive gravity provide explicit, calculable examples that exhibit the general interplay between superluminality, the well-posedness of the Cauchy problem, and strong coupling. For three particular classes of vacuum solutions, one of which is new to this work, we construct the conformal diagram for the characteristic surfaces on which isotropic stress-energy perturbations propagate. With one exception, all solutions necessarily possess spacelike characteristics, indicating perturbative superluminality. Foliating the spacetime with these surfaces gives a pathological frame where kinetic terms of the perturbations vanish, confusing the Hamiltonian counting of degrees of freedom. This frame dependence distinguishes the vanishing of kinetic terms from strong coupling of perturbations or an ill-posed Cauchy problem. We give examples where spacelike characteristics do and do not originate from a point where perturbation theory breaks down and where spacelike surfaces do or do not intersect all characteristics in the past light cone of a given observer. The global structure of spacetime also reveals issues that are unique to theories with two metrics: in all three classes of solutions, the Minkowski fiducial space fails to cover the entire de Sitter spacetime allowing worldlines of observers to end in finite proper time at determinant singularities. Characteristics run tangent to these surfaces requiring {\it ad hoc} rules to establish continuity across singularities.

preprint2015arXiv

Stability of cosmological solutions in extended quasidilaton massive gravity

We consider the stability of self-accelerating solutions to extended quasidilaton massive gravity in the presence of matter. By making a second or extended fiducial metric dynamical in this model, matter can cause it to evolve from a Lorentzian to Euclidean signature, triggering a ghost instability. We study this possibility with scalar field matter as it can model a wide range of cosmological expansion histories. For the $Λ$CDM expansion history, stability considerations substantially limit the available parameter space while for a kinetic energy dominated expansion, no choice of quasidilaton parameters is stable. More generally these results show that there is no mechanism intrinsic to the theory to forbid such pathologies from developing from stable initial conditions and that stability can only be guaranteed for particular choices for the matter configuration.

preprint2014arXiv

$νΛ$CDM: Neutrinos help reconcile Planck with the Local Universe

Current measurements of the low and high redshift Universe are in tension if we restrict ourselves to the standard six parameter model of flat $Λ$CDM. This tension has two parts. First, the Planck satellite data suggest a higher normalization of matter perturbations than local measurements of galaxy clusters. Second, the expansion rate of the Universe today, $H_0$, derived from local distance-redshift measurements is significantly higher than that inferred using the acoustic scale in galaxy surveys and the Planck data as a standard ruler. The addition of a sterile neutrino species changes the acoustic scale and brings the two into agreement; meanwhile, adding mass to the active neutrinos or to a sterile neutrino can suppress the growth of structure, bringing the cluster data into better concordance as well. For our fiducial dataset combination, with statistical errors for clusters, a model with a massive sterile neutrino shows 3.5$σ$ evidence for a non-zero mass and an even stronger rejection of the minimal model. A model with massive active neutrinos and a massless sterile neutrino is similarly preferred. An eV-scale sterile neutrino mass -- of interest for short baseline and reactor anomalies -- is well within the allowed range. We caution that 1) unknown astrophysical systematic errors in any of the data sets could weaken this conclusion, but they would need to be several times the known errors to eliminate the tensions entirely; 2) the results we find are at some variance with analyses that do not include cluster measurements; and 3) some tension remains among the datasets even when new neutrino physics is included.

preprint2014arXiv

Bounds on non-adiabatic evolution in single-field inflation

We examine the regime of validity of $N$-point spectra predictions of single field inflation models that invoke transient periods of non-adiabatic evolution. Such models generate oscillatory features in these spectra spanning frequencies up to the inverse time scale of the transient feature. To avoid strong coupling of fluctuations in these theories this scale must be at least $\sim 10^{-2}/c_s$ of the Hubble time during inflation, where $c_s$ is the inflaton sound speed. We show that, in such models, the signal-to-noise ratio of the bispectrum is bounded from above by that of the power spectrum, implying that searches for features due to non-adiabatic evolution are best focussed first on the latter.

preprint2014arXiv

Generalized Slow Roll for Tensors

The recent BICEP2 detection of degree scale CMB B-mode polarization, coupled with a deficit of observed power in large angle temperature anisotropy, suggest that the slow-roll parameter $ε_H$, the fractional variation in the Hubble rate per efold, is both relatively large and may evolve from an even larger value on scales greater than the horizon at recombination. The relatively large tensor contribution implied also requires finite matching features in the tensor power spectrum for any scalar power spectrum feature proposed to explain anomalies in the temperature data. We extend the generalized slow-roll approach for computing power spectra, appropriate for such models where the slow-roll parameters vary, to tensor features where scalar features are large. This approach also generalizes the tensor-scalar consistency relation to be between the ratio of tensor and scalar sources and features in the two power spectra. Features in the tensor spectrum are generically suppressed by $ε_H$ relative those in the scalar spectrum and by the smoothness of the Hubble rate, which must obey covariant conservation of energy, versus its derivatives. Their detection in near future CMB data would indicate a fast roll period of inflation where $ε_H$ approaches order unity, allowed but not required by inflationary explanations of temperature anomalies.

preprint2014arXiv

Neutrinos help reconcile Planck measurements with both Early and Local Universe

In light of the recent BICEP2 B-mode polarization detection, which implies a large inflationary tensor-to-scalar ratio r_{0.05}=0.2^{+0.07}_{-0.05}, we re-examine the evidence for an extra sterile massive neutrino, originally invoked to account for the tension between the cosmic microwave background (CMB) temperature power spectrum and local measurements of the expansion rate H0 and cosmological structure. With only the standard active neutrinos and power-law scalar spectra, this detection is in tension with the upper limit of r<0.11 (95% confidence) from the lack of a corresponding low multipole excess in the temperature anisotropy from gravitational waves. An extra sterile species with the same energy density as is needed to reconcile the CMB data with H0 measurements can also alleviate this new tension. By combining data from the Planck and ACT/SPT temperature spectra, WMAP9 polarization, H_0, baryon acoustic oscillation and local cluster abundance measurements with BICEP2 data, we find the joint evidence for a sterile massive neutrino increases to DeltaNeff=0.98\pm 0.26 for the effective number and ms= 0.52\pm 0.13 eV for the effective mass or 3.8 sigma and 4 sigma evidence respectively. We caution the reader that these results correspond to a joint statistical evidence and, in addition, astrophysical systematic errors in the clusters and H0 measurements, and small-scale CMB data could weaken our conclusions.

preprint2014arXiv

Power Spectrum Super-Sample Covariance

We provide a simple, unified approach to describing the impact of super-sample covariance, or beat coupling, on power spectrum estimation in a finite-volume survey. For a wide range of survey volumes, the sample variance that arises from modes that are larger than the survey dominates the covariance of power spectrum estimators for modes much smaller than the survey. The deeply nonlinear version of this effect is known as halo sample variance. We show that all variants are unified by the matter trispectrum of squeezed configurations and that such configurations obey a consistency relation which relates them to the response of the power spectrum to a change in the background density. Our method also applies to statistics that are based on radial projections of the density field such as weak lensing shear. While we use the halo model for an analytic description to expose the nature of the effect, the consistency description enables an accurate calibration of the full effect directly from simulations. It also suggests that super-sample covariance may be viewed as an additional interesting signal rather than excess noise.

preprint2014arXiv

Self-accelerating Massive Gravity: Covariant Perturbation Theory

We undertake a complete and covariant treatment for the quadratic Lagrangian of all of the degrees of freedom of massive gravity with a fixed flat fiducial metric for arbitrary massive gravity parameters around any isotropic self-accelerating background solution. Generically, 3 out of 4 Stuckelberg degrees of freedom propagate in addition to the usual 2 tensor degrees of freedom of general relativity. The complete kinetic structure typically is only revealed at an order in the graviton mass that is equivalently to retaining curvature terms in a locally flat expansion. These results resolve several apparent discrepancies in the literature where zero degrees of freedom propagate in either special cases or approximate treatments as well as decoupling limit analyses which attempt to count longitudinal degrees of freedom.

preprint2014arXiv

Super-Sample CMB Lensing

Lensing of the CMB by modes that are larger than the size of the survey dilates intrinsic scales in the temperature and polarization fields and coherently shifts their observed power spectra with respect to the ensemble or all-sky mean. The effect can be simply encapsulated as a contribution to the power spectrum covariance matrix in accordance with the lensing trispectrum or as an additional parameter, the mean convergence in the field, for parameter estimation. It should be included for upcoming surveys that precisely measure acoustic polarization features deep into the damping tail at multipoles of $\ell \gtrsim 1500$ with less than $10\%$ of sky. Its omission may lead to seemingly conflicting values for the angular scale of the sound horizon which may then provide erroneous cosmological parameters when compared to baryon acoustic oscillation measurements.

preprint2014arXiv

Super-Sample Covariance in Simulations

Using separate universe simulations, we accurately quantify super-sample covariance (SSC), the typically dominant sampling error for matter power spectrum estimators in a finite volume, which arises from the presence of super survey modes. By quantifying the power spectrum response to a background mode, this approach automatically captures the separate effects of beat coupling in the quasilinear regime, halo sample variance in the nonlinear regime and a new dilation effect which changes scales in the power spectrum coherently across the survey volume, including the baryon acoustic oscillation scale. It models these effects at typically the few percent level or better with a handful of small volume simulations for any survey geometry compared with directly using many thousands of survey volumes in a suite of large volume simulations. The stochasticity of the response is sufficiently small that in the quasilinear regime, SSC can be alternately included by fitting the mean density in the volume with these fixed templates in parameter estimation. We also test the halo model prescription and find agreement typically at better than the 10% level for the response.

preprint2014arXiv

Super-Sample Signal

When extracting cosmological information from power spectrum measurements, we must consider the impact of super-sample density fluctuations whose wavelengths are larger than the survey scale. These modes contribute to the mean density fluctuation $δ_b$ in the survey and change the power spectrum in the same way as a change in the cosmological background. They can be simply included in cosmological parameter estimation and forecasts by treating $δ_b$ as an additional cosmological parameter enabling efficient exploration of its impact. To test this approach, we consider here an idealized measurement of the matter power spectrum itself in the $Λ$CDM cosmology though our techniques can readily be extended to more observationally relevant statistics or other parameter spaces. Using sub-volumes of large-volume $N$-body simulations for power spectra measured with respect to either the global or local mean density, we verify that the minimum variance estimator of $δ_b$ is both unbiased and has the predicted variance. Parameter degeneracies arise since the response of the matter power spectrum to $δ_b$ and cosmological parameters share similar properties in changing the growth of structure and dilating the scale of features especially in the local case. For matter power spectrum measurements, these degeneracies can lead in certain cases to substantial error degradation and motivates future studies of specific cosmological observables such as galaxy clustering and weak lensing statistics with these techniques.

preprint2013arXiv

Neutrino physics from future weak lensing surveys

Given recent indications of additional neutrino species and cosmologically significant neutrino masses, we analyze their signatures in the weak lensing shear power spectrum. We find that a shear deficit in the 20-40% range or excess in the 20-60% range cannot be explained by variations in parameters of the flat LambdaCDM model that are allowed by current observations of the expansion history from Type Ia supernovae, baryon acoustic oscillations, and local measures of the Hubble constant H_0, coupled with observations of the cosmic microwave background from WMAP9 and the SPT 2500 square degree survey. Hence such a shear deficit or excess would indicate large masses or extra species, respectively, and we find this to be independent of the flatness assumption. We also discuss the robustness of these predictions to cosmic acceleration physics and the means by which shear degeneracies in joint variation of mass and species can be broken.

preprint2013arXiv

Non-Gaussian structure of the lensed CMB power spectra covariance matrix

Gravitational lensing of the Cosmic Microwave Background (CMB) encodes cosmological information in the observed anisotropies of temperature and polarization. Accurate extraction of this additional information requires precise modeling of the covariance matrix of the power spectra of observed CMB fields. We introduce a new analytical model to describe the non-Gaussian structure of this covariance matrix and display the importance of second-order terms that were previously neglected. When compared with direct numerical simulations our model captures parameter errors to better than a few percent for cases where the non-Gaussianity causes an order unity degradation in errors. We also provide a detailed comparison between the information content of lensed CMB power spectra and ideal reconstruction of the lensing potential. We illustrate the impact of the non-Gaussian terms in the power spectrum covariance by providing Fisher errors on the sum of the masses of the neutrinos, the dark energy equation of state, and the curvature of the Universe.

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.

preprint2013arXiv

Self-accelerating Massive Gravity: Bimetric Determinant Singularities

The existence of two metrics in massive gravity theories in principle allows solutions where there are singularities in new scalar invariants jointly constructed from them. These configurations occur when the two metrics differ substantially from each other, as in black hole and cosmological solutions. The simplest class of such singularities are determinant singularities. We investigate whether the dynamics of bimetric massive gravity -- where the second metric is allowed to evolve jointly with the spacetime metric -- can avoid these singularities. We investigate whether the dynamics of bimetric massive gravity, where the second metric is allowed to evolve jointly with the spacetime metric, averts the simplest class of such singularities, namely determinant singularities. We show that it is still possible to specify non-singular initial conditions that evolve to a determinant singularity. Determinant singularities are a feature of massive gravity of both fixed and dynamical metric type.

preprint2013arXiv

Self-accelerating Massive Gravity: How Zweibeins Walk through Determinant Singularities

The theory of massive gravity possesses ambiguities when the spacetime metric evolves far from the non-dynamical fiducial metric used to define it. We explicitly construct a spherically symmetric example case where the metric evolves to a coordinate-independent determinant singularity which does not exist in the initial conditions. Both the metric and the vierbein formulation of the theory are ill-defined at this point. In unitary gauge, the chart of the spacetime ends at this point and does not cover the full spacetime whereas the spherically symmetric vierbeins, or zweibeins, of the fiducial metric become non-invertible and do not describe a valid metric. Nonetheless it is possible to continuously join a zweibein solution on the other side of the singularity which picks one of the degenerate solutions of the metric square root. This continuous solution is not the choice conventionally made in the previous literature. We also show that the Stueckelberg equations of motion on the self-accelerating branch prevent solutions from evolving to a more pathological situation in which the spacetime vierbeins lack a crucial symmetry with the fiducial vierbeins and real square roots fail to exist.

preprint2013arXiv

Self-accelerating Massive Gravity: Time for Field Fluctuations

The ghost-free theory of massive gravity has exact solutions where the effective stress energy generated by the graviton mass term is a cosmological constant for any isotropic metric. Since they are exact, these solutions mimic a cosmological constant in the presence of any matter-induced isotropic metric perturbation. In the Stueckelberg formulation, this stress energy is carried entirely by the spatial Stueckelberg field. We show that any stress energy carried by fluctuations in the spatial field away from the exact solution always decays away in an expanding universe. However, the dynamics of the spatial Stueckelberg field perturbation depend on the background temporal Stueckelberg field, which is equivalent to the unitary gauge time coordinate. This dependence resolves an apparent conflict in the existing literature by showing that there is a special unitary time choice for which the field dynamics and energy density perturbations vanish identically. In general, the isotropic system has a single dynamical degree of freedom requiring two sets of initial data; however, only one of these initial data choices will affect the observable metric. Finally, we construct cosmological solutions with a well-defined perturbative initial value formulation and comment on alternate solutions that evolve to singularities.

preprint2013arXiv

The abundance of voids and the excursion set formalism

We present measurements of the number density of voids in the dark matter distribution from a series of N-body simulations of a ΛCDM cosmology. We define voids as spherical regions of ρ_v = 0.2ρ_m around density minima in order to relate our results to the predicted abundances using the excursion set formalism. Using a linear underdensity of δ_v = -2.7, from a spherical evolution model, we find that a volume conserving model, which does not conserve number density in the mapping from the linear to nonlinear regime, matches the measured abundance to within 16% for a range of void radii 1< r(Mpc/h)<15. This model fixes the volume fraction of the universe which is in voids and assumes that voids of a similar size merge as they expand by a factor of 1.7 to achieve a nonlinear density of ρ_v = 0.2ρ_m today. We find that the model of Sheth & van de Weygaert (2004) for the number density of voids greatly overpredicts the abundances over the same range of scales. We find that the volume conserving model works well at matching the number density of voids measured from the simulations at higher redshifts, z=0.5 and 1, as well as correctly predicting the abundances to within 25% in a simulation of a matter dominated Ω_m = 1 universe. We examine the abundance of voids in the halo distribution and find fewer small, r<10 Mpc/h, voids and many more large, r>10 Mpc/h, voids compared to the dark matter. These results indicate that voids identified in the halo or galaxy distribution are related to the underlying void distribution in the dark matter in a complicated way which merits further study if voids are to be used as a precision probe of cosmology.

preprint2012arXiv

Cosmic String constraints from WMAP and the South Pole Telescope

The predictions of the inflationary LCDM paradigm match today's high-precision measurements of the cosmic microwave background anisotropy extremely well. The same data put tight limits on other sources of anisotropy. Cosmic strings are a particularly interesting alternate source to constrain. Strings are topological defects, remnants of inflationary-era physics that persist after the big bang. They are formed in a variety of models of inflation, including string theory models such as brane inflation. We assume a "Nambu-Goto" model for strings, approximated by a collection of unconnected segments with zero width, and show that measurements of temperature anisotropy by the South Pole Telescope break a parameter degeneracy in the WMAP data, permitting us to place a strong upper limit on the possible string contribution to the CMB anisotropy: the power sourced by zero-width strings must be <1.75% (95% CL) of the total or the string tension Gmu <1.7x10^{-7}. These limits imply that the best hope for detecting strings in the CMB will come from B-mode polarization measurements at arcminute scales rather than the degree scale measurements pursued for gravitational wave detection.

preprint2012arXiv

Equivalence Principle Violation in Weakly Vainshtein-Screened Systems

Massive gravity, galileon and braneworld models that modify gravity to explain cosmic acceleration utilize the nonlinear field interactions of the Vainshtein mechanism to screen fifth forces in high density regimes. These source-dependent interactions cause apparent equivalence principle violations. In the weakly-screened regime violations can be especially prominent since the fifth forces are at near full strength. Since they can also be calculated perturbatively, we derive analytic solutions for illustrative cases: the motion of massive objects in compensated shells and voids and infall toward halos that are spherically symmetric. Using numerical techniques we show that these solutions are valid until the characteristic scale becomes comparable to the Vainshtein radius. We find a relative acceleration of more massive objects toward the center of a void and a reduction of the infall acceleration that increases with the mass ratio of the halos which can in principle be used to test the Vainshtein screening mechanism.

preprint2012arXiv

Fast Computation of Bispectrum Features with Generalized Slow Roll

We develop a fast technique based on the generalized slow roll (GSR) approach for computing the curvature bispectrum of inflationary models with features. We show that all triangle configurations can be expressed in terms of three simple integrals over the inflationary background with typical accuracy of better than ~20%. With a first order GSR approach the typical accuracy can be improved to better than the 5% level. We illustrate this technique with the step potential model that has been invoked to explain the WMAP temperature power spectrum glitches at ell ~ 20-40 and show that the maximum likelihood model falls short of observability by more than a factor of 100 in amplitude. We also explicitly demonstrate that the bispectrum consistency relation with the local slope of the power spectrum is satisfied for these models. In the GSR approach, the bispectrum arises from integrals of nearly the same function of the background slow-roll parameters as the power spectrum but with a stronger weight to the epoch before horizon crossing. Hence this technique enables reverse engineering of models with large bispectrum but small power spectrum features.

preprint2012arXiv

Fast Computation of First-Order Feature-Bispectrum Corrections

Features in the inflaton potential that are traversed in much less than an e-fold of the expansion can produce observably large non-Gaussianity. In these models first order corrections to the curvature mode function evolution induce effects at second order in the slow roll parameters that are generically greater than ~ 10% and can reach order unity for order unity power spectrum features. From a complete first order expression in generalized slow-roll, we devise a computationally efficient method that is as simple to evaluate as the leading order one and implements consistency relations in a controlled fashion. This expression matches direct numerical computation for step potential models of the dominant bispectrum configurations to better than 1% when features are small and 10% when features are order unity.

preprint2012arXiv

Non-Gaussianity from Step Features in the Inflationary Potential

We provide analytic solutions for the power spectrum and bispectrum of curvature fluctuations produced by a step feature in the inflaton potential, valid in the limit that the step is short and sharp. In this limit, the bispectrum is strongly scale dependent and its effective non-linearity attains a large oscillatory amplitude. The perturbations to the curvature power spectrum, on the other hand, remain a small component on top of the usual spectrum of fluctuations generated by slow roll. We utilize our analytic solutions to assess the observability of the predicted non-Gaussian signatures and show that, if present, only very sharp steps on scales larger than ~ 2 Gpc are likely to be able to be detected by Planck. Such features are not only consistent with WMAP7 data, but can also improve its likelihood by 2 Delta ln L ~ 12 for two extra parameters, the step location and height. If this improvement were due to a slow roll violating step as considered here, a bispectrum or corresponding polarization power spectrum detection would provide definitive checks as to its primordial origin.

preprint2012arXiv

Quantum Stability of Chameleon Field Theories

Chameleon scalar fields are dark energy candidates which suppress fifth forces in high density regions of the universe by becoming massive. We consider chameleon models as effective field theories and estimate quantum corrections to their potentials. Requiring that quantum corrections be small, so as to allow reliable predictions of fifth forces, leads to an upper bound $m < 0.0073 (ρ/ 10 {\rm g\, cm}^{-3})^{1/3}$eV for gravitational strength coupling whereas fifth force experiments place a lower bound of $m>0.0042$\,eV. An improvement of less than a factor of two in the range of fifth force experiments could test all classical chameleon field theories whose quantum corrections are well-controlled and couple to matter with nearly gravitational strength regardless of the specific form of the chameleon potential.

preprint2012arXiv

Self-accelerating Massive Gravity: Exact solutions for any isotropic matter distribution

We present an exact solution to the equations of massive gravity that displays cosmological constant-like behavior for any spherically symmetric distribution of matter, including arbitrary time dependence. On this solution, the new degrees of freedom from the massive graviton generate a cosmological constant-like contribution to stress-energy that does not interact directly with other matter sources. When the effective cosmological constant contribution dominates over other sources of stress energy the cosmological expansion self-accelerates, even when no other dark-energy-like ingredients are present. The new degrees of freedom introduced by giving the graviton the mass do not respond to arbitrarily large radial or homogeneous perturbations from other matter fields on this solution. We comment on possible implications of this result.

preprint2012arXiv

Testing dark energy paradigms with weak gravitational lensing

Any theory invoked to explain cosmic acceleration predicts consistency relations between the expansion history, structure growth, and all related observables. Currently there exist high-quality measurements of the expansion history from Type Ia supernovae, the cosmic microwave background temperature and polarization spectra, and baryon acoustic oscillations. We can use constraints from these datasets to predict what future probes of structure growth should observe. We apply this method to predict what range of cosmic shear power spectra would be expected if we lived in a LambdaCDM universe, with or without spatial curvature, and what results would be inconsistent and therefore falsify the model. Though predictions are relaxed if one allows for an arbitrary quintessence equation of state $-1\le w(z)\le 1$, we find that any observation that rules out LambdaCDM due to excess lensing will also rule out all quintessence models, with or without early dark energy. We further explore how uncertainties in the nonlinear matter power spectrum, e.g. from approximate fitting formulas such as Halofit, warm dark matter, or baryons, impact these limits.

preprint2011arXiv

Chameleon Halo Modeling in f(R) Gravity

We model the chameleon effect on cosmological statistics for the modified gravity f(R) model of cosmic acceleration. The chameleon effect, required to make the model compatible with local tests of gravity, reduces force enhancement as a function of the depth of the gravitational potential wells of collapsed structure and so is readily incorporated into a halo model by including parameters for the chameleon mass threshold and rapidity of transition. We show that the abundance of halos around the chameleon mass threshold are enhanced by both the merging from below and the lack of merging to larger masses. This property also controls the power spectrum in the nonlinear regime and we provide a description of the transition to the linear regime that is valid for a wide range of f(R) models.

preprint2011arXiv

Complete WMAP Constraints on Bandlimited Inflationary Features

Using a principal component (PC) basis that accommodates order unity features in the slow roll parameters as fine as 1/10 of a decade across more than 2 decades of the inflationary expansion, we test slow roll and single field inflation with the WMAP7 data. Detection of any non-zero component would represent a violation of ordinary slow roll and indicate a feature in the inflaton potential or sound speed. Although one component shows a deviation at the 98% CL, it cannot be considered statistically significant given the 20 components tested. The maximum likelihood PC parameters only improves 2DeltalnL by 17 for the 20 parameters associated with known glitches in the WMAP power spectrum at multipoles l<60. We make model-independent predictions for the matching glitches in the polarization spectrum that would test their inflationary origin. This complete analysis for bandlimited features in the source function of generalized slow roll can be used to constrain parameters of specific models of the inflaton potential without requiring a separate likelihood analysis for each choice. We illustrate its use by placing bounds on the height and width of a step-like feature in the potential proposed to explain the glitch at 20<l<40. Even allowing for the presence of features in the temperature spectrum, single field inflation makes sharp falsifiable predictions for the acoustic peaks in the polarization whose violation would require extra degrees of freedom.

preprint2011arXiv

Constraints on f(R) gravity from probing the large-scale structure

We study cosmological constraints on metric f(R) gravity models that are designed to reproduce the LCDM expansion history with modifications to gravity described by a supplementary cosmological freedom, the Compton wavelength parameter B_0. We conduct a Markov chain Monte Carlo analysis on the parameter space, utilizing the geometrical constraints from supernovae distances, the baryon acoustic oscillation distances, and the Hubble constant, along with all of the cosmic microwave background data, including the largest scales, its correlation with galaxies, and a probe of the relation between weak gravitational lensing and galaxy flows. The strongest constraints, however, are obtained through the inclusion of data from cluster abundance. Using all of the data, we infer a bound of B_0<0.0011 at the 95% C.L.

preprint2011arXiv

Generalized Slow Roll for Non-Canonical Kinetic Terms

We show that the generalized slow-roll approach for calculating the power spectrum where the inflationary slow roll parameters are neither small nor slowly varying can be readily extended to models with non-canonical kinetic terms in the inflaton action. For example, rapid sound speed variations can arise in DBI models with features in the warp factor leading to features in the power spectrum. Nonetheless there remains a single source function for deviations that is simply related to the power spectrum. Empirical constraints on this source function can be readily interpreted in the context of features in the inflaton potential or sound speed.

preprint2011arXiv

Simultaneous Falsification of LCDM and Quintessence with Massive, Distant Clusters

Observation of even a single massive cluster, especially at high redshift, can falsify the standard cosmological framework consisting of a cosmological constant and cold dark matter (LCDM) with Gaussian initial conditions by exposing an inconsistency between the well-measured expansion history and the growth of structure it predicts. Through a likelihood analysis of current cosmological data that constrain the expansion history, we show that the LCDM upper limits on the expected number of massive, distant clusters are nearly identical to limits predicted by all quintessence models where dark energy is a minimally coupled scalar field with a canonical kinetic term. We provide convenient fitting formulas for the confidence level at which the observation of a cluster of mass M at redshift z can falsify LCDM and quintessence given cosmological parameter uncertainties and sample variance, as well as for the expected number of such clusters in the light cone and the Eddington bias factor that must be applied to observed masses. By our conservative confidence criteria, which equivalently require masses 3 times larger than typically expected in surveys of a few hundred square degrees, none of the presently known clusters falsify these models. Various systematic errors, including uncertainties in the form of the mass function and differences between supernova light curve fitters, typically shift the exclusion curves by less than 10% in mass, making current statistical and systematic uncertainties in cluster mass determination the most critical factor in assessing falsification of LCDM and quintessence.

preprint2011arXiv

The Bispectrum of f(R) Cosmologies

In this paper we analyze a suite of cosmological simulations of modified gravitational action f(R) models, where cosmic acceleration is induced by a scalar field that acts as a fifth force on all forms of matter. In particular, we focus on the bispectrum of the dark matter density field on mildly non-linear scales. For models with the same initial power spectrum, the dark matter bispectrum shows significant differences for cases where the final dark matter power spectrum also differs. Given the different dependence on bias of the galaxy power spectrum and bispectrum, bispectrum measurements can close the loophole of galaxy bias hiding differences in the power spectrum. Alternatively, changes in the initial power spectrum can also hide differences. By constructing LCDM models with very similar final non-linear power spectra, we show that the differences in the bispectrum are reduced (<4%) and are comparable with differences in the imperfectly matched power spectra. These results indicate that the bispectrum depends mainly on the power spectrum and less sensitively on the gravitational signatures of the f(R) model. This weak dependence of the matter bispectrum on gravity makes it useful for breaking degeneracies associated with galaxy bias, even for models beyond general relativity.

preprint2010arXiv

Analytic Description of DGP Perturbations on All Scales

We develop analytic solutions for the linear evolution of metric perturbations in the DGP braneworld modified gravity scenario including near-horizon and superhorizon modes where solutions in the bulk are required. These solutions apply to both the self-accelerating and normal branch and elucidate the nature of coordinate singularities and initial data in the bulk as well as their effect on perturbation evolution on the brane. Even on superhorizon scales, the evolution of metric perturbations is no longer necessarily scale free due to multiple resonances in the bulk. Based on these analytic solutions, we devise convenient fitting functions for the evolution that bridge the various spatial and temporal regimes. Compared with a direct numerical integration of the bulk equations, the fits are accurate at the percent level and are sufficient for current and upcoming observational tests.

preprint2010arXiv

Cluster Abundance in f(R) Gravity Models

As one of the most powerful probes of cosmological structure formation, the abundance of massive galaxy clusters is a sensitive probe of modifications to gravity on cosmological scales. In this paper, we present results from N-body simulations of a general class of f(R) models, which self-consistently solve the non-linear field equation for the enhanced forces. Within this class we vary the amplitude of the field, which controls the range of the enhanced gravitational forces, both at the present epoch and as a function of redshift. Most models in the literature can be mapped onto the parameter space of this class. Focusing on the abundance of massive dark matter halos, we compare the simulation results to a simple spherical collapse model. Current constraints lie in the large-field regime, where the chameleon mechanism is not important. In this regime, the spherical collapse model works equally well for a wide range of models and can serve as a model-independent tool for placing constraints on f(R) gravity from cluster abundance. Using these results, we show how constraints from the observed local abundance of X-ray clusters on a specific f(R) model can be mapped onto other members of this general class of models.

preprint2010arXiv

CMB Constraints on Principal Components of the Inflaton Potential

We place functional constraints on the shape of the inflaton potential from the cosmic microwave background through a variant of the generalized slow roll approximation that allows large amplitude, rapidly changing deviations from scale-free conditions. Employing a principal component decomposition of the source function G'~3(V'/V)^2 - 2V''/V and keeping only those measured to better than 10% results in 5 nearly independent Gaussian constraints that maybe used to test any single-field inflationary model where such deviations are expected. The first component implies < 3% variations at the 100 Mpc scale. One component shows a 95% CL preference for deviations around the 300 Mpc scale at the ~10% level but the global significance is reduced considering the 5 components examined. This deviation also requires a change in the cold dark matter density which in a flat LCDM model is disfavored by current supernova and Hubble constant data and can be tested with future polarization or high multipole temperature data. Its impact resembles a local running of the tilt from multipoles 30-800 but is only marginally consistent with a constant running beyond this range. For this analysis, we have implemented a ~40x faster WMAP7 likelihood method which we have made publicly available.

preprint2010arXiv

Figures of merit for present and future dark energy probes

We compare current and forecasted constraints on dynamical dark energy models from Type Ia supernovae and the cosmic microwave background using figures of merit based on the volume of the allowed dark energy parameter space. For a two-parameter dark energy equation of state that varies linearly with the scale factor, and assuming a flat universe, the area of the error ellipse can be reduced by a factor of ~10 relative to current constraints by future space-based supernova data and CMB measurements from the Planck satellite. If the dark energy equation of state is described by a more general basis of principal components, the expected improvement in volume-based figures of merit is much greater. While the forecasted precision for any single parameter is only a factor of 2-5 smaller than current uncertainties, the constraints on dark energy models bounded by -1<w<1 improve for approximately 6 independent dark energy parameters resulting in a reduction of the total allowed volume of principal component parameter space by a factor of ~100. Typical quintessence models can be adequately described by just 2-3 of these parameters even given the precision of future data, leading to a more modest but still significant improvement. In addition to advances in supernova and CMB data, percent-level measurement of absolute distance and/or the expansion rate is required to ensure that dark energy constraints remain robust to variations in spatial curvature.

preprint2010arXiv

Observational Limits on Patchy Reionization: Implications for B-modes

The recent detection of secondary CMB anisotropy by the South Pole Telescope places a conservative bound on temperature fluctuations from the optical depth-modulated Doppler effect of T_{3000} < sqrt{13} microK at multipoles l~3000. This bound is the first empirical constraint on reionization optical depth fluctuations at arcminute scales, tau_{3000} = 0.001 T_{3000}/microK, implying that these fluctuations are no more than a few percent of the mean. Optical depth modulation of the quadrupole source to polarization generates B-modes that are correspondingly bounded as B_{3000} = 0.003 T_{3000}. The maximal extrapolation to the l~100 gravitational wave regime yields B_{100} = 0.1 T_{3000} and remains in excess of gravitational lensing if the effective comoving size of the ionizing regions is R > 80 Mpc. If patchy reionization is responsible for much of the observed arcminute scale temperature fluctuations, current bounds on B_{100} already require R < 200 Mpc and can be expected to improve rapidly. Frequency separation of thermal Sunyaev-Zel'dovich contributions to the measured secondary anisotropy would also substantially improve the limits on optical depth fluctuations and B-modes from reionization.

preprint2010arXiv

Spherical Collapse and the Halo Model in Braneworld Gravity

We present a detailed study of the collapse of a spherical perturbation in DGP braneworld gravity for the purpose of modeling simulation results for the halo mass function, bias and matter power spectrum. The presence of evolving modifications to the gravitational force in form of the scalar brane-bending mode lead to qualitative differences to the collapse in ordinary gravity. In particular, differences in the energetics of the collapse necessitate a new, generalized method for defining the virial radius which does not rely on strict energy conservation. These differences and techniques apply to smooth dark energy models with w unequal -1 as well. We also discuss the impact of the exterior of the perturbation on collapse quantities due to the lack of a Birkhoff theorem in DGP. The resulting predictions for the mass function, halo bias and power spectrum are in good overall agreement with DGP N-body simulations on both the self-accelerating and normal branch. In particular, the impact of the Vainshtein mechanism as measured in the full simulations is matched well. The model and techniques introduced here can serve as practical tools for placing consistent constraints on braneworld models using observations of large scale structure.

preprint2009arXiv

Cosmological Constraints on DGP Braneworld Gravity with Brane Tension

We perform a Markov Chain Monte Carlo analysis of the self-accelerating and normal branch of Dvali-Gabadadze-Porrati braneworld gravity. By adopting a parameterized post-Friedmann description of gravity, we utilize all of the cosmic microwave background data, including the largest scales, and its correlation with galaxies in addition to the geometrical constraints from supernovae distances and the Hubble constant. We find that on both branches brane tension or a cosmological constant is required at high significance with no evidence for the unique Dvali-Gabadadze-Porrati modifications. The cross-over scale must therefore be substantially greater than the Hubble scale H_0 r_c > 3 and 3.5 at the 95% CL with and without uncertainties from spatial curvature. With spatial curvature, the limit from the normal branch is substantially assisted by the galaxy cross-correlation which highlights its importance in constraining infrared modifications to gravity.

preprint2009arXiv

Generalized Slow Roll for Large Power Spectrum Features

We develop a variant of the generalized slow roll approach for calculating the curvature power spectrum that is well-suited for order unity deviations in power caused by sharp features in the inflaton potential. As an example, we show that predictions for a step function potential, which has been proposed to explain order unity glitches in the CMB temperature power spectrum at multipoles l=20-40, are accurate at the percent level. Our analysis shows that to good approximation there is a single source function that is responsible for observable features and that this function is simply related to the local slope and curvature of the inflaton potential. These properties should make the generalized slow roll approximation useful for inflation-model independent studies of features, both large and small, in the observable power spectra.

preprint2009arXiv

Hiding dark energy transitions at low redshift

We show that it is both observationally allowable and theoretically possible to have large fluctuations in the dark energy equation of state as long as they occur at ultra-low redshifts z<0.02. These fluctuations would masquerade as a local transition in the Hubble rate of a few percent or less and escape even future, high precision, high redshift measurements of the expansion history and structure. Scalar field models that exhibit this behavior have a sharp feature in the potential that the field traverses within a fraction of an e-fold of the present. The equation of state parameter can become arbitrarily large if a sharp dip or bump in the potential causes the kinetic and potential energy of the field to both be large and have opposite sign. While canonical scalar field models can decrease the expansion rate at low redshift, increasing the local expansion rate requires a non-canonical kinetic term for the scalar field.

preprint2009arXiv

Testable dark energy predictions from current data

Given a class of dark energy models, constraints from one set of cosmic acceleration observables make predictions for other observables. Here we present the allowed ranges for the expansion rate H(z), distances D(z), and the linear growth function G(z) (as well as other, derived growth observables) from the current combination of cosmological measurements of supernovae, the cosmic microwave background, baryon acoustic oscillations, and the Hubble constant. With a cosmological constant as the dark energy and assuming near-minimal neutrino masses, the growth function is already predicted to better than 2% precision at any redshift, with or without spatial curvature. Direct measurements of growth that match this precision offer the opportunity to stringently test and potentially rule out a cosmological constant. While predictions in the broader class of quintessence models are weaker, it is remarkable that they are typically within a factor of 2-3 of forecasts for future space-based supernovae and Planck CMB measurements. In particular, measurements of growth at any redshift, or the Hubble constant H_0, that exceed LambdaCDM predictions by substantially more than 2% would rule out not only a cosmological constant but also the whole quintessence class, with or without curvature and early dark energy. Barring additional systematic errors hiding in the data, such a discovery would require more exotic explanations of cosmic acceleration such as phantom dark energy, dark energy clustering, or modifications of gravity.

preprint2008arXiv

Challenges to the DGP Model from Horizon-Scale Growth and Geometry

We conduct a Markov Chain Monte Carlo study of the Dvali-Gabadadze-Porrati (DGP) self-accelerating braneworld scenario given the cosmic microwave background (CMB) anisotropy, supernovae and Hubble constant data by implementing an effective dark energy prescription for modified gravity into a standard Einstein-Boltzmann code. We find no way to alleviate the tension between distance measures and horizon scale growth in this model. Growth alterations due to perturbations propagating into the bulk appear as excess CMB anisotropy at the lowest multipoles. In a flat cosmology, the maximum likelihood DGP model is nominally a 5.3 sigma poorer fit than Lambda CDM. Curvature can reduce the tension between distance measures but only at the expense of exacerbating the problem with growth leading to a 4.8 sigma result that is dominated by the low multipole CMB temperature spectrum. While changing the initial conditions to reduce large scale power can flatten the temperature spectrum, this also suppresses the large angle polarization spectrum in violation of recent results from WMAP5. The failure of this model highlights the power of combining growth and distance measures in cosmology as a test of gravity on the largest scales.

preprint2008arXiv

Falsifying Paradigms for Cosmic Acceleration

Consistency relations between growth of structure and expansion history observables exist for any physical explanation of cosmic acceleration, be it a cosmological constant, scalar field quintessence, or a general component of dark energy that is smooth relative to dark matter on small scales. The high-quality supernova sample anticipated from an experiment like SNAP and CMB data expected from Planck thus make strong predictions for growth and expansion observables that additional observations can test and potentially falsify. We perform an MCMC likelihood exploration of the strength of these consistency relations based on a complete parametrization of dark energy behavior by principal components. For LCDM, future SN and CMB data make percent level predictions for growth and expansion observables. For quintessence, many of the predictions are still at a level of a few percent with most of the additional freedom coming from curvature and early dark energy. While such freedom is limited for quintessence where phantom equations of state are forbidden, it is larger in the smooth dark energy class. Nevertheless, even in this general class predictions relating growth measurements at different redshifts remain robust, although predictions for the instantaneous growth rate do not. Finally, if observations falsify the whole smooth dark energy class, new paradigms for cosmic acceleration such as modified gravity or interacting dark matter and dark energy would be required.

preprint2007arXiv

Cluster Mass Estimators from CMB Temperature and Polarization Lensing

Upcoming Sunyaev-Zel'dovich surveys are expected to return ~10^4 intermediate mass clusters at high redshift. Their average masses must be known to same accuracy as desired for the dark energy properties. Internal to the surveys, the CMB potentially provides a source for lensing mass measurements whose distance is precisely known and behind all clusters. We develop statistical mass estimators from 6 quadratic combinations of CMB temperature and polarization fields that can simultaneously recover large-scale structure and cluster mass profiles. The performance of these estimators on idealized NFW clusters suggests that surveys with a ~1' beam and 10uK' noise in uncontaminated temperature maps can make a ~10sigma detection, or equivalently a ~10% mass measurement for each 10^3 set of clusters. With internal or external acoustic scale E-polarization measurements, the ET cross correlation estimator can provide a stringent test for contaminants on a first detection at \~1/3 the significance. For surveys that reach below 3muK', the EB cross correlation estimator should provide the most precise measurements and potentially the strongest control over contaminants.

preprint2007arXiv

Model-independent constraints on reionization from large-scale CMB polarization

On large angular scales, the polarization of the CMB contains information about the evolution of the average ionization during the epoch of reionization. Interpretation of the polarization spectrum usually requires the assumption of a fixed functional form for the evolution, e.g. instantaneous reionization. We develop a model-independent method where a small set of principal components completely encapsulate the effects of reionization on the large-angle E-mode polarization for any reionization history within an adjustable range in redshift. Using Markov Chain Monte Carlo methods, we apply this approach to both the 3-year WMAP data and simulated future data. WMAP data constrain two principal components of the reionization history, approximately corresponding to the total optical depth and the difference between the contributions to the optical depth at high and low redshifts. The optical depth is consistent with the constraint found in previous analyses of WMAP data that assume instantaneous reionization, with only slightly larger uncertainty due to the expanded set of models. Using the principal component approach, WMAP data also place a 95% CL upper limit of 0.08 on the contribution to the optical depth from redshifts z>20. With improvements in polarization sensitivity and foreground modeling, approximately five of the principal components can ultimately be measured. Constraints on the principal components, which probe the entire reionization history, can test models of reionization, provide model-independent constraints on the optical depth, and detect signatures of high-redshift reionization.

preprint2006arXiv

Likelihood Methods for Cluster Dark Energy Surveys

Galaxy cluster counts at high redshift, binned into spatial pixels and binned into ranges in an observable proxy for mass, contain a wealth of information on both the dark energy equation of state and the mass selection function required to extract it. The likelihood of the number counts follows a Poisson distribution whose mean fluctuates with the large-scale structure of the universe. We develop a joint likelihood method that accounts for these distributions. Maximization of the likelihood over a theoretical model that includes both the cosmology and the observable-mass relations allows for a joint extraction of dark energy and cluster structural parameters.

preprint2005arXiv

Constraints on supernovae dimming from photon-pseudo scalar coupling

An alternative mechanism that dims high redshift supernovae without cosmic acceleration utilizes an oscillation of photons into a pseudo-scalar particle during transit. Since angular diameter distance measures are immune to the loss of photons, this ambiguity in interpretation can be resolved by combining CMB acoustic peak measurements with the recent baryon oscillation detection in galaxy power spectra. This combination excludes a non-accelerating dark energy species at the 4sigma level regardless of the level of the pseudo-scalar coupling. While solutions still exist with substantial non-cosmological dimming of supernovae, they may be tested with future improvement in baryon oscillation experiments.

preprint2005arXiv

Crossing the Phantom Divide: Dark Energy Internal Degrees of Freedom

Dark energy constraints have forced viable alternatives that differ substantially from a cosmological constant Lambda to have an equation of state w that evolves across the phantom divide set by Lambda. Naively, crossing this divide makes the dark energy gravitationally unstable, a problem that is typically finessed by unphysically ignoring the perturbations. While this procedure does not affect constraints near the favored cosmological constant model it can artificially enhance the confidence with which alternate models are rejected. Similar to the general problem of stability for w< 0, the solution lies in the internal degrees of freedom in the dark energy sector. We explicitly show how to construct a two scalar field model that crosses the phantom divide and mimics the single field behavior on either side to substantially better than 1% in all observables. It is representative of models where the internal degrees of freedom keep the dark energy smooth out to the horizon scale independently of the equation of state.

preprint2005arXiv

Self-Calibration of Cluster Dark Energy Studies: Observable-Mass Distribution

The exponential sensitivity of cluster number counts to the properties of the dark energy implies a comparable sensitivity to not only the mean but also the actual_distribution_ of an observable mass proxy given the true cluster mass. For example a 25% scatter in mass can provide a ~50% change in the number counts at z~2 for the upcoming SPT survey. Uncertainty in the scatter of this amount would degrade dark energy constraints to uninteresting levels. Given the shape of the actual mass function, the properties of the distribution may be internally monitored by the shape of the_observable_ mass function. An arbitrary evolution of the scatter of a mass-independent Gaussian distribution may be self-calibrated to allow a measurement of the dark energy equation of state of Delta w ~0.1. External constraints on the mass_variance_ of the distribution that are more accurate than Delta var < 0.01 at z~1 can further improve constraints by up to a factor of 2. More generally, cluster counts and their sample variance measured as a function of the observable provide internal consistency checks on the assumed form of the observable-mass distribution that will protect against misinterpretation of the dark energy constraints.

preprint2005arXiv

Weighing Neutrinos with Galaxy Cluster Surveys

Large future galaxy cluster surveys, combined with cosmic microwave background observations, can achieve a high sensitivity to the masses of cosmologically important neutrinos. We show that a weak lensing selected sample of ~100,000 clusters could tighten the current upper bound on the sum of masses of neutrino species by an order of magnitude, to a level of 0.03 eV. Since this statistical sensitivity is below the best existing lower limit on the mass of at least one neutrino species, a future detection is likely, provided that systematic errors can be controlled to a similar level.

preprint2004arXiv

A Low CMB Quadrupole from Dark Energy Isocurvature Perturbations

We explicate the origin of the temperature quadrupole in the adiabatic dark energy model and explore the mechanism by which scale invariant isocurvature dark energy perturbations can lead to its sharp suppression. The model requires anticorrelated curvature and isocurvature fluctuations and is favored by the WMAP data at about the 95% confidence level in a flat scale invariant model. In an inflationary context, the anticorrelation may be established if the curvature fluctuations originate from a variable decay rate of the inflaton; such models however tend to overpredict gravitational waves. This isocurvature model can in the future be distinguished from alternatives involving a reduction in large scale power or modifications to the sound speed of the dark energy through the polarization and its cross correlation with the temperature. The isocurvature model retains the same polarization fluctuations as its adiabatic counterpart but reduces the correlated temperature fluctuations. We present a pedagogical discussion of dark energy fluctuations in a quintessence and k-essence context in the Appendix.

preprint2004arXiv

Dark Energy Probes in Light of the CMB

CMB observables have largely fixed the expansion history of the universe in the deceleration regime and provided two self-calibrated absolute standards for dark energy studies: the sound horizon at recombination as a standard ruler and the amplitude of initial density fluctuations. We review these inferences and expose the testable assumptions about recombination and reionization that underly them. Fixing the deceleration regime with CMB observables, deviations in the distance and growth observables appear most strongly at z=0 implying that accurate calibration of local and CMB standards may be more important than redshift range or depth. The single most important complement to the CMB for measuring the dark energy equation of state at z~0.5 is a determination of the Hubble constant to better than a few percent. Counterintuitively, with fixed fractional distance errors and relative standards such as SNe, the Hubble constant measurement is best achieved in the high redshift deceleration regime. Degeneracies between the evolution and current value of the equation of state or between its value and spatial curvature can be broken if percent level measurement and calibration of distance standards can be made at intermediate redshifts or the growth function at any redshift in the acceleration regime. We compare several dark energy probes available to a wide and deep optical survey: baryon features in galaxy angular power spectra and the growth rate from galaxy-galaxy lensing, shear tomography and the cluster abundance.

preprint2004arXiv

Measuring Dark Energy Clustering with CMB-Galaxy Correlations

The integrated Sachs-Wolfe (ISW) effect in the cosmic microwave background (CMB) as measured through its correlation with galaxies provides a unique opportunity to study the dynamics of the dark energy through its large scale clustering properties. Ultimately, a deep all-sky galaxy survey out to z~2 can make a 10sigma or ~10% measurement of the correlation and limit ~3% changes in the gravitational potential or total density fluctuation due to dark energy clustering on the Gpc scale. A canonical single scalar field or quintessence model predicts that these clustering effects will appear on the horizon scale with a strength that reflects the evolution of the dark energy density. In terms of a constant equation of state, this would allow tests of the quintessence prediction for models where |1+w| > 0.05.

preprint2004arXiv

Weak Lensing of the CMB: Sampling Errors on B-Modes

The B modes generated by the lensing of CMB polarization are a primary target for the upcoming generation of experiments and can potentially constrain quantities such as the neutrino mass and dark energy equation of state. The net sample variance on the small scale B modes out to l=2000 exceeds Gaussian expectations by a factor of 10 reflecting the variance of the larger scale lenses that generate them. It manifests itself as highly correlated band powers with correlation coefficients approaching 70% for wide bands of Delta l/l \~0.25. It will double the total variance for experiments that achieve a sensitivity of approximately 4 uK-arcmin and a beam of several arcminutes or better. This non-Gaussianity must be taken into account in the analysis of experiments that go beyond first detection.

preprint2003arXiv

CMB Lensing Reconstruction on the Full Sky

Gravitational lensing of the microwave background by the intervening dark matter mainly arises from large-angle fluctuations in the projected gravitational potential and hence offers a unique opportunity to study the physics of the dark sector at large scales. Studies with surveys that cover greater than a percent of the sky will require techniques that incorporate the curvature of the sky. We lay the groundwork for these studies by deriving the full sky minimum variance quadratic estimators of the lensing potential from the CMB temperature and polarization fields. We also present a general technique for constructing these estimators, with harmonic space convolutions replaced by real space products, that is appropriate for both the full sky limit and the flat sky approximation. This also extends previous treatments to include estimators involving the temperature-polarization cross-correlation and should be useful for next generation experiments in which most of the additional information from polarization comes from this channel due to sensitivity limitations.

preprint2003arXiv

Model-Independent Reionization Observables in the CMB

We represent the reionization history of the universe as a free function in redshift and study the potential for its extraction from CMB polarization spectra. From a principal component analysis, we show that the ionization history information is contained in 5 modes, resembling low-order Fourier modes in redshift space. The amplitude of these modes represent a compact description of the observable properties of reionization in the CMB, easily predicted given a model for the ionization fraction. Measurement of these modes can ultimately constrain the total optical depth, or equivalently the initial amplitude of fluctuations to the 1% level regardless of the true model for reionization.

preprint2003arXiv

Self-Consistency and Calibration of Cluster Number Count Surveys for Dark Energy

Cluster number counts offer sensitive probes of the dark energy if and only if the_evolution_ of the cluster mass versus observable relation(s) is well calibrated. We investigate the potential for internal calibration by demanding consistency in the counts as a function of the observable. In the context of a constant dark energy equation of state, known initial fluctuation amplitude expected from the CMB, universal underlying mass function, and an idealized selection, we find that the ambiguity from the normalization of the mass-observable relationships, or an extrapolation of external mass-observable determinations from higher masses, can be largely eliminated with a sufficiently deep survey, even allowing for an arbitrary evolution. More generally, number counts as a function of both the redshift and the observable enable strong consistency tests on assumptions made in modelling the mass-observable relations and cosmology.

preprint2002arXiv

Benchmark Parameters for CMB Polarization Experiments

The recently detected polarization of the cosmic microwave background (CMB) holds the potential for revealing the physics of inflation and gravitationally mapping the large-scale structure of the universe, if so called B-mode signals below 10^{-7}, or tenths of a uK, can be reliably detected. We provide a language for describing systematic effects which distort the observed CMB temperature and polarization fields and so contaminate the B-modes. We identify 7 types of effects, described by 11 distortion fields, and show their association with known instrumental systematics such as common mode and differential gain fluctuations, line cross-coupling, pointing errors, and differential polarized beam effects. Because of aliasing from the small-scale structure in the CMB, even uncorrelated fluctuations in these effects can affect the large-scale B modes relevant to gravitational waves. Many of these problems are greatly reduced by having an instrumental beam that resolves the primary anisotropies (FWHM << 10'). To reach the ultimate goal of an inflationary energy scale of 3 \times 10^{15} GeV, polarization distortion fluctuations must be controlled at the 10^{-2}-10^{-3} level and temperature leakage to the 10^{-4}-10^{-3} level depending on effect. For example pointing errors must be controlled to 1.5'' rms for arcminute scale beams or a percent of the Gaussian beam width for larger beams; low spatial frequency differential gain fluctuations or line cross-coupling must be eliminated at the level of 10^{-4} rms.

preprint2002arXiv

CMB Temperature and Polarization Anisotropy Fundamentals

The tremendous experimental progress in cosmic microwave background (CMB) temperature and polarization anisotropy studies over the last few years has helped establish a standard paradigm for cosmology at intermediate epochs and has simultaneously raised questions regarding the physical processes at the two opposite ends of time. What is the physics behind the source of structure in the universe and the dark energy that is currently accelerating its expansion? We review the acoustic phenomenology that forms the cornerstone of the standard cosmological model and discuss internal consistency relations which lend credence to its interpretation. We discuss future milestones in the study of CMB anisotropy that have implications for inflationary and dark energy models. These include signatures of gravitational waves and gravitational lensing in the polarization fields and secondary temperature anisotropy from the transit of CMB photons across the large-scale structure of the universe.

preprint2002arXiv

Sample Variance Considerations for Cluster Surveys

We present a general statistical framework for describing the effect of sample variance in the number counts of virialized objects and examine its effect on cosmological parameter estimation. Specifically, we consider effects of sample variance on the power spectrum normalization and properties of dark energy extracted from current and future local and high-redshift samples of clusters. We show that for future surveys that probe ever lower cluster masses and temperatures, sample variance is generally comparable to or greater than shot noise and thus cannot be neglected in deriving precision cosmological constraints. For example, sample variance is usually more important than shot variance in constraints on the equation of state of the dark energy from z < 1 clusters. Although we found that effects of sample variance on the sigma_8-Omega_m constraints from the current flux and temperature limited X-ray surveys are not significant, they may be important for future studies utilizing the shape of the temperature function to break the sigma_8-Omega_m degeneracy. We also present numerical tests clarifying the definition of cluster mass employed in cosmological modelling and accurate fitting formula for the conversion between different definitions of halo mass (e.g., virial vs. fixed overdensity).

preprint2002arXiv

Second Order Corrections to Weak Lensing by Large-Scale Structure

We calculate corrections to the power spectrum predictions of weak lensing by large scale structure due to higher order effects in the gravitational potential. Using a perturbative approach to third order in transverse displacements, we calculate a second order correction to the angular power spectra of E and B mode shear and convergence resulting from dropping the so-called Born approximation, where one integrates along the unperturbed photon path. We also consider a correction to the power spectra from the coupling between lenses at different redshifts. Both effects generate B-mode shear and the latter also causes a net rotation of the background galaxy images. We show all these corrections are at least two orders of magnitude below the convergence or E-mode power and hence relevant only to future ultra high precision measurements. These analytical calculations are consistent with previous numerical estimates and validate the use of current large scale structure weak lensing predictions for cosmological studies and future use of B-modes as a monitor of systematic effects.

preprint2002arXiv

The Angular Trispectra of CMB Temperature and Polarization

We develop the formalism necessary to study four-point functions of the cosmic microwave background (CMB) temperature and polarization fields. We determine the general form of CMB trispectra, with the constraints imposed by the assumption of statistical isotropy of the CMB fields, and derive expressions for their estimators, as well as their Gaussian noise properties. We apply these techniques to initial non-Gaussianity of a form motivated by inflationary models. Due to the large number of four-point configurations, the sensitivity of the trispectra to initial non-Gaussianity approaches that of the temperature bispectrum at high multipole moment. These trispectra techniques will also be useful in the study of secondary anisotropies induced for example by the gravitational lensing of the CMB by the large scale structure of the universe.

preprint2001arXiv

Cosmic Microwave Background Anisotropies

Cosmic microwave background (CMB) temperature anisotropies have and will continue to revolutionize our understanding of cosmology. The recent discovery of the previously predicted acoustic peaks in the power spectrum has established a working cosmological model: a critical density universe consisting of mainly dark matter and dark energy, which formed its structure through gravitational instability from quantum fluctuations during an inflationary epoch. Future observations should test this model and measure its key cosmological parameters with unprecedented precision. The phenomenology and cosmological implications of the acoustic peaks are developed in detail. Beyond the peaks, the yet to be detected secondary anisotropies and polarization present opportunities to study the physics of inflation and the dark energy. The analysis techniques devised to extract cosmological information from voluminous CMB data sets are outlined, given their increasing importance in experimental cosmology as a whole.

preprint1999arXiv

A Model for Structure Formation Seeded by Gravitationally Produced Matter

This model assumes the baryons, radiation, three families of massless neutrinos, and cold dark matter were mutually thermalized before the baryon number was fixed, primeval curvature fluctuations were subdominant, and homogeneity was broken by scale-invariant fluctuations in a new dark matter component that behaves like a relativistic ideal fluid. The fluid behavior could follow if this new component were a single scalar field that interacts only with gravity and with itself by a pure quartic potential. The initial energy distribution could follow if this component were gravitationally produced by inflation. The power spectra of the present distributions of mass and radiation in this model are not inconsistent with the measurements but are sufficiently different from the adiabatic cold dark matter model to allow a sharp test in the near future.

preprint1998arXiv

An Isocurvature Mechanism for Structure Formation

We examine a novel mechanism for structure formation involving initial number density fluctuations between relativistic species, one of which then undergoes a temporary downward variation in its equation of state and generates superhorizon-scale density fluctuations. Isocurvature decaying dark matter models (iDDM) provide concrete examples. This mechanism solves the phenomenological problems of traditional isocurvature models, allowing iDDM models to fit the current CMB and large-scale structure data, while still providing novel behavior. We characterize the decaying dark matter and its decay products as a single component of ``generalized dark matter''. This simplifies calculations in decaying dark matter models and others that utilize this mechanism for structure formation.

preprint1998arXiv

Cosmic Complementarity: H_0 and Omega_m from Combining CMB Experiments and Redshift Surveys

We show that upcoming CMB satellite experiments and large redshift surveys can be used together to yield 5% determinations of H_0 and Omega_m, an order of magnitude improvement over CMB data alone. CMB anisotropies provide the sound horizon at recombination as a standard ruler. For reasonable baryon fractions, this scale is imprinted on the galaxy power spectrum as a series of spectral features. Measuring these features in redshift space determines the Hubble constant, which in turn yields Omega_m once combined with CMB data. Since the oscillations in both power spectra are frozen in at recombination, this test is insensitive to low-redshift cosmology.

preprint1998arXiv

The Structure of Structure Formation Theories

We study the general structure of models for structure formation, with applications to the reverse engineering of the model from observations. Through a careful accounting of the degrees of freedom in covariant gravitational instability theory, we show that the evolution of structure is completely specified by the stress history of the dark sector. The study of smooth, entropic, sonic, scalar anisotropic, vector anisotropic, and tensor anisotropic stresses reveals the origin, robustness, and uniqueness of specific model phenomenology. We construct useful and illustrative analytic solutions that cover cases with multiple species of differing equations of state relevant to the current generation of models, especially those with effectively smooth components. We present a simple case study of models with phenomenologies similar to that of a LambdaCDM model to highlight reverse-engineering issues. A critical-density universe dominated by a single type of dark matter with the appropriate stress history can mimic a LambdaCDM model exactly.

preprint1997arXiv

Baryonic Features in the Matter Transfer Function

We provide scaling relations and fitting formulae for adiabatic cold dark matter cosmologies that account for all baryon effects in the matter transfer function to better than 10% in the large-scale structure regime. They are based upon a physically well-motivated separation of the effects of acoustic oscillations, Compton drag, velocity overshoot, baryon infall, adiabatic damping, Silk damping, and cold-dark-matter growth suppression. We also find a simpler, more accurate, and better motivated form for the zero baryon transfer function than previous works. These descriptions are employed to quantify the amplitude and location of baryonic features in linear theory. While baryonic oscillations are prominent if the baryon fraction exceeds $Ω_0 h^2 + 0.2$, the main effect in more conventional cosmologies is a sharp suppression in the transfer function below the sound horizon. We provide a simple but accurate description of this effect and stress that it is not well approximated by a change in the shape parameter $Γ$.

preprint1995arXiv

COBE Constraints on Baryon Isocurvature Models

We consider Bayesian constraints on standard isocurvature baryon models from the slope and normalization of the anisotropy power spectrum detected by the COBE DMR experiment in their two year maps. In conjunction with either the amplitude of matter fluctuations $σ_8$ {\it or} its slope, all open models are ruled out at greater than 95\% confidence, whereas cosmological constant $Λ$ dominated models are constrained to be highly ionized. By including the COBE FIRAS 95\% confidence upper limit on spectral distortion, we further reduce the available parameter space for $Λ$ models by excluding these highly ionized models. These constraints define a single class of barely viable standard models which makes definite and testable predictions for degree scale anisotropies and large scale structure.

preprint1995arXiv

Concepts in CMB Anisotropy Formation

These lecture notes form a primer on the theory of cosmic microwave background (CMB) anisotropy formation. With emphasis on conceptual aspects rather than technical issues, we examine the physical foundations of anisotropy evolution in relativistic kinetic and perturbation theory as well as the manifestation of these principles in primary and secondary anisotropies. We discuss gauge choice and gauge invariance and their use in understanding the CMB. Acoustic, gravitational redshift and ionization effects have robust signatures in the CMB spectrum and may allow determination of classical cosmological parameters as well as reveal general distinctions between models for structure formation. We develop the tight and weak coupling approximations as analytic tools to help understand these effects and the robustness of their signatures.