🔭 astrophysics

Envelope Inflation and outflow Driven by Energy Deposition in Massive Stars

Using one-dimensional hydrodynamical simulations of a $70 \, \rm M_{\odot}$ star, this study demonstrates that impulsive energy deposition within the stellar envelope triggers significant expansion and, under specific conditions, drives strong outflows capable of unbinding substantial portions of the outer envelope.

Original authors: Bhawna Mukhija, Amit Kashi

Published 2026-02-03

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Bhawna Mukhija, Amit Kashi

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 a massive star as a giant, glowing balloon filled with hot gas. Usually, this balloon is in a delicate balance: the heat pushing outward is perfectly matched by the weight of the gas pulling inward. This is how stars stay stable for millions of years.

But what happens if you suddenly dump a massive amount of extra energy into the middle of that balloon? That is exactly what this paper investigates. The authors, Bhawna Mukhija and Amit Kashi, used a computer simulation to act like a "thermal bomb" inside a 70-sun-sized star to see how it reacts.

Here is a simple breakdown of their findings:

The Setup: The "Thermal Bomb"

The researchers didn't actually blow up a real star. Instead, they created a digital model of a massive star that had already finished burning its main fuel. They then decided to inject a specific amount of energy into a narrow layer deep inside the star (about 23 times the size of our Sun away from the center).

They wanted to see two different scenarios:

The "Slow Stretch" (Hydrostatic): What if the star reacts slowly, like a rubber band being stretched?
The "Explosive Pop" (Hydrodynamic): What if the star reacts instantly, like a balloon popping?

Scenario 1: The Slow Stretch (Case 1)

In the first test, the computer was set up to prevent the star from actually shooting material out into space.

The Result: The star got hot and puffed up. It expanded slightly, becoming a bit bigger and cooler (like a red giant).
The Analogy: Imagine blowing air into a balloon very slowly. The balloon gets bigger and the rubber stretches, but it doesn't burst. The energy just made the star "inflate" without losing any of its skin.

Scenario 2: The Explosive Pop (Case 2A)

In the second test, they turned off the "safety brakes" and let the physics run wild. They injected the same energy but allowed the star to react dynamically.

The Result: The energy was so intense that it overcame the star's own gravity. The outer layers of the star were "unbound"—meaning they were no longer held down by the star's pull.
The Analogy: This is like over-inflating a balloon until it finally bursts. A stream of gas (an outflow) shoots out into space. The star didn't just expand; it started spitting material away at high speeds.

The Surprising Twist: More Energy = Less Expansion?

The authors tested what would happen if they dumped even more energy into the star (Case 2C). You might think more energy equals a bigger explosion and a bigger star.

The Finding: Surprisingly, the star with the most energy ended up being smaller, hotter, and less luminous than the one with less energy.
The Analogy: Think of a campfire. If you throw a small log on it, the fire grows steadily. But if you throw a massive, explosive firework onto the fire, it blasts the heat and fuel away so violently that the fire actually shrinks and cools down momentarily because the energy was carried away by the explosion itself.
Why? The massive energy input caused such a violent, fast outflow that the star lost its heat and mass so quickly that it couldn't stay puffed up. The energy escaped with the flying debris before it could make the whole star expand.

The Shape of the Explosion

The researchers also changed how they delivered the energy:

Narrow Blast: If the energy was dumped in a very thin, tight layer, it created a sharp, quick push.
Wide Blast: If the energy was spread out over a wider area, it created a more sustained, heavier flow of material.
Double Peaks: When they used the highest energy levels, the outflow didn't just happen once; it had a "double peak," meaning the star seemed to eject material in two distinct waves, like a hiccup that happens twice.

The Bottom Line

This paper is a "stress test" for massive stars. It shows that when you dump energy into a star:

It can either just puff up (if the physics are calm).
It can blow material away (if the physics are dynamic).
If you dump too much energy too fast, the star actually loses its "puffiness" because the explosion carries the heat away too quickly.

The authors emphasize that this is a simplified, idealized experiment. They aren't saying this is exactly how every real star behaves, but rather that this is the fundamental "hydrodynamic response" (the basic physics of fluid movement) of a star's skin when it gets a sudden, massive energy injection. It helps scientists understand the mechanics behind the violent eruptions we see in the universe, like those from "Luminous Blue Variables" (stars that are known for throwing massive tantrums).

Technical Summary: Envelope Inflation and Outflow Driven by Energy Deposition in Massive Stars

Problem Statement
Evolved massive stars frequently exhibit high-mass ejections and outflows, phenomena observed in Luminous Blue Variables (LBVs) and pre-supernova (pre-SN) evolution. While mechanisms such as super-Eddington winds and binary interactions are proposed drivers, the specific hydrodynamic response of a stellar envelope to impulsive energy deposition remains not fully understood. Previous studies have explored continuous energy deposition or parametric energy injection but often lack a detailed investigation into the transition from hydrostatic expansion to dynamical outflow driven by localized, impulsive energy inputs. This work addresses the gap in understanding how a massive star's envelope responds to sudden, localized energy deposition, specifically quantifying the conditions under which this energy unbinds material to form an outflow versus simply inflating the envelope.

Methodology
The authors employ the 1D hydrodynamical stellar evolution code MESA to simulate a $70 M_\odot$ star in the post-main sequence phase. The study utilizes an idealized "thermal bomb" approach, where energy is deposited instantaneously within a narrow Lagrangian mass range (corresponding to a radial width of $\sim 0.5 R_\odot$ ) located at a critical radius where the envelope's binding energy equals the gravitational energy of a hypothetical companion ( $R_{crit} \approx 23.21 R_\odot$ ).

The simulations compare two distinct physical regimes:

Case 1 (Hydrostatic): The model includes the mlt++ convection treatment and neglects enhanced electron scattering. This setup suppresses outflow generation, allowing the study of quasi-hydrostatic envelope expansion.
Case 2 (Hydrodynamic): The mlt++ treatment is disabled, and electron scattering opacity is artificially enhanced by a factor of 100. This configuration enforces a radiation-pressure-dominated envelope to examine the dynamical response, allowing material to become energetically unbound.

The study systematically varies three parameters across different cases:

Case 2A: Baseline hydrodynamic run with a heating region width of $0.5 R_\odot$ .
Case 2B: Hydrodynamic run with an increased heating region width of $0.8 R_\odot$ .
Case 2C: Hydrodynamic run with a higher deposited energy ( $E_{dep}$ ) compared to Case 2A, while keeping the heating region width constant.

The simulations track the evolution of total energy, density, mass-loss rates, and kinetic energy over a duration of 1.2 years, with a focus on the transition from bound to unbound material.

Key Results

Hydrostatic vs. Hydrodynamic Response: In the hydrostatic case (Case 1), energy deposition causes the envelope to expand significantly (radius increases by a factor of $\sim 1.09$ ) but remains gravitationally bound. No outflow is generated; the total energy in the outer layers remains negative. Conversely, in the hydrodynamic cases (Cases 2A, 2B, 2C), the deposited energy is sufficient to unbind outer layers, leading to the formation of a sustained outflow.
Outflow Onset: In the hydrodynamic models, the outflow initiates at a radius of $R_{out} \approx 58.1 R_\odot$ , where the total energy becomes positive. This radius is distinct from the energy deposition site ( $23\text{--}23.5 R_\odot$ ), indicating that the energy propagates outward to unbind material at larger radii.
Effect of Heating Region Width (Case 2B): Increasing the width of the heating region from $0.5 R_\odot$ to $0.8 R_\odot$ results in a higher total unbound mass and a higher average mass-loss rate ( $\dot{M}$ ). However, the broader heating region leads to a slightly less extreme stellar expansion compared to the narrower, more localized deposition. The broader region energizes a larger mass range, creating a more sustained, mass-loaded outflow, whereas the narrow region creates a steeper pressure gradient that drives a more impulsive but less massive ejection.
Effect of Deposited Energy (Case 2C): Increasing the deposited energy triggers a more violent and rapid outflow. Counterintuitively, this results in a final stellar configuration that is somewhat hotter, less luminous, and less expanded than the lower-energy Case 2A. The authors attribute this to the rapid conversion of deposited energy into the kinetic energy of the ejecta, which removes energy from the system before it can drive prolonged envelope inflation.
Double-Peaked Behavior: In the high-energy Case 2C, the outflow profile exhibits a double-peaked structure, suggesting multiple ejection episodes during the burst. This feature is absent in the lower-energy or broader-width cases.
Mass Loss and Velocities: The simulations yield mass-loss rates reaching up to $\sim 10^{-2} M_\odot \text{yr}^{-1}$ and outflow velocities of $\sim 4.5 \times 10^3 \text{ km s}^{-1}$ . The fraction of deposited energy carried away by the ejecta decreases over time (e.g., from $\sim 24\%$ at $t=0.4$ yr to negligible levels at $t=1.2$ yr in Case 2A).

Significance and Claims
The paper claims to establish a "base case" for the generic hydrodynamic response of massive star envelopes to impulsive energy input. The primary contribution is the demonstration that the nature of the response (inflation vs. outflow) is highly sensitive to the treatment of convection and opacity, as well as the spatial distribution and magnitude of the energy deposition.

The authors emphasize that their models are idealized and do not aim to reproduce specific observed events (such as the Great Eruption of $\eta$ Carina) or the specific physical triggers (e.g., binary dynamics, wave transport) in a self-consistent manner. Instead, the study isolates the fundamental hydrodynamic behavior of the envelope. The results suggest that while energy deposition can drive strong outflows, the efficiency of mass ejection versus envelope inflation depends non-linearly on the deposition parameters. Specifically, higher energy inputs may lead to more rapid energy loss through ejecta, limiting the subsequent expansion of the star. The authors note that these findings complement existing studies on steady, luminosity-driven mass loss by highlighting the distinct outcomes of impulsive energy injection.