{"id":"869787cf-8758-42f8-ac09-4041b2a24173","arxiv_id":"2412.11545","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"Absolute physical parameters of ten W-type contact binaries are derived from simultaneous W-D fits to TMTS light curves and first-time LAMOST radial velocities.","lead":"Ten W UMa contact binaries are modeled with TMTS light curves and LAMOST spectra, yielding masses, radii, and luminosities for both components of each system. The catalog adds absolute parameters for six previously unstudied systems and flags two possible newly formed contact binaries.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Table 3 RV phases appear inconsistent with Table 5 ephemerides (e.g., J0305: 0.172 vs 0.216), so the W-D absolute masses may be systematically biased.","rationale":"The paper's main deliverable is Table 7, and the load-bearing input is the simultaneous W-D fit of photometry and RVs. The reader identified sparse phase coverage as the weak point. My independent check finds a sharper, internal problem: the RV phases printed in Table 3 do not match the ephemerides in Table 5. For J0305, JD 2458410.17413 with T0=2454085.459429 and P=0.246983 gives phase 0.172, not 0.216; for J0132 the mismatch is about 0.04 cycles on one night and about 0.01 on another, so it cannot be a simple zero-point shift. If those phases, or the JD/ephemeris pair used to compute them, entered the W-D fit, K1, K2, and Vγ are systematically displaced, and since M1 and M2 scale roughly with (K1+K2)^3/sin^3 i, the absolute masses and the evolutionary conclusions in Section 5 are not trustworthy. The quoted formal errors do not include this systematic. The proposed check is a direct recomputation and a re-fit; if the re-fit leaves Table 7 unchanged, the inconsistency is only typographic and the central claim stands. Until then, the reader's CONDITIONAL verdict is appropriate, with the condition sharpened to verifying the RV phase bookkeeping and re-running the fits with consistent ephemerides.","tokens_in":28800,"tokens_out":26412,"duration_ms":252383,"concrete_test":"Recompute every phase in Table 3 using the corrected epoch and period in Table 5, phase = frac((JD - T0) / P), and compare with the listed phase. For any target with |Δphase| > 0.005 cycle, re-run the W-D fit with the correctly recomputed phases and the same light curve, and compare q, i, M1, M2, and Vγ against Table 7. If the fitted parameters shift by more than the quoted formal errors, the reported absolute parameters and the Section 5 evolutionary classification are not robust; if they do not shift, the discrepancy is confined to the table and the analysis is unaffected.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that simultaneous W-D fits to TMTS light curves and LAMOST RVs yield reliable absolute masses and radii. A direct check of the published tables undermines this: recomputing the Table 3 phases from the Table 5 corrected epochs and periods gives systematic offsets. For J0132, BJD 2458450.02562 with T0=2458080.255492 and P=0.400938 gives phase 0.263, not the listed 0.30095; for J0305, BJD 2458410.17413 with T0=2454085.459429 and P=0.246983 gives phase 0.172, not the listed 0.216; similar offsets appear for J0047 and J0638. The offsets are not a single constant phase shift (e.g., J0132 differs by about 0.04 cycle on the first night but only about 0.01 on the second), so they cannot be absorbed by redefining the phase zero. If these phases, or the corresponding JDs with an inconsistent ephemeris, were used in the W-D fit, the RV semi-amplitudes K1, K2 and the systemic velocity Vγ are shifted. Because the masses scale roughly as (K1+K2)^3/sin^3 i, the absolute masses in Table 7 and all Section 5 evolutionary statements would be directly affected. The quoted 1-2 km/s RV errors and W-D formal errors do not capture this systematic. At minimum, the phase bookkeeping must be reconciled before the parameters of any of the ten systems are used as ground truth.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes ten W UMa-type contact binaries selected from the TMTS catalog and LAMOST medium-resolution spectra. Using simultaneous Wilson-Devinney fits to TMTS light curves and LAMOST radial velocities, the authors derive orbital parameters, mass ratios, inclinations, fill-out factors, and absolute masses, radii, and luminosities. They model the O'Connell effect with a dark spot, analyze O-C diagrams to obtain period-change rates and possible third-body/Applegate variations for three systems, measure H-alpha emission to identify active systems, and discuss the evolutionary states of the components. Four previously studied targets provide a literature cross-check for the mass ratios.","tokens_in":29291,"tokens_out":14251,"duration_ms":125952,"significance":"If the absolute parameters are reliable, this is a useful addition to the growing sample of contact binaries with spectroscopically determined masses and radii, and it provides the first such measurements for six of the ten systems. The paper deserves credit for using public survey data, providing machine-readable tables, and explicitly acknowledging degeneracies in the spot and third-body/Applegate interpretations. The consistency of the mass ratios with previous studies for J0132, J0305, J1300, and J1402 is a genuine sanity check. However, the central claim depends on radial-velocity phase bookkeeping and on sparse RV sampling, and the O-C analysis uses optimistic error assignments; these issues must be resolved before the absolute parameters can be used as ground truth.","major_comments":[{"comment":"The orbital phases in Table 3 are not reproducible from the corrected ephemerides listed in Table 5, which §4.1 states were used to convert BJD to phase. For example, J0132 at BJD 2458450.02562 with T0 = 2458080.255492 and P = 0.400938 gives (2458450.02562 − 2458080.255492)/0.400938 ≈ 922.26 cycles, i.e., phase ≈ 0.26, not the listed 0.30095; the second J0132 epoch differs by about 0.013 cycles from its listed value. For J0305, BJD 2458410.17413 with T0 = 2454085.459429 and P = 0.246983 gives phase ≈ 0.17, not the listed 0.21550. The offsets are not a constant phase shift, so they cannot be absorbed by redefining phase zero. Because the W-D simultaneous fit uses RV phases to determine K1, K2, and Vgamma, and hence q, a, and the absolute masses in Table 7, the authors must either correct Table 3 or explicitly document which ephemeris and which time columns were actually entered into the fit; otherwise the central parameter table is not reproducible and the quoted internal uncertainties understate the error.","section":"Table 3 vs Table 5; §4.1"},{"comment":"Several targets have very sparse radial-velocity coverage, and the paper does not quantify the resulting systematic error. Table 1 lists only three LAMOST MR exposures for J1300, J1402, and J2236, and the excerpted rows for J0047 contain three RV epochs all between phases 0.71 and 0.83 (Table 3). With no velocities near the opposite quadrature, K1, K2, and Vgamma are not independently pinned down; because the masses scale roughly as (K1+K2)^3/sin^3 i, small phase errors propagate directly into Table 7 and into the evolutionary statements in Section 5.3. Please show the complete RV table, state the number of usable epochs per target, and test robustness, for example by dropping one epoch per target or by fitting with the photometric mass ratio held fixed.","section":"Table 1 and Table 3; §4.1"},{"comment":"The treatment of timing uncertainties is not conservative enough to support the quoted precision. The Figure 3 caption states that errors not given in Table 4 are set to 0.001 d, yet many literature minima come from heterogeneous surveys with unknown systematics; this choice dominates the least-squares weighting and yields formal errors such as 7.00 ± 0.01 × 10^-8 d/yr for J0305 in Table 5. These period-change rates feed into the mass-transfer rates in Table 5, the timescale comparison in Table 8, and the third-body/Applegate discussion in Section 5.2. Please re-fit the O-C diagrams with realistic per-point errors or survey-by-survey weighting, and report the resulting range of dp/dt.","section":"Figure 3 caption; Table 5; §5.2"}],"minor_comments":[{"comment":"The paragraph labeled (iii) for J1300 contains several sentences about J0913 that appear to belong to a separate target description; please move them to the correct system.","section":"§2.2"},{"comment":"The text reports a periodic modulation of 10.14 ± 1.13 days for J1402, but Table 6 lists P3 = 17.9 ± 0.5 yr; please check the units and reconcile with the cited literature.","section":"§2.2(iv)"},{"comment":"For J0132 and J0305, the paper inverts the mass ratio before comparing with previous studies and then suggests that the temperature labels may be swapped; this ambiguity should be stated explicitly in a table note, because swapping T1 and T2 would change the W-type classification and the evolutionary interpretation.","section":"§5.1"},{"comment":"There are several typographical errors that should be corrected in a final pass: 'Chaina' in §2.1, 'exposrue' in Table 1, 'Normarlized flux' in Figure 7, 'photospere' in §4.1, and 'binaires' in §5.2.","section":"Throughout"},{"comment":"The typesetting of Equation (2) is garbled, especially the Irwin formula in the second line; please provide a cleanly formatted version with all defined symbols.","section":"Equation (2)"},{"comment":"The columns NLRS and NMRS are labeled 'exposure times' but the entries are integers that appear to be exposure counts; please rename these columns or clarify the units.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"This is a data-analysis paper whose value depends on careful bookkeeping. The phase/ephemeris inconsistency in Table 3 is the most serious issue; if the authors can show that the W-D fits actually used times and phases consistent with Table 5, or if they regenerate the fits with a corrected ephemeris, the central result may be salvageable. I would also ask the editor to ensure that the full RV table and the complete O-C minimum table are available in machine-readable form in the published version, since the printed excerpt is not sufficient to verify the fits."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read: this is a workmanlike contact-binary paper, fine for the TMTS series. It gives first RV-based absolute parameters for six of ten systems and re-derives four existing ones. The methodology follows the standard W-D + O-C template, and the four re-derived systems agree with the literature, which is a useful sanity check.\n\nThe real problem is in the tables: the RV phases in Table 3 don't match the corrected ephemerides in Table 5. I recomputed, e.g., J0305: BJD 2458410.17413 with T0=2454085.459429, P=0.246983 gives phase 0.172, not 0.21550. J0132 likewise is off by ~0.04 cycle on the first night. The offset isn't constant, so redefining the zero point doesn't fix it. If those phases were used in the fits, the K amplitudes and Vγ could be biased, and the masses scale as (K1+K2)^3, so Table 7 and Section 5 would be affected. If the fit was done with correct phases and the table is just wrong, that's still an error that has to be fixed before anyone trusts the parameters. The paper states the corrected ephemeris was used, but the published phases contradict it. This is the first thing I'd send back.\n\nOther soft spots: the RV phase coverage is genuinely thin for several targets (three points all on one side of the curve). The quoted 1-2 km/s errors won't capture the uncertainty from that sampling. Setting all literature O-C errors to 0.001 d is arbitrary; minor, but they should say what they did. The light curves aren't released, which is a shame for a survey paper. Also minor: the mass ratio comparison in Section 5.1 mixes M1/M2 and M2/M1 conventions without making the conversion explicit for two targets; the conclusion of consistency is probably right, but the presentation is sloppy.\n\nCredit where it's due: the six new systems are a real addition; the O-C work is comprehensive; the 'newly formed contact binary' inference for J0132 and J0913 is properly hedged as a possibility rather than a claim. The three-body/Applegate discussion is appropriately cautious. The self-citations to TMTS-V are fine; they provide the selection catalog.\n\nThis paper deserves a serious referee, but it needs a careful pass on the phase bookkeeping and error treatment before the absolute masses can be cited as ground truth. If those get fixed, it's a solid data paper for contact-binary specialists.","headline":"Solid, standard contact-binary parameter paper whose published RV phases don't match its own ephemerides; fix that before trusting the masses.","tokens_in":29896,"tokens_out":6833,"would_cite":true,"duration_ms":57659,"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 reports absolute masses, radii, and luminosities for ten W UMa contact binaries by fitting minute-cadence light curves and radial velocities together, and identifies two systems as likely newly formed contact binaries.","keywords":["W UMa contact binaries","absolute physical parameters","Wilson-Devinney modeling","radial velocities","O-C period analysis","magnetic activity","stellar evolution","close binary stars"],"falsifier":"Take new medium- or high-resolution spectra for these ten targets at orbital phases covering both velocity extrema; if the resulting semi-amplitudes or systemic velocities differ from the W-D solutions by more than the formal errors, the absolute parameters are biased. A single such re-measurement for J0132, with only a handful of RVs, would already test the claimed mass ratio.","tokens_in":28609,"feed_emoji":"🔭","tokens_out":12401,"duration_ms":103404,"temperature":0.7,"pith_summary":"The paper reports absolute physical parameters—masses, radii, luminosities, mass ratios, inclinations, and fill-out factors—for ten short-period W UMa contact binaries, obtained by fitting minute-cadence photometric light curves and radial-velocity measurements together in a Roche-geometry binary model. All ten systems turn out to be W-type shallow or medium contact binaries, and their unequal light-curve maxima (the O'Connell effect) can be reproduced by adding a dark spot on one component. Eclipse-timing analysis finds long-term orbital period changes in nine systems and cyclic variations in three, while spectral subtraction identifies strong chromospheric activity in four. In eight systems the currently more massive component is a main-sequence star and the less massive component has evolved beyond the main sequence; J0132 and J0913 instead sit near the detached/contact boundary and may be newly formed contact binaries. If these parameters hold, they enlarge the small sample of contact binaries with absolute masses and radii anchored by both photometry and spectroscopy.","feed_headline":"Ten contact binaries gain absolute masses and radii","feed_subtitle":"Simultaneous light-curve and radial-velocity fits mark two systems as newly formed contact binaries.","key_machinery":"The load-bearing machinery is the Wilson–Devinney (W-D) Roche-geometry model, which solves the phase-folded minute-cadence light curve and the measured radial velocities simultaneously. Its adjustable parameters include the orbital semimajor axis $a$, systemic velocity $V_\\gamma$, inclination $i$, mass ratio $q=M_2/M_1$, secondary temperature $T_2$, primary luminosity $L_1$, and common surface potential $\\Omega_1=\\Omega_2$; the fill-out factor is $f=(\\Omega-\\Omega_{\\rm in})/(\\Omega_{\\rm out}-\\Omega_{\\rm in})$. The evolutionary interpretation is carried by the empirical angular momentum–total mass relation $J_{\\rm orb}=1.24\\times10^{52}M_T^{3/5}P^{1/3}q(1+q)^{-2}$, which places J0132 and J0913 close to the detached/contact boundary, and by ZAMS/TAMS tracks from a binary stellar evolution code used to classify component evolutionary states.","core_discovery":"On the paper's own terms, the discovery is that simultaneous modeling of the ten systems' high-cadence light curves and their first measured radial velocities fixes the absolute physical parameters of each binary: component masses, radii, luminosities, mass ratio, inclination, fill-out factor, and separation. All ten are W-type contact binaries with shallow or medium fill-out, and a single dark spot accounts for the observed O'Connell effect. The O-C diagrams show long-term period growth or shrinkage in nine systems and periodic residuals in three that are attributed to possible third bodies or magnetic cycles. The evolutionary states split the sample: eight systems have a main-sequence primary and an evolved secondary, whereas J0132 has both components on the main sequence and J0913 has both above the terminal-age main sequence; from their low fill-out factors and position in the orbital angular momentum–total mass plane, the paper concludes that J0132 and J0913 are probably newly formed contact binaries that recently evolved from detached systems.","pith_inferences":["A natural extension is to apply the same simultaneous-fitting approach to other systems in the minute-cadence sample that currently lack radial velocities; a small set of such calibrators could train a machine-learning surrogate for the model and produce a large catalog of absolute parameters.","The number of radial-velocity epochs per system is very small—sometimes only three—so the quoted internal uncertainties should be treated as lower bounds; adding more spectra at previously unsampled orbital phases is the direct way to test whether the mass ratios are stable.","If J0132 and J0913 are indeed newly formed contact binaries, their light-curve asymmetries and spot configurations might change on year timescales, and their orbital periods might show correlated changes; tracking these systems over the next several years could distinguish spot activity from genuine period evolution."],"forward_implications":["If the fitted parameters are correct, these ten systems can be added to the small set of contact binaries whose absolute masses, radii, and luminosities rest on both photometry and spectroscopy, providing anchors for statistical studies of W-type binaries.","The measured mass ratios for the four previously studied targets agree with earlier photometric-only solutions, so the simultaneous method can be extended to other targets in the same survey sample.","For J0047, J0638, and J1402, the computed mass-transfer timescale is only 1–4% of the thermal timescale, so their decreasing periods are more plausibly driven by angular momentum loss than by conservative mass transfer.","If J0132, J1300, and J1402 host third bodies, J0132's companion would be substellar at about $0.066\\,M_\\odot$, while J1300 and J1402 would have M-dwarf companions; magnetic activity remains a viable alternative, especially for J1402.","The two candidate newly formed contact binaries, J0132 and J0913, should show the relationship between angular momentum loss and initial contact in their future period evolution."],"supporting_citations":[{"why":"Introduces the Roche-geometry light-curve model that the paper uses to solve photometry and RVs simultaneously.","marker":"Wilson & Devinney 1971"},{"why":"Extends the W-D code to the modern differential-correction form used in the fits.","marker":"Wilson 1990"},{"why":"Supplies the square-root limb-darkening coefficients adopted in the light-curve fits.","marker":"van Hamme 1993"},{"why":"Provides the radial-velocity zero-point calibration for the medium-resolution spectra from which the RVs are measured.","marker":"Zhang et al. 2021"},{"why":"Produces the machine-learned catalog of eclipsing binaries from which the ten targets were selected.","marker":"Guo et al. 2024"},{"why":"Supplies the ZAMS/TAMS tracks used to classify the components' evolutionary states.","marker":"Hurley et al. 2002"},{"why":"Gives the orbital angular momentum–total mass relation used to identify J0132 and J0913 as newly formed contact binaries.","marker":"Christopoulou & Papageorgiou 2013"},{"why":"Provides the detached/overcontact boundary line and comparison samples in the angular momentum–mass diagram.","marker":"Eker et al. 2006"},{"why":"Supplies the light-travel-time formula used to fit the cyclic O-C variations attributed to possible third bodies.","marker":"Irwin 1952"}],"fun_headline_variants":["Ten contact binaries get absolute masses and radii from joint data","Two newborn contact binaries identified among ten studied","Absolute parameters for ten contact binaries from two surveys","Joint light and velocity fits reveal ten contact binaries' true masses","High-cadence data give precise masses for ten contact binaries"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results stand or fall on the assumption that the sparse radial-velocity measurements—as few as three epochs, often near one quadrature—pin down the full velocity orbits and hence the component masses.","fun_headline_variants_meta":{"raw":{"variants":["Ten contact binaries get absolute masses and radii from joint data","Two newborn contact binaries identified among ten studied","Absolute parameters for ten contact binaries from two surveys","Joint light and velocity fits reveal ten contact binaries' true masses","High-cadence data give precise masses for ten contact binaries"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00151,"raw_usage":{"total_tokens":6101,"prompt_tokens":1039,"completion_tokens":5062,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":655,"completion_tokens_details":{"reasoning_tokens":4998}},"tokens_in":655,"tokens_out":5062,"duration_ms":34122,"temperature":1.0,"reasoning_tokens":4998,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:50:10.121655+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take new medium- or high-resolution spectra for these ten targets at orbital phases covering both velocity extrema; if the resulting semi-amplitudes or systemic velocities differ from the W-D solutions by more than the formal errors, the absolute parameters are biased. A single such re-measurement for J0132, with only a handful of RVs, would already test the claimed mass ratio.","supporting_citations":[],"review_version":1}