Constraining disk evolution prescriptions of planet population synthesis models with observed disk masses and accretion rates

While planets are commonly discovered around main-sequence stars, the processes leading to their formation are still far from being understood. Current planet population synthesis models, which aim to describe the planet formation process from the protoplanetary disk phase to the time exoplanets are observed, rely on prescriptions for the underlying properties of protoplanetary disks where planets form and evolve. The recent development in measuring disk masses and disk-star interaction properties, i.e., mass accretion rates, in large samples of young stellar objects demand a more careful comparison between the models and the data. We performed an initial critical assessment of the assumptions made by planet synthesis population models by looking at the relation between mass accretion rates and disk masses in the models and in the currently available data. We find that the currently used disk models predict mass accretion rate in line with what is measured, but with a much lower spread of values than observed. This difference is mainly because the models have a smaller spread of viscous timescales than what is needed to reproduce the observations. We also find an overabundance of weakly accreting disks in the models where giant planets have formed with respect to observations of typical disks. We suggest that either fewer giant planets have formed in reality or that the prescription for planet accretion predicts accretion on the planets that is too high. Finally, the comparison of the properties of transition disks with large cavities confirms that in many of these objects the observed accretion rates are higher than those predicted by the models. On the other hand, PDS70, a transition disk with two detected giant planets in the cavity, shows mass accretion rates well in line with model predictions

The New Generation Planetary Population Synthesis (NGPPS). IV. Planetary systems around low-mass stars

Context. Previous theoretical works on planet formation around low-mass stars have often been limited to large planets and individual systems. As current surveys routinely detect planets down to terrestrial size in these systems, models have shifted toward a more holistic approach that reflects their diverse architectures. Aims. Here, we investigate planet formation around low-mass stars and identify differences in the statistical distribution of modeled planets. We compare the synthetic planet populations to observed exoplanets and we discuss the identified trends. Methods. We used the Generation III Bern global model of planet formation and evolution to calculate synthetic populations, while varying the central star from Solar-like stars to ultra-late M dwarfs. This model includes planetary migration, N-body interactions between embryos, accretion of planetesimals and gas, and the long-term contraction and loss of the gaseous atmospheres. Results. We find that temperate, Earth-sized planets are most frequent around early M dwarfs (0.3 M⊙–0.5 M⊙) and that they are more rare for Solar-type stars and late M dwarfs. The planetary mass distribution does not linearly scale with the disk mass. The reason behind this is attributed to the emergence of giant planets for M⋆ ≥ 0.5 M⊙, which leads to the ejection of smaller planets. Given a linear scaling of the disk mass with stellar mass, the formation of Earth-like planets is limited by the available amount of solids for ultra-late M dwarfs. For M⋆ ≥ 0.3 M⊙, however, there is sufficient mass in the majority of systems, leading to a similar amount of Exo-Earths going from M to G dwarfs. In contrast, the number of super-Earths and larger planets increases monotonically with stellar mass. We further identify a regime of disk parameters that reproduces observed M-dwarf systems such as TRAPPIST-1. However, giant planets around late M dwarfs, such as GJ 3512b, only form when type I migration is substantially reduced. Conclusions. We are able to quantify the stellar mass dependence of multi-planet systems using global simulations of planet formation and evolution. The results fare well in comparison to current observational data and predict trends that can be tested with future observations

The New Generation Planetary Population Synthesis (NGPPS) VI. Introducing KOBE: Kepler Observes Bern Exoplanets

Context. Observations of exoplanets indicate the existence of several correlations in the architecture of planetary systems. Exoplanets within a system tend to be of similar size and mass, evenly spaced, and are often ordered in size and mass. Small planets are frequently packed in tight configurations, while large planets often have wider orbital spacing. Together, these correlations are called the peas in a pod trends in the architecture of planetary systems. Aims. In this paper these trends are investigated in theoretically simulated planetary systems and compared with observations. Whether these correlations emerge from astrophysical processes or the detection biases of the transit method is examined. Methods. Synthetic planetary system were simulated using the Generation III Bern Model. KOBE, a new computer code, simulates the geometrical limitations of the transit method and applies the detection biases and completeness of the Kepler survey. This allows simulated planetary systems to be compared with observations. Results. The architecture of synthetic planetary systems, observed via KOBE, show the peas in a pod trends in good agreement with observations. These correlations are also present in the theoretical underlying population, from the Bern Model, indicating that these trends are probably of astrophysical origin. Conclusions. The physical processes involved in planet formation are responsible for the emergence of evenly spaced planets with similar sizes and masses. The size–mass similarity trends are primordial and originate from the oligarchic growth of protoplanetary embryos and the uniform growth of planets at early times. Later stages in planet formation allows planets within a system to grow at different rates, thereby decreasing these correlations. The spacing and packing correlations are absent at early times and arise from dynamical interactions

The New Generation Planetary Population Synthesis (NGPPS). V. Predetermination of planet types in global core accretion models

Context. State-of-the-art planet formation models are now capable of accounting for the full spectrum of known planet types. This comes at the cost of an increasing complexity of the models, which calls into question whether established links between their initial conditions and the calculated planetary observables are preserved. Aims. In this paper, we take a data-driven approach to investigate the relations between clusters of synthetic planets with similar properties and their formation history. Methods. We trained a Gaussian mixture model on typical exoplanet observables computed by a global model of planet formation to identify clusters of similar planets. We then traced back the formation histories of the planets associated with them and pinpointed their differences. Using the cluster affiliation as labels, we trained a random forest classifier to predict planet species from properties of the originating protoplanetary disk. Results. Without presupposing any planet types, we identified four distinct classes in our synthetic population. They roughly correspond to the observed populations of (sub-)Neptunes, giant planets, and (super-)Earths, plus an additional unobserved class we denote as “icy cores”. These groups emerge already within the first 0.1 Myr of the formation phase and are predicted from disk properties with an overall accuracy of >90%. The most reliable predictors are the initial orbital distance of planetary nuclei and the total planetesimal mass available. Giant planets form only in a particular region of this parameter space that is in agreement with purely analytical predictions. Including N-body interactions between the planets decreases the predictability, especially for sub-Neptunes that frequently undergo giant collisions and turn into super-Earths. Conclusions. The processes covered by current core accretion models of planet formation are largely predictable and reproduce the known demographic features in the exoplanet population. The impact of gravitational interactions highlights the need for N-body integrators for realistic predictions of systems of low-mass planets

The Effect of Inefficient Accretion on Planetary Differentiation

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 non-accreted 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 Earth's 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.Comment: 21 pages, 2 tables, 7 figures. Published open access on PSJ: https://iopscience.iop.org/article/10.3847/PSJ/abf0a

Direct imaging and spectroscopy of exoplanets with the ELT/HARMONI high-contrast module

Combining high-contrast imaging with medium-resolution spectroscopy has been shown to significantly boost the direct detection of exoplanets. HARMONI, one of the first-light instruments to be mounted on ESO’s future extremely large telescope (ELT), will be equipped with a single-conjugated adaptive optics system to reach the diffraction limit of the ELT in the H and K bands, a high-contrast module dedicated to exoplanet imaging, and a medium-resolution (up to R = 17 000) optical and near-infrared integral field spectrograph. When combined, these systems will provide unprecedented contrast limits at separations between 50 and 400 mas. This paper is aimed at estimating the capabilities of the HARMONI high-contrast module for the direct detection of young giant exoplanets. We use an end-to-end model of the instrument to simulate high-contrast observations performed with HARMONI, based on realistic observing scenarios and conditions. We then analyze these data with the so-called “molecule mapping” technique combined with a matched-filter approach in order to disentangle companions from the host star and tellurics and to increase the signal-to-noise ratio (S/N) of the planetary signal. We detected planets above 5σ at contrasts up to 16 mag and separations down to 75 mas in several spectral configurations of the instrument. We show that molecule mapping allows for the detection of companions up to 2.5 mag fainter compared to state-of-the-art high-contrast imaging techniques based on angular differential imaging. We also demonstrate that the performance is not strongly affected by the spectral type of the host star and we show that we are able to reach close sensitivities for the best three quartiles of observing conditions at Armazones, which means that HARMONI could be used in near-critical observations during 60 to 70% of telescope time at the ELT. Finally, we simulated planets from population synthesis models to further explore the parameter space that HARMONI and its high-contrast module will open up and compare this to the current high-contrast instrumentation

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’1,2 (a region in mass–radius space that contains few planets) have proved to be particularly valuable in this regard. These planets include HD149026b3, which is thought to have an unusually massive core, and recent discoveries such as LTT9779b4 and NGTS-4b5, 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.6+2.7 Earth masses and a density of 5.2−0.8+0.7 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.9+0.8 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 formation6. 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

A giant exoplanet orbiting a very-low-mass star challenges planet formation models

Surveys have shown that super-Earth and Neptune-mass exoplanets are more frequent than gas giants around low-mass stars, as predicted by the core accretion theory of planet formation. We report the discovery of a giant planet around the very-low-mass star GJ 3512, as determined by optical and near-infrared radial-velocity observations. The planet has a minimum mass of 0.46 Jupiter masses, very high for such a small host star, and an eccentric 204-day orbit. Dynamical models show that the high eccentricity is most likely due to planet-planet interactions. We use simulations to demonstrate that the GJ 3512 planetary system challenges generally accepted formation theories, and that it puts constraints on the planet accretion and migration rates. Disk instabilities may be more efficient in forming planets than previously thought

Collision chains among the terrestrial planets. III. formation of the moon

In the canonical model of Moon formation, a Mars-sized protoplanet “Theia” collides with proto-Earth at close to their mutual escape velocity vesc and a common impact angle ∼45°. The “graze-and-merge” collision strands a fraction of Theia’s mantle into orbit, while Earth accretes most of Theia and its momentum. Simulations show that this produces a hot, high angular momentum, silicate-dominated protolunar system, in substantial agreement with lunar geology, geochemistry, and dynamics. However, a Moon that derives mostly from Theia’s mantle, as angular momentum dictates, is challenged by the fact that O, Ti, Cr, radiogenic W, and other elements are indistinguishable in Earth and lunar rocks. Moreover, the model requires an improbably low initial velocity. Here we develop a scenario for Moon formation that begins with a somewhat faster collision, when proto-Theia impacts proto-Earth at ∼ 1.2vesc, also around ∼45°. Instead of merging, the bodies come into violent contact for a half hour and their major components escape, a “hit-and-run” collision. N-body evolutions show that the “runner” often returns ∼0.1–1 Myr later for a second giant impact, closer to vesc; this produces a postimpact disk of ∼2–3 lunar masses in smoothed particle hydrodynamics simulations, with angular momentum comparable to canonical scenarios. The disk ends up substantially inclined, in most cases, because the terminal collision is randomly oriented to the first. Moreover, proto-Earth contributions to the protolunar disk are enhanced by the compounded mixing and greater energy of a collision chain. © 2021. The Author(s). Published by the American Astronomical Society.Open access journalThis item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at [email protected]

The New Generation Planetary Population Synthesis (NGPPS). III. Warm super-Earths and cold Jupiters: A weak occurrence correlation, but with a strong architecture-composition link

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. 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. 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. 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.Comment: 29 pages, 28 figures. Accepted for publication in Astronomy & Astrophysic