On the gravitational stratification of multi-fluid-multi-species plasma

This paper presents a numerical method for constructing multi-fluid-multi-species gravitational stratifications of solar atmospheric plasma that simultaneously satisfy hydrostatic and ionization equilibrium, thereby eliminating non-physical disturbances and numerical instabilities caused by initial in-equilibria in dynamic simulations.

The Gist

The Big Picture: The Solar Atmosphere is a Chaotic Crowd

Imagine the Sun's atmosphere not as a smooth, empty space, but as a massive, chaotic crowd of different types of people (particles) packed into a tall building.

• The Building: The Sun's atmosphere, stretching from the cool surface (photosphere) up to the super-hot outer layer (corona).
• The People: The "fluids" or particles. Some are heavy and slow (neutral atoms), some are light and fast (ions), and some are tiny and zippy (electrons).
• The Gravity: A constant force pulling everyone down toward the floor.
• The Heat: A force pushing everyone apart (pressure).

In the past, scientists tried to model this crowd by assuming everyone stood perfectly still in their own separate lines, balancing gravity and pressure independently. The authors of this paper say: "That doesn't work. People bump into each other, and that changes everything."

The Problem: The "Independent Lines" Mistake

For a long time, scientists modeled the Sun by treating different types of particles (like Hydrogen and Helium) as if they were in separate, non-communicating lines.

• The Old Way (Pure Hydrostatic Equilibrium): Imagine a line of heavy boxes and a line of light balloons. If you tell them to stand still against gravity, the heavy boxes will stack up tightly at the bottom, while the light balloons float way up high.
• The Reality: In the Sun, these particles are constantly bumping into each other (colliding). They aren't in separate lines; they are a mosh pit. The heavy boxes and light balloons are holding onto each other. If you treat them as separate, your model predicts that heavy elements (like Iron) will vanish from the lower atmosphere because they "should" sink too fast. But in reality, they stay mixed because the collisions hold them up.

If you start a computer simulation with the "separate lines" model, the computer gets confused. The particles realize they are in the wrong place, and the simulation crashes or produces nonsense because the "crowd" is trying to fix itself instantly.

The Solution: The "Team Huddle" Approach

The authors propose a new way to set up the starting position for their computer models. They call it Coupled Hydrostatic Equilibrium (cHE).

The Analogy: The Elevator Ride
Imagine this crowd of particles is inside an elevator that is slowly moving up a skyscraper.

1. The Goal: We want to know exactly where everyone is standing before the elevator starts shaking (before we add waves or explosions).
2. The Old Mistake: We told the heavy people to stand at the bottom and the light people to stand at the top, ignoring that they are holding hands.
3. The New Method: We assume that because they are bumping into each other so much, they act like one giant team. They share a "center of mass."
   • We calculate where the average person should be based on the total weight of the whole group.
   • Then, we let the collisions do the work. The heavy particles might try to sink, but the light particles push them back up. The light particles might try to float, but the heavy ones pull them down.
   • Result: Everyone finds a comfortable spot where the total group is balanced, even if individuals are drifting slightly up or down relative to each other.

How They Did It (The "Recipe")

The authors created a simple "recipe" (a numerical routine) to build this starting model:

1. Pick a Temperature: Decide how hot the "elevator" is at every floor.
2. Pick the Ionization: Decide how many people are "charged" (ions) vs. "neutral" (atoms) at each floor. This can be based on standard physics or complex simulations.
3. The Integration: Instead of solving a giant, impossible math puzzle all at once, they take tiny steps. They calculate the pressure at the bottom, move up a tiny bit, adjust for gravity and temperature, and repeat.
4. The Output: A perfectly balanced starting map where the "crowd" is stable, but individual particles are allowed to drift slightly (like people shifting their weight in a crowded room).

Why This Matters: The "Drift" is Real

One of the coolest findings is that in this new model, the particles do move, even when the Sun is "quiet."

• The Drift: Because the heavy particles want to sink and the light ones want to float, there is a constant, slow "drift" between them.
• The Analogy: Think of a crowd of people walking up a hill. The heavy backpacks (ions) want to slide down, but the people without backpacks (neutrals) are pushing them up. They end up walking at the same speed, but the backpacks are constantly slipping down the shoulders of the people carrying them.
• The Result: This "drift" is real physics. If you ignore it (like the old models did), you miss out on how the Sun actually behaves.

The Test: Waves and Wobbles

To prove their new model works, the authors ran simulations:

1. The Crash Test: They tried to start a simulation with the old "separate lines" model. It immediately became unstable and crashed because the particles were too far out of place.
2. The Smooth Ride: They used their new "team huddle" model. It was stable.
3. The Wave Test: They sent a wave (like a shout or a ripple) through the crowd.
   • In the old model, the heavy and light particles reacted differently and chaotically.
   • In the new model, the crowd reacted as a unit, but with the correct "drift" between the heavy and light particles. This allowed them to study how waves move through the Sun's atmosphere much more accurately.

The Takeaway

This paper gives scientists a better "starting line" for their solar simulations.

Instead of guessing where the particles should be and hoping they settle down (which often causes the simulation to break), they now have a method to set up the Sun's atmosphere exactly as it should be: a balanced, interacting crowd where heavy and light particles hold onto each other. This allows them to study solar flares, waves, and the heating of the Sun's corona with much greater accuracy, without the computer crashing due to "initial confusion."

In short: They stopped treating the Sun's atmosphere like separate lines of people and started treating it like a single, bumping, drifting crowd. And that makes all the difference.

Technical Summary

1. Problem Statement

The solar atmosphere is a complex, multi-layered environment ranging from the weakly ionized photosphere to the fully ionized corona. Modeling this region requires Multi-Fluid Multi-Species (MFMS) Magnetohydrodynamics (MHD) to account for the distinct behaviors of ions, neutrals, and electrons.

A critical challenge in setting up initial conditions for these models is establishing a gravitational stratification that satisfies two conditions simultaneously:

1. Hydrostatic Equilibrium (HE): The balance between pressure gradients and gravity.
2. Ionization Equilibrium (IE): The correct chemical composition (ionization fractions) dictated by local temperature and density.

The Core Conflict:

• Traditional Approach (Pure Hydrostatic Equilibrium - pHE): Assumes each fluid species is in independent hydrostatic equilibrium ($dP_i/dz = -\rho_i g$). Because different species have different masses and ionization states, they possess different pressure scale heights. This leads to ionization fractions that vary exponentially with height, often deviating by orders of magnitude from realistic atmospheric models (e.g., Model C7).
• Consequences: Using pHE initial conditions in simulations that include ionization/recombination physics leads to:
  • Non-physical initial imbalances.
  • Numerical instability.
  • Incorrect collisional frequencies (due to wrong densities).
  • Failure to model heavy elements correctly, as their scale heights differ significantly from hydrogen.

2. Methodology

The authors propose a method to construct a Coupled Hydrostatic Equilibrium (cHE) stratification.

Key Assumptions:

• In a static state (no high-frequency external driving), collisional interactions between fluids are sufficient to couple them, allowing the mixture to behave approximately as a single fluid regarding bulk motion.
• Ionization fractions can be pre-calculated based on either Statistical Equilibrium (SE) or Non-Equilibrium (NEQ) conditions derived from other models.

The Numerical Integration Routine:

1. Mean Atomic Mass: Instead of treating fluids independently, the authors define a mean atomic mass ($m_T$) for the bulk mixture:
   $$m_T = \frac{\sum N_i m_i}{\sum N_i}$$
   where $N_i$ and $m_i$ are the number density and mass of species $i$. Free electrons are included in the sum for mass calculation but assumed to have zero mass.
2. Total Pressure Equation: The hydrostatic balance is rewritten for the total plasma:
   $$\frac{dP}{dz} = -N_T m_T g = -\frac{P m_T g}{k_B T}$$
   where $P$ is the total thermal pressure and $N_T$ is the total number density.
3. Integration: A simple numerical integration routine (first-order left Riemann sum) is used to solve for pressure and density profiles step-by-step.
   • At each step, local ionization fractions (given as input) determine the partial number densities.
   • The total pressure is updated using the local mean mass and temperature.
   • This ensures the resulting stratification maintains the input ionization fractions while satisfying global hydrostatic balance.

Validation & Simulation:
The authors implemented this routine in the Ebysus code (a modular MFMS solver). They tested the resulting initial conditions by running 1D simulations:

• Static Tests: Running simulations without external driving forces to check for remnant disturbances and drift velocities.
• Dynamic Tests: Driving monochromatic Alfvén waves to observe the ponderomotive force and resulting plasma motions, comparing cHE vs. pHE initial conditions under various collisional regimes (weakly coupled, strongly coupled, and physical).

3. Key Contributions

• Novel Construction Method: Introduced a robust numerical routine to generate MFMS initial conditions that are simultaneously in hydrostatic and ionization equilibrium.
• Resolution of Scale Height Discrepancy: Demonstrated that treating fluids as a coupled bulk (cHE) prevents the unphysical exponential divergence of ionization fractions seen in pHE models.
• Handling of Heavy Elements: Showed that cHE allows for the inclusion of heavier elements (e.g., Iron, Neon) which would otherwise "disappear" or be incorrectly distributed in pHE models due to their different scale heights.
• Validation of Drift Velocities: Identified that in a cHE stratification, individual fluids possess non-zero drift velocities (to balance gravity and pressure gradients via collisions), while the total bulk velocity remains zero. This is a physically consistent state for a static, multi-fluid atmosphere.

4. Results

• Equilibrium Stability: Simulations initialized with cHE stratification remain stable over time. The density and ionization fractions change negligibly, confirming the initial state is close to true equilibrium.
• Failure of pHE: Simulations initialized with pHE stratification developed numerical instability when ionization/recombination was active, due to the massive initial imbalance in ionization fractions.
• Drift Velocities: In the cHE model, significant drift velocities appear between neutral and ionized fluids (especially Helium), particularly in the Transition Region where pressure gradients are steep. These drifts are physically real and necessary to maintain the coupled equilibrium.
• Ponderomotive Force Effects: When Alfvén waves were introduced:
  • cHE models correctly captured the decoupling of neutral helium from the ionized fluid, leading to distinct longitudinal motions driven by the ponderomotive force.
  • pHE models failed to reproduce these effects accurately and could not be reliably tested with ionization physics due to instability.
• NEQ Compatibility: The method successfully accommodated Non-Equilibrium (NEQ) ionization fractions extracted from 2.5D radiative MHD snapshots, proving its flexibility for studying dynamic atmospheric features like spicules.

5. Significance

This work provides a foundational tool for advancing solar atmospheric modeling:

• Improved Initial Conditions: It eliminates the "start-up shock" caused by initial non-equilibria, allowing researchers to study specific dynamics (waves, reconnection, chemical fractionation) without the results being dominated by numerical artifacts.
• Chemical Fractionation Studies: By correctly handling heavy elements and their scale heights, the cHE method enables accurate modeling of the First Ionization Potential (FIP) effect, a major unsolved problem in solar physics.
• Physical Realism: It acknowledges that in a static, collisional plasma, fluids are not independent; they drift relative to one another to maintain a coupled hydrostatic state. This bridges the gap between single-fluid approximations and complex multi-fluid dynamics.
• Future Applications: The method is applicable to 1D, 2D, and 3D models (with force-free magnetic fields) and can be integrated into codes like Ebysus and Bifrost to study the low solar atmosphere with higher fidelity.