A resolved and asymmetric ring of PAHs within the young circumstellar disk of IRS 48

TL;DR

This study uses high-resolution interferometry and radiative transfer modeling to reveal a resolved, asymmetric PAH ring within IRS 48's young, dust-depleted disk, suggesting a planetary influence and revising the system's age.

Contribution

It provides a new radiative transfer model incorporating PAHs, VSGs, and silicates, and revises the age and structure of IRS 48's disk with multi-epoch interferometric data.

Findings

Discovery of a resolved PAH and VSG ring at 11-26 AU.

Revised stellar age to 4 Myr, younger than previous estimates.

Identification of a potential planet at 40 AU influencing disk structure.

Abstract

For one decade, the spectral-type and age of the $\rho$ Oph object IRS-48 were subject to debates and mysteries. Modelling 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 its first 15AU, 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 SED 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 (42L$_{\odot}$) central-star with A$_v$=12.5mag and R$_v$=6.5 requires less near-infrared excess: the inner-most disk at $\approx$1AU is incompatible with the data. The revised stellar parameters place this system on a 4 Myr evolutionary track, 4 times younger than previous estimations, in better agreement with the surrounding $\rho$ Oph region and disk-lifetimes observations. The disk-structure converges to a classical-grains outer-disk from 55AU combined with a fully resolved VSG\&PAH-ring, at 11-26 AU. We find two over-luminosities in the PAH-ring at color-temperatures consistent with the radiative transfer simulations; one follows a Keplerian circular orbit at 14AU. We show a depletion of a factor $\approx$5 of classical dust grains compared to VSG\&PAH: the IRS-48 disk is nearly void of dust-grains in the first 55 AU. A 3.5M$_{Jup}$ planet on a 40AU orbit qualitatively explains the new disk-structure.

Tables

Table 1. Table 1 Main literature data of ρ 𝜌 \rho Oph IRS-48

Parameter	Value	Ref.
Stellar position	16h27m37${.}^s$18	B14
(J2000)	-24∘ 30′35.3″
Distance	120 pc	L08
Inclination	i=50∘	G07, B14
Systemic velocity	4.55 km.s-1	vdM13
Position angle	PA=96∘	G07, B14
Stellar type	A0${}_{-1}{}^{}{}_{}{}^{+4}$	B12
Stellar Temp.	9000$\pm$550 K	B12

Table 2. Table 2 Observations made with NaCo

Date	Band	Calibrators	Seq.
14 Mar. 2011	Ks(1,2)	Elia 2-35 & 2-37	4
14 Mar. 2011	Lp(1,2,3)	Elia 2-37(1,2)	4
		Hip 86311(3)
14 Mar. 2011	Mp(1,2)	Elia 2-15	2
1 Sept. 2011	Lp(1)	Elia 2-37	2
6 Mar. 2012	Lp(1)	Elia 2-37	3
25 Mar. 2013	Lp(1)	Elia 2-37 & 2-11	4

Table 3. Table 3 NaCo filters characteristics, interferometric 0.5 λ / D 𝜆 𝐷 \lambda/D criteria and field of view with respect to filters at central wavelengths, given a target distance of 120 parsec.

	Ks	Lp	Mp
Central wavelength ($\mu$m)	2.18	3.8	4.78
Bandwidth ($\mu$m)	0.35	0.62	0.59
Smallest Baseline, 1.78m
$0.5\lambda /D$ (AU)	15.2	26.6	33.5
Field of view (AU)	190	326	544
Longest Baseline, 6.43m
$0.5\lambda /D$ (AU)	4.2	7.4	9.2
Field of view (AU)	52.6	90.8	150

Table 4. Table 4 Previously published photometry and spectra of ρ 𝜌 \rho Oph IRS-48.

Wavelength(s)	Instrument	Reference
($\mu$m)
0.43, 0.64	NOMAD	Z05
0.65, 0.8	Hydra	E11
1.2, 1.6, 2.2	2MASS	C03
3.4, 4.6	WISE	W10
3.6, 4.5	Spitzer IRAC	vK09
5.9-36.9	Spitzer IRS	MC10
60-181	Herschel PACS	F13
450	ALMA	vdM13
850	SCUBA	A07
880, 1300	SMA	B12b

Table 5. 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 ) .

Parameter	old model	New model
Stellar parameters
Temperature	9250 K	9250 K*
Luminosity	14.3 L⊙	42 L⊙
Distance	120$pc$	120$pc$*
Rv	5.5	6.5
Av	11.5	12.9
Global disks parameters
Inclination	50∘	50∘ *
PA	96∘	96∘ *
Scale height	10 AU	14 AU
Ref. radius	60 AU	60 AU*
Inner-most disk
$R_{in}$ Inner	0.4 AU	–
$R_{out}$ Inner	1 AU	–
Flaring exp. $\beta$	1.3	–
Surface density exp $p$	-1	–
Grain-size power-law	-3.5	–
Dust Mass	8.0$\times {10}^{-12}$ M⊙	–
Silicate grains Mass	1%	–
Carb. grains Mass	99%	–
Grain sizes $a$	0.03 - 30 $\mu$m	–
VSP-ring
$R_{in}$ Inner	1 AU	11 AU
$R_{out}$ Inner	50 AU	26 AU
Flaring exp. $\beta$	1.3	0.6
Surface density exp $p$	-1	-0.1
Grain-size power-law	–	-4*
Dust Mass	8.0$\times {10}^{-10}$ M⊙	3.7$\times {10}^{-10}$ M⊙
VSG Mass	0%	20%
VSG sizes	–	10 - 30 Å *
Neutral PAH Mass	50%	60$\%$
Ionized PAH Mass	50%	20$\%$
PAH sizes	5 Å	4 - 10 Å *
Outer disk
$R_{in}$	55 AU	55 AU*
$R_{out}$	160 AU	250 AU
Flaring exp. $\beta$	1.3	0.67*
Surface density exp $p$	-1	-1*
Grain-size power-law	-3.5	-3.5*
Dust Mass	1.0$\times {10}^{-5}$ M⊙	9.0$\times {10}^{-6}$ M⊙
Silicate grains Mass	70$\%$	70$\%$*
Carb. grains Mass	30$\%$	30$\%$*
Grain sizes $a$	0.03 - 4000 $\mu$m	0.03 - 4000 $\mu$m*

Table 6. Table 6 New stellar model for ρ 𝜌 \rho Oph IRS-48 central star.

	This work	F15	B12
T (K)	9250	9500	9000
R (R⊙)	2.5(1)	1.8	1.4
L (L⊙)	42	23.6	14.3
Age (Myr)	4.1(1)	8	$>$15
M (M⊙)	2.5(1)	2.2	2.0
Av	12.9	12.0	11.5
Rv	6.5	(4.0)	5.5
Ak	(1.8)	1.5	(1.6)

Table 7. Table 7 Logarithm of the relative probability for each model and all datasets with respect to the null hypothesis of a single star without disk. The highest value the better. A negative value tells that the model performs worse than the null hypothesis. Bold values highlight the most probable model for all datasets. The relative probability values for different datasets cannot be compared between each-other. n is the number of independent data point (15 per acquisition sequence, see Table  2 ).

Model	NULL	1PS	2PS	3PS
Ks epoch 1, n=60
S	–	-1.16	4.68	-0.70
S+D	0.01	-1.36	4.64	-0.02
Lp epoch 1, n=60
S	–	84	178	166
S+D	-0.46	96	223	223
Mp epoch 1, n=30
S	–	12.4	12.2	4.6
S+D	-2.8	15.1	15.3	10.2
Lp epoch 2, n=30
S	–	32	72	65
S+D	-0.25	38	108	87
Lp epoch 3, n=45
S	–	52	113	131
S+D	-0.21	60	146	185
Lp epoch 4, n=60
S	–	51	176	227
S+D	0.06	59	215	210

Table 8. Table 8 Top Table : Overview of the quality of the T3 data fit with the favored model (2 point-sources and disk), for all datasets. The total resulting error ( ⟨ σ i , t ​ o ​ t ​ a ​ l ⟩ delimited-⟨⟩ subscript 𝜎 𝑖 𝑡 𝑜 𝑡 𝑎 𝑙 \langle\sigma_{i,total}\rangle ) defined by Eq. ( 4 ), composed of the maximum likelihood estimation of the additional error term ( σ s ​ y ​ s ​ t ​ e ​ m ​ a ​ t ​ i ​ c subscript 𝜎 𝑠 𝑦 𝑠 𝑡 𝑒 𝑚 𝑎 𝑡 𝑖 𝑐 \sigma_{systematic} ), and the average random errors of the data ( ⟨ σ i , r ​ a ​ n ​ d ​ o ​ m ⟩ delimited-⟨⟩ subscript 𝜎 𝑖 𝑟 𝑎 𝑛 𝑑 𝑜 𝑚 \langle\sigma_{i,random}\rangle ). The RMS of residuals remains close to unity while its average close to zero, meaning no or little bias. Bottom Table : Maximum likelihood estimation of the point-sources 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 datasets.

Data set	$⟨{\sigma }_{i,random}⟩$ (∘)	${\sigma }_{systematic}$ (∘)	$⟨{\sigma }_{i,total}⟩$ (∘)	$RMS\left(residual{s}_i\right)$	$⟨residual{s}_i⟩$
Ks epoch 1	$1.6$	$0.86\pm 0.17$	$1.82$	$0.97$	$0.35$
Lp epoch 1	$0.27$	$0.68\pm 0.17$	$0.73$	$1.01$	$0.20$
Mp epoch 1	$1.9$	$4.0\pm 0.6$	$4.5$	$1.00$	$0.10$
Lp epoch 2	$0.44$	$0.63\pm 0.23$	$0.77$	$0.97$	$-0.21$
Lp epoch 3	$0.25$	$0.86\pm 0.11$	$0.89$	$1.00$	$0.18$
Lp epoch 4	$0.26$	$0.97\pm 0.11$	$1.0$	$1.00$	$0.26$

Equations

$$\sum(r,z)=\sum_{0}~{}r^{p}\exp\Bigg{(}-\frac{1}{2}\bigg{(}\frac{z}{h}\bigg{)}^{2}\Bigg{)},$$

$$F_{\nu }=\Omega \times B_{\nu }(T_{dust})\times (1-{\exp }^{-{\tau }_{\nu }}).$$

$$\Big{(}1+\sum_{i=1}^{N}f_{i}\Big{)}~{}\mathcal{V}=1+\sum_{i=1}^{N}f_{i}~{}\exp\Big{(}-\textbf{j}\frac{2\pi}{\lambda}~{}(U~{}\delta_{i}~{}+~{}V~{}\alpha_{i})\Big{)},$$

$${\sigma }_{i,total}^2={\sigma }_{i,random}^2+{\sigma }_{systematic}^2$$

$$-2~{}ln\mathcal{L}=\sum_{i}\Bigg{[}\frac{\big{(}Data_{i}-Model_{i}\big{)}^{2}}{\sigma_{i,total}^{2}}-\ln\Big{(}\frac{2\pi}{\sigma_{i,total}^{2}}\Big{)}\Bigg{]},$$

$$BI{C}_{Model}=2(k\ln (n)-ln{\mathcal{L}}_{Model}{\mid }_{max}),$$

$$L{R}_{Model/null}$$

$$Residual{s}_i=\frac{Dat{a}_i-Mode{l}_i}{{\sigma }_{i,total}},$$

Full text

1. Introduction

The dust- and gas-rich disks surrounding numerous pre-main-sequence (PMS) stars are of key interest for unveiling how planetary system 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 magnitudes 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 infra-red 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 where the delicate balance of such phenomena is currently questioned. Also referred as $2MASS~{}J16273718-2430350$, $GY~{}304$ or $WLY~{}2-48$, it is an highly extinct A0${}_{-1}^{+4}$ star, part of the $\rho$ Oph cloud L1688; Table 1 sums up its main data.

First listed in Wilking et al. (1989), it shows a substantial far-infrared and millimeter excess, pointing towards a classification as 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 mysteries. While 8.6 through 11.3 $\mu$m images showed high and unresolved PAH fluxes, the 18.7 $\mu$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 $\mu$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*⊙* combined with an extinction of Av=11.5, positioning the object on a – suspiciously old – 15 Myr evolutionary track.

Subsequent ALMA and VLA observation at 0.44, 1.3, and 8.8 mm carried by van der Marel et al. (2013, 2016) unveiled a millimeter-grains asymmetry at a radius of 63 AU, in the southern extents of the disk. They postulate a planet of 10 MJup located at $\approx$18 AU triggering a vortex-shaped dust-trap. Further work on PAH 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 PAH, that 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 12CO and C17O lines measured by ALMA with a resolution of 30 AU. They proposed a slight gas depletion at $\approx$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 innermost-disk at $\approx$1 AU using SED modelling, although their fit suffers from over-luminous silicate features in the 9-18 $\mu$m window.

A number of mechanisms have been proposed to cause gaps in proto-planetary disks, including extensive grain growth (Dullemond & Dominik 2005), photo-evaporation (Clarke et al. 2001), binarity (Ireland & Kraus 2008), and 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-coagulatation theories (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 as large molecules: PAH, and in a lesser extent, Very Small Grains (VSG). It has been argued that PAH either may take 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 PAH in these processes is not conclusive and it is presently unclear how important these routes are. The question remains whether these very small particles take part in grain-coating, in the main grain-growth process (and their presence in many disks are to be explained by some replenishment processes), or if their intrinsic properties prevents them to grow (and their absence in many disks are to be explained by some other assimilation processes). Gorti & Hollenbach (2008) find that the abundance of PAH and small grains in disks is critical to the temperature profile. Indeed, they show that if some PAH are present, their contribution to disk heating via grain photoelectric emission can be larger by a factor of $\approx$2 than Xray heating at r$\approx$10 AU.

This paper reports the first direct detection of the full extents 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-environment observed from near-infrared to centimeter wavelengths, and sets limits on how much silicates grains – hence replenishment – is to be expected in the PAH and VSG ring.

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 Nasmyth Adaptive Optics System (NAOS) & 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-years period; these sets are detailed in Table 2. At that time, NaCo was integrated at the Nasmyth platform of the UT4 8.2m-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 ($>$16mag in the visible), we use the NaCo built-in infrared wavefront sensor for the Adaptive Optics (AO) correction.

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 $\mu$m, with $R_{average}\approx 800$. The Lp-band image is displayed in Figure 1 both with and without PSF-subtraction. One can easily see the East-West elongated nature of the circumstellar emission, at $\approx$100 milli-arcsecond (mas), as well as an asymmetry in the disk semi-major axis; the Western extension is $\approx$1.15 times brighter than the Eastern one. However, the $\approx$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 Å $\times$2.7 mas. The slit height is 86 mas; given the short scale extents of the disk and the high inclination ($\approx$50*∘*), we expect a large contamination of the disk at the stellar position (assumed to be the located at the peak intensity in the spectra). Similarly, given the full-width-half-maximum of the PSF at $\approx$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$\lambda$/D$\approx$115 mas. The PSF-reference image shows a very low stellar-flux in the null, $\approx$3% of the peak-flux; this highlights the high strehl ratio of the observation.

The spectra in Figure 1 shows 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$\lambda$/D for the null. One can see very strong PAH emission features at 3.3 $\mu$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 West-East 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 AU given the distance of 120 parsec, as shown in Table 3. Except for the longest baselines of ALMA ($\gtrsim$1 km), such resolution cannot be reached with any other instrument nor any other NaCo mode (coronographic, full-pupil, etc). Moreover, new extreme adaptive-optics (AO) systems such as SPHERE usually have their wavefront sensors in the visible and cannot close the AO loop on this extremely extinct target ($>$16mag 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 point-spread function (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}$$\approx50 mas at 4 $\mu$m). However, the interferometric field of view remains very limited, in practice of order 2.5\lambda/B_{min}$$\approx250 mas at 4 $\mu$m (although a somewhat lower sensitivity exist 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 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 7-hole non-redundant mask with baselines ranging from 1.78 to 6.43 meters, 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-band.

2.4. Spectral Energy Distribution

In addition to the multi-wavelength interferometric data listed above, we compiled SED fluxes for IRS-48, from the literature. All data are listed in Table 4 with the corresponding references. As can be seen, the photometry from literature has a broad coverage from visible-blue to mm wavelengths. No measurements in the ultraviolet exist because of the very high extinction of the source; measurements in cm-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 (1994) work in the visible. These extinction curves are parametrized with Rv (reddening slope) and Av (absorption in V band) such that A${}_{v}\equiv$ Rv**Extinction(Bmag*-Vmag). The initial values are taken from Brown et al. (2012a), Av=11.5 and Rv=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 Lp- or Mp-band (see Figures 3). 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 $\approx$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 to $\approx$2.7 for baselines $\gtrsim$3 meter.

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

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

3.2. New model: Fitting Strategy

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. Genetic algorithm allows for a very fast convergence for cases with a large number of dimensions where basic 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 modelling over these small grid-subsets allowed us to interpolates 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 a 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:

[TABLE]

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}$ with $\beta$ the flaring exponent ($\beta>0$, of order $\approx$1) and $h_{0}$ 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 $\approx$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 well the outer disk which mostly lies outside its field of view. Hence, we systematically adopt the simplest axi-symmetrical fitting solution consistent with the literature for this outer-disk, while we focus our fitting effort on the inner parts.

Similarly to the old model, we adopt a three disk-components model: two dusty inner- and outer-disks and an additional disk composed of Very Small Particles (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] $\mu$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, mid-way between Follette et al. (2015) work (9500 K) and Brown et al. (2012a) work (9000 K), who also estimated its spectral type to be A0${}_{-1}^{+4}$ (9000$\pm$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 $\approx$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 $\mu$m) and Mp-band (4.78 $\mu$m) – would immediately lead to a dramatic over-brightness of the main PAH emission features between 4 and 15 $\mu$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 $\approx$1 AU area from the star where this disk was usually needed to account for the 1-3 $\mu$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 (Rv* taken in $[4.0-10.0]$ and Av taken in $[10.0-15.0]$).

3.5. New model: Very Small Particles 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 $\approx$55 AU was not depleted of such grains (Geers et al. 2007a; Maaskant et al. 2014). This however means that the bright inner rim of the PAH disk must be fully resolved, i.e. further than $\approx$7 to 10 AU, in order to account favorably to the contrast ratio between unresolved and resolved flux.

The MIR SED does not show any bump peeking at 9.5 and 18.3 $\mu$m, characteristic of silicate grains smaller than $\approx$2 $\mu$m. It does instead show a clear deep at 9.8 $\mu$m, de facto eliminating any major silicate grains contribution to the MIR flux; refer to Section 4.5 for an study on classical thermal grains 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 very small particles, we follow the prescription of Natta et al. (1993) and use a mixture of carbonaceous VSG and PAH with a ionized/neutral fraction parameter, from Li & Draine (2001); Weingartner & Draine (2001); 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 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 (LTE). The quantum heating routines of MCFOST were benchmarked and validated by Camps et al. (2015).

4. Resulting Structure

This best solution is obtained for a scenario where the innermost 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, it improves dramatically the fit on the VIS2 data, which it explains entirely: Figures 3 show the comparison between the VIS2 data in Ks-, Lp- and Mp-band, and both models.

4.2. Central Star

The solution of a $L_{*}$=42 L*⊙* star with a 3.7$\times 10^{-10}$ M*⊙* 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 $\mu$m and 18 $\mu$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*⊙* and an effective temperature of 9250 K to find a larger stellar radius $R_{*}$=2.5 R*⊙, rather than 1.8 or 1.4 R⊙* 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 Av=12.9 and Rv=6.5. IRS-48 is 1.4 V-magnitude more extinct, with a more gray Rv than the extinction inferred in Brown et al. (2012a) (Av=11.5 and Rv=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 Ak=1.8 is consistent with that of $\approx$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.

The radiative modelling shows that there is no extinction due to the circumstellar disk; its $\approx$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 $\rho$ Oph.

The evolutionary tracks imply that such a star falls onto the $\approx$4 million years locus, with $M_{*}$=2.5 M*⊙*. 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 Particles 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 $\approx$13 $\mu$m, where the outer disks 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 $\approx$1$\times 10^{-17}$ kg.m*-3* (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-equlibrium“ and a temperature can be defined (i.e. low dispersion on the temperature probability distribution). It reaches particularly high values for neutral PAH grains, Tmax$\approx$1500 K and Tmin$\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: Tmax$\approx$1100 K and Tmin$\approx$400 K for ionized PAH and Tmax$\approx$550 K and Tmin$\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 mm-wavelengths observed with ALMA. We fitted the simplest centro-symmetric disk model to satisfy the FIR and mm-wavelength SED.

This centro-symmetric approach is however not void of meaning. Indeed, the outer disk is optically thin in its emission window (for wavelengths $\gtrsim$13 $\mu$m, $\tau$ in the mid-plane $<$10*-1*, 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 mm-wavelength SED, the outer-radius of the disk is found at $\approx$250 AU, similarly to Maaskant et al. (2014) (225 AU), although SED fitting is not very sentivite to this parameter. For that disk extend, the total mass of the outer-disk in the MCFOST model is 9$\times 10^{-6}$ M*⊙* (3 MEarth) for grains between 30nm and 4mm, with a grain-size power law of -3.5. This mass is similar although lower than Bruderer et al. (2014) (1.6$\times 10^{-6}$ M*⊙, 5.3 MEarth*) and 9 MEarth 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 685Ghz 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 heats the outer-disk more efficiently: as seen in Figure 6, the inner-rim temperature peaks at $T_{dust}$=148 K at 55 AU.

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

[TABLE]

When we integrate over the whole disk-model the optical thickness at $\nu$=685Ghz ($\lambda$=440 $\mu$m), we find $\tau_{685Ghz}$=0.19, closely consistent with the previous 0.14 found independently using ALMA flux. Using $M_{dust}=21\times\tau_{685Ghz}\sqrt{a_{max}}$ from Draine (2006); van der Marel et al. (2013) under the same hypothesis, we find Mdust=3 MEarth for $\tau$=0.14 and 4 MEarth for $\tau$=0.19, consistent with the outer-disk dust mass in our radiative transfer model (3 MEarth).

Figure 8 shows model-images generated in the mid-IR with the corresponding imaging data, at 4.78 $\mu$m (Mp-band), 8.6 and 18.7 $\mu$m. Although 8.6 and 18.7 $\mu$m images were not included into 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 $\mu$m image is mostly unresolved, with most flux arising from PAH in the inner 30 AU, while the 18.7 $\mu$m image is highly resolved, with most flux arising from the inner rim of the outer-disk located at $\approx$55 AU.

4.5. Classical Thermal Grains in the VSP-ring

Figure 9 compares intensity profiles along the semi-major axis of the observed and modelled disks at 18.7 $\mu$m. Several models were computed, where the amount of thermal grains (silicates (70%) & carbonaceous (30%)) with sizes between 0.3 and 300 $\mu$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 axi-symmetrical model cannot obviously reproduce the eastern asymmetry, it shows a distinctive dip at the star location corresponding closely to the 18.7 $\mu$m observed image. When the amount of classical thermal-grains increases, the inner-ring becomes brighter due to the efficiency of silicates grains to radiate at 18-19 $\mu$m. The central dip progressively disappears until the flux-profile is nearly flat, for a thermal-grains mass of $\approx$10 MVSP.

Note that the VIS2 data fitting was not affected by the adding of Classical Thermal Grains (CTG), which emit a negligible amount of flux in Ks-, Lp- and Mp-band. Indeed, the emission of such CTG is mostly visible through the emission “bumps” of the silicate grains, at $\approx$9.5 and 18 $\mu$m where the SED is dominated by the outer-disk. Until an addition of 5 MVSP worth of CTG in the VSP-ring, the fit to the SED do not show any variation. From 10 MVSP however, the model starts to clearly over-estimate the SED fluxes between $\approx$8 and 15 $\mu$m. Refer to Section 6.3 for a discussion.

5. Closure Phases 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 Lp-band, -18 and +15*∘* in Mp-band, and slightly less in Ks-band. This highlights 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. Modelling

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, 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 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 $\mathcal{V}$ at the (U,V) coordinates for N point-sources is expressed by:

[TABLE]

where $\lambda$ is the central wavelength of the band. $\delta_{i}$, $\alpha_{i}$ are the angular separation and position angle of the $i^{th}$ point-source with respect to the central star (projected along declination and right ascension axes), and $f_{i}$ is its 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, 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 non-linear form of the complex visibility expression to evaluate, see Eq. (3), makes it impossible to determine the number of degrees of freedom; the reduced chi-square can not 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:

[TABLE]

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

We use a Markov-Chain Monte-Carlo (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:

[TABLE]

where the index $i$ iterates over each (U,V) coordinate of a dataset. 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:

[TABLE]

where k is the number of free parameter in the model and n 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:

[TABLE]

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

5.3. Results

We compute the relative probability for all models on each of the six datasets, which we tabulate as log of Eq. (7) in Table 7. All six datasets 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 datasets 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 impact significantly 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 such that:

[TABLE]

where the index $i$ iterates over each (U,V) coordinate of a dataset. Table 8 shows that the additional systematic error $\sigma_{systematic}$ in Ks- and Lp-band is $\approx$0.8*∘, and reaches 4∘* in Mp-band. These uncalibrated systematic errors are consistent with previous performance and slightly worse SPHERE-SAM mode: $\approx$0.5*∘* (SPHERE SAM commissioning team, private conversation). This table also shows that on all datasets, the RMS of residuals calculated using $\sigma_{i,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. unbias).

The results for the maximum likelihood parameter estimation for the favored model with a disk and two 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 wavelength. The first point-source is $\approx$105 mas on the west of the star with a contrast ratio of $\approx$3.3 mag in L-band; the second point-source location is somewhat noisier due to its proximity to the star although still consistently found at $\approx$45 mas on the north-east with a contrast ratio of $\approx$2.5 mag in L-band. These locations correspond closely to the inner-rim of the VSP-disk.

5.4. Point-sources 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 (three point-source parameters and an additional diameter), Gaussian blobs (three point-source parameters and a semi-major axis) and two-dimensional Gaussian blobs (three point-source parameters, two semi-major axes and a position angle) instead of point-sources. However, the data was not able to constrain well the additional parameters: the point-sources are mostly unresolved (i.e. angular extension $\lesssim$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 black-body 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, 10000] 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 $\approx$1350 K, see Figure 11. However, one-sigma 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 modelling at the unprojected point source location, $\approx$12.5 AU; refer to Figure 6.

The second point-source shows a somewhat lower temperature, most likely of $\approx$1000 K. Here again, one-sigma errors-bars yield several hundreds Kelvin large: from 600 to 1300 K, refer to Figure 11. Due to its proximity to the star, at or below 0.5$\lambda/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 intrisic 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 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); 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 $\approx$14.2 AU, and a height of $\approx$3.5 AU (see below). This over-brightness location corresponds to the densest region of the VSP-ring in our radiative transfer modelling, as Figure 5 tells. However, the nature of the over-brightness is not clearly constrained. It may be 1) an orbiting 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 10 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 semi-major axis of $a$=14.2 AU. It lies flat 3.5 AU above the mid-plane (hence inclination $\approx$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 semi-major 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 un-projected Keplerian velocity of 12.7$\pm 1.9$ km.s*-1* in average over the four epochs (similar values are obtained for other intervals, between epochs 1-3 and 3-4); the total orbital period at this pace is $\approx$34 years. The theoretical Keplerian velocity at such distance and given the stellar mass is $\nu=\sqrt{G\times m_{*}/r}$=12.5 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 $\approx$1 Myr (Luhman & Rieke 1999) and the stellar age is widely decreased, even though a fair difference remains that cannot be explained to date, as it 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*⊙/yr by (Salyk et al. 2013; 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⊙* dust in the inner-most disk, and 8.0$\times 10^{-10}$ M*⊙* PAH until 50 AU), or our new model (3.7$\times 10^{-10}$ M*⊙* VSP at 11-26 AU). Moreover, the disk was found to be gas-depleted until 20 AU, which puts a very hard constrain on how much material is to be found in the first tens of AU, 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.median disk lifetimes was estimated to be between 2 and 3 Myr, with an exponential decay of 2.5 Myr characteristic-time (Mamajek09). Studies show that an inner disk is detected through NIR excess in 60 to 80% of stars younger than 1 Myr while no more than 10% of stars older than 10 Myr seem to possess one (Strom89).

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 $\approx$20% of 127 K- and M-stars are surrounded by optically thick disks while none of 30 F- and G-stars showed such disk for wavelength $\approx$16 $\mu$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$\geqslant$2 M*⊙* stars aged $\geqslant$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); van der Marel et al. (2016) find a $\approx$30 AU CO-ring with small radial extent, compatible with a Keplerian velocity for a 2.0 M*⊙* star with an inclination of 50*∘*. In addition, they find a gas-depleted region inside $\approx$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 AU. 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 configuration.

6.2.2 VSP-ring Inner Truncation, Binarity

Ghez et al. (1993); Leinert et al. (1993); 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 semi-major axis of order a few AU 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 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*⊙* and a mass of $\approx$2.2 M*⊙. This luminosity remains small compared to that of IRS-48 (45 L⊙*): any stellar-type cooler than K0 would remain a nearly invisible companion for the visible and infrared instruments, especially at $\approx$25 mas typical spatial separations.

An indirect way to detect late-type companions is to check if such A-type stars emit in Xray. Stelzer et al. (2006b, a) show that, given that A-type stars are not expected to have a corona and emit Xray, a high-energy emission from these stars is generally associated to late-type binarity. In the case of IRS-48, neither ROSAT and XMM archive data nor Chandra (Imanishi et al. 2001) showed Xray emission at this location. Although soft Xray would be absorbed by interstellar absorption, hard and extreme ones should be detected. Indeed, several very young stars (even Class I) are detected in Xray in the core-region of $\rho$ Oph by Ozawa et al. (2005). Given the completeness of several Xray studies as showed in their Figure 6, and the Av=12.9mag 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*⊙* (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 velocities field it imposes on its circumbinary disk. Brown et al. (2012a) find a 2 M*⊙* to explain a Keplerien 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*⊙* and below.

IRS-48 is most probably not a binary star. However, a low-mass companion with M$\lesssim$0.5 M*⊙* and L$\lesssim$0.01 L*⊙* would be qualitatively compatible with Xray 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 mm-asymmetry observed in the southern part of the disk, a 2 M*⊙* 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 MJup.

If we carry similar first-order Hill radius estimates with the new model disk-structure and a higher-mass star, we speculate a Mp$\approx$3.5 MJup planet at a radius of a_{p}$$\approx40 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 mass between 1 and 3 MJup, 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 classical thermal grains, 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 MVSP) worth of settled classical thermal grains 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 $\mu$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 semi-major axis of the disk. This result shows that the VSP-ring most probably does not contain more thermal grains than 5$\times$MVSP$\approx$2$\times 10^{-9}$ M*⊙*.

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 PAH 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 to 10 Å (20 to 100 carbon atoms); VSG span 10 to 20 Å (100 to 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 (CVSP/H). If we take updated values for PAH and VSG abundances from Tielens (2008) (their Table 2), we find CVSP/H=3.5%. Our model showing a VSP-to-dust mass fraction $>$20% in the VSP-ring between 11 and 26 AU clearly highlighs a depletion of a factor of $\approx$5-6 of classical dust grains up to 0.3mm compared to very small particles.

We are left with a problem 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 power index of $\approx$-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 likely are 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*-1*; 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 $\approx$1000 Å.

The orbiting velocity of the first point-source in the VSP-ring around IRS-48 is $\approx$12.7 km.s*-1*, similar to the Keplerian velocity at that distance for the 2.5 M*⊙* star. Moreover, A0 stars do not generally show strong winds or emit high-energy radiations. 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 classical thermal grains 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 W.m*-2* at 1 AU. The radial acceleration from radiation pressure produced on a particle at radius r is obtained with ap=W Seff c*-1*.r*-2*, with Seff the effective cross-section of the particle and c 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 (resp. ionized PAH, 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*⊙* shows a somewhat lower but still overwhelmingly strong radiation pressure compared to gravitational potential for these species: 16, 17 and 7 times larger in average for neutral, ionized PAH and VSG.

The fact that the VSP-ring was observed to be consistently similar at a radius between $\approx$11-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 presure – 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).