The Herschel view of GAS in Protoplanetary Systems (GASPS): First comparisons with a large grid of models

The Herschel GASPS key program is a survey of the gas phase of protoplanetary discs, targeting 240 objects which cover a large range of ages, spectral types, and disc properties. To interpret this large quantity of data and initiate self-consistent analyses of the gas and dust properties of protoplanetary discs, we have combined the capabilities of the radiative transfer code MCFOST with the gas thermal balance and chemistry code ProDiMo to compute a grid of ≈300  000 disc models (DENT). We present a comparison of the first Herschel/GASPS line and continuum data with the predictions from the DENT grid of models. Our objective is to test some of the main trends already identified in the DENT grid, as well as to define better empirical diagnostics to estimate the total gas mass of protoplanetary discs. Photospheric UV radiation appears to be the dominant gas-heating mechanism for Herbig stars, whereas UV excess and/or X-rays emission dominates for T Tauri stars. The DENT grid reveals the complexity in the analysis of far-IR lines and the difficulty to invert these observations into physical quantities. The combination of Herschel line observations with continuum data and/or with rotational lines in the (sub-)millimetre regime, in particular CO lines, is required for a detailed characterisation of the physical and chemical properties of circumstellar discs

Dust and Gas in the disc of HL Tauri: Surface density, dust settling, and dust-to-gas ratio

The recent ALMA observations of the disc surrounding HL Tau reveal a very complex dust spatial distribution. We present a radiative transfer model accounting for the observed gaps and bright rings as well as radial changes of the emissivity index. We find that the dust density is depleted by at least a factor 10 in the main gaps compared to the surrounding rings. Ring masses range from 10-100 M$_{\oplus}$ in dust, and, we find that each of the deepest gaps is consistent with the removal of up to 40 M$_{\oplus}$ of dust. If this material has accumulated into rocky bodies, these would be close to the point of runaway gas accretion. Our model indicates that the outermost ring is depleted in millimetre grains compared to the central rings. This suggests faster grain growth in the central regions and/or radial migration of the larger grains. The morphology of the gaps observed by ALMA - well separated and showing a high degree of contrast with the bright rings over all azimuths - indicates that the millimetre dust disc is geometrically thin (scale height $\approx$ 1 au at 100 au) and that a large amount of settling of large grains has already occurred. Assuming a standard dust settling model, we find that the observations are consistent with a turbulent viscosity coefficient of a few $10^{-4}$. We estimate the gas/dust ratio in this thin layer to be of the order of 5 if the initial ratio is 100. The HCO$^+$ and CO emission is consistent with gas in Keplerian motion around a 1.7 $M_\odot$ star at radii from $\leq 10 - 120\,$au.Comment: 12 pages, 8 figures, published in ApJ, same version as before but with slightly extended discussion on temperature and masses to account for literature published since initial submissio

Uncertainties in water chemistry in disks: An application to TW Hya

Context. This paper discusses the sensitivity of water lines to chemical processes and radiative transfer for the protoplanetary disk around TW Hya. The study focuses on the Herschel spectral range in the context of new line detections with the PACS instrument from the Gas in Protoplanetary Systems project (GASPS). Aims. The paper presents an overview of the chemistry in the main water reservoirs in the disk around TW Hya. It discusses the limitations in the interpretation of observed water line fluxes. Methods. ... (abbreviated) Results. We report new line detections of p-H2O (3_22-2_11) at 89.99 micron and CO J=18-17 at 144.78 micron for the disk around TW Hya. Disk modeling shows that the far-IR fine structure lines ([OI], [CII]) and molecular submm lines are very robust to uncertainties in the chemistry, while the water line fluxes can change by factors of a few. The water lines are optically thick, sub-thermally excited and can couple to the background continuum radiation field. The low-excitation water lines are also sensitive to uncertainties in the collision rates, e.g. with neutral hydrogen. The gas temperature plays an important role for the [OI] fine structure line fluxes, the water line fluxes originating from the inner disk as well as the high excitation CO, CH+ and OH lines. Conclusions. Due to their sensitivity on chemical input data and radiative transfer, water lines have to be used cautiously for understanding details of the disk structure. Water lines covering a wide range of excitation energies provide access to the various gas phase water reservoirs (inside and outside the snow line) in protoplanetary disks and thus provide important information on where gas-phase water is potentially located. Experimental and/or theoretical collision rates for H2O with atomic hydrogen are needed to diminish uncertainties from water line radiative transfer.Comment: accepted for publication in A&

Submillimetre dust polarisation and opacity in the HD163296 protoplanetary ring system

We present ALMA images of the sub-mm continuum polarisation and spectral index of the protoplanetary ringed disk HD163296. The polarisation fraction at 870{\mu}m is measured to be ~0.9% in the central core and generally increases with radius along the disk major axis. It peaks in the gaps between the dust rings, and the largest value (~4%) is found between rings 1 and 2. The polarisation vectors are aligned with the disk minor axis in the central core, but become more azimuthal in the gaps, twisting by up to +/-9degrees in the gap between rings 1 and 2. These general characteristics are consistent with a model of self-scattered radiation in the ringed structure, without requiring an additional dust alignment mechanism. The 870/1300{\mu}m dust spectral index exhibits minima in the centre and the inner rings, suggesting these regions have high optical depths. However, further refinement of the dust or the disk model at higher resolution is needed to reproduce simultaneously the observed degree of polarisation and the low spectral index.Comment: 5 pages +2 pages supplemental data. v2 - revised figures and final values; conclusions unchange

An upper limit on the mass of the circumplanetary disk for DH Tau b

DH Tau is a young ($\sim$1 Myr) classical T Tauri star. It is one of the few young PMS stars known to be associated with a planetary mass companion, DH Tau b, orbiting at large separation and detected by direct imaging. DH Tau b is thought to be accreting based on copious H${\alpha}$ emission and exhibits variable Paschen Beta emission. NOEMA observations at 230 GHz allow us to place constraints on the disk dust mass for both DH Tau b and the primary in a regime where the disks will appear optically thin. We estimate a disk dust mass for the primary, DH Tau A of $17.2\pm1.7\,M_{\oplus}$, which gives a disk-to-star mass ratio of 0.014 (assuming the usual Gas-to-Dust mass ratio of 100 in the disk). We find a conservative disk dust mass upper limit of 0.42$M_{\oplus}$ for DH Tau b, assuming that the disk temperature is dominated by irradiation from DH Tau b itself. Given the environment of the circumplanetary disk, variable illumination from the primary or the equilibrium temperature of the surrounding cloud would lead to even lower disk mass estimates. A MCFOST radiative transfer model including heating of the circumplanetary disk by DH Tau b and DH Tau A suggests that a mass averaged disk temperature of 22 K is more realistic, resulting in a dust disk mass upper limit of 0.09$M_{\oplus}$ for DH Tau b. We place DH Tau b in context with similar objects and discuss the consequences for planet formation models.Comment: accepted for publication in A

Physical characterisation of southern massive star-forming regions using Parkes NH3_3 observations

We have undertaken a Parkes ammonia spectral line study, in the lowest two inversion transitions, of southern massive star formation regions, including young massive candidate protostars, with the aim of characterising the earliest stages of massive star formation. 138 sources from the submillimetre continuum emission studies of Hill et al., were found to have robust (1,1) detections, including two sources with two velocity components, and 102 in the (2,2) transition. We determine the ammonia line properties of the sources: linewidth, flux density, kinetic temperature, NH$_3$ column density and opacity, and revisit our SED modelling procedure to derive the mass for 52 of the sources. By combining the continuum emission information with ammonia observations we substantially constrain the physical properties of the high-mass clumps. There is clear complementarity between ammonia and continuum observations for derivations of physical parameters. The MM-only class, identified in the continuum studies of Hill et al., display smaller sizes, mass and velocity dispersion and/or turbulence than star-forming clumps, suggesting a quiescent prestellar stage and/or the formation of less massive stars.Comment: 20 pages, 9 Figures, 1 appendix (to appear in full online only, a sample appendix in the paper); 7 tables. Accepted by MNRA

The Taurus Boundary of Stellar/Substellar (TBOSS) Survey II. Disk Masses from ALMA Continuum Observations

We report 885$\mu$m ALMA continuum flux densities for 24 Taurus members spanning the stellar/substellar boundary, with spectral types from M4 to M7.75. Of the 24 systems, 22 are detected at levels ranging from 1.0-55.6 mJy. The two non-detections are transition disks, though other transition disks in the sample are detected. Converting ALMA continuum measurements to masses using standard scaling laws and radiative transfer modeling yields dust mass estimates ranging from $\sim$0.3-20M$_{\oplus}$. The dust mass shows a declining trend with central object mass when combined with results from submillimeter surveys of more massive Taurus members. The substellar disks appear as part of a continuous sequence and not a distinct population. Compared to older Upper Sco members with similar masses across the substellar limit, the Taurus disks are brighter and more massive. Both Taurus and Upper Sco populations are consistent with an approximately linear relationship in $M_{dust}$ to $M_{star}$, although derived power-law slopes depend strongly upon choices of stellar evolutionary model and dust temperature relation. The median disk around early M-stars in Taurus contains a comparable amount of mass in small solids as the average amount of heavy elements in Kepler planetary systems on short-period orbits around M-dwarf stars, with an order of magnitude spread in disk dust mass about the median value. Assuming a gas:dust ratio of 100:1, only a small number of low-mass stars and brown dwarfs have a total disk mass amenable to giant planet formation, consistent with the low frequency of giant planets orbiting M-dwarfs.Comment: 41 pages and 32 figures, with all tables and appendices presented here in their entirety. Accepted for publication in AJ (November 26, 2017

Gas and dust in the Beta Pictoris Moving Group as seen by the Herschel Space Observatory

Context. Debris discs are thought to be formed through the collisional grinding of planetesimals, and can be considered as the outcome of planet formation. Understanding the properties of gas and dust in debris discs can help us to comprehend the architecture of extrasolar planetary systems. Herschel Space Observatory far-infrared (IR) photometry and spectroscopy have provided a valuable dataset for the study of debris discs gas and dust composition. This paper is part of a series of papers devoted to the study of Herschel PACS observations of young stellar associations. Aims. This work aims at studying the properties of discs in the Beta Pictoris Moving Group (BPMG) through far-IR PACS observations of dust and gas. Methods. We obtained Herschel-PACS far-IR photometric observations at 70, 100 and 160 microns of 19 BPMG members, together with spectroscopic observations of four of them. Spectroscopic observations were centred at 63.18 microns and 157 microns, aiming to detect [OI] and [CII] emission. We incorporated the new far-IR observations in the SED of BPMG members and fitted modified blackbody models to better characterise the dust content. Results. We have detected far-IR excess emission toward nine BPMG members, including the first detection of an IR excess toward HD 29391.The star HD 172555, shows [OI] emission, while HD 181296, shows [CII] emission, expanding the short list of debris discs with a gas detection. No debris disc in BPMG is detected in both [OI] and [CII]. The discs show dust temperatures in the range 55 to 264 K, with low dust masses (6.6*10^{-5} MEarth to 0.2 MEarth) and radii from blackbody models in the range 3 to 82 AU. All the objects with a gas detection are early spectral type stars with a hot dust component.Comment: 12 pages, 7 figures, 6 table

Mid - infrared interferometry of massive young stellar objects II Evidence for a circumstellar disk surrounding the Kleinmann - Wright object

The formation scenario for massive stars is still under discussion. To further constrain current theories, it is vital to spatially resolve the structures from which material accretes onto massive young stellar objects (MYSOs). Due to the small angular extent of MYSOs, one needs to overcome the limitations of conventional thermal infrared imaging, regarding spatial resolution, in order to get observational access to the inner structure of these objects.We employed mid - infrared interferometry, using the MIDI instrument on the ESO /VLTI, to investigate the Kleinmann - Wright Object, a massive young stellar object previously identified as a Herbig Be star precursor. Dispersed visibility curves in the N- band (8 - 13 {\mu}m) have been obtained at 5 interferometric baselines. We show that the mid - infrared emission region is resolved. A qualitative analysis of the data indicates a non - rotationally symmetric structure, e.g. the projection of an inclined disk. We employed extensive radiative transfer simulations based on spectral energy distribution fitting. Since SED - only fitting usually yields degenerate results, we first employed a statistical analysis of the parameters provided by the radiative transfer models. In addition, we compared the ten best - fitting self - consistent models to the interferometric observations. Our analysis of the Kleinmann - Wright Object suggests the existence of a circumstellar disk of 0.1M\odot at an intermediate inclination of 76\circ, while an additional dusty envelope is not necessary for fitting the data. Furthermore, we demonstrate that the combination of IR interferometry with radiative transfer simulations has the potential to resolve ambiguities arising from the analysis of spectral energy distributions alone.Comment: 12 pages, 22 figures accepted for publication in A&

Characteristics of small protoplanetary disc warps in kinematic observations

Stars and planetary system