{"id":"e7e436c3-d759-4295-86ac-111d629648d9","arxiv_id":"2412.20859","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Dark matter shells collapsing nonradially through the center of a primordial black hole binary accelerate its merger, boosting the merger rate by a factor of a few.","lead":"Primordial black hole pairs may merge faster than previously thought because streams of dark matter pass through the center of their halo and drain orbital energy. This raises the predicted gravitational-wave merger rate by a factor of a few and lowers the dark-matter fraction needed to explain LIGO/Virgo events.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The factor-of-few merger-rate enhancement rests on an unvalidated disk model for collisionless DM shell crossing; Eq. (7) and Eq. (26) need a phase-space or N-body check.","rationale":"The reader's weakest assumption is correct and is the same point I would put under stress. The paper's own Sections 7-8 acknowledge that the effect is computed in an approximate analytical model; no high-resolution simulation or code is offered. The central density rho_c in Eq. (7) is the only input that changes the merger rate relative to the no-halo calculation, and it comes from a geometric ansatz whose phase-space justification is only a citation to [21]. I do not see an internal contradiction in the algebra, and the direction of the effect (shells passing through the center increase hardening) is physically plausible; however, the quantitative claim 'several times' and the resulting f range are not robustly established. A purpose-built N-body experiment can settle whether the effective density at the binary is within a factor of order unity of Eq. (7). If it is not, the claimed reduction of f to 3-7 x 10^-4 would need revision. This supports, rather than overturns, the reader's CONDITIONAL verdict.","tokens_in":10041,"tokens_out":19183,"duration_ms":197656,"concrete_test":"Run a high-resolution collisionless N-body simulation of one 30 M_sun PBH binary with initial conditions at z_eq from the same tidal model: inflationary perturbations normalized as in Section 3 plus the nearest third PBH at the mean separation for f = 10^-4. Use enough particles per DM shell and small enough softening to resolve the first shell crossing around the binary; measure the DM density in a sphere of radius ~a around the binary and the induced d(1/a)/dt as functions of time. Compare these with Eq. (7) and Eq. (26). If the measured hardening rate differs by more than a factor of 2 from the disk-model prediction for f = 10^-4, the f ~ (3-7) x 10^-4 conclusion is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in Section 8, that f ~ (3-7) x 10^-4 suffices because the DM halo accelerates PBH mergers by 'several times,' depends directly on the central DM density rho_c(t) from Eq. (7) and on using that density in the hardening law, Eq. (26). Equation (7) is obtained from the disk ansatz in Section 2: a first-contracting spherical shell is treated as a thin uniform disk of radius r_t, with thickness set only by the spread in turnaround times (delta_r_c = v_c delta_t_c, Eq. (5)). For collisionless DM this is not a derived phase-space result. The first infall produces caustic structures whose density near the binary is not shown to be a homogeneous disk; the zero-angular-momentum locus in a tidal field may be a line or filament, and the velocity distribution at crossing is strongly anisotropic. Reference [21] is cited for the existence of zero-angular-momentum directions and for a figure, not for the disk density profile. This uncertainty is not a small order-unity refinement: rho_c enters d(1/a)/dt directly in Eq. (26), and through the semimajor-axis reduction it enters the merger probability nonlinearly. The coefficient H1 ~ 20 is itself taken from a stellar-dynamical fit for a Maxwellian velocity background, not from a monoenergetic cold stream; applying it to the disk flow adds an additional unquantified O(1) error. Since the paper provides no error bars and no code, the magnitude 'several times' is plausible but not yet anchored.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":10396,"tokens_out":6455,"duration_ms":68438,"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":[{"comment":"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).","section":"§2, Eqs. (3)-(7)"},{"comment":"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.","section":"§5, Eq. (26)"},{"comment":"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.","section":"Fig. 3 and §8"}],"minor_comments":[{"comment":"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...\".","section":"Abstract and §1"},{"comment":"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.","section":"§4, text after Eq. (20)"},{"comment":"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.","section":"Fig. 3"},{"comment":"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.","section":"Eqs. (18) and (23)"},{"comment":"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.","section":"§5, around Eq. (27)"},{"comment":"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.","section":"§3"}],"recommendation":"major_revision","confidential_remarks":"I see no issue with the novelty or integrity of the manuscript, and it is well within the scope of the journal. My concern is purely technical: the main quantitative claim rests on an unvalidated disk model and an extrapolated hardening coefficient. A revised version that includes a numerical or phase-space check of Eq. (7), a justification or sensitivity study of H1 in Eq. (26), and an uncertainty estimate for the merger-rate curves would make the factor-of-several claim convincing."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this paper is a legitimate order-of-magnitude estimate, not a breakthrough. The new piece is treating nonradial DM shell contraction as producing a persistent central flow that hardens PBH binaries, and quantifying it with a disk model. The authors are upfront that the disk is an ansatz, and they offer a conservative loss-cone estimate (Section 7) that doesn't depend on it. That internal cross-check is the best part of the paper: both routes give some halo enhancement, and the qualitative direction agrees with the simulations in [23].\n\nThe main soft spot is exactly what the stress-test flags: Eq. (7) for the central density comes from a thin uniform disk with thickness set only by turnaround-time spread. Collisionless shell crossing could produce caustics, filaments, or a cored profile, and the density near the binary could differ by more than the claimed factor of a few. The hardening law (26) also uses H1 ~ 20 from a Maxwellian stellar-dynamical fit, applied to a cold monoenergetic stream; that's an O(1) uncertainty on top. No error bars, no code. So the abstract's 'several times' is not firmly anchored.\n\nThat said, the central claim doesn't collapse. The authors made a reasonable analytic choice, they didn't hide its approximate nature, and they gave a second, more conservative route to the same qualitative conclusion. The f ~ (3-7) x 10^-4 range should be read as indicative, not precise.\n\nWho is this for? People working on PBH merger-rate predictions and LIGO/Virgo interpretations. It deserves a serious referee: the question is well-posed, the citation pattern is fair, and the main deficiency is a missing uncertainty budget and a missing validation of the disk model. I'd send it to review, with a request that the authors either compare against an N-body run that resolves the background flow or state clearly a regime where the disk ansatz breaks down.","headline":"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.","tokens_in":10954,"tokens_out":2216,"would_cite":true,"duration_ms":23469,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["primordial black holes","dark matter halos","binary black hole mergers","gravitational waves","LIGO/Virgo/KAGRA","dark matter shell contraction","loss cone","tidal forces"],"falsifier":"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.","tokens_in":9784,"feed_emoji":"🌌","tokens_out":5742,"duration_ms":52414,"temperature":0.7,"pith_summary":"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.","feed_headline":"Dark halo streams speed up black hole pair mergers severalfold","feed_subtitle":"A PBH abundance of 0.03–0.07 percent of dark matter matches LIGO/Virgo/KAGRA rates.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the fixed-point theorem guaranteeing directions of zero angular momentum in DM shells and the resulting disk-like first passage through the halo center, the geometric basis for the central density.","marker":"[21]"},{"why":"Provides the hardening law d(1/a)/dt = H₁ Gρ/v_c that converts the central DM density into orbital semimajor-axis decay.","marker":"[25]"},{"why":"Gives the earlier order-of-magnitude estimate that the third PBH directs DM particles into the loss cone of the pair, which this paper extends by modeling the halo growth around the third PBH and the shell passage.","marker":"[18]"},{"why":"Provides the numerical simulation result that PBH pair mergers are accelerated by DM halos, which this paper confirms and compares against, showing an additional enhancement from background DM passage.","marker":"[23]"},{"why":"Supplies the standard model for PBH pair formation, the probability distribution of pair separations, and the merger-rate statistics that the paper modifies by including the halo influence.","marker":"[9–11]"}],"fun_headline_variants":["Dark matter halos boost PBH merger rates","Halo flows accelerate black hole pair mergers","Nonradial halo infall speeds PBH coalescence","PBH mergers severalfold faster in DM halos"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Dark matter halos boost PBH merger rates","Halo flows accelerate black hole pair mergers","Nonradial halo infall speeds PBH coalescence","PBH mergers severalfold faster in DM halos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0002,"raw_usage":{"total_tokens":1353,"prompt_tokens":902,"completion_tokens":451,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":518,"completion_tokens_details":{"reasoning_tokens":402}},"tokens_in":518,"tokens_out":451,"duration_ms":5055,"temperature":1.0,"reasoning_tokens":402,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T23:08:55.934919+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the fixed-point theorem guaranteeing directions of zero angular momentum in DM shells and the resulting disk-like first passage through the halo center, the geometric basis for the central density."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the hardening law d(1/a)/dt = H₁ Gρ/v_c that converts the central DM density into orbital semimajor-axis decay."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the earlier order-of-magnitude estimate that the third PBH directs DM particles into the loss cone of the pair, which this paper extends by modeling the halo growth around the third PBH and the shell passage."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the numerical simulation result that PBH pair mergers are accelerated by DM halos, which this paper confirms and compares against, showing an additional enhancement from background DM passage."}],"review_version":1}