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.

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.

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

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.

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.

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.

He II lambda 4686 IN eta CARINAE: COLLAPSE OF THE WIND-WIND COLLISION REGION DURING PERIASTRON PASSAGE

The periodic spectroscopic events in eta Carinae are now well established and occur near the periastron passage of two massive stars in a very eccentric orbit. Several mechanisms have been proposed to explain the variations of different spectral features, such as an eclipse by the wind-wind collision (WWC) boundary, a shell ejection from the primary star or accretion of its wind onto the secondary. All of them have problems explaining all the observed phenomena. To better understand the nature of the cyclic events, we performed a dense monitoring of eta Carinae with five Southern telescopes during the 2009 low-excitation event, resulting in a set of data of unprecedented quality and sampling. The intrinsic luminosity of the He II lambda 4686 emission line (L similar to 310 L-circle dot) just before periastron reveals the presence of a very luminous transient source of extreme UV radiation emitted in the WWC region. Clumps in the primary's wind probably explain the flare-like behavior of both the X-ray and He II lambda 4686 light curves. After a short-lived minimum, He II lambda 4686 emission rises again to a new maximum, when X-rays are still absent or very weak. We interpret this as a collapse of the WWC onto the "surface" of the secondary star, switching off the hard X-ray source and diminishing the WWC shock cone. The recovery from this state is controlled by the momentum balance between the secondary's wind and the clumps in the primary's wind.

The young stellar cluster [DBS2003] 157 associated with the H ii region GAL 331.31-00.34

We report a study of the stellar content of the near-infrared (NIR) cluster [DBS2003] 157 embedded in the extended H ii region GAL 331.31-00.34, which is associated with the IRAS source 16085-5138. JHK photometry was carried out in order to identify potential ionizing candidates, and the follow-up NIR spectroscopy allowed the spectral classification of some sources, including two O-type stars. A combination of NIR photometry and spectroscopy data was used to obtain the distance of these two stars, with the method of spectroscopic parallax: IRS 298 (O6 V, 3.35 +/- 0.61 kpc) and IRS 339 (O9 V, 3.24 +/- 0.56 kpc). Adopting the average distance of 3.29 +/- 0.58 kpc and comparing the Lyman continuum luminosity of these stars with that required to account for the radio continuum flux of the H ii region, we conclude that these two stars are the ionizing sources of GAL 331.31-00.34. Young stellar objects (YSOs) were searched by using our NIR photometry and mid-infrared (MIR) data from the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) survey. The analysis of NIR and MIR colourcolour diagrams resulted in 47 YSO candidates. The GLIMPSE counterpart of IRAS 16085-5138, which presents IRAS colour indices compatible with an ultracompact H ii region, has been identified. The analysis of its spectral energy distribution between 2 and m revealed that this source shows a spectral index a= 3.6 between 2 and m, which is typical of a YSO immersed in a protostellar envelope. Lower limits to the bolometric luminosity and the mass of the embedded protostar have been estimated as L= 7.7 x 10(3) L? and M= 10 M?, respectively, which correspond to a B0B1 V zero-age main sequence star.

Color tunability of intense upconversion emission from Er3+-Yb3+ co-doped SiO2-Ta2O5 glass ceramic planar waveguides

This work reports on the infrared-to-visible CW frequency upconversion from planar waveguides based on Er3+-Yb3+-doped 100-xSiO(2)-xTa(2)O(5) obtained by a sol-gel process and deposited onto a SiO2-Si substrate by dip-coating. Surface morphology and optical parameters of the planar waveguides were analyzed by atomic force microscopy and the m-line technique. The influence of the composition on the electronic properties of the glass-ceramic films was followed by the band gap ranging from 4.35 to 4.51 eV upon modification of the Ta2O5 content. Intense green and red emissions were detected from the upconversion process for all the samples after excitation at 980 nm. The relative intensities of the emission bands around 550 nm and 665 nm, assigned to the H-2(11/2) -> I-4(15/2), S-4(3/2) -> I-4(15/2), and F-4(9/2) -> I-4(15/2) transitions, depended on the tantalum oxide content and the power of the laser source at 980 nm. The upconversion dynamics were investigated as a function of the Ta2O5 content and the number of photons involved in each emission process. Based on the upconversion emission spectra and 1931CIE chromaticity diagram, it is shown that color can be tailored by composition and pump power. The glass ceramic films are attractive materials for application in upconversion lasers and near infrared-to-visible upconverters in solar cells.

MHD numerical simulations of colliding winds in massive binary systems - I. Thermal versus non-thermal radio emission

In the past few decades detailed observations of radio and X-ray emission from massive binary systems revealed a whole new physics present in such systems. Both thermal and non-thermal components of this emission indicate that most of the radiation at these bands originates in shocks. O and B-type stars and WolfRayet (WR) stars present supersonic and massive winds that, when colliding, emit largely due to the freefree radiation. The non-thermal radio and X-ray emissions are due to synchrotron and inverse Compton processes, respectively. In this case, magnetic fields are expected to play an important role in the emission distribution. In the past few years the modelling of the freefree and synchrotron emissions from massive binary systems have been based on purely hydrodynamical simulations, and ad hoc assumptions regarding the distribution of magnetic energy and the field geometry. In this work we provide the first full magnetohydrodynamic numerical simulations of windwind collision in massive binary systems. We study the freefree emission characterizing its dependence on the stellar and orbital parameters. We also study self-consistently the evolution of the magnetic field at the shock region, obtaining also the synchrotron energy distribution integrated along different lines of sight. We show that the magnetic field in the shocks is larger than that obtained when the proportionality between B and the plasma density is assumed. Also, we show that the role of the synchrotron emission relative to the total radio emission has been underestimated.

Elemental abundances and classification of carbon-enhanced metal-poor stars

We present a detailed study of carbon-enhanced metal-poor (CEMP) stars, based on high-resolution spectroscopic observations of a sample of 18 stars. The stellar spectra for this sample were obtained at the 4.2 m William Herschel Telescope in 2001 and 2002, using the Utrecht Echelle Spectrograph, at a resolving power R similar to 52 000 and S/N similar to 40, covering the wavelength range lambda lambda 3700-5700 angstrom. The atmospheric parameters determined for this sample indicate temperatures ranging from 4750 K to 7100 K, log g from 1.5 to 4.3, and metallicities -3.0 <= [Fe/H]<=-1.7. Elemental abundances for C, Na, Mg, Sc, Ti, Cr, Cu, Zn, Sr, Y, Zr, Ba, La, Ce, Nd, Sm, Eu, Gd, Dy are determined. Abundances for an additional 109 stars were taken from the literature and combined with the data of our sample. The literature sample reveals a lack of reliable abundance estimates for species that might be associated with the r-process elements for about 67% of CEMP stars, preventing a complete understanding of this class of stars, since [Ba/Eu] ratios are used to classify them. Although eight stars in our observed sample are also found in the literature sample, Eu abundances or limits are determined for four of these stars for the first time. From the observed correlations between C, Ba, and Eu, we argue that the CEMP-r/s class has the same astronomical origin as CEMP-s stars, highlighting the need for a more complete understanding of Eu production.

The peculiar isolated neutron star in the Carina Nebula Deep XMM-Newton and ESO-VLT observations of 2XMM J104608.7-594306

While fewer in number than the dominant rotation-powered radio pulsar population, peculiar classes of isolated neutron stars (INSs) which include magnetars, the ROSAT-discovered "Magnificent Seven" (M7), rotating radio transients (RRATs), and central compact objects in supernova remnants (CCOs) - represent a key element in understanding the neutron star phenomenology. We report the results of an observational campaign to study the properties of the source 2XMM J104608.7-594306, a newly discovered thermally emitting INS. The evolutionary state of the neutron star is investigated by means of deep dedicated observations obtained with the XMM-Newton Observatory, the ESO Very Large Telescope, as well as publicly available gamma-ray data from the Fermi Space Telescope and the AGILE Mission. The observations confirm previous expectations and reveal a unique type of object. The source, which is likely within the Carina Nebula (N-H = 2.6x10(21) cm(-2)), has a spectrum that is both thermal and soft, with kT(infinity) = 135 eV. Non-thermal (magnetospheric) emission is not detected down to 1% (3 sigma, 0.1-12 keV) of the source luminosity. Significant deviations (absorption features) from a simple blackbody model are identified in the spectrum of the source around energies 0.6 keV and 1.35 keV. While the former deviation is likely related to a local oxygen overabundance in the Carina Nebula, the latter can only be accounted for by an additional spectral component, which is modelled as a Gaussian line in absorption with EW = 91 eV and sigma = 0.14 keV (1 sigma). Furthermore, the optical counterpart is fainter than m(V) = 27 (2 sigma) and no gamma-ray emission is significantly detected by either the Fermi or AGILE missions. Very interestingly, while these characteristics are remarkably similar to those of the M7 or the only RRAT so far detected in X-rays, which all have spin periods of a few seconds, we found intriguing evidence of very rapid rotation, P = 18.6ms, at the 4 sigma confidence level. We interpret these new results in the light of the observed properties of the currently known neutron star population, in particular those of standard rotation-powered pulsars, recycled objects, and CCOs. We find that none of these scenarios can satisfactorily explain the collective properties of 2XMM J104608.7-594306, although it may be related to the still poorly known class of Galactic anti-magnetars. Future XMM-Newton data, granted for the next cycle of observations (AO11), will help us to improve our current observational interpretation of the source, enabling us to significantly constrain the rate of pulsar spin down.