Engineering topics

Heerikhuisen, J.

Publications and source records attributed to Heerikhuisen, J..

At least 19 records

DOE OSTI2021

Turbulent Acceleration of Interstellar Pickup Ions at the Heliospheric Termination Shock Forms the Global ENA Spectrum

The heliospheric energetic neutral atom spectrum observed by the Interstellar Boundary Explorer (IBEX) reveals that the heliosheath proton distribution is consistent with a power law. The origin of the spectrum is likely from interstellar pickup ions (PUIs) accelerated at the heliospheric termination shock (HTS). We present an explanation of the proton spectrum origin using a test particle simulation of PUIs accelerated at the HTS. PUIs experience preferential heating by the motional electric field in the shock foot, but do not develop a power-law tail without the presence of turbulence at wavenumbers (k) close to the PUI gyroradius scale (R g ). Voyager 2 observations of the magnetic field downstream of the HTS indicate a moderate amount of turbulence at kR g ≅ 1, (δB /B 0 ) 2 @ 0.01, which we find to be sufficient for producing a downstream suprathermal PUI tail but not at intensities observed by IBEX. Within the shock ramp, however, Voyager observed the turbulence power at much smaller scales to be nearly 100 times stronger, suggesting the possibility of strong turbulence at the PUI gyroradius scale. We show that a proton distribution can develop a power law downstream of the HTS consistent with IBEX observations if (δB /B 0 ) 2 ≳ 0.1 at kR g ≅ 1 in the shock foot. Shock drift acceleration of PUIs by the motional electric field is aided by interactions with turbulence upstream of the shock overshoot. Steepening of the IBEX proton spectrum in directions farther from the heliospheric nose suggests the HTS compression ratio and/or turbulence power weakens near the heliotail.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Variability in the Position of the IBEX Ribbon over Nine Years: More Observational Evidence for a Secondary ENA Source

The ribbon of enhanced energetic neutral atom flux, discovered by the Interstellar Boundary Explorer (IBEX) in 2009, has redefined our understanding of the heliosphere's interaction with the local interstellar medium (LISM). Yet, its origin continues to be a topic of scientific debate. The ribbon is circular and traces the region where the putative LISM magnetic field (B LISM ) is perpendicular to the radial direction from the Sun. Using nine years of IBEX-Hi observations, we investigate the ribbon circularity and location as functions of time and energy. We provide updated locations of the ribbon center at five energy passbands (centered at 0.7, 1.1, 1.7, 2.7, and 4.3 keV) in ecliptic coordinates [longitude, latitude]: [217$_˙^°$dg41 ± 0$_˙^°$dg95, 44$_˙^°$dg36 ± 0$_˙^°$dg93], [219$_˙^°$dg72 ± 0$_˙^°$dg95, 41$_˙^°$dg50 ± 0$_˙^°$dg87], [220$_˙^°$dg51 ± 1$_˙^°$dg19, 39$_˙^°$dg96 ± 1$_˙^°$dg00], [218$_˙^°$dg08 ± 1$_˙^°$dg66, 38$_˙^°$dg44 ± 1$_˙^°$dg24], and [214$_˙^°$dg68 ± 1$_˙^°$dg48, 34$_˙^°$dg13 ± 1$_˙^°$dg19] respectively. The weighted mean center location over all energies and all years is [218$_˙^°$dg33 ± 0$_˙^°$dg68, 40$_˙^°$dg38 ± 0$_˙^°$dg88] and its radius is 74$_˙^°$dg81 ± 0$_˙^°$dg65. As viewed by IBEX at 1 au, we find that (1) the ribbon is stable over time, with distinct centers at each energy; (2) ribbon centers exhibit small temporal variations, likely caused by the solar wind (SW) speed and density variations; and (3) ribbon location in the sky appears to be driven by (i) the inherent alignment of the ribbon centers along the plane connecting the presumed B LISM and the heliospheric upwind direction, and (ii) the variable SW structure along the heliographic meridian, further emphasizing that the ribbon source is outside the heliosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Constraining the Evolution of the Proton Distribution Function in the Heliotail

We use Interstellar Boundary Explorer (IBEX) measurements of energetic neutral atoms (ENAs) to constrain the proton (mostly pickup ion, PUI) distribution in the heliotail. In our previous study, we solved the Parker transport equation and found that the velocity diffusion coefficient $D(v)$ for PUIs is approximately $D(v)$ ~ 1.1 × 10 -8 km 2 s -3 ($v/v_0$) 1.3 , assuming the initial proton distribution processed by the termination shock (TS), $f_{p,0}$, is a kappa distribution with kappa index $ₖ_{p,0}$ = 1.63. In this study, we test different forms for $f_{p,0}$. We find that if $f_{p,0}$ is kappa-distributed and $D(v)$ = $D_0(v/v_0)^{1.3}$, any kappa index in the range 1.5 < $κ_{p,0}$ < 10 is consistent with IBEX data if $D_0 ~$ 0.8–1.3 × 10 -8 km 2 s -3 . While the case where $D(v)$ ∝ $v^{1.3}$ yields ENA fluxes that appear to best reproduce IBEX data for any $κ_{p,0}$, it is possible for $D(v)$ to scale close to $~v^{2/3}$ or $~v^2$ within our uncertainties by changing $D_0$. We also show that an upstream PUI filled-shell distribution that is heated by a quasi-stationary TS, generating a downstream filled-shell with large cutoff speed, yields an excess of ENAs > 2 keV compared to IBEX. However, using a fully kinetic particle-in-cell simulation to process a PUI filled-shell across the TS yields ENA spectra consistent with IBEX , reinforcing the significance of self-consistent, preferential PUI heating and diffusion at the TS. Interestingly, an upstream PUI distribution inferred from the particle-in-cell simulation to reproduce Voyager 2 observations of the nose-ward TS is inconsistent with IBEX observations from the heliotail, suggesting differences in the upstream PUI distribution or TS properties.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Stochastic Acceleration of ~0.1–5 keV Pickup Ions in the Heliotail

We seek to understand the quantitative role of the dominant physical processes (charge-exchange, adiabatic heating, stochastic acceleration) governing the proton distribution in the heliotail using observations of hydrogen energetic neutral atoms (ENAs) from the Interstellar Boundary Explorer ( IBEX ). We solve the Parker transport equation for solar wind protons and pickup ions (PUIs) as they propagate from the termination shock (TS) down the heliotail, including charge-exchange between protons and neutral hydrogen atoms as source terms derived from an MHD-fluid and kinetic-neutral simulation of the heliosphere. We compute ENA fluxes at 1 au from the results of the proton transport model and compare them with IBEX observations. We find that, under the assumptions of our model, a stochastic acceleration process is needed to counteract the energy-dependent losses of ~0.1–5 keV PUIs from charge-exchange to reproduce IBEX data. The power-law velocity dependence of the diffusion coefficient (spectral index γ ) is limited to the range 0.67 < γ < 2, and the best fit to IBEX data appears close to γ ~ 1.25. The diffusion rate ~1.1 × 10 -8 km 2 s -3 ( v / v 0 ) 1.25 nearly balances the loss of ~0.1–5 keV PUIs by charge-exchange. Our analysis suggests that cyclotron resonance with two widely known incompressible MHD turbulence: namely, isotropic Kolmogorov and anisotropic Goldreich–Sridhar turbulence, as well as stochastic particle interactions with compressive waves are not by themselves the dominant diffusion mechanisms. However, some intermediate processes may be occurring due to the presence of PUIs.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2017

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

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 .

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Structure of the Heliotail from Interstellar Boundary Explorer Observations: Implications for the 11-year Solar Cycle and Pickup Ions in the Heliosheath

Interstellar Boundary Explorer ( IBEX ) measurements of energetic neutral atoms (ENAs) from the heliotail show a multi-lobe structure of ENA fluxes as a function of energy between ∼0.71 and 4.29 keV. Below ∼2 keV, there is a single structure of enhanced ENA fluxes centered near the downwind direction. Above ∼2 keV, this structure separates into two lobes, one north and one south of the solar equatorial plane. ENA flux from these two lobes can be interpreted as originating from the fast solar wind (SW) propagating through the inner heliosheath (IHS). Alternatively, a recently published model of the heliosphere suggests that the heliotail may split into a “croissant-like” shape, and that such a geometry could be responsible for the heliotail ENA feature. Here we present results from a time-dependent simulation of the heliosphere that produces a comet-like heliotail, and show that the 11-year solar cycle leads to the formation of ENA lobes with properties remarkably similar to those observed by IBEX . The ENA energy at which the north and south lobes appear suggests that the pickup ion (PUI) temperature in the slow SW of the IHS is ∼10{sup 7} K. Moreover, we demonstrate that the extinction of PUIs by charge-exchange is an essential process required to create the observed global ENA structure. While the shape and locations of the ENA lobes as a function of energy are well reproduced by PUIs that cross the termination shock, the results appear to be sensitive to the form of the distribution of PUIs injected in the IHS.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

LATITUDE, ENERGY, AND TIME VARIATIONS IN THE ENERGETIC NEUTRAL ATOM SPECTRAL INDICES MEASURED BY THE INTERSTELLAR BOUNDARY EXPLORER ( IBEX )

We investigate the latitude, energy, and temporal variations of the ∼0.5–6 keV energetic neutral atom (ENA) fluxes and spectral indices measured by the Interstellar Boundary Explorer ( IBEX ) outside of the Ribbon from 2009 to 2013. By combining all-sky maps for years 1–3 and years 4–5, we find that the >1.2 keV globally distributed (GD) ENA fluxes at all latitudes decrease by ∼25% from years 1–3 to years 4–5, but there is no change (<5%) in the corresponding spectral indices. We also show that the latitudinal profile and energy-dependence of the ≥1 keV ENA spectral indices outside the Ribbon exhibit no significant time differences between years 1–3 and 4–5, lending strong support for the notion that the highly organized, persistent energy dependence and latitudinal pattern of the ENA spectral indices during 2009–2014 are determined by that of the solar wind (SW) speed observed in the inner heliosphere during the deep solar minimum conditions of 2006–2010. In contrast, the ∼25% decrease in the >1.2 keV ENA fluxes at all latitudes from years 1–3 to years 4–5 occurs because the magnitude of the driving SW parameters (SW density or dynamic pressure) diminished significantly from 2006 to 2010. Based on the reconstructed latitudinal and temporal profiles of SW parameters from 2011 to 2014, i.e., during the rising phase of solar cycle 24, we suggest that the GD ENA fluxes in years 2014–2017, i.e., in maps 6–9, will either stabilize or increase and the latitudinal pattern and energy dependence of the corresponding spectral indices will be disrupted.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

THE IBEX RIBBON AND THE PICKUP ION RING STABILITY IN THE OUTER HELIOSHEATH. I. THEORY AND HYBRID SIMULATIONS

The nearly circular band of energetic neutral atom emission dominating the field of view of the Interplanetary Boundary Explorer ( IBEX ) satellite, is most commonly attributed to the effect of charge exchange of secondary pickup ions (PUIs) gyrating about the magnetic field in the outer heliosheath and the interstellar space beyond. Several models for the PUI dynamics of this mechanism have been proposed, each requiring either strong or weak scattering of the initial pitch angle. Conventional wisdom states that ring distributions tend to generate waves and scatter onto a shell on timescales too short for charge exchange to occur. We performed a careful study of ring and thin shell proton distribution stability using theoretical tools and hybrid plasma simulations. We show that the kinetic behavior of a freshly injected proton ring is a far more complicated process than previously thought. In the presence of a warm Maxwellian core, narrower rings could be more stable than broader toroidal distributions. The scattered rings possess a fine structure that can only be revealed using very large numbers of macroparticles in a simulation. It is demonstrated that a “stability gap” in ring temperature exists where the protons could retain large gyrating anisotropies for years, and the wave activity could remain below the level of the ambient magnetic fluctuations in interstellar space. In the directions away from the ribbon, however, a partial shell distribution is more likely to be unstable, leading to significant scattering into one hemisphere in velocity space. The process is accompanied by turbulence production, which is puzzling given the very low level of magnetic fluctuations measured in the outer heliosheath by Voyager 1 .

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

THE IBEX RIBBON AND THE PICKUP ION RING STABILITY IN THE OUTER HELIOSHEATH. II. MONTE-CARLO AND PARTICLE-IN-CELL MODEL RESULTS

The nearly circular ribbon of energetic neutral atom (ENA) emission discovered by NASA’s Interplanetary Boundary EXplorer satellite ( IBEX ), is most commonly attributed to the effect of charge exchange of secondary pickup ions (PUIs) gyrating about the magnetic field in the outer heliosheath (OHS) and the interstellar space beyond. The first paper in the series (Paper I) presented a theoretical analysis of the pickup process in the OHS and hybrid-kinetic simulations, revealing that the kinetic properties of freshly injected proton rings depend sensitively on the details of their velocity distribution. It was demonstrated that only rings that are not too narrow (parallel thermal spread above a few km s{sup −1}) and not too wide (parallel temperature smaller than the core plasma temperature) could remain stable for a period of time long enough to generate ribbon ENAs. This paper investigates the role of electron dynamics and the extra spatial degree of freedom in the ring ion scattering process with the help of two-dimensional full particle-in-cell (PIC) kinetic simulations. A good agreement is observed between ring evolution under unstable conditions in hybrid and PIC models, and the dominant modes are found to propagate parallel to the magnetic field. We also present more realistic ribbon PUI distributions generated using Monte Carlo simulations of atomic hydrogen in the global heliosphere and examine the effect of both the cold ring-like and the hot “halo” PUIs produced from heliosheath ENAs on the ring stability. It is shown that the second PUI population enhances the fluctuation growth rate, leading to faster isotropization of the solar-wind-derived ring ions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

GEOMETRY AND CHARACTERISTICS OF THE HELIOSHEATH REVEALED IN THE FIRST FIVE YEARS OF INTERSTELLAR BOUNDARY EXPLORER OBSERVATIONS

We investigate and interpret the geometry and characteristics of the inner heliosheath (IHS) plasma and their impact on the heliotail structure as observed in energetic neutral atom (ENA) maps acquired during the first 5 yr of Interstellar Boundary Explorer ( IBEX ) observations. In particular, IBEX observations of the heliotail reveal distinct, localized emission features (lobes) that provide a rich set of information about the properties and evolution of the heliosheath plasma downstream of the termination shock (TS). We analyze the geometry of the heliotail lobes and find that the plane intersecting the port and starboard heliotail lobe centers is ∼6° from the solar equatorial plane, and the plane intersecting the north and south heliotail lobe centers is ∼90° from the solar equatorial plane, both indicating strong correlation with the fast–slow solar wind asymmetry, and thus reflecting the structure of the IHS flow around the Sun. We also analyze the key parameters and processes that form and shape the port and starboard lobes, which are distinctly different from the north and south lobes. By comparing IBEX ENA observations with results from a simplistic flow model of the heliosphere and a multicomponent description for pickup ions (PUIs) in the IHS, we find that the port and starboard lobe formation is driven by a thin IHS, large nose–tail asymmetry of the distance to the TS (and consequently, a large nose–tail asymmetry of the relative abundance of PUIs at the TS) and the energy-dependent removal of PUIs by charge exchange in the IHS.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

LOCAL INTERSTELLAR MAGNETIC FIELD DETERMINED FROM THE INTERSTELLAR BOUNDARY EXPLORER RIBBON

The solar wind emanating from the Sun interacts with the local interstellar medium (LISM), forming the heliosphere. Hydrogen energetic neutral atoms (ENAs) produced by the solar-interstellar interaction carry important information about plasma properties from the boundaries of the heliosphere, and are currently being measured by NASA's Interstellar Boundary Explorer (IBEX). IBEX observations show the existence of a “ribbon” of intense ENA emission projecting a circle on the celestial sphere that is centered near the local interstellar magnetic field (ISMF) vector. Here we show that the source of the IBEX ribbon as a function of ENA energy outside the heliosphere, uniquely coupled to the draping of the ISMF around the heliopause, can be used to precisely determine the magnitude (2.93 ± 0.08 μG) and direction (227.°28 ± 0.°69, 34.°62 ± 0.°45 in ecliptic longitude and latitude) of the pristine ISMF far (∼1000 AU) from the Sun. We find that the ISMF vector is offset from the ribbon center by ∼8.°3 toward the direction of motion of the heliosphere through the LISM, and their vectors form a plane that is consistent with the direction of deflected interstellar neutral hydrogen, thought to be controlled by the ISMF. Our results yield draped ISMF properties close to that observed by Voyager 1, the only spacecraft to directly measure the ISMF close to the heliosphere, and give predictions of the pristine ISMF that Voyager 1 has yet to sample.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

THE HELIOTAIL

The heliotail is formed when the solar wind (SW) interacts with the local interstellar medium (LISM) and is shaped by the interstellar magnetic field (ISMF). While there are no spacecraft available to perform in situ measurements of the SW plasma and heliospheric magnetic field (HMF) in the heliotail, it is of importance for the interpretation of measurements of energetic neutral atom fluxes performed by Interstellar Boundary Explorer. It has been shown recently that the orientation of the heliotail in space and distortions of the unperturbed LISM caused by its presence may explain the anisotropy in the TeV cosmic ray flux detected in air shower observations. The SW flow in the heliotail is a mystery itself because it is strongly affected by charge exchange between the SW ions and interstellar neutral atoms. If the angle between the Sun’s magnetic and rotation axes is constant, the SW in the tail tends to be concentrated inside the HMF spirals deflected tailward. However, the twisted field soon becomes unstable and the reason for the SW collimation within a two-lobe structure vanishes. We demonstrate that kinetic treatment of the H atom transport becomes essential in this case for explaining the lobe absence further along the tail. We show that the heliotail flow is strongly affected by the solar cycle that eliminates artifacts, which is typical of solutions based on simplifying assumptions. The heliopause in the tail is subject to Kelvin–Helmholtz instability, while its orientation and shape are determined by the ISMF direction and strength.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

SIMULATIONS OF A DYNAMIC SOLAR CYCLE AND ITS EFFECTS ON THE INTERSTELLAR BOUNDARY EXPLORER RIBBON AND GLOBALLY DISTRIBUTED ENERGETIC NEUTRAL ATOM FLUX

Since 2009, observations by the Interstellar Boundary Explorer (IBEX) have vastly improved our understanding of the interaction between the solar wind (SW) and local interstellar medium through direct measurements of energetic neutral atoms (ENAs), which inform us about the heliospheric conditions that produced them. An enhanced feature of flux in the sky, the so-called IBEX ribbon, was not predicted by any global models before the first IBEX observations. A dominating theory of the origin of the ribbon, although still under debate, is a secondary charge-exchange process involving secondary ENAs originating from outside the heliopause. According to this mechanism, the evolution of the solar cycle should be visible in the ribbon flux. Therefore, in this paper we simulate a fully time-dependent ribbon flux, as well as globally distributed flux from the inner heliosheath (IHS), using time-dependent SW parameters from Sokół et al. as boundary conditions for our time-dependent heliosphere simulation. After post-processing the results to compute H ENA fluxes, our results show that the secondary ENA ribbon indeed should be time dependent, evolving with a period of approximately 11 yr, with differences depending on the energy and direction. Our results for the IHS flux show little periodic change with the 11 yr solar cycle, but rather with short-term fluctuations in the background plasma. While the secondary ENA mechanism appears to emulate several key characteristics of the observed IBEX ribbon, it appears that our simulation does not yet include all of the relevant physics that produces the observed ribbon.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Latitudinal and energy dependence of energetic neutral atom spectral indices measured by the Interstellar Boundary Explorer

We investigate the latitudinal and energy dependence of the globally distributed 0.5–6 keV energetic neutral atom (ENA) spectra measured by the Interstellar Boundary Explorer (IBEX) during the first 3 yrs of the mission. Our results are: (1) the ENA spectral indices at the two lowest energies (0.89 and 1.47 keV) exhibit no clear trend with ecliptic latitude θ, while those at ∼2.29 and ∼3.41 keV exhibit a clear latitudinal pattern; flatter spectra occur above 60° latitude and steeper spectra occur ±30° of the equator. (2) The latitudinal dependence of the spectral indices at different energies can be represented by the cosine function γ=a{sub 0}+a{sub 1}cos(a{sub 2}θ) with unique offsets, amplitudes, and phase angles; the higher energy ENA indices transition to successively larger amplitudes within ±45° of the equator. Our results confirm the previously reported latitudinal organization of the ENA spectra and their remarkable similarity to that of the solar wind (SW) speed observed by Ulysses in the inner heliosphere. While earlier studies showed that the ∼0.5–6 keV globally distributed ENA spectral indices could be represented as single power laws over much of the sky, our new results indicate that this is an over-simplification because the spectral indices have an energy and latitude dependence. This dependence is an important factor that must be taken into consideration by models and simulations that seek to map the IBEX ENA observations back to the latitudinal profile of the SW speed structure observed in the inner heliosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

κ -distributed protons in the solar wind and their charge-exchange coupling to energetic hydrogen

The interaction between the solar wind and the interstellar medium represents a collision between two plasma flows, resulting in a heliosphere with an extended tail. While the solar wind is mostly ionized material from the corona, the interstellar medium is only partially ionized. The ion and neutral populations are coupled through charge-exchange collisions that operate on length scales of tens to hundreds of astronomical units. About half the interstellar hydrogen flows into the heliosphere where it may charge-exchange with solar wind protons. This process gives rise to a nonthermal proton, known as a pickup ion, which joins the plasma. In this paper we investigate the effects of approximating the total ion distribution of the subsonic solar wind as a generalized Lorentzian, or κ distribution, using an MHD neutral code. We illustrate the effect different values of the κ parameter have on both the structure of the heliosphere and the energetic neutral atom flux at 1 AU. We find that using a κ distribution in our simulations yields levels of energetic neutral atom flux that are within a factor of about 2 or 3 over the IBEX-Hi range of energies from 0.5 to 6 keV. In conclusion, while the presence of a suprathermal tail in the proton distribution leads to the production of high-energy neutrals, the sharp decline in the charge-exchange cross section around 10 keV mitigates the enhanced transfer of energy from the ions to the neutrals that might otherwise be expected.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Turbulent Heating of the Core Solar Wind Protons in the Outer Heliosphere

A model of the interstellar pickup protons, slab component of the Alfvénic turbulence, and core solar wind (SW) protons is presented for r ≥ 1 AU along with a comparison with the Voyager 2 (V2) observations of the core SW temperature. The main conclusions of our study are summarized as follows. (1) A combined effect of the Wentzel-Kramers-Brillouin attenuation, large-scale driving, and pickup proton generated waves results in the turbulent energy sink in the region 1 AU < r < 10 AU, while energy is pumped in the turbulence beyond ˜ 10 AU. Without energy pumping the nonlinear energy cascade is suppressed for r < 10 AU. A similar situation takes place for the 2D turbulence. (2) For both the slab and 2D turbulent components at the heliocentric distances r < 10 AU, the energy source for the core SW protons due to the Alfvén wave dissipation is small. As a result, adiabatic cooling mostly controls the model SW temperature in this region, and the model temperature disagrees with the V2 observations in the region r < 20 AU. (3) The model core SW temperature agrees well with the V2 observations for r > 20 AU.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

The effect of new interstellar medium parameters on the heliosphere and energetic neutral atoms from the interstellar boundary

We present new results from three-dimensional simulations of the solar wind interaction with the local interstellar medium (LISM) using recent observations by NASA's Interstellar Boundary EXplorer (IBEX) mission estimates of the velocity and temperature of the LISM. We investigate four strengths of the LISM magnetic field, from 1 to 4 μG, and adjust the LISM proton and hydrogen densities so that the distance to the termination shock (TS) in the directions of the Voyager spacecraft is just below 90 AU, and the density of hydrogen at the TS is close to 0.09 cm{sup –3} in the nose direction. The orientation of the magnetic field is chosen to point toward the center of the ribbon of enhanced energetic neutral atom (ENA) flux seen in the IBEX data. Our simulations show that the plasma and neutral properties in the outer heliosheath vary considerably as a function of the LISM magnetic field strength. We also show that the heliotail points downwind in all cases, though its structure is strongly affected by the external magnetic field. Comparison and consistency between the simulated ENA flux and the circularity of the ribbon as measured by IBEX are most consistent with a LISM magnetic field strength aligned with the center of the ribbon and a magnitude in the range 2.5-3 μG.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Charge-exchange coupling between pickup ions across the heliopause and its effect on energetic neutral hydrogen flux

Pickup ions (PUIs) appear to play an integral role in the multi-component nature of the plasma in the interaction between the solar wind (SW) and local interstellar medium (LISM). Three-dimensional (3D) MHD simulations with a kinetic treatment for neutrals and PUIs are currently still not viable. In light of recent energetic neutral atom (ENA) observations by the Interstellar Boundary EXplorer, the purpose of this paper is to illustrate the complex coupling between PUIs across the heliopause (HP) as facilitated by ENAs using estimates of PUI properties extracted from a 3D MHD simulation of the SW-LISM interaction with kinetic neutrals. First, we improve upon the multi-component treatment of the inner heliosheath (IHS) plasma from Zank et al. by including the extinction of PUIs through charge-exchange. We find a significant amount of energy is transferred away from hot, termination shock-processed PUIs into a colder, 'freshly injected' PUI population. Second, we extend the multi-component approach to estimate ENA flux from the outer heliosheath (OHS), formed from charge-exchange between interstellar hydrogen atoms and energetic PUIs. These PUIs are formed from ENAs in the IHS that crossed the HP and experienced charge-exchange. Our estimates, based on plasma-neutral simulations of the SW-LISM interaction and a post-processing analysis of ENAs and PUIs, suggest the majority of flux visible at 1 AU from the front of the heliosphere, between ∼0.02 and 10 keV, originates from OHS PUIs, indicating strong coupling between the IHS and OHS plasmas through charge-exchange.

79 ASTRONOMY AND ASTROPHYSICS↗