Enhanced H{}_{\alpha} activity at periastron in the young and massive spectroscopic binary HD 200775 Based on observations made with the VEGA/CHARA instrument

Enhanced H activity at periastron in the young and massive spectroscopic binary HD 200775 thanks: Based on observations made with the VEGA/CHARA instrument

M. Benisty Institut d’Astrophysique et de Planétologie de Grenoble, CNRS-UJF UMR 5571, 414 rue de la Piscine, 38400 St Martin d’Hères, France Max Planck Institut für Astronomie, Konigstühl 17, 69117 Heidelberg, Germany    K. Perraut Institut d’Astrophysique et de Planétologie de Grenoble, CNRS-UJF UMR 5571, 414 rue de la Piscine, 38400 St Martin d’Hères, France    D. Mourard Laboratoire Lagrange, UMR 7293 UNS-CNRS-OCA, Boulevard de l’Observatoire, B.P. 4229 F, 06304 NICE Cedex 4, France    P. Stee Laboratoire Lagrange, UMR 7293 UNS-CNRS-OCA, Boulevard de l’Observatoire, B.P. 4229 F, 06304 NICE Cedex 4, France    G.H.R.A. Lima Institut d’Astrophysique et de Planétologie de Grenoble, CNRS-UJF UMR 5571, 414 rue de la Piscine, 38400 St Martin d’Hères, France    J.B. Le Bouquin Institut d’Astrophysique et de Planétologie de Grenoble, CNRS-UJF UMR 5571, 414 rue de la Piscine, 38400 St Martin d’Hères, France    M. Borges Fernandes Observatorio Nacional, Rua General José Cristino, 77, 20921-400, Sao Cristovao, Rio de Janeiro, Brazil    O. Chesneau Laboratoire Lagrange, UMR 7293 UNS-CNRS-OCA, Boulevard de l’Observatoire, B.P. 4229 F, 06304 NICE Cedex 4, France    N. Nardetto Laboratoire Lagrange, UMR 7293 UNS-CNRS-OCA, Boulevard de l’Observatoire, B.P. 4229 F, 06304 NICE Cedex 4, France    I. Tallon-Bosc Université de Lyon, 69003 Lyon, France; Université Lyon 1, Observatoire de Lyon, 9 avenue Charles André, 69230 Saint Genis Laval; CNRS, UMR 5574, Centre de Recherche Astrophysique de Lyon; Ecole Normale Supérieure, 69007 Lyon, France    H. McAlister Georgia State University, PO Box 3969, Atlanta GA 30302-3969, USA CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA    T. Ten Brummelaar CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA    S. Ridgway National Optical Astronomy Observatory, PO Box 26732, Tucson, AZ 85726, USA
   J. Sturmann CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA    L. Sturmann CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA    N. Turner CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA    C. Farrington CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA    P.J. Goldfinger CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA
Received 26 June 2012 / Accepted 7 May 2013
Key Words.:
Methods: observational Techniques: high angular resolution - Techniques: interferometric - Stars: binary (HD 200775) Stars: emission-line - Stars: circumstellar matter
offprints: Myriam.Benisty@obs.ujf-grenoble.fr
Abstract

Context:Young close binaries clear central cavities in their surrounding circumbinary disk from which the stars can still accrete material. This process takes place within the very first astronomical units, and is still not well constrained as the observational evidence has been gathered, until now, only by means of spectroscopy.

Aims:The young object HD 200775 (MWC 361) is a massive spectroscopic binary (separation of 15.9 mas, 5.0 AU), with uncertain classification (early/late Be), that shows a strong and variable H emission. We aim to study the mechanisms that produce the H line at the AU-scale.

Methods:Combining the radial velocity measurements and astrometric data available in the literature, we determined new orbital parameters. With the VEGA instrument on the CHARA array, we spatially and spectrally resolved the H emission of HD 200775 on a scale of a few milliarcseconds, at low and medium spectral resolutions (R1600 and 5000) over a full orbital period (3.6 years).

Results:We observe that the H equivalent width varies with the orbital phase, and increases close to periastron, as expected from theoretical models that predict an increase of the mass transfer from the circumbinary disk to the primary disk. In addition, using spectral visibilities and differential phases, we find marginal variations of the typical extent of the H emission (at 1 to 2 level) and location (at 1 to 5 level). The spatial extent of the H emission, as probed by a Gaussian FWHM, is minimum at the ascending node (0.670.20 mas, i.e., 0.220.06 AU), and more than doubles at periastron. In addition, the Gaussian photocenter is slightly displaced in the direction opposite to the secondary, ruling out the scenario in which all or most of the H emission is due to accretion onto the secondary. These findings, together with the wide H line profile, may be due to a non-spherical wind enhanced at periastron.

Conclusions:For the first time in a system of this kind, we spatially resolve the H line and estimate that it is emitted in a region larger than the one usually inferred in accretion processes. The H line could be emitted in a stellar or disk-wind, enhanced at periastron as a result of gravitational perturbation, after a period of increased mass accretion rate. Our results suggest a strong connection between accretion and ejection in these massive objects, consistent with the predictions for lower-mass close binaries.

1 Introduction

Herbig AeBe (HAeBe) stars are pre-main-sequence objects of intermediate mass, with spectral types from B to F. They are surrounded by protoplanetary disks of gas and dust, responsible for the observed excess emission from the infrared to the submillimeter (e.g., Alonso-Albi et al., 2009). Their spectra display many emission lines, which are signatures of accretion and outflow (Mundt & Ray, 1994). These objects are of particular interest as they lie between solar-mass young stars, thought to form through gravitational collapse of a molecular cloud, and the massive young stars for which the process of formation is still matter of debate.

Like their lower mass counterparts (e.g., Duchêne et al., 2007), a large fraction of Herbig AeBe stars is found in multiple systems (up to 6811%, Baines et al. (2006)). Once formed, close binary stars are thought to accrete mass from an envelope through a circumbinary disk. However, it is not clear what effect preferential mass accretion will have on the evolution of the two components. Observations have shown that a number of close T Tauri star binaries (e.g., DQ Tau, Basri et al., 1997) show enhanced emission line activity close to periastron, indicating that the accretion is non-axisymmetric. Numerical studies of young close binary systems have shown that an inner cavity forms inside the 2:1 resonance and that accretion streamers can still feed the stars inside the circumbinary disk, producing the observed periodic line changes (Artymowicz & Lubow, 1996; Günther & Kley, 2002; de Val-Borro et al., 2011). This interaction between the circumbinary disk and the stars occurs at (sub-)AU scales and, until now, only spatially unresolved observations have been published. Being able to directly probe the disk-binary interaction is crucial to test models of young binary evolution.

The object HD 200775 (MWC361) is a triple system consisting of a spectroscopic binary (SB) at 18 milliseconds of arc (mas) separation (Millan-Gabet et al., 2001), and a third companion at 6” (Li et al., 1994). Based on the analysis of the H line, a radial velocity (RV) orbit with a period of 1341 days was reported for the SB (Pogodin et al., 2004). From their analysis of 33 spectral features, Hernández et al. (2004) classified the SB as a Herbig Be star with a luminosity of 15000-L. However, because of the large uncertainty on the age of the system, the fundamental parameters of the individual sources remain highly uncertain, although it is likely that at least one of the two is an early Herbig Be star that dominates the spectrum.

Various terminologies and criteria have been used to distinguish one star from the other in the literature, generating some confusion. Using H-band interferometric measurements, Monnier et al. (2006) determined an astrometric orbit with a projected separation of 15.14  0.70 mas and an inclination of i658. They defined as primary target, the brightest near infrared (NIR) source, and found an H-band brightness ratio of 6.5  0.5. They modeled the NIR visibilities around the primary star, assuming that the secondary was unresolved in the NIR, i.e., it did not possess any extended disk at their angular resolution (4.3 mas). Using analytical models, they found a uniform disk diameter of 3.6  0.5 mas, i.e., 1.3  0.2 AU, at a distance of 360 pc.

HD 200775 was later observed during an extensive spectropolarimetric campaign (Alecian et al., 2008). Two individual components were found in the photospheric line profiles. The authors defined as the primary target the one that emits the sharper lines: the secondary was determined to be the one responsible for the broader, shallower lines. Based on the RV of the lines, the authors provided orbital parameters in overall agreement with Pogodin et al. (2004) and Monnier et al. (2006), and a mass ratio primary/secondary of 0.810.22, indicating that the star considered as secondary is in practice the most massive one. The H bisector velocities computed at a line intensity of 1.5 times the continuum were found to trace the RV of the secondary, suggesting that the line emission is dominated by this star. Their observations revealed a strongly inclined dipolar magnetic field (1000  150 G), found to have been stable for more than 2 years, and related to the object considered as the primary. Although the authors derive similar masses and effective temperatures (10 M and 18600 K, respectively), the discrepancy between the observational properties suggests that the two stellar components must have grown and evolved differently.

It seems very likely that the component defined as secondary in Alecian et al. (2008) is actually the one defined as primary in Monnier et al. (2006), the most massive component. The RV mass estimate and the dominant H activity is consistent with the presence of a circumstellar dusty disk as invoked in Monnier et al. (2006) around the most massive component. In this paper, we consider the primary to be the most massive component, i.e., the one that possesses a circumstellar dusty disk and dominates the H emission. The secondary is the least massive component that possesses a strong magnetic field.

Mid-infrared (MIR) images obtained with the Keck segment-tilting experiment indicated a large halo containing 45% of 10.7 m flux with a north-south elongation (Monnier et al., 2009) consistent with the orientation of the binary orbit measured by Monnier et al. (2006). This suggests that the halo is the remnant of a circumbinary disk. This was confirmed by Subaru MIR images that showed a diffuse emission with an elliptical shape, suggesting an inclined flared disk (i=54.51.2; Okamoto et al., 2009). The MIR emission extends up to 20 times the semi-major axis of the binary, which indicates a large gap in the system. Furthermore, the system lies in a large scale biconical cavity that has very likely been excavated by an extended bipolar outflow inclined by 70 (Fuente et al., 1998; Watt et al., 1986), a value also close to the orbital plane inclination. These results support the presence of a circumbinary disk in the same plane as the orbit.

The H and H emission lines show great variations over time, indicating changes in activity with a period of 3.68 years, in agreement with the binary period (Pogodin et al., 2004). In low states, the lines are double-peaked. In active states, the line intensities and equivalent widths (EW) increase, while their profiles show a complicated multi-component structure, including a doubling of the central absorption feature with a new, variable, blue-shifted component in addition to the pre-existing red-shifted one. The EW is found to be at its maximum right after periastron, indicating that the line activity is indeed related to the binarity.

The hot gas, responsible for the Hydrogen line, can be involved in accretion and ejection flows close to the source and can be used to probe the corresponding physical conditions. These phenomena, however, occur in a small region of a few AU around the star, corresponding to a few mas. With the recent advent of spectro-interferometric instruments, it has been possible to achieve such a resolution, and to spatially and spectrally resolve some of these lines. The first studies of the kind showed that the Br emission line was probably tracing winds or gas in a rotating disk (Malbet et al., 2007; Kraus et al., 2008). One of the challenging goals is to study the launching points of the jets and discriminate between the various theoretical models, X-wind (Shu et al., 1994) and disk-wind (Casse & Ferreira, 2000). Inspired by the magneto-centrifugal scenario for the acceleration of jets, Rousselet-Perraut et al. (2010) and Weigelt et al. (2011) provided realistic solutions to account for the line emission in disk winds.

From 2008 to 2011, we led an observing campaign over an entire orbit of HD 200775 with the optical spectro-interferometer VEGA installed at the CHARA Array. The paper is organized as follows: in Sect. 2 we describe the observations and the data processing. We present the new orbital solution, spectra, and interferometric observables in Sects. 3 and 4, and describe our modeling in Sect. 5. We discuss our results in Sect. 6 and conclude in Sect. 7.

Date UT HJD HA Bp PA Spectral Calibrator
(h) (2 450 000+) (h) (m) () resolution (HD number) (cm)
2008-07-29 10.26 4676.9271 1.86 26.7 -31.0 1600 204770 8-10 -0.012
2008-08-07 9.41 4685.8924 1.60 27.0 -28.0 1600 204770 7-8 -0.006
2008-11-24 4.13 4794.6722 5000 5 0.070
2009-07-26 11.86 5038.9944 1600 5-6 0.240
2009-08-02 8.42 5045.8507 1600 5-6 0.245
2009-08-26 8.46 5069.8228 1.86 26.8 -31.0 1600 197950 8 0.262
2010-08-01 7.00 5409.7917 -1.25 28.3 0.9 5000 197950 ; 204770 10-11 0.499
2010-09-14 6.95 5453.7833 5000 7 0.530
2010-09-17 8.07 5456.8361 5000 6-7 0.532
2011-07-28 9.61 5770.8896 5000 5 0.752
2011-07-29 7.77 5771.8236 5000 7 0.752
2011-08-28 9.63 5801.9007 3.11 25.1 -44.0 5000 197950 6-7 0.773
2011-10-16 5.61 5850.7340 2.31 26.2 -36.0 5000 197950 7-8 0.806
2011-10-21 3.62 5855.6507 5000 5-6 0.811
Table 1: Log of the observations. is the orbital phase (zero at periastron). Italics indicate that only spectra were retrieved.

2 Observations and data processing

2.1 VEGA observations

The data were collected at the CHARA array (ten Brummelaar et al., 2005), with the VEGA instrument (Mourard et al., 2009). Our datasets were obtained with two telescopes (S1S2) and cover the period from July 2008 to October 2011, i.e., an orbital period. The average projected baseline length (Bp) is about 27 m, and the average baseline position angles (PA) are given in Table 1. The angular resolution (/2Bp) of our observations is 2.5 mas. We therefore fully resolve the spectroscopic binary. The interferometric field of view is 2” in the slit direction and excludes the third object of the system (at 6”).

The first dataset was recorded with the lowest spectral resolution of VEGA (R1600, hereafter LR), as the target is at the instrumental sensitivity limit ( = 7.4). Instrumental improvements allowed us to later record data in medium spectral resolution (R5000, hereafter MR).

Each observation followed a calibrator-star-calibrator sequence, with 40 files of 1000 short exposures (15 ms) per observation. The calibrators were chosen to be close to the target both in distance and in spectral type, to be small enough at visible wavelengths, and to have an angular diameter known to an accuracy of a few percentage points. Using the SearchCal JMMC tool, we selected HD 204770 and HD 197950 (uniform disk diameters of 0.170.01 mas and 0.330.02 mas, respectively).

2.2 Data processing

Spectra: the spectra were extracted using a classical scheme of collapsing the 2D flux in one spectrum, calibrating the pixel-wavelength relation using a Thorium-Argon lamp, and normalizing the continuum by a polynomial fit. We used the H absorption lines of the calibrators to check the spectral calibration at medium resolution. The accuracy of the spectral calibration is 0.13 nm (i.e., 60 km/s) in MR, and 0.39 nm (i.e., 178 km/s) in LR.

Visibilities and differential phases: The standard routines of the VEGA data reduction pipeline were used (Mourard et al., 2009). We computed visibilities and differential phases in individual spectral channels, using the cross-spectrum method between two spectral channels [1] and [2] (for more details, see Berio et al., 1999). To reach a sufficient signal-to-noise ratio (SNR) (at least 1 photon per speckle, spectral channel and single exposure), we considered the spectral band [1] to be as wide as the entire spectral range (32 nm, i.e., a broad band measurement), and [2], to be 4 nm-wide. This method first led to the determination of the differential phase between [1] and [2], and of the product of the visibility amplitude V*V. The dataset was calibrated from the residual atmospheric piston and chromatic optical path difference with the model described in Mourard et al. (2009). were estimated using the integration of the spectral densities of the short exposures over the entire spectral range [1] and calibrated from the instrumental transfer function estimated on the calibrators. Using an estimate of V and of the transfer function, we deduced the calibrated visibility V for the narrow spectral channel [2]. By sliding the narrow spectral channel [2] with steps of 2 nm across the entire spectral range [1], we obtained a set of 13 visibilities V and differential phases, , with a final spectral resolution of 160. We considered equal to 0 in the continuum part of the spectrum, which means that the differential phases measured in the H line correspond to astrometric offsets along the baseline direction, with the photocenter in the continuum as a reference point.

The data reduction pipeline provided individual errors for each spectral channel, that account for photon noise only. To account for other sources of noise, we adopted conservative errors by considering the largest value between the rms in the continuum and the error computed by the pipeline. As the number of photons is much higher in the H spectral channel than in the continuum, we estimated the error in the H visibility by dividing the continuum error by . Finally, as HD 200775 is a faint target compared to the VEGA limiting magnitude, excellent weather conditions are required to obtain good quality data. As a consequence, only 6 observations over our 14 attempts have led to a sufficient SNR to produce interferometric observables.

Figure 1: Left: H VEGA spectra obtained at R=1600 (LR). Right: Additional VEGA spectra obtained at R=5000 (MR).
Figure 2: Equivalent widths of the H line over the orbit for all the spectra.
Figure 3: Variation of the line profile measured with VEGA at R5000, with the orbital phase. In each panel, the astrometric orbit is shown, with a full black circle indicating the corresponding orbital phase. The thick part is oriented towards the observer.
Date Detection V Detection
2008-07-29 -0.012 2.0 / 3.2 1.8
2008-08-07 -0.006 1.0 / 1.7 1.5
2009-08-26 0.262 - / 1.0 -
2010-08-01 0.499 - / - 3.4
2011-08-28 0.773 1.5 / 2.1 -
2011-10-16 0.806 1.0 / 1.8 4.6
Table 2: Detection level of the V and change in H. and are the errors in the continuum and in the line, respectively (computed as in Sect. 2.2).

3 Orbital parameters

Several determinations of the binary orbit exist in the literature. However, the SB2 radial velocities (Alecian et al., 2008) and resolved astrometric observations (Monnier et al., 2006) have never been fitted conjointly. Consequently, we revisit these works with the goal of determining the dynamical distance of the system and the individual masses of the components. We combined the available datasets following the formalism detailed in Le Bouquin et al. (2013) and performed a Levenberg-Marquardt least-square fit of the data. We are confident in our orbital determination, as the convergence toward a single and deep minimum is fast and robust for a wide range of initial guesses. The error bars on the parameters were obtained by bootstrapping. The best-fit parameters are provided in Table 3. We note that the dynamical distance is out of the confidence interval of the new reduction of the Hipparcos data ( pc, Van Leeuwen 2007) and the individual masses are smaller by a factor of 2 than the ones derived in Alecian et al. (2008).

Parameter Value Unit
JD-2400000
days
mas
km s
AU
pc
/
+
Table 3: Best fit orbital elements and related physical parameters

4 Spectro-interferometry

4.1 Spectroscopy

Figure 1 presents 14 different H spectra that show a drastic change of intensity with time as well as slight changes in the full width at half maximum (FWHM) of the line. With these measurements, we confirm that the equivalent width (EW) of the line (Fig. 2) varies with the orbital phase, and reaches its maximum close to the periastron. In the following, we refer to ’active phase or state’ when the system is close to the periastron, otherwise to ’quiescent phase or state’. In the quiescent state, we measure an EW of 35-40 Å, about twice as much as in the active phase which is in agreement with previous studies (Miroshnichenko et al., 1998; Alecian et al., 2008). For the sake of clarity, Fig. 3 shows six of the spectra obtained at medium spectral resolution, with a schematic of the system at the corresponding orbital epoch. The line is double peaked in the quiescent state with peak-velocities of 75-80 km/s, FWHM450 km/s and broad wings at very high velocities (800 km/s), as determined through Gaussian fitting of the profiles. The profile is almost symmetric close to the ascending node (), and slightly asymmetric as the system gets closer to the periastron (0.8), with more blue-shifted emission. At periastron (0), the line is single peaked, and shows a slight asymmetry with more redshifted emission. There is no variability in the broad wings. These results support the idea that the binarity is at the origin of the line profile and intensity variations, as expected and observed in other close binaries.

4.2 Interferometry

Figure 4 presents the visibilities and differential phases, as well as the corresponding spectra at a resolution of 1600.

Change in the spectral shape of the visibility. By comparing the interferometric quantities in the continuum to the ones in the spectral channel centred on H, we detect a marginal visibility drop in the line during the active states (0 and 0.8; Table 2). This means that at least at these dates, the bulk of the H emission is spatially resolved and extended.

Change in the spectral shape of the differential phases. We find that within large error bars, the differential phase in the line is consistent with zero for all measurements, except for 0.5 and 0.8. The differential phase signals can be translated into a photocenter displacement in the plane of the sky along the direction of the projected baseline. For the S1S2 baseline used for the observations, a positive phase corresponds to a photocenter displacement towards the south, and a negative phase corresponds to the displacement towards the north, as detailed in Mourard et al. (2012). Therefore, the measured non-zero differential phases trace for all epochs a photocenter displacement in the direction opposite to the secondary. As the measured differential phases are low (), the photocenter of the H emission is rather close to the continuum/binary photocenter. Detection levels are defined using the errors in the continuum and in the line ( and , respectively, computed as in Sect. 2.2) and are reported in Table 2.

Change in the absolute level of the visibility. We notice that the level of the continuum visibility varies with the orbital phase, even if the measurements were obtained at similar angular resolutions (2.5 mas). This behavior is expected since the continuum emission is due to the binary that is spatially resolved by the interferometer. Therefore, any change in binary separation, position angle and/or flux ratio strongly affects the visibilities. Close to the ascending node (2009 dataset), we measure a continuum visibility close to 1, corresponding to an unresolved emission along the baseline position angle, while in the active state, the continuum visibility can be as low as 0.76.

Date PA PA FR
(mas) (mas) () () (mas) (AU) (mas) (AU)
2008-07-29 -0.012 -1.372 -9.387 -171.7 -31.0 0.10 1.910.13 0.610.04 -0.100.03 -0.030.01
2008-08-07 -0.006 -1.712 -8.891 -169.1 -28.0 0.07 1.650.33 0.530.11 -0.190.03 -0.060.03
2009-08-26 0.262 -3.558 16.961 -11.8 -31.0 0.01 0.670.20 0.220.06 -0.030.04 -0.010.01
2010-08-01 0.499 3.766 15.567 13.6 0.9 0.29 1.410.17 0.450.05 0.00.10 0.00.03
2011-08-28 0.773 7.066 -4.305 121.3 -44.0 0.07 1.950.20 0.620.06 0.40.04 0.00.04
2011-10-16 0.806 6.510 -6.927 136.8 -36.0 0.99 1.850.20 0.590.06 -0.320.30 -0.100.10
Table 4: Orbital phase (), angular separation (, ), and position angle (PA) of the SB at the date of the interferometric observations. PA is the baseline position angle, FR the binary flux ratio in the continuum, and the H FWHM and displacements resulting from the best fit of a face-on Gaussian model, in mas and AU at 320 pc. The displacements are counted positive towards the secondary.
Figure 4: VEGA spectro-interferometric measurements over the orbit.

5 Results

The goal of the analysis is to derive the basic properties of the H emitting region. This requires a proper subtraction of the continuum emission that is only due to the binary stars. To do so, we first fit a binary model to the continuum visibilities, with fixed orbital parameters (Table 3), in order to determine the flux ratio (FR) between the two stars. This flux ratio will set our model in the continuum for both the visibilities and differential phases.

5.1 Binary flux ratio

In the spectral range of our observations, we assume that the only emitters in the continuum are the two stellar components. Their relative flux ratios, separations, and position angles have a strong impact on the visibility in the continuum. Assuming that the relative stellar fluxes vary over the orbit, we can use our continuum visibilities to retrieve the flux ratio at different epochs, assuming that the individual stars are unresolved, and using a simple analytical formula for the binary visibility.

We first compute the position angle and separation for each epoch, using the orbital parameters determined in Tab. 3. We then compute a large grid of 1000 binary models, and minimize a . Individual fits to each data set provide flux ratios from 0.07 to 0.29 (=1.5 at most). This suggests that the secondary is the dimmer object in the continuum. At the ascending node 0.26 where V1 and the binary is unresolved, all flux ratio values below 0.06 provide an equally good fit (=1.1), while for 0.81 all possible flux ratio values surprisingly provide a bad fit (=3.3). The second dataset shows a large intrinsic scatter in the continuum (rms10%), and seems to suffer from calibration problems that prevent us from determining a good absolute value for the visibility. Because we cannot determine a flux ratio for this measurement, we scale the average continuum visibilities, for this dataset only, to the binary model predictions with FR=0.10, which corresponds to the best-fit model assuming a constant flux ratio over the entire orbit (=3.8). The scaling factor is 25%, as shown in Fig. 7, together with the best fit of the binary model for each dataset. Values of flux ratios are given in Table 4. Because of the large error bars, we cannot determine whether the change in FR with the orbital period is significant and expect that further high SNR observations would answer this question.

5.2 Size and photocenter of the H emission

To model the H circumstellar emission, we consider the contribution of the binary to the measured visibilities and differential phases using the complex visibilities. Assuming that the binary is at the photocenter of the continuum emission (i.e., =0 in the continuum), we find that for each spectral channel

(1)

where the indices cont and refer to the continuum (binary) and circumstellar emission, is the measured differential phases in the spectral channel , and is the differential phase due to the H emitting region. The real and imaginary parts of Eq. 1 lead to two equations that can be solved for and (see, eg., Weigelt et al., 2007; Eisner et al., 2010).

With the best binary model, and the line to continuum ratio determined from the spectra after subtraction of the photospheric H absorption, we fit our model for the H emission to and , in the two spectral channels that contain most or all of the line. Considering the low quality of our datasets, we limit ourselves to a simplistic analytical approach with a face-on Gaussian model, which FWHM (considered as a typical size) and location can vary. We restrict the Gaussian displacement along the binary axis only, since our measurements have been obtained along a single baseline orientation at a time, which prevents us from determining a full 2-D position. The reference is taken to be the primary and displacements are counted positive towards the secondary. We vary the Gaussian FWHM from 0 to 6 mas, and since the sign of the differential phases indicates a displacement in the direction opposite to the secondary, we vary the Gaussian location from -4 to 0 mas, and minimize a to find the best model parameters.
Figure 5 shows the best-fit to the H visibilities and differential phases, and Fig. 6 presents the model parameters as they vary with the orbital phase (values are in Table 4). Our model fits the observations well, with a reduced of 2.7 (at most) in the line+continuum (1.1 in the H spectral channels only). Conservative error bars have been obtained by considering the extreme values of the visibilities and FR within their error bars.
We find that the spatial extent of the emission changes over the orbit. In the quiescent state, the Gaussian FWHM is 0.670.20 mas (i.e., 0.22 AU) and increases up to 1.950.20 mas (i.e., 0.62 AU), close to periastron. The displacements towards the opposite direction than the secondary appear to be quasi constant along the orbit, with absolute values close to 0.100.03 mas (i.e., 0.030.01 AU). The errors on the displacement are dominated by the errors on the binary flux ratio, hence are very large for 0.81.

Because our measurements only probe one direction at a time, we are not able to determine an inclination and a position angle for the H emitting region, and leave this for a detailed radiative transfer modeling combined with additional simultaneous multi-baseline interferometric observations. Thus, we would like to stress that the obtained parameter values are model-dependent and different geometries (e.g., ring, uniform disk) lead to slightly larger values of typical size, as expected with such simple prescriptions. In addition, the values of the displacement strongly depend on the binary flux ratio that set the phase values in the continuum. However, the variations of the H characteristic size along the orbit and the negative photocenter displacements still hold.

Figure 5: Best-fit continuum and H line model (dashed line). The first two lines give the visibilities, the following two, the differential phases. We recall that the visibilities for =0.806 have been scaled to the binary model with FR=0.10.
Figure 6: Best-fit model parameters, Gaussian FWHM (top) and displacement (bottom), over the orbit.

6 Discussion

6.1 Characteristic spatial extents in H 

It is useful to compare the H size estimates to typical radii in the close circumstellar regions probed by our measurements. If we consider the extreme values of luminosities found in the literature (3000 and 15000 L from Alecian et al. (2008); Hernández et al. (2004)), we derive an estimate of the location of the dust sublimation radius of the circumprimary disk to be between 3.7 and 8.4 AU (e.g., Dullemond & Monnier, 2010). These estimates are larger than the dust inner radius (1.73 mas, i.e., 0.55 AU at 320 pc) measured by Monnier et al. (2006). This apparent discrepancy can be explained with optically thick material inside the dust sublimation radius that efficiently shields the dust from the stellar light (Monnier & Millan-Gabet, 2002).

Considering the vsin(i) estimate and stellar properties from Alecian et al. (2008), we find that the corotation radius should be 0.300.15 AU. This estimate suffers from the large uncertainties on the stellar parameters inferred in the literature. Nonetheless, our Gaussian fit seems to indicate that the bulk of the H emission as described by a Gaussian model is most likely located between the corotation radius and the dust sublimation radius, and is more extended than the typical regions involved in accretion processes onto the star. In the case of a low mass star, the truncation radius is set by the interaction of the stellar magnetic field and the gaseous disk in a region very close to the star. In a massive star, although this is still a matter of debate, the accretion is thought to occur through a boundary layer.

Similar extents for Hydrogen line emission have been found in other Herbig Be stars, in the Br line, and interpreted as originating in a stellar or disk wind (e.g., Kraus et al., 2012; Weigelt et al., 2011; Benisty et al., 2010), in contrast with the findings of Eisner et al. (2010) on a survey of lower mass young stars that are consistent with accretion.

It also seems unlikely that the bulk of the H emission is due to a rotating disk, even in the quiescent state, where the H line profile shows a double peak at low velocities (75 km/s). If we consider a mass of 10 M, and the velocity field of a disk in Keplerian rotation (/sin(i)), the H double peak in the quiescent state could result in a rotating disk with a minimum outer radius of R1.2 AU, a much larger value than our estimates (Table 4). The disk region responsible for the broad H wings would be located very close to or inside the stellar surface (7.6 R).

6.2 Origin of the H burst

The results given in the previous section can be interpreted in the context of mass ejection. The extended H emission together with the increase of the H line intensity may trace a period of enhanced mass loss in a strong wind emitted by the primary. Such an event could follow a period of enhanced mass accretion from the circumbinary disk to the primary disk, triggered at periastron as predicted by numerical models. The photocenter displacements in the direction opposite to the secondary indicate that if a wind is responsible for the H emission, the emission should be enhanced along this direction. The wind therefore could not be fully spherical, as the consequent emission would be symmetric.

Figure 1 shows that the H broad line wings (800 km/s) are very stable, and can be due to outflowing matter, as such velocities are often seen in wind tracers. The additional emission at periastron is produced only at very low projected velocities, probably coming from a plane close to perpendicular to our line of sight. The lack of P Cygni absorption indicates that the wind is not self-absorbing along our line of sight. This would support a wind with a higher density towards low latitudes, close to the equatorial plane, in agreement with the model predictions of poorly collimated winds in massive stars (Vaidya et al., 2011). With a non-zero inclination, an outflow can also produce the double peaked profile observed in quiescent state, as shown in Stee & de Araujo (1994) for Be stars, and more recently in Weigelt et al. (2011) for a Herbig Be star of similar mass. However, complex line opacity effects can come into play, as suggested by Elitzur et al. (2012) who showed that it is in practice impossible to disentangle effects of kinematics and line opacity effects in double peaked lines and that proper radiative transfer is needed. In addition, we note that spectra obtained at high spectral resolution (Alecian et al., 2008) show that the line profile is extremely complex with variable features, which seem to indicate a combination of various variable mechanisms.

The additional emission at periastron could, in principle, come from either the primary, the secondary, or a combination of both. However, the low differential phases indicate that the photocenter of the H emitting region is close to the continuum photocenter, which is located close to the primary (since the binary best-fit flux ratio is 0.10). The displacement values found by our simple analytical model fitting therefore rule out a scenario in which all of the additional H emission at periastron comes from the secondary. It seems therefore very unlikely that accretion onto the secondary or wind-enhancement on the secondary is the key-mechanism at periastron or that both objects are equally responsible for the increase. On the contrary, because the displacements are in the opposite direction to the secondary, the burst of activity at periastron seems to be connected to the primary. Although the variation of photocenter with the orbital phase is marginal, a smaller displacement at the ascending node is consistent with a smaller emitting region and a lower line intensity. An almost-constant displacement at all other phases is also in agreement with the derivation of a similar extent at these dates. This is supported by recent numerical simulations of accretion flows in close binaries done at high resolution (de Val-Borro et al., 2011). These models suggested that accretion from the circumbinary disk onto the primary is favored at periastron. In that scenario, it is expected that most or all of the additional H emission in the active state comes from the primary star, as our data seem to indicate. We note, however, that with our single-baseline measurements, we are not able to disentangle more complex scenarios.

Whatever the mechanism favored to enhance the H emission at periastron, additional mass has to be fed into the innermost regions at these active stages. If mass is transferred through accretion streamers from the circumbinary disk to the primary disk, there should be a consequent time delay to allow matter to flow through the primary disk and reach the innermost AU. The viscous time is proportional to R/ ( c H/R), with R the radius, H the disk scale height, and the turbulence parameter (proportional to the viscosity; Shakura & Sunyaev (1973)). Considering a maximum outer disk radius for the circumprimary disk of 2 AU (because the binary separation at periastron is 4 AU, (Monnier et al., 2006)), =0.01, and H/R=0.05, the viscous time needed for the mass flow to propagate is of the order of a few ten thousand years, much larger than the binary period. An alternative scenario is that the secondary induces gravitational instabilities in the primary disk as it approaches. In fact, tidal effects could induce strong waves all over the disk spatial range and produce instabilities, which would in turn enhance the mass accretion in the inner AU and result in a higher mass loss. The exact physical process that leads to the H increase still remains to be understood.

The comparison of Fig. 2 and Fig. 6 seems to suggest that the variations in H EW and sizes are slightly shifted in time, as the size appears to increase faster than the absolute EW. This delay is difficult to interpret as we would expect a larger emitting area to produce a stronger H emission, but considering the large error bars on our size estimate and the sparse time sampling of the orbit, we cannot determine whether it is a real feature or not.

7 Summary and conclusion

In this paper, we report the first spatially and spectrally resolved observations of the H line emitting region in a close binary system using the VEGA spectrometer on the CHARA array. As previous authors have already pointed out (Pogodin et al., 2004; Alecian et al., 2008), we find that HD 200775 shows a variation of the H EW with the orbital phase. We use a simple Gaussian model located along the binary axis to fit the H visibilities and differential phases, and find marginal evidence for a change in spatial extent during the orbit that is minimum in quiescent state (close to the ascending node). The H Gaussian photocenter is at all epochs located in the direction opposite to the secondary, supporting a dominant emission by the primary.

We find typical H extents within the sublimation radius. These H spatial extents do not support an accretion origin for the bulk of the line emission, as it would involve smaller scales. Instead, we interpret these results in the context of an enhanced mass-loss event, triggered by the gravitational influence of the secondary. Such a scenario would suggest a strong connection between accretion and ejection in this Herbig Be star, as already noted in other massive objects (e.g., Ellerbroek et al., 2011; Benisty et al., 2010).

A dominant contribution to the H line by the primary is supported by our results and by the analysis of Alecian et al. (2008) and Monnier et al. (2006). However, although the stellar parameters are still highly uncertain, it is puzzling that only the primary has a circumstellar disk which accretes from the circumbinary disk. In particular, the kilo-Gauss magnetic field detected on the secondary should induce the formation of a gaseous disk which should in turn emit H. The circumstellar environments of the two objects must either have evolved differently, or the secondary could itself be an unknown spectroscopic binary of lower mass which would more easily explain the strong magnetic field, an uncommon feature in Herbig Be stars.

We find that the binary flux ratio in the continuum is 0.10 (from 0.07 to 0.29, Tab. 4). Because of large error bars on the visibilities, we are not able to retrieve a precise flux ratio at each orbital phase and to determine whether the change in H EW is correlated with some change in the amount of reprocessed and scattered light as observed in other young close binaries (van Boekel et al., 2010). We notice a discrepancy between our flux ratio estimates in the continuum and the estimate of Alecian et al. (2008), based on an analysis of spectral lines (0.67). These estimates are difficult to reconcile, unless one invokes variable extinction on one of the two objects, or if one of the two (most likely the magnetic secondary) is already chemically peculiar. Besides, the large uncertainties on the individual masses, the distance, and age of the system prevents the nature of the individual objects from being assessed precisely. As noted by Alecian et al. (2008), photometric observations of each component separately are necessary. In addition, a regular and full coverage of the orbit is required to answer a number of open questions, such as a possible time delay between periastron passage and maximum activity. More detailed modeling, as well as an intensive high SNR spectro-interferometric campaign with additional baselines at various orientations are also expected to properly constrain the complex interplay between accretion and ejection in such massive close binary systems.

Appendix A Appendix

Figure 7: Best-fit binary models in the continuum (dashed blue line). The red dotted line in the lower right panel is the prediction of a binary model with flux ratio of FR=0.10 to which the continuum visibilities measured for =0.806 were re-scaled.
Acknowledgements.
We acknowledge the anonymous referee for the constructive report that helped improve the manuscript. We thank E. Alecian, J. Bouvier, P. Garcia, A. Mérand, J.L. Monin and J. Monnier for fruitful discussions. We acknowledge all the great observers at the CHARA Array who helped to gather this large dataset. The CHARA Array is operated with support from the National Science Foundation through Grant AST-0908253, the W. M. Keck Foundation, the NASA Exoplanet Science Institute, and from Georgia State University. GHRAL acknowledges CAPES for the support. This research has made use of the Jean-Marie Mariotti Center SearchCal service 111Available at http://www.jmmc.fr/searchcal co-developped by FIZEAU and IPAG, and of CDS Astronomical Databases SIMBAD and VIZIER 222Available at http://cdsweb.u-strasbg.fr/.

References

  • Alecian et al. (2008) Alecian, E., Catala, C., Wade, G. A., et al. 2008, MNRAS, 385, 391
  • Alonso-Albi et al. (2009) Alonso-Albi, T., Fuente, A., Bachiller, R., et al. 2009, A&A, 497, 117
  • Artymowicz & Lubow (1996) Artymowicz, P. & Lubow, S. H. 1996, ApJ, 467, L77
  • Baines et al. (2006) Baines, D., Oudmaijer, R. D., Porter, J. M., & Pozzo, M. 2006, MNRAS, 367, 737
  • Basri et al. (1997) Basri, G., Johns-Krull, C. M., & Mathieu, R. D. 1997, AJ, 114, 781
  • Benisty et al. (2010) Benisty, M., Malbet, F., Dougados, C., et al. 2010, A&A, 517, L3
  • Berio et al. (1999) Berio, P., Mourard, D., Bonneau, D., et al. 1999, Journal of the Optical Society of America A, 16, 872
  • Casse & Ferreira (2000) Casse, F. & Ferreira, J. 2000, A&A, 353, 1115
  • de Val-Borro et al. (2011) de Val-Borro, M., Gahm, G. F., Stempels, H. C., & Pepliński, A. 2011, MNRAS, 413, 2679
  • Duchêne et al. (2007) Duchêne, G., Delgado-Donate, E., Haisch, Jr., K. E., Loinard, L., & Rodríguez, L. F. 2007, Protostars and Planets V, 379
  • Dullemond & Monnier (2010) Dullemond, C. P. & Monnier, J. D. 2010, ARA&A, 48, 205
  • Eisner et al. (2010) Eisner, J. A., Monnier, J. D., Woillez, J., et al. 2010, ApJ, 718, 774
  • Elitzur et al. (2012) Elitzur, M., Asensio Ramos, A., & Ceccarelli, C. 2012, MNRAS, 2604
  • Ellerbroek et al. (2011) Ellerbroek, L. E., Kaper, L., Bik, A., et al. 2011, ApJ, 732, L9
  • Fuente et al. (1998) Fuente, A., Martin-Pintado, J., Rodriguez-Franco, A., & Moriarty-Schieven, G. D. 1998, A&A, 339, 575
  • Günther & Kley (2002) Günther, R. & Kley, W. 2002, A&A, 387, 550
  • Hernández et al. (2004) Hernández, J., Calvet, N., Briceño, C., Hartmann, L., & Berlind, P. 2004, AJ, 127, 1682
  • Kraus et al. (2012) Kraus, S., Calvet, N., Hartmann, L., et al. 2012, ArXiv e-prints
  • Kraus et al. (2008) Kraus, S., Hofmann, K.-H., Benisty, M., et al. 2008, A&A, 489, 1157
  • Le Bouquin et al. (2013) Le Bouquin, J.-B., Beust, H., Duvert, G., et al. 2013, A&A, 551, A121
  • Li et al. (1994) Li, W., Evans, II, N. J., Harvey, P. M., & Colome, C. 1994, ApJ, 433, 199
  • Malbet et al. (2007) Malbet, F., Benisty, M., de Wit, W.-J., et al. 2007, A&A, 464, 43
  • Millan-Gabet et al. (2001) Millan-Gabet, R., Schloerb, F. P., & Traub, W. A. 2001, ApJ, 546, 358
  • Miroshnichenko et al. (1998) Miroshnichenko, A. S., Mulliss, C. L., Bjorkman, K. S., et al. 1998, PASP, 110, 883
  • Monnier et al. (2006) Monnier, J. D., Berger, J.-P., Millan-Gabet, R., et al. 2006, ApJ, 647, 444
  • Monnier & Millan-Gabet (2002) Monnier, J. D. & Millan-Gabet, R. 2002, ApJ, 579, 694
  • Monnier et al. (2009) Monnier, J. D., Tuthill, P. G., Ireland, M., et al. 2009, ApJ, 700, 491
  • Mourard et al. (2012) Mourard, D., Challouf, M., Ligi, R., et al. 2012, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 8445, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
  • Mourard et al. (2009) Mourard, D., Clausse, J. M., Marcotto, A., et al. 2009, A&A, 508, 1073
  • Mundt & Ray (1994) Mundt, R. & Ray, T. P. 1994, in Astronomical Society of the Pacific Conference Series, Vol. 62, The Nature and Evolutionary Status of Herbig Ae/Be Stars, ed. P. S. The, M. R. Perez, & E. P. J. van den Heuvel, 237
  • Okamoto et al. (2009) Okamoto, Y. K., Kataza, H., Honda, M., et al. 2009, ApJ, 706, 665
  • Pogodin et al. (2004) Pogodin, M. A., Miroshnichenko, A. S., Tarasov, A. E., et al. 2004, A&A, 417, 715
  • Rousselet-Perraut et al. (2010) Rousselet-Perraut, K., Benisty, M., Mourard, D., et al. 2010, A&A, 516, L1
  • Shakura & Sunyaev (1973) Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337
  • Shu et al. (1994) Shu, F., Najita, J., Ostriker, E., et al. 1994, ApJ, 429, 781
  • Stee & de Araujo (1994) Stee, P. & de Araujo, F. X. 1994, A&A, 292, 221
  • ten Brummelaar et al. (2005) ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., et al. 2005, ApJ, 628, 453
  • Vaidya et al. (2011) Vaidya, B., Fendt, C., Beuther, H., & Porth, O. 2011, ApJ, 742, 56
  • van Boekel et al. (2010) van Boekel, R., Juhász, A., Henning, T., et al. 2010, A&A, 517, A16
  • Watt et al. (1986) Watt, G. D., Burton, W. B., Choe, S.-U., & Liszt, H. S. 1986, A&A, 163, 194
  • Weigelt et al. (2011) Weigelt, G., Grinin, V. P., Groh, J. H., et al. 2011, A&A, 527, A103
  • Weigelt et al. (2007) Weigelt, G., Kraus, S., Driebe, T., et al. 2007, A&A, 464, 87
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
246522
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description