🔭 astrophysics

Identifying heating processes in simulations with an entropy-based scheme: A single jet episode in a galaxy cluster

This paper introduces and validates an entropy-based scheme to systematically disentangle heating mechanisms in galaxy cluster simulations, revealing that while shock heating dominates early jet evolution, light jets primarily heat the intracluster medium through turbulent dissipation at later stages, whereas denser jets continue to rely on shock dissipation.

Original authors: Moun Meenakshi, Rainer Weinberger, Christoph Pfrommer, Thomas Berlok

Published 2026-03-31

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Moun Meenakshi, Rainer Weinberger, Christoph Pfrommer, Thomas Berlok

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: The "Thermostat" Problem

Imagine a giant, glowing cloud of gas surrounding a galaxy cluster. This gas is supposed to be hot, but it keeps trying to cool down. If it cools too much, it would rain down on the center, forming billions of new stars and turning the galaxy into a chaotic mess. This is the "Cooling Flow Problem."

Scientists believe that supermassive black holes in the center of these galaxies act as a thermostat. They shoot out powerful jets of energy (like a cosmic blowtorch) to reheat the gas and stop it from cooling. But here's the mystery: How exactly does that jet heat the gas? Is it like a shockwave hitting a wall? Is it like stirring a pot of soup? Or is it like mixing hot and cold water?

This paper is like a detective story. The authors built a new "heat detector" to figure out exactly how these jets warm up the galaxy cluster.

The Detective Tool: The "Entropy Tag"

To solve the mystery, the scientists needed a way to track heat without getting confused by the chaos of the simulation.

Think of the gas in the cluster as a crowd of people.

The Old Way: Scientists tried to guess how much heat was added by looking at how fast the people were moving or how loud the noise was. It was messy and often wrong.
The New Way (This Paper): The authors gave every bit of gas a special invisible "Entropy Tag."
- If the gas is just moving around smoothly, the tag stays the same.
- If the gas gets heated up (friction, shocks, or mixing), the tag changes permanently.

By tracking how much these tags change, the scientists can calculate exactly how much heat was added and, crucially, what caused it.

The Experiment: Light vs. Heavy Jets

The researchers ran computer simulations of these "cosmic blowtorches" (jets) shooting into the gas cloud. They tested two main types of jets, like comparing a firehose to a water balloon:

The "Light" Jet (The Firehose): This jet is very light and puffy. It's like a firehose spraying water. It doesn't have much momentum, so it can't punch a straight hole through the gas. Instead, it puffs up huge, wide bubbles.
The "Heavy" Jet (The Water Balloon): This jet is dense and heavy. It's like a solid water balloon. It has a lot of momentum and punches a narrow, straight tunnel through the gas.

What They Found: Two Different Heating Styles

1. The "Heavy" Jet: The Shockwave Heater

The heavy, dense jet acts like a bulldozer.

How it heats: It pushes the gas out of the way so hard that it creates massive shockwaves (like a sonic boom).
The Analogy: Imagine running through a crowd. If you are heavy and fast, you shove people aside violently. The friction of that shove creates heat.
Result: The heavy jet heats the gas mostly by creating strong shockwaves right in front of it. However, because it punches a narrow hole, it doesn't heat up the center of the cluster very well. It just shoots past it.

2. The "Light" Jet: The Stirring Heater

The light, puffy jet acts like a spoon stirring a pot.

How it heats: Because it's light, it gets pushed around easily. It inflates huge bubbles that rise up like hot air balloons. As these bubbles rise, they swirl the gas around, creating turbulence (chaotic swirling).
The Analogy: Imagine stirring a pot of soup. You aren't hitting the soup with a hammer; you are creating swirls and eddies. The friction of the swirling liquid heats the soup.
Result: The light jet creates massive bubbles that displace a huge amount of gas in the center. Even after the jet stops shooting, the bubbles keep rising and swirling, continuing to heat the gas for millions of years.

The "After-Dinner" Effect

One of the coolest discoveries is what happens after the jet turns off.

Heavy Jets: Once the bulldozer stops, the heat stops almost immediately.
Light Jets: Even after the firehose turns off, the huge bubbles it created keep floating and swirling. They continue to "stir" the gas and keep it warm for a long time (tens of millions of years). This is crucial because it means the galaxy doesn't need a jet shooting 24/7 to stay warm; a few good bursts are enough.

The Verdict

The paper concludes that lighter jets are actually better at keeping the galaxy's core warm than heavy ones.

Heavy jets are like a sniper: precise, strong, but they miss the big picture.
Light jets are like a sprinkler: they cover a wide area, mix everything up, and keep the whole "soup" warm for a long time.

By using their new "Entropy Tag" method, the scientists finally proved that turbulence and mixing (stirring) are just as important as shocks (hitting) when it comes to keeping galaxy clusters from freezing over.

1. Problem Statement

The paper addresses the cooling flow problem in galaxy clusters. Many clusters possess "cool cores" characterized by high central gas densities and short cooling times, which theoretically should lead to massive inflows of gas and vigorous star formation. However, observations show these inflows are suppressed. The leading explanation is Active Galactic Nucleus (AGN) feedback, where relativistic jets inject energy into the Intracluster Medium (ICM) to offset radiative cooling.

While simulations confirm AGN jets can provide sufficient energy, the specific physical mechanisms by which jet energy couples to the ICM and converts into heat remain uncertain. Proposed mechanisms include:

Turbulent dissipation (numerical or physical).
Shock heating (dissipative and adiabatic).
Cosmic Ray (CR) heating via Alfvén wave damping.
Mixing between jet plasma and ambient gas.
Anisotropic Braginskii viscosity.

Existing studies often lack a unified framework to disentangle these channels, particularly in the presence of magnetic fields and CRs. This paper aims to develop and validate a robust numerical method to identify and quantify these distinct heating processes in Magnetohydrodynamic (MHD) simulations.

2. Methodology

The authors implemented a comprehensive entropy-based heating estimator within the moving-mesh code Arepo.

A. Theoretical Framework

The method relies on the thermodynamic relation between entropy generation and irreversible heating. The authors distinguish between:

Irreversible Heating: Processes that generate entropy (viscous dissipation, shock dissipation, CR damping).
Reversible Heating: Adiabatic compression/expansion (which changes temperature but not entropy).

The total heating rate density ( $H$ ) is decomposed into:

Turbulent/Mixing Heating ( $H_{tu}$ ): Derived from the difference between the entropy evolved via the energy equation and entropy evolved via a passive advection equation.
Shock Heating ( $H_{sh}$ ): Calculated using an in-situ shock finder (Schaal & Springel 2015) to apply Rankine-Hugoniot jump conditions, separating dissipative (entropy-generating) and adiabatic (compression) components.
CR Heating ( $H_{cr}$ ): Energy transfer from CRs to thermal plasma via Alfvén wave damping.
Braginskii Viscosity: Anisotropic viscous heating along magnetic field lines.

B. Numerical Implementation

Entropy Tracers: Two passive scalars are evolved: the entropic function $K = P\rho^{-\gamma}$ $K = P ρ^{- γ}$ and its logarithm $s = \log(K)$ $s = lo g (K)$ .
- $K_{egy}$ and $s_{egy}$ are calculated from the conservative hydrodynamic update (including all heating/cooling).
- $K_{adv}$ and $s_{adv}$ are evolved via pure advection (no heating).
- The difference ( $\Delta K$ or $\Delta s$ ) over a time step $\Delta t$ yields the heating rate.
Robustness Check: To avoid numerical artifacts (e.g., negative heating at contact discontinuities), the authors only accept heating values where both the $K$ -based and $s$ -based estimators yield positive results.
Shock Separation: Shock zones identified by the shock finder are explicitly excluded from the turbulent heating calculation to prevent double-counting.

C. Validation Tests

The method was validated against analytical solutions and known physical behaviors in several test cases:

Turbulent Heating: Forced subsonic turbulence in a periodic box.
Braginskii Viscosity: Standing magnetosonic waves.
CR Heating: Damping of Alfvén waves excited by CR streaming.
Shock Heating: 1D shock tubes for weak ( $M \approx 2$ ) and strong ( $M \approx 10$ ) shocks.
Adiabatic Heating: Standing sound waves.

3. Key Contributions

Development of a Multi-Mechanism Diagnostic: A unified framework to simultaneously track turbulent, shock, CR, and viscous heating in MHD simulations, including magnetic fields and CRs.
Validation of Entropy Tracers: Demonstrated that using dual entropy tracers ( $K$ and $s$ ) effectively filters out numerical noise while accurately capturing physical heating.
Systematic Study of Jet Density: Applied the method to simulate single-epoch AGN jet outbursts with varying jet-to-ambient density ratios ( $\eta$ ) but fixed power, isolating the effect of density on heating mechanisms.

4. Results

A. Validation Outcomes

The entropy-based estimator showed excellent agreement with analytical solutions for turbulent dissipation, Braginskii viscosity, CR heating, and shock heating across a range of Mach numbers. The shock finder accurately recovered Mach numbers and heating rates for both weak and strong shocks.

B. Jet Evolution and Heating Mechanisms

Simulations of three jet types (Very-Light, Light, Intermediate-Density) revealed distinct evolutionary paths:

Jet Morphology:
- Light Jets: Form wide, quasi-spherical bubbles that displace a large fraction of the core gas. They propagate slower and reach smaller distances.
- Intermediate-Dense Jets: Propagate efficiently to large distances, forming narrow cocoons with less displacement of the central gas.
Dominant Heating Channels:
- Early Evolution: Shock heating (both dissipative and adiabatic) dominates for all jet types, regardless of density.
- Late Evolution (Light Jets): As light jets inflate large bubbles, turbulent dissipation and mixing become the primary heating mechanisms. The forward shock weakens rapidly.
- Late Evolution (Intermediate-Dense Jets): These jets continue to dissipate energy primarily via shocks (sustaining strong bow shocks) with minimal turbulent contribution.
Spatial Distribution:
- Inside the Cocoon: Heating is dominated by turbulent dissipation and mixing, even for dense jets.
- Outside the Cocoon: Heating is dominated by shocks (for dense jets) or acoustic compressions (for light jets, where shocks have weakened to $M \approx 1$ ).
Cooling Suppression:
- Light Jets are more effective at suppressing cooling in the cluster core ( $r < 30$ kpc) because they displace more gas and heat a larger volume via turbulence.
- Intermediate Jets fail to offset core cooling efficiently due to their collimated nature.
- Residual Heating: Even after jets are switched off, light jets leave a residual heating signature in the core for tens of Myr due to the slow dissipation of turbulent energy in the inflated bubbles.

5. Significance

Resolving the Cooling Flow Problem: The study demonstrates that light, low-density jets are potentially more efficient at regulating cool cores than previously thought, provided their energy is thermalized via turbulence and mixing rather than just shocks.
Mechanism Identification: The paper provides a definitive tool for future simulations to distinguish between "shock heating" and "turbulent heating," which is crucial for understanding the long-term thermal evolution of clusters.
Future Directions: The authors note that while this study used fixed-power jets, realistic clusters involve self-regulated feedback. The developed methodology is ready to be applied to self-regulated jet simulations to understand how jet density and power evolve over cosmic time to maintain thermal balance.

In summary, this work establishes a rigorous numerical diagnostic for AGN feedback, revealing that the density of the jet dictates not just its morphology, but the fundamental physical mechanism (shocks vs. turbulence) by which it heats the galaxy cluster.