Toward more realistic analytic models of the heliotail: Incorporating magnetic flattening via distortion flows

Both physical arguments and simulations of the global heliosphere indicate that the tailward heliopause is flattened considerably in the direction perpendicular to both the incoming flow and the large-scale interstellar magnetic field. Despite this fact, all of the existing global analytical models of the outer heliosheath's magnetic field assume a circular cross section of the heliotail. To eliminate this inconsistency, we introduce a mathematical procedure by which any analytically or numerically given magnetic field can be deformed in such a way that the cross sections along the heliotail axis attain freely prescribed, spatially dependent values for their total area and aspect ratio. The distorting transformation of this method honors both the solenoidality condition and the stationary induction equation with respect to an accompanying flow field, provided that both constraints were already satisfied for the original magnetic and flow fields prior to the transformation. In order to obtain realistic values for the above parameters, we present the first quantitative analysis of the heliotail's overall distortion as seen in state-of-the-art three-dimensional hybrid MHD-kinetic simulations.Comment: 15 pages, 8 figures. Published in The Astrophysical Journa

Magnetic reconnection solutions in the presence of multiple nulls

It is known that exact analytic solutions can be constructed for incompressible magnetic reconnection in three space dimensions. In the case of an isolated X-point null, there are two types of reconnection solutions, namely, “spine” and “fan” models, which depend on the form of the X-point disturbance. However, such models cannot describe multiple null “separator” reconnection, for which there is independent observational evidence. Here we show that the spine formalism naturally extends to the case of multiple null fields. Solutions showing the characteristics of fan, spine, and separator are described, and a discussion is given of their energy dissipation properties. We demonstrate a family of multiple null, fast reconnection solutions and point out that the classical Sweet-Parker dissipation rate is the slowest that can be achieved with the present models

Hall current effects in dynamic magnetic reconnection solutions

The impact of Hall current contributions on flow driven planar magnetic merging solutions is discussed. The Hall current is important if the dimensionless Hall parameter (or normalized ion skin depth) satisfies cH>η where η is the inverse Lundquist number for the plasma. A dynamic analysis of the problem shows, however, that the Hall current initially manifests itself, not by modifying the planar reconnection field, but by inducing a non-reconnecting perpendicular "separator" component in the magnetic field. Only if the stronger condition c2/H > η is satisfied can Hall currents be expected to affect the planar merging. These analytic predictions are then tested by performing a series of numerical experiments in periodic geometry, using the full system of planar magnetohydrodynamic (MHD) equations. The numerical results confirm that the nature of the merging changes dramatically when the Hall coupling satisfies c2/H > η. In line with the analytic treatment of sheared reconnection, the coupling provided by the Hall term leads to the emergence of multiple current layers that can enhance the global Ohmic dissipation at the expense of the reconnection rate. However, the details of the dissipation depend critically on the symmetries of the simulation, and when the merging is "head-on" (i.e., comprises fourfold symmetry) the reconnection rate can be enhanced

Comparing various multi-component global heliosphere models

Modeling of the global heliosphere seeks to investigate the interaction of the solar wind with the partially ionized local interstellar medium. Models that treat neutral hydrogen self-consistently and in great detail, together with the plasma, but that neglect magnetic fields, constitute a sub-category within global heliospheric models. There are several different modeling strategies used for this sub-category in the literature. Differences and commonalities in the modeling results from different strategies are pointed out. Plasma-only models and fully self-consistent models from four research groups, for which the neutral species is modeled with either one, three, or four fluids, or else kinetically, are run with the same boundary parameters and equations. They are compared to each other with respect to the locations of key heliospheric boundary locations and with respect to the neutral hydrogen content throughout the heliosphere. In many respects, the models' predictions are similar. In particular, the locations of the termination shock agree to within 7% in the nose direction and to within 14% in the downwind direction. The nose locations of the heliopause agree to within 5%. The filtration of neutral hydrogen from the interstellar medium into the inner heliosphere, however, is model dependent, as are other neutral results including the hydrogen wall. These differences are closely linked to the strength of the interstellar bow shock. The comparison also underlines that it is critical to include neutral hydrogen into global heliospheric models.Comment: 10 pages, 4 figures, submitted to a special section at A&A of an ISSI team "Determination of the physical Hydrogen parameters of the LIC from within the Heliosphere

Properties of the Acceleration Regions in Several Loop-structured Solar Flares

Using {\em RHESSI} hard X-ray imaging spectroscopy observations, we analyze electron flux maps for a number of extended coronal loop flares. For each event, we fit a collisional model with an extended acceleration region to the observed variation of loop length with electron energy $E$, resulting in estimates of the plasma density in, and longitudinal extent of, the acceleration region. These quantities in turn allow inference of the number of particles within the acceleration region and hence the filling factor $f$ -- the ratio of the emitting volume to the volume that encompasses the emitting region(s). We obtain values of $f$ that lie mostly between 0.1 and 1.0; the (geometric) mean value is $f = 0.20 \times \div 3.9$, somewhat less than, but nevertheless consistent with, unity. Further, coupling information on the number of particles in the acceleration region with information on the total rate of acceleration of particles above a certain reference energy (obtained from spatially-integrated hard X-ray data) also allows inference of the specific acceleration rate (electron s$^{-1}$ per ambient electron above the chosen reference energy). We obtain a (geometric) mean value of the specific acceleration rate $\eta(20$ keV) $= (6.0 \times / \div 3.4) \times 10^{-3}$ electrons s$^{-1}$ per ambient electron; this value has implications both for the global electrodynamics associated with replenishment of the acceleration region and for the nature of the particle acceleration process

Exclusion of Tiny Interstellar Dust Grains from the Heliosphere

The distribution of interstellar dust grains (ISDG) observed in the Solar System depends on the nature of the interstellar medium-solar wind interaction. The charge of the grains couples them to the interstellar magnetic field (ISMF) resulting in some fraction of grains being excluded from the heliosphere while grains on the larger end of the size distribution, with gyroradii comparable to the size of the heliosphere, penetrate the termination shock. This results in a skewing the size distribution detected in the Solar System. We present new calculations of grain trajectories and the resultant grain density distribution for small ISDGs propagating through the heliosphere. We make use of detailed heliosphere model results, using three-dimensional (3-D) magnetohydrodynamic/kinetic models designed to match data on the shape of the termination shock and the relative deflection of interstellar neutral H and He flowing into the heliosphere. We find that the necessary inclination of the ISMF relative to the inflow direction results in an asymmetry in the distribution of the larger grains (0.1 micron) that penetrate the heliopause. Smaller grains (0.01 micron) are completely excluded from the Solar System at the heliopause.Comment: 5 pages, 5 figures, accepted for publication in the Solar Wind 12 conference proceeding

Solar Particle Acceleration at Reconnecting 3D Null Points

Context: The strong electric fields associated with magnetic reconnection in solar flares are a plausible mechanism to accelerate populations of high energy, non-thermal particles. One such reconnection scenario occurs at a 3D magnetic null point, where global plasma flows give rise to strong currents in the spine axis or fan plane. Aims: To understand the mechanism of charged particle energy gain in both the external drift region and the diffusion region associated with 3D magnetic reconnection. In doing so we evaluate the efficiency of resistive spine and fan models for particle acceleration, and find possible observables for each. Method: We use a full orbit test particle approach to study proton trajectories within electromagnetic fields that are exact solutions to the steady and incompressible magnetohydrodynamic equations. We study single particle trajectories and find energy spectra from many particle simulations. The scaling properties of the accelerated particles with respect to field and plasma parameters is investigated. Results: For fan reconnection, strong non-uniform electric drift streamlines can accelerate the bulk of the test particles. The highest energy gain is for particles that enter the current sheet, where an increasing "guide field" stabilises particles against ejection. The energy is only limited by the total electric potential energy difference across the fan current sheet. The spine model has both slow external electric drift speed and weak energy gain for particles reaching the current sheet. Conclusions: The electromagnetic fields of fan reconnection can accelerate protons to the high energies observed in solar flares, gaining up to 0.1 GeV for anomalous values of resistivity. However, the spine model, which gave a harder energy spectrum in the ideal case, is not an efficient accelerator after pressure constraints in the resistive model are included.Comment: 15 pages, 14 figures. Submitted to Astronomy and Astrophysic

Proton acceleration in analytic reconnecting current sheets

Particle acceleration provides an important signature for the magnetic collapse that accompanies a solar flare. Most particle acceleration studies, however, invoke magnetic and electric field models that are analytically convenient rather than solutions of the governing magnetohydrodynamic equations. In this paper a self-consistent magnetic reconnection solution is employed to investigate proton orbits, energy gains, and acceleration timescales for proton acceleration in solar flares. The magnetic field configuration is derived from the analytic reconnection solution of Craig and Henton. For the physically realistic case in which magnetic pressure of the current sheet is limited at small resistivities, the model contains a single free parameter that specifies the shear of the velocity field. It is shown that in the absence of losses, the field produces particle acceleration spectra characteristic of magnetic X-points. Specifically, the energy distribution approximates a power law ~ξ-3/2 nonrelativistically, but steepens slightly at the higher energies. Using realistic values of the “effective” resistivity, we obtain energies and acceleration times that fall within the range of observational data for proton acceleration in the solar corona