# The peculiar fast-rotating star 51 Ophiuchi probed by VEGA/CHARA

N. Jamialahmadi, P. Berio, A. Meilland, K. Perraut, D. Mourard, B. Lopez, P. Stee, N. Nardetto, B. Pichon, J. Clausse, et al.

► To cite this version:

N. Jamialahmadi, P. Berio, A. Meilland, K. Perraut, D. Mourard, et al. The peculiar fast-rotating star 51 Ophiuchi probed by VEGA/CHARA. Astronomy and Astrophysics - A&A, 2015, 579, pp.A81. 10.1051/0004-6361/201425473. [hal-03473862](https://hal.science/hal-03473862)

HAL Id: hal-03473862

https://hal.science/hal-03473862

Submitted on 13 Jan 2022

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
# The peculiar fast-rotating star 51 Ophiuchi probed by VEGA/CHARA

N. Jamiahmadi<fnref n="1" />,<fnref n="2" />, P. Berio<fnref n="2" />, A. Meilland<fnref n="2" />, K. Perraut<fnref n="3" />, D. Mourard<fnref n="2" />, B. Lopez<fnref n="2" />, P. Stee<fnref n="2" />, N. Nardetto<fnref n="2" />, B. Pichon<fnref n="2" />, J. M. Clausse<fnref n="2" />, A. Spang<fnref n="2" />, H. McAlister<fnref n="4" />,<fnref n="5" />, T. ten Brummelaar<fnref n="5" />, J. Sturmann<fnref n="5" />, L. Sturmann<fnref n="5" />, N. Turner<fnref n="5" />, C. Farrington<fnref n="5" />, N. Vargas<fnref n="5" />, and N. Scott<fnref n="5" />

Received 4 December 2014 / Accepted 19 May 2015

## ABSTRACT

*Context.* Stellar rotation is a key in our understanding of both mass-loss and evolution of intermediate and massive stars. It can lead to anisotropic mass-loss in the form of radiative wind or an excretion disk.

*Aims.* We wished to spatially resolve the photosphere and gaseous environment of 51 Oph, a peculiar star with a very high $v \sin i$ of 267 km s⁻<fnref n="1" /> and an evolutionary status that remains unsettled. It has been classified by different authors as a Herbig, a β Pic, or a classical Be star.

*Methods.* We used the VEGA visible beam combiner installed on the CHARA array that reaches a submilliarcsecond resolution. Observation were centered on the Hα emission line.

*Results.* We derived, for the first time, the extension and flattening of 51 Oph photosphere. We found a major axis of $\theta_{eq} = 8.08 \pm 0.70 R_\odot$ and a minor axis of $\theta_{pol} = 5.66 \pm 0.23 R_\odot$. This high photosphere distortion shows that the star is rotating close to its critical velocity. Finally, using spectro-interferometric measurements in the Hα line, we constrained the circumstellar environment geometry and kinematics and showed that the emission is produced in a $5.2 \pm 2 R_\star$ disk in Keplerian rotation.

*Conclusions.* From the visible point of view, 51 Oph presents all the features of a classical Be star: near critical-rotation and double-peaked Hα line in emission produced in a gaseous disk in Keplerian rotation. However, this does not explain the presence of dust as seen in the mid-infrared and millimeter spectra, and the evolutionary status of 51 Oph remains unsettled.

**Key words.** stars: emission-line, Be – stars: evolution – stars: rotation

## 1. Introduction

Many hot stars are observed to be rotating with equatorial velocities higher than 120 km s⁻<fnref n="1" /> (Abt & Morrell 1995; Abt et al. 2002). This fast stellar rotation can have strong effects on the observed stellar shape and the intensity distribution (gravity darkening). Some hot stars show emission lines in their spectra. These features indicate that a gaseous circumstellar environment surrounds these stars.

The study of the stellar photosphere and circumstellar environment of fast-rotating stars is limited by the lack of spatial resolution of single telescopes. An important and reliable way to extract this information is through long-baseline optical interferometry, which allows us to study the detailed stellar surface properties and the close-by environment through emission lines. Several rapid rotators have been studied using these techniques, including Altair, Vega, Achernar, Alderamin, Regulus, and Rasalhague (van Belle et al. 2001; Aufdenberg et al. 2006; Peterson et al. 2006; Domiciano de Souza et al. 2003; Monnier et al. 2007, respectively).

This paper is dedicated to the study of 51 Oph: a rapidly rotating star, $v \sin i = 267 \pm 5$ km s⁻<fnref n="1" /> (Dunkin et al. 1997a) with an age of 0.3 Myr and a mass of ~4 $M_\odot$ (Van den Ancker et al. 1998). This star appears to be a peculiar source in an

unusual transitional state. Compared with Herbig Ae/Be stars, van den Ancker et al. (2001) and Meeus et al. (2001) confirmed by analyzing ISO spectra that 51 Oph has more gas than other HAeBe stars. In addition, Malfait et al. (1998) and Leinert et al. (2004) found that the dusty disk of 51 Oph is much smaller (<2 AU) than the disk of most HAeBe stars.

Compared with β Pic-like stars, 51 Oph does not exhibit a far-infrared-excess bump, which is inconsistent with the presence of an outer dusty disk (Malfait et al. 1998). Thi et al. (2013) suggested a difference between 51 Oph and β Pic because the 51 Oph dust grains can be primordial, while the dust grains in an exozodiacal ring around β Pic have to be replenished. Furthermore, 51 Oph does exhibit an Hα line in emission, which is not the case of β Pic stars. Berthoud (2008) and Tatulli et al. (2008) concluded that 51 Oph is most likely a classical Be star surrounded by a compact gaseous disk. This scenario agrees well with the following properties reported in previous studies: high rotational velocity, most of the near-infrared continuum emission arising from the gaseous disk, and hot CO gas in a rotating Keplerian disk. However, the remaining presence of dust around 51 Oph is not supported by this scenario and suggests that this star follows an intriguing evolution scheme.

We here report new elements that lend credit to the hypothesis that 51 Oph is a Be star. Hereafter, we present observations of

<footnote n="1">¹ European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Munich, Germany
e-mail: njamiala@eso.org</footnote>

<footnote n="2">² Laboratoire Lagrange, UMR 7293, Univ. Nice Sophia-Antipolis, CNRS, Observatoire de la Côte d'Azur, 06304 Nice, France</footnote>

<footnote n="3">³ 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</footnote>

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

<footnote n="5">⁵ CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA</footnote>
<table><caption>**Table 1.** CHARA/VEGA observing log of 51 Oph.</caption><thead><tr><th colspan="2">Continuum</th><th colspan="2">H<sub>α</sub> line</th></tr><tr><th>Obs. time<br>(UTC)</th><th>Telescopes<br>(conf.)</th><th>Obs. time<br>(UTC)</th><th>Telescopes<br>(conf.)</th></tr></thead><tbody><tr><td>2014-07-04 05:40</td><td>E2-S2-W2</td><td>—</td><td>—</td></tr><tr><td>2014-07-04 06:15</td><td>E2-S2-W2</td><td>2013-05-25 07:12</td><td>E1-E2</td></tr><tr><td>2014-07-04 06:49</td><td>E2-S2-W2</td><td>2013-07-28 03:39</td><td>S1-S2 (1)</td></tr><tr><td>2014-07-08 05:53</td><td>E2-S2-W2</td><td>2013-07-28 05:43</td><td>S1-S2 (2)</td></tr><tr><td>2014-07-08 06:18</td><td>E2-S2-W2</td><td>2014-05-02 10:04</td><td>W1-W2 (1)</td></tr><tr><td>2014-07-08 06:43</td><td>E2-S2-W2</td><td>2014-05-02 10:47</td><td>W1-W2 (2)</td></tr></tbody></table>

**Notes.** For the continuum observations, the calibrators HD 150366 with the uniform disk angular diameter $\theta_{UD} = 0.28 \pm 0.02$ mas and $\theta_{UD} = 0.39 \pm 0.03$ mas for HD 163955 were used. For the observations in the Hα line, HD 170296 with $\theta_{UD} = 0.43 \pm 0.03$ mas was used. The predicted uniform disk angular diameter (in mas) in R band were derived from the JMMC SearchCal software (Bonneau et al. 2006).

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** (u, v) coverage obtained for the observations of 51 Oph in the continuum (left) and in the Hα line (right).</figcaption></figure>

51 Oph with a spectro-interferometer operating at optical wavelengths, the VEGA instrument (Mourard et al. 2009, 2011) installed at the CHARA Array (ten Brummelaar et al. 2005). This unique combination of spectrally and spatially resolved information allowed us to resolve the photosphere of the star and to probe the kinematics of the gaseous component producing the Hα line on sub-AU scales.

## 2. Observations and data processes

Observations were performed around the Hα line, between 630 and 670 nm, using the VEGA medium spectral resolution (MR) ($R = 5000$) to study the gaseous disk by combining three pairs of telescopes: W1W2, E1E2, and S1S2. The observation log is given in Table 1. The source 51 Oph was also observed in the continuum at 700 nm using the triplet E2S2W2 to study its photosphere. The (u, v) plane coverage corresponding to these data is presented in the left panel of Fig. 1 in continuum and in the right panel in the Hα line.

The data processing of VEGA is composed of two parts. The data in the continuum are processed using the power spectral method giving the squared visibilities (see Mourard et al. 2009, for details). For this analysis, we used a spectral bandwidth of 15 nm. The processing of the data in the spectral line is based on the cross-spectrum method, which provides differential visibilities and phases across the line. This method was applied with a spectral bandwidth of 0.2 nm.

## 3. Continuum emission: stellar photosphere

We used the LITpro2 model-fitting software for optical-to-infrared interferometric data developed by the Jean-Marie

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Normalized residuals between the continuum data and an elongated disk model (top) and a uniform disk model (bottom) are represented for three pairs of telescopes. The black points represent the residual average for each baseline.</figcaption></figure>

<table><caption>**Table 2.** Best-fit parameters obtained from the elliptical model for the stellar photosphere. $\chi_r^2$ is the total reduced $\chi^2$.</caption><thead><tr><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td>Major axis of θ<sub>eq</sub></td><td>0.6±0.05 mas</td></tr><tr><td>Minor axis of θ<sub>pol</sub></td><td>0.42±0.01 mas</td></tr><tr><td>Elongated ratio (θ<sub>eq</sub>/θ<sub>pol</sub>)</td><td>1.45 ± 0.12</td></tr><tr><td>Position angle (PA)</td><td>138 ± 3.9°</td></tr><tr><td>χ<sub>r</sub><sup>2</sup></td><td>1.52</td></tr></tbody></table>

Mariotti Center (JMMC) to interpret our data. Figure 2 shows that the elongated disk model is more consistent with the data than the uniform disk model for the stellar photosphere for the baselines S2-W2 and E2-S2: the mean of the residuals of S2-W2 and E2-S2 is not equal to zero for the uniform disk model. However, both models have nearly the same normalized residuals for the baseline W2-E2. According to Table 2, we found the major axis, $\theta_{eq}$, the minor axis, $\theta_{pol}$, the elongated ratio, $\theta_{eq}/\theta_{pol}$, and the major axis orientation, PA, using an elliptical model. Using HIPPARCOS 2007 data (van Leeuwen 2007) and calculating the distance from the parallax ($d = 124 \pm 4$ pc), we derived an equatorial radius $R_{eq} = 8.08 \pm 0.7$ R⊙ and a polar radius $R_{pol} = 5.66 \pm 0.23$ R⊙ from $\theta_{eq}$ and $\theta_{pol}$, respectively. We assumed the continuum flux to exclusively come from the central star. A comparison of the Hα emission line of 51 Oph with classical Be stars shows that the intensity level of 51 Oph is low (only 1.2 time above the continuum, see Fig. 3). The flux ratio between the disk and the central star is therefore lower than in classical Be stars. In classical Be stars, this flux ratio in the continuum is between 5% and 20% in the visible. We therefore assumed that the disk contribution in the continuum is negligible. The fit of an elliptical model for the stellar photosphere produces a reduced $\chi^2$ ($\chi_r^2$) of 1.52.
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Visibility and differential phase variations in the Hα emission line and in the surrounding continuum obtained from our medium spectral-resolution (MR) CHARA/VEGA observations (red). The best-fit kinematic model for this dataset is overplotted in blue. Right bottom panel: Hα intensity profile. Top right panel: elliptical model for the central star 51 Oph and our kinematic model for a purely rotating disk.</figcaption></figure>

# 4. Hα emission line

## 4.1. Qualitative analysis

The differential visibility and phase for the five baselines are plotted in Fig. 3. The visibility exhibits a drop in the emission line that is caused by a variation of the circumstellar environment extension and relative flux between the continuum and the line, as already explained in Meilland et al. (2007). The phase shows S-shaped variations in the line. This is characteristic of a rotating equatorial disk, as described in Meilland et al. (2012). For the longest baseline, W1W2 the disk is probably over-resolved, so that the differential phase deviates from the S shape, as explained in (Meilland et al. 2012). The VEGA spectrum of 51 Oph obtained during our observations exhibits an Hα emission line with a double-peak similar to the spectrum found by Dunkin et al. (1997b), which reflects the presence of the circumstellar disk. Hα profile wings clearly show that the disk is seen nearly edge-on.

## 4.2. Kinematic model

To model the wavelength dependence of the visibility, the differential phase and the intensity profile of Hα line emission, we used a simple kinematic model of an expanding and/or rotating thin equatorial disk<fnref n="1" />. This model has previously been used for some classical Be stars, for instance, by Meilland et al. (2011, 2012) and is described in detail by Delaa et al. (2011). Assuming no expansion, the model free-parameters can be classified into three categories:

1. the geometric parameters: stellar radius ($R_★$), distance ($d$), inclination angle ($i$), and disk major-axis position angle (PA),

2. the kinematic parameters: rotational velocity ($V_{rot}$) at the disk inner radius (i.e., photosphere) and the rotation ($\beta$) velocity laws, and

3. the disk emission line parameters: disk FWHM in the line (a) and the line equivalent width (EW).

<table><caption>**Table 3.** Best-fit parameters obtained from our simple kinematic model of the circumstellar disk.</caption><thead><tr><th>Parameter</th><th>Value</th><th>Remarks</th></tr></thead><tbody><tr><td colspan="3"><em>Geometric parameters</em></td></tr><tr><td>R<sub>eq</sub><sup>*</sup></td><td>8.08 R<sub>☉</sub></td><td>This work (Sect. 4)</td></tr><tr><td>d<sup>*</sup></td><td>124 pc</td><td>(van Leeuwen 2007)</td></tr><tr><td>i</td><td>80<sup>+7</sup><sub>-30</sub> deg</td><td>This work</td></tr><tr><td>PA</td><td>122<sup>+15</sup><sub>-25</sub> deg</td><td>This work</td></tr><tr><td colspan="3"><em>Kinematic parameter</em></td></tr><tr><td>V<sub>rot</sub><sup>*</sup></td><td>267 km s<sup>-1</sup></td><td>(Dunkin et al. 1997a)</td></tr><tr><td colspan="3"><em>H<sub>α</sub> line parameters</em></td></tr><tr><td>a</td><td>5.6 ± 2 R<sub>★</sub></td><td>This work</td></tr><tr><td>EW</td><td>5.0<sup>+2</sup><sub>-3</sub> Å</td><td>This work</td></tr><tr><td>χ<sup>2</sup><sub>r</sub></td><td>0.7</td><td></td></tr></tbody></table>

**Notes.** (*) Fixed parameters.

The star distance was calculated following van Leeuwen (2007). $R_{eq}$, the equator stellar radius, was derived from the fit of the elliptical model according to Sect. 3 (in solar radius). We assumed that the disk is directly connected to the stellar surface, therefore the rotational velocity ($V_{rot}$) should be equal to the stellar rotational velocity. We also assumed that the gaseous disk surrounding 51 Oph is Keplerian, so that we set $\beta$, the rotation velocity law, to -0.5. The other four parameters were left free.

For our target we computed several hundreds of models to constrain the free parameters, determined the uncertainties, and tried to detect any degeneracy or linked parameters. To reduce the number of computed models, we started with a qualitative estimation of the parameters from our interferometric data (PA, i, a, and EW) and explored the parameter space with decreasing steps to converge to the $\chi^2$ minimum ($\chi^2_r$). The final parameter values for the best-fit models are presented in Table 3. The corresponding differential visibilities and phases are overplotted in Fig. 3. The final reduced $\chi^2$ is 0.7.

<footnote n="1">¹ As input of our model, we used a synthetic absorption in $H_α$ line calculated with the TLUSTY code (Hubeny 1998).</footnote>
The inclination angle found for the gaseous disk in our study is compatible with the previous estimate of Thi et al. (2005), that is, $88^{\circ}_{-35} + 2^{\circ}_{+8}$ and that of Tatulli et al. (2008), $82^{\circ}_{-15} + 8$. Moreover, the disk position angle that we derived, that is, $122^{\circ}_{-25}$, is also compatible with that of Tatulli et al. (2008), $126^{\circ}_{-5} + 15$, and is likewise consistent with the angle found by Stark et al. (2009), $126^{\circ}_{-0.15} + 0.5$. Finally, the position angle of the gaseous disk in our study also agrees well with the position angle found for the major axis of the central star, that is, $138 \pm 3.9^{\circ}$, showing that the disk is in the equatorial plane.

## 5. Discussions

The first main result from our CHARA/VEGA observing campaign on 51 Oph is the measurement of the highly flattened stellar photosphere. Assuming that the photosphere distortion is due to the high rotational velocity of the star, the stellar rotational rate can be inferred. Assuming a Roche model and for an edge-on star, a $1.45 \pm 0.12$ distortion corresponds to an angular rotation rate of $\Omega/\Omega_c = 0.99 \pm 0.02$, or a linear one of $V/V_c = 0.96 \pm 0.1$. Another method for deriving the rotational rate is to use the mass and radius of the star and the measured $v \sin i$. The critical velocity is given by $V_c = \sqrt{GM_*/R_c} = \sqrt{2GM_*/(3R_p)}$, where the “e” and “p” subscripts correspond to the equatorial and polar values. Knowing the mass of 51 Oph, $4 M_\odot$ (van den Ancker et al. 1998) and $R_{eq}$, $8.08 R_\odot$, the calculated critical velocity is $308 \text{ km s}^{-1}$. Neglecting the stellar inclination angle, the rotational velocity of the star, that is, $267 \text{ km s}^{-1}$, is 86% of the critical one in linear velocity. As the star might not be seen fully edge-on, this represents a lower limit of the rotational rate. These two different methods give coherent results, and it is clear that 51 Oph is very close to critical rotation.

Assuming critical rotation for a star with a rotational velocity of $267 \text{ km s}^{-1}$ and a polar radius of $R_p = 5.66 \pm 0.23 R_\odot$, we can deduce a mass of $3.3 \pm 0.1 M_\odot$. This is 25% lower than the estimate of van den Ancker et al. (1998), who assumed pre-main sequence evolutionary tracks, that is, $4 M_\odot$, and 15% lower than the estimate of Meeus et al. (2001), who assumed post-main sequence evolutionary tracks, that is, $3.8 M_\odot$. Nevertheless, even if our measurement favors the post-main sequence hypothesis, the difference between the two estimated masses is too small to conclude on the evolutionary status of 51 Oph evolutionary.

The second important result of our study is the direct measurement of the polar radius. Considering the estimate of the effective temperature by Van den Ancker et al. (1998) of $10\,000$ K, it is too high for the star to be classified as a main-sequence star. The standard value for a B9.5V to A0V star is about $2.4-2.6 R_\odot$. Accordingly, 51 Oph should be classified as a giant star. This distance from the main sequence, which was already found by Ancker et al. (1998), was used to compute the age, assuming that 51 Oph was a pre-main sequence star still in the contraction phase. Nevertheless, as we stated above, this value is also compatible with a post-main sequence star.

The last main result from our CHARA/VEGA observing campaign concerns the circumstellar environment geometry and kinematics. We first note that the extension in the $H\alpha$ line, $0.20 \pm 0.07$ AU, is compatible with the one found by Tatulli et al. (2008) in the K-band CO lines using the VLTI/AMBER instrument, $0.15^{+0.07}_{-0.04}$ AU. This emission is also lower than the dust sublimation radius of 51 Oph computed using Eq. (1) from Dullemond & Monnier (2010), $1.17$ AU (i.e., with $T_{\text{eff}} = 10\,000$ K $T_{\text{rim}} = 1500$ K and $R_* = R_p$). With our kinematic model, we showed that the spectro-interferometric data in the

$H\alpha$ emission line are fully compatible with a geometrically thin disk in Keplerian rotation. We note that the differential phase variations cannot be reproduced with an expanding wind. The disk features are similar to those of the few classical Be stars studied with the same instrument (Meilland et al. 2011; Delaa et al. 2011), even if the extension in $H\alpha$ is slightly smaller in the case of 51 Oph. On the other hand, the emission could come from an accretion disk whose geometry and kinematics would not differ strongly from the one of excretion disk of classical Be stars. The main difference between the two models would be the radial velocity, positive for the excretion disk and negative for the accretion disk. However, this low radial velocity of about $0.1 \text{ to } 1 \text{ km s}^{-1}$ is negligible compared to the high rotational velocity in the disk and cannot be constrained in our data.

Finally, from the visible point of view, 51 Oph presents all the features of a classical Be star: critical rotation, double-peaked $H\alpha$ line in emission, and a circumstellar gas disk of a few stellar radii in Keplerian rotation. However, this does not explain the presence of dust as seen in the mid-infrared and millimeter range, and the evolutionary status of 51 Oph still remains unclear. Nevertheless, even if 51 Oph were a pre-main-sequence star, it would make a perfect progenitor for a classical Be star because it is very unlikely that it will lose enough angular momentum to significantly slow down during its on-going contraction phase.

**Acknowledgements.** The CHARA Array is operated with support from the National Science Foundation through grant AST-1211129, the W. M. Keck Foundation, the NASA Exoplanet Science Institute, and from Georgia State University. This research has made use of the SearchCal service of the Jean-Marie Mariotti Center and of CDS Astronomical Databases SIMBAD and VIZIER. N.J. warmly thanks all the VEGA observers that permitted the acquisition of this set of data. N.J. also acknowledges the Ph.D. financial support from the Observatoire de la Côte d'Azur and the PACA Region. We thank T. Lanz for useful discussions. This research has largely benefited from the advices of our colleague Olivier Chesneau, who passed away this Spring. The authors wish to acknowledge his help.

## References

* Abt, H. A., & Morrell, N. I. 1995, ApJS, 99, 135

* Abt, H. A., Levato, H., & Grosso, M. 2002, ApJ, 573, 359

* Aufdenberg, J. P., Mérand, A., Coudé du Foresto, V., et al. 2006, ApJ, 645, 664

* Berthoud, M. G. 2008, Ph.D. Thesis, Cornell University

* Bonneau, D., Clausse, J.-M., Delfosse, X., et al. 2006, A&A, 456, 789

* Delaa, O., Stee, P., Meilland, A., et al. 2011, A&A, 529, A87

* Domiciano de Souza, A., Kervella, P., Jankov, S., et al. 2003, A&A, 407, L47

* Dunkin, S. K., Barlow, M. J., & Ryan, S. G. 1997a, MNRAS, 286, 604

* Dunkin, S. K., Barlow, M. J., & Ryan, S. G. 1997b, MNRAS, 290, 165

* Hubeny, I. 1998, in 1997 Pacific Rim Conference on Stellar Astrophysics, eds.

* K. L. Chan, K. S. Cheng, & H. P. Singh, ASP Conf. Ser., 138, 139

* Malfait, K., Bogaert, E., & Waelkens, C. 1998, A&A, 331, 211

* Meeus, G., Waters, L. B. F. M., Bouwman, J., et al. 2001, A&A, 365, 476

* Meilland, A., Stee, P., Vannier, M., et al. 2007, A&A, 464, 59

* Meilland, A., Delaa, O., Stee, P., et al. 2011, A&A, 532, A80

* Meilland, A., Millour, F., Kanaan, S., et al. 2012, A&A, 538, A110

* Monnier, J. D., Zhao, M., Pedretti, E., et al. 2007, Science, 317, 342

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

* Mourard, D., Bério, P., Perraut, K., et al. 2011, A&A, 531, A110

* Peterson, D. M., Hummel, C. A., Pauls, T. A., et al. 2006, Nature, 440, 896

* Stark, C. C., Kuchner, M. J., Traub, W. A., et al. 2009, ApJ, 703, 1188

* Tatulli, E., Malbet, F., Ménard, F., et al. 2008, A&A, 489, 1151

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

* Thi, W.-F., van Dalen, B., Bik, A., & Waters, L. B. F. M. 2005, A&A, 430, L61

* Thi, W. F., Ménard, F., Meeus, G., et al. 2013, A&A, 557, A111

* van Belle, G., Ciardi, D. R., Thompson, R. R., Akeson, R. L., & Lada, E. A. 2001, BAAS, 33, 881

* van den Ancker, M. E., Meeus, G., Cami, J., Waters, L. B. F. M., & Waelkens, C. 2001, A&A, 369, L17

* van Leeuwen, F. 2007, A&A, 474, 653