Showing 1–50 of 76 results for author: Swisdak, M

Simulation Models for Exploring Magnetic Reconnection

Authors: Michael Shay, Subash Adhikari, Naoki Beesho, Joachim Birn, Jorg Buechner, Paul Cassak, Li-Jen Chen, Yuxi Chen, Giulia Cozzani, Jim Drake, Fan Guo, Michael Hesse, Neeraj Jain, Yann Pfau-Kempf, Yu Lin, Yi-Hsin Liu, Mitsuo Oka, Yuri A. Omelchenko, Minna Palmroth, Oreste Pezzi, Patricia H. Reiff, Marc Swisdak, Frank Toffoletto, Gabor Toth, Richard A. Wolf

Abstract: Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This chapter reviews many of these… ▽ More Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This chapter reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM). △ Less

Submitted 9 June, 2024; originally announced June 2024.

Comments: Chapter 5.2 of ISSI Book on Magnetic Reconnection, submitted to Space Science Reviews

Ohm's Law, the Reconnection Rate, and Energy Conversion in Collisionless Magnetic Reconnection

Authors: Yi-Hsin Liu, Michael Hesse, Kevin Genestreti, Rumi Nakamura, Jim Burch, Paul Cassak, Naoki Bessho, Jonathan Eastwood, Tai Phan, Marc Swisdak, Sergio Toledo-Redondo, Masahiro Hoshino, Cecilia Norgren, Hantao Ji, TKM Nakamura

Abstract: Magnetic reconnection is a ubiquitous plasma process that transforms magnetic energy into particle energy during eruptive events throughout the universe. Reconnection not only converts energy during solar flares and geomagnetic substorms that drive space weather near Earth, but it may also play critical roles in the high energy emissions from the magnetospheres of neutron stars and black holes. In… ▽ More Magnetic reconnection is a ubiquitous plasma process that transforms magnetic energy into particle energy during eruptive events throughout the universe. Reconnection not only converts energy during solar flares and geomagnetic substorms that drive space weather near Earth, but it may also play critical roles in the high energy emissions from the magnetospheres of neutron stars and black holes. In this review article, we focus on collisionless plasmas that are most relevant to reconnection in many space and astrophysical plasmas. Guided by first-principles kinetic simulations and spaceborne in-situ observations, we highlight the most recent progress in understanding this fundamental plasma process. We start by discussing the non-ideal electric field in the generalized Ohm's law that breaks the frozen-in flux condition in ideal magnetohydrodynamics and allows magnetic reconnection to occur. We point out that this same reconnection electric field also plays an important role in sustaining the current and pressure in the current sheet and then discuss the determination of its magnitude (i.e., the reconnection rate), based on force balance and energy conservation. This approach to determining the reconnection rate is applied to kinetic current sheets of a wide variety of magnetic geometries, parameters, and background conditions. We also briefly review the key diagnostics and modeling of energy conversion around the reconnection diffusion region, seeking insights from recently developed theories. Finally, future prospects and open questions are discussed. △ Less

Submitted 2 June, 2024; originally announced June 2024.

Comments: Submitted to Space Science Reviews. This is a review paper as an outcome of the 2022 Magnetic Reconnection Workshop in the International Space Science Institute

A Computational Model for Ion and Electron Energization during Macroscale Magnetic Reconnection

Authors: Zhiyu Yin, J. F. Drake, Marc Swisdak

Abstract: 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 exp… ▽ More 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én waves in a plasma with anisotropic electron and ion pressure are presented to benchmark the new model. △ Less

Submitted 20 May, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Comments: 10 pages, 2 figures

Constraints on the Alfvénicity of Switchbacks

Authors: O. V. Agapitov, J. F. Drake, M. Swisdak, K. -E. Choi, N. Raouafi

Abstract: 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 sp… ▽ More 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énicity by evaluating the distribution of the parameter defined as the ratio of the parallel to $Δ\vec{B}$ component of $Δ\vec{V}$ to $Δ\vec{V}_A = Δ\vec{B}/4πn_i m_i$ inside SBs, i.e. $α=V_{||}/|Δ\vec{V}_A|$ ( $α=Δ\vec{V}/|Δ\vec{V}_A|$ when $Δ\vec{V}~Δ\vec{B}$), which quantifies the deviation of the perturbation from an Alfvénic one. Based on Parker Solar Probe (PSP) observations, we show that $α$ inside SBs has systematically lower values than it has in the pristine solar wind: $α$ inside SBs observed during PSP Encounter 1 were distributed in a range of 0.2-0.9. The upper limit on $α$ is constrained by the requirement that the jump in velocity across the switchback boundary be less than the local Alfvén 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 $α=1$ (the Alfvénicity condition) and the highest possible $α$ 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én waves. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Wave Generation by Flare-Accelerated Ions and Implications for 3He Acceleration

Authors: A. Fitzmaurice, J. F. Drake, M. Swisdak

Abstract: 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 b… ▽ More 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. △ Less

Submitted 15 March, 2024; v1 submitted 20 November, 2023; originally announced November 2023.

Comments: 19 pages, 12 figures

Advanced methods for analyzing in-situ observations of magnetic reconnection

Authors: H. Hasegawa, M. R. Argall, N. Aunai, R. Bandyopadhyay, N. Bessho, I. J. Cohen, R. E. Denton, J. C. Dorelli, J. Egedal, S. A. Fuselier, P. Garnier, V. Genot, D. B. Graham, K. J. Hwang, Y. V. Khotyaintsev, D. B. Korovinskiy, B. Lavraud, Q. Lenouvel, T. C. Li, Y. -H. Liu, B. Michotte de Welle, T. K. M. Nakamura, D. S. Payne, S. M. Petrinec, Y. Qi , et al. (11 additional authors not shown)

Abstract: 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 a… ▽ More 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. △ Less

Submitted 24 June, 2024; v1 submitted 11 July, 2023; originally announced July 2023.

Comments: submitted to Space Science Reviews (116 pages, incl. 31 figures, 7 tables)

Wave generation and energetic electron scattering in solar flares

Authors: Hanqing Ma, James F. Drake, Marc Swisdak

Abstract: We conduct two-dimensional particle-in-cell simulations to investigate the scattering of electron heat flux by self-generated oblique electromagnetic waves. The heat flux is modeled as a bi-kappa distribution with a T_parallel > T_perp temperature anisotropy maintained by continuous injection at the boundaries. The anisotropic distribution excites oblique whistler waves and filamentary-like Weibel… ▽ More We conduct two-dimensional particle-in-cell simulations to investigate the scattering of electron heat flux by self-generated oblique electromagnetic waves. The heat flux is modeled as a bi-kappa distribution with a T_parallel > T_perp temperature anisotropy maintained by continuous injection at the boundaries. The anisotropic distribution excites oblique whistler waves and filamentary-like Weibel instabilities. Electron velocity distributions taken after the system has reached a steady state show that these in stabilities inhibit the heat flux and drive the total distributions towards isotropy. Electron trajectories in velocity space show a circular-like diffusion along constant energy surfaces in the wave frame. The key parameter controlling the scattering rate is the average speed, or drift speed vd, of the heat flux compared with the electron Alfven speed vAe, with higher drift speeds producing stronger fluctua tions and a more significant reduction of the heat flux. Reducing the density of the electrons carrying the heat flux by 50% does not significantly affect the scattering rate. A scaling law for the electron scattering rate versus vd/vAe is deduced from the simulations. The implications of these results for understanding energetic electron transport during solar flare energy release are discussed. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Development of a structured, turbulent solar wind as a result of interchange reconnection

Authors: J. F. Drake, S. D. Bale, M. Swisdak, N. E. Raouafi, M. Velli

Abstract: 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énic outflow and associated heating differ from previous models. Additional heating of the corona associated with Alfvén waves or other MHD turbulence,… ▽ More 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énic outflow and associated heating differ from previous models. Additional heating of the corona associated with Alfvén 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énic ejection and reconnection-driven heating yields a minimum value of the Alfvén 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énic 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énic flows that suggests these flows are the free-energy source of the Alfvénic 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. △ Less

Submitted 15 August, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Advancing Theory and Modeling Efforts in Heliophysics

Authors: Fan Guo, Spiro Antiochos, Paul Cassak, Bin Chen, Xiaohang Chen, Chuanfei Dong, Cooper Downs, Joe Giacalone, Colby C. Haggerty, Hantao Ji, Judith Karpen, James Klimchuk, Wen Li, Xiaocan Li, Mitsuo Oka, Katharine K. Reeves, Marc Swisdak, Weichao Tu

Abstract: Heliophysics theory and modeling build understanding from fundamental principles to motivate, interpret, and predict observations. Together with observational analysis, they constitute a comprehensive scientific program in heliophysics. As observations and data analysis become increasingly detailed, it is critical that theory and modeling develop more quantitative predictions and iterate with obse… ▽ More Heliophysics theory and modeling build understanding from fundamental principles to motivate, interpret, and predict observations. Together with observational analysis, they constitute a comprehensive scientific program in heliophysics. As observations and data analysis become increasingly detailed, it is critical that theory and modeling develop more quantitative predictions and iterate with observations. Advanced theory and modeling can inspire and greatly improve the design of new instruments and increase their chance of success. In addition, in order to build physics-based space weather forecast models, it is important to keep developing and testing new theories, and maintaining constant communications with theory and modeling. Maintaining a sustainable effort in theory and modeling is critically important to heliophysics. We recommend that all funding agencies join forces and consider expanding current and creating new theory and modeling programs--especially, 1. NASA should restore the HTMS program to its original support level to meet the critical needs of heliophysics science; 2. a Strategic Research Model program needs to be created to support model development for next-generation basic research codes; 3. new programs must be created for addressing mission-critical theory and modeling needs; and 4. enhanced programs are urgently required for training the next generation of theorists and modelers. △ Less

Submitted 8 September, 2022; originally announced September 2022.

Comments: White paper submitted to Heliophysics 2024 Decadal Survey

Interchange reconnection as the source of the fast solar wind within coronal holes

Authors: S. D. Bale, J. F. Drake, M. D. McManus, M. I. Desai, S. T. Badman, D. E. Larson, M. Swisdak, T. S. Horbury, N. E. Raouafi, T. Phan, M. Velli, D. J. McComas, C. M. S. Cohen, D. Mitchell, O. Panasenco, J. C. Kasper

Abstract: 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 coron… ▽ More 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. △ Less

Submitted 7 June, 2023; v1 submitted 16 August, 2022; originally announced August 2022.

Comments: 15 pages, 4 figures. Nature, 2023

Scaling of electron heating by magnetization during reconnection and applications to dipolarization fronts and super-hot solar flares

Authors: M. Hasan Barbhuiya, Paul. A. Cassak, Michael. A. Shay, Vadim Roytershteyn, Marc Swisdak, Amir Caspi, Andrei Runov, Haoming Liang

Abstract: Electron ring velocity space distributions have previously been seen in numerical simulations of magnetic reconnection exhausts and have been suggested to be caused by the magnetization of the electron outflow jet by the compressed reconnected magnetic fields [Shuster et al., ${\it Geophys.~Res.~Lett.}, {\bf 41}$, 5389 (2014)]. We present a theory of the dependence of the major and minor radii of… ▽ More Electron ring velocity space distributions have previously been seen in numerical simulations of magnetic reconnection exhausts and have been suggested to be caused by the magnetization of the electron outflow jet by the compressed reconnected magnetic fields [Shuster et al., ${\it Geophys.~Res.~Lett.}, {\bf 41}$, 5389 (2014)]. We present a theory of the dependence of the major and minor radii of the ring distributions solely in terms of upstream (lobe) plasma conditions, thereby allowing a prediction of the associated temperature and temperature anisotropy of the rings in terms of upstream parameters. We test the validity of the prediction using 2.5-dimensional particle-in-cell (PIC) simulations with varying upstream plasma density and temperature, finding excellent agreement between the predicted and simulated values. We confirm the Shuster et al. suggestion for the cause of the ring distributions, and also find that the ring distributions are located in a region marked by a plateau, or shoulder, in the reconnected magnetic field profile. The predictions of the temperature are consistent with observed electron temperatures in dipolarization fronts, and may provide an explanation for the generation of plasma with temperatures in the 10s of MK in super-hot solar flares. A possible extension of the model to dayside reconnection is discussed. Since ring distributions are known to excite whistler waves, the present results should be useful for quantifying the generation of whistler waves in reconnection exhausts. △ Less

Submitted 31 July, 2022; originally announced August 2022.

Comments: Accepted for publication in Journal of Geophysical Research: Space Physics

Journal ref: Journal of Geophysical Research: Space Physics, Vol. 127, Issue 8, e2022JA030610 (22pp); 2022 August 08

Correlated Spatiotemporal Evolution of Extreme-Ultraviolet Ribbons and Hard X-rays in a Solar Flare

Authors: Stephen J. Naus, Jiong Qiu, C. Richard DeVore, Spiro K. Antiochos, Joel T. Dahlin, James F. Drake, Marc Swisdak

Abstract: We analyze the structure and evolution of ribbons from the M7.3 SOL2014-04-18T13 flare using ultraviolet (UV) images from IRIS and SDO/AIA, magnetic data from SDO/HMI, hard X-ray (HXR) images from RHESSI, and light curves from Fermi/GBM, in order to infer properties of coronal magnetic reconnection. As the event progresses, two flare ribbons spread away from the magnetic polarity inversion line. T… ▽ More We analyze the structure and evolution of ribbons from the M7.3 SOL2014-04-18T13 flare using ultraviolet (UV) images from IRIS and SDO/AIA, magnetic data from SDO/HMI, hard X-ray (HXR) images from RHESSI, and light curves from Fermi/GBM, in order to infer properties of coronal magnetic reconnection. As the event progresses, two flare ribbons spread away from the magnetic polarity inversion line. The width of the newly brightened front along the extension of the ribbon is highly intermittent in both space and time, presumably reflecting non-uniformities in the structure and/or dynamics of the flare current sheet. Furthermore, the ribbon width grows most rapidly in regions exhibiting concentrated non-thermal HXR emission, with sharp increases slightly preceding the HXR bursts. The light curve of the ultraviolet emission matches the HXR light curve at photon energies above 25 keV. In other regions the ribbon-width evolution and light curves do not temporally correlate with the HXR emission. This indicates that the production of non-thermal electrons is highly non-uniform within the flare current sheet. Our results suggest a strong connection between the production of non-thermal electrons and the locally enhanced perpendicular extent of flare ribbon fronts, which in turn reflects inhomogeneous structure and/or reconnection dynamics of the current sheet. Despite this variability, the ribbon fronts remain nearly continuous, quasi-one-dimensional features. Thus, although the reconnecting coronal current sheets are highly structured, they remain quasi-two-dimensional and the magnetic energy release occurs systematically, rather than stochastically, through the volume of reconnecting magnetic flux. △ Less

Submitted 5 December, 2021; v1 submitted 30 September, 2021; originally announced September 2021.

Slow Shock Formation Upstream of Reconnecting Current Sheets

Authors: Harry Arnold, James Drake, Marc Swisdak, Fan Guo, Joel Dahlin, Qile Zhang

Abstract: The formation, development and impact of slow shocks in the upstream region of reconnecting current layers are explored. Slow shocks have been documented in the upstream region of magnetohydrodynamic (MHD) simulations of magnetic reconnection as well as in similar simulations with the {\it kglobal} kinetic macroscale simulation model. They are therefore a candidate mechanism for preheating the pla… ▽ More The formation, development and impact of slow shocks in the upstream region of reconnecting current layers are explored. Slow shocks have been documented in the upstream region of magnetohydrodynamic (MHD) simulations of magnetic reconnection as well as in similar simulations with the {\it kglobal} kinetic macroscale simulation model. They are therefore a candidate mechanism for preheating the plasma that is injected into the current layers that facilitate magnetic energy release in solar flares. Of particular interest is their potential role in producing the hot thermal component of electrons in flares. During multi-island reconnection, the formation and merging of flux ropes in the reconnecting current layer drives plasma flows and pressure disturbances in the upstream region. These pressure disturbances steepen into slow shocks that propagate along the reconnecting component of the magnetic field and satisfy the expected Rankine-Hugoniot jump conditions. Plasma heating arises from both compression across the shock and the parallel electric field that develops to maintain charge neutrality in a kinetic system. Shocks are weaker at lower plasma $β$, where shock steepening is slow. While these upstream slow shocks are intrinsic to the dynamics of multi-island reconnection, their contribution to electron heating remains relatively minor compared with that from Fermi reflection and the parallel electric fields that bound the reconnection outflow. △ Less

Submitted 9 December, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 9 pages and 5 figures

Flux rope merging and the structure of switchbacks in the solar wind

Authors: O. Agapitov, J. F. Drake, M. Swisdak, S. D. Bale, T. S. Horbury, J. C. Kasper, R. J. MacDowall, F. S. Mozer, T. D. Phan, M. Pulupa, N. E. Raouafi, M. Velli

Abstract: 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 ejectio… ▽ More 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. △ Less

Submitted 3 December, 2021; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 21 pages, 9 figures

The Reversibility Of Magnetic Reconnection

Authors: M. Xuan, M. Swisdak, J. F. Drake

Abstract: 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 simulati… ▽ More 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. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Electron Acceleration during Macroscale Magnetic Reconnection

Authors: Harry Arnold, James Drake, Marc Swisdak, Fan Guo, Joel Dahlin, Bin Chen, Gregory Fleishman, Lindsay Glesener, Eduard Kontar, Tai Phan, Chengcai Shen

Abstract: The first self-consistent simulations of electron acceleration during magnetic reconnection in a macroscale system are presented. Consistent with solar flare observations the spectra of energetic electrons take the form of power-laws that extend more than two decades in energy. The drive mechanism for these nonthermal electrons is Fermi reflection in growing and merging magnetic flux ropes. A stro… ▽ More The first self-consistent simulations of electron acceleration during magnetic reconnection in a macroscale system are presented. Consistent with solar flare observations the spectra of energetic electrons take the form of power-laws that extend more than two decades in energy. The drive mechanism for these nonthermal electrons is Fermi reflection in growing and merging magnetic flux ropes. A strong guide field is found to suppress the production of nonthermal electrons by weakening the Fermi drive mechanism. For a weak guide field the total energy content of nonthermal electrons dominates that of the hot thermal electrons even though their number density remains small. Our results are benchmarked with the hard x-ray, radio and extreme ultra-violet (EUV) observations of the X8.2-class solar flare on September 10, 2017. △ Less

Submitted 22 January, 2021; v1 submitted 2 November, 2020; originally announced November 2020.

Journal ref: Phys. Rev. Lett. 126, 135101 (2021)

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Authors: H. Ji, J. Karpen, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, A. Bhattacharjee, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, B. Chen, L. -J. Chen, Y. Chen, A. Chien, L. Comisso, D. Craig, J. Dahlin, W. Daughton, E. DeLuca , et al. (83 additional authors not shown)

Abstract: 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 interagen… ▽ More 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. △ Less

Submitted 16 September, 2020; originally announced September 2020.

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

Switchbacks as signatures of magnetic flux ropes generated by interchange reconnection in the corona

Authors: J. F. Drake, O. Agapitov, M. Swisdak, S. T. Badman, S. D. Bale, T. S. Horbury, Justin C. Kasper, R. J. MacDowall, F. S. Mozer, T. D. Phan, M. Pulupa, A. Szabo, M. Velli

Abstract: 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 rope… ▽ More 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. △ Less

Submitted 12 October, 2020; v1 submitted 11 September, 2020; originally announced September 2020.

Journal ref: A&A 650, A2 (2021)

Whistler-regulated MHD: Transport equations for electron thermal conduction in the high $β$ intracluster medium of galaxy clusters

Authors: J. F. Drake, C. Pfrommer, C. S. Reynolds, M. Ruszkowski, M. Swisdak, A. Einarsson, T. Thomas, A. B. Hassam, G. T. Roberg-Clark

Abstract: Transport equations for electron thermal energy in the high $β_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… ▽ More Transport equations for electron thermal energy in the high $β_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}/β_e\ll{v}_{te}$ dominate the turbulent spectrum, with $v_{te}$ the electron thermal speed and $β_e\gg1$ the ratio of electron thermal to magnetic pressure. We find: (1) temperature-gradient-driven whistlers dominate classical scattering when $L_c>L/β_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. △ Less

Submitted 27 August, 2021; v1 submitted 15 July, 2020; originally announced July 2020.

Comments: Accepted for publication in The Astrophysical Journal

Noise-Induced Magnetic Field Saturation in Kinetic Simulations

Authors: J. Juno, M. Swisdak, J. M. TenBarge, V. Skoutnev, A. Hakim

Abstract: 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 m… ▽ More 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. △ Less

Submitted 7 September, 2020; v1 submitted 15 April, 2020; originally announced April 2020.

Comments: 15 pages, 4 figures

Journal ref: J. Plasma Phys. 86 (2020) 175860401

Turbulence and Transport During Guide-Field Reconnection at the Magnetopause

Authors: L. Price, M Swisdak, J. F. Drake, D. B. Graham

Abstract: 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 f… ▽ More 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. △ Less

Submitted 4 October, 2019; originally announced October 2019.

Scattering of energetic electrons by heat-flux-driven whistlers in flares

Authors: G. T. Roberg-Clark, O. V. Agapitov, J. F. Drake, M. M. Swisdak

Abstract: 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 re… ▽ More 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. △ Less

Submitted 25 October, 2019; v1 submitted 18 August, 2019; originally announced August 2019.

Comments: Accepted for publication in The Astrophysical Journal

Large-Scale Parallel Electric Fields and Return Currents in a Global Simulation Model

Authors: Harry Arnold, James Drake, Marc Swisdak, Joel Dahlin

Abstract: A new computational model, kglobal, is being developed to explore energetic electron production via magnetic reconnection in macroscale systems. The model is based on the discovery that the production of energetic electrons during reconnection is controlled by Fermi reflection in large-scale magnetic fields and not by parallel electric fields localized in kinetic scale boundary layers. Thus, the m… ▽ More A new computational model, kglobal, is being developed to explore energetic electron production via magnetic reconnection in macroscale systems. The model is based on the discovery that the production of energetic electrons during reconnection is controlled by Fermi reflection in large-scale magnetic fields and not by parallel electric fields localized in kinetic scale boundary layers. Thus, the model eliminates these boundary layers. However, although the parallel electric fields that develop around the magnetic x-line and associated separatrices are not important in producing energetic electrons, there is a large scale electric field that kickstarts the heating of low-energy electrons and drives the cold-electron return current that accompanies escaping energetic electrons in open systems. This macroscale electric field is produced by magnetic-field-aligned gradients in the electron pressure. We have upgraded kglobal to include this large-scale electric field while maintaining energy conservation. The new model is tested by exploring the dynamics of electron acoustic modes which develop as a consequence of the presence of two electron species: hot kinetic and cold fluid electrons. Remarkably, the damping of electron acoustic modes is accurately captured by kglobal. Additionally, it has been established that kglobal correctly describes the dynamics of the interaction of the parallel electric field with escaping hot electrons through benchmarking simulations with the Particle-In-Cell (PIC) code p3d. △ Less

Submitted 18 October, 2019; v1 submitted 17 July, 2019; originally announced July 2019.

Comments: 7 pages, 3 figures, double columns

Journal ref: Physics of Plasmas 26, 102903 (2019)

Instabilities and Turbulence in Low-$β$ Guide Field Reconnection Exhausts with Kinetic Riemann Simulations

Authors: Qile Zhang, J. F. Drake, M. Swisdak

Abstract: The role of turbulence in low-$β$, 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… ▽ More The role of turbulence in low-$β$, 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én 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én speed and is expected to be stable in the low $β$ solar corona. The results suggest that in a guide field reconnection exhaust with $1 \gg β> 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. △ Less

Submitted 11 October, 2019; v1 submitted 6 July, 2019; originally announced July 2019.

Particle Heating and Energy Partition in Low-$β$ Guide Field Reconnection with Kinetic Riemann Simulations

Authors: Qile Zhang, J. F. Drake, M. Swisdak

Abstract: Kinetic Riemann simulations have been completed to explore particle heating during guide field reconnection in the low-$β$ 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… ▽ More Kinetic Riemann simulations have been completed to explore particle heating during guide field reconnection in the low-$β$ 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-$β$ 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. △ Less

Submitted 26 July, 2019; v1 submitted 29 April, 2019; originally announced April 2019.

Comments: Some figures' quality is downgraded due to the size limit of arXiv

Decomposition of Plasma Kinetic Entropy into Position and Velocity Space and the Use of Kinetic Entropy in Particle-in-Cell Simulations

Authors: Haoming Liang, Paul A. Cassak, Sergio Servidio, Michael A. Shay, James F. Drake, Marc Swisdak, Matt R. Argall, John C. Dorelli, Earl E. Scime, William H. Matthaeus, Vadim Roytershteyn, Gian Luca Delzanno

Abstract: We describe a systematic development of kinetic entropy as a diagnostic in fully kinetic particle-in-cell (PIC) simulations and use it to interpret plasma physics processes in heliospheric, planetary, and astrophysical systems. First, we calculate kinetic entropy in two forms -- the ``combinatorial'' form related to the logarithm of the number of microstates per macrostate and the ``continuous'' f… ▽ More We describe a systematic development of kinetic entropy as a diagnostic in fully kinetic particle-in-cell (PIC) simulations and use it to interpret plasma physics processes in heliospheric, planetary, and astrophysical systems. First, we calculate kinetic entropy in two forms -- the ``combinatorial'' form related to the logarithm of the number of microstates per macrostate and the ``continuous'' form related to $f \ln f$, where $f$ is the particle distribution function. We discuss the advantages and disadvantages of each and discuss subtleties about implementing them in PIC codes. Using collisionless PIC simulations that are two-dimensional in position space and three-dimensional in velocity space, we verify the implementation of the kinetic entropy diagnostics and discuss how to optimize numerical parameters to ensure accurate results. We show the total kinetic entropy is conserved to three percent in an optimized simulation of anti-parallel magnetic reconnection. Kinetic entropy can be decomposed into a sum of a position space entropy and a velocity space entropy, and we use this to investigate the nature of kinetic entropy transport during collisionless reconnection. We find the velocity space entropy of both electrons and ions increases in time due to plasma heating during magnetic reconnection, while the position space entropy decreases due to plasma compression. This project uses collisionless simulations, so it cannot address physical dissipation mechanisms; nonetheless, the infrastructure developed here should be useful for studies of collisional or weakly collisional heliospheric, planetary, and astrophysical systems. Beyond reconnection, the diagnostic is expected to be applicable to plasma turbulence and collisionless shocks. △ Less

Submitted 23 August, 2019; v1 submitted 7 February, 2019; originally announced February 2019.

Comments: 20 pages, 7 figures

Journal ref: Physics of Plasmas 26, 082903 (2019)

Transition from ion-coupled to electron-only reconnection: Basic physics and implications for plasma turbulence

Authors: P. Sharma Pyakurel, M. A. Shay, T. D. Phan, W. H. Matthaeus, J. F. Drake, J. M. TenBarge, C. C. Haggerty, K. Klein, P. A. Cassak, T. N. Parashar, M. Swisdak, A. Chasapis

Abstract: 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… ▽ More 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én 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. △ Less

Submitted 27 January, 2019; originally announced January 2019.

A computational model for exploring particle acceleration during reconnection in macro-scale systems

Authors: J. F. Drake, H. Arnold, M. Swisdak, J. T. Dahlin

Abstract: 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… ▽ More 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. △ Less

Submitted 21 December, 2018; v1 submitted 12 September, 2018; originally announced September 2018.

Wave generation and heat flux suppression in astrophysical plasma systems

Authors: G. T. Roberg-Clark, J. F. Drake, M. Swisdak, C. S. Reynolds

Abstract: Heat flux suppression in collisionless plasmas for a large range of plasma $β$ 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-$β$) and double-layer-dominated (low-$β$) heat flux suppression. Whistlers saturate at small amplitude in the low beta limit and are unable to e… ▽ More Heat flux suppression in collisionless plasmas for a large range of plasma $β$ 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-$β$) and double-layer-dominated (low-$β$) 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 $β$ explored. The range of plasma-$β$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. △ Less

Submitted 24 September, 2018; v1 submitted 12 July, 2018; originally announced July 2018.

Comments: 16 pages, 6 figures. Accepted for publication in Astrophys. J

Journal ref: G. T. Roberg-Clark et al 2018 ApJ 867 154

Characterizing Ion Flows Across a Dipolarization Jet

Authors: Harry Arnold, Marc Swisdak, James Drake

Abstract: The structure of dipolarization jets with finite width in the dawn-dusk direction relevant to magnetic reconnection in the Earth's magnetotail is explored with particle-in-cell simulations. We carry out Riemann simulations of the evolution of the jet in the dawn-dusk, north-south plane to investigate the dependence of the jet structure on the jet width in the dawn-dusk direction. We find that the… ▽ More The structure of dipolarization jets with finite width in the dawn-dusk direction relevant to magnetic reconnection in the Earth's magnetotail is explored with particle-in-cell simulations. We carry out Riemann simulations of the evolution of the jet in the dawn-dusk, north-south plane to investigate the dependence of the jet structure on the jet width in the dawn-dusk direction. We find that the magnetic field and Earth-directed ion flow structure depend on the dawn-dusk width. A reversal in the usual Hall magnetic field near the center of the current sheet on the dusk side of larger jets is observed. For small widths, the maximum velocity of the Earthward flow is significantly reduced below the theoretical limit of the upstream Alfvén speed. However, the ion flow speed approaches this limit once the width exceeds the ion Larmor radius based on the normal magnetic field, $B_z$. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Localized Oscillatory Dissipation in Magnetopause Reconnection

Authors: J. L. Burch, R. E. Ergun, P. A. Cassak, J. M. Webster, R. B. Torbert, B. L. Giles, J. C. Dorelli, A. C. Rager, K. -J. Hwang, T. D. Phan, K. J. Genestreti, R. C. Allen, L. -J. Chen, S. Wang, D. Gershman, O. Le Contel, C. T. Russell, R. J. Strangeway, F. D. Wilder, D. B. Graham, M. Hesse, J. F. Drake, M. Swisdak, L. M. Price, M. A. Shay , et al. (4 additional authors not shown)

Abstract: 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 ele… ▽ More 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). △ Less

Submitted 13 December, 2017; originally announced December 2017.

Localized and intense energy conversion in the diffusion region of asymmetric magnetic reconnection

Authors: M. Swisdak, J. F. Drake, L. Price, J. L. Burch, P. A. Cassak, T. -D. Phan

Abstract: 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… ▽ More 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. △ Less

Submitted 14 May, 2018; v1 submitted 12 October, 2017; originally announced October 2017.

Comments: Accepted version

Suppression of electron thermal conduction by whistler turbulence in a sustained thermal gradient

Authors: G. T. Roberg-Clark, J. F. Drake, C. S. Reynolds, M. Swisdak

Abstract: 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, $δB / B_{0} \simeq 1$, and resonantl… ▽ More 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, $δ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 $β$ astrophysical plasmas such as hot accretion flows and the intracluster medium of galaxy clusters. △ Less

Submitted 19 December, 2017; v1 submitted 31 August, 2017; originally announced September 2017.

Comments: Accepted for publication in Physical Review Letters. 5 pages, 5 figures. Supplementary material: 1 page, 1 figure

Journal ref: Phys. Rev. Lett. 120, 035101 (2018)

The role of three-dimensional transport in driving enhanced electron acceleration during magnetic reconnection

Authors: J. T. Dahlin, J. F. Drake, M. Swisdak

Abstract: 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énic outflows expel energetic electrons into flux ropes where they become trapped and disconnected from acce… ▽ More 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énic 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. △ Less

Submitted 1 June, 2017; originally announced June 2017.

Turbulence in Three-Dimensional Simulations of Magnetopause Reconnection

Authors: L. Price, M. Swisdak, J. F. Drake, J. L. Burch, P. A. Cassak, R. E. Ergun

Abstract: 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 charac… ▽ More 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ρ_e\sim(m_eT_e/m_iT_i)^{0.25}$ with $ρ_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{ρ_eρ_i}$ than the $ρ_e$ or $d_e$ scalings seen in 2-D reconnection simulations, and produces significant anomalous resistivity and viscosity in the electron diffusion region. △ Less

Submitted 3 May, 2018; v1 submitted 5 April, 2017; originally announced April 2017.

Comments: 11 pages, 10 figures

Journal ref: Journal of Geophysical Research: Space Physics, 122, 11086--11099 (2017)

The Twist of the Draped Interstellar Magnetic Field Ahead of the Heliopause: A Magnetic Reconnection Driven Rotational Discontinuity

Authors: M. Opher, J. F. Drake, M. Swisdak, B. Zieger, G. Toth

Abstract: 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 comp… ▽ More 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 ($δ=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 $β$ (the plasma to magnetic pressure) across the nose of HP is much greater than in the flanks because the heliosheath $β$ 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. △ Less

Submitted 20 February, 2017; originally announced February 2017.

Comments: 19 pages, 5 figures, submitted

Journal ref: Astrophysical Journal Letters 839, Number 1, 2017

The formation of magnetic depletions and flux annihilation due to reconnection in the heliosheath

Authors: J. F. Drake, M. Swisdak, M. Opher, J. D. Richardson

Abstract: 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, sp… ▽ More 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}. △ Less

Submitted 6 February, 2017; originally announced February 2017.

Particle-in-cell simulation study of the scaling of asymmetric magnetic reconnection with in-plane flow shear

Authors: C. E. Doss, P. A. Cassak, M. Swisdak

Abstract: 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 reco… ▽ More 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. △ Less

Submitted 19 July, 2016; originally announced July 2016.

Comments: 17 pages, 4 figures, submitted to Physics of Plasmas

Parallel electric fields are inefficient drivers of energetic electrons in magnetic reconnection

Authors: J. T. Dahlin, J. F. Drake, M. Swisdak

Abstract: 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 thi… ▽ More 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. △ Less

Submitted 20 December, 2016; v1 submitted 13 July, 2016; originally announced July 2016.

Suppression of electron thermal conduction in the high $β$ intracluster medium of galaxy clusters

Authors: G. T. Roberg-Clark, J. F. Drake, C. S. Reynolds, M. Swisdak

Abstract: 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… ▽ More 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. △ Less

Submitted 19 September, 2016; v1 submitted 14 June, 2016; originally announced June 2016.

Comments: 8 pages, 4 figures. To be published in ApJ Letters

The Effects of Turbulence on Three-Dimensional Magnetic Reconnection at the Magnetopause

Authors: L. Price, M. Swisdak, J. F. Drake, P. A. Cassak, J. T. Dahlin, R. E. Ergun

Abstract: 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 m… ▽ More 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. △ Less

Submitted 7 June, 2016; v1 submitted 27 April, 2016; originally announced April 2016.

Comments: Revised version accepted by GRL

Kinetic signatures of the region surrounding the X-line in asymmetric (magnetopause) reconnection

Authors: M. A. Shay, T. D. Phan, C. C. Haggerty, M. Fujimoto, J. F. Drake, K. Malakit, P. A. Cassak, M. Swisdak

Abstract: 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-lin… ▽ More 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. △ Less

Submitted 1 February, 2016; originally announced February 2016.

Comments: Submitted to Geophysical Research Letters

Quantifying Gyrotropy in Magnetic Reconnection

Authors: M. Swisdak

Abstract: 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-… ▽ More 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. △ Less

Submitted 17 December, 2015; v1 submitted 2 September, 2015; originally announced September 2015.

Comments: Under review. Comments welcome; Now accepted. New version incorporates reviewer comments

A model of the heliosphere with jets

Authors: J. F. Drake, M. Swisdak, M. Opher

Abstract: 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 pressur… ▽ More 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, $β=8π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. △ Less

Submitted 1 July, 2015; v1 submitted 6 May, 2015; originally announced May 2015.

Electron acceleration in three-dimensional magnetic reconnection with a guide field

Authors: J. T. Dahlin, J. F. Drake, M. Swisdak

Abstract: 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 acce… ▽ More 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. △ Less

Submitted 16 November, 2015; v1 submitted 7 March, 2015; originally announced March 2015.

Comments: 13 pages, 4 figures

Fast Magnetic Reconnection Due to Anisotropic Electron Pressure

Authors: P. A. Cassak, R. N. Baylor, R. L. Fermo, M. T. Beidler, M. A. Shay, M. Swisdak, J. F. Drake, H. Karimabadi

Abstract: 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 simulation… ▽ More 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. △ Less

Submitted 5 February, 2015; originally announced February 2015.

Comments: 5 pages, 2 figures, accepted for publication as a letter in Physics of Plasmas

Electron Heating During Magnetic Reconnection: A Simulation Scaling Study

Authors: M. A. Shay, C. C. Haggerty, T. D. Phan, J. F. Drake, P. A. Cassak, P. Wu, M. Oieroset, M. Swisdak, K. Malakit

Abstract: 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én speed $c_{Ar}$ based on the reconnecting magnetic field through the relation $ΔT_e = 0.033 \,m_i\,c_{Ar}^2$, where $ΔT_{e}$ is the increase in electron temperature. For the range… ▽ More 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én speed $c_{Ar}$ based on the reconnecting magnetic field through the relation $ΔT_e = 0.033 \,m_i\,c_{Ar}^2$, where $Δ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 $β$. 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 $ΔT_e$ was proportional to the inflowing Alfvén 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. △ Less

Submitted 5 October, 2014; originally announced October 2014.

Comments: 27 pages, 1 table, 7 figures. Submitted to Physics of Plasmas

Dynamics of double layers, ion acceleration and heat flux suppression during solar flares

Authors: T. C. Li, J. F. Drake, M. Swisdak

Abstract: 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 poten… ▽ More 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. △ Less

Submitted 7 August, 2014; originally announced August 2014.

Journal ref: ApJ, 793, 7 (2014)

The Mechanisms of Electron Heating and Acceleration during Magnetic Reconnection

Authors: J. T. Dahlin, J. F. Drake, M. Swisdak

Abstract: 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… ▽ More 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. △ Less

Submitted 4 June, 2014; v1 submitted 3 June, 2014; originally announced June 2014.

Comments: Submitted to Physics of Plasmas (19 pages, 13 figures)

The Onset of Ion Heating During Magnetic Reconnection with a Strong Guide Field

Authors: J. F. Drake, M. Swisdak

Abstract: 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 $α$ particles. Heating parallel to the local magnetic field during reconnection with a guide field is strongly reduced compared with the reconnection of anti-parallel… ▽ More 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 $α$ 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. △ Less

Submitted 30 April, 2014; originally announced April 2014.

Comments: arXiv admin note: text overlap with arXiv:1108.5750