Termination Shock Asymmetries as Seen by the Voyager Spacecraft: The Role of the Interstellar Magnetic Field and Neutral Hydrogen

We show that asymmetries of the termination shock due to the influence of the interstellar magnetic field (ISMF) are considerably smaller in the presence of neutral hydrogen atoms, which tend to symmetrize the heliopause, the termination shock, and the bow shock due to charge exchange with charged particles. This leads to a much stronger restriction on the ISMF direction and its strength. We demonstrate that in the presence of the interplanetary magnetic field the plane defined by the local interstellar medium (LISM) velocity and magnetic field vectors does not exactly coincide with the plane defined by the interstellar neutral helium and hydrogen velocity vectors in the supersonic solar wind region, which limits the accuracy of the inferred direction of the ISMF. We take into account the tilt of the LISM velocity vector with respect to the ecliptic plane and show that magnetic fields as strong as 3 μG or greater may be necessary to account for the observed asymmetry. Estimates are made of the longitudinal streaming anisotropy of energetic charged particles at the termination shock caused by the nonalignment of the interplanetary magnetic field with its surface. By investigating the behavior of interplanetary magnetic field lines that cross the Voyager 1 trajectory in the inner heliosheath, we estimate the length of the trajectory segment that is directly connected by these lines to the termination shock. A possible effect of the ISMF draping over the heliopause is discussed in connection with radio emission generated in the outer heliosheath

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

A Two-Dimensional, Self-Consistent Model of Galactic Cosmic Rays in the Heliosphere

We present initial results from our new two-dimensional (radius and latitude), self-consistent model of galactic cosmic rays in the heliosphere. We focus on the latitudinal variations in the solar wind flow caused by the energetic particles. Among other things our results show that the cosmic rays significantly modify the latitudinal structure of the solar wind flow downstream of the termination shock. Specifically, for A>0 (corresponding to the present solar minimum) the wind beyond the shock is driven towards the equator, resulting in a faster wind flow near the current sheet, while for A<0 the effect is reversed and the wind turns towards the pole, with a faster flow at high latitudes. We attribute this effect to the latitudinal gradients in the cosmic ray pressure, caused by drifts, that squeeze the flow towards the ecliptic plane or the pole, respectively.Comment: 10 pages, 4 Postscript figures, uses AAS LaTeX v4.0, to be published in The Astrophysical Journal Letter

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

Spectral Evolution of Anomalous Cosmic Rays at Voyager 1 beyond the Termination Shock

When the Voyager 1 spacecraft crossed the termination shock (TS) on 2004 December 16, the energy spectra of anomalous cosmic rays (ACRs) could not have been produced by steady-state diffusive shock acceleration. However, over the next few years, in the declining phase of the solar cycle, the spectra began to evolve into the expected power-law profile. Observations at the shock led to a broad range of alternative theories for ACR acceleration. In spite of that, in this work we show that the observations could be explained by assuming ACRs are accelerated at the TS. In this paper, we propose that the solar cycle had an important effect on the unrolling of the spectra in the heliosheath. To investigate the spectral evolution of ACRs, a magnetohydrodynamic background model with stationary solar-wind inner boundary conditions was used to model the transport of helium and oxygen ions. We used a backward-in-time stochastic integration technique where phase-space trajectories are integrated until the so-called "injection energy" is reached. Our simulation results were compared with Voyager 1 observations using three different diffusion models. It is shown that the spectral evolution of ACRs in the heliosheath at Voyager 1 could be explained by an increase in the source strength and an enhancement in diffusion as a result of a decrease of the turbulent correlation length in the declining phase of the solar cycle. At the same time, drift effects seem to have had a smaller effect on the evolution of the spectra

Cosmic rays beyond the boundary of the heliosphere

In August of 2012 the Voyager 1 space probe has left the solar-wind bubble of ionized gas we call the heliosphere and entered the denser and colder environment of the interstellar cloud surrounding the solar system. Energetic charged particles underwent dramatic changes past the heliopause: the heliospheric ions disappeared completely, while the galactic cosmic rays were for the first time measured in their unmodulated state. The interstellar medium turned out to be almost entirely devoid of turbulent magnetic fluctuations, therefore the transport of cosmic rays is governed by a large-scale geometry of the magnetic field. We discuss observations of heliospheric ions, including anomalous cosmic rays, near the heliopause transition, and propose interpretations of the measured intensities and pitch-angle distributions based on gradient drift in a weakly nonuniform magnetic field. The heliopause transition appears to be permeated by magnetic flux tubes connected to the interstellar space and facilitating particle escape. These flux tubes may be a product of interchange instability driven by a plasma pressure gradient across the heliopause. The curvature of magnetic field lines and the anti-sunward gradient in plasma kinetic pressure provide conditions favorable for an interchange. The two flux tube crossings by the spacecraft allowed an indirect measurement of the plasma radial velocity near the heliopause