Non-spherical effects on the mass function of Primordial Black Holes

Basic Information

• Paper ID: 2410.03451
• Title: Non-spherical effects on the mass function of Primordial Black Holes
• Authors: Albert Escrivà, Chul-Moon Yoo
• Classification: gr-qc astro-ph.CO hep-th
• Publication Date: October 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2410.03451

Abstract

This paper investigates the impact of non-spherical effects on the mass function of primordial black holes (PBHs) through 3+1 relativistic numerical simulations based on ellipticity-dependent critical thresholds. The study considers both the radiation equation of state w=P/ρ=1/3 and a softer equation of state w=1/10. Assuming curvature perturbations follow Gaussian statistics with a monochromatic power spectrum, peak theory is employed to examine the most probable ellipsoidal configurations. The findings demonstrate that while non-sphericity readily prevents black hole formation in the critical scaling regime, its impact on the mass function is relatively modest, reducing it by only a few factors, with negligible effects on the mass function distribution, including the power-law scaling in the low-mass region.

Research Background and Motivation

Core Problem

The fundamental question addressed by this research is: How do non-spherical effects influence the mass function of primordial black holes (PBHs)? Traditional studies predominantly assume spherical symmetry, whereas actual curvature perturbations are typically ellipsoidal in nature.

Significance

1. Dark Matter Candidate: PBHs represent an important dark matter candidate, particularly in the asteroid mass range (10^-15 - 10^-10 M⊙)
2. Theoretical Prediction Accuracy: Accurate predictions of PBH abundance are crucial for understanding the early universe and the nature of dark matter
3. Numerical Verification Requirements: The validity of the spherical symmetry assumption requires verification through non-spherical numerical simulations

Limitations of Existing Methods

• Most studies rely on spherical symmetry assumptions, lacking systematic investigation of non-spherical effects
• Only a limited number of non-spherical numerical simulation studies exist 1, 36
• Previous theoretical predictions 48 are based on heuristic models lacking numerical verification

Core Contributions

1. First Systematic Study: A comprehensive investigation combining 3+1 relativistic numerical simulations with peak theory to study non-spherical effects on the PBH mass function
2. Ellipticity-Dependent Threshold: Establishes a quantitative relationship between ellipticity e and critical threshold μc: ec(μ) = Ke(μ-μc,sp/μc,sp)^γe
3. Mass Function Correction: Demonstrates that non-spherical effects reduce the mass function by only 1-2 factors, with negligible impact on the overall distribution
4. Critical Scaling Preservation: Verifies that the critical scaling law M ∝ (μ-μc)^γ remains valid in the non-spherical case

Methodology Details

Task Definition

Investigate gravitational collapse of ellipsoidal curvature perturbations during the radiation-dominated era and compute the PBH mass function accounting for non-spherical effects.

Input: Ellipsoidal perturbation parameters (μ, e, p) Output: Corrected PBH mass function fPBH(M) Constraint: Super-horizon condition kp/(aHb) ≪ 1

Model Architecture

1. Ellipsoidal Perturbation Description

Employs the characteristic profile from peak theory:

ζ̄(x,y,z;μ) = ζ̄sp(r;μ) + 5μ/(2kp³r³)[3kpr cos(kpr) + (r³kp³-3)sin(kpr)] × [3e/r²(z²-y²) + p(1-3(x/r)²)]

Where:

• μ: peak amplitude
• e: ellipticity
• p: oblateness
• ζ̄sp(r;μ) = μ sinc(kpr): spherically symmetric profile

2. Ellipticity-Dependent Threshold

Critical ellipticity determined through numerical simulations:

ec(μ) = Ke((μ-μc,sp)/μc,sp)^γe

Parameter values:

• w=1/10: Ke≈0.532, γe≈0.413
• w=1/3: Ke≈0.563, γe≈0.448

3. Corrected Mass Scaling Law

Assumes the critical scaling law remains valid in the non-spherical case:

M = MH(μ)K(μ-μc(e))^γ

Where μc(e) is the function satisfying e = ec(μc(e)).

Technical Innovations

1. Ellipticity Probability Distribution: Utilizes peak theory to compute the conditional probability distribution of ellipticity and oblateness Pe,p(e,p|ν,ξ)
2. Mass Function Integration: Employs Jacobian transformation |J|(M,e) = ∂μ/∂lnM to convert peak amplitude to mass
3. Monochromatic Spectrum Approximation: In the γ→1 limit, simplifies integration using Dirac delta functions

Experimental Setup

Numerical Simulation Parameters

• Power Spectrum Peak: kp = 10^13.5 Mpc^-1
• Effective Degrees of Freedom: g* = 106.75
• Equation of State: w = 1/3 (radiation), w = 1/10 (soft fluid)
• Mass Range: Asteroid mass range

Evaluation Metrics

• PBH Mass Function: fPBH(M) = MnPBH(M)/ρDM
• Total PBH Fraction: f^tot_PBH = ∫fPBH(M)dlnM = 1

Comparison Methods

• Spherically Symmetric Case: Conventional calculations with e = p = 0
• Non-Spherical Case: Complete calculations accounting for ellipticity distribution

Experimental Results

Main Results

Numerical Parameters

• w=1/10: A≈1.242×10^-3, νc≈8.8, Mkp≈8.36×10^-15 M⊙
• w=1/3: A≈4.870×10^-3, νc≈8.7, Mkp≈1.22×10^-14 M⊙

Mass Function Comparison

From Figure 3, the following observations emerge:

1. Overall Suppression: Non-spherical effects reduce the mass function by approximately 1-2 factors
2. Shape Preservation: The overall shape of the mass function remains essentially unchanged
3. Power-Law Preservation: The power-law behavior in the low-mass tail fPBH(M) ∝ M^(1+1/γ(w)) is preserved

Key Findings

1. Concentrated Ellipticity Distribution: For large peak values (ν≫1), the ellipticity distribution concentrates at small values: ⟨e⟩≈6×10^-2, ⟨p⟩≈10^-3
2. Critical Region Dominance: The mass function is dominated by perturbations near the critical value μc(e)
3. Minimal Equation of State Dependence: Results for w=1/3 and w=1/10 are similar

Related Work

Major Research Directions

1. Spherically Symmetric PBH Formation: Extensive numerical studies 21-35 determine critical thresholds in the spherically symmetric case
2. Non-Spherical Studies: Limited works 1, 36, 44-46 investigate non-spherical collapse
3. Peak Theory Applications: 14, 57, 58 apply peak theory to PBH abundance calculations

Advantages of This Work

• First systematic study combining 3+1 numerical simulations with peak theory
• Provides quantitative relationships for ellipticity-dependent thresholds
• Quantifies the actual impact of non-spherical effects

Conclusions and Discussion

Main Conclusions

1. Limited Impact: Non-spherical effects have relatively modest impact on the PBH mass function (1-2 factors)
2. Validity of Spherical Approximation: The spherical symmetry assumption is highly accurate for calculating PBH abundance
3. Critical Scaling Preservation: The critical scaling law remains valid in the non-spherical case
4. Universality: Conclusions hold for both radiation and soft equations of state

Limitations

1. Monochromatic Spectrum Restriction: Only considers monochromatic power spectra; extension to broad spectra is needed
2. Gaussian Perturbation Assumption: Does not account for non-Gaussian perturbations
3. Specific Profile: Employs a specific profile from peak theory; other profiles may yield different results
4. Mass Scaling Assumption: Assumes the critical scaling law is unaffected by ellipticity

Future Directions

1. Investigate non-spherical effects for broad power spectra
2. Account for the influence of non-Gaussian perturbations
3. Verify mass scaling laws for large ellipticity cases
4. Examine effects of alternative profile functions

In-Depth Evaluation

Strengths

1. Methodological Innovation: First systematic combination of numerical simulations and peak theory for non-spherical effects
2. Technical Rigor: Employs complete 3+1 relativistic numerical simulations with reliable results
3. Theoretical Completeness: Provides a comprehensive theoretical framework from ellipticity thresholds to mass function
4. Practical Value: Validates the spherical approximation and provides important reference for PBH research

Weaknesses

1. Parameter Space Limitations: Considers only two equations of state and specific parameter ranges
2. Assumption Dependence: Key results depend on the assumption that critical scaling laws hold in the non-spherical case
3. Insufficient Numerical Verification: Does not directly verify mass scaling relationships in the non-spherical case

Impact

1. Theoretical Contribution: Provides important non-spherical corrections to PBH formation theory
2. Observational Significance: Affects constraints on PBHs as dark matter candidates
3. Methodological Reference: Provides methodological guidance for subsequent non-spherical studies

Applicable Scenarios

• Precise calculations of PBH abundance
• Investigation of dark matter formation mechanisms in the early universe
• Theoretical predictions of gravitational wave sources
• Cosmological parameter constraints

References

This paper cites 61 important references covering classical and frontier works in multiple fields including PBH formation, numerical relativity, and peak theory. Particularly noteworthy is the authors' companion paper 1, which provides the numerical simulation foundation for this research.