{"id":"36f65e49-dcd1-43cf-afe4-cf9178ccd8a1","arxiv_id":"2411.11459","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Combining Gaia, Kepler, and OGLE data shows bimodal mass ratios and eccentric orbits in compact hierarchical triples, tentatively linking them to disk-instability formation and age or metallicity effects.","lead":"The authors combine orbital catalogs from Gaia, Kepler, and OGLE to study how compact triple-star systems form, and they find signs that the outer stars' orbits and mass ratios depend on age and metallicity. The work is an early, approximate comparison that points to follow-up measurements rather than a settled result.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"GAIA/OGLE q3 bimodality may be an artifact of the fixed mbin=2 Msun assumption: q3 roughly scales as mbin^-1/3, so survey-dependent masses can create the two peaks.","rationale":"Good-faith reading: this is a short exploratory proceedings paper that combines existing CHT catalogs and attempts to read formation signatures from global e2 and q3 distributions. It is a useful synthesis, and the authors are transparent that the estimates are approximate. However, the load-bearing step is the construction of q3, because the bimodality and the young/old assignment are the main scientific output. The fixed mbin=2 Msun assumption is not a minor detail: the equation used to solve q3 depends on mbin nonlinearly, and the two peaks are supplied by two catalogs with likely different typical stellar masses. The paper offers no check that the bimodality is robust to varying mbin and no significance test. The reader's conditional verdict already captures this; my stress-test sharpens it by identifying the sign of the bias: under-massive Bulge binaries are pushed to low q3 and more massive GAIA hosts to high q3, precisely the observed split. Thus the central claim is not secure until this is tested. I do not see an internal contradiction in the paper's equations, and the exploration is honest about its limitations, so unconditional rejection would be too harsh; the correct route is to condition acceptance on the per-survey mass test.","tokens_in":4138,"tokens_out":17690,"duration_ms":179024,"concrete_test":"Take the spectroscopically characterized CHT subsample from Moharana et al. (2024) and measure the actual mbin distribution for each survey population (Kepler, OGLE, GAIA-NSS). Recompute the q3 histogram in Fig. 3 using per-survey median mbin instead of a fixed 2 Msun, and run a bootstrap test for whether the low-q3 OGLE peak and high-q3 GAIA peak remain separated at, say, >2 sigma. If the bimodality disappears or the peaks shift between surveys, the central DI+DI/two-population conclusion is an artifact of the assumed mass. A secondary check is to recompute GAIA q3 with the full K = (2*pi*a2*sin(i2)/P2)*(q/(1+q))/sqrt(1-e2^2) for systems with measured i2 and e2.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central evidence is the bimodal q3 distribution in Sec. 4.2/Fig. 3, interpreted as two populations of different ages and as support for sequential disk-instability formation. Every q3 value comes from Eq. 2 with a single assumed inner-binary mass, mbin = 2 Msun (Sec. 3). This is load-bearing because the mapping from mass function to q3 is nonlinear: for low q3, q3 roughly scales as (f/mbin)^(1/3), so a factor-of-two mass difference shifts q3 by about 20-25%. The low-q3 peak is dominated by the OGLE Bulge sample and the high-q3 peak by GAIA/younger systems. If the typical Bulge OGLE inner binary is less massive than 2 Msun, and the typical GAIA host is more massive, the fixed assumption will systematically push OGLE q3 down and GAIA q3 up, creating or strongly enhancing the very bimodality attributed to age. The GAIA values are also approximated by setting e2=0 and K=2*pi*a2/P2 (Eqs. 3-4), so the catalogs are not on the same physical footing. No error bars, bootstrap, or significance test supports the bimodality; Sec. 4.3 explicitly calls the q3 estimates approximate. Since the formation and age conclusions follow directly from the shape and survey decomposition of this distribution, the fixed-mbin assumption is the weakest point on which the central claim rests.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes orbital parameters of compact hierarchical triples (CHTs) from Kepler, OGLE, and Gaia DR3 NSS catalogs. It constructs distributions of tertiary eccentricity e2 and tertiary-to-inner-binary mass ratio q3, arguing that the e2 distribution shows survey-dependent trends attributable to metallicity and that the q3 distribution is bimodal, with a low-q3 peak dominated by old OGLE Bulge systems and a high-q3 peak from younger systems. These demographic features are interpreted as signatures of sequential disk-instability (DI+DI) formation, supplemented by tidal dissipation. The central empirical claims are the survey-decomposed e2 and q3 distributions in Figs. 2 and 3.","tokens_in":4470,"tokens_out":3039,"duration_ms":28973,"significance":"If the claims hold, the paper would provide an observationally grounded link between the formation mechanism of compact hierarchical triples and the demographic differences between Galactic populations, which is a valuable contribution given the recent growth of CHT catalogs. The paper draws on large, publicly available samples (110 Kepler CHTs, 177 OGLE CHTs, and 376 Gaia CHT candidates) and explicitly labels the approximations in the mass-ratio derivation. However, the central demographic conclusions rest on a single assumed inner-binary mass for all systems, on a simplified treatment of Gaia orbits, and on visual comparisons without statistical uncertainties. These limitations directly affect the bimodality and metallicity interpretations, so the significance of the results is conditional on addressing them.","major_comments":[{"comment":"The q3 distribution is constructed by assuming a single inner-binary mass, mbin = 2 Msun, for every system, and for Gaia systems by additionally setting e2 = 0 and approximating K as 2*pi*a2/P2. This is load-bearing because Eq. (2) maps the mass function to q3 nonlinearly; for q3 < 1, q3 scales approximately as (f/mbin)^(1/3), so a factor-of-two difference in the assumed mbin shifts q3 by roughly 20-25%. The low-q3 peak in Fig. 3 is dominated by the OGLE Bulge sample and the high-q3 peak by the Gaia sample, so if typical Bulge OGLE inner binaries are less massive than 2 Msun while typical Gaia hosts are more massive, the fixed assumption would systematically skew OGLE q3 downward and Gaia q3 upward, potentially creating or strongly enhancing the very bimodality that is interpreted as an age effect. The paper should either justify mbin per survey with empirical mass distributions, propagate a range of plausible mbin values, or explicitly show that the bimodality persists under such variations.","section":"Sec. 3, Eqs. (2)-(4)"},{"comment":"The claims of bimodality in q3 and of survey-dependent differences in e2 are based on visual inspection of histograms and cumulative distributions, with no error bars, bootstrap resampling, or significance tests. The sample sizes differ strongly across surveys (e.g., 177 OGLE systems vs. 45 robust Kepler solutions), so Poisson fluctuations and sample-selection effects could produce apparent structure that is not physical. At minimum, the authors should provide per-bin uncertainties and a statistical test (e.g., a bootstrap or a two-sample test) to assess whether the two q3 peaks are distinguishable and whether the e2 sub-flat trends are significant.","section":"Sec. 4.2, Fig. 3 and Sec. 4.1, Fig. 2"},{"comment":"The interpretation of the sub-flat e2 trend at e2 = 0.1-0.25 as indicating metal-poor CHTs, and the interpretation of the q3 bimodality as reflecting old versus young populations, are taken from the authors' own Moharana et al. (2024) and then applied to label the Kepler, OGLE, and Gaia samples without independent metallicity or age estimates for the present sample. This creates a circularity risk: the same calibration is used to interpret the new data and to support the formation scenario. The paper should either obtain or cite direct metallicity and age measurements for the individual systems, or clearly reframe these as tentative hypotheses rather than derived conclusions.","section":"Sec. 4.1 and Sec. 4.2"}],"minor_comments":[{"comment":"The author name \"K.G. He/suppress lminiak\" appears garbled; this should be corrected (presumably to K.G. Hełminiak).","section":"Author list"},{"comment":"The abstract contains received and accepted dates (May 1, 2020; July 28, 2020) that are inconsistent with the November 2024 submission date; these should be removed or updated.","section":"Abstract"},{"comment":"There is a grammatical error in the sentence \"This percentage increases are we go towards high-mass stars,\" which should read \"This percentage increases as we go towards high-mass stars.\"","section":"Introduction"},{"comment":"The statement \"our calculations are approximate and do not precisely estimate the q3\" is important caveat; it should be stated earlier, at the point where the q3 distribution is introduced (Sec. 4.2), so that readers are not misled by the apparent precision of the histograms in Fig. 3.","section":"Sec. 4.3, last paragraph"},{"comment":"The period distribution in Fig. 1 would benefit from a legend or explicit labeling of which histogram corresponds to which survey, as the current gray-scale differentiation is difficult to interpret.","section":"Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"This is a short proceedings-style contribution, and the authors are transparent about the approximate nature of the q3 estimates. The central claim, however, is currently supported only by a single-assumption mass-ratio derivation and visual distribution comparisons. The suggested revisions—propagating mass uncertainties, adding statistical tests, and softening the circular metallicity/age attributions—are within the scope of a revised manuscript and would materially strengthen the paper. I would encourage the editor to request a revision rather than reject, provided the authors can address the mass-ratio sensitivity concern."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is an honest, exploratory conference paper that brings together the GAIA-NSS, OGLE, and Kepler CHT samples and points out survey-to-survey differences in the e2 and q3 distributions. That joint presentation is genuinely new, and the cross-survey eccentricity comparison (Kepler sub-flat at low e2, OGLE intermediate behavior) is a real observation I haven't seen stated before. The OGLE dominance of the low-q3 peak is also new. The paper deserves credit for being transparent: it labels the q3 estimates 'ad-hoc' and says the calculations are approximate.\n\nBut the central interpretation—that the q3 bimodality traces age and the e2 shape traces metallicity—is not actually supported by the data as presented. The load-bearing assumption is a fixed inner-binary mass of 2 Msun for every system. Since q3 scales roughly as (f/mbin)^(1/3), a factor-of-two mass difference between surveys shifts q3 by 20-25%. If the typical OGLE Bulge inner binary is less massive than 2 Msun and the typical GAIA host is more massive, the fixed assumption pushes OGLE q3 down and GAIA q3 up, creating or strongly enhancing the very bimodality attributed to age. The GAIA values are also derived with e2=0 and K=2πa2/P2, so the two catalogs are not on the same physical footing. And the population labels (old vs young, metal-poor vs metal-rich) are imported from the authors' own earlier work and applied to the new samples without independent age or metallicity measurements. Finally, there are no error bars or significance tests anywhere; the bimodality and the sub-flat features are read off histograms. The paper itself admits the q3 estimates are approximate, so it isn't overclaiming, but the main conclusion is a hypothesis, not a measurement.\n\nThis is a conference proceedings, so the length constraints may explain some of the shortcuts. Still, the central claim could be tested cheaply: propagate a plausible range of mbin across surveys and see if the bimodality survives. That should have been done before submitting.\n\nWho is this for? People working on stellar multiplicity and triple formation who want a quick combined view of the current CHT samples. It deserves a serious referee rather than a desk reject, because it makes a concrete, falsifiable claim about formation channels and survey differences, and the underlying data are public. Recommendation: send to peer review, require the authors to test the mbin sensitivity, add at least basic bootstrap uncertainties, and tone down the age/metallicity interpretation unless they bring in external data.","headline":"A useful, honest conference paper that combines the three main CHT catalogs for the first time, but the central age/metallicity interpretation rests on an unexamined fixed-mass assumption and should be treated as a suggestion, not a result.","tokens_in":5000,"tokens_out":2333,"would_cite":true,"duration_ms":25392,"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":"The observed mass-ratio and eccentricity distributions of compact hierarchical triples favor sequential disk-instability formation, with a bimodal mass-ratio split between old and young systems.","keywords":["compact hierarchical triples","triple star formation","disk instability","eclipse timing variations","Gaia non-single stars","mass ratio distribution","eccentricity distribution","Galactic bulge"],"falsifier":"Recompute the q3 distribution using measured inner binary masses (from double-lined eclipsing binaries) and actual Gaia eccentricities for a subsample of CHTs: if the bimodality disappears or the low-q3 peak shifts, the two-population age interpretation is wrong. Alternatively, direct age measurements of the OGLE Bulge CHTs that show they are not old would falsify the age assignment.","tokens_in":3961,"feed_emoji":"⭐","tokens_out":7592,"duration_ms":64814,"temperature":0.7,"pith_summary":"Compact hierarchical triples are three-star systems whose outermost orbit is shorter than 1000 days. This paper combines eclipse-timing samples from Kepler and OGLE with Gaia astrometric orbits to reconstruct, for the first time across surveys, the outer mass-ratio and eccentricity distributions of these systems. It finds that the mass-ratio distribution is bimodal, with a low-q3 peak dominated by old Galactic Bulge systems from OGLE and a high-q3 peak from younger Gaia and Kepler systems. The eccentricity distribution sits above a flat distribution and differs between surveys in a way the authors tie to metallicity. The authors argue both distributions are consistent with sequential disk-instability formation plus later dynamical mass loss, and they identify the specific dynamical mechanisms (tidal dissipation, circumbinary disk interactions) that could drive the mass-ratio shift.","feed_headline":"Bimodal mass ratios in compact triple stars trace age","feed_subtitle":"Kepler, OGLE, and Gaia data show old Bulge triples at low outer mass ratios and young systems near unity.","key_machinery":"The key machinery is the transformation of each survey's orbital parameters into a common tertiary mass ratio q3 = mC/mbin using Kepler's mass function f(m3) = P2 $K^{3}$ / (2 pi G) (1 - $e2^{2}$)^{3/2}, converted into the cubic equation mbin $q3^{3}$ - f $q3^{2}$ - 2f q3 - f = 0 (Eq. 2). For ETV surveys f is measured directly; for Gaia, the paper approximates K = 2 pi a2 / P2 and sets e2 = 0, and for all surveys it assumes mbin = 2 solar masses. This common q3 scale is what allows the distributions to be compared across surveys, and its approximations carry the argument.","core_discovery":"The central claim is that the observed eccentricity and tertiary mass-ratio distributions of compact hierarchical triples are consistent with a single formation channel: sequential disk instability (DI+DI). The mass-ratio distribution is bimodal, with peaks at q3 about 0.2-0.35 and near q3 = 1; the low-q3 peak is populated almost entirely by the OGLE Bulge sample, which the authors interpret as old CHTs, while the high-q3 peak comes from younger GAIA/Kepler systems. The eccentricity distribution is right-shifted relative to flat, and the survey-by-survey differences at low eccentricity are interpreted as a metallicity signal. Taken together, the paper concludes that CHTs form through disk fragmentation in the circumbinary disk and then evolve through mass loss from tidal dissipation or interactions with the circumbinary accretion disk, which moves systems from high to low q3 over gigayears.","pith_inferences":["The mbin = 2 solar masses assumption is strong: if the true inner binary masses in the OGLE sample are systematically lower, the low-q3 peak would shift; measuring masses directly for a subset is a natural next step the paper itself calls for.","The low-q3 OGLE peak may be partly a selection effect of ETV surveys, which are less sensitive to low-mass tertiaries; completeness corrections could change the bimodality's height even if not its existence.","The inferred metallicity-eccentricity connection would be directly testable by measuring spectroscopic metallicities of Gaia CHTs, which currently enter only through the eccentricity-shape comparison.","With future mutual-inclination measurements, the tidal-dissipation versus circumbinary-disk-mass-loss mechanisms for shifting q3 could be distinguished."],"forward_implications":["The low-q3 peak (0.2-0.35) in the combined sample is dominated by OGLE Bulge systems, identifying them as old CHTs, while the near-unity q3 peak from Gaia and Kepler identifies younger systems.","The cumulative e2 distribution lies above a flat distribution at all eccentricities, so CHT outer orbits are dynamically excited beyond what random formation would produce.","Survey-to-survey differences in the low-eccentricity part of the e2 distribution track metallicity, with the Kepler sample apparently metal-poor.","Because stellar evolution cannot provide enough mass loss to move q3 from unity to about 0.2, tidal dissipation or circumbinary disk interactions must be acting, and this also explains the planar CHT population."],"supporting_citations":[{"why":"Supplies the first large Kepler CHT sample with ETV-derived outer orbits, providing the Kepler sub-sample.","marker":"Borkovits et al. (2016)"},{"why":"Supplies the OGLE Galactic Bulge CHT sample and its ETV parameters, the source of the low-q3 peak.","marker":"Hajdu et al. (2019)"},{"why":"Supplies the Gaia DR3 cross-matched CHT orbits used for the high-q3 sample.","marker":"Czavalinga et al. (2023)"},{"why":"Previous spectroscopy that proposed the age-metallicity interpretation and tidal dissipation mechanism this paper builds on.","marker":"Moharana et al. (2024)"},{"why":"Formulates the sequential disk-instability (DI+DI) formation scenario that the distributions are tested against.","marker":"Tokovinin (2017)"},{"why":"Provides the baseline multiplicity statistics and the flat eccentricity expectation for comparison.","marker":"Raghavan et al. (2010)"},{"why":"Offers the circumbinary disk interaction alternative for shifting q3, cited as an alternative explanation.","marker":"Borkovits et al. (2022)"}],"fun_headline_variants":["Bimodal mass ratios in compact triples track age","Old triples favor faint tertiaries, young ones equal","Compact triples form via disk instability alone","Triple star mass ratios betray birth age","Bulge triples skew light, young ones balanced"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole mass-ratio distribution assumes every inner binary has a total mass of 2 solar masses and that Gaia outer orbits can be treated as circular with K from the circular formula; if real binary masses vary between surveys, the two q3 peaks and the identification of the low-q3 peak with old Bulge systems would not be reliable.","fun_headline_variants_meta":{"raw":{"variants":["Bimodal mass ratios in compact triples track age","Old triples favor faint tertiaries, young ones equal","Compact triples form via disk instability alone","Triple star mass ratios betray birth age","Bulge triples skew light, young ones balanced"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000267,"raw_usage":{"total_tokens":1548,"prompt_tokens":809,"completion_tokens":739,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":425,"completion_tokens_details":{"reasoning_tokens":664}},"tokens_in":425,"tokens_out":739,"duration_ms":7150,"temperature":1.0,"reasoning_tokens":664,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:30:14.684521+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the q3 distribution using measured inner binary masses (from double-lined eclipsing binaries) and actual Gaia eccentricities for a subsample of CHTs: if the bimodality disappears or the low-q3 peak shifts, the two-population age interpretation is wrong. Alternatively, direct age measurements of the OGLE Bulge CHTs that show they are not old would falsify the age assignment.","supporting_citations":[],"review_version":1}