Component masses in stellar and substellar binaries from Gaia astrometry and photometry
Pith reviewed 2026-05-15 16:47 UTC · model grok-4.3
The pith
Gaia astrometry and photometry can determine individual component masses in unresolved binaries without follow-up observations.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using a likelihood approach, the individual component masses are obtained by sampling the posterior probability distribution over the stellar parameters, marginalizing over system age and metallicity, with the flux ratio inferred from mass-flux relations fitted from stellar isochrone models applied to the Gaia three-band photometry and the astrometric orbit.
What carries the argument
Mass-flux relation fitted from stellar isochrone models for each Gaia band, used to infer the unknown flux ratio from the total photometry.
If this is right
- Primary masses are determined with 10-20% precision (one-sigma posterior width) in 90% of cases for main sequence primaries within 300 pc.
- Secondary masses extending to planetary-mass objects are obtained, with half having better than 25% precision.
- Adding infrared photometry changes the mass estimates by less than 4%.
- Spectroscopic orbits from Gaia change the estimates by less than 1%.
- Interstellar extinction has little impact on the mass estimates for this nearby sample.
Where Pith is reading between the lines
- This method could be applied to future Gaia data releases to build a large catalog of binary masses for population studies.
- The approach might extend to other photometric surveys combined with astrometric orbits from other missions.
- Validation against known mass systems could refine the isochrone-based relations for low-mass objects.
Load-bearing premise
The mass-flux relations from stellar isochrone models accurately represent the true brightness-mass relationship for the observed binaries, including low-mass secondaries.
What would settle it
Obtaining independent mass measurements from radial-velocity monitoring or eclipsing binaries for a subset of these systems and comparing them to the Gaia-derived values would test the accuracy and precision.
read the original abstract
The masses of stars and planets can be measured dynamically in binary systems. For an unresolved binary, time series astrometry yields some orbital parameters, but it cannot provide the component masses, because we observe only the motion of the system's photocentre. However, as a star's luminosity is related to its mass, the observable photometry of both components together provides information on the system mass. Here we develop a method to determine the individual component masses of an unresolved binary using the astrometric orbit together with three-band photometry from Gaia. We use a mass-flux relation fitted from stellar isochrone models for each Gaia band to infer the unknown flux ratio. This enables our method to distinguish between near equal-mass, near equal-brightness stellar binaries and star-planet binaries, which otherwise have identical astrometric signatures. Using a likelihood approach, we sample the posterior probability distribution over the stellar parameters, marginalizing over system age and metallicity to get the individual masses. We apply this to 20 000 systems with a main sequence primary within 300 pc of the Sun using data from the Gaia data release 3 non-single star catalogue. Primary masses can be determined with a precision (one-sigma posterior width) of 10-20% in 90% of cases. Secondary masses, which extend down to planetary-mass objects, are less precise, although half are more than 25% precise. Interestingly, adding either infrared photometry or spectroscopic orbits from Gaia does not change the mass estimates much (less than 4% and 1% respectively). Interstellar extinction likewise has little impact for this sample. This work shows that reasonably precise masses can be obtained for stars and substellar objects using just the Gaia astrometry and photometry without need for extensive follow-up.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a method to determine individual component masses for unresolved binaries by combining Gaia astrometric orbital solutions with three-band photometry. Mass-flux relations are fitted from external stellar isochrone models to infer the unknown flux ratio; a likelihood framework then samples posteriors for the two masses while marginalizing over age and metallicity. The method is applied to ~20 000 systems with main-sequence primaries within 300 pc drawn from the Gaia DR3 non-single-star catalogue, yielding claimed 10–20 % precision on primary masses in 90 % of cases and >25 % precision on half the secondary masses, which extend into the planetary regime.
Significance. If the mass-flux relations prove accurate across the full mass range, the approach would enable dynamical-mass estimates for tens of thousands of binaries using only existing Gaia data, providing a scalable route to statistical studies of the stellar–substellar boundary and binary demographics without extensive follow-up.
major comments (3)
- [§3] §3 (mass-flux relation construction): the relations are taken directly from isochrone grids rather than fitted or validated against the Gaia sample itself; no quantitative residuals or systematic-error budget is reported for masses ≲0.1 M⊙, where atmospheric and deuterium-burning uncertainties are known to be largest. This assumption is load-bearing for the substellar-mass claims and the quoted 25 % precision.
- [§4.2] §4.2 (sample definition and main-sequence assumption): the selection requires a main-sequence primary, yet no explicit test (e.g., comparison with independent luminosity or temperature indicators) is shown to quantify contamination by pre-main-sequence or evolved stars; such contaminants would systematically bias the inferred secondary masses.
- [§5] §5 (precision and validation): posterior widths are presented as the precision metric, but the manuscript contains no external validation against systems with independent dynamical masses in the substellar regime; without such a test the claimed accuracy for planetary-mass secondaries remains unverified.
minor comments (2)
- [Abstract] The abstract states that infrared photometry or spectroscopic orbits change the masses by <4 % and <1 %, respectively; the corresponding quantitative comparison should be shown explicitly in a table or figure.
- [§2] Notation for the flux ratio and the three Gaia bands (G, G_BP, G_RP) should be defined once at first use and used consistently thereafter.
Simulated Author's Rebuttal
We thank the referee for their thorough review and constructive comments on our manuscript. We address each of the major comments below and outline the revisions we plan to make.
read point-by-point responses
-
Referee: [§3] §3 (mass-flux relation construction): the relations are taken directly from isochrone grids rather than fitted or validated against the Gaia sample itself; no quantitative residuals or systematic-error budget is reported for masses ≲0.1 M⊙, where atmospheric and deuterium-burning uncertainties are known to be largest. This assumption is load-bearing for the substellar-mass claims and the quoted 25 % precision.
Authors: We agree that providing a quantitative assessment of the mass-flux relations is important. In the revised manuscript, we will add an analysis of the residuals between the isochrone-derived mass-flux relations and literature mass-luminosity relations for low-mass stars and brown dwarfs. We will also include a discussion of systematic uncertainties in the substellar regime due to atmospheric models and deuterium burning. This will be presented as an appendix or extended section in §3. revision: yes
-
Referee: [§4.2] §4.2 (sample definition and main-sequence assumption): the selection requires a main-sequence primary, yet no explicit test (e.g., comparison with independent luminosity or temperature indicators) is shown to quantify contamination by pre-main-sequence or evolved stars; such contaminants would systematically bias the inferred secondary masses.
Authors: We will strengthen the sample selection section by including a quantitative estimate of potential contamination. Specifically, we will compare the positions of our sample in the Gaia color-magnitude diagram with isochrones for different ages and discuss the expected fraction of pre-main-sequence or evolved stars within 300 pc. If possible, we will cross-match with spectroscopic surveys for temperature validation. revision: yes
-
Referee: [§5] §5 (precision and validation): posterior widths are presented as the precision metric, but the manuscript contains no external validation against systems with independent dynamical masses in the substellar regime; without such a test the claimed accuracy for planetary-mass secondaries remains unverified.
Authors: We acknowledge the importance of external validation. However, there are currently very few unresolved binaries with Gaia astrometric solutions that also have independent dynamical mass measurements in the substellar or planetary-mass regime. We will add a section discussing the limited available validation samples, such as any known systems from radial velocity or other methods, and compare our inferred masses where possible. We will clarify that the reported precisions are based on posterior widths and highlight the need for future observational validation. revision: partial
- Lack of sufficient external validation data for planetary-mass secondaries
Circularity Check
No circularity: mass-flux relations drawn from external isochrone models
full rationale
The derivation combines Gaia astrometric orbits with three-band photometry to infer component masses. It adopts mass-flux relations that were fitted once from stellar isochrone models (external to the 20 000 Gaia systems) and then uses a likelihood to sample posteriors while marginalizing age and metallicity. No equation or step defines a quantity in terms of itself, renames a fitted parameter as a prediction, or relies on a load-bearing self-citation whose validity is internal to the paper. The central result therefore remains independent of the target sample and does not reduce to its inputs by construction.
Axiom & Free-Parameter Ledger
free parameters (1)
- mass-flux relation coefficients
axioms (2)
- domain assumption Mass-luminosity relation from isochrones holds for both components
- domain assumption All targets are main-sequence stars with the modeled properties
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.