The VLT-FLAMES Tarantula survey: XX. The nature of the X-ray bright emission-line star VFTS 399

Context: The stellar population of the 30 Doradus star-forming region in the Large Magellanic Cloud contains a subset of apparently single, rapidly rotating O-type stars. The physical processes leading to the formation of this cohort are currently uncertain. 

Aims: One member of this group, the late O-type star VFTS 399, is found to be unexpectedly X-ray bright for its bolometric luminosity-in this study we aim to determine its physical nature and the cause of this behaviour. 

Methods: To accomplish this we performed a time-resolved analysis of optical, infrared and X-ray observations. 

Results: We found VFTS 399 to be an aperiodic photometric variable with an apparent near-IR excess. Its optical spectrum demonstrates complex emission profiles in the lower Balmer series and select He i lines-taken together these suggest an OeBe classification. The highly variable X-ray luminosity is too great to be produced by a single star, while the hard, non-thermal nature suggests the presence of an accreting relativistic companion. Finally, the detection of periodic modulation of the X-ray lightcurve is most naturally explained under the assumption that the accretor is a neutron star. 

Conclusions: VFTS 399 appears to be the first high-mass X-ray binary identified within 30 Dor, sharing many observational characteristics with classical Be X-ray binaries. Comparison of the current properties of VFTS 399 to binary-evolution models suggests a progenitor mass 25 M_⊙ for the putative neutron star, which may host a magnetic field comparable in strength to those of magnetars. VFTS 399 is now the second member of the cohort of rapidly rotating "single" O-type stars in 30 Dor to show evidence of binary interaction resulting in spin-up, suggesting that this may be a viable evolutionary pathway for the formation of a subset of this stellar population.

DYNAMICAL EVOLUTION OF VISCOUS DISKS AROUND Be STARS. I. PHOTOMETRY

Be stars possess gaseous circumstellar disks that modify in many ways the spectrum of the central B star. Furthermore, they exhibit variability at several timescales and for a large number of observables. Putting the pieces together of this dynamical behavior is not an easy task and requires a detailed understanding of the physical processes that control the temporal evolution of the observables. There is an increasing body of evidence that suggests that Be disks are well described by standard alpha-disk theory. This paper is the first of a series that aims at studying the possibility of inferring several disk and stellar parameters through the follow-up of various observables. Here we study the temporal evolution of the disk density for different dynamical scenarios, including the disk buildup as a result of a long and steady mass injection from the star, the disk dissipation that occurs after mass injection is turned off, as well as scenarios in which active periods are followed by periods of quiescence. For those scenarios, we investigate the temporal evolution of continuum photometric observables using a three-dimensional non-LTE radiative transfer code. We show that light curves for different wavelengths are specific of a mass loss history, inclination angle, and alpha viscosity parameter. The diagnostic potential of those light curves is also discussed.

Variable rotational line broadening in the Be star Achernar

Aims. The main theoretical problem for the formation of a Keplerian disk around Be stars is how angular momentum is supplied from the star to the disk, even more so since Be stars probably rotate somewhat subcritically. For instance, nonradial pulsation may transport angular momentum to the stellar surface until (part of) this excess supports the disk-formation/replenishment. The nearby Be star Achernar is presently building a new disk and o ers an excellent opportunity to observe this process from relatively close-up. Methods. Spectra from various sources and epochs are scrutinized to identify the salient stellar parameters characterizing the disk life cycle as defined by H emission. The variable strength of the non-radial pulsation is confirmed, but does not a ect the other results. Results. For the first time it is demonstrated that the photospheric line width does vary in a Be star, by as much as v sin i . 35 km However, unlike assumptions in which a photospheric spin-up accumulates during the diskless phase and then is released into the disk as it is fed, the apparent photospheric spin-up is positively correlated with the appearance of H line emission. The photospheric line widths and circumstellar emission increase together, and the apparent stellar rotation declines to the value at quiescence after the H line emission becomes undetectable

Probing the properties of Be star discs with spectroastrometry and NLTE radiative transfer modelling: beta CMi

While the presence of discs around classical Be stars is well established, their origin is still uncertain. To understand what processes result in the creation of these discs and how angular momentum is transported within them, their physical properties must be constrained. This requires comparing high spatial and spectral resolution data with detailed radiative transfer modelling. We present a high spectral resolution, R similar to 80 000, sub-milliarcsecond precision, spectroastrometric study of the circumstellar disc around the Be star beta CMi. The data are confronted with 3D, non-local thermodynamic equilibrium radiative transfer calculations to directly constrain the properties of the disc. Furthermore, we compare the data to disc models featuring two velocity laws: Keplerian, the prediction of the viscous disc model, and angular momentum conserving rotation. It is shown that the observations of beta CMi can only be reproduced using Keplerian rotation. The agreement between the model and the observed spectral energy distribution, polarization and spectroastrometric signature of beta CMi confirms that the discs around Be stars are well modelled as viscous decretion discs.

THE FIRST DETERMINATION OF THE VISCOSITY PARAMETER IN THE CIRCUMSTELLAR DISK OF A Be STAR

Be stars possess gaseous circumstellar decretion disks, which are well described using standard alpha-disk theory. The Be star 28 CMa recently underwent a long outburst followed by a long period of quiescence, during which the disk dissipated. Here we present the first time-dependent models of the dissipation of a viscous decretion disk. By modeling the rate of decline of the V-band excess, we determine that the viscosity parameter alpha = 1.0 +/- 0.2, corresponding to a mass injection rate (M) over dot = (3.5 +/- 1.3) x 10(-8) M-circle dot yr(-1). Such a large value of a suggests that the origin of the turbulent viscosity is an instability in the disk whose growth is limited by shock dissipation. The mass injection rate is more than an order of magnitude larger than the wind mass-loss rate inferred from UV observations, implying that the mass injection mechanism most likely is not the stellar wind, but some other mechanism.

Stochastic gravito-inertial modes discovered by CoRoT in the hot Be star HD51452

Context. Be stars are rapidly rotating stars with a circumstellar decretion disk. They usually undergo pressure and/or gravity pulsation modes excited by the kappa-mechanism, i.e. an effect of the opacity of iron-peak elements in the envelope of the star. In the Milky Way, p-modes are observed in stars that are hotter than or equal to the B3 spectral type, while g-modes are observed at the B2 spectral type and cooler. Aims. We observed a B0IVe star, HD51452, with the high-precision, high-cadence photometric CoRoT satellite and high-resolution, ground-based HARPS and SOPHIE spectrographs to study its pulsations in great detail. We also used the lower resolution spectra available in the BeSS database. Methods. We analyzed the CoRoT and spectroscopic data with several methods: CLEAN-NG, FREQFIND, and a sliding window method. We also analyzed spectral quantities, such as the violet over red (V/R) emission variations, to obtain information about the variation in the circumstellar environment. We calculated a stellar structure model with the ESTER code to test the various interpretation of the results. Results. We detect 189 frequencies of variations in the CoRoT light curve in the range between 0 and 4.5 c d(-1). The main frequencies are also recovered in the spectroscopic data. In particular we find that HD51452 undergoes gravito-inertial modes that are not in the domain of those excited by the kappa-mechanism. We propose that these are stochastic modes excited in the convective zones and that at least some of them are a multiplet of r-modes (i.e. subinertial modes mainly driven by the Coriolis acceleration). Stochastically excited gravito-inertial modes had never been observed in any star, and theory predicted that their very low amplitudes would be undetectable even with CoRoT. We suggest that the amplitudes are enhanced in HD51452 because of the very rapid stellar rotation. In addition, we find that the amplitude variations of these modes are related to the occurrence of minor outbursts. Conclusions. Thanks to CoRoT data, we have detected a new kind of pulsations in HD51452, which are stochastically excited gravito-inertial modes, probably due to its very rapid rotation. These modes are probably also present in other rapidly rotating hot Be stars.

The relationship between gamma Cassiopeiae's X-ray emission and its circumstellar environment

gamma Cas is the prototypical classical Be star and is recently best known for its variable hard X-ray emission. To elucidate the reasons for this emission, we mounted a multiwavelength campaign in 2010 centered around four XMM-Newton observations. The observational techniques included long baseline optical interferometry (LBOI) from two instruments at CHARA, photometry carried out by an automated photometric telescope and H alpha observations. Because gamma Cas is also known to be in a binary, we measured radial velocities from the H alpha line and redetermined its period as 203.55 +/- 0.20 days and its eccentricity as near zero. The LBOI observations suggest that the star's decretion disk was axisymmetric in 2010, has an system inclination angle near 45 degrees, and a larger radius than previously reported. In addition, the Be star began an "outburst" at the beginning of our campaign, made visible by a brightening and reddening of the disk during our campaign and beyond. Our analyses of the new high resolution spectra disclosed many attributes also found from spectra obtained in 2001 (Chandra) and 2004 (XMM-Newton). As well as a dominant hot (approximate to 14 keV) thermal component, the familiar attributes included: (i) a fluorescent feature of Fe K even stronger than observed at previous times; (ii) strong lines of N VII and Ne XI lines indicative of overabundances; and (iii) a subsolar Fe abundance from K-shell lines but a solar abundance from L-shell ions. We also found that two absorption columns are required to fit the continuum. While the first one maintained its historical average of 1 x 10(21) cm(-2), the second was very large and doubled to 7.4 x 10(23) cm(-2) during our X-ray observations. Although we found no clear relation between this column density and orbital phase, it correlates well with the disk brightening and reddening both in the 2010 and earlier observations. Thus, the inference from this study is that much (perhaps all?) of the X-ray emission from this source originates behind matter ejected by gamma Cas into our line of sight.

New activity in the large circumstellar disk of the Be-shell star 48 Librae

Aims. Spectroscopic, polarimetric, and high spectral resolution interferometric data covering the period 1995-2011 are analyzed to document the transition into a new phase of circumstellar disk activity in the classical Be-shell star 48 Lib. The objective is to use this broad data set to additionally test disk oscillations as the basic underlying dynamical process. Methods. The long-term disk evolution is described using the V/R ratio of the violet and red emission components of H alpha and Br gamma, radial velocities and profiles of He I and optical metal shell lines, as well as multi-band BVRI polarimetry. Single-epoch broad-band and high-resolution interferometric visibilities and phases are discussed with respect to a classical disk model and the given baseline orientations. Results. Spectroscopic signatures of disk asymmetries in 48 Lib vanished in the late nineties but recovered some time between 2004 and 2007, as shown by a new large-amplitude and long-duration V/R cycle. Variations in the radial velocity and line profile of conventional shell lines correlate with the V/R behavior. They are shared by narrow absorption cores superimposed on otherwise seemingly photospheric He I lines, which may form in high-density gas at the inner disk close to the photosphere. Large radial velocity variations continued also during the V/R-quiet years, suggesting that V/R may not always be a good indicator of global density waves in the disk. The comparison of the polarization after the recovery of the V/R activity shows a slight increase, while the polarization angle has been constant for more than 20 years, placing tight limits on any 3-D precession or warping of the disk. The broad H-band interferometry gives a disk diameter of (1.72 +/- 0.2) mas (equivalent to 15 stellar radii), position angle of the disk (50 +/- 9)degrees and a relatively low disk flattening of 1.66 +/- 0.3. Within the errors the same disk position angle is derived from polarimetric observations and from photocenter shifts across Br gamma. The high-resolution interferometric visibility and phase profiles show a double or even multiple-component structure. A preliminary estimate based on the size of the Br gamma emitting region indicates a large diameter for the disk (tens of stellar radii). Overall, no serious contradiction between the observations and the disk-oscillation model could be construed.

Stochastic gravito-inertial modes discovered by CoRoT in the hot Be star HD51452

Context. Be stars are rapidly rotating stars with a circumstellar decretion disk. They usually undergo pressure and/or gravity pulsation modes excited by the κ-mechanism, i.e. an effect of the opacity of iron-peak elements in the envelope of the star. In the Milky Way, p-modes are observed in stars that are hotter than or equal to the B3 spectral type, while g-modes are observed at the B2 spectral type and cooler. Aims. We observed a B0IVe star, HD51452, with the high-precision, high-cadence photometric CoRoT satellite and high-resolution, ground-based HARPS and SOPHIE spectrographs to study its pulsations in great detail. We also used the lower resolution spectra available in the BeSS database. Methods. We analyzed the CoRoT and spectroscopic data with several methods: Clean-NG, FreqFind, and a sliding window method. We also analyzed spectral quantities, such as the violet over red (V/R) emission variations, to obtain information about the variation in the circumstellar environment. We calculated a stellar structure model with the ESTER code to test the various interpretation of the results. Results. We detect 189 frequencies of variations in the CoRoT light curve in the range between 0 and 4.5 c d−1. The main frequencies are also recovered in the spectroscopic data. In particular we find that HD51452 undergoes gravito-inertial modes that are not in the domain of those excited by the κ-mechanism. We propose that these are stochastic modes excited in the convective zones and that at least some of them are a multiplet of r-modes (i.e. subinertial modes mainly driven by the Coriolis acceleration). Stochastically excited gravito-inertial modes had never been observed in any star, and theory predicted that their very low amplitudes would be undetectable even with CoRoT. We suggest that the amplitudes are enhanced in HD51452 because of the very rapid stellar rotation. In addition, we find that the amplitude variations of these modes are related to the occurrence of minor outbursts. Conclusions. Thanks to CoRoT data, we have detected a new kind of pulsations in HD51452, which are stochastically excited gravito-inertial modes, probably due to its very rapid rotation. These modes are probably also present in other rapidly rotating hot Be stars.

Non-LTE Monte Carlo radiative transfer. II. Nonisothermal solutions for viscous Keplerian disks

We discuss the basic hydrodynamics that determines the density structure of the disks around hot stars. Observational evidence supports the idea that these disks are Keplerian (rotationally supported) gaseous disks. A popular scenario in the literature, which naturally leads to the formation of Keplerian disks, is the viscous decretion model. According to this scenario, the disks are hydrostatically supported in the vertical direction, while the radial structure is governed by the viscous transport. This suggests that the temperature is one primary factor that governs the disk density structure. In a previous study we demonstrated, using three-dimensional non-LTE Monte Carlo simulations, that viscous Keplerian disks can be highly nonisothermal. In this paper we build on our previous work and solve the full problem of the steady state nonisothermal viscous diffusion and vertical hydrostatic equilibrium. We find that the self-consistent solution departs significantly from the analytic isothermal density, with potentially large effects on the emergent spectrum. This implies that nonisothermal disk models must be used for a detailed modeling of Be star disks.

TOWARD UNDERSTANDING THE B[e] PHENOMENON. IV. MODELING OF IRAS 00470+6429

FS CMa type stars are a recently described group of objects with the B[e] phenomenon which exhibits strong emission-line spectra and strong IR excesses. In this paper, we report the first attempt for a detailed modeling of IRAS 00470+6429, for which we have the best set of observations. Our modeling is based on two key assumptions: the star has a main-sequence luminosity for its spectral type (B2) and the circumstellar (CS) envelope is bimodal, composed of a slowly outflowing disklike wind and a fast polar wind. Both outflows are assumed to be purely radial. We adopt a novel approach to describe the dust formation site in the wind that employs timescale arguments for grain condensation and a self-consistent solution for the dust destruction surface. With the above assumptions we were able to satisfactorily reproduce many observational properties of IRAS 00470+6429, including the Hi line profiles and the overall shape of the spectral energy distribution. Our adopted recipe for dust formation proved successful in reproducing the correct amount of dust formed in the CS envelope. Possible shortcomings of our model, as well as suggestions for future improvements, are discussed.

TOWARD UNDERSTANDING THE B[e] PHENOMENON. III. PROPERTIES OF THE OPTICAL COUNTERPART OF IRAS 00470+6429

FS CMa type stars are a group of Galactic objects with the B[e] phenomenon. They exhibit strong emission-line spectra and infrared excesses, which are most likely due to recently formed circumstellar dust. The group content and identification criteria were described in the first two papers of the series. In this paper we report our spectroscopic and photometric observations of the optical counterpart of IRAS 00470+6429 obtained in 2003-2008. The optical spectrum is dominated by emission lines, most of which have P Cyg type profiles. We detected significant brightness variations, which may include a regular component, and variable spectral line profiles in both shape and position. The presence of a weak Li I 6708 angstrom line in the spectrum suggests that the object is most likely a binary system with a B2-B3 spectral-type primary companion of a luminosity log L/L(circle dot) = 3.9 +/- 0.3 and a late-type secondary companion. We estimate a distance toward the object to be 2.0 +/- 0.3 kpc from the Sun.

A cyclic bipolar wind in the interacting binary V 393 Scorpii

V393 Scorpii is a double periodic variable characterized by a relatively stable non-orbital photometric cycle of 253 d. Mennickent et al. argue for the presence of a massive optically thick disc around the more massive B-type component and describe the evolutionary stage of the system. In this paper, we analyse the behaviour of the main spectroscopic optical lines during the long non-orbital photometric cycle. We study the radial velocity of the donor determining its orbital elements and find a small but significant orbital eccentricity (e = 0.04). The donor spectral features are modelled and removed from the spectrum at every observing epoch using the light-curve model given by Mennickent et al. We find that the line emission is larger during eclipses and mostly comes from a bipolar wind. We also find that the long cycle is explained in terms of a modulation of the wind strength; the wind has a larger line and continuum emissivity at the high state. We report the discovery of highly variable chromospheric emission in the donor, as revealed by the Doppler maps of the emission lines Mg II 4481 and C I 6588. We discuss notable and some novel spectroscopic features like discrete absorption components, especially visible at blue depressed O I 7773 absorption wings during the second half-cycle, Balmer double emission with V/R curves showing 'Z-type' and 'S-type' excursions around secondary and main eclipses, respectively, and H beta emission wings extending up to +/- 2000 km s(-1). We also discuss possible causes for these phenomena and for their modulations with the long cycle.

Astrophysical parameters of LS 2883 and implications for the PSR B1259-63 gamma-ray binary

Only a few binary systems with compact objects display TeV emission. The physical properties of the companion stars represent basic input for understanding the physical mechanisms behind the particle acceleration, emission, and absorption processes in these so-called gamma-ray binaries. Here we present high-resolution and high signal-to-noise optical spectra of LS 2883, the Be star forming a gamma-ray binary with the young non-accreting pulsar PSR B1259-63, showing it to rotate faster and be significantly earlier and more luminous than previously thought. Analysis of the interstellar lines suggests that the system is located at the same distance as (and thus is likely a member of) Cen OB1. Taking the distance to the association, d = 2.3 kpc, and a color excess of E(B – V) = 0.85 for LS 2883 results in MV ≈ –4.4. Because of fast rotation, LS 2883 is oblate (R eq sime 9.7 R ☉ and R pole sime 8.1 R ☉) and presents a temperature gradient (T eq≈ 27,500 K, log g eq = 3.7; T pole≈ 34,000 K, log g pole = 4.1). If the star did not rotate, it would have parameters corresponding to a late O-type star. We estimate its luminosity at log(L */L ☉) sime 4.79 and its mass at M * ≈ 30 M ☉. The mass function then implies an inclination of the binary system i orb ≈ 23°, slightly smaller than previous estimates. We discuss the implications of these new astrophysical parameters of LS 2883 for the production of high-energy and very high-energy gamma rays in the PSR B1259-63/LS 2883 gamma-ray binary system. In particular, the stellar properties are very important for prediction of the line-like bulk Comptonization component from the unshocked ultrarelativistic pulsar wind.

NGC 7419 as a template for red supergiant clusters

Context. The open cluster NGC 7419 is known to contain five red supergiants and a very high number of Be stars. However, there are conflicting reports about its age and distance that prevent a useful comparison with other clusters. Aims. We intend to obtain more accurate parameters for NGC 7419, using techniques different from those of previous authors, so that it may be used as a calibrator for more obscured clusters. Methods. We obtained Strömgren photometry of the open cluster NGC 7419, as well as classification spectroscopy of ~20 stars in the area. We then applied standard analysis and classification techniques. Results. We find a distance of 4 ± 0.4 kpc and an age of 14 ± 2 Myr for NGC 7419. The main-sequence turn-off is found at spectral type B1, in excellent agreement. We identify 179 B-type members, implying that there are more than 1200 M⊙ in B stars at present. Extrapolating this to lower masses indicates an initial cluster mass of between 7000 and 10 000 M⊙, depending on the shape of the initial mass function. We find a very high fraction (≈40%) of Be stars around the turn-off, but very few Be stars at lower masses. We also report for the first time a strong variability in the emission characteristics of Be stars. We verified that the parameters of the red supergiant members can be used to obtain accurate cluster parameters. Conclusions. NGC 7419 is sufficiently massive to serve as a testbed for theoretical predictions and as a template to compare more obscured clusters. The distribution of stars above the main-sequence turn-off is difficult to accommodate with current evolutionary tracks. Though the presence of five red supergiants is marginally consistent with theoretical expectations, the high number of Be stars and very low number of luminous evolved B stars hint at some unknown physical factor that is not considered in current synthesis models.