Comet 1P/Halley multifluid MHD model for the Giotto fly-by

The interaction of comets with the solar wind has been the focus of many studies including numerical modeling. We compare the results of our multifluid MHD simulation of comet 1P/Halley to data obtained during the flyby of the European Space Agency's Giotto spacecraft in 1986. The model solves the full set of MHD equations for the individual fluids representing the solar wind protons, the cometary light and heavy ions, and the electrons. The mass loading, charge-exchange, dissociative ion-electron recombination, and collisional interactions between the fluids are taken into account. The computational domain spans over several million kilometers, and the close vicinity of the comet is resolved to the details of the magnetic cavity. The model is validated by comparison to the corresponding Giotto observations obtained by the Ion Mass Spectrometer, the Neutral Mass Spectrometer, the Giotto magnetometer experiment, and the Johnstone Plasma Analyzer instrument. The model shows the formation of the bow shock, the ion pile-up, and the diamagnetic cavity and is able to reproduce the observed temperature differences between the pick-up ion populations and the solar wind protons. We give an overview of the global interaction of the comet with the solar wind and then show the effects of the Lorentz force interaction between the different plasma populations.

Plasma environment of a weak comet – Predictions for Comet 67P/Churyumov–Gerasimenko from multifluid-MHD and Hybrid models

The interaction of a comet with the solar wind undergoes various stages as the comet’s activity varies along its orbit. For a comet like 67P/Churyumov–Gerasimenko, the target comet of ESA’s Rosetta mission, the various features include the formation of a Mach cone, the bow shock, and close to perihelion even a diamagnetic cavity. There are different approaches to simulate this complex interplay between the solar wind and the comet’s extended neutral gas coma which include magnetohydrodynamics (MHD) and hybrid-type models. The first treats the plasma as fluids (one fluid in basic single fluid MHD) and the latter treats the ions as individual particles under the influence of the local electric and magnetic fields. The electrons are treated as a charge-neutralizing fluid in both cases. Given the different approaches both models yield different results, in particular for a low production rate comet. In this paper we will show that these differences can be reduced when using a multifluid instead of a single-fluid MHD model and increase the resolution of the Hybrid model. We will show that some major features obtained with a hybrid type approach like the gyration of the cometary heavy ions and the formation of the Mach cone can be partially reproduced with the multifluid-type model.

Molecular nitrogen in comet 67P/Churyumov-Gerasimenko indicates a low formation temperature

Molecular nitrogen (N2) is thought to have been the most abundant form of nitrogen in the protosolar nebula. It is the main N-bearing molecule in the atmospheres of Pluto and Triton and probably the main nitrogen reservoir from which the giant planets formed. Yet in comets, often considered the most primitive bodies in the solar system, N2 has not been detected. Here we report the direct in situ measurement of N2 in the Jupiter family comet 67P/Churyumov-Gerasimenko, made by the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis mass spectrometer aboard the Rosetta spacecraft. A N2/CO ratio of Embedded Image (2σ standard deviation of the sampled mean) corresponds to depletion by a factor of ~25.4 ± 8.9 as compared to the protosolar value. This depletion suggests that cometary grains formed at low-temperature conditions below ~30 kelvin.

Time variability and heterogeneity in the coma of 67P/Churyumov-Gerasimenko

Comets contain the best-preserved material from the beginning of our planetary system. Their nuclei and comae composition reveal clues about physical and chemical conditions during the early solar system when comets formed. ROSINA (Rosetta Orbiter Spectrometer for Ion and Neutral Analysis) onboard the Rosetta spacecraft has measured the coma composition of comet 67P/Churyumov-Gerasimenko with well-sampled time resolution per rotation. Measurements were made over many comet rotation periods and a wide range of latitudes. These measurements show large fluctuations in composition in a heterogeneous coma that has diurnal and possibly seasonal variations in the major outgassing species: water, carbon monoxide, and carbon dioxide. These results indicate a complex coma-nucleus relationship where seasonal variations may be driven by temperature differences just below the comet surface.

Inventory of the volatiles on comet 67P/Churyumov-Gerasimenko from Rosetta/ROSINA

Context. The ESA Rosetta spacecraft (S/C) is tracking comet 67P/Churyumov-Gerasimenko in close vicinity. This prolonged encounter enables studying the evolution of the volatile coma composition. Aims. Our work aims at comparing the diversity of the coma of 67P/Churyumov-Gerasimenko at large heliocentric distance to study the evolution of the comet during its passage around the Sun and at trying to classify it relative to other comets. Methods. We used the Double Focussing Mass Spectrometer (DFMS) of the ROSINA experiment on ESA’s Rosetta mission to determine relative abundances of major and minor volatile species. This study is restricted to species that have previously been detected elsewhere. Results. We detect almost all species currently known to be present in cometary coma with ROSINA DFMS. As DFMS measured the composition locally, we cannot derive a global abundance, but we compare measurements from the summer and the winter hemisphere with known abundances from other comets. Differences between relative abundances between summer and winter hemispheres are large, which points to a possible evolution of the cometary surface. This comet appears to be very rich in CO2 and ethane. Heavy oxygenated compounds such as ethylene glycol are underabundant at 3 AU, probably due to their high sublimation temperatures, but nevertheless, their presence proves that Kuiper belt comets also contain complex organic molecules.

Molecular oxygen in oort cloud comet 1P/Halley

Recently, the ROSINA mass spectrometer suite on board the European Space Agency's Rosetta spacecraft discovered an abundant amount of molecular oxygen, O2, in the coma of Jupiter family comet 67P/Churyumov-Gerasimenko of O2/H2O = 3.80 ± 0.85%. It could be shown that O2 is indeed a parent species and that the derived abundances point to a primordial origin. Crucial questions are whether the O2 abundance is peculiar to comet 67P/Churyumov-Gerasimenko or Jupiter family comets in general, and also whether Oort cloud comets such as comet 1P/Halley contain similar amounts of molecular oxygen. We investigated mass spectra obtained by the Neutral Mass Spectrometer instrument during the flyby by the European Space Agency's Giotto probe of comet 1P/Halley. Our investigation indicates that a production rate of O2 of 3.7 ± 1.7% with respect to water is indeed compatible with the obtained Halley data and therefore that O2 might be a rather common and abundant parent species.

Modeled Interaction of Comet 67P/Churyumov-Gerasimenko with the Solar Wind Inside 2 AU

Periodic comets move around the Sun on elliptical orbits. As such comet 67P/Churyumov-Gerasimenko (hereafter 67P) spends a portion of time in the inner solar system where it is exposed to increased solar insolation. Therefore given the change in heliocentric distance, in case of 67P from aphelion at 5.68 AU to perihelion at ~1.24 AU, the comet’s activity—the production of neutral gas and dust—undergoes significant variations. As a consequence, during the inbound portion, the mass loading of the solar wind increases and extends to larger spatial scales. This paper investigates how this interaction changes the character of the plasma environment of the comet by means of multifluid MHD simulations. The multifluid MHD model is capable of separating the dynamics of the solar wind ions and the pick-up ions created through photoionization and electron impact ionization in the coma of the comet. We show how two of the major boundaries, the bow shock and the diamagnetic cavity, form and develop as the comet moves through the inner solar system. Likewise for 67P, although most likely shifted back in time with respect to perihelion passage, this process is reversed on the outbound portion of the orbit. The presented model herein is able to reproduce some of the key features previously only accessible to particle-based models that take full account of the ions’ gyration. The results shown herein are in decent agreement to these hybrid-type kinetic simulations.

In situ plasma measurements of fragmented comet 73P Schwassmann–Wachmann 3

The interiors of comets contain some of the most pristine material in the solar system. Comet 73P/Schwassmann–Wachmann 3, discovered in 1930, is a Jupiter-family comet with a 5.34-year period. This comet split into 5 fragments in 1995 and disintegrated into nearly 70 major pieces in 2006. In 2006 May and June, recently ionized cometary particles originating from fragments including and surrounding some of these major objects were collected with the ACE/SWICS and Wind/STICS sensors. Due to a combination of the instrument characteristics and the close proximity of the fragments passing between those spacecraft and the Sun, unique measurements regarding the charge state composition and the elemental abundances of both cometary and heliospheric plasma were made during that time. The cometary material released from some of these fragments can be identified by the concentrations of water-group pickup ions having a mass-per-charge ratio of 16–18 amu e−1, indicating that while these fragments are small, they are still actively sublimating. We present an analysis of cometary composition, spatial distribution, and heliospheric interactions, with a focus on helium, C+/O+, and water-group ions.

The capabilities of ROSINA/DFMS to measure argon isotopes at comet 67P/Churyumov–Gerasimenko

Little is known about the noble gas abundances in comets. These highly volatile atoms are possible tracers of the history of cometary matter including the thermal evolution. They can help quantify the contribution of cometary impacts to terrestrial oceans and help elucidate on the formation history of comets and their role in the formation and evolution of planetary atmospheres. This paper focuses on argon and the capabilities to measure this noble gas with in situ mass spectrometry at comet 67P/Churyumov–Gerasimenko, the target of the European Space Agency׳s spacecraft Rosetta. Argon may have been detected by remote sensing in a single Oort cloud comet but to date nothing is known about the isotopic abundances of argon in comets. Furthermore, no detection of argon in a Jupiter-family comet has been reported. Comet 67P/Churyumov–Gerasimenko belongs to the group of Jupiter-family comets and originates most likely in the Kuiper belt. Onboard Rosetta is ROSINA/DFMS (Rosetta Orbiter Spectrometer for Ion and Neutral Analysis/Double Focusing Mass Spectrometer). DFMS has unprecedented mass resolution and high sensitivity and is designed to measure isotopic ratios including argon (Balsiger et al., 2007). Argon measurements using the DFMS lab model (identical to the flight model) demonstrate this capability. At very least, this mass spectrometer has the resolution and sensitivity to reduce the upper limit on the argon outgassing rate relative to water by more than three orders of magnitude (for 38Ar). Most likely, ROSINA/DFMS will provide the first detection of argon in a Jupiter-family comet together with the first determination of the ³⁶Ar/³⁸Ar ratio at a comet.

Sulfur Isotopic Ratios at 67P/Churyumov-Gerasimenko and Characterization of ROSINA-DFMS FM & FS

Comets are thought to be the most pristine bodies present in the Solar System. In consequence of spending the majority of their existence beyond 30 AU, their composition can give insights on the physical and chemical conditions during their formation. Since August 2014 the European Space Agency spacecraft Rosetta accompanies the Jupiter family comet 67P/Churyumov-Gerasimenko on its way to perihelion and beyond. In this study the isotope fractionation of 34S are reported in H2S, OCS, SO2, S2, and CS2 at 67P. In addition for the first time the isotope fractionation for 33S is presented for cometary volatiles. The ratio 32S/33S is given for H2S, SO2 and a tentative value is given for CS2. With a mean value of -50 ± 22‰ and -306 ± 31‰ for δ34S and δ33S respectively, H2S shows a significant depletion in both 34S and 33S. For SO2 the depletion is less distinct with δ34S and δ33S being -67 ± 40‰ and -130 ± 53‰, respectively. The strongest depletion is present for CS2 with -114 ± 21‰and -276 ± 55‰, respectively. For OCS and S2 only δ34S could be determined which is -252 ± 77‰ and -357 ± 145‰, respectively. A comparison with sulfur isotopic ratios measured in SiC grains revealed that both SiC grains and the five volatile species have similar sulfur isotopic ratios. However, it is beyond the scope of this work to investigate the possibility of a link between SiC grains and cometary ices. Nevertheless, mass-dependent or mass-independent fractionation due to photo dissociation can be ruled out as sole cause of the seen depletion of 33S and 34S. Furthermore, an upper limit of (9.64 ± 0.19)·10.4 for D/H in HDS has been determined. This value is about a factor two higher than D/H in H2O for the same comet reported by (Altwegg et al., 2015). Besides the investigation concerning isotopic ratios of sulfur bearing species in this work the calibration and characterization of ROSINA/DFMS has been continued. Here it is reported about the deviation of the mass scale for MCP/LEDA low resolution spectra and the calibration measurements performed in the laboratory. Furthermore the outcome of the attempt to describe the sensitivity of DFMS with an empirical function will be discussed. The last part of the characterization of DFMS is dedicated to determine the so-called individual pixel gain for the laboratory and the flight model. Moreover, correlation between the depletion’s manifestation of the MCP with respect to the applied voltages has been investigated for both models. It has been found that further measurements are needed to understand the manifestation of depletion at the laboratory model. For the model on board of Rosetta it could be shown that most of the present feature are due to the usage of the MCP and suggestions have been made in order to answer the remaining question considering the depletion of the MCP.

Cometary Isotopic Measurements

Isotopic ratios in comets provide keys for the understanding of the origin of cometary material, and the physical and chemical conditions in the early Solar Nebula. We review here measurements acquired on the D/H, N-14/N-15, O-16/O-18, C-12/C-13, and S-32/S-34 ratios in cometary grains and gases, and discuss their cosmogonic implications. The review includes analyses of potential cometary material available in collections on Earth, recent measurements achieved with the Herschel Space Observatory, large optical telescopes, and Rosetta, as well as recent results obtained from models of chemical-dynamical deuterium fractionation in the early solar nebula. Prospects for future measurements are presented.

New fully kinetic model for the study of electric potential, plasma, and dust above lunar landscapes

We have developed a new fully kinetic electrostatic simulation, HYBes, to study how the lunar landscape affects the electric potential and plasma distributions near the surface and the properties of lifted dust. The model embodies new techniques that can be used in various types of physical environments and situations. We demonstrate the applicability of the new model in a situation involving three charged particle species, which are solar wind electrons and protons, and lunar photoelectrons. Properties of dust are studied with test particle simulations by using the electric fields derived from the HYBes model. Simulations show the high importance of the plasma and the electric potential near the surface. For comparison, the electric potential gradients near the landscapes with feature sizes of the order of the Debye length are much larger than those near a flat surface at different solar zenith angles. Furthermore, dust test particle simulations indicate that the landscape relief influences the dust location over the surface. The study suggests that the local landscape has to be taken into account when the distributions of plasma and dust above lunar surface are studied. The HYBes model can be applied not only at the Moon but also on a wide range of airless planetary objects such as Mercury, other planetary moons, asteroids, and nonactive comets.

Gravitational slopes, geomorphology, and material strengths of the nucleus of comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Aims. We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov-Gerasimenko and provide constraints on the mechanical properties of the cometary material (tensile, shear, and compressive strengths). Methods. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface (Imhotep, Ash, Seth, Hathor, and Agilkia), using two shape models computed from OSIRIS images by the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive strengths using different surface morphologies (overhangs, collapsed structures, boulders, cliffs, and Philae's footprint) and mechanical considerations. Results. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain morphology: i) low-slope terrains (0-20 degrees) are covered by a fine material and contain a few large (>10 m) and isolated boulders; ii) intermediate-slope terrains (20-45 degrees) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from <1m to 10m for the majority of them; and iii) high-slope terrains (45-90 degrees) are cliffs that expose a consolidated material and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3-15 Pa (upper limit of 150 Pa), 4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (upper limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and tends to favor a formation of comets by the accrection of pebbles at low velocities.

A ~ 32-70 K formation temperature range for the ice grains agglomerated by Comet 67P/Churyumov-Gerasimenko

Grand Canonical Monte Carlo simulations are used to reproduce the N₂/CO ratio ranging between 1.7 x 10⁻³ and 1.6 x 10⁻² observed in situ in the Jupiter-family comet 67 P/Churyumov-Gerasimenko (67 P) by the ROSINA mass spectrometer on board the Rosetta spacecraft. By assuming that this body has been agglomerated from clathrates in the protosolar nebula (PSN), simulations are developed using elaborated interatomic potentials for investigating the temperature dependence of the trapping within a multiple-guest clathrate formed from a gas mixture of CO and N₂ in proportions corresponding to those expected for the PSN. By assuming that 67 P agglomerated from clathrates, our calculations suggest the cometary grains must have been formed at temperatures ranging between ~ 31.8 and 69.9 K in the PSN to match the N₂/CO ratio measured by the ROSINA mass spectrometer. The presence of clathrates in Jupiter-family comets could then explain the potential N₂ depletion (factor of up to ~ 87 compared to the protosolar value) measured in 67 P/Churyumov-Gerasimenko.

Comet 67P/Churyumov-Gerasimenko: Constraints on its origin from OSIRIS observations

Context. One of the main aims of the ESA Rosetta mission is to study the origin of the solar system by exploring comet 67P/Churyumov-Gerasimenko at close range. Aims. In this paper we discuss the origin and evolution of comet 67P/Churyumov-Gerasimenko in relation to that of comets in general and in the framework of current solar system formation models. Methods. We use data from the OSIRIS scientific cameras as basic constraints. In particular, we discuss the overall bi-lobate shape and the presence of key geological features, such as layers and fractures. We also treat the problem of collisional evolution of comet nuclei by a particle-in-a-box calculation for an estimate of the probability of survival for 67P/Churyumov-Gerasimenko during the early epochs of the solar system. Results. We argue that the two lobes of the 67P/Churyumov-Gerasimenko nucleus are derived from two distinct objects that have formed a contact binary via a gentle merger. The lobes are separate bodies, though sufficiently similar to have formed in the same environment. An estimate of the collisional rate in the primordial, trans-planetary disk shows that most comets of similar size to 67P/Churyumov-Gerasimenko are likely collisional fragments, although survival of primordial planetesimals cannot be excluded. Conclusions. A collisional origin of the contact binary is suggested, and the low bulk density of the aggregate and abundance of volatile species show that a very gentle merger must have occurred. We thus consider two main scenarios: the primordial accretion of planetesimals, and the re-accretion of fragments after an energetic impact onto a larger parent body. We point to the primordial signatures exhibited by 67P/Churyumov-Gerasimenko and other comet nuclei as critical tests of the collisional evolution.