Search

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.

We present a new hybrid kinetic model to simulate the response of plasmaspheric drainage plumes to impulsive interplanetary pressure pulses. Since particle distributions attending the interplanetary pulses and in the drainage plume are non-Maxwellian, wave-particle interactions play a crucial role in energy transport within and outside the plumes. Finite gyroradius effects become important in mass loading of the transmitted impulse with the drainage plume ions. A forward-reverse shock structure develops from the initial step-like transmitted shock. First results show that the impulse causes strong deformations in the global structure of the plume. The anisotropic ion velocity distribution functions at the impulse front and inside the plume help us determine energy transport via wave-particle interactions throughout the Earth’s inner magnetosphere.

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.

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.

We suggest a merging procedure for the Particle-Mesh Complex Particle Kinetic (PMCPK) method in case of inter-penetrating flow (multiple plasma beams). We examine the standard particle-in-cell (PIC) and the PMCPK methods in the case of particle acceleration by shock surfing for a wide range of the control numerical parameters. The plasma dynamics is described by a hybrid (particle-ion-fluid-electron) model. Note that one may need a mesh if modeling with the computation of an electromagnetic field. Our calculations use specified, time-independent electromagnetic fields for the shock, rather than self-consistently generated fields. While a particle-mesh method is a well-verified approach, the CPK method seems to be a good approach for multiscale modeling that includes multiple regions with various particle/fluid plasma behavior. However, the CPK method is still in need of a verification for studying the basic plasma phenomena: particle heating and acceleration by collisionless shocks, magnetic field reconnection, beam dynamics, etc.

A 2.5D numerical plasma model of the interaction of the solar wind (SW) with the Solar Probe Plus spacecraft (SPPSC) is presented. These results should be interpreted as a basic plasma model derived from the SW-interaction with the spacecraft (SC), which could have consequences for both plasma wave and electron plasma measurements on board the SC in the inner heliosphere. Compression waves and electric field jumps with amplitudes of about 1.5 V/m and (12-18) V/m were also observed. A strong polarization electric field was also observed in the wing of the plasma wake. However, 2.5D hybrid modeling did not show excitation of whistler/Alfven waves in the upstream connected with the bidirectional current closure that was observed in short-time 3D modeling SPPSC and near a tether in the ionosphere. The observed strong electromagnetic perturbations may be a crucial point in the electromagnetic measurements planned for the future Solar Probe Plus (SPP) mission. The results of modeling electromagnetic field perturbations in the SW due to shot noise in absence of SPPSC are also discussed.

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.

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 etal.,2007;Shematovichetal.,2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyro radius 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 O2 atmosphere is thermal with an extended cold population (Cassidyetal.,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.

A model is presented that describes the excitation of small-scale density perturbations and electromagnetic fields by standing shear Alfven waves in a Cartesian geometry. The model includes the effects of plasma betas and the Hall term effects. The characteristics of magnetospheric density cavities and the formation of the significant peak in density are discussed in the paper.

To be able to simulate the interaction of extrasolar planets with the stellar wind, a number of planetary parameters are required. Some of these (like planetary mass and radius) can be obtained directly from observational data. Other properties are not known very precisely. For example, up to now, there is no observation providing information on the strength of planetary magnetic moments. However, there is good reason to expect only very small magnetic moments for planets in very close orbits around their stars (like HD 209458 b and OGLE-TR-56 b). Thus, as a first step towards a more complete treatment, it seems reasonable to treat the interaction of the stellar wind with an unmagnetized planet. Calculations were performed for a nonconducting as well as for a weakly conducting planet. The interaction with the stellar wind and the resulting induced magnetosphere was simulated using a three dimensional hybrid code as well as in the drift-kinetic approximation. The effect of a interplanetary magnetic field oriented perpendicular to the incoming stellar wind was included. In the case of a weakly conducting body an asymmetrical Mach cone is formed, whereas for a nonconducting body no Mach cone is observed. These investigations will serve as the first step in the search for particular effects occurring at extrasolar planets, which could possibly lead to observable effects, e.g. radio emission. The results are also relevant for plasma structures near weakly conducting, unmagnetized bodies like the Earth's moon.

The interaction of the solar wind with weak comets, leading to the formation of the cometary magnetosphere with different types of structures in the solar wind and heavy ion plasmas, is simulated using a three-dimensional hybrid code. The simulations were performed for a wide range of gas production rates (6.2×10 26 mol / s ⩽Q⩽10 28 mol / s ) and for an interplanetary magnetic field which is perpendicular to the incoming solar wind. When Q 26 mol / s the cometary atmosphere forms a strong cycloid-type tail, whereas for Q⩾10 27 mol / s the cometary atmosphere forms a cone-type tail and structuring of the coma occurs. The results of these simulations may be applied to other weak massloading sources, e.g., dusty plasmas and cometary ion dynamics in the inner coma, AMPTE releases, and nonmagnetic bodies like Phobos, Deimos or even Pluto. Furthermore, the results presented here may be important for studying the dynamics of the ionized environment near the ‘Solar Probe’ spacecraft in the future.

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}

Bi-ion fluid simulations (Bogdanov et al., 1996; Sauer et al., 1995) have demonstrated several distinct features which occur in massloading zones associated with the interaction of the solar wind with comets, whose rate of gas production is small; specifically, bi-ion magneto-acoustic structuring and pulsations. Because of the limitated applicability of the bi-ion model, hybrid codes are used in order to check the results. The simulations were performed for gas production rates 5 × 1024s−1 ≤ Q ≤ 4 × 1027s−1 and for interplanetary magnetic field oriented perpendicular to simulation plane. In the case of Q < 1025s−1 the cometary atmosphere forms a strong cycloid-type tail. In case that Q ≥ 1025s−1, the cometary atmosphere forms a cone-type tail and structuring of the coma occurs. The results of these simulations apply to other weak solar wind massloading sources, e.g., the AMPTE releases, possibly Phobos and Deimos and possibly even Pluto. The results presented here may be important for studies of the dynamics of the ionized environment near the ‘Solar Probe’ spacecraft, in future.

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa moon-magnetosphere system with respect a 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; 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). In this paper we discuss two tasks: (1) the plasma wake structure dependence on the parameters of the upstream plasma and Europa's atmosphere (model I, cases (a) and (b) with a homogeneous Jovian magnetosphere field, an inductive magnetic dipole and high oceanic shell conductivity); and (2) estimation of the possible effect of an induced magnetic field arising from oceanic shell conductivity. This effect was estimated based on the difference between the observed and modeled magnetic fields (model II, case (c) with an inhomogeneous Jovian magnetosphere field, an inductive magnetic dipole and low oceanic shell conductivity)., 36 pages, 13 figures, paper was submitted to Planetary and Space Science in Nov. 2011, special volume "Outer Planets" but still under review

The paper reviews the studies of collective plasma processes responsible for anomalously fast exchange of energy between the plasma components, for a structure of the boundaries of the characteristic regions in the zone of the comet interaction with the solar wind, and for the anomalous transport across these boundaries. The position and structure of the cometary bow shock are obtained by numerical simulation using the particles in cells method. To study different plasma instabilities in the cometary atmospheres the simple analytic model of the cometary plasma environment is used that is based on the semikinetic description of the solar wind loading by cometary ions. The influence of a self-consistent electric field on the cometary plasma outflow is briefly discussed.

The paper considers the energy spectrum of particles accelerated near X-point of the magnetic field in the course of the tearing instability development in the current layer. The time dynamics of the spectrum at the linear and nonlinear stages of the instability is studied. The contribution of the particles inflowing into the vicinity of the magnetic X-point due to drift in the crossed electrical and magnetic fields was also taken into account. The dynamics of particle acceleration is numerically simulated and the results of the simulation are compared with the theoretical model. The analytical and numerical calculations are in reasonable agreement with the results of IMP satellites 7 and 8 energetic proton measurements in a wide energy interval.