# ALMA Observations of Elias 2-24: A Protoplanetary Disk with Multiple Gaps in the Ophiuchus Molecular Cloud

Lucas A. Cieza, Simon Casassus, Sebastian Pérez, Antonio Hales, Miguel Cárcamo, Megan Ansdell, Henning Avenhaus, Amelia Bayo, Gesa H. -M. Bertrang, Hector Cánovas, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Lucas A. Cieza, Simon Casassus, Sebastian Pérez, Antonio Hales, Miguel Cárcamo, et al. ALMA Observations of Elias 2-24: A Protoplanetary Disk with Multiple Gaps in the Ophiuchus Molecular Cloud. The Astrophysical Journal Letters, 2017, 851, 10.3847/2041-8213/aa9b7b. insu-03692462

HAL Id: insu-03692462

https://insu.hal.science/insu-03692462v1

Submitted on 28 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# ALMA Observations of Elias 2–24: A Protoplanetary Disk with Multiple Gaps in the Ophiuchus Molecular Cloud

Lucas A. Cieza<sup>1,2 Ⓓ</sup>, Simon Casassus<sup>2,3 Ⓓ</sup>, Sebastian Pérez<sup>2,3 Ⓓ</sup>, Antonio Hales<sup>4,5 Ⓓ</sup>, Miguel Cárcamo<sup>6</sup>, Megan Ansdell<sup>7 Ⓓ</sup>, Henning Avenhaus<sup>8</sup>, Amelia Bayo<sup>9,10 Ⓓ</sup>, Gesa H.-M. Bertrang<sup>2,3</sup>, Hector Cánovas<sup>11 Ⓓ</sup>, Valentin Christiaens<sup>2,3</sup>, William Dent<sup>4</sup>, Gabriel Ferrero<sup>12</sup>, Roberto Gamen<sup>12</sup>, Johan Olofsson<sup>9,10 Ⓓ</sup>, Santiago Orcajo<sup>12</sup>, Axel Osses<sup>13</sup>, Karla Peña-Ramirez<sup>14</sup>, David Principe<sup>15</sup>, Dary Ruíz-Rodríguez<sup>16 Ⓓ</sup>, Matthias R. Schreiber<sup>9,10</sup>, Gerrit van der Plas<sup>17</sup>, Jonathan P. Williams<sup>18 Ⓓ</sup>, and Alice Zurlo<sup>1,2 Ⓓ</sup>

<sup>1</sup> Facultad de Ingeniería y Ciencias, Núcleo de Astronomía, Universidad Diego Portales, Av. Ejercito 441, Santiago, Chile

<sup>2</sup> Millennium Nucleus “Protoplanetary Disks in ALMA Early Science,” Av. Ejercito 441, Santiago, Chile

<sup>3</sup> Departamento de Astronomía, Universidad de Chile, Casilla 36-D Santiago, Chile

<sup>4</sup> Joint ALMA Observatory, Alonso de Cordova 3107, Vitacura 763-0355, Santiago, Chile

<sup>5</sup> National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903-2475, USA

<sup>6</sup> Departamento de Ingeniería Informática, Universidad de Santiago de Chile, Av. Ecuador 3659, Santiago, Chile

<sup>7</sup> Department of Astronomy, University of California at Berkeley, Berkeley, CA 94720-3411, USA

<sup>8</sup> ETH Zurich, Institute for Particle Physics and Astrophysics, Wolfgang-Pauli-Strasse 27, CH-8093 Zurich, Switzerland

<sup>9</sup> Facultad de Ciencias, Instituto de Física y Astronomía, Universidad de Valparaíso, Av. Gran Bretaña 1111, 5030 Casilla, Valparaíso, Chile

<sup>10</sup> Millennium Nucleus for Planet Formation, Universidad de Valparaíso, Av. Gran Bretaña 1111, Valparaíso, Chile

<sup>11</sup> European Space Astronomy Centre (ESA), Camino Bajo del Castillo s/n, E-28692, Villanueva de la Cañada, Madrid, Spain

<sup>12</sup> Instituto de Astrofísica de La Plata y Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, Paseo del Bosque s/n, La Plata, Argentina

<sup>13</sup> El Departamento de Ingeniería Matemática, Facultad de Ciencias Físicas y Matemáticas, Universidad de Chile, Beauchef 851, Santiago, Chile

<sup>14</sup> Unidad de Astronomía de la Universidad de Antofagasta, Av. U. de Antofagasta, 02800 Antofagasta, Chile

<sup>15</sup> Massachusetts Institute of Technology, Kavli Institute for Astrophysics, Cambridge, MA, USA

<sup>16</sup> Chester F. Carlson Center for Imaging Science, Rochester Institute of Technology, Rochester, NY 14623-5603, USA

<sup>17</sup> Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France

<sup>18</sup> Institute for Astronomy, University of Hawaii at Manoa, Honolulu, HI 96822, USA

Received 2017 October 28; revised 2017 November 13; accepted 2017 November 16; published 2017 December 11

## Abstract

We present ALMA 1.3 mm continuum observations at $0.′′2$ (25 au) resolution of Elias 2–24, one of the largest and brightest protoplanetary disks in the Ophiuchus Molecular Cloud, and we report the presence of three partially resolved concentric gaps located at ~20, 52, and 87 au from the star. We perform radiative transfer modeling of the disk to constrain its surface density and temperature radial profile and place the disk structure in the context of mechanisms capable of forming narrow gaps such as condensation fronts and dynamical clearing by actively forming planets. In particular, we estimate the disk temperature at the locations of the gaps to be 23, 15, and 12 K (at 20, 52, and 87 au, respectively), very close to the expected snowlines of CO (23–28 K) and $N_2$ (12–15 K). Similarly, by assuming that the widths of the gaps correspond to $4–8×$ the Hill radii of forming planets (as suggested by numerical simulations), we estimate planet masses in the range of $0.2–1.5 M_{Jup}$, $1.0–8.0 M_{Jup}$, and $0.02–0.15 M_{Jup}$ for the inner, middle, and outer gap, respectively. Given the surface density profile of the disk, the amount of “missing mass” at the location of each one of these gaps (between 4 and 20 $M_{Jup}$) is more than sufficient to account for the formation of such planets.

**Key words:** circumstellar matter – planetary systems – protoplanetary disks – stars: individual (Elias 2–24) – techniques: interferometric

## 1. Introduction

Gas-rich circumstellar disks are the sites of planet formation. Since most of them are found still embedded in molecular clouds, only a handful of protoplanetary disks are located within 100 pc of Earth. Given their characteristic temperature and size ($T=20$ K and $r \lesssim 100$ au; Williams & Cieza 2011), resolving their thermal emission is best achieved by observations at (sub)millimeter wavelengths with sub-arcsecond angular resolution. Imaging protoplanetary disks in great detail has been one of the main scientific drivers for building the Atacama Large Millimeter/submillimeter Array (ALMA), which in recent years has revolutionized our view of protoplanetary disks thanks to its unprecedented resolution. Before ALMA, few structures were seen within these disks, mostly in the form of large central cavities, tens of au in radius (e.g., Brown et al. 2009; Andrews et al. 2010; Cieza et al. 2012). When observed at higher resolution and sensitivity,

protoplanetary disks show a variety of substructures such as narrow gaps (ALMA Partnership et al. 2015; Andrews et al. 2016; Isella et al. 2016), bright rings (Canovas et al. 2016; van der Plas et al. 2017), dust traps (Casassus et al. 2013; van der Marel et al. 2013; Kraus et al. 2017), spiral arms (Pérez et al. 2016), and sharp intensity breaks (Cieza et al. 2016). The origin of these structures and their role in disk evolution and planet formation is currently one of the major questions in the field.

Here, we present new Band 6 (1.3 mm) observations of Elias 2–24, a well-studied K5 star (Prato et al. 2003) hosting one of the largest and brightest protoplanetary disks in the nearby Ophiuchus Molecular Cloud (distance ~125 pc; Loinard et al. 2008). Andrews et al. (2010) observed this object with a $0.′′65 × 0.′′51$ beam at 880 µm and detected disk emission up to ~1" from the star with a total flux of 890 mJy. The disk is close to face-on ($i \sim 25$ deg) and has an estimated mass of
~0.12 M☉, while the central object is a heavily accreting T Tauri star (M˙* ~ 2 × 10⁻⁷ M☉ yr⁻¹; Natta et al. 2006) with a mass of ~1.0 M☉ and an estimated age of just 0.4 Myr (Siess et al. 2000; Andrews et al. 2010). We observed Elias 2–24 as part of the Ophiuchus Disk Survey Employing ALMA (ODISEA; L. Cieza et al. 2017, in preparation), a program studying 147 Ophiuchus objects at 0′′2 (25 au) resolution. The size, brightness, and orientation of its disk allow us to search for substructure in Elias 2–24 using what now can be considered modest resolution.

## 2. ALMA Observations and Data Analysis

### 2.1. 1.3 mm Observations

Elias 2–24 was observed by ALMA in Band 6 (230 GHz/ 1.3 mm) as part of the Cycle 4 program 2016.1.00545.S on 2017 July 13 and 14, in three different execution blocks. The precipitable water vapor (PWV) ranged from 1.1 to 1.9 mm in the three different observing sessions. During the observations, 42–45 of the 12 m ALMA antennas were used with baselines ranging from 16.7 to 2647.3 m. The ALMA correlator was configured with two spectral windows with 1.875 GHz bandwidths for continuum observations centered at 232.6 and 218.0 GHz. Also, three spectral windows were placed to cover the ¹²CO (2–1), ¹³CO (2–1), and C¹⁸O (2–1) transitions of carbon monoxide at 230.5380, 220.3987, and 219.5603 GHz, respectively. The first spectral window has a 0.04 km s⁻¹ spectral resolution, while the other two have a 0.08 km s⁻¹ resolution. J1517–2422 and J1733–1304 served as flux calibrators, while the quasars J1517–2422 and J1625–2527 were observed for bandpass and phase calibration, respectively. The total integration time on Elias 2–24 was 45 s.

### 2.2. Data Analysis

All the ODISEA data were calibrated using the Common Astronomy Software Applications package (CASA v4.2.1; McMullin et al. 2007) by the ALMA observatory, including the standard bandpass, phase, and amplitude calibrations; the offline Water Vapor Radiometer (WVR) calibration; and system temperature corrections. Online and nominal flagging, such as shadowed antennas and band edges, were applied for calibration. The observations from all 3 execution blocks were concatenated and processed together to increase the signal to noise and uv-coverage. We also used CASA to image all ODISEA data using the standard clean algorithm, with uniform weightings. The cleaning resulted in a synthesized beam of 0′′20 × 0′′25 and a continuum rms of 0.26 mJy beam⁻¹. After visual inspection, Elias 2–24 immediately stood out as one of the brightest and largest disks of the 147 objects in the survey and showed a clear ring pattern, with gaps at ~0′′4 and ~0′′7 (see Figure 1(a)). Applying one iteration of phase-only self-calibration reduced the noise level to 0.19 mJy beam⁻¹ for this object, consistent with the expected thermal noise. The CO line was detected with a peak emission of ~60 mJy beam⁻¹ in 0.25 km s⁻¹ channels and showed rotation broadly consistent with a Keplerian disk, but is not further discussed in this Letter. The ¹³CO and C¹⁸O lines were not detected.

Provided with a high dynamic range (>200) in the continuum, we investigated the super-resolution of the visibilities with non-parametric image synthesis. We used the uvmem package

(Casassus et al. 2006, 2015) to fit a model image to the data in a least χ² sense. In Figure 1(b), we show a deconvolved model image {$I_i$} using a measure of entropy regularization and the following objective function: $L = χ^2 - λS$, with $χ^2 = \sum_k ω_k ||V_k^o - V_k^m||^2$ and $S = \sum_i I_i \log(\frac{I_i + G^*(η+1)}{G})$. Here, $V^o$ and $V^m$ are the observed and model visibility data, $λ = 0.1$, $G = 30$, and $η = 1.0$. We also produce a pure χ² model (Figure 1(d)) with positivity regularization ($I_i > 0$), which allows for finer angular resolution at the expense of lowering the signal-to-noise ratio. The uncertainties in the model are calculated with Monte Carlo simulations, specifically 100 different injections of Gaussian random noise in the visibilities (see Cárcamo et al. 2017 for technical details regarding uvmen). These uvmem images have a pixel size of 0′′02 and a variable resolution (in the 0.1–0′′2 range, depending on the local signal to noise) that is a factor of ~2–3 higher than in the clean image with uniform weights. The clean and uvmem images are qualitatively very similar. The main difference is the depths of the two gaps, which are deeper at higher resolution, suggesting that these are partially resolved features. The integrated continuum flux is 345 mJy ± 35 mJy, where the error reflects the absolute flux calibration. An elliptical Gaussian fit to the model image indicates a position angle of 43.4° (east of north) and an inclination of 23.6°. Using this inclination and position angle, we create a higher signal-to-noise image averaging all pixels along concentric ellipses (Figure 1(f)).

Figure 2 shows a deprojected radial profile of this average image and a cut along the semimajor axis of the disk in the pure χ² uvmem image. This radial cut maximizes the resolution of the data and allows us to identify a possible third gap at 0′′18 thanks to the increased resolution of uvmen toward the center of the image. This third gap survives the different realization of noise injected to the visibilities and it is seen as a change in slope in the deprojected radial average. We thus consider it likely to be real. There is a small hint of a fourth gap at 0′′9, but it has a low significance and it is not further discussed in this Letter.

## 3. Radiative Transfer Modeling

In order to derive the mass, temperature, and surface density profiles of the Elias 2–24 disk, we perform radiative transfer modeling using the code RADMC-3D (Dullemond et al. 2012). While the disk has barely resolved substructures in the form of gaps, the large-scale structure can be approximated as a continuous disk. We adopt the same parameterization, stellar parameters, dust properties, and gas to dust mass ratio (i.e., 100) used by Andrews et al. (2010). In particular, the surface density profile of the disk, $Σ(r)$ is given by

$$ \Sigma(r) = \Sigma_c \left(\frac{r}{R_c}\right)^{-\gamma} \exp\left[-\left(\frac{r}{R_c}\right)^{2-\gamma}\right], \quad (1) $$

where $R_c$ is the disk’s characteristic radius. Similarly, the vertical scale height of the disk as a function of radius $h(r)$ is described as

$$ h(r) = h_{100} \left( \frac{r}{R_{100}} \right)^{\psi}, \quad (2) $$

where $h_{100}$ is the scale height at 100 au and $ψ$ defines how this scale height increases with radius, as expected for a flared disk. In this context, the structure of the disk can be fully described
<figure><img src="image_4.png" /><figcaption>**Figure 1.** The 1.3 mm *clean* image of Elias 2-24 with uniform weights (a) and the deconvolved model image using the *uvmem* package with the parameters described in the text (b). The best-fitting model (c) and the pure $\chi^2$ *uvmem* image (d). The radiative transfer residual (e), where the outer gaps are seen as negative features as they are not included in the modeling. The asymmetric structure close to the origin is likely to be due to the innermost gap, which is only marginally resolved. The elliptical average of the pure $\chi^2$ *uvmem* image (f).</figcaption></figure>

by five free parameters: $R_c$, $\Sigma_c$, $\gamma$, $h_{100}$, and $\psi$. To facilitate comparisons to other objects, $\Sigma_c$ can be replaced by $M_{\text{disk}}$ (gas + dust disk mass) by integrating Equation (1).

The parameter space {$M_{\text{disk}}, \gamma, R_c, h_{100}, \psi$} was explored using a Bayesian approach, which is described in more detail in S. Pérez et al. (2017, in preparation) and is based on the
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Radial brightness profile of Elias 2–24 at 1.3 mm. The gray line and red shading correspond to the *deprojected* radial average of the image (Figure 1(f)) and the errors in this average. The blue line corresponds to a cut along the semimajor axis of the pure $\chi^2$ uvmem image (Figure 1(d)) and maximizes the resolution of the data. In this case, the error in the profile (blue shading) is derived from 100 Monte Carlo realizations of noise injected to the visibilities. The vertical lines mark the location of the gaps, which correspond to distances of ~20, 52, and 87 au. The dotted line shows the midplane temperature profile of our continuous disk model with the scale on the right, while the horizontal lines indicate the approximate temperatures at the location of the gaps.</figcaption></figure>

Goodman & Weare’s Affine Invariant Markov Chain Monte Carlo (MCMC) ensemble sampler (Foreman-Mackey et al. 2013) and the publicly available PYTHON module *emcee*. The chains were initialized in a uniform distribution around the parameters reported in Andrews et al. (2010), and we used 100 walkers for 500 steps, which allowed the likelihood function to reach steady maximum values. The posterior distributions and best-fit parameters obtained from the MCMC exploration are shown in Figure 3. Our model parameters agree well (within $1\sigma-2\sigma$) with those derived by Andrews et al. (2010) using lower-resolution data. The only exception is the $\psi$ parameter, which is best constrained by the spectral energy distribution (not included in our modeling).

The comparison between the model and the observations is shown in Figures 1(a)–(e), where the deep gap at ~0.4 is clearly visible in the residual as well as hints of the gaps at ~0.18 and ~0.7. From our radiative transfer model, we can calculate the temperature of the disk midplane at the location of the gaps: ~23 K at 20 au, ~15 K at 52 au, and ~12 K at 87 au (see Figure 2). We can also estimate the mass of the “missing material” at each gap by integrating the surface density profile of the best-fit model over the width of the gap. Since the gaps are only partially resolved, there is a degeneracy between the width and the depth of the gaps. We estimate the widths of the gaps by assuming that they are fully evacuated and measuring their equivalent widths in the brightness profile along the major-axis (blue line in Figure 2). Gaussian fits to the gaps using the `splot` routine within the IRAF package `noao`. `onedspec` indicate that the gaps are located at 20 ± 3 au, 52 ± 2 au, and 87 ± 3 au and have widths of 6 ± 2 au, 28 ± 3 au, and 11 ± 4 au. Integrating the area under the curve between the edges of each gap, and/or using the deprojected radial average (gray line in Figure 2, in the case of the two outer gaps), produce consistent equivalent widths. From these gap widths, we derive “missing masses” of ~4 $M_{\text{Jup}}$, 20 $M_{\text{Jup}}$, and 10 $M_{\text{Jup}}$ for the gaps at 20, 52, and 87 au, respectively.

## 4. Discussion

ALMA observations at high resolution and sensitivity have revealed concentric gaps in several sources, including HL Tau (ALMA Partnership et al. 2015), TW Hydra (Andrews et al. 2016), and HD 163296 (Isella et al. 2016). However, the origin of these gaps and rings still remains to be established. Several potential explanations have been offered so far, including dynamical clearing by forming planets (Yen et al. 2016; Dong et al. 2017), snowlines (Zhang et al. 2015; Okuzumi et al. 2016), magnetohydrodynamic effects (Ruge et al. 2016; Flock et al. 2017), and viscous ring-instability (Dullemond & Penzlin 2017).

Distinguishing between the different potential explanations for concentric gaps in individual objects seems difficult and might require the measurements of magnetic fields, dust kinematics, more detailed modeling, and/or the direct detection of planets within the gaps. Deep high-contrast searches for planets in the TW Hydra system yield no detections (Ruane et al. 2017), but place limits of 1.2–2.5 $M_{\text{Jup}}$ for the mass of the putative planets at the location of the main gaps. However, these limits are not enough to rule out the planetary origin of the gaps as less massive planets might be responsible for them (e.g., Dong et al. 2017). The demographics of the gaps might also shed some light on their origin. For instance, if gaps are produced by snowlines, they must be ubiquitous and their location should correlate with the temperature and luminosity of the central source. On the other hand, if planets are responsible for these gaps, their location should not depend on disk temperature, but their widths and depths should correlate with the mass of the disk. In the case of our Ophiuchus survey (L. Cieza et al. 2017, in preparation), among 147 sources, only the disk around Elias 2–24 has the right combination of size, brightness, and orientation that allows the identification of three gaps. Deeper, higher-resolution images are needed to investigate the overall incidence and properties of concentric gaps in the Ophiuchus Molecular Cloud and other star-forming regions.

Meanwhile, we can place the results presented in Section 3 in the context of opacity gaps created by condensation fronts and the dynamical clearing by forming planets. Zhang et al. (2015) suggest that the most prominent gaps seen in HL Tau (at 13, 32, and 63 au) are due to the snowlines of water (~144 K) pure ammonia or ammonia hydrates (~84 K) and clathrate hydrates (~57 K). In the case of Elias 2–24, all these temperatures correspond to stellocentric distances <20 au, which are not resolved by our observations. The disk temperatures we derive for the location of the gaps (23, 15, and 12 K) are instead very close to the condensation fronts of other species: CO and $N_2$, which are expected to occur at temperatures of 23–28 K and 12–15 K (Mumma & Charnley 2011; Martín-Doménech et al. 2014). Interestingly, Andrews et al. (2016) report that the gaps seen at 22 and 37 au in TW Hydra are also located close to the expected snowlines of CO and $N_2$ according to their thermal model. However, the gaps in TW Hydra are narrower than those of Elias 2–24, both in terms of absolute size (1–6 au versus 6–28 au) and as a fraction of disk radius (<2%–8% versus 5%–20%).

To date, the closest analog to Elias 2–24 in terms of structure seems to be HD 163296 (Isella et al. 2016). The 1.3 mm images of both objects are remarkably similar, but HD 163296 is actually a scaled-up version (by a factor of two) of Elias 2–24 with respect to the stellar mass, the disk size, the location
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Posterior distributions of each of the five parameters used to model Elias 2–24 as a continuous disk. The vertical dashed lines represent the 16th, 50th, and 84th percentiles.</figcaption></figure>

of the gaps and their widths. HD 163296 has a dust disk 250 au in radius with three gaps at 60, 100, and 160 au. These gaps have estimated widths of 25, 22, and 45 au, respectively. Isella et al. (2016) argue that the CO snowline in HD 163296 could be anywhere between 50 and 180 au and therefore do not associate frostlines with any particular gap. Instead, they estimate the masses of the possible planets that could explain the observed gaps. They note that numerical simulations show that planets can open gaps that are 4–8 times the Hill radius (Wolf et al. 2007; Rosotti et al. 2016), $r_{Hill}$, which is given by

$$r_{\text{Hill}} \sim a \times \left( \frac{m_p}{3M_*} \right)^{1/3}, \quad (3)$$

where a is the semimajor axis of the planet's orbit, $m_p$ is the mass of the planet, and $M_*$ is the mass of the central star, $2.3 M_\odot$ in the case of HD 163296. Using this argument, they estimate that planets with masses between 0.05 and 0.5 $M_{\text{Jup}}$ are needed to explain the gaps in HD 163296. Applying the same approach to Elias 2–24, we derive planet masses in the range of 0.2–1.5 $M_{\text{Jup}}$, 1.0–8.0 $M_{\text{Jup}}$, and 0.02–0.15 $M_{\text{Jup}}$ for the gaps at 20, 52, and 87 au, respectively. These values assume the gap widths adopted in Section 3 (6, 28, and 11 au, respectively) and a stellar mass of $1.0 M_\odot$ (Andrews et al. 2016). The planet masses so derived are very modest compared to “missing masses,” also reported in Section 3, based on our continuous disk model (4–20 $M_{\text{Jup}}$). We thus conclude that there was more than enough mass at the location of the gaps (as an initial condition) to account for the formation of the putative planets. In this scenario, most of the “missing mass” would likely be pushed away by the tidal forces of the planet toward the edges of the gap, while only a small fraction of this mass would be accreted by each planet (Szulágyi et al. 2014; Dipierro & Laibe 2017). We note that the two explanations discussed above for the structure of Elias 2–24 (condensation fronts and dynamical clearing by forming planets) might be closely related, as one can imagine combined scenarios in which the conditions at the condensation fronts are responsible for the first steps of planet formation (the growth of pebbles and planetesimals) that eventually lead to large planets capable of dynamically clearing the gaps. However, opacity gaps produced by snowlines would only require the formation of centimeter-sized pebbles, while

dynamical clearing would imply the formation of giant planets at ~20–85 au by the age of the system, <1.0 Myr in the case of Elias 2–24 (Andrews et al. 2010).

## 5. Summary and Conclusions

We have obtained ALMA 1.3 mm observations of Elias 2–24 as part of the ODISEA program, which is a survey of 147 objects in Ophiuchus at ~25 au resolution. The disk shows a structure of three concentric gaps with estimated widths ranging from ~6 to 28 au. Radiative transfer modeling of the source indicates that the disk temperature at the location of the gaps is close to the expected snowlines of CO and $N_2$, consistent with claims that frostlines can result in gaps of dust opacity. The surface density profile of the disk is also consistent with formation of planets with masses similar to those of the giant planets in the solar system (from Neptune-mass to a few times the mass of Jupiter). However, other potential explanations still exist, including magnetohydrodynamic effects and disk instabilities. Elias 2–24 is one of the brightest and largest disks in the Ophiuchus Molecular Cloud, and it remains to be established whether smaller versions of these gaps (e.g., with locations and widths scaled down by disk radii and/or stellar luminosity) are typical of protoplanetary disks. Previous observations of TW Hydra at much higher resolution (~1 au) suggest that this might be the case.

This Letter uses the following ALMA data: ADS/JAO. ALMA #2016.1.00545.S. ALMA is a partnership of ESO, NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The NRAO is a facility of the NSF operated under cooperative agreement by Associated Universities, Inc. L.A.C., S.C., G.H.M.B., A.O., and A.Z. were supported by CONICYT-FONDECYT grant numbers 1171246, 1171624, 3170204, 1151512, and 3170657. L.A.C., S.C., S.P., A.B., G.H.M.B., J.O., M.R.S. and A.Z. acknowledge support from the Millennium Science Initiative (Chile) through grant RC130007 and the Núcleo Milenio de Formación Planetaria, NPF. This Letter use the Brelka cluster, financed by Fondequip project EQM140101. G.P. is supported by ANR (France) under #ANR-16-CE31-0013 (Planet-Forming-Disks).
ORCID iDs

Lucas A. Cieza https://orcid.org/0000-0002-2828-1153

Simon Casassus https://orcid.org/0000-0002-0433-9840

Sebastian Pérez https://orcid.org/0000-0003-2953-755X

Antonio Hales https://orcid.org/0000-0001-5073-2849

Megan Ansdell https://orcid.org/0000-0003-4142-9842

Amelia Bayo https://orcid.org/0000-0001-7868-7031

Hector Cánovas https://orcid.org/0000-0001-7668-8022

Johan Olofsson https://orcid.org/0000-0003-4475-3605

Dary Ruíz-Rodríguez https://orcid.org/0000-0003-3573-8163

Jonathan P. Williams https://orcid.org/0000-0001-5058-695X

Alice Zurlo https://orcid.org/0000-0002-5903-8316

References

ALMA Partnership, Brogan, C. L., Pérez, L. M., et al. 2015, ApJL, 808, L3

Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C., & Dullemond, C. P. 2010, ApJ, 723, 1241

Andrews, S. M., Wilner, D. J., Zhu, Z., et al. 2016, ApJL, 820, L40

Brown, J. M., Blake, G. A., Qi, C., et al. 2009, ApJ, 704, 496

Canovas, H., Caceres, C., Schreiber, M. R., et al. 2016, MNRAS, 458, L29

Cárcamo, M., Román, P., Casassus, S., Moral, V., & Rannou, F. R. 2017, arXiv:1703.02920

Casassus, S., Cabrera, G. F., Förster, F., et al. 2006, ApJ, 639, 951

Casassus, S., Marino, S., Pérez, S., et al. 2015, ApJ, 811, 92

Casassus, S., van der Plas, G., Perez, M. S., et al. 2013, Natur, 493, 191

Cieza, L. A., Casassus, S., Tobin, J., et al. 2016, Natur, 535, 258

Cieza, L. A., Mathews, G. S., Williams, J. P., et al. 2012, ApJ, 752, 75

Dipierro, G., & Laibe, G. 2017, MNRAS, 469, 1932

Dong, R., Li, S., Chiang, E., & Li, H. 2017, ApJ, 843, 127

Dullemond, C. P., Juhasz, A., Pohl, A., et al. 2012, RADMC-3D: A multi-purpose radiative transfer tool, Astrophysics Source Code Library, ascl:1202.015

Dullemond, C. P., & Penzlin, A. B. T. 2017, A&A, in press (arXiv:1709.00246)

Flock, M., Nelson, R. P., Turner, N. J., et al. 2017, ApJ, in press (arXiv:1710.06007)

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306

Isella, A., Guidi, G., Testi, L., et al. 2016, PhRvL, 117, 251101

Kraus, S., Kreplin, A., Fukugawa, M., et al. 2017, ApJL, 848, L11

Loinard, L., Torres, R. M., Mioduszewski, A. J., & Rodríguez, L. F. 2008, ApJL, 675, L29

Martín-Doménech, R., Muñoz Caro, G. M., Bueno, J., & Goesman, F. 2014, A&A, 564, A8

McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, adass XVI, 376, 127

Mumma, M. J., & Charnley, S. B. 2011, ARA&A, 49, 471

Natta, A., Testi, L., & Randich, S. 2006, A&A, 452, 245

Okuzumi, S., Momose, M., Sirono, S.-i., Kobayashi, H., & Tanaka, H. 2016, ApJ, 821, 82

Pérez, L. M., Carpenter, J. M., Andrews, S. M., et al. 2016, Sci, 353, 1519

Prato, L., Greene, T. P., & Simon, M. 2003, ApJ, 584, 853

Rosotti, G. P., Juhasz, A., Booth, R. A., & Clarke, C. J. 2016, MNRAS, 459, 2790

Ruane, G., Mawet, D., Kastner, J., et al. 2017, AJ, 154, 73

Ruge, J. P., Flock, M., Wolf, S., et al. 2016, A&A, 590, A17

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593

Szulágyi, J., Morbidelli, A., Crida, A., & Masset, F. 2014, ApJ, 782, 65

van der Marel, N., van Dishoeck, E. F., Bruderer, S., et al. 2013, Sci, 340, 1199

van der Plas, G., Wright, C. M., Ménard, F., et al. 2017, A&A, 597, A32

Williams, J. P., & Cieza, L. A. 2011, ARA&A, 49, 67

Wolf, S., Moro-Martín, A., & D'Angelo, G. 2007, P&SS, 55, 569

Yen, H.-W., Liu, H. B., Gu, P.-G., et al. 2016, ApJL, 820, L25

Zhang, K., Blake, G. A., & Bergin, E. A. 2015, ApJL, 806, L7