Assessment of detectability of neutral interstellar deuterium by IBEX observations

Context. The abundance of deuterium in the interstellar gas in front of the Sun gives insight into the processes of filtration of neutral interstellar species through the heliospheric interface and potentially into the chemical evolution of the Galactic gas. Aims: We investigate the possibility of detection of neutral interstellar deuterium at 1 AU from the Sun by direct sampling by the Interstellar Boundary Explorer (IBEX). Methods: Using both previous and the most recent determinations of the flow parameters of neutral gas in the local interstellar cloud (LIC) and an observation-based model of solar radiation pressure and ionization in the heliosphere, we simulated the flux of neutral interstellar D at IBEX for the actual measurement conditions. We assessed the number of interstellar D atom counts expected during the first three years of IBEX operation. We also simulated the observations expected during an epoch of high solar activity. In addition, we calculated the expected counts of D atoms from the thin terrestrial water layer covering the IBEX-Lo conversion surface, sputtered by neutral interstellar He atoms. Results: Most D counts registered by IBEX-Lo are expected to come from the water layer, exceeding the interstellar signal by 2 orders of magnitude. However, the sputtering should stop once the Earth leaves the portion of orbit traversed by interstellar He atoms. We identify seasons during the year when mostly the genuine interstellar D atoms are expected in the signal. During the first 3 years of IBEX operations about 2 detectable interstellar D atoms are expected. This number is comparable to the expected number of sputtered D atoms registered during the same time intervals. Conclusions: The most favorable conditions for the detection occur during low solar activity, in an interval including March and April each year. The detection chances could be improved by extending the instrument duty cycle, say, by making observations in the special deuterium mode of IBEX-Lo.

Theoretical models of planetary system formation: mass vs. semi-major axis

Context. Planet formation models have been developed during the past years to try to reproduce what has been observed of both the solar system and the extrasolar planets. Some of these models have partially succeeded, but they focus on massive planets and, for the sake of simplicity, exclude planets belonging to planetary systems. However, more and more planets are now found in planetary systems. This tendency, which is a result of radial velocity, transit, and direct imaging surveys, seems to be even more pronounced for low-mass planets. These new observations require improving planet formation models, including new physics, and considering the formation of systems. Aims: In a recent series of papers, we have presented some improvements in the physics of our models, focussing in particular on the internal structure of forming planets, and on the computation of the excitation state of planetesimals and their resulting accretion rate. In this paper, we focus on the concurrent effect of the formation of more than one planet in the same protoplanetary disc and show the effect, in terms of architecture and composition of this multiplicity. Methods: We used an N-body calculation including collision detection to compute the orbital evolution of a planetary system. Moreover, we describe the effect of competition for accretion of gas and solids, as well as the effect of gravitational interactions between planets. Results: We show that the masses and semi-major axes of planets are modified by both the effect of competition and gravitational interactions. We also present the effect of the assumed number of forming planets in the same system (a free parameter of the model), as well as the effect of the inclination and eccentricity damping. We find that the fraction of ejected planets increases from nearly 0 to 8% as we change the number of embryos we seed the system with from 2 to 20 planetary embryos. Moreover, our calculations show that, when considering planets more massive than ~5 M⊕, simulations with 10 or 20 planetary embryos statistically give the same results in terms of mass function and period distribution.

Planet formation models: the interplay with the planetesimal disc

Context. According to the sequential accretion model (or core-nucleated accretion model), giant planet formation is based first on the formation of a solid core which, when massive enough, can gravitationally bind gas from the nebula to form the envelope. The most critical part of the model is the formation time of the core: to trigger the accretion of gas, the core has to grow up to several Earth masses before the gas component of the protoplanetary disc dissipates. Aims: We calculate planetary formation models including a detailed description of the dynamics of the planetesimal disc, taking into account both gas drag and excitation of forming planets. Methods: We computed the formation of planets, considering the oligarchic regime for the growth of the solid core. Embryos growing in the disc stir their neighbour planetesimals, exciting their relative velocities, which makes accretion more difficult. Here we introduce a more realistic treatment for the evolution of planetesimals' relative velocities, which directly impact on the formation timescale. For this, we computed the excitation state of planetesimals, as a result of stirring by forming planets, and gas-solid interactions. Results: We find that the formation of giant planets is favoured by the accretion of small planetesimals, as their random velocities are more easily damped by the gas drag of the nebula. Moreover, the capture radius of a protoplanet with a (tiny) envelope is also larger for small planetesimals. However, planets migrate as a result of disc-planet angular momentum exchange, with important consequences for their survival: due to the slow growth of a protoplanet in the oligarchic regime, rapid inward type I migration has important implications on intermediate-mass planets that have not yet started their runaway accretion phase of gas. Most of these planets are lost in the central star. Surviving planets have masses either below 10 M⊕ or above several Jupiter masses. Conclusions: To form giant planets before the dissipation of the disc, small planetesimals (~0.1 km) have to be the major contributors of the solid accretion process. However, the combination of oligarchic growth and fast inward migration leads to the absence of intermediate-mass planets. Other processes must therefore be at work to explain the population of extrasolar planets that are presently known.

Evidence for distinct modes of solar activity

The Sun shows strong variability in its magnetic activity, from Grand minima to Grand maxima, but the nature of the variability is not fully understood, mostly because of the insufficient length of the directly observed solar activity records and of uncertainties related to long-term reconstructions. Here we present a new adjustment-free reconstruction of solar activity over three millennia and study its different modes. Methods. We present a new adjustment-free, physical reconstruction of solar activity over the past three millennia, using the latest verified carbon cycle, 14C production, and archeomagnetic field models. This great improvement allowed us to study different modes of solar activity at an unprecedented level of details. Results. The distribution of solar activity is clearly bi-modal, implying the existence of distinct modes of activity. The main regular activity mode corresponds to moderate activity that varies in a relatively narrow band between sunspot numbers 20 and 67. The existence of a separate Grand minimum mode with reduced solar activity, which cannot be explained by random fluctuations of the regular mode, is confirmed at a high confidence level. The possible existence of a separate Grand maximum mode is also suggested, but the statistics is too low to reach a confident conclusion. Conclusions. The Sun is shown to operate in distinct modes – a main general mode, a Grand minimum mode corresponding to an inactive Sun, and a possible Grand maximum mode corresponding to an unusually active Sun. These results provide important constraints for both dynamo models of Sun-like stars and investigations of possible solar influence on Earth’s climate.

The HARPS search for southern extra-solar planets

The meter-per-second precision achieved by today’s velocimeters enables us to search for 1−10 M⊕ planets in the habitable zone of cool stars. This paper reports on the detection of three planets orbiting GJ 163 (HIP 19394), a M3 dwarf monitored by our ESO/HARPS search for planets. We made use of the HARPS spectrograph to collect 150 radial velocities of GJ 163 over a period of eight years. We searched the radial-velocity time series for coherent signals and found five distinct periodic variabilities. We investigated the stellar activity and called into question the planetary interpretation for two signals. Before more data can be acquired we concluded that at least three planets are orbiting GJ 163. They have orbital periods of Pb = 8.632 ± 0.002, Pc = 25.63 ± 0.03, and Pd = 604 ± 8 days and minimum masses msini = 10.6 ± 0.6, 6.8 ± 0.9, and 29 ± 3 M⊕, respectively. We hold our interpretations for the two additional signals with periods P(e) = 19.4 and P(f) = 108 days. The inner pair presents an orbital period ratio of 2.97, but a dynamical analysis of the system shows that it lays outside the 3:1 mean motion resonance. The planet GJ 163c, in particular, is a super-Earth with an equilibrium temperature of Teq = (302 ± 10)(1 − A)1/4 K and may lie in the so-called habitable zone for albedo values (A = 0.34 − 0.89) moderately higher than that of Earth (A⊕ = 0.2−0.3).

The HARPS search for southern extra-solar planets

The vast diversity of planetary systems detected to date is defying our capability of understanding their formation and evolution. Well-defined volume-limited surveys are the best tool at our disposal to tackle the problem, via the acquisition of robust statistics of the orbital elements. We are using the HARPS spectrograph to conduct our survey of ≈850 nearby solar-type stars, and in the course of the past nine years we have monitored the radial velocity of HD 103774, HD 109271, and BD-061339. In this work we present the detection of five planets orbiting these stars, with msin   (i) between 0.6 and 7 Neptune masses, four of which are in two multiple systems, comprising one super-Earth and one planet within the habitable zone of a late-type dwarf. Although for strategic reasons we chose efficiency over precision in this survey, we have the capability to detect planets down to the Neptune and super-Earth mass range as well as multiple systems, provided that enough data points are made available.

From stellar nebula to planets: The refractory components

Context. To date, calculations of planet formation have mainly focused on dynamics, and only a few have considered the chemical composition of refractory elements and compounds in the planetary bodies. While many studies have been concentrating on the chemical composition of volatile compounds (such as H2O, CO, CO2) incorporated in planets, only a few have considered the refractory materials as well, although they are of great importance for the formation of rocky planets. Aims. We computed the abundance of refractory elements in planetary bodies formed in stellar systems with a solar chemical composition by combining models of chemical composition and planet formation. We also considered the formation of refractory organic compounds, which have been ignored in previous studies on this topic. Methods. We used the commercial software package HSC Chemistry to compute the condensation sequence and chemical composition of refractory minerals incorporated into planets. The problem of refractory organic material is approached with two distinct model calculations: the first considers that the fraction of atoms used in the formation of organic compounds is removed from the system (i.e., organic compounds are formed in the gas phase and are non-reactive); and the second assumes that organic compounds are formed by the reaction between different compounds that had previously condensed from the gas phase. Results. Results show that refractory material represents more than 50 wt % of the mass of solids accreted by the simulated planets with up to 30 wt % of the total mass composed of refractory organic compounds. Carbide and silicate abundances are consistent with C/O and Mg/Si elemental ratios of 0.5 and 1.02 for the Sun. Less than 1 wt % of carbides are present in the planets, and pyroxene and olivine are formed in similar quantities. The model predicts planets that are similar in composition to those of the solar system. Starting from a common initial nebula composition, it also shows that a wide variety of chemically different planets can form, which means that the differences in planetary compositions are due to differences in the planetary formation process. Conclusions. We show that a model in which refractory organic material is absent from the system is more compatible with observations. The use of a planet formation model is essential to form a wide diversity of planets in a consistent way.

The rotation state of 67P/Churyumov-Gerasimenko from approach observations with the OSIRIS cameras on Rosetta

Aims. Approach observations with the Optical, Spectroscopic, and Infrared Remote Imaging System (OSIRIS) experiment onboard Rosetta are used to determine the rotation period, the direction of the spin axis, and the state of rotation of comet 67P’s nucleus. Methods. Photometric time series of 67P have been acquired by OSIRIS since the post wake-up commissioning of the payload in March 2014. Fourier analysis and convex shape inversion methods have been applied to the Rosetta data as well to the available ground-based observations. Results. Evidence is found that the rotation rate of 67P has significantly changed near the time of its 2009 perihelion passage, probably due to sublimation-induced torque. We find that the sidereal rotation periods P1 = 12.76129 ± 0.00005 h and P2 = 12.4043 ± 0.0007 h for the apparitions before and after the 2009 perihelion, respectively, provide the best fit to the observations. No signs of multiple periodicity are found in the light curves down to the noise level, which implies that the comet is presently in a simple rotation state around its axis of largest moment of inertia. We derive a prograde rotation model with spin vector J2000 ecliptic coordinates λ = 65° ± 15°, β = + 59° ± 15°, corresponding to equatorial coordinates RA = 22°, Dec = + 76°. However, we find that the mirror solution, also prograde, at λ = 275° ± 15°, β = + 50° ± 15° (or RA = 274°, Dec = + 27°), is also possible at the same confidence level, due to the intrinsic ambiguity of the photometric problem for observations performed close to the ecliptic plane.

From planetesimals to planets: volatile molecules

Context. Solar and extrasolar planets are the subject of numerous studies aiming to determine their chemical composition and internal structure. In the case of extrasolar planets, the composition is important as it partly governs their potential habitability. Moreover, observational determination of chemical composition of planetary atmospheres are becoming available, especially for transiting planets. Aims. The present works aims at determining the chemical composition of planets formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data for models of planet formation and evolution, and future interpretation of chemical composition of solar and extrasolar planets. Methods. We have developed a model that computes the composition of ices in planets in different stellar systems with the use of models of ice and planetary formation. Results. We provide the chemical composition, ice/rock mass ratio and C:O molar ratio for planets in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, we produce a wide variety of planetary chemical compositions as a function of the mass of the disk and distance to the star. The volatile species incorporated in planets are mainly composed of H2O, CO, CO2, CH3OH, and NH3. Icy or ocean planets have systematically higher values of molecular abundances compared to giant and rocky planets. Gas giant planets are depleted in highly volatile molecules such as CH4, CO, and N2 compared to icy or ocean planets. The ice/rock mass ratio in icy or ocean and gas giant planets is, respectively, equal at maximum to 1.01 ± 0.33 and 0.8 ± 0.5, and is different from the usual assumptions made in planet formation models, which suggested this ratio to be 2–3. The C:O molar ratio in the atmosphere of gas giant planets is depleted by at least 30% compared to solar value.

From stellar nebula to planetesimals

Context. Solar and extrasolar comets and extrasolar planets are the subject of numerous studies in order to determine their chemical composition and internal structure. In the case of planetesimals, their compositions are important as they govern in part the composition of future planets. Aims. The present works aims at determining the chemical composition of icy planetesimals, believed to be similar to present day comets, formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data on chemical composition for models of planetesimals and comets, and models of planet formation and evolution. Methods. We have developed a model that calculates the composition of ices formed during the cooling of the stellar nebula. Coupled with a model of refractory element formation, it allows us to determine the chemical composition and mass ratio of ices to rocks in icy planetesimals throughout in the protoplanetary disc. Results. We provide relationships for ice line positions (for different volatile species) in the disc, and chemical compositions and mass ratios of ice relative to rock for icy planetesimals in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, a wide variety of chemical compositions of planetesimals were produced as a function of the mass of the disc and distance to the star. Ices incorporated in planetesimals are mainly composed of H2O, CO, CO2, CH3OH, and NH3. The ice/rock mass ratio is equal to 1 ± 0.5 in icy planetesimals following assumptions. This last value is in good agreement with observations of solar system comets, but remains lower than usual assumptions made in planet formation models, taking this ratio to be of 2–3.

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.

Composition-dependent outgassing of comet 67P/Churyumov-Gerasimenko from ROSINA/DFMS

Context. Early measurements of Rosetta’s target comet, 67P/Churyumov-Gerasimenko (67P), showed a strongly heterogeneous coma in H2O, CO, and CO2. Aims. The purpose of this work is to further investigate the coma heterogeneity of 67P, and to provide predictions for the near-perihelion outgassing profile based on the proposed explanations. Methods. Measurements of various minor volatile species by ROSINA/DFMS on board Rosetta are examined. The analysis focuses on the currently poorly illuminated winter (southern) hemisphere of 67P. Results. Coma heterogeneity is not limited to the major outgassing species. Minor species show better correlation with either H2O or CO2. The molecule CH4 shows a different diurnal pattern from all other analyzed species. Such features have implications for nucleus heterogeneity and thermal processing. Conclusions. Future analysis of additional volatiles and modeling the heterogeneity are required to better understand the observed coma profile.

Geomorphology of the Imhotep region on comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Context. Since August 2014, the OSIRIS Narrow Angle Camera (NAC) onboard the Rosetta spacecraft has acquired high spatial resolution images of the nucleus of comet 67P/Churyumov-Gerasimenko, down to the decimeter scale. This paper focuses on the Imhotep region, located on the largest lobe of the nucleus, near the equator. Aims. We map, inventory, and describe the geomorphology of the Imhotep region. We propose and discuss some processes to explain the formation and ongoing evolution of this region. Methods. We used OSIRIS NAC images, gravitational heights and slopes, and digital terrain models to map and measure the morphologies of Imhotep. Results. The Imhotep region presents a wide variety of terrains and morphologies: smooth and rocky terrains, bright areas, linear features, roundish features, and boulders. Gravity processes such as mass wasting and collapse play a significant role in the geomorphological evolution of this region. Cometary processes initiate erosion and are responsible for the formation of degassing conduits that are revealed by elevated roundish features on the surface. We also propose a scenario for the formation and evolution of the Imhotep region; this implies the presence of large primordial voids inside the nucleus, resulting from its formation process.