Search

We study the dynamics of the thermal O+ and H+ ions at Ganymede's magnetopause when Ganymede is inside and outside of the Jovian plasma sheet using a three-dimensional hybrid model of plasma (kinetic ions, fluid electrons). We present the global structure of the electric fields and power density (E ⋅ J) in the magnetosphere of Ganymede and show that the power density at the magnetopause is mainly positive and on average is +0.95 and +0.75 nW/m3 when Ganymede is inside and outside the Jovian plasma sheet, respectively, but locally it reaches over +20 nW/m3. Our kinetic simulations show that ion velocity distributions at the vicinity of the upstream magnetopause of Ganymede are highly non-Maxwellian. We investigate the energization of the ions interacting with the magnetopause and find that the energy of those particles on average increases by a factor of 8 and 30 for the O+ and H+ ions, respectively. The energy of these ions is mostly within 1–100 keV for both species after interaction with the magnetopause, but a few percentages reach to 0.1–1 MeV. Our kinetic simulations show that a small fraction ((Formula presented.) 25%) of the corotating Jovian plasma reach the magnetopause, but among those >50% cross the high-power density regions at the magnetopause and gain energy. Finally, we compare our simulation results with Galileo observations of Ganymede's magnetopause crossings (i.e., G8 and G28 flybys). There is an excellent agreement between our simulations and observations, particularly our simulations fully capture the size and structure of the magnetosphere.

Context: The plasma near the nucleus of a comet is subjected to an electric field to which a few different sources contribute: the convective electric field of the solar wind, the ambipolar electric field due to higher electron than ion speeds, and a polarisation field arising from the vastly different ion and electron trajectories. Aims: Our aim is to show how the ambipolar and polarisation electric fields arise and develop under the influence of space charge effects, and in doing so we paint a qualitative picture of the electric fields in the inner coma of a comet. Methods. We use an electrostatic particle-in-cell model to simulate a scaled-down comet, representing comet 67P/Churyumov-Gerasimenko with parameters corresponding to a 3.0 AU heliocentric distance. Results: We find that an ambipolar electric field develops early in the simulation and that this is soon followed by the emergence of a polarisation electric field, manifesting itself as an anti-sunward component prevalent in the region surrounding the centre of the comet. As plasma is removed from the inner coma in the direction of the convectional electric field of the solar wind, a density maximum develops on the opposite side of the centre of the comet. Conclusions: The ambipolar and polarisation electric fields both have a significant influence on the motion of cometary ions. This demonstrates the importance of space charge effects in comet plasma physics.

The bow shock is the first boundary the solar wind encounters as it approaches planets or comets. The Rosetta spacecraft was able to observe the formation of a bow shock by following comet 67P/Churyumov-Gerasimenko toward the Sun, through perihelion, and back outward again. The spacecraft crossed the newly formed bow shock several times during two periods a few months before and after perihelion; it observed an increase in magnetic field magnitude and oscillation amplitude, electron and proton heating at the shock, and the diminution of the solar wind further downstream. Rosetta observed a cometary bow shock in its infancy, a stage in its development not previously accessible to in situ measurements at comets and planets.

Aims. We aim at analytically modelling the solar wind proton trajectories during their interaction with a partially ionised cometary atmosphere, not in terms of bulk properties of the flow but in terms of single particle dynamics. Methods. We first derive a generalised gyromotion, in which the electric field is reduced to its motional component. Steady-state is assumed, and simplified models of the cometary density and of the electron fluid are used to express the force experienced by individual solar wind protons during the interaction. Results. A three-dimensional (3D) analytical expression of the gyration of two interacting plasma beams is obtained. Applying it to a comet case, the force on protons is always perpendicular to their velocity and has an amplitude proportional to 1/r2. The solar wind deflection is obtained at any point in space. The resulting picture presents a caustic of intersecting trajectories, and a circular region is found that is completely free of particles. The particles do not lose any kinetic energy and this absence of deceleration, together with the solar wind deflection pattern and the presence of a solar wind ion cavity, is in good agreement with the general results of the Rosetta mission. Conclusions. The qualitative match between the model and the in situ data highlights how dominant the motional electric field is throughout most of the interaction region for the solar wind proton dynamics. The model provides a simple general kinetic description of how momentum is transferred between these two collisionless plasmas. It also shows the potential of this semi-analytical model for a systematic quantitative comparison to the data.

The presence or absence of a magnetic field determines the nature of how a planet interacts with the solar wind and what paths are available for atmospheric escape. Magnetospheres form both around magnetised planets, such as Earth, and unmagnetised planets, like Mars and Venus, but it has been suggested that magnetised planets are better protected against atmospheric loss. However, the observed mass escape rates from these three planets are similar (in the approximate (0.5–2) kg s−1 range), putting this latter hypothesis into question. Modelling the effects of a planetary magnetic field on the major atmospheric escape processes, we show that the escape rate can be higher for magnetised planets over a wide range of magnetisations due to escape of ions through the polar caps and cusps. Therefore, contrary to what has previously been believed, magnetisation is not a sufficient condition for protecting a planet from atmospheric loss. Estimates of the atmospheric escape rates from exoplanets must therefore address all escape processes and their dependence on the planet’s magnetisation.

Aims: The lack of an upstream solar wind plasma monitor when a spacecraft is inside the highly dynamic magnetosphere of Mercury limits interpretations of observed magnetospheric phenomena and their correlations with upstream solar wind variations. Methods: We used AMITIS, a three-dimensional GPU-based hybrid model of plasma (particle ions and fluid electrons) to infer the solar wind dynamic pressure and Alfvén Mach number upstream of Mercury by comparing our simulation results with MESSENGER magnetic field observations inside the magnetosphere of Mercury. We selected a few orbits of MESSENGER that have been analysed and compared with hybrid simulations before. Then we ran a number of simulations for each orbit (~30–50 runs) and examined the effects of the upstream solar wind plasma variations on the magnetic fields observed along the trajectory of MESSENGER to find the best agreement between our simulations and observations. Results: We show that, on average, the solar wind dynamic pressure for the selected orbits is slightly lower than the typical estimated dynamic pressure near the orbit of Mercury. However, we show that there is a good agreement between our hybrid simulation results and MESSENGER observations for our estimated solar wind parameters. We also compare the solar wind dynamic pressure inferred from our model with those predicted previously by the WSA-ENLIL model upstream of Mercury, and discuss the agreements and disagreements between the two model predictions. We show that the magnetosphere of Mercury is highly dynamic and controlled by the solar wind plasma and interplanetary magnetic field. In addition, in agreement with previous observations, our simulations show that there are quasi-trapped particles and a partial ring current-like structure in the nightside magnetosphere of Mercury, more evident during a northward interplanetary magnetic field (IMF). We also use our simulations to examine the correlation between the solar wind dynamic pressure a

Aims. We wish to investigate the energy conversion between particles and electromagnetic fields and determine the location where it occurs in the plasma environment of comets. Methods. We used a hybrid plasma model that included photoionization, and we considered two cases of the solar extreme ultraviolet flux. Other parameters corresponded to the conditions of comet 67P/Churyumov-Gerasimenko at a heliocentric distance of 1.5 AU. Results. We find that a shock-like structure is formed upstream of the comet and acts as an electromagnetic generator, similar to the bow shock at Earth that slows down the solar wind. The Poynting flux transports electromagnetic energy toward the inner coma, where newly born cometary ions are accelerated. Upstream of the shock-like structure, we find local energy transfer from solar wind ions to cometary ions. We show that mass loading can be a local process with a direct transfer of energy, but also part of a dynamo system with electromagnetic generators and loads. Conclusions. The energization of cometary ions is governed by a dynamo system for weak ionization, but changes into a large conversion region with local transfer of energy directly from solar wind protons for high ionization.

We use data from the first two dayside seasons of the Magnetospheric Multiscale (MMS) mission to study current systems associated with quasi‐perpendicular bow shocks of generator type. We have analyzed 154 MMS bow shock crossings near the equatorial plane. We compute the current density during the crossings and conclude that the component perpendicular to the shock normal (J⊥) is consistent with a pileup of the interplanetary magnetic field (IMF) inside the magnetosheath. For predominantly southward IMF, we observe a component Jn parallel (antiparallel) to the normal for GSM Y> 0 (<0), and oppositely directed for northward IMF. This indicates current closure across the equatorial magnetosheath, and it is observed for IMF clock angles near 0∘ and 180∘. To our knowledge, these are the first observational evidence for bow shock current closure across the magnetosheath. Since we observe no clear signatures of |J⊥| decreasing toward large |Y| we suggest that the main region of current closure is further tailward, outside MMS probing region. For IMF clock angles near 90∘, we find indications of the current system being tilted toward the north‐south direction, obtaining a significant Jz component, and we suggest that the current closes off the equatorial plane at higher latitudes where the spacecraft are not probing. The observations are complicated for several reasons. For example, variations in the solar wind and the magnetospheric currents and loads affect the closure, and Jn is distributed over large regions, making it difficult to resolve inside the magnetosheath proper.

Alfvén wings are known to form when a conducting or mass-loading object slows down a flowing plasma in its vicinity. Alfvén wings are not expected to be generated when an inert moon such as Rhea interacts with Saturn's magnetosphere, where the plasma impacting the moon is absorbed and the magnetic flux passes unimpeded through the moon. However, in two close polar passes of Rhea, Cassini clearly observed magnetic field signatures consistent with Alfvén wings. In addition, observations from a high-inclination flyby (Distance > 100 R Rh ) of Rhea on 3 June 2010 showed that the Alfvén wings continue to propagate away from Rhea even at this large distance. We have performed three-dimensional hybrid simulations of Rhea's interaction with Saturn's magnetosphere which show that the wake refilling process generates a plasma density gradient directed in the direction of corotating plasma. The resulting plasma pressure gradient exerts a force directed toward Rhea and slows down the plasma streaming into the wake along field lines. As on the same field lines, outside of the wake, the plasma continues to move close to its full speed, this differential motion of plasma bends the magnetic flux tubes, generating Alfvén wings in the wake. The current system excited by the Alfvén wings transfers momentum to the wake plasma extracting it from plasma outside the wake. Our work demonstrates that Alfvén wings can be excited even when a moon does not possess a conducting exosphere.

Context. The ESA/Rosetta mission made it possible to monitor the plasma environment of a comet, from near aphelion to perihelion conditions. To understand the complex dynamics and plasma structures found at the comet, a modelling effort must be carried out in parallel. Aims. Firstly, we present a 3D hybrid model of the cometary plasma environment including photoionisation, solar wind charge exchange, and electron ionisation reactions; this model is used in stationary and dynamic conditions (mimicking the solar wind variations), and is thus especially adapted to a weakly outgassing comet such as 67P/Churyumov-Gerasimenko, the target of the ESA/Rosetta mission. Secondly, we use the model to study the respective effects of ionisation processes on the formation of the dayside macroscopic magnetic and density boundaries upstream of comet 67P in perihelion conditions at 1.3 AU. Thirdly, we explore and discuss the effects of these processes on the magnetic field line draping, ionisation rates, and composition in the context of the Rosetta mission. Methods. We used a new quasi-neutral hybrid model, originally designed for weakly magnetised planetary bodies, such as Venus, Mars, and Titan, and adapted here to comets. Ionisation processes were monitored individually and together following a probabilistic interaction scheme. Three-dimensional paraboloid fits of the bow shock surface, identified for a magnetosonic Mach number equal to 2, and of the cometopause surface, were performed for a more quantitative analysis. Results. We show that charge exchange and electron ionisation play a major role in the formation of a bow shock-like structure far upstream, while photoionisation is the main driver at and below the cometopause boundary, within 1000 km cometocentric distance. Charge exchange contributes to 42% of the total production rate in the simulation box, whereas production rates from electron ionisation and photoionisation reach 33% and 25%, respectively. We also discuss imp

The solar wind interaction with Ceres is studied for a high water vapor release from its surface using a hybrid model including photoionization. We use a water vapor production rate of 6 kg/s, thought to be due to subsurface sublimation, corresponding to a detection on 6 March 2013 by the Herschel Space Observatory. We present the general morphology of the plasma interactions, both close to Ceres and on a larger scale. Mass loading of water ions causes a magnetic pileup region in front of Ceres, where the solar wind deflects up to 15 ∘ and slows down by 15%. The global plasma interaction with Ceres is not greatly affected by the source location of water vapor nor on gravity, only on the production rate of water vapor. On a global scale, Ceres has a comet-like interaction with the solar wind with observable perturbations farther than 250 Ceres radii downstream of the body.

Here I study the “plasma interactions with icy bodies in the solar system”, that is, my quest to understand the fundamental processes that govern such interactions. By using numerical modelling combined with in situ observations, one can infer the internal structure of icy bodies and their plasma environments. After a broad overview of the laws governing space plasmas a more detailed part follows. This contains the method on how to model the interaction between space plasmas and icy bodies. Numerical modelling of space plasmas is applied to the icy bodies Callisto (a satellite of Jupiter), the dwarf planet Ceres (located in the asteroid main belt) and the comet 67P/Churyumov-Gerasimenko. The time-varying magnetic field of Jupiter induces currents inside the electrically conducting moon Callisto. These create magnetic field perturbations thought to be related to conducting subsurface oceans. The flow of plasma in the vicinity of Callisto is greatly affected by these magnetic field perturbations. By using a hybrid plasma solver, the interaction has been modelled when including magnetic induction and agrees well with magnetometer data from flybys (C3 and C9) made by the Galileo spacecraft. The magnetic field configuration allows an inflow of ions onto Callisto’s surface in the central wake. Plasma that hits the surface knocks away matter (sputtering) and creates Callisto’s tenuous atmosphere. A long term study of solar wind protons as seen by the Rosetta spacecraft was conducted as the comet 67P/Churyumov-Gerasimenko approached the Sun. Here, extreme ultraviolet radiation from the Sun ionizes the neutral water of the comet’s coma. Newly produced water ions get picked up by the solar wind flow, and forces the solar wind protons to deflect due to conservation of momentum. This effect of mass-loading increases steadily as the comet draws closer to the Sun. The solar wind is deflected, but does not lose much energy. Hybrid modelling of the solar wind interaction with the c, Här studerar jag “plasmaväxelverkan med isiga kroppar i solsystemet”, det vill säga, min strävan är att förstå de grundläggande processerna som styr sådana interaktioner. Genom att använda numerisk modellering i kombination med observationer på plats vid himlakropparna kan man förstå sig på deras interna strukturer och rymdmiljöer. Efter en bred översikt över de fysiska lagar som styr ett rymdplasma följer en mer detaljerad del. Denna innehåller metoder för hur man kan modellera växelverkan mellan rymdplasma och isiga kroppar. Numerisk modellering av rymdplasma appliceras på de isiga himlakropparna Callisto (en måne kring Jupiter), dvärgplaneten Ceres (lokaliserad i asteroidbältet mellan Mars och Jupiter) och kometen 67P/Churyumov-Gerasimenko. Det tidsvarierande magnetiska fältet kring Jupiter inducerar strömmar inuti den elektriskt ledande månen Callisto. Dessa strömmar skapar magnetfältsstörningar som tros vara relaterade till ett elektriskt ledande hav under Callistos yta. Plasmaflödet i närheten av Callisto påverkas i hög grad av dessa magnetfältsstörningar. Genom att använda en hybrid-plasma-lösare har växelverkan modellerats, där effekten av magnetisk induktion har inkluderats. Resultaten stämmer väl överens med magnetfältsdata från förbiflygningarna av Callisto (C3 och C9) som gjordes av den obemannade rymdfarkosten Galileo i dess bana kring Jupiter. Den magnetiska konfigurationen som uppstår möjliggör ett inflöde av laddade joner på Callistos baksida. Plasma som träffar ytan slår bort materia och skapar Callistos tunna atmosfär. En långtidsstudie av solvindsprotoner sett från rymdfarkosten Rosetta utfördes då kometen 67P/Churyumov-Gerasimenko närmade sig solen. Ultraviolett strålning från solen joniserar det neutrala vattnet i kometens koma (kometens atmosfär). Nyligt joniserade vattenmolekyler plockas upp av solvindsflödet och tvingar solvindsprotonernas banor att böjas av, så att rörelsemängden bevaras. Denna effekt ökar stadigt då kometen närmar sig solen

Context. The first long-term in-situ observation of the plasma environment in the vicinity of a comet, as provided by the European Rosetta spacecraft. Aims. Here we offer characterisation of the solar wind flow near 67P/Churyumov-Gerasimenko (67P) and its long term evolution during low nucleus activity. We also aim to quantify and interpret the deflection and deceleration of the flow expected from ionization of neutral cometary particles within the undisturbed solar wind. Methods. We have analysed in situ ion and magnetic field data and combined this with hybrid modeling of the interaction between the solar wind and the comet atmosphere. Results. The solar wind deflection is increasing with decreasing heliocentric distances, and exhibits very little deceleration. This is seen both in observations and in modeled solar wind protons. According to our model, energy and momentum are transferred from the solar wind to the coma in a single region, centered on the nucleus, with a size in the order of 1000 km. This interaction affects, over larger scales, the downstream modeled solar wind flow. The energy gained by the cometary ions is a small fraction of the energy available in the solar wind. Conclusions. The deflection of the solar wind is the strongest and clearest signature of the mass-loading for a small, low-activity comet, whereas there is little deceleration of the solar wind.

By using a hybrid plasma solver (ions as particles and electrons as a fluid), we have modeled the interaction between Callisto and Jupiter's magnetosphere for variable ambient plasma parameters. We compared the results with the magnetometer data from flybys (C3, C9, and C10) by the Galileo spacecraft. Modeling the interaction between Callisto and Jupiter's magnetosphere is important to establish the origin of the magnetic field perturbations observed by Galileo and thought to be related to a subsurface ocean. Using typical upstream magnetospheric plasma parameters and a magnetic dipole corresponding to the inductive response inside the moon, we show that the model results agree well with observations for the C3 and C9 flybys, but agrees poorly with the C10 flyby close to Callisto. The study does support the existence of a subsurface ocean at Callisto.