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.

To decode the information stored within a spectrum, detailed modelling of the physical state is required together with accurate radiative transfer solution schemes. In the analysis of stellar spectra, the numerical model often needs to account for high velocity outflows, multi-dimensional structures, and the effects of binary companions. Focusing now on binary systems, we present the BOSS-3D spectral synthesis code, which is capable of calculating synthetic line profiles for a variety of binary systems. Assuming the state of the circumstellar material to be known, the standard pz-geometry is extended by defining individual coordinate systems for each object. By embedding these coordinate systems within the observer’s frame, BOSS-3D automatically accounts for outflows or discs within both involved systems, and includes all Doppler shifts. Moreover, the code accounts for different length-scales, and thus could also be used to analyse transit-spectra of planetary atmospheres. As a first application of BOSS-3D, we model the phase-dependent line profiles for the enigmatic binary (or multiple) system LB-1.

A significant fraction of the stars near the tip of the AGB phase become regular or semi-regular (Mira-type, SRs) pulsators. However, some of these light curves have shown intriguing secondary minima or sharp dips with much longer periods. Although this phenomenon shows some resemblance with the R CrB variables, the light curve is generally symmetric before and after the dip, whereas in R CrB the luminosity recovers slower after its minimum. More recently, high-resolution ALMA CO observations revealed a spiral structure around some of these stars, which suggests the presence of a stellar or sub-stellar companion. In these cases, the long-term light curve minima could be caused by periodic eclipses of the primary by a spiral circumstellar structure, and the long-period would be related to the orbital period. In this paper we discuss the pros and cons of the various proposed scenarios for the long-term minima of pulsating AGB stars.

It is well established that late-type main-sequence (MS) stars display a relationship between X-ray activity and the Rossby number, Ro, the ratio of rotation period to the convective turnover time. This manifests itself as a saturated regime (where X-ray activity is constant) and an unsaturated regime (where X-ray activity anti-correlates with Rossby number). However, this relationship breaks down for the fastest rotators. We cross-correlated SuperWASP visually classified photometric light curves and All-Sky Automated Survey for Supernovae automatically classified photometric light curves with XMM-Newton X-ray observations to identify 3 178 stars displaying a photometrically defined rotational modulation in their light curve and corresponding X-ray observations. We fitted a power-law to characterise the rotation–activity relation of 900 MS stars. We identified that automatically classified rotationally modulated light curves are not as reliable as visually classified light curves for this work. We found a power-law index in the unsaturated regime of G- to M-type stars of $\beta=-1.84\,\pm\,0.18$ for the SuperWASP catalogue, in line with the canonical value of $\beta=-2$ . We find evidence of supersaturation in the fastest rotating K-type stars, with a power-law index of $\beta_{s}=1.42\pm0.26$ .

The growing number of black hole binary (BHB) mergers detected by the Laser Interferometer Gravitational-Wave Observatory have the potential to enable an unprecedented characterisation of the physical processes and astrophysical conditions that govern the formation of compact binaries. In this paper, we focus on investigating the dynamical formation of BHBs in dense star clusters through a state-of-the-art set of 58 direct N-body simulations with N $\leqslant200\,000$ particles which include stellar evolution, gravitational braking, orbital decay through gravitational radiation, and galactic tidal interactions. The simulations encompass a range of initial conditions representing typical young globular clusters, including the presence of primordial binaries. The systems are simulated for $\sim 12$ Gyr. The dataset yields 117 BHB gravitational wave (GW) events, with 97 binaries merging within their host cluster and 20 merging after having been ejected. Only 8% of all ejected BHBs merge within the age of the Universe. Systems in this merging subset tend to have smaller separations and larger eccentricities, as this combination of parameters results in greater emission of gravitational radiation. We confirm known trends from Monte Carlo simulations, such as the anti-correlation between the mass of the binary and age of the cluster. In addition, we highlight for the first time a difference at low values of the mass ratio distribution between in-cluster and ejected mergers. However, the results depend on assumptions on the strength of GW recoils, thus in-cluster mergers cannot be ruled out at a significant level of confidence. A more substantial catalogue of BHB mergers and a more extensive library of N-body simulations are needed to constrain the origin of the observed events.