How to disentangle geometry and mass-loss rate from AGB-star spectral energy distributions -- The case of EP Aqr

High-angular-resolution observations of asymptotic giant branch (AGB) stars often reveal non-spherical morphologies for the gas and dust envelopes. We aim to make a pilot study to quantify the impact of different geometries (spherically symmetric, spiral-shaped, and disc-shaped) of the dust component of AGB envelopes on spectral energy distributions (SEDs), mass estimates, and subsequent mass-loss rate (MLR) estimates. We also estimate the error made on the MLR if the SED is fitted by an inappropriate geometrical model. We use the 3D Monte-Carlo-based radiative-transfer code RADMC-3D to simulate emission from dusty envelopes with different geometries (but fixed spatial extension). We compare these predictions with each other, and with the SED of the AGB star EP Aqr that we use as a benchmark since its envelope is disc-like and known to harbour spiral arms, as seen in CO. The SEDs involving the most massive envelopes are those for which the different geometries have the largest impact, primarily on the silicate features at 10 and 18 um. These different shapes originate from large differences in optical depths. Massive spirals and discs appear akin to black bodies. Optically thick edge-on spirals and discs (with dust masses of 1e-4 and 1e-5 Msun) exhibit black-body SEDs that appear cooler than those from face-on structures and spheres of the same mass, while optically thick face-on distributions appear as warmer emission. We find that our more realistic models, combined spherical and spiral distributions, are 0.1 to 0.5 times less massive than spheres with similar SEDs. More extreme, less realistic scenarios give that spirals and discs are 0.01 to 0.05 times less massive than corresponding spheres. This means that adopting the wrong geometry for an AGB circumstellar envelope may result in a MLR that is incorrect by as much as 1 to 2 orders of magnitude when derived from SED fitting.

The flux-weighted gravity-luminosity relation of Galactic classical Cepheids

The flux-weighted gravity-luminosity relation (FWGLR) is investigated for a sample of 477 classical Cepheids (CCs), including stars that have been classified in the literature as such but are probably not. The luminosities are taken from the literature, based on the fitting of the spectral energy distributions (SEDs) assuming a certain distance and reddening. The flux-weighted gravity (FWG) is taken from gravity and effective temperature determinations in the literature based on high-resolution spectroscopy. There is a very good agreement between the theoretically predicted and observed FWG versus pulsation period relation that could serve in estimating the FWG (and $log g$) in spectroscopic studies with a precision of 0.1~dex. As was known in the literature, the theoretically predicted FWGLR relation for CCs is very tight and is not very sensitive to metallicity (at least for LMC and solar values), rotation rate, and crossing of the instability strip. The observed relation has a slightly different slope and shows more scatter (0.54~dex). This is due both to uncertainties in the distances and to the pulsation phase averaged FWG values. Data from future Gaia data releases should reduce these errors, and then the FWGLR could serve as a powerful tool in Cepheid studies.

The VMC Survey -- XXXVII. Pulsation periods of dust enshrouded AGB stars in the Magellanic Clouds

(abridged)Variability is a key property of stars on the asymptotic giant branch (AGB). Their pulsation period is related to the luminosity and mass-loss rate of the star. The long-period variables (LPVs) and Mira variables are the most prominent of all types of variability of evolved stars. The reddest, most obscured AGB stars are too faint in the optical and have eluded large variability surveys. Selection criteria are derived based on colour-colour and colour-magnitude diagrams from the combination of VISTA Magellanic Cloud (VMC) survey, Spitzer IRAC and AllWISE data. After eliminating LPVs with known periods shorter than 450 days, a sample of 1299 candidate obscured AGB stars is selected. K-band LCs are constructed combining the epoch photometry available in the VMC survey with literature data, analysed for variability and fitted with a single period sine curve to derive mean magnitudes, amplitudes and periods. A subset of 254 stars are either new variables, or known variables where the period we find is better determined than the literature value, or longer than 1000 days. The spectral energy distributions (SEDs) of these stars are fitted to a large number of templates. For this purpose the SEDs and Spitzer IRS spectra of some non-AGB stars are also fitted to have templates of the most likely contaminants in the sample. A sample of 217 likely LPVs is found. Thirty-four stars have periods longer than 1000 days although some of them have alternative shorter periods. Mass-loss rates and luminosities are estimated from the template fitting. Period-luminosity relations are presented for C- and O-rich Miras that appear to be extensions of relations derived in the literature for shorter periods. The fit for the C-stars is particularly well defined (with 182 objects) and reads Mbol = (-2.27 pm 0.20) log P + (1.45 pm 0.54)mag with an rms of 0.41 mag.

Analysing the spectral energy distributions of Galactic classical Cepheids

(abridged) Spectral energy distributions (SEDs) were constructed for a sample of 477 classical cepheids (CCs). The SEDs were fitted with a dust radiative transfer code. Four stars showed a large mid- or far-infrared excess and the fitting then included a dust component. These comprise the well-known case of RS Pup, and three stars that are (likely) Type-II cepheids (T2Cs), AU Peg, QQ Per, and FQ Lac. The remainder of the sample was fitted with a stellar photosphere to derive the best-fitting luminosity and effective temperature. Distance and reddening were taken from the literature. The stars were plotted in a Hertzsprung-Russell diagram and compared to evolutionary tracks for cepheids and theoretical instability strips. For the large majority of stars, the position in the HRD is consistent with the instability strip for a CC or T2C. About 5% of the stars are outliers in the sense that they are much hotter or cooler than expected. A comparison to effective temperatures derived from spectroscopy suggests in some cases that the photometrically derived temperature is not correct and that this is likely linked to an incorrectly adopted reddening. In this work the presence of a small NIR excess, as has been proposed in the literature for a few well-known cepheids, is investigated. Firstly, this was done by using a sample of about a dozen stars for which a mid-infrared spectrum is available. Secondly, the SEDs of all stars were fitted with a dust model to see if a statistically significant better fit could be obtained. The results were compared to recent work. Eight new candidates for exhibiting a NIR excess are proposed, solely based on the photometric SEDs. Obtaining mid-infrared spectra would be needed to confirm this excess. Finally, period-bolometric luminosity and period-radius relations are presented for samples of over 370 fundamental-mode CCs.

The VMC Survey - XXXIII. The tip of the red giant branch in the Magellanic Clouds

In this paper JK_s data from the VISTA Magellanic Cloud (VMC) survey are used to investigate the tip of the red giant branch (TRGB) as a distance indicator. A linear fit to recent theoretical models is used which reads M_{K_s} = -4.196 -2.013 (J-K_s), valid in the colour range 0.75 < (J-K_s) < 1.3 mag and in the 2MASS system. The observed TRGB is found based on a classical first-order and a second-order derivative filter applied to the binned luminosity function using the sharpened magnitude that takes the colour term into account. Extensive simulations are carried out to investigate any biases and errors in the derived distance modulus (DM). Based on these simulations criteria are established related to the number of stars per bin in the 0.5 magnitude range below the TRGB and related to the significance with which the peak in the filter response curve is determined such that the derived distances are unbiased. The DMs based on the second-order derivative filter are found to be more stable and are therefore adopted, although this requires twice as many stars per bin. The TRGB method is applied to specific lines-of-sight where independent distance estimates exist, based on detached eclipsing binaries in the LMC and SMC, classical Cepheids in the LMC, RR Lyrae stars in the SMC, and fields in the SMC where the star formation history (together with reddening and distance) has been derived from deep VMC data. The analysis shows that the theoretical calibration is consistent with the data, that the systematic error on the DM is approximately 0.045 mag, and that random errors of 0.015 mag are achievable. Reddening is an important element in deriving the distance: we find mean DMs ranging from 18.92 (for a typical E(B-V) of 0.15 mag) to 19.07 mag (E(B-V) about 0.04) for the SMC, and ranging from 18.48 (E(B-V) about 0.12 mag) to 18.57 mag (E(B-V) about 0.05) for the LMC.

The Cepheid period -- luminosity -- metallicity relation based on Gaia DR2 data

We use parallax data from the Gaia second data release (GDR2), combined with parallax data based on Hipparcos and HST data, to derive the period-luminosity-metallicity (PLZ) relation for Galactic classical cepheids (CCs) in the V,K, and Wesenheit WVK bands. An initial sample of 452 CCs are extracted from the literature with spectroscopically derived iron abundances. Reddening values, pulsation periods, and mean magnitudes are taken from the literature. Based on nine CCs with a goodness-of-fit (GOF) statistic <8 and with an accurate non-Gaia parallax, a parallax zero-point offset of -0.049 +- 0.018 mas is derived. Selecting a GOF statistic <8 removes about 40% of the sample most likely related due to binarity. Excluding first overtone and multi-mode cepheids and applying some other criteria reduces the sample to about 200 stars. The derived PL(Z) relations depend strongly on the parallax zero-point offset. The slope of the PL relation is found to be different from the relations in the LMC at the 3 sigma level. Fixing the slope to the value found in the LMC leads to a distance modulus (DM) to the LMC of order 18.7 mag, larger than the canonical distance. The canonical DM of around 18.5 mag would require a parallax zero-point offset of order $-0.1$ mas. Given the strong correlation between zero point, period and metallicity dependence of the PL relation, and the parallax zero-point offset there is no evidence for a metallicity term in the PLZ relation. The GDR2 release does not allow us to improve on the current distance scale based on CCs. The value of and the uncertainty on the parallax zero-point offset leads to uncertainties of order 0.15 mag on the distance scale. The parallax zero-point offset will need to be known at a level of 3 microas or better to have a 0.01 mag or smaller effect on the zero point of the PL relation and the DM to the LMC.

PACS and SPIRE range spectroscopy of cool, evolved stars

Context: At the end of their lives AGB stars are prolific producers of dust and gas. The details of this mass-loss process are still not understood very well. Herschel PACS and SPIRE spectra offer a unique way of investigating properties of AGB stars in general and the mass-loss process in particular. Methods: The HIPE software with the latest calibration is used to process the available PACS and SPIRE spectra of 40 evolved stars. The spectra are convolved with the response curves of the PACS and SPIRE bolometers and compared to the fluxes measured in imaging data of these sources. Custom software is used to identify lines in the spectra, and to determine the central wavelengths and line intensities. Standard molecular line databases are used to associate the observed lines. Because of the limited spectral resolution of the spectrometers several known lines are typically potential counterparts to any observed line. To help identifications the relative contributions in line intensity of the potential counterpart lines are listed for three characteristic temperatures based on LTE calculations and assuming optically thin emission. Result: The following data products are released: the reduced spectra, the lines that are measured in the spectra with wavelength, intensity, potential identifications, and the continuum spectra, i.e. the full spectra with all identified lines removed. As simple examples of how this data can be used in future studies we have fitted the continuum spectra with three power laws and find that the few OH/IR stars seem to have significantly steeper slopes than the other oxygen- and carbon-rich objects in the sample. As another example we constructed rotational diagrams for CO and fitted a two-component model to derive rotational temperatures.

The photodissociation of CO in circumstellar envelopes

Carbon monoxide is the most abundant molecule after H$_2$ and is important for chemistry in circumstellar envelopes around late-type stars. The size of the envelope is important when modelling low-J transition lines and deriving mass-loss rates from such lines. Now that ALMA is coming to full power the extent of the CO emitting region can be measured directly for nearby asymptotic giant branch (AGB) stars. In parallel, it has become obvious in the past few years that the strength of the interstellar radiation field (ISRF) can have a significant impact on the interpretation of the emission lines. In this paper an update and extension of the classical Mamon et al. (1988; ApJ 328, 797) paper is presented; these authors provided the CO abundance profile, described by two parameters, as a function of mass-loss rate and expansion velocity. Following recent work an improved numerical method and updated H$_2$ and CO shielding functions are used and a larger grid is calculated that covers more parameter space, including the strength of the ISRF. The effect of changing the photodissociation radius on the low-J CO line intensities is illustrated in two cases.

Dusty shells surrounding the carbon variables S Scuti and RT Capricorni

For the Mass-loss of Evolved StarS (MESS) programme, the unprecedented spatial resolution of the PACS photometer on board the Herschel space observatory was employed to map the dusty environments of asymptotic giant branch (AGB) and red supergiant (RSG) stars. Among the morphologically heterogeneous sample, a small fraction of targets is enclosed by spherically symmetric detached envelopes. Based on observations in the 70 {mu}m and 160 {mu}m wavelength bands, we investigated the surroundings of the two carbon semiregular variables S Sct and RT Cap, which both show evidence for a history of highly variable mass-loss. S Sct exhibits a bright, spherically symmetric detached shell, 138 in diameter and co-spatial with an already known CO structure. Moreover, weak emission is detected at the outskirts, where the morphology seems indicative of a mild shaping by interaction of the wind with the interstellar medium, which is also supported by the stellar space motion. Two shells are found around RT Cap that were not known so far in either dust emission or from molecular line observations. The inner shell with a diameter of 188 shows an almost immaculate spherical symmetry, while the outer ~5 structure is more irregularly shaped. MoD, a modification of the DUSTY radiative transfer code, was used to model the detached shells. Dust temperatures, shell dust masses, and mass-loss rates are derived for both targets.

MESS (Mass-loss of Evolved StarS), a Herschel Key Program

MESS (Mass-loss of Evolved StarS) is a Guaranteed Time Key Program that uses the PACS and SPIRE instruments on board the Herschel Space Observatory to observe a representative sample of evolved stars, that include asymptotic giant branch (AGB) and post-AGB stars, planetary nebulae and red supergiants, as well as luminous blue variables, Wolf-Rayet stars and supernova remnants. In total, of order 150 objects are observed in imaging and about 50 objects in spectroscopy. This paper describes the target selection and target list, and the observing strategy. Key science projects are described, and illustrated using results obtained during Herschels science demonstration phase. Aperture photometry is given for the 70 AGB and post-AGB stars observed up to October 17, 2010, which constitutes the largest single uniform database of far-IR and sub-mm fluxes for late-type stars.