No Arabic abstract
Spatial distribution and growth of dust in a clumpy protoplanetary disk subject to vigorous gravitational instability and fragmentation is studied numerically with sub-au resolution using the FEOSAD code. Hydrodynamics equations describing the evolution of self-gravitating and viscous protoplanetary disks in the thin-disk limit were modified to include a dust component consisting of two parts: sub-micron-sized dust and grown dust with a variable maximum radius. The conversion of small to grown dust, dust growth, friction of dust with gas, and dust self-gravity were also considered. We found that the disk appearance is notably time-variable with spiral arms, dusty rings, and clumps, constantly forming, evolving, and decaying. As a consequence, the total dust-to-gas mass ratio is highly non-homogeneous throughout the disk extent, showing order-of-magnitude local deviations from the canonical 1:100 value. Gravitationally bound clumps formed through gravitational fragmentation have a velocity pattern that deviates notably from the Keplerian rotation. Small dust is efficiently converted into grown dust in the clump interiors, reaching a maximum radius of several decimeters. Concurrently, grown dust drifts towards the clump center forming a massive compact central condensation (70-100 $M_oplus$). We argue that protoplanets may form in the interiors of inward migrating clumps before they disperse through the action of tidal torques. We foresee the formation of protoplanets at orbital distances of several tens of au with initial masses of gas and dust in the protoplanetary seed in the (0.25-1.6) $M_{rm Jup}$ and (1.0-5.5) $M_oplus$ limits, respectively. The final masses of gas and dust in the protoplanets may however be much higher due to accretion from surrounding massive metal-rich disks/envelopes.
ALMA observations of protoplanetary disks confirm earlier indications that there is a clear difference between the dust and gas radial extents. The origin of this difference is still debated, with both radial drift of the dust and optical depth effects suggested in the literature. In this work, the feedback of realistic dust particle distributions onto the gas chemistry and molecular emissivity is investigated, with a particular focus on CO isotopologues. The radial dust grain size distribution is determined using dust evolution models that include growth, fragmentation and radial drift. A new version of the code DALI is used to take into account how dust surface area and density influence the disk thermal structure, molecular abundances and excitation. The difference of dust and gas radial sizes is largely due to differences in the optical depth of CO lines and millimeter continuum, without the need to invoke radial drift. The effect of radial drift is primarily visible in the sharp outer edge of the continuum intensity profile. The gas outer radius probed by $^{12}$CO emission can easily differ by a factor of $sim 2$ between the models for a turbulent $alpha$ ranging between typical values. Grain growth and settling concur in thermally decoupling the gas and dust components, due to the low collision rate with large grains. As a result, the gas can be much colder than the dust at intermediate heights, reducing the CO excitation and emission, especially for low turbulence values. Also, due to disk mid-plane shadowing, a second CO thermal desorption (rather than photodesorption) front can occur in the warmer outer mid-plane disk. The models are compared to ALMA observations of HD 163296 as a test case. In order to reproduce the observed CO snowline of the system, a binding energy for CO typical of ice mixtures needs to be used rather than the lower pure CO value.
The central region of a circumstellar disk is difficult to resolve in global numerical simulations of collapsing cloud cores, but its effect on the evolution of the entire disk can be significant. We use numerical hydrodynamics simulations to model the long-term evolution of self-gravitating and viscous circumstellar disks in the thin-disk limit. Simulations start from the gravitational collapse of prestellar cores of 0.5--1.0~$M_odot$ and both gaseous and dusty subsystems were considered, including a model for dust growth. The inner unresolved 1.0 au of the disk is replaced with a central smart cell (CSC) -- a simplified model that simulates physical processes that may occur in this region. We found that the mass transport rate through the CSC has an appreciable effect on the evolution of the entire disk. Models with slow mass transport form more massive and warmer disks and they are more susceptible to gravitational instability and fragmentation, including a newly identified episodic mode of disk fragmentation in the T Tauri phase of disk evolution. Models with slow mass transport through the CSC feature episodic accretion and luminosity bursts in the early evolution, while models with fast transport are characterized by a steadily declining accretion rate with low-amplitude flickering. Dust grows to a larger, decimeter size in the slow transport models and efficiently drifts in the CSC, where it accumulates reaching the limit when streaming instability becomes operational. We argue that gravitational instability, together with streaming instability likely operating in the inner disk regions, constitute two concurrent planet-forming mechanisms, which may explain the observed diversity of exoplanetary orbits (Abridged).
Far-infrared and (sub)millimeter fluxes can be used to study dust in protoplanetary disks, the building blocks of planets. Here, we combine observations from the Herschel Space Observatory with ancillary data of 284 protoplanetary disks in the Taurus, Chamaeleon I, and Ophiuchus star-forming regions, covering from the optical to mm/cm wavelengths. We analyze their spectral indices as a function of wavelength and determine their (sub)millimeter slopes when possible. Most disks display observational evidence of grain growth, in agreement with previous studies. No correlation is found between other tracers of disk evolution and the millimeter spectral indices. A simple disk model is used to fit these sources, and we derive posterior distributions for the optical depth at 1.3 mm and 10 au, the disk temperature at this same radius, and the dust opacity spectral index. We find the fluxes at 70 microns to correlate strongly with disk temperatures at 10 au, as derived from these simple models. We find tentative evidence for spectral indices in Chamaeleon I being steeper than those of disks in Taurus/Ophiuchus, although more millimeter observations are needed to confirm this trend and identify its possible origin. Additionally, we determine the median spectral energy distribution of each region and find them to be similar across the entire wavelength range studied, possibly due to the large scatter in disk properties and morphologies.
We present 870 $mu$m ALMA observations of polarized dust emission toward the Class II protoplanetary disk IM Lup. We find that the orientation of the polarized emission is along the minor axis of the disk, and that the value of the polarization fraction increases steadily toward the center of the disk, reaching a peak value of ~1.1%. All of these characteristics are consistent with models of self-scattering of submillimeter-wave emission from an optically thin inclined disk. The distribution of the polarization position angles across the disk reveals that while the average orientation is along the minor axis, the polarization orientations show a significant spread in angles; this can also be explained by models of pure scattering. We compare the polarization with that of the Class I/II source HL Tau. A comparison of cuts of the polarization fraction across the major and minor axes of both sources reveals that IM Lup has a substantially higher polarization fraction than HL Tau toward the center of the disk. This enhanced polarization fraction could be due a number of factors, including higher optical depth in HL Tau, or scattering by larger dust grains in the more evolved IM Lup disk. However, models yield similar maximum grain sizes for both HL Tau (72 $mu$m) and IM Lup (61 $mu$m, this work). This reveals continued tension between grain-size estimates from scattering models and from models of the dust emission spectrum, which find that the bulk of the (unpolarized) emission in disks is most likely due to millimeter (or even centimeter) sized grains.
We estimate the mass loss rates of photoevaporative winds launched from the outer edge of protoplanetary discs impinged by an ambient radiation field. We focus on mild/moderate environments (the number of stars in the group/cluster is N ~ 50), and explore disc sizes ranging between 20 and 250 AU. We evaluate the steady-state structures of the photoevaporative winds by coupling temperature estimates obtained with a PDR code with 1D radial hydrodynamical equations. We also consider the impact of dust dragging and grain growth on the final mass loss rates. We find that these winds are much more significant than have been appreciated hitherto when grain growth is included in the modelling: in particular, mass loss rates > 1e-8 M_sun/yr are predicted even for modest background field strengths ( ~ 30 G_0) in the case of discs that extend to R > 150 AU. Grain growth significantly affects the final mass loss rates by reducing the average cross section at FUV wavelengths, and thus allowing a much more vigorous flow. The radial profiles of observable quantities (in particular surface density, temperature and velocity patterns) indicate that these winds have characteristic features that are now potentially observable with ALMA. In particular, such discs should have extended gaseous emission that is dust depleted in the outer regions, characterised by a non-Keplerian rotation curve, and with a radially increasing temperature gradient.