Primordial black hole mass function with mass gap

We analytically compute the power spectrum of primordial curvature perturbations in Type III hilltop inflation models under the slow-roll approximation. The model parameters are constrained using current Cosmic Microwave Background (CMB) data. The curvature perturbations that exit the horizon at small scales show sufficiently large amplitudes to produce primordial black holes (PBHs). We then consider the quantum one-loop corrections in these models from both the self-interaction of the inflaton and its interaction with the waterfall field. We show the loop corrections in both cases for 60 e-folds of inflation are negligible, ensuring the tree-level results are reliable within the chosen parameter regime.

The primordial black hole (PBH) is an effective candidate for dark matter. In this work, the PBH abundance $f$ is calculated in peak theory, with one or two perturbations in the inflaton potential. We construct an antisymmetric perturbation that can create a perfect plateau in the inflaton potential, leading inflation to the ultraslow-roll stage. During this stage, the power spectrum of primordial curvature perturbation is remarkably enhanced on small scales, generating abundant PBHs. The PBH abundance $f\ensuremath{\sim}0.1$ can be achieved in one or two typical mass windows at ${10}^{\ensuremath{-}17}\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$, ${10}^{\ensuremath{-}13}\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$, and $30\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$, without spoiling the nearly scale-invariant power spectrum on large scales. For comparison, $f$ is calculated in two approximate methods of peak theory (with different spectral moments) and also in the Press-Schechter theory. It is found that the Press-Schechter theory systematically underestimates $f$ by 2 or 3 orders of magnitude compared with peak theory.

We propose an inflation scenario with three independent stages of cold, warm and thermal inflation, respectively, driven by different scalar fields, motivated by the large number of such fields predicted in most extensions of the Standard Model. We show, in particular, that the intermediate period of warm inflation naturally leads to large density fluctuations on small scales, which can lead to primordial black hole formation in the mass window where they may account for all dark matter. This type of scenario yields a distinctive primordial black hole mass function with a mass gap, with the final period of thermal inflation diluting the abundance of very light black holes.

We consider the possibility that the majority of dark matter in our Universe consists of black holes of primordial origin. We determine the conditions under which such black holes may have originated from a single-field model of inflation characterized by a quartic polynomial potential. We also explore the effect of higher-dimensional operators. The large power spectrum of curvature perturbations that is needed for a large black hole abundance sources sizable second order tensor perturbations. The resulting stochastic background of primordial gravitational waves could be detected by the future space-based observatories LISA and DECIGO or—as long as we give up on the dark matter connection—by the ground-based Advanced LIGO-Virgo detector network.

Well-known causality arguments show that events occurring during or at the end of inflation, associated with reheating or preheating, could contribute a blue component to the spectrum of primordial curvature perturbations, with the dependence k3. We explore the possibility that they could be observably large in CMB, LSS, and Lyman-α data. We find that a k3 component with a cutoff at some maximum k can modestly improve the fits (Δχ2 = 2.0, 5.4) of the low multipoles (ℓ ~ 10–50) or the second peak (ℓ ~ 540) of the CMB angular spectrum when the three-year WMAP data are used. Moreover, the results from WMAP are consistent with the CBI, ACBAR, 2dFGRS, and SDSS data when they are included in the analysis. Including the SDSS galaxy clustering power spectrum, we find weak positive evidence for the k3 component at the level of Δχ2′ = 2.4, with the caveat that the nonlinear evolution of the power spectrum may not be properly treated in the presence of the k3 distortion. To investigate the high-k regime, we use the Lyman-α forest data (LUQAS, Croft et al., and SDSS Lyman-α); here we find evidence at the level Δχ2′ = 3.8. Considering that there are two additional free parameters in the model, the above results do not give a strong evidence for features; however, they show that surprisingly large bumps are not ruled out. We give constraints on the ratio between the k3 component and the nearly scale-invariant component, r3 < 1.5, over the range of wave numbers 2.3 × 10-3 Mpc -1 < k < 8.2 Mpc -1. We also discuss theoretical models which could lead to the k3 effect, including ordinary hybrid inflation and double D-term inflation models. We show that the well-motivated k3 component is also a good representative of the generic spikelike feature in the primordial perturbation power spectrum.

Constraints on inflation from cosmic microwave background and Lyman-α forest

Stellar evolution theories predict a gap in the black hole birth mass spectrum as the result of pair instability processes in the cores of massive stars. This gap, however, is not seen in the binary black hole masses inferred from gravitational wave data. One explanation is that black holes form dynamically in dense star clusters where smaller black holes merge to form more massive black holes, populating the mass gap. We show that this model predicts a distribution of the effective and precessing spin parameters, χ_{eff} and χ_{p}, within the mass gap that is insensitive to assumptions about black hole natal spins and other astrophysical parameters. We analyze the distribution of χ_{eff} as a function of primary mass for the black hole binaries in the third gravitational wave transient catalog. We infer the presence of a high mass and isotropically spinning population of black holes that is consistent with hierarchical formation in dense star clusters and a pair-instability mass gap with a lower edge at 44_{-4}^{+6}M_{⊙}. We compute a Bayes factor B>10^{4} relative to models that do not allow for a high mass population with a distinct χ_{eff} distribution. Upcoming data will enable us to tightly constrain the hierarchical formation hypothesis and refine our understanding of binary black hole formation.

Recent cosmological observations hint at a deviation from the simple power-law form of the primordial spectrum of curvature perturbations. In this paper we show that in the presence of a tensor component, a turn-over in the initial spectrum is preferred by current observations, and hence non-power-law models ought to be considered. For instance, for a power-law parameterisation with both a tensor component and running parameter, current data show a preference for a negative running at more than 2.5σ C.L. As a consequence of this deviation from a power-law, constraints on the tensor-to-scalar ratio r are slightly broader. We also present constraints on the inflationary parameters for a model-independent reconstruction and the Lasenby & Doran (LD) model. In particular, the constraints on the tensor-to-scalar ratio from the LD model are: rLD = 0.11±0.024. In addition to current data, we show expected constraints from Planck-like and CMB-Pol sensitivity experiments by using Markov-Chain-Monte-Carlo sampling chains. For all the models, we have included the Bayesian Evidence to perform a model selection analysis. The Bayes factor, using current observations, shows a strong preference for the LD model over the standard power-law parameterisation, and provides an insight into the accuracy of differentiating models through future surveys.

The properties of the recently discovered Higgs boson together with the absence of new physics at collider experiments allows us to speculate about consistently extending the Standard Model of particle physics all the way up to the Planck scale. In this context, the Standard Model Higgs non-minimally coupled to gravity could be responsible for the symmetry properties of the Universe at large scales and for the generation of the primordial spectrum of curvature perturbations seeding structure formation. We overview the minimalistic Higgs inflation scenario, its predictions, open issues and extensions and discuss its interplay with the possible metastability of the Standard Model vacuum.

We investigate the production of dark matter from metric perturbations produced during inflation or in standard cosmology. Perturbations break Weyl flatness of the Friedmann-Lemaître-Robertson-Walker metric, allowing conformally coupled fields-in particular, fermions studied here-to be produced even when their mass is negligible. Particle production can be computed by studying the Bogoliubov transformation induced by the stochastic background. For perturbations generated during inflation, we present a closed formula for the resulting abundance of particles that depends solely on the power spectrum of curvature perturbations at the end of inflation. This production mechanism can be dominant especially if the amplitude of curvature perturbations is enhanced for modes that exit the horizon toward the end of inflation. In the simplest scenario, the critical dark matter abundance is reproduced for M≳10^{6} GeV.

We report on the population properties of compact binary mergers inferred from gravitational-wave observations of these systems during the first three LIGO-Virgo observing runs. The Gravitational-Wave Transient Catalog 3 (GWTC-3) contains signals consistent with three classes of binary mergers: binary black hole, binary neutron star, and neutron star–black hole mergers. We infer the binary neutron star merger rate to be between 10 and 1700 Gpc−3 yr−1 and the neutron star–black hole merger rate to be between 7.8 and 140 Gpc−3 yr−1, assuming a constant rate density in the comoving frame and taking the union of 90% credible intervals for methods used in this work. We infer the binary black hole merger rate, allowing for evolution with redshift, to be between 17.9 and 44 Gpc−3 yr−1 at a fiducial redshift (z=0.2). The rate of binary black hole mergers is observed to increase with redshift at a rate proportional to (1+z)κ with κ=2.9−1.8+1.7 for z≲1. Using both binary neutron star and neutron star–black hole binaries, we obtain a broad, relatively flat neutron star mass distribution extending from 1.2−0.2+0.1 to 2.0−0.3+0.3M⊙. We confidently determine that the merger rate as a function of mass sharply declines after the expected maximum neutron star mass, but cannot yet confirm or rule out the existence of a lower mass gap between neutron stars and black holes. We also find the binary black hole mass distribution has localized over- and underdensities relative to a power-law distribution, with peaks emerging at chirp masses of 8.3−0.5+0.3 and 27.9−1.8+1.9M⊙. While we continue to find that the mass distribution of a binary’s more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above approximately 60M⊙, which would indicate the presence of a upper mass gap. Observed black hole spins are small, with half of spin magnitudes below χi≈0.25. While the majority of spins are preferentially aligned with the orbital angular momentum, we infer evidence of antialigned spins among the binary population. We observe an increase in spin magnitude for systems with more unequal-mass ratio. We also observe evidence of misalignment of spins relative to the orbital angular momentum. Published by the American Physical Society 2023

Cosmological parameter uncertainties are often stated assuming a particular model, neglecting the model uncertainty, even when Bayesian model selection is unable to identify a conclusive best model. Bayesian model averaging is a method for assessing parameter uncertainties in situations where there is also uncertainty in the underlying model. We apply model averaging to the estimation of the parameters associated with the primordial power spectra of curvature and tensor perturbations. We use CosmoNest and MultiNest to compute the model evidences and posteriors, using cosmic microwave data from WMAP, ACBAR, BOOMERanG, and CBI, plus large-scale structure data from the SDSS DR7. We find that the model-averaged 95% credible interval for the spectral index using all of the data is $0.940<{n}_{\mathrm{s}}<1.000$, where ${n}_{\mathrm{s}}$ is specified at a pivot scale $0.015\text{ }\text{ }{\mathrm{Mpc}}^{\ensuremath{-}1}$. For the tensors model averaging can tighten the credible upper limit, depending on prior assumptions.

Effects of features in the inflaton potential on the spectrum and bispectrum of primordial curvature perturbations

We present a modification of the Press–Schechter (PS) formalism to derive general mass functions for primordial black holes (PBHs), considering their formation as being associated with the amplitude of linear energy density fluctuations. To accommodate a wide range of physical relations between the linear and non-linear conditions for collapse, we introduce an additional parameter to the PS mechanism, and that the collapse occurs at either a given cosmic time, or as fluctuations enter the horizon. We study the case where fluctuations obey Gaussian statistics and follow a primordial power spectrum of broken power-law form with a blue spectral index for small scales. We use the observed abundance of supermassive black holes (SMBH) to constrain the extended mass functions taking into account dynamical friction. We further constrain the modified PS by developing a method for converting existing constraints on the PBH mass fraction, derived assuming monochromatic mass distributions for PBHs, into constraints applicable for extended PBH mass functions. We find that when considering well-established monochromatic constraints, there are regions in parameter space where all the dark matter can be made of PBHs. Of special interest is the region for the characteristic mass of the distribution ${\sim}10^2\, \mathrm{M}_\odot$, for a wide range of blue spectral indices in the scenario where PBHs form as they enter the horizon, where the linear threshold for collapse is of the order of the typical overdensities, as this is close to the black hole masses detected by LIGO, which are difficult to explain by stellar collapse.

Pair-instability (PI) is expected to open a gap in the mass spectrum of black holes (BHs) between ≈40–65 and ≈120 M⊙. The existence of the mass gap is currently being challenged by the detection of GW190521, with a primary component mass of $85^{+21}_{-14}$ M⊙. Here, we investigate the main uncertainties on the PI mass gap: the 12C(α, γ)16O reaction rate and the H-rich envelope collapse. With the standard 12C(α, γ)16O rate, the lower edge of the mass gap can be 70 M⊙ if we allow for the collapse of the residual H-rich envelope at metallicity Z ≤ 0.0003. Adopting the uncertainties given by the starlib database, for models computed with the 12C(α, γ)16O rate $-1\, \sigma$, we find that the PI mass gap ranges between ≈80 and ≈150 M⊙. Stars with MZAMS > 110 M⊙ may experience a deep dredge-up episode during the core helium-burning phase, that extracts matter from the core enriching the envelope. As a consequence of the He-core mass reduction, a star with MZAMS = 160 M⊙ may avoid the PI and produce a BH of 150 M⊙. In the $-2\, {}\sigma {}$ case, the PI mass gap ranges from 92 to 110 M⊙. Finally, in models computed with 12C(α, γ)16O $-3\, {}\sigma {}$, the mass gap is completely removed by the dredge-up effect. The onset of this dredge-up is particularly sensitive to the assumed model for convection and mixing. The combined effect of H-rich envelope collapse and low 12C(α, γ)16O rate can lead to the formation of BHs with masses consistent with the primary component of GW190521.