The role of ambipolar heating in the energy balance of solar prominences
Pith reviewed 2026-05-15 10:27 UTC · model grok-4.3
The pith
Ambipolar diffusion supplies part of the heating that balances radiative losses in solar prominence threads.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In one-dimensional models of solar prominences based on the Kippenhahn-Schlüter configuration with magnetic shear, ambipolar diffusion deposits energy from the non-force-free magnetic field into the plasma. Balancing this heating against radiative losses and thermal conduction produces stationary solutions with a cold, dense, partially ionized thread, a very thin prominence-corona transition region, and an extended hot fully ionized corona. Thread lengths consistent with observations are recovered for realistic values of central temperature, central pressure, magnetic field strength, and shear angle, while ambipolar diffusion simultaneously induces stationary flows associated with the down-d
What carries the argument
Ambipolar heating term in the energy equation, produced by the relative drift between neutrals and ions in the partially ionized plasma under the Lorentz force.
If this is right
- Ambipolar heating partly offsets radiative losses inside the cold thread.
- Gravitational drainage of neutrals produces stationary flows along the field in the partially ionized region.
- Thread length increases with stronger shear and weaker central pressure for fixed field strength.
- Realistic parameter choices yield thread lengths matching observed values.
- The effect must be retained when building multi-dimensional models.
Where Pith is reading between the lines
- Similar ambipolar heating could operate in other partially ionized structures such as spicules or chromospheric fibrils.
- The induced flows may influence the stability of the thread against Rayleigh-Taylor or other instabilities.
- Models that omit ambipolar diffusion may systematically underestimate the energy input needed to maintain observed prominence masses.
Load-bearing premise
The one-dimensional geometry with prescribed central parameters and shear captures the essential energy balance without additional heating mechanisms or multi-dimensional geometry.
What would settle it
High-resolution observations that find thread lengths independent of shear angle or that show no measurable stationary flows in the partially ionized core would falsify the predicted role of ambipolar diffusion.
read the original abstract
Solar prominence threads are typically located around magnetic dips, where cold and dense plasma is suspended against gravity in the hot corona thanks to the upward magnetic force. Because prominences are partially ionized, ambipolar diffusion can deposit part of the energy of their non-force-free magnetic field into the plasma. This ambipolar heating may therefore play a role in the energy balance of prominences. In this proof-of-concept work, we explore the effect of ambipolar diffusion in one-dimensional models that satisfy both mechanical equilibrium and energy balance. The magnetic configuration is based on the classic Kippenhahn-Schl\"uter model, incorporating a sheared magnetic field. The temperature profile along the magnetic field is computed numerically by balancing radiative losses, thermal conduction, and ambipolar heating. The resulting models consistently consist of a cold, dense, partially ionized thread with prominence core conditions, a very thin prominence-corona transition region, and an extended, hot, fully ionized region with coronal conditions. In addition to providing heating that partly compensates for radiative losses, ambipolar diffusion also gives rise to stationary flows associated with the gravitational drainage of neutrals in the partially ionized region. We investigate how the length of the cold threads depends on the central temperature, central pressure, magnetic field strength, and shear angle, and show that thread lengths compatible with observations are obtained for realistic values of these parameters. Therefore, we demonstrate that ambipolar diffusion plays a relevant role in this simple configuration, indicating that this effect should be incorporated into more elaborate multi-dimensional models and simulations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a proof-of-concept study using one-dimensional Kippenhahn-Schlüter models with a sheared magnetic field to investigate the role of ambipolar diffusion in the energy balance of solar prominences. By numerically solving for the temperature profile balancing radiative losses, thermal conduction, and ambipolar heating, the authors find that ambipolar diffusion contributes significantly to heating, enabling cold dense threads with observationally consistent lengths for realistic parameters, and inducing stationary flows from neutral drainage.
Significance. If the results hold, this work provides evidence that ambipolar heating is a relevant mechanism in prominence energy balance even in simplified geometries, offering a way to compensate radiative losses without additional ad-hoc heating and highlighting the need to include this effect in multi-dimensional models of the solar corona and prominences. The parameter study adds value by showing robustness for realistic inputs.
major comments (2)
- The description of the numerical solution of the energy-balance equation provides no convergence tests, grid-resolution sensitivity analysis, or error estimates on the computed temperature profiles and resulting thread lengths. This limits the strength of the quantitative claims that thread lengths match observations for realistic parameters.
- The reported thread lengths are obtained within the fixed 1D KS configuration with prescribed shear angle; no test or estimate is given of how the ambipolar heating term (proportional to the perpendicular current density) would change if the magnetic geometry were allowed to evolve self-consistently in 2D or 3D, where field-line curvature and cross-field currents could alter the local heating rate.
minor comments (2)
- The abstract states that ambipolar diffusion 'gives rise to stationary flows associated with the gravitational drainage of neutrals'; a short sentence explaining the origin of these flows from the momentum equation would improve accessibility.
- All symbols appearing in the energy-balance equation (radiative loss function, thermal conductivity, ambipolar diffusivity or heating rate) should be defined explicitly on first use, even if standard in the field.
Simulated Author's Rebuttal
We thank the referee for their constructive review and positive assessment of our proof-of-concept study. We have addressed the concerns about numerical validation by adding convergence tests and error estimates in the revised manuscript. We also clarify the scope and limitations of the 1D approach while maintaining that the results demonstrate the relevance of ambipolar heating even in this simplified geometry.
read point-by-point responses
-
Referee: The description of the numerical solution of the energy-balance equation provides no convergence tests, grid-resolution sensitivity analysis, or error estimates on the computed temperature profiles and resulting thread lengths. This limits the strength of the quantitative claims that thread lengths match observations for realistic parameters.
Authors: We agree that explicit convergence tests and error estimates would strengthen the quantitative aspects of our claims. In the revised manuscript, we have expanded the numerical methods section to include a description of the finite-difference scheme used to solve the energy-balance equation. We performed grid-resolution tests with 200, 500, and 1000 points along the field line, demonstrating that the cold-thread length converges to within 3% for resolutions of 500 points and finer. Residual errors in the energy balance are now quantified and remain below 1% throughout the domain for the parameter ranges explored. These additions support the robustness of the reported thread lengths for realistic inputs. revision: yes
-
Referee: The reported thread lengths are obtained within the fixed 1D KS configuration with prescribed shear angle; no test or estimate is given of how the ambipolar heating term (proportional to the perpendicular current density) would change if the magnetic geometry were allowed to evolve self-consistently in 2D or 3D, where field-line curvature and cross-field currents could alter the local heating rate.
Authors: We acknowledge that our models employ a fixed 1D Kippenhahn-Schlüter geometry with prescribed shear, so the ambipolar heating rate is computed from the static perpendicular current. Self-consistent 2D or 3D evolution would require time-dependent MHD simulations that allow the field to adjust under the influence of ambipolar diffusion and flows, which lies beyond the scope of this proof-of-concept paper. We have added a dedicated paragraph in the discussion section that explicitly states this limitation and notes that the heating term could be modified by evolving curvature and currents in higher dimensions. Our results nevertheless show that ambipolar heating is already sufficient to produce observationally consistent thread lengths in the simplest geometry that satisfies both force and energy balance, thereby motivating its inclusion in future multidimensional models. revision: partial
Circularity Check
No significant circularity in the derivation chain
full rationale
The paper solves the 1D mechanical equilibrium and energy balance equations numerically for a Kippenhahn-Schlüter configuration with prescribed shear, treating ambipolar heating as an explicit term in the energy equation alongside radiation and conduction. Thread lengths emerge directly from integrating the stationary balance for varied input parameters (central temperature, pressure, B strength, shear angle); no subset of data is fitted and then re-predicted, no self-citation supplies a load-bearing uniqueness theorem or ansatz, and no renaming of known results occurs. The central demonstration that ambipolar heating contributes to observed thread lengths is therefore a forward computation inside the stated 1D ansatz rather than a reduction to its own inputs.
Axiom & Free-Parameter Ledger
free parameters (4)
- central temperature
- central pressure
- magnetic field strength
- shear angle
axioms (2)
- domain assumption Mechanical equilibrium follows the classic Kippenhahn-Schlüter dipped-field configuration
- domain assumption Energy balance is achieved by equating radiative losses, thermal conduction, and ambipolar heating
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.