{"id":"8f9f24b8-a276-4cb9-83ab-f9b42dc530be","arxiv_id":"2607.06067","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":8,"one_line_summary":"An isolated hierarchical triple channel with chemically homogeneous evolution and triple common envelope can produce PISN mass-gap BBH mergers matching GW190706 at ~22% of the observed rate.","lead":"This paper proposes that some black-hole mergers detected by LIGO — where one black hole sits in the 'pair-instability mass gap' (45–130 solar masses) — can form via isolated triple-star systems. A tight inner binary undergoes chemically homogeneous evolution, a tertiary triggers a common envelope that shrinks the orbit, and the inner binary merges to produce a massive second-generation black hole.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The SCATTER triple-CE formalism (Eqs. 4–9) is the load-bearing assumption: post-TCE inner-orbit shrinkage (Eq. 9) sets whether tGW ≤ 1 yr (Eq. 18), and this prescription is unvalidated for triples.","rationale":"The reader correctly identifies the SCATTER formalism as the single most load-bearing assumption. I agree because: (1) SCATTER determines both the GW190706 track viability and the rate estimate — if the post-TCE inner orbit is misestimated, the inner BBH fails to merge promptly and no 2g+1g system forms; (2) the formalism was calibrated on binaries, not triples, and the extension to triples (Eq. 8, partitioning envelope mass via Roche-lobe geometry) is unvalidated against 3D hydrodynamics; (3) the rate calculation (Eq. 21) is linear in f_surv, which is effectively a placeholder for the unknown SCATTER success fraction — setting f_surv = 0.5 without measuring the actual fraction of systems satisfying tGW ≤ 1 yr makes the 22% claim unreliable by potentially an order of magnitude. The paper's other weaknesses (hand-tuned initial conditions for GW190706, neglect of delay-time distribution, q=1 assumption) are secondary: they affect precision but not the existence of the channel. The SCATTER concern is different because it could invalidate the channel entirely if the post-TCE orbit is too large. The paper is a legitimate hypothesis-generating study, and CONDITIONAL with MODERATE confidence is the appropriate verdict. The open-source TSE code provides a path to the Monte Carlo test I propose, but no reproduction artifacts are shipped with the paper. I do not adjust the verdict because the reader already captured the essential concern; my contribution is to specify the mechanism of failure (exponential sensitivity of Eq. 9 to M_interact_in) and the concrete test (Monte Carlo measurement of the tGW ≤ 1 yr success fraction).","tokens_in":22164,"tokens_out":8427,"duration_ms":483936,"concrete_test":"Run a Monte Carlo sampling of the initial parameter distributions (§III.B: M1∈[107,140]M☉, Pin∈[1.2,2.0]d, qin∈[0.9,1.0], Pout∈[400,1200]d, qout∈[0.4,0.5]) through the full TSE pipeline with the SCATTER formalism, and measure the fraction of systems satisfying tGW ≤ 1 yr (Eq. 18). If this fraction is ≪ 0.5, the assumed f_surv overestimates the rate. Additionally, vary the SCATTER parameters η, A, B, δ within their quoted uncertainties and recompute the post-TCE inner orbit for the GW190706 track; if tGW exceeds ~1 Myr for any parameter combination within 1σ, the channel is fragile.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader correctly identifies the SCATTER formalism as the critical weak point. The inner orbit shrinkage (Eq. 9) is exponentially sensitive to M_interact_in / M_binary × F(M1/M2), where M_interact_in (Eq. 8) is the envelope mass that interacts with the inner binary, estimated via Roche-lobe geometry (Eq. 7). The SCATTER formalism was calibrated on post-CE *binaries* (Di Stefano et al. 2023, ref [77]), not triples. The extension to triples assumes the envelope mass distribution between inner and outer orbits follows the same Roche-lobe partitioning, but no 3D hydrodynamic simulation of a triple CE phase has validated this. If the post-TCE inner orbit is larger than predicted by even a factor of ~2–3, tGW (Eq. 18) exceeds the helium-star lifetime of the tertiary (~0.1–1 Myr), the inner binary does not merge before the tertiary collapses, and no 2g+1g BBH system forms. This breaks both the GW190706 track and the rate estimate. The rate (Eq. 21) inherits this uncertainty because f_surv = 0.5 is applied uniformly without measuring what fraction of systems in the sampled parameter space (§III.B) actually satisfy tGW ≤ 1 yr. Since the rate is linear in P_sys × f_surv, and f_surv is effectively a placeholder for the unknown SCATTER success fraction, the 22% claim could be off by an order of magnitude in either direction.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"The manuscript proposes an isolated hierarchical triple stellar evolution channel to produce binary black hole (BBH) mergers with one component in the pair-instability supernova (PISN) mass gap. The channel combines tidally driven chemically homogeneous evolution (CHE) in a tight inner binary with a subsequent triple common envelope (TCE) phase computed via the SCATTER formalism. A specific evolutionary track is shown to reproduce the observed properties of GW190706 (m1≈72.8 M☉, m2≈39.3 M☉, χeff≈0.45). A volumetric merger rate of ~0.011 Gpc⁻³ yr⁻¹ is derived, claimed to account for ~22% of the empirical rate for this sub-population at z≈0.68.","tokens_in":22386,"tokens_out":3572,"duration_ms":201290,"significance":"The question of how black holes in the PISN mass gap form is actively debated, and an isolated triple channel is a legitimate and interesting contribution to this discussion. The paper provides a concrete, step-by-step evolutionary pathway and a falsifiable rate prediction. The combination of CHE (which naturally produces aligned-spin, equal-mass inner binaries with zero recoil) with a TCE-driven inspiral is a physically motivated mechanism. However, the significance of the rate claim is limited by the unvalidated nature of the SCATTER formalism for triples and the unjustified survival fraction, both of which are load-bearing for the quantitative prediction.","major_comments":[{"comment":"§II.B, Eqs. (4)–(9): The SCATTER formalism is the load-bearing physical ingredient of the entire channel. The post-TCE inner-orbit shrinkage (Eq. 9) determines whether t_GW ≤ 1 yr (Eq. 18), which is the criterion for a prompt inner-binary merger. SCATTER was calibrated on post-CE binaries (Di Stefano et al. 2023, ref [77]), not on triple systems. The extension to triples assumes that the envelope mass partitioning between inner and outer orbits follows the Roche-lobe approximation in Eq. (8), but no 3D hydrodynamic validation of this assumption is cited. If the post-TCE inner orbit is larger than predicted by a factor of ~2–3, t_GW will exceed the tertiary helium-star lifetime (~0.1–1 Myr), the inner binary will not merge promptly, and no 2g+1g system forms. The manuscript should explicitly acknowledge this uncertainty and, at minimum, provide a sensitivity analysis showing how the rate和","section":null},{"comment":"§III.C, Eq. (21): The survival fraction f_surv = 0.5 is introduced without justification ('Adopting a fiducial survival fraction of f_surv = 0.5 to account for dynamical disruption'). This parameter enters linearly into the final rate. More critically, it is unclear what f_surv is meant to capture: dynamical disruption, the fraction of systems satisfying t_GW ≤ 1 yr, or both? If it includes the latter, the rate calculation becomes partially circular, because P_sys is already computed from narrow parameter windows (§III.B) chosen because they produce the desired outcome, and f_surv then implicitly absorbs the unknown SCATTER success fraction. The paper should either compute this fraction from the simulations or provide a clear physical justification for the adopted value, with an exploration of how the rate scales with it.","section":null},{"comment":"§III.C, Eq. (21): The quantity ⟨M⟩ in Eq. (21) is never defined. It presumably represents the mean stellar mass for IMF normalization, but its value is not stated. Since the rate is linear in 1/⟨M⟩, this omission makes the rate calculation unreproducible. Please define ⟨M⟩ and state its value.","section":null},{"comment":"§III, §III.C: There is an internal inconsistency in the event count. §III states that after excluding events with negative χ_eff, 'the remaining sources, notably GW230824 and GW190706' are well-explained (implying N≈2). Yet in §III.C, Eq. (22), N_obs = 6 is used to compute the empirical rate R_obs ≈ 0.050 Gpc⁻³ yr⁻¹. Additionally, GW230824 does not appear in the initial sample of six events listed in §III (GW230107, GW230928, GW230820, GW190706, GW190620, GW170729). Please clarify which events are used for the rate baseline and reconcile the filtering described in §III with the N_obs = 6 used in Eq. (22).","section":null}],"minor_comments":[{"comment":"Figure 1 is difficult to parse: the orbital separation and mass labels overlap, and the time axis is non-linear with unmarked jumps (e.g., from 3.644 to 3.772 to 24.609 Myr). A cleaner version with a logarithmic time axis or separate panels would help.","section":null},{"comment":"§III.A: The statement that ω_spin/ω_crit ≈ 1.21 triggers CHE should cite the specific threshold criterion used (e.g., from Li et al. 2025, ref [71]) and note its uncertainty.","section":null},{"comment":"§III.A: The paper states the tertiary helium core 'also experiences violent episodic mass loss via PPISNe before undergoing a FSNe.' The pre-SN mass of the tertiary core and the resulting BH mass (39.3 M☉) should be traced explicitly, as is done for the inner binary components.","section":null},{"comment":"§III.B: The outer period range ΔP_out ∈ [400, 1200] days is stated, but the GW190706 fiducial system has a_out = 3600 R☉, which for ~109 M☉ total mass corresponds to P_out ≈ 1400 days — outside this range. Please reconcile.","section":null},{"comment":"Abstract and §IV: 'GW230107, GW230820, and GW230928' are cited as GWTC-4 events but no reference is given for these beyond the GWTC-4 catalog (ref [12]). If these are preliminary, this should be noted.","section":null},{"comment":"§II.E, Eq. (14): The SFRD parameters a=0.015, b=2.7, c=2.9, d=5.6 are attributed to Madau & Dickinson (2014) and Madau & Fragos (2017) but the values differ between those references. Please clarify which parameter set is used.","section":null},{"comment":"Several typographical issues: 'Merge' in the title should be 'Mergers'; 'A V AILABLE' in the Data Available section; inconsistent spacing in equations.","section":null}],"recommendation":"major_revision","confidential_remarks":"The central physical idea (CHE inner binary + TCE-driven inspiral producing mass-gap BBHs) is sound and interesting. The GW190706 track, while fine-tuned, is acceptable as a proof of concept. The main problem is that the rate prediction (the 22% claim) is not yet reliable: it depends on an unvalidated formalism (SCATTER for triples) and an unjustified placeholder (f_surv). The N_obs inconsistency suggests the empirical baseline was constructed hastily. I would encourage the authors to (a) add a sensitivity analysis on the SCATTER parameters, (b) justify or compute f_surv, and (c) fix the event-count inconsistency. If these are addressed, the paper could be acceptable for publication."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and constructive report. The referee raises four major points: (1) the SCATTER formalism lacks hydrodynamic validation for triples and a sensitivity analysis is needed; (2) the survival fraction f_surv=0.5 is unjustified and its physical meaning unclear; (3) ⟨M⟩ in Eq. (21) is undefined; (4) an internal inconsistency exists in the event count and event list between §III and §III.C. We agree that points (3) and (4) require correction, and that points (1) and (2) warrant additional discussion and sensitivity analysis. We do not claim that the current rate prediction is robust to order-unity uncertainties in the SCATTER formalism, and we will revise the manuscript to make this explicit.","responses":[{"response":"We agree that the SCATTER formalism has not been validated against 3D hydrodynamic simulations of triple common envelopes, and that this is a genuine uncertainty affecting the rate prediction. We will add an explicit caveat in §II.B acknowledging this limitation. Regarding the sensitivity analysis: the referee's concern is well-taken. The critical question is whether the post-TCE inner orbit a_in,f is small enough that t_GW ≤ 1 yr (Eq. 18). In our fiducial track, a_in,f ≈ 13.5 R☉ yields t_GW ≈ 0.05 yr, which is well below the 1 yr threshold. If a_in,f were larger by a factor of 2-3, t_GW would scale as a^4 (Eq. 18), giving t_GW ≈ 0.8-6.4 yr. At the lower end (factor ~2), the system would still satisfy t_GW ≤ 1 yr; at factor ~3, it would not. We will include a sensitivity analysis showing how the rate scales with a multiplicative factor on a_in,f, demonstrating that the channel survives for moderate (~2x) overestimates of the post-TCE inner orbit but fails for larger ones. We will also note that 3D hydrodynamic simulations of TCE phases are needed to resolve this uncertainty definitively.","revision_made":"yes","referee_comment":"SCATTER formalism is load-bearing, calibrated on post-CE binaries not triples, no 3D hydrodynamic validation cited. If post-TCE inner orbit is larger by factor 2-3, t_GW exceeds tertiary helium-star lifetime and no 2g+1g system forms. Should acknowledge uncertainty and provide sensitivity analysis."},{"response":"The referee is correct that f_surv = 0.5 is insufficiently justified and that its physical scope is ambiguous. To clarify: f_surv is intended to account for dynamical disruption of the triple system during the TCE phase and subsequent evolution — i.e., systems that are disrupted by natal kicks, dynamical instabilities, or envelope ejection that unbinds the inner binary. It is not intended to absorb the fraction of systems satisfying t_GW ≤ 1 yr, which is already implicitly encoded in the narrow parameter windows of §III.B. However, we acknowledge that this distinction is not clearly stated in the manuscript, and the referee is right that the current presentation risks circularity. We will revise the text to explicitly define f_surv as the fraction of systems surviving dynamical disruption (excluding the t_GW criterion), and we will add a rate scaling showing R_triple ∝ f_surv for f_surv ∈ [0.1, 1.0]. We note honestly that we cannot compute f_surv from first principles without a full population synthesis of TCE outcomes, which is beyond the scope of this paper. We will adopt f_surv = 0.5 as a fiducial value but present the rate as scaling linearly with this parameter.","revision_made":"yes","referee_comment":"f_surv = 0.5 introduced without justification. Unclear what it captures: dynamical disruption, fraction satisfying t_GW ≤ 1 yr, or both? If it includes the latter, rate calculation becomes partially circular since P_sys is already computed from narrow parameter windows. Should compute from simulations or provide physical justification with exploration of rate scaling."},{"response":"The referee is correct. ⟨M⟩ represents the mean stellar mass obtained by integrating the Kroupa IMF over the standard mass range [0.1, 150] M☉. Using the IMF P(M) ∝ M^{-2.3} for M > 0.5 M☉ (with the full Kroupa piecewise form below that), we compute ⟨M⟩ ≈ 0.35 M☉. We will add this definition and value to the manuscript, and we will show the explicit substitution so that the rate calculation is reproducible. We thank the referee for catching this omission.","revision_made":"yes","referee_comment":"⟨M⟩ in Eq. (21) is never defined. Presumably mean stellar mass for IMF normalization, but value not stated. Rate is linear in 1/⟨M⟩, making calculation unreproducible. Please define and state its value."},{"response":"The referee has identified a genuine inconsistency in our manuscript, and we thank them for catching it. The issue is twofold. First, GW230824 was erroneously mentioned in §III as a well-explained event; it does not appear in our initial sample of six events (GW230107, GW230928, GW230820, GW190706, GW190620, GW170729) and should not have been referenced. This is an error we will correct. Second, there is a logical inconsistency between the filtering described in §III (which excludes negative-χ_eff events, leaving fewer than six) and the use of N_obs = 6 in Eq. (22). The intent was to use all six events as the observational baseline for the empirical rate, while noting that our model specifically explains the subset with positive χ_eff (such as GW190706). However, this logic is not clearly presented. We will revise §III and §III.C to either (a) use N_obs = 6 as the total population baseline and explicitly state that our channel accounts for a subset, or (b) recompute R_obs using only the positive-χ_eff subset. We will adopt option (a) for consistency with the 22% comparison, but will clarify the distinction between the full baseline and the subset our model targets.","revision_made":"yes","referee_comment":"Internal inconsistency in event count. §III says after excluding negative χ_eff events, 'remaining sources, notably GW230824 and GW190706' are well-explained (implying N≈2). But §III.C Eq. (22) uses N_obs=6. Also GW230824 does not appear in the initial sample of six events listed in §III. Clarify which events are used for rate baseline and reconcile filtering with N_obs=6."}],"tokens_in":22232,"tokens_out":1488,"duration_ms":280490,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"Hi —","headline":"A new isolated triple channel for PISN mass-gap BBH mergers — plausible physics, but the rate estimate rests on an unvalidated common-envelope prescription and a hand-tuned example.","tokens_in":23387,"tokens_out":89,"would_cite":false,"duration_ms":98547,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Triple Stars Can Forge Black Holes in the 'Forbidden' Mass Gap","keywords":[],"falsifier":"If future GW catalogs show that PISN mass-gap events with high positive χ_eff are rarer than ~22% of the sub-population rate, or if their spin-orbit misalignments are inconsistent with the strictly coplanar, spin-aligned topology this channel requires, the triple-CHE-TCE pathway would be constrained to a smaller contribution or ruled out as the dominant channel.","tokens_in":22218,"feed_emoji":"⭐","tokens_out":1298,"duration_ms":210765,"temperature":0.7,"pith_summary":"Some LIGO gravitational-wave detections contain black holes weighing 45–130 solar masses — a range where pair-instability supernovae should destroy stars entirely, leaving no remnant. This paper proposes that isolated hierarchical triple stars (three stars orbiting each other in a stable configuration) can bypass this destruction. In the proposed channel, two stars in a very tight inner orbit undergo tidal synchronization that drives chemically homogeneous evolution (CHE), bypassing the giant-phase expansion that normally leads to problematic mass transfer. These stars collapse directly into black holes. When the third, outer star later expands into a giant, it engulfs the inner binary in a triple common envelope (TCE) phase, shrinking the inner orbit enough to force the two black holes to merge promptly. Because the inner black holes are equal-mass and spin-aligned, the merger produces no gravitational recoil kick, so the merged remnant stays bound to the third star's collapsing core. The result is a binary containing one massive second-generation black hole sitting squarely in the PISN mass gap, paired with a smaller first-generation black hole. The paper traces a specific evolutionary track reproducing GW190706 (primary ~72.8 solar masses, secondary ~39.3 solar masses, effective spin ~0.45) and computes a volumetric merger rate of ~0.011 Gpc⁻³ yr⁻¹ at z≈0.68, accounting for roughly 22% of the empirically inferred rate for this sub-population.","feed_headline":"Triple-Star Evolution Builds Black Holes in the 'Forbidden' Mass Gap","feed_subtitle":"A tidy chain of tidal locking, direct collapse, and triple common envelopes reproduces GW190706 and claims 22% of mass-gap merger events.","key_machinery":"Chemically homogeneous evolution (CHE) in the tidally locked inner binary; the triple common envelope (TCE) phase computed via the SCATTER angular-momentum-conservation formalism; the zero-recoil condition from equal-mass, spin-aligned inner mergers; and the direct-collapse (failed supernova) channel for forming first-generation black holes without natal kicks.","core_discovery":"The paper constructs a complete evolutionary chain from a coplanar hierarchical triple of ~109 solar-mass stars at low metallicity (Z=0.001) through CHE of the inner binary, direct collapse via failed supernovae, a TCE phase driven by the tertiary's giant expansion, and prompt inner-binary merger, arriving at a final BBH system whose masses (72.8 + 39.3 solar masses) and effective spin (χ_eff ≈ 0.45) match GW190706. The rate calculation, built from empirically motivated initial-parameter distributions and a triple-star fraction of 0.73, yields a birth probability P_sys ≈ 1.41×10⁻⁹ per massive star, translating to ~22% of the observed mass-gap merger rate.","pith_inferences":["The channel's predictability hinges on the inner binary having a mass ratio of exactly unity after early over-contact equilibration. If real systems retain q < 1, the gravitational recoil from the inner merger would be nonzero, potentially unbinding the system and reducing the rate well below 22%.","The SCATTER formalism is calibrated on post-CE binaries, not on triple systems with the extreme mass ratios and configurations studied here. If its orbital-shrinkage predictions are off by even a modest factor, the inner BBH may not merge promptly (t_GW > 1 yr threshold), breaking the chain.","The assumption of zero natal spin for 1g black holes (from efficient Tayler-Spruit angular momentum transport) is load-bearing for the χ_eff prediction. If 1g BHs retain even modest spins, the effective-spin signature would shift, potentially weakening the match with GW190706.","The delay-time approximation (neglecting delays beyond ~10 Myr) may underestimate or overestimate the rate at z≈0.68 depending on the true cosmic metallicity distribution at that redshift."],"forward_implications":["If this channel is real, a subset of PISN mass-gap BBH mergers should show high positive effective spins (~0.4–0.7) from second-generation remnant spin, distinguishable from cluster-formed hierarchical mergers which would show more isotropic spin orientations.","The channel predicts a correlation between mass-gap primary mass and effective spin: the 2g remnant spin of ~0.686 is inherited from the inner merger, so systems near the center of the gap should cluster around χ_eff ≈ 0.45 when paired with a non-spinning 1g companion.","Future space-based detectors could potentially detect the tertiary companion's gravitational signature or residual eccentricity in pre-merger systems, providing a direct test of the triple-origin hypothesis.","The 22% rate contribution leaves room for complementary channels (cluster dynamics, AGN disks) to produce the remaining mass-gap events, making this a partial rather than exclusive explanation."],"fun_headline_variants":["Triple-star model explains LIGO black holes in the mass gap","Hierarchical triple evolution reproduces GW190706 mass-gap black holes","Tidal locking in triples forms mass-gap binary black holes","Triple common envelopes merge black holes in the PISN mass gap","Triple-star channel bridges the pair-instability supernova mass gap"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The entire channel depends on the SCATTER triple common-envelope formalism, an angular-momentum-conservation prescription with fitted parameters (η, A, B, δ) calibrated on a limited sample of post-CE binaries. If this formalism misestimates the post-CE orbital shrinkage by even a factor of a few, the inner binary either fails to merge promptly or the system is disrupted, invalidating both the rate prediction and the GW190706 evolutionary track.","fun_headline_variants_meta":{"raw":{"variants":["Triple-star model explains LIGO black holes in the mass gap","Hierarchical triple evolution reproduces GW190706 mass-gap black holes","Tidal locking in triples forms mass-gap binary black holes","Triple common envelopes merge black holes in the PISN mass gap","Triple-star channel bridges the pair-instability supernova mass gap","Coplanar triple stars form LIGO's PISN mass gap black holes","Triple stellar evolution predicts 22% of mass-gap black hole mergers","Tidal locking in tight triples builds PISN mass gap black holes"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":1653,"prompt_tokens":637,"completion_tokens":1016,"prompt_tokens_details":null},"tokens_in":637,"tokens_out":1016,"duration_ms":72344,"temperature":1.0,"reasoning_tokens":842,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-08T17:26:38.583913+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If future GW catalogs show that PISN mass-gap events with high positive χ_eff are rarer than ~22% of the sub-population rate, or if their spin-orbit misalignments are inconsistent with the strictly coplanar, spin-aligned topology this channel requires, the triple-CHE-TCE pathway would be constrained to a smaller contribution or ruled out as the dominant channel.","supporting_citations":[],"review_version":1}