BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

M. Swisdak

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

ResearcherUnclaimed source record
physics.plasm-phphysics.space-phastro-ph.SRastro-ph.HEastro-phastro-ph.EPastro-ph.GAnlin.PSastro-ph.COastro-ph.IMnlin.AOphysics.comp-ph

Catalog footprint

What is connected

42works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

42 published item(s)

preprint2021arXiv

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

Transport equations for electron thermal energy in the high $β_e$ intracluster medium (ICM) are developed that include scattering from both classical collisions and self-generated whistler waves. The calculation employs an expansion of the kinetic electron equation along the ambient magnetic field in the limit of strong scattering and assumes whistler waves with low phase speeds $V_w\sim{v}_{te}/β_e\ll{v}_{te}$ dominate the turbulent spectrum, with $v_{te}$ the electron thermal speed and $β_e\gg1$ the ratio of electron thermal to magnetic pressure. We find: (1) temperature-gradient-driven whistlers dominate classical scattering when $L_c>L/β_e$, with $L_c$ the classical electron mean-free-path and $L$ the electron temperature scale length, and (2) in the whistler dominated regime the electron thermal flux is controlled by both advection at $V_w$ and a comparable diffusive term. The findings suggest whistlers limit electron heat flux over large regions of the ICM, including locations unstable to isobaric condensation. Consequences include: (1) the Field length decreases, extending the domain of thermal instability to smaller length-scales, (2) the heat flux temperature dependence changes from $T_e^{7/2}/L$ to $V_wnT_e\sim{T}_e^{1/2}$, (3) the magneto-thermal and heat-flux driven buoyancy instabilities are impaired or completely inhibited, and (4) sound waves in the ICM propagate greater distances, as inferred from observations. This description of thermal transport can be used in macroscale ICM models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Noise-Induced Magnetic Field Saturation in Kinetic Simulations

Monte Carlo methods are often employed to numerically integrate kinetic equations, such as the particle-in-cell method for the plasma kinetic equation, but these methods suffer from the introduction of counting noise to the solution. We report on a cautionary tale of counting noise modifying the nonlinear saturation of kinetic instabilities driven by unstable beams of plasma. We find a saturated magnetic field in under-resolved particle-in-cell simulations due to the sampling error in the current density. The noise-induced magnetic field is anomalous, as the magnetic field damps away in continuum kinetic and increased particle count particle-in-cell simulations. This modification of the saturated state has implications for a broad array of astrophysical phenomena beyond the simple plasma system considered here, and it stresses the care that must be taken when using particle methods for kinetic equations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

The role of turbulence in low-$β$, guide-field reconnection exhausts is explored in 2D reconnection and 2D and 3D Riemann simulations. The structure of the exhaust and associated turbulence is controlled by a pair of rotational discontinuities (RDs) at the exhaust boundary and a pair of slow shocks (SSs) that are generated by counterstreaming ions beams. In 2D the exhaust develops large-amplitude striations at the ion Larmor radius scale that are produced by electron-beam-driven ion cyclotron waves. The electron beams driving the instability are injected into the exhaust from one of the RDs. However, in 3D Riemann simulations, the additional dimension (in the out-of-plane direction) results in strong Buneman and electron-electron streaming instabilities at the RD which suppress electron beam formation and therefore the striations in the exhaust. The strength of the streaming instabilities at the RD are controlled by the ratio of the electron thermal speed to Alfvén speed, lower thermal speed being more unstable. In the 3D simulations an ion-ion streaming instability acts to partially thermalize the counterstreaming ion beams at the SSs. This instability is controlled by the ratio of the sound speed to Alfvén speed and is expected to be stable in the low $β$ solar corona. The results suggest that in a guide field reconnection exhaust with $1 \gg β> m_e/m_i$, the kinetic-scale turbulence that develops will be too weak to play a significant role in energy conversion and particle acceleration. Therefore, the energy conversion will be mostly controlled by laminar physics or multi-x-line reconnection.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

We present two-dimensional kinetic simulations, with a broad range of initial guide fields, that isolate the role of parallel electric fields ($E_\parallel$) in energetic electron production during collisionless magnetic reconnection. In the strong guide field regime, $E_\parallel$ drives essentially all of the electron energy gain, yet fails to generate an energetic component. We suggest that this is due to the weak energy scaling of particle acceleration from $E_\parallel$ compared to that of a Fermi-type mechanism responsible for energetic electron production in the weak guide-field regime. This result has important implications for energetic electron production in astrophysical systems and reconnection-driven dissipation in turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

We investigate magnetic reconnection in systems simultaneously containing asymmetric (anti-parallel) magnetic fields, asymmetric plasma densities and temperatures, and arbitrary in-plane bulk flow of plasma in the upstream regions. Such configurations are common in the high-latitudes of Earth's magnetopause and in tokamaks. We investigate the convection speed of the X-line, the scaling of the reconnection rate, and the condition for which the flow suppresses reconnection as a function of upstream flow speeds. We use two-dimensional particle-in-cell simulations to capture the mixing of plasma in the outflow regions better than is possible in fluid modeling. We perform simulations with asymmetric magnetic fields, simulations with asymmetric densities, and simulations with magnetopause-like parameters where both are asymmetric. For flow speeds below the predicted cutoff velocity, we find good scaling agreement with the theory presented in Doss et al., J.~Geophys.~Res., 120, 7748 (2015). Applications to planetary magnetospheres, tokamaks, and the solar wind are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

Understanding the thermodynamic state of the hot intracluster medium (ICM) in a galaxy cluster requires a knowledge of the plasma transport processes, especially thermal conduction. The basic physics of thermal conduction in plasmas with ICM-like conditions has yet to be elucidated, however. We use particle-in-cell simulations and analytic models to explore the dynamics of an ICM-like plasma (with small gyroradius, large mean-free-path, and strongly sub-dominant magnetic pressure) driven by the diffusive heat flux associated with thermal conduction. Lin- ear theory reveals that whistler waves are driven unstable electron heat flux, even when the heat flux is weak. The resonant interaction of electrons with these waves then plays a critical role in scattering electrons and suppressing the heat flux. In a 1D model where only whistler modes that are parallel to the magnetic field are captured, the only resonant electrons are moving in the opposite direction to the heat flux and the electron heat flux suppression is small. In 2D or more, oblique whistler modes also resonate with electrons moving in the direction of the heat flux. The overlap of resonances leads to effective symmetrization of the electron distribution function and a strong suppression of heat flux. The results suggest that thermal conduction in the ICM might be strongly suppressed, possibly to negligible levels.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

Two- and three-dimensional particle-in-cell simulations of a recent encounter of the Magnetospheric Multiscale Mission (MMS) with an electron diffusion region at the magnetopause are presented. While the two-dimensional simulation is laminar, turbulence develops at both the x-line and along the magnetic separatrices in the three-dimensional simulation. The turbulence is strong enough to make the magnetic field around the reconnection island chaotic and produces both anomalous resistivity and anomalous viscosity. Each contribute significantly to breaking the frozen-in condition in the electron diffusion region. A surprise is that the crescent-shaped features in velocity space seen both in MMS observations and in two-dimensional simulations survive, even in the turbulent environment of the three-dimensional system. This suggests that MMS's measurements of crescent distributions do not exclude the possibility that turbulence plays an important role in magnetopause reconnection.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A model of the heliosphere with jets

An analytic model of the heliosheath (HS) between the termination shock (TS) and the heliopause (HP) is developed in the limit in which the interstellar flow and magnetic field are neglected. The heliosphere in this limit is axisymmetric and the overall structure of the HS and HP are controlled by the solar magnetic field even in the limit in which the ratio of the plasma to magnetic field pressure, $β=8πP/B^2$, in the HS is large. The tension of the solar magnetic field produces a drop in the total pressure between the TS and the HP. This same pressure drop accelerates the plasma flow downstream of the TS into the North and South directions to form two collimated jets. The radii of these jets are controlled by the flow through the TS and the acceleration of this flow by the magnetic field -- a stronger solar magnetic field boosts the velocity of the jets and reduces the radii of the jets and the HP. Magnetohydrodynamic (MHD) simulations of the global helioshere embedded in a stationary interstellar medium match well with the analytic model. The results suggest that mechanisms that reduce the HS plasma pressure downstream of the TS can enhance the jet outflow velocity and reduce the HP radius to values more consistent with the Voyager 1 observations than in current global models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

Kinetic simulations of 3D collisionless magnetic reconnection with a guide field show a dramatic enhancement of energetic electron production when compared with 2D systems. In the 2D systems, electrons are trapped in magnetic islands that limit their energy gain, whereas in the 3D systems the filamentation of the current layer leads to a stochastic magnetic field that enables the electrons to access volume-filling acceleration regions. The dominant accelerator of the most energetic electrons is a Fermi-like mechanism associated with reflection of charged particles from contracting field lines.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantifying Gyrotropy in Magnetic Reconnection

A new scalar measure of the gyrotropy of a pressure tensor is defined. Previously suggested measures are shown to be incomplete by means of examples for which they give unphysical results. To demonstrate its usefulness as an indicator of magnetic topology, the new measure is calculated for electron data taken from numerical simulations of magnetic reconnection, shown to peak at separatrices and X-points, and compared to the other measures. The new diagnostic has potential uses in analyzing spacecraft observations and so a method for calculating it from measurements performed in an arbitrary coordinate system is derived.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coronal Electron Confinement by Double Layers

In observations of flare-heated electrons in the solar corona, a longstanding problem is the unexplained prolonged lifetime of the electrons compared to their transit time across the source. This suggests confinement. Recent particle-in-cell (PIC) simulations, which explored the transport of pre-accelerated hot electrons through ambient cold plasma, showed that the formation of a highly localized electrostatic potential drop, in the form of a double layer (DL), significantly inhibited the transport of hot electrons (T.C. Li, J.F. Drake, and M. Swisdak, 2012, ApJ, 757, 20). The effectiveness of confinement by a DL is linked to the strength of the DL as defined by its potential drop. In this work, we investigate the scaling of the DL strength with the hot electron temperature by PIC simulations, and find a linear scaling. We demonstrate that the strength is limited by the formation of parallel shocks. Based on this, we analytically determine the maximum DL strength, and find also a linear scaling with the hot electron temperature. The DL strength obtained from the analytic calculation is comparable to that from the simulations. At the maximum strength, the DL is capable of confining a significant fraction of hot electrons in the source.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

Observations of flare-heated electrons in the corona typically suggest confinement of electrons. The confinement mechanism, however, remains unclear. The transport of coronal hot electrons into ambient plasma was recently investigated by particle-in-cell (PIC) simulations. Electron transport was significantly suppressed by the formation of a highly localized, nonlinear electrostatic electric potential in the form of a double layer (DL). In this work large-scale PIC simulations are performed to explore the dynamics of DLs in larger systems where, instead of a single DL, multiple DLs are generated. The primary DL accelerates return current electrons, resulting in high velocity electron beams that interact with ambient ions. This forms a Buneman unstable system that spawns more DLs. Trapping of heated return current electrons between multiple DLs strongly suppresses electron transport. DLs also accelerate ambient ions and produce strong ion flows over an extended region. This clarifies the mechanism by which hot electrons in the corona couple to and accelerate ions to form the solar wind. These new dynamics in larger systems reveal a more likely picture of DL development and their impact on the ambient plasma in the solar corona. They are applicable to the preparation for in-situ coronal space missions like the Solar Probe Plus.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the 3-D structure and dissipation of reconnection-driven flow-bursts

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Mechanisms of Electron Heating and Acceleration during Magnetic Reconnection

The heating of electrons in collisionless magnetic reconnection is explored in particle-in-cell (PIC) simulations with non-zero guide fields so that electrons remain magnetized. In this regime electric fields parallel to B accelerate particles directly while those perpendicular to B do so through gradient-B and curvature drifts. The curvature drift drives parallel heating through Fermi reflection while the gradient B drift changes the perpendicular energy through betatron acceleration. We present simulations in which we evaluate each of these mechanisms in space and time in order to quantify their role in electron heating. For a case with a small guide field (20 % of the magnitude of the reconnecting component) the curvature drift is the dominant source of electron heating. However, for a larger guide field (equal to the magnitude of the reconnecting component) electron acceleration by the curvature drift is comparable to that of the parallel electric field. In both cases the heating by the gradient B drift is negligible in magnitude. Heating by the curvature-drift dominates in the outflow exhausts where bent field lines expand to relax their tension and is therefore distributed over a large area. In contrast, the parallel electric field is localized near X-lines. Acceleration by parallel electric fields may play a smaller role in large systems where the X-line occupies a vanishing fraction of the system. The curvature drift and the parallel electric field dominate the dynamics and drive parallel heating. A consequence is that the electron energy spectrum becomes extremely anisotropic at late time, which has important implications for quantifying the limits of electron acceleration due to synchrotron emission. An upper limit on electron energy gain that is substantially higher than earlier estimates is obtained by balancing reconnection drive with radiative loss.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

The onset of the acceleration of ions during magnetic reconnection is explored via particle-in-cell simulations in the limit of a strong ambient guide field that self-consistently and simultaneously follow the motions of protons and $α$ particles. Heating parallel to the local magnetic field during reconnection with a guide field is strongly reduced compared with the reconnection of anti-parallel magnetic fields. The dominant heating of thermal ions during guide field reconnection results from pickup behavior of ions during their entry into reconnection exhausts and dominantly produces heating perpendicular rather than parallel to the local magnetic field. Pickup behavior requires that the ion transit time across the exhaust boundary (with a transverse scale of the order of the ion sound Larmor radius) be short compared with the ion cyclotron period. This translates into a threshold in the strength of reconnecting magnetic field that favors the heating of ions with high mass-to-charge. A simulation with a broad initial current layer produces a reconnecting system in which the amplitude of the reconnecting magnetic field just upstream of the dissipation region increases with time as reconnection proceeds. The sharp onset of perpendicular heating when the pickup threshold is crossed is documented. A comparison of the time variation of the parallel and perpendicular ion heating with that predicted based on the strength of the reconnecting field establishes the scaling of ion heating with ambient parameters both below and above the pickup threshold. The relevance to observations of ion heating in the solar corona is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Porous, Layered Heliopause

The picture of the heliopause (HP) -- the boundary between the domains of the sun and the local interstellar medium (LISM) -- as a pristine interface with a large rotation in the magnetic field fails to describe recent Voyager 1 (V1) spacecraft data. Magnetohydrodynamic (MHD) simulations of the global heliosphere reveal that the rotation angle of the magnetic field across the HP at V1 is small. Particle-in-cell simulations, based on cuts through the MHD model at the location of V1, suggest that the sectored region of the heliosheath (HS) produces large-scale magnetic islands that reconnect with the interstellar magnetic field and mix LISM and HS plasma. Cuts across the simulation data reveal multiple, anti-correlated jumps in the number densities of LISM and HS particles at the magnetic separatrices of the islands, similar to those observed by V1. A model is presented, based on both the observations and simulation data, of the HP as a porous, multi-layered structure threaded by magnetic fields. This model further suggests that, contrary to the conclusions of recent papers, V1 has already crossed the HP.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Adiabatic Phase Mixing and Heating of Electrons in Buneman Turbulence

The nonlinear development of the strong Buneman instability and the associated fast electron heating in thin current layers with $Ω_e/ω_{pe} <1$ are explored. Phase mixing of the electrons in wave potential troughs and a rapid increase in temperature are observed during the saturation of the instability. We show that the motion of trapped electrons can be described using a Hamiltonian formalism in the adiabatic approximation. The process of separatrix crossing as electrons are trapped and de-trapped is irreversible and guarantees that the resulting electron energy gain is a true heating process.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Scaling of the growth rate of magnetic islands in the heliosheath

Current sheets thinner than the ion inertial length are unstable to the tearing instability and will develop magnetic islands that grow due to magnetic reconnection. We investigate whether the growth of magnetic islands in a current sheet can continue indefinitely, or in the case of the heliosheath until reaching a neighboring current sheet, and at what rate the islands grow. We investigate the development and growth of magnetic islands using a particle-in-cell code, starting from particle noise. Performing a scaling of the growth of magnetic islands versus the system size, we find that the growth rate is independent of the system size up to the largest simulation we were able to complete. The islands are able to continue growing as long as they merge with each other and maintain a high aspect ratio. Otherwise there is not enough magnetic tension to sustain reconnection. When applied to the sectored magnetic fields in the heliosheath, we show that the islands can continue growing until they reach the sector width and do so in much less time than it takes for the islands to convect through the heliosheath.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Suppression of energetic electron transport in flares by double layers

During flares and coronal mass ejections, energetic electrons from coronal sources typically have very long lifetimes compared to the transit times across the systems, suggesting confinement in the source region. Particle-in-cell simulations are carried out to explore the mechanisms of energetic electron transport from the corona to the chromosphere and possible confinement. We set up an initial system of pre-accelerated hot electrons in contact with ambient cold electrons along the local magnetic field, and let it evolve over time. Suppression of transport by a nonlinear, highly localized electrostatic electric field (in the form of a double layer) is observed after a short phase of free-streaming by hot electrons. The double layer (DL) emerges at the contact of the two electron populations. It is driven by an ion-electron streaming instability due to the drift of the back-streaming return current electrons interacting with the ions. The DL grows over time and supports a significant drop in temperature and hence reduces heat flux between the two regions that is sustained for the duration of the simulation. This study shows transport suppression begins when the energetic electrons start to propagate away from a coronal acceleration site. It also implies confinement of energetic electrons with kinetic energies less than the electrostatic energy of the DL for the DL lifetime, which is much longer than the electron transit time through the source region.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The effects of plasma beta and anisotropy instabilities on the dynamics of reconnecting magnetic fields in the heliosheath

The plasma β (the ratio of the plasma pressure to the magnetic pressure) of a system can have a large effect on its dynamics as high β enhances the effects of pressure anisotropies. We investigate the effects of β in a system of stacked current sheets that break up into magnetic islands due to magnetic reconnection. We find significant differences between β < 1 and β > 1. At low β growing magnetic islands are modestly elongated and become round as contraction releases magnetic stress and reduces magnetic energy. At high β the increase of the parallel pressure in contracting islands causes saturation of modestly elongated islands as island cores approach the marginal firehose condition. Only highly elongated islands reach finite size. The kinking associated with the Weibel and firehose instabilities prevents full contraction of these islands, leading to a final state of highly elongated islands in which further reconnection is suppressed. The results are directly relevant to reconnection in the sectored region of the heliosheath and possibly to saturation mechanisms of the magnetorotational instability in accretion flows.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The power-law spectra of energetic particles during multi-island magnetic reconnection

Power-law distributions are a near universal feature of energetic particle spectra in the heliosphere. Anomalous Cosmic Rays (ACRs), super-Alfvénic ions in the solar wind and the hardest energetic electron spectra in flares all have energy fluxes with power-laws that depend on energy $E$ approximately as $E^{-1.5}$. We present a new model of particle acceleration in systems with a bath of merging magnetic islands that self-consistently describes the development of velocity-space anisotropy parallel and perpendicular to the local magnetic field and includes the self-consistent feedback of pressure anisotropy on the merging dynamics. By including pitch-angle scattering we obtain an equation for the omni-directional particle distribution $f(v,t)$ that is solved in closed form to reveal $v^{-5}$ (corresponding to an energy flux varying as $E^{-1.5}$) as a near-universal solution as long as the characteristic acceleration time is short compared with the characteristic loss time. In such a state the total energy in the energetic particles reaches parity with the remaining magnetic free energy. More generally, the resulting transport equation can serve as the basis for calculating the distribution of energetic particles resulting from reconnection in large-scale inhomogeneous systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Is the magnetic field in the heliosheath laminar or a turbulent bath of bubbles?

All the current global models of the heliosphere are based on the assumption that the magnetic field in the heliosheath, in the region close to the heliopause is laminar. We argue that in that region the heliospheric magnetic field is not laminar but instead consists of magnetic bubbles. Recently, we proposed that the annihilation of the "sectored" magnetic field within the heliosheath as it is compressed on its approach to the heliopause produces the anomalous cosmic rays and also energetic electrons. As a product of the annihilation of the sectored magnetic field, densely-packed magnetic islands/bubbles are produced. These magnetic islands/bubbles will be convected with the ambient flows as the sector region is carried to higher latitudes filling the heliosheath. We further argue that the magnetic islands/bubbles will develop upstream within the heliosheath. As a result, the magnetic field in the heliosheath sector region will be disordered well upstream of the heliopause. We present a 3D MHD simulation with very high numerical resolution that captures the north-south boundaries of the sector region. We show that due to the high pressure of the interstellar magnetic field a north-south asymmetry develops such that the disordered sectored region fills a large portion of the northern part of the heliosphere with a smaller extension in the southern hemisphere. We suggest that this scenario is supported by the following changes that occur around 2008 and from 2009.16 onward: a) the sudden decrease in the intensity of low energy electrons detected by Voyager 2; b) a sharp reduction in the intensity of fluctuations of the radial flow; and c) the dramatic differences in intensity trends between GCRs at V1 and 2. We argue that these observations are a consequence of V2 leaving the sector region of disordered field during these periods and crossing into a region of unipolar laminar field.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The effects of strong temperature anisotropy on the kinetic structure of collisionless slow shocks and reconnection exhausts. Part I: PIC simulations

A 2-D Riemann problem is designed to study the development and dynamics of the slow shocks that are thought to form at the boundaries of reconnection exhausts. Simulations are carried out for varying ratios of normal magnetic field to the transverse upstream magnetic field (i.e., propagation angle with respect to the upstream magnetic field). When the angle is sufficiently oblique, the simulations reveal a large firehose-sense (P_parallel>P_perpendicular) temperature anisotropy in the downstream region, accompanied by a transition from a coplanar slow shock to a non-coplanar rotational mode. In the downstream region the firehose stability parameter epsilon=1-mu_0(P_parallel-P_perpendicular)/ B^2 tends to lock in to 0.25. This balance arises from the competition between counterstreaming ions, which drives epsilon down, and the scattering due to ion inertial scale waves, which are driven unstable by the downstream rotational wave. At very oblique propagating angles, 2-D turbulence also develops in the downstream region.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The effects of strong temperature anisotropy on the kinetic structure of collisionless slow shocks and reconnection exhausts. Part II: Theory

Simulations of collisionless oblique propagating slow shocks have revealed the existence of a transition associated with a critical temperature anisotropy epsilon=1-mu_0(P_parallel-P_perpendicular)/ B^2 = 0.25 (Liu, Drake and Swisdak (2011)). An explanation for this phenomenon is proposed here based on anisotropic fluid theory, in particular the Anisotropic Derivative Nonlinear-Schrodinger-Burgers equation, with an intuitive model of the energy closure for the downstream counter-streaming ions. The anisotropy value of 0.25 is significant because it is closely related to the degeneracy point of the slow and intermediate modes, and corresponds to the lower bound of the coplanar to non-coplanar transition that occurs inside a compound slow shock (SS)/rotational discontinuity (RD) wave. This work implies that it is a pair of compound SS/RD waves that bound the outflows in magnetic reconnection, instead of a pair of switch-off slow shocks as in Petschek's model. This fact might explain the rareness of in-situ observations of Petschek-reconnection-associated switch-off slow shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The structure of the magnetic reconnection exhaust boundary

The structure of shocks that form at the exhaust boundaries during collisionless reconnection of anti-parallel fields is studied using particle-in-cell (PIC) simulations and modeling based on the anisotropic magnetohydrodynamic equations. Large-scale PIC simulations of reconnection and companion Riemann simulations of shock development demonstrate that the pressure anisotropy produced by counterstreaming ions within the exhaust prevents the development of classical Petschek switch-off-slow shocks (SSS). The shock structure that does develop is controlled by the firehose stability parameter epsilon=1-mu_0(P_parallel-P_perpendicular)/ B^2 through its influence on the speed order of the intermediate and slow waves. Here P_parallel and P_perpendicular are the pressure parallel and perpendicular to the local magnetic field. The exhaust boundary is made up of a series of two shocks and a rotational wave. The first shock takes epsilon from unity upstream to a plateau of 0.25 downstream. The condition epsilon =0.25 is special because at this value the speeds of nonlinear slow and intermediate waves are degenerate. The second slow shock leaves epsilon=0.25 unchanged but further reduces the amplitude of the reconnecting magnetic field. Finally, in the core of the exhaust epsilon drops further and the transition is completed by a rotation of the reconnecting field into the out-of-plane direction. The acceleration of the exhaust takes place across the two slow shocks but not during the final rotation. The result is that the outflow speed falls below that expected from the Walen condition based on the asymptotic magnetic field. A simple analytic expression is given for the critical value of epsilon within the exhaust below which SSSs no longer bound the reconnection outflow.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Electron Holes and Heating in the Reconnection Dissipation Region

Using particle-in-cell simulations and kinetic theory, we explore the current-driven turbulence and associated electron heating in the dissipation region during 3D magnetic reconnection with a guide field. At late time the turbulence is dominated by the Buneman and lower hybrid instabilities. Both produce electron holes that co-exist but have very different propagation speeds. The associated scattering of electrons by the holes enhances electron heating in the dissipation region.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Vector Direction of the Interstellar Magnetic Field Outside the Heliosphere

We propose that magnetic reconnection at the heliopause only occurs where the interstellar magnetic field points nearly anti-parallel to the heliospheric field. By using large-scale magnetohydrodynamic (MHD) simulations of the heliosphere to provide the initial conditions for kinetic simulations of heliopause (HP) reconnection we show that the energetic pickup ions downstream from the solar wind termination shock induce large diamagnetic drifts in the reconnecting plasma and stabilize non-anti-parallel reconnection. With this constraint the MHD simulations can show where HP reconnection most likely occurs. We also suggest that reconnection triggers the 2-3 kHz radio bursts that emanate from near the HP. Requiring the burst locations to coincide with the loci of anti-parallel reconnection allows us to determine, for the first time, the vector direction of the local interstellar magnetic field. We find it to be oriented towards the southern solar magnetic pole.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Three-dimensional simulations of the orientation and structure of reconnection X-lines

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Statistical Model of Magnetic Islands in a Large Current Layer

We develop a statistical model describing the dynamics of magnetic islands in very large current layers that develop in space plasma. Two parameters characterize the island distribution: the flux contained in the island and the area it encloses. We derive an integro-differential evolution equation for this distribution function, based on rules that govern the small-scale generation of secondary islands, the rates of island growth, and island merging. Our numerical solutions of this equation produce island distributions relevant to the magnetosphere and corona. We also derive and analytically solve a differential equation for large islands that explicitly shows the role merging plays in island growth.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dissipation of the sectored heliospheric magnetic field near the heliopause: a mechanism for the generation of anomalous cosmic rays

The recent observations of the anomalous cosmic ray (ACR) energy spectrum as Voyagers 1 and 2 crossed the heliospheric termination shock have called into question the conventional shock source of these energetic particles. We suggest that the sectored heliospheric magnetic field, which results from the flapping of the heliospheric current sheet, piles up as it approaches the heliopause, narrowing the current sheets that separate the sectors and triggering the onset of collisionless magnetic reconnection. Particle-in-cell simulations reveal that most of the magnetic energy is released and most of this energy goes into energetic ions with significant but smaller amounts of energy going into electrons. The energy gain of the most energetic ions results from their reflection from the ends of contracting magnetic islands, a first order Fermi process. The energy gain of the ions in contracting islands increases their parallel (to the magnetic field ${\bf B}$) pressure $p_\parallel$ until the marginal firehose condition is reached, causing magnetic reconnection and associated particle acceleration to shut down. The model calls into question the strong scattering assumption used to derive the Parker transport equation and therefore the absence of first order Fermi acceleration in incompressible flows. A simple 1-D model for particle energy gain and loss is presented in which the feedback of the energetic particles on the reconnection drive is included. The ACR differential energy spectrum takes the form of a power law with a spectral index slightly above 1.5. The model has the potential to explain several key Voyager observations, including the similarities in the spectra of different ion species.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nonlinear Development of Streaming Instabilities In Strongly Magnetized Plasmas

The nonlinear development of streaming instabilities in the current layers formed during magnetic reconnection with a guide field is explored. Theory and 3-D particle-in-cell simulations reveal two distinct phases. First, the parallel Buneman instability grows and traps low velocity electrons. The remaining electrons then drive two forms of turbulence: the parallel electron-electron two-stream instability and the nearly-perpendicular lower hybrid instability. The high velocity electrons resonate with the turbulence and transfer momentum to the ions and low velocity electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The Orientation of the Reconnection X-line

We propose a criterion for identifying the orientation of the X-line when two regions of plasma with arbitrary densities, temperatures, and magnetic fields undergo reconnection. The X-line points in the direction that maximizes the (suitably-defined) Alfvén speed characterizing the reconnection outflow. For many situations a good approximation is that the X-line bisects the angle formed by the magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Simulation Study of a Positive Ionospheric Storm Phase Observed at Millstone Hill

Simulation results from the NRL ionospheric model SAMI2 indicate that the changes in the F-region over Millstone Hill during the geomagnetic storm beginning on 3 April 2004 are primarily due to the influence of a long-lasting eastward electric field, as was previously suggested by {\it Huang et al.} [2005]. A simulation of the storm day agrees well with the observational data and shows that, compared with the ionosphere of the previous quiet day, the peak electron density in the F-region (NmF2) increased by a factor of $\approx 2$, the altitude of the peak density (hmF2) rose by $\approx 80$ km, and the F-region electron temperature decreased by $\approx 1000$ K. Further simulations in which either the neutral atmosphere and winds or the electric field were replaced by their quiet day counterparts clearly suggest that the electric field played the dominant, although not exclusive, role in producing these effects.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access