Search

The first stars might have been fast rotators. This would have important consequences for their radiative, mechanical and chemical feedback. We discuss the impact of fast initial rotation on the evolution of massive Population III models and on their nitrogen and oxygen stellar yields. We explore the evolution of Population III stars with initial masses in the range of 9Msol < Mini < 120Msol starting with an initial rotation on the Zero Age Main Sequence equal to 70% of the critical one. We find that with the physics of rotation considered here, our rapidly-rotating Population III stellar models do not follow a homogeneous evolution. They lose very little mass in case mechanical winds are switched on when the surface rotation becomes equal or larger than the critical velocity. Impact on the ionising flux appears modest when compared to moderately-rotating models. Fast rotation favours, in models with initial masses above ~20Msol, the appearance of a very extended intermediate convective zone around the H-burning shell during the core He-burning phase. This shell has important consequences on the sizes of the He- and CO-cores and thus impacts the final fate of stars. Moreover, it has a strong impact on nucleosynthesis boosting the production of primary 14N. Fast initial rotation impacts significantly the chemical feedback of Population III stars. Observations of extremely metal-poor stars and/or starbursting regions are essential to provide constraints on the properties of the first stars., Comment: Accepted in A&A. Pages 17. Figures 15. Tables 2

A growing number of supernovae (SNe) are now known to exhibit evidence for significant interaction with a dense, pre-existing, circumstellar medium (CSM). SNe Ibn comprise one such class that can be characterised by both rapidly evolving light curves and persistent narrow He I lines. The origin of such a dense CSM in these systems remains a pressing question, specifically concerning the progenitor system and mass-loss mechanism. In this paper, we present multi-wavelength data of the Type Ibn SN 2020nxt, including $HST$/STIS ultraviolet spectra. We fit the data with recently updated CMFGEN models designed to handle configurations for SNe Ibn. The UV coverage yields strong constraints on the energetics and, when combined with the CMFGEN models, offer new insight on potential progenitor systems. We find the most successful model is a $\lesssim4 {\rm M}_\odot$ helium star that lost its $\sim 1\,{\rm M}_\odot$ He-rich envelope in the years preceding core collapse. We also consider viable alternatives, such as a He white dwarf merger. Ultimately, we conclude at least some SNe Ibn do not arise from single, massive ($>30 {\rm M}_\odot$) Wolf-Rayet-like stars., Comment: 17 pages, 13 figures, 1 table, submitted to MNRAS

Context: Grids of stellar models, computed with the same physical ingredients, allow one to study the impact of a given physics on a broad range of initial conditions and are a key ingredient for modeling the evolution of galaxies. Aims: We present a grid of single star models for masses between 0.8 and 120 $M_\odot$, with and without rotation for a mass fraction of heavy element Z=0.006, representative of the Large Magellanic Cloud (LMC). Methods: We used the Geneva stellar evolution code. The evolution was computed until the end of the central carbon-burning phase, the early asymptotic giant branch phase, or the core helium-flash for massive, intermediate, and low mass stars, respectively. Results: The outputs of the present stellar models are well framed by the outputs of the two grids obtained by our group for metallicities above and below the one considered here. The models of the present work provide a good fit to the nitrogen surface enrichments observed during the main sequence for stars in the LMC with initial masses around 15 $M_\odot$. They also reproduce the slope of the luminosity function of red supergiants of the LMC well, which is a feature that is sensitive to the time-averaged mass loss rate over the red supergiant phase. The most massive black hole that can be formed from the present models at Z=0.006 is around 55 $M_\odot$. No model in the range of mass considered will enter into the pair-instability supernova regime, while the minimal mass to enter the region of pair pulsation instability is around 60 $M_\odot$ for the rotating models and 85 $M_\odot$ for the nonrotating ones. Conclusions: The present models are of particular interest for comparisons with observations in the LMC and also in the outer regions of the Milky Way. We provide public access to numerical tables that can be used for computing interpolated tracks and for population synthesis studies., Comment: 11 pages, 11 figures, published in A&A

We present a grid of stellar models at super-solar metallicity (Z = 0.020) extending the previous grids of Geneva models at solar and sub-solar metallicities. A metallicity of Z = 0.020 was chosen to match that of the inner Galactic disk. A modest increase of 43% (=0.02/0.014) in metallicity compared to solar models means that the models evolve similarly to solar models but with slightly larger mass loss. Mass loss limits the final total masses of the super-solar models to 35 M$_\odot$ even for stars with initial masses much larger than 100 M$_\odot$. Mass loss is strong enough in stars above 20 M$_\odot$ for rotating stars (25 M$_\odot$ for non-rotating stars) to remove the entire hydrogen-rich envelope. Our models thus predict SNII below 20 M$_\odot$ for rotating stars (25 M$_\odot$ for non-rotating stars) and SNIb (possibly SNIc) above that. We computed both isochrones and synthetic clusters to compare our super-solar models to the Westerlund 1 (Wd1) massive young cluster. A synthetic cluster combining rotating and non-rotating models with an age spread between log10 (age/yr) = 6.7 and 7.0 is able to reproduce qualitatively the observed populations of WR, RSG and YSG stars in Wd1, in particular their simultaneous presence at log10(L/L$_\odot$) = 5-5.5. The quantitative agreement is imperfect and we discuss the likely causes: synthetic cluster parameters, binary interactions, mass loss and their related uncertainties. In particular, mass loss in the cool part of the HRD plays a key role., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS

The evolution of massive stars is affected by a variety of physical processes including convection, rotation, mass loss and binary interaction. Because these processes modify the internal chemical abundance profiles in multiple ways simultaneously, it can be challenging to determine which properties of the stellar interior are primarily driving the overall evolution. Building on previous work, we develop a new modelling approach called SNAPSHOT that allows us to isolate the key features of the internal abundance profile that drive the evolution of massive stars. Using our approach, we compute numerical stellar structure models in thermal equilibrium covering key phases of stellar evolution. For the main sequence, we demonstrate that models with the same mass and very similar surface properties can have different internal distributions of hydrogen and convective core masses. We discuss why massive stars expand after the main sequence and the fundamental reasons for why they become red, blue or yellow supergiants. For the post-main sequence, we demonstrate that small changes in the abundance profile can cause very large effects on the surface properties. We also discuss the effects that produce blue supergiants and the cause of blue loops. Our models show that massive stars with lower metallicity tend to be more compact due to the combined effect of lower CNO abundances in the burning regions and lower opacity in the envelope., Comment: Accepted in MNRAS

We investigate the low-luminosity supernova SN 2016bkv and its peculiar early-time interaction. For that, we compute radiative transfer models using the CMFGEN code. Because SN 2016bkv shows signs of interaction with material expelled by its progenitor, it offers a great opportunity to constrain the uncertain evolutionary channels leading to low-luminosity supernovae. Our models indicate that the progenitor had a mass-loss rate of (6.0 +- 2.0) x 1e-4 Msun/yr (assuming a velocity of 150 km/s). The surface abundances of the progenitor are consistent with solar contents of He and CNO. If SN 2016bkv's progenitor evolved as a single star, it was an odd red supergiant that did not undergo the expected dredge up for some reason. We propose that the progenitor more likely evolved through binary interaction. One possibility is that the primary star accreted unprocessed material from a companion and avoided further rotational and convective mixing until the SN explosion. Another possibility is a merger with a lower mass star, with the primary remaining with low N abundance until core collapse. Given the available merger models, we can only put a loose constraint on the pre-explosion mass around 10-20 Msun, with lower values being favored based on previous observational constraints from the nebular phase., Comment: 10 pages, 8 figures, MNRAS accepted

The first stars are thought to be one of the dominant sources of hydrogen reionization in the early Universe, with their high luminosities and surface temperatures expected to drive high ionizing photon production rates. In this work, we take our Geneva stellar evolution models of zero-metallicity stars and predict their production rates of photons capable to ionize H, He I and He II, based on a blackbody approximation. We present analytical fits in the range 1.7-500 solar masses. We then explore the impact of stellar initial mass, rotation, and convective overshooting for individual stars. We have found that ionizing photon production rates increase with increasing initial mass. For the rotational velocities considered we see changes of up to 25% to ionizing photons produced. This varies with initial mass and ionizing photon species and reflects changes to surface properties due to rotation. We have also found that higher convective overshooting increases ionizing photon production by approximately 20% for the change in overshooting considered here. For stellar populations, we explore how the production of ionizing photons varies as a function of the initial mass function (IMF) slope, and minimum and maximum initial masses. For a fixed population mass we have found changes of the order of 20-30% through varying the nature of the IMF. This work presents ionizing photon production predictions for the most up to date Geneva stellar evolution models of Population III stars, and provides insight into how key evolutionary parameters impact the contribution of the first stars to reionization., Comment: 20 pages, 18 figures. Accepted for publication in MNRAS

Using {\it Gaia} Early Data Release 3 (EDR3) parallaxes and Bayesian inference, we infer a parallax of the Westerlund 1 (Wd1) cluster. We find a parallax of $0.34\pm{0.05}$ mas corresponding to a distance of $2.8^{+0.7}_{-0.6}$ kpc. The new {\it Gaia} EDR3 distance is consistent with our previous result using {\it Gaia} DR2 parallaxes. This confirms that Wd1 is less massive and older than previously assumed. Compared to DR2, the EDR3 individual parallax uncertainties for each star decreased by 30\%. However, the aggregate parallax uncertainty for the cluster remained the same. This suggests that the uncertainty is dominated by systematics, which is possibly due to crowding, motions within the cluster, or motions due to binary orbits., Comment: 3 pages, 1 figure

Understanding the nature of the first stars is key to understanding the early universe. With new facilities such as JWST we may soon have the first observations of the earliest stellar populations, but to understand these observations we require detailed theoretical models. Here we compute a grid of stellar evolution models using the Geneva code with the aim to improve our understanding of the evolution of zero-metallicity stars, with particular interest in how rotation affects surface properties, interior structure, and metal enrichment. We produce a range of models of initial masses (Mini) from 1.7 Msun to 120 Msun, focusing on massive models of 9 Msun < Mini < 120 Msun. Our grid includes models with and without rotation, with rotating models having an initial velocity of 40% of the critical velocity. We find that rotation strongly impacts the evolution of the first stars, mainly through increased core size and stronger H-burning shells during core He-burning. Without radiative mass loss, angular momentum builds at the surface in rotating models, thus models of initial masses Mini > 60 Msun reach critical rotation on the main sequence and experience mass loss. We find that rotational mixing strongly affects metal enrichment, but does not always increase metal production as we see at higher metallicities. This is because rotation leads to an earlier CNO boost to the H shell during He-burning, which may hinder metal enrichment depending on initial mass and rotational velocity. Electronic tables of this new grid of Population III models are publicly available., Comment: 20 pages, 15 figures. Accepted for publication in MNRAS

GW190521 challenges our understanding of the late-stage evolution of massive stars and the effects of the pair-instability in particular. We discuss the possibility that stars at low or zero metallicity could retain most of their hydrogen envelope until the pre-supernova stage, avoid the pulsational pair-instability regime and produce a black hole with a mass in the mass gap by fallback. We present a series of new stellar evolution models at zero and low metallicity computed with the Geneva and MESA stellar evolution codes and compare to existing grids of models. Models with a metallicity in the range 0-0.0004 have three properties which favour higher BH masses as compared to higher metallicity models. These are (i) lower mass-loss rates during the post-MS phase, (ii) a more compact star disfavouring binary interaction and (iii) possible H-He shell interactions which lower the CO core mass. We conclude that it is possible that GW190521 may be the merger of black holes produced directly by massive stars from the first stellar generations. Our models indicate BH masses up to 70-75 Msun. Uncertainties related to convective mixing, mass loss, H-He shell interactions and pair-instability pulsations may increase this limit to ~85 Msun., Comment: 6 pages, 3 figures, accepted for publication in MNRAS Letters

We introduce SNAPSHOT, a technique to systematically compute stellar structure models in hydrostatic and thermal equilibrium based on 3 structural properties - core mass $M_{\rm core}$, envelope mass $M_{\rm env}$ and core composition. This approach allows us to connect these properties of stellar interiors to the luminosity and effective temperature $T_{\rm eff}$ in a more systematic way than with stellar evolution models. For MS models, we derive an analytical relationship between $M_{\rm core}$, $M_{\rm total}$ and central H abundance that can be used in rapid stellar evolution algorithms. Core-He burning models with $M_{\rm core}/M_{\rm total}$ from 0.2 to 0.8 have convective envelopes, low $T_{\rm eff}$ and will appear as red supergiants. For a given $M_{\rm core}$, they exhibit a small variation in luminosity (0.02 dex) and $T_{\rm eff}$ ($\sim 400\,\mathrm{K}$) over a wide range of $M_{\rm env}$ ($\sim 2 - 20\,\mathrm{M}_{\odot}$). This means that it is not possible to derive red supergiant masses from luminosities and $T_{\rm eff}$ alone. We derive the following relationship between $M_{\rm core}$ and the total luminosity of a red supergiant during core He burning: $\log M_{\rm core} \simeq 0.44 \log L/L_{\odot} - 1.38$. At $M_{\rm core}$/$M_{\rm total} \approx 0.2$, our models exhibit a bi-stability and jump from a RSG to a BSG structure. Our models with $M_{\rm core}/M_{\rm total} > 0.8$, which correspond to stripped stars produced by mass loss or binary interaction, show that $T_{\rm eff}$ has a strong dependence on $M_{\rm env}$, $M_{\rm core}$ and the core composition. We find the mass of one of these stripped stars in a binary system, HD 45166, to be less than its dynamical mass. When a large observational sample of stripped stars becomes available, our results can be used to constrain their masses and the physics of binary interaction., Comment: Accepted for publication in MNRAS

We investigate a suspected very massive star in one of the most metal-poor dwarf galaxies, PHL~293B. Excitingly, we find the sudden disappearance of the stellar signatures from our 2019 spectra, in particular the broad H lines with P~Cygni profiles that have been associated with a massive luminous blue variable (LBV) star. Such features are absent from our spectra obtained in 2019 with the ESPRESSO and X-shooter instruments of the ESO's VLT. We compute radiative transfer models using CMFGEN that fit the observed spectrum of the LBV and are consistent with ground-based and archival Hubble Space Telescope photometry. Our models show that during 2001--2011 the LBV had a luminosity $L_* = 2.5-3.5 \times 10^6 ~L_{\odot}$, a mass-loss rate $\dot{M} = 0.005-0.020 ~M_{\odot}$~yr$^{-1}$, a wind velocity of 1000~km~s$^{-1}$, and effective and stellar temperatures of $T_\mathrm{eff} = 6000-6800$~K and $T_\mathrm{*}=9500-15000$~K. These stellar properties indicate an eruptive state. We consider two main hypotheses for the absence of the broad emission components from the spectra obtained since 2011. One possibility is that we are seeing the end of an LBV eruption of a surviving star, with a mild drop in luminosity, a shift to hotter effective temperatures, and some dust obscuration. Alternatively, the LBV could have collapsed to a massive black hole without the production of a bright supernova., Comment: 8, pages, 7 figures, MNRAS accepted ; see also the ESO press release at: https://www.eso.org/public/news/eso2010/

In circumstellar gas, the complex organic molecule methanol has been found almost exclusively around young stellar objects, and is thus regarded as a signpost of recent star formation. Here we report the first probable detection of methanol around an evolved high-mass star, in the complex circumstellar environment around the Luminous Blue Variable $\eta$ Carinae, while using ALMA to investigate molecular cloud conditions traced by CO (2-1) in an orbit phase of the massive binary preceding the 2020 periastron. Favoring methanol over a $^{13}$CS alternative, the emission originates from hot ($T_{\rm{gas}} \simeq$ 700 K) material, $\sim$2$''$ (0.02 pc) across, centered on the dust-obscured binary in contrast to the CO which traces inner layers of the extended massive equatorial torus, and is accompanied by prominent absorption in a cooler ($T_{\rm{gas}} \simeq$ 110 K) layer of gas. We also report detections of water in $Herschel$/HIFI observations at 557 GHz and 988 GHz. The methanol abundance is several to 50 times higher than observed towards several lower mass stars, while water abundances are similar to those observed in cool, dense molecular clouds. The very high methanol:water abundance ratio in the core of $\eta$ Carinae may suggest methanol formation processes similar to Fischer-Tropsch-type catalytic reactions on dust grains. These observations prove that complex molecule formation can occur in the chemically evolved environments around massive stars in the end stages of their evolution, given sufficient gas densities and shielding conditions as may occur in material around massive interacting companions and merger remnants., Comment: Accepted by ApJ Letters

We show that it is not possible to determine the final mass $M_{\rm fin}$ of a red supergiant (RSG) at the pre-supernova (SN) stage from its luminosity $L$ and effective temperature $T_{\rm eff}$ alone. Using a grid of stellar models, we demonstrate that for a given value of $L$ and $T_{\rm eff}$, a RSG can have a range of $M_{\rm fin}$ as wide as 3 to $45~\mathrm{M}_{\odot}$. While the probability distribution within these limits is not flat, any individual determination of $M_{\rm fin}$ for a RSG will be degenerate. This makes it difficult to determine its evolutionary history and to map $M_{\rm fin}$ to an initial mass. Single stars produce a narrower range that is difficult to accurately determine without making strong assumptions about mass loss, convection, and rotation. Binaries would produce a wider range of RSG $M_{\rm fin}$. However, the final Helium core mass M$_{\rm He-core}$ is well determined by the final luminosity and we find $\log (\mathrm{M}_{\rm He-core}/M_{\odot}) = 0.659 \log (L/\mathrm{L}_{\odot}) -2.630$ Using this relationship, we derive M$_{\rm He-core}$ for directly imaged SN progenitors and one failed SN candidate. The value of $M_{\rm fin}$ for stripped star progenitors of SNe IIb is better constrained by $L$ and $T_{\rm eff}$ due to the dependence of $T_{\rm eff}$ on the envelope mass $M_{\rm env}$ for $M_{\rm env} \lesssim 1~$M$_{\odot}$. Given the initial mass function, our results apply to the majority of progenitors of core collapse SNe, failed SNe and direct collapse black holes., Comment: Re-submitted to MNRAS Letters after referee comments

We compute an extensive set of early-time spectra of supernovae interacting with circumstellar material using the radiative transfer code CMFGEN. Our models are applicable to events observed from 1 to a few days after explosion. Using these models, we constrain the progenitor and explosion properties of a sample of 17 observed interacting supernovae at early-times. Because massive stars have strong mass loss, these spectra provide valuable information about supernova progenitors, such as mass-loss rates, wind velocities, and surface abundances. We show that these events span a wide range of explosion and progenitor properties, exhibiting supernova luminosities in the 1e8 to 1e12 Lsun range, temperatures from 10000 to 60000 K, progenitor mass-loss rates from a few 1e-4 up to 1 Msun/yr, wind velocities from 100 to 800 km/s, and surface abundances from solar-like to H-depleted. Our results suggest that many progenitors of supernovae interacting with circumstellar material have significantly increased mass-loss before explosion compared to what massive stars show during the rest of their lifetimes. We also infer a lack of correlation between surface abundances and mass-loss rates. This may point to the pre-explosion mass-loss mechanism being independent of stellar mass. We find that the majority of these events have CNO-processed surface abundances. In the single star scenario this points to a preference towards high-mass RSGs as progenitors of interacting SNe, while binary evolution could impact this conclusion. Our models are publicly available and readily applicable to analyze results from ongoing and future large scale surveys such as the Zwicky Transient Factory., Comment: Accepted for publication in MNRAS. All the synthetic spectra presented in this article are publicly available and can be downloaded via WISeREP (https://wiserep.weizmann.ac.il/). 32 pages, 25 figures

We investigate the effects of mass loss during the main-sequence (MS) and post-MS phases of massive star evolution on black hole (BH) birth masses. We compute solar metallicity Geneva stellar evolution models of an 85 $M_{\odot}$ star with mass-loss rate ($\dot{M}$) prescriptions for MS and post-MS phases and analyze under which conditions such models could lead to very massive BHs. Based on the observational constraints for $\dot{M}$ of luminous stars, we discuss two possible scenarios that could produce massive BHs at high metallicity. First, if a massive BH progenitor evolves from the observed population of massive MS stars known as WNh stars, we show that its average post-MS mass-loss rate has to be less than $1\,\times10^{-5}\,M_{\odot}$/yr. However, this is lower than the typical observed mass-loss rates of luminous blue variables (LBV). Second, a massive BH progenitor could evolve from a yet undetected population of $80-85$ $M_{\odot}$ stars with strong surface magnetic fields, which could quench mass loss during the evolution. In this case, the average mass-loss rate during the post-MS LBV phase has to be less than $5\,\times10^{-5}\,M_{\odot}$/yr to produce 70 $M_{\odot}$ BHs. We suggest that LBVs that explode as SNe have large envelopes and small cores that could be prone to explosion, possibly evolving from binary interaction (either mergers or mass gainers that do not fully mix). Conversely, LBVs that directly collapse to BHs could have evolve from massive single stars or binary-star mergers that fully mix, possessing large cores that would favor BH formation., Comment: Accepted in ApJ

Massive stars are the drivers of the chemical evolution of dwarf galaxies. We review here the basics of massive star evolution and the specificities of stellar evolution in low-Z environment. We discuss nucleosynthetic aspects and what observations could constrain our view on the first generations of stars., Comment: 8 pages, 4 figures. Proceedings of IAUS 344 "Dwarf Galaxies: From the Deep Universe to the Present" held in Vienna (Austria) August 20-24 2018

The luminous blue variable (LBV) RMC143 is located in the outskirts of the 30~Doradus complex, a region rich with interstellar material and hot luminous stars. We report the $3\sigma$ sub-millimetre detection of its circumstellar nebula with ALMA. The observed morphology in the sub-millimetre is different than previously observed with HST and ATCA in the optical and centimetre wavelength regimes. The spectral energy distribution (SED) of RMC143 suggests that two emission mechanisms contribute to the sub-mm emission: optically thin bremsstrahlung and dust. Both the extinction map and the SED are consistent with a dusty massive nebula with a dust mass of $0.055\pm0.018~M_{\odot}$ (assuming $\kappa_{850}=1.7\rm\,cm^{2}\,g^{-1}$). To date, RMC143 has the most dusty LBV nebula observed in the Magellanic Clouds. We have also re-examined the LBV classification of RMC143 based on VLT/X-shooter spectra obtained in 2015/16 and a review of the publication record. The radiative transfer code CMFGEN is used to derive its fundamental stellar parameters. We find an effective temperature of $\sim 8500$~K, luminosity of log$(L/L_{\odot}) = 5.32$, and a relatively high mass-loss rate of $1.0 \times 10^{-5}~M_{\odot}$~yr$^{-1}$. The luminosity is much lower than previously thought, which implies that the current stellar mass of $\sim8~M_{\odot}$ is comparable to its nebular mass of $\sim 5.5~M_{\odot}$ (from an assumed gas-to-dust ratio of 100), suggesting that the star has lost a large fraction of its initial mass in past LBV eruptions or binary interactions. While the star may have been hotter in the past, it is currently not hot enough to ionize its circumstellar nebula. We propose that the nebula is ionized externally by the hot stars in the 30~Doradus star-forming region., Comment: Paper accepted by A&A on 09/05/2019 and in proof stage. Second comments by referee are included in this version

The effects of rotation on stellar evolution are particularly important at low metallicity, when mass loss by stellar winds diminishes and the surface enrichment due to rotational mixing becomes relatively more pronounced than at high metallicities. Here we investigate the impact of rotation and metallicity on stellar evolution. Using a similar physics as in our previous large grids of models at Z=0.002 and Z=0.014, we compute stellar evolution models with the Geneva code for rotating and nonrotating stars with initial masses (Mini) between 1.7 and 120 Msun and Z=0.0004 (1/35 solar). This is comparable to the metallicities of the most metal poor galaxies observed so far, such as I Zw 18. Concerning massive stars, both rotating and nonrotating models spend most of their core-helium burning phase with an effective temperature higher than 8000 K. Stars become red supergiants only at the end of their lifetimes, and few RSGs are expected. Our models predict very few to no classical Wolf-Rayet stars as a results of weak stellar winds at low metallicity. The most massive stars end their lifetimes as luminous blue supergiants or luminous blue variables, a feature that is not predicted by models with higher metallicities. Interestingly, due to the behavior of the intermediate convective zone, the mass domain of stars producing pair-instability supernovae is smaller at Z=0.0004 than at Z=0.002. We find that during the main sequence phase, the ratio between nitrogen and carbon abundances (N/C) remains unchanged for nonrotating models. However, N/C increases by factors of 10-20 in rotating models at the end of the MS. Cepheids coming from stars with Mini > 4-6 Msun are beyond the core helium burning phase and spend little time in the instability strip. Since they would evolve towards cooler effective temperatures, these Cepheids should show an increase of the pulsation period as a function of age., Comment: 19 pages, 18 figures, accepted for publication in A&A

Westerlund 1 (Wd1) is potentially the largest star cluster in the Galaxy. That designation critically depends upon the distance to the cluster, yet the cluster is highly obscured, making luminosity-based distance estimates difficult. Using {\it Gaia} Data Release 2 (DR2) parallaxes and Bayesian inference, we infer a parallax of $0.35^{+0.07}_{-0.06}$ mas corresponding to a distance of $2.6^{+0.6}_{-0.4}$ kpc. To leverage the combined statistics of all stars in the direction of Wd1, we derive the Bayesian model for a cluster of stars hidden among Galactic field stars; this model includes the parallax zero-point. Previous estimates for the distance to Wd1 ranged from 1.0 to 5.5 kpc, although values around 5 kpc have usually been adopted. The {\it Gaia} DR2 parallaxes reduce the uncertainty from a factor of 3 to 18\% and rules out the most often quoted value of 5 kpc with 99\% confidence. This new distance allows for more accurate mass and age determinations for the stars in Wd1. For example, the previously inferred initial mass at the main-sequence turn-off was around 40 M$_{\odot}$; the new {\it Gaia} DR2 distance shifts this down to about 22 M$_{\odot}$. This has important implications for our understanding of the late stages of stellar evolution, including the initial mass of the magnetar and the LBV in Wd1. Similarly, the new distance suggests that the total cluster mass is about four times lower than previously calculated., Comment: 14 pages, 10 figures

A large fraction of massive stars evolve in interacting binary systems, which dramatically modifies the outcome of stellar evolution. We investigated the properties of blue supergiants in binary systems and whether they are suitable for extragalactic distance determinations using the flux-weighted gravity luminosity relationship (FGLR). This is a relationship between the absolute bolometric magnitude $M_{\rm bol}$ and the spectroscopically determined flux-weighted gravity $g_{\rm F} = g/T^4_{\rm eff}$. We computed a grid of binary stellar evolution models with MESA and use the v2.1 BPASS models to examine whether they are compatible with the relatively small scatter shown by the observed relationship. Our models have initial primary masses of 9 - 30 $M_\odot$, initial orbital periods of 10 - 2511 days, mass ratio $q$ = 0.9, and metallicity $Z$ = 0.02. We find that the majority of primary stars that produce blue supergiant stages are consistent with the observed FGLR, with a small offset towards brighter bolometric magnitudes. In between 1\% - 24\% of cases, binary evolution may produce blue supergiants after a mass transfer episode, that lie below the observed FGLR. A small number of such stars have been found in extragalactic FGLR studies, suggesting that they may have evolved through binary interaction. Some models with shorter periods could resemble blue hypergiants and LBVs. We used CMFGEN radiative transfer models to investigate the effects of unresolved secondaries on diagnostics for $T_{\rm eff}$ and $g$, and the biases on the determination of interstellar reddening and $M_{\rm bol}$. We find that the effects are small and within the observed scatter, but could lead to a small overestimate of the luminosity, of $T_{\rm eff}$ and of $g$ for extreme cases. We conclude that the observed FGLR can, in principle, be well reproduced by close binary evolution models., Comment: 15 pages, 12 figures, accepted for publication in A&A

We examine parallaxes and distances for Galactic luminous blue variables (LBVs) in Gaia DR2. The sample includes 11 LBVs and 14 LBV candidates. For about half of the sample, DR2 distances are either similar to commonly adopted literature values, or the DR2 values have large uncertainties. For the rest, reliable DR2 distances differ significantly from values in the literature, and in most cases the Gaia DR2 distance is smaller. Two key results are that the S Doradus instability strip may not be as clearly defined as previously thought, and that there exists a population of LBVs at relatively low luminosities. LBVs seem to occupy a wide swath from the end of the main sequence at the blue edge to 8000 K at the red side, with a spread in luminosity reaching as low as log(L/Lsun)=4.5. The lower-luminosity group corresponds to effective single-star initial masses of 10-25 Msun, and includes objects that have been considered as confirmed LBVs. We discuss implications for LBVs including (1) their instability and origin in binary evolution, (2) connections to some supernova (SN) impostors such as the class of SN 2008S-like objects, and (3) LBVs that may be progenitors of SNe with dense circumstellar material across a wide initial mass range. Although some of the Gaia DR2 distances for LBVs have large uncertainty, this represents the most direct and consistent set of Galactic LBV distance estimates available in the literature., Comment: 20 pages, 5 figures, MNRAS accepted. Updated inclusion of excess astrometric noise and discussion of literature distances

In this paper we analyse the pre-explosion spectrum of SN2015bh by performing radiative transfer simulations using the CMFGEN code. This object has attracted significant attention due to its remarkable similarity to SN2009ip in both its pre- and post-explosion behaviour. They seem to belong to a class of events for which the fate as a genuine core-collapse supernova or a non-terminal explosion is still under debate. Our CMFGEN models suggest that the progenitor of SN2015bh had an effective temperature between 8700 and 10000 K, luminosity in the range ~ 1.8-4.74e6 Lsun, contained at least 25% H in mass at the surface, and half-solar Fe abundances. The results also show that the progenitor of SN 2015bh generated an extended wind with a mass-loss rate of ~ 6e-4 to 1.5e-3 Msun/yr and a velocity of 1000 km/s. We determined that the wind extended to at least 2.57e14 cm and lasted for at least 30 days prior to the observations, releasing 5e-5 Msun into the circumstellar medium. In analogy to 2009ip, we propose that this is the material that the explosive ejecta could interact at late epochs, perhaps producing observable signatures that can be probed with future observations. We conclude that the progenitor of SN 2015bh was most likely a warm luminous blue variable of at least 35 Msun before the explosion. Considering the high wind velocity, we cannot exclude the possibility that the progenitor was a Wolf-Rayet star that inflated just before the 2013 eruption, similar to HD5980 during its 1994 episode. If the star survived, late-time spectroscopy may reveal either a similar LBV or a Wolf-Rayet star, depending on the mass of the H envelope before the explosion. If the star exploded as a genuine SN, 2015bh would be a remarkable case of a successful explosion after black-hole formation in a star with a possible minimum mass 35 Msun at the pre-SN stage., Comment: 13 pages, 10 figures, accepted for publication in A&A

Eta Carinae, the closest, active, massive binary containing a highly unstable LBV, exhibits expanding, compressed wind shells, seen in emission, that are spatially and spectrally resolved by \hst/\stis. Starting in June 2009, these structures were mapped across its 5.54-year, highly elliptical, binary orbit to follow temporal changes in the light of [Fe~III] 4659\AA\ and [Fe~II] 4815\AA. The emissions trace portions of fossil wind-shells, that were formed by wind-wind interactions across each cycle. Over the high-ionization state, dense arcs, photo-ionized by far ultraviolet radiation from the hot secondary, are seen in [Fe~III]. Other arcs, ionized by mid-ultraviolet radiation from the primary star, are seen in [Fe II]. The [Fe III] structures tend to be interior to [Fe II] structures that trace extensive, less disturbed primary wind. During the brief periastron passage when the secondary plunges deep into the primary's extremely dense wind, on the far side of primary star, high-ionization [Fe III] structures fade and reappear in [Fe II]. Multiple fossil wind-structures were traced across the 5.7-year monitoring interval. The strong similarity of the expanding [Fe II] shells suggests that the wind and photo-ionization properties of the massive binary have not changed substantially from one orbit to the next over the past several orbital cycles. These observations trace structures that can be used to test three-dimensional hydrodynamical and radiative-transfer models of massive, interacting winds. They also provide a baseline for following future changes in Eta Carinae, especially of its winds and photoionization properties., Comment: 28 pages accepted MNRAS July 22 2016

The Galactic blue supergiant SBW1 with its circumstellar ring nebula represents the best known analog of the progenitor of SN 1987A. High-resolution imaging has shown H-alpha and IR structures arising in an ionized flow that partly fills the ring's interior. To constrain the influence of the stellar wind on this structure, we obtained an ultraviolet (UV) spectrum of the central star of SBW1 with the HST Cosmic Origins Spectrograph (COS). The UV spectrum shows none of the typical wind signatures, indicating a very low mass-loss rate. Radiative transfer models suggest an extremely low rate below 10$^{-10}$ Msun/yr, although we find that cooling timescales probably become comparable to or longer than the flow time below 10$^{-8}$ Msun/yr. We therefore adopt this latter value as a conservative upper limit. For the central star, the model yields $T_{\rm eff}$=21,000$\pm$1000 K, $L\simeq$5$\times$10$^4$ $L_{\odot}$, and roughly Solar composition except for enhanced N abundance. SBW1's very low mass-loss rate may hinder the wind's ability to shape the surrounding nebula. The very low mass-loss rate also impairs the wind's ability to shed angular momentum; the spin-down timescale for magnetic breaking is more than 500 times longer than the age of the ring. This, combined with the star's slow rotation rate, constrain merger scenarios to form ring nebulae. The mass-loss rate is at least 10 times lower than expected from mass-loss recipes, without any account of clumping. The physical explanation for why SBW1's wind is so weak presents an interesting mystery., Comment: 12 pages, 6 figs. submitted to MNRAS. comments welome

We have monitored the massive binary star $\eta$ Carinae with the CTIO/SMARTS 1.5~m telescope and CHIRON spectrograph from the previous apastron passage of the system through the recent 2014.6 periastron passage. Our monitoring has resulted in a large, homogeneous data set with an unprecedented time-sampling, spectral resolving power, and signal-to-noise. This allowed us to investigate temporal variability previously unexplored in the system and discover a kinematic structure in the P Cygni absorption troughs of neutral helium wind lines. The features observed occurred prior to the periastron passage and are seen as we look through the trailing arm of the wind-wind collision shock cone. We show that the bulk of the variability is repeatable across the last five periastron passages, and that the absorption occurs in the inner 230 AU of the system. In addition, we found an additional, high-velocity absorption component super-imposed on the P Cygni absorption troughs that has been previously un-observed in these lines, but which bears resemblance to the observations of the He~I $\lambda$10830 \AA\ feature across previous cycles. Through a comparison of the current smoothed particle hydrodynamical simulations, we show that the observed variations are likely caused by instabilities in the wind-wind collision region in our line-of-sight, coupled with stochastic variability related to clumping in the winds., Comment: accepted to MNRAS, 28 pages including 14 figures and a 9 page appendix

The evolved, massive highly eccentric binary system, eta Carinae, underwent a periastron passage in the summer of 2014. We obtained two coordinated X-ray observations with XMM-Newton and NuSTAR during the elevated X-ray flux state and just before the X-ray minimum flux state around this passage. These NuSTAR observations clearly detected X-ray emission associated with eta Car extending up to ~50 keV for the first time. The NuSTAR spectrum above 10 keV can be fit with the bremsstrahlung tail from a kT ~6 keV plasma. This temperature is Delta kT ~2 keV higher than those measured from the iron K emission line complex, if the shocked gas is in collisional ionization equilibrium. This result may suggest that the companion star's pre-shock wind velocity is underestimated. The NuSTAR observation near the X-ray minimum state showed a gradual decline in the X-ray emission by 40% at energies above 5 keV in a day, the largest rate of change of the X-ray flux yet observed in individual eta Car observations. The column density to the hardest emission component, NH ~1e24 cm-2, marked one of the highest values ever observed for eta Car, strongly suggesting increased obscuration of the wind-wind colliding X-ray emission by the thick primary stellar wind prior to superior conjunction. Neither observation detected the power-law component in the extremely hard band that INTEGRAL and Suzaku observed prior to 2011. If the non-detection by NuSTAR is caused by absorption, the power-law source must be small and located very near the WWC apex. Alternatively, it may be that the power-law source is not related to either eta Car or the GeV gamma-ray source., Comment: 24 pages, 6 figures, 2 tables

Wolf-Rayet (WR) stars, as they are advanced stages of the life of massive stars, provide a good test for various physical processes involved in the modelling of massive stars, such as rotation and mass loss. In this paper, we show the outputs of the latest grids of single massive stars computed with the Geneva stellar evolution code, and compare them with some observations. We present a short discussion on the shortcomings of single stars models and we also briefly discuss the impact of binarity on the WR populations., Comment: 3 pages, 1 figure, to appear in the proceedings of the international Wolf-Rayet stars workshop held in Potsdam (2015)

We present the first quantitative spectroscopic modeling of an early-time supernova that interacts with its progenitor wind. Using the radiative transfer code CMFGEN, we investigate the recently-reported 15.5 h post-explosion spectrum of the type IIb SN 2013cu. For the first time, we are able to directly measure the chemical abundances of a SN progenitor and find a relatively H-rich wind, with H and He abundances (by mass) of X=0.46 +- 0.2 and Y=0.52 +- 0.2, respectively. The wind is enhanced in N and depleted in C relative to solar values (mass fractions of 8.2e-3 and 1e-5). We obtain that a dense wind/circumstellar medium, with a mass-loss rate of Mdot= 3e-3 Msun/yr and wind velocity vwind=100 km/s, surrounds the star at the pre-SN stage. These values are lower than previous analytical estimates, although we find Mdot/vinf consistent with previous work. We also compute a CMFGEN model to constrain the progenitor spectral type and find that the high Mdot and low vwind imply that the star had an effective temperature of ~8000 K immediately before the SN explosion. Our models suggest that the progenitor was either an unstable luminous blue variable or a yellow hypergiant undergoing an eruptive phase, and rule out a WR star. We classify the post-explosion spectra at 15.5 h as XWN5(h) and advocate for the use of the prefix `X' (eXplosion) to avoid confusion between post-explosion, non-stellar spectra with those of massive stars. We show that the progenitor spectral type is significantly different than the early post-explosion spectral type owing to the huge differences in the ionization structure before and after the SN event. We find the following temporal evolution: LBV/YHG -> XWN5(h) -> SN IIb. Future early-time spectroscopy in the UV will give access to additional spectroscopic diagnostics and further constrain the properties of SN precursors, such as their metallicities., Comment: 4 pages, 3 figures, submitted to A&A, comments welcome

The well-studied Type IIn supernova (SN) 1998S is often dubbed the prototypical SN IIn, and it provides a unique opportunity to study its progenitor star from within as the SN lights up dense circumstellar material (CSM) launched from the progenitor. Here we present a Keck HIRES spectrum of SN 1998S taken within a few days after core collapse - both the earliest high-resolution ($\Delta\lambda < 1.0\,\AA$) spectrum published of a SN IIn and the earliest published spectrum of SN 1998S. Modern SN studies achieve impressively short turn-around times between SN detection and the first observed spectrum, but high-resolution spectra of very young supernovae are rare; the unique spectrum presented here provides a useful case study for observations of other young SN systems including SN 2013cu, which displayed a remarkably similar spectrum when very young. We examine the fully resolved emission-line profiles of SN 1998S, finding evidence for extreme mass loss from the progenitor at velocities much less than those characteristic of Wolf-Rayet (WR) stars. We model our high-resolution SN 1998S spectrum using the radiative-transfer code CMFGEN and explore the composition, density, and velocity gradients within the SN system. We find a mass-loss rate of $6.0 \times 10^{-3} M_{\odot}$ per year during the $\sim$15 yr before core collapse, while other studies indicate a much lower rate at earlier times (>15yr before core collapse). A comparison with a spectrum of SN 2013cu indicates many similarities, though SN 2013cu was of Type IIb - indicating that very different supernovae can arise from progenitors with extreme mass loss in the last few years before explosion., Comment: 33 pages, 10 figures

For the first time, the interior and spectroscopic evolution of a massive star is analyzed from the zero-age main sequence (ZAMS) to the pre-supernova (SN) stage. For this purpose, we combined stellar evolution models using the Geneva code and atmospheric models using CMFGEN. With our approach, we were able to produce observables, such as a synthetic high-resolution spectrum and photometry, aiding the comparison between evolution models and observed data. Here we analyze the evolution of a non-rotating 60 Msun star and its spectrum throughout its lifetime. Interestingly, the star has a supergiant appearance (luminosity class I) even at the ZAMS. We find the following evolutionary sequence of spectral types: O3 I (at the ZAMS), O4 I (middle of the H-core burning phase), B supergiant (BSG), B hypergiant (BHG), hot luminous blue variable (LBV; end of H-core burning), cool LBV (H-shell burning through the beginning of the He-core burning phase), rapid evolution through late WN and early WN, early WC (middle of He-core burning), and WO (end of He-core burning until core collapse). We find the following spectroscopic phase lifetimes: 3.22e6 yr for the O-type, 0.34e5 yr (BSG), 0.79e5 yr (BHG), 2.35e5 yr (LBV), 1.05e5 yr (WN), 2.57e4 yr (WC), and 3.80e4 yr (WO). Compared to previous studies, we find a much longer (shorter) duration for the early WN (late WN) phase, as well as a long-lived LBV phase. We show that LBVs arise naturally in single-star evolution models at the end of the MS when the mass-loss rate increases as a consequence of crossing the bistability limit. We discuss the evolution of the spectra, magnitudes, colors, and ionizing flux across the star's lifetime, and the way they are related to the evolution of the interior. [abridged], Comment: 23 pages, 10 figures, accepted for publication in Astronomy and Astrophysics

Recent work suggests that the mass-loss rate of the primary star (Eta A) in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present results from large- (r=1545au) and small- (r=155au) domain, 3D smoothed particle hydrodynamic (SPH) simulations of Eta Car's colliding winds for 3 Eta A mass-loss rates (2.4, 4.8, and 8.5 x 10^-4 M_sun/yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling, and radiative forces. We find that Eta A's mass-loss rate greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star (Eta B) switches from the adiabatic to the radiative-cooling regime during periastron passage. The SPH simulations together with 1D radiative transfer models of Eta A's spectra reveal that a factor of 2 or more drop in Eta A's mass-loss rate should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in Eta A's mass-loss rate was likely by a factor < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line-of-sight to the system lies close to the dense walls of the WWC zone. A modest decrease in Eta A's mass-loss rate may be responsible, but changes in the wind/stellar parameters of Eta B cannot yet be fully ruled out. We suggest observations during Eta Car's next periastron in 2014 to further test for decreases in Eta A's mass-loss rate. If Eta A's mass-loss rate is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009., Comment: 38 pages, 25 figures, 1 table. Accepted for publication in MNRAS

We discuss four questions dealing with massive star evolution. The first one is about the origin of slowly rotating, non-evolved, nitrogen rich stars. We propose that these stars may originate from initially fast rotating stars whose surface has been braked down. The second question is about the evolutionary status of alpha-Cygni variables. According to their pulsation properties, these stars should be post red supergiant stars. However, some stars at least present surface abundances indicating that they should be pre red supergiant stars. How to reconcile these two contradictory requirements? The third one concerns the various supernova types which are the end point of the evolution of stars with initial masses between 18 and 30 Msun, i.e. the most massive stars which go through a red supergiant phase during their lifetime. Do they produce types IIP, IIL, IIn, IIb or Ib supernovae or do they end without producing any SN event? Finally, we shall discuss reasons why so few progenitors of type Ibc supernovae have yet been detected?, Comment: 10 pages, 6 figures, conference, New advances in stellar physics: from microscopic to macroscopic processes, EAS publication series

We investigate the fundamental properties of core-collapse Supernova (SN) progenitors from single stars at solar metallicity. We combine Geneva stellar evolutionary models with initial masses of Mini=20-120 Msun with atmospheric/wind models using CMFGEN. We provide synthetic photometry and high-resolution spectra of hot stars at the pre-SN stage. For Mini=9-20 Msun, we supplement our analysis using publicly available MARCS model atmospheres of RSGs. We employ observational criteria of spectroscopic classification and find that massive stars, depending on Mini and rotation, end their lives as red supergiants (RSG), yellow hypergiants (YHG), luminous blue variables (LBV), and Wolf-Rayet (WR) stars of the WN and WO spectral types. For rotating models, we obtain the following types of SN progenitors: WO1-3 (Mini <= 32 Msun), WN10-11 (25 < Mini < 32 Msun), LBV (20 <= Mini < 25 Msun), G1 Ia+ (18 < Mini < 20 Msun), and RSGs (9 <= Mini <= 18 Msun). For non-rotating models, we find spectral types WO1-3 (Mini>40 Msun), WN7-8 (25 < Mini <= 40 Msun), WN11h/LBV (20 < Mini <= 25 Msun), and RSGs (9 <= Mini <= 20 Msun). Our rotating models indicate that SN IIP progenitors are all RSG, SN IIL/b progenitors are 56% LBVs and 44% YHGs, SN Ib progenitors are 96% WN10-11 and 4% WOs, and SN Ic progenitors are all WO stars. We find that not necessarily the most massive and luminous SN progenitors are the brighter ones in a given filter. We show that SN IIP progenitors (RSGs) are bright in the RIJHK_S filters and faint in the UB filters. SN IIL/b progenitors (LBVs and YHGs), and SN Ib progenitors (WNs) are relatively bright in optical/IR filters, while SN Ic progenitors (WOs) are faint in all optical filters. We argue that SN Ib and Ic progenitors from single stars should be undetectable in the available pre-explosion images with the current magnitude limits, in agreement with observational results.[abridged], Comment: 24 pages, 13 figures, accepted for publication in Astronomy and Astrophysics

(shortened) We provide a grid of single star models covering a mass range from 0.8 to 120 Msun with an initial metallicity Z = 0.002 with and without rotation. We discuss the impact of a change in the metallicity by comparing the current tracks with models computed with exactly the same physical ingredients but with a metallicity Z = 0.014 (solar). We show that the width of the main-sequence (MS) band in the upper part of the Hertzsprung-Russell diagram (HRD), for luminosity above log(L/Lsun) > 5.5, is very sensitive to rotational mixing. Strong mixing significantly reduces the MS width. We confirm, but here for the first time on the whole mass range, that surface enrichments are stronger at low metallicity provided that comparisons are made for equivalent initial mass, rotation and evolutionary stage. We show that the enhancement factor due to a lowering of the metallicity (all other factors kept constant) increases when the initial mass decreases. Present models predict an upper luminosity for the red supergiants (RSG) of log (L/Lsun) around 5.5 at Z = 0.002 in agreement with the observed upper limit of RSG in the Small Magellanic Cloud. We show that models using shear diffusion coefficient calibrated to reproduce the surface enrichments observed for MS B-type stars at Z = 0.014 can also reproduce the stronger enrichments observed at low metallicity. In the framework of the present models, we discuss the factors governing the timescale of the first crossing of the Hertzsprung gap after the MS phase. We show that any process favouring a deep localisation of the H-burning shell (steep gradient at the border of the H-burning convective core, low CNO content) and/or the low opacity of the H-rich envelope favour a blue position in the HRD for the whole or at least a significant fraction of the core He-burning phase., Comment: 17 pages, 15 figures, 2 tables. Accepted for publication in Astronomy and Astrophysics

Core-collapse supernova (SN) explosions mark the end of the tumultuous life of massive stars. Determining the nature of their progenitors is a crucial step towards understanding the properties of SNe. Until recently, no progenitor has been directly detected for SN of type Ibc, which are believed to come from massive stars that lose their hydrogen envelope through stellar winds and from binary systems where the companion has stripped the H envelope from the primary. Here we analyze recently reported observations of iPTF13bvn, which could possibly be the first detection of a SN Ib progenitor based on pre-explosion images. Very interestingly, the recently published Geneva models of single stars can reproduce the observed photometry of the progenitor candidate and its mass-loss rate, confirming a recently proposed scenario. We find that a single WR star with initial mass in the range 31-35 Msun fits the observed photometry of the progenitor of iPTF13bvn. The progenitor likely has a luminosity of log (L/Lsun)~5.55, surface temperature ~45000 K, and mass of ~10.9 Msun at the time of explosion. Our non-rotating 32 Msun model overestimates the derived radius of the progenitor, although this could likely be reconciled with a fine-tuned model of a more massive (between 40 and 50 Msun), hotter, and luminous progenitor. Our models indicate a very uncertain ejecta mass of ~8 Msun, which is higher than the average of the SN Ib ejecta mass that is derived from the lightcurve (2-4 Msun). This possibly high ejecta mass could produce detectable effects in the iPTF13bvn lightcurve and spectrum. If the candidate is indeed confirmed to be the progenitor, our results suggest that stars with relatively high initial masses (>30 Msun) can produce visible SN explosions at their deaths and do not collapse directly to a black hole., Comment: 4 pages, 2 figures, accepted for publication in A&A

Ideally, one would like to know which type of core-collapse SNe is produced by different progenitors and the channels of stellar evolution leading to these progenitors. These links have to be very well known to use the observed frequency of different types of SN events for probing the star formation rate and massive star evolution in different types of galaxies. We investigate the link between LBV as SN progenitors and the appearance of episodic radio light curve modulations of the SN event. We use the 20Msun and 25Msun models with rotation at solar metallicity, part of an extended grid of stellar models computed by the Geneva team. At their pre-SN stage, these two models have recently been shown to have spectra similar to those of LBV stars and possibly explode as Type IIb SNe. Based on the wind properties before the explosion, we derive the density structure of their circumstellar medium. This structure is used as input for computing the SN radio light curve. We find that the 20Msun model shows radio light curves with episodic luminosity modulations, similar to those observed in some Type IIb SNe. This occurs because the evolution of the 20Msun model terminates in a region of the HR diagram where radiative stellar winds present strong density variations, caused by the bistability limit. The 25Msun model, ending its evolution in a zone of the HR diagram where no change of the mass-loss rates is expected, presents no such modulations in its radio SN light curve. Our results reinforce the link between SN progenitors and LBV stars. We also confirm the existence of a physical mechanism for a single star to have episodic radio light curve modulations. In the case of the 25Msun progenitors, we do not obtain modulations in the radio light curve, but our models may miss some outbursting behavior in the late stages of massive stars., Comment: 5 pages, 3 figures, accepted by Astronomy & Astrophysics Letters

We review the significant role played by red supergiants (RSGs) in stellar populations, and some challenges and questions they raise for theoretical stellar evolution. We present how metallicity and rotation modify the way stars go to the red part of the Hertzsprung- Russell diagram or come back from it, and how RSGs might keep a trace of their main-sequence evolution. We compare theoretical popu- lation ratios with observed ones., Comment: 11 pages, 5 figures, Betelgeuse workshop, November 2012, Paris. To be published in the European Astronomical Society Publications Series, editors: Pierre Kervella, Thibaut Le Bertre & Guy Perrin

We discuss the physics and the evolution of a typical massive star passing through an evolutionary stage similar to that of Betelgeuse. After a brief introduction recalling various observed parameters of Betelgeuse, we discuss the Pre-Main-Sequence phase (PMS), the Main-Sequence (MS) phase, the physics governing the duration of the first crossing of the HR diagram, the red supergiant stage (RSG), the post-red supergiant phases and the final fate of solar metallicity stars with masses between 9 and 25 M$_\odot$. We examine the impact of different initial rotation and of various prescriptions for the mass loss rates during the red supergiant phase. We show that, whatever the initial rotation rate (chosen between 0 and 0.7$\times\upsilon_{\rm crit}$, $\upsilon_{\rm crit}$ being the surface equatorial velocity producing a centrifugal acceleration balancing exactly the gravity) and the mass loss rates during the RSG stage (varied between a standard value and 25 times that value), a 15 M$_\odot$ star always ends its lifetime as a RSG and explodes as a type II-P or II-L supernova., Comment: 12 pages, 5 figures, Betelgeuse workshop, November 2012, Paris. To be published in the European Astronomical Society Publications Series, editors: Pierre Kervella, Thibaut Le Bertre & Guy Perrin

Mass-loss rates are one of the most relevant parameters determining the evolution of massive stars. In particular, the rates at which the star loses mass during the red-supergiant (RSG) phase is the least constrained by the observations or theory. In this paper, we show how the mass loss during the RSG phase affects the later evolution of the star, as well as the final type of supernova towards which it leads. We also discuss some possibilities to discriminate between blue stars that went through a RSG phase and those which remained in the blue part of the Hertzsprung-Russell diagram., Comment: 9 pages, 3 figures, to be published in the proceedings of the Betelgeuse workshop held in Paris (November 2012) (EAS Publications Series)

Stars more massive than about 8 Msun end their lives as a Supernova (SN), an event of fundamental importance Universe-wide. Theoretically, these stars have been expected to be either at the red supergiant, blue supergiant, or Wolf-Rayet stage before the explosion. We performed coupled stellar evolution and atmospheric modeling of stars with initial masses between 20 Msun and 120 Msun. We found that the 20 Msun and 25 Msun rotating models, before exploding as SN, have spectra that do not resemble any of the aforementioned classes of massive stars. Rather, they have remarkable similarities with rare, unstable massive stars known as Luminous Blue Variables (LBV). While observations show that some SNe seem to have had LBVs as progenitors, no theoretical model had yet predicted that a star could explode at this stage. Our models provide theoretical support for relatively low-luminosity LBVs exploding as SN in the framework of single stellar evolution. This is a significant shift in paradigm, meaning that a fraction of LBVs could be the end stage of massive star evolution, rather than a transitory evolutionary phase. We suggest that type IIb SN could have LBV as progenitors, and a prime example could be SN 2008ax., Comment: 4 pages, 3 figures, accepted for publication in Astronomy and Astrophysics

The last part of SpS5 dealt with the circumstellar environment. Structures are indeed found around several types of massive stars, such as blue and red supergiants, as well as WRs and LBVs. As shown in the last years, the potential of IR for their study is twofold: first, IR can help discover many previously unknown nebulae, leading to the identification of new massive stars as their progenitors; second, IR can help characterize the nebular features. Current and new IR facilities thus pave the way to a better understanding of the feedback from massive stars., Comment: 10 pages, to be published in Highlights of Astronomy vol 16, 3rd part of the proceedings of special session #5 'IR view of massive stars' (ed. Y. Naze), a session which took place during the IAU general assembly 2012 ; see also http://www.gaphe.ulg.ac.be/IAU_XXVIII/prg.html

We analyze spatially resolved spectroscopic observations of the Eta Carinae binary system obtained with HST/STIS. Eta Car is enshrouded by the dusty Homunculus nebula, which scatters light emitted by the central binary and provides a unique opportunity to study a massive binary system from different vantage points. We investigate the latitudinal and azimuthal dependence of H$\alpha$ line profiles caused by the presence of a wind-wind collision (WWC) cavity created by the companion star. Using two-dimensional radiative transfer models, we find that the wind cavity can qualitatively explain the observed line profiles around apastron. Regions of the Homunculus which scatter light that propagated through the WWC cavity show weaker or no H alpha absorption. Regions scattering light that propagated through a significant portion of the primary wind show stronger P Cygni absorption. Our models overestimate the H alpha absorption formed in the primary wind, which we attribute to photoionization by the companion, not presently included in the models. We can qualitatively explain the latitudinal changes that occur during periastron, shedding light on the nature of Eta Car's spectroscopic events. Our models support the idea that during the brief period of time around periastron when the primary wind flows unimpeded toward the observer, H alpha absorption occurs in directions toward the central object and Homunculus SE pole, but not toward equatorial regions close to the Weigelt blobs. We suggest that observed latitudinal and azimuthal variations are dominated by the companion star via the WWC cavity, rather than by rapid rotation of the primary star., Comment: 5 pages, 4 figures. Accepted for publication in ApJL

We present 2D radiative transfer modeling of the Eta Carinae binary system accounting for the presence of a wind-wind collision (WWC) cavity carved in the optically-thick wind of the primary star. By comparing synthetic line profiles with HST/STIS spectra obtained near apastron, we show that the WWC cavity has a strong influence on multi-wavelength diagnostics. This influence is regulated by the modification of the optical depth in the continuum and spectral lines. We find that H-alpha, H-beta, and Fe II lines are the most affected by the WWC cavity, since they form over a large volume of the primary wind. These spectral lines depend on latitude and azimuth since, according to the orientation of the cavity, different velocity regions of a spectral line are affected. For 2D models with orientation corresponding to orbital inclination angle 110deg < i < 140deg and longitude of periastron 210deg < omega < 330deg, the blueshifted and zero-velocity regions of the line profiles are the most affected. These orbital orientations are required to simultaneously fit the UV and optical spectrum of Eta Car, for a half-opening angle of the cavity in the range 50-70deg. We find that the excess P-Cygni absorption seen in H-alpha, H-beta and optical Fe II lines in spherical models becomes much weaker or absent in the 2D models, in agreement with the observations. The observed UV spectrum of Eta Car, dominated by Fe II absorption lines, is superbly reproduced by our 2D cavity models. Small discrepancies still remain, as H-gamma and H-delta absorptions are overestimated by our models. We suggest that photoionization of the wind of the primary by the hot companion star is responsible for the weak absorption seen in these lines. Our CMFGEN models indicate that the primary star has a mass-loss rate of 8.5x10e-4 Msun/yr and wind terminal velocity of 420 km/s around the 2000 apastron., Comment: 20 pages, 14 figures, accepted for publication in MNRAS

We present a new model for the behavior of scattered time-dependent, asymmetric near-UV emission from the nearby ejecta of {\eta} Car. Using a 3-D hydrodynamical simulation of {\eta} Car's binary colliding winds, we show that the 3-D binary orientation derived by Madura et al. (2012) is capable of explaining the asymmetric near-UV variability observed in the Hubble Space Telescope Advanced Camera for Surveys/High Resolution Camera (HST ACS/HRC) F220W images of Smith et al. (2004b). Models assuming a binary orientation with i ~ 130 to 145 degrees, {\omega} ~ 230 to 315 degrees, PAz ~ 302 to 327 degrees are consistent with the observed F220W near-UV images. We find that the hot binary companion does not significantly contribute to the near-UV excess observed in the F220W images. Rather, we suggest that a bore-hole effect and the reduction of Fe II optical depths inside the wind-wind collision cavity carved in the extended photosphere of the primary star lead to the time-dependent directional illumination of circum-binary material as the companion moves about in its highly elliptical orbit., Comment: 14 pages, 4 figures, 1 table. Accepted for publication in ApJL

Eta Carinae is considered to be a massive colliding wind binary system with a highly eccentric (e \sim 0.9), 5.54-yr orbit. However, the companion star continues to evade direct detection as the primary dwarfs its emission at most wavelengths. Using three-dimensional (3-D) SPH simulations of Eta Car's colliding winds and radiative transfer codes, we are able to compute synthetic observables across multiple wavebands for comparison to the observations. The models show that the presence of a companion star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line profiles, as well as the ground-based photometric monitoring. Here, we focus on the Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in the dense primary wind, allowing increased escape of radiation from the hotter/deeper layers of the primary's extended wind photosphere. The results have important implications for interpretations of Eta Car's observables at multiple wavelengths., Comment: 5 pages, 4 figures, To be published in the proceedings of the meeting 'Four Decades of Research on Massive Stars' in honor of Tony Moffat, 11-15 July 2011, Saint-Michel-des-Saints, Quebec

We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission observed in Eta Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) simulations of Eta Car's binary colliding winds. Radiative transfer codes are used to generate synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA {\theta} that the orbital plane projection of the line-of-sight makes with the apastron side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously reproduce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38 degrees, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i \approx 130 to 145 degrees, {\theta} \approx -15 to +30 degrees, and an orbital axis projected on the sky at a PA \approx 302 to 327 degrees east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, Eta B, thus orbits clockwise on the sky and is on the observer's side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modeling to determine the stellar masses., Comment: 23 pages, 12 color figures, plus 2 online-only appendices (available in the /anc folder of the Source directory). Accepted for publication in MNRAS

We report new HST/STIS observations that map the high-ionization forbidden line emission in the inner arcsecond of Eta Car, the first that fully image the extended wind-wind interaction region of the massive colliding wind binary. These observations were obtained after the 2009.0 periastron at orbital phases 0.084, 0.163, and 0.323 of the 5.54-year spectroscopic cycle. We analyze the variations in brightness and morphology of the emission, and find that blue-shifted emission (-400 to -200 km s-1) is symmetric and elongated along the northeast-southwest axis, while the red-shifted emission (+100 to +200 km s-1) is asymmetric and extends to the north- northwest. Comparison to synthetic images generated from a 3-D dynamical model strengthens the 3-D orbital orientation found by Madura et al. (2011), with an inclination i \approx 138\degree, argument of periapsis {\omega} \approx 270\degree, and an orbital axis that is aligned at the same PA on the sky as the symmetry axis of the Homunculus, 312\degree. We discuss the potential that these and future mappings have for constraining the stellar parameters of the companion star and the long-term variability of the system., Comment: 10 pages, 4 figures, accepted Ap J Letters

Luminous Blue Variables (LBVs) oftentimes show double-troughed absorption in their strong Halpha lines, which are as yet not understood. Intriguingly, the feature has also been seen in the interacting supernova SN 2005gj, which was for this reason suggested to have an LBV progenitor. Our aims are to understand the double-troughed absorption feature in LBVs and investigate whether this phenomenon is related to wind variability. To this purpose, we perform time-dependent radiative transfer modeling using CMFGEN. We find that abrupt changes in the wind-terminal velocity - as expected from the bi-stability jump - are required to explain the double-troughed absorption profiles in LBVs. This strengthens scenarios that discuss the link between LBVs and SNe utilizing the progenitor's wind variability resulting from the bi-stability jump. We also discuss why the presence of double-troughed P-Cygni components may become an efficient tool to detect extra-galactic LBVs and how to analyze their mass-loss history on the basis of just one single epoch of spectral observations., Comment: 5 pages, 5 figures, accepted by A&A Letters

We show that the significantly different effective temperatures (Teff) achieved by the luminous blue variable AG Carinae during the consecutive visual minima of 1985-1990 (Teff= 22,800 K) and 2000-2001 (Teff=17,000 K) place the star on different sides of the bistability limit, which occurs in line-driven stellar winds around Teff ~ 21,000 K. Decisive evidence is provided by huge changes in the optical depth of the Lyman continuum in the inner wind as Teff changes during the S Dor cycle. These changes cause different Fe ionization structures in the inner wind. The bistability mechanism is also related to the different wind parameters during visual minima: the wind terminal velocity was 2-3 times higher and the mass-loss rate roughly two times smaller in 1985-1990 than in 2000-2003. We obtain a projected rotational velocity of 220 +- 50 km/s during 1985-1990 which, combined with the high luminosity (L=1.5x10^6 Lsun), puts AG Car extremely close to the Eddington limit modified by rotation (Omega-Gamma limit): for an inclination angle of 90 deg, Gamma_Omega > 1.0 for M< 60 Msun. Based on evolutionary models and mass budget, we obtain an initial mass of ~ 100 Msun and a current mass of 60-70 Msun for AG Car. Therefore, AG Car is close to, if not at, the Omega-Gamma limit during visual minimum. Assuming M=70 Msun, we find that Gamma_Omega decreases from 0.93 to 0.72 as AG Car expands toward visual maximum, suggesting that the star is not above the Eddington limit during maximum phases., Comment: 13 pages, 8 figures, accepted by ApJ