Primordial black holes may have formed in the radiative era of the early Universe from the collapse of large enough amplitude perturbations of the metric. These correspond to non linear energy density perturbations characterized by an amplitude larger than a certain threshold, measured when the perturbations reenter the cosmological horizon. The process of primordial black hole formation is studied here within spherical symmetry, using the gradient expansion approximation in the long wavelength limit, where the pressure gradients are small, and the initial perturbations are functions only of a time-independent curvature profile. In this regime it is possible to understand how the threshold for primordial black hole formation depends on the shape of the initial energy density profile, clarifying the relation between local and averaged measures of the perturbation amplitude. Although there is no universal threshold for primordial black hole formation, the averaged mass excess of the perturbation depends on the amplitude of the energy density peak, and it is possible to formulate a well-defined criterion to establish when a cosmological perturbation is able to form a black hole in terms of one of these two key quantities. This gives understanding of how the abundance of primordial black holes depends on the shape of the the inflationary power spectrum of cosmological perturbations.

Primordial black holes could have been formed in the early universe from sufficiently large cosmological perturbations re-entering the horizon when the Universe is still radiation dominated. These originate from the spectrum of curvature perturbations generated during inflation at small-scales. Because of the non-linear relation between the curvature perturbation $\zeta$ and the overdensity $\delta\rho$, the formation of the primordial black holes is affected by intrinsic non-Gaussianity even though the curvature perturbation is Gaussian. We investigate the impact of this non-Gaussianity on the critical threshold $\delta_c$ which measures the excess of mass of the perturbation, finding a relative change with respect to the value obtained using a linear relation between $\zeta$ and $\delta\rho$, of a few percent. This shows that the value of the critical threshold is rather robust against non-linearities. The same holds also when cosmologically interesting values of local primordial non-Gaussianity are added to the curvature perturbation.

The formation and abundance of primordial black holes (PBHs) arising from the curvature perturbation $\zeta$ is studied. The non-linear relation between $\zeta$ and the density contrast $\delta$ means that, even when $\zeta$ 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 $\mathcal{P}_\zeta$ must be a factor of $\sim2$ larger to produce the same number of PBHs as the linear model (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.

In the model where Primordial Black Holes (PBHs) form from large primordial curvature (C) perturbations, i.e., CPBHs, constraints on PBH abundance provide in principle constraints on the primordial curvature power spectrum. This connection however depends necessarily on the details of PBH formation mechanism. In this paper we provide, for the first time, constraints on the primordial curvature power spectrum from the latest limits on PBH abundance, taking into account all the steps from gravitational collapse in real space to PBH formation. In particular, we use results from numerical relativity simulations and peak theory to study the conditions for PBH formation for a range of perturbation shapes, including non-linearities, perturbation profile and a careful treatment of smoothing and filtering scales. We then obtain updated PBH formation conditions and translate that into primordial spectrum constraints for a wide range of shapes and abundances. These updated constraints cover a range of scales not probed by other cosmological observables. Our results show that the correct and accurate modelling of non-linearities, filtering and typical perturbation profile, is crucial for deriving meaningful cosmological implications.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

In this letter, combining peak theory and the numerical analysis of gravitational collapse in the radiation dominated era, we show that the abundance of primordial blacks holes, generated by an enhancement in the inflationary power spectrum, is extremely dependent on the shape of the peak. Given the amplitude of the power spectrum, we show that the density of primordial black holes generated from a narrow peak, is exponentially smaller than in the case of a broad peak. Specifically, for a top-hat profile of the power spectrum in Fourier space, we find that for having primordial black holes comprising all of the dark matter, one would only need a power spectrum amplitude an order of magnitude smaller than suggested previously whereas in the case of a narrow peak, one would instead need a much larger power spectrum amplitude, which in many cases would invalidate the perturbative analysis of cosmological perturbations. Finally, we show that, although critical collapse gives a broad mass spectrum, the density of primordial black holes formed is dominated by masses roughly equal to the cosmological horizon mass measured at horizon crossing.

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.

Following on after three previous papers discussing the formation of primordial black holes in the early universe during the radiation dominated era, we present here related results considering the theoretical possibility of having a fluid with a different sound speed. In practice we vary the parameter w of the equation of state p = we between 0 and 1, where e is the energy density, and p the pressure of the fluid. We consider also the possibility of having a cosmological constant interpreted as a vacuum energy component (corresponding to w = -1) revising the results presented in our first paper of this series, having now a more refined numerical code available. To perform our simulations we have fine tuned our previous computer code allowing us to observe that the nature of the critical collapse is conserved for all the different values of w considered, in agreement with previous studies of ideal collapses done in a different context from PBHs. As already done in our previous works, we start the simulations with initial supra-horizon scale perturbations of a type which could have come from inflation, characterized by a pure growing mode component.

This is a write-up of a talk given at the Opava RAGtime meeting in 2011, but it has been updated to include some subsequent related developments. The talk focused on discussion of some aspects of black hole and cosmological horizons under rather general circumstances, and on two different topics related to formation of cosmological structures at different epochs of the universe: virialization of cold dark matter during standard structure formation in the matter-dominated era, and primordial black hole formation during the radiative era.

Models including an energy transfer from CDM to DE are widely considered in the literature, namely to allow DE a significant high-z density. Strongly Coupled cosmologies assume a much larger coupling between DE and CDM, together with the presence of an uncoupled warm DM component, as the role of CDM is mostly restricted to radiative eras. This allows us to preserve small scale fluctuations even if the warm particle, possibly a sterile neutrino, is quite light, O(100 eV). Linear theory and numerical simulations show that these cosmologies agree with LCDM on supergalactic scales; e.g., CMB spectra are substantially identical. Simultaneously, simulations show that they significantly ease problems related to the properties of MW satellites and cores in dwarfs. SC cosmologies also open new perspectives on early black hole formation, and possibly lead towards unificating DE and inflationary scalar fields.

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.

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.

Results are presented 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 confirm the results of previous work concerning scaling-laws but note 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.