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

F. Zhang; J. Mart\'inez-Sykora; Q. M. Wargnier; V. H. Hansteen

arxiv: 2601.05321 · v3 · submitted 2026-01-08 · 🌌 astro-ph.SR · physics.plasm-ph

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

F. Zhang , J. Mart\'inez-Sykora , Q. M. Wargnier , V. H. Hansteen This is my paper

Pith reviewed 2026-05-16 15:37 UTC · model grok-4.3

classification 🌌 astro-ph.SR physics.plasm-ph

keywords gravitational stratificationmulti-fluid plasmasolar atmosphereionization equilibriumhydrostatic equilibriumnumerical integrationtransition region

0 comments

The pith

Numerical routine constructs multi-fluid gravitational stratifications in hydrostatic equilibrium for any ionization fractions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a numerical integration method to build initial conditions for multi-fluid models of the solar atmosphere that are both gravitationally stratified and in ionization equilibrium. This prevents starting simulations with imbalances that produce artificial motions or instabilities. The approach takes ionization fractions either from statistical equilibrium at given temperatures or from other models, then integrates density profiles for each species and fluid so the total pressure gradient balances gravity. If the method works, models can examine real dynamics such as waves or reconnection without transients caused by setup errors.

Core claim

A gravitational stratification constructed using the present numerical integration routine can be in hydrostatic equilibrium with any given ionization fractions of multi-fluid-multi-species plasmas. The routine sets the density profiles for each fluid and species while incorporating the specified ionization states, ensuring the sum of partial pressure gradients balances gravity. Without external driving, fluid decoupling appears in the transition region, but the total velocity of all fluids remains zero.

What carries the argument

Numerical integration routine that computes density and pressure profiles for each fluid and species to satisfy hydrostatic equilibrium given fixed ionization fractions.

Load-bearing premise

Collisional interactions between fluids are assumed sufficient to couple all fluids when no high-frequency external driving force is imposed.

What would settle it

A simulation started from the constructed stratification that develops non-zero total velocity or pressure imbalances in the absence of any external forces would show the method fails to maintain equilibrium.

read the original abstract

Context. The solar atmosphere is gravitationally stratified and consists of several layers at temperatures different by orders of magnitude. Consequently, the solar atmospheric plasma changes from weakly ionized in the photosphere, partially ionized in the chromosphere, to eventually fully ionized in the corona. However, integrating ionization and recombination processes into multi-fluid solar plasma models with gravitational stratification remains nontrivial. Aims. We intend to provide a method for constructing multi-fluid-multi-species gravitational stratification that satisfies ionization equilibrium and hydrostatic equilibrium at the same time, avoiding causing non-physical disturbances and numerical instability due to initial in-equilibria. Methods. We assume that collisional interactions between fluids are sufficient for coupling all fluids when there is no high-frequency external driving force imposed. Ionization fractions can be (I) calculated assuming ionization in statistical equilibrium at any given temperature, or (II) extracted from other atmospheric models. A simple numerical integration routine is then designed and used to construct multi-fluid-multi-species gravitational stratifications. Results. A gravitational stratification constructed using the present numerical integration routine can be in hydrostatic equilibrium with any given ionization fractions of multi-species plasmas. Meanwhile, without any dynamic driving force, fluid decoupling appears, particularly in the transition region of the constructed stratification, while the total velocity of all fluids remains zero. Conclusions. A gravitational stratification constructed using the present routine can be used in multi-fluid-multi-species models to study specific dynamics without being affected by initial in-equilibria.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper offers a practical numerical routine for initializing multi-fluid solar atmosphere models with consistent stratification and ionization balance, but the reported species decoupling without driving creates a tension with the equilibrium claim.

read the letter

This paper's main point is a numerical routine that constructs gravitational stratifications for multi-fluid multi-species plasmas while satisfying both ionization equilibrium and hydrostatic balance for given ionization fractions. It handles a practical problem in solar atmosphere modeling well. Getting initial conditions right in these simulations is tricky because of the big temperature jumps and changing ionization states, and this method provides a way to avoid starting with imbalances that trigger instabilities. The routine is described clearly as a simple integration that takes ionization fractions from equilibrium calculations or other models and builds the profiles accordingly. They confirm that the summed velocity stays at zero, which keeps the setup static overall. The soft spot comes from the results on fluid behavior. The paper finds that without external driving, fluids decouple in the transition region with non-zero relative velocities, even though collisions are assumed to couple them sufficiently. Since friction terms in the multi-fluid equations depend on those relative velocities, non-zero values should generate forces that prevent a steady equilibrium. This suggests the hydrostatic balance holds mainly for the total fluid rather than each species individually, which weakens the central claim somewhat. It's not a showstopper but needs clearer discussion on what equilibrium means here and how it affects subsequent dynamics. This is for solar physicists running multi-fluid simulations who need reliable starting points. The work shows honest engagement with the setup challenges, so it merits a serious referee to verify the numerics and perhaps add tests on residual accelerations. I recommend sending it for peer review.

Referee Report

2 major / 2 minor

Summary. The paper presents a numerical integration routine for constructing gravitationally stratified multi-fluid multi-species plasma profiles in the solar atmosphere. Given arbitrary ionization fractions (computed from statistical equilibrium at a prescribed temperature or taken from other models), the method integrates the hydrostatic balance equations under the assumption that collisional coupling is sufficient to keep all fluids together in the absence of high-frequency external driving. The resulting profiles are claimed to satisfy both ionization equilibrium and hydrostatic equilibrium simultaneously and are proposed as initial conditions for multi-fluid simulations to avoid spurious initial transients.

Significance. If the central claim is valid, the routine supplies a practical, parameter-free initialization procedure for multi-fluid solar-atmosphere models that incorporate partial ionization. This could reduce numerical artifacts in studies of chromospheric dynamics, wave propagation, or heating. The approach is grounded in a forward integration that takes ionization fractions as direct inputs rather than fitting parameters, which is a methodological strength.

major comments (2)

[Results] Results (as summarized in the abstract): the reported fluid decoupling in the transition region, with non-zero relative velocities between species while the total velocity remains zero, directly contradicts the hydrostatic-equilibrium claim. In the multi-fluid momentum equations, collisional friction/drag terms are proportional to relative velocity; non-zero v_rel therefore generates net forces that must produce acceleration, violating the time-independent, force-balanced state required for hydrostatic equilibrium. This inconsistency is load-bearing for the assertion that the routine produces equilibrium stratifications for arbitrary ionization fractions.
[Methods] Methods (assumption paragraph): the premise that 'collisional interactions between fluids are sufficient for coupling all fluids when there is no high-frequency external driving force imposed' is not reconciled with the observed decoupling. The numerical integration therefore cannot guarantee equilibrium unless the friction terms are shown to remain identically zero, which the results indicate they do not.

minor comments (2)

[Results] The manuscript supplies no quantitative diagnostics (residual force imbalance, convergence metrics, or comparison against analytic single-fluid limits) to demonstrate that hydrostatic equilibrium is actually achieved to within a stated tolerance.
[Methods] Notation for the multi-fluid momentum and continuity equations should be introduced explicitly with equation numbers before the integration routine is described.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for providing detailed comments that highlight important aspects of our approach. We respond to each major comment below and indicate the revisions we will make.

read point-by-point responses

Referee: [Results] Results (as summarized in the abstract): the reported fluid decoupling in the transition region, with non-zero relative velocities between species while the total velocity remains zero, directly contradicts the hydrostatic-equilibrium claim. In the multi-fluid momentum equations, collisional friction/drag terms are proportional to relative velocity; non-zero v_rel therefore generates net forces that must produce acceleration, violating the time-independent, force-balanced state required for hydrostatic equilibrium. This inconsistency is load-bearing for the assertion that the routine produces equilibrium stratifications for arbitrary ionization fractions.

Authors: We appreciate the referee bringing this potential inconsistency to our attention. In the multi-fluid system, the collisional drag terms between species are equal and opposite, so they cancel exactly in the total momentum equation. Thus, the total plasma can satisfy hydrostatic equilibrium (total pressure gradient balancing gravity with zero total velocity) independently of the relative velocities between species. The non-zero relative velocities in the transition region result from the sharp gradients in ionization fractions and partial pressures. While these would cause the individual fluids to relax over time, the profile provides a consistent initial state for the total system with no net force imbalance. We will revise the abstract, results, and discussion sections to clarify that hydrostatic equilibrium applies to the composite fluid and to quantify the magnitude of the relative velocities and associated relaxation timescales. revision: partial
Referee: [Methods] Methods (assumption paragraph): the premise that 'collisional interactions between fluids are sufficient for coupling all fluids when there is no high-frequency external driving force imposed' is not reconciled with the observed decoupling. The numerical integration therefore cannot guarantee equilibrium unless the friction terms are shown to remain identically zero, which the results indicate they do not.

Authors: The assumption is intended to justify treating the fluids as coupled for the purpose of constructing the gravitational stratification via integration of the hydrostatic balance. The decoupling observed in the results is localized to the transition region where ionization changes rapidly. We will update the methods section to explicitly reconcile this by explaining that the coupling assumption holds for the total momentum balance, while allowing for small relative drifts due to species-specific forces. We will also add a supplementary analysis demonstrating that the friction terms are balanced by differences in species pressure gradients in the individual equations, maintaining the overall equilibrium. revision: yes

Circularity Check

0 steps flagged

No circularity detected in the numerical construction of multi-fluid stratification

full rationale

The paper presents a forward numerical integration routine that accepts ionization fractions (either from statistical equilibrium or external models) as explicit inputs and computes density profiles to enforce hydrostatic balance under the stated collisional-coupling assumption. No parameter is fitted to a data subset and then relabeled as a prediction; the output stratification is produced directly by integrating the governing equations rather than by any self-referential definition or self-citation chain. The reported appearance of relative velocities in the transition region is an output of the multi-fluid momentum equations and does not retroactively make the construction circular. The derivation therefore remains self-contained and does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The method rests on standard domain assumptions of hydrostatic equilibrium and ionization equilibrium plus the key coupling assumption stated in the abstract. No free parameters or new entities are introduced; ionization fractions are treated as external inputs.

axioms (2)

domain assumption Collisional interactions between fluids are sufficient for coupling all fluids when there is no high-frequency external driving force imposed.
Explicitly stated as the basis for the multi-fluid treatment in the methods description.
domain assumption Ionization fractions can be calculated assuming statistical equilibrium at any given temperature or extracted from other atmospheric models.
Used as direct input to the numerical integration routine.

pith-pipeline@v0.9.0 · 5572 in / 1298 out tokens · 67553 ms · 2026-05-16T15:37:29.718897+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We assume that collisional interactions between fluids are sufficient for coupling all fluids when there is no high-frequency external driving force imposed... dP/dz = −NT mT g (Eq. 8); numerical integration P_{j+1} = P_j − (p_gT,j mT,j g_j / kB T_j) Δz_j (Eq. 11)
IndisputableMonolith/Foundation/AbsoluteFloorClosure.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

A gravitational stratification constructed using the present numerical integration routine can be in hydrostatic equilibrium with any given ionization fractions of multi-species plasmas.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.