Search

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., Comment: 13 pages, 4 figures

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., Comment: 5 pages, 2 figures, accepted for publication as a letter in Physics of Plasmas

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"., Comment: 19 pages, 5 figures; Astrophysical Journal Letters; in press

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\'en speed $c_{Ar}$ based on the reconnecting magnetic field through the relation $\Delta T_e = 0.033 \,m_i\,c_{Ar}^2$, where $\Delta 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 $\beta$. 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 $\Delta T_e$ was proportional to the inflowing Alfv\'en 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., Comment: 27 pages, 1 table, 7 figures. Submitted to Physics of Plasmas

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.

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., Comment: Submitted to Physics of Plasmas (19 pages, 13 figures)

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 $\alpha$ 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., Comment: arXiv admin note: text overlap with arXiv:1108.5750

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.

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.

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., Comment: 53 pages, 15 figures

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., Comment: 16 pages, 5 figures, accepted to Astrophysical Journal Letters in August, 2013. This version is the accepted version

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., Comment: 7 pages, 4 figures

The nonlinear development of the strong Buneman instability and the associated fast electron heating in thin current layers with $\Omega_e/\omega_{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., Comment: 4 figs, 11 pages, Accepted by Physics of Plasma

Power-law distributions are a near universal feature of energetic particle spectra in the heliosphere. Anomalous Cosmic Rays (ACRs), super-Alfv\'enic 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.

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.

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.

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., Comment: 13 pages, 5 figures

The plasma {\beta} (the ratio of the plasma pressure to the magnetic pressure) of a system can have a large effect on its dynamics as high {\beta} enhances the effects of pressure anisotropies. We investigate the effects of {\beta} in a system of stacked current sheets that break up into magnetic islands due to magnetic reconnection. We find significant differences between {\beta} < 1 and {\beta} > 1. At low {\beta} growing magnetic islands are modestly elongated and become round as contraction releases magnetic stress and reduces magnetic energy. At high {\beta} 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.

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., Comment: 18 pages, 10 figures

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\'en wave. For magnetotail parameters, the parallel propagation of this wave is super-Alfv\'enic (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., Comment: Submitted to PRL on 11/1/2010. Resubmitted on 4/5/2011

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., Comment: 36 pages, 15 figures, submitted to ApJ

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., Comment: 13 pages, 10 figures

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., Comment: Accepted by Physics of Plasmas

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., Comment: 14 pages, 5 figures, submitted to GRL

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., Comment: Submitted to ApJ; incorporates minor referee-suggested revisions

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., Comment: Submitted to ApJ; shortened abstract; degraded figure quality

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., Comment: 4 pages, 3 figures

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., Comment: Accepted for the Astrophysical Journal

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., Comment: Accepted by PRL

The mass symmetry between the two species in electron-positron (pair) plasmas has interesting consequences for collisionless magnetic reconnection because the Hall term, which plays a crucial role in supporting fast reconnection in electron-proton plasmas, vanishes. We perform kinetic simulations of pair reconnection in systems of various sizes, show that it remains fast, and identify the reason why this occurs. For sufficiently large systems a Weibel-like temperature anisotropy instability develops in the outflow from the X-point that causes the current layer to broaden and form a Petschek-like open outflow. We discuss the parameter regimes in which pair reconnection should be fast and the implications for astrophysical pair plasmas., Comment: Accepted by ApJ

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\'en speed characterizing the reconnection outflow. For many situations a good approximation is that the X-line bisects the angle formed by the magnetic fields.

Particle in cell (PIC) simulations of collisionless magnetic reconnection are presented that demonstrate that the electron dissipation region develops a distinct two-scale structure along the outflow direction. The length of the electron current layer is found to decrease with decreasing electron mass, approaching the ion inertial length for a proton-electron plasma. A surprise, however, is that the electrons form a high-velocity outflow jet that remains decoupled from the magnetic field and extends large distances downstream from the x-line. The rate of reconnection remains fast in very large systems, independent of boundary conditions and the mass of electrons., Comment: Submitted to Physical Review Letters, 4 pages, 4 figures

We present a model for the spontaneous onset of fast magnetic reconnection in a weakly collisional plasma, such as the solar corona. When a current layer of macroscopic width undergoes collisional (Sweet-Parker) reconnection, a narrow dissipation region forms around the X-line. This dissipation region naturally becomes narrower during the reconnection process as stronger magnetic fields are convected toward the X-line. When the dissipation region becomes thinner than the ion skin depth, resistive magnetohydrodynamics breaks down as the Hall effect becomes important and the Sweet-Parker solution ceases to exist. A transition to collisionless (Hall) reconnection ensues, increasing the reconnection rate by many orders of magnitude in a very short time. Predictions of the model are consistent with constraints set by observations of solar flares., Comment: 11 pages, 2 figures, submitted to Astrophysical Journal Letters

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., Comment: To be published in JGR

A catastrophe model for the onset of fast magnetic reconnection is presented that suggests why plasma systems with magnetic free energy remain apparently stable for long times and then suddenly release their energy. For a given set of plasma parameters there are generally two stable reconnection solutions: a slow (Sweet-Parker) solution and a fast (Alfv\'enic) Hall reconnection solution. Below a critical resistivity the slow solution disappears and fast reconnection dominates. Scaling arguments predicting the two solutions and the critical resistivity are confirmed with two-fluid simulations., Comment: 10 pages, 4 figures, submitted to Phys. Rev. Lett., Content changed, figure added

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., Comment: 10 pages, 10 figures; accepted by JGR; agu2001.cls and agu.bst included

Large objects, immersed in a homogeneous granular system, migrate when subjected to vibrations. Under certain conditions large heavy objects rise and similar light ones sink to the bottom. This is called reverse buoyancy. We report an experimental study of this singular behavior, for a large sphere immersed in a deep granular bed. A simple mechanism is proposed to describe the motion of a sphere, inside a vertically vibrated granular system. When reverse buoyancy is observed, the measured vertical velocity of the immersed object, as a function of its density, shows a simple behavior. With a one-dimensional mechanical model that takes into account a buoyancy force and the frictional drag, we obtain the rising velocity for heavy objects and the sinking rate for light ones. The model yields a very good qualitative and quantitative agreement with the experiment., Comment: 4 pages, 3 figures