# Enhanced H α activity at periastron in the young and massive spectroscopic binary HD 200775

M. Benisty, K. Perraut, D. Mourard, P. Stee, G. Lima, J. Le Bouquin, M.

Borges Fernandes, O. Chesneau, N. Nardetto, I. Tallon-Bosc, et al.

► To cite this version:

M. Benisty, K. Perraut, D. Mourard, P. Stee, G. Lima, et al., Enhanced H α activity at periastron in the young and massive spectroscopic binary HD 200775. Astronomy and Astrophysics - A&A, 2013, 555, pp.A113. 10.1051/0004-6361/201219893. hal-03473883

HAL Id: hal-03473883

https://hal.science/hal-03473883

Submitted on 17 Dec 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Enhanced $H_\alpha$ activity at periastron in the young and massive spectroscopic binary HD 200775*

M. Benisty<sup>1,2</sup>, K. Perraut<sup>1</sup>, D. Mourard<sup>3</sup>, P. Stee<sup>3</sup>, G. H. R. A. Lima<sup>1</sup>, J. B. Le Bouquin<sup>1</sup>, M. Borges Fernandes<sup>4</sup>, O. Chesneau<sup>3</sup>, N. Nardetto<sup>3</sup>, I. Tallon-Bosc<sup>5,6,7,8</sup>, H. McAlister<sup>9,10</sup>, T. Ten Brummelaar<sup>10</sup>, S. Ridgway<sup>11</sup>, J. Sturmann<sup>10</sup>, L. Sturmann<sup>10</sup>, N. Turner<sup>10</sup>, C. Farrington<sup>10</sup>, and P. J. Goldfinger<sup>10</sup>

Received 26 June 2012 / Accepted 7 May 2013

## ABSTRACT

**Context.** Young close binaries clear central cavities in their surrounding circumbinary disk from which the stellar objects can still accrete material. This process takes place within the first astronomical unit and is still not well constrained because the observational evidence has been gathered, until now, only by means of spectroscopy. Theoretical models for T Tauri stars in close binaries predict a variability of the hydrogen emission lines attributable to periodic changes in the accretion rates as the secondary approaches periastron. Whether a similar scenario applies to more massive objects is unclear, and still needs to be proven observationally.

**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_\alpha$ emission. We aim to study the mechanisms that produce the $H_\alpha$ line at the AU-scale, and their dependence on binarity.

**Methods.** Combining the radial velocity measurements and astrometric data available in the literature, we determined new orbital parameters and revised the distance to 320 ± 51 pc. With the VEGA instrument on the CHARA array, we spatially and spectrally resolved the $H_\alpha$ emission of HD 200775 on a scale of a few milliarcseconds, at low and medium spectral resolutions ($R \sim 1600$ and 5000). Our observations cover a single orbital period (~3.6 years). Spectra, spectral visibilities, and differential phases have been derived. A simple analytical model of a face-on Gaussian located along the binary axis was used to analyze the interferometric observables over the spectral range.

**Results.** We observe that the $H_\alpha$ 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_\alpha$ emission (at 1 to $2\sigma$ level) and location (at 1 to $5\sigma$ level). The spatial extent of the $H_\alpha$ emission, as probed by the Gaussian FWHM, is minimum at the ascending node (0.67 ± 0.20 mas, i.e., 0.22 ± 0.06 AU), and more than doubles at the 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_\alpha$ emission is due to accretion onto the secondary. This favors a scenario in which the primary is responsible for the enhanced $H_\alpha$ activity at periastron. These findings, together with the wide $H_\alpha$ 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_\alpha$ line and estimate that it is emitted in a region larger than the one usually inferred in accretion processes. The $H_\alpha$ 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.

**Key words.** instrumentation: interferometers – stars: variables: T Tauri, Herbig Ae/Be – stars: winds, outflows – binaries: spectroscopic

## 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

<footnote n="1"><sup>1</sup> 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
e-mail: Myriam.Benisty@obs.ujf-grenoble.fr</footnote>

<footnote n="2"><sup>2</sup> Max Planck Institut für Astronomie, Königstühl 17, 69117 Heidelberg, Germany</footnote>

<footnote n="3"><sup>3</sup> Laboratoire Lagrange, UMR 7293 UNS-CNRS-OCA, Boulevard de l'Observatoire, BP 4229 F, 06304 Nice Cedex 4, France</footnote>

<footnote n="4"><sup>4</sup> Observatorio Nacional, Rua General José Cristino, 77, 20921-400, Sao Cristovao, Rio de Janeiro, Brazil</footnote>

<footnote n="5"><sup>5</sup> Université de Lyon, 69003 Lyon, France</footnote>

<footnote n="6"><sup>6</sup> Université Lyon 1, Observatoire de Lyon, 9 avenue Charles André, 69230 Saint Genis Laval, France</footnote>

<footnote n="7"><sup>7</sup> CNRS, UMR 5574, Centre de Recherche Astrophysique de Lyon, France</footnote>

<footnote n="8"><sup>8</sup> École Normale Supérieure, 69007 Lyon, France</footnote>

<footnote n="9"><sup>9</sup> Georgia State University, PO Box 3969, Atlanta GA 30302-3969, USA</footnote>

<footnote n="10"><sup>10</sup> CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA</footnote>

<footnote n="11"><sup>11</sup> National Optical Astronomy Observatory, PO Box 26732, Tucson, AZ 85726, USA</footnote>

<footnote>* Based on observations made with the VEGA/CHARA instrument.</footnote>
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 68 ± 11%, *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 con-
sisting of a spectroscopic binary (SB) at ~18 milliseconds of arc
(mas) separation (Millan-Gabet et al. 2001), and a third compan-
ion 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 re-
ported 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 15 000-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 dis-
tinguish 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 pro-
jected separation of 15.14 ± 0.70 mas and an inclination of *i* ~
65° ± 8°. They defined as primary target, the brightest near in-
frared (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 resolu-
tion (~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<sup>+120</sup><sub>-70</sub> pc.

HD 200775 was later observed during an extensive spec-
tropolarimetric campaign (Alecian et al. 2008). Two individ-
ual components were found in the photospheric line profiles.
The authors defined as the primary target the one that emits
the sharper lines, and as the secondary 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.81 ± 0.22, indicating that the star con-
sidered 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, sug-
gesting 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 18 600 K, respectively), the discrepancy between
the observational properties suggests that the two stellar compo-
nents 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 com-
ponent, 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 mea-
sured 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 ellipt-
tical shape, suggesting an inclined flared disk (i = 54.5° ± 1.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 bicon-
ical 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 pro-
files show a complicated multi-component structure, including
a doubling of the central absorption feature with a new, vari-
able, 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 in-
volved 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 ad-
vent of spectro-interferometric instruments, it has been possible
to achieve such a resolution, and to spatially and spectrally re-
solve 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
<table><caption>**Table 1.** Log of the observations. $\phi$ is the orbital phase (zero at periastron).</caption><thead><tr><th>Date</th><th>UT (h)</th><th>HJD (2 450 000+)</th><th>HA (h)</th><th>Bp (m)</th><th>PA<sub>B</sub> (°)</th><th>Spectral resolution</th><th>Calibrator (HD number)</th><th>r<sub>0</sub> (cm)</th><th>φ</th></tr></thead><tbody><tr><td>2008-07-29</td><td>10.26</td><td>4676.9271</td><td>1.86</td><td>26.7</td><td>-31.0</td><td>1600</td><td>204770</td><td>8-10</td><td>-0.012</td></tr><tr><td>2008-08-07</td><td>9.41</td><td>4685.8924</td><td>1.60</td><td>27.0</td><td>-28.0</td><td>1600</td><td>204770</td><td>7-8</td><td>-0.006</td></tr><tr><td>2008-11-24</td><td>4.13</td><td>4794.6722</td><td>—</td><td>—</td><td>—</td><td>5000</td><td>—</td><td>5</td><td>0.070</td></tr><tr><td>2009-07-26</td><td>11.86</td><td>5038.9944</td><td>—</td><td>—</td><td>—</td><td>1600</td><td>—</td><td>5-6</td><td>0.240</td></tr><tr><td>2009-08-02</td><td>8.42</td><td>5045.8507</td><td>—</td><td>—</td><td>—</td><td>1600</td><td>—</td><td>5-6</td><td>0.245</td></tr><tr><td>2009-08-26</td><td>8.46</td><td>5069.8228</td><td>1.86</td><td>26.8</td><td>-31.0</td><td>1600</td><td>197950</td><td>8</td><td>0.262</td></tr><tr><td>2010-08-01</td><td>7.00</td><td>5409.7917</td><td>-1.25</td><td>28.3</td><td>0.9</td><td>5000</td><td>197950; 204770</td><td>10-11</td><td>0.499</td></tr><tr><td>2010-09-14</td><td>6.95</td><td>5453.7833</td><td>—</td><td>—</td><td>—</td><td>5000</td><td>—</td><td>7</td><td>0.530</td></tr><tr><td>2010-09-17</td><td>8.07</td><td>5456.8361</td><td>—</td><td>—</td><td>—</td><td>5000</td><td>—</td><td>6-7</td><td>0.532</td></tr><tr><td>2011-07-28</td><td>9.61</td><td>5770.8896</td><td>—</td><td>—</td><td>—</td><td>5000</td><td>—</td><td>5</td><td>0.752</td></tr><tr><td>2011-07-29</td><td>7.77</td><td>5771.8236</td><td>—</td><td>—</td><td>—</td><td>5000</td><td>—</td><td>7</td><td>0.752</td></tr><tr><td>2011-08-28</td><td>9.63</td><td>5801.9007</td><td>3.11</td><td>25.1</td><td>-44.0</td><td>5000</td><td>197950</td><td>6-7</td><td>0.773</td></tr><tr><td>2011-10-16</td><td>5.61</td><td>5850.7340</td><td>2.31</td><td>26.2</td><td>-36.0</td><td>5000</td><td>197950</td><td>7-8</td><td>0.806</td></tr><tr><td>2011-10-21</td><td>3.62</td><td>5855.6507</td><td>—</td><td>—</td><td>—</td><td>5000</td><td>—</td><td>5-6</td><td>0.811</td></tr></tbody></table>

**Notes.** Italics indicate that only spectra were retrieved.

describe our modeling in Sect. 5. We discuss our results in Sect. 6 and conclude in Sect. 7.

## 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<sub>B</sub>) 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 (R ~ 1600, hereafter LR), as the target is at the instrumental sensitivity limit ($m_V = 7.4$). Instrumental improvements allowed us to later record data in medium spectral resolution (R ~ 5000, 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.17 ± 0.01 mas and ~0.33 ± 0.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<sup>-1</sup>) in MR, and 0.39 nm (i.e., 178 km s<sup>-1</sup>) 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 (S/N; 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_1*V_2$. The dataset was calibrated from the residual atmospheric piston and chromatic optical path difference with the model described in Mourard et al. (2009). $V_1^2$ 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_1$ and of the transfer function, we deduced the calibrated visibility $V_2$ 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_2$ and differential phases, $\Delta\phi = \phi_2 - \phi_1$, with a final spectral resolution of ~160. We considered $\Delta\phi$ 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 $N_{Hα}$ 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 $\sqrt{N_{Hα}/N_{cont}}$. 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 S/N to produce interferometric observables.

## 3. Orbital parameters

Several determinations of the binary orbit exist in the literature. However, the SB2 RVs (Alelian et al. 2008) and resolved astrometric observations (Monnier et al. 2006) have never been fitted
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Left: $H_α$ VEGA spectra obtained at R = 1600 (LR). Right: additional VEGA spectra obtained at R = 5000 (MR).</figcaption></figure>

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 $\chi^2$ 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 (520 ± 150 pc, Van Leeuwen 2007) and the individual masses are smaller by a factor of 2 than the ones derived in Alecian et al. (2008).

# 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⁻¹, FWHM ~ 450 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 (φ ~ 0.5), 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

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Equivalent widths of the $H_α$ line over the orbit for all the spectra.</figcaption></figure>

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
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Variation of the line profile measured with VEGA at R ~ 5000, 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.</figcaption></figure>

<table><caption>**Table 2.** Detection level of the V and $\Delta\phi$ change in $H_α$.</caption><thead><tr><th>Date</th><th>φ</th><th>Detection V</th><th>Detection Δφ</th></tr></thead><tbody><tr><td>2008-07-29</td><td>-0.012</td><td>2.0σ<sub>cont</sub>/3.2σ<sub>Hα</sub></td><td>1.8σ<sub>cont</sub></td></tr><tr><td>2008-08-07</td><td>-0.006</td><td>1.0σ<sub>cont</sub>/1.7σ<sub>Hα</sub></td><td>1.5σ<sub>cont</sub></td></tr><tr><td>2009-08-26</td><td>0.262</td><td>-/1.0σ<sub>Hα</sub></td><td>—</td></tr><tr><td>2010-08-01</td><td>0.499</td><td>-/-</td><td>3.4σ<sub>cont</sub></td></tr><tr><td>2011-08-28</td><td>0.773</td><td>1.5σ<sub>cont</sub>/2.1σ<sub>Hα</sub></td><td>—</td></tr><tr><td>2011-10-16</td><td>0.806</td><td>1.0σ<sub>cont</sub>/1.8σ<sub>Hα</sub></td><td>4.6σ<sub>cont</sub></td></tr></tbody></table>

**Notes.** $\sigma_{cont}$ and $\sigma_{Hα}$ are the errors in the continuum and in the line, respectively (computed as in Sect. 2.2).

is consistent with zero for all measurements, except for $\phi \sim 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 ($\Delta\phi \le 20^\circ$), 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 ($\sigma_{cont}$ and $\sigma_{Hα}$, respectively, computed as in Sect. 2.2) and are reported in Table 2.

<table><caption>**Table 3.** Best fit orbital elements and related physical parameters.</caption><thead><tr><th>Parameter</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td>t<sub>p</sub></td><td>48 962 ± 55</td><td>JD-2 400 000</td></tr><tr><td>P</td><td>1433 ± 17</td><td>days</td></tr><tr><td>a<sub>app</sub></td><td>15.9 ± 0.7</td><td>mas</td></tr><tr><td>e</td><td>0.30 ± 0.02</td><td></td></tr><tr><td>Ω</td><td>-7.7 ± 6</td><td>deg</td></tr><tr><td>ω</td><td>-136 ± 1</td><td>deg</td></tr><tr><td>i</td><td>66 ± 7</td><td>deg</td></tr><tr><td>K<sub>a</sub></td><td>19.8 ± 0.5</td><td>km s<sup>-1</sup></td></tr><tr><td>K<sub>b</sub></td><td>16.1 ± 0.9</td><td>km s<sup>-1</sup></td></tr><tr><td>v<sub>0</sub></td><td>-8.1 ± 0.3</td><td>km s<sup>-1</sup></td></tr><tr><td>a<sub>phy</sub></td><td>5.0 ± 0.5</td><td>AU</td></tr><tr><td>d</td><td>320 ± 51</td><td>pc</td></tr><tr><td>M<sub>p</sub>/M<sub>s</sub></td><td>1.23 ± 0.08</td><td></td></tr><tr><td>M<sub>p</sub>+M<sub>s</sub></td><td>9.8 ± 3.6</td><td>M<sub>☉</sub></td></tr><tr><td>M<sub>p</sub></td><td>5.37 ± 1.9</td><td>M<sub>☉</sub></td></tr><tr><td>M<sub>s</sub></td><td>4.4 ± 1.7</td><td>M<sub>☉</sub></td></tr></tbody></table>

**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
<figure><img src="image_5.png" /><figcaption>Fig. 4. VEGA spectro-interferometric measurements over the orbit.</figcaption></figure>

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.
<table><caption>**Table 4.** Orbital and best-fit parameters.</caption><thead><tr><th>Date</th><th>φ</th><th>Δα<br>(mas)</th><th>Δδ<br>(mas)</th><th>PA<sub>bin</sub><br>(°)</th><th>PA<sub>B</sub><br>(°)</th><th>FR</th><th>Θ<sub>Hα</sub><br>(mas)</th><th>Θ<sub>Hα</sub><br>(AU)</th><th>ρ<sub>Hα</sub><br>(mas)</th><th>ρ<sub>Hα</sub><br>(AU)</th></tr></thead><tbody><tr><td>2008-07-29</td><td>-0.012</td><td>-1.372</td><td>-9.387</td><td>-171.7</td><td>-31.0</td><td>0.10<sup>+0.03</sup><sub>-0.04</sub></td><td>1.91 ± 0.13</td><td>0.61 ± 0.04</td><td>-0.10 ± 0.03</td><td>-0.03 ± 0.01</td></tr><tr><td>2008-08-07</td><td>-0.006</td><td>-1.712</td><td>-8.891</td><td>-169.1</td><td>-28.0</td><td>0.07<sup>+0.05</sup><sub>-0.06</sub></td><td>1.65 ± 0.33</td><td>0.53 ± 0.11</td><td>-0.19 ± 0.03</td><td>-0.06 ± 0.03</td></tr><tr><td>2009-08-26</td><td>0.262</td><td>-3.558</td><td>16.961</td><td>-11.8</td><td>-31.0</td><td>0.01<sup>+0.05</sup><sub>-0.01</sub></td><td>0.67 ± 0.20</td><td>0.22 ± 0.06</td><td>-0.03 ± 0.04</td><td>-0.01 ± 0.01</td></tr><tr><td>2010-08-01</td><td>0.499</td><td>3.766</td><td>15.567</td><td>13.6</td><td>0.9</td><td>0.29<sup>+0.18</sup><sub>-0.11</sub></td><td>1.41 ± 0.17</td><td>0.45 ± 0.05</td><td>0.0 ± 0.10</td><td>0.0 ± 0.03</td></tr><tr><td>2011-08-28</td><td>0.773</td><td>7.066</td><td>-4.305</td><td>121.3</td><td>-44.0</td><td>0.07<sup>+0.08</sup><sub>-0.06</sub></td><td>1.95 ± 0.20</td><td>0.62 ± 0.06</td><td>0.4 ± 0.04</td><td>0.0 ± 0.04</td></tr><tr><td>2011-10-16</td><td>0.806</td><td>6.510</td><td>-6.927</td><td>136.8</td><td>-36.0</td><td>0.99<sup>+0.01</sup><sub>-0.98</sub></td><td>1.85 ± 0.20</td><td>0.59 ± 0.06</td><td>-0.32 ± 0.30</td><td>-0.10 ± 0.10</td></tr></tbody></table>

**Notes.** φ is the orbital phase, (Δα, Δδ) the angular separation, PA<sub>bin</sub> the position angle of the SB at the date of the interferometric observations. PA<sub>B</sub> is the baseline position angle, FR the binary flux ratio in the continuum, Θ<sub>Hα</sub> and ρ<sub>Hα</sub> 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.

# 5. Results

The goal of the analysis is to derive the basic properties of the H<sub>α</sub> 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 Table 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<sup>+0.05</sup><sub>-0.06</sub> to 0.29<sup>+0.18</sup><sub>-0.11</sub> (χ² = 1.5 at most). This suggests that the secondary is the dimmer object in the continuum. At the ascending node φ ~ 0.26 where V ~ 1 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 (rms ~ 10%), 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. A.1, 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 S/N observations would answer this question.

## 5.2. Size and photocenter of the H<sub>α</sub> emission

To model the H<sub>α</sub> 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 k

$$V_k e^{i\Delta\phi_k} = \frac{F_{\text{cont},k} V_{\text{cont},k} + F_{\text{H}_\alpha,k} V_{\text{H}_\alpha,k} e^{i\Delta\phi_{\text{H}_\alpha,k}}}{F_{\text{cont},k} + F_{\text{H}_\alpha,k}}, \quad (1)$$

where the indices cont and H<sub>α</sub> refer to the continuum (binary) and circumstellar emission, Δφ<sub>k</sub> is the measured differential phases in the spectral channel k, and Δφ<sub>Hα</sub> is the differential phase due to the H<sub>α</sub> emitting region. The real and imaginary parts of Eq. (1) lead to two equations that can be solved for V<sub>Hα</sub> and Δφ<sub>Hα</sub> (see, e.g., 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<sub>α</sub> absorption, we fit our model for the H<sub>α</sub> emission to V<sub>Hα</sub> and Δφ<sub>Hα</sub>, 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 2D 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<sub>α</sub> 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<sub>α</sub> 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.67 ± 0.20 mas (i.e., ~0.22 AU) and increases up to ~1.95 ± 0.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.10 ± 0.03 mas (i.e., 0.03 ± 0.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
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Best-fit continuum and $H_\alpha$ line model (dashed line). The first two lines give the visibilities, the following two, the differential phases. We recall that the visibilities for $\phi = 0.806$ have been scaled to the binary model with FR = 0.10.</figcaption></figure>

angle for the $H_\alpha$ 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_\alpha$ characteristic size along the orbit and the negative photocenter displacements still hold.

# 6. Discussion

## 6.1. Characteristic spatial extents in $H_\alpha$

It is useful to compare the $H_\alpha$ 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 ~15 000 $L_\odot$) 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).

<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Best-fit model parameters, Gaussian FWHM (top) and displacement (bottom), over the orbit.</figcaption></figure>

Considering the $v\sin(i)$ estimate and stellar properties from Alecian et al. (2008), we find that the corotation radius should be ~0.30 ± 0.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_\alpha$ 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\gamma$ 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_\alpha$ emission is due to a rotating disk, even in the quiescent state, where the $H_\alpha$ line profile shows a double peak at low velocities (~75 km s⁻¹). If we consider a mass of ~10 $M_\odot$, and the velocity field of a disk in Keplerian rotation ($v = \sqrt{GM_*/R}/\sin(i)$), the $H_\alpha$ double peak in the quiescent state could result in a rotating disk with a minimum outer radius of $R \sim 1.2$ AU, a much larger value than our estimates (Table 4). The disk region responsible for the broad $H_\alpha$ wings would be located very close to or inside the stellar surface (~7.6 $R_\odot$).

## 6.2. Origin of the $H_\alpha$ burst

The results given in the previous section can be interpreted in the context of mass ejection. The extended $H_\alpha$ 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_s 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 Figs. 2 and 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, Table 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 re-processed 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 S/N 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

<figure><img src="image_8.png" /><figcaption>**Fig. A.1.** 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 $\phi = 0.806$ were re-scaled.</figcaption></figure>

*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. G.H.R.A.L. acknowledges CAPES for the support. This research has made use of the Jean-Marie Mariotti Center SearchCal service<fnref n="1" /> co-developped by FIZEAU and IPAG, and of CDS Astronomical Databases SIMBAD and VIZIER<fnref n="2" />.

## References

Alecian, E., Catala, C., Wade, G. A., et al. 2008, MNRAS, 385, 391

Alonso-Albi, T., Fuente, A., Bachiller, R., et al. 2009, A&A, 497, 117

Artymowicz, P., & Lubow, S. H. 1996, ApJ, 467, L77

Baines, D., Oudmaijer, R. D., Porter, J. M., & Pozzo, M. 2006, MNRAS, 367, 737

Basri, G., Johns-Krull, C. M., & Mathieu, R. D. 1997, AJ, 114, 781

Benisty, M., Malbet, F., Dougados, C., et al. 2010, A&A, 517, L3

Berio, P., Mourard, D., Bonneau, D., et al. 1999, J. Opt. Soc. Am. A, 16, 872

Casse, F., & Ferreira, J. 2000, A&A, 353, 1115

de Val-Borro, M., Gahm, G. F., Stempels, H. C., & Pepliński, A. 2011, MNRAS, 413, 2679

Duchêne, G., Delgado-Donate, E., Haisch, K. E., Jr., Loinard, L., & Rodríguez, L. F. 2007, Protostars and Planets V, 379

Dullemond, C. P., & Monnier, J. D. 2010, ARA&A, 48, 205

Eisner, J. A., Monnier, J. D., Woillez, J., et al. 2010, ApJ, 718, 774

Elitzur, M., Asensio Ramos, A., & Ceccarelli, C. 2012, MNRAS, 2604

Ellerbroek, L. E., Kaper, L., Bik, A., et al. 2011, ApJ, 732, L9

Fuente, A., Martin-Pintado, J., Rodriguez-Franco, A., & Moriarty-Schieven, G. D. 1998, A&A, 339, 575

Günther, R., & Kley, W. 2002, A&A, 387, 550

Hernández, J., Calvet, N., Briceño, C., Hartmann, L., & Berlind, P. 2004, AJ, 127, 1682

Kraus, S., Hofmann, K.-H., Benisty, M., et al. 2008, A&A, 489, 1157

Kraus, S., Calvet, N., Hartmann, L., et al. 2012, ApJ, 752, 11

Le Bouquin, J.-B., Beust, H., Duvert, G., et al. 2013, A&A, 551, A121

Li, W., Evans, N. J., II, Harvey, P. M., & Colome, C. 1994, ApJ, 433, 199

Malbet, F., Benisty, M., de Wit, W.-J., et al. 2007, A&A, 464, 43

Millan-Gabet, R., Schloerb, F. P., & Traub, W. A. 2001, ApJ, 546, 358

Miroshnichenko, A. S., Mulliss, C. L., Bjorkman, K. S., et al. 1998, PASP, 110, 883

Monnier, J. D., & Millan-Gabet, R. 2002, ApJ, 579, 694

Monnier, J. D., Berger, J.-P., Millan-Gabet, R., et al. 2006, ApJ, 647, 444

Monnier, J. D., Tuthill, P. G., Ireland, M., et al. 2009, ApJ, 700, 491

Mourard, D., Clausse, J. M., Marcotto, A., et al. 2009, A&A, 508, 1073

Mourard, D., Challouf, M., Ligi, R., et al. 2012, in SPIE Conf. Ser., 8445

Mundt, R., & Ray, T. P. 1994, In The Nature and Evolutionary Status of Herbig Ae/Be Stars, eds. P. S. The, M. R. Perez, & E. P. J. van den Heuvel, ASP Conf. Ser., 62, 237

Okamoto, Y. K., Kataza, H., Honda, M., et al. 2009, ApJ, 706, 665

Pogodin, M. A., Miroshnichenko, A. S., Tarasov, A. E., et al. 2004, A&A, 417, 715

Rousselet-Perraut, K., Benisty, M., Mourard, D., et al. 2010, A&A, 516, L1

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337

Shu, F., Najita, J., Ostriker, E., et al. 1994, ApJ, 429, 781

Stee, P., & de Araujo, F. X. 1994, A&A, 292, 221

ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., et al. 2005, ApJ, 628, 453

Vaidya, B., Fendt, C., Beuther, H., & Porth, O. 2011, ApJ, 742, 56

van Boekel, R., Juhász, A., Henning, T., et al. 2010, A&A, 517, A16

Watt, G. D., Burton, W. B., Choe, S.-U., & Liszt, H. S. 1986, A&A, 163, 194

Weigelt, G., Kraus, S., Driebe, T., et al. 2007, A&A, 464, 87

Weigelt, G., Grinin, V. P., Groh, J. H., et al. 2011, A&A, 527, A103

<footnote n="1">¹ http://www.jmmc.fr/searchcal</footnote>

<footnote n="2">² http://cdsweb.u-strasbg.fr/</footnote>