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.

Gas composition of the main volatile elements in protoplanetary discs and its implication for planet formation

Context. Direct observations of gaseous exoplanets reveal that their gas envelope has a higher C/O ratio than that of the host star (e.g., Wasp 12-b). This has been explained by considering that the gas phase of the disc could be inhomogeneous, exceeding the stellar C/O ratio in regions where these planets formed; but few studies have considered the drift of the gas and planet migration. Aims. We aim to derive the gas composition in planets through planet formation to evaluate if the formation of giant planets with an enriched C/O ratio is possible. The study focusses on the effects of different processes on the C/O ratio, such as the disc evolution, the drift of gas, and planet migration. Methods. We used our previous models for computing the chemical composition, together with a planet formation model, to which we added the composition and drift of the gas phase of the disc, which is composed of the main volatile species H2O, CO, CO2, NH3, N2, CH3OH, CH4, and H2S, H2 and He. The study focusses on the region where ice lines are present and influence the C/O ratio of the planets. Results. Modelling shows that the condensation of volatile species as a function of radial distance allows for C/O enrichment in specific parts of the protoplanetary disc of up to four times the solar value. This leads to the formation of planets that can be enriched in C/O in their envelope up to three times the solar value. Planet migration, gas phase evolution and disc irradiation enables the evolution of the initial C/O ratio that decreases in the outer part of the disc and increases in the inner part of the disc. The total C/O ratio of the planets is governed by the contribution of ices accreted, suggesting that high C/O ratios measured in planetary atmospheres are indicative of a lack of exchange of material between the core of a planet and its envelope or an observational bias. It also suggests that the observed C/O ratio is not representative of the total C/O ratio of the planet.

The influence of absorbed solar radiation by Saharan dust on hurricane genesis

To date, the radiative impact of dust and the Sahar an air layer (SAL) on North Atlantic hurricane activity is not yet known. According to previous studies, dust stabilizes the atmosphere due to absorption of solar radiation but thus shifts convection to regions more conducive for hurricane genesis. Here we analyze differences in hurricane genesis and frequency from ensemble sensitivity simulations with radiatively active and inactive dust in the aerosol-climate model ECHAM6-HAM. We investigate dust burden and other hurricane-related variables and determine their influence on disturbances which develop into hurricanes (developing disturbances, DDs) and those which do not (nondeveloping disturbances, NDDs). Dust and the SAL are found to potentially have both inhibiting and supporting influences on background conditions for hurricane genesis. A slight southward shift of DDs is determined when dust is active as well as a significant warming of the SAL, which leads to a strengthening of the vertical circulation associated with the SAL. The dust burden of DDs is smaller in active dust simulations compared to DDs in simulations with inactive dust, while NDDs contain more dust in active dust simulations. However, no significant influence of radiatively active dust on other variables in DDs and NDDs is found. Furthermore, no substantial change in the DD and NDD frequency due to the radiative effects of dust can be detected.

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.

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.

Evaluation of a Novel Screw Position in a Type III Distal Phalanx Fracture Model: An Ex Vivo Study

OBJECTIVE: Mechanical evaluation of a novel screw position used for repair in a type III distal phalanx fracture model and assessment of solar canal penetration (SCP). STUDY DESIGN: Experimental study. SAMPLE POPULATION: Disarticulated equine hooves (n = 24) and 24 isolated distal phalanges. METHODS: Hooves/distal phalanges cut in a sagittal plane were repaired with 1 of 2 different cortical screw placements in lag fashion. In group 1 (conventional screw placement), the screw was inserted halfway between the proximal border of the solar canal (SC) and the subchondral bone surface on a line parallel to the dorsal cortex, whereas in group 2, the screw was inserted more palmar/plantar, where a perpendicular line drawn from the group 1 position reached the palmar/plantar cortex. Construct strength was evaluated by 3-point bending to failure. SCP was assessed by CT imaging and macroscopically. RESULTS: Screws were significantly longer in group 2 and in forelimbs. Group 2 isolated distal phalanges had a significantly more rigid fixation compared with the conventional screw position (maximum point at failure 31%, bending stiffness 41% higher). Lumen reduction of the SC was observed in 13/52 specimens (all from group 2), of which 9 were forelimbs. CONCLUSIONS: More distal screw positioning compared with the conventionally recommended screw position for internal fixation of type III distal phalangeal fractures allows placement of a longer screw and renders a more rigid fracture fixation. The novel screw position, however, carries a higher risk of SCP

Dynamics of dissipative Bose-Einstein condensation

We resolve the real-time dynamics of a purely dissipative s=1/2 quantum spin or, equivalently, hard-core boson model on a hypercubic d-dimensional lattice. The considered quantum dissipative process drives the system to a totally symmetric macroscopic superposition in each of the S3 sectors. Different characteristic time scales are identified for the dynamics and we determine their finite-size scaling. We introduce the concept of cumulative entanglement distribution to quantify multiparticle entanglement and show that the considered protocol serves as an efficient method to prepare a macroscopically entangled Bose-Einstein condensate.

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.

3D-modeling of Mercury's solar wind sputtered surface-exosphere environment

The efficiency of sputtered refractory elements by H+ and He++ solar wind ions from Mercury's surface and their contribution to the exosphere are studied for various solar wind conditions. A 3D solar wind-planetary interaction hybrid model is used for the evaluation of precipitation maps of the sputter agents on Mercury's surface. By assuming a global mineralogical surface composition, the related sputter yields are calculated by means of the 2013 SRIM code and are coupled with a 3D exosphere model. Because of Mercury's magnetic field, for quiet and nominal solar wind conditions the plasma can only precipitate around the polar areas, while for extreme solar events (fast solar wind, coronal mass ejections, interplanetary magnetic clouds) the solar wind plasma has access to the entire dayside. In that case the release of particles form the planet's surface can result in an exosphere density increase of more than one order of magnitude. The corresponding escape rates are also about an order of magnitude higher. Moreover, the amount of He++ ions in the precipitating solar plasma flow enhances also the release of sputtered elements from the surface in the exosphere. A comparison of our model results with MESSENGER observations of sputtered Mg and Ca elements in the exosphere shows a reasonable quantitative agreement. (C) 2015 Elsevier Ltd. All rights reserved.

Origin and history of ureilitic material in the solar system: The view from asteroid 2008 TC3 and the Almahata Sitta meteorite

Asteroid 2008 TC3 (approximately 4m diameter) was tracked and studied in space for approximately 19h before it impacted Earth's atmosphere, shattering at 44-36km altitude. The recovered samples (>680 individual rocks) comprise the meteorite Almahata Sitta (AhS). Approximately 50-70% of these are ureilites (ultramafic achondrites). The rest are chondrites, mainly enstatite, ordinary, and Rumuruti types. The goal of this work is to understand how fragments of so many different types of parent bodies became mixed in the same asteroid. Almahata Sitta has been classified as a polymict ureilite with an anomalously high component of foreign clasts. However, we calculate that the mass of fallen material was 0.1% of the pre-atmospheric mass of the asteroid. Based on published data for the reflectance spectrum of the asteroid and laboratory spectra of the samples, we infer that the lost material was mostly ureilitic. Therefore, 2008 TC3 probably contained only a few percent nonureilitic materials, similar to other polymict ureilites except less well consolidated. From available data for the AhS meteorite fragments, we conclude that 2008 TC3 samples essentially the same range of types of ureilitic and nonureilitic materials as other polymict ureilites. We therefore suggest that the immediate parent of 2008 TC3 was the immediate parent of all ureilitic material sampled on Earth. We trace critical stages in the evolution of that material through solar system history. Based on various types of new modeling and re-evaluation of published data, we propose the following scenario. (1) The ureilite parent body (UPB) accreted 0.5-0.6Ma after formation of calcium-aluminum-rich inclusions (CAI), beyond the ice line (outer asteroid belt). Differentiation began approximately 1Ma after CAI. (2) The UPB was catastrophically disrupted by a major impact approximately 5Ma after CAI, with selective subsets of the fragments reassembling into daughter bodies. (3) Either the UPB (before breakup), or one of its daughters (after breakup), migrated to the inner belt due to scattering by massive embryos. (4) One daughter (after forming in or migrating to the inner belt) became the parent of 2008 TC3. It developed a regolith, mostly 3.8Ga ago. Clasts of enstatite, ordinary, and Rumuruti-type chondrites were implanted by low-velocity collisions. (5) Recently, the daughter was disrupted. Fragments were injected or drifted into Earth-crossing orbits. 2008 TC3 comes from outer layers of regolith, other polymict ureilites from deeper regolith, and main group ureilites from the interior of this body. In contrast to other models that have been proposed, this model invokes a stochastic history to explain the unique diversity of foreign materials in 2008 TC3 and other polymict ureilites.

The isotope composition of selenium in chondrites constrains the depletion mechanism of volatile elements in solar system materials

Solar nebula processes led to a depletion of volatile elements in different chondrite groups when compared to the bulk chemical composition of the solar system deduced from the Sun's photosphere. For moderately-volatile elements, this depletion primarily correlates with the element condensation temperature and is possibly caused by incomplete condensation from a hot solar nebula, evaporative loss from the precursor dust, and/or inherited from the interstellar medium. Element concentrations and interelement ratios of volatile elements do not provide a clear picture about responsible mechanisms. Here, the abundance and stable isotope composition of the moderately- to highly-volatile element Se are investigated in carbonaceous, ordinary, and enstatite chondrites to constrain the mechanism responsible for the depletion of volatile elements in planetary bodies of the inner solar system and to define a δ(82/78)Se value for the bulk solar system. The δ(82/78)Se of the studied chondrite falls are identical within their measurement uncertainties with a mean of −0.20±0.26‰ (2 s.d., n=14n=14, relative to NIST SRM 3149) despite Se abundance depletions of up to a factor of 2.5 with respect to the CI group. The absence of resolvable Se isotope fractionation rules out a kinetic Rayleigh-type incomplete condensation of Se from the hot solar nebula or partial kinetic evaporative loss on the precursor material and/or the parent bodies. The Se depletion, if acquired during partial condensation or evaporative loss, therefore must have occurred under near equilibrium conditions to prevent measurable isotope fractionation. Alternatively, the depletion and cooling of the nebula could have occurred simultaneously due to the continuous removal of gas and fine particles by the solar wind accompanied by the quantitative condensation of elements from the pre-depleted gas. In this scenario the condensation of elements does not require equilibrium conditions to avoid isotope fractionation. The results further suggest that the processes causing the high variability of Se concentrations and depletions in ordinary and enstatite chondrites did not involve any measurable isotope fractionation. Different degrees of element depletions and isotope fractionations of the moderately-volatile elements Zn, S, and Se in ordinary and enstatite chondrites indicate that their volatility is controlled by the thermal stabilities of their host phases and not by the condensation temperature under canonical nebular conditions.

Critical core mass for enriched envelopes: the role of H2O condensation

Context. Within the core accretion scenario of planetary formation, most simulations performed so far always assume the accreting envelope to have a solar composition. From the study of meteorite showers on Earth and numerical simulations, we know that planetesimals must undergo thermal ablation and disruption when crossing a protoplanetary envelope. Thus, once the protoplanet has acquired an atmosphere, not all planetesimals reach the core intact, i.e. the primordial envelope (mainly H and He) gets enriched in volatiles and silicates from the planetesimals. This change of envelope composition during the formation can have a significant effect on the final atmospheric composition and on the formation timescale of giant planets. Aims. We investigate the physical implications of considering the envelope enrichment of protoplanets due to the disruption of icy planetesimals during their way to the core. Particular focus is placed on the effect on the critical core mass for envelopes where condensation of water can occur. Methods. Internal structure models are numerically solved with the implementation of updated opacities for all ranges of metallicities and the software Chemical Equilibrium with Applications to compute the equation of state. This package computes the chemical equilibrium for an arbitrary mixture of gases and allows the condensation of some species, including water. This means that the latent heat of phase transitions is consistently incorporated in the total energy budget. Results. The critical core mass is found to decrease significantly when an enriched envelope composition is considered in the internal structure equations. A particularly strong reduction of the critical core mass is obtained for planets whose envelope metallicity is larger than Z approximate to 0.45 when the outer boundary conditions are suitable for condensation of water to occur in the top layers of the atmosphere. We show that this effect is qualitatively preserved even when the atmosphere is out of chemical equilibrium. Conclusions. Our results indicate that the effect of water condensation in the envelope of protoplanets can severely affect the critical core mass, and should be considered in future studies.