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

J. F. Drake

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

ResearcherUnclaimed source record
physics.space-phphysics.plasm-phastro-ph.SRastro-ph.EPastro-ph.HEastro-phastro-ph.GAnlin.PSastro-ph.COnlin.AO

Catalog footprint

What is connected

49works
10topics
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

49 published item(s)

preprint2021arXiv

A solar source of Alfvénic magnetic field switchbacks: {\em in situ} remnants of magnetic funnels on supergranulation scales

One of the striking observations from the Parker Solar Probe (PSP) spacecraft is the prevalence in the inner heliosphere of large amplitude, Alfvénic magnetic field reversals termed 'switchbacks'. These $δB_R/B \sim \mathcal{O}(1$) fluctuations occur on a range of timescales and in {\em patches} separated by intervals of quiet, radial magnetic field. We use measurements from PSP to demonstrate that patches of switchbacks are localized within the extensions of plasma structures originating at the base of the corona. These structures are characterized by an increase in alpha particle abundance, Mach number, plasma $β$ and pressure, and by depletions in the magnetic field magnitude and electron temperature. These intervals are in pressure-balance, implying stationary spatial structure, and the field depressions are consistent with overexpanded flux tubes. The structures are asymmetric in Carrington longitude with a steeper leading edge and a small ($\sim$1$^\circ$) edge of hotter plasma and enhanced magnetic field fluctuations. Some structures contain suprathermal ions to $\sim$85 keV that we argue are the energetic tail of the solar wind alpha population. The structures are separated in longitude by angular scales associated with supergranulation. This suggests that these switchbacks originate near the leading edge of the diverging magnetic field funnels associated with the network magnetic field - the primary wind sources. We propose an origin of the magnetic field switchbacks, hot plasma and suprathermals, alpha particles in interchange reconnection events just above the solar transition region and our measurements represent the extended regions of a turbulent outflow exhaust.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Suprathermal Ion Energy spectra and Anisotropies near the Heliospheric Current Sheet crossing observed by the Parker Solar Probe during Encounter 7

We present observations of >10-100 keV/nucleon suprathermal (ST) H, He, O, and Fe ions associated with crossings of the heliospheric current sheet (HCS) at radial distances <0.1 au from the Sun. Our key findings are: 1) very few heavy ions are detected during the 1st full crossing, the heavy ion intensities are reduced during the 2nd partial crossing and peak just after the 2nd crossing; 2) ion arrival times exhibit no velocity dispersion; 3) He pitch-angle distributions track the magnetic field polarity reversal and show up to ~10:1 anti-sunward, field-aligned flows and beams closer to the HCS that become nearly isotropic further from the HCS; 4) the He spectrum steepens either side of the HCS and the He, O, and Fe spectra exhibit power-laws of the form ~E^4-6; and 5) maximum energies EX increase with the ion's charge-to-mass (Q/M) ratio as EX/EH proportional to [(QX/MX)]^alpha where alpha~0.65-0.76, assuming that the average Q-states are similar to those measured in gradual and impulsive solar energetic particle events at 1 au. The absence of velocity dispersion in combination with strong field-aligned anisotropies closer to the HCS appears to rule out solar flares and near-sun coronal mass ejection-driven shocks. These new observations present challenges not only for mechanisms that employ direct parallel electric fields and organize maximum energies according to E/Q, but also for local diffusive and magnetic reconnection-driven acceleration models. Re-evaluation of our current understanding of the production and transport of energetic ions is necessary to understand this near-solar, current-sheet-associated population of ST ions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The development of a split-tail heliosphere and the role of non-ideal processes: a comparison of the BU and Moscow models

Global models of the heliosphere are critical tools used in the interpretation of heliospheric observations. There are several three-dimensional magnetohydrodynamic (MHD) heliospheric models that rely on different strategies and assumptions. Until now only one paper has compared global heliosphere models, but without magnetic field effects. We compare the results of two different MHD models, the BU and Moscow models. Both models use identical boundary conditions to compare how different numerical approaches and physical assumptions contribute to the heliospheric solution. Based on the different numerical treatments of discontinuities, the BU model allows for the presence of magnetic reconnection, while the Moscow model does not. Both models predict collimation of the solar outflow in the heliosheath by the solar magnetic field and produce a split-tail where the solar magnetic field confines the charged solar particles into distinct north and south columns that become lobes. In the BU model, the ISM flows between the two lobes at large distances due to MHD instabilities and reconnection. Reconnection in the BU model at the port flank affects the draping of the interstellar magnetic field in the immediate vicinity of the heliopause. Different draping in the models cause different ISM pressures, yielding different heliosheath thicknesses and boundary locations, with the largest effects at high latitudes. The BU model heliosheath is 15% thinner and the heliopause is 7% more inwards at the north pole relative to the Moscow model. These differences in the two plasma solutions may manifest themselves in energetic neutral atom measurements of the heliosphere.

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sunward propagating whistler waves collocated with localized magnetic field holes in the solar wind: Parker Solar Probe observations at 35.7 Sun radii

Observations by the Parker Solar Probe mission of the solar wind at about 35.7 solar radii reveal the existence of whistler wave packets with frequencies below 0.1 f/fce (20-80 Hz in the spacecraft frame). These waves often coincide with local minima of the magnetic field magnitude or with sudden deflections of the magnetic field that are called switchbacks. Their sunward propagation leads to a significant Doppler frequency downshift from 200-300 Hz to 20-80 Hz (from 0.2 f/fce to 0.5 f/fce). The polarization of these waves varies from quasi-parallel to significantly oblique with wave normal angles that are close to the resonance cone. Their peak amplitude can be as large as 2 to 4 nT. Such values represent approximately 10% of the background magnetic field, which is considerably more than what is observed at 1 a.u. Recent numerical studies show that such waves may potentially play a key role in breaking the heat flux and scattering the Strahl population of suprathermal electrons into a halo population.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Confinement of the Heliosheath Plasma by the Solar Magnetic Field as Revealed by Energetic Neutral Atom Simulations

Traditionally, the solar magnetic field has been considered to have a negligible effect in the outer regions of the heliosphere. Recent works have shown that the solar magnetic field may play a crucial role in collimating the plasma in the heliosheath. Interstellar Boundary Explorer (IBEX) observations of the heliotail indicated a latitudinal structure varying with energy in the energetic neutral atom (ENA) fluxes. At energies ~1 keV, the ENA fluxes show an enhancement at low latitudes and a deficit of ENAs near the poles. At energies >2.7 keV, ENA fluxes had a deficit within low latitudes, and lobes of higher ENA flux near the poles. This ENA structure was initially interpreted to be a result of the latitudinal profile of the solar wind during solar minimum. We extend the work of Kornbleuth et al. (2018) by using solar minimum-like conditions and the recently developed SHIELD model. The SHIELD model couples the magnetohydrodynamic (MHD) plasma solution with a kinetic description of neutral hydrogen. We show that while the latitudinal profile of the solar wind during solar minimum contributes to the lobes in ENA maps, the collimation by the solar magnetic field is important in creating and shaping the two high latitude lobes of enhanced ENA flux observed by IBEX. This is the first work to explore the effect of the changing solar magnetic field strength on ENA maps. Our findings suggest that IBEX is providing the first observational evidence of the collimation of the heliosheath plasma by the solar magnetic field.

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
preprint2019arXiv

Turbulence and Transport During Guide-Field Reconnection at the Magnetopause

We analyze the development and influence of turbulence in three-dimensional particle-in-cell simulations of guide-field magnetic reconnection at the magnetopause with parameters based on observations of an electron diffusion region by the Magnetospheric Multiscale (MMS) mission. Along the separatrices the turbulence is a variant of the lower hybrid drift instability (LHDI) that produces electric field fluctuations with amplitudes much greater than the reconnection electric field. The turbulence controls the scale length of the density and current profiles while enabling significant transport across the magnetopause despite the electrons remaining frozen-in to the magnetic field. Near the X-line the electrons are not frozen-in and the turbulence, which differs from the LHDI, makes a significant net contribution to the generalized Ohm's law through an anomalous viscosity. The characteristics of the turbulence and associated particle transport are consistent with fluctuation amplitudes in the MMS observations. However, for this event the simulations suggest that the MMS spacecraft were not close enough to the core of the electron diffusion region to identify the region where anomalous viscosity is important.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Universality of lower hybrid waves at Earth's magnetopause

Waves around the lower hybrid frequency are frequently observed at Earth's magnetopause, and readily reach very large amplitudes. Determining the properties of lower hybrid waves is crucial because they are thought to contribute to electron and ion heating, cross-field particle diffusion, anomalous resistivity, and energy transfer between electrons and ions. All these processes could play an important role in magnetic reconnection at the magnetopause and the evolution of the boundary layer. In this paper, the properties of lower hybrid waves at Earth's magnetopause are investigated using the Magnetospheric Multiscale (MMS) mission. For the first time, the properties of the waves are investigated using fields and direct particle measurements. The highest-resolution electron moments resolve the velocity and density fluctuations of lower hybrid waves, confirming that electrons remain approximately frozen in at lower hybrid wave frequencies. Using fields and particle moments the dispersion relation is constructed and the wave-normal angle is estimated to be close to $90^{\circ}$ to the background magnetic field. The waves are shown to have a finite parallel wave vector, suggesting that they can interact with parallel propagating electrons. The observed wave properties are shown to agree with theoretical predictions, the previously used single-spacecraft method, and four-spacecraft timing analyses. These results show that single-spacecraft methods can accurately determine lower hybrid wave properties.

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

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

Asymmetric magnetic reconnection with a flow shear and applications to the magnetopause

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.

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

Magnetized jets driven by the sun: the structure of the heliosphere revisited

The classic accepted view of the heliosphere is a quiescent, comet-like shape aligned in the direction of the Sun's travel through the interstellar medium (ISM) extending for 1000's of AUs (AU: astronomical unit). Here we show, based on magnetohydrodynamic (MHD) simulations, that the tension (hoop) force of the twisted magnetic field of the sun confines the solar wind plasma beyond the termination shock and drives jets to the North and South very much like astrophysical jets. These jets are deflected into the tail region by the motion of the Sun through the ISM similar to bent galactic jets moving through the intergalactic medium. The interstellar wind blows the two jets into the tail but is not strong enough to force the lobes into a single comet-like tail, as happens to some astrophysical jets (Morsony et al. 2013). Instead, the interstellar wind flows around the heliosphere and into equatorial region between the two jets. As in some astrophysical jets that are kink unstable (Porth et al. 2014) we show here that the heliospheric jets are turbulent (due to large-scale MHD instabilities and reconnection) and strongly mix the solar wind with the ISM beyond 400 AU. The resulting turbulence has important implications for particle acceleration in the heliosphere. The two-lobe structure is consistent with the energetic neutral atoms (ENAs) images of the heliotail from IBEX (McComas et al. 2013) where two lobes are visible in the North and South and the suggestion from the CASSINI (Krimigis et al. 2009, Dialynas et al. 2013) ENAs that the heliosphere is "tailless".

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 Interface Region Imaging Spectrograph (IRIS)

The Interface Region Imaging Spectrograph (IRIS) small explorer spacecraft provides simultaneous spectra and images of the photosphere, chromosphere, transition region, and corona with 0.33-0.4 arcsec spatial resolution, 2 s temporal resolution and 1 km/s velocity resolution over a field-of-view of up to 175 arcsec x 175 arcsec. IRIS was launched into a Sun-synchronous orbit on 27 June 2013 using a Pegasus-XL rocket and consists of a 19-cm UV telescope that feeds a slit-based dual-bandpass imaging spectrograph. IRIS obtains spectra in passbands from 1332-1358, 1389-1407 and 2783-2834 Angstrom including bright spectral lines formed in the chromosphere (Mg II h 2803 Angstrom and Mg II k 2796 Angstrom) and transition region (C II 1334/1335 Angstrom and Si IV 1394/1403 Angstrom). Slit-jaw images in four different passbands (C II 1330, Si IV 1400, Mg II k 2796 and Mg II wing 2830 Angstrom) can be taken simultaneously with spectral rasters that sample regions up to 130 arcsec x 175 arcsec at a variety of spatial samplings (from 0.33 arcsec and up). IRIS is sensitive to emission from plasma at temperatures between 5000 K and 10 MK and will advance our understanding of the flow of mass and energy through an interface region, formed by the chromosphere and transition region, between the photosphere and corona. This highly structured and dynamic region not only acts as the conduit of all mass and energy feeding into the corona and solar wind, it also requires an order of magnitude more energy to heat than the corona and solar wind combined. The IRIS investigation includes a strong numerical modeling component based on advanced radiative-MHD codes to facilitate interpretation of observations of this complex region. Approximately eight Gbytes of data (after compression) are acquired by IRIS each day and made available for unrestricted use within a few days of the observation.

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

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.

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

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.

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

Energetic Protons, Radionuclides and Magnetic Activity in Protostellar Disks

We calculate the location of the magnetically-inactive dead zone in the minimum-mass protosolar disk, under ionization scenarios including stellar X-rays, long- or short-lived radionuclide decay, and energetic protons arriving from the general interstellar medium, from a nearby supernova explosion, from the disk corona, or from the corona of the young star. The disk contains a dead zone in all scenarios except those with small dust grains removed and a fraction of the short-lived radionuclides remaining in the gas. All the cases without exception have an "undead zone" where intermediate resistivities prevent magneto-rotational turbulence while allowing shear-generated large-scale magnetic fields. The mass column in the undead zone is typically greater than the column in the turbulent surface layers. The results support the idea that the dead and undead zones are robust consequences of cold, dusty gas with mass columns exceeding 1000 g/cm^2.

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

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