In 2004, the Rosetta spacecraft was sent to comet 67P/Churyumov-Gerasimenko for the first ever long-term investigation of a comet. After its arrival in 2014, the spacecraft spent more than two years in immediate proximity to the comet. During these two years, the ROSINA Double Focusing Mass Spectrometer (DFMS) onboard Rosetta discovered a coma with an unexpectedly complex chemical composition that included many oxygenated molecules. Determining the exact cometary composition is an essential first step to understanding of the organic rich chemistry in star forming regions and protoplanetary disks that are ultimately conserved in cometary ices. In this study a joint approach of laboratory calibration and space data analysis was used to perform a detailed identification and quantification of CHO-compounds in the coma of 67P/Churyumov-Gerasimenko. The goal was to derive the CHO-compound abundances relative to water for masses up to 100 u. For this study, the May 2015 post-equinox period represent the best bulk abundances of comet 67P/Churyumov-Gerasimenko. A wide variety of CHO-compounds were discovered and their bulk abundances were derived. Finally, these results are compared to abundances of CHO-bearing molecules in other comets, obtained mostly from ground-based observations and modelling.

Cometary comae are generally depleted in nitrogen. The main carriers for volatile nitrogen in comets are NH3 and HCN. It is known that ammonia readily combines with many acids like e.g. HCN, HNCO, HCOOH, etc. encountered in the interstellar medium as well as in cometary ice to form ammonium salts (NH4+X-) at low temperatures. Ammonium salts, which can play a significant role in prebiotic chemistry, are hard to detect in space as they are unstable in the gas phase and their infrared signature is often hidden by thermal radiation or by e.g. OH in minerals. Here we report the presence of all possible sublimation products of five different ammonium salts at comet 67P/Churyumov-Gerasimenko measured by the ROSINA instrument on Rosetta. The relatively high sublimation temperatures of the salts leads to an apparent lack of volatile nitrogen in the coma. This then also explains the observed trend of higher NH3/H2O ratios with decreasing perihelion distances in comets.

This White Paper outlines the importance of addressing the fundamental science theme <<How are charged particles energized in space plasmas>> through a future ESA mission. The White Paper presents five compelling science questions related to particle energization by shocks, reconnection,waves and turbulence, jets and their combinations. Answering these questions requires resolving scale coupling, nonlinearity and nonstationarity, which cannot be done with existing multi-point observations. In situ measurements from a multi-point, multi-scale L-class plasma observatory consisting of at least 7 spacecraft covering fluid, ion and electron scales are needed. The plasma observatory will enable a paradigm shift in our comprehension of particle energization and space plasma physics in general, with very important impact on solar and astrophysical plasmas. It will be the next logical step following Cluster, THEMIS and MMS for the very large and active European space plasmas community. Being one of the cornerstone missions of the future ESA Voyage 2035-2050 science program, it would further strengthen the European scientific and technical leadership in this important field.

Comets are considered to be some of the most pristine and unprocessed solar system objects accessible to in-situ exploration. Investigating their molecular and elemental composition takes us on a journey back to the early period of our solar system and possibly even further. In this work, we deduce the bulk abundances of the major volatile species in comet 67P/Churyumov-Gerasimenko, the target of the European Space Agency's Rosetta mission. The basis are measurements obtained with the ROSINA instrument suite on board the Rosetta orbiter during a suitable period of high outgassing near perihelion. The results are combined with both gas and dust composition measurements published in the literature. This provides an integrated inventory of the major elements present in the nucleus of 67P/Churyumov-Gerasimenko. Similar to comet 1P/Halley, which was visited by ESA's Giotto spacecraft in 1986, comet 67P/Churyumov-Gerasimenko also shows near-solar abundances of oxygen and carbon, whereas hydrogen and nitrogen are depleted compared to solar. Still, the degree of devolatilization is lower than that of inner solar system objects, including meteorites and the Earth. This supports the idea that comets are among the most pristine objects in our solar system.

A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics.

Comets hold the key to the understanding of our solar system, its formation and its evolution, and to the fundamental plasma processes at work both in it and beyond it. A comet nucleus emits gas as it is heated by the sunlight. The gas forms the coma, where it is ionised, becomes a plasma and eventually interacts with the solar wind. Besides these neutral and ionised gases, the coma also contains dust grains, released from the comet nucleus. As a cometary atmosphere develops when the comet travels through the solar system, large-scale structures, such as the plasma boundaries, develop and disappear, while at planets such large-scale structures are only accessible in their fully grown, quasi-steady state. In situ measurements at comets enable us to learn both how such large-scale structures are formed or reformed and how small-scale processes in the plasma affect the formation and properties of these large scale structures. Furthermore, a comet goes through a wide range of parameter regimes during its life cycle, where either collisional processes, involving neutrals and charged particles, or collisionless processes are at play, and might even compete in complicated transitional regimes. Thus a comet presents a unique opportunity to study this parameter space, from an asteroid-like to a Mars- and Venus-like interaction. Fast flybys of comets have made many new discoveries, setting the stage for a multi-spacecraft mission to accompany a comet on its journey through the solar system. This white paper reviews the present-day knowledge of cometary plasmas, discusses the many questions that remain unanswered, and outlines a multi-spacecraft ESA mission to accompany a comet that will answer these questions by combining both multi-spacecraft observations and a rendezvous mission, and at the same time advance our understanding of fundamental plasma physics and its role in planetary systems.

The nucleus of the Jupiter-family comet 67P/Churyumov-Gerasimenko was discovered to be bi-lobate in shape when the European Space Agency spacecraft Rosetta first approached it in July 2014. The bi-lobate structure of the cometary nucleus has led to much discussion regarding the possible manner of its formation and on how the composition of each lobe might compare with that of the other. During its two-year-long mission from 2014 to 2016, Rosetta remained in close proximity to 67P/Churyumov-Gerasimenko, studying its coma and nucleus in situ. Based on lobe-specific measurements of HDO and H2O performed with the ROSINA DFMS mass spectrometer on board Rosetta, the Deuterium-to-Hydrogen ratios in water from the two lobes could be compared. No appreciable difference was observed, suggesting that both lobes formed in the same region and are homogeneous in their Deuterium-to-Hydrogen ratios.

The European Space Agency (ESA) spacecraft Rosetta accompanied the Jupiter-family comet (JFC) 67P / Churyumov-Gerasimenko for over two years along its trajectory through the inner solar system. Between 2014 and 2016, it performed almost continuous in-situ measurements of the comet's gaseous atmosphere in close proximity to its nucleus. In this study, the $^{16}$O/$^{18}$O ratio of H$_2$O in the coma of 67P / Churyumov-Gerasimenko, as measured by the ROSINA DFMS mass spectrometer on board Rosetta, was determined from the ratio of H${_2}^{16}$O / H${_2}^{18}$O and $^{16}$OH / $^{18}$OH. The value of 445 $\pm$ 35 represents an $\sim$ 11% enrichment of $^{18}$O compared with the terrestrial ratio of 498.7 $\pm$ 0.1. This cometary value is consistent with the comet containing primordial water, in accordance with leading self-shielding models. These models predict primordial water to be between 5% to 20% enriched in heavier oxygen isotopes compared to terrestrial water. Addendum: The $^{16}$O/$^{17}$O ratio of H$_2$O in the coma of 67P / Churyumov-Gerasimenko.

The ROSINA mass spectrometer DFMS on board ESA's Rosetta spacecraft detected the major isotopes of the noble gases argon, krypton, and xenon in the coma of comet 67P/Churyumov-Gerasimenko. Earlier, it has been shown that xenon exhibits an isotopic composition distinct from anywhere else in the solar system. However, argon isotopes, within error, were shown to be consistent with solar isotope abundances. This discrepancy suggested an additional exotic component of xenon in comet 67P/Churyumov-Gerasimenko. Here we show that also krypton exhibits an isotopic composition close to solar. Furthermore, we found a depletion compared to solar of argon with respect to krypton and of krypton with respect to xenon, which is a necessity to postulate an addition of exotic xenon in the comet.

Context: Compensated-current systems are established in response to hot ion beams in terrestrial foreshock regions, around supernova remnants, and in other space and astrophysical plasmas. Aims: We study a non-resonant reactive instability of Alfv\'en waves (AWs) propagating quasi-parallel to the background magnetic field $\mathbf{B}_{0}$ in such systems. Methods: The instability is investigated analytically in the framework of kinetic theory applied to the hydrogen plasmas penetrated by hot proton beams. Results: The instability arises at parallel wavenumbers $k_{z}$ that are sufficiently large to demagnetize the beam ions, $k_{z}V_{Tb}/\omega_{Bi}\gtrsim $ $1$ (here $V_{Tb}$ is the beam thermal speed along $\mathbf{B}_{0}$ and $\omega _{Bi}$ is the ion-cyclotron frequency). The Alfv\'en mode is then made unstable by the imbalance of perturbed currents carried by the magnetized background electrons and partially demagnetized beam ions. The destabilizing effects of the beam temperature and the temperature dependence of the instability threshold and growth rate are demonstrated for the first time. The beam temperature, density, and bulk speed are all destabilizing and can be combined in a single destabilizing factor $\alpha_{b}$ triggering the instability at {$\alpha _{b}>$ $\alpha _{b}^{\mathrm{thr}}$}, where the threshold varies in a narrow range $2.43\leq $ $\alpha _{b}^{\mathrm{thr}}\leq $ $4.87$. New analytical expressions for the instability growth rate and its boundary in the parameter space are obtained and can be directly compared with observations. Two applications to terrestrial foreshocks and foreshocks around supernova remnants are shortly discussed. In particular, our results suggest that the ions reflected by the shocks around supernova remnants can drive stronger instability than the cosmic rays.

Aims. To better understand the radiation environment in low Earth orbit (LEO), the analysis of in-situ observations of a variety of particles, at different atmospheric heights, and in a wide range of energies, is needed. Methods. We present an analysis of energetic particles, indirectly detected by the Large Yield RAdiometer (LYRA) instrument on board ESA's Project for On-board Autonomy 2 (PROBA2) satellite as background signal. Combining Energetic Particle Telescope (EPT) observations with LYRA data for an overlapping period of time, we identified these particles as electrons with an energy range of 2 to 8 MeV. Results. The observed events are strongly correlated to geo-magnetic activity and appear even during modest disturbances. They are also well confined geographically within the L=4-6 McIlwain zone, which makes it possible to identify their source. Conclusions. Although highly energetic particles are commonly perturbing data acquisition of space instruments, we show in this work that ultra-relativistic electrons with energies in the range of 2-8 MeV are detected only at high latitudes, while not present in the South Atlantic Anomaly region.

Context. The Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) was designed to measure the composition of the gas in the coma of comet 67P/Churyumov-Gerasimenko, the target of the European Space Agency's Rosetta mission. In addition to the volatiles, ROSINA measured refractories sputtered off the comet by the interaction of solar wind protons with the surface of the comet. Aims. The origin of different solar system materials is still heavily debated. Isotopic ratios can be used to distinguish between different reservoirs and investigate processes occurring during the formation of the solar system. Methods. ROSINA consisted of two mass spectrometers and a pressure sensor. In the ROSINA Double Focusing Mass Spectrometer (DFMS), the neutral gas of cometary origin was ionized and then deflected in an electric and a magnetic field that separated the ions based on their mass-to-charge ratio. The DFMS had a high mass resolution, dynamic range, and sensitivity that allowed detection of rare species and the known major volatiles. Results. We measured the relative abundance of all three stable silicon isotopes with the ROSINA instrument on board the Rosetta spacecraft. Furthermore, we measured $^{13}$C/$^{12}$C in C$_2$H$_4$, C$_2$H$_5$, and CO. The DFMS in situ measurements indicate that the average silicon isotopic composition shows depletion in the heavy isotopes $^{29}$Si and $^{30}$Si with respect to $^{28}$Si and solar abundances, while $^{13}$C to $^{12}$C is analytically indistinguishable from bulk planetary and meteorite compositions. Although the origin of the deficiency of the heavy silicon isotopes cannot be explained unambiguously, we discuss mechanisms that could have contributed to the measured depletion of the isotopes $^{29}$Si and $^{30}$Si.

Because of the high fraction of refractory material present in comets, the heat produced by the radiogenic decay of elements such as aluminium and iron can be high enough to induce the loss of ultravolatile species such as nitrogen, argon or carbon monoxide during their accretion phase in the protosolar nebula. Here, we investigate how heat generated by the radioactive decay of 26Al and 60Fe influences the formation of comet 67P/Churyumov-Gerasimenko, as a function of its accretion time and size of parent body. We use an existing thermal evolution model that includes various phase transitions, heat transfer in the ice-dust matrix, and gas diffusion throughout the porous material, based on thermodynamic parameters derived from Rosetta observations. Two possibilities are considered: either, to account for its bilobate shape, 67P/Churyumov-Gerasimenko was assembled from two primordial ~2 kilometer-sized planetesimals, or it resulted from the disruption of a larger parent body with a size corresponding to that of comet Hale-Bopp (~70 km). To fully preserve its volatile content, we find that either 67P/Churyumov-Gerasimenko's formation was delayed between ~2.2 and 7.7 Myr after that of Ca-Al-rich Inclusions in the protosolar nebula or the comet's accretion phase took place over the entire time interval, depending on the primordial size of its parent body and the composition of the icy material considered. Our calculations suggest that the formation of 67P/Churyumov-Gerasimenko is consistent with both its accretion from primordial building blocks formed in the nebula or from debris issued from the disruption of a Hale-Bopp-like body.

Compensated-current systems created by energetic ion beams are widespread in space and astrophysical plasmas. The well-known examples are foreshock regions in the solar wind and around supernova remnants. We found a new oblique Alfv\'enic instability driven by compensated currents flowing along the background magnetic field. Because of the vastly different electron and ion gyroradii, oblique Alfv\'enic perturbations react differently on the currents carried by the hot ion beams and the return electron currents. Ultimately, this difference leads to a non-resonant aperiodic instability at perpendicular wavelengths close to the beam ion gyroradius. The instability growth rate increases with increasing beam current and temperature. In the solar wind upstream of Earth's bow shock the instability growth time can drop below 10 proton cyclotron periods. Our results suggest that this instability can contribute to the turbulence and ion acceleration in space and astrophysical foreshocks.

We investigate electrostatic plasma instabilities of Farley-Buneman (FB) type driven by quasi-stationary neutral gas flows in the solar chromosphere. The role of these instabilities in the chromosphere is clarified. We find that the destabilizing ion thermal effect is highly reduced by the Coulomb collisions and can be ignored for the chromospheric FB-type instabilities. On the contrary, the destabilizing electron thermal effect is important and causes a significant reduction of the neutral drag velocity triggering the instability. The resulting threshold velocity is found as function of chromospheric height. Our results indicate that the FB type instabilities are still less efficient in the global chromospheric heating than the Joule dissipation of the currents driving these instabilities. This conclusion does not exclude the possibility that the FB type instabilities develop in the places where the cross-field currents overcome the threshold value and contribute to the heating locally. Typical length-scales of plasma density fluctuations produced by these instabilities are determined by the wavelengths of unstable modes, which are in the range $10-10^2$ cm in the lower chromosphere, and $10^2-10^3$ cm in the upper chromosphere. These results suggest that the decimetric radio waves undergoing scattering (scintillations) by these plasma irregularities can serve as a tool for remote probing of the solar chromosphere at different heights.

A weakly dispersive sub-range (WDR) of kinetic Alfv\'en turbulence is distinguished and investigated for the first time in the context of MHD/kinetic turbulence transition. We found perpendicular wavenumber spectra ~ k^{-3} and ~ k^{-4} formed in WDR by strong and weak turbulence of kinetic Alfv\'en waves (KAWs), respectively. These steep WDR spectra connect shallower spectra in the MHD and strongly dispersive KAW sub-ranges, which results in a specific double-kink (2-k) pattern often seen in observed turbulent spectra. The first kink occurs where MHD turbulence transforms into weakly dispersive KAW turbulence; the second one is between weakly and strongly dispersive KAW sub-ranges. Our analysis suggests that the partial turbulence dissipation due to amplitude-dependent super-adiabatic ion heating may occur in the vicinity of the first spectral kink. A threshold-like nature of this process results in a conditional selective dissipation affecting only largest over-threshold amplitudes and decreasing intermittency in the range below the first spectral kink. Several recent counter-intuitive observational findings can be explained by the selective dissipation coupled to the nonlinear interaction among weakly dispersive KAWs.