A comprehensive grid of massive binary evolution models for the Galaxy - Surface properties of post-mass transfer stars

This paper presents a comprehensive grid of over 38,000 massive binary evolution models using MESA to characterize the diverse surface properties, evolutionary stages, and chemical abundances of post-mass transfer stars, providing a crucial tool for distinguishing binary interaction products from single stars and improving supernova progenitor identification.

The Gist

Imagine the universe as a giant, bustling cosmic dance floor. For a long time, astronomers thought most of the "stars" (the dancers) were solo acts, spinning gracefully on their own. But we now know that most massive stars are actually dancing in pairs, holding hands (or rather, sharing gravity) in binary systems.

This paper is like a massive, detailed instruction manual for what happens when these star-dancers get too close and start interacting. The authors, a team of astrophysicists, ran over 38,000 computer simulations to see how these pairs evolve, specifically focusing on what happens after they swap mass (like one dancer giving their energy to the other).

Here is the breakdown of their findings, translated into everyday language:

1. The Setup: A Massive Simulation Lab

Think of the authors as chefs in a giant kitchen. Instead of cooking soup, they cooked 38,000 different star recipes.

• The Ingredients: They started with stars ranging from 5 to 100 times the mass of our Sun.
• The Recipe: They used a super-complex computer program (called MESA) that tracks everything: how fast they spin, how they lose mass like steam from a kettle, how they pull on each other, and even the nuclear chemistry happening deep inside them (tracking elements from Hydrogen all the way to Aluminum).
• The Goal: They wanted to know: "If two stars swap mass, how do they look and taste (chemically) compared to a star that stayed single?"

2. The Main Event: The "Mass Transfer" Dance

In a binary system, the bigger star (the "Donor") often expands as it ages. If it gets too close to its partner, it spills its outer layers onto the smaller star (the "Gainer").

• The Donor (The one losing mass): Imagine a balloon being deflated. As the Donor loses its outer skin, it gets stripped down.

  • The Surprise: Instead of just fading away, some of these stripped stars turn into Blue or Yellow Super Giants. They stay hot and bright for a long time.
  • The "LB-1" Mystery: There was a famous star system (LB-1) that confused astronomers. They thought it had a black hole, but this paper suggests it's actually a "stripped" star that lost its outer layers but is still burning helium in its core. The authors' models fit this perfectly.

• The Gainer (The one gaining mass): Imagine a star eating a massive buffet. It swells up and spins faster.

  • The Surprise: These "gluttonous" stars often become Blue or Yellow Super Giants too, but they are luminous and hotter than a single star of the same mass would be. They essentially "cheat" death by staying in the blue zone longer.

3. The "Chemical Fingerprint"

This is the most exciting part. When stars swap mass, they don't just change size; they change their chemical makeup.

• The Single Star: A solo star is like a quiet library. It changes slowly. If it spins fast, it mixes its own chemicals, bringing some "cooked" material (Nitrogen) to the surface, but it's a slow process.
• The Binary Star: A binary star is like a smoothie blender.
  • The Gainer: It swallows the Donor's "cooked" material (which is rich in Nitrogen and Helium). Suddenly, the Gainer's surface is a smoothie of Nitrogen and Helium. It also spins so fast it mixes everything up instantly.
  • The Donor: It gets stripped down to its core, revealing material that was deep inside, which is also rich in Nitrogen but poor in Hydrogen.
  • The Boron Clue: The authors found a special "detective" element: Boron. In single stars, Boron disappears slowly. In binary interactions, it gets wiped out almost instantly. If astronomers see a star with high Nitrogen but zero Boron, they know: "Aha! This star was in a binary system!"

4. The Grand Finale: Supernovae

When these stars die, they explode as supernovae. The paper shows that the "type" of explosion depends heavily on whether the star was a solo act or a dancer.

• Solo Stars: Usually explode as Red Super Giants (Type IIP).
• Binary Stars: Because they lost their hydrogen "skin" (envelope) during the dance, they explode differently.
  • Some become Type IIb (like the famous SN 1993J), which have very little hydrogen left.
  • Some become Type Ibc (stripped completely).
  • The authors' models predict that the "stripped" stars from binary systems are the perfect match for these observed explosions.

5. Why This Matters

Think of this paper as a Rosetta Stone for astronomers.

• Before: We saw a strange star and didn't know if it was a weird solo act or a binary survivor.
• Now: We have a massive grid of predictions. If we see a star with specific colors, temperatures, and chemical fingerprints (like high Nitrogen and low Boron), we can look up this "menu" and say, "Ah, this star must have been a mass-gainer in a binary system!"

In Summary:
The universe is full of star couples. When they interact, they don't just change their orbits; they fundamentally change their identity, their chemical taste, and how they die. This paper provides the ultimate guidebook to recognizing these "cosmic couples" and understanding the dramatic stories they tell through their light and chemistry.

Technical Summary

Here is a detailed technical summary of the paper "A comprehensive grid of massive binary evolution models for the Galaxy – Surface properties of post-mass transfer stars" by Jin et al. (2025).

1. Problem Statement

Massive stars frequently evolve in binary systems, where interactions such as mass transfer, tidal forces, and mergers significantly alter their evolutionary paths compared to single stars. While previous studies have focused on predicting the final fates of massive stars (e.g., supernovae, gravitational wave sources) or late-stage phenomena (e.g., Be X-ray binaries), there is a lack of comprehensive models covering earlier evolutionary stages (main sequence and post-main sequence) with detailed surface abundance predictions.

Specifically, existing grids often:

• Limit nuclear networks to helium and CNO abundances.
• Focus on subsets of interacting binaries rather than a full parameter space.
• Lack detailed predictions for the full range of hydrogen-burning nucleosynthesis products (from light elements to aluminum).

This gap hinders the ability to observationally distinguish binary interaction products (mass donors and mass gainers) from single stars using spectroscopic data, particularly for main-sequence and supergiant phases.

2. Methodology

The authors computed a comprehensive grid of 38,430 binary evolution models using the MESA (Modules for Experiments in Stellar Astrophysics) code (version 10398).

Key Physical Assumptions & Inputs:

• Stellar Physics:
  • Nuclear Network: A new, extensive network for hydrogen burning tracking stable isotopes up to Silicon and radioactive $^{26}$Al. It includes pp-chains, CNO-cycles, Ne-Na, and Mg-Al cycles, allowing the tracking of surface abundances from H to Al.
  • Convection & Mixing: Standard mixing-length theory ($\alpha_{MLT}=1.5$), mass-dependent overshooting (linearly increasing from 0.1 to 0.3), and thermohaline mixing.
  • Rotation: Initial rotation set to 20% of critical velocity. Rotational mixing implemented with diffusion coefficients ($f_c=1/30$, $f_\mu=0.1$) and the Spruit-Tayler dynamo.
  • Mass Loss: Combined prescriptions for hot (Vink et al.) and cool (Nieuwenhuijzen & de Jager) stars, plus Wolf-Rayet rates for stripped stars.
• Binary Physics:
  • Mass Transfer: Implicit schemes for Case A (core H-burning) to ensure Roche lobe containment; explicit schemes (Kolb & Ritter) for post-main sequence.
  • Accretion Efficiency: Determined self-consistently based on the accretor's rotation. If the accretor spins up to critical rotation, excess mass is lost from the system, carrying away specific orbital angular momentum.
  • Tides: Implemented following Detmers et al. (2008) with timescales based on Zahn (1977).
• Grid Parameters:
  • Primary Mass ($M_{1,i}$): 5 to 100 $M_\odot$.
  • Initial Mass Ratio ($q_i$): 0.1 to 0.95.
  • Initial Orbital Period ($P_i$): 0.3 to 5011 days (covering the full range of interacting systems).
  • Metallicity: Solar (Galactic) abundances (Asplund et al. 2021).

3. Key Contributions

• Comprehensive Grid: The largest and most detailed grid of massive binary models to date for the Galaxy, covering the full binary parameter space and evolutionary stages from ZAMS to pre-supernova.
• Extended Nucleosynthesis: Unlike previous grids limited to CNO/He, this work tracks the full suite of hydrogen-burning products (Li, Be, B, C, N, O, F, Ne, Na, Mg, Al, Si), providing a rich chemical fingerprint for binary interactions.
• Differentiation Strategy: The paper establishes specific surface abundance signatures (e.g., Boron depletion, Nitrogen enhancement, Isotopic ratios) that can distinguish post-mass transfer stars from single stars, even when the binary companion is no longer visible.

4. Key Results

A. Evolution in the Hertzsprung-Russell Diagram (HRD)

• Core Helium Burning (CHeB) Phase:
  • Mass Gainers: Many evolve as Blue (BSG) or Yellow Supergiants (YSGs) during CHeB, whereas single stars of the same mass would be Red Supergiants (RSGs). This occurs particularly in Case A mass transfer where significant mass is accreted, keeping the star hot.
  • Mass Donors: Exhibit a wide diversity. While many are stripped to become Wolf-Rayet stars or Helium stars, partially stripped donors remain as RSGs, YSGs, or BSGs. These partially stripped stars are often overluminous compared to single stars of the same mass.
• Supernova Progenitors:
  • Binary interactions produce progenitors for diverse supernova types (Type IIP, IIL, IIb, Ibc).
  • SN 1993J (Type IIb): The models successfully reproduce the progenitor's properties (YSG with a small H-envelope) as a partially stripped mass donor.
  • LB-1 System: The models suggest the stripped B-type star in LB-1 could be a core-helium-burning mass donor, rather than a star in a thermal contraction phase.

B. Surface Abundance Signatures

• Nitrogen (N) and Boron (B):
  • Mass Gainers: Show significantly stronger Nitrogen enhancement and Boron depletion compared to single stars of similar luminosity. This is due to the accretion of CNO-processed material and accretion-induced mixing.
  • Diagnostic Power: High N and low B are strong indicators of binary interaction, distinguishing them from single stars where such levels would require extremely high initial rotation velocities ($>400$ km/s), which are rare.
• Isotopic Ratios ($^{12}$C/$^{13}$C and $^{16}$O/$^{17}$O):
  • Mass Gainers: Exhibit lower $^{12}$C/$^{13}$C ratios than single stars due to contamination by CNO-processed material.
  • Mass Donors: Show a wide range of ratios depending on the degree of envelope stripping.
• Light Elements (Li, Be, B):
  • Mass gainers show distinct depletion patterns. For instance, some RSG mass gainers retain detectable Lithium, which is unexpected for single stars of that mass, potentially explaining anomalous observations in the Galaxy.

C. Comparison with Observations

• The models align well with observed progenitors of specific supernovae (e.g., SN 1993J, iPTF13bvn).
• The grid explains the existence of "runaway" stars with high N and low B (e.g., $\zeta$ Ophiuchi) as mass gainers.
• The models address the "RSG Problem" (the lack of observed high-luminosity RSG progenitors) by showing that binary interactions can produce high-luminosity RSGs that are actually mass gainers or partially stripped donors.

5. Significance and Future Applications

• Population Studies: The grid serves as a foundational tool for population synthesis, helping to estimate the occurrence rates of binary interaction products.
• Supernova Modeling: Provides a robust theoretical baseline for predicting the diversity of supernova progenitors and their pre-explosion properties (envelope mass, chemical composition).
• Observational Guidance: Offers specific chemical tracers (N, B, isotopic ratios) for astronomers to identify binary interaction products in large spectroscopic surveys (e.g., VFTS, IACOB, BLOeM).
• Physics Constraints: Future comparisons between these predictions and observations will help constrain uncertainties in stellar physics, such as mass transfer efficiency, rotational mixing, and angular momentum transport.

In conclusion, this work bridges the gap between theoretical binary evolution and observational spectroscopy, providing a critical resource for understanding the life cycles of massive stars in the Galaxy.