Researcher profile

Ravit Helled

Ravit Helled contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
31works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2025arXiv

Icy or rocky? Convective or stable? New interior models of Uranus and Neptune

We present a new framework for constructing agnostic and yet physical models for planetary interiors and apply it to Uranus and Neptune. Unlike previous research that either impose rigid assumptions or rely on simplified empirical profiles, our approach bridges both paradigms. Starting from randomly generated density profiles, we applied an iterative algorithm that converges towards models that simultaneously satisfy hydrostatic equilibrium, match the observed gravitational moments, and remain thermodynamically and compositionally consistent. The inferred interior models for Uranus and Neptune span a wide range of possible interior structures, in particular encompassing both water-dominated and rock-dominated configurations (rock-to-water mass ratios between 0.04-3.92 for Uranus and 0.20-1.78 for Neptune). All models contain convective regions with ionic water and have temperature-pressure profiles that remain above the demixing curves for hydrogen-helium-water mixtures. This offers both a plausible explanation for the observed non-dipolar magnetic fields and indicates that no hydrogen-helium-water demixing occurs. We find a higher H-He mass fraction in the outermost convection zones for Uranus (0.62-0.73) compared to Neptune (0.25-0.49) and that Uranus' magnetic field is likely generated deeper in the interior compared to Neptune. We infer upper limits of 0.69-0.74 (Uranus) versus 0.78-0.92 (Neptune) for the outer edges of the dynamo regions in units of normalised radii. Overall, our findings challenge the conventional classification of Uranus and Neptune as 'ice giants' and underscore the need for improved observational data or formation constraints to break compositional degeneracy.

preprint2022arXiv

A low accretion efficiency of planetesimals formed at planetary gap edges

Observations and models of giant planets indicate that such objects are enriched in heavy elements compared to solar abundances. The prevailing view is that giant planets accreted multiple Earth masses of heavy elements after the end of core formation. Such late solid enrichment is commonly explained by the accretion of planetesimals. Planetesimals are expected to form at the edges of planetary gaps, and here we address the question of whether these planetesimals can be accreted in large enough amounts to explain the inferred high heavy element contents of giant planets. We perform a series of N-body simulations of the dynamics of planetesimals and planets during the planetary growth phase, taking into account gas drag as well as the enhanced collision cross-section caused by the extended envelopes. We consider the growth of Jupiter and Saturn via gas accretion after reaching the pebble isolation mass and we include their migration in an evolving disk. We find that the accretion efficiency of planetesimals formed at planetary gap edges is very low: less than 10% of the formed planetesimals are accreted even in the most favorable cases, which in our model corresponds to a few Earth-masses. When planetesimals are assumed to form beyond the feeding zone of the planets, extending to a few Hill radii from a planet, accretion becomes negligible. Furthermore, we find that the accretion efficiency increases when the planetary migration distance is increased and that the efficiency does not increase when the planetesimal radii are decreased. Based on these results we conclude that it is difficult to explain the large heavy element content of giant planets with planetesimal accretion during the gas accretion phase. Alternative processes most likely are required, e.g. accretion of vapor deposited by drifting pebbles.

preprint2022arXiv

Empirical Structure Models of Uranus and Neptune

Uranus and Neptune are still poorly understood. Their gravitational fields, rotation periods, atmosphere dynamics, and internal structures are not well determined. In this paper we present empirical structure models of Uranus and Neptune where the density profiles are represented by polytropes. By using these models, that are set to fit the planetary gravity field, we predict the higher order gravitational coefficients $J_6$ and $J_8$ for various assumed rotation periods, wind depths, and uncertainty of the low-order harmonics. We show that faster rotation and/or deep winds favour centrally concentrated density distributions. We demonstrate that an accurate determination of $J_6$ or $J_8$ with a relative uncertainty no larger than $10\%$ could constrain wind depths of Uranus and Neptune. We also confirm that the Voyager rotation periods are inconsistent with the measured shapes of Uranus and Neptune. We next demonstrate that more accurate determination of the gravity field can significantly reduce the possible range of internal structures. Finally, we suggest that an accurate measurement of the moment of inertia of Uranus and Neptune with a relative uncertainty of $\sim1\%$ and $\sim0.1\%$, could constrain their rotation periods and depths of the winds, respectively.

preprint2022arXiv

Enrichment of Jupiter's atmosphere by late planetesimal bombardment

Jupiter's atmosphere is enriched with heavy elements by a factor of about 3 compared to proto-solar. The origin of this enrichment and whether it represent the bulk composition of the planetary envelope remain unknown. Internal structure models of Jupiter suggest that its envelope is separated from the deep interior and that the planet is not fully mixed. This implies that Jupiter's atmosphere was enriched with heavy elements just before the end of its formation. Such enrichment can be a result of late planetesimal accretion. However, in-situ Jupiter formation models suggest the decreasing accretion rate with increasing planetary mass, which cannot explain Jupiter's atmospheric enrichment. In this study, we model Jupiter's formation and show that an migration of proto-Jupiter from $\sim$ 20 AU to its current location can lead to a late planetesimal accretion and atmospheric enrichment. Late planetesimal accretion does not occur if proto-Jupiter migrates only a few AU. We suggest that if Jupiter's outermost layer is fully-mixed and is relatively thin (up to $\sim$ 20\% of its mass), such late accretion can explain its measured atmospheric composition. It is therefore possible that Jupiter underwent significant orbital migration followed by late planetesimal accretion.

preprint2022arXiv

Giant Planets from the Inside-Out

Giant planets acquire gas, ices and rocks during the early formation stages of planetary systems and thus inform us on the formation process itself. Proceeding from inside out, examining the connections between the deep interiors and the observable atmospheres, linking detailed measurements on giant planets in the solar system to the wealth of data on brown dwarfs and giant exoplanets, we aim to provide global constraints on interiors structure and composition for models of the formation of these planets. New developments after the Juno and Cassini missions point to both Jupiter and Saturn having strong compositional gradients and stable regions from the atmosphere to the deep interior. This is also the case of Uranus and Neptune, based on available, limited data on these planets. Giant exoplanets and brown dwarfs provide us with new opportunities to link atmospheric abundances to bulk, interior abundances \rev{and to link these abundances and isotopic ratios to formation scenarios. Analysing the wealth of data becoming available} will require new models accounting for the complexity of the planetary interiors and atmospheres

preprint2022arXiv

Jupiter's Temperature Structure: A Reassessment of the Voyager Radio Occultation Measurements

The thermal structure of planetary atmospheres is an essential input for predicting and retrieving the distribution of gases and aerosols, as well as the bulk chemical abundances. In the case of Jupiter, the temperature at a reference level - generally taken at 1 bar - serves as the anchor in models used to derive the planet's interior structure and composition. Most models assume the temperature measured by the Galileo probe (Seiff et al. 1998). However, those data correspond to a single location, an unusually clear, dry region, affected by local atmospheric dynamics. On the other hand, the Voyager radio occultation observations cover a wider range of latitudes, longitudes, and times (Lindal et al. 1981). The Voyager retrievals were based on atmospheric composition and radio refractivity data that require updating and were never properly tabulated: the few existing tabulations are incomplete and ambiguous. Here, we present a systematic electronic digitization of all available temperature profiles from Voyager, followed by their reanalysis, employing currently accepted values of the abundances and radio refractivities of atmospheric species. We find the corrected temperature at the 1 bar level to be up to 4 K greater than previously published values, i.e., 170.3{\pm}3.8 K at 12°S (Voyager 1 ingress) and 167.3{\pm}3.8 K at 0°N (Voyager 1 egress). This is to be compared with the Galileo probe value of 166.1{\pm}0.8 K at the edge of an unusual feature at 6.57°N. Altogether, this suggests that Jupiter's tropospheric temperatures may vary spatially by up to 7 K between 7°N and 12°S.

preprint2022arXiv

Possible in situ formation of Uranus and Neptune via Pebble Accretion

The origin of Uranus and Neptune is still unknown. In particular, it has been challenging for planet formation models to form the planets in their current radial distances within the expected lifetime of the solar nebula. In this paper, we simulate the in-situ formation of Uranus and Neptune via pebble accretion and show that both planets can form within ~ 3 Myr at their current locations, and have final compositions that are consistent with the heavy-element to H-He ratios predicted by structure models. We find that Uranus and Neptune could have been formed at their current locations. In several cases a few earth masses of heavy elements are missing, suggesting that Uranus and/or Neptune may have accreted $\sim$1 -- 3 M Earth's mass of heavy elements after their formation via planetesimal accretion and/or giant impacts.

preprint2022arXiv

Potential long-term habitable conditions on planets with primordial H-He atmospheres

Cold super-Earths which retain their primordial, H-He dominated atmosphere could have surfaces that are warm enough to host liquid water. This would be due to the collision induced absorption (CIA) of infra-red light by hydrogen, which increases with pressure. However, the long-term potential for habitability of such planet has not been explored yet. Here we investigate the duration of this potential exotic habitability by simulating planets of different core masses, envelope masses and semi-major axes. We find that terrestrial and super-Earth planets with masses of $\sim$1 - 10$M_{\oplus}$ can maintain temperate surface conditions up to 5 - 8 Gyr at radial distances larger than $\sim $2 AU. The required envelope masses are $\sim 10^{-4} \, M_{\oplus}$ (which is 2 orders of magnitude more massive than Earth's), but can be an order of magnitude smaller (when close-in) or larger (when far out). This result suggests that the concept of planetary habitability should be revisited and made more inclusive with respect to the classical definition.

preprint2022arXiv

Revelations on Jupiter's Formation, Evolution and Interior: Challenges from Juno Results

The Juno mission has revolutionized and challenged our understanding of Jupiter. As Juno transitioned to its extended mission, we review the major findings of Jupiter's internal structure relevant to understanding Jupiter's formation and evolution. Results from Juno's investigation of Jupiter's interior structure imply that the planet has compositional gradients and is accordingly non-adiabatic, with a complex internal structure. These new results imply that current models of Jupiter's formation and evolution require a revision. In this paper, we discuss potential formation and evolution paths that can lead to an internal structure model consistent with Juno data, and the constraints they provide. We note that standard core accretion formation models, including the heavy-element enrichment during planetary growth is consistent with an interior that is inhomogeneous with composition gradients in its deep interior. However, such formation models typically predict that this region, which could be interpreted as a primordial dilute core, is confined to about 10% of Jupiter's total mass. In contrast, structure models that fit Juno data imply that this region contains 30% of the mass or more. One way to explain the origin of this extended region is by invoking a relatively long (about 2 Myrs) formation phase where the growing planet accretes gas and planetesimals delaying the runaway gas accretion. Alternatively, Jupiter's fuzzy core could be a result of a giant impact or convection post-formation. These novel scenarios require somewhat special and specific conditions. Clarity on the plausibility of such conditions could come from future high-resolution observations of planet-forming regions around other stars, from the observed and modeled architectures of extrasolar systems with giant planets, and future Juno data obtained during its extended mission.

preprint2022arXiv

The rotation of planet-hosting stars

Understanding the distribution of angular momentum during the formation of planetary systems is a key topic in astrophysics. Data from the $\textit{Kepler}$ and $\textit{Gaia}$ missions allow to investigate whether stellar rotation is correlated with the presence of planets around Sun-like stars. Here, we perform a statistical analysis of the rotation period of 493 planet-hosting stars. These are matched to a control sample, without detected planets, with similar effective temperatures, masses, radii, metallicities, and ages. We find that planet-hosting stars rotate on average $1.63 \pm 0.40$ days slower. The difference in rotation is statistically significant both in samples including and not including planets confirmed by radial velocity follow-up observations. We also analyse the dependence of rotation distribution on various stellar and planetary properties. Our results could potentially be explained by planet detection biases depending on the rotation period of their host stars in both RV and transit methods. Alternatively, they could point to a physical link between the existence of planets and stellar rotation, emphasising the need to understand the role of angular momentum in the formation and evolution planetary systems.

preprint2022arXiv

Zonal winds of Uranus and Neptune: Gravitational harmonics, dynamic self-gravity, shape, and rotation

Uranus and Neptune exhibit fast surface zonal winds that can reach up to few hundred meters per second. Previous studies on zonal gravitational harmonics and Ohmic dissipation constraints suggest that the wind speeds diminish rapidly in relatively shallow depths within the planets. Through a case-by-case comparison between the missing dynamical gravitational harmonic $J^\prime_4$ from structure models, and with that expected from fluid perturbations, we put constraints on zonal wind decay in Uranus and Neptune. To this end, we generate polytropic empirical structure models of Uranus and Neptune using $4^{\rm th}$-order Theory of Figures (ToF) that leave hydrostatic $J_4$ as an open parameter. Allotting the missing dynamical contribution to density perturbations caused by zonal winds (and their dynamic self-gravity), we find that the maximum scale height of zonal winds are $\sim 2-3\%$ of the planetary radii for both planets. Allowing the models to have $J_2$ solutions in the $\pm 5 \times 10^{-6}$ range around the observed value has similar implications. The effect of self-gravity on $J^\prime_4$ is roughly a factor of ten lower than that of zonal winds, as expected. The decay scale heights are virtually insensitive to the proposed modifications to the bulk rotation periods of Uranus and Neptune in the literature. Additionally, we find that the dynamical density perturbations due to zonal winds have a measurable impact on the shape of the planet, and could potentially be used to infer wind decay and bulk rotation period via future observations.

preprint2021arXiv

A wide-orbit giant planet in the high-mass b Centauri binary system

Planet formation occurs around a wide range of stellar masses and stellar system architectures. An improved understanding of the formation process can be achieved by studying it across the full parameter space, particularly toward the extremes. Earlier studies of planets in close-in orbits around high-mass stars have revealed an increase in giant planet frequency with increasing stellar mass until a turnover point at 1.9 solar masses, above which the frequency rapidly decreases. This could potentially imply that planet formation is impeded around more massive stars, and that giant planets around stars exceeding 3 solar masses may be rare or non-existent. However, the methods used to detect planets in small orbits are insensitive to planets in wide orbits. Here we demonstrate the existence of a planet at 560 times the Sun-Earth distance from the 6-10 solar mass binary b Centauri through direct imaging. The planet-to-star mass ratio of 0.10-0.17% is similar to the Jupiter-Sun ratio, but the separation of the detected planet is ~100 times wider than that of Jupiter. Our results show that planets can reside in much more massive stellar systems than what would be expected from extrapolation of previous results. The planet is unlikely to have formed in-situ through the conventional core accretion mechanism, but might have formed elsewhere and arrived to its present location through dynamical interactions, or might have formed via gravitational instability.

preprint2021arXiv

Internal water storage capacity of terrestrial planets and the effect of hydration on the M-R relation

Understanding the chemical interactions between water and Mg-silicates or iron is essential to constrain the interiors of water-rich planets. Hydration effects have, however, been mostly neglected by the astrophysics community so far. As such effects are unlikely to have major impacts on theoretical mass-radius relations this is justified as long as the measurement uncertainties are large. However, upcoming missions, such as the PLATO mission (scheduled launch 2026), are envisaged to reach a precision of up to $\approx 3 \%$ and $\approx 10 \%$ for radii and masses, respectively. As a result, we may soon enter an area in exoplanetary research where various physical and chemical effects such as hydration can no longer be ignored. Our goal is to construct interior models for planets that include reliable prescriptions for hydration of the cores and the mantles. These models can be used to refine previous results for which hydration has been neglected and to guide future characterization of observed exoplanets. We have developed numerical tools to solve for the structure of multi-layered planets with variable boundary conditions and compositions. Here we consider three types of planets: dry interiors, hydrated interiors and dry interiors + surface ocean where the ocean mass fraction corresponds to the mass fraction of $\rm H_2 O$ equivalent in the hydrated case. We find H/OH storage capacities in the hydrated planets equivalent to $0-6 \rm \ wt\% \ \rm H_{2}O$ corresponding to up to $\approx 800 \rm \ km$ deep ocean layers. In the mass range $0.1 \leq M/M_\oplus \leq 3$ the effect of hydration on the total radius is found to be $\leq 2.5\%$ whereas the effect of differentiation into an isolated surface ocean is $\leq 5 \ \%$. Furthermore, we find that our results are very sensitive to the bulk composition.

preprint2021arXiv

Partially Diffusive Helium-Silica Compound in the Deep Interiors of Giant Planets

Helium is the second most abundant element in the universe, and together with silica, they are major components of giant planets. Exploring the reactivity and state of helium and silica under high pressure is of fundamental importance for developing and understanding of the evolution and internal structure of giant planets. Here, using first-principles calculations and crystal structure predictions, we identify four stable phases of a helium-silica compound with seven/eight-coordinated silicon atoms at pressure range of 600-4000 GPa, corresponding to the interior condition of the outer planets in the solar system. The density of HeSiO2 agrees with current structure models of the planets. This helium-silica compound exhibits a superionic-like helium diffusive state at the high pressure and high temperature conditions along the isentropes of Saturn, a metallic fluid state in Jupiter, and a solid state in the deep interiors of Uranus and Neptune. The reaction of helium and silica may lead to the erosion of the rocky core of giant planets and form a diluted core region. These results highlight the reactivity of helium under high pressure to form new compounds, and also provides evidence to help build more sophisticated interior models of giant planets.

preprint2021arXiv

The origin of the high metallicity of close-in giant exoplanets II The nature of the sweet spot for accretion

The composition of giant planets reflects their formation history. Planetesimal accretion during the phase of planetary migration could lead to the delivery of heavy elements into giant planets. In our previous paper (Shibata et al. 2020) we showed that planetesimal accretion during planetary migration occurs in a rather narrow region of the protoplanetary disk, which we refer as "the sweet spot for accretion". The goal of this paper is to reveal the nature of the sweet spot and investigate the role of the sweet spot in determining the composition of gas giant planets. We analytically derive the required conditions for the sweet spot. Then, we compare the derived equations with the numerical simulations. We find that the conditions required for the sweet spot can be expressed by the ratio of the gas damping timescale of the planetesimal orbits and the planetary migration timescale. If the planetary migration timescale depends on the surface density of disk gas inversely, the location of the sweet spot does not change with the disk evolution. The mass of planetesimals accreted by the planet depends on the amount of planetesimals that are shepherded by mean motion resonances. Our analysis suggests that tens Earth-mass of planetesimals can be shepherded into the sweet spot without planetesimal collisions. However, as more planetesimals are trapped into mean motion resonances, collisional cascade can lead to fragmentation of planetesimals. This could affect the location of the sweet spot and the population of small objects in planetary systems. We conclude that the composition of gas giant planets depends on whether the planets crossed the sweet spot during their formation. Constraining the metallicity of cold giant planets, that are expected to be outer than the sweet spot, would reveal key information for understanding the origin of heavy elements in giant planets.

preprint2020arXiv

A remnant planetary core in the hot-Neptune desert

The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to large uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary processes provide a route to understanding planetary interiors. Planets found in and near the typically barren hot-Neptune 'desert' (a region in mass-radius space that contains few planets) have proved to be particularly valuable in this regard. These planets include HD149026b, which is thought to have an unusually massive core, and recent discoveries such as LTT9779b and NGTS-4b, on which photoevaporation has removed a substantial part of their outer atmospheres. Here we report observations of the planet TOI-849b, which has a radius smaller than Neptune's but an anomalously large mass of $39.1^{+2.7}_{-2.6}$ Earth masses and a density of $5.2^{+0.7}_{-0.8}$ grams per cubic centimetre, similar to Earth's. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$ per cent of the total planetary mass. The planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it could have avoided substantial gas accretion, perhaps through gap opening or late formation. Although photoevaporation rates cannot account for the mass loss required to reduce a Jupiter-like gas giant, they can remove a small (a few Earth masses) hydrogen and helium envelope on timescales of several billion years, implying that any remaining atmosphere on TOI-849b is likely to be enriched by water or other volatiles from the planetary interior. We conclude that TOI-849b is the remnant core of a giant planet.

preprint2020arXiv

Constraining the Depth of the Winds on Uranus and Neptune via Ohmic Dissipation

Determining the depth of atmospheric winds in the outer planets of the Solar System is a key topic in planetary science. We provide constraints on these depths in Uranus and Neptune via the total induced Ohmic dissipation, due to the interaction of the zonal flows and the planetary magnetic fields. An upper bound can be placed on the induced dissipation via energy and entropy flux throughout the interior. The induced Ohmic dissipation is directly linked to the electrical conductivity profile of the materials involved in the flow. We present a method for calculating electrical conductivity profiles of ionically conducting hydrogen-helium-water mixtures under planetary conditions, using results from ab initio simulations. We apply this prescription on several ice giant interior structure models available in the literature, where all the heavy elements are represented by water. According to the energy (entropy) flux budget, the maximum penetration depth for Uranus lies above $0.93R_{\mathrm{\scriptscriptstyle{U}}}$ ($0.90R_{\mathrm{\scriptscriptstyle{U}}}$) and for Neptune above $0.95R_{\mathrm{\scriptscriptstyle{N}}}$ ($0.92R_{\mathrm{\scriptscriptstyle{N}}}$). These results for the penetration depths are upper bounds, and are consistent with previous estimates based on the contribution of the zonal winds to the gravity field. As expected, interior structure models with higher water abundance in the outer regions have also a higher electrical conductivity and therefore reach the Ohmic limit at shallower regions. Thus, our study shows that the likelihood of deep-seated winds on Uranus and Neptune drops significantly with the presence of water in the outer layers.

preprint2020arXiv

Could Uranus and Neptune form by collisions of planetary embryos?

The origin of Uranus and Neptune remains a challenge for planet formation models. A potential explanation is that the planets formed from a population of a few planetary embryos with masses of a few Earth masses which formed beyond Saturn's orbit and migrated inwards. These embryos can collide and merge to form Uranus and Neptune. In this work we revisit this formation scenario and study the outcomes of such collisions using 3D hydrodynamical simulations. We investigate under what conditions the perfect-merging assumption is appropriate, and infer the planets' final masses, obliquities and rotation periods, as well as the presence of proto-satellite disks. We find that the total bound mass and obliquities of the planets formed in our simulations generally agree with N-body simulations therefore validating the perfect-merging assumption. The inferred obliquities, however, are typically different from those of Uranus and Neptune, and can be roughly matched only in a few cases. In addition, we find that in most cases the planets formed in this scenario rotate faster than Uranus and Neptune, close to break-up speed, and have massive disks. We therefore conclude that forming Uranus and Neptune in this scenario is challenging, and further research is required. We suggest that future planet formation models should aim to explain the various physical properties of the planets such as their masses, compositions, obliquities, rotation rates and satellite systems.

preprint2020arXiv

Detailed Calculations of the Efficiency of Planetesimal Accretion in the Core-Accretion Model

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.

preprint2020arXiv

Ice Giant Systems: The Scientific Potential of Orbital Missions to Uranus and Neptune

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.

preprint2020arXiv

TESS Reveals a Short-period Sub-Neptune Sibling (HD 86226c) to a Known Long-period Giant Planet

The Transiting Exoplanet Survey Satellite mission was designed to find transiting planets around bright, nearby stars. Here we present the detection and mass measurement of a small, short-period ($\approx\,4$\,days) transiting planet around the bright ($V=7.9$), solar-type star HD 86226 (TOI-652, TIC 22221375), previously known to host a long-period ($\sim$1600 days) giant planet. HD 86226c (TOI-652.01) has a radius of $2.16\pm0.08$ $R_{\oplus}$ and a mass of 7.25$^{+1.19}_{-1.12}$ $M_{\oplus}$ based on archival and new radial velocity data. We also update the parameters of the longer-period, not-known-to-transit planet, and find it to be less eccentric and less massive than previously reported. The density of the transiting planet is $3.97$ g cm$^{-3}$, which is low enough to suggest that the planet has at least a small volatile envelope, but the mass fractions of rock, iron, and water are not well-constrained. Given the host star brightness, planet period, and location of the planet near both the ``radius gap'' and the ``hot Neptune desert'', HD 86226c is an interesting candidate for transmission spectroscopy to further refine its composition.

preprint2020arXiv

The Challenge of Forming a Fuzzy Core in Jupiter

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.

preprint2020arXiv

The Formation of Jupiter's Diluted Core by a Giant Impact

The Juno mission has provided an accurate determination of Jupiter's gravitational field, which has been used to obtain information about the planet's composition and internal structure. Several models of Jupiter's structure that fit the probe's data suggest that the planet has a diluted core, with a total heavy-element mass ranging from ten to a few tens of Earth masses (~5-15 % of the Jovian mass), and that heavy elements (elements other than H and He) are distributed within a region extending to nearly half of Jupiter's radius. Planet-formation models indicate that most heavy elements are accreted during the early stages of a planet's formation to create a relatively compact core and that almost no solids are accreted during subsequent runaway gas accretion. Jupiter's diluted core, combined with its possible high heavy-element enrichment, thus challenges standard planet-formation theory. A possible explanation is erosion of the initially compact heavy-element core, but the efficiency of such erosion is uncertain and depends on both the immiscibility of heavy materials in metallic hydrogen and on convective mixing as the planet evolves. Another mechanism that can explain this structure is planetesimal enrichment and vaporization during the formation process, although relevant models typically cannot produce an extended diluted core. Here we show that a sufficiently energetic head-on collision (giant impact) between a large planetary embryo and the proto-Jupiter could have shattered its primordial compact core and mixed the heavy elements with the inner envelope. Models of such a scenario lead to an internal structure that is consistent with a diluted core, persisting over billions of years. We suggest that collisions were common in the young Solar system and that a similar event may have also occurred for Saturn, contributing to the structural differences between Jupiter and Saturn.

preprint2020arXiv

TOI-824 b: A New Planet on the Lower Edge of the Hot Neptune Desert

We report the detection of a transiting hot Neptune exoplanet orbiting TOI-824 (SCR J1448-5735), a nearby (d = 64 pc) K4V star, using data from the \textit{Transiting Exoplanet Survey Satellite} (TESS). The newly discovered planet has a radius, $R_{\rm{p}}$ = 2.93 $\pm$ 0.20 R$_{\oplus}$, and an orbital period of 1.393 days. Radial velocity measurements using the Planet Finder Spectrograph (PFS) and the High Accuracy Radial velocity Planet Searcher (HARPS) spectrograph confirm the existence of the planet and we estimate its mass to be $M_{\rm{p}}$ = 18.47 $\pm$ 1.84 M$_{\oplus}$. The planet's mean density is $ρ_{\rm{p}}$ = 4.03$^{+0.98}_{-0.78}$ g cm$^{-3}$ making it more than twice as dense as Neptune. TOI-824 b's high equilibrium temperature makes the planet likely to have a cloud free atmosphere, and thus an excellent candidate for follow up atmospheric studies. The detectability of TOI-824 b's atmosphere from both ground and space is promising and could lead to the detailed characterization of the most irradiated, small planet at the edge of the hot Neptune desert that has retained its atmosphere to date.

preprint2020arXiv

TW Hya: an old protoplanetary disc revived by its planet

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.

preprint2020arXiv

Understanding dense hydrogen at planetary conditions

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.

preprint2020arXiv

Uranus and Neptune: Origin, Evolution and 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.

preprint2019arXiv

Bifurcation in the history of Uranus and Neptune: the role of giant impacts

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.

preprint2019arXiv

Jupiter's heavy-element enrichment expected from formation models

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.

preprint2019arXiv

Saturn's Probable Interior: An Exploration of Saturn's Potential Interior Density Structures

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.

preprint2019arXiv

The origin of the high metallicity of close-in giant exoplanets: Combined effect of the resonant and aerodynamic shepherding

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.