🔭 astrophysics

A Kinematic Study of Wolf-Rayet Stars at the Galactic Center I: Binary Candidates and Constraints on the Binary Fraction

This kinematic study of Wolf-Rayet stars within 0.5 pc of the Galactic center, utilizing the longest time-baseline (1994–2024) radial velocity data, identifies five binary candidates and infers a binary fraction of 0.56±0.18 (rising to 0.69±0.17 when combined with photometric data), a result consistent with both previous studies of young Galactic center stars and field Wolf-Rayet populations.

Original authors: Rory O. Bentley, Tuan Do, Andrea Ghez, Devin Chu, Anna Ciurlo, Abhimat K. Gautam, Zoë Haggard, Matthew W. Hosek Jr., Kelly Kosmo O'neil, Rebecca Lewis-Merrill, Gregory D. Martinez, Anna Pusack, Shoko

Published 2026-06-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Rory O. Bentley, Tuan Do, Andrea Ghez, Devin Chu, Anna Ciurlo, Abhimat K. Gautam, Zoë Haggard, Matthew W. Hosek Jr., Kelly Kosmo O'neil, Rebecca Lewis-Merrill, Gregory D. Martinez, Anna Pusack, Shoko Sakai, Jessica R. Lu, Mark R. Morris, Keith Matthews

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 Detective Story

Imagine the center of our galaxy, the Milky Way, as a bustling, crowded city square. In the very middle of this square sits a supermassive black hole (Sgr A*), acting like a giant, invisible whirlpool. Surrounding this whirlpool is a neighborhood of massive, short-lived stars called Wolf-Rayet (WR) stars. These are the "rock stars" of the galaxy: they are incredibly bright, huge, and burn through their fuel so fast they live fast and die young.

For a long time, astronomers knew these stars existed, but they didn't know a crucial detail about their social lives: Are they lonely, or are they in relationships?

This paper is a detective report that answers that question. The team of astronomers spent decades (from 1994 to 2024) watching these stars to see if they are part of binary systems—pairs of stars orbiting each other like a cosmic dance couple.

The Method: Listening to the Cosmic Dance

How do you tell if a star has a partner if you can't see the partner? You watch the star's movement.

Think of a star in a binary system like a person on a merry-go-round. If you are standing still and watching them, they don't just spin in place; they bob up and down and sway side-to-side as they orbit their partner.

The Clue: As the star moves toward us, its light gets slightly "bluer." As it moves away, the light gets slightly "redder." This is called a radial velocity change.
The Challenge: These stars are also blowing massive "winds" (streams of gas) that can wiggle their light, making it look like they are moving when they aren't. It's like trying to hear a whisper in a hurricane.

The researchers used the Keck and Gemini telescopes (the best "ears" in the galaxy) to listen to the light from 27 of these stars over a 30-year period. They measured the speed of the stars with extreme precision to see if the "wobble" was caused by a hidden partner or just the star's own windy atmosphere.

The Findings: Who is Dancing?

After analyzing the data, the team found:

Five Suspects: Out of the 27 stars they watched, 5 showed clear signs of having a partner.
- Three were already known to be couples (like IRS 16SW and IRS 16NE).
- Two were brand new discoveries identified in this paper: S8-181 and IRS 13E4.
The "Wobble" Proof: For the new candidates, the stars moved back and forth so dramatically that it couldn't be explained by their own winds or the gravity of the central black hole. It had to be a partner tugging on them.
- Analogy: Imagine a lighthouse beam. If the light just flickers because of a storm, that's one thing. But if the entire lighthouse is swinging back and forth on a giant pendulum, you know something heavy is attached to it. That's what they saw.

The Big Conclusion: The "Dating" Rate

The most important number in this paper is the Binary Fraction. This is the percentage of stars that are in relationships.

The Result: The team calculated that about 56% of these Wolf-Rayet stars near the galactic center are in binary pairs.
The Combined Result: When they combined their new data with previous studies of other young stars in the same area, the number goes up to about 69%.

What does this mean?
It means that in this extreme neighborhood near the black hole, being single is the exception, not the rule. Most of these massive stars are dancing in pairs.

Why Does This Matter?

The paper uses this "dating rate" to test theories about how stars are born in such a crazy place.

The Theory: Some scientists thought the black hole's gravity was so strong it would rip any star-forming clouds apart, making it impossible for stars (or pairs of stars) to form there.
The Reality: The fact that so many stars are in pairs suggests that star formation did happen there, likely in a swirling disk of gas that managed to overcome the black hole's pull. The high number of "couples" supports the idea that these stars formed in a chaotic, crowded disk rather than being scattered there individually.

Summary

Who: A team of astronomers led by Rory Bentley and Tuan Do.
What: They watched 27 massive stars near the center of our galaxy for 30 years.
How: They measured the stars' speeds to detect the "wobble" caused by invisible partners.
Result: They found 2 new star couples and confirmed that roughly two-thirds of these massive stars are in binary systems.
Takeaway: The center of our galaxy is a very social place where massive stars prefer to form and live in pairs, helping us understand how stars can survive near a supermassive black hole.

Technical Summary: A Kinematic Study of Wolf-Rayet Stars at the Galactic Center I: Binary Candidates and Constraints on the Binary Fraction

Problem Statement
The origin of the young stellar population at the Galactic Center, particularly the Wolf-Rayet (WR) stars located within 0.5 pc of the supermassive black hole Sgr A*, remains poorly understood. Standard models suggest that tidal shear from Sgr A* should prevent star formation at these distances, necessitating alternative formation scenarios such as in-situ formation within an accretion disk. The binary fraction of these stars serves as a critical diagnostic for testing these formation models; for instance, star formation in an accretion disk predicts a higher binary fraction due to high fragmentation rates. Furthermore, the binary fraction influences the dynamical evolution of the cluster through mergers and evaporation, and binary motion can bias kinematic measurements used to determine stellar membership in specific structures. While the binary fraction of the general young stellar population in the Galactic Center has been estimated via photometry (>71%), spectroscopic constraints specifically for the WR population have been limited by small sample sizes (typically ≤8 stars) and short time baselines.

Methodology
This study presents the longest time-baseline kinematic study (1994–2024) of WR stars in the Galactic Center to date.

Sample Assembly: The authors constructed a sample of 27 WR stars and candidates within 0.5 pc of Sgr A*, comprising 18 WN stars and 9 WC stars. This sample triples the number of WR stars included in previous spectroscopic binary searches in the region.
Observations: The study utilizes 110 new integral field unit observations taken between 2006 and 2024, primarily using the OSIRIS spectrograph on the W. M. Keck Observatory (101 observations) and the NIFS spectrograph on the Gemini North telescope (9 observations). These were combined with 70 previously published radial velocity ( $v_z$ ) measurements from the literature.
Data Reduction and Analysis: Spectra were reduced using standard pipelines (OSIRIS DRP and Nifty4Gemini). Radial velocities were derived using cross-correlation against empirical templates constructed from the observed spectra themselves to minimize sensitivity to specific atmospheric parameters. For the double-lined binary IRS 16SW, $v_z$ was determined via Gaussian fitting of the 2.112 µm He I absorption line.
Error Budget: Statistical uncertainties were derived via Monte Carlo simulations. Systematic uncertainties arising from wind-driven spectral variability were added in quadrature: 15 km/s for WN stars (including Ofpe/WN9 and hybrids) and 3 km/s for WC stars, based on prior characterizations of wind variability.
Binary Identification: Candidates were identified as stars exhibiting a $v_z$ variation between any two epochs with a significance greater than $3\sigma$ , after correcting for the orbital acceleration of the stars around Sgr A*.
Statistical Inference: To infer the intrinsic binary fraction ( $f_{bin}$ ), the authors performed 1000 simulations. These simulations randomly assigned binary status to the 27 stars based on a trial $f_{bin}$ , generated binary physical properties (period, mass ratio, eccentricity) based on distributions observed in field WR stars (Sana et al. 2012; Dsilva et al. 2022, 2023), and simulated $v_z$ measurements matching the observational cadence and uncertainties of the real data. The number of simulated binaries detected with $>3\sigma$ significance was compared to the actual number of detections to derive the posterior distribution of $f_{bin}$ .

Key Contributions and Results

New Binary Candidates: The study identifies two new binary candidates among the WR stars: S8-181 (WN8) and IRS 13E4 (WC9).
- S8-181: Shows a $v_z$ variation of $139 \pm 42$ km/s ( $3.4\sigma$ ) between 2022 and 2024 spectra. The variation exceeds wind-driven variability limits and gravitational acceleration from Sgr A*.
- IRS 13E4: Shows a $v_z$ variation of $71 \pm 15$ km/s ( $4.8\sigma$ ) over a 14-year baseline. The variation exceeds the expected wind variability for WC stars ( $\sim 10$ km/s) and Sgr A* acceleration.
Confirmed Binaries: The study confirms three previously known binaries: IRS 16SW, IRS 16NE, and S4-258.
Binary Fraction: Based on the detection of 5 binaries (3 known, 2 new) out of the 27 stars in the sample, and accounting for experimental sensitivity via simulation, the authors infer a spectroscopic binary fraction for Galactic Center WR stars of $f_{bin} = 0.56 \pm 0.18$ .
Combined Constraints: When combined with the photometric binary fraction results from Gautam et al. (2024) for the broader young stellar population, the study derives a combined binary fraction of $f_{bin} = 0.69 \pm 0.17$ .

Significance and Claims
The paper claims that the derived binary fraction of $0.56 \pm 0.18$ for Galactic Center WR stars is consistent with:

Previous photometric studies of young stars in the Galactic Center.
The binary fractions of field WR stars (specifically WN stars) and O-stars in young Milky Way clusters.
The suggestion that the binary fraction increases with radial distance from Sgr A*, as the measured fraction is significantly higher than the upper limit of $f_{bin} \le 0.17$ found for B-stars within 0.04 pc of the black hole (Chu et al. 2023).

The authors emphasize that while their result supports the high binary fraction predicted by accretion disk formation models, the inference relies on assumptions regarding the distribution of binary physical properties (period, mass ratio) which are not independently constrained for the Galactic Center population. The study also notes that dust-producing WC stars, often associated with colliding wind binaries, did not universally show velocity variations in this sample, leading the authors to conservatively exclude dust production as a sole indicator of binarity in their analysis. The work represents the largest spectroscopic sample of Galactic Center WR stars to date, providing a robust kinematic constraint on the binary population in this extreme environment.