The unidentified infrared (UIR) bands, whose carriers are thought to be organics, have been widely observed in various astrophysical environments. However, our knowledge of the detailed chemical composition and formation process of the carriers is still limited. We have synthesized laboratory organics named Quenched Nitrogen-included Carbonaceous Composite (QNCC) by quenching plasma produced from nitrogen gas and hydrocarbon solids. Infrared and X-ray analyses of QNCC showed that infrared properties of QNCC well reproduce the UIR bands observed in novae and amine structures contained in QNCC play an important role in the origin of the broad 8 m feature, which characterizes the UIR bands in novae. QNCC is at present the best laboratory analog of organic dust formed around dusty classical novae, which carries the UIR bands in novae via thermal emission process [Endo et al.(2021)].

Although red supergiants (RSGs) are observed to be undergoing vigorous mass loss, explaining the mechanism launching their winds has been a long-standing problem. Given the importance of mass loss to stellar evolution in this phase, this is a key uncertainty. In this contribution we present a recently published model (Kee et al. 2021) showing that turbulent pressure alone can extend the stellar atmosphere of an RSG to the degree that a wind is launched. This provides a fully analytic mass-loss prescription for RSGs. Moreover, utilising observationally inferred turbulent velocities for these objects, we find that this wind can carry an appropriate amount of mass to overall match observations. Intriguingly, when coupled to stellar evolution models the predicted mass-loss rates show that stars with initial masses above Mini∼17M⊙ may naturally evolve back to the blue and as such not end their lives as RSGs; this is also in overall good agreement with observations, here of Type II-P/L supernova progenitors. Moreover, since the proposed wind launching mechanism is not necessarily sensitive to metallicity, this could have important implications for stellar evolution predictions in low-metallicity environments.

A rich zoo of peculiar objects forms when Asymptotic Giant Branch (AGB) stars, undergo interactions in a binary system. For example, Barium (Ba) stars are main-sequence and red-giant stars that accreted mass from the outflows of a former AGB companion, which is now a dim white dwarf (WD). Their orbital properties can help us constrain AGB binary interaction mechanisms and their chemical abundances are a tracer of the nucleosynthesis processes that took place inside the former AGB star. The observational constraints concerning the orbital and stellar properties of Ba stars have increased in the past years, but important uncertainties remained concerning their WD companions. In this contribution, we used HD 76225 to demonstrate that by combining radial-velocity data with Hipparcos and Gaia astrometry, one can accurately constrain the orbital inclinations of these systems and obtain the absolute masses of these WDs, getting direct information about their AGB progenitors via initial-final mass relationships.

The link between hot and cool stellar outflows is shown to be critical for correctly predicting the masses of the most massive black holes (BHs) below the so-called pair-instability supernova (PISN) mass gap. Gravitational Wave (GW) event 190521 allegedly hosted an “impossibly” heavy BH of 85 M ⊙. Here we show how our increased knowledge of both metallicity Z and temperature dependent mass loss is critical for our evolutionary scenario of a low-Z blue supergiant (BSG) progenitor of an initially approx 100 M ⊙ star to work. We show using MESA stellar evolution modelling experiments that as long as we can keep such stars above 8000 K such low-Z BSGs can avoid strong winds, and keep a very large envelope mass intact before core collapse. This naturally leads to the Cosmic Time dependent maximum BH function below the PISN gap.

Massive stars are amongst the rarest but also most intriguing stars. Their extreme, magnetised stellar winds induce, by wind-ISM interaction, famous multi-wavelengths circumstellar gas nebulae of various morphologies, spanning from large-scale wind bubbles to stellar wind bow shocks, rings and bipolar shapes. We present two- and three-dimensional magneto-hydrodynamical (MHD) simulations of the circumstellar medium of such massive stars at different phase of their evolution. Particularly, we investigate the stability properties of 3D MHD bow shock nebulae around the runaway red supergiant stars IRC-10414 and Betelgeuse. Our results show that their astrospheres are stabilised by an organised, non-parallel ambient magnetic field. These findings suggest that Betelgeuse’s bar is of interstellar origin. Last, we explore the circular aspect of the young nebula around the Wolf-Rayet stars. It is found that Wolf-Rayet nebulae are not affected by the ISM gas distribution in which the stellar objects lie, even in the case of fast stellar motion: as testifies the ring-like surroundings of the Milky Way’s fastest Wolf-Rayet star, WR124. The morphology of these nebulae is tightly related to their pre-Wolf-Rayet wind geometry and to their phase evolution transition properties, which can generate bipolar shapes. We will further discuss their diffuse projected emission by means of radiative transfer calculations and show that the projected diffuse emission can appear as bipolar structures as in NGC6888.

The winds observed around asymptotic giant branch (AGB) stars are generally attributed to radiation pressure on dust, which is formed in the extended dynamical atmospheres of these pulsating, strongly convective stars. Current radiation-hydrodynamical models can explain many of the observed features, and they are on the brink of delivering a predictive theory of mass loss. This review summarizes recent results and ongoing work on winds of AGB stars, discussing critical ingredients of the driving mechanism, and first results of global 3D RHD star-and-wind-in-a-box simulations. With such models it becomes possible to follow the flow of matter, in full 3D geometry, all the way from the turbulent, pulsating interior of an AGB star, through its atmosphere and dust formation zone into the region where the wind is accelerated by radiation pressure on dust. Advanced instruments, which can resolve the stellar atmospheres, where the winds originate, provide essential data for testing the models.

The mechanisms driving mass loss from massive stars in late stages of their evolution is still very much unknown. Stellar evolution models indicate that the last stage before going supernova for many massive stars is the Wolf-Rayet (WR) phase, characterized by a strong, optically thick stellar wind. Stellar models show that these stars exceed the Eddington limit already in deep sub-surface layers around the so-called ‘iron-opacity’ bump, and so should launch a supersonic outflow from there. However, if the outward force does not suffice to accelerate the gas above the local escape speed, the initiated flow will stagnate and start raining down upon the stellar core. In previous, spherically symmetric, WR wind models, this has been circumvented by artificially increasing either clumping or the line force. Here, we present pioneering 3D time-dependent radiation-hydrodynamic simulations of WR winds. In these models, computed without any ad-hoc force enhancement, the stagnated flow leads to co-existing regions of up- and down-flows, which dynamically interact with each other to form a multi-dimensional and complex outflow. These density structures, and the resulting highly non-monotonic velocity field, can have important consequences for mass-loss rates and the interpretation of observed Wolf-Rayet spectra.

The post-main sequence evolutionary path of massive stars comprises various transition phases, in which the stars shed large amounts of material into their environments. Our studies focus on two of them: B[e] supergiants and yellow hypergiants, for which we investigate the structure and dynamics within their environments. We find that each B[e] supergiant is surrounded by a unique set of rings or arc-like structures. These structures are either stable over time or they display high variability, including expansion and dilution. In contrast, yellow hypergiants are embedded in multiple shells of gas and dust. These objects are famous for their outburst activity. Moreover, the dynamics in their extended atmospheres imply an enhanced pulsation activity prior to outburst. The physical mechanism(s) leading to episodic mass ejections in these two types of stars is still uncertain. We propose that strange-mode instabilities, excited in the inflated envelopes of these objects, play a significant role.

Classical Wolf-Rayet (WR) stars mark an important stage in the late evolution of massive stars. As hydrogen-poor massive stars, these objects have lost their outer layers, while still losing further mass through strong winds indicated by their prominent emission line spectra. Wolf-Rayet stars have been detected in a variety of different galaxies. Their strong winds are a major ingredient of stellar evolution and population synthesis models. Yet, a coherent theoretical picture of their strong mass-loss is only starting to emerge. In particular, the occurrence of WR stars as a function of metallicity (Z) is still far from being understood.

To uncover the nature of the complex and dense winds of Wolf-Rayet stars, we employ a new generation of model atmospheres including a consistent solution of the wind hydrodynamics in an expanding non-LTE situation. With this technique, we can dissect the ingredients driving the wind and predict the resulting mass-loss for hydrogen-depleted massive stars. Our modelling efforts reveal a complex picture with strong, non-linear dependencies on the luminosity-to-mass ratio and Z with a steep, but not totally abrupt onset for WR-type winds in helium stars. With our findings, we provide a theoretical motivation for a population of helium stars at low Z, which cannot be detected via WR-type spectral features. Our study of massive He-star atmosphere models yields the very first mass-loss recipe derived from first principles in this regime. Implementing our first findings in stellar evolution models, we demonstrate how traditional approaches tend to overpredict WR-type mass loss in the young Universe.

Interacting binaries within a common envelope, wherein the primary is a red giant are believed to result in a recently identified evolutionary class – the dusty post-RGB stars. Our SED modeling of eight post-RGBs in the LMC indicates the presence of geometrically thick disks with substantial opening angle in addition to the outer shells. We estimated the total dust mass (and gas mass assuming gas-to-dust ratio) in the disks and shells and set constraints on the dust grain compositions and sizes. The only known Galactic object of this class is the Boomerang nebula. Additionally, we present a DUSTY model of the Boomerang that can serve as a template for 2D modeling of the object using RADMC-3D. 2D modeling is essential to dissect the morphology of the spatially-unresolved post-RGBs in the LMC. These models may then be tested with future HST and ALMA imaging, together with JWST spectroscopy of these objects.

To look at propagating winds from evolved stars into the interstellar medium is to look at how they are sustained. To understand their origins, we must look to the circumstances that create them in the first instance. In this article, I examine the physical conditions under which pulsation-enhanced, dust-driven winds are first generated. These initial conditions can help constrain the late evolutionary stages of these stars and provide insight into the mechanisms that cause the mass loss itself.

Low- and intermediate mass stars experience a significant mass loss during the last phases of their evolution, which obscures them in a vast, dusty envelope. Although it has long been thought this envelope is generally spherically symmetric in shape, recent high-resolution observations find that most of these stars exhibit complex and asymmetrical morphologies, most likely resulting from binary interaction. In order to improve our understanding about these systems, theoretical studies are needed in the form of numerical simulations. Currently, a handful of simulations exist, albeit they mainly focus on the hydrodynamics of the outflow. Hence, we here present the pathway to more detailed and accurate modelling of companion-perturbed outflows with, by discussing the missing but crucial physical and chemical processes. With these state-of-the-art simulations we aim to make a direct comparison with observations to unveil the true identity on the embedded systems.

A recent analysis of a few carbon-oxygen white dwarfs in old open clusters of the Milky Way (MW) identified a kink in the initial-final mass relation (IFMR), located over a range of initial masses, 1.65 ≲ Mi/M⊙ ≲ 2.10, which unexpectedly interrupts the commonly assumed monotonic trend. The proposed interpretation links this observational fact to the formation of carbon stars and the modest outflows (with mass loss rate < 10−7 M⊙/yr) that are expected as long as the carbon excess remains too low to produce dust grains in sufficient amount. Under these conditions the mass of the carbon-oxygen core can grow more than is generally predicted by stellar models. We discuss these new findings also in light of a new systematic follow-up investigation, based on Gaia EDR3, of evolved giants (13 carbon stars, 3 S stars and 4 M stars) belonging to intermediate-age open clusters.

Some binary stars experience common envelope evolution, which is accompanied by drastic loss of angular momentum, mass, and orbital energy and which leaves behind close binaries often involving at least one white dwarf, neutron star, or black hole. The best studied phase of common envelope is the dynamical inspiral lasting few original orbital periods. We show theoretical interpretation of observations of V1309 Sco and AT2018bwo revealing that binaries undergo substantial prolonged mass loss before the dynamical event amounting up to few solar masses. This mass loss is concentrated in the orbital plane in the form of an outflow or a circumbinary disk. Collision between this slower mass loss and the subsequent faster dynamical ejection powers a bright red transient. The resulting radiative shock helps to shape the explosion remnant and provides a site of dust and molecule formation.

Understanding the nucleosynthetic origin of nitrogen and the evolution of the N/O ratio in the interstellar medium is crucial for a comprehensive picture of galaxy chemical evolution at high-redshift because most observational metallicity (O/H) estimates are implicitly dependent on the N/O ratio. The observed N/O at high-redshift shows an overall constancy with O/H, albeit with a large scatter. We show that these heretofore unexplained features can be explained by the pre-supernova wind yields from rotating massive stars (M≳10M⊙,ν/νcrit≳0.4). Our models naturally produce the observed N/O plateau, as well as the scatter at low O/H. We find the scatter to arise from varying star formation efficiency. However, the models that have supernovae dominated yields produce a poor fit to the observed N/O at low O/H. This peculiar abundance pattern at low O/H suggests that dwarf galaxies are most likely to be devoid of SNe yields and are primarily enriched by pre-supernova wind abundances.

The cosmic origin of fluorine is still under debate. Asymptotic giant branch (AGB) stars are among the few suggested candidates to efficiently synthesis F in our Galaxy, however their relative contribution is not clear. In this paper, we briefly review the theoretical studies from stellar yield models of the F synthesis and chemical equilibrium models of the F-containing molecules in the outflow around AGB stars. Previous detections of the F-bearing species towards AGB and post-AGB stars are also highlighted. We suggest that high-resolution ALMA observations of the AlF, one of the two main carriers of F in the outflow of AGB stars, can provide a reliable tracer of the F-budget in AGB stars. This will be helpful to quantify the role of AGB stars in the Galactic F budget.

Evolved massive stars are major cosmic engines, providing strong mechanical and radiative feedback on their host environment. They contribute to the enrichment of their environment through a strong stellar winds, still poorly understood. Wind physics across the life cycle of these stars is the key ingredient to accomplish a complete understanding of their evolution in the near and distant Universe. Nowadays, the development of the observational instruments is so advanced that the observations became very sensitive to the details of the stellar surface making possible to quantitatively study what happens on their surfaces and above where the stellar winds become dominant. Three-dimensional radiative hydrodynamics simulations of evolved stars are essential to a proper and quantitative analysis of these observations. This work presents how these simulations have been (and will be) crucial to prepare and interpret a multitude of observations and how they are important to achieve the knowledge of the mass-loss mechanism.

Most of a star’s mass is bound in a hydrostatic equilibrium in which pressure balances gravity. But if at some near-surface layer additional outward forces overcome gravity, this can transition to a supersonic, outflowing wind, with the sonic point, where the outward force cancels gravity, marking the division between hydrostatic atmosphere and wind outflow. This talk will review general issues with such transonic initiation of a stellar wind outflow, and how this helps set the wind mass loss rate. The main discussion contrasts the flow initiation in four prominent classes of steady-state winds: (1) the pressure-driven coronal wind of the sun and other cool stars; (2) line-driven winds from OB stars; (3) a two-stage initiation model for the much denser winds from Wolf-Rayet (WR) stars; and (4) the slow “overflow” mass loss from highly evolved giant stars. A follow on discussion briefly reviews eruptive mass loss, with particular focus on the giant eruption of η Carinae.

We present interferometric continuum and molecular line emission maps obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) of OH231.8+4.2, a well studied bipolar nebula around an asymptotic giant branch (AGB) star that is key to understand the remarkable changes in nebular morphology and kinematics during the short transition from the AGB to the Planetary Nebula (PN) phase. The excellent angular resolution of our maps (∼20 mas ∼30 AU) allows us to scrutinize the central nebular regions of OH231.8+4.2, which hold the clues to understanding how this iconic object assembled its complex nebular architecture. We report, for the first time in this object and others of its kind (i.e. pre-PNe with massive bipolar outflows), the discovery of a rotating circumbinary disk of radius ∼30 AU traced by NaCl, KCl, and H2O emission lines. The disk lies at the base of a young bipolar wind with signs of rotation as well. A compact spatially resolved dust disk is found perpendicular to the bipolar outflow. We also identify a point-like continuum source, which likely represents the central Mira star enshrouded by a ∼3 R* shell or disk of hot (∼1400 K) freshly formed dust. The point source is slightly off-centre from the disk centroid, enabling us for the first time to place constraints to the orbital separation of the central binary system.

This contribution focuses on a rare example of the class of post-Red Supergiants, IRAS 17163-3907, the central star of the Fried Egg nebula. In particular, we discuss some of our recently published results in detail. The inner parts of the circumstellar environment of this evolved massive star are probed at milli-arcsec resolution using VLTI’s GRAVITY instrument operating in the K-band (2 µm), while larger, arcsecond, scales are probed by VISIR diffraction limited images around 10 µm, supplemented by a complete Spectral Energy Distribution. The spectro-interferometric data cover important diagnostic lines (Brγ, Na I), which we are able to constrain spatially. Both the presence and size of the Na i doublet in emission has been traditionally challenging to explain towards other objects of this class. In this study we show that a two-zone model in Local Thermal Equilibrium can reproduce both the observed sizes and strengths of the emission lines observed in the K-band, without the need of a pseudo-photosphere. In addition, we find evidence for the presence of a third hot inner shell, and demonstrate that the star has undergone at least three mass-loss episodes over roughly the past century. To explain the properties of the observed non-steady mass-loss we explore pulsation-driven and line-driven mass-loss and introduce the bi-stability jump as a possible underlying mechanism to explain mass-loss towards Yellow Hypergiants.