arxiv: 2604.27496 · v1 · submitted 2026-04-30 · 🌌 astro-ph.CO · hep-ph· hep-th

Recognition: unknown

Primordial black hole dark matter from axion inflation

Gabriele Franciolini , Nadir Ijaz , Marco Peloso

Authors on Pith no claims yet

Pith reviewed 2026-05-07 09:01 UTC · model grok-4.3

classification 🌌 astro-ph.CO hep-phhep-th

keywords primordial black holesdark matteraxion inflationU(1) gauge fieldstochastic gravitational wavesLISAbackreaction

0 comments

The pith

Axion inflation with a coupled gauge field can produce primordial black holes that account for all dark matter in the asteroidal mass range.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The authors study how a U(1) gauge field coupled to the inflaton generates amplified perturbations during axion inflation, leading to the formation of primordial black holes. They improve prior calculations by staying within the homogeneous backreaction regime, using numerically computed gauge mode functions, and applying modern PBH abundance methods that include uncertainties in the density perturbation statistics. The central result is that these black holes can make up the full dark matter density for asteroid-scale masses, and this remains true even when the inflaton's gradient energy is only 0.01 to 0.1 percent of its kinetic energy. The same process necessarily produces a stochastic gravitational wave background whose amplitude is large enough for detection at LISA, which would also help distinguish the underlying statistics of the perturbations.

Core claim

PBHs can account for all of the dark matter in the asteroidal mass range, even when the inflaton gradient energy density is highly subdominant (10^{-4}--10^{-3} of the kinetic energy), supporting the validity of the backreaction scheme. This mechanism also unavoidably generates a stochastic gravitational wave background with an amplitude that will be measured at LISA and that will allow to indirectly discriminate between different statistics of δρ.

What carries the argument

The homogeneous backreaction regime applied to a U(1) gauge field with pseudo-scalar coupling to the inflaton, using numerically solved gauge mode functions together with updated PBH abundance calculations that fold in uncertainty from the statistics of δρ.

If this is right

Primordial black holes in the asteroidal mass window can constitute the entire dark matter density without needing additional fields.
The same inflationary dynamics generate a stochastic gravitational wave background whose amplitude and shape are fixed by the required PBH abundance.
Observations at LISA can indirectly test which statistics best describe the density perturbations produced by the gauge field.
The perturbative backreaction treatment remains consistent even when the inflaton gradient energy is orders of magnitude below the kinetic energy.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If the predicted gravitational wave signal is seen, it would tighten constraints on the allowed range of the gauge-field coupling strength during inflation.
Non-detection of the background at the expected level would force either a revision of the homogeneous backreaction assumption or a change in the modeling of perturbation statistics.
The same gauge-field amplification mechanism could be examined in other inflationary models that include vector fields to see whether asteroidal-mass PBH dark matter arises more generally.

Load-bearing premise

The homogeneous backreaction regime stays valid and the assumed statistics for the density perturbations correctly capture the uncertainties in calculating the black hole abundance.

What would settle it

A LISA measurement of the stochastic gravitational wave background whose amplitude or spectrum fails to match the value required to produce the observed dark matter density through this PBH channel.

Figures

Figures reproduced from arXiv: 2604.27496 by Gabriele Franciolini, Marco Peloso, Nadir Ijaz.

**Figure 1.** Figure 1: Potential of the model, Eq. (2.7), with parameters given in Eq. (2.12). Vertical lines denote the inflaton values at which different branches of the potential are connected. The inset shows the closeup of the initial branches, with the outmost left (red) vertical line representing the start of inflation. The first branch (barely visible in the main figure) is consistent with the CMB observations. • The thi… view at source ↗

**Figure 2.** Figure 2: Evolution of the inflaton (top-left panel), of the parameter ξ (top-right panel), of the Hubble rate (bottom-left panel), and of the equation of state (bottom-right panel), for the inflaton potential shown in view at source ↗

**Figure 3.** Figure 3: Vacuum (black dashed line) and sourced (red solid line) contribution to the power spectrum of the curvature perturbation ζ. The peak of the sourced contribution saturated the limit on PBH abundance (shown with a blue dashed segment) derived in the next section. The (blue) band at small wavenumbers indicate the CMB constraint on the power [1]. The left (respectively, right) panel assumes that the curvature … view at source ↗

**Figure 4.** Figure 4: Evolution of the energy densities for the choice of parameters leading to the power spectrum shown in the previous figure. The left (respectively, right) panel assumes that the curvature perturbation obeys Gaussian (respectively, χ 2 ) statistics. 4 PBH dark matter The PBHs of our interest form from the collapse of overdense regions at horizon re-entry [67]. The resulting PBH mass is controlled parametrica… view at source ↗

**Figure 5.** Figure 5: Left panel: PBH mass fraction obtained from the evolution described in the two previous figures, with black solid lines (receptively, red dashed lines) referring to the Gaussian (respectively, χ 2 ) case. The result is compared with current observational constraints on monochromatic PBH mass functions across different masses, taken from Ref. [78], showing that in both cases the PBHs obtained from this me… view at source ↗

**Figure 6.** Figure 6: Mass-weighted fraction of the universe that collapses to a PBH for any given mode, as a function of the amplitude of the power spectrum of the scalar curvature. We indicate the Gaussian and χ 2 limiting cases considered in the text. The radicand must be non-negative, which is equivalent to requiring the original overthreshold condition to hold. Switching the integration variables from dCg dζg to dMPBH dζg… view at source ↗

read the original abstract

We revisit the production of primordial black holes (PBHs) by a U(1) gauge field with a pseudo-scalar coupling to the inflaton. We improve upon the existing literature by working in the homogeneous backreaction regime with numerically computed gauge mode functions, adopting state-of-the-art PBH abundance calculations, and incorporating the uncertainty in the statistics of $\delta\rho$. We find that PBHs can account for all of the dark matter in the asteroidal mass range, even when the inflaton gradient energy density is highly subdominant ($10^{-4}$--$10^{-3}$ of the kinetic energy), supporting the validity of the backreaction scheme. This mechanism also unavoidably generates a stochastic gravitational wave background with an amplitude that will be measured at LISA and that will allow to indirectly discriminate between different statistics of $\delta \rho$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper finds PBHs can be all dark matter from axion inflation even with tiny gradient energy, but the homogeneous backreaction may not hold in the overdense patches that form them.

read the letter

The main result is that PBHs from a U(1) gauge field coupled to the inflaton can account for all dark matter in the asteroidal mass range, even when the inflaton gradient energy is only 10^{-4} to 10^{-3} of the kinetic energy. They reach this by solving the gauge mode functions numerically in the homogeneous backreaction regime, feeding the output into current PBH abundance formulas, and scanning the statistics of δρ to fold in uncertainty. A stochastic GW background comes along for free and sits in the LISA band, with amplitude that depends on which δρ statistics you pick. That gives a clear observational handle to test the setup later. The numerics and the updated abundance machinery are genuine improvements over earlier analytic approximations, and the subdominance claim is presented as a quantitative outcome of those upgrades. The LISA prediction is also a useful addition because it ties the mechanism to something measurable soon. The soft spot is the one the stress-test note flags. Gauge-field amplification depends exponentially on the local inflaton velocity and coupling strength. In the rare regions where δρ is large enough to collapse into a PBH, those local conditions can differ from the background, so the mode functions and backreaction there may not match the homogeneous calculation used for the global abundance. The paper varies the δρ statistics but does not appear to run a local-inhomogeneity check or estimate how much that shifts the tail. If the effect is order-one, the claimed abundance range could move. This is the sort of paper that belongs in the PBH-inflation and inflationary GW literature. Readers working on those topics will want the updated numerics and the LISA angle even if they end up disagreeing with the final abundance numbers. It is worth sending to peer review because the computational upgrades and the concrete GW signature give referees something concrete to evaluate, though the backreaction assumption in the relevant patches will need explicit discussion in revision.

Referee Report

2 major / 2 minor

Summary. The paper revisits PBH production in axion inflation coupled to a U(1) gauge field. It works in the homogeneous backreaction regime, numerically computes the gauge mode functions, adopts state-of-the-art PBH abundance formulas, and incorporates uncertainty in the statistics of δρ. The central claim is that PBHs can comprise all dark matter in the asteroidal mass window even when the inflaton gradient energy density is only 10^{-4}--10^{-3} of the kinetic energy, thereby supporting the backreaction scheme; the model also produces a SGWB detectable at LISA that can discriminate between δρ statistics.

Significance. If the result holds, the work supplies a concrete, testable PBH dark-matter scenario with an associated stochastic gravitational-wave background whose amplitude and shape are within LISA reach. Strengths include the use of numerically solved mode functions rather than analytic approximations, the adoption of modern abundance integrators, and the explicit propagation of statistical uncertainty in δρ; these elements reduce reliance on uncontrolled approximations common in earlier literature.

major comments (2)

[Section 3 (numerical mode functions and backreaction)] The homogeneous backreaction regime (used to compute the gauge mode functions that determine the δρ distribution): the global subdominance of gradient energy (10^{-4}--10^{-3}) does not automatically guarantee that the same regime holds inside the rare, high-δρ patches (δρ ≳ 0.1–0.3) that source PBHs. Because gauge-field amplification is exponentially sensitive to the local inflaton velocity, even a modest local enhancement of gradient or electric energy in those patches could alter the tail of the δρ PDF and therefore the predicted abundance. The manuscript does not provide a quantitative estimate or self-consistency check of this local feedback effect.
[Section 4 (PBH abundance and SGWB)] The mapping from the computed δρ statistics to the PBH mass function (Section 4): while state-of-the-art abundance tools are used and statistical uncertainty is varied, the input δρ distribution is still derived under the homogeneous-background assumption. Any correction to the high-δρ tail arising from the local-inhomogeneity issue raised above would propagate directly into the claimed “all of the dark matter” result and the LISA forecast.

minor comments (2)

[Figure 2] The caption of Figure 2 should explicitly state the range of δρ statistics explored and the precise definition of the “asteroidal mass window” used for the DM fraction.
[Section 3] A short paragraph comparing the numerical mode functions to the analytic WKB or slow-roll approximations employed in prior works would help readers assess the size of the improvement.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their careful reading, positive assessment of the manuscript's improvements, and constructive comments. We address the two major comments point by point below.

read point-by-point responses

Referee: [Section 3 (numerical mode functions and backreaction)] The homogeneous backreaction regime (used to compute the gauge mode functions that determine the δρ distribution): the global subdominance of gradient energy (10^{-4}--10^{-3}) does not automatically guarantee that the same regime holds inside the rare, high-δρ patches (δρ ≳ 0.1–0.3) that source PBHs. Because gauge-field amplification is exponentially sensitive to the local inflaton velocity, even a modest local enhancement of gradient or electric energy in those patches could alter the tail of the δρ PDF and therefore the predicted abundance. The manuscript does not provide a quantitative estimate or self-consistency check of this local feedback effect.

Authors: We agree that global subdominance alone does not rigorously prove the approximation holds locally in the rare, high-δρ patches. The gauge-field amplification depends exponentially on the inflaton velocity, so local velocity perturbations could in principle modify the tail. However, the velocity perturbation amplitude is controlled by the small gradient energy density; the fractional velocity fluctuation remains O(10^{-2}) even in the tail, limiting the change in the exponential growth factor. We will add a short paragraph in Section 3 that makes this scaling argument explicit and notes that a fully inhomogeneous simulation lies beyond present computational reach. This discussion clarifies the regime of validity without altering the central results. revision: partial
Referee: [Section 4 (PBH abundance and SGWB)] The mapping from the computed δρ statistics to the PBH mass function (Section 4): while state-of-the-art abundance tools are used and statistical uncertainty is varied, the input δρ distribution is still derived under the homogeneous-background assumption. Any correction to the high-δρ tail arising from the local-inhomogeneity issue raised above would propagate directly into the claimed “all of the dark matter” result and the LISA forecast.

Authors: The referee is correct that any modification to the high-δρ tail would affect the PBH abundance and the predicted SGWB. As explained in the response to the first comment, we expect the local correction to be perturbative given the smallness of the gradient energy. We will revise the discussion in Section 4 (and the conclusions) to explicitly state that the quoted abundance and LISA forecasts are obtained under the homogeneous backreaction assumption, while reiterating that the LISA signal shape remains a robust discriminator between different δρ statistics. No change to the numerical results is required. revision: partial

standing simulated objections not resolved

A fully self-consistent, inhomogeneous numerical evolution of the gauge-field modes inside the rare high-δρ patches is computationally prohibitive with current methods and cannot be performed within this work.

Circularity Check

0 steps flagged

No significant circularity: numerical mode computation and external PBH tools keep derivation self-contained

full rationale

The paper computes gauge mode functions numerically within the stated homogeneous backreaction regime, then inserts the resulting δρ statistics into independent, state-of-the-art PBH abundance formulas while varying the statistical assumptions. This chain does not reduce any output to a fitted parameter renamed as a prediction, nor does it rely on a self-citation that itself assumes the target result. The consistency check that PBH dark-matter abundance remains viable even for subdominant gradient energy (10^{-4}–10^{-3}) is an output of the calculation, not an input by definition. No load-bearing step matches any of the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review yields no extractable free parameters, axioms, or invented entities; paper likely inherits standard slow-roll inflation assumptions and PBH collapse criteria from cited literature.

pith-pipeline@v0.9.0 · 5444 in / 1084 out tokens · 38887 ms · 2026-05-07T09:01:45.001524+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

The Magnetic Origin of Primordial Black Holes: Ultralight PBHs and Secondary GWs
astro-ph.CO 2026-05 unverdicted novelty 5.0

Inflationary magnetic fields induce curvature perturbations that form ultralight PBHs, generating a stochastic GW background with model-specific features.

Reference graph

Works this paper leans on

80 extracted references · 76 canonical work pages · cited by 1 Pith paper · 4 internal anchors

[1]

Planck 2018 results. X. Constraints on inflation

Y. Akrami et al. (Planck), Astron. Astrophys.641, A10 (2020),1807.06211

work page internal anchor Pith review arXiv 2020
[2]

Byrnes, G

C. Byrnes, G. Franciolini, T. Harada, P. Pani, and M. Sasaki, eds.,Primordial Black Holes, Springer Series in Astrophysics and Cosmology (Springer, 2025), ISBN 978-981–978886-6, 978-981–978889-7, 978-981–978887-3

2025
[3]

B. Carr, K. Kohri, Y. Sendouda, and J. Yokoyama, Rept. Prog. Phys.84, 116902 (2021), 2002.12778

work page internal anchor Pith review arXiv 2021
[4]

Bartolo, V

N. Bartolo, V. De Luca, G. Franciolini, A. Lewis, M. Peloso, and A. Riotto, Phys. Rev. Lett. 122, 211301 (2019),1810.12218

work page arXiv 2019
[5]

Freese, J

K. Freese, J. A. Frieman, and A. V. Olinto, Phys. Rev. Lett.65, 3233 (1990). – 18 –

1990
[6]

M. S. Turner and L. M. Widrow, Phys. Rev. D37, 2743 (1988)

1988
[7]

W. D. Garretson, G. B. Field, and S. M. Carroll, Phys. Rev. D46, 5346 (1992), hep-ph/9209238

work page arXiv 1992
[8]

M. M. Anber and L. Sorbo, JCAP10, 018 (2006),astro-ph/0606534

work page arXiv 2006
[9]

Barnaby and M

N. Barnaby and M. Peloso, Phys. Rev. Lett.106, 181301 (2011),1011.1500

work page arXiv 2011
[10]

J. L. Cook and L. Sorbo, Phys. Rev. D85, 023534 (2012), [Erratum: Phys.Rev.D 86, 069901 (2012)],1109.0022

work page Pith review arXiv 2012
[11]

Barnaby, E

N. Barnaby, E. Pajer, and M. Peloso, Phys. Rev. D85, 023525 (2012),1110.3327

work page arXiv 2012
[12]

Domcke, M

V. Domcke, M. Pieroni, and P. Bin´ etruy, JCAP06, 031 (2016),1603.01287

work page arXiv 2016
[13]

Garcia-Bellido, A

J. Garcia-Bellido, A. Papageorgiou, M. Peloso, and L. Sorbo, JCAP01, 034 (2024), 2303.13425

work page arXiv 2024
[14]

¨Ozsoy, A

O. ¨Ozsoy, A. Papageorgiou, and M. Fasiello, JCAP12, 008 (2024),2405.14963

work page arXiv 2024
[15]

Gravitational waves from axion inflation in the gradient expansion formalism. Part I. Pure axion inflation

R. von Eckardstein, K. Schmitz, and O. Sobol (2025),2508.00798

work page internal anchor Pith review Pith/arXiv arXiv 2025
[16]

Teuscher, R

M. Teuscher, R. Durrer, K. Martineau, and A. Barrau, JCAP03, 043 (2026),2510.00869

work page arXiv 2026
[17]

Barbon, N

M. Barbon, N. Ijaz, and M. Peloso (2025),2510.17207

work page arXiv 2025
[18]

Linde, S

A. Linde, S. Mooij, and E. Pajer, Phys. Rev. D87, 103506 (2013),1212.1693

work page arXiv 2013
[19]

Axion inflation with gauge field production and primordial black holes

E. Bugaev and P. Klimai, Phys. Rev. D90, 103501 (2014),1312.7435

work page Pith review arXiv 2014
[20]

Garcia-Bellido, M

J. Garcia-Bellido, M. Peloso, and C. Unal, JCAP12, 031 (2016),1610.03763

work page arXiv 2016
[21]

Garcia-Bellido, M

J. Garcia-Bellido, M. Peloso, and C. Unal, JCAP09, 013 (2017),1707.02441

work page arXiv 2017
[22]

J. P. B. Almeida, N. Bernal, D. Bettoni, and J. Rubio, JCAP11, 009 (2020),2007.13776

work page arXiv 2020
[23]

¨Ozsoy and Z

O. ¨Ozsoy and Z. Lalak, JCAP01, 040 (2021),2008.07549

work page arXiv 2021
[24]

¨Ozsoy and G

O. ¨Ozsoy and G. Tasinato, Universe9, 203 (2023),2301.03600

work page arXiv 2023
[25]

Dimastrogiovanni, M

E. Dimastrogiovanni, M. Fasiello, and A. Papageorgiou, Phys. Rev. D110, 103542 (2024), 2403.13581

work page arXiv 2024
[26]

Domcke, V

V. Domcke, V. Guidetti, Y. Welling, and A. Westphal, JCAP09, 009 (2020),2002.02952

work page arXiv 2020
[27]

Peloso and L

M. Peloso and L. Sorbo, JCAP01, 038 (2023),2209.08131

work page arXiv 2023
[28]

Sobol, R

O. Sobol, R. von Eckardstein, E. Koch, S. Gurevich, U. Thiele, and K. Schmitz (2026), 2603.02570

work page arXiv 2026
[29]

Caravano, E

A. Caravano, E. Komatsu, K. D. Lozanov, and J. Weller, Phys. Rev. D105, 123530 (2022), 2110.10695

work page arXiv 2022
[30]

Caravano, E

A. Caravano, E. Komatsu, K. D. Lozanov, and J. Weller, Phys. Rev. D108, 043504 (2023), 2204.12874

work page arXiv 2023
[31]

D. G. Figueroa, J. Lizarraga, A. Urio, and J. Urrestilla, Phys. Rev. Lett.131, 151003 (2023), 2303.17436

work page arXiv 2023
[32]

Caravano and M

A. Caravano and M. Peloso, JCAP01, 104 (2025),2407.13405

work page arXiv 2025
[33]

Sharma, A

R. Sharma, A. Brandenburg, K. Subramanian, and A. Vikman, JCAP05, 079 (2025), 2411.04854

work page arXiv 2025
[34]

D. G. Figueroa, J. Lizarraga, N. Loayza, A. Urio, and J. Urrestilla, Phys. Rev. D111, 063545 (2025),2411.16368

work page arXiv 2025
[35]

Lizarraga, C

J. Lizarraga, C. L´ opez-Mediavilla, and A. Urio (2025),2505.19950

work page arXiv 2025
[36]

Jamieson, A

D. Jamieson, A. Caravano, and E. Komatsu (2025),2507.22285. – 19 –

work page arXiv 2025
[37]

Cheng, W

S.-L. Cheng, W. Lee, and K.-W. Ng, Phys. Rev. D93, 063510 (2016),1508.00251

work page arXiv 2016
[38]

Notari and K

A. Notari and K. Tywoniuk, JCAP12, 038 (2016),1608.06223

work page arXiv 2016
[39]

Dall’Agata, S

G. Dall’Agata, S. Gonz´ alez-Mart´ ın, A. Papageorgiou, and M. Peloso, JCAP08, 032 (2020), 1912.09950

work page arXiv 2020
[40]

He, K.-G

J.-F. He, K.-G. Zhang, C. Fu, and Z.-K. Guo, Phys. Rev. D111, 103525 (2025),2502.13158

work page arXiv 2025
[41]

O. O. Sobol, E. V. Gorbar, and S. I. Vilchinskii, Phys. Rev. D100, 063523 (2019),1907.10443

work page arXiv 2019
[42]

E. V. Gorbar, K. Schmitz, O. O. Sobol, and S. I. Vilchinskii, Phys. Rev. D104, 123504 (2021), 2109.01651

work page arXiv 2021
[43]

Durrer, O

R. Durrer, O. Sobol, and S. Vilchinskii, Phys. Rev. D108, 043540 (2023),2303.04583

work page arXiv 2023
[44]

von Eckardstein, M

R. von Eckardstein, M. Peloso, K. Schmitz, O. Sobol, and L. Sorbo, JHEP11, 183 (2023), 2309.04254

work page arXiv 2023
[45]

Domcke, Y

V. Domcke, Y. Ema, and S. Sandner, JCAP03, 019 (2024),2310.09186

work page arXiv 2024
[46]

Durrer, R

R. Durrer, R. von Eckardstein, D. Garg, K. Schmitz, O. Sobol, and S. Vilchinskii, Phys. Rev. D 110, 043533 (2024),2404.19694

work page arXiv 2024
[47]

Maleknejad and E

A. Maleknejad and E. Komatsu, JHEP05, 174 (2019),1808.09076

work page arXiv 2019
[48]

Ishiwata, E

K. Ishiwata, E. Komatsu, and I. Obata, JCAP03, 010 (2022),2111.14429

work page arXiv 2022
[49]

Iarygina, E

O. Iarygina, E. I. Sfakianakis, R. Sharma, and A. Brandenburg, JCAP04, 018 (2024), 2311.07557

work page arXiv 2024
[50]

Dimastrogiovanni, M

E. Dimastrogiovanni, M. Fasiello, A. Papageorgiou, and C. Z. Gatica (2025),2504.17750

work page arXiv 2025
[51]

Ishiwata and E

K. Ishiwata and E. Komatsu, JCAP03, 052 (2026),2512.10184

work page arXiv 2026
[52]

Cielo, E

M. Cielo, E. Dimastrogiovanni, M. Fasiello, and A. Papageorgiou (2026),2601.08622

work page arXiv 2026
[53]

Musco, Phys

I. Musco, Phys. Rev. D100, 123524 (2019),1809.02127

work page arXiv 2019
[54]

Musco, V

I. Musco, V. De Luca, G. Franciolini, and A. Riotto, Phys. Rev. D103, 063538 (2021), 2011.03014

work page arXiv 2021
[55]

Young, C

S. Young, C. T. Byrnes, and M. Sasaki, JCAP07, 045 (2014),1405.7023

work page arXiv 2014
[56]

Young, I

S. Young, I. Musco, and C. T. Byrnes, JCAP11, 012 (2019),1904.00984

work page arXiv 2019
[57]

De Luca, G

V. De Luca, G. Franciolini, A. Kehagias, M. Peloso, A. Riotto, and C. ¨Unal, JCAP07, 048 (2019),1904.00970

work page arXiv 2019
[58]

Ferrante, G

G. Ferrante, G. Franciolini, A. Iovino, Junior., and A. Urbano, Phys. Rev. D107, 043520 (2023),2211.01728

work page arXiv 2023
[59]

A. D. Gow, H. Assadullahi, J. H. P. Jackson, K. Koyama, V. Vennin, and D. Wands, EPL142, 49001 (2023),2211.08348

work page arXiv 2023
[60]

P. A. R. Ade et al. (Planck), Astron. Astrophys.594, A17 (2016),1502.01592

work page arXiv 2016
[61]

Adshead, J

P. Adshead, J. T. Giblin, M. Pieroni, and Z. J. Weiner, Phys. Rev. D101, 8 (2020), 1909.12842

work page arXiv 2020
[62]

Adshead, J

P. Adshead, J. T. Giblin, M. Pieroni, and Z. J. Weiner, Phys. Rev. Lett.124, 17 (2020), 1909.12843

work page arXiv 2020
[63]

A. R. Liddle and D. H. Lyth, Phys. Rept.231, 1 (1993),astro-ph/9303019

work page Pith review arXiv 1993
[64]

Planck 2018 results. VI. Cosmological parameters

N. Aghanim et al. (Planck), Astron. Astrophys.641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)],1807.06209

work page internal anchor Pith review arXiv 2020
[65]

Barnaby, R

N. Barnaby, R. Namba, and M. Peloso, JCAP04, 009 (2011),1102.4333. – 20 –

work page arXiv 2011
[66]

Riotto, ICTP Lect

A. Riotto, ICTP Lect. Notes Ser.14, 317 (2003),hep-ph/0210162

work page arXiv 2003
[67]

Suyama and C.-M

T. Suyama and C.-M. Yoo,Overall Picture: A Beginner’s Guide to Primordial Black Hole Formation(Springer Nature Singapore, Singapore, 2025), pp. 37–51, ISBN 978-981-97-8887-3, URLhttps://doi.org/10.1007/978-981-97-8887-3_2

work page doi:10.1007/978-981-97-8887-3_2 2025
[68]

Black hole formation in the Friedmann universe: Formulation and computation in numerical relativity

M. Shibata and M. Sasaki, Phys. Rev. D60, 084002 (1999),gr-qc/9905064

work page Pith review arXiv 1999
[69]

Musco,Numerical Simulations of Primordial Black Holes(Springer Nature Singapore, Singapore, 2025), pp

I. Musco,Numerical Simulations of Primordial Black Holes(Springer Nature Singapore, Singapore, 2025), pp. 93–112, ISBN 978-981-97-8887-3, URL https://doi.org/10.1007/978-981-97-8887-3_4

work page doi:10.1007/978-981-97-8887-3_4 2025
[70]

M. W. Choptuik, Phys. Rev. Lett.70, 9 (1993)

1993
[71]

C. R. Evans and J. S. Coleman, Phys. Rev. Lett.72, 1782 (1994),gr-qc/9402041

work page arXiv 1994
[72]

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

I. Musco, J. C. Miller, and A. G. Polnarev, Class. Quant. Grav.26, 235001 (2009),0811.1452

work page Pith review arXiv 2009
[73]

Escriv` a, Universe8, 66 (2022), 2111.12693

A. Escriv` a, Universe8, 66 (2022),2111.12693

work page arXiv 2022
[74]

The abundance of primordial black holes depends on the shape of the inflationary power spectrum

C. Germani and I. Musco, Phys. Rev. Lett.122, 141302 (2019),1805.04087

work page Pith review arXiv 2019
[75]

Escriv` a, C

A. Escriv` a, C. Germani, and R. K. Sheth, Phys. Rev. D101, 044022 (2020),1907.13311

work page arXiv 2020
[76]

Young, Int

S. Young, Int. J. Mod. Phys. D29, 2030002 (2019),1905.01230

work page arXiv 2019
[77]

Primordial Black Holes - Perspectives in Gravitational Wave Astronomy -

M. Sasaki, T. Suyama, T. Tanaka, and S. Yokoyama, Class. Quant. Grav.35, 063001 (2018), 1801.05235

work page Pith review arXiv 2018
[78]

B. Carr, A. J. Iovino, G. Perna, V. Vaskonen, and H. Veerm¨ ae (2026),2601.06024

work page arXiv 2026
[79]

Cosmological long-wavelength solutions and primordial black hole formation

T. Harada, C.-M. Yoo, T. Nakama, and Y. Koga, Phys. Rev. D91, 084057 (2015),1503.03934

work page Pith review arXiv 2015
[80]

Kehagias, I

A. Kehagias, I. Musco, and A. Riotto, JCAP12, 029 (2019),1906.07135. – 21 –

work page arXiv 2019