The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to major uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary pathways provide a new route to understanding planetary interiors. We present the discovery of TOI-849b, the remnant core of a giant planet, with a radius smaller than Neptune but an anomalously high mass $M_p=40.8^{+2.4}_{-2.5}M_{\oplus}$ and density of $5.5\pm0.8$gcm$^{-3}$, similar to the Earth. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$% of the total mass of the planet. TOI-849b transits a late G type star (T$_{\rm mag}=11.5$) with an orbital period of 18.4 hours, leading to an equilibrium temperature of 1800K. The planet's mass is larger than the theoretical threshold mass for runaway gas accretion. As such, the planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it avoided substantial gas accretion, perhaps through gap opening or late formation. Photoevaporation rates cannot provide the mass loss required to reduce a Jupiter-like gas giant, but can remove a few $M_\oplus$ hydrogen and helium envelope on timescales of several Gyr, implying that any remaining atmosphere is likely to be enriched by water or other volatiles from the planetary interior. TOI-849b represents a unique case where material from the primordial core is left over from formation and available to study.

Exoplanet characterization is one of the main foci of current exoplanetary science. For super-Earths and sub-Neptunes, we mostly rely on mass and radius measurements, which allow to derive the body's mean density and give a rough estimate of the planet's bulk composition. However, the determination of planetary interiors is a very challenging task. In addition to the uncertainty in the observed fundamental parameters, theoretical models are limited due to the degeneracy in determining the planetary composition. We aim to study several aspects that affect internal characterization of super-Earths and sub-Neptunes: observational uncertainties, location on the M-R diagram, impact of additional constraints as bulk abundances or irradiation, and model assumptions. We use a full probabilistic Bayesian inference analysis that accounts for observational and model uncertainties. We employ a Nested Sampling scheme to efficiently produce the posterior probability distributions for all the planetary structural parameter of interest. We include a structural model based on self-consistent thermodynamics of core, mantle, high-pressure ice, liquid water, and H-He envelope. Regarding the effect of mass and radius uncertainties on the determination of the internal structure, we find three different regimes: below the Earth-like composition line and above the pure-water composition line smaller observational uncertainties lead to better determination of the core and atmosphere mass respectively, and between them structure characterization only weakly depends on the observational uncertainties. We show that small variations in the temperature or entropy profiles lead to radius variations that are comparable to the observational uncertainty, suggesting that uncertainties linked to model assumptions can become more relevant to determine the internal structure than observational uncertainties.

Materials at high pressures and temperatures are of great interest for planetary science and astrophysics, warm dense matter physics, and inertial confinement fusion research. Planetary structure models rely on our understanding of the behaviour of elements (and their mixtures) at exotic conditions that do not exist on Earth, and at the same time planets serve as natural laboratories for studying materials at extreme conditions. In this review we discuss the connection between modelling planetary interiors and high-pressure physics of hydrogen and helium. First, we summarise key experiments for determining the equation of state and phase diagram of hydrogen and helium as well as state-of-the-art theoretical approaches. We next briefly review our current knowledge of the internal structures of the giant planets in the Solar System, Jupiter and Saturn, and the importance of high pressure physics to their characterisation.

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.

Uranus and Neptune, and their diverse satellite and ring systems, represent the least explored environments of our Solar System, and yet may provide the archetype for the most common outcome of planetary formation throughout our galaxy. Ice Giants will be the last remaining class of Solar System planet to have a dedicated orbital explorer, and international efforts are under way to realise such an ambitious mission in the coming decades. In 2019, the European Space Agency released a call for scientific themes for its strategic science planning process for the 2030s and 2040s, known as Voyage 2050. We used this opportunity to review our present-day knowledge of the Uranus and Neptune systems, producing a revised and updated set of scientific questions and motivations for their exploration. This review article describes how such a mission could explore their origins, ice-rich interiors, dynamic atmospheres, unique magnetospheres, and myriad icy satellites, to address questions at the heart of modern planetary science. These two worlds are superb examples of how planets with shared origins can exhibit remarkably different evolutionary paths: Neptune as the archetype for Ice Giants, whereas Uranus may be atypical. Exploring Uranus' natural satellites and Neptune's captured moon Triton could reveal how Ocean Worlds form and remain active, redefining the extent of the habitable zone in our Solar System. For these reasons and more, we advocate that an Ice Giant System explorer should become a strategic cornerstone mission within ESA's Voyage 2050 programme, in partnership with international collaborators, and targeting launch opportunities in the early 2030s.

Dark rings with bright rims are the indirect signposts of planets embedded in protoplanetary discs. In a recent first, an azimuthally elongated AU-scale blob, possibly a planet, was resolved with ALMA in TW Hya. The blob is at the edge of a cliff-like rollover in the dust disc rather than inside a dark ring. Here we build time-dependent models of TW Hya disc. We find that the classical paradigm cannot account for the morphology of the disc and the blob. We propose that ALMA-discovered blob hides a Neptune mass planet losing gas and dust. We show that radial drift of mm-sized dust particles naturally explains why the blob is located on the edge of the dust disc. Dust particles leaving the planet perform a characteristic U-turn relative to it, producing an azimuthally elongated blob-like emission feature. This scenario also explains why a 10 Myr old disc is so bright in dust continuum. Two scenarios for the dust-losing planet are presented. In the first, a dusty pre-runaway gas envelope of about 40 Earth mass Core Accretion planet is disrupted, e.g., as a result of a catastrophic encounter. In the second, a massive dusty pre-collapse gas giant planet formed by Gravitational Instability is disrupted by the energy released in its massive core. Future modelling may discriminate between these scenarios and allow us to study planet formation in an entirely new way -- by analysing the flows of dust and gas recently belonging to planets, informing us about the structure of pre-disruption planetary envelopes.

Recent structure models of Jupiter that match Juno gravity data suggest that the planet harbours an extended region in its deep interior that is enriched in heavy elements, a so-called dilute/fuzzy core. This finding raises the question of what possible formation pathways would lead to such a structure. We model Jupiter's formation and long-term evolution, starting at late-stage formation before runaway gas accretion. The formation scenarios we consider include both primordial composition gradients, as well as gradients that are built as proto-Jupiter rapidly acquires its gaseous envelope. We then follow Jupiter's evolution as it cools down and contracts, with a particular focus on the energy and material transport in the interior. We find that none of the scenarios we consider lead to a fuzzy core that is compatible with interior structure models. In all the cases, most of Jupiter's envelope becomes convective and fully mixed after at most a few million years. This is true even when we consider a case where the gas accretion leads to a cold planet, and large amounts of heavy elements are accreted. It is therefore concluded that it is very challenging to explain Jupiter's dilute core from standard formation models. We suggest that future work should consider more complex formation pathways as well as the modelling of additional physical processes that could lead to Jupiter's current-state internal structure.

There are still many open questions regarding the nature of Uranus and Neptune, the outermost planets in the Solar System. In this review we summarize the current-knowledge about Uranus and Neptune with a focus on their composition and internal structure, formation including potential subsequent giant impacts, and thermal evolution. We present key open questions and discuss the uncertainty in the internal structures of the planets due to the possibility of non-adiabatic and inhomogeneous interiors. We also provide the reasoning for improved observational constraints on their fundamental physical parameters such as their gravitational and magnetic fields, rotation rates, and deep atmospheric composition and temperature. Only this way will we be able to improve our understating of these planetary objects, and the many similar-sized objects orbiting other stars.

The Transiting Exoplanet Survey Satellite, TESS, is currently carrying out an all-sky search for small planets transiting bright stars. In the first year of the TESS survey, steady progress was made in achieving the mission's primary science goal of establishing bulk densities for 50 planets smaller than Neptune. During that year, TESS's observations were focused on the southern ecliptic hemisphere, resulting in the discovery of three mini-Neptunes orbiting the star TOI-125, a V=11.0 K0 dwarf. We present intensive HARPS radial velocity observations, yielding precise mass measurements for TOI-125b, TOI-125c and TOI-125d. TOI-125b has an orbital period of 4.65 days, a radius of $2.726 \pm 0.075 ~\mathrm{R_{\rm E}}$, a mass of $ 9.50 \pm 0.88 ~\mathrm{M_{\rm E}}$ and is near the 2:1 mean motion resonance with TOI-125c at 9.15 days. TOI-125c has a similar radius of $2.759 \pm 0.10 ~\mathrm{R_{\rm E}}$ and a mass of $ 6.63 \pm 0.99 ~\mathrm{M_{\rm E}}$, being the puffiest of the three planets. TOI-125d, has an orbital period of 19.98 days and a radius of $2.93 \pm 0.17~\mathrm{R_{\rm E}}$ and mass $13.6 \pm 1.2 ~\mathrm{M_{\rm E}}$. For TOI-125b and TOI-125d we find unusual high eccentricities of $0.19\pm 0.04$ and $0.17^{+0.08}_{-0.06}$, respectively. Our analysis also provides upper mass limits for the two low-SNR planet candidates in the system; for TOI-125.04 ($R_P=1.36 ~\mathrm{R_{\rm E}}$, $P=$0.53 days) we find a $2\sigma$ upper mass limit of $1.6~\mathrm{M_{\rm E}}$, whereas TOI-125.05 ( $R_P=4.2^{+2.4}_{-1.4} ~\mathrm{R_{\rm E}}$, $P=$ 13.28 days) is unlikely a viable planet candidate with upper mass limit $2.7~\mathrm{M_{\rm E}}$. We discuss the internal structure of the three confirmed planets, as well as dynamical stability and system architecture for this intriguing exoplanet system.

The low luminosity of Uranus is a long-standing challenge in planetary science. Simple adiabatic models are inconsistent with the measured luminosity, which indicates that Uranus is non-adiabatic because it has thermal boundary layers and/or conductive regions. A gradual composition distribution acts as a thermal boundary to suppress convection and slow down the internal cooling. Here we investigate whether composition gradients in the deep interior of Uranus can explain its low luminosity, the required composition gradient, and whether it is stable for convective mixing on a timescale of some billion years. We varied the primordial composition distribution and the initial energy budget of the planet, and chose the models that fit the currently measured properties (radius, luminosity, and moment of inertia) of Uranus. We present several alternative non-adiabatic internal structures that fit the Uranus measurements. We found that convective mixing is limited to the interior of Uranus, and a composition gradient is stable and sufficient to explain its current luminosity. As a result, the interior of Uranus might still be very hot, in spite of its low luminosity. The stable composition gradient also indicates that the current internal structure of Uranus is similar to its primordial structure. Moreover, we suggest that the initial energy content of Uranus cannot be greater than 20% of its formation (accretion) energy. We also find that an interior with a mixture of ice and rock, rather than separated ice and rock shells, is consistent with measurements, suggesting that Uranus might not be "differentiated". Our models can explain the luminosity of Uranus, and they are also consistent with its metal-rich atmosphere and with the predictions for the location where its magnetic field is generated.

The gravity field of a giant planet is typically our best window into its interior structure and composition. Through comparison of a model planet's calculated gravitational potential with the observed potential, inferences can be made about interior quantities, including possible composition and the existence of a core. Necessarily, a host of assumptions go into such calculations, making every inference about a giant planet's structure strongly model dependent. In this work we present a more general picture by setting Saturn's gravity field, as measured during the \emph{Cassini} Grand Finale, as a likelihood function driving a Markov-Chain Monte Carlo exploration of the possible interior density profiles. The result is a posterior distribution of the interior structure that is not tied to assumed composition, thermal state, or material equations of state. Constraints on interior structure derived in this Bayesian framework are necessarily less informative, but are also less biased and more general. These empirical and probabilistic constraints on the density structure are our main data product which we archive for continued analysis. We find that the outer half of Saturn's radius is relatively well constrained, and we interpret our findings as suggesting a significant metal enrichment, in line with atmospheric abundances from remote sensing. As expected, the inner half of Saturn's radius is less well-constrained by gravity, but we generally find solutions that include a significant density enhancement, which can be interpreted as a core, although this core is often lower in density and larger in radial extent than typically found by standard models. This is consistent with a dilute core and/or composition gradients.

The goal of this work is to investigate Jupiter's growth focusing on the amount of heavy elements accreted by the planet, and its comparison with recent structure models. Our model assumes an initial core growth dominated by pebble accretion, and a second growth phase that is characterized by a moderate accretion of both planetesimals and gas. The third phase is dominated by runaway gas accretion during which the planet becomes detached from the disk. The second and third phases are computed in detail, considering two different prescriptions for the planetesimal accretion and fits from hydrodynamical studies to compute the gas accretion in the detached phase. In order for Jupiter to consist of $\sim$20-40 $M_{\oplus}$ of heavy elements as suggested by structure models, we find that Jupiter's formation location is preferably at an orbital distance of $1\lesssim a \lesssim 10$ au once the accretion of planetesimals dominates. We find that Jupiter could accrete between $\sim$1 and $\sim$15 $M_{\oplus}$ of heavy elements during runaway gas accretion, depending on the assumed initial surface density of planetesimals and the prescription used to estimate the heavy-element accretion during the final stage of the planetary formation. This would yield an envelope metallicity of $\sim$0.5 to $\sim$3 times solar. By computing the solid accretion during the detached phase, we infer a planetary mass-metallicity ($M_P$-$M_Z$) relation of $M_Z \sim M_P^{2/5}$ when a gap in the planetesimal disk is created, and of $M_Z \sim M_P^{1/6}$ without a planetesimal gap. Our hybrid pebble-planetesimal model can account for Jupiter's bulk and atmospheric enrichment. The high bulk metallicity inferred for many giant exoplanets is difficult to explain from standard formation models. This might suggest a migration history for such highly enriched giant exoplanets and/or giant impacts after the disk's dispersal.

Despite many similarities, there are significant observed differences between Uranus and Neptune: while Uranus is tilted and has a regular set of satellites, suggesting their accretion from a disk, Neptune's moons are irregular and are captured objects. In addition, Neptune seems to have an internal heat source, while Uranus is in equilibrium with solar insulation. Finally, structure models based on gravity data suggest that Uranus is more centrally condensed than Neptune. We perform a large suite of high resolution SPH simulations to investigate whether these differences can be explained by giant impacts. For Uranus, we find that an oblique impact can tilt its spin axis and eject enough material to create a disk where the regular satellites are formed. Some of the disks are massive and extended enough, and consist of enough rocky material to explain the formation of Uranus' regular satellites. For Neptune, we investigate whether a head-on collision could mix the interior, and lead to an adiabatic temperature profile, which may explain its larger flux and higher moment of inertia value. We find that massive and dense projectiles can penetrate towards the centre and deposit mass and energy in the deep interior, leading to a less centrally concentrated interior for Neptune. We conclude that the dichotomy between the ice giants can be explained by violent impacts after their formation.

The masses and radii of exoplanets are fundamental quantities needed for their characterisation. Studying the different populations of exoplanets is important for understanding the demographics of the different planetary types, which can then be linked to planetary formation and evolution. We present an updated exoplanet catalog based on reliable, robust and as much as possible accurate mass and radius measurements of transiting planets up to 120 $M_{\oplus}$. The resulting mass-radius (M-R) diagram shows two distinct populations, corresponding to rocky and volatile-rich exoplanets which overlap in both mass and radius. The rocky exoplanet population shows a relatively small density variability and ends at mass of $\sim25 M_{\oplus}$, possibly indicating the maximum core mass that can be formed. We use the composition line of pure-water to separate the two populations, and infer two new empirical M-R relations based on this data: $ M = (0.9 \pm 0.06) \ R^{(3.45 \pm 0.12)}$ for the rocky population, and $ M = (1.74 \pm 0.38) \ R^{(1.58 \pm 0.10)}$ for the volatile-rich population. While our results for the two regimes are in agreement with previous studies, the new M-R relations better match the population in the transition-region from rocky to volatile-rich exoplanets, which correspond to a mass range of 5-25 $M_{\oplus}$ and a radius range of 2-3 $R_{\oplus}$.

Context.Recent studies suggest that many giant exoplanets are highly enriched with heavy elements compared to their host star andcontain several tens of Earth masses or more of heavy elements. Such enrichment is considered to have been brought by accretionof planetesimals in late formation stages. Previous dynamical simulations, however, show that planets are unable to collect so muchheavy elements throughin situplanetesimal accretion. Aims.We investigate whether a giant planet migrating inward can capture planetesimals efficiently to significantly increase its metal-licity. Methods.We performed orbital integrations of a migrating giant planet and planetesimals in a protoplanetary gas disc to infer theplanetesimal mass that is accreted by the planet. Results.We find that the two shepherding processes of mean motion resonances trapping and aerodynamic gas drag inhibit plan-etesimal capture of a migrating planet. However, the amplified libration allows the highly-excited planetesimals in the resonances toescape from the resonance trap and be accreted by the planet. Consequently, we show that a migrating giant planet captures planetes-imals with total mass of several tens of Earth masses, if the planet forms at a few tens of AU in a relatively massive disc. We alsofind that planetesimal capture occurs efficiently in a limited range of semi-major axis, and that the total captured planetesimal massincreases with increasing migration distances. Our results have important implications for understanding the relation between giantplanet metallicity and mass, as we suggest that it reflects the formation location of the planet, or more precisely, the location whererunaway gas accretion occurred. We also suggest the observed metal-rich close-in Jupiters migrated to their present locations fromafar, where they formed.

Direct imaging (DI) surveys suggest that gas giants beyond 20 AU are rare around FGK stars. However, it is not clear what this means for the formation frequency of Gravitational Instability (GI) protoplanets due to uncertainties in gap opening and migration efficiency. Here we combine state-of-the-art calculations of homogeneous planet contraction with a population synthesis code. We find DI constraints to be satisfied if protoplanet formation by GI occurs in tens of percent of systems if protoplanets `super migrate' to small separations. In contrast, GI may occur in only a few percent of systems if protoplanets remain stranded at wide orbits because their migration is `quenched' by efficient gap opening. We then use the frequency of massive giants in radial velocity surveys inside 5 AU to break this degeneracy - observations recently showed that this population does not correlate with the host star metallicity and is therefore suspected to have formed via GI followed by inward migration. We find that only the super-migration scenario can sufficiently explain this population whilst simultaneously satisfying the DI constraints and producing the right mass spectrum of planets inside 5 AU. If massive gas-giants inside 5 AU formed via GI, then our models imply that migration must be efficient and that the formation of GI protoplanets occurs in at least a tens of percent of systems.

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.

Probing the interiors of the gas giant planets in our Solar System is not an easy task. It requires a set of accurate measurements combined with theoretical models that are used to infer the planetary composition and its depth dependence. The masses of Jupiter and Saturn are 317.83 and 95.16 Earth masses, respectively, and since a few decades, we know that they mostly consist of hydrogen and helium. It is the mass of heavy elements (all elements heavier than helium) that is not well determined, as well as their distribution within the planets. While the heavy elements are not the dominating materials in Jupiter and Saturn they are the key for our understanding of their formation and evolution histories. The planetary internal structure is inferred from theoretical models that fit the available observational constraints by using theoretical equations of states (EOSs) for hydrogen, helium, their mixtures, and heavier elements (typically rocks and/or ices). However, there is no unique solution for the planetary structure and the results depend on the used EOSs and the model assumptions imposed by the modeler. Major model assumptions that can affect the derived internal structure include the number of layers, the heat transport mechanism within the planet (and its entropy), the nature of the core (compact vs. diluted), and the location (pressure) of separation between the two envelopes. Alternative structure models assume a less distinct division between the layers and/or a non-homogenous distribution of the heavy elements. Today, with accurate measurements of the gravitational fields of Jupiter and Saturn from the Juno and Cassini missions, structure models can be further constrained. At the same time, these measurements introduce new challenges for planetary modellers.

In the standard model for giant planet formation, the planetary growth begins with accretion of solids followed by a buildup of a gaseous atmosphere as more solids are accreted, and finally, by rapid accretion of gas. The interaction of the solids with the gaseous envelope determines the subsequent planetary growth and the final internal structure. In this work we simulate the interaction of planetesimals with a growing giant planet (proto-Jupiter) and investigate how different treatments of the planetesimal-envelope interaction affect the heavy-element distribution, and the inferred core mass. We consider various planetesimal sizes and compositions as well as different ablation and radiation efficiencies, and fragmentation models. We find that in most cases the core reaches a maximum mass of 2 Earth's Mass. We show that the value of the core's mass mainly depends on the assumed size and composition of the solids, while the heavy-element distribution is also affected by the fate of the accreted planetesimals (ablation/fragmentation). Fragmentation, which is found to be important for planetesimals > 1 km, typically leads to enrichment of the inner part of the envelope while ablation results in enrichment of the outer atmosphere. Finally, we present a semi-analytical prescription for deriving the heavy-element distribution in giant protoplanets.

The standard model for giant planet formation is based on the accretion of solids by a growing planetary embryo, followed by rapid gas accretion once the planet exceeds a so-called critical mass. The dominant size of the accreted solids (cm-size particles named pebbles or km to hundred km-size bodies named planetesimals) is, however, unknown. Recently, high-precision measurements of isotopes in meteorites provided evidence for the existence of two reservoirs in the early Solar System. These reservoirs remained separated from ~1 until ~ 3 Myr after the beginning of the Solar System's formation. This separation is interpreted as resulting from Jupiter growing and becoming a barrier for material transport. In this framework, Jupiter reached ~20 Earth masses within ~1 Myr and slowly grew to ~50 Earth masses in the subsequent 2 Myr before reaching its present-day mass. The evidence that Jupiter slowed down its growth after reaching 20 Earth masses for at least 2 Myr is puzzling because a planet of this mass is expected to trigger fast runaway gas accretion. Here, we use theoretical models to describe the conditions allowing for such a slow accretion and show that Jupiter grew in three distinct phases. First, rapid pebble accretion brought the major part of Jupiter's core mass. Second, slow planetesimal accretion provided the energy required to hinder runaway gas accretion during 2 Myr. Third, runaway gas accretion proceeded. Both pebbles and planetesimals therefore have an important role in Jupiter's formation.

In this study, we investigated the differences between four commonly-used exoplanet catalogs (exoplanet.eu; exoplanetarchive.ipac.caltech.edu; openexoplanetcatalogue.com; exoplanets.org) using a Kolmogorov-Smirnov (KS) test. We found a relatively good agreement in terms of the planetary parameters (mass, radius, period) and stellar properties (mass, temperature, metallicity), although a more careful analysis of the overlap and unique parts of each catalog revealed some differences. We quantified the statistical impact of these differences and their potential cause. We concluded that although statistical studies are unlikely to be significantly affected by the choice of catalog, it would be desirable to have one consistent catalog accepted by the general exoplanet community as a base for exoplanet statistics and comparison with theoretical predictions.

Constraining the planetary composition is essential for exoplanetary characterization. In this paper, we use a statistical analysis to determine the characteristic maximum (threshold) radii for various compositions for exoplanets with masses up to 25 Earth masses (M$_\oplus$). We confirm that most planets with radii larger than 1.6 Earth radius (R$_\oplus$) are not rocky, and must consist of lighter elements, as found by previous studies. We find that planets with radii above 2.6 R$_\oplus$ cannot be pure-water worlds, and must contain significant amounts of hydrogen and helium (H-He). We find that planets with radii larger than about 3 R$_\oplus$, 3.6 R$_\oplus$, and 4.3 R$_\oplus$ are expected to consist of 2%, 5% and 10% of H-He, respectively. We investigate the sensitivity of the results to the assumed internal structure, the planetary temperature and albedo, and the accuracy of the mass and radius determination. We show that the envelope's metallicity, the percentage of H-He and the distribution of the elements play a significant role in the determination of the threshold radius. Finally, we conclude that despite the degenerate nature of the problem, it is possible to put limits on the possible range of compositions for planets with well-measured mass and radius.

The origin of Mercury's high iron-to-rock ratio is still unknown. In this work we investigate Mercury's formation via giant impacts and consider the possibilities of a single giant impact, a hit-and-run, and multiple collisions in one theoretical framework. We study the standard collision parameters (impact velocity, mass ratio, impact parameter), along with the impactor's composition and the cooling of the target. It is found that the impactor's composition affects the iron distribution within the planet and the final mass of the target by up to 15\%, although the resulting mean iron fraction is similar. We suggest that an efficient giant impact requires to be head-on with high velocities, while in the hit-and-run case the impact can occur closer to the most probable collision angle (45$^{\circ}$). It is also shown that Mercury's current iron-to-rock ratio can be a result of multiple-collisions, with their exact number depending on the collision parameters. Mass loss is found to be more significant when the collisions are tight in time.

In this chapter we summarize current knowledge of the internal structure of giant planets. We concentrate on the importance of heavy elements and their role in determining the planetary composition and internal structure, in planet formation, and during the planetary long-term evolution. We briefly discuss how internal structure models are derived, present the possible structures of the outer planets in the Solar System, and summarise giant planet formation and evolution. Finally, we introduce giant exoplanets and discuss how they can be used to better understand giant planets as a class of planetary objects.

The study of extrasolar planets and of the Solar System provides complementary pieces of the mosaic represented by the process of planetary formation. Exoplanets are essential to fully grasp the huge diversity of outcomes that planetary formation and the subsequent evolution of the planetary systems can produce. The orbital and basic physical data we currently possess for the bulk of the exoplanetary population, however, do not provide enough information to break the intrinsic degeneracy of their histories, as different evolutionary tracks can result in the same final configurations. The lessons learned from the Solar System indicate us that the solution to this problem lies in the information contained in the composition of planets. The goal of the Atmospheric Remote-Sensing Infrared Exoplanet Large-survey (ARIEL), one of the three candidates as ESA M4 space mission, is to observe a large and diversified population of transiting planets around a range of host star types to collect information on their atmospheric composition. ARIEL will focus on warm and hot planets to take advantage of their well-mixed atmospheres, which should show minimal condensation and sequestration of high-Z materials and thus reveal their bulk composition across all main cosmochemical elements. In this work we will review the most outstanding open questions concerning the way planets form and the mechanisms that contribute to create habitable environments that the compositional information gathered by ARIEL will allow to tackle

The formation history of giant planets determines their primordial structure and consequent evolution. We simulate various formation paths of Jupiter to determine its primordial entropy, and find that a common outcome is for proto-Jupiter to have non-convective regions in its interior. We use planet formation models to calculate how the entropy and post-formation luminosity depend on model properties such as the solid accretion rate and opacity, and show that the gas accretion rate and its time evolution play a key role in determining the entropy profile. The predicted luminosity of Jupiter shortly after formation varies by a factor of $2$--$3$ for different choices of model parameters. We find that entropy gradients inside Jupiter persist for $\sim 10\ {\rm Myr}$ after formation. We suggest that these gradients should be considered together with heavy-element composition gradients when modeling Jupiter's evolution and internal structure.

Recent formation and structure models of Jupiter suggest that the planet can have composition gradients and not be fully convective (adiabatic). This possibility directly affects our understanding of Jupiter's bulk composition and origin. In this Letter we present Jupiter's evolution with a primordial structure consisting of a relatively steep heavy-element gradient of 40 Earth masses. We show that for a primordial structure with composition gradients, most of the mixing occurs in the outer part of the gradient during the early evolution (several 10^7 years), leading to an adiabatic outer envelope (60% of Jupiter's mass). We find that the composition gradient in the deep interior persists, suggesting that about 40% of Jupiter's mass can be non-adiabatic with a higher temperature than the one derived from Jupiter's atmospheric properties. The region that can potentially develop layered-convection in Jupiter today is estimated to be limited to about 10% of the mass.

We present a semi-analytical population synthesis model of protoplanetary clumps formed by disk instability at radial distances of 80 - 120 AU. Various clump density profiles, initial mass functions, protoplanetary disk models, stellar masses, and gap opening criteria are considered. When we use more realistic gap opening criteria, we find that gaps open only rarely which strongly affects clump survival rates and their physical properties (mass, radius and radial distance). The inferred surviving population is then shifted towards less massive clumps at smaller radial distances. We also find that populations of surviving clumps are very sensitive to the model assumptions and used parameters. Depending on the chosen parameters, the protoplanets occupy a mass range between 0.01 and 16 $M_{\text{J}}$ and may either orbit close to the central star or as far out as 75 AU, with a sweet spot at 10-30 AU for the massive ones. However, in all the cases we consider we find that massive giant planets at very large radial distances are rare in qualitative agreement with current direct imaging surveys. We conclude that caution should be taken in deriving population synthesis models as well as when comparing the models' results with observations.

Understanding planetary interiors is directly linked to our ability of simulating exotic quantum mechanical systems such as hydrogen (H) and hydrogen-helium (H-He) mixtures at high pressures and temperatures. Equations of State (EOSs) tables based on Density Functional Theory (DFT), are commonly used by planetary scientists, although this method allows only for a qualitative description of the phase diagram, due to an incomplete treatment of electronic interactions. Here we report Quantum Monte Carlo (QMC) molecular dynamics simulations of pure H and H-He mixture. We calculate the first QMC EOS at 6000 K for an H-He mixture of a proto-solar composition, and show the crucial influence of He on the H metallization pressure. Our results can be used to calibrate other EOS calculations and are very timely given the accurate determination of Jupiter's gravitational field from the NASA Juno mission and the effort to determine its structure.

Planets in the mass range from 2 to 15 MEarth are very diverse. Some of them have low densities, while others are very dense. By measuring the masses and radii, the mean densities, structure, and composition of the planets are constrained. These parameters also give us important information about their formation and evolution, and about possible processes for atmospheric loss.We determined the masses, radii, and mean densities for the two transiting planets orbiting K2-106. The inner planet has an ultra-short period of 0.57 days. The period of the outer planet is 13.3 days.Although the two planets have similar masses, their densities are very different. For K2-106b we derive Mb=8.36-0.94+0.96 MEarh, Rb=1.52+/-0.16 REarth, and a high density of 13.1-3.6+5.4 gcm-3. For K2-106c, we find Mc=5.8-3.0+3.3 MEarth, Rc=2.50-0.26+0.27 REarth and a relatively low density of 2.0-1.1+1.6 gcm-3.Since the system contains two planets of almost the same mass, but different distances from the host star, it is an excellent laboratory to study atmospheric escape. In agreement with the theory of atmospheric-loss processes, it is likely that the outer planet has a hydrogen-dominated atmosphere. The mass and radius of the inner planet is in agreement with theoretical models predicting an iron core containing 80+20-30% of its mass. Such a high metal content is surprising, particularly given that the star has an ordinary (solar) metal abundance. We discuss various possible formation scenarios for this unusual planet.

Mini-Neptunes seem to be common planets. In this work we investigate the possible formation histories and predicted occurrence rates of mini-Neptunes assuming the planets form beyond the iceline. We consider pebble and planetesimal accretion accounting for envelope enrichment and two different opacity conditions. We find that the formation of mini-Neptunes is a relatively frequent output when envelope enrichment by volatiles is included, and that there is a "sweet spot" for mini-Neptune formation with a relatively low solid accretion rate of ~10^{-6} Earth masses per year. This rate is typical for low/intermediate-mass protoplanetary disks and/or disks with low metallicities. With pebble accretion, envelope enrichment and high opacity favor the formation of mini-Neptunes, with more efficient formation at large semi-major axes (~30 AU) and low disk viscosity. For planetesimal accretion, such planets can form also without enrichment, with the opacity being a key aspect in the growth history and favorable formation location. Finally, we show that the formation of Neptune-like planets remains a challenge for planet formation theories.

The Juno spacecraft has measured Jupiter's low-order, even gravitational moments, $J_2$--$J_8$, to an unprecedented precision, providing important constraints on the density profile and core mass of the planet. Here we report on a selection of interior models based on ab initio computer simulations of hydrogen-helium mixtures. We demonstrate that a dilute core, expanded to a significant fraction of the planet's radius, is helpful in reconciling the calculated $J_n$ with Juno's observations. Although model predictions are strongly affected by the chosen equation of state, the prediction of an enrichment of $Z$ in the deep, metallic envelope over that in the shallow, molecular envelope holds. We estimate Jupiter's core to contain an 7--25 Earth mass of heavy elements. We discuss the current difficulties in reconciling measured $J_n$ with the equations of state, and with theory for formation and evolution of the planet.

Satellite tracking involves sending electromagnetic signals to Earth. Both the orbit of the spacecraft and the electromagnetic signals themselves are affected by the curvature of spacetime. The arrival time of the pulses is compared to the ticks of local clocks to reconstruct the orbital path of the satellite to high accuracy, and to implicitly measure general relativistic effects. In particular, Schwarzschild space curvature (static) and frame-dragging (stationary) due to the planet's spin affect the satellite's orbit. The dominant relativistic effect on the path of the signal photons is Shapiro delay due to static space curvature. We compute these effects for some current and proposed space missions, using a Hamiltonian formulation in four dimensions. For highly eccentric orbits, such as in the Juno mission and in the Cassini Grand Finale, the relativistic effects have a kick-like nature, which could be advantageous for detecting them if their signatures are properly modeled as functions of time. Frame-dragging appears, in principle, measurable by Juno and Cassini, though not by Galileo 5 and 6. Practical measurement would require disentangling frame-dragging from the Newtonian 'foreground' such as the gravitational quadrupole which has an impact on both the spacecraft's orbit and the signal propagation. The foreground problem remains to be solved.

Giant planets are thought to have cores in their deep interiors, and the division into a heavy-element core and hydrogen-helium envelope is applied in both formation and structure models. We show that the primordial internal structure depends on the planetary growth rate, in particular, the ratio of heavy elements accretion to gas accretion. For a wide range of likely conditions, this ratio is in one-to-one correspondence with the resulting post-accretion profile of heavy elements within the planet. This flux ratio depends sensitively on the assumed solid surface density in the surrounding nebula. We suggest that giant planets' cores might not be distinct from the envelope and includes some hydrogen and helium, and the deep interior can have a gradual heavy-element structure. Accordingly, Jupiter's core may not be well-defined. Accurate measurements of Jupiter's gravitational field by Juno could put constraints on Jupiter's core mass. However, as we suggest here, the definition of Jupiter's core is complex, and the core's physical properties (mass, density) depend on the actual definition of the core and on its growth history.

Today, with the large number of detected exoplanets and improved measurements, we can reach the next step of planetary characterization. Classifying different populations of planets is not only important for our understanding of the demographics of various planetary types in the galaxy, but also for our understanding of planet formation. We explore the nature of two regimes in the planetary mass-radius (M-R) relation. We suggest that the transition between the two regimes of "small" and "large" planets, occurs at a mass of 124 \pm 7, M_Earth and a radius of 12.1 \pm 0.5, R_Earth. Furthermore, the M-R relation is R \propto M^{0.55\pm 0.02} and R \propto M^{0.01\pm0.02} for small and large planets, respectively. We suggest that the location of the breakpoint is linked to the onset of electron degeneracy in hydrogen, and therefore, to the planetary bulk composition. Specifically, it is the characteristic minimal mass of a planet which consists of mostly hydrogen and helium, and therefore its M-R relation is determined by the equation of state of these materials. We compare the M-R relation from observational data with the one derived by population synthesis calculations and show that there is a good qualitative agreement between the two samples.

We aim to present a generalized Bayesian inference method for constraining interiors of super Earths and sub-Neptunes. Our methodology succeeds in quantifying the degeneracy and correlation of structural parameters for high dimensional parameter spaces. Specifically, we identify what constraints can be placed on composition and thickness of core, mantle, ice, ocean, and atmospheric layers given observations of mass, radius, and bulk refractory abundance constraints (Fe, Mg, Si) from observations of the host star's photospheric composition. We employed a full probabilistic Bayesian inference analysis that formally accounts for observational and model uncertainties. Using a Markov chain Monte Carlo technique, we computed joint and marginal posterior probability distributions for all structural parameters of interest. We included state-of-the-art structural models based on self-consistent thermodynamics of core, mantle, high-pressure ice, and liquid water. Furthermore, we tested and compared two different atmospheric models that are tailored for modeling thick and thin atmospheres, respectively. First, we validate our method against Neptune. Second, we apply it to synthetic exoplanets of fixed mass and determine the effect on interior structure and composition when (1) radius, (2) atmospheric model, (3) data uncertainties, (4) semi-major axes, (5) atmospheric composition (i.e., a priori assumption of enriched envelopes versus pure H/He envelopes), and (6) prior distributions are varied. Our main conclusions are: [...]

The composition of Jupiter and the primordial distribution of the heavy elements are determined by its formation history. As a result, in order to constrain the primordial internal structure of Jupiter the growth of the core and the deposition and settling of accreted planetesimals must be followed in detail. In this paper we determine the distribution of the heavy elements in proto-Jupiter and determine the mass and composition of the core. We find that while the outer envelope of proto-Jupiter is typically convective and has an homogeneous composition, the innermost regions have compositional gradients. In addition, the existence of heavy elements in the envelope leads to much higher internal temperatures (several times 10$^4$ K) than in the case of a hydrogen-helium envelope. The derived core mass depends on the actual definition of the core: if the core is defined as the region in which the heavy-element mass fraction is above some limit (say 0.5), then it can be much more massive ($\sim$ 15 M$_{\oplus}$) and more extended (10\% of the planet's radius) than in the case where the core is just the region with 100\% heavy elements. In the former case Jupiter's core also consists of hydrogen and helium. Our results should be taken into account when constructing internal structure models of Jupiter and when interpreting the upcoming data from the {\it Juno} (NASA) mission.

We review our current understanding of the interior structure and thermal evolution of Saturn, with a focus on recent results in the Cassini era. There has been important progress in understanding physical inputs, including equations of state of planetary materials and their mixtures, physical parameters like the gravity field and rotation rate, and constraints on Saturnian free oscillations. At the same time, new methods of calculation, including work on the gravity field of rotating fluid bodies, and the role of interior composition gradients, should help to better constrain the state of Saturn's interior, now and earlier in its history. However, a better appreciation of modeling uncertainties and degeneracies, along with a greater exploration of modeling phase space, still leave great uncertainties in our understanding of Saturn's interior. Further analysis of Cassini data sets, as well as precise gravity field measurements from the Cassini Grand Finale orbits, will further revolutionize our understanding of Saturn's interior over the next few years.

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 suggest the existence of a correlation between the planetary radius and orbital period for planets with radii smaller than 4 R_Earth. Using the Kepler data, we find a correlation coefficient of 0.5120, and suggest that the correlation is not caused solely by survey incompleteness. While the correlation coefficient could change depending on the statistical analysis, the statistical significance of the correlation is robust. Further analysis shows that the correlation originates from two contributing factors. One seems to be a power-law dependence between the two quantities for intermediate periods (3-100 days), and the other is a dearth of planets with radii larger than 2 R_Earth in short periods. This correlation may provide important constraints for small-planet formation theories and for understanding the dynamical evolution of planetary systems.

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.

The alignment of Saturn's magnetic pole with its rotation axis precludes the use of magnetic field measurements to determine its rotation period. The period was previously determined from radio measurements by the Voyager spacecraft to be 10h 39m 22.4s. When the Cassini spacecraft measured a period of 10h 47m 6s, which was additionally found to change between sequential measurements, it became clear that the radio period could not be used to determine the bulk planetary rotation period. Estimates based upon Saturn's measured wind fields have increased the uncertainty even more, giving numbers smaller than the Voyager rotation period, and at present Saturn's rotation period is thought to be between 10h 32m and 10h 47m, which is unsatisfactory for such a fundamental property. Here we report a period of 10h 32m 45s +- 46s, based upon an optimization approach using Saturn's measured gravitational field and limits on the observed shape and possible internal density profiles. Moreover, even when solely using the constraints from its gravitational field, the rotation period can be inferred with a precision of several minutes. To validate our method, we applied the same procedure to Jupiter and correctly recovered its well-known rotation period.

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.

In this paper we investigate the formation of Uranus and Neptune, according to the core-nucleated accretion model, considering formation locations ranging from 12 to 30 AU from the Sun, and with various disk solid-surface densities and core accretion rates. It is shown that in order to form Uranus-like and Neptune-like planets in terms of final mass and solid-to-gas ratio, very specific conditions are required. We also show that when recently proposed high solid accretion rates are assumed, along with solid surface densities about 10 times those in the minimum-mass solar nebula, the challenge in forming Uranus and Neptune at large radial distances is no longer the formation timescale, but is rather finding agreement with the final mass and composition of these planets. In fact, these conditions are more likely to lead to gas-giant planets. Scattering of planetesimals by the forming planetary core is found to be an important effect at the larger distances. Our study emphasizes how (even slightly) different conditions in the protoplanetary disk and the birth environment of the planetary embryos can lead to the formation of very different planets in terms of final masses and compositions (solid-to-gas ratios), which naturally explains the large diversity of intermediate-mass exoplanets.

Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases abundances and the precise measurement of the helium mixing ratio have only been made available through in situ measurements by the Galileo probe. This paper describes the main scientific goals to be addressed by the future in situ exploration of Saturn placing the Galileo probe exploration of Jupiter in a broader context and before the future probe exploration of the more remote ice giants. In situ exploration of Saturn's atmosphere addresses two broad themes that are discussed throughout this paper: first, the formation history of our solar system and second, the processes at play in planetary atmospheres. In this context, we detail the reasons why measurements of Saturn's bulk elemental and isotopic composition would place important constraints on the volatile reservoirs in the protosolar nebula. We also show that the in situ measurement of CO (or any other disequilibrium species that is depleted by reaction with water) in Saturn's upper troposphere would constrain its bulk O/H ratio. We highlight the key measurements required to distinguish competing theories to shed light on giant planet formation as a common process in planetary systems with potential applications to most extrasolar systems. In situ measurements of Saturn's stratospheric and tropospheric dynamics, chemistry and cloud-forming processes will provide access to phenomena unreachable to remote sensing studies. Different mission architectures are envisaged, which would benefit from strong international collaborations.

The JUNO mission to Jupiter is planned to measure the water abundance in Jupiter's atmosphere below the cloud layer. This measurement is important because it can be used to reveal valuable information on Jupiter's origin and its composition. In this paper we discuss the importance of this measurement, the challenges in its interpretation, and address how it can be connected to interior and formation models of Jupiter.

Giant planet formation in the core accretion (CA) paradigm is predicated by the formation of a core, assembled by the coagulation of grains and later by planetesimals within a protoplanetary disc. In contrast, in the disc instability paradigm, giant planet formation is believed to be independent of core formation: massive self-gravitating gas fragments cool radiatively and collapse as a whole. We show that giant planet formation in the disc instability model may be also enhanced by core formation for reasons physically very similar to the CA paradigm. In the model explored here, efficient grain sedimentation within an initial fragment (rather than the disc) leads to the formation of a core composed of heavy elements. We find that massive atmospheres form around cores and undergo collapse as a critical core mass is exceeded, analogous to CA theory. The critical mass of the core to initiate such a collapse depends on the fragment mass and metallicity, as well as core luminosity, but ranges from less than 1 to as much as $\sim80$ Earth masses. We therefore suggest that there are two channels for the collapse of a gaseous fragment to planetary scales within the disc instability model: (i) H$_2$ dissociative collapse of the entire gaseous clump, and (ii) core-assisted gas capture, as presented here. We suggest that the first of these two is favoured in metal-poor environments and for fragments more massive than $\sim 5-10$ Jupiter masses, whereas the second is favored in metal-rich environments and fragments of lower mass. [Abridged]

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

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.

The accurate determination of Saturn's gravitational coefficients by Cassini could provide tighter constrains on Saturn's internal structure. Also, occultation measurements provide important information on the planetary shape which is often not considered in structure models. In this paper we explore how wind velocities and internal rotation affect the planetary shape and the constraints on Saturn's interior. We show that within the geodetic approach (Lindal et al., 1985, ApJ, 90, 1136) the derived physical shape is insensitive to the assumed deep rotation. Saturn's re-derived equatorial and polar radii at 100 mbar are found to be 54,445 $\pm$10 km and 60,365$\pm$10 km, respectively. To determine Saturn's interior we use {\it 1 D} three-layer hydrostatic structure models, and present two approaches to include the constraints on the shape. These approaches, however, result in only small differences in Saturn's derived composition. The uncertainty in Saturn's rotation period is more significant: with Voyager's 10h39mns period, the derived mass of heavy elements in the envelope is 0-7 M$_{\oplus}$. With a rotation period of 10h32mns, this value becomes $<4$ $M_{\oplus}$, below the minimum mass inferred from spectroscopic measurements. Saturn's core mass is found to depend strongly on the pressure at which helium phase separation occurs, and is estimated to be 5-20 M$_{\oplus}$. Lower core masses are possible if the separation occurs deeper than 4 Mbars. We suggest that the analysis of Cassini's radio occultation measurements is crucial to test shape models and could lead to constraints on Saturn's rotation profile and departures from hydrostatic equilibrium.

Knowledge of Saturn's axial moment of inertia can provide valuable information on its internal structure. We suggest that Saturn's angular momentum may be determined by the Solstice Mission (Cassini XXM) by measuring Saturn's pole precession rate and the Lense-Thirring acceleration on the spacecraft, and therefore, put constraints on Saturn's moment of inertia. It is shown that Saturn's moment of inertia can change up to ~2% due to different core properties. However, a determination of Saturn's rotation rate is required to constrain its axial moment of inertia. A change of about seven minutes in rotation period leads to a similar uncertainty in the moment of inertia value as different core properties (mass, radius). A determination of Saturn's angular momentum and rotation period by the Solstice Mission could reveal important information on Saturn's internal structure, in particular, its core properties.

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.

Since the Voyager fly-bys of Uranus and Neptune, improved gravity field data have been derived from long-term observations of the planets' satellite motions, and modified shape and solid-body rotation periods were suggested. A faster rotation period (-40 min) for Uranus and a slower rotation period (+1h20) of Neptune compared to the Voyager data were found to minimize the dynamical heights and wind speeds. We apply the improved gravity data, the modified shape and rotation data, and the physical LM-R equation of state to compute adiabatic three-layer structure models, where rocks are confined to the core, and homogeneous thermal evolution models of Uranus and Neptune. We present the full range of structure models for both the Voyager and the modified shape and rotation data. In contrast to previous studies based solely on the Voyager data or on empirical EOS, we find that Uranus and Neptune may differ to an observationally significant level in their atmospheric heavy element mass fraction Z1 and nondimensional moment of inertia, nI. For Uranus, we find Z1 < 8% and nI=0.2224(1), while for Neptune Z1 < 65% and nI=0.2555(2) when applying the modified shape and rotation data, while for the unmodified data we compute Z1 < 17% and nI=0.230(1) for Uranus and Z1 < 54% and nI=0.2410(8) for Neptune. In each of these cases, solar metallicity models (Z1=0.015) are still possible. The cooling times obtained for each planet are similar to recent calculations with the Voyager rotation periods: Neptune's luminosity can be explained by assuming an adiabatic interior while Uranus cools far too slowly. More accurate determinations of these planets' gravity fields, shapes, rotation periods, atmospheric heavy element abundances, and intrinsic luminosities are essential for improving our understanding of the internal structure and evolution of icy planets.

We model the evolution of a Jupiter-mass protoplanet formed by the disk instability mechanism at various radial distances accounting for the presence of the disk. Using three different disk models, it is found that a newly-formed Jupiter-mass protoplanet at radial distance of $\lesssim$ 5-10 AU cannot undergo a dynamical collapse and evolve further to become a gravitational bound planet. We therefore conclude that {\it giant planets, if formed by the gravitational instability mechanism, must form and remain at large radial distances during the first $\sim$ 10$^5-10^6$ years of their evolution}. The minimum radial distances in which protoplanets of 1 Saturn-mass, 3 and 5 Jupiter-mass protoplanets can evolve using a disk model with $\dot{M}=10^{-6} M_{Sun}/yr$ and $\alpha=10^{-2}$ are found to be 12, 9, and 7 AU, respectively. The effect of gas accretion on the planetary evolution of a Jupiter-mass protoplanet is also investigated. It is shown that gas accretion can shorten the pre-collapse timescale substantially. Our study suggests that the timescale of the pre-collapse stage does not only depend on the planetary mass, but is greatly affected by the presence of the disk and efficient gas accretion.

The present OSS mission continues a long and bright tradition by associating the communities of fundamental physics and planetary sciences in a single mission with ambitious goals in both domains. OSS is an M-class mission to explore the Neptune system almost half a century after flyby of the Voyager 2 spacecraft. Several discoveries were made by Voyager 2, including the Great Dark Spot (which has now disappeared) and Triton's geysers. Voyager 2 revealed the dynamics of Neptune's atmosphere and found four rings and evidence of ring arcs above Neptune. Benefiting from a greatly improved instrumentation, it will result in a striking advance in the study of the farthest planet of the Solar System. Furthermore, OSS will provide a unique opportunity to visit a selected Kuiper Belt object subsequent to the passage of the Neptunian system. It will consolidate the hypothesis of the origin of Triton as a KBO captured by Neptune, and improve our knowledge on the formation of the Solar system. The probe will embark instruments allowing precise tracking of the probe during cruise. It allows to perform the best controlled experiment for testing, in deep space, the General Relativity, on which is based all the models of Solar system formation. OSS is proposed as an international cooperation between ESA and NASA, giving the capability for ESA to launch an M-class mission towards the farthest planet of the Solar system, and to a Kuiper Belt object. The proposed mission profile would allow to deliver a 500 kg class spacecraft. The design of the probe is mainly constrained by the deep space gravity test in order to minimise the perturbation of the accelerometer measurement.

We explore the change in Jupiter's normalized axial moment of inertia (NMOI) assuming that Jupiter undergoes core erosion. It is found that Jupiter's contraction combined with an erosion of 20 M_Earth from a primordial core of 30 M_Earth can change Jupiter's NMOI over time significantly. It is shown that Jupiter's NMOI could have changed from ~0.235 to ~0.264 throughout its evolution. We find that a NMOI value of ~0.235 as suggested by dynamical models (Ward & Canup, 2006, ApJ, 640, L91) could, in principle, be consistent with Jupiter's primordial internal structure. Low NMOI values, however, persist only for the first ~ 10^6 years of Jupiter's evolution. Re-evaluation of dynamical stability models as well as more sophisticated evolution models of Jupiter with core erosion seem to be required in order to provide more robust estimates for Jupiter's primordial NMOI.

The moment of inertia of a giant planet reveals important information about the planet's internal density structure and this information is not identical to that contained in the gravitational moments. The forthcoming Juno mission to Jupiter might determine Jupiter's normalized moment of inertia NMoI=C/MR^2 by measuring Jupiter's pole precession and the Lense-Thirring acceleration of the spacecraft (C is the axial moment of inertia, and M and R are Jupiter's mass and mean radius, respectively). We investigate the possible range of NMoI values for Jupiter based on its measured gravitational field using a simple core/envelope model of the planet assuming that J_2 and J_4 are perfectly known and are equal to their measured values. The model suggests that for fixed values of J_2 and J_4 a range of NMOI values between 0.2629 and 0.2645 can be found. The Radau-Darwin relation gives a NMoI value that is larger than the model values by less than 1%. A low NMoI of ~ 0.236, inferred from a dynamical model (Ward & Canup, 2006, ApJ, 640, L91) is inconsistent with this range, but the range is model dependent. Although we conclude that the NMoI is tightly constrained by the gravity coefficients, a measurement of Jupiter's NMoI to a few tenths of percent by Juno could provide an important constraint on Jupiter's internal structure. We carry out a simplified assessment of the error involved in Juno's possible determination of Jupiter's NMoI.

The exact solution for the shape and gravitational field of a rotating two-layer Maclaurin ellipsoid of revolution is compared with predictions of the theory of figures up to third order in the small rotational parameter of the theory of figures. An explicit formula is derived for the external gravitational coefficient $J_2$ of the exact solution. A new approach to the evaluation of the theory of figures based on numerical integration of ordinary differential equations is presented. The classical Radau-Darwin formula is found not to be valid for the rotational parameter \epsilon_2 = \Omega^2/(2\pi G\rho_2) >= 0.17 since the formula then predicts a surface eccentricity that is smaller than the eccentricity of the core-envelope boundary. Interface eccentricity must be smaller than surface eccentricity. In the formula for $\epsilon_2$, $\Omega$ is the angular velocity of the two-layer body, $\rho_2$ is the density of the outer layer, and G is the gravitational constant. For an envelope density of 3000 kg m^(-3) the failure of the Radau-Darwin formula corresponds to a rotation period of about 3 hr. Application of the exact solution and the theory of figures is made to models of Earth, Mars, Uranus, and Neptune. The two-layer model with constant densities in the layers can provide realistic approximations to terrestrial planets and icy outer planet satellites. The two-layer model needs to be generalized to allow for a continuous envelope (outer layer) radial density profile in order to realistically model a gas or ice giant planet.

Transit surveys combined with Doppler data have revealed a class of gas giant planets that are massive and highly enriched in heavy elements (e.g., HD149026b, GJ436b, and HAT-P-20b). It is tempting to consider these planets as validation of core accretion plus gas capture because it is often assumed that disk instability planets should be of nebular composition. We show in this paper, to the contrary, that gas giants that form by disk instability can have a variety of heavy element compositions, ranging from sub- to super-nebular values. High levels of enrichment can be achieved through one or multiple mechanisms, including enrichment at birth, planetesimal capture, and differentiation plus tidal stripping. As a result, the metallicity of an individual gas giant cannot be used to discriminate between gas giant formation modes.

The physical shape of a giant planet can reveal important information about its centrifugal potential, and therefore, its rotation. In this paper I investigate the response of Jupiter's shape to differential rotation on cylinders of various cylindrical radii using a simple equipotential theory. I find that both solid-body rotation (with System III rotation rate) and differential rotation on cylinders up to a latitude between 20 and 30 degrees are consistent with Jupiter's measured shape. Occultation measurements of Jupiter's shape could provide an independent method to constrain the depth of its zonal winds.

'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.

The two current models for giant planet formation are core accretion and disk instability. We discuss the core masses and overall planetary enrichment in heavy elements predicted by the two formation models, and show that both models could lead to a large range of final compositions. For example, both can form giant planets with nearly stellar compositions. However, low-mass giant planets, enriched in heavy elements compared to their host stars, are more easily explained by the core accretion model. The final structure of the planets, i.e., the distribution of heavy elements, is not firmly constrained in either formation model.

Giant protoplanets formed by gravitational instability in the outer regions of circumstellar disks go through an early phase of quasi-static contraction during which radii are large and internal temperatures are low. The main source of opacity in these objects is dust grains. We investigate two problems involving the effect of opacity on the evolution of planets of 3, 5, and 7 M_J. First, we pick three different overall metallicities for the planet and simply scale the opacity accordingly. We show that higher metallicity results in slower contraction as a result of higher opacity. It is found that the pre-collapse time scale is proportional to the metallicity. In this scenario, survival of giant planets formed by gravitational instability is predicted to be more likely around low-metallicity stars, since they evolve to the point of collapse to small size on shorter time scales. But metal-rich planets, as a result of longer contraction times, have the best opportunity to capture planetesimals and form heavy-element cores. Second, we investigate the effects of opacity reduction as a result of grain growth and settling, for the same three planetary masses and for three different values of overall metallicity. When these processes are included, the pre-collapse time scale is found to be of order 1000 years for the three masses, significantly shorter than the time scale calculated without these effects. In this case the time scale is found to be relatively insensitive to planetary mass and composition. However, the effects of planetary rotation and accretion of gas and dust, which could increase the timescale, are not included in the calculation. The short time scale we find would preclude metal enrichment by planetesimal capture, as well as heavy-element core formation, over a large range of planetary masses and metallicities.

Both Uranus and Neptune are thought to have strong zonal winds with velocities of several hundred meters per second. These wind velocities, however, assume solid-body rotation periods based on Voyager 2 measurements of periodic variations in the planets' radio signals and of fits to the planets' magnetic fields; 17.24h and 16.11h for Uranus and Neptune, respectively. The realization that the radio period of Saturn does not represent the planet's deep interior rotation and the complexity of the magnetic fields of Uranus and Neptune raise the possibility that the Voyager 2 radio and magnetic periods might not represent the deep interior rotation periods of the ice giants. Moreover, if there is deep differential rotation within Uranus and Neptune no single solid-body rotation period could characterize the bulk rotation of the planets. We use wind and shape data to investigate the rotation of Uranus and Neptune. The shapes (flattening) of the ice giants are not measured, but only inferred from atmospheric wind speeds and radio occultation measurements at a single latitude. The inferred oblateness values of Uranus and Neptune do not correspond to bodies rotating with the Voyager rotation periods. Minimization of wind velocities or dynamic heights of the 1 bar isosurfaces, constrained by the single occultation radii and gravitational coefficients of the planets, leads to solid-body rotation periods of ~16.58h for Uranus and ~17.46h for Neptune. Uranus might be rotating faster and Neptune slower than Voyager rotation speeds. We derive shapes for the planets based on these rotation rates. Wind velocities with respect to these rotation periods are essentially identical on Uranus and Neptune and wind speeds are slower than previously thought. Alternatively, if we interpret wind measurements in terms of differential rotation on cylinders there are essentially no residual atmospheric winds.

The final composition of giant planets formed as a result of gravitational instability in the disk gas depends on their ability to capture solid material (planetesimals) during their 'pre-collapse' stage, when they are extended and cold, and contracting quasi-statically. The duration of the pre-collapse stage is inversely proportional roughly to the square of the planetary mass, so massive protoplanets have shorter pre-collapse timescales and therefore limited opportunity for planetesimal capture. The available accretion time for protoplanets with masses of 3, 5, 7, and 10 Jupiter masses is found to be 7.82E4, 2.62E4, 1.17E4 and 5.67E3 years, respectively. The total mass that can be captured by the protoplanets depends on the planetary mass, planetesimal size, the radial distance of the protoplanet from the parent star, and the local solid surface density. We consider three radial distances, 24, 38, and 68 AU, similar to the radial distances of the planets in the system HR 8799, and estimate the mass of heavy elements that can be accreted. We find that for the planetary masses usually adopted for the HR 8799 system, the amount of heavy elements accreted by the planets is small, leaving them with nearly stellar compositions.

The determination of Saturn's atmospheric noble gas abundances are critical to understanding the formation and evolution of Saturn, and giant planets in general. These measurements can only be performed with an entry probe. A Saturn probe will address whether enhancement in heavy noble gases, as was found in Jupiter, are a general feature of giant planets, and their ratios will be a powerful constraint on how they form. The helium abundance will show the extent to which helium has phase separated from hydrogen in the planet's deep interior. Jupiter's striking neon depletion may also be tied to its helium depletion, and must be confirmed or refuted in Saturn. Together with Jupiter's measured atmospheric helium abundance, a consistent evolutionary theory for both planets, including "helium rain" will be possible. We will then be able to calibrate the theory of the evolution of all giant planets, including exoplanets. In addition, high pressure H/He mixtures under giant planet conditions are an important area of condensed matter physics that are beyond the realm of experiment.

Anderson & Schubert (2007, Science,317,1384) proposed that Saturn's rotation period can be ascertained by minimizing the dynamic heights of the 100 mbar isosurface with respect to the geoid; they derived a rotation period of 10h 32m 35s. We investigate the same approach for Jupiter to see if the Jovian rotation period is predicted by minimizing the dynamical heights of its isobaric (1 bar pressure level) surface using zonal wind data. A rotation period of 9h 54m 29s is found. Further, we investigate the minimization method by fitting Pioneer and Voyager occultation radii for both Jupiter and Saturn. Rotation periods of 9h 55m 30s and 10h 32m 35s are found to minimize the dynamical heights for Jupiter and Saturn, respectively. Though there is no dynamical principle requiring the minimization of the dynamical heights of an isobaric surface, the successful application of the method to Jupiter lends support to its relevance for Saturn. We derive Jupiter and Saturn rotation periods using equilibrium theory in which the solid-body rotation period (no winds) that gives the observed equatorial and polar radii at the 100 mbar level is found. Rotation periods of 9h 55m 20s and 10h 31m 49s are found for Jupiter and Saturn, respectively. We show that both Jupiter's and Saturn's shapes can be derived using solid-body rotation, suggesting that zonal winds have a minor effect on the planetary shape for both planets. The agreement in the values of Saturn's rotation period predicted by the different approaches supports the conclusion that the planet's period of rotation is about 10h 32m.

We calculate heavy element enrichment in a Jupiter-mass protoplanet formed by disk instability at various radial distances from the star, considering different disk masses and surface density distributions. Although the available mass for accretion increases with radial distance (a) for disk solid surface density (sigma) functions sigma=sigma_0*a^(-alpha) with alpha < 2, the accretion timescale is significantly longer at larger radial distances. Efficient accretion is limited to the first ~ 1E5 years of planetary evolution, when the planet is extended and before gap opening and type II migration take place. The accreted mass is calculated for disk masses of 0.01, 0.05 and 0.1 M_sun with alpha = 1/2, 1, and 3/2. We show that a Jupiter-mass protoplanet can accrete 1 to 110 M_earth of heavy elements, depending on the disk properties. Our results explain the large variation in heavy element enrichment found in extra-solar giant planets. Since higher disk surface density is found to lead to larger heavy element enrichment, our model results are consistent with the correlation between heavy element enrichment and stellar metallicity. Our calculations also suggest that Jupiter could have formed at a larger radial distance than its current location while still accreting the mass of heavy elements predicted by interior models. We conclude that in the disk instability model the final composition of a giant planet is strongly determined by its formation environment. The heavy element abundance of a giant planet does not discriminate between its origin by either disk instability or core accretion.

We present 'empirical' models (pressure vs. density) of Saturn's interior constrained by the gravitational coefficients J_2, J_4, and J_6 for different assumed rotation rates of the planet. The empirical pressure-density profile is interpreted in terms of a hydrogen and helium physical equation of state to deduce the hydrogen to helium ratio in Saturn and to constrain the depth dependence of helium and heavy element abundances. The planet's internal structure (pressure vs. density) and composition are found to be insensitive to the assumed rotation rate for periods between 10h:32m:35s and 10h:41m:35s. We find that helium is depleted in the upper envelope, while in the high pressure region (P >~ 1 Mbar) either the helium abundance or the concentration of heavier elements is significantly enhanced. Taking the ratio of hydrogen to helium in Saturn to be solar, we find that the maximum mass of heavy elements in Saturn's interior ranges from ~ 6 to 20 M_Earth. The empirical models of Saturn's interior yield a moment of inertia factor varying from 0.22271 to 0.22599 for rotation periods between 10h:32m:35s and 10h:41m:35s, respectively. A long-term precession rate of about 0.754" yr^{-1} is found to be consistent with the derived moment of inertia values and assumed rotation rates over the entire range of investigated rotation rates. This suggests that the long-term precession period of Saturn is somewhat shorter than the generally assumed value of 1.77 x 10^6 years inferred from modeling and observations.

Sedimentation rates of silicate grains in gas giant protoplanets formed by disk instability are calculated for protoplanetary masses between 1 M_Saturn to 10 M_Jupiter. Giant protoplanets with masses of 5 M_Jupiter or larger are found to be too hot for grain sedimentation to form a silicate core. Smaller protoplanets are cold enough to allow grain settling and core formation. Grain sedimentation and core formation occur in the low mass protoplanets because of their slow contraction rate and low internal temperature. It is predicted that massive giant planets will not have cores, while smaller planets will have small rocky cores whose masses depend on the planetary mass, the amount of solids within the body, and the disk environment. The protoplanets are found to be too hot to allow the existence of icy grains, and therefore the cores are predicted not to contain any ices. It is suggested that the atmospheres of low mass giant planets are depleted in refractory elements compared with the atmospheres of more massive planets. These predictions provide a test of the disk instability model of gas giant planet formation. The core masses of Jupiter and Saturn were found to be ~0.25 M_Earth and ~0.5 M_Earth, respectively. The core masses of Jupiter and Saturn can be substantially larger if planetesimal accretion is included. The final core mass will depend on planetesimal size, the time at which planetesimals are formed, and the size distribution of the material added to the protoplanet. Jupiter's core mass can vary from 2 to 12 M_Earth. Saturn's core mass is found to be ~8 M_Earth.

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 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.

We investigate a uniform sample of 113 light curves (LCs) of meteors collected at the Wise Observatory in November 2002 while observing the Leonid meteor shower. We use previously defined descriptors such as the skewness F and a recently defined pointedness parameter along with a number of other measurable or derived quantities to explore the parameter space in search of meaningful LC descriptors. We make extensive use of statistical techniques to reveal links among the variables and to understand their relative importance. In particular, we show that meteors with long-duration trails rise slowly to their maximal brightness and also decay slowly from there while showing milder flaring than other meteors. Early skewed LCs show a fast rise to the peak. We show that the duration of te luminous phase of the meteor is th emost important variable differentiating among the 2002 LCs. The skewness parameter F appears only as the 2nd or 3rd in explaining the LC variance. We suggest that the pointedness parameter P could possibly be useful to discriminate among meteors from different showers, or to compare observations and model predictions.