-
Quantum initial conditions for curved inflating universes
Authors:
Mary I. Letey,
Zakhar Shumaylov,
Fruzsina J. Agocs,
Will J. Handley,
Michael P. Hobson,
Anthony N. Lasenby
Abstract:
We discuss the challenges of motivating, constructing, and quantizing a canonically normalized inflationary perturbation in spatially curved universes. We show that this has historically proved challenging due to the interaction of nonadiabaticity with spatial curvature. We construct a novel curvature perturbation that is canonically normalized in the sense of its equation of motion and is unique…
▽ More
We discuss the challenges of motivating, constructing, and quantizing a canonically normalized inflationary perturbation in spatially curved universes. We show that this has historically proved challenging due to the interaction of nonadiabaticity with spatial curvature. We construct a novel curvature perturbation that is canonically normalized in the sense of its equation of motion and is unique up to a single scalar parameter. With this construction it becomes possible to set initial conditions invariant under canonical transformations, overcoming known ambiguities in the literature. This corrected quantization has potentially observational consequences via modifications to the primordial power spectrum at large angular scales, as well as theoretical implications for quantization procedures in curved cosmologies filled with a scalar field.
△ Less
Submitted 8 July, 2024; v1 submitted 30 November, 2022;
originally announced November 2022.
-
Finite inflation in curved space
Authors:
Lukas T. Hergt,
Fruzsina J. Agocs,
Will J. Handley,
Michael P. Hobson,
Anthony N. Lasenby
Abstract:
We investigate the effects of non-zero spatial curvature on cosmic inflation in the light of cosmic microwave background (CMB) anisotropy measurements from the Planck 2018 legacy release and from the 2015 observing season of BICEP2 and the Keck Array. Even a small percentage of non-zero curvature today would significantly limit the total number of e-folds of the scale factor during inflation, rend…
▽ More
We investigate the effects of non-zero spatial curvature on cosmic inflation in the light of cosmic microwave background (CMB) anisotropy measurements from the Planck 2018 legacy release and from the 2015 observing season of BICEP2 and the Keck Array. Even a small percentage of non-zero curvature today would significantly limit the total number of e-folds of the scale factor during inflation, rendering just-enough inflation scenarios with a kinetically dominated or fast-roll stage prior to slow-roll inflation more likely. Finite inflation leads to oscillations and a cutoff towards large scales in the primordial power spectrum and curvature pushes them into the CMB observable window. Using nested sampling, we carry out Bayesian parameter estimations and model comparisons taking into account constraints from reheating and horizon considerations. We confirm the preference of CMB data for closed universes with Bayesian odds of over $100:1$ and with a posterior on the curvature density parameter of $Ω_{K,0}=-0.051\pm0.017$ for a curvature extension of LCDM and $Ω_{K,0}=-0.031\pm0.014$ for Starobinsky inflation. Model comparisons of various inflation models give similar results as for flat universes with the Starobinsky model outperforming most other models.
△ Less
Submitted 15 May, 2022;
originally announced May 2022.
-
Simple and statistically sound recommendations for analysing physical theories
Authors:
Shehu S. AbdusSalam,
Fruzsina J. Agocs,
Benjamin C. Allanach,
Peter Athron,
Csaba Balázs,
Emanuele Bagnaschi,
Philip Bechtle,
Oliver Buchmueller,
Ankit Beniwal,
Jihyun Bhom,
Sanjay Bloor,
Torsten Bringmann,
Andy Buckley,
Anja Butter,
José Eliel Camargo-Molina,
Marcin Chrzaszcz,
Jan Conrad,
Jonathan M. Cornell,
Matthias Danninger,
Jorge de Blas,
Albert De Roeck,
Klaus Desch,
Matthew Dolan,
Herbert Dreiner,
Otto Eberhardt
, et al. (50 additional authors not shown)
Abstract:
Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by mul…
▽ More
Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283.
△ Less
Submitted 11 April, 2022; v1 submitted 17 December, 2020;
originally announced December 2020.
-
Dense output for highly oscillatory numerical solutions
Authors:
F. J. Agocs,
M. P. Hobson,
W. J. Handley,
A. N. Lasenby
Abstract:
We present a method to construct a continuous extension (otherwise known as dense output) for a numerical routine in the special case of the numerical solution being a scalar-valued function exhibiting rapid oscillations. Such cases call for numerical routines that make use of the known global behaviour of the solution, one example being methods using asymptotic expansions to forecast the solution…
▽ More
We present a method to construct a continuous extension (otherwise known as dense output) for a numerical routine in the special case of the numerical solution being a scalar-valued function exhibiting rapid oscillations. Such cases call for numerical routines that make use of the known global behaviour of the solution, one example being methods using asymptotic expansions to forecast the solution at each step of the independent variable. An example is oscode, numerical routine which uses the Wentzel-Kramers-Brillouin (WKB) approximation when the solution oscillates rapidly and otherwise behaves as a Runge-Kutta (RK) solver. Polynomial interpolation is not suitable for producing the solution at an arbitrary point mid-step, since efficient numerical methods based on the WKB approximation will step through multiple oscillations in a single step. Instead we construct the continuous solution by extending the numerical quadrature used in computing a WKB approximation of the solution with no additional evaluations of the differential equation or terms within, and provide an error estimate on this dense output. Finally, we draw attention to previous work on the continuous extension of Runge-Kutta formulae, and construct an extension to a RK method based on Gauss--Lobatto quadrature nodes, thus describing how to generate dense output from each of the methods underlying oscode.
△ Less
Submitted 9 July, 2020;
originally announced July 2020.
-
Quantum initial conditions for inflation and canonical invariance
Authors:
F. J. Agocs,
L. T. Hergt,
W. J. Handley,
A. N. Lasenby,
M. P. Hobson
Abstract:
We investigate the transformation of initial conditions for primordial curvature perturbations under two types of transformations of the associated action: simultaneous redefinition of time and the field to be quantised, and the addition of surface terms. The latter encompasses all canonical transformations, whilst the time- and field-redefinition is a distinct, non-canonical transformation since…
▽ More
We investigate the transformation of initial conditions for primordial curvature perturbations under two types of transformations of the associated action: simultaneous redefinition of time and the field to be quantised, and the addition of surface terms. The latter encompasses all canonical transformations, whilst the time- and field-redefinition is a distinct, non-canonical transformation since the initial and destination systems use different times. Actions related to each other via such transformations yield identical equations of motion and preserve the commutator structure. They further preserve the time-evolution of expectation values of quantum operators unless the vacuum state also changes under the transformation. These properties suggest that it is of interest to investigate vacuum prescriptions that also remain unchanged under canonical transformations. We find that initial conditions derived via minimising the vacuum expectation value of the Hamiltonian and those obtained using the Danielsson vacuum prescription are not invariant under these transformations, whereas those obtained by minimising the local energy density are. We derive the range of physically distinct initial conditions obtainable by Hamiltonian diagonalisation, and illustrate their effect on the scalar primordial power spectrum and the Cosmic Microwave Background under the just enough inflation model. We also generalise the analogy between the dynamics of a quantum scalar field on a curved, time-dependent spacetime and the gauge-invariant curvature perturbation. We argue that the invariance of the vacuum prescription obtained by minimising the renormalised stress--energy tensor should make it the preferred procedure for setting initial conditions for primordial perturbations. All other procedures reviewed in this work yield ambiguous initial conditions, which is problematic both in theory and practice.
△ Less
Submitted 27 May, 2020; v1 submitted 14 February, 2020;
originally announced February 2020.
-
An efficient method for solving highly oscillatory ordinary differential equations with applications to physical systems
Authors:
F. J. Agocs,
W. J. Handley,
A. N. Lasenby,
M. P. Hobson
Abstract:
We present a novel numerical routine (oscode) with a C++ and Python interface for the efficient solution of one-dimensional, second-order, ordinary differential equations with rapidly oscillating solutions. The method is based on a Runge-Kutta-like stepping procedure that makes use of the Wentzel-Kramers-Brillouin (WKB) approximation to skip regions of integration where the characteristic frequenc…
▽ More
We present a novel numerical routine (oscode) with a C++ and Python interface for the efficient solution of one-dimensional, second-order, ordinary differential equations with rapidly oscillating solutions. The method is based on a Runge-Kutta-like stepping procedure that makes use of the Wentzel-Kramers-Brillouin (WKB) approximation to skip regions of integration where the characteristic frequency varies slowly. In regions where this is not the case, the method is able to switch to a made-to-measure Runge-Kutta integrator that minimises the total number of function evaluations. We demonstrate the effectiveness of the method with example solutions of the Airy equation and an equation exhibiting a burst of oscillations, discussing the error properties of the method in detail. We then show the method applied to physical systems. First, the one-dimensional, time-independent Schrödinger equation is solved as part of a shooting method to search for the energy eigenvalues for a potential with quartic anharmonicity. Then, the method is used to solve the Mukhanov-Sasaki equation describing the evolution of cosmological perturbations, and the primordial power spectrum of the perturbations is computed in different cosmological scenarios. We compare the performance of our solver in calculating a primordial power spectrum of scalar perturbations to that of BINGO, an efficient code specifically designed for such applications.
△ Less
Submitted 13 December, 2019; v1 submitted 30 May, 2019;
originally announced June 2019.
