🔭 astrophysics

Population synthesis of Be X-ray binaries in the Small Magellanic Cloud: angular momentum recycling and stable mass transfer

This study utilizes binary population synthesis to demonstrate that the observed population of Be X-ray binaries in the Small Magellanic Cloud is best reproduced by models featuring stable mass transfer, angular momentum recycling via tides, low natal kicks, and a strong propeller effect that suppresses accretion in wide systems.

Original authors: Víctor López Oller, Boyuan Liu, Michela Mapelli, Stefano Rinaldi, Cecilia Sgalletta, Julia Bodensteiner, Giuliano Iorio, Rebekka Schupp

Published 2026-04-30

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Víctor López Oller, Boyuan Liu, Michela Mapelli, Stefano Rinaldi, Cecilia Sgalletta, Julia Bodensteiner, Giuliano Iorio, Rebekka Schupp

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

Imagine the Small Magellanic Cloud (SMC) as a nearby, quiet cosmic neighborhood. In this neighborhood, there is a specific type of "celestial couple" called a Be X-ray Binary. Think of this couple as a fast-spinning, hot blue star (the "Be" star) wearing a swirling, messy skirt of gas (a disk), paired with a tiny, invisible, super-dense partner (a neutron star) that is constantly trying to steal some of that gas.

Astronomers have found about 100 of these couples in the SMC. Because this neighborhood is so uniform and close, it's the perfect "laboratory" to figure out exactly how these couples form, dance, and interact.

The authors of this paper ran thousands of computer simulations to answer a big question: What are the specific rules of physics that allow these couples to form and survive exactly as we see them?

Here is the breakdown of their findings using simple analogies:

1. The Dance of the Couple (Mass Transfer)

When the blue star loses gas, the neutron star tries to eat it.

The Old Idea: Scientists used to think the neutron star was a very picky eater, only taking a tiny bite and spitting the rest out.
The New Finding: The paper suggests the neutron star is actually a moderate eater. It takes about 60% of the gas offered to it. If it eats too little, the blue star partner ends up too small and faint. If it eats too much, the blue star becomes too heavy and bright. The "Goldilocks" zone of 60% matches the observations perfectly.

2. The Spin and the Orbit (Angular Momentum)

As the neutron star eats gas, it starts spinning faster, like a figure skater pulling in their arms.

The Problem: If it spins too fast, it might fly apart.
The Solution: The paper found that when the neutron star gets close to spinning too fast, it doesn't just throw the extra "spin energy" away. Instead, it acts like a rechargeable battery. It takes that extra spin and pumps it back into the couple's orbit, pushing them slightly further apart.
Why it matters: This "recycling" prevents the couple from spiraling too close together and crashing. Without this recycling, the simulations produced too many couples that were too close together, which doesn't match what we see in the sky.

3. The "Propeller" Effect (The Cosmic Brake)

Sometimes, the neutron star spins so fast that its magnetic field acts like a giant propeller.

The Mechanism: Imagine a fan spinning so fast it blows away anything trying to land on it. If the neutron star is spinning fast enough, it slams the gas away before it can be eaten.
The Result: This is crucial for wide couples (those far apart). If the propeller is too strong, these wide couples never light up in X-rays and remain invisible. The paper found that this "brake" must be very effective. If we assume the brake is weak, the computer predicts too many faint, long-period couples that we simply don't see.

4. The "Birth Jump" (Natal Kicks)

When the neutron star is born from a supernova explosion, it often gets a "kick" that sends it flying.

The Finding: The paper shows that for these couples to survive, the kick must be gentle. Most of the neutron stars in these systems received a very small push (less than 100 km/s). If the kick were too hard, it would fling the neutron star away from its partner, breaking the couple up forever.

5. The Stability of the Relationship

Finally, the paper looked at how stable the gas transfer is.

The Verdict: The relationship must be stable. The gas must flow smoothly from the blue star to the neutron star. If the flow gets chaotic (a "common envelope" phase), the two stars would likely merge into one and the couple would disappear. The data strongly suggests that these couples avoid the chaotic phase entirely.

The Big Picture

The authors compared their computer models against the real list of 102 couples in the SMC. They found that the only models that worked were the ones where:

The neutron star eats a moderate amount of gas (60%).
Excess spin is recycled back into the orbit to keep the couple apart.
The "propeller" effect strongly blocks gas from falling onto the neutron star in wide systems.
The neutron star gets a gentle birth kick.
The gas transfer is smooth and stable, never chaotic.

In short, the universe seems to have a very specific recipe for making these cosmic couples, and this paper has finally figured out the exact ingredients and the right amount of mixing required to get the result we see in the sky.

1. Problem Statement

Be X-ray binaries (BeXRBs) are critical laboratories for understanding binary interaction physics, including mass transfer (MT), angular momentum (AM) transport, and supernova (SN) natal kicks. However, previous population synthesis studies have largely relied on qualitative comparisons between models and observations, leaving key physical parameters loosely constrained.

The authors aim to quantitatively identify the specific combination of binary-evolution parameters that simultaneously reproduces:

The total observed number of BeXRBs in the Small Magellanic Cloud (SMC).
The joint distribution of orbital periods ( $P_{\text{orb}}$ ) and optical $V$ -band magnitudes ( $m_V$ ).

The SMC is chosen as the testbed due to its low metallicity ( $Z \approx 0.0035$ ), nearly complete and homogeneous BeXRB catalog (102 candidates), and well-constrained star formation history, which minimizes degeneracies associated with metallicity-dependent processes.

2. Methodology

Population Synthesis Code

The study utilizes SEVN (a binary population synthesis code), extended to specifically model BeXRBs based on previous work by Liu et al. (2024). The simulations evolve $10^6$ binary systems using stellar tracks appropriate for SMC metallicity.

BeXRB Identification Criteria

A simulated system is classified as a BeXRB if it satisfies:

Components: An O/B main-sequence star ( $M > 3 M_\odot$ ) and a neutron star (NS).
Decretion Disk: The donor star rotates rapidly ( $\omega = \Omega_{\text{rot}}/\Omega_{\text{crit}} > \omega_{\text{min}}$ ) and the orbital period is $> 1$ day (to allow disk formation).
Roche Geometry: The donor does not fill its Roche lobe at periastron (distinguishing them from Roche-lobe overflow X-ray binaries).

Physical Parameters Explored

The authors performed an extensive grid analysis varying the following key physical processes:

Mass Transfer Efficiency ( $f_{\text{MT}}$ ): The fraction of donor mass accreted by the companion (ranging from 0.1 to 1.0).
Angular Momentum (AM) Loss: Three prescriptions for non-accreted mass: Jeans-mode (from donor), Isotropic Re-emission (from accretor), and Circumbinary Ring.
AM Feedback (Recycling): How surplus AM is handled when the accretor approaches critical rotation:
- Elimination (EAM): Surplus AM is ejected (default in SEVN).
- Recycling (RAM): Surplus AM is recycled back into the orbit via tidal coupling.
RLOF Stability:
- Modified Hurley: Stability depends on the mass ratio $q$ and evolutionary phase.
- Always-Stable (ASMT): MT is assumed stable regardless of $q$ for relevant phases.
Natal Kicks: NS kicks drawn from a Maxwellian distribution (rescaled by ejecta mass), with a focus on low-velocity regimes.
Propeller Effect: A mechanism where a rapidly rotating NS magnetosphere prevents accretion if the corotation radius exceeds the magnetospheric radius. The study tests strict suppression ( $f_{\text{prop}}=0$ ) vs. residual accretion.

Statistical Framework

Instead of visual comparison, the authors employ a likelihood-based statistical framework:

Kernel Density Estimation (KDE): Constructed for the synthetic populations to generate probability density functions (PDFs) for $P_{\text{orb}}$ and the joint $P_{\text{orb}} - m_V$ space.
Resampling: Random subsamples of the synthetic population are drawn to match the observed sample size ( $N=102$ ).
Likelihood Calculation: The likelihood of the real observed data under the synthetic PDF is computed, yielding a $p$ -value.
Number Count Constraint: Models are penalized if they predict fewer than 75 systems (outside the $2\sigma$ confidence interval of the observed count, accounting for SFR uncertainty).
Calibration: The minimum spin threshold ( $\omega_{\text{min}}$ ) is tuned for each model to match the observed count of 102 systems, then checked for physical plausibility.

3. Key Contributions

Quantitative Statistical Methodology: The paper moves beyond qualitative "eye-balling" of distributions, establishing a rigorous likelihood framework that combines distribution shapes with population counts.
Joint Constraints: By simultaneously fitting $P_{\text{orb}}$ and $m_V$ , the study breaks degeneracies between AM transport and mass accretion efficiency that single-distribution analyses cannot resolve.
AM Recycling Mechanism: The study provides strong evidence that angular momentum must be recycled back into the orbit (via tides) when the accretor spins up, rather than being ejected, to reproduce the observed orbital architecture.
Propeller Effect Sensitivity: It demonstrates that the observable population of wide, low-accretion BeXRBs is highly sensitive to the treatment of the propeller effect and detection thresholds.

4. Key Results

Best-Fitting Physical Parameters

The models that successfully reproduce the SMC BeXRB population (within the 98% confidence interval) share the following characteristics:

Mass Transfer Efficiency: Moderately non-conservative, with $f_{\text{MT}} \approx 0.6$ .
- Reasoning: Lower efficiencies ( $f_{\text{MT}} \sim 0.2$ ) produce Be stars that are too faint (low mass) to match the observed $m_V$ distribution. Higher efficiencies produce overly massive/luminous stars.
Angular Momentum Feedback: Recycling (RAM) is required. Surplus AM from the accretor must be transferred back to the orbit via tidal coupling.
- Reasoning: This prevents excessive orbital shrinkage during MT, reproducing the observed lack of very short-period systems ( $P_{\text{orb}} \lesssim 10$ days) without requiring extremely inefficient mass transfer.
RLOF Stability: Always-Stable MT (ASMT) is favored.
- Reasoning: Models allowing unstable MT (leading to Common Envelope phases) fail to produce the correct total number of systems, as CE evolution tends to merge binaries or destroy the BeXRB progenitor.
Natal Kicks: Low kicks are preferred ( $\lesssim 100$ km s $^{-1}$ ). Approximately 40% of successful systems experience effectively zero kick.
Propeller Effect: A strong suppression of accretion ( $f_{\text{prop}} \approx 0$ $f_{prop} \approx 0$ ) is necessary.
- Reasoning: Strict suppression prevents wide, low-accretion systems from becoming detectable X-ray sources, keeping the synthetic population consistent with the observed lack of very long-period BeXRBs. Allowing residual accretion creates a "tail" of undetected long-period systems.
Spin Threshold: The minimum spin for a star to be a Be star is found to be $\omega_{\text{min}} \approx 0.84 - 0.91$ , consistent with observational limits for classical Be stars.

Distribution Analysis

Orbital Period: The scarcity of short-period systems is reproduced not by low mass transfer efficiency, but by the AM recycling mechanism which widens the orbit during the spin-up phase.
Optical Magnitude: The $m_V$ distribution is the primary driver for constraining $f_{\text{MT}}$ . The joint analysis rules out models that fit the period distribution but fail the magnitude distribution (e.g., low $f_{\text{MT}}$ models).

5. Significance and Implications

Binary Evolution Physics: The results suggest that BeXRB progenitors evolve through stable mass transfer phases where angular momentum is efficiently recycled, challenging models that rely heavily on Common Envelope evolution or highly non-conservative mass loss to explain orbital periods.
Gravitational Wave Progenitors: Since BeXRBs are potential progenitors for compact object binaries (NS-NS or NS-BH), constraining their formation pathways (specifically the stability of MT and kick distributions) improves predictions for gravitational wave event rates.
Observational Biases: The study highlights that the observed absence of very long-period BeXRBs may be partly due to the propeller effect suppressing X-ray emission below detection thresholds, suggesting future deeper X-ray surveys could reveal a hidden population of wide binaries.
Metallicity Dependence: The success of these parameters in the low-metallicity SMC provides a baseline for understanding binary evolution in different environments, with implications for the Milky Way population (where the best SMC models also qualitatively reproduce the Galactic period distribution).

In conclusion, the paper establishes that the SMC BeXRB population is best explained by a scenario involving stable mass transfer, moderate accretion efficiency, angular momentum recycling, and low natal kicks, with the propeller effect playing a crucial role in defining the observable sample.