The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

We report the detection of the far-infrared (FIR) fine-structure line of singly ionised nitrogen, \Nplusa, within the peak epoch of galaxy assembly, from a strongly lensed galaxy, hereafter ``The Red Radio Ring''; the RRR, at z = 2.55. We combine new observations of the ground-state and mid-J transitions of CO (J$_{\rm up} =$ 1,5,8), and the FIR spectral energy distribution (SED), to explore the multi-phase interstellar medium (ISM) properties of the RRR. All line profiles suggest that the HII regions, traced by \Nplusa, and the (diffuse and dense) molecular gas, traced by the CO, are co-spatial when averaged over kpc-sized regions. Using its mid-IR-to-millimetre (mm) SED, we derive a non-negligible dust attenuation of the \Nplusa line emission. Assuming a uniform dust screen approximation results a mean molecular gas column density $> 10^{24}$\, cm$^{-2}$, with a molecular gas-to-dust mass ratio of 100. It is clear that dust attenuation corrections should be accounted for when studying FIR fine-structure lines in such systems. The attenuation corrected ratio of $L_{\rm NII205} / L_{\rm IR(8-1000\mu m)} = 2.7 \times 10^{-4}$ is consistent with the dispersion of local and $z >$ 4 SFGs. We find that the lower-limit, \Nplusa -based star-formation rate (SFR) is less than the IR-derived SFR by a factor of four. Finally, the dust SED, CO line SED and $L_{\rm NII205}$ line-to-IR luminosity ratio of the RRR is consistent with a starburst-powered ISM.

Wolf-Rayet star's winds can be so dense and so optically thick that the photosphere appears in the highly supersonic part of the outflow, veiling the underlying subsonic part of the star, and leaving the initial acceleration of the wind inaccessible to observations. We investigate the conditions and the structure of the subsonic part of the outflow of Galactic WR stars, in particular of the WNE subclass; our focus is on the conditions at the sonic point. We compute 1D hydrodynamic stellar structure models for massive helium stars adopting outer boundaries at the sonic point. We find that the outflows of our models are accelerated to supersonic velocities by the radiative force from opacity bumps either at temperatures of the order of 200kK by the Fe opacity bump or of the order of 50kK by the HeII opacity bump. For a given mass-loss rate, the conditions in the subsonic part of the outflow are independent from the detailed physical conditions in the supersonic part. The close proximity to the Eddington limit at the sonic point allows us to construct a Sonic HR diagram, relating the sonic point temperature to the L/M ratio and the stellar mass-loss rate, thereby constraining the sonic point conditions, the subsonic structure, and the stellar wind mass-loss rates from observations. The minimum mass-loss rate necessary to have the flow accelerated to supersonic velocities by the Fe opacity bump is derived. A comparison of the observed parameters of Galactic WNE stars to this minimum mass-loss rate indicates that their winds are launched to supersonic velocities by the radiation pressure arising from the Fe-bump. Conversely, models which do not show transonic flows from the Fe opacity bump form inflated envelopes. We derive an analytic criterion for the appearance of envelope inflation in the subphotospheric layers.

(abridged) The strong winds of Wolf-Rayet (WR) stars are important for the mechanical and chemical feedback of the most massive stars and determine whether they end their lives as neutron stars or black holes. In this work we investigate theoretically the mass-loss properties of H-free WR stars of the nitrogen sequence (WN stars). We connect stellar structure models for He stars with wind models for optically-thick winds and assess how both types of models can simultaneously fulfill their respective sonic-point conditions. Fixing the outer wind law and terminal wind velocity, we obtain unique solutions for the mass-loss rates of optically-thick, radiatively-driven winds of WR stars in the phase of core He-burning. The resulting mass-loss relations as a function of stellar parameters, agree well with previous empirical relations. Furthermore, we encounter stellar mass limits below which no continuous solutions exist. While these mass limits agree with observations of WR stars in the Galaxy, they are in conflict with observations in the LMC. While our results confirm in particular the slope of oft-used empirical mass-loss relations, they imply that only part of the observed WN population can be understood in the framework of the standard assumptions of a smooth transonic flow and compact stellar core. This means that alternative approaches, such as a clumped and inflated wind structure, or deviations from the diffusion limit at the sonic point may have to be invoked. Qualitatively, the existence of mass limits for the formation of WR-type winds may be relevant for the non-detection of low-mass WR stars in binary systems, which are believed to be progenitors of Type Ib/c supernovae. The sonic-point conditions derived in this work may provide a possibility to include optically-thick winds in stellar evolution models in a more physically motivated form than in current models.

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.

We exploit the recent discovery of pulsations in mixed-atmosphere (He/H), extremely low-mass white dwarf precursors (ELM proto-WDs) to test the proposition that rotational mixing is a fundamental process in the formation and evolution of low-mass helium core white dwarfs. Rotational mixing has been shown to be a mechanism able to compete efficiently against gravitational settling, thus accounting naturally for the presence of He, as well as traces of metals such as Mg and Ca, typically found in the atmospheres of ELM proto-WDs. Here we investigate whether rotational mixing can maintain a sufficient amount of He in the deeper driving region of the star, such that it can fuel, through HeII-HeIII ionization, the observed pulsations in this type of stars. Using state-of-the-art evolutionary models computed with MESA, we show that rotational mixing can indeed explain qualitatively the very existence and general properties of the known pulsating, mixed-atmosphere ELM proto-WDs. Moreover, such objects are very likely to pulsate again during their final WD cooling phase.

We aim to provide new empirical clues about macroturbulent spectral line broadening in O- and B-type stars to evaluate its physical origin. We use high-resolution spectra of ~430 stars with spectral types in the range O4-B9 (all luminosity classes). We characterize the line-broadening of adequate diagnostic metal lines using a combined FT and GOF technique. We perform a quantitative spectroscopic analysis of the whole sample using automatic tools coupled with a huge grid of FASTWIND models. We also incorporate quantitative information about line asymmetries to our observational description of the characteristics of the line-profiles, and present a comparison of the shape and type of line-profile variability found in a small sample of O stars and B supergiants with still undefined pulsational properties and B main sequence stars with variable line-profiles. We present a homogeneous and statistically significant overview of the (single snapshot) line-broadening properties of stars in the whole O and B star domain. We find empirical evidence of the existence of various types of non-rotational broadening agents acting in the realm of massive stars. Even though all of them could be quoted and quantified as a macroturbulent broadening from a practical point of view, their physical origin can be different. Contrarily to the early- to late-B dwarfs/giants, which present a mixture of cases in terms of line-profile shape and variability, the whole O-type and B supergiant domain (or, roughly speaking, stars with M_ZAMS > 15 M_sol) is fully dominated by stars with a remarkable non-rotational broadening component and very similar profiles (including type of variability). We provide some examples illustrating how this observational dataset can be used to evaluate scenarios aimed at explaining the existence of sources of non-rotational broadening in massive stars.

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.

A large number of extremely low-mass helium white dwarfs (ELM WDs) have been discovered in recent years. The majority of them are found in close binary systems suggesting they are formed either through a common-envelope phase or via stable mass transfer in a low-mass X-ray binary (LMXB) or a cataclysmic variable (CV) system. Here, we investigate the formation of these objects through the LMXB channel with emphasis on the proto-WD evolution in environments with different metallicities. We study, for the first time, the combined effects of rotational mixing and element diffusion (e.g. gravitational settling, thermal and chemical diffusion) on the evolution of proto-WDs and on the cooling properties of the resulting WDs. We present state-of-the-art binary stellar evolution models computed with MESA for metallicities between Z=0.0002 and Z=0.02, producing WDs with masses between 0.16-0.45 M$_{\odot}$. Our results confirm that element diffusion plays a significant role in the evolution of proto-WDs that experience hydrogen shell flashes. The occurrence of these flashes produces a clear dichotomy in the cooling timescales of ELM WDs, which has important consequences e.g. for the age determination of binary millisecond pulsars. Rotational mixing is found to counteract the effect of gravitational settling in the surface layers of young, bloated ELM proto-WDs and therefore plays a key role in determining their surface chemical abundances. We predict that these proto-WDs have helium-rich envelopes through a significant part of their lifetime, a crucial ingredient for understanding the newly observed ELM proto-WD pulsators. The hydrogen envelope at detachment, although small compared to the total mass of the WD, contains enough angular momentum such that the spin frequency of the resulting WD on the cooling track is well above the orbital frequency.

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.

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.

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.