Ultraslow radiative cooling lifetimes and adiabatic detachment energies for three astrochemically relevant anions, C$_n^-$ ($n=3-5$), are measured using the Double ElectroStatic Ion Ring ExpEriment (DESIREE) infrastructure at Stockholm University. DESIREE maintains a background pressure of $\approx$10$^{-14}$\,mbar and temperature of $\approx$13\,K, allowing storage of mass-selected ions for hours and providing conditions coined a "molecular cloud in a box". Here, we construct two-dimensional (2D) photodetachment spectra for the target anions by recording photodetachment signal as a function of irradiation wavelength and ion storage time (seconds to minute timescale). Ion cooling lifetimes, which are associated with infrared radiative emission, are extracted from the 2D photodetachment spectrum for each ion by tracking the disappearance of vibrational hot-band signal with ion storage time, giving $\frac{1}{e}$ cooling lifetimes of 3.1$\pm$0.1\,s (C$_3^-$), 6.8$\pm$0.5\,s (C$_4^-$) and 24$\pm$5\,s (C$_5^-$). Fits of the photodetachment spectra for cold ions, i.e. those stored for at least 30\,s, provides adiabatic detachment energies in good agreement with values from laser photoelectron spectroscopy. Ion cooling lifetimes are simulated using a Simple Harmonic Cascade model, finding good agreement with experiment and providing a mode-by-mode understanding of the radiative cooling properties. The 2D photodetachment strategy and radiative cooling modeling developed in this study could be applied to investigate the ultraslow cooling dynamics of wide range of molecular anions.

Substituted naphthalene anions (deprotonated 2-naphthol and 6-hydroxy-2-naphthoic acid) are spectroscopically probed in a tandem drift tube ion mobility spectrometer (IMS). Target anions are selected according to their drift speed through nitrogen buffer gas in the first IMS stage before being exposed to a pulse of tunable light that induces either photodissociation or electron photodetachment, which is conveniently monitored by scavenging the detached electrons with trace \ce{SF6} in the buffer gas. The photodetachment action spectrum of the 2-naphtholate anion exhibits a band system spanning 380-460\,nm~with a prominent series of peaks spaced by 440\,\cme, commencing at 458.5\,nm, and a set of weaker peaks near the electron detachment threshold corresponding to transitions to dipole-bound states. The two deprotomers of 6-hydroxy-2-naphthoic acid are separated and spectroscopically probed independently. The molecular anion formed from deprotonation of the hydroxy group possesses a photodetachment action spectrum similar to that of the 2-naphtholate anion with an onset at 470\,nm and a maximum at 420\,nm. Near threshold, photoreaction with \ce{SF6} is observed with displacement of an OH group by an F atom. In contrast, the anion formed from deprotonation of the carboxylic acid group features a photodissociation action spectrum, recorded on the \ce{CO2} loss channel, lying to much shorter wavelength with an onset at 360\,nm and maximum photoresponse at 325\,nm.

Dynamic model atmospheres of Mira stars predict variabilities in the photospheric radius and in atmospheric molecular layers which are not yet strongly constrained by observations. Here we measure the variability of the oxygen-rich Mira star R Peg. We used K band spectro-interferometry with a spectral resolution of about 4000 obtained at four epochs between post-maximum and minimum visual phases employing the newly available GRAVITY beam combiner at the VLTI. Our observations show a continuum radius that is anti-correlated with the visual lightcurve. Uniform disc (UD) angular diameters at a near-continuum wavelength of 2.25 mu are steadily increasing with values of 8.7 mas, 9.4 mas, 9.8 mas, and 9.9 mas at visual phases of 0.15, 0.36, 0,45, 0.53, respectively. UD diameters at a bandpass around 2.05 mu, dominated by water vapour, follow the near-continuum variability at larger UD diameters. UD diameters at the CO 2-0 bandhead, instead, are correlated with the visual lightcurve and anti-correlated with the near-continuum UD diameters. The observed anti-correlation between continuum radius and visual lightcurve is consistent with an earlier study of the oxygen-rich Mira S Lac, and with recent 1D CODEX dynamic model atmosphere predictions. The amplitude of the variation is comparable to the earlier observations of S Lac, and smaller than predicted by CODEX models. The wavelength-dependent visibility variations at our epochs can be reproduced by a set of CODEX models at model phases between 0.3 and 0.6. The anti-correlation of water vapour and CO contributions at our epochs suggests that these molecules undergo different processes along the stellar cycle. The newly available GRAVITY instrument is suited to conducting longer time series observations, which are needed to provide strong constraints on the model-predicted intra- and inter-cycle variability.

We study whether the condensation of silicate dust in Mira envelopes could be caused by cluster formation by the abundant SiO molecules. For a simplified model of the pulsational motions of matter in the the outer layers of a Mira variable which is guided by a numerical model for Mira pulsations, the equations of dust nucleation and growth are solved in the co-moving frame of a fixed mass element. It is assumed that seed particles form by clustering of SiO molecules. The calculation of the nucleation rate is based on the experimental data of Nuth and Donn (1982). The quantity of dust formed is calculated by a moment method and the calculation of radiation pressure on the dusty gas is based on a dirty silicate model. Dust nucleation occurs in the model at the upper culmination of the trajectory of a gas parcel where it stays for a considerable time at low temperatures while subsequent dust growth occurs during the descending part of the motion and continues after the next shock reversed motion. It is found that sufficient dust forms that radiation pressure exceeds gravitational pull of the stars such that the mass element is finally driven out of the star. Nucleation of dust particles by clustering of the abundant SiO molecules could be the mechanism that triggers silicate dust formation in Miras.

We obtained a total of 20 near-infrared K-band spectro-interferometric snapshot observations of the Mira variables o Cet, R Leo, R Aqr, X Hya, W Vel, and R Cnc with a spectral resolution of about 1500. We compared observed flux and visibility spectra with predictions by CODEX 1D dynamic model atmospheres and with azimuthally averaged intensities based on CO5BOLD 3D dynamic model atmospheres including convection. Our visibility data confirm the presence of spatially extended molecular atmospheres located above the continuum radii with large-scale inhomogeneities or clumps that contribute a few percent of the total flux. The detailed structure of the inhomogeneities or clumps show a variability on time scales of 3 months and above. Both modeling attempts provided satisfactory fits to our data. In particular, they are both consistent with the observed decrease in the visibility function at molecular bands of water vapor and CO, indicating a spatially extended molecular atmosphere. Observational variability phases are mostly consistent with those of the best-fit CODEX models, except for near-maximum phases, where data are better described by near-minimum models. Rosseland angular diameters derived from the model fits are broadly consistent between those based on the 1D and the 3D models and with earlier observations. We derived fundamental parameters including absolute radii, effective temperatures, and luminosities for our sources. Our results provide a first observational support for theoretical results that shocks induced by convection and pulsation in the 3D CO5BOLD models of AGB stars are roughly spherically expanding and of similar nature to those of self-excited pulsations in 1D CODEX models. Unlike for red supergiants, the pulsation- and shock-induced dynamics can levitate the molecular atmospheres of Mira variables to extensions that are consistent with observations.

The mass-loss process in Mira stars probably occurs in an asymmetric way where dust can form in inhomogeneous circumstellar molecular clumps. Following asymmetries along the pulsation cycle can give us clues about these mass-loss processes. We imaged the Mira star X Hya and its environnement at different epochs to follow the evolution of the morphology in the continuum and in the molecular bands. We observed X Hya with AMBER in J-H-K at low resolution at two epochs. We modelled squared visibilities with geometrical and physical models. We also present imaging reconstruction results obtained with MiRA and based on the physical a priori images. We report on the angular scale change of X Hya between the two epochs. 1D CODEX profiles allowed us to understand and model the spectral variation of squared visibilities and constrain the stellar parameters. Reconstructed model-dependent images enabled us to reproduce closure phase signals and the azimuthal dependence of squared visibilities. They show evidence for material inhomogeneities located in the immediate environment of the star.

We present the atmospheric structure and the fundamental parameters of three red supergiants, increasing the sample of RSGs observed by near-infrared spectro-interferometry. Additionally, we test possible mechanisms that may explain the large observed atmospheric extensions of RSGs. We carried out spectro-interferometric observations of 3 RSGs in the near-infrared K-band with the VLTI/AMBER instrument at medium spectral resolution. To comprehend the extended atmospheres, we compared our observational results to predictions by available hydrostatic PHOENIX, available 3-D convection, and new 1-D self-excited pulsation models of RSGs. Our near-infrared flux spectra are well reproduced by the PHOENIX model atmospheres. The continuum visibility values are consistent with a limb-darkened disk as predicted by the PHOENIX models, allowing us to determine the angular diameter and the fundamental parameters of our sources. Nonetheless, in the case of V602 Car and HD 95686, the PHOENIX model visibilities do not predict the large observed extensions of molecular layers, most remarkably in the CO bands. Likewise, the 3-D convection models and the 1-D pulsation models with typical parameters of RSGs lead to compact atmospheric structures as well, which are similar to the structure of the hydrostatic PHOENIX models. They can also not explain the observed decreases in the visibilities and thus the large atmospheric molecular extensions. The full sample of our RSGs indicates increasing observed atmospheric extensions with increasing luminosity and decreasing surface gravity, and no correlation with effective temperature or variability amplitude, which supports a scenario of radiative acceleration on Doppler-shifted molecular lines.

Aims. We study the effects of moderate deviations from solar abundances upon the atmospheric structure and colours of typical Mira variables. Methods. We present two model series of dynamical opacity-sampling models of Mira variables which have (1) 1 solar metallicity 3 and (2) "mild" S-type C/O abundance ratio ([C/O]=0.9) with typical Zr enhancement (solar +1.0). These series are compared to a previously studied solar-abundance series which has similar fundamental parameters (mass, luminosity, period, radius) that are close to those of o Cet. Results. Both series show noticeable effects of abundance upon stratifications and infrared colours but cycle-to-cycle differences mask these effects at most pulsation phases, with the exception of a narrow-water-filter colour near minimum phase.

We observed the AGB stars S Ori, GX Mon and R Cnc with the MIDI instrument at the VLTI. We compared the data to radiative transfer models of the dust shells, where the central stellar intensity profiles were described by dust-free dynamic model atmospheres. We used Al2O3 and warm silicate grains. Our S Ori and R Cnc data could be well described by an Al2O3 dust shell alone, and our GX Mon data by a mix of an Al2O3 and a silicate shell. The best-fit parameters for S Ori and R Cnc included photospheric angular diameters Theta(Phot) of 9.7+/-1.0mas and 12.3+/-1.0mas, optical depths tau(V)(Al2O3) of 1.5+/-0.5 and 1.35+/-0.2, and inner radii R(in) of 1.9+/-0.3R(Phot) and 2.2+/-0.3R(Phot), respectively. Best-fit parameters for GX Mon were Theta(Phot)=8.7+/-1.3mas, tau(V)(Al2O3)=1.9+/-0.6, R(in)(Al2O3)=2.1+/-0.3R(Phot), tau(V)(silicate)=3.2+/-0.5, and R(in)(silicate)=4.6+/-0.2R(Phot). Our model fits constrain the chemical composition and the inner boundary radii of the dust shells, as well as the photospheric angular diameters. Our interferometric results are consistent with Al2O3 grains condensing close to the stellar surface at about 2 stellar radii, co-located with the extended atmosphere and SiO maser emission, and warm silicate grains at larger distances of about 4--5 stellar radii. We verified that the number densities of aluminum can match that of the best-fit Al2O3 dust shell near the inner dust radius in sufficiently extended atmospheres, confirming that Al2O3 grains can be seed particles for the further dust condensation. Together with literature data of the mass-loss rates, our sample is consistent with a hypothesis that stars with low mass-loss rates form primarily dust that preserves the spectral properties of Al2O3, and stars with higher mass-loss rate form dust with properties of warm silicates.

Our ability to extract information from the spectra of stars depends on reliable models of stellar atmospheres and appropriate techniques for spectral synthesis. Various model codes and strategies for the analysis of stellar spectra are available today. We aim to compare the results of deriving stellar parameters using different atmosphere models and different analysis strategies. The focus is set on high-resolution spectroscopy of cool giant stars. Spectra representing four cool giant stars were made available to various groups and individuals working in the area of spectral synthesis, asking them to derive stellar parameters from the data provided. The results were discussed at a workshop in Vienna in 2010. Most of the major codes currently used in the astronomical community for analyses of stellar spectra were included in this experiment. We present the results from the different groups, as well as an additional experiment comparing the synthetic spectra produced by various codes for a given set of stellar parameters. Similarities and differences of the results are discussed. Several valid approaches to analyze a given spectrum of a star result in quite a wide range of solutions. The main causes for the differences in parameters derived by different groups seem to lie in the physical input data and in the details of the analysis method. This clearly shows how far from a definitive abundance analysis we still are.

Radio and infrared interferometry of SiO maser stars provide complementary information on the atmosphere and circumstellar environment at comparable spatial resolution. Here, we present the latest results on the atmospheric structure and the dust condensation region of AGB stars based on our recent infrared spectro-interferometric observations, which represent the environment of SiO masers. We discuss, as an example, new results from simultaneous VLTI and VLBA observations of the Mira variable AGB star R Cnc, including VLTI near- and mid-infrared interferometry, as well as VLBA observations of the SiO maser emission toward this source. We present preliminary results from a monitoring campaign of high-frequency SiO maser emission toward evolved stars obtained with the APEX telescope, which also serves as a precursor of ALMA images of the SiO emitting region. We speculate that large-scale long-period chaotic motion in the extended molecular atmosphere may be the physical reason for observed deviations from point symmetry of atmospheric molecular layers, and for the observed erratic variability of high-frequency SiO maser emission

We present 4 model series of the CODEX dynamical opacity-sampling models of Mira variables with solar abundances, designed to have parameters similar to $o$ Cet, R Leo and R Cas. We demonstrate that the CODEX models provide a clear physical basis for the molecular shell scenario used to explain interferometric observations of Mira variables. We show that these models generally provide a good match to photometry and interferometry at wavelengths between the near-infrared and the radio, and make the model outputs publicly available. These model also demonstrate that, in order to match visible and infrared observations, the Fe-poor silicate grains that form within 3 continuum radii must have small grain radii and therefore can not drive the winds from O-rich Mira variables.

We obtained 13 epochs of mid-infrared interferometry with the MIDI instrument at the VLTI between April 2004 and July 2007, covering pulsation phases 0.45-0.85 within four cycles. The data are modeled with a radiative transfer model of the dust shell where the central stellar intensity profile is described by a series of dust-free dynamic model atmospheres based on self-excited pulsation models. We examined two dust species, silicate and Al2O3 grains. We performed model simulations using variations in model phase and dust shell parameters to investigate the expected variability of our photometric and interferometric data. The observed visibility spectra do not show any indication of variations as a function of pulsation phase and cycle. The observed photometry spectra may indicate intracycle and cycle-to-cycle variations at the level of 1-2 standard deviations. The best-fitting model for our average pulsation phase of 0.64+/-0.15 includes the dynamic model atmosphere M21n (T_model=2550 K) with a photospheric angular diameter of 7.6+/-0.6 mas, and a silicate dust shell with an optical depth of 2.8+/-0.8, an inner radius of 4.1+/-0.7 R_Phot, and a power-law index of the density distribution of 2.6+/-0.3. The addition of an Al2O3 dust shell did not improve the model fit. The photospheric angular diameter corresponds to a radius of 520^+230_-140 R_sun and an effective temperature of ~ 2420+/-200 K. Our modeling simulations confirm that significant visibility variations are not expected for RR Aql at mid-infrared wavelengths within our uncertainties. We conclude that our RR Aql data can be described by a pulsating atmosphere surrounded by a silicate dust shell. The effects of the pulsation on the mid-infrared flux and visibility values are expected to be less than about 25% and 20%, respectively, and are too low to be detected within our measurement uncertainties.

We obtained K-band spectro-interferometric observations of the Miras R Cnc, X Hya, W Vel, and RW Vel with a spectral resolution of 1500 using the VLTI/AMBER instrument. We obtained concurrent JHKL photometry using the the Mk II instrument at the SAAO. Our sources have wavelength-dependent visibility values that are consistent with earlier low-resolution AMBER observations of S Ori and with the predictions of dynamic model atmosphere series based on self-excited pulsation models. The wavelength-dependent UD diameters show a minimum near the near-continuum bandpass at 2.25 um. They increase by up to 30% toward the H2O band at 2.0 um and by up to 70% at the CO bandheads. The dynamic model atmosphere series show a consistent wavelength-dependence, and their parameters such as the visual phase, effective temperature, and distances are consistent with independent estimates. The closure phases have significantly wavelength-dependent non-zero values indicating deviations from point symmetry. For example, the R Cnc closure phase is 110 degr in the 2.0 um H2O band, corresponding for instance to an additional unresolved spot contributing 3% of the total flux at a separation of ~4 mas. Our observations are consistent with the predictions of the latest dynamic model atmosphere series based on self-excited pulsation models. The wavelength-dependent radius variations are interpreted as the effect of molecular layers. The wavelength-dependent closure phase values are indicative of deviations from point symmetry at all wavelengths, thus a complex non-spherical stratification of the extended atmosphere. In particular, the significant deviation from point symmetry in the H2O band is interpreted as a signature on large scales of inhomogeneities or clumps in the water vapor layer. The observed inhomogeneities might be caused by pulsation- and shock-induced chaotic motion in the extended atmosphere.

We present recent studies using the near-infrared instrument AMBER of the VLT Interferometer (VLTI) to investigate the structure and shaping processes within the extended atmosphere of AGB stars. Spectrally resolved near-infrared AMBER observations of the Mira variable S Ori have revealed wavelength-dependent apparent angular sizes. These data were successfully compared to dynamic model atmospheres, which predict wavelength-dependent radii because of geometrically extended molecular layers. Most recently, AMBER closure phase measurements of several AGB stars have also revealed wavelength-dependent deviations from 0/180 deg., indicating deviations from point symmetry. The variation of closure phase with wavelength indicates a complex non-spherical stratification of the extended atmosphere, and may reveal whether observed asymmetries are located near the photosphere or in the outer molecular layers. Concurrent observations of SiO masers located within the extended molecular layers provide us with additional information on the morphology, conditions, and kinematics of this shell. These observations promise to provide us with new important insights into the shaping processes at work during the AGB phase. With improved imaging capabilities at the VLTI, we expect to extend the successful story of imaging studies of planetary nebulae to the photosphere and extended outer atmosphere of AGB stars.

We present an overview on our project to study the extended atmospheres and dust formation zones of Mira stars using coordinated observations with the Very Large Telescope Interferometer (VLTI), the Very Long Baseline Array (VLBA), and the Atacama Pathfinder Experiment (APEX). The data are interpreted using an approach of combining recent dynamic model atmospheres with a radiative transfer model of the dust shell, and combining the resulting model structure with a maser propagation model.

Aims. We aim to explore the photosphere of the very cool late-type star VX Sgr and in particular the existence and characterization of molecular layers above the continuum forming photosphere. Methods. We obtained interferometric observations with the VLTI/AMBER interferometer using the fringe tracker FINITO in the spectral domain 1.45-2.50 micron with a spectral resolution of about 35 and baselines ranging from 15 to 88 meters.We perform independent image reconstruction for different wavelength bins and fit the interferometric data with a geometrical toy model.We also compare the data to 1D dynamical models of Miras atmosphere and to 3D hydrodynamical simulations of red supergiant (RSG) and asymptotic giant branch (AGB) stars. Results. Reconstructed images and visibilities show a strong wavelength dependence. The H-band images display two bright spots whose positions are confirmed by the geometrical toy model. The inhomogeneities are qualitatively predicted by 3D simulations. At about 2,00 micron and in the region 2,35 - 2,50 micron, the photosphere appears extended and the radius is larger than in the H band. In this spectral region, the geometrical toy model locates a third bright spot outside the photosphere that can be a feature of the molecular layers. The wavelength dependence of the visibility can be qualitatively explained by 1D dynamical models of Mira atmospheres. The best-fitting photospheric models show a good match with the observed visibilities and give a photospheric diameter of theta = 8,82+-0,50 mas. The H2O molecule seems to be the dominant absorber in the molecular layers. Conclusions. We show that the atmosphere of VX Sgr rather resembles Mira/AGB star model atmospheres than RSG model atmospheres. In particular, we see molecular (water) layers that are typical for Mira stars.

We describe a combined dynamic atmosphere and maser propagation model of SiO maser emission in Mira variables. This model rectifies many of the defects of an earlier model of this type, particularly in relation to the infra-red (IR) radiation field generated by dust and various wavelength-dependent, optically thick layers. Modelled masers form in rings with radii consistent with those found in VLBI observations and with earlier models. This agreement requires the adoption of a radio photosphere of radius approximately twice that of the stellar photosphere, in agreement with observations. A radio photosphere of this size renders invisible certain maser sites with high amplification at low radii, and conceals high-velocity shocks, which are absent in radio continuum observations. The SiO masers are brightest at an optical phase of 0.1 to 0.25, which is consistent with observed phase-lags. Dust can have both mild and profound effects on the maser emission. Maser rings, a shock and the optically thick layer in the SiO pumping band at 8.13\micron appear to be closely associated in three out of four phase samples.

We present results from a spectro-interferometric study of the Miras o Cet, R Leo and W Hya obtained with the Keck Aperture Masking Experiment from 1998 Sep to 2002 Jul. The spectrally dispersed visibility data permit fitting with circularly symmetric brightness profiles such as a simple uniform disk. The stellar angular diameter obtained over up to ~ 450 spectral channels spaning the region 1.1-3.8 microns is presented. Use of a simple uniform disk brightness model facilitates comparison between epochs and with existing data and theoretical models. Strong size variations with wavelength were recorded for all stars, probing zones of H2O, CO, OH, and dust formation. Comparison with contemporaneous spectra extracted from our data show a strong anti-correlation between the observed angular diameter and flux. These variations consolidate the notion of a complex stellar atmosphere consisting of molecular shells with time-dependent densities and temperatures. Our findings are compared with existing data and pulsation models. The models were found to reproduce the functional form of the wavelength vs. angular diameter curve well, although some departures are noted in the 2.8-3.5 micron range.

We describe the Cool Opacity-sampling Dynamic EXtended (CODEX) atmosphere models of Mira variable stars, and examine in detail the physical and numerical approximations that go in to the model creation. The CODEX atmospheric models are obtained by computing the temperature and the chemical and radiative states of the atmospheric layers, assuming gas pressure and velocity profiles from Mira pulsation models, which extend from near the H-burning shell to the outer layers of the atmosphere. Although the code uses the approximation of Local Thermodynamic Equilibrium (LTE) and a grey approximation in the dynamical atmosphere code, many key observable quantities, such as infrared diameters and low-resolution spectra, are predicted robustly in spite of these approximations. We show that in visible light, radiation from Mira variables is dominated by fluorescence scattering processes, and that the LTE approximation likely under-predicts visible-band fluxes by a factor of two.

We have used near- and mid-infrared interferometry to investigate the pulsating atmosphere and the circumstellar environment of the Mira variable RR Aql. Observations were taken with the VLTI/AMBER (near infrared) and the VLTI/MIDI (mid infrared) instruments. We have obtained a total of 15 MIDI epochs between Apr 9, 2004 and Jul 28, 2007 covering 4 pulsation cycles and one AMBER epoch on Sep 9, 2006 at phase 2.82. This work is also part of an ongoing project of joint VLTI and VLBA observations to study the connection between stellar pulsation and the mass loss process. Here we present a comparison of the AMBER visibility data to a simple uniform disk model as well as to predictions by recent self-excited dynamic model atmospheres. The best fitting photospheric angular diameter of the model atmosphere at phase 2.82 is 9.9 +/- 2.4 mas.

Aims: We present J, H, K interferometry with a spectral resolution of 35 for the Mira variable S Orionis. We aim at measuring the diameter variation as a function of wavelength that is expected due to molecular layers lying above the continuum-forming photosphere. Methods: Visibility data of S Ori were obtained at phase 0.78 with the VLTI/AMBER instrument using the fringe tracker FINITO at 29 spectral channels between 1.29 and 2.32 mu. Apparent uniform disk (UD) diameters were computed for each spectral channel. In addition, the visibility data were directly compared to predictions by recent self-excited dynamic model atmospheres. Results: S Ori shows significant variations in the visibility values as a function of spectral channel that can only be described by a clear variation in the apparent angular size with wavelength. The closure phase values are close to zero at all spectral channels, indicating the absence of asymmetric intensity features. The apparent UD angular diameter is smallest at about 1.3 and 1.7 mu and increases by a factor of ~1.4 around 2.0 mu. The minimum UD angular diameter is 8.1 pm 0.5 mas, corresponding to ~420 R_sun. The S Ori visibility data and the apparent UD variations can be explained reasonably well by a dynamic atmosphere model that includes molecular layers. Conclusions: The measured visibility and UD diameter variations with wavelength resemble and generally confirm the predictions by recent dynamic model atmospheres. [abridged]

The angular diameters of six oxygen rich Mira-type long-period variables have been measured at various near-infrared (NIR) wavelengths using the aperture masking technique in an extensive observing program from 1997 Jan to 2004 Sep. These data sets span many pulsation cycles of the observed objects and represent the largest study of multi-wavelength, multi-epoch interferometric angular diameter measurements on Mira stars to date. The calibrated visibility data of o Cet, R Leo, R Cas, W Hya, chi Cyg and R Hya are fitted using a uniform disk brightness distribution model to facilitate comparison between epochs, wavelengths and with existing data and theoretical models. The variation of angular diameter as a function of wavelength and time are studied, and cyclic diameter variations are detected for all objects in our sample. These variations are believed to stem from time-dependent changes of density and temperature (and hence varying molecular opacities) in different layers of these stars. The similarities and differences in behaviour between these objects are analyzed and discussed in the context of existing theoretical models. Furthermore, we present time-dependent 3.08 micron angular diameter measurements, probing for the first time these zones of probable dust formation, which show unforeseen sizes and are consistently out of phase with other NIR layers shown in this study. The S-type Mira chi Cyg exhibits significantly different behaviour compared to the M-type Miras in this study.

We present the first multi-epoch study that includes concurrent mid-infrared and radio interferometry of an oxygen-rich Mira star. We obtained mid-infrared interferometry of S Ori with VLTI/MIDI at four epochs between December 2004 and December 2005. We concurrently observed v=1, J=1-0 (43.1 GHz), and v=2, J=1-0 (42.8 GHz) SiO maser emission toward S Ori with the VLBA at three epochs. The MIDI data are analyzed using self-excited dynamic model atmospheres including molecular layers, complemented by a radiative transfer model of the circumstellar dust shell. The VLBA data are reduced to the spatial structure and kinematics of the maser spots. The modeling of our MIDI data results in phase-dependent continuum photospheric angular diameters between about 7.9 mas (Phase 0.55) and 9.7 mas (Phase 1.16). The dust shell can best be modeled with Al2O3 grains using phase-dependent inner boundary radii between 1.8 and 2.4 photospheric radii. The dust shell appears to be more compact with greater optical depth near visual minimum, and more extended with lower optical depth after visual maximum. The ratios of the SiO maser ring radii to the photospheric radii are between about 1.9 and 2.4. The maser spots mark the region of the molecular atmospheric layers just beyond the steepest decrease in the mid-infrared model intensity profile. Their velocity structure indicates a radial gas expansion. Al2O3 dust grains and SiO maser spots form at relatively small radii of 1.8-2.4 photospheric radii. Our results suggest increased mass loss and dust formation close to the surface near the minimum visual phase, when Al2O3 dust grains are co-located with the molecular gas and the SiO maser shells, and a more expanded dust shell after visual maximum. Silicon does not appear to be bound in dust, as our data show no sign of silicate grains.

High dispersion near-infrared spectra have been taken of seven highly-evolved, variable, intermediate-mass (4-6 Msun) AGB stars in the LMC and SMC in order to look for C, N and O variations that are expected to arise from third dredge-up and hot-bottom burning. The pulsation of the objects has been modelled, yielding stellar masses, and spectral synthesis calculations have been performed in order to derive abundances from the observed spectra. For two stars, abundances of C, N, O, Na, Al, Ti, Sc and Fe were derived and compared with the abundances predicted by detailed AGB models. Both stars show very large N enhancements and C deficiencies. These results provide the first observational confirmation of the long-predicted production of primary nitrogen by the combination of third dredge-up and hot-bottom burning in intermediate-mass AGB stars. It was not possible to derive abundances for the remaining five stars: three were too cool to model, while another two had strong shocks in their atmospheres which caused strong emission to fill the line cores and made abundance determination impossible. The latter occurrence allows us to predict the pulsation phase interval during which observations should be made if successful abundance analysis is to be possible.

The formation of dust with temperature-dependent non-grey opacity is considered in a series of self-consistent model atmospheres at different phases of an O-rich Mira variable of mass 1.2 $M_\odot$. Photometric and interferometric properties of these models are predicted under different physical assumptions regarding the dust formation. The iron content of the initial silicate that forms and the availability of grain nuclei are found to be critical parameters that affect the observable properties. In particular, parameters were found where dust would form at 2-3 times the average continuum photospheric radius. This work provides a consistent physical explanation for the larger apparent size of Mira variables at wavelengths shorter than 1 $\mu$m than that predicted by dust free fundamental-mode pulsation models.

We present a comparison of dynamical model atmospheres with mid-infrared (~11 micron) interferometric and spectroscopic observations of the Mira variable o Cet. The dynamical model atmospheres of Mira variables pulsating in the fundamental mode can fairly explain, without assuming ad-hoc components, the seemingly contradictory mid-infrared spectroscopic and interferometric observations of o Cet: the 11 micron sizes measured in the bandpass without any salient spectral features are about twice as large as those measured in the near-infrared. Our calculations of synthetic spectra show that the strong absorption due to a number of optically thick H2O lines is filled in by the emission of these H2O lines originating in the geometrically extended layers, providing a possible physical explanation for the picture proposed by Ohnaka (2004a) based on a semi-empirical modeling. This filling-in effect results in rather featureless, continuum-like spectra in rough agreement with the observed high-resolution 11 micron spectra, although the models still predict the H2O lines to be more pronounced than the observations. The inverse P-Cyg profiles of some strong H2O lines observed in the 11 micron spectra can also be reasonably reproduced by our dynamical model atmospheres. The presence of the extended H2O layers manifests itself as mid-infrared angular diameters much larger than the continuum diameter. The 11 micron uniform-disk diameters predicted by our dynamical model atmospheres are in fair agreement with those observed with the Infrared Spatial Interferometer (ISI), but still somewhat smaller than the observed diameters. We discuss possible reasons for this discrepancy and problems with the current dynamical model atmospheres of Mira variables.

We present near-infrared K-band interferometric measurements of the Mira star R Leonis obtained in April 2001 and January 2002 with the VLTI, the commissioning instrument VINCI, and the two test siderostats. These epochs correspond to near-maximum stellar variability phases $\sim$ 0.08 and $\sim$ 1.02 (one cycle later), respectively. The April 2001 data cover a range of spatial frequencies (31--35 cycles/arcsecond) within the first lobe of the visibility function. These measurements indicate a center-to-limb intensity variation (CLV) that is clearly different from a uniform disk (UD) intensity profile. We show that these measured visibility values are consistent with predictions from recent self-excited dynamic Mira model atmospheres that include molecular shells close to continuum-forming layers. We derive high-precision Rosseland diameters of 28.5 +/- 0.4 mas and 26.2 +/- 0.8 mas for the April 2001 and January 2002 data, respectively. Together with literature estimates of the distance and the bolometric flux, these values correspond to linear radii of 350^{+50}_{-40}R$_\odot$ and 320^{+50}_{-40}R$_\odot$, and to effective temperatures of 2930 +/- 270K and 3080 +/- 310K, respectively.

Model-predicted and observed properties of the brightness distribution on M-type Mira disks are discussed. Fundamental issues of limb-darkening and diameter definition, of assigning observational data to diameter-type quantities and of interpreting such quantities in terms of model diameters are outlined. The influence of model properties upon interpretation of measured data is clarified. The dependence of the centre-to-limb variation (CLV) of intensity on wavelength, on stellar parameters and on variablity phase and cycle may be used for analyzing the geometrical and physical structure of the Mira atmosphere, for determining fundamental stellar parameters, and for investigating the quality of models. Desirable future observations include simultaneous observations in different spectral features at different phases and cycles, observation of the position of the shock front and observation of the time- and wavelength-dependence of deviations from spherical symmetry.

Models of M-type Miras with masses of 1 $M_\odot$ and 1.2 $M_\odot$, i.e. with envelope masses of about 0.4 $M_\odot$ and 0.6 $M_\odot$, have been constructed, and a comparison has been made of their observable properties. Geometric pulsation of continuum-forming layers is found to be little affected by the mass difference. The influence of molecular contamination of near-infrared continuum bandpasses upon interferometrically measured fit diameters ranges from undetectable to quite significant. Some pulsation cycles of the lower-mass model Mira show substantially stronger contamination than that found in any cycle of the higher-mass star. Observations which sample pulsation phase well and continuously are crucial for avoiding misinterpretations, because the assignment of absolute pulsation phases is inherently uncertain by at least 0.1 cycles, diameter changes may be strongly phase-dependent, and cycle-to-cycle variations may be substantial. In accord with expectations, we find that cycle-to-cycle variations that show up in light curves and in near-continuum diameters tend to be larger and more common in the low-mass models, leading to one possible way to discriminate mass. Two other methods, based on high-precision measurements of the pulsation amplitude and on derivation of pre-maximum effective temperatures from diameter measurements, are also discussed. High-layer features that may be strongly affected by mass are not well described by present dust-free models.

We present K-band commissioning observations of the Mira star prototype o Cet obtained at the ESO Very Large Telescope Interferometer (VLTI) with the VINCI instrument and two siderostats. The observations were carried out between 2001 October and December, in 2002 January and December, and in 2003 January. Rosseland angular radii are derived from the measured visibilities by fitting theoretical visibility functions obtained from center-to-limb intensity variations (CLVs) of Mira star models (Bessel et al. 1996, Hofmann et al. 1998, Tej et al. 2003). Using the derived Rosseland angular radii and the SEDs reconstructed from available photometric and spectrophotometric data, we find effective temperatures ranging from T_eff=3192 +/- 200 K at phase phi=0.13 to 2918 +/- 183 K at phi=0.26. Comparison of these Rosseland radii, effective temperatures, and the shape of the observed visibility functions with model predictions suggests that o Cet is a fundamental mode pulsator. Furthermore, we investigated the variation of visibility function and diameter with phase. The Rosseland angular diameter of o Cet increased from 28.9 +/- 0.3 mas (corresponding to a Rosseland radius of 332 +/- 38 R_sun for a distance of D=107 +/- 12 pc) at phi=0.13 to 34.9 +/- 0.4 mas (402 +/- 46 R_sun) at phi=0.4. The error of the Rosseland linear radius almost entirely results from the error of the parallax, since the error of the angular diameter is only approximately 1 %.

Pulsation models are crucial for the interpretation of spectrophotometric and interferometric observations of Mira variables. Comparing predicted and observed spectra is one way of establishing the validity of such models. In this paper, we focus on the models published between 1996 and 1998 by Bessell, Hofmann, Scholz and Wood. A few new model spectra are added, to improve available phase coverage. We compare the synthetic spectra with observed low resolution spectra of optically selected oxygen-rich Miras, over a range of optical and near-IR wavelengths that encompasses most of the stellar energy output. We investigate the overall energy distributions, and specific spectral features in the near-IR wavelength range. The agreement between the observed and the model-predicted properties is found to be reasonably good. However, there are discrepancies seen especially in various molecular bands. We find that different combinations of stellar parameters and pulsation phases often result in very similar model spectra. Therefore the problem of deriving parameters of a Mira variable from its spectrum has no unique solution. More advanced models than presently available, providing even better fits to the data and covering a wider range of parameters, would be needed to achieve better discrimination.

Dynamic models of M-type Mira variables predict the occurrence of water "shells", i.e. of zones of high H2O density and high H2O absorption inside the stellar atmosphere. The density, position and width of these shells is closely correlated with different types of two-component shapes of the intensity distribution on the disk in the H, K and L near-continuum bandpasses. We investigate these correlations and highlight the role of a spectrophotometric H2O index that warns against serious complications in diameter measurements in the case of substantial water contamination of the bandpass of observation. Simultaneous spectrophotometric and interferometric measurements may allow observers to estimate real continuum diameters more precisely.

We present K'-band observations of five Mira stars with the IOTA interferometer. The interferograms were obtained with the FLUOR fiber optics beam combiner, which provides high-accuracy visibility measurements in spite of time-variable atmospheric conditions. For the M-type Miras X Oph, R Aql, RU Her, R Ser, and the C-type Mira V CrB we derived the uniform-disk diameters 11.7mas, 10.9mas, 8.4mas, 8.1mas, and 7.9mas (+/- 0.3mas), respectively. Simultaneous photometric observations yielded the bolometric fluxes. The derived angular Rosseland radii and the bolometric fluxes allowed the determination of effective temperatures. For instance, the effective temperature of R Aql was determined to be 2970 +/- 110 K. A linear Rosseland radius for R Aql of (250 +100/-60) Rsun was derived from the angular Rosseland radius of 5.5mas +/- 0.2mas and the HIPPARCOS parallax of 4.73mas +/- 1.19mas. The observations were compared with theoretical Mira star models of Bessel et al. (1996) and Hofmann et al. (1998). The effective temperatures of the M-type Miras and the linear radius of R Aql indicate fundamental mode pulsation.

We present visibility measurements of the nearby Mira-like star R Doradus taken over a wide range of wavelengths (650-990 nm). The observations were made using MAPPIT (Masked APerture-Plane Interference Telescope), an interferometer operating at the 3.9-m Anglo-Australian Telescope. We used a slit to mask the telescope aperture and prism to disperse the interference pattern in wavelength. We observed in R Dor strong decreases in visibility within the TiO absorption bands. The results are in general agreement with theory but differ in detail, suggesting that further work is needed to refine the theoretical models.

We present visibility measurements of the nearby Mira-like star R Doradus taken over a wide range of wavelengths (650--990 nm). The observations were made using MAPPIT (Masked APerture-Plane Interference Telescope), an interferometer operating at the 3.9-m Anglo-Australian Telescope. We used a slit to mask the telescope aperture and prism to disperse the interference pattern in wavelength. We observed in R Dor strong decreases in visibility within the TiO absorption bands. The results are in general agreement with theory but differ in detail, suggesting that further work is needed to refine the theoretical models.

We have built an infrared beam combiner for the GI2T/REGAIN interferometer of the Observatoire de la Cote d'Azur. The beam combiner allows us to record spectrally dispersed Michelson interference fringes in the near-infrared J-, H- or K-bands. The beam combiner has the advantage that Michelson interferograms can simultaneously be recorded in about 128 different spectral channels. The tilt of the spectrally dispersed fringes is a measure of the instantaneous optical path difference. We present the optical design of the beam combiner and GI2T/REGAIN observations of the Mira star R Cas with this beam combiner in the spectral range of 2.00 micron - 2.18 micron (observations on 22 and 25 August 1999; variability phase 0.08; V-magnitude approx. 6; seven baselines between 12m and 24m; reference stars Vega and Beta Peg). The spectrograph of the beam combiner consists of an anamorphotic cylindrical lens system, an image plane slit, and a grism. A system of digital signal processors calculates the ensemble average power spectrum of the spectrally dispersed Michelson interferograms and the instantaneous optical path difference error in real time. From the observed R Cas visibilities at baselines 12.0m, 13.8m and 13.9m, a 2.1 micron uniform-disk diameter of 25.3mas +/-3.3mas was derived. The unusually high visibility values at baselines >16m show that the stellar surface of R Cas is more complex than previously assumed. The visibility values at baselines >16m can be explained by high-contrast surface structure on the stellar surface of R Cas or other types of unexpected center-to-limb variations. The R Cas observations were compared with theoretical Mira star models yielding a linear Rosseland radius of 276Rsun +/-66Rsun and an effective temperature of 2685K+/-238K for R Cas at phase 0.08.

We present K-band observations of five Mira stars with the IOTA interferometer. The interferograms were obtained with the FLUOR fiber optics beam combiner which provides high-accuracy visibility measurements in spite of time-variable atmospheric conditions. For the Mira stars X Oph, R Aql, RU Her, R Ser, and V CrB we derived the uniform-disk diameters 11.7 mas, 10.9 mas, 8.4 mas, 8.1 mas, and 7.9 mas (+/-0.3 mas), respectively. Simultaneous photometric observations yielded the bolometric fluxes. The derived angular Rosseland radii and the bolometric fluxes allowed the determination of effective temperatures. For instance, the effective temperature of R Aql was determined to be 3072 K +/- 161 K. A Rosseland radius for R Aql of 250 Rsun +/- 63 Rsun was derived from the angular Rosseland radius of 5.5 mas +/- 0.2 mas and the HIPPARCOS parallax of 4.73 mas +/- 1.19 mas. The observations were compared with theoretical Mira star models of Bessel, Scholz and Wood (1996) and Hofmann, Scholz and Wood (1998).

We discuss the wavelength dependence of angular diameters of M giants from an observational perspective. Observers cannot directly measure an optical-depth radius for a star, despite this being a common theoretical definition. Instead, they can use an interferometer to measure the square of the fringe visibility. We present new plots of the wavelength-dependent centre-to-limb variation (CLV) of intensity of the stellar disk as well as visibility for Mira and non-Mira M giant models. We use the terms ``CLV spectra'' and ``visibility spectra'' for these plots. We discuss a model-predicted extreme limb-darkening effect (also called the narrow-bright-core effect) in very strong TiO bands which can lead to a misinterpretation of the size of a star in these bands. We find no evidence as yet that this effect occurs in real stars. Our CLV spectra can explain the similarity in visibilities of R Dor (M8IIIe) that have been observed recently despite the use of two different passbands. We compare several observations with models and find the models generally under-estimate the observed variation in visibility with wavelength. We present CLV and visibility spectra for a model that is applicable to the M supergiant alpha Ori.