EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 585 (January 2016) A&A, 585 (2016) A100 Full HTML
Free Access
Issue
A&A
Volume 585, January 2016
Article Number A100
Number of page(s) 7
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201526691
Published online 23 December 2015

© ESO, 2015

1. Introduction

Young stellar objects (YSOs) such as T Tauri stars gather material from their surrounding disk via accretion. It is commonly believed that the observed brightness variability of YSOs are due to a temporal increase of the stellar accretion rate that only lasts a limited amount of time (Hartmann & Kenyon 1996; Mosoni et al. 2013). Various explanations for the temporal variations of accretion have been proposed: planet-disk interaction (Lodato & Clarke 2004; Ruge et al. 2014), thermal or gravitational instability (Pringle 1981; Hartmann & Kenyon 1996; Lodato & Rice 2004, 2005; Zhu et al. 2009), interactions between the disk and a central binary (Pfalzner 2008), and the impact of the magnetic field (D’Angelo & Spruit 2010). The flux variations can be traced from ultraviolet (UV) up to mid-infrared (MIR) wavelengths has been the subject of many photometric and spectroscopic studies in the past (Lorenzetti et al. 2009; Semkov et al. 2013; Antoniucci et al. 2014; Banzatti et al. 2014). With the rise of optical/infrared long baseline interferometry, one is potentially able to put stringent constraints on the region of the origin of the flux variations, even if not to spatially resolve it.

In this study, we investigate the variability of the protoplanetary disk of DR Tau, a classical T Tauri star. DR Tau is a member of the Taurus-Auriga molecular cloud at a distance of about 140 pc (Kenyon et al. 1994). It is known for its very strong, long- and short-term photometric and spectroscopic variability from the UV to the MIR (Kenyon & Hartmann 1995; Alencar et al. 2001). Therefore the obtained parameters of the surrounding disk and the star itself vary with every new study. For instance, Banzatti et al. (2014) showed that the accretion luminosity can change by a factor of about two to three within less than two months by analysing the UV excess. Besides, age and spectral type of DR Tau could only be limited to be less than 4 Myr and between K5 and K7, respectively (Muzerolle et al. 2003; Greaves 2004; Mohanty et al. 2005; Pontoppidan et al. 2011).

We present self-consistent disk models, simulated with the radiative transfer code MC3D (Wolf et al. 1999; Wolf 2003). These models reproduce highly spatially resolved interferometric data obtained with the MID-infrared Interferometric instrument (MIDI) at the VLTI at three different epochs. The measured visibilities show a clear change in the brightness distribution of the disk and, thus, we are able to determine the underlying variation of the dust density distribution. To reduce the degeneracies, a set of published photometric data is also used to model the spectral energy distribution (SED).

2. Mid-infrared interferometric observations and data reduction

The circumstellar disk of DR Tau was observed with MIDI (Leinert et al. 2003a,b; Morel et al. 2004) at the VLTI at three epochs. For all measurements the instrument was fed by two 8.2 m Unit Telescopes and the prism was used to disperse the light with a low resolution of about 30.

The data were reduced with the MIA+EWS package1. The visibilities shown were computed with MIA that is based on the analysis of the power spectrum. The visibility shown for the first night is the mean of the two calibrated measurements. The standard error is 0.07 on average. All visibilities have been calibrated by all the calibrators taken in the same night and with the same mode as the science target. The stability of the transfer function is reflected by the errors shown together with the visibilities. Calibrators with peculiar visibilities have been ignored and are indicated in the journal of observations; see Table 1. A comparison of the MIA results with the results obtained with EWS confirm the visibilities. For the last night the EWS visibility even tends to be lower up to about 0.1 at short N-band wavelengths. The EWS is based on a coherent analysis of the interferometric fringe signal (Jaffe 2004). The differential phases of our observations do not significantly differ from 0°. The maximum differential phase is at about 10° for E2 but strongly depends on the used calibrator and thus the uncertainties are very high. Besides, differential phases equal to zero do not necessarily mean that there are no asymmetries in the target because linear phase terms cannot be probed (Ratzka et al. 2009).

The data used for our analysis were taken in January 2005, October 2013, and December 2013. Additional measurements were made in October and November 2004 and in December 2013, but were excluded due to their low quality. The data set presented and modelled by Schegerer et al. (2009) has been re-reduced. The calibrated visibilities for the three epochs are shown in Fig. 1.

Table 1

Journal of our VLTI/MIDI observations.

thumbnail Fig. 1

Calibrated MIR visibility of DR Tau for the three different epochs January 2005 (E1), October 2013 (E2), and December 2013 (E3).
Open with DEXTER

3. Data analysis and disk modelling

The properties of our target DR Tau are listed in Table 2. Previous photometric measurements are given in Table 3.

Table 2

Stellar properties of DR Tau.

Table 3

Photometric data of DR Tau.

3.1. Temporal variability

Interferometry: the visibilities obtained at the first epoch in January 2005 (E12) differ significantly from those made at E3 in December 2013, although the projected baseline lengths (BL) and orientations are nearly the same. These observations indicate changes of the brightness distribution and, hence, the underlying density distribution in the innermost regions of the disk (≲4 AU). Thus, these two observations cannot be reproduced with a static disk model. Possible reasons could be a general brightness increase of the disk due to a variable accretion rate onto the central star or azimuthal brightness asymmetries that have significantly changed their position between both epochs of observation. Local brightness asymmetries on the scale of several dozens to hundreds of AU have already been observed in the case of various protoplanetary disks, e.g. GM Aur (Schneider et al. 2003), FS Tau (Hioki et al. 2011), IM Lupi (Pinte et al. 2008). However, in contrast to these observations only long-baseline interferometry in the optical/infrared provides the angular resolution so far to trace regions, so close to the central star that the orbital motion on the scale of months to years has to be taken into account in the data analysis. In the very near future, ALMA will provide a similar angular resolution in the submm/mm wavelength range.

On the other hand, the visibility obtained at E2 (October 2013) with the much longer BL of about 90 m is higher than the measurement with BL ≈ 60 m two months later (E3). A higher visibility means more compact structure in the case of simple structures. Thus, an observation with a longer BL, i.e. tracing higher spatial frequencies, result in a lower visibility, which is not the case here. Of course, we observed at two different baseline orientations, which at first glance might explain this counterintuitive scenario. A static, rotationally symmetric, inclined disk, observed at slightly different orientations can produce different visibilities if one observation is almost aligned with the minor or major axis of the disk. This difference becomes larger with higher orientation difference and higher disk inclination, however, our observations were made at an orientation difference of only 35° and there is no evidence of a disk inclination larger than 60° (Kenyon et al. 1994; Alencar et al. 2001; Brown et al. 2013; Eisner et al. 2014); see also Sect. 4. Thus, it could either be a hint of short-term variations at the same timescale as the UV variations observed by Banzatti et al. (2014) or for a non-rotationally symmetric disk, i.e. disk structures seen at different phases of their orbit. In summary, our interferometric observations indicate the presence of a non-axisymmetric structure of the disk and/or temporal variations of the disk structure in mid-infrared bright regions of the disk.

Photometry: for a moderately inclined disk, we expect to observe a silicate emission feature at ≈10 μm and a smooth SED without depressions or peaks. In contrast to the adjacent photometric data points, the Spitzer/IRAC data suggest a lack of flux at 8 μm; see Table 3. Moreover, this is the same wavelength region MIDI is sensitive to, and the non-stable visibility suggests a non-stable flux. Thus, it is possible that this photometric Spitzer observation was made during a more quiet phase of DR Tau. Because of this, we exclude the IRAC data point from our fitting procedure.

3.2. Model set-up

Earlier studies: many different values have been published for the orientation of the disk, described by the disk inclination i and the disk position angle PA. The inclination ranges from 9° (Pontoppidan et al. 2011) to 60° (Akeson et al. 2005) and the position angle is presumed to be between 75° (Eisner et al. 2014) and 180° (Andrews & Williams 2005). The disk is rather massive, Mdisk ≈ 0.01–0.1 M⊙ (Akeson et al. 2005; Robitaille et al. 2007; Schegerer et al. 2009; Ricci et al. 2010), and the inner radius has been determined to be less than 0.1 AU through analysing the emission of CO and H2O (Brown et al. 2013; Banzatti et al. 2014; Eisner et al. 2014).

Disk structure: Schegerer et al. (2009) also fit photometric and interferometric data for DR Tau but only for one epoch (E1). Their results are used as a starting point for our simulation, which also applies to the Ansatz of Shakura & Sunyaev (1973) for the density distribution. It is given by(1)where r and z are the usual cylindrical coordinates, R⋆ the stellar radius, and h(r) the scale height:(2)The two exponents α and β describe the radial density profile and disk flaring, respectively. Assuming a coupling of temperature and surface density, the number of free parameters decreases by one (Shakura & Sunyaev 1973):(3)This simple disk model has already been used successfully for fitting spatially resolved and unresolved multi-wavelength observations of several protoplanetary disks (Wolf et al. 2003; Ratzka et al. 2009; Sauter & Wolf 2011; Madlener et al. 2012; Gräfe et al. 2013; Liu et al. 2013).

Dust properties: the dust in our set-up is a mixture of 62.5% astronomical silicate and 37.5% graphite (12.5% parallel and 25% perpendicular) and has particle sizes between 5 nm and 250 nm. The grain size distribution is given by the MRN distribution (Mathis et al. 1977), i.e. follows a power law, n(s) ∝ s-3.5, where n·ds are the number of particles in the radius interval [s,s + ds].

Heating sources: the central star is treated as a black body with a fixed value for the total luminosity of Ltot = L⋆ + Lacc = 1.9 L⊙ (Muzerolle et al. 2003; Schegerer et al. 2009; Eisner et al. 2014). Further, we set the stellar effective temperature to T⋆ = 4050 K (Mohanty et al. 2005). The disk around DR Tau is actively accreting (Robitaille et al. 2007; Edwards et al. 2013; Eisner et al. 2014). Therefore, we also assume an additional accretion heating of the disk. To implement this, we use a very simple, static accretion model, by adding a second point-like, black-body source in the centre of our disk. The properties of this accretion source, i.e. temperature and truncation radius, are fixed in this study to Taccr = 8000 K and Rtrunc = 5 R⋆, respectively. These values were already successfully used to model other TTSs (Akeson et al. 2005; Schegerer et al. 2008, 2009).

To calculate the dust temperature distribution and the corresponding SED and mid-infrared visibilities, we use the radiative-transfer code MC3D (Wolf et al. 1999; Wolf 2003).

3.3. Modelling

Fitting procedure: in the standard fitting approach, the photometric and interferometric data were used simultaneously. The SED is a degenerate quantity and more than just one adequate density distribution could reproduce the observations. The visibility curves are then used to reduce these degeneracies, since they are only sensitive to small-scale structures in the interior of the disk. According to Sect. 3.1, this approach is not suitable for our observations. The visibilities obtained at the three epochs cannot be reproduced by a disk with only one parameter set. Thus, we independently fit all four sets, SED + three visibilities. Due to the set-up of our global model, a resulting parameter set based on a visibility is only valid for the innermost brightness distribution and does not represent the entire disk.

Parameter space and fitness quality: to explore the parameter space, we use an iteratively refined grid search. The boundaries of the parameter space and the overall number of used values for each parameter are shown in Table 4. To compare all individual models, we calculate χ2 for each resulting SED (4)where m is the number of wavelengths in Table 3 and σSED the photometric uncertainties. To evaluate the visibilities, we apply a similar approach. Because of the huge amount of time needed to calculate a radiation map with sufficient resolution, we do not calculate the visibilities and thus the χ2 value for every measured wavelength bin given by the spectral resolution of MIDI, but at seven wavelengths between 8 and 13 μm.

thumbnail Fig. 2

Best SED fit with corresponding MIR visibilities. The three visibility epochs are shown in different colours (black: first, red: second, blue: third). Observations: with uncertainties; model: solid lines.
Open with DEXTER

Table 4

Boundaries of parameter space.

Table 5

Model parameters for our best fits.

4. Results

In the upper part of Fig. 2, our best SED fit is shown. We excluded the Spitzer/IRAC data point at 8 μm from the fitting process (see Sect. 3.1). The properties of this disk are listed in the second column of Table 5. These values are in agreement with other DR Tau studies; cf. Banzatti et al. (2014), Eisner et al. (2014), Isella et al. (2009), Schegerer et al. (2009), Muzerolle et al. (2003). These parameters, however, cover a rather wide range, which could be directly related to the short- and long-term variations of the disk structure. Interestingly, the best-fit parameters result in mid-infrared visibilities, which are in agreement with those obtained in E3, although the visibilities were not considered during the SED fitting process (see lower part of Fig. 2).

The corresponding visibilities in Fig. 2 are plotted for a PA ≈ 50° from north to east for the major axis of the disk. As predicted in Sect. 3.1, the two visibilities with nearly the same set-up from January 2005 (E1) and December 2013 (E3) cannot be reproduced with the same disk parameter values. Although the PA is a free parameter here, the baseline orientation difference of just 10° is too small to produce this difference in the visibility. Because of the small inclination of the disk, the visibility is not very sensitive to changes in the PA. The calculated visibilities with BL = 60.6 m (blue solid line) vary only between 5% and 12% due to changes in the PA. This uncertainty is also confirmed by previous DR Tau studies: Pontoppidan et al. (2011) found a PA of about 0°, Isella et al. (2009) found 100°, Brown et al. (2013) published 140° and Akeson et al. (2005) derived PA = 160° ± 55°.

In Fig. 3 our best result for the visibility obtained at epoch E3 is shown. As you can see, both the visibility curve for E3 and the SED are very similar to those in Fig. 2. The major differences in both set-ups are only in the inner radius, which changed from 0.111 AU to 0.065 AU, and the outer radius, which changed from 200 AU to 90 AU. The larger outer radius of the first model results in more cold dust, which increases the flux in the far-IR and (sub-)mm regime and matches the photometry better at these wavelengths. Furthermore, MIDI measurements are not sensitive to the outer disk because of the limited FoV of the interferometer. The other parameters are the same or very similar for both models, especially inclination and flaring, and thus the radial density profile is almost unchanged (see Eq. (3)).

A reasonable scenario for the found differences between the visibilities of E1 and E3 could be fast changes in the accretion luminosity, as mentioned in Sect. 1. To investigate this, we used the disk set-up that reproduces the SED best, and set Laccr = 1.9 L⊙, resulting in a total luminosity of Ltot = 2.8 L⊙. As a consequence, the dust heating is stronger at all distances and disk layers. This leads to a SED that is higher at all wavelengths and to visibilities that drop by about 10% for all three baselines. If the accretion luminosity is decreased or even completely eliminated (Ltot = L⋆ = 0.9 L⊙), the visibilities increase up to 20% for the longest wavelengths. Therefore, the lack of accretion cannot reproduce our observations as the visibilities are still too low for E1 and E2, and hence a different disk set-up is needed.

thumbnail Fig. 3

Best fit for third epoch visibility with corresponding SED. The three visibility epochs are shown in different colours (black: first, red: second, blue: third). Observations: with uncertainties; model: solid lines.
Open with DEXTER

Our third parameter set reproduces the observations of E1 and E2 simultaneously, although they are almost nine years apart (see Fig. 4). If we compare the parameters of this best-fit model with those of the two previous models, we find that only the scale height and outer radius have changed and that flaring and inclination are again the same. This disk set-up is only valid for small stellar separations of r ≲ 10 AU. More distant regions are not covered by the MIR visibilities and hence not considered by the fitting process. Because of the huge outer radius of model 3, the density in the hot, inner region is lower. In contrast, the scale height is decreased by nearly 50%, which significantly increases the density in the midplane of the disk, resulting in a higher optical depth and a more compact emission. Hence, the visibility in the MIR is higher compared to the other two models and fits the observations. Owing to our purely global disk model, the SED in Fig. 4 does not reproduce the observations and is too low for nearly all wavelengths. With the much smaller scale height, the disk of model 3 is very compact, resulting in a higher optical depth. Thus, the dust heating is less efficient, resulting in a reduced flux.

thumbnail Fig. 4

Best fit for first and second epoch visibility with corresponding SED. The three visibility epochs are shown in different colours (black: first, red: second, blue: third). Observations: with uncertainties; model: solid lines.
Open with DEXTER

5. Conclusion

We presented MIR visibilities of the T Tauri star DR Tau obtained with MIDI/VLTI. We modelled the SED and the individual mid-infrared visibility curves at three epochs to constrain the underlying variation in the disk density structure. Our simulations suggest that the hot, inner regions of the disk were more compact in January 2005 (E1) than they were in December 2013 (E3), especially in the layers close to the midplane. Since these two observations were made with nearly the same baseline length and position angle, it is possible that the structural difference is due to a local density variation, such as a density clump observed at different stages of its orbit. If this is true, the relevant region can be approximated by assuming a Keplerian orbit. If this hypothetical clump revolved once or twice within these nine years, its orbit must be about 2.5 or 4 AU, respectively. For a stellar distance of 140 pc, this is exactly the resolution of MIDI with a BL ≈ 60 m. Since we cannot distinguish between PA = 0° and PA = 180°, we detect the same fringe pattern twice per orbit. Although it is very unlikely, it cannot be completely ruled out that we observed the same density structure at first and second epoch. Considering a circular orbit of 2.5 AU with a period of ≈4.5 years, in two months the clump would have moved ≈13°−14° on its orbit. Together with the PA difference of 35° between E2 and E3, it is also possible to explain the non-detection at E3.

However, we are aware of the fact that the photometric data set (see Table 3) is based on observations made at different epochs, i.e. at different stages of temporal evolution of the disk. Detailed studies over different timescales are necessary to achieve a better understanding of the temporal variations of protoplanetary disks and thus the underlying disk physics. With this pioneering study, we want to stress this problem, which has been commonly ignored in most disk modelling studies so far.

With the second generation instrument MATISSE on the VLTI, a much better uv-coverage will be possible than with MIDI (Lopez et al. 2014). The MATISSE instrument will measure up to six visibilities at once with four telescopes, which is crucial for a consistent view on the disk at one certain point in time. Moreover, studying the temporal evolution of substructures could constrain physical processes in the potential planet-forming region (Kley et al. 2009; Flock et al. 2015; Bitsch et al. 2015). With existing MIDI data, new observations will enable us to trace disk evolution on timescales of more than 15 years. Even though MATISSE will cover a much wider wavelength range than MIDI (L, M, and N band), it is mandatory to observe protoplanetary disks with other interferometric instruments in different wavelength regimes also. High-spatial-resolution observations in the NIR and (sub-)mm, such as obtained with PIONIER/VLTI, AMBER/VLTI, and ALMA, will allow us to constrain the density distribution for dust at different temperatures, i.e. different disk regions.

2

In the following, we denote the epochs 1, 2, and 3 as E1, E2, and E3, respectively.

Acknowledgments

This research was funded through the DFG grant: WO 857/13-1. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.

References

All Tables

Table 1

Journal of our VLTI/MIDI observations.

In the text
Table 2

Stellar properties of DR Tau.

In the text
Table 3

Photometric data of DR Tau.

In the text
Table 4

Boundaries of parameter space.

In the text
Table 5

Model parameters for our best fits.

In the text

All Figures

thumbnail Fig. 1

Calibrated MIR visibility of DR Tau for the three different epochs January 2005 (E1), October 2013 (E2), and December 2013 (E3).
Open with DEXTER
In the text
thumbnail Fig. 2

Best SED fit with corresponding MIR visibilities. The three visibility epochs are shown in different colours (black: first, red: second, blue: third). Observations: with uncertainties; model: solid lines.
Open with DEXTER
In the text
thumbnail Fig. 3

Best fit for third epoch visibility with corresponding SED. The three visibility epochs are shown in different colours (black: first, red: second, blue: third). Observations: with uncertainties; model: solid lines.
Open with DEXTER
In the text
thumbnail Fig. 4

Best fit for first and second epoch visibility with corresponding SED. The three visibility epochs are shown in different colours (black: first, red: second, blue: third). Observations: with uncertainties; model: solid lines.
Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.