{"id":"6aaecaeb-c2b6-4bc9-936b-fec728c57448","arxiv_id":"2607.21956","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"Combining three accretion disk states for binary massive black holes, the paper predicts a spectrum 'notch' and a mass-ratio evolution attractor at q~1e-3 for extreme-mass-ratio binaries.","lead":"Astrophysicists model how pairs of massive black holes accrete gas when one is much smaller than the other, computing the light they emit from infrared to X-rays. They find a universal 'notch' in the spectrum and predict that low-mass-ratio pairs evolve toward a fixed mass ratio of about one-thousandth, while more equal pairs evolve toward equality.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unvalidated interpolation in Eq. (7) sets the bifurcation threshold; test with sensitivity analysis or transitional-regime hydro sims.","rationale":"The paper's most novel and externally relevant result is the evolutionary bifurcation, which predicts a pile-up near q~1e-3 for low-q binaries and equal-mass mergers for high-q systems. The mechanism is purely the sign of (ṁ_s − ṁ_p), and the location of the zero crossings is derived from a fitting function whose central segment is an unvalidated guess. The reader flagged this as the weakest assumption; the paper's own §5.2 confirms the lack of dedicated transitional simulations. A sensitivity analysis is the minimal check that can falsify robustness without requiring new simulations. If the bifurcation boundary is stable under parameter variation, the claim is strengthened; if not, the conclusion must be downgraded to a speculative scenario. The SED predictions are less affected because they are evaluated at fixed q where the accretion modes are better constrained; this is why the concern is load-bearing specifically for the evolution claim, not for the overall framework.","tokens_in":31370,"tokens_out":7175,"duration_ms":67925,"concrete_test":"Perform a sensitivity analysis of Eq. (7): vary the interpolation parameters over conservatively wide ranges (e.g., q1 ∈ [0.005,0.03], q2 ∈ [0.0001,0.001], k ∈ [5,50]) and recompute the evolutionary tracks in Fig. 10 (bottom panel). If the critical initial q above which systems evolve to q→1 shifts by more than a factor of 2 from the nominal ~2.5×10^-3, or if the stable equilibrium at q~1e-3 disappears for any plausible parameter set, the bifurcation claim is not robust. This test uses the authors' own code and requires no new physics.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central evolutionary claim—bifurcation of mass-ratio evolution toward q~1e-3 or q→1—hinges on the sign of (ṁ_s − ṁ_p) in Eq. (10). The sign is governed by the joint fit λ(q) in Eq. (7), which splices a gap-case formula (Li et al. 2023) to a CBD formula (Lai & Muñoz 2023) via a hand-fitted logistic weight with parameters q1=0.012, q2=0.0002, k=20. The unstable fixed point (separatrix) q_crit ≈ 2.5×10^-3—the threshold above which systems evolve to q→1—lies inside the transition region 2×10^-4 < q < 0.012 where this interpolation has no simulation support. The paper itself states in §5.2: 'there are currently no dedicated simulations that concern the mass ratio in transitional state between the gap case and the CBD case.' Because the evolutionary tracks in Fig. 10 are computed with the joint fit across all q, the claim that initial q ≲ a few×10^-3 converge to q~1e-3 depends directly on the assumed shape of λ(q) in this unconstrained interval. A different, equally plausible interpolation could shift or eliminate the separatrix, changing which binaries become equal-mass and which stall at q~1e-3.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models accretion onto binary massive black holes (BHs) across mass ratios 10^-4 ≤ q ≤ 0.5 for a 10^8 M_sun binary with Eddington ratio 0.1. It combines a circumbinary disk and two mini-disks, assigning each disk one of three accretion states (SSD, HAF, or Slim) depending on the local accretion rate, and computes multi-component SEDs. The SED portion finds a characteristic 'notch' feature caused by the gap/cavity, plus distinctive mass-ratio-dependent broadband signatures including HAF X-ray emission and Slim-disk soft X-ray excess. The evolution portion derives dq/dt from the accretion-rate ratio λ(q)=Ṁ_s/Ṁ_p, using a joint fit that splices a gap-case prescription (Li et al. 2023) to a CBD prescription (Lai & Muñoz 2023), and claims that initial q ≲ a few × 10^-3 evolve toward an equilibrium q ~ 10^-3, while larger-q systems evolve toward q → 1 (Eqs. 10–11, Fig. 10).","tokens_in":31825,"tokens_out":2925,"duration_ms":33826,"significance":"If the evolution claim survives scrutiny, it has substantial implications: extreme- and intermediate-mass-ratio binaries would tend to merge near q ~ 10^-3, whereas comparable-mass systems would converge toward equality before merger, affecting LISA event-rate predictions and electromagnetic counterpart searches. The SED predictions are also valuable and more robust: they use standard disk models, show systematic variation with q and separation, and give falsifiable signatures (notch location scaling as a_sep^-0.76, HAF/Slim spectral state changes) that can be tested with UV/optical/X-ray facilities. A clear strength is that the paper is candid about its limitations, including explicitly stating in §5.2 that no dedicated simulations cover the transitional mass-ratio regime. The central weakness is that the claimed bifurcation and its threshold q_crit ≈ 2.5×10^-3 are direct products of the hand-fitted interpolation in Eq. (7), which is not supported by simulations in the relevant interval. This is the load-bearing point for the abstract's headline evolution result.","major_comments":[{"comment":"The evolutionary bifurcation is determined by the sign of (ṁ_s − ṁ_p) in Eq. (10), and the location of the fixed points is set by the joint fit λ_fit(q) in Eq. (7) with parameters q1=0.012, q2=0.0002, k=20. The attracting equilibrium q ≈ 1×10^-3 and the separatrix q_crit ≈ 2.5×10^-3 both lie inside the transition interval 2×10^-4 < q < 1.2×10^-2 where there are no dedicated simulations, as the paper itself states in §5.2: 'there are currently no dedicated simulations that concern the mass ratio in transitional state between the gap case and the CBD case.' The parameters q1, q2, k are hand-picked to make a smooth connection, with no error bars or sensitivity analysis. Because the claimed bifurcation in evolutionary pathways is the central new result, the manuscript needs either a dedicated sensitivity study over plausible alternative interpolations or new transitional-regime hydrodynamic","section":"§2.3.3, Eq. (7); §4, Fig. 10"},{"comment":"The paper's abstract and conclusions state that systems with q ≲ a few × 10^-3 evolve toward q ~ 10^-3, but the same section notes that the mass-growth timescale may exceed the gravitational-wave merger timescale. Concretely, the text states that q0 = 3×10^-3 requires about 1.4×10^8 yr to reach q = 6×10^-3, and Fig. 10 shows t_GW for a_sep = 2000 R_g,p0 and 1000 R_g,p0 can be comparable or shorter; §4 then cautions that 'most binary SMBHs with mass ratios q ≲ 10^-2 are likely to merge before they can evolve into systems with q = 10^-3 or large mass ratios q ∼ 1.' This is an important caveat that should appear in the abstract and conclusions, not only in the discussion. The evolutionary tracks in Fig. 10 need a direct quantitative comparison between the accretion-driven q-evolution timescale and t_GW/t_migration for representative separations, otherwise the headline claim overstates the a","section":"§4, Fig. 10; §5.2"},{"comment":"Both the gap-case formula (Eq. 4–5, from Li et al. 2023) and the CBD-case formula (Eq. 6, from Lai & Muñoz 2023) are calibrated by hydrodynamic simulations with disk aspect ratios h ≳ 0.03, whereas the present work applies them at h = 3×10^-3, the value appropriate for a cold AGN thin disk. The paper acknowledges this in §5.2 ('These simulations usually assume a SSD with an aspect ratio of h ≳ 0.03...'). Since ṁ_s and ṁ_p enter Eq. (10) linearly, this mismatch is not a cosmetic issue for the evolution calculation. The SED results are less sensitive because the qualitative disk-state assignments are robust, but the quantitative λ(q) and therefore the exact equilibrium values are not. A sensitivity analysis varying h (or at least showing the resulting spread in q_crit) would materially improve the paper.","section":"§2.3.1–2.3.3; §5.2"}],"minor_comments":[{"comment":"The abstract says 'q ≲ a few × 10^-3' while the conclusions say 'q ≲ 2.5×10^-3.' These are inconsistent; the quantitative threshold should be stated uniformly.","section":"Abstract and §6"},{"comment":"The functional form λ_fit is an inverse-weighted harmonic mean of λ_gap and λ_CBD. It would help readers to explain why this particular combination (rather than, say, a simple logistic interpolation in log λ) is chosen, and to state explicitly that the parameters are unconstrained.","section":"§2.3.3, Eq. (7)"},{"comment":"The text has several typos, e.g., 'T able' in Table 1, 'magneta' for 'magenta', 'binay' for 'binary', and 'the and dot-dashed' in the Fig. 10 caption. Please proofread carefully.","section":"§2.4"},{"comment":"The scaling ν_notch ∝ a_sep^-0.76 is quoted as a fit; it would be useful to state the expected range of validity (e.g., for which q and a_sep this power-law remains accurate), since the text notes the simple R^-3/4 scaling assumes q ≪ 1.","section":"§3.2, Fig. 6"},{"comment":"The 'Little Red Dots' discussion is interesting but the connection to the paper's binary notch is only qualitative. A brief mention in the conclusions would help emphasize the observational relevance.","section":"§5.1.2"}],"recommendation":"major_revision","confidential_remarks":"The paper has two distinct parts: the SED modeling is solid and worth publishing, but the evolutionary bifurcation claim is currently grounded in an unvalidated interpolation function. The authors have been transparent about this limitation, which is to their credit. My recommendation of major revision rather than rejection reflects the fact that the SED part is independently valuable and the evolution part could be made publishable with a proper sensitivity analysis or transitional-regime simulation support. If the authors cannot provide such support, the evolution claims should be substantially weakened in the abstract and conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this paper gives a plausible one-stop SED model for accreting binary massive black holes across q = 1e-4 to 0.5, combining SSD, hot accretion flow, and slim disk prescriptions for the two mini-disks plus a circumbinary disk. That combination is new relative to the mostly SSD-only binary SED work before it, and the qualitative trends—notch position, soft X-ray excess from a slim secondary at q ~ 1e-4, HAF contribution to X-rays at q ~ 0.1—are physically reasonable and worth having on the record. The authors are also honest: the limitations section explicitly says there are no dedicated simulations for the transitional mass-ratio regime.\n\nThe soft spot is the headline evolution claim. The bifurcation toward q ~ 1e-3 or q → 1 is controlled by the sign of (ṁ_s − ṁ_p), which in turn is controlled by the joint fit λ(q) in Eq. (7). That fit splices a gap-case formula to a CBD formula using a hand-picked logistic with q1 = 0.012, q2 = 0.0002, k = 20, and no error bars or sensitivity analysis. The critical separatix q_crit ≈ 2.5e-3 sits exactly in the unconstrained transition region. A different, equally plausible interpolation could move that threshold or erase the equilibrium entirely. The paper acknowledges this in §5.2, but the abstract and conclusions still present the bifurcation as a main result. So the equilibrium-q prediction is a model-dependent suggestion, not a robust finding.\n\nThat said, I would not call this a fatal flaw. The SED framework is useful even if the evolution part turns out to be wrong, and the authors did the right thing by flagging the missing simulations. The referee should ask them to add a sensitivity analysis around the interpolation parameters and to either soften the evolution claims or identify what observations could discriminate among interpolation choices.\n\nWho gets value from this: people hunting for binary SMBH candidates in UV/optical/X-ray surveys, and anyone modeling the electromagnetic counterparts to massive black hole mergers. I would not cite it for the q-evolution result, but I would cite it for the SED framework and the slim-secondary diagnostic. It should go to peer review—serious referees can push on the interpolation and the right parts are solid enough to warrant their time.","headline":"A useful unified SED framework for accreting binary massive BHs, but the mass-ratio evolution bifurcation is only as solid as an unvalidated interpolation in the gap–CBD transition regime.","tokens_in":755,"tokens_out":1365,"would_cite":true,"duration_ms":33712,"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":"Accreting binary black holes bifurcate: tiny-mass-ratio pairs settle toward q≈10^-3 while larger pairs evolve to equal mass.","keywords":["binary black holes","accretion disks","circumbinary disk","mass ratio evolution","spectral energy distribution","gap/cavity","hot accretion flow","gravitational wave sources"],"falsifier":"Run three-dimensional hydrodynamic simulations of an accreting binary with mass ratios q = 5 × 10^-4, 2 × 10^-3, 5 × 10^-3, and 0.01 in a thin disk (aspect ratio h ≈ 3 × 10^-3) and measure the time-averaged ratio Ṁ_s/Ṁ_p. If it is everywhere above or below one, or crosses once instead of twice, the q ≈ 10^-3 equilibrium and the claimed bifurcation do not hold. Alternatively, measure the mass-ratio distribution of mergers with a space-based gravitational-wave detector: no pile-up near q ≈ 10^-3 would contradict the attractor.","tokens_in":31283,"feed_emoji":"🕳️","tokens_out":5844,"duration_ms":55924,"temperature":0.7,"pith_summary":"The paper builds a self-consistent model of gas accretion onto a binary of massive black holes — an outer circumbinary disk plus two mini-disks, each of which can behave as a thin disk, a hot accretion flow, or a slim disk depending on the local accretion rate — and uses it to predict spectra and long-term evolution across mass ratios q = 10^-4 to 0.5. Its central claim is that the evolution of the mass ratio is set by the difference in Eddington-normalized accretion rates of the two holes: binaries with tiny initial mass ratios (q ≲ a few × 10^-3) are driven to a stable equilibrium at q ≈ 10^-3, while binaries with larger initial mass ratios grow toward equal masses before merger. The same model produces a universal spectral 'notch' from the near-infrared to ultraviolet caused by the cavity or gap the secondary carves in the disk, with signatures that vary systematically with q. A sympathetic reader would care because the predicted bifurcation changes what binaries should look like — both electromagnetically and in gravitational waves — and offers a way to identify candidate systems in existing and upcoming surveys.","feed_headline":"Binary black hole accretion drives two mass-ratio fates","feed_subtitle":"Extreme-ratio pairs settle at q≈10^-3 while larger pairs grow equal, reshaping merger predictions.","key_machinery":"The load-bearing object is the accretion-rate ratio λ(q) = Ṁ_s/Ṁ_p, joined across mass-ratio regimes by a weight-function interpolation (Eq. 7) with endpoints q1 = 0.012, q2 = 0.0002, sharpness k = 20, anchored to the gap-case model (accretion limited by the secondary's gravitational capture region for q ≲ 1.6 × 10^-3) and to the circumbinary-disk formula (≈ 0.5 + 4(1−q)/9) for q ≳ 0.04. This λ(q) feeds the evolutionary equation dq/dt = q(ṁ_s − ṁ_p)/τ_growth, whose two crossings of ṁ_s = ṁ_p define the stable equilibrium at q ≈ 10^-3 and the unstable crossing at q ≈ 2.5 × 10^-3. The same machinery routes the supplied total accretion rate into the two mini-disks, deciding which disk is a thin","core_discovery":"The paper's central claim is that an accreting binary of massive black holes has two evolutionary fates, not one. Writing the growth of the mass ratio q = M_secondary/M_primary as dq/dt = q(ṁ_s − ṁ_p)/τ_growth, where ṁ are accretion rates in Eddington units and τ_growth is the e-folding mass-growth time, the sign of (ṁ_s − ṁ_p) decides everything. The authors construct a smooth joint fit for Ṁ_s/Ṁ_p across the gap regime (q ≲ 1.6 × 10^-3), the circumbinary-disk regime (q ≳ 0.04), and the unmeasured transitional band between them. The fit crosses ṁ_s = ṁ_p twice, at q ≈ 10^-3 and q ≈ 2.5 × 10^-3, making q ≈ 10^-3 an attractor: systems starting below it tend toward q ≈ 10^-3, while systems sta","pith_inferences":["If the q ≈ 10^-3 attractor is real, the mass-ratio distribution of merging massive black holes observed by future space-based detectors should be bimodal, peaking near 10^-3 and 1; a flat distribution would falsify the interpolation.","The paper leaves the transition band 2 × 10^-4 ≲ q ≲ 0.012 without dedicated simulations; a direct simulation suite measuring λ(q) there at AGN-like aspect ratios would settle whether the attractor exists.","The notch mechanism offers a natural discriminator between binary-induced 'V-shaped' spectra and Balmer-break features in high-redshift compact-object populations, since the binary notch has no fixed rest-frame wavelength and should drift with orbital separation.","One can extend the model to higher total accretion rates or lower disk aspect ratios; if λ(q) changes shape there, the equilibrium mass ratio may shift, changing which binaries produce gravitational-wave counterparts."],"forward_implications":["Most extreme- and intermediate-mass-ratio binaries that survive to merger will do so near q ≈ 10^-3 rather than growing to comparable masses.","Binaries with initial q above roughly 2.5 × 10^-3 will tend to become equal-mass before merging, so their gravitational-wave chirp and electromagnetic signatures should be those of near-equal-mass systems.","A broad spectral depression from near-infrared through ultraviolet, with the notch frequency shifting as ν_notch ∝ a_sep^-0.76, is a generic sign of the binary cavity and can be searched for in high-redshift AGN samples.","Secondary-dominated systems with q ≈ 10^-4 can produce excess emission in X-rays, while systems with q ≈ 0.1 produce hot-accretion-flow emission from the primary's disk; broad-band infrared-to-X-ray observations can separate mass-ratio regimes.","Predicted periodic Doppler-boosted flux variations are up to about 20% for low-q systems, complementing the notch as an identification tool."],"fun_headline_variants":["Binary black holes: two mass-ratio fates","Accreting black hole binaries: settle or equalize","Black hole pair accretion drives q to 10^-3 or 1","Binary black hole evolution: two equilibrium paths"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The existence and location of the equilibrium at q ≈ 10^-3 rests on a hand-fitted interpolation of the accretion-rate ratio in the transition region q ≈ 2 × 10^-4 to 0.012, where the paper states no dedicated simulations exist; if the true λ(q) differs there, the attractor could disappear or move.","fun_headline_variants_meta":{"raw":{"variants":["Binary black holes: two mass-ratio fates","Accreting black hole binaries: settle or equalize","Black hole pair accretion drives q to 10^-3 or 1","Binary black hole evolution: two equilibrium paths"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000184,"raw_usage":{"total_tokens":1246,"prompt_tokens":924,"completion_tokens":322,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":668,"completion_tokens_details":{"reasoning_tokens":256}},"tokens_in":668,"tokens_out":322,"duration_ms":3739,"temperature":1.0,"reasoning_tokens":256,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T06:13:02.383273+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run three-dimensional hydrodynamic simulations of an accreting binary with mass ratios q = 5 × 10^-4, 2 × 10^-3, 5 × 10^-3, and 0.01 in a thin disk (aspect ratio h ≈ 3 × 10^-3) and measure the time-averaged ratio Ṁ_s/Ṁ_p. If it is everywhere above or below one, or crosses once instead of twice, the q ≈ 10^-3 equilibrium and the claimed bifurcation do not hold. Alternatively, measure the mass-ratio distribution of mergers with a space-based gravitational-wave detector: no pile-up near q ≈ 10^-3 would contradict the attractor.","supporting_citations":[],"review_version":1}