We present a dense model grid with tailored input chemical composition appropriate for the Large Magellanic Cloud. We use a one-dimensional hydrodynamic stellar evolution code, which accounts for rotation, transport of angular momentum by magnetic fi
elds, and stellar wind mass loss to compute our detailed models. We calculate stellar evolution models with initial masses of 70-500 Msun and with initial surface rotational velocities of 0-550 km/s, covering the core-hydrogen burning phase of evolution. We find our rapid rotators to be strongly influenced by rotationally induced mixing of helium, with quasi-chemically homogeneous evolution occurring for the fastest rotating models. Above 160 Msun, homogeneous evolution is also established through mass loss, producing pure helium stars at core hydrogen exhaustion independent of the initial rotation rate. Surface nitrogen enrichment is also found for slower rotators, even for stars that lose only a small fraction of their initial mass. For models above 150 MZAMS, and for models in the whole considered mass range later on, we find a considerable envelope inflation due to the proximity of these models to their Eddington limit. This leads to a maximum zero-age main sequence surface temperature of 56000 K, at 180 Msun, and to an evolution of stars in the mass range 50-100 Msun to the regime of luminous blue variables in the HR diagram with high internal Eddington factors. Inflation also leads to decreasing surface temperatures during the chemically homogeneous evolution of stars above 180 Msun. The cool surface temperatures due to the envelope inflation in our models lead to an enhanced mass loss, which prevents stars at LMC metallicity from evolving into pair-instability supernovae. The corresponding spin-down will also prevent very massive LMC stars to produce long-duration gamma-ray bursts, which might, however, originate from lower masses.
The distribution of stars in the Hertzsprung-Russell diagram narrates their evolutionary history and directly assesses their properties. Placing stars in this diagram however requires the knowledge of their distances and interstellar extinctions, whi
ch are often poorly known for Galactic stars. The spectroscopic Hertzsprung-Russell diagram (sHRD) tells similar evolutionary tales, but is independent of distance and extinction measurements. Based on spectroscopically derived effective temperatures and gravities of almost 600 stars, we derive for the first time the observational distribution of Galactic massive stars in the sHRD. While biases and statistical limitations in the data prevent detailed quantitative conclusions at this time, we see several clear qualitative trends. By comparing the observational sHRD with different state-of-the-art stellar evolutionary predictions, we conclude that convective core overshooting may be mass-dependent and, at high mass ($geq 15,M_odot$), stronger than previously thought. Furthermore, we find evidence for an empirical upper limit in the sHRD for stars with $T_{rm{eff}}$ between 10000 and 32000 K and, a strikingly large number of objects below this line. This over-density may be due to inflation expanding envelopes in massive main-sequence stars near the Eddington limit.
At least 5 per cent of the massive stars are moving supersonically through the interstellar medium (ISM) and are expected to produce a stellar wind bow shock. We explore how the mass loss and space velocity of massive runaway stars affect the morphol
ogy of their bow shocks. We run two-dimensional axisymmetric hydrodynamical simulations following the evolution of the circumstellar medium of these stars in the Galactic plane from the main sequence to the red supergiant phase. We find that thermal conduction is an important process governing the shape, size and structure of the bow shocks around hot stars, and that they have an optical luminosity mainly produced by forbidden lines, e.g. [OIII]. The Ha emission of the bow shocks around hot stars originates from near their contact discontinuity. The H$alpha$ emission of bow shocks around cool stars originates from their forward shock, and is too faint to be observed for the bow shocks that we simulate. The emission of optically-thin radiation mainly comes from the shocked ISM material. All bow shock models are brighter in the infrared, i.e. the infrared is the most appropriate waveband to search for bow shocks. Our study suggests that the infrared emission comes from near the contact discontinuity for bow shocks of hot stars and from the inner region of shocked wind for bow shocks around cool stars. We predict that, in the Galactic plane, the brightest, i.e. the most easily detectable bow shocks are produced by high-mass stars moving with small space velocities.
Recent discoveries of weak and fast optical transients raise the question of their origin. We investigate the minimum ejecta mass associated with core-collapse supernovae (SNe) of Type Ic. We show that mass transfer from a helium star to a compact co
mpanion can produce an ultra-stripped core which undergoes iron core collapse and leads to an extremely fast and faint SN Ic. In this Letter, a detailed example is presented in which the pre-SN stellar mass is barely above the Chandrasekhar limit, resulting in the ejection of only ~0.05-0.20 M_sun of material and the formation of a low-mass neutron star. We compute synthetic light curves of this case and demonstrate that SN 2005ek could be explained by our model. We estimate that the fraction of such ultra-stripped to all SNe could be as high as 0.001-0.01. Finally, we argue that the second explosion in some double neutron star systems (for example, the double pulsar PSR J0737-3039B) was likely associated with an ultra-stripped SN Ic.
Aims: Projected rotational velocities (vsini) have been estimated for 334 targets in the VLT-FLAMES Tarantula survey that do not manifest significant radial velocity variations and are not supergiants. They have spectral types from approximately O9.5
to B3. The estimates have been analysed to infer the underlying rotational velocity distribution, which is critical for understanding the evolution of massive stars. Methods: Projected rotational velocities were deduced from the Fourier transforms of spectral lines, with upper limits also being obtained from profile fitting. For the narrower lined stars, metal and non-diffuse helium lines were adopted, and for the broader lined stars, both non-diffuse and diffuse helium lines; the estimates obtained using the different sets of lines are in good agreement. The uncertainty in the mean estimates is typically 4% for most targets. The iterative deconvolution procedure of Lucy has been used to deduce the probability density distribution of the rotational velocities. Results: Projected rotational velocities range up to approximately 450 kms and show a bi-modal structure. This is also present in the inferred rotational velocity distribution with 25% of the sample having $0leq$ve$leq$100,kms and the high velocity component having ve$sim 250$,kms. There is no evidence from the spatial and radial velocity distributions of the two components that they represent either field and cluster populations or different episodes of star formation. Be-type stars have also been identified. Conclusions: The bi-modal rotational velocity distribution in our sample resembles that found for late-B and early-A type stars. While magnetic braking appears to be a possible mechanism for producing the low-velocity component, we can not rule out alternative explanations.
Rotation is thought to be a major factor in the evolution of massive stars, especially at low metallicity, with consequences for their chemical yields, ionizing flux and final fate. Determining the natal rotation-rate distribution of stars is of high
priority given its importance as a constraint on theories of massive star formation and as input for models of stellar populations in the local Universe and at high redshift. Recently, it has become clear that the majority of massive stars interact with a binary companion before they die. We investigate how this affects the distribution of rotation rates. For this purpose, we simulate a massive binary-star population typical for our Galaxy assuming continuous star formation. We find that, because of binary interaction, 20^+5_-10% of all massive main-sequence stars have projected rotational velocities in excess of 200km/s. We evaluate the effect of uncertain input distributions and physical processes and conclude that the main uncertainties are the mass transfer efficiency and the possible effect of magnetic braking, especially if magnetic fields are generated or amplified during mass accretion and stellar mergers. The fraction of rapid rotators we derive is similar to that observed. If indeed mass transfer and mergers are the main cause for rapid rotation in massive stars, little room remains for rapidly rotating stars that are born single. This implies that spin down during star formation is even more efficient than previously thought. In addition, this raises questions about the interpretation of the surface abundances of rapidly rotating stars as evidence for rotational mixing. Furthermore, our results allow for the possibility that all early-type Be stars result from binary interactions and suggest that evidence for rotation in explosions, such as long gamma-ray bursts, points to a binary origin.
We highlight the role of the light elements (Li, Be, B) in the evolution of massive single and binary stars, which is largely restricted to a diagnostic value, and foremost so for the element boron. However, we show that the boron surface abundance i
n massive early type stars contains key information about their foregoing evolution which is not obtainable otherwise. In particular, it allows to constrain internal mixing processes and potential previous mass transfer event for binary stars (even if the companion has disappeared). It may also help solving the mystery of the slowly rotating nitrogen-rich massive main sequence stars.
We study the convection zones in the outer envelope of hot massive stars which are caused by opacity peaks associated with iron and helium ionization. We determine the occurrence and properties of these convection zones as function of the stellar par
ameters. We then confront our results with observations of OB stars. A stellar evolution code is used to compute a grid of massive star models at different metallicities. In these models, the mixing length theory is used to characterize the envelope convection zones. We find the iron convection zone (FeCZ) to be more prominent for lower surface gravity, higher luminosity and higher initial metallicity. It is absent for luminosities below about $10^{3.2}Lsun$, $10^{3.9}Lsun$, and $10^{4.2}$Lsun$ for the Galaxy, LMC and SMC, respectively. We map the strength of the FeCZ on the Hertzsprung-Russell diagram for three metallicities, and compare this with the occurrence of observational phenomena in O stars: microturbulence, non-radial pulsations, wind clumping, and line profile variability. The confirmation of all three trends for the FeCZ as function of stellar parameters by empirical microturbulent velocities argues for a physical connection between sub-photospheric convective motions and small scale stochastic velocities in the photosphere of O- and B-type stars. We further suggest that clumping in the inner parts of the winds of OB stars could be caused by the same mechanism, and that magnetic fields produced in the FeCZ could appear at the surface of OB stars as diagnosed by discrete absorption components in ultraviolet absorption lines.
We have previously analysed the spectra of 135 early B-type stars in the LMC and found several groups of stars that have chemical compositions that conflict with the theory of rotational mixing. Here we extend this study to Galactic and SMC metallici
ties with the analysis of ~50 Galactic and ~100 SMC early B-type stars with rotational velocities up to ~300km/s. The surface nitrogen abundances are utilised as a probe of the mixing process. In the SMC, we find a population of slowly rotating nitrogen-rich stars amongst the early B type core-hydrogen burning stars, similar to the LMC. In the Galactic sample we find no significant enrichment amongst the core hydrogen-burning stars, which appears to be in contrast with the expectation from both rotating single-star and close binary evolution models. However, only a small number of the rapidly rotating stars have evolved enough to produce a significant nitrogen enrichment, and these may be analogous to the non-enriched rapid rotators previously found in the LMC sample. Finally, in each metallicity regime, a population of highly enriched supergiants is observed, which cannot be the immediate descendants of core-hydrogen burning stars. Their abundances are, however, compatible with them having gone through a previous red supergiant phase. Together, these observations paint a complex picture of the nitrogen enrichment in massive main sequence and supergiant stellar atmospheres, where age and binarity cause crucial effects. Whether rotational mixing is required to understand our results remains an open question at this time, but could be answered by identifying the true binary fraction in those groups of stars that do not agree with single-star evolutionary models (abridged).
One of the main uncertainties in evolutionary calculations of massive stars is the efficiency of internal mixing. It changes the chemical profile inside the star and can therefore affect the structure and further evolution. We demonstrate that ecli
psing binaries, in which the tides synchronize the rotation period of the stars and the orbital period, constitute a potentially strong test for the efficiency of rotational mixing. We present detailed stellar evolutionary models of massive binaries assuming the composition of the Small Magellanic Cloud. In these models we find enhancements in the surface nitrogen abundance of up to 0.6 dex.
