Source author record

Cyril Georgy

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

55works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

55 published item(s)

preprint2022arXiv

Grids of stellar models with rotation VI: Models from 0.8 to 120 $M_\odot$ at a metallicity Z = 0.006

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.

preprint2022arXiv

Grids of stellar models with rotation VII: Models from 0.8 to 300 M$_\odot$ at super-solar metallicity (Z = 0.020)

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.

preprint2022arXiv

The formation of the stripped envelope type II b Supernova progenitors: Rotation, Metallicity and Overshooting

Type IIb supernovae are believed to originate from core-collapse progenitors having kept only a very thin hydrogen envelope. We aim to explore how some physical factors, such as rotation, metallicity, overshooting, and the initial orbital period in binaries, significantly affect the Roche lobe overflow and the formation of type IIb supernovae. It is found that binaries are the main channel that capable of producing type typeIIb supernovae progenitors in the mass range for initial masses below 20 $M_{\odot}$. The formation of type IIb supernova progenitors is extremely sensitive to the initial orbital period. A less massive hydrogen indicates smaller radius and a higher effective temperatures, and vice versa. Binary systems with initial periods between 300 and 720 days produce type IIb progenitors that are a red supergiant. Those with an initial period between 50 and 300 days produce yellow supergiant progenitors and those with initial periods shorter than 50 days, blue supergiant progenitors. Both rapid rotation and larger overshooting can enlarge the carbon-oxygen core mass and lead to higher core temperature and lower central density at the pre-collapse phase. They are also beneficial to surface nitrogen enrichment but restrict the efficiency of the first dredge-up. SN IIb progenitors with low metallicity have smaller hydrogen envelope masses and radii than the high metallicity counterparts. Ultra-stripped binary models have systematically higher core mass fraction $\rm ^{12}C$ left, which has important influence on the compactness of type IIb progenitors.

preprint2021arXiv

Close binary evolution based on Gaia DR2: the origin of late WC-type Wolf-Rayet stars with low luminosity

The observed late-type WC Wolf-Rayet stars (WC7-9) with low luminosity below $\rm \log L/L_{\odot} < 5.4$ in the HR diagram cannot be reproduced satisfactorily by the evolutionary track of single stars. The mass transfer due to Roche lobe overflow drastically modifies the internal structure and surface compositions of two components. Therefore, binaries provide a very promising evolutionary channel to produce these WC stars.

preprint2021arXiv

Is GW190521 the merger of black holes from the first stellar generations?

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.

preprint2020arXiv

How stellar rotation shapes the colour magnitude diagram of the massive intermediate-age star cluster NGC 1846

We present a detailed study of stellar rotation in the massive 1.5 Gyr old cluster NGC 1846 in the Large Magellanic Cloud. Similar to other clusters at this age, NGC 1846 shows an extended main sequence turn-off (eMSTO), and previous photometric studies have suggested it could be bimodal. In this study, we use MUSE integral-field spectroscopy to measure the projected rotational velocities (vsini) of around 1400 stars across the eMSTO and along the upper main sequence of NGC 1846. We measure vsini values up to ~250 km/s and find a clear relation between the vsini of a star and its location across the eMSTO. Closer inspection of the distribution of rotation rates reveals evidence for a bimodal distribution, with the fast rotators centred around vsini = 140 km/s and the slow rotators centred around vsini = 60 km/s. We further observe a lack of fast rotating stars along the photometric binary sequence of NGC 1846, confirming results from the field that suggest that tidal interactions in binary systems can spin down stars. However, we do not detect a significant difference in the binary fractions of the fast and slowly rotating sub-populations. Finally, we report on the serendipitous discovery of a planetary nebula associated with NGC 1846.

preprint2020arXiv

Massive black holes regulated by luminous blue variable mass loss and magnetic fields

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.

preprint2020arXiv

Relative Importance of Convective Uncertainties in Massive Stars

In this work, we investigate the impact of uncertainties due to convective boundary mixing (CBM), commonly called `overshoot', namely the boundary location and the amount of mixing at the convective boundary, on stellar structure and evolution. For this we calculated two grids of stellar evolution models with the MESA code, each with the Ledoux and the Schwarzschild boundary criterion, and vary the amount of CBM. We calculate each grid with the initial masses $15$, $20$ and $25\,\rm{M}_\odot$. We present the stellar structure of the models during the hydrogen and helium burning phases. In the latter, we examine the impact on the nucleosynthesis. We find a broadening of the main-sequence with more CBM, which is more in agreement with observations. Furthermore during the core hydrogen burning phase there is a convergence of the convective boundary location due to CBM. The uncertainties of the intermediate convective zone remove this convergence. The behaviour of this convective zone strongly affects the surface evolution of the model, i.e. how fast it evolves red-wards. The amount of CBM impacts the size of the convective cores and the nucleosynthesis, e.g. the $^{12}$C to $^{16}$O ratio and the weak s-process. Lastly, we determine the uncertainty that the range of parameter values investigated introduce and we find differences of up to $70\%$ for the core masses and the total mass of the star.

preprint2020arXiv

SNAPSHOT: Connections between Internal and Surface Properties of Massive Stars

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.

preprint2019arXiv

The Luminosity Function of Red Supergiants in M31

The mass-loss rates of red supergiant stars (RSGs) are poorly constrained by direct measurements, and yet the subsequent evolution of these stars depends critically on how much mass is lost during the RSG phase. In 2012 the Geneva evolutionary group updated their mass-loss prescription for RSGs with the result that a 20 solar mass star now loses 10x more mass during the RSG phase than in the older models. Thus, higher mass RSGs evolve back through a second yellow supergiant phase rather than exploding as Type II-P supernovae, in accord with recent observations (the so-called "RSG Problem"). Still, even much larger mass-loss rates during the RSG phase cannot be ruled out by direct measurements of their current dust-production rates, as such mass-loss is episodic. Here we test the models by deriving a luminosity function for RSGs in the nearby spiral galaxy M31 which is sensitive to the total mass loss during the RSG phase. We carefully separate RSGs from asymptotic giant branch stars in the color-magnitude diagram following the recent method exploited by Yang and collaborators in their Small Magellanic Cloud studies. Comparing our resulting luminosity function to that predicted by the evolutionary models shows that the new prescription for RSG mass-loss does an excellent job of matching the observations, and we can readily rule out significantly larger values.

preprint2018arXiv

Relative Importance of Convective Uncertainties

Convection plays a key role in the evolution of stars due to energy transport and mixing of composition. Despite its importance, this process is still not well understood. One longstanding conundrum in all 1D stellar evolution codes is the treatment of convective boundaries. In this study we compare two convective uncertainties, the boundary location (Ledoux versus Schwarzschild) and the amount of extra mixing, and their impact on the early evolution of massive stars. With increasing convective boundary mixing (CBM), we find a convergence of the two different boundary locations, a decreasing blue to red super giant ratio and a reduced importance of semiconvection.

preprint2016arXiv

Linking 1D Evolutionary to 3D Hydrodynamical Simulations of Massive Stars

Stellar evolution models of massive stars are important for many areas of astrophysics, for example nucleosynthesis yields, supernova progenitor models and understanding physics under extreme conditions. Turbulence occurs in stars primarily due to nuclear burning at different mass coordinates within the star. The understanding and correct treatment of turbulence and turbulent mixing at convective boundaries in stellar models has been studied for decades but still lacks a definitive solution. This paper presents initial results of a study on convective boundary mixing (CBM) in massive stars. The 'stiffness' of a convective boundary can be quantified using the bulk Richardson number ($\textrm{Ri}_B$), the ratio of the potential energy for restoration of the boundary to the kinetic energy of turbulent eddies. A 'stiff' boundary ($\textrm{Ri}_B \sim 10^4$) will suppress CBM, whereas in the opposite case a 'soft' boundary ($\textrm{Ri}_B \sim 10$) will be more susceptible to CBM. One of the key results obtained so far is that lower convective boundaries (closer to the centre) of nuclear burning shells are 'stiffer' than the corresponding upper boundaries, implying limited CBM at lower shell boundaries. This is in agreement with 3D hydrodynamic simulations carried out by Meakin and Arnett [The Astrophysical Journal 667:448-475, 2007]. This result also has implications for new CBM prescriptions in massive stars as well as for nuclear burning flame front propagation in Super-Asymptotic Giant Branch stars and also the onset of novae.

preprint2016arXiv

Massive Star Evolution: What we do (not) know

The modelling of massive star evolution is a complex task, and is very sensitive to the way physical processes (such as convection, rotation, mass loss, etc.) are included in stellar evolution code. Moreover, the very high observed fraction of binary systems among massive stars makes the comparison with observations difficult. In this paper, we focus on discussing the uncertainties linked to the modelling of convection and rotation in single massive stars.

preprint2016arXiv

On the Effect of Rotation on Populations of Classical Cepheids II. Pulsation Analysis for Metallicities 0.014, 0.006, and 0.002

Classical Cepheid variable stars are high-sensitivity probes of stellar evolution and fundamental tracers of cosmic distances. While rotational mixing significantly affects the evolution of Cepheid progenitors (intermediate-mass stars), the impact of the resulting changes in stellar structure and composition on Cepheids on their pulsational properties is hitherto unknown. Here we present the first detailed pulsational instability analysis of stellar evolution models that include the effects of rotation, for both fundamental mode and first overtone pulsation. We employ Geneva evolution models spanning a three-dimensional grid in mass (1.7 - 15 $M_\odot$), metallicity (Z = 0.014, 0.006, 0.002), and rotation (non-rotating, average & fast rotation). We determine (1) hot and cool instability strip (IS) boundaries taking into account the coupling between convection and pulsation, (2) pulsation periods, and (3) rates of period change. We investigate relations between period and (a) luminosity, (b) age, (c) radius, (d) temperature, (e) rate of period change, (f) mass, (g) the flux-weighted gravity-luminosity relation (FWGLR). We confront all predictions aside from those for age with observations, finding generally excellent agreement. We tabulate period-luminosity relations (PLRs) for several photometric pass-bands and investigate how the finite IS width, different IS crossings, metallicity, and rotation affect PLRs. We show that a Wesenheit index based on H, V, and I photometry is expected to have the smallest intrinsic PLR dispersion. We confirm that rotation resolves the Cepheid mass discrepancy. Period-age relations depend significantly on rotation (rotation increases Cepheid ages), offering a straightforward explanation for evolved stars in binary systems that cannot be matched by conventional isochrones assuming a single age. Finally, we show that Cepheids obey a tight FWGLR.

preprint2016arXiv

The advanced stages of stellar evolution: impact of mass loss, rotation, and link with B[e] stars

In this paper, we discuss some consequences of rotation and mass loss on the evolved stages of massive star evolution. The physical reasons of the time evolution of the surface velocity are explained, and then we show how the late-time evolution of massive stars are impacted in combination with the effects of mass loss. The most interesting result is that in some cases, a massive star can have a blue-red-blue evolution, opening the possibility that Blue Supergiants are composed by two distinct populations of stars: one just leaving the main sequence and crossing the HRD for the first time, and the other one evolving back to the blue side of the HRD after a Red Supergiant phase. We discuss a few possible observational tests that can allow to distinguish these two populations, and how supergiant B[e] stars fit in this context.

preprint2015arXiv

$s$-process production in rotating massive stars at solar and low metallicities

Rotation was shown to have a strong impact on the structure and light element nucleosynthesis in massive stars. In particular, models including rotation can reproduce the primary nitrogen observed in halo extremely metal-poor (EMP) stars. Additional exploratory models showed that rotation may enhance $s$-process production at low metallicity. Here we present a large grid of massive star models including rotation and a full $s$-process network to study the impact of rotation on the weak $s$-process. We explore the possibility of producing significant amounts of elements beyond the strontium peak, which is where the weak $s$-process usually stops. We used the Geneva stellar evolution code coupled to an enlarged reaction network with 737 nuclear species up to bismuth to calculate $15-40\,\text{M}_\odot$ models at four metallicities ($Z = 0.014,10^{-3}$, $10^{-5}$, and $10^{-7}$) from the main sequence up to the end of oxygen burning. We confirm that rotation-induced mixing between the convective H-shell and He-core enables an important production of primary $^{14}$N and $^{22}$Ne and $s$-process at low metallicity. At low metallicity, even though the production is still limited by the initial number of iron seeds, rotation enhances the $s$-process production, even for isotopes heavier than strontium, by increasing the neutron to seed ratio. The increase in this ratio is a direct consequence of the primary production of $^{22}$Ne. Despite nuclear uncertainties affecting the $s$-process production and stellar uncertainties affecting the rotation-induced mixing, our results show a robust production of $s$ process at low metallicity when rotation is taken into account. Considering models with a distribution of initial rotation rates enables to reproduce the observed large range of the [Sr/Ba] ratios in (carbon-enhanced and normal) EMP stars.

preprint2015arXiv

Impact of rotation on stellar models

After a brief recall of the main impacts of stellar rotation on the structure and the evolution of stars, four topics are addressed: 1) the links between magnetic fields and rotation; 2) the impact of rotation on the age determination of clusters; 3) the exchanges of angular momentum between the orbit of a planet and the star due to tides; 4) the impact of rotation on the early chemical evolution of the Milky Way and the origin of the Carbon-Enhanced-Metal-Poor stars.

preprint2015arXiv

Mass loss of red supergiants: a key ingredient for the final evolution of massive stars

Mass-loss rates during the red supergiant phase are very poorly constrained from an observational or theoretical point of view. However, they can be very high, and make a massive star lose a lot of mass during this phase, influencing considerably the final evolution of the star: will it end as a red supergiant? Will it evolve bluewards by removing its hydrogen-rich envelope? In this paper, we briefly summarise the effects of this mass loss and of the related uncertainties, particularly on the population of blue supergiant stars.

preprint2015arXiv

The flux-weighted gravity-luminosity relationship of blue supergiant stars as a constraint for stellar evolution

(abridged) The flux-weighted gravity-luminosity relationship (FGLR) of blue supergiant stars (BSG) links their absolute magnitude to the spectroscopically determined flux-weighted gravity log g = Teff^4. BSG are the brightest stars in the universe at visual light and the application of the FGLR has become a powerful tool to determine extragalactic distances. Observationally, the FGLR is a tight relationship with only small scatter. It is, therefore, ideal to be used as a constraint for stellar evolution models. The goal of this work is to investigate whether stellar evolution can reproduce the observed FGLR and to develop an improved foundation of the FGLR as an extragalactic distance indicator. We use different grids of stellar models for initial masses between 9 and 40 Msun, for metallicities between Z = 0.002 and 0.014, with and without rotation, computed with various mass loss rates during the red supergiant phase. For each of these models we discuss the details of post-main sequence evolution and construct theoretical FGLRs by means of population synthesis models which we then compare with the observed FGLR. In general, the stellar evolution model FGLRs agree reasonably well with the observed one. There are, however, differences between the models, in particular with regard to the shape and width (scatter) in the flux-weighted gravity-luminosity plane. The best agreement is obtained with models which include the effects of rotation and assume that the large majority, if not all the observed BSG evolve towards the red supergiant phase and only a few are evolving back from this stage. The effects of metallicity on the shape and scatter of the FGLR are small. (...)

preprint2015arXiv

Wolf-Rayet stars as an evolved stage of stellar life

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.

preprint2014arXiv

Code dependencies of pre-supernova evolution and nucleosynthesis in massive stars: Evolution to the end of core helium burning

Massive stars are key sources of radiative, kinetic, and chemical feedback in the universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes---GENEC, KEPLER and MESA---to compute non-rotating stellar models of $15~\mathrm{M}_\odot$, $20~\mathrm{M}_\odot$, and $25~\mathrm{M}_\odot$ and compare their nucleosynthesis. We follow the evolution from the main sequence until the end of core helium burning. The GENEC and KEPLER models hold physics assumptions used in large grids of published models. The MESA code was set up to use convective core overshooting such that the CO core masses are consistent with those obtained by GENEC. For all models, full nucleosynthesis is computed using the NuGrid post-processing tool MPPNP. We find that the surface abundances predicted by the models are in reasonable agreement. In the helium core, the standard deviation of the elemental overproduction factors for Fe to Mo is less than $30\,\%$---smaller than the impact of the present nuclear physics uncertainties. For our three initial masses, the three stellar evolution codes yield consistent results. Differences in key properties of the models, e.g., helium and CO core masses and the time spent as a red supergiant, are traced back to the treatment of convection and, to a lesser extent, mass loss. The mixing processes in stars remain the key uncertainty in stellar modelling. Better constrained prescriptions are thus necessary to improve the predictive power of stellar evolution models.

preprint2014arXiv

Evolution of surface CNO abundances in massive stars

The nitrogen to carbon (N/C) and nitrogen to oxygen (N/O) ratios are the most sensitive quantities to mixing in stellar interiors of intermediate and massive stars. We further investigate the theoretical properties of these ratios as well as put in context recent observational results obtained by the VLT-FLAMES Survey of massive stars in the Galaxy and the Magellanic Clouds. We consider analytical relations and numerical models of stellar evolution as well as our own stellar atmosphere models, and we critically re-investigate observed spectra. On the theoretical side, the N/C vs N/O plot shows little dependence on the initial stellar masses, rotation velocities, and nature of the mixing processes up to relative enrichment of N/O by a factor of about four, thus this plot constitutes an ideal quality test for observational results. The comparison between the FLAMES Survey and theoretical values shows overall agreement, despite the observational scatter of the published results. The existence of some mixing of CNO products is clearly confirmed, however the accuracy of the data is not sufficient for allowing a test of the significant differences between different models of rotating stars and the Geneva models. We discuss reasons (for the most part due to observational bias) why part of the observational data points should not be considered for this comparison. When these observational data points are not considered, the scatter is reduced. Finally, the N/C vs N/O plot potentially offers a powerful way for discriminating blue supergiants before the red supergiant stage from those after it. Also, red supergiants of similar low velocities may exhibit different N enrichments, depending on their initial rotation during the main-sequence phase.

preprint2014arXiv

Linking 1D Stellar Evolution to 3D Hydrodynamical Simulations

In this contribution we present initial results of a study on convective boundary mixing (CBM) in massive stellar models using the GENEVA stellar evolution code. Before undertaking costly 3D hydrodynamic simulations, it is important to study the general properties of convective boundaries, such as the: composition jump; pressure gradient; and `stiffness'. Models for a 15Mo star were computed. We found that for convective shells above the core, the lower (in radius or mass) boundaries are `stiffer' according to the bulk Richardson number than the relative upper (Schwarzschild) boundaries. Thus, we expect reduced CBM at the lower boundaries in comparison to the upper. This has implications on flame front propagation and the onset of novae.

preprint2014arXiv

Populations of rotating stars III. SYCLIST, the new Geneva Population Synthesis code

(abridged) We present a new tool for building synthetic colour-magnitude diagrams of coeval stellar populations. We study, from a theoretical point of view, the impact of axial rotation of stars on various observed properties of single-aged stellar populations: magnitude at the turnoff, photometric properties of evolved stars, surface velocities, surface abundances, and the impact of rotation on the age determination of clusters by an isochrone fitting. Stellar models for different initial masses, metallicities, and zero-age main sequence (ZAMS) rotational velocities are used for building interpolated stellar tracks, isochrones, and synthetic clusters for various ages and metallicities. The synthetic populations account for the effects of the initial distribution of the rotational velocities on the ZAMS, the impact of the inclination angle and the effects of gravity and limb darkening, unresolved binaries and photometric errors. Interpolated tracks, isochrones, and synthetic clusters can be computed through a public web interface. For clusters with a metallicity in the range [0.002, 0.014] and an age between 30 Myr and 1 Gyr, the fraction of fast rotators on the main sequence (MS) band is the largest just below the turnoff. This is a natural consequence of the increase in the MS lifetime due to rotation. The fraction of fast rotators one magnitude below the turnoff also increases with the age of the cluster between 30 Myr and 1 Gyr. There is an increase in the fraction of enriched stars when the metallicity decreases. We show that the use of isochrones computed from rotating stellar models with an initial rotation that is representative of the average initial rotation of the stars in clusters provides a reasonable estimate of the age, even though stars in a real cluster did not start their evolution with an identical initial rotation.

preprint2014arXiv

Rotation and the Cepheid Mass Discrepancy

We recently showed that rotation significantly affects most observable Cepheid quantities, and that rotation, in combination with the evolutionary status of the star, can resolve the long-standing Cepheid mass discrepancy problem. We therefore provide a brief overview of our results regarding the problem of Cepheid masses. We also briefly mention the impact of rotation on the Cepheid period-luminosity(-color) relation, which is crucial for determining extragalactic distances, and thus for calibrating the Hubble constant.

preprint2014arXiv

Stability boundaries for massive stars in the sHR diagram

Stability boundaries of radial pulsations in massive stars are compared with positions of variable and non-variable blue-supergiants in the spectroscopic HR (sHR) diagram (Langer & Kudritzki 2014), whose vertical axis is $4\log T_{\rm eff}-\log g \,(= \log L/M)$. Observational data indicate that variables tend to have higher $L/M$ than non-variables in agreement with the theoretical prediction. However, many variable blue-supergiants are found to have values of $L/M$ below the theoretical stability boundary; i.e., surface gravities seem to be too high by around 0.2-0.3 dex.

preprint2014arXiv

The evolution of massive stars and their spectra I. A non-rotating 60 Msun star from the zero-age main sequence to the pre-supernova stage

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]

preprint2013arXiv

Circum-stellar medium around rotating massive stars at solar metallicity

Aims. Observations show nebulae around some massive stars but not around others. If observed, their chemical composition is far from homogeneous. Our goal is to put these observational features into the context of the evolution of massive stars and their circumstellar medium (CSM) and, more generally, to quantify the role of massive stars for the chemical and dynamical evolution of the ISM. Methods. Using the A-MAZE code, we perform 2d-axisymmetric hydrodynamical simulations of the evolution of the CSM, shaped by stellar winds, for a whole grid of massive stellar models from 15 to 120 Msun and following the stellar evolution from the zero-age main-sequence to the time of supernova explosion. In addition to the usual quantities, we also follow five chemical species: H, He, C, N, and O. Results. We show how various quantities evolve as a function of time: size of the bubble, position of the wind termination shock, chemical composition of the bubble, etc. The chemical composition of the bubble changes considerably compared to the initial composition, particularly during the red-supergiant (RSG) and Wolf-Rayet (WR) phases. In some extreme cases, the inner region of the bubble can be completely depleted in hydrogen and nitrogen, and is mainly composed of carbon, helium and oxygen. We argue why the bubble typically expands at a lower rate than predicted by self-similarity theory. In particular, the size of the bubble is very sensitive to the density of the ISM, decreasing by a factor of around 2.5 for each additional dex in ISM density. The bubble size also decreases with the metallicity of the central star, as low-metallicity stars have weaker winds. Our models qualitatively fit the observations of WR ejecta nebulae.

preprint2013arXiv

Evolution and fate of very massive stars

There is observational evidence that supports the existence of Very Massive Stars in the local universe. First, very massive stars (Mini<=320 M) have been observed in the Large Magellanic Cloud . Second, there are observed SNe that bear the characteristics of Pair Creation Supernovae which have very massive stars as progenitors. The most promising candidate to date is SN2007bi. In order to investigate the evolution and fate of nearby very massive stars, we calculated a new grid of models for such objects, for solar, LMC and SMC metallicities, which covers the initial mass range from 120 to 500M. Both rotating and non-rotating models were calculated using the Geneva stellar evolution code and evolved until at least the end of helium burning and for most models until oxygen burning. Since very massive stars have very large convective cores during the Main-Sequence phase, their evolution is not so much affected by rotational mixing, but more by mass loss through stellar winds. Their evolution is never far from a homogeneous evolution even without rotational mixing. All the VMS, at all the metallicities studied here, end their life as WC(WO) type stars. At solar metallicity, none of our models is expected to explode as a PCSN. At the metallicity of the LMC, only stars more massive than 300 M are expected to explode as PCSNe. At the SMC metallicity, the mass range for the PCSN progenitors is much larger and comprises stars with initial masses between about 100 and 290 M . All VMS stars in the metallicity range studied here produce either a type Ib or a type Ic SN. We estimate that the progenitor of SN2007bi, assuming a SMC metallicity, had an initial mass between 160 and 175 M . None of models presented in this grid produce GRBs or magnetars. They lose too much angular momentum by mass loss or avoid the formation of a BH by producing a completely disruptive PCSN.

preprint2013arXiv

Evolution of blue supergiants and αCygni variables; Puzzling CNO surface abundances

A massive star can enter the blue supergiant region either evolving directly from the main-sequence, or evolving from a previous red supergiant stage. The fractions of the blue supergiants having different histories depend on the internal mixing and mass-loss during the red supergiant stage. We study the possibility to use diagnostics based on stellar pulsation to discriminate blue supergiants having different evolution histories. For this purpose we have studied the pulsation property of massive star models calculated with the Geneva stellar evolution code for initial masses ranging from 8 to 50 M$_\odot$ with a solar metallicity of $Z=0.014$. We have found that radial pulsations are excited in the blue-supergiant region only in the models that had been red-supergiants before. This would provide us with a useful mean to diagnose the history of evolution of each blue-supergiant. At a given effective temperature, much more nonradial pulsations are excited in the model after the red-supergiant stage than in the model evolving towards the red-supergiant. The properties of radial and nonradial pulsations in blue supergiants are discussed. Predicted periods are compared with period ranges observed in some α-Cygni variables in the Galaxy and NGC 300. We have found that blue supergiant models after the red-supergiant stage roughly agree with observed period ranges in most cases. However, we are left with the puzzle that the predicted surface N/C and N/O ratios seem to be too high compared with those of Deneb and Rigel.

preprint2013arXiv

Four open questions in massive star evolution

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?

preprint2013arXiv

Fundamental properties of core-collapse Supernova and GRB progenitors: predicting the look of massive stars before death

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]

preprint2013arXiv

Grids of stellar models with rotation - III. Models from 0.8 to 120 Msun at a metallicity Z = 0.002

(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.

preprint2013arXiv

How the mass-loss rates of red-supergiants determine the fate of massive stars ?

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.

preprint2013arXiv

Populations of rotating stars. - I. Models from 1.7 to 15 Msun at Z = 0.014, 0.006, and 0.002 with Ω/Ωcrit between 0 and 1

B-type stars are known to rotate at various velocities, including very fast rotators near the critical velocity as the Be stars. In this paper, we provide stellar models covering the mass range between 1.7 to 15 Msun, which includes the typical mass of known Be stars, at Z = 0.014, 0.006, and 0.002 and for an extended range of initial velocities on the zero-age main sequence. We used the Geneva stellar-evolution code, including the effects of shellular rotation, with a numerical treatment that has been improved so the code can precisely track the variation in the angular momentum content of the star as it changes under the influence of radiative winds and/or mechanical mass loss. We discuss the impact of the initial rotation rate on the tracks in the Hertzsprung-Russell diagram, the main-sequence (MS) lifetimes, the evolution of the surface rotation and abundances, as well as on the ejected masses of various isotopes. Among the new results obtained from the present grid we find that 1) fast-rotating stars with initial masses around 1.7 Msun present at the beginning of the core hydrogen-burning phase quite small convective cores with respect to their slowly rotating counterparts. This fact may be interesting to keep in mind in the framework of the asteroseismic studies of such stars. 2) The contrast between the core and surface angular velocity is higher in slower rotating stars. The values presently obtained are in agreement with the very few values obtained for B-type stars from asteroseismology. 3) At Z = 0.002, the stars in the mass range of 1.7 to 3 Msun with a mean velocity on the MS of the order of 150 km/s show N/H enhancement superior to 0.2 dex at mid-MS, and superior to 0.4 dex at the end of the MS phase. At solar metallicity the corresponding values are below 0.2 dex at any time in the MS.

preprint2013arXiv

Progenitors of supernova Ibc: a single Wolf-Rayet star as the possible progenitor of the SN Ib iPTF13bvn

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.

preprint2013arXiv

Red supergiants and stellar evolution

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.

preprint2013arXiv

Strange mode instability for micro-variations in Luminous Blue Variables

If a massive star has lost significant mass during its red-supergiant stage, it would return to blue region in the HR diagram and spend a part of the core-He burning stage as a blue supergiant having a luminosity to mass ratio (L/M) considerably larger than about 10^4 (in solar units); the duration depends on the degree of internal mixing and on the metallicity. Then, various stellar pulsations are excited by enhanced κ-mechanism and strange mode instability. Assuming these pulsations to be responsible for (at least some of) the quasi-periodic light and radial-velocity variations in αCygni variables including luminous blue variables (LBVs; or S Dor variables), we can predict masses and surface compositions for these variables, and compare them with observed ones to constrain the evolutionary models. We discuss radial pulsations excited in evolutionary models of an initial mass of 40 M_\odot with solar metallicity of Z=0.014, and compare them to micro-variations in the two Galactic LBVs, HR Car and HD 160529. We have found that these stars should has lost more than half of the initial mass and their surface CNO abundances should be significantly modified from the original ones showing partial H-burning products.

preprint2013arXiv

The past and future evolution of a star like Betelgeuse

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.

preprint2013arXiv

The puzzle of the CNO abundances of α Cygni variables resolved by the Ledoux criterion?

Recent stellar evolution computations show that the blue supergiant (BSG) stars could come from two distinct populations: a first group arising from massive stars that just left the main sequence (MS) and are crossing the Hertzsprung-Russell diagram (HRD) towards the red supergiant (RSG) branch, and a second group coming from stars that have lost considerable amount of mass during the RSG stage and are crossing the HRD for a second time towards the blue region. Due to very different luminosity-to- mass ratio, only stars from the second group are expected to have excited pulsations observable at the surface. In a previous work, we have shown that our models were able to reproduce the pulsational properties of BSGs. However, these models failed to reproduce the surface chemical composition of stars evolving back from a RSG phase. In this paper, we show how the use of the Ledoux criterion instead of the Schwarzschild one for convection allows to significantly improve the agreement with the observed chemical composition, while keeping the agreement with the pulsation periods. This gives some support to the Ledoux criterion.

preprint2012arXiv

Grids of stellar models with rotation II. WR populations and supernovae/GRB progenitors at Z = 0.014

We used a recent grid of stellar models computed with and without rotation to make predictions concerning the WR populations and the frequency of different types of core-collapse SNe. Current rotating models were checked to provide good fits to the following features: solar luminosity and radius at the solar age, main-sequence width, red-giant and red-supergiant (RSG) positions in the HRD, surface abundances, and rotational velocities. Rotating stellar models predict that about half of the observed WR stars and at least half of the type Ibc SNe may be produced through the single-star evolution channel. Rotation increases the duration of the WNL and WNC phases, while reducing those of the WNE and WC phases, as was already shown in previous works. Rotation increases the frequency of type Ic SNe. The upper mass limit for type II-P SNe is \sim 19.0 MSun for the non rotating models and \sim 16.8 MSun for the rotating ones. Both values agree with observations. Moreover, present rotating models provide a very good fit to the progenitor of SN 2008ax. We discuss future directions of research for further improving the agreement between the models and the observations. We conclude that the mass-loss rates in the WNL and RSG phases are probably underestimated at present. We show that up to an initial mass of 40 M\odot, a surface magnetic field inferior to about 200 G may be sufficient to produce some braking. Much lower values are needed at the red supergiant stage. We suggest that the presence/absence of any magnetic braking effect may play a key role in questions regarding rotation rates of young pulsars and the evolution leading to LGRBs.

preprint2012arXiv

On the Eddington limit and WR Stars

We examine some properties of stars evolving close to the classical Eddington limit for electron-scattering opacity, when these stars maintain a chemically homogeneous structure as a result of mixing and/or mass loss. We consider analytical relations and models computed with the Geneva code. Homologous, chemically homogeneous stars evolving with a constant Eddington factor obey a relation of the form mu^2 M = const. This applies, for example, to Wolf-Rayet (WR) stars in stages without hydrogen. The value of the constant may depend on the metallicity, initial mass, evolutionary stage, and physical processes included in the considered homologous evolutionary sequence. An average value of the constant between 20 and 40 in solar units is consistent with the masses of Galactic WR stars.

preprint2012arXiv

The Wolf-Rayet Content of M31

Wolf-Rayet stars are evolved massive stars, and the relative number of WC-type and WN-type WRs should vary with the metallicity of the host galaxy, providing a sensitive test of stellar evolutionary theory. However, past studies of the WR content of M31 have been biased towards detecting WC stars, as their emission line signatures are much stronger than those of WNs. Here we present the results of a survey covering all of M31's optical disk (2.2 deg^2), with sufficient sensitivity to detect the weaker-lined WN-types. We identify 107 newly found WR stars, mostly of WN-type. This brings the total number of spectroscopically confirmed WRs in M31 to 154, a number we argue is complete to ~95%, except in regions of unusually high reddening. This number is consistent with what we expect from the integrated Halpha luminosity compared to that of M33. The majority of these WRs formed in OB associations around the Population I ring, although 5% are truly isolated. Both the relative number of WC to WN-type stars as well as the WC subtype distribution suggest that most WRs exist in environments with higher-than-solar metallicities, which is consistent with studies of M31's metallicity. Although the WC-to-WN ratio we find for M31 is much lower than that found by previous studies, it is still higher than what the Geneva evolutionary models predict. This may suggest that Roche-lobe overflow produces the excess of WC stars observed at high metallicity, or that the assumed rotational velocities in the models are too high.

preprint2011arXiv

Grids of stellar models with rotation - I. Models from 0.8 to 120 Msun at solar metallicity (Z = 0.014)

[abridged] Many topical astrophysical research areas, such as the properties of planet host stars, the nature of the progenitors of different types of supernovae and gamma ray bursts, and the evolution of galaxies, require complete and homogeneous sets of stellar models at different metallicities in order to be studied during the whole of cosmic history. We present here a first set of models for solar metallicity, where the effects of rotation are accounted for in a homogeneous way. We computed a grid of 48 different stellar evolutionary tracks, both rotating and non-rotating, at Z=0.014, spanning a wide mass range from 0.8 to 120 Msun. For each of the stellar masses considered, electronic tables provide data for 400 stages along the evolutionary track and at each stage, a set of 43 physical data are given. These grids thus provide an extensive and detailed data basis for comparisons with the observations. The rotating models start on the ZAMS with a rotation rate Vini/Vcrit=0.4. The evolution is computed until the end of the central carbon-burning phase, the early AGB phase, or the core helium-flash for, respectively, the massive, intermediate, and both low and very low mass stars. The initial abundances are those deduced by Asplund and collaborators, which best fit the observed abundances of massive stars in the solar neighbourhood. We update both the opacities and nuclear reaction rates, and introduce new prescriptions for the mass-loss rates as stars approach the Eddington and/or the critical velocity. We account for both atomic diffusion and magnetic braking in our low-mass star models. [...]

preprint2011arXiv

Red Supergiants, Luminous Blue Variables and Wolf-Rayet stars: the single massive star perspective

We discuss, in the context of the single star scenario, the nature of the progenitors of Red Supergiants (RSG), of Luminous Blue Variables (LBV) and of Wolf-Rayet (WR) stars. These three different populations correspond to evolved phases of Main Sequence (MS) OB stars. Axial rotation and mass loss have a great influence on massive star evolution in general and more specifically on the durations of these different phases. Moderate rotation and mass loss, during the MS phase, favor the evolution towards the RSG stage. Fast rotation and strong mass loss during the MS phase, in contrast, prevent the star from becoming a RSG and allow the star to pass directly from the OB star phase into the WR phase. Mass loss during the RSG stage may make the star evolve back in the blue part of the HR diagram. We argue that such an evolution may be more common than presently accounted for in stellar models. This might be the reason for the lack of type IIP SNe with RSG progenitors having initial masses between 18 and 30 M$_\odot$. The LBVs do appear as a possible transition phase between O and WR stars or between WNL and WNE stars. Fast rotation and/or strong mass loss during the Main-Sequence phase prevent the formation of LBV stars. The mechanisms driving the very strong ejections shown by LBV stars are still unknown. (abridged)

preprint2011arXiv

Simulating the circum-stellar environment of supernova and GRB progenitors by combining stellar evolution models and hydrodynamical code

The medium around massive stars is strongly shaped by the stellar winds. Those winds depend on various stellar parameters (effective temperature, luminosity, chemical composition, rotation, ...), which are varying as a function of the time. Using the wind properties obtained by classical stellar evolution code allows for the multi- D hydrodynamical simulation of the circum-stellar medium accounting for the time variations of the wind. We present here the preliminary results of the simulation of the medium around a fast rotating star. Comparing the results of such simulations with the observed properties of the circum-stellar medium will allow to better understand the interactions of the stellar winds with the interstellar medium, and could give hints on the past history of the mass loss around massive stars, providing constraints for stellar models.

preprint2011arXiv

Yellow supergiants as supernova progenitors: an indication of strong mass loss for red supergiants?

The increasing observed number of supernova events allows for finding ever more frequently the progenitor star in archive images. In a few cases, the progenitor star is a yellow supergiant star. The estimated position in the Hertzsprung-Russell diagram of these stars is not compatible with the theoretical tracks of classical single star models. According to several authors, the mass-loss rates during the red supergiant phase could be underestimated. We study the impact of an increase of these mass-loss rates on the position of 12 to 15 M\odot stars at the end of their nuclear life, in order to reconcile the theoretical tracks with the observed yellow supergiant progenitors. We perform calculations of 12 to 15 M\odot rotating stellar models using the Geneva stellar evolution code. To account for the uncertainties in the mass-loss rates during the RSG phase, we increase the mass-loss rate of the star (between 3 and 10 times the standard one) during that phase and compare the evolution of stars undergoing such high mass-loss rates with models computed with the standard mass-loss prescription. We show that the final position of the models in the Hertzsprung-Russell diagram depends on the mass loss they undergo during the red supergiant phase. With an increased mass-loss rate, we find that some models end their nuclear life at positions that are compatible with the observed position of several supernova progenitors. We conclude that an increased mass-loss rate (whom physical mechanism still need to be clarified) allows single star models to reproduce simultaneously the estimated position in the HRD of the YSG SN progenitors, as well as the SN type.

preprint2010arXiv

Are C-rich ultra iron poor stars also He-rich?

(abridged) The three most iron poor stars presently known ([Fe/H] equal to -5.96, -5.4 and -4.75) are carbon-rich, they are called C-Rich Ultra-Metal Poor Stars (CRUMPS). The origin of their peculiar surface abundances is not understood. We propose a synthetic view of the different models so far proposed to explain the peculiar abundances observed at the surface of the CRUMP stars. We deduce some expected trends based on nucleosynthetic arguments and look for signatures allowing to discriminate among models. We discuss the conditions for having CRUMP stars which are He-rich, i.e. with a mass fraction of helium greater than 0.30 and up to 0.60. We discuss the chemical composition of stars made of interstellar medium mixed with wind material of very metal poor massive stars, with wind plus supernova ejecta and with material extracted from the envelope of early AGB stars. Rotating and non-rotating models are considered. CRUMP stars formed from wind material of massive stars mixed with small amounts of pristine interstellar medium are He-rich (helium mass fraction between 0.30 and 0.60), Li-depleted and present low $^{12}$C/$^{13}$C ratios (inferior to 10 in number). Such He-rich stars, if discovered, would confirm that the most metal poor CRUMPs formed from essentially pure wind/envelope material. They would provide the most direct way to probe the nucleosynthetic outputs of the first generations of stars. We show that rotation is a key ingredient to explain the abundance patterns of CRUMPS stars.

preprint2010arXiv

Constraints on rotational mixing from surface evolution of light elements in massive stars

Context. Light elements and nitrogen surface abundances together can constrain the mixing efficiencies in massive stars on the main sequence, because moderate mixing at the surface leads to a depletion of light elements but only later to an enrichment in nitrogen. Aims. We want to test the rotational mixing prescriptions included in the Geneva stellar evolution code (GENEC) by following the evolution of surface abundances of light isotopes in massive stars. Methods. The GENEC is a 1D code containing sophisticated prescriptions for rotational mixing. We implemented an extended reaction network into this code including the light elements Li, Be and B, which allowed us to perform calculations testing the rotation induced mixing. Results. We followed 9, 12 and 15 solar mass models with rotation from the zero age main sequence up to the end of He burning. The calculations show the expected behaviour with faster depletion of light isotopes for faster rotating stars and more massive stars. Conclusions. We find that the mixing prescriptions used in the present rotating models for massive single stars can account for most of the observations. However the uncertainties are quite large making it hard to draw a firm conclusion on the mixing scenario.

preprint2010arXiv

Effects of anisotropic winds on massive stars evolution

Whenever stars are rotating very fast (Omega/Omega_crit > 0.7, with Omega_crit the Keplerian angular velocity of the star accounting for its deformation) radiative stellar winds are enhanced in polar regions. This theoretical prediction is now confirmed by interferometric observations of fast rotating stars.} Polar winds remove less angular momentum than spherical winds and thus allow the star to keep more angular momentum. We quantitatively assess the importance of this effect. First we use a semi-analytical approach to estimate the variation of the angular momentum loss when the rotation parameter increases. Then we compute complete 9 MSun stellar models at very high angular velocities (starting on the ZAMS with Omega/Omega_crit = 0.8 and reaching the critical velocity during the Main Sequence) with and without radiative wind anisotropies. When wind anisotropies are accounted for, the angular momentum loss rate is reduced by less than 4% for Omega/Omega_crit < 0.9 with respect to the case of spherical winds. The reduction amounts to at most 30% when the star is rotating near the critical velocity. These values result from two counteracting effects: on the one hand polar winds reduce the loss of angular momentum, on the other hand, surface deformations imply that the mass which is lost at high co-latitude is lost at a larger distance from the rotational axis and thus removes more angular momentum. In contrast with previous studies, which neglected surface deformations, we show that the radiative wind anisotropies have a relatively modest effect on the evolution of the angular momentum content of fast rotating stars.

preprint2010arXiv

Equatorial mass loss from Be stars

Be stars are thought to be fast rotating stars surrounded by an equatorial disc. The formation, structure and evolution of the disc are still not well understood. In the frame of single star models, it is expected that the surface of an initially fast rotating star can reach its keplerian velocity (critical velocity). The Geneva stellar evolution code has been recently improved, in order to obtain some estimates of the total mass loss and of the mechanical mass loss rates in the equatorial disc during the whole critical rotation phase. We present here the first results of the computation of a grid of fast rotating B stars evolving towards the Be phase, and discuss the first estimates we obtained.

preprint2010arXiv

Massive stellar models: rotational evolution, metallicity effects

The Be star phenomenon is related to fast rotation, although the cause of this fast rotation is not yet clearly established. The basic effects of fast rotation on the stellar structure are reviewed: oblateness, mixing, anisotropic winds. The processes governing the evolution of the equatorial velocity of a single star (transport mechanisms and mass loss) are presented, as well as their metallicity dependence. The theoretical results are compared to observations of B and Be stars in the Galaxy and the Magellanic Clouds.

preprint2009arXiv

Evolution of Massive Stars along the Cosmic History

Massive stars are "cosmic engines" (cf the title of the IAU Symposium 250). They drive the photometric and chemical evolution of galaxies, inject energy and momentum through stellar winds and supernova explosions, they modify in this way the physical state of the interstellar gas and have an impact on star formation. The evolution of massive stars depends sensitively on the metallicity which has an impact on the intensity of the line driven stellar winds and on rotational mixing. We can distinguish four metallicity regimes: 1.- the Pop III regime $0 \le Z < \sim 10^{-10}$; 2.- The low metallicity regime $10^{-10} \le Z < 0.001$; 3.- The near solar metallicity regime $0.001 \le Z < 0.020$; 4.- The high metallicity regime $0.020 \le Z$. In each of these metallicity ranges, some specific physical processes occur. In this review we shall discuss these physical processes and their consequences for nucleosynthesis and the massive star populations in galaxies. We shall mainly focus on the effects of axial rotation and mass loss by line driven winds, although of course other processes like binarity, magnetic fields, transport processes by internal waves may also play important roles.

preprint2009arXiv

Nucleosynthesis in Rotating massive stars and Abundances in the Early Galaxy

We discuss three effects of axial rotation at low metallicity. The first one is the mixing of the chemical species which is predicted to be more efficient in low metallicity environments. A consequence is the production of important quantities of primary $^{14}$N, $^{13}$C, $^{22}$Ne and a strong impact on the nucleosynthesis of the {\it s}-process elements. The second effect is a consequence of the first. Strong mixing makes possible the apparition at the surface of important quantities of CNO elements. This increases the opacity of the outer layers and may trigger important mass loss by line driven winds. The third effect is the fact that, during the main-sequence phase, stars, at very low metallicity, reach more easily than their more metal rich counterparts, the critical velocity\footnote{The critical velocity is the surface equatorial velocity such that the centrifugal acceleration compensates for the local gravity.}. We discuss the respective importance of these three effects as a function of the metallicity. We show the consequences for the early chemical evolution of the galactic halo and for explaining the CEMP stars. We conclude that rotation is probably a key feature which contributes in an important way to shape the evolution of the first stellar generations in the Universe.

preprint2008arXiv

Evolution and chemical and dynamical effects of high-mass stars

We review general characteristics of massive stars, present the main observable constraints that stellar models should reproduce. We discuss the impact of massive star nucleosynthesis on the early phases of the chemical evolution of the Milky Way (MW). We show that rotating models can account for the important primary nitrogen production needed at low metallicity. Interestingly such rotating models can also better account for other features as the variation with the metallicity of the C/O ratio. Damped Lyman Alpha (DLA) systems present similar characteristics as the halo of the MW for what concern the N/O and C/O ratios. Although in DLAs, the star formation history might be quite different from that of the halo, in these systems also, rotating stars (both massive and intermediate) probably play an important role for explaining these features. The production of primary nitrogen is accompanied by an overproduction of other elements as $^{13}$C, $^{22}$Ne and s-process elements. We show also how the observed variation with the metallicity of the number ratio of type Ibc to type II supernovae may be a consequence of the metallicity dependence of the line-driven stellar winds.