An official website of the United States government
Abstract Models fail to reproduce observations of the coldest parts of the Sun’s atmosphere, where interactions between multiple ionized and neutral species prevent an accurate MHD representation. This paper argues that a meter-scale electrostatic plasma instability develops in these regions and causes heating. We refer to this instability as the Thermal Farley–Buneman Instability (TFBI). Using parameters from a 2.5D radiative MHD Bifrost simulation, we show that the TFBI develops in many of the colder regions in the chromosphere. This paper also presents the first multifluid simulation of the TFBI and validates this new result by demonstrating close agreement with theory during the linear regime. The simulation eventually develops turbulence, and we characterize the resulting wave-driven heating, plasma transport, and turbulent motions. These results all contend that the effects of the TFBI contribute to the discrepancies between solar observations and radiative MHD models.
more »
« less
This content will become publicly available on June 3, 2026
Multifluid and Kinetic 2D and 3D Simulations of Thermal Farley–Buneman Instability Turbulence in the Solar Chromosphere
Models currently fail to reproduce observations of the coldest regions in the Sun’s atmosphere, though recent work suggests the thermal Farley–Buneman instability (TFBI) may play a critical role. This meter-scale, electrostatic, multifluid plasma instability causes turbulence and heating in the coldest regions of the solar chromosphere. This paper describes how TFBI turbulence and heating varies across multifluid 2D, kinetic 2D, and kinetic 3D simulations. It also presents the first 3D simulations of the TFBI. We find that multifluid and kinetic 2D simulations produce similar results overall, despite using vastly different approaches. Additionally, our kinetic 3D simulations produce a similar or somewhat larger amount of heating compared to 2D, as contributions from the parallel electric field account for only (13 ± 2.5)% of the total turbulent heating in 3D.
more »
« less
- Award ID(s):
- 2532187
- PAR ID:
- 10598154
- Publisher / Repository:
- American Astronomical Society
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 986
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 23
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We present high-resolution 2D and 3D simulations of magnetized decaying turbulence in relativistic, resistive magnetohydrodynamics. The simulations show dynamic formation of large-scale intermittent long-lived current sheets being disrupted into plasmoid chains by the tearing instability. These current sheets are locations of enhanced magnetic-field dissipation and heating of the plasma. We find magnetic energy spectra ∝k−3/2, together with strongly pronounced dynamic alignment of Elsässer fields and of velocity and magnetic fields, for strong guide-field turbulence, whereas we retrieve spectra ∝k−5/3for the case of a weak guide-field.more » « less
-
Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair PlasmasAbstract The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.more » « less
-
Plasmoid instability accelerates reconnection in collisional plasmas by transforming a laminar reconnection layer into numerous plasmoids connected by secondary current sheets in two dimensions (2D) and by fostering self-generated turbulent reconnection in three dimensions (3D). In large-scale astrophysical and space systems, plasmoid instability likely initiates in the collisional regime but may transition into the collisionless regime as the fragmentation of the current sheet progresses toward kinetic scales. Hall magnetohydrodynamics (MHD) models are widely regarded as a simplified yet effective representation of the transition from collisional to collisionless reconnection. However, plasmoid instability in 2D Hall MHD simulations often leads to a single-X-line reconnection configuration, which significantly differs from fully kinetic particle-in-cell simulation results. This study shows that single-X-line reconnection is less likely to occur in 3D compared to 2D. Moreover, depending on the Lundquist number and the ratio between the system size and the kinetic scale, Hall MHD can also realize 3D self-generated turbulent reconnection. We analyze the features of the self-generated turbulent state, including the energy power spectra and the scale dependence of turbulent eddy anisotropy.more » « less
-
Abstract The proton–alpha drift instability is a possible mechanism of the alpha-particle deceleration and the resulting proton heating in the solar wind. We present hybrid numerical simulations of this instability with particle-in-cell ions and a quasi-neutralizing electron fluid for typical conditions at 1 au. For the parameters used in this paper, we find that fast magnetosonic unstable modes propagate only in the direction opposite to the alpha-particle drift and do not produce the perpendicular proton heating necessary to accelerate the solar wind. Alfvén modes propagate in both directions and heat the protons perpendicularly to the mean magnetic field. Despite being driven by the alpha temperature anisotropy, the Alfvén instability also extracts the energy from the bulk motion of the alpha particles. In the solar wind, the instabilities operate in a turbulent ambient medium. We show that the turbulence suppresses the Alfvén instability but the perpendicular proton heating persists. Unlike a static nonuniform background, the turbulence does not invert the sense of the proton heating associated with the fast magnetosonic instability and it remains preferentially parallel.more » « less
