Outer Heliosheath Research Articles

Constraining the Properties of the Multicomponent Local Interstellar Medium: MHD-kinetic Modeling Validated by Voyager and New Horizons Data

We introduce the first solar-cycle simulations from our 3D, global MHD-plasma/kinetic-neutrals model, where both hydrogen and helium atoms are treated kinetically, while electrons and helium ions are described as individual fluids. Using Voyager/PWS observations of electron density up to 160 au from the Sun for validation of several different global models, we conclude that the current estimates for the proton density in the local interstellar medium (LISM) need a revision. Our findings indicate that the commonly accepted value of 0.054 cm−3 may need to be increased to values exceeding 0.07 cm−3. We also show how different assumptions regarding the proton velocity distribution function in the outer heliosheath may affect the global solution. A new feature revealed by our simulations is that the helium ion flow may be significantly compressed and heated in the heliotail at heliocentric distances exceeding ∼400 au. Additionally, we identify a Kelvin–Helmholtz instability at the boundary of the slow and fast solar wind in the inner heliosheath, which acts as a driver of turbulence in the heliotail. These results are crucial for inferring the properties of the LISM and of the global heliosphere structure.

Interstellar Neutral Hydrogen in the Heliosphere: New Horizons Observations in the Context of Models

Interstellar neutral (ISN) hydrogen is the most abundant species in the outer heliosheath and the very local interstellar medium (VLISM). Charge-exchange collisions in the outer heliosheath result in filtration, reducing the ISN hydrogen density inside the heliosphere. Additionally, these atoms are intensively ionized close to the Sun, resulting in a substantial reduction of their density within a few astronomical units from the Sun. The products of this ionization—pickup ions (PUIs)—are detected by charged particle detectors. The Solar Wind Around Pluto instrument on New Horizons provides, for the first time, PUI observations from the distant heliosphere. We analyze the observations collected between 22 and 52 au from the Sun to find the ISN hydrogen density profile and compare the results with predictions from global heliosphere models. We conclude that the density profile derived from the observations is inconsistent with steady-state model predictions. This discrepancy is not explained by time variations close to the Sun and thus may be related to the temporal evolution of the outer boundaries or VLISM conditions. Furthermore, we show that the cold and hot models of ISN hydrogen distribution are not a good approximation closer to the termination shock. Therefore, we recommend a new fiduciary point based on the available New Horizons observations at 40 au from the Sun, at ecliptic direction (285.°62, 1.°94), where the ISN hydrogen density is 0.11 cm−3. The continued operation of New Horizons should give better insight into the source of the discussed discrepancy.

The unexpected role of heliospheric boundaries in facilitating interstellar dust penetration at 1–5 AU

Aims. Interstellar dust (ISD) particles penetrate the heliosphere because of the relative motion of the local interstellar cloud and the Sun. The penetrated particles pass through the heliospheric interface, that is, the region in which solar wind and interstellar plasma interact. As a result, the ISD flow is modified after the passage through this region under the influence of electromagnetic force. The main goal of this work is to show how the heliospheric interface affects the distribution of ISD particles near the Sun. Methods. We have developed a Monte Carlo model of the ISD distribution in the heliosphere. It first takes the effects of the heliospheric interface and the rotating heliospheric current sheet into account. The effects of the heliospheric interface were probed using a global heliospheric model. Results. The computation results show that the heliospheric interface strongly influences the distribution of relatively small (radius a = 150 − 250 nm) astronomical silicates. The unexpected finding is that the heliospheric interface facilitates the penetration of a = 150 nm particles at small heliocentric distances and, particularly, to the Ulysses orbit (1 − 5 AU). We demonstrate that the deflection of ISD particles in the outer heliosheath is the principal mechanism that causes the effects of the heliospheric interface on the distribution near the Sun. The computations with different heliospheric models show that the distribution near the Sun is sensitive to the plasma parameters in the pristine local interstellar medium. Thus, we demonstrated that being measured near the Sun, the ISD may serve as a new independent diagnostics of the local interstellar medium and the heliospheric boundaries.

Magnetic Trapping of Galactic Cosmic Rays in the Outer Heliosheath and Their Preferential Entry into the Heliosphere

This paper examines the geometry of interstellar magnetic field lines close to the boundary of the heliosphere in the direction of the unperturbed local interstellar magnetic field, where the field lines are spread apart by the heliopause (HP). Such field parting establishes a region of weaker magnetic field of about 300 au in size in the northern hemisphere that acts as a giant magnetic trap affecting the propagation of galactic cosmic rays (GCRs). The choice of an analytic model of the magnetic field in the very local interstellar medium allows us to qualitatively study the resulting magnetic field draping pattern while avoiding unphysical dissipation across the HP-impeding numerical magnetohydrodynamic (MHD) models. We investigate GCR transport in the region exterior to the heliosphere, including the magnetic trap, subject to guiding center drifts, pitch angle scattering, and perpendicular diffusion. The transport coefficients were derived from Voyager 1 observations of magnetic turbulence in the VLISM. Our results predict a ring current of energetic ions drifting around the interior of the magnetic trap. It is also demonstrated that GCRs cross the HP for the first time preferentially through a crescent-shaped region between the magnetic trap and the upwind direction. The paper includes results of MHD modeling of the heliosphere that provide the coordinates of the center of the magnetic trap in ecliptic coordinates. In addition to the heliosphere, we examine several extreme field draping configurations that could describe the astrospheres of other stars.

The Effect of Angular Scattering Imposed by Charge Exchange and Elastic Collisions on Interstellar Neutral Hydrogen Atoms

Angular scattering (AS) in charge exchange and elastic collisions between interstellar ions and neutral (ISN) atoms has been assumed to be negligible in previous studies. Here, we investigate the momentum transfer associated with the AS of H atoms using Monte Carlo calculations to simulate their transport through the outer heliosheath. We consider two cases where charge exchange and elastic collisions between ISN H atoms and protons occur with and without AS in the outer heliosheath. We show that considering AS decelerates and heats primary ISN H, reducing the effect of selective charge exchange in the outer heliosheath. Secondary ISN H atoms, on the other hand, are not significantly affected by AS. We then simulate the transport of ISN H atoms inside the heliosphere to simulate count rates observed in the lowest energy bin of IBEX-Lo. We study the effect of radiation pressure on the ISN H measurements for the cases with and without AS and compare them with our previous findings. We find an effective radiation parameter (μ eff, which represents force associated with radiation pressure relative to gravity) for the years 2009–2018 based on the longitudinal shift of the ISN H signal. The two cases with and without AS reproduce the longitudinal shift in accordance with variations in solar activity, in agreement with our previous results, and they result in similar values for the μ eff, which is ∼22%–23% larger than the value found based on solar Lyα profile observations.

The Impact of Pickup Ion Thermal Spread on Pickup Ion Ring-beam-driven Instabilities and Scattering in the Outer Heliosheath

The present study investigates the unstable waves driven by the pickup ions of ring-beam distributions with various pickup angles and pickup ion temperatures in the outer heliosheath, using both linear instability analysis and hybrid simulations. While previous papers have generally assumed specific pickup ion temperatures, this study takes a more comprehensive approach by examining a wide range of pickup ion temperatures that may occur in the outer heliosheath. Our one-dimensional simulations demonstrate that regardless of the initial pickup ion temperature, the pitch-angle scattering of the pickup ions at small pickup angles saturates before the pickup ions can reach the hemisphere of negative parallel velocities with respect to the background magnetic field in velocity space, while at a 90° pickup angle, the pickup ions of ring distributions attain a significant level of isotropy. In contrast, the two-dimensional simulation results show that at all pickup angles, increasing the pickup ion temperature only slightly reduces the pickup ion scattering rate, but does not prevent the pickup ions from reaching the hemisphere of negative parallel velocities. Overall, the results do not align with the requirements of the earlier proposed scenarios for the secondary energetic neutral atom (ENA) mechanism of the Interstellar Boundary Explorer ENA ribbon, which demand either weak pickup ion scattering in the outer heliosheath or at least incomplete pickup ion scattering in the off-ribbon directions.

Science Opportunities for IMAP-Lo Observations of Interstellar Neutral Helium, Neon, and Oxygen during a Maximum of Solar Activity

Direct-sampling observations of interstellar neutral (ISN) species and their secondary populations give information about the physical state of the local interstellar medium and processes occurring in the outer heliosheath. Such observations are performed from Earth’s orbit by the IBEX-Lo experiment on board the Interstellar Boundary Explorer (IBEX) mission. IBEX ISN viewing is restricted to directions close to perpendicular to the Earth–Sun line, which limits the observations of interstellar species to several months during the year. A greatly improved data set will be possible for the upcoming Interstellar Mapping and Acceleration Probe (IMAP) mission due to a novel concept of putting the IMAP ISN detector, IMAP-Lo, on a pivot platform that varies the angle of observation relative to the Sun–Earth line and the detector boresight. Here, we suggest a 2 yr scenario for varying the viewing angle in such a way that all the necessary atom components can be observed sufficiently well to achieve the science goals of the nominal IMAP mission. This scenario facilitates, among others, removal of the correlation of the inflow parameters of interstellar gas, unambiguous analysis of the primary and secondary populations of interstellar helium (He), neon (Ne), and oxygen (O), and determination of the ionization rates of He and Ne free of possible calibration bias. The scheme is operationally simple, provides good counting statistics, and synergizes observations of interstellar species and heliospheric energetic neutral atoms.

Interstellar Conditions Deduced from Interstellar Neutral Helium Observed by IBEX and Global Heliosphere Modeling

In situ observations of interstellar neutral (ISN) helium atoms by the IBEX-Lo instrument on board the Interstellar Boundary Explorer (IBEX) mission are used to determine the velocity and temperature of the pristine very local interstellar medium (VLISM). Most ISN helium atoms penetrating the heliosphere, known as the primary population, originate in the pristine VLISM. As the primary atoms travel through the outer heliosheath, they charge exchange with He+ ions in slowed and compressed plasma, creating the secondary population. With more than 2.4 million ISN helium atoms being sampled by IBEX during ISN seasons 2009–2020, we compare the observations with the predictions of a parameterized model of ISN helium transport in the heliosphere. We account for the filtration of ISN helium atoms at the heliospheric boundaries by charge-exchange and elastic collisions. We examine the sensitivity of the ISN helium fluxes to the interstellar conditions described by the pristine VLISM velocity, temperature, magnetic field, and composition. We show that comprehensive modeling of the filtration processes is critical for interpreting ISN helium observations, as the change in the derived VLISM conditions exceeds the statistical uncertainties when accounting for these effects. The pristine VLISM parameters found by this analysis are the flow speed (26.6 km s−1), inflow direction in ecliptic coordinates (255.°7, 5.°04), temperature (7350 K), and B − V plane inclination to the ecliptic plane (53.°7). The derived pristine VLISM He+ density is 9.7 × 10−3 cm−3. Additionally, we show a strong correlation between the interstellar plasma density and magnetic field strength deduced from these observations.

The Direction of the Flow of Interstellar Neutral H Based on Photometric Observations from SOHO/SWAN

Interstellar neutral hydrogen flows into the heliosphere as a mixture of the primary and secondary populations from two somewhat different directions due to splitting occurring in the magnetized outer heliosheath. The direction of the inflow of interstellar neutral H observed in the inner heliosphere, confronted with that of the unperturbed flow of interstellar neutral helium, is important for understanding the geometry of the distortion of the heliosphere from axial symmetry. It is also needed for facilitating remote-sensing studies of the solar wind structure based on observations of the helioglow, such as those presently performed by SOHO/SWAN, and in the near future by IMAP/GLOWS. In the past, the only means to measure the direction of the flow of interstellar hydrogen were spectroscopic observations of the helioglow. Here, we propose a new method to determine this parameter based on a long series of photometric observations of the helioglow. The method is based on purely geometric considerations and does not depend on any model and absolute calibration of the measurements. We apply this method to sky maps of the helioglow available from the SOHO/SWAN experiment and derive the mean longitude of the flow of interstellar hydrogen. We obtain 253.°1 ± 2.°8, which is in perfect agreement with the previously obtained results based on spectroscopic observations.

Oblique Electromagnetic Instabilities Driven by Pickup Ion Ring-beam Distributions in the Outer Heliosheath. I. Linear Instability Analysis

The energetic neutral atom ribbon observed by the Interstellar Boundary Explorer spacecraft is believed to originate from the pickup ions in the outer heliosheath. The outer heliosheath pickup ions generally have a ring-beam velocity distribution at a certain pickup angle, α, the angle at which these ions are picked up by the interstellar magnetic field. The pickup ion ring-beam distributions can drive unstable waves of different propagation angles with respect to the background interstellar magnetic field, θ. Previous studies of the outer heliosheath pickup ion dynamics were mainly focused on ring-like pickup ion distributions with α ≈ 90° and/or the parallel- and antiparallel-propagating unstable waves (θ = 0° and 180°). The present study carries out linear kinetic instability analysis to investigate both the parallel and oblique unstable modes (0° ≤ θ ≤ 180°) driven by ring-beam pickup ion distributions of different pickup angles between 0° and 90°. Our linear instability analysis reveals that ring-beam pickup ion distributions can excite oblique mirror waves as well as parallel/quasi-parallel and oblique right- and left-helicity waves. The maximum growth rate among all the instabilities belongs to the parallel-propagating left-helicity waves at most pickup angles. Furthermore, the evolution of the unstable mirror waves by varying pickup angle indicates that as the pickup angle increases, the maximum growth rate of the mirror modes increases, while its propagation angle decreases.

Oblique Electromagnetic Instabilities Driven by Pickup Ion Ring-beam Distributions in the Outer Heliosheath. II. Hybrid Simulations

Hybrid simulations are carried out to study the instabilities driven by the outer heliosheath pickup ions of ring-beam velocity distributions. Following the linear instability analysis in the companion paper, the electromagnetic instabilities at various propagation directions (0° ≤ θ ≤ 180°) are studied for pickup ion ring-beam distributions of different pickup angles (0° ≤ α ≤ 90°). Previous studies of the ring-beam pickup ion distributions have been limited to the parallel- and antiparallel-propagating modes (θ = 0° and 180°). In this paper, several one-dimensional hybrid simulations are performed to investigate the most unstable modes of different propagation angles predicted by the linear instability analysis. Moreover, two-dimensional hybrid simulations are conducted to examine the simultaneous excitation of the parallel, antiparallel, and oblique unstable modes and their role in the scattering of the outer heliosheath pickup ions. Our recent one-dimensional hybrid and particle-in-cell simulations with simulation domains along the background magnetic field demonstrated that at small pickup angles, the outer heliosheath pickup ions experience an incomplete scattering limited to the hemisphere of positive parallel velocities with respect to the background magnetic field. In contrast, the present two-dimensional hybrid simulation results show that when both parallel/antiparallel and oblique unstable modes are acting together, the pickup ions can be scattered into the hemisphere of negative parallel velocities and approach a near-isotropic velocity distribution. This is in conflict with the spatial retention scenario of the energetic neutral atom ribbon observed by the Interstellar Boundary Explorer, which requires pickup ions at small pickup angles to stream away from the ribbon directions.

Whence the Interstellar Magnetic Field Shaping the Heliosphere?

Abstract Measurements of starlight polarized by aligned interstellar dust grains are used to probe the relation between the orientation of the ambient interstellar magnetic field (ISMF) and the ISMF traced by the ribbons of energetic neutral atoms discovered by the Interstellar Boundary Explorer spacecraft. We utilize polarization data, many acquired specifically for this study, to trace the configuration of the ISMF within 40 pc. A statistical analysis yields a best-fit ISMF orientation, B magpol, aligned with Galactic coordinates ℓ = 42°, b = 49°. Further analysis shows the ISMF is more orderly for “downfield” stars located over 90° from B magpol. The data subset of downfield stars yields an orientation for the nearby ISMF at ecliptic coordinates λ, β ≈ 219° ± 15°, 43° ± 9° (Galactic coordinates l, b ≈ 40°, 56°, ±17°). This best-fit ISMF orientation from polarization data is close to the field direction obtained from ribbon models. This agreement suggests that the ISMF shaping the heliosphere belongs to an extended ordered magnetic field. Extended filamentary structures are found throughout the sky. A previously discovered filament traversing the heliosphere nose region, “Filament A,” extends over 300° of the sky, and crosses the upwind direction of interstellar dust flowing into the heliosphere. Filament A overlaps the locations of the Voyager kilohertz emissions, three quasar intraday variables, cosmic microwave background (CMB) components, and the inflow direction of interstellar grains sampled by Ulysses and Galileo. These features are likely located in the upstream outer heliosheath where ISMF drapes over the heliosphere, suggesting Filament A coincides with a dusty magnetized plasma. A filament 55° long is aligned with a possible shock interface between local interstellar clouds. A dark spot in the CMB is seen within 5° of the filament and within 10° of the downfield ISMF direction. Two large magnetic arcs are centered on the directions of the heliotail. The overlap between CMB components and the aligned dust grains forming Filament A indicates the configuration of dust entrained in the ISMF interacting with the heliosphere provides a measurable foreground to the CMB.

Particle-in-cell simulations of high-frequency waves driven by pickup ion ring-beam distributions in the outer heliosheath

ABSTRACT Scattering of pickup ion ring-beam distributions in the outer heliosheath is a fundamental element in the spatial retention scenario of the energetic neutral atom (ENA) ribbon observed by the Interstellar Boundary EXplorer (IBEX). According to our earlier linear instability analysis, pickup ion ring-beam distributions trigger magnetic field-aligned, right-hand polarized unstable waves in two separate frequency ranges which are near and far above the proton cyclotron frequency, respectively. We have performed hybrid simulations to study the unstable waves near the proton cyclotron frequency. However, the high-frequency waves well above the proton cyclotron frequency are beyond the reach of hybrid simulations. In this paper, particle-in-cell simulations are carried out to study the parallel- and antiparallel-propagating high-frequency waves excited by the outer heliosheath pickup ions at different pickup angles as well as the scattering of the pickup ions by the waves excited. In the early stages of the simulations, the results confirm the excitation of the parallel-propagating, right-hand polarized high-frequency waves as predicted by the earlier linear analysis. Later in the simulations, enhanced antiparallel-propagating modes also emerge. Furthermore, the evolution of the pickup ion ring-beam distributions of the selected pickup angles reveals that the high-frequency waves do not significantly contribute to the pickup ion scattering. These results are favourable regarding the plausibility of the spatial retention scenario of the IBEX ENA ribbon.

Hybrid simulations of the ring-beam instabilities driven by the pickup ions in the outer heliosheath

ABSTRACT One-dimensional hybrid simulations are performed to investigate the plasma instabilities driven by the pickup ions in the outer heliosheath, which are believed to have a ring-beam velocity distribution in general. The modes propagating parallel and antiparallel to the background magnetic field are studied for pickup ion ring-beam distributions of different pickup angles, following the linear kinetic instability analysis presented in our previous work. The simulation results first confirm the unstable modes predicted by the linear analysis and further reveal how the pickup ions are scattered by the excited waves. While the maximum growth rate of the low-frequency, left-helicity waves occurs at the pickup angle of 82° (consistent with the linear analysis results), the saturation level of the fluctuating magnetic field is the highest around 45° pickup angle. For the pickup ion parameters examined, the pitch-angle spread of the pickup ions at the end of the simulations decreases with increasing pickup angle. Our study goes beyond previous studies of the pickup-ion-driven instabilities in the outer heliosheath which mainly focused on ring-like distributions of around 90° pickup angle. With the caveats that the pickup ion parameters are in the Alfvén cyclotron stability gap and that the one-dimensional simulations only include parallel and antiparallel-propagating modes, our results show that the scattering of pickup ions at small pickup angles is limited to one hemisphere so they can drift away from the ribbon direction as required by the spatial retention scenario of the energetic neutral atom ribbon observed by the Interstellar Boundary EXplorer.

Stability analysis of the pickup ion ring-beam distributions in the outer heliosheath

ABSTRACT The electromagnetic instabilities driven by the pickup ions in the outer heliosheath play a crucial role in interpreting the energetic neutral atom (ENA) ribbon observed by the Interstellar Boundary EXplorer (IBEX). Previous studies on the stability of the outer heliosheath pickup ions focused on ring-like velocity distributions with pickup angles equal or close to 90°. In this study, we investigate the instabilities driven by a ring-beam velocity distribution of the pickup ions in the outer heliosheath using linear instability analysis. The modes propagating along (both parallel and antiparallel to) the background magnetic field are examined, and the emphasis is on how the instabilities change with the pickup angle, α, varying from zero to 90°. The parallel thermal spread of the pickup ions is chosen to make the plasma parameters lie in the Alfvén cyclotron ‘stability gap’ to exclude the well-studied Alfvén cyclotron instability. Our linear instability analysis reveals that the unstable modes all have left helicity but may propagate parallel or antiparallel to the background magnetic field. At small α, the parallel-propagating modes are unstable in two separate wavenumber (and frequency) ranges. The two unstable ranges merge into one when α increases to 82°. In addition, the antiparallel-propagating modes become unstable when α is above 5°, but the parallel-propagating modes have larger growth rates until α becomes larger than αc = 82.8°. The growth rates of the parallel-propagating modes quickly drop below zero when α goes beyond αc, while the antiparallel-propagating modes remain unstable until α approaches 90°.

Variations in the Pickup Ion Density Structure in Response to the Growth of the Kelvin–Helmholtz Instability along the Heliopause

Abstract Features of the response of pickup ions (PUIs) to the Kelvin–Helmholtz instability (KHI) on the heliopause (HP) are examined by means of two-dimensional hybrid simulations. We assume the supersonic neutral solar wind as the source of PUIs gyrating about the magnetic field in the outer heliosheath. These PUIs become energetic neutral atoms (ENAs) via charge exchange with interstellar hydrogen, and a portion of these ENAs are detected by spacecraft such as the Interstellar Boundary Explorer (IBEX). To evaluate the possibility of identifying the KHI on the HP from ENA observations, we assume that an imprint of the KHI may be displayed in spatial and temporal variations in the observed ENA profile. As an alternative to ENA, the column density of PUIs integrated across the HP is calculated. The KH-inducing vortex forces not only background protons but also PUIs to roll up deep in the inner heliosheath. The KH vortex also results in the emission of magnetosonic pulses that sweep PUIs in the outer heliosheath and lead to their local confinement. These effects elongate the PUIs’ spatial distribution in the direction normal to the HP. The appearance of the local confined structure in the PUIs’ column density is consequently confirmed, and this feature can be confirmed as the KHI evolution. Although the simulation cannot be quantitatively compared with the observations currently available because its resolution is too low, we expect that the derived properties will be useful for diagnosing the nature of HP fluctuation in future missions.

Slowdown and Heating of Interstellar Neutral Helium by Elastic Collisions Beyond the Heliopause.

Direct sampling of interstellar neutral (ISN) atoms close to the Sun enables studies of the very local interstellar medium (VLISM) around the heliosphere. The primary population of ISN helium atoms has, until now, been assumed to reflect the pristine VLISM conditions at the heliopause. Consequently, the atoms observed at 1 au by the Interstellar Boundary Explorer (IBEX) were used to determine the VLISM temperature and velocity relative to the Sun, without accounting for elastic collisions with other species outside the heliopause. Here, we evaluate the effect of these collisions on the primary ISN helium population. We follow trajectories of helium atoms and track their collisions with slowed plasma and interstellar hydrogen atoms ahead of the heliopause. Atoms typically collide a few times in the outer heliosheath, and only ~1.5% of the atoms are not scattered at all. We use calculated differential cross sections to randomly choose scattering angles in these collisions. We estimate that the resulting primary ISN helium atoms at the heliopause are slowed down by ~0.45 km s-1 and heated by ~1100 K compared to the pristine VLISM. The resulting velocity distribution is asymmetric and shows an extended tail in the antisunward direction. Accounting for this change in the parameters derived from IBEX observations gives the Sun's relative speed of 25.85 km s-1 and temperature of 6400 K in the pristine VLISM. Finally, this paper serves as a source of the differential cross sections for elastic collisions with helium atoms.

Mirror Instability Driven by Pickup Ions in the Outer Heliosheath

Abstract The stability of the pickup ions in the outer heliosheath is a vital factor in the generation of the energetic neutral atom (ENA) ribbon observed by the Interstellar Boundary EXplorer according to the secondary ENA mechanism. Most previous studies of the pickup ion stability assumed simple, idealized velocity distributions of the pickup ions and focused on the parallel-propagating modes only. This paper takes a more realistic multicomponent pickup ion velocity distribution given by the global modeling of neutral atoms in the heliosphere and investigates the role of the oblique mirror waves, in addition to the parallel modes. Both linear kinetic instability analysis and hybrid simulations are performed. In contrast to a recent investigation using the same distribution that demonstrated the growth of parallel waves but reported an insignificant contribution of oblique modes, our study reveals substantial growth of the oblique mirror instability. The oblique mirror modes and the parallel/quasi-parallel ion cyclotron waves grow simultaneously with different growth rates. The pickup ion scattering by two types of excited waves together is more pronounced than by either type of wave alone. More importantly, our two-dimensional hybrid simulation results demonstrate that the development of the mirror instability not only produces its own pickup ion scattering, but also leads to the growth of extra ion cyclotron waves (in a quasi-linear manner), which further enhances the pickup ion scattering. The results suggest an important role of the mirror modes that should not be ignored in the stability study of the outer heliosheath pickup ions.

Probing the thermodynamic conditions of the heliosheath plasma by shock wave propagation

Context. The pressure equilibrium between the inner heliosheath and the outer heliosheath (referred to as the local interstellar medium) is an eminent theoretical and practical problem; theoretical, because the relevant pressure carriers have to be identified, and practical, because data must be gathered in order to confirm such a pressure equilibrium. The problem is closely connected with the stability of the heliopause, that is, of the tangential discontinuity between these two counterflowing media, and is of utmost importance for understanding the stability of the whole circumsolar plasma structure. Aims. In this paper we analyze the thermodynamic conditions of the multi-fluid plasma between the solar wind termination shock and the heliopause determining the total heliosheath pressure. We look into this problem from a theoretical standpoint and revisit theoretical descriptions of the solar wind plasma after its passage over the solar wind termination shock, thereafter forming the subsonic heliosheath region. Methods. Hereby we take into account the 3D magnetohydrodynamics shock conditions and the resulting 3D temperature structure of the downstream plasma flow. We use a kind of seismological procedure to probe the heliosheath plasma by inquiring into the propagation conditions of traveling shock wave perturbations in this predetermined 3D heliosheath plasma structure. We discuss the fact that the front geometry of such a traveling shock wave most probably does not remain spherical, if it was to begin with, due to asymmetric shock propagation conditions. In contrast, the wave front is likely to become strongly deformed into an upwind bulge. Results. Concerning the plasma pressure, in addition to solar wind and pick-up proton pressures, we have to take into account the solar wind electron pressure which as a surprise turns out to be of comparable magnitude. As a consequence, the characteristic propagation speed of the traveling shock wave in the weakly magnetized heliosheath plasma is given as a mixed speed expressed by the sound speeds of the protons and the electrons. We describe local low-energy proton density signatures that can be found in Voyager-2 proton data as a consequence of traveling shock wave passages and show that the total local plasma pressure can be directly derived from them.

Heliospheric Lyman-α Absorption Toward Voyager 2

Motivated by recent Voyager 2 (V2) measurements of surprisingly high plasma temperatures outside the heliopause, we study the heliospheric hydrogen Lyman-a absorption observed by the Hubble Space Telescope for the closest available line of sight to the V2 direction, 5 Pav, which is 9° from V2. The absorption is well reproduced by a typical global heliospheric model that predicts outer heliosheath temperatures of T«20, 000 K. Arbitrarily increasing these temperatures by a factor of two to be more consistent with the initial T=30, 000-50, 000 K values measured by V2 leads to significant overprediction of absorption. Thus, the high plasma temperatures first encountered by V2 near the heliopause cannot be present throughout the outer heliosheath.