arxiv: 2604.12492 · v1 · submitted 2026-04-14 · 🌌 astro-ph.SR

Recognition: unknown

The μ Herculis system solved after nearly three centuries

Marcus L. Marcussen , Mikkel N. Lund , Frank Grundahl , Daniel Huber , Emil Knudstrup , Adam L. Kraus , Christoph Baranec , Pere L. Pall\'e

show 6 more authors

Trent J. Dupuy Guillaume Huber James Ou Zach Werber Ruihan Zhang Reed Riddle

Authors on Pith no claims yet

Pith reviewed 2026-05-10 15:54 UTC · model grok-4.3

classification 🌌 astro-ph.SR

keywords μ Herculisquadruple star systemdynamical massesasteroseismologyhierarchical orbitsastrometryradial velocities

0 comments

The pith

Precise dynamical masses have been determined for all four stars in the μ Herculis quadruple system from a joint fit to three centuries of observations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

μ Herculis is a bright nearby quadruple star whose brightest member shows solar-like oscillations, so model-independent masses are needed to turn the system into a benchmark for asteroseismology. The authors combine radial velocities with relative and absolute astrometry from Hipparcos, Gaia DR3, and ground-based catalogues spanning nearly three centuries. They fit the inner Aa-Ab and B-C orbits, the wide A-BC orbit, and the centre-of-mass motion simultaneously. This produces sub-percent precision masses for every component and a reconciled system parallax. A sympathetic reader cares because these masses anchor tests of stellar structure and evolution without relying on theoretical models.

Core claim

Through a simultaneous forward model of the inner Aa-Ab and B-C Keplerian orbits, the wide A-BC orbit, and the sky motion and parallax of the total centre of mass, constrained by radial velocities and astrometry spanning nearly three centuries including the 2023 Aa-Ab periastron passage, the component masses are determined to sub-percent precision as M_Aa = 1.134 ± 0.007 M⊙, M_Ab = 0.2286 ± 0.0006 M⊙, M_C = 0.445 ± 0.005 M⊙, and M_B = 0.417 ± 0.005 M⊙, together with a system parallax of ϖ_CM = 120.069 ± 0.089 mas.

What carries the argument

The joint forward-modelling framework that simultaneously constrains the three Keplerian orbits and the centre-of-mass motion using radial velocities plus absolute and relative astrometry.

If this is right

The mass of the oscillating star Aa is now known independently of stellar models to better than one percent.
All orbital elements of the hierarchical quadruple system are fully determined.
The reconciled parallax improves the distance and luminosity of the entire system.
These masses supply direct calibration points for stellar evolution models across 0.2 to 1.1 solar masses.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The masses can anchor improvements to asteroseismic scaling relations applied to other solar-like oscillators.
Continued monitoring could test whether any non-Keplerian effects appear over longer baselines.
The same joint-fitting approach may resolve other long-period multiple systems that have remained unsolved for centuries.

Load-bearing premise

The three orbits are purely Keplerian with no significant perturbations from unseen companions or relativistic effects over the 300-year baseline, and the joint fit correctly separates the center-of-mass motion from the internal orbits.

What would settle it

A future high-precision astrometric or radial-velocity measurement that shows statistically significant residuals from the predicted orbital paths or masses.

Figures

Figures reproduced from arXiv: 2604.12492 by Adam L. Kraus, Christoph Baranec, Daniel Huber, Emil Knudstrup, Frank Grundahl, Guillaume Huber, James Ou, Marcus L. Marcussen, Mikkel N. Lund, Pere L. Pall\'e, Reed Riddle, Ruihan Zhang, Trent J. Dupuy, Zach Werber.

**Figure 2.** Figure 2: Sky-plane trajectory and residuals for the absolute-astrometry solution. [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗

**Figure 3.** Figure 3: Relative astrometry of the µ Her A subsystem. The plot is centred on the primary µ Her Aa. Data points are colour-coded by epoch, progressing from purple to yellow. To illustrate the constraints provided by relative astrometry alone, we display 600 random draws from the corresponding posterior distribution (red traces). The solid black lines trace the orbit derived from our full joint fit, with the median … view at source ↗

**Figure 4.** Figure 4: Relative astrometry of the BC subsystem. The plot is centred on the primary of that system, [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗

**Figure 5.** Figure 5: Relative astrometry of the outer A–BC system. The plot is centred on [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗

**Figure 6.** Figure 6: Radial-velocity data and model fits for µ Her Aa. Top: Radial velocities versus time. The points are coloured by instrument (see legend). Note that SONG RVs are nightly medians. Green curves represent random draws from the posterior of the RV-only solution, while black curves depict the joint solution; the median joint model is overlaid in white. The dashed lines trace the long-term trends: the green line … view at source ↗

**Figure 7.** Figure 7: Comparison of the primary mass, MAa, derived in this work (vertical grey band, showing the 1σ credibility interval) against literature values. The markers represent estimates from non-seismic spectroscopy (Fuhrmann 1998), stellar modelling (Yang & Meng 2010; Li et al. 2019), and asteroseismic analyses (Grundahl et al. 2017; Gupta et al. 2025). With the advent of higher-precision seismic constraints, the in… view at source ↗

read the original abstract

$\mu$ Herculis is a bright, nearby quadruple system. Its brightest member, $\mu$ Her Aa, displays solar-like oscillations, establishing the system as a crucial benchmark for asteroseismology, provided that its mass can be determined independently of stellar models. We aim to resolve the full hierarchical architecture of the system and determine precise, model-independent dynamical masses for all four components (Aa, Ab, B, and C), along with a consistent astrometric solution for the system's centre of mass. We performed a joint fit of radial velocities, relative astrometry and absolute astrometry from \textit{Hipparcos}, \textit{Gaia} DR3, and ground-based catalogues, spanning nearly three centuries. Our forward-modelling framework simultaneously constrains the Keplerian orbits of the inner Aa--Ab and B--C subsystems, the wide A--BC orbit, and the sky motion and parallax of the total centre of mass. Leveraging several complementary datasets and the decisive 2023 periastron passage of the Aa--Ab pair, we precisely determine all orbital parameters and obtain sub-percent precision on the component masses: $M_{\rm Aa} = 1.134 \pm 0.007\,M_{\odot}$, $M_{\rm Ab} = 0.2286 \pm 0.0006\,M_{\odot}$, $M_{\rm C} = 0.445 \pm 0.005\,M_{\odot}$, and $M_{\rm B} = 0.417 \pm 0.005\,M_{\odot}$. We derive a system parallax of $\varpi_{\rm CM} = 120.069 \pm 0.089\,\mathrm{mas}$ that reconciles and improves upon the individual \textit{Hipparcos} and \textit{Gaia} DR3 values.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

μ Herculis now has the first sub-percent dynamical masses for all four components from a joint fit that folds in the 2023 periastron.

read the letter

The main result is a full hierarchical orbital solution for the μ Herculis quadruple with dynamical masses at the sub-percent level: 1.134 ± 0.007 solar masses for the oscillating primary, plus tight values for the other three stars. This comes from one forward model that fits radial velocities, relative astrometry, and absolute positions from Hipparcos, Gaia, and ground data across nearly 300 years. The recent periastron passage supplies the leverage that tightens the inner Aa-Ab orbit and lets the masses follow directly from the fitted elements via standard Keplerian relations. The parallax they derive also sits cleanly between the earlier Hipparcos and Gaia values. That combination is what makes the work useful as a benchmark for asteroseismic scaling relations. The modeling is straightforward and the circularity is low because the masses are computed from the orbital parameters rather than assumed inside the fit. The data volume and the new observation give the central claims real support. The wide A-BC orbit and center-of-mass motion are fitted together, and the quoted uncertainties look consistent with the constraints available. The stress-test worry about degeneracy between the long-period wide orbit and linear CM motion has some basis given the sparse high-precision absolute points, but the joint data set appears to separate them at the level claimed. The purely Keplerian assumption over three centuries is standard and reasonable here, though small perturbations would be hard to detect. No code or data release is mentioned, which is a small practical limitation for a paper built on public catalogs. This paper is for people who need model-independent masses for solar-like oscillators or hierarchical systems. It is a careful application of existing tools to a high-value target rather than a new technique. The thinking is clear and the literature is engaged properly. I would bring it to a reading group on stellar benchmarks and send it to peer review.

Referee Report

0 major / 2 minor

Summary. The manuscript presents a solution to the μ Herculis quadruple system by conducting a joint forward-model fit to radial velocity measurements, relative astrometry, and absolute astrometry from Hipparcos, Gaia DR3, and ground-based catalogues spanning nearly three centuries. This allows determination of all orbital parameters for the Aa-Ab, B-C, and wide A-BC orbits, yielding dynamical masses with sub-percent precision and a consistent center-of-mass parallax that reconciles previous measurements.

Significance. If the results hold, this provides important model-independent masses for a bright star exhibiting solar-like oscillations, establishing it as a benchmark for asteroseismology. The approach is strengthened by incorporating the 2023 periastron observation of the Aa-Ab pair, which is decisive for constraining the inner orbit. The simultaneous fit of multiple datasets and hierarchical levels is a strength, as is the low circularity in deriving masses from orbital elements using Newtonian gravity. The degeneracy concern between the wide orbit and center-of-mass motion does not appear to land given the data constraints and quoted uncertainties.

minor comments (2)

The abstract could briefly note the total number of data points or the degrees of freedom in the joint fit to provide context for the achieved precisions.
Figure captions for the astrometric plots would benefit from explicit labels indicating the time baselines covered by each dataset.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of the manuscript, accurate summary of the joint orbital fit, and recommendation to accept. The significance statement correctly identifies the value of the model-independent masses for asteroseismology and the role of the 2023 periastron data.

Circularity Check

0 steps flagged

No circularity: dynamical masses derived from independent data via standard Keplerian relations

full rationale

The paper conducts a joint forward-model fit of independent observational datasets (radial velocities, relative astrometry, and absolute astrometry from Hipparcos/Gaia/ground catalogs over ~300 years) to constrain Keplerian orbital elements for the inner Aa-Ab and B-C subsystems, the wide A-BC orbit, and the center-of-mass parallax/proper motion. Component masses are subsequently computed from the fitted periods and semi-major axes using Newtonian gravity (Kepler's third law). This chain does not reduce to self-definition, fitted-input-as-prediction, or self-citation load-bearing steps, as the input measurements are external and the mass formulas are standard physics applied post-fit. No ansatzes, uniqueness theorems, or renamings are invoked for the core result. The derivation remains self-contained.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The masses follow from applying Kepler's third law to the fitted orbital elements of three nested Keplerian orbits; the only free parameters are the orbital elements themselves, which are constrained by the data rather than chosen ad hoc.

free parameters (1)

orbital elements (period, eccentricity, argument of periastron, inclination, node, semi-major axis) for Aa-Ab, B-C, andA
These are fitted simultaneously to the combined astrometric and RV time series; their values determine the masses via Kepler's law.

axioms (2)

domain assumption All orbits are purely Keplerian (Newtonian two-body motion with no significant third-body perturbations or relativistic corrections over the baseline).
Invoked in the forward-modeling framework to allow analytic orbit integration.
domain assumption The center-of-mass proper motion and parallax can be separated from the internal orbital motions.
Required for the joint fit of absolute astrometry.

pith-pipeline@v0.9.0 · 5695 in / 1561 out tokens · 46315 ms · 2026-05-10T15:54:45.015991+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

3 extracted references · 1 canonical work pages

[1]

Abt, H. A. 1973, ApJS, 26, 365 Adams, W. S. & Joy, A. H. 1923, ApJ, 57, 149 Andersen, M. F., Handberg, R., Weiss, E., et al. 2019, PASP, 131, 045003 Argelander, F. W. A. 1903, Eds Marcus and Weber’s Verlag, 0 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33 Baines, E. K., Armstrong, J. T., Schmitt, H. R., et al. 2018, ...

1973
[2]

M., Michalik, D., Brandt, T

Vol.1: Introduction and explanatory tables; Right Ascension 6h - 12h; Right Ascension 12h - 18h; Vol.5: Right Ascension 18h - 24h Brandt, G. M., Michalik, D., Brandt, T. D., et al. 2021, AJ, 162, 230 Campbell, W. W. 1913, Lick Observatory Bulletin, 229, 113 Campbell, W. W. 1928, Publications of Lick Observatory, 16, 1 Chaplin, W. J., Cegla, H. M., Watson,...

2021
[3]

∆α∗ orb,0(t) ∆δorb,0(t) # = β− a1 a

Article number, page 15 of 22 A&A proofs:manuscript no. main Fabricius, C., Luri, X., Arenou, F., et al. 2021, A&A, 649, A5 Feng, F., Butler, R. P., Jones, H. R. A., et al. 2021, MNRAS, 507, 2856 Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 Fouqu´ e, P., Moutou, C., Malo, L., et al. 2018, MNRAS, 475, 1960 Fricke, W., Kopff...

work page arXiv 2021