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

The scattering of electrons by heat-flux-driven whistler waves is explored with a particle-in-cell (PIC) simulation relevant to the transport of energetic electrons in flares. The simulation is initiated with a large heat flux that is produced using a kappa distribution of electrons with positive velocity and a cold return current beam. This system represents energetic electrons escaping from a reconnection-driven energy release site. This heat flux system drives large amplitude oblique whistler waves propagating both along and against the heat flux, as well as electron acoustic waves. While the waves are dominantly driven by the low energy electrons, including the cold return current beam, the energetic electrons resonate with and are scattered by the whistlers on time scales of the order of a hundred electron cyclotron times. Peak whistler amplitudes of $\tilde{B} / B_{0} \sim 0.125$ and angles of $\sim 60 \degree$ with respect to the background magnetic field are observed. Electron perpendicular energy is increased while the field-aligned electron heat flux is suppressed. The resulting scattering mean-free-paths of energetic electrons are small compared with the typical scale size of energy release sites in flares, which might lead to the effective confinement of energetic electrons that is required for the production of very energetic particles.

The role of turbulence in low-$\beta$, 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\'en 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\'en speed and is expected to be stable in the low $\beta$ solar corona. The results suggest that in a guide field reconnection exhaust with $1 \gg \beta > 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.

Kinetic Riemann simulations have been completed to explore particle heating during guide field reconnection in the low-$\beta$ environment of the inner heliosphere and the solar corona. The reconnection exhaust is bounded by two rotational discontinuities (RD) and two slow shocks (SS) form within the exhaust as in magnetohydrodynamic (MHD) models. At the RDs, ions are accelerated by the magnetic field tension to drive the reconnection outflow as well as flows in the out-of-plane direction. The out-of-plane flows stream toward the midplane and meet to drive the SSs. The SSs differ greatly from those in the MHD model. The turbulence at the shock fronts and both upstream and downstream is weak so the shocks are laminar and produce little dissipation. Downstream of the SSs the counterstreaming ion beams lead to higher density, which leads to a positive potential between the SSs that acts to confine the downstream electrons to maintain charge neutrality. The potential accelerates electrons from upstream of the SSs to downstream region and traps a small fraction but only modestly increases the downstream electron temperature above the upstream value. In the low-$\beta$ limit the released magnetic energy is split between bulk flow and ion heating with little energy going to electrons. That the model does not produce strong electron heating nor an energetic electron component suggests that other mechanisms, such as multiple x-line reconnection, are required to explain energetic electron production in large flares. The model can be tested with the expected data from the Parker Solar Probe.

A new computational model is presented suitable for exploring the self-consistent production of energetic electrons during magnetic reconnection in macroscale systems. The equations are based on the recent discovery that parallel electric fields are ineffective drivers of energetic particles during reconnection so that the kinetic scales which control the development of such fields can be ordered out of the equations. The resulting equations consist of a magnetohydrodynamic (MHD) backbone with the energetic component represented by macro-particles described by the guiding center equations. Crucially, the energetic component feeds back on the MHD equations so that the total energy of the MHD fluid and the energetic particles is conserved. The equations correctly describe the firehose instability, whose dynamics plays a key role in throttling reconnection and in controlling the spectra of energetic particles. The results of early tests of the model, including the propagation of Alfven waves in a system with pressure anisotropy and the growth of firehose modes, establish that the basic algorithm is stable and produces reliable physics results in preparation for further benchmarking with particle-in-cell models of reconnection.

Heat flux suppression in collisionless plasmas for a large range of plasma $\beta$ is explored using two-dimensional particle-in-cell simulations with a strong, sustained thermal gradient. We find that a transition takes place between whistler-dominated (high-$\beta$) and double-layer-dominated (low-$\beta$) heat flux suppression. Whistlers saturate at small amplitude in the low beta limit and are unable to effectively suppress the heat flux. Electrostatic double layers suppress the heat flux to a mostly constant factor of the free streaming value once this transition happens. The double layer physics is an example of ion-electron coupling and occurs on a scale of roughly the electron Debye length. The scaling of ion heating associated with the various heat flux driven instabilities is explored over the full range of $\beta$ explored. The range of plasma-$\beta$s studied in this work makes it relevant to the dynamics of a large variety of astrophysical plasmas, including the intracluster medium of galaxy clusters, hot accretion flows, stellar and accretion disk coronae, and the solar wind.

The dynamics of weakly magnetized collisionless plasmas in the presence of an imposed temperature gradient along an ambient magnetic field is explored with particle-in-cell simulations and modeling. Two thermal reservoirs at different temperatures drive an electron heat flux that destabilizes off-angle whistler-type modes. The whistlers grow to large amplitude, $\delta B / B_{0} \simeq 1$, and resonantly scatter the electrons, significantly reducing the heat flux. A surprise is that the resulting steady state heat flux is largely independent of the thermal gradient. The rate of thermal conduction is instead controlled by the finite propagation speed of the whistlers, which act as mobile scattering centers that convect the thermal energy of the hot reservoir. The results are relevant to thermal transport in high $\beta$ astrophysical plasmas such as hot accretion flows and the intracluster medium of galaxy clusters.

Data from the NASA Magnetospheric Multiscale (MMS) mission are used to investigate asymmetric magnetic reconnection at the dayside boundary between the Earth's magnetosphere and the solar wind (the magnetopause). High-resolution measurements of plasmas, electric and magnetic fields, and waves are used to identify highly localized (~15 electron Debye lengths) standing wave structures with large electric-field amplitudes (up to 100 mV/m). These wave structures are associated with spatially oscillatory dissipation, which appears as alternatingly positive and negative values of J dot E (dissipation). For small guide magnetic fields the wave structures occur in the electron stagnation region at the magnetosphere edge of the EDR. For larger guide fields the structures also occur near the reconnection x-line. This difference is explained in terms of channels for the out-of-plane current (agyrotropic electrons at the stagnation point and guide-field-aligned electrons at the x-line).

Magnetic reconnection is an important driver of energetic particles in many astrophysical phenomena. Using kinetic particle-in-cell (PIC) simulations, we explore the impact of three-dimensional reconnection dynamics on the efficiency of particle acceleration. In two-dimensional systems, Alfv\'enic outflows expel energetic electrons into flux ropes where they become trapped and disconnected from acceleration regions. However, in three-dimensional systems these flux ropes develop axial structure that enables particles to leak out and return to acceleration regions. This requires a finite guide field so that particles may move quickly along the flux rope axis. We show that greatest energetic electron production occurs when the guide field is of the same order as the reconnecting component: large enough to facilitate strong transport, but not so large as to throttle the dominant Fermi mechanism responsible for efficient electron acceleration. This suggests a natural explanation for the envelope of electron acceleration during the impulsive phase of eruptive flares.

Based on the difference between the orientation of the interstellar $B_{ISM}$ and the solar magnetic fields, there was an expectation that the magnetic field direction would rotate dramatically across the heliopause (HP). However, the Voyager 1 spacecraft measured very little rotation across the HP. Previously we showed that the $B_{ISM}$ twists as it approaches the HP and acquires a strong T component (East-West). Here we establish that reconnection in the eastern flank of the heliosphere is responsible for the twist. On the eastern flank the solar magnetic field has twisted into the positive N direction and reconnects with the Southward pointing component of the $B_{ISM}$. Reconnection drives a rotational discontinuity (RD) that twists the $B_{ISM}$ into the -T direction and propagates upstream in the interstellar medium towards the nose. The consequence is that the N component of $B_{ISM}$ is reduced in a finite width band upstream of the HP. Voyager 1 currently measures angles ($\delta=sin^{-1}(B_{N}/B)$) close to solar values. We present MHD simulations to support this scenario, suppressing reconnection in the nose region while allowing it in the flanks, consistent with recent ideas about reconnection suppression from diamagnetic drifts. The jump in plasma $\beta$ (the plasma to magnetic pressure) across the nose of HP is much greater than in the flanks because the heliosheath $\beta$ is greater there than in the flanks. Large-scale reconnection is therefore suppressed in the nose but not at the flanks. Simulation data suggest that $B_{ISM}$ will return to its pristine value $10-15~AU$ past the HP.

The misalignment of the solar rotation axis and the magnetic axis of the Sun produces a periodic reversal of the Parker spiral magnetic field and the sectored solar wind. The compression of the sectors is expected to lead to reconnection in the heliosheath (HS). We present particle-in-cell simulations of the sectored HS that reflect the plasma environment along the Voyager 1 and 2 trajectories, specifically including unequal positive and negative azimuthal magnetic flux as seen in the Voyager data \citep{Burlaga03}. Reconnection proceeds on individual current sheets until islands on adjacent current layers merge. At late time bands of the dominant flux survive, separated by bands of deep magnetic field depletion. The ambient plasma pressure supports the strong magnetic pressure variation so that pressure is anti-correlated with magnetic field strength. There is little variation in the magnetic field direction across the boundaries of the magnetic depressions. At irregular intervals within the magnetic depressions are long-lived pairs of magnetic islands where the magnetic field direction reverses so that spacecraft data would reveal sharp magnetic field depressions with only occasional crossings with jumps in magnetic field direction. This is typical of the magnetic field data from the Voyager spacecraft \citep{Burlaga11,Burlaga16}. Voyager 2 data reveals that fluctuations in the density and magnetic field strength are anti-correlated in the sector zone as expected from reconnection but not in unipolar regions. The consequence of the annihilation of subdominant flux is a sharp reduction in the "number of sectors" and a loss in magnetic flux as documented from the Voyager 1 magnetic field and flow data \citep{Richardson13}.

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.

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, $\beta=8\pi 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.

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.

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.

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.

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

A procedure for generating random variates from a relativistic Maxwellian distribution with arbitrary temperature and drift velocity is presented. The algorithm is based on the rejection method and can be used to initialize particle velocities in kinetic simulations of plasmas and gases.

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.

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.

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.

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.

The acceleration of solar flare ions during magnetic reconnection is explored via particle-in-cell simulations that self-consistently follow the motions of both protons and $\alpha$ particles. We demonstrate that the dominant ion heating during reconnection with a guide field (a magnetic component perpendicular to the reconnection plane) results from pickup behavior during the entry into reconnection exhausts. In contrast with anti-parallel reconnection, the temperature increment is dominantly transverse, rather than parallel, to the local magnetic field. The comparison of protons and alphas reveals a mass-to-charge ($M/Q$) threshold in pickup behavior that favors heating of high $M/Q$ ions over protons, which is consistent with impulsive flare observations.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

The nature of collisionless reconnection in a three-species plasma composed of a heavy species, protons, and electrons is examined. Besides the usual two length scales present in two-species reconnection, there are two additional larger length scales in the system: one associated with a "heavy whistler" which produces a large scale quadrupolar out-of-plane magnetic field, and one associated with the "heavy Alfven" wave which can slow the outflow speed and thus the reconnection rate. The consequences for reconnection in the magnetotail with an O+ population present are discussed.

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.

We describe an approach for finding the eigenfrequencies of solar acoustic modes (p modes) in a convective envelope in the WKB limit. This approximation restricts us to examining the effects of fluid motions which are large compared to the mode wavelength, but allows us to treat the three-dimensional mode as a localized ray. The method of adiabatic switching is then used to investigate the frequency shifts resulting from simple perturbations to a polytropic model of the convection zone as well as from two basic models of a convective cell. We find that although solely depth-dependent perturbations can give frequency shifts which are first order in the strength of the perturbation, models of convective cells generate downward frequency shifts which are second order in the perturbation strength. These results may have implications for resolving the differences between eigenfrequencies derived from solar models and those found from helioseismic observations.