skip to main content

Title: Maximum temperatures in evolving protoplanetary discs and composition of planetary building blocks
ABSTRACT The maximum temperature and radial temperature profile in a protoplanetary disc are important for the condensation of different elements in the disc. We simulate the evolution of a set of protoplanetary discs from the collapse of their progenitor molecular cloud cores as well as the dust decoupling within the discs as they evolve. We show how the initial properties of the cloud cores affect the thermal history of the protoplanetary discs using a simple viscous disc model. Our results show that the maximum mid-plane temperature in the disc occurs within 0.5 au. It increases with the initial cloud temperature and decreases with its angular velocity and the viscosity of the disc. From the observed properties of the molecular cloud cores, we find the median value of the maximum temperature is around 1250 K, with roughly 90 per cent of them being less than 1500 K – a value that is lower than the 50 per cent condensation temperatures of most refractory elements. Therefore, only cloud cores with high initial temperatures or low-angular velocities and/or low viscosities within the planet-forming discs will result in refractory-rich planetesimals. To reproduce the volatile depletion pattern of CM, CO, and CV chondrites and the terrestrial planets in Solar system, one more » must either have rare properties of the initial molecular cloud cores like high core temperature, or other sources of energy to heat the disc to sufficiently high temperatures. Alternatively, the volatile depletion observed in these chondrites may be inherited from the progenitor molecular cloud. « less
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
5254 to 5262
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Partial condensation of dust from the Solar nebula is likely responsible for the diverse chemical compositions of chondrites and rocky planets/planetesimals in the inner Solar system. We present a forward physical–chemical model of a protoplanetary disc to predict the chemical compositions of planetary building blocks that may form from such a disc. Our model includes the physical evolution of the disc and the condensation, partial advection, and decoupling of the dust within it. The chemical composition of the condensate changes with time and radius. We compare the results of two dust condensation models: one where an element condenses when the mid-plane temperature in the disc is lower than the 50 per cent condensation temperature ($\rm T_{50}$) of that element and the other where the condensation of the dust is calculated by a Gibbs free energy minimization technique assuming chemical equilibrium at local disc temperature and pressure. The results of two models are generally consistent with some systematic differences of ∼10 per cent depending upon the radial distance and an element’s condensation temperature. Both models predict compositions similar to CM, CO, and CV chondrites provided that the decoupling time-scale of the dust is of the order of the evolution time-scale of the disc ormore »longer. If the decoupling time-scale is too short, the composition deviates significantly from the measured values. These models may contribute to our understanding of the chemical compositions of chondrites, and ultimately the terrestrial planets in the Solar system, and may constrain the potential chemical compositions of rocky exoplanets.« less
  2. Comparing compositional models of the terrestrial planets provides insights into physicochemical processes that produced planet-scale similarities and differences. The widely accepted compositional model for Mars assumes Mn and more refractory elements are in CI chondrite proportions in the planet, including Fe, Mg, and Si, which along with O make up >90% of the mass of Mars. However, recent improvements in our understandings on the composition of the solar photosphere and meteorites challenge the use of CI chondrite as an analog of Mars. Here we present an alternative model composition for Mars that avoids such an assumption and is based on data from Martian meteorites and spacecraft observations. Our modeling method was previously applied to predict the Earth’s composition. The model establishes the absolute abundances of refractory lithophile elements in the bulk silicate Mars (BSM) at 2.26 times higher than that in CI carbonaceous chondrites. Relative to this chondritic composition, Mars has a systematic depletion in moderately volatile lithophile elements as a function of their condensation temperatures. Given this finding, we constrain the abundances of siderophile and chalcophile elements in the bulkMars and its core. The Martian volatility trend is consistent with <7 wt% S in its core, which is significantlymore »lower than that assumed in most core models (i.e., >10 wt% S). Furthermore, the occurrence of ringwoodite at the Martian core-mantle boundary might have contributed to the partitioning of O and H into the Martian core.« less
  3. Abstract Prestellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures (<10 K) allow the formation of thick icy dust mantles, which will be partially preserved in future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particular CO, are heavily depleted in prestellar cores due to the efficient molecular freeze-out onto the surface of cold dust grains. However, N-bearing species such as NH 3 and, in particular, its deuterated isotopologues appear to maintain high abundances where CO molecules are mainly in the solid phase. Thanks to ALMA, we present here the first clear observational evidence of NH 2 D freeze-out toward the L1544 prestellar core, suggestive of the presence of a “complete depletion zone” within a ≃1800 au radius, in agreement with astrochemical prestellar core model predictions. Our state-of-the-art chemical model coupled with a non-LTE radiative transfer code demonstrates that NH 2 D becomes mainly incorporated in icy mantles in the central 2000 au and starts freezing out already at ≃7000 au. Radiative transfer effects within the prestellar core cause the NH 2 D(1 11 − 1 01 ) emission tomore »appear centrally concentrated, with a flattened distribution within the central ≃3000 au, unlike the 1.3 mm dust continuum emission, which shows a clear peak within the central ≃1800 au. This prevented NH 2 D freeze-out from being detected in previous observations, where the central 1000 au cannot be spatially resolved.« less
  4. ABSTRACT Discs are essential to the formation of both stars and planets, but how they form in magnetized molecular cloud cores remains debated. This work focuses on how the disc formation is affected by turbulence and ambipolar diffusion (AD), both separately and in combination, with an emphasis on the protostellar mass accretion phase of star formation. We find that a relatively strong, sonic turbulence on the core scale strongly warps but does not completely disrupt the well-known magnetically induced flattened pseudo-disc that dominates the inner protostellar accretion flow in the laminar case, in agreement with previous work. The turbulence enables the formation of a relatively large disc at early times with or without AD, but such a disc remains strongly magnetized and does not persist to the end of our simulation unless a relatively strong AD is also present. The AD-enabled discs in laminar simulations tend to fragment gravitationally. The disc fragmentation is suppressed by initial turbulence. The AD facilitates the disc formation and survival by reducing the field strength in the circumstellar region through magnetic flux redistribution and by making the field lines there less pinched azimuthally, especially at late times. We conclude that turbulence and AD complement eachmore »other in promoting disc formation. The discs formed in our simulations inherit a rather strong magnetic field from its parental core, with a typical plasma-β of order a few tens or smaller, which is 2–3 orders of magnitude lower than the values commonly adopted in magnetohydrodynamic simulations of protoplanetary discs. To resolve this potential tension, longer term simulations of disc formation and evolution with increasingly more realistic physics are needed.« less

    Despite many methods developed to find young massive planets in protoplanetary discs, it is challenging to directly detect low-mass planets that are embedded in discs. On the other hand, the core-accretion theory suggests that there could be a large population of embedded low-mass young planets at the Kelvin-Helmholtz (KH) contraction phase. We adopt both 1D models and 3D simulations to calculate the envelopes around low-mass cores (several to tens of M⊕) with different luminosities, and derive their thermal fluxes at radio wavelengths. We find that, when the background disc is optically thin at radio wavelengths, radio observations can see through the disc and probe the denser envelope within the planet’s Hill sphere. When the optically thin disc is observed with the resolution reaching one disc scale height, the radio thermal flux from the planetary envelope around a 10 M⊕ core is more than 10 per cent higher than the flux from the background disc. The emitting region can be extended and elongated. Finally, our model suggests that the au-scale clump at 52 au in the TW Hydrae disc revealed by ALMA is consistent with the envelope of an embedded 10–20 M⊕ planet, which can explain the detected flux, the spectral index dip,more »and the tentative spirals. The observation is also consistent with the planet undergoing pebble accretion. Future ALMA and ngVLA observations may directly reveal more such low-mass planets, enabling us to study core growth and even reconstruct the planet formation history using the embedded ‘protoplanet’ population.

    « less