Local Interstellar Medium Plasma Research Articles

PIC Simulation of a Shock Tube: Implications for Wave Transmission in the Heliospheric Boundary Region

Abstract A shock tube problem is solved numerically by using one-dimensional full particle-in-cell simulations under the condition that a relatively tenuous and weakly magnetized plasma is continuously pushed by a relatively dense and strongly magnetized plasma having supersonic relative velocity. A forward and a reverse shock and a contact discontinuity are self-consistently reproduced. The spatial width of the contact discontinuity increases as the angle between the discontinuity normal and ambient magnetic field decreases. The inner structure of the discontinuity shows different profiles between magnetic field and plasma density, or pressure, which is caused by a non-MHD effect of the local plasma. The region between the two shocks is turbulent. The fluctuations in the relatively dense plasma are compressible and propagating away from the contact discontinuity, although the fluctuations in the relatively tenuous plasma contain both compressible and incompressible components. The source of the compressible fluctuations in the relatively dense plasma is in the relatively tenuous plasma. Only compressible fast mode fluctuations generated in the relatively tenuous plasma are transmitted through the contact discontinuity and propagate in the relatively dense plasma. These fast mode fluctuations are steepened when passing the contact discontinuity. This wave steepening and probably other effects may cause the broadening of the wave spectrum in the very local interstellar medium plasma. The results are discussed in the context of the heliospheric boundary region or heliopause.

Three-dimensional Features of the Outer Heliosphere Due to Coupling between the Interstellar and Heliospheric Magnetic Field. V. The Bow Wave, Heliospheric Boundary Layer, Instabilities, and Magnetic Reconnection

Abstract The heliosphere is formed due to interaction between the solar wind (SW) and local interstellar medium (LISM). The shape and position of the heliospheric boundary, the heliopause, in space depend on the parameters of interacting plasma flows. The interplay between the asymmetrizing effect of the interstellar magnetic field and charge exchange between ions and neutral atoms plays an important role in the SW–LISM interaction. By performing three-dimensional, MHD plasma/kinetic neutral atom simulations, we determine the width of the outer heliosheath—the LISM plasma region affected by the presence of the heliosphere—and analyze quantitatively the distributions in front of the heliopause. It is shown that charge exchange modifies the LISM plasma to such extent that the contribution of a shock transition to the total variation of plasma parameters becomes small even if the LISM velocity exceeds the fast magnetosonic speed in the unperturbed medium. By performing adaptive mesh refinement simulations, we show that a distinct boundary layer of decreased plasma density and enhanced magnetic field should be observed on the interstellar side of the heliopause. We show that this behavior is in agreement with the plasma oscillations of increasing frequency observed by the plasma wave instrument onboard Voyager 1. We also demonstrate that Voyager observations in the inner heliosheath between the heliospheric termination shock and the heliopause are consistent with dissipation of the heliospheric magnetic field. The choice of LISM parameters in this analysis is based on the simulations that fit observations of energetic neutral atoms performed by Interstellar Boundary Explorer.

ADVECTIVE TRANSPORT OF INTERSTELLAR PLASMA INTO THE HELIOSPHERE ACROSS THE RECONNECTING HELIOPAUSE

We discuss results of magnetohydrodynamical model simulations of plasma dynamics in the proximity of the heliopause (HP). The model is shown to fit details of the magnetic field variations observed by Voyager 1 spacecraft during the transition from the heliosphere to the local interstellar medium (LISM). We propose an interpretation of magnetic field structures observed by Voyager 1 in terms of fine-scale physical processes. Our simulations reveal an effective transport mechanism of relatively dense LISM plasma across the reconnecting HP into the heliosphere. The mechanism is associated with annihilation of magnetic sectors in the heliospheric plasma near the HP.

HELIOSPHERIC STRUCTURE: THE BOW WAVE AND THE HYDROGEN WALL

Recent IBEX observations indicate that the local interstellar medium (LISM) flow speed is less than previously thought (23.2 km s–1 rather than 26 km s–1). Reasonable LISM plasma parameters indicate that the LISM flow may be either marginally super-fast magnetosonic or sub-fast magnetosonic. This raises two challenging questions: (1) Can a LISM model that is barely super-fast or sub-fast magnetosonic account for Lyα observations that rely critically on the additional absorption provided by the hydrogen wall (H-wall)? and (2) If the LISM flow is weakly super-fast magnetosonic, does the transition assume the form of a traditional shock or does neutral hydrogen (H) mediate shock dissipation and hence structure through charge exchange? Both questions are addressed using three three-dimensional self-consistently coupled magnetohydrodynamic plasma—kinetic H models with different LISM magnetic field strengths (2, 3, and 4 μG) as well as plasma and neutral H number densities. The 2 and 3 μG models are fast magnetosonic far upwind of the heliopause whereas the 4 μG model is fully subsonic. The 2 μG model admits a broad (~50-75 AU) bow-shock-like structure. The 3 μG model has a smooth super-fast-sub-fast magnetosonic transition that resembles a very broad, ~200 AU thick, bow wave. A theoretical analysis shows that the transition from a super-fast to a sub-fast magnetosonic downstream state is due to the charge exchange of fast neutral H and hot neutral H created in the supersonic solar wind and hot inner heliosheath, respectively. For both the 2 μG and the 3 μG models, the super-fast magnetosonic LISM flow passes through a critical point located where the fast magnetosonic Mach number M = 1 and Qe = γ/(γ – 1)UQm , where Qe and Qm are the plasma energy and momentum source terms due to charge exchange, U is the LISM flow speed, and γ is the plasma adiabatic index. Because the Mach number is only barely super-fast magnetosonic in the 3 μG case, the hot and fast neutral H can completely mediate the transition and impose a charge exchange length scale on the structure, making the solar-wind-LISM interaction effectively bow-shock-free. The charge exchange of fast and hot heliospheric neutral H therefore provides a primary dissipation mechanism at the weak heliospheric bow shock, in some cases effectively creating a one-shock heliosphere (i.e., a heliospheric termination shock only). Both super-fast magnetosonic models produce a sizeable H-wall. We find that (1) a sub-fast magnetosonic LISM flow cannot model the observed Lyα absorption profiles along the four sightlines considered (α Cen, 36 Oph, DK UMa, and χ1 Ori—upwind, sidewind, and downwind respectively); (2) both the super-fast magnetosonic models can account for the Lyα observations, with possibly the bow-shock-free 3 μG model being slightly favored. Subject to further modeling and comparison against further lines of sight, we conclude with the tantalizing possibility that IBEX may have discovered a class of interstellar shocks mediated by neutral H.

Heliopause imaging in EUV: Oxygen O+ ion 83.4‐nm resonance line emission

We explore the possibility of remote, from 1 AU, study of the heliopause by an observer outside the geocorona. We argue that the heliopause, a boundary that separates the solar wind and the galactic plasma of the local interstellar medium (LISM), can be imaged by detecting solar extreme ultraviolet (EUV) radiation reflected by interstellar ions. Such EUV imaging would map the heliopause and provide important insight into its three‐dimensional structure and the LISM parameters as well. We consider heliopause mapping in the oxygen O+ ion resonance line (83.4 nm); imaging in the helium He+ ion line (30.4 nm) will be considered in a future article. We show that the expected heliopause brightness map at 83.4 nm is essentially different from that of the foreground glow of the solar wind O+ pickup ions. The interstellar plasma glow is brighter in the upwind (with respect to the interstellar wind) direction, while the pickup ion glow dominates in the downwind direction. The spectral characteristics of the radiation scattered by the LISM plasma and by the pickup ions are significantly different. The all‐sky images at 83.4 nm are highly sensitive to the ionization state of the LISM and would allow one to probe the asymmetry of the interstellar magnetic field. We briefly discuss the experimental requirements to heliopause EUV mapping, which would require 3 orders of magnitude improvement in instrumentation sensitivity. This is a challenging but not impossible task.

Effect of local interstellar medium hydrogen fractional ionization on the distant solar wind and interface region

The paper investigates the influence of the local interstellar medium (LISM) proton number density on the gasdynamics of LISM / solar wind interaction. It uses the self‐consistent axisymmetric model of Baranov and Malama (1993) for a constant value of the LISM neutral hydrogen number density (nH∞ = 0.14 cm−3). It is shown that the flow structure is asymmetric relative to the upwind and downwind directions. Numerical results are presented for four values of the proton number density in the LISM: np∞ = 0.005; 0.025; 0.1 and 0.28 cm−3. The results are given along the axis of symmetry only (in upwind and downwind directions). Velocity, number density, and temperature distributions of the plasma component as well as of the H atoms of the LISM are presented. It is shown that the influence of H atoms, penetrating from the LISM into the solar system on the solar wind may be strong enough to be observed by the “Voyager” spacecraft. The results show that the number density of the LISM's plasma component as well as of the solar wind has very large gradients in the interface region (between the bow and termination shocks). This effect may be very important for the interpretation of the kHz radiation detected by “Voyager”. Moreover, the asymmetric (in upwind and downwind directions) distributions of H atom number density, bulk velocity, and temperature may be observed on the basis of the scattered solar Lyman alpha radiation data.