{"id":"5f53b8fd-99b8-4214-9029-707e5bdd5e76","arxiv_id":"2608.02960","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Using CHARA interferometry and radial velocities, the authors measured the visual orbits, dynamical masses (0.5% precision), and distances of two F-type binaries, HD 210763 and HD 221950.","lead":"Astronomers measured the full 3D orbits of two close pairs of stars, obtaining masses precise to 0.5% and model-independent distances. The results give stellar evolution models two clean test points from binaries that behave like single stars.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 0.5% mass and 0.2% distance claims rest on an interferometric model that assumes two unresolved point sources and a single H/K flux ratio; finite stellar diameters or a small band-dependent flux-ratio difference would bias the astrometry below the quoted statistical errors.","rationale":"The paper's core claims are the 0.5% masses and 0.2% distances. The most obvious apparent defects on reading are the HD 210763 distance mismatch (Section 4 gives 94.19 pc, while Table 5 lists 93.62 pc even though the table's alpha = 3.738 mas and a = 0.352 AU imply 94.17 pc) and the HD 221950 age averaging (MIST 3.5 Gyr and BaSTI 4.4 Gyr, combined after adopting [Fe/H] = -0.3 instead of the measured -0.5). I do not make either of these the primary attack: the distance value is likely a typographical error that does not propagate into the mass solution, and the age is a secondary claim whose uncertainty is already visible in the paper's own discussion. The primary risk to the headline '0.5%/0.2%' numbers is the interferometric binary model. The bootstrap error analysis propagates only statistical scatter from the fitted data; it cannot include the model bias from assuming unresolved components and a single H/K flux ratio. The proposed refit with finite angular diameters and band-separated flux ratios is a decisive and inexpensive test: if positions and alpha are stable, the concern is resolved; if not, the quoted uncertainties need a systematic term. This reinforces the reader's conditional verdict without moving it, hence the verdict should remain unchanged.","tokens_in":14811,"tokens_out":14924,"duration_ms":147668,"concrete_test":"Refit the 2024-08-04 epoch of HD 221950 (MJD 60526.4299) with the Section 3 grid-search method twice: first with the current unresolved-point-source and single-flux-ratio model, and second with the SED-based uniform-disk angular diameters fixed (0.18 mas primary, 0.15 mas secondary) and allowing independent f2/f1 for the MIRC-X and MYSTIC data. Then re-run the global orbit fit using the revised relative positions and compare alpha, i, and the derived distance and masses with Table 5. If the recovered position shifts by more than about 0.01 mas, or alpha by more than 0.5%, the quoted uncertainties are underestimated and a systematic term must be added; if the positions agree within the error ellipses, the concern is resolved. Repeat the same check for HD 210763 to confirm the VLTI agreement is not masking a bias.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims inherit their precision from the relative astrometry, and the least secure step is the binary model in Section 3: each system is treated as two unresolved point sources, with one flux ratio f2/f1 fitted jointly to MIRC-X (H-band) and MYSTIC (K-band) data. The paper justifies this by saying the component temperatures are very similar, but the justification is incomplete. The SED-based angular diameters quoted in Section 7 are 0.15 and 0.09 mas for HD 210763 (0.18 and 0.15 mas for HD 221950). At the longest CHARA baselines, with B = 331 m and lambda = 1.65 um, a 0.15 mas uniform disk has visibility V about 0.98, a roughly 2% deficit in the primary's visibility; the model assumes V = 1. A similar-order bias can enter if the true H- and K-band flux ratios differ by even a few tenths of a percent, an effect not included in the bootstrap uncertainties in Table 5. The fitted flux ratios in Tables 3 and 4 scatter by 0.02-0.03 from the H-band-only values (e.g., 0.36 versus 0.38 for HD 210763; 0.68-0.71 versus 0.71 for HD 221950), larger than the approximately 0.003 difference expected from the 100 K Teff difference. If the resulting astrometric bias is comparable to the quoted error ellipses, roughly 0.01-0.02 mas, it translates to a 0.3-0.5% shift in alpha and a corresponding shift in distance. For HD 210763 the independent VLTI orbit provides an external check, but for HD 221950 there is no such check, so the model systematic is the main unresolved risk to the headline precision claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents new CHARA/MIRC-X and MYSTIC interferometric visual orbits for two double-lined spectroscopic binaries, HD 210763 and HD 221950, combined with APO/ARCES and CTIO/CHIRON radial velocities. The authors jointly fit the astrometric and spectroscopic orbits to derive component masses, orbital parallax distances, and full three-dimensional orbital elements. They then disentangle the spectra to estimate effective temperatures and metallicities, fit SEDs to obtain stellar radii and luminosities, and compare the resulting parameters to MIST and BaSTI evolutionary tracks to infer system ages. The headline results are masses of 1.748 and 1.493 solar masses for HD 210763 and 1.098 and 1.031 solar masses for HD 221950, orbital-parallax distances of about 94 pc and 32.75 pc, and ages of about 1.6 Gyr and 3.76 Gyr.","tokens_in":15180,"tokens_out":7027,"duration_ms":65401,"significance":"If the quoted uncertainties are reliable, this is a valuable contribution to the empirical calibration of stellar masses and model-independent distances for non-eclipsing binaries beyond the tidal-circularization period. The HD 210763 orbit is independently confirmed by the VLTI solution of Gallenne et al. (2023), and the orbital-parallax distances agree with Gaia DR3 while improving the HD 221950 distance precision by an order of magnitude. The age analysis provides a useful test of evolutionary-model predictions, especially the end-of-main-sequence position of the HD 210763 primary. The paper is transparent about the use of bootstrapping for orbital uncertainties and about the inconsistencies between MIST and BaSTI ages. However, the central precision claims rest on modeling assumptions—point-source components and a single H/K flux ratio—that are not fully quantified, and at least one internally inconsistent distance value must be corrected before the results can be taken at face value.","major_comments":[{"comment":"The distance to HD 210763 is quoted as 94.19 +/- 0.21 pc in Section 4 but as 93.62 +/- 0.21 pc in Table 5. This is not a rounding difference, and the distance is a headline 0.2%-precision result that is also used as an input to the SED radius fit in Section 5.2 and compared to Gaia DR3 (93.81 +/- 0.26 pc). The paper must adopt a single value, propagate it consistently through the radii and luminosities, and explain the discrepancy.","section":"Section 4 / Table 5"},{"comment":"The astrometric model treats each component as an unresolved point source and fits one flux ratio f2/f1 simultaneously to MIRC-X and MYSTIC data. The angular diameters quoted in Section 7 (0.15 and 0.09 mas for HD 210763; 0.18 and 0.15 mas for HD 221950) imply uniform-disk visibility deficits of roughly 1-4% at lambda = 1.65 microns and the longest CHARA baselines, so the point-source assumption is not obviously negligible at the claimed precision. In addition, the fitted flux ratios in Tables 3 and 4 vary by 0.02-0.03 between nights, larger than the ~0.003 difference expected from the ~100 K temperature difference. These effects are not included in the bootstrap uncertainties of Table 5. I request a quantitative estimate of the induced astrometric bias, for example by fitting with resolved disk diameters or with band-dependent flux ratios, or at least a demonstration that the bias is below the ~0.01 mas level needed to support the 0.5% mass and 0.2% distance claims. This is especially important for HD 221950, which lacks the external VLTI check available for HD 210763.","section":"Section 3 / Tables 3 and 4"},{"comment":"The quoted age of HD 221950, 3.76 +/- 0.38 Gyr, is obtained after abandoning the measured [Fe/H] = -0.5 +/- 0.10 from Section 5.1 and instead computing MIST and BaSTI tracks at [Fe/H] = -0.3. The paper notes this inconsistency but does not propagate the metallicity ambiguity into the age uncertainty. Because the age comparison is a central conclusion, the authors should either quote ages for the full allowed [Fe/H] range or provide a quantitative justification for the -0.3 choice based on abundance-scale systematics. As written, the age error bar reflects only model scatter, not the dominant systematic.","section":"Section 6"}],"minor_comments":[{"comment":"The secondary mass of HD 210763 is 1.492 +/- 0.006 solar masses in the abstract but 1.493 +/- 0.006 solar masses in Section 4 and Table 5; these should be made consistent.","section":"Abstract / Section 4 / Table 5"},{"comment":"Position angles such as 398.9, 407.9, and 416.5 degrees exceed 360 degrees; wrapping these values to the 0-360 degree range or explicitly stating that they are unwrapped would improve clarity.","section":"Table 3"},{"comment":"The caption labels the BaSTI model ages as 1.67 Gyr without overshooting and 1.50 Gyr with overshooting, but the text in Section 6 states the opposite assignment; the caption and text should agree.","section":"Figure 5 caption"},{"comment":"The SED analysis uses the H-band flux ratio from CHARA, but Tables 3 and 4 show night-to-night scatter in f2/f1; please state explicitly which value was adopted and whether its uncertainty was propagated into R1 and R2.","section":"Section 5.2"},{"comment":"The spelling 'BasTI' appears in the Figure 5 caption while the text uses 'BaSTI'; please standardize.","section":"Figure 5 / Section 6"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid, incremental contribution to an established series, and the external validations (VLTI for HD 210763, Gaia for both distances) are genuine strengths. The main blockers are the internal distance inconsistency and the unquantified point-source/flux-ratio modeling systematic, both of which bear directly on the headline precision claims. These should be fixable in revision without changing the scope of the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is Paper V in a reliable series, and it delivers what the series promises—two precise dynamical masses, one first visual orbit (HD 221950), and a distance that beats Gaia by 10x. The HD 210763 orbit reproduces the independent VLTI solution, which is good evidence the pipeline works. If you work on stellar fundamental parameters or tests of evolutionary models, this is useful.\n\nWhat's genuinely new: the first visual orbit of HD 221950 and its component masses (1.098 and 1.031 Msun). The HD 210763 masses (1.748, 1.492 Msun) agree with Gallenne et al. 2023. The age anchor for HD 210763 at 1.6 Gyr looks tight because the primary is at the terminal-age main sequence. The RUWE comparison across the series is a sensible figure.\n\nThe soft spots are not fatal but they need attention. Section 4 gives the HD 210763 distance as 94.19 ± 0.21 pc, while Table 5 lists 93.62 ± 0.21 pc. The table's own alpha and a imply 94.19 pc, so the table entry looks like a typo. Doesn't change the masses, but it's a consistency error a referee will catch. Second, the age for HD 221950 is a mess: the measured [Fe/H] is -0.5, but the models need -0.3, which the authors admit. They then quote a weighted average of MIST (3.5±0.5) and BaSTI (4.4±0.2) as 3.76±0.38; that weighted average is actually 4.3±0.2, so the arithmetic is wrong. The age isn't central to the paper's main claim, but it should be fixed.\n\nThe bigger concern is the interferometric model. The paper assumes two unresolved point sources and a single flux ratio for both H and K bands. The quoted angular diameters are 0.15 mas for the HD 210763 primary and up to 0.18 mas for HD 221950; at the longest CHARA baseline this means visibilities around 0.98, not 1. That is a subtle but real bias in the relative astrometry, and the fitted flux ratios scatter by 2-3% across nights, larger than the few-tenths-percent expected from the 100 K temperature difference. The authors don't propagate this into the quoted 0.5% mass errors. For HD 210763 the agreement with VLTI gives an external check; for HD 221950 there is none, so its headline precision may be optimistic. This is a modeling systematic, not a flaw in the data.\n\nBottom line: this is a legitimate, careful observational paper that deserves a serious referee. It needs a minor revision to fix the table typo, the age arithmetic, and to add a paragraph quantifying the unresolved-source/flux-ratio systematic. I'd send it to review.","headline":"A solid incremental binary-mass paper with one genuinely new orbit; the headline precision needs a caveat about the unresolved-source model and two internal numerical inconsistencies.","tokens_in":15921,"tokens_out":3370,"would_cite":true,"duration_ms":27016,"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":"Combining interferometric and radial-velocity orbits, this paper measures the dynamical masses of HD 210763 and HD 221950 to within 0.5% and their distances to within 0.2%, then uses the resulting stellar parameters to test evolutionary…","keywords":["binaries: spectroscopic","binaries: visual","stars: fundamental parameters","long baseline interferometry","radial velocities","dynamical masses","orbital parallax","stellar evolution models"],"falsifier":"Fit the H- and K-band flux ratios separately to the MIRC-X and MYSTIC data instead of forcing one value: if the resulting flux ratios disagree beyond the fitted uncertainties, or if direct angular-diameter measurements at CHARA's longest baselines show either component to be non-zero in size, the point-source equal-flux-ratio assumption fails and the quoted inclination, distance, and masses would shift.","tokens_in":14530,"feed_emoji":"🔭","tokens_out":5652,"duration_ms":46203,"temperature":0.7,"pith_summary":"Using long-baseline interferometry at the CHARA Array and high-resolution spectroscopy with APO and CTIO, this paper measures the full three-dimensional orbits of two F-type double-lined spectroscopic binaries, HD 210763 and HD 221950. The combined astrometric and radial-velocity fit delivers component masses to 0.5% precision and distances to 0.2%, giving $M_1=1.748\\pm0.008\\,M_\\odot$ and $M_2=1.492\\pm0.006\\,M_\\odot$ for HD 210763, and $M_1=1.098\\pm0.006\\,M_\\odot$ and $M_2=1.031\\pm0.005\\,M_\\odot$ for HD 221950. Temperatures and radii from spectral disentangling and spectral energy distribution analysis allow the authors to compare both systems with MIST and BaSTI evolutionary models, yielding ages of $1.6\\pm0.1$ Gyr and $3.76\\pm0.38$ Gyr. Because these systems have periods well beyond the tidal circularization limit, they serve as relatively clean tests of single-star evolution, free of the tidal distortion and eclipses that complicate shorter-period binaries.","feed_headline":"Orbits nail two binary stars' masses to 0.5%","feed_subtitle":"Merging CHARA astrometry with radial velocities also gives a distance ten times sharper than Gaia's for one system.","key_machinery":"The load-bearing object is the combined visual-plus-spectroscopic orbit solution, obtained through the Schaefer et al. (2016) grid-search procedure: interferometric squared visibilities and closure phases from the MIRC-X and MYSTIC beam combiners are fitted to recover the secondary star's relative position and error ellipse at each epoch, and those astrometric points are then fitted simultaneously with TODCOR-derived radial velocities to solve for the full Keplerian orbital elements. Bootstrapping distributions of those elements convert directly into mass and distance uncertainties. This mechanism turns a handful of milliarcsecond position measurements into sub-percent stellar masses.","core_discovery":"The central discovery is that pairing milliarcsecond astrometry of the secondary star's position with radial-velocity curves yields a complete three-dimensional orbit for each binary, from which component masses follow directly: for HD 210763, $M_1 = 1.748\\pm0.008\\,M_\\odot$ and $M_2 = 1.493\\pm0.006\\,M_\\odot$; for HD 221950, $M_1 = 1.098\\pm0.006\\,M_\\odot$ and $M_2 = 1.031\\pm0.005\\,M_\\odot$. The same fit gives orbital-parallax distances of $94.19\\pm0.21$ pc and $32.75\\pm0.07$ pc for the two targets, the latter improving on the Gaia DR3 parallax by an order of magnitude. With masses fixed, the paper derives effective temperatures near 6400 K and radii ($R_1 = 2.96\\pm0.11\\,R_\\odot$ and $R_2 = 1.81\\pm0.07\\,R_\\odot$ for HD 210763; $R_1 = 1.30\\pm0.08\\,R_\\odot$ and $R_2 = 1.09\\pm0.07\\,R_\\odot$ for HD 221950). Comparing these to evolutionary tracks, the MIST models match both systems at single ages, while BaSTI models succeed only for the lower-mass system, exposing differences in convective core overshooting and abundance scales.","pith_inferences":["Fitting the H- and K-band flux ratios separately rather than forcing a single value would test whether the assumed equal-flux-ratio approximation biases inclination and distance at the quoted precision; for systems with a larger temperature difference between components, this modeling choice probably becomes the dominant systematic.","The mismatch between the SpecMatch-derived iron abundance for HD 221950 and the value the evolutionary models need could reflect an abundance-scale offset; a high-resolution abundance analysis would let the age estimate be recalculated on a consistent metallicity scale.","With even longer baselines at CHARA, direct angular-diameter measurements could reduce radius uncertainties from roughly 6% toward 1%, turning these two systems into full calibrators of stellar model radii rather than mass-only tests."],"forward_implications":["The 0.5%-precision masses join the set of dynamical masses for non-eclipsing binaries, giving stellar evolution models a test that is free of tidal distortion and eclipse geometry.","Orbital-parallax distances independently confirm Gaia's trigonometric parallaxes, and for HD 221950 improve the precision by a factor of ten, providing a check on systematic biases such as the RUWE indicator.","If the MIST models are right, each of the two systems is coeval at a single age, and HD 210763's primary, sitting at the main-sequence turnoff, pins the system age to $\\pm0.1$ Gyr.","The same observational method transfers directly to other long-period spectroscopic binaries: each new orbit yields masses and distances without requiring eclipses, extending precision stellar astrophysics to wider systems."],"supporting_citations":[{"why":"Supplies the grid-search and combined orbit-fitting method used to derive the relative astrometry and full orbital solutions.","marker":"Schaefer et al. 2016"},{"why":"Prior spectroscopic orbit of HD 210763, used for RV templates, comparison, and literature radial velocities.","marker":"Fekel et al. 2011"},{"why":"Prior spectroscopic orbit and metallicity estimate for HD 221950, used for templates and comparison.","marker":"Tomkin & Fekel 2008"},{"why":"First visual orbit of HD 210763 with VLTI, used to check consistency between CHARA and VLTI results.","marker":"Gallenne et al. 2023"},{"why":"Gaia DR3 trigonometric parallaxes that the orbital parallax distances are compared against and confirmed by.","marker":"Gaia Collaboration et al. 2023"},{"why":"The MIST evolutionary models used to estimate system ages and test convective core overshooting.","marker":"Dotter 2016; Choi et al. 2016"},{"why":"The BaSTI-IAC evolutionary models used as an independent comparison for ages and overshooting prescriptions.","marker":"Hidalgo et al. 2018"},{"why":"Nominal solar values adopted when converting orbital elements into component masses and distances.","marker":"Prša et al. 2016"}],"fun_headline_variants":["CHARA orbits yield 0.5% masses for two binaries","Visual orbits and RVs pin masses to 0.5%","Long-period binaries get precise masses and ages","Two binaries weighed to 0.5% via CHARA orbits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The binaries are modeled as two unresolved point sources with a single flux ratio applied to both the H- and K-band data; if the true flux ratios differ between bands, or if either star is slightly resolved at the longest CHARA baselines, the recovered relative astrometry, inclination, and distance would be biased outside the quoted uncertainties.","fun_headline_variants_meta":{"raw":{"variants":["CHARA orbits yield 0.5% masses for two binaries","Visual orbits and RVs pin masses to 0.5%","Long-period binaries get precise masses and ages","Two binaries weighed to 0.5% via CHARA orbits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000328,"raw_usage":{"total_tokens":1938,"prompt_tokens":1156,"completion_tokens":782,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":772,"completion_tokens_details":{"reasoning_tokens":712}},"tokens_in":772,"tokens_out":782,"duration_ms":20658,"temperature":1.0,"reasoning_tokens":712,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:54:51.149465+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit the H- and K-band flux ratios separately to the MIRC-X and MYSTIC data instead of forcing one value: if the resulting flux ratios disagree beyond the fitted uncertainties, or if direct angular-diameter measurements at CHARA's longest baselines show either component to be non-zero in size, the point-source equal-flux-ratio assumption fails and the quoted inclination, distance, and masses would shift.","supporting_citations":[{"cited_title":"C., Tomkin, J., Williamson, M","cited_arxiv_id":null,"evidence_quote":"Prior spectroscopic orbit of HD 210763, used for RV templates, comparison, and literature radial velocities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior spectroscopic orbit and metallicity estimate for HD 221950, used for templates and comparison."}],"review_version":2}