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

We present a new boundary condition for simulations of magnetic reconnection based on the topology of a Klein bottle. When applicable, the new condition is computationally cheaper than fully periodic boundary conditions, reconnects more flux than systems with conducting boundaries, and does not require assumptions about regions external to the simulation as is necessary for open boundaries. The new condition reproduces the expected features of reconnection, but cannot be straightforwardly applied in systems with asymmetric upstream plasmas., Comment: 7 pages, 8 figures, submitted to Physics of Plasmas

This short article highlights the unsolved problems of magnetic reconnection in collisionless plasma. The advanced in-situ plasma measurements and simulations enabled scientists to gain a novel understanding of magnetic reconnection. Still, outstanding questions remain on the complex dynamics and structures in the diffusion region, on the cross-scale and regional couplings, on the onset of magnetic reconnection, and on the details of energetics. Future directions of the magnetic reconnection research in terms of new observations, new simulations and interdisciplinary approaches are discussed., Comment: Submitted to Space Science Reviews. This is a review paper and is an outcome of the International Space Science Institute (ISSI) Workshop on Magnetic Reconnection: Explosive Energy Conversion in Space Plasmas held at June 27 - July 1, 2022

Switchbacks (SBs) are localized structures in the solar wind containing deflections of the magnetic field direction relative to the background solar wind magnetic field. The amplitudes of the magnetic field deflection angles for different SBs vary from ~40 to ~160-170 degrees. Alignment of the perturbations of the magnetic field and the bulk solar wind velocity is observed inside SBs and causes spiky enhancements of the radial bulk velocity inside SBs. We have investigated the deviations of SB perturbations from Alfv\'enicity by evaluating the distribution of the parameter defined as the ratio of the parallel to $\Delta\vec{B}$ component of $\Delta\vec{V}$ to $\Delta\vec{V}_A = \Delta\vec{B}/4\pi n_i m_i$ inside SBs, i.e. $\alpha=V_{||}/|\Delta\vec{V}_A|$ ( $\alpha=\Delta\vec{V}/|\Delta\vec{V}_A|$ when $\Delta\vec{V}~\Delta\vec{B}$), which quantifies the deviation of the perturbation from an Alfv\'enic one. Based on Parker Solar Probe (PSP) observations, we show that $\alpha$ inside SBs has systematically lower values than it has in the pristine solar wind: $\alpha$ inside SBs observed during PSP Encounter 1 were distributed in a range of 0.2-0.9. The upper limit on $\alpha$ is constrained by the requirement that the jump in velocity across the switchback boundary be less than the local Alfv\'en speed. This prevents the onset of shear flow instabilities. The consequence is that the perturbation of the proton bulk velocity in SBs with deflection greater than 60 degrees cannot reach $\alpha=1$ (the Alfv\'enicity condition) and the highest possible $\alpha$ for a SB with the full reversal of B is 0.5. These results have consequences for the interpretation of switchbacks as large amplitude Alfv\'en waves.

The waves generated by high-energy proton and alpha particles streaming from solar flares into regions of colder plasma are explored using particle-in-cell simulations. Initial distribution functions for the protons and alphas consist of two populations: an energetic, streaming population represented by an anisotropic ($T_{\parallel} > T_{\perp}$), one-sided kappa function and a cold, Maxwellian background population. The anisotropies and non-zero heat fluxes of these distributions destabilize oblique waves with a range of frequencies below the proton cyclotron frequency. These waves scatter particles out of the tails of the initial distributions along constant energy surfaces in the wave frame. Overlap of the nonlinear resonance widths allows particles to scatter into near isotropic distributions by the end of the simulations. The dynamics of $^3$He are explored using test particles. Their temperatures can increase by a factor of nearly 20. Propagation of such waves into regions above and below the flare site can lead to heating and transport of $^3$He into the flare acceleration region. The amount of heated $^3$He that will be driven into the flare site is proportional to the wave energy. Using values from our simulations, we show that the abundance of $^3$He driven into the acceleration region should approach that of $^4$He in the corona. Therefore, waves driven by energetic ions produced in flares are a strong candidate to drive the enhancements of $^3$He observed in impulsive flares., Comment: 19 pages, 12 figures

There is ample evidence for magnetic reconnection in the solar system, but it is a nontrivial task to visualize, to determine the proper approaches and frames to study, and in turn to elucidate the physical processes at work in reconnection regions from in-situ measurements of plasma particles and electromagnetic fields. Here an overview is given of a variety of single- and multi-spacecraft data analysis techniques that are key to revealing the context of in-situ observations of magnetic reconnection in space and for detecting and analyzing the diffusion regions where ions and/or electrons are demagnetized. We focus on recent advances in the era of the Magnetospheric Multiscale mission, which has made electron-scale, multi-point measurements of magnetic reconnection in and around Earth's magnetosphere., Comment: submitted to Space Science Reviews (116 pages, incl. 31 figures, 7 tables)

Using observations from the Solar Dynamics Observatory's Atmosphere Imaging Assembly and the Ramaty High Energy Solar Spectroscopic Imager, we present novel measurements of the shear of post-reconnection flare loops (PRFLs) in SOL20141218T21:40 and study its evolution with respect to magnetic reconnection and flare emission. Two quasi-parallel ribbons form adjacent to the magnetic polarity inversion line (PIL), spreading in time first parallel to the PIL and then mostly in a perpendicular direction. We measure magnetic reconnection rate from the ribbon evolution, and also the shear angle of a large number of PRFLs observed in extreme ultraviolet passbands ($\lesssim$1 MK). For the first time, the shear angle measurements are conducted using several complementary techniques allowing for a cross-validation of the results. In this flare, the total reconnection rate is much enhanced before a sharp increase of the hard X-ray emission, and the median shear decreases from 60$^\circ$-70$^\circ$ to 20$^\circ$, on a time scale of ten minutes. We find a correlation between the shear-modulated total reconnection rate and the non-thermal electron flux. These results confirm the strong-to-weak shear evolution suggested in previous observational studies and reproduced in numerical models, and also confirm that, in this flare, reconnection is not an efficient producer of energetic non-thermal electrons during the first ten minutes when the strongly sheared PRFLs are formed. We conclude that an intermediate shear angle, $\le 40^\circ$, is needed for efficient particle acceleration via reconnection, and we propose a theoretical interpretation., Comment: 19 pages, 10 figures

The role of interchange reconnection as a drive mechanism for the solar wind is explored by solving the global magnetic-field-aligned equations describing wind acceleration. Boundary conditions in the low corona, including a reconnection-driven Alfv\'enic outflow and associated heating differ from previous models. Additional heating of the corona associated with Alfv\'en waves or other MHD turbulence, which has been the foundation of many earlier models, is neglected. For this simplified model a sufficient condition for interchange reconnection to overcome gravity to drive the wind is derived. The combination of Alfv\'enic ejection and reconnection-driven heating yields a minimum value of the Alfv\'en speed of the order of 350-400$km/s$ that is required to drive the wind. Recent evidence based on Parker Solar Probe (PSP) observations suggests that this threshold is typically exceeded in the coronal holes that are the source regions of the fast wind. On the other hand, since reconnection in the coronal environment is predicted to have a bursty character, the magnitude of reconnection outflows can be highly variable. The consequence is a highly non-uniform wind in which in some regions the velocity increases sharply to super-Alfv\'enic values while in adjacent regions the formation of an asymptotic wind fails. A simple model is constructed to describe the turbulent mixing of these highly-sheared super-Alfv\'enic flows that suggests these flows are the free-energy source of the Alfv\'enic turbulence and associated switchbacks that have been documented in the PSP data in the near coronal environment. The global wind profiles are presented and benchmarked with Parker Solar Probe (PSP) observations at 12 solar radii.

The fast solar wind that fills the heliosphere originates from deep within regions of open magnetic field on the Sun called coronal holes. The energy source responsible for accelerating the plasma to high speeds is widely debated, however there is evidence that it is ultimately magnetic in nature with candidate mechanisms including wave heating^(1,2) and interchange reconnection^(3,4,5). The coronal magnetic field near the solar surface is structured on scales associated with supergranulation convection cells, where descending flows create intense fields. The energy density in these network magnetic field bundles is a likely candidate as an energy source of the wind. Here we report measurements of fast solar wind streams from the Parker Solar Probe (PSP) spacecraft^6 which provides strong evidence for the interchange reconnection mechanism. We show that supergranulation structure at the coronal hole base remains imprinted in the near-Sun solar wind resulting in asymmetric patches of magnetic 'switchbacks'^(7,8) and bursty wind streams with power law-like energetic ion spectra to beyond 100 keV. Computer simulations of interchange reconnection support key features of the observations, including the ion spectra. Important characteristics of interchange reconnection in the low corona are inferred from the data including that the reconnection is collisionless and that the energy release rate is sufficient to power the fast wind. In this scenario, open magnetic flux undergoes continuous reconnection and the wind is driven both by the resulting plasma pressure and the radial Alfvenic flow bursts., Comment: 15 pages, 4 figures. Nature, 2023

A major discovery of Parker Solar Probe (PSP) was the presence of large numbers of localized increases in the radial solar wind speed and associated sharp deflections of the magnetic field - switchbacks (SB). A possible generation mechanism of SBs is through magnetic reconnection between open and closed magnetic flux near the solar surface, termed interchange reconnection that leads to the ejection of flux ropes (FR) into the solar wind. Observations also suggest that SBs undergo merging, consistent with a FR picture of these structures. The role of FRs merging in controlling the structure of SB in the solar wind is explored through direct observations, through analytic analysis, and numerical simulations. Analytic analysis reveals key features of the structure of FR and their scaling with the heliocentric distance R that are consistent with observations and that reveal the critical role of merging in controlling the SB structure. FR merging is shown to be energetically favorable to reduce the strength of the wrapping magnetic field and drive the observed elongation of SBs. A further consequence is the resulting dominance of the axial magnetic field within SBs that leads to the characteristic sharp rotation of the magnetic field into the axial direction at the SB boundary that is revealed in observations. Finally, the radial scaling of the SB area in the FR model of SBs suggests that the observational probability of SB identification should be insensitive to R, which is consistent with the most recent statistical analysis of SB observations from PSP., Comment: 21 pages, 9 figures

The reversibility of the transfer of energy from the magnetic field to the surrounding plasma during magnetic reconnection is examined. Trajectories of test particles in an analytic model of the fields demonstrate that irreversibility is associated with separatrix crossings and regions of weaker magnetic field. Inclusion of a guide field increases the degree of reversibility. Full kinetic simulations with a particle-in-cell code support these results and demonstrate that while time-reversed simulations at first "un-reconnect", they eventually evolve into a reconnecting state.

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

The structure of magnetic flux ropes injected into the solar wind during reconnection in the coronal atmosphere is explored with particle-in-cell simulations and compared with in situ measurements of magnetic "switchbacks" from the Parker Solar Probe. We suggest that multi-x-line reconnection between open and closed flux in the corona injects flux ropes into the solar wind and that these flux ropes convect outward over long distances before eroding due to reconnection. Simulations that explore the magnetic structure of flux ropes in the solar wind reproduce the following key features of the switchback observations: a rapid rotation of the radial magnetic field into the transverse direction, which is a consequence of reconnection with a strong guide field; and the potential to reverse the radial field component. The potential implication of the injection of large numbers of flux ropes in the coronal atmosphere for understanding the generation of the solar wind is discussed.

Transport equations for electron thermal energy in the high $\beta_e$ intracluster medium (ICM) are developed that include scattering from both classical collisions and self-generated whistler waves. The calculation employs an expansion of the kinetic electron equation along the ambient magnetic field in the limit of strong scattering and assumes whistler waves with low phase speeds $V_w\sim{v}_{te}/\beta_e\ll{v}_{te}$ dominate the turbulent spectrum, with $v_{te}$ the electron thermal speed and $\beta_e\gg1$ the ratio of electron thermal to magnetic pressure. We find: (1) temperature-gradient-driven whistlers dominate classical scattering when $L_c>L/\beta_e$, with $L_c$ the classical electron mean-free-path and $L$ the electron temperature scale length, and (2) in the whistler dominated regime the electron thermal flux is controlled by both advection at $V_w$ and a comparable diffusive term. The findings suggest whistlers limit electron heat flux over large regions of the ICM, including locations unstable to isobaric condensation. Consequences include: (1) the Field length decreases, extending the domain of thermal instability to smaller length-scales, (2) the heat flux temperature dependence changes from $T_e^{7/2}/L$ to $V_wnT_e\sim{T}_e^{1/2}$, (3) the magneto-thermal and heat-flux driven buoyancy instabilities are impaired or completely inhibited, and (4) sound waves in the ICM propagate greater distances, as inferred from observations. This description of thermal transport can be used in macroscale ICM models., Comment: Accepted for publication in The Astrophysical Journal

Monte Carlo methods are often employed to numerically integrate kinetic equations, such as the particle-in-cell method for the plasma kinetic equation, but these methods suffer from the introduction of counting noise to the solution. We report on a cautionary tale of counting noise modifying the nonlinear saturation of kinetic instabilities driven by unstable beams of plasma. We find a saturated magnetic field in under-resolved particle-in-cell simulations due to the sampling error in the current density. The noise-induced magnetic field is anomalous, as the magnetic field damps away in continuum kinetic and increased particle count particle-in-cell simulations. This modification of the saturated state has implications for a broad array of astrophysical phenomena beyond the simple plasma system considered here, and it stresses the care that must be taken when using particle methods for kinetic equations., Comment: 15 pages, 4 figures

We analyze the development and influence of turbulence in three-dimensional particle-in-cell simulations of guide-field magnetic reconnection at the magnetopause with parameters based on observations of an electron diffusion region by the Magnetospheric Multiscale (MMS) mission. Along the separatrices the turbulence is a variant of the lower hybrid drift instability (LHDI) that produces electric field fluctuations with amplitudes much greater than the reconnection electric field. The turbulence controls the scale length of the density and current profiles while enabling significant transport across the magnetopause despite the electrons remaining frozen-in to the magnetic field. Near the X-line the electrons are not frozen-in and the turbulence, which differs from the LHDI, makes a significant net contribution to the generalized Ohm's law through an anomalous viscosity. The characteristics of the turbulence and associated particle transport are consistent with fluctuation amplitudes in the MMS observations. However, for this event the simulations suggest that the MMS spacecraft were not close enough to the core of the electron diffusion region to identify the region where anomalous viscosity is important.

The scattering of electrons by heat-flux-driven whistler waves is explored with a particle-in-cell (PIC) simulation relevant to the transport of energetic electrons in flares. The simulation is initiated with a large heat flux that is produced using a kappa distribution of electrons with positive velocity and a cold return current beam. This system represents energetic electrons escaping from a reconnection-driven energy release site. This heat flux system drives large amplitude oblique whistler waves propagating both along and against the heat flux, as well as electron acoustic waves. While the waves are dominantly driven by the low energy electrons, including the cold return current beam, the energetic electrons resonate with and are scattered by the whistlers on time scales of the order of a hundred electron cyclotron times. Peak whistler amplitudes of $\tilde{B} / B_{0} \sim 0.125$ and angles of $\sim 60 \degree$ with respect to the background magnetic field are observed. Electron perpendicular energy is increased while the field-aligned electron heat flux is suppressed. The resulting scattering mean-free-paths of energetic electrons are small compared with the typical scale size of energy release sites in flares, which might lead to the effective confinement of energetic electrons that is required for the production of very energetic particles., Comment: Accepted for publication in The Astrophysical Journal

The role of turbulence in low-$\beta$, guide-field reconnection exhausts is explored in 2D reconnection and 2D and 3D Riemann simulations. The structure of the exhaust and associated turbulence is controlled by a pair of rotational discontinuities (RDs) at the exhaust boundary and a pair of slow shocks (SSs) that are generated by counterstreaming ions beams. In 2D the exhaust develops large-amplitude striations at the ion Larmor radius scale that are produced by electron-beam-driven ion cyclotron waves. The electron beams driving the instability are injected into the exhaust from one of the RDs. However, in 3D Riemann simulations, the additional dimension (in the out-of-plane direction) results in strong Buneman and electron-electron streaming instabilities at the RD which suppress electron beam formation and therefore the striations in the exhaust. The strength of the streaming instabilities at the RD are controlled by the ratio of the electron thermal speed to Alfv\'en speed, lower thermal speed being more unstable. In the 3D simulations an ion-ion streaming instability acts to partially thermalize the counterstreaming ion beams at the SSs. This instability is controlled by the ratio of the sound speed to Alfv\'en speed and is expected to be stable in the low $\beta$ solar corona. The results suggest that in a guide field reconnection exhaust with $1 \gg \beta > m_e/m_i$, the kinetic-scale turbulence that develops will be too weak to play a significant role in energy conversion and particle acceleration. Therefore, the energy conversion will be mostly controlled by laminar physics or multi-x-line reconnection.

Kinetic Riemann simulations have been completed to explore particle heating during guide field reconnection in the low-$\beta$ environment of the inner heliosphere and the solar corona. The reconnection exhaust is bounded by two rotational discontinuities (RD) and two slow shocks (SS) form within the exhaust as in magnetohydrodynamic (MHD) models. At the RDs, ions are accelerated by the magnetic field tension to drive the reconnection outflow as well as flows in the out-of-plane direction. The out-of-plane flows stream toward the midplane and meet to drive the SSs. The SSs differ greatly from those in the MHD model. The turbulence at the shock fronts and both upstream and downstream is weak so the shocks are laminar and produce little dissipation. Downstream of the SSs the counterstreaming ion beams lead to higher density, which leads to a positive potential between the SSs that acts to confine the downstream electrons to maintain charge neutrality. The potential accelerates electrons from upstream of the SSs to downstream region and traps a small fraction but only modestly increases the downstream electron temperature above the upstream value. In the low-$\beta$ limit the released magnetic energy is split between bulk flow and ion heating with little energy going to electrons. That the model does not produce strong electron heating nor an energetic electron component suggests that other mechanisms, such as multiple x-line reconnection, are required to explain energetic electron production in large flares. The model can be tested with the expected data from the Parker Solar Probe., Comment: Some figures' quality is downgraded due to the size limit of arXiv

Using 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 Larmor radius, the ions do not respond to the reconnection dynamics leading to ''electron-only'' reconnection with very large quasi-steady reconnection rates. 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 Alfv\'en 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.

A new computational model is presented suitable for exploring the self-consistent production of energetic electrons during magnetic reconnection in macroscale systems. The equations are based on the recent discovery that parallel electric fields are ineffective drivers of energetic particles during reconnection so that the kinetic scales which control the development of such fields can be ordered out of the equations. The resulting equations consist of a magnetohydrodynamic (MHD) backbone with the energetic component represented by macro-particles described by the guiding center equations. Crucially, the energetic component feeds back on the MHD equations so that the total energy of the MHD fluid and the energetic particles is conserved. The equations correctly describe the firehose instability, whose dynamics plays a key role in throttling reconnection and in controlling the spectra of energetic particles. The results of early tests of the model, including the propagation of Alfven waves in a system with pressure anisotropy and the growth of firehose modes, establish that the basic algorithm is stable and produces reliable physics results in preparation for further benchmarking with particle-in-cell models of reconnection.

Heat flux suppression in collisionless plasmas for a large range of plasma $\beta$ is explored using two-dimensional particle-in-cell simulations with a strong, sustained thermal gradient. We find that a transition takes place between whistler-dominated (high-$\beta$) and double-layer-dominated (low-$\beta$) heat flux suppression. Whistlers saturate at small amplitude in the low beta limit and are unable to effectively suppress the heat flux. Electrostatic double layers suppress the heat flux to a mostly constant factor of the free streaming value once this transition happens. The double layer physics is an example of ion-electron coupling and occurs on a scale of roughly the electron Debye length. The scaling of ion heating associated with the various heat flux driven instabilities is explored over the full range of $\beta$ explored. The range of plasma-$\beta$s studied in this work makes it relevant to the dynamics of a large variety of astrophysical plasmas, including the intracluster medium of galaxy clusters, hot accretion flows, stellar and accretion disk coronae, and the solar wind., Comment: 16 pages, 6 figures. Accepted for publication in Astrophys. J

Data from the NASA Magnetospheric Multiscale (MMS) mission are used to investigate asymmetric magnetic reconnection at the dayside boundary between the Earth's magnetosphere and the solar wind (the magnetopause). High-resolution measurements of plasmas, electric and magnetic fields, and waves are used to identify highly localized (~15 electron Debye lengths) standing wave structures with large electric-field amplitudes (up to 100 mV/m). These wave structures are associated with spatially oscillatory dissipation, which appears as alternatingly positive and negative values of J dot E (dissipation). For small guide magnetic fields the wave structures occur in the electron stagnation region at the magnetosphere edge of the EDR. For larger guide fields the structures also occur near the reconnection x-line. This difference is explained in terms of channels for the out-of-plane current (agyrotropic electrons at the stagnation point and guide-field-aligned electrons at the x-line).

We analyze a high-resolution simulation of magnetopause reconnection observed by the Magnetospheric Multiscale (MMS) mission and explain the occurrence of strongly localized dissipation with an amplitude more than an order of magnitude larger than expected. Unlike symmetric reconnection, wherein reconnection of the ambient reversed magnetic field drives the dissipation, we find the annihilation of the self-generated, out-of-plane (Hall) magnetic field plays the dominant role. Electrons flow along the magnetosheath separatrices, converge in the diffusion region, and jet past the X-point into the magnetosphere. The resulting accumulation of negative charge generates intense parallel electric fields that eject electrons along the magnetospheric separatrices and produce field-aligned beams. Many of these features match MMS observations., Comment: Accepted version

The dynamics of weakly magnetized collisionless plasmas in the presence of an imposed temperature gradient along an ambient magnetic field is explored with particle-in-cell simulations and modeling. Two thermal reservoirs at different temperatures drive an electron heat flux that destabilizes off-angle whistler-type modes. The whistlers grow to large amplitude, $\delta B / B_{0} \simeq 1$, and resonantly scatter the electrons, significantly reducing the heat flux. A surprise is that the resulting steady state heat flux is largely independent of the thermal gradient. The rate of thermal conduction is instead controlled by the finite propagation speed of the whistlers, which act as mobile scattering centers that convect the thermal energy of the hot reservoir. The results are relevant to thermal transport in high $\beta$ astrophysical plasmas such as hot accretion flows and the intracluster medium of galaxy clusters., Comment: Accepted for publication in Physical Review Letters. 5 pages, 5 figures. Supplementary material: 1 page, 1 figure

Magnetic reconnection is an important driver of energetic particles in many astrophysical phenomena. Using kinetic particle-in-cell (PIC) simulations, we explore the impact of three-dimensional reconnection dynamics on the efficiency of particle acceleration. In two-dimensional systems, Alfv\'enic outflows expel energetic electrons into flux ropes where they become trapped and disconnected from acceleration regions. However, in three-dimensional systems these flux ropes develop axial structure that enables particles to leak out and return to acceleration regions. This requires a finite guide field so that particles may move quickly along the flux rope axis. We show that greatest energetic electron production occurs when the guide field is of the same order as the reconnecting component: large enough to facilitate strong transport, but not so large as to throttle the dominant Fermi mechanism responsible for efficient electron acceleration. This suggests a natural explanation for the envelope of electron acceleration during the impulsive phase of eruptive flares.

We present detailed analysis of the turbulence observed in three-dimensional particle-in-cell simulations of magnetic reconnection at the magnetopause. The parameters are representative of an electron diffusion region encounter of the Magnetospheric Multiscale (MMS) mission. The turbulence is found to develop around both the magnetic X line and separatrices, is electromagnetic in nature, is characterized by a wave vector $k$ given by $k\rho_e\sim(m_eT_e/m_iT_i)^{0.25}$ with $\rho_e$ the electron Larmor radius, and appears to have the ion pressure gradient as its source of free energy. Taken together, these results suggest the instability is a variant of the lower hybrid drift instability. The turbulence produces electric field fluctuations in the out-of-plane direction (the direction of the reconnection electric field) with an amplitude of around $\pm 10$~mV/m, which is much greater than the reconnection electric field of around $0.1$~mV/m. Such large values of the out-of-plane electric field have been identified in the MMS data. The turbulence in the simulations controls the scale lengths of the density profile and current layers in asymmetric reconnection, driving them closer to $\sqrt{\rho_e\rho_i}$ than the $\rho_e$ or $d_e$ scalings seen in 2-D reconnection simulations, and produces significant anomalous resistivity and viscosity in the electron diffusion region., Comment: 11 pages, 10 figures

Based on the difference between the orientation of the interstellar $B_{ISM}$ and the solar magnetic fields, there was an expectation that the magnetic field direction would rotate dramatically across the heliopause (HP). However, the Voyager 1 spacecraft measured very little rotation across the HP. Previously we showed that the $B_{ISM}$ twists as it approaches the HP and acquires a strong T component (East-West). Here we establish that reconnection in the eastern flank of the heliosphere is responsible for the twist. On the eastern flank the solar magnetic field has twisted into the positive N direction and reconnects with the Southward pointing component of the $B_{ISM}$. Reconnection drives a rotational discontinuity (RD) that twists the $B_{ISM}$ into the -T direction and propagates upstream in the interstellar medium towards the nose. The consequence is that the N component of $B_{ISM}$ is reduced in a finite width band upstream of the HP. Voyager 1 currently measures angles ($\delta=sin^{-1}(B_{N}/B)$) close to solar values. We present MHD simulations to support this scenario, suppressing reconnection in the nose region while allowing it in the flanks, consistent with recent ideas about reconnection suppression from diamagnetic drifts. The jump in plasma $\beta$ (the plasma to magnetic pressure) across the nose of HP is much greater than in the flanks because the heliosheath $\beta$ is greater there than in the flanks. Large-scale reconnection is therefore suppressed in the nose but not at the flanks. Simulation data suggest that $B_{ISM}$ will return to its pristine value $10-15~AU$ past the HP., Comment: 19 pages, 5 figures, submitted

The misalignment of the solar rotation axis and the magnetic axis of the Sun produces a periodic reversal of the Parker spiral magnetic field and the sectored solar wind. The compression of the sectors is expected to lead to reconnection in the heliosheath (HS). We present particle-in-cell simulations of the sectored HS that reflect the plasma environment along the Voyager 1 and 2 trajectories, specifically including unequal positive and negative azimuthal magnetic flux as seen in the Voyager data \citep{Burlaga03}. Reconnection proceeds on individual current sheets until islands on adjacent current layers merge. At late time bands of the dominant flux survive, separated by bands of deep magnetic field depletion. The ambient plasma pressure supports the strong magnetic pressure variation so that pressure is anti-correlated with magnetic field strength. There is little variation in the magnetic field direction across the boundaries of the magnetic depressions. At irregular intervals within the magnetic depressions are long-lived pairs of magnetic islands where the magnetic field direction reverses so that spacecraft data would reveal sharp magnetic field depressions with only occasional crossings with jumps in magnetic field direction. This is typical of the magnetic field data from the Voyager spacecraft \citep{Burlaga11,Burlaga16}. Voyager 2 data reveals that fluctuations in the density and magnetic field strength are anti-correlated in the sector zone as expected from reconnection but not in unipolar regions. The consequence of the annihilation of subdominant flux is a sharp reduction in the "number of sectors" and a loss in magnetic flux as documented from the Voyager 1 magnetic field and flow data \citep{Richardson13}.

We investigate magnetic reconnection in systems simultaneously containing asymmetric (anti-parallel) magnetic fields, asymmetric plasma densities and temperatures, and arbitrary in-plane bulk flow of plasma in the upstream regions. Such configurations are common in the high-latitudes of Earth's magnetopause and in tokamaks. We investigate the convection speed of the X-line, the scaling of the reconnection rate, and the condition for which the flow suppresses reconnection as a function of upstream flow speeds. We use two-dimensional particle-in-cell simulations to capture the mixing of plasma in the outflow regions better than is possible in fluid modeling. We perform simulations with asymmetric magnetic fields, simulations with asymmetric densities, and simulations with magnetopause-like parameters where both are asymmetric. For flow speeds below the predicted cutoff velocity, we find good scaling agreement with the theory presented in Doss et al., J.~Geophys.~Res., 120, 7748 (2015). Applications to planetary magnetospheres, tokamaks, and the solar wind are discussed., Comment: 17 pages, 4 figures, submitted to Physics of Plasmas

We present two-dimensional kinetic simulations, with a broad range of initial guide fields, that isolate the role of parallel electric fields ($E_\parallel$) in energetic electron production during collisionless magnetic reconnection. In the strong guide field regime, $E_\parallel$ drives essentially all of the electron energy gain, yet fails to generate an energetic component. We suggest that this is due to the weak energy scaling of particle acceleration from $E_\parallel$ compared to that of a Fermi-type mechanism responsible for energetic electron production in the weak guide-field regime. This result has important implications for energetic electron production in astrophysical systems and reconnection-driven dissipation in turbulence.

Understanding the thermodynamic state of the hot intracluster medium (ICM) in a galaxy cluster requires a knowledge of the plasma transport processes, especially thermal conduction. The basic physics of thermal conduction in plasmas with ICM-like conditions has yet to be elucidated, however. We use particle-in-cell simulations and analytic models to explore the dynamics of an ICM-like plasma (with small gyroradius, large mean-free-path, and strongly sub-dominant magnetic pressure) driven by the diffusive heat flux associated with thermal conduction. Lin- ear theory reveals that whistler waves are driven unstable electron heat flux, even when the heat flux is weak. The resonant interaction of electrons with these waves then plays a critical role in scattering electrons and suppressing the heat flux. In a 1D model where only whistler modes that are parallel to the magnetic field are captured, the only resonant electrons are moving in the opposite direction to the heat flux and the electron heat flux suppression is small. In 2D or more, oblique whistler modes also resonate with electrons moving in the direction of the heat flux. The overlap of resonances leads to effective symmetrization of the electron distribution function and a strong suppression of heat flux. The results suggest that thermal conduction in the ICM might be strongly suppressed, possibly to negligible levels., Comment: 8 pages, 4 figures. To be published in ApJ Letters

Two- and three-dimensional particle-in-cell simulations of a recent encounter of the Magnetospheric Multiscale Mission (MMS) with an electron diffusion region at the magnetopause are presented. While the two-dimensional simulation is laminar, turbulence develops at both the x-line and along the magnetic separatrices in the three-dimensional simulation. The turbulence is strong enough to make the magnetic field around the reconnection island chaotic and produces both anomalous resistivity and anomalous viscosity. Each contribute significantly to breaking the frozen-in condition in the electron diffusion region. A surprise is that the crescent-shaped features in velocity space seen both in MMS observations and in two-dimensional simulations survive, even in the turbulent environment of the three-dimensional system. This suggests that MMS's measurements of crescent distributions do not exclude the possibility that turbulence plays an important role in magnetopause reconnection., Comment: Revised version accepted by GRL

Kinetic particle-in-cell simulations are used to identify signatures of the electron diffusion region (EDR) and its surroundings during asymmetric magnetic reconnection. A "shoulder" in the sunward pointing normal electric field (EN > 0) at the reconnection magnetic field reversal is a good indicator of the EDR, and is caused by magnetosheath electron meandering orbits in the vicinity of the x-line. Earthward of the X-line, electrons accelerated by EN form strong currents and crescent-shaped distribution functions in the plane perpendicular to B. Just downstream of the X-line, parallel electric fields create field-aligned crescent electron distribution functions. In the immediate upstream magnetosheath, magnetic field strength, plasma density, and perpendicular electron temperatures are lower than the asymptotic state. In the magnetosphere inflow region, magnetosheath ions intrude resulting in an Earthward pointing electric field and parallel heating of magnetospheric particles. Many of the above properties persist with a guide field of at least unity., Comment: Submitted to Geophysical Research Letters

A new scalar measure of the gyrotropy of a pressure tensor is defined. Previously suggested measures are shown to be incomplete by means of examples for which they give unphysical results. To demonstrate its usefulness as an indicator of magnetic topology, the new measure is calculated for electron data taken from numerical simulations of magnetic reconnection, shown to peak at separatrices and X-points, and compared to the other measures. The new diagnostic has potential uses in analyzing spacecraft observations and so a method for calculating it from measurements performed in an arbitrary coordinate system is derived., Comment: Under review. Comments welcome; Now accepted. New version incorporates reviewer comments

An analytic model of the heliosheath (HS) between the termination shock (TS) and the heliopause (HP) is developed in the limit in which the interstellar flow and magnetic field are neglected. The heliosphere in this limit is axisymmetric and the overall structure of the HS and HP are controlled by the solar magnetic field even in the limit in which the ratio of the plasma to magnetic field pressure, $\beta=8\pi P/B^2$, in the HS is large. The tension of the solar magnetic field produces a drop in the total pressure between the TS and the HP. This same pressure drop accelerates the plasma flow downstream of the TS into the North and South directions to form two collimated jets. The radii of these jets are controlled by the flow through the TS and the acceleration of this flow by the magnetic field -- a stronger solar magnetic field boosts the velocity of the jets and reduces the radii of the jets and the HP. Magnetohydrodynamic (MHD) simulations of the global helioshere embedded in a stationary interstellar medium match well with the analytic model. The results suggest that mechanisms that reduce the HS plasma pressure downstream of the TS can enhance the jet outflow velocity and reduce the HP radius to values more consistent with the Voyager 1 observations than in current global models.

Kinetic simulations of 3D collisionless magnetic reconnection with a guide field show a dramatic enhancement of energetic electron production when compared with 2D systems. In the 2D systems, electrons are trapped in magnetic islands that limit their energy gain, whereas in the 3D systems the filamentation of the current layer leads to a stochastic magnetic field that enables the electrons to access volume-filling acceleration regions. The dominant accelerator of the most energetic electrons is a Fermi-like mechanism associated with reflection of charged particles from contracting field lines., Comment: 13 pages, 4 figures

A new regime of fast magnetic reconnection with an out-of-plane (guide) magnetic field is reported in which the key role is played by an electron pressure anisotropy described by the Chew-Goldberger-Low gyrotropic equations of state in the generalized Ohm's law, which even dominates the Hall term. A description of the physical cause of this behavior is provided and two-dimensional fluid simulations are used to confirm the results. The electron pressure anisotropy causes the out-of-plane magnetic field to develop a quadrupole structure of opposite polarity to the Hall magnetic field and gives rise to dispersive waves. In addition to being important for understanding what causes reconnection to be fast, this mechanism should dominate in plasmas with low plasma beta and a high in-plane plasma beta with electron temperature comparable to or larger than ion temperature, so it could be relevant in the solar wind and some tokamaks., Comment: 5 pages, 2 figures, accepted for publication as a letter in Physics of Plasmas

Electron bulk heating during magnetic reconnection with symmetric inflow conditions is examined using kinetic particle-in-cell (PIC) simulations. The degree of electron heating is well correlated with the inflowing Alfv\'en speed $c_{Ar}$ based on the reconnecting magnetic field through the relation $\Delta T_e = 0.033 \,m_i\,c_{Ar}^2$, where $\Delta T_{e}$ is the increase in electron temperature. For the range of simulations performed, the heating shows almost no correlation with inflow total temperature $T_{tot} = T_i + T_e$ or plasma $\beta$. An out-of-plane (guide) magnetic field of similar magnitude to the reconnecting field does not affect the total heating, but it does quench perpendicular heating, with almost all heating being in the parallel direction. These results are qualitatively consistent with a recent statistical survey of electron heating in the dayside magnetopause, which also found that $\Delta T_e$ was proportional to the inflowing Alfv\'en speed. The net electron heating varies very little with distance downstream of the x-line. The simulations show at most a very weak dependence of electron heating on the ion to electron mass ratio. In the antiparallel reconnection case, the largely parallel heating is eventually isotropized downstream due a scattering mechanism such as stochastic particle motion or instabilities. The study highlights key properties that must be satisfied by an electron heating mechanism: (1) Preferential heating in the parallel direction; (2) Heating proportional to $m_i\,c_{Ar}^2$; (3) At most a weak dependence on electron mass; and (4) An exhaust electron temperature that varies little with distance from the x-line., Comment: 27 pages, 1 table, 7 figures. Submitted to Physics of Plasmas

Observations of flare-heated electrons in the corona typically suggest confinement of electrons. The confinement mechanism, however, remains unclear. The transport of coronal hot electrons into ambient plasma was recently investigated by particle-in-cell (PIC) simulations. Electron transport was significantly suppressed by the formation of a highly localized, nonlinear electrostatic electric potential in the form of a double layer (DL). In this work large-scale PIC simulations are performed to explore the dynamics of DLs in larger systems where, instead of a single DL, multiple DLs are generated. The primary DL accelerates return current electrons, resulting in high velocity electron beams that interact with ambient ions. This forms a Buneman unstable system that spawns more DLs. Trapping of heated return current electrons between multiple DLs strongly suppresses electron transport. DLs also accelerate ambient ions and produce strong ion flows over an extended region. This clarifies the mechanism by which hot electrons in the corona couple to and accelerate ions to form the solar wind. These new dynamics in larger systems reveal a more likely picture of DL development and their impact on the ambient plasma in the solar corona. They are applicable to the preparation for in-situ coronal space missions like the Solar Probe Plus.

The heating of electrons in collisionless magnetic reconnection is explored in particle-in-cell (PIC) simulations with non-zero guide fields so that electrons remain magnetized. In this regime electric fields parallel to B accelerate particles directly while those perpendicular to B do so through gradient-B and curvature drifts. The curvature drift drives parallel heating through Fermi reflection while the gradient B drift changes the perpendicular energy through betatron acceleration. We present simulations in which we evaluate each of these mechanisms in space and time in order to quantify their role in electron heating. For a case with a small guide field (20 % of the magnitude of the reconnecting component) the curvature drift is the dominant source of electron heating. However, for a larger guide field (equal to the magnitude of the reconnecting component) electron acceleration by the curvature drift is comparable to that of the parallel electric field. In both cases the heating by the gradient B drift is negligible in magnitude. Heating by the curvature-drift dominates in the outflow exhausts where bent field lines expand to relax their tension and is therefore distributed over a large area. In contrast, the parallel electric field is localized near X-lines. Acceleration by parallel electric fields may play a smaller role in large systems where the X-line occupies a vanishing fraction of the system. The curvature drift and the parallel electric field dominate the dynamics and drive parallel heating. A consequence is that the electron energy spectrum becomes extremely anisotropic at late time, which has important implications for quantifying the limits of electron acceleration due to synchrotron emission. An upper limit on electron energy gain that is substantially higher than earlier estimates is obtained by balancing reconnection drive with radiative loss., Comment: Submitted to Physics of Plasmas (19 pages, 13 figures)

The onset of the acceleration of ions during magnetic reconnection is explored via particle-in-cell simulations in the limit of a strong ambient guide field that self-consistently and simultaneously follow the motions of protons and $\alpha$ particles. Heating parallel to the local magnetic field during reconnection with a guide field is strongly reduced compared with the reconnection of anti-parallel magnetic fields. The dominant heating of thermal ions during guide field reconnection results from pickup behavior of ions during their entry into reconnection exhausts and dominantly produces heating perpendicular rather than parallel to the local magnetic field. Pickup behavior requires that the ion transit time across the exhaust boundary (with a transverse scale of the order of the ion sound Larmor radius) be short compared with the ion cyclotron period. This translates into a threshold in the strength of reconnecting magnetic field that favors the heating of ions with high mass-to-charge. A simulation with a broad initial current layer produces a reconnecting system in which the amplitude of the reconnecting magnetic field just upstream of the dissipation region increases with time as reconnection proceeds. The sharp onset of perpendicular heating when the pickup threshold is crossed is documented. A comparison of the time variation of the parallel and perpendicular ion heating with that predicted based on the strength of the reconnecting field establishes the scaling of ion heating with ambient parameters both below and above the pickup threshold. The relevance to observations of ion heating in the solar corona is discussed., Comment: arXiv admin note: text overlap with arXiv:1108.5750