-
The Importance of Heat Flux in Quasi-Parallel Collisionless Shocks
Authors:
Colby C. Haggerty,
Damiano Caprioli,
Paul A. Cassak,
M. Hasan Barbhuiya,
Lynn Wilson III,
Drew Turner
Abstract:
Collisionless plasma shocks are a common feature of many space and astrophysical systems and are sources of high-energy particles and non-thermal emission, channeling as much as 20\% of the shock's energy into non-thermal particles. The generation and acceleration of these non-thermal particles have been extensively studied, however, how these particles feed back on the shock hydrodynamics has not…
▽ More
Collisionless plasma shocks are a common feature of many space and astrophysical systems and are sources of high-energy particles and non-thermal emission, channeling as much as 20\% of the shock's energy into non-thermal particles. The generation and acceleration of these non-thermal particles have been extensively studied, however, how these particles feed back on the shock hydrodynamics has not been fully treated. This work presents the results of self-consistent hybrid particle-in-cell simulations that show the effect of self-generated non-thermal particle populations on the nature of collisionless, quasi-parallel shocks. They contribute to a significant heat flux density upstream of the shock. Non-thermal particles downstream of the shock leak into the upstream region, taking energy away from the shock. This increases the compression ratio, slows the shock down, and flattens the non-thermal population's spectral index for lower Mach number shocks. We incorporate this into a revised theory for the Rankine-Hugoniot jump conditions that include this effect and it shows excellent agreement with simulations. The results have the potential to explain discrepancies between predictions and observations in a wide range of systems, such as inaccuracies of predictions of arrival times of coronal mass ejections and the conflicting radio and x-ray observations of intracluster shocks.
These effects will likely need to be included in fluid modeling to accurately predict shock evolution.
△ Less
Submitted 23 July, 2023;
originally announced July 2023.
-
Quantifying Energy Conversion in Higher Order Phase Space Density Moments in Plasmas
Authors:
Paul A. Cassak,
M. Hasan Barbhuiya,
Haoming Liang,
Matthew R. Argall
Abstract:
Weakly collisional and collisionless plasmas are typically far from local thermodynamic equilibrium (LTE), and understanding energy conversion in such systems is a forefront research problem. The standard approach is to investigate changes in internal (thermal) energy and density, but this omits energy conversion that changes any higher order moments of the phase space density. In this study, we c…
▽ More
Weakly collisional and collisionless plasmas are typically far from local thermodynamic equilibrium (LTE), and understanding energy conversion in such systems is a forefront research problem. The standard approach is to investigate changes in internal (thermal) energy and density, but this omits energy conversion that changes any higher order moments of the phase space density. In this study, we calculate from first principles the energy conversion associated with all higher moments of the phase space density for systems not in LTE. Particle-in-cell simulations of collisionless magnetic reconnection reveal that energy conversion associated with higher order moments can be locally significant. The results may be useful in numerous plasma settings, such as reconnection, turbulence, shocks, and wave-particle interactions in heliospheric, planetary, and astrophysical plasmas.
△ Less
Submitted 1 June, 2023;
originally announced June 2023.
-
Scaling of electron heating by magnetization during reconnection and applications to dipolarization fronts and super-hot solar flares
Authors:
M. Hasan Barbhuiya,
Paul. A. Cassak,
Michael. A. Shay,
Vadim Roytershteyn,
Marc Swisdak,
Amir Caspi,
Andrei Runov,
Haoming Liang
Abstract:
Electron ring velocity space distributions have previously been seen in numerical simulations of magnetic reconnection exhausts and have been suggested to be caused by the magnetization of the electron outflow jet by the compressed reconnected magnetic fields [Shuster et al., ${\it Geophys.~Res.~Lett.}, {\bf 41}$, 5389 (2014)]. We present a theory of the dependence of the major and minor radii of…
▽ More
Electron ring velocity space distributions have previously been seen in numerical simulations of magnetic reconnection exhausts and have been suggested to be caused by the magnetization of the electron outflow jet by the compressed reconnected magnetic fields [Shuster et al., ${\it Geophys.~Res.~Lett.}, {\bf 41}$, 5389 (2014)]. We present a theory of the dependence of the major and minor radii of the ring distributions solely in terms of upstream (lobe) plasma conditions, thereby allowing a prediction of the associated temperature and temperature anisotropy of the rings in terms of upstream parameters. We test the validity of the prediction using 2.5-dimensional particle-in-cell (PIC) simulations with varying upstream plasma density and temperature, finding excellent agreement between the predicted and simulated values. We confirm the Shuster et al. suggestion for the cause of the ring distributions, and also find that the ring distributions are located in a region marked by a plateau, or shoulder, in the reconnected magnetic field profile. The predictions of the temperature are consistent with observed electron temperatures in dipolarization fronts, and may provide an explanation for the generation of plasma with temperatures in the 10s of MK in super-hot solar flares. A possible extension of the model to dayside reconnection is discussed. Since ring distributions are known to excite whistler waves, the present results should be useful for quantifying the generation of whistler waves in reconnection exhausts.
△ Less
Submitted 31 July, 2022;
originally announced August 2022.
-
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.
-
Dissipation measures in weakly-collisional plasmas
Authors:
O. Pezzi,
H. Liang,
J. L. Juno,
P. A. Cassak,
C. L. Vasconez,
L. Sorriso-Valvo,
D. Perrone,
S. Servidio,
V. Roytershteyn,
J. M. TenBarge,
W. H. Matthaeus
Abstract:
The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. W…
▽ More
The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell VPIC, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov-Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thus exhibiting local correlation. The distribution function-based measures show a broader width compared to energy-based proxies, suggesting that energy transfer is co-localized at coherent structures, but can affect the particle distribution function in wider regions. The effect of interparticle collisions on these parameters is finally discussed.
△ Less
Submitted 26 June, 2021; v1 submitted 3 January, 2021;
originally announced January 2021.
-
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.
-
Elongation of Flare Ribbons
Authors:
Jiong Qiu,
Dana W. Longcope,
Paul A. Cassak,
Eric R. Priest
Abstract:
We present an analysis of the apparent elongation motion of flare ribbons along the polarity inversion line (PIL) as well as the shear of flare loops in several two-ribbon flares. Flare ribbons and loops spread along the PIL at a speed ranging from a few to a hundred km/s. The shear measured from conjugate foot-points is consistent with the measurement from flare loops, and both show the decrease…
▽ More
We present an analysis of the apparent elongation motion of flare ribbons along the polarity inversion line (PIL) as well as the shear of flare loops in several two-ribbon flares. Flare ribbons and loops spread along the PIL at a speed ranging from a few to a hundred km/s. The shear measured from conjugate foot-points is consistent with the measurement from flare loops, and both show the decrease of shear toward a potential field as a flare evolves and ribbons and loops spread along the PIL. Flares exhibiting fast bi-directional elongation appear to have a strong shear, which may indicate a large magnetic guide field relative to the reconnection field in the coronal current sheet. We discuss how the analysis of ribbon motion could help infer properties in the corona where reconnection takes place.
△ Less
Submitted 8 July, 2017;
originally announced July 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.
-
Asymmetric magnetic reconnection with a flow shear and applications to the magnetopause
Authors:
C. E. Doss,
C. M. Komar,
P. A. Cassak,
F. D. Wilder,
S. Eriksson,
J. F. Drake
Abstract:
We perform a theoretical and numerical study of anti-parallel 2D magnetic reconnection with asymmetries in the density and reconnecting magnetic field strength in addition to a bulk flow shear across the reconnection site in the plane of the reconnecting fields, which commonly occurs at planetary magnetospheres. We predict the speed at which an isolated X-line is convected by the flow, the reconne…
▽ More
We perform a theoretical and numerical study of anti-parallel 2D magnetic reconnection with asymmetries in the density and reconnecting magnetic field strength in addition to a bulk flow shear across the reconnection site in the plane of the reconnecting fields, which commonly occurs at planetary magnetospheres. We predict the speed at which an isolated X-line is convected by the flow, the reconnection rate, and the critical flow speed at which reconnection no longer takes place for arbitrary reconnecting magnetic field strengths, densities, and upstream flow speeds, and confirm the results with two-fluid numerical simulations. The predictions and simulation results counter the prevailing model of reconnection at Earth's dayside magnetopause which says reconnection occurs with a stationary X-line for sub-Alfvenic magnetosheath flow, reconnection occurs but the X-line convects for magnetosheath flows between the Alfven speed and double the Alfven speed, and reconnection does not occur for magnetosheath flows greater than double the Alfven speed. We find that X-line motion is governed by momentum conservation from the upstream flows, which are weighted differently in asymmetric systems, so the X-line convects for generic conditions including sub-Alfvenic upstream speeds. For the reconnection rate, while the cutoff condition for symmetric reconnection is that the difference in flows on the two sides of the reconnection site is twice the Alfven speed, we find asymmetries cause the cutoff speed for asymmetric reconnection to be higher than twice the asymmetric form of the Alfven speed. The results compare favorably with an observation of reconnection at Earth's polar cusps during a period of northward interplanetary magnetic field, where reconnection occurs despite the magnetosheath flow speed being more than twice the magnetosheath Alfven speed, the previously proposed suppression condition.
△ Less
Submitted 4 September, 2015;
originally announced September 2015.
-
Comparative Analysis of Dayside Reconnection Models in Global Magnetosphere Simulations
Authors:
C. M. Komar,
R. L. Fermo,
P. A. Cassak
Abstract:
We test and compare a number of existing models predicting the location of magnetic reconnection at Earth's dayside magnetopause for various solar wind conditions. We employ robust image processing techniques to determine the locations where each model predicts reconnection to occur. The predictions are then compared to the magnetic separators, the magnetic field lines separating different magneti…
▽ More
We test and compare a number of existing models predicting the location of magnetic reconnection at Earth's dayside magnetopause for various solar wind conditions. We employ robust image processing techniques to determine the locations where each model predicts reconnection to occur. The predictions are then compared to the magnetic separators, the magnetic field lines separating different magnetic topologies. The predictions are tested in distinct high-resolution simulations with interplanetary magnetic field (IMF) clock angles ranging from 30 to 165 degrees in global magnetohydrodynamic simulations using the three-dimensional Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) code with a uniform resistivity, although the described techniques can be generally applied to any self-consistent magnetosphere code. Additional simulations are carried out to test location model dependence on IMF strength and dipole tilt. We find that most of the models match large portions of the magnetic separators when the IMF has a southward component, with the models saying reconnection occurs where the local reconnection rate and reconnection outflow speed are maximized performing best. When the IMF has a northward component, none of the models tested faithfully map the entire magnetic separator, but the maximum magnetic shear model is the best at mapping the separator in the cusp region where reconnection has been observed. Predictions for some models with northward IMF orientations improve after accounting for plasma flow shear parallel to the reconnecting components of the magnetic fields. Implications for observations are discussed.
△ Less
Submitted 23 April, 2015;
originally announced April 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.
-
On the 3-D structure and dissipation of reconnection-driven flow-bursts
Authors:
J. F. Drake,
M. Swisdak,
P. A. Cassak,
T. D. Phan
Abstract:
The structure of magnetic reconnection-driven outflows and their dissipation are explored with large-scale, 3-D particle-in-cell (PIC) simulations. Outflow jets resulting from 3-D reconnection with a finite length x-line form fronts as they propagate into the downstream medium. A large pressure increase ahead of this ``reconnection jet front'' (RJF), due to reflected and transmitted ions, slows th…
▽ More
The structure of magnetic reconnection-driven outflows and their dissipation are explored with large-scale, 3-D particle-in-cell (PIC) simulations. Outflow jets resulting from 3-D reconnection with a finite length x-line form fronts as they propagate into the downstream medium. A large pressure increase ahead of this ``reconnection jet front'' (RJF), due to reflected and transmitted ions, slows the front so that its velocity is well below the velocity of the ambient ions in the core of the jet. As a result, the RJF slows and diverts the high-speed flow into the direction perpendicular to the reconnection plane. The consequence is that the RJF acts as a thermalization site for the ion bulk flow and contributes significantly to the dissipation of magnetic energy during reconnection even though the outflow jet is subsonic. This behavior has no counterpart in 2-D reconnection. A simple analytic model predicts the front velocity and the fraction of the ion bulk flow energy that is dissipated.
△ Less
Submitted 18 March, 2014; v1 submitted 27 January, 2014;
originally announced January 2014.
-
Tracing magnetic separators and their dependence on IMF clock angle in global magnetospheric simulations
Authors:
C. M. Komar,
P. A. Cassak,
J. C. Dorelli,
A. Glocer,
M. M. Kuznetsova
Abstract:
A new, efficient, and highly accurate method for tracing magnetic separators in global magnetospheric simulations with arbitrary clock angle is presented. The technique is to begin at a magnetic null and iteratively march along the separator by finding where four magnetic topologies meet on a spherical surface. The technique is verified using exact solutions for separators resulting from an analyt…
▽ More
A new, efficient, and highly accurate method for tracing magnetic separators in global magnetospheric simulations with arbitrary clock angle is presented. The technique is to begin at a magnetic null and iteratively march along the separator by finding where four magnetic topologies meet on a spherical surface. The technique is verified using exact solutions for separators resulting from an analytic magnetic field model that superposes dipolar and uniform magnetic fields. Global resistive magnetohydrodynamic simulations are performed using the three-dimensional BATS-R-US code with a uniform resistivity, in eight distinct simulations with interplanetary magnetic field (IMF) clock angles ranging from 0 (parallel) to 180 degrees (anti-parallel). Magnetic nulls and separators are found in the simulations, and it is shown that separators traced here are accurate for any clock angle, unlike the last closed field line on the Sun-Earth line that fails for southward IMF. Trends in magnetic null locations and the structure of magnetic separators as a function of clock angle are presented and compared with those from the analytic field model. There are many qualitative similarities between the two models, but quantitative differences are also noted. Dependence on solar wind density is briefly investigated.
△ Less
Submitted 3 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.
-
Estimates of Densities and Filling Factors from a Cooling Time Analysis of Solar Microflares Observed with RHESSI
Authors:
R. N. Baylor,
P. A. Cassak,
S. Christe,
I. G. Hannah,
Säm Krucker,
D. J. Mullan,
M. A. Shay,
H. S. Hudson,
R. P. Lin
Abstract:
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 tha…
▽ More
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).
△ Less
Submitted 20 July, 2011;
originally announced July 2011.
-
Three-dimensional simulations of the orientation and structure of reconnection X-lines
Authors:
R. Schreier,
M. Swisdak,
J. F. Drake,
P. A. Cassak
Abstract:
This work employs Hall magnetohydrodynamic (MHD) simulations to study the X-lines formed during the reconnection of magnetic fields with differing strengths and orientations embedded in plasmas of differing densities. Although random initial perturbations trigger the growth of X-lines with many orientations, at late time a few robust X-lines sharing an orientation reasonably consistent with the di…
▽ More
This work employs Hall magnetohydrodynamic (MHD) simulations to study the X-lines formed during the reconnection of magnetic fields with differing strengths and orientations embedded in plasmas of differing densities. Although random initial perturbations trigger the growth of X-lines with many orientations, at late time a few robust X-lines sharing an orientation reasonably consistent with the direction that maximizes the outflow speed, as predicted by Swisdak and Drake [Geophys. Res. Lett., 34, L11106, (2007)], dominate the system. The existence of reconnection in the geometry examined here contradicts the suggestion of Sonnerup [J. Geophys. Res., 79, 1546 (1974)] that reconnection occurs in a plane normal to the equilibrium current. At late time the growth of the X-lines stagnates, leaving them shorter than the simulation domain.
△ Less
Submitted 8 September, 2010;
originally announced September 2010.
-
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.