Search

We study linear scalar perturbations in single-field models of inflation featuring a non-attractor phase. These models lead to a peak in the curvature power spectrum that may result in the formation of primordial black holes. We develop a transfer-matrix formalism, analogous to the S-matrix program in quantum-field theory, that maps perturbations throughout the transitory phase. At scales smaller than the peak, the power spectrum features damped oscillations, and the duration of the transition sets the scale at which power-law damping switches to exponential damping. At scales larger than the peak, we demonstrate that a dip appears in the power spectrum if and only if the inflaton's velocity does not flip sign. We show that the amplitude at the dip always scales as the inverse square-rooted amplitude of the peak, and comment on the physical consequences of this universal relationship. We also test the robustness of our results with a few toy models and interpret them with an intuitive mechanical analogy., Comment: 31 pages without appendices (total 40 pages), 6 figures

We introduce a novel framework to implement stochastic inflation on stochastic trees, modelling the inflationary expansion as a branching process. Combined with the $\delta N$ formalism, this allows us to generate real-space maps of the curvature perturbation that fully capture quantum diffusion and its non-perturbative backreaction during inflation. Unlike lattice methods, trees do not proceed on a fixed background since new spacetime units emerge dynamically as trees unfold, naturally incorporating metric fluctuations. The recursive structure of stochastic trees also offers remarkable numerical efficiency, and we develop the FOrtran Recursive Exploration of Stochastic Trees ($\texttt{FOREST}$) tool and demonstrate its performance. We show how primordial black holes blossom at unbalanced nodes of the trees, and how their mass distribution can be obtained while automatically accounting for the "cloud-in-cloud" effect. In the "quantum-well" toy model, we find broad mass distributions, with mild power laws terminated by exponential tails. We finally compare our results with existing approximations in the literature and discuss several prospects., Comment: 39 pages, 13 figures

We simulate the distribution of very rare, large excursions in the primordial density field produced in models of inflation in the very early universe which include a strong enhancement of the power spectrum. The stochastic $\delta \mathcal{N}$ formalism is used to identify the probability distribution for the primordial curvature perturbation with the first-passage-time distribution, $P(\delta \mathcal{N})$, and we compare our stochastic results with those obtained in the classical $\delta \mathcal{N}$ approach. We extend the PyFPT numerical code to simulate the full 2D phase space, and apply importance sampling which allows very rare fluctuations to be simulated in $\mathcal{O}(10)$ minutes on a single CPU, where previous direct simulations required supercomputers. We demonstrate that the stochastic noise due to quantum fluctuations after a sudden transition to ultra-slow roll can be accurately modelled using an analytical Bessel-function ansatz to identify the homogeneous growing mode. The stochastic noise found in this way is a function of the field value only. This enables us to coarse grain the inflation field at the Hubble scale and include non-linear, stochastic evolution on all super-Hubble length scales., Comment: 25 pages plus appendix. 8 figures. To be submitted to JCAP

We show how to compute the purity and entanglement entropy for quantum fields in a systematic perturbative expansion. To that end, we generalize the in-in formalism to non-unitary dynamics (i.e. accounting for the presence of an environment) and to the calculation of quantum information measures, which are not observables in the usual sense. This allows us to reduce the problem to one involving standard correlation functions, and to organize their computation in a diagrammatic expansion for which we construct the corresponding Feynman rules. As an illustration, we apply the formalism to a cosmological setting inspired by the effective field theory of inflation. We find that at late times, non-linear loop corrections share the same time behavior as the linear contribution, and only yield a slight redressing of the purity. In particular, when the environment is heavy compared to the Hubble scale, the phenomenon of recoherence previously encountered is robust to the class of non-linear extensions considered. Bridging the gap between perturbative quantum field theory and open quantum systems paves the way to a better understanding of renormalization and resummation in open effective field theories. It also enables a more systematic exploration of quantum information properties in field theoretic settings., Comment: 43 pages without appendices (66 pages in total), 14 figures, matches published version in JCAP

We show that the simplest, and currently favoured, theoretical realizations of cosmic inflation yield a sharp prediction for the running of the spectral index $\alpha_\mathrm{S}$. Using latest cosmological data, we compute its marginalized posterior probability distribution over the space of nearly 300 models of single-field slow-roll inflation. The most probable value is $\alpha_\mathrm{S}=-6.3 \times 10^{-4}$, lying within the $98\%$ credible interval $-1.8 \times 10^{-3}< \alpha_\mathrm{S}< -9.1 \times 10^{-5}$. Within the landscape of all the proposed slow-roll inflationary models, positive values for the running are therefore disfavoured at more than three-sigma., Comment: 6 pages, 4 figures, uses epl2. Misprints corrected, matches published version

The capability of Cosmic Inflation to explain the latest Cosmic Microwave Background and Baryonic Acoustic Oscillation data is assessed by performing Bayesian model comparison within the landscape of nearly three-hundred models of single-field slow-roll inflation. We present the first Bayesian data analysis based on the third-order slow-roll primordial power spectra. In particular, the fourth Hubble-flow function $\epsilon_4$ remains unbounded while the third function verifies, at two-sigma, $\epsilon_{3}\in[-0.4,0.5]$, which is perfectly compatible with the slow-roll predictions for the running of the spectral index. We also observe some residual excess of $B$-modes within the BICEP/Keck data favoring, at a non-statistically significant level, non-vanishing primordial tensor modes: $\log(\epsilon_{1}) > -3.9$, at $68\%$ confidence level. Then, for 283 models of single-field inflation, we compute the Bayesian evidence, the Bayesian dimensionality and the marginalized posteriors of all the models' parameters, including the ones associated with the reheating era. The average information gain on the reheating parameter $R_\mathrm{reh}$ reaches $1.3 \pm 0.18$ bits, which is more than a factor two improvement compared to the first Planck data release. As such, inflationary model predictions cannot meet data accuracy without specifying, or marginalizing over, the reheating kinematics. We also find that more than $40\%$ of the scenarios are now strongly disfavored, which shows that the constraining power of cosmological data is winning against the increase of the number of proposed models. In addition, about $20\%$ of all models have evidences within the most probable region and are all favored according to the Jeffreys' scale of Bayesian evidences., Comment: 79 pages, 37 figures, uses jcappub. Discussions added, matches published version

We compute the rate with which unobserved fields decohere other fields to which they couple, both in flat space and in de Sitter space, for spectator scalar fields prepared in their standard adiabatic vacuum. The process is very efficient in de Sitter space once the modes in question pass outside the Hubble scale, displaying the tell-tale phenomenon of secular growth that indicates the breakdown of perturbative methods on a time scale parameterically long compared with the Hubble time. We show how to match the perturbative evolution valid at early times onto a late-time Lindblad evolution whose domain of validity extends to much later times, thereby allowing a reliable resummation of the perturbative result beyond the perturbative regime. Super-Hubble modes turn out to be dominantly decohered by unobserved modes that are themselves also super-Hubble. Although our calculation is done for spectator fields, if applied to curvature perturbations during inflation our observations here could close a potential loophole in recent calculations of the late-time purity of the observable primordial fluctuations., Comment: 37 pages without appendices (59 pages in total), 7 figures; matches published version in JCAP

Quantum diffusion describes the inflow of vacuum quantum fluctuations as they get amplified by gravitational instability, and stretched to large distances during inflation. In this picture, the dynamics of the universe's expansion becomes stochastic, and the statistics of the curvature perturbation is encoded in the distribution of the duration of inflation. This provides a non-perturbative framework to study cosmological fluctuations during inflation, which is well-suited to the case of primordial black holes since they originate from large fluctuations. We show that standard, perturbative expectations for the primordial black hole abundance can be significantly modified by quantum-diffusion effects, and we identify a few open challenges., Comment: Invited chapter to the book "Primordial Black Holes'', Springer 2024, ed. Chris Byrnes, Gabriele Franciolini, Tomohiro Harada, Paolo Pani, Misao Sasaki. 22 pages, 3 figures, accepted version

We study how large fluctuations are spatially correlated in the presence of quantum diffusion during inflation. This is done by computing real-space correlation functions in the stochastic-$\delta N$ formalism. We first derive an exact description of physical distances as measured by a local observer at the end of inflation, improving on previous works. Our approach is based on recursive algorithmic methods that consistently include volume-weighting effects. We then propose a "large-volume'' approximation under which calculations can be done using first-passage time analysis only, and from which a new formula for the power spectrum in stochastic inflation is derived. We then study the full two-point statistics of the curvature perturbation. Due to the presence of exponential tails, we find that the joint distribution of large fluctuations is of the form $P(\zeta_{R_1}, \zeta_{R_2}) = F(R_1,R_2,r) P(\zeta_{R_1})P( \zeta_{R_2})$, where $\zeta_{R_1}$ and $\zeta_{R_2}$ denote the curvature perturbation coarse-grained at radii $R_1$ and $R_2$, around two spatial points distant by $r$. This implies that, on the tail, the reduced correlation function, defined as $P(\zeta_{R_1}>\zeta_{\rm{c}}, \zeta_{R_2}>\zeta_{\rm{c}})/[P(\zeta_{R_1}>\zeta_{\rm{c}}) P(\zeta_{R_2}>\zeta_{\rm{c}})]-1$, is independent of the threshold value $\zeta_{\rm{c}}$. This contrasts with Gaussian statistics where the same quantity strongly decays with $\zeta_{\rm{c}}$, and shows the existence of a universal clustering profile for all structures forming in the exponential tails. Structures forming in the intermediate (i.e. not yet exponential) tails may feature different, model-dependent behaviours., Comment: 52 pages, 14 figures. Discussion around Eq.(3.6) expanded, a few minor changes and typos fixed. Some references added. Matches published version in JCAP

If dark matter is made of QCD axions, its abundance is determined by the vacuum expectation value acquired by the axion field during inflation. The axion is usually assumed to follow the equilibrium distribution arising from quantum diffusion during inflation. This leads to the so-called stochastic window under which the QCD axion can make up all the dark matter. It is characterised by $10^{10.4}\mathrm{GeV}\leq f\leq 10^{17.2}\mathrm{GeV}$ and $H_{\mathrm{end}}>10^{-2.2}\mathrm{GeV}$, where $f$ is the axion decay constant and $H_{\mathrm{end}}$ is the Hubble expansion rate at the end of inflation. However, in realistic inflationary potentials, we show that the axion never reaches the equilibrium distribution at the end of inflation. This is because the relaxation time of the axion is much larger than the typical time scale over which $H$ varies during inflation. As a consequence, the axion acquires a quasi-flat distribution as long as it remains light during inflation. This leads us to reassessing the stochastic axion window, and we find that $ 10^{10.3}\mathrm{GeV}\leq f\leq 10^{14.1}\mathrm{GeV}$ and $H_{\mathrm{end}}>10^{-13.8}\mathrm{GeV}$., Comment: 28 pages, 8 figures, minor changes, matches the published version in JCAP

The separate-universe approach gives an intuitive way to understand the evolution of cosmological perturbations in the long-wavelength limit. It uses solutions of the spatially-homogeneous equations of motion to model the evolution of the inhomogeneous universe on large scales. We show that the separate-universe approach fails on a finite range of super-Hubble scales at a sudden transition from slow roll to ultra-slow roll during inflation in the very early universe. Such transitions are a feature of inflation models giving a large enhancement in the primordial power spectrum on small scales, necessary to produce primordial black holes after inflation. We show that the separate-universe approach still works in a piece-wise fashion, before and after the transition, but spatial gradients on finite scales require a discontinuity in the homogeneous solution at the transition. We discuss the implications for the $\delta N$ formalism and stochastic inflation, which employ the separate-universe approximation., Comment: 24 pages plus appendices. 9 figures. v2: to appear in JCAP

In the recent years, primordial black holes (PBHs) have emerged as one of the most interesting and hotly debated topics in cosmology. Among other possibilities, PBHs could explain both some of the signals from binary black hole mergers observed in gravitational wave detectors and an important component of the dark matter in the Universe. Significant progress has been achieved both on the theory side and from the point of view of observations, including new models and more accurate calculations of PBH formation, evolution, clustering, merger rates, as well as new astrophysical and cosmological probes. In this work, we review, analyse and combine the latest developments in order to perform end-to-end calculations of the various gravitational wave signatures of PBHs. Different ways to distinguish PBHs from stellar black holes are emphasized. Finally, we discuss their detectability with LISA, the first planned gravitational-wave observatory in space., Comment: 161 pages, 47 figures, comments welcome

Most high-energy constructions that realise a phase of cosmic inflation contain many degrees of freedom. Yet, cosmological observations are all consistent with single-field embeddings. We show how volume selection effects explain this apparent paradox. Because of quantum diffusion, different regions of space inflate by different amounts. In regions that inflate most, and eventually dominate the volume of the universe, a generic mechanism is unveiled that diverts the inflationary dynamics towards single-field attractors. The formalism of constrained stochastic inflation is developed to this end., Comment: 5 pages and 2 figures without appendices (total 17 pages, 7 figures)

Cosmological perturbation theory is an example of a gauge theory, where gauge transformations correspond to changes in the space-time coordinate system. To determine physical quantities, one is free to introduce gauge conditions (\ie to work with specific space-time coordinates), and such conditions are often used to simplify technical aspects of the calculation or to facilitate the interpretation of the physical degrees of freedom. Some of the prescriptions introduced in the literature are known to fix the gauge only partially, but it is commonly assumed that the remaining gauge degrees of freedom can be fixed somehow. In this work, we show that this is not necessarily the case, and that some of these gauges are indeed pathological. We derive a systematic procedure to determine whether a gauge is pathological or not, and to complete partially-fixed gauges into healthy gauges when this is possible. In this approach, the Lagrange multipliers (\ie the perturbed lapse and shift in the ADM formalism) cannot appear in the off-shell definition of the gauges, they necessarily arise as on-shell consequences of the gauge conditions. As illustrative applications, we propose an alternative, non-pathological formulation of the synchronous gauge, and we show that the uniform-expansion gauge (as well as any gauge ensuring vanishing lapse perturbations) can hardly be made healthy. Our methodology also allows us to construct all gauge-invariant variables. We further show that our non-pathological criterion for gauges is also the one that ensures Dirac brackets to be properly defined. This allows cosmological perturbations to be quantised in a gauge-fixed way. We finally discuss possible generalisations of our formalism., Comment: 40 pages without appendices, total 69 pages, 1 figure. Matches published version in JCAP. Link with Dirac-bracket quantisation introduced around Eq. (6.7). Abstract modified accordingly

The aim of this paper is to highlight the challenges and potential gains surrounding a coherent description of physics from the high-energy scales of inflation down to the lower energy scales probed in particle-physics experiments. As an example, we revisit the way inflation can be realised within an effective Minimal Supersymmetric Standard Model (eMSSM), in which the $LLe$ and $udd$ flat directions are lifted by the combined effect of soft-supersymmetric-breaking masses already present in the MSSM, together with the addition of effective non-renormalizable operators. We clarify some features of the model and address the question of the one-loop Renormalization Group improvement of the inflationary potential, discussing its impact on the fine-tuning of the model. We also compare the parameter space that is compatible with current observations (in particular the amplitude, $A_{\scriptscriptstyle{\mathrm{S}}}$, and the spectral index, $n_{\scriptscriptstyle{\mathrm{S}}}$, of the primordial cosmological fluctuations) at tree level and at one loop, and discuss the role of reheating. Finally we perform combined fits of particle and cosmological observables (mainly $A_{\scriptscriptstyle{\mathrm{S}}}$, $n_{\scriptscriptstyle{\mathrm{S}}}$, the Higgs mass, and the cold-dark-matter energy density) with the one-loop inflationary potential applied to some examples of dark-matter annihilation channels (Higgs-funnel, Higgsinos and A-funnel), and discuss the status of the ensuing MSSM spectra with respect to the LHC searches., Comment: 38 pages, 7 figures. Includes updated bibliography and minor additions in the text. Matches the published version in PRD

We revisit how super-Hubble cosmological fluctuations induce, at any time in the cosmic history, a non-vanishing spatial curvature of the local background metric. The random nature of these fluctuations promotes the curvature density parameter to a stochastic quantity for which we derive novel non-perturbative expressions for its mean, variance, higher moments and full probability distribution. For scale-invariant Gaussian perturbations, such as those favored by cosmological observations, we find that the most probable value for the curvature density parameter $\Omega_\mathrm{K}$ today is $-10^{-9}$, that its mean is $+10^{-9}$, both being overwhelmed by a standard deviation of order $10^{-5}$. We then discuss how these numbers would be affected by the presence of large super-Hubble non-Gaussianities, or, if inflation lasted for a very long time. In particular, we find that substantial values of $\Omega_\mathrm{K}$ are obtained if inflation lasts for more than a billion e-folds., Comment: 7 pages, 2 figures, uses RevTex. Misprints corrected, references added, matches published version

Primordial black holes (PBH) may form from large cosmological perturbations, produced during inflation when the inflaton's velocity is sufficiently slowed down. This usually requires very flat regions in the inflationary potential. In this paper we investigate another possibility, namely that the inflaton climbs up its potential. When it turns back, its velocity crosses zero, which triggers a short phase of ``uphill inflation'' during which cosmological perturbations grow at a very fast rate. This naturally occurs in double-well potentials if the width of the well is close to the Planck scale. We include the effect of quantum diffusion in this scenario, which plays a crucial role, by means of the stochastic-$\delta N$ formalism. We find that ultra-light black holes are produced with very high abundances, which do not depend on the energy scale at which uphill inflation occurs, and which suffer from substantially less fine tuning than in alternative PBH-production models. They are such that PBHs later drive a phase of PBH domination., Comment: matches published version in JCAP. 25 pages plus appendices (total 33 pages), 12 figures

Despite being created through a fundamentally quantum-mechanical process, cosmological structures have not yet revealed any sign of genuine quantum correlations. Among the obstructions to the direct detection of quantum signatures in cosmology, environmental-induced decoherence is arguably one of the most inevitable. Yet, we discover a mechanism of quantum recoherence for the adiabatic perturbations when they couple to an entropic sector. After a transient phase of decoherence, a turning point is reached, recoherence proceeds and adiabatic perturbations exhibit a large amount of self-coherence at late-time. This result is also understood by means of a non-Markovian master equation, which reduces to Wilsonian effective-field theory in the unitary limit. This allows us to critically assess the validity of open-quantum-system methods in cosmology and to highlight that re(de)coherence from linear interactions has no flat-space analogue., Comment: 5 + 12 pages, 4 + 2 figures, matches published version in EPL

We compute the rate with which super-Hubble cosmological fluctuations are decohered during inflation, by their gravitational interactions with unobserved shorter-wavelength scalar and tensor modes. We do so using Open Effective Field Theory methods, that remain under control at the late times of observational interest, contrary to perturbative calculations. Our result is minimal in the sense that it only incorporates the self-interactions predicted by General Relativity in single-clock models (additional interaction channels should only speed up decoherence). We find that decoherence is both suppressed by the first slow-roll parameter and by the energy density during inflation in Planckian units, but that it is enhanced by the volume comprised within the scale of interest, in Hubble units. This implies that, for the scales probed in the Cosmic Microwave Background, decoherence is effective as soon as inflation proceeds above $\sim 5\times 10^{9}$ GeV. Alternatively, if inflation proceeds at GUT scale decoherence is incomplete only for the scales crossing out the Hubble radius in the last ~ 13 e-folds, of inflation. We also compute how short-wavelength scalar modes decohere primordial tensor perturbations, finding a faster rate unsuppressed by slow-roll parameters. Identifying the parametric dependence of decoherence, and the rate at which it proceeds, helps suggest ways to look for quantum effects., Comment: 31 pages + appendices, 7 figures

Measuring the quantumness of a system can be done with a variety of methods. In this article we compare different criteria, namely quantum discord, Bell inequality violation and non-separability, for systems placed in a Gaussian state. When the state is pure, these criteria are equivalent, while we find that they do not necessarily coincide when decoherence takes place. Finally, we prove that these criteria are essentially controlled by the semi-minor axis of the ellipse representing the state's Wigner function in phase space., Comment: 8 pages without appendix (total 19 pages), 5 figures. Matches published version in EPL

We present a non-perturbative method for calculating the abundance of primordial black holes given an arbitrary one-point probability distribution function for the primordial curvature perturbation, $P(\zeta)$. A non-perturbative method is essential when considering non-Gaussianities that cannot be treated using a conventional perturbative expansion. To determine the full statistics of the density field, we relate $\zeta$ to a Gaussian field by equating the cumulative distribution functions. We consider two examples: a specific local-type non-Gaussian distribution arising from ultra slow roll models, and a general piecewise model for $P(\zeta)$ with an exponential tail. We demonstrate that the enhancement of primordial black hole formation is due to the intermediate regime, rather than the far tail. We also show that non-Gaussianity can have a significant impact on the shape of the primordial black hole mass distribution., Comment: 9 pages, 2 figures; Version 3: matches published version

If the inflaton gets trapped in a local minimum of its potential shortly before the end of inflation, it escapes by building up quantum fluctuations in a process known as stochastic tunnelling. In this work we study cosmological fluctuations produced in such a scenario, and how likely they are to form Primordial Black Holes (PBHs). This is done by using the stochastic-$\delta N$ formalism, which allows us to reconstruct the highly non-Gaussian tails of the distribution function of the number of $e$-folds spent in the false-vacuum state. We explore two different toy models, both analytically and numerically, in order to identify which properties do or do not depend on the details of the false-vacuum profile. We find that when the potential barrier is small enough compared to its width, $\Delta V/V < \Delta\phi^2/M_{\text{Pl}}^2$, the potential can be approximated as being flat between its two local extrema, so results previously obtained in a "flat quantum well'' apply. Otherwise, when $\Delta V/V < V/M_{\text{Pl}}^4$, the PBH abundance depends exponentially on the height of the potential barrier, and when $\Delta V/V > V/M_{\text{Pl}}^4$ it depends super-exponentially ( i.e. as the exponential of an exponential) on the barrier height. In that later case PBHs are massively produced. This allows us to quantify how much flat inflection points need to be fine-tuned. In a deep false vacuum, we also find that slow-roll violations are typically encountered unless the potential is close to linear. This motivates further investigations to generalise our approach to non-slow-roll setups., Comment: 25 pages without appendix ( 36 pages in total ), 13 figures. Matches the published version in JCAP

Master equations are commonly employed in cosmology to model the effect of additional degrees of freedom, treated as an "environment", onto a given "system". However, they rely on assumptions that are not necessarily satisfied in cosmology, where the environment may be out of equilibrium and the background is dynamical. In this work, we apply the master-equation program to a model that is exactly solvable, and which consists of two linearly coupled scalar fields evolving on a cosmological background. The light field plays the role of the system and the heavy field is the environment. By comparing the exact solution to the output of the master equation, we can critically assess its performance. We find that the master equation exhibits a set of "spurious" terms that explicitly depend on the initial conditions, and which arise as a consequence of working on a dynamical background. Although they cancel out in the perturbative limit of the theory (i.e. at leading orders in the interaction strength), they spoil resummation. However, when those terms are removed, the master equation performs impressively well to reproduce the power spectra and the amount of the decoherence of the light field, even in the strongly decohered regime. We conclude that master equations are able to perform late-time resummation, even though the system is far from the Markovian limit, provided spurious contributions are suppressed., Comment: 25 pages without appendices (48 pages in total), 5 figures, matches published version in EPJC

It is generally assumed within the standard cosmological model that initial density perturbations are Gaussian at all scales. However, primordial quantum diffusion unavoidably generates non-Gaussian, exponential tails in the distribution of inflationary perturbations. These exponential tails have direct consequences for the formation of collapsed structures in the universe, as has been studied in the context of primordial black holes. We show that these tails also affect the very-large-scale structures, making heavy clusters like "El Gordo", or large voids like the one associated with the cosmic microwave background cold spot, more probable. We compute the halo mass function and cluster abundance as a function of redshift in the presence of exponential tails. We find that quantum diffusion generically enlarges the number of heavy clusters and depletes subhalos, an effect that cannot be captured by the famed $f_{\mathrm{NL}}$ corrections. These late-universe signatures could thus be fingerprints of quantum dynamics during inflation that should be incorporated in $N$-body simulations and checked against astrophysical data., Comment: 5 pages without appendices (total 8 pages), 3 figures, matches published version in Phys. Rev. Letters

We show how importance sampling can be used to reconstruct the statistics of rare cosmological fluctuations in stochastic inflation. We have developed a publicly available package, PyFPT, that solves the first-passage time problem of generic one-dimensional Langevin processes. In the stochastic-$\delta N$ formalism, these are related to the curvature perturbation at the end of inflation. We apply this method to quadratic inflation, where the existence of semi-analytical results allows us to benchmark our approach. We find excellent agreement within the estimated statistical error, both in the drift- and diffusion-dominated regimes. The computation takes at most a few hours on a single CPU, and can reach probability values corresponding to less than one Hubble patch per observable universe at the end of inflation. With direct sampling, this would take more than the age of the universe to simulate even with the best current supercomputers. As an application, we study how the presence of large-field boundaries might affect the tail of the probability distribution. We also find that non-perturbative deviations from Gaussianity are not always of the simple exponential type., Comment: 21+13 pages, 11 figures. Github: https://github.com/Jacks0nJ/PyFPT. Version accepted for publication in JCAP

The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational wave observations by LISA to probe the universe.

Bell-inequality violations reveal the presence of quantum correlations between two particles that have interacted and then separated. Their generalisation to quantum fields is necessary to study a number of field-theoretic setups, such as cosmological density fluctuations. In this work, we show how Bell operators can be constructed for quantum fields in real space, and for Gaussian states we compute their expectation value in terms of the field power spectra. We then apply our formalism to a scalar field in de-Sitter space-time. We find that, in spite of the tremendous production of entangled particles with opposite wave momenta on large scales, Bell inequalities are not violated in real space. The reason is that, when considering measurements of a field at two distinct locations in real space, one implicitly traces over the configuration of the field at every other location, leading to a mixed bipartite system. This "effective decoherence" effect is responsible for the erasure of quantum features, and casts some doubts on our ability to reveal the quantum origin of cosmological structures. We finally discuss these results in the light of quantum discord., Comment: 22 pages without appendices (31 pages in total), 10 figures. Comments and references added. Matches version published in JCAP

In quantum information theory, quantum discord has been proposed as a tool to characterise the presence of "quantum correlations" between the subparts of a given system. Whether a system behaves quantum-mechanically or classically is believed to be impacted by the phenomenon of decoherence, which originates from the unavoidable interaction between this system and an environment. Generically, decoherence is associated with a decrease of the state purity, i.e. a transition from a pure to a mixed state. In this paper, we investigate how quantum discord is modified by this quantum-to-classical transition. This study is carried out on systems described by quadratic Hamiltonians and Gaussian states, with generalised squeezing parameters. A generic parametrisation is also introduced to describe the way the system is partitioned into two subsystems. We find that the evolution of quantum discord in presence of an environment is a competition between the growth of the squeezing amplitude and the decrease of the state purity. In phase space, this corresponds to whether the semi-minor axis of the Wigner ellipse increases or decreases, which has a clear geometrical interpretation. Finally, these considerations are applied to primordial cosmological perturbations, thus allowing us to investigate how large-scale structures in our universe, which are believed to arise from quantum fluctuations, can exhibit classical properties., Comment: 35 pages without appendices (total 57 pages), 10 figures, typo corrected in the labels of figure 3

We present a generic framework to compute the one-point statistics of cosmological perturbations, when coarse-grained at an arbitrary scale $R$, in the presence of quantum diffusion. Making use of the stochastic-$\delta N$ formalism, we show how it can be related to the statistics of the amount of expansion realised until the scale $R$ crosses out the Hubble radius. This leads us to explicit formulae for the probability density function (PDF) of the curvature perturbation, the comoving density contrast, and the compaction function. We then apply our formalism to the calculation of the mass distribution of primordial black holes produced in a single-field model containing a "quantum well" (i.e. an exactly flat region in the potential). We confirm that the PDFs feature heavy, exponential tails, with an additional cubic suppression in the case of the curvature perturbation. The large-mass end of the mass distribution is shown to be mostly driven by stochastic-contamination effects, which produce black holes more massive than those naively expected. This work bridges the final gap between the stochastic-inflation formalism and the calculation of the mass distribution of astrophysical objects such as primordial black holes, and opens up various prospects that we finally discuss., Comment: 39 pages, 6 figures; v2 matches the published version in JCAP; v3 corrects typos; v4 corrects typos

The separate-universe approach provides an effective description of cosmological perturbations at large scales, where the universe can be described by an ensemble of independent, locally homogeneous and isotropic patches. By reducing the phase space to homogeneous and isotropic degrees of freedom, it greatly simplifies the analysis of large-scale fluctuations. It is also a prerequisite for the stochastic-inflation formalism. In this work, we formulate the separate-universe approach in the Hamiltonian formalism, which allows us to analyse the full phase-space structure of the perturbations. Such a phase-space description is indeed required in dynamical regimes which do not benefit from a background attractor, as well as to investigate quantum properties of cosmological perturbations. We find that the separate-universe approach always succeeds in reproducing the same phase-space dynamics for homogeneous and isotropic degrees of freedom as the full cosmological perturbation theory, provided that the wavelength of the modes under consideration are larger than some lower bound that we derive. We also compare the separate-universe approach and cosmological perturbation theory at the level of the gauge-matching procedure, where the agreement is not always guaranteed and requires specific matching prescriptions that we present., Comment: Discussion around equation (4.12) expanded, few minor changes (main conclusions unchanged)

We compute the entanglement entropy, mutual information and quantum discord of the Cosmic Microwave Background (CMB) fluctuations in real space. To that end, we first show that measurements of these fluctuations at two distinct spatial locations can be described by a bipartite, continuous Gaussian system. This leads to explicit formulas for the mutual information and the quantum discord in terms of the Fourier-space power spectra of the curvature perturbation. We then find that quantum entanglement, that builds up in Fourier space between opposite wave momenta as an effect of quantum squeezing, is transferred to real space. In particular, both the mutual information and quantum discord, which decay as the fourth power of the distance between the two measurements in flat space time, asymptotes a constant in cosmological backgrounds. At the scales probed in the CMB however, they are highly suppressed, while they can reach order-one values at much smaller scales, where primordial black holes could have formed., Comment: companion of arXiv:2106.14575, 37 pages without appendices (total 46 pages), 7 figures, typos corrected in equations (A.8), (A.10), (A.13), (A.15)

We introduce a new method permitting the analytical determination of entanglement entropy (and related quantities) between configurations of a quantum field, which is either free or in interaction with a classical source, at two distinct spatial locations. We show how such a setup can be described by a bipartite, continuous Gaussian system. This allows us to derive explicit and exact formulas for the entanglement entropy, the mutual information and the quantum discord, solely in terms of the Fourier-space power spectra of the field. This contrasts with previous studies, which mostly rely on numerical considerations. As an illustration, we apply our formalism to massless fields in flat space, where exact expressions are derived that only involve the ratio between the size of the regions over which the field is coarse-grained, and the distance between these regions. In particular, we recover the well-known fact that mutual information decays as the fourth power of this ratio at large distances, as previously observed in numerical works. Our method leads to the first analytical derivation of this result, and to an exact formula that also applies to arbitrary distances. Finally, we determine the quantum discord and find that it identically vanishes (unless coarse-graining is performed over smeared spheres, in which case it obeys the same suppression at large distance as mutual information)., Comment: companion of arXiv:2106.15100, matches published version in Phys. Rev. D

We construct the four-mode squeezed states and study their physical properties. These states describe two linearly-coupled quantum scalar fields, which makes them physically relevant in various contexts such as cosmology. They are shown to generalise the usual two-mode squeezed states of single-field systems, with additional transfers of quanta between the fields. To build them in the Fock space, we use the symplectic structure of the phase space. For this reason, we first present a pedagogical analysis of the symplectic group $\mathrm{Sp}(4,\mathbb{R})$ and its Lie algebra, from which we construct the four-mode squeezed states and discuss their structure. We also study the reduced single-field system obtained by tracing out one of the two fields. This procedure being easier in the phase space, it motivates the use of the Wigner function which we introduce as an alternative description of the state. It allows us to discuss environmental effects in the case of linear interactions. In particular, we find that there is always a range of interaction coupling for which decoherence occurs without substantially affecting the power spectra (hence the observables) of the system., Comment: 38 pages without appendices (total 49 pages), matches published version in EPJC

The Continuous Spontaneous Localisation (CSL) theory in the cosmological context is subject to uncertainties related to the choice of the collapse operator. In this paper, we constrain its form based on generic arguments. We show that, if the collapse operator is even in the field variables, it is unable to induce the collapse of the wavefunction. Instead, if it is odd, we find that only linear operators are such that the outcomes are distributed according to Gaussian statistics, as required by measurements of the cosmic microwave background. We discuss implications of these results for previously proposed collapse operators. We conclude that the cosmological CSL collapse operator should be linear in the field variables., Comment: typo in the threading terminology corrected

We consider the effect of quantum diffusion on the dynamics of the inflaton during a period of ultra-slow-roll inflation. We extend the stochastic-$\delta\mathcal{N}$ formalism to the ultra-slow-roll regime and show how this system can be solved analytically in both the classical-drift and quantum-diffusion dominated limits. By deriving the characteristic function, we are able to construct the full probability distribution function for the primordial density field. In the diffusion-dominated limit, we recover an exponential tail for the probability distribution, as found previously in slow-roll inflation. To complement these analytical techniques, we present numerical results found both by very large numbers of simulations of the Langevin equations, and through a new, more efficient approach based on iterative Volterra integrals. We illustrate these techniques with two examples of potentials that exhibit an ultra-slow-roll phase leading to the possible production of primordial black holes., Comment: 35 pages without appendices (total 43 pages), 9 figures, matches the version published in JCAP

We compute the power spectrum of curvature perturbations in stochastic inflation. This combines the distribution of first crossing times through the end-of-inflation surface, which has been previously studied, with the distribution of the fields value at the time when a given scale crosses out the Hubble radius during inflation, which we show how to compute. This allows the stochastic-$\delta N$ formalism to make concrete contact with observations. As an application, we study how quantum diffusion at small scales (arising e.g. in models leading to primordial black holes) affects the large-scale perturbations observed in the cosmic microwave background. We find that even if those sets of scales are well separated, large effects can arise from the distortion of the classical relationship between field values and wavenumbers brought about by quantum diffusion near the end of inflation. This shows that cosmic microwave background measurements can set explicit constraints on the entire inflationary potential down to the end of inflation., Comment: References added, matches published version in JCAP

We calculate the mass distribution of Primordial Black Holes (PBHs) produced during metric preheating. After inflation, the oscillations of the inflaton at the bottom of its potential source a parametric resonant instability for small-scale scalar perturbations, that may collapse into black holes. After reviewing in a pedagogical way different techniques that have been developed in the literature to compute mass distributions of PBHs, we focus on the excursion-set approach. We derive a Volterra integral equation that is free of a singularity sometimes encountered, and apply it to the case of metric preheating. We find that if the energy density at which the instability stops, $\rho_\Gamma$, is sufficiently smaller than the one at which inflation ends, $\rho_\mathrm{end}$, namely if $\rho_\Gamma^{1/4}/\rho_\mathrm{end}^{1/4}< 10^{-5}(\rho_\mathrm{end}^{1/4}/10^{16}\mathrm{GeV})^{3/2}$, then PBHs dominate the universe content at the end of the oscillatory phase. This confirms the previous analysis of arXiv:1907.04236 . By properly accounting for the "cloud-in-cloud" mechanism, we find that the mass distribution is more suppressed at low masses than previously thought, and peaks several orders of magnitude above the Hubble mass at the end of inflation. The peak mass ranges from $10$ g to stellar masses, giving rise to different possible cosmological effects that we discuss., Comment: Matches the published version on JCAP

Ultra-light primordial black holes, with masses $m_\mathrm{PBH}<10^9\mathrm{g}$, evaporate before big-bang nucleosynthesis and can therefore not be directly constrained. They can however be so abundant that they dominate the universe content for a transient period (before reheating the universe via Hawking evaporation). If this happens, they support large cosmological fluctuations at small scales, which in turn induce the production of gravitational waves through second-order effects. Contrary to the primordial black holes, those gravitational waves survive after evaporation, and can therefore be used to constrain such scenarios. In this work, we show that for induced gravitational waves not to lead to a backreaction problem, the relative abundance of black holes at formation, denoted $ \Omega_\mathrm{PBH,f} $, should be such that $ \Omega_\mathrm{PBH,f} <10^{-4}(m_\mathrm{PBH}/10^9\mathrm{g})^{-1/4}$. In particular, scenarios where primordial black holes dominate right upon their formation time are all excluded (given that $m_\mathrm{PBH}>10\, \mathrm{g}$ for inflation to proceed at $\rho^{1/4}<10^{16}\mathrm{GeV}$). This sets the first constraints on ultra-light primordial black holes., Comment: matches the published version in JCAP

Our recent letter "Cosmic Microwave Background Constraints Cast a Shadow On Continuous Spontaneous Localization Models" [arXiv:1906.04405] has recently been criticised in [G. R. Bengochea, G. Leon, P. Pearle, and D. Sudarsky, arXiv:2008.05285, see also arXiv:2006.05313]. In this reply, we explain why the arguments presented in those articles are either incorrect or a confirmation of the robustness of our results., Comment: typo in the threading terminology corrected

During inflation, vacuum quantum fluctuations are amplified and stretched to astrophysical distances. They give rise to fluctuations in the cosmic microwave background (CMB) temperature and polarisation, and to large-scale structures in our universe. They can also trigger the formation of primordial black holes (PBHs). Such objects could provide the progenitors of the recently detected black-hole mergers, and constitute part or all of the dark matter. Their observation would give invaluable access to parts of the inflationary sector that are unconstrained by the CMB. Since PBHs require large inhomogeneities to form, they are produced in scenarios where quantum fluctuations substantially modify the dynamics of the universe. In this habilitation thesis, this "backreaction" effect is investigated using the stochastic inflation formalism, an effective theory for the long-wavelengths of quantum fields during inflation, which can be described in a classical but stochastic way once the small wavelengths have been integrated out. It describes an inflating background that gets randomly corrected by the vacuum quantum fluctuations as they get stretched to large distances. After a brief review of the stochastic inflation formalism, we explain how it can be combined with standard techniques of cosmological perturbation theory (the $\delta N$ formalism) to provide the full probability density function of curvature perturbations in the presence of non-perturbative quantum diffusion. These results are then applied to PBHs, where we show that quantum diffusion can change the expected abundance by several orders of magnitude. Finally, since inflationary models giving rise to cosmologically relevant PBHs often feature violations of slow roll, the stochastic-$\delta N$ formalism is generalised to non slow-roll dynamics. We conclude by highlighting several research directions that remain to be explored., Comment: This is the manuscript of my habilitation thesis, which was defended on the 30th of June 2020. It draws results from arXiv:1506.04732, arXiv:1604.04502, arXiv:1604.06017, arXiv:1703.00447, arXiv:1707.00537, arXiv:1806.09553, arXiv:1905.06300, arXiv:1910.01916 and arXiv:1912.05399. 169 pages, 50 figures. a very few minor typos corrected in v2

Bipartite temporal Bell inequalities are similar to the usual Bell inequalities except that, instead of changing the direction of the polariser at each measurement, one changes the time at which the measurement is being performed. By doing so, one is able to test for realism and locality, but relying on position measurements only. This is particularly useful in experimental setups where the momentum direction cannot be probed (such as in cosmology for instance). We study these bipartite temporal Bell inequalities for continuous systems placed in two-mode squeezed states, and find some regions in parameter space where they are indeed violated. We highlight the role played by the rotation angle, which is one of the three parameters characterising a two-mode squeezed state (the other two being the squeezing amplitude and the squeezing angle). In single-time measurements, it only determines the overall phase of the wavefunction and can therefore be discarded, but in multiple-time measurements, its time dynamics becomes relevant and crucially determines when bipartite temporal Bell inequalities can be violated. Our study opens up the possibility of new experimental designs for the observation of Bell inequality violations., Comment: matches published version

Contracting cosmologies are known to be flawed with a shear instability, where the contribution from the anisotropic stress to the overall energy density grows as $a^{-6}$, with $a$ the scale factor. Classically, whether or not this contribution becomes important before the bounce depends on its initial value, which can always be sufficiently fine tuned to make it irrelevant. However, vacuum quantum fluctuations inevitably provide a non-vanishing source of anisotropic stress. In this work, we compute the minimum amount of shear that is obtained if one assumes that it vanishes initially, but lets quantum fluctuations build it up. In practice, we consider a massless test scalar field, and describe its quantum fluctuations by means of the stochastic "inflation" (though here applied to a contracting phase) formalism. We find that, if the equation-of-state parameter of the contraction satisfies $w>-1/9$, regardless of when the contracting phase is initiated, the time at which the shear becomes sizeable is always when the Hubble scale approaches the Planck mass (which is also where the bounce is expected to take place). However, if $w<-1/9$, the shear backreaction becomes important much earlier, at a point that depends on the overall amount of contraction., Comment: matches the published version in EPJC

At the end of inflation, the inflaton oscillates at the bottom of its potential and these oscillations trigger a parametric instability for scalar fluctuations with wavelength $\lambda$ comprised in the instability band $(3H m)^{-1/2} <\lambda < H^{-1}$, where $H$ is the Hubble parameter and $m$ the curvature of the potential at its minimum. This "metric preheating" instability, which proceeds in the narrow resonance regime, leads to various interesting phenomena such as early structure formation, production of gravitational waves and formation of primordial black holes. In this work we study its fate in the presence of interactions with additional degrees of freedom, in the form of perturbative decay of the inflaton into a perfect fluid. Indeed, in order to ensure a complete transition from inflation to the radiation-dominated era, metric preheating must be considered together with perturbative reheating. We find that the decay of the inflaton does not alter the instability structure until the fluid dominates the universe content. As an application, we discuss the impact of the inflaton decay on the production of primordial black holes from the instability. We stress the difference between scalar field and perfect fluid fluctuations and explain why usual results concerning the formation of primordial black holes from perfect fluid inhomogeneities cannot be used, clarifying some recent statements made in the literature., Comment: 27 pages, 7 figures, matches published version

Attempts to apply quantum collapse theories to Cosmology and cosmic inflation are reviewed. These attempts are motivated by the fact that the theory of cosmological perturbations of quantum-mechanical origin suffers from the single outcome problem, which is a modern incarnation of the quantum measurement problem, and that collapse models can provide a solution to these issues. Since inflationary predictions can be very accurately tested by cosmological data, this also leads to constraints on collapse models. These constraints are derived in the case of Continuous Spontaneous Localization (CSL) and are shown to be of unprecedented efficiency., Comment: minor typo in equation (30) corrected, conclusions and results unchanged

The curvature perturbations produced during an early era of inflation are known to have quasi-Gaussian distribution functions close to their maximum, where they are well constrained by measurements of the cosmic microwave background anisotropies and of the large-scale structures. In contrast, the tails of these distributions are poorly known, although this part is the relevant one for rare, extreme objects such as primordial black holes. We show that these tails are highly non-Gaussian, and cannot be described with standard non-Gaussian expansions, that are designed to approximate the distributions close to their maximum only. Using the stochastic-$\delta N$ formalism, we develop a generic framework to compute the tails, which are found to have an exponential, rather than Gaussian, decay. These exponential tails are inevitable, and do not require any non-minimal feature as they simply result from the quantum diffusion of the inflaton field along its potential. We apply our formalism to a few relevant single-field, slow-roll inflationary potentials, where our analytical treatment is confirmed by comparison with numerical results. We discuss the implications for the expected abundance of primordial black holes in these models, and highlight that it can differ from standard results by several orders of magnitude. In particular, we find that potentials with an inflection point overproduce primordial black holes, unless slow roll is violated., Comment: 34 pages, 14 figures. Justification of the characteristic function expansion in section 2.2 amended; definition of a new scalar product in equation (2.21); figure 13 is corrected. Our results are unchanged

Canonical transformations are ubiquitous in Hamiltonian mechanics, since they not only describe the fundamental invariance of the theory under phase-space reparameterisations, but also generate the dynamics of the system. In the first part of this work we study the symplectic structure associated with linear canonical transformations. After reviewing salient mathematical properties of the symplectic group in a pedagogical way, we introduce the squeezing formalism, and show how any linear dynamics can be cast in terms of an invariant representation. In the second part, we apply these results to the case of cosmological perturbations, and focus on scalar field fluctuations during inflation. We show that different canonical variables select out different vacuum states, and that this leaves an ambiguity in observational predictions if initial conditions are set at a finite time in the past. We also discuss how the effectiveness of the quantum-to-classical transition of cosmological perturbations depends on the set of canonical variables used to describe them., Comment: 82 pages, 2 figures, matches published version

After the end of inflation, the inflaton field oscillates around a local minimum of its potential and decays into ordinary matter. These oscillations trigger a resonant instability for cosmological perturbations with wavelengths that exit the Hubble radius close to the end of inflation. In this paper, we study the formation of Primordial Black Holes (PBHs) at these enhanced scales. We find that the production mechanism can be so efficient that PBHs subsequently dominate the content of the universe and reheating proceeds from their evaporation. Observational constraints on the PBH abundance also restrict the duration of the resonant instability phase, leading to tight limits on the reheating temperature that we derive. We conclude that the production of PBHs during reheating is a generic and inevitable property of the simplest inflationary models, and does not require any fine tuning of the inflationary potential., Comment: 29 pages plus appendices (total 42 pages), 14 figures, matches published version

The Continuous Spontaneous Localisation (CSL) model solves the measurement problem of standard quantum mechanics, by coupling the mass density of a quantum system to a white-noise field. Since the mass density is not uniquely defined in general relativity, this model is ambiguous when applied to cosmology. We however show that most natural choices of the density contrast already make current measurements of the cosmic microwave background incompatible with other laboratory experiments., Comment: Minor typo in equation (52) corrected (other results unchanged)