N. Markova

Bright OB stars in the Galaxy - III. Constraints on the radial stratification of the clumping factor in hot star winds from a combined H

Context. Recent results strongly challenge the canonical picture of massive star winds: various evidence indicates that currently accepted mass-loss rates, u M, may need to be revised downwards, by factors extending to one magnitude or even more. This is because the most commonly used mass-loss diagnostics are affected by “clumping” (small-scale density inhomogeneities), influencing our interpretation of observed spectra and fluxes. Aims. Such downward revisions would have dramatic consequences for the evolution of, and feedback from, massive stars, and thus robust determinations of the clumping properties and mass-loss rates are urgently needed. We present a first attempt concerning this objective, by means of constraining the radial stratification of the so-called clumping factor. Methods. To this end, we have analyzed a sample of 19 Galactic O-type supergiants/giants, by combining our own and archival data for Hα, IR, mm and radio fluxes, and using approximate methods, calibrated to more sophisticated models. Clumping has been included into our analysis in the “conventional” way, by assuming the inter-clump matter to be void. Because (almost) all our diagnostics depends on the square of density, we cannot derive absolute clumping factors, but only factors normalized to a certain minimum. Results. This minimum was usually found to be located in the outermost, radio-emitting region, i.e., the radio mass-loss rates are the lowest ones, compared to u M derived from Hα and the IR. The radio rates agree well with those predicted by theory, but are only upper limits, due to unknown clumping in the outer wind. Hα turned out to be a useful tool to derive the clumping properties inside r < 3 ... 5R� . Our most important result concerns a (physical) difference between denser and thinner winds: for denser winds, the innermost region is more strongly clumped than the outermost one (with a normalized clumping factor of 4.1 ± 1.4), whereas thinner winds have similar clumping properties in the inner and outer regions. Conclusions. Our findings are compared with theoretical predictions, and the implications are discussed in detail, by assuming different scenarios regarding the still unknown clumping properties of the outer wind.

\mathsf{_{\alpha}}

Context. B-type supergiants represent an important phase in the evolution of massive stars. Reliable estimates of their stellar and wind parameters, however, are scarce, especially at mid and late spectral subtypes. Aims. We apply the NLTE atmosphere code FASTWIND to perform a spectroscopic study of a small sample of Galactic B-supergiants from B0 to B9. By means of the resulting data and incorporating additional datasets from alternative studies, we investigate the properties of OB-supergiants and compare our findings with theoretical predictions. Methods. Stellar and wind parameters of our sample stars are determined by line profile fitting, based on synthetic profiles, a Fourier technique to investigate the individual contributions of stellar rotation and “macro-turbulence” and an adequate approach to determine the Si abundances in parallel with micro-turbulent velocities. Results. Due to the combined effects of line- and wind-blanketing, the temperature scale of Galactic B-supergiants needs to be revised downward, by 10 to 20%, the latter value being appropriate for stronger winds. Compared to theoretical predictions, the wind properties of OB-supergiants indicate a number of discrepancies. In fair accordance with recent results, our sample indicates a gradual decrease in v∞ over the bi-stability region, where the limits of this region are located at lower Teff than those predicted. Introducing �

, IR and radio analysis

Context. The 30 Doradus (30 Dor) region of the Large Magellanic Cloud, also known as the Tarantula nebula, is the nearest starburst region. It contains the richest population of massive stars in the Local Group, and it is thus the best possible laboratory to investigate open questions on the formation and evolution of massive stars. Aims. Using ground-based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Dor. The sample is large enough to obtain statistically significant information and to search for variations among subpopulations – in terms of spectral type, luminosity class, and spatial location – in the field of view. Methods. We measured projected rotational velocities, ν_esini, by means of a Fourier transform method and a profile fitting method applied to a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, P(ν_e), of the equatorial rotational velocity, ν_e. Results. The distribution of νesini shows a two-component structure: a peak around 80 kms^(-1) and a high-velocity tail extending up to ~600 kms^(-1). This structure is also present in the inferred distribution P(ν_e) with around 80% of the sample having 0 < ν_e ≤ 300 kms^(-1) and the other 20% distributed in the high-velocity region. The presence of the low-velocity peak is consistent with what has been found in other studies for late O- and early B-type stars. Conclusions. Most of the stars in our sample rotate with a rate less than 20% of their break-up velocity. For the bulk of the sample, mass loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. If massive-star formation results in stars rotating at birth with a large portion of their break-up velocities, an alternative braking mechanism, possibly magnetic fields, is thus required to explain the present-day rotational properties of the O-type stars in 30 Dor. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Even though we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. That the high-velocity tail may be populated primarily (and perhaps exclusively) by post-binary interaction products has important implications for the evolutionary origin of systems that produce gamma-ray bursts.

Bright OB stars in the Galaxy IV. Stellar and wind parameters of early to late B supergiants

Detailed spectral classifications are presented for 352 O-B0 stars in the VLT-FLAMES Tarantula Survey ESO Large Programme, of which 213 O-type are judged of sufficiently high quality for further morphological analysis. Among them, six subcategories of special interest are distinguished. (1) Several new examples of the earliest spectral types O2-O3 have been found, while a previously known example has been determined to belong to the nitrogen-rich ON2 class. (2) A group of extremely rapidly rotating main-sequence objects has been isolated, including the largest vsini values known, the spatial and radial-velocity distributions of which suggest ejection from the two principal ionizing clusters NGC 2070 and NGC 2060. (3) Several new examples of the evolved, rapidly rotating Onfp class show similar evidence, although at least some of them are spectroscopic binaries. (4) No fewer than 48 members of the Vz category, hypothesized to be on or near the zero-age main sequence, are found in this sample; in contrast to the rapid rotators, they are strongly concentrated to the ionizing clusters and a newly recognized region of current and recent star formation to the north, supporting their interpretation as very young objects, as do their relatively faint absolute magnitudes. (5) A surprisingly large fraction of the main-sequence spectra belong to the recently recognized V((fc)) class, with C iii emission lines of similar strength to the usual N iii in V((f)) spectra, although a comparable number of the latter are also present, as well as six objects with very high-quality data but no trace of either emission feature, presenting new challenges to physical interpretations. (6) Two mid-O Vz and three late-O giant/supergiant spectra with morphologically enhanced nitrogen lines have been detected. Absolute visual magnitudes have been derived for each star with individual extinction laws, and composite Hertzsprung-Russell diagrams provide evidence of the multiple generations present in this field. Spectroscopic binaries, resolved visual multiples, and possible associations with X-ray sources are noted. Astrophysical and dynamical analyses of this unique dataset underway will provide new insights into the evolution of massive stars and starburst clusters.

The VLT-FLAMES Tarantula Survey - XII. Rotational velocities of the single O-type stars

We investigate the line-profile variability (lpv) of Hα for a large sample of O-type supergiants (15 objects between O4 and O9.7), in an objective, statistically rigorous manner. We employed the Temporal Variance Spectrum (TVS) analysis, developed for the case of photospheric absorption lines and modified by us to take into account the effects of wind emission. By means of a comparative analysis we place constraints on the properties of this variability - quantified in terms of a mean and a newly defined fractional amplitude of deviations - as a function of stellar and wind parameters. The results of our analysis show that all the stars in the sample show evidence of significant lpv in Hα, mostly dominated by processes in the wind. The variations occur between zero and 0.3 v∞ (i.e., below ∼1.5 R� ), in good agreement with results from similar studies. A comparison between the observations and corresponding line-profile simulations indicates that for stars with intermediate wind densities the properties of the Hα variability can be explained by simple models consisting of coherent or broken shells (blobs) uniformly distributed over the wind volume, with an intrinsic scatter in the maximum density contrast of about a factor of two. For stars at lower and higher wind densities, on the other hand, we found certain inconsistencies between the observations and our predictions, most importantly concerning the mean amplitude and the symmetry properties of the TVS. This disagreement might be explained by the presence of coherent large-scale structures, partly confined in a volume close to the star. Interpreted in terms of a variable mass-loss rate, the observed variations of Hα indicate changes of ±4% with respect to the mean value of u M for stars with stronger winds and of ± 16% for stars with weaker winds. The effect of these variations on the corresponding wind momenta is rather insignificant (less than 0.16 dex), increasing only the local scatter without affecting the Wind Momentum Luminosity Relationship.

The VLT-FLAMES Tarantula Survey. XIV. The O-type stellar content of 30 Doradus

Detailed studies of resolved young massive star clusters are necessary to determine their dynamical state and evaluate the importance of gas expulsion and early cluster evolution. In an effort to gain insight into the dynamical state of the young massive cluster R136 and obtain the first measurement of its velocity dispersion, we analyse multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud obtained as part of the VLT-FLAMES Tarantula Survey. Following a quantitative assessment of the variability, we use the radial velocities of non-variable sources to place an upper limit of 6 km s −1 on the line-of-sight velocity dispersion of stars within a projected distance of 5 pc from the centre of the cluster. After accounting for the contributions of undetected binaries and measurement errors through Monte Carlo simulations, we conclude that the true velocity dispersion is likely between 4a nd 5 km s −1 given a range of standard assumptions about the binary distribution. This result is consistent with what is expected if the cluster is in virial equilibrium, suggesting that gas expulsion has not altered its dynamics. We find that the velocity dispersion would be ∼25 km s −1 if binaries were not identified and rejected, confirming the importance of the multi-epoch strategy and the risk

Bright OB stars in the Galaxy II. Wind variability in O supergiants as traced by Hα

Context. Model atmosphere analyses have been previously undertaken for both Galactic and extragalactic B-type supergiants. By contrast, little attention has been given to a comparison of the properties of single supergiants and those that are members of multiple systems. Aims. Atmospheric parameters and nitrogen abundances have been estimated for all the B-type supergiants identified in the VLTFLAMES Tarantula survey. These include both single targets and binary candidates. The results have been analysed to investigate the role of binarity in the evolutionary history of supergiants. Methods. TLUSTY non-LTE (local thermodynamic equilibrium) model atmosphere calculations have been used to determine atmospheric parameters and nitrogen abundances for 34 single and 18 binary supergiants. E ective temperatures were deduced using the silicon balance technique, complemented by the helium ionisation in the hotter spectra. Surface gravities were estimated using Balmer line profiles and microturbulent velocities deduced using the silicon spectrum. Nitrogen abundances or upper limits were estimated from the Nii spectrum. The e ects of a flux contribution from an unseen secondary were considered for the binary sample. Results. We present the first systematic study of the incidence of binarity for a sample of B-type supergiants across the theoretical terminal age main sequence (TAMS). To account for the distribution of e ective temperatures of the B-type supergiants it may be necessary to extend the TAMS to lower temperatures. This is also consistent with the derived distribution of mass discrepancies, projected rotational velocities and nitrogen abundances, provided that stars cooler than this temperature are post-red supergiant objects. For all the supergiants in the Tarantula and in a previous FLAMES survey, the majority have small projected rotational velocities. The distribution peaks at about 50 km s 1 with 65% in the range 30 km s 1 ve sini 60 km s 1 . About ten per cent have larger ve sini ( 100 km s 1 ), but surprisingly these show little or no nitrogen enhancement. All the cooler supergiants have low projected rotational velocities of 70km s 1 and high nitrogen abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. Additionally, there are a lack of cooler binaries, possibly reflecting the small sample sizes. Single-star evolutionary models, which include rotation, can account for all of the nitrogen enhancement in both the single and binary samples. The detailed distribution of nitrogen abundances in the single and binary samples may be di erent, possibly reflecting di erences in their evolutionary history. Conclusions. The first comparative study of single and binary B-type supergiants has revealed that the main sequence may be significantly wider than previously assumed, extending to Te = 20 000 K. Some marginal di erences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolutionary status of blue supergiants.

The VLT-FLAMES Tarantula Survey. VII. A low velocity dispersion for the young massive cluster R136

Context. O Vz stars, a subclass of O-type dwarfs characterized by having He ii λ4686 stronger in absorption than any other helium line in their blue-violet spectra, have been suggested to be on or near the zero-age main sequence (ZAMS). If their youth were confirmed, they would be key objects with which to advance our knowledge of the physical properties of massive stars in the early stages of their lives. Aims. We test the hypothesis of O Vz stars being at a different (younger) evolutionary stage than are normal O-type dwarfs. Methods. We have performed the first comprehensive quantitative spectroscopic analysis of a statistically meaningful sample of O Vz and O V stars in the same star-forming region, exploiting the large number of O Vz stars identified by the VLT-FLAMES Tarantula Survey in the 30 Doradus region of the Large Magellanic Cloud (LMC). We obtained the stellar and wind parameters of 38 O Vz stars (and a control sample of 46 O V stars) using the FASTWIND stellar atmosphere code and the IACOB-GBAT, a grid-based tool developed for automated quantitative analysis of optical spectra of O stars. In the framework of a differential study, we compared the physical and evolutionary properties of both samples, locating the stars in the log g vs. log Teff, log Q vs. log Teff, and log L/L⊙ vs. log Teff diagrams. We also investigated the predictions of the FASTWIND code regarding the O Vz phenomenon. Results. We find a differential distribution of objects in terms of effective temperature, with O Vz stars dominant at intermediate values. The O Vz stars in 30 Doradus tend to be younger (i.e., closer to the ZAMS) and less luminous, and they have weaker winds than the O V stars, but we also find examples with ages of 2−4 Myr and with luminosities and winds that are similar to those of normal O dwarfs. Moreover, the O Vz stars do not appear to have higher gravities than the O V stars. In addition to effective temperature and wind strength, our FASTWIND predictions indicate how important it is to take other stellar parameters (gravity and projected rotational velocity) into account for correctly interpreting the O Vz phenomenon. Conclusions. In general, the O Vz stars appear to be on or very close to the ZAMS, but there are some examples where the Vz classification does not necessarily imply extreme youth. In particular, the presence of O Vz stars in our sample at more evolved phases than expected is likely a consequence of modest O-star winds owing to the low-metallicity environment of the LMC.

The VLT-FLAMES Tarantula Survey: XIX. B-type supergiants: Atmospheric parameters and nitrogen abundances to investigate the role of binarity and the width of the main sequence

Although it has important ramifications for both the formation of star clusters and their subsequent dynamical evolution, rotation remains a largely unexplored characteristic of young star clusters (few Myr). Using multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud obtained as part of the VLT-FLAMES Tarantula Survey, we search for rotation of the young massive cluster R136. From the radial velocities of 36 apparently single O-type stars within a projected radius of 10 pc from the centre of the cluster, we find evidence, at the 95% confidence level, for rotation of the cluster as a whole. We use a maximum likelihood method to fit simple rotation curves to our data and find a typical rotational velocity of ∼ 3k m s −1 . When compared to the low velocity dispersion of R136, our result suggests that star clusters may form with at least ∼20% of the kinetic energy in rotation.

The VLT-FLAMES Tarantula Survey XIII: On the nature of O Vz stars in 30 Doradus

Rotation is of key importance for the evolution of hot massive stars, however, the rotational velocities of these stars are difficult to determine. Based on our own data for 31 Galactic O stars and incorporating similar data for 86 OB supergiants from the literature, we aim at investigating the properties of rotational and extra line-broadening as a function of stellar parameters and at testing model predictions about the evolution of stellar rotation. Fundamental stellar parameters were determined by means of the code FASTWIND. Projected rotational and extra broadening velocities originate from a combined Ft + GOF method. Model calculations published previously were used to estimate the initial evolutionary masses. The sample O stars with Minit > 50 Msun rotate with less that 26% of their break-up velocity, and they also lack objects with v sin i 35 Msun on the hotter side of the bi-stability jump, the observed and predicted rotational rates agree quite well; for those on the cooler side of the jump, the measured velocities are systematically higher than the predicted ones. In general, the derived extra broadening velocities decrease toward cooler Teff, whilst for later evolutionary phases they appear, at the same v sin i, higher for high-mass stars than for low-mass ones. None of the sample stars shows extra broadening velocities higher than 110 km/s. For the majority of the more massive stars, extra broadening either dominates or is in strong competition with rotation. Conclusions: For OB stars of solar metallicity, extra broadening is important and has to be accounted for in the analysis. When appearing at or close to the zero-age main sequence, most of the single and more massive stars rotate slower than previously thought. Model predictions for the evolution of rotation in hot massive stars may need to be updated.