No Arabic abstract
We assess the ionising effect of low energy protostellar cosmic rays in protoplanetary disks around a young solar mass star for a wide range of disk parameters. We assume a source of low energy cosmic rays located close to the young star which travel diffusively through the protoplanetary disk. We use observationally inferred values from nearby star-forming regions for the total disk mass and the radial density profile. We investigate the influence of varying the disk mass within the observed scatter for a solar mass star. We find that for a large range of disk masses and density profiles that protoplanetary disks are optically thin to low energy ($sim$3 GeV) cosmic rays. At $Rsim10$au, for all of the disks that we consider ($M_mathrm{disk}=6.0times10^{-4} - 2.4times 10^{-2}M_odot$), the ionisation rate due to low energy stellar cosmic rays is larger than that expected from unmodulated galactic cosmic rays. This is in contrast to our previous results which assumed a much denser disk which may be appropriate for a more embedded source. At $Rsim70$au, the ionisation rate due to stellar cosmic rays dominates in $sim$50% of the disks. These are the less massive disks with less steep density profiles. At this radius there is at least an order of magnitude difference in the ionisation rate between the least and most massive disk that we consider. Our results indicate, for a wide range of disk masses, that low energy stellar cosmic rays provide an important source of ionisation at the disk midplane at large radii ($sim$70au).
Protoplanetary disks around young stars are the sites of planet formation. While the dust mass can be estimated using standard methods, determining the gas mass - and thus the amount of material available to form giant planets - has proven to be very difficult. Hydrogen deuteride (HD) is a promising alternative to the commonly-used gas mass tracer, CO. We aim to examine the robustness of HD as tracer of the disk gas mass, specifically the effect of gas mass on the HD FIR emission and its sensitivity to the vertical structure. Deuterium chemistry reactions relevant for HD were implemented in the thermochemical code DALI and models were run for a range of disk masses and vertical structures. The HD J=1-0 line intensity depends directly on the gas mass through a sublinear power law relation with a slope of ~0.8. Assuming no prior knowledge about the vertical structure of a disk and using only the HD 1-0 flux, gas masses can be estimated to within a factor of 2 for low mass disks (M$_{rm disk} < 10^{-3}$ M$_odot$). For more massive disks, this uncertainty increases to more than an order of magnitude. Adding the HD 2-1 line or independent information about the vertical structure can reduce this uncertainty to a factor of ~3 for all disk masses. For TW Hya, using the radial and vertical structure from Kama et al. 2016b the observations constrain the gas mass to $6cdot10^{-3}$ M$_odot$ < M$_{rm disk} < 9cdot10^{-3}$ M$_odot$. Future observations require a 5$sigma$ sensitivity of $1.8cdot10^{-20}$ W m$^{-2}$ ($2.5cdot10^{-20}$ W m$^{-2}$) and a spectral resolving power R > 300 (1000) to detect HD 1-0 (HD 2-1) for all disk masses above $10^{-5}$ M$_odot$ with a line-to-continuum ratio > 0.01. These results show that HD can be used as an independent gas mass tracer with a relatively low uncertainty and should be considered as an important science goal for future FIR missions.
The volatile composition of a planet is determined by the inventory of gas and ice in the parent disk. The volatile chemistry in the disk is expected to evolve over time, though this evolution is poorly constrained observationally. We present ALMA observations of C18O, C2H, and the isotopologues H13CN, HC15N, and DCN towards five Class 0/I disk candidates. Combined with a sample of fourteen Class II disks presented in Bergner et al. (2019b), this data set offers a view of volatile chemical evolution over the disk lifetime. Our estimates of C18O abundances are consistent with a rapid depletion of CO in the first ~0.5-1 Myr of the disk lifetime. We do not see evidence that C2H and HCN formation are enhanced by CO depletion, possibly because the gas is already quite under-abundant in CO. Further CO depletion may actually hinder their production by limiting the gas-phase carbon supply. The embedded sources show several chemical differences compared to the Class II stage, which seem to arise from shielding of radiation by the envelope (impacting C2H formation and HC15N fractionation) and sublimation of ices from infalling material (impacting HCN and C18O abundances). Such chemical differences between Class 0/I and Class II sources may affect the volatile composition of planet-forming material at different stages in the disk lifetime.
Both observations of arc-like structures and luminosity bursts of stars > 1 Myr in age indicate that at least some stars undergo late infall events. We investigate scenarios of replenishing the mass reservoir around a star via capturing and infalling events of cloudlets. We carry out altogether 24 three-dimensional hydrodynamical simulations of cloudlet encounters with a Herbig star of mass 2.5 solar mass using the moving mesh code AREPO. To account for the two possibilities of a star or a cloudlet traveling through the interstellar medium (ISM), we put either the star or the cloudlet at rest with respect to the background gas. For absent cooling in the adiabatic runs, almost none of the cloudlet gas is captured due to high thermal pressure. However, second-generation disks easily form when accounting for cooling of the gas. The disk radii range from several 100 au to about 1000 au and associated arc-like structures up to 10 000 au in length form around the star for runs with and without stellar irradiation. Consistent with angular momentum conservation, the arcs and disks are larger for larger impact parameters. Accounting for turbulence in the cloudlet only mildly changes the model outcome. In the case of the star being at rest with the background gas, the disk formation and mass replenishment process is more pronounced and the associated arc-shaped streamers are longer-lived. The results of our models confirm that late encounter events lead to the formation of transitional disks associated with arc-shaped structures such as observed for AB Aurigae or HD 100546. In addition, we find that second-generation disks and their associated filamentary arms are longer lived (>100 000 yrs) in infall events, when the star is at rest with the background gas.
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.
We report FUV, optical, and NIR observations of three T Tauri stars in the Orion OB1b subassociation with H$alpha$ equivalent widths consistent with low or absent accretion and various degrees of excess flux in the mid-infrared. We aim to search for evidence of gas in the inner disk in HST ACS/SBC spectra, and to probe the accretion flows onto the star using H$alpha$ and He I $lambda$10830 in spectra obtained at the Magellan and SOAR telescopes. At the critical age of 5 Myr, the targets are at different stages of disk evolution. One of our targets is clearly accreting, as shown by redshifted absorption at free-fall velocities in the He I line and wide wings in H$alpha$; however, a marginal detection of FUV H$_2$ suggests that little gas is present in the inner disk, although the spectral energy distribution indicates that small dust still remains close to the star. Another target is surrounded by a transitional disk, with an inner cavity in which little sub-micron dust remains. Still, the inner disk shows substantial amounts of gas, accreting onto the star at a probably low, but uncertain rate. The third target lacks both a He I line or FUV emission, consistent with no accretion or inner gas disk; its very weak IR excess is consistent with a debris disk. Different processes occurring in targets with ages close to the disk dispersal time suggest that the end of accretion phase is reached in diverse ways.