Engineering Papers⌕ Search

Engineering topics

Wands, David

Publications and source records attributed to Wands, David.

At least 19 records

The stochastic spectator

We study the stochastic distribution of spectator fields predicted in different slow-roll inflation backgrounds. Spectator fields have a negligible energy density during inflation but may play an important dynamical role later, even giving rise to primordial density perturbations within our observational horizon today. During de-Sitter expansion there is an equilibrium solution for the spectator field which is often used to estimate the stochastic distribution during slow-roll inflation. However slow roll only requires that the Hubble rate varies slowly compared to the Hubble time, while the time taken for the stochastic distribution to evolve to the de-Sitter equilibrium solution can be much longer than a Hubble time. We study both chaotic (monomial) and plateau inflaton potentials, with quadratic, quartic and axionic spectator fields. We give an adiabaticity condition for the spectator field distribution to relax to the de-Sitter equilibrium, and find that the de-Sitter approximation is never a reliable estimate for the typical distribution at the end of inflation for a quadratic spectator during monomial inflation. The existence of an adiabatic regime at early times can erase the dependence on initial conditions of the final distribution of field values. In these cases, spectator fields acquire sub-Planckian expectation values. Otherwise spectator fields may acquire much larger field displacements than suggested by the de-Sitter equilibrium solution. We quantify the information about initial conditions that can be obtained from the final field distribution. Our results may have important consequences for the viability of spectator models for the origin of structure, such as the simplest curvaton models.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum diffusion during inflation and primordial black holes

We calculate the full probability density function (PDF) of inflationary curvature perturbations, even in the presence of large quantum backreaction. Making use of the stochastic-δ N formalism, two complementary methods are developed, one based on solving an ordinary differential equation for the characteristic function of the PDF, and the other based on solving a heat equation for the PDF directly. In the classical limit where quantum diffusion is small, we develop an expansion scheme that not only recovers the standard Gaussian PDF at leading order, but also allows us to calculate the first non-Gaussian corrections to the usual result. In the opposite limit where quantum diffusion is large, we find that the PDF is given by an elliptic theta function, which is fully characterised by the ratio between the squared width and height (in Planck mass units) of the region where stochastic effects dominate. We then apply these results to the calculation of the mass fraction of primordial black holes from inflation, and show that no more than ∼ 1 e -fold can be spent in regions of the potential dominated by quantum diffusion. We explain how this requirement constrains inflationary potentials with two examples.

79 ASTRONOMY AND ASTROPHYSICS↗

Relativistic initial conditions for N-body simulations

Initial conditions for (Newtonian) cosmological N-body simulations are usually set by re-scaling the present-day power spectrum obtained from linear (relativistic) Boltzmann codes to the desired initial redshift of the simulation. This back-scaling method can account for the effect of inhomogeneous residual thermal radiation at early times, which is absent in the Newtonian simulations. We analyse this procedure from a fully relativistic perspective, employing the recently-proposed Newtonian motion gauge framework. We find that N-body simulations for ΛCDM cosmology starting from back-scaled initial conditions can be self-consistently embedded in a relativistic space-time with first-order metric potentials calculated using a linear Boltzmann code. This space-time coincides with a simple ''N-body gauge'' for z < 50 for all observable modes. Care must be taken, however, when simulating non-standard cosmologies. As an example, we analyse the back-scaling method in a cosmology with decaying dark matter, and show that metric perturbations become large at early times in the back-scaling approach, indicating a breakdown of the perturbative description. We suggest a suitable ''forwards approach' for such cases.

79 ASTRONOMY AND ASTROPHYSICS↗

Relativistic interpretation of Newtonian simulations for cosmic structure formation

The standard numerical tools for studying non-linear collapse of matter are Newtonian N -body simulations. Previous work has shown that these simulations are in accordance with General Relativity (GR) up to first order in perturbation theory, provided that the effects from radiation can be neglected. In this paper we show that the present day matter density receives more than 1% corrections from radiation on large scales if Newtonian simulations are initialised before z =50. We provide a relativistic framework in which unmodified Newtonian simulations are compatible with linear GR even in the presence of radiation. Our idea is to use GR perturbation theory to keep track of the evolution of relativistic species and the relativistic space-time consistent with the Newtonian trajectories computed in N -body simulations. If metric potentials are sufficiently small, they can be computed using a first-order Einstein–Boltzmann code such as CLASS. We make this idea rigorous by defining a class of GR gauges, the Newtonian motion gauges, which are defined such that matter particles follow Newtonian trajectories. We construct a simple example of a relativistic space-time within which unmodified Newtonian simulations can be interpreted.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining curvatonic reheating

We derive the first systematic observational constraints on reheating in models of inflation where an additional light scalar field contributes to primordial density perturbations and affects the expansion history during reheating. This encompasses the original curvaton model but also covers a larger class of scenarios. We find that, compared to the single-field case, lower values of the energy density at the end of inflation and of the reheating temperature are preferred when an additional scalar field is introduced. For instance, if inflation is driven by a quartic potential, which is one of the most favoured models when a light scalar field is added, the upper bound T {sub reh} < 5 × 10{sup 4} GeV on the reheating temperature T {sub reh} is derived, and the implications of this value on post-inflationary physics are discussed. The information gained about reheating is also quantified and it is found that it remains modest in plateau inflation (though still larger than in the single-field version of the model) but can become substantial in quartic inflation. The role played by the vev of the additional scalar field at the end of inflation is highlighted, and opens interesting possibilities for exploring stochastic inflation effects that could determine its distribution.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiple fields in stochastic inflation

Stochastic effects in multi-field inflationary scenarios are investigated. A hierarchy of diffusion equations is derived, the solutions of which yield moments of the numbers of inflationary e-folds. Solving the resulting partial differential equations in multi-dimensional field space is more challenging than the single-field case. A few tractable examples are discussed, which show that the number of fields is, in general, a critical parameter. When more than two fields are present for instance, the probability to explore arbitrarily large-field regions of the potential, otherwise inaccessible to single-field dynamics, becomes non-zero. In some configurations, this gives rise to an infinite mean number of e-folds, regardless of the initial conditions. Another difference with respect to single-field scenarios is that multi-field stochastic effects can be large even at sub-Planckian energy. This opens interesting new possibilities for probing quantum effects in inflationary dynamics, since the moments of the numbers of e-folds can be used to calculate the distribution of primordial density perturbations in the stochastic-δN formalism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Galaxy bispectrum, primordial non-Gaussianity and redshift space distortions

Measurements of the non-Gaussianity of the primordial density field have the power to considerably improve our understanding of the physics of inflation. Indeed, if we can increase the precision of current measurements by an order of magnitude, a null-detection would rule out many classes of scenarios for generating primordial fluctuations. Large-scale galaxy redshift surveys represent experiments that hold the promise to realise this goal. Thus, we model the galaxy bispectrum and forecast the accuracy with which it will probe the parameter f {sub NL}, which represents the degree of primordial local-type non Gaussianity. Specifically, we address the problem of modelling redshift space distortions (RSD) in the tree-level galaxy bispectrum including f {sub NL}. We find novel contributions associated with RSD, with the characteristic large scale amplification induced by local-type non-Gaussianity. These RSD effects must be properly accounted for in order to obtain un-biased measurements of f {sub NL} from the galaxy bispectrum. We propose an analytic template for the monopole which can be used to fit against data on large scales, extending models used in the recent measurements. Finally, we perform idealised forecasts on σ {sub f} {sub N{sub L}}—the accuracy of the determination of local non-linear parameter f {sub NL}—from measurements of the galaxy bispectrum. Our findings suggest that current surveys can in principle provide f {sub NL} constraints competitive with Planck , and future surveys could improve them further.

79 ASTRONOMY AND ASTROPHYSICS↗

The intrinsic matter bispectrum in ΛCDM

We present a fully relativistic calculation of the matter bispectrum at second order in cosmological perturbation theory assuming a Gaussian primordial curvature perturbation. For the first time we perform a full numerical integration of the bispectrum for both baryons and cold dark matter using the second-order Einstein-Boltzmann code, SONG. We review previous analytical results and provide an improved analytic approximation for the second-order kernel in Poisson gauge which incorporates Newtonian nonlinear evolution, relativistic initial conditions, the effect of radiation at early times and the cosmological constant at late times. Our improved kernel provides a percent level fit to the full numerical result at late times for most configurations, including both equilateral shapes and the squeezed limit. We show that baryon acoustic oscillations leave an imprint in the matter bispectrum, making a significant impact on squeezed shapes.

79 ASTRONOMY AND ASTROPHYSICS↗

Inflation with an extra light scalar field after Planck

Bayesian inference techniques are used to investigate situations where an additional light scalar field is present during inflation and reheating. This includes (but is not limited to) curvaton-type models. We design a numerical pipeline where ≅ 200 inflaton setups × 10 reheating scenarios = 2000 models are implemented and we present the results for a few prototypical potentials. We find that single-field models are remarkably robust under the introduction of light scalar degrees of freedom. Models that are ruled out at the single-field level are not improved in general, because good values of the spectral index and the tensor-to-scalar ratio can only be obtained for very fine-tuned values of the extra field parameters and/or when large non-Gaussianities are produced. The only exception is quartic large-field inflation, so that the best models after Planck are of two kinds: plateau potentials, regardless of whether an extra field is added or not, and quartic large-field inflation with an extra light scalar field, in some specific reheating scenarios. Using Bayesian complexity, we also find that more parameters are constrained for the models we study than for their single-field versions. This is because the added parameters not only contribute to the reheating kinematics but also to the cosmological perturbations themselves, to which the added field contributes. The interplay between these two effects lead to a suppression of degeneracies that is responsible for having more constrained parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Encyclopædia curvatonis

We investigate whether the predictions of single-field models of inflation are robust under the introduction of additional scalar degrees of freedom, and whether these extra fields change the potentials for which the data show the strongest preference. We study the situation where an extra light scalar field contributes both to the total curvature perturbations and to the reheating kinematic properties. Ten reheating scenarios are identified, and all necessary formulas allowing a systematic computation of the predictions for this class of models are derived. They are implemented in the public library ASPIC, which contains more than 75 single-field potentials. This paves the way for a forthcoming full Bayesian analysis of the problem. A few representative examples are displayed and discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-local bias in the halo bispectrum with primordial non-Gaussianity

Primordial non-Gaussianity can lead to a scale-dependent bias in the density of collapsed halos relative to the underlying matter density. The galaxy power spectrum already provides constraints on local-type primordial non-Gaussianity complementary those from the cosmic microwave background (CMB), while the bispectrum contains additional shape information and has the potential to outperform CMB constraints in future. We develop the bias model for the halo density contrast in the presence of local-type primordial non-Gaussianity, deriving a bivariate expansion up to second order in terms of the local linear matter density contrast and the local gravitational potential in Lagrangian coordinates. Nonlinear evolution of the matter density introduces a non-local tidal term in the halo model. Furthermore, the presence of local-type non-Gaussianity in the Lagrangian frame leads to a novel non-local convective term in the Eulerian frame, that is proportional to the displacement field when going beyond the spherical collapse approximation. We use an extended Press-Schechter approach to evaluate the halo mass function and thus the halo bispectrum. We show that including these non-local terms in the halo bispectra can lead to corrections of up to 25% for some configurations, on large scales or at high redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

CMB hemispherical asymmetry from non-linear isocurvature perturbations

We investigate whether non-adiabatic perturbations from inflation could produce an asymmetric distribution of temperature anisotropies on large angular scales in the cosmic microwave background (CMB). We use a generalised non-linear δ N formalism to calculate the non-Gaussianity of the primordial density and isocurvature perturbations due to the presence of non-adiabatic, but approximately scale-invariant field fluctuations during multi-field inflation. This local-type non-Gaussianity leads to a correlation between very long wavelength inhomogeneities, larger than our observable horizon, and smaller scale fluctuations in the radiation and matter density. Matter isocurvature perturbations contribute primarily to low CMB multipoles and hence can lead to a hemispherical asymmetry on large angular scales, with negligible asymmetry on smaller scales. In curvaton models, where the matter isocurvature perturbation is partly correlated with the primordial density perturbation, we are unable to obtain a significant asymmetry on large angular scales while respecting current observational constraints on the observed quadrupole. However in the axion model, where the matter isocurvature and primordial density perturbations are uncorrelated, we find it may be possible to obtain a significant asymmetry due to isocurvature modes on large angular scales. Such an isocurvature origin for the hemispherical asymmetry would naturally give rise to a distinctive asymmetry in the CMB polarisation on large scales.

79 ASTRONOMY AND ASTROPHYSICS↗

Generalised tensor fluctuations and inflation

Using an effective field theory approach to inflation, we examine novel properties of the spectrum of inflationary tensor fluctuations, that arise when breaking some of the symmetries or requirements usually imposed on the dynamics of perturbations. During single-clock inflation, time-reparameterization invariance is broken by a time-dependent cosmological background. In order to explore more general scenarios, we consider the possibility that spatial diffeomorphism invariance is also broken by effective mass terms or by derivative operators for the metric fluctuations in the Lagrangian. We investigate the cosmological consequences of the breaking of spatial diffeomorphisms, focussing on operators that affect the power spectrum of fluctuations. We identify the operators for tensor fluctuations that can provide a blue spectrum without violating the null energy condition, and operators for scalar fluctuations that lead to non-conservation of the comoving curvature perturbation on superhorizon scales even in single-clock inflation. In the last part of our work, we also examine the consequences of operators containing more than two spatial derivatives, discussing how they affect the sound speed of tensor fluctuations, and showing that they can mimic some of the interesting effects of symmetry breaking operators, even in scenarios that preserve spatial diffeomorphism invariance.

79 ASTRONOMY AND ASTROPHYSICS↗

EINSTEIN'S SIGNATURE IN COSMOLOGICAL LARGE-SCALE STRUCTURE

We show how the nonlinearity of general relativity generates a characteristic nonGaussian signal in cosmological large-scale structure that we calculate at all perturbative orders in a large-scale limit. Newtonian gravity and general relativity provide complementary theoretical frameworks for modeling large-scale structure in ΛCDM cosmology; a relativistic approach is essential to determine initial conditions, which can then be used in Newtonian simulations studying the nonlinear evolution of the matter density. Most inflationary models in the very early universe predict an almost Gaussian distribution for the primordial metric perturbation, ζ. However, we argue that it is the Ricci curvature of comoving-orthogonal spatial hypersurfaces, R, that drives structure formation at large scales. We show how the nonlinear relation between the spatial curvature, R, and the metric perturbation, ζ, translates into a specific nonGaussian contribution to the initial comoving matter density that we calculate for the simple case of an initially Gaussian ζ. Our analysis shows the nonlinear signature of Einstein's gravity in large-scale structure.

79 ASTRONOMY AND ASTROPHYSICS↗

The intrinsic B-mode polarisation of the Cosmic Microwave Background

We estimate the B-polarisation induced in the Cosmic Microwave Background by the non-linear evolution of density perturbations. Using the second-order Boltzmann code SONG, our analysis incorporates, for the first time, all physical effects at recombination. We also include novel contributions from the redshift part of the Boltzmann equation and from the bolometric definition of the temperature in the presence of polarisation. The remaining line-of-sight terms (lensing and time-delay) have previously been studied and must be calculated non-perturbatively. The intrinsic B-mode polarisation is present independent of the initial conditions and might contaminate the signal from primordial gravitational waves. We find this contamination to be comparable to a primordial tensor-to-scalar ratio of r ≅ 10{sup −7} at the angular scale ℓ ≅ 100, where the primordial signal peaks, and r ≅ 5 × 10{sup −5} at ℓ ≅ 700, where the intrinsic signal peaks. Therefore, we conclude that the intrinsic B-polarisation from second-order effects is not likely to contaminate future searches of primordial gravitational waves.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-Gaussian initial conditions in ΛCDM: Newtonian, relativistic, and primordial contributions

The goal of the present paper is to set initial conditions for structure formation at nonlinear order, consistent with general relativity, while also allowing for primordial non-Gaussianity. We use the nonlinear continuity and Raychaudhuri equations, which together with the nonlinear energy constraint, determine the evolution of the matter density fluctuation in general relativity. We solve this equations at first and second order in a perturbative expansion, recovering and extending previous results derived in the matter-dominated limit and in the Newtonian regime. We present a second-order solution for the comoving density contrast in a ΛCDM universe, identifying nonlinear contributions coming from the Newtonian growing mode, primordial non-Gaussianity and intrinsic non-Gaussianity, due to the essential nonlinearity of the relativistic constraint equations. We discuss the application of these results to initial conditions in N-body simulations, showing that relativistic corrections mimic a non-zero nonlinear parameter f {sub NL}.

79 ASTRONOMY AND ASTROPHYSICS↗