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

Dispersive Fast Magnetosonic Waves and Shock‐Driven Compressible Turbulence in the Inner Heliosheath

The solar wind in the inner heliosheath beyond the termination shock (TS) is a nonequilibrium collisionless plasma consisting of thermal solar wind ions, suprathermal pickup ions, and electrons. In such multi‐ion plasma, two fast magnetosonic wave modes exist, the low‐frequency fast mode and the high‐frequency fast mode. Both fast modes are dispersive on fluid and ion scales, which results in nonlinear dispersive shock waves. We present high‐resolution three‐fluid simulations of the TS and the inner heliosheath up to a few astronomical units (AU) downstream of the TS. We show that downstream propagating nonlinear fast magnetosonic waves grow until they steepen into shocklets, overturn, and start to propagate backward in the frame of the downstream propagating wave. The counterpropagating nonlinear waves result in 2‐D fast magnetosonic turbulence, which is driven by the ion‐ion hybrid resonance instability. Energy is transferred from small scales to large scales in the inverse cascade range, and enstrophy is transferred from large scales to small scales in the direct cascade range. We validate our three‐fluid simulations with in situ high‐resolution Voyager 2 magnetic field observations in the inner heliosheath. Our simulations reproduce the observed magnetic turbulence spectrum with a spectral slope of −5/3 in frequency domain. However, the fluid‐scale turbulence spectrum is not a Kolmogorov spectrum in wave number domain because Taylor’s hypothesis breaks down in the inner heliosheath. The magnetic structure functions of the simulated and observed turbulence follow the Kolmogorov‐Kraichnan scaling, which implies self‐similarity.Key PointsNonlinear dispersive fast magnetosonic waves produce 2‐D compressible turbulence downstream of the termination shockTaylor’s hypothesis breaks down in the subfast magnetosonic solar wind in the inner heliosheathThe magnetic turbulence spectrum observed by Voyager 2 in the inner heliosheath is reproduced by self‐consistent three‐fluid MHD simulationPeer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/163373/2/jgra56004_am.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/163373/1/jgra56004.pd

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

Galactic Cosmic Rays Modulation in the Vicinity of Corotating Interaction Regions: Observations During the Last Two Solar Minima

Corotating interaction regions (CIRs) are responsible for short-term recurrent cosmic-ray modulation, prominent near solar minima. Using the OMNI data sets for two periods of low solar activity near the beginning and end of solar cycle 24, superposed epoch analysis was performed on the solar wind plasma features for 53 and 43 events during periods 2007–2008 and 2017–2018, respectively. Turbulent properties of the solar wind were studied using the variance method for each CIR. Power spectra have been constructed for overlapped subintervals in the vicinity of stream interfaces (SIs). Using measured correlation lengths and turbulent energies, parallel and perpendicular diffusion mean free paths for cosmic-ray ions have been inferred based on two distinct theoretical formulations. For the two periods with opposite solar polarities, our results show that unlike solar wind speed, magnetic field strength, flow pressure, and proton density are relatively higher during the latest period. Increased turbulent energy and reduced parallel transport coefficients of energetic particles at the SIs are observed. The diffusion coefficients follow the same trends during both periods. The perpendicular diffusion starts increasing nearly a day before SIs and is higher in the fast wind. Superposed epoch analysis is performed on the >120 MeV proton count rate obtained from the CRIS instrument on board the ACE spacecraft for the same events. The recorded proton rates have peaks half a day before a SI and reach their minimum more than a day after a SI and have a high anticorrelation with the perpendicular diffusion coefficient