Stages Of Planet Formation Research Articles

Measuring the 34S and 33S Isotopic Ratios of Volatile Sulfur during Planet Formation

Stable isotopic ratios constitute powerful tools for unraveling the thermal and irradiation history of volatiles. In particular, we can use our knowledge of the isotopic fractionation processes active during the various stages of star, disk, and planet formation to infer the origins of different volatiles with measured isotopic patterns in our own solar system. Observations of planet-forming disks with the Atacama Large Millimeter/submillimeter Array (ALMA) now readily detect the heavier isotopologues of C, O, and N, while the isotopologue abundances and isotopic fractionation mechanisms of sulfur species are less well understood. Using ALMA observations of the SO and SO2 isotopologues in the nearby, molecule-rich disk around the young star Oph-IRS 48 we present the first constraints on the combined 32S/34S and 32S/33S isotope ratios in a planet-forming disk. Given that these isotopologues likely originate in relatively warm gas (>50 K), like most other Oph-IRS 48 volatiles, SO is depleted in heavy sulfur, while SO2 is enriched compared to solar system values. However, we cannot completely rule out a cooler gas reservoir, which would put the SO sulfur ratios more in line with comets and other solar system bodies. We also constrain the S18O/SO ratio and find the limit to be consistent with solar system values given a temperature of 60 K. Together these observations show that we should not assume solar isotopic values for disk sulfur reservoirs, but additional observations are needed to determine the chemical origin of the abundant SO in this disk, inform on what isotopic fractionation mechanism(s) are at play, and aid in unraveling the history of the sulfur budget during the different stages of planet formation.

Revisiting Collisional Dust Growth in Class 0/I Protostellar Disks: Sweep-up Can Convert a Few 10 M ⊕ of Dust into Kilogram Pebbles in 0.1 Myr

Recent observations suggest that the first stages of planet formation likely take place in the Class 0/I phase of young stellar object evolution, when the star and the disk are still embedded in an infalling envelope. In this study we perform grain coagulation calculations to investigate the very first stage of planet formation, the collisional growth of dust grains, in Class 0/I disks. We find that the slow increase in grain mass by high-velocity collision with much smaller grains (“sweep-up”) allows ∼50 M ⊕ of grains to grow well beyond the fragmentation barrier into ∼kilogram pebbles by the end of Class 0/I (0.1 Myr). We analyze the linear growth and saturation of sweep-up to understand our results quantitatively, and test whether the sweep-up outcome is sensitive to disk parameters and details of the grain coagulation model. The sweep-up pebble population could be important for planet formation, because they are less well-coupled to the gas (compared to the main population below the fragmentation barrier) and therefore more favorable to known mechanisms of dust clump formation (which initiate planetesimal formation). It also contains enough mass to form all planet cores, based on observational estimates of the planet mass budget. Our findings motivate future studies of grain growth and planetesimal formation in Class 0/I disks, including the subsequent evolution of this sweep-up population.

Constraining giant planet formation with synthetic ALMA images of the Solar System’s natal protoplanetary disk

New ALMA observations of protoplanetary disks allow us to probe planet formation in other planetary systems, giving us new constraints on planet formation processes. Meanwhile, studies of our own Solar System rely on constraints derived in a completely different way. However, it is still unclear what features the Solar System protoplanetary disk could have produced during its gas phase. By running 2D isothermal hydro-simulations used as inputs for a dust evolution model, we derive synthetic images at millimeter wavelengths using the radiative transfer code RADMC3D. We find that the embedded multiple giant planets strongly perturb the radial gas velocities of the disk. These velocity perturbations create traffic jams in the dust, producing over-densities different from the ones created by pressure traps and located away from the planets’ positions in the disk. By deriving the images atλ= 1.3 mm from these dust distributions, we show that very high resolution observations are needed to distinguish the most important features expected in the inner part (<15 AU) of the disk. The traffic jams, observable with a high resolution, further blur the link between the number of gaps and rings in disks and the number of embedded planets. We additionally show that a system capable of producing eccentric planets by scattering events that match the eccentricity distributions in observed exoplanets does not automatically produce bright outer rings at large radii in the disk. This means that high resolution observations of disks of various sizes are needed to distinguish between different giant planet formation scenarios during the disk phase, where the giants form either in the outer regions of the disks or in the inner regions. In the second scenario, the disks do not present planet-related features at large radii. Finally, we find that, even when the dust temperature is determined self-consistently, the dust masses derived observationally might be off by up to a factor of ten compared to the dust contained in our simulations due to the creation of optically thick regions. Our study clearly shows that in addition to the constraints from exoplanets and the Solar System, ALMA has the power to constrain different stages of planet formation already during the first few million years, which corresponds to the gas disk phase.

The search for gas in debris discs: ALMA detection of CO gas in HD 36546

ABSTRACT Debris discs represent the last stages of planet formation and as such are expected to be depleted of primordial gas. None the less, in the last few years the presence of cold gas has been reported in ∼20 debris discs from far-infrared to (sub-)mm observations and hot gas has been observed in the optical spectra of debris discs for decades. While the origin of this gas is still uncertain, most pieces of evidence point towards a secondary origin, as a result of collisions and evaporation of small bodies in the disc. In this paper, we present ALMA observations aimed at the detection of CO gas in a sample of eight debris discs with optical gas detections. We report the detection of CO (12CO and 13CO) gas in HD 36546, the brightest and youngest disc in our sample, and provide upper limits to the presence of gas in the remaining seven discs.

Mars’ atmospheric neon suggests volatile-rich primitive mantle

Martian atmospheric neon (Ne) has been detected by Viking and also found as trapped gas in Martian meteorites, though its abundance and isotopic composition have not been well determined. Because the timescale of Ne loss via atmospheric escape estimated from recent measurements with MAVEN is short (0.6–1 × 108 years), the abundance and isotope composition of Martian atmospheric Ne reflect recent atmospheric gas supply mostly from volcanic degassing. Thus, it can serve as a probe for the volatile content of the interior. Here we show that the tentatively-informed atmospheric Ne abundance suggests recent active volcanism and the mantle being richer in Ne than Earth’s mantle today by more than a factor of 5–80. The estimated mantle Ne abundance requires efficient solar nebular gas capture or accretion of Ne-rich materials such as solar-wind-implanted dust in the planet formation stage, both of which provide important constraints on the abundance of other volatile elements in the interior and the accretion history of Mars. More precise determination of atmospheric Ne abundance and isotopic composition by in situ analysis or Mars sample return is crucial for distinguishing the possible origins of Ne.

Growth of aggregates with liquid-like ice shells in protoplanetary discs

ABSTRACT During the first stages of planet formation, the collision growth of dust aggregates in protoplanetary discs (PPDs) is interrupted at the bouncing barrier. Dust aggregates coated by different species of ice turn out to be helpful to shift the bouncing barrier towards larger sizes due to their enhanced sticking properties in some cases. A rarely noticed fact is that H2O ice and H2O–CH3OH–NH3 ice behave liquid-like when UV irradiated within a PPD-similar environment. Dust aggregates coated by these ice species might be damped in collisions due to the liquid-like ice shell, which would effectively result in an increase of the sticking velocity. In this work, collisions of dust aggregates covered by the liquid-like H2O–CH3OH–NH3 ice shell are considered numerically. The coefficient of restitution and the sticking velocity are calculated for different thicknesses of the ice shell. The simulation predicts that an ice-shell thickness of few microns would be sufficient to allow the growth of cm-sized clusters in the PPD.

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.

On the Period-Ratio-Mass-Ratio Correlation of Extra-Solar Multi-Planet Systems

In order to gain possible hints for planet formation from the current data of known extra-solar planets, the period-ratios and mass-ratios of adjacent planet pairs in multi-planet systems are determined. A moderate period-ratio-mass-ratio correlation is found to have a correlation coefficient r=0.5779 with 99% confidence interval (0.464, 0.672). In contrast, for non-adjacent planet pairs, the correlation coefficient is r=0.2820 with 99% confidence interval (0.133,0.419). Our results reveal the imprint of planet-planet interactions of the adjacent planet pairs in a certain fraction of the multi-planet systems during the stage of planet formation.

Far and extreme ultraviolet radiation fields and consequent disc destruction in star-forming regions

ABSTRACT The first stages of planet formation usually occur when the host star is still in a (relatively) dense star-forming region, where the effects of the external environment may be important for understanding the outcome of the planet formation process. In particular, star-forming regions that contain massive stars have strong far-ultraviolet (FUV) and extreme ultraviolet (EUV) radiation fields, which can induce mass-loss from protoplanetary discs due to photoevaporation. In this paper, we present a parameter-space study of the expected FUV and EUV fields in N-body simulations of star-forming regions with a range of initial conditions. We then use recently published models to determine the mass-loss due to photoevaporation from protoplanetary discs. In particular, we focus on the effects of changing the initial degree of spatial structure and initial virial ratio in the star-forming regions, as well as the initial stellar density. We find that the FUV fields in star-forming regions are much higher than in the interstellar medium, even when the regions have stellar densities as low as in the Galactic field, due to the presence of intermediate-mass, and massive, stars (>5 M⊙). These strong radiation fields lead to the destruction of the gas component in protoplanetary discs within 1 Myr, implying that gas giant planets must either form extremely rapidly (<1 Myr), or that they exclusively form in star-forming regions like Taurus, which contain no intermediate-mass or massive stars. The latter scenario is in direct tension with meteoritic evidence from the Solar system that suggests the Sun and its protoplanetary disc was born in close proximity to massive stars.

Erosion of planetesimals by gas flow

The first stages of planet formation take place in protoplanetary disks that are largely made up of gas. Understanding how the gas affects planetesimals in the protoplanetary disk is therefore essential. In this paper, we discuss whether or not gas flow can erode planetesimals. We estimated how much shear stress is exerted onto the planetesimal surface by the gas as a function of disk and planetesimal properties. To determine whether erosion can take place, we compared this with previous measurements of the critical stress that a pebble-pile planetesimal can withstand before erosion begins. If erosion took place, we estimated the erosion time of the affected planetesimals. We also illustrated our estimates with two-dimensional numerical simulations of flows around planetesimals using the lattice Boltzmann method. We find that the wall shear stress can overcome the critical stress of planetesimals in an eccentric orbit within the innermost regions of the disk. The high eccentricities needed to reach erosive stresses could be the result of shepherding by migrating planets. We also find that if a planetesimal erodes, it does so on short timescales. For planetesimals residing outside of 1 au, we find that they are mainly safe from erosion, even in the case of highly eccentric orbits.

Clues to late accretion from Venus’s atmosphere

Whether Earth’s water was delivered early or late in its formation is debated. The composition of Venus’s atmosphere may indicate that late accretion, the final stage of planet formation, delivered little water to the terrestrial planets.

Volatile depletion in planet-forming disks

Newly born stars are surrounded by gas and dust with a flattened axisymmetric distribution termed protoplanetary disk, in which planets are formed. Observations of these objects are necessary for understanding the formation and early evolution of stars and planets, and for revealing the composition of the raw material from which planets are made. Numerical models can extract important parameters from the observational data, including the gas and dust mass of the disk. These parameters are used as input for further modeling, e.g., to calculate the chemical composition of the disk. A consistent thermochemical model should be able to reproduce the abundances of different species in the disk. However, this good wish has been challenged for many disks: models over-predict the emission line intensity of some species; namely, they are depleted (with respect to expectations from canonical models). In this review we show how this disparity indicates that dust evolution has significant effects on gas chemistry, and may indicate the earliest stages of planet formation.

Planets in Binaries: Formation and Dynamical Evolution

Binary systems are very common among field stars, yet the vast majority of known exoplanets have been detected around single stars. While this relatively small number of planets in binaries is probably partly due to strong observational biases, there is, however, statistical evidence that planets are indeed less frequent in binaries with separations smaller than 100 au, strongly suggesting that the presence of a close-in companion star has an adverse effect on planet formation. It is indeed possible for the gravitational pull of the second star to affect all the different stages of planet formation, from proto-planetary disk formation to dust accumulation into planetesimals, to the accretion of these planetesimals into large planetary embryos and, eventually, the final growth of these embryos into planets. For the crucial planetesimal-accretion phase, the complex coupling between dynamical perturbations from the binary and friction due to gas in the proto-planetary disk suggests that planetesimal accretion might be hampered due to increased, accretion-hostile impact velocities. Likewise, the interplay between the binary’s secular perturbations and mean motion resonances lead to unstable regions, where not only planet formation is inhibited, but where a massive body would be ejected from the system on a hyperbolic orbit. The amplitude of these two main effects is different for S- and P-type planets, so that a comparison between the two populations might outline the influence of the companion star on the planet formation process. Unfortunately, at present the two populations (circumstellar or circumbinary) are not known equally well and different biases and uncertainties prevent a quantitative comparison. We also highlight the long-term dynamical evolution of both S and P-type systems and focus on how these different evolutions influence the final architecture of planetary systems in binaries.

Chemical evolution of planetary materials in a dynamic solar nebula

AbstractAs observational facilities improve, providing new insights into the chemistry occurring in protoplanetary disks, it is important to develop more complete pictures of the processes that shapes the chemical evolution of materials during this stage of planet formation. Here we describe how primitive meteorites in our own Solar System can provide insights into the processes that shaped planetary materials early in their evolution around the Sun. In particular, we show how this leads us to expect protoplanetary disks to be very dynamic objects and what modeling and laboratory studies are needed to provide a more complete picture for the early chemical evolution that occurs for planetary systems.

A circumbinary protoplanetary disk in a polar configuration

Nearly all young stars are initially surrounded by `protoplanetary' discs of gas and dust, and in the case of single stars at least 30\% of these discs go on to form planets. The process of protoplanetary disc formation can result in initial misalignments, where the disc orbital plane is different to the stellar equator in single star systems, or to the binary orbital plane in systems with two stars. A quirk of the dynamics means that initially misaligned `circumbinary' discs -- those that surround two stars -- are predicted to evolve to one of two possible stable configurations, one where the disc and binary orbital planes are coplanar, and one where they are perpendicular (a `polar' configuration). Prior work has found coplanar circumbinary discs, but no polar examples were known until now. Here we report the first discovery of a protoplanetary circumbinary disc in the polar configuration, supporting the predictions that such discs should exist. The disc shows some characteristics that are similar to discs around single stars, and that are attributed to dust growth. Thus, the first stages of planet formation appear able to proceed in polar circumbinary discs.

Close-in Super-Earths: The first and the last stages of planet formation in an MRI-accreting disc

We explore in situ formation and subsequent evolution of close-in super-Earths and mini-Neptunes. We adopt a steady-state inner protoplanetary gas disc structure that arises from viscous accretion due to the magneto-rotational instability (MRI). We consider the evolution of dust in the inner disc, including growth, radial drift and fragmentation, and find that dust particles that radially drift into the inner disc fragment severely due to the MRI-induced turbulence. This result has two consequences: (1) radial drift of grains within the inner disc is quenched, leading to an enhancement of dust in the inner regions which scales as dust-to-gas-mass-flux-ratio at ~1 AU; (2) however, despite this enhancement, planetesimal formation is impeded by the small grain size. Nevertheless, assuming that planetary cores are present in the inner disc, we then investigate the accretion of atmospheres onto cores and their subsequent photoevaporation. We then compare our results to the observed exoplanet mass-radius relationship. We find that: (1) the low gas surface densities and high temperatures in the inner disc reduce gas accretion onto cores compared to the minimum mass solar nebula, preventing the cores from growing into hot Jupiters, in agreement with the data; (2) however, our predicted envelope masses are still typically larger than observed ones. Finally, we sketch a qualitative picture of how grains may grow and planetesimals form in the inner disc if grain effects on the ionization levels and the MRI and the back-reaction of the dust on the gas (both neglected in our calculations) are accounted for.

Towards Realistic Understandings of Gas Dynamics in Protoplanetary Disks

AbstractThe gas dynamics of protoplanetary disks (PPDs) plays a crucial role in almost all stages of planet formation, yet it is far from being well understood largely due to the complex interplay among various microphysical processes. Primarily, PPD gas dynamics is likely governed by magnetic fields, and their coupling with the weakly ionized gas is described by non-ideal magnetohydrodynamic (MHD) effects. Incorporating these effects, I will present the first fully global simulations of PPDs that include the most realistic disk microphysics. Accretion and disk evolution is primarily driven by magnetized disk winds with significant mass loss comparable to accretion rate. The overall disk gas dynamics strongly depends on the polarity of large-scale poloidal magnetic field threading the disk owing to the Hall effect. The flow structure in the disk is highly unconventional with major implications on planet formation.

Theoretical calculation of equilibrium Mg isotope fractionation between silicate melt and its vapor

Isotope fractionation during the evaporation of silicate melt and condensation of vapor has been widely used to explain various isotope signals observed in lunar soils, cosmic spherules, calcium–aluminum-rich inclusions, and bulk compositions of planetary materials. During evaporation and condensation, the equilibrium isotope fractionation factor (α) between high-temperature silicate melt and vapor is a fundamental parameter that can constrain the melt’s isotopic compositions. However, equilibrium α is difficult to calibrate experimentally. Here we used Mg as an example and calculated equilibrium Mg isotope fractionation in MgSiO3 and Mg2SiO4 melt–vapor systems based on first-principles molecular dynamics and the high-temperature approximation of the Bigeleisen–Mayer equation. We found that, at 2500 K, δ25Mg values in the MgSiO3 and Mg2SiO4 melts were 0.141 ± 0.004 and 0.143 ± 0.003‰ more positive than in their respective vapors. The corresponding δ26Mg values were 0.270 ± 0.008 and 0.274 ± 0.006‰ more positive than in vapors, respectively. The general $$\alpha - T$$ equations describing the equilibrium Mg α in MgSiO3 and Mg2SiO4 melt–vapor systems were: $$\alpha_{{{\text{Mg}}\left( {\text{l}} \right) - {\text{Mg}}\left( {\text{g}} \right)}} = 1 + \frac{{5.264 \times 10^{5} }}{{T^{2} }}\left( {\frac{1}{m} - \frac{1}{{m^{\prime}}}} \right)$$ and $$\alpha_{{{\text{Mg}}\left( {\text{l}} \right) - {\text{Mg}}\left( {\text{g}} \right)}} = 1 + \frac{{5.340 \times 10^{5} }}{{T^{2} }}\left( {\frac{1}{m} - \frac{1}{{m^{\prime}}}} \right)$$ , respectively, where m is the mass of light isotope 24Mg and m′ is the mass of the heavier isotope, 25Mg or 26Mg. These results offer a necessary parameter for mechanistic understanding of Mg isotope fractionation during evaporation and condensation that commonly occurs during the early stages of planetary formation and evolution.

On the Role of Dissolved Gases in the Atmosphere Retention of Low-mass Low-density Planets

Abstract Low-mass low-density planets discovered by Kepler in the super-Earth mass regime typically have large radii for their inferred masses, implying the presence of H2–He atmospheres. These planets are vulnerable to atmospheric mass loss due to heating by the parent star’s XUV flux. Models coupling atmospheric mass loss with thermal evolution predicted a bimodal distribution of planetary radii, which has gained observational support. However, a key component that has been ignored in previous studies is the dissolution of these gases into the molten core of rock and iron that constitute most of their mass. Such planets have high temperatures (>2000 K) and pressures (∼kbars) at the core-envelope boundary, ensuring a molten surface and a subsurface reservoir of hydrogen that can be 5–10 times larger than the atmosphere. This study bridges this gap by coupling the thermal evolution of the planet and the mass loss of the atmosphere with the thermodynamic equilibrium between the dissolved H2 and the atmospheric H2 (Henry’s law). Dissolution in the interior allows a planet to build a larger hydrogen repository during the planet formation stage. We show that the dissolved hydrogen outgasses to buffer atmospheric mass loss. The slow cooling of the planet also leads to outgassing because solubility decreases with decreasing temperature. Dissolution of hydrogen in the interior therefore increases the atmosphere retention ability of super-Earths. The study highlights the importance of including the temperature- and pressure-dependent solubility of gases in magma oceans and coupling outgassing to planetary evolution models.

Growth and evolution of satellites in a Jovian massive disc

The formation of satellite systems in circum-planetary discs is considered to be similar to the formation of rocky planets in a proto-planetary disc, especially Super-Earths. Thus, it is possible to use systems with large satellites to test formation theories that are also applicable to extrasolar planets. Furthermore, a better understanding of the origin of satellites might yield important information about the environment near the growing planet during the last stages of planet formation. In this work we investigate the formation and migration of the Jovian satellites through N-body simulations. We simulated a massive, static, low viscosity, circum-planetary disc in agreement with the minimum mass sub-nebula model prescriptions for its total mass. In hydrodynamic simulations we found no signs of gaps, therefore type II migration is not expected. Hence, we used analytic prescriptions for type I migration, eccentricity and inclination damping, and performed N-body simulations with damping forces added. Detailed parameter studies showed that the number of final satellites is strong influenced by the initial distribution of embryos, the disc temperature, and the initial gas density profile. For steeper initial density profiles it is possible to form systems with multiple satellites in resonance while a flatter profile favours the formation of satellites close to the region of the Galilean satellites. We show that the formation of massive satellites such as Ganymede and Callisto can be achieved for hotter discs with an aspect ratio of H/r ~ 0.15 for which the ice line was located around 30 R_J.