Pith. sign in

REVIEW 3 major objections 6 minor 33 references

Mergers of Binary Primordial Black Holes in Evolving Dark Matter Halos

T0 review · 3 major / 6 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Dark matter halos that assemble around pairs of primordial black holes continuously feed dark matter through the halo center, extracting orbital energy and raising the merger rate by a factor of several.

desk verdict A plausible order-of-magnitude estimate that background DM streaming through halo centers hardens PBH binaries, but the factor-of-few rests on an unvalidated disk ansatz and needs a phase-space or N-body check. read the letter →

arxiv 2412.20859 v1 pith:7VARITHA submitted 2024-12-30 astro-ph.CO

V.D. Stasenko , Yu.N. Eroshenko This is my paper · ORCID
classification astro-ph.CO
keywords primordialblackholesdarkmatterhalosbinaryholemergersgravitationalwavesLIGO/Virgo/KAGRAshellcontractionlossconetidalforces
topics Dark Matter
open problems Dark Matter
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper argues that the dark matter halo forming around a pair of primordial black holes does not stay static: as shells of dark matter contract, a continuous stream of DM passes through the halo center, where it scatters off the black holes and drains orbital energy. That extra energy loss hardens the binary and shortens its merger time. For primordial black holes of about 30 solar masses, including this halo feedback raises the predicted merger rate by several times compared with calculations that ignore the halo. Consequently, a primordial black hole abundance of only about f ≈ (3–7) × 10⁻⁴ of the dark matter would be enough to explain the LIGO/Virgo/KAGRA merger rates.

What carries the argument

The disk ansatz: each infalling dark matter shell is represented at first core passage as a thin uniform disk of radius r_t passing through the halo center, with central density ρ_c(t) = (π/8) ρ̄(t̃) (r_s(t̃)/r_t)². The radius r_t is set by a quadrature sum of tidal displacements from inflationary perturbations, the third nearest PBH, and the pair's radius of influence, with r_t capped by the shell stopping radius. This disk density, inserted into the Quinlan hardening rate d(1/a)/dt = H₁ Gρ/v_c, is what converts the continuous central stream into an orbital energy loss that accelerates the merger.

What would settle it

Run a high-resolution N-body simulation of a 30 M☉ PBH pair embedded in a cosmological DM halo, with enough mass resolution to track smooth background DM shells (not just discrete clumps), and measure the DM density at the halo center and the binary's semimajor-axis decay rate. If the measured central density is much lower than Eq. (7) at redshifts z ≈ 500–3000, or if the hardening rate is not several times higher than in the no-halo case, the disk ansatz and the factor-of-several enhancement would be refuted.

Watch

Extended reading notes

Core claim

The central claim is that nonspherical, nonradial contraction of dark matter shells produces a nonzero, slowly varying dark matter density at the center of the halo surrounding a PBH pair. Because every shell has a direction of zero angular momentum, a shell collapses through the halo center before dispersing, and because new shells continuously detach from the cosmological expansion, the pair is perpetually immersed in a central DM flow. The paper derives this central density as ρ_c(t) = (π/8) ρ̄(t̃) (r_s(t̃)/r_t)², with the disk radius r_t determined by tidal forces from inflationary density perturbations and from the third nearest PBH, and combines it with the Quinlan hardening law to compute accelerated orbital decay. Including both the disk-like shell passage and a more conservative loss-cone estimate raises the merger rate several times above the no-halo case, confirming the qualitative conclusion of earlier numerical simulations but with a larger enhancement.

Load-bearing premise

The calculation assumes that a first-contracting dark matter shell can be represented as a thin uniform disk of radius r_t passing through the halo center with a specific central density; this disk picture is an analytic approximation, not derived from a full phase-space or N-body treatment, so the true central density and the resulting orbital hardening could differ by more than the claimed factor of a few.

Editorial extensions

If this is right

  • The PBH fraction needed to match LIGO/Virgo/KAGRA rates drops to roughly (3–7) × 10⁻⁴ when the halo-fed central DM stream is included.
  • The predicted merger rate with halo feedback (solid curve in Fig. 3) is several times higher than the no-halo rate for f ≲ 10⁻², with the enhancement most significant where inflationary tidal perturbations dominate over the third-PBH effect.
  • The effect operates only while DM shells are still detaching from the cosmological expansion; it ceases once halo growth stops, so the enhancement is tied to the assembly epoch of the halo.
  • A conservative loss-cone calculation that ignores the one-dimensional shell passage still yields a non-negligible enhancement, bounding the expected rate if the disk picture is not fully realized.
  • The authors note that additional factors, such as PBH clustering, could raise the admissible PBH fraction beyond the quoted range, so the rate estimate is a lower bound within this model.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the central DM stream is strongest during halo assembly, the redshift distribution of PBH mergers should show an enhancement at high redshifts; a dedicated search for a high-z merger excess could test this mechanism indirectly.
  • If the disk ansatz fails — for example, if shell crossing produces caustics, a cored density profile, or a strongly velocity-dependent flow — the factor-of-several enhancement could shrink or vanish, so the quantitative result hinges on the accuracy of the disk representation.
  • The same central-stream mechanism should apply to any binary embedded in an assembling dark matter halo, not only PBH pairs; binary black holes of intermediate mass in early halos could show a comparable orbital-hardening effect, testable with high-resolution cosmological simulations that resolve smooth background DM flow.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper studies how the dark-matter halo around a primordial-black-hole pair affects the pair's orbital evolution and merger rate. Its central mechanism is that a collapsing DM shell is not spherical: along directions of zero angular momentum the shell passes through the halo center, so a continuously renewed DM flow of density rho_c(t) (Eq. 7) is present at the pair. The authors compute rho_c from tidal perturbations produced by inflationary density fluctuations and by the third PBH, use the hardening law Eq. (26) for the semimajor axis, and then compute the merger rate including gravitational-radiation decay. They find that the halo increases the rate "by several times" and that a PBH fraction f ~ (3-7)x10^-4 suffices to match the LIGO/Virgo/KAGRA rate.

Significance. If the mechanism is real, the paper identifies a previously under-appreciated smooth-DM contribution to PBH-binary hardening and gives a direct route to lower inferred PBH abundances. The analytic model is transparent and includes a conservative alternative estimate based on the loss cone, and the comparison with the simulation of [23] is a useful anchor. However, the headline factor-of-several depends on two unvalidated ingredients: the thin-disk ansatz for collisionless shell crossing and the use of a Maxwellian stellar-dynamical hardening coefficient for a cold stream. The paper is therefore a plausible order-of-magnitude estimate rather than an established prediction until those ingredients are tested. There is no code or detailed reproducibility information, but the formulas are explicit enough that the calculation could be reconstructed.

major comments (3)
  1. [§2, Eqs. (3)-(7)] The central density rho_c(t) is derived from the thin-disk ansatz: a first-contracting shell is represented as a uniform disk of radius r_t and thickness delta_r_c, with delta_r_c/delta_r_s = 9 pi / 8. For collisionless dark matter this is not a derived phase-space result. Reference [21] is cited for the existence of zero-angular-momentum directions and for a figure, but not for a uniform disk density profile. First infall can produce caustics, filaments, or a cored flow, and a velocity-dependent flow would change rho_c by more than the claimed factor of a few. Since rho_c enters linearly in Eq. (26) and nonlinearly in the merger probability, the quantitative claim of the paper requires an explicit phase-space or N-body check of Eq. (7).
  2. [§5, Eq. (26)] The hardening coefficient H1 ~ 20 is taken from Quinlan (1996), a stellar-dynamical fit obtained for a Maxwellian velocity distribution of incoming stars. In this paper the DM crossing the binary is a cold, coherent stream at one velocity v_c. The hardening rate for such a stream need not equal the Maxwellian-averaged rate, and the resulting O(1) error is unquantified. Because the factor-of-several enhancement in Fig. 3 is directly proportional to the hardening rate in Eq. (26), the authors should either derive H1 for their monoenergetic stream geometry or show that the final rate is insensitive to H1 over a plausible range.
  3. [Fig. 3 and §8] The conclusion that f ~ (3-7)x10^-4 is enough to explain the observed merger rate is based on comparing the solid and dotted curves in Fig. 3, which differ by only a factor of 2-3, and no error bars or sensitivity tests are given. Equations (18), (23), (25), and (26) each carry order-of-magnitude uncertainties, and the curves are close. The authors should propagate uncertainties through the rate calculation, for example by varying H1, the shell thickness ratio delta_r_c/delta_r_s, and the quadrature combination in Eq. (25), before the quoted f range can be regarded as robust. The paper's own characterization of the model as an "approximate analytical model" in §1 and "for our estimate" in §2 is not carried through to the quantitative claim in §8.
minor comments (6)
  1. [Abstract and §1] The phrase "around pair of primordial black holes" should be "around a pair of primordial black holes"; also, "each shell upon the first contraction passes through the halo center" is stated as a fact before the disk approximation is introduced, so it would be clearer to say "is modeled as passing...".
  2. [§4, text after Eq. (20)] The sentence beginning "the characteristic comoving distance to the third PBH is y_ch ~ xbar/2" is duplicated and left incomplete; it should be rewritten as a single sentence.
  3. [Fig. 3] The horizontal axis is labeled "log f" with tick values from -4 to -2; please specify that the base is 10 and describe in the caption what the dashed and dash-dotted curves represent.
  4. [Eqs. (18) and (23)] The function kappa(delta_eq) is not derived anywhere; since the text quotes kappa -> 2.9 as delta_eq -> 0, a brief derivation or plot of kappa would help the reader reproduce the figures.
  5. [§5, around Eq. (27)] The quantity K1 is left unspecified and the reference to [25] is not enough to make Eq. (27) self-contained; because the eccentricity change is subsequently neglected, an explicit statement of the relevant K1 form would be useful.
  6. [§3] The paper does not state the cosmological parameters used for rho_eq, z_eq, and the Planck-normalized power spectrum; specifying them would improve reproducibility of the numerical integration.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the central DM-halo enhancement is computed from independent tidal and hardening inputs; the disk ansatz is an approximation, not a self-referential fit.

full rationale

No load-bearing circular step was found. The claimed merger-rate enhancement is assembled from independent inputs. The central dark-matter density rho_c(t) in Eq. (7) is computed from the shell-contraction geometry together with the tidal radii r_t,in and r_t,3 obtained from the inflationary perturbation spectrum and the third-PBH separation in Eqs. (18) and (23); it is not fitted to the final merger rate. The hardening law d(1/a)/dt = H1 G rho / v_c in Eq. (26) is imported from Quinlan's numerical scattering experiments, and the no-halo baseline in Eq. (32) is the separate standard formula of Refs. [9-11]. The fraction f is scanned rather than fitted to the LIGO/Virgo/KAGRA band, so the conclusion that f ~ (3-7) x 10^-4 suffices is a derived prediction rather than a re-statement of an input. The only self-citations, Refs. [19,20], appear in the introduction and conclusions as contextual remarks about additional PBH-clustering effects; they do not carry the central derivation. The disk/shell-crossing approximation in Section 2 is a genuine modeling assumption whose quantitative accuracy is not demonstrated, and applying the H1 coefficient outside its Maxwellian-averaged context introduces possible error; however, these are correctness or calibration risks, not cases where a prediction equals its input by construction. The paper explicitly defers full validation to future high-resolution simulations, which further confirms that the present claim is an extrapolation rather than a circular restatement.

Assumptions & free parameters 4 free parameters · 6 assumptions · 0 invented entities

The central rate enhancement rests on a chain of standard cosmological perturbation and binary-hardening results, plus one ad hoc disk model for shell crossing. The ad hoc disk model is the most consequential item: it converts a one-dimensional focus into a density enhancement of order (rs/rt)^2. The free parameters are order-unity coefficients inherited from prior numerical fits, not tuned to the target rate; MPBH and f are physical inputs, with f scanned. No new particles, forces, or conserved quantities are introduced.

free parameters (4)
  • H1 (binary hardening coefficient) = ~20 (from Quinlan 1996)
    Appears in Eq (26) for the semimajor-axis decay induced by passing DM shells; taken from numerical experiments in [25], not derived here, and directly scales the halo-induced hardening rate.
  • alpha, beta (semimajor and semiminor axis coefficients) = ~1
    Appear in Eq (30) for the initial binary orbit from pair-formation theory [9-11]; order-unity values are assumed and enter the merger-time formula.
  • eta (loss-cone energy extraction efficiency) = ~1
    Used in the conservative estimate in Section 7 for the energy extracted per DM particle in the loss cone; not determined in the paper.
  • MPBH (PBH mass) = 30 M_sun
    Fixed as the representative LIGO/Virgo mass; the paper does not scan over mass, so the reported rates are specific to this value.
assumptions (6)
  • domain assumption Spherical collapse solution around a PBH: DM shells follow the parametric equations (1)-(2) with delta_eq = 2 MPBH/M.
    This is the standard spherical top-hat collapse model adapted to a point PBH; it fixes the turn-around radius, collapse time, and shell density used throughout Section 2.
  • domain assumption For every contracting DM shell there is at least one zero-angular-momentum direction, and this direction passes through the halo center, so shells cross the center continuously.
    Inherited from [21] via the fixed point theorem and the accompanying figure; the central-density calculation in Eq (7) collapses without this assumption.
  • ad hoc to paper A crossing shell can be represented as a thin uniform disk of transverse radius rt and normal thickness delta_rc, with density given by Eq (3) and delta_rc/delta_rs = 9 pi / 8.
    This disk ansatz is the paper's key modeling step; it is justified qualitatively by Fig 2 of [21], not derived from the full phase-space distribution.
  • domain assumption Tidal perturbations from inflationary density fluctuations and the third PBH act independently, with transverse displacement obtained by integrating Eq (12) along the unperturbed radial trajectory; total rt is obtained from quadrature in Eq (25).
    Linear perturbation theory is used, and the min and quadrature combination is an approximation; no nonlinear corrections or coupling between tidal sources are included.
  • domain assumption The third nearest PBH is distributed with flat probability (31), and the pair's initial semimajor and semiminor axes are given by Eq (30) with alpha and beta of order unity.
    Standard pair-formation statistics from [9-11]; used to convert the per-binary lifetime into the merger rate in Section 6.
  • domain assumption The DM halo mass around the third PBH grows as MH(z) = (3/2)(2/(3 pi))^(2/3) ((1 + z_eq)/(1 + z)) MPBH (Eq 22).
    Taken from self-similar halo growth; enters the tidal force from the third PBH and the conservative loss-cone estimate.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Mergers of Binary Primordial Black Holes in Evolving Dark Matter Halos." pith.science (2026). https://pith.science/paper/7VARITHA

@misc{pith2026241220859,
  author       = {Pith},
  title        = {Pith review of: Mergers of Binary Primordial Black Holes in Evolving Dark Matter Halos},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7VARITHA}},
  note         = {Machine review of arXiv:2412.20859}
}
abstract

The influence of a dark matter halo around pair of primordial black holes on their orbit evolution and the black hole merger rate is considered. Because of the nonspherical (nonradial) contraction of DM shells, each shell upon the first contraction passes through the halo center in the direction of the radius vector corresponding to zero angular momentum. Since the shell contraction is a continuous process, at each instant of time there is a nonzero dark matter density at the halo center. This density is determined by the influence of the tidal gravitational forces from inflationary density perturbations and from other primordial black holes. The scattering of dark matter particles by a pair of black holes leads to a loss of the energy of its orbital motion and to an accelerated pair merger. In the case of primordial black holes with masses $\sim30M_\odot$, the black hole merger rate in the presence of a dark matter halo is several times higher than that without such a halo.

Figures

Figures reproduced from arXiv: 2412.20859 by the authors.

Figure 1
Figure 1. The characteristic DM shell radius rt determined by the combined effect of the tidal forces from the inflationary perturbations and the third nearest PBH with respect to the shell stopping radius rs at redshift z. The curves correspond to the PBH mass MPBH = 30M⊙ and (from top to bottom) f = 10−1 , f = 10−2 , f = 10−3 , and f = 10−4 [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. The density at the center of the DM halo sur [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. The rate of gravitational bursts from the mergers o [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

33 extracted references · 33 canonical work pages

  1. [23]

    Eroshenko Yu.N., Astronomy Letters 42, 347 (2016)

    work page 2016

  2. [21]

    Carr B., Kuhnel F., Annual Review of Nuclear and Par- ticle Science 70, 355 (2020)

    work page 2020

  3. [1]

    A multitude of studies of the possible role of PBHs in various astrophysi- cal processes and in cosmology have also been performed

    INTRODUCTION The idea about the formation of primordial black holes (PBHs) in the early Universe was proposed almost 60 years ago by [1–3], and several models of their forma- tion have been proposed since then [3–8]. A multitude of studies of the possible role of PBHs in various astrophysi- cal processes and in cosmology have also been performed. As examp...

    work page

  4. [2]

    eq” and the quantities at the present time t0 by the subscript “0

    THE DENSITY AT THE CENTER OF THE HALO SURROUNDING A PBH P AIR A region with a density ∼ ρeq and a radius ∼ (3MPBH/(4πρeq))1/ 3 is formed around a PBH pair near the transition time teq of the Universe to the dust-like stage of evolution (hereafter we mark the quantities at the time teq by the subscript “eq” and the quantities at the present time t0 by the ...

    work page

  5. [3]

    In this section we will find the part of rt associated with the influence of inflationary DM density perturbations

    INFLUENCE OF THE TIDAL FORCES PRODUCED BY INFLATIONAR Y PER TURBATIONS The quantity rt can contain the contribution from var- ious tidal processes. In this section we will find the part of rt associated with the influence of inflationary DM density perturbations. A DM halo is formed around a PBH; the DM shells successively pass from expansion to contraction....

    work page

  6. [4]

    radius of influence

    INFLUENCE OF THE THIRD PBH The angular momentum of the PBH pair is determined by the third nearest PBH. However, this third PBH in- fluences the motion of the DM shells. Near teq the DM halo mass is still small and the inflationary perturbations did not grow. Therefore, one might expect that at larger f the third PBH will make a relatively large contribu- t...

    work page 2000

  7. [5]

    EVOLUTION OF THE PBH P AIR ORBIT Using a numerical experiment, [25] obtained the law of evolution of the semimajor axis of a black hole pair un- der the influence of incoming stars with the same velocity and then performed an averaging over the velocities un- der the assumption of their Maxwellian distribution. In our case, it is necessary to use the expre...

    work page

  8. [6]

    Denote the comoving separations be- tween the components of the PBH pair and between the center of mass of the pair and the third PBH by x and y, respectively

    THE PBH MERGER RATE Following the approach of [9–11], let us write out the basic relations required to calculate the statistics of PBH mergers in pairs by taking into account the influence of the DM halo. Denote the comoving separations be- tween the components of the PBH pair and between the center of mass of the pair and the third PBH by x and y, respect...

    work page

Show all 33 references
  1. [7]

    Particles with a low an- gular momentum can pass through the orbit of the PBH pair, causing the energy of its orbital motion to be lost

    A CONSER V ATIVE ESTIMATE OF THE MERGER RATE In this section, for comparison, we will give a more con- servative estimate of the merger rate without including the one-dimensional DM shell contraction, but by taking into account the angular momentum distribution of DM particles...

  2. [8]

    The DM parti- cles in the halo forming around the PBH pair acquire an angular momentum under the influence of tidal grav- itational forces

    CONCLUSIONS In this paper we considered the influence of the DM flow continuously passing through the center of a DM halo on the evolution of a PBH pair. The DM parti- cles in the halo forming around the PBH pair acquire an angular momentum under the influence of tidal grav- itat...

  3. [9]

    and Novikov I.D., Sov

    Zel’dovich Ya.B. and Novikov I.D., Sov. Astron. 10, 602 (1967)

    1967

  4. [10]

    Hawking S., MNRAS 15, 75 (1971)

    1971

  5. [11]

    Khlopov M.Yu., Polnarev A.G., Phys. Lett. B 97 , 383 (1980)

    1980

  6. [12]

    Berezin V.A., Kuzmin V.A., Tkachev I.I., Phys. Lett. B 120, 91 (1983)

    1983

  7. [13]

    Jedamzik K., Phys. Rev. D 55, R5871(R) (1997)

    1997

  8. [14]

    Dolgov A., Silk J., Phys. Rev. D 47, 4244 (1993)

    1993

  9. [15]

    Cos- mol

    Rubin S.G., Khlopov M.Yu., Sakharov A.S., Grav. Cos- mol. S 6, 51 (2000)

    2000

  10. [16]

    Rubin S.G., Sakharov A.S., Khlopov M.Yu., JETP 92, 921 (2001)

    2001

  11. [17]

    Nakamura T., Sasaki M., Tanaka T., Thorne K.S., The Astrophys. J. 487, L139 (1997)

    1997

  12. [18]

    Ioka K., Chiba T., Tanaka T., Nakamura T., Phys. Rev. D 58, 063003 (1998)

    1998

  13. [19]

    Sasaki M., Suyama T., Tanaka T., Yokoyama S., Phys. Rev. Lett. 117, 061101 (2016)

    2016

  14. [20]

    Dolgov A. D. et al., Journal of Cosmology and Astropar- ticle Physics 12, 017 (2020)

    2020

  15. [22]

    Phys.: Conf

    Eroshenko Yu.N., J. Phys.: Conf. Ser. 1051, 012010 (2018)

    2018

  16. [24]

    Dokuchaev V.I., Eroshenko Y.N., Astronomy Letters 27, 759 (2001)

    2001

  17. [25]

    Dokuchaev V.I., Eroshenko Y.N., Astronomical & Astro- physical Transactions 22, 727 (2003)

    2003

  18. [26]

    Hayasaki K., Takahashi K., Sendouda Y., Nagataki S., Publications of the Astronomical Society of Japan 68, 66 (2016)

    2016

  19. [27]

    Stasenko V., Belotsky K., MNRAS 526, 4308 (2023)

    2023

  20. [28]

    Stasenko V., Phys. Rev. D 109, 123546 (2024)

    2024

  21. [29]

    Sikivie P., Tkachev I.I., Wang Y., Phys. Rev. D 56, 1863 (1997)

    1997

  22. [30]

    Tkachev M.V., Pilipenko S.V., Yepes G., Mon. Not. R. Astron. Soc. 499, 4854 (2020)

    2020

  23. [31]

    Pilipenko S., Tkachev M., Ivanov P., Phys. Rev. D 105, 123504 (2022)

    2022

  24. [32]

    Berezinsky V., Dokuchaev V., Eroshenko Y., Phys. Rev. D 68, 103003 (2003)

    2003

  25. [33]

    1, 35 (1996)

    Quinlan G.D., New Astron. 1, 35 (1996)

    1996

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.