Hot, massive stars are known to host unstable, radiation-driven outflowing winds, giving rise to dense clumps of material which severely affect the diagnostic techniques used to derive wind properties of massive stars. Most of the current diagnostic models account for wind inhomogeneities by assuming a one-component medium consisting of optically thin clumps, and maintaining a smooth velocity-field. However, this neglects important light-leakage effects through porous channels in-between the clumps. These light-leakage effects have recently been incorporated in the stellar atmosphere modelling code FASTWIND, and here we will present quantitative mass-loss results from a combined Ultraviolet-Optical wind analysis of O-supergiants in the Galaxy. Using a genetic-algorithm fitting-approach, we systematically investigate the impact the wind physics has on derived stellar and wind parameters, and how this depends on metallicity and spectral type. We compare our findings with earlier results (which do not take into account such light-leakage effects), to standard mass-loss rates usually included in evolution model studies of massive stars, and with theoretical predictions of clumping properties. We will also present the first systematic empirical constraints on the new wind parameters, associated with light-leakage, and compare these with theoretical predictions.

I describe (i) our recent updates on first star formation, with particular emphasis on their binaries, (ii) formation of low-metallicity stars and the transition of their initial mass functions with metal enrichment, and finally (iii) formation of supermassive stars from slightly metal-enriched gas by the newly found super-competitive accretion channel.

Despite the important role mass loss in the red supergiant phase plays in controlling stellar evolution and massive stars’ final supernova fates, a theoretical explanation of the mechanism driving this mass loss has been elusive. In this contribution we present a recent breakthrough [Kee et∼al., 2021] showing that turbulent pressure alone is sufficient to markedly extend the atmospheres of red supergiants and allow a wind to be launched. The resulting theory provides a fully analytic prescription for red supergiant mass-loss rates. Moreover, the theoretical mass-loss rates computed from observationally inferred turbulent velocities are in overall good agreement with observationally inferred red supergiant mass loss. A particularly interesting aspect of this theory is that it is not sensitive to metallicity, providing important implications for stellar evolution and the so-called “red-supergiant problem” for supernova progenitors in various environments.

Very metal-poor massive stars hold the key to interpret high-redshift star-forming galaxies and the early reionization epoch, but also contemporary events such as gravitational waves. To study these objects in resolved environments, we need to resort to dwarf irregular galaxies far from the potential wells of M31 and the Milky Way, and therefore distant. While the archives, recently boosted by the ULLYSES and XSHOOTU programs, store a healthy dataset of massive stars in the Milky Way and the Magellanic Clouds, the number of observed targets with poorer metal content than the SMC (1/5 Zȯ) is dramatically small. This paper reviews the state of observations of very metal-poor massive stars, assessing what can be realistically learned about their physics and evolution with current instrumentation, and arguing whether or not near-future facilities can remedy the gaps in the knowledge that remain.

Rotation plays an important role in the structure and evolution of massive stars. It leads to deviation from spherical symmetry for very fast rotating stars, mixing in otherwise unmixed radiative regions and generally increased mass loss. In addition, magnetic fields interact with rotation and lead to significant transport of angular momentum. In this article, we review the various rotational and magnetic instabilities present in massive stars and their implementation in one-dimensional stellar evolution codes. We then focus on their impact on the evolution of single rotating stars. Finally, we compare rotating models to observations and discuss ways to disentangle between various uncertainties.

Low metallicity (Z) massive stars are among the main feedback agents in the early Universe and in present-day blue dwarf galaxies. The nearby star-forming SMC galaxy offers conditions which resemble those at redshift z∼2 i.e. where modern galaxies formed and star formation peaked. Here we present the recent results about the nature of the eclipsing O-type binary in the SMC, AzV 476, to gain insights on the properties of massive stars and binaries at earlier cosmic epochs. We find that the primary has surprisingly low mass while being much brighter and hotter than the secondary. To place the measured stellar properties in the evolutionary context we modeled the system and confirm that AzV 476 is a post-interaction binary with the primary already being core helium (He) burning, while still having a hydrogen-rich (H-rich) envelope. These results constrain massive binary evolutionary scenarios and guide the searches of stripped stars in low-Z environments.

Luminous Blue Variable stars (LBVs) are rare and enigmatic. Often cited as evolutionary stages in the single-star evolution, the idea that binary evolution produces the LBV state was already considered, 30 years ago.

It is now commonly accepted that a significant part of massive stars are born in multiple systems. One aspect that also emerged is that massive stars have on average at least two companions, i.e. they are triples. This immediately implies that a number of LBVs should have evolved as part of multiple systems.

While some LBVs are confirmed as binaries, different methods were used to derive their multiplicity, with different results. We report on a systematic search for multiplicity using spectroscopy, interferometry in a sample of 20 LBVs. Spectroscopy provides us with a bias-corrected binary fraction of $\[62_{ - 24}^{ + 38}\]$%, and a percentage of 50–70% is found from interferometry. This has a high impact on the way that these objects might be formed.

The temperature independent part of the Humphreys-Davidson (HD) limit sets the boundary for evolutionary channels of massive stars that either end their lives as red supergiants (RSGs) or as the hotter blue supergiants (BSGs) and Wolf-Rayet stars. Recent downward revision of most luminous RSGs the Galaxy below log(L / L⊙) ≈ 5.5, more in line with the Magellanic Clouds, might hint towards a metallicity (Z)-independent HD limit. We present MESA single star models in the 15-40 M⊙ range and study the different Z-dependent processes that could potentially affect the location of the upper luminosity limit of RSGs.

Since massive stars form preferentially as members of close binary systems, we use dense grids of detailed binary evolution models to explore how binary evolution shapes the main-sequence morphology of young star clusters. We propose that binary mergers might be the origin of the blue main sequence stars in young star clusters. Our results imply that stars may either form by accretion, or through a binary merger, and that both paths lead to distinctly different spins, magnetic fields, and stellar mass distributions.

In this contribution, we explore the question on the formation of multiple massive stellar systems via disk fragmentation with the help of the highest-resolution simulations to date of a fragmenting disk in the context of massive star formation. The simulations start from a collapsing cloud of 200 solar masses, followed by the formation of an accretion disk that develops spiral arms and fragments. Due to the high resolution of our grid, we are able to self-consistently form the fragments without the need for a subgrid module such as sink particles. We track the formed fragments into the first stages of companion formation, which allows us to give an estimate of the multiplicity of the final system due to disk fragmentation. We find in total around ∼6 fragments, some at orbits of ∼ 1000 au, and some close (possibly spectroscopic) companions.

MWC 656 has been reported as classical Be star with a black hole companion. Revisited spectral variability properties render this unlikely, with a hot subdwarf more probable.

Exploring the low-mass end of the companion mass function around massive stars is of crucial importance to constrain massive star formation theories. We present a high-contrast imaging study of 20 O- and early B-type stars in the Scorpius OB1 association. From the analysis of VLT/SPHERE data, we identify a total of 789 sources. The data probe the brown dwarf regime around massive stars, resulting in the discovery of large-separation multiple systems with mass-ratios as low as 0.001 (comparable to Jupiter-Sun mass-ratio).

One significant difficulty in reliable quantification of the rates of mass-loss from hot, massive stars lies in uncertainties associated with quantifying temporal and spatial variability within stellar winds. The consequences of low-metallicity conditions for wind structure also merit continued investigation. We present initial results from ULLYSES data with the aim of identifying structure within the stellar winds of early B type supergiants with sub-solar metallicities in the Large and Small Magellanic Clouds. We demonstrate how single-epoch ULLYSES data can be used to investigate significant wind structure for these stars.

We investigate the physical properties of dust in the environment of three core-collapse supernovae (CCSNe) through mid-infrared (mid-IR) spectral energy distribution (SED) modeling (both analytical and numerical methods) and interpret our results within a Bayesian framework. We provide evidence that the observed late-time mid-IR excess of the SNe can be described by dust models. We conclude that in case of various types of SNe, numerical dust models with a shell-like geometry can be reconciled with analytical models, regarding the essential properties of dust grains.

Mass loss through stellar winds plays a dominant role in the evolution of massive stars. Very massive stars (VMSs, > 100Mȯ) display Wolf-Rayet spectral morphologies (WNh) whilst on the main-sequence. Bestenlehner (2020) extended the elegant and widely used stellar wind theory by Castor, Abbott & Klein (1975) from the optically thin (O star) to the optically thick main-sequence (WNh) wind regime. The new mass-loss description is able to explain the empirical mass-loss dependence on the Eddington parameter and is suitable for incorporation into stellar evolution models for massive and very massive stars. The prescription can be calibrated with the transition mass-loss rate defined in Vink & Gräfener (2012). Based on the stellar sample presented in Bestenlehner et al. (2014) we derive a mass-loss recipe for the Large Magellanic Cloud using the new theoretical mass-loss prescription of Bestenlehner (2020).

Adequate stellar atmosphere models are prerequisite to derive robust stellar parameters from spectroscopic analyses. I will briefly review recent results obtained with the Potdam Wolf-Rayet (PoWR) model atmosphere code, which is applicable to all types of hot stars. Using multi-wavelength observations including the UV, we analyzed large samples of massive stars at various metallicities, gaining important insights on their cosmic role and the feedback to their environment.

A recent extension of PoWR allows to compose the model atmosphere from two zones. A rapidly rotating star, e.g., might possess a cooler equatorial region with a slow wind, and two polar cones with higher photospheric temperature and fast wind. For two examples of rapidly rotating O-type stars, we demonstrate that such model can reproduce wind-line profiles which otherwise would stay inconsistent. Fast rotation, which prevails in particular at low metallicities, thus might bias empirically derived parameters, having implications for feedback as well as for angular-momentum losses of SN and GRB progenitors.

Upcoming large-scale spectroscopic surveys such as WEAVE and 4MOST will provide thousands of spectra of massive stars, which need to be analysed in an efficient and homogeneous way. Studies on massive stars are usually based on samples of a few hundred objects which pushes current spectroscopic analysis tools to their limits because visual inspection is necessary to verify the spectroscopic fit.

The novel spectroscopic analysis pipeline takes advantage of the statistics that large samples provide, and determines the model error to account for imperfections in stellar atmosphere codes due to simplified, wrong or missing physics. Considering observational plus model uncertainties improve spectroscopic fits. The pipeline utilises the entire spectrum rather than selected diagnostic lines allowing a wider range of temperature from B to early O stars to be analysed. A small fraction of stars like peculiar, contaminated or spectroscopic binaries require visual inspection, which are identified through their larger uncertainties.

B-type supergiants show enormous potential as resourceful tools to address a wide range of astrophysical questions concerning stellar atmospheres, stellar and galactic evolution and the cosmic distance scale. For the purposes of a comprehensive analysis of these objects we test a hybrid non-LTE approach – line-blanketed model atmospheres computed under the assumptions of local thermodynamic equilibrium (LTE) in combination with non-LTE line-formation calculations. An observational sample of 14 Galactic B-type supergiants with masses below about 30 Mȯ is investigated on the basis of high-resolution Echelle spectra. The results of this analysis – atmospheric and fundamental stellar parameters, the characterisation of the interstellar sightlines to the objects, as well as derived spectroscopic distances and multi-species abundances – are subjected to multiple tests of consistency.

There is evidence that some red supergiants (RSGs) experience phases of episodic mass-loss. These episodes yield more extreme mass-loss rates, further stripping the envelope of the RSG, significantly affecting the further evolution towards the final collapse of the star. Mass lost through RSG outbursts/superwinds will flow outwards and form dust further out from the surface and this dust may be detected and modelled. Here, we aim to derive the surface properties and estimate the global properties of Mid-IR bright RSGs in the Magellanic Clouds. These properties will then be compared to evolutionary predictions and used for future spectral energy distribution fitting studies to measure the mass-loss rates from present circumstellar dust.

In this paper, we present a glimpse of our observations of two Wolf-Rayet (WR) nebulae, NGC2359 and NGC6888 obtained with the SITELLE imaging Fourier transform spectrograph. The data are of unprecedented spatial coverage and cover a broad wavelength range.