🔭 astrophysics

EWOCS-V: Is Wd1-72 a recent post-interaction WR+O binary?

This paper proposes that the unique nebular morphology of the Wd1-72 binary in Westerlund 1 is best explained by a hydrodynamic simulation of a recent (~10 kyr) post-Roche Lobe Overflow event, suggesting the system is a hydrogen-free, multi-episode mass-transfer binary that serves as a key benchmark for studying the evolution of gravitational-wave progenitors.

Original authors: C. J. K. Larkin, J. Mackey, H. Jin, A. A. C. Sander, B. Reville, K. Anastasopoulou, M. Andersen, A. Bayo, J. J. Drake, E. K. Grebel, M. G. Guarcello, T. J. Haworth, V. M. Kalari, R. R. Lefever, F. Naj

Published 2026-03-31

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: C. J. K. Larkin, J. Mackey, H. Jin, A. A. C. Sander, B. Reville, K. Anastasopoulou, M. Andersen, A. Bayo, J. J. Drake, E. K. Grebel, M. G. Guarcello, T. J. Haworth, V. M. Kalari, R. R. Lefever, F. Najarro, B. W. Ritchie, E. Sabbi

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: A Cosmic "Aftermath" Mystery

Imagine you are looking at a massive, crowded city (the Westerlund 1 star cluster). In the middle of this city, there is a very famous, powerful celebrity star called Wd1-72. This star is a "Wolf-Rayet" star, which is essentially a massive, aging star that has blown off its outer layers and is now burning very hot and bright.

Recently, astronomers pointed the James Webb Space Telescope (JWST) at this star and saw something weird. Instead of a clean, empty space around it, there was a messy, structured cloud of gas and dust. It looked like two different things happening at once:

Comet Tails: Some gas was streaming away from the star, looking like the tails of comets.
Spherical Droplets: Other gas was clumped together in round, ball-like shapes, floating closer to the center of the cluster.

The big question was: How did this happen? Was it a violent explosion from the star itself? Or was it the result of a messy interaction with a neighbor?

The Theory: A Cosmic "Haircut" Gone Wrong

The authors of this paper propose a specific story: Wd1-72 is a star that just got a very messy haircut from its partner.

Here is the analogy:

The Binary System: Wd1-72 isn't alone; it's in a tight dance with another massive star (an O-type supergiant). They are so close they are practically hugging.
The "Roche Lobe Overflow" (RLOF): Imagine the two stars are holding hands. If they get too close, the "gravity hand" of the heavier star starts pulling the outer skin (atmosphere) off the other star. This is called Roche Lobe Overflow. It's like one star is stealing a massive amount of hair from the other.
The "Non-Conservative" Mess: Usually, when stars steal mass, they keep it. But in this case, the authors suggest the mass was thrown away (non-conservative). It was like the star tried to give the hair to its partner, but the partner was too full, so the hair was flung out into space in a chaotic spray.

The Simulation: Recreating the Mess

To test this idea, the scientists built a digital movie (hydrodynamic simulation) of what happens when this "hair spray" meets the environment.

The Wind: The cluster is full of powerful winds blowing from other stars, acting like a strong breeze in a stadium.
The Collision: When the "hair spray" (the stolen mass) comes out, the strong cluster wind hits it.
The Result:
- The "Comet Tails": The wind catches the loose material and drags it away, stretching it out into long tails, just like a kite tail or a comet.
- The "Droplets": Where the wind hits the material head-on (coming from the direction of the cluster center), it creates a pressure zone. This pressure squashes the gas into round, ball-like droplets, similar to how water beads up on a windshield.

The computer simulation looked almost exactly like the photos taken by the JWST. This suggests the "messy haircut" theory is likely correct.

Why This Matters: The "Second Time" Clue

There is a twist in the story. The star Wd1-72 has no hydrogen left on its surface.

The Analogy: If you shave a head once, you still have some stubble. To get a completely bald head, you usually have to shave it multiple times, or shave it very aggressively.
The Science: The authors argue that a single "haircut" (mass transfer event) wouldn't strip the star clean enough to make it a Wolf-Rayet star this quickly. This suggests Wd1-72 has been through this messy interaction more than once. It's a veteran of cosmic violence.

Why Should We Care?

This isn't just about one weird star. It helps us understand the "family trees" of the universe.

Gravitational Waves: When massive stars die, they can turn into black holes. If two black holes merge, they create ripples in space-time called gravitational waves.
The Connection: To get two black holes close enough to merge, they usually need to go through these messy mass-transfer phases first. By studying Wd1-72, we are essentially looking at a "training ground" for the systems that will eventually create those gravitational wave signals.

Summary

The paper argues that the strange, comet-like clouds around the star Wd1-72 are the leftover debris from a recent, violent mass-transfer event between two close stars. The star was "stripped" of its outer layers, and the cluster's wind shaped that debris into the beautiful, complex patterns we see today. It's a snapshot of a star system in the middle of a dramatic transformation, potentially setting the stage for a future black hole merger.

1. Problem Statement

The evolutionary origin of Wolf-Rayet (WR) stars at Solar metallicity remains debated. While single-star evolution suggests massive O stars strip their hydrogen envelopes via steady winds or Luminous Blue Variable (LBV) eruptions, observed binary period distributions of WR subtypes challenge this. Specifically, WN (nitrogen-rich) binaries often have short periods, while WC (carbon-rich) binaries have long periods, implying different formation pathways.

The study focuses on Wd1-72, a WN7b + O supergiant binary system in the young massive cluster Westerlund 1 (Wd 1). JWST observations reveal complex nebulosity around Wd1-72:

Upstream (towards the cluster core): Quasi-spherical droplets.
Downstream (away from the core): Cometary tails pointing away from the star.
This morphology is unexpected for a cluster of this age, as feedback should have dispersed natal gas. The authors investigate whether this structure results from a recent, intense mass-loss event via non-conservative Roche Lobe Overflow (RLOF) interacting with the cluster wind, rather than a giant LBV eruption.

2. Methodology

The authors employed a two-pronged approach combining stellar evolution modeling and hydrodynamic simulations.

A. Stellar Evolution Modeling (MESA)

Code: Used the Modules for Experiments in Stellar Astrophysics (MESA, v10398).
Track Selection: Identified a specific binary evolution track (Solar metallicity, initial masses $M_1 = 28.18 M_\odot$ , $M_2 = 19.73 M_\odot$ , initial period 3.98 d) consistent with Wd1-72's age, luminosity, and period.
Mass-Loss Adjustment: To match the observed mass-loss rate ( $\dot{M} \approx 9 \times 10^{-5} M_\odot \text{yr}^{-1}$ ) of Wd1-72, the authors artificially increased the primary star's wind mass-loss rate by a factor of 4.7 immediately following the second RLOF event.
Outcome: This adjusted track predicts the system reaches a surface hydrogen abundance of $\sim 1\%$ only $\sim 20$ kyr post-RLOF, supporting a scenario of multiple RLOF episodes (required to strip H completely before the droplets disperse).

B. Hydrodynamic Simulations (PION)

Code: Used the magneto-hydrodynamic code PION.
Setup: 2D cylindrical domain ( $r, z$ ) with static mesh refinement (finest resolution $\Delta x \approx 1.22 \times 10^{15}$ cm).
Physics:
- Cluster Wind: Constant velocity ( $1500 \text{ km s}^{-1}$ ), density ( $1 \times 10^{-24} \text{ g cm}^{-3}$ ), and pressure.
- Source Term: Injected RLOF material from the binary (velocity $20 \text{ km s}^{-1}$ ) into the cluster wind.
- Radiation: Included radiative heating/cooling and full photoionization by the cluster's EUV field; magnetic fields were neglected.
- Dust: Post-processed with TORUS to calculate thermal dust emission at $11 \mu\text{m}$ , assuming a 1% dust mass fraction in RLOF material.

3. Key Results

A. Morphological Reproduction
The hydrodynamic simulation successfully qualitatively reproduced the observed JWST/MIRI F1130W morphology:

Upstream ( $+z$ ): RLOF material forms Rayleigh-Taylor (RT) "clumps" that appear as quasi-spheroidal droplets due to the stagnation region where the ejecta collides with the cluster wind.
Downstream ( $-z$ ): The cluster wind sweeps the material, stretching the RT clumps into cometary tails pointing away from the star.
Fragmentation: The fragmentation of the RLOF shell into droplets is driven by RT instabilities as the WR wind expands into the slower RLOF ejecta.

B. Physical Constraints & Timescales

Timescale: The simulation suggests the system is only $\sim 10$ kyr post-RLOF. This short timescale is necessary because the droplets disperse rapidly ( $\lesssim 50$ kyr).
Hydrogen Abundance: The model requires the system to be in a second (or subsequent) RLOF episode to explain the hydrogen-free (WN7) nature of the star, as a single RLOF event would likely leave residual hydrogen if the system is this young.
Mass-Loss Rate: The required mass-loss rate for the observed morphology is $\sim \text{few} \times 10^{-4} M_\odot \text{yr}^{-1}$ . This is incompatible with a giant LBV eruption (which typically requires $\sim 0.1 M_\odot \text{yr}^{-1}$ ), strongly favoring the RLOF scenario.

C. Discrepancies

Spatial Scale: The simulated droplets appear roughly twice as far from the star as observed. This is attributed to the 2D nature of the simulation (fragmentation is more efficient in 3D) and resolution limits.
Tail Distribution: The simulation produces uniform tails, whereas observations show a non-uniform distribution (sparser in the upper part). This may be due to the assumption of spherical symmetry in the outflow or unknown system orientation relative to the cluster wind.

4. Key Contributions

First Hydrodynamic Model of Post-RLOF WR Nebulae: Demonstrated that non-conservative RLOF mass-loss, fragmented by a WR wind and swept by a cluster wind, can reproduce the specific "droplet + cometary tail" morphology seen in Wd1-72.
Constraint on LBV vs. Binary Interaction: Provided physical arguments (ram-pressure balance and required $\dot{M}$ ) to rule out a giant LBV eruption as the origin of the nebula, reinforcing the binary interaction hypothesis.
Evolutionary Benchmark: Identified Wd1-72 as a potential "fossil" of a recent, unstable mass-transfer event, occurring only $\sim 10$ kyr ago.
Dust Modeling: Estimated dust emission properties, confirming that the temperature gradient (warmer near the cluster core, cooler downstream) matches JWST observations.

5. Significance

Binary Evolution: If Wd1-72 is indeed a recent post-RLOF system, it provides empirical evidence that unstable mass transfer occurs in short-period WN binaries. This challenges single-star formation models and highlights the critical role of binary interactions in WR evolution.
Gravitational Wave Progenitors: Understanding mass transfer efficiency and angular momentum loss in these systems is crucial for modeling the formation of double compact object binaries (black hole/neutron star mergers). Wd1-72 offers a rare, observable snapshot of the immediate aftermath of mass transfer.
Future Observations: The paper outlines a roadmap for follow-up, including integral field spectroscopy to measure velocities and composition, and multi-epoch imaging to measure proper motions. These observations will constrain the transition timescale from sgB[e] to WR and the cluster wind velocity.

In conclusion, the paper posits that Wd1-72 is a young, hydrogen-depleted binary that recently underwent a violent mass-transfer event, creating a unique laboratory for studying the physics of non-conservative mass loss and the early stages of binary evolution leading to gravitational wave sources.