Small-scale Current Sheets Research Articles

Small-scale Current Sheets and Associated Switchback Activity in the Inner Heliosphere

Several long-standing theories postulate that turbulent dissipation can heat solar wind protons in situ. Turbulent dissipation can occur via current sheets, which are small-scale structures embedded in the solar wind magnetic field. This study examines the role that switchbacks—intermediate-scale reversals in the interplanetary magnetic field—may play in heating the solar wind by generating current sheets. We explore this possible relationship by analyzing the characteristics of current sheets within and around switchback regions. Previous studies investigated current sheet properties during Parker Solar Probe's first solar encounter, analyzed current sheets using a wide range of statistics, and explored trends that switchbacks follow with radial distance from the Sun. The present study builds on these works by analyzing the distribution and maximum values of solar wind current sheets using the Partial Variance of Increments method and focusing on how these properties correlate with the presence of switchbacks to better understand how switchbacks contribute to current sheet activity. We conclude that there are no increased current sheet populations observed within and around switchbacks, with most current sheets being observed outside switchbacks. We find a consistent distribution of current sheets regardless of whether there is concurrent switchback activity. We also observe that current sheets follow a uniform occurrence rate with increased distance from the Sun, while switchback regions significantly evolve with larger radial distances. Our findings suggest that local turbulence may be responsible for generating solar wind current sheets and does so with the same efficiency inside and outside of switchback regions.

Coronal energy release by MHD avalanches

Context. Magnetohydrodynamic (MHD) instabilities, such as the kink instability, can trigger the chaotic fragmentation of a twisted magnetic flux tube into small-scale current sheets that dissipate as aperiodic impulsive heating events. In turn, the instability could propagate as an avalanche to nearby flux tubes and lead to a nanoflare storm. Our previous work was devoted to related 3D MHD numerical modeling, which included a stratified atmosphere from the solar chromosphere to the corona, tapering magnetic field, and solar gravity for curved loops with the thermal structure modelled by plasma thermal conduction, along with optically thin radiation and anomalous resistivity for 50 Mm flux tubes. Aims. Using 3D MHD modeling, this work addresses predictions for the extreme-ultraviolet (EUV) imaging spectroscopy of such structure and evolution of a loop, with an average temperature of 2–2.5 MK in the solar corona. We set a particular focus on the forthcoming MUSE mission, as derived from the 3D MHD modeling. Methods. From the output of the numerical simulations, we synthesized the intensities, Doppler shifts, and non-thermal line broadening in 3 EUV spectral lines in the MUSE passbands: Fe IX 171 Å, Fe XV 284 Å, and Fe XIX 108 Å, emitted by ∼1 MK, ∼2 MK, and ∼10 MK plasma, respectively. These data were detectable by MUSE, according to the MUSE expected pixel size, temporal resolution, and temperature response functions. We provide maps showing different view angles (front and top) and realistic spectra. Finally, we discuss the relevant evolutionary processes from the perspective of possible observations. Reults. We find that the MUSE observations might be able to detect the fine structure determined by tube fragmentation. In particular, the Fe IX line is mostly emitted at the loop footpoints, where we might be able to track the motions that drive the magnetic stressing and detect the upward motion of evaporating plasma from the chromosphere. In Fe XV, we might see the bulk of the loop with increasing intensity, with alternating filamentary Doppler and non-thermal components in the front view, along with more defined spots in the topward view. The Fe XIX line is very faint within the chosen simulation parameters; thus, any transient brightening around the loop apex may possibly be emphasized by the folding of sheet-like structures, mainly at the boundary of unstable tubes. Conclusions. In conclusion, we show that coronal loop observations with MUSE can pinpoint some crucial features of MHD-modeled ignition processes, such as the related dynamics, helping to identify the heating processes.

INFLUENCE OF A UNIFORM MAGNETIC FIELD ON THE GENERATION OF STRONG SMALL-SCALE MAGNETIC FIELDS DURING THE INJECTION OF A PLASMA WITH HOT ELECTRONS INTO AN INHOMOGENEOUS COLD PLASMA LAYER

We carry out a numerical modeling of plasma injection with hot electrons into a thin layer of cold plasma in the presence of an external magnetic field. We show that the latter can significantly affect the emerging small-scale current filaments and sheets, even if it does not magnetize the particles and does not change the overall dynamics of the redistribution of the total plasma density in the process of injection. The effect observed depends on the orientation of the external magnetic field that is parallel to the plane that bounds the cold plasma layer, if the injection occurs from a narrow strip lying in this plane. In this situation, which corresponds to the ablation of a flat target by a femtosecond laser beam using cylindrical focusing, we study the evolution of the characteristic structures of the formed small-scale magnetic field. It is established that its generation is associated with instabilities of the anisotropic velocity distribution of electrons, and its value can be many times greater than the value of the external magnetic field.

Electron magnetic reconnection without ion coupling in Earth's turbulent magnetosheath.

Magnetic reconnection in current sheets is a magnetic-to-particle energy conversion process that is fundamental to many space and laboratory plasma systems. In the standard model of reconnection, this process occurs in a minuscule electron-scale diffusion region1,2. On larger scales, ions couple to the newly reconnected magnetic-field lines and are ejected away from the diffusion region in the form of bi-directional ion jets at the ion Alfvén speed3-5. Much of the energy conversion occurs in spatially extended ion exhausts downstream of the diffusion region 6 . In turbulent plasmas, which contain a large number of small-scale current sheets, reconnection has long been suggested to have a major role in the dissipation of turbulent energy at kinetic scales7-11. However, evidence for reconnection plasma jetting in small-scale turbulent plasmas has so far been lacking. Here we report observations made in Earth's turbulent magnetosheath region (downstream of the bow shock) of an electron-scale current sheet in which diverging bi-directional super-ion-Alfvénic electron jets, parallel electric fields and enhanced magnetic-to-particle energy conversion were detected. Contrary to the standard model of reconnection, the thin reconnecting current sheet was not embedded in a wider ion-scale current layer and no ion jets were detected. Observations of this and other similar, but unidirectional, electron jet events without signatures of ion reconnection reveal a form of reconnection that can drive turbulent energy transfer and dissipation in electron-scale current sheets without ion coupling.

TURBULENCE-GENERATED PROTON-SCALE STRUCTURES IN THE TERRESTRIAL MAGNETOSHEATH

ABSTRACT Recent results of numerical magnetohydrodynamic simulations suggest that in collisionless space plasmas, turbulence can spontaneously generate thin current sheets. These coherent structures can partially explain the intermittency and the non-homogenous distribution of localized plasma heating in turbulence. In this Letter, Cluster multi-point observations are used to investigate the distribution of magnetic field discontinuities and the associated small-scale current sheets in the terrestrial magnetosheath downstream of a quasi-parallel bow shock. It is shown experimentally, for the first time, that the strongest turbulence-generated current sheets occupy the long tails of probability distribution functions associated with extremal values of magnetic field partial derivatives. During the analyzed one-hour time interval, about a hundred strong discontinuities, possibly proton-scale current sheets, were observed.

DIFFUSIVE SHOCK ACCELERATION AND RECONNECTION ACCELERATION PROCESSES

Shock waves, as shown by simulations and observations, can generate high levels of downstream vortical turbulence, including magnetic islands. We consider a combination of diffusive shock acceleration (DSA) and downstream magnetic-island-reconnection-related processes as an energization mechanism for charged particles. Observations of electron and ion distributions downstream of interplanetary shocks and the heliospheric termination shock (HTS) are frequently inconsistent with the predictions of classical DSA. We utilize a recently developed transport theory for charged particles propagating diffusively in a turbulent region filled with contracting and reconnecting plasmoids and small-scale current sheets. Particle energization associated with the anti-reconnection electric field, a consequence of magnetic island merging, and magnetic island contraction, are considered. For the former only, we find that (i) the spectrum is a hard power law in particle speed, and (ii) the downstream solution is constant. For downstream plasmoid contraction only, (i) the accelerated spectrum is a hard power law in particle speed; (ii) the particle intensity for a given energy peaks downstream of the shock, and the distance to the peak location increases with increasing particle energy, and (iii) the particle intensity amplification for a particular particle energy, f(x,c/c_0)/f(0,c/c_0), is not 1, as predicted by DSA, but increases with increasing particle energy. The general solution combines both the reconnection-induced electric field and plasmoid contraction. The observed energetic particle intensity profile observed by Voyager 2 downstream of the HTS appears to support a particle acceleration mechanism that combines both DSA and magnetic-island-reconnection-related processes.

Particle acceleration by combined diffusive shock acceleration and downstream multiple magnetic island acceleration

As a consequence of the evolutionary conditions [28; 29], shock waves can generate high levels of downstream vortical turbulence. Simulations [32-34] and observations [30; 31] support the idea that downstream magnetic islands (also called plasmoids or flux ropes) result from the interaction of shocks with upstream turbulence. Zank et al. [18] speculated that a combination of diffusive shock acceleration (DSA) and downstream reconnection-related effects associated with the dynamical evolution of a “sea of magnetic islands” would result in the energization of charged particles. Here, we utilize the transport theory [18; 19] for charged particles propagating diffusively in a turbulent region filled with contracting and reconnecting plasmoids and small-scale current sheets to investigate a combined DSA and downstream multiple magnetic island charged particle acceleration mechanism. We consider separately the effects of the anti-reconnection electric field that is a consequence of magnetic island merging [17], and magnetic island contraction [14]. For the merging plasmoid reconnection- induced electric field only, we find i) that the particle spectrum is a power law in particle speed, flatter than that derived from conventional DSA theory, and ii) that the solution is constant downstream of the shock. For downstream plasmoid contraction only, we find that i) the accelerated particle spectrum is a power law in particle speed, flatter than that derived from conventional DSA theory; ii) for a given energy, the particle intensity peaks downstream of the shock, and the peak location occurs further downstream of the shock with increasing particle energy, and iii) the particle intensity amplification for a particular particle energy, f(x, c/c0)/f(0, c/c0), is not 1, as predicted by DSA theory, but increases with increasing particle energy. These predictions can be tested against observations of electrons and ions accelerated at interplanetary shocks and the heliospheric termination shock.

ENERGY DISSIPATION IN MAGNETIC NULL POINTS AT KINETIC SCALES

We use kinetic particle-in-cell and magnetohydrodynamic simulations supported by an observational dataset to investigate magnetic reconnection in clusters of null points in space plasma. The magnetic configuration under investigation is driven by fast adiabatic flux rope compression that dissipates almost half of the initial magnetic field energy. In this phase powerful currents are excited producing secondary instabilities, and the system is brought into a state of `intermittent turbulence' within a few ion gyro-periods. Reconnection events are distributed all over the simulation domain and energy dissipation is rather volume-filling. Numerous spiral null points interconnected via their spines form null lines embedded into magnetic flux ropes; null point pairs demonstrate the signatures of torsional spine reconnection. However, energy dissipation mainly happens in the shear layers formed by adjacent flux ropes with oppositely directed currents. In these regions radial null pairs are spontaneously emerging and vanishing, associated with electron streams and small-scale current sheets. The number of spiral nulls in the simulation outweighs the number of radial nulls by a factor of 5\---10, in accordance with Cluster observations in the Earth's magnetosheath. Twisted magnetic fields with embedded spiral null points might indicate the regions of major energy dissipation for future space missions such as Magnetospheric Multiscale Mission (MMS).

A dynamical model of plasma turbulence in the solar wind.

A dynamical approach, rather than the usual statistical approach, is taken to explore the physical mechanisms underlying the nonlinear transfer of energy, the damping of the turbulent fluctuations, and the development of coherent structures in kinetic plasma turbulence. It is argued that the linear and nonlinear dynamics of Alfvén waves are responsible, at a very fundamental level, for some of the key qualitative features of plasma turbulence that distinguish it from hydrodynamic turbulence, including the anisotropic cascade of energy and the development of current sheets at small scales. The first dynamical model of kinetic turbulence in the weakly collisional solar wind plasma that combines self-consistently the physics of Alfvén waves with the development of small-scale current sheets is presented and its physical implications are discussed. This model leads to a simplified perspective on the nature of turbulence in a weakly collisional plasma: the nonlinear interactions responsible for the turbulent cascade of energy and the formation of current sheets are essentially fluid in nature, while the collisionless damping of the turbulent fluctuations and the energy injection by kinetic instabilities are essentially kinetic in nature.

COLLISIONLESS DAMPING AT ELECTRON SCALES IN SOLAR WIND TURBULENCE

The dissipation of turbulence in the weakly collisional solar wind plasma is governed by unknown kinetic mechanisms. Two candidates have been suggested to play an important role in the dissipation, collisionless damping via wave-particle interactions and dissipation in small-scale current sheets. High resolution spacecraft measurements of the turbulent magnetic energy spectrum provide important constraints on the dissipation mechanism. The limitations of popular fluid and hybrid numerical schemes for simulation of the dissipation of solar wind turbulence are discussed, and instead a three-dimensional kinetic approach is recommended. We present a three-dimensional nonlinear gyrokinetic simulation of solar wind turbulence at electron scales that quantitatively reproduces the exponential form of the turbulent magnetic energy spectrum measured in the solar wind. A weakened cascade model that accounts for nonlocal interactions and collisionless Landau damping also quantitatively agrees with the observed exponential form. These results establish that a turbulent cascade of kinetic Alfven waves that is terminated by collisionless Landau damping is sufficient to explain the observed magnetic energy spectrum in the dissipation range of solar wind turbulence.

CURRENT SHEETS AND COLLISIONLESS DAMPING IN KINETIC PLASMA TURBULENCE

We present the first study of the formation and dissipation of current sheets at electron scales in a wave-driven, weakly collisional, 3D kinetic turbulence simulation. We investigate the relative importance of dissipation associated with collisionless damping via resonant wave-particle interactions versus dissipation in small-scale current sheets in weakly collisional plasma turbulence. Current sheets form self-consistently from the wave-driven turbulence, and their filling fraction is well correlated to the electron heating rate. However, the weakly collisional nature of the simulation necessarily implies that the current sheets are not significantly dissipated via Ohmic dissipation. Rather, collisionless damping via the Landau resonance with the electrons is sufficient to account for the measured heating as a function of scale in the simulation, without the need for significant Ohmic dissipation. This finding suggests the possibility that the dissipation of the current sheets is governed by resonant wave-particle interactions and that the locations of current sheets correspond spatially to regions of enhanced heating.

KINETIC PLASMA TURBULENCE IN THE FAST SOLAR WIND MEASURED BYCLUSTER

The k-filtering technique and wave polarization analysis are applied to Cluster magnetic field data to study plasma turbulence at the scale of the ion gyroradius in the fast solar wind. Waves are found propagating in directions nearly perpendicular to the background magnetic field at such scales. The frequencies of these waves in the solar wind frame are much smaller than the proton gyro-frequency. After the wave vector ${\bf k}$ is determined at each spacecraft frequency $f_{sc}$, wave polarization property is analyzed in the plane perpendicular to ${\bf k}$. Magnetic fluctuations have $\delta B_\perp>\delta B_\parallel$ (here the $\parallel$ and $\perp$ refer to the background magnetic field ${\bf B}_0$). The wave magnetic field has right-handed polarization at propagation angles $\theta_{\bf kB}<90^\circ$ and $>90^\circ$. The magnetic field in the plane perpendicular to ${\bf B}_0$ however has no clear sense of a dominant polarization but local rotations. We discuss the merits and limitations of linear kinetic Aflv\'en waves (KAWs) and coherent Alfv\'en vortices in the interpretation of the data. We suggest that the fast solar wind turbulence may be populated with KAWs, small scale current sheets and Alfv\'en vortices at ion kinetic scales.

SPONTANEOUS CURRENT-LAYER FRAGMENTATION AND CASCADING RECONNECTION IN SOLAR FLARES. I. MODEL AND ANALYSIS

Magnetic reconnection is commonly considered as a mechanism of solar (eruptive) flares. A deeper study of this scenario reveals, however, a number of open issues. Among them is the fundamental question, how the magnetic energy is transferred from large, accumulation scales to plasma scales where its actual dissipation takes place. In order to investigate this transfer over a broad range of scales we address this question by means of a high-resolution MHD simulation. The simulation results indicate that the magnetic-energy transfer to small scales is realized via a cascade of consecutive smaller and smaller flux-ropes (plasmoids), in analogy with the vortex-tube cascade in (incompressible) fluid dynamics. Both tearing and (driven) "fragmenting coalescence" processes are equally important for the consecutive fragmentation of the magnetic field (and associated current density) to smaller elements. At the later stages a dynamic balance between tearing and coalescence processes reveals a steady (power-law) scaling typical for cascading processes. It is shown that cascading reconnection also addresses other open issues in solar flare research such as the duality between the regular large-scale picture of (eruptive) flares and the observed signatures of fragmented (chaotic) energy release, as well as the huge number of accelerated particles. Indeed, spontaneous current-layer fragmentation and formation of multiple channelised dissipative/acceleration regions embedded in the current layer appears to be intrinsic to the cascading process. The multiple small-scale current sheets may also facilitate the acceleration of a large number of particles. The structure, distribution and dynamics of the embedded potential acceleration regions in a current layer fragmented by cascading reconnection are studied and discussed.

SHEAR PHOTOSPHERIC FORCING AND THE ORIGIN OF TURBULENCE IN CORONAL LOOPS

We present a series of numerical simulations aimed at understanding the nature and origin of turbulence in coronal loops in the framework of the Parker model for coronal heating. A coronal loop is studied via reduced magnetohydrodynamic (MHD) simulations in Cartesian geometry. A uniform and strong magnetic field threads the volume between the two photospheric planes, where a velocity field in the form of a one-dimensional shear flow pattern is present. Initially, the magnetic field that develops in the coronal loop is a simple map of the photospheric velocity field. This initial configuration is unstable to a multiple tearing instability that develops islands with X and O points in the plane orthogonal to the axial field. Once the nonlinear stage sets in the system evolution is characterized by a regime of MHD turbulence dominated by magnetic energy. A well-developed power law in energy spectra is observed and the magnetic field never returns to the simple initial state mapping the photospheric flow. The formation of X and O points in the planes orthogonal to the axial field allows the continued and repeated formation and dissipation of small-scale current sheets where the plasma is heated. We conclude that the observed turbulent dynamics are not induced by the complexity of the pattern that the magnetic field-line footpoints follow but they rather stem from the inherent nonlinear nature of the system.

THE HEATING OF TEST PARTICLES IN NUMERICAL SIMULATIONS OF ALFVÉNIC TURBULENCE

We study the heating of charged test particles in three-dimensional numerical simulations of weakly compressible magnetohydrodynamic (MHD) turbulence (``Alfvenic turbulence''); these results are relevant to particle heating and acceleration in the solar wind, solar flares, accretion disks onto black holes, and other astrophysics and heliospheric environments. The physics of particle heating depends on whether the gyrofrequency of a particle is comparable to the frequency of a turbulent fluctuation that is resolved on the computational domain. Particles with these frequencies nearly equal undergo strong perpendicular heating (relative to the local magnetic field) and pitch angle scattering. By contrast, particles with large gyrofrequency undergo strong parallel heating. Simulations with a finite resistivity produce additional parallel heating due to parallel electric fields in small-scale current sheets. Many of our results are consistent with linear theory predictions for the particle heating produced by the Alfven and slow magnetosonic waves that make up Alfvenic turbulence. However, in contrast to linear theory predictions, energy exchange is not dominated by discrete resonances between particles and waves; instead, the resonances are substantially ``broadened.'' We discuss the implications of our results for solar and astrophysics problems, in particular the thermodynamics of the near-Earth solar wind. We conclude that Alfvenic turbulence produces significant parallel heating via the interaction between particles and magnetic field compressions (``slow waves''). However, on scales above the proton Larmor radius, Alfvenic turbulence does not produce significant perpendicular heating of protons or minor ions.

Explosive reconnection in magnetars

X-ray activity of Anomalous X-ray Pulsars and Soft Gamma-Ray Repeaters may result from the heating of their magnetic corona by direct currents dissipated by magnetic reconnection. We investigate the possibility that X-ray flares and bursts observed from AXPs and SGRs result from magnetospheric reconnection events initiated by development of tearing mode in magnetically-dominated relativistic plasma. We formulate equations of resistive force-free electrodynamics, discuss its relation to ideal electrodynamics, and give examples of both ideal and resistive equilibria. Resistive force-free current layers are unstable toward the development of small-scale current sheets where resistive effects become important. Thin current sheets are found to be unstable due to the development of resistive force-free tearing mode. The growth rate of tearing mode is intermediate between the short \Alfven time scale $ \tau_A$ and a long resistive time scale $\tau_R$: $\Gamma \sim 1/(\tau_R \tau_A)^{1/2}$, similar to the case of non-relativistic non-force-free plasma. We propose that growth of tearing mode is related to the typical rise time of flares, $\sim 10$ msec. Finally, we discuss how reconnection may explain other magnetar phenomena and ways to test the model.

Stationary magnetospheric convection on November 24, 1981. 2. Small-scale structures in the dayside cusp/cleft

Abstract. A case study of the dayside cusp/cleft region during an interval of stationary magnetospheric convection (SMC) on November, 24, 1981 is presented, based on detailed measurements made by the AUREOL-3 satellite. Layered small-scale field-aligned current sheets, or loops, superimposed to a narrow V-shaped ion dispersion structure, were observed just equatorward from the region of the "cusp proper". The equatorward sheet was accompanied by a very intense and short (less than 1 s) ion intensity spike at 100 eV. No major differences were noted of the characteristics of the LLBL, or "boundary cusp", and plasma mantle precipitation during this SMC period from those typical of the cusp/cleft region for similar IMF conditions. Simultaneous NOAA-6 and NOAA-7 measurements described in Despirak et al. were used to estimate the average extent of the "cusp proper" (defined by dispersed precipitating ions with the energy flux exceeding 10-3 erg cm-2 s-1) during the SMC period, as ~0.73° ILAT width, 2.6-3.4 h in MLT, and thus the recently merged magnetic flux, 0.54-0.70 × 107 Wb. This, together with the average drift velocity across the cusp at the convection throat, ~0.5 km s-1, allowed to evaluate the cusp merging contribution to the total cross-polar cap potential difference, ~33.8-43.8 kV. It amounts to a quite significant part of the total cross-polar cap potential difference evaluated from other data. A "shutter" scenario is suggested for the ion beam injection/penetration through the stagnant plasma region in the outer cusp to explain the pulsating nature of the particle injections in the low- and medium-altitude cusp region.Key words. Magnetospheric physics (current systems; magnetopause · cusp · and boundary layers; solar wind-magnetosphere interactions).

Quiet auroral arcs: Ionosphere effect of magnetospheric convection stratification

A detailed study of the mechanism of electromagnetic stratification of the large-scale stationary magnetospheric convection due to a friction of the convective flow in the ionosphere layer was performed. Magnetosphere-ionosphere interaction was taken into account by means of the effective boundary conditions on the ionosphere top and bottom boundaries including the actual height profile of charge particles velocity in the ionosphere. It has been shown that the magnetospheric convection is stratified into small-scale current sheets which are respective in the linear approximation to an oblique Alfvén wave. The dispersion equation was deduced for the Alfvén mode and its solution obtained determining the space-time scales and the increment of instability. The maximum increment is realized for the disturbances stretched along the convection velocity that is correspondent to the actual orientation of the auroral arcs. In the conditions of rapid growth of Alfvén velocity above the maximum of the ionosphere F layer, it was shown that small-scale disturbances with the transverse scales l ⊥ ≲ 1 km are localized at the altitudes up to several thousand kilometers whereas the large-scale stratification penetrate into the equatorial plane of the magnetosphere. A mechanism is proposed to intensify the parallel electric field acting at that stratification stage when the field-aligned currents in the Alfvén wave are sufficient to form abnormal resistance along geomagnetic lines of force.