Source author record

Richard Easther

Richard Easther 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

54works
9topics
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

54 published item(s)

preprint2023arXiv

Constraining Cosmic Inflation with Observations: Prospects for 2030

The ability to test and constrain theories of cosmic inflation will advance substantially over the next decade. Key data sources include cosmic microwave background (CMB) measurements and observations of the distribution of matter at low-redshift from optical, near-infrared, and 21cm intensity surveys. A positive detection of a CMB B-mode consistent with a primordial stochastic gravitational wave background (SGWB) is widely viewed as a smoking gun for an inflationary phase. Still, a null result does not exclude inflation. However, in a significant class of inflationary scenarios, a low SGWB amplitude is correlated with a more significant running, $α_s$, in the primordial density perturbations than is seen with the simplest inflationary potentials. With this motivation, we forecast the precision with which the spectral index $n_{\rm{s}}$ and $α_{\rm{s}}$ can be constrained by currently envisaged observations, including CMB (Simons Observatory, CMB-S4 and LiteBIRD), optical/near infra-red (DESI and SPHEREx), and 21cm intensity mapping (Tianlai and CHIME) surveys. We identify optimal combinations of datasets for constraining the running and show that they may yield additional and informative constraints on the overall inflationary parameter space if the SGWB remains undetected.

preprint2022arXiv

Dynamical Friction From Ultralight Dark Matter

We simulate the gravitational dynamics of a massive object interacting with Ultralight / Fuzzy Dark Matter (ULDM/FDM), non-relativistic quantum matter described by the Schrodinger-Poisson equation. We first consider a point mass moving in a uniform background, and then a supermassive black hole (SMBH) moving within a ULDM soliton. After replicating simple dynamical friction scenarios to verify our numerical strategies, we demonstrate that the wake induced by a moving mass in a uniform medium may undergo gravitational collapse that dramatically increases the drag force, albeit in a scenario unlikely to be encountered astrophysically. We broadly confirm simple estimates of dynamical friction timescales for a black hole at the center of a halo but see that a large moving point mass excites coherent "breathing modes" in a ULDM soliton. These can lead to "stone skipping" trajectories for point masses which do not sink uniformly toward the center of the soliton, as well as stochastic motion near the center itself. These effects will add complexity to SMBH-ULDM interactions and to SMBH mergers in a ULDM universe.

preprint2022arXiv

Gravitational Collapse in the Post-Inflationary Universe

The Universe may pass through an effectively matter-dominated epoch between inflation and Big Bang Nucleosynthesis during which gravitationally bound structures can form on subhorizon scales. In particular, the inflaton field can collapse into inflaton halos, forming "large scale" structure in the very early universe. We combine N-body simulations with high-resolution zoom-in regions in which the non-relativistic Schrödinger-Poisson equations are used to resolve the detailed, wave-like structure of inflaton halos. Solitonic cores form inside them, matching structure formation simulations with axion-like particles in the late-time universe. We denote these objects \textit{inflaton stars}, by analogy with boson stars. Based on a semi-analytic formalism we compute their overall mass distribution which shows that some regions will reach overdensities of $10^{15}$ if the early matter-dominated epoch lasts for 20 $e$-folds. The radii of the most massive inflaton stars can shrink below the Schwarzschild radius, suggesting that they could form primordial black holes prior to thermalization.

preprint2022arXiv

New Horizons for Fundamental Physics with LISA

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas.

preprint2022arXiv

Schrödinger-Poisson Solitons: Perturbation Theory

Self-gravitating quantum matter may exist in a wide range of cosmological and astrophysical settings from the very early universe through to present-day boson stars. Such quantum matter arises in a number of different theories, including the Peccei-Quinn axion and UltraLight (ULDM) or Fuzzy (FDM) dark matter scenarios. We consider the dynamical evolution of perturbations to the spherically symmetric soliton, the ground state solution to the Schrödinger-Poisson system common to all these scenarios. We construct the eigenstates of the Schrödinger equation, holding the gravitational potential fixed to its ground state value. We see that the eigenstates qualitatively capture the properties seen in full ULDM simulations, including the soliton "breathing" mode, the random walk of the soliton center, and quadrupolar distortions of the soliton. We then show that the time-evolution of the gravitational potential and its impact on the perturbations can be well described within the framework of time-dependent perturbation theory. Applying our formalism to a synthetic ULDM halo reveals considerable mixing of eigenstates, even though the overall density profile is relatively stable. Our results provide a new analytic approach to understanding the evolution of these systems as well as possibilities for faster approximate simulations.

preprint2020arXiv

Inflaton Clusters and Inflaton Stars

In a broad class of scenarios, inflation is followed by an extended era of matter-dominated expansion during which the inflaton condensate is nonrelativistic on subhorizon scales. During this phase density perturbations grow to the point of nonlinearity and collapse into bound structures. This epoch strongly resembles structure formation with ultra-light axion-like particles. This parallel permits us to adapt results from studies of cosmological structure formation to describe the nonlinear dynamics of this post-inflationary epoch. We show that the inflaton condensate fragments into "inflaton clusters", analogues of axion dark matter halos in present-day cosmology. Moreover, solitonic objects or "inflaton stars" can form inside these clusters, leading to density contrasts as large as $10^6$ in the post-inflationary universe.

preprint2020arXiv

Robust posterior inference when statistically emulating forward simulations

Scientific analyses often rely on slow, but accurate forward models for observable data conditioned on known model parameters. While various emulation schemes exist to approximate these slow calculations, these approaches are only safe if the approximations are well understood and controlled. This workshop submission reviews and updates a previously published method, which has been used in cosmological simulations, to (1) train an emulator while simultaneously estimating posterior probabilities with MCMC and (2) explicitly propagate the emulation error into errors on the posterior probabilities for model parameters. We demonstrate how these techniques can be applied to quickly estimate posterior distributions for parameters of the $Λ$CDM cosmology model, while also gauging the robustness of the emulator approximation.

preprint2020arXiv

The Core-Cusp Problem Revisited: ULDM vs. CDM

The core-cusp problem is a widely cited motivation for the exploration of dark matter models beyond standard CDM. One such alternative is ULDM; extremely light scalar particles exhibiting wavelike properties on kiloparsec scales. Astrophysically realistic ULDM halos are expected to consist of inner solitonic cores embedded in NFW-like outer halos. The presence of the solitonic core suggests that ULDM may resolve the core-cusp discrepancy associated with pure NFW halos without recourse to baryonic physics. However, it has been demonstrated that the density of ULDM halos can exceed those of comparable NFW configurations at some radii and halo masses, apparently exacerbating the problem rather than solving it. This situation arises because, although solitonic cores are flat at their centres, they obey an inverse mass-radius scaling relationship. Meanwhile, the mass of the inner soliton increases with the total halo mass, and therefore the inner core becomes more peaked at large halo masses. We describe a parameterisation of the radial density profiles of ULDM halos that allows for environmental variability of the core-halo mass relation in order to investigate this issue in more detail. For halos up to 10^12 solar masses we find feasible ULDM profiles for which the central density is lower than their NFW counterparts at astrophysically accessible radii. However, comparisons to observed profiles do not strongly favour either option; both give reasonable fits to subsets of the data for some parameter choices. Consequently, we find that robust tests of the core-cusp problem in ULDM will require more comprehensive observational data and simulations that include baryonic feedback.

preprint2020arXiv

The Distribution of Vacua in Random Landscape Potentials

Landscape cosmology posits the existence of a convoluted, multidimensional, scalar potential -- the "landscape" -- with vast numbers of metastable minima. Random matrices and random functions in many dimensions provide toy models of the landscape, allowing the exploration of conceptual issues associated with these scenarios. We compute the relative number and slopes of minima as a function of the vacuum energy $Λ$ in an $N$-dimensional Gaussian random potential, quantifying the associated probability density, $p(Λ)$. After normalisations $p(Λ)$ depends only on the dimensionality $N$ and a single free parameter $γ$, which is related to the power spectrum of the random function. For a Gaussian landscape with a Gaussian power spectrum, the fraction of positive minima shrinks super-exponentially with $N$; at $N=100$, $p(Λ>0) \approx 10^{-1197}$. Likewise, typical eigenvalues of the Hessian matrices reveal that the flattest approaches to typical minima grow flatter with $N$, while the ratio of the slopes of the two flattest directions grows with $N$. We discuss the implications of these results for both swampland and conventional anthropic constraints on landscape cosmologies. In particular, for parameter values when positive minima are extremely rare, the flattest approaches to minima where $Λ\approx 0$ are much flatter than for typical minima, increasingly the viability of quintessence solutions.

preprint2019arXiv

Lighting the Dark: The Evolution of the Post-Inflationary Universe

In simple inflationary cosmological scenarios the near-exponential growth can be followed by a long period in which the Universe is dominated by the oscillating inflaton condensate. The condensate is initially almost homogeneous, but perturbations grow gravitationally, eventually fragmenting the condensate if it is not disrupted more quickly by resonance or prompt reheating. We show that the gravitational fragmentation of the condensate is well-described by the Schrödinger-Poisson equations and use numerical solutions to show that large overdensities form quickly after the onset of nonlinearity. This is the first exploration of this phase of nonlinear dynamics in the very early universe, which can affect the detailed form of the inflationary power spectrum and the dark matter fraction when the dark sector is directly coupled to the inflaton.

preprint2016arXiv

Counting Vacua in Random Landscapes

It is speculated that the correct theory of fundamental physics includes a large landscape of states, which can be described as a potential which is a function of N scalar fields and some number of discrete variables. The properties of such a landscape are crucial in determining key cosmological parameters including the dark energy density, the stability of the vacuum, the naturalness of inflation and the properties of the resulting perturbations, and the likelihood of bubble nucleation events. We codify an approach to landscape cosmology based on specifications of the overall form of the landscape potential and illustrate this approach with a detailed analysis of the properties of N-dimensional Gaussian random landscapes. We clarify the correlations between the different matrix elements of the Hessian at the stationary points of the potential. We show that these potentials generically contain a large number of minima. More generally, these results elucidate how random function theory is of central importance to this approach to landscape cosmology, yielding results that differ substantially from those obtained by treating the matrix elements of the Hessian as independent random variables.

preprint2016arXiv

Designing and testing inflationary models with Bayesian networks

Even simple inflationary scenarios have many free parameters. Beyond the variables appearing in the inflationary action, these include dynamical initial conditions, the number of fields, and couplings to other sectors. These quantities are often ignored but cosmological observables can depend on the unknown parameters. We use Bayesian networks to account for a large set of inflationary parameters, deriving generative models for the primordial spectra that are conditioned on a hierarchical set of prior probabilities describing the initial conditions, reheating physics, and other free parameters. We use $N_f$--quadratic inflation as an illustrative example, finding that the number of $e$-folds $N_*$ between horizon exit for the pivot scale and the end of inflation is typically the most important parameter, even when the number of fields, their masses and initial conditions are unknown, along with possible conditional dependencies between these parameters.

preprint2016arXiv

Testing for New Physics: Neutrinos and the Primordial Power Spectrum

We test the sensitivity of neutrino parameter constraints from combinations of CMB and LSS data sets to the assumed form of the primordial power spectrum (PPS) using Bayesian model selection. Significantly, none of the tested combinations, including recent high-precision local measurements of $\mathrm{H}_0$ and cluster abundances, indicate a signal for massive neutrinos or extra relativistic degrees of freedom. For PPS models with a large, but fixed number of degrees of freedom, neutrino parameter constraints do not change significantly if the location of any features in the PPS are allowed to vary, although neutrino constraints are more sensitive to PPS features if they are known a priori to exist at fixed intervals in $\log k$. Although there is no support for a non-standard neutrino sector from constraints on both neutrino mass and relativistic energy density, we see surprisingly strong evidence for features in the PPS when it is constrained with data from Planck 2015, SZ cluster counts, and recent high-precision local measurements of $\mathrm{H}_0$. Conversely combining Planck with matter power spectrum and BAO measurements yields a much weaker constraint. Given that this result is sensitive to the choice of data this tension between SZ cluster counts, Planck and $\mathrm{H}_0$ measurements is likely an indication of unmodeled systematic bias that mimics PPS features, rather than new physics in the PPS or neutrino sector.

preprint2016arXiv

Ultracompact minihalos as probes of inflationary cosmology

Cosmological inflation generates primordial density perturbations on all scales, including those far too small to contribute to the cosmic microwave background. At these scales, isolated ultracompact minihalos of dark matter can form well before standard structure formation, if the perturbations have sufficient amplitude. Minihalos affect pulsar timing data and are potentially bright sources of gamma rays. The resulting constraints significantly extend the observable window of inflation in the presence of cold dark matter, coupling two of the key problems in modern cosmology.

preprint2015arXiv

Learn-As-You-Go Acceleration of Cosmological Parameter Estimates

Cosmological analyses can be accelerated by approximating slow calculations using a training set, which is either precomputed or generated dynamically. However, this approach is only safe if the approximations are well understood and controlled. This paper surveys issues associated with the use of machine-learning based emulation strategies for accelerating cosmological parameter estimation. We describe a learn-as-you-go algorithm that is implemented in the Cosmo++ code and (1) trains the emulator while simultaneously estimating posterior probabilities; (2) identifies unreliable estimates, computing the exact numerical likelihoods if necessary; and (3) progressively learns and updates the error model as the calculation progresses. We explicitly describe and model the emulation error and show how this can be propagated into the posterior probabilities. We apply these techniques to the Planck likelihood and the calculation of $Λ$CDM posterior probabilities. The computation is significantly accelerated without a pre-defined training set and uncertainties in the posterior probabilities are subdominant to statistical fluctuations. We have obtained a speedup factor of $6.5$ for Metropolis-Hastings and $3.5$ for nested sampling. Finally, we discuss the general requirements for a credible error model and show how to update them on-the-fly.

preprint2015arXiv

Signatures of the Very Early Universe: Inflation, Spatial Curvature and Large Scale Anomalies

A short inflationary phase may not erase all traces of the primordial universe. Associated observables include both spatial curvature and "anomalies" in the microwave background or large scale structure. The present curvature $Ω_{K,0}$ reflects the initial curvature, $Ω_{K,\mathrm{start}}$, and the angular size of anomalies depends on $k_\mathrm{start}$, the comoving horizon size at the onset of inflation. We estimate posteriors for $Ω_{K,\mathrm{start}}$ and $k_\mathrm{start}$ using current data and simulations, and show that if either quantity is measured to have a non-zero value, both are likely to be observable. Mappings from $Ω_{K,\mathrm{start}}$ and $k_\mathrm{start}$ to present-day observables depend strongly on the primordial equation of state; $Ω_{K,0}$ spans ten orders of magnitude for a given $Ω_{K,\mathrm{start}}$ while a simple and general relationship connects $Ω_{K,0}$ and $k_\mathrm{start}$. We show that current bounds on $Ω_{K,0}$ imply that if $k_\mathrm{start}$ is measurable, the curvature was already small when inflation began. Finally, since the energy density changes slowly during inflation, primordial gravitational wave constraints require that a short inflationary phase is preceded by a nontrivial pre-inflationary phase with critical implications for the expected value of $Ω_{K,\mathrm{start}}$.

preprint2014arXiv

Gravitational wave consistency relations for multifield inflation

We study the tensor spectral index $n_t$ and the tensor-to-scalar ratio $r$ in the simplest multifield extension to single-field, slow-roll inflation models. We show that multifield models with potentials $V \sim \sum_i λ_i |ϕ_i|^p$ have different predictions for $n_t/r$ than single-field models, even when all the couplings are equal $λ_i=λ_j$, due to the probabilistic nature of the fields' initial values. We analyze well-motivated prior probabilities for the $λ_i$ and initial conditions to make detailed predictions for the marginalized probability distribution of $n_t/r$. With $\mathcal O(100)$ fields and $p>3/4$, we find that $n_t/r$ differs from the single-field result of $n_t/r=-1/8$ at the 5$σ$ level. This gives a novel and testable prediction for the simplest multifield inflation models.

preprint2014arXiv

Inflating an Inhomogeneous Universe

While cosmological inflation can erase primordial inhomogeneities, it is possible that inflation may not begin in a significantly inhomogeneous universe. This issue is particularly pressing in multifield scenarios, where even the homogeneous dynamics may depend sensitively on the initial configuration. This paper presents an initial survey of the onset of inflation in multifield models, via qualitative lattice-based simulations that do not include local gravitational backreaction. Using hybrid inflation as a test model, our results suggest that small subhorizon inhomogeneities do play a key role in determining whether inflation begins in multifield scenarios. Interestingly, some configurations which do not inflate in the homogeneous limit "succeed" after inhomogeneity is included, while other initial configurations which inflate in the homogeneous limit "fail" when inhomogeneity is added.

preprint2014arXiv

Large Scale Anomalies in the Microwave Background: Causation and Correlation

Most treatments of large scale anomalies in the microwave sky are a posteriori, with unquantified look-elsewhere effects. We contrast these with physical models of specific inhomogeneities in the early universe which then generate apparent anomalies. Physical models predict correlations between candidate anomalies, as well as the corresponding signals in polarization and large scale structure, reducing the impact of cosmic variance. We compute the apparent spatial curvature associated with large-scale inhomogeneities and show that it is typically small, allowing for a self-consistent analysis. Focussing on a single large plane wave inhomogeneity, we show this can contribute to low-l mode alignment and odd-even asymmetry in the power spectra and the best fit inhomogeneity accounts for a significant part of the claimed odd-even asymmetry. We argue that this approach can be generalized to provide a more quantitative assessment of potential large scale anomalies in the universe.

preprint2014arXiv

MultiModeCode: An efficient numerical solver for multifield inflation

We present MultiModeCode, a Fortran 95/2000 package for the numerical exploration of multifield inflation models. This program facilitates efficient Monte Carlo sampling of prior probabilities for inflationary model parameters and initial conditions and is the first publicly available code that can efficiently generate large sample-sets for inflation models with $\mathcal O(100)$ fields. The code numerically solves the equations of motion for the background and first-order perturbations of multi-field inflation models with canonical kinetic terms and arbitrary potentials, providing the adiabatic, isocurvature, and tensor power spectra at the end of inflation. For models with sum-separable potentials MultiModeCode also computes the slow-roll prediction via the $δN$ formalism for easy model exploration and validation. We pay particular attention to the isocurvature perturbations as the system approaches the adiabatic limit, showing how to avoid numerical instabilities that affect some other approaches to this problem. We demonstrate the use of MultiModeCode by exploring a few toy models. Finally, we give a concise review of multifield perturbation theory and a user's manual for the program.

preprint2014arXiv

Simple predictions from multifield inflationary models

We explore whether multifield inflationary models make unambiguous predictions for fundamental cosmological observables. Focusing on $N$-quadratic inflation, we numerically evaluate the full perturbation equations for models with 2, 3, and $\mathcal{O}(100)$ fields, using several distinct methods for specifying the initial values of the background fields. All scenarios are highly predictive, with the probability distribution functions of the cosmological observables becoming more sharply peaked as $N$ increases. For $N=100$ fields, 95% of our Monte Carlo samples fall in the ranges $n_s \in (0.9455,0.9534)$; $α\in (-9.741,-7.047)\times 10^{-4}$; $r\in(0.1445,0.1449)$; and $r_\mathrm{iso} \in (0.02137,3.510)\times 10^{-3}$ for the spectral index, running, tensor-to-scalar ratio, and isocurvature-to-adiabatic ratio, respectively. The expected amplitude of isocurvature perturbations grows with $N$, raising the possibility that many-field models may be sensitive to post-inflationary physics and suggesting new avenues for testing these scenarios.

preprint2014arXiv

The Knotted Sky I: Planck constraints on the primordial power spectrum

Using the temperature data from Planck we search for departures from a power-law primordial power spectrum, employing Bayesian model-selection and posterior probabilities. We parametrize the spectrum with $n$ knots located at arbitrary values of $\log{k}$, with both linear and cubic splines. This formulation recovers both slow modulations and sharp transitions in the primordial spectrum. The power spectrum is well-fit by a featureless, power-law at wavenumbers $k>10^{-3} \, \mathrm{Mpc}^{-1}$. A modulated primordial spectrum yields a better fit relative to $Λ$CDM at large scales, but there is no strong evidence for a departure from a power-law spectrum. Moreover, using simulated maps we show that a local feature at $k \sim 10^{-3} \, \mathrm{Mpc}^{-1}$ can mimic the suppression of large-scale power. With multi-knot spectra we see only small changes in the posterior distributions for the other free parameters in the standard $Λ$CDM universe. Lastly, we investigate whether the hemispherical power asymmetry is explained by independent features in the primordial power spectrum in each ecliptic hemisphere, but find no significant differences between them.

preprint2014arXiv

The Knotted Sky II: Does BICEP2 require a nontrivial primordial power spectrum?

An inflationary gravitational wave background consistent with BICEP2 is difficult to reconcile with a simple power-law spectrum of primordial scalar perturbations. Tensor modes contribute to the temperature anisotropies at multipoles with $l\lesssim 100$, and this effect --- together with a prior on the form of the scalar perturbations --- was the source of previous bounds on the tensor-to-scalar ratio. We compute Bayesian evidence for combined fits to BICEP2 and Planck for three nontrivial primordial spectra: a) a running spectral index, b) a cutoff at fixed wavenumber, and c) a spectrum described by a linear spline with a single internal knot. We find no evidence for a cutoff, weak evidence for a running index, and significant evidence for a "broken" spectrum. Taken at face-value, the BICEP2 results require two new inflationary parameters in order to describe both the broken scale invariance in the perturbation spectrum and the observed tensor-to-scalar ratio. Alternatively, this tension may be resolved by additional data and more detailed analyses.

preprint2013arXiv

Constraining Monodromy Inflation

We use cosmic microwave background (CMB) data from the 9-year WMAP release to derive constraints on monodromy inflation, which is characterized by a linear inflaton potential with a periodic modulation. We identify two possible periodic modulations that significantly improve the fit, but it is unclear whether this improvement is associated with a "signal" or with the scatter in the measured angular power spectrum. The unmodulated potential is mildly favored by standard model selection criteria. A modulated inflationary potential can generate substantial primordial non-Gaussianity, of a specific and characteristic form. For the best-fit parameters to the WMAP angular power spectrum, the corresponding non-Gaussianity might be detectable in upcoming CMB data, allowing nontrivial consistency checks on the predictions of a modulated inflationary potential.

preprint2013arXiv

Gravitational Waves from Oscillon Preheating

Oscillons are long-lived, localized excitations of nonlinear scalar fields which may be copiously produced during preheating after inflation, leading to a possible oscillon-dominated phase in the early Universe. For example, this can happen after axion monodromy inflation, on which we run our simulations. We investigate the stochastic gravitational wave background associated with an oscillon-dominated phase. An isolated oscillon is spherically symmetric and does not radiate gravitational waves, and we show that the flux of gravitational radiation generated between oscillons is also small. However, a significant stochastic gravitational wave background may be generated during preheating itself (i.e, when oscillons are forming), and in this case the characteristic size of the oscillons is imprinted on the gravitational wave power spectrum, which has multiple, distinct peaks.

preprint2013arXiv

Initial conditions and sampling for multifield inflation

We investigate the initial conditions problem for multifield inflation. In these scenarios the pre-inflationary dynamics can be chaotic, increasing the sensitivity of the onset of inflation to the initial data even in the homogeneous limit. To analyze physically equivalent scenarios we compare initial conditions at fixed energy. This ensures that each trajectory is counted once and only once, since the energy density decreases monotonically. We present a full analysis of hybrid inflation that reveals a greater degree of long range order in the set of "successful" initial conditions than was previously apparent. We explore the effective smoothing scale for the fractal set of successful initial conditions induced by the finite duration of the pre-inflationary phase. The role of the prior information used to specify the initial data is discussed in terms of Bayesian sampling.

preprint2013arXiv

Planck Constraints on Monodromy Inflation

We use data from the nominal Planck mission to constrain modulations in the primordial power spectrum associated with monodromy inflation. The largest improvement in fit relative to the unmodulated model has Δχ^2~10 and we find no evidence for a primordial signal, in contrast to a previous analysis of the WMAP9 dataset, for which Δχ^2~20. The Planck and WMAP9 results are broadly consistent on angular scales where they are expected to agree as far as best-fit values are concerned. However, even on these scales the significance of the signal is reduced in Planck relative to WMAP, and is consistent with a fit to the noise associated with cosmic variance. Our results motivate both a detailed comparison between the two experiments and a more careful study of the theoretical predictions of monodromy inflation.

preprint2013arXiv

Supersymmetry, Nonthermal Dark Matter and Precision Cosmology

Within the Minimal Supersymmetric Standard Model (MSSM), LHC bounds suggest that scalar superpartner masses are far above the electroweak scale. Given a high superpartner mass, nonthermal dark matter is a viable alternative to WIMP dark matter generated via freezeout. In the presence of moduli fields nonthermal dark matter production is associated with a long matter dominated phase, modifying the spectral index and primordial tensor amplitude relative to those in a thermalized primordial universe. Nonthermal dark matter can have a higher self-interaction cross-section than its thermal counterpart, enhancing astrophysical bounds on its annihilation signals. We constrain the contributions to the neutralino mass from the bino, wino and higgsino using existing astrophysical bounds and direct detection experiments for models with nonthermal neutralino dark matter. Using these constraints we quantify the expected change to inflationary observables resulting from the nonthermal phase.

preprint2012arXiv

Bayesian Analysis of Inflation II: Model Selection and Constraints on Reheating

We discuss the model selection problem for inflationary cosmology. We couple ModeCode, a publicly-available numerical solver for the primordial perturbation spectra, to the nested sampler MultiNest, in order to efficiently compute Bayesian evidence. Particular attention is paid to the specification of physically realistic priors, including the parametrization of the post-inflationary expansion and associated thermalization scale. It is confirmed that while present-day data tightly constrains the properties of the power spectrum, it cannot usefully distinguish between the members of a large class of simple inflationary models. We also compute evidence using a simulated Planck likelihood, showing that while Planck will have more power than WMAP to discriminate between inflationary models, it will not definitively address the inflationary model selection problem on its own. However, Planck will place very tight constraints on any model with more than one observationally-distinct inflationary regime -- e.g. the large- and small-field limits of the hilltop inflation model -- and put useful limits on different reheating scenarios for a given model.

preprint2012arXiv

Bayesian Analysis of Inflation III: Slow Roll Reconstruction Using Model Selection

We implement Slow Roll Reconstruction -- an optimal solution to the inverse problem for inflationary cosmology -- within ModeCode, a publicly available solver for the inflationary dynamics. We obtain up-to-date constraints on the reconstructed inflationary potential, derived from the WMAP 7-year dataset and South Pole Telescope observations, combined with large scale structure data derived from SDSS Data Release 7. Using ModeCode in conjunction with the MultiNest sampler, we compute Bayesian evidence for the reconstructed potential at each order in the truncated slow roll hierarchy. We find that the data are well-described by the first two slow roll parameters, εand η, and that there is no need to include a nontrivial ξparameter.

preprint2011arXiv

Bayesian Analysis of Inflation: Parameter Estimation for Single Field Models

Future astrophysical datasets promise to strengthen constraints on models of inflation, and extracting these constraints requires methods and tools commensurate with the quality of the data. In this paper we describe ModeCode, a new, publicly available code that computes the primordial scalar and tensor power spectra for single field inflationary models. ModeCode solves the inflationary mode equations numerically, avoiding the slow roll approximation. It is interfaced with CAMB and CosmoMC to compute cosmic microwave background angular power spectra and perform likelihood analysis and parameter estimation. ModeCode is easily extendable to additional models of inflation, and future updates will include Bayesian model comparison. Errors from ModeCode contribute negligibly to the error budget for analyses of data from Planck or other next generation experiments. We constrain representative single field models (phi^n with n=2/3, 1, 2, and 4, natural inflation, and "hilltop" inflation) using current data, and provide forecasts for Planck. From current data, we obtain weak but nontrivial limits on the post-inflationary physics, which is a significant source of uncertainty in the predictions of inflationary models, while we find that Planck will dramatically improve these constraints. In particular, Planck will link the inflationary dynamics with the post-inflationary growth of the horizon, and thus begin to probe the "primordial dark ages" between TeV and GUT scale energies.

preprint2011arXiv

Caching and Interpolated Likelihoods: Accelerating Cosmological Monte Carlo Markov Chains

We describe a novel approach to accelerating Monte Carlo Markov Chains. Our focus is cosmological parameter estimation, but the algorithm is applicable to any problem for which the likelihood surface is a smooth function of the free parameters and computationally expensive to evaluate. We generate a high-order interpolating polynomial for the log-likelihood using the first points gathered by the Markov chains as a training set. This polynomial then accurately computes the majority of the likelihoods needed in the latter parts of the chains. We implement a simple version of this algorithm as a patch (InterpMC) to CosmoMC and show that it accelerates parameter estimatation by a factor of between two and four for well-converged chains. The current code is primarily intended as a "proof of concept", and we argue that there is considerable room for further performance gains. Unlike other approaches to accelerating parameter fits, we make no use of precomputed training sets or special choices of variables, and InterpMC is almost entirely transparent to the user.

preprint2011arXiv

Constraining holographic inflation with WMAP

In a class of recently proposed models, the early universe is strongly coupled and described holographically by a three-dimensional, weakly coupled, super-renormalizable quantum field theory. This scenario leads to a power spectrum of scalar perturbations that differs from the usual empirical LCDM form and the predictions of generic models of single field, slow roll inflation. This spectrum is characterized by two parameters: an amplitude, and a parameter g related to the coupling constant of the dual theory. We estimate these parameters, using WMAP and other astrophysical data. We compute Bayesian evidence for both the holographic model and standard LCDM and find that their difference is not significant, although LCDM provides a somewhat better fit to the data. However, it appears that Planck will permit a definitive test of this holographic scenario.

preprint2011arXiv

Inflation and the Scale Dependent Spectral Index: Prospects and Strategies

We consider the running of the spectral index as a probe of both inflation itself, and of the overall evolution of the very early universe. Surveying a collection of simple single field inflationary models, we confirm that the magnitude of the running is relatively consistent, unlike the tensor amplitude, which varies by orders of magnitude. Given this target, we confirm that the running is potentially detectable by future large scale structure or 21 cm observations, but that only the most futuristic measurements can distinguish between these models on the basis of their running. For any specified inflationary scenario, the combination of the running index and unknown post-inflationary expansion history induces a theoretical uncertainty in the predicted value of the spectral index. This effect can easily dominate the statistical uncertainty with which Planck and its successors are expected to measure the spectral index. More positively, upcoming cosmological experiments thus provide an intriguing probe of physics between TeV and GUT scales by constraining the reheating history associated with any specified inflationary model, opening a window into the "primordial dark age" that follows the end of inflation.

preprint2011arXiv

Oscillons After Inflation

Oscillons are massive, long-lived, localized excitations of a scalar field. We show that in a large class of well-motivated single-field models, inflation is followed by self-resonance, leading to copious oscillon generation and a lengthy period of oscillon domination. These models are characterized by an inflaton potential which has a quadratic minimum and is shallower than quadratic away from the minimum. This set includes both string monodromy models and a class of supergravity inspired scenarios, and is in good agreement with the current central values of the concordance cosmology parameters. We assume that the inflaton is weakly coupled to other fields, so as not to quickly drain energy from the oscillons or prevent them from forming. An oscillon-dominated universe has a greatly enhanced primordial power spectrum on very small scales relative to that seen with a quadratic potential, possibly leading to novel gravitational effects in the early universe.

preprint2010arXiv

Delayed Reheating and the Breakdown of Coherent Oscillations

We analyze the evolution of the perturbations in the inflaton field and metric following the end of inflation. We present accurate analytic approximations for the perturbations, showing that the coherent oscillations of the post-inflationary condensate necessarily break down long before any current phenomenological constraints require the universe to become radiation dominated. Further, the breakdown occurs on length-scales equivalent to the comoving post-inflationary horizon size. This work has implications for both the inflationary "matching" problem, and the possible generation of a stochastic gravitational wave background in the post-inflationary universe.

preprint2010arXiv

Inflaton Fragmentation and Oscillon Formation in Three Dimensions

Analytical arguments suggest that a large class of scalar field potentials permit the existence of oscillons -- pseudo-stable, non-topological solitons -- in three spatial dimensions. In this paper we numerically explore oscillon solutions in three dimensions. We confirm the existence of these field configurations as solutions to the Klein-Gorden equation in an expanding background, and verify the predictions of Amin and Shirokoff for the characteristics of individual oscillons for their model. Further, we demonstrate that significant numbers of oscillons can be generated via fragmentation of the inflaton condensate, consistent with the analysis of Amin. These emergent oscillons can easily dominate the post-inflationary universe. Finally, both analytic and numerical results suggest that oscillons are stable on timescales longer than the post-inflationary Hubble time. Consequently, the post-inflationary universe can contain an effective matter-dominated phase, during which it is dominated by localized concentrations of scalar field matter.

preprint2010arXiv

PSpectRe: A Pseudo-Spectral Code for (P)reheating

PSpectRe is a C++ program that uses Fourier-space pseudo-spectral methods to evolve interacting scalar fields in an expanding universe. PSpectRe is optimized for the analysis of parametric resonance in the post-inflationary universe, and provides an alternative to finite differencing codes, such as Defrost and LatticeEasy. PSpectRe has both second- (Velocity-Verlet) and fourth-order (Runge-Kutta) time integrators. Given the same number of spatial points and/or momentum modes, PSpectRe is not significantly slower than finite differencing codes, despite the need for multiple Fourier transforms at each timestep, and exhibits excellent energy conservation. Further, by computing the post-resonance equation of state, we show that in some circumstances PSpectRe obtains reliable results while using substantially fewer points than a finite differencing code. PSpectRe is designed to be easily extended to other problems in early-universe cosmology, including the generation of gravitational waves during phase transitions and pre-inflationary bubble collisions. Specific applications of this code will be pursued in future work.

preprint2008arXiv

Primordial Black Holes, Eternal Inflation, and the Inflationary Parameter Space after WMAP5

We consider constraints on inflation driven by a single, minimally coupled scalar field in the light of the WMAP5 dataset, as well as ACBAR and the SuperNova Legacy Survey. We use the Slow Roll Reconstruction algorithm to derive optimal constraints on the inflationary parameter space. The scale dependence in the slope of the scalar spectrum permitted by WMAP5 is large enough to lead to viable models where the small scale perturbations have a substantial amplitude when extrapolated to the end of inflation. We find that excluding parameter values which would cause the overproduction of primordial black holes or even the onset of eternal inflation leads to potentially significant constraints on the slow roll parameters. Finally, we present a more sophisticated approach to including priors based on the total duration of inflation, and discuss the resulting restrictions on the inflationary parameter space.

preprint2007arXiv

Large Non-Gaussianities in Single Field Inflation

We compute the 3-point correlation function for a general model of inflation driven by a single, minimally coupled scalar field. Our approach is based on the numerical evaluation of both the perturbation equations and the integrals which contribute to the 3-point function. Consequently, we can analyze models where the potential has a "feature", in the vicinity of which the slow roll parameters may take on large, transient values. This introduces both scale and shape dependent non-Gaussianities into the primordial perturbations. As an example of our methodology, we examine the ``step'' potentials which have been invoked to improve the fit to the glitch in the $<TT>$ $C_l$ for $l \sim 30$, present in both the one and three year WMAP data sets. We show that for the typical parameter values, the non-Gaussianities associated with the step are far larger than those in standard slow roll inflation, and may even be within reach of a next generation CMB experiment such as Planck. More generally, we use this example to explain that while adding features to potential can improve the fit to the 2-point function, these are generically associated with a greatly enhanced signal at the 3-point level. Moreover, this 3-point signal will have a very nontrivial shape and scale dependence, which is correlated with the form of the 2-point function, and may thus lead to a consistency check on the models of inflation with non-smooth potentials.

preprint2007arXiv

Non-Gaussianities in Multi-field Inflation

We compute the amplitude of the non-Gaussianities in inflationary models with multiple, uncoupled scalar fields. This calculation thus applies to all models of assisted inflation, including N-flation, where inflation is driven by multiple axion fields arising from shift symmetries in a flux stabilized string vacuum. The non-Gaussianities are associated with nonlinear evolution of the field (and density) perturbations, characterized by the parameter $f_{NL}$. We derive a general expression for the nonlinear parameter, incorporating the evolution of perturbations after horizon-crossing. This is valid for arbitrary separable potentials during slow roll. To develop an intuitive understanding of this system and to demonstrate the applicability of the formalism we examine several cases with quadratic potentials: two-field models with a wide range of mass ratios, and a general N-field model with a narrow mass spectrum. We uncover that $f_{NL}$ is suppressed as the number of e-foldings grows, and that this suppression is increased in models with a broad spectrum of masses. On the other hand, we find no enhancement to $f_{NL}$ that increases with the number of fields. We thus conclude that the production of a large non-Gaussian signal in multi-field models of inflation is very unlikely as long as fields are slowly rolling and potentials are of simple, quadratic form. Finally, we compute a spectrum for the scalar spectral index that incorporates the nonlinear corrections to the fields' evolution.

preprint2006arXiv

Counting Pockets with World Lines in Eternal Inflation

We consider the long standing puzzle of how to obtain meaningful probabilities in eternal inflation. We demonstrate a new algorithm to compute the probability distribution of pocket universe types, given a multivacua inflationary potential. The computed probability distribution is finite and manifestly gauge-independent. We argue that in some scenarios this technique can be applied to disfavor some eternally inflating potentials.

preprint2006arXiv

Implications of a Running Spectral Index for Slow Roll Inflation

We analyze the weak (2 sigma) evidence for a running spectral index seen in the three-year WMAP dataset and its implications for single field, slow roll inflation. We assume that the running is comparable to the central value found from the WMAP data analysis, and use the Hubble Slow Roll formalism to follow the evolution of the slow roll parameters. For all parameter choices consistent with a large, negative running, single field, slow roll inflation lasts less than 30 efolds after CMB scales leave the horizon. Thus, a definitive observation of a large negative running would imply that any inflationary phase requires multiple fields or the breakdown of slow roll. Alternatively, if single field, slow roll inflation is sources the primordial fluctuations, we can expect the observed running to move much closer to zero as the CMB is measured more accurately at small angular scales.

preprint2006arXiv

Random Matrices and the Spectrum of N-flation

N-flation is a promising embedding of inflation in string theory in which many string axions combine to drive inflation. We characterize the dynamics of a general N-flation model with non-degenerate axion masses. Although the precise mass of a single axion depends on compactification details in a complicated way, the distribution of masses can be computed with very limited knowledge of microscopics: the shape of the mass distribution is an emergent property. We use random matrix theory to show that a typical N-flation model has a spectrum of masses distributed according to the Marchenko-Pastur law. This distribution depends on a single parameter, the number of axions divided by the dimension of the moduli space. We use this result to describe the inflationary dynamics and phenomenology of a general N-flation model. We produce an ensemble of models and use numerical integration to track the axions' evolution and the resulting scalar power spectrum. For realistic initial conditions, the power spectrum is considerably more red than in single-field $m^2ϕ^2$ inflation. We conclude that random matrix models of N-flation are surprisingly tractable and have a rich phenomenology that differs in testable ways from that of single-field $m^2ϕ^2$ inflation.

preprint2006arXiv

Recovering the Inflationary Potential and Primordial Power Spectrum With a Slow Roll Prior: Methodology and Application to WMAP 3 Year Data

We introduce a new method for applying an inflationary prior to a cosmological dataset that includes relations between observables at arbitrary order in the slow roll expansion. The process is based on the inflationary flow equations, and the slow roll parameters appear explicitly in the cosmological parameter set. We contrast our method to other ways of imposing an inflationary prior on a cosmological dataset, and argue that this method is ideal for use with heterogeneous datasets. In particular, it would be well suited to exploiting any direct detection of fundamental tensor modes by a BBO-style mission. To demonstrate the practical use of this method we apply it to the WMAPI+All dataset, and the newly released WMAPII dataset on its own and together with the SDSS data. We find that all basic classes of single field inflationary models are still allowed at the 1-2sigma level, but the overall parameter space is sharply constrained. In particular, we find evidence that the combination of WMAPII+SDSS is sensitive to effects arising from terms that are quadratic in the two leading-order slow roll parameters.

preprint2006arXiv

Slow Roll Reconstruction: Constraints on Inflation from the 3 Year WMAP Dataset

We study the constraints on the inflationary parameter space derived from the 3 year WMAP dataset using ``slow roll reconstruction'', using the SDSS galaxy power spectrum to gain further leverage where appropriate. This approach inserts the inflationary slow roll parameters directly into a Monte Carlo Markov chain estimate of the cosmological parameters, and uses the inflationary flow hierarchy to compute the parameters' scale-dependence. We work with the first three parameters (epsilon, eta and xi) and pay close attention to the possibility that the 3 year WMAP dataset contains evidence for a ``running'' spectral index, which is dominated by the xi term. Mirroring the WMAP team's analysis we find that the permitted distribution of xi is broad, and centered away from zero. However, when we require that inflationary parameters yield at least 30 additional e-folds of inflation after the largest observable scales leave the horizon, the bounds on xi tighten dramatically. We make use of the absence of an explicit pivot scale in the slow roll reconstruction formalism to determine the dependence of the computed parameter distributions on the pivot. We show that the choice of pivot has a significant effect on the inferred constraints on the inflationary variables, and the spectral index and running derived from them. Finally, we argue that the next round of cosmological data can be expected to place very stringent constraints on the region of parameter space open to single field models of slow roll inflation.

preprint2005arXiv

Boundary Effective Field Theory and Trans-Planckian Perturbations: Astrophysical Implications

We contrast two approaches to calculating trans-Planckian corrections to the inflationary perturbation spectrum: the New Physics Hypersurface [NPH] model, in which modes are normalized when their physical wavelength first exceeds a critical value, and the Boundary Effective Field Theory [BEFT] approach, where the initial conditions for all modes are set at the same time, and modified by higher dimensional operators enumerated via an effective field theory calculation. We show that these two approaches -- as currently implemented -- lead to radically different expectations for the trans-Planckian corrections to the CMB and emphasize that in the BEFT formalism we expect the perturbation spectrum to be dominated by quantum gravity corrections for all scales shorter than some critical value. Conversely, in the NPH case the quantum effects only dominate the longest modes that are typically much larger than the present horizon size. Furthermore, the onset of the breakdown in the standard inflationary perturbation calculation predicted by the BEFT formalism is likely to be associated with a feature in the perturbation spectrum, and we discuss the observational signatures of this feature in both CMB and large scale structure observations. Finally, we discuss possible modifications to both calculational frameworks that would resolve the contradictions identified here.

preprint2005arXiv

Observing Trans-Planckian Signatures in the Cosmic Microwave Background

We examine the constraints cosmological observations can place on any trans-Planckian corrections to the primordial spectrum of perturbations underlying the anisotropies in the Cosmic Microwave Background. We focus on models of trans-Planckian physics which lead to a modulated primordial spectrum. Rather than looking at a generic modulated spectrum, our calculations are based on a specific model, and are intended as a case study for the sort of constraints one could hope to apply on a well-motivated model of trans-Planckian physics. We present results for three different approaches -- a grid search in a subset of the overall parameter space, a Fisher matrix estimate of the likely error ellipses, and a Monte Carlo Markov Chain fit to a simulated CMB sky. As was seen in previous analyses, the likelihood space has multiple peaks, and we show that their distribution can be reproduced via a simple semi-analytic argument. All three methods lead to broadly similar results. We vary 10 cosmological parameters (including two related to the trans-Planckian terms) and show that the amplitude of the tensor perturbations is directly correlated with the detectability of any trans-Planckian modulation. We argue that this is likely to be true for any trans-Planckian modulation in the paradigm of slow-roll inflation. For the specific case we consider, we conclude that if the tensor to scalar ratio, $r \sim 0.15$, the ratio between the inflationary Hubble scale $H$, and the scale of new physics $M$ has to be on the order of 0.004 if the modulation is detectable at the 2$σ$ level. For a lower value of $r$, the bound on $H/M$ becomes looser.

preprint2005arXiv

The Hubble Slow Roll Expansion for Multi Field Inflation

We examine the dynamics of inflation driven by multiple, interacting scalar fields and derive a multi field version of the Hubble slow roll expansion. We show that the properties of this expansion naturally generalize those of the single field case. We present an analogous hierarchy of slow roll parameters, and derive the system of ``flow'' equations that describes their evolution, and show that when this system is truncated at finite order, it can always be solved exactly. Lastly, we express the scalar and tensor perturbation spectra in terms of the slow roll parameters.

preprint2004arXiv

String windings in the early universe

We study string dynamics in the early universe. Our motivation is the proposal of Brandenberger and Vafa, that string winding modes may play a key role in decompactifying three spatial dimensions. We model the universe as a homogeneous but anisotropic 9-torus filled with a gas of excited strings. We adopt initial conditions which fix the dilaton and the volume of the torus, but otherwise assume all states are equally likely. We study the evolution of the system both analytically and numerically to determine the late-time behavior. We find that, although dynamical evolution can indeed lead to three large spatial dimensions, such an outcome is not statistically favored.

preprint2002arXiv

Cosmological String Gas on Orbifolds

It has long been known that strings wound around incontractible cycles can play a vital role in cosmology. In particular, in a spacetime with toroidal spatial hypersurfaces, the dynamics of the winding modes may help yield three large spatial dimensions. However, toroidal compactifications are phenomenologically unrealistic. In this paper we therefore take a first step toward extending these cosmological considerations to $D$-dimensional toroidal orbifolds. We use numerical simulation to study the timescales over which "pseudo-wound" strings unwind on these orbifolds with trivial fundamental group. We show that pseudo-wound strings can persist for many ``Hubble times'' in some of these spaces, suggesting that they may affect the dynamics in the same way as genuinely wound strings. We also outline some possible extensions that include higher-dimensional wrapped branes.

preprint2001arXiv

Inflation as a Probe of Short Distance Physics

We show that a string-inspired Planck scale modification of general relativity can have observable cosmological effects. Specifically, we present a complete analysis of the inflationary perturbation spectrum produced by a phenomenological Lagrangian that has a standard form on large scales but incorporates a string-inspired short distance cutoff, and find a deviation from the standard result. We use the de Sitter calculation as the basis of a qualitative analysis of other inflationary backgrounds, arguing that in these cases the cutoff could have a more pronounced effect, changing the shape of the spectrum. Moreover, the computational approach developed here can be used to provide unambiguous calculations of the perturbation spectrum in other heuristic models that modify trans-Planckian physics and thereby determine their impact on the inflationary perturbation spectrum. Finally, we argue that this model may provide an exception to constraints, recently proposed by Tanaka and Starobinsky, on the ability of Planck-scale physics to modify the cosmological spectrum.

preprint1999arXiv

Preheating and the Einstein Field Equations

We inaugurate a framework for studying preheating and parametric resonance after inflation without resorting to any approximations, either in gravitational perturbation theory or in the evolution of the field(s). We do this by numerically solving the Einstein field equations in the post-inflationary universe. In this letter we show how to compare our results to those of gauge invariant perturbation theory. We then verify Finelli and Brandenberger's analysis (hep-ph/9809490) of super-horizon modes in $m^2ϕ^2$ inflation, showing that they are not amplified by resonant effects. Lastly, we make a preliminary survey of the nonlinear couplings between modes, which will be important in models where the primordial metric perturbations undergo parametric amplification.

preprint1993arXiv

Exact Superstring Motivated Cosmological Models

We present a number of new, exact scalar field cosmologies where the potential consists of two or more exponential terms. Such potentials are motivated by supergravity or superstring models formulated in higher dimensional spacetimes that have been compactified to (3+1)-dimensions. We have found solutions in both curved and flat Robertson Walker spacetimes. These models have a diverse range of properties and often possess several distinct phases, with a smooth transition between ordinary and inflationary expansion. While exponential potentials typically produce powerlaw inflation, we find models where the inflationary period contains eras of both powerlaw and exponential growth.