In discussing open question in the field of massive stars, I consider their evolution from birth to death. After touching upon massive star formation, which may be bi-modal and not lead to a zero-age main sequence at the highest masses, I consider the consequences of massive stars being close to their Eddington limit. Then, when discussing the effects of a binary companion, I highlight the importance of massive Algols and contact binaries for understanding the consequences of mass transfer, and the role of binaries in forming Wolf-Rayet stars. Finally, a discussion on pair instability supernovae and of superluminous supernovae is provided.

An overview is provided of the scientific goals of the Magellanic Cloud component of the STScI Directors Discretionary UV initiative ULLYSES, together with the complementary spectroscopic survey XShootU (VLT/Xshooter) and other ancillary datasets. Together, ULLYSES and XShootU permit the first comprehensive, homogeneous study of wind densities and velocities in metal-poor massive stars, plus UV/optical spectroscopic libraries for population synthesis models and a large number of interstellar sight-lines towards the Magellanic Clouds.

Rest-frame far-ultraviolet spectra are fundamental to our understanding of star-forming galaxies, providing a unique window on massive star populations (MSs), chemical evolution, feedback processes, and reionization. JWST is ushering in a new era, pushing the FUV frontier beyond z=10. The success of such endeavors hinges on a comprehensive understanding of the MSs and gas conditions that power the observed spectra. The COS Legacy Archive Spectroscopic SurveY (CLASSY) is a powerful and promising solution providing high-quality, high-resolution FUV spectra of 45 nearby star-forming galaxies. The spectra contain a suite of features that simultaneously characterize the MSs that populate metal-poor galaxies, physical properties of large-scale outflows, and chemical abundance patterns. The CLASSY sample is consistent with the z 0 mass-metallicity relationship and spans 1.5 dex in metallicity. These unique properties make CLASSY the benchmark training set for studies of MSs in star-forming galaxies both across cosmic time and connecting metal-poor to metal-rich populations.

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.

We analyzed archival HST and IUE ultraviolet spectra of 29 nearby star-forming galaxies. The range of aperture sizes permits studies of the galaxy properties over pc to kpc scales. We measured line strengths and spectral energy distributions over the 1200 – 1300 Åwavelength range and established trends with galaxy properties. Updated oxygen abundances were measured from ancillary optical data. Star-formation rates and internal dust attenuations were derived from comparison with synthesis models. The interstellar absorption lines are heavily saturated, yet scale with oxygen abundance. We interpret this as due to macroscopic velocities arising in a turbulent ISM and large-scale outflows. The stellar-wind lines also scale with oxygen abundance. As these lines are shaped by mass loss, which is driven by the Fe abundance, we can study the α-element/Fe ratio in these galaxies.

Local Group (LG) very metal-poor massive stars are the best proxy for the First Stars of the Universe, fundamental to model the early evolution of the first galaxies, and key to unravelling new evolutionary pathways restricted to low metallicities, such as chemically homogeneous evolution. Yet, due to the great leap in distance required to reach metal-poor galaxies of the LG and vicinity, no comprehensive spectroscopic studies have been carried out at sub-$$$$ metallicities so far.

We focus on the massive star population of the 1/10Ȯ galaxy Sextans A. After five observing campaigns at the 10.4-m Gran Telescopio Canarias (GTC), we have assembled a spectroscopic catalogue of more than 150 OB stars. This catalogue will be fundamental to test stellar evolution at very low metallicity, to detect the first binary systems at 1/10Ȯ and to unveil the most recent star formation history of this galaxy.

Super star cluster (SSC) A1 in starburst galaxy NGC 3125 has the strongest broad He II λ1640 emission line ever observed in the nearby Universe and constitutes an important template for interpreting observations of galaxies that are located out to a redshift of z∼3. We use observations of SSC A1 obtained with the Cosmic Origins Spectrograph (COS) on board of the Hubble Space Telescope (HST) in order to check if there is a contribution of nebular emission to the He II line. In addition, we compare the COS G130M + G160M observations of A1 (1150 – 1750∘A) to the latest Charlot & Bruzual population synthesis models, which account for Very Massive Stars (VMS) of up to 300 Mȯ. A model with Z = 0.008 and age = 2.4 Myr provides a very reasonable fit to the C III λ1175, N V λ1240, C IV λ1550, He II λ1640, and N IV λ1718 stellar-wind features, although the O V λ1371 line is not well reproduced. Overall, our results show the great improvement of stellar evolution and population synthesis models over the past decade, and in particular, the improved formulation of stellar mass loss rates.

We probe how common extremely rapid rotation is among massive stars in the early universe by measuring the OBe star fraction in nearby metal-poor dwarf galaxies. We apply a new method that uses broad-band photometry to measure the galaxy-wide OBe star fractions in the Magellanic Clouds and three more distant, more metal-poor dwarf galaxies. We find OBe star fractions of ∼20% in the Large Magallanic Cloud (0.5Zȯ), and ∼30% in the Small Magellanic Cloud (0.2Zȯ) as well as in the so-far unexplored metallicity range 0.1 Z/Zȯ < 0.2 occupied by the other three dwarf galaxies. Our results imply that extremely rapid rotation is common among massive stars in metal-poor environments such as the early universe.

The formation of multiples has seen some significant progress over the past years mainly due to the advent and the expansion of high-angular resolution facilities. Star-forming regions are the laboratories where massive stars can be caught right after their formation phase. Still, the observational constraints and the properties of young multiple systems are poorly documented. These proceedings contain recent results about the multiplicity properties of six young O-type stars in the M17 star-forming region, observed by the means of near-IR interferometric observations, which have provided insight into the origin of massive close binaries in a cluster environment.

At the time of this meeting, the latest Gaia data release is EDR3, published on 3 December 2020, but the next one, DR3, will appear soon, on 13 June 2022. This contribution describes, on the one hand, Gaia EDR3 results on massive stars and young stellar clusters, placing special emphasis on how a correct treatment of the astrometric and photometric calibration yields results that are simultaneously precise and accurate. On the other hand, it gives a brief description of the exciting results we can expect from Gaia DR3.

We use the RIOTS4 sample of SMC field OB stars to determine the origin of massive runaways in this low-metallicity galaxy using Gaia proper motions, together with stellar masses obtained from RIOTS4 data. These data allow us to estimate the relative contributions of stars accelerated by the dynamical ejection vs binary supernova mechanisms, since dynamical ejection favors faster, more massive runaways, while SN ejection favors the opposite trend. In addition, we use the frequencies of classical OBe stars, high-mass X-ray binaries, and non-compact binaries to discriminate between the mechanisms. Our results show that the dynamical mechanism dominates by a factor of 2 – 3. This also implies a significant contribution from two-step acceleration that occurs when dynamically ejected binaries are followed by SN kicks. We update our published quantitative results from Gaia DR2 proper motions with new data from DR3.

Nebular Heii emission implies the presence of energetic photons (E≽54 eV). Despite the great deal of effort dedicated to understanding Heii ionization, its origin has remained mysterious, particularly in metal-deficient star-forming galaxies. Unfolding Heii-emitting, metal-poor starbursts at z∼0 can yield insight into the powerful ionization processes occurring in the primordial universe. Here we present a study on the origin of the extended nebular Heii emission in SBS 0335-052E, one of the most metal-poor (Z ∼ 3% Z⊙ Heii-emitter starbursts known locally. Based on optical VLT/MUSE spectroscopic and Chandra X-ray observations, and current stellar models we found that the Heii-ionization budget of SBS 0335-052E can only be produced by peculiar, nearly metal-free ionizing stars (called here “PopIII-like” stars) with a top-heavy initial mass function. This result is in line with recent simulations for PopIII star formation down to z=0.

The Gaia-ESO Survey (GES) is a large public spectroscopic survey that has collected spectra of about 100,000 stars. The survey provides not only the reduced spectra, but also the radial velocities, stellar parameters and surface abundances resulting from the analysis of the spectra. We present the work of the groups that analysed the spectra of the hottest stars in that Survey. The large temperature range that is covered (Teff = 7,000 to 50,000 K) requires the use of different analysis codes by the different groups. Eight groups each analysed part of the data, with significant overlap that allowed cross-checks. In total 17,693 spectra of 6,462 stars were analysed, most of them in 37 open star clusters. The homogenisation of all this information led to stellar parameters for 5,584 stars. Abundances for at least one of the elements He, C, N, O, Ne, Mg, Al, Si and Sc were determined for 292 stars. The GES hot star data, as well as the Survey data in general, will be of considerable use in future studies of stellar evolution and open clusters.

NGC3603 is one of the youngest massive clusters of the Milky Way which uniquely enables studying the interplay between massive star feedback and the surrounding interstellar medium. Yet, a deep infrared (IR) view of the cluster is missing. We present guaranteed time observations of NGC3603 consisting of near infrared spectroscopy taken with VLT-KMOS. This data set will provide a first, rather complete IR census.

Empirical constraints are master keys for testing theoretical evolutionary model predictions. In massive stars, the region in the Hertzsprung-Russell diagram (HRD) in which the Blue Supergiants (BSGs) are located sets several important constraints to the models, and in particular to the theoretical end of the main sequence (MS). So far, we are missing from a full quantitative spectroscopic analysis (QSA) of a sample of BSGs large enough to be statistically significant and without observational biases. We present results from a QSA of a sample of ∼700 Galactic BSGs for which we have high-resolution multi-epoch optical spectra.

The Alicante Survey of MAssive Stars in Hii Regions (A-SMASHeR) is aimed at finding the ratio of massive stars that are born in isolation. We present LIRIS/WHT images and EMIR/GTC spectra of the massive stellar content in A-SMASHeR regions. Our preliminary analysis yields ∼20% of regions hosting relatively (or truly) isolated massive stars.

The Intensity Interferometry technique consists of measuring the spatial coherence (visibility) of an object via its intensity fluctuations over a sufficient range of telescope separations (baselines). This allows us to study the size, shape and morphology of stars with an unprecedented resolution. Cherenkov telescopes have a set of characteristics that coincidentally allow for Intensity Interferometry observations: very large reflective surfaces, sensitivity to individual photons, temporal resolution of nanoseconds and the fact that they come in groups of several telescopes. In the recent years, the MAGIC Collaboration has developed a deadtime-free Intensity Interferometry setup for its two 17 m diameter Cherenkov telescopes that includes a 4-channel GPU-based real-time correlator, 410–430 nm filters and new ways of splitting its primary mirrors into submirrors using Active Mirror Control (AMC). With this setup, MAGIC can operate as a long-baseline optical interferometer in the baseline range 40–90 m, which translates into angular resolutions of 0.5-1 mas. Additionally, thanks to its AMC, it can simultaneously measure the zero-baseline correlation or, by splitting into submirrors, access shorter baselines under 17 m in multiple u-v plane orientations. The best candidates to observe with this technique are relatively small and bright stars, in other words, massive stars (O, B and A types). We will present the science cases that are currently being proposed for this setup, as well as the prospects for the future of the system and technique, like the possibility of large-scale implementation with CTA.

We use the MIDE3700 code to find effective temperatures (Teff) and surface gravities (log g) via the Barbier Chalonge Divan (BCD) method for 222 B-type stars in the ESO archive in preparation for their inclusion as an extension to the X-shooter Spectral Library (XSL). We find agreement of Δ log Teff ∼0.1σ and Δ log g ∼0.25σ of our results with a sample of literature stars. We populate a previously bare region of the XSL Kiel diagram in the ranges 9000 ≤ Teff ≤ 23000 K and 2.8 ≤ log g ≤ 4.0 dex, and thereby extend the lower age limit of XSL stellar population models by up to a factor ∼10 at [Fe/H] = −1.2 dex, and by a factor ∼2 at Solar metallicity.

The aim of this survey is the homogeneous characterization of a large sample of H ii regions with active star formation in order to detect observational trends supporting the two main models of massive star formation.

We present results from our ongoing infrared spectroscopic studies of the massive stellar content at the Center of the Milky Way (GC) and across the obscured Galactic disk. Together with the full characterization of these clusters, we seek to obtain a present day metallicity 2-D map of the inner Galaxy and characterize the influence on the bar in the chemical evolution. We will also constrain the clusters IMFs, infer the presence of possible top-heavy recent star formation histories and test massive star formation channels: clusters vs isolation.