Source author record

J. Puls

J. Puls appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

44works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2022arXiv

Metallicity Estimation of Core-Collapse Supernova HII Regions in Galaxies within 30 Mpc

This work presents measurements of the local HII environment metallicities of core-collapse supernovae (SNe) within a luminosity distance of 30 Mpc. 76 targets were observed at the Isaac Newton Telescope and environment metallicities could be measured for 65 targets using the N2 and O3N2 strong emission line method. The cumulative distribution functions (CDFs) of the environment metallicities of Type Ib and Ic SNe tend to higher metallicity than Type IIP, however Type Ic are also present at lower metallicities whereas Type Ib are not. The Type Ib frequency distribution is narrower (standard deviation $\sim$0.06 dex) than the Ic and IIP distributions ($\sim$0.15 dex) giving some evidence for a significant fraction of single massive progenitor stars; the low metallicity of Type Ic suggests a significant fraction of compact binary progenitors. However, both the Kolmogorov-Smirnov test and the Anderson-Darling test indicate no statistical significance for a difference in the local metallicities of the three SN types. Monte-Carlo simulations reveal a strong sensitivity of these tests to the uncertainties of the derived metallicities. Given the uncertainties of the strong emission methods, the applicability of the tests seems limited. We extended our analysis with the data of the Type Ib/Ic/IIP SN sample from Galbany et al. (2018). The CDFs created with their sample confirm our CDFs very well. The statistical tests, combining our sample and the Galbany et al. (2018) sample, indicate a significant difference between Type Ib and Type IIP with <5% probability that they are drawn from the same parent population.

preprint2022arXiv

New predictions for radiation-driven, steady-state mass-loss and wind-momentum from hot, massive stars III. Updated mass-loss rates for stellar evolution

Massive stars lose a large fraction of their mass to radiation-driven winds throughout their entire life. These outflows impact both the life and death of these stars and their surroundings. Theoretical mass-loss rates of hot, massive stars are derived to be used in applications such as stellar evolution. The behaviour of these rates in the OB-star regime is analysed, and their effects on massive-star evolution predictions is studied. Dynamically-consistent models are computed by solving the spherically symmetric, steady-state equation-of-motion for a large grid of hot, massive stars. The radiative acceleration is derived from non-local thermodynamic equilibrium radiative transfer in the co-moving frame. The resulting mass-loss rates are used to derive a simple scaling recipe with stellar parameters and to evaluate some first impacts upon evolution tracks. We provide a new prescription for steady-state, radiation-driven mass-loss from hot, massive stars depending on their fundamental parameters. The rates for O-stars are lower by about a factor ~3 than the rates typically used in previous stellar-evolution calculations, where differences generally decrease with increasing luminosity and temperature. For cooler B-giants/supergiants we find larger discrepancies. This arises because we do not find any systematic increase in mass-loss rates below the so-called bi-stability region; indeed, our results do not show any sign of a significant bi-stability jump within the parameter range covered by the grid. Due to the lower mass-loss rates we find that envelopes are not easily stripped by means of standard steady-state winds, making it difficult to create classical Wolf-Rayet stars via this channel. However, a remaining key uncertainty regarding these predictions concerns unsteady mass loss for very high-luminosity stars close to the Eddington limit as well as the impact of non-line-driven winds.

preprint2022arXiv

The nature of the Cygnus extreme B-supergiant 2MASS J20395358+4222505

2MASS J20395358+4222505 is an obscured early B supergiant near the massive OB star association Cyg OB2. Despite its bright infrared magnitude (K$_{s}$=5.82) it has remained largely ignored because of its dim optical magnitude (B=16.63, V=13.68). In a previous paper we classified it as a highly reddened, potentially extremely luminous, early B-type supergiant. We obtained its spectrum in the U, B and R spectral bands during commissioning observations with the instrument MEGARA@GTC. It displays a particularly strong H$α$ emission for its spectral type, B1 Ia. The star seems to be in an intermediate phase between super- and hypergiant, a group that it will probably join in the near (astronomical) future. We observe a radial velocity difference between individual observations and determine the stellar parameters, obtaining T$_{eff}$ = 24000 K, logg$_{c}$= 2.88 $\pm$ 0.15. The rotational velocity found is large for a B-supergiant, vsini= 110 $\pm$ 25 km s$^{-1}$. The abundance pattern is consistent with solar, with a mild C underabundance (based on a single line). Assuming that J20395358+4222505 is at the distance of Cyg OB2 we derive the radius from infrared photometry, finding R= 41.2 $\pm$ 4.0 R$_{\odot}$, log(L/L$_{\odot}$)= 5.71 $\pm$ 0.04 and a spectroscopic mass of 46.5 $\pm$ 15.0 M$_{\odot}$. The clumped mass-loss rate (clumping factor 10) is very high for the spectral type, $\dot{M}$ = 2.4x10$^{-6}$ M$_{\odot}$ a$^{-1}$. The high rotational velocity and mass-loss rate place the star at the hot side of the bi-stability jump. Together with the nearly solar CNO abundance pattern, they may also point to evolution in a binary system, J20395358+4222505 being the initial secondary.

preprint2021arXiv

Upper Mass-Loss Limits and Clumping in the Intermediate and Outer Wind Regions of OB stars

We probe the radial clumping stratification of OB stars in the intermediate and outer wind regions (r>~2 R*) to derive upper limits for mass-loss rates, and compare to current mass-loss implementation. Together with archival multi-wavelength data, our new far-infrared continuum observations for a sample of 25 OB stars (including 13 B Supergiants) uniquely constrain the clumping properties of the intermediate wind region. We derive the minimum radial stratification of the clumping factor through the stellar wind, fclmin(r), and the corresponding maximum mass-loss rate, Mdotmax, normalising clumping factors to the outermost wind region (clfar=1). The clumping degree for r>~2 R* decreases or stays constant with increasing radius for almost the whole sample. There is a dependence on luminosity class and spectral type at the intermediate region relative to the outer ones: O Supergiants (OSGs) present a factor 2 larger clumping factors than B Supergiants (BSGs). The maximum clumping of roughly 1/3 of the OB Supergiants occurs close to the wind base (r<~2 R*) and then decreases monotonically. This contrasts with the more frequent case where the lowermost clumping increases towards a maximum, and needs to be addressed by theoretical models. Additionally, the estimated Mdotmax for BSGs is at least one order of magnitude lower than theoretical values, whereas for OSGs our results and predictions agree within errors. Assuming values of clfar=4-9 from hydrodynamical models would imply a reduction of mass-loss rates included in stellar evolution models by a factor 2-3 for OSGs and by factors 6-200 for BSGs below the first bi-stability jump. This implies large reductions of mass-loss rates applied in evolution-models for BSGs, independently of the actual clumping properties of these winds, and a thorough re-investigation of BSG mass-loss rates and their effects on stellar evolution.

preprint2020arXiv

Properties of OB star-black hole systems derived from detailed binary evolution models

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

preprint2020arXiv

The VLT-FLAMES Tarantula Survey XXXII. Low-luminosity late O-type stars -- classification, main physical parameters, and silicon abundances

Analysis of late O-type stars observed in the Large Magellanic Cloud (LMC) by the VLT-FLAMES Tarantula Survey (VFTS) revealed a discrepancy between the physical properties estimated from model-atmosphere analysis and those expected from their morphological classifications. Here we revisit the analysis of 32 of these puzzling objects using new hydrogen-helium-silicon FASTWIND models and a different fitting approach to re-evaluate their physical properties. Our new analysis confirms that these stars indeed have properties that are typical of late O-type dwarfs. We also present the first estimates of silicon abundances for O-type stars in the 30 Dor clusters NGC 2060 and NGC 2070, with a weighted mean abundance for our sample of 7.05 +/- 0.03. Our values are about 0.20 dex lower than those previously derived for B-type stars in the LMC clusters N 11 and NGC 2004 using TLUSTY models. Various possibilities (e.g. differences in the analysis methods, effects of microturbulence, and real differences between stars in different clusters) were considered to account for these results. We also used our grid of FASTWIND models to reassess the impact of using the Galactic classification criteria for late O-type stars in the LMC by scrutinising their sensitivity to different stellar properties. At the cool edge of the O star regime the HeII 4686/HeI 4713 ratio used to assign luminosity class for Galactic stars can mimic giants or bright giants in the LMC, even for objects with high gravities (log_g > 4.0 dex). We argue that this line ratio is not a reliable luminosity diagnostic for late O-type stars in the LMC, and that the SiIV 4989/HeI4026 ratio is more robust for these types.

preprint2019arXiv

New predictions for radiation-driven, steady-state mass-loss and wind-momentum from hot, massive stars. I. Method and first results

[Abridged] Context: Radiation-driven mass loss plays a key role in the life-cycles of massive stars. However, basic predictions of such mass loss still suffer from significant quantitative uncertainties. Aims: We develop new radiation-driven, steady-state wind models for massive stars with hot surfaces, suitable for quantitative predictions of global parameters like mass-loss and wind-momentum rates. Methods: The simulations presented here are based on a self-consistent, iterative grid-solution to the spherically symmetric, steady-state equation of motion, using full NLTE radiative transfer solutions in the co-moving frame to derive the radiative acceleration. We do not rely on any distribution functions or parametrization for computation of the line force responsible for the wind driving. Results: In this first paper, we present models representing two prototypical O-stars in the Galaxy, one with a higher stellar mass M/Msun=59 and luminosity log10(L/Lsun)= 5.87 and one with a lower M/Msun= 27 and log10(L/Lsun)= 5.1. For these simulations, basic predictions for global mass-loss rates and velocity laws are given, and the influence from additional parameters like wind clumping and microturbulent speeds discussed. A key result is that our mass-loss rates are significantly lower than those predicted by the mass-loss recipes normally included in models of massive-star evolution. Conclusions: Our results support previous suggestions that Galactic O-star mass-loss rates may be overestimated in present-day stellar evolution models, and that new rates thus might be needed. Indeed, future papers in this series will incorporate our new models into such simulations of stellar evolution, extending the very first simulations presented here toward larger grids covering a range of metallicities, B supergiants across the bistability jump, and possibly also Wolf-Rayet stars.

preprint2018arXiv

Response to comment on "An excess of massive stars in the local 30 Doradus starburst"

Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradus.

preprint2016arXiv

Atmospheric NLTE-Models for the Spectroscopic Analysis of Blue Stars with Winds. III. X-ray emission from wind-embedded shocks

X-rays/EUV radiation emitted from wind-embedded shocks in hot, massive stars can affect the ionization balance in their outer atmospheres, and can be the mechanism responsible for the production of highly ionized species. To allow for these processes in the context of spectral analysis, we have implemented such emission into our unified, NLTE model atmosphere/spectrum synthesis code FASTWIND. The shock structure and corresponding emission is calculated as a function of user-supplied parameters. We account for a temperature and density stratification inside the post-shock cooling zones, calculated for radiative and adiabatic cooling in the inner and outer wind, respectively. The high-energy absorption of the cool wind is considered by adding important K-shell opacities, and corresponding Auger ionization rates have been included into the NLTE network. We tested and verified our implementation carefully against corresponding results from various alternative model atmosphere codes, and studied the effects from shock emission for important ions from He, C, N, O, Si, and P. Surprisingly, dielectronic recombination turned out to play an essential role for the ionization balance of OIV/OV around Teff = 45,000 K. Finally, we investigated the behavior of the mass absorption coefficient, kappa_nu(r), important in the context of X-ray line formation in massive star winds. In almost all considered cases, direct ionization is of major influence, and Auger ionization significantly affects only NVI and OVI. The approximation of a radially constant kappa_nu is justified for r > 1.2 Rstar and lambda < 18 A, and also for many models at longer wavelengths. To estimate the actual value of this quantity, however, the HeII opacities need to be calculated from detailed NLTE modeling, at least for wavelengths longer than 18 to 20 A, and information on the individual CNO abundances has to be present.

preprint2016arXiv

The IACOB project: III. New observational clues to understand macroturbulent broadening in massive O- and B-type stars

We aim to provide new empirical clues about macroturbulent spectral line broadening in O- and B-type stars to evaluate its physical origin. We use high-resolution spectra of ~430 stars with spectral types in the range O4-B9 (all luminosity classes). We characterize the line-broadening of adequate diagnostic metal lines using a combined FT and GOF technique. We perform a quantitative spectroscopic analysis of the whole sample using automatic tools coupled with a huge grid of FASTWIND models. We also incorporate quantitative information about line asymmetries to our observational description of the characteristics of the line-profiles, and present a comparison of the shape and type of line-profile variability found in a small sample of O stars and B supergiants with still undefined pulsational properties and B main sequence stars with variable line-profiles. We present a homogeneous and statistically significant overview of the (single snapshot) line-broadening properties of stars in the whole O and B star domain. We find empirical evidence of the existence of various types of non-rotational broadening agents acting in the realm of massive stars. Even though all of them could be quoted and quantified as a macroturbulent broadening from a practical point of view, their physical origin can be different. Contrarily to the early- to late-B dwarfs/giants, which present a mixture of cases in terms of line-profile shape and variability, the whole O-type and B supergiant domain (or, roughly speaking, stars with M_ZAMS > 15 M_sol) is fully dominated by stars with a remarkable non-rotational broadening component and very similar profiles (including type of variability). We provide some examples illustrating how this observational dataset can be used to evaluate scenarios aimed at explaining the existence of sources of non-rotational broadening in massive stars.

preprint2016arXiv

The R136 star cluster dissected with Hubble Space Telescope/STIS. I. Far-ultraviolet spectroscopic census and the origin of HeII 1640 in young star clusters

We introduce a HST/STIS stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS/MAMA spectroscopy of R136 using 17 contiguous 52x0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90% of the 57 sources brighter than m_F555W = 16.0 mag within a radius of 0.5 parsec of R136a1, plus 8 additional nearby sources including R136b (O4\,If/WN8). We measure wind velocities for 52 early-type stars from CIV 1548-51, including 16 O2-3 stars. For the first time we spectroscopically classify all Weigelt & Baier members of R136a, which comprise three WN5 stars (a1-a3), two O supergiants (a5-a6) and three early O dwarfs (a4, a7, a8). A complete Hertzsprung-Russell diagram for the most massive O stars in R136 is provided, from which we obtain a cluster age of 1.5+0.3_-0.7 Myr. In addition, we discuss the integrated ultraviolet spectrum of R136, and highlight the central role played by the most luminous stars in producing the prominent HeII 1640 emission line. This emission is totally dominated by very massive stars with initial masses above ~100 Msun. The presence of strong HeII 1640 emission in the integrated light of very young star clusters (e.g A1 in NGC 3125) favours an initial mass function extending well beyond a conventional upper limit of 100 Msun. We include montages of ultraviolet spectroscopy for LMC O stars in the Appendix. Future studies in this series will focus on optical STIS/CCD medium resolution observations.

preprint2015arXiv

The ephemeris, orbital decay, and masses of 10 eclipsing HMXBs

We take advantage of more than 10 years of monitoring of the eclipsing HMXB systems LMC X-4, Cen X-3, 4U 1700-377, 4U 1538-522, SMC X-1, IGR J18027-2016, Vela X-1, IGR J17252-3616, XTE J1855-026, and OAO 1657-415 with the ASM on-board RXTE and ISGRI on-board INTEGRAL to update their ephemeris. These results are used to refine previous measurements of the orbital period decay of all sources (where available) and provide the first accurate values of the apsidal advance in Vela X-1 and 4U 1538-522. Updated values for the masses of the neutron stars hosted in the ten HMXBs are also provided, as well as the long-term lightcurves folded on the sources best determined orbital parameters. These lightcurves reveal complex eclipse ingresses and egresses, that are understood mostly as being due to the presence of accretion wakes. The results reported in this paper constitute a database to be used for population and evolutionary studies of HMXBs, as well as theoretical modelling of long-term accretion in wind-fed X-ray binaries.

preprint2015arXiv

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

TLUSTY non-LTE model atmosphere calculations have been used to determine atmospheric parameters and nitrogen (N) abundances for 34 single and 18 binary B-type supergiants (BSGs). The effects of flux contribution from an unseen secondary were considered for the binary sample. We present the first systematic study of the incidence of binarity for a sample of BSGs across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the BSGs it may be necessary to extend the TAMS to lower temperatures. This is consistent with the derived distribution of mass discrepancies, projected rotational velocities (vsini) and N abundances, provided that stars cooler than this temperature are post RSG objects. For the BSGs in the Tarantula and previous FLAMES surveys, most have small vsini. About 10% have larger vsini (>100 km/s) but surprisingly these show little or no N enhancement. All the cooler BSGs have low vsini of <70km/s and high N abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. A lack of cool binaries, possibly reflects the small sample size. Single star evolutionary models, which include rotation, can account for the N enhancement in both the single and binary samples. The detailed distribution of N abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. The first comparative study of single and binary BSGs has revealed that the main sequence may be significantly wider than previously assumed, extending to Teff=20000K. Some marginal differences 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 evolution of BSGs. Full abstract in paper

preprint2014arXiv

Mass loss from inhomogeneous hot star winds III. An effective-opacity formalism for line radiative transfer in accelerating, clumped two-component media, and first results on theory and diagnostics

[Abridged] We develop and benchmark a fast and easy-to-use effective-opacity formalism for line and continuum radiative transfer in an accelerating two-component clumpy medium. The formalism bridges the limits of optically thin and thick clumps, and is here used to i) design a simple vorosity-modified Sobolev with exact integration (vmSEI) method for analyzing UV wind resonance lines in hot, massive stars, and ii) derive simple correction factors to the line force driving the outflows of such stars. We show that (for a given ionization factor) UV resonance doublets may be used to analytically predict the upward corrections in empirically inferred mass-loss rates associated with porosity in velocity space (a.k.a. velocity-porosity, or vorosity), but that severe solution degeneracies exist. For an inter-clump density set to 1 % of the mean density, we for O and B supergiants derive upward empirical mass-loss corrections of typically factors of either ~5 or ~50, depending on which of the two applicable solutions is chosen. Overall, our results indicate this solution dichotomy severely limits the use of UV resonance lines as direct mass-loss indicators of clumped hot stellar winds. We next apply the effective-opacity formalism to the standard CAK theory of line-driven winds. By analytic and numerical hydrodynamics calculations, we show that in cases where vorosity is important at the critical point setting the mass-loss rate, the reduced line-force leads to a lower theoretical mass loss, by a factor scaling with the normalized velocity filling factor fvel. On the other hand, if vorosity is important only above this critical point, the predicted mass loss is not affected, but the wind terminal speed is reduced. This shows that porosity in velocity space can have a significant impact not only on the diagnostics, but also on the dynamics and theory of radiatively driven winds.

preprint2014arXiv

Physics of mass loss in massive stars

We review potential mass-loss mechanisms in the various evolutionary stages of massive stars, from the well-known line-driven winds of O-stars and BA-supergiants to the less-understood winds from Red Supergiants. We discuss optically thick winds from Wolf-Rayet stars and Very Massive Stars, and the hypothesis of porosity-moderated, continuum-driven mass loss from stars formally exceeding the Eddington limit, which might explain the giant outbursts from Luminous Blue Variables. We finish this review with a glance on the impact of rapid rotation, magnetic fields and small-scale inhomogeneities in line-driven winds.

preprint2014arXiv

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, $v_e\sin i$, of $\sim$330 O-type objects, i.e. $\sim$210 spectroscopic single stars and $\sim$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 $v_e\sin i$ distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100\,$\rm{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\,$\rm{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.

preprint2014arXiv

The IACOB project: II. On the scatter of O-dwarf spectral type -- effective temperature calibrations

We are now in an era of large spectroscopic surveys of OB-type stars. Quantitative spectroscopic analysis of these modern datasets is enabling us to review the physical properties of blue massive stars with robust samples, not only revisiting mean properties and general trends, but also incorporating information about the effects of second-order parameters. We investigate the spectral type -- effective temperature (SpT - Teff) calibration for O-type dwarfs, and its claimed dependence on metallicity, using statistically-meaningful samples of stars extracted from the IACOB and VFTS surveys. We perform a homogeneous differential spectroscopic analysis of 33 Galactic and 53 LMC O~dwarfs (spanning spectral types of O4 -- O9.7) using the IACOB-GBAT package, a $χ^2$-fitting algorithm based on a large pre-computed grid of FASTWIND models, and standard techniques for the hydrogen/helium analysis of O-type stars. We compare the estimated effective temperatures and gravities as a function of (internally consistent) spectral classifications. While the general trend is that the temperature of a star increases with earlier spectral types and decreasing metallicity, we show that the large range of gravities found for O-type dwarfs -- spaning up to 0.45-0.50 dex in some spectral bins -- plays a critical role on the dependence of the effective temperature calibrations as a function of spectral type and metallicity. This result warns us about the use of SpT - Teff calibrations for O-dwarfs which ignore the effects of gravity, and highlights the risks of employing calibrations based on small samples. The effects of this scatter in gravities (evolutionary status) for O-type dwarfs should be included in future recipes which employ SpT - Teff calibrations.

preprint2013arXiv

On magnetic inhibition of photospheric macro-turbulence generated in the iron-bump opacity zone of O-stars

Massive, hot OB-stars show clear evidence of strong macroscopic broadening (in addition to rotation) in their photospheric spectral lines. This paper examines the occurrence of such "macro-turbulence" in slowly rotating O-stars with strong, organised surface magnetic fields. Focusing on the CIV 5811A line, we find evidence for significant macro-turbulent broadening in all stars except NGC1624-2, which also has (by far) the strongest magnetic field. Instead, the very sharp CIV lines in NGC1624-2 are dominated by magnetic Zeeman broadening, from which we estimate a dipolar field of approximately 20 kG. By contrast, magnetic broadening is negligible in the other stars (due to their weaker field strengths, on order 1 kG), and their CIV profiles are typically very broad and similar to corresponding lines observed in non-magnetic O-stars. Quantifying this by an isotropic, Gaussian macro-turbulence, we derive vmac = 2.2 (+- 0.9/2.2) km/s for NGC1624-2, and vmac = 20-65 km/s for the rest of the magnetic sample. We use these observational results to test the hypothesis that the field can stabilise the atmosphere and suppress the generation of macro-turbulence down to stellar layers where the magnetic pressure PB and the gas pressure Pg are comparable. Using a simple grey atmosphere to estimate the temperature T0 at which PB = Pg, we find that T0 > Teff for all investigated magnetic stars, but that T0 reaches the ~160000 K layers associated with the iron opacity-bump in hot stars only for NGC1624-2. This is consistent with the view that the responsible physical mechanism for photospheric O-star macro-turbulence may be stellar gravity-mode oscillations excited by sub-surface convection zones, and suggests that a sufficiently strong magnetic field can suppress such iron-bump generated convection and associated pulsational excitation.

preprint2013arXiv

Rotational properties of the O-type star population in the Tarantula region

The 30 Doradus (30\,Dor) region in the Large Magellanic Cloud (also known as the Tarantula Nebula) is the nearest massive starburst region, containing the richest sample of massive stars in the Local Group. It is the best possible laboratory to investigate aspects of the formation and evolution of massive stars. Here, we focus on rotation which is a key parameter in the evolution of these objects. We establish the projected rotational velocity, $v_{e}\sin i$, distribution of an unprecedented sample of 216 radial velocity constant ($\rm{ΔRV\, \leq\, 20 \,km s^{-1}}$) O-type stars in 30\,Dor observed in the framework of the VLT-FLAMES Tarantula Survey (VFTS). The distribution of $v_{e}\sin i$ shows a two-component structure: a peak around 80 $\rm{km s^{-1}}$ and a high-velocity tail extending up to $\sim$600 $\rm{km s^{-1}}$. Around 75% of the sample has 0 $\leq\, v_{e}\sin i \leq$ 200 $\rm{km s^{-1}}$ with the other 25% distributed in the high-velocity tail. The presence of the low-velocity peak is consistent with that found in other studies of late-O and early-B stars. The high-velocity tail is compatible with expectations from binary interaction synthesis models and may be predominantly populated by post-binary interaction, spun-up, objects and mergers. This may have important implications for the nature of progenitors of long-duration gamma ray bursts.

preprint2013arXiv

Spectroscopic and physical parameters of Galactic O-type stars. II. Observational constraints on projected rotational and extra broadening velocities as a function of fundamental parameters and stellar evolution

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 < 50 km/s. For the stars with Minit > 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.

preprint2013arXiv

The rotation rates of massive stars: How slow are the slow ones?

Context: Rotation plays a key role in the life cycles of stars with masses above 8 Msun. Hence, accurate knowledge of the rotation rates of such massive stars is critical for understanding their properties and for constraining models of their evolution. Aims: This paper investigates the reliability of current methods used to derive projected rotation speeds v sin i from line-broadening signatures in the photospheric spectra of massive stars, focusing on stars that are not rapidly rotating. Methods: We use slowly rotating magnetic O-stars with well-determined rotation periods to test the Fourier transform (FT) and goodness-of-fit (GOF) methods typically used to infer projected rotation rates of massive stars. Results: For our two magnetic test stars with measured rotation periods longer than one year, i.e., with v sin i < 1 km/s, we derive v sin i ~ 40-50 km/s from both the FT and GOF methods. These severe overestimates are most likely caused by an insufficient treatment of the competing broadening mechanisms referred to as microturbulence and macroturbulence. Conclusions: These findings warn us not to rely uncritically on results from current standard techniques to derive projected rotation speeds of massive stars in the presence of significant additional line broadening, at least when v sin i <~ 50 km/s. This may, for example, be crucial for i) determining the statistical distribution of observed rotation rates of massive stars, ii) interpreting the evolutionary status and spin-down histories of rotationally braked B-supergiants, and iii) explaining the deficiency of observed O-stars with spectroscopically inferred v sin i ~ 0 km/s. Further investigations of potential shortcomings of the above techniques are presently under way.

preprint2013arXiv

The VLT-FLAMES Tarantula Survey - XI. A census of the hot luminous stars and their feedback in 30 Doradus

We compile the first comprehensive census of hot luminous stars in the 30 Doradus (30 Dor) star forming region of the LMC. The census extends to a radius of 10arcmin (150pc) from the central cluster, R136. Stars were selected photometrically and combined with the latest spectral types. 1145 candidate hot luminous stars were identified of which >700 were considered genuine early type stars that contribute to feedback. We assess the spectroscopic completeness to be 85% in outer regions (>5pc) but fall to 35% in the vicinity of R136, giving a total of 500 hot luminous stars with spectroscopy. Stellar calibrations and models were used to obtain their physical parameters before integrated values were compared to global observations and the population synthesis code, Starburst99. The 31 W-R and Of/WN stars made large contributions to the total ionising and wind luminosities of ~40% and ~50%, respectively. Stars with Minit>100Msun also showed high contributions to the global feedback, ~25% in both cases. Such massive stars are not accounted for by the current Starburst99 code, which underestimated the ionising and wind luminosities of R136 by factors of ~2 and ~9, respectively. The census inferred a SFR of 0.073+/-0.04Msun/yr for 30 Dor, typically higher than results from popular SFR calibrations. However, it remained consistent with a far-UV luminosity tracer and a combined Halpha and mid-infrared tracer, but only after correcting for Halpha extinction. The global ionising output exceeded measurements from the associated gas and dust, suggesting ~6(+55/-6)% of ionising photons escape the region. When studying very luminous star forming regions, it is therefore essential to include the most massive stars to ensure a reliable energy budget. If 30 Dor is typical of other large star forming regions, estimates of the SFR will be underpredicted if this escape fraction is not accounted for.(abridged)

preprint2013arXiv

The VLT-FLAMES Tarantula Survey XII. Rotational velocities of the single O-type 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 Doradus (30 Dor). Methods. We measured projected rotational velocities, \vrot, by means of a Fourier transform method and a profile fitting method applied on a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, $\rm{P(\veq)}$, of the equatorial rotational velocity, \veq. Results. The distribution of \vrot\ shows a two-component structure: a peak around 80 \kms\ and a high-velocity tail extending up to $\sim$600 \kms. This structure is also present in the inferred distribution $\rm{P(\veq)}$ with around 80% of the sample having 0 $<$ \veq\, $\leq\, 300$ \kms\ and the other 20% distributed in the high-velocity region. 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. 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. Despite the fact that we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. The fact that the high-velocity tail may be preferentially (and perhaps even exclusively), populated by post-binary interaction products, has important implications for the evolutionary origin of systems that produce gamma-ray bursts.

preprint2013arXiv

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

AIMS. We test the hypothesis of O Vz stars (characterized by having HeII4686 stronger in absorption than other He lines in their blue-violet spectra) being at a younger evolutionary stage than are normal O-type dwarfs. METHODS. We have performed a quantitative spectroscopic analysis of a sample of 38 O Vz and 46 O V 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 both samples using the FASTWIND stellar atmosphere code and the IACOB-GBAT grid-based automatic tool. In the framework of a differential study, we compared the physical and evolutionary properties of both samples, regarding Teff, logg, logQ and logL. We also investigated the predictions of the FASTWIND code about 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 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.

preprint2012arXiv

A peculiar Of star in the Local Group galaxy IC 1613

Context. Results from the theory of radiatively driven winds are incorporated in stellar evolutionary and population synthesis models, and used in our interpretation of the observations of the deep Universe. Yet, the theory has been confirmed only until Small Magellanic Cloud metallicities. Analyses of O-stars at lower metallicities are needed to prove the theory. Aims. We have observed GHV-62024, an O6.5 IIIf star in the low-metallicity galaxy IC1613. According to a previous preliminary analysis this star could challenge the radiatively driven wind theory at low metallicities. Methods. Our observations were obtained with VIMOS at VLT, at R~2000 and were analysed using the latest version of the model atmosphere code FASTWIND, which includes N III Results. We obtain the stellar parameters and conclude that the star follows the average wind momentum-luminosity relationship (WLR) expected for its metallicity, but with a high value for the exponent of the wind velocity law, beta. We suggest that this high value may be reached because GHV-62024 could be a fast rotator seen at a low inclination angle. While the derived beta value does not change by adopting a lower wind terminal velocity, a wrong $V_\infty$ has a clear impact on the position of the star in the WLR diagram. The N and He abundances are very high, consistent with strong CNO mixing that could have been caused by the fast rotation, although we cannot discard a different origin. We find again the well-known mass-discrepancy. Conclusions. We conclude that the star follows the WLR expected for its metallicity. The results are consistent with GHV-62024 being a fast rotator seen close to pole-on, contaminated at the surface with CNO products and with a wind structure altered by the fast rotation without modifying the global WLR. We suggest that this could be a general property of fast rotators.

preprint2012arXiv

Nitrogen line spectroscopy in O-stars -- III. The earliest O-stars

This is the third paper in a series aiming at the analysis of nitrogen abundances in O-type stars, to enable further constraints on the early evolution of massive stars. We provide first theoretical predictions for the NIV4058/NIII4640 emission line ratio in dependence of various parameters, and confront them with results from the analysis of a sample of early LMC/SMC O-stars. Stellar and wind parameters are determined by line profile fitting of H/He/N lines, exploiting the helium and nitrogen ionization balance. Corresponding synthetic spectra are calculated using the NLTE atmospheric code FASTWIND. Though there is a monotonic relationship between the emission line ratio and Teff, all other parameters being equal, theoretical predictions indicate additional dependencies, most notably, on the nitrogen abundance. These basic predictions are confirmed by results from atmospheric code CMFGEN. The effective temperatures for the earliest O-stars, inferred from the nitrogen ionization balance, are partly considerably hotter than indicated by previous studies. Consistent with earlier results, effective temperatures increase from supergiants to dwarfs for all spectral types in the LMC. The relation between observed NIV4058/NIII4640 emission line ratio and Teff, for a given luminosity class, turned out to be quite monotonic for our sample stars, and fairly consistent with our model predictions. The scatter within a spectral sub-type is mainly produced by abundance effects. Our findings suggest that the Walborn et al. (2002) classification scheme is able to provide a meaningful relation between spectral type and Teff, provided that it is possible to discriminate for the luminosity class. This might be difficult to achieve in low-Z environments such as the SMC, owing to rather low wind-strengths. According to our predictions, the major bias of the classification scheme is due to nitrogen content.

preprint2011arXiv

L-band spectroscopy of Galactic OB-stars

Context. Mass-loss, occurring through radiation driven supersonic winds, is a key issue throughout the evolution of massive stars. Two outstanding problems are currently challenging the theory of radiation-driven winds: wind clumping and the weak-wind problem. Aims. We seek to obtain accurate mass-loss rates of OB stars at different evolutionary stages to constrain the impact of both problems in our current understanding of massive star winds. Methods. We perform a multi-wavelength quantitative analysis of a sample of ten Galactic OB-stars by means of the atmospheric code CMFGEN, with special emphasis on the L-band window. A detailed investigation is carried out on the potential of Brα and Pfγ as mass-loss and clumping diagnostics. Results. For objects with dense winds, Brα samples the intermediate wind while Pfγ maps the inner one. In combination with other indicators (UV, Hα, Brγ) these lines enable us to constrain the wind clumping structure and to obtain "true" mass-loss rates. For objects with weak winds, Brα emerges as a reliable diagnostic tool to constrain the mass-loss rates. The emission component at the line Doppler-core superimposed on the rather shallow Stark absorption wings reacts very sensitively to mass loss already at very low mass-loss values. On the other hand, the line wings display similar sensitivity to mass loss as Hα, the classical optical mass loss diagnostics. Conclusions. Our investigation reveals the great diagnostic potential of L-band spectroscopy to derive clumping properties and mass-loss rates of hot star winds. We are confident that Brα will become the primary diagnostic tool to measure very low mass-loss rates with unprecedented accuracy

preprint2011arXiv

Mass loss from inhomogeneous hot star winds II. Constraints from a combined optical/UV study

Mass-loss rates currently in use for hot, massive stars have recently been seriously questioned, mainly because of the effects of wind clumping. We investigate the impact of clumping on diagnostic ultraviolet resonance and optical recombination lines. Optically thick clumps, a non-void interclump medium, and a non-monotonic velocity field are all accounted for in a single model. We used 2D and 3D stochastic and radiation-hydrodynamic (RH) wind models, constructed by assembling 1D snapshots in radially independent slices. To compute synthetic spectra, we developed and used detailed radiative transfer codes for both recombination lines (solving the "formal integral") and resonance lines (using a Monte-Carlo approach). In addition, we propose an analytic method to model these lines in clumpy winds, which does not rely on optically thin clumping. Results: Synthetic spectra calculated directly from current RH wind models of the line-driven instability are unable to in parallel reproduce strategic optical and ultraviolet lines for the Galactic O-supergiant LCep. Using our stochastic wind models, we obtain consistent fits essentially by increasing the clumping in the inner wind. A mass-loss rate is derived that is approximately two times lower than predicted by the line-driven wind theory, but much higher than the corresponding rate derived from spectra when assuming optically thin clumps. Our analytic formulation for line formation is used to demonstrate the potential impact of optically thick clumping in weak-winded stars and to confirm recent results that resonance doublets may be used as tracers of wind structure and optically thick clumping. (Abridged)

preprint2011arXiv

Nitrogen line spectroscopy in O-stars -- II. Surface nitrogen abundances for O-stars in the Large Magellanic Cloud

This is the second paper in a series aiming at the analysis of nitrogen abundances in O-type stars, to enable further constraints on the early evolution of massive stars. We investigate the NIV lambda4058 emission line formation, provide nitrogen abundances for a substantial O-star sample in the LMC, and compare our findings with recent predictions from stellar evolutionary models. Stellar and wind parameters are determined by line profile fitting of hydrogen, helium and nitrogen lines, based on synthetic spectra calculated by FASTWIND. We derive nitrogen abundances for 20 O- and 5 B-stars, by analyzing all nitrogen lines present in the available optical spectra. The dominating process responsible for emission at NIV lambda4058 in O-stars is the strong depopulation of the lower level of the transition, which increases as a function of Mdot. Unlike the NIII triplet emission, resonance lines do not play a role for typical mass-loss rates and below. The bulk of our sample O-stars seems to be strongly nitrogen-enriched, and a clear correlation of nitrogen and helium enrichment is found. By comparing the nitrogen abundances as a function of vsini ('Hunter-plot') with tailored evolutionary calculations, we identify a considerable number of highly enriched objects at low rotation. Due to the low initial abundance, the detection of strong Nitrogen enrichment in the bulk of O-stars indicates that efficient mixing takes place already during the very early phases of stellar evolution of LMC O-stars.

preprint2011arXiv

Spectroscopic and physical parameters of Galactic O-type stars. I. Effects of rotation and spectral resolving power in the spectral classification of dwarfs and giants

The modern-era spectral classification of O-stars relies on either the Walborn or the Conti-Mathys scheme. Since both of these approaches have been developed using low-quality photographic data, their application to high-quality digital data might not be straightforward and be hampered by problems and complications that have not yet been appreciated. Using high-resolution spectra obtained with the ESO/MPG 2.2\,m telescope in La Silla and following the premises of the Walborn and Conti classification schemes, we determined the spectral types and luminosity classes of 19 Galactic O-type stars and compared them to those attributed by Walborn and Mathys based on low-quality data. Our analysis reveals that the morphological spectral types assigned using high-resolution data are systematically later (by up to 1.5 subtypes) then those attributed by Walborn. By means of line-profile simulations, we show that part of this discrepancy is more likely caused by the combined effect of stellar rotation and high spectral resolution on the depth of helium lines used as spectral type indicators. In addition, we demonstrate that at least for narrow-lined stars the "rotational effect" does not disappear when the high-resolution spectra are degraded to the resolution of the Walborn standards. We also find evidence of a systematic difference between our high-resolution quantitative spectral types and those assigned by Mathys. Rotation and spectral resolution are important third parameters in the spectral classification of O-type stars. To obtain reliable spectral classes within the Walborn approach, the unknown and the standard spectra must be compared at the same resolution and \vsini. Owing to resolution effects, the Conti approach might also need to be updated.

preprint2011arXiv

The IACOB project: A grid-based automatic tool for the quantitative spectroscopic analysis of O-stars

We present the IACOB grid-based automatic tool for the quantitative spectroscopic analysis of O-stars. The tool consists of an extensive grid of FASTWIND models, and a variety of programs implemented in IDL to handle the observations, perform the automatic analysis, and visualize the results. The tool provides a fast and objective way to determine the stellar parameters and the associated uncertainties of large samples of O-type stars within a reasonable computational time.

preprint2011arXiv

The VLT-FLAMES Tarantula Survey I: Introduction and observational overview

The VLT-FLAMES Tarantula Survey (VFTS) is an ESO Large Programme that has obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Here we introduce our scientific motivations and give an overview of the survey targets, including optical and near-infrared photometry and comprehensive details of the data reduction. One of the principal objectives was to detect massive binary systems via variations in their radial velocities, thus shaping the multi-epoch observing strategy. Spectral classifications are given for the massive emission-line stars observed by the survey, including the discovery of a new Wolf-Rayet star (VFTS 682, classified as WN5h), 2' to the northeast of R136. To illustrate the diversity of objects encompassed by the survey, we investigate the spectral properties of sixteen targets identified by Gruendl & Chu from Spitzer photometry as candidate young stellar objects or stars with notable mid-infrared excesses. Detailed spectral classification and quantitative analysis of the O- and B-type stars in the VFTS sample, paying particular attention to the effects of rotational mixing and binarity, will be presented in a series of future articles to address fundamental questions in both stellar and cluster evolution.

preprint2010arXiv

Are the Stellar Winds in IC1613 stronger than expected?

In this poster we present the results of our analyses of three early massive stars in IC 1613, whose spectra have been observed with VIMOS and analyzed with CMFGEN and FASTWIND. One of the targets resulted a possible LBV and the other two are Of stars with unexpectedly strong winds. The Of stars seem to be strongly contaminated by CNO products. Our preliminary results may represent a challenge for the theory of stellar atmospheres, but they still have to be confirmed by the analysis of more objects and a more complete coverage of the parameter space.

preprint2010arXiv

Is macroturbulence in OB Sgs related to pulsations?

As part of a long term observational project, we are investigating the macroturbulent broadening in O and B supergiants (Sgs) and its possible connection with spectroscopic variability phenomena and stellar oscillations. We present the first results of our project, namely firm observational evidence for a strong correlation between the extra broadening and photospheric line-profile variations in a sample of 13 Sgs with spectral types ranging from O9.5 to B8.

preprint2010arXiv

Is macroturbulent broadening in OB Supergiants related to pulsations?

The spectrum of O and B Supergiants is known to be affected by an important extra line-broadening (usually called macroturbulence) that adds to stellar rotation. Recent analysis of high resolution spectra has shown that the interpretation of this line-broadening as a consequence of large-scale turbulent motions would imply highly super-sonic velocity fields, making this scenario quite improbable. Stellar oscillations have been proposed as a likely alternative explanation. We present first encouraging results of an observational project aimed at investigating the $macroturbulent$ broadening in O and B Supergiants, and its possible connection with spectroscopic variability phenomena and stellar oscillations: a) all the studied B Supergiants show line profile variations, quantified by means of the first (< v >) and third velocity (< v^3 >) moments of the lines, b) there is a strong correlation between the peak-to-peak amplitudes of the < v > and < v^3 > variability and the size of the extra-broadening.

preprint2010arXiv

Mass loss from inhomogeneous hot star winds I. Resonance line formation in 2D models

Small-scale clumping in the winds of hot, massive stars is conventionally included in spectral analyses by assuming optically thin clumps, a void inter-clump medium, and a smooth velocity field. To reconcile investigations of different diagnostics within such models, a highly clumped wind with very low mass-loss rates needs to be invoked. Particularly, unsaturated UV resonance lines seem to indicate rates an order of magnitude (or even more) lower than previously accepted values. We investigate resonance line formation in inhomogeneous hot star winds with non-monotonic velocity fields by means of 2D stochastic and pseudo-2D radiation-hydrodynamic wind models. A Monte-Carlo radiative transfer code is presented and used to produce synthetic line spectra. Results: The optically thin clumping limit is only valid for very weak lines. For intermediate strong lines, the velocity spans of the clumps are of central importance. Current hydrodynamical models predict spans that are too large to reproduce observed profiles unless a very low mass-loss rate is invoked. By simulating lower spans in 2D stochastic models, the profile strengths become drastically reduced, and are consistent with higher mass-loss rates. To simultaneously meet the constraints from strong lines, the inter-clump medium must be non-void. A first comparison to the observed PV doublet in the O6 supergiant lam Cep confirms that a stochastic 2D model reproduces observations with a mass-loss rate roughly ten times higher than that derived from the same lines but assuming optically thin clumping. Tentatively this may resolve discrepancies between theoretical predictions, evolutionary constraints, and recent derived mass-loss rates, and suggests a re-investigation of the structure predicted by current hydrodynamical models.

preprint2010arXiv

Observational evidence for a correlation between macroturbulent broadening and line-profile variations in OB Supergiants

The spectra of O and B supergiants are known to be affected by a significant form of extra line broadening (usually referred to as macroturbulence) in addition to that produced by stellar rotation. Recent analyses of high resolution spectra have shown that the interpretation of this line broadening as a consequence of large scale turbulent motions would imply highly supersonic velocity fields in photospheric regions, making this scenario quite improbable. Stellar oscillations have been proposed as a likely alternative explanation. As part of a long term observational project, we are investigating the macroturbulent broadening in O and B supergiants and its possible connection with spectroscopic variability phenomena and stellar oscillations. In this letter, we present the first encouraging results of our project, namely firm observational evidence for a strong correlation between the extra broadening and photospheric line-profile variations in a sample of 13 supergiants with spectral types ranging from O9.5 to B8.

preprint2009arXiv

Collective pulsational velocity broadening due to gravity modes as a physical explanation for macroturbulence in hot massive stars

We aimed at finding a physical explanation for the occurrence of macroturbulence in the atmospheres of hot massive stars, a phenomenon found in observations since more than a decade but yet unexplained. We computed time series of line profiles for evolved massive stars broadened by rotation and by hundreds of low-amplitude nonradial gravity-mode pulsations which are predicted to be excited for evolved massive stars. In general, line profiles based on macrotubulent broadening can mimic those subject to pulsational broadening. In several cases, though, good fits require macroturbulent velocities that pass the speed of sound for realistic pulsation amplitudes. Moreover, we find that the rotation velocity can be seriously underestimated by using a simple parameter description for macroturbulence rather than an appropriate pulsational model description to fit the line profiles. We conclude that macroturbulence is a likely signature of the collective effect of pulsations. We provide line diagnostics and their typical values to decide whether or not pulsational broadening is present in observed line profiles, as well as a procedure to avoid an inaccurate estimation of the rotation velocity.

preprint2009arXiv

Spectroscopic determination of the fundamental parameters of 66 B-type stars in the field-of-view of the CoRoT satellite

We aim to determine the fundamental parameters of a sample of B stars with apparent visual magnitudes below 8 in the field-of-view of the CoRoT space mission, from high-resolution spectroscopy. We developed an automatic procedure for the spectroscopic analysis of B-type stars with winds, based on an extensive grid of FASTWIND model atmospheres. We use the equivalent widths and/or the line profile shapes of continuum normalized hydrogen, helium and silicon line profiles to determine the fundamental properties of these stars in an automated way. After thorough tests, both on synthetic datasets and on very high-quality, high-resolution spectra of B stars for which we already had accurate values of their physical properties from alternative analyses, we applied our method to 66 B-type stars contained in the ground-based archive of the CoRoT space mission. We discuss the statistical properties of the sample and compare them with those predicted by evolutionary models of B stars. Our spectroscopic results provide a valuable starting point for any future seismic modelling of the stars, should they be observed by CoRoT.

preprint2009arXiv

The VLT-FLAMES Tarantula Survey

The Tarantula Survey is an ambitious ESO Large Programme that has obtained multi-epoch spectroscopy of over 1,000 massive stars in the 30 Doradus region of the Large Magellanic Cloud. Here we introduce the scientific motivations of the survey and give an overview of the observational sample. Ultimately, quantitative analysis of every star, paying particular attention to the effects of rotational mixing and binarity, will be used to address fundamental questions in both stellar and cluster evolution.

preprint2005arXiv

A Medium Resolution Near-Infrared Spectral Atlas of O and Early B Stars

We present intermediate resolution (R ~ 8,000 - 12,000) high signal-to-noise H- and K-band spectroscopy of a sample of 37 optically visible stars, ranging in spectral type from O3 to B3 and representing most luminosity classes. Spectra of this quality can be used to constrain the temperature, luminosity and general wind properties of OB stars, when used in conjunction with sophisticated atmospheric model codes. Most important is the need for moderately high resolutions (R > 5000) and very high signal-to-noise (S/N > 150) spectra for a meaningful profile analysis. When using near-infrared spectra for a classification system, moderately high signal-to-noise (S/N ~ 100) is still required, though the resolution can be relaxed to just a thousand or two. In the appendix we provide a set of very high quality near-infrared spectra of Brackett lines in six early-A dwarfs. These can be used to aid in the modeling and removal of such lines when early-A dwarfs are used for telluric spectroscopic standards.