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 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)$.

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 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.

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.

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.

We explore the possibility that Dark Matter is the lightest hadron made of two stable color octet Dirac fermions ${\cal Q}$. The cosmological DM abundance is reproduced for $M_{\cal Q}\approx 12.5$ TeV, compatibly with direct searches (the Rayleigh cross section, suppressed by $1/M_{\cal Q}^6$, is close to present bounds), indirect searches (enhanced by ${\cal Q}{\cal Q}+\bar{\cal Q}\bar{\cal Q}\to {\cal Q}\bar{\cal Q}+{\cal Q}\bar{\cal Q}$ recombination), and with collider searches (where ${\cal Q}$ manifests as tracks, pair produced via QCD). Hybrid hadrons, made of $\cal Q$ and of SM quarks and gluons, have large QCD cross sections, and do not reach underground detectors. Their cosmological abundance is $10^5$ times smaller than DM, such that their unusual signals seem compatible with bounds. Those in the Earth and stars sank to their centers; the Earth crust and meteorites later accumulate a secondary abundance, although their present abundance depends on nuclear and geological properties that we cannot compute from first principles.

The results of the analysis of the data collected with the NEMO Phase-2 tower, deployed at 3500 m depth about 80 km off-shore Capo Passero (Italy), are presented. Cherenkov photons detected with the photomultipliers tubes were used to reconstruct the tracks of atmospheric muons. Their zenith-angle distribution was measured and the results compared with Monte Carlo simulations. An evaluation of the systematic effects due to uncertainties on environmental and detector parameters is also included. The associated depth intensity relation was evaluated and compared with previous measurements and theoretical predictions. With the present analysis, the muon depth intensity relation has been measured up to 13 km of water equivalent.