Free Access
Issue
A&A
Volume 634, February 2020
Article Number A88
Number of page(s) 12
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201937124
Published online 13 February 2020

© ESO 2020

1 Introduction

The elusive total gas mass of a protoplanetary disc is relevant for planet formation, dust dynamics, and for testing disc evolution models. As a result of difficulties in observing H2, Mgas has been robustly measured in only three cases (Bergin et al. 2013; McClure et al. 2016). In this work, we use Herschel archival data to constrain Mgas in a sample of 15 Herbig Ae/Be discs, and we determine the mass of HD 163296 to within a factor of a few.

The gas mass is dominated by H2, which has a large energy spacing between its lowest rotational levels (para-H2 J = 2–0, ΔE = 512 K) and lacks a dipole moment. As such, H2 is not emissive at the 10–100 K temperatures typical for discs. Dust continuum emission at millimetre wavelengths is often used to estimate Mgas. Gas and dust are linked through a mass ratio, canonically Δg/d = 100 for solar-composition material below ~103 K (e.g. Lodders 2003). While dust emission is easy to detect, the different dust and gas evolution as well as uncertain opacity values limit its reliability in measuring Mgas. The most precise Mgas measurements to date are from hydrogen deuteride (HD) rotational lines. The relative abundance of this deuterated isotopologue of H2 is set by the local absolute atomic ratio, DH  = (2.0 ± 0.1) × 10−5 (Prodanović et al. 2010), and is minimally affected by disc chemistry (Trapman et al. 2017). As the J = 1 rotational level is at EkB = 128.5 K, HD emits from warm gas (Tgas ≈ 30–50 K; Bergin et al. 2013; Trapman et al. 2017). This is sufficient to constrain the total Mgas, especially ifthe temperature structure is constrained via other observables. The HD J = 1–0 line at 112 μm is however impossible to observe from the ground owing to atmospheric absorption and requires air- or space-borne telescopes.

After the pioneering HD 1–0 detection in TW Hya (Bergin et al. 2013; Trapman et al. 2017), facilitated by the PACS spectrometer (Poglitsch et al. 2010)on the Herschel Space Observatory (Pilbratt et al. 2010), further detections were only made in DM Tau and GM Aur (McClure et al. 2016) before the instrument expired. The masses of these T Tauri discs are Mgas = (6−9) × 10−3, (1− 4.7) × 10−2, and (2.5−20.4) × 10−2 M, respectively. An upper limit Mgas ≤ 8× 10−2 M was obtained for the Herbig Ae/Be system HD 100546 (Kama et al. 2016, revised down from the published value because of a mistake in the D/H ratio).

In this work, we use the 2D physical-chemical code DALI (Bruderer et al. 2012; Bruderer 2013) to constrain Mgas in 15 discs by analysing Herschel Space Observatory archival data covering the HD 1–0 and 2–1 lines. The data and models are discussed in Sects. 2 and 3, respectively. In Sect. 4, we explorethe disc mass constraints, with a focus on HD 163296, and discuss the potential for gravitational instability. In Sect. 5, we compare the mass of discs, stars, and planetary systems for stars over 1.4 M. We also discuss future observations of HD with SOFIA/HIRMES (Richards et al. 2018) and SPICA/SAFARI (Nakagawa et al. 2014; Audley et al. 2018).

Table 1

HD line flux upper limits (3σ) for the sample.

2 Observations and sample

We use archival data from the Herschel Space Observatory (Pilbratt et al. 2010) key programme DIGIT (PI: N.J. Evans), which targeted 30 protoplanetary discs with the PACS (Poglitsch et al. 2010) instrument at 50–210 μm. Detected gaseous species in this data are presented in Fedele et al. (2013) and Meeus et al. (2013). We analyse upper limits on HD J = 1–0 and 2–1 lines at 112 and 56 μm for the 15 Herbig Ae/Be discs in the sample. Because of the intrinsically higher luminosity of their host stars (~ 10 –100 L), these discs are warmer and brighter in continuum and line emission than those around T Tauri stars. This enables tighter constraints for discs at equivalent distance.

We selected discs around stars of spectral type mid-F to late-B, including well-known targets such as HD 100546 and HD 163296. HD 50138 was excluded as it is likely an evolved star (Ellerbroek et al. 2015), and HD 35187 was excluded because it is a binary of two intermediate-mass stars and not directly comparable to our model grid. The data are spectrally unresolved, with δv ≈ 100 km s−1 (λδλ = 3000) at the shortest wavelengths (51 μm), while expected linewidths are ≤ 10 km s−1. Exposure times ranged from 4356 to 8884 s. The system parameters and 3σ line flux upper limits are given in Table 1.

We obtained flux limits for the HD transitions from the 1σ noise reported for the nearest lines of other molecules from Fedele et al. (2013): OH 2Π1∕2 J = 9∕2–7∕2+ at 55.89 μm for the 56 μm line and OH 2Π3∕2 J = 5∕2–3∕2+ at 119.23 μm for the 112 μm line. With a typical 1σ uncertainty of 5 × 10−18 W m−2 at 112 μm and 2 × 10−17W m−2 at 56 μm, neither of the HD lines is detected in the targets, individually or stacked. For comparison, the HD 1–0 detections in Bergin et al. (2013) and McClure et al. (2016) have respective uncertainties of roughly 7 × 10−19 W m−2 and 5 × 10−19 W m−2, which illustrates the difference between those targeted, deep integrations and the survey-type observations analysed in this work.

The discs fall into two categories: cold (flat, group II in the Meeus classification, Meeus et al. 2001) and warm (flaring, group I). This characterises the shape of the radial optically thick surface, where starlight is effectively absorbed. Starlight impinges at a shallow angle on flat discs, and heating is inefficient compared to that above the same mid-plane location in a flaring disc. In addition, among the Herbig Ae/Be systems flaring, group I discs have resolved cavities or gaps 10–100 au scales in their millimetre dust emission (Maaskant et al. 2013; Kama et al. 2015).

3 Modelling

3.1 DALI

To determine the behaviour of the HD 1–0 line and 1.3 mm continuum flux as a function of disc structure parameters, we ran a grid of models with the 2D physical-chemical disc code DALI (Bruderer et al. 2012; Bruderer 2013). The surface density is parametrised following the viscous accretion disc formalism (Lynden-Bell & Pringle 1974; Hartmann et al. 1998): (1)

where Σc is the surface density at the characteristic radius Rc and γ the power-law index, which is generally 1. Assuming an isothermal structure in hydrostatic equilibrium, the vertical structure is given by a Gaussian density distribution (Kenyon & Hartmann 1987): (2)

In this equation , ψ is the flaring index, and hc is the disc opening angle at Rc.

A population of small grains (0.005–1 μm), which have a mass fraction fsmall, follows the gas density distribution given in Eq. (2). A second population, consisting of large grains (1 μm–1 mm), has a mass fraction flarge. The scale height of this second population is χh, where χ ∈ (0, 1) is the settling parameter.

For the dust opacities of both small and large grain populations we assume a standard interstellar composition following Weingartner & Draine (2001), in line with Bruderer (2013). The absorption coefficient for the small (large) grains is 29.9 cm2 g−1 (30.0 cm2 g−1) at 112 μm and 154 cm2 g−1 (46.3 cm2 g−1) at 56 μm.

First, the radiation field and dust temperature are determined from Monte Carlo radiative transfer. Next, the gas temperature (heating-cooling balance) and chemical composition (steady-state) are solved for iteratively. Raytracing then yields simulated line and continuum observations.

HD chemical network versus fixed abundance

The HD abundance (HD/H2) can be prescribed as a constant or obtained from solving a chemical reaction network. In the parametric approach, the HD abundance is determined by the local D/H ratio, which for the local interstellar medium (ISM; within ≈ 2 kpc) is measured to be (DH)ISM = (2.0 ± 0.1) × 10−5 (Prodanović et al. 2010). Assuming all deuterium is in HD, this gives HD/H2 = 4 × 10−5.

A more refined approach is to calculate the HD abundance using a reaction network which includes deuterium. Trapman et al. (2017) extends the standard DALI chemical network (originally based on the UMIST06 database Woodall et al. 2007) to include the species HD, D, HD+, and D+. These authors include HD formation on dust and ion-exchange reactions in addition to HD self-shielding. The details of the implementation are described in Sect. 2.3 of Trapman et al. (2017).

Using the chemical network approach, we find that all of the available deuterium is locked up in HD for the vast majority of the disc, and the parametric abundance of HD/H2 = 4 × 10−5 is appropriate to use. The network produces less HD in only two regions: the uppermost layers of the disc where HD is photodissociated and in a thin intermediate layer, where the HD abundance is decreased by a factor of ~ 2. Tests determined that neither of these significantly affects the disc-integrated HD line flux. Given the very close match between the two approaches, we opt for simplicity and fix the HD/H2 ratio at 4 × 10−5.

Part of our analysis below involves modelling CO rotational lines. Owing to processes such as chemical conversion and freeze-out, the gas-phase total abundance of C and O nuclei can be more than a factor of ten below nominal (e.g. Favre et al. 2013; Kama et al. 2016), which makes CO-based mass estimates highly uncertain. We refer to the reduction of gas-phase C and O nuclei below their total values with the term depletion, and the phenomenon can be included in our modelling as a reduction of the total amount of volatile C or O input into a given DALI model. This is relevant for Sect. 4.3, in particular.

Table 2

DALI model grid parameters.

3.2 Model grid

To investigate the range of disc properties constrained by the Herschel upper limits on the HD 1–0 line, we ran a grid of Herbig Ae/Be disc models covering a wide range of parameters, summarised in Table 2. The disc gas masses are Mgas = 10−3, 10−2, and 10−1 M. Dust mass is defined by the gas-to-dust mass ratio that has values Δg/d = 10, 50, 100, and 300 and ranges from Mdust = 3 × 10−6 to 10−2 M. The shape of the stellar spectrum, including UV excess, is based on HD 100564 from Bruderer et al. (2012). The spectrum is scaled to the total stellar luminosity, L ∈ [10, 50, 115] L. This covers the sources in our sample, as given in Table 1. In total we ran 2304 models; parameters are given in Table 2. Our fiducial model has hc = 0.15, Rc = 50 au, Δgd = 100, flarge = 0.95, χ = 0.2, and L = 10 L.

Figure 1 shows the HD 1–0 emitting regions and disc mass outline for models representing extremes in flaring (Ψ = 0.0 and hc = 0.05 for flat, and Ψ = 0.3 and hc = 0.15 for flared), radial extent (Rc = 50 and 125 au), and total disc mass. From the figure it is clear that the flared disc (ψ = 0.3, hc = 0.15), shown in red, has a much large emitting region than the flat disc (Ψ = 0.0, hc = 0.05), shown in blue. In both cases the HD 1–0 emission originates from the warm layer above the midplane.

thumbnail Fig. 1

HD1–0 line emitting regions in our flat/cold (blue) and flared/warm (red) disc models. Solid contours contain the middle 75% of cumulative line emission. Dashed lines indicate gas number density iso-contours for ngas = 106 cm−3, acting as a disc outline.

4 Results

In Fig. 2, we show the HD J = 1–0 flux as a function of Mgas and 1.3 mm continuum flux. The warm, flaring, group I discs and cold, flat, group II discs are highlighted separately for clarity.

4.1 Parameter dependences in the grid

Dependences of the HD 1–0 line and 1.3 mm continuum flux on the main model parameters are shown in Fig. 3. The HD line flux depends linearly on Mgas, which has only a marginal effect on the dust emission. For a fixed Mgas, a 1 dex increase in Mdust leads to a factor 6.7 lower HD and 2.5 higher continuum flux. The flaring structure of the disc has the largest influence, as the HD line flux increases by a factor of 26 when the flaring parameter Ψ goes from 0 (height is linear with radius, inefficient heating) to 0.3 (very flared and efficiently heated). The Meeus group corresponds to the flaring structure (group I discs are flared, II flat).

A near-linear dependence of HD line flux on Mgas arises because the HD line emission in the models is vertically limited by the dust optical depth τ at 112 μm out to ≈ 100 au radii, beyond which the surface density drops rapidly. Thus the HD contribution from the gas above and radially outside the dust scales linearly with the total gas mass. Dust emission, to first order, is optically thin at 1.3 mm, and thus scales linearly with the total dust mass. Again because the dust optical depth dominates at the 112 μm wavelength of HD 1–0, increasing the dust mass in a given column lifts the vertical τ(112 μm) = 1 surface, hiding a larger fraction of the HD molecules.

4.2 Constraints on Mgas across the sample

A comparison of the HD upper limits from Herschel with our DALI model grid (Figs. 2 and 4) places an upper limit of approximately Mgas ≤ 0.1 M for the discs in our sample. Among the flared, group I discs (Fig. 2, upper row), we find Mgas < 0.02–0.03 M for IRS 48, HD 36112, HD 100453, and HD 135344B, while among the flat, group II discs HD 163296 has a limit at < 0.1 M.

Source-specific models can tighten the mass limit for individual discs. We ran a small grid of models for HD 163296, in which we have a strong HD upper limit and a wide comparison range of indirect gas mass estimates from the literature based on various isotopologues of CO.

Table 3

Adopted model for HD 163296.

4.3 HD 163296

We constrain the gas mass in the HD 163296 disc to Mgas ≤ 0.067 M (Fig. 5). Given that the disc-integrated dust mass in our model is 6.7 × 10−4 M, this constrains the gas-to-dust ratio to Δg/d ≤ 100 and has implications for the gas-phase volatile abundances, which we discuss below. This source-specific model matches the continuum spectral energy distribution,12CO rotational ladder and isotopologue lines, and several other key volatile species. The full details of this modelling are outside the scope of this paper and will be published separately; below we focus on the main outcomes of the continuum, CO, and HD modelling.

HD 163296 is one of the largest known discs with a CO J = 3–2 gas emission radius of 540 au (Rosenfeld et al. 2013). Fitting of CO and 850 μm continuum emission, observed by ALMA, with a tapered surface density power law yielded γ = 0.9 and Rc = 125 au (Tilling et al. 2012; de Gregorio-Monsalvo et al. 2013). We model HD 163296 with the stellar spectrum from the ProDiMo project (Woitke et al. 2019), fixing the shape of the dust surface density profile to the above parameters and varying the gas mass. To satisfy the radial profile of CO 3–2 emission simultenaously with the spectral energy distribution, we find the density profile flaring index in Eq. (2) is around ψ =0.05, consistent with the range 0.019–0.066 found by Tilling et al. (2012). The morpology of the 12 CO 3–2 channel maps, in which both the near and far side of the disc can be seen, suggest HD 163296 is more flared (ψ ≈ 0.12, de Gregorio-Monsalvo et al. 2013)than our model (ψ = 0.05). However, these two ψ-s differ in physical meaning: the CO-based ψ measures the observed shape of the CO-emitting surface, while the disc structure parameter ψ characterises the shape of the total gas mass distribution (see Eq. (2)).

Our model, which hits the HD upper limit, reproduces the observed dust emission across the far-infrared and sub-millimetre wavelengths as well as various spatially resolved and unresolved emission lines of 12 CO and its isotopologues, and has a gas-to-dust ratio Δg/d = 100.

Most previous estimates of the HD 163296 gas mass relied on low-J emission lines of CO isotopologues, and used a range of modelling approaches from generic model grids to tailored modelling with physical-chemical codes. Those Mgas estimates range from 8 × 10−3 to 5.8 × 10−1 M (Isella et al. 2007; Williams & Best 2014; Boneberg et al. 2016; Miotello et al. 2016; Williams & McPartland 2016; Powell et al. 2019; Woitke et al. 2019; Booth et al. 2019). The mass obtained from the most optically thin isotopologue among these, 13 C17O, was 2.1 × 10−1M (Booth et al. 2019).

Above, we assumed an undepleted solar abundance for elemental gas-phase carbon and oxygen. Our model matching the HD upper limit overproduces the low-J line fluxes of CO isotopologues by a factor of a few. Since the rarer isotopologues are progressively more optically thin, we can reproduce their line fluxes by decreasing the gas-phase elemental carbon and oxygen abundance proportionately to the flux mismatch. Since the millimetre-wave dust emission and HD upper limit constrain the gas-to-dust mass ratio to be ≤ 100, we can combine the above considerations to arrive at the following three distinct hypotheses for HD 163296.

  • 1.

    Mgas is just sufficiently below our upper limit of 6.7 × 10−2 M for HD not to be detected. If so, then as the dust mass is fixed, it follows from our models that Δg/d = 100 and that total gas-phase elemental C and O are depleted by up to a factor of a few.

  • 2.

    Mgas is a factor of a few below our HD limit, and the total gas-phase elemental C and O abundances are not depleted with respect to their interstellar values. If so, the implication is that Δg/d ≈ 20–50. This relative depletion of gas over dust is supported by the hydrostatic MCMax modelling of the SED and low-J 12CO, 13CO, and C18O lines by Boneberg et al. (2016), whose best models had 9.2 < Δg/d < 18. It is also consistent with the inner disc value Δg/d ≈ 55, measured using accretion onto the central star by Kama et al. (2015, their Fig. 2).

  • 3.

    Mgas is far below our upper limit. In this hypothesis, the total C and O abundance in the gas must be enhanced above the interstellar baseline to still match the optically thin CO isotopologues. This would be the first known case of C and O enhancement, however the inner disc composition analysis by Kama et al. (2015) does not show evidence for a strong enhancement of gas-phase volatile elements over total hydrogen.

Thus Δg/d > 100 is ruled out by the HD 1–0 upper limit for HD 163296, independently of assumptions about the precise abundance of gas-phase volatiles.

The abundance of volatile elements in the HD 163296 disc may be depleted or enhanced by up to a factor of a few, depending on the true value of Mgas and on the somewhat uncertain underlying number abundance ratios of 12CO and its various isotopologues. We note that even with the flat, cold disc structure of HD 163296, our Δg/d = 100 model somewhat over-predicts the CO emission outside of ~ 100 au for an undepleted elemental carbon abundance (C/H = 1.35 × 10−4). A more flared surface would aggravate this over-prediction, while globally reducing the elemental C under-predicts the CO 3–2 inside ~100 au. This may indicate that any depletion of gas-phase volatile elemental C and O, reflected in the CO abundance in the warm molecular layer, is restricted to the region beyond the CO snowline, which has been observed to be at ≈ 90 au (Qi et al. 2015). The same conclusion was recently reached by Zhang et al. (2019) through an analysis of spatially resolved C18 O data, which yielded a factor of ten depletion of gas-phase CO outside the CO snowline.

thumbnail Fig. 2

Distance-normalised 3σ upper limits on HD 112 μm line flux forthe disc sample (black lines and circles) compared with our grid of DALI disc models (coloured crosses). Highlighted crosses show the HD 112 μm line flux ofour fiducial model. The top panels show the group I sources compared to models with flaring angle ψ = 0.3. The bottom panels show the group II sources compared to models with ψ = 0.0. Left: models are separated based on gas mass. Right: HD 1–0 upper limitsset against 1.3 mm continuum fluxes for both observations and models.

thumbnail Fig. 3

HD 1–0 line and dust continuum flux dependences on disc and stellar parameters.

thumbnail Fig. 4

HD 1–0 line flux vs. stellar luminosity. Observed stellar luminosities taken from Table 1. Model stellar luminosities are given a small offset for clarity. Highlighted crosses show our fiducial model (hc = 0.15, Rc = 50 au, Δgd = 100, flarge = 0.95, χ = 0.2).

thumbnail Fig. 5

Comparing the HD 163296 specific models to the HD 1–0 upper limit (Fedele et al. 2013). All models have a dust mass Mdust = 6.6 × 10−4 M (Table 3). The red bar shows the range of gas masses inferred from CO in the literature.

4.4 HD 100546

HD 100546 was previously modelled with DALI by Bruderer et al. (2012) who determined the radial and vertical structure of the disc mainly from CO lines and continuum emission. A refined version of this modelling effort included the Herschel HD upper limits, the C0 and C2 H fluxes, and the spatially resolved CO 3–2 emission, constraining the gas mass to 8.1 × 10−3Mgas ≤ 2.4 × 10−1 M (Kama et al. 2016). The highest mass model had Δg/d = 300, with a dust mass anchored by the continuum spectral energy distribution. Owing to a factor of four error in the D abundance used in that model, we revise those numbers to ≲ 100 and thus Mgas ≲ 0.08 M from the Kama et al. (2016) model. This is about a factor of two stronger than the constraint from our general model grid, so in Fig. 6 we adopt Mgas ≲ 0.08 M.

4.5 Other individual discs

HD 97048

Hosts a massive dust disc, Mdust ≃ 6.7 × 10−4 M (Walsh et al. 2016), so it is likely the gas mass is also high. The disc surface is highly flared (Ψ = 0.5−0.73; see e.g. Lagage et al. 2006; Walsh et al. 2016; Ginski et al. 2016; van der Plas et al. 2019). This exceeds the largest Ψ in our general grid, but we note again that the CO-surface Ψ and the density structure Ψ differ in physical meaning. From our grid we find Mgas ≤ 9.4 × 10−2 Mg/d ≤ 200).

HD 104237

For this disc, Hales et al. (2014) determined Mdust = 4× 10−4 M, which assuming Δg/d = 100 implies a total mass Mgas = 4× 10−2 M. This is consistent with our upper limit from HD 1–0, which yields an upper limit of Δg/d ≤ 300 (Fig. 2).

HD 36112 (MWC 758)

Based on millimetre continuum interferometry, Guilloteau et al. (2011) inferred a disc mass of (1.1 ± 0.2) × 10−2 M. Our analysis of the 1.3 mm continuum flux and the HD 1 – 0 upper limit matches both data points for Δg/d ≈ 100 and a disc mass of order 10−2 M. A substantially lower gas mass would imply a very low Δg/d mass ratio.

thumbnail Fig. 6

Mass of selected discs and planets around B9 to F5 type stars. Vertical lines show the cumulative mass of each planetary system; dots highlight planets from the most massive at bottom. Disc gas mass upper limits from HD lines are taken from this work (HD 163296) and from Kama et al. (2016, HD 100546). For HD 163296, the range of CO isotopologue-based disc mass estimates is shown by a light blue bar (8 × 10−3 to 5.8 × 10−1 M; references in text). Also shown are the stellar mass divided by 10 and age; the mass limit for a gravitationally unstable disc (dashed line); an extrapolated dust-based disc mass range (dotted lines; Pascucci et al. 2016); and a population density colour map for planets around B9 to F5 type stars (data retrieved from exoplanets.org on 2019.07.16; bins contain from bottom to top 7, 6, and 1 planet). See text for individual planet and stellar mass references.

HD 31293 (AB Aurigae)

From 1.3 mm continuum observations performed using the SMA, Andrews et al. (2013) inferred a dust mass of (1.56 ± 0.09) × 10−4 M, implying Mgas = 1.56 × 10−2 M assuming Δg/d = 100. The high upper limit of HD 1–0 for this source does not allow us to put any meaningful constraints on the gas mass based on HD.

HD 135344B

It has been modelled by van der Marel et al. (2016) to determine the physical structure. Using ALMA observations of 13 CO J = 3–2, C18 O J = 3–2, 12 CO J = 6–5 and dust 690 GHz continuum, these authors determined a gas mass Mgas = 1.5 × 10−2 M. We ran models based on their physical structure and find the resulting HD 1–0 flux to be inagreement with the upper limit (see Fig. A.1).

HD 142527

Modelling interferometric 880 μm continuum and 13 CO 3–2 and C18 O 3–2 line observations, Boehler et al. (2017) determine a dust mass of 1.5 × 10−3 M and a gas mass of 5.7 × 10−3 M (see also Muto et al. 2015). This gives 3 ≤ Δg/d ≤ 5 and suggests the gas is either strongly depleted in elemental C and O, or dissipating entirely. Because of the loose HD 1–0 upper limit for this source, we cannot provide an independent check of the low Δg/d derived from CO.

HD 179218

From the integrated 1.3 mm flux Mannings & Sargent (2000) inferred a dust mass of (1.5 ± 0.15) × 10−4 M, implying Mgas = 1.5 × 10−2 M assuming Δg/d = 100. Again the HD 1–0 upper limit provides no meaningful constraint on the gas mass.

HD 100453

Based on millimetre continuum interfermotric observations, van der Plas et al. (2019) infer a dust mass of 6.7 × 10−5 M. By comparing the 13CO 2–1 and C18O 2–1 to the disc model grid in Williams & Best (2014), they determine a gas mass of (1 − 3×) × 10−3 M. Combining both disc masses implies a gas-to-dust mass ratio of Δg/d 15−45. From our analysis of the 1.3 mm continuum flux and the HD 1–0 upper limit we constrain gas mass to Mgas ≤ 10−2 M and the gas-to-dust mass ratio Δg/d ≤ 300. Both constraints are in agreement with the results of van der Plas et al. (2019).

HD 169142

From interferometric 1.3 mm continuum and 12CO 2–1, 13CO 2–1, and C18 O 2–1 line observations, Panić et al. (2008) derived a dust mass of 2.16 × 10−4 M and a gas mass of (0.6−3.0) × 10−2 M. Fedele et al. (2017) find similar disc masses based on higher resolution observations. Constraints based on our analysis of the 1.3 mm continuum flux and the HD 1–0 upper limit put the gas mass at Mgas ≤ 4 × 10−2 M and Δg/d ≤ 300, both of which are in good agreement with previous results.

Oph IRS48 (WLY 2-48)

van der Marel et al. (2016) modelled the resolved 440 μm continuum and 13CO 6–5 and C18O 6–5 line observations. These authors derive a dust mass of 1.5 × 10−5 M and a gas mass of 5.5 × 10−4 M, giving a gas-to-dust mass ratio of Δg/d ≈ 37. Constraints from our analysis of the HD 1–0 line flux and 1.3 mm continuum give Mgas ≾ 10−2 M and Δg/d ≾ 300. These upper limits agree with previous results.

4.6 The gravitational stability of the discs

Constraints on Mgas allow us to test whether the discs in our sample are currently gravitationally stable. Gravitational instability, leading to spirals or fragmentation, occurs in disc regions which are dense and cold and have low orbital shearing on the timescale of the instability (i.e. at large radii). This is quantified with the Toomre Q parameter, (Toomre 1964), which simplifies to (3)

following Kimura & Tsuribe (2012). If Q < 1, the disc fragments. For 1 < Q < 2, the disc is marginally stable, developing transient spirals and clumps, while for Q > 2 it is stable against gravitational collapse. Assuming a surface density profile and MdiscMgas, we obtain (4)

Our most massive disc models have Mgas = 0.1 M. Taking a characteristic radius r0 = 100 au, we find Q = 1.5, which is marginally stable. The discs for which we have the weakest upper limits relative to the massive disc models – HD 142527, HD 144668, HD 179218, and HD 31293 – may potentially be gravitationally unstable within the limits of the Herschel HD data. For the rest of the sample, a gravitationally unstable Mgas is effectively ruled out, i.e. they are most likely stable.

Dust dips, rings, or cavities may locally affect the temperature structure of the gas and thus, through the sound speed, the local Q in a disc (QTkin0.5, i.e. a weak dependence). In general, a lower dust surface density leads to more efficient external heating and thus more stability. Inside a local dust enhancement, the temperature may drop somewhat, but if the region is already very optically thick, the effect on Tkin is minor. We therefore have not considered the effect of such Σdust perturbations in this paper.

5 Discussion

5.1 Mass of discs, stars, and planets

Intermediate-mass stars (spectral types B9 to F5, masses 1.5–3 M) host some of the best-studied protoplanetary discs and high-mass planetary systems. Several Herbig Ae/Be protoplanetary discs have also yielded detections of protoplanet candidates. This presents an opportunity to investigate equivalent planetary systems at different stages of evolution.

In Fig. 6, we compare the disc gas mass with the host star and the candidate protoplanets in the disc. We show two Herbig Ae/Be systems with strong mass limits, HD 163296 (Mgas ≤ 0.067 M; this work) and HD 100546 (Mgas ≤ 0.08 M; Kama et al. 2016). Much of the work on embedded protoplanet candidates quoted below is very new. There are large, at least a factor of two to ten, uncertainties behind the planet mass estimates below, in particular for those inferred from dust gaps where the α viscosity parameter plays a role. We adopt middle-ground values from the literature to begin a discussion on comparing disc and embedded planet masses.

For HD 163296, our HD-based upper limit rules out a large fraction of the wide range of CO isotopologue based Mgas estimates from the literature. Of those still possible, the lowest is Mgas = 8 × 10−3 M. The presence of five giant planets has been inferred from dust gaps and gas kinematics: at 10 au with a mass (0.53 ± 0.18) MJup for αvisc = 10−4 to 10−3 (Zhang et al. 2018); at 48 au with 0.46 MJup (Isella et al. 2016; Liu et al. 2018); at 86 au with (1 ± 0.5) MJup (Liu et al. 2018; Teague et al. 2018); at 145 au with 1.3 MJup (Liu et al. 2018; Teague et al. 2018); and at 260 au with 2 MJup (Pinte et al. 2018). Using the HD- and CO-based Mgas limits, and taking the combined mass of all published protoplanets in this disc as ≈ 5 MJup, we find the HD 163296 disc has converted 10–40% of its mass into giant planets.

For HD100546, the planet masses were constrained to be ≈ 10 MJup at 10 au and ~ 10 MJup at 70 au by Pinilla et al. (2015). The mass of the outer planet could be < 5 MJup (>15 MJup) if it formed very early (late), so we adopt 10 MJup. The HD-based Mgas upper limit and the combined mass of the candidate planets yield a lower limit on the disc-to-planet mass conversion efficiency of ≳ 30%.

Such high disc-to-planet mass conversion efficiencies combined with the presence of several gas giants per star raise the question of whether the planets formed through gravitational instability. Adding the Mgas upper limit and combined mass of proposed planets in either disc gives a result close to 0.1 M. This is approximately at the gravitationally unstable limit, so such a formation pathway may be feasible even with the current total mass in the system, although the local Toomre Q varies with radius and may leave the outer disc still far from instability (e.g. Booth et al. 2019).

We also show in Fig. 6 three somewhat older stars of similar mass (HD 95086, β Pic, and HR 8799) and their planets; standard disc mass estimates for stars of 1.5 and 3 M based on Mdust relations from Pascucci et al. (2016) and scaled up with Δg/d = 100; and a shaded log-scale histogram of the mass distribution of known planets around early-type stars1. Stellar masses are from the Gaia DR2 analysis by Vioque et al. (2018), and from David & Hillenbrand (2015, β Pic) and Stassun et al. (2018, HD 95086). Planet masses for individual systems are plotted as cumulative bars with the highest mass planet at the base. We compiled planet data from Teague et al. (2018), Pinte et al. (2018, 2019), Pinilla et al. (2015), Liu et al. (2018), Zhang et al. (2018), Rameau et al. (2013a,b), De Rosa et al. (2016), and Marois et al. (2008, 2010). Individual stellar masses are from Rhee et al. (2007), David & Hillenbrand (2015), Stassun et al. (2018), and Vioque et al. (2018).

The two HD-based disc Mgas limits in Fig. 6 exceed the combined mass of planets around HR 8799, the most massive known planetary system, by a factor of only three. The disc mass limits are also only a factor of three above combined mass of candidate protoplanets in the HD 100546 disc. Either A-type star discs can, in some cases, convert 10% or more of their mass into giant planets, or these planetary systems formed at a very early stage, perhaps while the central protostar and massive initial disc were still heavily accreting from the protostellar envelope in which they were embedded. The mass distribution of giant planets around main-sequence A and B stars (Fig. 6) is strongly skewed towards lower masses, suggesting that such extreme mass conversion events are either rare, or that the high-mass planetary systems are not stable on timescales beyond a few times 10 Myr.

5.2 Observing HD in Herbig discs with SOFIA/HIRMES, SPICA/SAFARI and Origins Space Telescope

In the coming years, several facilities will or may become available for observing HD rotational lines. The HIRMES instrument for SOFIA is currently undergoing commissioning and is due to be delivered at the end of 2020 (Richards et al. 2018). HIRMES will have a high spectral resolution ofR ~ 100 000, allowing us for the first time to spectrally resolve the HD 1–0 line. The sensitivity of HIRMES will be similar to Herschel/PACS. Our models suggest some Herbig Ae/Be discs will be detectable with this instrument, assuming the necessary hours per source are available.

Figure 7 shows the detectability of our disc models with a 10 h SOFIA/HIRMES observation, assuming a distance of 150 pc. Of the flat models (group II discs), only the most massive (Mgas ~ 0.1 M) around stars with the highest stellar luminosity (L* ≥ 50 L) are detectable. Among the flared models (group I), a larger fraction of discs is observable. All of the disc models Mgas = 0.1 M where Δgd > 10 should be detectable in 10 h with SOFIA/HIRMES. For those discs with Mgas = 0.01 M, all systems with L* = 125 L and most systems with L* = 50 L are detectable. To maximise the chance of success, future SOFIA/HIRMES observations should select group I sources with high stellar luminosity.

Based on the stellar luminosities in Table 1 there are four group I sources that match these criteria best for SOFIA/HIRMES to detect the HD 1–0 line: HD 31293 (AB Aur), HD 100546, HD 179218 and HD 97048. For these sources a 10 h observation with SOFIA/HIRMES would improve the current upper limits by a factor 3–10 and constrain the gas-to-dust mass ratio to Δg/d ≤ 50−100 if the sources remain undetected.

Beyond SOFIA/HIRMES there are two proposed space missions focussing on far-infrared observations: SPICA/SAFARI and Origins Space Telescope. The SPICA mission is one of the competitors for ESA’s M5 opportunity, with a resolving power R ~ 3000 and a 5σ 1 h sensitivity of 1.3 × 10−19 W m−2 at 112 μm (Audley et al. 2018). The Origin Space Telescope is a NASA mission concept. It would have high spectral resolution (R ~ 43 000) and sensitivity (~ 1 × 10−20 W m−2 in 1 h) at 112 μm (Bonato et al. 2019). Hydrogen deuteride in all Herbig Ae/Be discs, and many T Tauri discs, within ~ 200 pc will be detectable with these missions. However, both still require final approval and would only become available at the end of the 2020 s at the earliest. If approved, these missions would be an enormous step forward in planet-forming disc studies.

thumbnail Fig. 7

Observability of Group Ia,Ib (left) and Group IIa (right) models with SOFIA HIRMES. Coloured disc models are detectable (≥5σ) with a 10 h integration. Dark red dashed line shows the SPICA/SAFARI 1 h detection limit. The Origins Space Telescope 1 h detection limit (~1 × 10−20 W m−2) lies below the limits of the figure. The fluxes are calculated for a distance of 150 pc.

6 Conclusions

We have studied the significance, in terms of total gas mass, of Herschel Space Observatory upper limits on HD 1–0 line emission for 15 individual protoplanetary discs. We find an overall gas mass upper limit of Mgas ≤ 0.1 M for most of the discs studied. None of the discs are very likely to be strongly gravitationally unstable, although the constraints for HD 142527, HD 144668, HD 179218, and HD 31293 (AB Aur) are weak enough to allow for the possibility.

The HD 163296 disc mass is Mgas ≤ 6.7 × 10−2 M, based on the HD 1–0 upper limit. The CO-based literature lower limit is Mgas = 8× 10−3 M, which is contingent on the true level of gas-phase volatile depletion. The gas-to-dust ratio is thus 12 ≤ Δg/d ≤ 100, indicating gas dissipation may be proceeding faster than dust removal in this disc. This is consistent with Δg/d=55 inferred from the accretion-contaminated photosphere of the central star (Kama et al. 2015).

Comparing the HD 163296 and HD 100546 Mgas constraints with their protoplanet candidates and the HR 8799 giant planet system, we find that at least some Herbig Ae/Be discs convert the equivalent of 10–40% of their present-day mass into giant planets.

Near-future SOFIA/HIRMES observations will probe the mass of flaring discs and large flat discs around A-type stars within ≈ 150 pc with ≳ 10 h integrations. The SPICA/SAFARI mission will be crucial for larger sample studies of Mgas in discs. The Origins Space Telescope, if approved, would further revolutionise the field.

Acknowledgements

We thank Sebastiaan Krijt for useful discussions. M.K. gratefully acknowledges funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie Fellowship grant agreement No 753799. L.T. is supported by NWO grant 614.001.352. D.F. acknowledges financial support provided by the Italian Ministry of Education, Universities and Research, project SIR (RBSI14ZRHR). A.M. and C.C. gratefully acknowledge funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 823823, (RISE DUSTBUSTERS), and A.M. funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)– Ref no. FOR2634/1 ER685/11-1. E.A.B. gratefully acknowledges support from NASA via grant NNX16AB48G-XRP. All figures were generated with the PYTHON-based package MATPLOTLIB (Hunter 2007).

Appendix A HD 1–0 fluxes for HD 135344B

Based on the HD 135344B source-specific model from van der Marel et al. (2015, 2016), we ran a series of ten models, varying the disc gas mass between 3.75 × 10−3 M and 3 × 10−2 M. Figure A.1 compares the HD 1–0 line fluxes of these models to the observed upper limit (Table 1). From the CO isotopologue observations van der Marel et al. (2016) infer Mgas = 1.5 × 10−2 M. This gas mass is in agreement with the gas mass upper limit inferred from HD 1–0, Mgas ≤ 2.3 × 10−2 M. We note that both gas masses are much lower than 0.1 M, making it highly unlikely that HD 135344B is gravitationally unstable (Sect. 4.6).

thumbnail Fig. A.1

Comparing the HD 135344B specific models from van der Marel et al. (2016) to HD 1–0 upper limit (Fedele et al. 2013). All models have a dust mass Mdust = 1.3 × 10−4 M (cf. Table 3 in van der Marel et al. 2016). The red circle shows the gas mass inferred from CO by van der Marel et al. (2016).

Appendix B HD 2–1 upper limits versus the model fluxes

thumbnail Fig. B.1

Upper limits on HD 56 μm line flux for the sample of Herbig Ae/Be disc systems (black lines) compared with our grid of DALI disc models (crosses). The top panels show the group I sources compared to models with flaring angle ψ = 0.3. The bottom panels show the group II sources compared to models with ψ = 0.0. Left: modelsare separated based on gas mass. Right: HD 2–1 upper limits set against 1.3 mm continuum fluxes for both observationsand models.

Appendix C HD 1–0 line versus 1.3 mm continuum fluxes, showing gas-to-dust ratios and stellar luminosities

thumbnail Fig. C.1

HD 1–0 line flux vs. 1.3 mm continuum fluxes for both observations and models. The panels are similar to right panels of Fig. 2, but also show the model gas-to-dust mass ratios (marker shape) and stellar luminosities (marker size).

References

  1. Andrews, S. M., Rosenfeld, K. A., Kraus, A. L., & Wilner, D. J. 2013, ApJ, 771, 129 [NASA ADS] [CrossRef] [Google Scholar]
  2. Audley, M. D., de Lange, G., Gao, J.-R., et al. 2018, SPIE Conf. Ser., 10708, 107080K [Google Scholar]
  3. Bergin, E. A., Cleeves, L. I., Gorti, U., et al. 2013, Nature, 493, 644 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  4. Boehler, Y., Weaver, E., Isella, A., et al. 2017, ApJ, 840, 60 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bonato, M., De Zotti, G., Leisawitz, D., et al. 2019, PASA, 36, e017 [NASA ADS] [CrossRef] [Google Scholar]
  6. Boneberg, D. M., Panić, O., Haworth, T. J., Clarke, C. J., & Min, M. 2016, MNRAS, 461, 385 [NASA ADS] [CrossRef] [Google Scholar]
  7. Booth, A. S., Walsh, C., Ilee, J. D., et al. 2019, ApJL, 882, L31 [NASA ADS] [CrossRef] [Google Scholar]
  8. Brown, A., Vallenari, A., Prusti, T., et al. 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Bruderer, S. 2013, A&A, 559, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bruderer, S., van Dishoeck, E. F., Doty, S. D., & Herczeg, G. J. 2012, A&A, 541, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. David, T. J., & Hillenbrand, L. A. 2015, ApJ, 804, 146 [NASA ADS] [CrossRef] [Google Scholar]
  12. de Gregorio-Monsalvo, I., Ménard, F., Dent, W., et al. 2013, A&A, 557, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. De Rosa, R. J., Rameau, J., Patience, J., et al. 2016, ApJ, 824, 121 [NASA ADS] [CrossRef] [Google Scholar]
  14. Ellerbroek, L. E., Benisty, M., Kraus, S., et al. 2015, A&A, 573, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Fairlamb, J. R., Oudmaijer, R. D., Mendigutía, I., Ilee, J. D., & van den Ancker, M. E. 2015, MNRAS, 453, 976 [NASA ADS] [CrossRef] [Google Scholar]
  16. Favre, C., Cleeves, L. I., Bergin, E. A., Qi, C., & Blake, G. A. 2013, ApJ, 776, L38 [NASA ADS] [CrossRef] [Google Scholar]
  17. Fedele, D., Bruderer, S., van Dishoeck, E. F., et al. 2013, A&A, 559, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Fedele, D., Carney, M., Hogerheijde, M. R., et al. 2017, A&A, 600, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Folsom, C. P., Bagnulo, S., Wade, G. A., et al. 2012, MNRAS, 422, 2072 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  20. Ginski, C., Stolker, T., Pinilla, P., et al. 2016, A&A, 595, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Guilloteau, S., Dutrey, A., Piétu, V., & Boehler, Y. 2011, A&A, 529, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Hales, A. S., De Gregorio-Monsalvo, I., Montesinos, B., et al. 2014, AJ, 148, 47 [NASA ADS] [CrossRef] [Google Scholar]
  23. Hartmann, L., Calvet, N., Gullbring, E., & D’Alessio, P. 1998, ApJ, 495, 385 [NASA ADS] [CrossRef] [Google Scholar]
  24. Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [Google Scholar]
  25. Isella, A., Testi, L., Natta, A., et al. 2007, A&A, 469, 213 [CrossRef] [EDP Sciences] [Google Scholar]
  26. Isella, A., Guidi, G., Testi, L., et al. 2016, Phys. Rev. Lett., 117, 251101 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kama, M., Folsom, C. P., & Pinilla, P. 2015, A&A, 582, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Kama, M., Bruderer, S., Carney, M., et al. 2016, A&A, 588, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Kenyon, S. J., & Hartmann, L. 1987, ApJ, 323, 714 [NASA ADS] [CrossRef] [Google Scholar]
  30. Kimura, S. S., & Tsuribe, T. 2012, PASJ, 64, 116 [NASA ADS] [CrossRef] [Google Scholar]
  31. Lagage, P.-O., Doucet, C., Pantin, E., et al. 2006, Science, 314, 621 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  32. Liu, S.-F., Jin, S., Li, S., Isella, A., & Li, H. 2018, ApJ, 857, 87 [NASA ADS] [CrossRef] [Google Scholar]
  33. Lodders, K. 2003, ApJ, 591, 1220 [NASA ADS] [CrossRef] [Google Scholar]
  34. Lynden-Bell, D., & Pringle, J. E. 1974, MNRAS, 168, 603 [NASA ADS] [CrossRef] [Google Scholar]
  35. Maaskant, K. M., Honda, M., Waters, L. B. F. M., et al. 2013, A&A, 555, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Maaskant, K. M., Min, M., Waters, L. B. F. M., & Tielens, A. G. G. M. 2014, A&A, 563, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Mannings, V., & Sargent, A. I. 2000, ApJ, 529, 391 [NASA ADS] [CrossRef] [Google Scholar]
  38. Marois, C., Macintosh, B., Barman, T., et al. 2008, Science, 322, 1348 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  39. Marois, C., Zuckerman, B., Konopacky, Q. M., Macintosh, B., & Barman, T. 2010, Nature, 468, 1080 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  40. McClure, M. K., Bergin, E. A., Cleeves, L. I., et al. 2016, ApJ, 831, 167 [NASA ADS] [CrossRef] [Google Scholar]
  41. Meeus, G., Waters, L. B. F. M., Bouwman, J., et al. 2001, A&A, 365, 476 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Meeus, G., Salyk, C., Bruderer, S., et al. 2013, A&A, 559, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Miotello, A., van Dishoeck, E. F., Kama, M., & Bruderer, S. 2016, A&A, 594, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Muto, T., Tsukagoshi, T., Momose, M., et al. 2015, PASJ, 67, 122 [NASA ADS] [CrossRef] [Google Scholar]
  45. Nakagawa, T., Shibai, H., Onaka, T., et al. 2014, SPIE, 9143, 91431I [Google Scholar]
  46. Panić, O., Hogerheijde, M. R., Wilner, D., & Qi, C. 2008, A&A, 491, 219 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Pascucci, I., Testi, L., Herczeg, G. J., et al. 2016, ApJ, 831, 125 [NASA ADS] [CrossRef] [Google Scholar]
  48. Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
  49. Pinilla, P., Birnstiel, T., & Walsh, C. 2015, A&A, 580, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Pinte, C., Price, D., Ménard, F., et al. 2018, ApJ, 860, L13 [NASA ADS] [CrossRef] [Google Scholar]
  51. Pinte, C., van der Plas, G., Menard, F., et al. 2019, Nat. Astron., 3, 1109 [Google Scholar]
  52. Poglitsch, A., Waelkens, C., Geis, N., et al. 2010, A&A, 518, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Powell, D., Murray-Clay, R., Pérez, L. M., Schlichting, H. E., & Rosenthal, M. 2019, ApJ, 878, 116 [NASA ADS] [CrossRef] [Google Scholar]
  54. Prodanović, T., Steigman, G., & Fields, B. D. 2010, MNRAS, 406, 1108 [NASA ADS] [Google Scholar]
  55. Qi, C., Öberg, K. I., Andrews, S. M., et al. 2015, ApJ, 813, 128 [NASA ADS] [CrossRef] [Google Scholar]
  56. Rameau, J., Chauvin, G., Lagrange, A. M., et al. 2013a, ApJ, 772, L15 [NASA ADS] [CrossRef] [Google Scholar]
  57. Rameau, J., Chauvin, G., Lagrange, A. M., et al. 2013b, ApJ, 779, L26 [NASA ADS] [CrossRef] [Google Scholar]
  58. Rhee, J. H., Song, I., Zuckerman, B., & McElwain, M. 2007, ApJ, 660, 1556 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  59. Richards, S. N., Moseley, S. H., Stacey, G., et al. 2018, J. Astron. Instrum., 7, 1840015 [CrossRef] [Google Scholar]
  60. Rosenfeld, K. A., Andrews, S. M., Hughes, A. M., Wilner, D. J., & Qi, C. 2013, ApJ, 774, 16 [NASA ADS] [CrossRef] [Google Scholar]
  61. Salyk, C., Herczeg, G. J., Brown, J. M., et al. 2013, ApJ, 769, 21 [NASA ADS] [CrossRef] [Google Scholar]
  62. Stassun, K. G., Oelkers, R. J., Pepper, J., et al. 2018, AJ, 156, 102 [NASA ADS] [CrossRef] [Google Scholar]
  63. Teague, R., Bae, J., Bergin, E. A., Birnstiel, T., & Foreman-Mackey, D. 2018, ApJ, 860, L12 [NASA ADS] [CrossRef] [Google Scholar]
  64. Tilling, I., Woitke, P., Meeus, G., et al. 2012, A&A, 538, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Toomre, A. 1964, ApJ, 139, 1217 [Google Scholar]
  66. Trapman, L., Miotello, A., Kama, M., van Dishoeck, E. F., & Bruderer, S. 2017, A&A, 605, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. van der Marel,N., van Dishoeck, E. F., Bruderer, S., Pérez, L., & Isella, A. 2015, A&A, 579, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. van der Marel,N., van Dishoeck, E. F., Bruderer, S., et al. 2016, A&A, 585, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. van der Plas, G., Ménard, F., Gonzalez, J. F., et al. 2019, A&A, 624, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Vioque, M., Oudmaijer, R. D., Baines, D., Mendigutía, I., & Pérez-Martínez, R. 2018, A&A, 620, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Walsh, C., Juhász, A., Meeus, G., et al. 2016, ApJ, 831, 200 [NASA ADS] [CrossRef] [Google Scholar]
  72. Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296 [NASA ADS] [CrossRef] [Google Scholar]
  73. Williams, J. P., & Best, W. M. J. 2014, ApJ, 788, 59 [NASA ADS] [CrossRef] [Google Scholar]
  74. Williams, J. P., & McPartland, C. 2016, ApJ, 830, 32 [NASA ADS] [CrossRef] [Google Scholar]
  75. Woitke, P., Kamp, I., Antonellini, S., et al. 2019, PASP, 131, 064301 [NASA ADS] [CrossRef] [Google Scholar]
  76. Woodall, J., Agúndez, M., Markwick-Kemper, A. J., & Millar, T. J. 2007, A&A, 466, 1197 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Zhang, S., Zhu, Z., Huang, J., et al. 2018, ApJ, 869, L47 [NASA ADS] [CrossRef] [Google Scholar]
  78. Zhang, K., Bergin, E. A., Schwarz, K., Krijt, S., & Ciesla, F. 2019, ApJ, 883, 98 [NASA ADS] [CrossRef] [Google Scholar]

1

Planets retrieved from exoplanets.org on 2019.07.16.

All Tables

Table 1

HD line flux upper limits (3σ) for the sample.

Table 2

DALI model grid parameters.

Table 3

Adopted model for HD 163296.

All Figures

thumbnail Fig. 1

HD1–0 line emitting regions in our flat/cold (blue) and flared/warm (red) disc models. Solid contours contain the middle 75% of cumulative line emission. Dashed lines indicate gas number density iso-contours for ngas = 106 cm−3, acting as a disc outline.

In the text
thumbnail Fig. 2

Distance-normalised 3σ upper limits on HD 112 μm line flux forthe disc sample (black lines and circles) compared with our grid of DALI disc models (coloured crosses). Highlighted crosses show the HD 112 μm line flux ofour fiducial model. The top panels show the group I sources compared to models with flaring angle ψ = 0.3. The bottom panels show the group II sources compared to models with ψ = 0.0. Left: models are separated based on gas mass. Right: HD 1–0 upper limitsset against 1.3 mm continuum fluxes for both observations and models.

In the text
thumbnail Fig. 3

HD 1–0 line and dust continuum flux dependences on disc and stellar parameters.

In the text
thumbnail Fig. 4

HD 1–0 line flux vs. stellar luminosity. Observed stellar luminosities taken from Table 1. Model stellar luminosities are given a small offset for clarity. Highlighted crosses show our fiducial model (hc = 0.15, Rc = 50 au, Δgd = 100, flarge = 0.95, χ = 0.2).

In the text
thumbnail Fig. 5

Comparing the HD 163296 specific models to the HD 1–0 upper limit (Fedele et al. 2013). All models have a dust mass Mdust = 6.6 × 10−4 M (Table 3). The red bar shows the range of gas masses inferred from CO in the literature.

In the text
thumbnail Fig. 6

Mass of selected discs and planets around B9 to F5 type stars. Vertical lines show the cumulative mass of each planetary system; dots highlight planets from the most massive at bottom. Disc gas mass upper limits from HD lines are taken from this work (HD 163296) and from Kama et al. (2016, HD 100546). For HD 163296, the range of CO isotopologue-based disc mass estimates is shown by a light blue bar (8 × 10−3 to 5.8 × 10−1 M; references in text). Also shown are the stellar mass divided by 10 and age; the mass limit for a gravitationally unstable disc (dashed line); an extrapolated dust-based disc mass range (dotted lines; Pascucci et al. 2016); and a population density colour map for planets around B9 to F5 type stars (data retrieved from exoplanets.org on 2019.07.16; bins contain from bottom to top 7, 6, and 1 planet). See text for individual planet and stellar mass references.

In the text
thumbnail Fig. 7

Observability of Group Ia,Ib (left) and Group IIa (right) models with SOFIA HIRMES. Coloured disc models are detectable (≥5σ) with a 10 h integration. Dark red dashed line shows the SPICA/SAFARI 1 h detection limit. The Origins Space Telescope 1 h detection limit (~1 × 10−20 W m−2) lies below the limits of the figure. The fluxes are calculated for a distance of 150 pc.

In the text
thumbnail Fig. A.1

Comparing the HD 135344B specific models from van der Marel et al. (2016) to HD 1–0 upper limit (Fedele et al. 2013). All models have a dust mass Mdust = 1.3 × 10−4 M (cf. Table 3 in van der Marel et al. 2016). The red circle shows the gas mass inferred from CO by van der Marel et al. (2016).

In the text
thumbnail Fig. B.1

Upper limits on HD 56 μm line flux for the sample of Herbig Ae/Be disc systems (black lines) compared with our grid of DALI disc models (crosses). The top panels show the group I sources compared to models with flaring angle ψ = 0.3. The bottom panels show the group II sources compared to models with ψ = 0.0. Left: modelsare separated based on gas mass. Right: HD 2–1 upper limits set against 1.3 mm continuum fluxes for both observationsand models.

In the text
thumbnail Fig. C.1

HD 1–0 line flux vs. 1.3 mm continuum fluxes for both observations and models. The panels are similar to right panels of Fig. 2, but also show the model gas-to-dust mass ratios (marker shape) and stellar luminosities (marker size).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.