ALMA survey of Class II protoplanetary disks in Corona Australis: a young region with low disk masses

In recent years, the disk populations in a number of young star-forming regions have been surveyed with ALMA. Understanding the disk properties and their correlation with those of the central star is critical to understand planet formation. In particular, a decrease of the average measured disk dust mass with the age of the region has been observed. We conducted high-sensitivity continuum ALMA observations of 43 Class II young stellar objects in CrA at 1.3 mm (230 GHz). The typical spatial resolution is 0.3". The continuum fluxes are used to estimate the dust masses of the disks, and a survival analysis is performed to estimate the average dust mass. We also obtained new VLT/X-Shooter spectra for 12 of the objects in our sample. 24 disks are detected, and stringent limits have been put on the average dust mass of the non-detections. Accounting for the upper limits, the average disk mass in CrA is $6\pm3\,\rm M_\oplus$, significantly lower than that of disks in other young (1-3 Myr) star forming regions (e.g. Lupus) and appears consistent with the 5-10 Myr old Upper Sco. The position of the stars in our sample on the HR diagram, however, seems to confirm that that CrA has age similar to Lupus. Neither external photoevaporation nor a lower than usual stellar mass distribution can explain the low disk masses. On the other hand, a low-mass disk population could be explained if the disks are small, which could happen if the parent cloud has a low temperature or intrinsic angular momentum, or if the the angular momentum of the cloud is removed by some physical mechanism such as magnetic braking. In order to fully explain and understand the dust mass distribution of protoplanetary disks and their evolution, it may also be necessary to take into consideration the initial conditions of star and disk formation process, which may vary from region to region, and affect planet formation.


Introduction
Planets form in protoplanetary disks around young stars, and the way these disks evolve also impacts what kind of planetary system will be formed (Morbidelli & Raymond 2016). The evolution of the disk mass with time is one of the key ingredients of planetary synthesis models (Benz et al. 2014). For a long time infrared telescopes (e.g., Spitzer) have shown how the inner regions of disks dissipate on a timescale of ∼3-5 Myr (Haisch et al. 2001;Hernández et al. 2007;Fedele et al. 2010;Bell et al. 2013).
Only recently, however, we have been able to measure the bulk disk mass for statistically significant samples of disks, thanks to the high sensitivity of the Atacama Large Millimeter/submillimeter Array (ALMA). Pre-ALMA surveys of disk Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme ID 299.C-5048 and 0101.C-0893 masses were restricted to the northern hemisphere Taurus, Ophiuchus and Orion Nebula Cluster regions (Andrews & Williams 2005;Andrews et al. 2009Andrews et al. , 2013Eisner et al. 2008;Mann & Williams 2010). In the first years of operations of ALMA this has changed dramatically: hundreds of disks have been surveyed to determine the disk population in the ∼1-3 Myr old Lupus, Chamaeleon I, Orion Nebula Cluster, Ophiuchus, IC348 and Taurus regions (Ansdell et al. 2016;Pascucci et al. 2016;Eisner et al. 2018;Cieza et al. 2019;Ruíz-Rodríguez et al. 2018;Long et al. 2018), in the ∼3-5 Myr old σ-Orionis region (Ansdell et al. 2017), and in the older ∼5-10 Myr Upper Scorpius association (Barenfeld et al. 2016). These surveys have shown that the typical mass of protoplanetary disks decreases with the age of the region, in line with the observations that the inner regions of disks are dissipated within ∼ 3-5 Myr, similar to the dissipation time scale measured in the infrared. A positive correlation between disk and stellar mass was also found, and a steepening of its slope with time was identified (Ansdell et al. 2016(Ansdell et al. , 2017 . This is consistent with the result that massive planets form and are found preferentially around more massive stars (e.g. Bonfils et al. 2013;Alibert et al. 2011). Finally, the steepening of the relation with time is explained with more efficient radial drift around low mass stars (Pascucci et al. 2016), and it suggests that a significant portion of the planet formation process, especially around low mass stars, must happen in the first ∼1-2 Myr, when enough material to form planets is still available in disks Manara et al. 2018). Studying the evolution of the M disk − M relation in as many different environments as possible is therefore critical for understanding how the planet formation process is affected by the mass of the central stars.
We present here a survey of the Class II disks in the Corona Australis star forming region (CrA). Located at an average distance of about 154 pc (Gaia Collaboration et al. 2018;Dzib et al. 2018), the CrA molecular cloud complex is one of the nearest star-forming regions (see review in Neuhäuser & Forbrich 2008). It has been the target of many infrared surveys, the most recent being the Gould Belt (GB) Spitzer Legacy program presented in Peterson et al. (2011). At the center of the CrA region is located the Coronet cluster, which is a region of young embedded objects in the vicinity of R CrA (Herbig Ae star, Neuhäuser et al. 2000), on which many of the previous studies have focused. All studies agree in assigning to the Coronet an age < 3 Myr (e.g. Meyer & Wilking 2009;Sicilia-Aguilar et al. 2011). However, there are also some indications of a more evolved population (e.g. Neuhäuser et al. 2000;Peterson et al. 2011;Sicilia-Aguilar et al. 2011). A deep, sub-mm wavelength survey of the disk population in the region can help to further understand the formation and evolutionary history of CrA.
We therefore use ALMA to conduct a high-sensitivity millimeter wavelength survey of all the known Class II sources in CrA and compare the results with other regions surveyed todate. In Sec. 2 the sample is described, while the ALMA observations are detailed in Sec. 3. We also describe there new VLT/X-Shooter observations to determine the stellar charachteristics. The continuum millimeter measurements, their conversion to dust masses and a comparison with other star-forming regions is presented in Sec. 4. Our findings are interpreted in the context of disk evolution in Sec. 5. Finally, the work is summarized in Sec. 6.   (2014), (12) Patten (1998) of 116 YSOs, 14 of which are classified as Class I, 5 as Flat Spectrum (FS), 43 as Class II and 54 as Class III. The Infrared Class was determined by calculating the spectral slope α over the widest possible range of IR wavelengths as follows: where λ is the wavelength and F λ the flux at λ. Sources with α ≥ 0.3 are classified as Class I; FS have −0.3 ≤ α < 0.3; Class II have −1.6 ≤ α < −0.3; sources with α < −1.6 are Class III (Evans et al. 2009;Peterson et al. 2011). Fig. 1 shows the spatial distribution of the sources and their classification on top of the Herschel 250 µm map of the molecular Cloud.
Our sample includes all the Class II sources from the Peterson et al. (2011) catalogue. Two of them  were later identified as background, evolved stars based on parallax measurements with Gaia (Gaia Collaboration et al. 2016;Lindegren et al. 2018;Luri et al. 2018;Gaia Collaboration et al. 2018) and on our VLT/X-Shooter spectra (see Sec. 3.2). We then checked our sample against the more recently published Article number, page 3 of 20  Dunham et al. (2015) in which the Spitzer data are re-analysed and the spectral slopes re-calculated. We find broad agreement between the classification in Peterson et al. (2011) and Dunham et al. (2015), except for a few very marginal cases at the boundaries of classes.
Our final sample contains 41 targets, two of which are clearly resolved binaries (S CrA and CrA-45). Of the 43 targeted disks, 24 are detected with ALMA. The spectral type (SpT) was known for only 26 of the stars from the literature. We obtained VLT/X-Shooter spectra for 11 of the remaining targets, and derived their properties as explained in Sec. 4.1.
The basic stellar properties for the CrA sample are given in Table 1, the distribution of SpTs is shown in Fig. 2, while the millimeter observations, flux densities, and calculated disk masses are presented in Table 3.  Table 2. 12 CO (2-1), 13 CO (2-1) and C 18 O (2-1) transitions were also targeted with our spectral setup, but no clear detection was found because of strong foreground contamination. SO (6-5) and SiO  lines were also covered and not detected.

Observations
The data were manually calibrated using the CASA v5.1.1 software package (McMullin et al. 2007) . The gain calibrator for the first epoch of observations was faint. To yield reasonably high signal-to-noise (S/N) ratios when deriving the gain phase solutions, the phase offsets among spectral windows were first solved using the passband calibration scan. After applying the phase offsets solution, the gain phase solution was then derived by combining all spectral windows. The calibration of the other two epochs of observations followed the standard procedure of ALMA quality assurance (i.e., QA2). The bootstrapped flux values of the calibrator quasar J1924-2914 were consistent with the SMA Calibrator list 1 (Gurwell et al. 2007) to ∼10%. After calibration, we fit the continuum baseline and subtract it from the spectral line data, using the CASA task uvcontsub.
The continuum data imaging was performed with multifrequency synthesis (MFS) imaging of the continuum data using the CASA-clean task, and correcting for the primary beam. By jointly imaging all three epochs of data, for each target source field, the achieved continuum root-mean-square (RMS) noise level is ∼0.15 mJy beam −1 , and the synthesized beam is θ maj × θ min =0 . 33×0 . 31 (P.A.=67 • ), corresponding to a spatial resolution of ∼ 50 au at d = 154 pc. The imaged detections are presented in Fig. 3.
It is important to note that because of an error when setting the observation coordinates, the decimal places of the target RAs have been trimmed: this results in an offset of the sources of up to 15 east of the phase center: as a consequence, our images had to be primary beam corrected. The images in Fig. 3 have therefore been re-centered using the best-fit positions in Tab. 3.

VLT/X-Shooter observations
The spectroscopic follow-up observations for the 13 targets with missing spectral type information were carried out in Pr.Id. 299.C-5048 (PI Manara) and Pr.Id. 0101.C-0893 (PI Cazzoletti) with the VLT/X-Shooter spectrograph (Vernet et al. 2011). This instrument covers the wavelength range from ∼300 nm to ∼2500 nm simultaneously, dividing the spectrum in three arms, the UVB (λλ ∼ 300-550 nm), the VIS (λλ ∼ 500-1050 nm), and the NIR (λλ ∼ 1000-2500 nm). All targets were observed both with a narrow slit -1.0 in the UVB, 0.9 in the VIS and NIR arms -leading to R∼9000 and ∼10000, respectively, and a wide slit of 5.0 used to obtain an accurate flux calibration of the spectra. The log of the observations is reported in Table B.2. The spectra of all the observed targets are detected in the NIR arm, while only 5 targets are bright enough and not extincted too much to be detected also in the UVB arm.
The reduction of the data was performed using the ESO X-Shooter pipeline 2.9.3 (Modigliani et al. 2010). The pipeline performs the typical reduction steps, such as flat fielding, bias subtraction, order extraction and combination, rectification, wavelength calibration, flux calibration using standard stars observed in the same night. We extracted the 1D spectra from the 2D images produced by the pipeline using IRAF and then removed telluric absorption lines in the VIS and NIR arms using telluric standard stars observed close in time and airmass (see e.g., Alcalá et al. 2014). The S/N of the spectra at different wavelengths is reported in

Stellar properties
The spectral type for the targets were obtained from the literature (see Tab. 1) or from the VLT/X-Shooter spectra. The procedure used for the analysis of the X-Shooter spectra was as follows. First, we corrected the spectra for extinction using the values from the literature (Dunham et al. 2015;Sicilia-Aguilar et al. 2008 and the reddening law by Cardelli et al. (1989) with R V =3.1, as suggested by Sicilia-Aguilar et al. (2008). Then,   we calculated the values of a number of spectral indices both at wavelengths in the VIS and the NIR arms, taken by those calibrated by Riddick et al. (2007), Jeffries et al. (2007), and Herczeg & Hillenbrand (2014), as in Manara et al. (2017), and by Testi et al. (2001), as in Manara et al. (2013a). The spectral types derived from these indices are presented in Tab. B.1 in Appendix B. The spectral indices in the VIS arms are more reliable, and we select the spectral type from these indices when available. The observed spectra along with a template of the relative Spectral Types are presented in Fig. B.1.
The spectral types are converted in effective temperatures (T eff ) using the relation by Herczeg & Hillenbrand (2014). Stellar luminosity (L ) is obtained from the reddening-corrected J-band magnitudes and using the bolometric correction from Herczeg & Hillenbrand (2014), assuming for all the target the average distance of 154 pc calculated by Dzib et al. (2018). With this information, we have been able to plot our data on the HR diagram ( Fig. 6) and to estimate the stellar masses (M ) for all the targets using the evolutionary tracks by Baraffe et al. (2015) for M < 1.4M and Siess et al. (2000) for higher M and ages younger than 1 Myr . The stellar parameters for the targets are reported in Tab

mm continuum emission
Among the 41 targets, 20 of them show a clear (≥ 4σ) detection within a 1 radius from the nominal Spitzer location from Peterson et al. (2011). In addition, CrA-42 and T CrA show a ∼ 36σ and a ∼ 22σ detection respectively at a slightly larger distance from their nominal Spitzer positions (1 .05 for CrA-42 and 1 .34 T CrA), and are also regarded as detections. S CrA is a known binary (Reipurth & Zinnecker 1993;Ghez et al. 1997;Takami et al. 2003), and we detected millimeter emission associated with both binary components. CrA-45 is also identified as a binary. The total number of detections is therefore 24 out of the 43 targeted disks, so the detection rate is ∼ 56% .
None of the disks show clear substructures, no transition disk with cavities with radius > 25 au are found and all of them appear to be unresolved or marginally resolved: a Gaussian is therefore fitted to the detected sources (two Gaussians for the binaries) in the image plane using the imfit task in CASA. The task returns the total flux-density F 1.3 mm of the source along with the statistical uncertainty, the FWHM along the semi-major (a maj ) and semi-minor (a min ) axis and the position angle (PA). The results of the fit are shown in Table 3 2 The right ascension offset (∆α) and the declination offset (∆δ) with respect to the Spitzer coordinates is also shown. The rms noise for the non-detections was calculated using the imstat task within a 1 radius centered at the Spitzer coordinates; for the detection, it was calculated in an annular region centered on the source and with inner and outer radii equal to 2 and 4 , respectively.
In order to constrain the average flux density of individually undetected sources, a stacking analysis was also performed. The images were centered at their Spitzer coordinates (Table 1) and then stacked. Even after the stacking, no detection was found and an average rms noise is 0.017 mJy beam −1 , corresponding to a 3σ upper limit of 0.051 mJy is found assuming unresolved disks. However, it should be noted that the average offset between the disks and the Spitzer positions, measured on the detections, are < ∆α >= −0.13 and < ∆δ >= 0.47 : it is therefore possible that the undetected sources did not overlap during the stacking, and that the upper limit is actually higher than that quoted.

Dust masses
Assuming that the observed sub-millimeter emission is optically thin and isothermal, the relation between the emitting dust mass (M dust ) and the observed continuum flux at frequency ν (F ν ) is as follows (Hildebrand 1983): where d is the distance of the object, F ν is measured the fluxdensity, B ν (T dust ) is the Planck function for a given dust temperature T dust and κ ν is the dust opacity at frequency ν. To make the comparison with previous surveys easier, for the dust opacity κ ν we follow the same approach of Ansdell et al. (2016), assuming κ ν = 10 cm 2 g −1 at 1000 GHz (Beckwith et al. 1990) and scaling it to our frequency using β = 1. The adopted value is therefore κ ν = 2.3 cm 2 g −1 at ν = 230 GHz (1.3 mm). In the right-hand side of Eq. 2, the distance d is measured in pc and the flux density 2005), rather than the T dust = 25 K×(L * /L ) 0.25 relation based on two-dimensional continuum radiative transfer by Andrews et al. (2013) and used in other works (e.g. Law et al. 2017). We adopt this simplified approach with a single grain opacity and temperature for all the disks in the sample following the approach of Ansdell et al. (2016) and to facilitate the comparison with other star-forming regions (see Sec. 4.4). Moreover, it should be noted that no dependence of the average dust temperature on the stellar parameters was found with the more detailed modelling by Tazzari et al. (2017) for the Lupus disks.
The dust masses of the disks in our sample are presented in Tab. 3, along with the relative uncertainty calculated from the flux uncertainty. Only 3 disks out of 24 detections have a dust content ≥ 10 M ⊕ 3 and large enough to form the cores of giant planets in the future. However, it is still possible that a similar amount of dust mass is hidden at the inner few region due to very high optical depth (e.g. Zhu et al. 2010;Liu et al. 2017;Vorobyov et al. 2018). Also note that very recent high angular resolution ALMA and VLA observations of disks are revealing that an important amount of dust is located in dense regions such as rings (e.g. Andrews et al. 2018), which are optically thick at wavelengths around 1 mm (Dullemond et al. 2018). When optically thin emission is detected, higher masses are estimated (Carrasco-González et al. 2016).
The stacking of the non-detections gives an average 3σ upper limit corresponding to 0.036 M ⊕ , about 3 Lunar masses.

Comparison with other regions
The surveys of nearby star forming regions over the last years have shown growing evidence of a decrease in the mass of the disks with age, reflecting dust growth and disk dispersal. Ansdell et al. (2016Ansdell et al. ( , 2017 found consistent results, calculating the highest average mass in the youngest regions (1)(2)(3), and the lowest for the the oldest Upper Sco association (5-10 Myrs). The 2-3 Myrs old IC348 is the only exception, showing an average dust mass of only 4 ± 1 M ⊕ (between the average σ Orionis and that of Upper Sco) despite its young age. This can be explained by the low-mass stellar population in the region (Ruíz-Rodríguez et al. 2018) (also see Tab. 4).
The same analysis was done here for CrA. The dust masses are uniformly calculated following the approach used by Ansdell et al. (2016), namely using Eq. 2 with the continuum fluxes (or the 3 σ upper limits) from our ALMA data or from the literature, assuming a uniform T = 20 K, and inputting the frequency of the observation for each specific dataset. The distances assumed for each region are listed in Tab. 4. For the Upper Sco region, only the disks classified as "full", "evolved" and "transitional" from the Barenfeld et al. (2016) sample are included, while the "debris" and Class III YSOs, which likely represent a separate evolutionary stage, are excluded. Finally, in order to facilitate the comparison with the other samples, in this analysis we only include the disks around stars with masses above the brown-dwarf limit (M ≥ 0.1 M ). The Kaplan-Meier estimator from the lifelines 4 and ASURV (Lavalley et al. 1992) packages were then used to estimate the cumulative mass distribution and to calculate the average dust mass and its uncertainty while properly accounting for the upper limits by using wellestablished techniques for left-censored data sets. Fig. 4 presents the results accounting for the upper limits given by the non-detections. With an average dust mass of 6 ± 3 M ⊕ , the distribution of the CrA disks appear closer to that of the old Upper Sco region rather than to those of the younger systems.

M disk − M relation
A clear correlation between the dust mass of disks and the mass of the central star has been identified across all protoplanetary disk populations surveyed (Pascucci et al. 2016;Ansdell et al. 2017). This finding highlights how the disk properties are affected by the central star, and is consistent with the correlation between frequency of giant planets and mass of the host star, both from the observational and theoretical points of view (Alibert et al. 2011;Bonfils et al. 2013). Moreover, the slope of this relation has been observed to steepen with time, with the young Taurus, Lupus and Chemeleon I regions (∼ 1 − 3 Myr) having slopes similar to each other and shallower than that found for the disks in the Upper Sco association ( low mass of the disk population found in Sec. 4.4. We derive the M dust −M relation using the same linear-regression Bayesian approach followed by Ansdell et al. (2017) and presented by Kelly (2007) 5 . Unlike other linear regression methods, this approach is capable of simultaneously accounting for the uncertainties in both the measurements of M dust and M , of the intrinsic scatter of the data and of the disk non detections, which result in upperlimits on the disk masses. Note that the SpT, and therefore the stellar mass, is missing for 5 of our targets: for these objects the stellar mass is randomly drawn from the stellar mass distribution of the entire sample. In particular, 4 of the objects with unknown SpT are also not detected with ALMA, while the other one (CrA-42) shows a clear detection of a disk at mm-wavelengths. For the 4 non detections, the stellar mass is therefore randomly drawn among the masses of the stars with non-detected disks, while the mass of CrA-42 is drawn from those showing a detection with ALMA. This uncertainty is also taken into account in the Bayesian approach we adopt by performing 100 different draws. In our fit, a standard uncertainty of 20% of M on the stellar mass is assumed (Alcalá et al. 2017;Manara et al. 2017), while the uncertainties shown in Tab. 3 were used for the M dust values. Finally, it should be noted that only 1 out of 89 sources in Lupus was a Herbig Ae/Be star, while Upper Sco did not include any Herbig. We therefore decided not to include T CrA and TY CrA in the fit , for which the M dust − M relation might not hold.
The best fit relation we find is then plotted in Fig. 5 in dark red, along with a subsample of all the models in the chains to show the uncertainty. As in the other surveys, we also find a correlation, where the best-fit model has a slope β = 2.32 ± 0.77 and intercept α = 1.29 ± 0.60. This regression intercept is lower that that of other regions, as a consequence of the low disk masses found in the region. The uncertainties of the best-fit parameters reflect the large scatter in the data and the low number statistics.
In order to test that no strong bias was introduced by our procedure, we also run the fit described above without any random draw, finding consistent results.  Barenfeld et al. (2016). The evolutionary tracks for different stellar masses and the relative isochrones from Baraffe et al. (2015) are also plotted for reference. The isochrones refer (from top to bottom) to the 1 Myr, 2Myr, 5 Myr and 10 Myr isochrone. The coloured solid lines show the approximate median value of the luminosity at each temperature.

Is CrA old?
The observed low disk dust masses suggest that the CrA objects targeted in our survey may have an age comparable to that of the Upper Sco association, rather than to the young Lupus region. Unlike CrA, however, Upper Sco shows no presence of Class 0 or Class I sources, as expected for a 5 − 10 Myr region (Dunham et al. 2015). Moreover, most studies agree in assigning Corona Australis an age < 3 Myr (e.g. Meyer & Wilking 2009;Nisini et al. 2005;Sicilia-Aguilar et al. 2008. On the other hand, most of these studies focused only on the Coronet cluster, a small region extending ∼ 1 pc around the R CrA YSO, and where most of the young embedded Class 0 and Class I sources are located (see Fig. 1). The hypothesis that the large scale YSO population of the whole CrA cloud also includes a population of older objects therefore cannot be entirely ruled out. Some evidence of an additional older population has already been presented in previous studies. Neuhäuser et al. (2000) for example identify two classical T Tauri stars located outside the main cloud with an age of ∼ 10 Myr using ROSAT data. In addition, Peterson et al. (2011) perform a clustering analysis of the 116 YSOs in their sample, identifying a single core (corresponding to the Coronet) and a more extended population of PMS stars showing an age gradient west of the Coronet. They also observe that in the central core, the ratio Class II/ Class I=1.8, while the same ratio is Class II/ Class I=2.3 when all the objects in the sample are considered, again hinting toward a younger population inside the Coronet. In order to further test if the Class II population in our sample indeed includes an older population, we have placed them on the HR diagram, by using the spectral types listed in Tab. 1 and by deriving effective temperatures and bolometric corrections using the relationships in Herczeg & Hillenbrand (2014) and tables in Herczeg & Hillenbrand (2015), respectively. The obtained diagram is presented in Fig. 6. For comparison, the Upper Sco and Lupus objects are also plotted. In contrast with what Fig. 4 suggests, the HR diagram supports the scenario of a young CrA cluster with an age more consistent to that of Lupus than to Upper Sco.
In order to make this conclusion evident, the median values of the bolometric luminosities for each temperature are also shown (solid coloured lines in Fig. 6). The indicative age of the cluster is the isochrone closer to those median values: these lines also suggest that CrA is younger than Upper Sco. However, a more extended spectral classification for a larger number of objects in CrA would be needed to fully test this older-population scenario.

Is CrA young?
If the whole CrA is coeval with an age of 1 − 3 Myr, some other mechanism has to be invoked to explain the low observed mm fluxes. For example, these fluxes could be due to low metallicity. However, James et al. (2006) determined metallicities for three T Tauri stars in CrA, finding them to be only slightly sub-solar, and not low enough to explain our obesrvations.
External photo-evaporation is also known to play an important role in the disk mass evolution (Facchini et al. 2016;Winter et al. 2018a), and evidence of it occurring has been found in σ Ori (Maucó et al. 2016;Ansdell et al. 2017), where a clear correlation between disk mass and distance from the central Herbig O9V star has been observed and in the Orion Nebula Cluster (Mann & Williams 2010;Eisner et al. 2018). However, in CrA no correlation between the mass of the disks (or the disk detection rate) and the distance from the brightest star (R CrA) is found. Moreover, in σ Ori external photo-evaporation has been shown to affect disks up to 2 pc away from the Herbig star, where the geometrically diluted far-ultraviolet (FUV) flux reached a value of ∼ 2000 G 0 . The spectral type of R CrA is still uncertain, ranging from F5 (e.g. Garcia Lopez et al. 2006) to B8 (e.g. Hamaguchi et al. 2005). Even in the latter case, assuming a typical FUV luminosity for a B8 star of L FUV ∼ 10 L (Antonellini et al. 2015) and accounting for geometric dilution, we find that the FUV flux would drop to ∼ 1 G 0 in the first inner pc from R CrA, thus ruling-out external photo-evaporation as an explanation. Also, this calculation neglects dust absorption, which is probably very effective in the Coronet cluster around R CrA.
Because of the M disk − M relation presented in Sec. 4.5, it is also possible that a system dominated by low-mass stars shows a low-mass disk population, regardless of its age, as in the case of IC348 (Ruíz-Rodríguez et al. 2018). It is therefore important, when comparing disk dust masses from different regions, to verify that they have the same stellar mass distribution. In order to do this, we employ a Monte Carlo (MC) approach similar to that used by Andrews et al. (2013). We first normalize the stellar populations by defining stellar mass bins and randomly drawing the same number of sources in each bin from the reference sample (CrA) and from a comparison sample (Lupus, Chamaeleon I or Upper Sco). We then perform a two-sample logrank test for censored datasets between the disk dust masses of the two samples, to test the probability (p φ value) that the two samples are randomly drawn from the same parent population. A low p φ value  Andrews et al. (2013). p φ is the probability that the synthetic population drawn from the comparison sample (Lupus and Upper Sco) and the reference sample come from the same parent population. f (< p φ ) is the cumulative distribution for p φ resulting from the logrank two-sample test for censored datasets after 10 4 MC iterations.
indicates that the difference in disk masses cannot only be ascribed to different stellar populations and that some other factor, such as disk evolution and the age of the system, must play a role. This process is repeated 10 4 times, and the results are used to create the cumulative distributions shown in Fig. 7. When using Upper Sco as a comparison sample, we find a median p φ value of 0.53, while the median p φ for Lupus is only 0.004. The conclusion is that even when accounting for the M dust −M relation, the disk dust mass distribution of CrA appears to be statistically different from that of Lupus, while it is significantly more similar to that from that of Upper Sco. Therefore, the comparably low masses of the protoplanetary disks in CrA cannot be explained in terms of the low stellar masses.
Another way a disk can lose part of its mass is via tidal interaction with other stars (e.g. Clarke & Pringle 1993;Pfalzner et al. 2005). This mechanism is, however, only effective in much denser environments than CrA (e.g. Winter et al. 2018b). In principle, it is possible to imagine that at very early stages most of the stars were located in a dense region (e.g. the Coronet) where they interacted violently before being ejected. However, the very low velocity dispersion of the stars in the cluster makes this scenario very unlikely (Neuhäuser et al. 2000). Tidal interaction can be effective in removing dust mass from a disk even in later stages when the disk is in a binary system (e.g. Artymowicz & Lubow 1994), as proposed to explain the low mm flux of some objects in Taurus by Long et al. (2018). A higher than usual binary fraction could therefore explain the low disk masses observed in CrA. However, Ghez et al. (1997) show that the binary fraction of CrA is indistinguishable from those of Lupus and Chameleon I.
Finally it is possible that the low mass distribution observed today is a consequence of a population of disks that has formed with a low mass from the very beginning. For example, the disk formation efficiency in a cloud with mass M 0 depends on the sound speed c s and on the solid body rotation rate Ω 0 , where we have defined the disk formation efficiency as the fraction of M 0 that is in the disk at the end of the collapse stage, or as the ratio between M disk /M at that time (Cassen & Moosman 1981;Tere-bey et al. 1984). In particular, clouds with higher c s and Ω 0 (i.e. warmer or more turbulent) will form more massive disks (also see Appendix A in Visser et al. 2009). Therefore, a cold parent cloud or one with low intrinsic angular momentum Ω 0 , will form disks with a lower mass, and with a lower M disk /M as observed in CrA. Consistently, observations of dense cloud cores in the CrA cloud show line-widths lower than in other regions (Tachihara et al. 2002). Moreover, because of the smaller circularization radius, the formed disks would alse be smaller (e.g. Dullemond et al. 2006) and potentially mostly optically thick, thus hiding an even larger fraction of the mass . Alternatively, small and optically thick disks could result from magnetic braking of the disks by means of the magnetic field threading the disk and the surrounding molecular cloud at the formation stage (e.g. Mellon & Li 2008;Herczeg & Hillenbrand 2014;Krumholz et al. 2013). The same scenario was proposed by Maury et al. (2019) to explain the low occurrence of large (> 60 AU) Class 0 disks in the CALYPSO sample.
Such scenarios, although not testable with the present dataset, are consistent with the low disk mass distribution and with the low intercept of the M disk − M in CrA and are not in contradiction with the young age of the stellar popultion. If the parent cloud initial conditions are indeed responsible for the low masses observed, this would be an additional critical aspect to be considered when studying planet formation and evolution. Since the conditions at the epoch of disk formation can be different in each star-forming region, proper modelling is required to assert to which extent they can affect the initial disk mass distribution, the subsequent disk evolution, planet formation and planetary populations.
Observationally, this could be tested by observing the mass of disks around Class 0 and Class I objects in CrA: if the disks are born with a low-mass, the disk mass distribution even at these younger stages should be significantly lower than in other regions.

Conclusion
We presented the first ALMA survey of 43 Class II protoplanetary disks in the Corona Australis nearby (d = 160 pc) star forming region, in order to measure their dust content and understand how it scales with the stellar properties. The ultimate goal was to test if the relations between disk properties, age of the stellar population found in other surveys also hold for this region.
1. The average mm fluxes from the disks in CrA is low. This in turn converts into a low disk mass distribution. Even though our observations are able to constrain dust masses down to ∼ 0.2 M ⊕ , the detection rate is only 56%. Moreover, we find that only 3 disks in our sample have a dust mass ≥ 10 M ⊕ and thus sufficient mass to form giant planet cores. 2. We obtained VLT/X-Shooter spectra for 8 objects with previously unknown spectral type, and derived their stellar physical properties. 3. Despite the apparent young age of the CrA stellar population, we find that the dust mass distribution of the disks in CrA is much lower than that of the Lupus young star forming region which shares a similar age, while it appears to be consistent with that in the 5-10 Myr old Upper Sco association. The correlation between disk dust mass M dust and stellar mass M previously identified in all other surveyed star forming regions is confirmed. However, because of the low mass of the disks in our sample we find a much lower intercept. The large scatter of the data points does not allow the slope of the relation to be well constrained for CrA.
4. Since most of the age estimates of the CrA regions are based on the population of the compact Coronet cluster, a possible explanation for the low disk masses might be in principle that CrA also hosts an old population of disks, consistently with previous observations. The position of the objects of our sample on the HR diagram, however, seems to support the idea of a mostly coeval, young population. 5. Low disk masses in a young star forming region can be explained by external photo-evaporation (as in the case of σ Ori) or by a low stellar mass population (as in IC348). With our analysis, we can rule out both these scenarios for CrA. Tidal interaction between different members of CrA, stripping material from the disks, as well as close binaries can also be ruled out. 6. We suggest that initial conditions may play a crucial role in setting the initial disk mass distribution and its subsequent evolution. Small disks with low mass can originate from a cloud with very low turbulence or sound speed, or can alternatively result from disk magnetic braking. It is therefore important to better study the impact of initial conditions on the disk properties, especially if planet formation occurs even before 1 Myr age, as the recent results from Tychoniec et al. (2018) and Manara et al. (2018) suggest.
Future surveys including younger Class 0 and I objects in CrA and other star forming regions will help testing wether or not initial conditions play a critical role in shaping the physical properties of circumstellar disks.   Uncertain estimate of SpT due to the low S/N of the spectra. Table B.2: Night log and basic information on the spectra. In Column 1 is the name of the source, in Column 2 the date and time of the observations, in Column 3-5 the exposure times, in Column 6-8 the slit widths, in Column 9-11 the S/N measured at the indicated wavelengths, in Column 12-13 we show whether or not the H α and Li lines have been detected.