A&A 484, L13-L16 (2008)
DOI: 10.1051/0004-6361:200809765

LETTER TO THE EDITOR

The close-in companion of the fast rotating Be star Achernar[*],[*]

P. Kervella1 - A. Domiciano de Souza2 - Ph. Bendjoya2


1 - LESIA, Observatoire de Paris, CNRS UMR 8109, UPMC, Université Paris Diderot, 5 place Jules Janssen, 92195 Meudon, France
2 - Lab. H. Fizeau, CNRS UMR 6525, Univ. de Nice-Sophia Antipolis, Observatoire de la Côte d'Azur, 06108 Nice Cedex 2, France

Received 11 March 2008 / Accepted 16 April 2008

Abstract
Context. The Be stars are massive dwarf or subgiant stars that present temporary emission lines in their spectrum, and particularly in the H${\alpha}$ line. The mechanism triggering these Be episodes is currently unknown, but binarity could play an important role.
Aims. Previous observations with the VLT/VISIR instrument (Kervella & Domiciano de Souza 2007, A&A, 474, L49) revealed a faint companion to Achernar, the brightest Be star in the sky. The present observations are intended to characterize the physical nature of this object.
Methods. We obtained near-IR images and an H-band spectrum of Achernar B using the VLT/NACO adaptive optics systems.
Results. Our images clearly show the displacement of Achernar B over a portion of its orbit around Achernar A. Although there are not enough data to derive the orbital parameters, they indicate a period of about 15 yr. The projected angular separation of the two objects in December 2007 was less than 0.15 $\hbox{$^{\prime\prime}$ }$, or 6.7 AU at the distance of Achernar.
Conclusions. From its flux distribution in the near- and thermal-infared, Achernar B is most likely an A1V-A3V star. Its orbital period appears similar to the observed pseudo-periodicity of the Be phenomenon of Achernar. This indicates that an interaction between A and B at periastron could be the trigger of the Be episodes.

Key words: stars: individual: Achernar - techniques: high angular resolution - stars: emission-line, Be - stars: binaries: close

1 Introduction

As the brightest (mV=0.46) and nearest Be star in the sky, Achernar (${\alpha}$ Eri, HD 10144) has been the focus of a lot of interest over the past decades. Its very fast rotation velocity $v \sin i$ is estimated between 220 to 270 km s-1 and its effective temperature between 15 000 to 20 000 K (Vinicius et al. 2006). Achernar was chosen as the subject of the first VLTI observations, which revealed its extraordinarily distorted interferometric profile (Domiciano de Souza et al. 2003). Different possibilities have been proposed recently to explain the exceptionally high flattening ratio of the photosphere of the star (Jackson et al. 2004; Carciofi et al. 2008). Further interferometric observations have revealed the presence of the stellar wind emitted by the overheated poles of the star (Kervella & Domiciano de Souza 2006), resulting in a slight revision of its flattening ratio. A model of the envelope of Achernar has recently been presented by Meilland (2007). Last year, we discovered a close-in faint companion to Achernar, from diffraction-limited thermal IR imaging with VLT/VISIR (Kervella & Domiciano de Souza 2007). The present Letter reports the follow-up adaptive optics observations in the near-IR domain to characterize this companion, hereafter referred to as Achernar B.

   
2 Observations

2.1 Imaging

We observed Achernar at several epochs in the second half of 2007 using the Nasmyth Adaptive Optics System (NAOS, Rousset et al. 2003) of the Very Large Telescope (VLT), coupled to the CONICA infrared camera (Lenzen et al. 1998), abbreviated as NACO. Table A.1 gives the list of the observations of Achernar and the standard star $\delta $ Phe (HD 9362). We selected the smallest available pixel scale of 13.26 $\pm$ 0.03 mas/pix (Masciadri et al. 2003), giving a field of view of 13.6 $\hbox{$^{\prime\prime}$ }$ $\times$ 13.6 $\hbox{$^{\prime\prime}$ }$. Due to the brightness of Achernar, we employed narrow-band filters at wavelengths 1.094 $\pm$ 0.015, 1.644 $\pm$ 0.018, and 2.166 $\pm$ $0.023~\mu$m (hereafter abbreviated as 1.09, 1.64, and 2.17) together with a neutral density filter (labeled ``ND2_short''), with a transmission of about 1.5%. The raw images were processed using the Yorick[*] and IRAF[*] software packages in a standard way, except that we did not subtract the negligible sky background. Examples of the images of Achernar A and B and $\delta $ Phe are presented in Fig. 1.


  \begin{figure}
\par\includegraphics[width=4cm,clip]{09765fig1a.eps}\includegraphics[width=4cm,clip]{09765fig1b.eps}\end{figure} Figure 1: NACO images of Achernar A and B ( left) and $\delta $ Phe ( right) in the 2.17 $\mu $m filter obtained on 22 December 2007. The arrows indicate the position of Achernar B. The field of view is $1\hbox {$^{\prime \prime }$ }$, and the grey scale is normalized for both images to the central part of the PSF.
Open with DEXTER

2.1.1 Astrometry

To obtain the position of B, we first prepared the images by subtracting a 180$^\circ $ rotated version of each image to itself. This removes the contribution of Achernar A and mainly leaves a pair of positive and negative images of B, corrected from the PSF wing leaks from A (Fig. 2). We then directly measured the position of Achernar B relative to A using a Gaussian fit procedure on these images. Although Achernar B is clearly visible at all three wavelengths, the Gaussian fit converged essentially in the K band where the high Strehl ratio (compared to the 1.09 $\mu $m images in particular) results in a better separation of the two objets. The formal fitting error is below $\pm$0.1 pixel, but we estimate the true uncertainty to $\pm$0.5 pixel (0.066 $\hbox{$^{\prime\prime}$ }$), due to the presence of residual speckles from Achernar A close to B. The pixel scale and detector orientation introduce negligible systematic uncertainties (Masciadri et al. 2003; Chauvin et al. 2005).


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{09765fig2.eps}\end{figure} Figure 2: Result of the subtraction of a 180$^\circ $ rotated version of the image of the Achernar system presented in Fig. 1 ( left) from itself. The ``B+'' and ``B-'' arrows point at the positions of the positive and negative images of Achernar B, and the ``A res.'' arrow indicates the A subtraction residuals. The aperture used for photometry is shown in the lower left corner, and the field of view is $1\hbox {$^{\prime \prime }$ }$ on each side.
Open with DEXTER

Table 1: Position of Achernar B relative to Achernar A.

Between epochs 2006.760 (Kervella & Domiciano 2007) and 2007.975, we measured an apparent displacement of B relative to A of $\rho_{\rm B} = 161$ $\pm$ 10 mas along an azimuth of $\alpha = +113^\circ$ (Fig. 3 and Table 1). According to the Hipparcos catalog (ESA 1997), the amplitude of Achernar A's proper motion on the sky over 1.2 yr is $\rho_{\rm pm} = 116$ mas along an azimuth of $\alpha_{\rm pm} = +115^\circ$. If B was a background source, we would expect an apparent displacement opposite in azimuth to A's proper motion. As B is clearly comoving with A, we can rule out the possibility that it is a background source. As a remark, the parallactic oscillation of the position of Achernar A is small ( $\pi = 22.68$ $\pm$ 0.57 mas) compared to the observed displacements, but the presence of B could have affected the Hipparcos proper motion and parallax measurement of Achernar. The other epochs listed in Table A.1 give astrometric positions compatible with the trajectory shown in Fig. 3, although with lower accuracy due to poor seeing conditions.


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{09765fig3.eps}\end{figure} Figure 3: Position of Achernar B relative to A for four epochs. The open square indicates the VISIR observation (Kervella & Domiciano 2007), and the dots represent the new NACO epochs. The dashed curve is a simple quadratic fit through the data points intended to guide the eye. The segment over Achernar A indicates its projected rotation axis and polar wind as measured by Kervella & Domiciano (2006). The apparent angular sizes of Achernar A and its polar wind are approximately represented to scale.
Open with DEXTER

2.1.2 Photometry

For the photometric calibration of our images, we observed a standard star immediately before or after Achernar, $\delta $ Phe (HD 9362, G9III) in the same narrow-band filters. It was chosen in the Cohen et al. (1999) catalogue of spectrophotometric standards for infrared wavelengths. Due to the small separation of the Achernar A-B pair, PSF fitting is not possible for measuring the photometry of B. We therefore proceeded in three steps:

1.
We measured classical aperture photometry of Achernar relative to $\delta $ Phe, using a large aperture of 1.2 $\hbox{$^{\prime\prime}$ }$ in radius.

2.
As for the astrometry measurement, we subtracted a 180$^\circ $ rotated version of each image to itself to remove the contribution of Achernar A and leave a positive and a negative image of B, along with limited residuals from A (Fig. 2).

3.
We obtained aperture photometry of B from the subtracted images using an aperture radius of 0.040 $\hbox{$^{\prime\prime}$ }$ (shown in Fig. 2).
This procedure has ensured that we obtain properly referenced photometry for Achernar A (the influence of B is negligible at step 1), that we can transfer to B at step 3. Step 2 is necessary because B is located inside the wings of the PSF of A. This subtraction is efficient thanks to the good circular symmetry of the PSF produced by NACO. At step 3, we measure the photometry of A and B over the same aperture and within the same images. This allows us to avoid the problem of the variable Strehl ratio of adaptive optics images. The conversion of the measured narrow-band magnitudes to standard $\it JHK$ band magnitudes requires taking the position of the quasi-monochromatic wavelengths within the bands and the shape of the observed spectra into account. Starting from the narrow-band fluxes measured on Achernar and $\delta $ Phe, we used the Pickles (1998)[*] reference spectra corresponding to their spectral types (B3V and G9III) to recover the corresponding broadband flux ratio. The filter profiles were taken from Bessell & Brett (1988). From the October-December 2007 observations, we derived magnitudes of: $m_J({\rm A}) = 0.58_{\pm 0.14}$, $m_H({\rm A})= 0.80_{\pm 0.11}$, $m_K({\rm A}) = 0.81_{\pm 0.11}$. These values generally agree with the $\it JHK$ magnitudes of Achernar from the 2MASS (Skrutskie et al. 2006; $0.82_{\pm 0.25}$, $0.87_{\pm 0.32}$, $0.88_{\pm 0.33}$) and Ducati (2002; 0.79, 0.86, 0.88) catalogs, although systematically brighter by 0.1-0.2 mag, particularly in the J band. This may be due our narrow-band filters corresponding to the emission lines sometimes present in the spectrum of Achernar. As we may slightly overestimate the brightness of Achernar, we quadratically added an uncertainty of 0.1 mag.

From the magnitudes of Achernar A, the small-aperture photometry obtained at step 3 gives the following broadband magnitudes for Achernar B: $m_J({\rm B}) = 4.23_{\pm 0.36}$, $m_H({\rm B})= 4.29_{\pm 0.22}$, $m_K({\rm B}) = 4.61_{\pm 0.27}$. The magnitude differences measured between A and B in the three narrow-band filters are: $\Delta m_{1.09} = 3.64_{\pm 0.47}$, $\Delta m_{1.64} = 3.50_{\pm 0.19}$, $\Delta m_{2.17} = 3.80_{\pm 0.25}$, corresponding to an average contrast of $\approx$30 between the two stars in the near-infrared. These magnitudes of B are averages over the October-December 2007 observations (epoch 2007.8) and the error bars contain the statistical uncertainty and the dispersion of all the measurements. The magnitudes at the other epochs (list of the observations in Table A.1) show a marginally significant brightening of B by $-0.7_{\pm 0.4}$ mag in the K band between 2007.488 and 2007.975 as it approaches A.

2.2 Spectroscopy

For the spectroscopic observations, we used the S27_3_SH mode of NACO, featuring a slit width of 86 mas on the sky, a spectral resolution of 1500 over the H band (dispersion of 0.34 nm/pixel), and an angular pixel scale of 27 mas/pixel. As shown in Fig. 4, the slit was centered on Achernar B. Although most of the light fed into the slit comes from the companion, it is not excluded that Achernar A may contribute with part of the H band spectrum. Details on the 3 recorded ABBA sequences are given in Table A.1. Observations were performed using the classical IR NACO sequence where the spectrum is recorded with the target positioned at two different slit positions, i.e. the ABBA observing sequence. For the data reduction we applied two different procedures: (1) the classical spectral reduction procedure using the IRAF packages for dark and flat field corrections, and (2) the A-B, B-A procedure. Because the sky background is negligible and the target is rather bright, both procedures give identical results. In the last step, IRAF was used to obtain 1D spectra from weighted averages of 2D spectra and to perform the wavelength calibration. Figure 5 shows the continuum-normalized average spectrum derived from our 12 observations. The atmospheric absorption lines were corrected using the transmission by Lord (1992), but some residuals are visible in particular around $\lambda=1.66~\mu$m (CH4 and water vapor lines). The absorption lines from the Brackett transitions 11 to 23 of hydrogen are clearly visible in the spectrum and marked by vertical lines.


  \begin{figure}
\par\includegraphics[width=7.9cm,clip]{09765fig4.eps}\end{figure} Figure 4: Position of the NACO slit relative to Achernar A and B.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=16.6cm,clip]{09765fig5.eps}\end{figure} Figure 5: Average spectrum of Achernar B (solid curve) compared to the H band spectrum of HR 2763 ($\lambda $ Gem, A3V, dashed curve) from Ranade et al. (2004). The positions of the Brackett transitions 11-4 to 23-4 of hydrogen are marked with vertical lines.
Open with DEXTER

   
3 Discussion

The Hipparcos parallax of Achernar ( $\pi = 22.68$ $\pm$ 0.57 mas; ESA 1997) corresponds to a distance modulus of $\mu=3.22$ $\pm$ 0.06 mag. The absolute magnitudes of B in the $\it JHK$ and thermal infrared N band (average of the PAH1 and PAH2 magnitudes, see Appendix A) are therefore: $M_J = 1.00_{\pm 0.37}$, $M_H = 1.07_{\pm 0.22}$, $M_K = 1.39_{\pm 0.28}$, and $M_N = 1.88_{\pm 0.12}$. The $\it JHK$ absolute magnitudes suggest a spectral type around A1V, slightly fainter than Vega and slightly brighter than Sirius, which have respective magnitudes of $M_{JHK} \approx 0.55$ and 1.56. The N band absolute magnitude of Achernar B is also very similar to that of Sirius ( MN = 1.62). It thus does not appear to present an infrared excess that would betray a Vega-like dusty envelope. Because Vega itself presents a significant excess in this band ( MN = 0.55) due to the presence of dust, it is not comparable.


  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{09765fig6.eps}\end{figure} Figure 6: Cross-correlation between our spectrum of B and reference spectra from Ranade et al. (2004), Meyer et al. (1998), and Pickles (1998).
Open with DEXTER

Our NACO H band spectrum of B presents a lot of similarity with early A-type star spectra, as shown in Fig. 5. This is confirmed by the computation of linear Pearson cross-correlation coefficients between our normalized spectrum and the reference H band spectra from Ranade et al. (2004)[*] and Pickles (1998)[*]. As shown in Fig. 6, the spectral types A0V-A4V give the best correlation with our observation. As shown by Meyer et al. (1998)[*], the equivalent width (EW) of the 11-4 Brackett line at $1.681~\mu$m offers a possibility to test the effective temperature, although it is not monotonic. An integration normalized to the pseudo-continuum between 1.670 and $1.690~\mu$m gives an equivalent width of 0.675 nm (=2.4 cm-1). From Fig. 7 in Meyer et al. (1998), this gives $\log T_{\rm eff} \approx 3.9{-}4.1$, also compatible with an early A-type star. Such a star has an approximate mass of $\approx$ $2~M_{\odot}$ (Kervella et al. 2003).

Although we cannot derive the full parameters of Achernar B's orbit from our limited astrometry, the combination of its estimated mass ($\approx$$2~M_{\odot}$), the mass of A ($\approx$ $6.7~M_\odot$; Vinicius et al. 2006, see also Harmanec 1988) and the maximum apparent A-B separation ( $r_{\rm AB} \approx 12.3$ AU) allow us to roughly estimate its period. From Kepler's third law and by assuming an elliptic orbit, we have T2 = a3/M where a is the semi-major axis in AU, T the period in years, and M the total mass in $M_\odot$. Assuming that a is equal to the observed maximum separation of 12.3 AU, we obtain a minimum period of T = 14 to 15 yr.

4 Conclusion

From our photometry and spectroscopy, Achernar B is most probably an A1V-A3V star. Our data are currently insufficient for deriving its full orbit, but its minimum period is $\approx$15 yr. Its fast orbital motion should allow the derivation of reliable parameters within a few years. The periodic approach of the companion could be the cause of the observed $\approx$15 yr pseudo-periodicity of the Be episodes of Achernar (Vinicius et al. 2006). The passage of B at periastron within a few AUs of A could extract material from the equator of A, where the effective gravity is very low. The next such passage should happen around 2010 (Meilland 2007). Achernar appears similar to the B0.2IVe star $\delta $ Sco, which has a 1.5 mag fainter companion on a highly excentric 10.6 yr orbit (Bedding 1993; Miroshnichenko et al. 2001). This suggests that the presence of companions around Be stars should be examined carefully, as it may play a key role in triggering the Be phenomenon.

Acknowledgements
This research made use of the SIMBAD and VIZIER databases at the CDS, Strasbourg (France), and NASA's Astrophysics Data System Bibliographic Services. We also received the support of PHASE, the high angular resolution partnership between ONERA, Observatoire de Paris, the CNRS, and University Denis Diderot Paris 7.

References

 

  
Online Material

   
Appendix A: Thermal infrared imaging at 11.25 $\mu $m

Kervella & Domiciano de Souza (2007) did not detect Achernar B in the 11.25 $\mu $m images (PAH2 band) of the system obtained with the VLT/VISIR instrument in BURST mode. A re-analysis of these observations has shown that the reason for this non-detection was the introduction of noise during the normalization and subtraction of the PSF reference star ($\delta $ Phe). For our VISIR program on the B[e] star MWC300 (Domiciano de Souza et al. 2008), we obtained BURST mode images of another PSF calibrator ($\eta $ Ser) 3.4 h before the observations of Achernar (under similar seeing conditions). Since $\eta $ Ser is brighter than $\delta $ Phe, we could obtain a cleaner PSF-subtracted image of Achernar B, as presented in Fig. A.1. From this image, aperture photometry over a 0.22 $\hbox{$^{\prime\prime}$ }$ diameter gives a flux ratio of 1.74% between Achernar B and A (the ratio of peak intensities is 2.14%). From the 16.8 Jy absolute flux derived by Kervella & Domiciano de Souza (2007) for A, the contribution from B is therefore 0.3 Jy at 11.25 $\mu $m. This value is comparable to the 0.4 Jy flux derived by Kervella & Domiciano de Souza (2007) at 8.59 $\mu $m.


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{09765fig7.eps}\end{figure} Figure A.1: Residual of the subtraction of the VISIR PAH2 image of $\eta $ Ser from the image of Achernar in the same band. The contribution of Achernar B is clearly visible. The arrow indicates the rotation axis of Achernar A.
Open with DEXTER

Table A.1: Log of the NACO observations of Achernar and $\delta $ Phe.



Copyright ESO 2008