The Science of Exoplanets and Their Systems

A scientific forum on “The Future Science of Exoplanets and Their Systems,” sponsored by Europlanet* and the International Space Science Institute (ISSI)† and co-organized by the Center for Space and Habitability (CSH)‡ of the University of Bern, was held during December 5 and 6, 2012, in Bern, Switzerland. It gathered 24 well-known specialists in exoplanetary, Solar System, and stellar science to discuss the future of the fast-expanding field of exoplanetary research, which now has nearly 1000 objects to analyze and compare and will develop even more quickly over the coming years. The forum discussions included a review of current observational knowledge, efforts for exoplanetary atmosphere characterization and their formation, water formation, atmospheric evolution, habitability aspects, and our understanding of how exoplanets interact with their stellar and galactic environment throughout their history. Several important and timely research areas of focus for further research efforts in the field were identified by the forum participants. These scientific topics are related to the origin and formation of water and its delivery to planetary bodies and the role of the disk in relation to planet formation, including constraints from observations as well as star-planet interaction processes and their consequences for atmosphere-magnetosphere environments, evolution, and habitability. The relevance of these research areas is outlined in this report, and possible themes for future ISSI workshops are identified that may be proposed by the international research community over the coming 2–3 years.

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.

The influence of recent major crater impacts on the surrounding surfaces of (21) Lutetia

We present 3-D simulations of impacts into Asteroid 21 Lutetia, the subject of a fly-by by the European Space Agency's Rosetta mission to Comet 67P/Churyumov-Gerasimenko. Using a 3-D shape model of the asteroid, impacts of sizes sufficient to reproduce the observed craters in Lutetia's North Polar Crater Cluster (NPCC) as observed by the OSIRIS experiment have been simulated using the Smoothed Particle Hydrodynamics technique. The asteroid itself has been modelled both as a homogeneous body and as a body with an iron core.

Micrometer-sized ice particles for planetary-science experiments – II. Bidirectional reflectance

We have measured the bidirectional reflectance of spherical micrometer-sized water-ice particles in the visible spectral range over a wide range of incidence and emission angles. The small ice spheres were produced by spraying fine water droplets directly into liquid nitrogen. The resulting mean particle radii are 1.47 + 0.96 - 0.58 μm. Such a material shares many properties with ice in comets and at the surface of icy satellites. Measurements show that the fresh sample material is highly backscattering, contrasting with natural terrestrial snow and frost. The formation of agglomerates of particles during the sample production results in a noticeable variability of the photometric properties of the samples in their initial state. We have also observed significant temporal evolutions of the scattering behavior of the samples, shifting towards more forward scattering within some tens of hours, resulting most likely from sintering processes. All reflectance data are fitted by the Hapke photometric model (1993 and 2002 formulation) with a one/two/three-parameter Henyey-Greenstein phase function and the resulting Hapke parameters are provided. These parameters can be used to compare laboratory results with the observed photometric behaviors of astronomical objects. We show, in particular, that the optical properties of the fresh micrometer-sized ice samples can be used to reproduce the predominant backscattering in the phase curves of Enceladus and Europa.

Observations of the northern seasonal polar cap on Mars II: HiRISE photometric analysis of evolution of northern polar dunes in spring

We present an overview of our analyses of HiRISE observations of spring evolution of selected dune areas of the north polar erg. The north polar erg is covered annually by seasonal volatile ice layer, a mixture of CO2 and H2O with mineral dust contamination. In spring, this layer sublimes creating visually enigmatic phenomena, e.g. dark and bright fan-shaped deposits, dark–bright–dark bandings, dark down-slope streaks, and seasonal polygonal cracks. Similar phenomena in southern polar areas are believed to be related to the specific process of solid-state greenhouse effect. In the north, it is currently unclear if the solid-state greenhouse effect is able to explain all the observed phenomena especially because the increased influence of H2O on the time scales of this process has not yet been quantified. HiRISE observations of our selected locations show that the ground exhibits a temporal behaviour similar to the one observed in the southern polar areas: a brightening phase starting close to the spring equinox with a subsequent darkening towards summer solstice. The resolution of HiRISE enabled us to study dunes and substrate individually and even distinguish between different developments on windward and slip face sides of single dunes. Differences in the seasonal evolution between steep slip faces and flatter substrate and windward sides of dunes have been identified and compared to CRISM data of CO2 and H2O distributions on dunes. We also observe small scale dark blotches that appear in early observations and tend to sustain a low reflectivity throughout the spring. These blotches can be regarded as the analogue of dark fan deposits in southern polar areas, leading us to the conclusion that both martian polar areas follow similar spring evolutions.

Observations of the northern seasonal polar cap on Mars III: CRISM/HiRISE observations of spring sublimation

We analyze a series of targeted CRISM and HiRISE observations of seven regions of interest at high latitudes in the Northern polar regions of Mars. These data allow us to investigate the temporal evolution of the composition of the seasonal ice cap during spring, with a special emphasis on peculiar phenomena occurring in the dune fields and in the vicinity of the scarps of the North Polar Layered Deposits (NPLDs). The strength of the spectral signature of CO2 ice continuously decreases during spring whereas the one of H2O ice first shows a strong increase until Ls = 50°. This evolution is consistent with a scenario previously established from analysis of OMEGA data, in which a thin layer of pure H2O ice progressively develops at the surface of the volatile layer. During early spring (Ls < 10°), widespread jet activity is observed by HiRISE while strong spectral signatures of CO2 ice are detected by CRISM. Later, around Ls = 20-40°, activity concentrates at the dune fields where CRISM also detects a spectral enrichment in CO2 ice, consistent with "Kieffer's model" (Kieffer, H.H. [2007]. J. Geophys. Res. 112, E08005. doi:10.1029/2006JE002816) for jet activity. Effects of wind are prominent across the dune fields and seem to strongly influence the sublimation of the volatile layer. Strong winds blowing down the scarps could also be responsible for the significant spatial and temporal variability of the surface ice composition observed close to the NPLD.

Coupling of LMS with a fs-laser ablation ion source: elemental and isotope composition measurements

Mass spectrometric analysis of elemental and isotopic compositions of several NIST standards is performed by a miniature laser ablation/ionisation reflectron-type time-of-flight mass spectrometer (LMS) using a fs-laser ablation ion source (775 nm, 190 fs, 1 kHz). The results of the mass spectrometric studies indicate that in a defined range of laser irradiance (fluence) and for a certain number of accumulations of single laser shot spectra, the measurements of isotope abundances can be conducted with a measurement accuracy at the per mill level and at the per cent level for isotope concentrations higher and lower than 100 ppm, respectively. Also the elemental analysis can be performed with a good accuracy. The LMS instrument combined with a fs-laser ablation ion source exhibits similar detection efficiency for both metallic and non-metallic elements. Relative sensitivity coefficients were determined and found to be close to one, which is of considerable importance for the development of standard-less instruments. Negligible thermal effects, sample damage and excellent characteristics of the fs-laser beam are thought to be the main reason for substantial improvement of the instrumental performance compared to other laser ablation mass spectrometers.

Performance evaluation of a miniature laser ablation time-of-flight mass spectrometer designed for in situ investigations in planetary space research

Key performance features of a miniature laser ablation time-of-flight mass spectrometer designed for in situ investigations of the chemical composition of planetary surfaces are presented. This mass spectrometer is well suited for elemental and isotopic analysis of raw solid materials with high sensitivity and high spatial resolution. In this study, ultraviolet laser radiation with irradiances suitable for ablation (< 1 GW/cm2) is used to achieve stable ion formation and low sample consumption. In comparison to our previous laser ablation studies at infrared wavelengths, several improvements to the experimental setup have been made, which allow accurate control over the experimental conditions and good reproducibility of measurements. Current performance evaluations indicate significant improvements to several instrumental figures of merit. Calibration of the mass scale is performed within a mass accuracy (Δm/m) in the range of 100 ppm, and a typical mass resolution (m/Δm) ~600 is achieved at the lead mass peaks. At lower laser irradiances, the mass resolution is better, about (m/Δm) ~900 for lead, and limited by the laser pulse duration of 3 ns. The effective dynamic range of the instrument was enhanced from about 6 decades determined in previous study up to more than 8 decades at present. Current studies show high sensitivity in detection of both metallic and non-metallic elements. Their abundance down to tens of ppb can be measured together with their isotopic patterns. Due to strict control of the experimental parameters, e.g. laser characteristics, ion-optical parameters and sample position, by computer control, measurements can be performed with high reproducibility. Copyright © 2012 John Wiley & Sons, Ltd.

Coupled W–Os–Pt isotope systematics in IVB iron meteorites: In situ neutron dosimetry for W isotope chronology

Tungsten isotope compositions of magmatic iron meteorites yield ages of differentiation that are within ±2 Ma of the formation of CAIs, with the exception of IVB irons that plot to systematically less radiogenic compositions yielding erroneously old ages. Secondary neutron capture due to galactic cosmic ray (GCR) irradiation is known to lower the ε182W of iron meteorites, adequate correction of which requires a measure of neutron dosage which has not been available, thus far. The W, Os and Pt isotope systematics of 12 of the 13 known IVB iron meteorites were determined by MC-ICP-MS (W, Os, Pt) and TIMS (Os). On the same dissolutions that yield precise ε182W, stable Os and Pt isotopes were determined as in situ neutron dosimeters for empirical correction of the ubiquitous cosmic-ray induced burn-out of 182W in iron meteorites. The W isotope data reveal a main cluster with ε182W of ∼−3.6, but a much larger range than observed in previous studies including irons (Weaver Mountains and Warburton Range) that show essentially no cosmogenic effect on their ε182W. The IVB data exhibits resolvable negative anomalies in ε189Os (−0.6ε) and complementary ε190Os anomalies (+0.4ε) in Tlacotepec due to neutron capture on 189Os which has approximately the same neutron capture cross section as 182W, and captures neutrons to produce 190Os. The least irradiated IVB iron, Warburton Range, has ε189Os and ε190Os identical to terrestrial values. Similarly, Pt isotopes, which are presented as ε192Pt, ε194Pt and ε196Pt range from +4.4ε to +53ε, +1.54ε to −0.32ε and +0.73ε to −0.20ε, respectively, also identify Tlacotepec and Dumont as the most GCR-damaged samples. In W–Os and W–Pt isotope space, the correlated isotope data back-project toward a 0-epsilon value of ε192Pt, ε189Os and ε190Os from which a pre-GCR irradiation ε182W of −3.42±0.09 (2σ) is derived. This pre-GCR irradiation ε182W is within uncertainty of the currently accepted CAI initial ε182W. The Pt and Os isotope correlations in the IVB irons are in good agreement with a nuclear model for spherical irons undergoing GCR spallation, although this model over-predicts the change of ε182W by ∼2×, indicating a need for better W neutron capture cross section determinations. A nucleosynthetic effect in ε184W in these irons of −0.14±0.08 is confirmed, consistent with the presence of Mo and Ru isotope anomalies in IVB irons. The lack of a non-GCR Os isotope anomaly in these irons requires more complex explanations for the production of W, Ru and Mo anomalies than nebular heterogeneity in the distribution of s-process to r-process nuclides.

Neutron capture on Pt isotopes in iron meteorites and the Hf–W chronology of core formation in planetesimals

The short-lived 182Hf–182W isotope system can provide powerful constraints on the timescales of planetary core formation, but its application to iron meteorites is hampered by neutron capture reactions on W isotopes resulting from exposure to galactic cosmic rays. Here we show that Pt isotopes in magmatic iron meteorites are also affected by capture of (epi)thermal neutrons and that the Pt isotope variations are correlated with variations in 182W/184W. This makes Pt isotopes a sensitive neutron dosimeter for correcting cosmic ray-induced W isotope shifts. The pre-exposure 182W/184W derived from the Pt–W isotope correlations of the IID, IVA and IVB iron meteorites are higher than most previous estimates and are more radiogenic than the initial 182W/184W of Ca–Al-rich inclusions (CAI). The Hf–W model ages for core formation range from +1.6±1.0 million years (Ma; for the IVA irons) to +2.7±1.3 Ma after CAI formation (for the IID irons), indicating that there was a time gap of at least ∼1 Ma between CAI formation and metal segregation in the parent bodies of some iron meteorites. From the Hf–W ages a time limit of <1.5–2 Ma after CAI formation can be inferred for the accretion of the IID, IVA and IVB iron meteorite parent bodies, consistent with earlier conclusions that the accretion of differentiated planetesimals predated that of most chondrite parent bodies.

Cosmogenic production rates and recoil loss effects in micrometeorites and interplanetary dust particles

We present a purely physical model to determine cosmogenic production rates for noble gases and radionuclides in micrometeorites (MMs) and interplanetary dust particles (IDPs) by solar cosmic-rays (SCR) and galactic cosmic-rays (GCR) fully considering recoil loss effects. Our model is based on various nuclear model codes to calculate recoil cross sections, recoil ranges, and finally the percentages of the cosmogenic nuclides that are lost as a function of grain size, chemical composition of the grain, and the spectral distribution of the projectiles. The main advantage of our new model compared with earlier approaches is that we consider the entire SCR particle spectrum up to 240 MeV and not only single energy points. Recoil losses for GCR-produced nuclides are assumed to be equal to recoil losses for SCR-produced nuclides. Combining the model predictions with Poynting-Robertson orbital lifetimes, we calculate cosmic-ray exposure ages for recently studied MMs, cosmic spherules, and IDPs. The ages for MMs and the cosmic-spherule are in the range <2.2–233 Ma, which corresponds, according to the Poynting-Robertson drag, to orbital distances in the range 4.0–34 AU. For two IDPs, we determine exposure ages of longer than 900 Ma, which corresponds to orbital distances larger than 150 AU. The orbital distance in the range 4–6 AU for one MM and the cosmic spherule indicate an origin either in the asteroid belt or release from comets coming either from the Kuiper Belt or the Oort Cloud. Three of the studied MMs have orbital distances in the range 23–34 AU, clearly indicating a cometary origin, either from short-period comets from the Kuiper Belt or from the Oort Cloud. The two IDPs have orbital distances of more than 150 AU, indicating an origin from Oort Cloud comets.