{"id":"759938b6-2e1d-47bd-a5e2-4a6f7f4ed015","arxiv_id":"2607.07813","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.5,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Cosmological MBHBs coalesce in ~1 Gyr with high eccentricities; scaling relations from 30 Griffin re-simulations link dynamical-friction, hardening and total times to galaxy and orbital properties.","lead":"High-resolution N-body re-simulations of IllustrisTNG galaxy mergers show massive black hole binaries typically coalesce in about 1 Gyr after highly eccentric encounters. The derived scaling relations supply practical subgrid recipes that bring predicted merger times into line with pulsar-timing-array constraints.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"The ~1 Gyr coalescence peak rests on extrapolating 4-parameter linear fits from N=30 gas-poor major mergers to a larger TNG sample whose selection and resolution differ.","rationale":"The reader correctly isolates the representativeness of the N=30 gas-poor major-merger suite as the weakest link supporting the strongest claim. The short coalescence times are a genuine, PTA-relevant result obtained with realistic cosmological orbits and high-resolution N-body evolution; the scaling relations themselves are carefully constructed and cross-validated. The concern is therefore not that the claim is false, but that its quantitative peak values (0.67 Gyr / 1.0 Gyr) rest on an extrapolation whose stability has not been demonstrated. A leave-one-out or bootstrap re-fit test of the kind proposed above would settle whether that extrapolation is secure. Because the paper already flags the limitation and the overall methodology is sound, the appropriate verdict remains CONDITIONAL; no stronger rejection is warranted.","tokens_in":34119,"tokens_out":718,"duration_ms":8252,"concrete_test":"Re-fit relations (c) and (q) on 25 randomly chosen members of the 30-run sample, apply the new coefficients to the same 174 TNG mergers, and recompute the KDE peaks of Δt_df and t_coal (Fig. 13). If either peak shifts by more than the quoted intrinsic scatter (~0.1–0.15 dex) relative to the published 0.67/1.0 Gyr values, the extrapolated PTA-consistency claim is not robust to the limited training sample.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim (abstract, §4.5, Fig. 13) is that applying relations (c) and (q) to 174 TNG major mergers yields galaxy-merger times peaking at ~0.67 Gyr and total BH coalescence times peaking at ~1.0 Gyr, implying efficient evolution consistent with PTA constraints. Those relations were fitted exclusively to the 30 high-resolution Griffin runs (§2.1.2, Table 3). The 30-run sample is restricted to gas mass <30 % of stellar mass inside 2 r_hm, q_★>1/4, z≤1 and chirp mass >10^8 M_⊙; the 174-run catalogue relaxes only the stellar-mass floor and still excludes gas-rich systems. Relation (q) itself has R^{2}=0.86, Δ=0.11 and LOOCV RMSE=0.37 dex on the training set, and the hardening piece (relation j) drops to R^{2}=0.45 once binary-orbit parameters are unavailable. Because dynamical-friction time is well captured by a_0 but the hardening contribution still carries residual stochasticity in e_b (§4.4.4, §5.2), any systematic offset between the gas-poor training set and the broader TNG population (or between TNG300 resolution and the fitted profiles) shifts the location of the 1 Gyr peak that is advertised as PTA-consistent. The paper acknowledges the N=30 and gas-cut limitations (§5.3) but does not quantify how much the peak would move under a different training subset or under inclusion of the systems that fail the gas cut.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper extracts 30 major, gas-poor galaxy mergers from IllustrisTNG (z≤1, q★>1/4, chirp mass >10^8 M⊙), re-samples them at high resolution with Agama, and evolves them with the Griffin N-body code through dynamical friction and stellar hardening. A semi-analytic Peters integration with MCMC uncertainty on the hardening rate then carries each binary to coalescence. Cosmological encounters and the resulting binaries are found to be highly eccentric (e0 peak ~0.94, eb peak ~0.91, ePTA peak ~0.992). Empirical log-linear scaling relations (Table 3) are derived for dynamical-friction, hardening and total coalescence times via exhaustive subset selection, VIF multicollinearity cuts and LOOCV. Applying the most practical relations (c and q) to a larger TNG sample of 174 mergers yields galaxy-merger times peaking at ~0.67 Gyr and total BH coalescence times peaking at ~1.0 Gyr, presented as consistent with PTA constraints on efficient binary evolution.","tokens_in":34622,"tokens_out":1224,"duration_ms":11866,"significance":"If the short coalescence-time distribution is robust, the work supplies a concrete, observationally motivated subgrid prescription that can replace the common “prompt merger” or purely analytic hardening recipes used in cosmological GWB forecasts. The use of genuine cosmological initial conditions (highly eccentric, DM-stripped secondaries) rather than idealised equal-halo setups is a clear advance over Holley-Bockelmann et al. (2025). Strengths include the careful mass-refinement and mass-dependent softening scheme, the MCMC treatment of s(t), the VIF/LOOCV selection protocol, and the explicit practical guide in §5.4. The eccentricity statistics at PTA entry are also of direct interest for waveform modelling.","major_comments":[{"comment":"The central claim of §4.5 and Fig. 13 (peaks at ~0.67 Gyr and ~1.0 Gyr) rests on extrapolating relations (c) and (q) fitted exclusively to the N=30 gas-poor training set. Relation (q) has R²=0.86 and LOOCV RMSE=0.37 dex; the pure-global hardening relation (j) drops to R²=0.45. The paper acknowledges the gas cut and sample size in §5.3 but does not quantify how the peak location shifts under leave-k-out subsets, bootstrap resampling of the 30 runs, or inclusion of systems that fail the 30 % gas cut. A short sensitivity test is needed before the PTA-consistency statement can be regarded as secure.","section":null},{"comment":"§4.4.4 and Table 3: once binary-orbit parameters are unavailable, the hardening timescale is only weakly constrained (relation j, R²=0.45, Δ=0.38). Because dynamical friction dominates the total time, the coalescence relations remain usable, but the residual stochasticity in eb (explicitly discussed in §5.2) is not folded into the Gaussian mixture of Fig. 13. Propagating that extra scatter (or at least quoting an enlarged uncertainty band) would strengthen the cosmological application.","section":null},{"comment":"§2.1.2 and §5.3: the training sample is restricted to major stellar mergers with gas mass <30 % of stellar mass inside 2 rhm. The larger catalogue of 174 systems still excludes gas-rich mergers. Given that PTA sources may include a non-negligible gas-rich fraction, the claim that the ~1 Gyr peak is representative of the full PTA-relevant population needs either a quantitative bound or a clearer statement of the restricted domain of applicability.","section":null}],"minor_comments":[{"comment":"Fig. 12 panels are dense; adding a one-to-one residual histogram or a colour-coding by TNG volume would help the reader assess residual trends.","section":null},{"comment":"Table 3: the units of each coefficient are not always obvious from the log10 expressions; a short column or footnote listing the physical units of every predictor would improve usability for subgrid implementers.","section":null},{"comment":"§2.2.5: the switch from εBH,df=10 pc to εBH,bin=2 pc after the third pericentre is sensible, but a one-sentence justification that the choice does not affect the measured hardening rate would be welcome.","section":null},{"comment":"Appendix B: the mild offset of the TNG black holes below the Reines & Volonteri relation is noted but not discussed; a brief remark on whether this biases the chirp-mass selection would be useful.","section":null},{"comment":"A few typographical inconsistencies remain (e.g., “parametres” vs “parameters”, occasional missing spaces around ~).","section":null}],"recommendation":"minor_revision","confidential_remarks":"The work is solid and timely for MNRAS. The N=30 limitation is real but the authors are transparent about it; a short sensitivity appendix would convert the paper from “useful” to “authoritative” for subgrid use. No novelty or citation concerns."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The new piece is the set of multi-parameter scalings (Table 3) that map cosmological orbital parameters and host properties onto dynamical-friction, hardening and total coalescence times, derived from 30 high-resolution Griffin re-simulations of IllustrisTNG major mergers. Previous relations (Holley-Bockelmann et al. 2025) used idealised single-galaxy setups with fixed e=0.5; here the encounters are highly eccentric (peak e0~0.94), dark-matter stripped, and drawn from actual cosmological trees. That is a genuine step forward for sub-grid recipes.\n\nWhat they do well is careful. Mass refinement, mass-dependent softening, FMM force control, MCMC hardening-rate fits, VIF cuts and LOOCV are all executed properly. Eccentricity distributions at formation and PTA entry (peak ePTA~0.99, residence ~130 Myr) are cleanly measured. Dynamical friction is tightly predicted by a0 (R2 up to 0.99); total coalescence remains usable even with only global quantities (R2 0.85–0.88). Applying relations (c) and (q) to 174 TNG mergers gives the advertised peaks (~0.67 Gyr galaxy merger, ~1.0 Gyr total BH coalescence). The application step is out-of-sample relative to the fit, so the circularity burden is low.\n\nThe soft spot is exactly the one the stress-test flags, and it is real but not fatal. The 30-run training set is gas-poor, major-merger only, z≤1, chirp mass >10^8 M☉. Relation (q) has R2=0.86 and LOOCV RMSE 0.37 dex; the pure-hardening relation without binary-orbit parameters drops to R2=0.45. Residual stochasticity in eb is acknowledged. The paper is open about N=30 and the gas cut (§5.3) and does not claim universality, but it also does not quantify how far the 1 Gyr peak would move if the training subset changed or gas-rich systems were included. That is a limitation of scope, not a load-bearing flaw in the math or the data that exist.\n\nThis is for people who build PTA background models or write sub-grid BH merger recipes. The methods and the relations themselves are solid enough that a serious editor should send it to referees; the caveats are already written and the central claim holds for the population they studied. I would cite the eccentricity distributions and the usable global-parameter relations. Bring it to reading group if the group cares about PTA or cosmological BH dynamics.","headline":"Solid, usable scaling relations from realistic cosmological ICs; the ~1 Gyr PTA-friendly peak is real for the gas-poor major-merger population they actually trained on, but rests on N=30 linear fits that should not be over-sold as universal.","tokens_in":35162,"tokens_out":664,"would_cite":true,"duration_ms":8423,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Cosmological black-hole binaries typically merge in about 1 Gyr, short enough to match pulsar-timing signals.","keywords":["massive black hole binaries","gravitational wave background","pulsar timing arrays","dynamical friction","stellar hardening","scaling relations","galaxy mergers","IllustrisTNG"],"falsifier":"A larger suite of high-resolution re-simulations that includes gas-rich or minor mergers and finds systematically longer coalescence times (several Gyr or more) would falsify the claim that the typical delay is only ~1 Gyr.","tokens_in":35039,"feed_emoji":"🌌","tokens_out":584,"duration_ms":6018,"temperature":0.7,"pith_summary":"Pulsar timing arrays have detected a nanohertz gravitational-wave background whose amplitude sits a little high compared with standard models of massive black-hole binaries. The missing ingredient is how long those binaries take to merge once their host galaxies collide. This paper takes 30 major galaxy mergers straight from a cosmological simulation, re-simulates them at high resolution, and follows the black holes all the way to coalescence. The encounters and the resulting binaries are almost always highly eccentric. From the runs the authors extract simple scaling relations that link dynamical-friction, hardening and total merger times to host-galaxy and orbital properties. When the relations are applied back to the full cosmological merger catalogue, galaxy mergers finish in roughly 0.7 Gyr and black holes coalesce in roughly 1 Gyr. Those short delays mean binaries evolve efficiently enough to produce the observed gravitational-wave background.","feed_headline":"Black-hole binaries merge in about 1 Gyr","feed_subtitle":"Short delays from cosmological initial conditions match the nanohertz gravitational-wave background","key_machinery":"Empirical power-law scaling relations (Table 3) that express dynamical-friction, hardening and total coalescence times as linear combinations of a few host-galaxy and orbital parameters (chiefly the encounter semi-major axis, eccentricity and central densities or masses).","core_discovery":"Applying the scaling relations derived from high-resolution re-simulations of 30 IllustrisTNG major mergers yields a galaxy-merger time distribution that peaks near 0.67 Gyr and a total black-hole coalescence time that peaks near 1.0 Gyr, implying efficient binary evolution consistent with current pulsar-timing-array constraints.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["MBH binaries coalesce in ~1 Gyr from cosmological starts","Galaxy mergers peak at 0.7 Gyr, BH coalescence at 1 Gyr","Scaling relations yield ~1 Gyr black-hole merger times","Efficient MBHB evolution completes in about 1 Gyr","Cosmological BH binaries merge within ~1 Gyr total"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The 30 carefully selected gas-poor major mergers are representative enough that the fitted linear relations can be extrapolated to the whole cosmological merger population.","fun_headline_variants_meta":{"raw":{"variants":["MBH binaries coalesce in ~1 Gyr from cosmological starts","Galaxy mergers peak at 0.7 Gyr, BH coalescence at 1 Gyr","Scaling relations yield ~1 Gyr black-hole merger times","Efficient MBHB evolution completes in about 1 Gyr","Cosmological BH binaries merge within ~1 Gyr total"]},"model":"grok-4.5","effort":"low","cost_usd":0.004584,"raw_usage":{"total_tokens":1384,"prompt_tokens":834,"num_sources_used":0,"completion_tokens":94,"cost_in_usd_ticks":45840000,"prompt_tokens_details":{"text_tokens":834,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":456,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":834,"tokens_out":94,"duration_ms":5866,"temperature":1.0,"reasoning_tokens":456,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-10T17:34:44.477799+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A larger suite of high-resolution re-simulations that includes gas-rich or minor mergers and finds systematically longer coalescence times (several Gyr or more) would falsify the claim that the typical delay is only ~1 Gyr.","supporting_citations":[],"review_version":1}