Primordial black hole formation in the early universe: critical behaviour and self-similarity

Following on after three previous papers discussing the formation of primordial black holes during the radiative era of the early universe, we present here a further investigation of the critical nature of the process involved, aimed at making contact with some of the basic underlying ideas from the literature on critical collapse. We focus on the intermediate state, which we have found appearing in cases with perturbations close to the critical limit, and examine the connection between this and the similarity solutions which play a fundamental role in the standard picture of critical collapse. We have derived a set of self-similar equations for the null-slicing form of the metric which we are using for our numerical calculations, and have then compared the results obtained by integrating these with the ones coming from our simulations for collapse of cosmological perturbations within an expanding universe. We find that the similarity solution is asymptotically approached in a region which grows to cover both the contracting matter and part of the semi-void which forms outside it. Our main interest is in the situation relevant for primordial black hole formation in the radiative era of the early universe, where the relation between the pressure $p$ and the energy density $e$ can be reasonably approximated by an expression of the form $p = we$ with $w=1/3$. However, we have also looked at other values of $w$, both because these have been considered in previous literature and also because they can be helpful for giving further insight into situations relevant for primordial black hole formation. As in our previous work, we have started our simulations with initial supra-horizon scale perturbations of a type which could have come from inflation.Comment: 23 pages, 8 figures, new abstract, submitted to Classical and Quantum Gravity. This new version of the paper has been completely rewritten with respect the previous one, with several changes and substantial additional wor

Causal Nature and Dynamics of Trapping Horizons in Black Hole Collapse

In calculations of gravitational collapse to form black holes, trapping horizons (foliated by marginally trapped surfaces) make their first appearance either within the collapsing matter or where it joins on to a vacuum exterior. Those which then move outwards with respect to the matter have been proposed for use in defining black holes, replacing the global concept of an "event horizon" which has some serious drawbacks for practical applications. We here present results from a study of the properties of both outgoing and ingoing trapping horizons, assuming strict spherical symmetry throughout. We have investigated their causal nature (i.e. whether they are spacelike, timelike or null), making contact with the Misner-Sharp- Hernandez formalism, which has often been used for numerical calculations of spherical collapse. We follow two different approaches, one using a geometrical quantity related to expansions of null geodesic congruences, and the other using the horizon velocity measured with respect to the collapsing matter. After an introduction to these concepts, we then implement them within numerical simulations of stellar collapse, revisiting pioneering calculations from the 1960s where some features of the emergence and subsequent behaviour of trapping horizons could already be seen. Our presentation here is aimed firmly at "real world" applications of interest to astrophysicists and includes the effects of pressure, which may be important for the asymptotic behaviour of the ingoing horizon.Comment: 33 pages, 11 figure

Primordial black hole formation during the QCD phase transition: threshold, mass distribution and abundance

Primordial black hole (PBH) formation during cosmic phase transitions and annihilation periods, such as the QCD transition or the $e^+e^-$-annihilation, is known to be particularly efficient due to a softening of the equation of state. We present a detailed numerical study of PBH formation during the QCD epoch in order to derive an accurate PBH mass function. We also briefly consider PBH formation during the $e^+e^-$-annihilation epoch. Our investigation confirms that, for nearly scale-invariant spectra, PBH abundances on the QCD scale are enhanced by a factor $\sim 10^3$ compared to a purely radiation dominated Universe. For a power spectrum producing an (almost) scale-invariant PBH mass function outside of the transition, we find a peak mass of $M_{\rm pbh}\approx 1.9 M_{\odot}$ of which a fraction $f\approx 1.5\times 10^{-2}$ of the PBHs have a mass of $M_{\rm pbh} > 10 M_{\odot}$, possibly contributing to the LIGO-Virgo black hole merger detections. We point out that the physics of PBH formation during the $e^+e^-$-annihilation epoch is more complex as it is very close to the epoch of neutrino decoupling. We argue that neutrinos free-streaming out of overdense regions may actually hinder PBH formation.Comment: 20 pages, 12 Figure

Primordial black hole formation and abundance: contribution from the non-linear relation between the density and curvature perturbation

International audienceThe formation and abundance of primordial black holes (PBHs) arising from the curvature perturbation ζ is studied. The non-linear relation between ζ and the density contrast δ means that, even when ζ has an exactly Gaussian distribution, significant non-Gaussianities affecting PBH formation must be considered. Numerical simulations are used to investigate the critical value and the mass of PBHs which form, and peaks theory is used to calculate the mass fraction of the universe collapsing to form PBHs at the time of formation. A formalism to calculate the total present day PBH abundance and mass function is also derived. It is found that the abundance of PBHs is very sensitive to the non-linear effects, and that the power spectrum &calP;ζ must be a factor of &calO; (2) larger to produce the same number of PBHs as if using the linear relation between ζ and δ (where the exact value depends on the critical value for a region to collapse and form a PBH). This also means that the derived constraints on the small-scale power spectrum from constraints on the abundance of PBHs are weaker by the same factor

Primordial black hole formation in the radiative era: investigation of the critical nature of the collapse

Following on after two previous papers discussing the formation of primordial black holes in the early universe, we present here results from an in-depth investigation of the extent to which primordial black hole formation in the radiative era can be considered as an example of the critical collapse phenomenon. We focus on initial supra-horizon-scale perturbations of a type which could have come from inflation, with only a growing component and no decaying component. In order to study perturbations with amplitudes extremely close to the supposed critical limit, we have modified our previous computer code with the introduction of an adaptive mesh refinement scheme. This has allowed us to follow black hole formation from perturbations whose amplitudes are up to eight orders of magnitude closer to the threshold than we could do before. We find that scaling-law behaviour continues down to the smallest black hole masses that we are able to follow and we see no evidence of shock production such as has been reported in some previous studies and which led there to a breaking of the scaling-law behaviour at small black-hole masses. We attribute this difference to the different initial conditions used. In addition to the scaling law, we also present other features of the results which are characteristic of critical collapse in this context.Comment: 21 pages, 7 figures, the present version is updated with some changes and two new appendix. Accepted for pubblication in Classical and Quantum Gravit

Primordial black holes and their gravitational-wave signatures

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 welcom

Cosmology with the Laser Interferometer Space Antenna

254 pags:, 44 figs.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.This work is partly supported by: A.G. Leventis Foundation; Academy of Finland Grants 328958 and 345070; Alexander S. Onassis Foundation, Scholarship ID: FZO 059-1/2018-2019; Amaldi Research Center funded by the MIUR program “Dipartimento di Eccellenza” (CUP: B81I18001170001); ASI Grants No. 2016-24-H.0 and No. 2016-24-H.1-2018; Atracción de Talento Grant 2019-T1/TIC-15784; Atracción de Talento contract no. 2019-T1/TIC-13177 granted by the Comunidad de Madrid; Ayuda ‘Beatriz Galindo Senior’ by the Spanish ‘Ministerio de Universidades’, Grant BG20/00228; Basque Government Grant (IT-979-16); Belgian Francqui Foundation; Centre national d’Etudes spatiales; Ben Gurion University Kreitman Fellowship, and the Israel Academy of Sciences and Humanities (IASH) & Council for Higher Education (CHE) Excellence Fellowship Program for International Postdoctoral Researchers; Centro de Excelencia Severo Ochoa Program SEV-2016-0597; CERCA program of the Generalitat de Catalunya; Cluster of Excellence “Precision Physics, Fundamental Interactions, and Structure of Matter” (PRISMA? EXC 2118/1); Comunidad de Madrid, Contrato de Atracción de Talento 2017-T1/TIC-5520; Czech Science Foundation GAČR, Grant No. 21-16583M; Delta ITP consortium; Department of Energy under Grant No. DE-SC0008541, DE-SC0009919 and DESC0019195; Deutsche Forschungsgemeinschaft (DFG), Project ID 438947057; Deutsche Forschungsgemeinschaft under Germany’s Excellence Strategy - EXC 2121 Quantum Universe - 390833306; European Structural and Investment Funds and the Czech Ministry of Education, Youth and Sports (Project CoGraDS - CZ.02.1.01/0.0/0.0/15 003/0000437); European Union’s H2020 ERC Consolidator Grant “GRavity from Astrophysical to Microscopic Scales” (Grant No. GRAMS-815673); European Union’s H2020 ERC, Starting Grant Agreement No. DarkGRA-757480; European Union’s Horizon 2020 programme under the Marie Sklodowska-Curie Grant Agreement 860881 (ITN HIDDeN); European Union’s Horizon 2020 Research and Innovation Programme Grant No. 796961, “AxiBAU” (K.S.); European Union’s Horizon 2020 Research Council grant 724659 MassiveCosmo ERC-2016-COG; FCT through national funds (PTDC/FIS-PAR/31938/2017) and through project “BEYLA – BEYond LAmbda” with Ref. Number PTDC/FIS-AST/0054/2021; FEDER-Fundo Europeu de Desenvolvimento Regional through COMPETE2020 - Programa Operacional Competitividade e Internacionalização (POCI-01-0145- FEDER-031938) and research Grants UIDB/04434/2020 and UIDP/04434/2020; Fondation CFM pour la Recherche in France; Foundation for Education and European Culture in Greece; French ANR project MMUniverse (ANR-19-CE31-0020); FRIA Grant No.1.E.070.19F of the Belgian Fund for Research, F.R. S.-FNRS Fundação para a Ciência e a Tecnologia (FCT) through Contract No. DL 57/2016/CP1364/ CT0001; Fundação para a Ciência e a Tecnologia (FCT) through Grants UIDB/04434/2020, UIDP/04434/ 2020, PTDC/FIS-OUT/29048/2017, CERN/FIS-PAR/0037/2019 and “CosmoTests – Cosmological tests of gravity theories beyond General Relativity” CEECIND/00017/2018; Generalitat Valenciana Grant PROMETEO/2021/083; Grant No. 758792, project GEODESI; Government of Canada through the Department of Innovation, Science and Economic Development and Province of Ontario through the Ministry of Colleges and Universities; Grants-in-Aid for JSPS Overseas Research Fellow (No. 201960698); I?D Grant PID2020-118159GB-C41 of the Spanish Ministry of Science and Innovation; INFN iniziativa specifica TEONGRAV; Israel Science Foundation (Grant No. 2562/20); Japan Society for the Promotion of Science (JSPS) KAKENHI Grant Nos. 20H01899 and 20H05853; IFT Centro de Excelencia Severo Ochoa Grant SEV-2; Kavli Foundation and its founder Fred Kavli; Minerva Foundation; Ministerio de Ciencia e Innovacion Grant PID2020-113644GB-I00; NASA Grant 80NSSC19K0318; NASA Hubble Fellowship grants No. HST-HF2-51452.001-A awarded by the Space Telescope Science Institute with NASA contract NAS5-26555; Netherlands Organisation for Science and Research (NWO) Grant Number 680-91-119; new faculty seed start-up grant of the Indian Institute of Science, Bangalore, the Core Research Grant CRG/2018/002200 of the Science and Engineering; NSF Grants PHY-1820675, PHY-2006645 and PHY-2011997; Polish National Science Center Grant 2018/31/D/ ST2/02048; Polish National Agency for Academic Exchange within the Polish Returns Programme under Agreement PPN/PPO/2020/1/00013/U/00001; Pró-Reitoria de Pesquisa of Universidade Federal de Minas Gerais (UFMG) under Grant No. 28359; Ramón y Cajal Fellowship contract RYC-2017-23493; Research Project PGC2018-094773-B-C32 [MINECO-FEDER]; Research Project PGC2018-094773-B-C32 [MINECO-FEDER]; ROMFORSK Grant Project. No. 302640; Royal Society Grant URF/R1/180009 and ERC StG 949572: SHADE; Shota Rustaveli National Science Foundation (SRNSF) of Georgia (Grant FR/18-1462); Simons Foundation/SFARI 560536; SNSF Ambizione grant; SNSF professorship Grant (No. 170547); Spanish MINECO’s “Centro de Excelencia Severo Ochoa” Programme Grants SEV-2016- 0597 and PID2019-110058GB-C22; Spanish Ministry MCIU/AEI/FEDER Grant (PGC2018-094626-BC21); Spanish Ministry of Science and Innovation (PID2020-115845GB-I00/AEI/10.13039/ 501100011033); Spanish Proyectos de I?D via Grant PGC2018-096646-A-I00; STFC Consolidated Grant ST/T000732/1; STFC Consolidated Grants ST/P000762/1 and ST/T000791/1; STFC Grant ST/ S000550/1; STFC Grant ST/T000813/1; STFC Grants ST/P000762/1 and ST/T000791/1; STFC under the research Grant ST/P000258/1; Swiss National Science Foundation (SNSF), project The Non-Gaussian Universe and Cosmological Symmetries, Project Number: 200020-178787; Swiss National Science Foundation Professorship Grants No. 170547 and No. 191957; SwissMap National Center for Competence in Research; “The Dark Universe: A Synergic Multi-messenger Approach” Number 2017X7X85K under the MIUR program PRIN 2017; UK Space Agency; UKSA Flagship Project, Euclid.Peer reviewe

Formation of primordial black holes

This thesis presents results from general relativistic numerical computations of primordial black-hole formation during the radiation-dominated era of the universe. Growing-mode perturbations are specified within the linear regime and their subsequent evolution is followed as they become nonlinear. We use a spherically symmetric Lagrangian code and study both super-critical perturbations, which go on to produce black holes, and sub-critical perturbations, for which the overdensity eventually disperses into the background medium. For super-critical perturbations, we revisit the results of previous work concerning scaling-laws, noting that the threshold amplitude for a perturbation to lead to black-hole formation is substantially reduced when the initial conditions are taken to represent purely growing modes. For sub-critical cases, where an initial collapse is followed by a subsequent re-expansion, strong compressions and rarefactions are seen for perturbation amplitudes near to the threshold. We have also investigated the effect of including a significant component of vacuum energy and have calculated the resulting changes in the threshold and in the slope of the scaling law. The specification of the growing-mode perturbations in the above work is approximate and in the later part of the thesis, we introduce a more sophisticated and elegant formulation in terms of curvature perturbations. This allows a direct connection to be made with the spectrum of perturbations coming from inflation and also, using this, we find that there is no longer evidence of shock production in connection with primordial black hole formation. Introducing adaptive mesh refinement into our code, we are able to follow black hole formation nearer to the critical limit and find evidence suggesting that scaling laws may continue down to very small n1asses, in contrast with previous suggestions in the literature