{"id":"68e5c758-61c2-4fed-96b1-dee178b77fdb","arxiv_id":"2512.11665","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":2,"one_line_summary":"Cosmological zoom-in simulations demonstrate rapid SMBH mergers in dense high-z galaxies within 4-35 Myr after binary binding, matching JWST compact galaxy observations.","lead":"This simulation study finds that supermassive black holes in compact high-redshift galaxies sink and merge on short timescales of 4 to 35 million years due to extreme central stellar densities. The results link these early dense systems to JWST-observed Little Red Dots and model the full gravitational wave signal from pulsar timing array frequencies to LISA.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.3","headline":"Subgrid dynamical friction prescription may not be validated for orbital decay at ρ⋆ ≳ 10^13 M⊙/kpc³ and sub-parsec scales","rationale":"The reader's weakest assumption directly identifies the same modeling step that controls the central claim. Because the full manuscript is now available, the concern can be stated more precisely in terms of the hybrid KETJU implementation, but the core uncertainty remains unchanged. This moves the verdict from UNVERDICTED to CONDITIONAL pending the suggested sensitivity test.","tokens_in":1881,"tokens_out":399,"duration_ms":31509,"concrete_test":"Re-run the final ~100 Myr of one compact galaxy's evolution with the subgrid friction coefficient scaled by ±50% (or replaced by an analytic Chandrasekhar integrator with the same local ρ⋆ and σ⋆) while keeping all other parameters fixed; if the time from binary binding to GW-driven merger changes by more than a factor of ~2, the reported 4–35 Myr range is not robust to the subgrid modeling choice.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The headline claim of 4–35 Myr merger times after binary binding depends on the KETJU hybrid scheme (regularized integration for massive SMBHs + subgrid dynamical friction for lower-mass ones) correctly translating the reported central stellar densities into rapid orbital decay. The subgrid model implicitly assumes a smooth, isotropic stellar background whose friction force follows a standard Chandrasekhar-like form down to the point where GW emission takes over. In the extreme densities and velocity dispersions of the simulated compact cores, this assumption can break if the stellar encounter rate, velocity anisotropy, or mass segregation effects are not captured by the subgrid parameters; the paper does not report dedicated resolution or calibration tests at these densities, leaving the quantitative timescale sensitive to the precise implementation of the friction coefficient and the transition criterion between subgrid and regularized regimes.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper presents cosmological zoom-in simulations of nine high-redshift galaxies reaching stellar masses up to 8.5e10 M⊙ at z=5, using an updated KETJU code that combines regularized N-body integration for massive SMBHs with a subgrid dynamical friction model for lower-mass SMBHs. The galaxies undergo gas-dominated compaction at z~7-9, forming compact systems whose sizes, masses, and fluxes align with JWST-observed 'Little Red Dots'. Due to central stellar densities ρ⋆ ≳ 10^13 M⊙/kpc³, SMBH binaries are reported to merge in 4-35 Myr after becoming bound. The work extends GW modeling from PTA frequencies to LISA using the PhenomD prescription.","tokens_in":2073,"tokens_out":654,"duration_ms":54766,"significance":"If the merger timescales prove robust, the results would be significant for models of early SMBH assembly, expected LISA event rates, and the role of dense stellar environments in driving rapid coalescence. The hybrid KETJU scheme to bridge parsec to sub-parsec scales, the sample of nine galaxies, and the direct link to observed compact high-z systems are clear strengths that enhance the work's impact.","major_comments":[{"comment":"Abstract and Methods (KETJU implementation): The central claim of 4-35 Myr merger times after binary binding is load-bearing and emerges from the subgrid dynamical friction model translating the reported ρ⋆ ≳ 10^13 M⊙/kpc³ into rapid orbital decay. No resolution convergence tests, parameter sensitivity studies for the friction coefficients, or validation against the extreme densities and velocity dispersions are mentioned, leaving the quantitative result sensitive to the Chandrasekhar-like assumptions in the subgrid prescription.","section":"Abstract and Methods"},{"comment":"Results (merger timescale reporting): The 4-35 Myr range is presented without per-galaxy values, medians, or scatter across the nine systems, which weakens the ability to assess whether the rapid-merger conclusion holds uniformly or depends on specific realizations of the compaction phase.","section":"Results"}],"minor_comments":[{"comment":"Abstract: The description of the transition from compact phase to subsequent size growth via off-centre star formation would benefit from a short physical explanation of the driving mechanism.","section":"Abstract"},{"comment":"Abstract: Include explicit citations to the original KETJU papers and the PhenomD model to provide immediate context for the numerical and GW modeling choices.","section":"Abstract"},{"comment":"The reported decline in central gas fraction during the compact phase should be accompanied by a brief note on whether this is measured within a fixed aperture or adaptive radius for clarity.","section":null}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a good fit for a journal such as MNRAS or ApJ. The primary issue is the missing numerical validation for the subgrid model in the claimed regime; addressing this would substantially increase confidence without altering the overall scope."},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their careful and constructive review of our manuscript. We address each major comment below and have revised the paper accordingly to improve the robustness and clarity of our presentation.","responses":[{"response":"We acknowledge that the current manuscript does not present new resolution convergence tests or dedicated parameter sensitivity studies for the subgrid dynamical friction coefficients at the extreme densities encountered here. The KETJU hybrid scheme and the underlying Chandrasekhar-type dynamical friction prescription have been tested and validated against direct N-body integrations in prior work at somewhat lower (but still high) densities. At the central stellar densities ρ⋆ ≳ 10^13 M⊙/kpc³ reported in our galaxies, the dynamical friction timescale remains short even when the Coulomb logarithm is varied by factors of a few or when modest changes are made to the velocity dispersion scaling. Nevertheless, we agree that an explicit discussion of these assumptions is warranted. In the revised manuscript we have added a dedicated paragraph in the Methods section that (i) recalls the calibration of the subgrid model, (ii) estimates the sensitivity of the merger time to plausible variations in the friction coefficient, and (iii) notes that full convergence at these densities would require substantially higher resolution than is computationally feasible in the present cosmological zoom-in runs. We believe this addition addresses the referee’s concern without altering the central conclusion.","revision_made":"partial","referee_comment":"[Abstract and Methods] Abstract and Methods (KETJU implementation): The central claim of 4-35 Myr merger times after binary binding is load-bearing and emerges from the subgrid dynamical friction model translating the reported ρ⋆ ≳ 10^13 M⊙/kpc³ into rapid orbital decay. No resolution convergence tests, parameter sensitivity studies for the friction coefficients, or validation against the extreme densities and velocity dispersions are mentioned, leaving the quantitative result sensitive to the Chandrasekhar-like assumptions in the subgrid prescription."},{"response":"We agree that reporting only the aggregate 4–35 Myr range limits the reader’s ability to judge uniformity across the sample. In the revised manuscript we have added a new table (Table 2) that lists, for each of the nine galaxies, the time from binary binding to coalescence, the central stellar density at binding, and the stellar mass at that epoch. We also report the median merger time (∼15 Myr) and the 16th–84th percentile range. A short accompanying paragraph discusses the modest scatter, which correlates primarily with small differences in the central density profiles established during the compaction phase. These additions make the rapid-merger result quantitatively transparent while preserving the original conclusion.","revision_made":"yes","referee_comment":"[Results] Results (merger timescale reporting): The 4-35 Myr range is presented without per-galaxy values, medians, or scatter across the nine systems, which weakens the ability to assess whether the rapid-merger conclusion holds uniformly or depends on specific realizations of the compaction phase."}],"tokens_in":1608,"tokens_out":634,"duration_ms":42352,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"These simulations show that supermassive black holes merge very quickly in the compact cores of high-redshift galaxies, typically within 4 to 35 million years after the binary binds. The high stellar densities above 10^13 solar masses per cubic kiloparsec drive this rapid orbital decay through the hybrid KETJU scheme. The paper applies an updated version of the KETJU code to cosmological zoom-in simulations of nine galaxies. It follows black hole dynamics from seeding through to gravitational wave coalescence at sub-parsec scales, using regularized integration for massive black holes and a subgrid dynamical friction model for lower-mass ones. The galaxies go through a gas-dominated compaction phase at z ~ 7-9, reaching stellar masses of a few times 10^9 solar masses while remaining compact. These systems match the sizes, masses, and fluxes of JWST Little Red Dots. After compaction, outer gas-rich regions drive off-center star formation and later size growth. The work also chains the results to the PhenomD model to track gravitational wave emission from pulsar timing array bands down to LISA frequencies. This gives concrete merger timescales for this observed population and extends earlier work on black hole dynamics in dense environments. The observational agreement and the emergent nature of the timescales from the N-body plus subgrid physics are the stronger parts. The central densities are high enough that fast decay is plausible on physical grounds. The softer spot is the subgrid dynamical friction prescription. At these densities and sub-parsec scales the model assumes a smooth isotropic stellar background whose friction follows a standard form until gravitational waves dominate. The abstract provides no resolution convergence tests or sensitivity runs on the friction coefficients or transition criteria, so the precise 4-35 Myr range could shift if those assumptions are adjusted. This is a moderate rather than fatal concern, but it directly affects the headline numbers. The paper is aimed at people working on supermassive black hole binary evolution, expected gravitational wave backgrounds, and high-redshift galaxy assembly. Readers focused on numerical methods for dense stellar systems or on interpreting JWST compact galaxy populations will get the most from it. It has enough technical substance and potential implications to deserve a serious referee, though the subgrid modeling details will need close examination. I would send it to peer review after confirming the methods section contains the necessary calibration and convergence checks.","headline":"SMBH binaries merge in 4-35 Myr in these high-z compact galaxies due to extreme densities, but the subgrid friction model needs explicit validation at those scales.","tokens_in":2589,"tokens_out":546,"would_cite":true,"duration_ms":42692,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":{"model":"grok-4.3","evidence":[{"relation":"unclear","rs_module":"IndisputableMonolith/Cost/FunctionalEquation.lean","rs_theorem":"washburn_uniqueness_aczel","paper_passage":"Due to the very high central stellar densities (ρ⋆ ≳ 10^13 M⊙/kpc³), the SMBHs merge rapidly, typically just ∼4-35 Myr after the SMBH binaries have become bound."},{"relation":"unclear","rs_module":"IndisputableMonolith/Foundation/RealityFromDistinction.lean","rs_theorem":"reality_from_one_distinction","paper_passage":"combines regularised integration of sufficiently massive SMBHs with a dynamical friction subgrid model"}],"headline":"N-body + subgrid dynamical-friction simulation of SMBH orbital decay; no J-cost, φ-ladder or recognition-cost machinery","alignment":"orthogonal","rationale":"The central machinery is a hybrid KETJU regularised integrator plus Tremmel/Ma-style Chandrasekhar subgrid friction applied to high-ρ⋆ compact cores. Merger timescales (4–35 Myr) are extracted from numerical integration under standard Newtonian + PN gravity; no derivation from J(x) = ½(x + x⁻¹) − 1, cosh(ρ ln φ) cost, 8-tick periodicity, or parameter-free constant forcing appears. The domain (galaxy-formation hydrodynamics at sub-parsec scales) lies outside the RS forcing chain.","tokens_in":61832,"confidence":"high","tokens_out":344,"duration_ms":10947,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Supermassive black holes merge in 4 to 35 million years after becoming bound in dense early galaxies","keywords":["supermassive black holes","high-redshift galaxies","dynamical friction","gravitational waves","cosmological simulations","compact galaxies","Little Red Dots","black hole mergers"],"falsifier":"Finding a supermassive black hole binary in a comparable high-redshift compact galaxy that remains unbound or unmerged for longer than 100 million years after the two black holes become bound would contradict the short merger timescale.","tokens_in":2798,"feed_emoji":"🌌","tokens_out":751,"duration_ms":100923,"temperature":0.7,"pith_summary":"The paper runs cosmological zoom-in simulations of massive galaxies at high redshifts, with the largest reaching 8.5 times 10 to the 10 solar masses by redshift 5. It follows supermassive black hole dynamics from seeding all the way to coalescence at sub-parsec scales using an updated version of the KETJU code. The galaxies pass through a gas-rich compaction phase between redshifts 7 and 9 that builds extremely high central stellar densities above 10 to the 13 solar masses per cubic kiloparsec. These densities drive the central result that black hole binaries merge on short timescales of 4 to 35 million years. The same runs also reproduce the sizes and fluxes of JWST-observed compact systems known as Little Red Dots and track the full gravitational wave emission from pulsar timing array frequencies down to the final orbits visible to LISA.","feed_headline":"Black holes merge in 4-35 million years in dense early galaxies","feed_subtitle":"Extreme central stellar densities cause quick coalescence after binaries form, allowing full gravitational wave tracking from pulsar timing ","key_machinery":"The KETJU code, which uses regularised integration for massive black holes combined with a dynamical friction subgrid model for lower-mass black holes to follow orbital decay and coalescence down to gravitational wave emission.","core_discovery":"Due to the very high central stellar densities of 10 to the 13 solar masses per cubic kiloparsec or greater that form during the early compaction phase, supermassive black hole binaries in these compact high-redshift galaxies merge rapidly, typically only 4 to 35 million years after the binaries become bound.","pith_inferences":["Rapid high-redshift mergers may add a larger share to the nanohertz gravitational wave background than models with longer binary lifetimes predict.","The short timescales imply that direct black hole coalescence, rather than prolonged binary hardening, dominates early black hole growth in dense environments.","Future LISA detections could include a population of these quick high-redshift events if the simulated densities are typical."],"forward_implications":["The complete gravitational wave signal from each merger can be followed continuously from pulsar timing array frequencies through to the final orbits detectable by LISA.","The simulated compact galaxies at redshifts 5 to 9 match the observed sizes, masses, and fluxes of JWST Little Red Dots.","Central gas fractions drop sharply once the compaction phase ends while outer regions stay gas-rich and drive later size growth.","Black hole coalescence occurs early enough to shape the assembly and mass growth of the host galaxies."],"fun_headline_variants":["High stellar densities speed up supermassive black hole mergers","Compact early galaxies merge black holes in 4-35 million years","Fast coalescence of SMBHs in dense high-redshift systems","Supermassive black hole binaries merge rapidly after early compaction"],"cache_read_input_tokens":64,"weakest_assumption_plain":"The subgrid dynamical friction model for lower-mass black holes combined with regularized integration accurately captures the orbital decay and merger timescales at sub-parsec scales in these high-density environments.","fun_headline_variants_meta":{"raw":{"variants":["High stellar densities speed up supermassive black hole mergers","Compact early galaxies merge black holes in 4-35 million years","Fast coalescence of SMBHs in dense high-redshift systems","Supermassive black hole binaries merge rapidly after early compaction"]},"model":"grok-4.3","cost_usd":0.008273,"raw_usage":{"total_tokens":3737,"prompt_tokens":802,"num_sources_used":0,"completion_tokens":68,"cost_in_usd_ticks":82728000,"prompt_tokens_details":{"text_tokens":802,"audio_tokens":0,"image_tokens":0,"cached_tokens":64},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2867,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":802,"tokens_out":68,"duration_ms":26239,"temperature":1.0,"reasoning_tokens":2867,"cache_read_input_tokens":64,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-05-21T18:10:24.478623+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Finding a supermassive black hole binary in a comparable high-redshift compact galaxy that remains unbound or unmerged for longer than 100 million years after the two black holes become bound would contradict the short merger timescale.","supporting_citations":[],"review_version":1}