Search

Alongside magnetic reconnection, turbulence is another fundamental nonlinear plasma phenomenon that plays a key role in energy transport and conversion in space and astrophysical plasmas. From a numerical, theoretical, and observational point of view there is a long history of exploring the interplay between these two phenomena in space plasma environments; however, recent high-resolution, multi-spacecraft observations have ushered in a new era of understanding this complex topic. The interplay between reconnection and turbulence is both complex and multifaceted, and can be viewed through a number of different interrelated lenses - including turbulence acting to generate current sheets that undergo magnetic reconnection (turbulence-driven reconnection), magnetic reconnection driving turbulent dynamics in an environment (reconnection-driven turbulence) or acting as an intermediate step in the excitation of turbulence, and the random diffusive/dispersive nature of magnetic field lines embedded in turbulent fluctuations enabling so-called stochastic reconnection. In this paper, we review the current state of knowledge on these different facets of the interplay between turbulence and reconnection in the context of collisionless plasmas, such as those found in many near-Earth astrophysical environments, from a theoretical, numerical, and observational perspective. Particular focus is given to several key regions in Earth's magnetosphere - Earth's magnetosheath, magnetotail, and Kelvin-Helmholtz vortices on the magnetopause flanks - where NASA's Magnetospheric Multiscale mission has been providing new insights on the topic., Comment: Chapter 2.2 of ISSI Book on Magnetic Reconnection, submitted to Space Science Reviews

How ions are energized and heated is a fundamental problem in the study of energy dissipation in magnetized plasmas. In particular, the heating of heavy ions (including ${}^{4}\mathrm{He}^{2+}$, ${}^{3}\mathrm{He}^{2+}$ and others) has been a constant concern for understanding the microphysics of impulsive solar flares. In this article, via two-dimensional hybrid-kinetic Particle-in-Cell simulations, we study the heating of Helium ions (${}^{4}\mathrm{He}^{2+}$) by turbulence driven by cascading waves launched at large scales from the left-handed polarized Helium ion cyclotron wave branch of a multi-ion plasma composed of electrons, protons, and Helium ions. We find significant parallel (to the background magnetic field) heating for both Helium ions and protons due to the formation of beams and plateaus in their velocity distribution functions along the background magnetic field. The heating of Helium ions in the direction perpendicular to the magnetic field starts with a lower rate than that in the parallel direction, but overtakes the parallel heating after a few hundreds of the proton gyro-periods due to cyclotron resonances with mainly obliquely propagating waves induced by the cascade of injected Helium ion cyclotron waves at large scales. There is however little evidence for proton heating in the perpendicular direction due to the absence of left-handed polarized cyclotron waves near the proton cyclotron frequency. Our results are useful for understanding the preferential heating of ${}^{3}\mathrm{He}$ and other heavy ions in the ${}^{3}\mathrm{He}$-rich solar energetic particle events, in which Helium ions play a crucial role as a species of background ions regulating the kinetic plasma behavior., Comment: accepted by ApJ

Solar Orbiter is the first space mission observing the solar plasma both in situ and remotely, from a close distance, in and out of the ecliptic. The ultimate goal is to understand how the Sun produces and controls the heliosphere, filling the Solar System and driving the planetary environments. With six remote-sensing and four in-situ instrument suites, the coordination and planning of the operations are essential to address the following four top-level science questions: (1) What drives the solar wind and where does the coronal magnetic field originate? (2) How do solar transients drive heliospheric variability? (3) How do solar eruptions produce energetic particle radiation that fills the heliosphere? (4) How does the solar dynamo work and drive connections between the Sun and the heliosphere? Maximising the mission's science return requires considering the characteristics of each orbit, including the relative position of the spacecraft to Earth (affecting downlink rates), trajectory events (such as gravitational assist manoeuvres), and the phase of the solar activity cycle. Furthermore, since each orbit's science telemetry will be downloaded over the course of the following orbit, science operations must be planned at mission level, rather than at the level of individual orbits. It is important to explore the way in which those science questions are translated into an actual plan of observations that fits into the mission, thus ensuring that no opportunities are missed. First, the overarching goals are broken down into specific, answerable questions along with the required observations and the so-called Science Activity Plan (SAP) is developed to achieve this. The SAP groups objectives that require similar observations into Solar Orbiter Observing Plans (SOOPs), resulting in a strategic, top-level view of the optimal opportunities for science observations during the mission lifetime., Comment: 20 pages, 1 figure, accepted by Astronomy & Astrophysics

Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events., Comment: 4 pages (including a title page), white paper submitted to Helio2050 workshop at https://www.hou.usra.edu/meetings/helio2050/. arXiv admin note: substantial text overlap with arXiv:2004.00079

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process., Comment: 6 pages, 1 figure, white paper submitted to both Plasma 2020 and Astro 2020 Decadal Surveys

We analyzed properties of waves excited by mildly relativistic electron beams propagating along magnetic field with a ring-shape perpendicular momentum distribution in neutral and current-free solar coronal plasmas. These plasmas are subject to both the beam and the electron cyclotron maser (ECM) instabilities driven by the positive momentum gradient of the ring-beam electron distribution in the directions parallel and perpendicular to the ambient magnetic field, respectively. To explore the related kinetic processes self-consistently, 2.5-dimensional fully kinetic particle-in-cell (PIC) simulations were carried out. To quantify excited wave properties in different coronal conditions, we investigated the dependence of their energy and polarization on the ring-beam electron density and magnetic field. In general, electrostatic waves dominate the energetics of waves and nonlinear waves are ubiquitous. In weakly magnetized plasmas, where the electron cyclotron frequency $\omega_{ce}$ is lower than the electron plasma frequency $\omega_{pe}$, it is difficult to produce escaping electromagnetic waves with frequency $\omega > \omega_{pe}$ and small refractive index $|c k / \omega| < 1$ ($k$ and $c$ are the wavenumber and the light speed, respectively). Highly polarized and anisotropic escaping electromagnetic waves can, however, be effectively excited in strongly magnetized plasmas with $\omega_{ce}/\omega_{pe} \geq 1$. The anisotropy of the energy, circular polarization degree (CPD), and spectrogram of these escaping electromagnetic waves strongly depend on the number density ratio of the ring-beam electrons to the background electrons. In particular, their CPDs can vary from left-handed to right-handed with the decrease of the ring-beam density, which may explain some observed properties of solar radio bursts (e.g., radio spikes) from the solar corona., Comment: ApJ, accepted

In space and astrophysical plasmas, like in planetary magnetospheres, as that of Mercury,energetic electrons are often found near current sheets (CSs), which hints at electron acceleration by magnetic reconnection. Unfortunately, electron acceleration by reconnection is not well understood, yet. In particular, acceleration by turbulent plasmoid reconnection. We have investigated electron acceleration by turbulent plasmoid reconnection, described by MHD simulations, via test particle calculations. In order to avoid resolving all relevant turbulence scales down to the dissipation scales, a mean-field turbulence model is used to describe the turbulence of sub-grid scales (SGS) and their effects via a turbulent electromotive force (EMF). The mean-field model describes the turbulent EMF as a function of the mean values of current density, vorticity, magnetic field as well as of the energy, cross-helicity and residual helicity of the turbulence. We found that, mainly around X-points of turbulent reconnection, strongly enhanced localized EMFs most efficiently accelerated electrons and caused the formation of power-law spectra. Magnetic-field-aligned EMFs, caused by the turbulence, dominate the electron acceleration process. Scaling the acceleration processes to parameters of the Hermean magnetotail, electron energies up to 60 keV can be reached by turbulent plasmoid reconnection through the thermal plasma., Comment: 2018PhPl...25d2904Z

We report a newly found two-stage mechanism of electron acceleration near X-lines of 3D collisionless guide-field magnetic reconnection in the non-relativistic regime typical, e.g., for stellar coronae. We found that after electrons are first pre-accelerated during the linear growth of reconnection, they become additionally accelerated in the course of the nonlinear stage of 3D guide-field magnetic reconnection. This additional acceleration is due to the filamentation of electric and magnetic fields caused by streaming instabilities. In addition to enhanced parallel electric fields, the filamentation leads to additional curvature-driven electron acceleration in the guide-field direction. As a result, part of the the accelerated electron spectra becomes a power law with a spectral index of $\sim-1.6$ near the X-line. This second stage of acceleration due to nonlinear reconnection is relevant for the production of energetic electrons in, e.g., thin current sheets of stellar coronae., Comment: 15 pages. 11 figures. Revised version. Published in The Astrophysical Journal

Although turbulence has been conjectured to be important for magnetic reconnection, still very little is known about its role in collisionless plasmas. Previous attempts to quantify the effect of turbulence on reconnection usually prescribed Alfv\'enic or other low-frequency fluctuations or investigated collisionless kinetic effects in just two-dimensional configurations and antiparallel magnetic fields. In view of this, we analyzed the kinetic turbulence self-generated by three-dimensional guide-field reconnection through force-free current sheets in frequency and wavenumber spaces, utilizing 3D particle-in-cell code numerical simulations. Our investigations reveal reconnection rates and kinetic turbulence with features similar to those obtained by current in-situ spacecraft observations of MMS as well as in the laboratory reconnection experiments MRX, VTF and \textsc{Vineta}-II. In particular we found that the kinetic turbulence developing in the course of 3D guide-field reconnection exhibits a broadband power-law spectrum extending beyond the lower-hybrid frequency and up to the electron frequencies. In the frequency space the spectral index of the turbulence appeared to be close to -2.8 at the reconnection X-line. In the wavenumber space it also becomes -2.8 as soon as the normalized reconnection rate reaches 0.1. The broadband kinetic turbulence is mainly due to current-streaming and electron-flow-shear instabilities excited in the sufficiently thin current sheets of kinetic reconnection. The growth of the kinetic turbulence corresponds to high reconnection rates which exceed those of fast laminar, non-turbulent reconnection., Comment: 13 pages, 10 figures. Revised version. Published in Physical Review E

We present a new hybrid algorithm implemented in the code CHIEF (Code Hybrid with Inertial Electron Fluid) for simulations of electron-ion plasmas. The algorithm treats the ions kinetically, modeled by the Particle-in-Cell (PiC) method, and electrons as an inertial fluid, modeled by electron fluid equations without any of the approximations used in most of the other hybrid codes with an inertial electron fluid. This kind of code is appropriate to model a large variety of quasineutral plasma phenomena where the electron inertia and/or ion kinetic effects are relevant. We present here the governing equations of the model, how these are discretized and implemented numerically, as well as six test problems to validate our numerical approach. Our chosen test problems, where the electron inertia and ion kinetic effects play the essential role, are: 0) Excitation of parallel eigenmodes to check numerical convergence and stability, 1) parallel (to a background magnetic field) propagating electromagnetic waves, 2) perpendicular propagating electrostatic waves (ion Bernstein modes), 3) ion beam right-hand instability (resonant and non-resonant), 4) ion Landau damping, 5) ion firehose instability, and 6) 2D oblique ion firehose instability. Our results reproduce successfully the predictions of linear and non-linear theory for all these problems, validating our code. All properties of this hybrid code make it ideal to study multi-scale phenomena between electron and ion scales such as collisionless shocks, magnetic reconnection and kinetic plasma turbulence in the dissipation range above the electron scales., Comment: 57 pages, 19 figures. Revised to match with the version published in Computer Physics Communications

Transport in collisionless plasmas is usually called anomalous, being due to the interaction between the particles and the self-generated turbulence by their collective interactions. Because of its relevance for astrophysical and space plasmas, we explore the excitation of turbulence in current sheets prone to component- or guide-field reconnection, a process not well understood, yet. We analyze the anomalous transport properties by using 2.5D Particle-in-Cell (PiC) code simulations. We split off the mean, slow variation (in contrast to the fast turbulent fluctuations) of the macroscopic observables and determine the main transport terms of the generalized Ohm's law. We verify our findings by comparing with the independently determined slowing-down rate of the macroscopic currents and with the transport terms obtained by the first order correlations of the turbulent fluctuations. We find that the turbulence is most intense in the "low density" separatrix region of guide-field reconnection. It is excited by streaming instabilities, it is mainly electrostatic, "patchy" in space and so is the associated quasi-periodic anomalous transport. The remaining irreversible anomalous resistivity can be parametrized by an effective collision frequency ranging from the local ion-cyclotron to the lower-hybrid frequency. The contributions to the parallel and the perpendicular (to the magnetic field) components of the slowly varying, DC-electric fields, balanced by the turbulence, are similar. This anomalous electric field is, however, smaller than the contributions of the off-diagonal pressure and electron inertia terms of the Ohm's law. This result can now be verified by in-situ measurements of the turbulence, in and around the magnetic reconnection regions of the Earth's magnetosphere by the multi-spacecraft mission MMS and in laboratory experiments like MRX and VINETA-II., Comment: 21 pages, 21 figures. Revised to match with the version published in Physics of Plasmas. An abridged version of the abstract is shown here

Non-Maxwellian electron velocity space distribution functions (EVDF) are useful signatures of plasma conditions and non-local consequences of collisionless magnetic reconnection. In the past, EVDFs were obtained mainly for antiparallel reconnection and under the influence of weak guide-fields in the direction perpendicular to the reconnection plane. EVDFs are, however, not well known, yet, for oblique (or component-) reconnection in dependence on stronger guide-magnetic fields and for the exhaust (outflow) region of reconnection away from the diffusion region. In view of the multi-spacecraft Magnetospheric Multiscale Mission (MMS), we derived the non-Maxwellian EVDFs of collisionless magnetic reconnection in dependence on the guide-field strength $b_g$ from small ($b_g\approx0$) to very strong ($b_g=8$) guide-fields, taking into account the feedback of the self-generated turbulence. For this sake, we carried out 2.5D fully-kinetic Particle-in-Cell simulations using the ACRONYM code. We obtained anisotropic EVDFs and electron beams propagating along the separatrices as well as in the exhaust region of reconnection. The beams are anisotropic with a higher temperature in the direction perpendicular rather than parallel to the local magnetic field. The beams propagate in the direction opposite to the background electrons and cause instabilities. We also obtained the guide-field dependence of the relative electron-beam drift speed, threshold and properties of the resulting streaming instabilities including the strongly non-linear saturation of the self-generated plasma turbulence. This turbulence and its non-linear feedback cause non-adiabatic parallel electron acceleration and EVDFs well beyond the limits of the quasi-linear approximation, producing phase space holes and an isotropizing pitch-angle scattering., Comment: 21 pages, 8 figures. Revised to match with the version published in Physics of Plasmas. An abridged version of the abstract is shown here

Aims: We investigate the electron acceleration in convective electric fields of cascading magnetic reconnection in a flaring solar corona and show the resulting hard X-ray (HXR) radiation spectra caused by Bremsstrahlung for the coronal source. Methods: We perform test particle calculation of electron motions in the framework of a guiding center approximation. The electromagnetic fields and their derivatives along electron trajectories are obtained by linearly interpolating the results of high-resolution adaptive mesh refinement (AMR) MHD simulations of cascading magnetic reconnection. Hard X-ray (HXR) spectra are calculated using an optically thin Bremsstrahlung model. Results: Magnetic gradients and curvatures in cascading reconnection current sheet accelerate electrons: trapped in magnetic islands, precipitating to the chromosphere and ejected into the interplanetary space. The final location of an electron is determined by its initial position, pitch angle and velocity. These initial conditions also influence electron acceleration efficiency. Most of electrons have enhanced perpendicular energy. Trapped electrons are considered to cause the observed bright spots along coronal mass ejection CME-trailing current sheets as well as the flare loop-top HXR emissions., Comment: submitted to A&A

In this work, we compare gyrokinetic (GK) and fully kinetic Particle-in-Cell (PIC) simulations of magnetic reconnection in the limit of strong guide field. In particular, we analyze the limits of applicability of the GK plasma model compared to a fully kinetic description of force free current sheets for finite guide fields ($b_g$). Here we report the first part of an extended comparison, focusing on the macroscopic effects of the electron flows. For a low beta plasma ($\beta_i=0.01$), it is shown that both plasma models develop magnetic reconnection with similar features in the secondary magnetic islands if a sufficiently high guide field ($b_g\gtrsim 30$) is imposed in the kinetic PIC simulations. Outside of these regions, in the separatrices close to the X points, the convergence between both plasma descriptions is less restrictive ($b_g\gtrsim 5$). Kinetic PIC simulations using guide fields $b_g \lesssim 30$ reveal secondary magnetic islands with a core magnetic field and less energetic flows inside of them in comparison to the GK or kinetic PIC runs with stronger guide fields. We find that these processes are mostly due to an initial shear flow absent in the GK initialization and negligible in the kinetic PIC high guide field regime, in addition to fast outflows on the order of the ion thermal speed that violate the GK ordering. Since secondary magnetic islands appear after the reconnection peak time, a kinetic PIC/GK comparison is more accurate in the linear phase of magnetic reconnection. For a high beta plasma ($\beta_i=1.0$) where reconnection rates and fluctuations levels are reduced, similar processes happen in the secondary magnetic islands in the fully kinetic description, but requiring much lower guide fields ($b_g\lesssim 3$)., Comment: 18 pages, 13 figures. Revised to match with the published version in Physics of Plasmas

In this work, we investigate the influence of the anisotropic heating on the spontaneous instability and evolution of thin Harris-type collisionless current sheets, embedded in antiparallel magnetic fields. In particular, we explore the influence of the macroparticle shape-function using a 2D version of the PIC code ACRONYM. We also investigate the role of the numerical collisionality due to the finite number of macroparticles in PIC codes. It is shown that it is appropriate to choose higher order shape functions of the macroparticles compared to a larger number of macroparticles per cell. This allows to estimate better the anisotropic electron heating due to the collisions of macroparticles in a PIC code. Temperature anisotropies can stabilize the tearing mode instability and trigger additional current sheet instabilities. We found a good agreement between the analytically derived threshold for the stabilization of the anisotropic tearing mode and other instabilities, either spontaneously developing or initially triggered ones. Numerical effects causing anisotropic heating at electron time scales, become especially important for higher mass ratios (above $m_i/m_e = 180$). If numerical effects are carefully taken into account, one can recover the theoretical estimated linear growth rates of the tearing instability of thin isotropic collisionless current sheets, also for higher mass ratios., Comment: 26 pages, 30 figures

The numerical simulation of turbulence and flows in almost ideal, large-Reynolds-number astrophysical plasmas motivates the implementation of almost conservative MHD computer codes. They should efficiently calculate, use highly parallelized schemes scaling well with large numbers of CPU cores, allows to obtain a high grid resolution over large simulation domains and which can easily be adapted to new computer architectures as well as to new initial and boundary conditions, allow modular extensions. The new massively parallel simulation code GOEMHD3 enables efficient and fast simulations of almost ideal, large-Reynolds-number astrophysical plasma flows, well resolved and on huge grids covering large domains. Its abilities are validated by major tests of ideal and weakly dissipative plasma phenomena. The high resolution ($2048^3$ grid points) simulation of a large part of the solar corona above an observed active region proved the excellent parallel scalability of the code using more than 30.000 processor cores., Comment: The revised version

Context. The solar corona is heated to high temperatures of the order of 10^{6} K. The coronal energy budget and specifically possible mechanisms of coronal heating (wave, DC-electric fields, ..) are poorly understood. This is particularly true as far as the formation of X-ray bright points (BPs) is concerned. Aims. Investigation of the energy budget with emphasis on the relative role and contribution of adiabatic compression versus current dissipation to the formation of coronal BPs. Methods. Three-dimensional resistive MHD simulation starts with the extrapolation of the observed magnetic field from SOHO/MDI magnetograms, which are associated with a BP observed on 19 December 2006 by Hinode. The initial radially non-uniform plasma density and temperature distribution is in accordance with an equilibrium model of chromosphere and corona. The plasma motion is included in the model as a source of energy for coronal heating. Results. Investigation of the energy conversion due to Lorentz force, pressure gradient force and Ohmic current dissipation for this bright point shows the minor effect of Joule heating in comparison to the work done by pressure gradient force in increasing the thermal energy by adiabatic compression. Especially at the time when the temperature enhancement above the bright point starts to form, compressional effects are quite dominant over the direct Joule heating. Conclusions. Choosing non-realistic high resistivity in compressible MHD models for simulation of solar corona can lead to unphysical consequences for the energy balance analysis, especially when local thermal energy enhancements are being considered., Comment: 10 pages, 12 figures

We present a newly developed approach to solar coronal magnetic field extrapolation from vector magnetograms, based on the Principle of Minimum Dissipation Rate (MDR). The MDR system was derived from a variational problem that is more suitable for an open and externally driven system, like the solar corona. The resulting magnetic field equation is more general than force-free. Its solution can be expressed as the superposition of two linear (constant-$\alpha$) force-free fields (LFFFs) with distinct $\alpha$ parameters, and one potential field. Thus the original extrapolation problem is decomposed into three LFFF extrapolations, utilizing boundary data. The full MDR-based approach requires two layers of vector magnetograph measurements on solar surface, while a slightly modified practical approach only requires one. We test both approaches against 3D MHD simulation data in a finite volume. Both yield quantitatively good results. The errors in the magnetic energy estimate are within a few percents. In particular, the main features of relatively strong perpendicular current density structures, representative of the non-force freeness of the solution, are well recovered., Comment: submitted to ApJ