Evolution Of Massive Stars Research Articles

The IACOB project

Context. Blue supergiants (BSGs) are key objects for understanding the evolution of massive stars, which play a crucial role in the evolution of galaxies. However, discrepancies between theoretical predictions and empirical observations have opened up important questions yet to be answered. Studying statistically significant and unbiased samples of these objects can help to improve the situation. Aims. We perform a homogeneous and comprehensive quantitative spectroscopic analysis of a large sample of Galactic luminous blue stars (a majority of which are BSGs) from the IACOB spectroscopic database, providing crucial parameters to refine and improve theoretical evolutionary models. Methods. We derived the projected rotational velocity (υ sin i) and macroturbulent broadening (υmac) using IACOB-BROAD, which combines Fourier transform and line-profile fitting techniques. We compared high-quality optical spectra with state-of-the-art simulations of massive star atmospheres computed with the FASTWIND code. This comparison allowed us to derive effective temperatures (Teff), surface gravities (log 𝑔), microturbulences (ξ), surface abundances of silicon and helium, and to assess the relevance of stellar winds through a wind-strength parameter (log Q). Results. We provide estimates and associated uncertainties of the above-mentioned quantities for the largest sample of Galactic luminous O9 to B5 stars spectroscopically analyzed to date, comprising 527 targets. We find a clear drop in the relative number of stars at Teff ≈ 21 kK, coinciding with a scarcity of fast rotating stars below that temperature. We speculate that this feature (roughly corresponding to B2 spectral type) might be roughly delineating the location of the empirical terminal-age main sequence in the mass range between 15 and 85 M⊙. By investigating the main characteristics of the υ sin i distribution of O stars and BSGs as a function of Teff, we propose that an efficient mechanism transporting angular momentum from the stellar core to the surface might be operating along the main sequence in the high-mass domain. We find correlations between ξ,υmac and the spectroscopic luminosity 𝓛 (defined as Teff4 / g). We also find that no more than 20% of the stars in our sample have atmospheres clearly enriched in helium, and suggest that the origin of this specific subsample might be in binary evolution. We do not find clear empirical evidence of an increase in the wind strength over the wind bi-stability region toward lower Teff.

Fits for the Convective Envelope Mass in Massive Stars

We explore the evolution of massive stars (>8 M ⊙) with 1D models and present analytical fits to the masses and binding energies of the convective portions of their envelopes. These fits are given as functions of total mass, metallicity, and surface temperature (used as a proxy for evolutionary phase). They enable the application of the two-stage common envelope (CE) formalism in rapid binary population synthesis frameworks. We estimate that the degree of orbital hardening following CE ejection spans 6 orders of magnitude and is a very strong function of the accretor mass, and, to a lesser extent, donor evolutionary phase.

Exploring extreme brightness variations in blue supergiant MACHO 80.7443.1718: Evidence for companion-driven enhanced mass loss

Context. The evolution of massive stars is dominated by interactions within binary and multiple systems. In order to accurately model this evolution, it is necessary to investigate all possible forms of an interaction in binary systems that may affect the evolution of the components. One of the “laboratories” plausible for this kind of investigation is the massive eccentric binary system MACHO 80.7443.1718 (ExtEV), which exhibits an exceptionally large amplitude of light variability close to the periastron passage of its 32.8-day orbit. Aims. We examine whether the light variability of ExtEV can be explained by a wind-wind collision (WWC) binary system model. We also critically review other models proposed to explain the light curve of ExtEV. Methods. We conducted an analysis of (i) the broadband multicolor photometry of ExtEV spanning a wide range of wavelengths from the ultraviolet to near-infrared, (ii) the time-series space photometry from the Transiting Exoplanet Survey Satellite (TESS), (iii) ground-based Johnson UBV photometry, and (iv) time-series high-resolution spectroscopy. To derive the parameters of the primary component of the system, we fit the spectral energy distribution (SED) and calculated evolutionary models of massive stars that included mass loss. Using radial-velocity data, we determined the spectroscopic parameters of the system. We also fit an analytical model of light variations to the TESS light curve of ExtEV. Results. The ExtEV system exhibits an infrared excess, indicating an increased mass-loss rate. The system does not match the characteristics of B[e] stars, however. We rule out the possibility of the presence of a Keplerian disk around the primary component. We also argue that the scenario with periodic Roche-lobe overflow at periastron may not be consistent with the observations of ExtEV. Analysis of the SED suggests that the primary component has a radius of about 30 R⊙ and a luminosity of ∼6.6 × 105 L⊙. With the analysis of the radial-velocity data, we refine the orbital parameters of ExtEV and find evidence for the presence of a tertiary component in the system. Using evolutionary models we demonstrate that the primary component’s mass is between 25 and 45 M⊙. We successfully reproduced the light curve of ExtEV with our analytical model, showing that the dominant processes shaping its light curve can be attributed to the atmospheric eclipse and light scattered in the WWC cone. We also estimate the primary’s mass loss rate due to stellar wind for 4.5 × 10−5 M⊙ yr−1. Conclusions. ExtEV is most likely not an extreme eccentric ellipsoidal variable, but rather an exceptional WWC binary system. The mass loss rate we derived exceeds theoretical predictions by up to two orders of magnitude. This implies that the wind in the system is likely enhanced by tidal interactions, rotation, and possibly also tidally excited oscillations. Therefore, ExtEV represents a rare evolutionary phase of a binary system that may help to understand the role of companion-driven enhanced mass loss in the evolution of massive binary systems.

Investigating episodic mass loss in evolved massive stars

Context. Episodic mass-loss events such as giant eruptions in luminous blue variables or pre-supernova eruptions in red supergiants drastically alter the evolutionary path of a massive star, resulting in a rich and complex circumstellar environment and IR excess. However, the incidence of these events, and hence their importance in massive star evolution, remains unknown. Aims. The ASSESS project (Episodic Mass Loss in Evolved Massive Stars: Key to Understanding the Explosive early Universe) aims to determine the role of episodic mass loss in the evolution of massive stars. As a first step, we constructed a catalog of spectroscopically identified dusty, evolved massive stars in ten southern galaxies for which Spitzer point-source catalogs are available. The resulting catalog may be used to identify stars that have undergone an episodic mass-loss event. The target galaxies span a range of metallicities, Z = 0.06–1.6 Z⊙, allowing for the investigation of a potential metallicity dependence. Methods. We conducted multi-object spectroscopy of dusty massive star candidates in ten target galaxies using the Very Large Telescope. We obtained 763 spectra from WLM, NGC 55, NGC 247, NGC 253, NGC 300, NGC 1313, NGC 3109, Sextans A, M83, and NGC 7793. The targets were selected using their Spitzer photometry, by prioritizing targets with a strong IR excess, which indicates the presence of hot dust. We determined a spectral classification for each target. Additionally, we used archival images from the Hubble Space Telescope (HST), available for 150 of our targets, to provide a visual classification for 80 targets, as a star, cluster, or galaxy. Results. We provide a catalog of 541 spectroscopically classified sources that includes 185 massive stars, of which 154 are newly classified massive stars. The catalog contains 129 red supergiants, 27 blue supergiants, 10 yellow supergiants, 4 luminous blue variable candidates, 7 supergiant B[e] stars, and 8 emission-line objects. Evidence for circumstellar dust is found in 24% of these massive stars, based on their IR colors. We report a success rate of 28% for identifying massive stars from our observed spectra, while the average success rate of our priority system for selecting evolved massive stars is 36%. Additionally, the catalog contains 21 background galaxies (including active galactic nuclei and quasars), 10 carbon stars, and 99 H II regions. We measured the line ratios [N II]/Hα and [S II]/Hα for 76 H II regions and 36 other spectra with nebular emission lines, thereby identifying eight sources with shocked emission. Conclusions. We present the largest catalog of evolved massive stars and in particular of red supergiants in nearby galaxies at low Z beyond the Local Group. The brightest and reddest of these are candidates for episodic mass loss. The fraction of dusty massive stars observed with respect to the initial selection is ~30%. We expect this catalog to trigger follow-up studies and pave the way for a comprehensive study of the eruptive late stages of massive star evolution in the era of the James Webb Space Telescope and the new survey telescopes (e.g., the Euclid mission, Nancy Grace Roman Space Telescope, and Vera C. Rubin Observatory).

β Cephei Pulsators in Eclipsing Binaries Observed with TESS

The combined strength of asteroseismology and empirical stellar basic parameter determinations for in-depth asteroseismic analysis of massive pulsators in eclipsing binaries shows great potential for treating the challenging and mysterious discrepancies between observations and models of stellar structure and the evolution of massive stars. This paper compiles a comprehensive list of massive pulsators in eclipsing binary systems observed with TESS. The TESS light curves and discrete Fourier transforms of a sample of 8055 stars of spectral type B0–B3 were examined for eclipses and stellar pulsations, and the ephemerides of the resulting subsample of massive pulsators in eclipsing binaries were computed. This subsample was also crossmatched with existing catalogs of massive pulsators. Until now, fewer than 30 β Cephei pulsators in eclipsing binaries have been reported in the literature. Here we announce a total of 78 pulsators of the β Cephei type in eclipsing binaries, 59 of which are new discoveries. Forty-three are recognized as definite, and 35 are candidate pulsators. Our sample of pulsating massive stars in eclipsing binaries allows for future asteroseismic modeling to better understand the internal mixing profile and to resolve the mass discrepancy in massive stars. We have already started follow-up work on some of the most interesting candidates.

Nuclear structure properties and weak interaction rates of even–even Fe isotopes

Nuclear structure properties and weak interaction rates of neutron-rich even–even iron (Fe) isotopes (A = 50 – 70) are investigated using the Interacting Boson Model-1 (IBM-1) and the proton–neutron Quasiparticle Random Phase Approximation (pn-QRPA) model. The IBM-1 is used for the calculation of energy levels and the B(E2) values of neutron-rich Fe isotopes. Later their geometry was predicted within the potential energy formalism of the IBM-1 model. Weak interaction rates on neutron-rich nuclei are needed for the modeling and simulation of presupernova evolution of massive stars. In the current study, we investigate the possible effect of nuclear deformation on stellar rates of even–even Fe isotopes. The pn-QRPA model is applied to calculate the weak interaction rates of selected Fe isotopes using three different values of deformation parameter. It is noted that, in general, bigger deformation values led to smaller total strength and larger centroid values of the resulting Gamow–Teller distributions. This later translated to smaller computed weak interaction rates. The current finding warrants further investigation before it may be generalized. The reported stellar rates are up to 4 orders of magnitude smaller than previous calculations and may bear astrophysical significance.

Investigating the equation-of-state, stability and mass–radius relationship of anisotropic and massive neutron stars embedded in f(R,T) modified gravity

In this study, our main focus is to investigate the mass–radius relation and several important properties of massive neutron stars to realize the nature, behavior and evolution of these kinds of compact objects at present time. Also, we want to understand the equation-of-state of the core nuclear matter precisely with their stable equilibrium configuration. We have chosen a few massive binary pulsars and investigated on maximum attainable mass and lowest possible radius of them. We have considered Einstein–Hilbert action as [Formula: see text], where [Formula: see text] is the Ricci scalar and [Formula: see text], the trace of energy–momentum tensor with [Formula: see text] as the coupling parameter and also used modified TOV equations. The interior space-time of the spherical neutron star is matched to the exterior Schwarzschild line element at the surface of the star. The [Formula: see text]–[Formula: see text] curve can predict the maximum achievable mass is about [Formula: see text] with lowest possible radius of around [Formula: see text] km for the massive compact stars under stable equilibrium. We have obtained a clear picture of structural evolution of massive neutron stars through accelerating space-time and can put some constraints on several quantities related to them. It is depicted from our present investigation that all the derived outcomes are compatible with physically adopted regimes which reveal the viability of our current model in the context of [Formula: see text] modified gravity.

The Physics of Core-Collapse Supernovae: Explosion Mechanism and Explosive Nucleosynthesis

Recent developments in multi-dimensional simulations of core-collapse supernovae have considerably improved our understanding of this complex phenomenon. In addition to that, one-dimensional (1D) studies have been employed to study the explosion mechanism and its causal connection to the pre-collapse structure of the star, as well as to explore the vast parameter space of supernovae. Nonetheless, many uncertainties still affect the late stages of the evolution of massive stars, their collapse, and the subsequent shock propagation. In this review, we will briefly summarize the state-of-the-art of both 1D and 3D simulations and how they can be employed to study the evolution of massive stars, supernova explosions, and shock propagation, focusing on the uncertainties that affect each of these phases. Finally, we will illustrate the typical nucleosynthesis products that emerge from the explosion.

A Spectroscopic Hunt for Post-red Supergiants in the Large Magellanic Cloud. I. Preliminary Results

Yellow supergiants (YSGs) are rare and poorly understood, and studying them is critical to constraining massive star evolution. We obtained flux-calibrated Magellan Inamori Kyocera Echelle high-resolution spectra of 40 YSGs in the Large Magellanic Cloud (LMC); this sample likely contains post-red supergiants. Fitting these data with ATLAS9 model atmospheres, we determined fundamental parameters for these stars. We measure the first spectroscopic luminosities for YSGs above 20 M ⊙, providing us a novel probe of the luminosity-to-mass ratio. Many stars in our sample appear to have anomalously high surface gravities, despite being confirmed LMC supergiants. We manually inspected our data finding evidence for binary companions and ongoing mass loss. Our work demonstrates the valuable role of high-resolution spectroscopy in interpreting the evolutionary status of cool supergiants.

Re-Examination of the Effect of Pairing Gaps on Gamow–Teller Strength Distributions and β-Decay Rates

β-decay is one of the key factors for understanding the r-process and evolution of massive stars. The Gamow–Teller (GT) transitions drive the β-decay process. We employ the proton–neutron quasiparticle random phase approximation (pn-QRPA) model to calculate terrestrial and stellar β-decay rates for 50 top-ranked nuclei possessing astrophysical significance according to a recent survey. The model parameters of the pn-QRPA model affect the predicted results of β-decay. The current study investigates the effect of nucleon–nucleon pairing gaps on charge-changing transitions and the associated β decay rates. Three different values of pairing gaps, namely TF, 3TF, and 5TF, were used in our investigation. It was concluded that both GT strength distributions and half-lives are sensitive to pairing gap values. The 3TF pairing gap scheme, in our chosen nuclear model, resulted in the best prediction with around 80% of the calculated half-lives within a factor 10 of the measured ones. The 3TF pairing scheme also led to the calculation of the biggest β-decay rates in stellar matter.

Observational Properties of 155 β Cephei Pulsating Variable Stars

β Cephei pulsating variable (BCEP) stars are the most massive pulsating variable stars in the main sequence, exhibiting both p- and g-mode pulsations. In this study, we identified 155 BCEP stars or candidates using data from TESS and Gaia, of which 83 were first confirmed as BCEP stars. They have visual magnitudes ranging from 8 to 12 mag and effective temperatures between approximately 20,000 and 30,000 K, while the parallaxes of most targets are between 0.2 and 0.6 mas. The study indicates that these BCEP stars have pulsation periods ranging from 0.06 to 0.31 days, with amplitudes ranging from 0.1 to 55.8 mmag in the TESS band. Additionally, the number of BCEP stars increases as the pulsation amplitude decreases. These targets align with the distribution region of BCEP stars in the luminosity–period (L–P) and temperature–period diagrams. We have updated the L–P relation of BCEP stars. The Hertzsprung–Russell diagram indicates that these targets are in the main-sequence evolutionary phase, with masses ranging from 7 to 20 M ⊙ and luminosities between 2800 and 71,000 L ⊙. They are almost in the theoretical instability region of BCEP stars but as previously reported, this region at the low-mass end (red) is not filled. The distribution of the pulsation constant indicates that the dominant pulsation periods of BCEP stars consist mainly of low-order p-mode pulsations with a high proportion of radial fundamental modes. These BCEP stars are excellent objects for enhancing our understanding of the structure and evolution of massive stars through asteroseismology.

The dragon-II simulations – III. Compact binary mergers in clusters with up to 1 million stars: mass, spin, eccentricity, merger rate, and pair instability supernovae rate

ABSTRACT Compact binary mergers forming in star clusters may exhibit distinctive features that can be used to identify them among observed gravitational-wave sources. Such features likely depend on the host cluster structure and the physics of massive star evolution. Here, we dissect the population of compact binary mergers in the dragon-II simulation data base, a suite of 19 direct N-body models representing dense star clusters with up to 106 stars and $\lt 33~{{\ \rm per\ cent}}$ of stars in primordial binaries. We find a substantial population of black hole binary (BBH) mergers, some of them involving an intermediate-mass BH (IMBH), and a handful mergers involving a stellar BH and either a neutron star (NS) or a white dwarf (WD). Primordial binary mergers, $\sim 30~{{\ \rm per\ cent}}$ of the whole population, dominate ejected mergers. Dynamical mergers, instead, dominate the population of in-cluster mergers and are systematically heavier than primordial ones. Around 20 per cent of dragon-II mergers are eccentric in the Laser Interferometer Space Antenna (LISA) band and 5 per cent in the LIGO band. We infer a mean cosmic merger rate of $\mathcal {R}\sim 30(4.4)(1.2)$ yr−1 Gpc−3 for BBHs, NS–BH, and WD–BH binary mergers, respectively, and discuss the prospects for multimessenger detection of WD–BH binaries with LISA. We model the rate of pair-instability supernovae (PISNe) in star clusters and find that surveys with a limiting magnitude mbol = 25 can detect ∼1–15 yr−1 PISNe. Comparing these estimates with future observations could help to pin down the impact of massive star evolution on the mass spectrum of compact stellar objects in star clusters.

The drivers of massive star evolution

An elusive population of stripped binary stars is finally revealed

Stellar-wind variability in Cygnus X-1 from high-resolution excess variance spectroscopy with Chandra

Context. Stellar winds of massive stars are known to be driven by line absorption of UV photons, a mechanism that is prone to instabilities, causing the wind to be clumpy. The clumpy structure hampers wind mass-loss estimates, limiting our understanding of massive star evolution. The wind structure also impacts accretion in high-mass X-ray binary (HMXB) systems. Aims. We aim to analyse the wavelength-dependent variability of X-ray absorption in the wind to study its structure. Such an approach is possible in HMXBs, where the compact object serves as an X-ray backlight. We probe different parts of the wind by analysing data taken at superior and inferior conjunctions. Methods. We applied excess variance spectroscopy to study the wavelength-dependent soft (2–14 Å) X-ray variability of the HMXB Cygnus X-1 in the hard spectral state. Excess variance spectroscopy quantifies the variability of an object above the statistical noise as a function of wavelength, which allows us to study the variability of individual spectral lines. This technique was applied to high-resolution gratings spectra provided by Chandra, accounting for various systematic effects. The frequency dependence is investigated by changing the time binning. Results. The strong orbital phase dependence we observe in the excess variance is consistent with column-density variations predicted by a simple model for a clumpy wind. We identify spikes of increased variability with spectral features found by previous spectroscopic analyses of the same data set, most notably from silicon in over-dense clumps in the wind. In the silicon line region, the variability power is redistributed towards lower frequencies, hinting at increased line variability in large clumps. In prospect of the microcalorimetry missions that are scheduled to launch within the next decade, excess variance spectra present a promising approach to constraining the wind structure, especially if accompanied by models that consider changing ionisation.

A probable Keplerian disk feeding an optically revealed massive young star.

The canonical picture of star formation involves disk-mediated accretion, with Keplerian accretion disks and associated bipolar jets primarily observed in nearby, low-mass young stellar objects (YSOs). Recently, rotating gaseous structures and Keplerian disks have been detected around several massive (M > 8 M⊙) YSOs (MYSOs)1-4, including several disk-jet systems5-7. All the known MYSO systems are in the Milky Way, and all are embedded in their natal material. Here we report the detection of a rotating gaseous structure around an extragalactic MYSO in the Large Magellanic Cloud. The gas motion indicates that there is a radial flow of material falling from larger scales onto a central disk-like structure. The latter exhibits signs of Keplerian rotation, so that there is a rotating toroid feeding an accretion disk and thus the growth of the central star. The system is in almost all aspects comparable to Milky Way high-mass YSOs accreting gas from a Keplerian disk. The key difference between this source and its Galactic counterparts is that it is optically revealed rather than being deeply embedded in its natal material as is expected of such a massive young star. We suggest that this is the consequence of the star having formed in a low-metallicity and low-dust content environment. Thus, these results provide important constraints for models of the formation and evolution of massive stars and their circumstellar disks.

On the plerionic rectangular supernova remnants of static progenitors

ABSTRACT Pulsar wind nebulae are a possible final stage of the circumstellar evolution of massive stars, where a fast-rotating, magnetized neutron star produces a powerful wind that interacts with the supernova ejecta. The shape of these so-called plerionic supernova remnants is influenced by the distribution of circumstellar matter at the time of the explosion, itself impacted by the magnetic field of the ambient medium, responsible for the expansion of the circumstellar bubble of the progenitor star. To understand the effects of magnetization on the circumstellar medium and resulting pulsar nebulae, we conduct 2D magnetohydrodynamic simulations. Our models explore the impact of the interstellar medium’s (ISM) magnetic field on the morphology of a supernova remnant and pulsar wind nebula that develop in the circumstellar medium of massive star progenitor in the warm phase of the Milky Way’s ISM. Our simulations reveal that the jet-like structures formed on both sides perpendicularly to the equatorial plane of the pulsar, creating complex radio synthetic synchrotron emissions. This morphology is characterized by a rectangular-like remnant, which is typical of the circumstellar medium of massive stars in a magnetized medium, along with the appearance of a spinning top structure within the projected rectangle. We suggest that this mechanism may be partially responsible for the complex morphologies observed in pulsar wind nebulae that do not conform to the typical torus/jet or bow shock/tail shapes observed in most cases.

Evolutionary tracks of massive stars with different rotation and metallicity in neutrino H–R diagram

ABSTRACT Neutrino losses play a crucial role in the evolution of massive stars. We study the neutrino luminosity of stars ranging from 20 to 90 $\rm M_{\odot }$ from Zero Age Main Sequence to Fe Core Collapse (FeCC) with different rotation and metallicity in a neutrino Hertzsprung–Russell diagram. In our simulations, we consider $\rm \omega /\omega _{crit}= 0$ and 0.7 to represent non-rotation and high rotation, respectively, and set the metallicities to 0.014, 0.001, and 0.0001. During hydrogen burning stages, neutrino luminosity primarily originates from CNO cycle, and increases with higher stellar mass while decreasing with increasing metallicity. For the high metallicity models (Z = 0.014) during the helium burning stage, the reduction of the hydrogen envelope caused by a larger mass loss rate leads to a gradual decrease in neutrino luminosity. The rapid rotation results in extra mixing inside massive stars, which increases the neutrino luminosity during main sequence, while decreases the neutrino luminosity during helium burning phase. Simultaneously, the rapid rotation also increases CO core mass, which enhances the neutrino luminosity during C and O burning phase. We also investigate the effect of neutrino magnetic moment (NMM) on the massive stars. We find that the energy loss caused by the NMM does not have effects on the evolutionary destiny of massive stars, and it does not significant change the compactness at the time of FeCC.

Exploring the Red Supergiant wind kink

Red supergiants (RSG) are key objects in studying the evolution of massive stars and their endpoints, but uncertainties related to their underlying mass-loss mechanism have stood in the way of an appropriate framework for massive star evolution thus far. In this work, we analyse a recently uncovered empirical mass-loss “kink” feature and we highlight its similarity to hot star radiation-driven wind models and observations at the optically thin-to-thick transition point. We motivate a new RSG mass-loss prescription that depends on the Eddington factor, Γ, (including both a steep luminosity, L, dependence and an inverse steep mass, Mcur, dependence). We subsequently implement this new RSG mass-loss prescription in the stellar evolution code MESA. We find that our physically motivated mass-loss behaviour naturally reproduces the Humphreys-Davidson limit without the need for any ad hoc tweaks. It also resolves the RSG supernova “problem”. We argue that a universal behaviour that is seen for radiation-driven winds across the HR diagram, independent of the exact source of opacity, is a key feature of the evolution of the most massive stars.

Left Ringing: Betelgeuse Illuminates the Connection between Convective Outbursts, Mode Switching, and Mass Ejection in Red Supergiants

Betelgeuse, the nearest red supergiant, dimmed to an unprecedented level in early 2020. The star emerged from this “Great Dimming” episode with its typical, roughly 400 days pulsation cycle halved, and a new dominant period of around 200 days. The dimming event has been attributed to a surface mass ejection, in which rising material drove shocks through the stellar atmosphere and expelled some material, partially obscuring the star as it formed molecules and dust. In this paper, we use hydrodynamic simulations to reveal the connections between Betelgeuse's vigorously convective envelope, the surface mass ejection, and the pulsation mode switching that ensued. An anomalously hot convective plume, generated rarely but naturally in the star's turbulent envelope, can rise and break free from the surface, powering an upwelling that becomes the surface mass ejection. The rising plume also breaks the phase coherence of the star's pulsation, causing the surface to keep expanding even as the deeper layers contract. This drives a switch from the 400 days fundamental mode of pulsation, in which the whole star expands and contracts synchronously, to the 200 days first overtone, where a radial node separates the interior and exterior of the envelope moving in opposite phase. We predict that the star's convective motions will damp the overtone oscillation and Betelgeuse will return to its previous, 400 days fundamental-mode pulsation in the next 5–10 yr. With its resolved surface and unprecedentedly detailed characterization, Betelgeuse opens a window to episodic surface mass ejection in the late-stage evolution of massive stars.

A Metallicity Dependence on the Occurrence of Core-collapse Supernovae

Core-collapse supernovae (CCSNe) are widely accepted to be caused by the explosive death of massive stars with initial masses ≳8 M ⊙. There is, however, a comparatively poor understanding of how properties of the progenitors—mass, metallicity, multiplicity, rotation, etc.—manifest in the resultant CCSN population. Here, we present a minimally biased sample of nearby CCSNe from the All-Sky Automated Survey for Supernovae survey whose host galaxies were observed with integral-field spectroscopy using MUSE at the Very Large Telescope. This data set allows us to analyze the explosion sites of CCSNe within the context of global star formation properties across the host galaxies. We show that the CCSN explosion site oxygen abundance distribution is offset to lower values than the overall H ii region abundance distribution within the host galaxies. We further split the sample at dex and show that within the subsample of low-metallicity host galaxies, the CCSNe unbiasedly trace the star formation with respect to oxygen abundance, while for the subsample of higher-metallicity host galaxies, they preferentially occur in lower-abundance star-forming regions. We estimate the occurrence of CCSNe as a function of oxygen abundance per unit star formation and show that there is a strong decrease as abundance increases. Such a strong and quantified metallicity dependence on CCSN production has not been shown before. Finally, we discuss possible explanations for our result and show that each of these has strong implications not only for our understanding of CCSNe and massive star evolution but also for star formation and galaxy evolution.