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.

We highlight a few results from ALMA Band 6 observations of high-rotational transition number (J) SiO masers towards oxygen-rich AGB and red supergiant stars carried out as part of the ATOMIUM Large Programme in 2018–2020. A search for a relationship between mass-loss rates and flux-weighted mean angular distances of maser components was inconclusive, as linear regression models for the 28SiO v=1 J =5–4 and J =6–5 transitions were inconsistent. Supplementary APEX observations towards the ATOMIUM AGB stars also suggest variability at different stellar pulsation phases.

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.

This contribution presents recent advances in identifying the stellar upper mass limit using simulations of UV radiative feedback during the star formation process. Generally, due to computational costs and a focus on au to parsec scales, simulations of massive star formation do not trace the flow of material to distances closer than a few au from the forming star. However, UV line-acceleration acts directly on accreting material in the sub-au circumstellar region, thereby efficiently ablating the surface layers off the protostellar disk. For stars on the order of a few hundred solar masses, this disk destruction rate exceeds the accretion rate, destroying the disk faster than it is replenished, and setting a maximum stellar mass as a function of metallicity that can be attained by single star formation channels.

Debris disks are the natural by-products of the star and planet formation processes. Since the 1980’s several thousands of debris disks have been detected, and the presence of a disk is inferred by the detection of excess emission over the photospheric emission. This thermal emission arises from (micron-sized or slightly bigger) dust grains heated by the central star. However, in the vast majority of cases, these observations are not spatially resolving the radial distribution of the dust, resulting in strong degeneracies in the modeling approach (radial distance vs minimum grain size mostly). Therefore the properties of the small dust grains remained largely unconstrained until the arrival of high angular resolution instruments, especially at optical and near-infrared wavelengths. In these proceeding some of the main results are presented that have been obtained over the past few years on the properties of small dust grains in debris disks, and it is discussed how laboratory experiments contributed to put those results in context.

UV radiation plays a critical role in the chemistry of circumstellar envelopes (CSEs) around evolved stars on the asymptotic giant branch (AGB). However, the effects of all potential sources of UV radiation have not been included in models. We present preliminary results of adding an internal source of UV to the CSE chemistry and predict large enhancements of atomic and ionic species arising from photo-destruction of parent species. Observations of atomic carbon towards the UV-bright AGB star o Cet are consistent with the modelling. In addition, we calculate the precise depth dependence of the CO photodissociation rate in an expanding CSE. We incorporate this within a chemical network active in the outflows of AGB stars, which includes 933 species and 15182 reactions. Our results show that the CO envelope size is about 30% smaller at half abundance than the most commonly used radius reported by Mamon et al. (1988).

We derive azimuthally-averaged surface-brightness profiles of 16 AGB stars in the far-IR and sub-mm with the aim of studying the resolved historic mass loss in the extended circumstellar envelope. The PSF-subtracted extended component fluxes were found to be ∼40% of the total source flux. By fitting SEDs at each radial point we derive the dust temperature, column density and spectral index of emissivity via Bayesian inference. The measured dust-to-gas ratios were somewhat consistent with canonical values however with a large scatter.

A long-standing problem of the general paradigm of low-mass star formation is the “luminosity problem”: protostars are less luminous than theoretically predicted. One possible solution is that the accretion process is episodic. FU Orionis-type stars (FUors) are thought to be the visible examples for objects in the high accretion state and it is still debated what physical mechanism triggers the phenomenon. For many of these objects their disk properties are still largely unknown so we conducted a deep, high spatial resolution (down to 20 au) ALMA Band 6 (1.3 mm) dust continuum survey of a sub-sample of known FUors. Here we present preliminary results of our survey, including the mass, size and spectral slope of each disk.

The Korean VLBI Network (KVN) is a unique millimeter VLBI system which is consisted of three 21 m telescopes with relatively short baselines. We present the preliminary results of simultaneous monitoring observations of the 22.2 GHz H2O and 43.1/42.8/86.2/129.3 GHz SiO masers based on the KVN Key Science Project (KSP). We obtained the astrometrically registered maps of the H2O and SiO masers toward nine evolved stars using the source frequency phase referencing method (SFPR). The SFPR maps of the H2O and SiO masers enabled us to investigate the spatial structure and kinematics from the SiO to H2O maser regions including the development of an outward motion from the ring-like or elliptical structures of SiO masers to the asymmetric structures of the 22.2 GHz H2O maser features. In particular, the 86.2/129.3 GHz SiO (v=1, J=2–1 and J=3–2) masers were clearly imaged toward several objects for the first time. The SiO v=1, J=3–2 maser shows different distributions compared to those of the SiO v=1, 2, J=1–0 and v=1, J=2–1 masers implying a different physical condition.

The radio non-detection of two Type Ia supernovae (SNe) SN 2011fe and SN 2014J has been modeled considering synchrotron radiation from shock accelerated electrons in the SN shock fronts. With 10% each of the bulk kinetic energy in electric and magnetic fields, a very low density of the medium around both the SNe has been estimated from the null detection of radio emission, around 1 and 4 years after the explosion of SNe 2014J and 2011fe, respectively. Keeping the fraction of energy in electrons fixed at 10%, a medium with particle density ~ 1cm−3 is found when 1% of the post shock energy is in magnetic fields. In case of a wind medium, the former predicts the mass loss rate Ṁ to be <10−9M⊙ yr−1, and the latter gives an upper limit ~10−9M⊙ yr−1, for wind velocity of 100 kms−1, for both the SNe. The tenuous media obtained from this study favor the double degenerate as well as a spin up/down model for both SNe 2011fe and 2014J.

Eta Carinae is one of the most massive observable binaries. Yet determination of its orbital and physical parameters is hampered by obscuring winds. However the effects of the strong, colliding winds changes with phase due to the high orbital eccentricity. We wanted to improve measures of the orbital parameters and to determine the mechanisms that produce the relatively brief, phase-locked minimum as detected throughout the electromagnetic spectrum. We conducted intense monitoring of the He ii λ4686 line in η Carinae for 10 months in the year 2014, gathering ~300 high S/N spectra with ground- and space-based telescopes. We also used published spectra at the FOS4 SE polar region of the Homunculus, which views the minimum from a different direction. We used a model in which the He ii λ4686 emission is produced by two mechanisms: a) one linked to the intensity of the wind-wind collision which occurs along the whole orbit and is proportional to the inverse square of the separation between the companion stars; and b) the other produced by the ‘bore hole’ effect which occurs at phases across the periastron passage. The opacity (computed from 3D SPH simulations) as convolved with the emission reproduces the behavior of equivalent widths both for direct and reflected light. Our main results are: a) a demonstration that the He ii λ4686 light curve is exquisitely repeatable from cycle to cycle, contrary to previous claims for large changes; b) an accurate determination of the longitude of periastron, indicating that the secondary star is ‘behind’ the primary at periastron, a dispute extended over the past decade; c) a determination of the time of periastron passage, at ~4 days after the onset of the deep light curve minimum; and d) show that the minimum is simultaneous for observers at different lines of sight, indicating that it is not caused by an eclipse of the secondary star, but rather by the immersion of the wind-wind collision interior to the inner wind of the primary.

The evolution of massive stars encompasses short-lived transition phases in which mass-loss is more enhanced and usually eruptive. A complex environment, combining atomic, molecular and dust regions, is formed around these stars. In particular, the circumstellar environment of B[e] Supergiants is not well understood. To address that, we have initiated a campaign to investigate their environments for a sample of Galactic and Magellanic Cloud sources. Using high-resolution optical and near-infrared spectra (MPG-ESO/FEROS, GEMINI/Phoenix and VLT/CRIRES, respectively), we examine a set of emission features ([OI], [CaII], CO bandheads) to trace the physical conditions and kinematics in their formation regions. We find that the B[e] Supergiants are surrounded by a series of rings of different temperatures and densities, a probable result of previous mass-loss events. In many cases the CO forms very close to the star, while we notice also an alternate mixing of densities and temperatures (which give rise to the different emission features) along the equatorial plane.

This article presents an overview of the published results for planetary nebulae based on images and spectroscopy from the PACS, SPIRE, and HIFI instruments on board the Herschel satellite.

I look briefly at protostellar environments and at what information can be gained from the millimetre waveband. A short overview of the theory of protostellar environments and evidence of disks will be followed by a more detailed review of recent publications of observations of protostellar halos/disks, specifically in 3 mm.

Using data from the WISE All-Sky Survey, we have found >100 new infrared excess sources around main-sequence Hipparcos stars within 75 pc. Our empirical calibration of WISE photospheric colors and removal of non-trivial false-positive sources are responsible for the high confidence (>99.5%) of detections, while our corrections to saturated W1 and W2 photometry have for the first time allowed us to search for new infrared excess sources around bright field stars in WISE. The careful calibration and filtering of the WISE data have allowed us to probe excess fluxes down to roughly 8% of the photospheric emission at 22μm around saturated stars in WISE. We expect that the increased sensitivity of our survey will not only aid in understanding the evolution of debris disks, but will also benefit future studies using WISE.

Debris disks are usually thought to be gas-poor, the gas being dissipated by accretion or evaporation during the protoplanetary phase. HD141569A is a 5 Myr old star harboring a famous debris disk, with multiple rings and spiral features. I present here the first PdBI maps of the 12CO(2−1), 13CO(2−1) gas and dust emission at 1.3 mm in this disk. The analysis reveals there is still a large amount of (primordial) gas extending out to 250 AU, i.e. inside the rings observed in scattered light. HD141569A is thus a hybrid disk with a huge debris component, where dust has evolved and is produced by collisions, with a large remnant reservoir of gas.

Using 22 μm data from the Wide Field Infrared Survey Explorer (WISE), we have completed a sensitive all-sky survey for debris disks in Hipparcos and Tycho catalog stars within 120 pc. This warm excess emission traces material in the circumstellar region likely to host terrestrial planets. Several hundred previously unknown debris disk candidates were identified. We are currently performing follow-up observations to characterize the stars, companions, and circumstellar material in these systems with a variety of facilities including Keck, Herschel, and HST. Thirteen WISE debris disks have been observed to date using HST/STIS coronagraphy. Five of these disks have been detected in scattered light. One is a large and highly asymmetric edge-on disk which appears to be both warped and bifurcated.

We compare in a systematic way spectrometric, photometric and mid-infrared (VLTI/MIDI) interferometric measurements with different types of model atmospheres. Self-consistent dynamic model atmospheres in particular were used to interpret in a consistent way the dynamic behavior of gas and dust. The results underline how the joint use of different kind of observations, as photometry, spectroscopy and interferometry, is essential to understand the atmospheres of pulsating C-rich AGB stars. The sample of C-rich stars discussed in this work provides crucial constraints for the atmospheric structure.

A rare subclass of Type Ia supernovae (SNe Ia) shows evidence of strong interaction with a hydrogen-rich circumstellar medium (CSM); these objects are referred to as SNe Ia-CSM. PTF11kx began life as a SN Ia, but after a month it began to show indications of significant interaction with its CSM. This well-studied object solidified the connection between SNe Ia-CSM and more typical SNe Ia, despite their spectral similarity to Type IIn SNe (which likely come from massive star progenitors, as opposed to the white dwarf progenitors for the SNe Ia-CSM). The spectra of all ~20 known SNe Ia-CSM are dominated by Hα emission (with widths of ~2000 km s−1) and exhibit large Hα/Hβ intensity ratios; moreover, they have an almost complete lack of He I emission (see left panel of Figure 1). They also show possible evidence of dust formation through a decrease in the red wing of Hα 75–100 days past maximum brightness. The absolute magnitudes of SNe Ia-CSM are found to be -21.3 mag ≤ MR ≤ −19 mag (see right panel of Figure 1), and they also show ultraviolet emission at early times and strong infrared emission at late times (but no detected radio or X-ray emission). Finally, the host galaxies of SNe Ia-CSM imply that these objects come from a relatively young stellar population.

Starlight polarization provides insight into the physical mechanisms in and around the source as well as its geometry, whether or not the source is resolved. In this talk we will review mechanisms that polarize light in stellar envelopes. The observations and modeling can be used to probe the physics of the circumstellar environment as well as its relation to the ambient interstellar environment.