Search

The linear stability of waves driven by ion beams produced during solar flare energy release are explored to assess their role in driving abundance enhancements in minority species such as $^3$He and in controlling, through pitch-angle scattering, proton/alpha confinement during energy release. The Arbitrary Linear Plasma Solver (ALPS) is used to solve the linear dispersion relation for a population of energetic, reconnection-accelerated protons streaming into a cold background plasma. We assume equal densities of the two populations, using an anisotropic ($T_\parallel/T_\perp = 10$), one-sided kappa distribution for the energetic streaming population and a cold Maxwellian for the background. We find two unstable modes with parallel propagation: a right-handed wave with a frequency of the order of the proton cyclotron frequency ($\Omega_{cp}$) and a left-handed, lower frequency mode. We also find highly oblique modes with frequencies below $\Omega_{cp}$ that are unstable for higher beam energies. Through resonant interactions, all three modes will contribute to the scattering of the high-energy protons, thereby limiting their transport out of the flare-acceleration region. The higher-frequency oblique mode, which can be characterized as a kinetic Alfv\'en wave, will preferentially heat $^3$He, making it a good candidate for the driver of the abundance enhancements commonly observed for this species in impulsive events., Comment: 19 pages, 9 figures, submitted to Physics of Plasmas

We report observations of direct evidence of energetic protons being accelerated above ~400 keV within the reconnection exhaust of a heliospheric current sheet (HCS) crossing by NASA's Parker Solar Probe (PSP) at a distance of ~16.25 solar radii (Rs) from the Sun. Inside the extended exhaust, both the reconnection-generated plasma jets and the accelerated protons propagated toward the Sun, unambiguously establishing their origin from HCS reconnection sites located beyond PSP. Within the core of the exhaust, PSP detected stably trapped energetic protons up to ~400 keV, which is approximately 1000 times greater than the available magnetic energy per particle. The differential energy spectrum of the accelerated protons behaved as a pure power-law with spectral index of about -5. Supporting simulations using the kglobal model suggest that the trapping and acceleration of protons up to ~400 keV in the reconnection exhaust is likely facilitated by merging magnetic islands with a guide field between ~0.2-0.3 of the reconnecting magnetic field, consistent with the observations. These new results, enabled by PSP's proximity to the Sun, demonstrate that magnetic reconnection in the HCS is a significant new source of energetic particles in the near-Sun solar wind. The discovery of in-situ particle acceleration via magnetic reconnection at the HCS provides valuable insights into this fundamental process which frequently converts the large magnetic field energy density in the near-Sun plasma environment and may be responsible for heating the sun's atmosphere, accelerating the solar wind, and energizing charged particles to extremely high energies in solar flares., Comment: 24 pages, 4 figures

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

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

A set of equations are developed that extend the macroscale magnetic reconnection simulation model kglobal to include particle ions. The extension from earlier versions of kglobal, which included only particle electrons, requires the inclusion of the inertia of particle ions in the fluid momentum equation, which was not required in the electron-only model. The new equations will facilitate the exploration of the simultaneous non-thermal energization of ions and electrons during magnetic reconnection in macroscale systems. Numerical tests of the propagation of Alfv\'en waves in a plasma with anisotropic electron and ion pressure are presented to benchmark the new model., Comment: 10 pages, 2 figures

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

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.

A variety of energy sources, ranging from dynamic processes like magnetic reconnection and waves to quasi-steady terms like the plasma pressure, may contribute to the acceleration of the solar wind. We utilize a combination of charged particle and magnetic field observations from the Parker Solar Probe (PSP) to attempt to quantify the steady-state contribution of the proton pressure, the electric potential, and the wave energy to the solar wind proton acceleration observed by PSP between 13.3 and ~100 solar radii (RS). The proton pressure provides a natural kinematic driver of the outflow. The ambipolar electric potential acts to couple the electron pressure to the protons, providing another definite proton acceleration term. Fluctuations and waves, while inherently dynamic, can act as an additional effective steady-state pressure term. To analyze the contributions of these terms, we utilize radial binning of single-point PSP measurements, as well as repeated crossings of the same stream at different distances on individual PSP orbits (i.e. "fast radial scans"). In agreement with previous work, we find that the electric potential contains sufficient energy to fully explain the acceleration of the slower wind streams. On the other hand, we find that the wave pressure plays an increasingly important role in the faster wind streams. The combination of these terms can explain the continuing acceleration of both slow and fast wind streams beyond 13.3 RS.

We present EUV solar observations showing evidence for omnipresent jetting activity driven by small-scale magnetic reconnection at the base of the solar corona. We argue that the physical mechanism that heats and drives the solar wind at its source is ubiquitous magnetic reconnection in the form of small-scale jetting activity (i.e., a.k.a. jetlets). This jetting activity, like the solar wind and the heating of the coronal plasma, are ubiquitous regardless of the solar cycle phase. Each event arises from small-scale reconnection of opposite polarity magnetic fields producing a short-lived jet of hot plasma and Alfv\'en waves into the corona. The discrete nature of these jetlet events leads to intermittent outflows from the corona, which homogenize as they propagate away from the Sun and form the solar wind. This discovery establishes the importance of small-scale magnetic reconnection in solar and stellar atmospheres in understanding ubiquitous phenomena such as coronal heating and solar wind acceleration. Based on previous analyses linking the switchbacks to the magnetic network, we also argue that these new observations might provide the link between the magnetic activity at the base of the corona and the switchback solar wind phenomenon. These new observations need to be put in the bigger picture of the role of magnetic reconnection and the diverse form of jetting in the solar atmosphere., Comment: 10 pages, 4 figures

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

We present observations of >10-100 keV/nucleon suprathermal (ST) H, He, O, and Fe ions associated with crossings of the heliospheric current sheet (HCS) at radial distances <0.1 au from the Sun. Our key findings are: 1) very few heavy ions are detected during the 1st full crossing, the heavy ion intensities are reduced during the 2nd partial crossing and peak just after the 2nd crossing; 2) ion arrival times exhibit no velocity dispersion; 3) He pitch-angle distributions track the magnetic field polarity reversal and show up to ~10:1 anti-sunward, field-aligned flows and beams closer to the HCS that become nearly isotropic further from the HCS; 4) the He spectrum steepens either side of the HCS and the He, O, and Fe spectra exhibit power-laws of the form ~E^4-6; and 5) maximum energies EX increase with the ion's charge-to-mass (Q/M) ratio as EX/EH proportional to [(QX/MX)]^alpha where alpha~0.65-0.76, assuming that the average Q-states are similar to those measured in gradual and impulsive solar energetic particle events at 1 au. The absence of velocity dispersion in combination with strong field-aligned anisotropies closer to the HCS appears to rule out solar flares and near-sun coronal mass ejection-driven shocks. These new observations present challenges not only for mechanisms that employ direct parallel electric fields and organize maximum energies according to E/Q, but also for local diffusive and magnetic reconnection-driven acceleration models. Re-evaluation of our current understanding of the production and transport of energetic ions is necessary to understand this near-solar, current-sheet-associated population of ST ions., Comment: 4 Figures, 2 Tables

Global models of the heliosphere are critical tools used in the interpretation of heliospheric observations. There are several three-dimensional magnetohydrodynamic (MHD) heliospheric models that rely on different strategies and assumptions. Until now only one paper has compared global heliosphere models, but without magnetic field effects. We compare the results of two different MHD models, the BU and Moscow models. Both models use identical boundary conditions to compare how different numerical approaches and physical assumptions contribute to the heliospheric solution. Based on the different numerical treatments of discontinuities, the BU model allows for the presence of magnetic reconnection, while the Moscow model does not. Both models predict collimation of the solar outflow in the heliosheath by the solar magnetic field and produce a split-tail where the solar magnetic field confines the charged solar particles into distinct north and south columns that become lobes. In the BU model, the ISM flows between the two lobes at large distances due to MHD instabilities and reconnection. Reconnection in the BU model at the port flank affects the draping of the interstellar magnetic field in the immediate vicinity of the heliopause. Different draping in the models cause different ISM pressures, yielding different heliosheath thicknesses and boundary locations, with the largest effects at high latitudes. The BU model heliosheath is 15% thinner and the heliopause is 7% more inwards at the north pole relative to the Moscow model. These differences in the two plasma solutions may manifest themselves in energetic neutral atom measurements of the heliosphere., Comment: 2 tables, 12 figures, accepted to ApJ

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

One of the striking observations from the Parker Solar Probe (PSP) spacecraft is the prevalence in the inner heliosphere of large amplitude, Alfv\'enic magnetic field reversals termed 'switchbacks'. These $\delta B_R/B \sim \mathcal{O}(1$) fluctuations occur on a range of timescales and in {\em patches} separated by intervals of quiet, radial magnetic field. We use measurements from PSP to demonstrate that patches of switchbacks are localized within the extensions of plasma structures originating at the base of the corona. These structures are characterized by an increase in alpha particle abundance, Mach number, plasma $\beta$ and pressure, and by depletions in the magnetic field magnitude and electron temperature. These intervals are in pressure-balance, implying stationary spatial structure, and the field depressions are consistent with overexpanded flux tubes. The structures are asymmetric in Carrington longitude with a steeper leading edge and a small ($\sim$1$^\circ$) edge of hotter plasma and enhanced magnetic field fluctuations. Some structures contain suprathermal ions to $\sim$85 keV that we argue are the energetic tail of the solar wind alpha population. The structures are separated in longitude by angular scales associated with supergranulation. This suggests that these switchbacks originate near the leading edge of the diverging magnetic field funnels associated with the network magnetic field - the primary wind sources. We propose an origin of the magnetic field switchbacks, hot plasma and suprathermals, alpha particles in interchange reconnection events just above the solar transition region and our measurements represent the extended regions of a turbulent outflow exhaust., Comment: 18 pages, 9 figures, submitted to Astrophys. J

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.

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

Traditionally, the solar magnetic field has been considered to have a negligible effect in the outer regions of the heliosphere. Recent works have shown that the solar magnetic field may play a crucial role in collimating the plasma in the heliosheath. Interstellar Boundary Explorer (IBEX) observations of the heliotail indicated a latitudinal structure varying with energy in the energetic neutral atom (ENA) fluxes. At energies ~1 keV, the ENA fluxes show an enhancement at low latitudes and a deficit of ENAs near the poles. At energies >2.7 keV, ENA fluxes had a deficit within low latitudes, and lobes of higher ENA flux near the poles. This ENA structure was initially interpreted to be a result of the latitudinal profile of the solar wind during solar minimum. We extend the work of Kornbleuth et al. (2018) by using solar minimum-like conditions and the recently developed SHIELD model. The SHIELD model couples the magnetohydrodynamic (MHD) plasma solution with a kinetic description of neutral hydrogen. We show that while the latitudinal profile of the solar wind during solar minimum contributes to the lobes in ENA maps, the collimation by the solar magnetic field is important in creating and shaping the two high latitude lobes of enhanced ENA flux observed by IBEX. This is the first work to explore the effect of the changing solar magnetic field strength on ENA maps. Our findings suggest that IBEX is providing the first observational evidence of the collimation of the heliosheath plasma by the solar magnetic field., Comment: Accepted for publication in ApJL

Observations by the Parker Solar Probe mission of the solar wind at about 35.7 solar radii reveal the existence of whistler wave packets with frequencies below 0.1 f/fce (20-80 Hz in the spacecraft frame). These waves often coincide with local minima of the magnetic field magnitude or with sudden deflections of the magnetic field that are called switchbacks. Their sunward propagation leads to a significant Doppler frequency downshift from 200-300 Hz to 20-80 Hz (from 0.2 f/fce to 0.5 f/fce). The polarization of these waves varies from quasi-parallel to significantly oblique with wave normal angles that are close to the resonance cone. Their peak amplitude can be as large as 2 to 4 nT. Such values represent approximately 10% of the background magnetic field, which is considerably more than what is observed at 1 a.u. Recent numerical studies show that such waves may potentially play a key role in breaking the heat flux and scattering the Strahl population of suprathermal electrons into a halo population.

Magnetic reconnection in current sheets converts magnetic energy into particle energy. The process may play an important role in the acceleration and heating of the solar wind close to the Sun. Observations from Parker Solar Probe provide a new opportunity to study this problem, as it measures the solar wind at unprecedented close distances to the Sun. During the 1st orbit, PSP encountered a large number of current sheets in the solar wind through perihelion at 35.7 solar radii. We performed a comprehensive survey of these current sheets and found evidence for 21 reconnection exhausts. These exhausts were observed in heliospheric current sheets, coronal mass ejections, and regular solar wind. However, we find that the majority of current sheets encountered around perihelion, where the magnetic field was strongest and plasma beta was lowest, were Alfv\'enic structures associated with bursty radial jets and these current sheets did not appear to be undergoing local reconnection. We examined conditions around current sheets to address why some current sheets reconnected, while others did not. A key difference appears to be the degree of plasma velocity shear across the current sheets: The median velocity shear for the 21 reconnection exhausts was 24% of the Alfv\'en velocity shear, whereas the median shear across 43 Alfv\'enic current sheets examined was 71% of the Alfv\'en velocity shear. This finding could suggest that large, albeit sub-Alfv\'enic, velocity shears suppress reconnection. An alternative interpretation is that the Alfv\'enic current sheets are isolated rotational discontinuities which do not undergo local reconnection., Comment: In Press (accepted by ApJS on 2019-11-08) This paper is part of the Parker Solar Probe ApJS Special Issue Current citation: Phan, T. D., S. D. Bale, J. P. Eastwood, et al. (2019), Parker Solar Probe In-Situ Observations of Magnetic Reconnection Exhausts During Encounter 1, ApJS., in press, doi: 10.3847/1538-4365/ab55ee

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.

Waves around the lower hybrid frequency are frequently observed at Earth's magnetopause, and readily reach very large amplitudes. Determining the properties of lower hybrid waves is crucial because they are thought to contribute to electron and ion heating, cross-field particle diffusion, anomalous resistivity, and energy transfer between electrons and ions. All these processes could play an important role in magnetic reconnection at the magnetopause and the evolution of the boundary layer. In this paper, the properties of lower hybrid waves at Earth's magnetopause are investigated using the Magnetospheric Multiscale (MMS) mission. For the first time, the properties of the waves are investigated using fields and direct particle measurements. The highest-resolution electron moments resolve the velocity and density fluctuations of lower hybrid waves, confirming that electrons remain approximately frozen in at lower hybrid wave frequencies. Using fields and particle moments the dispersion relation is constructed and the wave-normal angle is estimated to be close to $90^{\circ}$ to the background magnetic field. The waves are shown to have a finite parallel wave vector, suggesting that they can interact with parallel propagating electrons. The observed wave properties are shown to agree with theoretical predictions, the previously used single-spacecraft method, and four-spacecraft timing analyses. These results show that single-spacecraft methods can accurately determine lower hybrid wave properties., Comment: 43 pages, 18 figures

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.

Magnetic reconnection is an energy conversion process important in many astrophysical contexts including the Earth's magnetosphere, where the process can be investigated in-situ. Here we present the first encounter of a reconnection site by NASA's Magnetospheric Multiscale (MMS) spacecraft in the magnetotail, where reconnection involves symmetric inflow conditions. The unprecedented electron-scale plasma measurements revealed (1) super-Alfvenic electron jets reaching 20,000 km/s, (2) electron meandering motion and acceleration by the electric field, producing multiple crescent-shaped structures, (3) spatial dimensions of the electron diffusion region implying a reconnection rate of 0.1-0.2. The well-structured multiple layers of electron populations indicate that, despite the presence of turbulence near the reconnection site, the key electron dynamics appears to be largely laminar., Comment: 4 pages, 3 figures, and supplementary material

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 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

We perform a theoretical and numerical study of anti-parallel 2D magnetic reconnection with asymmetries in the density and reconnecting magnetic field strength in addition to a bulk flow shear across the reconnection site in the plane of the reconnecting fields, which commonly occurs at planetary magnetospheres. We predict the speed at which an isolated X-line is convected by the flow, the reconnection rate, and the critical flow speed at which reconnection no longer takes place for arbitrary reconnecting magnetic field strengths, densities, and upstream flow speeds, and confirm the results with two-fluid numerical simulations. The predictions and simulation results counter the prevailing model of reconnection at Earth's dayside magnetopause which says reconnection occurs with a stationary X-line for sub-Alfvenic magnetosheath flow, reconnection occurs but the X-line convects for magnetosheath flows between the Alfven speed and double the Alfven speed, and reconnection does not occur for magnetosheath flows greater than double the Alfven speed. We find that X-line motion is governed by momentum conservation from the upstream flows, which are weighted differently in asymmetric systems, so the X-line convects for generic conditions including sub-Alfvenic upstream speeds. For the reconnection rate, while the cutoff condition for symmetric reconnection is that the difference in flows on the two sides of the reconnection site is twice the Alfven speed, we find asymmetries cause the cutoff speed for asymmetric reconnection to be higher than twice the asymmetric form of the Alfven speed. The results compare favorably with an observation of reconnection at Earth's polar cusps during a period of northward interplanetary magnetic field, where reconnection occurs despite the magnetosheath flow speed being more than twice the magnetosheath Alfven speed, the previously proposed suppression condition., Comment: 46 pages, 7 figures, abstract abridged here, accepted to Journal of Geophysical Research - Space Physics

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.

The physical processes that control the partition of released magnetic energy between electrons and ions during reconnection is explored through particle-in-cell simulations and analytical techniques. We demonstrate that the development of a large-scale parallel electric field and its associated potential controls the relative heating of electrons and ions. The potential develops to restrain heated exhaust electrons and enhances their heating by confining electrons in the region where magnetic energy is released. Simultaneously the potential slows ions entering the exhaust below the Alfv\'enic speed expected from the traditional counterstreaming picture of ion heating. Unexpectedly, the magnitude of the potential and therefore the relative partition of energy between electrons and ions is not a constant but rather depends on the upstream parameters and specifically the upstream electron normalized temperature (electron beta). These findings suggest that the fraction of magnetic energy converted into the total thermal energy may be independent of upstream parameters.

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

The classic accepted view of the heliosphere is a quiescent, comet-like shape aligned in the direction of the Sun's travel through the interstellar medium (ISM) extending for 1000's of AUs (AU: astronomical unit). Here we show, based on magnetohydrodynamic (MHD) simulations, that the tension (hoop) force of the twisted magnetic field of the sun confines the solar wind plasma beyond the termination shock and drives jets to the North and South very much like astrophysical jets. These jets are deflected into the tail region by the motion of the Sun through the ISM similar to bent galactic jets moving through the intergalactic medium. The interstellar wind blows the two jets into the tail but is not strong enough to force the lobes into a single comet-like tail, as happens to some astrophysical jets (Morsony et al. 2013). Instead, the interstellar wind flows around the heliosphere and into equatorial region between the two jets. As in some astrophysical jets that are kink unstable (Porth et al. 2014) we show here that the heliospheric jets are turbulent (due to large-scale MHD instabilities and reconnection) and strongly mix the solar wind with the ISM beyond 400 AU. The resulting turbulence has important implications for particle acceleration in the heliosphere. The two-lobe structure is consistent with the energetic neutral atoms (ENAs) images of the heliotail from IBEX (McComas et al. 2013) where two lobes are visible in the North and South and the suggestion from the CASSINI (Krimigis et al. 2009, Dialynas et al. 2013) ENAs that the heliosphere is "tailless"., Comment: 19 pages, 5 figures; Astrophysical Journal Letters; in press

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.