The impact of accretion disc winds on the optical spectra of cataclysmic variables

Many high-state non-magnetic cataclysmic variables (CVs) exhibit blueshifted absorption or P-Cygni profiles associated with ultraviolet (UV) resonance lines. These features imply the existence of powerful accretion disc winds in CVs. Here, we use our Monte Carlo ionization and radiative transfer code to investigate whether disc wind models that produce realistic UV line profiles are also likely to generate observationally significant recombination line and continuum emission in the optical waveband. We also test whether outflows may be responsible for the single-peaked emission line profiles often seen in high-state CVs and for the weakness of the Balmer absorption edge (relative to simple models of optically thick accretion discs). We find that a standard disc wind model that is successful in reproducing the UV spectra of CVs also leaves a noticeable imprint on the optical spectrum, particularly for systems viewed at high inclination. The strongest optical wind-formed recombination lines are H alpha and He ii lambda 4686. We demonstrate that a higher density outflow model produces all the expected H and He lines and produces a recombination continuum that can fill in the Balmer jump at high inclinations. This model displays reasonable verisimilitude with the optical spectrum of RW Trianguli. No single-peaked emission is seen, although we observe a narrowing of the double-peaked emission lines from the base of the wind. Finally, we show that even denser models can produce a single-peaked H alpha line. On the basis of our results, we suggest that winds can modify, and perhaps even dominate, the line and continuum emission from CVs.

Self-consistent modelling of line-driven hot-star winds with Monte Carlo radiation hydrodynamics

Radiative pressure exerted by line interactions is a prominent driver of outflows in astrophysical systems, being at work in the outflows emerging from hot stars or from the accretion discs of cataclysmic variables, massive young stars and active galactic nuclei. In this work, a new radiation hydrodynamical approach to model line-driven hot-star winds is presented. By coupling a Monte Carlo radiative transfer scheme with a finite volume fluid dynamical method, line-driven mass outflows may be modelled self-consistently, benefiting from the advantages of Monte Carlo techniques in treating multiline effects, such as multiple scatterings, and in dealing with arbitrary multidimensional configurations. In this work, we introduce our approach in detail by highlighting the key numerical techniques and verifying their operation in a number of simplified applications, specifically in a series of self-consistent, one-dimensional, Sobolev-type, hot-star wind calculations. The utility and accuracy of our approach are demonstrated by comparing the obtained results with the predictions of various formulations of the so-called CAK theory and by confronting the calculations with modern sophisticated techniques of predicting the wind structure. Using these calculations, we also point out some useful diagnostic capabilities our approach provides. Finally, we discuss some of the current limitations of our method, some possible extensions and potential future applications.

Self-regulated Shocks in Massive Star Binary Systems

In an early-type, massive star binary system, X-ray bright shocks result from the powerful collision of stellar winds driven by radiation pressure on spectral line transitions. We examine the influence of the X-rays from the wind-wind collision shocks on the radiative driving of the stellar winds using steady-state models that include a parameterized line force with X-ray ionization dependence. Our primary result is that X-ray radiation from the shocks inhibits wind acceleration and can lead to a lower pre-shock velocity, and a correspondingly lower shocked plasma temperature, yet the intrinsic X-ray luminosity of the shocks, L X, remains largely unaltered, with the exception of a modest increase at small binary separations. Due to the feedback loop between the ionizing X-rays from the shocks and the wind driving, we term this scenario as self-regulated shocks. This effect is found to greatly increase the range of binary separations at which a wind-photosphere collision is likely to occur in systems where the momenta of the two winds are significantly different. Furthermore, the excessive levels of X-ray ionization close to the shocks completely suppress the line force, and we suggest that this may render radiative braking less effective. Comparisons of model results against observations reveal reasonable agreement in terms of log (L X/L bol). The inclusion of self-regulated shocks improves the match for kT values in roughly equal wind momenta systems, but there is a systematic offset for systems with unequal wind momenta (if considered to be a wind-photosphere collision).

Chemical abundances and winds of massive stars in M31: a B-type supergiant and a WC star in OB 10

We present high quality spectroscopic data for two massive stars in the OB 10 association of M31, OB 10-64 (B0 la) and OB 10-WRI (WC6). Medium resolution spectra of both stars were obtained using the ISIS spectrograph on the William Herschel Telescope. This is supplemented with Hubble Space Telescope STIS UV spectroscopy and Keck I HIRES data for OB 10-64. A non- local thermodynamic equilibrium (LTE) model atmosphere and abundance analysis for OB 10-64 is presented, indicating that this star has similar photospheric CNO, Mg and Si abundances to solar neighbourhood massive stars. A wind analysis of this early B-type supergiant reveals a mass-loss rate of (M)over dot = 1.6 x 10(-6) M-circle dot yr(-1), and v(infinity) = 1650 km s(-1). The corresponding wind momentum is in good agreement with the wind momentum-luminosity relationship found for Galactic early-B supergiants. Observations of OB 10-WRI are analysed using a non-LTE, line-blanketed code, to reveal approximate stellar parameters of log L/L-circle dot similar to 5.7, T-* - 75 kK, v(infinity) similar to 3000 km s(-1), (M)over dot/(M-circle dot yr(-1)) similar to 10(-4.3) adopting a clumped wind with a filling factor of 10 per cent. Quantitative comparisons are made with the Galactic WC6 star HD 92809 (WR23) revealing that OB 10-WR1 is 0.4 dex more luminous, though it has a much lower C/He ratio (similar to0.1 versus 0.3 for HD 92809). Our study represents the first detailed, chemical model atmosphere analysis for either a B-type supergiant or a Wolf- Rayet (WR) star in Andromeda, and shows the potential of how such studies can provide new information on the chemical evolution of galaxies and the evolution of massive stars in the local Universe.

The VLT-FLAMES survey of massive stars: stellar parameters and rotational velocities in NGC 3293, NGC 4755 and NGC 6611

An analysis is presented of VLT-FLAMES spectroscopy for three Galactic clusters, NGC3293, NGC4755 and NGC6611. Non-LTE model atmosphere calculations have been used to estimate effective temperatures (from either the helium spectrum or the silicon ionization equilibrium) and gravities (from the hydrogen spectrum). Projected rotational velocities have been deduced from the helium spectrum (for fast and moderate rotators) or the metal line spectrum (for slow rotators). The origin of the low gravity estimates for apparently near main sequence objects is discussed and is related to the stellar rotational velocity. The atmospheric parameters have been used to estimate cluster distances (which are generally in good agreement with previous determinations) and these have been used to estimate stellar luminosities and evolutionary masses. The observed Hertzsprung-Russell diagrams are compared with theoretical predictions and some discrepancies including differences in the main sequence luminosities are discussed. Cluster ages have been deduced and evidence for non-coeval star formation is found for all three of the clusters. Projected rotational velocities for targets in the older clusters, NGC3293 and NGC4755, have been found to be systematically larger than those for the field, confirming recent results in other similar age clusters. The distribution of projected rotational velocities are consistent with a Gaussian distribution of intrinsic rotational velocities. For the relatively unevolved targets in the older clusters, NGC3293 and NGC4755, the peak of the velocity distribution would be 250 km s(-1) with a full-width-half-maximum of approximately 180 km s(-1). For NGC6611, the sample size is relatively small but implies a lower mean rotational velocity. This may be evidence for the spin-down effect due to angular momentum loss through stellar winds, although our results are consistent with those found for very young high mass stars. For all three clusters we deduce present day mass functions with Gamma-values in the range of -1.5 to -1.8, which are similar to other young stellar clusters in the Milky Way.

The VLT-FLAMES survey of massive stars: mass loss and rotation of early-type stars in the SMC

We have studied the optical spectra of a sample of 31 O- and early B-type stars in the Small Magellanic Cloud, 21 of which are associated with the young massive cluster NGC 346. Stellar parameters are determined using an automated fitting method (Mokiem et al. 2005, A&A, 441, 711), which combines the stellar atmosphere code FASTWIND (Puls et al. 2005, A&A, 435, 669) with the genetic algorithm based optimisation routine PIKAIA (Charbonneau 1995, ApJS, 101, 309). Comparison with predictions of stellar evolution that account for stellar rotation does not result in a unique age, though most stars are best represented by an age of 1-3 Myr. The automated method allows for a detailed determination of the projected rotational velocities. The present day v(r) sin i distribution of the 21 dwarf stars in our sample is consistent with an underlying rotational velocity (v(r)) distribution that can be characterised by a mean velocity of about 160-190 km s(-1) and an effective half width of 100-150 km s(-1). The vr distribution must include a small percentage of slowly rotating stars. If predictions of the time evolution of the equatorial velocity for massive stars within the environment of the SMC are correct (Maeder & Meynet 2001, A&A, 373, 555), the young age of the cluster implies that this underlying distribution is representative for the initial rotational velocity distribution. The location in the Hertzsprung-Russell diagram of the stars showing helium enrichment is in qualitative agreement with evolutionary tracks accounting for rotation, but not for those ignoring vr. The mass loss rates of the SMC objects having luminosities of log L-star/L-circle dot greater than or similar to 5.4 are in excellent agreement with predictions by Vink et al. (2001, A&A, 369, 574). However, for lower luminosity stars the winds are too weak to determine. M accurately from the optical spectrum. Three targets were classified as Vz stars, two of which are located close to the theoretical zero-age main sequence. Three lower luminosity targets that were not classified as Vz stars are also found to lie near the ZAMS. We argue that this is related to a temperature effect inhibiting cooler from displaying the spectral features required for the Vz luminosity class.

The VLT-FLAMES survey of massive stars: atmospheric parameters and rotational velocity distributions for B-type stars in the Magellanic Clouds

Aims.We aim to provide the atmospheric parameters and rotational velocities for a large sample of O- and early B-type stars, analysed in a homogeneous and consistent manner, for use in constraining theoretical models. Methods: Atmospheric parameters, stellar masses, and rotational velocities have been estimated for approximately 250 early B-type stars in the Large (LMC) and Small (SMC) Magellanic Clouds from high-resolution VLT-FLAMES data using the non-LTE TLUSTY model atmosphere code. This data set has been supplemented with our previous analyses of some 50 O-type stars (Mokiem et al. 2006, 2007) and 100 narrow-lined early B-type stars (Hunter et al. 2006; Trundle et al. 2007) from the same survey, providing a sample of ~400 early-type objects. Results: Comparison of the rotational velocities with evolutionary tracks suggests that the end of core hydrogen burning occurs later than currently predicted and we argue for an extension of the evolutionary tracks. We also show that the large number of the luminous blue supergiants observed in the fields are unlikely to have directly evolved from main-sequence massive O-type stars as neither their low rotational velocities nor their position on the H-R diagram are predicted. We suggest that blue loops or mass-transfer binary systems may populate the blue supergiant regime. By comparing the rotational velocity distributions of the Magellanic Cloud stars to a similar Galactic sample, we find that (at 3s confidence level) massive stars (above 8 M?) in the SMC rotate faster than those in the solar neighbourhood. However there appears to be no significant difference between the rotational velocity distributions in the Galaxy and the LMC. We find that the v sin i distributions in the SMC and LMC can modelled with an intrinsic rotational velocity distribution that is a Gaussian peaking at 175 km s-1 (SMC) and 100 km s-1 (LMC) with a 1/e half width of 150 km s-1. We find that in NGC 346 in the SMC, the 10-25 M? main-sequence stars appear to rotate faster than their higher mass counterparts. It is not expected that O-type stars spin down significantly through angular momentum loss via stellar winds at SMC metallicity, hence this could be a reflection of mass dependent birth spin rates. Recently Yoon et al. (2006) have determined rates of GRBs by modelling rapidly rotating massive star progenitors. Our measured rotational velocity distribution for the 10-25 M? stars is peaked at slightly higher velocities than they assume, supporting the idea that GRBs could come from rapid rotators with initial masses as low as 14 M? at low metallicities.

The empirical metallicity dependence of the mass-loss rate of O- and early B-type stars

We present a comprehensive study of the observational dependence of the mass-loss rate in stationary stellar winds of hot massive stars on the metal content of their atmospheres. The metal content of stars in the Magellanic Clouds is discussed, and a critical assessment is given of state-of-the-art mass-loss determinations of OB stars in these two satellite systems and the Milky-Way. Assuming a power-law dependence of mass loss on metal content,. M. Z(m), and adopting a theoretical relation between the terminal flow velocity and metal content, v(infinity). Z(0.13) (Leitherer et al. 1992, ApJ, 401, 596), we find m = 0.83 +/- 0.16 for non-clumped outflows from an analysis of the wind momentum luminosity relation (WLR) for stars more luminous than 105.2 L circle dot. Within the errors, this result is in agreement with the prediction m = 0.69 +/- 0.10 by Vink et al. (2001, A& A, 369, 574). Absolute empirical values for the mass loss, based on Ha and ultraviolet (UV) wind lines, are found to be a factor of two higher than predictions in this high luminosity regime. If this difference is attributed to inhomogeneities in the wind, and this clumping does not impact the predictions, this would imply that luminous O and early-B stars have clumping factors in their Ha and UV line forming regions of about a factor of four. For lower luminosity stars, the winds are so weak that their strengths can generally no longer be derived from optical spectral lines (essentially Ha) and one must currently rely on the analysis of UV lines. We confirm that in this low-luminosity domain the observed Galactic WLR is found to be much steeper than expected from theory (although the specific sample is rather small), leading to a discrepancy between UV mass-loss rates and the predictions by a factor 100 at luminosities of L similar to 10(4.75) L circle dot, the origin of which is unknown. We emphasize that even if the current mass-loss rates of hot luminous stars would turn out to be overestimated as a result of wind clumping, but the degree of clumping would be rather independent of metallicity, the scalings derived in this study are expected to remain correct.

Blue luminous stars in nearby galaxies: Quantitative spectral analysis of M33 B-type supergiant stars

We present the detailed spectral analysis of a sample of M33 B-type supergiant stars, aimed at the determination of their fundamental parameters and chemical composition. The analysis is based on a grid of non-LTE metal line-blanketed model atmospheres including the effects of stellar winds and spherical extension computed with the code FASTWIND. Surface abundance ratios of C, N, and O are used to discuss the chemical evolutionary status of each individual star. The comparison of observed stellar properties with theoretical predictions of massive star evolutionary models shows good agreement within the uncertainties of the analysis. The spatial distribution of the sample allows us to investigate the existence of radial abundance gradients in the disk of M33. The comparison of stellar and H II region O abundances ( based on direct determinations of the electron temperature of the nebulae) shows good agreement. Using a simple linear radial representation, the stellar oxygen abundances result in a gradient of -0.0145 +/- 0.005 dex arcmin(-1) (or -0.06 +/- 0.02 dex kpc(-1)) up to a distance equal to similar to 1.1 times the isophotal radius of the galaxy. A more complex representation cannot be completely discarded by our stellar sample. The stellar Mg and Si abundances follow the trend displayed by O abundances, although with shallower gradients. These differences in gradient slope cannot be explained at this point. The derived abundances of the three alpha-elements yield solar metallicity in the central regions of the disk of M33. A comparison with recent planetary nebula data from Magrini and coworkers indicates that the disk of M33 has not suffered from a significant O enrichment in the last 3 Gyr.

Rotational velocities of single and binary O-type stars in the Tarantula Nebula

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, υ e sin i, of ~330 O-type objects, i.e. ~210 spectroscopic single stars and ~110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30 Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the υ e sin i distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100 km s^-1. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300 km s^-1, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part - and could potentially be completely due - to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences.

High-resolution Keck I spectroscopy of Galactic halo post-asymptotic giant branch stars

Absolute and differential abundance analyses have been performed from high-resolution, high signal-to-noise ratio optical (Keck I) spectra for three evolved Galactic halo stars, namely PG 1704 + 222, HD 341617 and LSIV -0401. Their derived atmospheric parameters indicate that all three objects are undergoing a post-asymptotic giant branch (post-AGB) phase of evolution. A differential abundance analysis reveals HD 341617 as having a mild carbon deficiency of 0.74 dex, possibly due to the star having evolved off the AGB before the onset of the third dredge-up. Although such carbon underabundances are typical of hot post-AGB objects, the same trend is not observed in PG 1704 + 222, where the carbon abundance is found to be consistent with those derived for nitrogen and oxygen. Hence, a dredge-up scenario need not be invoked to explain the chemical composition of PG 1704 + 222. For LSIV -0401 no iron deficiency is apparent relative to magnesium and silicon, and hence a gas- dust separation event in the AGB progenitor need not be invoked for this star.

CaII K interstellar observations towards early-type disc and halo stars

We present high-resolution (R = lambda/Deltalambda similar to 40 000) Ca II K interstellar observations (lambda(air) = 3933.66Angstrom) towards 88 mainly B-type stars, of which 74 are taken from the Edinburgh-Cape or Palomar-Green surveys, and 81 have > 25degrees. The majority of the data come from previously existing spectroscopy, although also included are 18 new observations of stars with echelle spectra taken with UVES on the Very Large Telescope UT2 (Kueyen). Some 49 of the sample stars have distance estimates above the Galactic plane (z) greater than or equal to 1 kpc, and are thus good probes of the halo interstellar medium. Of the 362 interstellar Ca K components that we detect, 75 (21 per cent) have absolute values of their LSR velocity values exceeding 40 km s(-1). In terms of the deviation velocity for the sightlines with distance estimates, 46/273 (17 per cent) of components have velocity values exceeding those predicted by standard Galactic rotation by more than 40 km s(-1). Combining this data set with previous observations, we find that the median value of the reduced equivalent width (REW) of stars with z greater than or equal to 1 kpc (EW x sin ) is similar to 115 mAngstrom (n = 80), similar to that observed in extragalactic sightlines by Bowen. Using data of all z distances, the REW at infinity is found to be similar to 130 mAngstrom, with the scaleheight (1) of the Ca II K column density distribution being;z 800 pc (n = 196) and reduced column density at infinity of log[N(Ca II K) cm(-2)] similar to 12.24. This implies that similar to30 per cent of Ca II K absorption occurs at distances exceeding similar to1 kpc. For nine sightlines, with distance exceeding 1 kpc and with a companion object within 5degrees, we find that all but two have values of Ca II reduced equivalent width the same to within similar to20 per cent, when the REW of the nearest object is extrapolated to the distance of the further of the pair, and assuming 1 = 800 pc. For 29 of our sightlines with z greater than or equal to 1 kpc and a H I detection from the Leiden-Dwingeloo survey (beamsize of 0.5degrees), we find log(N(Ca II K)IN(H I)) ranging from -7.4 to - 8.4. Values of the Ca II K abundance relative to neutral hydrogen (log[N(Ca II K) cm(-2)] - log[N(H I) cm(-2)]) are found to be more than similar to0.5 dex higher in stars with distances exceeding approximate to100 pc, when compared with the (log[N(Ca II K) cm(-2)] -log[N(H-tot) cm(-2)]) values found in nearby sightlines such as those in Wakker & Mathis (2000). Finally, stellar Ca II K equivalent widths of the sample are determined for 26 objects.

CaII K interstellar observations towards early-type disc and halo stars - distances to intermediate- and high-velocity clouds

We compare existing high spectral resolution (R = lambda/Deltalambda similar to 40 000) Ca II Kobservations (lambda(air) = 3933.66 Angstrom) towards 88 mainly B-type stars, and new observations taken using the Intermediate dispersion Spectrograph and Imaging System (ISIS) on the William Herschel Telescope at R similar to 10 000 towards three stars taken from the Palomar-Green Survey, with 21-cm HI emission-line profiles, in order to search for optical absorption towards known intermediate- and high-velocity cloud complexes. Given certain assumptions, limits to the gas phase abundance of Ca II are estimated for the cloud components. We use the data to derive the following distances from the Galactic plane (z). (i) Tentative lower z-height limits of 2800 and 4100 pc towards complex C using lack of absorption in the spectra of HD341617 and PG 0855 + 294, respectively. (ii) A weak lower z-height of 1400 pc towards complex WA-WB using lack of absorption in EC 09470-1433 and a weak lower limit of 2470 pc using lack of absorption in EC 09452-1403. (iii) An upper z- height of 2470 pc towards a southern intermediate- velocity cloud (IVC) with v(LSR) = -55 km s(-1) using PG 2351 + 198. (iv) Detection of a possible IVC in Ca II absorption at v(LSR) = +52 km s(-1) using EC 20104-2944. No associated HI in emission is detected. At this position, normal Galactic rotation predicts velocities of up to similar to+ 25 km s(-1). The detection puts an upper z-height of 1860 pc to the cloud. (v) Tentative HI and Ca II K detections towards an IVC at similar to+70 km s(-1) in the direction of high-velocity cloud (HVC) complex WE, sightline EC 06387-8045, indicating that the IVC may be at a z-height lower than 1770 pc. (vi) Detection of Ca II K absorption in the spectrum of PG 0855 + 294 in the direction of IV20, indicating that this IVC has a z-height smaller than 4100 pc. (vii) A weak lower z-height of 4300 pc towards a small HVC with v(LSR) = +115 km s(-1) at l, b = 200degrees, + 52degrees, using lack of absorption in the Ca II K spectrum of PG 0955 + 291.