-
Scale Filtering Analysis of Kinetic Reconnection and its Associated Turbulence
Authors:
Subash Adhikari,
Yan Yang,
William H. Matthaeus,
Paul A. Cassak,
Tulasi N. Parashar,
Michael A. Shay
Abstract:
Previously, using an incompressible von Kármán-Howarth formalism, the behavior of cross-scale energy transfer in magnetic reconnection and turbulence was found to be essentially identical to each other, independent of an external magnetic (guide) field, in the inertial and energy-containing ranges (Adhikari et al., Phys. Plasmas 30, 082904, 2023). However, this description did not account for the…
▽ More
Previously, using an incompressible von Kármán-Howarth formalism, the behavior of cross-scale energy transfer in magnetic reconnection and turbulence was found to be essentially identical to each other, independent of an external magnetic (guide) field, in the inertial and energy-containing ranges (Adhikari et al., Phys. Plasmas 30, 082904, 2023). However, this description did not account for the energy transfer in the dissipation range for kinetic plasmas. In this letter, we adopt a scale-filtering approach to investigate this previously unaccounted-for energy transfer channel in reconnection. Using kinetic particle-in-cell (PIC) simulations of antiparallel and component reconnection, we show that the pressure-strain (PS) interaction becomes important at scales smaller than the ion inertial length, where the nonlinear energy transfer term drops off. Also, the presence of a guide field makes a significant difference in the morphology of the scale-filtered energy transfer. These results are consistent with kinetic turbulence simulations, suggesting that the pressure strain interaction is the dominant energy transfer channel between electron scales and ion scales.
△ Less
Submitted 25 October, 2023;
originally announced October 2023.
-
Effective Viscosity, Resistivity, and Reynolds Number in Weakly Collisional Plasma Turbulence
Authors:
Yan Yang,
William H. Matthaeus,
Sean Oughton,
Riddhi Bandyopadhyay,
Francesco Pecora,
Tulasi N. Parashar,
Vadim Roytershteyn,
Alexandros Chasapis,
Michael A. Shay
Abstract:
We examine dissipation and energy conversion in weakly collisional plasma turbulence, employing in situ observations from the Magnetospheric Multiscale (MMS) mission and kinetic Particle-in-Cell (PIC) simulations of proton-electron plasma. A previous result indicated the presence of viscous-like and resistive-like scaling of average energy conversion rates -- analogous to scalings characteristic o…
▽ More
We examine dissipation and energy conversion in weakly collisional plasma turbulence, employing in situ observations from the Magnetospheric Multiscale (MMS) mission and kinetic Particle-in-Cell (PIC) simulations of proton-electron plasma. A previous result indicated the presence of viscous-like and resistive-like scaling of average energy conversion rates -- analogous to scalings characteristic of collisional systems. This allows for extraction of collisional-like coefficients of effective viscosity and resistivity, and thus also determination of effective Reynolds numbers based on these coefficients. The effective Reynolds number, as a measure of the available bandwidth for turbulence to populate various scales, links macro turbulence properties with kinetic plasma properties in a novel way.
△ Less
Submitted 5 September, 2023;
originally announced September 2023.
-
Statistics of Pressure Fluctuations in Turbulent Kinetic Plasmas
Authors:
Subash Adhikari,
William H. Matthaeus,
Tulasi N. Parashar,
Michael A. Shay,
Paul A. Cassak
Abstract:
In this study we explore the statistics of pressure fluctuations in kinetic collisionless turbulence. A 2.5D kinetic particle-in-cell (PIC) simulation of decaying turbulence is used to investigate pressure balance via the evolution of thermal and magnetic pressure in a plasma with beta of order unity. We also discuss the behavior of thermal, magnetic and total pressure structure functions and thei…
▽ More
In this study we explore the statistics of pressure fluctuations in kinetic collisionless turbulence. A 2.5D kinetic particle-in-cell (PIC) simulation of decaying turbulence is used to investigate pressure balance via the evolution of thermal and magnetic pressure in a plasma with beta of order unity. We also discuss the behavior of thermal, magnetic and total pressure structure functions and their corresponding wavenumber spectra. The total pressure spectrum exhibits a slope of -7/3 extending for about a decade in the ion-inertial range. In contrast, shallower -5/3 spectra are characteristic of the magnetic pressure and thermal pressure. The steeper total pressure spectrum is a consequence of cancellation caused by density-magnetic field magnitude anticorrelation. Further, we evaluate higher order total pressure structure functions in an effort to discuss intermittency and compare the power exponents with higher order structure functions of velocity and magnetic fluctuations. Finally, applications to astrophysical systems are also discussed.
△ Less
Submitted 3 May, 2023;
originally announced May 2023.
-
Three-Dimensional Magnetic Reconnection Spreading in Current Sheets of Non-Uniform Thickness
Authors:
Milton Arencibia,
P. A. Cassak,
M. A. Shay,
Jiong Qiu,
Steven M. Petrinec,
Haoming Liang
Abstract:
Magnetic reconnection in naturally occurring and laboratory settings often begins locally and elongates, or spreads, in the direction perpendicular to the plane of reconnection. Previous work has largely focused on current sheets with a uniform thickness, for which the predicted spreading speed for anti-parallel reconnection is the local speed of the current carriers. We derive a scaling theory of…
▽ More
Magnetic reconnection in naturally occurring and laboratory settings often begins locally and elongates, or spreads, in the direction perpendicular to the plane of reconnection. Previous work has largely focused on current sheets with a uniform thickness, for which the predicted spreading speed for anti-parallel reconnection is the local speed of the current carriers. We derive a scaling theory of three-dimensional (3D) spreading of collisionless anti-parallel reconnection in a current sheet with its thickness varying in the out-of-plane direction, both for spreading from a thinner to thicker region and a thicker to thinner region. We derive an expression for calculating the time it takes for spreading to occur for a current sheet with a given profile of its thickness. A key result is that when reconnection spreads from a thinner to a thicker region, the spreading speed in the thicker region is slower than both the Alfvén speed and the speed of the local current carriers by a factor of the ratio of thin to thick current sheet thicknesses. This is important because magnetospheric and solar observations have previously measured the spreading speed to be slower than previously predicted, so the present mechanism might explain this feature. We confirm the theory via a parametric study using 3D two-fluid numerical simulations. We use the prediction to calculate the time scale for reconnection spreading in Earth's magnetotail during geomagnetic activity. The results are also potentially important for understanding reconnection spreading in solar flares and the dayside magnetopause of Earth and other planets.
△ Less
Submitted 3 March, 2023;
originally announced March 2023.
-
Effect of a guide field on the turbulence like properties of magnetic reconnection
Authors:
Subash Adhikari,
Michael A. Shay,
Tulasi N. Parashar,
William H. Matthaeus,
Prayash S. Pyakurel,
Julia E. Stawarz,
Jonathan P. Eastwood
Abstract:
The effect of an external guide field on the turbulence-like properties of magnetic reconnection is studied using five different 2.5D kinetic particle-in-cell (PIC) simulations. The magnetic energy spectrum is found to exhibit a slope of approximately -5/3 in the inertial range, independent of the guide field. On the contrary, the electric field spectrum, in the inertial range steepens more with t…
▽ More
The effect of an external guide field on the turbulence-like properties of magnetic reconnection is studied using five different 2.5D kinetic particle-in-cell (PIC) simulations. The magnetic energy spectrum is found to exhibit a slope of approximately -5/3 in the inertial range, independent of the guide field. On the contrary, the electric field spectrum, in the inertial range steepens more with the guide field and approaches a slope of -5/3. In addition, spectral analysis of the different terms of the generalized Ohm's law is performed and found to be consistent with PIC simulations of turbulence and MMS observations. Finally, guide field effect on the energy transfer behavior is examined using von-Kármán Howarth (vKH) equation based on incompressible Hall-MHD. The general characteristics of the vKH equation with constant rate of energy transfer in the inertial range, is consistent in all the simulations. This suggests that the qualitative behavior of energy spectrum, and energy transfer in reconnection is similar to that of turbulence, indicating that reconnection fundamentally involves an energy cascade.
△ Less
Submitted 19 February, 2023;
originally announced February 2023.
-
Strategies for determining the cascade rate in MHD turbulence: isotropy, anisotropy, and spacecraft sampling
Authors:
Yanwen Wang,
Rohit Chhiber,
Subash Adhikari,
Yan Yang,
Riddhi Bandyopadhyay,
Michael A. Shay,
Sean Oughton,
William H. Matthaeus,
Manuel E. Cuesta
Abstract:
``Exact'' laws for evaluating cascade rates, tracing back to the Kolmogorov ``4/5'' law, have been extended to many systems of interest including magnetohydrodynamics (MHD), and compressible flows of the magnetofluid and ordinary fluid types. It is understood that implementations may be limited by the quantity of available data and by the lack of turbulence symmetry. Assessment of the accuracy and…
▽ More
``Exact'' laws for evaluating cascade rates, tracing back to the Kolmogorov ``4/5'' law, have been extended to many systems of interest including magnetohydrodynamics (MHD), and compressible flows of the magnetofluid and ordinary fluid types. It is understood that implementations may be limited by the quantity of available data and by the lack of turbulence symmetry. Assessment of the accuracy and feasibility of such ``third-order'' (or Yaglom) relations is most effectively accomplished by examining the von Karman-Howarth equation in increment form, a framework from which the third-order laws are derived as asymptotic approximations. Using this approach, we examine the context of third-order laws for incompressible MHD in some detail. The simplest versions rely on the assumption of isotropy and the presence of a well-defined inertial range, while related procedures generalize the same idea to arbitrary rotational symmetries. Conditions for obtaining correct and accurate values of the dissipation rate from these laws based on several sampling and fitting strategies are investigated using results from simulations. The questions we address are of particular relevance to sampling of solar wind turbulence by one or more spacecraft.
△ Less
Submitted 2 September, 2022; v1 submitted 31 August, 2022;
originally announced September 2022.
-
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.
-
Energy Dissipation in Turbulent Reconnection
Authors:
R. Bandyopadhyay,
A. Chasapis,
W. H. Matthaeus,
T. N. Parashar,
C. C. Haggerty,
M. A. Shay,
D. J. Gershman,
B. L. Giles,
J. L. Burch
Abstract:
We study the nature of pressure-strain interaction at reconnection sites, detected by NASA's Magnetospheric Multiscale (MMS) Mission. We employ data from a series of published case studies, including a large-scale reconnection event at the magnetopause, three small-scale reconnection events at the magnetosheath current sheets, and one example of the recently discovered electron-only reconnection.…
▽ More
We study the nature of pressure-strain interaction at reconnection sites, detected by NASA's Magnetospheric Multiscale (MMS) Mission. We employ data from a series of published case studies, including a large-scale reconnection event at the magnetopause, three small-scale reconnection events at the magnetosheath current sheets, and one example of the recently discovered electron-only reconnection. In all instances, we find that the pressure-strain shows signature of conversion into (or from) internal energy at the reconnection site. The electron heating rate is larger than the ion heating rate and the compressive heating is dominant over the incompressive heating rate in all cases considered. The magnitude of thermal energy conversion rate is close to the electromagnetic energy conversion rate in the reconnection region. Although in most cases the pressure-strain interaction indicates that the particle internal energy is increasing, in one case the internal energy is decreasing. These observations indicate that the pressure-strain interaction can be used as an independent measure of energy conversion and dynamics in reconnection regions, in particular independent of measures based on the electromagnetic work. Finally, we explore a selected reconnection site in a turbulent Particle-in-Cell (PIC) simulation which further supports the observational results.
△ Less
Submitted 4 November, 2021;
originally announced November 2021.
-
Scaling Theory of 3D Magnetic Reconnection Spreading
Authors:
Milton Arencibia,
Paul A. Cassak,
Michael A. Shay,
Eric R. Priest
Abstract:
We develop a first-principles scaling theory of the spreading of three-dimensional (3D) magnetic reconnection of finite extent in the out of plane direction. This theory addresses systems with or without an out of plane (guide) magnetic field, and with or without Hall physics. The theory reproduces known spreading speeds and directions with and without guide fields, unifying previous knowledge in…
▽ More
We develop a first-principles scaling theory of the spreading of three-dimensional (3D) magnetic reconnection of finite extent in the out of plane direction. This theory addresses systems with or without an out of plane (guide) magnetic field, and with or without Hall physics. The theory reproduces known spreading speeds and directions with and without guide fields, unifying previous knowledge in a single theory. New results include: (1) Reconnection spreads in a particular direction if an x-line is induced at the interface between reconnecting and non-reconnecting regions, which is controlled by the out of plane gradient of the electric field in the outflow direction. (2) The spreading mechanism for anti-parallel collisionless reconnection is convection, as is known, but for guide field reconnection it is magnetic field bending. We confirm the theory using 3D two-fluid and resistive-magnetohydrodynamics simulations. (3) The theory explains why anti-parallel reconnection in resistive-magnetohydrodynamics does not spread. (4) The simulation domain aspect ratio, associated with the free magnetic energy, influences whether reconnection spreads or convects with a fixed x-line length. (5) We perform a simulation initiating anti-parallel collisionless reconnection with a pressure pulse instead of a magnetic perturbation, finding spreading is unchanged rather than spreading at the magnetosonic speed as previously suggested. The results provide a theoretical framework for understanding spreading beyond systems studied here, and are important for applications including two-ribbon solar flares and reconnection in Earth's magnetosphere.
△ Less
Submitted 31 July, 2021;
originally announced August 2021.
-
Magnetic reconnection as an energy cascade process
Authors:
S. Adhikari,
T. N. Parashar,
M. A. Shay,
W. H. Matthaeus,
P. Sharma Pyakurel,
S. Fordin,
J. E. Stawarz,
J. P. Eastwood
Abstract:
Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrody…
▽ More
Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady cascade rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. These results provide explicit evidence that reconnection itself is fundamentally an energy cascade process.
△ Less
Submitted 24 April, 2021;
originally announced April 2021.
-
Parker Solar Probe In-Situ Observations of Magnetic Reconnection Exhausts During Encounter 1
Authors:
T. D. Phan,
S. D. Bale,
J. P. Eastwood,
B. Lavraud,
J. F. Drake,
M. Oieroset,
M. A. Shay,
M. Pulupa,
M. Stevens,
R. J. MacDowall,
A. W. Case,
D. Larson,
J. Kasper,
P. Whittlesey,
A. Szabo,
K. E. Korreck,
J. W. Bonnell,
T. Dudok de Wit,
K. Goetz,
P. R. Harvey,
T. S. Horbury,
R. Livi,
D. Malaspina,
K. Paulson,
N. E. Raouafi
, et al. (1 additional authors not shown)
Abstract:
Magnetic reconnection in current sheets converts magnetic energy into particle energy. The process may play an important role in the acceleration and heating of the solar wind close to the Sun. Observations from Parker Solar Probe provide a new opportunity to study this problem, as it measures the solar wind at unprecedented close distances to the Sun. During the 1st orbit, PSP encountered a large…
▽ More
Magnetic reconnection in current sheets converts magnetic energy into particle energy. The process may play an important role in the acceleration and heating of the solar wind close to the Sun. Observations from Parker Solar Probe provide a new opportunity to study this problem, as it measures the solar wind at unprecedented close distances to the Sun. During the 1st orbit, PSP encountered a large number of current sheets in the solar wind through perihelion at 35.7 solar radii. We performed a comprehensive survey of these current sheets and found evidence for 21 reconnection exhausts. These exhausts were observed in heliospheric current sheets, coronal mass ejections, and regular solar wind. However, we find that the majority of current sheets encountered around perihelion, where the magnetic field was strongest and plasma beta was lowest, were Alfvénic structures associated with bursty radial jets and these current sheets did not appear to be undergoing local reconnection. We examined conditions around current sheets to address why some current sheets reconnected, while others did not. A key difference appears to be the degree of plasma velocity shear across the current sheets: The median velocity shear for the 21 reconnection exhausts was 24% of the Alfvén velocity shear, whereas the median shear across 43 Alfvénic current sheets examined was 71% of the Alfvén velocity shear. This finding could suggest that large, albeit sub-Alfvénic, velocity shears suppress reconnection. An alternative interpretation is that the Alfvénic current sheets are isolated rotational discontinuities which do not undergo local reconnection.
△ Less
Submitted 16 January, 2020;
originally announced January 2020.
-
Reconnection from a Turbulence Perspective
Authors:
S. Adhikari,
M. A. Shay,
T. N. Parashar,
P. Sharma Pyakurel,
W. H. Matthaeus,
D. Godzieba,
J. E. Stawarz,
J. P. Eastwood,
J. T. Dahlin
Abstract:
The spectral properties associated with laminar, anti-parallel reconnection are examined using a 2.5D kinetic particle in cell (PIC) simulation. Both the reconnection rate and the energy spectrum exhibit three distinct phases: an initiation phase where the reconnection rate grows, a quasi-steady phase, and a declining phase where both the reconnection rate and the energy spectrum decrease. During…
▽ More
The spectral properties associated with laminar, anti-parallel reconnection are examined using a 2.5D kinetic particle in cell (PIC) simulation. Both the reconnection rate and the energy spectrum exhibit three distinct phases: an initiation phase where the reconnection rate grows, a quasi-steady phase, and a declining phase where both the reconnection rate and the energy spectrum decrease. During the steady phase, the energy spectrum exhibits approximately a double power-law behavior, with a slope near -5/3 at wavenumbers smaller than the inverse ion inertial length, and a slope steeper than -8/3 for larger wavenumbers up to the inverse electron inertial length. This behavior is consistent with a Kolmogorov energy cascade and implies that laminar reconnection may fundamentally be an energy cascade process. Consistent with this idea is that the reconnection rate exhibits a rough correlation with the energy spectrum at wave numbers near the inverse ion inertial length. The 2D spectrum is strongly anisotropic with most energy associated with the wave vector direction normal to the current sheet. Reconnection acts to isotropize the energy spectrum, reducing the Shebalin angle from an initial value of 70 degrees to about 48 degrees (nearly isotropic) by the end of the simulation.
△ Less
Submitted 16 March, 2020; v1 submitted 14 November, 2019;
originally announced November 2019.
-
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.
-
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.
-
Electron-Scale Dynamics of the Diffusion Region during Symmetric Magnetic Reconnection in Space
Authors:
R. B. Torbert,
J. L. Burch,
T. D. Phan,
M. Hesse,
M. R. Argall,
J. Shuster,
R. E. Ergun,
L. Alm,
R. Nakamura,
K. Genestreti,
D. J. Gershman,
W. R. Paterson,
D. L. Turner,
I. Cohen,
B. L. Giles,
C. J. Pollock,
S. Wang,
L. -J. Chen,
Julia Stawarz,
J. P. Eastwood,
K. - J. Hwang,
C. Farrugia,
I. Dors,
H. Vaith,
C. Mouikis
, et al. (24 additional authors not shown)
Abstract:
Magnetic reconnection is an energy conversion process important in many astrophysical contexts including the Earth's magnetosphere, where the process can be investigated in-situ. Here we present the first encounter of a reconnection site by NASA's Magnetospheric Multiscale (MMS) spacecraft in the magnetotail, where reconnection involves symmetric inflow conditions. The unprecedented electron-scale…
▽ More
Magnetic reconnection is an energy conversion process important in many astrophysical contexts including the Earth's magnetosphere, where the process can be investigated in-situ. Here we present the first encounter of a reconnection site by NASA's Magnetospheric Multiscale (MMS) spacecraft in the magnetotail, where reconnection involves symmetric inflow conditions. The unprecedented electron-scale plasma measurements revealed (1) super-Alfvenic electron jets reaching 20,000 km/s, (2) electron meandering motion and acceleration by the electric field, producing multiple crescent-shaped structures, (3) spatial dimensions of the electron diffusion region implying a reconnection rate of 0.1-0.2. The well-structured multiple layers of electron populations indicate that, despite the presence of turbulence near the reconnection site, the key electron dynamics appears to be largely laminar.
△ Less
Submitted 18 September, 2018;
originally announced September 2018.
-
Kinetic range spectral features of cross-helicity using MMS
Authors:
Tulasi N. Parashar,
Alexandros Chasapis,
Riddhi Bandyopadhyay,
Rohit Chhiber,
W. H. Matthaeus,
B. Maruca,
M. A. Shay,
J. L. Burch,
T. E. Moore,
B. L. Giles,
D. J. Gershman,
C. J. Pollock,
R. B. Torbert,
C. T. Russell,
R. J. Strangeway,
Vadim Roytershteyn
Abstract:
We study spectral features of ion velocity and magnetic field correlations in the solar wind and in the magnetosheath using data from the Magnetospheric Multi-Scale (MMS) spacecraft. High resolution MMS observations enable the study of transition of these correlations between their magnetofluid character at larger scales into the sub-proton kinetic range, previously unstudied in spacecraft data. C…
▽ More
We study spectral features of ion velocity and magnetic field correlations in the solar wind and in the magnetosheath using data from the Magnetospheric Multi-Scale (MMS) spacecraft. High resolution MMS observations enable the study of transition of these correlations between their magnetofluid character at larger scales into the sub-proton kinetic range, previously unstudied in spacecraft data. Cross-helicity, angular alignment and energy partitioning is examined over a suit- able range of scales, employing measurements based on the Taylor frozen-in approximation as well as direct two-spacecraft correlation measurements. The results demonstrate signatures of alignment at large scales. As kinetic scales are approached, the alignment between v and b is destroyed by demagnetization of protons.
△ Less
Submitted 6 September, 2018;
originally announced September 2018.
-
Dependence of Kinetic Plasma Turbulence on Plasma beta
Authors:
Tulasi N. Parashar,
William H. Matthaeus,
Michael A Shay
Abstract:
We study the effects of plasma \b{eta} (ratio of plasma pressure to magnetic pressure) on the evolution of kinetic plasma turbulence using fully kinetic particle-in-cell simulations of decaying turbulence. We find that the plasma \b{eta} systematically affects spectra, measures of intermittency, decay rates of turbulence fluctuations, and partitioning over different channels of energy exchange Mor…
▽ More
We study the effects of plasma \b{eta} (ratio of plasma pressure to magnetic pressure) on the evolution of kinetic plasma turbulence using fully kinetic particle-in-cell simulations of decaying turbulence. We find that the plasma \b{eta} systematically affects spectra, measures of intermittency, decay rates of turbulence fluctuations, and partitioning over different channels of energy exchange More specifically, an increase in plasma \b{eta} leads to greater total heating, with proton heating preferentially more than electrons. Implications for achieving magnetosheath like temperature ratios are discussed.
△ Less
Submitted 30 July, 2018;
originally announced July 2018.
-
MMS Observations of Beta-Dependent Constraints on Ion Temperature-Anisotropy in Earth's Magnetosheath
Authors:
Bennett A. Maruca,
A. Chasapis,
S. P. Gary,
R. Bandyopadhyay,
R. Chhiber,
T. N. Parashar,
W. H. Matthaeus,
M. A. Shay,
J. L. Burch,
T. E. Moore,
C. J. Pollock,
B. J. Giles,
W. R. Paterson,
J. Dorelli,
D. J. Gershman,
R. B. Torbert,
C. T. Russell,
R. J. Strangeway
Abstract:
Protons (ionized hydrogen) in the solar wind frequently exhibit distinct temperatures ($T_{\perp p}$ and $T_{\parallel p}$) perpendicular and parallel to the plasma's background magnetic-field. Numerous prior studies of the interplanetary solar-wind have shown that, as plasma beta ($β_{\parallel p}$) increases, a narrower range of temperature-anisotropy (…
▽ More
Protons (ionized hydrogen) in the solar wind frequently exhibit distinct temperatures ($T_{\perp p}$ and $T_{\parallel p}$) perpendicular and parallel to the plasma's background magnetic-field. Numerous prior studies of the interplanetary solar-wind have shown that, as plasma beta ($β_{\parallel p}$) increases, a narrower range of temperature-anisotropy ($R_p\equiv T_{\perp p}\,/\,T_{\parallel p}$) values is observed. Conventionally, this effect has been ascribed to the actions of kinetic microinstabilities. This study is the first to use data from the Magnetospheric Multiscale Mission (MMS) to explore such $β_{\parallel p}$-dependent limits on $R_p$ in Earth's magnetosheath. The distribution of these data across the $(β_{\parallel p},R_p)$-plane reveals limits on both $R_p>1$ and $R_p<1$. Linear Vlasov theory is used to compute contours of constant growth-rate for the ion-cyclotron, mirror, parallel-firehose, and oblique-firehose instabilities. These instability thresholds closely align with the contours of the data distribution, which suggests a strong association of instabilities with extremes of ion temperature anisotropy in the magnetosheath. The potential for instabilities to regulate temperature anisotropy is discussed.
△ Less
Submitted 22 June, 2018;
originally announced June 2018.
-
Incompressive Energy Transfer in the Earth's Magnetosheath: Magnetospheric Multiscale Observations
Authors:
Riddhi Bandyopadhyay,
A. Chasapis,
R. Chhiber,
T. N. Parashar,
W. H. Matthaeus,
M. A. Shay,
B. A. Maruca,
J. L. Burch,
T. E. Moore,
C. J. Pollock,
B. L. Giles,
W. R. Paterson,
J. Dorelli,
D. J. Gershman,
R. B. Torbert,
C. T. Russell,
R. J. Strangeway
Abstract:
Using observational data from the \emph{Magnetospheric Multiscale} (MMS) Mission in the Earth's magnetosheath, we estimate the energy cascade rate using different techniques within the framework of incompressible magnetohydrodynamic (MHD) turbulence. At the energy containing scale, the energy budget is controlled by the von Kármán decay law. Inertial range cascade is estimated by fitting a linear…
▽ More
Using observational data from the \emph{Magnetospheric Multiscale} (MMS) Mission in the Earth's magnetosheath, we estimate the energy cascade rate using different techniques within the framework of incompressible magnetohydrodynamic (MHD) turbulence. At the energy containing scale, the energy budget is controlled by the von Kármán decay law. Inertial range cascade is estimated by fitting a linear scaling to the mixed third-order structure function. Finally, we use a multi-spacecraft technique to estimate the Kolmogorov-Yaglom-like cascade rate in the kinetic range, well below the ion inertial length scale. We find that the inertial range cascade rate is almost equal to the one predicted by the von Kármán law at the energy containing scale, while the cascade rate evaluated at the kinetic scale is somewhat lower, as anticipated in theory~\citep{Yang2017PoP}. Further, in agreement with a recent study~\citep{Hadid2018PRL}, we find that the incompressive cascade rate in the Earth's magnetosheath is about $1000$ times larger than the cascade rate in the pristine solar wind.
△ Less
Submitted 29 August, 2018; v1 submitted 11 June, 2018;
originally announced June 2018.
-
The Reduction of Magnetic Reconnection Outflow Jets to Sub-Alfvénic Speeds
Authors:
Colby C. Haggerty,
Michael A. Shay,
Alexandros Chasapis,
Tai D. Phan,
James F. Drake,
Kittipat Malakit,
Paul A. Cassak,
Rungployphan Kieokaew
Abstract:
The outflow velocity of jets produced by collisionless magnetic reconnection is shown to be reduced by the ion exhaust temperature in simulations and observations. We derive a scaling relationship for the outflow velocity based on the upstream Alfvén speed and the parallel ion exhaust temperature, which is verified in kinetic simulations and observations. The outflow speed reduction is shown to be…
▽ More
The outflow velocity of jets produced by collisionless magnetic reconnection is shown to be reduced by the ion exhaust temperature in simulations and observations. We derive a scaling relationship for the outflow velocity based on the upstream Alfvén speed and the parallel ion exhaust temperature, which is verified in kinetic simulations and observations. The outflow speed reduction is shown to be due to the firehose instability criterion, and so for large enough guide fields this effect is suppressed and the outflow speed reaches the upstream Alfvén speed based on the reconnecting component of the magnetic field.
△ Less
Submitted 25 April, 2018;
originally announced April 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.
-
On the collisionless asymmetric magnetic reconnection rate
Authors:
Yi-Hsin Liu,
M. Hesse,
P. A. Cassak,
M. A. Shay,
S. Wang,
L. -J. Chen
Abstract:
A prediction of the steady-state reconnection electric field in asymmetric reconnection is obtained by maximizing the reconnection rate as a function of the opening angle made by the upstream magnetic field on the weak magnetic field (magnetosheath) side. The prediction is within a factor of two of the widely examined asymmetric reconnection model [Cassak and Shay, Phys. Plasmas 14, 102114, 2007]…
▽ More
A prediction of the steady-state reconnection electric field in asymmetric reconnection is obtained by maximizing the reconnection rate as a function of the opening angle made by the upstream magnetic field on the weak magnetic field (magnetosheath) side. The prediction is within a factor of two of the widely examined asymmetric reconnection model [Cassak and Shay, Phys. Plasmas 14, 102114, 2007] in the collisionless limit, and they scale the same over a wide parameter regime. The previous model had the effective aspect ratio of the diffusion region as a free parameter, which simulations and observations suggest is on the order of 0.1, but the present model has no free parameters. In conjunction with the symmetric case [Liu et al., Phys. Rev. Lett. 118, 085101, 2017], this work further suggests that this nearly universal number 0.1, essentially the normalized fast reconnection rate, is a geometrical factor arising from maximizing the reconnection rate within magnetohydrodynamic (MHD)-scale constraints.
△ Less
Submitted 17 November, 2017;
originally announced November 2017.
-
A Review of the 0.1 Reconnection Rate Problem
Authors:
P. A. Cassak,
Y. -H. Liu,
M. A. Shay
Abstract:
A long-standing problem in magnetic reconnection is to explain why it tends to proceed at or below a normalized rate of 0.1. This article gives a review of observational and numerical evidence for this rate and discusses recent theoretical work addressing this problem. Some remaining open questions are summarized.
A long-standing problem in magnetic reconnection is to explain why it tends to proceed at or below a normalized rate of 0.1. This article gives a review of observational and numerical evidence for this rate and discusses recent theoretical work addressing this problem. Some remaining open questions are summarized.
△ Less
Submitted 11 August, 2017;
originally announced August 2017.
-
Exploring the Statistics of Magnetic Reconnection X-points in Kinetic Particle-in-Cell (PIC) Turbulence
Authors:
Colby C. Haggerty,
Tulasi N. Parashar,
William H. Matthaeus,
Michael A. Shay,
Yan Yang,
Minping Wan,
Penny Wu,
Sergio Servidio
Abstract:
Magnetic reconnection is a ubiquitous phenomenon in turbulent plasmas. It is an important part of the turbulent dynamics and heating of space and astrophysical plasmas. We examine the statistics of magnetic reconnection using a quantitative local analysis of the magnetic vector potential, previously used in magnetohydrodynamics simulations, and now generalized to fully kinetic PIC simulations. Dif…
▽ More
Magnetic reconnection is a ubiquitous phenomenon in turbulent plasmas. It is an important part of the turbulent dynamics and heating of space and astrophysical plasmas. We examine the statistics of magnetic reconnection using a quantitative local analysis of the magnetic vector potential, previously used in magnetohydrodynamics simulations, and now generalized to fully kinetic PIC simulations. Different ways of reducing the particle noise for analysis purposes including multiple smoothing techniques are explored. We find that a Fourier filter applied at the Debye scale is an optimal choice for analyzing PIC data. Finlay, we find a broader distribution of normalized reconnection rates compared to the MHD limit with rates as large as 0.5 but with an average of approximately 0.1.
△ Less
Submitted 15 June, 2017;
originally announced June 2017.
-
Why does steady-state magnetic reconnection have a maximum local rate of order 0.1?
Authors:
Yi-Hsin Liu,
M. Hesse,
F. Guo,
W. Daughton,
H. Li,
P. A. Cassak,
M. A. Shay
Abstract:
Simulations suggest collisionless steady-state magnetic reconnection of Harris-type current sheets proceeds with a rate of order 0.1, independent of dissipation mechanism. We argue this long-standing puzzle is a result of constraints at the magnetohydrodynamic (MHD) scale. We perform a scaling analysis of the reconnection rate as a function of the opening angle made by the upstream magnetic fields…
▽ More
Simulations suggest collisionless steady-state magnetic reconnection of Harris-type current sheets proceeds with a rate of order 0.1, independent of dissipation mechanism. We argue this long-standing puzzle is a result of constraints at the magnetohydrodynamic (MHD) scale. We perform a scaling analysis of the reconnection rate as a function of the opening angle made by the upstream magnetic fields, finding a maximum reconnection rate close to 0.2. The predictions compare favorably to particle-in-cell simulations of relativistic electron-positron and non-relativistic electron-proton reconnection. The fact that simulated reconnection rates are close to the predicted maximum suggests reconnection proceeds near the most efficient state allowed at the MHD-scale. The rate near the maximum is relatively insensitive to the opening angle, potentially explaining why reconnection has a similar fast rate in differing models.
△ Less
Submitted 23 November, 2016;
originally announced November 2016.
-
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.
-
Reduction of Ion Heating During Magnetic Reconnection by Large-Scale Effective Potentials
Authors:
C. C. Haggerty,
M. A. Shay,
J. F. Drake,
T. D. Phan,
C. T. McHugh
Abstract:
The physical processes that control the partition of released magnetic energy between electrons and ions during reconnection is explored through particle-in-cell simulations and analytical techniques. We demonstrate that the development of a large-scale parallel electric field and its associated potential controls the relative heating of electrons and ions. The potential develops to restrain heate…
▽ More
The physical processes that control the partition of released magnetic energy between electrons and ions during reconnection is explored through particle-in-cell simulations and analytical techniques. We demonstrate that the development of a large-scale parallel electric field and its associated potential controls the relative heating of electrons and ions. The potential develops to restrain heated exhaust electrons and enhances their heating by confining electrons in the region where magnetic energy is released. Simultaneously the potential slows ions entering the exhaust below the Alfvénic speed expected from the traditional counterstreaming picture of ion heating. Unexpectedly, the magnitude of the potential and therefore the relative partition of energy between electrons and ions is not a constant but rather depends on the upstream parameters and specifically the upstream electron normalized temperature (electron beta). These findings suggest that the fraction of magnetic energy converted into the total thermal energy may be independent of upstream parameters.
△ Less
Submitted 2 March, 2016; v1 submitted 8 April, 2015;
originally announced April 2015.
-
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.
-
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.
-
Nonlinear and Linear Timescales near Kinetic Scales in Solar Wind Turbulence
Authors:
W. H. Matthaeus,
S. Oughton,
K. T. Osman,
S. Servidio,
M. Wan,
S. P. Gary,
M. A. Shay,
F. Valentini,
V. Roytershteyn,
H. Karimabadi,
S. C. Chapman
Abstract:
The application of linear kinetic treatments to plasma waves, damping, and instability requires favorable inequalities between the associated linear timescales and timescales for nonlinear (e.g., turbulence) evolution. In the solar wind these two types of timescales may be directly compared using standard Kolmogorov-style analysis and observational data. The estimated local nonlinear magnetohydrod…
▽ More
The application of linear kinetic treatments to plasma waves, damping, and instability requires favorable inequalities between the associated linear timescales and timescales for nonlinear (e.g., turbulence) evolution. In the solar wind these two types of timescales may be directly compared using standard Kolmogorov-style analysis and observational data. The estimated local nonlinear magnetohydrodynamic cascade times, evaluated as relevant kinetic scales are approached, remain slower than the cyclotron period, but comparable to, or faster than, the typical timescales of instabilities, anisotropic waves, and wave damping. The variation with length scale of the turbulence timescales is supported by observations and simulations. On this basis the use of linear theory - which assumes constant parameters to calculate the associated kinetic rates - may be questioned. It is suggested that the product of proton gyrofrequency and nonlinear time at the ion gyroscales provides a simple measure of turbulence influence on proton kinetic behavior.
△ Less
Submitted 25 April, 2014;
originally announced April 2014.
-
A New Electric Field in Asymmetric Magnetic Reconnection
Authors:
Kittipat Malakit,
Michael A. Shay,
Paul A. Cassak,
David Ruffolo
Abstract:
We present a theory and numerical evidence for the existence of a previously unexplored in-plane electric field in collisionless asymmetric magnetic reconnection. This electric field, dubbed the "Larmor electric field," is associated with finite Larmor radius effects and is distinct from the known Hall electric field. Potentially, it could be an important indicator for the upcoming Magnetospheric…
▽ More
We present a theory and numerical evidence for the existence of a previously unexplored in-plane electric field in collisionless asymmetric magnetic reconnection. This electric field, dubbed the "Larmor electric field," is associated with finite Larmor radius effects and is distinct from the known Hall electric field. Potentially, it could be an important indicator for the upcoming Magnetospheric Multiscale (MMS) mission to locate reconnection sites as we expect it to appear on the magnetospheric side, pointing Earthward, at the dayside magnetopause reconnection site.
△ Less
Submitted 11 September, 2013;
originally announced September 2013.
-
On the Cause of Supra-Arcade Downflows in Solar Flares
Authors:
P. A. Cassak,
J. F. Drake,
J. T. Gosling,
T. -D. Phan,
M. A. Shay,
L. S. Shepherd
Abstract:
A model of supra-arcade downflows (SADs), dark low density regions also known as tadpoles that propagate sunward during solar flares, is presented. It is argued that the regions of low density are flow channels carved by sunward-directed outflow jets from reconnection. The solar corona is stratified, so the flare site is populated by a lower density plasma than that in the underlying arcade. As th…
▽ More
A model of supra-arcade downflows (SADs), dark low density regions also known as tadpoles that propagate sunward during solar flares, is presented. It is argued that the regions of low density are flow channels carved by sunward-directed outflow jets from reconnection. The solar corona is stratified, so the flare site is populated by a lower density plasma than that in the underlying arcade. As the jets penetrate the arcade, they carve out regions of depleted plasma density which appear as SADs. The present interpretation differs from previous models in that reconnection is localized in space but not in time. Reconnection is continuous in time to explain why SADs are not filled in from behind as they would if they were caused by isolated descending flux tubes or the wakes behind them due to temporally bursty reconnection. Reconnection is localized in space because outflow jets in standard two-dimensional reconnection models expand in the normal (inflow) direction with distance from the reconnection site, which would not produce thin SADs as seen in observations. On the contrary, outflow jets in spatially localized three-dimensional reconnection with an out-of-plane (guide) magnetic field expand primarily in the out-of-plane direction and remain collimated in the normal direction, which is consistent with observed SADs being thin. Two-dimensional proof-of-principle simulations of reconnection with an out-of-plane (guide) magnetic field confirm the creation of SAD-like depletion regions and the necessity of density stratification. Three-dimensional simulations confirm that localized reconnection remains collimated.
△ Less
Submitted 2 September, 2013; v1 submitted 15 July, 2013;
originally announced July 2013.
-
von Kármán energy decay and heating of protons and electrons in a kinetic turbulent plasma
Authors:
P. Wu,
M. Wan,
W. H. Matthaeus,
M. A. Shay,
M. Swisdak
Abstract:
Decay in time of undriven weakly collisional kinetic plasma turbulence in systems large compared to the ion kinetic scales is investigated using fully electromagnetic particle-in-cell simulations initiated with transverse flow and magnetic disturbances, constant density, and a strong guide field. The observed energy decay is consistent with the von Kármán hypothesis of similarity decay, in a formu…
▽ More
Decay in time of undriven weakly collisional kinetic plasma turbulence in systems large compared to the ion kinetic scales is investigated using fully electromagnetic particle-in-cell simulations initiated with transverse flow and magnetic disturbances, constant density, and a strong guide field. The observed energy decay is consistent with the von Kármán hypothesis of similarity decay, in a formulation adapted to magnetohydrodyamics (MHD). Kinetic dissipation occurs at small scales, but the overall rate is apparently controlled by large scale dynamics. At small turbulence amplitude the electrons are preferentially heated. At larger amplitudes proton heating is the dominant effect. In the solar wind and corona the protons are typically hotter, suggesting that these natural systems are in large amplitude turbulence regime.
△ Less
Submitted 30 August, 2013; v1 submitted 18 June, 2013;
originally announced June 2013.
-
Super-Alfvénic propagation of reconnection signatures and Poynting flux during substorms
Authors:
M. A. Shay,
J. F. Drake,
J. P. Eastwood,
T. D. Phan
Abstract:
The propagation of reconnection signatures and their associated energy are examined using kinetic particle-in-cell simulations and Cluster satellite observations. It is found that the quadrupolar out-of-plane magnetic field near the separatrices is associated with a kinetic Alfvén wave. For magnetotail parameters, the parallel propagation of this wave is super-Alfvénic (V_parallel ~ 1500 - 5500 km…
▽ More
The propagation of reconnection signatures and their associated energy are examined using kinetic particle-in-cell simulations and Cluster satellite observations. It is found that the quadrupolar out-of-plane magnetic field near the separatrices is associated with a kinetic Alfvén wave. For magnetotail parameters, the parallel propagation of this wave is super-Alfvénic (V_parallel ~ 1500 - 5500 km/s) and generates substantial Poynting flux (S ~ 10^-5 - 10^-4 W/m^2) consistent with Cluster observations of magnetic reconnection. This Poynting flux substantially exceeds that due to frozen-in ion bulk outflows and is sufficient to generate white light aurora in the Earth's ionosphere.
△ Less
Submitted 5 April, 2011;
originally announced April 2011.
-
Kinetic dissipation and anisotropic heating in a turbulent collisionless plasma
Authors:
T. N. Parashar,
M. A. Shay,
P. A. Cassak,
W. H. Matthaeus
Abstract:
The kinetic evolution of the Orszag-Tang vortex is studied using collisionless hybrid simulations. In the magnetohydrodynamic regime this vortex leads rapidly to broadband turbulence. Significant differences from MHD arise at small scales, where the fluid scale energy dissipates into heat almost exclusively through the magnetic field because the protons are decoupled from the magnetic field. Alt…
▽ More
The kinetic evolution of the Orszag-Tang vortex is studied using collisionless hybrid simulations. In the magnetohydrodynamic regime this vortex leads rapidly to broadband turbulence. Significant differences from MHD arise at small scales, where the fluid scale energy dissipates into heat almost exclusively through the magnetic field because the protons are decoupled from the magnetic field. Although cyclotron resonance is absent, the protons heat preferentially in the plane perpendicular to the mean field, as in the corona and solar wind. Effective transport coefficients are calculated.
△ Less
Submitted 20 June, 2008; v1 submitted 29 December, 2007;
originally announced January 2008.
-
From Solar and Stellar Flares to Coronal Heating: Theory and Observations of How Magnetic Reconnection Regulates Coronal Conditions
Authors:
P. A. Cassak,
D. J. Mullan,
M. A. Shay
Abstract:
There is currently no explanation of why the corona has the temperature and density it has. We present a model which explains how the dynamics of magnetic reconnection regulates the conditions in the corona. A bifurcation in magnetic reconnection at a critical state enforces an upper bound on the coronal temperature for a given density. We present observational evidence from 107 flares in 37 sun…
▽ More
There is currently no explanation of why the corona has the temperature and density it has. We present a model which explains how the dynamics of magnetic reconnection regulates the conditions in the corona. A bifurcation in magnetic reconnection at a critical state enforces an upper bound on the coronal temperature for a given density. We present observational evidence from 107 flares in 37 sun-like stars that stellar coronae are near this critical state. The model may be important to self-organized criticality models of the solar corona.
△ Less
Submitted 5 February, 2008; v1 submitted 17 October, 2007;
originally announced October 2007.
-
Two-scale structure of the electron dissipation region during collisionless magnetic reconnection
Authors:
M. A. Shay,
J. F. Drake,
M. Swisdak
Abstract:
Particle in cell (PIC) simulations of collisionless magnetic reconnection are presented that demonstrate that the electron dissipation region develops a distinct two-scale structure along the outflow direction. The length of the electron current layer is found to decrease with decreasing electron mass, approaching the ion inertial length for a proton-electron plasma. A surprise, however, is that…
▽ More
Particle in cell (PIC) simulations of collisionless magnetic reconnection are presented that demonstrate that the electron dissipation region develops a distinct two-scale structure along the outflow direction. The length of the electron current layer is found to decrease with decreasing electron mass, approaching the ion inertial length for a proton-electron plasma. A surprise, however, is that the electrons form a high-velocity outflow jet that remains decoupled from the magnetic field and extends large distances downstream from the x-line. The rate of reconnection remains fast in very large systems, independent of boundary conditions and the mass of electrons.
△ Less
Submitted 5 April, 2007;
originally announced April 2007.
-
A Model for Spontaneous Onset of Fast Magnetic Reconnection
Authors:
P. A. Cassak,
J. F. Drake,
M. A. Shay
Abstract:
We present a model for the spontaneous onset of fast magnetic reconnection in a weakly collisional plasma, such as the solar corona. When a current layer of macroscopic width undergoes collisional (Sweet-Parker) reconnection, a narrow dissipation region forms around the X-line. This dissipation region naturally becomes narrower during the reconnection process as stronger magnetic fields are conv…
▽ More
We present a model for the spontaneous onset of fast magnetic reconnection in a weakly collisional plasma, such as the solar corona. When a current layer of macroscopic width undergoes collisional (Sweet-Parker) reconnection, a narrow dissipation region forms around the X-line. This dissipation region naturally becomes narrower during the reconnection process as stronger magnetic fields are convected toward the X-line. When the dissipation region becomes thinner than the ion skin depth, resistive magnetohydrodynamics breaks down as the Hall effect becomes important and the Sweet-Parker solution ceases to exist. A transition to collisionless (Hall) reconnection ensues, increasing the reconnection rate by many orders of magnitude in a very short time. Predictions of the model are consistent with constraints set by observations of solar flares.
△ Less
Submitted 31 March, 2006;
originally announced April 2006.
-
The Transition from Anti-Parallel to Component Magnetic Reconnection
Authors:
M. Swisdak,
J. F. Drake,
M. A. Shay,
J. G. McIlhargey
Abstract:
We study the transition between anti-parallel and component collisionless magnetic reconnection with 2D particle-in-cell simulations. The primary finding is that a guide field \approx 0.1 times as strong as the asymptotic reconnecting field -- roughly the field strength at which the electron Larmor radius is comparable to the width of the electron current layer -- is sufficient to magnetize the…
▽ More
We study the transition between anti-parallel and component collisionless magnetic reconnection with 2D particle-in-cell simulations. The primary finding is that a guide field \approx 0.1 times as strong as the asymptotic reconnecting field -- roughly the field strength at which the electron Larmor radius is comparable to the width of the electron current layer -- is sufficient to magnetize the electrons in the vicinity of the x-line, thus causing significant changes to the structure of the electron dissipation region. This implies that great care should be exercised before concluding that magnetospheric reconnection is antiparallel. We also find that even for such weak guide fields strong inward-flowing electron beams form in the vicinity of the magnetic separatrices and Buneman-unstable distribution functions arise at the x-line itself. As in the calculations of {\it Hesse et al.} [2002] and {\it Yin and Winske} [2003], the non-gyrotropic elements of the electron pressure tensor play the dominant role in decoupling the electrons from the magnetic field at the x-line, regardless of the magnitude of the guide field and the associated strong variations in the pressure tensor's spatial structure. Despite these changes, and consistent with previous work, the reconnection rate does not vary appreciably with the strength of the guide field as it changes between 0 and a value equal to the asymptotic reversed field.
△ Less
Submitted 8 March, 2005;
originally announced March 2005.
-
A catastrophe model for fast magnetic reconnection onset
Authors:
P. A. Cassak,
M. A. Shay,
J. F. Drake
Abstract:
A catastrophe model for the onset of fast magnetic reconnection is presented that suggests why plasma systems with magnetic free energy remain apparently stable for long times and then suddenly release their energy. For a given set of plasma parameters there are generally two stable reconnection solutions: a slow (Sweet-Parker) solution and a fast (Alfvénic) Hall reconnection solution. Below a c…
▽ More
A catastrophe model for the onset of fast magnetic reconnection is presented that suggests why plasma systems with magnetic free energy remain apparently stable for long times and then suddenly release their energy. For a given set of plasma parameters there are generally two stable reconnection solutions: a slow (Sweet-Parker) solution and a fast (Alfvénic) Hall reconnection solution. Below a critical resistivity the slow solution disappears and fast reconnection dominates. Scaling arguments predicting the two solutions and the critical resistivity are confirmed with two-fluid simulations.
△ Less
Submitted 9 June, 2005; v1 submitted 31 January, 2005;
originally announced February 2005.
-
Three species collisionless reconnection: Effect of O+ on magnetotail reconnection
Authors:
M. A. Shay,
M. Swisdak
Abstract:
The nature of collisionless reconnection in a three-species plasma composed of a heavy species, protons, and electrons is examined. Besides the usual two length scales present in two-species reconnection, there are two additional larger length scales in the system: one associated with a "heavy whistler" which produces a large scale quadrupolar out-of-plane magnetic field, and one associated with…
▽ More
The nature of collisionless reconnection in a three-species plasma composed of a heavy species, protons, and electrons is examined. Besides the usual two length scales present in two-species reconnection, there are two additional larger length scales in the system: one associated with a "heavy whistler" which produces a large scale quadrupolar out-of-plane magnetic field, and one associated with the "heavy Alfven" wave which can slow the outflow speed and thus the reconnection rate. The consequences for reconnection in the magnetotail with an O+ population present are discussed.
△ Less
Submitted 13 September, 2004; v1 submitted 15 June, 2004;
originally announced June 2004.
-
Diamagnetic Suppression of Component Magnetic Reconnection at the Magnetopause
Authors:
M. Swisdak,
B. N. Rogers,
J. F. Drake,
M. A. Shay
Abstract:
We present particle-in-cell simulations of collisionless magnetic reconnection in a system (like the magnetopause) with a large density asymmetry across the current layer. In the presence of an ambient component of the magnetic field perpendicular to the reconnection plane the gradient creates a diamagnetic drift that advects the X-line with the electron diamagnetic velocity. When the relative d…
▽ More
We present particle-in-cell simulations of collisionless magnetic reconnection in a system (like the magnetopause) with a large density asymmetry across the current layer. In the presence of an ambient component of the magnetic field perpendicular to the reconnection plane the gradient creates a diamagnetic drift that advects the X-line with the electron diamagnetic velocity. When the relative drift between the ions and electrons is of the order the Alfven speed the large scale outflows from the X-line necessary for fast reconnection cannot develop and the reconnection is suppressed. We discuss how these effects vary with both the plasma beta and the shear angle of the reconnecting field and discuss observational evidence for diamagnetic stabilization at the magnetopause.
△ Less
Submitted 7 May, 2003;
originally announced May 2003.