Issue 
A&A
Volume 591, July 2016



Article Number  L3  
Number of page(s)  6  
Section  Letters  
DOI  https://doi.org/10.1051/00046361/201628549  
Published online  06 June 2016 
Evidence for a correlation between mass accretion rates onto young stars and the mass of their protoplanetary disks
^{1}
Scientific Support Office, Directorate of Science, European Space Research
and Technology Centre (ESA/ESTEC), Keplerlaan 1,
2201 AZ
Noordwijk,
The Netherlands
email:
cmanara@cosmos.esa.int
^{2}
Institute of Astronomy, University of Cambridge,
Madingley Road, Cambridge
CB30HA,
UK
^{3}
European Southern Observatory, KarlSchwarzschildStr. 2, 85748
Garching bei München,
Germany
^{4}
INAF/Osservatorio Astrofisico di Arcetri, Largo E. Fermi
5, 50125
Firenze,
Italy
^{5}
Excellence Cluster Universe, Boltzmannstr. 2, 85748
Garching bei München,
Germany
^{6}
School of Cosmic Physics, Dublin Institute for Advanced
Studies, 31 Fitzwilliams Place,
2
Dublin,
Ireland
^{7}
INAF/Osservatorio Astronomico di Capodimonte, Salita
Moiariello, 16
80131
Napoli,
Italy
^{8}
Institute for Astronomy, University of Hawai’i at
Mänoa, Honolulu,
HI,
USA
^{9}
Leiden Observatory, Leiden University,
PO Box 9513, 2300 RA
Leiden, The
Netherlands
^{10}
California Institute of Technology, 1200 East California Blvd, Pasadena, CA
91125,
USA
^{11}
University of Hawaii, Department of Physics and
Astronomy, 2505
Correa Rd. Honolulu, HI
96822,
USA
^{12}
Observatório Nacional/MCTI, 20921400, Rio de Janeiro, Brazil
^{13}
MaxPlankInstitut für Extraterrestrische Physik,
Giessenbachstraße 1,
85748
Garching,
Germany
Received: 18 March 2016
Accepted: 18 May 2016
A relation between the mass accretion rate onto the central young star and the mass of the surrounding protoplanetary disk has long been theoretically predicted and observationally sought. For the first time, we have accurately and homogeneously determined the photospheric parameters, mass accretion rate, and disk mass for an essentially complete sample of young stars with disks in the Lupus clouds. Our work combines the results of surveys conducted with VLT/XShooter and ALMA. With this dataset we are able to test a basic prediction of viscous accretion theory, the existence of a linear relation between the mass accretion rate onto the central star and the total disk mass. We find a correlation between the mass accretion rate and the disk dust mass, with a ratio that is roughly consistent with the expected viscous timescale when assuming an interstellar medium gastodust ratio. This confirms that mass accretion rates are related to the properties of the outer disk. We find no correlation between mass accretion rates and the disk mass measured by CO isotopologues emission lines, possibly owing to the small number of measured disk gas masses. This suggests that the mmsized dust mass better traces the total disk mass and that masses derived from CO may be underestimated, at least in some cases.
Key words: accretion, accretion disks / protoplanetary disks / stars: premain sequence / stars: variables: T Tauri, Herbig Ae/Be
© ESO, 2016
1. Introduction
The evolution of a protoplanetary disk significantly influences the planetary system that is formed. The final mass distribution of planets resembles the evolution of the surface density of gas in the disk (e.g., Thommes et al. 2008) and more massive disks lead to systems with more massive planets (Mordasini et al. 2012). The evolution of the disk structure is mainly driven by processes happening in the disk, such as dust evolution (Testi et al. 2014), and by interaction between the disk and the central star through viscous accretion and winds (Alexander et al. 2014).
In the context of viscously evolving protoplanetary disks, the mass accretion rate onto the central star (Ṁ_{acc}) and the mass of the disk (M_{disk}) should be directly correlated (e.g., Eq. (7) of Hartmann et al. 1998). The ratio between these quantities is related to the viscous timescale (t_{ν}) at the outer radius of the disk (R_{out}) and the assumptions about the disk viscous properties. Overall, it is expected that Ṁ_{acc}~M_{disk}/t_{ν}(R_{out}) with a coefficient of order unity (e.g., Jones et al. 2012). In disks that evolved viscously, the ratio M_{disk}/Ṁ_{acc} must be comparable to the age of the system independent of the initial conditions and value of the viscosity parameter α. The tight correlation between Ṁ_{acc} and stellar masses (Ṁ_{acc}∝M_{⋆}^{1.8}, e.g., Muzerolle et al. 2003; Natta et al. 2006; Alcalá et al. 2014; Manara et al. 2016) was also explained by Dullemond et al. (2006) as a consequence of the initial rotation rate of the cores where disks formed. They predict a strong dependence of M_{disk} on and a tight correlation of M_{disk} with Ṁ_{acc} as a consequence of viscous evolution. This theoretical relation between M_{disk} and Ṁ_{acc} has been empirically investigated, but previous studies were unable to find any significant correlation (e.g., Andrews et al. 2010; Ricci et al. 2010).
In this Letter we present a study of an almost complete and homogeneous dataset of young stars in the ~1–3 Myr old Lupus starforming region (Comerón 2008, d = 150–200 pc). We collected Ṁ_{acc} measured from ultraviolet (UV) excess with the VLT/XShooter spectrograph, and M_{disk} measured both from submm continuum and CO line emission with Atacama Large Millimeter/submillimeter Array (ALMA). We look for correlations between Ṁ_{acc} and M_{disk}, as predicted by viscous theory.
2. Data sample
The sample analyzed here includes Class II and transition disk (TD) young stellar objects (YSOs) with 0.1 < M_{⋆}/M_{⊙} < 2.2, thus including only the TTauri stars of the ALMA sample. Both the ALMA and XShooter surveys are complete at the ~95% level. In total, there are 66 objects with ALMA and XShooter data available. The list of targets included in the analysis is reported in Table A.1.
The ALMA data are presented by Ansdell et al. (2016, hereafter AW16). The setting includes continuum emission at 335.8 GHz (890 μm) at a resolution of ~0.34′′×0.28′′ (~25 × 20 AU radius at 150 pc) and the ^{13}CO and C^{18}O 3–2 transitions. From the continuum emission, detected for 54 of the targets included here, AW16 derive disk dust mass (M_{disk,dust}) using typical assumptions of a single dust grain opacity κ(890 μm) = 3.37 cm^{2}/g and a single dust temperature T_{dust} = 20 K. From the CO emission lines, AW16 derive the disk gas mass (M_{disk,gas}) for 29 disks in our sample, but for 22 of these the lower bound of the acceptable values of M_{disk,gas} is unconstrained because of the lack of C^{18}O detection. Upper limits are calculated for the other 37 targets.
We obtain Ṁ_{acc} from the XShooter spectra (Alcalá et al. 2014, and in prep.). Briefly, the stellar and accretion parameters are derived by finding the best fit among a grid of models including photospheric templates, a slab model for the accretion spectrum, and reddening (Manara et al. 2013). We use the UV excess as a main tracer of accretion and the broad wavelength range covered by XShooter (λλ ~ 330−2500) to constrain both the spectral type of the target and the extinction. Among the objects discussed here, 57 have Ṁ_{acc} derived from XShooter measurements, while 5 have an accretion rate compatible with chromospheric noise (non accretors), and 4 targets are observed edgeon, thus their Ṁ_{acc} are underestimated. The evolutionary models by Siess et al. (2000) are used to determine M_{⋆} and, thus, Ṁ_{acc}.
Finally, the sample includes several resolved binaries. All of these binaries have separations ≳2″ and are indicated in Table A.1.
3. Results
3.1. Disk dust mass
Fig. 1 Logarithm of Ṁ_{acc} vs. logarithm of M_{disk,dust}. Green filled squares are used for measured values, open squares for edgeon objects, and downward pointing open triangle for objects with accretion compatible with chromospheric noise. Transition disks are indicated with a circle. We show fit results obtained using the Bayesian fitting procedure by Kelly (2007), which considers errors on both axes and is only applied to detected targets. The assumed best fit is represented with a red solid line, while the light red lines are a subsample of the results of some chains. The best fitting with this procedure overlaps with the leastsquares best fit. 
The values of M_{disk,dust} derived by AW16 (see Sect. 2) are a measure of the bulk dust mass in the disk. Figure 1 shows the values of Ṁ_{acc} measured for our targets as a function of M_{disk,dust}^{1}.
We first search for a correlation between the two quantities running a leastsquares linear regression on the targets with both M_{disk,dust} and Ṁ_{acc} measurements and find a moderate correlation with r = 0.53 and a twosided pvalue of 1.5 × 10^{4} for the null hypothesis that the slope of this correlation is zero. The best fit obtained with this method has a slope of 0.7 and a standard deviation of the fit of 0.2. Then, we compute the linear regression coefficients using the fully Bayesian method by Kelly (2007)^{2}, which allows us to include uncertainties on both axes in the fitting procedure. Uniform priors are used for the linear regression coefficients. We include some single chain results in Fig. 1, as well as the best fit obtained with this method, which has the same slope and intercept of the leastsquares fit relation. We adopt the median of the results of the chains as best fit values. We refer to Appendix B for the corner plots with the posterior analysis results. The best fit obtained with this method has a slope of 0.7 ± 0.2, a standard deviation of 0.4 ± 0.1, and a correlation coefficient of 0.56 ± 0.12. We also verified that the two quantities are still correlated when upper limits on M_{disk,dust} are properly considered using the same tool. The correlation coefficient increases to 0.7 ± 0.1, while the slope is larger (1.2 ± 0.2) but compatible with that obtained using detections only. The same slope is obtained including upper limits and using the emmethod and buckleyjames method in ASURV. However, the slope estimated when including upper limits is not well constrained (Kelly 2007) and should be considered with caution. We then find a probability lower than 10^{4} of nocorrelation using the Cox hazard test for censored data in ASURV (Lavalley et al. 1992) including upper limits on M_{disk,dust}.
Fig. 2 Logarithm of L_{acc} vs. logarithm of continuum emission normalized to a distance of 150 pc. Symbols are as in Fig. 1. 
We show the dependence of the accretion luminosity (L_{acc}) on the submm continuum flux normalized to a distance of 150 pc in Fig. 2 to confirm that the correlation is not induced by the conversion from L_{acc} to Ṁ_{acc}. Indeed, L_{acc} is directly measured from the spectra, while the conversion to Ṁ_{acc} depends on M_{⋆}, which is derived from evolutionary models. A correlation is found with r = 0.6, a slope of 0.8 ± 0.2, and a standard deviation of 0.5 ± 0.1.
We then test for the robustness of the correlation, given our assumptions to convert the continuum emission in M_{disk,dust}. First, we assumed a single disk opacity and gastodust ratio for all disks. To test whether a random variation of these parameters would affect our results, we perform the same statistical tests on the same targets after randomly displacing the values of M_{disk,dust} within a uniform distribution with size ±1 dex and centered on the measured value to mimic the uncertainties. We perform this test ten times, and the correlation is still present in nine out of ten realizations. Then, we test the effect of our assumption of a single T_{dust} by modifying our M_{disk,dust} values assuming (Andrews et al. 2010). The correlation becomes less robust, but still significant (r = 0.3, pvalue = 0.03). We conclude that there is a statistically significant relation between the logarithm of Ṁ_{acc} and the logarithm of M_{disk,dust}. This relation has a slope slightly smaller than unity.
The location of TDs in Fig. 4 is also highlighted. All but one of the TDs are found to be below the bestfit relation in agreement with, for example, Najita et al. (2015). This suggests they have either lower Ṁ_{acc}, or larger disk mass, or a different gastodust ratio, than typical full disks.
3.2. Disk gas mass
Fig. 3 Logarithm of Ṁ_{acc} vs. logarithm of the disk mass derived from CO emission. Symbols are as in Fig. 1. No correlation is found between these quantities. 
Fig. 4 Logarithm of Ṁ_{acc} vs. logarithm of M_{disk}= 100 ·M_{disk,dust}. Symbols are as in Fig. 1. Also here, the best fit with the procedure by Kelly (2007) overlaps with the leastsquares best fit. The dashed lines represent different ratios of M_{disk}/Ṁ_{acc}, as labeled. 
The lower detection rates of CO lines than continuum emission (AW16) implies that we only measure M_{disk,gas} for a few objects. Figure 3 reports the Ṁ_{acc} vs. M_{disk,gas} plot.
We perform the same statistical tests as for the Ṁ_{acc} – M_{disk,dust} relation. We find no correlation between the logarithm of Ṁ_{acc} and the logarithm of M_{disk,gas} using the leastsquares linear regression on the targets with measured M_{disk,gas} (r = 0.2, pvalue = 0.3). When considering uncertainties on the measurements we find a value for the correlation coefficient of 0.5, and thus we find no correlation. We obtain the same statistically insignificant value for the correlation coefficient, which points toward no correlation, even when we include the upper limits. Finally, the Cox hazard test for censored data gives a probability of 0.25 that the two quantities are not correlated. We then conclude that we do not detect any correlation between these two quantities. The large number of upper limits compared to detection is probably a limiting factor in studying this relation and the large error bars of the measurements are another limiting factor. Deeper ALMA surveys of CO emission in protoplanetary disks are needed to further study this relation. In Fig. 3, the TDs are mixed with full disks.
4. Discussion
As mentioned in the Introduction, viscous evolution theory predicts that Ṁ_{acc}∝M_{disk}/t_{ν}(R_{out}). The evolution of the surface density of the disk (Σ) can be analytically described provided that the viscosity (ν) is known (e.g., Pringle 1981; Lodato 2008).
As described by Jones et al. (2012), a similarity solution (M_{disk} ∝ t^{− σ}) is reached for times much larger than the viscous timescales under the simple assumptions ν ∝ R^{n} or ν ∝ Σ^{m}R^{n}, where R is the disk radius. By differentiating this solution, one obtains Ṁ_{acc} ∝ σt^{− (1 + σ)} and thus it is possible to define the “viscous disk age” as t_{disk} = M_{disk}/Ṁ_{acc} = t/σ. Measurements of the decline of Ṁ_{acc} with time suggest that σ ~ 0.5 (e.g., Hartmann et al. 1998; SiciliaAguilar et al. 2010). For viscously evolving disks, this implies that the age of the objects should be within a factor ~2 of the ratio M_{disk}/Ṁ_{acc}. Jones et al. (2012) have also shown that more complex assumptions on the disk viscosity, such as different values for the α viscosity description, lead to the same asymptotic behavior with M_{disk}/Ṁ_{acc} ratios usually larger than the age of the objects by a factor 2–3, but always less than 10.
Other processes happening during disk evolution, such as a layered accretion, photoevaporation, and even planet formation, all lead to very similar values of M_{disk}/Ṁ_{acc} at late times, and these values are always higher than the age of the object (Jones et al. 2012). Thus, a disk that evolved only from internal processes has a M_{disk}/Ṁ_{acc} ratio similar or larger than its age regardless of the assumption on the disk viscosity. The only means by which a disk might have an M_{disk}/Ṁ_{acc} ratio smaller than its age is if the disk is externally truncated (Rosotti et al., in prep.).
We compare our results with these theoretical expectations by showing in Figs. 3 and 4 the M_{disk}/Ṁ_{acc} ratios for three different values of M_{disk}/Ṁ_{acc} = 0.1 Myr, 1 Myr, and 3 Myr. We assume that the total disk mass (M_{disk}) is M_{disk} = M_{disk,gas} in Fig. 3, while M_{disk} = 100·M_{disk,dust} in Fig. 4. Indeed, to convert M_{disk,dust} to M_{disk} one needs to know the gastodust ratio. We assume an interstellar medium value of 100 for the gastodust ratio, as is commonly done (e.g., Andrews et al. 2010; Ricci et al. 2010). If the gastodust ratio has no dependence on M_{disk}, this has no impact on the correlation between Ṁ_{acc} and M_{disk} but is instructive for the discussion. The typical age of Lupus targets is ~1–3 Myr with a spread of 1–2 Myr (e.g., Alcalá et al. 2014).
The location of the targets in Fig. 4 is in general agreement with the aforementioned theoretical expectations. Most of the targets (60%) have positions between or compatible with the 1 and 3 Myr lines. However, several of the targets in Fig. 3 do not match the expectations from viscous evolution theory as they lie above the M_{disk}/Ṁ_{acc} = 1 Myr line.
The lack of correlation between Ṁ_{acc} and M_{disk,gas} is contrary to expectations from viscous evolution theory. When assuming M_{disk} = 100 · M_{disk,dust}, however, we find a correlation between M_{disk} and Ṁ_{acc} and also M_{disk}/Ṁ_{acc} ratios that are compatible with expectations from theory. Thus, we are inclined to conclude that the total disk mass M_{disk} ∝ M_{disk,dust}, as with this assumption the correlation is present. This in turn suggests that M_{disk,gas} measured from CO emission is possibly lower than the total M_{disk}, at least for the more massive disks. A possible explanation for this might be that carbon is processed in more complex molecules (e.g., Bergin et al. 2014; Kama et al. 2016) or that more detailed modeling of CO lines is needed, but this discussion is out of the scope of this paper.
The slope of the observed correlation between Ṁ_{acc} and M_{disk}, as measured from dust emission, is consistent with being linear, as expected if all disks evolve viscously, however, we cannot exclude that it is actually shallower. The exact slope can be derived with a better handle on the uncertainty in the M_{disk} estimate, namely the gastodust ratio, disk grain opacity, disk temperature, and their dependence on the stellar properties. More constraints on these values are awaited from future ALMA survey of disks with higher sensitivity, multiple band observations, and targeting several molecules in order to better determine the chemical properties of the disks. The interest in further constraining this slope is related to the fact that this relation can tell us what evolutionary processes dominate at different stellar masses.
5. Conclusions
In this Letter we compared the most complete and homogeneous datasets of properties of young stars and their disks to date. We used accretion rates onto the central star determined from UV excess with the VLT/XShooter spectrograph and disk masses from both submm continuum and CO line emission measured by ALMA.
We detected a statistically significant correlation between Ṁ_{acc} and M_{disk} with a slope that is slightly smaller than 1. This is found when assuming that the total disk mass is proportional to the disk dust mass, but not when using the disk gas mass. The latter result could be due to large uncertainties in M_{disk,gas} estimate and low number statistics. For this reason, deeper surveys of gas emission in disks are needed. When measuring M_{disk} from dust emission, transitional disks are found to have either a smaller Ṁ_{acc} or a larger M_{disk} than full disks.
We compared the observed M_{disk} derived from dust emission and Ṁ_{acc} with basic predictions from viscous evolution theory and we found a good agreement with the expected M_{disk}/Ṁ_{acc} ratios for our targets.
Future studies should look for the M_{disk} – Ṁ_{acc} correlation for objects with different ages and in different environments.
Acknowledgments
We thank Cathie Clarke and Phil Armitage for insightful discussions. We thank the anonymous referee for insightful comments that helped to improve the presentation of the results. C.F.M. gratefully acknowledges an ESA Research Fellowship. G.R. is supported by the DISCSIM project, grant agreement 341137 funded by the European Research Council under ERC2013ADG. A.N. would like to acknowledge funding from Science Foundation Ireland (Grant 13/ERC/I2907). Leiden is supported by the European Union AERC grant 291141 CHEMPLAN, by the Netherlands Research School for Astronomy (NOVA), and by grant 614.001.352 from the Netherlands Organization for Scientific Research (NWO). JPW and MA were supported by NSF and NASA grants AST1208911 and NNX15AC92G.
References
 Alcalá, J. M., Natta, A., Manara, C. F., et al. 2014, A&A, 561, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Alexander, R., Pascucci, I., Andrews, S., Armitage, P., & Cieza, L. 2014, Protostars and Planets VI, 475 [Google Scholar]
 Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C., & Dullemond, C. P. 2010, ApJ, 723, 1241 [NASA ADS] [CrossRef] [Google Scholar]
 Ansdell, M., Williams, J. P., van der Marel, N., et al. 2016, ArXiv eprints [arXiv:1604.05719] [Google Scholar]
 Bergin, E. A., Cleeves, L. I., Crockett, N., & Blake, G. A. 2014, Faraday Discussions, 168, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Comerón, F. 2008, Handbook of Star Forming Regions, 5, 295 [Google Scholar]
 Dullemond, C. P., Natta, A., & Testi, L. 2006, ApJ, 645, L69 [NASA ADS] [CrossRef] [Google Scholar]
 Hartmann, L., Calvet, N., Gullbring, E., & D’Alessio, P. 1998, ApJ, 495, 385 [NASA ADS] [CrossRef] [Google Scholar]
 Jones, M. G., Pringle, J. E., & Alexander, R. D. 2012, MNRAS, 419, 925 [NASA ADS] [CrossRef] [Google Scholar]
 Kama, M., Bruderer, S., Carney, M., et al. 2016, A&A, 588, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kelly, B. C. 2007, ApJ, 665, 1489 [NASA ADS] [CrossRef] [Google Scholar]
 Lavalley, M., Isobe, T., & Feigelson, E. 1992, Astronomical Data Analysis Software and Systems I, 25, 245 [Google Scholar]
 Lodato, G. 2008, New Astron. Rev., 52, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Manara, C. F., Beccari, G., Da Rio, N., et al. 2013, A&A, 558, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Manara, C. F., Fedele, D., Herczeg, G. J., & Teixeira, P. S. 2016, A&A, 585, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mordasini, C., Alibert, Y., Benz, W., Klahr, H., & Henning, T. 2012, A&A, 541, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Muzerolle, J., Hillenbrand, L., Calvet, N., Briceño, C., & Hartmann, L. 2003, ApJ, 592, 266 [NASA ADS] [CrossRef] [Google Scholar]
 Natta, A., Testi, L., & Randich, S. 2006, A&A, 452, 245 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Najita, J. R., Andrews, S. M., & Muzerolle, J. 2015, MNRAS, 450, 3559 [NASA ADS] [CrossRef] [Google Scholar]
 Pringle, J. E. 1981, ARA&A, 19, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Ricci, L., Testi, L., Natta, A., et al. 2010, A&A, 512, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 SiciliaAguilar, A., Henning, T., & Hartmann, L. W. 2010, ApJ, 710, 597 [NASA ADS] [CrossRef] [Google Scholar]
 Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 [Google Scholar]
 Testi, L., Birnstiel, T., Ricci, L., et al. 2014, Protostars and Planets VI, 339 [Google Scholar]
 Thommes, E. W., Matsumura, S., & Rasio, F. A. 2008, Science, 321, 814 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
Appendix A: Sample
Targets included in the analysis.
Appendix B: Corner plots of fit
We show in Fig. B.1 the corner plots showing the intercept, slope, standard deviation, and correlation coefficient for the fit
Fig. B.1 Corner plot of the log Ṁ_{acc} vs. log(100·M_{disk,dust}) fit done using the linmix routine (Kelly 2007). 
of the log Ṁ_{acc} vs. log(100 · M_{disk,dust}) relation (see Sect. 3.1 and Fig. 4). This is done using linmix (Kelly 2007), only data with detection on both axes, and including the errors on both axes.
All Tables
All Figures
Fig. 1 Logarithm of Ṁ_{acc} vs. logarithm of M_{disk,dust}. Green filled squares are used for measured values, open squares for edgeon objects, and downward pointing open triangle for objects with accretion compatible with chromospheric noise. Transition disks are indicated with a circle. We show fit results obtained using the Bayesian fitting procedure by Kelly (2007), which considers errors on both axes and is only applied to detected targets. The assumed best fit is represented with a red solid line, while the light red lines are a subsample of the results of some chains. The best fitting with this procedure overlaps with the leastsquares best fit. 

In the text 
Fig. 2 Logarithm of L_{acc} vs. logarithm of continuum emission normalized to a distance of 150 pc. Symbols are as in Fig. 1. 

In the text 
Fig. 3 Logarithm of Ṁ_{acc} vs. logarithm of the disk mass derived from CO emission. Symbols are as in Fig. 1. No correlation is found between these quantities. 

In the text 
Fig. 4 Logarithm of Ṁ_{acc} vs. logarithm of M_{disk}= 100 ·M_{disk,dust}. Symbols are as in Fig. 1. Also here, the best fit with the procedure by Kelly (2007) overlaps with the leastsquares best fit. The dashed lines represent different ratios of M_{disk}/Ṁ_{acc}, as labeled. 

In the text 
Fig. B.1 Corner plot of the log Ṁ_{acc} vs. log(100·M_{disk,dust}) fit done using the linmix routine (Kelly 2007). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.