Search

The cluster R136 in the giant star-forming region 30 Doradus in the Large Magellanic Cloud (LMC) offers a unique opportunity to resolve a stellar population in a starburst-like environment. We obtain the near-infrared to ultraviolet extinction towards 50 stars in the core of R136, employing the `extinction without standards' method. To assure good fits over the full wavelength range, we combine and modify existing extinction laws. We detect a strong spatial gradient in the extinction properties across the core of R136, coinciding with a gradient in density of cold gas that is part of a molecular cloud lying northeast of the cluster. In line with previous measurements of R136 and the 30 Doradus region, we obtain a high total-to-relative extinction ($R_V = 4.38 \pm 0.87$). However, the high values of $R_V$ are accompanied by relatively strong extinction in the ultraviolet, contrary to what is observed for Galactic sightlines. The relatively strong ultraviolet extinction suggests that the properties of the dust towards R136 differ from those in the Milky Way. For $R_{V} \sim 4.4$, about three times fewer ultraviolet photons can escape from the ambient dust environment relative to the canonical Galactic value of $R_{V} \sim 3.1$ at the same $A_{V}$. Therefore, if dust in the R136 star-bursting environment is characteristic for cosmologically distant star-bursting regions, the escape fraction of ultraviolet photons from such regions is overestimated by a factor of three relative to the standard Milky Way assumption for the total-to-selective extinction. Furthermore, a comparison with average curves tailored to other regions of the LMC shows that large differences in ultraviolet extinction exist within this galaxy. Further investigation is required in order to decipher whether or not there is a relation between $R_V$ and ultraviolet extinction in the LMC., Accepted for publication in A&A (18 pages, 12 figures)

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\v{c}ka & Kub\'at (2018) matches best the observed mass-loss rates., Comment: Accepted for publication in A&A; Appendix I will not be included in the published version

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only a few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray-faint binary in the Large Magellanic Cloud. With an orbital period of 10.4 d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5σ confidence level. The minimum companion mass implies that it is a black hole. No other X-ray-quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the cosmos., This research has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (T.S. and H.S., grant agreement 772225: MULTIPLES). T.S. acknowledges support from the European Union’s Horizon 2020 programme under Marie Skłodowska-Curie grant agreement 101024605. This work is based on observations collected at the ESO under programme IDs 182.D-0222, 090.D-0323 and 092.D-0136. L.M. thanks the ESA and the Belgian Federal Science Policy Office (BELSPO) for their support in the framework of the PRODEX programme. P. Marchant acknowledges support from FWO junior postdoctoral fellowship 12ZY520N. We thank J. Hillier for making the CMFGEN code available. C.H. acknowledges support from the KU Leuven Research Council (grant C16/17/007: MAESTRO). P.A.C. acknowledges support from UK Science and Technology Facilities Council research grant ST/V000853/1. V.H.-B. acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC) through grant RGPIN-2020-05990. M. Gieles acknowledges support from the Ministry of Science and Innovation through a Europa Excelencia grant (EUR2020-112157). The PoWR code, as well as the associated post-processing and visualization tool WRplot, has been developed under the guidance of W.-R. Hamann with substantial contributions from L. Koesterke, G. Gräfener, A. Sander, T.S. and other co-workers and students. This work has received funding from the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement 945806). It is also supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy EXC 2181/1-390900948 (the Heidelberg STRUCTURES Excellence Cluster). S.T. acknowledges support from the Netherlands Research Council NWO (grants VENI 639.041.645, VIDI 203.061). A.H. and D.J.L. acknowledge support by the Spanish MCI through grant PGC-2018-0913741-B-C22 and the Severo Ochoa Programme through CEX2019-000920-S. J.M.A. acknowledges support from the Spanish Government Ministerio de Ciencia, Innovación y Universidades through grant PGC2018-095-049-B-C22.

Stars strongly impact their environment, and shape structures on all scales throughout the universe, in a process known as “feedback.” Due to the complexity of both stellar evolution and the physics of larger astrophysical structures, there remain many unanswered questions about how feedback operates and what we can learn about stars by studying their imprint on the wider universe. In this white paper, we summarize discussions from the Lorentz Center meeting “Bringing Stellar Evolution and Feedback Together” in 2022 April and identify key areas where further dialog can bring about radical changes in how we view the relationship between stars and the universe they live in.

A pioneering study showed that the fine structure in the luminosity function (LF) of young star clusters contains information about the evolutionary stage (age) and composition of the stellar population. The notable features include the H-peak, which is the result of the onset of hydrogen burning turning pre-main sequence stars into main sequence stars. The feature moves toward the faint end of the LF, and eventually disappears as the population evolves. Another detectable feature is the Wielen dip, a dip at M_V ~ 7 mag in the LF first identified in 1974 for stars in the solar environment. Later studies also identified this feature in the LF of star clusters. The Wielen dip is caused by the increased importance of H- opacity in a certain range of low-mass stars. We studied the detailed structure in the luminosity function using the data from Gaia DR2 and PARSEC stellar evolution models with the aim to further our understanding of young stellar populations. We analyzed the astrometric properties of stars in the solar neighborhood (< 20 pc) and in various relatively nearby (< 400 pc) young (< 50 Myr) open clusters and OB associations, and compare the features in the luminosity function with those generated by PARSEC models. The Wielen dip is confirmed in the LF of all the populations, including the solar neighborhood, at M_G ~7 mag. The H-peak is present in the LF of the field stars in the solar neighborhood. It likely signals that the population is mixed with a significant number of stars younger than 100 Myr. The H-peak is found in the LF of young open clusters and OB associations, and its location varies with age. Our observations with Gaia DR2 confirm the evolution of the H-peak from 5 Myr up to 47 Myr. The fine structure in the luminosity function in young stellar populations can be used to estimate their age., 16 pages, A&A in press

The majority of close massive binary stars with initial periods of a few days experience a contact phase, in which both stars overflow their Roche lobes simultaneously. We perform the first dedicated study of the evolution of massive contact binaries and provide a comprehensive prediction of their observed properties. We compute 2790 detailed binary models for the Large and Small Magellanic Clouds each, assuming a conservative mass transfer. The initial parameter space for both grids span total masses from 20 to 80$\,\textrm{M}_{\odot}$, orbital periods of 0.6 to 2 days and mass ratios of 0.6 to 1.0. We find that models that remain in contact over nuclear timescales evolve towards equal masses, echoing the mass ratios of their observed counterparts. Ultimately, the fate of our nuclear-timescale models is to merge on the main sequence. Our predicted period-mass ratio distributions of O-type contact binaries are similar for both galaxies, and we expect 10 such systems together in both Magellanic Clouds. While we can largely reproduce the observed distribution, we over-estimate the population of equal-mass contact binaries. This situation is somewhat remedied if we also account for binaries that are approaching contact. Our theoretical distributions work particularly well for contact binaries with periods $, Comment: 23 pages, 11 figures, accepted for publication in MNRAS

Observations of massive stars in young open clusters (< ~8 Myr) have shown that a majority of them are in binary systems, most of which will interact during their life. Populations of massive stars older than ~20 Myr allow us to probe the outcome of such interactions after many systems have experienced mass and angular momentum transfer. Using multi-epoch integral-field spectroscopy, we investigate the multiplicity properties of the massive-star population in NGC 330 (~40 Myr) in the Small Magellanic Cloud to search for imprints of stellar evolution on the multiplicity properties. From six epochs of VLT/MUSE observations supported by adaptive optics we extract spectra and measure radial velocities for stars brighter than F814W = 19. We identify single-lined spectroscopic binaries through significant RV variability as well as double-lined spectroscopic binaries, and quantify the observational biases for binary detection. The observed spectroscopic binary fraction is 13.2+/-2.0 %. Considering period and mass ratio ranges from log(P)=0.15-3.5, and q = 0.1-1.0, and a representative set of orbital parameter distributions, we find a bias-corrected close binary fraction of 34 +8 -7 %. This seems to decline for the fainter stars, which indicates either that the close binary fraction drops in the B-type domain, or that the period distribution becomes more heavily weighted towards longer orbital periods. Both fractions vary strongly in different regions of the color-magnitude diagram which probably reveals the imprint of the binary history of different groups of stars. We provide the first homogeneous RV study of a large sample of B-type stars at a low metallicity. The overall bias-corrected close binary fraction of B stars in NGC 330 is lower than the one reported for younger Galactic and LMC clusters. More data are needed to establish whether this result from an age or a metallicty effect., 19 page (incl. appendix), 10 figures, 4 tables, accepted for publication in A&A

Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we observed the M-type asymptotic giant branch (AGB) star R Hydrae with ALMA. The morphology of the inner wind of R Hya, which has a known companion at ~3500 au, was determined from maps of CO and SiO obtained at high angular resolution. A map of the CO emission reveals a multi-layered structure consisting of a large elliptical feature at an angular scale of ~10'' that is oriented along the north-south axis. The wind morphology within the elliptical feature is dominated by two hollow bubbles. The bubbles are on opposite sides of the AGB star and lie along an axis with a position angle of ~115 deg. Both bubbles are offset from the central star, and their appearance in the SiO channel maps indicates that they might be shock waves travelling through the AGB wind. An estimate of the dynamical age of the bubbles yields an age of the order of 100 yr, which is in agreement with the previously proposed elapsed time since the star last underwent a thermal pulse. When the CO and SiO emission is examined on subarcsecond angular scales, there is evidence for an inclined, differentially rotating equatorial density enhancement, strongly suggesting the presence of a second nearby companion. The position angle of the major axis of this disc is ~70 deg in the plane of the sky. We tentatively estimate that a lower limit on the mass of the nearby companion is ~0.65 Msol on the basis of the highest measured speeds in the disc and the location of its inner rim at ~6 au from the AGB star., 21 pages, 23 figures

In this paper we make predictions for the behaviour of wind bubbles around young massive stars using analytic theory. We do this in order to determine why there is a discrepancy between theoretical models that predict that winds should play a secondary role to photoionisation in the dynamics of HII regions, and observations of young HII regions that seem to suggest a driving role for winds. In particular, regions such as M42 in Orion have neutral hydrogen shells, suggesting that the ionising radiation is trapped closer to the star. We first derive formulae for wind bubble evolution in non-uniform density fields, focusing on singular isothermal sphere density fields with a power law index of -2. We find that a classical "Weaver"-like expansion velocity becomes constant in such a density distribution. We then calculate the structure of the photoionised shell around such wind bubbles, and determine at what point the mass in the shell cannot absorb all of the ionising photons emitted by the star, causing an "overflow" of ionising radiation. We also estimate perturbations from cooling, gravity, magnetic fields and instabilities, all of which we argue are secondary effects for the conditions studied here. Our wind-driven model provides a consistent explanation for the behaviour of M42 to within the errors given by observational studies. We find that in relatively denser molecular cloud environments \around single young stellar sources, champagne flows are unlikely until the wind shell breaks up due to turbulence or clumping in the cloud., Comment: 17 pages, 10 figures, published in MNRAS

We present a spectroscopic analysis of the most rapidly rotating stars currently known, VFTS 102 ($v_{e} \sin i = 649 \pm 52$ km s$^{-1}$; O9: Vnnne+) and VFTS 285 ($v_{e} \sin i = 610 \pm 41$ km s$^{-1}$; O7.5: Vnnn), both members of the 30 Dor complex in the Large Magellanic Cloud. This study is based on high resolution ultraviolet spectra from HST/COS and optical spectra from VLT X-shooter plus archival VLT GIRAFFE spectra. We utilize numerical simulations of their photospheres, rotationally distorted shape, and gravity darkening to calculate model spectral line profiles and predicted monochromatic absolute fluxes. We use a guided grid search to investigate parameters that yield best fits for the observed features and fluxes. These fits produce estimates of the physical parameters for these stars (plus a Galactic counterpart, $\zeta$ Oph) including the equatorial rotational velocity, inclination, radius, mass, gravity, temperature, and reddening. We find that both stars appear to be radial velocity constant. VFTS 102 is rotating at critical velocity, has a modest He enrichment, and appears to share the motion of the nearby OB association LH 99. These properties suggest that the star was spun up through a close binary merger. VFTS 285 is rotating at $95\%$ of critical velocity, has a strong He enrichment, and is moving away from the R136 cluster at the center of 30 Dor. It is mostly likely a runaway star ejected by a supernova explosion that released the components of the natal binary system., Comment: 28 pages, 12 figures, 4 tables

The nebular circumstellar environments of cool evolved stars are known to harbour a rich morphological complexity of gaseous structures on different length scales. A large part of these density structures are thought to be brought about by the interaction of the stellar wind with a close companion. The S-type asymptotic giant branch star Pi1 Gruis, which has a known companion at ~440 au and is thought to harbour a second, closer-by (, 19 pages, 20 figures

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 $\gamma \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., Comment: Accepted for publication in MNRAS, 21 pages, 14 figures, 2 tables plus 33 pages of supplementary material

The majority of massive stars, the progenitors of core-collapse supernovae (SNe), are found in close binary systems. Zapartas et al. (2019) modeled the fraction of hydrogen-rich, Type II SN progenitors which have their evolution affected by mass exchange with their companion, finding this to be between 1/3 and 1/2 for most assumptions. Here we study in more depth the impact of this binary history of Type II SN progenitors on their final pre-SN core mass distribution, using population synthesis simulations. We find that binary star progenitors of Type II SNe typically end their life with a larger core mass than they would have had if they had lived in isolation, because they gained mass or merged with a companion before explosion. The combination of the diverse binary evolutionary paths typically lead to a marginally shallower final core mass distribution. Discussing our results in the context of the red supergiant problem, i.e., the reported lack of detected high luminosity progenitors, we conclude that binary evolution does not seem to significantly affect the issue. This conclusion is quite robust against our variations in the assumptions of binary physics. We also predict that inferring the initial masses of Type II SN progenitors from "age-dating" its surrounding environment systematically yields lower masses compared to methods that probe the pre-SN core mass or luminosity. A robust discrepancy between the inferred initial masses of a SN progenitor from those different techniques could indicate an evolutionary history of binary mass accretion or merging., Published in Astronomy & Astrophysics, Volume 645

Winds from young massive stars contribute a large amount of energy to their host molecular clouds. This has consequences for the dynamics and observable structure of star-forming clouds. In this paper, we present radiative magnetohydrodynamic simulations of turbulent molecular clouds that form individual stars of 30, 60 and 120 solar masses emitting winds and ultraviolet radiation following realistic stellar evolution tracks. We find that winds contribute to the total radial momentum carried by the expanding nebula around the star at 10 % of the level of photoionisation feedback, and have only a small effect on the radial expansion of the nebula. Radiation pressure is largely negligible in the systems studied here. The 3D geometry and evolution of wind bubbles is highly aspherical and chaotic, characterised by fast-moving "chimneys" and thermally-driven "plumes". These plumes can sometimes become disconnected from the stellar source due to dense gas flows in the cloud. Our results compare favourably with the findings of relevant simulations, analytic models and observations in the literature while demonstrating the need for full 3D simulations including stellar winds. However, more targeted simulations are needed to better understand results from observational studies., Comment: 18 pages, 11 figures, accepted to MNRAS

Hydrogen-rich supernovae, known as Type II (SNe II), are the most common class of explosions observed following the collapse of the core of massive stars. We use analytical estimates and population synthesis simulations to assess the fraction of SNe II progenitors that are expected to have exchanged mass with a companion prior to explosion. We estimate that 1/3 to 1/2 of SN II progenitors have a history of mass exchange with a binary companion before exploding. The dominant binary channels leading to SN II progenitors involve the merger of binary stars. Mergers are expected to produce a diversity of SN II progenitor characteristics, depending on the evolutionary timing and properties of the merger. Alternatively, SN II progenitors from interacting binaries may have accreted mass from their companion, and subsequently been ejected from the binary system after their companion exploded. We show that the overall fraction of SN II progenitors that are predicted to have experienced binary interaction is robust against the main physical uncertainties in our models. However, the relative importance of different binary evolutionary channels is affected by changing physical assumptions. We further discuss ways in which binarity might contribute to the observed diversity of SNe II by considering potential observational signatures arising from each binary channel. For supernovae which have a substantial H-rich envelope at explosion (i.e., excluding Type IIb SNe), a surviving non-compact companion would typically indicate that the supernova progenitor star was in a wide, non-interacting binary. We argue that a significant fraction of even Type II-P SNe are expected to have gained mass from a companion prior to explosion., The paper has been accepted for publication in Astronomy and Astrophysics

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., 24 pages, 4 tables, 12 figures, accepted for publication in ApJ

High-resolution interferometric observations of the circumstellar environments of AGB stars show a variety of morphologies. Guided by the unusual carbon monoxide line profile of the AGB star EP Aquarii, we have observed its circumstellar environment with ALMA band 6 in cycle 4. We describe the morphological complexity of the CO, SiO and SO2 molecular emission. The CO emission exhibits the characteristics of a bi-conical wind with a bright nearly face-on spiral feature around the systemic velocity. This is the first convincing detection of a spiral morphology in an O-rich wind. Based on the offsets of the centers of the two bi-conical wind hemispheres, we deduce the position angle of the inclination axis to be ~150 deg measured counterclockwise from north. Based on the velocity width of the spiral signature we estimate the inclination angle of the system to be between 4 deg and 18 deg. The central emission zone exhibits a morphology that resembles simulations modelling the spiral-inducing wind-Roche-lobe-overflow mechanism. Though the spiral may be a (companion-induced) density enhancement in the stellar outflow, the extremely narrow width of the spiral signature in velocity space suggests that it may be a hydrodynamical perturbation in a face-on differentially rotating disk. The SiO emission does not show the spiral, but exhibits a local emission void approximately 0.5'' west of the continuum brightness peak. We hypothesise that this may be a local environment caused by the presence of a (stellar) companion with a mass of at most 0.1 Msol, based on its non-detection in the continuum..., 20 pages, 16 Figures

A spectral scan of the circumstellar environment of the asymptotic giant branch (AGB) star R~Doradus was taken with ALMA in cycle 2 at frequencies between 335 and 362 GHz and with a spatial resolution of $\sim$150 milliarcseconds. Many molecular lines show a spatial offset between the blue and red shifted emission in the innermost regions of the wind. The position-velocity diagrams of this feature, in combination with previous SPHERE data and theoretical work point towards the presence of a compact differentially rotating disk, orientated nearly edge-on. We model the $^{\rm 28}$SiO ($v=1,~J=8\to7$) emission with a disk model. We estimate the disk mass and angular momentum to be $3 \times 10^{-6}$ Solar masses and $5 \times 10^{40}\ {\rm m^2 kg/s}$. The latter presents an `angular momentum problem' that may be solved by assuming that the disk is the result of wind-companion interactions with a companion of at least 2.5 earth masses, located at 6 AU, the tentatively determined location of the disk's inner rim. An isolated clump of emission is detected to the south-east with a velocity that is high compared to the previously determined terminal velocity of the wind. Its position and mean velocity suggest that it may be associated with a companion planet, located at the disk's inner rim., Comment: 11 pages, 7 Figures

The circumstellar environment of the AGB star L$_{\rm 2}$ Puppis was observed with ALMA in cycle 3, with a resolution of $15 \times 18 \rm\ mas$. The molecular emission shows a differentially rotating disk, inclined to a nearly edge-on position. In the first paper in this series (paper I) the molecular emission was analysed to accurately deduce the motion of the gas in the equatorial regions of the disk. In this work we model the optically thick $^{\rm 12}$CO $J=$3$-$2 and the optically thin $^{\rm 13}$CO $J=$3$-$2 rotational transition to constrain the physical conditions in the disk. To realise this effort we make use of the 3D NLTE radiative transfer code {\tt LIME}. The temperature structure and velocity structure show a high degree of complexity, both radially and vertically. The radial H$_{\rm 2}$ density profile in the disk plane is characterised by a power law with a slope of $-3.1$. We find a $^{\rm 12}$CO over $^{\rm 13}$CO abundance ratio of 10 inside the disk. Finally, estimations of the angular momentum in the disk surpass the expected available angular momentum of the star, strongly supporting the indirect detection of a compact binary companion reported in paper I. We estimate the mass of the companion to be around 1 Jupiter mass., 21 pages; 19 figures

The carriers of the diffuse interstellar bands (DIBs) are largely unidentified molecules ubiquitously present in the interstellar medium (ISM). After decades of study, two strong and possibly three weak near-infrared DIBs have recently been attributed to the C60+ fullerene based on observational and laboratory measurements. There is great promise for the identification of the over 400 other known DIBs, as this result could provide chemical hints towards other possible carriers. In an effort to systematically study the properties of the DIB carriers, we have initiated a new large-scale observational survey: the ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES). The main objective is to build on and extend existing DIB surveys to make a major step forward in characterising the physical and chemical conditions for a statistically significant sample of interstellar lines-of-sight, with the goal to reverse-engineer key molecular properties of the DIB carriers. EDIBLES is a filler Large Programme using the Ultraviolet and Visual Echelle Spectrograph at the Very Large Telescope at Paranal, Chile. It is designed to provide an observationally unbiased view of the presence and behaviour of the DIBs towards early-spectral-type stars whose lines-of-sight probe the diffuse-to-translucent ISM. Such a complete dataset will provide a deep census of the atomic and molecular content, physical conditions, chemical abundances and elemental depletion levels for each sightline. Achieving these goals requires a homogeneous set of high-quality data in terms of resolution (R ~ 70000 -- 100000), sensitivity (S/N up to 1000 per resolution element), and spectral coverage (305--1042 nm), as well as a large sample size (100+ sightlines). In this first paper the goals, objectives and methodology of the EDIBLES programme are described and an initial assessment of the data is provided., 29 pages (including 14 pages appendix), 15 figures, accepted by A&A

International audience

We report on a comprehensive study of the wind properties of 115 O- and early B-type stars in the Galaxy and the Large Magellanic Clouds. This work is part of the VLT/FLAMES Survey of Massive Stars. The data is used to construct the empirical dependence of the mass-loss in stellar winds on the metal content of their atmospheres. The metal content of early-type stars in the Magellanic Clouds is discussed. Assuming a power-law dependence of mass loss on metal content, Ṁ ∝ Zm, we find m = 0.83 ± 0.16 from an analysis of the wind momentum luminosity relation (Mokiem et al. 2007b). This result is in good agreement with the prediction m = 0.69 ± 0.10 by Vink et al. (2001). Though the scaling agrees, the absolute empirical value of mass loss is found to be a factor of two higher than predictions. This may be explained by a modest amount of clumping in the outflows of the objects studied.

Context. Recent high-resolution observations have shown that stellar winds harbour complexities that strongly deviate from spherical symmetry, which generally is assumed as standard wind model. One such morphology is the Archimedean spiral, which is generally believed to be formed by binary interactions, as has been directly observed in multiple sources.Aims. We seek to investigate the manifestation in the observables of spiral structures embedded in the spherical outflows of cool stars. We aim to provide an intuitive bedrock with which upcoming ALMA data can be compared and interpreted.Methods. By means of an extended parameter study, we modelled rotational CO emission from the stellar outflow of asymptotic giant branch stars. To this end, we developed a simplified analytical parametrised description of a 3D spiral structure. This model is embedded into a spherical wind and fed into the 3D radiative transfer code LIME, which produces 3D intensity maps throughout velocity space. Subsequently, we investigated the spectral signature of rotational transitions of CO in the models, as well as the spatial aspect of this emission by means of wide-slit position-velocity (PV) diagrams. Additionally, we quantified the potential for misinterpreting the 3D data in a 1D context. Finally, we simulated ALMA observations to explore the effect of interferometric noise and artefacts on the emission signatures.Results. The spectral signatures of the CO rotational transition v = 0J = 3 − 2 are very efficient at concealing the dual nature of the outflow. Only a select few parameter combinations allow for the spectral lines to disclose the presence of the spiral structure. If the spiral cannot be distinguished from the spherical signal, this might result in an incorrect interpretation in a 1D context. Consequently, erroneous mass-loss rates would be calculated. The magnitude of these errors is mainly confined to a factor of a few, but in extreme cases can exceed an order of magnitude. CO transitions of different rotationally excited levels show a characteristical evolution in their line shape that can be brought about by an embedded spiral structure. However, if spatial information on the source is also available, the use of wide-slit PV diagrams systematically expose the embedded spiral. The PV diagrams also readily provide most of the geometrical and physical properties of the spiral-harbouring wind. Simulations of ALMA observations prove that the choice of antenna configuration is strongly dependent on the geometrical properties of the spiral. We conclude that exploratory endeavours should observe the object of interest with a range of different maximum-baseline configurations.

We present radiation-driven wind models for Luminous Blue Variables (LBVs) and predict their mass-loss rates. We study the effects of lower masses and modified abundances in comparison to the normal OB supergiants, and we find that the main difference in mass loss is due to the lower masses of LBVs. In addition, we find that the increase in helium abundance changes the mass-loss properties by small amounts (up to about 0.2 dex in log Mdot), while CNO processing is relatively unimportant for the mass-loss rate. A comparison between our mass loss predictions and the observations is performed for four relatively well-studied LBVs. The comparison shows that (i) the winds of LBVs are driven by radiation pressure on spectral lines, (ii) the variable mass loss behaviour of LBVs during their S Doradus-type variation cycles is explained by changes in the line driving efficiency, notably due to the recombination/ionisation of Fe IV/III and Fe III/II, and finally, (iii) the winds of LBVs can be used to derive their masses, as exemplified by the case of AG Car, for which we derive a present-day mass of 35 Msun., Comment: 12 pages; A&A accepted

Powerful telescopes equipped with multi-fibre or integral field spectrographs combined with detailed models of stellar atmospheres and automated fitting techniques allow for the analysis of large number of stars. These datasets contain a wealth of information that require new analysis techniques to bridge the gap between observations and stellar evolution models. To that end, we develop BONNSAI (BONN Stellar Astrophysics Interface), a Bayesian statistical method, that is capable of comparing all available observables simultaneously to stellar models while taking observed uncertainties and prior knowledge such as initial mass functions and distributions of stellar rotational velocities into account. BONNSAI can be used to (1) determine probability distributions of fundamental stellar parameters such as initial masses and stellar ages from complex datasets, (2) predict stellar parameters that were not yet observationally determined and (3) test stellar models to further advance our understanding of stellar evolution. An important aspect of BONNSAI is that it singles out stars that cannot be reproduced by stellar models through $\chi^{2}$ hypothesis tests and posterior predictive checks. BONNSAI can be used with any set of stellar models and currently supports massive main-sequence single star models of Milky Way and Large and Small Magellanic Cloud composition. We apply our new method to mock stars to demonstrate its functionality and capabilities. In a first application, we use BONNSAI to test the stellar models of Brott et al. (2011a) by comparing the stellar ages inferred for the primary and secondary stars of eclipsing Milky Way binaries. Ages are determined from dynamical masses and radii that are known to better than 3%. We find that the stellar models reproduce the Milky Way binaries well. BONNSAI is available through a web-interface at http://www.astro.uni-bonn.de/stars/bonnsai., Comment: Accepted for publication in A&A; 15 pages, 10 figures, 4 tables; BONNSAI is available through a web-interface at http://www.astro.uni-bonn.de/stars/bonnsai

We have calculated a grid of massive star wind models and mass-loss rates for a wide range of metal abundances between 1/100 and 10 Z/Zsun. The calculation of this grid completes the Vink et al. (2000) mass-loss recipe with an additional parameter Z. We have found that the exponent of the power law dependence of mass loss vs. metallicity is constant in the range between 1/30 and 3 Z/Zsun. The mass-loss rate scales as Mdot \propto Z^0.85 Vinf^p with p = -1.23 for stars with Teff \ga 25000 K, and p = -1.60 for the B supergiants with Teff \la 25000 K. Taking also into account the metallicity dependence of Vinf, using the power law dependence Vinf \propto Z^0.13 from Leitherer et al. (1992), the overall result of mass loss as a function of metallicity can be represented by Mdot \propto Z^0.69 for stars with Teff \ga 25000 K, and Mdot \propto Z^0.64 for B supergiants with Teff \la 25000 K. Our mass-loss predictions are successful in explaining the observed mass-loss rates for Galactic and Small Magellanic Cloud O-type stars, as well as in predicting the observed Galactic bi-stability jump. Hence, we believe that our predictions are reliable and suggest that our mass-loss recipe be used in future evolutionary calculations of massive stars at different metal abundance. A computer routine to calculate mass loss is publicly available., Accepted by A&A, 16 pages, (email: jvink@ic.ac.uk)

Empirical ionization fractions of C IV ,N V ,S iIV, and empirical ionization plus excitation fractions of C III* and N IV* in the winds of 34 O stars and one B star have been derived. We combine the mass-loss rates derived from radio measurements and Ha with the line —tting of ultraviolet resonance lines and subordinate lines using the Sobolev plus exact integration (SEI) method. The dependence of the empiri- cal ionization fractions, SqT, on the stellar eUective temperature and on the mean wind density is dis- cussed. This sets constraints for the models of ionization in the winds of hot stars. The ionization and excitation fractions can be expressed in terms of an empirical radiation temperature. This radiation tem- perature scales with and we derive empirical relations for as a function of The radiation T eff , T rad T eff . temperatures are on the order of 0.5¨0.9 with signi—cant diUerences between the ions. The derived T eff , relations between the ionization fractions and the stellar parameters have an uncertainty of 0.2 dex for Si IV ,N V, and C III*, and about 0.26 dex for N IV*. For C IV, we can only derive an expression for the mean ionization fraction in the wind if the mass-loss rate is small, yr~1, because the C IV M 0\ 10~6 M _ lines are usually saturated for higher mass-loss rates. The resulting expressions for can be used to T rad derive the mass-loss rates from studies of ultraviolet P Cygni pro—les in the range of stellar parameters studied here: 30,000 K, and ( 7.5 K ( T eff ( 50,500 5.2 ( log L /L * ( 6.4, M _ ( log M0 ( (4.6 M _ yr~1. An accuracy of about a factor of 2 or better can be reached, depending on the lines that are used and the accuracy of the line —ts and the stellar parameters. The Si IV lines give the most reliable mass- loss rates, because the abundance is about the same for all O stars, the lines saturate only for high mass-loss rates, the doublet lines only partly overlap, and the mass-loss rate is proportional to the square root of the column density. The radiation temperature of N V shows a surprisingly strict relation with with a scatter of only The mass-loss rate cannot be derived from the N V T eff , *T rad /T eff \ 0.01. lines, because the column density of the N V ions in the wind is independent of A consistency check M0 . and a test of the method for the stars HD 14749 and HD 190429 show that the mass-loss rate derived from the UV lines with the ionization fractions of this paper agree very well with the mass-loss rate derived from new radio —ux measurements. Subject headings: stars: atmospheresstars: early-typestars: mass lossultraviolet: stars

We report on a detailed spectroscopic study of 12 very massive and luminous stars (M 35M?) in the core of the compact cluster R136a, near the center of the 30 Doradus complex. The three brightest stars of the cluster, R136a1, R136a2, and R136a3, have been investigated earlier by de Koter, Heap, & Hubeny. Low-resolution spectra (

We report on a systematic study of the most massive stars, in which we analyzed the spectra of four very luminous stars in the Large Magellanic Cloud. The stars lie in the 30 Doradus complex, three of which are located in the core of the compact cluster, R136a (R136a1, R136a3, and R136a5), and the fourth (Melnick 42), located about 8'' north of R136a. Low-resolution spectra (

This paper represents a collective effort to provide an extensive electronic data base useful for the interpretation of the spectra and evolution of galaxies. A broad variety of empirical and theoretical data is discussed here, and the data are made fully available in the AAS CD-ROM Series, Vol. 7. Several empirical stellar libraries are part of this data base. They cover the ultraviolet spectral range observed with IUE, optical data from different ground-based telescopes, and ground-based infrared data. Spectral type coverage depends on the wavelength, but is mostly complete for types O to M and luminosity classes V to I. A large metallicity range is covered as well. Theoretical libraries of selected spectral indices of cool stars and of stellar continuum fluxes in the temperature range 2000-50,000 K as well as Wolf-Rayet energy distributions are presented. Several libraries of star clusters and early-type galaxies have been selected for this data base. We discuss an extensive set of empirical spectral templates covering the wavelength region from 1200 to 9800 Angstrom, as well as narrow-band line indices in a large number of passbands. Bench-mark spectra of nearby galaxies for model tests are included as well. We compiled numerous evolutionary models and isochrones for stars of all mass ranges of interest, wide metallicity range, and for all evolutionary phases, including the pre-main-sequence phase. The majority of the models have been computed by the Geneva and Padova groups. Evolutionary synthesis models computed by several independent groups are made available. They can be applied to old and young systems, and are optimized with respect to different aspects of input physics. The model predictions include stellar (colors, magnitudes, absorption features) and nebular (emission-line fluxes) properties. Finally, we present models of ionized gas to be used for the interpretation of active galactic nuclei and young star-forming galaxies. The community is encouraged to make use of this electronic data base and to perform a critical comparison between the individual datasets.

We present a new empirical prescription for the mass-loss rates of hydrogen-free Wolf-Rayet stars based on results of detailed spectral analyses of WC and WO stars. Compared to the prescription of Nugis & Lamers (2000), Ṁ is less sensitive to the surface helium abundance, implying a stronger mass loss at the late stages of Wolf-Rayet evolution. The winds of hydrogen-free WN stars have a strong metallicity dependence, while those of WC and WO stars have a very weak metallicity dependence.

The Tarantula Nebula (30 Dor) is a spectacular star-forming region in the Large Magellanic Cloud, seen through gas in the Galactic Disc and Halo. Diffuse Interstellar Bands offer a unique probe of the diffuse, cool-warm gas in these regions. The aim is to use DIBs as diagnostics of the local interstellar conditions, whilst at the same time deriving properties of the yet-unknown carriers. Spectra of over 800 early-type stars from the VLT Flames Tarantula Survey (VFTS) were analysed. Maps were created, separately, for the Galactic and LMC absorption in the DIBs at 4428 and 6614 Ang and - in a smaller region near the central cluster R136 - neutral sodium (Na I D); we also measured the DIBs at 5780 and 5797 Ang. The maps show strong 4428 and 6614 Ang DIBs in the quiescent cloud complex to the south of 30 Dor but weak absorption in the harsher environments to the north (bubbles) and near the OB associations. The Na maps show at least five kinematic components in the LMC and a shell-like structure surrounding R136, and small-scale structure in the Milky Way. The strengths of the 4428, 5780, 5797 and 6614 Ang DIBs are correlated, also with Na absorption and visual extinction. The strong 4428 Ang DIB is present already at low Na column density but the 6614, 5780 and 5797 Ang DIBs start to be detectable at subsequently larger Na column densities. The relative strength of the 5780 and 5797 Ang DIBs clearly confirm the Tarantula Nebula and Galactic high-latitude gas to represent a harsh radiation environment. The resilience of the 4428 Ang DIB suggests its carrier is large, compact and neutral. Structure is detected in the distribution of cool-warm gas on scales between one and >100 pc in the LMC and as little as 0.01 pc in the Sun's vicinity. Stellar winds from the central cluster R136 have created an expanding shell; some infalling gas is also detected, reminiscent of a galactic "fountain"., Accepted for publication in Astronomy and Astrophysics

We present observations of the embedded massive young stellar object (YSO) candidate 08576nr292, obtained with X-shooter and SINFONI on the ESO Very Large Telescope (VLT). The flux-calibrated, medium-resolution X-shooter spectrum (300 - 2500 nm) includes over 300 emission lines, but no (photospheric) absorption lines and is consistent with a reddened disk spectrum. Among the emission lines are three hydrogen series and helium lines, both permitted and forbidden metal lines, and CO first-overtone emission. A representative sample of lines with different morphologies is presented. The H alpha and Ca II triplet lines are very strong, with profiles indicative of outflow and - possibly - infall, usually observed in accreting stars. These lines include a blue-shifted absorption component at ~ -125 km/s. The He I and metal-line profiles are double-peaked, with a likely origin in a circumstellar disk. The forbidden lines, associated with outflow, have a single blue-shifted emission component centered at -125 km/s, coinciding with the absorption components in H alpha and Ca II. SINFONI H- and K-band integral-field spectroscopy of the cluster environment demonstrates that the [Fe II] emission is produced by a jet originating at the location of 08576nr292. Because the spectral type of the central object cannot be determined, its mass remains uncertain. We argue that 08576nr292 is an intermediate-mass YSO with a high accretion rate Mdot_acc ~ 10^-5 - 10^-6 Msun/yr. These observations demonstrate the potential of X-shooter and SINFONI to study in great detail an accretion disk-jet system, rarely seen around the more massive YSOs., Comment: 6 pages, 4 figures, 1 appendix (11 pages), accepted for publication in ApJ Letters

We have performed a pilot study of mass loss predictions for late-type Wolf-Rayet (WR) stars as a function of metal abundance, over a range between 10^{-5} < (Z/Zsun) < 10. We find that the winds of nitrogen-rich Wolf-Rayet stars are dominated by iron lines, with a dependence of mass loss on Z similar to that of massive OB stars. For more evolved, carbon-rich, WR stars the wind strength is also found to be dependent on the Fe abundance, so that they depend on the chemical environment of the host galaxy, but with a mass loss metallicity dependence that is less steep than for OB stars. Our finding that WR mass loss is Z-dependent is a new one, with important consequences for black hole formation and X-ray population studies in external galaxies. A further finding of our study is that the Z dependence of C-rich WR stars becomes weaker at metallicities below Z/Zsun < 1/10, and mass loss no longer declines once the metal abundance drops below (Z/Zsun) = 10^{-3}. This is the result of an increased importance of radiative driving by intermediate mass elements, such as carbon. In combination with rapid rotation and/or proximity to the Eddington limit -- likely to be relevant for massive Population III stars -- this effect may indicate a role for mass loss in the appearance and evolution of these objects, as well as a potential role for stellar winds in enriching the intergalactic medium of the early Universe., Accepted by Astronomy & Astrophysics (11 pages)

We discuss an interesting feature of the distribution of luminous blue variables on the H-R diagram, and we propose a connection with the bistability jump in the winds of early-type supergiants. There appears to be a deficiency of quiescent LBVs on the S Dor instability strip at luminosities between log L/Lsun = 5.6 and 5.8. The upper boundary, is also where the temperature-dependent S Dor instability strip intersects the bistability jump at about 21,000 K. Due to increased opacity, winds of early-type supergiants are slower and denser on the cool side of the bistability jump, and we postulate that this may trigger optically-thick winds that inhibit quiescent LBVs from residing there. We conduct numerical simulations of radiation-driven winds for a range of temperatures, masses, and velocity laws at log L/Lsun=5.7 to see what effect the bistability jump should have. We find that for relatively low stellar masses the increase in wind density at the bistability jump leads to the formation of a modest to strong pseudo photosphere -- enough to make an early B-type star appear as a yellow hypergiant. Thus, the proposed mechanism will be most relevant for LBVs that are post-red supergiants. Yellow hypergiants like IRC+10420 and rho Cas occupy the same luminosity range as the ``missing'' LBVs, and show apparent temperature variations at constant luminosity. If these yellow hypergiants do eventually become Wolf-Rayet stars, we speculate that they may skip the normal LBV phase, at least as far as their apparent positions on the HR diagram are concerned., Comment: 20 pages, 4 figs, accepted by ApJ

The installation of the Corrective Optics Space Telescope Axial Replacement (COSTAR) Instrument on the Hubble Space Telescope (HST) makes it possible to observe stars in very crowded regions with high spatial and spectral purity. To demonstrate this capability, we have used the Goddard High Resolution Spectrograph (GHRS) to obtain spectra of two stars in the dense center of the 30 Doradus ionizing cluster: R136a5, and its nearest neighbor, R136a2, only 0.17 sec away. R136a5 is shown to ben an O3f/WN star, while R136a2 is a WN4-w star. From both Wide Field/Planetary Camera (WFPC) photometry and GHRS, spectroscopy we estimate the following properties of R136a5: T(sub eff) = 42,500 K, R = 16.4 solar radius, L(sub bol) = 8 x 10(exp 5) solar luminosity, and M approx. equals 50 solar mass -- all indicating that, despite its spectral type, R136a5, as indicated by the strength of He II lambda 1640 emission. The observed mass-loss rate, dot-M = 1.8 x 10(exp -5) solar mass/yr, is an order of magnitude higher than is assumed by current stellar evolutionary models. We argue that this high rate of mass loss will alter drastically the evolutionary path of R136a5. If so, evolutionary models for massive stars require substantial revision.

We present a theoretical analysis of the Hubble Space Telescope (HST) Goddard High Resolution Spectrograph (GHRS) spectrum of R136a5, an O3fWN star in the R136a cluster in 30 Doradus. Using non-LTE extended and expanding model atmospheres, we find a surprisingly high mass-loss rate dot-M = 1.8 +/- 0.5 x 10(exp -5) solar mass/yr and an effective temperature T(sub eff) = 42.5 +/- 2.5 kK. With the observed visual magnitude, this implies a radius R(sub *) = 17 +/- 2 solar radius and a solar luminosity L = 8.5 +/- 1 x 10(exp 5) solar luminosity. We discuss possible sources of the uncertainties in the derived stellar parameters and conclude that the effective temperature may be underestimated; however, if so, the value of the mass-loss rate would not be affected.