Formation Of super-Earths Research Articles

Buoyancy torques prevent low-mass planets from stalling in low-turbulence radiative discs

ABSTRACT Low-mass planets migrating inwards in laminar protoplanetary discs (PPDs) experience a dynamical corotation torque (DCT), which is expected to slow down migration to a stall. However, baroclinic effects can reduce or even reverse this effect, leading to rapid inward migration. In the radiatively inefficient inner disc, one such mechanism is the buoyancy response of the disc to an embedded planet. Recent work has suggested that radiative cooling can quench this response, but for parameters that are not necessarily representative of the inner regions of PPDs. We perform global 3D inviscid radiation hydrodynamics simulations of planet–disc interaction to investigate the effect of radiative cooling on the buoyancy-driven torque in a more realistic disc model. We find that the buoyancy response exerts a negative DCT – albeit partially damped due to radiative cooling – resulting in sustained, rapid inward migration. Models that adopt a local cooling prescription significantly overestimate the impact of the buoyancy response, highlighting the importance of a realistic treatment of radiation transport that includes radiative diffusion. Our results suggest that low-mass planets should migrate inwards faster than has been previously expected in radiative discs, with implications for the formation and orbital distribution of super-Earths and sub-Neptunes at intermediate distances from their host stars, unless additional physical processes that can slow down migration are considered.

Which stars can form planets: Planetesimal formation at low metallicities

The diversity of exoplanets has been linked to the disc environment in which they form, where the host star metallicity and the formation pathways play a crucial role. In the context of the core accretion paradigm, the initial stages of planet formation require the growth of dust material from micrometre-sized to planetesimal-sized bodies before core accretion can kick in. Although numerous studies have been conducted on planetesimal formation, it is still poorly understood how this process takes place in low-metallicity stellar environments. In this work, we explore how planetesimals are formed in stellar environments with primarily low metallicities. We performed global 1D viscous disc evolution simulations, including the growth of dust particles and the evaporation and condensation of chemical species at ice lines. We followed the formation of planetesimals during disc evolution and tested different metallicities, disc sizes, and turbulent viscosity strengths. We find that at solar and sub-solar metallicities, there is a significant increase in the midplane dust-to-gas mass ratios at the ice lines, but this leads to planetesimal formation only at the water--ice line. In our simulations Fe/H $= -0.6$ is the lower limit of metallicity for planetesimal formation where a few Earth masses of planetesimals can form. Our results further show that for such extreme disc environments, large discs are more conducive than small discs for forming large amounts of planetesimals at a fixed metallicity because the pebble flux can be maintained for a longer time, resulting in a longer and more efficient planetesimal formation phase. At lower metallicities, planetesimal formation is less supported in quiescent discs compared to turbulent discs, which produce larger amounts of planetesimals, because the pebble flux can be maintained for a longer time. The amount of planetesimals formed at sub-solar metallicities in our simulations places a limit on core sizes that could potentially result only in the formation of super-Earths.

Formation of rocky super-earths from a narrow ring of planetesimals

The formation of super-Earths, the most abundant planets in the Galaxy, remains elusive. These planets have masses that typically exceed that of the Earth by a factor of a few, appear to be predominantly rocky, although often surrounded by H/He atmospheres, and frequently occur in multiples. Moreover, planets that encircle the same star tend to have similar masses and radii, whereas those belonging to different systems exhibit remarkable overall diversity. Here we advance a theoretical picture for rocky planet formation that satisfies the aforementioned constraints: building upon recent work, which has demonstrated that planetesimals can form rapidly at discrete locations in the disk, we propose that super-Earths originate inside rings of silicate-rich planetesimals at approximately ~1 au. Within the context of this picture, we show that planets grow primarily through pairwise collisions among rocky planetesimals until they achieve terminal masses that are regulated by isolation and orbital migration. We quantify our model with numerical simulations and demonstrate that our synthetic planetary systems bear a close resemblance to compact, multi-resonant progenitors of the observed population of short-period extrasolar planets. A planetary origin model that forms exoplanets from a narrow ring of silicate material at a stellocentric distance of 1 au is able to explain the physical properties of super-Earths and reproduce the ‘peas in a pod’ pattern of uniformity within planetary architecture.

Formation of super-Earths in icy dead zones around low-mass stars

ABSTRACT While giant planet occurrence rates increase with stellar mass, occurrence rates of close-in super-Earths decrease. This is in contradiction to the expectation that the total mass of the planets in a system scale with the protoplanetary disc mass and hence the stellar mass. Since the snow line plays an important role in the planet formation process, we examine differences in the temperature structure of protoplanetary gas discs around stars of different mass. Protoplanetary discs likely contain a dead zone at the mid-plane that is sufficiently cold and dense for the magneto-rotational instability to be suppressed. As material builds up, the outer parts of the dead zone may be heated by self-gravity. The temperature in the disc can be below the snow line temperature far from the star and in the inner parts of a dead zone. The inner icy region has a larger radial extent around smaller mass stars. The increased mass of solid icy material may allow for the in situ formation of larger and more numerous planets close to a low-mass star. Super-Earths that form in the inner icy region may have a composition that includes a significant fraction of volatiles.

The New Generation Planetary Population Synthesis (NGPPS)

Context. Recent observational findings have suggested a positive correlation between the occurrence rates of inner super-Earths and outer giant planets. These results raise the question of whether this trend can be reproduced and explained by planet formation theory. Aims. Here, we investigate the properties of inner super-Earths and outer giant planets that form according to a core accretion scenario. We study the mutual relations between these planet species in synthetic planetary systems and compare them to the observed exoplanet population. Methods. We invoked the Generation 3 Bern model of planet formation and evolution to simulate 1000 multi-planet systems. We then confronted these synthetic systems with the observed sample, taking into account the detection bias that distorts the observed demographics. Results. The formation of warm super-Earths and cold Jupiters in the same system is enhanced compared to the individual appearances, although it is weaker than what has been proposed through observations. We attribute the discrepancy to warm and dynamically active giant planets that frequently disrupt the inner systems, particularly in high-metallicity environments. In general, a joint occurrence of the two planet types requires intermediate solid reservoirs in the originating protoplanetary disk. Furthermore, we find differences in the volatile content of planets in different system architectures and predict that high-density super-Earths are more likely to host an outer giant. This correlation can be tested observationally.

The ‘breaking the chains’ migration model for super-Earth formation: the effect of collisional fragmentation

ABSTRACT Planets between 1 and 4 R⊕ (Earth radius) with orbital periods <100 d are strikingly common. The migration model proposes that super-Earths migrate inwards and pile up at the disc inner edge in chains of mean motion resonances. After gas disc dispersal, simulations show that super-Earth’s gravitational interactions can naturally break their resonant configuration leading to a late phase of giant impacts. The instability phase is key to matching the orbital spacing of observed systems. Yet, most previous simulations have modelled collisions as perfect accretion events, ignoring fragmentation. In this work, we investigate the impact of imperfect accretion on the ‘breaking the chains’ scenario. We performed N-body simulations starting from distributions of planetary embryos and modelling the effects of pebble accretion and migration in the gas disc. Our simulations also follow the long-term dynamical evolution of super-Earths after the gas disc dissipation. We compared the results of simulations where collisions are treated as perfect merging events with those where imperfect accretion and fragmentation are allowed. We concluded that the perfect accretion is a suitable approximation in this regime, from a dynamical point of view. Although fragmentation events are common, only ∼10 per cent of the system mass is fragmented during a typical ‘late instability phase’, with fragments being mostly reacreted by surviving planets. This limited total mass in fragments proved to be insufficient to alter qualitatively the final system dynamical configuration – e.g. promote strong dynamical friction or residual migration – compared to simulations where fragmentation is neglected.

Linking Outer Disk Pebble Dynamics and Gaps to Inner Disk Water Enrichment

Abstract Millimeter continuum imaging of protoplanetary disks reveals the distribution of solid particles and the presence of substructures (gaps and rings) beyond 5–10 au, while infrared (IR) spectra provide access to abundances of gaseous species at smaller disk radii. Building on recent observational findings of an anti-correlation between the inner disk water luminosity and outer dust disk radius, we aim here at investigating the dynamics of icy solids that drift from the outer disk and sublimate their ice inside the snow line, enriching the water vapor that is observed in the IR. We use a volatile-inclusive disk evolution model to explore a range of conditions (gap location, particle size, disk mass, and α viscosity) under which gaps in the outer disk efficiently block the inward drift of icy solids. We find that inner disk vapor enrichment is highly sensitive to the location of a disk gap, yielding for each particle size a radial “sweet spot” that reduces the inner disk vapor enrichment to a minimum. For pebbles of 1–10 mm in size, which carry the most mass, this sweet spot is at 7–15 au, suggesting that inner gaps may have a key role in reducing ice delivery to the inner disk and may not allow the formation of Earths and super-Earths. This highlights the importance of observationally determining the presence and properties of inner gaps in disks. Finally, we argue that the inner water vapor abundance can be used as a proxy for estimating the pebble drift efficiency and mass flux entering the inner disk.

Why Do M Dwarfs Have More Transiting Planets?

Abstract We propose a planet formation scenario to explain the elevated occurrence rates of transiting planets around M dwarfs compared to Sun-like stars discovered by Kepler. We use a pebble drift and accretion model to simulate the growth of planet cores inside and outside of the snow line. A smaller pebble size interior to the snow line delays the growth of super-Earths, allowing giant planet cores in the outer disk to form first. When those giant planets reach pebble isolation mass they cut off the flow of pebbles to the inner disk and prevent the formation of close-in super-Earths. We apply this model to stars with masses between 0.1 and 2 M ⊙ and for a range of initial disk masses. We find that the masses of hot super-Earths and of cold giant planets are anticorrelated. The fraction of our simulations that form hot super-Earths is higher around lower-mass stars and matches the exoplanet occurrence rates from Kepler. The fraction of simulations forming cold giant planets is consistent with the stellar mass dependence from radial-velocity surveys. A key testable prediction of the pebble accretion hypothesis is that the occurrence rates of super-Earths should decrease again for M dwarfs near the substellar boundary like Trappist-1.

Formation of planetary systems by pebble accretion and migration

At least 30% of main sequence stars host planets with sizes of between 1 and 4 Earth radii and orbital periods of less than 100 days. We useN-body simulations including a model for gas-assisted pebble accretion and disk–planet tidal interaction to study the formation of super-Earth systems. We show that the integrated pebble mass reservoir creates a bifurcation between hot super-Earths or hot-Neptunes (≲15M⊕) and super-massive planetary cores potentially able to become gas giant planets (≳15M⊕). Simulations with moderate pebble fluxes grow multiple super-Earth-mass planets that migrate inwards and pile up at the inner edge of the disk forming long resonant chains. We follow the long-term dynamical evolution of these systems and use the period ratio distribution of observed planet-pairs to constrain our model. Up to ~95% of resonant chains become dynamically unstable after the gas disk dispersal, leading to a phase of late collisions that breaks the original resonant configurations. Our simulations naturally match observations when they produce a dominant fraction (≳95%) of unstable systems with a sprinkling (≲5%) of stable resonant chains (the Trappist-1 system represents one such example). Our results demonstrate that super-Earth systems are inherently multiple (N≥ 2) and that the observed excess of single-planet transits is a consequence of the mutual inclinations excited by the planet–planet instability. In simulations in which planetary seeds are initially distributed in the inner and outer disk, close-in super-Earths are systematically ice rich. This contrasts with the interpretation that most super-Earths are rocky based on bulk-density measurements of super-Earths and photo-evaporation modeling of their bimodal radius distribution. We investigate the conditions needed to form rocky super-Earths. The formation of rocky super-Earths requires special circumstances, such as far more efficient planetesimal formation well inside the snow line, or much faster planetary growth by pebble accretion in the inner disk. Intriguingly, the necessary conditions to match the bulk of hot super-Earths are at odds with the conditions needed to match the Solar System.

The growth of super-Earths

The conditions in the protoplanetary disk are determinant for the various planet formation mechanisms. We present a framework that combines self-consistent disk structures with the calculations of the growth rates of planetary embryos via pebble accretion, in order to study the formation of super-Earths. We first perform 2D hydrodynamical simulations of the inner disks, considering a grain size distribution with multiple chemical species and their corresponding size and composition dependent opacities. The resulting aspect ratios are almost constant with orbital distance, resulting in radially constant pebble isolation masses, the mass where pebble accretion stops. This supports the “peas-in-a-pod” constraint from the Kepler observations. The derived pebble sizes are used to calculate the growth rates of planetary embryos via pebble accretion. Disks with low levels of turbulence (expressed through the α-viscosity) and/or high dust fragmentation velocities allow larger particles, hence lead to lower pebble isolation masses, and the contrary. At the same time, small pebble sizes lead to low accretion rates. We find that there is a trade-off between the pebble isolation mass and the growth timescale; the best set of parameters is an α-viscosity of 10−3 and a dust fragmentation velocity of 10 m s−1, mainly for an initial gas surface density (at 1 AU) greater than 1000 g cm−2. A self-consistent treatment between the disk structures and the pebble sizes is thus of crucial importance for planet formation simulations.

Six transiting planets and a chain of Laplace resonances in TOI-178

Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this context, TOI-178 has been the subject of particular attention since the first TESS observations hinted at the possible presence of a near 2:3:3 resonant chain. Here we report the results of observations from CHEOPS, ESPRESSO, NGTS, and SPECULOOS with the aim of deciphering the peculiar orbital architecture of the system. We show that TOI-178 harbours at least six planets in the super-Earth to mini-Neptune regimes, with radii ranging from 1.152−0.070+0.073 to 2.87−0.13+0.14 Earth radii and periods of 1.91, 3.24, 6.56, 9.96, 15.23, and 20.71 days. All planets but the innermost one form a 2:4:6:9:12 chain of Laplace resonances, and the planetary densities show important variations from planet to planet, jumping from 1.02−0.23+0.28 to 0.177−0.061+0.055 times the Earth’s density between planets c and d. Using Bayesian interior structure retrieval models, we show that the amount of gas in the planets does not vary in a monotonous way, contrary to what one would expect from simple formation and evolution models and unlike other known systems in a chain of Laplace resonances. The brightness of TOI-178 (H = 8.76 mag, J = 9.37 mag, V = 11.95 mag) allows for a precise characterisation of its orbital architecture as well as of the physical nature of the six presently known transiting planets it harbours. The peculiar orbital configuration and the diversity in average density among the planets in the system will enable the study of interior planetary structures and atmospheric evolution, providing important clues on the formation of super-Earths and mini-Neptunes.

Rapid formation of super-Earths around low-mass stars

ABSTRACT NASA’s TESS mission is expected to discover hundreds of M dwarf planets. However, few studies focus on how planets form around low-mass stars. We aim to better characterize the formation process of M dwarf planets to fill this gap and aid in the interpretation of TESS results. We use ten sets of N-body planet formation simulations that vary in whether a gas disc is present, initial range of embryo semimajor axes, and initial solid surface density profile. Each simulation begins with 147 equal-mass embryos around a 0.2 solar mass star and runs for 100 Myr. We find that planets form rapidly, with most collisions occurring within the first 1 Myr. The presence of a gas disc reduces the final number of planets relative to a gas-free environment and causes planets to migrate inward. We find that roughly a quarter of planetary systems experience their final giant impact inside the gas disc, suggesting that some super-Earths may be able to reaccrete an extended gaseous envelope after their final giant impact, though these may be affected by additional processes such as photoevaporation. In addition, we find that the final distribution of planets does not retain a memory of the slope of the initial surface density profile, regardless of whether or not a gas disc is present. Thus, our results suggest that present-day observations are unlikely to provide sufficient information to accurately reverse-engineer the initial distribution of solids.

Composition of super-Earths, super-Mercuries, and their host stars

Because of their common origin, it was assumed that the composition of planet building blocks should, to a first order, correlate with stellar atmospheric composition, especially for refractory elements. In fact, information on the relative abundance of refractory and major rock-forming elements such as Fe, Mg, Si has been commonly used to improve interior estimates for terrestrial planets. Recently Adibekyan et al. (2021) presented evidence of a tight chemical link between rocky planets and their host stars. In this study we add six recently discovered exoplanets to the sample of Adibekyan et al and re-evaluate their findings in light of these new data. We confirm that i) iron-mass fraction of rocky exoplanets correlates (but not a 1:1 relationship) with the composition of their host stars, ii) on average the iron-mass fraction of planets is higher than that of the primordial f star iron , iii) super-Mercuries are formed in disks with high iron content. Based on these results we conclude that disk-chemistry and planet formation processes play an important role in the composition, formation, and evolution of super-Earths and super-Mercuries.

Linking planetary embryo formation to planetesimal formation

Context.The accretion of pebbles onto planetary cores has been widely studied in recent years and is found to be a highly effective mechanism for planetary growth. While most studies assume planetary cores as an initial condition in their simulation, the question of the manner, location, and time at which these cores form is often neglected.Aims.We study the effect of pebble accretion during the formation phase and subsequent evolution of planetary embryos in the early stages of circumstellar disk evolution. In doing so, we aim to quantify the timescales and local dependence of planetary embryo formation based on the solid evolution of the disk.Methods.We connected a one-dimensional two-population model for solid evolution and pebble-flux-regulated planetesimal formation to theN-body code LIPAD. We focused on the growth of planetesimals with an initial size of 100 km in diameter by planetesimal collisions and pebble accretion for the first one million years of a viscously evolving disk. We compared 18 differentN-body simulations in which we varied the total planetesimal mass after one million years, the surface density profile of the planetesimal disk, the radial pebble flux, and the possibility of pebble accretion.Results.Pebble accretion leads to the formation of fewer but substantially more massive embryos. The area of possible embryo formation is weakly affected by the accretion of pebbles, and the innermost embryos tend to form slightly earlier than in simulations in which pebble accretion is neglected.Conclusions.Pebble accretion strongly enhances the formation of super-Earths in the terrestrial planet region, but it does not enhance the formation of embryos at larger distances.

Carbon, isotopic ratio 12C/13C, and nitrogen in solar twins: constraints for the chemical evolution of the local disc

ABSTRACT Abundances of light elements in dwarf stars of different ages are important constraints for stellar yields, Galactic chemical evolution and exoplanet chemical composition studies. We have measured C and N abundances and 12C/13C ratios for a sample of 63 solar twins spanning a wide range in age, based on spectral synthesis of a comprehensive list of CH A-X and CN B-X features using HARPS spectra. The analysis of 55 thin disc solar twins confirms the dependences of [C/Fe] and [N/Fe] on [Fe/H]. [N/Fe] is investigated as a function of [Fe/H] and age for the first time for these stars. Our derived correlation [C/Fe]–age agrees with works for solar-type stars and solar twins, but the [N/Fe]–age correlation does not. The relations [C,N/Fe]–[Fe/H] and [C,N/Fe]–age for the solar twins lay under solar. 12C/13C is found correlated with [Fe/H] and seems to have decreased along the evolution of the local thin disc. Predictions from chemical evolution models for the solar vicinity corroborate the relations [C,N/Fe]–[Fe/H], 12C/13C–age, and [N/O]–[O/H], but do not for the 12C/13C–[Fe/H] and [C/O]–[O/H] relations. The N/O ratio in the Sun is placed at the high end of the homogeneous distribution of solar twins, which suggests uniformity in the N–O budget for the formation of icy planetesimals, watery superearths, and giant planets. C and N had different nucleosynthetic origins along the thin disc evolution, as shown by the relations of [C/N], [C/O], and [N/O] against [O/H] and age. [C/N] and [C/O] are particularly observed increasing in time for solar twins younger than the Sun.

Unified Simulations of Planetary Formation and Atmospheric Evolution. II. Rapid Disk Clearing by Photoevaporation Yields Low-mass Super-Earth Atmospheres

Abstract Super-Earths possess low-mass H2/He atmospheres (typically less than 10% by mass). However, the origins of super-Earth atmospheres have not yet been ascertained. We investigate the role of rapid disk clearing by photoevaporation during the formation of super-Earths and their atmospheres. We perform unified simulations of super-Earth formation and atmospheric evolution in evolving disks that consider both photoevaporative winds and magnetically driven disk winds. For the growth mode of planetary cores, we consider two cases in which planetary embryos grow with and without pebble accretion. Our main findings are summarized as follows. (i) The time span of atmospheric accretion is shortened by rapid disk dissipation due to photoevaporation, which prevents super-Earth cores from accreting massive atmospheres. (ii) Even if planetary cores grow rapidly by embryo accretion in the case without pebble accretion, the onset of runaway gas accretion is delayed because the isolation mass for embryo accretion is small. Together with rapid disk clearing, the accretion of massive atmospheres can be avoided. (iii) After rapid disk clearing, a number of high-eccentricity embryos can remain in outer orbits. Thereafter, such embryos may collide with the super-Earths, leading to efficient impact erosion of accreted atmospheres. Therefore, we find that super-Earths with low-mass H2/He atmospheres are naturally produced by N-body simulations that consider realistic disk evolution.

Formation of planetary populations – III. Core composition and atmospheric evaporation

ABSTRACT The exoplanet mass–radius diagram reveals that super-Earths display a wide range of radii, and therefore mean densities, at a given mass. Using planet population synthesis models, we explore the key physical factors that shape this distribution: planets’ solid core compositions, and their atmospheric structure. For the former, we use equilibrium disc chemistry models to track accreted minerals on to planetary cores throughout the formation. For the latter, we track gas accretion during the formation and consider photoevaporation-driven atmospheric mass-loss to determine what portion of accreted gas escapes after the disc phase. We find that atmospheric stripping of Neptunes and sub-Saturns at small orbital radii (≲0.1 au) plays a key role in the formation of short-period super-Earths. Core compositions are strongly influenced by the trap in which they formed. We also find a separation between Earth-like planet compositions at small orbital radii ≲0.5 au and ice-rich planets (up to 50 per cent by mass) at larger orbits ∼1 au. This corresponds well with the Earth-like mean densities inferred from the observed position of the low-mass planet radius valley at small orbital periods. Our model produces planet radii comparable to observations at masses ∼1–3 M⊕. At larger masses, planets’ accreted gas significantly increases their radii to be larger than most of the observed data. While photoevaporation, affecting planets at small orbital radii ≲0.1 au, reduces a subset of these planets’ radii and improves our comparison, most planets in our computed populations are unaffected due to low-far ultraviolet fluxes as they form at larger separations.

The origins of nearly coplanar, non-resonant systems of close-in super-Earths

ABSTRACT Some systems of close-in ‘super-Earths’ contain five or more planets on non-resonant but compact and nearly coplanar orbits. The Kepler-11 system is an iconic representative of this class of system. It is challenging to explain their origins given that planet–disc interactions are thought to be essential to maintain such a high degree of coplanarity, yet these same interactions invariably cause planets to migrate into chains of mean motion resonances. Here, we mine a large data set of dynamical simulations of super-Earth formation by migration. These simulations match the observed period ratio distribution as long as the vast majority of planet pairs in resonance become dynamically unstable. When instabilities take place resonances are broken during a late phase of giant impacts, and typical surviving systems have planet pairs with significant mutual orbital inclinations. However, a subset of our unstable simulations matches the Kepler-11 system in terms of coplanarity, compactness, planet-multiplicity, and non-resonant state. This subset has dynamical instability phases typically much shorter than ordinary systems. Unstable systems may keep a high degree of coplanarity post-instability if planets collide at very low orbital inclinations (≲1○) or if collisions promote efficient damping of orbital inclinations. If planetary scattering during the instability takes place at low orbital inclinations (i ≲ 1○), orbital inclinations are barely increased by encounters before planets collide. When planetary scattering pumps orbital inclinations to higher values (≳1○) planets tend to collide at higher mutual orbital inclinations, but depending on the geometry of collisions mergers’ orbital inclinations may be efficiently damped. Each of these formation pathways can produce analogues to the Kepler-11 system.

Unified Simulations of Planetary Formation and Atmospheric Evolution: Effects of Pebble Accretion, Giant Impacts, and Stellar Irradiation on Super-Earth Formation

Abstract A substantial number of super-Earths have been discovered, and atmospheres of transiting super-Earths have also been observed by transmission spectroscopy. Several lines of observational evidence indicate that most super-Earths do not possess massive H2/He atmospheres. However, accretion and retention of less massive atmospheres on super-Earths challenge planet formation theory. We consider the following three mechanisms: (i) envelope heating by pebble accretion, (ii) mass loss during giant impacts, and (iii) atmospheric loss by stellar X-ray and EUV photoevaporation. We investigate whether these mechanisms influence the amount of the atmospheres that form around super-Earths. We develop a code combining an N-body simulation of pebble-driven planetary formation and an atmospheric evolution simulation. We demonstrate that the observed orbital properties of super-Earths are well reproduced by the results of our simulations. However, (i) heating by pebble accretion ceases prior to disk dispersal, (ii) the frequency of giant impact events is too low to sculpt massive atmospheres, and (iii) many super-Earths having H2/He atmospheres of ≳10 wt% survive against stellar irradiation for 1 Gyr. Therefore, it is likely that other mechanisms, such as suppression of gas accretion, are required to explain less massive atmospheres (≲10 wt%) of super-Earths.

Forming Diverse Super-Earth Systems In Situ

Abstract Super-Earths and mini-Neptunes exhibit great diversity in their compositional and orbital properties. Their bulk densities span a large range, from those dense enough to be purely rocky to those needing a substantial contribution from volatiles to their volumes. Their orbital configurations range from compact, circular multitransiting systems like Kepler-11 to systems like our solar system’s terrestrial planets with wider spacings and modest but significant eccentricities and mutual inclinations. Here we investigate whether a continuum of formation conditions resulting from variation in the amount of solids available in the inner disk can account for the diversity of orbital and compositional properties observed for super-Earths, including the apparent dichotomy between single and multitransiting systems. We simulate in situ formation of super-Earths via giant impacts and compare to the observed Kepler sample. We find that intrinsic variations among disks in the amount of solids available for in situ formation can account for the orbital and compositional diversity observed among Kepler’s transiting planets. Our simulations can account for the planets’ distributions of orbital period ratios, transit duration ratios, and transit multiplicity; higher eccentricities for planets in single transiting systems than for those in multitransiting systems; smaller eccentricities for larger planets; scatter in the mass–radius relation, including lower densities for planets with masses measured with transit timing variations instead of with radial velocity; and similarities in planets’ sizes and spacings within each system. Our findings support the theory that variation among super-Earth and mini-Neptune properties is primarily locked in by different in situ formation conditions, rather than arising stochastically through subsequent evolution.