Search

Massive stars are predominantly born in stellar associations or clusters. Their radiation fields, stellar winds, and supernovae strongly impact their local environment. In the first few million years of a cluster's life, massive stars are dynamically ejected running away from the cluster at high speed. However, the production rate of dynamically ejected runaways is poorly constrained. Here we report on a sample of 55 massive runaway stars ejected from the young cluster R136 in the Large Magellanic Cloud. Astrometric analysis with Gaia reveals two channels of dynamically ejected runaways. The first channel ejects massive stars in all directions and is consistent with dynamical interactions during and after the birth of R136. The second channel launches stars in a preferred direction and may be related to a cluster interaction. We find that 23-33% of the most luminous stars initially born in R136 are runaways. Model predictions have significantly underestimated the dynamical escape fraction of massive stars. Consequently, their role in shaping and heating the interstellar and galactic medium, along with their role in driving galactic outflows, is far more important than previously thought., Comment: Published in Nature on 2024 October 09. 39 pages, 17 figures

High-resolution observations reveal that the outflows of evolved low- and intermediate-mass stars harbour complex morphological structures that are linked to the presence of one or multiple companions. Hydrodynamical simulations provide a way to study the impact of a companion on the shaping of the AGB outflow. Using smoothed particle hydrodynamic (SPH) simulations of a mass losing AGB star with a binary companion, we study the impact of including HI atomic line cooling on the flow morphology and on the properties of accretion disks that form around the companion. We perform high-resolution 3D SPH simulations of the interaction of a solar-mass companion with the outflow of an AGB star using different wind velocities and eccentricities. We compare the models properties computed with and without the inclusion of HI cooling. The inclusion of HI cooling produces a large decrease in the temperature up to one order of magnitude in the region closely surrounding the companion star. As a consequence, morphological irregularities and relatively energetic outflows that were obtained without HI cooling no longer appear. In case of an eccentric orbit and low wind velocity, the morphologies are still highly asymmetric, but the same structures recur at every orbital period, making the morphology more regular. Flared accretion disks with a (sub)-Keplerian velocity profile form around the companion in all models, provided the accretion radius is small enough. The disks have radial sizes ranging from about 0.4 to 0.9 au. For the considered wind velocities, mass accretion onto the companion is up to 2 times higher than predicted by the BHL rate, and ranges between 0.04-0.21 times the AGB wind mass loss rate. The lower the wind velocity at the location of the companion, the larger and the more massive the disk, and the higher the mass accretion efficiency., Comment: 18 pages (32 including appendices), 18 figures (33 including apendices) For associated input files, final output dumps and 3 movies, see https://zenodo.org/records/13224218

The interplay between T Tauri stars and their circumstellar disks, and how this impacts the onset of planet formation has yet to be established. We studied a seemingly old T Tauri star, PDS 111, and its disk. We analyzed optical, infrared, and sub-millimeter observations obtained with VLT/X-shooter, Mercator/HERMES, TESS, VLT/SPHERE, and ALMA, providing a new view on PDS 111 and its protoplanetary disk. The multi-epoch spectroscopy yields photospheric lines to classify the star, and emission lines to study variability in the hot inner disk and to determine the mass-accretion rate. The SPHERE and ALMA observations are used to characterize the dust distribution of the small and large grains, respectively. PDS 111 is a weak-line T Tauri star with spectral type G2, exhibits strong H$\alpha$ variability and with a low mass-accretion rate of $1-5\times10^{-10}$\,M$_{\odot}$\,yr$^{-1}$. We measured an age of the system of 15.9$^{+1.7}_{-3.7}$ Myr using pre-main sequence tracks. The SPHERE observations show a strongly flaring disk with an asymmetric substructure. The ALMA observations reveal a 30 au cavity in the dust continuum emission with a low contrast asymmetry in the South-West of the disk and a dust disk mass of 45.8\,$M_\oplus$. The $^{12}$CO radial extension is at least three times larger than that of the dust emission. Although the measured age is younger than suggested in literature, PDS 111 still seems relatively old; this provides insight into disk properties at an advanced stage of pre-main sequence evolution. The characteristics of this disk are very similar to its younger counterparts: strongly flaring, an average disk mass, a typical radial extent of the disk gas and dust, and the presence of common substructures. This suggests that disk evolution has not significantly changed the disk properties. These results show similarities with the "Peter Pan disks" around M-dwarfs., Comment: 23 pages, 23 figures, accepted by A&A (abstract shortened)

In this paper, we explore the significant, non-linear impact that stellar winds have on H ii regions. We perform a parameter study using three-dimensional radiative magnetohydrodynamic simulations of wind and ultraviolet radiation feedback from a 35 Msun star formed self-consistently in a turbulent, self-gravitating cloud, similar to the Orion Nebula (M42) and its main ionizing source Theta 1 Ori C. Stellar winds suppress early radiative feedback by trapping ionizing radiation in the shell around the wind bubble. Rapid breakouts of warm photoionized gas ('champagne flows') still occur if the star forms close to the edge of the cloud. The impact of wind bubbles can be enhanced if we detect and remove numerical overcooling caused by shocks crossing grid cells. However, the majority of the energy in the wind bubble is still lost to turbulent mixing between the wind bubble and the gas around it. These results begin to converge if the spatial resolution at the wind bubble interface is increased by refining the grid on pressure gradients. Wind bubbles form a thin chimney close to the star, which then expands outwards as an extended plume once the wind bubble breaks out of the dense core the star formed in, allowing them to expand faster than a spherical wind bubble. We also find wind bubbles mixing completely with the photoionized gas when the H ii region breaks out of the cloud as a champagne flow, a process we term 'hot champagne'., Comment: 19 pages, 17 figures, published in MNRAS. Sorry it took so long to upload here, I got distracted

The fraction of massive stars in young stellar clusters is of importance as they are the dominant sources of both mechanical and radiative feedback, strongly influencing the thermal and dynamical state of their birth environments. It turns out that a significant fraction of massive stars escape from their parent cluster via dynamical interactions of single stars and/or multiple stellar systems. M 17 is the nearest giant H II region hosting a very young and massive cluster: NGC 6618. Our aim is to identify stars brighter than G < 21 mag that belong to NGC 6618, including the (massive) stars that may have escaped since its formation, and to determine the cluster distance and age. We identified 42 members of NGC 6618 of which eight have a spectral type of O, with a mean distance of 1675 pc and a transversal velocity dispersion of about 3 km/s , and a radial velocity dispersion of 6 km/s. Another ten O stars are associated with NGC 6618, but they cannot be classified as members due to poor astrometry or high extinction. We have also identified six O star runaways. The relative transverse velocity of these runaways ranges from 10 to 70 km/s and their kinematic age ranges from about 100 to 750 kyr. Given the already established young age of NGC 6618 (< 1 Myr), this implies that massive stars are being ejected from the cluster already directly after or during the cluster formation process. When constructing the initial mass function, one has to take into account the massive stars that have already escaped from the cluster, that is, about 30% of the O stars of the original population of NGC 6618. The trajectories of the O runaways can be traced back to the central 0.25 pc region of NGC 6618. The good agreement between the evolutionary and kinematic age of the runaways implies that the latter provides an independent way to estimate (a lower limit to) the age of the cluster., Comment: to be published in A&A

We present the first results of the eXtreme UV Environments (XUE) James Webb Space Telescope (JWST) program, that focuses on the characterization of planet forming disks in massive star forming regions. These regions are likely representative of the environment in which most planetary systems formed. Understanding the impact of environment on planet formation is critical in order to gain insights into the diversity of the observed exoplanet populations. XUE targets 15 disks in three areas of NGC 6357, which hosts numerous massive OB stars, among which some of the most massive stars in our Galaxy. Thanks to JWST we can, for the first time, study the effect of external irradiation on the inner ($< 10$ au), terrestrial-planet forming regions of proto-planetary disks. In this study, we report on the detection of abundant water, CO, CO$_2$, HCN and C$_2$H$_2$ in the inner few au of XUE 1, a highly irradiated disk in NGC 6357. In addition, small, partially crystalline silicate dust is present at the disk surface. The derived column densities, the oxygen-dominated gas-phase chemistry, and the presence of silicate dust are surprisingly similar to those found in inner disks located in nearby, relatively isolated low-mass star-forming regions. Our findings imply that the inner regions of highly irradiated disks can retain similar physical and chemical conditions as disks in low-mass star-forming regions, thus broadening the range of environments with similar conditions for inner disk rocky planet formation to the most extreme star-forming regions in our Galaxy., Comment: Accepted for publication in ApJ Letters. 20 pages, 7 figures

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., Comment: Accepted for publication in A&A (18 pages, 12 figures)

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 April 2022, and identify key areas where further dialogue can bring about radical changes in how we view the relationship between stars and the universe they live in., Comment: Accepted to the Publications of the Astronomical Society of the Pacific

In the first few Myr the massive stars dynamically interact, produce runaways and affect the initial binary population. Observing and interpreting the dynamics of young massive clusters is key to our understanding of the star formation process and predicting the outcome of stellar evolution. We have studied NGC6611 in the Eagle Nebula (M16), a young massive cluster hosting 19 O stars. We used Gaia EDR3 data to determine the membership, age, cluster dynamics and the kinematics of the massive stars including runaways. The membership analysis yields 137 members located at a mean distance of 1706 $\pm$ 7 pc. The colour - absolute magnitude diagram reveals a blue and a red population of pre-main-sequence stars, consistent with two distinct populations of stars. In line with earlier studies, the youngest population has a mean extinction $A_V$ = 3.6 $\pm$ 0.1 mag and an age = 1.3 $\pm$ 0.2 Myr, while the older population of stars has a mean extinction $A_V$ = 2.0 $\pm$ 0.1 mag and an age = 7.5 $\pm$ 0.4 Myr. The latter population is more spatially extended than the younger generation of stars. We argue that most of the OB stars belong to the younger population. We identify 8 runaways originating from the center of NGC6611, consistent with the dynamical ejection scenario. We show that ~ 50% of the O stars have velocities comparable to or greater than the escape velocity. These O stars can be traced back to the center of NGC6611 with kinematic ages ranging from 0 to 2 Myr. This suggests that dynamical interactions played an important role in the early evolution of NGC6611, which is surprising considering the low current stellar density. Comparing this to simulations of young massive clusters, the required initial radius of 0.1-0.5 pc is not consistent with that of NGC6611. The O stars could have initially formed in wide binaries and possibly harden through dynamical interactions., Comment: replaced with accepted version to A&A. 28 pages, 15 figures

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 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-sigma 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., Comment: Accepted to Nature Astronomy, 64 pages, 15 figures, 2 tables; ESO press release: https://www.eso.org/public/news/eso2210/; Nat Asr paper URL: https://www.nature.com/articles/s41550-022-01730-y

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

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

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

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., Comment: 16 pages, A&A in press

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., Comment: 21 pages, 23 figures

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 $<$2 days and total masses $\lessapprox45\,\textrm{M}_{\odot}$. We expect stellar winds, non-conservative mass transfer and envelope inflation to have played a role in the formation of the more massive and longer-period contact binaries., Comment: 23 pages, 11 figures, accepted for publication in MNRAS

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 (<10 au) companion, was observed with the Atacama Large Millimeter/submillimeter Array as part of the ATOMIUM Large programme. In this work, the brightest CO, SiO, and HCN molecular line transitions are analysed. The continuum map shows two maxima, separated by 0.04'' (6 au). The CO data unambiguously reveal that Pi1 Gru's circumstellar environment harbours an inclined, radially outflowing, equatorial density enhancement. It contains a spiral structure at an angle of 38+/-3 deg with the line-of-sight. The HCN emission in the inner wind reveals a clockwise spiral, with a dynamical crossing time of the spiral arms consistent with a companion at a distance of 0.04'' from the asymptotic giant branch star, which is in agreement with the position of the secondary continuum peak. The inner wind dynamics imply a large acceleration region, consistent with a beta-law power of ~6. The CO emission suggests that the spiral is approximately Archimedean within 5'', beyond which this trend breaks down as the succession of the spiral arms becomes less periodic. The SiO emission at scales smaller than 0.5'' exhibits signatures of gas in rotation, which is found to fit the expected behaviour of gas in the wind-companion interaction zone. An investigation of SiO maser emission reveals what could be a stream of gas accelerating from the surface of the AGB star to the companion. Using these dynamics, we have tentatively derived an upper limit on the companion mass to be ~1.1 Msol., Comment: 19 pages, 20 figures

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

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., Comment: Published in Astronomy & Astrophysics, Volume 645

Rapid rotation in massive stars imposes a latitudinal variation in the mass loss from radiatively driven winds that can lead to enhanced mass loss at the poles (with little angular momentum loss) and/or equator (with maximal angular momentum loss). Here we present an examination of the stellar wind lines of the two O-type stars with the fastest known equatorial velocities, VFTS 102 ($V\sin i = 610 \pm 30$ km/s; O9:Vnnne+) and VFTS 285 ($V\sin i = 609 \pm 29$ km/s; O7.5 Vnnn) in the Large Magellanic Cloud. Ultraviolet spectra of both stars were obtained with the Hubble Space Telescope Cosmic Origins Spectrograph. The spectrum of VFTS 285 displays a fast outflow in N V and a much slower wind in Si IV, and we argue that there is a two-wind regime in which mass loss is strong at the poles (fast and tenuous wind) but dominant at the equator (slow and dense winds). These ions and wind lines are not present in the spectrum of the cooler star VFTS 102, but the double-peaked H$\alpha$ emission in its spectrum implies equatorial mass loss into a circumstellar disk. The results suggest that in the fastest rotating O-stars, most mass is lost as an equatorial outflow, promoting angular momentum loss that contributes to a spin down over time., Comment: Accepted for publication in ApJ. 17 pages, 5 figures

Unravelling the composition and characteristics of gas and dust lost by asymptotic giant branch (AGB) stars is important as these stars play a vital role in the chemical life cycle of galaxies. The general hypothesis of their mass loss mechanism is a combination of stellar pulsations and radiative pressure on dust grains. However, current models simplify dust formation, which starts as a microscopic phase transition called nucleation. Various nucleation theories exist, yet all assume chemical equilibrium, growth restricted by monomers, and commonly use macroscopic properties for a microscopic process. Such simplifications for initial dust formation can have large repercussions on the type, amount, and formation time of dust. By abandoning equilibrium assumptions, discarding growth restrictions, and using quantum mechanical properties, we have constructed and investigated an improved nucleation theory in AGB wind conditions for four dust candidates, TiO$_2$, MgO, SiO and Al$_2$O$_3$. This paper reports the viability of these candidates as first dust precursors and reveals implications of simplified nucleation theories. Monomer restricted growth underpredicts large clusters at low temperatures and overpredicts formation times. Assuming the candidates are present, Al$_2$O$_3$ is the favoured precursor due to its rapid growth at the highest considered temperatures. However, when considering an initially atomic chemical mixture, only TiO$_2$-clusters form. Still, we believe Al$_2$O$_3$ to be the prime candidate due to substantial physical evidence in presolar grains, observations of dust around AGB stars at high temperatures, and its ability to form at high temperatures and expect the missing link to be insufficient quantitative data of Al-reactions., Comment: Accepted for publication in MNRAS. 19 pages (68 incl. appendix)

Current burning issues in stellar physics, for both hot and cool stars, concern their magnetism. In hot stars, stable magnetic fields of fossil origin impact their stellar structure and circumstellar environment, with a likely major role in stellar evolution. However, this role is complex and thus poorly understood as of today. It needs to be quantified with high-resolution UV spectropolarimetric measurements. In cool stars, UV spectropolarimetry would provide access to the structure and magnetic field of the very dynamic upper stellar atmosphere, providing key data for new progress to be made on the role of magnetic fields in heating the upper atmospheres, launching stellar winds, and more generally in the interaction of cool stars with their environment (circumstellar disk, planets) along their whole evolution. UV spectropolarimetry is proposed on missions of various sizes and scopes, from POLLUX on the 15-m telescope LUVOIR to the Arago M-size mission dedicated to UV spectropolarimetry., Comment: 25 pages, 11 figures, white paper submitted to ESA Voyage 2050 call

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., Comment: The paper has been accepted for publication in Astronomy and Astrophysics

We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena., Comment: 60 pages, 22 figures, minor updates

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

The material lost through stellar winds of Asymptotic Giant Branch (AGB) stars is one of the main contributors to the chemical enrichment of galaxies. The general hypothesis of the mass loss mechanism of AGB winds is a combination of stellar pulsations and radiative pressure on dust grains, yet current models still suffer from limitations. Among others, they assume chemical equilibrium of the gas, which may not be justified due to rapid local dynamical changes in the wind. This is important as it is the chemical composition that regulates the thermal structure of the wind, the creation of dust grains in the wind, and ultimately the mass loss by the wind. Using a self-consistent hydrochemical model, we investigated how non-equilibrium chemistry affects the dynamics of the wind. This paper compares a hydrodynamical and a hydrochemical dust-free wind, with focus on the chemical heating and cooling processes. No sustainable wind arises in a purely hydrodynamical model with physically reasonable pulsations. Moreover, temperatures are too high for dust formation to happen, rendering radiative pressure on grains impossible. A hydrochemical wind is even harder to initiate due to efficient chemical cooling. However, temperatures are sufficiently low in dense regions for dust formation to take place. These regions occur close to the star, which is needed for radiation pressure on dust to sufficiently aid in creating a wind. Extending this model self-consistently with dust formation and evolution, and including radiation pressure, will help to understand the mass loss by AGB winds., Comment: Accepted for publication in MNRAS. 30 pages (incl. Appendix), 19 figures

CubeSats are routinely used for low-cost photometry from space. Space-borne spectroscopy, however, is still the exclusive domain of much larger platforms. Key astrophysical questions in e.g. stellar physics and exoplanet research require uninterrupted spectral monitoring from space over weeks or months. Such monitoring of individual sources is unfortunately not affordable with these large platforms. With CUBESPEC we plan to offer the astronomical community a low-cost CubeSat solution for near-UV/optical/near-IR spectroscopy that enables this type of observations. CUBESPEC is a generic spectrograph that can be configured with minimal hardware changes to deliver both low resolution (R=100) with very large spectral coverage (200-1000nm), as well as high resolution (R=30,000) over a selected wavelength range. It is built around an off-axis Cassegrain telescope and a slit spectrograph with configurable dispersion elements. CUBESPEC will use a compact attitude determination and control system for coarse pointing of the entire spacecraft, supplemented with a fine-guidance system using a fast steering mirror to center the source on the spectrograph slit and to cancel out satellite jitter. An extremely compact optical design allows us to house this instrument in a 6U CubeSat with a volume of only 10x20x30cm$^{3}$, while preserving a maximized entrance pupil of ca. 9x19cm$^{2}$. In this contribution, we give an overview of the CUBESPEC project, discuss its most relevant science cases, and present the design of the instrument., Comment: SPIE Astronomical Telescopes + Instrumentation conference proceedings, Austin 2018

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..., Comment: 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

Equatorial density enhancements (EDEs) are a very common astronomical phenomenon. Studies of the circumstellar environments (CSE) of young stellar objects and of evolved stars have shown that these objects often possess these features. These are believed to originate from different mechanisms, ranging from binary interactions to the gravitational collapse of interstellar material. Quantifying the effect of the presence of this type of EDE on the observables is essential for a correct interpretation of high-resolution data. We seek to investigate the manifestation in the observables of a circumstellar EDE, to assess which properties can be constrained, and to provide an intuitive bedrock on which to compare and interpret upcoming high-resolution data (e.g. \emph{ALMA} data) using 3D models. We develop a simplified analytical parametrised description of a 3D EDE, with a possible substructure such as warps, gaps, and spiral instabilities. In addition, different velocity fields (Keplerian, radial, super-Keplerian, sub-Keplerian and rigid rotation) are considered. The effect of a bipolar outflow is also investigated. The geometrical models are fed into the 3D radiative transfer code {\tt LIME}, that produces 3D intensity maps throughout velocity space. We investigate the spectral signature of the $J$=3$-$2 up to $J$=7$-$6 rotational transitions of CO in the models, as well as the spatial aspect of this emission by means of channel maps, wide-slit position-velocity (PV) diagrams, stereograms, and spectral lines. Additionally, we discuss methods of constraining the geometry of the EDE, the inclination, the mass-contrast between the EDE and the bipolar outflow, and the global velocity field. Finally, we simulated \emph{ALMA} observations...

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., Comment: 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., Comment: 29 pages (including 14 pages appendix), 15 figures, accepted by A&A

The formation process of massive stars is still poorly understood. Massive young stellar objects (mYSOs) are deeply embedded in their parental clouds, they are rare and thus typically distant, and their reddened spectra usually preclude the determination of their photospheric parameters. M17 is one of the best studied HII regions in the sky, is relatively nearby, and hosts a young stellar population. With X-shooter on the ESO Very Large Telescope we have obtained optical to near-infrared spectra of candidate mYSOs, identified by Hanson et al. (1997), and a few OB stars in this region. The large wavelength coverage enables a detailed spectroscopic analysis of their photospheres and circumstellar disks. We confirm the pre-main sequence (PMS) nature of six of the stars and characterise the O stars. The PMS stars have radii consistent with being contracting towards the main sequence and are surrounded by a remnant accretion disk. The observed infrared excess and the (double-peaked) emission lines provide the opportunity to measure structured velocity profiles in the disks. We compare the observed properties of this unique sample of young massive stars with evolutionary tracks of massive protostars by Hosokawa & Omukai (2009), and propose that these mYSOs near the western edge of the HII region are on their way to become main-sequence stars ($\sim 6 - 20$ $M_{\odot}$) after having undergone high mass-accretion rates (${\dot{M}_{\rm acc}} \sim 10^{-4} - 10^{-3}$ $M_{\odot}$ $\rm yr^{-1}$). Their spin distribution upon arrival at the zero age main sequence (ZAMS) is consistent with that observed for young B stars, assuming conservation of angular momentum and homologous contraction., Comment: Accepted for publication in A&A. Appendixes A and B have been truncated due to size limitations, the full version will be available on A&A

Recent high-resolution observations have shown stellar winds to harbour complexities which strongly deviate from spherical symmetry, generally assumed as standard wind model. One such morphology is the archimedean spiral, generally believed to be formed by binary interactions, which has been directly observed in multiple sources. 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. By means of an extended parameter study, we model rotational CO emission from the stellar outflow of asymptotic giant branch stars. To this end, we develop 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 investigate the spectral signature of rotational transitions of CO of the models, as well as the spatial aspect of this emission by means of wide-slit PV diagrams. Additionally, the potential for misinterpretation of the 3D data in a 1D context is quantified. Finally, we simulate ALMA observations to explore the impact of interefrometric noise and artifacts on the emission signatures. The spectral signatures of the CO rotational transition v=0 J=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. The inability to disentangle the spiral from the spherical signal can result in an incorrect interpretation in a 1D context. Consequently, erroneous mass loss rates would be calculated...

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible., Comment: Significantly expanded and updated version of previous ELT-MOS White Paper, so there is some textual overlap with arXiv:1303.0029

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

Massive stars rapidly change their masses through strong stellar winds and mass transfer in binary systems. We show that such mass changes leave characteristic signatures in stellar mass functions of young star clusters which can be used to infer their ages and to identify products of binary evolution. We model the observed present day mass functions of the young Galactic Arches and Quintuplet star clusters using our rapid binary evolution code. We find that shaping of the mass function by stellar wind mass loss allows us to determine the cluster ages to 3.5$\pm$0.7 Myr and 4.8$\pm$1.1 Myr, respectively. Exploiting the effects of binary mass exchange on the cluster mass function, we find that the most massive stars in both clusters are rejuvenated products of binary mass transfer, i.e. the massive counterpart of classical blue straggler stars. This resolves the problem of an apparent age spread among the most luminous stars exceeding the expected duration of star formation in these clusters. We perform Monte Carlo simulations to probe stochastic sampling, which support the idea of the most massive stars being rejuvenated binary products. We find that the most massive star is expected to be a binary product after 1.0$\pm$0.7 Myr in Arches and after 1.7$\pm$1.0 Myr in Quintuplet. Today, the most massive 9$\pm$3 stars in Arches and 8$\pm$3 in Quintuplet are expected to be such objects. Our findings have strong implications for the stellar upper mass limit and solve the discrepancy between the claimed 150 $\mathrm{M}_\odot$ limit and observations of fours stars with initial masses of 165-320 $\mathrm{M}_\odot$ in R136 and of SN 2007bi, which is thought to be a pair-instability supernova from an initial 250 $\mathrm{M}_\odot$ star. Using the stellar population of R136, we revise the upper mass limit to values in the range 200-500 $\mathrm{M}_\odot$., Comment: 13 pages (+4 pages appendix), 8 figures, 2 tables; accepted for publication in ApJ

Globular clusters contain many stars with surface abundance patterns indicating contributions from hydrogen burning products, as seen in the anti-correlated elemental abundances of e.g. sodium and oxygen, and magnesium and aluminium. Multiple generations of stars can explain this phenomenon, with the second generation forming from a mixture of pristine gas and ejecta from the first generation. We show that massive binary stars may be a source of much of the material that makes this second generation of stars. Mass transfer in binaries is often non-conservative and the ejected matter moves slowly enough that it can remain inside a globular cluster and remain available for subsequent star formation. Recent studies show that there are more short-period massive binaries than previously thought, hence also more stars that interact and eject nuclear-processed material., Comment: Contribution to the proceedings of "Reading the book of globular clusters with the lens of stellar evolution", Rome, 26-28 November 2012

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"., Comment: Accepted for publication in Astronomy and Astrophysics

As part of the VLT/X-shooter science verification, we obtained the first optical medium-resolution spectrum of a previously identified bright O-type object in NGC55, an LMC-like galaxy at a distance of \sim2.0 Mpc. Based on the stellar and nebular spectrum, we investigate the nature and evolutionary status of the central object(s) and its influence on the surrounding interstellar medium. We conclude that the source, NGC55_C1_31, is a composite object, likely a stellar cluster, which contains one or several hot (T_eff \simeq 50000 K) WN stars with a high mass-loss rate (\sim3 \times 10^{-5} M_\odot yr^{-1}) and a helium-rich composition (N_He/N_H = 0.8). The visual flux is dominated by OB-type (super)giant stars with T_eff \sim< 35000 K, solar helium abundance (N_He/N_H = 0.1), and mass-loss rate \sim2 \times 10^{-6} M_\odot yr^{-1}. The surrounding H II region has an electron density n_e < 10^2 cm^{-3} and an electron temperature T(OIII) \simeq 11500 \pm 600 K. The oxygen abundance of this region is [O/H] = 8.18 \pm 0.03 which corresponds to Z = 0.31 \pm 0.04 Z_\odot. We observed no significant gradients in T(OIII), n_e or [O/H] on a scale of 73 pc extending in four directions from the ionising source. The properties of the HII region can be reproduced by a CLOUDY model which uses the central cluster as ionising source, thus providing a self-consistent interpretation of the data. We also report on the serendipitous discovery of HeII nebular emission associated with the nearby source NGC55_C2_35, a feature usually associated with strong X-ray sources., Comment: 12 pages, 10 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society; the definitive version will be available at wwww.blackwell-synergy.com

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

Rotational mixing in massive stars is a widely applied concept, with far reaching consequences for stellar evolution. Nitrogen surface abundances for a large and homogeneous sample of massive B-type stars in the LMC were obtained by the VLT-FLAMES Survey of Massive Stars. This sample is the first covering a broad range of projected stellar rotational velocities, with a large enough sample of high quality data to allow for a statistically significant analysis. We use the sample to provide the first rigorous test of the theory of rotational mixing in massive stars. We calculated a grid of stellar evolution models, using the FLAMES sample to calibrate some of the uncertain mixing processes. We developed a new population-synthesis code, which uses this grid to simulate a large population of stars with masses, ages and rotational velocity distributions consistent with those from the FLAMES sample. The synthesized population is then filtered by the selection effects in the observed sample, to enable a direct comparison between the empirical results and theoretical predictions. Our simulations reproduce the fraction of stars without significant nitrogen enrichment. The predicted number of rapid rotators with enhanced nitrogen is about twice as large as found observationally. Furthermore, a group of stars consisting of slowly rotating, nitrogen-enriched objects and another consisting of rapidly rotating un-enriched objects can not be reproduced by our single-star population synthesis. Additional physical processes appear to be required to understand the population of massive main-sequence stars from the FLAMES sample.We discuss the possible role of binary stars and magnetic fields in the interpretation of our results. We find that the population of slowly rotating nitrogen-enriched stars is unlikely produced via mass transfer and subsequent tidal spin-down in close binary systems, Comment: Accepted for publication in A&A

We present a dense grid of evolutionary tracks and isochrones of rotating massive main-sequence stars. We provide three grids with different initial compositions tailored to compare with early OB stars in the Small and Large Magellanic Clouds and in the Galaxy. Each grid covers masses ranging from 5 to 60 Msun and initial rotation rates between 0 and about 600 km/s. To calibrate our models we used the results of the VLT-FLAMES Survey of Massive Stars. We determine the amount of convective overshooting by using the observed drop in rotation rates for stars with surface gravities log g < 3.2 to determine the width of the main sequence. We calibrate the efficiency of rotationally induced mixing using the nitrogen abundance determinations for B stars in the Large Magellanic cloud. We describe and provide evolutionary tracks and the evolution of the central and surface abundances. In particular, we discuss the occurrence of quasi-chemically homogeneous evolution, i.e. the severe effects of efficient mixing of the stellar interior found for the most massive fast rotators. We provide a detailed set of isochrones for rotating stars. Rotation as an initial parameter leads to a degeneracy between the age and the mass of massive main sequence stars if determined from its observed location in the Hertzsprung-Russell diagram. We show that the consideration of surface abundances can resolve this degeneracy., Comment: Accepted for publication in A&A

We report the results of observations of ten rotational transitions of water vapor toward the carbon-rich AGB (asymptotic giant branch) star IRC+10216 (CW Leonis), carried out with Herschel's HIFI instrument. Each transition was securely detected by means of observations using the dual beam switch mode of HIFI. The measured line ratios imply that water vapor is present in the inner outflow at small distances (few x 1.E+14 cm) from the star, confirming recent results reported by Decin et al. from observations with Herschel's PACS and SPIRE instruments. This finding definitively rules out the hypothesis that the observed water results from the vaporization of small icy objects in circular orbits. The origin of water within the dense C-rich envelope of IRC+10216 remains poorly understood. We derive upper limits on the H2-17O/H2-16O and H2-18O/H2-16O isotopic abundance ratios of ~ 5.E-3 (3 sigma), providing additional constraints on models for the origin of the water vapor in IRC+10216., Comment: Accepted for publication in the Astrophysical Journal Letters

Massive stars and supernovae (SNe) have a huge impact on their environment. Despite their importance, a comprehensive knowledge of which massive stars produce which SNe is hitherto lacking. We use a Monte Carlo method to predict the mass-loss rates of massive stars in the Hertzsprung-Russell Diagram (HRD) covering all phases from the OB main sequence, the unstable Luminous Blue Variable (LBV) stage, to the final Wolf-Rayet (WR) phase. Although WR produce their own metals, a strong dependence of the mass-loss rate on the initial iron abundance is found at sub-solar metallicities (1/10 -- 1/100 solar). This may present a viable mechanism to prevent the loss of angular momentum by stellar winds, which could inhibit GRBs occurring at solar metallicities -- providing a significant boost to the collapsar model. Furthermore, we discuss recently reported quasi-sinusoidal modulations in the radio lightcurves of SNe 2001ig and 2003bg. We show that both the sinusoidal behaviour and the recurrence timescale of these modulations are consistent with the predicted mass-loss behaviour of LBVs. We discuss potential ramifications for the ``Conti'' scenario for massive star evolution., Comment: 8 pages, 6 figures, conference talk. To be published in the proceedings of "Mass loss from stars and the evolution of stellar clusters"

We report on observations of circumstellar disks around young stars that have been obtained with the MIDI instrument, which is mounted on the VLT Interferometer and operates in the 10 micrometer atmospheric window. The maximum spatial resolution of 5 milli-arcsec corresponds to sub-AU scales at the distance to nearby star formation regions. Thus, we can study the disks on the spatial scales at which important processes occur, such as accretion, dust processing, and planet formation. The main results obtained so far can be summarized as follows: 1. The measured interferometric visibilities are in good qualitative agreement with those predicted by models of circumstellar disks. In particular, a predicted correlation between the strength of the far-infrared excess and the spatial structure of the disk is confirmed by direct measurements; 2. In several objects strong evidence for deviations from circular symmetry is present, indicating that an inclined disk is indeed the dominant component seen in the mid-infrared; 3. The dust properties are not uniform over the disk, but are instead a strong function of distance to the central star. The dust in the innermost disk regions is observed to be more ``processed'' than the dust further out, both in Herbig Ae star disks and in those around T-Tauri stars., Comment: 19 pages, 12 figures, invited talk at SPIE conference "Astronomical telescopes and instrumentation", Orlando, may 2006