We present the results of KVN Key Science Program (KSP) for evolved stars, which was launched in 2014. The first phase of KSP ended in June 2020 and the second phase started in October 2020. The goal of KSP is to study the physical characteristics of the evolved stars by observing the spatial distribution and temporal variability of the stellar masers at four frequency-bands (K, Q, W and D bands). The 22 GHz H2O maser is usually observed from the outer part of circumstellar envelopes compared to the 43, 86, 129 GHz SiO masers, thus the kinematic links between these regions can be studied by the multi-frequency simultaneous observations of KSP along the stellar pulsation cycles. This eventually enable us to study the enormous mass-loss rate of evolved stars, and the accumulated results from KSP are expected to shed light on the study of the late stage of the stellar evolution.

Eight planetary nebulae (PNe) are known to emit OH and/or H2O masers, but there is no report of an SiO maser in this type of objects. We present a search for SiO masers in 16 confirmed and candidate PNe, carried out with the Australia Telescope Compact Array. We found no evidence of association between SiO masers and PNe in our data. Previous detections of thermal SiO emission in PNe show that these molecules can be present in gas phase in this type of objects, but it is not yet clear whether they can be found where the physical conditions are appropriate for maser pumping. We suggest that the best candidates for a new search would be PNe showing high-velocity outflows.

Intense mass loss through cool, low-velocity winds is a defining characteristic of low-to-intermediate mass stars during the asymptotic giant branch (AGB) evolutionary stage. Such winds return up ∼80% of the initial stellar mass to the interstellar medium and play a major role in enriching it with dust and heavy elements. A challenge to understanding the physics underlying AGB mass loss is its dependence on an interplay between complex and highly dynamic processes, including pulsations, convective flows, shocks, magnetic fields, and opacity changes resulting from dust and molecule formation. I highlight some examples of recent advances in our understanding of late-stage stellar mass loss that are emerging from radio and (sub)millimeter observations, with a particular focus on those that resolve the surfaces and extended atmospheres of evolved stars in space, time, and frequency.

Maser effects occur in recombination lines when the plasma departs from local thermodynamic equilibrium (LTE). Its consequence is not as dramatic as that found in molecular masers, and therefore it is more difficult to recognize. Besides, it occurs in compact high density regions, and its lines fall at millimeter and submillimeter wavelengths, only recently available with good angular resolution. This review will focus on the theoretical aspects of maser recombination lines and on the recent detection of these masers in different astronomical objects.

We have investigated the spectral evolutions of H2O and SiO masers associated with 12 “water fountain” sources in our FLASHING (Finest Legacy Acquisitions of SiO-/H2O-maser Ignitions by Nobeyama Generation) project. Our monitoring observations have been conducted using the Nobeyama 45 m telescope every 2 weeks–2 months since 2018 December except during summer seasons. We have found new extremely high velocity H2O maser components, breaking the records of jet speeds in this type of sources. Systematic line-of-sight velocity drifts of the H2O maser spectral peaks have also been found, indicating acceleration of the entrained material hosting the masers around the jet. Moreover, by comparing with previous spectral data, we can find decadal growths/decays of H2O maser emission. Possible periodic variations of the maser spectra are further being inspected in order to explore the periodicity of the central stellar system (a pulsating star or a binary). Thus we expect to see the real-time evolution/devolutions of the water fountains over decades.

We report new detections of SiO ν = 1 and ν = 2 J = 1 → 0 masers in the “water fountain” source IRAS 16552-3050, which was observed with the Nobeyama 45 m telescope from March 2021 to April 2023. Water fountains are evolved stars whose H2O maser spectra trace high-velocity outflows of >100 Km s−. This is the second known case of SiO masers in a water fountain, after their prototypical source, W 43A. The line-of-sight velocity of the SiO masers are blue-shifted by ∼25 km s−1 from the systemic velocity. This velocity offset imply that the SiO masers are associated with nozzle structure formed by a jet penetrating the circumstellar envelope, and that new gas blobs of the jet erupted recently. Thus, the SiO masers imply this star to be in a new evolutionary stage.

In this contribution, I briefly review the long-term evolution of the solar wind (its mass-loss rate), including the evolution of observed properties that are intimately linked to the solar wind (rotation, magnetism and activity). I also briefly discuss implications of the evolution of the solar wind on the evolving Earth. I argue that studying exoplanetary systems could open up new avenues for progress to be made in our understanding of the evolution of the solar wind.

Wolf-Rayet stars are regarded as candidates for progenitors of core-collapse supernovae, and they are expected to be progenitors of long gamma-ray bursts. These types of stars are considered to be fast rotators. Their high rotation speed breaks the sphericity of the star and leads to an axisymmetric wind density structure. In such a case, the electron scattering takes place in a nonspherical environment, and as a result, we might expect an intrinsic polarization. We present a 2.5D radiation hydrodynamic stellar wind model of these stars. The model simulations account for the deformation of the stellar surface due to rotation, gravity darkening, and nonradial forces. We computed the polarization from the density variable of the hydrodynamic model, derived the upper limit of rotational velocities, and found no conflict with the previous studies of Wolf-Rayet stars.

During the evolution of massive stars, their properties change significantly. But stellar parameters of massive stars have the biggest uncertainties in stellar astrophysics, specifically in the post-main sequence stages where blue supergiant stars are located. These stars experience mass loss events during their evolution which are supposed to be related to strange-modes instabilities. In this work, we explore the stability of oscillation modes in massive stars for different masses, and a range of mass-loss rates, with the aim to provide clues about the connection between strange modes instabilities and mass-loss events.

The role of mass loss from massive stars, especially episodic mass loss, is one of the outstanding open questions facing stellar evolution theory. Multiple lines of evidence are pointing to violent, episodic mass-loss events being responsible for removing a large part of the massive stellar envelope, especially in low-metallicity galaxies. The ERC ASSESS project aims to determine whether episodic mass loss is a dominant process in the evolution of the most massive stars by conducting the first extensive, multi-wavelength survey of evolved massive stars in the nearby Universe. The project hinges on the fact that mass-losing stars form dust and are bright in the mid-infrared. We aim to investigate the properties of evolved targets in nearby galaxies and estimate the amount of ejected mass, which will constrain evolutionary models. In this work we present some of our first observational results from the galaxies NGC 6822 and IC 10 obtained with OSIRIS (GTC).

UV wind line variability in OB stars appears to be universal. We review the evidence that the variability is due to large, dense, optically thick structures rooted in or near the photosphere. Using repeated observations and a simple model we translate observed profile variations into optical depth variations and, consequently, variations in measured mass loss rates. Although global rates may be stable, measured rates vary. Consequently, profile variations infer how mass loss rates determined from UV wind lines vary. These variations quantify the intrinsic error inherent in any mass loss rate derived from a single observation. These derived rates can differ by factors of 3 or more. Our results also imply that rates from non-simultaneous observations (such as UV and ground based data) need not agree. Finally, we use our results to examine the nature of the structures responsible for the variability.

Hot, massive stars are known to host unstable, radiation-driven outflowing winds, giving rise to dense clumps of material which severely affect the diagnostic techniques used to derive wind properties of massive stars. Most of the current diagnostic models account for wind inhomogeneities by assuming a one-component medium consisting of optically thin clumps, and maintaining a smooth velocity-field. However, this neglects important light-leakage effects through porous channels in-between the clumps. These light-leakage effects have recently been incorporated in the stellar atmosphere modelling code FASTWIND, and here we will present quantitative mass-loss results from a combined Ultraviolet-Optical wind analysis of O-supergiants in the Galaxy. Using a genetic-algorithm fitting-approach, we systematically investigate the impact the wind physics has on derived stellar and wind parameters, and how this depends on metallicity and spectral type. We compare our findings with earlier results (which do not take into account such light-leakage effects), to standard mass-loss rates usually included in evolution model studies of massive stars, and with theoretical predictions of clumping properties. We will also present the first systematic empirical constraints on the new wind parameters, associated with light-leakage, and compare these with theoretical predictions.

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.

Mass loss through stellar winds plays a dominant role in the evolution of massive stars. Very massive stars (VMSs, > 100Mȯ) display Wolf-Rayet spectral morphologies (WNh) whilst on the main-sequence. Bestenlehner (2020) extended the elegant and widely used stellar wind theory by Castor, Abbott & Klein (1975) from the optically thin (O star) to the optically thick main-sequence (WNh) wind regime. The new mass-loss description is able to explain the empirical mass-loss dependence on the Eddington parameter and is suitable for incorporation into stellar evolution models for massive and very massive stars. The prescription can be calibrated with the transition mass-loss rate defined in Vink & Gräfener (2012). Based on the stellar sample presented in Bestenlehner et al. (2014) we derive a mass-loss recipe for the Large Magellanic Cloud using the new theoretical mass-loss prescription of Bestenlehner (2020).

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.

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.

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 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.

We discovered a nebula in the low-metallicity (Z=0.1 ) nearby ( Mpc) dwarf galaxy NGC 4068, which reveals broad Hα line profile and unusual emission line fluxes in its spectrum. The object also shows significant nitrogen overabundance, not typical for metal-poor environment. We assumed that the nebula could be ionized by an evolved massive star with , Wolf-Rayet or Blue Supergiant, and built models of the nebula using Cloudy and CMFGEN codes. Our models successfully reproduce the optical emission spectrum of the object, including the peculiar [SII]/[NII] ratio and the presence of the HeII λ4686 line.

Spectroscopy can decode the radiation from stars in an appropriate way and derive many properties of different stellar objects. In this work we seek to derive simultaneously stellar and wind parameters of massive stars. To model the data we use the radiative transport code Fastwind with the hydrodynamic solutions derived using our stationary code Hydwind as input, instead of the β-law. Then, ISOSCELES, our grid of stellar atmosphere and hydrodynamic models of massive stars, is used to derive the physical properties of the observed spectra through spectral line fittings. This quantitative spectroscopic analysis provide an estimation about the line–force parameters, whose theoretical calculations are complex. In addition, we expect to confirm that the hydrodynamic δ-slow solutions, describe quite reliable the radiation line-driven winds of A and late B supergiant stars and, at the same time, explain disagreements between observational data and theoretical models for the Wind–Momentum Luminosity Relationship (WLR).