Search

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 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 small, we expect that the derived properties will be useful for diagnosing the nature of HP fluctuation in future missions., Comment: 13 pages, 7 figures. Accepted for publication in ApJ

The study of multiscale pickup ion phase-mixing in the lunar plasma wake with a hybrid model is the main subject of our investigation in this paper. Photoionization and charge exchange of protons with the lunar exosphere are the ionization processes included in our model. The computational model includes the self-consistent dynamics of the light ( H + or H 2 + and He + ), and heavy ( Na + ) pickup ions. The electrons are considered as a fluid. The lunar interior is considered as a weakly conducting body. In this paper we considered for the first time the cumulative effect of heavy neutrals in the lunar exosphere (e.g., Al, Ar), an effect which was simulated with one species of Na + but with a tenfold increase in total production rates. We find that various species produce various types of plasma tail in the lunar plasma wake. Specifically, Na + and He + pickup ions form a cycloid-like tail, whereas the H + or H 2 + pickup ions form a tail with a high density core and saw-like periodic structures in the flank region. The length of these structures varies from 1.5 R M to 3.3 R M depending on the value of gyroradius for H + or H 2 + pickup ions. The light pickup ions produce more symmetrical jump in the density and magnetic field at the Mach cone which is mainly controlled by the conductivity of the interior, an effect previously unappreciated. Although other pickup ion species had little effect on the nature of the interaction of the Moon with the solar wind, the global structure of the lunar tail in these simulations appeared quite different when the H 2 + production rate was high.

Pickup ion detection at Titan is challenging because ion cyclotron waves are rarely detected in the vicinity of the moon. In this work, signatures left by freshly produced pick up heavy ions (m/q ∼ 16 to m/q ∼ 28) as detected in the plasma data by the CAPS/IMS instrument on-board Cassini are analyzed. In order to discern whether these correspond to ions of exospheric origin, one of the flybys during which the reported signatures were observed is investigated in detail. For this purpose, ion composition data from time of flight measurements and test particle simulations to constrain the ions' origin are used. After being validated, the detection method is applied to all the flybys for which the CAPS/IMS instrument gathered valid data, constraining the region around the moon where the signatures are observed. The results reveal an escape region located in the anti-Saturn direction as expected from the nominal corotation electric field direction. These findings provide new constraints for the area of freshly produced pickup ion escape, giving an approximate escape rate of 3.3−2+3×1023ions·s−1.

A Dione O 2 and CO 2 exosphere of similar composition and density to Rhea's is confirmed by Cassini spacecraft Ion Neutral Mass Spectrometer (INMS) flyby data. INMS results from three Dione and two Rhea flybys show exospheric spatial and temporal variability indicative of seasonal exospheres, modulated by winter polar gas adsorption and desorption at the equinoxes. Cassini Plasma Spectrometer (CAPS) pickup ion fluxes also show exospheric structure and evolution at Rhea consistent with INMS, after taking into consideration the anticipated charge exchange, electron impact, and photo-ionization rates. Data-model comparisons show the exospheric evolution to be consistent with polar frost diffusion into the surface regolith, which limits surface exposure and loss of the winter frost cap by sputtering. Implied O 2 source rates of ∼45(7) × 10 21 s −1 at Dione(Rhea) are ∼50(300) times less than expected from known O 2 radiolysis yields from ion-irradiated pure water ice measured in the laboratory, ruling out secondary sputtering as a major exospheric contributor, and implying a nanometer scale surface refractory lag layer consisting of concentrated carbonaceous impurities. We estimate ∼30:1(2:1) relative O 2 :CO 2 source rates at Dione(Rhea), consistent with a stoichiometric bulk composition below the lag layer of 0.01(0.13) C atoms per H 2 O molecule, deriving from endogenic constituents, implanted micrometeoritic organics, and (in particular at Dione) exogenous H 2 O delivery by E-ring grains. Impact deposition, gardening and vaporization may thereby control the global O 2 source rates by fresh H 2 O ice exposure to surface radiolysis and trapped oxidant ejection.

We report fortuitous observations of low-energy lunar pickup ion fluxes near the Moon while in the solar wind by the Lunar Dust Experiment (LDEX) on board the Lunar Atmosphere and Dust Environment Explorer (LADEE). We describe the method of observation and the empirical calibration of the instrument for ion observations. LDEX observes several trends in the exospheric ion production rate, including a scale height of approximately 100 km, a positive, linear correlation with solar wind flux, and evidence of a slight enhancement near 7–8 h local time. We compare the LDEX observations to both LADEE Neutral Mass Spectrometer ion mode observations and theoretical models. The LDEX data are best fit by total exospheric ion production rates of ≈6 × 103 m−3 s−1 with dominant contributions from Al+, CO+, and Ar+, although the LDEX data suggest that the aluminum neutral density and corresponding ion production rate are lower than predicted by recent models.

We investigate a generation mechanism of O+ ion beams observed above the Martian bow shock by analyzing ion velocity distribution functions (VDFs) measured by the Suprathermal and Thermal Ion Composition instrument on the Mars Atmosphere and Volatile Evolution (MAVEN) spacecraft. In the solar wind near Mars, MAVEN often observes energetic O+ ion beams (similar to 10 keV or higher). Accompanied with the O+ ion beam events, we sometimes observe characteristic ion VDFs in the magnetosheath: a partial ring distribution. The partial ring distribution corresponds to pickup ions with a finite initial velocity (i.e., not newborn pickup ions), and its phase space density is much smaller than that of local pickup O+ ions of the magnetosheath. Thus, the partial ring distribution is most likely produced by the reflection of pickup O+ ions precipitating from the upstream solar wind below the bow shock. After being injected into the magnetosheath from the solar wind, the precipitating O+ ions are subject to the significantly enhanced magnetic field in this region and start to gyrate around the guiding center of the plasma frame in the magnetosheath. Consequently, a part of precipitating O+ ions are reflected back to the solar wind, generating O+ beams in the solar wind. The beams direct quasi-sunward near the subsolar region but have large angle with respect to the sunward direction at high solar zenith angles (>50 degrees). The reflected O+ beams are accelerated by the convection electric field of the solar wind and may escape Mars.

Non-thermal pickup ions (PUIs) are created in the solar wind (SW) by charge-exchange between SW ions (SWIs) and slow interstellar neutral atoms. It has long been theorized, but not directly observed, that PUIs should be preferentially heated at quasi-perpendicular shocks compared to thermal SWIs. We present in situ observations of interstellar hydrogen (H+) PUIs at an interplanetary shock by the New Horizons' Solar Wind Around Pluto (SWAP) instrument at ~34 au from the Sun. At this shock, H+ PUIs are only a few percent of the total proton density but contain most of the internal particle pressure. A gradual reduction in SW flow speed and simultaneous heating of H+ SWIs is observed ahead of the shock, suggesting an upstream energetic particle pressure gradient. H+ SWIs lose ~85% of their energy flux across the shock and H+ PUIs are preferentially heated. Moreover, a PUI tail is observed downstream of the shock, such that the energy flux of all H+ PUIs is approximately six times that of H+ SWIs. We find that H+ PUIs, including their suprathermal tail, contain almost half of the total downstream energy flux in the shock frame., 11 pages, 4 figures, accepted for publication in Physical Review Letters

The Rosetta Ion and Electron Sensor (IES) has been measuring solar wind ions intermittently since exiting from hibernation in May 2014. On 19 August, when Rosetta was ~80 km from the comet 67P/Churyumov-Gerasimenko, which was ~3.5 AU from the Sun, IES began to see ions at its lowest energy range, ~4–10 eV. We identify these as ions created from neutral species emitted by the comet nucleus, photoionized by solar UV radiation in the neighborhood of the Rosetta spacecraft (S/C), and attracted by the small negative potential of the S/C resulting from the population of thermal electrons. Later, IES began to see higher-energy ions that we identify as having been picked up and accelerated by the solar wind. IES continues to measure changes in the solar wind and the development of the pickup ion structure.

Interstellar neutrals penetrating into the inner heliosphere are ionized by photoionization, charge exchange with solar wind ions, and electron impact ionization. These processes comprise the first step in the evolution of interstellar pickup ion (PUI) distributions. Typically, PUI distributions have been described in terms of velocity distribution functions that cool adiabatically under solar wind expansion, with a cooling index of 3/2. Recently, the cooling index has been determined experimentally in observations of He PUI distributions with Advanced Composition Explorer (ACE)/Solar Wind Ion Composition Spectrometer and found to vary substantially over the solar cycle. The experimental determination of the cooling index depends on the knowledge of the ionization rates and their spatial variation. Usually, ionization rates increase with 1/r2 as neutral particles approach the Sun, which is not exactly true for electron impact ionization, because the electron temperature increases with decreasing distance from the Sun due to the complexity of its distributions and different radial gradients in temperature. This different dependence on distance may become important in the study of the evolution of PUI distributions and is suspected as one of the potential reasons for the observed variation of the cooling index. Therefore, we investigate in this paper the impact of electron ionization on the variability of the cooling index. We find that the deviation of the electron ionization rate from the canonical 1/r2 behavior of other ionization processes plays only a minor role.

The microstructure of the heliospheric termination shock and the accompanied local acceleration processes of both ions and electrons are investigated by utilizing one-dimensional full particle-in-cell simulations for a variety of parameters. The relative pickup ion density is assumed to be 20–30%. The magnetic field and the shock potential profiles exhibit significant differences, since the former mostly reflects the dynamics of solar wind ions, whereas the latter is mainly sustained by the bulk motion of the reflected pickup ions in the extended foot. The discrepancy between the magnetic field profile and the potential profile increases with Alfven Mach number. Most of the downstream thermal energy is gained by the pickup ions, while some heating of the solar wind ions and electrons occurs through the modified two-stream instability excited in the extended foot. Self-reformation can occur when the relative pickup ion density is 20% but is blurred when it becomes as large as 30%. Reformation is also suppressed if the local solar wind ion temperature in the extended foot is high, which can either be due to heating by the modified two-stream instability or is already determined by the solar wind temperature far upstream. In all runs presented in this study no evidence for shock surfing acceleration of pickup ions could be found. Nonthermal particle acceleration occurs for oblique shocks. Electron (pickup ion) shock drift acceleration is evidenced when the shock angle is below 80° (60°).

In this report we discuss the global plasma environment of the TA flyby from the perspective of 3D hybrid modeling. In our model the background, pickup, and ionospheric ions are considered as particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. Our modeling shows that mass loading of the background plasma ( H + , O + ) by pickup ions H 2 + , CH 4 + and N 2 + differs from the T9 encounter simulations when O + ions are not introduced into the background plasma. In our hybrid modeling we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M =28 amu ions that were generated inside the ionosphere. Titan׳s interior is considered as a weakly conducting body. Special attention has been paid to comparing the simulated pickup ion density distribution with CAPS-ELS and with RPWS LP observations by the Cassini–Huygens spacecraft along the TA trajectory. Our modeling shows an asymmetry of the ion density distribution and the magnetic field, including the formation of Alfven wing-like structures.

In this report we discuss the self-consistent dynamics of pickup ions in the solar wind flow around the lunar-like object. In our model the solar wind and pickup ions are considered as a particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. The Moon will be chosen as a basic object for our modeling. The current modeling shows that mass loading by pickup ions H + , H 2 + , He + , and Na + may be very important in the global dynamics of the solar wind around the Moon. In our hybrid modeling we use exponential profiles for the exospheric components. The Moon is considered as a weakly conducting body. Special attention will be paid to comparing the modeling pickup ion velocity distribution with ARTEMIS observations. Our modeling shows an asymmetry of the Mach cone due to mass loading, the upstream flow density distribution and the magnetic field. The pickup ions form an asymmetrical plasma tails that may disturb the lunar plasma wake.

[1] We use Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun (ARTEMIS) measurements of lunar exospheric pickup ions in the terrestrial magnetotail lobes combined with a particle-tracing model to constrain the source species and distributions of the lunar neutral exosphere. These pickup ions, generated by photoionization of neutral species while the Moon is in the magnetotail lobes, undergo acceleration from both the magnetotail convection electric field and the lunar surface photoelectric field, giving rise to distinct pickup ion flux, pitch angle, and energy distributions. By simulating the behavior of lunar pickup ions in the magnetotail lobes and the response of the twin ARTEMIS probes under various ambient conditions, we can constrain several physical quantities associated with these observations, including the source ion production rate and the magnetotail convection velocity (and hence, electric field). Using the model-derived source ion production rate and established photoionization rates, we present upper limits on the density of several species potentially in the lunar exosphere. In certain cases, these limits are lower than those previously reported. We also present evidence that the lunar exosphere is displaced toward the lunar dawnside while in the terrestrial magnetotail based on fits to the observed pickup ion distributions.

In this report we discuss the self-consistent dynamics of pickup ions in the solar wind flow around the lunar-like object. In our model the solar wind and pickup ions are considered as a particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. The Moon will be chosen as a basic object for our modeling. The current modeling shows that mass loading by pickup ions Na + and He + may be very important in the global dynamics of the solar wind around the Moon. In our hybrid modeling we use exponential profiles for the exospheric components. The Moon is considered as a weakly conducting body. Special attention will be paid to comparing the modeling pickup ion velocity distribution with ARTEMIS observations. Our modeling shows an asymmetry of the Mach cone due to mass loading, the upstream flow density distribution and the magnetic field. The pickup ions form an asymmetrical plasma tails that may disturb the lunar plasma wake.

Global dynamics of ionized and neutral gases in the environment of Titan plays an important role in the interaction of Saturn s magnetosphere with Titan. Several hybrid simulations of this problem have already been done (Brecht et al., 2000; Kallio et al., 2004; Modolo et al., 2007a; Simon et al., 2007a, 2007b; Modolo and Chanteur, 2008). Observational data from CAPS for the T9 encounter (Sittler et al., 2009) indicates an absence of O(+) heavy ions in the upstream that change the models of interaction which were discussed in current publications (Kallio et al., 2004; Modolo et al., 2007a; Simon et al., 2007a, 2007b; Ma et al., 2007; Szego et al., 2007). Further analysis of the CAPS data shows very low density or even an absence of H(+) ions in upstream. In this paper we discuss two models of the interaction of Saturn s magnetosphere with Titan: (A) high density of H(+) ions in the upstream flow (0.1/cu cm), and (B) low density of H(+) ions in the upstream flow (0.02/cu cm). The hybrid model employs a fluid description for electrons and neutrals, whereas a particle approach is used for ions. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. The model atmosphere includes exospheric H(+), H(2+), N(2+)and CH(4+) pickup ion production as well as an immobile background ionosphere and a shell distribution for active ionospheric ions (M(sub i)=28 amu). The hybrid model allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of Alfven wing-like structures. The results of the ion dynamics in Titan s environment are compared with Cassini T9 encounter data (CAPS).

In this report we discuss the ion velocity distribution dynamics from the 3D hybrid simulation. In our model the background, pickup, and ionospheric ions are considered as a particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. The current simulation shows that mass loading by pickup ions H(+); H2(+), CH4(+) and N2(+) is stronger than in the previous simulations when O+ ions are introduced into the background plasma. In our hybrid simulations we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M = 28 amu ions that were generated inside the ionosphere. The moon is considered as a weakly conducting body. Special attention will be paid to comparing the simulated pickup ion velocity distribution with CAPS T9 observations. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of the Alfve n wing-like structures. The simulation also shows that the ring-like velocity distribution for pickup ions relaxes to a Maxwellian core and a shell-like halo.

Titan has a very extensive atmosphere and exosphere which interact strongly with the corotating magnetospheric plasma. Some of the new pickup ions created in the vicinity of the exosphere will re-impact the upper atmosphere causing additional energy input. Previous investigations of the atmospheric collisional interaction of the pickup ions have generally assumed the atmosphere to be spherically symmetric. From the Voyager and Cassini observations, we know that the magnetic field configuration and plasma flow field is highly asymmetric. To study the possible spatial variation in the pickup ion influx, we have employed the plasma data of the three dimensional MHD simulation of Kopp and Ip [Kopp, A., Ip, W.-H. Asymmetric mass loading effect at Titan’s ionosphere. J. Geophys. Res. 106, 8323–8332, 2001] to compute the trajectories of the pickup ions. In this work, we calculate the ion influx and energy deposit into Titan’s exobase for the H2+, CH4+ and N2+ pickup ions separately. The model results of four different Titan’s orbital locations are also presented.

The Cassini plasma spectrometer (CAPS) instrument made measurements of Titan's plasma environment when the Cassini Orbiter flew through the moon's plasma wake October 26, 2004 (flyby TA). Initial CAPS ion and electron measurements from this encounter will be compared with measurements made by the Voyager 1 plasma science instrument (PLS). The comparisons will be used to evaluate previous interpretations and predictions of the Titan plasma environment that have been made using PLS measurements. The plasma wake trajectories of flyby TA and Voyager 1 are similar because they occurred when Titan was near Saturn's local noon. These similarities make possible direct, meaningful comparisons between the various plasma wake measurements. They lead to the following: (A) The light and heavy ions, H+and N+/O+, were observed by PLS in Saturn's magnetosphere in the vicinity of Titan while the higher mass resolution of CAPS yielded H+ and H2+as the light constituents and O+/CH4+ as the heavy ions. (B) Finite gyroradius effects were apparent in PLS and CAPS measurements of ambient O+ ions as a result of their absorption by Titan's extended atmosphere. (C) The principal pickup ions inferred from both PLS and CAPS measurements are H+, H2+, N+, CH4+ and N2+. (D) The inference that heavy pickup ions, observed by PLS, were in narrow beam distributions was empirically established by the CAPS measurements. (E) Slowing down of the ambient plasma due to pickup ion mass loading was observed by both instruments on the anti-Saturn side of Titan. (F) Strong mass loading just outside the ionotail by a heavy ion such as N2+ is apparent in PLS and CAPS measurements. (G) Except for the expected differences due to the differing trajectories, the magnitudes and structures of the electron densities and temperatures observed by both instruments are similar. The high-energy electron bite-out observed by PLS in the magnetotail is consistent with that observed by CAPS.

A two-dimensional kinetic model calculation for the water group species (H2O, H2, O2, OH, O, H) in Europa's atmosphere is undertaken to determine its basic compositional structure, gas escape rates, and velocity distribution information to initialize neutral cloud model calculations for the most important gas tori. The dominant atmospheric species is O2 at low altitudes and H2 at higher altitudes with average day–night column densities of 4.5 × 10 14 and 7.7 × 10 13 cm −2 , respectively. H2 forms the most important gas torus with an escape rate of ∼ 2 × 10 27 s −1 followed by O with an escape rate of ∼ 5 × 10 26 s −1 , created primarily as exothermic O products from O2 dissociation by magnetospheric electrons. The circumplanetary distributions of H2 and O are highly peaked about the satellite location and asymmetrically distributed near Europa's orbit about Jupiter, have substantial forward clouds extending radially inward to Io's orbit, and have spatially integrated cloud populations of 4.2 × 10 33 molecules for H2 and 4.0 × 10 32 atoms for O that are larger than their corresponding populations in Europa's local atmosphere by a factor of ∼200 and ∼1000, respectively. The cloud population for H2 is a factor of ∼3 times larger than that for the combined cloud population of Io's O and S neutral clouds and provides the dominant neutral population beyond the so-called ramp region at 7.4–7.8 R J in the plasma torus. The calculated brightness of Europa's O cloud on the sky plane is very dim at the sub-Rayleigh level. The H2 and O tori provide a new source of europagenic molecular and atomic pickup ions for the thermal plasma and introduce a neutral barrier in which new plasma sinks are created for the cooler iogenic plasma as it is transported radially outward and in which new sinks are created to alter the population and pitch angle distribution of the energetic plasma as it is transported radially inward. The europagenic instantaneous pickup ion rates are peaked at Europa's orbit, dominate the iogenic pickup ion rates beyond the ramp region, and introduce new secondary plasma source peaks in the solution of the plasma transport problem. The H2 torus is identified as the unknown Europa gas torus that creates both the observed loss of energetic H+ ions at Europa's orbit and the corresponding measured ENA production rate for H.

The Ultraviolet Coronograph on board SOHO (UVCS) has observed the 58.4 nm glow of the interplanetary He focusing cone at regular intervals since 1996. The intensity decrease with time already observed during the first two years (Michels et al. 2002) has dramatically amplified during the solar activity increase. Intensities seem to reach a plateau in 2001. Using a model of the cone emission which takes into account both photoionization and electron impact ionization of neutral helium we show that the photoionization increase alone cannot explain the observed intensity drop. Data can be fitted if at minimum activity the electron impact ionization rate is the solar cycle average rate predicted by Rucinski & Fahr (1989), and if this rate is increased by a factor of about 3.5 between 1996 and 2001. Assuming the Rucinski and Fahr radial dependence, such high electron impact rates create averaged He + pickup ion (PUI) fluxes which may reach 50% of the fluxes of ions born after photoionization, or 35% of PUI total fluxes, as far as 1 AU from the Sun. In slow and dense solar wind enhancements, in particular in the presence of strong suprathermal tails, PUIs from electron impact could be the dominant species. This could explain a fraction of the observed correlation between He + and H + pickups and anticorrelation of He + fluxes with solar wind velocity.

Maintaining the stability of pickup ions in the outer heliosheath is a critical element for the secondary energetic neutral atom (ENA) mechanism, a theory put forth to explain the nearly annular band of ENA emission observed by the Interstellar Boundary EXplorer. A recent study showed that a pickup ion ring can remain stable to the Alfven/ion cyclotron (AC) instability at propagation parallel to the background magnetic field when the parallel thermal spread of the ring is comparable to that of a background population. This study investigates the potential role that the mirror or ion Bernstein (IB) instabilities can play in the stability of pickup ions when conditions are such that the AC instability is suppressed. Linear Vlasov theory predicts relatively fast mirror and IB instability growth even though AC instability growth is suppressed. For a few such cases, two-dimensional hybrid and macroscopic quasi-linear simulations are carried out to examine how the unstable mirror and IB modes evolve and affect the pickup ion ring beyond the linear theory picture. For the parameters used, the mirror mode dominates initially and leads to a rapid parallel heating of the pickup ions in excess of the parallel temperature of the background protons. The heated pickup ions subsequently trigger onset of the AC mode, which grows sufficiently large to be the dominant pitch angle scattering agent after the mirror mode has decayed away. The present results indicate that the pickup ion stability needed may not be guaranteed once the mirror and IB instabilities are taken into account.

A large number of numerical experiments has been performed in order to study the interaction of interstellar pickup protons and helium ions with quasiparallel collisionless shocks. The shocks are modeled by a one-dimensional hybrid simulation method which treats the ions as macroparticles and the electrons as a massless fluid. Solar wind alpha particles and pickup protons are included self-consistently. In addition, the particle splitting method is used for the solar wind ions so that the distribution function can be followed over more than 10 orders of magnitude. A large part of the pickup ion distribution is reflected; the reflection efficiency is very high, and can reach in cases where the pickup ion density is low as much as 50%–60%. The reflection efficiency is almost independent of magnetic field–shock normal angle. This indicates that magnetic mirroring is unimportant and does not lead to larger reflection efficiencies. The reflection efficiency of pickup protons rapidly decreases when the pickup i...

The effects of Alfven cyclotron waves generated due to the formation in the outer heliosphere of pickup ions on the transport coefficients of low-energy electrons is investigated here. To this end, parallel mean free path (MFP) expressions are derived from quasilinear theory, employing the damping model of dynamical turbulence. These are then used as inputs for existing expressions for the perpendicular MFP and turbulence-reduced drift coefficient. Using outputs generated by a two-component turbulence transport model, the resulting diffusion coefficients are compared with those derived using a more typically assumed turbulence spectral form, which neglects the effects of pickup ion-generated waves. It is found that the inclusion of pickup ion effects greatly leads to considerable reductions in the parallel and perpendicular MFPs of 1–10 MeV electrons beyond ~10 au, which are argued to have significant consequences for studies of the transport of these particles.

In this work, we discuss the properties of the ion cyclotron waves observed by the Galileo spacecraft in the Io torus on December 7, 1995. These waves are generated by anisotropies in the pickup ion distributions present in the region. To zeroth order the waves grow at the sulfur dioxide gyrofrequency, are left-hand polarized and propagate along the magnetic field. A more detailed study reveals that in some regions the waves grow over a finite bandwidth, with frequencies near the SO + and S + gyrofrequencies, propagate at an angle to the field, and are elliptically polarized. We perform a kinetic dispersion analysis and show that while the parallel propagating SO 2 + cyclotron mode has the largest growth, oblique cyclotron waves can also grow with a significant rate in the multicomponent plasma near Io. We find that decreasing the sulfur dioxide pickup ion density, and including a SO + pickup component can lead to instabilities generated by the sulfur monoxide and sulfur components in the plasma, and to the observed bandwidth of frequencies. Therefore, it is possible that in addition to the waves generated by SO 2 + pickup ions, in some regions Galileo observed ion cyclotron waves generated by anisotropic distributions of SO + and S + ions. When Voyager 1 flew by the torus in 1979 the spacecraft detected the Alfven wing but no cyclotron waves were observed in the data. We evaluate the group and phase velocities of ion cyclotron waves and find that both are smaller than the Alfven velocity. Then it is possible that the cyclotron waves can get carried away before they reach Voyager.

Flowing plasma interactions with small bodies in the solar system are frequently complicated by the significant size of the ion gyroradii relative to the scale size of the body and the interaction region. This situation often applies to both the flowing ion populations and the atmospheric pickup ions added to the flow in the vicinity of the body. Titan represents an important example of this type of system and is scheduled to be probed in some detail during the Cassini mission tour of Saturn. Initial results from global hybrid (particle ions, fluid electrons) numerical simulations of Saturn's magnetospheric plasma encountering Titan's ionosphere reveal the complexity of the interaction introduced by the ion kinetics. The scale of the interaction region is dominated by the heavy ion gyroradii of the ambient and pickup ions rather than the size of Titan and is found to depend on the mass loading of the magnetospheric flow encountering Titan. Comparison with data from the Voyager 1 flyby suggests that ion kinetics as well as the relative orientations of external plasma ram and pickup ion production regions affect what was observed.

Ion pickup processes invariably produce anisotropic ion distributions and hence can give rise to wave generation. This paper discusses the waves due to the Io pickup source in Jupiter's magnetosphere. In the near-Io torus, the Galileo spacecraft observed ion cyclotron waves near the sulfur dioxide gyrofrequency, which grow due to the absence of a thermalized background component of these molecular ions in a torus plasma of predominantly dissociated species. These wave observations allow limits on pickup source rates to be estimated. On the edge of the Io wake, mirror mode waves were seen where the pickup ion contribution to the pressure anisotropy is sufficient to overcome the instability threshold. The mirror mode can dominate over the individual ion cyclotron modes in a multi-species plasma. The mass density in the torus cannot increase indefinitely and the plasma must be radially transported through the Jovian magnetosphere. Wave processes may be associated with any ion distribution anisotropies that arise during their radial transport, and with unsteady reconnection processes that may facilitate their eventual loss from the magnetotail.

Multiscale hybrid kinetic simulations of low-{beta}{sub p} supercritical shocks demonstrate that pickup ions may be strongly accelerated by shock surfing, also known as multiply reflected ion acceleration. {copyright} {ital 1999} {ital The American Physical Society}

The time-dependent cosmic ray modulation due to solar activity and global merged interaction regions is studied with a self-consistent model of the interaction of a time-varying solar wind, including pickup ions, with cosmic rays. After a study of the spatial evolution of global merged interaction regions, their different effects on local as well as global modulation of both galactic and anomalous cosmic rays are analyzed. It is found that the long-term modulation cannot be explained alone with successive step decreases in the cosmic ray intensities caused by global merged interaction regions, but that an additional long-term global variation in spatial diffusion is required. The analysis for a solar wind with pickup ions reveals that there occurs a significant upstream expansion of global merged interaction regions accompanied by a decrease of their peak density enhancement. Both effects lead to an immediate recovery of the cosmic ray intensities after a step decrease. Furthermore, the calculated response of the solar wind termination shock to the passage of global merged interaction regions shows minor changes in the structure and position so that the shock remains in a gas-dynamic (pickup ion) dominated state; that is, it is never cosmic ray dominated. This investigation also shows that the decay of global merged interaction regions is particularly fast during their interaction with the termination shock and in the region beyond, so they probably cannot be interpreted as the downstream modulation barriers discussed in the literature.

We have performed a number of numerical experiments to study the interaction of interstellar pickup protons and helium ions with quasi-parallel stationary (bow) shocks and interplanetary traveling shocks. The collisionless shocks are modeled by the hybrid simulation method which treats the ions as macroparticles and the electrons as a massless background fluid. Solar wind protons, alpha particles, pickup H+ and He+ are included self-consistently. The pickup ion distributions are modeled as spherical shells in velocity space convected with the solar wind. Pickup ions are rather efficiently reflected at both types of shocks. The reflection coefficient, measured as the ratio of incident to reflected ions, is more than one order of magnitude larger than the reflection coefficients for either solar wind protons or alpha particles, and is of the order of 30%. We have found that at high Mach number shocks pickup ions can reach many tens to over hundred times the plasma ramming energy in essentially one shock encounter. The reflected pickup ion distributions have a pancake-type shape in velocity space.

On the basis of recent observations that pickup ions are inefficiently scattered through 90° pitch angle, a pickup ion transport model has been developed which assumes that pickup ion distributions are uniform in the so-called forward and reverse hemispheres, containing ions with pitch angles less than 90° and greater than 90°, respectively. A simple, fast numerical scheme is then described which solves the transport equations. Several examples are solved and comparisons are performed between the numerical solutions and previously used analytical solutions. Results from the model are then compared to pickup ion spectra obtained with the solar wind ion composition spectrometer on the Ulysses satellite.

We present a simple model of anisotropic pickup ion distributions on a flux tube containing both a radial segment and segments oriented perpendicular to the radial direction. The flux tube is carried away from the Sun by a constant speed solar wind. We model the pitch-angle scattering of these ions by a hemispherical approximation, which assumes that the scattering is instantaneous except through the pitch angle of 90°, where the scattering is finite. The model includes a radially dependent ionization of pickup ions, with the newly created ions concentrated in the sunward hemisphere within the radial segment but distributed equally between the hemispheres in the perpendicular segments. The pickup ions stream along the magnetic field and decelerate in the expanding solar wind, but we assume that the radial segment is short enough that the effects of adiabatic focusing may be neglected. For weak scattering through 90°, the distributions exhibit nonlocal properties due to ion streaming along the field from regions where conditions are different than observed locally. We find that this model can yield a depletion of antisunward ions in the radial segment of the flux tube, such as has been observed, but the streaming of ions into the radial region limits the size of this effect unless the radial segment is very long.

During the Galileo-Io flyby, nearly field-aligned, left-hand circularly polarized ion cyclotron waves were observed in a band near the sulfur-dioxide ion gyrofrequency. We have performed a warm plasma dispersion analysis using nominal Io torus composition ratios, pickup ion “ring”-type velocity-space distributions, and a thermal background plasma of typical torus temperature. Analysis shows that the SO2+ wave is dominant, particularly as the ring begins to broaden. The observed spectral peak indicates that the ring distribution of pickup SO2+ is highly unstable and generates waves, but there is unlikely to be sufficient time for ions to fully thermalize before dissociation occurs. Assuming wave-particle scattering of ions towards a “bispherical” shell-type distribution, a free energy analysis and comparison with observed wave amplitudes suggests that the ions are not scattered far from the ring and/or that the SO2+ composition ratio in the Io torus falls off considerably from Io wake values.

Pickup ions are expected to play a fundamental role in determining the global structure of the heliosphere since they dominate the internal energy of the outer solar wind beyond the ionization cavity. The dynamical effect of pickup ions on the propagation characteristics of interplanetary shocks is addressed here. It is found that (1) forward and reverse shocks in a heliosphere mediated by pickup ions have higher speeds than expected of shocks in an adiabatic solar wind without pickup ions; (2) identical initial conditions evolve differently in a pickup ion mediated heliosphere compared to an adiabatic heliosphere, producing more complicated velocity, density and temperature structure, and often with multiple forward-reverse shock pairs; (3) pickup ions act to weaken and dissipate shocks in the outer heliosphere although the velocity jumps can remain large; (4) density fluctuations in the outer pickup ion mediated heliosphere do not correlate with velocity fluctuations as well as they do in the inner heliosphere, and (5) the amplitudes of density fluctuations in a pickup ion mediated heliosphere are considerably smaller than those found in an adiabatic heliosphere without pickup ions.

According to one possible picture of the Titan atmosphere interaction with Saturn's magnetosphere, newborn exospheric pickup ions may exhibit a marked asymmetry with respect to the corotation wake due to their large gyroradii in the local ∼5 nT magnetic field. These finite gyroradius effects should persist until the ions are scattered by wave-particle interactions (a scale length that is presently unknown). To the extent that simple test-particle-like motion is maintained, the combination of the small (on the scale of most ion gyroradii) size of Titan and the small, steady local magnetic field at ∼20 Saturn Radii provides a natural ion mass spectrometer configuration. The resulting spatial dispersion of the pickup ion trajectories can be used for both analyzing the composition of the Titan exosphere and determining the ionization rates of its various constituents. Possible evidence of this ion mass spectrometer effect may be visible as a magnetic field perturbation in the Voyager 1 Titan flyby data.

Pickup ions measured deep inside the heliosphere open a new way to determine the absolute atomic density of a number of elements and isotopes in the local interstellar cloud (LIC). We derive the atomic abundance of hydrogen and the two isotopes of helium from the velocity and spatial distributions of interstellar pickup protons and ionized helium measured with the Solar Wind Ion Composition Spectrometer (SWICS) on the Ulysses spacecraft between ~2 and ~5 AU. The atomic hydrogen density near the termination shock derived from interstellar pickup ion measurements is 0.115±0.025 cm-3, and the atomic H/He ratio from these observations is found to be 7.7 ± 1.3 in the outer heliosphere. Comparing this value with the standard universal H/He ratio of ~10 we conclude that filtration of hydrogen is small and that the ionization fraction of hydrogen in the LIC is low.

Data from the Shuttle Potential and Return Electron Experiment (SPREE) flown as part of the Tethered Satellite System (TSS 1) are used to determine the detailed characteristics of beam arc distributions of pickup ions due to molecules outgassed or ejected from the shuttle. These ion distributions are only detected near the plane perpendicular to the magnetic field direction in an angular range of ±45° about the minimum angle to ram. Their flux is largest when the angle between this plane and the shuttle ram direction is smallest. Generally, ion spectra peak in the range 19 to 25 eV at the minimum angle between the perpendicular plane and the ram direction. The peak energy decreases smoothly as this angle increases. Weak fluxes are measured above the peak, to energies as high as 150 eV. Within the SPREE energy range, two-dimensional distribution functions of beam arc ions in the perpendicular plane have teardrop shapes, symmetric about the minimum angle to shuttle ram with deep minima in the centers. Variations in the peak energies of differential number fluxes agree with collisionless trajectory analysis, assuming that the ions are H2O+ and allowing for different initial velocities before charge exchange. The lowest densities for beam arc ions occur during periods of purely residual outgassing from the shuttle. Ion densities increase by a factor of 5 during waste water dumps. The highest densities occur during operation of the flash evaporator system when the pickup ion densities in daylight can exceed 105 ions cm−3, about 30% of the estimated ambient plasma density. We also present a nonlinear numerical analysis to study the stability of beam arc generated plasma waves and explain electrostatic spectra measured during previous shuttle flights.

Over the last 10 years, the experimental basis for the study of the local interstellar medium has been substantially enhanced by the direct detection of interstellar pickup ions and of interstellar neutral helium within the heliosphere. Pickup ions can be studied for a wide range of interstellar species. However, currently the accuracy of the method to determine the parameters of the interstellar medium, namely neutral density, temperature and relative velocity, is hampered by two problems: (1) In most cases the crucial ionization rates are not available from simultaneous measurements and (2) the transport of the pickup ions in the interplanetary medium substantially modifies the measured spatial distribution of the ions. In this study we will discuss how the enhanced capabilities of the instrumentation on SOHO and Cluster in combination with ongoing efforts to model the pickup ion distributions will lead to a significant improvement over the coming years.

We investigate the properties of a bispherical shell distribution of pickup ions when dispersion of the resonant waves is taken into account. We also extend the method of Johnstone et al. [1991] and Huddleston and Johnstone [1992] for determining the spectra of the quasi-linear self-generated waves to include dispersion. Our specific results assume that all waves propagate along the average magnetic field, and that the waves satisfy the cold electron-proton dispersion relation. The major effect of including dispersion in the analysis is that the particles in a bispherical distribution retain more of their initial energy than in the nondispersive treatment, so the resulting wave spectra are less intense than predicted when dispersion is neglected. The dispersive analysis also encounters the well-known gap in the cyclotron resonance with fast-mode waves which is absent in the nondispersive picture. We present detailed results for pickup protons in the two extreme configurations of magnetic field perpendicular to, and parallel to, the solar wind flow for a range of values of VA/VSW. We find that when pickup protons are scattered to a bispherical distribution with VA/VSW = 0.1, typical of conditions in the solar wind, the discrepancy in the total energy density of the self-generated waves between the dispersive and nondispersive results amounts to 18% in the perpendicular configuration. This discrepancy falls between 9% and 73% in the parallel configuration, depending on the extent of the particle transport through the resonance gap. These discrepancies are substantially larger for larger values of VA/VSW. These results indicate that dispersive effects can be important in the correct modeling of the quasi-linear resonant wave-particle interaction between pickup ions and self-generated waves in the solar wind.

Pickup ions in a ring velocity distribution are unstable to several kinetic plasma instabilities. At large heliocentric distances where the overall plasma β (ratio of kinetic to magnetic energy) is dominated by the energy density of interstellar pickup ions and pickup is perpendicular to the interplanetary magnetic field, the dominant of these is the Alfven ion cyclotron instability (AIC). We demonstrate by hybrid particle simulation that, for conditions where the solar wind β is low, AIC driven by the pickup ions couples to the solar wind. The result is perpendicular heating, leading to an anisotropic solar wind distribution. This process may contribute to enhanced solar wind temperatures at large heliocentric distances and may allow for indirect measurement of interstellar pickup ions.

Pickup ions are created when interstellar neutral atoms resonantly exchange charge with the solar wind (SW) ions, especially in the supersonic part of the wind, where they carry most of the plasma pressure. Here we present numerical simulation results of the 3D heliospheric interface treating pickup ions as a separate proton fluid. To satisfy the fundamental conservation laws, we solve the system of equations describing the flow of the mixture of electrons, thermal protons, and pickup ions. To find the density and pressure of pickup ions behind the termination shock, we employ simple boundary conditions that take into account the \emph{Voyager} observations that showed that the decrease in the kinetic energy of the mixture at the termination shock predominantly contributed to the increase in the pressure of pickup ions. We show that this model adequately describes the flow of the plasma mixture and results in a noticeable decrease in the heliosheath width., Comment: Accepted for publication in Journal of Physics, Conference Series

We present results from hybrid simulations (kinetic ion/fluid electron) of the interaction of interstellar pickup ions with collisionless shocks. Since cross-field transport is unphysically suppressed in the one-dimensional geometry used here, an ad hoc scattering algorithm is used to model this effect. This is a necessary step to accelerate the pickup ions from their initial low energies at quasi-perpendicular shocks to the high energies which are often observed associated with traveling interplanetary shocks by Ulysses.

[1] Observations of interstellar pickup ions inside the heliosphere provide an indirect method to access information on the surrounding interstellar medium. The so-called pickup ion focusing cone and pickup ion crescent, which show an imprint of the related longitudinal distribution of interstellar neutrals in form of two overabundances on the down- and upwind side of the sun, are both believed to be aligned along the inflow vector of the interstellar medium. By finding their longitudinal positions, we can give an accurate value for the inflow direction λISM of interstellar matter. For that we performed an epoch analysis of interstellar pickup ions measured by the PLAsma and SupraThermal Ion Composition instrument (PLASTIC) on the Solar TErrestrial RElations Observatory mission (STEREO) and were able to reveal in situ the longitudinal distribution of interstellar He+, O+, and Ne+ pickup ions in the ecliptic plane at 1 AU. The previously accepted values for the inflow direction of interstellar matter in ecliptic longitude, as obtained with Ulysses/GAS (λ = 75.4° ± 0.5°), Prognoz 6 (λ = 74.5° ± 1°), and ACE/SWICS (λ = 74.43° ± 0.33°), are currently debated, especially in view of recent results from the Interstellar Boundary Explorer (IBEX) mission that show an inflow direction of interstellar neutral helium of λ = 79° + 3.0°(−3.5°). Four years of data collected with PLASTIC aboard STEREO A provided statistics sufficient not only to obtain values for the inflow direction of interstellar helium (λCone = 77.4° ± 1.9° and λCrescent = 80.4° ± 5.4°, deduced from an analysis of the He+ focusing cone and crescent, respectively) but also to derive values for the inflow direction of interstellar neon (λCone = 77.4° ± 5.0° and λCrescent = 79.7° ± 2.6°) and oxygen (λCrescent = 78.9° ± 3.1°). Although our values for He+, O+, and Ne+ are consistent with results from ACE, Ulysses, and Prognoz 6, considering the statistical and systematic uncertainties (except λNe,Crescent), they are systematically larger than the previously accepted values of 74.99 ± 0.55° and show a better agreement with the values from IBEX.

ARTEMIS observes pickup ions around the Moon, at distances of up to 20,000 km from the surface. The observed ions form a plume with a narrow spatial and angular extent, generally seen in a single energy/angle bin of the ESA instrument. Though ARTEMIS has no mass resolution capability, we can utilize the analytically describable characteristics of pickup ion trajectories to constrain the possible ion masses that can reach the spacecraft at the observation location in the correct energy/angle bin. We find that most of the observations are consistent with a mass range of approx. 20-45 amu, with a smaller fraction consistent with higher masses, and very few consistent with masses below 15 amu. With the assumption that the highest fluxes of pickup ions come from near the surface, the observations favor mass ranges of approx. 20-24 and approx. 36-40 amu. Although many of the observations have properties consistent with a surface or near-surface release of ions, some do not, suggesting that at least some of the observed ions have an exospheric source. Of all the proposed sources for ions and neutrals about the Moon, the pickup ion flux measured by ARTEMIS correlates best with the solar wind proton flux, indicating that sputtering plays a key role in either directly producing ions from the surface, or producing neutrals that subsequently become ionized.

[1] Observations upstream of Mars show strong and persistent plasma wave activity near the local proton cyclotron frequency. These waves are thought to be generated by unstable populations of newborn ions produced by ionization of the Martian hydrogen exosphere. It is generally believed that the amplitudes of the waves are related to the energy of the local pickup ion population, such that the observed spatial and temporal distribution of the wave amplitudes could be used as a diagnostic of local exospheric structure and loss rate. To understand the relationship between wave amplitudes and ion pickup rate, we carry out a parametric study using the 1D hybrid simulation (kinetic ions, fluid electrons) for varying pickup rates and pickup geometry for the Martian planetary environment. Our results indicate that, contrary to expectations, the observed local ion cyclotron wave amplitudes cannot be easily related to local ion pickup rates because the growth time of the instability is long compared to the transit time of the waves past Mars. Additionally, we find the amount of pickup ion energy lost to wave growth is not a constant over all injection rates and wave amplitudes are likely affected by propagation through regions of spatially non-uniform ion production. Because the waves in the upstream region are in a state of growth rather than saturated, the simulations suggest the pickup ion distributions are likely only partially scattered.

We outline the disagreement between the theoretical expectation of rapid isotropization for interstellar pickup ions in the solar wind and the observations of frequent anisotropy. This disagreement is related to the further problem of explaining the scarcity of the wave enhancements predicted to accompany the isotropization. We describe a computational investigation of the kinetic processes involved, and our plans to resolve these discrepancies. We present results of a preliminary calculation to illustrate the capabilities of this technique.

Pickup ions formed from ionized neutral exospheres in flowing plasmas have phase space distributions that reflect their source's spatial distributions. Phase space distributions of the ions are derived from the Vlasov equation with a delta function source using three.dimensional neutral exospheres. The ExB drift produced by plasma motion picks up the ions while the effects of magnetic field draping, mass loading, wave particle scattering, and Coulomb collisions near a planetary body are ignored. Previously, one.dimensional exospheres were treated, resulting in closed form pickup ion distributions that explicitly depend on the ratio rg/H, where rg is the ion gyroradius and H is the neutral scale height at the exobase. In general, the pickup ion distributions, based on three.dimensional neutral exospheres, cannot be written in closed form, but can be computed numerically. They continue to reflect their source's spatial distributions in an implicit way. These ion distributions and their moments are applied to several bodies, including He(+) and Na(+) at the Moon, H(+2) and CH(+4) at Titan, and H+ at Venus. The best places to use these distributions are upstream of the Moon's surface, the ionopause of Titan, and the bow shock of Venus.

[1] The Solar Wind Ion Detectors (SWIDs) on the Chang'E-1 spacecraft, while orbiting the Moon, occasionally observed two continuous flux peaks with energies not exceeding 8 and 4 times that of the prevailing solar wind proton energy. These form parallel curves (PCs) with an energy ratio of 2 in the energy-time spectrogram. The fluxes of the two curves are comparable, around 10−5 ∼ 10− 4 of the solar wind flux. The pitch angle distribution of PC particles is concentrated around 90°. The velocity space distribution of PC particles shows distinct double-ring feature, suggesting the existence of a pickup ion species with m/q = 2. Pickup ions from local interstellar medium, the inner sources and the lunar exosphere are investigated. We conclude that this observation may be the first in situ evidence for H2+ ions in the lunar exosphere, thus providing new insights on the evolution and fate of solar wind hydrogen in the solar system.

The cometary ion pickup process along the sun-comet line at Comet Halley is investigated using a quasi-linear diffusion model including both pitch angle and energy diffusion, adiabatic compression, and convective motion with the solar wind flow. The model results are compared with energetic ion distributions observed by instruments on board the Giotto spacecraft. The observed power spectrum index of magnetic turbulence (gamma) is 2-2.5. The present simulation shows that when gamma was 2, the calculated proton distributions were much more isotropic than the observed ones. The numerical solutions of the quasi-linear diffusion equations show that the isotropization of the pickup ion distribution, particularly at the pickup velocity, is not complete even close to the bow shock. Given the observed turbulence level, quasi-linear theory yields pickup ion energy distributions that agree with the observed ones quite well and easily produces energetic ions with energies up to hundreds of keV.