{"id":"e1e64bcf-928e-45e4-ae41-5978fb73bb40","arxiv_id":"2506.15989","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"New photometric solutions and absolute parameters are presented for 84 TESS contact binaries, including 15 candidate merger systems.","lead":"This paper analyzes TESS light curves of 84 totally eclipsing contact binary stars and derives their physical parameters using the PHOEBE code. It identifies 15 low-mass ratio, deep contact binaries that may be on the path to merging into single stars.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Absolute parameters rest on one empirical P–a relation whose scatter is not propagated and whose applicability to this sample is untested; this directly undermines the merger-candidate identification.","rationale":"The reader's weakest assumption correctly identified the reliance on the Li et al. (2021b) period–semi-major-axis relation as the most fragile link. I concur that this is the single most load-bearing concern because all absolute parameters in Table 3 and the subsequent evolutionary and merger-candidate discussion in Section 4 are built on it. The paper's own text flags systematic temperature errors as the reason primaries sit below the ZAMS in the M–L diagram, which reinforces that the absolute scale is not independently verified. The spot-model degeneracy and the lack of code are real but secondary; they affect the interpretation of the O'Connell effect and reproducibility, not the main classification of deep/moderate/shallow contact or the merger-candidate count as directly. The comparison with historical mass ratios for 44 targets provides partial validation of q, but not of the absolute scale. A targeted RV check or scatter-propagation test is the most direct way to settle whether the absolute parameters—and therefore the evolutionary status claims—are reliable. Thus the conditional verdict is appropriate; no verdict change is needed.","tokens_in":19711,"tokens_out":1812,"duration_ms":19622,"concrete_test":"Obtain or simulate RV mass ratios and semi-major axes for at least 10–15 of the 84 targets (e.g., from the literature where available, or new moderate-resolution spectroscopy) and compare the spectroscopically derived a and M_tot to those obtained from Eq. (3). Additionally, recompute Table 3 including the reported scatter in Eq. (3) (approximately 0.02 dex in log a) propagated through Kepler's third law; if the mass uncertainties increase by more than a factor of 3–5, the quoted precision of the absolute parameters is misleading and the evolutionary conclusions in Section 4 require revision.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim depends on converting the 84 photometric solutions into absolute masses, radii, and luminosities via Eq. (3): log a = 0.864 log P + 0.783, fitted to 168 contact binaries in Li et al. (2021b) that had RV and light-curve solutions. The transformation then uses Kepler's third law to extract total mass M = a^3/P^2 (scaled), and individual masses from the fitted mass ratio q. The load-bearing weakness is that Eq. (3) is an average relation with intrinsic scatter, yet no scatter is propagated into Table 3 masses, radii, or luminosities. The quoted MCMC uncertainties in Table 3 (0.02–0.07 M_sun) reflect only the photometric fitting, not the calibration uncertainty. Moreover, Eq. (3) was derived mostly from shorter-period or different metallicity/evolutionary-state systems; if those systems are not representative of the 84 TESS targets, the absolute masses inherit a systematic offset. Since the merger-candidate classification in Section 4 depends on q and f (from photometry) but the evolutionary discussion, J_orb, and initial masses depend on absolute M, R, L, a, a bias in a directly shifts the inferred masses and the positions in the M–L and M–R diagrams. The paper also notes in Section 4 that primary stars fall below ZAMS in the M–L diagram, attributing this to systematic temperature errors—an admission that luminosity is affected by an acknowledged systematic. If the same temperature systematics also affect the adopted T_mean used to set T1, they enter L and hence the evolution conclusions. The 15 merger candidates are selected by q<0.25 and f>50%; their absolute parameters are then used to claim low J_orb as supporting evidence. A concrete check: recompute Table 3 with the scatter in Eq. (3) propagated, and separately test whether the 84 targets follow the fitted relation by checking Gaia parallax-based luminosity against the photometric L or by obtaining a few RV orbits.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents a homogeneous photometric study of 84 totally eclipsing contact binaries selected from the TESS eclipsing binary catalog of Prša et al. (2022). The authors model the TESS light curves with the PHOEBE code and MCMC, obtaining mass ratios, inclinations, fillout factors, temperature ratios, and spot parameters. They classify the sample into 18 deep, 39 moderate, and 27 shallow contact binaries, identify 43 systems showing the O'Connell effect, and single out 15 low-mass ratio deep contact binaries as potential merger candidates. Absolute masses, radii, and luminosities are derived by converting each orbital period to a semi-major axis via an empirical period–semi-major axis relation and applying Kepler's third law. The evolutionary states are discussed using mass–luminosity and mass–radius diagrams, and the photometric solutions are compared with historical studies for 44 targets.","tokens_in":19997,"tokens_out":4585,"duration_ms":49697,"significance":"If the stated results hold, the paper provides a useful, uniform sample of photometric orbital parameters for 84 TESS contact binaries and a candidate list of 15 late-stage merger systems. Strengths include the use of a standard, well-tested modeling code with MCMC uncertainty estimates, a documented and reproducible target-selection procedure, and a broad comparison with previously published solutions, including a re-analysis of two radial velocity curves. The principal weakness is that all absolute parameters—masses, radii, luminosities, orbital angular momenta, and initial masses—rely on a single empirical period–semi-major axis relation whose calibration scatter is not propagated and whose applicability to the TESS sample is not validated. This limits the weight that can be placed on the evolutionary conclusions, although the directly fitted photometric quantities and the merger-candidate classification based on q and f are not affected by this issue.","major_comments":[{"comment":"The absolute masses, radii, and luminosities in Table 3 are obtained by converting each period to a semi-major axis with the relation log a = 0.864 log P + 0.783 from Li et al. (2021b), but the calibration scatter of this relation is not propagated into any of the derived quantities. The quoted uncertainties (e.g., 0.02–0.07 M_sun on M1) reflect only the PHOEBE/MCMC light-curve fitting. Because the evolutionary conclusions in §4 (M-L and M-R positions, initial masses, J_orb) are built on these absolute values, the reader cannot tell how much of the claimed precision is real. Please propagate the slope/intercept uncertainties and the intrinsic scatter of the P-a relation, or at least quantify the resulting systematic uncertainty and add it to the quoted errors.","section":"§4, Eq. (3), Table 3"},{"comment":"The paper states that most primary stars fall below the ZAMS and attributes this to systematic temperature errors that affect luminosity. Since the luminosities also depend on the P-a relation and on the adopted T_mean values, the systematic error budget is not fully characterized. A concrete test, such as varying T_mean by ±100 K and recomputing the positions in the M-L diagram, would show whether the ZAMS offset is robust, or the authors should soften the evolutionary-status statements accordingly.","section":"§4, M-L diagram"},{"comment":"The P-a relation was fitted by Li et al. (2021b), a paper with overlapping authorship, using systems with radial-velocity measurements. Its applicability to the present TESS sample is not tested: the comparison in Table 4 covers only q, i, T2/T1, L2/L1, and f, not absolute parameters, and only two of the 84 targets have re-analyzed radial velocities (TIC 267043786 and TIC 207174531). Please validate the relation on the subset of targets with available RV-based absolute parameters, or explicitly state the assumption and add a systematic caveat to all absolute parameters in Table 3 and to the derived evolutionary discussions.","section":"§4, Eq. (3)"},{"comment":"The division into A-subtype and W-subtype contact binaries is used to interpret the M-L, M-R, and J_orb diagrams, but the operational criterion by which each of the 84 targets is assigned to a subtype is never defined. It is not clear whether the subtype is based on the temperature ratio, the mass ratio, the light-curve morphology, or the relative temperatures of the components. Please state the definition explicitly before drawing evolutionary conclusions from the two subgroups.","section":"§3, Table 2, §4"}],"minor_comments":[{"comment":"Equation (1) contains an extra opening parenthesis and reads awkwardly; also, the radius ratio k is used in the equation before it is defined. Please define k = r2/r1 immediately before the equation and clean up the parentheses.","section":"§3, Eq. (1)"},{"comment":"The sentence 'By examining the 168 contact binaries ... we derived the relationship between the semi-major axis a and the period P' attributes the fit to the present authors, but the relation is taken from Yu et al. (2022) and Li et al. (2021b). Please correct the attribution.","section":"§4, first paragraph"},{"comment":"The phrase 'Phoebe was used to calculate the absolute parameters' is vague: the absolute parameters are computed from the semi-major axis via Kepler's third law and the light-curve-derived q, r1, r2, and T2/T1, not by PHOEBE itself. Please describe the derivation more precisely.","section":"§4, text after Eq. (3)"},{"comment":"The log J_orb column is quoted without uncertainties. Either state explicitly that no error propagation was performed for this quantity, or provide uncertainties that include the P-a calibration scatter.","section":"Table 3"},{"comment":"Some spot temperature entries have clearly erroneous error ranges, e.g., TIC 117978580 lists 0.954^{+0.955}_{-1.058}, which is internally inconsistent. Please check and correct all entries in Table A1.","section":"Table A1"},{"comment":"The selection criterion 'the light curves do not change with time' is described only qualitatively. Please state whether stability was assessed by a quantitative metric (e.g., inspection of time-resolved light curves, scatter statistics) or was purely visual, so the reader can judge the reproducibility of the sample selection.","section":"§2.1, fourth criterion"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's absolute-parameter section relies heavily on the empirical P-a relation of Li et al. (2021b), which shares authorship with this paper. This is not in itself improper, but the relation's scatter is not propagated and its applicability is not demonstrated, so the evolutionary conclusions are substantially weaker than the photometric catalog itself. The referee recommends that the editor ask the authors to either validate the relation with existing RV data or reframe the absolute-parameter discussion as exploratory, with explicit systematic caveats."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a useful, honest sample paper. The authors have extracted PHOEBE+MCMC light-curve solutions for 84 totally eclipsing contact binaries from TESS, 40 of which they say have no prior published study. The catalog itself — q, i, fillout, T2/T1, spot parameters, and the 15 low-mass-ratio deep-contact merger candidates — is the kind of thing workers in this subfield will actually use. I give them credit for the historical comparison: 44 targets have prior studies, they check consistency, and they re-fit two radial-velocity curves that disagreed, which is more diligence than you often see.\n\nThe main soft spot is the absolute-parameter calibration. All masses, radii, and luminosities in Table 3 come from a single empirical P–a relation (log a = 0.864 log P + 0.783, from Li et al. 2021b with overlapping authorship). The scatter in that relation isn't propagated, so the quoted MCMC uncertainties of 0.02–0.07 M_sun are artificially small. The paper even admits that primary stars fall below ZAMS in the M–L diagram due to systematic temperature errors affecting luminosity — yet those luminosities feed the evolutionary-status and initial-mass arguments. This is a genuine weakness, but it is fixable: propagate the relation's scatter, and check a subset against Gaia parallaxes or a few new RVs.\n\nTwo smaller items. The spot model is explicitly degenerate, and the authors acknowledge it, but spot parameters are still listed as results; fine, as long as they're presented as a reference interpretation, which they more or less are. And no code is provided, though machine-readable tables are present. That's not disqualifying but I'd ask for the fitting tables.\n\nIs the central result undermined? Not as much as the stress-test note fears. The 15 merger candidates are selected by q<0.25 and f>50% directly from photometry, not from absolute parameters, so the candidate list probably survives even if the J_orb supporting argument is weakened. The absolute masses and evolutionary claims are the soft part, and the paper would be stronger if the authors said so more plainly.\n\nOverall: a competent, niche sample study that deserves a serious referee. The flaws are addressable, not fatal. I'd send it to review and ask for the P–a scatter propagation and a Gaia cross-check.","headline":"A useful homogeneous catalog of 84 TESS contact binaries, with a real but fixable weakness in how absolute parameters are scaled from a period-separation relation.","tokens_in":20646,"tokens_out":2857,"would_cite":true,"duration_ms":33443,"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 paper derives homogeneous orbital and physical parameters for 84 TESS totally eclipsing contact binaries and identifies 15 low-mass-ratio deep systems as plausible pre-merger contact binaries.","keywords":["contact binaries","eclipsing binaries","TESS","light-curve modeling","mass ratio","merger candidates","O'Connell effect","orbital angular momentum"],"falsifier":"Take double-lined radial-velocity spectra for a subset of the 84 systems, especially the 15 low-mass-ratio merger candidates, and compare the dynamical masses and semi-major axes with Table 3. If the measured masses disagree with the period-derived values by more than the quoted uncertainties in a systematic way, the period-semi-major axis calibration is biased and the absolute parameters, evolutionary states, and merger-candidate status would need revision.","tokens_in":19486,"feed_emoji":"⭐","tokens_out":8659,"duration_ms":94987,"temperature":0.7,"pith_summary":"Using TESS space-telescope light curves, the paper analyzes 84 totally eclipsing contact binary systems whose light curves are stable over time, and it produces a uniform set of orbital and physical parameters for them. A contact binary is a close pair whose two stars overflow their Roche lobes and share a common envelope; a total eclipse means each component passes fully behind the other. The derived parameter set yields 18 deep, 39 moderate, and 27 shallow contact binaries, with 43 systems showing the O'Connell effect (unequal maxima attributed to star-spots). Fifteen systems are deep contact binaries with low mass ratio, and the paper identifies them as candidate contact binaries that may merge into single stars. Absolute masses, radii, and luminosities are then derived through an empirical period-semi-major axis relation, and the systems' evolutionary states are traced on mass-radius and mass-luminosity diagrams.","feed_headline":"15 of 84 contact binaries may be merger-bound","feed_subtitle":"TESS light curves yield uniform masses and radii for the whole sample, flagging late-stage merger candidates.","key_machinery":"The analysis runs on a chain of standard tools: TESS 2-minute PDCSAP light curves are phase-folded and binned to 200 normal points, then fitted with a physics-based eclipsing-binary light-curve model driven by a Markov Chain Monte Carlo sampler, with starting values from a grid search or a genetic algorithm. Asymmetric light curves are handled by adding cool star-spots, with spot parameters included in the fit. The step that turns orbital parameters into physical ones is the empirical relation $\\log a = 0.864\\log P + 0.783$, calibrated on 168 contact binaries with both radial-velocity and light-curve solutions; it converts each orbital period to a semi-major axis, and Kepler's third law then gives masses, radii, and luminosities. Evolutionary status is read from mass-radius and mass-luminosity diagrams with zero-age and terminal-age main-sequence lines, while orbital angular momentum and initial masses come from published formulas.","core_discovery":"The paper's central claim is that a homogeneous re-analysis of 84 TESS totally eclipsing contact binaries gives reliable orbital and absolute parameters, and that 15 low-mass-ratio deep systems among them are plausible contact-binary merger candidates. The 15 candidates satisfy $q<0.25$ and $f>50\\%$, and as a group they sit near the low-angular-momentum edge of the $J_{\\rm orb}$-$M_T$ diagram, which the paper interprets as late-stage evolution toward a single rapidly rotating star. The paper also claims that A-subtype and W-subtype contact binaries follow different evolutionary paths: at a given mass, A-subtype components have larger radii and higher luminosities and lower orbital angular momentum, placing them at a later evolutionary stage. These conclusions follow from photometric solutions alone, with absolute parameters depending on a period-semi-major axis calibration rather than on direct radial-velocity masses.","pith_inferences":["A radial-velocity campaign on even a few of the 84 systems would directly test the period-semi-major axis calibration; without it, the absolute masses rest on a single empirical relation and could carry a systematic bias.","The 15 merger candidates can be checked in later TESS sectors for period decrease, eclipse-timing variations, or a brightening that would signal the onset of coalescence.","The star-spot interpretation of the O'Connell effect is degenerate; multi-band photometry or spectroscopic surface imaging could distinguish cool spots from hot accretion-related regions and alter the derived inclinations.","If the period-semi-major axis relation is internally tied to the same kind of systems, comparison with Gaia parallax-based luminosities and radii would provide an independent check of the absolute parameters."],"forward_implications":["The 84-system parameter set gives a homogeneous sample for studying how contact binaries evolve from shallow to deep contact and eventually merge.","The 15 low-mass-ratio deep systems are concrete targets for follow-up period and eclipse-timing monitoring; if they are pre-merger, their periods should show measurable changes.","The systematic differences between A-subtype and W-subtype components in luminosity, radius, and orbital angular momentum support separate evolutionary tracks for the two subtypes.","The comparison with 44 previously studied targets indicates that mass ratios from machine-learning pipelines need per-target verification, while fillout factors from older studies are generally the least consistent parameter."],"supporting_citations":[{"why":"provides the TESS eclipsing-binary catalog from which the 84 targets were selected.","marker":"Prša et al. (2022)"},{"why":"supplies the 168 contact binaries with radial-velocity and light-curve solutions used to calibrate the period-semi-major axis relation.","marker":"Li et al. (2021b)"},{"why":"is the eclipsing-binary modelling code used for the photometric solutions.","marker":"Prša et al. (2016)"},{"why":"is the machine-learning historical study compared against for 27 of the targets.","marker":"Ding et al. (2023)"},{"why":"provides Gaia DR3 astrometry and temperatures used for cross-matching and primary-star temperatures.","marker":"Gaia Collaboration et al. (2021)"},{"why":"provides the formula used to estimate initial masses of the components.","marker":"Yildiz & Doğan (2013)"},{"why":"provides the ZAMS and TAMS lines for the mass-radius and mass-luminosity diagrams.","marker":"Hurley et al. (2002)"},{"why":"provides the orbital angular momentum formula used in the evolutionary comparison.","marker":"Eker et al. (2006)"}],"fun_headline_variants":["15 contact binaries near merger in TESS sample","TESS survey flags 15 merging contact binaries","Low-mass contact binaries may merge, TESS shows","Two evolutionary paths found for contact binaries","84 contact binaries reanalyzed with TESS light curves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Every absolute mass, radius, and luminosity in this paper comes from a single empirical relation between orbital period and semi-major axis that was calibrated on a different set of contact binaries, and the paper does not test that relation with radial velocities for the 84 targets.","fun_headline_variants_meta":{"raw":{"variants":["15 contact binaries near merger in TESS sample","TESS survey flags 15 merging contact binaries","Low-mass contact binaries may merge, TESS shows","Two evolutionary paths found for contact binaries","84 contact binaries reanalyzed with TESS light curves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000301,"raw_usage":{"total_tokens":1760,"prompt_tokens":993,"completion_tokens":767,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":609,"completion_tokens_details":{"reasoning_tokens":695}},"tokens_in":609,"tokens_out":767,"duration_ms":8888,"temperature":1.0,"reasoning_tokens":695,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T23:44:36.774461+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take double-lined radial-velocity spectra for a subset of the 84 systems, especially the 15 low-mass-ratio merger candidates, and compare the dynamical masses and semi-major axes with Table 3. If the measured masses disagree with the period-derived values by more than the quoted uncertainties in a systematic way, the period-semi-major axis calibration is biased and the absolute parameters, evolutionary states, and merger-candidate status would need revision.","supporting_citations":[],"review_version":1}