Context. Red supergiant stars possess surface features and extended molecular atmospheres. Photospheric convection may be a crucial factor of the levitation of the outer atmospheric layers. However, the mechanism responsible is still poorly understood. Aims. We image the stellar surface of V602 Carinae (V602 Car) to constrain the morphology and contrast of the surface features and of the extended atmospheric layers. Methods. We observed V602 Car with the Very Large Telescope Interferometer (VLTI) PIONIER instrument (1.53-1.78 $\mathrm{\mu}$m) between May and July 2016, and April and July 2019 with different telescope configurations. We compared the image reconstructions with 81 temporal snapshots of 3D radiative-hydrodynamics (RHD) CO$^5$BOLD simulations in terms of contrast and morphology, using the Structural Similarity Index. Results. The interferometric data are compatible with an overall spherical disk of angular diameter 4.4$\pm$0.2 mas, and an extended molecular layer. In 2016, the reconstructed image reveals a bright arc-like feature toward the northern rim of the photospheric surface. In 2019, an arc-like feature is seen at a different orientation and a new peak of emission is detected on the opposite side. The contrasts of the reconstructed surface images are 11%$\pm$2% and 9%$\pm$2% for 2016 and 2019, respectively. The morphology and contrast of the two images are consistent with 3D RHD simulations, within our achieved spatial resolution and dynamic range. The extended molecular layer contributes 10--13% of the total flux with an angular diameter of 6--8 mas. It is present but not clearly visible in the reconstructed images because it is close to the limits of the achieved dynamic range. The presence of the molecular layer is not reproduced by the 3D RHD simulations. Conclusions...

The AGB star $\pi^{1}$ Gruis has a known companion (at a separation of ~400 AU) which cannot explain the strong deviations from the spherical symmetry of the CSE. Recently, hydrodynamic simulations of mass transfer in closer binary systems have successfully reproduced the spiral-shaped CSEs found around a handful of sources. There is growing evidence for an even closer, undetected companion complicating the case of $\pi^{1}$ Gruis further. The improved spatial resolution allows for the investigation of the complex circumstellar morphology and the search for imprints on the CSE of the third component. We have observed the 12CO J=3-2 line emission from $\pi^{1}$ Gruis using both the compact and extended array of Atacama Large Millimeter/submillimeter Array (ALMA). The interferometric data has furthermore been combined with data from the ALMA total power (TP) array. The imaged brightness distribution has been used to constrain a non-local, non-LTE 3D radiative transfer model of the CSE. The high-angular resolution ALMA data have revealed the first example of a source on the AGB where both a faster bipolar outflow and a spiral pattern along the orbital plane can be seen in the gas envelope. The spiral can be traced in the low- to intermediate velocity, from 13 to 25 km s$^{-1}$, equatorial torus. The largest spiral-arm separation is $\approx$5".5 and consistent with a companion with an orbital period of $\approx$330 yrs and a separation of less than 70 AU. The kinematics of the bipolar outflow is consistent with it being created during a mass-loss eruption where the mass-loss rate from the system increased by at least a factor of 5 during 10-15 yrs. The spiral pattern is the result of an undetected companion. The bipolar outflow is the result of a rather recent mass-loss eruption event.

Cool evolved stars are the main source of chemical enrichment of the interstellar medium. Understanding their mass loss offers a unique opportunity to study the cycle of matter. We discuss interferometric studies and their comparison to latest state-of-the-art dynamic model atmospheres. They show broad agreement for asymptotic giant branch stars. For red supergiants, however, current models cannot explain observed extensions by far, pointing to missing physical processes in their models, and uncertainties in our general understanding of mass loss. We present ongoing imaging and time-series observations that may provide the strongest constraint and may help to identify missing dynamic processes. VLTI studies will remain the highest spatial resolution observations at ESO into the ELT era, complemented by ALMA observations. We discuss crucial improvements in both instrumental and operational areas for surface imaging of cool evolved stars in the era of the ELT.

The purpose of this work is to present an extensive grid of dynamical atmosphere and wind models for M-type AGB stars, covering a wide range of relevant stellar parameters. We used the DARWIN code, which includes frequency-dependent radiation-hydrodynamics and a time-dependent description of dust condensation and evaporation, to simulate the dynamical atmosphere. The wind-driving mechanism is photon scattering on submicron-sized Mg$_2$SiO$_4$ grains. The grid consists of $\sim4000$ models, with luminosities from $L_\star=890\,{\mathrm{L}}_\odot$ to $L_\star=40000\,{\mathrm{L}}_\odot$ and effective temperatures from 2200K to 3400K. For the first time different current stellar masses are explored with M-type DARWIN models, ranging from 0.75M$_\odot$ to 3M$_\odot$. The modelling results are radial atmospheric structures, dynamical properties such as mass-loss rates and wind velocities, and dust properties (e.g. grain sizes, dust-to-gas ratios, and degree of condensed Si). We find that the mass-loss rates of the models correlate strongly with luminosity. They also correlate with the ratio $L_*/M_*$: increasing $L_*/M_*$ by an order of magnitude increases the mass-loss rates by about three orders of magnitude, which may naturally create a superwind regime in evolution models. There is, however, no discernible trend of mass-loss rate with effective temperature, in contrast to what is found for C-type AGB stars. We also find that the mass-loss rates level off at luminosities higher than $\sim14000\,{\mathrm{L}}_\odot$, and consequently at pulsation periods longer than $\sim800$ days. The final grain radii range from 0.25 micron to 0.6 micron. The amount of condensed Si is typically between 10% and 40%, with gas-to-dust mass ratios between 500 and 4000.

Models of the chemical evolution of the interstellar medium, galaxies, and the Universe rely on our understanding of the amounts and chemical composition of the material returned by stars and supernovae. Stellar yields are obtained from stellar-evolution models, which currently lack predictive prescriptions of stellar mass loss, although it significantly affects stellar lifetimes, nucleosynthesis, and chemical ejecta. Galaxy properties are derived from observations of the integrated light of bright member stars. Stars in the late stages of their evolution are among the infrared-brightest objects in galaxies. An unrealistic treatment of the mass-loss process introduces significant uncertainties in galaxy properties derived from their integrated light. We describe current efforts and future needs and opportunities to characterize AGB outflows: driving mechanisms, outflow rates, underlying fundamental physical and chemical processes such as dust grain formation, and dependency of these on metallicity.

The thermally-pulsing asymptotic giant branch (TP-AGB) experienced by low- and intermediate-mass stars is one of the most uncertain phases of stellar evolution and the models need to be calibrated with the aid of observations. To this purpose, we couple high-quality observations of resolved stars in the Small Magellanic Cloud (SMC) with detailed stellar population synthesis simulations computed with the TRILEGAL code. The strength of our approach relies on the detailed spatially-resolved star formation history of the SMC, derived from the deep near-infrared photometry of the VISTA survey of the Magellanic Clouds, as well as on the capability to quickly and accurately explore a wide variety of parameters and effects with the COLIBRI code for the TP-AGB evolution. Adopting a well-characterized set of observations -- star counts and luminosity functions -- we set up a calibration cycle along which we iteratively change a few key parameters of the TP-AGB models until we eventually reach a good fit to the observations. Our work leads to identify two best-fitting models that mainly differ in the efficiencies of the third dredge-up and mass loss in TP-AGB stars with initial masses larger than about 3 M$_{\odot}$. On the basis of these calibrated models we provide a full characterization of the TP-AGB stellar population in the SMC in terms of stellar parameters (initial masses, C/O ratios, carbon excess, mass-loss rates). Extensive tables of isochrones including these improved models are publicly available.

Dust grains forming in the extended atmospheres of AGB stars are critical for the heavy mass loss of these cool luminous giants, as they provide radiative acceleration for the stellar winds. Characteristic mid-IR spectral features indicate that the grains consist mainly of silicates and corundum. The latter species seems to form in a narrow zone within about 2 stellar radii, preceding the condensation of silicate dust, which triggers the outflow. Recent high-angular-resolution observations show clumpy, variable dust clouds at these distances. We explore possible causes for the formation of inhomogeneous dust layers, using 3D dynamical simulations. We modeled the outer convective envelope and the dust-forming atmosphere of an M-type AGB star with the CO5BOLD radiation-hydrodynamics code. The simulations account for frequency-dependent gas opacities, and include a time-dependent description of grain growth and evaporation for corundum and olivine-type silicates. In the inner, gravitationally bound, and corundum-dominated layers of the circumstellar envelope, a patchy distribution of the dust emerges naturally, due to atmospheric shock waves that are generated by large-scale convective flows and pulsations. The formation of silicate dust at somewhat larger distances probably indicates the outer limit of the gravitationally bound layers. The current models do not describe wind acceleration, but the cloud formation mechanism should also work for stars with outflows. Timescales of atmospheric dynamics and grain growth are similar to observed values. In spherical averages of dust densities the variable 3D morphology manifests itself as cycle-to-cycle variations. Grain growth in the wake of large-scale non-spherical shock waves, generated by convection and pulsations, is a likely mechanism for producing the observed clumpy dust clouds, and for explaining their physical and dynamical properties.

The wind-driving mechanism of asymptotic giant branch (AGB) stars is commonly attributed to a two-step process: first, gas in the stellar atmosphere is levitated by shockwaves caused by stellar pulsation, then accelerated outwards by radiative pressure on newly formed dust, inducing a wind. Dynamical modelling of such winds usually assumes a spherically symmetric star. We explore the potential consequences of complex stellar surface structures, as predicted by three-dimensional (3D) star-in-a-box modelling of M-type AGB stars, on the resulting wind properties with the aim to improve the current wind models. Two different modelling approaches are used; the CO$^5$BOLD 3D star-in-a-box code to simulate the convective, pulsating interior and lower atmosphere of the star, and the DARWIN one-dimensional (1D) code to describe the dynamical atmosphere where the wind is accelerated. The gas dynamics of the inner atmosphere region at distances of $R\sim1-2R_\star$, which both modelling approaches simulate, are compared. Dynamical properties and luminosity variations derived from CO$^5$BOLD interior models are used as input for the inner boundary in DARWIN wind models in order to emulate the effects of giant convection cells and pulsation, and explore their influence on the dynamical properties. The CO$^5$BOLD models are inherently anisotropic, with non-uniform shock fronts and varying luminosity amplitudes, in contrast to the spherically symmetrical DARWIN wind models. DARWIN wind models with CO$^5$BOLD-derived inner boundary conditions produced wind velocities and mass-loss rates comparable to the standard DARWIN models, however the winds show large density variations on time-scales of 10-20 years.

The S-type asymptotic giant branch (AGB) star $\pi^{1}$ Gruis has a known companion at a separation of $\approx$400 AU. The envelope structure, including an equatorial torus and a fast bipolar outflow, is rarely seen in the AGB phase and is particularly unexpected in such a wide binary system. Therefore a second, closer companion has been suggested, but the evidence is not conclusive. The new ALMA $^{12}$CO and $^{13}$CO $J$=3-2 data, together with previously published $^{12}$CO $J$=2-1 data from the Submillimeter Array (SMA), and the $^{12}$CO $J$=5-4 and $J$=9-8 lines observed with Herschel/Heterodyne Instrument for the Far-Infrared (HIFI), is modeled with the 3D non-LTE radiative transfer code SHAPEMOL. The data analysis clearly confirms the torus-bipolar structure. The 3D model of the CSE that satisfactorily reproduces the data consists of three kinematic components: a radially expanding torus with velocity slowly increasing from 8 to 13 km s$^{-1}$ along the equator plane; a radially expanding component at the center with a constant velocity of 14 km s$^{-1}$; and a fast, bipolar outflow with velocity proportionally increasing from 14 km s$^{-1}$ at the base up to 100 km s$^{-1}$ at the tip, following a linear radial dependence. The results are used to estimate an average mass-loss rate during the creation of the torus of 7.7$\times$10$^{-7}$ M$_{\odot}$ yr$^{-1}$. The total mass and linear momentum of the fast outflow are estimated at 7.3$\times$10$^{-4}$ M$_{\odot}$ and 9.6$\times$10$^{37}$ g cm s$^{-1}$, respectively. The momentum of the outflow is in excess (by a factor of about 20) of what could be generated by radiation pressure alone, in agreement with recent findings for more evolved sources. The best-fit model also suggests a $^{12}$CO/$^{13}$CO abundance ratio of 50. Possible shaping scenarios for the gas envelope are discussed

The heavy mass loss observed in evolved asymptotic giant branch stars is usually attributed to a two-step process: atmospheric levitation by pulsation-induced shock waves, followed by radiative acceleration of newly formed dust grains. Detailed wind models suggest that the outflows of M-type AGB stars may be triggered by photon scattering on Fe-free silicates with grain sizes of about 0.1 - 1 $\mu$m. Due to the low grain temperature, these Fe-free silicates can condense close to the star, but they do not produce the characteristic mid-IR features that are often observed in M-type AGB stars. However, it is probable that the silicate grains are gradually enriched with Fe as they move away from the star, to a degree where the grain temperature stays below the sublimation temperature, but is high enough to produce emission features. We investigate whether differences in grain temperature in the inner wind region, which are related to changes in the grain composition, can be detected with current interferometric techniques, in order to put constraints on the wind mechanism. To investigate this we use radial structures of the atmosphere and wind of an M-type AGB star, produced with the 1D radiation-hydrodynamical code DARWIN. The spectral energy distribution is found to be a poor indicator of different temperature profiles and therefore is not a good tool for distinguishing different scenarios of changing grain composition. However, spatially resolved interferometric observations have promising potential. They show signatures even for Fe-free silicates (found at 2-3 stellar radii), in contrast to the spectral energy distribution. Observations with baselines that probe spatial scales of about 4 stellar radii and beyond are suitable for tracing changes in grain composition, since this is where effects of Fe enrichment should be found.

Wind-driving in asymptotic giant branch (AGB) stars is commonly attributed to a two-step process. First, matter in the stellar atmosphere is levitated by shock waves, induced by stellar pulsation, and second, this matter is accelerated by radiation pressure on dust, resulting in a wind. In dynamical atmosphere and wind models the effects of the stellar pulsation are often simulated by a simplistic prescription at the inner boundary. We test a sample of dynamical models for M-type AGB stars, for which we kept the stellar parameters fixed to values characteristic of a typical Mira variable but varied the inner boundary condition. The aim was to evaluate the effect on the resulting atmosphere structure and wind properties. The results of the models are compared to observed mass-loss rates and wind velocities, photometry, and radial velocity curves, and to results from 1D radial pulsation models. Dynamical atmosphere models are calculated, using the DARWIN code for different combinations of photospheric velocities and luminosity variations. The inner boundary is changed by introducing an offset between maximum expansion of the stellar surface and the luminosity and/or by using an asymmetric shape for the luminosity variation. Models that resulted in realistic wind velocities and mass-loss rates, when compared to observations, also produced realistic photometric variations. For the models to also reproduce the characteristic radial velocity curve present in Mira stars (derived from CO $\Delta v = 3$ lines), an overall phase shift of 0.2 between the maxima of the luminosity and radial variation had to be introduced. We find that a group of models with different boundary conditions (29 models, including the model with standard boundary conditions) results in realistic velocities and mass-loss rates, and in photometric variations.

Winds of asymptotic giant branch (AGB) stars are commonly assumed to be driven by radiative acceleration of dust grains. For M-type AGB stars, the nature of the wind-driving dust species has been a matter of intense debate. A proposed source of the radiation pressure triggering the outflows is photon scattering on Fe-free silicate grains. This wind-driving mechanism requires grain radii of about 0.1 - 1 micron in order to make the dust particles efficient at scattering radiation around the stellar flux maximum. Grain size is therefore an important parameter for understanding the physics behind the winds of M-type AGB stars. We seek to investigate the diagnostic potential of scattered polarized light for determining dust grain sizes. We have developed a new tool for computing synthetic images of scattered light in dust and gas shells around AGB stars, which can be applied to detailed models of dynamical atmospheres and dust-driven winds. We present maps of polarized light using dynamical models computed with the DARWIN code. The synthetic images clearly show that the intensity of the polarized light, the position of the inner edge of the dust shell, and the size of the dust grains near the inner edge are all changing with the luminosity phase. Non-spherical structures in the dust shells can also have an impact on the polarized light. We simulate this effect by combining different pulsation phases into a single 3D structure before computing synthetic images. An asymmetry of the circumstellar envelope can create a net polarization, which can be used as diagnostics for the grain size. The ratio between the size of the scattering particles and the observed wavelength determines at what wavelengths net polarization switches direction. If observed, this can be used to constrain average particle sizes.

Context: Observations of asymptotic giant branch (AGB) stars with increasing spatial resolution reveal new layers of complexity of atmospheric processes on a variety of scales. Aim: To analyze the physical mechanisms that cause asymmetries and surface structures in observed images, we use detailed 3D dynamical simulations of AGB stars; these simulations self-consistently describe convection and pulsations. Methods: We used the CO5BOLD radiation-hydrodynamics code to produce an exploratory grid of global "star-in-a-box" models of the outer convective envelope and the inner atmosphere of AGB stars to study convection, pulsations, and shock waves and their dependence on stellar and numerical parameters. Results: The model dynamics are governed by the interaction of long-lasting giant convection cells, short-lived surface granules, and strong, radial, fundamental-mode pulsations. Radial pulsations and shorter wavelength, traveling, acoustic waves induce shocks on various scales in the atmosphere. Convection, waves, and shocks all contribute to the dynamical pressure and, thus, to an increase of the stellar radius and to a levitation of material into layers where dust can form. Consequently, the resulting relation of pulsation period and stellar radius is shifted toward larger radii compared to that of non-linear 1D models. The dependence of pulsation period on luminosity agrees well with observed relations. The interaction of the pulsation mode with the non-stationary convective flow causes occasional amplitude changes and phase shifts. The regularity of the pulsations decreases with decreasing gravity as the relative size of convection cells increases. The model stars do not have a well-defined surface. Instead, the light is emitted from a very extended inhomogeneous atmosphere with a complex dynamic pattern of high-contrast features.

We present near-infrared interferometry of the carbon-rich asymptotic giant branch (AGB) star R Sculptoris. The visibility data indicate a broadly circular resolved stellar disk with a complex substructure. The observed AMBER squared visibility values show drops at the positions of CO and CN bands, indicating that these lines form in extended layers above the photosphere. The AMBER visibility values are best fit by a model without a wind. The PIONIER data are consistent with the same model. We obtain a Rosseland angular diameter of 8.9+-0.3 mas, corresponding to a Rosseland radius of 355+-55 Rsun, an effective temperature of 2640+-80 K, and a luminosity of log L/Lsun=3.74+-0.18. These parameters match evolutionary tracks of initial mass 1.5+-0.5 Msun and current mass 1.3+-0.7 Msun. The reconstructed PIONIER images exhibit a complex structure within the stellar disk including a dominant bright spot located at the western part of the stellar disk. The spot has an H-band peak intensity of 40% to 60% above the average intensity of the limb-darkening-corrected stellar disk. The contrast between the minimum and maximum intensity on the stellar disk is about 1:2.5. Our observations are broadly consistent with predictions by dynamic atmosphere and wind models, although models with wind appear to have a circumstellar envelope that is too extended compared to our observations. The detected complex structure within the stellar disk is most likely caused by giant convection cells, resulting in large-scale shock fronts, and their effects on clumpy molecule and dust formation seen against the photosphere at distances of 2-3 stellar radii.

The recent progress in high-spatial-resolution techniques, spanning wavelengths from the visual to the radio regime, is leading to new valuable insights into the complex dynamical atmospheres of Asymptotic Giant Branch (AGB) stars and their wind forming regions. Striking examples are images of asymmetries and inhomogeneities in the photospheric and dust-forming layers which vary on time-scales of months. These features are probably related to large-scale convective flows predicted by 3D star-in-a-box models. Furthermore, high-resolution observations make it possible to measure dust condensation distances, and they give information about the chemical composition and sizes of dust grains in the close vicinity of cool giants. These are essential constraints for building realistic models of wind acceleration and developing a predictive theory of mass loss for AGB stars, which is a crucial ingredient of stellar and galactic chemical evolution models.

This paper presents a summary of four invited and twelve contributed presentations on asymptotic giant branch stars and red supergiants, given over the course of two afternoon splinter sessions at the 19th Cool Stars Workshop. It highlights both recent observations and recent theory, with some emphasis on high spatial resolution, over a wide range of wavelengths. Topics covered include 3D models, convection, binary interactions, mass loss, dust formation and magnetic fields.

High spatial resolution techniques have given valuable insights into the mass loss mechanism of AGB stars, which presumably involves a combination of atmospheric levitation by pulsation-induced shock waves and radiation pressure on dust. Observations indicate that Al$_2$O$_3$ condenses at distances of about 2 stellar radii or less, prior to the formation of silicates. Al$_2$O$_3$ grains are therefore prime candidates for producing the scattered light observed in the close vicinity of several M-type AGB stars, and they may be seed particles for the condensation of silicates at lower temperatures. We have constructed a new generation of Dynamic Atmosphere & Radiation-driven Wind models based on Implicit Numerics (DARWIN), including a time-dependent treatment of grain growth & evaporation for both Al$_2$O$_3$ and Fe-free silicates (Mg$_2$SiO$_4$). The equations describing these dust species are solved in the framework of a frequency-dependent radiation-hydrodynamical model for the atmosphere & wind structure, taking pulsation-induced shock waves and periodic luminosity variations into account. Condensation of Al$_2$O$_3$ at the close distances and in the high concentrations implied by observations requires high transparency of the grains in the visual and near-IR region to avoid destruction by radiative heating. For solar abundances, radiation pressure due to Al$_2$O$_3$ is too low to drive a wind. Nevertheless, this dust species may have indirect effects on mass loss. The formation of composite grains with an Al$_2$O$_3$ core and a silicate mantle can give grain growth a head start, increasing both mass loss rates and wind velocities. Furthermore, our experimental core-mantle grain models lead to variations of visual and near-IR colors during a pulsation cycle which are in excellent agreement with observations.

Winds of AGB stars are thought to be driven by a combination of pulsation-induced shock waves and radiation pressure on dust. In dynamic atmosphere and wind models, the stellar pulsation is often simulated by prescribing a simple sinusoidal variation in velocity and luminosity at the inner boundary of the model atmosphere. We experiment with different forms of the luminosity variation in order to assess the effects on the wind velocity and mass-loss rate, when progressing from the simple sinusoidal recipe towards more realistic descriptions. Using state-of-the-art dynamical models of C-rich AGB stars, a range of different asymmetric shapes of the luminosity variation and a range of phase shifts of the luminosity variation relative to the radial variation are tested. These tests are performed on two stellar atmosphere models. The first model has dust condensation and, as a consequence, a stellar wind is triggered, while the second model lacks both dust and wind. The first model with dust and stellar wind is very sensitive to moderate changes in the luminosity variation. There is a complex relationship between the luminosity minimum, and dust condensation: changing the phase corresponding to minimum luminosity can either increase or decrease mass-loss rate and wind velocity. The luminosity maximum dominates the radiative pressure on the dust, which in turn, is important for driving the wind. These effects of changed luminosity variation are coupled with the dust formation. In contrast there is very little change to the structure of the model without dust. Changing the luminosity variation, both by introducing a phase shift and by modifying the shape, influences wind velocity and the mass-loss rate. To improve wind models it would probably be desirable to extract boundary conditions from 3D dynamical interior models or stellar pulsation models.

Dust jets, i.e. fuzzy collimated streams of cometary material arising from the nucleus, have been observed in-situ on all comets since the Giotto mission flew by comet 1P/Halley in 1986. Yet their formation mechanism remains unknown. Several solutions have been proposed, from localized physical mechanisms on the surface/sub-surface (see review in Belton (2010)) to purely dynamical processes involving the focusing of gas flows by the local topography (Crifo et al. 2002). While the latter seems to be responsible for the larger features, high resolution imagery has shown that broad streams are composed of many smaller features (a few meters wide) that connect directly to the nucleus surface. We monitored these jets at high resolution and over several months to understand what are the physical processes driving their formation, and how this affects the surface. Using many images of the same areas with different viewing angles, we performed a 3-dimensional reconstruction of collimated jets, and linked them precisely to their sources on the nucleus. Results.We show here observational evidence that the Northern hemisphere jets of comet 67P arise from areas with sharp topographic changes and describe the physical processes involved. We propose a model in which active cliffs are the main source of jet-like features, and therefore the regions eroding the fastest on comets. We suggest that this is a common mechanism taking place on all comets.

The winds of cool luminous AGB stars are commonly assumed to be driven by radiative acceleration of dust grains which form in the extended atmospheres produced by pulsation-induced shock waves. The dust particles gain momentum by absorption or scattering of stellar photons, and they drag along the surrounding gas particles through collisions, triggering an outflow. This scenario, here referred to as Pulsation-Enhanced Dust-DRiven Outflow (PEDDRO), has passed a range of critical observational tests as models have developed from empirical and qualitative to increasingly self-consistent and quantitative. A reliable theory of mass loss is an essential piece in the bigger picture of stellar and galactic chemical evolution, and central for determining the contribution of AGB stars to the dust budget of galaxies. In this review, I discuss the current understanding of wind acceleration and indicate areas where further efforts by theorists and observers are needed.

Stellar winds observed in asymptotic giant branch (AGB) stars are usually attributed to a combination of stellar pulsations and radiation pressure on dust. Shock waves triggered by pulsations propagate through the atmosphere, compressing the gas and lifting it to cooler regions, which create favourable conditions for grain growth. If sufficient radiative acceleration is exerted on the newly formed grains through absorption or scattering of stellar photons, an outflow can be triggered. Strong candidates for wind-driving dust species in M-type AGB stars are magnesium silicates (Mg$_2$SiO$_4$ and MgSiO$_3$). Such grains can form close to the stellar surface, they consist of abundant materials and, if they grow to sizes comparable to the wavelength of the stellar flux maximum, they experience strong acceleration by photon scattering. We use a frequency-dependent radiation-hydrodynamics code with a detailed description for the growth of Mg$_2$SiO$_4$ grains to calculate the first extensive set of time-dependent wind models for M-type AGB stars. The resulting wind properties, visual and near-IR photometry and mid-IR spectra are compared with observations.We show that the models can produce outflows for a wide range of stellar parameters. We also demonstrate that they reproduce observed mass-loss rates and wind velocities, as well as visual and near-IR photometry. However, the current models do not show the characteristic silicate features at 10 and 18 $\mu$m as a result of the cool temperature of Mg$_2$SiO$_4$ grains in the wind. Including a small amount of Fe in the grains further out in the circumstellar envelope will increase the grain temperature and result in pronounced silicate features, without significantly affecting the photometry in the visual and near-IR wavelength regions.

Asymptotic giant branch (AGB) stars play a key role in the enrichment of galaxies with heavy elements. Due to their large amplitude variability, the measurement of elemental abundances is a highly challenging task that has not been solved in a satisfactory way yet. Following our previous work we use hydrostatic and dynamical model atmospheres to simulate observed high-resolution near-infrared spectra of 12 variable and non-variable red giants in the globular cluster 47 Tuc. The 47 Tuc red giants are independently well-characterized in important parameters (mass, metallicity, luminosity). The principal aim was to compare synthetic spectra based on the dynamical models with observational spectra of 47 Tuc variables. Assuming that the abundances are unchanged on the upper giant branch in these low-mass stars, our goal is to estimate the impact of atmospheric dynamics on the abundance determination. We present new measurements of the C/O and 12C/13C ratio for 5 non-variable red giants in 47Tuc. The equivalent widths measured for our 7 variable stars strongly differ from the non-variable stars and cannot be reproduced by either hydrostatic or dynamical model atmospheres. Nevertheless, the dynamical models fit the observed spectra of long-period variables much better than any hydrostatic model. For some spectral features, the variations in the line intensities predicted by dynamical models over a pulsation cycle give similar values as a sequence of hydrostatic models with varying temperature and constant surface gravity.

The evolution and spectral properties of stars on the AGB are significantly affected by mass loss through dusty stellar winds. Dynamic atmosphere and wind models are an essential tool for studying these evolved stars, both individually and as members of stellar populations, to understand their contribution to the integrated light and chemical evolution of galaxies. This paper is part of a series testing state-of-the-art atmosphere and wind models of carbon stars against observations, and making them available for use in various theoretical and observational studies. We have computed low-resolution spectra and photometry (in the wavelength range 0.35-25 mu) for a grid of 540 dynamic models with stellar parameters typical of solar-metallicity C-rich AGB stars and with a range of pulsation amplitudes. The models cover the dynamic atmosphere and dusty outflow (if present), assuming spherical symmetry, and taking opacities of gas-phase species and dust grains consistently into account. To characterize the time-dependent dynamic and photometric behaviour of the models in a concise way we defined a number of classes for models with and without winds. Comparisons with observed data in general show a quite good agreement for example regarding mass-loss rates vs. (J-K) colours or K magnitudes vs. (J-K) colours. Some exceptions from the good overall agreement, however, are found and attributed to the range of input parameters (e.g. relatively high carbon excesses) or intrinsic model assumptions (e.g. small particle limit for grain opacities). While current results indicate that some changes in model assumptions and parameter ranges should be made in the future to bring certain synthetic observables into better agreement with observations, it seems unlikely that these pending improvements will significantly affect the mass-loss rates of the models.

The heavy mass loss observed in evolved asymptotic giant branch (AGB) stars is usually attributed to a two-stage process: atmospheric levitation by pulsation-induced shock waves, followed by radiative acceleration of newly formed dust grains. The dust transfers momentum to the surrounding gas through collisions and thereby triggers a general outflow. Radiation-hydrodynamical models of M-type AGB stars suggest that these winds can be driven by photon scattering -- in contrast to absorption -- on Fe-free silicate grains of sizes 0.1--1\,$\mu$m. In this paper we study photometric constraints for wind-driving dust species in M-type AGB stars, as part of an ongoing effort to identify likely candidates among the grain materials observed in circumstellar envelopes. To investigate the scenario of stellar winds driven by photon scattering on dust, and to explore how different optical and chemical properties of wind-driving dust species affect photometry we focus on two sets of dynamical models atmospheres: (i) models using a detailed description for the growth of Mg$_2$SiO$_4$ grains, taking into account both scattering and absorption cross-sections when calculating the radiative acceleration, and (ii) models using a parameterized dust description, constructed to represent different chemical and optical dust properties. By comparing synthetic photometry from these two sets of models to observations of M-type AGB stars we can provide constraints on the properties of wind-driving dust species. Photometry from wind models with a detailed description for the growth of Mg$_2$SiO$_4$ grains reproduces well both the values and the time-dependent behavior of observations of M-type AGB stars, providing further support for the scenario of winds driven by photon scattering on dust.

We study the circumstellar environment of the M-type AGB star RT Vir using mid-infrared high spatial resolution observations from the ESO-VLTI focal instrument MIDI. The aim of this study is to provide observational constraints on theoretical prediction that the winds of M-type AGB objects can be driven by photon scattering on iron-free silicate grains located in the close environment (about 2 to 3 stellar radii) of the star. We interpreted spectro-interferometric data, first using wavelength-dependent geometric models. We then used a self-consistent dynamic model atmosphere containing a time-dependent description of grain growth for pure forsterite dust particles to reproduce the photometric, spectrometric, and interferometric measurements of RT Vir. Since the hydrodynamic computation needs stellar parameters as input, a considerable effort was first made to determine these parameters. MIDI differential phases reveal the presence of an asymmetry in the stellar vicinity. Results from the geometrical modeling give us clues to the presence of aluminum and silicate dust in the close circumstellar environment (< ~5 stellar radii). Comparison between spectro-interferometric data and a self-consistent dust-driven wind model reveals that silicate dust has to be present in the region between 2 to 3 stellar radii to reproduce the 59 and 63 m baseline visibility measurements around 9.8 micron. This gives additional observational evidence in favor of winds driven by photon scattering on iron-free silicate grains located in the close vicinity of an M-type star. However, other sources of opacity are clearly missing to reproduce the 10-13 micron visibility measurements for all baselines. This study is a first attempt to understand the wind mechanism of M-type AGB stars by comparing photometric, spectrometric, and interferometric measurements with state-of-the-art, self-consistent dust-driven wind models. The agreement of the dynamic model atmosphere with interferometric measurements in the 8-10 micron spectral region gives additional observational evidence that the winds of M-type stars can be driven by photon scattering on iron-free silicate grains. Finally, a larger statistical study and progress in advanced self-consistent 3D modeling are still required to solve the remaining problems.

The purpose of this corrigendum is to point out that a handful of models in the original paper were computed with faulty initial structures. Using exactly the same modelling methods we have recomputed the faulty models with new initial structures. The new results slightly changes some of the trends in the wind properties with stellar parameters, but the overall effects are small. The conclusions are not affected.

This work is part of an ongoing effort aiming at identifying the actual wind-drivers among the dust species observed in circumstellar envelopes. In particular, we focus on the interplay between a strong stellar radiation field and the dust formation process. To identify critical properties of potential wind-driving dust species we use detailed radiation-hydrodynamical models which include a parameterized dust description, complemented by simple analytical estimates to help with the physical interpretation of the numerical results. The adopted dust description is constructed to mimic different chemical and optical dust properties in order to systematically study the effects of a realistic radiation field on the second stage of the mass loss mechanism. We see distinct trends in which combinations of optical and chemical dust properties are needed to trigger an outflow. Dust species with a low condensation temperature and a NIR absorption coefficient that decreases strongly with wavelength will not condense close enough to the stellar surface to be considered as potential wind-drivers. Our models confirm that metallic iron and Fe-bearing silicates are not viable as wind-drivers due to their near-infrared optical properties and resulting large condensation distances. TiO2 is also excluded as a wind-driver due to the low abundance of Ti. Other species, such a SiO2 and Al2O3, are less clear-cut cases due to uncertainties in the optical and chemical data and further work is needed. A strong candidate is Mg2SiO4 with grain sizes of 0.1-1 micron, where scattering contributes significantly to the radiative acceleration, as suggested by earlier theoretical work and supported by recent observations.

[Abridged] In this paper we explore grain size effects on wind properties of carbon stars, using a generalized description of radiative cross sections valid for particles of arbitrary sizes. The purpose of the study is to investigate under which circumstances the small particle limit (SPL) may give acceptable results, and to quantify the possible errors that may occur when it does not hold. The time-dependent description of grain growth in our detailed radiation-hydrodynamical models gives information about dust particle radii in every layer at every instant of time. These grain radii are used for computing opacities and determining the radiative acceleration of the dust-gas mixture. It is shown that in the critical cases the effect of the generalized description of dust opacities can be significant, resulting in more intense mass loss and higher wind velocities compared to models using SPL opacities. For well-developed winds, however, grain size effects on mass loss rates and wind velocities are found to be small. Both groups of models tend towards lower degrees of dust condensation compared to corresponding SPL models, due to a self-regulating feedback between grain growth and radiative acceleration. Consequently, the "dust-loss rates" are lower in the models with the generalized treatment of grain opacities. We conclude that our previous results on mass loss rates obtained with SPL opacities are reliable within a wide region of stellar parameter space, except for critical cases close to thresholds of dust-driven outflows where SPL models will tend to underestimate the mass loss rates and wind velocities.

By using self-consistent dynamic model atmospheres which simulate pulsation-enhanced dust-driven winds of AGB stars we studied in detail the influence of (i) pulsations of the stellar interiors, and (ii) the development of dusty stellar winds on the spectral appearance of long period variables with carbon-rich atmospheric chemistry. While the pulsations lead to large-amplitude photometric variability, the dusty envelopes cause pronounced circumstellar reddening. Based on one selected dynamical model which is representative of C-type Mira variables with intermediate mass loss rates, we calculated synthetic spectra and photometry for standard broad-band filters from the visual to the near-infrared. Our modelling allows to investigate in detail the substantial effect of circumstellar dust on the resultant photometry. The pronounced absorption of amorphous carbon dust grains leads to colour indices which are significantly redder than the corresponding ones based on hydrostatic dust-free models. Only if we account for this circumstellar reddening we get synthetic colours that are comparable to observations of evolved AGB stars. The photometric variations of the dynamical model were compared to observed lightcurves of the C-type Mira RU_Vir which appears to be quite similar to the model. We found good agreement concerning the principal behaviour of the BVRIJHKL lightcurves and also quantitatively fitting details. The analysed model is able to reproduce the variations of RU_Vir and other Miras in (J-H) vs. (H-K) diagrams throughout the light cycle. Contrasting the model photometry with observational data for a variety of galactic C-rich giants in such colour-colour diagrams proved that the chosen atmospheric model fits well into a sequence of objects with increasing mass loss rates, i.e., redder colour indices.

There are strong observational indications that the dense slow winds of cool luminous AGB stars are driven by radiative pressure on dust grains which form in the extended atmospheres resulting from pulsation-induced shocks. For carbon stars, detailed models of outflows driven by amorphous carbon grains show good agreement with observations. Some still existing discrepancies may be due to a simplified treatment of cooling in shocks, drift of the grains relative to the gas, or effects of giant convection cells or dust-induced pattern formation. For stars with C/O < 1, recent models indicate that absorption by silicate dust is probably insuffcient to drive their winds. A possible alternative is scattering by Fe-free silicate grains with radii of a few tenths of a micron. In this scenario one should expect less circumstellar reddening for M- and S-type AGB stars than for C-stars with comparable stellar parameters and mass loss rates.

It is commonly believed that winds of cool giants in their late evolutionary stages are driven by radiative pressure on dust grains, but the actual grain species responsible for driving winds of M-type AGB stars are still a matter of debate (see H\"ofner, this volume). The optical properties of the dust particles responsible for initiating the outflow affect the stellar spectra in two ways: (i) indirectly, through their influence on the dynamical structure of the atmosphere/envelope and the resulting molecular features, and (ii) directly, by changes of the spectral energy distribution due to absorption and scattering on dust grains. The qualitative differences in the energy distributions of C-type and M-type AGB stars in the visual and near-infrared regions suggest that the dust particles in oxygen rich atmospheres are relatively transparent to radiation. By using detailed dynamical models of gas and radiation combined with a simple description for the dust opacity (which can be adjusted to mimic different wavelength dependencies and condensation temperatures) and also by adjusting the fraction of the opacity that is treated as true absorption, we investigate which dust properties produce synthetic photometry consistent with observations. The goal of this study is to narrow down the possible dust species that may be driving the winds in M-type AGB stars.

Mass loss is a fundamental, observationally well-established feature of AGB stars but many aspects of this process still remain to be understood. To date, self-consistent dynamical models of dust-driven winds reproducing the observed mass-loss rates seem only possible for M-type stars if the grains in the close circumstellar environment grow to larger sizes than previously assumed. In order to study the grain-size distribution where the mass loss is initiated, high-spatial-resolution interferometry observations are necessary. We have observed two M-type stars using the VLTI/MIDI instrument to constrain the dust-grain sizes through modeling the 10-micron silicate feature. Complementary observations are scheduled and we will present preliminary results.

We study the circumstellar environment of the carbon-rich star R Scl using the near- and mid-infrared high spatial resolution observations from the ESO-VLTI instruments VINCI and MIDI. These observations aim at increasing our knowledge of the dynamic processes in play within the very close circumstellar environment where the mass loss of AGB stars is initiated. Data are interpreted using a self-consistent dynamic model. Interferometric observations do not show any significant variability effect at the 16 m baseline between phases 0.17 and 0.23 in the K band, and for both the 15 m baseline between phases 0.66 and 0.97 and the 31 m baseline between phases 0.90 and 0.97 in the N band. We find fairly good agreement between the dynamic model and the spectrophotometric data from 0.4 to 25 $\mu$m. The model agrees well with the time-dependent flux data at 8.5 $\mu$m, whereas it is too faint at 11.3 and 12.5 $\mu$m. The VINCI visibilities are reproduced well, meaning that the extension of the model is suitable in the K-band. In the mid-infrared, the model has the proper extension to reveal molecular structures of C2H2 and HCN located above the stellar photosphere. However, the windless model used is not able to reproduce the more extended and dense dusty environment. Among the different explanations for the discrepancy between the model and the measurements, the strong nonequilibrium process of dust formation is one of the most probable. The complete dynamic coupling of gas and dust and the approximation of grain opacities with the small-particle limit in the dynamic calculation could also contribute to the difference between the model and the data.

We have computed wind models with time-dependent dust formation and grain-size dependent opacities, where (1) the problem is simplified by assuming a fixed dust-grain size, and where (2) the radiation pressure efficiency is approximated using grain sizes based on various means of the actual grain size distribution. It is shown that in critical cases, the effect of grain sizes can be significant. For well-developed winds, however, the effects on the mass-loss rate and the wind speed are small.

(abbreviated) Measuring the surface abundances of AGB stars is an important tool for studying the effects of nucleosynthesis and mixing in the interior of low- to intermediate mass stars during their final evolutionary phases. The atmospheres of AGB stars can be strongly affected by stellar pulsation and the development of a stellar wind, though, and the abundance determination of these objects should therefore be based on dynamic model atmospheres. We investigate the effects of stellar pulsation and mass loss on the appearance of selected spectral features (line profiles, line intensities) and on the derived elemental abundances by performing a systematic comparison of hydrostatic and dynamic model atmospheres. High-resolution synthetic spectra in the near infrared range were calculated based on two dynamic model atmospheres (at various phases during the pulsation cycle) as well as a grid of hydrostatic COMARCS models. Equivalent widths of a selection of atomic and molecular lines were derived in both cases and compared with each other. In the case of the dynamic models, the equivalent widths of all investigated features vary over the pulsation cycle. A consistent reproduction of the derived variations with a set of hydrostatic models is not possible, but several individual phases and spectral features can be reproduced well with the help of specific hydrostatic atmospheric models. In addition, we show that the variations in equivalent width that we found on the basis of the adopted dynamic model atmospheres agree qualitatively with observational results for the Mira R Cas over its light cycle. The findings of our modelling form a starting point to deal with the problem of abundance determination in strongly dynamic AGB stars (i.e., long-period variables).

We argue that the energy injection of pulsations may be of greater importance to the mass-loss rate of AGB stars than metallicity, and that the mass-loss trend with metallicity is not as simple as sometimes assumed. Using our detailed radiation hydrodynamical models that include dust formation, we illustrate the effects of pulsation energy on wind properties. We find that the mass-loss rate scales with the kinetic energy input by pulsations as long as a dust-saturated wind does not occur, and all other stellar parameters are kept constant. This includes the absolute abundance of condensible carbon (not bound in CO), which is more relevant than keeping the C/O-ratio constant when comparing stars of different metallicity. The pressure and temperature gradients in the atmospheres of stars, become steeper and flatter, respectively, when the metallicity is reduced, while the radius where the atmosphere becomes opaque is typically associated with a higher gas pressure. This effect can be compensated for by adjusting the velocity amplitude of the variable inner boundary (piston), which is used to simulate the effects of pulsation, to obtain models with comparable kinetic-energy input. Hence, it is more relevant to compare models with similar energy-injections than of similar velocity amplitude. Since there is no evidence for weaker pulsations in low-metallicity AGB stars, we conclude that it is unlikely that low-metallicity C-stars have a lower mass-loss rate, than their more metal-rich counterparts with similar stellar parameters, as long as they have a comparable amount of condensible carbon.

The stability and conservation properties of a recently proposed polymerization model are studied. The achiral (racemic) solution is linearly unstable once the relevant control parameter (here the fidelity of the catalyst) exceeds a critical value. The growth rate is calculated for different fidelity parameters and cross-inhibition rates. A chirality parameter is defined and shown to be conserved by the nonlinear terms of the model. Finally, a truncated version of the model is used to derive a set of two ordinary differential equations and it is argued that these equations are more realistic than those used in earlier models of that form.