Search

Explaining the origin of supermassive black holes via a primordial origin is severely challenged by the tight spectral distortion constraints on the amplitude of the primordial perturbations. Following the first calculation of how the $\mu$ constraints are modified by non-Gaussianity in a companion paper, we here make the first robust constraints on primordial black hole formation under large non-Gaussianity. Even the infinite $f_{\rm NL}$ limit is insufficiently non-Gaussian but much higher-order non-Gaussianity of the form ${\cal R}={\cal R}_{\rm G}^5$ may allow the formation of any mass primordial black hole without conflicting with distortion constraints. We caution that such extreme models face other challenges., Comment: 24 pages, 6 figures. v2: Minor typos corrected and references added to match the published version

We examine the shape of the primordial black hole mass distribution arising from a peak in the primordial power spectrum. In light of improvements to the modelling, we revisit the claim that the effects of critical collapse produce a distribution that is not described by the commonly assumed lognormal, showing that this conclusion remains valid, particularly for narrow peaks where the shape of the mass distribution is insensitive to the peak properties and critical collapse determines a minimum width. We propose some alternative models that may better describe the shape, both for the narrow peak case and for much broader peaks where the effect of the peak shape is significant. We highlight the skew-lognormal and a generalised model motivated by the physics of critical collapse as the best of these possible alternatives. These models can be used as an accurate and fast approximation to the numerically calculated mass distribution, allowing for efficient implementation in an MCMC analysis. We advocate the use of one of these two models instead of the lognormal with sufficiently accurate data, such as future LIGO-Virgo observations, or when considering strongly mass dependent constraints on the PBH abundance., Comment: 15 pages + appendices, 5 figures; v3: additional discussion of mass distribution insensitivity to narrow peak shape, matches published version

We conduct a thorough Bayesian analysis of the possibility that the black hole merger events seen in gravitational waves are primordial black hole (PBH) mergers. Using the latest merger rate models for PBH binaries drawn from a lognormal mass function we compute posterior parameter constraints and Bayesian evidences using data from the first two observing runs of LIGO-Virgo. We account for theoretical uncertainty due to possible disruption of the binary by surrounding PBHs, which can suppress the merger rate significantly. We also consider simple astrophysically motivated models and find that these are favoured decisively over the PBH scenario, quantified by the Bayesian evidence ratio. Paying careful attention to the influence of the parameter priors and the quality of the model fits, we show that the evidence ratios can be understood by comparing the predicted chirp mass distribution to that of the data. We identify the posterior predictive distribution of chirp mass as a vital tool for discriminating between models. A model in which all mergers are PBH binaries is strongly disfavoured compared with astrophysical models, in part due to the over-prediction of heavy systems having $\mathcal{M}_{{\rm chirp}} \gtrsim 40 \, M_\odot$ and positive skewness over the range of observed masses which does not match the observations. We find that the fit is not significantly improved by adding a maximum mass cut-off, a bimodal mass function, or imposing that PBH binaries form at late times. We argue that a successful PBH model must either modify the lognormal shape of the initial mass function significantly or abandon the hypothesis that all observed merging binaries are primordial. We develop and apply techniques for analysing PBH models with gravitational wave data which will be necessary for robust statistical inference as the gravitational wave source sample size increases., Comment: 29+8 pages, 18 figures. Minor edits to match version published in Physical Review D

We compare primordial black hole (PBH) constraints on the power spectrum and mass distributions using the traditional Press Schechter formalism, peaks theory, and a recently developed version of peaks theory relevant to PBHs. We show that, provided the PBH formation criteria and the power spectrum smoothing are treated consistently, the constraints only vary by $\sim$10\% between methods (a difference that will become increasingly important with better data). Our robust constraints from PBHs take into account the effects of critical collapse, the non-linear relation between $\zeta$ and $\delta$, and the shift from the PBH mass to the power spectrum peak scale. We show that these constraints are remarkably similar to the pulsar timing array (PTA) constraints impacting the black hole masses detected by the LIGO and Virgo, but that the $\mu$-distortion constraints rule out supermassive black hole (SMBH) formation and potentially even the much lighter mass range of $\sim$(1-100) $\mathrm{M}_\odot$ that LIGO/Virgo probes., Comment: 24 pages + Appendices, 8 figures; Matches published version

We consider a disformal coupling between Standard Model matter and a cubic Galileon scalar sector, assumed to be a relict of some other physics that solves the cosmological constant problem rather than a solution in its own right. This allows the energy density carried by the Galileon scalar to be sufficiently small that it evades stringent constraints from the integrated Sachs-Wolfe effect, which otherwise rules out the cubic Galileon theory. Although the model with disformal coupling does not exhibit screening, we show there is a `screening-like' phenomenon in which the energy density carried by the Galileon scalar is suppressed during matter domination when the quadratic and cubic Galileon operators are both relevant and the quadratic sector has a stable kinetic term. We obtain the explicit 3+1 form of Maxwell's equations in the presence of the disformal coupling, and the wave equations that govern electromagnetic waves. The disformal coupling is known to generate a small mass that modifies their velocity of propagation. We use the WKB approximation to study electromagnetic waves in this theory and show that, despite remarkable recent constraints from the LIGO/Virgo observatories that restrict the difference in propagation velocity between electromagnetic and gravitational radiation to roughly 1 part in $10^{15}$, the disformal coupling is too weak to be constrained by events such as GW170817 or by the dispersion of electromagnetic radiation at different wavelengths., Comment: 18 pages, 1 figure. v2: Stability analysis added. This is an author-created, un-copyedited version of an article published in JCAP. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at 10.1088/1475-7516/2021/10/085

It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years., Comment: 67 pages, 18 figures. v2: Discussion and references added. Accepted for publication in The Open Journal of Astrophysics

We have calculated the detectable merger rate of primordial black holes, as a function of the redshift, as well as the binary's mass ratio, total mass and chirp mass (observables that have not previously been explored in great detail for PBHs). We consider both the current and design sensitivity of LIGO and five different primordial black hole mass functions, as well as showing a comparison to a predicted astrophysical black hole merger rate. We show that the empirical preference for nearly equal-mass binaries in current LIGO/Virgo data can be consistent with a PBH hypothesis once observational selection effects are taken into account. However, current data do exclude some PBH mass distributions, and future data may be able to rule out the possibility that all observed BH mergers had a primordial origin., Comment: 23 pages + appendices, 17 figures; v2: Minor changes, clarification of probabilities in table 1 supersedes published version

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress., Comment: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

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., Comment: 27 pages without appendices (total 35 pages), 6 figures, 1 table, matches published version in JCAP

We reinspect the calculation for the mass fraction of primordial black holes (PBHs) which are formed from primordial perturbations, finding that performing the calculation using the comoving curvature perturbation $\mathcal{R}_{c}$ in the standard way vastly overestimates the number of PBHs, by many orders of magnitude. This is because PBHs form shortly after horizon entry, meaning modes significantly larger than the PBH are unobservable and should not affect whether a PBH forms or not - this important effect is not taken into account by smoothing the distribution in the standard fashion. We discuss alternative methods and argue that the density contrast, $\Delta$, should be used instead as super-horizon modes are damped by a factor $k^{2}$. We make a comparison between using a Press-Schechter approach and peaks theory, finding that the two are in close agreement in the region of interest. We also investigate the effect of varying the spectral index, and the running of the spectral index, on the abundance of primordial black holes., Comment: 16 pages, 7 figures. Version 2: updated to match published version and include journal reference. Version 3 and 4: minor corrections, conclusions unchanged

We compute the scale dependence of fNL for models of multi-field inflation, allowing for an arbitrary field space metric. We show that, in addition to multi-field effects and self interactions, the curved field space metric provides another source of scale dependence, which arises from the field-space Riemann curvature tensor and its derivatives. The scale dependence may be detectable within the near future if the amplitude of fNL is not too far from the current observational bounds., Comment: (v1) 9 pages, 1 figure; (v2) 10 pages, minor updates including the title, to appear in Physics Letters B

We study the effects of perturbative reheating on the evolution of the curvature perturbation \zeta, in two-field inflation models. We use numerical methods to explore the sensitivity of f_NL, n_s and r to the reheating process, and present simple qualitative arguments to explain our results. In general, if a large non-Gaussian signal exists at the start of reheating, it will remain non zero at the end of reheating. Unless all isocurvature modes have completely decayed before the start of reheating, we find that the non-linearity parameter, f_NL, can be sensitive to the reheating timescale, and that this dependence is most appreciable for `runaway' inflationary potentials that only have a minimum in one direction. For potentials with a minimum in both directions, f_NL can also be sensitive to reheating if a mild hierarchy exists between the decay rates of each field. Within the class of models studied, we find that the spectral index n_s, is fairly insensitive to large changes in the field decay rates, indicating that n_s is a more robust inflationary observable, unlike the non-linearity parameter f_NL. Our results imply that the statistics of \zeta, especially f_NL, can only be reliably used to discriminate between models of two-field inflation if the physics of reheating are properly accounted for., Comment: 25 pages, 8 figures. Comments welcome

Primordial Black Holes (PBH's) can form in the early Universe from the collapse of large density fluctuations. Tight observational limits on their abundance constrain the amplitude of the primordial fluctuations on very small scales which can not otherwise be constrained, with PBH's only forming from the extremely rare large fluctuations. The number of PBH's formed is therefore sensitive to small changes in the shape of the tail of the fluctuation distribution, which itself depends on the amount of non-Gaussianity present. We study, for the first time, how quadratic and cubic local non-Gaussianity of arbitrary size (parameterised by f_nl and g_nl respectively) affects the PBH abundance and the resulting constraints on the amplitude of the fluctuations on very small scales. Intriguingly we find that even non-linearity parameters of order unity have a significant impact on the PBH abundance. The sign of the non-Gaussianity is particularly important, with the constraint on the allowed fluctuation amplitude tightening by an order of magnitude as f_nl changes from just -0.5 to 0.5. We find that if PBH's are observed in the future, then regardless of the amplitude of the fluctuations, non-negligible negative f_nl would be ruled out. Finally we show that g_nl can have an even larger effect on the number of PBH's formed than f_nl., Comment: 9 pages, 5 figures, v2: version to appear in Phys. Rev. D with minor changes, v3: typos corrected (including factor of 1/2 in erfc prefactor), no changes to results

We outline the expected constraints on non-Gaussianity from the cosmic microwave background (CMB) with current and future experiments, focusing on both the third (f_{NL}) and fourth-order (g_{NL} and \tau_{NL}) amplitudes of the local configuration or non-Gaussianity. The experimental focus is the skewness (two-to-one) and kurtosis (two-to-two and three-to-one) power spectra from weighted maps. In adition to a measurement of \tau_{NL} and g_{NL} with WMAP 5-year data, our study provides the first forecasts for future constraints on g_{NL}. We describe how these statistics can be corrected for the mask and cut-sky through a window function, bypassing the need to compute linear terms that were introduced for the previous-generation non-Gaussianity statistics, such as the skewness estimator. We discus the ratio A_{NL} = \tau_{NL}/(6f_{NL}/5)^2 as an additional test of single-field inflationary models and discuss the physical significance of each statistic. Using these estimators with WMAP 5-Year V+W-band data out to l_{max}=600 we constrain the cubic order non-Gaussianity parameters \tau_{NL}, and g_{NL} and find -7.4 < g_{NL}/10^5 < 8.2 and -0.6 < \tau_{NL}/10^4 < 3.3 improving the previous COBE-based limit on \tau_{NL} < 10^8 nearly four orders of magnitude with WMAP., Comment: 15 pages. 14 figures

We study the generation of non-Gaussianity in models of hybrid inflation with two inflaton fields, (2-brid inflation). We analyse the region in the parameter and the initial condition space where a large non-Gaussianity may be generated during slow-roll inflation which is generally characterised by a large f_NL, tau_NL and a small g_NL. For certain parameter values we can satisfy tau_NL>>f_NL^2. The bispectrum is of the local type but may have a significant scale dependence. We show that the loop corrections to the power spectrum and bispectrum are suppressed during inflation, if one assume that the fields follow a classical background trajectory. We also include the effect of the waterfall field, which can lead to a significant change in the observables after the waterfall field is destabilised, depending on the couplings between the waterfall and inflaton fields., Comment: 16 pages, 6 figures; v2: comments and references added, typos corrected, matches published version

We study the level of primordial non-Gaussianity in slow-roll two-field inflation. Using an analytic formula for the nonlinear parameter f_nl in the case of a sum or product separable potential, we find that it is possible to generate significant non-Gaussianity even during slow-roll inflation with Gaussian perturbations at Hubble exit. In this paper we give the general conditions to obtain large non-Gaussianity and calculate the level of fine-tuning required to obtain this. We present explicit models in which the non-Gaussianity at the end of inflation can exceed the current observational bound of |f_nl|<100., Comment: 16 pages, 6 figures, 1 table, v2: typos corrected and references added, matches version accepted by JCAP

We present Feynman type diagrams for calculating the n-point function of the primordial curvature perturbation in terms of scalar field perturbations during inflation. The diagrams can be used to evaluate the corresponding terms in the n-point function at tree level or any required loop level. Rules are presented for drawing the diagrams and writing down the corresponding terms in real space and Fourier space. We show that vertices can be renormalised to automatically account for diagrams with dressed vertices. We apply these rules to calculate the primordial power spectrum up to two loops, the bispectrum including loop corrections, and the trispectrum., Comment: 17 pages, 13 figures. v2: Comments and references added, v3: Introduction expanded, subsection on evaluating loop diagrams added, minor errors corrected, references added