{"id":"801e8f01-99f1-4718-8f5a-35f58637f623","arxiv_id":"1908.02719","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"HD 81357 is a 33.77-day semi-detached binary at about 63 degrees inclination, with a 3.36 solar-mass B star and a 0.34 solar-mass Roche-lobe-filling cool companion in the late stages of mass exchange.","lead":"Astronomers measured the little-known binary star HD 81357 in detail for the first time, finding a hot B star orbited by a cool, bloated, low-mass star that is spilling gas onto it. The system appears to be a late-stage mass-exchange binary, a useful test case for theories of how binary stars trade mass.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Table 6 is not self-consistent: the quoted a=63.01 R_sun contradicts Kepler's third law for P=33.774 d and M1+M2=3.70 M_sun, which require a≈68 R_sun; the reported physical elements rest on an inconsistent PHOEBE solution.","rationale":"The paper is a careful multi-observatory study, and the binary nature is well supported by anti-phased radial velocities and a passband-dependent double-wave light curve. The central quantitative claim, however, rests on the PHOEBE solution reported in Table 6. A reader cannot verify that solution because the quoted semi-major axis violates Kepler's third law when combined with the quoted period and masses. This is not a disagreement with an outside model; it is an internal inconsistency in the headline parameters. The same table contains the abstract's overstated R2 uncertainty (13.97±0.05 versus 14.0±0.7), and the text concedes that chi^2 varies by less than 3% over broad parameter ranges, so the formal errors are not informative. The reader's weakest assumption about degenerate light-curve modeling is related, but the Kepler violation is the sharper and more easily settled problem: if a truly is 63.01 R_sun, then either the masses or the period are wrong; if a should be about 68 R_sun, the Roche-lobe-filling story actually becomes more consistent, and the paper needs a correction but not a change of the central interpretation. I therefore do not propose changing the CONDITIONAL verdict, but the condition should explicitly include verifying and correcting the self-consistency of Table 6 before the quoted physical elements are used.","tokens_in":24047,"tokens_out":24277,"duration_ms":239521,"concrete_test":"Recompute the semi-major axis from the quoted period and masses with Kepler's third law: a = 4.2067 (M1+M2)^(1/3) P^(2/3) R_sun = 68.0 R_sun, and compare with Table 6's a=63.01±0.09. Additionally verify the relation a sin i = P(K1+K2)/(2π) using the adopted K1, K2, and i. Then rerun the PHOEBE solution with a as a free parameter or fixed to the Keplerian value, and check whether the resulting masses, radii, and detached-versus-semi-detached comparison change. This single check settles whether the reported physical elements are internally consistent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing basis for the claimed physical elements is the PHOEBE solution in Table 6 (Section 3.6). As published, that solution is internally inconsistent. With P=33.77445 d, M1=3.36 M_sun, and M2=0.34 M_sun, Kepler's third law gives a=4.2067 (M1+M2)^(1/3) P^(2/3) R_sun = 68.0 R_sun, yet Table 6 quotes a=63.01±0.09 R_sun. The quoted a corresponds instead to M_total=2.95 M_sun, a ~20% deficit. It also disagrees with the adopted RVs: K1+K2=89.5 km/s and i=63 deg give a sin i = P(K1+K2)/(2π) = 59.8 R_sun, i.e. a=67.1 R_sun. The error is not marginal: the reported a is roughly 5-7% too small. This affects the Roche-lobe geometry: with a=63.0, the secondary's lobe radius is about 13.0 R_sun, so R2=14.0 overfills by 7%; with the correct a=68.0, R_L2≈14.1 R_sun and R2 is consistent with filling its lobe. The paper's own statement at the end of Section 3.6 that the total chi^2 changes by less than 3% over a large range of inclinations and mass ratios further indicates that the formal errors in Table 6, and the abstract's R2=13.97±0.05, do not represent the true uncertainty. If the reported a is a misprojected or stale quantity, the physical elements must be recomputed; if it is not, the quoted masses, radii, and Roche-lobe filling conclusion are not supported by a self-consistent model.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a multi-technique study of the little-known Be-type binary HD 81357, based on new and archival spectroscopy and photometry. The authors infer a double-lined spectroscopic binary and ellipsoidal variable with an orbital period of 33.77445 d, an inclination near 63°, component masses of about 3.36 and 0.34 M⊙, radii of about 3.9 and 14.0 R⊙, and a mass ratio near 10. They interpret the system as a semi-detached binary in the final slow phase of mass exchange after mass-ratio reversal, and they support this with MESA evolutionary models that reproduce the observed masses, period, and secondary radius for a set of adopted initial conditions.","tokens_in":24447,"tokens_out":5065,"duration_ms":52514,"significance":"If the derived parameters are correct, HD 81357 is a valuable new representative of the rare class of hot emission-line binaries with a Roche-lobe-filling cool secondary, and the paper's careful multi-pronged approach—independent RV methods, MCMC-based effective-temperature estimation, two independent light-curve codes, a Gaia parallax constraint, and an explicit consistency check between photometric and spectroscopic radii—sets a high methodological standard. The paper is also refreshingly honest about the degeneracies it encounters, particularly the template-dependent hot-component RVs and the flat χ² landscape. However, the central quantitative claims rest on a PHOEBE solution that is internally inconsistent with Kepler's third law, and the quoted error bars do not reflect the admitted degeneracy; both issues must be resolved before the physical elements can be accepted.","major_comments":[{"comment":"The published solution is not self-consistent. With P=33.77445 d, M1=3.36 M⊙, M2=0.34 M⊙, Kepler's third law requires a=4.2067 (M_total)^(1/3) P^(2/3) R⊙ ≈ 68.0 R⊙, whereas Table 6 quotes a=63.01±0.09 R⊙. The RV amplitudes give the same answer: K1+K2=89.5 km/s and i=63° imply a sin i = P(K1+K2)/(2π) ≈ 59.8 R⊙, hence a≈67.1 R⊙. The quoted a corresponds instead to a total mass of about 2.95 M⊙, a ~20% deficit. This is not marginal: it changes the Roche-lobe geometry materially, with R_L2≈13.0 R⊙ for a=63.0 R⊙ (implying an overfill of ~7%) versus R_L2≈14.1 R⊙ for a=68.0 R⊙ (implying nominal lobe filling). The physical elements, the semi-detached conclusion, and the evolutionary comparison in Section 5 therefore rest on a solution that violates the adopted orbital constraints. The authors must recompute the PHOEBE solution with the Kepler relation enforced (or demonstrate that the quoted a and masses are not the ones actually used), and re-derive all derived quantities that depend on a, including the Roche-lobe radius and the MESA comparison.","section":"Section 3.6, Table 6"},{"comment":"The formal errors quoted for the masses and radii are not a realistic measure of the solution's uncertainty. The text states that over a large range of inclinations and a tolerable range of mass ratios the total χ² changes by less than 3%, and Section 3.3.2 concedes that the mass ratio cannot be uniquely determined from the RVs (K1 varies from 4.78 to 11.09 km/s among the three KOREL templates). The Table 6 errors are derived from the covariance matrix only and take none of this degeneracy into account. Yet the abstract and tables report, for example, M1=3.36±0.15 M⊙ and M2=0.34±0.04 M⊙, which will be read as tight constraints. The paper should either quantify the true parameter correlations with a grid-based or MCMC marginalization, or explicitly present the parameters as a representative solution within a broad, unquantified degeneracy. As published, the error budget is understated.","section":"Section 3.6, Abstract and Table 6"}],"minor_comments":[{"comment":"The abstract gives the secondary radius as 13.97±0.05 R⊙, but Table 6 lists R2=14.0±0.7 R⊙. Please reconcile the value and, in particular, the uncertainty, which differs by more than an order of magnitude.","section":"Abstract vs. Table 6"},{"comment":"The caption refers to an observed semimajor axis a=63.95 R⊙, while Table 6 quotes a=63.01±0.09 R⊙. These two numbers should be made consistent, and the source of the difference explained.","section":"Figure 9 caption"},{"comment":"Equation (1) is difficult to parse as typeset; the relationship between the synchronicity parameter F1, the orbital period, v1 sin i, and the equatorial radius should be written with explicit parentheses and units so that the reader can reproduce the iterative update.","section":"Equation (1)"},{"comment":"There are several typos: \"re-caculated\" in Section 3.6, \"absortion\" in Section 4, and \"shorward\" in Section 4. A careful proofreading pass is recommended.","section":"General text"}],"recommendation":"major_revision","confidential_remarks":"The internal Kepler inconsistency in Table 6 is the kind of issue that must be fixed before publication; the central interpretation may survive a recomputation, but the current numerical values cannot be cited. The manuscript also appears to be part of a longer series, and the authors have a track record of careful work, which makes me inclined to request a revision rather than reject, but the load-bearing numbers need to be regenerated and re-validated against the Gaia parallax and the RV constraints."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: this is the first detailed study of HD 81357, and it earns a real referee. The data are broad — red, blue, IR spectroscopy, Hvar, Hipparcos, ASAS-SN photometry — and the analysis is genuinely multi-pronged: KOREL disentangling, asTODCOR RVs, PYTERPOL temperatures, PHOEBE and BINSYN light-curve solutions, a Gaia parallax consistency check. The paper is honest about the degeneracies in an ellipsoidal variable; it does not hide that the mass ratio is not uniquely pinned by RVs alone.\n\nThat said, there is a problem in Table 6 that needs to be fixed before the physical elements are usable. With P=33.77445 d and M1+M2=3.70 M_sun, Kepler's third law gives a≈68 R_sun, not the quoted 63.01±0.09. The K1+K2 sum from the same table gives a similar value (about 67 R_sun at i=63°). So the published solution is internally inconsistent. This is not a cosmetic typo: the Roche-lobe radius of the secondary changes by about 1 R_sun between a=63 and a=68, which is exactly the threshold for whether R2=14.0 fills its lobe. The paper's semi-detached interpretation hinges on that.\n\nThe abstract's R2=13.97±0.05 is also inconsistent with Table 6's 14.0±0.7, and the paper's own statement that the total chi^2 changes by less than 3% over a large range of inclinations and mass ratios shows that the formal errors understate the true uncertainty. These are addressable in revision, but they are not just presentational.\n\nThe MESA section is a tuned inverse search and is best read as an illustration that the observed state can be reached by mass transfer, not as a determination of the initial conditions. The paper itself is appropriately cautious about that.\n\nBottom line: if the Table 6 inconsistency is resolved — meaning the PHOEBE solution is re-derived with a self-consistent a, or the quoted a is corrected and the elements recomputed — this becomes a valuable new anchor point for intermediate-mass interacting binaries near the end of mass exchange. I would send it to review, with a request to fix the internal consistency and to replace the formal errors with something that reflects the real degeneracy. It deserves a serious referee; the raw material is good and the caveats are mostly stated.","headline":"A careful, data-rich first study of a newly identified Be binary, but the quoted semimajor axis contradicts Kepler's third law given the paper's own masses and period, and the formal errors are unrealistic.","tokens_in":25131,"tokens_out":3007,"would_cite":true,"duration_ms":30400,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"HD 81357 is a 33.77-day double-lined binary with a 10:1 mass ratio and a Roche-lobe-filling cool star, observed in the final slow phase of mass exchange.","keywords":["Stars: close","Stars: binaries: spectroscopic","Stars: emission-line, Be","Stars: fundamental parameters","Stars: individual: HD 81357","ellipsoidal variable","Roche-lobe overflow","mass transfer"],"falsifier":"A clean radial-velocity curve of the hot component, obtained from lines free of circumstellar emission (for example in the ultraviolet), would directly test the mass ratio: if its semi-amplitude K1 differs significantly from the adopted value near 8 km/s, the masses 3.36 and 0.34 solar masses would have to be revised. Alternatively, high-angular-resolution observations that resolved the cool star's disk could test the central geometrical premise that its radius equals its Roche-lobe radius.","tokens_in":23806,"feed_emoji":"🔭","tokens_out":8424,"duration_ms":78966,"temperature":0.7,"pith_summary":"This paper tries to establish what the little-studied Be star HD 81357 actually is: not a single star but a 33.77-day double-lined spectroscopic binary and ellipsoidal variable. From new spectra, archival and new photometry, and a very accurate trigonometric parallax, the authors derive component masses of 3.36 and 0.34 solar masses, radii of 3.9 and about 14 solar radii, an inclination near 63 degrees, and a mass ratio of about 10. They argue that the cool, low-mass star fills its Roche lobe and that the system is seen in the final slow phase of mass exchange after the mass-ratio reversal, with the hot star as the mass gainer. A sympathetic reader would care because such systems are rare benchmarks for testing the theory of mass transfer in close binaries, and because the paper shows how a precise parallax can break the usual degeneracies of ellipsoidal variables.","feed_headline":"HD 81357 is a 10-to-1 binary in its final mass-exchange phase","feed_subtitle":"Spectra and light curves yield a 3.36 and 0.34 solar-mass pair with a cool star filling its Roche lobe.","key_machinery":"The load-bearing machinery is the simultaneous model of the ellipsoidal light curve and the radial-velocity curve of the cool component, anchored by a very accurate trigonometric parallax. The ellipsoidal variations carry the geometric information: their amplitude and its decrease from red to blue passbands require the cool star to fill its Roche lobe, and the parallax converts the measured flux of the hot component into a physical radius, which in turn fixes the scale of the orbit. Effective temperatures of both stars are first pinned down by an automated comparison of observed and synthetic blue spectra; the light-curve and orbital model then yields masses, radii, inclination, and the mass ratio. A binary stellar evolution code provides the time dimension, showing that conservative mass exchange from an initially more massive star produces the present-day masses, period, and the inflated cool star.","core_discovery":"On the paper's own terms, the central discovery is that HD 81357 is a new member of the class of hot emission-line binaries with a Roche-lobe-filling cool secondary, and that it is caught near the end of its mass-transfer phase. The double-wave ellipsoidal light variations, the clean sinusoidal radial-velocity curve of the cool star with semi-amplitude 81.8 km/s, and the broad, partly circumstellar-veiled lines of the hot star together define a binary with period 33.77445 days and a circular orbit. A combined fit to all photometry and radial velocities, constrained by an accurate parallax, gives masses 3.36 and 0.34 solar masses, radii 3.9 and 14.0 solar radii, and mass ratio 10.0. Because detached models cannot reproduce the light-curve amplitudes, the authors conclude that the cool star fills its Roche lobe. Evolutionary modelling then shows that a binary starting with roughly 1.5 and 2.2 solar masses in a 2.4-day orbit evolves conservatively into exactly such a system, so HD 81357 is interpreted as being observed in the final slow phase of mass exchange after the mass-ratio reversal.","pith_inferences":["The same recipe of using a precise parallax to convert ellipsoidal light-curve amplitudes into physical radii could break the mass-ratio degeneracy in other emission-line binaries with cool Roche-lobe-filling stars, where spectral disentangling alone leaves a wide range of acceptable mass ratios.","If the system is truly near the end of mass transfer, its orbital period should be changing only very slowly; monitoring the times of the ellipsoidal light-curve minima over decades could test this prediction.","The model predicts specific passband-dependent ellipsoidal amplitudes; space-based ultraviolet photometry, where the hot star dominates, would provide an independent check of the inclination and of the cool star's contribution.","The presence of long-term cyclic variations in some similar systems but not in this one may indicate that those cycles require an actively sustained mass-transfer rate or an extended outer disk, not merely the existence of a Roche-lobe-filling donor."],"forward_implications":["HD 81357 becomes a well-measured example of a Be binary in the late stage of Algol-like mass exchange, useful for calibrating the end points of mass-transfer theory.","The derived masses, radii, and temperatures give direct constraints on the structure of a low-mass star with an inflated envelope after most of its mass has been transferred.","Because the cool star fills its Roche lobe, the ellipsoidal light curve is a monitor of the system's present geometry; future photometry can check for slow changes in the amplitudes.","The absence of cyclic light variations on timescales an order of magnitude longer than the orbital period, in contrast with similar binaries, supports the suggestion that such cycles are not present at the end of the mass-exchange phase.","The H-alpha emission follows the hot component's orbital motion with a phase shift of about five days, indicating an asymmetric distribution of circumstellar matter around the mass-gaining star."],"supporting_citations":[{"why":"Reported the 33.8-day periodicity in metallic-line radial velocities and Hipparcos photometry, establishing the binary nature that this paper verifies and measures.","marker":"Koubský et al. (2012)"},{"why":"The predecessor study of BR CMi that supplies the iterative analysis strategy and the finding that a wide range of mass ratios fits the disentangled spectra.","marker":"Harmanec et al. (2015)"},{"why":"Demonstrates the 2-D cross-correlation approach for hot binaries with circumstellar matter and provides the AU Mon comparison used in the discussion.","marker":"Desmet et al. (2010)"},{"why":"Describes the light-curve and radial-velocity modelling code used for the combined photometric and orbital solution.","marker":"Prša & Zwitter (2005, 2006)"},{"why":"Provides the accurate parallax that the paper uses to fix the physical radius scale and restrict the mass-ratio range.","marker":"Gaia Collaboration et al. (2018)"},{"why":"The stellar evolution code that models conservative mass exchange and reproduces the observed masses, period, and secondary radius.","marker":"Paxton et al. (2015)"},{"why":"Describes the synthetic-spectrum interpolation tool used to determine the component effective temperatures and projected rotational velocities.","marker":"Nemravová et al. (2016)"},{"why":"Establishes the 2-D cross-correlation technique used to measure the hot component's radial velocities from the blue spectra.","marker":"Zucker & Mazeh (1994)"}],"fun_headline_variants":["HD 81357: a 10-to-1 binary caught in late mass exchange","Binary with 10:1 mass ratio enters final mass-exchange phase","Roche-filling cool star reveals 10:1 binary's late mass transfer","HD 81357: mass-ratio 10 binary after reversal, final slow phase","Ellipsoidal binary with 10:1 mass ratio near mass-exchange end"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result rests on the assumption that the ellipsoidal variability comes entirely from a cool star that exactly fills its Roche lobe, rotates synchronously, and obeys the adopted gravity-darkening and albedo laws; because the model's total chi-squared changes by less than three percent over a wide range of inclinations and mass ratios, a different acceptable set of these assumptions could shift the inferred masses and inclination.","fun_headline_variants_meta":{"raw":{"variants":["HD 81357: a 10-to-1 binary caught in late mass exchange","Binary with 10:1 mass ratio enters final mass-exchange phase","Roche-filling cool star reveals 10:1 binary's late mass transfer","HD 81357: mass-ratio 10 binary after reversal, final slow phase","Ellipsoidal binary with 10:1 mass ratio near mass-exchange end"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000813,"raw_usage":{"total_tokens":3668,"prompt_tokens":1152,"completion_tokens":2516,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":768,"completion_tokens_details":{"reasoning_tokens":2422}},"tokens_in":768,"tokens_out":2516,"duration_ms":18662,"temperature":1.0,"reasoning_tokens":2422,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:37:42.519158+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A clean radial-velocity curve of the hot component, obtained from lines free of circumstellar emission (for example in the ultraviolet), would directly test the mass ratio: if its semi-amplitude K1 differs significantly from the adopted value near 8 km/s, the masses 3.36 and 0.34 solar masses would have to be revised. Alternatively, high-angular-resolution observations that resolved the cool star's disk could test the central geometrical premise that its radius equals its Roche-lobe radius.","supporting_citations":[{"cited_title":"2010, , 401, 418","cited_arxiv_id":null,"evidence_quote":"Demonstrates the 2-D cross-correlation approach for hot binaries with circumstellar matter and provides the AU Mon comparison used in the discussion."}],"review_version":1}