Planet Formation Models Research Articles

Color Dependence of Planetary Transit Depths due to Large Starspots and Dust Clouds: Application to PTFO 8-8695

Abstract Although thousands of exoplanets have now been discovered, there is still a significant lack of observations of young planets only a few Myr old. Thus there is little direct evidence available to differentiate between various models of planet formation. The detection of planets of this age would provide much-needed data that could help constrain the planet formation process. To explore what transit observations of such planets may look like, we model the effects of large starspots and dust clouds on the depths of exoplanet transits across multiple wavelengths. We apply this model to the candidate planet PTFO 8-8695b, whose depths vary significantly across optical and infrared wavelengths. Our model shows that, while large starspots can significantly increase the color dependence of planetary transits, a combination of starspots and a large cloud surrounding the planet is required to reproduce the observed transit depths across four wavelengths.

Molecules with ALMA at Planet-forming Scales (MAPS). I. Program Overview and Highlights

Abstract Planets form and obtain their compositions in dust- and gas-rich disks around young stars, and the outcome of this process is intimately linked to the disk chemical properties. The distributions of molecules across disks regulate the elemental compositions of planets, including C/N/O/S ratios and metallicity (O/H and C/H), as well as access to water and prebiotically relevant organics. Emission from molecules also encodes information on disk ionization levels, temperature structures, kinematics, and gas surface densities, which are all key ingredients of disk evolution and planet formation models. The Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program was designed to expand our understanding of the chemistry of planet formation by exploring disk chemical structures down to 10 au scales. The MAPS program focuses on five disks—around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480—in which dust substructures are detected and planet formation appears to be ongoing. We observed these disks in four spectral setups, which together cover ∼50 lines from over 20 different species. This paper introduces the Astrophysical Journal Supplement’s MAPS Special Issue by presenting an overview of the program motivation, disk sample, observational details, and calibration strategy. We also highlight key results, including discoveries of links between dust, gas, and chemical substructures, large reservoirs of nitriles and other organics in the inner disk regions, and elevated C/O ratios across most disks. We discuss how this collection of results is reshaping our view of the chemistry of planet formation.

Limits on the mass and initial entropy of 51 Eri b from Gaia EDR3 astrometry

ABSTRACT 51 Eri b is one of the only young planets consistent with a wide range of possible initial entropy states, including the cold-start scenario associated with some models of planet formation by core accretion. The most direct way to constrain the initial entropy of a planet is by measuring its luminosity and mass at a sufficiently young age that the initial conditions still matter. We present the tightest upper limit on 51 Eri b’s mass yet (M < 11 MJup at 2σ) using a cross-calibration of Hipparcos and Gaia EDR3 astrometry and the orbit-fitting code orvara. We also reassess its luminosity using a direct, photometric approach, finding $\log (\rm{L_{\rm bol}}/\rm{\mathrm{L}_{\odot }}) = -5.5\pm 0.2$ dex. Combining this luminosity with the 24 ± 3 Myr age of the β Pic moving group, of which 51 Eri is a member, we derive mass distributions from a grid of evolutionary models that spans a wide range of initial entropies. We find that 51 Eri b is inconsistent with the coldest-start scenarios, requiring an initial entropy of >8 kB baryon−1 at 97 per cent confidence. This result represents the first observational constraint on the initial entropy of a potentially cold-start planet, and it continues the trend of dynamical masses for directly imaged planets pointing to warm- or hot-start formation scenarios.

The Occurrence-weighted Median Planets Discovered by Transit Surveys Orbiting Solar-type Stars and Their Implications for Planet Formation and Evolution

Abstract Since planet occurrence and primordial atmospheric retention probability increase with period, the occurrence-weighted median planets discovered by transit surveys may bear little resemblance to the low-occurrence, short-period planets sculpted by atmospheric escape ordinarily used to calibrate mass–radius relations and planet formation models. An occurrence-weighted mass–radius relation for the low-mass planets discovered so far by transit surveys orbiting solar-type stars requires both occurrence-weighted median Earth-mass and Neptune-mass planets to have a few percent of their masses in hydrogen/helium (H/He) atmospheres. Unlike the Earth that finished forming long after the protosolar nebula was dissipated, these occurrence-weighted median Earth-mass planets must have formed early in their systems’ histories. The existence of significant H/He atmospheres around Earth-mass planets confirms an important prediction of the core-accretion model of planet formation. It also implies core masses M c in the range 2 M ⊕ ≲ M c ≲ 8 M ⊕ that can retain their primordial atmospheres. If atmospheric escape is driven by photoevaporation due to extreme-ultraviolet (EUV) flux, then our observation requires a reduction in the fraction of incident EUV flux converted into work usually assumed in photoevaporation models. If atmospheric escape is core driven, then the occurrence-weighted median Earth-mass planets must have large Bond albedos. In contrast to Uranus and Neptune that have at least 10% of their masses in H/He atmospheres, these occurrence-weighted median Neptune-mass planets are H/He poor. The implication is that they experienced collisions or formed in much shorter-lived and/or hotter parts of their parent protoplanetary disks than Uranus and Neptune’s formation location in the protosolar nebula.

Zodiacal exoplanets in time (ZEIT) XII: a directly imaged planetary-mass companion to a young Taurus M dwarf star

ABSTRACT We report the discovery of a resolved (0.9 arcsec) substellar companion to a member of the 1–5 Myr Taurus star-forming region. The host star (2M0437) is a single mid-M type (Teff ≈ 3100 K) dwarf with a position, space motion, and colour–magnitude that support Taurus membership, and possible affiliation with a ∼2.5-Myr-old subgroup. A comparison with stellar models suggests a 2–5 Myr age and a mass of 0.15–0.18M⊙. Although K2 detected quasi-periodic dimming from close-in circumstellar dust, the star lacks detectable excess infrared emission from a circumstellar disc and its H α emission is not commensurate with accretion. Astrometry based on 3 yr of AO imaging shows that the companion (2M0437b) is comoving, while photometry of two other sources at larger separation indicates that they are likely heavily reddened background stars. A comparison of the luminosity of 2M0437b with models suggests a mass of 3–5 MJUP, well below the deuterium burning limit, and an effective temperature of 1400–1500 K, characteristic of a late L spectral type. The H − K colour is redder than the typical L dwarf, but comparable to other directly detected young planets, e.g. those around HR 8799. The discovery of a super-Jupiter around a very young, very low-mass star challenges models of planet formation by either core accretion (which requires time) or disc instability (which requires mass). We also detected a second, comoving, widely separated (75 arcsec) object that appears to be a heavily extincted star. This is certainly a fellow member of this Taurus subgroup and statistically likely to be a bound companion.

A unique hot Jupiter spectral sequence with evidence for compositional diversity

The emergent spectra of close-in, giant exoplanets ("hot Jupiters") are expected to be distinct from those of self-luminous objects with similar effective temperatures because hot Jupiters are primarily heated from above by their host stars rather than internally from the release of energy from their formation. Theoretical models predict a continuum of dayside spectra for hot Jupiters as a function of irradiation level, with the coolest planets having absorption features in their spectra, intermediate-temperature planets having emission features due to thermal inversions, and the hottest planets having blackbody-like spectra due to molecular dissociation and continuum opacity from the H- ion. Absorption and emission features have been detected in the spectra of a number of individual hot Jupiters, and population-level trends have been observed in photometric measurements. However, there has been no unified, population-level study of the thermal emission spectra of hot Jupiters such as has been done for cooler brown dwarfs and transmission spectra of hot Jupiters. Here we show that hot Jupiter secondary eclipse spectra centered around a water absorption band at 1.4 microns follow a common trend in water feature strength with temperature. The observed trend is broadly consistent with model predictions for how the thermal structures of solar-composition planets vary with irradiation level. Nevertheless, the ensemble of planets exhibits some degree of scatter around the mean trend for solar composition planets. The spread can be accounted for if the planets have modest variations in metallicity and/or elemental abundance ratios, which is expected from planet formation models. (abridged abstract)

Carbon monoxide gas produced by a giant impact in the inner region of a young system.

Models of terrestrial planet formation predict that the final stages of planetary assembly-lasting tens of millions of years beyond the dispersal of young protoplanetary disks-are dominated by planetary collisions. It is through these giant impacts that planets like the young Earth grow to their final mass and achieve long-term stable orbital configurations1. A key prediction is that these impacts produce debris. So far, the most compelling observational evidence for post-impact debris comes from the planetary system around the nearby 23-million-year-old A-type star HD 172555. This system shows large amounts of fine dust with an unusually steep size distribution and atypical dust composition, previously attributed to either a hypervelocity impact2,3 or a massive asteroid belt4. Here we report the spectrally resolved detection of a carbon monoxide gas ring co-orbiting with dusty debris around HD 172555 between about six and nine astronomical units-a region analogous to the outer terrestrial planet region of our Solar System. Taken together, the dust and carbon monoxide detections favour a giant impact between large, volatile-rich bodies. This suggests that planetary-scale collisions, analogous to the Moon-forming impact, can release large amounts of gas as well as debris, and that this gas is observable, providing a window into the composition of young planets.

Global 3D radiation hydrodynamic simulations of proto-Jupiter’s convective envelope

ABSTRACT The core accretion model of giant planet formation has been challenged by the discovery of recycling flows between the planetary envelope and the disc that can slow or stall envelope accretion. We carry out 3D radiation hydrodynamic simulations with an updated opacity compilation to model the proto-Jupiter’s envelope. To isolate the 3D effects of convection and recycling, we simulate both isolated spherical envelopes and envelopes embedded in discs. The envelopes are heated at given rates to achieve steady states, enabling comparisons with 1D models. We vary envelope properties to obtain both radiative and convective solutions. Using a passive scalar, we observe significant mass recycling on the orbital time-scale. For a radiative envelope, recycling can only penetrate from the disc surface until ∼0.1–0.2 planetary Hill radii, while for a convective envelope, the convective motion can ‘dredge up’ the deeper part of the envelope so that the entire convective envelope is recycled efficiently. This recycling, however, has only limited effects on the envelopes’ thermal structure. The radiative envelope embedded in the disc has identical structure as the isolated envelope. The convective envelope has a slightly higher density when it is embedded in the disc. We introduce a modified 1D approach which can fully reproduce our 3D simulations. With our updated opacity and 1D model, we recompute Jupiter’s envelope accretion with a 10 M⊕ core, and the time-scale to runaway accretion is shorter than the disc lifetime as in prior studies. Finally, we discuss the implications of the efficient recycling on the observed chemical abundances of the planetary atmosphere (especially for super-Earths and mini-Neptunes).

The CARMENES search for exoplanets around M dwarfs

Context. The CARMENES exoplanet survey of M dwarfs has obtained more than 18 000 spectra of 329 nearby M dwarfs over the past five years as part of its guaranteed time observations (GTO) program. Aims. We determine planet occurrence rates with the 71 stars from the GTO program for which we have more than 50 observations. Methods. We use injection-and-retrieval experiments on the radial-velocity time series to measure detection probabilities. We include 27 planets in 21 planetary systems in our analysis. Results. We find 0.06−0.03+0.04 giant planets (100 M⊕ < Mpl sin i < 1000 M⊕) per star in periods of up to 1000 d, but due to a selection bias this number could be up to a factor of five lower in the whole 329-star sample. The upper limit for hot Jupiters (orbital period of less than 10 d) is 0.03 planets per star, while the occurrence rate of planets with intermediate masses (10 M⊕ < Mpl sin i < 100 M⊕) is 0.18−0.05+0.07 planets per star. Less massive planets with 1 M⊕ < Mpl sin i < 10 M⊕ are very abundant, with an estimated rate of 1.32−0.31+0.33 planets per star for periods of up to 100 d. When considering only late M dwarfs with masses M⋆ < 0.34 M⊙, planets more massive than 10 M⊕ become rare. Instead, low-mass planets with periods shorter than 10 d are significantly overabundant. Conclusions. For orbital periods shorter than 100 d, our results confirm the known stellar mass dependences from the Kepler survey: M dwarfs host fewer giant planets and at least two times more planets with Mpl sin i < 10 M⊕ than G-type stars. In contrast to previous results, planets around our sample of very low-mass stars have a higher occurrence rate in short-period orbits of less than 10 d. Our results demonstrate the need to take into account host star masses in planet formation models.

The importance of thermal torques on the migration of planets growing by pebble accretion

ABSTRACT A key process in planet formation is the exchange of angular momentum between a growing planet and the protoplanetary disc, which makes the planet migrate through the disc. Several works show that in general low-mass and intermediate-mass planets migrate towards the central star, unless corotation torques become dominant. Recently, a new kind of torque, called the thermal torque, was proposed as a new source that can generate outward migration of low-mass planets. While the Lindblad and corotation torques depend mostly on the properties of the protoplanetary disc and on the planet mass, the thermal torque depends also on the luminosity of the planet, arising mainly from the accretion of solids. Thus, the accretion of solids plays an important role not only in the formation of the planet but also in its migration process. In a previous work, we evaluated the thermal torque effects on planetary growth and migration mainly in the planetesimal accretion paradigm. In this new work, we study the role of the thermal torque within the pebble accretion paradigm. Computations are carried out consistently in the framework of a global model of planet formation that includes disc evolution, dust growth and evolution, and pebble formation. We also incorporate updated prescriptions of the thermal torque derived from high-resolution hydrodynamical simulations. Our simulations show that the thermal torque generates extended regions of outward migration in low-viscosity discs. This has a significant impact in the formation of the planets.

Planetesimal Dynamics in the Presence of a Giant Planet

Abstract Standard models of planet formation explain how planets form in axisymmetric, unperturbed disks in single star systems. However, it is possible that giant planets could have already formed when other planetary embryos start to grow. We investigate the dynamics of planetesimals under the perturbation of a giant planet in a gaseous disk. Our aim is to understand the effect of the planet’s perturbation on the formation of giant planet cores outside the orbit of the planet. We calculate the orbital evolution of planetesimals ranging from 1013 to 1020 g, with a Jupiter-mass planet located at 5.2 au. We find orbital alignment of planetesimals distributed at ≃9–15 au, except for the mean motion resonance (MMR) locations. The degree of alignment increases with increasing distance from the planet and decreasing planetesimal mass. Aligned orbits lead to low encounter velocity and thus faster growth. The typical velocity dispersion for identical-mass planetesimals is ∼  ( 10 ) m s − 1 except for the MMR locations. The relative velocity decreases with increasing distance from the planet and decreasing mass ratio of planetesimals. When the eccentricity vectors of planetesimals reach equilibrium under the gas drag and secular perturbation, the relative velocity becomes lower when the masses of two planetesimals are both at the larger end of the mass spectrum. Our results show that with a giant planet embedded in the disk, the growth of another planetary core outside the planet’s orbit might be accelerated in certain locations.

A Tidal Origin for a Three-body Resonance in Kepler-221

Abstract Over the course of the last two decades, traditional models of planet formation have been repeatedly challenged by the emerging census of extrasolar planets. Key among them is the orbital architecture problem: while standard models of orbital migration predict resonant orbits for short-period objects, most planets do not appear to lie in orbital resonances. Here, we show that the four-planet system Kepler-221, not previously recognized to have active orbital resonances, has a three-body commensurability relation unique within the Kepler sample. Using a suite of numerical experiments as well as a perturbative analysis, we demonstrate that this system likely began as a resonant chain and proceeded to undergo large-scale divergence away from resonance, under the action of tidal dissipation. Our results further indicate that obliquity tides, driven by a secular spin–orbit resonance and mutual inclination, are an excellent candidate for driving this orbital divergence, and that the high tidal luminosity may also explain the anomalous size of planet b, which lies within the Fulton radius gap.

What Can Meteorites Tell Us About the Formation of Jupiter?

AbstractGas giants like Jupiter are a fundamental component of planetary systems, but how they formed has been uncertain. Here we discuss how paleomagnetic records in meteorites of the solar nebula may tell us about Jupiter's final growth stage. We suggest that under certain testable assumptions, the meteorite data indicate that proto‐Jupiter grew from a mass of ∼50 Earth masses (M⨁) at >3.46 million years (Ma) after solar system formation to its final mass of 318 M⊕ over just <0.5 Ma. This rapid acceleration is consistent with a key prediction of the core accretion model for giant planet formation.

The Effect of Inefficient Accretion on Planetary Differentiation

Abstract Pairwise collisions between terrestrial embryos are the dominant means of accretion during the last stage of planet formation. Hence, their realistic treatment in N-body studies is critical to accurately model the formation of terrestrial planets and to develop interpretations of telescopic and spacecraft observations. In this work, we compare the effects of two collision prescriptions on the core−mantle differentiation of terrestrial planets: a model in which collisions are always completely accretionary (“perfect merging”), and a more realistic model based on neural networks that has been trained on hydrodynamical simulations of giant impacts. The latter model is able to predict the loss of mass due to imperfect accretion and the evolution of nonaccreted projectiles in hit-and-run collisions. We find that the results of the neural network model feature a wider range of final core mass fractions and metal−silicate equilibration pressures, temperatures, and oxygen fugacities than the assumption of perfect merging. When used to model collisions in N-body studies of terrestrial planet formation, the two models provide similar answers for planets more massive than ≈0.1 M ⊕ (Earth masses). For less massive final bodies, however, the inefficient-accretion model predicts a higher degree of compositional diversity. This phenomenon is not reflected in planet formation models of the solar system that use perfect merging to determine collisional outcomes. Our findings confirm the role of giant impacts as important drivers of planetary diversity and encourage a realistic implementation of inefficient accretion in future accretion studies.

Understanding planet formation using microgravity experiments

In 2018, images were released of a planet being formed around the star PDS 70, offering a tantalizing glimpse into how planets come into being. However, many questions remain about how dust evolves into planets, and astrophysical observations are unable to provide all the answers. It is therefore necessary to perform experiments to reveal key details and, to avoid unwanted effects from the Earth's gravitational pull, it is often necessary to perform such experiments in microgravity platforms. This Review sketches current models of planet formation and describes the experiments needed to test the models.

K2-138 g: Spitzer Spots a Sixth Planet for the Citizen Science System

Abstract K2 greatly extended Kepler’s ability to find new planets, but it was typically limited to identifying transiting planets with orbital periods below 40 days. While analyzing K2 data through the Exoplanet Explorers project, citizen scientists helped discover one super-Earth and four sub-Neptune sized planets in the relatively bright (V = 12.21, K = 10.3) K2-138 system, all which orbit near 3:2 mean-motion resonances. The K2 light curve showed two additional transit events consistent with a sixth planet. Using Spitzer photometry, we validate the sixth planet’s orbital period of 41.966 ± 0.006 days and measure a radius of , solidifying K2-138 as the K2 system with the most currently known planets. There is a sizeable gap between the outer two planets, since the fifth planet in the system, K2-138 f, orbits at 12.76 days. We explore the possibility of additional nontransiting planets in the gap between f and g. Due to the relative brightness of the K2-138 host star, and the near resonance of the inner planets, K2-138 could be a key benchmark system for both radial velocity and transit-timing variation mass measurements, and indeed radial velocity masses for the inner four planets have already been obtained. With its five sub-Neptunes and one super-Earth, the K2-138 system provides a unique test bed for comparative atmospheric studies of warm to temperate planets of similar size, dynamical studies of near-resonant planets, and models of planet formation and migration.

A New Sample of Warm Extreme Debris Disks from the ALLWISE Catalog

Abstract Extreme debris disks (EDDs) are rare systems with peculiarly large amounts of warm dust that may stem from recent giant impacts between planetary embryos during the final phases of terrestrial planet growth. Here we report on the identification and characterization of six new EDDs. These disks surround F5-G9 type main-sequence stars with ages >100 Myr, have dust temperatures higher than 300 K, and fractional luminosities between 0.01 and 0.07. Using time-domain photometric data at 3.4 and 4.6 μm from the WISE all-sky surveys, we conclude that four of these disks exhibited variable mid-infrared (IR) emission between 2010 and 2019. Analyzing the sample of all known EDDs, now expanded to 17 objects, we find that 14 of them showed changes at 3–5 μm over the past decade, suggesting that mid-IR variability is an inherent characteristic of EDDs. We also report that wide-orbit pairs are significantly more common in EDD systems than in the normal stellar population. While current models of rocky planet formation predict that the majority of giant collisions occur in the first 100 Myr, we find that the sample of EDDs is dominated by systems older than this age. This raises the possibility that the era of giant impacts may be longer than we think, or that some other mechanism(s) can also produce EDDs. We examine a scenario where the observed warm dust stems from the disruption and/or collisions of comets delivered from an outer reservoir into the inner regions, and explore what role the wide companions could play in this process.

Four Jovian planets around low-luminosity giant stars observed by the EXPRESS and PPPS

We report the discovery of planetary companions orbiting four low-luminosity giant stars with M⋆ between 1.04 and 1.39 M⊙. All four host stars have been independently observed by the EXoPlanets aRound Evolved StarS (EXPRESS) program and the Pan-Pacific Planet Search (PPPS). The companion signals were revealed by multi-epoch precision radial velocities obtained in nearly a decade. The planetary companions exhibit orbital periods between ~1.2 and 7.1 yr, minimum masses of mpsin i ~ 1.8–3.7 MJ, and eccentricities between 0.08 and 0.42. With these four new systems, we have detected planetary companions to 11 out of the 37 giant stars that are common targets in the EXPRESS and PPPS. After excluding four compact binaries from the common sample, we obtained a fraction of giant planets (mp ≳ 1– 2 MJ) orbiting within 5 AU from their parent star of f = 33.3−7.1+9.0%. This fraction is slightly higher than but consistent at the 1σ level with previous results obtained by different radial velocity surveys. Finally, this value is substantially higher than the fraction predicted by planet formation models of gas giants around stars more massive than the Sun.

Assessment of Cr Isotopic Heterogeneities of Volatile-rich Asteroids Based on Multiple Planet Formation Models

Abstract Describing the comprehensive evolutionary scenario for asteroids is key to explaining the various physical processes of the solar system. Bulk-scale carbonaceous chondrites (CCs) possibly record the primordial information associated with the formation processes of their parent bodies. In this study, we tried to estimate the relative formation region of volatile-rich asteroids by utilizing the nucleosynthetic Cr isotopic variation (54Cr/52Cr) in bulk-scale CCs. Numerical calculations were conducted to track the temporal evolution of isotopically different (solar and presolar) dust and 54Cr/52Cr values for mixed materials with disk radius. First, we found that isotopic heterogeneities in CC formation regions were preserved with a weak turbulence setting that increased the timescales of the advection and diffusion in the disk. Second, we assessed the effects of gaps formed by giant planets. Finally, the distance from the injected supernovae and Cr isotopic compositions of the presolar grains were investigated in terms of the estimated formation region of CCs. In our results, a plausible formation region of four types of CCs could be obtained with the supernova from approximately 2 pc and typical Cr isotopic compositions of presolar grains. Among the parent bodies of CCs (i.e., volatile-rich asteroids), B-type asteroids formed in the outermost region, which is inconsistent with the present population, showing that D-type asteroids are generally located beyond most of the C-complex asteroids. Both the initial and present orbits of asteroids might be explained by the scatter attributed to the inward-outward migration of Jupiter and Saturn.

TESS Discovery of a Super-Earth and Three Sub-Neptunes Hosted by the Bright, Sun-like Star HD 108236

Abstract We report the discovery and validation of four extrasolar planets hosted by the nearby, bright, Sun-like (G3V) star HD 108236 using data from the Transiting Exoplanet Survey Satellite (TESS). We present transit photometry, reconnaissance, and precise Doppler spectroscopy, as well as high-resolution imaging, to validate the planetary nature of the objects transiting HD 108236, also known as the TESS Object of Interest (TOI) 1233. The innermost planet is a possibly rocky super-Earth with a period of days and has a radius of 1.586 ± 0.098 R ⊕. The outer planets are sub-Neptunes, with potential gaseous envelopes, having radii of R ⊕, 2.72 ± 0.11 R ⊕, and R ⊕ and periods of days, days, and days, respectively. With V and K s magnitudes of 9.2 and 7.6, respectively, the bright host star makes the transiting planets favorable targets for mass measurements and, potentially, for atmospheric characterization via transmission spectroscopy. HD 108236 is the brightest Sun-like star in the visual (V) band known to host four or more transiting exoplanets. The discovered planets span a broad range of planetary radii and equilibrium temperatures and share a common history of insolation from a Sun-like star (R ⋆ = 0.888 ± 0.017 R ⊙, T eff = 5730 ± 50 K), making HD 108236 an exciting, opportune cosmic laboratory for testing models of planet formation and evolution.