Mira variables share essential characteristics: High visual amplitude, periods of hundreds of days, red colors (spectral types M, S, and C), and the presence of emission lines at some phases. They are fundamental mode pulsators, with progenitor masses ranging from >1 to several solar masses. In this review, we summarize what is known from modeling and observational studies, including recent measurements from optical and IR interferometry, and studies involving large samples of stars particularly in the Magellanic Clouds. While we have a good idea of how these stars fit into the big picture of stellar evolution, many important details remain to be settled by a combination of more ambitious models and new observational techniques. Carrying on observations of bright Mira variables will be essential for interpreting observations of large numbers of fainter sources as well as for assessing the completeness and accuracy of the models.

We present the initial-final mass relation derived from 10 white dwarfs in wide binaries that consist of a main sequence star and a white dwarf. The temperature and gravity of each white dwarf was measured by fitting theoretical model atmospheres to the observed spectrum using a $\chi^{2}$ fitting algorithm. The cooling time and mass was obtained using theoretical cooling tracks. The total age of each binary was estimated from the chromospheric activity of its main sequence component to an uncertainty of about 0.17 dex in log \textit{t} The difference between the total age and white dwarf cooling time is taken as the main sequence lifetime of each white dwarf. The initial mass of each white dwarf was then determined using stellar evolution tracks with a corresponding metallicity derived from spectra of their main sequence companions, thus yielding the initial-final mass relation. Most of the initial masses of the white dwarf components are between 1 - 2 M$_{\odot}$. Our results suggest a correlation between the metallicity of a white dwarf's progenitor and the amount of post-main-sequence mass loss it experiences - at least among progenitors with masses in the range of 1 - 2 M$_{\odot}$. A comparison of our observations to theoretical models suggests that low mass stars preferentially lose mass on the red giant branch.

We report imaging observations of the symbotic long-period Mira variable R Aquarii (R Aqr) at near-infrared and radio wavelengths. The near-infrared observations were made with the IOTA imaging interferometer in three narrow-band filters centered at 1.51, 1.64, and 1.78 $\mu$m, which sample mainly water, continuum, and water features, respectively. Our near-infrared fringe visibility and closure phase data are analyzed using three models. (a) A uniform disk model with wavelength-dependent sizes fails to fit the visibility data, and is inconsistent with the closure phase data. (b) A three- component model, comprising a Mira star, water shell, and an off-axis point source, provide a good fit to all data. (c) A model generated by a constrained image reconstruction analysis provides more insight, suggesting that the water shell is highly non-uniform, i.e., clumpy. The VLBA observations of SiO masers in the outer molecular envelope show evidence of turbulence, with jet-like features containing velocity gradients.

Mass loss rate formulae are derived from observations or from suites of models. For theoretical models, the following have all been identified as factors greatly influencing the atmospheric structure and mass loss rates: Pulsation with piston amplitude scaling appropriately with stellar L; dust nucleation and growth, with radiation pressure and grain-gas interactions and appropriate dependence on temperature and density; non-grey opacity with at least 51 frequency samples; non-LTE and departures from radiative equilibrium in the compressed and expanding flows; and non-equilibrium processes affecting the composition (grain formation; molecular chemistry). No one set of models yet includes all the factors known to be important. In fact, it is very difficult to construct a model that can simultaneously include these factors and be useful for computing spectra. Therefore, although theoretical model grids are needed to separate the effects of M,L,R and/or $T_{\mathrm{eff}}$ or Z on the mass loss rates, these models must be carefully checked against observations. Getting the right order of magnitude for the mass loss rate is only the first step in such a comparison, and is not sufficient to determine whether the mass loss formula is correct. However, there are observables that do test the validity of mass loss formulae as they depend directly on $d\log \dot M/d\log L$, $d\log \dot M/d\log R$, or $d\log \dot M/d\log P$.

We have measured non-zero closure phases for about 29% of our sample of 56 nearby Asymptotic Giant Branch (AGB) stars, using the 3-telescope Infrared Optical Telescope Array (IOTA) interferometer at near-infrared wavelengths (H band) and with angular resolutions in the range 5-10 milliarcseconds. These nonzero closure phases can only be generated by asymmetric brightness distributions of the target stars or their surroundings. We discuss how these results were obtained, and how they might be interpreted in terms of structures on or near the target stars. We also report measured angular sizes and hypothesize that most Mira stars would show detectable asymmetry if observed with adequate angular resolution.

Stellar evolution theory predicts that asymptotic giant branch stars undergo a series of short thermal pulses that significantly change their luminosity and mass on timescales of hundreds to thousands of years. Secular changes in these stars resulting from thermal pulses can be detected as measurable changes in period if the star is undergoing Mira pulsations. The American Association of Variable Star Observers (AAVSO) International Database currently contains visual data for over 1500 Mira variables. Light curves for these stars span nearly a century in some cases, making it possible to study the secular evolution of the pulsation behavior on these timescales. In this paper, we present the results of our study of period change in 547 Mira variables using data from the AAVSO. We find non-zero rates of period change, dlnP/dt, at the 2-sigma significance level in 57 of the 547 stars, at the 3-sigma level in 21 stars, and at the level of 6-sigma or greater in eight of the 547. The latter eight stars have been previously noted in the literature, and our derived rates of period changes largely agree with published values. The largest and most statistically significant dlnP/dt are consistent with the rates of period change expected during thermal pulses on the AGB. A number of other stars exhibit non-monotonic period changes on decades-long timescales, the cause of which is not yet known.

Numerical models of the extended atmospheres of long period variable or Mira stars have shown that their winds have a very simple, power law structure when averaged over the pulsation cycle. This structure is stable and robust despite the pulsational wave disturbances, and appears to be strongly self-regulated. Observational studies support these conclusions. The models also show that dust-free winds are nearly adiabatic, with little heating or cooling. The classical, steady, adiabatic wind solution to the hydrodynamic equations fails to account for an extensive region of nearly constant outflow velocity. We investigate analytic solutions which include the effects of wave pressure, heating, and the resulting entropy changes. Wave pressure is represented by a term like that in the Reynolds turbulence equation for the mean velocity. Although the pressure from individual waves is modest, the waves are likely the primary agent of self-regulation of dust-free winds. In models of dusty winds, the gas variables also adopt a power law dependence on radius. Heating is required at all radii to maintain this flow, and grain heating and heat transfer to the gas are significant. Both hydrodynamic and gas/grain thermal feedbacks transform the flow towards self-regulated forms. (Abridged)

We present numerical hydrodynamical modeling of the effects of a giant planet or brown dwarf companion orbiting within the extended atmosphere and wind formation zone of an approximately solar-mass Mira variable star. The large-scale, time-dependent accretion flows within the radially oscillating and outflowing circumstellar gas around Miras are related to Bondi-Hoyle-Lyttleton flows, but have not, to our knowledge, been previously modelled. The new models presented in this paper illustrate the changes in accretion and wake dynamics as the companion mass is varied over a range from 10 to 50 Jupiter masses ($M_J$). The character of the accretion onto the companion changes greatly as the companion mass is increased. At the lowest companion masses considered here, a low continuous rate of mass accretion is punctuated by large, nearly periodic bursts of accretion. When the companion mass is large, the mass accretion has both a continuous part, and a rapidly varying, nearly stochastic part. These trends can be understood as the result of the interplay between the shocks and radial oscillations in the circumstellar gas, and the wake flow behind the companion. Models with accretion bursts may produce observable optical brightenings, and may affect SiO maser emission. (Abridged)

We present numerical hydrodynamical models of the effects of planets or brown dwarfs orbiting within the extended atmosphere and wind formation zone of Mira variables. We find time-dependent wake dynamics and episodic accretion phenomena which may give rise to observable optical events and affect SiO maser emission.

Short time scale variations have been reported for a few Miras, including sudden 0.2 magnitude or more brightening in the visual lasting from a few hours to a few days. Are these flashes real? Are they hot? What are they? Almost all the natural time scales for variation in the atmosphere and wind of a Mira variable are months to years long. What could cause such short-term variations? One intriguing possibility is that the variations are associated with the interaction of a Jovian planet with the time-dependent outflow of a Mira wind. Here we discuss observable features of such an interaction in terms of order-of-magnitude estimates and phenomenology, and make one clear prediction requiring observational followup. Future work needs to include both theoretical calculations and the development of systematic methods for searching for such events.