Protoplanetary Disks Research Articles

Study of the Spectral Energy Distribution Model of the Protoplanetary Disk Geometric around the Brown Dwarf CFHT-BD-Tau 4

The accretion circumstellar disk of young stars and the Brown dwarf plays an essential role in the formation and evaluation of the planet. Our main work in this paper is to investigate the geometrical shape model for the protoplanetary disk around one of the Brown Dwarfs. The photometric measurements for the brown dwarf CFHT-BD-Tau 4 were extracted from the Vizier archive. We used a numerical simulation to build a model of the spectral energy distribution of our target CFHT-BD-Tau 4. The spectral energy distribution model was fitted with observational data for the brown dwarf CFHT-BD-Tau 4. A transitional disk has been assumed around CFHT-BD-Tau 4. We obtained physical properties of the two disks and the size of the gap between them by fitting the SED. The gap in the protoplanetary disk proves that a planet formation process occurred around the Brown dwarf.

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 We compare for the first time the spatial location of the emission in the and 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 emission is co-spatial with the disk emission, whereas that of is emitted both on the disk surface and on the ionization front of the proplyds. We show for the first time that the 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 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 and 6300 lines are well detected in all these targets, the total detection rate is $ in the case of the line. This number increases substantially to a $ detection rate in sigma -Orionis, a region with higher UV radiation than in low-mass SFRs, but lower than in the ONC. The spatial location of the 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.

Presolar Grains as Probes of Supernova Nucleosynthesis.

We provide an overview of the isotopic signatures of presolar supernova grains, specifically focusing on 44Ti-containing grains with robustly inferred supernova origins and their implications for nucleosynthesis and mixing mechanisms in supernovae. Recent technique advancements have enabled the differentiation between radiogenic (from 44Ti decay) and nonradiogenic 44Ca excesses in presolar grains, made possible by enhanced spatial resolution of Ca-Ti isotope analyses with the Cameca NanoSIMS (Nano-scale Secondary Ion Mass Spectrometer) instrument. Within the context of presolar supernova grain data, we discuss (i) the production of 44Ti in supernovae and the impact of interstellar medium heterogeneities on the galactic chemical evolution of 44Ca/40Ca, (ii) the nucleosynthesis processes of neutron bursts and explosive H-burning in Type II supernovae, and (iii) challenges in identifying the progenitor supernovae for 54Cr-rich presolar nanospinel grains. Drawing on constraints and insights derived from presolar supernova grain data, we also provide an overview of our current understanding of the roles played by various supernova types- including Type II, Type Ia, and electron capture supernovae- in accounting for the diverse array of nucleosynthetic isotopic variations identified in bulk meteorites and meteoritic components. We briefly overview the potential mechanisms that have been proposed to explain these nucleosynthetic variations by describing the transport and distribution of presolar dust carriers in the protoplanetary disk. We highlight existing controversies in the interpretation of presolar grain data and meteoritic nucleosynthetic isotopic variations, while also outlining potential directions for future research.

Population of giant planets around B stars from the first part of the BEAST survey

Exoplanets form from circumstellar protoplanetary disks whose fundamental properties (notably their extent, composition, mass, temperature, and lifetime) depend on the host star properties, such as their mass and luminosity. B stars are among the most massive stars and their protoplanetary disks test extreme conditions for exoplanet formation. This paper investigates the frequency of giant planet companions around young B stars (median age of 16 Myr) in the Scorpius-Centaurus (Sco-Cen) association, the closest association containing a large population of B stars. We systematically searched for massive exoplanets with the high-contrast direct imaging instrument SPHERE using the data from the BEAST survey, which targets a homogeneous sample of young B stars from the wide Sco-Cen association. We derived accurate detection limits in the case of non-detections. We found evidence in previous papers for two substellar companions around 42 stars. The masses of these companions are straddling the sim 13 Jupiter mass deuterium burning limit, but their mass ratio with respect to their host star is close to that of Jupiter. We derived a frequency of such massive planetary-mass companions around B stars of $11_ $<!PCT!>, accounting for the survey sensitivity. The discoveries of substellar companions b and B happened after only a few stars in the survey had been observed, raising the possibility that massive Jovian planets might be common around B stars. However, our statistical analysis shows that the occurrence rate of such planets is similar around B stars and around solar-type stars of a similar age, while B-star companions exhibit low mass ratios and a larger semi-major axis.

Correction to: Early stages of gap opening by planets in protoplanetary discs

Correction to: Early stages of gap opening by planets in protoplanetary discs

Linking the primordial composition of planet building disks to the present-day composition of rocky exoplanets

The composition of rocky planets is strongly driven by the primordial materials in the protoplanetary disk, which can be inferred from the abundances of the host star. Understanding this compositional link is crucial for characterizing exoplanets. We aim to investigate the relationship between the compositions of low-mass planets and their host stars. We determined the primordial compositions of host stars using high-precision present-day stellar abundances and stellar evolutionary models. These primordial abundances were then input into a stoichiometric model to estimate the composition of planet-building blocks. Additionally, we employed a three-component planetary interior model (core, mantle, and water in different phases) to estimate planetary compositions based only on their radius and mass. We find that although stellar abundances vary over time, relevant abundance ratios such as Fe/Mg remain relatively constant during the main sequence evolution for low temperature stars. A strong correlation is found between the iron-to-silicate mass fraction of protoplanetary disks and planets, while no significant correlation was observed for water mass fractions. The Fe/Mg ratio varies significantly between planets and their stars, indicating substantial disk-driven compositional diversity, and this ratio also correlates with planetary radius. While stellar abundances, as a proxy of the composition of protoplanetary disk, provide a baseline for planetary composition, significant deviations arise due to complex disk processes, challenging the assumption of a direct, one-to-one elemental relationship between stars and their planets.

Measurement of interstellar extinction for classical T Tauri stars using far-UV H_2 line fluxes

Understanding the interstellar and potentially circumstellar extinction in the sight lines of classical T Tauri stars is an important ingredient for constructing reliable spectral energy distributions, which catalyze protoplanetary disk chemistry, for example. Therefore, some attempts of measuring $A_ V $ toward individual stars have been made using partly different wavelength regimes and different underlying assumptions. We used strong lines of Lyalpha fluorescent H$_ and derived the extinction based on the assumption of optically thin transitions. We investigated a sample of 72 classical T Tauri stars observed with the Hubble Space Telescope in the framework of the ULLYSES program. We computed $A_V$ and $R_V$ values for the 34 objects with sufficient data quality and an additionally $A_V$ value for the canonical $R_V=3.1$ value. Our results agree largely with values obtained from optical data. Moreover, we confirm the degeneracy between $A_ V $ and $R_ V $ and present possibilities to break this. Finally, we discuss whether the assumption of optical thin lines is valid.

Dust mass in protoplanetary disks with porous dust opacities

Atacama Large Millimeter/submillimeter Array surveys have suggested that protoplanetary disks are not massive enough to form the known exoplanet population, based on the assumption that the millimeter continuum emission is optically thin. In this work, we investigate how the mass determination is influenced when the porosity of dust grains is considered in radiative transfer models. The results show that disks with porous dust opacities yield similar dust temperatures, but systematically lower millimeter fluxes, as compared to disks that incorporate compact dust grains. Moreover, we have recalibrated the relation between dust temperature and stellar luminosity for a wide range of stellar parameters. We also calculated the dust masses of a large sample of disks using the traditionally analytic approach. The median dust mass from our calculation is about six times higher than the literature result, and this is mostly driven by the different opacities of porous and compact grains. A comparison of the cumulative distribution function between disk dust masses and exoplanet masses shows that the median exoplanet mass is about two times lower than the median dust mass when grains are assumed to be porous and there are no exoplanetary systems with masses higher than the most massive disks. Our analysis suggests that adopting porous dust opacities may alleviate the mass budget problem for planet formation. As an example illustrating the combined effects of optical depth and porous dust opacities on the mass estimation, we conducted new IRAM/NIKA-2 observations toward the IRAS\,04370+2559 disk and performed a detailed radiative transfer modeling of the spectral energy distribution (SED). The best-fit dust mass is roughly 100 times higher than the value given by a traditionally analytic calculation. Future spatially resolved observations at various wavelengths are required to better constrain the dust mass.

Evidence for Magnetically‐Driven Accretion in the Distal Solar System

AbstractPaleomagnetic measurements of meteorites indicate that magnetic fields existed in the inner solar nebula capable of driving accretion at rates similar to those observed for young stellar objects with protoplanetary disks. However, the field strength in the solar system beyond ∼7 astronomical units (AU) and its role in accretion remain poorly constrained. Returned samples from asteroid (162173) Ryugu offer the possibility of determining the nebular field intensity in this distal region. Here, we report paleomagnetic studies of three Ryugu particles which reveal that alteration occurred in the presence of a null or relatively weak (<15.8 μT) field within 3 million years (Ma) after solar system formation. This resolves previously contrasting reports that Ryugu's parent body experienced alteration in the presence of a strong (>80 μT) magnetic field and weak or null field (<3 μT). In addition, we re‐examine previous paleomagnetic and Mn‐Cr chronometry studies of three other distally‐sourced meteorites, Tagish Lake, Tarda, and Wisconsin Range 91600, which measured paleointensities of <0.9, <1.7 and 5.1 ± 4.5 μT respectively. While it was previously unclear whether these records were acquired while the nebula was present, our re‐analysis suggests that their records are sufficiently old (i.e., <3.5 Ma after solar system formation) to be nebular in origin. Collectively, these data demonstrate that the distal solar system nebular field, while faint, was likely still strong enough to drive accretion at rates like those observed in the inner solar system.

REsolved ALMA and SMA Observations of Nearby Stars. REASONS: A population of 74 resolved planetesimal belts at millimetre wavelengths

Planetesimal belts are ubiquitous around nearby stars, and their spatial properties hold crucial information for planetesimal and planet formation models. We present resolved dust observations of 74 planetary systems as part of the REsolved ALMA and SMA Observations of Nearby Stars (REASONS) survey and archival reanalysis. We uniformly modelled interferometric visibilities for the entire sample to obtain the basic spatial properties of each belt, and combined these with constraints from multi-wavelength photometry. We report key findings from a first exploration of this legacy dataset: (1) Belt dust masses are depleted over time in a radially dependent way, with dust being depleted faster in smaller belts, as predicted by collisional evolution. (2) Most belts are broad discs rather than narrow rings, with much broader fractional widths than rings in protoplanetary discs. We link broad belts to either unresolved substructure or broad planetesimal discs produced if protoplanetary rings migrate. (3) The vertical aspect ratios ($h=H/R$) of 24 belts indicate orbital inclinations of sim 1-20$^ circ $, implying relative particle velocities of sim 0.1-4 km/s, and no clear evolution of heights with system age. This could be explained by early stirring within the belt by large bodies (with sizes of at least sim 140 km to the size of the Moon), by inheritance of inclinations from the protoplanetary disc stage, or by a diversity in evolutionary pathways and gravitational stirring mechanisms. We release the REASONS legacy multidimensional sample of millimetre-resolved belts to the community as a valuable tool for follow-up multi-wavelength observations and population modelling studies.

Planet migration in windy discs

ABSTRACT Accretion of protoplanetary discs (PPDs) could be driven by magnetohydrodynamic disc winds rather than turbulent viscosity. With a dynamical prescription for angular momentum transport induced by disc winds, we perform 2D simulations of PPDs to systematically investigate the rate and direction of planet migration in a windy disc. We find that the the strength of disc winds influences the corotation region similarly to the ‘desaturation’ effect of viscosity. The magnitude and direction of torque depend sensitively on the hierarchy between the radial advection time-scale across the horseshoe due to disc wind $\tau _{\rm dw}$, the horseshoe libration time-scale $\tau _{\rm lib}$ and U-turn time-scale $\tau _{\rm U-turn}$. Initially, as wind strength increases and the advection time-scale shortens, a non-linear horseshoe drag emerges when $\tau _{\rm dw} \lesssim \tau _{\rm lib}$, which tends to drive strong outward migration. Subsequently, the drag becomes linear and planets typically still migrate inward when $\tau _{\rm dw} \lesssim \tau _{\rm U-turn} \sim \tau _{\rm lib}h$, where h is the disc aspect ratio. For a planet with mass ratio of ${\sim} 10^{-5}$, the zone of outward migration sandwiched between inner and outer inward migration zones corresponds to $\sim$10–100 au in a PPD with accretion rates between $10^{-8}$ and $10^{-7}\, \mathrm{ M}_\odot \text{yr}^{-1}$.

Effect of Time-varying X-Ray Emission from Stellar Flares on the Ionization of Protoplanetary Disks

X-rays have significant impacts on cold, weakly ionized protoplanetary disks by increasing the ionization rate and driving chemical reactions. Stellar flares are explosions that emit intense X-rays and are the unique source of hard X-rays with an energy of ≳10 keV in the protoplanetary disk systems. Hard X-rays should be carefully taken into account in models as they can reach the disk midplane as a result of scattering in the disk atmospheres. However, previous models are insufficient to predict the hard X-ray spectra because of simplifications in flare models. We develop a model of X-ray spectra of stellar flares based on observations and flare theories. The flare temperature and nonthermal electron emissions are modeled as functions of flare energy, which allows us to better predict the hard X-ray photon flux than before. Using our X-ray model, we conduct radiative transfer calculations to investigate the impact of flare hard X-rays on disk ionization, with a particular focus on the protoplanetary disk around a T Tauri star. We demonstrate that for a flare with an energy of 1035 erg, X-ray photons with ≳5 keV increase the ionization rates more than galactic cosmic rays down to z ≈ 0.1R. The contribution of flare X-rays to the ionization at the midplane depends on the disk parameters such as disk mass and dust settling. We also find that the 10 yr averaged X-rays from multiple flares could certainly contribute to the ionization. These results emphasize the importance of stellar flares on the disk evolution.

Dust ring and gap formation by gas flow induced by low-mass planets embedded in protoplanetary disks. II. Time-dependent model

The observed dust rings and gaps in protoplanetary disks could be imprints of forming planets. Even low-mass planets in the 1--10 Earth-mass regime, which have not yet carved deep gas gaps, can generate such dust rings and gaps by driving a radially-outward gas flow, as shown in previous work. However, understanding the creation and evolution of these dust structures is challenging due to dust drift and diffusion, requiring an approach beyond previous steady state models. Here we investigate the time evolution of the dust surface density influenced by the planet-induced gas flow based on post-processing three-dimensional hydrodynamical simulations. We find that planets larger than a dimensionless thermal mass of $m=0.05$ corresponding to $0.3$ Earth mass at 1 au or 1.7 Earth masses at 10 au generate dust rings and gaps provided that solids have small Stokes numbers St and that the disk midplane is weakly turbulent ($ diff As dust particles pile up outside the orbit of the planet, the interior gap expands with time when the advective flux dominates over diffusion. Dust gap depths range from a factor of a few to several orders of magnitude, depending on planet mass and the level of midplane particle diffusion. We constructed a semi-analytic model describing the width of the dust ring and gap, and then compared it with the observational data. We find that up to 65<!PCT!> of the observed wide-orbit gaps could be explained as resulting from the presence of a low-mass planet, assuming diff and St However, it is more challenging to explain the observed wide rings, which in our model would require the presence of a population of small particles ($ Further work is needed to explore the role of pebble fragmentation, planet migration, and the effect of multiple planets.

High stability of charged particle clusters in protoplanetary disks

Context. The initial particle growth in protoplanetary disks is limited by a bouncing barrier at submillimeter wavelengths. Bouncing leads to tribocharging and the electrostatic attraction of tribocharged aggregates may eventually draw them into large clusters. A charge- mediated growth phase allows for the formation of larger entities, namely, clusters of aggregates that are more prone to further particle concentrations, such as the streaming instability. Aims. We aim to quantify the strength of the electrostatic forces. Methods. In laboratory experiments, we used an acoustic trap to levitate small aggregates of tribocharged submm grains. These aggregates spin up within the trap until they lose grains. Thus, we used the centrifugal force as a measure of the local force. Results. Grains are regularly bound strongly to their neighbors. In comparison, the force at ejection can be stronger than the attractive scattering forces of the trap and can therefore be several orders of magnitude larger than expected. We note that these forces are long- ranging, compared to van der Waals forces. Thus, charged aggregates are much more stable than uncharged ones. Conclusions. Particle aggregates in disks might grow to centimeter clusters or larger as tribocharging increases the effective binding forces. This allows for hydrodynamic concentration and planetesimal formation to eventually take place throughout a wide part of the disk.

Light Echoes of Protoplanetary Disks

Light echoes offer a means of studying protoplanetary disks, including their geometry and composition, even when they are not spatially resolved. We present a test of this approach applied specifically to optically thick, geometrically flared disks around active stars. Here we adopt stellar parameters of an active M dwarf to calculate light echoes for disks and rings with radii that would produce time delays consistent with TESS short cadence (about 2 minutes) time bins. Our results show successful fits to disk parameters, highlighting the potential effectiveness of this method in the search for protoplanetary disks.

Mind the Kinematics Simulation of Planet–Disk Interactions: Time Evolution and Numerical Resolution

Planet–disk interactions can produce kinematic signatures in protoplanetary disks. While recent observations have detected non-Keplerian gas motions in disks, their origins are still being debated. To explore this, we conduct 3D hydrodynamic simulations using the code FARGO3D to study nonaxisymmetric kinematic perturbations at two scale heights induced by Jovian planets in protoplanetary disks, followed by examinations of detectable signals in synthetic CO emission line observations at millimeter wavelengths. We advocate for using residual velocity or channel maps, generated by subtracting an azimuthally averaged background of the disk, to identify planet-induced kinematic perturbations. We investigate the effects of two basic simulation parameters, simulation duration and numerical resolution, on the simulation results. Our findings suggest that a short simulation (e.g., 100 orbits) is insufficient to establish a steady velocity pattern given our chosen viscosity (α = 10−3) and displays plenty of fluctuations on an orbital timescale. Such transient features could be detected in observations. By contrast, a long simulation (e.g., 1000 orbits) is required to reach steady state in kinematic structures. At 1000 orbits, the strongest detectable velocity structures are found in the spiral wakes close to the planet. Through numerical convergence tests, we find hydrodynamics results converge in spiral regions at a resolution of 14 cells per disk scale height or higher. Meanwhile, synthetic observations produced from hydrodynamic simulations at different resolutions are indistinguishable with 0.″1 angular resolution and 10 hr of integration time on Atacama Large Millimeter/submillimeter Array.

Dust Drift Timescales in Protoplanetary Disks at the Cusp of Gravitational Instability

Millimeter emitting dust grains have sizes that make them susceptible to drift in protoplanetary disks due to the difference between their orbital speed and that of the gas. The characteristic drift timescale depends on the surface density of the gas. By comparing disk radius measurements from Atacama Large Millimeter/submillimeter Array CO and continuum observations at millimeter wavelengths, the gas surface density profile and dust drift time can be self-consistently determined. We find that profiles which match the measured dust mass have very short drift timescales, an order of magnitude or more shorter than the stellar age, whereas profiles for disks that are on the cusp of gravitational instability, defined via the minimum value of the Toomre parameter, Qmin∼1−2 , have drift timescales comparable to the stellar lifetime. This holds for disks with masses of dust ≳5 M ⊕ across a range of absolute ages from less than 1 Myr to over 10 Myr. The inferred disk masses scale with stellar mass as Mdisk≈M*/5Qmin . This interpretation of the gas and dust disk sizes simultaneously solves two long standing issues regarding the dust lifetime and exoplanet mass budget, and suggests that we consider millimeter wavelength observations as a window into an underlying population of particles with a wide size distribution in secular evolution with a massive planetesimal disk.

Shadows Wreak Havoc in Transition Disks

We demonstrate that shadows cast on a protoplanetary disk can drive it eccentric. Stellar irradiation dominates heating across much of these disks, so an uneven illumination can have interesting dynamical effects. Here, we focus on transition disks. We carry out 3D Athena++ simulations, using a constant thermal relaxation time to describe the disk’s response to changing stellar illumination. We find that an asymmetric shadow, a feature commonly observed in real disks, perturbs the radial pressure gradient and distorts the fluid streamlines into a set of twisted ellipses. Interactions between these streamlines have a range of consequences. For a narrow ring, an asymmetric shadow can sharply truncate its inner edge, possibly explaining the steep density drop-offs observed in some disks and obviating the need for massive perturbers. For a wide ring, such a shadow can dismantle it into two (or possibly more) eccentric rings. These rings continuously exert torque on each other and drive gas accretion at a healthy rate, even in the absence of disk viscosity. Signatures of such twisted eccentric rings may have already been observed as, e.g., twisted velocity maps inside gas cavities. We advocate for more targeted observations and for a better understanding on the origin of such shadows.

ALMA Observations of Proper Motions of the Dust Clumps in the Protoplanetary Disk MWC 758

To study the dust dynamics in the dust-trapping vortices in the protoplanetary disk around MWC 758, we analyzed the 1.3 mm continuum images of the MWC 758 disk obtained with the Atacama Large Millimeter/submillimeter Array in 2017 and 2021. We detect proper motions of 22 mas and 24 mas in the two dust clumps at radii of 0.″32 and 0.″54 in the disk on the plane of the sky, respectively. On the assumption that the dust clumps are located in the disk midplane, the velocities of the observed proper motions along the azimuthal direction of the inner and outer dust clumps are sub- and super-Keplerian, respectively, and both have angular velocities corresponding to the Keplerian angular velocity at a radius of 0.″46 ± 0.″04. This deviation from the Keplerian motion is not expected in the conventional theory of vortices formed by the Rossby wave instability. The observed non-Keplerian proper motions of the dust clumps are unlikely due to the disk warp and eccentricity, nor are they associated with any predicted planets. The two dust clumps are likely spatially coincident with the infrared spirals. In addition, we detect the changes in the intensity profiles of the dust clumps over the 4 yr span. Therefore, we suggest that the observed proper motions are possibly due to changes in the density distributions in the dust clumps caused by their interaction with the spirals in the disk.

Effects of Dust Coagulation on Streaming Instability

Streaming instability (SI) in dust has long been thought to be a promising process in triggering planetesimal formation in the protoplanetary disks (PPDs). In this study, we present the first numerical investigation that models the SI in the vertically stratified disk together with the dust coagulation process. Our simulations reveal that, even with the initially small dust sizes, because dust coagulation promotes dust size growth, SI can eventually still be triggered. Specifically, dust coagulation, limited only by dust fragmentation, broadens the parameter boundaries obtained from previous SI studies using single dust species. We describe the various stages of dust dynamics along with their size evolution and explore the impact of different dust fragmentation velocities. Implications of these results for realistic PPDs are also discussed.