© P. Pinilla et al. 2021

1 Introduction

The current exoplanet population shows a large diversity of properties, which may have an equivalent in the properties and initial conditions of the parental protoplanetary disks. In the context of exoplanets, some trends have already emerged with stellar mass, such as giant planets being more frequent around more massive and more metal-rich stars, but sub-Neptunes being more frequent around lower mass stars (e.g., Santos et al. 2000; Johnson et al. 2010; Mulders et al. 2015). In the context of protoplanetary disks, surveys at millimeter wavelengths have shown that dust disk masses also increase with stellar mass (e.g., Andrews et al. 2013; Pascucci et al. 2016; Ansdell et al. 2017). These two observational results are consistent with the predictions of core-accretion models for the formation of planets, in which giant planet formation is relatively inefficient in disks with low mass around very low-mass stars (e.g., Payne & Lodato 2007; Liu et al. 2020). Alternatively, giant planets around M-dwarfs may form when the disk is gravitationally unstable (Mercer & Stamatellos 2020).

Recent observations of protoplanetary disks at high angular resolution and sensitivity, in particular with the Atacama Large Millimeter/submillimeter Array (ALMA), revealed that nearly all large disks (≳50 au in dust continuum emission) host substructures (Andrews et al. 2018; Long et al. 2018), with rings and gaps being the most common type. If these substructures are created by planets, most of these planets have masses and orbital radii that are not detectable by any of the methods currently used to detect exoplanets (Bae et al. 2018; Lodato et al. 2019).

However, most of the current ALMA programs attempt to study the prevalence and diversity of these structures around T Tauri and Herbig AeBe stars with different properties; little is known about the existence of substructures in disks around very low-mass stars (VLMSs, ≲0.1–0.2 M_⊙) or brown dwarfs (BDs) because observations of these disks are challenging as they are fainter and smaller (e.g., Rilinger et al. 2019; Sanchis et al. 2020). Recently, Kurtovic et al. (2021) presented a sample of six disks around VLMSs in the Taurus star forming region observed at high resolution with ALMA, demonstrating that these disks follow the expected relations for disk dust masses versus stellar mass ($M_{\text{dust}}\propto M_{\star }^{1.1}$) and disk dust radii versus disk luminosity ($R_{\text{dust}}\propto L_{\text{mm}}^{0.54}$) as seen in their counterparts (more massive and larger disks, e.g., Andrews et al. 2013; Pascucci et al. 2016; Hendler et al. 2020). These recent observations are also unveiling that substructures in the low-mass star regime could be as common as in T-Tauri and Herbig disks (50% of the sample revealed substructures at a resolution of 0.1″, although the sample was selected to target the brightest disks around VLMSs). These observations suggest that structures may be universal and they question the origin of the observed structures; for example, whether or not these rings and gaps may be formed by (giant) planets independent of stellar mass.

The spectral index, which is the slope of the spectral energy distribution (SED) at the millimeter wavelength, can be interpreted in terms of grains size with low values (≲3.0) corresponding to large grains (e.g., Draine 2006). Disks around VLMSs and BDs are known to host millimeter (mm)-sized particles, as evidenced by their low millimeter spectral indices (e.g., Ricci et al. 2014; Pinilla et al. 2017). These observations challenge current models of grain growth in protoplanetary disks, where radial drift velocities of dust particles towards the star are higher in VLMSs and BDs than around higher mass young stars (Pinilla et al. 2013; Zhu et al. 2018). To retain enough particles in VLMS disks and explain low spectral index values, multiple pressure bumps of high amplitude are required. This would translate into multiple rings in millimeter observations with infinite angular resolution. Vortices can also trap particles (e.g., Barge & Sommeria 1995; Klahr & Bodenheimer 2006; Varnière & Tagger 2006; Ataiee et al. 2013), in which case observationswould reveal one or more nonaxisymmetric features. Nonetheless, from the observational point of view, our current understanding of substructures around VLMSs is very limited, with only two clear cavities detected in CIDA 1 (Pinilla et al. 2018a) and MHO 6 (Kurtovic et al. 2021) and a hint of rings and gaps in the disk around 2MASS04334465 + 2615005 (Kurtovic et al. 2021).

CIDA 1 (2MASS J04141760 + 2806096) is a young VLMS in the Taurus star forming region that is close to the hydrogen burning limit. It is known to be surrounded by a very luminous (for its mass) disk at millimeter emission. For this reason it was selected as a candidate to measure the spectral index and constrain the presence of large dust grains. ALMA Cycle 0 observations at 0.89 and 3.4 mm revealed a shallow integrated spectral index (~ 2) and relatively large disk (Ricci et al. 2014), suggesting the presence in the disk of a substantial population of millimeter-sized grains. ALMA Cycle 3 observations of CIDA 1 revealed for the first time a dust-depleted inner cavity around a VLMS, with a ring with maximum brightness at ~20 au (Pinilla et al. 2018a), which is a typical feature of the transition disk (TD) population. The observed structure in CIDA 1 may suggest that cavity opening processes are common in millimeter-bright disks independent of the central star mass, including VLMSs (Pinilla et al. 2018b, 2020; Sinclair et al. 2020).

The mechanisms that are most likely to form the cavity in the CIDA 1 disk include the interaction of the disk with a massive embedded planet. However, core-accretion models predict that it is very challenging to create the required pressure bump with an embedded planet under the physical conditions of CIDA 1, because a large fraction of the disk mass (5–45%, where higher values correspond to a disk with high viscosity) is needed to form such a massive planet in the first place. One alternative is that the initial disk mass is high enough such that the disk is gravitationally unstable and capable of forming massive planets (Mercer & Stamatellos 2020). In the case of CIDA 1, a Saturn-mass planet is the minimum mass required to open a gap (with low viscosity) and to trap millimeter-sized particles, but the predicted integrated spectral indices for this scenario are too high compared to observations (Pinilla et al. 2017). This contrasts with T-Tauri and Herbig disks, where a single massive planet is capable of opening the cavity, and trapping and retaining enough mm-grains to explain low spectral indices (Pinilla et al. 2014). As a consequence, it is expected that strong and multiple pressure bumps are present in VLMS disks with low spectral indices.

Nevertheless, in previous ALMA observations the cavity of CIDA1 was marginally resolved and the nature of this structure was unclear. It therefore remains unclear as to whether or not the cavity is empty of millimeter-sized particles, whether or not the ring is composed of multiple structures, and whether or not there are variations of the spectral index throughout the disk that indicate grain growth in particular regions. In this paper, we present observations with ALMA of CIDA 1 in Band 7 (0.9 mm) and Band 4 (2.1 mm) with a high angular resolution at the two wavelengths (~ 0.05″ × 0.034″) in order to potentially resolve small-scale features and resolve possible variations of the spectral index throughout the disk. In addition, X-shooter observations on the Very Large Telescope (VLT) were acquired to get better constraints of the CIDA 1 accretion rate, which can help to test models of cavity formation, such as photoevaportation (e.g., Ercolano & Pascucci 2017). This paper is organized as follows. Section 2 summarizes the observations and data reduction. Section 3 presents the results of our X-shooter and ALMA observations, and the analysis of the dust morphology in the visibility plane. Section 4 discusses our observations in the context of dust evolution models and planet–disk interaction. Finally, in Sect. 5 we discuss the results of our observations, and in Sect. 6 we present our conclusions.

2 Observations

2.1 ALMA observations

The data sets studied in this work include ALMA observations at high angular resolution (≈ 0.05 × 0.034″) of CIDA1 at 0.9 and 2.1 mm wavelengths, under the ALMA project 2018.1.00536.S (PI: A. Natta). The 0.9 mm data were obtained by ALMA Band 7, and the correlator was configured to observe four spectral windows: two covered dust continuum emission centered at 344.0 and 333.25 GHz, and two centered at the molecular lines ¹²CO (J : 3−2) and ¹³CO (J : 3−2). The frequency resolution of the channels is 976.562 kHz for continuum and ¹³CO spectral windows, while it was 244.141 kHz for ¹²CO (approximately 1 and 0.2 km s⁻¹, respectively). For the imaging and analysis, we combined all available ALMA Band 7 observations, thus including the archival data from 2015.1.00934.S (PI: L. Ricci) and 2016.1.01511.S (PI: J. Patience) already published in Pinilla et al. (2018a) and Kurtovic et al. (2021).

Observations at 2.1 mm were carried out with ALMA Band 4. The correlator was configured to observe four spectral windows: three were located to observe dust continuum emission, centered at 134.4, 136.2, and 148.2 GHz, while the remaining one was centered at the molecular line CS J:3−2. The frequency resolution of the continuum channels was 976.562 kHz, while the CS line was observed at 31 kHz (approximately 2.4 and 0.1 km s⁻¹, respectively). The details of all the data used in this paper are summarized in Table 1.

After ALMA standard pipeline calibration, a series of preparation steps have to be applied to combine and self-calibrate the datasets. Using CASA 5.6.2, we extract the dust continuum emission by flagging the channels located at ± 25 km s⁻¹ from each targeted spectral line. The remaining channels from all spectral windows are averaged into 125 MHz channels. To avoid differences in the weight calculation through different ALMA cycles, we apply the task statwt to recalculate the visibility weighting according to the observed scatter. Next, we aligned the datasets and compared the flux calibration of each individual execution to ensure they were consistent. In Band 7, we found and corrected a 10% discrepancy between the observation from project 2015.1.00934.S and all the other observations, and so we rescaled this compact configuration data to match the others. In Band 4, we corrected a discrepancy of 19% between the first observation of the extended baselines and all others.

To enhance the signal-to-noise ratio (S/N), self-calibration was performed in two stages. First, we self-calibrated the data sets we classified as having “compact” baseline configuration in Table 1, and then we combined these with the “extended” baselines and self-calibrated again. We used a Briggs robust parameter of 0.5 for the imaging of the self-calibration process. The combined data sets from Band 4 had a peak S/N of ≈ 23 before self-calibration, thus we only applied one phase and one amplitude calibration, using the whole integration time as the solution interval.

In order to reduce the data volume for the visibility analysis, we averaged the continuum emission into channels of 250 MHz in width and 30 s time-binning. We used the central frequency of each binned channel to convert the uv-coordinates into wavelength units.

All the dust continuum calibration steps, including the centroid shifting, flux calibration, and self-calibration tables, were then applied to the molecular line emission channels. The continuum emission was subtracted using the uvcontsub task. We generate images of 0.5 km s⁻¹ in velocity width, with a robust parameter of 1. To enhance the S/N, we chose to apply a uv-tapering with a Gaussian of 0.03″ on the ¹²CO, and 0.04″ on the ¹³CO. We used thepackage bettermoments (Teague & Foreman-Mackey 2018) to create the kinematic map. This package fits a quadratic function to find the peak intensity of the line emission in each pixel, and the velocity associated to it.

2.2 X-shooter observations

We acquired one broad-band flux-calibrated spectrum of CIDA 1 with the X-shooter instrument on the VLT. The observations were taken for Pr.Id. 105.2061.001 (PI: A. Natta), on the night of October 30, 2020. Sky conditions were clear (CLR), seeing at zenith was ~1′′. The spectrum was observed using slit width 1.0′′ –0.4′′–0.4′′ in the three arms (UVB, VIS, NIR) of the X-shooter spectrograph, respectively. These spectra were reduced using the X-shooter pipeline (Modigliani et al. 2010) v3.5.0 run through the ESOReflex tool. This pipeline performs the main reduction steps, including flat fielding, correction for bias and dark, wavelength calibration, rectification, and extraction of the 1D spectrum. The latter is then corrected for telluric absorption lines using the molecfit tool, version 3.0.3 (Smette et al. 2015). The absolute flux calibration of the spectra is then performed by scaling the spectra to the one obtained in the same observing block using a wider set of 5.0′′ slits. The flux rescaling factors are 1.2, 3.7, and 3.5 in the three arms, in line with the fact that the slit widths in the three arms are smaller than the seeing at the corresponding wavelengths.

3 Results from ALMA and X-shooter observations

3.1 Images, dust disk mass, and spectral index

Figure 1 shows the final images of CIDA 1 of the ALMA observations in Band 7 (0.9 mm) and in Band4 (2.1 mm) and the details of the images are summarized in Table 2. The fluxes from this table are obtained from the images assuming all the emission within a circular area of 0.5″ from the center. The images of the ¹²CO and ¹³CO are shown in Fig. 2; these were obtained by processing the channel maps with the Python package bettermoments (Teague & Foreman-Mackey 2018), and clipping the emission at 3σ. The emission of CS was not detected in Band 4, although an increase in the noise level is observed close to the central velocity channel. We imaged the channels around the expected location of the CS line with a velocity resolution of 1.0 km s⁻¹. Three moment 0 images were generated from these channel maps, spanning the velocities from 1–4, 5–8, and 9–12 km s⁻¹. When the rms is measured in a circle of 1′′ centered on CIDA 1, the moment 0 image spanning 5–8 km s⁻¹ has a rms of 1.56 mJy beam⁻¹, while the other two moment 0 images have an rms of 1.04 mJy beam⁻¹.

The dust continuum images at both wavelengths reveal unresolved emission from the inner disk, a large gap,and a ring-like structure peaking at around 20 au. We use the distance to CIDA 1 of 137.5 ± 1.1 pc derived by Bailer-Jones et al. (2018) using Gaia DR2 (Gaia Collaboration 2018). Figure 3 shows in linear and logarithmic scale the azimuthally averaged radial intensity profiles of the deprojected continuum images at each wavelength (using the inclination and position angle report in Table 3). For this figure, each profile is normalized to the peak of the ring and the images are previously restored with the same beam before deprojection (0.05′′ × 0.034′′). For the deprojection, we use the inclination (inc) and position angle (PA) from Table 3, which are obtained from the visibility analysis explained in Sect. 3.2. The error bars include the standard deviation of each elliptical bin divided by the square root of the number of beams spanning full azimuthal angle at each radial bin. In this figure, it is possible to see that emission at both wavelengths peaks at the ring location (20 au) and there is a very bright inner disk. The ¹²CO peak intensity map also shows a cavity that remains unresolved with a size of ~ 7 au along the major axis of the disk. Another interesting result from the ¹²CO peak intensity map is that this emission is more extended than the dust continuum emission, as found in Kurtovic et al. (2021). The analysis of the ¹²CO and ¹³CO will be presented in a future paper by Curone et al. (in prep.).

We calculate the optical depth of the peak of the continuum ring, assuming τ = −ln[1 − T_brightness∕T_physical], with T_brightness and T_physical being the brightness and physical temperature respectively. We obtain the brightness temperature from the full black-body Planck function. We used a radially constant physical temperature of 20 K and we find that the optical depth at the peak is τ_peak,B7 = 0.48 and τ_peak,B4 = 0.28 for Band 7 and Band 4, respectively. We also calculate τ assuming a luminosity-dependent relation for the temperature from van der Plas et al. (2016), that is, $T=22{(L_{\star }/L_{\odot })}^{0.16}$ K, in which case we obtain a physical temperature of 14.7 K and τ_peak,B7 = 0.74 and τ_peak,B4 = 0.41. Using our radiative transfer calculations presented in Sect. 4 for CIDA 1, the temperature expected at the midplane and at the peak location (~20 au) is around 30 K, in which case τ_peak,B7 = 0.29 and τ_peak,B4 = 0.18. We therefore assume that most of the emission is optically thin to calculate the dust disk mass. These values for the optical depths fall in the same range as reported by the DSHARP sample (Dullemond et al. 2018) for more massive disks. One potential explanation for this range of the optical depth inside the ring is planetesimal formation in this region (Stammler et al. 2019).

Under the assumption of optically thin emission, we calculate the dust disk mass as $M_{\text{dust}}\simeq \frac{d^2{F}_{\nu }}{{\kappa }_{\nu }{B}_{\nu }(T)}$, where d is the distance to the source, F_ν is the total flux at a given wavelength, and B_ν is the blackbody surface brightness at a given temperature (Hildebrand 1983). Taking a mass absorption coefficient (κ_ν) at a given frequency as κ_ν = 2.3 cm² g${}^{-1}\times {(\nu /230\text{\hspace{0.17em}GHz})}^{0.4}$ (which is typical for a population of large millimeter grains, Beckwith et al. 1990; Andrews et al. 2013, which is compatible within 2σ with the spectral index we infer for CIDA 1 as shown later), and a dust temperature of 20.0 K or 14.7 K, the dust disk mass is 7.4 M_Earth or 12.0 M_Earth, respectively,when taking the total flux from the image in Band 7. For Band 4, the dust disk mass is 7.8 M_Earth or 11.4 M_Earth for 20 or 14.7 K, respectively. Therefore, the inferred dust disk mass is consistent between the two wavelengths with a large uncertainty that depends on the assumed temperature. In addition, κ_ν is a very unknown quantity in protoplanetary disks, and for direct comparison we take the same κ_ν assumed to calculate the dust disk mass of large populations of protoplanetary disks (e.g., Pascucci et al. 2016). Recent work from Lin et al. (2020) about inferring (sub)millimeter dust opacities from resolved multi-wavelength continuum observations of edge-on disks provides support to this widely used opacity assumption from Beckwith et al. (1990).

With the total flux from the images at the two wavelengths, we obtained the spatially integrated spectral index as α_mm = log(F_B7∕F_B4)∕log(ν₇∕ν₄) = 2.0 ± 0.2, where the uncertainty includes the rms of the observations and 10% uncertainty from the flux calibration. This value is in agreement with the values found in Ricci et al. (2014). The very low fluxes within the gap do not allow us to obtain reliable values of the spectral index in this region (see right panel of Fig. 3). When taking the emission either from a circled area that is centered at the location of CIDA 1 and that is as big as the beam or an area that encircles the ring, the integrated spectral index in each region is around 2.0. The right panel of Fig. 3 shows the radial profile of the spectral index, which shows a slight decrease inside the ring, but the differences are not significant (<1σ). Within the uncertainties, the spectral index seems to be relatively constant throughout the disk, as found by Tazzari et al. (2020) in a sample of disks in the Lupus star formation region.

Fig. 1 ALMA Band 7 (0.9 mm) and Band 4 (2.1 mm) continuum images of the CIDA 1 disk. The synthesized beam is shown in the lower left corner of each image. The details of the images are summarized in Table 2.

Fig. 2 Moment 0 and moment 1 of the 12CO and 13CO. The upper right bar corresponds to 10 au size. The dashed circles (projected with the same inclination and PA as CIDA 1, Table 3) are 1.0″ and 0.6″ in size from the center for the 12CO and 13CO, respectively. The details of the images are summarized in Table 2.

Fig. 3 Azimuthally averaged radial intensity profiles of the deprojected continuum images (lineal scale on the left and logarithmic scale in the center) at each wavelength. Each profile is normalized to the peak. The profiles are obtained from images that are restored with the same beam (50 mas × 34 mas). The error bars include the standard deviation of each elliptical bin divided by the square root of the number of beams spanning the full azimuthal angle at each radial bin. The right panel shows the radial profile of the spectral index taking into account the errors on the radial intensity profiles in addition to 10% of the error due to flux calibration.

3.2 Dust morphology from visibility fitting

We describe the emission on both ALMA bands with axisymmetric parametric models, and for comparison we also use the Gaussian Processes modeler Python package frank (Jennings et al. 2020). While a single centrally peaked Gaussian is enough to describe the inner emission of CIDA 1, for the ring emission we tried several different combinations of Gaussian profiles and/or radially broken Gaussian profiles (Gaussians with different widths to each side of the peak). After these tests, we find that the best parametric models for the ALMA Band 4 (2.1 mm) and Band 7 (0.9 mm) emission are given by: $$I_{\text{2}.\text{1\hspace{0.05em}mm}}(r)={10}^{f_0}g(r,{\sigma }_0)+{10}^{f_1}{g}_{\text{bg}}(r-r_1,{\sigma }_{1i},{\sigma }_{1o})\text{,}$$(1) $$I_{\text{0}.\text{9\hspace{0.05em}mm}}(r)=\text{\hspace{0.17em} \&}{10}^{f_0}g(r,{\sigma }_0)+{10}^{f_1}{g}_{\text{bg}}(r-r_1,{\sigma }_{1i},{\sigma }_{1o})\text{nn \&}+{10}^{f_2}g(r,{\sigma }_2)$$(2)

where g is a Gaussian profile, and g_bg is a broken Gaussian, described by: $$g(x,\sigma )=\exp \left(-\frac{x^2}{2{\sigma }^2}\right)\text{,}$$(3) $$g_{\text{bg}}(x,\sigma )=g(x,{\sigma }_i)\text{\hspace{0.17em}for }x\le 0.\text{}g(x,{\sigma }_o)\text{\hspace{0.17em}for }x>0..$$(4)

For each model, the visibilities are obtained by optimizing the model profile with a spatial offset (δ_RA, δ_Dec), an inclination and PA, which are used to deproject the observational data. Therefore, each model has four extra free parameters in addition to those that describe the intensity profile. The Fourier Transforms to obtain the models visibilities and the χ² calculation are computed with the galario Python package (Tazzari et al. 2018) using a pixel size of 0.2 mas.

We sample the posterior probability distribution with a Markov chain Monte Carlo (MCMC) routine based on the emcee Python package (Foreman-Mackey et al. 2013). We used a flat prior probability distribution over a wide parameter range, such that the walkers would only be initially restricted by geometric considerations (inc ∈ [0, 90], PA ∈ [0, 180], σ ≥0). We ran more than 250 000 steps after convergence to find the most likely parameter set for each model, as well as taking the 16th and 84th percentile for the error bars. Our results are shown in Table 3.

Figure 4 shows the real part of the visibilities after centering and deprojecting the data in Band 7 and in Band 4 vs. the best fit model that minimizes χ² from galario of the continuum data. In addition, this figure includes the intensity profile from the best fit in Band 7 and Band 4 (orange and blue line, respectively), with the 5000 randomly selected chains after convergence overlaid, and the right panel shows the spectral index radial profile obtained from the intensity profiles. The intensity profiles from the best fit reveal that the ring in Band 7 seems to have a shoulder after the peak of emission that is not visible in the Band 4 data. It is worth noting that a fit in Band 7 without this additional ring in the outer disk (e.g., the model used for Band 4) produced significant residuals, suggesting the presence of a further ring. In addition, our best fits show that the inner disk is as bright as or brighter than the peak brightness of the ring. The gap seems to be shallower in Band 4, but this result needs to be taken with care because the S/N of Band 4 is around half that of Band 7, and to confirm these differences between Band 7 and Band 4, deeper observations are required in Band 4. The residuals of the Band 7 fit hint at a nonaxisymmetric structure of the ring, and deeper observations at this wavelength are also needed to check for this potential asymmetry.

From the best intensity profiles, we calculate the radial profile of the spectral index presented in the right panel of Fig. 4, which shows the same trend as those seen in the images (right panel of Fig. 3), where there is a hint of a slight decrease inside the ring and overall a low value of the spectral index in the entire disk (~2) with a high decrease in the gap. However, the spectral index there is not trustworthy because of very low fluxes and large uncertainties as shown in the right panel of Fig. 3. We also check if there is any evidence from the visibilities that the spectral index varies by taking the ratio of the visibilities, and this ratio is constant within the uncertainties.

Figure 5 shows the observations versus model images from our best fit with galario (before and after convolution)and the residuals, for both Band 7 and Band 4. Our residuals do not show any clear sign of additional substructures, either radialor azimuthal at the current resolution and sensitivity of the observations.

Fig. 4 From left to right: real part of the visibilities after centering and deprojecting the data in Band 7 and in Band 4 vs. the best-fit model from galario of the continuum data. Intensity profile from the best fit in Band 7 and Band 4 (orange and blue line, respectively), with 5000 randomly selected chains after convergence overlaid. Right panel: spectral index radial profile obtained from the intensity profiles.

3.3 Stellar properties from the X-shooter spectrum

The spectrum of CIDA 1 is rich in permitted emission lines and shows a prominent Balmer jump. We fitted the spectrum following the procedure described by Manara et al. (2013). Briefly, we explore a parameter space including several photospheric templates of young stars, a set of slab models to reproduce the accretion emission continuum spectra, and increasing extinction value until the sum of the photospheric template and the accretion spectrum better reproduce the de-reddened spectrum of CIDA 1. In this case, the solution is highly degenerate between templates with spectral types between M4.5 and M6.5, and extinction values A_V ~ 2−4 mag. We therefore use an independent method to constrain the best values for A_V, and we measure the emission line fluxes for a number of permitted emission lines in the whole X-shooter spectrum. We then convert the measured line luminosity in accretion luminosity using the relations by Alcalá et al. (2017), and repeat this exercise by increasing the amount of A_V from 0 to 4 mag in steps of 0.25 mag. The accretion luminosity derived from the lines shows the minimum dispersion with A_V ~ 3.5 mag (see Fig. 6).

Constraining the values of A_V between 3 and 3.8 mag, we obtain a best fit from the X-shooter spectrum following the evolutionary tracks from Baraffe et al. (2015), with SpT M5, A_V = 3.5 mag, resulting in a value of stellar luminosity L_⋆ = 0.15 ± 0.03 L_⊙, accretion luminosity log (L_acc∕L_⊙) = −1.34 ± 0.3, a stellar mass M_⋆ = 0.1–0.2 M_⊙, and a mass accretion rate Ṁ_acc = 1.4 × 10⁻⁸ M_⊙ yr⁻¹. The best fit is shown in Fig. 7. This value of mass accretion rate is in line with the results of Herczeg & Hillenbrand (2008), and is higher than the one we reported in Pinilla et al. (2018a). The addition of the infrared emission lines in the analysis allows us to better constrain the extinction value from the emission lines.

We note that CIDA 1, despite being a strong accretor, shows very little emission in forbidden lines, with only the [OI]λ6300 Å line detected in the spectrum. The [OI] emission may be a tracer of some winds in the inner part of the disk (Simon et al. 2016; McGinnis et al. 2018; Banzatti et al. 2019). These authors find that [OI] lines in transition disks are narrow, like the one we see in CIDA 1, but it is unclear in this case whether it originates from photoevaporation (e.g., Ercolano & Owen 2016) or a magneto-hydrodynamical (MHD) wind (e.g., Banzatti et al. 2019).

Fig. 5 Observations vs. model images from our best fit with galario (before and after convolution) and the residuals. Top panels correspond to Band 7 and the bottom panels correspond to Band 4.

Fig. 6 Accretion luminosity derived from the luminosity of various emission lines, as labeled, using the relations by Alcalá et al. (2017), and assuming AV = 3.5 mag. The blue line is the best-fit accretion luminosity obtained fitting the whole X-shooter spectrum of CIDA 1. The red line is the mean value of Lacc obtained from the luminosity of the emission lines.

Fig. 7 Fit of the spectrum of CIDA 1 obtained here. The dereddened spectrum of the target (red) is reproduced by the sum of a photospheric template (green) and a slab model to reproduce the accretion shock emission (cyan).

4 Comparison with dust evolution models

In this section, we compare our observations with dust evolution models in combination with radiative transfer calculations. Our models are in the context of giant planets opening a gap in the disk and potentially creating a cavity to explain the observed properties of transition disks.

A pressure bump at the outer edge of a planetary gap can trap millimeter-sized particles and can create a dust-depleted cavity as seen at millimeter wavelength depending on the mass of the planet. The planet mass also determines whether or not an inner disk survives after million-year timescales. In a subset of transition disks, near-infrared (NIR) excess has been detected, which is interpreted as an optically thick dusty inner disk located within the first astronomical units (e.g., Espaillat et al. 2010). Pinilla et al. (2016) showed that this NIR excess emission can be explained at any time of evolution when the embedded planet is not more massive than ~ 1−5 M_Jup around a solar-type star (with a disk viscosity of 10⁻³). The gap formed by such a planet allows dust particles to pass throughout the gap and refill the inner disk with dust particles. In CIDA 1, there is a very bright inner emission in both Band 7 and Band 4, which can help to constrain the mass and position of the potential planet creating the large gap.

Under the assumption that it is one single planet creating the large gap, Pinilla et al. (2018a) analytically calculated the minimum mass of the planet needed to open a gap in the gas surface density under the conditions of CIDA 1 and create a ring at 20 au. Our current observations show that the minimum of this gap is around 8 au. Assuming this location for the position of the potential planet, the minimum mass of the planet to open a gap in a disk with α_visc = 10⁻³ is ~ 1.2 M_Saturn (or a planet-to-star mass ratio of 3.2 × 10⁻³). When assuming such a planet at 8 au, the pressure maximum is located at ~ 13 au. Hence either the planet must be more massive to create a wider gap and move the pressure maximum further out, or the planet should be located further out. In the latter case, the minimum of the gap seen in millimeter emission would not correspond to the location of the planet. This is a typical assumption in most hydrodynamical simulations that aim to explain some observational properties of protoplanetary disks, where it is assumed that the minimum of the dust emission co-locates with the minimum of the gas surface density, i.e., the local minimum of a gap opened by a planet.

When forcing the planet to be located at 8 au, a planet-to-star mass ratio of 2.3 × 10⁻² is needed to have the pressure maximum at ~20 au, where the peak of the emission is located. This corresponds to a planet mass of ~ 2.4M_Jup for CIDA 1. Such a planet is unlikely to form by core accretion around such a low-mass star and in a disk of low mass (Liu et al. 2020). As an experiment, we performed dust evolution simulations assuming a gap created by a ~ 2.4M_Jup located at 8 au. The shape of the gap is obtained analytically using the minimum mass to open a gap criterion by Crida et al. (2006) and following Pinilla et al. (2015). For this simulation, we assumed an initial dust-to-gas ratio of 1/100 where all the grains are initially micron-sized particles. We used the dust evolution package DustPy (Stammer & Birnstiel, in prep.). The simulations include the growth, fragmentation, and erosion of particles. The fragmentation velocity of the particles is 10 m s⁻¹ in the entire disk. We assumed α_visc = 10⁻³ and a disk mass of 5.5 M_Jup around a stellar object of 0.1 M_⊙ mass (i.e., M_disk∕M_CIDA1 ~ 0.05). For more details of the dust evolution models see Birnstiel et al. (2010).

Figure 8 shows the dust density distribution as a function of grain size and distance from the star for the case of a ~ 2.4M_Jup located at 8 au and three different times of evolution (0.1, 0.5, and 1 Myr). The white line corresponds to the particle size for which the Stokes number (St) is equal to unity ($\text{St\hspace{0.17em}=\hspace{0.17em}}\frac{\text{a}{\rho }_{\text{s}}}{{\Sigma }_{\text{g}}}\frac{\pi }{\text{2}}$), and therefore it also represents the shape of the gas surface density. These results show how the inner disk is depleted, because in this case the gap formed by such a massive planet is so deep that it blocks all dust particles, and thus there is no replenishment of dust particles from the outer part to the inner part of the disk. As a result, an empty cavity is created, and already at 0.1 Myr of evolution almost all the grains that were initially within the inner edge of the gap have grown and drifted towards the star.

This experiment demonstrates that such a massive planet is unlikely to exist in this disk, not only because its formation challenges current models of planet formation, but also because the new observational evidence of the inner disk at millimeter wavelengths excludes this possibility.

As a potential solution, we performed an experiment moving the location of the planet to 12 au, which allows us to have a lower mass planet to form a pressure maximum at 20 au. This situation relaxes the constraint of having the planet sitting at the minimum of the gap at 8 au. The planet-to-star mass ratio in this case is ~ 5 × 10⁻³, which corresponds to ~0.5 M_Jup around a 0.1 M_⊙ stellar object.

Figure 9 shows the dust density distribution for the case of a 0.5 M_Jup at 12 au at 0.1, 0.5, and 1.0 Myr of evolution. The dust density distribution of the dust evolution simulations are then combined with radiative transfer calculations using the Monte Carlo code RADMC-3D¹ as explainedin previous works (e.g., Pohl et al. 2016). Figure 10 shows the synthetic images at each time of evolution and for both wavelengths (without convolution to clearly see the differences). The dust density distribution shows that in this case the gap formed by this planet allows dust particles from the outer disk to move throughout the gap and replenish the inner disk. This inner dust-rich disk survives even after 1 Myr. With time, the dust density increases inside the ring because particles that are growing in the outer parts of the disk move towards the pressure maximum, which leads to higher intensity differences between the ring and the inner disk with time.

The synthetic images shown in Fig. 10 are normalized to the peak of the emission and they are in logarithmic scale. As seen in these images and in Fig. 8, there is always an inner emission. The outer tail is only visible in the images at 0.1 Myr of evolution and at 0.5 Myr and 1.0 Myr only the inner emission and the outer ring peaking at 20 au remain. Figure 11 shows the comparison of the deprojected radial profile from Fig. 3 and the models assuming a 0.5 M_Jup planet at 12 au after deprojecting and convolving with a beam similar to the one from observations. This comparison of the radial profiles illustrates how the models at 0.1 Myr nicely reproduce the contrast between the inner disk and the peak of the ring, and the width and depth of the gap in Band 7. The minimum of the gap (~ 8 au) is also in good agreement, even though the location of the planet is set to 12 au. In Band 4, the depth of the gap is higher than observed, although our visibility analysis suggests that the gap is shallower, in better agreement with these models. The main discrepancy between the radial profile at 0.1 Myr and the observations is the bright and large tail in the outer disk that is not seen in the radialprofile of the observations. At later times of evolution, while the outer tail vanishes, the inner disk also becomes less bright with respect to the peak. The discrepancy in the inner disk brightness between the observations and models is more noticeable for Band 4. Therefore, these models cannot explain the very bright inner disk at 1 My as observed in the absence of a mechanism helping to slow down the radial drift of the particles in the inner disk. One alternative is that this planet gradually grows and reaches its final mass in the last 0.1 Myr, such that the inner disk may still be as bright as observed.

To summarize these results of dust evolution models and comparison with observations, a very massive planet (as massive or more massive than Jupiter) is excluded because the gap would filter grains of all sizes and the inner disk would vanish very quickly. A 0.5 M_Jup allows replenishment of dust from the outer disk, which leads to good agreement between the contrast of the inner disk and the ring at early times of evolution (0.1 Myr). At longer times, it is challenging to keep this bright inner disk.

Fig. 8 Particle density distribution after 0.1, 0.5, and 1 Myr of evolution as a function of the distance from the star and grain size when a ~2.4 MJup planet is located at 8 au from the star.

Fig. 9 Particle density distribution after 0.1, 0.5, and 1 Myr of evolution as a function of distance from the star and grain size. The white solid line represents a Stokes number of unity. This model assumes the influence of a 0.5MJup planet located at 12 au in the disk.

5 Discussion

The high accretion rate and the large cavity of the disk of CIDA 1 exclude models of photo-evaporation (Ercolano & Pascucci 2017) as a potential explanation for the origin of the cavity. Moreover, the new detection of the inner disk with ALMAprovides further support to exclusion of photoevaporation as the mechanism to open this cavity. A possible alternative is the dust trapping that is expected at the outer edge of a dead zone (Flock et al. 2015). However, as discussed in Pinilla et al. (2018a) most of the nonideal MHD models have focused on T-Tauri or Herbig starsand little is known about how large a dead zone can be around a VLMS such as CIDA 1. In the remainder of the discussion, we focus on the possibility that the cavity is opened by a giant planet as in the previous section.

5.1 Inner disk

One of the main results from our observations is the detection of an unresolved inner disk in CIDA 1. Recent observations with ALMA revealed unresolved inner disk emission in several transition disks (e.g., T Cha, SR 24S, HD 100546, GG Tau, HD 100453, Hendler et al. 2018; Pinilla et al. 2019; Pérez et al. 2020; Phuong et al. 2020; Rosotti et al. 2020, respectively), see also Francis & van der Marel (2020).

As discussed in the previous section, the existence of these inner disks can provide upper limits for the potential embedded planets in cavities, because if the planet-to-star mass ratio is very high (≳ 5 × 10⁻³ − 1 × 10⁻²), the inner disk vanishes before 1 Myr. In the case of CIDA 1, the inner disk is very bright. Our visibility fits suggest that the ring peak is only ~30% of the peak of the inner disk emission (in Band 7) and in addition our spectral index analysis suggests that the inner disk and outer ring share similar values, which implies that both regions are rich in millimeter-sized particles assuming optically thin emission. This assumption of optically thin emission is probably inaccurate in the very hot inner parts of the disk (we obtained 160 K at 1 au from our radiative transfer calculations), in which case the spectral index does not provide any information about the grain size.

In our dust evolution models we have difficulty simultaneously explaining the inner disk (in particular its brightness) and the formation of the ring by one single massive planet. A less massive planet can help to have less dust filtration at the outer edge of the gap and hence a brighter inner disk with a higher replenishment of dust from the outer disk; but in such a case the pressure bump and the location of the ring would not be as far from the star as what is observed (at 20 au).

Our results also suggest that there may be a physical mechanism that helps to keep the inner disk brighter on million-year timescales. In our models, the grains that pass the gap and reach the inner disk are small dust particles (micron-sized particles), which then grow again and quickly drift towards the star. Hence, it is necessary to stop or slow down the radial drift of these dust particles in the inner disk to keep it bright. To slow down the radial drift of particles in the inner part of disks, Pinilla et al. (2016) considered the effect of the snow line in their dust evolution models. At the location of the snow line, small particles can be recreated due to the lower fragmentation velocities that are expected from dry grain particles compared to icy dust particles (see however recent results from laboratory experiments of dust collisions; Gundlach et al. 2018; Musiolik & Wurm 2019; Steinpilz et al. 2019). The inclusion of variation of the fragmentation velocity at the snow line (see also Gárate et al. 2020) can explain the NIR excess from the SEDs seen in a subset of transition disks (pre-transition disks). However, in the simulations from Pinilla et al. (2016), this snow-line effect does not produce a significant inner disk at submillimeter emission (see their Fig. A.1).

Another possibility is to trap the particles at the inner edge of a dead zone (e.g., Dzyurkevich et al. 2010; Hu et al. 2016; Ueda et al. 2019). Very close to the star where the temperatures are similar to the dust sublimation temperatures, the magnetorotational instability (MRI) is expected to operate due to thermal ionization. The inner edge of a dead zone has been proposed to be an excellent location for trapping of pebbles and for the formation of close-in planets. For the luminosity of CIDA 1, it is expected that the dead-zone inner edge is at about 0.01 au, similar to the case of Trappist 1 and the location of the closest-in planet in that system (see Fig. 8 from Flock et al. 2019). However, if the grains are too small near the inner edge of the dead zone, trapping is not expected (Jankovic et al. 2019), which can be the case inside the snow line.

Observations of exoplanets suggest that cold Jupiters may be more common when they are accompanied by a closer-in super-Earth (e.g., Zhu & Wu 2018; Bryan et al. 2019). If the large cavities observed in transition disks are the product of cold Jupiters embedded in the disk, it is possible that the inner disks remain due to the presence of an inner pressure bump created by a super-Earth, but detailed modeling is required to test this idea.

Alternatively, in the context of planets creating the cavity, it is possible that there are multiple planets whose gaps overlap and create a wider but shallower gap (Duffell & Dong 2015). In this case, it may be possible to still have a bright inner disk over longer timescales of evolution, an idea that needs confirmation with hydrodynamical simulations of multiple planets, dust evolution models, and radiative transfer calculations.

Fig. 10 Synthetic observations at 0.9 mm (top) and 2.1 mm (bottom) from dust evolution models and radiative transfer calculations after 0.1 Myr, 0.5 Myr, and 1.0 Myr of evolution for the case of one 0.5 MJup planet at 12 au creating a pressure bump at ~20 au.

Fig. 11 Comparison of the deprojected radial profile from Fig. 3 and the models from Fig. 9 (0.5 MJup planet at 12 au) after deprojecting and convolving with a beam similar to the one from observations.

5.2 Multiple rings or tail?

Our results from the visibility fits indicate that there is outer emission outside the main ring of the CIDA 1 disk. Our galario fits are in favor of an extra Gaussian on top of the main ring to reproduce that outer emission, which is only required in Band 7 and is negligible in the best model of the fit of the Band 4 observations. These results suggest that there is an outer shoulder in Band 7 or that the ring is actually a composition of two (or more) unresolved rings.

As discussed in the previous section, we expect a tail of emission in the outer disk from dust evolution models when a planet is present, which is significant compared to the ring at early times of evolution. The tail in the dust density distribution does not disappear at longer times, but because the ring becomes brighter due to the drift of millimeter-sized particles in the outer disk, the emission of the tail becomes negligible with time. This emission in the models also has a different nature from what our analysis of the visibilities suggests. In the dust evolution models, it looks more like extended emission towards the outer disk (see Figs. 11 and 10), while the best-fit models from visibility fitting using Band 7 data suggest that it looks more like an unresolved ring.

Recent observations with ALMA of several transition disks showed that what was known to be a wide ring around a cavity is actually a composition of several narrow rings (e.g., Pérez et al. 2020; Facchini et al. 2020). A couple of ideas have been proposed to explain these narrow rings. Pérez et al. (2020) suggested that a single migrating low-mass planet (10 M_⊕) in between these rings can explain the multiple rings and their separation, but this planet cannot explain the formation of the cavity itself.Facchini et al. (2020) demonstrated that the appearance of one ring or multiple rings depends on the assumptionsof the thermodynamics in hydrodynamical simulations where a planet is assumed to be embedded inside the cavity. Alternatively, low values of the disk viscosity, both in hydrodynamical simulations and dust evolution models, can change the appearance of the disk that is hosting a giant planet, spanning appearances that suggest a compact disk all the way to those suggesting multiple rings and gaps (e.g., de Juan Ovelar et al. 2016; Bae et al. 2018).

In the case of CIDA 1, it remains unclear as to whether or not this outer ring is a composition of multiple structures, and deeper observations at higher angular resolution are required to solve this problem and to explore potential scenarios.

5.3 Spectral index

Our observations suggest low average values of the spectral index of ~ 2. According to our calculations, the emission along the ring is mostly optically thin, which suggests that this low spectral index is due to a large number of millimeter-sized particles in the ring (the optical depth at the ring location from our radiative transfer models is much lower than unity). However, from the theoretical point of view of dust evolution models, it is challenging to reproduce these low values of the spectral index in disks around low-mass stars. This results from a combination of different factors. For example, radial drift velocities are higher around VLMSs than around T-tauri stars. This influences the spectral index in two different ways: (i) the millimeter-sized particles are lost towards the star on shorter timescales and (ii) drift mainly sets the maximum grain size of the particles (Pinilla et al. 2013). This last element together with these disks being lower in mass leads to a maximum grain size that is at most 1 cm in the inner disk (and usually decreases outwards), even when pressure bumps are present (see Fig. 9). This problem of having a low growth barrier, in particular in disks around VLMSs, is worse when taking into account the recent results from laboratory experiments that suggest fragmentation velocities of icy particles of 1 m s⁻¹ or even less (e.g., Musiolik & Wurm 2019).

Dust evolution simulations of disks around VLMSs, with one planet creating a bump and trapping the grains in a single ring, fail in reproducing very low spectral index values as observed in CIDA 1. Figure 12 shows the spectral index obtained from the intensity profiles from models that assume one single planet of 0.5 M_Jup at 12 au (Sect. 4). It is important to note that these models already include the full treatment of the dust scattering, which is known to decrease the spectral index values (Liu 2019).

In this figure, it is possible to see how the spectral index is expected to decrease inside the ring to values as low as 2.5 and to values of ~ 2.6 when integrating in the entire disk. At 1 Myr of evolution, 20% of the dust has been lost due to drift in the dust evolution model that considers a 0.5 M_Jup at 12 au, but a very small fraction of the remaining mass is in millimeter-sized particles. The dust mass in particles larger than 0.1 mm is only 0.3 M_⊕, while the dust mass in the smaller grains is 13.7M_⊕, which means that only 2% of the mass is in large grains (≥0.1 mm). From simple radiative transfer calculations, a way to mitigate this problem is to change the fraction of small versus large grains, which already provides values of the spectral index as observed (i.e., ~ 2). To have this switch in our dust evolution models, we would need multiple and strong pressure bumps in the entire disk that allow to accumulate a lot of millimeter-sized particles in the whole disk (as in Pinilla et al. 2013). However, our high-angular-resolution observationsonly suggest the presence of one single particle trap and therefore it remains unclear as to how we can explain the low spectral index of CIDA 1. This will require further investigation with dust evolution models to search for the set of parameters that allows us to have the large majority of the dust in large particles, including fragmentation velocities and gas and dust diffusion parameters.

Fig. 12 Spectral index profile obtained from the intensity profiles from models that assume one single planet of 0.5 MJup at 12 au.

6 Conclusions

In this paper we present new ALMA and X-shooter observations of the disk around CIDA 1 (0.1–0.2 M_⊙), which is one of the very few known disks that hosts a large cavity (20 au in size) around a very low mass star. Our ALMA observations include Band 7 (0.9 mm) and Band 4 (2.1 mm) both with a resolution of ~ 0.05″ × 0.034″. Our main findings can be summarized as follows.

• The X-shooter observations confirm the accretion rate of CIDA 1 of Ṁ_acc = 1.4 × 10⁻⁸ M_⊙ yr⁻¹ for its stellar mass. This high value of Ṁ_acc and the large cavity of 20 au exclude photo-evaporation as a potential origin of the observed structures.

• The new ALMA observations reveal an unresolved inner disk in the continuum emission at the two wavelengths. This inner disk is very bright when compared with the emission of the main ring, which peaks at 20 au. The discovery of this inner disk helps to constrain models for the formation of the observed cavity.

• Our experiments with dust evolution models that include a massive planet in the disk reveal an inner disk that is fainter than observations suggest for this system, demonstrating that a massive planet (≳0.5 M_Jup) is unlikely to exist in this disk. Our models of dust evolution with a 0.5 M_Jup planet at 12 au can explain the contrast of the inner disk and ring only at 0.1 Myr of evolution. After 1 Myr the inner disk becomes fainter because a nonphysical mechanism stops the radial drift of the particles in the inner disk. A potential alternative is that the planet around CIDA 1 has gradually grown and reached its final mass in the last 0.1 Myr.

• There are additional possibilities to explain a long-lived inner disk as observed with ALMA for CIDA 1. These include a traffic jam expected at the snow line due to variation of fragmentation velocities of the particles, the inner edge of dead zones, and/or multiple planets inside the cavity; e.g., a super-Earth in the very inner part of the disk. The effects of the inclusion of these scenarios in planet–disk interaction and dust evolution models remain to be explored.

• Besides the inner disk and the ring-like structures, there is no further substructure found in CIDA 1 at the current resolution and sensitivity. There is a hint of a shoulder beyond the ring-like structure from our fit of the visibilities in Band 7. It is possible that the ring is composed of two or more narrow rings as observed in other transition disks, and higher angular resolution observations are needed to verify this possibility.

• There is weak evidence of a decrease in the spectral index inside the ring, both from observations and from our fits of the visibilities. The spatially integrated spectral index is very low (~2), which isvery challenging to explain in dust evolution models. Better agreement between models and observations is obtained if most of the dust inside the ring of the models remains as millimeter-sized particles, which calls into questionour understanding of the different physical processes, such as growth, fragmentation, and diffusion inside pressure bumps.

Acknowledgements

P.P. and N.T.K. acknowledge support provided by the Alexander von Humboldt Foundation in the framework of the Sofja Kovalevskaja Award endowed by the Federal Ministry of Education and Research. This work was partly supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - Ref no. FOR 2634/1 TE 1024/1-1. M.T. has been supported by the UK Science and Technology research Council (STFC) via the consolidated grant ST/S000623/1, and by the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 823823 (RISE DUSTBUSTERS project). S.M.S. acknowledges funding from the European Research Council (ERC) under the European Unions Horizon 2020 research and innovation programme under grant agreement No. 714769, the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy - EXC-2094 - 390783311, and the support from the DFG Research Unit “Transition Disks” (FOR 2634/1). This paper makes use of the following ALMA data: ADS/JAO.ALMA#2018.1.00536.S, ADS/JAO.ALMA#2015.1.00934.S, and ADS/JAO.ALMA#2016.1.01511.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ.

References

All Tables

All Figures

	Fig. 1 ALMA Band 7 (0.9 mm) and Band 4 (2.1 mm) continuum images of the CIDA 1 disk. The synthesized beam is shown in the lower left corner of each image. The details of the images are summarized in Table 2.
In the text

	Fig. 2 Moment 0 and moment 1 of the 12CO and 13CO. The upper right bar corresponds to 10 au size. The dashed circles (projected with the same inclination and PA as CIDA 1, Table 3) are 1.0″ and 0.6″ in size from the center for the 12CO and 13CO, respectively. The details of the images are summarized in Table 2.
In the text

Fig. 3 Azimuthally averaged radial intensity profiles of the deprojected continuum images (lineal scale on the left and logarithmic scale in the center) at each wavelength. Each profile is normalized to the peak. The profiles are obtained from images that are restored with the same beam (50 mas × 34 mas). The error bars include the standard deviation of each elliptical bin divided by the square root of the number of beams spanning the full azimuthal angle at each radial bin. The right panel shows the radial profile of the spectral index taking into account the errors on the radial intensity profiles in addition to 10% of the error due to flux calibration.

In the text

Fig. 4 From left to right: real part of the visibilities after centering and deprojecting the data in Band 7 and in Band 4 vs. the best-fit model from galario of the continuum data. Intensity profile from the best fit in Band 7 and Band 4 (orange and blue line, respectively), with 5000 randomly selected chains after convergence overlaid. Right panel: spectral index radial profile obtained from the intensity profiles.

In the text

	Fig. 5 Observations vs. model images from our best fit with galario (before and after convolution) and the residuals. Top panels correspond to Band 7 and the bottom panels correspond to Band 4.
In the text

	Fig. 6 Accretion luminosity derived from the luminosity of various emission lines, as labeled, using the relations by Alcalá et al. (2017), and assuming AV = 3.5 mag. The blue line is the best-fit accretion luminosity obtained fitting the whole X-shooter spectrum of CIDA 1. The red line is the mean value of Lacc obtained from the luminosity of the emission lines.
In the text

	Fig. 7 Fit of the spectrum of CIDA 1 obtained here. The dereddened spectrum of the target (red) is reproduced by the sum of a photospheric template (green) and a slab model to reproduce the accretion shock emission (cyan).
In the text

	Fig. 8 Particle density distribution after 0.1, 0.5, and 1 Myr of evolution as a function of the distance from the star and grain size when a ~2.4 MJup planet is located at 8 au from the star.
In the text

	Fig. 9 Particle density distribution after 0.1, 0.5, and 1 Myr of evolution as a function of distance from the star and grain size. The white solid line represents a Stokes number of unity. This model assumes the influence of a 0.5MJup planet located at 12 au in the disk.
In the text

	Fig. 10 Synthetic observations at 0.9 mm (top) and 2.1 mm (bottom) from dust evolution models and radiative transfer calculations after 0.1 Myr, 0.5 Myr, and 1.0 Myr of evolution for the case of one 0.5 MJup planet at 12 au creating a pressure bump at ~20 au.
In the text

	Fig. 11 Comparison of the deprojected radial profile from Fig. 3 and the models from Fig. 9 (0.5 MJup planet at 12 au) after deprojecting and convolving with a beam similar to the one from observations.
In the text

	Fig. 12 Spectral index profile obtained from the intensity profiles from models that assume one single planet of 0.5 MJup at 12 au.
In the text