pith. sign in

arxiv: 2604.25027 · v1 · submitted 2026-04-27 · 🌌 astro-ph.SR

The Unstable Chromosphere: Effects of the Thermal Farley-Buneman Instability Across a Broad Range of Solar Chromospheric Conditions

Samuel Evans , Meers Oppenheim , Juan Mart\'inez-Sykora , Alexander Green This is my paper

Pith reviewed 2026-05-07 17:56 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords fieldtfbidrivingneutralinstabilityrelativesolaracross
0
0 comments X

The pith

Multifluid plasma simulations find that the Thermal Farley-Buneman Instability boosts neutral heating by (43±7)% on average, increases charged-species temperatures proportionally to driving field strength, and rotates current density components in the solar chromosphere.

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

The solar chromosphere is a thin, relatively cool layer in the sun's atmosphere where standard models often fail to match real observations. This paper uses computer simulations of a plasma process called the Thermal Farley-Buneman Instability. In this instability, flows of neutral gas create electric fields that push charged particles into small-scale turbulent motions and extra heating. The team ran 33 simulations covering a wide set of temperatures, densities, and magnetic field strengths typical of the chromosphere. They measured how the instability changes temperatures, gas motions, electric currents, and heating rates. The results show that the extra heating of neutral gas averages 43 percent with a standard deviation of 7 percent. Most simulations show at least a 29 percent increase. The instability also makes the average current flow more aligned with the driving electric field while reducing the perpendicular Hall component. Temperature rises in the charged particles scale directly with how far the driving field exceeds the theoretical minimum needed for the instability to grow.

Core claim

TFBI-driven turbulence increases neutral heating rates by (43±7)% on average, with 85% (28/33) of simulations having mean increase more than 29%, and a flat line of best fit suggesting zero correlation with driving field. The TFBI also causes average current density to rotate toward the driving field, with the Pederson/ambipolar component increasing by up to roughly 60% while the Hall component decreases by up to roughly 80%.

Load-bearing premise

The multifluid plasma simulation framework and the selected range of chromospheric parameters accurately capture the dominant physics without significant interference from other instabilities or unmodeled effects; the theoretical TFBI threshold field is taken as given from prior analytic work.

Figures

Figures reproduced from arXiv: 2604.25027 by Alexander Green, Juan Mart\'inez-Sykora, Meers Oppenheim, Samuel Evans.

Figure 1
Figure 1. Figure 1: Density perturbations, temperatures, and velocities from a typical TFBI suite simulation. Panels A (top four plots) show electron density perturbations relative to the mean, across the 2D simulation box at four snapshots in time. The remaining panels show values versus time for electrons (blue), HII (orange), (C+N+O+Ne)II (green), (Na+Mg+Al+Si+S)II (red), and (K+Ca+Cr+Fe+Ni)II (purple). Panel B shows the l… view at source ↗
Figure 2
Figure 2. Figure 2: Neutral heating rate (∂Tn/∂t) versus time, for the typical TFBI suite simulation shown in view at source ↗
Figure 3
Figure 3. Figure 3: Relative turbulent heating of each charged species (∆T (turb) s /T(0) s ) versus relative surplus driving field ((E(0)/Ethr) − 1), across the TFBI simulation suite. Each point represents the value computed from a single simulation, with error bars indicating one standard deviation, analogously to the computations from Section 3.1. Text in the top left corner of each panel indicates the plotted species, as … view at source ↗
Figure 4
Figure 4. Figure 4: Turbulent motions of each charged species (stddev(ns|⃗us|)/mean(ns)) versus relative surplus driving field ((E(0)/Ethr) − 1), across the TFBI simulation suite. Dot￾dashed lines near the bottom of each panel indicate “background turbulent motions” (due to particle noise) for each species. Electrons are excluded here due to significant background turbulent motions. Formatting is otherwise similar to view at source ↗
Figure 5
Figure 5. Figure 5: Turbulence-driven relative increase to inferred neutral heating rate from equa￾tion (2) versus relative surplus driving field ((E(0)/Ethr) − 1), across the TFBI simulation suite. Similarly to view at source ↗
Figure 6
Figure 6. Figure 6: Electron density perturbations relative to the mean, and all charged species’ average temperatures throughout simulation (3,2,3,2,4)F. Temperature plot has the same formatting as panel (C) of view at source ↗
Figure 7
Figure 7. Figure 7: Electron density perturbations relative to the mean, and all charged species’ average temperatures throughout simulation (4,3,5,2,4)G, with the same plot formatting as in view at source ↗
Figure 8
Figure 8. Figure 8: Density perturbations relative to the mean, at selected snapshots of simulation (3,1,4,4,6)E. From subplot to subplot, time increases from left to right (as indicated by labels on top), while each row shows a different charged species (as indicated by labels on the right) view at source ↗
Figure 9
Figure 9. Figure 9: Electron density perturbations relative to the mean (top), and all charged species’ average temperatures (bottom) throughout the (2,1,4,3,6) hybrid simulation. Tem￾perature plot has the same formatting as Panel (C) of view at source ↗
read the original abstract

In the coldest regions of the solar atmosphere, lingering discrepancies between models and observations may be caused by the Thermal Farley-Buneman Instability (TFBI). This meter-scale, electrostatic, collisional, multifluid plasma instability converts energy from neutral flows into turbulent motions and heating. In the neutral frame of reference, these neutral flows manifest as an electric field which can drive the TFBI. In this work, we simulate the TFBI across a broad range of solar chromospheric conditions. We find clear proportionality between between TFBI-driven relative temperature increases of charged species ($\Delta T_{s}^{\text{(turb)}} / T_{s}^{(0)}$) and driving electric field strength relative to the theoretical threshold field required for TFBI growth. We also discover a correlation between relative driving field strength and turbulent motions. Additionally, the TFBI consistently causes average current density to rotate towards the driving field, with Pederson/ambipolar component ($\vec{J} \cdot \hat{E}^{(0)}$) increasing by up to roughly 60% while the Hall component ($\vec{J} \cdot \hat{E}^{(0)} \times \hat{B}$) decreases in magnitude by up to roughly 80%. Meanwhile, we find TFBI-driven turbulence increases neutral heating rates by $(43\pm7)\%$ on average, with 85% (28/33) of simulations having mean increase more than 29%, and a flat line of best fit suggesting zero correlation with driving field.

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.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Central claims rest on the validity of the multifluid collisional plasma model for chromospheric conditions and the accuracy of the analytic TFBI growth threshold used to normalize driving fields; no new entities are postulated.

axioms (1)
  • domain assumption The theoretical threshold electric field for onset of the Thermal Farley-Buneman Instability is known accurately from prior linear theory and applies directly to the simulated nonlinear regime.
    Used to define relative driving field strength for all reported proportionality and correlation results.

pith-pipeline@v0.9.0 · 5587 in / 1469 out tokens · 100076 ms · 2026-05-07T17:56:35.440363+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

4 extracted references · 4 canonical work pages

  1. [1]

    Birdsall, C. K. 1991, IEEE Transactions on Plasma Science, 19, 65, doi: 10.1109/27.106800 Brooks, D. H., Ugarte-Urra, I., & Warren, H. P. 2015, Nature Communications, 6, 5947, doi: 10.1038/ncomms6947 Butcher, J. C. 2016, Numerical Methods for Ordinary Differential Equations (Newark, UNITED KINGDOM: John Wiley & Sons, Incorporated). http://ebookcentral.pro...

    work page doi:10.1109/27.106800 1991

  2. [2]

    A publicly available simulation of an enhanced network region of the Sun

    1146/annurev-astro-081817-052044 Carlsson, M., Hansteen, V. H., Gudiksen, B. V., Leenaarts, J., & De Pontieu, B. 2016, A&A, 585, A4, doi: 10.1051/0004-6361/201527226 Centeno, R., Blanco Rodr´ ıguez, J., Del Toro Iniesta, J. C., et al. 2017, ApJS, 229, 3, doi: 10.3847/1538-4365/229/1/3 Chintzoglou, G., De Pontieu, B., Mart´ ınez-Sykora, J., et al. 2021, Ap...

    work page doi:10.1051/0004-6361/201527226 2016

  3. [3]

    but here including a thermal equation for electrons rather than using an isothermal or adiabatic equation of state. 32 The electron thermal equation in hybrid EPPIC is given by: ∂Te ∂t =−⃗ ue · ∇Te − 2 3 Te∇ ·⃗ ue + 2 3 µe,nνe,n|⃗ ue|2 −δ e,nνe,n(Te −T n) (D1) whereµ e,n is the electron reduced mass,ν e,n is the electron neutral collision fre- quency, and...

    work page 2004

  4. [4]

    The electron temperature is then used in the usual field solve as described in (Young et al

    to update the electron temperature using equation (D1). The electron temperature is then used in the usual field solve as described in (Young et al. 2017). Figure 9.Electron density perturbations relative to the mean (top), and all charged species’ average temperatures (bottom) throughout the (2,1,4,3,6) hybrid simulation. Tem- perature plot has the same ...

    work page 2017