Pith. Sign in

REVIEW

The EBLM project - XX. Absolute dynamical masses of five low-mass eclipsing binaries, and an evaluation of mass estimation methods

Solar and Stellar Astrophysics · Submitted 13 Oct 2026 · Published version · 166 references

Pith's verdict: Five new EBLM absolute masses are solid and worth having; the SB1 bias mechanism is plausible, but the 'almost solely' attribution outruns the evidence — the missing R1 residual check is one figure away.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 14, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.18339 (version 1)
Pith IDHBQ7QHJ6, Pith's permanent record of this paper and its review
CitationAdam T. Stevenson, Alexandre Santerne, Alix V. Freckelton, Amaury H. M. J. Triaud, Daniel Sebastian, David V. Martin, Neda Heidari, Nicola J. Miller, Pierre F. L. Maxted, Thomas A. Baycroft, Yasmin T. Davis. "The EBLM project - XX. Absolute dynamical masses of five low-mass eclipsing binaries, and an evaluation of mass estimation methods." arXiv:2610.18339 (2026). https://pith.science/paper/HBQ7QHJ6
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Extracting the hidden M dwarf spectrum in five eclipsing binaries yields direct stellar masses and shows that single-lined mass estimates run systematically high because the companion's light inflates the isochrone primary radius.

The reading

This paper reports the detection of the spectral signature of the M dwarf secondary star in five eclipsing binaries from the EBLM project, whose secondary stars contribute less than one per cent of the optical light. The detections, made with high-resolution cross-correlation spectroscopy, convert the systems from single-lined to double-lined binaries and yield model-independent dynamical masses for both stars. Comparing those masses with the standard single-lined method—radial velocities plus eclipse photometry plus an isochrone-derived primary radius—the paper finds a systematic overestimation of about 8 per cent, and traces its source to the isochrone fitting: the companion's light shifts the measured photometric magnitudes and inflates the recovered primary radius. An injection-recovery simulation quantifies the inflation as a function of the binary mass ratio, and a mass-ratio-dependent correction removes the discrepancy. If the correction holds up, the far larger single-lined EBLM sample can be calibrated to anchor mass-radius-luminosity relations for the smallest and most common stars in the Galaxy.

What carries the argument

The load-bearing machinery is the pairing of a high-resolution cross-correlation spectroscopy pipeline with a model-dependent comparison. On the dynamical side, singular-value decomposition removes the primary's spectrum chunk by chunk, the K-focusing method builds a cross-correlation map in the secondary's rest frame by stepping over trial semi-amplitudes $K_2$, and the Saltire model—a double-Gaussian cross-correlation shape with side-lobes, fitted through the same focusing process—extracts $K_2$, with systematic errors estimated by time-splitting the data or bootstrapping. On the single-lined side, the mass function $f = M_2^3/(M_1+M_2)^2 = P K_1^3 \sin^3 i\,(1-e^2)^{3/2}/(2\pi G)$ is solved together with the Keplerian relation $M_1/R_1^3 = 4\pi^2(a/R_1)^3/(G P^2)$, which requires the isochrone primary radius $R_1$ as its one external input. The bias argument rides on that input: the injection-recovery grid measures how much the companion's light inflates the recovered $R_1$ as a function of mass ratio, and the fitted correction $f(q) = A\exp(-B q) + C$ with $(A,B,C) = (1.2786, -3.6498, -1.5167)$ is the device that lets the full single-lined sample be recalibrated.

Extended reading notes

Core claim

The central claim is that single-lined mass estimates for EBLM binaries run systematically high because the unresolved M dwarf companion's light changes the spectral energy distribution used in isochrone fitting, so the fitted primary radius—the single external input the single-lined method needs—comes out too large. Over a grid of simulated binaries covering primary masses from 0.65 to 1.5 solar masses and secondaries from 0.1 to 0.65, fitting a single-star model to MIST isochrones recovers primary radii systematically above the input truth, with the overestimate growing with mass ratio $q$ and reaching several per cent for the most massive, most luminous secondaries. The astrometric parallax zero-point offset and terminal-age bias are examined and rejected as the dominant cause. Feeding the corrected radii back through the single-lined method brings the resulting masses into scatter around zero relative to the HRCCS dynamical masses, which the authors read as showing that the discrepancy between the two approaches can be almost solely attributed to single-star isochrone modelling. The newly measured secondary masses are among the most accurate and precise fundamental parameters available in the M dwarf regime.

Load-bearing premise

The quantitative attribution of the single-lined mass bias to isochrone radius inflation rests on the simulated injection-recovery grid—one stellar evolution grid at a fixed age of 1 Gyr, distance of 300 pc, and extinction of $A_V = 0.2$—being representative of real EBLM systems, and on the visually chosen exponential form of the correction being the right function; the authors themselves state the correction is preliminary.

Editorial extensions

If this is right

  • The full EBLM single-lined sample—including the faintest, lowest-mass secondaries that visible-light HRCCS cannot reach—can be recalibrated with the mass-ratio-dependent radius correction to deliver accurate M dwarf masses, extending the sample far beyond the handful of systems accessible to direct detection.
  • The newly double-lined systems roughly double the number of precisely characterised, well-detached (period above 10 days) low-mass binaries in the reference catalogue of detached eclipsing binaries, a region of parameter space that is otherwise nearly empty.
  • For radius-inflation studies, the single-lined method remains usable: because the mass and radius posteriors are strongly correlated, a biased mass moves the inferred secondary along the theoretical isochrone slope rather than away from it.
  • Fundamental effective temperatures derived from the spectral energy distributions of these systems agree with spectroscopic temperatures within $1\sigma$ (RMS scatter 171 K), giving benchmark stars for calibrating stellar atmosphere models and empirical relations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the bias mechanism is driven by the secondary's fractional light contribution, the same isochrone-radius inflation should affect any single-star fit to an unresolved binary—including exoplanet host stars with stellar companions—so the correction has a natural testbed in transiting-planet systems with known faint companions.
  • The authors' robustness checks (varying metallicity, age, and distance leaves the fit unchanged) suggest the effect is geometric in the spectral energy distribution rather than specific to one evolution grid; if so, the same fractional radius overestimate should appear in single-star classification of binary blends in large photometric surveys.
  • The wavelength dependence of the mechanism is testable in the near-infrared, where M dwarfs are brighter relative to their primaries: the radius bias should grow with the secondary's fractional light, even though the improved contrast is exactly what makes near-infrared HRCCS detections of the lowest-mass secondaries feasible.
Original abstract

Eclipsing binaries consisting of a solar-type primary and an M dwarf secondary are key laboratories to study objects at the bottom of the main sequence. Fundamental parameters of such stars have a direct impact upon exoplanetary studies, as searches for temperate Earth-sized planets focus on M dwarf hosts. Eclipsing binaries are essential to measure these parameters. However for stellar masses $<$ 0.35 M$_\odot$ , few double-lined eclipsing systems are known. Thus, stellar evolution models and empirical relations have had to rely on single-lined eclipsing binaries, for instance those from the EBLM project (Eclipsing Binary - Low Mass). Using high resolution spectroscopy to extract the weak signal of the secondary (flux ratio $<$ 1 per cent in the optical), we transform five more systems from single to double-lined, deriving absolute masses for both stars. We then compare methods typically employed to estimate stellar masses for field stars or single-lined binaries with the absolute dynamical masses we obtain, to assess their accuracy. Our results for SB1 and SB2 analyses are mostly consistent within 1-$\sigma$, though we identify a small systematic overestimation of SB1 masses derived from a combination of radial velocities and eclipse photometry. We show that the bias originates from primary radius overestimation during isochronal analysis, wherein the secondary light reduces photometric magnitudes and impacts the recovered radii by modifying the spectral energy distribution. With proper calibration, we can ensure the accuracy of using the full EBLM SB1 sample to constrain mass-radius-luminosity relationships for the smallest, yet most abundant, stars in the Galaxy

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-14: 1 request, 7,725 input tokens.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.79 · human participants 0.02 · controlled comparison 0.25 · randomized assignment 0.01 · states quantitative result 0.79 · compares with prior results 0.95 · mentions code or data availability 0.03.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

Eclipsing binaries consisting of a solar-type primary and an M dwarf secondary are key laboratories to study objects at the bottom of the main sequence. Fundamental parameters of such stars have a direct impact upon exoplanetary studies, as searches for temperate Earth-sized planets focus on M dwarf hosts. Eclipsing binaries are essential to measure these parameters. However for stellar masses $<$ 0.35 M$_\odot$ , few double-lined eclipsing systems are known. Thus, stellar evolution models and empirical relations have had to rely on single-lined eclipsing binaries, for instance those from the EBLM project (Eclipsing Binary - Low Mass). Using high resolution spectroscopy to extract the weak signal of the secondary (flux ratio $<$ 1 per cent in the optical), we transform five more systems from single to double-lined, deriving absolute masses for both stars. We then compare methods typically employed to estimate stellar masses for field stars or single-lined binaries with the absolute dynamical masses we obtain, to assess their accuracy. Our results for SB1 and SB2 analyses are mostly consistent within 1-$\sigma$, though we identify a small systematic overestimation of SB1 masses derived from a combination of radial velocities and eclipse photometry. We show that the bias originates from primary radius overestimation during isochronal analysis, wherein the secondary light reduces photometric magnitudes and impacts the recovered radii by modifying the spectral energy distribution. With proper calibration, we can ensure the accuracy of using the full EBLM SB1 sample to constrain mass-radius-luminosity relationships for the smallest, yet most abundant, stars in the Galaxy

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Eclipsing binaries consisting of a solar-type primary and an M dwarf secondary are key laboratories to study objects at the bottom of the main sequence. background 1.00 0.06 0.02 asserted 1.00
s1 Fundamental parameters of such stars have a direct impact upon exoplanetary studies, as searches for temperate Earth-sized planets focus on M dwarf hosts. background 0.99 0.05 0.02 asserted 1.00
s2 Eclipsing binaries are essential to measure these parameters. background 1.00 0.05 0.03 asserted 1.00
s3 However for stellar masses $<$ 0.35 M$_\odot$ , few double-lined eclipsing systems are known. gap 0.99 0.10 0.02 asserted 1.00
s4 Thus, stellar evolution models and empirical relations have had to rely on single-lined eclipsing binaries, for instance those from the EBLM project (Eclipsing Binary - Low Mass). background 0.70 0.07 0.03 asserted 1.00
s5 Using high resolution spectroscopy to extract the weak signal of the secondary (flux ratio $<$ 1 per cent in the optical), we transform five more systems from single to double-lined, deriving absolute masses for both stars. method 0.55 0.97 0.08 asserted 1.00
s6 We then compare methods typically employed to estimate stellar masses for field stars or single-lined binaries with the absolute dynamical masses we obtain, to assess their accuracy. contribution 0.54 0.92 0.07 asserted 0.81
s7 Our results for SB1 and SB2 analyses are mostly consistent within 1-$\sigma$, though we identify a small systematic overestimation of SB1 masses derived from a combination of radial velocities and eclipse photometry. result 1.00 0.96 0.33 asserted 0.87
s8 We show that the bias originates from primary radius overestimation during isochronal analysis, wherein the secondary light reduces photometric magnitudes and impacts the recovered radii by modifying the spectral energy distribution. interpretation 0.83 0.89 0.15 asserted 1.00
s9 With proper calibration, we can ensure the accuracy of using the full EBLM SB1 sample to constrain mass-radius-luminosity relationships for the smallest, yet most abundant, stars in the Galaxy interpretation 0.85 0.56 0.59 hedged possibility 0.73

What kind of paper the classifier reports this is: empirical study 1.00 · evidence from observation or experiment 0.97.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Singular value decomposition and reconstruction for a given spec￾tral chunk. The principal components (PCs) are identified from all available spectra for a particular target. An example spectral chunk is shown in the top panel of this plot. The first few PCs of this spectrum are reconstructed from the identified common features of the primary star. Removing more than 2 or 3 PCs, the reconstructed ‘cleaned’ spectrum no longer shows the absorption features present in the input spectrum. Each spectrum/component is shifted arbitrarily for visibility, and components beyond the first are…
Figure 2. HRCCS detections of targets listed in
Figure 3. An example light curve fit, shown for EBLM J0608-59. The TESS SPOC data are plotted as purple points, with empty circles representing binned data. Top panels: phase-folded primary eclipses and residuals for the best-fitting model. Bottom panels: phase-folded secondary eclipses and residials for the best-fitting model. In both plots, 50 random posterior samples are displayed (cyan curves). 𝑇eff and [Fe/H] are used with the Gaia DR3 parallax, 2MASS J, H, and K magnitudes (Skrutskie et al. 2006), and AllWISE W1, W2, and W3 magnitudes (Cutri et al. 2021) to estimate the primary radii. No
Figure 4. Figures comparing primary masses from different methods. Top: primary masses of the eight systems included in our sample (
Figure 5. These are presented as absolute distance offsets (upper panel) and percentage offsets (lower panel) from the corrected par￾allax distance. For most systems, both distances are in agreement to within 1𝜎, with overlapping error bars. J0608-59 (BEBOP-1) is only consistent to slightly over 1𝜎. Nevertheless, the parallax zero point offset in Gaia DR3 does represent a systematic bias towards larger distances for the targets in our small sample. The magnitude varies from star to star, though the sign of the correction required is always constant. The re-interpolated isochrones, using…
Figure 6. Radii calculated from isochrones using corrected Gaia DR3 par￾allaxes are plotted along 𝑦 = 0. Originally determined radii (from e.g. Freck￾elton et al. 2024 and Section 3) are plotted as a percentage overestimation of the corrected value. The largest impact is observed for J0608-59, where the radius was ∼ 1 per cent too large. et al. (2018) investigated this by simulating combined binary spectra, mimicking those obtained through large scale surveys like APOGEE, GALAH, and LAMOST. They found that best-fitting effective tem￾peratures can be systematically under-estimated in systems…
Figure 7. A comparison between SB1 results when using either original or parallax-zero-point-corrected primary radii. The values are plotted as percentage deviation from the Saltire SB2 results. The updated SB1 masses are in a different colour, and indicated with square markers. For most systems, the difference is almost negligible, and we see that this does not account for the identified systematic bias in results from this method. J1928-38 shows a larger modification (albeit fully within 1𝜎 uncertainties) that brings the quoted median mass more into alignment with the SB2 measurement. The…
Figure 8. Figures showing the radius difference between the recovered SingleStarModel isochrone for our simulated binary systems, and the ground truth primary radius from the MIST track used to generate the system parameters. Top panel: as a function of secondary mass. Bottom panel: as a function of mass ratio. The amount of radius over-estimation increases in a consistent manner as a function of 𝑞. The over-estimation is similar no matter the mass of the simulated primary (indicated by line colours), though is slightly stronger for the least massive primaries – a consequence of cooler stars…
Figure 9. Similar to Figures 4 and A2, we now include the SB1 mass estimates when adjusting the input primary isochrone radius based on the over-estimation caused by using a single star model (light blue points). The adjusted values are far closer to the SB2 measurements, with typical scatter around zero difference. Top: Primary Masses. Bottom: Secondary masses. Both panels are ordered from left to right by increasing primary mass. sample of targets. Inspecting the CCFs, we see that the secondary is not visible above the noise level in an individual CCF, in agreement with SNRsec < 1. The…
Figure 10. Left: The joint 𝑀2, 𝑅2 posterior from the SB1 method (adjusted for 𝑅1 over-estimation), with a band over-plotted indicating the secondary mass from an SB2 analysis. Contours are drawn at the 1, 2, and 3𝜎 levels. One can clearly see that the marginalised 1𝜎 medians and uncertainties do not accurately represent such a correlated posterior distribution. Right: A re-calculation of the secondary radius using observables from the RV + photometric modelling and the total system mass from Saltire as the one piece of assumed information. The resulting distribution lies ever-so-slightly…
Figure 11. Mass–radius diagrams showing 1-𝜎 contours for the eight systems compared in this work. The theoretical BHAC15 5 Gyr isochrone is plotted in black (Baraffe et al. 2015). Blue ellipses show the Saltire results, with radii calculated by combining total mass 𝑀 with RV + Photometry posteriors. The original RV + Photometry results are shown as grey contours, and the pink ellipses show the results after correcting the MIST primary radius (Section 4.4.3). stars more massive than 1.3 M⊙ or with signs of a tertiary compan￾ion. We have chosen to use 10 d as an approximate, though perhaps…
Figure 12. Masses of DEBCat systems, where binaries with period below 10 d have been greyed-out/represented with cross markers. Three EBLM binaries compared in this work are already in DEBCat (J0608-59, J0723+79, Kepler￾16), as is EBLM J0113+31, all indicated by the circled points. With Saltire SB2 masses overlaid in larger dark blue diamond markers, one can see that some previously included systems are not exactly centred, and have changed slightly. The additional systems aid in expanding coverage of parameter space for low-mass stars. N.B: all points include error bars, yet most are smaller…
Figure 13. A comparison between the fundamental effective temperatures for 7 out of 8 stars in our sample, and 𝑇eff determined via spectroscopic methods (see Section 3.3, and Freckelton et al. 2024). The dashed line shows 1:1 agreement. Blue circles are those where we have calculated the fundamental 𝑇eff in this work, and black squares are those previously published (Maxted et al. 2022a, 2025) (2025) compared two sources of 𝑇eff, derived from spectral observa￾tions of a sample of bright single stars, and found an RMS difference of 101 K. This is fully consistent with the typical limiting…

Figures reproduced from arXiv:2610.18339 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.