Kinematics of Wolf-Rayet Stars in the LMC: Clues to Subtype Origins

Caden Burkhardt; Fiona Han; Mathieu Renzo; M. S. Oey; Natalia Ivanova

arxiv: 2603.10826 · v3 · pith:RI2ZFIE3new · submitted 2026-03-11 · 🌌 astro-ph.SR

Kinematics of Wolf-Rayet Stars in the LMC: Clues to Subtype Origins

Caden Burkhardt , Fiona Han , M. S. Oey , Natalia Ivanova , Mathieu Renzo This is my paper

Pith reviewed 2026-05-21 11:56 UTC · model grok-4.3

classification 🌌 astro-ph.SR

keywords Wolf-Rayet starsLarge Magellanic Cloudstellar kinematicsbinary evolutionproper motionsGaia astrometrymassive starsenvelope stripping

0 comments

The pith

Binary WC stars in the LMC show higher median transverse velocities and luminosities than single WC stars, indicating distinct formation channels.

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

The paper uses Gaia DR3 astrometry to measure transverse proper motions of Wolf-Rayet stars in the Large Magellanic Cloud and compares velocity distributions across subtypes. Very massive stars and some WNL populations display both slow and fast movers consistent with dynamical ejections on short timescales. For WC stars, binaries reach a median transverse velocity of 54 km/s and higher luminosities, while singles sit at 38 km/s; this pattern implies that single WC stars are not evolved descendants of the binary population. The authors conclude that mass transfer strips the envelopes of binary WC stars, whereas single WC stars arise through a separate process such as explosive mergers onto the shell-burning layer. Single and binary WNE stars both appear as ejected populations with single-peaked velocity distributions, pointing to yet another ejection mechanism.

Core claim

Binary WC stars appear faster with a median transverse velocity of 54 km/s and higher luminosities than single WC stars at 38 km/s. This difference indicates that single WC stars are not descendants of the binary population. Binaries are therefore stripped by mass transfer, while single WC stars originate from another process. The high velocities of binary WC stars align with predictions that lower-mass clusters produce fast dynamical ejections. Single WC and WN3/O3 stars share similarly high velocities around 38 km/s, possibly tied to their lower masses.

What carries the argument

Transverse velocity distributions separated by subtype and by single versus binary status for WC stars, derived from Gaia DR3 proper motions.

If this is right

Single WC stars must arise from an envelope-stripping channel other than binary mass transfer.
Binary WC stars are consistent with dynamical ejection from lower-mass clusters.
Dynamical ejections likely dominate the kinematics of very massive stars and classical WNL populations as well.
Both single and binary WNE stars form ejected populations with a common ejection mechanism.

Where Pith is reading between the lines

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

Models of massive-star envelope stripping will need separate pathways for single and binary WC stars to match the observed velocity split.
Extending the same Gaia-based velocity analysis to other nearby galaxies could test whether the single-binary distinction is universal.
The link between binary WC velocities and lower-mass clusters offers a direct test for cluster-dynamics simulations.

Load-bearing premise

The measured transverse velocities directly trace distinct ejection or formation mechanisms without significant contamination from measurement uncertainties, projection effects, or incomplete subtype and binary classification.

What would settle it

A larger sample with full three-dimensional space velocities that shows statistically indistinguishable velocity distributions for single and binary WC stars would falsify the claim of separate origins.

read the original abstract

We measure transverse proper motion velocities of LMC Wolf-Rayet (WR) stars using Gaia DR3 astrometry. The combined velocity distribution of WNh, O If*/WN, and WNL very massive stars ($>100\ M_\odot$; VMS) shows both slow, unejected objects ($v_\perp < 10$ $\rm km\ s^{-1}$) and stars dominated by fast, runaway velocities ($v_\perp > 24$ $\rm km\ s^{-1}$). This supports expectations that VMS ages are comparable to the dynamical ejection timescale ($\sim1.5$ Myr). These kinematics share similarities with those of lower-luminosity, classical WNh, O If*/WN, and WNL stars, as well as the SMC field OB stars, suggesting that dynamical ejections may also dominate these populations. In contrast, both single and binary WNE stars are ejected populations that show single-peaked velocity distributions, suggesting a different ejection mechanism(s). We speculate that single WNE stars might result from explosive mergers onto the shell-burning layer, thereby stripping the H envelope. Binary WC stars appear to be faster (median $v_\perp = 54$ $\rm km\ s^{-1}$) and have higher luminosities than singles (median $v_\perp = 38$ $\rm km\ s^{-1}$), suggesting that single WC stars are not descendants of the binaries. Thus, the binaries are probably stripped by mass transfer, while the WC singles likely originate from another process. The high velocities of binary WC stars are consistent with some predictions that lower mass clusters generate fast dynamical ejections. Single WC and WN3/O3 stars have ambiguous kinematics, but both show high $v_\perp$ (median $\sim 38$ $\rm km\ s^{-1}$), possibly linked to their lower masses.

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

Gaia DR3 velocities show binary WC stars faster than singles in the LMC, but thin statistics leave the separate-origin claim on shaky ground.

read the letter

The main takeaway is that binary WC stars show a higher median transverse velocity (54 km/s) than single WC stars (38 km/s), which the authors read as evidence that the singles do not descend from the binary population and must come from another channel such as dynamical ejection or mergers. They also note mixed slow and fast components in the very massive star group, consistent with dynamical ejection on Myr timescales, and single-peaked fast distributions for WNE stars. This is the concrete new content from the Gaia application. The work does well by breaking the sample into subtypes and reporting these median contrasts plus the luminosity offset for the WC stars. It turns an existing LMC WR catalog into a kinematic dataset that was not previously quantified this way. The data handling itself looks like a straightforward use of proper motions at LMC distance. The soft spots sit in the central WC claim. The medians are given without sample sizes per bin, uncertainty estimates, or any test for whether the 16 km/s difference is significant. LMC WR samples are typically modest, Gaia errors are a few km/s, and projection plus possible misclassification of single versus binary status can move medians noticeably. The stress-test note is right on this point: without those checks the inference that singles are not binary descendants rests on an untested contrast. The rest of the subtype comparisons are more qualitative and less load-bearing. This paper is for people working on massive-star evolution, cluster dynamics, and the channels that produce different WR subtypes. A reader who needs updated observational constraints from Gaia on ejection versus stripping would find the velocity numbers useful as a starting point. It deserves a serious referee because the astrometric application is new and the topic is relevant, even if the statistical presentation and robustness checks need work. I would send it to peer review rather than desk reject.

Referee Report

2 major / 2 minor

Summary. The manuscript measures transverse proper motions of Wolf-Rayet stars in the LMC with Gaia DR3 astrometry. It reports that the combined velocity distribution of WNh, O If*/WN, and WNL very massive stars shows both slow and fast runaway components consistent with dynamical ejection on ~1.5 Myr timescales, while single and binary WNE stars exhibit single-peaked ejected distributions. For WC stars, binary systems display a higher median transverse velocity (54 km s^{-1}) and higher luminosities than single WC stars (38 km s^{-1}), leading the authors to conclude that single WC stars are not descendants of the binary channel and that binaries are stripped by mass transfer while singles arise via a separate process (possibly explosive mergers). Single WC and WN3/O3 stars show similarly high velocities (~38 km s^{-1}).

Significance. If the reported median velocity contrast and its interpretation survive detailed statistical scrutiny, the work supplies direct kinematic evidence distinguishing binary mass-transfer stripping from other envelope-removal channels for WC stars at LMC metallicity. This would tighten constraints on massive-star evolution models and on the relative contributions of dynamical ejection versus binary interaction in low-metallicity environments.

major comments (2)

[Abstract and §3] Abstract and §3 (velocity results): the central claim that binary WC stars (median v_⊥ = 54 km s^{-1}) are not the progenitors of single WC stars (median v_⊥ = 38 km s^{-1}) rests on the reported median difference, yet no sample sizes, bootstrap uncertainties, Kolmogorov-Smirnov or Mann-Whitney p-values, or outlier-robustness tests are provided. With LMC WC samples typically modest, even a few high-velocity objects or Gaia proper-motion errors of a few km s^{-1} can shift the median by tens of percent; without these statistics the contrast cannot be shown to be load-bearing for the distinct-channel conclusion.
[§4] §4 (interpretation): the inference that the velocity offset directly traces distinct ejection mechanisms assumes negligible contamination from projection effects, incomplete binary/single classification, and distance uncertainties. No quantitative assessment of these systematics (e.g., Monte-Carlo projection of 3-D velocities or sensitivity to adopted LMC distance) is presented, weakening the link between observed v_⊥ distributions and the proposed formation scenarios.

minor comments (2)

[Figures 2-4] Figure captions and text should explicitly state the number of stars contributing to each median and the velocity-error floor adopted from Gaia DR3 at LMC distance.
[§2] Notation for transverse velocity (v_⊥) is used consistently, but the conversion from proper motion to km s^{-1} (including the precise LMC distance and solar-motion correction) should be given once in §2.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive feedback on our manuscript. Their comments have prompted us to enhance the statistical presentation and systematic analysis. Below we provide point-by-point responses to the major comments.

read point-by-point responses

Referee: [Abstract and §3] Abstract and §3 (velocity results): the central claim that binary WC stars (median v_⊥ = 54 km s^{-1}) are not the progenitors of single WC stars (median v_⊥ = 38 km s^{-1}) rests on the reported median difference, yet no sample sizes, bootstrap uncertainties, Kolmogorov-Smirnov or Mann-Whitney p-values, or outlier-robustness tests are provided. With LMC WC samples typically modest, even a few high-velocity objects or Gaia proper-motion errors of a few km s^{-1} can shift the median by tens of percent; without these statistics the contrast cannot be shown to be load-bearing for the distinct-channel conclusion.

Authors: We agree with the referee that additional statistical details are necessary to robustly support the claimed velocity difference between binary and single WC stars. In the revised version, we will report the sample sizes for each group, provide bootstrap-derived uncertainties on the median velocities, and include the results of a two-sample Kolmogorov-Smirnov test (or Mann-Whitney U test) to quantify the significance of the difference. We will also perform outlier-robustness checks by recomputing medians after removing the highest-velocity objects and discuss the impact of Gaia proper motion uncertainties. These changes will be incorporated into Section 3 and referenced in the abstract. revision: yes
Referee: [§4] §4 (interpretation): the inference that the velocity offset directly traces distinct ejection mechanisms assumes negligible contamination from projection effects, incomplete binary/single classification, and distance uncertainties. No quantitative assessment of these systematics (e.g., Monte-Carlo projection of 3-D velocities or sensitivity to adopted LMC distance) is presented, weakening the link between observed v_⊥ distributions and the proposed formation scenarios.

Authors: We acknowledge the importance of addressing potential systematics in our interpretation. For projection effects, we will add a Monte Carlo analysis assuming random orientations to estimate the distribution of 3D velocities from the observed transverse components. Regarding distance uncertainties, we will test the sensitivity of our velocity calculations to variations in the adopted LMC distance (e.g., ±2% and ±5%). For binary/single classification, we will expand the discussion to include the reliability of current classifications based on available spectroscopy and note any potential for contamination, perhaps with a simple sensitivity test. These quantitative assessments will be added to Section 4 to better link the observations to the formation scenarios. revision: yes

Circularity Check

0 steps flagged

No circularity: direct Gaia astrometry measurements and median comparisons are self-contained observational results.

full rationale

This is an observational kinematics study that derives transverse velocities directly from Gaia DR3 proper motions and parallaxes for LMC WR stars. Median v_perp values are computed from the measured data for binary vs. single WC subsamples without any parameter fitting, model ansatz, or self-referential prediction that reduces to the same inputs by construction. The interpretive suggestion regarding distinct formation channels is an inference from the observed contrasts, not a derivation that loops back to fitted quantities or self-citations. No equations, uniqueness theorems, or load-bearing self-citations appear in the provided text; the analysis relies on external astrometric data and is therefore self-contained against independent benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on standard assumptions about LMC distance and the conversion of proper motion to transverse velocity; no free parameters are fitted to the WR data itself and no new entities are postulated.

axioms (1)

domain assumption LMC distance and geometry allow direct conversion of Gaia proper motions to transverse velocities without large systematic errors.
Invoked when reporting v_perp values in km/s from astrometric data.

pith-pipeline@v0.9.0 · 5892 in / 1287 out tokens · 38616 ms · 2026-05-21T11:56:58.395954+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Binary WC stars appear to be faster (median v⊥ = 54 km s−1) and have higher luminosities than singles (median v⊥ = 38 km s−1), suggesting that single WC stars are not descendants of the binaries.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We speculate that single WNE stars might result from explosive mergers onto the shell-burning layer

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.