Recently it has been found that models of massive stars reach the Eddington limit in their interior, which leads to dilute extended envelopes. We perform a comparative study of the envelope properties of massive stars at different metallicities, with the aim to establish the impact of the stellar metallicity on the effect of envelope inflation. We analyse published grids of core-hydrogen burning massive star models computed with metallicities appropriate for massive stars in the Milky Way, the LMC and the SMC, the very metal poor dwarf galaxy I Zwicky 18, and for metal-free chemical composition. Stellar models of all the investigated metallicities reach and exceed the Eddington limit in their interior, aided by the opacity peaks of iron, helium and hydrogen, and consequently develop inflated envelopes. Envelope inflation leads to a redward bending of the zero-age main sequence and a broadening of the main sequence band in the upper part of the Hertzsprung-Russell diagram. We derive the limiting L/M-values as function of the stellar surface temperature above which inflation occurs, and find them to be larger for lower metallicity. While Galactic models show inflation above ~29 Msun, the corresponding mass limit for Population III stars is ~150 Msun. While the masses of the inflated envelopes are generally small, we find that they can reach 1-100 Msun in models with effective temperatures below ~8000 K, with higher masses reached by models of lower metallicity. Envelope inflation is expected to occur in sufficiently massive stars at all metallicities, and is expected to lead to rapidly growing pulsations, high macroturbulent velocities, and might well be related to the unexplained variability observed in Luminous Blue Variables like S Doradus and Eta Carina.

Macroturbulence, introduced as a fudge to reproduce the width and shape of stellar absorption lines, reflects gas motions in stellar atmospheres. While in cool stars, it is thought to be caused by convection zones immediately beneath the stellar surface, the origin of macroturbulence in hot stars is still under discussion. Recent works established a correlation between the turbulent-to-total pressure ratio inside the envelope of stellar models and the macroturbulent velocities observed in corresponding Galactic stars. To probe this connection further, we evaluated the turbulent pressure that arises in the envelope convective zones of stellar models in the mass range 1-125 Msun based on the mixing-length theory and computed for metallicities of the Large and Small Magellanic Cloud. We find that the turbulent pressure contributions in models with these metallicities located in the hot high-luminosity part of the Hertzsprung-Russel (HR) diagram is lower than in similar models with solar metallicity, whereas the turbulent pressure in low-metallicity models populating the cool part of the HR-diagram is not reduced. Based on our models, we find that the currently available observations of hot massive stars in the Magellanic Clouds appear to support a connection between macroturbulence and the turbulent pressure in stellar envelopes. Multidimensional simulations of sub-surface convection zones and a larger number of high-quality observations are necessary to test this idea more rigorously.

The envelopes of stars near the Eddington limit are prone to various instabilities. A high Eddington factor in connection with the Fe opacity peak leads to convective instability, and a corresponding envelope inflation may induce pulsational instability. Here, we investigate the occurrence and consequences of both instabilities in models of Wolf-Rayet stars. We determine the convective velocities in the sub-surface convective zones to estimate the amplitude of the turbulent velocity at the base of the wind that potentially leads to the formation of small-scale wind structures, as observed in several WR stars. We also investigate the effect of mass loss on the pulsations of our models. We approximated solar metallicity WR stars by models of mass-losing helium stars, and we characterized the properties of convection in the envelope adopting the standard MLT. Our results show the occurrence of sub-surface convective regions in all studied models. Small surface velocity amplitudes are predicted for models with masses below 10Msun. For models with M>10Msun, the surface velocity amplitudes are of the order of 10km/s. Moreover we find the occurrence of pulsations for stars in the mass range 9-14Msun, while mass loss appears to stabilize the more massive WR stars. We confront our results with observationally derived line variabilities of 17 WN stars. The data suggest variability to occur for stars above 10Msun, which is increasing linearly with mass above this value, in agreement with our results. We further find some of our models to be unstable to radial pulsations, and predict local magnetic fields of the order of hundreds of Gauss in WR stars more massive than 10Msun. Our study relates the surface velocity fluctuations induced by sub-surface convection to the formation of clumping in the inner part of the wind. From this mechanism, we expect a stronger variability in more massive WR stars.

A significant fraction of the envelope of low- and intermediate-mass stars is unstable to convection, leading to sub-surface turbulent motion. Here, we consider and include the effects of turbulence pressure in our stellar evolution calculations. In search of an observational signature, we compare the fractional contribution of turbulent pressure to the observed macroturbulent velocities in stars at different evolutionary stages. We find a strong correlation between the two quantities, similar to what was previously found for massive OB stars. We therefore argue that turbulent pressure fluctuations of finite amplitude may excite high-order, high-angular degree stellar oscillations, which manifest themselves at the surface an additional broadening of the spectral lines, i.e., macroturbulence, across most of the HR diagram. When considering the locations in the HR diagram where we expect high-order oscillations to be excited by stochastic turbulent pressure fluctuations, we find a close match with the observational $\gamma$ Doradus instability strip, which indeed contains high-order, non-radial pulsators. We suggest that turbulent pressure fluctuations on a percentual level may contribute to the $\gamma$ Dor phenomenon, calling for more detailed theoretical modelling in this direction.

Massive, luminous stars reaching the Eddington limit in their interiors develop very dilute, extended envelopes. This effect is called envelope inflation. If the progenitors of Type Ib/c supernovae, which are believed to be Wolf-Rayet (WR) stars, have inflated envelopes then the shock breakout signals diffuse in them and can extend their rise times significantly. We show that our inflated, hydrogen-free, WR stellar models with a radius of ~Rsun can have shock breakout signals longer than ~60 s. The puzzlingly long shock breakout signal observed in the Type Ib SN 2008D can be explained by an inflated progenitor envelope, and more such events might argue in favour of existence of inflated envelopes in general.

The major mass fraction of the envelope of hot luminous stars is radiatively stable. However, the partial ionisation of hydrogen, helium and iron gives rise to extended sub-surface convection zones in all of them. In this work, we investigate the effect of the pressure induced by the turbulent motion in these zones based on the mixing length theory, and search for observable consequences. We find that the turbulent pressure fraction can amount up to ~5% in OB supergiants, and to ~30% in cooler supergiants. The resulting structural changes are, however, not significantly affecting the evolutionary tracks compared to previous calculations. Instead, a comparison of macroturbulent velocities derived from high quality spectra of OB stars with the turbulent pressure fraction obtained in corresponding stellar models reveals a strong correlation of these two quantities. We discuss a possible physical connection, and conclude that turbulent pressure fluctuations may drive high-order oscillations, which - as conjectured earlier - manifest themselves as macroturbulence in the photospheres of hot luminous stars.

Context. Oxygen sequence Wolf-Rayet (WO) stars represent a very rare stage in the evolution of massive stars. Their spectra show strong emission lines of helium-burning products, in particular highly ionized carbon and oxygen. The properties of WO stars can be used to provide unique constraints on the (post-)helium burning evolution of massive stars, as well as their remaining lifetime and the expected properties of their supernovae. Aims. We aim to homogeneously analyse the currently known presumed-single WO stars to obtain the key stellar and outflow properties and to constrain their evolutionary state. Methods. We use the line-blanketed non-local thermal equilibrium atmosphere code cmfgen to model the X-Shooter spectra of the WO stars and deduce the atmospheric parameters. We calculate dedicated evolutionary models to determine the evolutionary state of the stars. Results. The WO stars have extremely high temperatures that range from 150 kK to 210 kK, and very low surface helium mass fractions that range from 44% down to 14%. Their properties can be reproduced by evolutionary models with helium zero-age main sequence masses of $M_{\mathrm{He, ini}} = 15-25 M_{\odot}$ that exhibit a fairly strong (on the order of a few times $10^{-5} M_{\odot} \mathrm{yr}^{-1}$), homogeneous ($f_\mathrm{c} > 0.3$) stellar wind. Conclusions. WO stars represent the final evolutionary stage of stars with estimated initial masses of $M_{\mathrm{ini}} = 40-60 M_{\odot}$. They are post core-helium burning and predicted to explode as type Ic supernovae within a few thousand years.

Massive rotating single stars with an initial metal composition appropriate for the dwarf galaxy I Zw 18 ([Fe/H]=$-$1.7) are modelled during hydrogen burning for initial masses of 9-300 M$_{\odot}$ and rotational velocities of 0-900 km s$^{-1}$. Internal mixing processes in these models were calibrated based on an observed sample of OB-type stars in the Magellanic Clouds. Even moderately fast rotators, which may be abundant at this metallicity, are found to undergo efficient mixing induced by rotation resulting in quasi chemically-homogeneous evolution. These homogeneously-evolving models reach effective temperatures of up to 90 kK during core hydrogen burning. This, together with their moderate mass-loss rates, make them Transparent Wind Ultraviolet INtense stars (TWUIN star), and their expected numbers might explain the observed HeII ionizing photon flux in I Zw 18 and other low-metallicity HeII galaxies. Our slowly rotating stars above $\sim$80 M$_{\odot}$ evolve into late B- to M-type supergiants during core hydrogen burning, with visual magnitudes up to 19$^{\mathrm{m}}$ at the distance of I Zw 18. Both types of stars, TWUIN stars and luminous late-type supergiants, are only predicted at low metallicity. Massive star evolution at low metallicity is shown to differ qualitatively from that in metal-rich environments. Our grid can be used to interpret observations of local star-forming dwarf galaxies and high-redshift galaxies, as well as the metal-poor components of our Milky Way and its globular clusters.

The evolution of massive stars even on the main sequence is not yet well understood. Due to the steep mass-luminosity relation, massive main sequence stars become very luminous. This brings their envelopes very close to the Eddington limit. We are analysing stellar evolutionary models in which the Eddington limit is reached and exceeded, and explore the rich diversity of physical phenomena which take place in their envelopes, and investigate their observational consequences. We use the grids of detailed stellar models by Brott et al. (2011) and Koehler et al. (2015), to investigate the envelope properties of core hydrogen burning massive stars. We find that at the stellar surface, the Eddington limit is almost never reached, even for stars up to 500 Msun. When an appropriate Eddington limit is defined locally in the stellar envelope, most stars more massive than 40 Msun actually exceed this limit, in particular in the partial ionization zones of iron, helium or hydrogen. While most models adjust their structure such that the local Eddington limit is exceeded at most by a few per cent, our most extreme models do so by a factor of more than seven. We find that the local violation of the Eddington limit has severe consequences for the envelope structure, as it leads to envelope inflation, convection, density inversions and possibly to pulsations. We find that all models with luminosities higher than 4*10^5 Lsun i.e. stars above 40 Msun show inflation, with a radius increase of up to a factor of about 40. We find that the hot edge of the S Dor variability region coincides with a line beyond which our models are inflated by more than a factor of two, indicating a possible connection between S Dor variability and inflation. Furthermore, our coolest models show highly inflated envelopes with masses of up to several solar masses, and appear to be candidates to produce major LBV eruptions.

Stars close to the Eddington luminosity can have large low-density inflated envelopes. We show that the rise times of shock breakout signals from supernovae can be extended significantly if supernova progenitors have an inflated stellar envelope. If the shock breakout occurs in such inflated envelopes, the shock breakout signals diffuse in them, and their rise time can be significantly extended. Then, the rise times of the shock breakout signals are dominated by the diffusion time in the inflated envelope rather than the light-crossing time of the progenitors. We show that our inflated Wolf-Rayet star models whose radii are of the order of the solar radius can have shock breakout signals which are longer than ~100 sec. The existence of inflated envelopes in Wolf-Rayet supernova progenitors may be related to the mysterious long shock breakout signal observed in Type Ib SN 2008D. Extended shock breakout signals may provide evidence for the existence of inflated stellar envelopes and can be used to constrain the physical properties of these enigmatic structures.

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 fields, 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.

Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs), formed via type Ib/c core-collapse supernovae (SNe), which have been spun up to high rotation rates via accretion from a companion star in a low-mass X-ray binary (LMXB). In an alternative formation channel, NSs are produced via the accretion-induced collapse (AIC) of a massive white dwarf (WD) in a close binary. Here we investigate binary evolution leading to AIC and examine if NSs formed in this way can subsequently be recycled to form MSPs and, if so, how they can observationally be distinguished from pulsars formed via the standard core-collapse SN channel in terms of their masses, spins, orbital periods and space velocities. Numerical calculations with a detailed stellar evolution code were used for the first time to study the combined pre- and post-AIC evolution of close binaries. We investigated the mass transfer onto a massive WD in 240 systems with three different types of non-degenerate donor stars: main-sequence stars, red giants, and helium stars. When the WD is able to accrete sufficient mass (depending on the mass-transfer rate and the duration of the accretion phase) we assumed it collapses to form a NS and we studied the dynamical effects of this implosion on the binary orbit. Subsequently, we followed the mass-transfer epoch which resumes once the donor star refills its Roche lobe and calculated the continued LMXB evolution until the end. We demonstrate that the final properties of these MSPs are, in general, remarkably similar to those of MSPs formed via the standard core-collapse SN channel. However, the resultant MSPs created via the AIC channel preferentially form in certain orbital period intervals. Finally, we discuss the link between AIC and young NSs in globular clusters. Our calculations are also applicable to progenitor binaries of SNe Ia under certain conditions. [Abridged]