🔭 astrophysics

Reevaluating thermal instability in a uniform plasma: an extended analysis of instability domains

This study reevaluates thermal instability in a uniform, non-magnetic plasma by extending the dispersion analysis of Field (1965) and Waters & Proga (2019) to clarify instability domains across all wavelength regimes, correct misconceptions regarding adiabatic perturbations, and refine the classification of thermal modes based on the Field length.

Original authors: Varsha Felsy, Ramon Oliver, Jaume Terradas, Amanda Stricklan, Timothy Waters, James A. Klimchuk

Published 2026-03-03

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Varsha Felsy, Ramon Oliver, Jaume Terradas, Amanda Stricklan, Timothy Waters, James A. Klimchuk

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 universe is filled with a giant, invisible soup of super-hot gas called plasma. This soup exists in stars, between galaxies, and even in our own Sun's atmosphere. Usually, this soup is in a delicate balance: it's being heated up by energy sources (like magnetic fields) and cooled down by radiating light away.

This paper is like a weather report for that cosmic soup. It asks a simple question: What happens when a little bit of this soup gets a tiny bit colder or hotter than its neighbors? Does it stay that way, or does it spiral out of control?

Here is the breakdown of the research using everyday analogies:

1. The Big Idea: The "Snowball Effect"

The scientists are studying Thermal Instability. Think of it like a snowball rolling down a hill.

Stable: If you push a snowball slightly, it might wobble but stop.
Unstable: If the hill is steep enough, a tiny nudge makes the snowball roll faster, gathering more snow, getting bigger, and rolling uncontrollably.

In space, if a patch of hot gas cools down slightly and starts radiating more heat away as it gets colder, it creates a runaway effect. It cools down faster, gets denser, and eventually collapses into cold, dense clumps (like coronal rain falling on the Sun or cold clouds forming in space).

2. The Old Map vs. The New GPS

For decades, scientists used a map created by a physicist named G. B. Field in 1965 to predict where these "snowballs" would form. Later, in 2019, a team led by Waters and Proga tried to simplify this map into a single, easy-to-read chart based on one number (let's call it the "R-Score").

The Problem: The 2019 map was a bit like a GPS that worked perfectly on a straight highway but got confused when you hit a bumpy road. It ignored a crucial factor: Thermal Conduction.

The Analogy: Imagine trying to keep a cup of coffee hot. If you have a lid (conduction), the heat spreads out and keeps the whole cup warm. Without a lid, the center gets cold while the edges stay hot. The 2019 map mostly ignored the "lid."

The New Paper's Job: The authors (Felsy, Oliver, et al.) took that 2019 map and added the "lid" back in. They wanted to see if the map still worked when heat could flow around freely.

3. The Three "Modes" of the Soup

When the soup gets disturbed, it reacts in three main ways, like a drum being hit:

The Thermal Mode (The Slow Sneeze): The gas just gets hotter or colder without moving much. It's a slow, steady change.
The Acoustic Mode (The Sound Wave): The gas ripples back and forth like a sound wave. It's fast and bouncy.
The "Catastrophic Cooling" Mode: A specific type of slow sneeze where the gas just gives up, cools down instantly, and collapses.

The paper maps out exactly when the soup chooses to be a "sneeze" (thermal) or a "ripple" (acoustic), and when it decides to be stable (calm) or unstable (chaotic).

4. The "Field Length" (The Size of the Spoon)

The most important discovery in this paper is about the Field Length.

The Analogy: Imagine you are stirring a pot of soup with a giant spoon.
- If you stir with a tiny spoon (short wavelengths), the heat spreads out so fast that it stops any cold spots from forming. The soup stays smooth.
- If you use a giant spoon (long wavelengths), the heat can't spread fast enough to save the cold spot. The cold spot grows, and the instability takes over.

The authors found that the old 2019 map only works if your "spoon" is small enough. If the spoon gets too big (meaning the heat can't conduct fast enough to stop the instability), the map breaks down, and the rules change.

5. The "Coronal Rain" Mystery

The paper ends with a specific application to the Sun. There is a phenomenon called Coronal Rain, where hot plasma in the Sun's atmosphere cools down and falls like rain.

Some recent theories suggested this happens because of a specific "fastest-growing" instability (the Catastrophic Cooling mode).
The Twist: The authors checked the actual physics of the Sun's atmosphere and found that, for the Sun, the "R-Score" is usually too high for that specific "fastest-growing" mode to be the winner.
The Verdict: The Sun's rain might not be caused by the "fastest" instability everyone thought it was. It's likely caused by a different, more complex version of the instability that peaks at a specific size, not at the very largest size.

Summary: Why Should You Care?

This paper is like rewriting the instruction manual for a cosmic engine.

It fixes errors in a popular guide used by astronomers.
It explains why some parts of the universe form cold clouds while others stay hot.
It helps us understand why the Sun rains "plasma" and how stars and galaxies evolve over time.

In short: They took a complex mathematical puzzle, added a missing piece (heat conduction), and showed us that the picture of the universe is a little more complicated—and a little more interesting—than we thought.

1. Problem Statement

Thermal instability is a fundamental mechanism driving the formation of multi-phase structures (e.g., cold clouds, filaments, coronal rain) in astrophysical plasmas. While the foundational theory was established by G. B. Field (1965) and recently revisited by T. Waters & D. Proga (2019), several gaps and misconceptions remain:

Classification Limitations: Field's original instability diagrams were complex and dependent on multiple wavenumber parameters. Waters & Proga (2019) simplified this using a single parameter $R$ , but their classification assumed negligible thermal conduction ( $\lambda_F \to 0$ ) and contained errors regarding critical wavelength transitions.
Misconceptions on Perturbation Types: There is a common belief that purely isobaric (constant pressure) or isochoric (constant density) perturbations exist across all scales. The authors aim to clarify the existence of these modes, particularly in the limits of very short and very long wavelengths.
Catastrophic Cooling: The recent proposal of "catastrophic cooling" as a distinct instability mechanism (Waters & Stricklan 2025) requires re-evaluation within the standard linear thermal instability framework.
Solar Application: Understanding the specific instability regimes relevant to solar coronal loops and the formation of coronal rain.

2. Methodology

The study employs a linear stability analysis of an infinite, homogeneous, non-magnetic hydrodynamic plasma in thermal equilibrium.

Governing Equations: The authors utilize the linearized equations of gas dynamics, incorporating a general heat-loss function $L = \Lambda - H$ (radiative cooling minus heating) and Spitzer's thermal conductivity ( $\kappa \propto T^{5/2}$ ).
Dispersion Relation: They analyze the cubic dispersion relation derived by Field (1965) and reformulated by Waters & Proga (2019). The relation is expressed in terms of the growth rate $n$ and wavenumber $k$ .
Dimensionless Parameters:
- $N_{\rho0}$ and $N_{p0}$ : Dimensionless forms of the isochoric and isobaric instability criteria.
- $R = N_{p0} / (\gamma N_{\rho0})$ : A key parameter determining the mode structure.
- $\lambda_F$ : The Field length (characteristic scale where cooling balances conduction).
- $\lambda_{th}$ : The thermal wavelength ( $c_s t_{cool}$ ).
Analytical & Numerical Approach:
- Asymptotic Analysis: Derivation of analytical solutions for the limits $\lambda \to 0$ (short wavelength) and $\lambda \to \infty$ (long wavelength).
- Critical Wavelengths: Solving the dispersion relation to find critical wavenumbers where mode types (thermal vs. acoustic) transition.
- Parameter Space Exploration: Systematic plotting of growth rates ( $n_R$ ) and frequencies ( $n_I$ ) across different values of $R$ , $N_{\rho0}$ , $N_{p0}$ , and the ratio $\lambda_F/\lambda_{th}$ .

3. Key Contributions

Refined Instability Classification (Table 1): The authors correct and expand the classification of instability domains proposed by Waters & Proga (2019). They introduce a comprehensive table that accounts for:
- The inclusion of thermal conduction (finite $\lambda_F$ ).
- The stabilization of modes at short wavelengths.
- The merging of previously distinct regimes (specifically $1/9 < R < 1/3$ ).
Correction of Critical Wavelengths: They demonstrate that the critical wavelengths ( $\lambda_c, \lambda_\pm$ ) separating thermal and acoustic modes depend on $\lambda_F$ . The simplified formulas from Waters & Proga are only valid when $\lambda_F \ll \lambda_{th}$ .
Clarification of Perturbation Nature:
- Isochoric Modes: Confirmed that purely isochoric perturbations (zero density/velocity perturbation) exist in the limits $\lambda \to 0$ and $\lambda \to \infty$ .
- Isobaric Modes: Debunked the existence of purely isobaric perturbations in a conducting plasma. While modes can be almost isobaric at intermediate wavelengths, true isobaric perturbations only exist if thermal conduction is zero.
Re-evaluation of "Catastrophic Cooling": The authors show that the "catastrophic cooling mode" proposed by Waters & Stricklan (2025) is mathematically identical to the standard isochoric thermal mode in the $\lambda \to \infty$ limit, rather than a distinct physical mechanism.

4. Key Results

Role of Parameter $R$ : The behavior of the system is primarily dictated by $R$ $R$ .
- $R < 0$ : The system is unstable in either isochoric or isobaric conditions (but not both). One thermal mode is unstable; acoustic modes are stable.
- $0 < R < 1/9$ : The system is unstable in all three criteria (isochoric, isobaric, isentropic). Two critical wavelengths ( $\lambda_-, \lambda_+$ ) exist where acoustic modes convert to thermal modes and back.
- $1/9 < R < 1$ : No critical wavelengths exist for mode conversion; the system transitions directly from stable to unstable acoustic/thermal modes as wavelength increases.
- $R > 1$ : The system is isentropically stable. Acoustic modes are stable; thermal instability depends on specific criteria.
Thermal Conduction Effects:
- For $\lambda_F \lesssim \lambda_{th}$ , the classification in Table 1 holds well.
- For $\lambda_F \gtrsim \lambda_{th}$ , the simple classification based on $R$ breaks down, and the number of critical wavenumbers can increase (up to three), making the instability domains complex and non-uniform.
Short vs. Long Wavelengths:
- $\lambda \to 0$ : All modes are stabilized by thermal conduction. The thermal mode becomes isochoric (damped temperature/pressure perturbations, no density/velocity). Acoustic modes become damped isothermal waves.
- $\lambda \to \infty$ : Thermal conduction is negligible. The system supports three thermal modes if $R < 0$ or one thermal + two acoustic modes if $R > 0$ .
Solar Coronal Application:
- Using CHIANTI atomic data for coronal loops, the authors find that for typical active region temperatures (1–3 MK), $R > 1$ .
- This implies that coronal rain formation is not driven by the "catastrophic cooling" mode (which peaks at $\lambda \to \infty$ ). Instead, the thermal mode achieves its maximum growth rate at a finite wavelength, suggesting a different mechanism for condensation formation than previously hypothesized for these specific conditions.

5. Significance

Theoretical Rigor: The paper provides the most complete analytical description of thermal instability in a uniform medium to date, correcting errors in recent literature and clarifying the mathematical limits of perturbation types.
Astrophysical Modeling: The updated classification (Table 1) and growth rate diagrams (Figures 4–10) serve as essential tools for interpreting numerical simulations of astrophysical plasmas, from the interstellar medium to galaxy clusters.
Solar Physics Impact: By demonstrating that $R > 1$ for typical coronal conditions, the study challenges the applicability of the "catastrophic cooling" hypothesis to solar coronal rain, suggesting that finite-wavelength thermal instability is the dominant driver.
Future Directions: The authors outline a path for future research into non-linear saturation regimes and the classification of instability in magnetized plasmas, which are critical for understanding laboratory and space plasmas.

In summary, this work re-establishes the foundational understanding of thermal instability, corrects recent simplifications, and provides a robust framework for analyzing plasma stability across diverse astrophysical environments.