# A Resolved and Asymmetric Ring of PAHs within the Young Circumstellar Disk of IRS 48

Guillaume Schworer, Sylvestre Lacour, Nuria Huélamo, Christophe Pinte, Gael Chauvin, Vincent Coudé Du Foresto, David Ehrenreich, Julien Girard, Peter Tuthill

► To cite this version:

Guillaume Schworer, Sylvestre Lacour, Nuria Huélamo, Christophe Pinte, Gael Chauvin, et al. A Resolved and Asymmetric Ring of PAHs within the Young Circumstellar Disk of IRS 48. The Astrophysical Journal, 2017, 842 (2), pp.77. 10.3847/1538-4357/aa74b7. obspm-02190411

HAL Id: obspm-02190411

https://hal-obspm.ccsd.cnrs.fr/obspm-02190411v1

Submitted on 27 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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A Resolved and Asymmetric Ring of PAHs within the Young Circumstellar Disk of IRS 48

Guillaume Schworer¹,², Sylvestre Lacour¹,³, Nuria Huélamo⁴, Christophe Pinte⁵,⁶, Gaël Chauvin⁵, Vincent Coudé du Foresto¹, David Ehrenreich⁷, Julien Girard⁸, and Peter Tuthill²

¹ LESIA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Université, UPMC Univ. Paris 06, Univ. Paris Diderot, Sorbonne Paris Cité, France

²Sydney Institute for Astronomy, School of Physics, The University of Sydney, NSW 2006, Australia

³Cavendish Laboratory, University of Cambridge, JJ Thomson Avenue, Cambridge CB3 0HE, UK

⁴Dpto. Astrofísica, Centro de Astrobiología (INTA-CSIC), ESAC Campus, P.O. Box 78, E-28691, Villanueva de la Cañada, Spain

⁵Univ. Grenoble Alpes, IPAG, F-38000 Grenoble, France CNRS, IPAG, F-38000 Grenoble, France

⁶UMI-FCA, CNRS/INSU France (UMI 3386), and Departamento de Astronomía, Universidad de Chile, Casilla 36-D Santiago, Chile

⁷Observatoire de l'Université de Genève, 51 chemin des Maillettes, 1290 Versoux, Switzerland

⁸European Southern Observatory, Alonso de Cordova 3107, Casilla 19001 Vitacura, Santiago 19, Chile

Received 2016 September 28; revised 2017 May 4; accepted 2017 May 20; published 2017 June 15

## Abstract

For one decade, the spectral type and age of the ρ Oph object IRS-48 were subject to debate and mystery. Modeling its disk with mid-infrared to millimeter observations led to various explanations to account for the complex intricacy of dust holes and gas-depleted regions. We present multi-epoch high-angular-resolution interferometric near-infrared data of spatially resolved emissions in the first 15 au of IRS-48, known to have very strong polycyclic aromatic hydrocarbon (PAH) emissions within this dust-depleted region. We make use of new Sparse-Aperture-Masking data to instruct a revised radiative-transfer model, where spectral energy distribution fluxes and interferometry are jointly fitted. Neutral and ionized PAH, very small grains (VSG), and classical silicates are incorporated into the model; new stellar parameters and extinction laws are explored. A bright (42 $L_{⊙}$) and hence large (2.5 $R_{⊙}$) central star with $A_{V}$ = 12.5 mag and $R_{V}$ = 6.5 requires less near-infrared excess: the inner-most disk at ≈1 au is incompatible with the interferometric data. The revised stellar parameters place this system on a 4 Myr evolutionary track, four times younger than the previous estimations, which is in better agreement with the surrounding ρ Oph region and disk-lifetime observations. The disk-structure solution converges to a classical-grain outer disk from 55 au combined with an unsettled and fully resolved VSG and PAH ring, between 11 and 26 au. We find two overluminosities in the PAH ring at color-temperatures consistent with the radiative transfer simulations; one follows a Keplerian circular orbit at 14 au. We show a depletion of a factor of ≈5 of classical dust grains up to 0.3 mm compared to very small particles: the IRS-48 disk is nearly void of dust grains in the first 55 au. A 3.5 $M_{Jup}$ planet on a 40 au orbit can qualitatively explain the new disk structure.

**Key words:** planetary systems – protoplanetary disks – radiative transfer – stars: individual (IRS-48) – stars: pre-main sequence

## 1. Introduction

The dust- and gas-rich disks surrounding numerous pre-main-sequence (PMS) stars are of key interest for unveiling how planetary systems are formed; they are the initial conditions for planetary formation. Protoplanetary disks have a rich structure, with different physics playing a role in different regions of the disk. The dynamic ranges involved span two to five orders of magnitude on spatial scales, orbital times, temperatures, and much more in dust or gas densities. The extreme dynamic ranges involved in the structure and composition of these objects mean that very different observational techniques have to be combined together to probe their various regions.

While spatially resolved observations give a direct localization of emitting regions, allowing precise disambiguation of structures in the spatial domain, a somewhat similar approach exists in the spectral domain through the analysis of the relative brightness distribution of parts (or all of) the disk. The discovery of a medium-age sub-class of young circumstellar disks, known as transition disks, characterized by a distinctive dip in their infrared spectral energy distribution (SED) suggests that a partially evacuated gap exists in the inner region of the protoplanetary disk (e.g., Calvet et al. 2002). The profound

implications for studies of planetary formation have become increasingly apparent with the confirmation of the disk-gap architecture by sub-millimeter measurements and optical interferometry (e.g., Andrews et al. 2011).

IRS-48 is a spectacular candidate for which the delicate balance of such phenomena is currently questioned. Also referred as 2MASS J16273718-2430350, *GY* 304, or *WLY* 2-48, it is a highly extinct $A_{O}^{+4}_{-1}$ star, part of the ρ Oph cloud L1688; Table 1 sums up its main data.

First listed in Wilking et al. (1989), IRS-48 shows a substantial far-infrared and millimeter excess, pointing toward a classification as a transition disk. Since the first resolved images obtained by Geers et al. (2007a) in the near- and mid-infrared, and the first substantial evidence of Polycyclic Aromatic Hydrocarbon (PAH) emission lines in Geers et al. (2007b), IRS-48 has been an object of interest and mystery. While 8.6 through 11.3 µm images showed high and unresolved PAH fluxes, the 18.7 µm image uncovered a purely thermal and asymmetric outer disk from a radius of 55 au. Using the VLT-CRIRES spectrograph to observe the 4.7 µm CO fundamental rovibrational band, Brown et al. (2012a) imaged a 30 au thin ring. Additionally, they solve the spectral type puzzle of the central star and set a luminosity of 14.3 $L_{⊙}$.
<table><caption>**Table 1**
Main Literature Data of ρ Oph IRS-48</caption><thead><tr><th>Parameter</th><th>Value</th><th>Reference</th></tr></thead><tbody><tr><td>Stellar position<br>(J2000)</td><td>16<sup>h</sup>27<sup>m</sup>37<sup>s</sup> 18<br>-24°30'35"3</td><td>B14</td></tr><tr><td>Distance</td><td>120 pc</td><td>L08</td></tr><tr><td>Inclination</td><td>i = 50°</td><td>G07, B14</td></tr><tr><td>Systemic velocity</td><td>4.55 km s<sup>-1</sup></td><td>vdM13</td></tr><tr><td>Position angle</td><td>PA = 96°</td><td>G07, B14</td></tr><tr><td>Stellar type</td><td>A0<sup>+4</sup><sub>-1</sub></td><td>B12</td></tr><tr><td>Stellar Temp.</td><td>9000 ± 550 K</td><td>B12</td></tr></tbody></table>

**Note.** B14—Bruderer et al. (2014), L08—Loinard et al. (2008), B12—Brown et al. (2012a), vdM13—van der Marel et al. (2013).

combined with an extinction of $A_v = 11.5$, positioning the object on a—suspiciously old—15 Myr evolutionary track.

Subsequent ALMA and VLA observations at 0.44, 1.3, and 8.8 mm carried out by van der Marel et al. (2013, 2016) unveiled a millimeter-grain asymmetry at a radius of 63 au, in the southern extents of the disk. They postulate a planet of $10 M_{Jup}$ located at ≈18 au triggering a vortex-shaped dust trap. Further work on PAHs by Maaskant et al. (2014), using unresolved spectral energy distribution data, revealed that the near- and mid-infrared spectrum of IRS-48 is dominated by a mixture of neutral and ionized PAHs, which they postulate to arise from an extended region between the inner-most disk at 1 au to the outer disk at 55 au.

Most recent studies by Bruderer et al. (2014) showed that a single gas depletion in the inner disk could not explain the flux profiles of $^{12}$CO and C$^{17}$O lines measured by ALMA with a resolution of 30 au. They proposed a slight gas depletion at ≈20–50 au in addition to a complete depletion inside 20 au. Follette et al. (2015) report the first reflected light H- and Ks-band direct images of the outer 55 au disk, and question the abundance of the inner-most disk at ≈1 au using SED modeling, though their fit suffers from over-luminous silicate features in the 9–18 µm window.

A number of mechanisms have been proposed to cause gaps in protoplanetary disks, including extensive grain growth (Dullemond & Dominik 2005), photo-evaporation (Clarke et al. 2001), binarity (Ireland & Kraus 2008), and a tidal barrier created by dynamical interaction with low-mass disk objects (e.g., Bryden et al. 1999). These different mechanisms can be distinguished by studying the distribution of the gas and dust within the gaps: a large (stellar) companion or photo-evaporation would almost completely evacuate the inner regions while a less massive planetary companion would allow gas and small dust grains to exist within its orbit (Lubow et al. 1999). Furthermore, in the latter case, the measurement of the size and distribution of this material would allow the orbit and mass of the planetary companion to be constrained.

Key grain-coagulation theories (see Williams & Cieza 2011 for a review) are able to explain growth from submicrons to millimeters, and then hectometer to planets. Several possible solutions have been proposed to account for the formation of bodies up to meter-sizes, such as turbulent vortices (Heng & Kenyon 2010). Another mystery regarding grain-coagulation exists for the smallest grains, that can arguably be considered to be large molecules: PAHs, and in a lesser extent, very small grains (VSGs). It has been argued that PAHs may either take

<table><caption>**Table 2**
Observations Made with NaCo</caption><thead><tr><th>Date</th><th>Band</th><th>Calibrators</th><th>Seq.</th></tr></thead><tbody><tr><td>2011 Mar 14</td><td>Ks<sup>a,b</sup></td><td>Elia 2–35 and 2–37</td><td>4</td></tr><tr><td>2011 Mar 14</td><td>Lp<sup>a,b,c</sup></td><td>Elia 2–37<sup>a,b</sup><br>Hip 86311<sup>c</sup></td><td>4</td></tr><tr><td>2011 Mar 14</td><td>Mp<sup>a,b</sup></td><td>Elia 2–15</td><td>2</td></tr><tr><td>2011 Sep 1</td><td>Lp<sup>a</sup></td><td>Elia 2–37</td><td>2</td></tr><tr><td>2012 Mar 6</td><td>Lp<sup>a</sup></td><td>Elia 2–37</td><td>3</td></tr><tr><td>2013 Mar 25</td><td>Lp<sup>a</sup></td><td>Elia 2–37 and 2–11</td><td>4</td></tr></tbody></table>

**Notes.** A calibration sequence (“Seq.” in the headline, also called “bracket”) is a unitary interferometric measurement following the observation pattern “Calibrator–Science Target–Calibrator” Single-frame exposition time is 0.4 s for both SAM and imaging modes.

a SAM.

b Full-pupil image.

c Spectra.

part in the dust coagulation process (Dullemond et al. 2007) or are replenished by mixing processes in the disk (Siebenmorgen & Krügel 2010). However, observational evidence of the role of PAHs in these processes is not conclusive and it is presently unclear how important these routes are. The question remains whether these very small particles (VSPs) take part in grain coating, in the main grain-growth process (and their presence in many disks is to be explained by some replenishment processes), or if their intrinsic properties prevent them from growing (and their absence in many disks is to be explained by some other assimilation processes). Gorti & Hollenbach (2008) find that the abundance of PAHs and small grains in disks is critical to the temperature profile. Indeed, they show that if some PAHs are present, their contribution to disk heating via grain photoelectric emission can be larger by a factor of ≈2 than X-ray heating at $r \approx 10$ au.

This paper reports the first direct detection of the full extent of a PAH ring in a young circumstellar disk, presents a revised model for the IRS-48 object to explain the rich and complex dust- and gas-environments observed from near-infrared to centimeter wavelengths, and sets limits on how many silicates grains—hence replenishment—are to be expected in the PAH and VSG rings.

## 2. Observations and Data Reduction

### 2.1. Observation Strategy

In this work, we make use of new data sets acquired using the instrument NaCo (abbreviation for the Nasmyth Adaptive Optics System (NAOS) and the Near-infrared Imager and Spectrograph (CONICA)) commissioned on the Very Large Telescope (VLT) at the Paranal Observatory of the European Southern Observatory (ESO). This data was acquired in different near-infrared bands and at four epochs over a two-year period; these sets are detailed in Table 2. At that time, NaCo was integrated at the Nasmyth platform of the UT4 8.2 m telescope. We used this instrument in Sparse-Aperture Masking (SAM) and full-pupil imaging modes. For both modes, given the large reddening of the target (>16 mag in the visible), we use the NaCo built-in infrared wavefront sensor for the Adaptive Optics (AO) correction.
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Full-pupil image in Lp-band (top left), and PSF-subtracted (top right) in normalized flux (peak of the un-subtracted image). The 86 mas slit height is overlaid in red, while the black dashed lines lying at ≈115 mas illustrate the first null of the central star, at 1.22λ/D. The displayed spectra show the disk-spectra taken at the east and west stellar-nulls, divided by the spectra taken at the peak intensity. Prominent PAH emission is seen from the disk at 3.3 µm, and Brα from the star at 4.05 µm. The same east-west asymmetry is seen as the one observed on the image.</figcaption></figure>

## 2.2. Direct Imaging and Spectroscopy

At epoch 1, full-pupil images were acquired in Ks-, Lp-, and Mp-bands, as well as spectra between 2.87 and 4.11 µm, with $R_{average} \approx 800$. The Lp-band image is displayed in Figure 1 both with and without point-spread function (PSF) subtraction. One can easily see the east-west elongated nature of the circumstellar emission, at ≈100 milli-arcsecond (mas), as well as an asymmetry in the disk semimajor axis; the western extension is ≈1.15 times brighter than the eastern one. However, the ≈100 mas spatial scale is similar to the resolution criterion $\lambda/D = 96$ mas at these wavelengths, which enormously complexifies the use of these images.

The slit of the spectrograph was orientated along the east-west direction as displayed in the latter figure. The pixel scale of the detector is constant over the whole spectra: 1.6 Å × 2.7 mas. The slit height is 86 mas; given the short scale extents of the disk and the high inclination (≈50°), we expect a large contamination of the disk at the stellar position (assumed to be located at the peak intensity in the spectra). Similarly, given the full width at half maximum of the PSF at ≈50 mas, we expect a large contamination of the star in the disk. To minimize the stellar contamination on the disk, we make use of the spectra at the first null of the star, 1.22λ/D ≈ 115 mas. The PSF-reference image shows a very low stellar flux in the null,

≈3% of the peak flux; this highlights the high strehl ratio of the observation.

The spectra in Figure 1 show the flux of the disk—taken at the first stellar null—relative to the peak flux of the spectra. The extraction aperture is three pixels for the peak, and one pixel interpolated at 1.22λ/D for the null. One can see very strong PAH emission features at 3.3 µm, characteristic of neutral PAH: in this emission line, the disk-to-peak contrast jumps by nearly a factor of two, compared to the disk-to-peak continuum contrast. Note that this emission line is not recorded in the interferometric data or full-pupil images given that it is outside the Lp-band filter bandpass, see Figure 1. Although both east and west spectra are very similar in shape, the Western disk-extension is constantly—at all spectrum wavelengths—1.15 times brighter than its eastern region. The same east-west asymmetry was measured on the full-pupil image.

## 2.3. Near-infrared Sparse-aperture Masking

The most interesting data sets are SAM data in that they offer an angular resolution down to a few tens of mas, equivalent to a few astronomical units given the distance of 120 pc, as shown in Table 3. Except for the longest baselines of ALMA (≥1 km), such resolution cannot be reached with any other instrument nor any other NaCo mode (coronographic, full-pupil, etc). Moreover, new extreme AO systems such as SPHERE usually
<table><caption>**Table 3**
NaCo Filter Characteristics, Interferometric $0.5 \lambda/D$ Criteria, and Field of View with Respect to Filters at Central Wavelengths, Given a Target Distance of 120 pc</caption><thead><tr><th></th><th>Ks</th><th>Lp</th><th>Mp</th></tr></thead><tbody><tr><td>Central wavelength (μm)</td><td>2.18</td><td>3.8</td><td>4.78</td></tr><tr><td>Bandwidth (μm)</td><td>0.35</td><td>0.62</td><td>0.59</td></tr><tr><td></td><td colspan="3">Smallest Baseline, 1.78 m</td></tr><tr><td>0.5 λ/D (au)</td><td>15.2</td><td>26.6</td><td>33.5</td></tr><tr><td>Field of view (au)</td><td>190</td><td>326</td><td>544</td></tr><tr><td></td><td colspan="3">Longest Baseline, 6.43 m</td></tr><tr><td>0.5 λ/D (au)</td><td>4.2</td><td>7.4</td><td>9.2</td></tr><tr><td>Field of view (au)</td><td>52.6</td><td>90.8</td><td>150</td></tr></tbody></table>

have their wavefront sensors in the visible and cannot close the AO loop on this extremely extinct target (>16 mag in the visible), which means that NaCo-SAM is the only instrument able to observe IRS-48 at such resolution in visible or in the infrared (IR).

SAM, also called Non-Redundant Aperture Masking (NRM) is an observational technique in which a single mirror is transformed into a multi-pupil Fizeau interferometer. A mask placed in the pupil plane of the telescope transforms the usual Airy PSF into a pattern made of multiple fringes. The technique consists in measuring the amplitude and phase of these interferometric fringes so that information can be retrieved down to a fraction of a resolution element of the longest baseline of the mask (typically $0.5\lambda/B_{max} \approx 50$ mas at 4 µm). However, the interferometric field of view remains very limited, in practice, of the order of $2.5\lambda/B_{min} \approx 250$ mas at 4 µm (though a somewhat lower sensitivity exists until $\approx 6\lambda/B_{min}$). Using closure phases (T3)—an estimator robust to atmospheric turbulence and optical aberrations (Jennison 1958; Monnier 2000)—allows high-contrast capability without the inevitable “blind-spot” of other techniques such as coronography.

The observation strategy for the SAM mode was designed to intertwine the science target IRS-48 between two blocks on an interferometric calibrator, here Elia 2–35, 2–37, or 2–15, as detailed in Table 2. The mask used was a seven-hole non-redundant mask with baselines ranging from 1.78 to 6.43 m, allowing the simultaneous measurement of 21 visibility-squared (VIS2) and 35 closure phases (T3, among which 15 are independent). In full-pupil imaging mode, we used the same calibrators to have PSF references.

The data was processed using the SAMP pipeline (Lacour et al. 2011). It includes sky subtraction, bad pixel subtraction and fringe fitting. T3 and VIS2 measurements are calibrated by the two calibrators for each bracket.

Calibration to PSF reference help reduce systematic errors in the data, but hardly corrects all of them. A fitting mechanism to account for that is specifically developed, refer to Section 5. The error bars on the data are obtained from statistical measurements on the data, which account for random errors. T3 errors they average $\approx 0.5$ in Lp-band and $\approx 1.7$ in Ks and M bands.

## 2.4. Spectral Energy Distribution

In addition to the multi-wavelength interferometric data listed above, we compiled SED fluxes for IRS-48 from the

<table><caption>**Table 4**
Previously Published Photometry and Spectra of $\rho$ Oph IRS-48</caption><thead><tr><th>Wavelength(s) (µm)</th><th>Instrument</th><th>Reference</th></tr></thead><tbody><tr><td>0.43, 0.64</td><td>NOMAD</td><td>Z05</td></tr><tr><td>0.65, 0.8</td><td>Hydra</td><td>E11</td></tr><tr><td>1.2, 1.6, 2.2</td><td>2MASS</td><td>C03</td></tr><tr><td>3.4, 4.6</td><td>WISE</td><td>W10</td></tr><tr><td>3.6, 4.5</td><td>Spitzer IRAC</td><td>vK09</td></tr><tr><td>5.9–36.9</td><td>Spitzer IRS</td><td>MC10</td></tr><tr><td>60–181</td><td>Herschel PACS</td><td>F13</td></tr><tr><td>450</td><td>ALMA</td><td>vdM13</td></tr><tr><td>850</td><td>SCUBA</td><td>A07</td></tr><tr><td>880, 1300</td><td>SMA</td><td>B12b</td></tr></tbody></table>

**Note.** Z05—Zacharias et al. (2005), E11—Erickson et al. (2011), C03—Cutri et al. (2003), W10—Wright et al. (2010), vK09—van Kempen et al. (2009), MC10—McClure et al. (2010), F13—Fedele et al. (2013), A07—Andrews & Williams (2007), B12b—Brown et al. (2012b), vdM13—van der Marel et al. (2013).

literature. All data are listed in Table 4 with the corresponding references. As can be seen, the photometry from the literature has a broad coverage from visible-blue to millimeter wavelengths. No measurements in the ultraviolet exist because of the very high extinction of the source; measurements in centimeter wavelengths are discarded given that they trace cold and large-scale structures outside the region of sensitivity and interest of our observations. For the fitting exercise, the errors on photometry fluxes were assumed to be of 5% if not specified in the literature.

The correction of the interstellar extinction was made using Cardelli et al. (1989), improved by O’Donnell’s (1994) work in the visible. These extinction curves are parametrized with $R_v$ (reddening slope) and $A_v$ (absorption in V band) such that $A_v \equiv R_v \cdot \text{Extinction}(B_{mag}-V_{mag})$. The initial values are taken from Brown et al. (2012a), $A_v = 11.5$ and $R_v = 5.5$.

# 3. Radiative Transfer Model

## 3.1. Old Model

We take over the disk structure presented in Bruderer et al. (2014). Given the prominence of PAH in the SED of IRS-48, we merge into this disk model the results of Maaskant et al. (2014) on ionized PAH in several young disks. We obtain a disk structure in three parts: (1) an inner-most dusty disk between 0.4 and 1 au, (2) an outer dusty disk between 55 and 160 au, and (3) a PAH disk in between. For the sake of simplicity, we call this disk model the old model. Its main parameters can be seen in Table 5. This model represents the literature knowledge of IRS-48.

We use MCFOST, a parallelized 3D radiative transfer code based on a Monte-Carlo method (Pinte et al. 2006), able to output SED data and monochromatic images when given an input parameter file. After a few minor adjustments, the fit to the SED literature data is remarkable. However, this model or any small variation of it fails to explain the VIS2 data in either the Lp- or Mp-band. This highlights the degeneracy of SED-fitting exercises.

Indeed, the inner dust disk between 0.4 and 1 au emits very strongly in NIR, while the PAH grains further out are much dimmer. This creates a very high unresolved flux (inner disk and star) compared to the resolved flux from the PAH in the
<table><caption>**Table 5**
Model Parameters of the Radiative Transfer of IRS-48, Comparison between the Best-fit for the New Model and the Old Model as a Merge of Bruderer et al. (2014) and Maaskant et al. (2014)</caption><thead><tr><th>Parameter</th><th>Old Model</th><th>New Model</th></tr></thead><tbody><tr><td></td><td colspan="2" style="text-align: center;">Stellar parameters</td></tr><tr><td>Temperature</td><td>9250 K</td><td>9250 K<sup>a</sup></td></tr><tr><td>Luminosity</td><td>14.3 L<sub>☉</sub></td><td>42 L<sub>☉</sub></td></tr><tr><td>Distance</td><td>120 pc</td><td>120 pc<sup>a</sup></td></tr><tr><td>R<sub>v</sub></td><td>5.5</td><td>6.5</td></tr><tr><td>A<sub>v</sub></td><td>11.5</td><td>12.9</td></tr><tr><td></td><td colspan="2" style="text-align: center;">Global disks parameters</td></tr><tr><td>Inclination</td><td>50°</td><td>50°<sup>a</sup></td></tr><tr><td>PA</td><td>96°</td><td>96°<sup>a</sup></td></tr><tr><td>Scale height</td><td>10 au</td><td>14 au</td></tr><tr><td>Ref. radius</td><td>60 au</td><td>60 au<sup>a</sup></td></tr><tr><td></td><td colspan="2" style="text-align: center;">Inner-most disk</td></tr><tr><td>R<sub>in</sub> Inner</td><td>0.4 au</td><td>...</td></tr><tr><td>R<sub>out</sub> Inner</td><td>1 au</td><td>...</td></tr><tr><td>Flaring exp. β</td><td>1.3</td><td>...</td></tr><tr><td>Surface density exp p</td><td>-1</td><td>...</td></tr><tr><td>Grain-sizepower law</td><td>-3.5</td><td>...</td></tr><tr><td>Dust Mass</td><td>8.0 × 10<sup>-12</sup> M<sub>☉</sub></td><td>...</td></tr><tr><td>Silicate grains Mass</td><td>1%</td><td>...</td></tr><tr><td>Carb. grains Mass</td><td>99%</td><td>...</td></tr><tr><td>Grain sizes a</td><td>0.03–30 µm</td><td>...</td></tr><tr><td></td><td colspan="2" style="text-align: center;">VSP ring</td></tr><tr><td>R<sub>in</sub> Inner</td><td>1 au</td><td>11 au</td></tr><tr><td>R<sub>out</sub> Inner</td><td>50 au</td><td>26 au</td></tr><tr><td>Flaring exp. β</td><td>1.3</td><td>0.6</td></tr><tr><td>Surface density exp p</td><td>-1</td><td>-0.1</td></tr><tr><td>Grain-sizepower law</td><td>...</td><td>-4<sup>a</sup></td></tr><tr><td>Dust Mass</td><td>8.0 × 10<sup>-10</sup> M<sub>☉</sub></td><td>3.7 × 10<sup>-10</sup> M<sub>☉</sub></td></tr><tr><td>VSG Mass</td><td>0%</td><td>20%</td></tr><tr><td>VSG sizes</td><td>...</td><td>10–30 Å<sup>a</sup></td></tr><tr><td>Neutral PAH Mass</td><td>50%</td><td>60%</td></tr><tr><td>Ionized PAH Mass</td><td>50%</td><td>20%</td></tr><tr><td>PAH sizes</td><td>5 Å</td><td>4–10 Å<sup>a</sup></td></tr><tr><td></td><td colspan="2" style="text-align: center;">Outer disk</td></tr><tr><td>R<sub>in</sub></td><td>55 au</td><td>55 au<sup>a</sup></td></tr><tr><td>R<sub>out</sub></td><td>160 au</td><td>250 au</td></tr><tr><td>Flaring exp. β</td><td>1.3</td><td>0.67<sup>a</sup></td></tr><tr><td>Surface density exp p</td><td>-1</td><td>-1<sup>a</sup></td></tr><tr><td>Grain-sizepower law</td><td>-3.5</td><td>-3.5<sup>a</sup></td></tr><tr><td>Dust Mass</td><td>1.0 × 10<sup>-5</sup> M<sub>☉</sub></td><td>9.0 × 10<sup>-6</sup> M<sub>☉</sub></td></tr><tr><td>Silicate grains Mass</td><td>70%</td><td>70%<sup>a</sup></td></tr><tr><td>Carb. grains Mass</td><td>30%</td><td>30%<sup>a</sup></td></tr><tr><td>Grain sizes a</td><td>0.03–4000 µm</td><td>0.03–4000 µm<sup>a</sup></td></tr></tbody></table>

**Note.**

≈9–55 au region. Hence, the inner regions of the old model disk appears spatially too small, mostly unresolved and the resulting visibilities are close to unity: much larger than the VIS2 data at all baselines. Given the difference of the VIS2 values, the resolved-flux depletion is estimated as ≈2.7 for baselines ≳3 m.

The only solution to make the old model fit the VIS2 is to incorporate a very high flux in the ≈9–35 au region—and in this region only—to make up for the contrast ratio of ≈2 in Lp-band between unresolved and resolved flux.

This implies fundamentally changing the structure, and hence the nature, of the disk around IRS-48.

## 3.2. New Model: Fitting Strategy for VIS2

In our data, VIS2 measurements are not excessively sensitive to asymmetries, essentially because of their relatively large error bars. Hence, they encode a point-symmetrical and smooth disk structure and inform on high-level morphology. We make use of the VIS2 data in Lp-band from all epochs and the Mp-band taken at epoch 1, as well as the SED fluxes from the literature (refer to Table 4). Ks-band from epoch 1 are discarded given that their large uncertainties would not set meaningful constraints on the disk structure. The Lp-band VIS2 data does not show time-variation along the four epochs.

We decided to perform a manual initial exploration of the parameters to identify the most impacting ones. We refine the fit later on with both a genetic algorithm and low-dimensional grids. The genetic algorithm allows for a very fast convergence for cases with a large number of dimensions where basic Markov chain Monte Carlo (MCMC) would fail. Grid calculations only come in a second time and ensure a full exploration of a small subset of the entire parameter space. Verifying smoothness for the outputs of the modeling over these small grid-subsets allowed us to interpolate the results and thus minimize the processing time for these grids. The final convergence was achieved by eye.

For the purpose of measuring visibilities on these MCFOST images, we choose their pixel resolution in “mas per pixel” to display the full interferometric field of view in a constant width of 1001 pixels. Also, in order to account for bandwidth smearing, we generate 15 monochromatic images linearly spaced in wavelength over the bandwidth of each filter (see Table 3) and combine them into a unique polychromatic image.

We perform an SED-fitting “inside-out,” or “in increasing wavelengths,” so that constraints and parameters obtained from fitting interior structures can remain (mostly) unchanged when fitting structures located further away from the star.

## 3.3. New Model: Initial Configuration and Assumptions

We assume an axial-symmetry perpendicular to the plane of the disk as well as plane symmetry above and below the plane of the disk. We adopt the standard flared disk prescription (Shakura & Sunyaev 1973), described in cylindrical coordinates (r, z) such that the density profile is given by

$$ \sum(r, z) = \sum_{0} r^{p} \exp\left(-\frac{1}{2}\left(\frac{z}{h}\right)^{2}\right), \quad (1) $$

where $\sum_0$ is a density-normalization constant, p is the surface density exponent, and h is the disk scale height; h varies with radius as $h = h_0(r/r_0)^\beta$, where $\beta$ is the flaring exponent ($\beta > 0$, of order ≈1) and $h_0$ is the scale height at radius $r_0$.

The SED shape shows that IRS-48 consists of—at least—two main disk-components: an outer and extended disk beyond ≈55 au radius, and a closer one, that emits most of the PAH light, as found by Maaskant et al. (2014). Our VIS2 data does not constrain the outer disk very well, which mostly lies outside its field of view. Hence, we systematically adopt the simplest axisymmetrical fitting solution consistent with the literature for this outer disk, while we focus our fitting effort on the inner parts.

<footnote><sup>a</sup> Set-value parameter.</footnote>
<table><caption>**Table 6**
New Stellar Model for ρ Oph IRS-48 Central Star</caption><thead><tr><th></th><th>This Work</th><th>F15</th><th>B12</th></tr></thead><tbody><tr><td>T (K)</td><td>9250</td><td>9500</td><td>9000</td></tr><tr><td>R (R<sub>☉</sub>)</td><td>2.5<sup>a</sup></td><td>1.8</td><td>1.4</td></tr><tr><td>L (L<sub>☉</sub>)</td><td>42</td><td>23.6</td><td>14.3</td></tr><tr><td>Age (Myr)</td><td>4.1<sup>a</sup></td><td>8</td><td>&gt;15</td></tr><tr><td>M (M<sub>☉</sub>)</td><td>2.5<sup>a</sup></td><td>2.2</td><td>2.0</td></tr><tr><td>A<sub>v</sub></td><td>12.9</td><td>12.0</td><td>11.5</td></tr><tr><td>R<sub>v</sub></td><td>6.5</td><td>(4.0)</td><td>5.5</td></tr><tr><td>A<sub>k</sub></td><td>(1.8)</td><td>1.5</td><td>(1.6)</td></tr></tbody></table>

Notes. F15—Follette et al. (2015), B12—Brown et al. (2012a).

<sup>a</sup> According to Siess et al. (2000) YSO evolutionary models.

Similarly to the old model, we adopt a three-disk-component model: two dusty inner- and outer-disks and an additional disk composed of VSP in between. For each of them, the inner radius, outer radius, and total dust mass are free parameters. As an initial step, the parameters “scale height” and “flaring exponent” are shared between all three disks. Other less-impacting parameters, such as the surface density exponent, the grain minimum and maximum sizes, the ratio between carbonaceous (from Li & Greenberg 1997) and silicate (from Draine & Lee 1984) grains, the ratio between ionized and neutral PAH, and the grain-size power-law index or the disk inclination are left-aside for the initial fitting exercise and set to standard values from IRS-48 and disk literature. The flaring exponent = 0.67 (van der Marel et al. 2013), the inclination = 50° (Geers et al. 2007a; van der Marel et al. 2013; Bruderer et al. 2014), the surface density exponent = -1 (Andrews et al. 2013), the inner-rim radius = 55 au (Geers et al. 2007a), the silicate over carbonaceous grains mass ratio = 70/30 (prescriptions span 85/15 in Wood et al. 2002; 45/55 in Kim et al. 1994; 80/20 in Maaskant et al. 2014) and the min/max grain sizes = [0.03, 4000] µm (van der Marel et al. 2013) of the outer disk are set parameters, see Table 5.

The effective temperature of the central star is set to 9250 K, midway between Follette et al.’s (2015) work (9500 K) and Brown et al.’s (2012a) work (9000 K), who also estimated its spectral type to be A0$^{+4}_{-1}$ (9000 ± 550 K), see Table 6.

### 3.4. New Model: Inner-most Region and Star

VIS2 data gives stringent constraints on both the flux ratio between the unresolved flux (inner-most disk and star) and resolved flux (further than ≈9 au). To correct for the resolved-flux depletion of about a factor of 2.7 compared to the unresolved flux, one can either (1) decrease the inner-most disk opacity to decrease the unresolved flux or (2) increase the PAH disk opacity further than 9 au to increase the resolved flux. This later solution—increasing PAH emission in Lp band (3.8 µm) and Mp band (4.78 µm)—would immediately lead to a dramatic over-brightness of the main PAH emission features between 4 and 15 µm, with no (or few) degrees of freedom to compensate for it (i.e., structural parameters such as flaring exponent, scale height, etc).

Decreasing the disk brightness within the resolution angle of the instrument would automatically lead to an under-fitting of the NIR flux. This can easily be corrected by increasing the brightness of the star, de facto decreasing the amount of NIR excess arising from dust needed in the fit. Indeed, as discussed by Follette et al. (2015), a hotter and especially larger star questions the existence of a third inner-most disk in the ≈1 au

area from the star where this disk was usually needed to account for the 1–3 µm IR excess. We hence decided to release fitting constraints by allowing the stellar luminosity to vary between 10 and 100 $L_{☉}$. This implied adding two additional free parameters to tune the extinction of the target in order to fit the photometric fluxes measured in the visible and NIR ($R_{v}$ taken in [4.0–10.0] and $A_{v}$ taken in [10.0–15.0]).

### 3.5. New Model: Very Small Particle Disk

A morphological fit of a simple two-dimensional Gaussian shows that the characteristic location of the resolved flux should lie in a region between 10 and 20 au. The PAH disk is a natural candidate to account for this flux, since previous work showed that the cavity up to ≈55 au was not depleted of such grains (Geers et al. 2007a; Maaskant et al. 2014). However, this means that the bright inner rim of the PAH disk must be fully resolved, i.e., further than ≈7–10 au, in order to account favorably for the contrast ratio between unresolved and resolved flux.

The MIR SED does not show any bump peeking at 9.5 and 18.3 µm, characteristic of silicate grains smaller than ≈2 µm. It does instead show a clear deep at 9.8 µm, de facto eliminating any major silicate grains contribution to the MIR flux; refer to Section 4.5 for a study on CTG in the inner regions. Hence, we decide not to incorporate any small or classical grains in the cavity between the inner and outer disks.

To model the VSP, we follow the prescription of Natta et al. (1993) and use a mixture of carbonaceous VSG and PAH with an ionized/neutral fraction parameter, from Li & Draine (2001), Weingartner & Draine (2001), and Draine & Li (2007). Following the classification of Tielens (2008; Table 2), we set $a_{min} = 4$ Å and $a_{max} = 10$ Å for PAH (ionized and neutral), $a_{min} = 10$ Å and $a_{max} = 30$ Å for VSG (referred to as PAH clusters and VSG in their table). The grain size power-law index for PAH and VSG is set to -4; only the relative abundances of these three species are left as free parameters for fitting. Finally, we assume that these VSP are quantum heated by UV stellar-flux, i.e., outside local thermodynamic equilibrium. The quantum heating routines of MCFOST were benchmarked and validated by Camps et al. (2015).

## 4. Resulting Structure

The best solution is obtained for a scenario where the inner-most disk (0.4–1 au) has such a low mass that it does not contribute to IRS-48 flux in any wavelength. The new IRS-48 model is hence composed of two disk elements: a VSP-only disk between 11 and 26 au and an outer disk from 55 au, as portrayed in the Lp-band model image (see Figure 2). The star has a higher luminosity $L_{*} = 42 L_{☉}$, leading to a larger radius $R_{*} = 2.5 R_{☉}$ than previously reported. The main parameters of the solution are listed in Table 5. Note again that the VIS2 data is mostly blind to the outer disk; the solution found represents the simplest centro-symmetrical structure that satisfied the SED fluxes.

### 4.1. General Results

Figure 4 shows the SED for both models with their respective interstellar extinction. Our new model provides the same fit quality as the old model. It matches very closely all the main features of the SED: the photosphere of the star, the NIR excess, the PAH emission features, and the outer disk. Most importantly,
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Lp-band model image normalized to the stellar flux, obtained from the radiative transfer on our model. This image was obtained out of 15 sub-images generated at wavelengths linearly spread in the Lp bandpass. The inner disk between 11 and 26 au appears to be much brighter in the Lp band than in the outer disk from 55 au. Two cavities appear, one inside at 11 au and another one between 26 and 55 au.</figcaption></figure>

it dramatically improves the fit on the VIS2 data, which it explains entirely: Figure 3 shows the comparison between the VIS2 data in Ks, Lp, and Mp bands, and both models.

## 4.2. Central Star

The solution of a $L_* = 42 L_\odot$ star with a $3.7 \times 10^{-10} M_\odot$ mass VSP ring and no inner-most disk (0.4–1 au) is a complex and precise balance between the fit to the SED and to the VIS2 data. Indeed, the NIR SED requires either a brighter star than the old model, or a similar star with an inner-most disk to account for a large NIR excess. The VIS2 data, however, pinpoints a precise contrast ratio between the resolved VSP ring and the unresolved star and inner-most disk altogether. Adding even the smallest inner-most disk mass creates an excessive NIR emission, which must be compensated with a large increase of the VSP disk-mass to fit the VIS2 data. Such an increase, however, dramatically over-estimates the PAH emission features between 5.5 and 18 µm. This strongly pushes to a solution with no inner-most disk inside the VSP ring.

We follow the isochrones of Siess et al. (2000) young stellar object (YSO) evolutionary models with a luminosity increased from 14 to 42 $L_\odot$ and an effective temperature of 9250 K to find a larger stellar radius $R_* = 2.5 R_\odot$, rather than 1.8 or 1.4 $R_\odot$ as previously reported. The main stellar parameters are listed in Table 6, with a comparison to the previous literature values.

To compensate for the excess stellar luminosity, interstellar extinction was re-evaluated to $A_v = 12.9$ and $R_v = 6.5$. IRS-48 is 1.4 V-magnitude more extinct, with a more gray $R_v$ than the extinction inferred in Brown et al. (2012a; $A_v = 11.5$ and $R_v = 5.5$). This confirms and reinforces the presence of substantial large dust grains in the line of sight (Kim et al. 1994; Indebetouw et al. 2005), typical of star-forming regions such as Ophiuchus. An equivalent value of $A_k = 1.8$ is consistent with that of ≈1.6 read on the low spatial-resolution extinction map computed by Lombardi et al. (2008); their estimate for that Ophiuchus-core region suffers, however, from a high discrepancy, most probably due to the unknown depth of their targets into this dense molecular cloud.

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Diagram showing the Ks-band (top), Lp-band (middle), and Mp-band (bottom) VIS2 data at all epochs, the calculated squared visibility from the old model image (green) and those of our model (blue) as a function of the unprojected baselines in meters. Most residuals for our model span ±1, while they exceed +4 for the old model. Although the Ks-band data is much noisier, our model still does provide an improvement over the old model.</figcaption></figure>

The radiative modeling shows that there is no extinction due to the circumstellar disk; its ≈50° inclination is not large enough to extenuate the stellar flux like (nearly) edge-on disks. All extinction is due to interstellar extinction, most of which arises from the local molecular cloud of ρ Oph.
<figure><img src="image_5.png" /><figcaption>**Figure 4.** SED of IRS-48 obtained from the old model—or the literature model—(green line) and our model (blue line). Both explain the photosphere of the star (though with different extinction parameters), the NIR excess and the outer disk SED-bump centered on 20 µm.</figcaption></figure>

The evolutionary tracks imply that such a star falls onto the ≈4 million year locus, with $M_* = 2.5 M_\odot$. This is much younger than the 8 or 15 Myr previously inferred, partly solving the evolutionary puzzle on the presence of a disk. A discussion on the age of IRS-48 is to be found in Section 6.1.

### 4.3. Very Small Particle Ring

This new model consists of a smooth centro-symmetrical and optically thin VSP-only ring between 11 and 26 au. It is responsible for all of the IR-excess until ≈13 µm, where the outer disk starts dominating the dust-emission. In the fitting exercise, the VIS2 data constrained mostly the structure of the VSP ring, while the MIR-SED PAH emission features constrained the relative abundance of the PAH grains and VSG. The VSP density tops at ≈1 × 10⁻¹⁷ kg m⁻³ (see Figure 5).

The resulting disk temperature is shown in Figure 6 for neutral PAH grains. Because of the close proximity to the star, the time between two successive absorptions is much lower than the relaxation time of the PAH grains. In this case, MCFOST treats quantum heating as “quasi-equilibrium” and a temperature can be defined (i.e., low dispersion on the temperature probability distribution). It reaches particularly high values for neutral PAH grains, $T_{max} \approx 1500$ K and $T_{min} \approx 600$ K. Ionized PAH grains and carbonaceous VSG follow a very similar temperature distribution, with lower temperatures due to their difference in electrical charge, mass, and structure: $T_{max} \approx 1100$ K and $T_{min} \approx 400$ K for ionized PAH and $T_{max} \approx 550$ K and $T_{min} \approx 250$ K for VSG.

### 4.4. Outer Disk

As previously highlighted, the interferometric data is mostly blind to the outer disk, and is not sensitive in any case to the southern asymmetry in millimeter wavelengths observed with ALMA. We fitted the simplest centro-symmetric disk model to satisfy the FIR and millimeter-wavelength SED.

<figure><img src="image_6.png" /><figcaption>**Figure 5.** Map showing the PAH dust density [kg m⁻³] in log10 scale, in a cut perpendicular to the plane of the disk. Isocontours on the map are shown at the levels of the color-bar ticks. The densest region of the VSP ring spans from 11 to 20 au with an average height of ≈5 au above the mid-plane. The ring shows a smooth and nearly flat density structure with increasing radius due to the high surface density exponent of −0.1.</figcaption></figure>

This centro-symmetric approach is, however, not void of meaning. Indeed, the outer disk is optically thin in its emission window (for wavelengths ≳13 µm, τ in the mid-plane <10⁻¹, see Figure 7). Hence the precise azimuthal location of grains (azimuthally asymmetric or not) only plays a secondary role in fitting the SED.

In order to satisfy the millimeter-wavelength SED, the outer radius of the disk is found at ≈250 au, similarly to Maaskant et al. (2014; 225 au), though SED fitting is not very sensitive to this parameter. For that disk extend, the total mass of the outer disk in the MCFOST model is 9 × 10⁻⁶ $M_\odot$ (3 $M_{Earth}$) for grains between 30 and 4 mm, with a grain-size power law of −3.5. This mass is similar, though lower than Bruderer et al.
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Map showing the temperatures (Kelvin) of the neutral PAH dust (left) and of the outer-disk grains (right), in a cut perpendicular to the plane of the disk. Iso-surfaces on the map are shown at the levels of the color-bar ticks. Compared to neutral PAH, ionized PAH grains and carbonaceous VSG follow a very similar temperature distribution, but with lower temperature bounds: $T_{max} \approx 1100$ K and $T_{min} \approx 400$ K for ionized PAH and $T_{max} \approx 550$ K and $T_{min} \approx 250$ K for VSG.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 7.** Optical thickness $\tau$ taken in the mid-plane as a function of wavelength for the VSP ring and the outer disk. One recognizes the PAH lines from the VSP ring and the two silicate bumps of the outer disk at 9.2 and 18 µm. The VSP ring becomes extremely thin for wavelengths above a micron.</figcaption></figure>

(2014; $1.6 \times 10^{-6} M_{\odot}$, $5.3 M_{\text{Earth}}$) and $9 M_{\text{Earth}}$ found by van der Marel et al. (2013) for the same grain-size distribution and bounds. This is due to the fact that the later mass estimate obtained from 685 Ghz ALMA flux was biased as they used a lower temperature, $T = 60$ K at 60 au, than what is currently observed in our model. Indeed, our star is brighter, hence it heats the outer disk more efficiently: as seen in Figure 6, the inner-rim temperature peaks at $T_{\text{dust}} = 148$ K at 55 au.

We update their dust mass calculations with the values of our model to find $\tau_{685 \text{ GHz}} = 0.14$ for $T = 140$ K at 60 au, using their expression

$$F_{\nu} = \Omega \times B_{\nu}(T_{\text{dust}}) \times (1 - \exp^{-\tau_{\nu}}). \quad (2)$$

When we integrate over the whole disk model the optical thickness at $\nu = 685$ GHz ($\lambda = 440$ µm), we find $\tau_{685 \text{ GHz}} = 0.19$, closely consistent with the previous 0.14 found independently using ALMA flux. Using $M_{\text{dust}} = 21 \times \tau_{685 \text{ GHz}} \sqrt{a_{\text{max}}}$ from

Draine (2006) and van der Marel et al. (2013) under the same hypothesis, we find $M_{\text{dust}} = 3 M_{\text{Earth}}$ for $\tau = 0.14$ and $4 M_{\text{Earth}}$ for $\tau = 0.19$, consistent with the outer disk dust mass in our radiative transfer model ($3 M_{\text{Earth}}$).

Figure 8 shows model images generated in the mid-IR with the corresponding imaging data, at 4.78 µm (Mp-band), 8.6, and 18.7 µm. Although 8.6 and 18.7 µm images were not included in the fit—only SED fluxes were fitted at these wavelengths—our model can reproduce the observed structure accurately. Geers et al. (2007a) noted that the 8.6 µm image is mostly unresolved, with most flux arising from PAH in the inner 30 au, while the 18.7 µm image is highly resolved, with most flux arising from the inner rim of the outer disk located at $\approx 55$ au.

### 4.5. CTG in the VSP Ring

Figure 9 compares intensity profiles along the semimajor axis of the observed and modeled disks at 18.7 µm. Several models were computed, where the amount of thermal grains (silicates (70%) and carbonaceous (30%)) with sizes between 0.3 and 300 µm (mass power index of -3.5), are introduced in the VSP ring as settled material (one-fifth of the VSP-ring scale height). Although our axisymmetrical model cannot obviously reproduce the eastern asymmetry, it shows a distinctive dip at the star location corresponding closely to the 18.7 µm observed image. When the amount of classical thermal grains (CTG) increases, the inner-ring becomes brighter due to the efficiency of silicate grains to radiate at 18–19 µm. The central dip progressively disappears until the flux-profile is nearly flat, for a thermal grain mass of $\approx 10 M_{\text{VSP}}$.

Note that the VIS2 data fitting was not affected by the adding of CTGs, which emit a negligible amount of flux in the Ks, Lp, and Mp bands. Indeed, the emission of such CTGs is mostly visible through the emission “bumps” of the silicate grains, at $\approx 9.5$ and 18 µm, where the SED is dominated by the outer disk. Until an addition of 5 $M_{\text{VSP}}$ worth of CTGs in the VSP ring, the fit to the SED do not show any variation. From 10 $M_{\text{VSP}}$, however, the model starts to clearly overestimate the
<figure><img src="image_9.png" /><figcaption>**Figure 8.** Synthetic-images (left and middle columns) and data images (right column) at wavelengths 4.78 µm (Mp-band, top row), 8.6 µm (middle row), and 18.7 µm (bottom row). Left-column images are shown in arbitrary flux (the central star being removed). Synthetic-images in the middle column are convolved with an 8.2 m telescope PSF and shown in the color scale of the observed data images seen in the right column. Middle and right column color scales are fluxes per square arcsec, normalized to the total flux of the disk, i.e., the fraction of the total disk surface brightness. Data images are taken from this work (Mp-band) and Geers et al. (2007a; 8.6, 18.7 µm). Isocontours are overlaid at the color-bar tick values.</figcaption></figure>

SED fluxes between ≈8 and 15 µm. Refer to Section 6.3 for a discussion.

## 5. Closure Phase Morphology

T3 are especially robust to atmospheric turbulence and optical aberrations. They are, however, in essence only sensitive to asymmetries: a point-symmetrical structure would lead to null T3 signal. As a consequence, an equal-brightness binary star would not produce any T3 signal either.

In the data, T3 at all epochs vary between −12° and +9° in the Lp band, −18° and +15° in the Mp band, and slightly less in the Ks band. This highlights the fact that strong asymmetries exist in the IRS-48 disk at all wavelengths and epochs. T3 calculated on the radiative transfer model image in all bands merely reach 0.5; a disk alone does not explain the T3 data.

### 5.1. Modeling

We fit N point sources to the T3 data in order to highlight the main locations of the asymmetries in the disk. Such simplistic models allow analytical, and hence fast, fitting to the data.

We make use of six analytical models: from one to three point sources in addition to the central star, either with or without the disk modeled with MCFOST. The central star is set as a reference to a flux of 1 and in the middle of the field. Each point source is described with three parameters: (1) angular separation to the central star within [20, 300] mas, (2) position angle (east of north) within [0, 360]°, (3) relative flux within $[\frac{1}{200}, 1]$ stellar flux.

The disk is scaled and has its flux normalized with respect to that of the central star. It is treated as a “background image” on which T3 are measured at the (U, V) coordinates of the data. Complex visibilities from each unitary model object are added
<figure><img src="image_10.png" /><figcaption>**Figure 9.** Cuts along the semimajor axis of the data image (thick black line) at 18.7 µm and of the synthetic images at the same wavelength and computed for an increasing amount of classical thermal grains (CTGs) in the PAH ring (0, 1, 5, and 10 $M_{VSP}$ worth of CTGs). Each flux profile is normalized to the total flux of the disk in each image.</figcaption></figure>

together and T3 values for the global model are obtained from the phases of the complex visibilities at the (U, V) coordinates of each T3 triangle.

The complex visibility $V$ at the (U, V) coordinates for N point sources is expressed by

$$ \left(1 + \sum_{i=1}^{N} f_i\right) V = 1 + \sum_{i=1}^{N} f_i \exp\left(-j \frac{2\pi}{\lambda} (U \delta_i + V \alpha_i)\right), \quad (3) $$

where $\lambda$ is the central wavelength of the band. $\delta_i$, $\alpha_i$ are the angular separation and position angle of the ith point source with respect to the central star (projected along declination and right ascension axes), and $f_i$ is the flux ratio.

The error estimation of interferometric data is extremely challenging, due to the intrinsic nature of the measurements that go through heavy data reduction. Systematic errors are calibrated using several PSF references (thanks to the observation pattern “Calibrator–Science Target–Calibrator”), while random errors are estimated from the scatter on the data. However, possibly large systematic errors remain even after calibration. The reason for this is a variation of experimental conditions between the measurements on the PSF reference(s) and the science target: either turbulence (i.e., AO correction), airmass, position on the detector, or position angle (i.e., angle-dependent aberrations in the optics).

## 5.2. Fitting

A thorough treatment of the variance-covariance matrix of T3 requires at least as many independent measurements of each target as independent T3 (Ireland 2013), which is in practice not feasible. Instead, the variance-covariance matrix can be modeled (Kraus et al. 2008), or an additional error term can be added in quadrature (i.e., sum of independent variances) to the error estimation on the data and tuned to normalize the reduced chi-square distribution of the fitting (Hinkley et al. 2011). However, the nonlinear form of the complex visibility expression to evaluate, see Equation (3), makes it impossible to determine the number of degrees of freedom; the reduced chi-square cannot be used to compare models (Andrae et al. 2010). Rather than tuning it using reduced

chi-square, we choose to fit the additional variance term $\sigma_{systematic}^2$ according to a maximum-likelihood estimation. We define the total variance $\sigma_{i,total}^2$ to be

$$ \sigma_{i,total}^2 = \sigma_{i,random}^2 + \sigma_{systematic}^2 \quad (4) $$

in order to account for error correlation and uncalibrated systematic errors of the T3.

We use an MCMC simulation using the emcee pure-Python implementation of Goodman & Weare’s Affine Invariant Ensemble sampler Foreman-Mackey et al. (2013) to maximize the full log-likelihood expression $\ln \mathcal{L}$ defined by

$$ -2 \ln \mathcal{L} = \sum_i \left[ \frac{(\text{Data}_i - \text{Model}_i)^2}{\sigma_{i,total}^2} - \ln \left( \frac{2\pi}{\sigma_{i,total}^2} \right) \right], \quad (5) $$

where the index i iterates over each (U, V) coordinate of a data set. While the first term under the sum favors high values of $\sigma_{i,total}$ (hence higher values of the fitted additional variance $\sigma_{systematic}$), the second term penalizes the increase of $\sigma_{systematic}$. The balance between these terms ensures the convergence of the estimation.

In order to compare models, we calculate the Bayesian information criteria (BIC) at the maximum likelihood of each model. The BIC is defined by Schwarz (1978), such as

$$ \text{BIC}_{\text{Model}} = 2(k \ln(n) - \ln \mathcal{L}_{\text{Model}|\text{max}}), \quad (6) $$

where k is the number of free parameters in the model and n is the sample size. The presence of k weighted with the sample size ensures that the more simple models will be preferred (to a certain extent) over more complex models. We estimate the relative probability of all models with respect to the null hypothesis of a unique central star without any disk using likelihood ratio tests (LR-Test) in its most general formalism:

$$ \begin{aligned} \text{LR}_{\text{Model/null}} &= \exp\left(\frac{\text{BIC}_{\text{null}} - \text{BIC}_{\text{Model}}}{2}\right) \\ &= \frac{1}{n^{3N}} \frac{\mathcal{L}_{\text{Model}|\text{max}}}{\mathcal{L}_{\text{null}}}, \end{aligned} \quad (7) $$

where N is the number of point sources in the model.

## 5.3. Results

We compute the relative probability for all models on each of the six data sets, which we tabulate as log of Equation (7) in Table 7. All six data sets except Lp-band epoch 3 favor the 2 point-source model rather than the 0, 1, or 3 point-source models, and all models favor the presence of a disk over a disk-less star (with the exception of Ks-band epoch 1, where the relative probabilities of a star with and without a disk are nearly identical). Since the fits on all T3 data sets were performed independently, there is a very clear statistical significance of a model that comprises a disk and two additional point sources. Note that given the relatively low flux from the point sources compared to that of the central star and the disk altogether, the additional two point sources do not significantly impact the fit of the radiative transfer disk model on the VIS2 data.

Finding the most probable model among several does not, however, ensure that the fit, hence the model, is satisfactory. We define residuals to be the normalized error of the model
<table><caption>**Table 7**
Logarithm of the Relative Probability for Each Model and All Data Sets with Respect to the Null Hypothesis of a Single Star Without a Disk</caption><thead><tr><th>Model</th><th>NULL</th><th>1PS</th><th>2PS</th><th>3PS</th></tr></thead><tbody><tr><td></td><td colspan="4">Ks epoch 1, n = 60</td></tr><tr><td>S</td><td>...</td><td>-1.16</td><td><strong>4.68</strong></td><td>-0.70</td></tr><tr><td>S + D</td><td>0.01</td><td>-1.36</td><td><strong>4.64</strong></td><td>-0.02</td></tr><tr><td></td><td colspan="4">Lp epoch 1, n = 60</td></tr><tr><td>S</td><td>...</td><td>84</td><td>178</td><td>166</td></tr><tr><td>S + D</td><td>-0.46</td><td>96</td><td><strong>223</strong></td><td>223</td></tr><tr><td></td><td colspan="4">Mp epoch 1, n = 30</td></tr><tr><td>S</td><td>...</td><td>12.4</td><td>12.2</td><td>4.6</td></tr><tr><td>S + D</td><td>-2.8</td><td>15.1</td><td><strong>15.3</strong></td><td>10.2</td></tr><tr><td></td><td colspan="4">Lp epoch 2, n = 30</td></tr><tr><td>S</td><td>...</td><td>32</td><td>72</td><td>65</td></tr><tr><td>S + D</td><td>-0.25</td><td>38</td><td><strong>108</strong></td><td>87</td></tr><tr><td></td><td colspan="4">Lp epoch 3, n = 45</td></tr><tr><td>S</td><td>...</td><td>52</td><td>113</td><td>131</td></tr><tr><td>S + D</td><td>-0.21</td><td>60</td><td>146</td><td><strong>185</strong></td></tr><tr><td></td><td colspan="4">Lp epoch 4, n = 60</td></tr><tr><td>S</td><td>...</td><td>51</td><td>176</td><td>227</td></tr><tr><td>S + D</td><td>0.06</td><td>59</td><td><strong>215</strong></td><td>210</td></tr></tbody></table>

**Note.** The higher the value the better. A negative value indicates that the model performs worse than the null hypothesis. Bold values highlight the most probable model for all data sets. The relative probability values for different data sets cannot be compared between each other. *n* is the number of independent data points (15 per acquisition sequence, see Table 2). NULL ≡ null hypothesis, i.e., single star without disk, D ≡ Disk model image in the corresponding band, PS ≡ point source model object.

such that

$$ \text{Residuals}_i = \frac{\text{Data}_i - \text{Model}_i}{\sigma_{i,\text{total}}}, \qquad (8) $$

where the index *i* iterates over each (U, V) coordinate of a data set. Table 8 shows that the additional systematic error $\sigma_{systematic}$ in Ks and Lp bands is ≈0.8, and reaches 4° in the Mp band. These uncalibrated systematic errors are consistent with previous performance and slightly worse SPHERE-SAM mode: ≈0.5 (SPHERE SAM commissioning team, private communication). This table also shows that on all data sets, the rms of residuals calculated using $\sigma_{i,\text{total}}$ is very close to unity, meaning this model does not under- or over-fit the T3 data. Note that the fact that rms values of residuals is unity is a pure consequence of the likelihood maximization of the additional variance parameter, not an a posteriori normalization. The average of residuals is only a fraction of unity, which denotes only a slight deviation from normal residuals (i.e., unbiased).

The results for the maximum-likelihood parameter estimation for the favored model with a disk and 2 point source are tabulated in Table 8. Although point sources were free to take any positions and fluxes, the parameter estimation locates both point sources at a consistent spot at all epochs and wavelengths. The first point source is ≈105 mas on the west of the star with a contrast ratio of ≈3.3 mag in the L band; the second point-source location is somewhat noisier due to its proximity to the star, though still consistently found at ≈45 mas on the northeast with a contrast ratio of ≈2.5 mag in the L band.

These locations correspond closely to the inner rim of the VSP disk.

### 5.4. Point-source Characterization

Fitting point sources to T3 data is equivalent to finding the photo-center of asymmetrically emitting regions; it does not constrain their size nor shape. Additional fits were carried using uniform disks (3 point-source parameters and an additional diameter), Gaussian blobs (3 point-source parameters and a semimajor axis) and two-dimensional Gaussian blobs (3 point-source parameters, 2 semimajor axes, and a position angle) instead of point sources. However, the data was not able to constrain well the additional parameters well: the point sources are mostly unresolved (i.e., angular extension ≲3 au).

We determine the color-temperature of both point sources at epoch 1, using their contrast ratios in Ks, Lp, and Mp bands. We assume a blackbody behavior constrained by its effective temperature and its radiating surface, assumed circular and defined by its radius. We use *emcee* to maximize the likelihood of the temperature parameter *T* taken in [100, 10,000] K and the logarithm of the radius parameter *R* taken in [0.1, 40] R⊙.

The maximum-likelihood color temperature found for the first point source is ≈1350 K, see Figure 10. However, 1σ error bars yield several hundreds of Kelvin because of large error bars on the point-source contrast ratios obtained with T3 data. Still, this temperature corresponds closely to PAH-temperatures found with radiative transfer modeling at the unprojected point-source location, ≈12.5 au; refer to Figure 6.

The second point source shows a somewhat lower temperature, most likely of ≈1000 K. Here again, 1σ error bars yield several hundreds of Kelvin large: from 600 to 1300 K, refer to Figure 10. Due to its proximity to the star, at or below 0.5λ/D, the angular separation of the second point source is particularly degenerate with its contrast ratio. However, the position angle remains mostly unaffected and shows no significant angular movement with epochs.

Discrepancy on the epoch 1 locations of the second point source between the different bands may be due to several factors, given that fitting T3 data only gives the location of the photo-center of an asymmetrically bright region. First, this location may change with wavelength due to the intrinsic cause of the asymmetry: a clump of segregated grains, an arm of material showing decreasing temperatures as it spirals out from the star, or visual effects of optical depth, opacity, and shadows. Second, the data—especially the sample size and the error bars—allowed to fit only two point sources. This second point source may, however, well be the (poorly defined) photo-center of a more complex system of asymmetries, unresolved by the data.

The first point source in the Lp band shows a clear rotational movement through different epochs. It is unlikely related to the (U, V) field rotation, since it shows much smaller variations over the four epochs, and not monotonic. This movement is compatible with a circular orbit in the VSP disk and in the same direction as the gas kinematic found by Brown et al. (2012a) and Bruderer et al. (2014; with south as the near side). If we assume a circular orbit, the origin of the over-brightness is located at a radius of ≈14.2 au, and a height of ≈3.5 au (see below). This over-brightness location corresponds to the densest region of the VSP ring in our radiative transfer modeling, as Figure 5 tells. However, the nature of the over-brightness is not clearly constrained. It may be (1) an orbiting
<table><caption>**Table 8**
Top Table: Overview of the Quality of the T3 Data Fit with the Favored Model (2 Point Sources and Disk), for all Data Sets</caption><thead><tr><th>Data Set</th><th>⟨σ<sub>i,random</sub>⟩ (°)</th><th>σ<sub>systematic</sub> (°)</th><th>⟨σ<sub>i,total</sub>⟩ (°)</th><th>rms(Residuals<sub>i</sub>)</th><th>⟨Residuals<sub>i</sub>⟩</th></tr></thead><tbody><tr><td>Ks epoch 1</td><td>1.6</td><td>0.86 ± 0.17</td><td>1.82</td><td>0.97</td><td>0.35</td></tr><tr><td>Lp epoch 1</td><td>0.27</td><td>0.68 ± 0.17</td><td>0.73</td><td>1.01</td><td>0.20</td></tr><tr><td>Mp epoch 1</td><td>1.9</td><td>4.0 ± 0.6</td><td>4.5</td><td>1.00</td><td>0.10</td></tr><tr><td>Lp epoch 2</td><td>0.44</td><td>0.63 ± 0.23</td><td>0.77</td><td>0.97</td><td>-0.21</td></tr><tr><td>Lp epoch 3</td><td>0.25</td><td>0.86 ± 0.11</td><td>0.89</td><td>1.00</td><td>0.18</td></tr><tr><td>Lp epoch 4</td><td>0.26</td><td>0.97 ± 0.11</td><td>1.0</td><td>1.00</td><td>0.26</td></tr></tbody></table><table><thead><tr><th rowspan="2">Data Set</th><th colspan="3">First Point-source</th><th colspan="3">Second Point-source</th></tr><tr><th>Sep. (mas)</th><th>PA (°)</th><th>Δmag</th><th>Sep. (mas)</th><th>PA (°)</th><th>Δmag</th></tr></thead><tbody><tr><td>Ks epoch 1</td><td>101 ± 4</td><td>281 ± 4</td><td>4.9 ± 0.3</td><td>69 ± 5</td><td>16 ± 8</td><td>4.8 ± 0.3</td></tr><tr><td>Lp epoch 1</td><td>105 ± 2.5</td><td>270.3 ± 1.1</td><td>3.35 ± 0.12</td><td>41<sup>+9</sup><sub>-3</sub></td><td>41 ± 4</td><td>2.7 ± 0.5</td></tr><tr><td>Mp epoch 1</td><td>89 ± 11</td><td>269 ± 4</td><td>2.9 ± 0.3</td><td>39 ± 9</td><td>52 ± 12</td><td>1.8<sup>+1.0</sup><sub>-0.6</sub></td></tr><tr><td>Lp epoch 2</td><td>111 ± 2.5</td><td>275.1 ± 1.3</td><td>3.20 ± 0.20</td><td>34<sup>+15</sup><sub>-2</sub></td><td>47 ± 4</td><td>1.8<sup>+1.7</sup><sub>-1.0</sub></td></tr><tr><td>Lp epoch 3</td><td>109 ± 3</td><td>276.5 ± 1.3</td><td>3.40 ± 0.05</td><td>39<sup>+10</sup><sub>-2</sub></td><td>44 ± 4</td><td>2.3<sup>+0.7</sup><sub>-0.5</sub></td></tr><tr><td>Lp epoch 4</td><td>119 ± 4</td><td>284.6 ± 1.1</td><td>3.25 ± 0.06</td><td>50 ± 4</td><td>43.3 ± 1.1</td><td>2.72 ± 0.22</td></tr></tbody></table>

**Note.** The total resulting error (⟨σᵢ,total⟩) defined by Equation (4), composed of the maximum-likelihood estimation of the additional error term (σsystematic), and the average random errors of the data (⟨σᵢ,random⟩). The rms of residuals remains close to unity, while its average remains close to zero, meaning no or little bias. Bottom table: maximum-likelihood estimation of the point-source parameters for the most probable model, consisting of two point sources in addition to a central star with the disk found in the radiative transfer VIS2 data fitting, for all data sets.

<figure><img src="image_11.png" /><figcaption>**Figure 10.** One- and two-dimensional histograms showing the variance, covariance, and maximum likelihood (blue lines) of R and T parameters obtained from the fit. First point source is on the left, the second one is on the right. Dashed lines on one-dimensional histograms represent the 16th, 50th (i.e., median), and 84th percentiles, while isocontours on two-dimensional histograms represent 0.5, 1, 1.5, and 2σ areas. The y-axes scaling on one-dimensional histogram is arbitrary; it shows the relative probability distribution function.</figcaption></figure>

over-density of VSP material, or (2) the central star orbiting a center of mass, shifted from the center of the disk due to a massive companion.

Figure 11 shows the Lp-band point-source location in the disk as well as two circular-orbit solutions. The first solution is an orbit with a semimajor axis of *a* = 14.2 au. It lies flat 3.5 au above the mid-plane (hence inclination ≈0°). This suggests that either the clump of material producing the over-brightness is mostly pressure-supported or that only the top layer of it is

detected. This latter option, however, implies that the center of mass and the photo-center of the clump are not co-located, inconsistent with the high transparency of this VSP ring (see Figure 7). The second solution is an inclined orbit, *i* = 14°, with a similar semimajor axis, *a* = 14.6 au; this implies a gravity-dominated clump of material, orbiting the star.

In both orbital solutions, the over-brightness moves at an unprojected Keplerian velocity of 12.7 ± 1.9 km s⁻¹ in average over the four epochs (similar values are obtained for other
<figure><img src="image_12.png" /><figcaption>**Figure 11.** Locations of the first point source in Lp band through the four epochs (noted 1–4). Red and yellow ellipsoids are 1σ and 2σ uncertainties for these locations. The background image represents the Lp-band model-image of the radiative transfer disk. It is shown in linear scale, the star was subtracted. Purple, red, and blue lines represent a circular orbit of semimajor axis *a* = 14.2 au for a body located, respectively, 3.5 au above the mid-plane, in the mid-plane, and 3.5 au below the mid-plane. The black line represents a 14° inclined circular orbit with the same *a*.</figcaption></figure>

intervals, between epochs 1–3 and 3–4); the total orbital period at this pace is ≈34 years. The theoretical *Keplerian velocity* at such a distance and given the stellar mass is $v = \sqrt{G \times m_*/r} = 12.5 \text{ km s}^{-1}$.

# 6. Discussion

## 6.1. Age of IRS-48

Finding a brighter star, hence younger (4 Myr), partly solves an evolutionary puzzle on this target.

First, the discrepancy between the molecular Ophiuchus cloud age of ≈1 Myr (Luhman & Rieke 1999) and the stellar age is widely decreased, even though a fair difference remains that cannot be explained to date because this is the case for all class II disks in Ophiuchus. IRS-48 might have been part of an earlier star formation episode in Ophiuchus, as postulated by Brown et al. (2012a).

Second, the disk of our model is entirely dust-depleted up to 55 au, except for a thin VSP disk between 11 and 26 au. This implies that the evolutionary stage of the disk is probably more advanced than what was previously thought. Accretion was estimated at $\approx 10^{-8.5} M_\odot \text{ yr}^{-1}$ by Salyk et al. (2013) and Follette et al. (2015) using $Pa\beta$ and $Br\gamma$ lines. This non-negligible accretion rate raises the question of the origin of such material and does not seem consistent with either the old model (where one finds $8.0 \times 10^{-12} M_\odot$ dust in the inner-most disk, and $8.0 \times 10^{-10} M_\odot$ PAH until 50 au), or our new model ($3.7 \times 10^{-10} M_\odot$ VSP at 11–26 au). Moreover, the disk was found to be gas-depleted until 20 au, which puts a very hard constraint on how much material is to be found in the first tens of astronomical units, available for accretion. This means that some active replenishment processes take place to sustain accretion, or that IRS-48 is on the verge of becoming a debris disk.

Finally, a younger age for IRS-48 together with a later disk-evolution stage is consistent with the fact that massive stars erode their disk at a larger pace than solar-type stars. Indeed, earlier studies on Upper Scorpius OB association (Carpenter

et al. 2006; Williams & Cieza 2011) showed that ≈20% of 127 K- and M-stars are surrounded by optically thick disks, while none of the 30 F- and G-stars showed such disks for a wavelength ≈16 µm, implying that primordial circumstellar disks are consumed faster by a massive star. This result is confirmed by Ribas et al. (2015) who show that only 2% of $M \ge 2 M_\odot$ stars aged ≥3 Myr show a protoplanetary disk.

## 6.2. Disk Morphology, Stellar Binarity, and Planetary Formation

### 6.2.1. 30 au CO Gas Ring

Brown et al. (2012a), Bruderer et al. (2014), and van der Marel et al. (2016) find a ≈30 au CO ring with a small radial extent, compatible with a Keplerian velocity for a $2.0 M_\odot$ star with an inclination of 50°. In addition, they find a gas-depleted region inside ≈20–25 au.

Although our data is blind to gas—consequently our model does not incorporate any gas structure—the new IRS-48 model and the later disk-evolution stage postulated remain qualitatively compatible with these previous findings. Indeed, a photo-evaporation-dominated disk would show a clear depletion of gas and dust in the inner tens of astronomical units. The VSP ring from 11 to 26 au acts as a very efficient shield of UV-radiation, effectively protecting gas further out from aggressive photo-evaporation. A thin ring of gas at 30 au—and potentially as close as 26 au—appears to be a consistent consequence of such a configuration.

### 6.2.2. VSP Ring Inner Truncation, Binarity

Ghez et al. (1993), Leinert et al. (1993), and Reipurth & Zinnecker (1993) detect an over-abundance by a factor of two of binaries among the young stars in the Taurus region when compared to the results for the main sequence (MS) stars. Although further studies with larger populations slightly lower this over-abundance, it appears clear that YSOs have a binarity over-abundance compared to MS stars, to a degree that is, however, still under discussion as underlined by Ratzka et al. (2005).

Our data resolves a somewhat sharp truncation of the inner rim at 11 au. Besides photo-evaporation and accretion, a stellar companion with a semimajor axis of the order of a few astronomical units would reproduce the VSP-ring truncation at 11 au (Artymowicz & Lubow 1994). To the knowledge of the authors, no observation already performed in the visible or infrared, had the sensitivity to resolve such a binary system. According to the PMS evolutionary models of Siess et al. (2000), a K0-type star (5400 K) of 4 Myr has a luminosity of about $7.5 L_\odot$ and a mass of ≈2.2 $M_\odot$. This luminosity remains small compared to that of IRS-48 (45 $L_\odot$): any stellar type cooler than K0 would remain a nearly invisible companion for the visible and infrared instruments, especially at ≈25 mas typical spatial separations.

An indirect way to detect late-type companions is to check if such A-type stars emit in X-ray. Stelzer et al. (2006a, 2006b) show that, given that A-type stars are not expected to have a corona and emit X-ray, a high-energy emission from these stars is generally associated with late-type binarity. In the case of IRS-48, neither ROSAT and XMM archive data nor Chandra (Imanishi et al. 2001) showed X-ray emission at this location. Although soft X-ray would be absorbed by interstellar absorption, hard and extreme ones should be detected. Indeed,
<figure><img src="image_13.png" /><figcaption>**Figure 12.** Schematic of the IRS-48 disk as presented in the new model. The outer disk is from Geers et al. (2007a), the CO emission is from Brown et al. (2012a), the millimeter-asymmetry from van der Marel et al. (2013), the inner 20 au gas depletion from Bruderer et al. (2014). The proposed gas dip from 35 to 50 au is proposed from Bruderer et al. (2014) findings on a first gas depletion between 20 and 55 au.</figcaption></figure>

several very young stars (even Class I) are detected in X-ray in the core region of ρ Oph by Ozawa et al. (2005). Given the completeness of several X-ray studies as showed in their Figure 6, and the $A_v = 12.9$ mag for IRS-48, a 4 Myr old companion with $T_{eff} \gtrsim 3000$ K (M6 star and brighter) should already have been detected around IRS-48.

Also, a heavy companion, $\gtrsim 1 M_\odot$ (corresponding to a star brighter than $\approx$K6 stellar type on the 4 Myr isochrone) would reveal its presence through its mass, and the Keplerian velocity field it imposes on its circumbinary disk. Brown et al. (2012a) find a $2 M_\odot$ to explain a Keplerian rotation of a gas ring at 30 au; this mass represents a lower limit given that the gas ring must have a sub-Keplerian velocity. IRS-48 does not likely have a bright and massive secondary star, and one can set a total binary mass at $\approx 4 M_\odot$ and below.

IRS-48 is most probably not a binary star. However, a low-mass companion with $M \lesssim 0.5 M_\odot$ and $L \lesssim 0.01 L_\odot$ would be qualitatively compatible with X-ray non-detection, CO 30 au rotation speed and with our model.

To date, the inner-truncatin at 11 au is still to be explained.

### 6.2.3. 26–55 au Cavity

Both the presence of the outer-rim of the VSP ring and the inner-rim of the outer disk require the existence of an explanatory phenomenon. Based on the millimeter asymmetry observed in the southern part of the disk, a $2 M_\odot$ star and Hill radius estimates, van der Marel et al. (2013) suggested the presence of a planet at 17–20 au with $M \gtrsim 10 M_{Jup}$.

If we carry similar first-order Hill radius estimates with the new model disk structure and a higher-mass star, we speculate a $M_p \approx 3.5 M_{Jup}$ planet at a radius of $a_p \approx 40$ au. Indeed, Dodson-Robinson & Salyk (2011) find that the outer edge of the gas cavity created by a planet is expected at $\approx 5$ Hill radii $r_H$, defined by $3 \times (r_H/a_p)^3 = M_p/M_*$. Pinilla et al. (2012) find that for planets with masses between 1 and $3 M_{Jup}$, dust dynamically accumulates at $7 r_H$ from their location. With the later planetary parameters, we are able to obtain the dust bump at 63 au (van der Marel et al. 2013), together with a VSP-ring

bump at 17 au, and gas densities peaking at 26 and 56 au (all radii given as distances from the star). These radii correspond closely to the structure of our model; a schematic picture is shown on Figure 12.

### 6.3. PAH and VSG Evolution

Dullemond et al. (2007) suggested that PAH and VSG do not settle well compared to CTG, and found that, under conditions of low turbulence, strong PAH emission features can “hide” a substantial amount of settled thermal silicate grains.

We were indeed able to add up to five times the mass of the VSP ring ($5 M_{VSP}$) worth of settled CTG in the VSP ring, without affecting the VIS2 data nor the SED fits. However, adding such grains in the close vicinity of the star affects the MIR emission in such a way that the 18.7 µm model image does not display the central dip with two bright spots at $\pm 55$ au from the star like the data image does (Figure 8). Instead, a contiguously bright region lies along the whole semimajor axis of the disk. This result shows that the VSP ring most probably does not contain more thermal grains than $5 \times M_{VSP} \approx 2 \times 10^{-9} M_\odot$.

The presence of PAH in many circumstellar disks reveals that in these objects the dust coagulation process was apparently not effective enough to remove the PAH from these disks, or that some other process continuously replenishes PAH grains (Dullemond et al. 2007). Out of theoretical and laboratory works, Jochims et al. (1994) showed that the largest and most regular PAH species are extremely stable to destructive radiations, i.e., photo-evaporation. Indeed, large PAHs of 20–30 carbon atoms will preferably relax through PAH emission in bands rather than photo-fragmentation. As a comparison, one of the most stable PAH species, circumcoronene C54H18, has a typical size of $4.87$ Å. PAH used in our model span $4–10$ Å (20–100 carbon atoms); VSG span $10–20$ Å (100–1000 carbon atoms), see Tielens (2008).

Geers et al. (2006) find a PAH-to-dust mass fraction of 6% based on the abundance of $5 \times 10^{-5}$ carbon atoms locked in PAH molecules per H nuclei ($C_{VSP}/H$). If we take updated values for PAH and VSG abundances from Tielens (2008; their
Table 2), we find $C_{VSP}/H = 3.5\%$. Our model showing a VSP-to-dust mass fraction >20% in the VSP ring between 11 and 26 au clearly highlights a depletion of a factor of ≈5–6 of classical dust grains up to 0.3 mm compared to VSP.

We are left with the problem of where one shall study the disappearance rate of VSP through grain growth, radiation pressure, or inward drift, versus the rate of some replenishment process, e.g., the collisional destruction of larger grains or inward drift from an outer reservoir.

There is presently little experimental data on the size distribution resulting from sub-micrometer grain–grain collisions. Tielens (2008) report that theoretical studies on grain–grain collisions result in a mass distribution with a power index of ≈−3.3. Although the minimum size is not constrained in this later study, given the layered structure of graphitic materials, it is expected that the smallest fragments are likely small, two-dimensional structures (e.g., PAH-like molecular species). Further processing by atomic reactions, as well as UV processing, may then quickly transform these species into compact PAH (Tielens 2008). In the frame of interstellar medium shocks, Jones et al. (1996) calculate that the catastrophic destruction (i.e., complete destruction of the grains, rather than cratering) of a 1000 Å grain by a 50 Å grain needs a critical velocity of approximately 75 km s⁻¹; larger grains need larger collisional energy. They also conclude that interstellar grain-growth processes shall be much more efficient than catastrophic destruction in order to explain the observation that most of the interstellar dust comes in grains larger than ≈1000 Å.

The orbiting velocity of the first point source in the VSP ring around IRS-48 is ≈12.7 km s⁻¹, similar to the Keplerian velocity at that distance for the 2.5 $M_\odot$ star. Moreover, A0 stars do not generally show strong winds or emit high-energy radiation. It seems very unlikely that IRS-48 possesses a high catastrophic collision rate within its VSP ring at 11–26 au. It is therefore expected that PAH and VSG cannot be replenished easily through grain–grain collisions in the VSP ring, in spite of a possible co-located reservoir of CTG up to five times the mass of VSP.

The radiation pressure estimate of the star can be achieved by calculating the stellar constant $W \approx 45,000 \text{ W m}^{-2}$ at 1 au. The radial acceleration from radiation pressure produced on a particle at radius $r$ is obtained with $a_p = W S_{\text{eff}} c^{-1} r^{-2}$, where $S_{\text{eff}}$ is the effective cross-section of the particle and $c$ is the speed of light in vacuum. It is found to be 39 (resp. 42, 17) times larger than the gravitational potential ($G M_* r^{-2}$) on the neutral PAH particles (respectively, ionized PAH and VSG). This probably shows that the inward drift of such small particles is very low. The previous IRS-48 stellar model with radius $1.8 R_\odot$ shows a somewhat lower but still overwhelmingly strong radiation pressure compared to gravitational potential for these species: 16, 17, and 7 times larger on average for neutral, ionized PAH and VSG.

The fact that the VSP ring was observed to be consistently similar at a radius between ≈11 and 26 au through the four epochs of observation spanning two years time, highlights that the VSP ring is governed by a much more complex balance of physical processes. In order to remain observationally in place, it must either be (1) constantly replenished by the outer-reservoir outside 55 au—possibly channeling triggered by the 40 au companion—or (2) kept in place as-is through (2a) gravitational shepherding from additional companions, or (2b)

coating of large grains—unaffected by radiation pressure—by PAH molecules.

We wish to thank Ewine van Dishoeck for her insight on SMA data, Vincent Geers for providing the VISIR NIR images, and Daniel Rouan for our talks on PAH and VSG. We also thank the Paranal staff for their support during the observations. This work was supported by the French National Agency for Research (ANR-13-JS05-0005) and the European Research Council (ERC-STG-639248). Based on observations collected at the European Southern Observatory (ESO) during runs 086.C-0497(A), 087.C-0450(B), 088.C-0527(A), and 089.C-0721(A).

## References

Alexander, R. D., Clarke, C. J., & Pringle, J. E. 2006, MNRAS, 369, 216

Andrae, R., Schulze-Hartung, T., & Melchior, P. 2010, arXiv:1012.3754

Andrews, S. M., Rosenfeld, K. A., Kraus, A. L., & Wilner, D. J. 2013, ApJ, 771, 129

Andrews, S. M., Rosenfeld, K. A., Wilner, D. J., & Bremer, M. 2011, ApJL, 742, L5

Andrews, S. M., & Williams, J. P. 2007, ApJ, 659, 705

Artymowicz, P., & Lubow, S. H. 1994, ApJ, 421, 651

Brown, J. M., Herczeg, G. J., Pontoppidan, K. M., & van Dishoeck, E. F. 2012a, ApJ, 744, 116

Brown, J. M., Rosenfeld, K. A., Andrews, S. M., Wilner, D. J., & van Dishoeck, E. F. 2012b, ApJL, 758, L30

Bruderer, S., van der Marel, N., van Dishoeck, E. F., & van Kempen, T. A. 2014, A&A, 562, A26

Bryden, G., Chen, X., Lin, D. N. C., Nelson, R. P., & Papaloizou, J. C. B. 1999, ApJ, 514, 344

Calvet, N., D'Alessio, P., Hartmann, L., et al. 2002, ApJ, 568, 1008

Camps, P., Misselt, K., Bianchi, S., et al. 2015, A&A, 580, A87

Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245

Carmona, A., Pinte, C., Thi, W. F., et al. 2014, A&A, 567, A51

Carpenter, J. M., Mamajek, E. E., Hillenbrand, L. A., & Meyer, M. R. 2006, ApJL, 651, L49

Cieza, L., Padgett, D. L., Stapelfeldt, K. R., et al. 2007, ApJ, 667, 308

Clarke, C. J., Gendrin, A., & Sotomayor, M. 2001, MNRAS, 328, 485

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, yCat, 2246, 0

Dodson-Robinson, S. E., & Salyk, J. 2011, ApJ, 738, 131

Draine, B. T. 2006, ApJ, 636, 1114

Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89

Draine, B. T., & Li, A. 2007, ApJ, 657, 810

Dullemond, C. P., & Dominik, C. 2005, A&A, 434, 971

Dullemond, C. P., Henning, T., Visser, R., et al. 2007, A&A, 473, 457

Erickson, K. L., Wilking, B. A., Meyer, M. R., Robinson, J. G., & Stephenson, L. N. 2011, AJ, 142, 140

Fedele, D., Bruderer, S., van Dishoeck, E. F., et al. 2013, A&A, 559, A77

Follette, K. B., Grady, C. A., Swearingen, J. R., et al. 2015, ApJ, 798, 132

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

Geers, V. C., Augereau, J.-C., Pontoppidan, K. M., et al. 2006, A&A, 459, 545

Geers, V. C., Pontoppidan, K. M., van Dishoeck, E. F., et al. 2007a, A&A, 469, L35

Geers, V. C., van Dishoeck, E. F., Visser, R., et al. 2007b, A&A, 476, 279

Ghez, A. M., Neugebauer, G., & Matthews, K. 1993, AJ, 106, 2005

Gorti, U., & Hollenbach, D. 2008, ApJ, 683, 287

Heng, K., & Kenyon, S. J. 2010, MNRAS, 408, 1476

Hinkley, S., Carpenter, J. M., Ireland, M. J., & Kraus, A. L. 2011, ApJL, 730, L21

Imanishi, K., Koyama, K., & Tsuboi, Y. 2001, ApJ, 557, 747

Indebetouw, R., Mathis, J. S., Babler, B. L., et al. 2005, ApJ, 619, 931

Ireland, M. J. 2013, MNRAS, 433, 1718

Ireland, M. J., & Kraus, A. L. 2008, ApJL, 678, L59

Jennison, R. C. 1958, MNRAS, 118, 276

Jochims, H. W., Ruhl, E., Baumgartel, H., Tobita, S., & Leach, S. 1994, ApJ, 420, 307

Jones, A. P., Tielens, A. G. G. M., & Hollenbach, D. J. 1996, ApJ, 469, 740

Kim, S.-H., Martin, P. G., & Hendry, P. D. 1994, ApJ, 422, 164

Kraus, A. L., Ireland, M. J., Martinache, F., & Lloyd, J. P. 2008, ApJ, 679, 762

Lacour, S., Tuthill, P., Amico, P., et al. 2011, A&A, 532, A72

Leinert, C., Zinnecker, H., Weitzel, N., et al. 1993, A&A, 278, 129

Li, A., & Draine, B. T. 2001, ApJ, 554, 778
Li, A., & Greenberg, J. M. 1997, A&A, 323, 566

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

Lombardi, M., Lada, C. J., & Alves, J. 2008, A&A, 489, 143

Lubow, S. H., Seibert, M., & Artymowicz, P. 1999, ApJ, 526, 1001

Luhman, K. L., & Rieke, G. H. 1999, ApJ, 525, 440

Maaskant, K. M., Min, M., Waters, L. B. F. M., & Tielens, A. G. G. M. 2014, A&A, 563, A78

McClure, M. K., Furlan, E., Manoj, P., et al. 2010, ApJS, 188, 75

Monnier, J. D. 2000, in Principles of Long Baseline Stellar Interferometry, ed. P. R. Lawson (Pasadena, CA: JPL), 203

Natta, A., Prusti, T., & Krugel, E. 1993, A&A, 275, 527

O'Donnell, J. E. 1994, ApJ, 422, 158

Owen, J. E., Ercolano, B., & Clarke, C. J. 2011, MNRAS, 412, 13

Ozawa, H., Grosso, N., & Montmerle, T. 2005, A&A, 429, 963

Pinilla, P., Benisty, M., & Birnstiel, T. 2012, A&A, 545, A81

Pinte, C., Ménard, F., Duchêne, G., & Bastien, P. 2006, A&A, 459, 797

Ratzka, T., Köhler, R., & Leinert, C. 2005, A&A, 437, 611

Reipurth, B., & Zinnecker, H. 1993, A&A, 278, 81

Ribas, Á., Bouy, H., & Merín, B. 2015, A&A, 576, A52

Salyk, C., Herczeg, G. J., Brown, J. M., et al. 2013, ApJ, 769, 21

Schwarz, G. 1978, AnSta, 6, 461

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337

Siebenmorgen, R., & Krügel, E. 2010, A&A, 511, A6

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

Skrutskie, M. F., Dutkevitch, D., Strom, S. E., et al. 1990, AJ, 99, 1187

Stelzer, B., Huélamo, N., Micela, G., & Hubrig, S. 2006a, A&A, 452, 1001

Stelzer, B., Micela, G., Hamaguchi, K., & Schmitt, J. H. M. M. 2006b, A&A, 457, 223

Strom, K. M., Strom, S. E., Edwards, S., Cabrit, S., & Skrutskie, M. F. 1989, AJ, 97, 1451

Tielens, A. G. G. M. 2008, ARA&A, 46, 289

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

van der Marel, N., van Dishoeck, E. F., Bruderer, S., et al. 2016, A&A, 585, A58

van Kempen, T. A., van Dishoeck, E. F., Salter, D. M., et al. 2009, A&A, 498, 167

Weidenschilling, S. J. 1977, MNRAS, 180, 57

Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296

Wilking, B. A., Lada, C. J., & Young, E. T. 1989, ApJ, 340, 823

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

Wood, K., Wolff, M. J., Bjorkman, J. E., & Whitney, B. 2002, ApJ, 564, 887

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

Zacharias, N., Monet, D. G., Levine, S. E., et al. 2005, yCat, 1297, 0