1. Introduction
Red Supergiants (RSGs) are cool, evolved massive stars with initial masses of 8–30 M
$_\unicode{x2299}$, representing the final evolutionary stage before they end their lives in a Type II supernova (SN) explosion (Ekström et al. Reference Ekström, Georgy and Eggenberger2012; Levesque Reference Levesque2017). However, the apparent lack of luminous RSGs detected as SN progenitors has sparked an ongoing debate over the fate of these stars, known as the so-called “RSG problem”. While several studies support its existence (Li et al. Reference Li, Van Dyk and Filippenko2006; Smartt et al. Reference Smartt, Eldridge, Crockett and Maund2009; Smartt Reference Smartt2015; Fang et al. Reference Fang, Moriya and Maeda2025), others argue that this problem is not statistically significant and may stem from the methodologies used to derive SN progenitor properties (Davies and Beasor Reference Davies and Beasor2020; Beasor et al. Reference Beasor, Smith and Jencson2024; Strotjohann et al. Reference Strotjohann, Ofek and Gal-Yam2024). Alternative explanations include the failed-SN scenario, involving direct collapse into a black hole (De et al. Reference De, MacLeod and Jencson2024), or post-RSG evolution (i.e., Yellow Hypergiant (YHG) phase; de Jager Reference de Jager1998; Gordon and Humphreys Reference Gordon and Humphreys2019) in which the star loses its outer envelope. Nevertheless, recent empirical studies (e.g., Beasor et al. Reference Beasor, Davies and Smith2020; Decin et al. Reference Decin, Richards, Marchant and Sana2024; Antoniadis et al. Reference Antoniadis, Bonanos and de Wit2024) have derived weak winds that are insufficient to produce post-RSG objects. The detection of vast episodic mass-loss events in luminous RSGs (e.g., NML Cyg and VY CMa; Richards et al. Reference Richards, Yates and Cohen1996; Singh et al. Reference Singh, Richards, Humphreys, Decin and Ziurys2023) provides an alternative mechanism to the steady winds. The physical mechanism driving these episodic events remains unknown, and their impact on stellar evolution is unclear, as the significance of episodic mass loss varies greatly among RSGs (Humphreys and Jones Reference Humphreys and Jones2022). Eruptions and large mass ejections due to instabilities (Glatzel and Kiriakidis Reference Glatzel and Kiriakidis1993) are expected at luminosities close to the Humphreys-Davidson limit (
$\log(L/{L_\unicode{x2299}})=5.5-5.8$; Humphreys and Davidson Reference Humphreys and Davidson1979; Davies and Beasor Reference Davies and Beasor2020; McDonald et al. Reference McDonald, Davies and Beasor2022). Alternatively, binary interactions can trigger episodic mass loss (Landri and Pejcha Reference Landri and Pejcha2024; MacLeod et al. Reference MacLeod, Blunt and De Rosa2024). In general, binarity adds further complexity to the fate of RSGs as it can drive stellar evolution by initiating a common-envelope phase (CE), stripping their envelope, or leading to a merger (Sana et al. Reference Sana, de Mink and de Koter2012; Ivanova et al. Reference Ivanova, Justham and Chen2013). The scarcity of very luminous RSGs, the presence of large amounts of circumstellar material, and their inaccurate distances (Humphreys Reference Humphreys2019) hinder accurate measurements of their stellar properties and whether they are single or binary systems, complicating the assessment of their evolutionary stage and ultimate fate.
WOH G64 (IRAS 04553-6825) (Westerlund et al. Reference Westerlund, Olander and Hedin1981) stands out as the most extreme RSG in the Large Magellanic Cloud (LMC), on the brink of the luminosity limit (
$\log(L/L_\unicode{x2299})=5.45$; Ohnaka et al. Reference Ohnaka, Driebe, Hofmann, Weigelt and Wittkowski2008), being one of the coolest (M5-7e; Elias et al. Reference Elias, Frogel and Schwering1986; Levesque et al. Reference Levesque, Massey, Plez and Olsen2009), largest (
$R=1540$ R
$_\unicode{x2299}$; Levesque et al. Reference Levesque, Massey, Plez and Olsen2009) and the record-holder in terms of mass-loss rate (
$\dot{M}\gt10^{-4}$ M
$_\unicode{x2299}$ yr
$^{-1}$; Beasor and Smith Reference Beasor and Smith2022; Antoniadis et al. Reference Antoniadis, Bonanos and de Wit2024). Interferometric observations in the near and mid-IR revealed elongated emission, which was attributed to a dusty, disk-like structure of 3-9 M
$_\unicode{x2299}$ (Ohnaka et al. Reference Ohnaka, Driebe, Hofmann, Weigelt and Wittkowski2008; Ohnaka et al. Reference Ohnaka, Hofmann and Weigelt2024). It is also an OH/IR star with strong multi-component maser emission, possibly due to an external disk, two expanding dust shells, or bipolar outflows (van Loon et al. Reference van Loon, Zijlstra, Bujarrabal and Nyman2001; Marshall et al. Reference Marshall, van Loon and Matsuura2004).
2. Results
In this talk, we present the results from Munoz-Sanchez et al. (Reference Munoz-Sanchez, Kalitsounaki and de Wit2024). We combined new observations with archival photometric and spectroscopic data to analyze WOH G64.
The optical light curve of WOH G64, spanning from 1992 to the present, clearly shows two distinct phases, separated by a transition in 2014 (Fig. 1). WOH G64 exhibits a sinusoidal pattern from the 1990s until 2013, but between mid-2013 and mid-2014, the star became significantly bluer, with a color change of
$\Delta(V-I) \approx 1.8$ mag. The change was caused exclusively by a sudden brightening in the V-band within a year, while the I-band remained approximately constant. Since 2014, WOH G64 has exhibited irregular variability, including a recent fading of 2 mag over 10 months, reaching its faintest observed I-band state to date. Mid-infrared (IR) brightness before and after the transition does not show significant changes, except for three brightenings of
$\Delta W1 \approx 0.5$ mag with a potential periodicity of
$\sim$ 750 days.
Optical light curves of WOH G64 spanning
$\gt$30 years. Each dataset is labeled with the corresponding survey and filter. The shadowed patch between 2013.5 and 2014.5 constrains the epoch of the transition. The vertical lines indicate the epochs of spectroscopy and the corresponding optical spectroscopic classifications reported for WOH G64.

Figure 1 Long description
Panel: A line graph showing the optical light curves of WOH G64 over a period of more than 30 years. The horizontal axis represents the Modified Julian Date (MJD) ranging from 1993 to 2025. The vertical axis represents the magnitude (mag) ranging from 12 to 20. Multiple datasets are labeled with corresponding surveys and filters: MACHO R_KC and V_KC, OGLE I and V, RP, G, ATLAS O, BP, and C. The shadowed patch between 2013.5 and 2014.5 indicates the epoch of the transition. Vertical lines mark the epochs of spectroscopy and corresponding optical spectroscopic classifications reported for WOH G64.
An extraordinary change in the optical spectral features accompanied the dramatic photometric transition in 2014 (Fig. 2). In 2007, WOH G64 exhibits spectral features characteristic of an M6 Ie-type RSG (Solf Reference Solf1978), including deep TiO bands and strong Ca <sc>ii</sc> triplet absorption lines. Remarkably, we detect strong double-peaked emission from forbidden emission lines (e.g., [N <sc>ii</sc>] and [S <sc>ii</sc>]), strangely seen in RSGs. Multi-epoch spectroscopy after the transition revealed the most extreme change ever seen in the optical spectral features of a RSG. In 2016 and 2021, the optical spectral range is dominated by asymmetric emission lines from permitted (e.g., Fe <sc>i</sc>, H <sc>i</sc>, and He <sc>i</sc>) and forbidden transitions (e.g., [Ca <sc>ii</sc>], [Fe <sc>ii</sc>], [N <sc>ii</sc>], and [S <sc>ii</sc>]) leading to a B[e]-star classification (Lamers et al. Reference Lamers, Zickgraf, de Winter, Houziaux and J.1998; Kraus Reference Kraus2019). Simultaneous near-IR spectroscopy in 2016 revealed a composite spectrum, with the H band characterized by neutral atomic absorption lines typically seen in late-G to early-K stars (Rayner et al. Reference Rayner, Cushing and Vacca2009). Fig. 3 presents a comparison of the H-band spectrum of WOH G64 with that of a supergiant B[e] star (sgB[e]), two Yellow Supergiants (YSGs), and a RSG. If WOH G64 was a sgB[e], its full spectral range would consistently resemble that of a hot B-type star embedded in low-density circumstellar material. However, the H-band features observed in WOH G64 closely resemble those of YSGs, suggesting an effective temperature of
$T_{\mathrm{eff}} \sim 4700$ K (Fig. 3). The absence of typical RSG features in the H-band further supports the spectral transition reported in the optical. The composite spectral features of a warm star in the near-IR and a B[e] star in the optical indicate that WOH G64 is a massive symbiotic binary (Lamers et al. Reference Lamers, Zickgraf, de Winter, Houziaux and J.1998), where the main RSG has evolved into a YHG.
Comparison of the multi-epoch optical spectroscopy of WOH G64. The main spectral features are indicated for each epoch. The spectra are scaled to unity between 8700-8800 Å, but the UVES spectrum is shown with an offset.

Comparison of the H band of WOH G64 from X-Shooter (blue) in 2016 with the APOGEE spectra of four supergiants of the LMC, indicating their T
$_{\mathrm{eff}}$ from APOGEE.

3. Discussion
Binary interactions in a symbiotic system provide a compelling framework to explain the observations of WOH G64, including its asymmetric morphology and surrounding torus (Ohnaka et al. Reference Ohnaka, Driebe, Hofmann, Weigelt and Wittkowski2008; Ohnaka et al. Reference Ohnaka, Hofmann and Weigelt2024), the He <sc>I</sc> emission detected after 2016, likely produced by a hot companion (e.g., Cowley Reference Cowley1969) or an accretion disk around it (e.g., Beristain et al. Reference Beristain, Edwards and Kwan2001), the double-peaked forbidden emission tracing bipolar outflows from the accretor during mass transfer (Soker Reference Soker2015), and the radial velocity shifts seen in several lines.
Binary interactions could therefore have driven the observed transition. Mass transfer is expected in RSGs within wide binaries, potentially leading to a CE phase (Paczyński and Sienkiewicz Reference Paczyński and Sienkiewicz1972; Ercolino et al. Reference Ercolino, Jin, Langer and Dessart2024). However, predicting the outcome or duration of CE evolution remains challenging, as it depends sensitively on stellar parameters and involves complex, multi-timescale 3D processes that are difficult to model (Ivanova et al. Reference Ivanova, Justham and Chen2013). WOH G64 has long exhibited an unusually cool and extended atmosphere compared to typical RSGs (Levesque Reference Levesque2017), which may have been a pseudo-atmosphere formed during a CE phase. During the CE evolution, the system could enter a slow, self-regulated phase where orbital decay energy is transported outward and radiated efficiently (Clayton et al. Reference Clayton, Podsiadlowski, Ivanova and Justham2017). This process may drive large-amplitude, recurrent pulsations that eject shells every few decades, yielding an average mass-loss rate of
$\dot{M}\sim10^{-3}$ M
$_\unicode{x2299}$ yr
$^{-1}$, comparable to the estimations for WOH G64 (
$\dot{M}\gt10^{-4}$ M
$_\unicode{x2299}$ yr
$^{-1}$; Antoniadis et al. Reference Antoniadis, Bonanos and de Wit2024).
Beyond binary interactions, the transition could be driven by single-stellar physics. The luminosity of WOH G64 is close to the Humphreys-Davidson limit where significant eruptions are expected (Glatzel and Kiriakidis Reference Glatzel and Kiriakidis1993). For instance, Var A, one of the most luminous stars in M33 (Hubble and Sandage Reference Hubble and Sandage1953), underwent an eruption lasting over 45 yr and serves as a notable example (Humphreys et al. Reference Humphreys, Jones and Gehrz1987; Humphreys et al. Reference Humphreys, Jones and Polomski2006). During the eruption, Var A exhibited late-M spectral features due to its optically thick wind. When the eruption ended, Var A returned to its YHG state, exhibiting F-type features along with both forbidden and permitted emission lines. Similar to Var A, WOH G64 might have been a YHG undergoing an eruption for more than 30 yr, until it returned to a quiescent state in 2014 when the false pseudo-atmosphere dissipated. These two cases are outstanding because other YHGs exhibit typically short-lived eruptions (e.g., Koumpia et al. Reference Koumpia, Oudmaijer and Graham2020). However, the dramatic optical changes in WOH G64, fast transition, and higher obscuration suggest an even more extreme event than in Var A.
Our findings challenge current models of massive star evolution by underscoring the critical role of binary interactions and episodic eruptions in shaping the fate of evolved massive stars. In particular, WOH G64 offers valuable insight into the late evolutionary stages of the most luminous RSGs, and the long-standing “RSG problem”. Recent studies of the extreme RSGs VY CMa and NML Cyg have uncovered potential evidence of companions (Smith et al. Reference Smith, Beasor and Close2025; De Beck et al. Reference De Beck, Andrews, Quintana-Lacaci and Vlemmings2025). Determining whether binarity is responsible for their extreme properties and their vast episodic mass-loss events is crucial for understanding the final evolutionary path of massive stars (Marchant and Bodensteiner Reference Marchant and Bodensteiner2024). In the case of WOH G64, its future remains uncertain due to limited constraints on the system’s physical and orbital parameters. Continued spectroscopic and photometric monitoring of this unique system will be essential to constrain its binary properties, uncover the mechanism behind its dramatic transition, and predict its ultimate fate.



>

eff