Context . We use a global 3D hybrid plasma model to investigate the interaction between Mercury’s magnetosphere and the solar wind for the second BepiColombo swingby, evaluate magnetospheric regions, and study the typical energy profile of protons. Aims . The objective of this study is to gain a better understanding of solar wind entry and analyze simulated plasma data along a trajectory using BepiColombo swingby 2 conditions, with the goal of enhancing our comprehension of measurement data and potentially providing forecasts for future swingbys. Methods . To model Mercury’s plasma environment, we used the hybrid code AIKEF and developed a method to extract the particle (ion) data in order to compute the proton energy spectrum along the trajectory of BepiColombo during its second Mercury swingby on June 23, 2022. We evaluate magnetopause and bow shock stand-off distances under average upstream solar wind conditions with the Interplanetary Magnetic Field (IMF) condition derived from the BepiColombo magnetic field measurements during the second Mercury swingby. Results . We found that the magnetosheath on the quasi-perpendicular (dusk) side of the bow shock is thicker than that on the quasi-parallel (dawn) side, where a foreshock is formed. Multiple plasma populations can be extracted from our modeled energy spectra that assist in identifying magnetospheric regions. We observed protons of solar wind origin entering Mercury’s magnetosphere. Their energies range from a few electron volts in the magnetosphere up to 10 keV in the magnetosheath.

The induced magnetotails on Mars and Venus are considered to arise through the interplanetary magnetic field (IMF) draping around the planet and the solar wind deceleration due to the mass loading effect. They have very similar structures as that on Earth, two magnetic lobes of opposite radial magnetic fields and a plasma sheet in between. However, the orientation and geometry of the induced magnetotails are controlled by the IMF, not the planetary intrinsic magnetic field. In this study, we present another characteristic of the induced magnetotails on Mars and Venus with the observations of MAVEN and Venus Express. It is found that the magnetic flux in the induced magnetotails on Mars and Venus are inhomogeneous. There is more magnetic flux in the +E hemisphere than -E hemisphere. The magnetic flux is observed to transport gradually from the +E hemisphere to the -E hemisphere along the magnetotail. The magnetotail magnetic flux transport seems to be faster on Mars than that at Venus. Based on these observations, we suggest that the finite gyro-radius effect of the planetary ions that are picked up by the solar wind is responsible to the magnetic flux inhomogeneity and transport in the induced magnetotails. The role of the magnetic pressure gradient in the magnetotail will be discussed.

Magnetic fields inconsistent with draped interplanetary magnetic fields and crustal fields have been observed on Mars. Considering the discovery of a global looping magnetic field around the Venusian magnetotail and the similarities in the solar wind interactions between Mars and Venus, we use Mars Atmosphere and Volatile Evolution observations to investigate the global looping field on Mars and its formation mechanism. It is found that a global looping field also exists on Mars; therefore, this type of global looping field is a common feature of unmagnetized planetary bodies with ionospheres, and therefore should also exist on Titan and near- Sun comets. The comparison of the looping fields on Mars and Venus shows that the looping field is stronger on Mars. Solar wind azimuthal flows around the magnetotail toward the -E magnetotail polar region (X-MSE < 0, Y-MSE = 0, Z(MSE) < - 1R(M)) are observed. We illustrate that the looping field can be formed by bending the draped field lines with these azimuthal flows, and that these azimuthal flows are associated with heavy ion plumes along the +E direction that are expected to be stronger on Mars than Venus. The current system associated with the looping field and its possible connection with the nightside ionosphere formations and ion escapes on Mars and Venus are discussed.

We model the emission of energetic neutral atoms (ENAs) that are generated by the interaction between energetic ions from Saturn's magnetosphere and neutrals from the upper atmosphere of the giant planet's largest moon Titan. The trajectories of the parent ions and the resulting ENA emission morphology are highly sensitive to the electromagnetic field configuration near the moon. We therefore compare the ENA emission pattern for spatially homogeneous fields to the emission obtained from a magnetohydrodynamic (MHD) and a hybrid (kinetic ions, fluid electrons) model of Titan's magnetospheric interaction, by computing the trajectories of several billion energetic test particles. While the MHD model takes into account the draping of the magnetic field lines around Titan, the hybrid approach also considers the significant asymmetries in the electromagnetic fields due to the large gyroradii of pick-up ions from Titan's ionosphere. In all three models, the upstream parameters correspond to the conditions during Cassini's TA flyby of Titan. The shape, magnitude, and location of the ENA emission maxima vary considerably between these three field configurations. The magnetic pile-up region at Titan's ramside deflects a large number of the energetic parent ions, thereby reducing the ENA flux. However, the draped magnetic field lines in Titan's lobes rotate the gyration planes of the incident energetic ions, thereby facilitating the observable ENA production. Overall, the ENA flux calculated for the MHD model is weaker than the emission obtained for the electromagnetic fields from the hybrid code. In addition, we systematically investigate the dependency of the ENA emission morphology on the energy of the parent ions and on the upstream magnetic field strength.

Using hybrid simulations and analytical calculations, we investigate the observable magnetic perturbations during the 12 planned Callisto flybys of the JUpiter ICy moons Explorer mission. During four of these encounters, Callisto will be embedded within Jupiter's magnetospheric current sheet. In these cases, Callisto's Alfven wings and ramside magnetic field pileup will partially obscure any magnetic signatures associated with induction in a conducting layer at the moon, thereby severely complicating attempts to further constrain properties of a possible subsurface ocean. During one of these flybys, the plasma interaction will even generate magnetic signatures that are qualitatively similar to an induced field from the moon's interior. In this case, highly accurate measurements of the upstream flow parameters and Callisto's ionosphere are required to disentangle the induction signal from plasma effects. For the remaining eight encounters, Callisto's plasma interaction is expected to be sufficiently weak for an unobstructed observation of the moon's inductive signature.

We apply a hybrid (kinetic ions and fluid electrons) simulation model to study Pluto's plasma environment during the New Horizons encounter on 14 July 2015. We show that Pluto's plasma interaction is dominated by significant north-south asymmetries, driven by large pickup ion gyroradii on the order of 200 Pluto radii. The transition region from the ambient solar wind to the population of plutogenic ions (called the "Plutopause") also shows considerable asymmetries that cannot be explained by a fluid picture. Since the New Horizons spacecraft does not carry a magnetometer, we use our model to estimate the strength and direction of the interplanetary magnetic field (IMF) at the time of the flyby by comparing output from the hybrid simulation to the plasma signatures observed during the New Horizons encounter. We find that an IMF strength of at least 0.24 nT is required to generate the observed plasma signatures. An IMF orientation either parallel or antiparallel to Pluto's orbital motion is able to explain the observed plasma densities and velocities along the New Horizons trajectory. Our simulations are able to quantitatively reproduce all key features of the plasma observations, specifically the gradual slowing of the solar wind, as well as the location and thickness of the Plutopause and bow shock.

We apply a combination of analytical modeling, hybrid simulations, and data analysis techniques to provide a comprehensive study of magnetometer data from four Galileo flybys of Callisto (C21, C22, C23, and C30) that have never been discussed in the literature before. Callisto's distance to the center of Jupiter's magnetospheric current sheet varied considerably from flyby to flyby. Therefore, the relative strength of the magnetic field perturbations due to Callisto's plasma interaction with Jupiter's magnetosphere and induction within Callisto's subsurface ocean drastically changed as well. During C21, a strong magnetic field perturbation along the corotation direction was detected in Callisto's geometric plasma shadow. This enhancement can be explained with Callisto's steady state plasma interaction only, if the upstream flow possessed a nonnegligible component away from Jupiter. During C22, Galileo only grazed Callisto's Alfvén wings which were elevated out of the flyby plane due to the ambient magnetospheric field orientation. During C23, the combination of an inclined flyby trajectory and finite gyroradius effects caused Callisto's observed Alfvén wings to be slightly asymmetric between both hemispheres. During C30, a discontinuity with a surface normal pointed toward Jupiter was detected within Callisto's geometric plasma shadow, similar to the earlier C10 flyby. Due to strong plasma interaction and an unfavorable flyby geometry (C21), a large closest approach altitude (C22), or weak inducing field (C23 and C30), no discernible induction signatures were observed during these four flybys. Based on data from all available Galileo flybys, we determine requirements on future flyby geometries that must be satisfied for an identification of Callisto's subsurface ocean in magnetometer data.

We apply a combination of data analysis and hybrid modeling to study Callisto's interaction with Jupiter's magnetosphere during the Galileo C10 flyby on 17 September 1997. This encounter took place while Callisto was located near the center of Jupiter's current sheet. Therefore, induction in Callisto's subsurface ocean and magnetospheric field line draping around the moon's ionosphere both made nonnegligible contributions to the observed magnetic perturbations. The induction signal during C10 was obscured by plasma currents to a significant degree, in contrast to previously studied Callisto flybys. Our analysis reveals that at large distances to Callisto, its magnetic environment was dominated by field line draping, leading to the formation of Alfven wings. Closer to the surface and in Callisto's wake, Galileo encountered a quasi-dipolar core region that was partially shielded from the plasma interaction and was dominated by the induced field. When exiting this core region, the spacecraft crossed a rotational discontinuity where the magnetic field vector rotated by approximately 50 degrees. The hybrid model is able to quantitatively explain numerous key features of the observed magnetic signatures, especially the transitions between draping- and dipole-dominated regimes along the C10 trajectory. The model also reproduces the electron number density enhancement by 3-4 orders of magnitude detected in Callisto's wake, requiring a substantial ionosphere to surround the moon during C10. For flybys with nonnegligible plasma currents, comprehensive knowledge of the incident flow conditions and properties of Callisto's atmosphere is required to refine existing constraints on the subsurface ocean (conductivity, thickness, and depth) based on magnetic field data. These findings are highly relevant for the upcoming JUpiter ICy moon Explorer (JUICE) mission, which will include multiple Callisto flybys.

Abstract By applying a hybrid (kinetic ions and fluid electrons) simulation code, we study the plasma environment of Saturn's largest moon Titan during Cassini's T96 flyby on 1 December 2013. The T96 encounter marks the only observed event of the entire Cassini mission where Titan was located in the supersonic solar wind in front of Saturn's bow shock. Our simulations can quantitatively reproduce the key features of Cassini magnetic field and electron density observations during this encounter. We demonstrate that the large‐scale features of Titan's induced magnetosphere during T96 can be described in terms of a steady state interaction with a high‐pressure solar wind flow. About 40 min before the encounter, Cassini observed a rotation of the incident solar wind magnetic field by almost 90°. We provide strong evidence that this rotation left a bundle of fossilized magnetic field lines in Titan's ionosphere that was subsequently detected by the spacecraft.

We study how the local electromagnetic disturbances introduced by Titan affect the ionization rates of the atmosphere. For this, we model the precipitation of energetic particles, specifically hydrogen and oxygen ions with energies between 1keV and 1MeV, into Titan׳s exobase for the specific magnetospheric configuration of the T9 flyby. For the study, a particle tracing software package is used which consists of an integration of the single particle Lorentz force equation using a 4th order Runge–Kutta numerical method. For the electromagnetic disturbances, the output of the A.I.K.E.F. hybrid code (kinetic ions, fluid electrons) is used, allowing the possibility of analyzing the disturbances and asymmetries in the access of energetic particles originated by their large gyroradii. By combining these methods, 2D maps showing the access of each set of particles were produced. We show that the access of different particles is largely dominated by their gyroradii, with the complexity of the maps increasing with decreasing gyroradius, due to the larger effect that local disturbances introduced by the presence of the moon have in the trajectory of the particles with lower energies. We also show that for particles with gyroradii much larger than the moon׳s radius, simpler descriptions of the electromagnetic environment can reproduce similar results to those obtained when using the full hybrid simulation description, with simple north–south fields being sufficient to reproduce the hybrid code results for O+ ions with energies larger than 10keV but not enough to reproduce those for H+ ions at any of the energies covered in the present study. Finally, by combining the maps created with upstream plasma flow measurements by the MIMI/CHEMS instrument, we are able to estimate normalized fluxes arriving at different selected positions of the moon׳s exobase. We then use these fluxes to calculate energy deposition and non-dissociative N2 ionization rates for precipitating O+ and H+ ions and find differences in the ion production rates of up to almost 80% at the selected positions. All these results combined show that the electromagnetic field disturbances present in the vicinity of Titan significantly affect the contribution of energetic ions to local ionization profiles.

The interaction between Callisto's atmosphere and ionosphere with the surrounding magnetospheric environment is analyzed by applying a hybrid simulation code, in which the ions are treated as particles and the electrons are treated as a fluid. Callisto is unique among the Galilean satellites in its interaction with the ambient magnetospheric plasma as the gyroradii of the impinging plasma and pickup ions are large compared to the size of the moon. A kinetic representation of the ions is therefore mandatory to adequately describe the resulting asymmetries in the electromagnetic fields and the deflection of the plasma flow near Callisto. Multiple model runs are performed at various distances of the moon to the center of Jupiter's magnetospheric current sheet, with differing angles between the corotational plasma flow and the ionizing solar radiation. When Callisto is embedded in the Jovian current sheet, magnetic perturbations due to the plasma interaction are more than twice the strength of the background field and may therefore obscure any magnetic signal generated via induction in a subsurface ocean. The magnetic field perturbations generated by Callisto's ionospheric interaction are very similar at different orbital positions of the moon, demonstrating that local time is only of minor importance when disentangling magnetic signals generated by the magnetosphere‐ionosphere interaction from those driven by induction. Our simulations also suggest that deflection of the magnetospheric plasma around the moon cannot alone explain the density enhancement of 2 orders of magnitude measured in Callisto's wake during Galileo flybys. However, through inclusion of an ionosphere surrounding Callisto, modeled densities in the wake are consistent with in situ measurements.