Jovian plasma torus interaction with Europa: 3D hybrid kinetic simulation. First results

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa moon-magnetosphere system with respect to variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo orbiter mission, and for planning flyby and orbital measurements (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa ( [Cassidy et al., 2007] and [Shematovich et al., 2005]). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyroradius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions). Non-thermal distributions of upstream plasma will be addressed in future work. Photoionization, electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider two models for background plasma: (a) with O⁺⁺ ions; (b) with O⁺⁺ and S⁺⁺ ions. The majority of O₂ atmosphere is thermal with an extended cold population (Cassidy et al., 2007). A few first simulations already include an induced magnetic dipole; however, several important effects of induced magnetic fields arising from oceanic shell conductivity will be addressed in later work.     

Cross-scale coupling at a perpendicular collisionless shock

A full particle simulation study is carried out on a perpendicular collisionless shock with a relatively low Alfven Mach number (M_A = 5). Recent self-consistent hybrid and full particle simulations have demonstrated ion kinetics are essential for the non-stationarity of perpendicular collisionless shocks, which means that physical processes due to ion kinetics modify the shock jump condition for fluid plasmas. This is a cross-scale coupling between fluid dynamics and ion kinetics. On the other hand, it is not easy to study cross-scale coupling of electron kinetics with ion kinetics or fluid dynamics, because it is a heavy task to conduct large-scale full particle simulations of collisionless shocks. In the present study, we have performed a two-dimensional (2D) electromagnetic full particle simulation with a “shock-rest-frame model”. The simulation domain is taken to be larger than the ion inertial length in order to include full kinetics of both electrons and ions. The present simulation result has confirmed the transition of shock structures from the cyclic self-reformation to the quasi-stationary shock front. During the transition, electrons and ions are thermalized in the direction parallel to the shock magnetic field. Ions are thermalized by low-frequency electromagnetic waves (or rippled structures) excited by strong ion temperature anisotropy at the shock foot, while electrons are thermalized by high-frequency electromagnetic waves (or whistler mode waves) excited by electron temperature anisotropy at the shock overshoot. Ion acoustic waves are also excited at the shock overshoot where the electron parallel temperature becomes higher than the ion parallel temperature. We expect that ion acoustic waves are responsible for parallel diffusion of both electrons and ions, and that a cross-scale coupling between an ion-scale mesoscopic instability and an electron-scale microscopic instability is important for structures and dynamics of a collisionless perpendicular shock.     

The variation of Io's auroral footprint brightness with the location of Io in the plasma torus

Ultraviolet and near-infrared observations of auroral emissions from the footprint of Io's magnetic Flux Tube (IFT) mapping to Jupiter's ionosphere have been interpreted via a combination of the unipolar inductor model [Goldreich, P., Lynden-Bell, D., 1969. Astrophys. J. 156, 59-78] and the multiply-reflected Alfvén wave model [Belcher, J.W., 1987. Science 238, 170-176]. While both models successfully explain the general nature of the auroral footprint and corotational wake, and both predict the presence of multiple footprints, the details of the interaction near Io are complicated [Saur, J., Neubauer, F.M., Connerney, J.E.P., Zarka, P., Kivelson, M.G., 2004. In: Bagenal, F., Dowling, T.E., McKinnon, W.B. (Eds.), Jupiter: The Planet, Satellites and Magnetosphere. Cambridge University Press, Cambridge, UK, pp. 537-560; Kivelson, M.G., Bagenal, F., Kurth, W.S., Neubauer, F.M., Paranicas, C., Saur, J., 2004. In: Bagenal, F., Dowling, T.E., McKinnon, W.B. (Eds.), Jupiter: The Planet, Satellites and Magnetosphere. Cambridge University Press, Cambridge, UK, pp. 513-536]. The auroral footprint brightness is believed to be a good remote indicator of the strength of the interaction near Io, indicating the energy and current strength linking Io with Jupiter's ionosphere. The brightness may also depend in part on local auroral acceleration processes near Jupiter. The relative importance of different physical processes in this interaction can be tested as Jupiter's rotation and Io's orbital motion shift Jupiter's magnetic centrifugal equator past Io, leading to longitudinal variations in the plasma density near Io and functionally different variations in the local field strength near Jupiter where the auroral emissions are produced. Initial HST WFPC2 observations found a high degree of variability in the footprint brightness with time, and some evidence for systematic variations with longitude [Clarke, J.T., Ben Jaffel, L., Gérard, J.-C., 1998. J. Geophys. Res. 103, 20217-20236], however the data were not of sufficient quality to determine functional relationships. In this paper we report the results from a second, more thorough study, using a series of higher resolution and sensitivity HST STIS observations and a model for the center to limb dependence of the optically thin auroral emission brightness based on measurements of the auroral curtain emission distribution with altitude. A search for correlations between numerous parameters has revealed a strong dependence between Io's position in the plasma torus and the resulting footprint brightness that persists over several years of observations. The local magnetic field strength near Jupiter (i.e. the size of the loss cone) and the expected north/south asymmetry in auroral brightness related to the path of currents generated near Io through the plasma torus en route to Jupiter appear to be less important than the total plasma density near Io. This is consistent with the near-Io interaction being dominated by collisions of corotating plasma and mass pickup, a long-standing view which has been subject to considerable debate. The brightness of the auroral footprint emissions, however, does not appear to be proportional to the incident plasma density or energy, and the interpretation of this result will require detailed modeling of the interaction near Io.     

The kinetic nature of turbulence at short scales in the solar wind

The analysis of the transition from the large-scale fluid regime to the short-scale kinetic range of wavelengths in the development of the turbulent cascade of energy is nowadays subject of fervent discussion in the space plasmas scientific community. We make use of Hybrid Vlasov-Maxwell simulations where the full kinetic dynamics of ions is taken into account, while electrons are treated as a fluid. We investigate the development of turbulence in the solar wind, in 1D-3V phase space configuration and in the frequency range across the ion cyclotron frequency. These simulations allow for the analysis of the role of kinetic effects in the short-scale region of the energy spectra in the direction parallel to the background magnetic field. Our numerical results show the presence of a significant electrostatic activity at small wavelengths, triggered by the resonant interaction of ions with longitudinal waves. Our model does not allow to take into account the evolution of the turbulent spectra in the plane perpendicular to the ambient field, due to limited dimensionality in phase space. On the other hand, this model permits to isolate and study the possibility of transferring the electromagnetic large-scale energy on the small-scale kinetic electrostatic component of the spectrum. Peculiar features observed in the spacecraft data in the solar wind are qualitatively reproduced within the hybrid-Vlasov model, such as the generation of perpendicular temperature anisotropy and accelerated longitudinal beams of ions in the distribution of particle velocities as well as the appearance of a marked peak of electrostatic activity in the short-scale termination of the turbulent spectra.     

Electron transport and precipitation at Mercury during the MESSENGER flybys: Implications for electron-stimulated desorption

To examine electron transport, energization, and precipitation in Mercury's magnetosphere, a hybrid simulation study has been carried out that follows electron trajectories within the global magnetospheric electric and magnetic field configuration of Mercury. We report analysis for two solar-wind parameter conditions corresponding to the first two MESSENGER Mercury flybys on January 14, 2008, and October 6, 2008, which occurred for similar solar wind speed and density but contrasting interplanetary magnetic field (IMF) directions. During the first flyby the IMF had a northward component, while during the second flyby the IMF was southward. Electron trajectories are traced in the fields of global hybrid simulations for the two flybys. Some solar wind electrons follow complex trajectories at or near where dayside reconnection occurs and enter the magnetosphere at these locations. The entry locations depend on the IMF orientation (north or south). As the electrons move through the entry regions they can be energized as they execute non-adiabatic (demagnetized) motion. Some electrons become magnetically trapped and drift around the planet with energies on the order of 1–10 keV. The highest energy of electrons anywhere in the magnetosphere is about 25 keV, consistent with the absence of high-energy (>35 keV) electrons observed during either MESSENGER flyby. Once within the magnetosphere, a fraction of the electrons precipitates at the planetary surface with fluxes on the order of 10⁹ cm⁻² s⁻¹ and with energies of hundreds of eV. This finding has important implications for the viability of electron-stimulated desorption (ESD) as a mechanism for contributing to the formation of the exosphere and heavy ion cloud around Mercury. From laboratory estimates of ESD ion yields, a calculated ion production rate due to ESD at Mercury is found to be on par with ion sputtering yields.     

Radiation Induced Dynamics of the Galactic Halo

We show results of numerical simulations of a three component plasma consisting of electrons, ions and dust with external gravitation and radiation fields. We perform simulation runs, starting from an analytic halo equilibrium, balancing pressure, gravitational, and radiative forces. Within these the equilibrium is perturbed by the radiation of a typical OB-star association. The perturbation has a total energy input of 10⁷ L _⊙ and a duration of 30 Myrs. After switching off the perturbation, the simulations are continued to further investigate the dynamics induced. We start with a self consistent one-fluid MHD model without background magnetic field and show for an asymmetric case that the system approaches a new equilibrium after switching on the perturbation. Later it relaxes into the starting configuration again, when the additional radiation is turned off. We then show, first by including a disk-parallel magnetic field and then by redoing the simulations with a full three-fluid code, the influence of magnetic fields and species separation on the plasma dynamics. With our computations we demonstrate that these features can be important for the explanation of the structures of galactic halos and large scale mass flows. This revised version was published online in July 2006 with corrections to the Cover Date.     

First lunar wake passage of ARTEMIS: Discrimination of wake effects and solar wind fluctuations by 3D hybrid simulations

The spacecraft P1 of the new ARTEMIS (Acceleration, Reconnection, Turbulence, and Electrodynamics of the Moon's Interaction with the Sun) mission passed the lunar wake for the first time on February 13, 2010. We present magnetic field and plasma data of this event and results of 3D hybrid simulations. As the solar wind magnetic field was highly dynamic during the passage, a simulation with stationary solar wind input cannot distinguish whether distortions were caused by these solar wind variations or by the lunar wake; therefore, a dynamic real-time simulation of the flyby has been performed. The input values of this simulation are taken from NASA OMNI data and adapted to the P1 data, resulting in a good agreement between simulation and measurements. Combined with the stationary simulation showing non-transient lunar wake structures, a separation of solar wind and wake effects is achieved. An anisotropy in the magnitude of the plasma bulk flow velocity caused by a non-vanishing magnetic field component parallel to the solar wind flow and perturbations created by counterstreaming ions in the lunar wake are observed in data and simulations. The simulations help to interpret the data granting us the opportunity to examine the entire lunar plasma environment and, thus, extending the possibilities of measurements alone: A comparison of a simulation cross section to theoretical predictions of MHD wave propagation shows that all three basic MHD modes are present in the lunar wake and that their expansion governs the lunar wake refilling process.     

Magnetic field fossilization and tail reconfiguration in Titan's plasma environment during a magnetopause passage: 3D adaptive hybrid code simulations

We present a hybrid simulation study (kinetic ions, fluid electrons) of Titan's plasma interaction during an excursion of this moon from Saturn's magnetosphere into its magnetosheath, as observed for the first time during Cassini's T32 flyby on 13 June 2007. In contrast to earlier simulations of Titan's plasma environment under non-stationary upstream conditions, our model considers a difference in the flow directions of magnetospheric and magnetosheath plasma. Two complementary scenarios are investigated, with the flow directions of the impinging magnetospheric/magnetosheath plasmas being (A) antiparallel and (B) parallel. In both cases, our simulations show that due to the drastically reduced convection speed in the slow and dense heavy ion plasma near Titan, the satellite carries a bundle of “fossilized” magnetic field lines from the magnetosphere in the magnetosheath. Furthermore, the passage through Saturn's magnetopause goes along with a disruption of Titan's pick-up tail. Although the tail is not detached from the satellite, large clouds of heavy ion plasma are stripped of its outer flank, featuring a wave-like pattern. Whereas in case (B) under parallel flow conditions there is only a small retardation of about 5 min between the passage of Titan through the magnetopause and the reconfiguration of the pick-up tail, the tail reconfiguration in the case (A) scenario is completed not until 25 min after the magnetopause passage. The lifetime of fossil fields in the moon's ionosphere is approximately 25 min, regardless of whether parallel or antiparallel flow conditions are applied.     

Multi-fluid model of comet 1P/Halley

A new three-dimensional magnetohydrodynamic model of the coma of a comet has been developed and applied to simulations of a Halley-class coma using the solar-wind conditions of the Giotto flyby of Halley in 1986. The code developed for high-performance parallel processing computers, combines the high spatial resolution of smaller than 1 km grid spacing near the nucleus, with a large computational domain that enables structures nearly 10 million km down the comet tail to be modeled. Ions, neutrals, and electrons are considered as separate interacting fluids. Significant physical processes treated by the model include both photo and electron impact ionization of neutrals, recombination of ions, charge exchange between solar-wind ions and cometary neutrals, and frictional interactions between the three fluids considered in the model. A variety of plasma structures and physical parameters that are the output of this model are compared with relevant Giotto data from the 1986 Halley flyby.     

Origin of Mercury’s double magnetopause: 3D hybrid simulation study with A.I.K.E.F.

During the first and second Mercury flyby the MESSENGER spacecraft detected a dawn side double-current sheet inside the Hermean magnetosphere that was labeled the “double magnetopause” (Slavin, J.A. et al. [2008]. Science 321, 85). This double current sheet confines a region of decreased magnetic field that is referred to as Mercury’s “dayside boundary layer” (Anderson, M., Slavin, J., Horth, H. [2011]. Planet. Space Sci.). Up to the present day the double current sheet, the boundary layer and the key processes leading to their formation are not well understood. In order to advance the understanding of this region we have carried out self-consistent plasma simulations of the Hermean magnetosphere by means of the hybrid simulation code A.I.K.E.F. (Müller, J., Simon, S., Motschmann, U., Schüle, J., Glassmeier, K., Pringle, G.J. [2011]. Comput. Phys. Commun. 182, 946–966). Magnetic field and plasma results are in excellent agreement with the MESSENGER observations. In contrast to former speculations our results prove this double current sheet may exist in a pure solar wind hydrogen plasma, i.e. in the absence of any exospheric ions like sodium. Both currents are similar in orientation but the outer is stronger in intensity. While the outer current sheet can be considered the “classical” magnetopause, the inner current sheet between the magnetopause and Mercury’s surface reveals to be sustained by a diamagnetic current that originates from proton pressure gradients at Mercury’s inner magnetosphere. The pressure gradients in turn exist due to protons that are trapped on closed magnetic field lines and mirrored between north and south pole. Both, the dayside and nightside diamagnetic decreases that have been observed during the MESSENGER mission show to be direct consequences of this diamagnetic current that we label Mercury’s “boundary-layer-current“.     

Multicomponent theory of buoyancy instabilities in magnetized plasmas: the case of magnetic field parallel to gravity

We investigate electromagnetic buoyancy instabilities of the electron-ion plasma with the heat flux based on not the magnetohydrodynamic (MHD) equations, but using the multicomponent plasma approach when the momentum equations are solved for each species. We consider a geometry in which the background magnetic field, gravity, and stratification are directed along one axis. The nonzero background electron thermal flux is taken into account. Collisions between electrons and ions are included in the momentum equations. No simplifications usual for the one-fluid MHD-approach in studying these instabilities are used. We derive a simple dispersion relation, which shows that the thermal flux perturbation generally stabilizes an instability for the geometry under consideration. This result contradicts to conclusion obtained in the MHD-approach. We show that the reason of this contradiction is the simplified assumptions used in the MHD analysis of buoyancy instabilities and the role of the longitudinal electric field perturbation which is not captured by the ideal MHD equations. Our dispersion relation also shows that the medium with the electron thermal flux can be unstable, if the temperature gradients of ions and electrons have the opposite signs. The results obtained can be applied to the weakly collisional magnetized plasma objects in laboratory and astrophysics.     

Simulation studies of electron acceleration by ion ring distributions in solar flares

Using a 2 1/2-D fully relativistic electromagnetic particle-in-cell code (PIC) we have investigated a potential electron acceleration mechanism in solar flares. The free energy is provided by ions which have a ring velocity distribution about the magnetic field direction. Ion rings may be produced by perpendicular shocks, which could in turn be generated by the super-Alfvénic motion of magnetic flux tubes emerging from the photosphere or by coronal mass ejections (CMEs). Such ion distributions are known to be unstable to the generation of lower hybrid waves, which have phase velocities in excess of the electron thermal speed parallel to the field and can, therefore, resonantly accelerate electrons in that direction. The simulations show the transfer of perpendicular ion energy to energetic electrons via lower hybrid wave turbulence. With plausible ion ring velocities, the process can account for the observationally inferred fluxes and energies of non-thermal electrons during the impulsive phase of flares. Our results also show electrostatic wave generation close to the plasma frequency: we suggest that this is due to a bump-in-tail instability of the electron distribution.     

The equilibrium of atmospheric sodium

We now have four examples of planetary objects with detectable sodium (and potassium) in their atmospheres—Earth, Io, Mercury and the Moon. After a summary of the observational data, this survey discusses proposed sources and sinks. It appears that Io's surface material is rich in frozen SO₂, but with around 1% of some sodium compound. The Io plasma torus contains ions of S, O and Na, also with at least one molecular ion containing Na. In turn, impact by these ions probably sustains the torus, as well as an extended neutral corona. A primary source for the Earth, Mercury and the Moon is meteoroidal bombardment; at Mercury and perhaps the Moon it may be supplemented by degassing of atoms from the regolith. Photoionization is important everywhere, although hot electrons are dominant at Io.     

Ion cyclotron waves at Io: implications for the temporal variation of Io's atmosphere

When the flowing torus plasma encounters the upper atmosphere of Jupiter's moon, Io, newly created ions are rapidly accelerated by the motional electric field. Many of these ions are reneutralized and form a spray of fast neutrals that travel far away from Io before being reionized by photoionization and impact. These ions, now far from Io, are unstable to the generation of ion cyclotron waves. These waves in turn act as a mass spectrometer allowing Galileo magnetic measurements to be used to probe the composition of the atmosphere of Io and how it varies in time and in space. We now have six Galileo passes by Io on which we have measurements with sufficient cadence to examine the ion cyclotron waves. One of these passes, on Galileo's 32nd orbit has not been discussed previously. These passes provide sufficient observations to begin to distinguish the sources of variability. We find that while the atmosphere of Io varies temporally throughout the mission, it also has a spatial variation in composition at any instant of time.     

Observations of suprathermal electrons in Mercury's magnetosphere during the three MESSENGER flybys

In 2008 the MESSENGER spacecraft made the first direct observation of Mercury's magnetosphere in the more than 30 years since the Mariner 10 encounters. During MESSENGER's first flyby on 14 January 2008, the interplanetary magnetic field (IMF) was northward immediately prior to and following MESSENGER's equatorial passage through this small magnetosphere. The Energetic Particle Spectrometer (EPS), one of two sensors on the Energetic Particle and Plasma Spectrometer instrument that responds to electrons from ∼35 keV to 1 MeV and ions from ∼35 keV to 2.75 MeV, saw no increases in particle intensity above instrumental background (∼5 particles/cm²/sr/s/keV at 45 keV) at any time during the probe's magnetospheric passage. During MESSENGER's second flyby on 6 October 2008, there was a steady southward IMF, and intense reconnection was observed between the planet's magnetic field and the IMF. However, once again EPS did not observe bursts of energetic particles similar to those reported by Mariner 10 from its March 1974 encounter. On 29 September 2009, MESSENGER flew by Mercury for the third and final time before orbit insertion in March 2011. Although a spacecraft safe-hold event stopped science measurements prior to the outbound portion of the flyby, all instruments recorded full observations until a few minutes before the closest approach. In particular, the MESSENGER Magnetometer documented several substorm-like signatures of extreme loading of Mercury's magnetotail, but again EPS measured no energetic ions or electrons above instrument background during the inbound portion of the flyby. MESSENGER's X-Ray Spectrometer (XRS) nonetheless observed photons resulting from low-energy (∼10 keV) electrons impinging on its detectors during each of the three flybys. We infer that suprathermal plasma electrons below the EPS energy threshold caused the bremsstrahlung seen by XRS. In this paper, we summarize the energetic particle observations made by EPS and XRS during MESSENGER's three Mercury flybys, and we revisit the observations reported by Mariner 10 in the context of these new results.     

Monte Carlo modeling of Io’s [OI] 6300 Å and [SII] 6716 Å auroral emission in eclipse

We present a Monte Carlo (MC) model of [OI] 6300 Å and [SII] 6716 Å emission from Io entering eclipse. The simulation accounts for the 3-D distribution of SO₂, O, SO, S, and O₂ in Io’s atmosphere, several volcanic plumes, and the magnetic field around Io. Thermal electrons from the jovian plasma torus are input along the simulation domain boundaries and move along the magnetic field lines distorted by Io, occasionally participating in collisions with neutrals. We find that the atmospheric asymmetry resulting from varying degrees of atmospheric collapse across Io (due to eclipse ingress) and the presence of volcanoes contributes significantly to the unique morphology of the [OI] 6300 Å emission. The [OI] radiation lifetime of ∼134 s limits the emission to regions that have a sufficiently low neutral density so that intermolecular collisions are rare. We find that at low altitudes (typically <40 km) and in volcanic plumes (Pele, Prometheus, etc.) the number density is large enough (>4 × 10⁹ cm⁻³) to collisionally quench nearly all (>95%) of the excited oxygen for reasonable quenching efficiencies. Upstream (relative to the plasma flow), Io’s perturbation of the jovian magnetic field mirrors electrons with high pitch angles, while downstream collisions can trap the electrons. This magnetic field perturbation is one of the main physical mechanisms that results in the upstream/downstream brightness asymmetry in [OI] emission seen in the observation by Trauger et al. (Trauger, J.T., Stapelfeldt, K.R., Ballester, G.E., Clarke, J.I., 1997. HST observations of [OI] emissions from Io in eclipse. AAS-DPS Abstract (1997DPS29.1802T)). There are two other main causes for the observed brightness asymmetry. First, the observation’s viewing geometry of the wake spot crosses the dayside atmosphere and therefore the wake’s observational field of view includes higher oxygen column density than the upstream side. Second, the phased entry into eclipse results in less atmospheric collapse and thus higher collisional quenching on the upstream side relative to the wake. We compute a location (both in altitude and latitude) for the intense wake emission feature that agrees reasonably well with this observation. Furthermore, the peak intensity of the simulated wake feature is less than that observed by a factor of ∼3, most likely because our model does not include direct dissociation-excitation of SO₂ and SO. We find that the latitudinal location of the emission feature depends not so much on the tilt of the magnetic field as on the relative north/south flux tube depletion that occurs due to Io’s changing magnetic latitude in the plasma torus. From 1-D simulations, we also find that the intensity of [SII] 6716 and 6731 Å emission is much weaker than that of [OI] even if the [SII] excitation cross section is 10³ times larger than excitation to [OI]. This is because the density of S⁺ is much less than that of O and because the Einstein-A coefficient of the [SII] emission is a factor of ∼10 smaller than that of [OI].     

A Hall MHD model to study energetic charged particle events produced during fluctuations in the interplanetary magnetic field intensity

Energetic charged particles, which are often observed in solar active regions, may be also produced in interplanetary space due to the decoupling of ions and electrons in plasma. The Hall term in general Ohm's law is generally thought to be responsible for the decoupling of electrons and ions in plasma during magnetic reconnection. In this paper, a Hall MHD model is developed to study energetic charged particle events produced during fluctuations in the interplanetary magnetic field intensity. Two energetic charged particle events are used to test this model. It is concluded that the Hall effect does not only play the important role in the process of magnetic reconnection, but also in energetic charged particle events produced during fluctuations in the interplanetary magnetic field intensity.