No Arabic abstract
Computing the flow from externally FUV irradiated protoplanetary discs requires solving complicated and expensive photodissociation physics iteratively in conjunction with hydrodynamics. Previous studies have therefore been limited to 1D models of this process. In this paper we compare 2D-axisymmetric models of externally photoevaporating discs with their 1D analogues, finding that mass loss rates are consistent to within a factor four. The mass loss rates in 2D are higher, in part because half of the mass loss comes from the disc surface (which 1D models neglect). 1D mass loss rates used as the basis for disc viscous evolutionary calculations are hence expected to be conservative. We study the anatomy of externally driven winds including the streamline morphology, kinematic, thermal and chemical structure. A key difference between the 1D and 2D models is in the chemical abundances. For instance in the 2D models CO can be dissociated at smaller radial distances from the disc outer edge than in 1D calculations because gas is photodissociated by radiation along trajectories that are assumed infinitely optically thick in 1D models. Multidimensional models will hence be critical for predicting observable signatures of environmentally photoevaporating protoplanetary discs.
There is growing theoretical and observational evidence that protoplanetary disc evolution may be significantly affected by the canonical levels of far ultraviolet (FUV) radiation found in a star forming environment, leading to substantial stripping of material from the disc outer edge even in the absence of nearby massive stars. In this paper we perform the first full radiation hydrodynamic simulations of the flow from the outer rim of protoplanetary discs externally irradiated by such intermediate strength FUV fields, including direct modelling of the photon dominated region (PDR) which is required to accurately compute the thermal properties. We find excellent agreement between our models and the semi-analytic models of Facchini et al. (2016) for the profile of the flow itself, as well as the mass loss rate and location of their critical radius. This both validates their results (which differed significantly from prior semi-analytic estimates) and our new numerical method, the latter of which can now be applied to elements of the problem that the semi--analytic approaches are incapable of modelling. We also obtain the composition of the flow, but given the simple geometry of our models we can only hint at some diagnostics for future observations of externally irradiated discs at this stage. We also discuss the potential for these models as benchmarks for future photochemical-dynamical codes.
We present an open access grid of 3930 calculations of externally evaporating protoplanetary discs. This spans a range of disc sizes (1-400AU), disc masses, UV field strengths (10-10$^4$G$_0$) and stellar masses (0.05-1.9M$_odot$). The grid is publicly available for download, and offers a means of cheaply including external photoevaporation in disc evolutionary calculations. It can also be queried using an online tool for quick estimates of instantaneous mass loss rates (e.g for convenient evaluation of real observed systems). The `FRIED grid itself illustrates that for discs around stars $leq0.3$M$_odot$ external photoevaporation is effective down to small radii ($<50$AU) down to UV fields at least as weak as 10G$_0$. At the other end of the scale, in a $10^4$G$_0$ environment photoevaporation is effective down to 1AU even for stellar masses at least as high as 1.9M$_odot$. We also illustrate in which regimes CO survives in the photoevaporative outflow for significant mass loss rates; marking a system a good candidate to detect external photoevaporation in weak-intermediate UV environments through sub-Keplerian rotation. Finally we make illustrative mass loss rate estimates for discs in Taurus based on the Guilloteau et al. (2011) star-disc parameters, finding that around half are expected to have both significant mass loss and retain CO in the photoevaporative outflow.
The stellar birth environment can significantly shorten protoplanetary disc (PPD) lifetimes due to the influence of stellar feedback mechanisms. The degree to which these mechanisms suppress the time and mass available for planet formation is dependent on the local far-ultraviolet (FUV) field strength, stellar density, and ISM properties. In this work, we present the first theoretical framework quantifying the distribution of PPD dispersal time-scales as a function of parameters that describe the galactic environment. We calculate the probability density function for FUV flux and stellar density in the solar neighbourhood. In agreement with previous studies, we find that external photoevaporation is the dominant environment-related factor influencing local stellar populations after the embedded phase. Applying our general prescription to the Central Molecular Zone of the Milky Way (i.e. the central ~250 pc), we predict that 90% of PPDs in the region are destroyed within 1 Myr of the dispersal of the parent molecular cloud. Even in such dense environments, we find that external photoevaporation is the dominant disc depletion mechanism over dynamical encounters between stars. PPDs around low-mass stars are particularly sensitive to FUV-induced mass loss, due to a shallower gravitational potential. For stars of mass ~1 $M_odot$, the solar neighbourhood lies at approximately the highest gas surface density for which PPD dispersal is still relatively unaffected by external FUV photons, with a median PPD dispersal timescale of ~4 Myr. We highlight the key questions to be addressed to further contextualise the significance of the local galactic environment for planet formation.
We present the first ALMA survey of protoplanetary discs at 3 mm, targeting 36 young stellar objects in the Lupus star-forming region with deep observations (sensitivity 20-50 microJy/beam) at ~0.35 resolution (~50 au). Building on previous ALMA surveys at 0.89 and 1.3 mm that observed the complete sample of Class II discs in Lupus at a comparable resolution, we aim to assess the level of grain growth in the relatively young Lupus region. We measure 3 mm integrated fluxes, from which we derive disc-averaged 1-3 mm spectral indices. We find that the mean spectral index of the observed Lupus discs is $alpha_mathrm{1-3 mm}=2.23pm0.06$, in all cases $alpha_mathrm{1-3 mm}<3.0$, with a tendency for larger spectral indices in the brightest discs and in transition discs. Furthermore, we find that the distribution of spectral indices in Lupus discs is statistically indistinguishable from that of the Taurus and Ophiuchus star-forming regions. Assuming the emission is optically thin, the low values $alpha_mathrm{1-3 mm}leq 2.5$ measured for most discs can be interpreted with the presence of grains larger than 1 mm. The observations of the faint discs in the sample can be explained without invoking the presence of large grains, namely through a mixture of optically thin and optically thick emission from small grains. However, the bright (and typically large) discs do inescapably require the presence of millimeter-sized grains in order to have realistic masses. Based on a disc mass argument, our results challenge previous claims that the presence of optically thick sub-structures may be a universal explanation for the empirical millimeter size-luminosity correlation observed at 0.89 mm.
Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared scattered light and (sub-) millimeter continuum and CO line emission. Using representative prograde and retrograde encounters for direct comparison, we post-process hydrodynamical simulations with radiative transfer methods to generate a catalogue of multi-wavelength observations. This provides a reference to identify flybys in recent near-infrared and sub-millimetre observations (e.g., RW Aur, AS 205, HV Tau & DO Tau, FU Ori, V2775 Ori, and Z CMa).