Umeå University's logo

We examine the origin of electrons in a weakly outgassing comet, using Rosetta mission data and a 3D collisional model of electrons at a comet. We have calculated a new data set of electron-impact ionization (EII) frequency throughout the Rosetta escort phase, with measurements of the Rosetta Plasma Consortium's Ion and Electron Sensor (RPC/IES). The EII frequency is evaluated in 15-min intervals and compared to other Rosetta data sets. EII is the dominant source of electrons at 67P away from perihelion and is highly variable (by up to three orders of magnitude). Around perihelion, EII is much less variable and less efficient than photoionization at Rosetta. Several drivers of the EII frequency are identified, including magnetic field strength and the outgassing rate. Energetic electrons are correlated to the Rosetta-upstream solar wind potential difference, confirming that the ionizing electrons are solar wind electrons accelerated by an ambipolar field. The collisional test particle model incorporates a spherically symmetric, pure water coma and all the relevant electron-neutral collision processes. Electric and magnetic fields are stationary model inputs, and are computed using a fully kinetic, collision-less Particle-in-Cell simulation. Collisional electrons are modelled at outgassing rates of Q = 1026 s-1 and Q = 1.5 × 1027 s-1. Secondary electrons are the dominant population within a weakly outgassing comet. These are produced by collisions of solar wind electrons with the neutral coma. The implications of large ion flow speed estimates at Rosetta, away from perihelion, are discussed in relation to multi-instrument studies and the new results of the EII frequency obtained in this study.

We present fast algorithms to automatically estimate the statistical position of the bow shock from spacecraft data, using existing analytical two-dimensional (2D) and three-dimensional (3D) models of the shock surface. We derive expressions of the standoff distances in 2D and 3D and of the normal to the bow shock at any given point on it. Two simple bow shock detection algorithms are constructed, one solely based on a geometrical predictor from existing models, the other using this predicted position to further refine it with the help of magnetometer data, an instrument flown on many planetary missions. Both empirical techniques are applicable to any planetary environment with a defined shock structure. Applied to the Martian environment and the NASA/MAVEN mission, the predicted shock position is on average within 0.15 planetary radius Rp of the bow shock crossing. Using the predictor-corrector algorithm, this estimate is further refined to within a few minutes of the true crossing (≈0.05Rp). Between 2014 and 2021, we detect 14,929 clear bow shock crossings, predominantly quasi-perpendicular. Thanks to 2D conic and 3D quadratic fits, we investigate the variability of the shock surface with respect to Mars Years (MY), solar longitude (Ls), and solar EUV flux levels. Although asymmetry in Y and Z Mars Solar Orbital coordinates is on average small, we show that for MY32 and MY35, Ls = [135°−225°] and high solar flux, it can become particularly noticeable, and is superimposed to the usual North-South asymmetry due in part to the presence of crustal magnetic fields.

Context. The observation of pits at the surface of comets offers the opportunity to take a glimpse into the properties and the mechanisms that shape a nucleus through cometary activity. If the origin of these pits is still a matter of debate, multiple studies have recently suggested that known phase transitions (such as volatile sublimation or amorphous water ice crystallization) alone could not have carved these morphological features on the surface of 67P/Churyumov-Gerasimenko (hereafter 67P).

Aims. We want to understand how the progressive modification of 67P' s surface due to cometary activity might have affected the characteristics of pits and alcoves. In particular, we aim to understand whether signatures of the formation mechanism of these surface morphological features can still be identified.

Methods. To quantify the amount of erosion sustained at the surface of 67P since it arrived on its currently observed orbit, we selected 380 facets of a medium-resolution shape model of the nucleus, sampling 30 pits and alcoves across the surface. We computed the surface energy balance with a high temporal resolution, including shadowing and self-heating contributions. We then applied a thermal evolution model to assess the amount of erosion sustained after ten orbital revolutions under current illumination conditions.

Results. We find that the maximum erosion sustained after ten orbital revolutions is on the order of 80 m, for facets located in the southern hemisphere. We thus confirm that progressive erosion cannot form pits and alcoves, as local erosion is much lower than their observed depth and diameter. We find that plateaus tend to erode more than bottoms, especially for the deepest depressions, and that some differential erosion can affect their morphology. As a general rule, our results suggest that sharp morphological features tend to be erased by progressive erosion.

Conclusions. This study supports the assumption that deep circular pits, such as Seth_01, are the least processed morphological features at the surface of 67P, or the best preserved since their formation.

Context: Amongst the different features and boundaries encountered around comets, one remains of particular interest to the plasma community: the diamagnetic cavity. Crossed for the first time at 1P/Halley during the Giotto flyby in 1986 and later met more than 700 times by the ESA Rosetta spacecraft around Comet 67P/Churyumov-Gerasimenko, this region, almost free of any magnetic field, surrounds nuclei of active comets. However, previous observations and modelling of this part of the coma have not yet provided a definitive answer as to the origin of such a cavity and on its border, the diamagnetic cavity boundary layer.

Aims: We investigate which forces and equilibrium might be at play and balance the magnetic pressure at this boundary down to the spatial and temporal scales of the electrons in the 1D collisionless case. In addition, we scrutinise assumptions made in magneto-hydrodynamic and hybrid simulations of this environment and check for their validity.

Methods: We simulated this region at the electron scale by means of 1D3V particle-in-cell simulations and SMILEI code.

Results: Across this layer, depending on the magnetic field strength, the electric field is governed by different equilibria, with a thin double-layer forming ahead. In addition, we show that the electron distribution function departs from Maxwellian and/or gyrotropic distributions and that electrons do not behave adiabatically. We demonstrate the need to investigate this region at the electron scale in depth with fully kinetic simulations.

In the White Paper, submitted in response to the European Space Agency (ESA) Voyage 2050 Call, we present the importance of advancing our knowledge of plasma-neutral gas interactions, and of deepening our understanding of the partially ionized environments that are ubiquitous in the upper atmospheres of planets and moons, and elsewhere in space. In future space missions, the above task requires addressing the following fundamental questions: (A) How and by how much do plasma-neutral gas interactions influence the re-distribution of externally provided energy to the composing species? (B) How and by how much do plasma-neutral gas interactions contribute toward the growth of heavy complex molecules and biomolecules? Answering these questions is an absolute prerequisite for addressing the long-standing questions of atmospheric escape, the origin of biomolecules, and their role in the evolution of planets, moons, or comets, under the influence of energy sources in the form of electromagnetic and corpuscular radiation, because low-energy ion-neutral cross-sections in space cannot be reproduced quantitatively in laboratories for conditions of satisfying, particularly, (1) low-temperatures, (2) tenuous or strong gradients or layered media, and (3) in low-gravity plasma. Measurements with a minimum core instrument package (< 15 kg) can be used to perform such investigations in many different conditions and should be included in all deep-space missions. These investigations, if specific ranges of background parameters are considered, can also be pursued for Earth, Mars, and Venus.

The environment of a comet is a fascinating and unique laboratory to study plasma processes and the formation of structures such as shocks and discontinuities from electron scales to ion scales and above. The European Space Agency's Rosetta mission collected data for more than two years, from the rendezvous with comet 67P/Churyumov-Gerasimenko in August 2014 until the final touch-down of the spacecraft end of September 2016. This escort phase spanned a large arc of the comet's orbit around the Sun, including its perihelion and corresponding to heliocentric distances between 3.8 AU and 1.24 AU. The length of the active mission together with this span in heliocentric and cometocentric distances make the Rosetta data set unique and much richer than sets obtained with previous cometary probes. Here, we review the results from the Rosetta mission that pertain to the plasma environment. We detail all known sources and losses of the plasma and typical processes within it. The findings from in-situ plasma measurements are complemented by remote observations of emissions from the plasma. Overviews of the methods and instruments used in the study are given as well as a short review of the Rosetta mission. The long duration of the Rosetta mission provides the opportunity to better understand how the importance of these processes changes depending on parameters like the outgassing rate and the solar wind conditions. We discuss how the shape and existence of large scale structures depend on these parameters and how the plasma within different regions of the plasma environment can be characterised. We end with a non-exhaustive list of still open questions, as well as suggestions on how to answer them in the future.

The flow direction of low-energy ions around comet 67P/Churyumov–Gerasimenko has previously been difficult to constrain due to the influence of the spacecraft potential. The Ion Composition Analyzer of the Rosetta Plasma Consortium (RPC-ICA) on Rosetta measured the distribution function of positive ions with energies down to just a few eV/q throughout the escort phase ofthe mission. Unfortunately, the substantial negative spacecraft potential distorted the directional information of the low-energy data. In this work, we present the flow directions of low-energy ions around comet 67P, corrected for the spacecraft potential using Particle-In-Cell simulation results. We focus on the region in and around the diamagnetic cavity, where low-energy ions are especially important for the dynamics. We separate between slightly accelerated ‘burst’ features and a more constant ‘band’ of low-energy ions visible in the data. The ‘bursts’ are flowing radially outwards from the nucleus with an antisunward component while the ‘band’ is predominantly streaming back towards the comet. This provides evidence of counter-streaming ions, which has implications for the overall expansion velocity of the ions. The backstreaming ions are present also at times when the diamagnetic cavity was not detected, indicating that the process accelerating the ions back towards the comet is not connected to the cavity boundary.

Ion acoustic waves were observed between 15 and 30 km from the centre of comet 67P/Churyumov–Gerasimenko by the Rosetta spacecraft during its close flyby on 28 March 2015. There are two electron populations: one cold at kBTe≈0.2 eV and one warm at kBTe≈2 eV. The ions are dominated by a cold (a few hundredths of electronvolt) distribution of water group ions with a bulk speed of (3–3.7) km s−1. A warm kBTe≈6 eV ion population, which also is present, has no influence on the ion acoustic waves due to its low density of only 0.25 % of the plasma density. Near closest approach the propagation direction was within 50∘ from the direction of the bulk velocity. The waves, which in the plasma frame appear below the ion plasma frequency fpi≈2 kHz, are Doppler-shifted to the spacecraft frame where they cover a frequency range up to approximately 4 kHz. The waves are detected in a region of space where the magnetic field is piled up and draped around the inner part of the ionised coma. Estimates of the current associated with the magnetic field gradient as observed by Rosetta are used as input to calculations of dispersion relations for current-driven ion acoustic waves, using kinetic theory. Agreement between theory and observations is obtained for electron and ion distributions with the properties described above. The wave power decreases over cometocentric distances from 24 to 30 km. The main difference between the plasma at closest approach and in the region where the waves are decaying is the absence of a significant current in the latter. Wave observations and theory combined supplement the particle measurements that are difficult at low energies and complicated by spacecraft charging.

Aims: We aim to determine whether dissociative excitation of cometary neutrals by electron impact is the major source of far-ultraviolet (FUV) emissions at comet 67P/Churyumov-Gerasimenko in the southern hemisphere at large heliocentric distances, both during quiet conditions and impacts of corotating interaction regions observed in the summer of 2016.

Methods: We combined multiple datasets from the Rosetta mission through a multi-instrument analysis to complete the first forward modelling of FUV emissions in the southern hemisphere of comet 67P and compared modelled brightnesses to observations with the Alice FUV imaging spectrograph. We modelled the brightness of OI1356, OI1304, Lyman-β, CI1657, and CII1335 emissions, which are associated with the dissociation products of the four major neutral species in the coma: CO₂, H₂O, CO, and O₂. The suprathermal electron population was probed by the Ion and Electron Sensor of the Rosetta Plasma Consortium and the neutral column density was constrained by several instruments: the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA), the Microwave Instrument for the Rosetta Orbiter and the Visual InfraRed Thermal Imaging Spectrometer.

Results: The modelled and observed brightnesses of the FUV emission lines agree closely when viewing nadir and dissociative excitation by electron impact is shown to be the dominant source of emissions away from perihelion. The CII1335 emissions are shown to be consistent with the volume mixing ratio of CO derived from ROSINA. When viewing the limb during the impacts of corotating interaction regions, the model reproduces brightnesses of OI1356 and CI1657 well, but resonance scattering in the extended coma may contribute significantly to the observed Lyman-β and OI1304 emissions. The correlation between variations in the suprathermal electron flux and the observed FUV line brightnesses when viewing the comet's limb suggests electrons are accelerated on large scales and that they originate in the solar wind. This means that the FUV emissions are auroral in nature.

Context: The Rosetta spacecraft escorted Comet 67P/Churyumov-Gerasimenko for 2 yr along its journey through the Solar System between 3.8 and 1.24 au. Thanks to the high resolution mass spectrometer on board Rosetta, the detailed ion composition within a coma has been accurately assessed in situ for the very first time.

Aims: Previous cometary missions, such as Giotto, did not have the instrumental capabilities to identify the exact nature of the plasma in a coma because the mass resolution of the spectrometers onboard was too low to separate ion species with similar masses. In contrast, the Double Focusing Mass Spectrometer (DFMS), part of the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis on board Rosetta (ROSINA), with its high mass resolution mode, outperformed all of them, revealing the diversity of cometary ions.

Methods: We calibrated and analysed the set of spectra acquired by DFMS in ion mode from October 2014 to April 2016. In particular, we focused on the range from 13–39 u q⁻¹. The high mass resolution of DFMS allows for accurate identifications of ions with quasi-similar masses, separating ¹³C⁺ from CH⁺, for instance.

Results: We confirm the presence in situ of predicted cations at comets, such as CH_m⁺ (m = 1−4), H_nO⁺ (n = 1−3), O⁺, Na⁺, and several ionised and protonated molecules. Prior to Rosetta, only a fraction of them had been confirmed from Earth-based observations. In addition, we report for the first time the unambiguous presence of a molecular dication in the gas envelope of a Solar System body, namely CO₂⁺⁺.