Planet Formation Process Research Articles

Characterization of Ring Substructures in the Protoplanetary Disk of HD 169142 from Multiwavelength Atacama Large Millimeter/submillimeter Array Observations

Abstract We present a detailed multiwavelength characterization of the multi-ring disk of HD 169142. We report new Atacama Large Millimeter/submillimeter Array (ALMA) observations at 3 mm and analyze them together with archival 0.89 and 1.3 mm data. Our observations resolve three out of the four rings in the disk previously seen in high-resolution ALMA data. A simple parametric model is used to estimate the radial profile of the dust optical depth, temperature, density, and particle size distribution. We find that the multiple ring features of the disk are produced by annular accumulations of large particles, probably associated with gas pressure bumps. Our model indicates that the maximum dust grain size in the rings is ∼1 cm, with slightly flatter power-law size distributions than the interstellar medium-like size distribution (p ∼ 3.5) found in the gaps. In particular, the inner ring (∼26 au) is associated with a strong and narrow buildup of dust particles that could harbor the necessary conditions to trigger the streaming instability. According to our analysis, the snowlines of the most important volatiles do not coincide with the observed substructures. We explore different ring formation mechanisms and find that planet–disk interactions are the most likely scenario to explain the main features of HD 169142. Overall, our multiwavelength analysis provides some of the first unambiguous evidence of the presence of radial dust traps in the rings of HD 169142. A similar analysis in a larger sample of disks could provide key insights on the impact that disk substructures have on the dust evolution and planet formation processes.

Evidence for an additional planet in the β Pictoris system

With its imaged debris disk of dust, its evaporating exocomets, and an imaged giant planet, the young (~23 Myr) β Pictoris system is a unique proxy for detailed studies of planet formation processes as well as planet–disk interactions. Here, we study ten years of European Southern Observatory/High Accuracy Radial Velocity Planet Searcher (HARPS) high-resolution spectroscopic data of β Pictoris. After removing the radial velocity (RV) signals arising from the δ Scuti pulsations of the star, a ~1,200-d periodic signal remains, which, within our current knowledge, we can only attribute to a second planet in the system. The β Pic c mass is about nine times the mass of Jupiter; it orbits at ~2.7 au on an eccentric (e ~ 0.24) orbit. More RV data are needed to obtain more precise estimates of the properties of β Pic c. The current modelling of the planet’s properties and the dynamic of the whole system has to be reinvestigated in light of this detection. Radial velocity data of the young β Pictoris system acquired by HARPS and spanning 15 years show evidence of β Pic c, a gas giant of ~9 Jupiter masses orbiting on an eccentric orbit at ~2.4 au from the star, near the theoretical snowline. Both β Pic b and c, located close to the star, may have formed in situ by core accretion.

Determination of the isotopic composition of tungsten using MC-ICP-MS

The processes of planetary accretion or formation of the Earth and other celestial objects can be studied by using the 182Hf-182W chronometer which requires precise measurements of tungsten isotope ratios. Many comparative measurements for the isotopic composition of tungsten have been performed using either thermal ionization mass spectrometry (TIMS) or multi-collector inductively coupled plasma mass spectrometry (MC-ICP-MS). Yet, calibrated measurements of tungsten isotope ratios, and, in turn, isotopic abundances and atomic weight, are still lacking. In this study, we report the first independent measurements of all tungsten isotope ratios in five commercial tungsten reagents, including the new NRC candidate isotopic reference material WOLF-1 by MC-ICP-MS with use of the-state-of-the-art optimized regression mass bias correction model and NIST SRM 989 isotopic rhenium as calibrator.

Pairwise tidal equilibrium states and the architecture of extrasolar planetary systems

ABSTRACT Current observations indicate that the planet formation process often produces multiple planet systems with nearly circular orbits, regular spacing, a narrow range of inclination angles, and similar planetary masses of order mp ∼ 10 M⊕. Motivated by the observational sample, this paper determines the tidal equilibrium states for this class of extrasolar planetary systems. We start by considering two-planet systems with fixed orbital spacing and variable mass ratios. The basic conjecture explored in this paper is that the planet formation process will act to distribute planetary masses in order to achieve a minimum energy state. The resulting minimum energy configuration – subject to the constraint of constant angular momentum – corresponds to circular orbits confined to a plane, with nearly equal planetary masses (as observed). We then generalize the treatment to include multiple planet systems, where each adjacent pair of planets attains its (local) tidal equilibrium state. The properties of observed planetary systems are close to those expected from this pairwise equilibrium configuration. In contrast, observed systems do not reside in a global minimum energy state. Both the equilibrium states of this paper and observed multiplanet systems, with planets of nearly equal mass on regularly spaced orbits, have an effective surface density of the form σ ∝ r−2, much steeper than most disc models.

Pebble-driven planet formation for TRAPPIST-1 and other compact systems

Recently, seven Earth-sized planets were discovered around the M-dwarf star TRAPPIST-1. Thanks to transit-timing variations, the masses and therefore the bulk densities of the planets have been constrained, suggesting that all TRAPPIST-1 planets are consistent with water mass fractions on the order of 10%. These water fractions, as well as the similar planet masses within the system, constitute strong constraints on the origins of the TRAPPIST-1 system. In a previous work, we outlined a pebble-driven formation scenario. In this paper we investigate this formation scenario in more detail. We used a Lagrangian smooth-particle method to model the growth and drift of pebbles and the conversion of pebbles to planetesimals through the streaming instability. We used the N-body code MERCURY to follow the composition of planetesimals as they grow into protoplanets by merging and accreting pebbles. This code is adapted to account for pebble accretion, type-I migration, and gas drag. In this way, we modelled the entire planet formation process (pertaining to planet masses and compositions, not dynamical configuration). We find that planetesimals form in a single, early phase of streaming instability. The initially narrow annulus of planetesimals outside the snowline quickly broadens due to scattering. Our simulation results confirm that this formation pathway indeed leads to similarly-sized planets and is highly efficient in turning pebbles into planets. Our results suggest that the innermost planets in the TRAPPIST-1 system grew mostly by planetesimal accretion at an early time, whereas the outermost planets were initially scattered outwards and grew mostly by pebble accretion. The water content of planets resulting from our simulations is on the order of 10%, and our results predict a “V-shaped” trend in the planet water fraction with orbital distance: from relatively high (innermost planets) to relatively low (intermediate planets) to relatively high (outermost planets).

A sub-Neptune exoplanet with a low-metallicity methane-depleted atmosphere and Mie-scattering clouds

With no analogues in the Solar System, the discovery of thousands of exoplanets with masses and radii intermediate between Earth and Neptune was one of the big surprises of exoplanet science. These super-Earths and sub-Neptunes likely represent the most common outcome of planet formation. Mass and radius measurements indicate a diversity in bulk composition much wider than for gas giants; however, direct spectroscopic detections of molecular absorption and constraints on the gas mixing ratios have largely remained limited to planets more massive than Neptune. Here, we analyze a combined Hubble/Spitzer Space Telescope dataset of 12 transits and 20 eclipses of the sub-Neptune GJ 3470 b, whose mass of 12.6 $M_\oplus$ places it near the half-way point between previously studied exo-Neptunes (22-23 $M_\oplus$) and exoplanets known to have rocky densities (7 $M_\oplus$). Obtained over many years, our data set provides a robust detection of water absorption (>5$\sigma$) and a thermal emission detection from the lowest irradiated planet to date. We reveal a low-metallicity, hydrogen-dominated atmosphere similar to a gas giant, but strongly depleted in methane gas. The low, near-solar metallicity (O/H=0.2-18) sets important constraints on the potential planet formation processes at low masses as well as the subsequent accretion of solids. The low methane abundance indicates that methane is destroyed much more efficiently than previously predicted, suggesting that the CH$_4$/CO transition curve has to be revisited for close-in planets. Finally, we also find a sharp drop in the cloud opacity at 2-3 $\mu$m characteristic of Mie scattering, which enables narrow constraints on the cloud particle size and makes GJ 3470b a keystone target for mid-IR characterization with JWST.

Fluorine and chlorine fractionation during magma ocean crystallization: Constraints on the origin of the non-chondritic F/Cl ratio of the Earth

Previous studies have reported that the relative abundances of volatile elements in the silicate Earth are non-chondritic. The abundance and distribution of volatile elements in terrestrial planets would have been predominantly controlled by planetary formation processes, including core-mantle separation, magma ocean crystallization, and volatility-dependent high-temperature fractionation. Thus, the current abundance patterns of volatile elements in the silicate fraction of terrestrial planets are the key to understanding the accretional history of terrestrial volatiles and the chemical differentiation of terrestrial planets. Although the origin of the non-chondritic ratios of volatile elements in terrestrial planets has been previously studied, it is still a matter of debate. In this study, we focused on the super-chondritic F/Cl ratio of the bulk silicate Earth and experimentally investigated the silicate mineral-melt partitioning of fluorine and chlorine at pressures from 18 GPa to 25 GPa. Our experimental results show that fluorine is moderately compatible with mantle minerals, whereas chlorine is highly incompatible. These results support the formation of a solid mantle with high F/Cl ratios, and a residual magma ocean and steam atmosphere with low F/Cl ratios during magma ocean crystallization. Thus, the F/Cl ratio in the residual solid parts of terrestrial planets would have become relatively enriched following escape of volatile elements from the planetary surface into outer space. This model is consistent with the collisional erosion hypothesis of primordial crusts and atmospheres, and current observations on the abundance and distribution of terrestrial fluorine and chlorine.

Explosion gravitation field algorithm with dust sampling for unconstrained optimization

Gravitation Field Algorithm (GFA) is a novel optimization algorithm derived from the Solar Nebular Disk Model (SNDM) in astronomy and inspired by the formation process of planets. Although it has achieved good performance when solving many unconstrained optimization problems, which demonstrated its promising application potential in many real-world problems, GFA still has much room for improvement, especially when it comes to the accuracy and efficiency of the algorithm.In this research, an improved GFA algorithm called Explosion Gravitation Field Algorithm (EGFA) is proposed for unconstrained optimization problems, with the introduction of two strategies: Dust Sampling (DS) and Explosion Operation. The task of DS is to locate the space that contains the optimal solution(s) by initializing the dust population randomly in the search space; while the Explosion Operator is to improve the accuracy of solutions and decrease the probability of the algorithm falling into local optima by generating the new population around the center dust to replace the original population.A comparison of experimental results on six classical unconstrained benchmark problems with different dimensions demonstrates that the proposed EGFA outperforms the original GFA and several classical metaheuristic optimization algorithms, such as Genetic Algorithm (GA) and Particle Swarm Optimization (PSO), in terms of accuracy and efficiency in lower dimensions. Additionally, the comparison of results on three real datasets indicate that EGFA performs better than the original GFA and k-means for solving clustering problems.

Thermal torque effects on the migration of growing low-mass planets

As planets grow the exchange of angular momentum with the gaseous component of the protoplanetary disc produces a net torque resulting in a variation of the semi-major axis of the planet. For low-mass planets not able to open a gap in the gaseous disc this regime is known as type I migration. Pioneer works studied this mechanism in isothermal discs finding fast inward type I migration rates that were unable to reproduce the observed properties of extrasolar planets. In the last years, several improvements have been made in order to extend the study of type I migration rates to non-isothermal discs. Moreover, it was recently shown that if the planet's luminosity due to solid accretion is taken into account, inward migration could be slowed down and even reversed. In this work, we study the planet formation process incorporating, and comparing, updated type I migration rates for non-isothermal discs and the role of planet's luminosity over such rates. We find that the latter can have important effects on planetary evolution, producing a significant outward migration for the growing planets.

Multiple Populations of Extrasolar Gas Giants

Abstract There are two planetary formation scenarios: core accretion and gravitational disk instability. Based on the fact that gaseous objects are preferentially observed around metal-rich host stars, most extrasolar gaseous objects discovered to date are thought to have been formed by core accretion. Here, we present 569 samples of gaseous planets and brown dwarfs found in 485 planetary systems that span three mass regimes with boundary values at 4 and 25 Jupiter-mass masses through performing cluster analyses of these samples regarding the host-star metallicity, after minimizing the impact of the selection effect of radial-velocity measurement on the cluster analysis. The larger mass is thought to be the upper mass limit of the objects that were formed during the planetary formation processes. In contrast, the lower mass limit appears to reflect the difference between planetary formation processes around early-type and G-type stars; disk instability plays a greater role in the planetary formation process around early-type stars than that around G-type stars. Populations with masses between 4 and 25 Jupiter masses that orbit early-type stars comprise planets formed not only via the core-accretion process but also via gravitational disk instability because the population preferentially orbits metal-poor stars or is independent of the host-star metallicity. Therefore, it is essential to have a hybrid scenario for the planetary formation of the diverse systems.

Formation of short-period planets by disc migration

ABSTRACT Protoplanetary discs are thought to be truncated at orbital periods of around 10 d. Therefore, the origin of rocky short-period planets with P < 10 d is a puzzle. We propose that many of these planets may form through the Type-I migration of planets locked into a chain of mutual mean motion resonances. We ran N-body simulations of planetary embryos embedded in a protoplanetary disc. The embryos experienced gravitational scatterings, collisions, disc torques, and dampening of orbital eccentricity and inclination. We then modelled Kepler observations of these planets using a forward model of both the transit probability and the detection efficiency of the Kepler pipeline. We found that planets become locked into long chains of mean motion resonances that migrate in unison. When the chain reaches the edge of the disc, the inner planets are pushed past the edge due to the disc torques acting on the planets farther out in the chain. Our simulated systems successfully reproduce the observed period distribution of short-period Kepler planets between 1 and 2 R⊕. However, we obtain fewer closely packed short-period planets than in the Kepler sample. Our results provide valuable insight into the planet formation process, and suggests that resonance locks, migration, and dynamical instabilities play important roles in the formation and evolution of close-in small exoplanets.

Dust in debris disk: Observations and laboratory experiments

AbstractDebris disks are the natural by-products of the star and planet formation processes. Since the 1980’s several thousands of debris disks have been detected, and the presence of a disk is inferred by the detection of excess emission over the photospheric emission. This thermal emission arises from (micron-sized or slightly bigger) dust grains heated by the central star. However, in the vast majority of cases, these observations are not spatially resolving the radial distribution of the dust, resulting in strong degeneracies in the modeling approach (radial distance vs minimum grain size mostly). Therefore the properties of the small dust grains remained largely unconstrained until the arrival of high angular resolution instruments, especially at optical and near-infrared wavelengths. In these proceeding some of the main results are presented that have been obtained over the past few years on the properties of small dust grains in debris disks, and it is discussed how laboratory experiments contributed to put those results in context.

Connecting substellar and stellar formation: the role of the host star’s metallicity

Context.Most of our current understanding of the planet formation mechanism is based on the planet metallicity correlation derived mostly from solar-type stars harbouring gas-giant planets.Aims.To achieve a more extensive grasp on the substellar formation process, we aim to analyse in terms of their metallicity a diverse sample of stars (in terms of mass and spectral type) covering the whole range of possible outcomes of the planet formation process (from planetesimals to brown dwarfs and low-mass binaries).Methods.Our methodology is based on the use of high-precision stellar parameters derived by our own group in previous works from high-resolution spectra by using the iron ionisation and equilibrium conditions. All values were derived in an homogeneous way, except for the M dwarfs where a methodology based on the use of pseudo equivalent widths of spectral features was used.Results.Our results show that as the mass of the substellar companion increases the metallicity of the host star tends to lower values. The same trend is maintained when analysing stars with low-mass stellar companions and a tendency towards a wide range of host star’s metallicity is found for systems with low-mass planets. We also confirm that more massive planets tend to orbit around more massive stars.Conclusions.The core-accretion formation mechanism for planet formation achieves its maximum efficiency for planets with masses in the range 0.2–2MJup. Substellar objects with higher masses have higher probabilities of being formed as stars. Low-mass planets and planetesimals might be formed by core-accretion even around low-metallicity stars.

Atmospheres on Nonsynchronized Eccentric-tilted Exoplanets. I. Dynamical Regimes

Abstract Relatively long-period, nonsynchronized planets—such as warm Jupiters—potentially retain the primordial rotation, eccentricity, and obliquity that might encapsulate information on planetary climate and formation processes. To date, there has not been a systematic study on climate patterns on these planets that will significantly influence their observations. Here we investigate the atmospheric dynamics of nonsynchronized, fast-rotating exoplanets across various radiative timescales, eccentricities, and obliquities using a shallow water model. The dynamical pattern can be demarcated into five regimes in terms of radiative timescale τ rad and obliquity θ. An atmosphere with τ rad shorter than a planetary day usually exhibits a strong day–night temperature contrast and a day-to-night flow pattern. In the intermediate τ rad regime between a planetary day and a year, the atmosphere is dominated by steady temperature and eastward jet patterns for θ ≤ 18° but shows a strong seasonal variation for θ ≥ 18° because the polar region undergoes an intense heating at around the summer solstice. If τ rad is larger than a year, seasonal variation is very weak. In this regime, eastward jets are developed for θ ≤ 54° and westward jets are developed for θ ≥ 54°. These dynamical regimes are also applicable to planets in eccentric orbits. The large effects of exoplanetary obliquities on circulation patterns might offer observational signatures, which will be investigated in Paper II of this study.

Resolved Imaging of the AR Puppis Circumbinary Disk*

Abstract Circumbinary disks are common around post-asymptotic giant branch (post-AGB) stars with a stellar companion on orbital timescales of a few 100 to few 1000 days. The presence of a disk is usually inferred from the system’s spectral energy distribution and confirmed, for a sub-sample, by interferometric observations. We used the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument on the Very Large Telescope to obtain extreme adaptive optics assisted scattered light images of the post-AGB binary system AR Puppis. Data have been obtained in the V, I, and H bands. Our observations have produced the first resolved images of AR Puppis’s circumbinary disk and confirm its edge-on orientation. In our high-angular-resolution and high-dynamic-range images we identify several structural components such as a dark mid-plane, the disk surface, and arc-like features. We discuss the nature of these components and use complementary photometric monitoring to relate them to the orbital phase of the binary system. Because the star is completely obscured by the disk at visible wavelengths, we conclude that the long-term photometric variability of the system must be caused by variable scattering, not extinction, of starlight by the disk over the binary orbit. Finally, we discuss how the short disk lifetimes and fast evolution of the host stars compared to the ages at which protoplanetary disks are typically observed make systems like AR Puppis valuable extreme laboratories to study circumstellar disk evolution and constrain the timescale of dust grain growth during the planet formation process.

A water budget dichotomy of rocky protoplanets from 26Al-heating

In contrast to the water-poor inner solar system planets, stochasticity during planetary formation and order of magnitude deviations in exoplanet volatile contents suggest that rocky worlds engulfed in thick volatile ice layers are the dominant family of terrestrial analogues among the extrasolar planet population. However, the distribution of compositionally Earth-like planets remains insufficiently constrained, and it is not clear whether the solar system is a statistical outlier or can be explained by more general planetary formation processes. Here we employ numerical models of planet formation, evolution, and interior structure, to show that a planet's bulk water fraction and radius are anti-correlated with initial $^{26}$Al levels in the planetesimal-based accretion framework. The heat generated by this short-lived radionuclide rapidly dehydrates planetesimals prior to accretion onto larger protoplanets and yields a system-wide correlation of planet bulk abundances, which, for instance, can explain the lack of a clear orbital trend in the water budgets of the TRAPPIST-1 planets. Qualitatively, our models suggest two main scenarios of planetary systems' formation: high-$^{26}$Al systems, like our solar system, form small, water-depleted planets, whereas those devoid of $^{26}$Al predominantly form ocean worlds, where the mean planet radii between both scenarios deviate by up to about 10%.

A giant impact as the likely origin of different twins in the Kepler-107 exoplanet system

Measures of exoplanet bulk densities indicate that small exoplanets with radius less than 3 Earth radii ($R_\oplus$) range from low-density sub-Neptunes containing volatile elements to higher density rocky planets with Earth-like or iron-rich (Mercury-like) compositions. Such astonishing diversity in observed small exoplanet compositions may be the product of different initial conditions of the planet-formation process and/or different evolutionary paths that altered the planetary properties after formation. Planet evolution may be especially affected by either photoevaporative mass loss induced by high stellar X-ray and extreme ultraviolet (XUV) flux or giant impacts. Although there is some evidence for the former, there are no unambiguous findings so far about the occurrence of giant impacts in an exoplanet system. Here, we characterize the two innermost planets of the compact and near-resonant system Kepler-107. We show that they have nearly identical radii (about $1.5-1.6~R_\oplus$), but the outer planet Kepler-107c is more than twice as dense (about $12.6~\rm g\,cm^{-3}$) as the innermost Kepler-107b (about $5.3~\rm g\,cm^{-3}$). In consequence, Kepler-107c must have a larger iron core fraction than Kepler-107b. This imbalance cannot be explained by the stellar XUV irradiation, which would conversely make the more-irradiated and less-massive planet Kepler-107b denser than Kepler-107c. Instead, the dissimilar densities are consistent with a giant impact event on Kepler-107c that would have stripped off part of its silicate mantle. This hypothesis is supported by theoretical predictions from collisional mantle stripping, which match the mass and radius of Kepler-107c.

Ring structure in the MWC 480 disk revealed by ALMA

Gap-like structures in protoplanetary disks are likely related to planet formation processes. In this paper, we present and analyze high-resolution (0.17′′× 0.11′′) 1.3 mm ALMA continuum observations of the protoplanetary disk around the Herbig Ae star MWC 480. Our observations show for the first time a gap centered at ~74 au with a width of ~23 au, surrounded by a bright ring centered at ~98 au from the central star. Detailed radiative transfer modeling of the ALMA image and the broadband spectral energy distribution is used to constrain the surface density profile and structural parameters of the disk. If the width of the gap corresponds to 4–8 times the Hill radius of a single forming planet, then the putative planet would have a mass of 0.4–3 MJ. We test this prediction by performing global three-dimensional smoothed particle hydrodynamic gas/dust simulations of disks hosting a migrating and accreting planet. We find that the dust emission across the disk is consistent with the presence of an embedded planet with a mass of ~2.3 MJ at an orbital radius of ~78 au. Given the surface density of the best-fit radiative transfer model, the amount of depleted mass in the gap is higher than the mass of the putative planet, which satisfies the basic condition for the formation of such a planet.

Excitation and Depletion of the Asteroid Belt in the Early Instability Scenario

Abstract Containing only a few percentages of the mass of the moon, the current asteroid belt is around three to four orders of magnitude smaller than its primordial mass inferred from disk models. Yet dynamical studies have shown that the asteroid belt could not have been depleted by more than about an order of magnitude over the past ∼4 Gyr. The remainder of the mass loss must have taken place during an earlier phase of the solar system’s evolution. An orbital instability in the outer solar system occurring during the process of terrestrial planet formation can reproduce the broad characteristics of the inner solar system. Here, we test the viability of this model within the constraints of the main belt’s low present-day mass and orbital structure. Although previous studies modeled asteroids as massless test particles because of limited computing power, our work uses graphics processing unit acceleration to model a fully self-gravitating asteroid belt. We find that depletion in the main belt is related to the giant planets’ exact evolution within the orbital instability. Simulations that produce the closest matches to the giant planets’ current orbits deplete the main belt by two to three orders of magnitude. These simulated asteroid belts are also good matches to the actual asteroid belt in terms of their radial mixing and broad orbital structure.

Debris Disks among Kepler Solar Rotational Analog Stars

Abstract Observations of circumstellar disks provide a powerful tool for our understanding of planetary system dynamics. Analogs to the solar system asteroid belts, debris disks result from the collision of the remaining solid material of the planet formation process. Even if the presence of disks is now reported for hundreds of stars, its detection around stars similar to the Sun is still very sparse. We report the results of a search for debris disks around Kepler stars with surface physical parameters close to solar values, including rotation period, using observations by the Wide-field infrared Survey Explorer. From the entire sample of Kepler stars, 881 targets were identified with these parameters and only six of them (KIC 1868785, 7267949, 7435796, 10533222, 11352643, and KIC 11666436) show unambiguous infrared excess, for which we determined as debris disk physical parameters. Interestingly, the present study reveals traces of debris disks much more massive and brighter than the solar system zodiacal dust, probably resulting from recent violent collisional events, orbiting stars with ages around the solar values.