Source author record

M. A. Shay

M. A. Shay appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2025arXiv

Factors Controlling the Statistics of Magnetic Reconnection in MHD Turbulence

We study the statistics of dynamical quantities associated with magnetic reconnection events embedded in a sea of strong background magnetohydrodynamic (MHD) turbulence using direct numerical simulations. We focus on the relationship of the reconnection properties to the statistics of global turbulent fields. For the first time, we show that the distribution in turbulence of reconnection rates (determined by upstream fields) is strongly correlated with the magnitude of the global turbulent magnetic field at the correlation scale. The average reconnection rates, and associated dissipation rates, during turbulence are thus much larger than predicted by using turbulent magnetic field fluctuation amplitudes at the dissipation or kinetic scales. Magnetic reconnection may therefore be playing a major role in energy dissipation in astrophysical and heliospheric turbulence.

preprint2021arXiv

Magnetic reconnection as an energy cascade process

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.

preprint2020arXiv

Parker Solar Probe In-Situ Observations of Magnetic Reconnection Exhausts During Encounter 1

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.

preprint2020arXiv

Reconnection from a Turbulence Perspective

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.

preprint2019arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

Fast Magnetic Reconnection Due to Anisotropic Electron Pressure

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.

preprint2014arXiv

Electron Heating During Magnetic Reconnection: A Simulation Scaling Study

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.

preprint2014arXiv

Nonlinear and Linear Timescales near Kinetic Scales in Solar Wind Turbulence

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.

preprint2013arXiv

On the Cause of Supra-Arcade Downflows in Solar Flares

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.

preprint2013arXiv

von Kármán energy decay and heating of protons and electrons in a kinetic turbulent plasma

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.

preprint2011arXiv

Estimates of Densities and Filling Factors from a Cooling Time Analysis of Solar Microflares Observed with RHESSI

We use more than 4,500 microflares from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) microflare data set (Christe et al., 2008, Ap. J., 677, 1385) to estimate electron densities and volumetric filling factors of microflare loops using a cooling time analysis. We show that if the filling factor is assumed to be unity, the calculated conductive cooling times are much shorter than the observed flare decay times, which in turn are much shorter than the calculated radiative cooling times. This is likely unphysical, but the contradic- tion can be resolved by assuming the radiative and conductive cooling times are comparable, which is valid when the flare loop temperature is a maximum and when external heating can be ignored. We find that resultant radiative and con- ductive cooling times are comparable to observed decay times, which has been used as an assumption in some previous studies. The inferred electron densities have a mean value of 10^11.6 cm^-3 and filling factors have a mean of 10^-3.7. The filling factors are lower and densities are higher than previous estimates for large flares, but are similar to those found for two microflares by Moore et al. (Ap. J., 526, 505, 1999).

preprint2011arXiv

Super-Alfvénic propagation of reconnection signatures and Poynting flux during substorms

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.

preprint2008arXiv

Kinetic dissipation and anisotropic heating in a turbulent collisionless plasma

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.

preprint2005arXiv

The Transition from Anti-Parallel to Component Magnetic Reconnection

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.

preprint2003arXiv

Diamagnetic Suppression of Component Magnetic Reconnection at the Magnetopause

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.