An official website of the United States government
ABSTRACT The coupling state between ions and neutrals in the interstellar medium plays a key role in the dynamics of magnetohydrodynamic (MHD) turbulence, but is challenging to study numerically. In this work, we investigate the damping of MHD turbulence in a partially ionized medium using 3D two-fluid (ions + neutrals) simulations generated with the athenak code. Specifically, we examine the velocity, density, and magnetic field statistics of the two-fluid MHD turbulence in different regimes of neutral-ion coupling. Our results demonstrate that when ions and neutrals are strongly coupled, the velocity statistics resemble those of single-fluid MHD turbulence. Both the velocity structures and kinetic energy spectra of ions and neutrals are similar, while their density structures can be significantly different. With an excess of small-scale sharp density fluctuations in ions, the density spectrum in ions is shallower than that of neutrals. When ions and neutrals are weakly coupled, the turbulence in ions is more severely damped due to the ion-neutral collisional friction than that in neutrals, resulting in a steep kinetic energy spectrum and density spectrum in ions compared to the Kolmogorov spectrum. We also find that the magnetic energy spectrum basically follows the shape of the kinetic energy spectrum of ions, irrespective of the coupling regime. In addition, we find large density fluctuations in ions and neutrals and thus spatially inhomogeneous ionization fractions. As a result, the neutral-ion decoupling and damping of MHD turbulence take place over a range of length-scales.
more »
« less
Multi-fluid and kinetic models of partially ionized magnetic reconnection
Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomenon in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet–Parker scaling, where neutrals are unimportant, or a fully coupled Sweet–Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfvén speed, vA⋆. These differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain νniL/vA⋆≳1, where νni is the neutral–ion collision frequency and L is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow νniL/vA⋆<1, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.
more »
« less
- Award ID(s):
- 1931388
- PAR ID:
- 10662701
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- Physics of Plasmas
- Volume:
- 28
- Issue:
- 4
- ISSN:
- 1070-664X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Using 2.5 dimensional kinetic particle-in-cell (PIC) simulations, we simulate reconnection conditions appropriate for the magnetosheath and solar wind, i.e., plasma beta (ratio of gas pressure to magnetic pressure) greater than 1 and low magnetic shear (strong guide field). Changing the simulation domain size, we find that the ion response varies greatly. For reconnecting regions with scales comparable to the ion inertial length, the ions do not respond to the reconnection dynamics leading to “electron-only” reconnection with very large quasi-steady reconnection rates. Note that in these simulations the ion Larmor radius is comparable to the ion inertial length. The transition to more traditional “ion-coupled” reconnection is gradual as the reconnection domain size increases, with the ions becoming frozen-in in the exhaust when the magnetic island width in the normal direction reaches many ion inertial lengths. During this transition, the quasi-steady reconnection rate decreases until the ions are fully coupled, ultimately reaching an asymptotic value. The scaling of the ion outflow velocity with exhaust width during this electron-only to ion-coupled transition is found to be consistent with a theoretical model of a newly reconnected field line. In order to have a fully frozen-in ion exhaust with ion flows comparable to the reconnection Alfven speed, an exhaust width of at least several ion inertial lengths is needed. In turbulent systems with reconnection occurring between magnetic bubbles associated with fluctuations, using geometric arguments we estimate that fully ion-coupled reconnection requires magnetic bubble length scales of at least several tens of ion inertial lengths.more » « less
-
Magnetic reconnection has been intensively studied in fully ionized plasmas. However, plasmas are often partially ionized in astrophysical environments. The interactions between the neutral particles and ionized plasmas might strongly affect the reconnection mechanisms. We review magnetic reconnection in partially ionized plasmas in different environments from theoretical, numerical, observational and experimental points of view. We focus on mechanisms which make magnetic reconnection fast enough to compare with observations, especially on the reconnection events in the low solar atmosphere. The heating mechanisms and the related observational evidence of the reconnection process in the partially ionized low solar atmosphere are also discussed. We describe magnetic reconnection in weakly ionized astrophysical environments, including the interstellar medium and protostellar discs. We present recent achievements about fast reconnection in laboratory experiments for partially ionized plasmas.more » « less
-
Abstract Electron-only magnetic reconnection was first detected by the Magnetospheric Multiscale (MMS) mission in Earth’s turbulent magnetosheath. Its prevalence in kinetic-scale turbulence has attracted great interest in heliophysics, but also revealed a great challenge in identifying it in turbulence, where electron flows are often complex. The magnetic flux transport (MFT) method is an innovative method to identify active reconnection in numerical simulations and in situ observations of turbulent plasmas. Here we extend this method to distinguish between electron-only and ion-coupled reconnection. The coupling of magnetic field motion with plasma flows in the diffusion regions sets distinct scales in the MFT velocity. While both forms of reconnection satisfy the MFT signature for active reconnection as MFT inflows and outflows at an X-line, the specific electron-only MFT signature is only an electron-scale MFT outflow along the current sheet normal direction, whereas the specific ion-coupled signature is a two-scale, outer-ion-and-inner-electron-scale MFT outflow in the electron diffusion region, which evolves into a single ion-scale in the ion diffusion region. These signatures are verified in a simulation of gyrokinetic turbulence. The dependence of the MFT outflow on the distance downstream from the X-lines also agrees well with the framework of magnetic field–plasma flow coupling. The new MFT signatures provide a clear and reliable tool for investigating electron-only reconnection in turbulence, independent of the development of electron outflows. They are directly applicable to kinetic and fluid simulations, and have potential application to observations of diffusion region crossings by spacecraft missions such as MMS.more » « less
-
Abstract Nonequilibrium ionization (NEI) is essentially required for astrophysical plasma diagnostics once the plasma status departs from the assumption of ionization equilibrium. In this work, we perform fast NEI calculations combined with magnetohydrodynamic (MHD) simulations and analyze the ionization properties of a Petschek-type magnetic reconnection current sheet during solar eruptions. Our simulation reveals Petschek-type slow-mode shocks in the classical Spitzer thermal conduction models and conduction flux-limitation situations. The results show that under-ionized features can be commonly found in shocked reconnection outflows and thermal halo regions outside the shocks. The departure from equilibrium ionization strongly depends on plasma density. In addition, this departure is sensitive to the observable target temperature: the high-temperature iron ions are strongly affected by the effects of NEI. The under-ionization also affects the synthetic SDO/AIA intensities, which indicates that the reconstructed hot reconnection current sheet structure may be significantly underestimated either for temperature or apparent width. We also perform an MHD-NEI analysis on the reconnection current sheet in the classical solar flare geometry. Finally, we show the potential reversal between the under-ionized and over-ionized states at the lower tip of reconnection current sheets where the downward outflow collides with closed magnetic loops, which can strongly affect multiple SDO/AIA band ratios along the reconnection current sheet.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1063/5.0039860
