skip to main content


Title: Theory of the magnetothermal instability in coronal plasma flows

The theory of the magnetothermal instability (MTI) [D. A. Tidman and R. A. Shanny, Phys. Fluids 17, 1207 (1974)] is revisited through the lens of the stability of uniform systems. The linear stability analysis includes flow advection and Nernst transport. The instability criteria derived distinguish between the convective and the absolute nature of the perturbation growth. It is proven that, in the region where the Nernst and plasma blowoff velocities cancel, the MTI can be absolute and wave-packet perturbations grow in situ. This instability is mediated by the internal feedback between the Biermann battery and Righi–Leduc terms. The analysis is extended to derive the dispersion relation for short-wavelength perturbations developing in nonuniform profiles with the application to coronal plasmas. It is found that the condition for MTI requires the net B-field convection velocity to be small at the isothermal sonic section, and the plasma conditions in this section govern the dynamics of the instability. Analysis of hydro-equivalent implosions suggests that unstable perturbations undergo more e-foldings of growth in larger-size targets.

 
more » « less
Award ID(s):
2123496
NSF-PAR ID:
10371889
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Institute of Physics
Date Published:
Journal Name:
Physics of Plasmas
Volume:
29
Issue:
9
ISSN:
1070-664X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This paper presents a multi-instrument observational analysis of the equatorial plasma bubbles (EPBs) variation over the American sector during a geomagnetically quiet time period of 07–10 December 2019. The day-to-day variability of EPBs and their underlying drivers are investigated through coordinately utilizing the Global-scale Observations of Limb and Disk (GOLD) ultraviolet images, the Ionospheric Connection Explorer (ICON) in-situ and remote sensing data, the global navigation satellite system (GNSS) total electron content (TEC) observations, as well as ionosonde measurements. The main results are as follows: 1) The postsunset EPBs’ intensity exhibited a large day-to-day variation in the same UT intervals, which was fairly noticeable in the evening of December 07, yet considerably suppressed on December 08 and 09, and then dramatically revived and enhanced on December 10. 2) The postsunset linear Rayleigh-Taylor instability growth rate exhibited a different variation pattern. It had a relatively modest peak value on December 07 and 08, yet a larger peak value on December 09 and 10. There was a 2-h time lag of the growth rate peak time in the evening of December 09 from other nights. This analysis did not show an exact one-to-one relationship between the peak growth rate and the observed EPBs intensity. 3) The EPBs’ day-to-day variation has a better agreement with that of traveling ionospheric disturbances and atmospheric gravity waves signatures, which exhibited relatively strong wavelike perturbations preceding/accompanying the observed EPBs on December 07 and 10 yet relatively weak fluctuations on December 08 and 09. These coordinate observations indicate that the initial wavelike seeding perturbations associated with AGWs, together with the catalyzing factor of the instability growth rate, collectively played important roles to modulate the day-to-day variation of EPBs. A strong seeding perturbation could effectively compensate for a moderate strength of Rayleigh-Taylor instability growth rate and therefore their combined effect could facilitate EPB development. Lacking proper seeding perturbations would make it a more inefficient process for the development of EPBs, especially with a delayed peak value of Rayleigh-Taylor instability growth rate. 
    more » « less
  2. ABSTRACT

    Active Galactic Nuclei (AGN) are believed to provide the energy that prevents runaway cooling of gas in the cores of galaxy clusters. However, how this energy is transported and thermalized throughout the Intracluster Medium (ICM) remains unclear. In recent work, we showed that streaming cosmic rays (CRs) destabilize sound waves in dilute ICM plasmas. Here, we show that CR streaming in the presence of gravity also destabilizes a pressure-balanced wave. We term this new instability the CR buoyancy instability (CRBI). In stark contrast to standard results without CRs, the pressure-balanced mode is highly compressible at short wavelengths due to CR streaming. Maximal growth rates are of order (pc/pg)β1/2ωff, where pc/pg is the ratio of CR pressure to thermal gas pressure, β is the ratio of thermal to magnetic pressure, and ωff is the free-fall frequency. The CRBI operates alongside buoyancy instabilities driven by background heat fluxes, i.e. the heat-flux-driven buoyancy instability (HBI) and the magneto-thermal instability (MTI). When the thermal mean free path lmfp is ≪ the gas scale height H, the HBI/MTI set the growth rate on large scales, while the CRBI sets the growth rate on small scales. Conversely, when lmfp ∼ H and (pc/pg)β1/2 ≳ 1, CRBI growth rates exceed HBI/MTI growth rates even on large scales. Our results suggest that CR-driven instabilities may be partially responsible for the sound waves/weak shocks and turbulence observed in galaxy clusters. CR-driven instabilities generated near radio bubbles may also play an important role redistributing AGN energy throughout clusters.

     
    more » « less
  3. ABSTRACT

    Relativistic radiation-mediated shocks are likely formed in prodigious cosmic explosions. The structure and emission of such shocks are regulated by copious production of electron–positron pairs inside the shock-transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the non-linear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced-transport equation for the particles via pitch-angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to non-adiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of non-thermal distributions, arising from the microturbulence, can affect the observed shock-breakout signal in ways unaccounted for by current single-fluid models.

     
    more » « less
  4. A newly developed three-dimensional electrostatic fluid model solving continuity and current closure equations aims to study phenomena that generate ionospheric turbulence. The model is spatially discretized using a pseudo-spectral method with full Fourier basis functions and evolved in time using a four-stage, fourth-order Runge Kutta method. The 3D numerical model is used here to investigate the behavior and evolution of ionospheric plasma clouds. This problem has historically been used to study the processes governing the evolution of the irregularities in the F region of the ionosphere. It has been shown that these artificial clouds can become unstable and structure rapidly (i.e., cascade to smaller scales transverse to the ambient magnetic field). The primary mechanism which causes this structuring of ionospheric clouds is the E×B, or the gradient drift instability (GDI). The persistence and scale sizes of the resulting structures cannot be fully explained by a two-dimensional model. Therefore, we suggest here that the inclusion of three-dimensional effects is key to a successful interpretation of mid-latitude irregularities, as well as a prerequisite for a credible simulation of these processes. We investigate the results of 2D and 3D nonlinear simulations of the GDI and secondary Kelvin–Helmholtz instability (KHI) in plasma clouds for three different regimes: highly collisional (≈200 km), collisional (≈300 km), and inertial (≈450 km). The inclusion of inertial effects permits the growth of the secondary KHI. For the three different regimes, the overall evolution of structuring of plasma cloud occurs on longer timescales in 3D simulations. The inclusion of three-dimensional effects, in particular, the ambipolar potential in the current closure equation, introduces an azimuthal “twist“ about the axis of the cloud (i.e., the magnetic field B). This azimuthal “twist” is observed in the purely collisional regime, and it causes the perturbations to have a non-flute-like character (k‖≠0). However, for the 3D inertial simulations, the cloud rapidly diffuses to a state in which the sheared azimuthal flow is substantially reduced; subsequently, the cloud becomes unstable and structures, by retaining the flute-like character of the perturbations (k‖=0).

     
    more » « less
  5. Abstract

    Thermal effects in the ionospheric plasma instabilities including the Farley‐Buneman instability (FBI), the gradient‐drift instability (GDI), and the ion‐thermal instability (ITI) are analyzed, focusing on analytic analysis of the previously obtained general expression of the combined instability growth rate. It is shown that thermal effects lead to generally nonmonotonic behavior of the growth rate with the background electric fieldor, equivalently with the plasma drift speed, in contrast with the conventional quadratic (FBI) or linear (GDI) dependence. The threshold electric field for FBI is demonstrated to exist at all altitudes, although at higher altitudes and/or at longer wavelength the threshold fields become too high to be observed in the ionosphere. The GDI growth rate is shown to be significantly modified by the thermal effects at long wavelengths. In the absence of thermal effects, the growth rate is proportional to a cosine of the flow angle, the angle between the differential plasma drift and the wavevector. Thermal effects result in an additional phase shift that maximizes at an altitude around 120 km, using representative high‐latitude ionospheric parameters, and in an increase in the proportionality coefficient by up to a factor of 2. The study highlights the usefulness of an analytic treatment for revealing additional insights into the instability behavior in the broad range of ionospheric parameters that may remain hidden using only numerical treatment.

     
    more » « less