Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.gao.spb.ru/english/personal/sfg/eng/articles/2011/1101.3695v1.pdf Äàòà èçìåíåíèÿ: Fri Jan 27 15:45:01 2012 Äàòà èíäåêñèðîâàíèÿ: Tue Oct 2 03:39:36 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: m 80

Astronomy & Astrophysics manuscript no. 15676 January 20, 2011

c ESO 2011

VLTI/AMBER spectro-interferometry of the Herbig Be star MWC 297 with spectral resolution 12 000
G. Weigelt1 , V. P. Grinin2,3 , J. H. Groh1 , K.-H. Hofmann1 , S. Kraus4 , A. S. Miroshnichenko5 , D. Schertl1 , L. V. Tamb ovtseva2 , M. Benisty6 , T. Drieb e7 , S. Lagarde8 , F. Malb et9 , A. Meilland1,8 , R. Petrov8 , and E. Tatulli9
1 2

arXiv:1101.3695v1 [astro-ph.SR] 19 Jan 2011

3 4 5 6 7 8 9

Max-Planck-Institut fur Radioastronomie, Auf dem Hugel 69, D-53121 Bonn, Germany ¨ ¨ Pulkovo Observatory of RAS, Pulkovskoe shosse 65, St. Petersburg 196140, Russia V.V. Sob olev Astronomical Institute, St. Petersburg University, St. Petersburg, Russia Department of Astronomy, University of Michigan, 500 Church Street, Ann Arb or, MI 48109, USA University of North Carolina at Greensb oro, P.O. Box 26170, Greensb oro, NC 27402, USA INAF-Osservatorio Astrofisico di Arcetri, Istituto Nazionale di Astrofisica, Largo E. Fermi 5, I-50125 Firenze, Italy Deutsches Zentrum fur Luft- und Raumfahrt e.V., K¨nigswinterer Str. 522-524, D-53227 Bonn, Germany ¨ o Lab oratoire H. Fizeau, University of Nice Sophia Antip olis, CNRS, OCA, Parc Valrose, 06108 Nice, France Lab oratoire d'Astrophysique de Grenoble, UMR 5571 Universit´ Joseph Fourier/CNRS, BP 53, F-38041 Grenoble e Cedex 9, France

Received / Accepted ABSTRACT Context. Circumstellar disks and outflows play a fundamental role in star formation. Infrared sp ectro-interferometry allows the inner accretion-ejection region to b e resolved. Aims. We study the disk and Br -emitting region of MWC 297 with high spatial and sp ectral resolution and compare our observations with disk-wind models. Methods. We measured interferometric visibilities, wavelength-differential phases, and closure phases of MWC 297 with a sp ectral resolution of 12 000. To interpret our MWC 297 observations, we employed disk-wind models. Results. The measured continuum visibilities confirm previous results that the continuum-emitting region of MWC 297 is remarkably compact. We derive a continuum ring-fit radius of 2.2 mas (0.56 AU at a distance of 250 p c), which is 5.4 times smaller than the 3 AU dust sublimation radius exp ected for silicate grains (in the absence of radiation-shielding material). The strongly wavelength-dep endent and asymmetric Br -emitting region is more extended (2.7 times) than the continuum-emitting region. At the center of the Br line, we derive a Gaussian fit radius of 6.3 mas HWHM (1.6 AU). To interpret the observations, we employ a magneto-centrifugally driven disk-wind model consisting of an accretion disk, which emits the observed continuum radiation, and a disk wind, which emits the Br line. The calculated wavelength-dep endent model intensity distributions and Br line profiles are compared with the observations (i.e., K -band sp ectrum, visibilities, differential phases, and closure phases). The closest fitting model predicts a continuumemitting disk with an inner radius of 0.3 AU and a disk wind ejection region with an inner radius of 0.5 AU (17.5 stellar radii). We obtain a disk-wind half-op ening angle (the angle b etween the rotation axis and the innermost streamline of the disk wind) of 80 , which is larger than in T Tau models, and a disk inclination angle of 20 (i.e., almost p ole-on). Conclusions. Our observations with a sp ectral resolution of 12 000 allow us to study the AU-scale environment of MWC 297 in 10 different sp ectral channels across the Br emission line. We show that the K -band flux, visibilities, and remarkably strong phases can b e explained by the employed magneto-centrifugally driven disk wind model. Key words. Stars: individual: MWC 297, Stars: pre-main sequence, Stars: winds, outflows, Stars: circumstellar matter, Techniques: interferometric, Techniques: sp ectroscopic

1. Intro duction
Circumstellar disks and outflows play a fundamental role in star formation. Testing theoretical accretion-ejection models is very challenging since important processes occur on AU and sub-AU scales, which cannot be resolved by conventional imaging. However, infrared spectro-interferometric observations with milli-arcsecond spatial resolution and
Send offprint requests to : G. Weigelt email: weigelt@mpifr.de Based on observations made with ESO telescop es at Paranal Observatory under programme ID 081.D-0230(A).

high spectral resolution allow unprecedented studies of the inner accretion region, which hosts fascinating astrophysical ob jects such as inner gaseous accretion disks, dust disks, disk winds, jets, etc. Interferometry in the near- and midinfrared is able to probe the gas and dust distribution on these scales and constrain their physical properties. These observations allow us to study the nature of the accretion process and the launching of jets and winds. With a spectral type of B1.5V, a mass of 10 M , and a distance of 250 ± 50 pc (Drew et al. 1997), MWC 297 is one of the nearest massive young stars. The geometry of the disk of MWC 297 has been investigated by various techniques
1


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

at near-infrared (Millan-Gabet et al. 2001; Eisner et al. 2004a; Monnier et al. 2006; Malbet et al. 2007; Acke et al. 2008; Kraus et al. 2008a), mid-infrared (Acke et al. 2008), and radio wavelengths (Mannings & Sargent 1997; Drew et al. 1997). Interestingly, several interferometric observations of MWC 297 (Eisner et al. 2004b; Monnier et al. 2005; Malbet et al. 2007; Acke et al. 2008; Kraus et al. 2008a) have shown that the compact characteristic size of the near-infrared continuum-emitting region is several times smaller than the theoretically expected 3 AU radius of the dust sublimation rim (the expected radius if there is no shielding of the stellar radiation; Monnier et al. 2005). The estimates of the inclination angle of MWC 297's circumstellar disk have been highly controversial. Drew et al. (1997) argued for an almost edge-on disk orientation. Moderate inclination angles of the circumstellar disk of MWC 297 were found with infrared interferometric observations. For instance, Malbet et al. (2007) obtained VLTI/AMBER observations with a medium spectral resolution of 1500 and measured the wavelength dependence of the visibility across the Br line. Malbet et al. (2007) interpreted their spectro-interferometric measurements of MWC 297 with an anisotropic stellar-wind model (Stee & de Araujo 1994; Stee et al. 1995) and found an inclination of 20. Acke et al. (2008) obtained both nearinfrared AMBER H - and K -band observations and midinfrared MIDI observations and concluded that the inclination is smaller than 40 . In this paper, we present the first spectrointerferometric measurements of the inner circumstellar environment of MWC 297 with AMBER's high spectral resolution mode (R = 12 000). This high spectral resolution allows us to spectrally resolve the Br emission line and measure the interferometric observables (spectrum, visibilities, wavelength-differential phases, and closure phases) in 10 different spectral channels across the line. In Sect. 2, we present the AMBER observations. In Sect. 3, we derive characteristic continuum and Br line sizes of the circumstellar environment using simple geometric models. In Sect. 4, we present radiative transfer modeling of the disk and disk wind and compare the K -band flux, visibilities, and phases with the model predictions in the continuum and several spectral channels across the Br emission line. In the Appendix, we discuss the disk-wind modeling, spectroscopy of MWC 297, wavelength calibration, and the use of AMBER's Beam Commutation Device.

velocity (LSR) [km/s] -100 0 100 5.0 4.0 3.0 Flux 2.0 1.0 1.0 0.9 0.8 0.7 Visibility 0.6 0.5 0.4 0.3 0.2 0.1 0.0 160 120 Diff. Phase [deg] 80 40 0 -40 -80 -120 -160 160 Closure Phase [deg] 120 80 40 0 -40 -80 -120 -160 PBL 14.0m, PA 68o PBL 28.0m, PA 68 o PBL 42.1m, PA 68
o

2. Interferometric VLTI/AMBER/FINITO observations with sp ectral resolution of 12 000
We observed MWC 297 on 6 April 2008 with ESO's Very Large Telescope Interferometer (VLTI) and its AMBER beam combiner instrument in the course of the AMBERGTO program 081.D-0230(A). AMBER (Petrov et al. 2007) allows interferometric three-telescope observations with low, medium, or high spectral resolution. The observations are described in Table 1. For these observations, the E0-G0-H0 array of the 1.8 m auxiliary telescopes (ATs) and AMBER's high resolution mode (HR mode; R = 12 000) were employed. To obtain interferograms with a high signalto-noise ratio (SNR), we used the fringe tracker FINITO
2

2.163 2.164 2.165 2.166 2.167 2.168 2.169 Wavelength [µm]

Fig. 1. AMBER observations of MWC 297 with spectral resolution of 12 000: (from top to bottom) wavelength dependence of flux, visibilities, wavelength-differential phases (for better visibility, the differential phases of the shortest and longest baselines are shifted by +40 and -40, respectively), and closure phase observed at pro jected baselines of 14.0, 28.0, and 42.1 m along the position angle (PA) of 68.0 on the sky. The wavelength scale at the bottom is the directly observed one, i.e. without heliocentric or LSR correction (see text). However, the radial velocity scale at the top, gives spectrum, visibilities, and phases as a function of the Br Doppler shift in the Local Standard of Rest (LSR) frame. The dashed vertical line indicates the centroid Br vacuum wavelength (2.1661 µm) in the LSR frame.


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

and a detector integration time of 8.0 seconds per interferogram. With excellent atmospheric conditions (seeing 0.4­0.6), the FINITO instrument in the co-phasing mode allowed us to stabilize the atmospheric and mechanically induced fringe motion. FINITO records temporally modulated H -band interferograms with a high frequency, from which the atmospheric phase shift is determined and then compensated for using the VLTI delay lines. For data reduction, we used AMBER's amd lib data reduction software (release 2.2), which employs the P2VM algorithm (Tatulli et al. 2007) to derive wavelengthdependent visibilities, wavelength-differential phases, and closure phases. In addition, the spectrum can be extracted as an additional observable. Along with the science observation of MWC 297, we observed the interferometric calibrator star HD 175 583, which allowed us to calibrate the atmospheric transfer function. We recorded 100 target interferograms (plus 40 calibrator interferograms) without AMBER's beam commutation device (BCD OUT data) and 100 target interferograms plus 40 calibrator interferograms with the BCD (BCD IN data). The BCD allows the study of systematic instrumental effects (see Appendix C). A fraction of the interferograms were of lower quality because of the lower FINITO performance. Therefore, we selected the 50% of the BCD OUT data and 50% of the BCD IN data with the highest fringe SNR among both the target and the calibrator interferograms to improve the visibility calibration (Tatulli et al. 2007). From this data, we derived the spectrum, visibilities, wavelength-differential phases, and closure phases, observed along the position angle (PA) of 68.0 on the sky (see Fig. 1). The differential phases and closure phases derived from both independent data sets (BCD OUT and BCD IN) are shown in Fig. C.1 of Appendix C. Subtraction of the closure phases obtained with and without BCD cancels potential instrumental phase drifts. Figure C.2 presents the derived closure phase and a fit of a constant to the continuum closure phases. The obtained averaged continuum (2.147­2.194 µm) closure phase is -1.28 ± 0.42 . With the high spectral resolution of R = 12 000, it is possible to spectrally resolve the Br line and measure the visibilities and phases in 10 different spectral channels across the Doppler-broadened Br line. The visibilities show a strong wavelength dependence of the size of the lineemitting region. The wavelength-differential phases and the closure phase indicate remarkably strong photocenter shifts and asymmetries, respectively, of the circumstellar environment in many spectral channels across the Br line. We show in Sect. 4 that this strong wavelength dependence and the large values of the phases can be explained by disk-wind models. The wavelength calibration of the AMBER data in Fig. 1 was accomplished using the many telluric lines present in the region 2.15­2.19 µm, as described in the Appendix B. We estimate an uncertainty in the wavelength calibration of 3 km s-1 . The spectrum in Fig. 1 (top) was normalized to the nearby continuum. The wavelength scale at the bottom is shown without heliocentric or Local Standard of Rest (LSR) correction. However, the radial velocity scale at the top of Fig. 1 is calibrated to give spectrum, visibilities, and phases as a function of the Br Doppler shift in the Local Standard of Rest (LSR) frame (16.5 km s-1 has to be subtracted to convert the LSR values to the heliocentric frame).

The obtained Br line profile has an equivalent width of -17.2 °. Using AMBER's HR mode, we were able to reA solve the Br line and determined a FWHM of 60 km s-1 . By performing a Gaussian fitting of the line profile, a centroid LSR velocity of 13.8 ± 3 km s-1 was obtained. No P-Cygni absorption or double-peaked emission is seen. The FWHM of the Br line is compatible with the value of 53 km s-1 found by Drew et al. (1997) for Br. On the other hand, the centroid velocity of Br is more similar to the one measured for H and H (16 and 19 km s-1 , respectively; Drew et al. 1997) than to that of Br (3 km s-1 , Drew et al. 1997).

3. Comparison of the observations with geometric mo dels
3.1. Characteristic size of the continuum-emitting region To measure the characteristic size of the circumstellar environment in the K -band continuum, we fitted geometric models to the visibilities obtained along the PA of 68.0 on the sky (see Table 2). As pointed out by Acke et al. (2008), who observed MWC 297 with the same VLTI E0-G0-H0 array, one-component uniform-disk, ring, and Gaussian disk models are unable to fit the shape of the continuum visibility. Therefore, we fitted a two-component Gaussian plus a point source with a 10% flux contribution (derived from the SED; as in Acke et al. 2008) to the wavelength-averaged continuum visibilities. Our observations are averages over the continuum wavelengths between 2.147 and 2.194 µm (except the line region) because of the narrower wavelength range of AMBER's HR mode. Figure 2 (left) shows that the fitted two-component Gaussian consists of a dominant compact Gaussian with a FWHM diameter of 4.63 ± 0.21 mas (corresponding to 1.16 ± 0.05 AU at a distance of 250 pc) in addition to an extended Gaussian halo with a FWHM diameter of 22 mas (see Table 2). Owing to the lack of data at low spatial frequencies, the size and shape of the extended halo component is not well constrained. The flux ratio of the dominant compact component to the extended halo component is 4.6. Our FWHM diameter of 4.6 mas obtained for the dominant compact Gaussian component is in good agreement with the diameter of 4.3 mas reported by Acke et al. (2008). To verify MWC 297's location in the size-luminosity relation (where ring fits are used; Monnier et al. 2005; Kraus et al. 2008a), we determined its ring-fit radius for a ring thickness of 20% of the radius. Figure 2 (right) shows the fit obtained of a compact ring, an extended Gaussian halo, and a 10% point source to the wavelengthaveraged (2.147­2.194 µm) continuum visibilities. The fitted FWHM diameter (inner ring diameter) of the dominant compact model ring is 4.47 ± 0.20 mas (corresponding to 1.12 ± 0.05 AU for a distance of 250 pc; see Table 2). The measured continuum ring radius of 0.56 AU is 5.4 times smaller than the 3 AU dust sublimation radius of MWC 297 predicted by the size-luminosity relation, if there is no shielding of the stellar radiation (Monnier et al. 2005; Kraus et al. 2008a). This surprisingly compact size was already discussed by several other authors (Eisner et al. 2004a; Monnier et al. 2005, 2006; Malbet et al. 2007; Kraus et al. 2008a; Acke et al. 2008). The correct interpretation of the compact continuum3


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

Table 1. Log of the VLTI/AMBER/FINITO observations of MWC 297.
Date Time [UT] Start End 08:14 08:57 AT array Sp ectral mo d e H R - K- F
d

Wavelength Rrange [µm] 2.147­2.194

DITa [s] 8

b Nt

Calibrator

c Nc

Uniform-disk diameter of the calibrator [mas] 1.02e

2008 Apr. 06

E0-G0-H0

200

HD 175 583

80

Notes ­ a Detector integration time p er interferogram; b Numb er of MWC 297 interferograms; c Numb er of calibrator interferograms; d HR mode in the K band using the fringe tracker FINITO; e Uniform-disk diameter taken from M´rand et al. e (2005).

1

MWC297, continuum visibility at 2.160 micron Two-Gaussian fit: FWHM 4.63 mas & 21.95 mas

1

MWC297, continuum visibility at 2.160 micron Ring diam. 4.47 mas & Gaussian FWHM 19.37 mas

0.8

0.8

visibility

0.6

visibility 0 20 40 60 80 100 120 140 160

0.6

0.4

0.4

0.2

0.2

0 cycles/arcsec

0 0 20 40 60 80 100 120 140 160 cycles/arcsec

Fig. 2. Left: Fit of a two-component Gaussian plus a point source with a 10% flux contribution to the wavelengthaveraged (2.147­2.194 µm) K -band continuum visibilities near Br . Right: Fit of a compact ring plus Gaussian halo plus a point source with a 10% flux contribution to the K -band continuum visibilities.
1 14 12 10 dBr [mas] 8 6 4 2 0 2.1654 2.1656 2.1658 [µm] 2.166 2.1662 flux [arbitrary units] 0.8 dBr = 12.61 ± 0.96 mas

5 4 3 2

visibility

0.6

0.4

0.2 1 0 2.1664 0 0 20 40 60 cycles/arcsec 80 100 120

Fig. 3. Left: Gaussian FWHM fit diameters of the Br linethe center and the line wings) across the Br emission line compensated visibilities. Therefore, they are the diameters of of both line plus continuum. The AMBER spectrum of MWC fit to the continuum-compensated AMBER visibilities at the emitted region of MWC 297 and several other Herbig Be stars remains disputed. Possible explanations include the absorption of the stellar light by an inner gas disk allowing dust to survive closer to the star (e.g., Millan-Gabet et al. 2001; Monnier et al. 2005; Malbet et al. 2007), emission from an optically thick inner gas disk (e.g., Monnier et al. 2005; Kraus et al. 2008b), and emission from special types of dust grains that can survive at higher temperatures (e.g., Benisty et al. 2010).

emitting region for several individual spectral channels (in (red bullets). The diameters are derived from continuumthe Br line-emitting region and not of the emission region 297 is shown for comparison (green curve). Right: Gaussian peak intensity of the Br emission line. 3.2. The characteristic size of the Br -emitting region The high spectral resolution of the AMBER data allows us to measure the characteristic size of the Br lineemitting region for several individual spectral channels across the Br emission line (see Fig. 3 left). These diameters were determined by fitting Gaussians to the continuumcompensated Br visibilities (which are the visibilities free from contributions from the underlying continuum within the line region). The right panel of Fig. 3 shows such a Gaussian fit for the central wavelength of the Br emission line.

4


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

Table 2. Geometric model fit radii of MWC 297 in the continuum and the center of the Br emission line.
Wavelength range Continuum Continuum Br center
a a

comp onent C1 (dominant compact core) Gaussian Ring (20%) Gaussian
b

diameter of C1 4.63 ± 0.21 mas 4.47 ± 0.20 mas 12.6 ± 0.75 mas

comp onent C2 (halo) Gaussian Gaussian -

diameter of C2 > 22 mas > 19 mas c c

flux ratio (C1/C2) 4.6 3.6 -

d

Notes ­ a average over continuum wavelengths b etween 2.147 and 2.194 µm (except the line region). b Inner radius of a ring model with a ring width of 20% of the inner ring radius. c Due to the lack of data at small spatial frequencies, the size of the extended halo comp onent is not well constrained by our measurements. Therefore, we can only give a rough estimate: 22(19) ± 5 mas or larger, but smaller than the AMBER AT FOV of 250 mas. d Continuum-comp ensated diameter (see text).

The continuum-compensated visibilities required for the size determination of the line-emitting region discussed above were calculated in the following way. Within the wavelength region of Br , the measured visibility has two constituents, a pure line-emitting component and a continuum-emitting one, the second of which includes emission of both the circumstellar environment and the unresolved star. The emission line visibility VBr can be written as (Weigelt et al. 2007) FBr VBr = = |Ftot Vtot |2 + |Fc Vc |2 - 2 Ftot Vtot Fc Vc · cos, (1) Fig. 4. Sketch of the MWC 297 disk-wind mo del adopted in this paper.

where Vtot (Ftot ) denotes the measured total visibility (flux) in the Br line, Vc (Fc ) is the measured visibility (flux) in the continuum, and describes the measured wavelengthdifferential phase within the Br line. ¿From the AMBER spectrum (see Fig. 1), we find Ftot = 4.5 and FBr = 3.5âFc at the emission line center. We derived a large continuum-compensated radius (HWHM) of 6.3 ± 0.4 mas (1.6 ± 0.1 AU) for the lineemitting region at the peak wavelength of the emission line, which is 3 times larger than the 0.56 AU HWHM radius of the compact continuum-emitting Gaussian component. The smaller line-continuum size ratio obtained by Kraus et al. (2008a) can be explained by the lower spectral resolution of R = 1500, which leads to an averaging over different diameters corresponding to different wavelengths. Figure 3 also shows that the size of the line-emitting region is smaller in the wings of the Br emission line than at the line center. This wavelength dependence is one of several observational results that have to be explained by models, such as the disk-wind model presented in the next section.

topology of the magnetic field in the disk: magnetic field lines in the wind-launching region are inclined with respect to the disk plane by less than 60 (Blandford & Payne 1982). The radiation pressure of hot stars probably changes the tra jectories of the gas streams, leading to a flatter disk wind (de Kool & Begelman 1995; Drew et al. 1998; Everett et al. 2001). Figure 4 shows a sketch of the disk-wind model that we consider. It assumes that both an accretion disk and a disk wind contribute to the observed emission. 4.1. Geometry and kinematics of the disk wind To model the extended disk wind, we follow an approach used by Shlosman & Vitello (1993) in their studies of disk winds in cataclysmic variables. This approach provides a simple parametrization of a disk wind, which has geometrical and kinematical properties similar to those used by Kurosawa et al. (2006) for T Tauri stars (see Fig. A.1 in the Appendix). Here, we use a similar description of the kinematical parameters of the wind where according to the conservation of angular momentum, the tangential velocity u decreases along a streamline from the Keplerian velocity at the base of streamline. The radial velocity v increases along the streamline from the initial value v0 to v (see Eqs. A.1 and A.2 and Fig. A.2 in Appendix A), where v0 and v are the model parameters. The local mass-loss rate per unit area of the disk, mw , is a function of the cylindric radius . To describe mw , we use a simple power-law mw ( )
-

4. Comparison of the observations with a mo del of a magneto-centrifugally driven disk wind
We now present a disk-wind model and compare the observed spectrum, visibilities, and phases of MWC 297 with the model predictions. In this disk-wind model, mass ejection is driven by magneto-centrifugal forces. The theory of these winds was originally proposed by Blandford & Payne (1982) and has been elaborated by many authors, mainly for T Tauri stars (see K¨nigl & Pudritz 2000, o Pudritz & Banerjee 2005, Ferreira 2007, and references therein). An important property of the launching mechanism in magneto-centrifugally driven disk winds is the

,

(2)
5


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

where is a free parameter. The other model parameter is the total mass-loss rate Mw = 2

1 N

mw ( ) 2 d ,

(3)

where 1 and N are the inner and outer boundaries of the disk-wind ejection region, and the factor two accounts for the mass loss from both sides of the disk surface. 4.2. The parameters of the accretion disk In classical accretion-disk models, the surface density increases toward the star as ( )
-p

,

(4)

where p = 1­1.5 (see, e.g., Calvet et al. 2000). At a mass accretion rate of Macc 10-7 M yr-1 , this disk is optically thick (due to the gas opacity) between the dust evaporation radius and either the stellar magnetosphere (Muzerolle et al. 2004) or the stellar surface. For MWC 297, the value of Macc is fairly uncertain (see Malbet et al. 2007). If the surface density in the accretion disk of MWC 297 changes as -1 , then this disk will be optically thick up to the stellar surface. However, if the accretion disk extended to the innermost vicinity of the star, the disk luminosity at near-infrared wavelengths would exceed the observed one. This is the main reason why an inner disk gap is required. As in Acke et al. (2008), we used a blackbody approximation with a power-law distribution of the disk temperature Td Td (R) = Td (Rin )(R/Rin ) , (5)

where Rin is the inner radius of the disk, and is a parameter. The parameter is equal to -3/4 in the two important special cases when a flat disk is just heated by the stellar radiation (passive disk; Adams & Shu 1986) or radiates due to the viscous dissipation (Lynden-Bell & Pringle 1974). However, this disk (as well as other simple power-law disk models) fails to reproduce the shape of the observed continuum visibilities (Acke et al. 2008). More sophisticated models can be obtained with two-component temperaturegradient models with slightly different temperature gradients (see below). 4.3. Model calculations For simplicity, it is assumed that the disk wind contains only hydrogen atoms. For the calculations of the ionization state and the number densities of the atomic levels, we adopted the numerical codes developed by Grinin & Mitskevich (1990) and Tambovtseva et al. (2001) for moving media. These codes are based on the Sobolev (1960) approximation in combination with the exact integration of the line intensities (see Appendix A). This method takes into account the radiative coupling in the local environment of each point caused by multiple scattering. The stellar radiation of MWC 297 is described by a model atmosphere (Kurucz 1979) with an effective temperature Teff = 25 000 K, log g = 4, and the stellar radius R =
6

6.1 R as in Drew et al. (1997). The calculations were performed for an isothermal wind. We adopted electron temperature Te values in the range 8000­10 000 K, which are typical of the disk winds of young stars (see, e.g., Kurosawa et al. 2006). In our models, we adopted the fixed parameters stellar mass M = 10 M and a distance of 250 pc (Drew et al. 1997). Figure 5 presents examples of the intensity distributions (pro jected on the plane of the sky) of our best-fit continuum-disk and disk-wind model number 5 (see Tables A.1 and 3). The inclination angle of the model (i.e., angle between the polar axis and the viewing direction) is i = 20 , i.e., almost pole-on. Model images at 15 different wavelengths (Doppler shifts) across the Br line profile are shown. The images are highly asymmetric and wavelength-dependent. Their wavelength-dependent photocenter shift (with respect to the continuum image) and asymmetric shape change across the line, resulting in wavelength-dependent visibilities, differential phases, and closure phases. This asymmetry of the model images is the result of the Doppler shift of the radiation due to the gas motion in the emitting region. The main part of the Br line emission is formed near the base of the disk wind where the dominant motion is Keplerian rotation. For the calculation of the model images shown in Fig. 5, a clockwise motion of the disk wind was assumed. When comparing with the observations, we calculated the model images for both clockwise and anti-clockwise motion. Examples of theoretical Br line profiles of model 5 for different inclination angles i between the polar disk axis and the viewing direction (10, 20, 40, and 60 ) are shown in Fig. 6. They suggest that a small inclination of the circumstellar disk (20 ± 10) is needed to explain the observed line profiles (Fig. 1). In our model calculations, we calculated more than 100 different models to examine the parameter space described in Table 3. Five examples of these models are listed in Table A.1. Closest agreement is obtained for a two-component disk model (with two temperature-gradient regimes) consisting of two zones from Rin = 0.3 AU to RS = 0.9 AU and from RS = 0.9 AU to Rout = 3 AU separated at Rs = 0.9 AU with inner and outer radii Rin = 0.3 AU and Rout = 3 AU, respectively, and a power-law temperature distribution in each zone given by Zone 1: Td = Tin (R/Rin )1 ; T
in

= 1800 K; 1 = -0.5;

Zone 2: Td =Tin (Rs /Rin )1 (R/Rs )2 ; 2 = -0.33. If we choose, for example, a stronger temperature gradient, then the effective average disk size becomes smaller leading to model visibilities that are higher than those observed. The calculations show that a similar agreement between model and observations can also be obtained if the disk parameters vary within the ranges given by Rin = 0.4 ± 0.1 AU, Td (Rin ) = 1700 ± 100 K, and 1 = 0.5 ± 0.05 up to R = 0.9 AU and 2 = 0.4 ± 0.05 for R 1 AU. The calculations of the intensity distributions of the Br -emitting region and corresponding interferometric model quantities (i.e., visibilities, wavelength-differential phases, and closure phases of the model images) show that they are very sensitive to the geometrical and kinematical parameters of the disk-wind model. Furthermore, they also depend on the disk inclination. We used the following ap-


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

0
25 Y/R 0 -25 -25
25 Y/R 0 -25

0

25

-1 0

10

25

-2 0

20

0 -25
12000

Fig. 6. Normalized Br line profiles of the disk-wind model 5 from Table A.1 at the inclination angles i = 10 (dotted line), 20 (solid), 40 (dashed), and 60 (dashed-dotted; high spectral resolution, not degraded to the spectral resolution of AMBER). Table 3. Range of parameter variations for our continuumdisk plus disk-wind model calculations
10000

Y/R

25

- 30

30

Y/R

0 -25

8000

25 Y/R 0 -25

-4 0

40
6000

Parameters Disk: Rin Rout Rs 1 2 Tin Disk wind: 1 N f 1 Mw

Range 0.25­3 AU (8.8­105 R ) 1­5 AU (35­175 R ) 0.85­1.25 AU (30­44 R ) -0.4­ -0.75 -0.34­ -0.4 1400­2000 K 0.1­3 AU (3.5­105 R ) 0.5­5.7 AU (17.5­200 R ) -1­ 5 0.5­3 0.3­2 10 ­80 10-9 ­10-6 M yr-1

Model 5 0.3 AU (10.5 R ) 3 AU (105 R ) 0.9 AU (31.5 R ) -0.5 -0.33 1800 K 0.5 AU (17.5 R ) 1 AU (35 R ) 2 0.5­3 1 80 10-7 M yr-1

4000

25 Y/R 0 -25

-50

50
2000

0

25 Y/R 0 -25

-60

60

-100
25 Y/R 0 -25 -25 0 X/R 25 -25 0 X/R

100

25

Fig. 5. Intensity distributions of our best-fit disk-wind model 5 (i.e., intensity distribution of the continuum disk plus the disk wind; the central star is not shown; see Tables 3 and A.1) at the center of the Br line (v = 0 km s-1 ) and at 14 other velocities (the labels give the velocity in km s-1 ). For the calculation of the model images in this figure, a clockwise motion of the disk wind was assumed. Therefore, in the blue-shifted images (left panels), mainly the disk regions on the right hand side of the star are bright. The radius of the inner edge of the disk wind ejection region (i.e., radius of the inner hole) is 1 = 17.5 R (0.3 AU). The inclination angle (angle between the polar axis and the viewing direction) of the model is i = 20 (i.e., almost pole-on). The colors represent the intensity in erg ster-1 s-1 °-1 cm-2 . In these images, AMBER's A spectral resolution of 12 000 was modeled, as described in Sect. A.3.

proach to find models that can approximately reproduce the observations: 1) To select the most appropriate the disk parameters, we calculate the continuum intensity of the disk, Id , at wavelengths near Br and compare it with the stellar intensity I . A good solution has to reproduce the observed ratio Id /I , the observed K ­band flux, and the observed visibilities. 2) We determine the velocity and density distribution along each streamline corresponding to the wind geometry and kinematics parameters (see Sect. 4 and Appendix A). 3) Using the velocity and density distribution, and choosing the temperature law along each streamline (we use a constant electron temperature of 8000 K), we solve the equations of the statistical equilibrium for the hydrogen atoms and compute the population of the hydrogen atomic levels and the ionization degree along each streamline. 4) We then calculate the Br line profile and compare it with the observed line profile. The solution is good if the theoretical line profile is in good agreement with the intensity and shape of the observed one. 5) For each disk-continuum and disk-wind model (see Table 3), we compute the two-dimensional intensity distribution map for several inclination angles. We calculate all
7


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

model intensity distributions for 15 different wavelengths across the Br line and the continuum. The spectral resolution of AMBER of R = 12 000 is modeled in the final model intensity distributions (see discussion on the importance of modeling AMBER's spectral resolution by convolution in Section A.3 and Fig. A.3). 6) Finally, we compute visibilities, wavelengthdifferential phases, and closure phases of the obtained 2D model intensity distributions, which are needed for the comparison with the observations. 4.4. Comparison of the models with the observations Using this method, we found that model 5 (see Tables A.1 and 3) provides the closest agreement with the AMBER data (see Figs. 7, A.4, and A.5) and the K -band continuum flux, as discussed below. The inner radius of the diskwind ejection region of model 5 is 0.5 AU (1 = 17.5 r ) and the half-opening angle of the disk wind is 80 . The disk-wind half-opening angle is defined as the angle between the rotation axis and the innermost wind streamline (see Fig. A.1). Figure 7 shows flux, visibilities, wavelengthdifferential phases, and closure phases of model 5 for the same interferometric baselines as the AMBER observations presented in Fig. 1. The spectral energy distribution of model 5 predicts that log F [erg cm-2 s-1 ] = -7.13 in the K -band continuum near Br for a distance of 250 pc, which is in good agreement with the dereddened SED of MWC 297 from Acke et al. (2008) and Alonso-Albi et al. (2009). The model ratio Id /I = 10 also agrees with the observations. Since the position angle of MWC 297's pro jected polar disk axis on the sky is unknown, the interferometric observables of the best-fit model 5 (spectrum, visibilities, wavelength-differential phases, and closure phases) were calculated for all PAs in steps of 10 (see Sect. A.3). Figure 7 compares the observations (left panel) with two of the models that are approximately in agreement with all observables (clockwise motion, i = 20 , middle: PA = 65 , right: PA = 300 ). A similar agreement can also be obtained for some smaller inclination angles (see Fig. 8). The detailed PA dependence is presented in Sect. A.3 and Figs. A.4 and A.5, which show the following results. For models with clockwise motion of the disk wind and i 20 (i.e., almost pole-on), we obtained the closest agreement between model and observations for the PAs of approximately 65, 120, 245, and 300 (±5 ; see Fig. A.4). For the models with anti-clockwise motion and i 20 , we obtained the best agreement for the PAs of approximately 15, 70, 195, and 250 (see Fig. A.5).

5. Discussion
Our disk-wind model allows us to approximately reproduce all interferometric observables, including the remarkably strong differential and closure phases measured in our AMBER/VLTI observations of MWC 297 and to draw several conclusions: 1) The AMBER continuum visibilities confirm previous results (Monnier et al. 2005; Malbet et al. 2007; Acke et al. 2008; Kraus et al. 2008a) that the continuum-emitting region is remarkably compact. We used geometric models consisting of a stellar point source, an extended halo, and
8

a dominant compact Gaussian or ring to characterize its size. For the compact ring component, we obtained an inner ring-fit radius of 2.2 mas (0.56 AU). This compact ring radius is about 5.4 times smaller than the 3 AU dust sublimation radius expected for silicate grains, if there is no radiation-shielding material between the star and the dust rim. 2) The Br emission line region is far more extended than the compact continuum-emitting region. The strong and wavelength-dependent closure phase reveals that the line-emitting region is very asymmetric along our measured interferometric position angle and the asymmetry dramatically changes across the Br line. At the center of the Br line, we derived a Gaussian fit HWHM radius of 6.8 mas (1.6 AU), which is 2.7 times larger than the compact 0.56 AU continuum Gaussian. The diameters of the Br -line-emitting region shown in Fig. 3 are derived from continuum-compensated visibilities. Therefore, they are the diameters of the Br -line-emitting region (not of the emission region of both line and continuum). 3) To interpret our AMBER observations, we employed a magneto-centrifugally driven disk-wind model. This model consists of an accretion disk, which emits the observed continuum radiation, and a disk wind, which emits the Br line. The disk wind starts at a distance of 0.5 AU from the star, which is further out than the typical value of 0.07 AU suggested for the inner radius of T Tauri star disk winds (Safier 1993; Ferreira 2007). The required large half-opening angle (80 ) of the wind in our model is much larger than the typical value (30­45) considered for T Tauri stars (see, e.g., Safier 1993; Ferreira 2007). Both of these features of the model are likely to be caused by the strong influence of the radiation pressure of the star (Drew et al. 1998). The model images of the disk wind are highly asymmetric and their photocenter shift and shape strongly change across the line profile, resulting in a strong wavelength dependence of the visibilities, wavelength-differential phases, and closure phase. This asymmetry of the model images is a result of the Doppler shift of the radiation caused by the gas motion in the emitting region. The Doppler shift strongly depends on the position within the disk and the disk inclination angle. The properties of the inner region of the circumstellar environment of MWC 297 are remarkable in many respects. In the NIR, it radiates at temperatures of Td 2000 K at a distance of only 0.5 AU from the star. If the disk is optically thick between the star and the sublimation zone (i.e., the disk is able to shield the radiation from the cental star), the radius of this zone, Rsub , will be smaller than for the optically thin limit. Using the model of optically thick and geometrically thin disks (Adams & Shu 1986) heated by the central star, one obtains Rsub 0.8 AU for MWC 297, which is comparable to the inner disk radius of the best-fit model 5. The mass-loss rate of the disk wind of the best-fit model 5 is 10-7 M yr-1 (Table A.1). The theory of magneto centrifugal disk winds predicts a typical ratio Mw /Macc on the order of 0.1 (K¨nigl & Pudritz 2000). In this case, o the mass accretion rate in the disk of MWC 297, Macc , is -6 -1 10 M yr . According to Muzerolle et al. (2004), such a disk is optically thick between the star and the sublimation zone because of the gas opacity.


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

velocity (LSR) [km/s] -100 0 100 5.0 4.0 3.0 Flux Flux 2.0 1.0 1.0 0.9 0.8 0.7 Visibility Visibility 0.6 0.5 0.4 0.3 0.2 0.1 0.0 160 120 Diff. Phase [deg] 80 40 0 -40 -80 -120 -160 160 Closure Phase [deg] Closure Phase [deg] 120 80 40 0 -40 -80 -120 -160 2.163 2.164 2.165 2.166 2.167 2.168 2.169 Wavelength [µm] PBL 14.0m, PA PBL 28.0m, PA PBL 42.1m, PA 68
o 68 o 68o

velocity (LSR) [km/s] -100 0 100 5.0 PA 65 deg 4.0 3.0 2.0 1.0 1.0 0.9 0.8 0.7 0.5 0.4 0.3 0.2 0.1 0.0 160 120 Diff. Phase [deg] 40 0 -40 -80 -120 -160 160 80 40 0 -40 -80 -120 -160 2.163 2.164 2.165 2.166 2.167 2.168 2.169 Wavelength [µm] Closure Phase [deg] 120 Diff. Phase [deg] 80 Visibility 0.6 Flux 4.0 3.0 2.0 1.0 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 160 120 80 40 0 -40 -80 -120 -160 160 120 80 40 0 -40 -80 -120 -160 5.0

velocity (LSR) [km/s] -100 0 100 PA 300 deg

2.163 2.164 2.165 2.166 2.167 2.168 2.169 Wavelength [µm]

Fig. 7. Comparison of the observations (left panel; see Fig. 1) with the corresponding model quantities of model 5 for two different PAs of the model image (middle and right panel; the model quantities are calculated for AMBER's spectral resolution of 12 000, as discussed in Sect. A.3 and Fig. A.3). The middle and right panels show the dependence of the interferometric observables (spectrum, visibilities, wavelength-differential phases, and closure phases) of our best-fit diskwind model 5 (disk-wind emitting region, continuum accretion disk, plus central star; Tables 3 and A.1) on the wavelength across the Br line for an inclination angle of i = 20 , clock-wise motion of the disk wind, and for two different PAs of the pro jected disk polar axis on the sky: 65 (middle) and 300 (right). Several other PAs that are also approximately in agreement with the observations are discussed in the text and Figs. A.4 and A.5. The detailed dependence of the interferometric observables on the PA is presented in Figs. A.4 and A.5. What type of matter radiates in this compact inner disk inside the dust sublimation radius? We assume that it is a mixture of the warm, mostly neutral molecular gas plus refractory (e.g. graphite) grains. Similar properties of the inner compact material were found by Benisty et al. (2010) for the HAeBe star HD 163296. Interaction of the dust grains with the stellar radiation possibly plays an important role in the acceleration of dust and gas and the formation of the disk wind. As shown above, an is needed to explain the be filled with material semi-transparent in the inner gap in the disk of MWC 297 observations. This gap region may that is either fully transparent or infrared (e.g., Tannirkulam et al. 2008). The formation of this gap could be a result of binarity (Artymowicz & Lubow 1994). Another possibility for gap formation is the interaction of stellar radiation and wind with the inner disk. The strong radiation pressure and the stellar wind can blow away the disk atmosphere (Drew et al. 1998). As a result, the disk can partially or fully dissipate in the vicinity of the star. 4) The critical rotation velocity vcrit = (2GM /(3R ) (Maeder & Meynet 2000) of MWC 297 would be 450 km s-1 for R = 6 R and M = 10 M , as adopted in our paper. For our determined inclination angle i = 20 ± 10 , the maximum value of vrot sin i would be 150 ± 80 km s-1 for vrot < vcrit . This is in clear contradiction
9


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

velocity (LSR) [km/s] -100 0 100 5.0 PA 300 deg 4.0 3.0 Flux Flux 2.0 1.0 1.0 0.9 0.8 0.7 Visibility Visibility 0.6 0.5 0.4 0.3 0.2 0.1 0.0 160 120 Diff. Phase [deg] Diff. Phase [deg] 80 40 0 -40 -80 -120 -160 160 Closure Phase [deg] Closure Phase [deg] 120 80 40 0 -40 -80 -120 -160 2.163 2.164 2.165 2.166 2.167 2.168 2.169 Wavelength [µm] 4.0 3.0 2.0 1.0 1.0 0.9 0.8 0.7 0.5 0.4 0.3 0.2 0.1 0.0 160 120 40 0 -40 -80 -120 -160 160 80 40 0 -40 -80 -120 -160 120 80 5.0

velocity (LSR) [km/s] -100 0 100 5.0 PA 300 deg 4.0 3.0 Flux 2.0 1.0 1.0 0.9 0.8 0.7 Visibility 0.6 0.5 0.4 0.3 0.2 0.1 0.0 160 120 Diff. Phase [deg] 80 40 0 -40 -80 -120 -160 160 Closure Phase [deg] 120 80 40 0 -40 -80 -120 -160

velocity (LSR) [km/s] -100 0 100 PA 300 deg

0.6

2.163 2.164 2.165 2.166 2.167 2.168 2.169 Wavelength [µm]

2.163 2.164 2.165 2.166 2.167 2.168 2.169 Wavelength [µm]

Fig. 8. Dependence of the interferometric quantities (visibilities, wavelength-differencial phases, and closure phases) of our best-fit disk-wind model 5 (see Table A.1) on the inclination angle i: 10, 15, and 20 from left to right for the PA 300 of the disk (see Fig. 7). In contrast to Fig. 7, the Br line profiles are given without convolution with the instrumental profile of the AMBER spectrograph for clearer visibility of the theoretical line profiles. with the value of vrot sin i = 350 ± 50 km s-1 obtained by Drew et al. (1997) for MWC 297. What is the cause of this discrepancy? Is it our poor knowledge of either mass or radius, the inclination angle of 20 ± 10 , or the observed vrot sin i = 350 ± 50 km s-1 ? The value of i is constrained not only by our modeling but also by the recent interferometric studies reported by Malbet et al. (2007) and Acke et al. (2008). Nevertheless, further studies of i will be required to constrain i more reliably. Furthermore, it cannot be excluded that the value of vrot sin i = 350 ± 50 km s-1 is overestimated, since it is difficult to derive the exact rotational velocity from the width of the observed absorption lines because of the medium spectral resolution employed and the quite low SNR of the Drew et al. spectrum caused by the relatively faint magnitude of V = 12.2. Moreover, the absorption lines in MWC 297 might be: a) contaminated by emission from outflowing material, as is often the case in massive young stars, or b) formed at a certain distance from the photosphere, where the broadening of the spectral line would be dominated by the outflow kinematics. Nevertheless, this question is far
10

from being solved, and we encourage independent determinations of vrot sin i using, for instance, high-resolution spectroscopy and state-of-the-art atmospheric modeling. Future multi-wavelength observations of MWC 297 using high spectral and spatial resolutions, combined with nonLTE radiative transfer modeling of the stellar atmosphere, will be essential to improve our knowledge of the stellar parameters such as the rotational velocity. 5) According to our data, the Br line has a low redshift in the LSR frame (but not in the heliocentric system), in contrast to the model spectrum. However, it is not yet known whether the Br - emitting region is really redshifted with respect to the star, since the small observed redshift could simply be caused by the (unknown) radial velocity of the star (see Appendix D). The model of the emitting region of MWC 297 proposed in this paper needs to be improved. In particular, we need to investigate the dynamics of the magneto-centrifugal disk wind in the presence of the strong radiation field of the star. It is also necessary to investigate the effect on the thermal


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

regime of the wind of both stellar radiation and ambipolar diffusion. Additional information about the properties of the emitting region can be obtained by analyzing other hydrogen lines, first of all the H line. The observations by Drew et al. (1997) show that in the spectrum of MWC 297, this line is broader than the Br line and has symmetric, extended wings up to 400­500 km s-1 . The Br line is also much narrower than H. Our preliminary test calculations show that this difference can be the result of opacity effects ( (H) (Br )) in addition to Stark broadening. A more detailed discussion of the H line properties as well as those of other Balmer lines will be presented in a forthcoming paper.

We are grateful to Rene Oudmaijer and Bram Acke for sharing their spectroscopic data. Finally, we thank the referee for his comments that helped to improve the manuscript.

App endix A: Description of the mo del calculations
A.1. Disk wind For the description of the disk wind, we use a coordinate system (l, , ) centered at point S located on the rotation axis (see the sketch in Fig. A.1). In MHD models of the disk wind, the inclination angle of the first streamline with respect to the disk plane is typically assumed to be around 60 (Blandford & Payne 1982). However, in the case of Herbig Be stars, the strong stellar radiation pressure possibly bends down the streamlines and makes the disk wind flatter. Therefore, in our models we use a smaller inclination angle for the innermost streamline (e.g., between 10 and 30 ). This effect influences the characteristic size of the wind-launching region. As mentioned in Section 4, the tangential velocity of the wind changes along streamlines as u( ) = uK (i ) ( /i )-1 , (A.1)

6. Conclusions
We have presented VLTI/AMBER observations of the Herbig Be star MWC 297 with a spectral resolution of 12 000. This high spectral resolution has allowed us to study the structure and kinematics of the inner sub-AU accretion-ejection region in approximately 10 different spectral channels (radial velocities) across the resolved Br emission line and compare the observations with disk-wind models. The K -band continuum observations confirm that MWC 297's continuum emitting region is remarkably compact. We measured a ring-fit radius of 2.2 mas (0.56 AU) for the dominant compact component. At all Doppler velocities across the Br line, the line-emitting region is both asymmetric (as indicated by the strong closure phase) and more extended than the compact continuum-emitting region. At the center of the Br line, we derived a Gaussian continuum-compensated radius (HWHM) of 6.3 ± 0.4 mas (1.6 ± 0.1 AU) for the line-emitting region, which is 3 times larger than the 0.58 AU HWHM of the compact continuum-emitting Gaussian component. To interpret our AMBER observations, we employed a magneto-centrifugally driven disk-wind model consisting of an accretion disk, which emits the observed continuum radiation, and a disk wind, which emits the Br line. We have demonstrated that all observables, i.e., the observed Br line profile, K -band continuum flux, visibilities, wavelength-differential phases, and remarkably strong closure phases, can simultaneously be reproduced by our disk-wind model. The modeling also shows that disk-wind models developed for T Tauri stars can be used to interpret Herbig Be star observations if suitable modifications are applied. The modification has to take into account the higher radiation pressure from hot stars, which can change the inclination angle of the streamlines of the disk wind and cause a flatter disk wind than in T Tauri stars. The disk wind of MWC 297 starts at a distance of 0.5 AU (17.5 stellar radii) from the star, which is a larger distance than the typical value of 0.07 AU predicted for T Tauri stars.
Acknow ledgements. We thank the ESO VLTI team on Paranal for the excellent collaboration. The data presented here were reduced using the publicly available data-reduction software package amd lib kindly provided by the Jean-Marie Mariotti Center (http://www.jmmc.fr/data processing amber.htm). The telluric spectra used in this work for spectral calibration of the AMBER data were created from data that was kindly made available by the NSO/Kitt Peak Observatory. A.M., V.G., and L.T. thank the Max-PlanckSociety for support during their stay in Bonn. This publication makes use of the SIMBAD database operated at CDS, Strasbourg, France.

where = l sin is the distance of the point (l, ) from the rotation axis, i = li sin , and uK (i ) = (G M /i )1/2 at the point (i ) at the base of streamline i. The radial velocity v increases along the streamlines as v (l) = v0 + (v - v0 ) (1 - li /l) , (A.2)

where v0 and v are the initial and terminal values of the radial velocity and is a parameter. We adopt v = f uK (i ), where uK (i ) is the Keplerian velocity at distance i from the disk axis and f is the scale factor of the asymptotic terminal velocity to the local Keplerian velocity at the wind emerging point (see Kurosawa et al. 2006). The parameter f can change with distance because the magnetic field in the accretion disk decreases with . An example of the velocity field of one of the models considered below is shown in Fig. A.2. The kinematical model described above is typical of disk-wind models adopted for different astrophysical ob jects with accretion disks including T Tauri stars (Kurosawa et al. 2006), cataclysmic variables (Shlosman & Vitello 1993), and even active galactic nuclei (Murray & Chiang 1997). In these models, the gas streams that start in the innermost regions of the disk-wind launching region have the highest radial velocities. These kinematical models are characterized by non-local radiative coupling between distant parts of the moving media (see, e.g., Grachev & Grinin 1975; Rybicki & Hummer 1978), which complicate the radiative transfer problem. In our calculations, we neglect the influence of the non-local radiative coupling since its effect on the radiative excitation in the flattened emitting regions (such as accretion disks or disk winds with large opening angles) is small. For the calculation of the density distribution in the disk wind, we use the mass-loss rate per unit solid angle, Mw ()/4 . One can show that Mw () and the mass-loss rate per unit area of the disk, mw , are related by Mw () = 4 d2 mw ( ) cos-3 , (A.3)
11


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER











Fig. A.1. Sketch of the geometry of the disk-wind model adopted.

Fig. A.2. Left: Radial velocity of the outflow along the different streamlines as a function of l; numbers 1-6 denote the number of the streamline. Right: The same for the rotational velocity; the dashed line shows the Keplerian velocity in the disk. where d is the distance between point S and the star. The total mass-loss rate is then Mw =

1 N

Mw () sin d.

(A.4)

is the integrated line opacity in the considered spectral line, (, r) is the line optical depth at point r and for the frequency in the direction to an observer, and ni is the number density of atoms in the i-th state
z
max

In this case, we can write the continuity equation in its usual form 4 (l, )v (l, )l2 = Mw (). (A.5) Using this equation, one can calculate the distribution of the number density at each point in the disk wind. A.2. Calculation of the model intensity distributions For the calculations of the 2-D intensity distributions (maps) of the emitting region, we use the coordinate system (x, y , z ) centered on the star. The (x, y ) plane coincides with the sky plane, the x axis is the intersection of the disk plane with the sky plane, and the z axis is parallel to the line-of-sight. In this case, the intensity of radiation at a frequency within a spectral line is
z
max

(, r) =
z

(r )( -

0

vz (r ) )dz , c

(A.8)

Iw (, x, y ) =
z
min

S (r)( -

0

vz (r) )e c

- (,r)

(r)dz , (A.6)

where r is a vector, |r| = (x2 + y 2 + z 2 )1/2 , vz (r) is the pro jection of the velocity at point r on the line-of-sight, and S is a source function for a transition between energy levels i and j S (r) = 2 h c2
3

nj (r) gi -1 ni (r) gj

-1

,

(A.7)

where is the profile function normalized to unity, r is the vector with coordinates (x, y , z ), and z changes from the current value z in Eq. A.6 to the outer boundary value z m ax . The optical depth of the disk wind in the line frequencies in the direction toward the star is ( ) = (, r = 0). We assume that there is a complete redistribution of the line frequencies in the reference frame of the atom and use the Doppler profile in the calculations of the Br and the Voigt profile for the H line (with the same constants as in Kurosawa et al. 2006). The calculations of the ionization state and the number densities of the atomic levels were performed in a cylinder with a radius rc up to 300 R and a height hc = 600 R divided into grid cells in l, coordinates. In each cell, we solved the equations of the statistical equilibrium for the 15level hydrogen atom + continuum, taking into account both collision and radiative processes of excitation and ionization (see Grinin & Mitskevich 1990 for more details). We assumed that the distribution of the atomic sub-levels follows their statistical weights. Johnson collision rates (Johnson 1972) were used for all transitions except 1-2; for the latter, we used those from Scholz et al. (1990).

12


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

velocity (LSR) [km/s] -100 0 100 5.0 PA 300 deg 4.0 3.0 Flux 2.0 1.0 1.0 0.9 0.8 0.7 Visibility 0.5 0.4 0.3 0.2 0.1 0.0 160 120 Diff. Phase [deg] 40 0 -40 -80 -120 -160 160 Closure Phase [deg] 80 40 0 -40 -80 -120 -160 2.163 2.164 2.165 2.166 2.167 2.168 2.169 Wavelength [µm] Closure Phase [deg] 120 Diff. Phase [deg] 80 Visibility 0.6 Flux 4.0 3.0 2.0 1.0 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 160 120 80 40 0 -40 -80 -120 -160 160 120 80 40 0 -40 -80 -120 -160 5.0

velocity (LSR) [km/s] -100 0 100 PA 300 deg

I ( )exp(- ( )) to this sum. Our calculations have shown that the contribution of the free-bound and free-free radiation of the disk wind at the wavelength of the Br line is small compared with the disk radiation. To model observations with a spectral resolution lower than the theoretical one, the intensity maps have to be convolved with the instrumental profile
+

Ic (, x, y ) =
-

I ( , x, y )i( - ) d ,

(A.11)

where i( ) is the normalized instrumental profile of AMBER. In our case, the spectral resolution is quite high (25 km s-1 ). The calculations of the theoretical maps (Fig. 5) and the interferometric model quantities (Figs. 7, A.4, and A.5) were made using Eq. A.11. A comparison of the results obtained with and without spectral convolution shows (see section A.3) that noticeable differences in the interferometric quantities exist even at this high spectral resolution. The flux at distance D from the ob ject is an integral over the entire area in the x,y -plane of the emitting region F =
Wavelength [µm]

1 D2

I (, x, y )dx dy .
A

(A.12)

2.163 2.164 2.165 2.166 2.167 2.168 2.169

To summarize, the model assumptions are: ­ The disk wind is launched from the accretion disk surface in the range from distance 1 (the starting point of the innermost streamline) to N (starting point of the outermost streamline). The half-opening angle (which is the angle between the rotation axis and the innermost streamline) is 1 (see Fig. A1). ­ The wind is assumed to be isothermal with an electron temperature of 8000 K (for all models in Table A.1). ­ The disk wind possesses both radial and tangential (rotation) motion. ­ The number density distribution is calculated from the mass continuity equation for each streamline using the local mass-loss rate per unit area. ­ The atomic level populations were computed in the Sobolev approximation; a 15-level hydrogen atom has been considered. The 3­D radiative problem (exact integration) has been solved in a cylinder with a radius of 300 R and a height 600 R . A.3. Interferometric model quantities

Fig. A.3. Comparison of the theoretical interferometric quantities of model 5 obtained with (left) and without (right) modeling of AMBER's spectral resolution of 12 000. In both cases, position and inclination angles are identical (PA = 300 , i = 20 ). The main sources of continuum radiation near the Br line are: i) the circumstellar disk, ii) the stellar atmosphere, and iii) the free-free and free-bound radiation of the disk wind. If the radiation of the star (I ) intersects the disk wind in the line of sight, it has to be corrected for the absorption in the wind and we observe I exp(- ( )). In the adopted coordinate system, we obtain for the observed disk intensity distribution Id (, x, y ) = B (Td (R)) e
-w (,x,y )

,

(A.9)

where R = (x2 + y 2 / cos2 i)1/2 , w (, x, y ) is the optical thickness of the disk wind at frequency between an observer and the point in the disk with the coordinates x = x and y = y / cos i, and i is the inclination angle (i = 0 corresponds to the pole-on orientation). It is assumed that the disk is flat and optically thick; each elementary area of the disk radiates as a black body with the local temperature Td ( ). The total intensity map of the ob ject at frequency is the sum I (, x, y ) = Iw (, x, y ) + Id (, x, y ), (A.10)

Table A.1. Five examples of the computed disk-wind models.
Model 1 2 3 4 5 Mw 10-8 M /y r 3 3 1.2 0.6 10 w deg. 70 10 80 80 80 1 r 20 30 25 10 17.5 N r 80 200 200 50 35 0 1 1 0 2 f 3 2 2 * *

where Iw (, x, y ) is the intensity of the disk wind at the point x, y at an emission line frequency , and Id is the disk intensity distribution. At the central part of the map where x2 + y 2 r , we have to add the intensity of the stellar radiation

Note: * In models 4 and 5, f is different for different streamlines (see Sect. A.1)

13


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 5.0 PA 0 deg 4.0 3.0 Flux 2.0 1.0 1.0 0.9 0.8 0.7 Visibility 0.6 0.5 0.4 0.3 0.2 0.1 0.0 160 120 Diff. Phase [deg] 80 40 0 -40 -80 -120 -160 160 120 Closure Phase [deg] 80 40 0 -40 -80 -120 -160 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] PA 10 deg PA 20 deg PA 30 deg PA 40 deg PA 50 deg PA 60 deg PA 70 deg PA 80 deg

velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 5.0 PA 90 deg 4.0 3.0 Flux 2.0 1.0 1.0 0.9 0.8 0.7 Visibility 0.6 0.5 0.4 0.3 0.2 0.1 0.0 160 120 Diff. Phase [deg] 80 40 0 -40 -80 -120 -160 160 120 Closure Phase [deg] 80 40 0 -40 -80 -120 -160 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] PA 100 deg PA 110 deg PA 120 deg PA 130 deg PA 140 deg PA 150 deg PA 160 deg PA 170 deg

Fig. A.4. Dependence of the visibilities and phases of model 5 (see Tables 3 and A.1) on both the wavelength across the Br line and the PA of the model on the sky for an inclination angle of i = 20 and clockwise motion of the disk wind (the corresponding anti-clockwise plots are shown in Fig. A.5). Top: PAs from 0 to 80 in steps of 10 from left to right; Bottom: PAs = 90 to 170 . Closest agreement between observation and model is obtained for PAs of 65, 120, 245, a4 d 300 , as discussed in Sect. A.3. n 1


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 5.0 PA 0 deg 4.0 3.0 Flux 2.0 1.0 1.0 0.9 0.8 0.7 Visibility 0.6 0.5 0.4 0.3 0.2 0.1 0.0 160 120 Diff. Phase [deg] 80 40 0 -40 -80 -120 -160 160 120 Closure Phase [deg] 80 40 0 -40 -80 -120 -160 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] PA 10 deg PA 20 deg PA 30 deg PA 40 deg PA 50 deg PA 60 deg PA 70 deg PA 80 deg

velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] velocity (LSR) [km/s] -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 -100 0 100 5.0 PA 90 deg 4.0 3.0 Flux 2.0 1.0 1.0 0.9 0.8 0.7 Visibility 0.6 0.5 0.4 0.3 0.2 0.1 0.0 160 120 Diff. Phase [deg] 80 40 0 -40 -80 -120 -160 160 120 Closure Phase [deg] 80 40 0 -40 -80 -120 -160 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] 2.165 2.166 Wavelength [µm] PA 100 deg PA 110 deg PA 120 deg PA 130 deg PA 140 deg PA 150 deg PA 160 deg PA 170 deg

Fig. A.5. Same as Fig. A.4, but for anti-clockwise motion. The best agreement between model and observations is obtained for the PAs of approximately 15, 70, 195, and 250 (see figure caption A.4).

15


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

Fig. A.6. Disk-wind model profiles of the Br line for different inclination angles (from top to bottom: 20, 40, and 60 ). From left to right, the profiles correspond to models 1, 2, 3, and 4 in Table A.1 (model 5 is discussed in Fig. 6). Using the theoretical intensity distribution maps of the emitting regions, one can calculate the interferometric quantities of the model: spectrum, visibilities, wavelengthdifferential phases, and closure phases for the same interferometric baseline lengths and position angles as the AMBER observations. Examples of interferometric quantities of disk-wind model 5 (see Table A.1) are shown in Figs. 7, 8, A.3, A.4, and A.5. In Figure A.3, we show the interferometric model quantities of the Br line calculated for model 5 with and without convolution of the theoretical intensity maps with the spectral AMBER profile corresponding to spectral resolution 12 000 (see Eq. A.11). One can see that in both cases, the visibilities inside the line profiles are very similar for all three baselines. However, the differential phases for the largest baseline (42 m) and the closure phases are quite different. In Figures A.4 and A.5, we show the detailed dependence of the visibilities and phases of model 5 (see Tables 3 and A.1) on both the wavelength across the Br line and the PA of the model on the sky for both an inclination angle of i = 20 and clockwise (Fig. A.4) and anti-clockwise (Fig. A.5) motion of the disk wind (the PA is counted from north to east with north up and east to the left). In Fig. A.4, the closest agreement between observation and model is obtained for PAs of 65 and 120 . The different sign of the phase at the PA = 120 does not disagree with the AMBER observations since the sign is unknown. The plots for the PAs of 180 to 350 are not shown since they give the same results as the PAs from 0 to 170 , except an opposite sign of the differential and closure phases, since a 180 rotation leads to a conjugate complex Fourier transform. Therefore, the PAs of 65 + 180 = 245 and 120 + 180 = 300 agree with the observations in addition to the PAs of 65 and 120 . For anti-clockwise motion (Fig. A.5), the closest agreement be16

tween model and observations is obtained for the PAs of approximately 15, 70, 195, and 250 . A.4. Br line as a function of the model parameters Examples of the theoretical Br line profiles calculated for the models listed in Table A.1 are presented in Fig. A.6. It shows that even for large inclination angles i, the profiles are quite narrow. This is due to the large radius of the disk-wind launching region. Except for a small range of i values close to pole-on view (i 20 ), all orientations of the disk result in a double-peaked Br line. In observations, the double-peaked structure can, of course, appear single-lined if the line is narrow and the spectral resolution is too low. A single-peaked Br emission line can also be obtained if a weak, narrow emission component is added to the model. Such emission can arise at large distances from the star, as discussed in the following Sect. A.4.1. A.4.1. Low-velocity component of the disk wind A low-velocity outflow from circumstellar disks can be caused by the heating of the upper layers of the disk atmosphere by EUV (h > 13.6 eV) radiation from the central star. As a result, the thermal velocity of atoms and ions in the outer part of the disk can exceed the escape velocity. This is the region where the disk can evaporate (Hollenbach 2 et al. 1994). The inner radius of this zone is rg = GM /vs , where vs is the sound speed. For ionized gas and M = 10 M , rg is 7â104/Te 70 AU (Gorti & Hollenbach 2009). This value exceeds the radius of the AMBER FOV (120 mas or 30 AU at the distance of MWC 297), meaning that the photo-evaporated disk wind cannot contribute to the emission observed with VLTI. However, in the extended region of the disk between the outer radius of the


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

MHD wind ejection region (1 AU) and rg , the EUV photons ionize the upper layers of the disk atmosphere and produce a mini-H II region. This part of the disk does not evaporate, but can radiate in the Br and other lines. This emission is broadened by the Keplerian rotation and thermal motion of atoms. Its contribution to the total emission may be small, but potentially enough to fill in the small depression in the cental part of the double-peaked Br line profile formed by the disk wind. We simulated this profile (Fig. A.7) by changing the parameter f from f = 3 at 1 to 0.5 at N . A more detailed modeling of this narrow emission is beyond the scope of this paper and will be presented in a separate one.

Fig. A.8. Dependence of the Br line profile on screening caused by an optically thick, dusty disk wind. The solid line represents the unscreened line profile of model 2 (see Table A.1) with f = 5, while the long-dashed (Hd = 100 R ; = 20 ) and dot-dashed (Hd = 130 R , = 25 ) lines show models with screening. For all models shown, the inclination angle is i = 70 . See text for further details.

App endix C: Beam commutation device
To optimize the calibration accuracy, we employed the AMBER Beam Commutation Device (BCD, Petrov et al. 2007), as described in Sect. 2. This device is located very early in the optical path of the AMBER instrument and allows us to exchange two of the three telescope beams within a few seconds. Performing this beam commutation between two subsequent exposures, allows one to correct potential instrumental drifts. The insertion of the BCD causes a sign change in one of the three differential phases and the CP signal (Millour et al. 2008). Our BCD OUT/BCD IN data sets precisely showed these expected sign changes. The differential phases and closure phases derived from both independent data sets (BCD OUT/BCD IN) are shown in Fig. C.1. Subtraction of the closure phases obtained with and without BCD cancels potential instrumental phase drifts. Fig. C.2 presents the derived closure phase and a fit of a constant to the continuum closure phases. The obtained averaged continuum closure phase is -1.28 ± 0.42.

Fig. A.7. Br line profile of model 4 with constant (a) and outward decreasing (b) parameter f (i = 20 ; see Section A.4.1 for details).

A.4.2. Screening effects As noted above, the optically thick part of the disk may screen some parts of the disk wind. This could lead to a red-shift of emission lines (depending on i, see Fig. A.8). The calculations show that this effect strongly depends on the opening angle 1 of the disk wind. At large i, the approaching part of the disk wind can be screened from an observer by optically thick layers of the flared disk. We simulated this screening effect by introducing a screen with a sharp upper boundary at some distance from the center. The angular size of the screen seen from the stellar position is the parameter that we introduced. The calculations show (see Fig. A.8) that for some (depending on the disk-wind parameters) and i 70 , the screening effect leads to a red-shifted emission line. It is likely that this mechanism is not important in our case, since the observations suggest that MWC 297's inclination angle is small (10 to 20 ).

App endix D: Optical sp ectroscopy
MWC 297 has been rarely observed spectroscopically mostly because of its relative faintness (V 12.2 mag) and high reddening (B - V 2.1 mag). Only two spectra that cover almost the entire optical range, although with a low resolution of 0.8 °, have ever been published (Drew et al. A 1997; Andrillat & Jaschek 1998). Higher-resolution spectroscopy was obtained for a few regions, mainly centered on strong emission or absorption lines (He i 5876 °, Na i A D-lines, [O i] 6300 °, and H, Zickgraf 2003; Acke et al. A 2005). A summary of the published optical spectroscopic data is presented in Table D.1. A detailed study of the ob ject's spectrum is beyond the scope of this paper. Thus, we only focus on the details that are relevant to the modeling of our interferometric results. One important issue for this purpose is the geometry of the circumstellar environment. The strongest emission-line profiles exhibit a double-peaked structure, which can only be revealed at a resolving power of R 30 000 (see Zickgraf 2003; Acke et al. 2005). This prompted Zickgraf (2003) to suggest that MWC 297 might be surrounded by a disk that
17

App endix B: Wavelength calibration
Using the telluric spectrum provided by NOAO as a reference (see Fig. B.1), standard IRAF routines and specific IDL routines developed by JHG were used to find a solution to the pixel-to-wavelength calibration of the AMBER spectrum (Fig. B.2). The estimated uncertainty in the wavelength calibration amounts to 2.5 km s-1 (Fig. B.1).


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

vLSR [km/s] -1000 1.0 normalized raw spectrum [arbitrary units] 0 1000

0.8

1.0 0.9 atmospheric transparancy atmospheric transparancy

0.6 0.8 0.4 0.7 0.6 0.2 0.5 0.0 0.4 2.145 2.15 2.155 2.16 2.165 2.17 2.175 2.18 2.185 2.19 2.195 [µm] vLSR [km/s] -100 1.0 0 100

normalized raw spectrum [arbitrary units]

0.8

1.0 0.9

0.6 0.8 0.4 0.7 0.6 0.2 0.5 0.0 0.4 2.16 2.161 2.162 2.163 2.164 2.165 2.166 2.167 2.168 2.169 2.17 [µm]

Fig. B.1. Wavelength calibration of the AMBER data. Top: Raw AMBER spectrum of MWC 297 (solid blue line) compared to a telluric spectrum from the Kitt Peak Observatory. The green line shows the original Kitt Peak spectrum with a spectral resolution of 60 000, while the red line displays the telluric spectrum that has been convolved to the same spectral resolution as the AMBER data (R = 12 000). For clarity, the red curve has been shifted down by 0.4. As the figure reveals, after the wavelength calibration, the AMBER raw spectrum of MWC 297 matches several telluric features. Bottom: Same as top panel, but just for the wavelength region around the Br line. is viewed at an intermediate inclination with respect to the line-of-sight. Another crucial parameter is the systemic (stel lar) radial velocity. Since this had never been carefully discussed, we critically examined possible methods to derive this quantity and performed an independent study of the published spectra. Absorption lines of the photospheric origin are only seen in the spectra obtained by Drew et al. (1997). They include a few He i lines blueward of 5800 ° and those A of C iii and N iii in the 4500­4700 ° range. All of them are A broad, which was interpreted to be caused by a high stellar rotation rate (Drew et al. 1997). Additionally, the residual intensities of the carbon and nitrogen lines are within 5 % of the continuum, hampering high accuracy measurements even at the SNR of 120 achieved in the spectrum. The average LSR radial velocity, which we derived from five He i absorption lines at 4009, 4026, 4143, 4471, and 4713 °, is A
18

Fig. B.2. Top: Calibration between the pixel and wavelength values of the VLTI/AMBER observations of MWC 297. A least-square polynomial fit of the data is also shown (red line). Bottom: Residuals between the measured position of the telluric line spectrum obtained from the wavelength calibration of the AMBER data and the expected position from the reference telluric spectrum provided by NOAO. The root-mean square value is of the order o f 2 .5 k m s - 1 . -7 ± 7 km s-1 . However, He i lines in the spectra of B-type emission-line stars are known to be strongly affected by the circumstellar material in the dense regions close to the star. Even pure absorption lines can exhibit variations that result in the observed shifts from the photospheric positions (see Israelian et al. 1996). Therefore, these lines cannot be considered as a reliable source of the systemic velocity. The most clearly recognized non-helium photospheric line is the C iii 4650.16 ° line, which however, forms a blend with two A N iii lines at 4640.64 and 4641.90 °. This line gives a more A negative LSR radial velocity of -29 km s-1 . The accuracy of this value is at least of the order of 15 km s-1 . The line can also be affected by the circumstellar material, which also ensures that this measurement is unreliable. An alternative method for deriving the systemic velocity can be used if photospheric lines are not seen in the spectrum or there are no data of sufficient quality for these lines. It employs centroid radial velocities of emission lines, such as those of Fe ii. It has been successfully tested on complex ob jects, such as LBVs (e.g., Humphreys et al. 1989) and stars that exhibit the B[e] phenomenon (e.g., Miroshnichenko et al. 2002). The resulting systemic velocity is used to determine a kinematical distance and lumi-


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER

Table D.1. Summary of published optical spectroscopic observations of MWC 297
Date 1987 Sept. 9 1987 Sept. 13 1992 Sept. 9 1994 June 20­22 1997 1997 2002 2002 July July Apr. June 21, 23 20, 24, 25 1 25 , ° A 6535­6600 5868­5903 6290­6335 6290­6759 3876­5198 5840­8809 4119­5734 6297­8911 6279­6321 5635­6637 R / 23 000 45 000 8 000a 8 000a 5 000b 125 000 30 000 +9 ± 1e +10 ± 2f V km s- Ref.
1

+6 ± 1c +10 ± 8d

Fig. C.1. Differential phases and closure phases derived from the data obtained without (left, BCD OUT) and with (right, BCD IN) the beam commutation device. The shown signals are not yet corrected using the calibrator data. After we corrected the sign changes introduced by the BCD device, the two data sets yield the shown signals, which are nearly identical and demonstrate the high accuracy of the measured phase signals.
Closure Phase: (BCD-OUT - BCD-IN)/2
average continuum [2.147µm, 2.194µm]: -1.28° +-0.42°

1 1 1 2 3 3 4 4 5 5

20.0 Closure phase [deg]

0.0

-20.0

Column 1: observing date; Column 2: observed sp ectral range; Column 3: average resolving p ower; Column 4: average LSR radial velocity (see explanations b elow and in the text) ; Column 5: reference to the rep orting pap er. a ­ two-pixel resolution 0.8 ° A b ­ 1.3 ° resolution at = 6000 ° A A c ­ from the [O i] 6300 ° line A d ­ from emission lines of Fe i i 4923 and 5018 ° [Fe i i] 5159 A °, Si i i 5041 and 5056 °, and Si i i 5958 and 5979 ° A A A e ­ from the [O i] 6300 and 6363 ° lines A f ­ result published by Acke et al. (2008) References: 1 ­ Zickgraf (2003); 2 ­ Vieira et al. (2003); 3 ­ Drew et al. (1997); 4 ­ Andrillat & Jaschek (1998); 5 ­ Acke et al. (2005)

-40.0

-60.0 2.163 2.164 2.165 2.166 [µm] 2.167 2.168 2.169

Fig. C.2. BCD-calibrated closure phase signal obtained by subtracting the calibrated closure phases obtained from the BCD IN and BCD OUT data. The blue line is the fit of a constant to the continuum (2.147­2.194µm) closure phase outside the Br line. nosity of the ob ject that are in general good agreement with those found by other methods (e.g., spectroscopic parallax, interstellar extinction versus distance relationship). The method seems to work for any type of centrally symmetric circumstellar envelopes. Moreover, it has been shown to trace the central star motion even in binary Be stars (e.g., Harmanec et al. 2000), whose disks may contain asymmetric density perturbations (see Okazaki 1991). Drew et al. (1997) comment that the centroid LSR radial velocities of the H, H , and the He i 5876 ° lines are A located between +2 and +20 km s-1 . Other measurements, including those we made in the optical spectra and which we have at our disposal, are listed in Table D.1. The average value is +8 ± 1 km s-1 , similar to +10 ± 2 km s-1 from Acke et al. (2008). We note that the formal statistical uncertainty of these velocities may be larger, because the number of lines used is relatively small. In addition, a systematic error due to, for instance, a slight asymmetry in the line profiles (see those shown by Zickgraf 2003) cannot be excluded. Therefore, we assume that the real uncertainty in the systemic velocity estimated with this method might be as large as 5 km s-1 .

To summarize the above discussion, we conclude that the positions of the emission lines suggest a systemic radial velocity of MWC 297 of +8 ± 5 km s-1 . In contrast, the He i absorption lines give -7 km s-1 with an even larger uncertainty, because of their possible displacement from the photospheric position due to processing by the circumstellar material. The value derived from the emission lines seems to be more reliable because it is based on data from several spectra obtained over a period of 15 years. Nevertheless, it still needs to be confirmed using high-resolution data of more lines.

References
Acke, B., van den Ancker, M. E., & Dullemond, C. P. 2005, A&A, 436, 209 Acke, B., Verhoelst, T., van den Ancker, M. E., et al. 2008, A&A, 485, 209 Adams, F. C. & Shu, F. H. 1986, ApJ, 308, 836 Alonso-Albi, T., Fuente, A., Bachiller, R., et al. 2009, A&A, 497, 117 Andrillat, Y. & Jaschek, C. 1998, A&AS, 131, 479 Artymowicz, P. & Lubow, S. H. 1994, in Circumstellar Dust Disks and Planet Formation, ed. R. Ferlet & A. Vidal-Madjar, 339 Benisty, M., Natta, A., Isella, A., et al. 2010, A&A, 511, A74+ Blandford, R. D. & Payne, D. G. 1982, MNRAS, 199, 883 Calvet, N., Hartmann, L., & Strom, S. E. 2000, Protostars and Planets IV, 377 de Kool, M. & Begelman, M. C. 1995, ApJ, 455, 448 Drew, J. E., Busfield, G., Hoare, M. G., et al. 1997, MNRAS, 286, 538 Drew, J. E., Proga, D., & Stone, J. M. 1998, MNRAS, 296, L6 Eisner, J. A., Lane, B. F., Hillenbrand, L. A., Akeson, R. L., & Sargent, A. I. 2004a, ApJ, 613, 1049 Eisner, J. A., Lane, B. F., Hillenbrand, L. A., Akeson, R. L., & Sargent, A. I. 2004b, ApJ, 613, 1049

19


G. Weigelt et al.: Sp ectro-interferometry of MWC 297 with VLTI/AMBER
Everett, J. E., K¨nigl, A., & Kartje, J. F. 2001, in Astronomical o Society of the Pacific Conference Series, Vol. 224, Probing the Physics of Active Galactic Nuclei, ed. B. M. Peterson, R. W. Pogge, & R. S. Polidan, 441­+ Ferreira, J. 2007, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 723, Lecture Notes in Physics, Berlin Springer Verlag, ed. J. Ferreira, C. Dougados, & E. Whelan, 181 Gorti, U. & Hollenbach, D. 2009, ApJ, 690, 1539 Grachev, S. I. & Grinin, V. P. 1975, Astrophysics, 11, 20 Grinin, V. P. & Mitskevich, A. S. 1990, Astrofizika, 32, 383 Harmanec, P., Habuda, P., Stefl, S., et al. 2000, A&A, 364, L85 Humphreys, R. M., Lamers, H. J. G. M., Hoekzema, N., & Cassatella, A. 1989, A&A, 218, L17 Israelian, G., Friedjung, M., Graham, J., et al. 1996, A&A, 311, 643 Johnson, L. C. 1972, ApJ, 174, 227 K¨nigl, A. & Pudritz, R. E. 2000, Protostars and Planets IV, 759 o Kraus, S., Hofmann, K., Benisty, M., et al. 2008a, A&A, 489, 1157 Kraus, S., Preibisch, T., & Ohnaka, K. 2008b, ApJ, 676, 490 Kurosawa, R., Harries, T. J., & Symington, N. H. 2006, MNRAS, 370, 580 Kurucz, R. L. 1979, ApJS, 40, 1 Lynden-Bell, D. & Pringle, J. E. 1974, MNRAS, 168, 603 Maeder, A. & Meynet, G. 2000, A&A, 361, 159 Malbet, F., Benisty, M., de Wit, W.-J., et al. 2007, A&A, 464, 43 Mannings, V. & Sargent, A. I. 1997, ApJ, 490, 792 M´rand, A., Bord´, P., & Coud´ du Foresto, V. 2005, A&A, 433, 1155 e e e Millan-Gabet, R., Schloerb, F. P., & Traub, W. A. 2001, ApJ, 546, 358 Millour, F., Petrov, R. G., Vannier, M., & Kraus, S. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7013, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series Miroshnichenko, A. S., Klochkova, V. G., Bjorkman, K. S., & Panchuk, V. E. 2002, A&A, 390, 627 Monnier, J. D., Berger, J.-P., Millan-Gabet, R., et al. 2006, ApJ, 647, 444 Monnier, J. D., Millan-Gabet, R., Billmeier, R., et al. 2005, ApJ, 624, 832 Murray, N. & Chiang, J. 1997, ApJ, 474, 91 Muzerolle, J., D'Alessio, P., Calvet, N., & Hartmann, L. 2004, ApJ, 617, 406 Okazaki, A. T. 1991, PASJ, 43, 75 Petrov, R. G., Malbet, F., Weigelt, G., et al. 2007, A&A, 464, 1 Pudritz, R. E. & Banerjee, R. 2005, in IAU Symposium, Vol. 227, Massive Star Birth: A Crossroads of Astrophysics, ed. R. Cesaroni, M. Felli, E. Churchwell, & M. Walmsley, 163 Rybicki, G. B. & Hummer, D. G. 1978, ApJ, 219, 654 Safier, P. N. 1993, ApJ, 408, 115 Scholz, T. T., Walters, H. R. J., Burke, P. J., & Scott, M. P. 1990, MNRAS, 242, 692 Shlosman, I. & Vitello, P. 1993, ApJ, 409, 372 Sobolev, V. V. 1960, Moving envelopes of stars (Cambridge: Harvard University Press). Stee, P. & de Araujo, F. X. 1994, A&A, 292, 221 Stee, P., de Araujo, F. X., Vakili, F., et al. 1995, A&A, 300, 219 Tambovtseva, L. V., Grinin, V. P., Rodgers, B., & Kozlova, O. V. 2001, Astronomy Reports, 45, 442 Tannirkulam, A., Monnier, J. D., Harries, T. J., et al. 2008, ApJ, 689, 513 Tatulli, E., Millour, F., Chelli, A., et al. 2007, A&A, 464, 29 Vieira, S. L. A., Corradi, W. J. B., Alencar, S. H. P., et al. 2003, AJ, 126, 2971 Weigelt, G., Kraus, S., Driebe, T., et al. 2007, A&A, 464, 87 Zickgraf, F.-J. 2003, A&A, 408, 257

20