We report the first detections of OH$^+$ emission in planetary nebulae (PNe). As part of an imaging and spectroscopy survey of 11 PNe in the far-IR using the PACS and SPIRE instruments aboard the Herschel Space Observatory, we performed a line survey in these PNe over the entire spectral range between 51 and 672$\mu$m to look for new detections. OH$^+$ rotational emission lines at 152.99, 290.20, 308.48, and 329.77$\mu$m were detected in the spectra of three planetary nebulae: NGC 6445, NGC 6720, and NGC 6781. Excitation temperatures and column densities derived from these lines are in the range of 27 to 47 K and 2$\times$10$^{10}$ to 4 $\times$10$^{11}$ cm$^{-2}$, respectively. In PNe, the OH+ rotational line emission appears to be produced in the photodissociation region (PDR) in these objects. The emission of OH+ is observed only in PNe with hot central stars (T$_{eff}$ > 100000 K), suggesting that high-energy photons may play a role in the OH+ formation and its line excitation in these objects, as it seems to be the case for ultraluminous galaxies.

This is the first of a series of investigations into far-IR characteristics of 11 planetary nebulae (PNs) under the Herschel Space Observatory Open Time 1 program, Herschel Planetary Nebula Survey (HerPlaNS). Using the HerPlaNS data set, we look into the PN energetics and variations of the physical conditions within the target nebulae. In the present work, we provide an overview of the survey, data acquisition and processing, and resulting data products. We perform (1) PACS/SPIRE broadband imaging to determine the spatial distribution of the cold dust component in the target PNs and (2) PACS/SPIRE spectral-energy-distribution (SED) and line spectroscopy to determine the spatial distribution of the gas component in the target PNs. For the case of NGC 6781, the broadband maps confirm the nearly pole-on barrel structure of the amorphous carbon-richdust shell and the surrounding halo having temperatures of 26-40 K. The PACS/SPIRE multi-position spectra show spatial variations of far-IR lines that reflect the physical stratification of the nebula. We demonstrate that spatially-resolved far-IR line diagnostics yield the (T_e, n_e) profiles, from which distributions of ionized, atomic, and molecular gases can be determined. Direct comparison of the dust and gas column mass maps constrained by the HerPlaNS data allows to construct an empirical gas-to-dust mass ratio map, which shows a range of ratios with the median of 195+-110. The present analysis yields estimates of the total mass of the shell to be 0.86 M_sun, consisting of 0.54 M_sun of ionized gas, 0.12 M_sun of atomic gas, 0.2 M_sun of molecular gas, and 4 x 10^-3 M_sun of dust grains. These estimates also suggest that the central star of about 1.5 M_sun initial mass is terminating its PN evolution onto the white dwarf cooling track.

Context. The interaction between stellar winds and the interstellar medium (ISM) can create complex bow shocks. The photometers on board the Herschel Space Observatory are ideally suited to studying the morphologies of these bow shocks. Aims. We aim to study the circumstellar environment and wind-ISM interaction of the nearest red supergiant, Betelgeuse. Methods. Herschel PACS images at 70, 100, and 160 micron and SPIRE images at 250, 350, and 500 micron were obtained by scanning the region around Betelgeuse. These data were complemented with ultraviolet GALEX data, near-infrared WISE data, and radio 21 cm GALFA-HI data. The observational properties of the bow shock structure were deduced from the data and compared with hydrodynamical simulations. Results. The infrared Herschel images of the environment around Betelgeuse are spectacular, showing the occurrence of multiple arcs at 6-7 arcmin from the central target and the presence of a linear bar at 9 arcmin. Remarkably, no large-scale instabilities are seen in the outer arcs and linear bar. The dust temperature in the outer arcs varies between 40 and 140 K, with the linear bar having the same colour temperature as the arcs. The inner envelope shows clear evidence of a non-homogeneous clumpy structure (beyond 15 arcsec), probably related to the giant convection cells of the outer atmosphere. The non-homogeneous distribution of the material even persists until the collision with the ISM. A strong variation in brightness of the inner clumps at a radius of 2 arcmin suggests a drastic change in mean gas and dust density some 32 000 yr ago. Using hydrodynamical simulations, we try to explain the observed morphology of the bow shock around Betelgeuse. Conclusions: [abbreviated]

CW Leo has been observed six times between October 2009 and June 2012 with the SPIRE instrument on board the Herschel satellite. Variability has been detected in the flux emitted by the central star with a period of 639 \pm 4 days, in good agreement with determinations in the literature. Variability is also detected in the bow shock around CW Leo that had previously been detected in the ultraviolet and Herschel PACS/SPIRE data. Although difficult to prove directly, our working hypothesis is that this variability is directly related to that of the central star. In this case, fitting a sine curve with the period fixed to 639 days results in a time-lag in the variability between bow shock and the central star of 402 \pm 37 days. The orientation of the bow shock relative to the plane of the sky is unknown (but see below). For an inclination angle of zero degrees, the observed time-lag translates into a distance to CW Leo of 130 \pm 13 pc, and for non-zero inclination angles the distance is smaller. Fitting the shape of the bow shock with an analytical model (Wilkin 1996), the effect of the inclination angle on the distance may be estimated. Making the additional assumption that the relative peculiar velocity between the interstellar medium (ISM) and CW Leo is determined entirely by the star space velocity with respect to the local standard of rest (i.e. a stationary ISM), the inclination angle is found to be (-33.3 \pm 0.8) degrees based on the observed proper motion and radial velocity. Using the Wilkin model, our current best estimate of the distance to CW Leo is 123 \pm 14 pc. For a distance of 123 pc, we derive a mean luminosity of 7790 \pm 150 Lsol (internal error).

We present Herschel Space Observatory photometric observations of the unique, long-period eclipsing binary star Epsilon Aurigae. Its extended spectral energy distribution is consistent with our previously published cool (550 K) dust disk model. We also present an archival infrared spectral energy distribution of the side of the disk facing the bright F-type star in the binary, which is consistent with a warmer (1150 K) disk model. The lack of strong molecular emission features in the Herschel bands suggests that the disk has a low gas-to-dust ratio. The spectral energy distribution and Herschel images imply that the 250 GHz radio detection reported by Altenhoff et al. is likely contaminated by infrared-bright, extended background emission associated with a nearby nebular region and should be considered an upper limit to the true flux density of Epsilon Aur.

S-type AGB stars are thought to be in the transitional phase between M-type and C-type AGB stars. Because of their peculiar chemical composition, one may expect a strong influence of the stellar C/O ratio on the molecular chemistry and the mineralogy of the circumstellar dust. In this paper, we present a large sample of 87 intrinsic galactic S-type AGB stars, observed at infrared wavelengths with the Spitzer Space Telescope, and supplemented with ground-based optical data. On the one hand, we derive the stellar parameters from the optical spectroscopy and photometry, using a grid of model atmospheres. On the other, we decompose the infrared spectra to quantify the flux-contributions from the different dust species. Finally, we compare the independently determined stellar parameters and dust properties. For the stars without significant dust emission, we detect a strict relation between the presence of SiS absorption in the Spitzer spectra and the C/O ratio of the stellar atmosphere. These absorption bands can thus be used as an additional diagnostic for the C/O ratio. For stars with significant dust emission, we define three groups, based on the relative contribution of certain dust species to the infrared flux. We find a strong link between group-membership and C/O ratio. We show that these groups can be explained by assuming that the dust-condensation can be cut short before silicates are produced, while the remaining free atoms and molecules can then form the observed magnesium sulfides or the carriers of the unidentified 13 and 20 micron features. Finally, we present the detection of emission features attributed to molecules and dust characteristic to C-type stars, such as molecular SiS, hydrocarbons and magnesium sulfide grains. We show that we often detect magnesium sulfides together with molecular SiS and we propose that it is formed by a reaction of SiS molecules with Mg.

Far-infrared Herschel/PACS images at 70 and 160 micron of a sample of 78 Galactic evolved stars are used to study the (dust) emission structures, originating from stellar wind-ISM interaction. In addition, two-fluid hydrodynamical simulations of the coupled gas and dust in wind-ISM interactions are used to compare with the observations. Four distinct classes of wind-ISM interaction (i.e. "fermata", "eyes", "irregular", and "rings") are identified and basic parameters affecting the morphology are discussed. We detect bow shocks for ~40% of the sample and detached rings for ~20%. De-projected stand-off distances (R_0) -- defined as the distance between the central star and the nearest point of the interaction region -- of the detected bow shocks ("fermata" and "eyes") are derived from the PACS images and compared to previous results, model predictions and the simulations. All observed bow shocks have stand-off distances smaller than 1 pc. Observed and theoretical stand-off distances are used together to independently derive the local ISM density. Both theoretical (analytical) models and hydrodynamical simulations give stand-off distances for adopted stellar properties that are in good agreement with the measured de-projected stand-off distance of wind-ISM bow shocks. The possible detection of a bow shock -- for the distance limited sample -- appears to be governed by its physical size as set roughly by the stand-off distance. In particular the star's peculiar space velocity and the density of the ISM appear decisive in detecting emission from bow shocks or detached rings. Tentatively, the "eyes" class objects are associated to (visual) binaries, while the "rings" generally appear not to occur for M-type stars, only for C or S-type objects that have experienced a thermal pulse.

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 Herschel's 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.

Detached circumstellar dust shells are detected around three carbon variables using Herschel-PACS. Two of them are already known on the basis of their thermal CO emission and two are visible as extensions in IRAS imaging data. By model fits to the new data sets, physical sizes, expansion timescales, dust temperatures, and more are deduced. A comparison with existing molecular CO material shows a high degree of correlation for TT Cyg and U Ant but a few distinct differences with other observables are also found.

Herschel PACS and SPIRE images have been obtained over a 30'x30' area around the well-known carbon star CW Leo (IRC +10 216). An extended structure is found in an incomplete arc of ~22' diameter, which is cospatial with the termination shock due to interaction with the interstellar medium (ISM) as defined by Sahai & Chronopoulos from ultraviolet GALEX images. Fluxes are derived in the 70, 160, 250, 350, and 550 um bands in the region where the interaction with the ISM takes place, and this can be fitted with a modified black body with a temperature of 25+-3 K. Using the published proper motion and radial velocity for the star, we derive a heliocentric space motion of 25.1 km/s. Using the PACS and SPIRE data and the analytical formula of the bow shock structure, we infer a de-projected standoff distance of the bow shock of R0 = (8.0+-0.3)x10^17 cm. We also derive a relative velocity of the star with respect to the ISM of (106.6+-8.7)/sqrt(n_ISM) km/s, where n_ISM is the number density of the local ISM.

As a result of the variability survey in Chi Persei and NGC6910, the number of Beta Cep stars that are members of these two open clusters is increased to twenty stars, nine in NGC6910 and eleven in Chi Persei. We compare pulsational properties, in particular the frequency spectra, of Beta Cep stars in both clusters and explain the differences in terms of the global parameters of the clusters. We also indicate that the more complicated pattern of the variability among B type stars in Chi Persei is very likely caused by higher rotational velocities of stars in this cluster. We conclude that the sample of pulsating stars in the two open clusters constitutes a very good starting point for the ensemble asteroseismology of Beta Cep-type stars and maybe also for other B-type pulsators.

During their evolution, Asymptotic Giant Branch (AGB) stars experience a high mass-loss which leads to the formation of a Circumstellar Envelope (CSE) of dust and gas. The mass-loss process is the most important phenomenon during this evolutionary stage. In order to understand it, it is important to study the physical parameters of the CSE. The emission of the CSE in the (sub)millimetre range is dominated by the dust continuum. This means that (sub)millimetre observations are a key tool in tracing the dust and improving our knowledge of the mass-loss process. In this paper we analyse new sub-millimetre observations of 9 evolved stars in order to constrain the CSE physical parameters. The data were taken by the new APEX bolometer LABoCa. The fluxes at 870 um are derived and the extended emission is investigated. We compute the spectral energy distribution (SEDs) using a 1D radiative transfer code, DUSTY which we compared to literature data. Grain properties are calculated using both spherical grains distribution and a Continuous Distribution of Ellipsoids (CDE) and a comparison between the two is drawn. Synthetic surface brightness maps have been derived from the modelling and were compared to the LABoCa brightness maps. We detected the presence of extended emission around four stars. Physical parameters of the circumstellar envelope are derived from SED modelling, such as the dust chemical composition, the dust condensation temperature and the total mass of the envelope. It proves difficult however to fit the SED and the intensity profile simultaneously.

CONTEXT: Recent progress in the seismic interpretation of field beta Cep stars has resulted in improvements of the physics in the stellar structure and evolution models of massive stars. Further asteroseismic constraints can be obtained from studying ensembles of stars in a young open cluster, which all have similar age, distance and chemical composition. AIMS: To improve our comprehension of the beta Cep stars, we studied the young open cluster NGC 884 to discover new B-type pulsators, besides the two known beta Cep stars, and other variable stars. METHODS: An extensive multi-site campaign was set up to gather accurate CCD photometry time series in four filters (U, B, V, I) of a field of NGC884. Fifteen different instruments collected almost 77500 CCD images in 1286 hours. The images were calibrated and reduced to transform the CCD frames into interpretable differential light curves. Various variability indicators and frequency analyses were applied to detect variable stars in the field. Absolute photometry was taken to deduce some general cluster and stellar properties. RESULTS: We achieved an accuracy for the brightest stars of 5.7 mmag in V, 6.9 mmag in B, 5.0 mmag in I and 5.3 mmag in U. The noise level in the amplitude spectra is 50 micromag in the V band. Our campaign confirms the previously known pulsators, and we report more than one hundred new multi- and mono-periodic B-, A- and F-type stars. Their interpretation in terms of classical instability domains is not straightforward, pointing to imperfections in theoretical instability computations. In addition, we have discovered six new eclipsing binaries and four candidates as well as other irregular variable stars in the observed field.

Accurate photometric CoRoT space observations of a secondary seismological target, HD 174884, led to the discovery that this star is an astrophysically important double-lined eclipsing spectroscopic binary in an eccentric orbit (e of about 0.3), unusual for its short (3.65705d) orbital period. The high eccentricity, coupled with the orientation of the binary orbit in space, explains the very unusual observed light curve with strongly unequal primary and secondary eclipses having the depth ratio of 1-to-100 in the CoRoT 'seismo' passband. Without the high accuracy of the CoRoT photometry, the secondary eclipse, 1.5 mmag deep, would have gone unnoticed. A spectroscopic follow-up program provided 45 high dispersion spectra. The analysis of the CoRoT light curve was performed with an adapted version of PHOEBE that supports CoRoT passbands. The final solution was obtained by simultaneous fitting of the light and the radial velocity curves. Individual star spectra were derived by spectrum disentangling. The uncertainties of the fit were derived by bootstrap resampling and the solution uniqueness was tested by heuristic scanning. The results provide a consistent picture of the system composed of two late B stars. The Fourier analysis of the light curve fit residuals yields two components, with orbital frequency multiples and an amplitude of about 0.1 mmag, which are tentatively interpreted as tidally induced pulsations. An extensive comparison with theoretical models is carried out by means of the Levenberg-Marquardt minimization technique and the discrepancy between models and the derived parameters is discussed. The best fitting models yield a young system age of 125 million years which is consistent with the eccentric orbit and synchronous component rotation at periastron.

While the first binary post-AGB stars were serendipitously discovered, the distinct characteristics of their Spectral Energy Distribution (SED) allowed us to launch a more systematic search for binaries. We selected post-AGB objects which show a broad dust excess often starting already at H or K, pointing to the presence of a gravitationally bound dusty disc in the system. We started a very extensive multi-wavelength study of those systems and here we report on our radial velocity and photometric monitoring results for six stars of early F type, which are pulsators of small amplitude. To determine the radial velocity of low signal-to-noise time-series, we constructed dedicated auto-correlation masks. The radial velocity variations were subjected to detailed analysis to differentiate between pulsational variability and variability due to orbital motion. Finally orbital minimalisation was performed to constrain the orbital elements. All of the six objects are binaries, with orbital periods ranging from 120 to 1800 days. Five systems have non-circular orbits. The mass functions range from 0.004 to 0.57 solar mass and the companions are likely unevolved objects of (very) low initial mass. We argue that these binaries must have been subject to severe binary interaction when the primary was a cool supergiant. Although the origin of the circumstellar disc is not well understood, the disc is generally believed to be formed during this strong interaction phase. The eccentric orbits of these highly evolved objects remain poorly understood. With the measured orbits and mass functions we conclude that the circumbinary discs seem to have a major impact on the evolution of a significant fraction of binary systems.

Like other starburst galaxies, M82 hosts compact, massive young star clusters that are interesting both in their own right and as benchmarks for population synthesis models. Can spectral synthesis models at resolutions around 1000 adequately reproduce the near-IR spectral features and the energy distribution of these clusters between 0.8 and 2.4 microns? How do the derived cluster properties compare with previous results from optical studies? We analyse the spectra of 5 massive clusters in M82, using data acquired with the spectrograph SpeX on the InfraRed Telescope Facility (NASA/IRTF) and a new population synthesis tool with a highly improved near-IR extension, based on a recent collection of empirical and theoretical spectra of red supergiant stars. We obtain excellent fits across the near-IR with models at quasi-solar metallicity and a solar neighbourhood extinction law. Spectroscopy breaks a strong degeneracy between age and extinction in the near-IR colours in the red supergiant-dominated phase of evolution. The estimated near-IR ages cluster between 9 and 30 Myr, i.e. the ages at which the molecular bands due to luminous red supergiants are strongest in the current models. They do not always agree with optical spectroscopic ages. Adding optical data sometimes leads to the rejection of the solar neighbourhood extinction law. This is not surprising considering small-scale structure around the clusters, but it has no significant effect on the near-IR based spectroscopic ages. [abridged]

We provide a grid of PHOENIX spectra of red giant and supergiant stars, that extend through optical and near-IR wavelengths. For the first time, models are also provided with modified surface abundances of C, N and O, as a step towards accounting for the changes known to occur due to convective dredge-up (and to be enhanced in the case of rotation). The aims are (i) to assess how well current models reproduce observed spectra, (ii) to quantify the effects of the abundance changes on the spectra, and (iii) to determine how these changes affect estimates of fundamental stellar parameters. Observed giant star spectra can be fitted very well at solar metallicity down to about 3400K. Modified surface abundances are preferred in only a minority of cases for luminosity class II, possibly indicating mixing in excess of standard first dredge-up. Supergiant stars show a larger variety of near-IR spectra, and good fits are currently obtained for about one third of the observations only. Modified surface abundances help reproducing strong CN bands, but do not suffice to resolve all the difficulties. The effect of the abundance changes on the estimated Teff depends on the wavelength range of observation and can amount several 100K. Reasons for the remaining discrepancies are discussed.

We report on recent progress in the modelling of the near-IR spectra of young stellar populations, i.e. populations in which red supergiants (RSGs) are dominant. First, we discuss the determination of fundamental parameters of RSGs using fits to their near-IR spectra with new PHOENIX model spectra; RSG-specific surface abundances are accounted for and effects of the microturbulence parameter are explored. New population synthesis predictions are then described and, as an example, it is shown that the spectra of young star clusters in M82 can be reproduced very well from 0.5 to 2.4 micrometers. We warn of remaining uncertainties in cluster ages.