Aims. Our goal is to study the different morphologies associated to the interaction of the stellar winds of AGB stars and red supergiants with the interstellar medium (ISM) to follow the fate of the circumstellar matter injected into the interstellar medium. Methods. Far-infrared Herschel/PACS images at 70 and 160 μm of a sample of 78 Galactic evolved stars are used to study the (dust) emission structures developing out of stellar wind-ISM interaction. In addition, two-fluid hydrodynamical simulations of the coupled gas and dust in wind-ISM interactions are used for comparison with the observations. Results. 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%. The total dust and gas mass inferred from the observed infrared emission is similar to the stellar mass loss over a period of a few thousand years, while in most cases it is less than the total ISM mass potentially swept-up by the wind-ISM interaction. 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. Conclusions. 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. In most cases the derived ISM densities concur with those typical of the warm neutral and ionised gas in the Galaxy, though some cases point towards the presence of cold diffuse clouds. Tentatively, the "eyes" class objects are associated to (visual) binaries, while the "rings" generally do not appear to occur for M-type stars, only for C or S-type objects that have experienced a thermal pulse.

Context. The influence of binarity on the late stages of stellar evolution remains an open issue. Aims. 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 exhibit a broad dust excess starting either at H or K, pointing to the presence of a gravitationally bound dusty disc in the system. We initiated an extensive multiwavelength study of those systems and here report on our radial velocity and photometric monitoring results for six stars of early F type, which are pulsators of small amplitude. Methods. To determine the radial velocity of low signal-to-noise ratio time-series data, we constructed dedicated autocorrelation masks based on high signal-to-noise ratio spectra, used in our published chemical studies. The radial velocity variations were analysed in detail to differentiate between pulsational variability and variability caused by orbital motion. When available, the photometric monitoring data were used to complement the time series of radial velocity data and to establish the nature of the pulsation. Finally, orbital minimalisation was performed to constrain the orbital elements. Results. 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 M and the companions are probably unevolved objects of (very) low initial mass. We argue that these binaries must have evolved through a phase of strong 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 have formed during this strong interaction phase. The eccentric orbits of these highly evolved objects remain poorly understood. In one object, the line-of-sight grazes the edge of the puffed-up inner rim of the disc. Conclusions. These results corroborate our earlier statement that evolved objects in binary stars create a Keplerian dusty circumbinary disc. 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. Tables 4-6 are only available in electronic form at

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 μm and SPIRE images at 250, 350, and 500 μm 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 from the central target and the presence of a linear bar at ∼9 . 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 ), 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 suggests a drastic change in mean gas and dust density ∼32 000 yr ago. Using hydrodynamical simulations, we try to explain the observed morphology of the bow shock around Betelgeuse. Conclusions. Different hypotheses, based on observational and theoretical constraints, are formulated to explain the origin of the multiple arcs and the linear bar and the fact that no large-scale instabilities are visible in the bow shock region. We infer that the two main ingredients for explaining these phenomena are a non-homogeneous mass-loss process and the influence of the Galactic magnetic field. The hydrodynamical simulations show that a warm interstellar medium, reflecting a warm neutral or partially ionized medium, or a higher temperature in the shocked wind also prevent the growth of strong instabilities. The linear bar is probably an interstellar structure illuminated by Betelgeuse itself.

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 ∼ 0.3), unusual for its short 3. d 65705 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 a simultaneous fitting of the light and the radial velocity curves. Individual star spectra were obtained by spectrum disentangling. The uncertainties of the fit were achieved 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 ∼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 the 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.

Context. It remains difficult to interpret the near-IR emission of young stellar populations. One main reason is our incomplete understanding of the spectra of luminous red stars. Aims. This work provides a grid of theoretical spectra of red giant and supergiant stars, that extends 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 that occur due to convective dredge-up in red supergiants or may occur at earlier evolutionary stages in the case of rotation. The aims are (i) to assess how well current models reproduce observed spectra, in particular in the near-IR; (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. Methods. Spectra are computed with the model atmosphere code PHOENIX and compared with a homogeneous set of observations. Although the empirical spectra have a resolution of only λ/∆λ ∼ 1000, we emphasize that models must be calculated at high spectral resolution in order to reproduce the shapes of line blends and molecular bands. Results. Giant star spectra of class III can be fitted extremely well at solar metallicity down to ∼3400 K, where difficulties appear in the modelling of near-IR H 2 O and TiO absorption bands. Luminous giants of class II can be fitted well too, with modified surface abundances preferred in a minority of cases, 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 the difficulties. The effect of the abundance changes on the estimated T eff depends on the wavelength range of observation and can amount several 100 K. Conclusions. While theoretical spectra for giant stars are becoming very satisfactory, red supergiants require further work. The model grid must be extended, in particular to larger micro-turbulent velocities. Some observed spectra may call for models with even lower gravities than explored here (and therefore probably stellar winds), and/or with more extreme abundances than predicted by standard non-rotating evolution models. Non-static atmospheres models should also be envisaged.

Context. Like other starburst galaxies, M 82 hosts compact, massive (>5 × 10 5 M ) young star clusters that are interesting both in their own right and as benchmarks for population synthesis models. Aims. In addition to assessing or reassessing the properties of some of the brightest near-IR sources in M 82, this paper addresses the following questions. Can population synthesis models at λ/δλ 750 adequately reproduce the near-IR spectral features and the energy distribution of these clusters between 0.8 and 2.4 μm? How do the derived cluster properties compare with previous results from optical studies? Methods. We analyse the spectra of 5 massive clusters in M 82, 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. Results. 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. Conclusions. The observed IR-bright clusters are part of the most recent episode of extended star formation in M 82. The near-IR study of clusters that are too faint for optical observation adds important elements to the age distribution of massive clusters in dusty starbursts. Further joint optical and near-IR spectroscopic studies will provide strong constraints on the uncertain physics of massive stars on which population synthesis models rest.

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 ± 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 ± 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 ± 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, ApJ, 459, L31), 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 ± 0.8) • 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 ± 14 pc. For a distance of 123 pc, we derive a mean luminosity of 7790 ± 150 L (internal error).

Herschel PACS and SPIRE images have been obtained over a 30 × 30 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 μm 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 −1 . 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 R 0 = (8.0 ± 0.3) × 10 17 cm. We also derive a relative velocity of the star with respect to the ISM of (106.6 ± 8.7)/ √ n ISM km s −1 , where n ISM is the number density of the local ISM.