* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download Astronomy Astrophysics Circumstellar emission in Be/X-ray binaries of the Magellanic
Aries (constellation) wikipedia , lookup
Corona Borealis wikipedia , lookup
X-ray astronomy satellite wikipedia , lookup
Circumstellar habitable zone wikipedia , lookup
Canis Minor wikipedia , lookup
Auriga (constellation) wikipedia , lookup
Constellation wikipedia , lookup
Gamma-ray burst wikipedia , lookup
Cassiopeia (constellation) wikipedia , lookup
Corona Australis wikipedia , lookup
Perseus (constellation) wikipedia , lookup
Aquarius (constellation) wikipedia , lookup
Cygnus (constellation) wikipedia , lookup
Observational astronomy wikipedia , lookup
International Ultraviolet Explorer wikipedia , lookup
Timeline of astronomy wikipedia , lookup
X-ray astronomy wikipedia , lookup
Malmquist bias wikipedia , lookup
Astrophysical X-ray source wikipedia , lookup
Star catalogue wikipedia , lookup
Stellar evolution wikipedia , lookup
H II region wikipedia , lookup
Corvus (constellation) wikipedia , lookup
High-velocity cloud wikipedia , lookup
Cosmic distance ladder wikipedia , lookup
Stellar classification wikipedia , lookup
Astronomy & Astrophysics A&A 539, A114 (2012) DOI: 10.1051/0004-6361/201117738 c ESO 2012 Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way M. S. Riquelme1 , J. M. Torrejón1,2 , and I. Negueruela1 1 2 Departamento de Física, Ingeniería de Sistemas y Teoría de la Señal, Escuela Politécnica Superior, Universidad de Alicante, Spain e-mail: [email protected] Instituto de Física Aplicada a las Ciencias y las Tecnologías, Universidad de Alicante, 03690 Alicante, Spain Received 21 July 2011 / Accepted 20 December 2011 ABSTRACT Aims. We study the optical and near-infrared colour excesses produced by circumstellar emission in a sample of Be/X-ray binaries. Our main goals are exploring whether previously published relations, valid for isolated Be stars, are applicable to Be/X-ray binaries and computing the distance to these systems after correcting for the effects of the circumstellar contamination. Methods. Simultaneous U BVRI photometry and spectra in the 3500−7000 Å spectral range were obtained for 11 optical counterparts to Be/X-ray binaries in the LMC, 5 in the SMC and 12 in the Milky Way. As a measure of the amount of circumstellar emission we used the Hα equivalent width corrected for photospheric absorption. Results. We find a linear relationship between the strength of the Hα emission line and the component of E(B − V) originating from the circumstellar disk. This relationship is valid for stars with emission lines weaker than EW ≈ −15 Å. Beyond this point, the circumstellar contribution to E(B − V) saturates at a value ≈0.17 mag. A similar relationship is found for the (V − I) near infrared colour excess, albeit with a steeper slope and saturation level. The circumstellar excess in (B − V) is found to be about five times higher for Be/X-ray binaries than for isolated Be stars with the same equivalent width EW(Hα ), implying significant differences in the physical properties of their circumstellar envelopes. The distance to Be/X-ray binaries (with non-shell Be star companions) can only be correctly estimated by taking into account the excess emission in the V band produced by free-free and free-bound transitions in the circumstellar envelope. We provide a simple method to determine the distances that includes this effect. Key words. circumstellar matter – X-rays: binaries – stars: early-type 1. Introduction Be/X-ray binaries comprise the major subclass within the high mass X-ray binaries, and the majority of X-ray pulsars belong to this class (see Reig 2011, for a comprehensive recent review). They are composed of a neutron star (NS) that accretes circumstellar material from the disk around a massive Be star. Be stars in X-ray binaries have spectral types restricted to a fairly narrow interval, O9 to B3 and luminosity classes V to III (Negueruela 1998). These stars are surrounded by a circumstellar envelope that profoundly alters the emerging spectrum. Spectra of Be stars show some essential characteristics. 1) The Balmer lines are affected by the emission that is produced by bound-bound transitions in the envelope, which fills in the photospheric lines or drives them into emission. 2) Free-free transitions in the envelope (thermal Bremsstrahlung) produce an excess of photons at long wavelengths, compared to a B star of the same spectral type. This effect is referred to as infrared excess and increases with wavelength. 3) Thermal Bremsstrahlung also makes the star appear brighter, compared to a B star of the same luminosity class. Other mechanisms can help to modify even further the brightness of the system. Thomson scattering in the envelope, for example, could add photons into the observer’s line of sight, which would additionally increase the brightness of the star. Partially based on observations collected at the European Southern Observatory, La Silla, Chile (66.D-0292; 074.D-0529) and the Nordic Optical Telescope. Owing to these effects, the direct application of the photometric and spectroscopic techniques for calculating the physical parameters of a star leads to incorrect results. In particular, the distances to some Be/X-ray binaries, and therefore their X-ray luminosities, could be underestimated, as we will see. All these effects are directly connected with the envelope and, therefore, a relationship between them and the overall amount of emission is expected. In this paper, we present a method to correctly treat these effects by means of simple correlations between the circumstellar excess in (B − V) and (V − I) colours and the amount of emission, as measured by the Hα equivalent width. To this goal, we have observed a sample of Be/X-ray binaries in the Magellanic Clouds (LMC and SMC) and in the Milky Way galaxy (MWG). Such relations had been established so far for isolated Be stars, that is to say, Be stars not known to be in binaries (Dachs et al. 1988, for bright field stars, Fabregat & Torrejón 1998, for Be stars in open clusters). Nevertheless, we now know that the circumstellar envelope of Be stars in binaries is substantially modified by the presence of the accreting compact object. Circumstellar envelopes in binary systems tend to be geometrically more compact and optically thicker (Reig et al. 1997; Okazaki & Negueruela 2001). It has been shown in the past years by direct interferometric results that the disk is indeed truncated in all cases studied (Tycner et al. 2006; Gies et al. 2007; Gründstrom & Gies 2006). There is evidence that the circumstellar reddening given by the calibrations for isolated Article published by EDP Sciences A114, page 1 of 13 A&A 539, A114 (2012) Be stars is too small when applied to Be/X-ray binaries for a given amount of emission (Zamanov et al. 2001). Therefore, the question arises whether these relationships can be safely applied as well to Be/X-ray binaries. In this paper we aim to answer this question. 2. Observations 2.1. Sample selection and instrumentation In order to be able to decouple the interstellar reddening and extinction from the circumstellar reddening and brightening1 we need a set of programme stars with well-known distances and interstellar reddenings as well as good determinations of the spectral types. In principle we could choose Be/X-ray binaries in well-studied open clusters and follow an approach similar to the one used by Fabregat & Torrejón (1998). Unfortunately, such stars do not seem to exist or are very scarce. Therefore we used stars in the Large Magellanic Cloud (LMC) whose distance and reddening are satisfactorily studied by independent methods. Furthermore, the LMC presents a population distribution that is similar to the Milky Way Galaxy (MWG; Antoniou et al. 2009) as well as a similar metallicity. The sample is almost complete because we have observed 11 out of 15 BeX known in the LMC. As a control sample we included five BeX in the Small Magellanic Cloud (SMC), chosen randomly amongst those with already determined and published spectral types. Finally we have chosen a sample of BeX in the MWG covering a wide range of emission strengths and spectral types. All 28 observed stars have spectral types already determined and published (the references are given in Table 4). This is mandatory, because the parameters of the underlying B star will be used as input for the calibration. Because Be stars are variable even on timescales of hours, we obtained simultaneous U BVRI Bessel photometry and Balmer spectroscopy in three different runs using the following telescopes and instrumentation: – The Danish Faint Object Spectrograph and Camera (DFOSC), attached to the Danish 1.54 m telescope, located at the ESO site in La Silla (Chile) during a run on 2001 February, 8–10. – The Andalucian Faint Object Spectrograph and Camera (ALFOSC), installed on the 2.5 m Nordic Optical Telescope (NOT) in La Palma (Spain) during a run on 2001 July, 2–5. – The ESO Multi-Mode Instrument (EMMI), at the 3.6 m New Technology Telescope (NTT), also located in La Silla, on 2004 December, 14 & 15. These instruments allow the sequential acquisition of images and spectra, with virtually no time delay, which is mandatory for reducing the uncertainty produced by short-term variability. 2.2. Photometry The photometric images were taken sequentially with the Bessell U BVRI (broadband filters). In some runs only BVR photometry is available. For data reduction we used the Starlink suite2 . Bias subtraction and flatfielding were carried out on all frames using the ccdpack package (Draper et al. 2000). Photometry was performed using the PSF fitting optimum extraction method implemented in gaia (Draper & Gray 2000). 1 2 For classical Be stars. See discussion below for Be shell stars. Now available at http://starlink.jach.hawaii.edu/ A114, page 2 of 13 Calibration of our instrumental magnitudes was carefully determined from eight different standards fields, taken from the list of Landolt (1995): Lan 95, Rubin 152 and PG0942 for the first run, Lan 107−602, Lan 110−230 and Lan 113−163 for the second and Lan 96 and Lan 97 for the third run, respectively. All these fields were followed throughout the night at a wide range of airmasses. The standardization process was carried out using the following set of equations: U − u = ξu + Au (U − B) + Ku X + Ku X (U − B) B − b = ξb + Ab (B − V) + Kb X + Kb X (B − V) V − v = ξv + Av (B − V) + Kv X + Kv X (B − V) R − r = ξr + Ar (V − R) + Kr X + Kr X (V − R) I − i = ξi + Ai (R − I) + Ki X + Ki X (R − I), (1) where U and u are the standard and instrumental magnitudes, ξu is the zero point, Ku and Ku are extinction coefficients, X is the airmass and Au is a constant. Similar notation is used for the remaining filters. Final calibration equations were applied to all programme stars, resulting in the magnitudes quoted in Table 1. Internal mean errors were judged from the scatter between various measurements of the PSF fitting and photon-counting statistics. These were always less than 0.01 mag. After the standardization process, we finally assigned an error to our calibration given by the rms of the residuals between listed and calculated magnitudes for the standard stars. These errors are quoted in Table 2. 2.3. Spectroscopy Following each photometric sequence, medium resolution spectroscopy was performed in the spectral classification and Balmer Hα regions. The details of the several instrumental setups used are given in Table 3. Slits of 1 were used in all runs. Spectroscopic data were reduced using Starlink packages ccdpack (Draper et al. 2000) and figaro (Shortridge et al. 1997) using standard techniques. After correcting the images for bias and flatfield, spectra were extracted from the object and the sky. This last spectrum was subtracted from that of the object and the resulting spectrum was wavelength-calibrated. For that purpose He-Ar calibration lamp spectra were taken following each programme star observation. Integration times were chosen to obtain signal to noise ratios (S/N) higher than 150. Finally, several spectrophotometric standard stars were observed each night to allow for absolute flux calibration. In Fig. 1, left panel, we show the spectrum of the programme star KS 1947+300 as an example. Note the strong Hα emission and reddening. All spectra were subsequently analysed using the Starlink dipso package (Howarth et al. 1998). In order to measure the observed equivalent width of the Hα emission line, we normalized the continuum to unity, as shown in Fig. 1, right panel, which allows for a more precise assessment of the wings extension. All lines showed a full width at half maximum (FWHM) much higher than the spectral resolution of the different instrumental setups. Therefore, the impact of spectral resolution on the measured EW is negligible. Errors in the equivalent width measurements were estimated averaging over several adopted continua as well as through several assessments of the wings extension. The process of sky subtraction can introduce non negligible errors around the Hα line because there is often strong and spatially variable diffuse Hα emission in the regions where the Be/X-ray binaries lie in the Magellanic Clouds. In order to asses and minimize these effects we used a narrow 1 slit, which minimizes the sky contribution. Next we extracted sky spectra from several positions along M. S. Riquelme et al.: Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way Table 1. U BVRI photometry for BeX stars in LMC, SMC and the Milky Way are listed in Cols. 3 through 7. Galaxy LMC SMC MWG Star CAL E CAL 9 RX J0516.0−6916 RX J0516.0−6916 RX J0520.5−6932 RX J0529.8−6556 EXO 0531.1−6609 EXO 0531.1−6609 RX J0531.5−6518 RX J0531.5−6518 RX J0535.0−6700 1A 0538−66 1A 0538−66 SAX J0544.1−7100 H 0544−66 RX J0051.9−7311 RX J0052.1−7319 RX J0059.2−7138 RX J0101.0−7206 RX J0101.3−7211 4U 0115+634 3A 0726−260 GRO J1008−57 LS 1698 1A 1118−61 GX 304−1 SS 397 XTE 1946+274 KS 1947+300 Cep X-4 4U 2206+54 SAX J2239.3+6116 U 12.89 13.56 ··· ··· 13.31 13.54 ··· ··· ··· ··· ··· 13.64 ··· ··· ··· ··· ··· ··· ··· ··· 16.96 11.32 17.21 11.34 12.63 16.46 13.97 19.05 14.58 15.20 9.27 16.19 B 14.03 14.72 14.99 15.06 14.46 14.43 14.74 14.76 15.93 15.89 14.86 14.68 15.12 15.35 15.29 14.40 14.86 14.04 15.68 15.47 16.96 12.12 16.74 11.87 12.95 15.63 13.78 18.63 15.01 15.47 10.01 16.06 V 14.17 14.81 15.12 15.15 14.46 14.59 14.85 14.82 15.99 15.99 14.94 14.91 15.25 15.23 15.49 14.42 14.86 14.07 15.72 15.53 15.54 11.61 15.08 11.40 12.09 13.94 12.30 16.89 14.08 14.45 9.67 14.55 R 14.10 14.77 15.15 15.19 14.40 14.50 14.81 14.80 16.05 16.01 14.91 14.93 15.29 15.18 15.58 14.36 14.80 14.06 15.60 15.54 ··· 11.12 13.77 10.95 11.26 12.69 11.49 15.65 13.48 13.78 9.59 13.60 I 14.10 ··· ··· ··· 14.23 14.56 ··· ··· ··· ··· ··· 15.16 ··· ··· ··· ··· ··· ··· ··· ··· ··· 10.88 12.68 10.74 10.69 11.67 10.49 14.78 12.78 13.01 9.39 12.74 Run Danish Danish NTT14 NTT15 Danish Danish NTT14 NTT15 NTT14 NTT15 NTT14 Danish NTT14 NTT14 NTT14 NTT15 NTT15 NTT15 NTT15 NTT15 NOT Danish Danish Danish Danish Danish NOT NOT NOT NOT NOT NOT Notes. U-band and I-band photometry was not performed during the NTT campaign. Table 2. U BVRI photometry errors for the different runs. U 0.03 ··· ··· 0.03 B 0.02 0.03 0.02 0.03 V 0.02 0.03 0.02 0.01 R 0.02 0.03 0.02 0.02 I 0.03 ··· ··· 0.06 Run Danish NTT14 NTT15 NOT the slit, including at least two positions close to the target spectrum, and checked for possible variations. The background spectra appeared to be quite stable. Finally, they were averaged and the resulting spectrum subtracted from the that of the star. Final errors were calculated adding the spectrum and model contributions quadratically. The estimated equivalent width errors in the model absorption line is ±0.1 Å. The errors in the measured emission line were always lower than 1 Å for the stars in the Clouds and lower than 2 Å for stars in the MWG. Final errors are given in Tables 5 and 8. 3. Circumstellar emission and reddening 3.1. Colour excess As pointed out in the introduction, material surrounding the massive star produces three main effects in the spectrum. In this section we are concerned with the first two. The total excess of photons in the optical is the sum of two components: the circumstellar excess, E cs (B − V), entirely caused by the envelope, and Table 3. Grism characteristics for each instrumental setup. Instrument DFOSC ALFOSC EMMI Grism # 3 5 3 5 5 λc (Å) 3900 6500 4300 7000 5300 Dispersion (Å/pix) 2.55 3.3 2.3 3.1 0.83 Δλ (Å) 5.5 7.5 6.2 8.4 4 Range (Å) 3300–6400 5200–10 200 3200–6700 5000–10 250 3800–7020 Notes. Central wavelengths λc correspond to the maximum efficiency. the interstellar excess, E is (B − V), produced by the interstellar material lying between the system and the observer. Interstellar material consists mainly of gas and dust. The gas absorbs radiation at very specific wavelengths, leaving absorption lines in the spectrum of the star. The dust, however, scatters the star light over a continuum of wavelenghts. The scatter cross section varies as σ ∼ λ−4 . Because blue photons are more easily scattered out of the star light beam than the red photons, the spectrum of the star becomes reddened. On the other hand, the thermal free electrons in the envelope scatter off ions via free-free transitions, thereby emitting a continuum spectrum. This emission spectrum increases its yield towards long wavelengths and is very strong in the infrared. The net effect over the observed spectrum will then be: E total (B − V) = E cs (B − V) + E is (B − V). (2) A114, page 3 of 13 A&A 539, A114 (2012) 1.6 Flux / 10-14 ( erg / cm -2 s -1 A-1 ) 1.5 1.4 1.3 1.2 1.0 0.8 0.7 0.6 0.4 0.2 4500 5000 5500 6000 ο 6500 7000 7500 λ (Α) Fig. 1. Left panel: flux-calibrated spectra of KS 1947+300. Note the strong Hα emission line at 6500 Å and the reddening. Right panel: zoom of the emission line normalized to continuum. The Kurucz atmosphere model with the corresponding Hα is also shown with an offset for a better comparison. The wavelength dependence of these two excesses, or reddenings, is completely different: the interstellar diminishes towards the red, while the circumstellar increases towards the red. But their effects upon the spectrum are indistinguishable. In order to separate these two contributions, we need a set of Be/X-ray binaries with known distances and reddenings. As explained before, no Be/X-ray binaries in open clusters have been detected so far. Therefore we observed stars in the LMC and the SMC because the interstellar reddening is already sufficiently known a priori by independent methods. Equation (2) can be explicitly split into two contributions: ∗ E (B − V) = (B − V) − (B − V)0 cs (3) and E is (B − V) = (B − V)obs − (B − V)∗ , (4) where (B − V)∗ is the colour of the star+envelope system (before interstellar reddening). Therefore the total excess, as defined previously in Eq. (2), can be obtained as E total (B − V) = (B − V)obs − (B − V)0 , (5) where (B − V)obs is the observed colour computed from Table 1 and (B − V)0 is the intrinsic colour for the corresponding spectral type (Table 4) as given by Schmidt & Kaler (1982). For some stars the spectral type is still uncertain. The spectral lines are filled in by emission and this renders the classification quite uncertain in extreme cases. When two spectral types are possible, calculations were made for both, and we chose the value that matched the distance and interstellar reddening of the Clouds correctly. For example, in the case of 1A 0538−66, the luminosity class quoted in the literature is III. However, a class V fits into the expected value for the distance to the LMC much better. Therefore we chose this second value. Intrinsic (B − V) colours, effective temperatures and log g were taken from Schmidt & Kaler (1982). Intrinsic (V − I) colours were taken from Ducati et al. (2001). For consistency, the spectral type deduced from our spectra was compared with the published values whenever possible. A good agreement was found. Owing to the limited spectral resolution of our data we will use the published values. All these parameters are listed in Table 4. The errors for the (B − V)obs were calculated adding the contributions of both filters quadratically. This error was considered A114, page 4 of 13 as the error for the total excess, which is equivalent to assuming that the error in the theoretical values of (B−V) for every spectral type is negligible. Of course this is not strictly correct because we always have an observational error of about ± one subtype in the spectral type determination. Fortunately, as can be seen in Table 4, the spectral types are concentrated in a fairly narrow interval, from B0 to B1. This covers a range in theoretical colours of about Δ(B − V)0 = 0.03, which, when added quadratically to the photometric colour errors, yields a final uncertainty of ±0.06. 3.2. E total (B − V ) vs. W env (Hα ) Next we compared our spectrum with the stellar atmosphere models of Kurucz (1979), for the effective temperatures and gravities given in Table 4. These models represent the spectrum of the underlying star without reddening and without emission. Then we measured the equivalent widths of the Hα line, both in the observed spectrum, W spc (Hα ), normally in emission, and in the model, W mdl (Hα ), always in absorption (see right panel Fig. 1). We normalized the spectra to a continuum of unity because this allows a much better assessment of the true extension of the line wings. Finally we subtracted them to obtain the true emission from the envelope W env (Hα ) corrected for photospheric absorption: W env (Hα ) = W spc (Hα ) − W mdl (Hα ). (6) This value is quoted in Table 5 for LMC and SMC stars and in Table 8 for MWG stars. As explained before, the spectral type range of the targets ranges from B0 to B1. The corresponding equivalent widths cover a range of Δ(EWmodel ) 0.1 Å for the same spectral types (Kurucz 1979). Added quadratically to the observational spectroscopic errors, this results in a final uncertainty of less than 1.5 Å except for the MWG source XTE J1946+274, which is 2.0 Å. In Fig. 2 we display E total (B − V) vs. W env (Hα ) for all our programme stars. Milky Way binaries are represented by (red) squares, LMC binaries by (blue) asterisks and SMC binaries by (green) triangles. Several conclusions can be drawn from inspection of the plot: 1. total excesses for Milky Way Galaxy stars are of the order of 10 times higher than for stars in the Magellanic Clouds; M. S. Riquelme et al.: Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way Table 4. Spectral type in the Morgan-Keenan system and its corresponding parameters T eff , (B − V)0 , MV and log g, taken from Schmidt & Kaler (1982). Star CAL E CAL 9 RX J0516.0−6916 RX J0520.5−6932 RX J0529.8−6556 EXO 0531.1−6609 RX J0531.5−6518 RX J0535.0−6700 1A 0538−66 SAX J0544.1−7100 H 0544−66 RX J0051.9−7311 RX J0052.1−7319 RX J0059.2−7138 RX J0101.0−7206 RX J0101.3−7211 4U 0115+634 3A 0726−260 GRO J1008−57 LS 1698 1A 1118−61 GX 304−1 SS 397 XTE J1946+274 KS 1947+300 Cep X-4 4U 2206+54 SAX J2239.3+6116 MK B0 Va B0 Va B0.5 IVb 09 Va B0.5 Va B0.7 Va B2 Va B0 Va B0.5 Va B0 Vc B0 Va B9.5 Vd B0 Vd B1 IIId B1 V d B0.5 Vd B0.2 V e O8.5 V f B0 III/Ve B0 III/Vg O9.5 Vh B0.7 Vi B1 Vm B0–1 V j B0 Vk B1–2 Vl O9.5 Vm B0 Vi T (K) 30 000 30 000 27 000 33 000 28 000 27 000 22 000 30 000 28 000 30 000 30 000 31 000 30 000 24 000 25 000 28 000 29 000 34 000 29 000/30 000 29 000/30 000 31 000 27 000 25 000 30 000/25 000 30 000 25 000/22 000 31 000 30 000 (B − V)0 −0.30 −0.30 −0.28 −0.31 −0.28 −0.28 −0.24 −0.30 −0.28 −0.30 −0.30 −0.31 −0.30 −0.27 −0.27 −0.28 −0.29 −0.32 −0.30 −0.30 −0.31 −0.28 −0.27 −0.30/−0.27 −0.30 −0.27/−0.24 −0.31 −0.30 MV −4.0 −4.0 −4.2 −4.5 −3.5 −3.4 −2.4 −4.0 −3.5 −4.0 −4.0 −4.3 −4.0 −4.4 −3.2 −3.5 −3.8 −4.7 −5.1/−4.0 −5.1/−4.0 −4.3 −3.4 −3.2 −4.0/ − 3.2 −4.0 −3.2/ − 2.4 −4.3 −4.0 log g 4.0 4.0 3.5 4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 3.5 4.0 4.0 4.0 4.0 3.5/4.0 3.5/4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 4.0 W mdl (Hα ) 2.7 2.7 2.4 2.0 3.0 3.2 4.2 2.7 3.0 2.7 2.7 2.5 2.7 2.8 3.6 3.0 2.8 1.9 2.0/2.7 2.0/2.7 2.5 3.2 3.6 2.7/3.6 2.7 3.6/4.2 2.5 2.7 Notes. (a) Negueruela & Coe (2002), (b) Liu et al. (2000), (c) Coe et al. (2001), (d) Mc Bride et al. (2008), (e) Negueruela & Okazaki (2001), (f) Negueruela et al. (1996), (g) Reig & Roche (1999), (h) Janot-Pacheco et al. (1981), (i) Ziolkowski (2002), ( j) Wilson et al. (2003), (k) Negueruela et al. (2003), (l) Bonnet-Bidaud & Mouchet (1998), (m) Status as BeX still unknown. 2.5 Point 1 can be explained by the MWG Be/X-ray binaries location on the Galactic Plane, where the amount of interstellar material is large, while LMC and SMC sources are observed at high galactic latitudes and suffer from very low interstellar extinction. Etot(B−V) 2 Since the E is (B − V) is essentially constant for the LMC and SMC sources (but see below), points 1 and 2 strongly suggest that the lack of correlation for the MWG stars is entirely attributable to the interstellar reddening term. The relationship between E cs (B − V) and the amount of emission for MWG stars that should otherwise be apparent, is smeared out because E is (B − V) is totally independent for each star. In view of this, we will henceforth assume that the relation observed for LMC and SMC stars also holds for the MWG. 1.5 1 0.5 0 0 10 20 30 40 50 −W 60 70 80 o env (H )(A) α total Fig. 2. Relationship between E (B−V) and the W env (Hα ) for the Milky Way (red squares), LMC (blue asterisks) and SMC (green triangles) programme stars. Some stars have error bars of the order of the the symbol size. 2. the LMC and SMC stars occupy a narrow band in this plot, indicating a relationship. No apparent relation is observed for the MWG stars; 3. the SMC sources tend to have much stronger emissions than the LMC sources. Point 3 is surprising. Indeed, there are no strong Hα emitters within our LMC sample while there is only one mild emitter within our SMC sample. A dichotomy seems to exist around −W(Hα ) ≈ 15 Å. The MWG, however, displays the whole range of emissions. This could be caused by an observational bias owing to low number statistics. In order to explore this point further we also plotted the non simultaneous photometric and spectroscopic observations of Coe et al. (2005, open triangles) and Antoniou et al. (2009, open circles) in Fig. 3, left panel. Evidently, they occupy the same locus on average, although with a larger scatter. Few mild emitters seem to be present in the Coe et al. (2005) and Antoniou et al. (2009) samples. But because the number of sources is still very small, we dit not proceed with this question here. A114, page 5 of 13 A&A 539, A114 (2012) 0.3 0.25 0.5 0.2 E (B−V) 0.15 0.3 0.1 cs tot E (B−V) 0.4 0.05 0.2 0 0.1 −0.05 0 0 10 20 30 40 50 −W 60 70 80 o env −0.1 0 5 10 env (H )(A) −W α 15 20 o (Hα)(A) Fig. 3. Left panel: E total (B − V) vs. W env (Hα ), for LMC (blue circles) SMC (green triangles). Open triangles belong to SMC stars data from Coe et al. (2005) and spectral types from McBride et al. (2008). Open circles are data from Antoniou et al. (2009). The best fit is plotted for all stars with envelope emission lowerer than 15 Å. Right panel: E cs (B − V) vs. W env (Hα ) for stars with emission lower than 15 Å. Symbols as in left panel. The bottom line links two points that correspond to two different observations of 1A 0538−66. 3.3. E cs (B − V ) vs. W env (Hα ) If we represent only stars in the Clouds, (see Fig. 3, left panel), it is possible to appreciate some degree of correlation between both parameters for emission strengths −W env (Hα ) below ∼15 Å. Beyond this point, however, the relationship saturates. Then, as a first approximation, we performed a least-squares fit for stars with emission weaker than 15 Å and obtained E total (B − V) = −(0.011 ± 0.003)W env(Hα ) + (0.08 ± 0.03) r = 0.67, (7) where r represents the correlation coefficient. The zero point, 0.08 ± 0.02 mag, would correspond to the total excess without any emission at all, that is to say, the interstellar excess. This is, indeed, the average reddening of the Clouds if one takes the corresponding values in Table 5. This is encouraging. One of the stars from the LMC appears displaced upwards from the rest. This star, SAX J0544.1−7100, is in the southern part of the LMC where the interstellar reddening is higher (it is well known that a gradient exists in the LMC interstellar excess that increases towards the south, while the reddening across the SMC face seems to be more uniform). Indeed, according to Cusumano (1998), this source seems to have an even greater excess which would reach 0.32 mag. Consequently, we subtracted the individual E is (B − V) for each star. This excess has two contributions. On one hand there is the foreground interstellar reddening towards the Clouds. For this purpose we used the reddening maps of Schwering & Israel (1991). On the other hand there is the interstellar reddening within the Clouds. In this case we used the surveys performed by Zaritsky et al. (2004) for the LMC and Zaritsky et al. (2002) for the SMC. In this manner, we calculated the individual circumstellar excess for each star by means of E cs (B − V) = E total (B − V) − E is(B − V). Surprisingly with this approach, some of the resulting circumstellar excesses turned out to be negative, which means that we were over-correcting for reddening. Furthermore, the scatter in the relationship increases. As stated previously, the independent term in Eq. (7) is already very close to the interstellar value to the LMC. There are several reasons why introducing the intracloud reddening does not improve the correlation. As stated by Zaritsky et al. (2002, 2004), there are a number of caveats when using the estimated extinctions. For A114, page 6 of 13 example, the extinction values are only reliable if a number of stars are used to compute the average value. However, for the positions of our sources there is always a very high dispersion around the average value, sometimes ∼100%. Furthermore there seem to be effects due to the position of an object along the line-of-sight. All these effects contribute to increase the scatter in the correlations. Notwithstanding, the slope of the obtained correlations was similar within the errors. Finally, we used also the more recent and accurate reddening maps of the Clouds by Haschke et al. (2011). These maps have been obtained with two independent methods involving the red clump in colour magnitude diagrams and the absolute magnitudes of RR Lyrae stars, respectively. The obtained reddenings are less than 0.03 mag for the vast majority of our sources, much lower than the foreground value and well within the photometric accuracy of our measurements. Unfortunately, not all our sources are covered by the Haschke et al. (2011) survey. Therefore, we applied only the foreground correction. We present the E cs (B − V) vs. W env (Hα ) in Fig. 3, right panel. Obviously, the scatter of points is reduced and the degree of correlation is improved. We now obtain the following relation: E cs (B − V) = −(0.011 ± 0.003)W env(Hα ) r = 0.7. (8) The zero point in this equation should be strictly null, since both parameters are pure envelope contributions. Indeed, this is less than the photometric errors and, therefore, compatible with zero. On the other hand, the slope of the relation remains almost unchanged, indicating that subtraction of individual reddenings does not introduce any artificial trend. Finally, although there are a few points with EW(Hα ) ≥ −15Å, it is quite obvious that a saturation occurs around 0.17 ± 0.03 mag. The non contemporaneous data from Coe et al. (2005) and Antoniou et al. (2009) also support this conclusion. In Fig. 3, left panel, there are several stars with total excesses well above the mean. The majority of those stars indeed present intracloud reddening values of the order of ∼0.23 mag. In these data sets the photometry and the spectroscopy are not simultaneous. Therefore we will use below only stars within our sample. Taking into account the above discussion, we propose as the M. S. Riquelme et al.: Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way 2.0 1.8 0.25 1.6 Normalized Flux 0.3 Etot(B−V) 0.2 0.15 1.4 1.2 0.1 1.0 0.05 0.8 0.6 0 0.4 −0.05 6480 6500 6520 6540 6560 6580 6600 6620 6640 6660 o 0 10 20 30 40 50 Wavelength (A) 60 o −Wenv(Hα)(A) Fig. 5. Hα line shown by 1A0538−66 in two different epochs. Fig. 4. Correlation for Be/X-ray binaries (filled blue dots) compared with the similar correlation for isolated Be stars (red crosses). definitive calibration the following equations: E cs (B − V) = −(0.011 ± 0.003)W env (Hα ) if − W env (Hα ) < 15 Å ∼0.17 ± 0.03 if − W env (Hα ) ≥ 15 Å. (9) Comparing the first of these last equations with that valid for the isolated Be stars E cs (B − V) = −0.002W(Hα) (10) given by Dachs et al. (1988) (Eq. (39) of that reference), we see that the amount of circumstellar reddening in E cs (B − V) is about five times higher for Be stars in binaries with the same amount of emission3 , that is to say, the same value of W(Hα ). In the extensive study performed by Dachs et al. (1988) on a sample of isolated Be stars, the presence of saturation can also be seen in all infrared colours (Figs. 1 to 5 in that reference). However, this is not observed in the optical colour (Fig. 10 in that reference). The circumstellar excess in isolated Be stars, however, is always less than 0.1 mag even for the stronger emissions, while it is much higher for BeX stars (see Fig. 4). This lack of saturation in the excess at optical wavelengths is also confirmed by Fabregat & Torrejón (1998) using Strömgren photometry of (isolated) Be stars in open clusters (Fig. 3 of that reference). Therefore, there must exist a fundamental difference in the physical properties of circumstellar envelopes of Be stars in Be/X-ray binaries with respect to their isolated counterparts because the saturation is present here at higher photon energies. Finally, 1A 0538−66 was observed twice in different runs. In Fig. 5 we present the spectra of this source in two different states. During the second run the source showed an Hα line in absorption with an equivalent width consistent with the expected photospheric value, that is to say, the Be star was without envelope. The corresponding points are connected by a dotted line in Fig. 3. As can be seen a) it correctly sets the zero point of the calibration (within the photometric errors) and b) it follows the slope of the calibration almost exactly. This behaviour is encouraging and suggests that the calibration will be giving the corresponding value for each individual star and not only the statistical trend. 3 Although these authors used equivalent widths not corrected for photospheric absorption. 3.4. E total (V − I ) vs. W env (Hα ) As stated before, the nature of the circumstellar excess is such that it increases with wavelength. For heavily obscured systems, it is often useful to use the infrared colour excess because it diminishes with increasing wavelength. One of the indices that we can obtain from our sample is (R − I). However, the Hα line is located within the wavelength range of the R broadband filter. This results in contaminated values of R, rendering its use as a measure of the slope of the continuum very problematic. In view of this, we instead explored the relationship between (V − I) and the amount of emission. We applied the same procedure as described above for (B − V) to this near-infrared colour (V − I). The total excess is now defined as E total (V − I) = (V − I)obs − (V − I)0 , (11) where the intrinsic colour (V −I)0 is obtained from Ducati (2001) (listed here in Table 4) and (V − I)obs is obtained from our photometry (Table 1). In Fig. 6 we display the (V − I) total excesses, calculated in this way (Table 6), for all programme stars with measured I magnitude. Unfortunately, we were able to use only four stars from the LMC sample and none from the SMC (because no I is available). For the SMC we therefore added data from Coe et al. (2005). Using only LMC stars, we obtain the following equation: E total (V − I) = −(0.040 ± 0.005)W env(Hα ) + (0.19 ± 0.04). (12) Stars with the higher levels of emission, −W env (Hα ) > 15 Å, clearly present a saturation around 0.60 ± 0.10 mag. This can now be compared with the analogous relation for (B − V) (Eq. (8)). Evidently the infrared correlation presents a higher slope (about four times higher), as expected. This is consistent with the fact that the circumstellar excess increases towards the red. This statement is based on a very small number of data points and it is only suggestive. Again, the zero point of the relation should be the total excess for zero emission, that is to say, the interstellar reddening towards the LMC. From Eq. (12) this turns out to be E is (V − I) = 0.19 ± 0.04. In a totally independent way, we can also obtain this value using the relation between interstellar excesses at different wavelengths given by Fitzpatrick (1999), is is (V − I) = 1.6ELMC (B − V), ELMC (13) A114, page 7 of 13 A&A 539, A114 (2012) Table 5. Total, interstellar and circumstellar excesses for Magellanic Cloud programme stars as well as their Hα emissions and photosphere corrected Hα emissions. LMC SMC CAL E CAL 9 RX J0516.0−6916 RX J0516.0−6916 RX J0520.5−6932 RX J0529.8−6556 EXO 0531.1−6609 EXO 0531.1−6609 RX J0531.5−6518 RX J0531.5−6518 RX J0535.0−6700 1A 0538−66 1A 0538−66 SAX J0544.1−7100 H 0544−66 RX J0051.9−7311 RX J0052.1−7319 RX J0059.2−7138 RX J0101.0−7206 RX J0101.3−7211 E total (B − V) 0.16 ± 0.03 0.21 ± 0.03 0.15 ± 0.04 0.19 ± 0.03 0.30 ± 0.03 0.12 ± 0.03 0.17 ± 0.04 0.22 ± 0.03 0.18 ± 0.04 0.14 ± 0.03 0.22 ± 0.04 0.05 ± 0.03 0.15 ± 0.04 0.42 ± 0.04 0.10 ± 0.04 0.28 ± 0.03 0.30 ± 0.03 0.24 ± 0.03 0.23 ± 0.03 0.22 ± 0.03 E is (B − V) 0.09 0.12 0.11 0.11 0.08 0.07 0.08 0.08 0.06 0.06 0.10 0.09 0.09 0.32 0.09 0.07 0.07 0.10 0.10 0.10 E cs (B − V) 0.07 ± 0.03 0.09 ± 0.03 0.04 ± 0.04 0.08 ± 0.03 0.22 ± 0.03 0.05 ± 0.03 0.09 ± 0.04 0.14 ± 0.03 0.12 ± 0.04 0.08 ± 0.03 0.12 ± 0.04 0.00∗ ±0.03 0.06 ± 0.04 0.10 ± 0.04 0.01 ± 0.04 0.21 ± 0.03 0.23 ± 0.03 0.14 ± 0.03 0.13 ± 0.03 0.12 ± 0.03 W spc (Hα ) −5.6 −10.4 −3.0 −3.2 −11.0 −2.3 −6.6 −6.7 −10.2 −10.1 −5.2 2.3 −10.5 −5.3 −3.0 −14.6 −30.6 −23.2 −56.0 −10.7 W env (Hα ) −8.3 −13.1 −5.4∗∗ −5.6∗∗ −13.0 −5.3 −9.8 −9.9 −14.4 −14.3 −7.9 −0.7 −13.5∗∗ −8.0 −5.7 −17.1 −33.3 −26.0 −59.6 −13.7 ±δ(W env (Hα )) 0.2 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 1.0 0.3 0.5 0.1 0.9 0.5 0.5 0.5 1.0 0.2 Notes. (∗) In this particular case (lack of emission), E is (B − V) is greater than E total (B − V), so that E cs (B − V) = −0.04 becomes negative. This is unphysical. It simply reflects the photometric uncertainty in the zero point. In consequence, we set this point to zero. (∗∗) Shell line. along with the corresponding values of E is (B − V) from Table 5. The average value of E(V − I) derived from Eq. (13) turns out to be 0.13 ± 0.06, which is compatible with the previously obtained value of 0.19 ± 0.04 within the errors. Taking into account the previous discussion and a mean insterstellar excess of 0.13 for the SMC (obtained by applying Eq. (13) to the E is (B−V) of SMC stars taken from the Schwering & Israel 1991, reddening maps), we propose the following equation: −(0.040 ± 0.005)W env (Hα ) if − W env (Hα ) < 15 Å E cs (V − I) = env ∼0.47 ± 0.12 if − W (Hα ) ≥ 15 Å, Table 6. (V − I) excesses for LMC and Milky Way stars for those stars with photometry in I filter. Star LMC CAL E CAL 9 RX J0520.5−6932 RX J0529.8−6556 1A 0538−66 E tot (V − I) E cs (V − I) E is (V − I) 0.51 ··· 0.69 0.46 0.18 0.33 0.52 0.52 0.21 0.03 0.18 ··· 0.17 0.25 0.15 2.23 2.55/2.53 1.74 1.86/1.84 0.73 2.25 ··· 1.20 2.84 1.10 1.85 2.70 0.50 0.50 0.78 0.50 0.26 0.55 0.07 0.36 0.71 0.27 0.50 0.25 1.73 2.05/1.03 0.96 1.36/1.34 0.47 1.70 ··· 0.84 2.13 0.83 1.35 2.45 MWG (14) where only data from our sample were used. Again, since this is based on a very small number of data points, the result is only suggestive. The values of E cs (V − I) and E is (V − I) so obtained are shown in Table 6. These equations are based on a quite small and heterogeneous sample though. The values of the calculated circumstellar excesses in (V − I) must be taken, therefore, with due caution. Finally we note that the total excesses of MWG and LMC stars displayed in Fig. 6 are closer than those displayed in Fig. 2 for (B − V). This can be understood if we consider that total excesses in (V − I) will now have a higher contribution from the circumstellar component compared to the interstellar component (compare Eqs. (8) and (12)). 4. V absorption and emission Once E is (B − V) is known, we can estimate the V extinction by using AV = 3.1E is (B − V). This expression assumes a standard reddening law, which will not always hold (Fitzpatrick 1999). But it will almost exactly match the low reddenings of the Clouds. The value thus calculated, however, cannot be used to compute the distance modulus by direct comparison with the A114, page 8 of 13 SS 397 XTE 1946+274 KS 1947+300 X Cep X-4 4U 2206+54 SAX J2239.3+6116 4U 0115+634 3A 0726−260 GRO J1008−57 LS 1698 1A 1118−61 GX 304−1 Notes. The uncertainty is less than 0.06 mag for total excess and 0.10 mag for circumstellar and interstellar excesses. absolute magnitude for the corresponding spectral type of the star. The reason, of course, is that the V magnitude is also affected by the envelope. The classical Be stars brighten with increasing emission, that is to say, they are brighter when redder (Zorec & Briot 1991; Fabregat & Torrejón 1998). However, shell Be stars4 show the opposite trend (Harmanec 1983). Two stars in 4 These are believed to be Be stars viewed edge-on so that the circumstellar envelope self-absorbs the circumstellar emission. M. S. Riquelme et al.: Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way Replacing E cs (B − V) by Eq. (10) and EVcs by Eq. (16), we obtain Rcs = 1.7 for isolated Be stars. Assuming now that this factor also holds for Be/X-ray binaries, we would have 3 2.5 EVcs = 1.7(−0.011W env(Hα )), Etot(V−I) 2 1.5 1 0.5 0 0 10 20 30 o 40 −Wenv(H )(A) 50 60 α Fig. 6. E (V − I) vs. W (Hα ), for LMC (squares), Milky Way stars (triangles). Crosses belong to the SMC data of Coe et al. (2005) and spectral types from McBride et al. (2008). The line represents the relationship for LMC stars only. Some stars have error bars of the order of the the symbol size. tot env our LMC sample have shown the shell phenomenon in the spectrum5 in the past, RXJ 0516.0−6916 and 1A0538−66. We will restrict ourselves here to the classical Be stars and therefore will not consider these two stars in the following discussion, even though they presented non-shell spectra during our campaign. In principle, we could follow the same approach as for the colour excesses. We know the total extinction in V by means of Atot V = Vobs − V0 = V − V0 , where V0 can be derived from the known luminosity class and distance to the LMC and SMC. The net effect over the V magnitude is a combination of two contributions: Atotal = AisV + EVcs . V which, according to these authors, accounts for the excess emission from a circumstellar envelope caused by free-free and freebound transitions. In order to find a suitable relationship for Be/X-ray binaries, we considered a relationship of the form 5 where we made use of our calibration (Eq. (9)). As stated in the introduction, and in agreement with the results presented in previous sections, there is evidence that the physical properties of circumstellar envelopes in Be/X-ray binaries differ significantly from those of isolated Be stars. Therefore the previous assumption might not hold. However, as will be shown below in Sect. 5, only the above relation yields the correct values for the independently known distances to our sources, both in the Clouds and the MWG. For the sake of comparison, we recast the previous equation in the form env W (Hα ) EVcs ≈ −0.6 · (19) −30 Å Comparing this equation with Eq. (16) we see that the excess in V would be six times higher for Be stars in binaries. On the other hand, for stars with −W env (Hα ) > 15 Å, the saturation regime holds and E cs (V) = −1.7 · 0.17 = −0.3, (20) independent of the amount of emission. Therefore, we propose the following calibration, valid for Be/X-ray binaries, under the assumption that the factor Rcs is the same for isolated Be stars and Be/X-ray binaries: env 3.1E is(B − V) − 0.6 W−30(HÅα ) if − W env (Hα ) < 15 Å tot AV = if − W env (Hα ) ≥ 15 Å. 3.1E is(B − V) − 0.3 (21) (15) The two terms in Eq. (15) have opposite signs, positive for interstellar absorption and negative for the circumstellar emission term. The final brightness of the star depends on each particular combination of interstellar extinction and circumstellar brightening. Since the interstellar term is known a priori we can isolate the EVcs and compare it with the known values of −W env (Hα ), in the same way as we did with the (B − V) excess. Unfortunately, this procedure does not give any observable trend within our programme stars in the Clouds. In view of this, we will explore whether the value given by Dachs et al. (1988) (Eq. (40) in that reference) for isolated Be stars is valid for our objects. The circumstellar emission in V is related to the emission-line strength through the equation env W (Hα ) EVcs ≈ −0.1 , (16) −30 Å EVcs = Rcs E cs (B − V). (18) (17) Central absorption in the Hα emission line below the continuum level. 5. Distance modulus Next we evaluate the goodness (self-consistency) of the correction proposed in the preceding section by computing the distance to the Magellanic Clouds, using our programme stars, with three different methods: 1. The method commonly found in the literature, using the total excess AV = 3.1 E total (B − V). This method implicitly (and certainly incorrectly) assumes that the total-to-selective absorption in the V band has exactly the same expression (and sign) for the interstellar and circumstellar terms. 2. We compute the distance using only the correctly calculated interstellar value AV = 3.1 E is (B − V). This method would completely neglect the brightening of the Be star by the circumstellar emission. 3. Finally, we use the V excess determined by mean of Eq. (21), which correctly takes into account both effects, but assumes that the E cs (B − V) contribution to total emission valid for isolated Be stars also holds for Be/X-ray binary stars. Now, for each one of these methods, we computed the distance modulus with the commonly used expression: DM = V − MV − AV , (22) where V is the observed value taken from Table 1, MV is the absolute magnitude for the corresponding spectral type and luminosity classs and AV as proposed in the preceding section. The A114, page 9 of 13 A&A 539, A114 (2012) Table 7. Distance modulus for Magellanic Clouds stars calculated with the three methods explained in the text. Galaxy LMC mean σ SMC Star CAL E CAL 9 RX J0520.5−6932 RX J0529.8−6556 EXO 0531.1−6609 EXO 0531.1−6609 RX J0531.5−6518 RX J0531.5−6518 RX J0535.0−6700 SAX J0544.1−7100 RX J0051.9−7311 RX J0052.1−7319 RX J0059.2−7138 RX J0101.0−7206 RX J0101.3−7211 mean σ DM1 17.67 18.16 18.00 17.92 17.72 17.54 17.93 18.06 18.26 17.93 18.0 0.2 17.82 17.93 17.73 18.21 18.55 18.1 0.3 DM2 17.96 18.61 18.44 18.10 18.06 17.88 18.42 18.54 18.53 18.20 18.3 0.3 18.44 18.55 18.35 18.83 19.02 18.6 0.3 DM3 18.12 18.87 18.73 18.20 18.25 18.07 18.71 18.83 18.69 18.36 18.5 0.3 18.65 18.76 18.55 19.03 19.29 18.9 0.3 resulting individual values, mean values and standard deviation for LMC and SMC are quoted in Table 7, where DM1 corresponds to the first method, DM2 to the second and DM3 to the third. Obviously DM1 differs significantly from the other methods, while DM2 and DM3 are compatible within the errors. The next step was to compare the results obtained by all three methods with the distance to the Clouds obtained independently with dedicated surveys. Storm et al. (2004) obtained 18.48 ± 0.07 for LMC and 18.88 ± 0.13 for SMC based on metallicity-dependent Cepheids calibrations. These values excellently agree with those derived through other independent methods such as large samples of eclipsing binaries, which deliver DM = 18.91±0.03±0.1 for the SMC (Hilditch et al. 2005). There are two important caveats. On one hand these are average distances. Distance determinations from several eclipsing binaries in the LMC suggest that there is a considerable depth along the line of sight (Fitzpatrick et al. 2003). The SMC is also known to be very extended along the line of sight, with many estimates well in excess of 10 kpc (Mathewson et al. 1986; Crowl et al. 2001). In a recent survey Subramaniam & Subramaniam (2009) found that, indeed, the SMC has a greater depth than the LMC. More important, these authors have pointed out the possible existence of a bulge in the optical center of the SMC where both the extinction and the depth are larger than in outer areas. Within our SMC sample, two stars show DM larger than the average value of the Cloud. These are located in the outer areas. The other three, in turn, show DM values lower than the Cloud average. These are located near the optical centre of the Cloud. This can be understood if we consider the increased optical depth of the bulge, which would prevent us from observing sources behind it, much in the same way as the MWG bulge does. In turn, the sources located far from the optical centre of the Cloud can be observed at greater distances. This possible explanation is based on a very small number of sources, however, and is, therefore, very tentative. These effects can explain the scatter of σ = 0.3 mag obtained within our small sample and why the direct determination of the V excess, described at the beginning of the previous section, does not work. On the other hand, the luminosity classes of many stars in our sample have been determined without taking A114, page 10 of 13 into account the possible brightening of the Be star caused by the circumstellar emission. That is to say, a star classified as class III could actually be class V. Only two class III stars are included in our sample, one in the SMC and one in the LMC. The difference in absolute magnitude from luminosity class III to class V for the spectral types concerned here is ΔMV 0.4. This is compatible with the standard deviation given in Table 7. The proposed correction was based on the assumption that the ratio of circumstellar emission in V with respect to (B − V) is Rcs = 1.7, obtained from isolated Be stars. In order to check the validity of this assumption, we varied this factor in steps of 0.05 and recalculated the distance to the Clouds. Any departure from the assumed value resulted in a significant worsening of the computed distance. Furthermore, as will be shown in the next section, only this factor yields the consistent results for Milky Way sources whose distances have been computed by completely independent methods. As can be seen in Table 7, Method 1 (the one commonly used) gives completely inaccurate values. The second and third methods significantly improve the distance to the Clouds and they are compatible within the errors. However, it is only the third one which delivers the correct mean value. Therefore we propose method three as the the best one and we will use it to compute the distance to MWG binaries in the next section. 6. Milky Way stars Even though the extrapolation to the MWG of calibrations based on stars in other galaxies might not be entirely accurate, we will assume in the following that such an assumption is tenable. The fact that no Be/X-ray binaries are known to exist in any MWG open cluster leaves, however, little margin to assume otherwise if we want to move forward in this matter. In any case, this is a fundamental limitation that affects any calibration based on stars in the Clouds. It has been shown that the lower metallicity of the SMC probably has implications on its higher abundance of Be stars (Martayan et al. 2010) and on the onset of the Be phenomenon. However, once the Be star is formed, there seem to be no fundamental differences between Be stars in the SMC, the LMC and the MWG. Therefore we applied Eqs. (9) and 21 to the MWG sources of our sample. First we computed the circumstellar excesses in (B − V) from the W env (Hα ) from Eq. (9), as explained in previous sections. Then we subtracted this excess from the total and obtained the interstellar excess through Eq. (2). This excess, along with the photosphere-corrected circumstellar emission, allows to estimate Atotal by means of Eq. (21). V Finally we computed the distance through the commonly used equation: d = 10(V−MV −AV +5)/5 . tot (23) In Table 8 we present the results for the MWG programme stars. For the shell stars no distances were computed because Eqs. (21) are not applicable. For the stars with a remaining uncertainty on their spectral type, both possible distances are given. We now compare the computed distances with previous determinations. In Table 9 we show the distances obtained with the proposed calibration to those compiled from the literature. Evidently, our distances exceed those previously published. Throughout explicitly excluding the brightening of the source caused by circumstellar emission leads to a underestimation of the distance. In two cases, our distance determinations are much higher than those quoted in Raguzova & Popov (2005). In both M. S. Riquelme et al.: Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way Table 8. Parameters obtained for Milky Way stars. Shell stars have neither AV nor distance because they do not satisfy Eq. (19). Star E total (B − V) ± 0.03 W spc (Hα ) 4U 0115+63 1.71 1.0 3A 0726−260 0.83 −7.0 GRO J1008−57 1.96 −15.0 LS 1698 0.77 −4.0∗ 1A 1118−61 1.17 −68.0 GX 304-1 1.97 −3.0∗ SS 397 1.75 −38.0 XTE J1946+274 2.04/2.01 −42.0 KS 1947+300 1.23 −16.8 Cep X-4 1.29/1.26 −54.0 4U 2206+54 0.65 −4.0 SAX J2239.3+6116 1.81 −11.0 W env (Hα ) −1.8 ± 0.1 −8.9 ± 1.0 −17.0/−17.7 ± 0.5 −6.0/−6.7 ± 0.2∗ −70.5 ± 1.0 −6.2 ± 0.2∗ −41.6 ± 1.0 −44.7/−45.6 ± 2.0 −19.5 ± 0.5 −57.6/−58.2 ± 1.0 −6.5 ± 0.5 −13.7 ± 1.0 E cs (B − V) E is (B − V) Atotal ± 0.31 d (kpc) V 0.02 ± 0.01 1.69 5.2 7.0 ± 0.3 0.10 ± 0.03 0.73 ± 0.04 2.1 7.0 ± 0.5 0.17 ± 0.04 1.79 ± 0.05 5.2 9.7 ± 0.8/5.8 ± 0.5 0.07 ± 0.01 0.70 ··· ··· 0.17 ± 0.05 1.00 ± 0.06 2.8 5.2 ± 0.9 0.07 ± 0.01 1.90 ± 0.03 ··· ··· 0.17 ± 0.05 1.58 ± 0.06 4.6 1.5 ± 0.2 0.17 ± 0.05 1.87/1.84 ± 0.06 5.5/5.4 12.0 ± 1.7/8.7 ± 1.2 0.17 ± 0.04 1.06 ± 0.05 3.0 10.4 ± 0.9 0.17 ± 0.05 1.12/1.09 ± 0.06 3.2/3.1 7.9 ± 1.2/5.9 ± 0.9 0.07 ± 0.02 0.58 ± 0.04 1.7 2.9 ± 0.2 0.15 ± 0.04 1.66 ± 0.05 4.9 5.4 ± 0.5 Notes. (∗) Shell line. Table 9. Distances for MWG stars calculated using our proposed method compared to those previously published in the literature. Source KS1947+300 4U2206+54 SAXJ2239.3+61 3A0726-26 1A1118-61 Cep X-4 GROJ1008-57 XTE 1946+27 4U0115+63 d (this work) (kpc) 10.4 ± 0.9 2.9 ± 0.2 5.4 ± 0.5 7.0 ± 0.5 5.2 ± 0.9 7.9 ± 1.2/5.9 ± 0.9 9.7 ± 0.8/5.8 ± 0.5 12.0 ± 1.7/8.7 ± 1.2 7.0 ± 0.3 d (previous) (kpc) 10 2.5 4.4 6.0 5.0 3.8 2.0 9.5 ± 2.9 7-8 Ref. RP05 RP05 RP05 RP05 RP05 RP05 RP05 Wi03 NO01 Notes. RP05 (Raguzova & Popov 2005). Wi03 (Wilson et al. 2003). NO01 (Negueruela & Okazaki 2001). XTE 1946+274 1A 1118-61 KS 1947+300 Cep X-4 4U 2206+54 SS397 GRO J1008-57 cases, though, the discrepancy comes from the assumed parameters and has nothing to do with the method presented here. For example GRO J1008−57: the distance quoted in Raguzova & Popov (2005) was taken from a work by Haberl & Sasaki (2005). These last authors used, in turn, results from Macomb et al. (1994). These authors estimated a V = 13.9 ± 0.2 mag for the counterpart using a cross calibration of archive plates. Our more accurate observed value is V = 15.08 ± 0.02, which accounts for the difference. Our derived distance, however, agrees well with the estimate made by Petre & Gehrels (1994) for luminosity class V. Remarkably, this distance puts the source right on the outer Sagittarius arm and, together with the argument by Petre & Gehrels (1994), would argue in favour of a mainsequence counterpart. The case of 4U0115+63 is a fortunate one because this system was observed by Negueruela & Okazaki (2001) in a diskless state. These authors favour a 7 kpc distance and warn that this could be an upper limit since there could be still some residual emission present. Considering that we observed the source with significant Hα emission, the agreement is excellent, and constitutes an independent check of the validity of this work. Likewise, the distance to XTE J1946+27 has been estimated using completely independent methods based on other energy ranges. Wilson et al. (2003) estimated a distance of 9.5 ± 2.9 kpc based on a correlation between measured spin-up rate and X-ray flux. This value excellently agrees with ours. Figure 7 shows the distribution of our sources on the Galactic Plane. As can be seen, and as is expected, the binaries trace the SAX J2239+6116 4U 0115+63 3A 0726-260 Fig. 7. Sketch of the Milky Way (based on NASA Spitzer satellite data) showing the position of the Be/X-ray binaries of our sample. Sources with two possible spectral types are connected by a line. arms of the MWG, following the location of the OB star populations. Sources such as SS397, GRO J1008−57 (as V class) and 1A 1118−61 belong to the outer Sagittarius Arm. The group of SAX J2239.3+61, 4U 2206+54, 3A 0726−26, GRO J1008−57 (as III class), Cep X−4 (as B2) y XTE J1946+274 (as B1), trace the Perseus arm. Finally, XTE J1946+274 (as B0), KS1947+300, Cep X−4 (as B1) and 4U 0115+63, trace the outer Norma arm. 7. Summary and discussion We summarize here the main conclusions and equations: – The total reddening in the Milky Way sources is always much higher than for Magellanic Cloud sources, an effect attributable to the interstellar component. A114, page 11 of 13 A&A 539, A114 (2012) – The circumstellar contribution to the (B − V) colour can be computed by means of the equation E cs (B − V) = −0.011W env (Hα ) if − Wenv (Hα ) < 15 Å ∼ 0.17 if − Wenv (Hα ) > 15 Å. – This result confirms that E cs (B − V) is about five times higher for Be stars in binaries than for isolated Be stars with the same amount of emission (same −W env (Hα )). As we pointed out in the introduction, the circumstellar envelopes of Be stars in binaries are geometrically smaller and optically thicker than the envelopes in isolated Be stars. Therefore it is plausible that the electron density in the envelopes of the former is higher than in the latter. This has been pointed out by Zamanov et al. (2001). The circumstellar continuum emission in turn will also be higher in the former because the Bremsstrahlung emission coefficient jν ∼ γne ni is proporcional to the electron density, as is the the circumstellar reddening. – The study of the near-infrared colour (V − I) suggested a similar behaviour albeit with a steeper slope and saturation level. This agrees with the fact that circumstellar reddening increases towards the red (see for example Dachs et al. 1988). This conclusion is based on a very small number of data points and therefore must be regarded as tentative. – In order to compute the distance to a Be/X-ray binary (non shell Be star) in the MWG, we must include the circumstellar contribution to the observed V magnitude by means of Atotal = V env 3.1E is (B − V) − 0.6 W−30(HÅα ) if − W env (Hα ) < 15 Å if − W env (Hα ) ≥ 15 Å. 3.1E is (B − V) − 0.3 The first term is the commonly used correction for interstellar absorption while the second one takes into account the over luminosity of Be stars. Neglecting the increase in V brightness caused by circumstellar emission results in an underestimation of the distance (hence luminosity) for the Be/X-ray binaries. The importance of this effect depends on the combination of interstellar extinction and circumstellar emission for a particular source. It shows that the maximum value for the overluminosity for BeX with strong Hα emission is 0.3 mag. This last relationship was obtained assuming that the ratio Rcs = E cs (V)/E cs (B − V) = 1.7 deduced from isolated Be stars also holds for Be/X-ray binaries. However, this assumption seems to be justified because any departure from this ratio yields unacceptable values of the independently known distances to BeX stars in the Magellanic Clouds and the Milky Way Galaxy. Appendix A: Distance errors We calculate the errors in the distance for Milky Way stars, σd , showed in Table 8, using error propagation. We start with d = 10(V−MV −AV +5)/5 . Then we differentiate Eq. (A.1) σV σAtotal σd V = + ln 10, d 5 5 (A.1) (A.2) and solve for σd , σd = d (σV + σAtotal ) ln 10, V 5 A114, page 12 of 13 (A.3) comes from Eqs. (21) where σV is taken from Table 2 and σAtotal V for non-saturated stars by means of σAtotal (nosat) = 3.1σ(E is(B − V)) + V 0.6 σ(W env (Hα )), 30 (A.4) while for saturated stars the errors are σAtotal (sat) = 3.1 σ(E is (B − V)). V (A.5) In Eq. (A.4), the first term is σ(E (B − V)) = is σ(E total (B − V)) 2 + (σ(E cs (B − V)))2 , (A.6) where σ(E total (B − V)) is taken from Table 2 and σ(E cs (B − V)) comes from Eq. (9) in this way: σ(E cs (B − V)) = 0.011σ(W env(Hα )) + 0.0017W env(Hα ). (A.7) The first term in this equation takes into account the error when measuring W env (Hα ) in the spectrum, listed in Table 8. The second term corresponds to the error in the slope of the circumstellar excess calibration (Eq. (9)). Acknowledgements. This research has been partially supported by the Spanish Ministerio de Ciencia e Innovación, MICINN, under grant AYA2010-15431. The authors gratefully acknowledge the anonymous referee, whose comments and suggestions improved the content of the paper. References Antoniou, V., Hatzidimitriou, D., Zezas, A., & Reig P. 2009, ApJ, 707, 1080 Bonnet-Bidaud, J. M., & Mouchet, M. 1998, A&A, 332, L9 Coe, M. J., Negueruela, I., Buckley, D. A. H., Haigh, N. J., & Laycock, S. G. T. 2001, MNRAS, 324, 623 Coe, M. J., Edge, W. R. T., Galache, J. L., & McBride, V. A. 2005, MNRAS, 356, 502 Crowl, H. H., Sarajedini, A., Piatti, A. E., et al. 2001, AJ, 122, 220 Cusumano, G., Israel, G. L., Mannucci, F., et al. 1998, A&A, 337, 772 Dachs, J., Engels, D., & Kiehling, R. 1988, A&A, 194, 167 Draper, P. W., & Gray, N. 2000, Starlink User Note 214.7, R.A.L. Draper, P. W., Taylor, M., & Allan, A. 2000, Starlink User Note 139.12, R.A.L. Ducati, J. R., Bevilacqua, C. M., Rembold, S. B., & Ribeiro, D. 2001, ApJ, 558, 309 Fabregat, J., & Torrejón, J. M. 1998, A&A 332, 643 Fernández Barba D. 2005, Ph.D. Thesis, University of Barcelona, Dpto. Astronomy & Meteorology, Area de conocimiento: ciencias experimentales y matematicas Fitzpatrick, E. L. 1999, PASP, 111, 63 Fitzpatrick, E. L., Ribas, I., Guinan, E. F., et al. 2003, ApJ, 587, 685 Gies, D. R., Bagnuolo, W. G., Baines, E. K., et al. 2007, ApJ, 654, 527 Grundstrom, E. D., & Gies, D. R. 2006, ApJ, 651, L53 Harmanec, P. 1983, Hvar Obs. Bull., 7, 55 Haschke, R., Grebel, E. K., & Duffau, S. 2011, AJ, 141, 158 Hilditch, R. W., Howarth, I. D., & Harries, T. J. 2005, MNRAS, 357, 304 Howarth, I., Murray, J., Mills, D., & Berry, D. S. 1998, Starlink User Note 50.21, R.A.L. Janot-Pacheco, E., Ilovaisky, S. A., & Chevalier, C. 1981, A&A, 99, 274 Kurucz, R. L. 1979, ApJS, 40, 1 Landolt, A. U. 1995, VizieR On-line Data Catalog: II/118 Liu, Q. Z., van Paradijs, J., & van den Heuvel, E. P. J. 2000, A&AS, 147, 25 Macomb, D. J., Schrader, C. R., & Schultz, A. B. 1994, ApJ, 437, 845 Martayan, C., Baade, D., & Fabregat, J. 2010, A&A, 509, 11 Mathewson, D. S., Ford, V. L., & Visvanathan, N. 1986, ApJ, 301, 664 McBride, V. A., Coe, M. J., Negueruela, I., Schurch, M. P. E., & McGowan, K. E. 2008, MNRAS, 388, 1198 Negueruela, I. 1998, A&A, 338, 505 Negueruela, I., & Coe, M. J. 2002, A&A, 385, 517 Negueruela, I., & Okazaki, A. T. 2001, A&A, 369, 108 Negueruela, I., & Reig, P. 2001, A&A, 371, 1056 Negueruela, I., Roche, P., Buckley, D. A. H., et al. 1996, A&A, 315, 160 Negueruela, I., Israel, G. L., Marco, A., et al. 2003, A&A, 397, 739 M. S. Riquelme et al.: Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way Okazaki, A. T., & Negueruela, I. 2001, A&A, 337, 161 Petre, R., & Gehrels, N. 1994, A&A 282, L33 Raguzova, N. V., & Popov, S. B. 2005, Astron. Astrophys. Trans., 24, 151 Reig, P. 2011, Ap&SS 332, 1 Reig, P., & Roche, P. 1999, MNRAS, 306, 100R Reig, P., Fabregat, J., & Coe, M. 1997, A&A, 322, 193 Schmidt-Kaler, T. H. 1982, Physical parameters of the stars, in LandoltBornstein New Series, Vol. 2b, Astronomy and Astrophysics – stars and star clusters, ed. K. Schaifers, & H. H. Voigt (New York: Springer, Verlag) Schwering, P. B. W., & Israel, F. P. 1991, A&A, 246, 231 Shortridge, K., Meyerdicks, H., Currie, M., et al. 1997, Stalink User Note 86.15 R.A.L. Storm, J., Carney, B. W., Gieren, W. P., et al. 2004, A&A, 415, 531 Subramaniam, S., & Subramaniam, A. 2009, A&A, 496, 399S Torrejón, J. M., Negueruela, I., & Riquelme, M. S. 2005, in Active OB-stars: Laboratories for Stellar and Circumstellar Physics, Sapporo, Japan, 361, 503 Tycner, C., Gilnreath, G. C., Zavala, R. T., et al. 2006, AJ, 131, 2710 Vilchez, J. M., & Esteban, C. 1996, MNRAS, 280, 720 Wilson, C. A., Finger, M. H., Coe, M. J., & Negueruela, I. 2003, ApJ, 584, 996 Zamanov, R. K., Reig, P., Martí, J., et al. 2001, A&A, 367, 884 Zaritsky, D., Harris, J., Thompson, I. B., et al. 2002, AJ, 123, 855 Zaritsky, D., Harris, J., Thompson, I. B., et al. 2004, AJ, 128, 1606 Ziolkowski, J. 2002, Mem. Soc. Astron. Ital., 73, 1038 Zorec, J., & Briot, D. 1991, A&A, 245, 150 A114, page 13 of 13