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.

Direct observations of the products of binary interactions are sparse, yet they provide important insights on the outcome of the interaction and the physics at play. Young and intermediate-age star clusters are the ideal tool to search for, and characterize such interaction products and allow for a detailed comparison to theoretical predictions. We here report on integral field spectroscopy obtained with MUSE for several such clusters in the Magellanic Clouds.

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 HII environment metallicities of 65 supernovae (SNe), obtained with INT/IDS, have been determined using the N2 and O3N2 strong emission line methods. Resulting cumulative distribution functions reveal a narrower distribution for Ib SNe (standard deviation σ ∼ 0.06 dex) compared to Ic and IIP distributions (σ ∼ 0.15 dex). This narrow distribution of Ib SNe is confirmed with an extended dataset using data from Galbany et al. (2018). Statistical tests confirm a statistically significant difference between the Ib and II-P metallicity distributions with < 5% probability that they result from the same progenitors. This narrow distribution suggests a lack of Type Ib SNe in low metallicity environments and points towards single star progenitors for these Type Ib SNe, rather than binaries. It also suggests that single massive stars at low metallicity are not commonly able to produce helium-rich Type Ib supernovae.

We present the results obtained using spectroscopic data taken with the intermediate-resolution Multi Unit Spectroscopic Explorer (MUSE) of B and A-type supergiants and bright giants in the Sculptor Group galaxy NGC 300. For our analysis, a hybrid local thermodynamic equilibrium (LTE) line-blanketing+non-LTE method was used to improve the previously published results for the same data. In addition, we present some further applications of this work, which includes extending the flux-weighted gravity luminosity relationship (FGLR), a distance determination method for supergiants. This pioneering work opens up a new window to explore this relation, and also demonstrates the enormous potential of integral field spectroscopy (IFS) for extragalactic quantitative stellar studies.

The empirical upper limit to Red Supergiant (RSG) luminosity, known as the Humphreys-Davidson (HD) limit, has been commonly explained as being caused by the stripping of stellar envelopes by metallicity-dependent, line-driven winds. As such, the theoretical expectation is that the HD limit should be higher at lower metallicity, where weaker mass-loss rates mean that higher initial masses are required for an envelope to be stripped. In this work, we test this prediction by measuring the luminosity function of RSGs in M31 and comparing to those in the LMC and SMC. We find that $\[\log ({L_{{\rm{m}}ax}}/{L_ \odot }) = 5.53 \pm 0.03\]$ in M31 (Z ≳ Z⊙), consistent with the limit found for both the LMC (Z ∼ 0.5 Z⊙) and SMC (Z ∼ 0.25 Z⊙), while the RSG luminosity distributions in these 3 galaxies are consistent to within 1σ. We therefore find no evidence for a metallicity dependence on both the HD limit and the RSG luminosity function, and conclude that line-driven winds on the main sequence are not the cause of the HD limit.

Stellar mergers produce more massive, rejuvenated (strongly magnetic) stars, with potentially peculiar properties, and can be detected as luminous red novae. Using a grid of detailed 1D binary evolution models, we aim to determine which binary systems are likely to merge and at what evolutionary stage. This will tell us more about the merger products, and might help us understand some of the trends found in observed single- and multiple-star populations.

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.

Long gamma-ray bursts (LGRBs) and superluminous supernovae (SLSNe) are expected to result from massive star deaths. However to date, there has been no direct observational measurement of their cloud collapse timescales nor progenitor lifetimes to help constrain their mass. Our analyses of z 2 LGRB afterglow spectra and Hubble Space Telescope images find a higher fraction of host galaxies that are interacting, have a close companion, and/or may have experienced a recent galaxy ‘fly by’ as compared to the general z 2 galaxy population. A smaller set of z 2 SLSNe suggests a similar result. Under the hypothesis that galaxy interactions induce cloud collapse and star formation near their closest approach, we explore measurements of the host and companion galaxy velocities and separations at the time of the LRGB/SLSN event as a direct physical means to measure the timescale of cloud collapse plus progenitor star lifetime.

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.

We report recent ESPaDOnS and HARPSpol spectropolarimetric observations from our ongoing magnetic survey of the brightest twenty-five classical Cepheids. Stokes V magnetic signatures are detected in eight of fifteen targets observed to date. The Stokes V profiles show a diversity of morphologies with weak associated longitudinal field measurements of order 1 G. Many of the Stokes V profiles are difficult to interpret in the context of the normal Zeeman effect. They consist of approximately unipolar single or double lobe(s) of positive or negative circular polarization. We hypothesize that these unusual signatures are due to the Zeeman effect modified by atmospheric velocity or magnetic field gradients. In contrast, the Stokes V profiles of Polaris and MY Pup appear qualitatively similar to the complex magnetic signatures of non-pulsating cool supergiants, possibly due to the low pulsation amplitudes of these two stars.

A fundamental question for theories of massive star formation is whether OB stars can form in isolation. We assess the contribution of any in-situ OB star formation by using 210 field OB stars in the Small Magellanic Cloud (SMC) from the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4). We search for tiny, sparse clusters around our target OB stars using cluster-finding algorithms. Employing statistical tests, we compare these observations with random-field data sets. We find that ∼5% of our target fields do show evidence of higher central stellar densities, implying the presence of small clusters. This frequency of small clusters is low and within errors, it is also consistent with the field OB population being composed entirely of runaway and walkaway stars. Assuming this small cluster fraction is real, it implies that some OB stars may form in highly isolated conditions. The low frequency could be caused by these clusters evaporating on a short timescale. However, another interpretation is that the low fraction of small clusters is observed because these form rarely, or not at all, implying a higher cluster lower-mass limit and generally consistent with a relationship between maximum stellar mass (mmax) and the cluster mass (Mcl).

In 2020, HR 6819 was reported to be a triple system containing the closest black hole to Earth. However, these results were contested, with an alternative explanation of a post-interaction binary suggested. Using new integral field spectroscopic and interferometric data, we have been able to determine the true nature of this exotic source.

When a star is rapidly rotating, it deviates from spherical symmetry causing non-uniform distributions of the surface gravity and temperature across the surface. These three-dimensional effects lead to an inclination dependence of many observable spectroscopic parameters, however this is often neglected when analyzing rapidly rotating systems. Using spamms, we generate synthetic spectra that account for the 3D geometry of the system and fit them with 1D models to investigate how much the 3D effects can change the derived stellar parameters. We show that these 3D effects can lead to observed temperature differences of thousands of kelvin for the same star viewed at different inclinations, and a systematic underestimation of the helium abundance.