Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open questions include the shape and extent of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries. In this contribution, we present the results of a large and systematic analysis of a sample of over 80 single-lined O-type spectroscopic binaries (SB1s) in the Milky Way and in the Large Magellanic Cloud (LMC). We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries.

Over the last two decades there have been considerable advances in modelling the spectra of massive stars and supernovae (SNe). Despite this progress, there are still numerous uncertainties that affect the accuracy of models. For massive stars, convection, instabilities, clumping, and our inability to model stellar winds self-consistently likely introduce systematic errors into our analyses. For SNe, and particularly for core-collapse SNe, departures from spherical symmetry strongly affect observed spectra and need to be taken into account. There are also issues with clumping, and mixing processes (both in the progenitor and the SN explosion) that need to be resolved. For both massive stars and SNe, the accuracy and availability of atomic data continues to be an ongoing issue influencing analyses.

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.

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.

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.

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.

The jittering jets explosion mechanism of core collapse supernovae proposes that stochastic convective motion, amplified by instabilities, results in intermittent accretion disk that launches jittering jets that explode the star. We conduct one-dimensional simulations of a wide range of stellar masses and show that the convective motion in the pre-collapse stellar core, when scaled to corresponding three-dimensional simulations, supply sufficient angular momentum fluctuation to form the intermittent accretion disk the launches the jittering jets. The resulting neutron star masses are consistent with observations.

Supernova properties in radio strongly depend on their circumstellar environment and they are an important probe to investigate the mass loss of supernova progenitors. Recently, core-collapse supernova observations in radio have been assembled and the rise time and peak luminosity distribution of core-collapse supernovae in radio has been obtained. In this talk, we will discuss the constraints on the mass-loss prescriptions of red supergiants obtained from the assembled radio properties of Type II supernovae. We take a couple of mass-loss prescriptions for red supergiants, calculate the rise time and peak luminosity distribution based on them, and compare the results with the observed distribution. We found that the widely spread radio rise time and peak luminosity distribution of Type II supernovae can only be explained by mass-loss prescriptions having strong dependence on the luminosity. Red supergiant mass-loss prescriptions should have steep luminosity dependence in the supernova progenitor range.

When a supernova shockwave launched deep inside the star exits the surface, it probes the circumstellar medium established by prior mass loss from the pre supernova star. The bright electromagnetic display accompanying the shock breakout is influenced by the properties of the star and scripts the history of the stellar mass loss. We investigate with MESA and STELLA codes the radiative display resulting from a set of progenitors that we evolved to core collapse. We simulate with different internal convective overshoot and compositional mixing and two sets of mass loss schema, one the standard “Dutch” scheme and another, an enhanced, episodic mass loss at a late stage. Shock breakout from the star shows double peaked bolometric light curves for the Dutch wind, as well as high velocity ejecta accelerated during shock breakout. We contrast the breakout flash from an optically thick CSM with that of the rarified medium.

The nuclear equation-of-state (EOS) describing newly formed proto-neutron stars (PNSs) in core-collapse supernovae (CCSNe) is yet uncertain, and varying its prescription affects multimessenger signatures in CCSN simulations. Focusing on the gravitational wave (GW) signal, we demonstrate the effect of varying parameter values in the EOS. We conclude that an especially important parameter is the effective mass of nucleons which affect thermal properties and subsequently the PNS compactness, regulating the GW signal in both amplitude and frequency. By radially decomposing the GW emission, we show where in the PNS the GWs originate from.

The combination of strong magnetic fields and fast rotation is often invoked as a characteristic of the central engine for outstanding sources such as GRBs, hypernovae, and superluminous supernovae. However, the actual properties of the magnetic field during the collapse of the stellar progenitor are very uncertain, since they depend on the evolution of the star and can be affected by complex dynamo processes occurring in the central proto-neutron star. Using 3D relativistic MHD models we show that higher-order multipolar fields can lead to the onset of a supernova, although they tend to produce less energetic explosions and less collimated outflows. Quadrupolar fields efficiently extract angular momentum from the central core, but the rotational energy is partly stored in the equatorial regions, rather than powering up the polar outflows. Finally, our results show a strong magnetic quenching of the hydrodynamic non-axisymmetric instabilities that are associated to the emission of GWs.

The convection-enhanced neutrino-driven supernova engine’s success in explaining a myriad of supernova properties has set it as the standard engine behind supernova. However, due to the success of rotationally-powered engines in explaining astrophysical transients like gamma-ray bursts, these engines have been revived as possible drivers of normal supernovae, competing with this standard engine. In this paper, these competing engines, and the constraints placed by compact remnant observations on these engines, are reviewed. We find that, with these constraints, such rotationally-powered engines can explain less than 1% of the current supernova remnants. In addition, we find that the remnant mass distribution can be used to constrain properties of the convection-enhanced neutrino-driven engine, helping astronomers understand the nature of convection in this engine.

Motivated by their role as the direct or indirect source of many of the elements in the Universe, numerical modeling of core collapse supernovae began more than five decades ago. Progress toward ascertaining the explosion mechanism(s) has been realized through increasingly sophisticated models, as physics and dimensionality have been added, as physics and numerical modeling have improved, and as the leading computational resources available to modelers have become far more capable. The past five to ten years have witnessed the emergence of a consensus across the core collapse supernova modeling community that had not existed in the four decades prior. For the majority of progenitors – i.e., slowly rotating progenitors – the efficacy of the delayed shock mechanism, where the stalled supernova shock wave is revived by neutrino heating by neutrinos emanating from the proto-neutron star, has been demonstrated by all core collapse supernova modeling groups, across progenitor mass and metallicity. With this momentum, and now with a far deeper understanding of the dynamics of these events, the path forward is clear. While much progress has been made, much work remains to be done, but at this time we have every reason to be optimistic we are on track to answer one of the most important outstanding questions in astrophysics: How do massive stars end their lives?

We have detected 27 new supernova remnants (SNRs) using a new data release of the GLEAM survey from the Murchison Widefield Array telescope, including the lowest surface brightness SNR ever detected, G 0.1 – 9.7. Our method uses spectral fitting to the radio continuum to derive spectral indices for 26/27 candidates, and our low-frequency observations probe a steeper spectrum population than previously discovered. None of the candidates have coincident WISE mid-IR emission, further showing that the emission is non-thermal. Using pulsar associations we derive physical properties for six candidate SNRs, finding G 0.1 – 9.7 may be younger than 10 kyr. Sixty per cent of the candidates subtend areas larger than 0.2 deg2 on the sky, compared to < 25% of previously detected SNRs. We also make the first detection of two SNRs in the Galactic longitude range 220°–240°.

We examined the latest data release from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey covering 345° < l < 60° and 180° < l < 240°, using these data and that of the Widefield Infrared Survey Explorer to follow up proposed candidate Supernova Remnant (SNR) from other sources. Of the 101 candidates proposed in the region, we are able to definitively confirm ten as SNRs, tentatively confirm two as SNRs, and reclassify five as H ii regions. A further two are detectable in our images but difficult to classify; the remaining 82 are undetectable in these data. We also investigated the 18 unclassified Multi-Array Galactic Plane Imaging Survey (MAGPIS) candidate SNRs, newly confirming three as SNRs, reclassifying two as H ii regions, and exploring the unusual spectra and morphology of two others.

Theories that attempt to unify the four fundamental interactions and alternative theories of gravity predict time and/or spatial variation of the fundamental constants of nature. Different versions of these theories predict different behaviours for these variations. As a consequence, experimental and observational bounds are an important tool to check the validity of such proposals. In this paper, we review constraints on the possible variation of the fundamental constants from astronomical observations and geophysical experiments designed to test the constancy of the fundamental constants of nature over different timescales. We also focus on the limits that can be obtained from white dwarfs, which can constrain the variation of the constants with the gravitational potential.

SN1991bg-like supernovae are a distinct subclass of thermonuclear Type Ia supernovae (SNe Ia). Their spectral and photometric peculiarities indicate that their progenitors and explosion mechanisms differ from ‘normal’ SNe Ia. One method of determining information about supernova progenitors we cannot directly observe is to observe the stellar population adjacent to the apparent supernova explosion site to infer the distribution of stellar population ages and metallicities. We obtain integral field observations and analyse the spectra extracted from regions of projected radius $\sim\,\text{kpc}$ about the apparent SN explosion site for 11 91bg-like SNe in both early- and late-type galaxies. We utilise full-spectrum spectral fitting to determine the ages and metallicities of the stellar population within the aperture. We find that the majority of the stellar populations that hosted 91bg-like supernovae have little recent star formation. The ages of the stellar populations suggest that that 91bg-like SN progenitors explode after delay times of >6 Gyr, much longer than the typical delay time of normal SNe Ia, which peaks at $\sim$1 Gyr.

Supernovae provide environments with strong links to laboratory astrophysics. Diverse physical processes spanning from hot gas and semi-relativistic particles down to cold dusty clumps require extensive atomic data and understanding of processes across different physical regimes. The current status of modelling and analyzing supernova spectra is reviewed, with focus on recent results for diagnosing the production of oxygen and nickel.

We present results of a supernova lightcurve population synthesis, predicting the range of possible supernova lightcurves arising from a population of progenitor stars that include interacting binary systems. We show that the known diversity of supernova lightcurves can be interpreted as arising from binary interactions. Given detailed models of the progenitor stars, we are able to the determine what parameters within these stars determine the shape of their supernova lightcurve. The primary factors are the mass of supernova ejecta and the mass of hydrogen in the final progenitor. We find that there is a continuum of lightcurve behaviour from type IIP, IIL, to IIb supernovae related to the range of hydrogen and ejecta masses. Most type IIb supernovae arise from a relatively narrow range of initial masses from 10 to 15 M⊙. We also find a few distinct lightcurves that are the result of stellar mergers.