We present WASP-43b climate simulations with deep wind jets (down to 700~bar) that are linked to retrograde (westward) flow at the equatorial day side for $p<0.1$~bar. Retrograde flow inhibits efficient eastward heat transport and naturally explains the small hotspot shift and large day-night-side gradient of WASP-43b ($P_{\text{orb}}=P_{\text{rot}}=0.8135$~days) observed with Spitzer. We find that deep wind jets are mainly associated with very fast rotations ($P_{\text{rot}}=P_{\text{orb}}\leq 1.5$~days) which correspond to the Rhines length smaller than $2$ planetary radii. We also diagnose wave activity that likely gives rise to deviations from superrotation. Further, we show that we can achieve full steady state in our climate simulations by imposing a deep forcing regime for $p>10$~bar: convergence time scale $\tau_{\text{conv}}=10^6-10^8$~s to a common adiabat, as well as linear drag at depth ($p\geq 200$~bar), which mimics to first order magnetic drag. Lower boundary stability and the deep forcing assumptions were also tested with climate simulations for HD~209458b ($P_{\text{orb}}=P_{\text{rot}}=3.5$~days). HD~209458b simulations always show shallow wind jets (never deeper than 100~bar) and unperturbed superrotation. If we impose a fast rotation ($P_{\text{orb}}=P_{\text{rot}}=0.8135$~days), also the HD~209458b-like simulation shows equatorial retrograde flow at the day side. We conclude that the placement of the lower boundary at $p=200$~bar is justified for slow rotators like HD~209458b, but we suggest that it has to be placed deeper for fast-rotating, dense hot Jupiters ($P_{\text{orb}}\leq 1.5$~days) like WASP-43b. Our study highlights that the deep atmosphere may have a strong influence on the observable atmospheric flow in some hot Jupiters.

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.

The large brightness variation in the observed lightcurve of Oumuamua is probably related to its shape, i.e., to the ratio between its longest axis and its shortest axis (aspect ratio). Several approaches found the aspect ratio of Oumuamua to be unusually elongated. Moreover, the spin axis orientation has to be almost perpendicular to the observer in order to obtain such an extreme lightcurve, a configuration which is unlikely. However, interstellar Oumuamua may have different surface properties than we know in our solar system. Therefore, in this work we widen the parameter space for surface properties beyond the asteroid-like models and study its effect on the lightcurve of Oumuamua. We calculate reflection from a rotating ellipsoidal object for four models: Lambertian reflection, specular reflection, single scattering diffusive and backscatter. We then calculate the probability to obtain a lightcurve ratio larger than the observed, as a function of the object aspect ratio, assuming an isotopic spin orientation distribution. We find the elongation of Oumuamua to be less extreme for the Lambertian and specular reflection models. Consequently, the probability to observe the lightcurve ratio of Oumuamua given its unknown spin axis orientation is larger for those models. We conclude that different surface reflection properties may suggest alternatives to the extreme shape of Oumuamua, relieving the need for complicated formation scenario, extreme albedo variation, or unnatural origin. Although the models suggested here are for ideal ellipsoidal shape and ideal reflection method, the results emphasize the importance of surface properties for the derived aspect ratio.

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 presence of rocky exoplanets with a large refractory carbon inventory is predicted by chemical evolution models of protoplanetary disks of stars with photospheric C/O >0.65, and by models studying the radial transport of refractory carbon. High-pressure high-temperature laboratory experiments show that most of the carbon in these exoplanets differentiates into a graphite outer shell. Our aim is to evaluate the effects of a graphite outer shell on the thermal evolution of rocky exoplanets containing a metallic core and a silicate mantle. We implement a parameterized model of mantle convection to determine the thermal evolution of rocky exoplanets with graphite layer thicknesses up to 1000 km. We find that, due to the high thermal conductivity of graphite, conduction is the dominant heat transport mechanism in a graphite layer for the long-term evolution (>200 Myr). The conductive graphite shell essentially behaves like a stagnant lid with a fixed thickness. Models of Kepler-37b (Mercury-size) and a Mars-sized exoplanet show that a planet with a graphite lid cools faster than a planet with a silicate lid, and a planet without a stagnant lid cools the fastest. A graphite lid needs to be approximately ten times thicker than a corresponding silicate lid in order to produce similar thermal evolution.

Context: Planet formation by pebble accretion is an alternative to planetesimal-driven core accretion. In this scenario, planets grow by accreting cm-to-m-sized pebbles instead of km-sized planetesimals. One of the main differences with planetesimal-driven core accretion is the increased thermal ablation experienced by pebbles. This provides early enrichment to the planet's envelope, which changes the process of core growth. Aims: We aim to predict core masses and envelope compositions of planets that form by pebble accretion and compare mass deposition of pebbles to planetesimals. Methods: We model the early growth of a proto-planet by calculating the structure of its envelope, taking into account the fate of impacting pebbles or planetesimals. The region where high-Z material can exist in vapor form is determined by the vapor pressure. We include enrichment effects by locally modifying the mean molecular weight. Results: In the pebble case, three phases of core growth can be identified. In the first phase, pebbles impact the core without significant ablation. During the second phase, ablation becomes increasingly severe. A layer of high-Z vapor starts to form around the core that absorbs a small fraction of the ablated mass. The rest either rains out to the core or mixes outwards instead, slowing core growth. In the third phase, the vapor inner region expands outwards, absorbing an increasing fraction of the ablated material as vapor. Rainout ends before the core mass reaches 0.6 M_Earth, terminating direct core growth. In the case of icy H2O pebbles, this happens before 0.1 M_Earth. Conclusions: Our results indicate that pebble accretion can directly form rocky cores up to only 0.6 M_Earth, and is unable to form similarly sized icy cores. Subsequent core growth can proceed indirectly when the planet cools, provided it is able to retain its high-Z material.

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.

Sub-Neptune planets are a very common type of planets. They are inferred to harbour a primordial (H/He) envelope, on top of a (rocky) core, which dominates the mass. Here, we investigate the long-term consequences of the core properties on the planet mass-radius relation. We consider the role of various core energy sources resulting from core formation, its differentiation, its solidification (latent heat), core contraction and radioactive decay. We divide the evolution of the rocky core into three phases: the formation phase, which sets the initial conditions, the magma ocean phase, characterized by rapid heat transport, and the solid state phase, where cooling is inefficient. We find that for typical sub-Neptune planets of ~2-10 Earth masses and envelope mass fractions of 0.5-10% the magma ocean phase lasts several Gyrs, much longer than for terrestrial planets. The magma ocean phase effectively erases any signs of the initial core thermodynamic state. After solidification, the reduced heat flux from the rocky core causes a significant drop in the rocky core surface temperature, but its effect on the planet radius is limited. In the long run, radioactive heating is the most significant core energy source in our model. Overall, the long term radius uncertainty by core thermal effects is up to 15%.

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

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.

Sub-Neptune planets are very common in our galaxy and show a large diversity in their mass-radius relation. In sub-Neptunes most of the planet mass is in the rocky part (hereafter core) which is surrounded by a modest hydrogen-helium envelope. As a result, the total initial heat content of such a planet is dominated by that of the core. Nonetheless, most studies contend that the core cooling will only have a minor effect on the radius evolution of the gaseous envelope, because the core's cooling is in sync with the envelope, i.e., most of the initial heat is released early on timescales of about 10-100 Myr. In this Letter we examine the importance of the core cooling rate for the thermal evolution of the envelope. Thus, we relax the early core cooling assumption and present a model where the core is characterized by two parameters: the initial temperature and the cooling time. We find that core cooling can significantly enhance the radius of the planet when it operates on a timescale similar to the observed age, i.e. several Gyr. Consequently, the interpretation of sub-Neptunes' mass-radius observations depends on the assumed core thermal properties and the uncertainty therein. The degeneracy of composition and core thermal properties can be reduced by obtaining better estimates of the planet ages (in addition to their radii and masses) as envisioned by future observations.

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.

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

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

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.