We present the results of a magnitude-limited spectroscopic survey of Galactic Wolf-Rayet stars with the HERMES spectrograph mounted on the Mercator telescope. Using cross-correlation to measure radial velocities, we measured the observed binary fractions of the Galactic carbon- (WC) and nitrogen-rich (WN) Wolf-Rayet stars to be and . We used Monte-Carlo simulations with a Bayesian framework to derive the intrinsic multiplicity properties and found and . We find that the majority of WN binaries reside in short-period systems, similar to O stars. However, the orbital period distribution of the Galactic WC population peaks at 5000 d, a discrepancy that challenges our current understanding of binary evolution in Wolf-Rayet stars.

B supergiants (BSGs) lie on the cool end of line-driven wind regime, such that the study of their atmospheres can help us to understand the physics of line-driven winds. So far key features of their spectra, especially in the UV region, could not be reproduced consistently with atmosphere models. This represents a significant gap in our knowledge of their physical properties and behavior, which is particularly striking for BSGs on the cool side of the Bi-Stability Jump (cooler than B1). To address this problem, we analysed a sample of Galactic cool BSGs, with sufficient UV and optical coverage. None of our targets are detected in X-rays with only upper limits existing for some of them.

. We present UVIT/Astrosat UV photometry of the RSG population of the Small Cloud galaxy (SMC). As RSGs are extremely faint in the far-UV, these observations directly probe potential companion stars. From a sample of 861 SMC RSGs, we find 88 have detections at far-UV wavelengths: a clear signature of binarity. Stellar parameters are determined for both components, which allows us to study - for the first time - the mass-ratio (q) distribution of RSG binary systems. We find a flat mass-ratio distribution best describes the observations up to MRSG ∼15M⊙. We account for our main observing bias (i.e. the limiting magnitude of the UVIT survey) to determine the intrinsic RSG binary fraction of 18.8 ± 1.5 %, for mass-ratios in the range 0.3.<q<1.0 and orbital periods approximately in the range 3<log P[days]<8.

The evolutionary link between Red Supergiants and Luminous Blue variables is interesting, but still poorly understood. We present the results of a study of the Galactic candidate luminous blue variable Wray 15-906, revealed via the detection of its infrared circumstellar shell (of ≍2 pc in diameter) with the Wide-field Infrared Survey Explorer (WISE) and the Herschel Space Observatory. Using the stellar atmosphere code CMFGEN and the Gaia parallax, we found that Wray 15-906 is a relatively low-luminosity, log(L/Lȯ) ≍ 5.4, star with a temperature of 25±2 kK. In the framework of single star evolution, the obtained results suggest that Wray 15-906 is a post-red supergiant star with an initial mass of ≍ 25Mȯ and that before exploding as a supernova it could transform for a short time into a WN11h star. The presence of a shell with a mass 2.9±0.5Me indicates that Wray 15-906 has suffered substantial mass loss in the recent past.

The young open cluster NGC 6231 hosts a rich population of O-type binary stars. We study several of these eccentric short-period massive eclipsing binaries and assess their fundamental parameters. The properties of these systems make them interesting targets to study tidally induced apsidal motion. The analysis of apsidal motion offers a powerful means to obtain information about the internal structure of the stars. Indeed, since the rate of apsidal motion in a binary system is proportional to the internal structure constants of the stars composing it, its value gives direct insight into the internal structure and evolutionary state of these stars. Stellar evolution models are constructed based on the observationally-determined fundamental parameters and a theoretical rate of apsidal motion is inferred. The results are striking: Adopting standard stellar evolution models yields a theoretical rate of apsidal motion much larger than the observational value. This discrepancy results from the standard models predicting too low an efficiency of internal mixing and thus too homogenous stars in terms of density. By enforcing the theoretical rates of apsidal motion to match the observational values, enhanced mixing is required, through a large overshooting parameter and/or additional turbulent/rotational mixing. Our analysis leads to the conclusion that the chemically mixed cores in those massive stars must be more extended than anticipated from standard models.

We present our measurements of the amplitude of photometric and spectroscopic variability due to clumping in the wind of Wolf-Rayet (WR) stars. Photometric variability was assessed using TESS light-curves, while spectroscopic variations were obtained from almost 20 years of monitoring of nearly 100 classical (presumably single) stars. Our results show an apparent dependence of the variability amplitude with the stars’ surface temperature and/or terminal velocity. Our interpretation is that it supports the idea that the dominating driver of the clumps in WR winds is a sub-surface convection region.

Mass loss is a key property to understand stellar evolution and in particular for low-metallicity environments. Our knowledge has improved dramatically over the last decades both for single and binary evolutionary models. However, episodic mass loss although definitely present observationally, is not included in the models, while its role is currently undetermined. A major hindrance is the lack of large enough samples of classified stars. We attempted to address this by applying an ensemble machine-learning approach using color indices (from IR/Spitzer and optical/Pan-STARRS photometry) as features and combining the probabilities from three different algorithms. We trained on M31 and M33 sources with known spectral classification, which we grouped into Blue/Yellow/Red/B[e] Supergiants, Luminous Blue Variables, classical Wolf-Rayet and background galaxies/AGNs. We then applied the classifier to about one million Spitzer point sources from 25 nearby galaxies, spanning a range of metallicites (). Equipped with spectral classifications we investigated the occurrence of these populations with metallicity.

We present low-frequency (0.40 – 1.25 GHz) radio observations and modeling of a Fast Blue Optical Transient (FBOT), AT2018cow [Nayana & Chandra(2021)]. Our data are best modeled as an inhomogeneous synchrotron emitting region expanding into an ionized circumstellar medium. We estimate the mass-loss rate of the progenitor star and shock parameters at multiple epochs post-explosion and find that the progenitor has gone through an enhanced phase of mass-loss close to its end-of-life.

We present results from a recent study of the spin rate properties of a sample of more than 400 Galactic O-type stars surveyed by the IACOB and OWN projects. We combine vsini, Teff, and logg estimates with information about the spectroscopic binarity status for 285 of the stars in the sample, and provide a renewed overview about how the empirical distribution of projected rotational velocities in the O-star domain depends on mass, evolutionary and binary status. The obtained distributions are then compared with predictions of state-of-the-art population synthesis simulations including binary interaction, and used to provide hints about the initial velocity distribution of stars with masses in the range 15-80 M⊙.

We compare pre-supernova observations with synthetic photometry from stellar evolution models to infer the progenitor properties of the seven known progenitors of Type Ib and IIb supernovae. Our results are roughly consistent with a hydrogen mass threshold of for a Type II appearance.

Recent observations of Type II supernovae have revealed that their red-supergiant progenitors lose a significant amount of mass during the last years of their evolution. However, because it is difficult to discover supernovae within days of explosion, the diversity of mass loss in red supergiants has not yet been fully mapped. This talk presented the case of SN 2021yja, which was serendipitously imaged within hours of explosion and observed with a sub-day cadence during its rise to peak. From the exceptionally long plateau period and the high nickel mass, we infer a relatively massive red-supergiant progenitor star. However, archival imaging from the Hubble Space Telescope places a stringent upper limit of on its progenitor mass. We discuss these conflicting constraints in the context of the larger sample of exploding red supergiants. Our analysis helps illuminate the poorly understood mechanism(s) behind red-supergiant mass loss.

The recent generation of dedicated wide-field, high-cadence sky-surveys have overwhelmed discovery statistics for all manner of extra-galactic transients, and uncovered new phenomena seemingly linked to the demise of massive stars. For the more established classes of transients, such as core-collapse supernovae, surges in discoveries are allowing true population studies to provide quantitative constraints not only on the explosion properties, but also on the progenitor populations. Crucially, such population insights are benefiting from creation of samples of transients constructed with largely unbiased methods for discovery and characterisation. Surrounding these discoveries are increasing samples of extreme transients that do not fit the standard core-collapse paradigm - requiring the invocation of exotic progenitor stars and placing demands on the stellar evolution of such systems. Here I will provide a high-level observationally-driven overview of recent results related to massive stellar transients.

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 present in-progress resolution test and parameter space studies for very massive stars using MESA, showcasing current MESA version convergence studies.

We present the analysis of the dust properties of the Wolf-Rayet nebulae M 1-67 and RCW 58 around the WN8h stars WR 124 and WR 40, respectively. Modeling with the photoionization code Cloudy shows that in both nebulae the IR spectral energy distributions and ionized gas properties can be reproduced by a dust shell consisting of two populations of dust grains. Furthermore, taking into account the initial mass, the morphology and the kinematics of the nebulae we propose M 1-67 and RCW 58, together with their progenitor stars, as the first observational evidences of post-common envelope evolution in nebulae around massive stars.

We have collected a database of more than 43,000 spectra of atmospheres of massive stars. These spectra have been generated with the CMFGEN code of Hillier & Miller (1998) by systematically varying stellar parameters: effective temperature, luminosity, metallicity and mass loss rate for stars from 9 to 120 solar masses (Zsargó et al. 2020) In this work we present a web-based platform for accessing the database. The platform allows an online comparison between an observed and a synthetic spectrum to quickly assess the stellar and wind parameters. The platform will be available without cost to the astronomical community and will be hosted on servers shared between Mexican Universities.

A century of study has characterized Plaskett’s Star (HD 47129) as an evolved, massive, short-period, equal mass O+O binary system. The discovery of a magnetic field in the broad-line component by Grunhut et al. (2013) renewed interest in the study of this system and led to its establishment as the most rapidly rotating magnetic O-type star. Grunhut et al. (2021) observed the circular polarization signatures of the magnetic star to exhibit no radial velocity variations while the narrow-line star demonstrates radial velocity variations consistent with the established orbital period. This has raised fundamental questions about the architecture of this system and the nature of the magnetic star which have led to a major shift in our understanding of HD 47129.

Whether it be due to rapid rotation or binary interactions, deviations from spherical symmetry are common in massive stars. These deviations from spherical symmetry are known to cause non-uniform distributions of various parameters across the surface including temperature, which can drive internal mixing processes within the envelopes of these massive stars. Despite how common these 3D distortions are, they are often neglected in spectroscopic analyses. We present a new spectral analysis code called spamms (Spectroscopic PAtch Model for Massive Stars) specifically designed to analyze non-spherical systems. We discuss how the code works and discuss its assumptions. Furthermore, we demonstrate how spamms can be applied to a variety of different types of systems and we show how it can model 3D effects in a way that current analysis techniques are not able to.

Our knowledge of massive star evolution is limited by uncertainties linked with multi-dimensional processes taking place in stellar interiors. Important examples are convective boundary mixing (CBM) and entrainment, which are implemented in 1D stellar evolution models assuming simplified prescriptions. 3D hydrodynamics models can improve these prescriptions by studying realistic multi-D processes for a short timerange (minutes or hours). In these proceedings, we present results coming from a new set of high-resolution hydrodynamics simulations of a neon-burning shell in a massive star, and discuss how the entrainment law can be calibrated from 3D models and then used to improve 1D stellar evolution prescriptions.

Especially in the upper Hertzsprung-Russell diagram, where stellar physics is least understood, obtaining model independent masses is of great value. Spectroscopic binaries that are also resolved astrometrically are an excellent alternative to eclipsing double-lined spectroscopic binaries where dynamical masses can be measured. 9 Sgr is such a massive binary. However, its characterization is troubled by conflicting conclusions from the spectroscopic analysis on the one hand and the interferometric one on the other hand. In this work, we attempt to resolve this tension by applying a novel approach to spectral disentangling of the spectroscopic data to constrain better the mass of 9 Sgr.