arXiv:1704.05675v1 [astro-ph.SR] 19 Apr 2017 utseisteen(.. a okle l 04.Multi-e th 2004). of al. distribution et spatial Boekel van the (e.g., on therein and species 2015), dust Leine al. (e.g., et disks Menu Telescope abo protoplanetary 2004; of Large information structure Very of and the size wealth the at a provided 2003) (VLTI) al. Interferometer et Leinert (MIDI, Interferometric yo ment Mid-infrared of The disks (YSOs). planet-forming objects requ of stellar resolution regions angular innermost the the achieve resolve to possibility the interfero t fers infrared for Long-baseline conditions formation. initial planet the trial aroun determines disks circumstellar stars of sequence part main inner the of structure The Introduction 1. 07(I .Msn) 9.-00(I h aza,ad092.C and Ratzka), Th. (PI: Ratzka). 090.C-0040 Th. (PI: Mosoni), L. prog (PI: the under 1007 (Chile) Observatory Paranal at terferometer etme 5 2018 25, September Astronomy ⋆ i-nrrditreoercvraiiyo GTu impli Tau: DG of variability interferometric Mid-infrared ae nosrain aewt h S eyLreTelescope Large Very ESO the with made observations on Based .Varga J. stars eevd accepted Received; xlie ydsymtra ihaoeteds ln,whose words. plane, Key disk the above high material dusty by explained ftemdifae-mtigrgo erae rm1 from changed decreased feature region silicate mid-infrared-emitting the u the of highly of strength is The feature disk. spectral outer disk the outer and inner the between ( Conclusions. Results. Methods. Aims. Context. ne ik eprtr neso nteds ufc n a and surface disk the on inversion temperature a disk, inner msin euepyial oiae aitv rnfrm transfer radiative motivated tem physically by use spectra We mid-infrared mod emission. the disk model geometric also We a disk. fit We the 2014. about and 2011 mi between the epochs six of variability at temporal the study We processing. dust air t i-nrrdtemleiso ssrnl time strongly is emission thermal mid-infrared its havior: temporarily. 8 7 6 5 3 4 1 2 r ako tr nZei,064 Zselic, in Stars of Park nttt fAtooy aige od abig B OHA, CB3 Cambridge Road, Madingley Astronomy, of Institute ednOsraoy ednUiest,NesBhwg2 23 2, Bohrweg Niels University, Leiden Observatory, Leiden nttt o hoeia srpyis edlegUnive Heidelberg Astrophysics, Theoretical for Institute nttt o Physics for Institute a lnkIsiuefrAtooy öisul1,D-69117 17, Königstuhl Astronomy, for Institute Planck Max e-mail: Ear Hungary and Budapest, Astronomy for Centre Research Observatory, Konkoly nttu orSernud,K evn eetjela 2 Celestijnenlaan Leuven, KU Sterrenkunde, voor Instituut > 1 − elo o h esnfrteseta aiblt thg sp high at variability spectral the for reason the for look We & u pcr hwacytlieslct etr nemissio in feature silicate crystalline a show spectra au) 3 h ne ik( disk inner The 1 GTui o-aspemi eunesa,woesrnl a strongly whose star, sequence pre-main low-mass a is Tau DG [email protected] nrrditreoer a ptal eov h thermal the resolve spatially can interferometry Infrared .É Gabányi É. K. , Astrophysics Dullemond o h rgno h bopinw ics orpsil exp possible four discuss we absorption the of origin the For rtpaeaydss–sas r-ansqec tr:in stars: – sequence pre-main stars: – disks protoplanetary / GM AIGa,Uiest fGa,Uiesttpaz5 Universitätsplatz Graz, of University Graz, NAWI IGAM, r 5 aucitn.main no. manuscript h Henning Th. , < 1 1 / rz,H77 slcifld Hungary Zselickisfalud, H-7477 hrsz., 2 .Ábrahám P. , − u pcr xii 10 a exhibit spectra au) 3 h ne-ikstructure inner-disk the 4 .Ja W. , 1 .Chen L. , as088.C- rams ff . 5t 0 to 15 er of- metry e 6 vral,ee unn h 10 the turning even -variable, te al. et rt .Juhász A. , Instru- rdto ired pre- d erres- ymr hnafco fto ewe 01ad21 h half-li the 2014 and 2011 Between two. of factor a than more by -nrrditreoercsga,osre ihteVLTI the with observed signal, interferometric d-infrared µ st,Abr-eel-tas ,D610Hiebr,Ger Heidelberg, D-69120 2, Albert-Ueberle-Strasse rsity, dln oitrrttemdifae nefrmti spe interferometric mid-infrared the interpret to odeling 1 -0086 poch uulaogTTuisas h i-nrrdvraiiyis variability mid-infrared The stars. Tauri T among nusual ABSTRACT bopinfauerltdt mrhu iiaegrains. silicate amorphous to related feature absorption m .Kóspál Á. , iaindinrgoer.Teslct msini h out the in emission silicate The geometry. inner misaligned . ung lt tigt hrceietepol n iedependence time and profile the characterize to fitting plate au. 7 0,30,Lue,Belgium Leuven, 3001, 00D, In- ascncag ihtm,psil u otruec nted the in turbulence to due possibly time, with change can mass ut 3C edn h Netherlands The Leiden, CA 33 e edleg Germany Heidelberg, ,smlrt h pcr fcmtHl-op h tiigdi striking The Hale-Bopp. comet of spectra the to similar n, lt h bevditreoercsga ooti spatial obtain to signal interferometric observed the to el hSine,HnainAaeyo cecs ..Bx6,H-1 67, Box P.O. Sciences, of Academy Hungarian Sciences, th ta eouinada utpeepochs. multiple at and resolution atial msino h icmtla ik logvn informatio giving also disk, circumstellar the of emission UK niciainageof angle inclination an ihasgicn crto ae( rate accretion significant a with tri h arssa omn ein tadsac fabout of luminos- distance stellar a a with at 0 1987), of region, Thaddeus ity & forming (Ungerechts star pc Taurus 140 the in star technique this far. but so disk, explored inner physical little the been temporal of the dynamics the of changes and study cesses explore enabling to structure, used disk the be could sequences interferometric n hs bevtos hyfuda ue aiso 70 of radius outer an found they observations, these ing GTuwt h obndAryfrRsac nMillimeter- in 1 Research at for (CARMA) Array ref Astronomy and Combined observati wave 2014, the performed al. with (2010) et to Tau Maurri al. extending DG et 2008; shocks Isella al. bow et therein). and (Güdel erences au knots 1500 with least jet at studied well a has yr − 7 1 .Moór A. , GTu(e i.1 salwms 6 r-ansequence pre-main K6V low-mass a is 1) Fig. (see Tau DG ht hz20;Wie&Hlebad20) talso It 2004). Hillenbrand & White 2001; Ghez & White , ceigpoolntr ikehbt ofreimtcb enigmatic so-far a exhibits disk protoplanetary ccreting iiul GTu–tcnqe:itreoerc–infrared: – interferometric techniques: – Tau DG dividual: 1 , aain:acl bcrn neoe naceinheated accretion an envelope, obscuring cold a lanations: . 4 L 9 .Menu J. , ⊙ Pla&Salr20) ti lsia ar star Tauri T classical a is It 2002). Stahler & (Palla / I 00 rz Austria Graz, 8010, II, 1 .Mosoni L. , µ iiaefauefo msint absorption to emission from feature silicate m 2 h Ratzka Th. , ⋆ i = 1 28 , 8 ◦ n .Sipos N. and , ± M 3 ˙ . .vnBoekel van R. , mada 2 at and mm 3 10 = ◦ 4 ril ubr ae1o 16 of 1 page number, Article ikpsto nl of angle position disk a , . 6 × ain for cations ctra. / IIinstrument MIDI 10 many h ue disk outer The − rds a be can disk er oiae by dominated ftesilicate the of 1 8 information − . h radius ght m Model- mm. 7 7
isk. 4 ff c . about n 4 .P. C. , erence S 2018 ESO × 525 10 − e- n of ons − 0au, 80 7 pro- M has in ⊙ - A&A proofs: manuscript no. main
ing interferometric observations of DG Tau with the VLTI/MIDI 30 using UTs in several epochs. All observations were performed in HIGH-SENS mode with the low-resolution (λ/∆λ 30) prism (Chesneau 2007; Ratzka et al. 2007). The measurements≈ were 15 taken at various projected baseline lengths (B) and position an- gles (φB), as summarized in Table 1. The values of the position angles fall into two narrow intervals, one between 30◦ and 45◦ 0 DG Tau and one 80◦. The projected baseline lengths range from 33m to 89 m.∼ Fig. 3 (panel a) shows the distribution of the baselines in the uv-plane. Since most observations were taken with similar position angles but with different baseline lengths (U1-U2 and -15 U1-U3 configurations), we are able to sample the disk at 30◦ to 45◦ on different spatial scales. Additionally, we have high-
Arc Seconds resolution spatial information at a direction of 80 (U2-U4 ∼ ◦ -30 configuration). In all epochs, the obtained data sets consist of the N-band (7.5 13 µm) low-resolution spectrum (hereafter total spec- trum,−F ), and the interferometric spectrum (hereafter corre- -45 tot,ν DG Tau B lated spectrum, Fcorr,ν) of the target, both in the same wavelength range with the same spectral resolution. Visibility is defined as 2000 au V = Fcorr,ν/Ftot,ν. There is one exception, on 2014 January 15: -60 Correlated spectra were taken in three consecutive time slots, but only two total spectra were recorded. The majority of our 30 15 0 -15 -30 -45 observations have good quality. The observations taken on 2012 Arc Seconds February 4 had to be repeated the following night; according to Center: R.A. 04 27 04.69 Dec +26 06 15.7 the observation log the errors were too high. Therefore we dis- card this observation from the further analysis. On 2012 Febru- Fig. 1. DG Tau (at 0,0) as seen by the Sloan Digital Sky Survey (SDSS, ary 6, the total spectrum was not of adequate quality, thus it is York et al. 2000) in optical wavelengths. The red dot in the center is also excluded. 1′′ wide, which is the effective field of view of the VLTI/MIDI. An accompanying object, DG Tau B is located at the lower right corner, marked by a cross. This image is a false-color composite created by us 2.2. Data reduction using u, g, r, i, z images. For the interferometric data reduction we used the Expert Work Station (EWS) package 2.01, which is based on a coherent lin- φ = 135◦ 21◦, and a disk mass between 0.009 and 0.07 M ear averaging method (Chesneau 2007; Burtscher et al. 2012), depending± on the model parameters. ⊙ and is one of the standard tools for processing MIDI data. DG Tau is known to have dramatic variability in its 10 µm EWS routines were called from a python wrapper developed by silicate feature (e.g., Bary et al. 2009, and references therein). A Menu et al. (2015). unique aspect of this variability is that the silicate featurewas de- We applied the “direct flux calibration” method as described tected in absorption as well as emission, and during a few epochs in Burtscher et al. (2012) instead of visibility calibration. In this it could not be detected at all. Althoughseveral models have been way one can avoid the usage of the less accurate total spectrum proposed (e.g., Bary et al. 2009; Sitko et al. 2008), the origin of measurements when calibrating the correlated spectrum of the variability is not yet understood. target. The higher uncertainty level of the total spectrum is due Here we present a study of the variability of the silicate fea- to the strong and variable background in the MIDI wavelength ture in DG Tau, based on our multi-epoch interferometric data range. When observing the interferometric signal, the incoherent set. Our aim is to investigate the spatial distribution of dust background noise cancels out. However, to subtract the back- species in the disk, the variation of the silicate feature at different ground in the total spectrum chopping was used, which may spatial scales, and the temporal variability of the inner disk. In leave residual sky signal. Sect. 2 we describe the observations and the data reduction. In The calibrated correlated spectrum of the target can be ob- Sect. 3 we show the results of the interferometric modeling, es- tained from the observed (raw) correlated spectrum of the tar- timate disk size for each epoch, and model mid-infrared spectra get by dividing it with a transfer function, which following by template fitting. In Sect. 4 we delineate and discuss several Menu et al. (2015), can be written as: physical processes which can explain our findings. Finally, in Ccal Sect. 5, we summarize our results. T = raw , (1) corr FcalVcal cal 2. Observations and data reduction where Craw is the observed (raw) correlated spectrum of the cal- ibrator source, Fcal is the known total spectrum of the calibra- 2.1. Observations tor, and Vcal is the visibility of the calibrator (calculated from its known diameter). Vcal is very close to unity, since typically The VLTI/MIDI instrument (Leinert et al. 2003) combines the point-like, unresolved sources are chosen as MIDI calibrators. signal from two 8.2 m Unit Telescopes (UT) or from two 1.8 m To determine T , it is preferable to use all calibrators observed Auxiliary Telescopes (AT). It provides spectrally resolved inter- corr ferometric data in the wavelength range of 7.5 13 µm. Between 1 The EWS package can be obtained from: 2011 and 2014, we performed a dedicated program− by conduct- home.strw.leidenuniv.nl/~nevec/MIDI/index.html
Article number, page 2 of 16 J. Varga et al.: Mid-infrared variability of DG Tau with VLTI/MIDI
Table 1. Overview of VLTI/MIDI observations of DG Tau. B is the projected baseline length and φB is the projected position angle of the baseline (measured from North through East). The resolution is the approximate diameter of the beam at 10.7 µm, converted to physical scale. In the last two columns we list the name of the calibrators and the time of their measurement.
Target Calibrator
Date and time Telescopes B φB Resolution Seeing Airmass Name Time (UTC) (m)(◦) (au) (′′) (UTC) 2011-10-1006:37 U1-U2 33 30 5.7 0.8 1.7 HD27482 06:51 2011-10-1007:43 U1-U3 76 45 2.5 0.8 1.6 HD27482 07:29 2011-12-1301:53 U1-U3 56 36 3.4 0.8a 1.8 HD27482 01:32 2011-12-1302:59 U1-U2 37 34 5.1 0.7 1.6 HD27482 02:38 HD27482 03:19 2012-02-04 01:02b U1-U2 45 40 4.2 1.5 1.6 HD27482 00:42 HD27482 01:23 2012-02-0500:46 U1-U2 44 39 4.3 1.1 1.6 HD27482 00:29 HD27482 01:02 2012-02-06 00:39c U1-U3 83 46 2.3 0.9 1.6 HD27482 00:22 2012-11-0305:07 U1-U3 64 41 2.9 0.9 1.7 HD27482 04:29 2014-01-1402:39 U1-U3 88 45 2.1 0.8 1.7 HD20644 01:17 HD25604 02:27 2014-01-1501:47 U2-U4 89 83 2.1 1.5 1.6 HD20644 00:42 HD25604 02:43 HD20644 02:03 2014-01-1502:17 U2-U4 89 80 2.1 0.9 1.6 HD20644 00:42 HD25604 02:43 HD20644 02:03 2014-01-1502:22 U2-U4 89 80 2.1 0.7 1.6 HD20644 00:42 HD25604 02:43 HD20644 02:03
Notes. (a) The seeing for this epoch was missing from the observation data. We indicate here an interpolated value. (b) According to the ESO observing log, this observation had to be repeated the following day because of its poor quality. Therefore we do not use these data in our work. (c) The total flux observation was not of adequate quality so it was discarded. on the same baseline as our target in a given night. Only in the surements; and 2012 November, when only one total spectrum last epochs could we use more than one calibrator. In the previ- was measured. Because of the atmospheric ozone absorption fea- ous epochs we used 1 2 measurements of the same calibrator ture, data points between 9.4 µm and 10.0 µm have significantly − (Table 1). Tcorr is a time- and airmass-dependent quantity. Us- higher uncertainties than in other parts of the total spectrum, ing the routines of Menu et al. (2015) the time-dependence of therefore we decided to discard points of the total spectra in this Tcorr was taken into account by applying linear interpolation. To wavelength regime from the analysis. correct for airmass (X) dependency, Tcorr values were multiplied In the case of correlated spectra, we average only those ob- by the correction factor exp( fλX), where fλ is the wavelength servations that have the same baseline configuration. On 2014 dependent atmospheric extinction. Calibration of the total spec- January 15 three consecutive correlated spectra were observed trum of our target was conducted in a similar way as the corre- with the same baseline length, and position angles within a few lated spectrum, but using Vcal 1 in Eq. (1). In this case the ≡ degrees, so we averaged these measurements. Thus we ended up transfer function is the atmospheric transparency. with nine independent correlated spectra (see Table 4). The outputs of the data reduction pipeline are the calibrated total spectrum and the calibrated correlated spectrum of the tar- 2.3. Calibration quality get. Visibilities are also calculated as the correlated spectrum divided by the total spectrum. The correlated spectrum can be MIDI spectra are known to have calibration biases interpreted as the spectrum of the very inner part of an object. In (Burtscheretal. 2012). The uncertainties on the spectra the case of DG Tau this means that we have spectral information provided by the EWS seem to be much larger than the variance with a resolution of 2 6 au, depending on the baseline length − of the neighboring spectral data points (see Fig. 7 for example). (see Table 1). This suggests that the main error source is systematic, due to The accuracy of the total spectra is not adequate to detect calibration issues, that is, the absolute level of the spectrum variations from one night to the next, therefore we averaged the can only be calibrated with a higher uncertainty, although the calibrated total spectra taken on consecutive nights to increase spectral shape is more reliable. The calibration biases arise from the signal-to-noise ratio. This resulted in five total spectra: 2011 flux losses in the MIDI instrument (Burtscher et al. 2012). October, 2011 December, and 2012 February, obtained by aver- To verify the calibration quality, we reduced calibrators of aging two measurements; 2014 January, by averaging three mea- DG Tau as if they were science targets, using other calibrators
Article number, page 3 of 16 A&A proofs: manuscript no. main to calibrate them. A calibrator is supposed to be non-variable, and compact, so it remains mostly unresolved by the interferom- 4.0 eter. From Table 1 we can see that all observations from 2011 and 2012 were done using the same source (HD 27482) for cal- ibration. Data from 2014 were calibrated with two other stars 3.5 (HD 20644, HD 25604). We obtained five correlated spectra, (Jy) 3.0 and three total spectra for HD 27482. The standard deviation m µ of these spectra at 10.7 µmis0.4 Jy (7%) both for the correlated and total flux density. HD 20644 was observed three times on 2.5 corr, 10.7
two consecutive days. The standard deviation for the correlated F flux densities (again at 10.7 µm) is 1.7 Jy (10%), while for the 2.0 total flux densities one finds 2.0 Jy (12%). We take these standard deviations as the errors of the abso- 1.5 lute calibration. The error bars on the calibrator spectra calcu- lated by the reduction pipeline are relatively similar: 7 10% for 1.0 the correlated spectra and 10 14% for the total spectra.− Thus 0 20 40 60 80 B (m) we conclude that the uncertainties− provided by the EWS are reli- eff able. In the case of DG Tau, typical measurementerrors are 6% ( 0.15 Jy at 10.7 µm) for the correlated spectra and 10 ∼20% ∼ − 0 100 200 300 400 500 600 700 800 (0.5 1.0Jyat10.7 µm) for the total spectra. Observation date (days since 2011−10−10) − Fig. 2. 10.7 µm correlated flux densities as a function of the effective 3. Results baseline length (Beff, see text) for DG Tau, showing the result of the model fitting (blue curve). The symbols are color coded for observation 3.1. Interferometric modeling date. The visibilities and correlated spectra contain spatial information about our target. Because of the sparse uv-sampling and lack of Ftot,ν only includes the flux from the disk. At mid-infrared wave- closure phases, we cannot obtain an image of the disk, instead lengths the stellar photospheric flux is negligible. We estimate we can try to fit simple disk models to the interferometric data. that the stellar flux is at most 3% of the disk flux, also taking Note that we use correlated spectra for the fitting, instead of visi- into account the possible contribution from accretion hot spots bilities, because the latter have much larger measurement uncer- on the stellar surface. tainties. At a specific wavelength, we have correlated flux densities Theaim ofthe modelfitting is to measurethe size of the mid- measured as a function of the projected baseline length. The infrared emitting region of the DG Tau disk. We fit the measured brightness profile of the model disk should be transformed in the correlated flux densities with a physically motivated geometric uv-space. This transformation is known as a Hankel transform, disk model also used by Menu et al. (2015). For the disk incli- which is essentially the two-dimensional Fourier transformofa nation and position angle we adopted the values of Isella et al. circularly symmetric function, with the following formula: (2010) (Sect. 1). The model geometry is a thin, flat disk, which starts at the Rout rBν (T (r)) J0 (2πrBeff/ (λd)) dr dust sublimation radius and extends out to Rout = 300 au, where Rsub Fcorr,ν (Beff) = Ftot,ν R . (6) the mid-IR radiation is negligible. The temperature drops as a Rout rBν (T (r)) dr Rsub power-law, and the disk emits thermal radiation: R Iν (r) = τνBν (T (r)) , (2) Here J0 is the 0th order Bessel function, d is the distance from the source and Beff is the effective projected baseline length, cor- where τν is the optical depth, Bν is the Planck-function and the rected for the inclination (i) and position angle (φ) of the disk as temperature (T) as a function of radius (r) is given as follows: q r − T (r) = Tsub . (3) α = atan2 (sin (φB φ) , cos i cos (φB φ)) , (7) Rsub ! − −
Here Tsub is the dust sublimation temperature, fixed at 1500 K. Bu, ff = B cos α cos φ B cos i sin α sin φ, (8) Rsub is the sublimation radius, calculated from the known lumi- e − nosity (L⋆) of the central star:
1/2 Bv,eff = B cos α sin φ + B cos i sin α cos φ, (9) L⋆ Rsub = . (4) 4πσT 4 sub B = B2 + B2 . (10) With L = 0.9 L , the sublimation radius is R = 0.07 au. eff u,eff v,eff ⋆ sub q The model has two⊙ free parameters: q, which is the power-law Fig. 2 shows a fit to our data at λ = 10.7 µm, obtained by exponent of the temperature profile, and τ , the optical depth. ν minimizing the reduced χ2 using the measurement uncertainties The latter can be related to the total flux density of the object provided by the data reduction pipeline. The exponent q deter- (F ), which is also measured by MIDI, as follows: tot,ν mines how rapidly the disk fades outwards, thus it can be used as a measure of the compactness of the emission. The best-fit value Rout for q is 0.54 0.02, which closely matches what one would ex- F = πrτ B (T (r)) r. tot,ν 2 ν ν d (5) pect for the± optically thin emission from a warm disk surface ZRsub Article number, page 4 of 16 J. Varga et al.: Mid-infrared variability of DG Tau with VLTI/MIDI
Table 2. Results of the interferometric modeling, for three wavelengths, baselines on a given night. However, the total spectrum measure- calculated from the correlated spectra. ments (corresponding to the source at zero baseline length) can be used to determine sizes for each epoch separately. Such cal- culations are generally less accurate than correlated flux mea- λ Ftot,ν q re re surements at multiple baselines, but still enable us to study time (µm) (Jy) (mas) (au) variability of the disk size. 8.0 3.3 0.1 0.58 0.03 2.5 0.3 0.35 0.05 We utilize here almost the same model as in Sect. 3.1, but λ = . µ 10.7 4.0 ± 0.2 0.54 ± 0.02 5.2 ± 0.4 0.73 ± 0.10 we use visibilities measured at 10 7 m instead of corre- lated flux densities. The resulting disk sizes are shown in Table 3 13.0 5.8 ± 0.4 0.53 ± 0.02 8.2 ± 0.9 1.15 ± 0.19 ± ± ± ± and in Fig. 3 (panel b). It is a polar plot showing the measured sizes at specific baseline position angles. Note that the values in the table are corrected for inclination, but the plotted values in Table 3. Results of the interferometric modeling calculated from the the figure are not. Assuming that there was no significant change visibilities for each epoch (for λ = 10.7 µm). in the brightness distribution from one night to another, we can make some constraints on the disk inclination using the size val- ues from the last two epochs (shown in purple in Fig. 3 panel V V V Epoch φB q re re b). The gray ellipse on the plot shows the shape of the projected (◦) (mas) (au) disk from mm observations(Isella et al. 2010, arbitrarily scaled). 2011Oct 37 0.50 0.02 7.7 1.3 1.08 0.20 From the figure we can readily see that our mid-IR shape is con- 2011Dec 35 0.49 ± 0.02 8.2 ± 1.7 1.15 ± 0.28 sistent with the millimetric shape within the uncertainties. 2012Feb 39 0.51 ± 0.03 6.7 ± 1.3 0.94 ± 0.22 With the sizes determined for each epoch, we can now look 2012Nov 41 0.52 ± 0.02 6.3 ± 0.6 0.89 ± 0.13 for the signs of temporal size variations. Although the values 2014Jan 45 0.54 ± 0.02 5.2 ± 0.7 0.73 ± 0.12 have relatively large error bars, there seems to be a decreasing ± ± ± trend in the disk size: between 2011 and 2014 the half-light ra- 2014Jan 80 0.53 0.02 5.4 0.6 0.76 0.13 V ± ± ± dius decreased from re = 1.15 0.28 au to 0.73 0.12 au. We performed a linear fit to the sizes± as a function of± time, and the resulting rate of change is 0.14 0.07 au/yr,which is a 2σ re- sult for variability. In parallel,− total± flux also seems to delineate a layer. Following Menu et al. (2015) we define re, a half-light ra- dius: decreasing trend. The size at the last epoch is consistent with the average size derived from the correlated spectra in the previous F re tot,ν = section. This means the whole disk can be described with a sin- 2πrIν (r) dr. (11) V 2 ZRsub gle temperature component. At earlier epochs, however, re > re, which suggest that the outer disk may have a shallower temper- From the best fit in Fig. 2, the resulting half-light radius at ature profile. This may imply an ongoing structural rearrange- 10.7 µm is re = 0.73 0.10 au (see Table 2 for results for other ment of the dusty disk material. Millan-Gabet et al. (2016) ob- wavelengths). The disk± appears more extended with increasing served DG Tau with the Keck Interferometer in 2010, and got wavelength (Table 2). This is because the bulk of the radiation at 1.15 0.02 au for the mid-infrared (N-band) size by fitting a larger wavelengths comes from a wider, cooler area with a larger Gaussian± model with an inner gap. This value fits the trend we average distance from the star (Schegerer et al. 2008). observed. We note, however, that their value may not be directly Due to the high accretion rate, mentioned in Sect. 1, the ac- comparable to our sizes due to the different fitting techniques. cretion luminosity of DG Tau can be also significant. Thus the From the interferometric modeling in Sect. 3.1 we have seen, total luminosity of the source can be higher than the luminosity that the difference of individual correlated flux densities from the of the stellar photosphere, which can move Rsub outwards. The best-fit model is in the range of 10 20%. This implies that real variable accretion rates yield 1 10 L accretion luminosities, variability of the half light radius− derived from the correlated − ⊙ and Rsub will be in the range of 0.1 0.2 au. There is a 40% flux densities should be in the same order. However, indicated increase in q with 10 L total luminosity,− with respect to the by the half-light radii from the visibilities (i.e., including the to- photosphere-only case. With⊙ a modest accretion luminosity of tal flux density) we can observe somewhat larger size variations 1L the change in q is only 10%. However, the resulting re val- (up to 30%) among the different epochs, indicating that the total ⊙ ues are completely insensitive to these changes in Rsub, changing spectrum is much more variable than the correlated spectrum.In by only 4% at 10.7 µm. the next section we show a detailed analysis of the mid-infrared The scatter of the correlated flux densities along the fitted spectral behavior of DG Tau. curve in Fig. 2 is 10 20%. It means that a single non-variable disk model can constrain− the average size of the mid-infrared emitting region. However, the correlated flux densities measured 3.3. Mid-infrared spectral behavior of DG Tau at the two shortest baselines (Beff < 40 m) show 20% flux dif- DG Tau is famous for its highly variable silicate feature in its ∼ ference over two months, thus hinting at temporal variability. In mid-infraredspectrum (e.g., Woodward et al. 2004; Kóspál et al. the following section, we elaborate on this possible size variabil- 2012). As early as in 1982 and 1983, Cohen & Witteborn (1985) ity. observed the mid-infrared spectrum of DG Tau and found that it showed the silicate feature in absorption. Later, Wooden et al. 3.2. Size variability (2000) reported that it showed silicate in emission in 1996 (Wooden et al. 2000). At that time, amorphous and crystalline In the previous section, we deduced the average disk size using forms of olivine and pyroxene were identified, while in the fol- all the correlated spectra. Nevertheless, we cannot determine in- lowing year (in 1997 September) the mineralogy was domi- dividualsizes for each epochdue to the limited range of observed nated by Mg-rich crystalline olivine. In 1998, the silicate fea-
Article number, page 5 of 16 A&A proofs: manuscript no. main
0 0 a) b) 30 330 30 330
60 300 60 300
100 1.5 80 60 U2U4 1 40 20 0.5 90 270 90 270
U1U2
120 U1U3 240 120 240
150 210 150 210 180 180
uv¡ plane (m) Image plane (au)
0 200 400 600 800 ¢ Observation date (days since 2011 10 10) Fig. 3. a) uv-plot of DG Tau MIDI observations. The tracks correspond to specific telescope configurations. b) Disk half-light radii, determined from visibilities for each epoch. The symbols are color-coded for observation date. Position angles are measured from North to East (counterclockwise). The radial scale of the plot is given in au. ture turned into absorption, and a drop in visual magnitude was observed prior to this event (Wooden et al. 2000). Although with different shape, the silicate feature was also in absorption 11 in 2001, and this was simultaneous with an overall brightness SWS 1997 Sep 18 10 increase in the 3 5 µm wavelength range. The feature dis- TIMMI2 2002 Dec 25 appeared by 2001− October (Woodward et al. 2004). A Spitzer 9 spectrum in 2004 (Kóspálet al. 2012) indicates a weak sili- 8 cate emission. According to Sitko et al. (2008), in 2006, DG 7 Tau again experienced an outburst in the 3 5 µm range and 6 − the silicate feature appeared in emission in its spectrum again. 5
Bary et al. (2009) presented nine Spitzer observations of DG Tau Flux density (Jy) 4 obtained between 2004 and 2007. The silicate feature showed 3 variations over month- and year-long timescales; shorter varia- 2 tions on day- to week-long timescales were not detected. The spectra were dominated by crystalline forsterite in emission in 1 0 all epochs. Recently, Millan-Gabet et al. (2016) reported mid- 2 3 4 5 6 7 8 9 10 11 12 13 infrared and near-infrared measurements of DG Tau performed Wavelength (micron) in 2010 September, October, and November with the Keck In- terferometer. At that time, the spectrum did not show a silicate Fig. 4. DG Tau mid-infrared spectra taken by the SWS on 1997 Septem- emission feature. ber 18 (Frieswijk et al. 2007) and by TIMMI2 on 2002 December 25 Checking the literature, we identified two additional observa- (Przygodda 2004). From the latter the wavelength range affected by the atmospheric ozone feature ( 9 10 µm) is excluded. tions of the silicate feature of DG Tau (see Fig. 4, Frieswijk et al. ∼ − 2007; Przygodda 2004). On 1997 September 18, a spectrum taken by the Short Wavelength Spectrometer (SWS) on-board to the 2001 October data. This might be indicative of a longer the Infrared Space Observatory2 showed weak emission. This is featureless period. consistent with the report of Wooden et al. (2000) on the detec- Using interferometric observations, it is possible to narrow tion of a similar silicate emission feature in 1997 September. The down the spatial regions responsible for the silicate emission other observation taken by the TIMMI2 instrument at La Silla (e.g., van Boekel et al. 2004). The mid-infrared measurements Observatory on 2002 December 25 showed no feature, similarly of Millan-Gabet et al. (2016) already suggested a difference be- tween the correlated and total spectra in 2010. While DG Tau 2 Data were downloaded from the ISO archive did not show emission in the total spectrum, it showed indi- http://iso.esac.esa.int/ida/ cation of silicate in absorption in the correlated spectrum. Our
Article number, page 6 of 16 J. Varga et al.: Mid-infrared variability of DG Tau with VLTI/MIDI
4
3
2
1 2011 Oct 33m 30deg 2011 Dec 37m 34deg 2012 Feb 44m 40deg <2.9 au <2.6 au <2.2 au
4
3
2
Flux density (Jy) 1 2011 Oct 76m 45deg 2011 Dec 56m 36deg 2012 Nov 64m 41deg <1.3 au <1.6 au <1.5 au
4
3
2
1 2012 Feb 83m 46deg 2014 Jan 88m 45deg 2014 Jan 89m 80deg <1.2 au <1.1 au <1.1 au 0 8 9 10 11 12 8 9 10 11 12 8 9 10 11 12 13 Wavelength (micron)
Fig. 5. Correlated spectra (black filled circles) and best fit curves (solid red line) of DG Tau at different projected baselines and in different epochs. The rows are categorized by approximate size of the observed baseline (top: short baseline, middle: intermediate baseline, bottom: long baseline). Gray shading indicates the level of total calibration uncertainties, while black error bars represent the random point-to-point uncertainties (see Appendix A for more details on error calculation).
VLTI/MIDI observations strengthen this notion. The correlated Table 4. Results of the fit to the correlated spectra. and total spectra of DG Tau have dramatically different shape. In the following we attempt to decompose the correlated spectra Epoch B bcorr acorr ccorr (Sect. 3.3.1; dominated by the inner few au region of the cir- (m) (Jy) (Jy/µm) (Jy) cumstellar disk) and the uncorrelated spectra (Sect. 3.3.2; the difference between the total and correlated spectra, which thus 2011Oct 33 3.2 0.2 0.25 0.02 0.29 0.02 76 2.6 ± 0.2 0.07 ± 0.01 0.41 ± 0.03 describes the outer, large-scale part of the disk) of DG Tau into a ± ± ± linear continuum and the silicate feature. Using the resulting pa- 2011Dec 56 3.3 0.2 0.12 0.01 0.42 0.03 rameter values we then study the changes in the silicate feature 37 4.0 ± 0.2 0.24 ± 0.01 0.36 ± 0.02 in different epochs and at different resolutions. ± ± ± 2012Feb 44 3.3 0.2 0.18 0.01 0.43 0.08 83 2.2 ± 0.2 0.031± 0.004 0.20 ± 0.02 2012Nov 64 3.3 ± 0.2 0.09 ±0.01 0.48 ± 0.03 3.3.1. Spectra of the inner disk ± ± ± 2014Jan 88 2.3 0.1 0.006 0.002 0.30 0.04 89 2.2 ± 0.2 −0.05 ±0.02 0.27 ± 0.08 The correlated spectra, shown in Fig. 5, are dominated by the in- ± − ± ± ner, compact region of the disk, whose size is determined by the resolution of the interferometer for the given baseline (see Ta- ble 1). All of them are very smooth. The main trend is a decrease upto 9.5 µm and a flat or increasing part at longerwavelengths. This shape∼ suggests that we see the silicate feature in absorption To describe the correlated spectra in Fig. 5, we model them here, and the profile and the wavelength of the minimum indi- as a combination of a linear continuum and an amorphous sili- cates amorphous grains. cate absorption. The spectra, Fcorr(λ) were fitted with the follow-
Article number, page 7 of 16 A&A proofs: manuscript no. main
Table 5. Results of the fit to the uncorrelated spectra. 0.18
Epoch B buncorr auncorr cuncorr 0.16 (m) (Jy) (Jy/µm) (Jy) 2011Oct 33 1.4 0.6 0.3 0.1 0.3 0.1 0.14 76 2.0 ± 0.6 0.5 ± 0.1 0.5 ± 0.1 ± ± ± corr 0.12 2011Dec 56 2.3 1.0 0.6 0.2 0.9 0.2 /b 37 1.6 ± 1.0 0.5 ± 0.2 0.8 ± 0.2 corr
c ± ± ± 0.10 2012Feb 44 0.9 0.6 0.1 0.1 0.7 0.2 83 2.2 ± 0.5 0.4 ± 0.1 0.3 ± 0.1 ± ± ± 0.08 2012Nov 64 1.8 0.4 0.34 0.04 0.6 0.1 ± ± ± 2014Jan 88 1.2 0.4 0.1 0.1 0.8 0.1 0.06 89 1.4 ± 0.4 0.2 ± 0.1 0.8 ± 0.1 30 40 50 60 70 80 90 ± ± ± B (m) B > 80 m do not fit in this trend. This may be a hint that the ab- sorption comes from a ring-shaped region with inner and outer Fig. 6. Amplitude of the silicate absorption feature (ccorr) normalized radii of 1 au and 3 au, respectively (these values refer to the for the 10.7 µm continuum (bcorr) as a function of the projected baseline ∼ ∼ length (B). Symbols are color-coded for observation date (color-coding resolution of the interferometer at the corresponding baselines). is the same as in Fig. 2). The presence of a general trend in Fig. 6 suggests that the main variations in the normalized absorption are related to the differ- ent spatial resolutions. Thus there is no indication of temporal ing formula: variability within the inner regions.
Fcorr(λ) = acorr (λ 10.7 µm) + bcorr ccorr G(λ), (12) · − − · 3.3.2. Spectra of the outer disk where a , b , and c are wavelength independent factors corr corr corr The uncorrelated spectrum is the emission component that is and G(λ) is the normalized spectrum of interstellar amorphous spatially resolved out by the interferometer, thus does not appear silicate grains measured by Kemper et al. (2004) towards the in the correlated spectrum. It can be computed by subtracting the Galactic Center. The correlated spectra are fitted in the wave- correlated spectrum from the total spectrum. In the case of our length range between 8 µm and 12.7 µm. We sampled a grid of VLTI/MIDI observations uncorrelated spectra are dominated by the parameter values of a , b , and c . The fits were com- corr corr corr the outer parts of the 10 µm emitting region of the circumstellar pared using the obtained χ2 values. We note that the parameters disk, outside a radiusof 1 3 au (dependingon the resolution, see were fitted simultaneously, making it possible to determine a re- Table 1). In contrast to the− previous subsection, the uncorrelated liable continuum level, even though the broad absorption feature spectra show no absorption; rather we see an emission feature would not allow to constrain the continuum alone. The best-fit (see Fig. 7), similar to the observations of Bary et al. (2009), for parameter values corresponding to the lowest χ2 are listed in Ta- example. We note that Bary et al. (2009) did not have spatially ble 4. In Appendix A we describe how uncertainties of the pa- resolved observations, that is, they studied total spectra. rameters were calculated. In Fig. 5 the best fitting F (λ) curves corr In order to describe the shape of the uncorrelated spectra, are displayed. Apart from the shortest investigated wavelengths we have to choose a spectral template for the silicate emis- (8 9 µm), the resulting curves are good representations of the sion feature. Bary et al. (2009) obtained high-resolution Spitzer measurements.− spectra of the silicate feature, indicating the presence of crys- The derived absolute flux density levels at 10.7 µm (bcorr in Table 4) are mainly determined by the projected baseline talline forsterite. Previous studies (e.g., Ábrahám et al. 2009) length as shown in Fig 2. The slope of the linear continuum showed that the silicate feature in the spectrum of comets can be very similar to that of protoplanetary disks. Therefore we chose (acorr) monotonically decreases with increasing projected base- line length. This is expected, since with better resolution (longer the continuum-subtracted, normalized mid-infared spectrum of projected baseline lengths) we observe the more compact re- cometHale-Bopp (Crovisier et al. 1997, and see redcurve in Fig. gions, which are dominated by the hotter dust emitting at shorter 8) as a template to fit the silicate emission feature. This spectrum wavelengths. At the longest baselines (83m in 2012 February also shows the presence of crystalline silicate grains. Since the and 88m and 89m in 2014 January), the derived continuum is general shape of the Hale-Bopp spectrum matches well the spec- almost completely flat. Since both a and b have a depen- trum of DG Tau, it was unnecessary to introduce any correction corr corr ff dence on B, the two parameters correlate: their correlation coef- for possible temperature di erences between the comet and the ficient is 0.87. disk. We model the spectra by the following formula: The average value of the amplitude of the silicate absorption Funcorr(λ) = auncorr (λ 10.7 µm) + buncorr + cuncorr H(λ), (13) is c = 0.35Jy with a standard deviation of σc, = 0.09 Jy. · − · h corri corr In Fig. 6, we show the amplitude of the absorption normalized with parameters auncorr (slope) and buncorr (10.7 µm flux density), for the 10.7 µm continuum (ccorr/bcorr) as a function of the pro- and the silicate emission template (H(λ)) from the spectrum of jected baseline length (B). There is a clear trend of increas- comet Hale-Bopp with an amplitude cuncorr. ing ccorr/bcorr with B, which might indicate that the absorption The derivation of the best-fit parameters and their uncer- mostly originates from a confined region in the inner part of tainties was done in the same way as for the correlated spectra the dusty disk. It is interesting to note, that ccorr/bcorr values at (see Sect. 3.3.1). We note that the uncertainties of the measured
Article number, page 8 of 16 J. Varga et al.: Mid-infrared variability of DG Tau with VLTI/MIDI
5 2011 Oct 33m 30deg 2011 Dec 37m 34deg 2012 Feb 44m 40deg 4
3
2
1
5 2011 Oct 76m 45deg 2011 Dec 56m 36deg 2012 Nov 64m 41deg 4
3
2 Flux density (Jy) 1
5 2012 Feb 83m 46deg 2014 Jan 88m 45deg 2014 Jan 89m 80deg 4
3
2
1
8 9 10 11 12 8 9 10 11 12 8 9 10 11 12 Wavelength (micron)
Fig. 7. Uncorrelated spectra (black symbols) of DG Tau at different projected baselines and at different epochs. The part affected by the telluric ozone feature, between 9.4 µm and 10 µm is not shown. Red line represents the best-fit curve using the silicate spectrum of Hale-Bopp (Eq. 13). The rows are categorized by approximate size of the observed baseline (top: short baseline, middle: intermediate baseline, bottom: long baseline). Gray shading indicates the level of total calibration uncertainties, while black error bars represent the random point-to-point uncertainties (see Appendix A for more details on error calculation).
points are much larger for the uncorrelated spectra than for the Changes in cuncorr can also be traced in Fig. 8, where we correlated spectra. This is caused by the larger uncertainties of show the continuum subtracted emission features from our MIDI the total spectra due to the high systematic background residuals data (blue curves), and from earlier Spitzer measurements of (see Sect. 2.3). In the case of the correlated spectra the back- Bary et al. (2009) (gray curves). In order to be able to compare ground noise almost completely cancels out since it is not co- the Spitzer and MIDI spectra, we only took into accountthe com- herent (Chesneau 2007). Hence, the correlated spectra are only mon wavelength range (8 13 µm). Additionally, the continua limited by the statistical photon noise. Additionally, the wave- in the Spitzer spectra were− calculated in the same way as in the length range affected by the atmospheric ozone feature (between MIDI uncorrelated spectra by fitting them with Eq. 13. (We note 9.4 µm and 10 µm) was disregarded in the fitting process. The that in general the wider wavelength range where Spitzer op- resulting best-fit parameter values are given in Table 5. The best erated allows a more precise continuum fitting, as was shown fit curves are overplotted on the data in Fig. 7. The fits repro- in Bary et al. 2009). Spitzer observed the total spectra, so to be duce the measured data points reasonably well, and the silicate able to directly compare these data with the uncorrelated MIDI emission feature is clearly detected at least at the 3σ level in all spectra, a correction for the inner disk absorption is also needed. epochs. It is quite remarkable that the spectrum of a comet is Thus we subtracted an absorption profile with an amplitude of very similar to the spectrum of a protoplanetary disk. 0.25 Jy from the Spitzer spectra. The amplitude was chosen to match with the MIDI observations with the corresponding base- lines. In both the Spitzer and MIDI datasets a large and small The average value of the amplitude of the silicate emission is silicate feature is shown. The amplitudes and shapes of the sil- c = 0.63Jy with a standarddeviationof σc, = 0.22 Jy, h uncorri uncorr icate features are relatively similar. We also show in Fig. 8 the which is larger than the value for the absorption in the correlated spectrum derived from the comet Hale-Bopp multiplied by 0.8, spectra (σc,corr) bya factorof 2.5, which means more significant corresponding to best-fit value obtained for the January MIDI variability. If we take the relative scatter (i.e., σc/ c ), then we uncorrelated spectrum (red line). still see that the variations in the emission are abouth i one third higher than the variations in the absorption.
Article number, page 9 of 16 A&A proofs: manuscript no. main
1.2
0.8 1.0
0.6 0.8 0.4 (Jy) 0.2 0.6
0 uncorr c
Flux density (Jy) 0.4 -0.2 Spitzer 2005 Oct 9 Spitzer 2006 March 17 -0.4 Comet Hale-Bopp VLTI/MIDI 2012 Feb 83m 0.2 -0.6 VLTI/MIDI 2014 Jan 88m 8.5 9 9.5 10 10.5 11 11.5 12 12.5 13 0.0 Wavelength (micron) 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
auncorr (Jy/µm) Fig. 8. Emission feature in the spectra of DG Tau. Dashed and solid blue lines show the continuum subtracted uncorrelated spectra of DG Tau Fig. 10. Amplitude of the silicate emission feature (cuncorr) plotted as in two epochs. Dashed and solid gray lines show the continuum sub- a function of the uncorrelated continuum slope (auncorr). Symbols are tracted spectra measured by Spitzer in 2005 October 9 and 2006 March color-coded for observation date (color-coding is the same as in Fig. 2). 17, respectively (Bary et al. 2009), corrected for an average absorption (see text for details). The red line represents the emission feature in the spectrum of comet Hale-Bopp, scaled to match the solid blue curve.
2011 Dec had the lowest value in 2011 October, but it almost tripled two 2011 Oct 2012 Feb 2012 Nov 2014 Jan months later. At the longest baselines we see the same amount of change from 2012 February to 2014 January. To deduce the spatial variability, we can compare observations taken at the 1.0 same epochs with different baseline configurations. According 80 to Fig. 9, we do not see any trend in the silicate emission feature 0.8 as a function of the baseline length. We note however, that the 70 uncorrelated spectra are much less sensitive on the baseline than the correlated spectra, because the uncorrelated emission comes
(Jy) 0.6 60 from a much larger area. uncorr c 0.4 50 Bary et al. (2009) presented a figure of the silicate emission amplitude as a function of the [6.0] [13.5] µm color index, based on their Spitzer data. They found− a weak correlation be- 0.2 40 tween the two quantities. Instead of this color index we can use the continuum slope auncorr to explore this relation. In Fig. 10 we 0.0 30 plot cuncorr versus the mid-infrared spectral slope auncorr, fromour 0 200 400 600 800 B (m) uncorrelated MIDI spectra. The correlation coefficient between Observation date (days since 2011-10-10) auncorr and cuncorr is ρa,c,uncorr = 0.04, thus we find no correlation between the two parameters. − Fig. 9. Amplitude of the silicate emission feature (cuncorr) plotted as a function of observation date. Symbols are color-coded for the projected Previously, in the correlated spectra we have seen a correla- baseline (B). Notice the large temporal variations of cuncorr. The gray stripe indicates the range of the emission feature strength obtained from tion between the continuum slope (acorr) and the 10.7 µm flux radiative transfer modeling, see Sect. 4.4 for details. density (bcorr). This relation is a consequence of the different spatial scales probed by the interferometer: with a longer base- line, the correlated spectrum comes from a smaller region, thus To disentangle temporal and spatial variability in interfero- having less absolute flux densities. At the same time, we see metric observations one has to consider observations taken at the hotter inner parts of the disk having bluer color and less similar projected baseline lengths at different epochs, and obser- steep spectral slope. Interestingly, we observe a similar rela- vations taken at different projected baseline lengths at the same tion in the uncorrelated spectra between auncorr and buncorr, with epoch. In Fig. 9 we plot the amplitude of the silicate emission ρa,b,uncorr = 0.86 (the same value as in the case of acorr and bcorr). feature (cuncorr) as a function of time, color-coded for projected This means that when the outer disk gets brighter, its spectral baseline length. Comparing the observations taken at similar an- slope gets steeper (corresponding to a lower average temper- gular resolutions in 2012 February and 2014 January at a pro- ature). But neither auncorr nor buncorr exhibit significant depen- jected baseline length of 85m (purple points), and in 2011 dence on the baseline length. Thus we attribute the scatter in October and December at∼ a projected baseline length of 33m these parameters to temporal variations that cause changes in the ∼ (red points) we find significant variability in the amplitude of continuum in a way that preserves the relation between auncorr the emission feature over time. At the shortest baselines, cuncorr and buncorr.
Article number, page 10 of 16 J. Varga et al.: Mid-infrared variability of DG Tau with VLTI/MIDI
4. Discussion Table 6. Parameters used in our radiative transfer model. Those in italics were kept fixed during the modeling. 4.1. Origin of the mid-infrared variability While an increasing number of studies demonstrated that System Parameters mid-infrared variability is common among young stars (e.g., Distance (d) 140pc Kóspál et al. 2012; Bary et al. 2009), DG Tau exhibits unusually Inclination (i) 28 ◦ large amplitudes up to a factor of two in the 10 µm silicate fea- Stellar Parameters ture. Moreover, the feature may even turn into absorption, which Temperature (T ) 4300K is very rarely observed in other young stars. Bary et al. (2009) Luminosity (L )∗ 1.0L developed a two-temperature two-slab model to explain the vari- Mass (M )∗ 1.0M⊙ ∗ ⊙ ability of the silicate emission feature, based on the method of Visual extinction (AV) 1.6 Sargent et al. (2009). This model includes two emitting layers on Selective extinction coefficient 3.1 the disk surface, with the top, intervening layer being cooler, en- (RV) abling self-absorption. The authors were able to reproduce their Gaseous Accretion Disk multi-epoch data set with this model by fitting the temperatures Inner radius (Rg ) 0.01au ing and dust masses. An important finding is that the intervening Outer radius (Rout) 0.07au 7 1 cooler dust layer has the same composition as the emitting re- Accretion rate (M˙ ) 10− M yr− gion, therefore the variations do not indicate any change in the Dusty Disk ⊙ dust composition. Bary et al. (2009) proposed three scenarios to d Inner radius (Rin) 0.07au explain the observed silicate variability on month- and week- d Outer radius (Rout) 100au long timescales: a) variable disk-shadowing by a warped puffed- d Inner scale height (hin) 0.04 up inner rim, b) turbulent mixing, and c) disk winds, capable of Power-law index for scale height 0.228 lifting dusty material above the disk. (ph) One of the main results of our MIDI observations is that the Power-law index for surface den- 0 silicate emission and absorption are detected in spatially distinct sity (pΣ) regions in the disk. Thus, the variability observed in the total 2 Mass (Mdisk) 3.4 10− M spectrum of DG Tau can be explained as the combination of an Gas-to-dust mass ratio × 100 ⊙ almost constant absorption towards the inner disk and a variable Dust Grain Size Distribution silicate emission from the outer disk outside 1 3au. If the emis- 2 3 − Diameter range (a) [10− , 10 ] µm sion feature completely fills up the absorption dip, no feature is Power-law index (p) 3.5 seen in the total spectrum (e.g., in the TIMMI2 observation in Dust species: − 2002, see Fig. 4), while stronger emission gives rise to the sili- amorphous carbon 50 % cate feature usually observed. astronomical silicate 50 % The other important finding of our analysis is that the shape Cloud Component of the absorption feature in the correlated spectra is character- c Inner radius (Rin) 5au istic of amorphous grains, but the emission feature in the un- c Outer radius (Rout) 7.5au correlated spectra indicates the presence of crystalline grains. Center height (hc ) 0.6 This is an unexpected result, as it seems to be in contradiction center Vertical size (∆hcloud) 0.6 with the usual picture of protoplanetary disks, where thermal an- 7 Mass (Mcloud) 3.7 10− , M nealing causes the crystallization of initially amorphous grains ⊙ 7.4× 10 7 M (Colangeli et al. 2003) in the inner disk, near the sublimation ra- − Gas-to-dust mass ratio × 100 ⊙ dius, while in the outer disk the grains are largely left unpro- cessed (van Boekel et al. 2004). In the case of DG Tau, however, we see that the highly processed dust lies in the outer regions of fective temperature, the mid-infrared emitting disk (outside 1 au). In the following we attempt to explain the spatially distinct silicate emission and ˙ 4 3GM M g 1/2 absorption features. T (r) = ∗ [1 R /r ], (14) 8πσr3 in − where G is the gravitational constant, σ is the Stefan-Boltzmann 4.2. Radiative transfer modeling constant, M is the stellar mass, M˙ is the accretion rate and Rg is ∗ in To evaluate the different scenarios we use the radiative transfer the inner radius of the gaseous accretion disk. For the dusty disk code RADMC-3D3 (Dullemond 2012). We assume a flared, ax- we assume a power-law radial distribution for its surface density isymmetric circumstellar disk geometry, where the inner radius Σ, is set at the sublimation radius. DG Tau has a considerableaccre- pΣ r d d tion rate, thus we include a gaseousaccretion disk inside the dust Σ (r) = Σin , R < r < R , (15) Rd in out sublimation radius. We also take into consideration the accretion in heating in the mid-plane of the dusty disk, thus it is heated by d both the accretion and the illumination from the central engine. where Σin is the surface density at the inner radius Rin, pΣ is the d For computing the spectral energy distribution (SED) of the in- power-law exponent, and Rout is the outer radius of the dusty ner gaseous disk, we adopt the Shakura-Sunyaev disk descrip- disk. We assume the vertical distribution to be a Gaussian func- tion (Shakura & Sunyaev 1973), where a geometrically thin, op- tion, and the resulting density structure of the disk is tically thick disk emits black-body radiation with a surface ef- 1 z2 3 ρ (r, z) = Σ (r) exp , (16) http://www.ita.uni-heidelberg.de/~dullemond/software/radmc-3d/ √2πH −2H2 ! Article number, page 11 of 16 A&A proofs: manuscript no. main where ρ(r, z) denotes the dust density as a function of r and the -8 height z above the mid-plane, and H is the scale height. We as- 10 sume that the dependence of H on r is also a power-law, 10-9 ph ) 1
H (r) d r − s h (r) = h , (17) in d 2 ≡ r R − -10 in 10 cm d where h (r) is the dimensionless scale height, h is the dimen- erg in ( -11 10 Modeled SED sionless scale height at the inner radius, and ph is the power-law