Uranus and Neptune are commonly considered ice giants, and it is often assumed that, in addition to a solar mix of hydrogen and helium, they contain roughly twice as much water as rock. This classical picture has led to successful models of their internal structure and has been understood to be compatible with the composition of the solar nebula during their formation (Reynolds and Summers 1965; Podolak and Cameron 1974; Podolak and Reynolds 1984; Podolak et al. 1995; Nettelmann et al. 2013). However, the dominance of water has been recently questioned (Teanby et al. 2020; Helled and Fortney 2020; Podolak et al. 2022). Planetesimals in the outer solar system are composed mainly of refractory materials, leading to an inconsistency between the icy composition of Uranus and Neptune and the ice-poor planetesimals they accreted during formation (Podolak et al. 2022). Here we elaborate on this problem, and propose a new potential solution. We show that chemical reactions between planetesimals dominated by organic-rich refractory materials and the hydrogen in gaseous atmospheres of protoplanets can form large amounts of methane 'ice'. Uranus and Neptune could thus be compatible with having accreted refractory-dominated planetesimals, while still remaining icy. Using random statistical computer models for a wide parameter space, we show that the resulting methane-rich internal composition could be a natural solution, giving a good match to the size, mass and moment of inertia of Uranus and Neptune, whereas rock-rich models appear to only work if a rocky interior is heavily mixed with hydrogen. Our model predicts a lower than solar hydrogen to helium ratio, which can be tested. We conclude that Uranus, Neptune and similar exoplanets could be methane-rich, and discuss why Jupiter and Saturn cannot.

Continuing our initiative on advancing the calculations of planetesimal accretion in the core-accretion model, we present here the results of our recent study of the contributions of planetesimals around and beyond the orbit of Saturn. In our first two papers, where our focus was on the effects of the Sun and Saturn, the initial distribution of planetesimals was limited to the regions around the accretion zone of a growing Jupiter. In this paper, we expanded that distribution to regions beyond the accretion zone of Saturn. We integrated the orbits of a large ensemble of planetesimals and studied the rate of their capture by the growing proto-Jupiter. In order to be consistent with our previous studies, we did not consider the effect of the nebular gas. Results demonstrated that the exterior planetesimals, especially those beyond 8 au, have only slight contributions to the growth and metallicity of the growing Jupiter. The final mass and composition of this planet is mainly due to the planetesimals inside and around its accretion zone. Our study shows that although the rate of capture varies slightly by the size and composition of planetesimals, in general, the final results are independent of the size and material of these bodies. Results also pointed to a new finding: the rate of accretion follows the same trend as that of the evolution of Jupiter’s envelope, with the largest accretion occurring during the envelope’s collapse. We present details of our analysis and discuss the implications of its results.

Although carbon monoxide (CO) is an abundant molecule and may have great importance for planetary interiors, measurements of its properties are difficult due to its extreme volatility. We calculate the equation of state for CO over a range of temperature and density that is applicable to the conditions in planetary interiors. Previous experimental and theoretical studies cover only a limited temperature-density range. Our calculations match these early results well, but now cover the full range of relevance. The method of calculation is based on the general-purpose quotidian equation of state described by More et al. (1988), which is here used in order to generate a freely downloadable look-up table to be used by the community.

We consider planets composed of water ice and rock, located far from a central star. In an earlier study, computing the growth of planets by continuous accretion, we found that a large fraction of the ice evaporates upon accretion, creating a water vapor atmosphere. Here we consider accretion as a discrete series of planetesimal impacts (of order 10^8), at the same time-dependent accretion rate, and investigate the fate of the vapor, as a result of its interaction with the accreting planetesimals. We find that a large fraction of the vapor escapes. The remaining fraction may form an outer layer of ice after the termination of accretion and cooling of the surface. The escaped water mass may significantly alter the ice-to-rock ratio of the planet. We investigate the effect of different choices of parameters such as the ice-to-rock ratio, the planetesimal size distribution, and the impact velocities. We find that the planetesimal size distribution has a negligible effect and explain why. By contrast, the ice-to-rock ratio and impact velocities affect the fraction of retained water masses considerably.

It is reasonable to assume that the structure of a planet and the interior distribution of its components are determined by its formation history. We thus follow the growth of a planet from a small embryo through its subsequent evolution. We estimate the accretion rate range based on a protoplanetary disk model at a large-enough distance from the central star for water ice to be a major component. We assume the accreted material to be a mixture of silicate rock and ice, with no H–He envelope, as the accretion timescale is much longer than the time required for the nebular gas to dissipate. We adopt a thermal evolution model that includes accretional heating, radioactive energy release, and separation of ice and rock. Taking the Safronov parameter and the ice-to-rock ratio as free parameters, we compute growth and evolutionary sequences for different parameter combinations, for 4.6 Gyr. We find the final structure to depend significantly on both parameters. Low initial ice-to-rock ratios and high accretion rates, each resulting in an increased heating rate, lead to the formation of extended rocky cores, while the opposite conditions leave the composition almost unchanged and result in relatively low internal temperatures. When rocky cores form, the ice-rich outer mantles still contain rock mixed with the ice. We find that a considerable fraction of the ice evaporates upon accretion, depending on parameters, and assume it is lost, thus the final surface composition and bulk density of the planet do not necessarily reflect the protoplanetary disk composition.

As part of our ongoing initiative on accurately calculating the accretion rate of planetesimals in the core-accretion model, we demonstrated in a recent article that when the calculations include the gravitational force of the Sun (the original core-accretion model did not include solar gravity), results change considerably [ApJ, 899:45]. In this paper, we have advanced our previous study by including the effect of Saturn. To maintain focus on the effect of this planet, and in order to be consistent with previous studies, we did not include the effect of the nebular gas. Results demonstrated that as expected, Saturn's perturbation decreases the rate of accretion by scattering many planetesimals out of Jupiter's accretion zone. It also increases the velocities with which planetesimals encounter the envelope, which in agreement with our previous findings, enhances their break-up due to the ram-pressure. Results also show that, because the effect of Saturn in scattering of planetesimals increases with its mass, this planet might not have played a significant role in the accretion of planetesimals by proto-Jupiter during the early stage of its growth. Finally, the late accretion of planetesimals, as obtained in our previous study, appears in our new results as well, implying that combined with the rapid in-fall of the gas, it can result in the mixing of material in the outer regions of the envelope which may explain the enhancement of the envelope's high-Z material.

The interest in the structure of ice-rich planetary bodies, in particular the differentiation between ice and rock, has grown due to the discovery of Kuiper Belt objects and exoplanets. We thus carry out a parameter study for a range of planetary masses M , yielding radii 50 ≲ R ≲ 3000 km, and for rock to ice mass ratios between 0.25 and 4, evolving them for 4.5 Gyr in a cold environment, to obtain the present structure. We use a thermal evolution model that allows for liquid and vapor flow in a porous medium, solving mass and energy conservation equations under hydrostatic equilibrium for a spherical body in orbit around a central star. The model includes the effect of pressure on porosity and on the melting temperature, heating by long-lived radioactive isotopes, and temperature-dependent serpentinization and dehydration. We obtain the boundary in parameter space (size, rock content) between bodies that differentiate, forming a rocky core, and those which remain undifferentiated: small bodies, bodies with a low rock content, and the largest bodies considered, which develop high internal pressures and barely attain the melting temperature. The final differentiated structure comprises a rocky core, an ice-rich mantle, and a thin dense crust below the surface. We obtain and discuss the bulk density–radius relationship. The effect of a very cold environment is investigated, and we find that at an ambient temperature of ∼20 K, small bodies preserve the ice in amorphous form to the present.

Modeling the interior of a planet is difficult because the small number of measured parameters is insufficient to constrain the many variables involved in describing the interior structure and composition. One solution is to invoke additional constraints based on arguments about how the planet formed. However, a planet’s actual structure and composition may hold clues to its formation which would be lost if this structure were not allowed by the initial assumptions. It is therefore interesting to explore the space of allowable compositions and structures in order to better understand which cosmogonic constraints are absolutely necessary. To this end, we describe a code for generating random, monotonic, density distributions, ρ(r), that fit a given mass, radius, and moment of inertia. Integrating the equation of hydrostatic equilibrium gives the pressure, P (r), at each point in the body. We then provide three algorithms for generating a monotonic temperature distribution, T (r), and an associated composition that is consistent with the ρ−P relation, and realistic equations of state. We apply this code to Uranus as a proof of concept, and show that the ratio of rock to water cannot be much larger than 2.

Solar-type young stellar objects undergo periodic, energetic outbursts that appear to be the result of enhanced mass accretion driven by the gravitational instability of their disks. Such FU Orionis outbursts may have profound consequences for the earliest solids in a protoplanetary disk, namely the refractory inclusions containing abundant calcium and aluminum (CAIs). We present models of the orbital evolution of centimeter-radius particles representing large CAIs in marginally gravitationally unstable disks. The hydrodynamical evolution of the disks is calculated with a fully three-dimensional code, including compressional heating and cooling in the beta cooling approximation. The particles are initially distributed uniformly throughout the disk, which extends from 1 to 10 au around a solar-mass protostar, but within ∼100 yr the particles are concentrated by gas drag into regions surrounding the spiral arms and rings formed by the gas disk. The particles settle down toward the disk midplane, only to be lofted repeatedly upward by shock fronts. Large-scale radial transport both outward and inward occurs, with significant numbers of particles reaching the outer disk (∼10 au) and surviving for considerably longer times than would be the case in a quiescent disk with gas pressure monotonically decreasing with distance from the protostar. Individual particles experience wide ranges of disk temperatures during their journeys, ranging from 60 K in the outer disk to nearly 2000 K in spiral features. Future work will consider the implications for CAI rims of the thermochemical processing experienced during FU Orionis outbursts.

It has been a common assumption of interior models that the outer planets of our solar system are convective, and that the internal temperature distributions are therefore adiabatic. This assumption is also often applied to exoplanets. However, if a large portion of the thermal flux can be transferred by conduction, or if convection is inhibited, the thermal profile could be substantially different and would therefore affect the inferred planetary composition. Here we investigate how the assumption of non-adiabatic temperature profiles in Uranus and Neptune affects their internal structures and compositions. We use a set of plausible temperature profiles together with density profiles that match the measured gravitational fields to derive the planets' compositions. We find that the inferred compositions of both Uranus and Neptune are quite sensitive to the assumed thermal profile in the outer layers, but relatively insensitive to the thermal profile in the central, high pressure region. The overall value of the heavy element mass fraction, Z, for these planets is between 0.8 and 0.9. Finally, we suggest that large parts of Uranus' interior might be conductive, a conclusion that is consistent with Uranus dynamo models and a hot central inner region.

On 12 November 2014, during the descent of the Rosetta lander Philae to the surface of comet 67P/Churyumov-Gerasimenlco the Dust Impact Monitor (DIM) on board Philae recorded an impact of a cometary dust impact of a cometary dust particle at 2.4 km from the comet surface (5 km from the nucleus' barycentre). In this work, we report further experiments that support the identification of this particle. We use aerogel as a comet analog material to characterise the properties of this particle. Our experiments show that this particle has a radius of 0.9 mm, a low density of 0.25 g/cm(3) and a high porosity close to 90%. The particle likely moved at near 4 m/s with respect to the comet. (C) 2017 Elsevier Inc. All rights reserved.

Abstract Modern observational techniques are still not powerful enough to directly view planet formation, and so it is necessary to rely on theory. However, observations do give two important clues to the formation process. The first is that the most primitive form of material in interstellar space exists as a dilute gas. Some of this gas is unstable against gravitational collapse, and begins to contract. Because the angular momentum of the gas is not zero, it contracts along the spin axis, but remains extended in the plane perpendicular to that axis, so that a disk is formed. Viscous processes in the disk carry most of the mass into the center where a star eventually forms. In the process, almost as a by-product, a planetary system is formed as well. The second clue is the time required. Young stars are indeed observed to have gas disks, composed mostly of hydrogen and helium, surrounding them, and observations tell us that these disks dissipate after about 5 to 10 million years. If planets like Jupiter and Saturn, which are very rich in hydrogen and helium, are to form in such a disk, they must accrete their gas within 5 million years of the time of the formation of the disk. Any formation scenario one proposes must produce Jupiter in that time, although the terrestrial planets, which don’t contain significant amounts of hydrogen and helium, could have taken longer to build. Modern estimates for the formation time of the Earth are of the order of 100 million years. To date there are two main candidate theories for producing Jupiter-like planets. The core accretion (CA) scenario supposes that any solid materials in the disk slowly coagulate into protoplanetary cores with progressively larger masses. If the core remains small enough it won’t have a strong enough gravitational force to attract gas from the surrounding disk, and the result will be a terrestrial planet. If the core grows large enough (of the order of ten Earth masses), and the disk has not yet dissipated, then the planetary embryo can attract gas from the surrounding disk and grow to be a gas giant. If the disk dissipates before the process is complete, the result will be an object like Uranus or Neptune, which has a small, but significant, complement of hydrogen and helium. The main question is whether the protoplanetary core can grow large enough before the disk dissipates. A second scenario is the disk instability (DI) scenario. This scenario posits that the disk itself is unstable and tends to develop regions of higher than normal density. Such regions collapse under their own gravity to form Jupiter-mass protoplanets. In the DI scenario a Jupiter-mass clump of gas can form—in several hundred years which will eventually contract into a gas giant planet. The difficulty here is to bring the disk to a condition where such instabilities will form. Now that we have discovered nearly 3000 planetary systems, there will be numerous examples against which to test these scenarios.

We examine the possibility that icy super-Earth mass planets, formed over long timescales (0.1-1 Gyr) at large distances (similar to 200-1000 au) from their host stars, will develop massive H-rich atmospheres. Within the interior of these planets, high pressure converts CH4 into ethane, butane, or diamond and releases H-2. Using simplified models that capture the basic physics of the internal structure, we show that the physical properties of the atmosphere depend on the outflux of H-2 from the mantle. When this outflux is less than or similar to 10(10) molec cm(-2) s(-1), the outgassed atmosphere has a. base pressure of less than or similar to 1 bar. Larger outflows result in a substantial atmosphere where the base pressure may approach 10(3) - 10(4) bar. For any pressure, the mean density of these planets, 2.4-3. g. cm(-3), is much larger than the mean density of Uranus and Neptune, 1.3-1.6. g. cm(-3). Thus, observations can distinguish between a Planet Nine with a primordial H/He-rich atmosphere accreted from the protosolar nebula and one with an atmosphere outgassed from the core.

We consider super-Earth sized planets which have a water mass fraction large enough to form an external mantle composed of high-pressure water-ice polymorphs and also lack a substantial H/He atmosphere. We consider such planets in their habitable zone, so that their outermost condensed mantle is a global, deep, liquid ocean. For these ocean planets, we investigate potential internal reservoirs of CO2, the amount of CO2 dissolved in the ocean for the various saturation conditions encountered, and the ocean-atmosphere exchange flux of CO2. We find that, in a steady state, the abundance of CO2 in the atmosphere has two possible states. When wind-driven circulation is the dominant CO2 exchange mechanism, an atmosphere of tens of bars of CO2 results, where the exact value depends on the subtropical ocean surface temperature and the deep ocean temperature. When sea-ice formation, acting on these planets as a CO2 deposition mechanism, is the dominant exchange mechanism, an atmosphere of a few bars of CO2 is established. The exact value depends on the subpolar surface temperature. Our results suggest the possibility of a negative feedback mechanism, unique to water planets, where a reduction in the subpolar temperature drives more CO2 into the atmosphere to increase the greenhouse effect.

We consider super-Earth sized planets which have a water mass fraction that is large enough to form an external mantle composed of high pressure water ice polymorphs and that lack a substantial H/He atmosphere. We consider such planets in their habitable zone so that their outermost condensed mantle is a global deep liquid ocean. For these ocean planets we investigate potential internal reservoirs of CO2; the amount of CO2 dissolved in the ocean for the various saturation conditions encountered, and the ocean-atmosphere exchange flux of CO2. We find that in steady state the abundance of CO2 in the atmosphere has two possible states. When the wind-driven circulation is the dominant CO2 exchange mechanism, an atmosphere of tens of bars of CO2 results, where the exact value depends on the subtropical ocean surface temperature and the deep ocean temperature. When sea-ice formation, acting on these planets as a CO2 deposition mechanism, is the dominant exchange mechanism, an atmosphere of a few bars of CO2 is established. The exact value depends on the subpolar surface temperature. Our results suggest the possibility of a negative feedback mechanism, unique to water planets, where a reduction in the subpolar temperature drives more CO2 into the atmosphere to increase the greenhouse effect.

The internal structure of gas giant planets may be more complex than the commonly assumed core-envelope structure with an adiabatic temperature profile. Different primordial internal structures as well as various physical processes can lead to non-homogenous compositional distributions. A non-homogenous internal structure has a significant impact on the thermal evolution and final structure of the planets. In this paper, we present alternative structure and evolution models for Jupiter and Saturn allowing for non-adiabatic primordial structures and the mixing of heavy elements by convection as these planets evolve. We present the evolution of the planets accounting for various initial composition gradients, and in the case of Saturn, include the formation of a helium-rich region as a result of helium rain. We investigate the stability of regions with composition gradients against convection, and find that the helium shell in Saturn remains stable and does not mix with the rest of the envelope. In other cases, convection mixes the planetary interior despite the existence of compositional gradients, leading to the enrichment of the envelope with heavy elements. We show that non-adiabatic structures (and cooling histories) for both Jupiter and Saturn are feasible. The interior temperatures in that case are much higher than those for standard adiabatic models. We conclude that the internal structure is directly linked to the formation and evolution history of the planet. These alternative internal structures of Jupiter and Saturn should be considered when interpreting the upcoming Juno and Cassini data.

We present a simple model for gas and dust flow from 67P/Churyumov-Gerasimenko that can be used to understand the grain impact observed by the DIM instrument on Philae (Krueger et al. 2015) We show how model results when applied to the GIADA measurements (Rotundi et al. 2015; Della Corte et al. 2015) can be used, in conjunction with the results found by the MIRO (Schloerb et al. 2015) and VIRTIS (De Sanctis et al. 2015) instruments to infer surface properties such as surface temperature and surface ice fraction.

The internal structure of gas giant planets may be more complex than the commonly assumed core-envelope structure with an adiabatic temperature profile. Different primordial internal structures as well as various physical processes can lead to non-homogenous compositional distributions. A non-homogenous internal structure has a significant impact on the thermal evolution and final structure of the planets. Here we present alternative structure and evolution models for Jupiter and Saturn allowing for both adiabatic and non-adiabatic evolution. In convective regions we calculate the mixing of heavy elements by convection, as these planets evolve. We present the thermal and structural evolution of the planets accounting for various initial composition gradients, and in the case of Saturn, include the formation of a helium-rich region as a result of helium rain. We investigate the stability of the structure against convection, and find that the helium shell in Saturn remains stable and does not mix with the rest of the envelope. In other cases, convection mixes the planetary interior despite the existence of compositional gradients, leading to enrichment of the envelope with metals. We show that non-adiabatic structures (and cooling histories) for both Jupiter and Saturn are feasible, and the interior temperatures in that case are much higher that for standard adiabatic models. Moreover, we show that non-adiabatic evolution can suggest more than one mechanism to explain the current structures, including Saturn's high luminosity. We conclude that the internal structure is directly linked to the formation and evolution history of the planet. These alternative internal structures of Jupiter and Saturn should be considered when interpreting the upcoming Juno and Cassini data.

Knowledge of both the mass and radius of an exoplanet allows us to estimate its mean density, and therefore its composition. Exoplanets seem to fill a very large parameter space in terms of mass and composition, and unlike the solar-system's planets, exoplanets also have intermediate masses (similar to 5-50 M-circle plus) with various densities. In this Letter, we investigate the behavior of the mass-radius relation for methane (CH4) planets and show that when methane planets are massive enough (M-p greater than or similar to 15 M-circle plus), the methane can dissociate and lead to a differentiated planet with a carbon core, a methane envelope, and a hydrogen atmosphere. The contribution of a rocky core to the behavior of a CH4 planet is also considered. We also develop interior models for several detected intermediate-mass planets that could, in principle, be methane/methane-rich planets. The example of methane planets emphasizes the complexity of the mass-radius relation and the challenge involved in uniquely inferring the planetary composition.

The Rosetta lander Philae successfully landed on the nucleus of comet 67P/Churyumov-Gerasimenko on 12 November 2014. Philae carries the Dust Impact Monitor (DIM) on board, which is part of the Surface Electric Sounding and Acoustic Monitoring Experiment (SESAME). DIM employs piezoelectric PZT sensors to detect impacts by sub-millimeter and millimeter-sized ice and dust particles that are emitted from the nucleus and transported into the cometary coma. The DIM sensor measures dynamical data like flux and the directionality of the impacting particles. Mass and speed of the particles can be constrained assuming density and elastic particle properties. DIM was operated during three mission phases of Philae at the comet: (1) Before Philae's separation from Rosetta at distances of about 9.6 km, 11.8 km, and 25.3 km from the nucleus barycenter. In this mission phase particles released from the nucleus on radial trajectories remained undetectable because of significant obscuration by the structures of Rosetta, and no dust particles were indeed detected. (2) During Philae's descent to its nominal landing site Agilkia, DIM detected one approximately millimeter-sized particle at a distance of 5.0 km from the nucleus' barycenter, corresponding to an altitude of 2.4 km from the surface. This is the closest ever dust detection at a cometary nucleus by a dedicated in-situ dust detector. (3) At Philae's final landing site, Abydos, DIM detected no dust impact which may be due to low cometary activity in the vicinity of Philae, or due to shading by obstacles close to Philae, or both. Laboratory calibration experiments showed that the material properties of the detected particle are compatible with a porous particle having a bulk density of approximately $250\, \mathrm{kg\,m^{-3}}$. The particle could have been lifted off from the comet's surface by sublimating water ice.