🔭 astrophysics

Radiatively Cooled Binary Mass Transfer: Flow Structure, Luminosities, and L2 Outflows Across Mass Transfer Rates

This paper presents hydrodynamical simulations of radiatively cooled binary mass transfer that reveal how high mass transfer rates ( $\gtrsim 10^{-3} M_\odot$ /yr) drive significant L2 outflows and generate luminous optical transients ( $\sim 10^5 L_\odot$ ), whereas lower rates result in mostly conservative accretion.

Original authors: Peter Scherbak, Wenbin Lu, Jim Fuller

Published 2026-05-19

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Peter Scherbak, Wenbin Lu, Jim Fuller

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 two stars dancing in a tight embrace, locked in a gravitational waltz. One is a massive, bloated "donor" star, and the other is a smaller "accretor" star. Sometimes, the donor star gets so big that it spills its outer layers of gas over onto its partner. This is called Mass Transfer.

In the past, scientists thought this spillage was like a gentle stream of water filling a bucket: the smaller star catches everything, and the system stays balanced. But this new study by Scherbak, Lu, and Fuller suggests that when the spill is really fast, things get messy, chaotic, and spectacularly bright.

Here is what they found, explained through simple analogies:

1. The "Overflowing Bathtub" Effect

Think of the smaller star (the accretor) as a bathtub. When the donor star pours gas into it, the water level rises.

Slow Spill: If the donor pours slowly, the bathtub catches almost all the water. The system is "conservative" (nothing is wasted).
Fast Spill: If the donor pours water in at a massive rate (like a firehose), the bathtub overflows. The water doesn't just splash out; it builds up a huge, swirling wave that crashes over the edge of the tub and shoots out into the room.

In the stars, this "room" is space. When the spill is fast enough, the gas builds up a thick, puffy disk around the accretor until it spills over a specific "leak" in the gravitational field called the L2 point. This creates a massive, donut-shaped wind of gas (a Circumbinary Outflow) that escapes the system entirely.

2. The "Speed Limit" of the Leak

The researchers ran computer simulations to see what happens at different "pouring speeds" (Mass Transfer Rates).

The Tipping Point: They found a critical speed. If the donor pours slower than about 1,000 times the mass of the Sun per year, the smaller star catches most of the gas. But if the pouring speed jumps higher (above 1,000 solar masses per year), the system suddenly becomes "non-conservative." The smaller star can't hold it all, and a huge amount of gas escapes.
The Escape Route: The escaping gas doesn't just fly off randomly. It flows out through a specific "door" (the L2 point) on the far side of the smaller star. It's like water finding the lowest point in a dam and rushing out.

3. The "Spinning Top" and the "Brake"

When gas escapes, it carries something very important: Angular Momentum (the spin energy that keeps the two stars orbiting each other).

The Metaphor: Imagine two ice skaters holding hands and spinning. If one skater throws a heavy backpack away, they spin faster or slower depending on how they throw it.
The Finding: The gas escaping through the L2 point carries away a huge amount of spin energy—much more than if it just fell onto the smaller star. This acts like a powerful brake on the binary system, likely causing the two stars to spiral closer together over time.

4. The "Glowing Fireworks"

The paper also looked at how bright these events would be.

Friction and Heat: As the gas swirls around the smaller star and crashes into the escaping stream, it heats up due to friction (shocks), just like rubbing your hands together creates heat.
The Light Show:
- Slow Spills: The glow comes mostly from the gas swirling right next to the smaller star. It's bright, but manageable.
- Fast Spills: When the spill is massive, the escaping wind itself becomes incredibly hot and bright. The gas in the escaping "donut" wind glows as brightly as the gas near the star.
- The Colors: Depending on how fast the gas is pouring, the light changes color. Fast spills glow in ultraviolet and blue light (very hot), while slower spills glow in red and infrared (cooler).

5. Why This Matters (According to the Paper)

The authors note that while we can't see these specific "fast spill" events easily right now, they likely happen in nature.

Real-World Examples: They mention systems like SS 433 and W Serpentis, which seem to be spewing out gas in ways that match their simulations.
Supernova Clues: They suggest that before a massive star explodes as a supernova, it might go through a phase of rapid mass transfer, creating a cloud of gas around it. When the star finally explodes, that gas interacts with the blast, changing how we see the explosion.

Summary

In short, this paper uses computer simulations to show that when a massive star dumps gas onto a partner too quickly, the partner can't catch it all. Instead, the gas builds up and shoots out through a specific gap, carrying away the system's spin energy and creating a massive, glowing cloud of hot gas. The faster the dump, the more gas escapes, and the brighter the cosmic fireworks become.

Technical Summary: Radiatively Cooled Binary Mass Transfer

Problem Statement
Mass transfer (MT) in binary star systems is a critical process in stellar evolution, preceding the formation of gravitational wave sources and core-collapse supernovae. A major uncertainty in modeling these systems is the degree to which MT is conservative—specifically, whether the accretor retains all transferred material or if mass escapes the system, carrying away angular momentum (AM). Previous work (e.g., Scherbak et al. 2025, "Paper I") utilized hydrodynamical simulations of rapid, stable MT under adiabatic assumptions, finding that material flows out through the outer Lagrange point ( $L_2$ ), forming circumbinary outflows (CBOs). However, Paper I neglected radiative cooling, which is expected to become significant at moderate MT rates ( $\lesssim 10^{-2} M_\odot \text{yr}^{-1}$ ). The lack of cooling in prior models left open questions regarding the transition between conservative and non-conservative MT regimes, the specific AM carried by escaping mass, and the observational properties (luminosity and temperature) of these transients.

Methodology
The authors perform 3D Newtonian hydrodynamical simulations using the PLUTO code to model stable MT between a $10 M_\odot$ donor star and a $5 M_\odot$ point-mass accretor in a fixed circular orbit. Unlike Paper I, this study incorporates an approximate radiative cooling model.

Cooling Prescription: The cooling power per volume is estimated assuming photons escape perpendicular to the orbital plane. The model interpolates opacity tables (solar metallicity) from MESA and solves for gas temperature by balancing gas and radiation pressure. The cooling term acts to reduce internal energy, subject to safety floors to prevent unphysical values (e.g., a temperature floor of 2000 K in high-density regions and a pressure floor in low-density regions).
Simulation Grid: The authors vary the orbital separation $a$ to probe a range of physical MT rates ( $\dot{M}_{\text{phys}}$ ) from $\sim 10^{-5}$ to $\sim 10^{-1} M_\odot \text{yr}^{-1}$ . Four primary simulations are conducted: lowest_mdot ( $2.0 \times 10^{-5}$ ), mid_mdot ( $2.2 \times 10^{-3}$ ), high_mdot ( $2.3 \times 10^{-2}$ ), and highest_mdot ( $1.5 \times 10^{-1}$ ).
Constraints: The simulations do not include explicit viscosity, viscous heating, or detailed radiative transport. The accretor is represented by a Plummer softening potential. The donor star is driven to overflow its Roche lobe via heat injection on a Kelvin–Helmholtz timescale.

Key Results

Flow Structure and Conservativeness:
- Mass flows from the donor into a disk around the accretor. For high MT rates ( $\gtrsim 10^{-3} M_\odot \text{yr}^{-1}$ ), significant equatorially concentrated outflows occur through $L_2$ .
- The fraction of mass retained by the accretor, $\beta$ , decreases sharply as the MT rate increases. For $\dot{M} \gtrsim 10^{-2} M_\odot \text{yr}^{-1}$ , MT is highly non-conservative ( $\beta \approx 0$ ). For lower rates ( $\dot{M} \lesssim 10^{-3} M_\odot \text{yr}^{-1}$ ), MT is mostly conservative ( $\beta \gtrsim 0.9$ ), though a small amount of material still escapes via $L_2$ .
- The transition to highly non-conservative MT occurs approximately when the accretion luminosity at the outer edge of the accretor's disk exceeds the Eddington limit. The authors note their threshold for non-conservative behavior is higher (by a factor of $\sim 2-5$ ) than semianalytic predictions by Lu et al. (2023).
Angular Momentum and Energetics:
- In all cases, the specific AM ( $h_{\text{loss}}$ ) of the outflowing gas asymptotes to a value approximately equal to that of the $L_2$ point ( $h_{L2}$ ), regardless of the MT rate. This implies that even small outflows efficiently extract AM from the binary.
- The asymptotic radial velocity of the outflow increases as the MT rate decreases (from $\sim 10\%$ of orbital velocity at high rates to $\sim 50\%$ at lower rates). However, in physical units, velocities remain in the range of $30-60 \text{ km s}^{-1}$ for most simulations.
- Bernoulli parameters indicate that outflows at high MT rates are bound, while those at lower rates are marginally unbound.
Luminosities and Observational Appearance:
- Gas cooling luminosity ( $L$ ) and temperature ( $T$ ) increase with MT rate.
- At low MT rates, cooling is dominated by the accretor's disk. At high rates ( $\sim 10^{-1} M_\odot \text{yr}^{-1}$ ), the CBO contributes significantly, doubling the total luminosity.
- Peak cooling luminosities reach $\sim 10^{39} \text{ erg s}^{-1}$ for the highest MT rate and $\sim 10^{38} \text{ erg s}^{-1}$ for the high MT rate.
- Effective temperatures ( $T_{\text{eff}}$ ) range from $\sim 5 \times 10^4 \text{ K}$ (UV) in the accretor's disk for high rates to $\sim 10^2 \text{ K}$ (mid-infrared) in the CBO for lower rates. The authors caution that dust formation is possible but difficult to predict definitively within their simulation domain, as temperatures required for condensation may only be reached at radii larger than their grid extends.

Significance and Claims
The paper claims to delineate the regime where radiative cooling significantly alters the dynamics of stable MT compared to adiabatic models. By systematically varying MT rates, the authors identify a critical threshold ( $\sim 10^{-3}$ to $10^{-2} M_\odot \text{yr}^{-1}$ ) where the system transitions from conservative to highly non-conservative, driven by the super-Eddington nature of the outer disk.

The authors assert that their results provide necessary inputs for population synthesis calculations investigating the stability of MT and the orbital evolution of binary compact objects prior to merger. They highlight that the specific AM carried by $L_2$ outflows is consistently high, which will impact the long-term evolution of the binary orbit. Furthermore, the calculated luminosities and temperatures offer predictions for the observational appearance of transients associated with rapid MT, such as those potentially seen in systems like SS 433, W Serpentis, or V1309 Sco. The authors modestly note that their results are limited to a single mass ratio and solar metallicity, and that future work is required to model dust formation at larger radii and to incorporate viscous accretion physics.