Showing 1–50 of 126 results for author: Puls, J

Empirical mass-loss rates and clumping properties of O-type stars in the LMC

Authors: C. Hawcroft, L. Mahy, H. Sana, J. O. Sundqvist, M. Abdul-Masih, S. A. Brands, L. Decin, A. deKoter, J. Puls

Abstract: We constrain wind parameters of a sample of 18 O-type stars in the LMC, through analysis with stellar atmosphere and wind models including the effects of optically thick clumping. This allows us to determine the most accurate spectroscopic mass-loss and wind structure properties of massive stars at sub-solar metallicity to date and gain insight into the impact of metallicity on massive stellar win… ▽ More We constrain wind parameters of a sample of 18 O-type stars in the LMC, through analysis with stellar atmosphere and wind models including the effects of optically thick clumping. This allows us to determine the most accurate spectroscopic mass-loss and wind structure properties of massive stars at sub-solar metallicity to date and gain insight into the impact of metallicity on massive stellar winds. Combining high signal to noise (S/N) ratio spectroscopy in the UV and optical gives us access to diagnostics of multiple different physical processes in the stellar wind. We produce synthetic spectra using the stellar atmosphere modelling code FASTWIND, and reproduce the observed spectra using a genetic algorithm based fitting technique. We empirically constrain 15 physical parameters associated with the stellar and wind properties, including temperature, surface gravity, surface abundances, rotation, macroturbulence and wind parameters. We find, on average, mass-loss rates a factor of 4-5 lower than those predicted by Vink et al. 2001, in good agreement with predictions from Bjorklund et al. 2021, and the best agreement with those from Krticka et al. 2018. In the 'weak-wind' regime we find mass-loss rates orders of magnitude below any theoretical predictions. We find a positive correlation of clumping factors (fcl) with effective temperature with an average fcl = 14 +- 8 for the full sample. Above 38 kK an average 46 +- 24% of the wind velocity span is covered by clumps and the interclump density is 10-30% of the mean wind. Below an effective temperature of roughly 38 kK there must be additional light leakage for supergiants. For dwarf stars at low temperatures there is a statistical preference for very low clump velocity spans, however it is unclear if this can be physically motivated as there are no clearly observable wind signatures in UV diagnostics. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 12 pages with 6 figures, plus 16 pages and 18 figures in appendix. Accepted for publication in A&A

X-Shooting ULLYSES: Massive Stars at low metallicity -- IV. Spectral analysis methods and exemplary results for O stars

Authors: A. A. C. Sander, J. -C. Bouret, M. Bernini-Peron, J. Puls, F. Backs, S. R. Berlanas, J. M. Bestenlehner, S. A. Brands, A. Herrero, F. Martins, O. Maryeva, D. Pauli, V. Ramachandran, P. A. Crowther, V. M. A. Gómez-González, A. C. Gormaz-Matamala, W. -R. Hamann, D. J. Hillier, R. Kuiper, C. J. K. Larkin, R. R. Lefever, A. Mehner, F. Najarro, L. M. Oskinova, E. C. Schösser , et al. (4 additional authors not shown)

Abstract: CONTEXT: The spectral analysis of hot, massive stars is a fundamental astrophysical method to obtain their intrinsic properties and their feedback. Quantitative spectroscopy for hot, massive stars requires detailed numerical modeling of the atmosphere and an iterative treatment to obtain the best solution within a given framework. AIMS: We present an overview of different techniques for the quanti… ▽ More CONTEXT: The spectral analysis of hot, massive stars is a fundamental astrophysical method to obtain their intrinsic properties and their feedback. Quantitative spectroscopy for hot, massive stars requires detailed numerical modeling of the atmosphere and an iterative treatment to obtain the best solution within a given framework. AIMS: We present an overview of different techniques for the quantitative spectroscopy of hot stars employed within the X-Shooting ULLYSES collaboration, from grid-based approaches to tailored fits. By performing a blind test, we gain an overview about the similarities and differences of the resulting parameters. Our study aims to provide an overview of the parameter spread caused by different approaches. METHODS: For three different stars from the sample (SMC O5 star AzV 377, LMC O7 star Sk -69 50, and LMC O9 star Sk -66 171), we employ different atmosphere codes (CMFGEN, Fastwind, PoWR) and strategies to determine their best-fitting model. For our analyses, UV and optical spectra are used to derive the properties with some methods relying purely on optical data for comparison. To determine the overall spectral energy distribution, we further employ additional photometry from the literature. RESULTS: Effective temperatures for each of three sample stars agree within 3 kK while the differences in log g can be up to 0.2 dex. Luminosity differences of up to 0.1 dex result from different reddening assumptions, which seem to be larger for the methods employing a genetic algorithm. All sample stars are nitrogen-enriched. CONCLUSIONS: We find a reasonable agreement between the different methods. Tailored fitting tends to be able to minimize discrepancies obtained with more course or automatized treatments. UV spectral data is essential for the determination of realistic wind parameters. For one target (Sk -69 50), we find clear indications of an evolved status. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 18+15 pages, 21+4 figures, under review at A&A, condensed abstract

X-Shooting ULLYSES: Massive stars at low metallicity -- V. Effect of metallicity on surface abundances of O stars

Authors: F. Martins, J. -C. Bouret, D. J. Hillier, S. A. Brands, P. A. Crowther, A. Herrero, F. Najarro, D. Pauli, J. Puls, V. Ramachandran, A. A. C. Sander, J. S. Vink, the XshootU collaboration

Abstract: Massive stars rotate faster, on average, than lower mass stars. Stellar rotation triggers hydrodynamical instabilities which transport angular momentum and chemical species from the core to the surface. Models of high-mass stars that include these processes predict that chemical mixing is stronger at lower metallicity. We aim to test this prediction by comparing the surface abundances of massive s… ▽ More Massive stars rotate faster, on average, than lower mass stars. Stellar rotation triggers hydrodynamical instabilities which transport angular momentum and chemical species from the core to the surface. Models of high-mass stars that include these processes predict that chemical mixing is stronger at lower metallicity. We aim to test this prediction by comparing the surface abundances of massive stars at different metallicities. We performed a spectroscopic analysis of single O stars in the Magellanic Clouds (MCs) based on the ULLYSES and XshootU surveys. We determined the fundamental parameters and helium, carbon, nitrogen, and oxygen surface abundances of 17 LMC and 17 SMC non-supergiant O6-9.5 stars. We complemented these determinations by literature results for additional MCs and also Galactic stars to increase the sample size and metallicity coverage. We investigated the differences in the surface chemical enrichment at different metallicities and compared them with predictions of three sets of evolutionary models. Surface abundances are consistent with CNO-cycle nucleosynthesis. The maximum surface nitrogen enrichment is stronger in MC stars than in Galactic stars. Nitrogen enrichment is also observed in stars with higher surface gravities in the SMC than in the Galaxy. This trend is predicted by models that incorporate chemical transport caused by stellar rotation. The distributions of projected rotational velocities in our samples are likely biased towards slow rotators. A metallicity dependence of surface abundances is demonstrated. The analysis of larger samples with an unbiased distribution of projected rotational velocities is required to better constrain the treatment of chemical mixing and angular momentum transport in massive single and binary stars. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 15 pages + appendix. Accepted in Astronomy & Astrophysics

Towards early-type eclipsing binaries as extragalactic milestones: III. Physical properties of the O-type eclipsing binary OGLE LMC-ECL-21568 in a quadruple system

Authors: Mónica Taormina, R. -P. Kudritzki, B. Pilecki, G. Pietrzyński, I. B. Thompson, J. Puls, M. Górski, B. Zgirski, D. Graczyk, W. Gieren, G. Hajdu

Abstract: We present the results from a complex study of an eclipsing O-type binary (Aa+Ab) with the orbital period $P_{A}=3.2254367$ days, that forms part of a higher-order multiple system in a configuration (A+B)+C. We derived masses of the Aa+Ab binary $M_{1}= 19.02 \pm 0.12 \,M_\odot$, $M_{2}= 17.50 \pm 0.13 \,M_\odot$, radii $R_{1}= 7.70 \pm 0.05 \,R_\odot$, $R_{2}= 6.64 \pm 0.06 \,R_\odot$, and temper… ▽ More We present the results from a complex study of an eclipsing O-type binary (Aa+Ab) with the orbital period $P_{A}=3.2254367$ days, that forms part of a higher-order multiple system in a configuration (A+B)+C. We derived masses of the Aa+Ab binary $M_{1}= 19.02 \pm 0.12 \,M_\odot$, $M_{2}= 17.50 \pm 0.13 \,M_\odot$, radii $R_{1}= 7.70 \pm 0.05 \,R_\odot$, $R_{2}= 6.64 \pm 0.06 \,R_\odot$, and temperatures $T_1 = 34250 \pm 500 $ K, $T_2 = 33750 \pm 500 $ K. From the analysis of radial velocities, we found a spectroscopic orbit of A in the outer A+B system with $P_{A+B}=195.8$ days ($P_{A+B}/P_{A}\approx 61$). In the O-C analysis, we confirmed this orbit and found another component orbiting the A+B system with $P_{AB+C}=2550$ days ($P_{AB+C}\,/P_{A+B}\approx 13$). From the total mass of the inner binary and its outer orbit, we estimated the mass of the third object, $M_B \gtrsim 10.7 M_\odot$. From the light-travel time effect fit to the O-C data, we obtained the limit for the mass of the fourth component, $M_C \gtrsim 7.3 M_\odot$. These extra components contribute to about 20% to 30% (increasing with wavelength) of the total system light. From the comparison of model spectra with the multiband photometry, we derived a distance modulus of 18.59 $\pm$ 0.06 mag, a reddening of 0.16 $\pm$ 0.02 mag, and an $R_V$ of $3.2$. This work is part of our ongoing project, which aims to calibrate the surface brightness-color relation for early-type stars. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 15 pages, 14 figures, 3 tables, accepted for publication in The Astrophysical Journal

The IACOB project X. Large-scale quantitative spectroscopic analysis of Galactic luminous blue stars

Authors: Abel de Burgos, Sergio Simón-Díaz, Miguel A. Urbaneja, Joachim Puls

Abstract: Blue supergiants (BSGs) are key objects for understanding the evolution of massive stars. However, discrepancies between theoretical predictions and empirical observations have opened up important questions yet to be answered. Studying statistically significant and unbiased samples of BSGs can help to improve the situation. We aim to perform a homogeneous and comprehensive quantitative spectroscop… ▽ More Blue supergiants (BSGs) are key objects for understanding the evolution of massive stars. However, discrepancies between theoretical predictions and empirical observations have opened up important questions yet to be answered. Studying statistically significant and unbiased samples of BSGs can help to improve the situation. We aim to perform a homogeneous and comprehensive quantitative spectroscopic analysis of a large sample of Galactic luminous blue stars (being the majority BSGs) from the IACOB spectroscopic database. We derive the projected rotational velocity ($v\sin i$) and macroturbulent broadening ($v_{\rm mac}$) using IACOB-BROAD. We used FASTWIND computations to derive effective temperatures ($T_{\rm eff}$), surface gravities, microturbulences ($ξ$), Si and He surface abundances, and the wind-strength parameter. We provide estimates of these quantities for the largest sample of Galactic luminous O9-B5 stars spectroscopically analyzed to date, with 527 targets. We find a drop in the relative number of stars at ~21 kK, coinciding with a scarcity of fast rotating stars below that temperature. We speculate that this feature might be roughly delineating the location of the Terminal-Age-Main-Sequence in the 15-85M$_\odot$ range. By investigating the main characteristics of the $v\sin i$ distribution as a function of $T_{\rm eff}$, we propose that an efficient mechanism transporting angular momentum from the stellar core to the surface might be operating along the main sequence. We find correlations between $ξ$, $v_{\rm mac}$ and the spectroscopic luminosity. We also find that no more than 20% of the stars in our sample have atmospheres clearly enriched in He, and suggest that the origin of this specific sub-sample might be in binary evolution. We do not find clear empirical evidence of an increase in the wind-strength over the wind bi-stability region towards lower temperatures. △ Less

Submitted 25 March, 2024; v1 submitted 30 November, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy and Astrophysics, 22 pages, 20 figures

To clump or not to clump The impact of wind inhomogeneities on the optical and NIR spectroscopic analysis of massive OB stars

Authors: K. Rübke, A. Herrero, J. Puls

Abstract: Winds of massive stars have density inhomogeneities (clumping) that may affect the formation of spectral lines in different ways, depending on their formation region. Most of previous and current spectroscopic analyses have been performed in the optical or ultraviolet domain. However, massive stars are often hidden behind dense clouds rendering near-infrared observations necessary. Our objective i… ▽ More Winds of massive stars have density inhomogeneities (clumping) that may affect the formation of spectral lines in different ways, depending on their formation region. Most of previous and current spectroscopic analyses have been performed in the optical or ultraviolet domain. However, massive stars are often hidden behind dense clouds rendering near-infrared observations necessary. Our objective is to investigate whether a spectroscopic analysis using either optical or infrared observations results in the same stellar parameters with comparable accuracy, and whether clumping affects them in different ways. We analyzed optical and near-infrared observations of a set of massive O stars with spectral types O4-O9.5 and all luminosity classes. We obtain similar stellar parameters in the optical and the infrared, although with larger uncertainties in the near-infrared, both with and without clumping, albeit with some individual deviating cases. We find that the inclusion of clumping improves the fit to H$_α$ or HeII 4686 in the optical for supergiants, as well as that of Br$_γ$ in the near-infrared, but it sometimes worsens the fit to HeII 2.18$μ$m. Globally, there are no significant differences when using the clumping laws tested in this work. The infrared can be used for spectroscopic analyses, giving similar parameters as from the optical, though with larger uncertainties. The best fits to different lines are obtained with different (linear) clumping laws, indicating that the wind structure may be more complex than adopted in the present work. No clumping law results in a better global fit, or improves the consistency between optical and infrared stellar parameters. Our work shows that the optical and infrared lines are not sufficient to break the dichotomy between the mass-loss rate and clumping factor. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Journal ref: A&A 679, A19 (2023)

X-Shooting ULLYSES: Massive stars at low metallicity. III. Terminal wind speeds of ULLYSES massive stars

Authors: C. Hawcroft, H. Sana, L. Mahy, J. O. Sundqvist, A. de Koter, P. A. Crowther, J. M. Bestenlehner, S. A. Brands, A. David-Uraz, L. Decin, C. Erba, M. Garcia, W. -R. Hamann, A. Herrero, R. Ignace, N. D. Kee, B. Kubátová, R. Lefever, A. Moffat, F. Najarro, L. Oskinova, D. Pauli, R. Prinja, J. Puls, A. A. C. Sander , et al. (4 additional authors not shown)

Abstract: The winds of massive stars have an impact on stellar evolution and on the surrounding medium. The maximum speed reached by these outflows, the terminal wind speed, is a global wind parameter and an essential input for models of stellar atmospheres and feedback. With the arrival of the ULLYSES programme, a legacy UV spectroscopic survey with HST, we have the opportunity to quantify the wind speeds… ▽ More The winds of massive stars have an impact on stellar evolution and on the surrounding medium. The maximum speed reached by these outflows, the terminal wind speed, is a global wind parameter and an essential input for models of stellar atmospheres and feedback. With the arrival of the ULLYSES programme, a legacy UV spectroscopic survey with HST, we have the opportunity to quantify the wind speeds of massive stars at sub-solar metallicity (in the Large and Small Magellanic Clouds, 0.5Z and 0.2Z) at an unprecedented scale. We empirically quantify the wind speeds of a large sample of OB stars, including supergiants, giants, and dwarfs at sub-solar metallicity. Using these measurements, we investigate trends of terminal wind speed with a number of fundamental stellar parameters, namely effective temperature, metallicity, and surface escape velocity. We empirically determined the terminal wind speed for a sample of 149 OB stars in the Magellanic Clouds either by directly measuring the maximum velocity shift of the absorption component of the Civ 1548-1550 line profile, or by fitting synthetic spectra produced using the Sobolev with exact integration method. Stellar parameters were either collected from the literature, obtained using spectral-type calibrations, or predicted from evolutionary models. We find strong trends of terminal wind speed with effective temperature and surface escape speed when the wind is strong enough to cause a saturated P Cygni profile in Civ 1548-1550. We find evidence for a metallicity dependence on the terminal wind speed proportional to Z^0.22+-0.03 when we compared our results to previous Galactic studies. Our results suggest that effective temperature rather than surface escape speed should be used as a straightforward empirical prediction of terminal wind speed and that the observed metallicity dependence is steeper than suggested by earlier works. △ Less

Submitted 24 March, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: 21 pages, 16 figures, 8 tables. Accepted in A&A

A model of anisotropic winds from rotating stars for evolutionary calculations

Authors: Ben Hastings, Norbert Langer, Joachim Puls

Abstract: Context: The surface properties of rotating stars can vary from pole to equator, resulting in anisotropic stellar winds which are not included in the currently available evolutionary models. Aims: We develop a formalism to describe the mass and angular momentum loss of rotating stars which takes into account both the varying surface properties and distortion due to rotation. Methods: Adopting… ▽ More Context: The surface properties of rotating stars can vary from pole to equator, resulting in anisotropic stellar winds which are not included in the currently available evolutionary models. Aims: We develop a formalism to describe the mass and angular momentum loss of rotating stars which takes into account both the varying surface properties and distortion due to rotation. Methods: Adopting the mass-loss recipe for non-rotating stars, we assigned to each point on the surface of a rotating star an equivalent non-rotating star, for which the surface mass flux is given by the recipe. The global mass-loss and angular momentum loss rates are then given by integrating over the deformed stellar surface as appropriate. Evolutionary models were computed and our prescription is compared to the currently used simple mass-loss enhancement recipes for rotating stars. Results: We find that mass-loss rates are largely insensitive to rotation for models not affected by the bi-stability jump. For those affected by the bi-stability jump, the increase in mass-loss rates with respect to time is smoothed. As our prescription considers the variation of physical conditions over the stellar surface, the region affected by the bi-stability jump is able to grow gradually instead of the whole star suddenly being affected. Conclusion: We have provided an easy to implement and flexible, yet physically meaningful prescription for calculating mass and angular momentum loss rates of rotating stars in a one-dimensional stellar evolution code which compares favourably to more physically comprehensive models. The implementation of our scheme in the stellar evolution code MESA is available online: https://zenodo.org/record/7437006 △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A. 12 pages

Journal ref: A&A 672, A60 (2023)

Spin-down and reduced mass loss in early-type stars with large-scale magnetic fields

Authors: Z. Keszthelyi, A. de Koter, Y. Götberg, G. Meynet, S. A. Brands, V. Petit, M. Carrington, A. David-Uraz, S. T. Geen, C. Georgy, R. Hirschi, J. Puls, K. J. Ramalatswa, M. E. Shultz, A. ud-Doula

Abstract: Magnetism can greatly impact the evolution of stars. In some stars with OBA spectral types there is direct evidence via the Zeeman effect for stable, large-scale magnetospheres, which lead to the spin-down of the stellar surface and reduced mass loss. So far, a comprehensive grid of stellar structure and evolution models accounting for these effects was lacking. For this reason, we computed and st… ▽ More Magnetism can greatly impact the evolution of stars. In some stars with OBA spectral types there is direct evidence via the Zeeman effect for stable, large-scale magnetospheres, which lead to the spin-down of the stellar surface and reduced mass loss. So far, a comprehensive grid of stellar structure and evolution models accounting for these effects was lacking. For this reason, we computed and studied models with two magnetic braking and two chemical mixing schemes in three metallicity environments with the MESA software instrument. We find notable differences between the subgrids, which affects the model predictions and thus the detailed characterisation of stars. We are able to quantify the impact of magnetic fields in terms of preventing quasi-chemically homogeneous evolution and producing slowly-rotating, nitrogen-enriched ("Group 2") stars. Our model grid is fully open access and open source. △ Less

Submitted 13 November, 2022; originally announced November 2022.

Comments: to appear in "Winds of Stars and Exoplanets" Proceedings IAUS 370, 2022, eds.: A.A. Vidotto, L. Fossati, J.S. Vink

The effects of surface fossil magnetic fields on massive star evolution: IV. Grids of models at Solar, LMC, and SMC metallicities

Authors: Z. Keszthelyi, A. de Koter, Y. Götberg, G. Meynet, S. A. Brands, V. Petit, M. Carrington, A. David-Uraz, S. T. Geen, C. Georgy, R. Hirschi, J. Puls, K. J. Ramalatswa, M. E. Shultz, A. ud-Doula

Abstract: Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the MESA software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to… ▽ More Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the MESA software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3-60 M$_\odot$), surface equatorial magnetic field strength (0-50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly-rotating, nitrogen-enriched ("Group 2") stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is always less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this grid of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of i) revisiting the derived stellar parameters of known magnetic stars, and ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes. △ Less

Submitted 28 October, 2022; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: Accepted for publication in MNRAS. A full reproduction package is available on Zenodo at https://doi.org/10.5281/zenodo.7069766

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

Authors: R. Björklund, J. O. Sundqvist, S. M. Singh, J. Puls, F. Najarro

Abstract: 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-… ▽ More 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. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: Accepted for publication in A&A

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

Authors: R. Ganss, J. L. Pledger, A. E. Sansom, P. A. James, J. Puls, S. M. Habergham-Mawson

Abstract: 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 T… ▽ More 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. △ Less

Submitted 28 March, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 19 pages, 10 figures, 10 tables, accepted for publication in MNRAS, minor changes at editor's request

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 512, Issue 1, May 2022, Pages 1541-1556

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

Authors: A. Herrero, S. R. Berlanas, A. Gil de Paz, F. Comerón, J. Puls, S. Ramírez Alegría, M. García, D. J. Lennon, F. Najarro, S. Simón-Díaz, M. A. Urbaneja, J. Gallego, E. Carrasco, J. Iglesias, R. Cedazo, M. L. García Vargas, A. Castillo-Morales, S. Pascual, N. Cardiel, A. Pérez-Calpena, P. Gómez-Alvarez, I. Martínez-Delgado

Abstract: 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 a… ▽ More 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. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS, 13 pages, 9 figures

The R136 star cluster dissected with Hubble Space Telescope/STIS. III. The most massive stars and their clumped winds

Authors: Sarah A. Brands, Alex de Koter, Joachim M. Bestenlehner, Paul A. Crowther, Jon O. Sundqvist, Joachim Puls, Saida M. Caballero-Nieves, Michael Abdul-Masih, Florian A. Driessen, Miriam García, Sam Geen, Götz Gräfener, Calum Hawcroft, Lex Kaper, Zsolt Keszthelyi, Norbert Langer, Hugues Sana, Fabian R. N. Schneider, Tomer Shenar, Jorick S. Vink

Abstract: Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind propert… ▽ More Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind properties of the massive stars in R136, in particular the parameters related to wind clumping. Methods: We simultaneously analyse optical and UV spectroscopy of 53 O-type and 3 WNh-stars using the FASTWIND model atmosphere code and a genetic algorithm. The models account for optically thick clumps and effects related to porosity and velocity-porosity, as well as a non-void interclump medium. Results: We obtain stellar parameters, surface abundances, mass-loss rates, terminal velocities and clumping characteristics and compare these to theoretical predictions and evolutionary models. The clumping properties include the density of the interclump medium and the velocity-porosity of the wind. For the first time, these characteristics are systematically measured for a wide range of effective temperatures and luminosities. Conclusions: We confirm a cluster age of 1.0-2.5 Myr and derive an initial stellar mass of $\geq 250 {\rm M}_\odot$ for the most massive star in our sample, R136a1. The winds of our sample stars are highly clumped, with an average clumping factor of $f_{\rm cl}=29\pm15$. We find tentative trends in the wind-structure parameters as a function of mass-loss rate, suggesting that the winds of stars with higher mass-loss rates are less clumped. We compare several theoretical predictions to the observed mass-loss rates and terminal velocities and find that none satisfactorily reproduces both quantities. The prescription of Krtička & Kubát (2018) matches best the observed mass-loss rates. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: Accepted for publication in A&A; Appendix I will not be included in the published version

Journal ref: A&A 663, A36 (2022)

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

Authors: M. M. Rubio-Díez, J. O. Sundqvist, F. Najarro, A. Traficante, J. Puls, L. Calzoletti, D. Figer

Abstract: 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 i… ▽ More 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. △ Less

Submitted 26 August, 2021; originally announced August 2021.

Comments: 23 pages + 12 appendices pages, 18 figures, accepted in Astronomy & Astrophysics

Journal ref: A&A 658, A61 (2022)

Empirical mass-loss rates and clumping properties of Galactic early-type O supergiants

Authors: C. Hawcroft, H. Sana, L. Mahy, J. O. Sundqvist, M. Abdul-Masih, J. C. Bouret, S. A. Brands, A. de Koter, F. A. Driessen, J. Puls

Abstract: We investigate the impact of optically thick clumping on stellar wind diagnostics in O supergiants and constrain wind parameters associated with porosity in velocity space. This is the first time the effects of optically thick clumping have been investigated for a sample of massive hot stars, using models including a full optically thick clumping description. We re-analyse spectroscopic observatio… ▽ More We investigate the impact of optically thick clumping on stellar wind diagnostics in O supergiants and constrain wind parameters associated with porosity in velocity space. This is the first time the effects of optically thick clumping have been investigated for a sample of massive hot stars, using models including a full optically thick clumping description. We re-analyse spectroscopic observations of a sample of eight O supergiants. Using a genetic algorithm wrapper around the NLTE atmosphere code FASTWIND we obtain simultaneous fits to optical and UV spectra and determine photospheric and wind properties and surface abundances. We provide empirical constraints on a number of wind parameters including the clumping factors, mass-loss rates and terminal wind velocities. Additionally, we establish the first systematic empirical constraints on velocity filling factors and interclump densities. These parameters describe clump distribution in velocity-space and density of the interclump medium in physical-space, respectively. We observe a mass-loss rate reduction of a factor of 3.6 compared to theoretical predictions from Vink et al. (2000), and mass-loss rates within a factor 1.4 of predictions from Björklund et al. (2021). We confirm that including optically thick clumping allows simultaneous fitting of recombination lines and resonance lines, including the unsaturated UV phosphorus lines (Pv 1118-1128), without reducing the phosphorus abundance. We find that, on average, half of the wind velocity field is covered by dense clumps. We also find that these clumps are 25 times denser than the average wind, and that the interclump medium is 3-10 times less dense than the mean wind. The former result agrees well with theoretical predictions, the latter suggests that lateral filling-in of radially compressed gas might be critical for setting the scale of the rarefied interclump matter. △ Less

Submitted 18 August, 2021; originally announced August 2021.

Comments: Accepted for publication in A&A; 35 pages, 20 figures

Journal ref: A&A 655, A67 (2021)

Constraining the Overcontact Phase in Massive Binary Evolution I. Mixing in V382 Cyg, VFTS 352, and OGLE SMC-SC10 108086

Authors: Michael Abdul-Masih, Hugues Sana, Calum Hawcroft, Leonardo A. Almeida, Sarah A. Brands, Selma E. de Mink, Stephen Justham, Norbert Langer, Laurent Mahy, Pablo Marchant, Athira Menon, Joachim Puls, Jon Sundqvist

Abstract: As potential progenitors of several exotic phenomena including gravitational wave sources, magnetic stars, and Be stars, close massive binary systems probe a crucial area of the parameter space in massive star evolution. Despite the importance of these systems, large uncertainties regarding the nature and efficiency of the internal mixing mechanisms still exist. In this work, we aim to provide rob… ▽ More As potential progenitors of several exotic phenomena including gravitational wave sources, magnetic stars, and Be stars, close massive binary systems probe a crucial area of the parameter space in massive star evolution. Despite the importance of these systems, large uncertainties regarding the nature and efficiency of the internal mixing mechanisms still exist. In this work, we aim to provide robust observational constraints on the internal mixing processes by spectroscopically analyzing a sample of three massive overcontact binaries at different metallicities. Using optical phase-resolved spectroscopic data, we perform an atmosphere analysis using more traditional 1D techniques and using state-of-the-art 3D techniques. We compare and contrast the assumptions and results of each technique and investigate how the assumptions affect the final derived atmospheric parameters. We find that in all three cases, both components of system are highly overluminous indicating either efficient internal mixing of helium or previous non-conservative mass transfer. However, we do not find strong evidence of helium or CNO surface abundance changes usually associated with mixing. Additionally, we find that in unequal mass systems, the measured effective temperature and luminosity of the less massive component places it very close to the more massive component on the Hertzsprung-Russell diagram. These results were obtained independently using both of the techniques mentioned above, which suggests that these measurements are robust. The observed discrepancies between the temperature and the surface abundance measurements when compared to theoretical expectations indicate that unaccounted for additional physical mechanisms may be at play. △ Less

Submitted 16 April, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: Accepted for publication in A&A; 28 pages, 20 figures

Mapping the core of the Tarantula Nebula with VLT-MUSE II. The spectroscopic Hertzsprung-Russell diagram of OB stars in NGC 2070

Authors: N. Castro, P. A. Crowther, C. J. Evans, J. S. Vink, J. Puls, A. Herrero, M. Garcia, F. J. Selman, M. M. Roth, S. Simón-Díaz

Abstract: We present the spectroscopic analysis of 333 OB-type stars extracted from VLT-MUSE observations of the central 30 x 30 pc of NGC 2070 in the Tarantula Nebula on the Large Magellanic Cloud, the majority of which are analysed for the the first time. The distribution of stars in the spectroscopic Hertzsprung-Russell diagram (sHRD) shows 281 stars in the main sequence. We find two groups in the main s… ▽ More We present the spectroscopic analysis of 333 OB-type stars extracted from VLT-MUSE observations of the central 30 x 30 pc of NGC 2070 in the Tarantula Nebula on the Large Magellanic Cloud, the majority of which are analysed for the the first time. The distribution of stars in the spectroscopic Hertzsprung-Russell diagram (sHRD) shows 281 stars in the main sequence. We find two groups in the main sequence, with estimated ages of 2.1$\pm$0.8 and 6.2$\pm$2 Myr. A subgroup of 52 stars is apparently beyond the main sequence phase, which we consider to be due to emission-type objects and/or significant nebular contamination affecting the analysis. As in previous studies, stellar masses derived from the sHRD are systematically larger than those obtained from the conventional HRD, with the differences being largest for the most massive stars. Additionally, we do not find any trend between the estimated projected rotational velocity and evolution in the sHRD. The projected rotational velocity distribution presents a tail of fast rotators that resembles findings in the wider population of 30 Doradus. We use published spectral types to calibrate the HeI$λ$4921/HeII$λ$5411 equivalent-width ratio as a classification diagnostic for early-type main sequence stars when the classical blue-visible region is not observed. Our model-atmosphere analyses demonstrate that the resulting calibration is well correlated with effective temperature. △ Less

Submitted 5 February, 2021; originally announced February 2021.

Comments: 19 pages, 12 figures. Accepted in A&A

Journal ref: A&A 648, A65 (2021)

Atmospheric NLTE models for the spectroscopic analysis of blue stars with winds. V. Complete comoving frame transfer, and updated modeling of X-ray emission

Authors: J. Puls, F. Najarro, J. O. Sundqvist, K. Sen

Abstract: Context. Obtaining precise stellar and wind properties and abundance patterns of massive stars is crucial to understanding their nature and interactions with their environments, as well as to constrain their evolutionary paths and end-products. Aims. To enable higher versatility and precision of the complete UV to optical range, we improve our NLTE atmosphere and spectrum synthesis code FASTWIND… ▽ More Context. Obtaining precise stellar and wind properties and abundance patterns of massive stars is crucial to understanding their nature and interactions with their environments, as well as to constrain their evolutionary paths and end-products. Aims. To enable higher versatility and precision of the complete UV to optical range, we improve our NLTE atmosphere and spectrum synthesis code FASTWIND. We also aim to obtain an advanced description of X-ray emission from wind-embedded shocks. Methods. We include a detailed comoving frame radiative transfer for the essential frequency range, but still apply methods that enable low turnaround times. Results. In most cases, our new results agree excellently with those from the alternative code CMFGEN, both regarding the total radiative acceleration, strategic optical lines, and the UV-range. The agreement regarding NIII4634-4640-4642 has improved, though there are still certain discrepancies, mostly related to line overlap effects in the extreme ultraviolet. In the UV range of our coolest models, we find differences in the predicted depression of the pseudo-continuum, most pronounced around Ly_alpha. The comparison between our new and previous FASTWIND version reveals an almost perfect agreement, except for NV 4603-4619. Using an improved, depth-dependent description for the filling factors of hot, X-ray emitting material, we confirm previous analytic scaling relations with our numerical models. Conclusions. We warn against uncritically relying on transitions, which are strongly affected by direct or indirect line-overlap effects. The predicted UV-continuum depression for the coolest grid-models needs to be checked, both observationally, and regarding the underlying atomic data. Wind lines from "super-ionized" ions such as OVI can, in principle, be used to constrain the distribution of wind-embedded shocks. △ Less

Submitted 4 November, 2020; originally announced November 2020.

Comments: This paper is dedicated to the late Adi Pauldrach, who left us much too early

Journal ref: Astronomy & Astrophysics (2020) 642, A172

The R136 star cluster dissected with Hubble Space Telescope/STIS. II. Physical properties of the most massive stars in R136

Authors: Joachim M. Bestenlehner, Paul A. Crowther, Saida M. Caballero-Nieves, Fabian R. N. Schneider, Sergio Simon-Diaz, Sarah A. Brands, Alex de Koter, Goetz Graefener, Artemio Herrero, Norbert Langer, Daniel J. Lennon, Jesus Maiz Apellaniz, Joachim Puls, Jorick S. Vink

Abstract: We present an optical analysis of 55 members of R136, the central cluster in the Tarantula Nebula of the Large Magellanic Cloud. Our sample was observed with STIS aboard the Hubble Space Telescope, is complete down to about 40\,$M_{\odot}$, and includes 7 very massive stars with masses over 100\,$M_{\odot}$. We performed a spectroscopic analysis to derive their physical properties. Using evolution… ▽ More We present an optical analysis of 55 members of R136, the central cluster in the Tarantula Nebula of the Large Magellanic Cloud. Our sample was observed with STIS aboard the Hubble Space Telescope, is complete down to about 40\,$M_{\odot}$, and includes 7 very massive stars with masses over 100\,$M_{\odot}$. We performed a spectroscopic analysis to derive their physical properties. Using evolutionary models we find that the initial mass function (IMF) of massive stars in R136 is suggestive of being top-heavy with a power-law exponent $γ\approx 2 \pm 0.3$, but steeper exponents cannot be excluded. The age of R136 lies between 1 and 2\,Myr with a median age of around 1.6\,Myr. Stars more luminous than $\log L/L_{\odot} = 6.3$ are helium enriched and their evolution is dominated by mass loss, but rotational mixing or some other form of mixing could be still required to explain the helium composition at the surface. Stars more massive than 40\,$M_{\odot}$ have larger spectroscopic than evolutionary masses. The slope of the wind-luminosity relation assuming unclumped stellar winds is $2.41\pm0.13$ which is steeper than usually obtained ($\sim 1.8$). The ionising ($\log Q_0\,[{\rm ph/s}] = 51.4$) and mechanical ($\log L_{\rm SW}\,[{\rm erg/s}] = 39.1$) output of R136 is dominated by the most massive stars ($>100\,M_{\odot}$). R136 contributes around a quarter of the ionising flux and around a fifth of the mechanical feedback to the overall budget of the Tarantula Nebula. For a census of massive stars of the Tarantula Nebula region we combined our results with the VLT-FLAMES Tarantula Survey plus other spectroscopic studies. We observe a lack of evolved Wolf-Rayet stars and luminous blue and red supergiants. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS, 21 pages, 14 figures, 2 tables plus 33 pages of supplementary material

New predictions for radiation-driven, steady-state mass-loss and wind-momentum from hot, massive stars II. A grid of O-type stars in the Galaxy and the Magellanic Clouds

Authors: R. Björklund, J. O. Sundqvist, J. Puls, F. Najarro

Abstract: Reliable predictions of mass-loss rates are important for massive-star evolution computations. We aim to provide predictions for mass-loss rates and wind-momentum rates of O-type stars, carefully studying the behaviour of these winds as functions of stellar parameters like luminosity and metallicity. We use newly developed steady-state models of radiation-driven winds to compute the global propert… ▽ More Reliable predictions of mass-loss rates are important for massive-star evolution computations. We aim to provide predictions for mass-loss rates and wind-momentum rates of O-type stars, carefully studying the behaviour of these winds as functions of stellar parameters like luminosity and metallicity. We use newly developed steady-state models of radiation-driven winds to compute the global properties of a grid of O-stars. The self-consistent models are calculated by means of an iterative solution to the equation of motion using full NLTE radiative transfer in the co-moving frame to compute the radiative acceleration. In order to study winds in different galactic environments, the grid covers main-sequence stars, giants and supergiants in the Galaxy and both Magellanic Clouds. We find a strong dependence of mass-loss on both luminosity and metallicity. Mean values across the grid are $\dot{M}\sim L_{\ast}^{2.2}$ and $\dot{M}\sim Z_{\ast}^{0.95}$, however we also find a somewhat stronger dependence on metallicity for lower luminosities. Similarly, the mass loss-luminosity relation is somewhat steeper for the SMC than for the Galaxy. In addition, the computed rates are systematically lower (by a factor 2 and more) than those commonly used in stellar-evolution calculations. Overall, our results agree well with observations in the Galaxy that account properly for wind-clumping, with empirical $\dot{M}$ vs. $Z_\ast$ scaling relations, and with observations of O-dwarfs in the SMC. Our results provide simple fit relations for mass-loss rates and wind momenta of massive O-stars stars as functions of luminosity and metallicity, valid in the range $T_{\rm eff} = 28000 - 45000$\,K. Due to the systematically lower $\dot{M}$, our new models suggest that new rates might be needed in evolution simulations of massive stars. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: Accepted for publication in A&A. 16 pages, 13 figures

Journal ref: A&A 648, A36 (2021)

Spectroscopic patch model for massive stars using PHOEBE II and FASTWIND

Authors: Michael Abdul-Masih, Hugues Sana, Kyle E. Conroy, Jon Sundqvist, Andrej Prša, Angela Kochoska, Joachim Puls

Abstract: Massive stars play an important role in the mechanical and chemical evolution of galaxies. Understanding the internal processes of these stars is vital to our understanding of their evolution and eventual end products. Deformations from spherical geometry are common for massive stars; however, the tools that are currently available for the study of these systems are almost exclusively one-dimensio… ▽ More Massive stars play an important role in the mechanical and chemical evolution of galaxies. Understanding the internal processes of these stars is vital to our understanding of their evolution and eventual end products. Deformations from spherical geometry are common for massive stars; however, the tools that are currently available for the study of these systems are almost exclusively one-dimensional. We present a new spectroscopic analysis tool tailored for massive stars that deviate from spherical symmetry. This code (entitled SPAMMS) is a spectroscopic patch model that takes the three-dimensional surface geometry of the system into account to produce spectral profiles at given phases and orientations. In using the Wilson-Devinney-like code PHOEBE in combination with the non local thermodynamic equilibrium (NLTE) radiative transfer code FASTWIND, we created a three-dimensional mesh that represents the surface geometry of our system and we assigned FASTWIND emergent intensity line profiles to each mesh triangle, which take the local parameters such as temperature, surface gravity, and radius into account. These line profiles were then integrated across the visible surface, where their flux contribution and radial velocity are taken into account, thus returning a final line profile for the visible surface of the system at a given phase. We demonstrate that SPAMMS can accurately reproduce the morphology of observed spectral line profiles for overcontact systems. Additionally, we show how line profiles of rapidly-rotating single stars differ when taking rotational distortion into account, and the effects that these can have on the determined parameters. Finally, we demonstrate the code's ability to reproduce the Rossiter-Mclaughlin and Struve-Sahade effects. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: 10 pages, 10 figures, accepted for publication in A&A

A search for strong magnetic fields in massive and very massive stars in the Magellanic Clouds

Authors: S. Bagnulo, G. A. Wade, Y. Naze, J. H. Grunhut, M. E. Shultz, D. J. Asher, P. A. Crowther, C. J. Evans, A. David-Uraz, I. D. Howarth, N. Morrell, M. S. Munoz, C. Neiner, J. Puls, M. K. Szymanski, J. S. Vink

Abstract: Despite their rarity, massive stars dominate the ecology of galaxies via their strong, radiatively-driven winds throughout their lives and as supernovae in their deaths. However, their evolution and subsequent impact on their environment can be significantly affected by the presence of a magnetic field. While recent studies indicate that about 7% of OB stars in the Milky Way host strong, stable, o… ▽ More Despite their rarity, massive stars dominate the ecology of galaxies via their strong, radiatively-driven winds throughout their lives and as supernovae in their deaths. However, their evolution and subsequent impact on their environment can be significantly affected by the presence of a magnetic field. While recent studies indicate that about 7% of OB stars in the Milky Way host strong, stable, organised (fossil) magnetic fields at their surfaces, little is known about the fields of very massive stars, nor the magnetic properties of stars outside our Galaxy. We aim to continue searching for strong magnetic fields in a diverse set of massive and very massive stars (VMS) in the Large and Small Magellanic Clouds (LMC/SMC), and we evaluate the overall capability of FORS2 to usefully search for and detect stellar magnetic fields in extra-galactic environments. We have obtained FORS2 spectropolarimetry of a sample of 41 stars, which principally consist of spectral types B, O, Of/WN, WNh, and classical WR stars in the LMC and SMC. Four of our targets are Of?p stars; one of them was just recently discovered. Each spectrum was analysed to infer the longitudinal magnetic field. No magnetic fields were formally detected in our study, although Bayesian statistical considerations suggest that the Of?p star SMC159-2 is magnetic with a dipolar field of the order of 2.4 to 4.4kG. In addition, our first constraints of magnetic fields in VMS provide interesting insights into the formation of the most massive stars in the Universe. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Comments: Accepted by A&A

Journal ref: A&A, 635, A163 (15p, 2020)

Towards early-type eclipsing binaries as extragalactic milestones: II. NLTE spectral analysis and stellar parameters of the detached O-type system OGLE-LMC-ECL-06782 in the LMC

Authors: Mónica Taormina, Rolf-Peter Kudritzki, Joachim Puls, Bogumił Pilecki, Eva Sextl, G. Pietrzyński, Miguel A. Urbaneja, Wolfgang Gieren

Abstract: We combine the NLTE spectral analysis of the detached O-type eclipsing binary OGLE-LMC-ECL-06782 with the analysis of the radial velocity curve and light curve to measure an independent distance to the LMC. In our spectral analysis we study composite spectra of the system at quadrature and use the information from radial velocity and light curve about stellar gravities, radii and component flux ra… ▽ More We combine the NLTE spectral analysis of the detached O-type eclipsing binary OGLE-LMC-ECL-06782 with the analysis of the radial velocity curve and light curve to measure an independent distance to the LMC. In our spectral analysis we study composite spectra of the system at quadrature and use the information from radial velocity and light curve about stellar gravities, radii and component flux ratio to derive effective temperature, reddening, extinction and intrinsic surface brightness. We obtain a distance modulus to the LMC of m - M = 18.53 +/- 0.04 mag. This value is 0.05 mag larger than the precision distance obtained recently from the analysis of a large sample of detached, long period late spectral type eclipsing binaries but agrees within the margin of the uncertainties. We also determine the surface brightnesses of the system components and find good agreement with the published surface brightness color relationship. A comparison of the observed stellar parameters with the prediction of stellar evolution based on the MESA stellar evolution code shows reasonable agreement, but requires a reduction of the internal angular momentum transport to match the observed rotational velocities. △ Less

Submitted 14 January, 2020; originally announced January 2020.

Comments: 18 pages, 16 figures, 3 tables, Accepted for publication in ApJ

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

Authors: N. Markova, J. Puls, P. L. Dufton, D. J. Lennon, C. J. Evans, A. de Koter, O. H. Ramirez-Agudelo, H. Sana, J. S. Vink

Abstract: 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… ▽ More 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. △ Less

Submitted 6 January, 2020; originally announced January 2020.

Comments: 18 pages, 10 figures accepted for publication in A&A

Journal ref: A&A 634, A16 (2020)

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

Authors: N. Langer, C. Schürmann, K. Stoll, P. Marchant, D. J. Lennon, L. Mahy, S. E. de Mink, M. Quast, W. Riedel, H. Sana, P. Schneider, A. Schootemeijer, Chen Wang, L. A. Almeida, J. M. Bestenlehner, J. Bodensteiner, N. Castro, S. Clark, P. A. Crowther, P. Dufton, C. J. Evans, L. Fossati, G. Gräfener, L. Grassitelli, N. Grin , et al. (16 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 10 April, 2020; v1 submitted 20 December, 2019; originally announced December 2019.

Comments: 19 pages, 14 figures, Astronomy and Astrophysics, in press

Journal ref: A&A 638, A39 (2020)

A 3D short-characteristics method for continuum and line scattering problems in the winds of hot stars

Authors: L. Hennicker, J. Puls, N. D. Kee, J. O. Sundqvist

Abstract: Context: Knowledge about hot, massive stars is usually inferred from quantitative spectroscopy. To analyse non-spherical phenomena, the existing 1D codes must be extended to higher dimensions, and corresponding tools need to be developed. Aims: We present a 3D radiative transfer code that is capable of calculating continuum and line scattering problems in the winds of hot stars. By considering sph… ▽ More Context: Knowledge about hot, massive stars is usually inferred from quantitative spectroscopy. To analyse non-spherical phenomena, the existing 1D codes must be extended to higher dimensions, and corresponding tools need to be developed. Aims: We present a 3D radiative transfer code that is capable of calculating continuum and line scattering problems in the winds of hot stars. By considering spherically symmetric test models, we discuss potential error sources, and indicate advantages and disadvantages by comparing different solution methods. Further, we analyse the UV resonance line formation in the winds of rapidly rotating O stars. Methods: We consider both a (simplified) continuum model including scattering and thermal sources, and a UV resonance line transition approximated by a two-level-atom. We applied the short-characteristics (SC) method, using linear or monotonic Bézier interpolations, to solve the equation of radiative transfer on a non-uniform Cartesian grid. To calculate scattering dominated problems, our solution method is supplemented by an accelerated $Λ$-iteration scheme. Results: For the spherical test models, the mean relative error of the source function is on the $5-20\,\%$ level, depending on the applied interpolation technique and the complexity of the considered model. All calculated line profiles are in excellent agreement with corresponding 1D solutions. The predicted line profiles from fast rotating stars show a distinct behaviour as a function of rotational speed and inclination. This behaviour is tightly coupled to the wind structure and the description of gravity darkening and stellar surface distortion. Conclusions: Our SC methods are ready to be used for quantitative analyses of UV resonance line profiles. When calculating optically thick continua, both SC methods give reliable results, in contrast to the alternative finite-volume method. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Journal ref: A&A 633, A16 (2020)

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

Authors: J. O. Sundqvist, R. Björklund, J. Puls, F. Najarro

Abstract: [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 rat… ▽ More [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. △ Less

Submitted 15 October, 2019; originally announced October 2019.

Comments: 13 pages, 10 figures, accepted for publication in Astronomy and Astrophysics

Modelling the photometric variability of magnetic massive stars with the Analytical Dynamical Magnetosphere model

Authors: M. S. Munoz, G. A. Wade, Y. Nazé, J. Puls, S. Bagnulo, M. K. Szymański

Abstract: In this paper, we investigate the photometric variability of magnetic O-type stars. Such stars possess oblique, predominantly dipolar magnetic fields that confine their winds roughly axisymmetrically about the magnetic equator, thus forming a magnetosphere. We interpret their photometric variability as phase-dependent magnetospheric occultations. For massive star winds dominated by electron scatte… ▽ More In this paper, we investigate the photometric variability of magnetic O-type stars. Such stars possess oblique, predominantly dipolar magnetic fields that confine their winds roughly axisymmetrically about the magnetic equator, thus forming a magnetosphere. We interpret their photometric variability as phase-dependent magnetospheric occultations. For massive star winds dominated by electron scattering opacity in the optical and NIR, we can compute synthetic light curves from simply knowing the magnetosphere's mass density distribution. We exploit the newly-developed Analytical Dynamical Magnetosphere model (ADM) in order to obtain the predicted circumstellar density structures of magnetic O-type stars. The simplicity in our light curve synthesis model allows us to readily conduct a parameter space study. For validation purposes, we first apply our algorithm to HD 191612, the prototypical Of?p star. Next, we attempt to model the photometric variability of the Of?p-type stars identified in the Magellanic Clouds using OGLE photometry. We evaluate the compatibility of the ADM predictions with the observed photometric variations, and discuss the magnetic field properties that are implied by our modelling. △ Less

Submitted 13 October, 2019; originally announced October 2019.

Comments: 17 pages, 11 figures

Journal ref: MNRAS, 492, 1199-1213 (2020)

Clues on the Origin and Evolution of Massive Contact Binaries: Atmosphere Analysis of VFTS 352

Authors: Michael Abdul-Masih, Hugues Sana, Jon Sundqvist, Laurent Mahy, Athira Menon, Leonardo A. Almeida, Alex De Koter, Selma E. de Mink, Stephen Justham, Norbert Langer, Joachim Puls, Tomer Shenar, Frank Tramper

Abstract: The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved… ▽ More The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: 24 pages, 4 tables, 12 figures, accepted for publication in ApJ

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

Authors: F. R. N. Schneider, H. Sana, C. J. Evans, J. M. Bestenlehner, N. Castro, L. Fossati, G. Gräfener, N. Langer, O. H. Ramírez-Agudelo, C. Sabín-Sanjulián, S. Simón-Díaz, F. Tramper, P. A. Crowther, A. de Koter, S. E. de Mink, P. L. Dufton, M. Garcia, M. Gieles, V. Hénault-Brunet, A. Herrero, R. G. Izzard, V. Kalari, D. J. Lennon, J. Maíz Apellániz, N. Markova , et al. (7 additional authors not shown)

Abstract: 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 Doradu… ▽ More 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. △ Less

Submitted 25 July, 2018; originally announced July 2018.

Comments: Authors' version of a response to a technical comment published in Science; 8 pages, 1 figure

3D radiative transfer: Continuum and line scattering in non-spherical winds from OB stars

Authors: L. Hennicker, J. Puls, N. D. Kee, J. O. Sundqvist

Abstract: Context: State of the art quantitative spectroscopy of OB-stars compares synthetic spectra (calculated by means of 1D, spherically symmetric computer codes) with observations. Certain stellar atmospheres, however, show strong deviations from spherical symmetry, and need to be treated in 3D. Aims: We present a newly developed 3D radiative transfer code, tailored to the solution of the radiation fie… ▽ More Context: State of the art quantitative spectroscopy of OB-stars compares synthetic spectra (calculated by means of 1D, spherically symmetric computer codes) with observations. Certain stellar atmospheres, however, show strong deviations from spherical symmetry, and need to be treated in 3D. Aims: We present a newly developed 3D radiative transfer code, tailored to the solution of the radiation field in rapidly expanding stellar atmospheres. We apply our code to the continuum transfer in wind-ablation models, and to the UV resonance line formation in magnetic winds. Methods: We have used a 3D finite-volume method for the solution of the equation of radiative transfer, to study continuum- and line-scattering problems. Convergence has been accelerated by a non-local approximate Lambda-iteration scheme. Particular emphasis has been put on careful (spherically symmetric) test cases. Results: Typical errors of the source functions, when compared to 1D solutions, are of the order of 10-20 %, and increase for optically thick continua. In circumstellar discs, the radiation temperatures in the (optically thin) transition region from wind to disc are quite similar to corresponding values in the wind. For MHD simulations of dynamical magnetospheres, the line profiles, calculated with our 3D code, agree well with previous solutions using a 3D-SEI method. When compared with profiles resulting from the `analytic dynamical magnetosphere' (ADM) model, significant differences become apparent. Conclusions: Due to similar radiation temperatures in the wind and the transition region to the disc, the same line-strength distribution can be applied within radiation hydrodynamic calculations for circumstellar discs in `accreting high-mass stars'. To properly describe the UV line formation in dynamical magnetospheres, the ADM model needs to be further developed, at least in a large part of the outer wind. △ Less

Submitted 21 June, 2018; originally announced June 2018.

Journal ref: A&A 616, A140 (2018)

Atmospheric NLTE-models for the spectroscopic analysis of blue stars with winds. IV. Porosity in physical and velocity space

Authors: J. O. Sundqvist, J. Puls

Abstract: [Abridged] Clumping in the radiation-driven winds of hot, massive stars affects the derivation of synthetic observables across the electromagnetic spectrum. We implement a formalism for treating wind clumping - in particular the light-leakage effects associated with a medium that is porous in physical and velocity space - into the global (photosphere+wind) NLTE model atmosphere code FASTWIND. We a… ▽ More [Abridged] Clumping in the radiation-driven winds of hot, massive stars affects the derivation of synthetic observables across the electromagnetic spectrum. We implement a formalism for treating wind clumping - in particular the light-leakage effects associated with a medium that is porous in physical and velocity space - into the global (photosphere+wind) NLTE model atmosphere code FASTWIND. We assume a stochastic, two-component wind consisting of a mixture of optically thick and thin clumps embedded in a rarefied inter-clump medium. We account fully for the reductions in opacity associated with porosity in physical and velocity-space, and for the well-known effect that opacities depending on rho^2 are higher in clumpy winds than in smooth ones of equal mass-loss rate. By formulating our method in terms of suitable mean and effective opacities for the clumpy wind, we are able to compute models with the same speed (~15 min. on a modern laptop) as in previous code-generations. Some first, generic results of the new models include: i) Confirming earlier results that velocity-space porosity is critical for analysis of UV wind lines in O-stars; ii) for the optical Halpha line, optically thick clumping effects are small for O-stars, but potentially very important for late B and A-supergiants; iii) spatial porosity is a marginal effect for absorption of high-energy X-rays in O-stars, as long as the mean-free path between clumps are kept at realistic values; iv) porosity is negligible at typical O-star radio-photosphere radii; v) regarding the wind ionization balance, a general trend is that increased rates of recombination in simulations with optically thin clumps lead to overall lower degrees of ionization than in corresponding smooth models, but that this effect now is counteracted by the increased levels of light-leakage associated with porosity in physical and velocity space. △ Less

Submitted 28 May, 2018; originally announced May 2018.

Comments: 12 pages, 5 figures, accepted for publication in Astronomy & Astrophysics

Spectroscopic and physical parameters of Galactic O-type stars. III. Mass discrepancy and rotational mixing

Authors: Nevena Markova, Joachim Puls, Norbert Langer

Abstract: Massive stars play a key role in the evolution of the Universe. Our goal is to compare observed and predicted properties of single Galactic O stars to identify and constrain uncertain physical parameters and processes in stellar evolution and atmosphere models. We used a sample of 53 objects with spectral types from O3 to O9.7. For 30 of these, we determined the main photospheric and wind paramete… ▽ More Massive stars play a key role in the evolution of the Universe. Our goal is to compare observed and predicted properties of single Galactic O stars to identify and constrain uncertain physical parameters and processes in stellar evolution and atmosphere models. We used a sample of 53 objects with spectral types from O3 to O9.7. For 30 of these, we determined the main photospheric and wind parameters, using optical spectroscopy and applying the FASTWIND code. For the remaining objects, literature data, obtained by means of the CMFGEN code, were used instead. The properties of our sample were compared to published predictions based on two grids evolution models that include rotationally induced mixing. Within each luminosity class, we find a close correlation of N surface abundance and luminosity, and a stronger N enrichment in more massive and evolved O stars. Additionally, a correlation of the surface nitrogen and helium abundances is observed. The large number of nitrogen-enriched stars above ~30 solar masses argues for rotationally induced mixing as the most likely explanation. However, none of the considered models can match the observed trends correctly, especially in the high mass regime. We confirm mass discrepancy for objects in the low mass O-star regime. We conclude that the rotationally induced mixing of helium to the stellar surface is too strong in some of the models. We also suggest that present inadequacies of the models to represent the N enrichment in more massive stars with relatively slow rotation might be related to problematic efficiencies of rotational mixing. We are left with a picture in which invoking binarity and magnetic fields is required to achieve a more complete agreement of the observed surface properties of a population of massive main- sequence stars with corresponding evolutionary models. △ Less

Submitted 9 March, 2018; originally announced March 2018.

Comments: 30 pages, 16 figures, 3 appendices

Journal ref: A&A 613, A12 (2018)

The Of?p stars of the Magellanic Clouds: Are they strongly magnetic?

Authors: Melissa Munoz, Gregg Wade, Yaël Nazé, Stefano Bagnulo, Joachim Puls

Abstract: All known Galactic Of?p stars have been shown to host strong, organized, magnetic fields. Recently, five Of?p stars have been discovered in the Magellanic Clouds. They posses photometric \citep{Naze} and spectroscopic \citep{Walborn} variability compatible with the Oblique Rotator Model (ORM). However, their magnetic fields have yet to be directly detected. We have developed an algorithm allowing… ▽ More All known Galactic Of?p stars have been shown to host strong, organized, magnetic fields. Recently, five Of?p stars have been discovered in the Magellanic Clouds. They posses photometric \citep{Naze} and spectroscopic \citep{Walborn} variability compatible with the Oblique Rotator Model (ORM). However, their magnetic fields have yet to be directly detected. We have developed an algorithm allowing for the synthesis of photometric observables based on the Analytic Dynamical Magnetosphere (ADM) model of \citet{Owocki}. We apply our model to OGLE photometry in order to constrain their magnetic geometries and surface dipole strengths. We predict that the field strengths for some of these candidate extra-Galactic magnetic stars may be within the detection limits of the FORS2 instrument. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: 5 pages, 2 figures

An excess of massive stars in the local 30 Doradus starburst

Authors: F. R. N. Schneider, H. Sana, C. J. Evans, J. M. Bestenlehner, N. Castro, L. Fossati, G. Gräfener, N. Langer, O. H. Ramírez-Agudelo, C. Sabín-Sanjulián, S. Simón-Díaz, F. Tramper, P. A. Crowther, A. de Koter, S. E. de Mink, P. L. Dufton, M. Garcia, M. Gieles, V. Hénault-Brunet, A. Herrero, R. G. Izzard, V. Kalari, D. J. Lennon, J. Maíz Apellániz, N. Markova , et al. (7 additional authors not shown)

Abstract: The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation s… ▽ More The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation started about $8\,\mathrm{Myr}$ ago and the star-formation rate seems to have declined in the last $1\,\mathrm{Myr}$. The IMF is densely sampled up to $200\,\mathrm{M}_\odot$ and contains $32\pm12\%$ more stars above $30\,\mathrm{M}_\odot$ than predicted by a standard Salpeter IMF. In the mass range $15-200\,\mathrm{M}_\odot$, the IMF power-law exponent is $1.90^{+0.37}_{-0.26}$, shallower than the Salpeter value of 2.35. △ Less

Submitted 9 January, 2018; originally announced January 2018.

Comments: Authors' version of paper published in Science at http://science.sciencemag.org/content/359/6371/69 ; 15 pages main text (3 figures), 47 pages supplementary materials (10 figures, 3 tables)

Journal ref: Science 359, 69 (2018)

The IACOB project. V. Spectroscopic parameters of the O-type stars in the modern grid of standards for spectral classification

Authors: G. Holgado, S. Simón-Díaz, R. H. Barbá, J. Puls, A. Herrero, N. Castro, M. Garcia, J. Maíz Apellániz, I. Negueruela, C. Sabín-Sanjulián

Abstract: The IACOB and OWN surveys are two ambitious complementary observational projects which have made available a large multi-epoch spectroscopic database of optical high resolution spectra of Galactic massive O-type stars. As a first step in the study of the full sample of (more than 350) O stars surveyed by the IACOB/OWN projects, we have performed the quantitative spectroscopic analysis of a subsamp… ▽ More The IACOB and OWN surveys are two ambitious complementary observational projects which have made available a large multi-epoch spectroscopic database of optical high resolution spectra of Galactic massive O-type stars. As a first step in the study of the full sample of (more than 350) O stars surveyed by the IACOB/OWN projects, we have performed the quantitative spectroscopic analysis of a subsample of 128 stars included in the modern grid of O-type standards for spectral classification. We use semi-automatized tools to determine the set of spectroscopic parameters that can be obtained from the optical spectrum of O-type stars. We also benefit from the multi-epoch character of the surveys to perform a spectroscopic variability study of the sample, accounting for spectroscopic binarity and variability of the main wind diagnostic lines. We provide a general overview of the stellar and wind parameters of this reference sample, and updated recipes for the SpT\,--\,Teff/log g calibrations for Galactic O-type stars. We evaluate our semi-automatized analysis strategy with $\sim$40 stars from the literature, and find a good agreement. The agreement between the synthetic spectra associated with fastwind best fitting models and the observed spectra is good for most targets, but 46 stars present a particular behavior of the wind diagnostic lines that cannot be reproduced by our grid of spherically symmetric unclumped models. These are potential targets of interest for more detailed investigations of clumpy winds and/or the existence of additional circumstellar components. Last, our variability study has led to the detection of signatures of spectroscopic binarity in 27\% of the stars and small amplitude radial velocity variations in the photospheric lines of another 30\%. Additionally, 31\% of the investigated stars show variability in the wind diagnostic lines. △ Less

Submitted 6 June, 2018; v1 submitted 27 November, 2017; originally announced November 2017.

Comments: 20 pages, 18 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 613, A65 (2018)

2D wind clumping in hot, massive stars from hydrodynamical line-driven instability simulations using a pseudo-planar approach

Authors: J. O. Sundqvist, S. P. Owocki, J. Puls

Abstract: Context: Clumping in the radiation-driven winds of hot, massive stars arises naturally due to the strong, intrinsic instability of line-driving (the `LDI'). But LDI wind models have so far mostly been limited to 1D, mainly because of severe computational challenges regarding calculation of the multi-dimensional radiation force. Aims: To simulate and examine the dynamics and multi-dimensional natur… ▽ More Context: Clumping in the radiation-driven winds of hot, massive stars arises naturally due to the strong, intrinsic instability of line-driving (the `LDI'). But LDI wind models have so far mostly been limited to 1D, mainly because of severe computational challenges regarding calculation of the multi-dimensional radiation force. Aims: To simulate and examine the dynamics and multi-dimensional nature of wind structure resulting from the LDI. Methods: We introduce a `pseudo-planar', `box-in-a-wind' method that allows us to efficiently compute the line-force in the radial and lateral directions, and then use this approach to carry out 2D radiation-hydrodynamical simulations of the time-dependent wind. Results: Our 2D simulations show that the LDI first manifests itself by mimicking the typical shell-structure seen in 1D models, but how these shells then quickly break up into complex 2D density and velocity structures, characterized by small-scale density `clumps' embedded in larger regions of fast and rarefied gas. Key results of the simulations are that density-variations in the well-developed wind statistically are quite isotropic and that characteristic length-scales are small; a typical clump size is ~0.01R at 2R, thus resulting also in rather low typical clump-masses ~10^17 g. Overall, our results agree well with the theoretical expectation that the characteristic scale for LDI-generated wind-structure is of order the Sobolev length. We further confirm some earlier results that lateral `filling-in' of radially compressed gas leads to somewhat lower clumping factors in 2D simulations than in comparable 1D models. We conclude by discussing an extension of our method toward rotating LDI wind models that exhibit an intriguing combination of large- and small-scale structure extending down to the wind base. △ Less

Submitted 21 October, 2017; originally announced October 2017.

Comments: 9 pages, 7 figures + 1 Appendix with 1 figure. Recommended for publication in A&A

Journal ref: A&A 611, A17 (2018)

Carbon line formation and spectroscopy in O-type stars

Authors: Luiz P. Carneiro, J. Puls, T. L. Hoffmann

Abstract: The determination of chemical abundances constitutes a fundamental requirement for obtaining a complete picture of a star. Particularly in massive stars, CNO abundances are of prime interest, due to the nuclear CNO-cycle and various mixing processes which bring these elements to the surface. We aim at enabling a reliable carbon spectroscopy for our unified NLTE atmosphere code FASTWIND. We dev… ▽ More The determination of chemical abundances constitutes a fundamental requirement for obtaining a complete picture of a star. Particularly in massive stars, CNO abundances are of prime interest, due to the nuclear CNO-cycle and various mixing processes which bring these elements to the surface. We aim at enabling a reliable carbon spectroscopy for our unified NLTE atmosphere code FASTWIND. We develop a new carbon model atom including CII/III/IV/V, and discuss problems related to carbon spectroscopy in O-type stars. We describe different tests to examine the reliability of our implementation, and investigate which mechanisms influence the carbon ionization balance. By comparing with high-resolution spectra from six O-type stars, we check in how far observational constraints can be reproduced by our new carbon line synthesis. Carbon lines are even more sensitive to a variation of temperature, gravity, and mass-loss rate, than hydrogen/helium lines. We are able to reproduce most of the observed lines from our stellar sample, and to estimate those specific carbon abundances which bring the lines from different ions into agreement. For hot dwarfs and supergiants earlier than O7, X-rays from wind-embedded shocks can impact the synthesized line strengths, particularly for CIV, potentially affecting the abundance determination. We have demonstrated our capability to derive realistic carbon abundances by means of FASTWIND, using our recently developed model atom. We found that complex effects can have a strong influence on the carbon ionization balance in hot stars. For a further understanding, the UV range needs to be explored as well. By means of detailed nitrogen and oxygen model atoms available to use, we will be able to perform a complete CNO abundance analysis for larger samples of massive stars, and to provide constraints on corresponding evolutionary models and aspects. △ Less

Submitted 30 August, 2017; v1 submitted 27 August, 2017; originally announced August 2017.

Comments: 22 pages, 16 figures, 6 tables

Journal ref: A&A 615, A4 (2018)

Stellar Winds in Massive X-ray Binaries

Authors: Peter Kretschmar, Silvia Martínez-Núñez, Enrico Bozzo, Lidia M. Oskinova, Joachim Puls, Lara Sidoli, Jon Olof Sundqvist, Pere Blay, Maurizio Falanga, Felix Fürst, Angel Gímenez-García, Ingo Kreykenbohm, Matthias Kühnel, Andreas Sander, José Miguel Torrejón, Jörn Wilms, Philipp Podsiadlowski, Antonios Manousakis

Abstract: Strong winds from massive stars are a topic of interest to a wide range of astrophysical fields. In High-Mass X-ray Binaries the presence of an accreting compact object on the one side allows to infer wind parameters from studies of the varying properties of the emitted X-rays; but on the other side the accretor's gravity and ionizing radiation can strongly influence the wind flow. Based on a coll… ▽ More Strong winds from massive stars are a topic of interest to a wide range of astrophysical fields. In High-Mass X-ray Binaries the presence of an accreting compact object on the one side allows to infer wind parameters from studies of the varying properties of the emitted X-rays; but on the other side the accretor's gravity and ionizing radiation can strongly influence the wind flow. Based on a collaborative effort of astronomers both from the stellar wind and the X-ray community, this presentation attempts to review our current state of knowledge and indicate avenues for future progress. △ Less

Submitted 22 February, 2017; originally announced February 2017.

Comments: 4 pages, 1 figure, to be published in the Proceedings of the IAU Symposium No. 329 "The lives and death-throes of massive stars"

The VLT-FLAMES Tarantula Survey XXVI: Properties of the O-dwarf population in 30 Doradus

Authors: C. Sabín-Sanjulián, S. Simón-Díaz, A. Herrero, J. Puls, F. R. N. Schneider, C. J. Evans, M. Garcia, F. Najarro, I. Brott, N. Castro, P. A. Crowther, A. de Koter, S. E. de Mink, G. Gräfener, N. J. Grin, G. Holgado, N. Langer, D. J. Lennon, J. Maíz Apellániz, O. H. Ramírez-Agudelo, H. Sana, W. D. Taylor, J. S. Vink, N. R. Walborn

Abstract: The VLT-FLAMES Tarantula Survey has observed hundreds of O-type stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). We study the properties of 105 apparently single O-type dwarfs. To determine stellar and wind parameters, we used the IACOB-GBAT package, an automatic procedure based on a large grid of atmospheric models calculated with the FASTWIND code. In addition to classical tec… ▽ More The VLT-FLAMES Tarantula Survey has observed hundreds of O-type stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). We study the properties of 105 apparently single O-type dwarfs. To determine stellar and wind parameters, we used the IACOB-GBAT package, an automatic procedure based on a large grid of atmospheric models calculated with the FASTWIND code. In addition to classical techniques, we applied the Bayesian BONNSAI tool to estimate evolutionary masses. We provide a new calibration of effective temperature vs. spectral type for O-type dwarfs in the LMC, based on our homogeneous analysis of the largest sample of such objects to date and including all spectral subtypes. Good agreement with previous results is found, although the sampling at the earliest subtypes could be improved. Rotation rates and helium abundances are studied in an evolutionary context. We find that most of the rapid rotators (vsini higher than 300 km/s ) in our sample have masses below 25 MSun and intermediate rotation-corrected gravities (log gc between 3.9 and 4.1). Such rapid rotators are scarce at higher gravities (i.e. younger ages) and absent at lower gravities (larger ages). This is not expected from theoretical evolutionary models, and does not appear to be due to a selection bias in our sample. We compare the estimated evolutionary and spectroscopic masses, finding a trend that the former is higher for masses below 20 MSun. This can be explained as a consequence of limiting our sample to the O-type stars, and we see no compelling evidence for a systematic mass discrepancy. For most of the stars in the sample we were unable to estimate the wind-strength parameter (hence mass-loss rates) reliably, particularly for objects with luminosity lower than logL/LSun about 5.1. Ultraviolet spectroscopy is needed to undertake a detailed investigation of the wind properties of these dwarfs. △ Less

Submitted 15 February, 2017; originally announced February 2017.

Comments: 23 pages, 16 figures. Accepted for publication in A&A

Journal ref: A&A 601, A79 (2017)

Towards a unified view of inhomogeneous stellar winds in isolated supergiant stars and supergiant high mass X-ray binaries

Authors: Silvia Martínez-Núñez, Peter Kretschmar, Enrico Bozzo, Lidia M. Oskinova, Joachim Puls, Lara Sidoli, Jon Olof Sundqvist, Pere Blay, Maurizio Falanga, Felix Fürst, Ángel Gímenez-García, Ingo Kreykenbohm, Matthias Kühnel, Andreas Sander, José Miguel Torrejón, Jörn Wilms

Abstract: Massive stars, at least $\sim$ 10 times more massive than the Sun, have two key properties that make them the main drivers of evolution of star clusters, galaxies, and the Universe as a whole. On the one hand, the outer layers of massive stars are so hot that they produce most of the ionizing ultraviolet radiation of galaxies; in fact, the first massive stars helped to re-ionize the Universe after… ▽ More Massive stars, at least $\sim$ 10 times more massive than the Sun, have two key properties that make them the main drivers of evolution of star clusters, galaxies, and the Universe as a whole. On the one hand, the outer layers of massive stars are so hot that they produce most of the ionizing ultraviolet radiation of galaxies; in fact, the first massive stars helped to re-ionize the Universe after its Dark Ages. Another important property of massive stars are the strong stellar winds and outflows they produce. This mass loss, and finally the explosion of a massive star as a supernova or a gamma-ray burst, provide a significant input of mechanical and radiative energy into the interstellar space. These two properties together make massive stars one of the most important cosmic engines: they trigger the star formation and enrich the interstellar medium with heavy elements, that ultimately leads to formation of Earth-like rocky planets and the development of complex life. The study of massive star winds is thus a truly multidisciplinary field and has a wide impact on different areas of astronomy. [...] This detailed review summarises the current knowledge on the theory and observations of winds from massive stars, as well as on observations and accretion processes in wind-fed high mass X-ray binaries. The aim is to combine in the near future all available theoretical diagnostics and observational measurements to achieve a unified picture of massive star winds in isolated objects and in binary systems. △ Less

Submitted 10 March, 2017; v1 submitted 30 January, 2017; originally announced January 2017.

Comments: Published at Journal of Space Science Reviews, Springer

The VLT-FLAMES Tarantula Survey XXIV. Stellar properties of the O-type giants and supergiants in 30 Doradus

Authors: O. H. Ramírez-Agudelo, H. Sana, A. de Koter, F. Tramper, N. J. Grin, F. R. N. Schneider, N. Langer, J. Puls, N. Markova, J. M. Bestenlehner, N. Castro, P. A. Crowther, C. J. Evans, M. García, G. Gräfener, A. Herrero, B. van Kempen, D. J. Lennon, J. Maíz Apellániz, F. Najarro, C. Sabín-Sanjulián, S. Simón-Díaz, W. D. Taylor, J. S. Vink

Abstract: The Tarantula region in the Large Magellanic Cloud contains the richest population of spatially resolved massive O-type stars known so far. This unmatched sample offers an opportunity to test models describing their main-sequence evolution and mass-loss properties. Using ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to determine stell… ▽ More The Tarantula region in the Large Magellanic Cloud contains the richest population of spatially resolved massive O-type stars known so far. This unmatched sample offers an opportunity to test models describing their main-sequence evolution and mass-loss properties. Using ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to determine stellar, photospheric and wind properties of 72 presumably single O-type giants, bright giants and supergiants and to confront them with predictions of stellar evolution and of line-driven mass-loss theories. We apply an automated method for quantitative spectroscopic analysis of O stars combining the non-LTE stellar atmosphere model {\sc fastwind} with the genetic fitting algorithm {\sc pikaia} to determine the following stellar properties: effective temperature, surface gravity, mass-loss rate, helium abundance, and projected rotational velocity. We present empirical effective temperature versus spectral subtype calibrations at LMC-metallicity for giants and supergiants. In the spectroscopic and classical Hertzsprung-Russell diagrams, our sample O stars are found to occupy the region predicted to be the core hydrogen-burning phase by Brott et al. (2011) and Köhler et al. (2015). Except for five stars, the helium abundance of our sample stars is in agreement with the initial LMC composition. The aforementioned five stars present moderate projected rotational velocities (i.e., $v_{\mathrm{e}}\,\sin\,i\,<\,200\,\mathrm{km\,s^{-1}}$) and hence do not agree with current predictions of rotational mixing in main-sequence stars. Adopting theoretical results for the wind velocity law, we find modified wind momenta for LMC stars that are $\sim$0.3 dex higher than earlier results. [Due to the limitation of characters, the abstract appearing here is slightly shorter than that in the PDF file.] △ Less

Submitted 7 March, 2018; v1 submitted 17 January, 2017; originally announced January 2017.

Comments: 39 pages, 12 figures in the main text, 12 figures in the appendix. Accepted for publication in A&A. Tables C1-C5 will be available at the CDS. Appendix E has been truncated due to size limitations, the full version will be available on A&A

Journal ref: A&A 600, A81 (2017)

Magnetic massive stars as progenitors of "heavy" stellar-mass black holes

Authors: V. Petit, Z. Keszthelyi, R. MacInnis, D. H. Cohen, R. H. D. Townsend, G. A. Wade, S. L. Thomas, S. P. Owocki, J. Puls, J. A. ud-Doula

Abstract: The groundbreaking detection of gravitational waves produced by the inspiralling and coalescence of the black hole (BH) binary GW150914 confirms the existence of "heavy" stellar-mass BHs with masses >25 Msun. Initial modelling of the system by Abbott et al. (2016a) supposes that the formation of black holes with such large masses from the evolution of single massive stars is only feasible if the w… ▽ More The groundbreaking detection of gravitational waves produced by the inspiralling and coalescence of the black hole (BH) binary GW150914 confirms the existence of "heavy" stellar-mass BHs with masses >25 Msun. Initial modelling of the system by Abbott et al. (2016a) supposes that the formation of black holes with such large masses from the evolution of single massive stars is only feasible if the wind mass-loss rates of the progenitors were greatly reduced relative to the mass-loss rates of massive stars in the Galaxy, concluding that heavy BHs must form in low-metallicity (Z < 0.25-0.5 Zsun) environments. However, strong surface magnetic fields also provide a powerful mechanism for modifying mass loss and rotation of massive stars, independent of environmental metallicity (ud-Doula & Owocki 2002; ud-Doula et al. 2008). In this paper we explore the hypothesis that some heavy BHs, with masses >25 Msun such as those inferred to compose GW150914, could be the natural end-point of evolution of magnetic massive stars in a solar-metallicity environment. Using the MESA code, we developed a new grid of single, non-rotating, solar metallicity evolutionary models for initial ZAMS masses from 40-80 Msun that include, for the first time, the quenching of the mass loss due to a realistic dipolar surface magnetic field. The new models predict TAMS masses that are significantly greater than those from equivalent non-magnetic models, reducing the total mass lost by a strongly magnetized 80 Msun star during its main sequence evolution by 20 Msun. This corresponds approximately to the mass loss reduction expected from an environment with metallicity Z = 1/30 Zsun. △ Less

Submitted 27 November, 2016; originally announced November 2016.

Comments: 10 pages, 7 figures, 1 table. Accepted for publication in MNRAS

Modeling the early evolution of massive OB stars with an experimental wind routine

Authors: Zsolt Keszthelyi, Joachim Puls, Gregg Wade

Abstract: Stellar evolution models of massive stars are very sensitive to the adopted mass-loss scheme. The magnitude and evolution of mass-loss rates significantly affect the main sequence evolution, and the properties of post-main sequence objects, including their rotational velocities. Driven by potential discrepancies between theoretically predicted and observationally derived mass-loss rates in the OB… ▽ More Stellar evolution models of massive stars are very sensitive to the adopted mass-loss scheme. The magnitude and evolution of mass-loss rates significantly affect the main sequence evolution, and the properties of post-main sequence objects, including their rotational velocities. Driven by potential discrepancies between theoretically predicted and observationally derived mass-loss rates in the OB star range, we particularly aim to investigate the response to mass-loss rates that are lower than currently adopted, in parallel with the mass-loss behavior at the "first" bi-stability jump. We perform 1D hydrodynamical model calculations of single $20 - 60 \, M_{\odot}$ Galactic ($Z = 0.014$) stars where the effects of stellar winds are already significant during the main sequence phase. We develop an experimental wind routine to examine the behavior and response of the models under the influence of different mass-loss rates. This observationally guided, simple and flexible wind routine is not a new mass-loss description but a useful tool based on the Wind-momentum Luminosity Relation and other scaling relations, and provides a meaningful base for various tests and comparisons. The main result of this study indicates a dichotomy when accounting for currently debated problems regarding mass-loss rates of hot massive stars. In a fully diffusive approach, and for commonly adopted initial rotational velocities, lower mass-loss rates than theoretically predicted require to invoke an additional source of angular momentum loss (either due to bi-stability braking, or yet unidentified) to brake down surface rotational velocities. [...] △ Less

Submitted 14 November, 2016; v1 submitted 15 October, 2016; originally announced October 2016.

Comments: accepted by Astronomy and Astrophysics

Journal ref: A&A 598, A4 (2017)

The VLT-FLAMES Tarantula Survey XXV. Surface nitrogen abundances of O-type giants and supergiants

Authors: N. J. Grin, O. H. Ramirez-Agudelo, A. de Koter, H. Sana, J. Puls, I. Brott, P. A. Crowther, P. L. Dufton, C. J. Evans, G. Graefener, A. Herrero, N. Langer, D. J. Lennon, J. Th. van Loon, N. Markova, S. E. de Mink, F. Najarro, F. R. N. Schneider, W. D. Taylor, F. Tramper, J. S. Vink, N. R. Walborn

Abstract: Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Using a sample of 72 presumably single O-type giants to supergiants observed in the context of the VLT-FLAMES Tarantula Survey (VFTS), we aim to investigate rotational mixing in evolve… ▽ More Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Using a sample of 72 presumably single O-type giants to supergiants observed in the context of the VLT-FLAMES Tarantula Survey (VFTS), we aim to investigate rotational mixing in evolved core-hydrogen burning stars initially more massive than $15\,M_\odot$ by analysing their surface nitrogen abundances. Using stellar and wind properties derived in a previous VFTS study, we constrained the nitrogen abundance by fitting the equivalent widths of relatively strong lines that are sensitive to changes in the abundance of this element. Given the quality of the data, we constrained the nitrogen abundance in 38 cases; for 34 stars only upper limits could be derived, which includes almost all stars rotating at $v_\mathrm{e}\sin i >200\,\mathrm{km s^{-1}}$. We analysed the nitrogen abundance as a function of projected rotation rate $v_\mathrm{e}\sin i$ and confronted it with predictions of rotational mixing. The upper limits on the nitrogen abundance of the rapidly rotating stars are not in apparent violation with theoretical expectations. However, we found a group of N-enhanced slowly-spinning stars that is not in accordance with predictions of rotational mixing in single stars. Among O-type stars with (rotation-corrected) gravities less than $\log\,g_c = 3.75$ this group constitutes 30$-$40 percent of the population. We found a correlation between nitrogen and helium abundance which is consistent with expectations, suggesting that, whatever the mechanism that brings N to the surface, it displays CNO-processed material. △ Less

Submitted 1 September, 2016; originally announced September 2016.

Journal ref: A&A 600, A82 (2017)

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

Authors: S. Simón-Díaz, M. Godart, N. Castro, A. Herrero, C. Aerts, J. Puls, J. Telting, L. Grassitelli

Abstract: 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 spectroscopi… ▽ More 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. △ Less

Submitted 19 August, 2016; originally announced August 2016.

Comments: 17 pages, 11 figures, 3 tables. Accepted for publication in A&A

Journal ref: A&A 597, A22 (2017)

The SILCC project: III. Regulation of star formation and outflows by stellar winds and supernovae

Authors: A. Gatto, S. Walch, T. Naab, P. Girichidis, R. Wünsch, S. C. O. Glover, R. S. Klessen, P. C. Clark, T. Peters, D. Derigs, C. Baczynski, J. Puls

Abstract: We study the impact of stellar winds and supernovae on the multi-phase interstellar medium using three-dimensional hydrodynamical simulations carried out with FLASH. The selected galactic disc region has a size of (500 pc)$^2$ x $\pm$ 5 kpc and a gas surface density of 10 M$_{\odot}$/pc$^2$. The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse h… ▽ More We study the impact of stellar winds and supernovae on the multi-phase interstellar medium using three-dimensional hydrodynamical simulations carried out with FLASH. The selected galactic disc region has a size of (500 pc)$^2$ x $\pm$ 5 kpc and a gas surface density of 10 M$_{\odot}$/pc$^2$. The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse heating, sink particles representing star clusters, stellar winds from these clusters which combine the winds from indi- vidual massive stars by following their evolution tracks, and subsequent supernova explosions. Dust and gas (self-)shielding is followed to compute the chemical state of the gas with a chemical network. We find that stellar winds can regulate star (cluster) formation. Since the winds suppress the accretion of fresh gas soon after the cluster has formed, they lead to clusters which have lower average masses (10$^2$ - 10$^{4.3}$ M$_{\odot}$) and form on shorter timescales (10$^{-3}$ - 10 Myr). In particular we find an anti-correlation of cluster mass and accretion time scale. Without winds the star clusters easily grow to larger masses for ~5 Myr until the first supernova explodes. Overall the most massive stars provide the most wind energy input, while objects beginning their evolution as B-type stars contribute most of the supernova energy input. A significant outflow from the disk (mass loading $\gtrsim$ 1 at 1 kpc) can be launched by thermal gas pressure if more than 50% of the volume near the disc mid-plane can be heated to T > 3x10$^5$ K. Stellar winds alone cannot create a hot volume-filling phase. The models which are in best agreement with observed star formation rates drive either no outflows or weak outflows. △ Less

Submitted 16 June, 2016; originally announced June 2016.

Comments: 23 pages; submitted to MNRAS

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

Authors: Paul A. Crowther, S. M. Caballero-Nieves, K. A. Bostroem, J. Maiz Apellaniz, F. R. N. Schneider, N. R. Walborn, C. R. Angus, I. Brott, A. Bonanos, A. de Koter, S. E. de Mink, C. J. Evans, G. Grafener, A. Herrero, I. D. Howarth, N. Langer, D. J. Lennon, J. Puls, H. Sana, J. S. Vink

Abstract: 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… ▽ More 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. △ Less

Submitted 16 March, 2016; originally announced March 2016.

Comments: 20 pages plus four Appendices providing LMC UV O spectral templates, UV spectral atlas in R136, wind velocities of LMC O stars and photometry of additional R136 sources

Journal ref: MNRAS 458, 624-659 (2016)

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

Authors: Luiz P. Carneiro, J. Puls, J. O. Sundqvist, T. L. Hoffmann

Abstract: 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. Th… ▽ More 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. △ Less

Submitted 3 March, 2016; originally announced March 2016.

Comments: accepted by A&A

Journal ref: A&A 590, A88 (2016)