The LIGO and Virgo Interferometers have so far provided 11 gravitational-wave (GW) observations of black-hole binaries. Similar detections are bound to become very frequent in the near future. With the current and upcoming wealth of data, it is possible to confront specific formation models with observations. We investigate here whether current data are compatible with the hypothesis that LIGO/Virgo black holes are of primordial origin. We compute in detail the mass and spin distributions of primordial black holes (PBHs), their merger rates, the stochastic background of unresolved coalescences, and confront them with current data from the first two observational runs, also including the recently discovered GW190412. We compute the best-fit values for the parameters of the PBH mass distribution at formation that are compatible with current GW data. In all cases, the maximum fraction of PBHs in dark matter is constrained by these observations to be $f_{\text{PBH}}\approx {\rm few}\times 10^{-3}$. We discuss the predictions of the PBH scenario that can be directly tested as new data become available. In the most likely formation scenarios where PBHs are born with negligible spin, the fact that at least one of the components of GW190412 is moderately spinning is incompatible with a primordial origin for this event, unless accretion or hierarchical mergers are significant. In the absence of accretion, current non-GW constraints already exclude that LIGO/Virgo events are all of primordial origin, whereas in the presence of accretion the GW bounds on the PBH abundance are the most stringent ones in the relevant mass range. A strong phase of accretion during the cosmic history would favour mass ratios close to unity, and a redshift-dependent correlation between high masses, high spins and nearly-equal mass binaries, with the secondary component spinning faster than the primary.

We consider the constraints on the fraction of dark matter in the universe in the form of primordial black holes taking into account the crucial role of accretion which may change both their mass and mass function. We show that accretion may drastically weaken the constraints at the present epoch for primordial black holes with masses larger than a few solar masses.

Primordial black holes in the mass range of ground-based gravitational-wave detectors can comprise a significant fraction of the dark matter. Mass and spin measurements from coalescences can be used to distinguish between an astrophysical or a primordial origin of the binary black holes. In standard scenarios the spin of primordial black holes is very small at formation. However, the mass and spin can evolve through the cosmic history due to accretion. We show that the mass and spin of primordial black holes are correlated in a redshift-dependent fashion, in particular primordial black holes with masses below ${\cal O}(30)M_\odot$ are likely non-spinning at any redshift, whereas heavier black holes can be nearly extremal up to redshift $z\sim10$. The dependence of the mass and spin distributions on the redshift can be probed with future detectors such as the Einstein Telescope. The mass and spin evolution affect the gravitational waveform parameters, in particular the distribution of the final mass and spin of the merger remnant, and that of the effective spin of the binary. We argue that, compared to the astrophysical-formation scenario, a primordial origin of black hole binaries might better explain the spin distribution of merger events detected by LIGO-Virgo, in which the effective spin parameter of the binary is compatible to zero except possibly for few high-mass events. Upcoming results from LIGO-Virgo third observation run might reinforce or weaken these predictions.

An observable stochastic background of gravitational waves is generated whenever primordial black holes are created in the early universe thanks to a small-scale enhancement of the curvature perturbation. We calculate the anisotropies and non-Gaussianity of such stochastic gravitational waves background which receive two contributions, the first at formation time and the second due to propagation effects. The former contribution can be generated if the distribution of the curvature perturbation is characterized by a local and scale-invariant shape of non-Gaussianity. Under such an assumption, we conclude that a sizeable magnitude of anisotropy and non-Gaussianity in the gravitational waves would suggest that primordial black holes may not comply the totality of the dark matter.

We elaborate on the mass function of primordial black holes in the case in which the power spectrum of the curvature perturbation is broad. For the case of a broad and flat spectrum, we argue that such a mass function is peaked at the smallest primordial black mass which can be formed and possesses a tail decaying like $M^{-3/2}$, where $M$ is the mass of the primordial black hole.

We study the spin of primordial black holes produced by the collapse of large inhomogeneities in the early universe. Since such primordial black holes originate from peaks, that is, from maxima of the local overdensity, we resort to peak theory to obtain the probability distribution of the spin at formation. We show that the spin is a first-order effect in perturbation theory: it results from the action of first-order tidal gravitational fields generating first-order torques upon horizon-crossing, and from the asphericity of the collapsing object. Assuming an ellipsoidal shape, the typical value of the dimensionless parameter $a_{\rm s}=S/G_N M^2$, where $S$ is the spin and $M$ is the mass of the primordial black hole, is about $ (\Omega_{\rm dm}/\pi) \sigma_\delta\sqrt{1-\gamma^2}$. Here, $\sigma^2_\delta$ is the variance of the overdensity at horizon crossing, $\Omega_{\rm dm}$ measures the current abundance of the dark matter and the parameter $\gamma$ is a measure of the width of the power spectrum giving rise to primordial black holes. One has $\gamma=1$ for monochromatic spectra. For these narrow spectra, the suppression arises because the velocity shear, which is strongly correlated with the inertia tensor, tends to align with the principal axis frame of the collapsing object. Typical values of $a_{\rm s}$ are at the percent level.

We discuss a novel mechanism to generate gravitational waves in the early universe. A standard way to produce primordial black holes is to enhance at small-scales the overdensity perturbations generated during inflation. The latter, upon horizon re-entry, collapse into black holes. They must be sizeable enough and are therefore associated to rare peaks. There are however less sizeable and much less rare overdensity peaks which do not end up forming primordial black holes and have a non-spherical shape. Upon collapse, they possess a time-dependent non-vanishing mass quadrupole which gives rise to the generation of gravitational waves. By their nature, such gravitational waves are complementary to those sourced at second-order by the very same scalar perturbations responsible for the formation of the primordial black holes. Their amplitude is nevertheless typically about two orders of magnitude smaller and therefore hardly measurable.

We study the formation of primordial black holes when they are generated by the collapse of large overdensities in the early universe. Since the density contrast is related to the comoving curvature perturbation by a nonlinear relation, the overdensity statistics is unavoidably non-Gaussian. We show that the abundance of primordial black holes at formation may not be captured by a perturbative approach which retains the first few cumulants of the non-Gaussian probability distribution. We provide two techniques to calculate the non-Gaussian abundance of primordial black holes at formation, one based on peak theory and the other on threshold statistics. Our results show that the unavoidable non-Gaussian nature of the inhomogeneities in the energy density makes it harder to generate PBHs. We provide simple (semi-)analytical expressions to calculate the non-Gaussian abundances of the primordial black holes and show that for both narrow and broad power spectra the gaussian case from threshold statistics is reproduced by increasing the amplitude of the power spectrum by a factor ${\cal O}(2\div 3)$.

The detection of primordial gravitational waves would not only have extraordinary implications for our understanding of early cosmology, but would also give non-trivial constraints on Standard Model parameters, under the assumption that no new physics enters below the Higgs instability scale. We study the resulting bounds on the top quark mass and the strong coupling constant, discussing their theoretical uncertainties and their robustness against changes in other parameters.

Tensor non-Gaussianities are a key ingredient to test the symmetries and the presence of higher spin fields during the inflationary epoch. Indeed, the shape of the three point correlator of the graviton is totally fixed by the symmetries of the de Sitter stage and, in the case of parity conservation, gets contributions only from the ordinary gravity action plus a higher derivative term called the (Weyl)$^3$ action. We discuss current and future bounds on the three point tensor contribution from the (Weyl)$^3$ term using cosmic microwave background (CMB) bispectra. Our results indicate that forthcoming experiments, such as LiteBIRD, CMB-S4 and CORE, will detect the presence of the (Weyl)$^3$ term if $M_p^4 L^4 \sim 10^{17} r^{-4}$, where $L$ parametrizes the strength of the (Weyl)$^3$ term and $r$ is the tensor-to-scalar ratio, which corresponds to $L\gtrsim 3.2 \times 10^5 M_p^{-1}$, while the current upper limit is $M_p^4 L^4 = (1.1 \pm 4.0) \times 10^{19} r^{-4}$ (68%CL).

The idea that primordial black holes (PBHs) can comprise most of the dark matter of the universe has recently reacquired a lot of momentum. Observational constraints, however, rule out this possibility for most of the PBH masses, with a notable exception around $10^{-12} M_\odot$. These light PBHs may be originated when a sizeable comoving curvature perturbation generated during inflation re-enters the horizon during the radiation phase. During such a stage, it is unavoidable that gravitational waves (GWs) are generated. Since their source is quadratic in the curvature perturbations, these GWs are generated fully non-Gaussian. Their frequency today is about the mHz, which is exactly the range where the LISA mission has the maximum of its sensitivity. This is certainly an impressive coincidence. We show that this scenario of PBHs as dark matter can be tested by LISA by measuring the GW two-point correlator. On the other hand, we show that the short observation time (as compared to the age of the universe) and propagation effects of the GWs across the perturbed universe from the production point to the LISA detector suppress the bispectrum to an unobservable level. This suppression is completely general and not specific to our model.

A common mechanism to form primordial black holes in the early universe is by enhancing at small-scales the scalar perturbations generated during inflation. If these fluctuations have a large enough amplitude, they may collapse into primordial black holes upon horizon re-entry. Such primordial black holes may comprise the totality of the dark matter. We offer some considerations about the formation and clustering of primordial black holes when the scalar perturbations are characterised by a broad spectrum. Using the excursion set method, as well as the supreme statistics, we show that the cloud-in-cloud phenomenon, for which small mass primordial black holes may be absorbed by bigger mass ones, is basically absent. This is due to the fact that the formation of a primordial black hole is an extremely rare event. We also show that, from the point of view of mass distribution, broad and narrow spectra give similar results in the sense that the mass distribution is tilted towards a single mass. Furthermore, we argue that primordial black holes from broad spectra are not clustered at formation, their distribution is dominantly Poissonian.

There has recently been renewed interest in the possibility that the dark matter in the universe consists of primordial black holes (PBHs). Current observational constraints leave only a few PBH mass ranges for this possibility. One of them is around $10^{-12} M_\odot$. If PBHs with this mass are formed due to an enhanced scalar-perturbation amplitude, their formation is inevitably accompanied by the generation of gravitational waves (GWs) with frequency peaked in the mHz range, precisely around the maximum sensitivity of the LISA mission. We show that, if these primordial black holes are the dark matter, LISA will be able to detect the associated GW power spectrum. Although the GW source signal is intrinsically non-Gaussian, the signal measured by LISA is a sum of the signal from a large number of independent sources suppressing the non-Gaussianity at detection to an unobservable level. We also discuss the effect of the GW propagation in the perturbed universe. PBH dark matter generically leads to a detectable, purely isotropic, Gaussian and unpolarised GW signal, a prediction that is testable with LISA.

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.

We explore non-Gaussian features of a massless spin-two field in the Vasiliev theory of higher-spin gravity. The theory contains an infinite tower of interacting gauge fields with increasing spin, and admits four-dimensional asymptotically de Sitter configurations. Using a recent proposal for calculating late-time quantum correlations in Vasiliev theory, we provide an exact formula for the tensor non-Gaussianities of the massless spin-two graviton field. By general symmetry considerations, we relate our result to that produced by a tree-level calculation in a gravitational theory containing an Einstein term and a term cubic in the Weyl tensor. The relative coefficient between the two terms is calculated explicitly, exhibiting a significant contribution from the Weyl cubed term. We discuss potential cosmological implications of our results.

The final ringdown phase in a coalescence process is a valuable laboratory to test General Relativity and potentially constrain additional degrees of freedom in the gravitational sector. We introduce here an effective description for perturbations around spherically symmetric spacetimes in the context of scalar-tensor theories, which we apply to study quasi-normal modes for black holes with scalar hair. We derive the equations of motion governing the dynamics of both the polar and the axial modes in terms of the coefficients of the effective theory. Assuming the deviation of the background from Schwarzschild is small, we use the WKB method to introduce the notion of "light ring expansion". This approximation is analogous to the slow-roll expansion used for inflation, and it allows us to express the quasinormal mode spectrum in terms of a small number of parameters. This work is a first step in describing, in a model independent way, how the scalar hair can affect the ringdown stage and leave signatures on the emitted gravitational wave signal. Potential signatures include the shifting of the quasi-normal spectrum, the breaking of isospectrality between polar and axial modes, and the existence of scalar radiation.

The presence of massive particles with spin during inflation induces distinct signatures on correlation functions of primordial curvature fluctuations. In particular, the bispectrum of primordial perturbations obtains an angular dependence determined by the spin of the particle, which can be used to set constraints on the presence of such particles. If these particles are long-lived on super-Hubble scales, as is the case for example for partially massless particles, their imprint on correlation functions of curvature perturbations would be unsuppressed. In this paper, we make a forecast for how well such angular dependence can be constrained by the upcoming EUCLID spectroscopic survey via the measurement of galaxy bispectrum.

Primordial black holes as dark matter may be generated in single-field models of inflation thanks to the enhancement at small scales of the comoving curvature perturbation. This mechanism requires leaving the slow-roll phase to enter a non-attractor phase during which the inflaton travels across a plateau and its velocity drops down exponentially. We argue that quantum diffusion has a significant impact on the primordial black hole mass fraction making the classical standard prediction not trustable.

Primordial black holes may owe their origin to the small-scale enhancement of the comoving curvature perturbation generated during inflation. Their mass fraction at formation is markedly sensitive to possible non-Gaussianities in such large, but rare fluctuations. We discuss a path-integral formulation which provides the exact mass fraction of primordial black holes at formation in the presence of non-Gaussianity. Through a couple of classes of models, one based on single-field inflation and the other on spectator fields, we show that restricting to a Gaussian statistics may lead to severe inaccuracies in the estimate of the mass fraction as well as on the clustering properties of the primordial black holes.

Inflation may provide a suitable collider to probe physics at very high energies. In this paper we investigate the impact on the CMB bispectrum of higher spin fields which are long-lived on super-Hubble scales, e.g. partially massless higher spin fields. We show that distinctive statistical anisotropic signals on the CMB three-point correlator are induced and we investigate their detectability.

If there exist higher-spin particles during inflation which are light compared to the Hubble rate, they may leave distinct statistical anisotropic imprints on the correlators involving scalar and graviton fluctuations. We characterise such signatures using the dS/CFT$_3$ correspondence and the operator product expansion techniques. In particular, we obtain generic results for the case of partially massless higher-spin states.