{"id":"02e9a84d-608d-4d8a-a998-51a808194215","arxiv_id":"2508.00774","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Massive primordial black holes with strong gas-dark matter streaming can trigger direct-collapse black hole formation in their wake, producing early supermassive black hole binaries.","lead":"This paper simulates a new way supermassive black hole pairs could form in the early universe: a heavy primordial black hole acts as a seed, and its radiation plus the relative motion of gas and dark matter help a second, smaller black hole form nearby. If the channel is real, it gives JWST, ALMA, and future space-based gravitational wave detectors concrete systems to search for.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The q~0.1 binary prediction rests on an unresolved single-progenitor assumption; sub-pc fragmentation could replace the DCBH with a cluster and erase the claimed mass ratio.","rationale":"The reader's weakest assumption identifies exactly the load-bearing issue: the paper's headline binary mass ratio and separation follow from the unresolved assumption that the sub-pc inflow feeds a single protostar. I agree with that assessment. The paper is internally consistent and honestly flags the limitation in Section 3.2, so this is an unresolved condition rather than a contradiction. Because the condition is not checked by the present simulations, the conditional verdict is appropriate; no shift to accept or reject is warranted on this basis alone. A secondary issue noted in the reader's rationale—the abstract and conclusions quote sustained inflow rates of 0.01-0.1 Msun/yr while Section 3.2 reports sustained values of only ~1e-3 Msun/yr with peaks above ~1e-2—reinforces the need for caution but is less central than the fragmentation question. The proposed high-resolution zoom-in test would settle whether the single-progenitor assumption holds; until then, the q~0.1 and ~10 pc predictions should be treated as conditional.","tokens_in":22712,"tokens_out":10435,"duration_ms":144867,"concrete_test":"Re-simulate the PBH LW str fd005 initial conditions at the moment the first sink forms, using a zoom-in AMR or SPH run that resolves <0.01 pc and includes H2/H- chemistry, the same LW feedback prescription, and sink particles with accretion. Count distinct sinks and track the mass of the most massive object over 5 Myr; if multiple sinks form and no single object reaches ~1e4-1e5 Msun, or the accretion rate onto the primary sink falls below ~0.01 Msun/yr, the single-progenitor assumption fails and q~0.1 is not supported. A cheaper first check is to compute the enclosed specific angular momentum and Toomre Q profile below 1 pc in the existing snapshots; Q<1 at r<0.5 pc would indicate fragmentation is likely.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative prediction—a secondary black hole of ~1e5 Msun and binary mass ratio q~O(0.1)—requires that the gas crossing into the sink region feed a single protostar. The simulations resolve only scales above roughly 0.5 pc physical (softening 0.01 h^-1 kpc, physical ~0.5 pc at z~17), and Section 3.2 explicitly states that sub-pc turbulence, radiative transfer, and angular momentum transport are unresolved and could instead produce a rotationally supported disk or multiple cores. If fragmentation-induced starvation occurs, the cloud forms a cluster rather than a single supermassive star; even if a massive black hole later assembles via collisions, its final mass and the binary mass ratio are no longer the simulated m_col ~5e4-1e5 Msun and q~0.1. The abstract and conclusions present q~0.1 and ~10 pc separations as outcomes, not merely as one possibility, so this unresolved scale is load-bearing. The paper cites competitive accretion and stellar collisions as possible rescue mechanisms, but those mechanisms themselves require unresolved sub-pc physics, so they broaden the uncertainty rather than close it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses cosmological hydrodynamical simulations of an isolated ~10^6 Msun primordial black hole to argue that accretion-driven Lyman-Werner feedback, combined with baryon-dark matter streaming velocities v_bchi >= 0.8 sigma_bchi, suppresses H2 cooling and triggers the collapse of atomically-cooling gas clouds into ~10^5 Msun direct-collapse black hole seeds at z ~ 20-10. The authors identify dense collapsing cores in runs with streaming 0.8-1.6 sigma, infer high gas infall rates, and interpret the ratio of the collapsing-core mass to the PBH mass as forming SMBH binaries with q ~ 0.05-0.1 and separations ~10 pc. The paper then discusses potential electromagnetic and gravitational-wave signatures, including links to Little Red Dots and LISA/TianQin sources.","tokens_in":22965,"tokens_out":5120,"duration_ms":60520,"significance":"If the proposed channel holds, it offers a novel in situ route to high-redshift SMBH binaries, with concrete, falsifiable predictions for mass ratios, separations, and multi-messenger observability. The manuscript is transparent about many of its limitations, builds on a well-established simulation code, and makes its initial-condition generator publicly available. However, the central quantitative predictions rest on several under-tested or unresolved ingredients: the sustained inflow-rate claim is weaker than the abstract states, the single-progenitor assumption is unresolved at sub-pc scales, and key parameters and thresholds are tested with only a single realization or a single fiducial value. The significance is therefore conditional on closing these gaps.","major_comments":[{"comment":"The abstract and Section 5 state that the collapsing clouds exhibit sustained inflow rates Mdot_infall >= 0.01-0.1 Msun/yr, which is the canonical DCBH threshold, but the body of Section 3.2 and Figure 4 report rates that consistently exceed ~1e-3 Msun/yr with peaks frequently surpassing ~1e-2 Msun/yr. The sustained rate is therefore an order of magnitude below the quoted critical range; only peaks reach it. This mismatch is load-bearing because the 'ideal conditions for DCBH formation' claim is the bridge from simulated clouds to supermassive stars, and the abstract currently overstates what the simulations show.","section":"§3.2, Fig. 4, Abstract"},{"comment":"The binary mass ratio q ~ 0.05-0.1 is obtained by equating the simulated collapsing-cloud mass with the mass of a single secondary black hole, under the explicit assumption that the inflow feeds a single protostar. The same section acknowledges that sub-pc processes such as turbulence, radiative transfer, and angular momentum transport are unresolved and could instead produce a rotationally supported disk or fragmentation into multiple cores; if fragmentation-induced starvation operates, the secondary mass and q are not the simulated m_col. Because Section 4 presents q ~ 0.1 and ~10 pc separations as baseline predictions, this unresolved scale is load-bearing for the central claim rather than a minor caveat.","section":"§3.2 and §4"},{"comment":"The ad hoc veto on sink formation at z>200 is introduced to avoid numerical artifacts, but no test is shown demonstrating that high-redshift dense structures are indeed spurious rather than physical collapse events. Since the proposed channel targets z ~ 20-10, the veto may not change the final conclusions, yet the paper should justify it with a resolution study or by following one vetoed event; otherwise the collapse criterion is not consistently applied.","section":"§2.3"},{"comment":"The threshold v_bchi >= 0.8 sigma_bchi is inferred from a single 0.4 sigma run without collapse and a single 0.8 sigma run with collapse. The paper itself notes the stochastic nature of collapse and the lack of a clear trend with streaming amplitude. With one realization per velocity, the apparent threshold could be sampling noise; multiple realizations per streaming value, or an analytic criterion, are needed to support a 'critical regime' claim.","section":"§2.1, Table 1, §3.1"},{"comment":"All successful collapse runs use the fiducial thermal feedback coupling epsilon_r = 0.005; the runs that vary epsilon_r (0.05, 0.005, 0.0005) have no streaming and do not collapse, while the streaming runs hold epsilon_r fixed. Since epsilon_r is a free parameter and the paper identifies it as critical in prior work, the central result is not demonstrated to be robust to epsilon_r in the regime where collapse actually occurs.","section":"Table 1, §2.2"}],"minor_comments":[{"comment":"Please define all symbols in the equation caption and state explicitly whether J_LW is an unshielded or self-shielded intensity; the text later applies a shielding factor, but Eq. (4) appears to give unshielded values.","section":"Equation (4)"},{"comment":"The blue star marking the collapsing cloud and the black dot marking the PBH are very small; larger markers or a zoomed inset would make the ~10 pc separation and the relative velocity vectors clearer.","section":"Figure 2"},{"comment":"The reference list contains two distinct 2025 entries by Zhang, Liu, and Bromm (one Zenodo, one ApJ); please use year-letter suffixes and disambiguate in the text so the reader can tell which prior work is being cited.","section":"References"},{"comment":"The entries without collapse use '-' for z_col, but the footnote explaining symbols does not define '-'; please add a short definition for clarity.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript presents a creative and potentially important scenario, but the abstract overstates the simulation results relative to the body, and several load-bearing assumptions (single-progenitor collapse, epsilon_r sensitivity, single-realization streaming threshold) need to be addressed before the quantitative predictions can be accepted. I do not recommend rejection because the central idea is sound enough to warrant further work, but the required revisions are substantial."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe genuinely new thing here is the specific in situ route: a ~1e6 Msun PBH with accretion LW feedback plus baryon-DM streaming that exceeds ~0.8 sigma can make an atomic-cooling cloud collapse in the PBH wake, giving a secondary ~1e5 Msun seed at ~10 pc separation and q~0.1 at z~20-10. That combination is not in the earlier PBH seeding literature, including the authors' own 2025 paper. It is a real idea worth engaging.\n\nThe paper does several things well. The simulations are internally consistent, the PBH initial conditions and accretion/LW treatment are built on documented prior work with public code, and the authors repeatedly flag their resolution limits. They also honestly describe the alternative fate: sub-pc fragmentation into a cluster, with competitive accretion or collisions as possible rescue mechanisms. The observational discussion (JWST/ALMA signatures, LISA/TianQin relevance) is clearly separated from the simulation results.\n\nThe soft spots are in proportion. The load-bearing detail is the unresolved sub-pc fate: the q~0.1 result essentially assumes one central protostar gets all the inflow. The 1-10 pc resolution cannot tell you whether the gas instead forms a disk or fragments. The stress-test note is correct that the rescue mechanisms themselves require unresolved physics, so the uncertainty is broad, not narrow. The abstract/body mismatch is real: the abstract states sustained inflow of 0.01-0.1 Msun/yr and q~0.1 as outcomes, while Section 3.2 presents these as conditional on a single-progenitor assumption. That should be fixed. Also minor: only one realization per streaming velocity, a coarse 0.4 vs 0.8 sigma threshold, a hand-tuned epsilon_r=0.005, and an artificial z>200 sink veto. None of these individually kills the channel, but collectively they mean the claim is a plausible proof-of-concept, not a validated prediction.\n\nThis paper deserves a serious referee. The channel is novel, the simulations are reproducible enough to check, and the authors have identified the key physics (streaming-enhanced wake collapse plus PBH feedback) that could matter for early SMBH binaries. I would send it to review with a request that the authors soften the abstract, clearly label q and separation as resolution-dependent estimates, and ideally add a second realization or a convergence test. I would not cite the q~0.1 value as established until higher-resolution work confirms it.\n\nRecommendation: engage, send to peer review, require the overclaims be trimmed.","headline":"A plausible PBH-catalyzed DCBH binary channel with honest caveats; the q~0.1 prediction is a resolution-limited inference, not a secure result.","tokens_in":23503,"tokens_out":646,"would_cite":false,"duration_ms":10790,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Simulations show a 10^6-solar-mass primordial black hole can trigger the collapse of a gas cloud into a 10^5-solar-mass companion, forming a supermassive black hole binary at redshift 20-10.","keywords":["primordial black holes","direct-collapse black holes","supermassive black hole binaries","Lyman-Werner radiation","baryon-dark matter streaming","high-redshift galaxy formation","gravitational waves","Little Red Dots"],"falsifier":"A radiation-hydrodynamic simulation that resolves the sub-parsec collapse and finds the atomically-cooling cloud fragments into many low-mass protostars, starving any single protostar of accretion, would remove the DCBH seed and with it the predicted $q\\sim0.1$ binary.","tokens_in":22486,"feed_emoji":"🌌","tokens_out":11245,"duration_ms":110835,"temperature":0.7,"pith_summary":"This paper argues that a single massive primordial black hole can build a supermassive black hole binary on its own, without a galaxy merger. In hydrodynamical simulations, a $\\sim10^6\\,M_\\odot$ PBH accretes gas and radiates Lyman-Werner photons, which suppress molecular hydrogen cooling and force the surrounding gas to cool through atomic hydrogen. When the gas is also moving with a streaming velocity of at least $0.8\\sigma_{b\\chi}$ relative to dark matter, dense atomically-cooling clouds form in the PBH's wake, collapse on scales of about 10 pc, and grow by sustained infall at $\\gtrsim0.01-0.1\\,M_\\odot\\,\\mathrm{yr}^{-1}$. The authors identify these clouds as sites where $\\sim10^5\\,M_\\odot$ direct-collapse black hole seeds can form at $z\\sim20-10$, leaving a binary with mass ratio $q\\sim0.1$ around the PBH. If correct, this gives a natural in situ origin for early SMBH binaries and yields concrete targets for JWST, ALMA, LISA, and TianQin.","feed_headline":"Primordial black holes can seed supermassive black hole binaries","feed_subtitle":"Simulations show a primordial black hole's wake can collapse gas into a second black hole, forming a pair JWST and LISA could detect.","key_machinery":"The mechanism is wake-collapse around an accreting PBH. The load-bearing pieces are the Lyman-Werner intensity fitting formula (Eq. 4), derived from ADAF and thin-disk spectral models, which dissociates H$_2$ and H$^-$ and shifts cooling to atomic hydrogen; the Bondi-Hoyle accretion scheme with a radiative-efficiency interpolation between advection-dominated and thin-disk regimes; the baryon-dark matter streaming velocity offset that displaces the gas center of mass from the PBH; and a sink-particle criterion ($n_H\\gtrsim10^6\\,\\mathrm{cm}^{-3}$, $t_{\\mathrm{survive}}\\gtrsim t_{\\mathrm{ff}}$) that tags collapsing cloud gas. The physical identity carrying the argument is the Jeans-rate scaling $\\dot{M}_{\\mathrm{infall}}\\sim c_s^3/G$, which turns atomic-cooling temperatures ($\\sim5000-10^4$ K) into super-Eddington inflow rates of $0.01-0.1\\,M_\\odot\\,\\mathrm{yr}^{-1}$, above the threshold for forming direct-collapse black hole seeds.","core_discovery":"The central claim is that a PBH-DCBH pair naturally forms when three conditions coincide: a $\\sim10^6\\,M_\\odot$ PBH accreting at a few percent of Eddington, Lyman-Werner radiation from its accretion flow that keeps H$_2$ abundance low, and baryon-dark matter streaming of $v_{b\\chi}\\gtrsim0.8\\sigma_{b\\chi}$. Under those conditions, the simulations produce a gravitationally unstable, atomically-cooling cloud of $\\sim10^5\\,M_\\odot$ in the PBH wake, offset from the PBH by roughly 10 pc, with inflow rates high enough to feed a bloated supermassive star that collapses into a direct-collapse black hole. The resulting systems start with mass ratio $q\\sim O(0.1)$ and separation of about 10 pc at redshift $z\\sim20-10$. The paper does not simulate the final seed collapse; it establishes the conditions and argues, from infall rates and cloud masses, that DCBH formation is the expected outcome.","pith_inferences":["One untested corollary is that this channel would also produce off-nuclear massive black holes—secondaries not at galaxy centers—which could show up as offset X-ray sources or as high-redshift electromagnetic counterparts to gravitational waves.","The same wake-collapse logic may extend to other PBH masses: the dimensionless condition is that the PBH's Lyman-Werner bubble suppresses H$_2$ cooling while the streaming wake provides a seed overdensity, so lighter PBHs with faster accretion could be probed with the same machinery.","A statistical prediction follows from the single-seed assumption: if this channel dominates, the secondary seed masses should cluster near $\\sim10^5\\,M_\\odot$ with little scatter, whereas fragmentation would produce a broad low-mass distribution.","The authors' simulation volume represents a tiny effective PBH abundance; folding the per-PBH collapse condition into a full PBH mass function would give a cosmic merger rate that could be compared directly with the LISA stochastic background."],"forward_implications":["SMBH binaries can assemble in situ before galactic mergers, giving a plausible route to the massive black holes JWST sees at $z\\gtrsim7$ without requiring heavy seeds from other channels.","The predicted $q\\sim0.1$, $\\sim10$ pc binaries are direct search templates for ALMA and for dual-AGN signatures in Little Red Dots.","Rapid accretion onto the secondary after formation can drive the mass ratio toward unity, so these systems may end up as strong millihertz gravitational-wave sources for LISA and TianQin.","The collapse outcome is stochastic in streaming velocity, implying a formation rate set by the high-velocity tail of the streaming distribution rather than by the mean.","Even if fragmentation prevents a single DCBH, the resulting dense star cluster would still produce observable extreme-mass-ratio inspirals, tidal disruption events, and intermediate-mass black hole mergers."],"supporting_citations":[{"why":"supplies the gizmo meshless-finite-mass hydrodynamics solver used for all simulations in this paper.","marker":"P. F. Hopkins 2015"},{"why":"provides the PBH initial-condition generator, the baseline CDM and PBH runs, and the streaming-velocity prescription this paper extends to Lyman-Werner feedback.","marker":"S. Zhang et al. 2025"},{"why":"provides the ADAF and thin-disk spectra integrated to obtain the Lyman-Werner intensity fitting formula in Eq. (4).","marker":"V. Takhistov et al. 2022"},{"why":"supplies the PBH accretion-feedback framework and prior evidence for PBH-seeded early structure formation.","marker":"B. Liu et al. 2022"},{"why":"provides the BH accretion kernel model that sets how gas is captured around the PBH.","marker":"M. Tremmel et al. 2017"},{"why":"supplies the radiative-efficiency interpolation between advection-dominated and thin-disk accretion used in feedback energy injection.","marker":"A. Negri & M. Volonteri 2017"},{"why":"sets the 'zone of no return' density threshold for inescapable collapse that motivates the sink-particle criterion.","marker":"K. Inayoshi et al. 2014"},{"why":"provides the local H2 self-shielding factor used to compute photodissociation rates from the Lyman-Werner field.","marker":"J. Wolcott-Green et al. 2011"}],"fun_headline_variants":["PBH wakes collapse gas into companion black holes","Primordial black holes catalyze early SMBH binaries","Supermassive binaries born from PBH accretion wakes","PBH seed black hole pairs in the early universe","Early black hole binaries from PBH-catalyzed collapse"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The channel assumes the infalling gas feeds a single central protostar, because the simulations do not resolve the sub-parsec turbulence, radiative transfer, and angular-momentum transport that decide whether the cloud forms one supermassive star or fragments into many cores.","fun_headline_variants_meta":{"raw":{"variants":["PBH wakes collapse gas into companion black holes","Primordial black holes catalyze early SMBH binaries","Supermassive binaries born from PBH accretion wakes","PBH seed black hole pairs in the early universe","Early black hole binaries from PBH-catalyzed collapse"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000261,"raw_usage":{"total_tokens":1694,"prompt_tokens":1150,"completion_tokens":544,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":766,"completion_tokens_details":{"reasoning_tokens":466}},"tokens_in":766,"tokens_out":544,"duration_ms":7217,"temperature":1.0,"reasoning_tokens":466,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T05:56:42.072990+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A radiation-hydrodynamic simulation that resolves the sub-parsec collapse and finds the atomically-cooling cloud fragments into many low-mass protostars, starving any single protostar of accretion, would remove the DCBH seed and with it the predicted $q\\sim0.1$ binary.","supporting_citations":[{"cited_title":"B., et al","cited_arxiv_id":null,"evidence_quote":"provides the ADAF and thin-disk spectra integrated to obtain the Lyman-Werner intensity fitting formula in Eq. (4)."},{"cited_title":"2017, MNRAS, 467, 3475, doi: 10.1093/mnras/stx362","cited_arxiv_id":null,"evidence_quote":"supplies the radiative-efficiency interpolation between advection-dominated and thin-disk accretion used in feedback energy injection."}],"review_version":1}