We present results of a detailed study of the rate of the accretion of planetesimals by a growing proto-Jupiter in the core-accretion model. Using a newly developed code, we accurately combine a detailed three-body trajectory calculation with gas drag experienced during the passage of planetesimals in the protoplanet's envelope. We find that the motion of planetesimals is excited to the extent that encounters with the proto-planetary envelope become so fast that ram pressure breaks up the planetesimals in most encounters. As a result, the accretion rate is largely independent of the planetesimal size and composition. For the case we explored of a planet forming at 5.2 AU from the Sun in a disk with a solid surface density of 6 g/cm^2 (Lozovsky et al. 2017) the accretion rate we compute differs in several respects from that assumed by those authors. We find that only 4-5 M_Earth is accreted in the first 1.5x10^6 years before the onset of rapid gas accretion. Most of the mass, some 10 M_Earth, is accreted simultaneously with this rapid gas accretion. In addition, we find that the mass accretion rate remains small, but non-zero for at least a million years after this point, and an additional 0.3-0.4 M_Earth is accreted during that time. This late accretion, together with a rapid infall of gas could lead to the accreted material being mixed throughout the outer regions, and may account for the enhancement of high-Z material in Jupiter's envelope.

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.

We examine the possibility that icy super-Earth mass planets, formed over long time scales (0.1--1~Gyr) at large distances ($\sim$ 200--1000~AU) from their host stars, will develop massive H-rich atmospheres. Within the interior of these planets, high pressure converts CH$_4$ into ethane, butane, or diamond and releases H$_2$. Using simplified models which 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 $\lesssim 10^{10}$ [molec cm$^{-2}$ s$^{-1}$], the outgassed atmosphere has base pressure $\lesssim$ 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 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 that 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 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.

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 (~ 5 - 50 M_Earth) with various densities. In this letter, we investigate the behavior of the Mass-Radius relation for methane (CH_4) planets and show that when methane planets are massive enough (Mp >~ 15 M_Earth), 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 CH_4 planet is considered as well. 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.

We study the evolution of planetesimals in evolved gaseous disks, which orbit a solar-mass star and harbor a Jupiter-mass planet at a_p~5AU. The gas dynamics is modeled with a three-dimensional hydrodynamics code that employes nested-grids and achieves a resolution of one Jupiter's radius in the circumplanetary disk. The code models solids as individual particles. Planetesimals are subjected to gravitational forces by the star and the planet, drag force by the gas, disruption via ram pressure, and mass loss through ablation. The mass evolution of solids is calculated self-consistently with their temperature, velocity, and position. We consider icy and icy/rocky bodies of radius 0.1-100km, initially deployed on orbits around the star within a few Hill radii (Rhill) of the planet's orbit. Planetesimals are scattered inward, outward, and toward disk regions of radius r>>a_p. Scattering can relocate significant amounts of solids, provided that regions |r-a_p|~ 3Rhill are replenished with planetesimals. Scattered bodies can be temporarily captured on planetocentric orbits. Ablation consumes nearly all solids at gas temperatures > ~220K. Super-keplerian rotation around and beyond the outer edge of the gas gap can segregate < ~0.1km bodies, producing solid gap edges at size-dependent radial locations. Capture, break-up, and ablation of solids result in a dust-laden circumplanetary disk with low surface densities of km-size planetesimals, implying relatively long timescales for satellite formation. After a giant planet acquires most of its mass, accretion of solids is unlikely to alter significantly its heavy-element content. The luminosity generated by solids' accretion can be of a similar order of magnitude to the contraction luminosity.

The primordial internal structures of gas giant planets are unknown. Often giant planets are modeled under the assumption that they are adiabatic, convective, and homogeneously mixed, but this is not necessarily correct. In this work, we present the first self-consistent calculation of convective transport of both heat and material as the planets evolve. We examine how planetary evolution depends on the initial composition and its distribution, whether the internal structure changes with time, and if so, how it affects the evolution. We consider various primordial distributions, different compositions, and different mixing efficiencies and follow the distribution of heavy elements in a Jupiter-mass planet as it evolves. We show that a heavy-element core cannot be eroded by convection if there is a sharp compositional change at the core-envelope boundary. If the heavy elements are initially distributed within the planet according to some compositional gradient, mixing occurs in the outer regions resulting in a compositionally homogeneous outer envelope. Mixing of heavy materials that are injected in a convective gaseous envelope are found to mix efficiently. Our work demonstrates that the primordial internal structure of a giant planet plays a substantial role in determining its long-term evolution and that giant planets can have non-adiabatic interiors. These results emphasize the importance of coupling formation, evolution, and internal structure models of giant planets self-consistently.

We study the transport of methane in the external water envelopes surrounding water-rich super-Earths and estimate its outgassing into the atmosphere. We investigate the influence of methane on the thermodynamics and mechanics of the water mantle. We find that including methane in the water matrix introduces a new phase (filled ice) resulting in hotter planetary interiors. This effect renders the super-ionic and reticulating phases accessible to relatively low mass planets lacking a H/He atmosphere. We model the thermal and structural profile of the planetary crust and discuss five possible crustal regimes. The formation of methane clathrate in the subsurface is shown to inhibit the formation of a subterranean ocean. This effect results in increased stresses on the lithosphere making modes of ice plate tectonics possible. The dynamics of the tectonic plates are analysed. We derive overturn and resurfacing time scales as well as the melt fraction underneath spreading centers. Ice mantle dynamics is found to be important for assessing the composition of the atmosphere. We formulate the relation between the outgassing flux of methane and the tectonic mode dynamics. We give numerical estimates for the global outgassing rate of methane into the atmosphere.

The large number of detected giant exoplanets offers the opportunity to improve our understanding of the formation mechanism, evolution, and interior structure of gas giant planets. The two main models for giant planet formation are core accretion and disk instability. There are substantial differences between these formation models, including formation timescale, favorable formation location, ideal disk properties for planetary formation, early evolution, planetary composition, etc. First, we summarize the two models including their substantial differences, advantages, and disadvantages, and suggest how theoretical models should be connected to available (and future) data. We next summarize current knowledge of the internal structures of solar- and extrasolar- giant planets. Finally, we suggest the next steps to be taken in giant planet exploration.

We model the evolution of planets with various masses and compositions. We investigate the effects of the composition and its depth dependence on the long-term evolution of the planets. The effects of opacity and stellar irradiation are also considered. It is shown that the change in radius due to various compositions can be significantly smaller than the change in radius caused by the opacity. Irradiation also affects the planetary contraction but is found to be less important than the opacity effects. We suggest that the mass-radius relationship used for characterization of observed extrasolar planets should be taken with great caution since different physical conditions can result in very different mass-radius relationships.

Whether volatiles can be entrapped in a background matrix composing planetary envelopes and be dragged via convection to the surface is a key question in understanding atmospheric fluxes, cycles and composition. In this paper we consider super-Earths with an extensive water mantle (i.e. water planets), and the possibility of entrapment of methane in their extensive water ice envelopes. We adopt the theory developed by van der Waals & Platteeuw (1959) for modelling solid solutions, often used for modelling clathrate hydrates, and modify it in order to estimate the thermodynamic stability field of a new phase, called methane filled ice Ih. We find that in comparison to water ice VII the filled ice Ih structure may be stable not only at the high pressures but also at the high temperatures expected at the core-water mantle transition boundary of water planets.

The internal structures and compositions of Uranus and Neptune are not well constrained due to the uncertainty in rotation period and flattening, as well as the relatively large error bars on the gravitational coefficients. While Uranus and Neptune are similar in mass and radius, they differ in other physical properties such as thermal emission, obliquity, and inferred atmospheric enrichment. In this letter we consider the uncertainty in the planetary rotation periods, show that rotation periods more consistent with the measured oblateness imply that Uranus and Neptune have different internal structures, and speculate on the source of that difference. We conclude that Uranus and Neptune might have very different structures and/or compositions despite their similar masses and radii. We point out that understanding these differences can have important implications for our view of the formation and evolution of Uranus and Neptune as well as intermediate-mass extra-solar planets in general.

The solar nebula is thought to have undergone a number of episodes of FU Orionis outbursts during its early evolution. We present here the first calculations of the trajectories of particles in a marginally gravitationally unstable solar nebula during an FU Orionis outburst, which show that 0.1 to 10 cm-sized particles traverse radial distances of 10 AU or more, inward and outward, in less than 200 yrs, exposing the particles to temperatures from $\sim$ 60 K to $\sim$ 1500 K. Such trajectories can thus account for the discovery of refractory particles in comets. Refractory particles should acquire Wark-Lovering-like rims as they leave the highest temperature regions of the disk, and these rims should have significant variations in their stable oxygen isotope ratios. Particles are likely to be heavily modified or destroyed if they pass within 1 AU of the Sun, and so are only likely to survive if they formed in the final few FU Orionis outbursts, or were transported to the outer reaches of the solar system. Calcium, aluminum-rich inclusions (CAIs) from primitive meteorites are the oldest known solar system objects and have a very narrow age range. Most CAIs may have formed at the end of the FU Orionis outbursts phase, with an age range reflecting the period between the last few outbursts.

By using a hydrodynamic atmospheric escape mechanism \citep{amit09} we show how the unusually high mass density of Quaoar could have been predicted (constrained), without any knowledge of a binary companion. We suggest an explanation of the recent spectroscopic observations of Orcus and Charon \citep{delsanti10,cook07}. We present a simple relation between the detection of certain volatile ices and the body mass density and diameter. As a test case we implement the relations on the KBO 2003 AZ$_{84}$ and give constraints on its mass density. We also present a method of relating the latitude-dependence of hydrodynamic gas escape to the internal structure of a rapidly rotating body and apply it to Haumea.

We investigate the evolution of grains composed of an ice shell surrounding an olivine core as they pass through a spiral shock in a protoplanetary disk. We use published three-dimensional radiation-hydrodynamics simulations of massive self-gravitating protoplanetary disks to extract the thermodynamics of spiral shocks in the region between 10 and 20 AU from the central star. As the density wave passes, it heats the grains, causing them to lose their ice shell and resulting in a lowering of the grain opacity. In addition, since grains of different sizes will have slightly different temperatures, there will be a migration of ice from the hotter grains to the cooler ones. The rate of migration depends on the temperature of the background gas, so a grain distribution that is effectively stable for low temperatures, can undergo an irreversible change in opacity if the gas is temporarily heated to above $\sim 150$\,K. We find that the opacity can drop more, and at a significantly faster rate throughout the spiral shocks relative to the prediction of standard dust grains models adopted in hydrodynamical calculations of protoplanetary disks. This would lead to faster gas cooling within spiral arms. We discuss the implications of our results on the susceptibility of disk fragmentation into sub-stellar objects at distances of a few tens of astronomical units.

'Empirical' models (pressure vs. density) of Uranus and Neptune interiors constrained by the gravitational coefficients J_2, J_4, the planetary radii and masses, and Voyager solid-body rotation periods are presented. The empirical pressure-density profiles are then interpreted in terms of physical equations of state of hydrogen, helium, ice (H_2O), and rock (SiO_2) to test the physical plausibility of the models. The compositions of Uranus and Neptune are found to be similar with somewhat different distributions of the high-Z material. The big difference between the two planets is that Neptune requires a non-solar envelope while Uranus is best matched with a solar composition envelope. Our analysis suggests that the heavier elements in both Uranus' and Neptune's interior might increase gradually towards the planetary centers. Indeed it is possible to fit the gravitational moments without sharp compositional transitions.

Numerical simulations, based on the core-nucleated accretion model, are presented for the formation of Jupiter at 5.2 AU in 3 primordial disks with three different assumed values of the surface density of solid particles. The grain opacities in the envelope of the protoplanet are computed using a detailed model that includes settling and coagulation of grains and that incorporates a recalculation of the grain size distribution at each point in time and space. We generally find lower opacities than the 2% of interstellar values used in previous calculations [Hubickyj, O., Bodenheimer, P., Lissauer, J. J., 2005. Icarus 179, 415--431; Lissauer, J. J., Hubickyj, O., D'Angelo, G., Bodenheimer, P., 2009. Icarus 199, 338-350]. These lower opacities result in more rapid heat loss from and more rapid contraction of the protoplanetary envelope. For a given surface density of solids, the new calculations result in a substantial speedup in formation time as compared with those previous calculations. Formation times are calculated to be 1.0, 1.9, and 4.0 Myr, and solid core masses are found to be 16.8, 8.9, and 4.7 Earth masses, for solid surface densities, sigma, of 10, 6, and 4 grams per squared centimeter, respectively. For sigma=10 and sigma=6 g/cm^2, respectively, these formation times are reduced by more than 50% and more than 80% compared with those in a previously published calculation with the old approximation to the opacity.

The snow line in a gas disk is defined as the distance from the star beyond which the water ice is stable against evaporation. Since oxygen is the most abundant element after hydrogen and helium, the presence of ice grains can have important consequences for disk evolution. However, determining the position of the snow line is not simple. I discuss some of the important processes that affect the position of the snow line.

We present a calculation of the sedimentation of grains in a giant gaseous protoplanet such as that resulting from a disk instability of the type envisioned by Boss (1998). Boss (1998) has suggested that such protoplanets would form cores through the settling of small grains. We have tested this suggestion by following the sedimentation of small silicate grains as the protoplanet contracts and evolves. We find that during the course of the initial contraction of the protoplanet, which lasts some $4\times 10^5$ years, even very small (> 1 micron) silicate grains can sediment to create a core both for convective and non-convective envelopes, although the sedimentation time is substantially longer if the envelope is convective, and grains are allowed to be carried back up into the envelope by convection. Grains composed of organic material will mostly be evaporated before they get to the core region, while water ice grains will be completely evaporated. These results suggest that if giant planets are formed via the gravitational instability mechanism, a small heavy element core can be formed due to sedimentation of grains, but it will be composed almost entirely of refractory material. Including planetesimal capture, we find core masses between 1 and 10 M$_{\oplus}$, and a total high-Z enhancement of ~40 M$_{\oplus}$. The refractories in the envelope will be mostly water vapor and organic residuals.

We have computed the size distribution of silicate grains in the outer radiative region of the envelope of a protoplanet evolving according to the scenario of Pollack et al. (1996). Our computation includes grain growth due to Brownian motion and overtake of smaller grains by larger ones. We also include the input of new grains due to the breakup of planetesimals in the atmosphere. We follow the procedure of Podolak (2003), but have speeded it up significantly. This allows us to test the sensitivity of the code to various parameters. We have also made a more careful estimate of the resulting grain opacity. We find that the grain opacity is of the order of $10^{-2}\ \mathrm{cm^2 g^{-1}}$ throughout most of the outer radiative zone as Hubickyj et al. (2005) assumed for their low opacity case, but near the outer edge of the envelope, the opacity can increase to $\sim{1} \mathrm{cm^2 g^{-1}}$. We discuss the effect of this on the evolution of the models.

We follow the contraction and evolution of a typical Jupiter-mass clump created by the disk instability mechanism, and compute the rate of planetesimal capture during this evolution. We show that such a clump has a slow contraction phase lasting ~3x10^5 years. By following the trajectories of planetesimals as they pass through the envelope of the protoplanet, we compute the cross-section for planetesimal capture at all stages of the protoplanet's evolution. We show that the protoplanet can capture a large fraction of the solid material in its feeding zone, which will lead to an enrichment of the protoplanet in heavy elements. The exact amount of this enrichment depends upon, but is not very sensitive to the size and random speed of the planetesimals.

In a protoplanetary disk, the inner edge of the region where the temperature falls below the condensation temperature of water is referred to as the 'snow line'. Outside the snow line, water ice increases the surface density of solids by a factor of 4. The mass of the fastest growing planetesimal (the 'isolation mass') scales as the surface density to the 3/2 power. It is thought that ice-enhanced surface densities are required to make the cores of the gas giants (Jupiter and Saturn) before the disk gas dissipates. Observations of the Solar System's asteroid belt suggest that the snow line occurred near 2.7 AU. In this paper we revisit the theoretical determination of the snow line. In a minimum-mass disk characterized by conventional opacities and a mass accretion rate of 10^-8 solar masses per year, the snow line lies at 1.6-1.8 AU, just past the orbit of Mars. The minimum-mass disk, with a mass of 0.02 solar, has a life time of 2 million years with the assumed accretion rate. Moving the snow line past 2.7 AU requires that we increase the disk opacity, accretion rate, and/or disk mass by factors ranging up to an order of magnitude above our assumed baseline values.