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.

The atmospheres and winds of C-type AGB stars are modeled with the 1D spherically symmetric radiation-hydrodynamical code DARWIN. To explore the metallicity-dependence of mass loss we calculate model grids at three different chemical abundances. Since carbon may be dredged up during the thermal pulses as AGB stars evolve, we keep the carbon abundance as a free parameter. The models in these three different grids all have a current mass of one solar mass; effective temperatures of 2600K, 2800K, 3000K, or 3200K; and stellar luminosities equal to log(L/Lsun)=3.70, 3.85, or 4.00. The models show that mass loss in carbon stars is facilitated by high luminosities, low effective temperatures, and a high carbon excess (C-O) at both solar and subsolar metallicities. Similar combinations of effective temperature, luminosity, and carbon excess produce outflows at both solar and subsolar metallicities. There are no large systematic differences in the mass-loss rates and wind velocities produced by these wind models with respect to metallicity, nor any systematic difference concerning the distribution of grain sizes or how much carbon is condensed into dust. DARWIN models at subsolar metallicity have approximately 15% lower mass-loss rates compared to DARWIN models at solar metallicity with the same stellar parameters and carbon excess. For both solar and subsolar environments typical grain sizes range between 0.1 and 0.5 micron, the degree of condensed carbon varies between 5% and 40%, and the gas-to-dust ratios between 500 and 10000. Furthermore, stellar evolution models can use the mass-loss rates calculated from DARWIN models at solar metallicity when modeling the AGB phase at subsolar metallicities if carbon excess is used as the critical abundance parameter instead of the C/O ratio.

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 symbiotic star R Aqr is part of a small sample of binary AGB stars observed with ALMA. The sample stars are: R Aqr, Mira, W Aql, and pi1 Gru. The sample covers a range in binary separation and wind properties, where R Aqr is the source with the smallest separation. The R Aqr binary pair consists of an M-type AGB star and a white dwarf at a separation of 45 mas, equivalent to about 10 AU at 218 pc. The aim of the ALMA study is to investigate the dependence of the wind shaping on the binary separation and to provide constraints for hydrodynamical binary interaction models. R Aqr is particularly interesting as the source with the smallest separation and a complex circumstellar environment that is strongly affected by the interaction between the two stars and by the high-energy radiation resulting from this interaction and from the hot white dwarf companion. The CO(J=3-2) line emission has been observed at 0.5" spatial resolution. The CO envelope around the binary pair is marginally resolved, showing what appears to be a rather complex distribution. The outer radius of the CO emitting region is estimated from the data and found to be about a factor of 10 larger than previously thought. This implies an average mass-loss rate during the past 100 yr of 2x10-7 Msun/yr, a factor of 45 less than previous estimates. The channel maps are presented and the molecular gas distribution is discussed and set into the context of what was previously known about the system. Additional molecular line emission detected within the bandwidth covered by the ALMA observations is also presented. Because of the limited extent of the emission, firm conclusions about the dynamical evolution of the system will have to wait for higher spatial resolution observations. However, the data presented here support the assumption that the mass-loss rate from the Mira star strongly varies and is focused on the orbital plane.

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.

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.

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.