Low-beta Plasma Research Articles

Simulation of radio-frequency heating and fast-ion generation in Wendelstein 7-X

The next scientific operation phase of Wendelstein 7-X (W7-X) is scheduled to begin in late autumn of 2022 and will, for the first time, include experiments in which the ICRH (ion-cyclotron-resonance heating) antenna will be used.In addition to heating the plasma, this system will generate fast ions and thus offers a new way to assess fast-ion confinement in a stellarator such as W7-X. The first plasmas that will be used for the upcoming ICRH operation will be Helium-4 plasmas with a small Hydrogen minority on the order of about 10%. In tokamaks such plasmas typically offer good power absorption and are thus considered a safe way for gaining first experiences with the new antenna in W7-X.This assessment is confirmed by the SCENIC simulations carried out in this contribution that use profiles foreseen for the upcoming campaign as input. The simulations are carried out in the standard configuration of W7-X in low-beta (0.3% ≲ 〈β〉 ≲ 1%) plasmas. However, also scans over minority concentration and background-plasma density are performed. We find that the power absorbed by the Hydrogen minority directly from the radio-frequency wave is typically (provided that the minority concentration is not too high) on the order of about 90% with the rest going to the electrons. Very little power goes to the Helium-4 ions. Under the present simulation conditions only fast-ion energies up to about E ≈ 50 keV can be reached.Combining SCENIC and ASCOT simulations enables us to track lost particles through the scrape-off-layer to the 3D wall of W7-X and to compute wall loads caused by ICRH. The results show that the wall loads that can be expected from ICRH under the first operating conditions are benign.

What are the fundamental modes of energy transfer and partitioning in the coupled Magnetosphere-Ionosphere system?

The fundamental processes responsible for energy exchange between large-scale electromagnetic fields and plasma are well understood theoretically, but in practice these theories have not been tested. These processes are ubiquitous in all plasmas, especially at the interface between high and low beta plasmas in planetary magnetospheres and other magnetic environments. Although such boundaries pervade the plasma Universe, the processes responsible for the release of the stored magnetic and thermal plasma energy have not been fully identified and the importance of the relative impact of each process is unknown. Despite advances in understanding energy release through the conversion of magnetic to kinetic energy in magnetic reconnection, how the extreme pressures in the regions between stretched and more relaxed field lines in the transition region are balanced and released through adiabatic convection of plasma and fields is still a mystery. Recent theoretical advances and the predictions of large-scale instabilities must be tested. In essence, the processes responsible remain poorly understood and the problem unresolved. The aim of the White Paper submitted to ESA’s Voyage 2050 call, and the contents of this paper, is to highlight three outstanding open science questions that are of clear international interest: (i) the interplay of local and global plasma physics processes: (ii) the partitioning during energy conversion between electromagnetic and plasma energy: and (iii) what processes drive the coupling between low and high beta plasmas. We present a discussion of the new measurements and technological advances required from current state-of-the-art, and several candidate mission profiles with which these international high-priority science goals could be significantly advanced.

Parametric Decay of Alfvénic Wave Packets in Nonperiodic Low-beta Plasmas

Abstract The parametric decay of finite-size Alfvén waves in nonperiodic low-beta plasmas is investigated using one-dimensional (1D) hybrid simulations. Compared with the usual small periodic system, a wave packet in a large system under the absorption boundary condition shows different decay dynamics, including reduced energy transfer, localized density cavitation, and ion heating. The resulting Alfvén wave dynamics are influenced by several factors relating to this instability, including the growth rate, central wave frequency, and unstable bandwidth. A final steady state of the wave packet may be achieved when the instability does not have enough time to develop within the residual packet, and the packet size shows well-defined scaling dependencies on the growth rate, wave amplitude, and plasma beta. Under the proper conditions, enhanced secondary decay can also be excited in the form of a narrow, amplified wave packet. These results may help to interpret laboratory and spacecraft observations of Alfvén waves, and to refine our understanding of the associated energy transport and ion heating.

A new locally divergence-free WLS-ENO scheme based on the positivity-preserving finite volume method for ideal MHD equations

In this paper, the WLS-ENO (Weighted-Least-Squares based Essentially Non-Oscillatory) reconstruction is modified to maintain the conservation of the cell average values. Furthermore, the divergence-free constraint is combined with the conservative WLS-ENO reconstruction, which can make the magnetic field locally divergence-free. The main merit of the proposed reconstruction scheme is that it can keep both the divergence-free constraint and ENO property for the magnetic field without using any limiter. We apply the scheme to the simulations of ideal MHD equations within the framework of a positivity-preserving finite volume method. The convergence, the magnetic field divergence error and the capability for low plasma-beta of the scheme are tested by some MHD benchmark problems.

High-fidelity kinetic modeling of instabilities and gyromotion physics in nonuniform low-beta plasmas

A fourth-order accurate continuum kinetic Vlasov solver and a systematic method for constructing customizable kinetic equilibria are demonstrated to be powerful tools for the study of nonuniform collisionless low-beta plasmas. The noise-free methodology is applied to investigate two gradient-driven instabilities in 4D (x,y,vx,vy) phase space: the Kelvin–Helmholtz instability and the lower hybrid drift instability. Nonuniform two-species configurations where ion gyroradii are comparable to gradient scale lengths are explored. The approach sheds light on the evolution of the pressure tensor in Kelvin–Helmholtz instabilities and demonstrates that the associated stress tensor deviates significantly from the gyroviscous stress tensor. Even at high magnetization, first-order approximations to finite-gyromotion physics are shown to be inadequate for the Kelvin–Helmholtz instability, as shear scales evolve to become on par with gyromotion scales. The methodology facilitates exploring transport and energy partitioning properties associated with lower hybrid drift instabilities in low-beta plasma configurations. Distribution function features are captured in detail, including the formation of local extrema in the vicinity of particle-wave resonances. The approach enables detailed targeted investigations and advances kinetic simulation capability for plasmas in which gyromotion plays an important role.

Ion flow driven by low frequency Alfvén waves in a low-beta plasma

In a low-beta plasma, the ion flow in the parallel direction along the background magnetic field is investigated when ions are accelerated by low-frequency polarized Alfvén waves with the finite amplitude propagating along the magnetic field due to nonresonant interactions. The results indicate that the magnitude of the ion flow is closely related to not only the amplitude but also the frequency and the polarization of the wave, which is verified by a test particle simulation. The ion flows driven by the nondispersive and dispersive Alfvén waves are proportional to vAα21±ε2 and vAα21±ε3/2, respectively, where α is the ratio of the magnetic field component of the Alfvén wave to the background magnetic field, ε is the ratio of the wave frequency to the ion cyclotron frequency, and the positive sign corresponds to the right-handed Alfvén wave and the negative sign corresponds to the left-handed Alfvén wave, respectively. For a wave with finite frequency, the ion flow is different from different polarized waves in low beta plasmas. In particular, the saturation value of ion flow has a minimum threshold for the left-handed polarized Alfvén wave and a maximum threshold for the right-handed polarized Alfvén wave. If the frequency is less than 1/10 of the ion cyclotron frequency, the Alfvén wave can be seen as propagating in nondispersive medium. When the frequency of the Alfvén wave is far less than the frequency of the ion cyclotron, the ion flows driven by the left-handed and right-handed polarized Alfvén waves with and without wave dispersion tend to be the same.

Threaded-field-line Model for the Low Solar Corona Powered by the Alfvén Wave Turbulence

Abstract We present an updated global model of the solar corona, including the transition region. We simulate the realistic three-dimensional (3D) magnetic field using the data from the photospheric magnetic field measurements and assume the magnetohydrodynamic (MHD) Alfvén wave turbulence and its nonlinear dissipation to be the only source for heating the coronal plasma and driving the solar wind. In closed-field regions, the dissipation efficiency in a balanced turbulence is enhanced. In the coronal holes, we account for a reflection of the outward-propagating waves, which is accompanied by the generation of weaker counterpropagating waves. The nonlinear cascade rate degrades in strongly imbalanced turbulence, thus resulting in colder coronal holes. The distinctive feature of the presented model is the description of the low corona as almost-steady-state low-beta plasma motion and heat flux transfer along the magnetic field lines. We trace the magnetic field lines through each grid point of the lower boundary of the global corona model, chosen at some heliocentric distance, R = R b ∼ 1.1R ⊙, well above the transition region. One can readily solve the plasma parameters along the magnetic field line from 1D equations for the plasma motion and heat transport together with the Alfvén wave propagation, which adequately describe the physics within the heliocentric distance range R ⊙ < R < R b , in the low solar corona. By interfacing this threaded-field-line model with the full MHD global corona model at r = R b , we find the global solution and achieve a faster-than-real-time performance of the model on ∼200 cores.

Field Localization and Density Cavitation in Low-Beta Plasmas

In the present paper, filamentous structure formation, associated turbulent spectrum, and density cavity formation phenomena have been investigated for low-β plasma β≪me/mi applicable to the auroral region. A set of dimensionless equations governing the dynamics of three dimensionally propagating inertial Alfvén wave (3D-IAW) and perpendicularly propagating magnetosonic wave (PMSW) has been developed. Ponderomotive force due to 3D-IAW has been included in the dynamics of the PMSW. Numerical simulation has been performed to study the nonlinear coupling of these two waves. From the obtained results, we found that the field intensity localization takes place which may further lead to the additional dissipation/turbulence process for particle heating and acceleration in space plasma. The associated turbulent spectrum is obtained with scaling nearly k−4.28 at smaller scales (in the dissipation range). Relevance of the obtained results with the observations reported by various spacecrafts such as Hawkeye and Heos 2 has been discussed. Also, density fluctuations (depletion) of ∼0.10 n0 are calculated, which are consistent with the FAST spacecraft observation reported.

An Analytical Treatment for Particle Acceleration at Shocks inside Coronal Mass Ejections near 1 au

Abstract We present an analytical treatment for time-dependent diffusive shock acceleration at shocks inside magnetic clouds (MCs) observed near 1 au. The model includes the effects of (i) spatial diffusion of test particles upstream and downstream of the shock, (ii) proton advection with the plasma inside MCs, (iii) a reflecting boundary at distance L upstream of the shock to mimic the boundary of the MCs, and (iv) particle leakage out of the system at a constant rate, possibly through open field lines introduced by magnetic reconnection between the closed field lines of the MC and open field lines in the corona or heliosphere. The analysis reveals that the mean time for accelerating particles from p 0 to p is naturally reduced if the MC characteristic length is much smaller than the spatial diffusion length of energetic protons upstream of the shock. However, because most shocks inside MCs observed at 1 au are located in the back half of the MC, the time that the shock has propagated into the MCs is not sufficient to cause significant SEP enhancement—even with a reflecting boundary—if particles are only injected from the low-beta plasma inside MCs. To cause large SEP enhancements inside the shock–MC structure, magnetic reconnection at the back MC is essential to allow particles energized by the shock prior to its interaction with the MC to enter the MC. These particles consequently become the seed energetic protons that are reaccelerated at the shock inside MC.

Ion Jets Within Current Sheets in the Martian Magnetosphere

AbstractMagnetic reconnection is often suggested to be implicated in a variety of dynamic phenomena observed in the complex magnetosphere of Mars. We present a first global survey of ion jets within current sheets, which are thought of as a primary product of magnetic reconnection, with Mars Atmosphere and Volatile EvolutioN (MAVEN) data obtained on the dayside and nightside in the Martian magnetosphere. We develop a fully automated algorithm to efficiently and reliably identify ion flow velocity deviations during current sheet crossings recorded by MAVEN. The obtained statistical results show that (i) sunward and anti‐sunward ion jets embedded within current sheets are commonly and widely observed in the dayside and nightside magnetosphere of Mars, (ii) measured magnetic field profiles and inferred magnetic topology during the jet crossings are generally consistent with those expected for reconnecting current sheets, (iii) the jet occurrence is apparently independent of upstream driver conditions, and (iv) most of the identified current sheets are thin with a roughly ion scale half‐thickness and are embedded in low beta plasma, suggesting that they are seemingly capable of reconnection based on the known onset conditions. The widespread distributions and common detection of ion jets, as well as the ubiquitous presence of thin, low‐beta current sheets, imply that magnetic reconnection could operate frequently in many locations around Mars, thereby playing an important role in the dynamics of the Martian magnetosphere for most, if not all, of the solar wind conditions in the present epoch.

The Origin of Switchbacks in the Solar Corona: Linear Theory

Abstract The origin, structure, and propagation characteristics of a switchback are compelling questions posed by Parker Solar Probe (PSP) observations of velocity spikes and magnetic field reversals. By assuming interchange reconnection between coronal loop and open magnetic field, we show that this results in the generation of upward (into the heliosphere) and downward complex structures propagating at the fast magnetosonic speed (i.e., the Alfvén speed in the low plasma beta corona) that can have an arbitrary radial magnetic field deflection, including “S-shaped.” We derive the evolution equation for the switchback radial magnetic field as it propagates through the inhomogeneous supersonic solar corona. An analytic solution for arbitrary initial conditions is used to investigate the properties of a switchback propagating from launch ∼6 to ∼35 R ⊙ where PSP observed switchbacks during its first encounter. We provide a detailed comparison to an example event, showing that the magnetic field and plasma solutions are in accord with PSP observations. For a simple single switchback, the model predicts either a single or a double-humped structure; the former corresponding to PSP observing either the main body or the flanks of the switchback. The clustering of switchbacks and their sometimes complicated structure may be due to the formation of multiple closely spaced switchbacks created by interchange reconnection with numerous open and loop magnetic field lines over a short period. We show that their evolution yields a complex, aggregated group of switchbacks that includes “sheaths” with large-amplitude radial magnetic field and velocity fluctuations.

Tearing mode analysis in electron magnetohydrodynamics with pressure gradient

A two-dimensional three-component (2D/3C) electron magnetohydrodynamic (EMHD) model is implemented to investigate the linear behavior of collisionless tearing modes in slab geometry. Owing to nonuniformity of thermal pressure and plasma density, the electron diamagnetic drift and Biermann battery effects are involved. The linear structures, growth rate, and real frequency are analyzed with a thin current sheet in the electron inertia scale. The ratio of the electron current to the total current in equilibrium can notably promote the growth of the tearing mode in EMHD. More numerical results then show that the effect of the pressure gradient on the tearing mode is dependent on the plasma beta, stabilizing the mode in a low enough beta limit but destabilizing it with the higher beta. The frequency of the mode caused by the pressure gradient is found to be increasing with it. The Biermann battery effect slightly stabilizes the tearing mode in low beta plasma but is indicated to be significant in much higher beta conditions.

Two-fluid and kinetic transport physics of Kelvin–Helmholtz instabilities in nonuniform low-beta plasmas

Hall-magnetohydrodynamic (Hall-MHD) theory, two-fluid simulations, and kinetic simulations are used to investigate the cross-field transport properties of Kelvin–Helmholtz instabilities in nonuniform low-beta collisionless plasmas. Hall-MHD analysis shows how the linear properties of the instability are modified by density gradients and magnetization. High-order accurate two-fluid and kinetic simulations, with complete dynamics of finite-mass electrons and ions, are applied to a suite of parameter cases to systematically assess the effects of diamagnetic drift, magnetization, charge separation, and finite Larmor motion. Initialization of exact two-species kinetic equilibria facilitates the study of isolated physical effects and enables detailed cross-comparisons between two-fluid and kinetic simulations, including for cases where ion gyroradii are comparable to gradient scale lengths. For nonuniform plasmas with significant space charge, the results of two-fluid and kinetic simulations are found to disagree with Hall-MHD predictions. Kelvin–Helmholtz instability growth rates, per unit shear, are shown to be smaller when ion diamagnetic drift and E × B drift are parallel and larger when the two drifts are antiparallel. The effect is attributed to polarization drift in the shear layer, which leads to redistribution of charge, alters the electric field that drives plasma advection, and consequently modifies growth rates. Instability-induced mass transport for different parameters is characterized in terms of the flux across the shear layer and a simplified diffusion model. Distribution functions from kinetic simulations are shown to deviate substantially from Maxwellian reconstructions, indicating the importance of kinetic physics during the nonlinear phase of the instability.

Dynamic charging effect on electromagnetic wave in low $\beta $ plasmas

Low beta plasma system is chosen for the study of low frequency electromagnetic Shear Alfven waves exhibiting dynamical charging of micron/sub micron size dust particles. A semi classical Quantum Hydrodynamic Model (QHD) is employed which includes the quantum effects for plasma electrons in the terms of statistical degenerate pressure, Landau quantization, exchange-correlation potential, and tunneling potential. The low beta plasma characteristics provides the opportunity for the usage of two potential theory in order to derive a complex equation of dispersion of the shear Alfven wave in quantum dusty magneto plasmas. Analytical simplification of dispersion equation provides the damping rate of the Shear Alfven wave which is verified graphically for a typical set of parameters. It is noticed that the damping rate depends upon the quantum properties of Landau quantization, exchange-correlation, and tunneling potential in addition to the dust radius which modifies the floating potential of the dust particles. This work may seek its applications in the astrophysical objects in general, and in the semiconductor laboratory plasmas for particular dust size and magnetic field.

Onset of Magnetic Explosion in Solar Coronal Jets in Quiet Regions on the Central Disk

Abstract We examine the initiation of 10 coronal jet eruptions in quiet regions on the central disk, thereby avoiding near-limb spicule-forest obscuration of the slow-rise onset of the minifilament eruption. From the Solar Dynamics Observatory/Atmospheric Imaging Assembly 171 Å 12 s cadence movie of each eruption, we (1) find and compare the start times of the minifilament's slow rise, the jet-base bright point, the jet-base-interior brightening, and the jet spire, and (2) measure the minifilament's speed at the start and end of its slow rise. From (a) these data, (b) prior observations showing that each eruption was triggered by magnetic flux cancelation under the minifilament, and (c) the breakout-reconnection current sheet observed in one eruption, we confirm that quiet-region jet-making minifilament eruptions are miniature versions of CME-making filament eruptions, and surmise that in most quiet-region jets: (1) the eruption starts before runaway reconnection starts, (2) runaway reconnection does not start until the slow-rise speed is at least ∼1 km s−1, and (3) at and before eruption onset, there is no current sheet of appreciable extent. We therefore expect that (i) many CME-making filament eruptions are triggered by flux cancelation under the filament, (ii) emerging bipoles seldom, if ever, directly drive jet production because the emergence is seldom, if ever, fast enough, and (iii) at a separatrix or quasi-separatrix in any astrophysical setting of a magnetic field in low-beta plasma, a current sheet of appreciable extent can be built only dynamically by a magnetohydrodynamic convulsion of the field, not by quasi-static gradual converging of the field.

Entropy modes in multi-component plasmas confined by a dipole field

In this paper, the collisionless electrostatic instability of the entropy mode is investigated linearly in multi-component plasmas confined by a dipole magnetic field, which commonly exists in space environments, as well as laboratory plasmas, such as Collisionless Terrella eXperiment, Levitated Dipole eXperiment, and Dipole Research EXperiment. We focus on characteristics of the entropy mode driven by the plasma density or/and the temperature gradient at low plasma beta (=8πP0/B02). The theoretical analysis of this work agrees with the calculated results qualitatively. It is indicated that the peak growth rate of the instability is in the regime of k⊥ρi ∼ 1, and the entropy mode tends to be more stable as the percentage of the heavy ion increases. For multi-component plasmas, each component has the entropy mode feature of its own, instead of simply averaging all the components. While for the electron with an isotropic temperature, the use of weighted harmonic average can be a good approximation for simplification.

Thermal instabilities: Fragmentation and field misalignment of filament fine structure

Context. Prominences show a surprising amount of fine structure and it is widely believed that their threads, as seen in Hα observations, provide indirect information concerning magnetic field topology. Both prominence and coronal rain condensations most likely originate from thermal instabilities in the solar corona. It is still not understood how non-linear instability evolution shapes the observed fine structure of prominences. Investigating this requires multidimensional, high-resolution simulations to resolve all emerging substructure in great detail. Aims. We investigate the spontaneous emergence and evolution of fine structure in high-density condensations formed through the process of thermal instability under typical solar coronal conditions. Our study reveals intricate multidimensional processes that occur through in situ condensations in a representative coronal volume in a low plasma beta regime. Methods. We quantified slow wave eigenfunctions used as perturbations and discuss under which conditions the thermal mode is unstable when anisotropic thermal conduction effects are included. We performed 2D and 3D numerical simulations of interacting slow magnetohydrodynamic (MHD) wave modes when all relevant non-adiabatic effects are included. Multiple levels of adaptive mesh refinement achieve a high resolution near regions with high density, thereby resolving any emerging fine structure automatically. Our study employs a local periodic coronal region traversed by damped slow waves inspired by the presence of such waves observed in actual coronal magnetic structures. Results. We show that the interaction of multiple slow MHD wave modes in a regime unstable to the thermal mode leads to thermal instability. This initially forms pancake-like structures almost orthogonal to the local magnetic field, while low-pressure induced inflows of matter generate rebound shocks. This is succeeded by the rapid disruption of these pancake sheets through thin-shell instabilities evolving naturally from minute ram pressure imbalances. This eventually creates high-density blobs accompanied by thread-like features from shear flow effects. The further evolution of the blobs follows the magnetic field lines, such that a dynamical realignment with the background magnetic field appears. However, the emerging thread-like features are not at all field-aligned, implying only a very weak link between fine structure orientation and magnetic field topology. Conclusions. As seen in our synthetic Hα views, threads formed by non-linear thermal instability evolution do not strictly outline magnetic field structure and this finding has far-reaching implications for field topology interpretations based on Hα observations.

Suppression of electrostatic micro-instabilities in maximum-J stellarators

We demonstrate favorable stability properties of maximum-J stellarators, exemplified by the Wendelstein 7-X (W7-X) device, in scenarios with low plasma beta. A large number of electrostatic linear gyrokinetic simulations are conducted to scan the relevant parameter space for different configurations, resulting in stability maps that account for the key micro-instabilities thought to drive turbulent transport. These maps exhibit a ‘stability valley’ in the region where the normalized ion temperature gradient is roughly equal to the normalized density gradient. In this valley, the electrostatic instabilities are partly suppressed thanks to the maximum-J property of the W7-X field. This property varies across different W7-X configurations, and this measurable difference is demonstrated to affect the size of the stability valley. Finally, the impact of the isotope effect and collisions on the valley is examined.

Local Regimes of Turbulence in 3D Magnetic Reconnection

Abstract The process of magnetic reconnection when studied in nature or when modeled in 3D simulations differs in one key way from the standard 2D paradigmatic cartoon: it is accompanied by many fluctuations in the electromagnetic fields and plasma properties. We developed a diagnostics to study the spectrum of fluctuations in the various regions around a reconnection site. We define the regions in terms of the local value of the flux function that determines the distance from the reconnection site, with positive values in the outflow and negative values in the inflow. We find that fluctuations belong to two very different regimes depending on the local plasma beta (defined as the ratio of plasma and magnetic pressures). The first regime develops in the reconnection outflows where beta is high and it is characterized by a strong link between plasma and electromagnetic fluctuations, leading to momentum and energy exchanges via anomalous viscosity and resistivity. But there is a second, low-beta regime: it develops in the inflow and in the region around the separatrix surfaces, including the reconnection electron diffusion region itself. It is remarkable that this low-beta plasma, where the magnetic pressure dominates, remains laminar even though the electromagnetic fields are turbulent.

Whistler instability stimulated by the suprathermal electrons present in space plasmas

In the absence of efficient collisions, deviations from thermal equilibrium of plasma particle distributions are controlled by the self-generated instabilities. The whistler instability is a notorious example, usually responsible for the regulation of electron temperature anisotropy $A = T_{\perp}/T_\parallel>$ (with $\perp, \parallel$ respective to the magnetic field direction) observed in space plasmas, e.g., solar wind and planetary magnetospheres. Suprathermal electrons present in these environments change the plasma dispersion and stability properties, with expected consequences on the kinetic instabilities and the resulting fluctuations, which, in turn, scatter the electrons and reduce their anisotropy. In order to capture these mutual effects we use a quasilinear kinetic approach and PIC simulations, which provide a comprehensive characterization of the whistler instability under the influence of suprathermal electrons. Analysis is performed for a large variety of plasma conditions, ranging from low-beta plasmas encountered in outer corona or planetary magnetospheres to a high-beta solar wind characteristic to large heliospheric distances. Enhanced by the suprathermal electrons, whistler fluctuations stimulate the relaxation of temperature anisotropy, and this influence of suprathermals increases with plasma beta parameter.