-
Robust posterior inference when statistically emulating forward simulations
Authors:
Grigor Aslanyan,
Richard Easther,
Nathan Musoke,
Layne C. Price
Abstract:
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 cosmolo…
▽ More
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.
△ Less
Submitted 24 April, 2020;
originally announced April 2020.
-
Estimating Cosmological Parameters from the Dark Matter Distribution
Authors:
Siamak Ravanbakhsh,
Junier Oliva,
Sebastien Fromenteau,
Layne C. Price,
Shirley Ho,
Jeff Schneider,
Barnabas Poczos
Abstract:
A grand challenge of the 21st century cosmology is to accurately estimate the cosmological parameters of our Universe. A major approach to estimating the cosmological parameters is to use the large-scale matter distribution of the Universe. Galaxy surveys provide the means to map out cosmic large-scale structure in three dimensions. Information about galaxy locations is typically summarized in a "…
▽ More
A grand challenge of the 21st century cosmology is to accurately estimate the cosmological parameters of our Universe. A major approach to estimating the cosmological parameters is to use the large-scale matter distribution of the Universe. Galaxy surveys provide the means to map out cosmic large-scale structure in three dimensions. Information about galaxy locations is typically summarized in a "single" function of scale, such as the galaxy correlation function or power-spectrum. We show that it is possible to estimate these cosmological parameters directly from the distribution of matter. This paper presents the application of deep 3D convolutional networks to volumetric representation of dark-matter simulations as well as the results obtained using a recently proposed distribution regression framework, showing that machine learning techniques are comparable to, and can sometimes outperform, maximum-likelihood point estimates using "cosmological models". This opens the way to estimating the parameters of our Universe with higher accuracy.
△ Less
Submitted 6 November, 2017;
originally announced November 2017.
-
Effect of reheating on predictions following multiple-field inflation
Authors:
Selim C. Hotinli,
Jonathan Frazer,
Andrew H. Jaffe,
Joel Meyers,
Layne C. Price,
Ewan R. M. Tarrant
Abstract:
We study the sensitivity of cosmological observables to the reheating phase following inflation driven by many scalar fields. We describe a method which allows semi-analytic treatment of the impact of perturbative reheating on cosmological perturbations using the sudden decay approximation. Focusing on $\mathcal{N}$-quadratic inflation, we show how the scalar spectral index and tensor-to-scalar ra…
▽ More
We study the sensitivity of cosmological observables to the reheating phase following inflation driven by many scalar fields. We describe a method which allows semi-analytic treatment of the impact of perturbative reheating on cosmological perturbations using the sudden decay approximation. Focusing on $\mathcal{N}$-quadratic inflation, we show how the scalar spectral index and tensor-to-scalar ratio are affected by the rates at which the scalar fields decay into radiation. We find that for certain choices of decay rates, reheating following multiple-field inflation can have a significant impact on the prediction of cosmological observables.
△ Less
Submitted 23 January, 2018; v1 submitted 24 October, 2017;
originally announced October 2017.
-
The Spectrum of the Axion Dark Sector
Authors:
Matthew J. Stott,
David J. E. Marsh,
Chakrit Pongkitivanichkul,
Layne C. Price,
Bobby S. Acharya
Abstract:
Axions arise in many theoretical extensions of the Standard Model of particle physics, in particular the "string axiverse". If the axion masses, $m_a$, and (effective) decay constants, $f_a$, lie in specific ranges, then axions contribute to the cosmological dark matter and dark energy densities. We compute the background cosmological (quasi-)observables for models with a large number of axion fie…
▽ More
Axions arise in many theoretical extensions of the Standard Model of particle physics, in particular the "string axiverse". If the axion masses, $m_a$, and (effective) decay constants, $f_a$, lie in specific ranges, then axions contribute to the cosmological dark matter and dark energy densities. We compute the background cosmological (quasi-)observables for models with a large number of axion fields, $n_{\rm ax}\sim \mathcal{O}(10-100)$, with the masses and decay constants drawn from statistical distributions. This reduces the number of parameters from $2n_{\rm ax}$ to a small number of "hyperparameters". We consider a number of distributions, from those motivated purely by statistical considerations, to those where the structure is specified according to a class of M-theory models. Using Bayesian methods we are able to constrain the hyperparameters of the distributions. In some cases the hyperparameters can be related to string theory, e.g. constraining the number ratio of axions to moduli, or the typical decay constant scale needed to provide the correct relic densities. Our methodology incorporates the use of both random matrix theory and Bayesian networks.
△ Less
Submitted 10 June, 2017;
originally announced June 2017.
-
Testing for New Physics: Neutrinos and the Primordial Power Spectrum
Authors:
Nicolas Canac,
Grigor Aslanyan,
Kevork N. Abazajian,
Richard Easther,
Layne C. Price
Abstract:
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 degre…
▽ More
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.
△ Less
Submitted 20 September, 2016; v1 submitted 9 June, 2016;
originally announced June 2016.
-
Reconstructing the redshift evolution of escaped ionizing flux from early galaxies with Planck and HST observations
Authors:
Layne C. Price,
Hy Trac,
Renyue Cen
Abstract:
While galaxies at $6 \lesssim z \lesssim 10$ are believed to dominate the epoch of cosmic reionization, the escape fraction of ionizing flux $f_\mathrm{esc}$ and the photon production rate $\dot n_γ$ from these galaxies must vary with redshift to simultaneously match CMB and low-redshift observations. We constrain $f_\mathrm{esc}(z)$ and $\dot n_γ(z)$ with Planck 2015 measurements of the Thomson o…
▽ More
While galaxies at $6 \lesssim z \lesssim 10$ are believed to dominate the epoch of cosmic reionization, the escape fraction of ionizing flux $f_\mathrm{esc}$ and the photon production rate $\dot n_γ$ from these galaxies must vary with redshift to simultaneously match CMB and low-redshift observations. We constrain $f_\mathrm{esc}(z)$ and $\dot n_γ(z)$ with Planck 2015 measurements of the Thomson optical depth $τ$, recent low multipole E-mode polarization measurements from Planck 2016, SDSS BAO data, and $3 \lesssim z \lesssim 10$ galaxy observations. We compare different galaxy luminosity functions that are calibrated to HST observations, using both parametric and non-parametric statistical methods that marginalize over the effective clumping factor $C_\mathrm{HII}$, the LyC production efficiency $ξ_\mathrm{ion}$, and the time-evolution of the UV limiting magnitude $dM_\mathrm{SF}/dz$. Using a power-law model, we find $f_\mathrm{esc} \lesssim 0.5$ at $z=8$ with slope $β\gtrsim 2.0$ at $68\%$ confidence with little dependence on the galaxy luminosity function or data, although there is non-negligible probability for no redshift evolution $β\sim 0$ or small escape fraction $f_\mathrm{esc} \sim 10^{-2}$. A non-parametric form for $f_\mathrm{esc}(z)$ evolves significantly with redshift, yielding $f_\mathrm{esc} \sim 0.2, 0.3, 0.6$ at $z=6,9,12$, respectively. However, a model-independent reconstruction of $\dot n_γ(z)$ predicts a suppressed escaped photon production rate at $z=9$ for the latest Planck data compared to the other models, implying a quicker period of reionization. We find evidence for redshift evolution in the limiting magnitude of the galaxy luminosity function for empirical models of the galaxy luminosity function.
△ Less
Submitted 12 May, 2016;
originally announced May 2016.
-
The Velocity Distribution Function of Galaxy Clusters as a Cosmological Probe
Authors:
M. Ntampaka,
H. Trac,
J. Cisewski,
L. C. Price
Abstract:
We present a new approach for quantifying the abundance of galaxy clusters and constraining cosmological parameters using dynamical measurements. In the standard method, galaxy line-of-sight (LOS) velocities, $v$, or velocity dispersions are used to infer cluster masses, $M$, in order to quantify the halo mass function (HMF), $dn(M)/d\log(M)$, which is strongly affected by mass measurement errors.…
▽ More
We present a new approach for quantifying the abundance of galaxy clusters and constraining cosmological parameters using dynamical measurements. In the standard method, galaxy line-of-sight (LOS) velocities, $v$, or velocity dispersions are used to infer cluster masses, $M$, in order to quantify the halo mass function (HMF), $dn(M)/d\log(M)$, which is strongly affected by mass measurement errors. In our new method, the probability distribution of velocities for each cluster in the sample are summed to create a new statistic called the velocity distribution function (VDF), $dn(v)/dv$. The VDF can be measured more directly and precisely than the HMF and it can also be robustly predicted with cosmological simulations which capture the dynamics of subhalos or galaxies. We apply these two methods to mock cluster catalogs and forecast the bias and constraints on the matter density parameter $Ω_m$ and the amplitude of matter fluctuations $σ_8$ in flat $Λ$CDM cosmologies. For an example observation of 200 massive clusters, the VDF with (without) velocity errors constrains the parameter combination $σ_8Ω_m^{0.29\ (0.29)} = 0.587 \pm 0.011\ (0.583 \pm 0.011)$ and shows only minor bias. However, the HMF with dynamical mass errors is biased to low $Ω_m$ and high $σ_8$ and the fiducial model lies well outside of the forecast constraints, prior to accounting for Eddington bias. When the VDF is combined with constraints from the cosmic microwave background (CMB), the degeneracy between cosmological parameters can be significantly reduced. Upcoming spectroscopic surveys that probe larger volumes and fainter magnitudes will provide a larger number of clusters for applying the VDF as a cosmological probe.
△ Less
Submitted 20 October, 2016; v1 submitted 4 February, 2016;
originally announced February 2016.
-
Ultracompact minihalos as probes of inflationary cosmology
Authors:
Grigor Aslanyan,
Layne C. Price,
Jenni Adams,
Torsten Bringmann,
Hamish A. Clark,
Richard Easther,
Geraint F. Lewis,
Pat Scott
Abstract:
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 ga…
▽ More
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.
△ Less
Submitted 20 October, 2016; v1 submitted 14 December, 2015;
originally announced December 2015.
-
Designing and testing inflationary models with Bayesian networks
Authors:
Layne C. Price,
Hiranya V. Peiris,
Jonathan Frazer,
Richard Easther
Abstract:
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,…
▽ More
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.
△ Less
Submitted 29 January, 2016; v1 submitted 30 October, 2015;
originally announced November 2015.
-
Power spectrum oscillations from Planck-suppressed operators in effective field theory motivated monodromy inflation
Authors:
Layne C. Price
Abstract:
We consider a phenomenological model of inflation where the inflaton is the phase of a complex scalar field $Φ$. Planck-suppressed operators of $\mathcal O(f^5/M_\mathrm{pl})$ modify the geometry of the vev $\langle Φ\rangle$ at first order in the decay constant $f$, which adds a first order periodic term to the definition of the canonically normalized inflaton $φ$. This correction to the inflaton…
▽ More
We consider a phenomenological model of inflation where the inflaton is the phase of a complex scalar field $Φ$. Planck-suppressed operators of $\mathcal O(f^5/M_\mathrm{pl})$ modify the geometry of the vev $\langle Φ\rangle$ at first order in the decay constant $f$, which adds a first order periodic term to the definition of the canonically normalized inflaton $φ$. This correction to the inflaton induces a fixed number of extra oscillatory terms in the potential $V \sim θ^p$. We derive the same result in a toy scenario where the vacuum $\langle Φ\rangle$ is an ellipse with an arbitrarily large eccentricity. These extra oscillations change the form of the power spectrum as a function of scale $k$ and provide a possible mechanism for differentiating EFT-motivated inflation from models where the angular shift symmetry is a gauge symmetry.
△ Less
Submitted 18 November, 2015; v1 submitted 29 July, 2015;
originally announced July 2015.
-
Learn-As-You-Go Acceleration of Cosmological Parameter Estimates
Authors:
Grigor Aslanyan,
Richard Easther,
Layne C. Price
Abstract:
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 de…
▽ More
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.
△ Less
Submitted 1 September, 2015; v1 submitted 2 June, 2015;
originally announced June 2015.
-
MultiModeCode: An efficient numerical solver for multifield inflation
Authors:
Layne C. Price,
Jonathan Frazer,
Jiajun Xu,
Hiranya V. Peiris,
Richard Easther
Abstract:
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…
▽ More
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.
△ Less
Submitted 2 October, 2014;
originally announced October 2014.
-
Gravitational wave consistency relations for multifield inflation
Authors:
Layne C. Price,
Hiranya V. Peiris,
Jonathan Frazer,
Richard Easther
Abstract:
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' initia…
▽ More
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.
△ Less
Submitted 20 December, 2014; v1 submitted 8 September, 2014;
originally announced September 2014.
-
Inflating an Inhomogeneous Universe
Authors:
Richard Easther,
Layne C. Price,
Javier Rasero
Abstract:
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 qualit…
▽ More
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.
△ Less
Submitted 19 August, 2014; v1 submitted 11 June, 2014;
originally announced June 2014.
-
The Knotted Sky II: Does BICEP2 require a nontrivial primordial power spectrum?
Authors:
Kevork N. Abazajian,
Grigor Aslanyan,
Richard Easther,
Layne C. Price
Abstract:
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…
▽ More
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.
△ Less
Submitted 26 August, 2014; v1 submitted 24 March, 2014;
originally announced March 2014.
-
The Knotted Sky I: Planck constraints on the primordial power spectrum
Authors:
Grigor Aslanyan,
Layne C. Price,
Kevork N. Abazajian,
Richard Easther
Abstract:
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…
▽ More
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.
△ Less
Submitted 26 August, 2014; v1 submitted 24 March, 2014;
originally announced March 2014.
-
Simple predictions from multifield inflationary models
Authors:
Richard Easther,
Jonathan Frazer,
Hiranya V. Peiris,
Layne C. Price
Abstract:
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…
▽ More
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.
△ Less
Submitted 8 April, 2014; v1 submitted 14 December, 2013;
originally announced December 2013.
-
Initial conditions and sampling for multifield inflation
Authors:
Richard Easther,
Layne C. Price
Abstract:
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, sinc…
▽ More
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.
△ Less
Submitted 17 July, 2013; v1 submitted 15 April, 2013;
originally announced April 2013.