Primordial Black Hole Formation in $f(R)=R+\alpha R^2$ Gravity: Perturbative and Non-Perturbative Analysis

This paper investigates primordial black hole formation in quadratic $f(R)$ gravity by combining a first-order perturbative analysis around General Relativity with a non-perturbative Einstein-frame numerical study of the scalaron field to determine the critical overdensity threshold for collapse.

The Gist

The Big Picture: Gravity's "Extra Gear"

Imagine General Relativity (our current best theory of gravity) as a standard car engine. It works perfectly for driving on normal roads (like planets orbiting stars). But the authors of this paper are asking: What happens if we add a turbocharger?

In this study, the "turbocharger" is a specific mathematical tweak to gravity called $f(R) = R + \alpha R^2$.

• $R$ represents the curvature of space (how bent space is).
• $\alpha$ is a tiny knob that controls how strong the "turbo" is.
• When space is flat or gently curved, the turbo does nothing, and gravity acts like normal.
• But when space gets extremely bent (like right before a black hole forms), the turbo kicks in, changing how gravity behaves.

The paper investigates what happens when a giant cloud of dust collapses under its own weight to form a black hole, specifically looking at how this "turbo" changes the process.

---

Part 1: The "Dust Cloud" Experiment (Perturbative Analysis)

The researchers first looked at a simplified scenario: a collapsing cloud of "dust" (matter with no pressure, like a pile of sand). They treated the "turbo" as a very small addition to normal gravity to see the first-order effects.

The Analogy: The Race to the Finish Line
Imagine two runners starting a race to collapse a cloud into a black hole:

1. Runner A (Normal Gravity/GR): Runs at a steady, predictable pace.
2. Runner B (Modified Gravity): Has a tiny bit of extra energy (the $\alpha$ term).

The Finding:
The paper found that Runner B finishes faster.

• The "turbo" makes the cloud collapse quicker than it would in normal gravity.
• Because the cloud shrinks faster, the "finish line" (the event horizon, the point of no return) is crossed sooner.
• The Consequence: If you need a certain amount of "push" (density) to start a collapse, and the gravity is stronger/faster, you actually need less push to get the job done. The paper suggests this means it becomes easier to form black holes in this theory.

The Twist (The Radiation Case):
The researchers also tried this with a cloud of "radiation" (like light or hot gas) instead of dust.

• The Result: In this specific simplified model, the "turbo" didn't work at all for the radiation cloud. The curvature of space in a radiation-dominated universe is different, and the math showed that the extra term canceled out.
• The Takeaway: To see the "turbo" effect on radiation, you can't use simple math; you need to look at the messy, complex, real-world chaos (non-linear effects).

---

Part 2: The "Hidden Engine" (Non-Perturbative Analysis)

Since the simple math had limits, the authors switched to a different way of looking at the problem, called the Einstein Frame.

The Analogy: Changing the Camera Angle
Imagine you are watching a movie of a car crash.

• The first method was like watching from a distance, trying to guess what happened by looking at the smoke.
• The second method (Einstein Frame) is like putting a camera inside the engine.

In this view, the "turbo" isn't just a tweak to gravity; it reveals a hidden particle called the scalaron.

• Think of the scalaron as a spring-loaded weight attached to the universe.
• When the universe is calm, the spring is relaxed.
• When the universe gets squeezed (like during the formation of a black hole), the spring gets compressed and pushes back, changing the dynamics of the collapse.

The authors wrote down a complete set of rules (equations) describing how this spring (scalaron) moves alongside the collapsing cloud. They didn't solve these equations with a computer in this paper, but they provided the blueprint so others can do it. This blueprint allows scientists to calculate exactly how much easier it is to form a black hole under these extreme conditions.

---

Part 3: What Does This Mean for the Universe? (Observational Constraints)

If this "turbo" makes black holes form too easily, we should see a lot more of them than we do.

The Analogy: The Goldilocks Zone

• If the "turbo" is too weak, we don't see the effect.
• If the "turbo" is too strong, we would have a universe filled with black holes, which would mess up the cosmic microwave background (the afterglow of the Big Bang) and the light from distant stars.
• The Paper's Conclusion: By looking at how many black holes we actually see (or don't see), we can put a limit on how big the "turbo" knob ($\alpha$) can be.
• The paper suggests that if the "turbo" is too strong, it would create too many black holes, violating what we observe. Therefore, the value of $\alpha$ must be very small, or it must behave differently in the very early universe compared to today.

Summary of Key Points

1. Faster Collapse: In the presence of this specific gravity tweak, dust clouds collapse faster than in normal gravity.
2. Easier Black Hole Formation: Because collapse is faster, the threshold (the minimum density needed) to create a black hole is likely lower.
3. Radiation is Tricky: In a simple model, radiation doesn't show this effect, meaning the real physics is more complex and requires advanced computer simulations.
4. The Blueprint: The authors provided the mathematical "blueprint" (ODE system) for the "hidden spring" (scalaron) so future scientists can run the numbers to predict exactly how many black holes should exist.
5. Real-World Check: Observations of the universe (like the lack of too many black holes) tell us that this "gravity turbo" cannot be too powerful, or it would have created a universe that doesn't look like ours.

What the paper does NOT do:

• It does not claim to have found a new type of black hole.
• It does not provide a final, exact number for how many black holes exist.
• It does not apply this to medical technology or everyday life; it is strictly about the physics of the early universe and black holes.

Technical Summary

Technical Summary: Primordial Black Hole Formation in $f(R) = R + \alpha R^2$ Gravity

Problem Statement
The study addresses the formation of Primordial Black Holes (PBHs) within the context of the quadratic $f(R)$ gravity model, specifically the Starobinsky model defined by $f(R) = R + \alpha R^2$. While General Relativity (GR) provides a well-established framework for PBH formation—where overdense regions collapse if the density contrast $\delta$ exceeds a critical threshold $\delta_c \approx 0.4\text{--}0.5$ during the radiation-dominated era—the dynamics of gravitational collapse may be significantly altered in modified gravity theories. The central problem is to determine how the curvature correction term $\alpha R^2$ modifies the collapse dynamics of overdense regions, the timing of horizon formation, and the resulting critical threshold $\delta_c$ for PBH production. The authors aim to bridge the gap between perturbative corrections near the GR limit and the fully non-perturbative dynamics required in high-curvature regimes, such as those near the end of inflation.

Methodology
The paper employs a dual approach, combining perturbative analytic methods with a non-perturbative formulation in the Einstein frame:

1. Perturbative Expansion around GR:

   • The authors expand the field equations of the $f(R) = R + \alpha R^2$ model to first order in the small parameter $\alpha$.
   • They analyze a collapsing, spatially flat Friedmann-Lemaître-Robertson-Walker (FLRW) interior filled with dust ($p=0$) and radiation ($p=\rho/3$).
   • The metric is expanded as $g_{\mu\nu} = g^{(0)}_{\mu\nu} + \alpha h_{\mu\nu}$, where $g^{(0)}_{\mu\nu}$ is the standard GR solution. The linearized field equations are derived to compute first-order corrections to the scale factor $a(t)$, the Hubble parameter $H(t)$, and the stellar radius.
   • This homogeneous treatment is explicitly noted as a simplified proxy to isolate curvature effects on background dynamics, rather than a direct first-principles calculation of the inhomogeneous PBH threshold.

2. Einstein-Frame Non-Perturbative Formulation:

   • To access the high-curvature regime where perturbative methods break down, the theory is reformulated in the Einstein frame via a conformal transformation.
   • In this frame, the model is equivalent to GR plus a canonical scalar field (the scalaron, $\phi$) with a Starobinsky potential $V(\phi) = \frac{1}{8\alpha}(1 - e^{-\sqrt{2/3}\phi})^2$.
   • The authors derive a complete system of Ordinary Differential Equations (ODEs) governing both the cosmological background and the evolution of a separate, closed FLRW patch representing an overdense region.
   • This system includes the scalaron evolution, Hubble friction, and dissipation terms, allowing for the numerical integration required to determine the critical overdensity $\delta_c(k)$ near the end of inflation.

Key Results

• Dust Collapse Dynamics: In the perturbative regime ($\alpha R \ll 1$), the $R^2$ correction accelerates the gravitational collapse of a dust-dominated FLRW interior. The corrected scale factor is found to be $a(t) = a_0(t)[1 - \alpha u^{-2}]$, where $u = t_0 - t$. This implies that for $\alpha > 0$, the physical radius of the collapsing cloud shrinks faster than in GR, leading to an earlier horizon formation time.
• Radiation Background Limitation: For a flat, radiation-dominated background, the Ricci scalar $R^{(0)}$ vanishes identically in the GR limit. Consequently, the correction tensor $K_{\mu\nu}$, which depends on $R$ and its derivatives, vanishes at first order. The authors conclude that within this specific homogeneous perturbative framework, the radiation-dominated background receives no nontrivial first-order correction. Any modification to PBH formation in the radiation era must therefore arise from nonlinear or non-perturbative effects.
• PBH Threshold Trend: Based on the accelerated collapse in the dust case, the authors infer a qualitative trend: the modification of collapse dynamics suggests a reduction in the effective critical overdensity threshold, $\delta_c^{(\alpha)} \lesssim \delta_c^{GR}$. While Eq. (47) provides a heuristic relation linking the shift in collapse time to $\delta_c$, the authors emphasize that this is not a precise quantitative derivation of the threshold, which requires full nonlinear inhomogeneous simulations.
• Critical Exponent Modification: The study proposes a leading-order estimate for the modification of the critical scaling exponent $\gamma$ in the mass relation $M_{PBH} \propto (\delta - \delta_c)^\gamma$. The curvature correction is estimated to reduce the exponent, $\gamma^{(\alpha)} < \gamma_{GR}$, implying that PBHs formed slightly above the threshold would be lighter, sharpening the mass function toward low masses.
• Observational Constraints: The paper translates the potential enhancement of PBH abundance (due to a reduced $\delta_c$) into constraints on the parameter $\alpha$. Using the Press-Schechter formalism, the authors note that even a small reduction in $\delta_c$ leads to an exponential increase in the PBH mass fraction $\beta(M)$. Comparing this with observational limits from LIGO/Virgo, microlensing surveys, CMB anisotropies, and evaporation constraints, they derive an order-of-magnitude bound $\alpha \lesssim 10^{-2} M_{Pl}^{-2}$. This is contrasted with the much larger $\alpha$ required for Starobinsky inflation ($\alpha \sim 10^9 M_{Pl}^{-2}$), suggesting that the parameter governing inflation may differ from the effective parameter in the high-curvature collapse regime.

Significance and Claims
The paper claims to provide a "complete analytic and semi-analytic study" of PBH formation in the Starobinsky model. Its primary significance lies in:

1. Explicit Analytic Corrections: Deriving the first-order corrections to the scale factor, Hubble parameter, and horizon formation time for collapsing dust, demonstrating explicitly that $\alpha > 0$ accelerates collapse.
2. Methodological Framework: Providing the full ODE system in the Einstein frame that allows for the future quantitative determination of $\delta_c(k)$ across all scales, bridging the gap between perturbative insights and non-perturbative numerical requirements.
3. Qualitative Insight into Thresholds: Offering a heuristic link between curvature corrections and the PBH formation threshold, suggesting that modified gravity can enhance PBH production in high-curvature regimes.
4. Distinction of Regimes: Clarifying that while dust collapse is modified at first order, the radiation-dominated background is not, highlighting the necessity of non-perturbative treatments for realistic PBH formation scenarios in the early universe.

The authors maintain a modest tone regarding their quantitative claims, repeatedly stating that the perturbative results provide "qualitative insight" and "heuristic indications" rather than precise values for $\delta_c$, which require the nonlinear numerical integration of the derived ODE system. The work positions itself as a foundational step connecting modified gravity dynamics to compact object formation, consistent with existing literature on scalar-tensor PBH studies while adding specific analytic derivations for the $R^2$ model.