Evolution Of Protoplanetary Disks Research Articles

On the Formation of Planets in the Milky Way’s Thick Disk

Abstract Exoplanet demographic surveys have revealed that close-in (≲1 au) small planets orbiting stars in the Milky Way’s thick disk are ∼50% less abundant than those orbiting stars in the Galactic thin disk. One key difference between the two stellar populations is the time at which they emerged: thick-disk stars are the likely product of cosmic noon (redshift z ∼ 2), an era characterized by high star formation rate, massive and dense molecular clouds, and strong supersonic turbulence. Solving for the background radiation field in these early star-forming regions, we demonstrate that protoplanetary disks at cosmic noon experienced radiation fields up to ∼7 orders of magnitude more intense than in solar neighborhood conditions. Coupling the radiation field to a one-dimensional protoplanetary disk evolution model, we find that external UV photoevaporation destroys protoplanetary disks in just ∼0.2–0.5 Myr, limiting the timescale over which planets can assemble. Disk temperatures exceed the sublimation temperatures of common volatile species for ≳Myr timescales, predicting more spatial homogeneity in gas chemical composition. Our calculations imply that the deficit in planet occurrence around thick-disk stars should be even more pronounced for giant planets, particularly those at wide orbital separations, predicting a higher rocky-to-giant planet ratio in the Galactic thick disk versus thin disk.

There is no disk mass budget problem of planet formation

The inferred dust masses from Class II protoplanetary disk observations are lower than or equal to the masses of the observed exoplanet systems. This poses the question of how planets form if their natal environments do not contain enough mass. This hypothesis has entered the literature as the ``mass budget problem'' of planet formation. We utilized numerical simulations of planet formation via pebble and gas accretion, including migration, in a viscously evolving protoplanetary disk, while tracing the time evolution of the dust mass. As expected, we found that the presence of a giant planet in the disk can influence the evolution of the disk itself and prevent rapid dust mass loss by trapping the dust outside its orbit. Early formation is crucial for giant planet formation, as we found in our previous work; therefore, our findings strengthen the hypothesis that planet formation has already occurred or is ongoing in Class II disks. More importantly, we find that the optically thin dust mass significantly underestimates the total dust mass in the presence of a dust-trapping deep gap. We also show that the beam convolution smears out the feature from a deep gap, especially if the planet forms in the inner disk. Such hidden dust mass, along with early planet formation, could be the answer to the hypothetical mass budget problem.

A tell-tale tracer for externally irradiated protoplanetary disks: Comparing the [C I] 8727 Å line and ALMA observations in proplyds

The evolution of protoplanetary disks in regions with massive OB stars is influenced by externally driven winds that deplete the outer parts of these disks. The winds have previously been studied via forbidden oxygen emission lines, which also arise in isolated disks in low-mass star-forming regions (SFRs) with weak external UV fields in photoevaporative or magnetic (internal) disk winds. It is crucial to determine how to disentangle external winds from internal ones. Here, we report a proxy for unambiguously identifying externally driven winds with a forbidden line of neutral atomic carbon, [C I] 8727 Å. We compare for the first time the spatial location of the emission in the [O I] 5577 Å, [O I] 6300 Å, and [C I] 8727 Å lines traced by VLT/MUSE-NFM with the ALMA Band 7 continuum disk emission in a sample of 12 proplyds in the Orion Nebula Cluster (ONC). We confirm that the [O I] 5577 Å emission is co-spatial with the disk emission, whereas that of [O I] 6300 Å is emitted both on the disk surface and on the ionization front of the proplyds. We show for the first time that the [C I] 8727 Å line is also co-spatial with the disk surface in proplyds, as seen in the MUSE and ALMA data comparison. The peak emission is compatible with the stellar location in all cases, apart from one target with high relative inclination with respect to the ionizing radiation, where the peak emission is located at the disk edge in the direction of the ionizing radiation. To verify whether the [C I] 8727 Å line is detected in regions where external photoevaporation is not expected, we examined VLT/X-Shooter spectra for young stars in low-mass SFRs. Although the [O I] 5577 Å and 6300 Å lines are well detected in all these targets, the total detection rate is ≪10% in the case of the [C I] 8727 Å line. This number increases substantially to a ∼40% detection rate in σ-Orionis, a region with higher UV radiation than in low-mass SFRs, but lower than in the ONC. The spatial location of the [C I] 8727 Å line emission and the lack of its detection in isolated disks in low-mass SFRs strongly suggest that this line is a tell-tale tracer of externally driven photoevaporative winds, which agrees with recent excitation models.

Changing disc compositions via internal photoevaporation

The chemical evolution of protoplanetary discs is a complex process that is not fully understood. Several factors influence the final spatial distribution of atoms and molecules in the disc. One such factor is the inward drift and evaporation of volatile-rich pebbles that can enrich the inner disc with vapour. In particular, the inner disc is first enriched with evaporating water-ice, resulting in a low C/O ratio, before carbon-rich gas from the outer disc – originating from the evaporation of CO, CO2, and CH4 ice – is transported viscously inwards, elevating the C/O ratio again. However, it is unclear how internal photoevaporation – which carries away gas and opens gaps in the disc that can block inward drifting pebbles – affects the chemical composition of the disc. Our goal is to study how and to what extent internal photoevaporation and the subsequent opening of gaps influence the chemical evolution of protoplanetary discs around solar-like stars (M* = 1 M⊙), where we specifically focus on the C/O ratio and the water content. To carry out our simulations, we use a semi-analytical 1D disc model. The code chemcomp includes viscous evolution and heating, pebble growth and drift, pebble evaporation and condensation, as well as a simple chemical partitioning model for the disc. We show that internal photoevaporation plays a major role in the evolution of protoplanetary discs and their chemical composition: As photoevaporation opens a gap, inward drifting pebbles are stopped and can no longer contribute to the volatile content in the gas. In addition, volatile-rich gas from the outer disc, originating from evaporated CO, CO2, or CH4 ice, is carried away by the photoevaporative winds. Consequently, the C/O ratio in the inner disc remains low. In contrast, gaps opened by giant planets still allow the gas to pass, resulting in an elevated C/O ratio in the inner disc, similar to the evolution of viscous discs without internal photoevaporation. This opens the possibility to distinguish observationally between these two scenarios when measuring the C/O ratio, implying that we can infer the root cause of deep gap structures when observing protoplanetary discs. In the case of a clear separation of the disc by photoevaporation, we additionally find an elevated water content in the inner disc, because the water vapour and ice undergo a cycle of evaporation and recondensation, preventing the inward accretion of water onto the star, in contrast to the situation for hydrogen and helium. We conclude that it is very difficult to achieve supersolar C/O ratios in the inner parts of protoplanetary discs when taking internal photoevaporation into account. This indicates the potential importance of photoevaporation for understanding the chemical evolution of these discs and the planets forming in them.

TriPoD: Tri-Population size distributions for Dust evolution

Context. Dust coagulation and fragmentation impact the structure and evolution of protoplanetary disks and set the initial conditions for planet formation. Dust grains dominate the opacities, they determine the cooling times of the gas via thermal accommodation in collisions, they influence the ionization state of the gas, and the available grain surface area is an important parameter for the chemistry in protoplanetary disks. Therefore, dust evolution is an effect that should not be ignored in numerical studies of protoplanetary disks. Available dust coagulation models are, however, too computationally expensive to be implemented in large-scale hydrodynamic simulations. This limits detailed numerical studies of protoplanetary disks, including these effects, mostly to one-dimensional models. Aims. We aim to develop a simple – yet accurate – dust coagulation model that can be easily implemented in hydrodynamic simulations of protoplanetary disks. Our model shall not significantly increase the computational cost of simulations and provide information about the local grain size distribution. Methods. The local dust size distributions are assumed to be truncated power laws. Such distributions can be fully characterized by only two dust fluids (large and small grains) and a maximum particle size, truncating the power law. We compare our model to state- of-the-art dust coagulation simulations and calibrate it to achieve a good fit with these sophisticated numerical methods. Results. Running various parameter studies, we achieved a good fit between our simplified three-parameter model and DustPy, a state-of-the-art dust coagulation software. Conclusions. We present TriPoD, a sub-grid dust coagulation model for the PLUTO code. With TriPoD, we can perform twodimensional, vertically integrated dust coagulation simulations on top of a hydrodynamic simulation. Studying the dust distributions in two-dimensional vortices and planet-disk systems is thus made possible.

Formation and evolution of a protoplanetary disk: Combining observations, simulations, and cosmochemical constraints

The formation and evolution of protoplanetary disks remains elusive. We have numerous astronomical observations of young stellar objects of different ages with their envelopes and/or disks. Moreover in the last decade there has been tremendous progress in numerical simulations of star and disk formation. New simulations use realistic equations of state for the gas and treat the interaction of matter and the magnetic field with the full set of nonideal magnetohydrodynamic (MHD) equations. However, it is still not fully clear how a disk forms and whether it happens from inside-out or outside-in. Open questions remain regarding where material is accreted onto the disk and comes from, how dust evolves in disks, and the timescales of appearance of disk's structures. These unknowns limit our understanding of how planetesimals and planets form and evolve. We attempted to reconstruct the evolutionary history of the protosolar disk, guided by the large amount of cosmochemical constraints derived from the study of meteorites, while using astronomical observations and numerical simulations as a guide to pinpointing plausible scenarios. Our approach is highly interdisciplinary and we do not present new observations or simulations in this work. Instead, we combine, in an original manner, a large number of published results concerning young stellar objects observations, and numerical simulations, along with the chemical, isotopic and petrological nature of meteorites. We have achieved a plausible and coherent view of the evolution of the protosolar disk that is consistent with cosmochemical constraints and compatible with observations of other protoplanetary disks and sophisticated numerical simulations. The evidence that high-temperature condensates, namely, calcium-aluminum inclusions (CAIs) and amoeboid olivine aggregates (AOAs), formed near the protosun before being transported to the outer disk can be explained in two ways: there could have either been an early phase of vigorous radial spreading of the disk that occurred or fast transport of these condensates from the vicinity of the protosun toward large disk radii via the protostellar outflow. The assumption that the material accreted toward the end of the infall phase was isotopically distinct allows us to explain the observed dichotomy in nucleosynthetic isotopic anomalies of meteorites. It leads us toward intriguing predictions on the possible isotopic composition of refractory elements in comets. At a later time, when the infall of material waned, the disk started to evolve as an accretion disk. Initially, dust drifted inward, shrinking the radius of the dust component to $ 45$ au, probably about to about half of the width of the gas component. Next, structures must have emerged, producing a series of pressure maxima in the disk, which trapped the dust on Myr timescales. This allowed planetesimals to form at radically distinct times without significantly changing any of the isotopic properties. We also conclude that there was no late accretion of material onto the disk via streamers. The disk disappeared at about 5 My, as indicated by paleomagnetic data in meteorites. The evolution of the protosolar disk seems to have been quite typical in terms of size, lifetime, and dust behavior. This suggests that the peculiarities of the Solar System with respect to extrasolar planetary systems probably originate from the chaotic nature of planet formation and not from the properties of the parental disk itself.

JWST study of the DG Tau B disk-wind candidate

Context. The origin of outflows and their impact on protoplanetary disk evolution and planet-formation processes are still crucial open questions. DG Tau B is a Class I protostar associated with a structured disk and a rotating conical CO outflow and was recently identified by ALMA as one of the best CO disk wind candidate. It is therefore a perfect target for studying these questions. Aims. We aim to map and study any outflow component intermediate between the axial jet and the CO outflow in order to constrain the origin (irradiated or shocked disk wind or swept-up material) of the redshifted molecular outflows in DG Tau B. Methods. We analyzed observations obtained with James Webb Space Telescope NIRSpec-IFU and NIRCam, supplemented with IFU data from the SINFONI/VLT instrument. We investigated the morphology, kinematics, and excitation conditions of the ro-vibrational H2 emission lines and their relation with the atomic jet and CO outflow. We focus our analysis on the redshifted outflow lobe. Results. We observe a global layered structure of the redshifted outflows in DG Tau B, with the atomic jet inside the H2 cavity, which in turn is nested inside the CO conical outflow. We also find temperature, velocity, and collimation to be increasing towards the axis of the flow. The redshifted H2 emission traces a narrow conical cavity (semi-opening angle of 9.4°) that is wider than the axial jet but nested just inside the CO outflow. Both the jet and the H2 cavity originate from the innermost regions of the disk (r0 < 6 au). The redshifted H2 cavity flows with a constant vertical velocity of Vz = 22.5 ± 0.8 km s−1, twice faster than the conical CO flow. The excitation conditions imply a hot H2 gas (Tex ≃ 2200 K) with an average mass flux of Ṁ(H2) = 3 × 10−11 M⊙ yr−1, which is significantly lower than the jet and CO values. Conclusions. The global layered H2-CO structure in temperature, velocity, and collimation in the DG Tau B redshifted lobe is consistent with a magneto-hydrodynamic disk wind scenario. The hot H2 could trace the inner, dense photodissociation layer in the wind. An H2 launching region at disk radii of 0.2−0.4 au combined with a large ejection efficiency (ξ ≃ 1) would account for the mass flux and kinematics. Alternatively, the near-IR ro-vibrational H2 could be emitted in the interaction layer driven by successive jet bow shocks into an outer disk wind or envelope. Further constraints on both scenarios will be obtained from the analysis of MIRI observations.

A question of personalities: evolution of viscous and wind-driven protoplanetary discs in the presence of dead zones

A question of personalities: evolution of viscous and wind-driven protoplanetary discs in the presence of dead zones

Does Metallicity Affect the Protoplanetary Disk Fraction? Answers from the Outer Milky Way

The role of metallicity in shaping protoplanetary disk evolution remains poorly comprehended. This study analyzes the disk fraction of 10 young (0.9–2.1 Myr) and low-metallicity (0.34–0.83 Z ⊙) clusters located in the outer Milky Way with Galactocentric distances between 10 and 13 kpc. Using JHK data obtained from UKIDSS, the calculated disk fraction values for low-mass stars (0.2–2 M ⊙) ranged from 42% to 7%. To enhance the statistical reliability of our analysis, eight additional low-metallicity clusters are sourced from previous studies with metallicity range 0.25–0.85 Z ⊙ along with our sample, resulting in a total of 18 regions with low metallicity. We find that low-metallicity clusters exhibit on average a 2.6 ± 0.2 times lower disk fraction compared to solar-metallicity clusters in all the age bins we have. Within the age range we can probe, our study does not find evidence of faster disk decay in subsolar-metallicity regions compared to solar-metallicity regions. Furthermore, we observe a positive correlation between cluster disk fraction and metallicity for two different age groups of 0.3–1.4 and 1.4–2.5 Myr. We emphasize that both cluster age and metallicity significantly affect the fraction of stars with evidence of inner disks.

Can near-to-mid Infrared Spectral Energy Distribution Quantitatively Trace Protoplanetary Disk Evolution?

Infrared (IR) spectral energy distribution (SED) is the major tracer of protoplanetary disks. It was recently proposed to use the near-to-mid IR (or K-24) SED slope α defined between 2 and 24μm as a potential quantitative tracer of disk age. We critically examine the viability of this idea and confront it with additional statistics of IR luminosities and SED shapes. We point out that, because the statistical properties of most of the complicated physical factors involved in disk evolution are still poorly understood in a quantitative sense, the only viable way is to assume them to be random so that an idealized “average disk” can be defined, which allows the α histogram to trace its age. We confirm that the statistics of the zeroth order (luminosity), first order (slope α), and second order characteristics (concavity) of the observed K-24 SEDs indeed carry useful information upon the evolutionary processes of the “average disk”. We also stress that intrinsic diversities in K-24 SED shapes and luminosities are always large at the level of individual stars so that the application of the evolutionary path of the “average disk” to individual stars must be done with care. The data of most curves in plots are provided on GitHub (Disk-age package https://github.com/starage/disk-age/).

Co-evolution of dust grains and protoplanetary disks. II. Structure and evolution of protoplanetary disks: An analytical approach

Abstract In our previous study (Tsukamoto et al. 2023b, PASJ, 75, 835), we investigated the formation and early evolution of protoplanetary disks with 3D non-ideal magnetohydrodynamics simulations considering dust growth, and found that the modified equations of the conventional steady accretion disk model that consider magnetic braking, dust growth, and ambipolar diffusion reproduce the disk structure (such as density and vertical magnetic field) obtained from simulations very well. In this paper, as a sequel to our previous study, we analytically investigate the structure and evolution of protoplanetary disks corresponding to Class 0/I young stellar objects using the modified steady accretion disk model combining an analytical model of envelope accretion. We estimate that the disk radius is several astronomical units at the disk formation epoch and increases to several hundred astronomical units at the end of the accretion phase. The disk mass is estimated to be $0.01 \lesssim M_{\rm disk} \lesssim 0.1 \, M_\odot$ for a disk with a radius of several tens of astronomical units and a mass accretion rate of $\dot{M}_{\rm disk} \sim 10^{-6} \, M_\odot \,\, {\rm yr^{-1}}$. These estimates seems to be consistent with recent observations. We also found that, with typical disk ionization rates (ζ ≳ 10−19 s−1) and a moderate mass accretion rate ($\dot{M}_{\rm disk}\gtrsim 10^{-8} \, M_\odot \,\, {\rm yr^{-1}}$), magnetorotational instability is suppressed in the disk because of low plasma β and efficient ambipolar diffusion. We argue that the radial profile of specific angular momentum (or rotational velocity) at the disk outer edge should be continuously connected to that of the envelope if the disk evolves by magnetic braking, and should be discontinuous if the disk evolves by an internal angular momentum transport process such as gravitational instability or magnetorotational instability. Future detailed observations of the specific angular momentum profile around the disk outer edge are important for understanding the angular momentum transport mechanism of protoplanetary disks.

The interplay between forming planets and photoevaporating discs

Context. Disc winds and planet–disc interactions are two crucial mechanisms that define the structure, evolution, and dispersal of protoplanetary discs. While winds are capable of removing material from discs, eventually leading to their dispersal, massive planets can shape their disc by creating sub-structures such as gaps and spiral arms. Aims. We studied the interplay between an X-ray photoevaporative disc wind and the sub-structures generated due to planet–disc interactions to determine how their mutual interactions affect the disc’s and the planet’s evolution. Methods. We performed 3D hydrodynamic simulations of viscous discs (α = 6.9 × 10−4) that host a Jupiter-like planet and undergo X-ray photoevaporation. We traced the gas flows within the disc and wind and measured the rate of accretion onto the planet, as well as the gravitational torque that is acting on it. Results. Our results show that the planetary gap removes the wind’s pressure support, allowing wind material to fall back into the gap. This opens new pathways for material from the inner disc (and part of the outer disc) to be redistributed through the wind towards the gap. Consequently, the gap becomes shallower and the flow of mass across the gap in both directions is significantly increased, as is the planet’s mass-accretion rate (by factors of ≈5 and ≈2, respectively). Moreover, the wind-driven redistribution results in a denser inner disc and a less dense outer disc, which, combined with the recycling of a significant portion of the inner wind, leads to longer lifetimes for the inner disc, contrary to the expectation in a planet-induced photoevaporation scenario that has been proposed in the past.

Evolution of the relation between the mass accretion rate and the stellar and disk mass from brown dwarfs to stars

The time evolution of the dependence of the mass accretion rate with the stellar mass and the disk mass represents a fundamental way to understand the evolution of protoplanetary disks and the formation of planets. In this work, we present observations with X-shooter of 26 Class II very low-mass stars (< 0.2 M⊙) and brown dwarfs in the Ophiuchus, Chamaeleon-I, and Upper Scorpius star-forming regions. These new observations extend the measurement of the mass accretion rate down to spectral type (SpT) M9 (∼0.02 M⊙) in Ophiuchus and Chamaeleon-I and add 11 very-low-mass stars to the sample of objects previously studied with broadband spectroscopy in Upper Scorpius. We obtained the spectral type and extinction, as well as the physical parameters of the sources. We used the intensity of various emission lines in the spectra of these sources to derive the accretion luminosity and mass accretion rates for the entire sample. Combining these new observations with data from the literature, we compare relations between accretion and stellar and disk properties of four different star-forming regions with different ages: Ophiuchus (∼1 Myr), Lupus (∼2 Myr), Chamaeleon-I (∼3 Myr), and Upper Scorpius (5−12 Myr). We find the slopes of the accretion relationships (L* − Lacc, M∗ − Ṁacc) to steepen in the 1−3 Myr age range (i.e., between Ophiuchus, Lupus, and Chamaeleon-I) and that both relationships may be better described with a single power law. We find that previous claims for a double power-law behavior of the M∗ − Ṁacc relationship may have been triggered by the use of a different SpT–Teff scale. We also find the relationship between the protoplanetary disk mass and the mass accretion rate of the stellar population to steepen with time down to the age of Upper Scorpius. Overall, we observe hints of a faster evolution into low accretion rates of low-mass stars and brown dwarfs. At the same time, we also find that brown dwarfs present higher Mdisk/Ṁacc ratios (i.e., longer accretion depletion timescales) than stars in Ophiuchus, Lupus, and Cha-I. This apparently contradictory result may imply that the evolution of protoplanetary disks around brown dwarfs may be different than what is seen in the stellar regime.

Thermal processing of primordial pebbles in evolving protoplanetary disks

During protoplanetary disk formation, dust grains located in the outer disk retain their pristine icy composition, while solids in the inner stellar-heated disk undergo volatile loss. This process may have left a fossil record in Solar System material, showing different nucleosynthetic imprints that have been attributed to different degrees of thermal processing. However, it remains unclear how a large mass fraction of thermally processed inner-disk pebbles is produced and how these grains are subsequently transported throughout the disk. In this work, we numerically investigate the evolution in time of a two-component pebble disk consisting of both pristine pebbles and those that underwent ice sublimation. We find that stellar outbursts exceeding 1000 times the solar luminosity are efficient in thermally altering, through ice sublimation, a large mass fraction of pebbles (around 80%). After the establishment of this initial radial dust composition gradient throughout the disk, the subsequent mixing and inward drift of pristine outer-disk pebbles alter the inner disk bulk composition from processed to more unprocessed in time. Therefore, if processed pebbles without ice mantles have an isotopic composition similar to ureilite meteorites from the inner Solar System, inner-disk minor bodies forming from the early pebble flux (<1 Myr) will be isotopically ureilite-like, while later-formed bodies will be increasingly admixed with the signature of the lateincoming, CI chondrite-like unprocessed pebbles. This appears to be largely consistent with the trend seen between the accretion age of different meteoric classes and their different stable isotope composition anomalies (in μ54Cr, μ48Ca, μ30Si, and μ58Ni), but further work may be needed to explain the role of isotopically anomalous refractory inclusions and anomaly trends in other elements. Our findings further support an early thermal processing of ice mantles via stellar outbursts that are common around young Sun-like stars.

The Dynamic, Chimeric Inner Disk of PDS 70

Transition disks, with inner regions depleted in dust and gas, could represent later stages of protoplanetary disk evolution when newly formed planets are emerging. The PDS 70 system has attracted particular interest because of the presence of two giant planets in orbits at tens of astronomical units within the inner disk cavity, at least one of which is itself accreting. However, the region around PDS 70 most relevant to understanding the planet populations revealed by exoplanet surveys of middle-aged stars is the inner disk, which is the dominant source of the system’s excess infrared emission but only marginally resolved by the Atacama Large Millimeter/submillimeter Array. Here we present and analyze time-series optical and infrared photometry and spectroscopy that reveal the inner disk to be dynamic on timescales of days to years, with occultation by submicron dust dimming the star at optical wavelengths, and 3–5 μm emission varying due to changes in disk structure. Remarkably, the infrared emission from the innermost region (nearly) disappears for ∼1 yr. We model the spectral energy distribution of the system and its time variation with a flattened warm (T ≲ 600 K) disk and a hotter (1200 K) dust that could represent an inner rim or wall. The high dust-to-gas ratio of the inner disk, relative to material accreting from the outer disk, means that the former could be a chimera consisting of depleted disk gas that is subsequently enriched with dust and volatiles produced by collisions and evaporation of planetesimals in the inner zone.

Isotopic evolution of the inner solar system revealed by size-dependent oxygen isotopic variations in chondrules

The systematic isotopic difference between refractory inclusions and chondrules, particularly for oxygen, has long indicated an isotopic evolution of the solar protoplanetary disk. However, it remains underconstrained whether such evolution continued during chondrule formation. Intrigued by past reports of the size-dependent oxygen isotopic compositions of chondrules in ordinary chondrites (OC), we analyzed type I olivine-rich chondrules of various sizes in two LL3 chondrites. Although most chondrules show positive Δ17O values comparable to bulk ordinary chondrites, a population of smaller (less than about 0.1 mm2 in cross-section), including many isolated olivine grains (sensu lato), are 16O-enriched (with Δ17O values down to −13.2 ‰). Literature data allow the same observation for R chondrites. All sub-TFL type I chondrules (i.e., Δ17O < 0) chondrules have Mg# > 97 mol% while the supra-TFL type I chondrule olivines extend to the formal boundary with type II chondrules (i.e., Mg# = 90 mol%). The sub-TFL chondrules are likely genetically related to isotopically similar aluminum-rich chondrules described in the literature. They therefore must have formed earlier than most OC and R chondrules when the inner disk was still sub-TFL. This interpretation is supported by the presence of similarly 16O-rich relict grains in supra-TFL OC and R chondrules that must be remains of this incompletely recycled precursor material. The non-carbonaceous reservoir was thus still evolving isotopically towards 16O-poorer composition when chondrule formation began, whether by mixing with outer disk material, late accretion streamers and/or an increase in the solid/gas ratio due to magnetothermal disk winds.

The Effect of Ultraviolet Photon Pumping of H2 in Dust-deficient Protoplanetary Disks

We perform radiation hydrodynamics simulations to study the structure and evolution of a photoevaporating protoplanetary disk. Ultraviolet and X-ray radiation from the host star heats the disk surface, where H2 pumping also operates efficiently. We run a set of simulations in which we varied the number of dust grains or the dust-to-gas mass ratio, which determines the relative importance between photoelectric heating and H2 pumping. We show that H2 pumping and X-ray heating contribute more strongly to the mass loss of the disk when the dust-to-gas mass ratio is D≤10−3 . The disk mass-loss rate decreases with a lower dust amount, but remains around 10−10−11 M ⊙yr−1. In these dust-deficient disks, H2 pumping enhances photoevaporation from the inner disk region and shapes the disk mass-loss profile. We thus argue that the late-stage disk evolution is affected by the ultraviolet H2 pumping effect. The mass-loss rates derived from our simulations can be used in the study of long-term disk evolution.

Dust growth and pebble formation in the initial stages of protoplanetary disk evolution

Aims. The initial stages of planet formation may start concurrently with the formation of a gas-dust protoplanetary disk. This makes the study of the earliest stages of protoplanetary disk formation crucially important. Here we focus on dust growth and pebble formation in a protoplanetary disk that is still accreting from a parental cloud core. Methods. We have developed an original three-dimensional numerical hydrodynamics code, which computes the collapse of rotating clouds and disk formation on nested meshes using a novel hybrid Coarray Fortran-OpenMP approach for distributed and shared memory parallelization. Dust dynamics and growth are also included in the simulations. Results. We found that the dust growth from ~1 µm to 1–10 mm already occurs in the initial few thousand years of disk evolution but the Stokes number hardly exceeds 0.1 because of higher disk densities and temperatures compared to the minimum mass Solar nebular. The ratio of the dust-to-gas vertical scale heights remains rather modest, 0.2–0.5, which may be explained by the perturbing action of spiral arms that develop in the disk soon after its formation. The dust-to-gas mass ratio in the disk midplane is highly nonhomogeneous throughout the disk extent and is in general enhanced by a factor of several compared to the fiducial 1:100 value. Low St hinders strong dust accumulation in the spiral arms compared to the rest of the disk and the nonsteady nature of the spirals is also an obstacle. The spatial distribution of pebbles in the disk midplane exhibits a highly nonhomogeneous and patchy character. The total mass of pebbles in the disk increases with time and reaches a few tens of Earth masses after a few tens of thousand years of disk evolution. Conclusions. We found that protoplanetary disks with an age ≤20 kyr can possess notable amounts of pebbles and feature dust-togas density enhancements in the disk midplane. Hence, these young disks can already be ripe for the planet formation process to start. Multidimensional numerical models of disk formation that consider the coevolution of gas and dust including dust growth are important to improve our understanding of planet formation.

Planet formation around intermediate-mass stars

Context. The study of protoplanetary disc evolution and theories of planet formation has predominantly concentrated on solar- (and low-) mass stars since they host the majority of confirmed exoplanets. Nevertheless, the confirmation of numerous planets orbiting stars more massive than the Sun (up to ~3 M⊙) has sparked considerable interest in understanding the mechanisms involved in their formation, and thus in the evolution of their hosting protoplanetary discs. Aims. We aim to improve our knowledge of the evolution of the gaseous component of protoplanetary discs around intermediate-mass stars and to set the stage for future studies of planet formation around them. Methods. We study the long-term evolution of protoplanetary discs affected by viscous accretion and photoevaporation by X-ray and far-ultraviolet (FUV) photons from the central star around stars in the range of 1–3 M⊙, considering the effects of stellar evolution and solving the vertical structure equations of the disc. We explore the effect of different values of the viscosity parameter and the initial mass of the disc. Results. We find that the evolutionary pathway of protoplanetary disc dispersal due to photoevaporation depends on the stellar mass. Our simulations reveal four distinct evolutionary pathways for the gas component not reported before that are a consequence of stellar evolution and that likely have a substantial impact on the dust evolution, and thus on planet formation. As the stellar mass increases from one solar mass to ~1.5–2 M⊙, the evolution of the disc changes from the conventional inside-out clearing, in which X-ray photoevaporation generates inner holes, to a homogeneous disc evolution scenario where both inner and outer discs formed after a gap is opened by photoevaporation vanish over a similar timescale. As the stellar mass continues to increase, reaching ~2–3 M⊙, we identify a distinct pathway that we refer to as revenant disc evolution. In this scenario, the inner and outer discs reconnect after the gap opened. For the largest masses, we observe outside-in disc dispersal, in which the outer disc dissipates first due to a stronger FUV photoevaporation rate. Revenant disc evolution stands out as it is capable of extending the disc lifespan. Otherwise, the disc dispersal timescale decreases with increasing stellar mass except for low-viscosity discs.

Analytical solutions for the evolution of MHD wind-driven accretion discs

ABSTRACT We present new analytical solutions for the evolution of protoplanetary discs (PPDs) where magnetohydrodynamic (MHD) wind-driven processes dominate. Our study uses a 1D model which incorporates equations detailing angular momentum extraction by MHD winds and mass-loss rates. Our solutions demonstrate that the disc retains its initial state during the early phases; however, it rapidly evolves towards a self-similar state in the later stages of disc evolution. The total disc mass undergoes a continuous decline over time, with a particularly rapid reduction occurring beyond a certain critical time threshold. This gradual decrease in mass is influenced by the wind parameters and the initial surface density of the disc. In the MHD wind-dominated regime, we show that the disc’s lifespan correlates positively with the magnetic lever arm up to a certain threshold, irrespective of the initial disc size. PPDs with a larger magnetic lever arm are found to maintain significantly higher total disc mass over extended periods compared to their counterparts. The mass ejection-to-accretion ratio increases in efficient wind scenarios and is further amplified by a steeper initial surface density profile. Our analysis also reveals varied evolutionary trajectories in the plane of accretion rate and total disc mass, influenced by magnetic parameters and initial disc size. In scenarios with efficient MHD winds, discs with bigger sizes have extended operation time for mechanisms governing planet formation.