The search for extrasolar planets has driven rapid advances in instrumentation, resulting in cameras such as SPHERE at the VLT, GPI at Gemini South and SCExAO at Subaru, capable of achieving very high contrast ($\sim10^{6}$) around bright stars with small inner working angles ($\sim 0\farcs{1}$). The optimal exploitation of data from these instruments depends on the availability of easy-to-use software to process and analyse their data products. We present a pure-python pipeline, {\sc precision}, which provides fast, memory-efficient reduction of data from the SPHERE/IRDIS near-infrared imager, and can be readily extended to other instruments. We apply {\sc precision} to observations of the extreme red supergiant VX~Sgr, the inner outflow of which is revealed to host complex, asymmetric structure in the near-IR. In addition, optical polarimetric imaging reveals clear extended polarised emission on $\sim0.5^{\prime\prime}$ scales which varies significantly with azimuth, confirming the asymmetry. While not conclusive, this could suggest that the ejecta are confined to a disc or torus, which we are viewing nearly face on, although other non-spherical or clumpy configurations remain possible. VX~Sgr has no known companions, making such a geometry difficult to explain, as there is no obvious source of angular momentum in the system.

We report on the succesful search for CO (2-1) and (3-2) emission associated with OH/IR stars in the Galactic Bulge. We observed a sample of eight extremely red AGB stars with the APEX telescope and detected seven. The sources were selected at sufficient high galactic latitude to avoid interference by interstellar CO, which hampered previous studies of inner galaxy stars. To study the nature of our sample and the mass loss we constructed the SEDs from photometric data and Spitzer IRS spectroscopy. In a first step we apply radiative transfer modelling to fit the SEDs and obtain luminosities and dust mass loss rates (MLR). Through dynamical modelling we then retrieve the total MLR and the gas-to-dust ratios. We derived variability periods of our stars. The luminosities range between approximately 4000 and 5500 Lsun and periods are below 700 days. The total MLR ranges between 1E-5 and 1E-4 Msun/yr. Comparison with evolutionary models shows that the progenitor mass is approximately 1.5 Msun, similar to the Bulge Miras and are of intermediate age (3 Gyr). The gas-to-dust ratios are between 100 and 400 and are similar to what is found for OH/IR stars in the galactic Disk. One star, IRAS 17347-2319, has a very short period of approximately 300 days which may be decreasing further. It may belong to a class of Mira variables with a sudden change in period as observed in some Galactic objects. It would be the first example of an OH/IR star in this class and deserves further follow-up observations.

We examine polycyclic aromatic hydrocarbon (PAH), dust and atomic/molecular emission toward the Galactic bulge using Spitzer Space Telescope observations of four fields: C32, C35, OGLE and NGC 6522. These fields are approximately centered on (l, b) = (0.0{\deg}, 1.0{\deg}), (0.0{\deg}, -1.0{\deg}), (0.4{\deg}, -2.1{\deg}) and (1.0{\deg}, -3.8{\deg}), respectively. Far-infrared photometric observations complement the Spitzer/IRS spectroscopic data and are used to construct spectral energy distributions. We find that the dust and PAH emission are exceptionally similar between C32 and C35 overall, in part explained due to their locations---they reside on or near boundaries of a 7 Myr-old Galactic outflow event and are partly shock-heated. Within the C32 and C35 fields, we identify a region of elevated H{\alpha} emission that is coincident with elevated fine-structure and [O IV] line emission and weak PAH feature strengths. We are likely tracing a transition zone of the outflow into the nascent environment. PAH abundances in these fields are slightly depressed relative to typical ISM values. In the OGLE and NGC 6522 fields, we observe weak features on a continuum dominated by zodiacal dust. SED fitting indicates that thermal dust grains in C32 and C35 have comparable temperatures to those of diffuse, high-latitude cirrus clouds. Little variability is detected in the PAH properties between C32 and C35, indicating that a stable population of PAHs dominates the overall spectral appearance. In fact, their PAH features are exceptionally similar to that of the M82 superwind, emphasizing that we are probing a local Galactic wind environment.

The radiation of stars heats dust grains in the diffuse interstellar medium and in star-forming regions in galaxies. Modelling this interaction provides information on dust in galaxies, a vital ingredient for their evolution. It is not straightforward to identify the stellar populations heating the dust, and to link attenuation to emission on a sub-galactic scale. Radiative transfer models are able to simulate this dust-starlight interaction in a realistic, three-dimensional setting. We investigate the dust heating mechanisms on a local and global galactic scale, using the Andromeda galaxy (M31) as our laboratory. We perform a series of panchromatic radiative transfer simulations of Andromeda with our code SKIRT. The high inclination angle of M31 complicates the 3D modelling and causes projection effects. However, the observed morphology and flux density are reproduced fairly well from UV to sub-millimeter wavelengths. Our model reveals a realistic attenuation curve, compatible with previous, observational estimates. We find that the dust in M31 is mainly (91 % of the absorbed luminosity) heated by the evolved stellar populations. The bright bulge produces a strong radiation field and induces non-local heating up to the main star-forming ring at 10 kpc. The relative contribution of unevolved stellar populations to the dust heating varies strongly with wavelength and with galactocentric distance.The dust heating fraction of unevolved stellar populations correlates strongly with NUV-r colour and specific star formation rate. These two related parameters are promising probes for the dust heating sources at a local scale.

The young star beta Pictoris is well known for its dusty debris disk, produced through the grinding down by collisions of planetesimals, kilometre-sized bodies in orbit around the star. In addition to dust, small amounts of gas are also known to orbit the star, likely the result from vaporisation of violently colliding dust grains. The disk is seen edge on and from previous absorption spectroscopy we know that the gas is very rich in carbon relative to other elements. The oxygen content has been more difficult to assess, however, with early estimates finding very little oxygen in the gas at a C/O ratio 20x higher than the cosmic value. A C/O ratio that high is difficult to explain and would have far-reaching consequences for planet formation. Here we report on observations by the far-infrared space telescope Herschel, using PACS, of emission lines from ionised carbon and neutral oxygen. The detected emission from C+ is consistent with that previously reported being observed by the HIFI instrument on Herschel, while the emission from O is hard to explain without assuming a higher-density region in the disk, perhaps in the shape of a clump or a dense torus, required to sufficiently excite the O atoms. A possible scenario is that the C/O gas is produced by the same process responsible for the CO clump recently observed by ALMA in the disk, and that the re-distribution of the gas takes longer than previously assumed. A more detailed estimate of the C/O ratio and the mass of O will have to await better constraints on the C/O gas spatial distribution.

Context. The recent detection of warm H$_2$O vapor emission from the outflows of carbon-rich asymptotic giant branch (AGB) stars challenges the current understanding of circumstellar chemistry. Two mechanisms have been invoked to explain warm H$_2$O vapor formation. In the first, periodic shocks passing through the medium immediately above the stellar surface lead to H$_2$O formation. In the second, penetration of ultraviolet interstellar radiation through a clumpy circumstellar medium leads to the formation of H$_2$O molecules in the intermediate wind. Aims. We aim to determine the properties of H$_2$O emission for a sample of 18 carbon-rich AGB stars and subsequently constrain which of the above mechanisms provides the most likely warm H$_2$O formation pathway. Methods, Results, and Conclusions. See paper.

Massive stars live short lives, losing large amounts of mass through their stellar wind. Their mass is a key factor determining how and when they explode as supernovae, enriching the interstellar medium with heavy elements and dust. During the red supergiant phase, mass-loss rates increase prodigiously, but the driving mechanism has proven elusive. Here we present high-contrast optical polarimetric-imaging observations of the extreme red supergiant VY Canis Majoris and its clumpy, dusty, mass-loss envelope, using the new extreme-adaptive-optics instrument SPHERE at the VLT. These observations allow us to make the first direct and unambiguous detection of submicron dust grains in the ejecta; we derive an average grain radius $\sim$ 0.5 $\mu$m, 50 times larger than in the diffuse ISM, large enough to receive significant radiation pressure by photon scattering. We find evidence for varying grain sizes throughout the ejecta, highlighting the dynamical nature of the envelope. Grains with 0.5 $\mu$m sizes are likely to reach a safe distance from the eventual explosion of VY Canis Majoris; hence it may inject upwards of 10$^{-2}$ M$_\odot$ of dust into the ISM.

We present far-infrared and submillimeter maps from the Herschel Space Observatory and the James Clerk Maxwell Telescope of the debris disk host star AU Microscopii. Disk emission is detected at 70, 160, 250, 350, 450, 500 and 850 micron. The disk is resolved at 70, 160 and 450 micron. In addition to the planetesimal belt, we detect thermal emission from AU Mic's halo for the first time. In contrast to the scattered light images, no asymmetries are evident in the disk. The fractional luminosity of the disk is $3.9 \times 10^{-4}$ and its mm-grain dust mass is 0.01 MEarth (+/- 20%). We create a simple spatial model that reconciles the disk SED as a blackbody of 53 +/- 2 K (a composite of 39 and 50 K components) and the presence of small (non-blackbody) grains which populate the extended halo. The best fit model is consistent with the "birth ring" model explored in earlier works, i.e., an edge-on dust belt extending from 8.8-40 AU, but with an additional halo component with an $r^{-1.5}$ surface density profile extending to the limits of sensitivity (140 AU). We confirm that AU Mic does not exert enough radiation force to blow out grains. For stellar mass loss rates of 10-100x solar, compact (zero porosity) grains can only be removed if they are very small, consistently with previous work, if the porosity is 0.9, then grains approaching 0.1 micron can be removed via corpuscular forces (i.e., the stellar wind).

ALMA observations show a non-detection of carbon monoxide around the four most luminous asymptotic giant branch (AGB) stars in the globular cluster 47 Tucanae. Stellar evolution models and star counts show that the mass-loss rates from these stars should be ~1.2-3.5 x 10^-7 solar masses per year. We would naively expect such stars to be detectable at this distance (4.5 kpc). By modelling the ultraviolet radiation field from post-AGB stars and white dwarfs in 47 Tuc, we conclude CO should be dissociated abnormally close to the stars. We estimate that the CO envelopes will be truncated at a few hundred stellar radii from their host stars and that the line intensities are about two orders of magnitude below our current detection limits. The truncation of CO envelopes should be important for AGB stars in dense clusters. Observing the CO (3-2) and higher transitions and targeting stars far from the centres of clusters should result in the detections needed to measure the outflow velocities from these stars.

In this article, we describe the MIRI Imager module (MIRIM), which provides broad-band imaging in the 5 - 27 microns wavelength range for the James Webb Space Telescope. The imager has a 0"11 pixel scale and a total unobstructed view of 74"x113". The remainder of its nominal 113"x113" field is occupied by the coronagraphs and the low resolution spectrometer. We present the instrument optical and mechanical design. We show that the test data, as measured during the test campaigns undertaken at CEA-Saclay, at the Rutherford Appleton Laboratory, and at the NASA Goddard Space Flight Center, indicate that the instrument complies with its design requirements and goals. We also discuss the operational requirements (multiple dithers and exposures) needed for optimal scientific utilization of the MIRIM.

An analysis of high-resolution near-infrared spectra of a sample of 45 asymptotic giant branch (AGB) stars towards the Galactic bulge is presented. The sample consists of two subsamples, a larger one in the inner and intermediate bulge, and a smaller one in the outer bulge. The data are analysed with the help of hydrostatic model atmospheres and spectral synthesis. We derive the radial velocity of all stars, and the atmospheric chemical mix ([Fe/H], C/O, $^{12}$C/$^{13}$C, Al, Si, Ti, and Y) where possible. Our ability to model the spectra is mainly limited by the (in)completeness of atomic and molecular line lists, at least for temperatures down to $T_{\rm eff}\approx3100$ K. We find that the subsample in the inner and intermediate bulge is quite homogeneous, with a slightly sub-solar mean metallicity and only few stars with super-solar metallicity, in agreement with previous studies of non-variable M-type giants in the bulge. All sample stars are oxygen-rich, C/O$<$1.0. The C/O and carbon isotopic ratios suggest that third dredge-up (3DUP) is absent among the sample stars, except for two stars in the outer bulge that are known to contain technetium. These stars are also more metal-poor than the stars in the intermediate or inner bulge. Current stellar masses are determined from linear pulsation models. The masses, metallicities and 3DUP behaviour are compared to AGB evolutionary models. We conclude that these models are partly in conflict with our observations. Furthermore, we conclude that the stars in the inner and intermediate bulge belong to a more metal-rich population that follows bar-like kinematics, whereas the stars in the outer bulge belong to the metal-poor, spheroidal bulge population.

Aim: The late stages of stellar evolution are mainly governed by the mass of the stars. Low- and intermediate-mass stars lose copious amounts of mass during the asymptotic giant branch (AGB) which obscure the central star making it difficult to study the stellar spectra and determine the stellar mass. In this study, we present observational data that can be used to determine lower limits to the stellar mass. Method: Spectra of nine heavily reddened AGB stars taken by the Herschel Space Observatory display numerous molecular emission lines. The strongest emission lines are due to H2O. We search for the presence of isotopologues of H2O in these objects. Result: We detected the 16O and 17O isotopologues of water in these stars, but lines due to H2^{18}O are absent. The lack of 18O is predicted by a scenario where the star has undergone hot-bottom burning which preferentially destroys 18O relative to 16O and 17O. From stellar evolution calculations, this process is thought to occur when the stellar mass is above 5 Msun for solar metallicity. Hence, observations of different isotopologues of H2O can be used to help determine the lower limit to the initial stellar mass. Conclusion: From our observations, we deduce that these extreme OH/IR stars are intermediate-mass stars with masses of >= 5 Msun. Their high mass-loss rates of ~ 1.0e-4 Msun/yr may affect the enrichment of the interstellar medium and the overall chemical evolution of our Galaxy.

Context: AGB stars lose a large percentage of their mass in a dust-driven wind. This creates a circumstellar envelope, which can be studied through thermal dust emission and molecular emission lines. In the case of high mass-loss rates, this study is complicated by the high optical depths and the intricate coupling between gas and dust radiative transfer characteristics. An important aspect of the physics of gas-dust interactions is the strong influence of dust on the excitation of several molecules, including H2O. Aims: The dust and gas content of the envelope surrounding the high mass-loss rate OH/IR star OH 127.8+0.0, as traced by Herschel observations, is studied, with a focus on the H2O content and the dust-to-gas ratio. We report detecting a large number of H2O vapor emission lines up to J = 9 in the Herschel data, for which we present the measured line strengths. Methods: The treatments of both gas and dust species are combined using two numerical radiative transfer codes. The method is illustrated for both low and high mass-loss-rate sources. Specifically, we discuss different ways of assessing the dust-to-gas ratio: 1) from the dust thermal emission spectrum and the CO molecular gas line strengths; 2) from the momentum transfer from dust to gas and the measured gas terminal velocity; and 3) from the determination of the required amount of dust to reproduce H2O lines for a given H2O vapor abundance. These three diagnostics probe different zones of the outflow, for the first time allowing an investigation of a possible radial dependence of the dust-to-gas ratio. Results: ... Continued in article.

We study the grain properties and location of the forsterite crystals in the circumstellar environment of the pre-planetary nebula (PPN) IRAS 17150-3224 in order to learn more about the as yet poorly understood evolutionary phase prior to the PPN. We use the best-fit model for IRAS 17150-3224 of Meixner et al. (2002) and add forsterite to this model. We investigate different spatial distributions and grain sizes of the forsterite crystals in the circumstellar environment. We compare the spectral bands of forsterite in the mid-infrared and at 69 micrometre in radiative transport models to those in ISO-SWS and Herschel/PACS observations. We can reproduce the non-detection of the mid-infrared bands and the detection of the 69 micrometre feature with models where the forsterite is distributed in the whole outflow, in the superwind region, or in the AGB-wind region emitted previous to the superwind, but we cannot discriminate between these three models. To reproduce the observed spectral bands with these three models, the forsterite crystals need to be dominated by a grain size population of 2 micrometre up to 6 micrometre. We hypothesise that the large forsterite crystals were formed after the superwind phase of IRAS 17150-3224, where the star developed an as yet unknown hyperwind with an extremely high mass-loss rate (10^-3 Msol/yr). The high densities of such a hyperwind could be responsible for the efficient grain growth of both amorphous and crystalline dust in the outflow. Several mechanisms are discussed that might explain the lower-limit of 2 micrometre found for the forsterite grains, but none are satisfactory. Among the mechanisms explored is a possible selection effect due to radiation pressure based on photon scattering on micron-sized grains.

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

The evolution of low- and intermediate-mass stars on the asymptotic giant branch (AGB) is mainly controlled by the rate at which these stars lose mass in a stellar wind. Understanding the driving mechanism and strength of the stellar winds of AGB stars and the processes enriching their surfaces with products of nucleosynthesis are paramount to constraining AGB evolution and predicting the chemical evolution of galaxies. In a previous paper we have constrained the structure of the outflowing envelope of W Hya using spectral lines of the $^{12}$CO molecule. Here we broaden this study by modelling an extensive set of H$_{2}$O and $^{28}$SiO lines observed by the three instruments on board Herschel using a state-of-the-art molecular excitation and radiative transfer code. The oxygen isotopic ratios and the $^{28}$SiO abundance profile can be connected to the initial stellar mass and to crucial aspects of dust formation at the base of the stellar wind, respectively. The modelling of H$_{2}$O and $^{28}$SiO confirms the properties of the envelope model of W Hya derived from $^{12}$CO lines. We find an H$_2$O ortho-to-para ratio of 2.5\,$^{+2.5}_{-1.0}$, consistent with what is expected for an AGB wind. The O$^{16}$/O$^{17}$ ratio indicates that W Hya has an initial mass of about 1.5 M$_\odot$. Although the ortho- and para-H$_{2}$O lines observed by HIFI appear to trace gas of slightly different physical properties, a turbulence velocity of $0.7\pm0.1$ km s$^{-1}$ fits the HIFI lines of both spin isomers and those of $^{28}$SiO well. The ortho- and para-H$_2^{16}$O and $^{28}$SiO abundances relative to H$_{2}$ are $(6^{+3}_{-2}) \times 10^{-4}$, $(3^{+2}_{-1}) \times 10^{-4}$, and $(3.3\pm 0.8)\times 10^{-5}$, respectively. Assuming a solar silicon-to-carbon ratio, the $^{28}$SiO line emission model is consistent with about one-third of the silicon atoms being locked up in dust particles.

$\tau$ Ceti is a nearby, mature G-type star very similar to our Sun, with a massive Kuiper Belt analogue (Greaves et al. 2004) and possible multiplanet system (Tuomi et al. 2013) that has been compared to our Solar System. We present Herschel Space Observatory images of the debris disk, finding the disk is resolved at 70 and 160 microns, and marginally resolved at 250 microns. The Herschel images and infrared photometry from the literature are best modelled using a wide dust annulus with an inner edge between 1-10 AU and an outer edge at ~55 AU, inclined from face-on by 35$\pm$10 degrees, and with no significant azimuthal structure. We model the proposed tightly-packed planetary system of five super-Earths and find that the innermost dynamically stable disk orbits are consistent with the inner edge found by the observations. The photometric modelling, however, cannot rule out a disk inner edge as close to the star as 1 AU, though larger distances produce a better fit to the data. Dynamical modelling shows that the 5 planet system is stable with the addition of a Neptune or smaller mass planet on an orbit outside 5 AU, where the Tuomi et al. analysis would not have detected a planet of this mass.

The imprints of dust-starlight interactions are visible in scaling relations between stellar characteristics, star formation parameters and dust properties. We aim to examine dust scaling relations on a sub-kpc resolution in the Andromeda galaxy (M31) by comparing the properties on a local and global scale to other galaxies of the local universe. New Herschel observations are combined with available data from GALEX, SDSS, WISE and Spitzer to construct a dataset covering UV to submm wavelengths. We work at the resolution of the SPIRE $500\; \mu$m beam, with pixels corresponding to physical regions of 137 x 608 pc in the galaxy's disk. A panchromatic spectral energy distribution was modelled for each pixel and several dust scaling relations are investigated. We find, on a sub-kpc scale, strong correlations between $M_d/M_\star$ and NUV-r, and between $M_d/M_\star$ and $\mu_\star$ (the stellar mass surface density). Striking similarities with corresponding relations based on integrated galaxies are found. We decompose M31 in four macro-regions based on their FIR morphology; the bulge, inner disk, star forming ring and the outer disk. All regions closely follow the galaxy-scale average trends. The specific star formation characteristics we derive for these macro-regions give strong hints of an inside-out formation of the bulge-disk geometry, as well as an internal downsizing process. However, within each macro-region, a great diversity in individual micro-regions is found. Furthermore, we confirm that dust in the bulge of M31 is heated only by the old stellar populations. In general, the local dust scaling relations indicate that the dust content in M31 is maintained by a subtle interplay of past and present star formation. The similarity with galaxy-based relations strongly suggests that they are in situ correlations, with underlying processes that must be local in nature. (Abriged)

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.

We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67-72 $\mu$m. This wavelength range covers the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{2-2x}\text{Fe}_{(2x)}\text{SiO}_{4}$). The width and wavelength position of this band are sensitive to the temperature and composition of the crystalline olivine. Our sample covers a wide range of objects: from high mass-loss rate AGB stars (OH/IR stars, $\dot M \ge 10^{-5}$ M$_\odot$/yr), through post-AGB stars with and without circumbinary disks, to planetary nebulae and even a few massive evolved stars. The goal of this study is to exploit the spectral properties of the 69 $\mu$m band to determine the composition and temperature of the crystalline olivine. Since the objects cover a range of evolutionary phases, we study the physical and chemical properties in this range of physical environments. We fit the 69 $\mu$m band and use its width and position to probe the composition and temperature of the crystalline olivine. For 27 sources in the sample, we detected the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{(2-2x)}\text{Fe}_{(2x)}\text{SiO}_{4}$). The 69 $\mu$m band shows that all the sources produce pure forsterite grains containing no iron in their lattice structure. The temperature of the crystalline olivine as indicated by the 69 $\mu$m band, shows that on average the temperature of the crystalline olivine is highest in the group of OH/IR stars and the post-AGB stars with confirmed Keplerian disks. The temperature is lower for the other post-AGB stars and lowest for the planetary nebulae. A couple of the detected 69 $\mu$m bands are broader than those of pure magnesium-rich crystalline olivine, which we show can be due to a temperature gradient in the circumstellar environment of these stars. continued...

Context: The dusty debris disk around the $\sim$20 Myr old main-sequence A-star {\beta} Pictoris is known to contain gas. Evidence points towards a secondary origin of the gas as opposed to being a direct remnant form the initial protoplanetary disk, although the dominant gas production mechanism is so far not identified. The origin of the observed overabundance of C and O compared to solar abundances of metallic elements, e.g. Na and Fe, is also unclear. Aims: Our goal is to constrain the spatial distribution of C in the disk, and thereby the gas origin and its abundance pattern. Methods: We used the HIFI instrument onboard Herschel to observe and spectrally resolve CII emission at 158 $\mu$m from the {\beta} Pic debris disk. Assuming Keplerian rotation, we use the spectrally resolved line profile to constrain the spatial distribution of the gas. Results: We show that most of the gas is located around $\sim$100 AU or beyond. We estimate a total C gas mass of $1.3\times10^{-2}$ M$_\oplus$. The data suggest that more gas is located on the southwest side of the disk than on the northeast side. The data are consistent with the hypothesis of a well-mixed gas (constant C/Fe ratio throughout the disk). Assuming instead a spatial profile expected from a simplified accretion disk model, we found it to give a significantly worse fit to the observations. Conclusions: Since the bulk of the gas is found outside 30 AU, we argue that the cometary objects known as "falling evaporating bodies" are unlikely to be the dominant source of gas; production from grain-grain collisions or photodesorption seems more likely. The incompatibility of the observations with a simplified accretion disk model could favour a preferential depletion explanation for the overabundance of C and O. More stringent constraints on the spatial distribution will be available from ALMA observations of CI at 609 $\mu$m.

Asymptotic giant branch (AGB) stars lose their envelopes by means of a stellar wind whose driving mechanism is not understood well. Characterizing the composition and thermal and dynamical structure of the outflow provides constraints that are essential for understanding AGB evolution, including the rate of mass loss and isotopic ratios. We modeled the CO emission from the wind of the low mass-loss rate oxygen-rich AGB star W Hya using data obtained by the HIFI, PACS, and SPIRE instruments onboard the Herschel Space Observatory and ground-based telescopes. $^{12}$CO and $^{13}$CO lines are used to constrain the intrinsic $^{12}$C/$^{13}$C ratio from resolved HIFI lines. The acceleration of the outflow up to about 5.5 km/s is quite slow and can be represented by a beta-type velocity law with index 5. Beyond this point, acceleration up the terminal velocity of 7 km/s is faster. Using the J=10-9, 9-8, and 6-5 transitions, we find an intrinsic $^{12}$C/$^{13}$C ratio of $18\pm10$ for W Hya, where the error bar is mostly due to uncertainties in the $^{12}$CO abundance and the stellar flux around 4.6 $\mu$m. To match the low-excitation CO lines, these molecules need to be photo-dissociated at about 500 stellar radii. The radial dust emission intensity profile measured by PACS images at 70 $\mu$m shows substantially stronger emission than our model predicts beyond 20 arcsec. The initial slow acceleration of the wind implies inefficient wind driving in the lower part of the envelope. The final injection of momentum in the wind might be the result of an increase in the opacity thanks to the late condensation of dust species. The derived intrinsic isotopologue ratio for W Hya is consistent with values set by the first dredge-up and suggestive of an initial mass of 2 M$_\odot$ or more. However, the uncertainty in the main-sequence mass derived based on this isotopologic ratio is large.

Spectra of OH/IR stars show prominent spectral bands of crystalline olivine (Mg$_{(2-2x)}$Fe$_{(2x)}$SiO$_{4}$). To learn more about the timescale of the outflows of OH/IR stars, we study the spectral band of crystalline olivine at 69 {\mu}m. The 69 {\mu}m band is of interest because its width and peak wavelength position are sensitive to the grain temperature and to the exact composition of the crystalline olivine. With Herschel/PACS, we observed the 69 {\mu}m band in the outflow of 14 OH/IR stars. By comparing the crystalline olivine features of our sample with those of model spectra, we determined the size of the outflow and its crystalline olivine abundance. The temperature indicated by the observed 69 {\mu}m bands can only be reproduced by models with a geometrically compact superwind ($R_{\rm{SW}}\lesssim$ 2500 AU = 1400 R$_{*}$).This means that the superwind started less than 1200 years ago (assuming an outflow velocity of 10 km/s). The small amount of mass lost in one superwind and the high progenitor mass of the OH/IR stars introduce a mass loss and thus evolutionary problem for these objects, which has not yet been understood.

We present an analysis of the far-infrared and submillimetre molecular emission line spectrum of the luminous M-supergiant VY CMa, observed with the SPIRE and PACS spectrometers aboard the Herschel Space Observatory. Over 260 emission lines were detected in the 190-650-micron SPIRE FTS spectra, with one-third of the observed lines being attributable to H2O. Other detected species include CO, 13CO, H2^18O, SiO, HCN, SO, SO2, CS, H2S, and NH3. Our model fits to the observed 12CO and 13CO line intensities yield a 12C/13C ratio of 5.6+-1.8, consistent with measurements of this ratio for other M supergiants, but significantly lower than previously estimated for VY CMa from observations of lower-J lines. The spectral line energy distribution for twenty SiO rotational lines shows two temperature components: a hot component at 1000 K, which we attribute to the stellar atmosphere and inner wind, plus a cooler ~200 K component, which we attribute to an origin in the outer circumstellar envelope. We fit the line fluxes of 12CO, 13CO, H2O and SiO, using the SMMOL non-LTE line transfer code, with a mass-loss rate of 1.85x10^-4 Msun yr^-1 between 9 R* and 350 R*. To fit the observed line fluxes of 12CO, 13CO, H2O and SiO with SMMOL non-LTE line radiative transfer code, along with a mass-loss rate of 1.85x10^-4 Msun yr^-1. To fit the high rotational lines of CO and H2O, the model required a rather flat temperature distribution inside the dust condensation radius, attributed to the high H2O opacity. Beyond the dust condensation radius the gas temperature is fitted best by an r^-0.5 radial dependence, consistent with the coolant lines becoming optically thin. Our H2O emission line fits are consistent with an ortho:para ratio of 3 in the outflow.

This paper provides an overview of the PACS photometer flux calibration concept, in particular for the principal observation mode, the scan map. The absolute flux calibration is tied to the photospheric models of five fiducial stellar standards (alpha Boo, alpha Cet, alpha Tau, beta And, gamma Dra). The data processing steps to arrive at a consistent and homogeneous calibration are outlined. In the current state the relative photometric accuracy is around 2% in all bands. Starting from the present calibration status, the characterization and correction for instrumental effects affecting the relative calibration accuracy is described and an outlook for the final achievable calibration numbers is given. After including all the correction for the instrumental effects, the relative photometric calibration accuracy (repeatability) will be as good as 0.5% in the blue and green band and 2% in the red band. This excellent calibration starts to reveal possible inconsistencies between the models of the K-type and the M-type stellar calibrators. The absolute calibration accuracy is therefore mainly limited by the 5% uncertainty of the celestial standard models in all three bands. The PACS bolometer response was extremely stable over the entire Herschel mission and a single, time-independent response calibration file is sufficient for the processing and calibration of the science observations. The dedicated measurements of the internal calibration sources were needed only to characterize secondary effects. No aging effects of the bolometer or the filters have been found. Also, we found no signs of filter leaks. The PACS photometric system is very well characterized with a constant energy spectrum nu*Fnu = lambda*Flambda = const as a reference. Colour corrections for a wide range of sources SEDs are determined and tabulated.

We report spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1'' of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations.

Context. Edge-on spiral galaxies with prominent dust lanes provide us with an excellent opportunity to study the distribution and properties of the dust within them. The HEROES project was set up to observe a sample of seven large edge-on galaxies across various wavelengths for this investigation. Aims. Within this first paper, we present the Herschel observations and perform a qualitative and quantitative analysis on them, and we derive some global properties of the far infrared and submillimetre emission. Methods. We determine horizontal and vertical profiles from the Herschel observations of the galaxies in the sample and describe the morphology. Modified black-body fits to the global fluxes, measured using aperture photometry, result in dust temperatures and dust masses. The latter values are compared to those that are derived from radiative transfer models taken from the literature. Results. On the whole, our Herschel flux measurements agree well with archival values. We find that the exponential horizontal dust distribution model often used in the literature generally provides a good description of the observed horizontal profiles. Three out of the seven galaxies show signatures of extended vertical emission at 100 and 160 {\mu}m at the 5{\sigma} level, but in two of these it is probably due to deviations from an exactly edge-on orientation. Only for NGC 4013, a galaxy in which vertically extended dust has already been detected in optical images, we can detect vertically extended dust, and the derived scaleheight agrees with the value estimated through radiative transfer modelling. Our analysis hints at a correlation between the dust scaleheight and its degree of clumpiness, which we infer from the difference between the dust masses as calculated from modelling of optical data and from fitting the spectral energy distribution of Herschel datapoints.

The Mass loss of Evolved StarS (MESS) sample offers a selection of 78 Asymptotic Giant Branch (AGB) stars and Red Supergiants (RSGs) observed with the PACS photometer on-board Herschel at 70 and 160 {\mu}m. For most of these objects, the dusty AGB wind differs from spherical symmetry and the wind shape can be subdivided into four classes. In the present paper we concentrate on the influence of a companion on the morphology of the stellar wind. Literature was searched to find binaries in the MESS sample and these are subsequently linked to their wind-morphology class to assert that the binaries are not distributed equally among the classes. In the second part of the paper we concentrate on the circumstellar environment of the two prominent objects R Aqr and W Aql. Each shows a characteristic signature of a companion interaction with the stellar wind. For the symbiotic star R Aqr, PACS revealed two perfectly opposing arms which in part reflect the previously observed ring-shaped nebula in the optical. However, from the far-IR there is evidence that the emitting region is elliptical rather than circular. The outline of the wind of W Aql seems to follow a large Archimedean spiral formed by the orbit of the companion but also shows strong indications of an interaction with the interstellar medium. The nature of the companion of W Aql was investigated and the magnitude of the orbital period supports the size of the spiral outline.

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]

Some planetary systems harbour debris disks containing planetesimals such as asteroids and comets. Collisions between such bodies produce small dust particles, the spectral features of which reveal their composition and, hence, that of their parent bodies. A measurement of the composition of olivine crystals has been done for the protoplanetary disk HD100546 and for olivine crystals in the warm inner parts of planetary systems. The latter compares well with the iron-rich olivine in asteroids (x<0.29). In the cold outskirts of the Beta Pictoris system, an analogue to the young Solar System, olivine crystals were detected but their composition remained undetermined, leaving unknown how the composition of the bulk of Solar System cometary olivine grains compares with that of extrasolar comets. Here we report the detection of the 69-micrometre-wavelength band of olivine crystals in the spectrum of Beta Pictoris. Because the disk is optically thin, we can associate the crystals with an extrasolar proto-Kuiper belt a distance of 15-45 astronomical units from the star (one astronomical unit is the Sun-Earth distance), determine their magnesium-rich composition (x=0.01+/-0.001) and show that they make up 3.6+/-1.0 per cent of the total dust mass. These values are strikingly similar to those for the dust emitted by the most primitive comets in the Solar System, even though Beta Pictoris is more massive and more luminous and has a different planetary system architecture.

We present new Herschel/PACS images at 70, 100, and 160 micron of the well-known, nearby, carbon-rich asymptotic giant branch star IRC+10216 revealing multiple dust shells in its circumstellar envelope. For the first time, dust shells (or arcs) are detected until 320 arcsec. The almost spherical shells are non-concentric and have an angular extent between 40 deg and 200 deg. The shells have a typical width of 5 arcsec - 8 arcsec, and the shell separation varies in the range of 10 arcsec - 35 arcsec, corresponding to 500-1700 yr. Local density variations within one arc are visible. The shell/intershell density contrast is typically 4, and the arcs contain some 50% more dust mass than the smooth envelope. The observed (nested) arcs record the mass-loss history over the past 16 000 yr, but Rayleigh-Taylor and Kelvin-Helmholtz instabilities in the turbulent astropause and astrosheath will erase any signature of the mass-loss history for at least the first 200 000 yr of mass loss. Accounting for the bowshock structure, the envelope mass around IRC+10216 contains >2Msun of gas and dust mass. It is argued that the origin of the shells is related to non-isotropic mass-loss events and clumpy dust formation.

We have obtained Herschel images at five wavelengths from 100 to 500 micron of a ~5.5x2.5 degree area centred on the local galaxy M31 (Andromeda), our nearest neighbour spiral galaxy, as part of the Herschel guaranteed time project "HELGA". The main goals of HELGA are to study the characteristics of the extended dust emission, focusing on larger scales than studied in previous observations of Andromeda at an increased spatial resolution, and the obscured star formation. In this paper we present data reduction and Herschel maps, and provide a description of the far-infrared morphology, comparing it with features seen at other wavelengths. We use high--resolution maps of the atomic hydrogen, fully covering our fields, to identify dust emission features that can be associated to M31 with confidence, distinguishing them from emission coming from the foreground Galactic cirrus. Thanks to the very large extension of our maps we detect, for the first time at far-infrared wavelengths, three arc-like structures extending out to ~21, ~26 and ~31 kpc respectively, in the south-western part of M31. The presence of these features, hosting ~2.2e6 Msol of dust, is safely confirmed by their detection in HI maps. Overall, we estimate a total dust mass of ~5.8e7 Msol, about 78% of which is contained in the two main ring-like structures at 10 and 15 kpc, at an average temperature of 16.5 K. We find that the gas-to-dust ratio declines exponentially as a function of the galacto-centric distance, in agreement with the known metallicity gradient, with values ranging from 66 in the nucleus to ~275 in the outermost region. [Abridged]

We present an analysis of the dust and gas in Andromeda, using Herschel images sampling the entire far-infrared peak. We fit a modified-blackbody model to ~4000 quasi-independent pixels with spatial resolution of ~140pc and find that a variable dust-emissivity index (beta) is required to fit the data. We find no significant long-wavelength excess above this model suggesting there is no cold dust component. We show that the gas-to-dust ratio varies radially, increasing from ~20 in the center to ~70 in the star-forming ring at 10kpc, consistent with the metallicity gradient. In the 10kpc ring the average beta is ~1.9, in good agreement with values determined for the Milky Way (MW). However, in contrast to the MW, we find significant radial variations in beta, which increases from 1.9 at 10kpc to ~2.5 at a radius of 3.1kpc and then decreases to 1.7 in the center. The dust temperature is fairly constant in the 10kpc ring (ranging from 17-20K), but increases strongly in the bulge to ~30K. Within 3.1kpc we find the dust temperature is highly correlated with the 3.6 micron flux, suggesting the general stellar population in the bulge is the dominant source of dust heating there. At larger radii, there is a weak correlation between the star formation rate and dust temperature. We find no evidence for 'dark gas' in M31 in contrast to recent results for the MW. Finally, we obtained an estimate of the CO X-factor by minimising the dispersion in the gas-to-dust ratio, obtaining a value of (1.9+/-0.4)x10^20 cm^-2 [K kms^-1]^-1.

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 report on the discovery of an infrared emission band in the Spitzer spectrum of the S-type AGB star NP Aurigae that is caused by TiO molecules in the circumstellar environment. We modelled the observed emission to derive the temperature of the TiO molecules (\approx 600 K), an upper limit on the column density (\approx 10^17.25 cm^{-2}) and a lower limit on the spatial extent of the layer that contains these molecules. (\approx 4.6 stellar radii). This is the first time that this TiO emission band is observed. A search for similar emission features in the sample of S-type stars yielded two additional candidates. However, owing to the additional dust emission, the identification is less stringent. By comparing the stellar characteristics of NP Aur to those of the other stars in our sample, we find that all stars with TiO emission show large-amplitude pulsations, s-process enrichment, and a low C/O ratio. These characteristics might be necessary requirements for a star to show TiO in emission, but they are not sufficient.

Fomalhaut is a young, nearby star that is suspected to harbor an infant planetary system, interspersed with one or more belts of dusty debris. We present far-infrared images obtained with the Herschel Space Observatory with an angular resolution between 5.7 and 36.7 arcsec at wavelengths between 70 and 500 micrometer. The images show the main debris belt in great detail. Even at high spatial resolution, the belt appears smooth. The region in between the belt and the central star is not devoid of material; thermal emission is observed here as well. Also at the location of the star, excess emission is detected. We use a dynamical model together with radiative-transfer tools to derive the parameters of the debris disk. We include detailed models of the interaction of the dust grains with radiation, for both the radiation pressure and the temperature determination. Comparing these models to the spatially resolved temperature information contained in the images allows us to place strong constraints on the presence of grains that will be blown out of the system by radiation pressure. We use this to derive the dynamical parameters of the system. The appearance of the belt points towards a remarkably active system in which dust grains are produced at a very high rate by a collisional cascade in a narrow region filled with dynamically excited planetesimals. Dust particles with sizes below the blow-out size are abundantly present. The equivalent of 2000 one-km-sized comets are destroyed every day, out of a cometary reservoir amounting to 110 Earth masses. From comparison of their scattering and thermal properties, we find evidence that the dust grains are fluffy aggregates, which indicates a cometary origin. The excess emission at the location of the star may be produced by hot dust with a range of temperatures, but may also be due to gaseous free-free emission from a stellar wind.

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.

We obtained high-resolution near-IR spectra of 45 AGB stars located in the Galactic bulge. The aim of the project is to determine key elemental abundances in these stars to help constrain the formation history of the bulge. A further aim is to link the photospheric abundances to the dust species found in the winds of the stars. Here we present a progress report of the analysis of the spectra.

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.

Herschel's PACS instrument observed the environment of the binary system Mira Ceti in the 70 and 160 micron bands. These images reveal bright structures shaped as five broken arcs and fainter filaments in the ejected material of Mira's primary star. The overall shape of the IR emission around Mira deviates significantly from the expected alignment with Mira's exceptionally high space velocity. The observed broken arcs are neither connected to each other nor are they of a circular shape; they stretch over angular ranges of 80 to 100 degrees. By comparing Herschel and GALEX data, we found evidence for the disruption of the IR arcs by the fast outflow visible in both Halpha and the far UV. Radial intensity profiles are derived, which place the arcs at distances of 6-85" (550 - 8000 AU) from the binary. Mira's IR environment appears to be shaped by the complex interaction of Mira's wind with its companion, the bipolar jet, and the ISM.

The asymptotic giant branch (AGB) stars X Her and TX Psc have been imaged at 70 and 160 microns with the PACS instrument onboard the Herschel satellite, as part of the large MESS (Mass loss of Evolved StarS) Guaranteed Time Key Program. The images reveal an axisymmetric extended structure with its axis oriented along the space motion of the stars. This extended structure is very likely to be shaped by the interaction of the wind ejected by the AGB star with the surrounding interstellar medium (ISM). As predicted by numerical simulations, the detailed structure of the wind-ISM interface depends upon the relative velocity between star+wind and the ISM, which is large for these two stars (108 and 55 km/s for X Her and TX Psc, respectively). In both cases, there is a compact blob upstream whose origin is not fully elucidated, but that could be the signature of some instability in the wind-ISM shock. Deconvolved images of X Her and TX Psc reveal several discrete structures along the outermost filaments, which could be Kelvin-Helmholtz vortices. Finally, TX Psc is surrounded by an almost circular ring (the signature of the termination shock?) that contrasts with the outer, more structured filaments. A similar inner circular structure seems to be present in X Her as well, albeit less clearly.

In 2001, the discovery of circumstellar water vapour around the ageing carbon star IRC+10216 was announced. This detection challenged the current understanding of chemistry in old stars, since water vapour was predicted to be absent in carbon-rich stars. Several explanations for the occurrence of water vapour were postulated, including the vaporization of icy bodies (comets or dwarf planets) in orbit around the star, grain surface reactions, and photochemistry in the outer circumstellar envelope. However, the only water line detected so far from one carbon-rich evolved star can not discriminate, by itself, between the different mechanisms proposed. Here we report on the detection by the Herschel satellite of dozens of water vapour lines in the far-infrared and sub-millimetre spectrum of IRC+10216, including some high-excitation lines with energies corresponding to ~1000 K. The emission of these high-excitation water lines can only be explained if water vapour is present in the warm inner region of the envelope. A plausible explanation for the formation of warm water vapour appears to be the penetration of ultraviolet (UV) photons deep into a clumpy circumstellar envelope. This mechanism triggers also the formation of other molecules such as ammonia, whose observed abundances are much higher than hitherto predicted.

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.

Observations of high-excitation molecular emission lines can greatly increase our understanding of AGB winds, as they trace the innermost regions of the circumstellar envelope. The PACS spectrometer on-board the Herschel Space Telescope, provides for the first time the spectral resolution and sensitivity necessary to trace these lines. We report on the first modelling efforts of a PACS spectral scan for the OH/IR star V669 Cas. Central to our methodology is the consistent treatment of both dust and gas by using a line radiative transfer and a continuum radiative transfer code conjointly. Water emission lines are found to be extremely sensitive to the dust-to-gas ratio, emphasizing the need of consistent modelling for dust and gas.

In this article we present the detection of the 69 {\mu}m band of the crystalline olivine forsterite within the MESS key program of Herschel. We determine the temperature of the forsterite grains by fitting the 69 {\mu}m band.

We present a diagnostic tool to determine the forsterite abundance of the dust ejected by AGB stars. Our method is based on a comparison between the observed strength of spectral bands of forsterite and model calculations. We show that the 11.3 {\mu}m forsterite band is a robust indicator of the forsterite abundance of the current mass-loss period for AGB stars with an optically thick dust shell. The 33.6 {\mu}m band of forsterite is sensitive to changes in the density and the geometry of the emitting dust shell, and so a less robust indicator. We apply this method to six high mass-loss rate AGB stars, showing that AGB stars can have forsterite abundances of 12% by mass and higher, which is more than the previously found maximum abundance of 5%.

METIS is a mid-infrared instrument proposed for the European Extremely Large Telescope (E-ELT). It is designed to provide imaging and spectroscopic capabilities in the 3 - 14 micron region up to a spectral resolution of 100000. One of the novel concepts of METIS is that of a high-resolution integral field spectrograph (IFS) for a diffraction-limited mid-IR instrument. While this concept has many scientific and operational advantages over a long-slit spectrograph, one drawback is that the spectral resolution changes over the field of view. This has an impact on the procedures to correct for telluric absorption lines imprinted on the science spectra. They are a major obstacle in the quest to maximize spectral fidelity, the ability to distinguish a weak spectral feature from the continuum. The classical technique of division by a standard star spectrum, observed in a single IFS spaxel, cannot simply be applied to all spaxels, because the spectral resolution changes from spaxel to spaxel. Here we present and discuss possible techniques of telluric line correction of METIS IFS spectra, including the application of synthetic model spectra of telluric transmission, to maximize spectral fidelity.

In this paper we will discuss the images of Planetary Nebulae that have recently been obtained with PACS and SPIRE on board the Herschel satellite. This comprises results for NGC 650 (the little Dumbbell nebula), NGC 6853 (the Dumbbell nebula), and NGC 7293 (the Helix nebula).

We report on the initial analysis of Herschel/HIFI carbon monoxide (CO) observations of the Martian atmosphere performed between 11 and 16 April 2010. We selected the (7-6) rotational transitions of the isotopes ^{13}CO at 771 GHz and C^{18}O at 768 GHz in order to retrieve the mean vertical profile of temperature and the mean volume mixing ratio of carbon monoxide. The derived temperature profile agrees within less than 5 K with general circulation model (GCM) predictions up to an altitude of 45 km, however, show about 12-15 K lower values at 60 km. The CO mixing ratio was determined as 980 \pm 150 ppm, in agreement with the 900 ppm derived from Herschel/SPIRE observations in November 2009.

We report on an initial analysis of Herschel/HIFI observations of hydrogen chloride (HCl), hydrogen peroxide (H_2O_2), and molecular oxygen (O_2) in the Martian atmosphere performed on 13 and 16 April 2010 (L_s ~ 77{\deg}). We derived a constant volume mixing ratio of 1400 +/- 120 ppm for O_2 and determined upper limits of 200 ppt for HCl and 2 ppb for H_2O_2. Radiative transfer model calculations indicate that the vertical profile of O_2 may not be constant. Photochemical models determine the lowest values of H_2O_2 to be around L_s ~ 75{\deg} but overestimate the volume mixing ratio compared to our measurements.

The high spectral resolution and sensitivity of Herschel/HIFI allows for the detection of multiple rotational water lines and accurate determinations of water production rates in comets. In this letter we present HIFI observations of the fundamental 110-101 (557 GHz) ortho and 111-000 (1113 GHz) para rotational transitions of water in comet 81P/Wild 2 acquired in February 2010. We mapped the extent of the water line emission with five point scans. Line profiles are computed using excitation models which include excitation by collisions with electrons and neutrals and solar infrared radiation. We derive a mean water production rate of $1.0 \times 10^{28}$ molecules s$^{-1}$ at a heliocentric distance of 1.61 AU about 20 days before perihelion, in agreement with production rates measured from the ground using observations of the 18-cm OH lines. Furthermore, we constrain the electron density profile and gas kinetic temperature, and estimate the coma expansion velocity by fitting the water line shapes.

We present results of performance modelling for METIS, the Mid-infrared European Extremely Large Telescope (E-ELT) Imager and Spectrograph. Designed by a consortium of NOVA (Netherlands), UK Astronomy Technology Centre (UK), MPIA Heidelberg (Germany), CEA Saclay (France) and KU Leuven (Belgium), METIS will cover the atmospheric windows in L, M and N-band and will offer imaging, medium-resolution slit spectroscopy (R~1000-3000) and high-resolution integral field spectroscopy (R~100,000). Our model uses a detailed set of input parameters for site characteristics and atmospheric profiles, optical design, thermal background and the most up-to-date IR detector specifications. We show that METIS will bring an orders-of-magnitude level improvement in sensitivity and resolution over current ground-based IR facilities, bringing mid-IR sensitivities to the micro-Jansky regime. As the only proposed E-ELT instrument to cover this entire spectral region, and the only mid-IR high-resolution integral field unit planned on the ground or in space, METIS will open up a huge discovery space in IR astronomy in the next decade.

We report on the initial analysis of a Herschel/PACS full range spectrum of Neptune, covering the 51-220 micrometer range with a mean resolving power of ~ 3000, and complemented by a dedicated observation of CH4 at 120 micrometers. Numerous spectral features due to HD (R(0) and R(1)), H2O, CH4, and CO are present, but so far no new species have been found. Our results indicate that (i) Neptune's mean thermal profile is warmer by ~ 3 K than inferred from the Voyager radio-occultation; (ii) the D/H mixing ratio is (4.5+/-1) X 10**-5, confirming the enrichment of Neptune in deuterium over the protosolar value (~ 2.1 X 10**-5); (iii) the CH4 mixing ratio in the mid stratosphere is (1.5+/-0.2) X 10**-3, and CH4 appears to decrease in the lower stratosphere at a rate consistent with local saturation, in agreement with the scenario of CH4 stratospheric injection from Neptune's warm south polar region; (iv) the H2O stratospheric column is (2.1+/-0.5) X 10**14 cm-2 but its vertical distribution is still to be determined, so the H2O external flux remains uncertain by over an order of magnitude; and (v) the CO stratospheric abundance is about twice the tropospheric value, confirming the dual origin of CO suspected from ground-based millimeter/submillimeter observations.

We have obtained the first continuous disk averaged spectrum of Mars from 450 to 1550 Ghz using the Herschel-SPIRE Fourier Transform Spectrometer. The spectrum was obtained at a constant resolution of 1.4 GHz across the whole band. The flux from the planet is such that the instrument was operated in "bright source" mode to prevent saturation of the detectors. This was the first successful use of this mode and in this work we describe the method used for observing Mars together with a detailed discussion of the data reduction techniques required to calibrate the spectrum. We discuss the calibration accuracy obtained and describe the first comparison with surface and atmospheric models. In addition to a direct photometric measurement of the planet the spectrum contains the characteristic transitions of 12CO from J 5-4 to J 13-12 as well as numerous H2O transitions. Together these allow the comparison to global atmospheric models allowing the mean mixing ratios of water and 12CO to be investigated. We find that it is possible to match the observed depth of the absorption features in the spectrum with a fixed water mixing ratio of 1 x 10-4 and a 12CO mixing ratio of 9 x 10-4

We summarize the first results from the Gould Belt survey, obtained toward the Aquila Rift and Polaris Flare regions during the 'science demonstration phase' of Herschel. Our 70-500 micron images taken in parallel mode with the SPIRE and PACS cameras reveal a wealth of filamentary structure, as well as numerous dense cores embedded in the filaments. Between ~ 350 and 500 prestellar cores and ~ 45-60 Class 0 protostars can be identified in the Aquila field, while ~ unbound starless cores and no protostars are observed in the Polaris field. The prestellar core mass function (CMF) derived for the Aquila region bears a strong resemblance to the stellar initial mass function (IMF), already confirming the close connection between the CMF and the IMF with much better statistics than earlier studies. Comparing and contrasting our Herschel results in Aquila and Polaris, we propose an observationally-driven scenario for core formation according to which complex networks of long, thin filaments form first within molecular clouds, and then the densest filaments fragment into a number of prestellar cores via gravitational instability.

The interstellar medium is enriched primarily by matter ejected from evolved low and intermediate mass stars. The outflows from these stars create a circumstellar envelope in which a rich gas-phase and dust-nucleation chemistry takes place. We observed the nearest carbon-rich evolved star, IRC+10216, using the PACS (55-210 {\mu}m) and SPIRE (194-672 {\mu}m) spectrometers on board Herschel. We find several tens of lines from SiS and SiO, including lines from the v=1 vibrational level. For SiS these transitions range up to J=124-123, corresponding to energies around 6700K, while the highest detectable transition is J=90-89 for SiO, which corresponds to an energy around 8400K. Both species trace the dust formation zone of IRC+10216, and the broad energy ranges involved in their detected transitions permit us to derive the physical properties of the gas and the particular zone in which each species has been formed. This allows us to check the accuracy of chemical thermodynamical equilibrium models and the suggested depletion of SiS and SiO due to accretion onto dust grains.

We report on the detection of anhydrous hydrochloric acid (hydrogen chlorine, HCl) in the carbon-rich star IRC+10216 using the spectroscopic facilities onboard the Herschel satellite. Lines from J=1-0 up to J=7-6 have been detected. From the observed intensities, we conclude that HCl is produced in the innermost layers of the circumstellar envelope with an abundance relative to H2 of 5x10^-8 and extends until the molecules reach its photodissociation zone. Upper limits to the column densities of AlH, MgH, CaH, CuH, KH, NaH, FeH, and other diatomic hydrides have also been obtained.

We obtained Herschel PACS and SPIRE images of the thermal emission of the debris disk around the A5V star {\beta} Pic. The disk is well resolved in the PACS filters at 70, 100, and 160 {\mu}m. The surface brightness profiles between 70 and 160 {\mu}m show no significant asymmetries along the disk, and are compatible with 90% of the emission between 70 and 160 {\mu}m originating in a region closer than 200 AU to the star. Although only marginally resolving the debris disk, the maps obtained in the SPIRE 250 - 500 {\mu}m filters provide full-disk photometry, completing the SED over a few octaves in wavelength that had been previously inaccessible. The small far-infrared spectral index ({\beta} = 0.34) indicates that the grain size distribution in the inner disk (<200AU) is inconsistent with a local collisional equilibrium. The size distribution is either modified by non-equilibrium effects, or exhibits a wavy pattern, caused by an under-abundance of impactors which have been removed by radiation pressure.

The Photodetector Array Camera and Spectrometer (PACS) is one of the three science instruments on ESA's far infrared and submillimetre observatory. It employs two Ge:Ga photoconductor arrays (stressed and unstressed) with 16x25 pixels, each, and two filled silicon bolometer arrays with 16x32 and 32x64 pixels, respectively, to perform integral-field spectroscopy and imaging photometry in the 60-210\mu\ m wavelength regime. In photometry mode, it simultaneously images two bands, 60-85\mu\ m or 85-125\mu\m and 125-210\mu\ m, over a field of view of ~1.75'x3.5', with close to Nyquist beam sampling in each band. In spectroscopy mode, it images a field of 47"x47", resolved into 5x5 pixels, with an instantaneous spectral coverage of ~1500km/s and a spectral resolution of ~175km/s. We summarise the design of the instrument, describe observing modes, calibration, and data analysis methods, and present our current assessment of the in-orbit performance of the instrument based on the Performance Verification tests. PACS is fully operational, and the achieved performance is close to or better than the pre-launch predictions.

We present five band imaging of the Vega debris disc obtained using the Herschel Space Observatory. These data span a wavelength range of 70-500 um with full-width half-maximum angular resolutions of 5.6-36.9". The disc is well resolved in all bands, with the ring structure visible at 70 and 160 um. Radial profiles of the disc surface brightness are produced, and a disc radius of 11" (~ 85 AU) is determined. The disc is seen to have a smooth structure thoughout the entire wavelength range, suggesting that the disc is in a steady state, rather than being an ephemeral structure caused by the recent collision of two large planetesimals.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

High-resolution far-infrared and sub-millimetre spectroscopy of water lines is an important tool to understand the physical and chemical properties of cometary atmospheres. We present observations of several rotational ortho- and para-water transitions in comet C/2008 Q3 (Garradd) performed with HIFI on Herschel. These observations have provided the first detection of the 2_{12}-1_{01} (1669 GHz) ortho and 1_{11}-0_{00} (1113 GHz) para transitions of water in a cometary spectrum. In addition, the ground-state transition 1_{10}-1_{01} at 557 GHz is detected and mapped. By detecting several water lines quasi-simultaneously and mapping their emission we can constrain the excitation parameters in the coma. Synthetic line profiles are computed using excitation models which include excitation by collisions, solar infrared radiation, and radiation trapping. We obtain the gas kinetic temperature, constrain the electron density profile, and estimate the coma expansion velocity by analyzing the map and line shapes. We derive water production rates of 1.7-2.8 x 10^{28} s^{-1} over the range r_h = 1.83-1.85 AU.

With a luminosity > 10^5 Lsun and a mass-loss rate of about 2.10-4 Msun/yr, the red supergiant VY CMa truly is a spectacular object. Because of its extreme evolutionary state, it could explode as supernova any time. Studying its circumstellar material, into which the supernova blast will run, provides interesting constraints on supernova explosions and on the rich chemistry taking place in such complex circumstellar envelopes. We have obtained spectroscopy of VYCMa over the full wavelength range offered by the PACS and SPIRE instruments of Herschel, i.e. 55 to 672 micron. The observations show the spectral fingerprints of more than 900 spectral lines, of which more than half belong to water. In total, we have identified 13 different molecules and some of their isotopologues. A first analysis shows that water is abundantly present, with an ortho-to-para ratio as low as 1.3:1, and that chemical non-equilibrium processes determine the abundance fractions in the inner envelope.

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.

Comet C/2006 W3 (Christensen) was observed in November 2009 at 3.3 AU from the Sun with Herschel. The PACS instrument acquired images of the dust coma in 70- and 160-micrometers filters, and spectra covering several H2O rotational lines. Spectra in the range 450-1550 GHz were acquired with SPIRE. The comet emission continuum from 70 to 672 micrometers was measured, but no lines were detected. The spectral energy distribution indicates thermal emission from large particles and provides a measure of the size distribution index and dust production rate. The upper limit to the water production rate is compared to the production rates of other species (CO, CH3OH, HCN, H2S, OH) measured with the IRAM 30-m and Nancay telescopes. The coma is found to be strongly enriched in species more volatile than water, in comparison to comets observed closer to the Sun. The CO to H2O production rate ratio exceeds 220%. The dust to gas production rate ratio is on the order of 1.

Herschel PACS and SPIRE images have been obtained of NGC 6720 (the Ring Nebula). This is an evolved planetary nebula with a central star that is currently on the cooling track, due to which the outer parts of the nebula are recombining. From the PACS and SPIRE images we conclude that there is a striking resemblance between the dust distribution and the H2 emission, which appears to be observational evidence that H2 forms on grain surfaces. We have developed a photoionization model of the nebula with the Cloudy code which we used to determine the physical conditions of the dust and investigate possible formation scenarios for the H2. We conclude that the most plausible scenario is that the H2 resides in high density knots which were formed after the recombination of the gas started when the central star entered the cooling track. Hydrodynamical instabilities due to the unusually low temperature of the recombining gas are proposed as a mechanism for forming the knots. H2 formation in the knots is expected to be substantial after the central star underwent a strong drop in luminosity about one to two thousand years ago, and may still be ongoing at this moment, depending on the density of the knots and the properties of the grains in the knots.

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.

Polycyclic aromatic hydrocarbons (PAHs) produce characteristic infrared emission bands that have been observed in a wide range of astrophysical environments, where carbonaceous material is subjected to ultraviolet (UV) radiation. Although PAHs are expected to form in carbon-rich AGB stars, they have up to now only been observed in binary systems where a hot companion provides a hard radiation field. In this letter, we present low-resolution infrared spectra of four S-type AGB stars, selected from a sample of 90 S-type AGB stars observed with the infrared spectrograph aboard the Spitzer satellite. The spectra of these four stars show the typical infrared features of PAH molecules. We confirm the correlation between the temperature of the central star and the centroid wavelength of the 7.9 {\mu}m feature, present in a wide variety of stars spanning a temperature range from 3 000 to 12 000 K. Three of four sources presented in this paper extend this relation towards lower temperatures. We argue that the mixture of hydrocarbons we see in these S-stars has a rich aliphatic component. The fourth star, BZ CMa, deviates from this correlation. Based on the similarity with the evolved binary TU Tau, we predict that BZ CMa has a hot companion as well.

Aims. We present a diagnostic tool to determine the abundance of the crystalline silicate forsterite in AGB stars surrounded by a thick shell of silicate dust. Using six infrared spectra of high mass-loss oxygen rich AGB stars we obtain the forsterite abundance of their dust shells. Methods. We use a monte carlo radiative transfer code to calculate infrared spectra of dust enshrouded AGB stars. We vary the dust composition, mass-loss rate and outer radius. We focus on the strength of the 11.3 and the 33.6 \mu m forsterite bands, that probe the most recent (11.3 \mu m) and older (33.6 \mu m) mass-loss history of the star. Simple diagnostic diagrams are derived, allowing direct comparison to observed band strengths. Results. Our analysis shows that the 11.3 \mu m forsterite band is a robust indicator for the forsterite abundance of the current mass-loss period for AGB stars with an optically thick dust shell. The 33.6 \mu m band of forsterite is sensitive to changes in the density and the geometry of the emitting dust shell, and so a less robust indicator. Applying our method to six high mass-loss rate AGB stars shows that AGB stars can have forsterite abundances of 12% by mass and higher, which is more than the previously found maximum abundance of 5%.

Aims: We aim at measuring mass-loss rates and the luminosities of a statistically large sample of Galactic bulge stars at several galactocentric radii. The sensitivity of previous infrared surveys of the bulge has been rather limited, thus fundamental questions for late stellar evolution, such as the stage at which substantial mass-loss begins on the red giant branch and its dependence on fundamental stellar properties, remain unanswered. We aim at providing evidence and answers to these questions. Methods: To this end, we observed seven 15 times 15 arcmin^2 fields in the nuclear bulge and its vicinity with unprecedented sensitivity using the IRAC and MIPS imaging instruments on-board the Spitzer Space Telescope. In each of the fields, tens of thousands of point sources were detected. Results: In the first paper based on this data set, we present the observations, data reduction, the final catalogue of sources, and a detailed comparison to previous mid-IR surveys of the Galactic bulge, as well as to theoretical isochrones. We find in general good agreement with other surveys and the isochrones, supporting the high quality of our catalogue.

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.

Context: Among late-type red giants, an interesting change occurs in the structure of the outer atmospheric layers as one moves to later spectral types in the Hertzsprung-Russell diagram: a chromosphere is always present, but the coronal emission diminishes and a cool massive wind steps in. Aims: Where most studies have focussed on short-wavelength observations, this article explores the influence of the chromosphere and the wind on long-wavelength photometric measurements. Methods: The observational spectral energy distributions are compared with the theoretical predictions of the MARCS atmosphere models for a sample of 9 K- and M-giants. The discrepancies found are explained using basic models for flux emission originating from a chromosphere or an ionized wind. Results: For 7 out of 9 sample stars, a clear flux excess is detected at (sub)millimeter and/or centimeter wavelengths. The precise start of the excess depends upon the star under consideration. The flux at wavelengths shorter than about 1 mm is most likely dominated by an optically thick chromosphere, where an optically thick ionized wind is the main flux contributor at longer wavelengths. Conclusions: Although the optical to mid-infrared spectrum of the studied K- and M-giants is well represented by a radiative equilibrium atmospheric model, the presence of a chromosphere and/or ionized stellar wind at higher altitudes dominates the spectrum in the (sub)millimeter and centimeter wavelength ranges. The presence of a flux excess also has implications on the role of these stars as fiducial spectrophotometric calibrators in the (sub)millimeter and centimeter wavelength range.

We present spectroscopic observations from the {\it Spitzer Space Telescope} of six carbon-rich AGB stars in the Sagittarius Dwarf Spheroidal Galaxy (Sgr dSph) and two foreground Galactic carbon stars. The band strengths of the observed C$_2$H$_2$ and SiC features are very similar to those observed in Galactic AGB stars. The metallicities are estimated from an empirical relation between the acetylene optical depth and the strength of the SiC feature. The metallicities are higher than those of the LMC, and close to Galactic values. While the high metallicity could imply an age of around 1 Gyr, for the dusty AGB stars, the pulsation periods suggest ages in excess of 2 or 3 Gyr. We fit the spectra of the observed stars using the DUSTY radiative transfer model and determine their dust mass-loss rates to be in the range 1.0--3.3$\times 10^{-8} $M$_{\odot}$yr$^{-1}$. The two Galactic foreground carbon-rich AGB stars are located at the far side of the solar circle, beyond the Galactic Centre. One of these two stars show the strongest SiC feature in our present Local Group sample.

Asymptotic Giant Branch variables are found to obey period-luminosity relations in the mid-IR similar to those seen at K_S (2.14 microns), even at 24 microns where emission from circumstellar dust is expected to be dominant. Their loci in the M, logP diagrams are essentially the same for the LMC and for NGC6522 in spite of different ages and metallicities. There is no systematic trend of slope with wavelength. The offsets of the apparent magnitude vs. logP relations imply a difference between the two fields of 3.8 in distance modulus. The colours of the variables confirm that a principal period with log P > 1.75 is a necessary condition for detectable mass-loss. At the longest observed wavelength, 24 microns, many semi-regular variables have dust shells comparable in luminosity to those around Miras. There is a clear bifurcation in LMC colour-magnitude diagrams involving 24 micron magnitudes.

METIS, the Mid-infrared ELT Imager and Spectrograph (formerly called MIDIR), is a proposed instrument for the European Extremely Large Telescope (E-ELT), currently undergoing a phase-A study. The study is carried out within the framework of the ESO-sponsored E-ELT instrumentation studies. METIS will be designed to cover the E-ELT science needs at wavelengths longward of 3um, where the thermal background requires different operating schemes. In this paper we discuss the main science drivers from which the instrument baseline has been derived. Specific emphasis has been given to observations that require very high spatial and spectral resolution, which can only be achieved with a ground-based ELT. We also discuss the challenging aspects of background suppression techniques, adaptive optics in the mid-IR, and telescope site considerations. The METIS instrument baseline includes imaging and spectroscopy at the atmospheric L, M, and N bands with a possible extension to Q band imaging. Both coronagraphy and polarimetry are also being considered. However, we note that the concept is still not yet fully consolidated. The METIS studies are being performed by an international consortium with institutes from the Netherlands, Germany, France, United Kingdom, and Belgium.

It is generally acknowledged that the mass loss of Asymptotic Giant Branch (AGB) stars undergoes variations on different time scales. We constructed models for the dust envelopes for a sample of AGB stars to assess whether mass-loss variations influence the spectral energy distribution (SED). To constrain the variability, extra observations at millimetre wavelengths (1.2 mm) were acquired. From the analysis of the dust models, two indications for the presence of mass-loss variations can be found, being (1) a dust temperature at the inner boundary of the dust envelope that is far below the dust condensation temperature and (2) an altered density distribution with respect to rho(r) \propto r^{-2} resulting from a constant mass-loss rate. For 5 out of the 18 studied sources a two-component model of the envelope is required, consisting of an inner region with a constant mass-loss rate and an outer region with a less steep density distribution. For one source an outer region with a steeper density distribution was found. Moreover, in a search for time variability in our data set at 1.2 mm, we found that WX Psc shows a large relative time variation of 34% which might partially be caused by variable molecular line emission.

We present Spitzer Space telescope spectroscopic observations of 14 carbon-rich AGB stars in the Small Magellanic Cloud. SiC dust is seen in most of the carbon-rich stars but it is weak compared to LMC stars. The SiC feature is strong only for stars with significant dust excess, opposite to what is observed for Galactic stars. We argue that in the SMC, SiC forms at lower temperature than graphite dust, whereas the reverse situation occurs in the Galaxy where SiC condenses at higher temperatures and forms first. Dust input into the interstellar medium by AGB stars consists mostly of carbonaceous dust, with little SiC or silicate dust. Only the two coolest stars show a 30-micron band due to MgS dust. We suggest that this is due to the fact that, in the SMC, mass-losing AGB stars generally have low circumstellar (dust) optical depth and therefore effective heating of dust by the central star does not allow temperatures below the 650 K necessary for MgS to exist as a solid. Gas phase C$_2$H$_2$ bands are stronger in the SMC than in the LMC or Galaxy. This is attributed to an increasing C/O ratio at low metallicity. We present a colour-colour diagram based on Spitzer IRAC and MIPS colours to discriminate between O- and C-rich stars. We show that AGB stars in the SMC become carbon stars early in the thermal-pulsing AGB evolution, and remain optically visible for $\sim 6 \times 10^5$ yr. For the LMC, this lifetime is $\sim 3 \times 10^5$ yr. The superwind phase traced with Spitzer lasts for $\sim 10^4$ yr. Spitzer spectra of a K supergiant and a compact HII region are also given.

To study the nature of Bulge AGB stars and in particular their circumstellar dust, we have analysed mid-infrared spectra obtained with the ISOCAM CVF spectrometer in three Bulge fields. The ISOCAM 5-16.5 micron CVF spectra were obtained as part of the ISOGAL infrared survey of the inner Galaxy. A classification of the shape of the 10 micron dust feature was made for each case. The spectra of the individual sources were modelled using a radiative transfer model. Different combinations of amorphous silicates and aluminium-oxide dust were used in the modelling. Spectra were obtained for 29 sources of which 26 are likely to be Bulge AGB stars. Our modelling shows that the stars suffer mass loss rates in the range of 10^{-8} - 5 x 10^{-7} Msun / yr, which is at the low end of the mass-loss rates experienced on the Thermally Pulsing AGB. The luminosities range from 1,700 to 7,700 Lsun as expected for a population of AGB stars with Minit of 1.5 - 2Msun. In agreement with the condensation sequence scenario, we find that the dust content is dominated by Al_2O_3 grains in this sample of low mass-loss rate stars.

We investigate the molecular bands in carbon-rich AGB stars in the Large Magellanic Cloud (LMC), using the InfraRed Spectrograph (IRS) on board the Spitzer Space Telescope (SST) over the 5--38 micron range. All 26 low-resolution spectra show acetylene (C2H2) bands at 7 and 14 micron. The hydrogen cyanide (HCN) bands at these wavelengths are very weak or absent. This is consistent with low nitrogen abundances in the LMC. The observed 14 micron C2H2 band is reasonably reproduced by an excitation temperature of 500 K. There is no clear dilution of the 14 micron band by circumstellar dust emission. This 14 micron band originates from molecular gas in the circumstellar envelope in these high mass-loss rate stars, in agreement with previous findings for Galactic stars. The C2H2,column density, derived from the 13.7 micron band, shows a gas mass-loss rate in the range 3x10^-6 to 5x10^{-5} Msun yr-1. This is comparable with the total mass-loss rate of these stars estimated from the spectral energy distribution. Additionally, we compare the line strengths of the 13.7 micron C2H2 band of our LMC sample with those of a Galactic sample. Despite the low metallicity of the LMC, there is no clear difference in the C2H2 abundance among LMC and Galactic stars. This reflects the effect of the 3rd dredge-up bringing self-produced carbon to the surface, leading to high C/O ratios at low metallicity.

Context. Galactic Centre (GC) OH/IR stars can be, based on the expansion velocities of their circumstellar shells, divided into two groups which are kinematically different and therefore are believed to have evolved from different stellar populations. Aims. To study the metallicity distribution of the OH/IR stars population in the GC on basis of a theoretical relation between EW(Na), EW(Ca) and EW(CO) and the metallicity. Methods. For 70 OH/IR stars in the GC, we obtained near-IR spectra. The equivalent line widths of NaI, CaI, 12CO(2,0) and the curvature of the spectrum around 1.6 micron due to water absorption are determined. Results. The near-IR spectrum of OH/IR stars is influenced by several physical processes. OH/IR stars are variable stars suffering high mass-loss rates. The dust that is formed around the stars strongly influences the near-IR spectra and reduces the equivalent line widths of NaI, CaI. A similar effect is caused by the water content in the outer atmosphere of the OH/IR star. Because of these effects, it is not possible with our low resolution near-infrared spectroscopy to determine the metallicities of these stars.

We present a Spitzer Space Telescope spectroscopic survey of mass-losing carbon stars (and one oxygen-rich star) in the Large Magellanic Cloud. The spectra cover the wavelength range 5--38 micron. They show varying combinations of dust continuum, dust emission features (SiC, MgS) and molecular absorption bands (C2H2, HCN). A set of four narrow bands, dubbed the Manchester system, is used to define the infrared continuum for dusty carbon stars. The relations between the continuum colours and the strength of the dust and molecular features are studied, and are compared to Galactic stars of similar colours. The circumstellar 7-micron C2H2 band is found to be stronger at lower metallicity, from a comparison of stars in the Galaxy, the LMC and the SMC. This is explained by dredge-up of carbon, causing higher C/O ratios at low metallicity (less O). A possible 10-micron absorption feature seen in our spectra may be due to C3. This band has also been identified with interstellar silicate or silicon-nitrite dust. The line-to-continuum ratio of theSiC and MgS bands shows some indication of being lower at low metallicity. The MgS band is only seen at dust temperatures below 600 K. We discuss the selection of carbon versus oxygen-rich AGB stars using the J-K vs. K-A colours, and show that these colours are relatively insensitive to chemical type. Metal-poor carbon stars form amorphous carbon dust from self-produced carbon. The formation efficiency of oxygen-rich dust depends more strongly on metallicity. We suggest that in lower-metallicity environments, the dust input into the Interstellar Medium by AGB stars is efficient but may be strongly biassed towards carbonaceous dust, as compare to the Galaxy.

We present the first thermal-infrared spectra of an extra-galactic Young Stellar Object (YSO), IRAS05328-6827 in the HII region LHA 120-N148 in the Large Magellanic Cloud. The observed and modelled spectral energy distribution reveals a massive YSO, M=20 Msun, which is heavily-embedded and probably still accreting. The reduced dust content as a consequence of the lower metallicity of the LMC allows a unique view into this object, and together with a high C/O ratio may be responsible for the observed low abundance of water ice and relatively high abundances of methanol and CO2 ices.

(abridged) The 222000 I-band light curves of variable stars detected by the OGLE-II survey in the direction of the Galactic Bulge have been fitted and have also been correlated with the DENIS and 2MASS databases. Results are presented for 2691 objects with I-band semi-amplitude larger than 0.45 magnitude, corresponding to classical Mira variables. The Mira period distribution of 6 fields at similar longitude but spanning latitudes from -1.2 to -5.8 degrees are statistically indistinguisable indicating similar populations with initial masses of 1.5-2 \msol (corresponding to ages of 1-3 Gyr). A field at similar longitude at b = -0.05 from Glass et al. (2001) does show a significantly different period distribution, indicating the presence of a younger population of 2.5-3 \msol and ages below 1 Gyr. The K-band period-luminosity relation is presented for the whole sample, and for sub-fields. The zero point depends on Galactic longitude. Simulations are carried out to show that the observed dependence of the zero point with $l$, and the number of stars per field are naturally explained using the model of disk and bulge stars of Binney et al. (1997), for a viewing angle of 43 +/- 17 degrees. A comparison is made with similar objects in the Magellanic Clouds, studied in a previous paper. The slope of the PL-relation in the Bulge and the MCs agree within the errorbars. A distance to the Galactic Centre between 8.6 and 9.0 kpc is derived (for an assumed DM of 18.50 to the LMC).

We investigated the physical and thermal properties of 65 Cybele}, one of the largest main-belt asteroids. Based on published and recently obtained thermal infrared observations, including ISO measurements, we derived through thermophysical modelling (TPM) a size of 302x290x232 km (+/- 4 %) and an geometric visible albedo of 0.050+/-0.005. Our model of a regolith covered surface with low thermal inertia and "default" roughness describes the wavelengths and phase angle dependent thermal aspects very well. Before/after opposition effect and beaming behaviour can be explained in that way. We found a constant emissivity of 0.9 at wavelengths up to about 100 micron and lower values towards the submillimetre range, indicating a grain size distribution dominated by 200 micron particle sizes. The spectroscopic analysis revealed an emissivity increase between 8.0 and 9.5 micron. We compared this emissivity behaviour with the Christiansen features of carbonaceous chondrite meteorites, but a conclusive identification was not possible. A comparison between the Standard Thermal Model (STM) and the applied TPM clearly demonstrates the limitations and problems of the STM for the analysis of multi-epoch and -wavelengths observations. While the TPM produced a unique diameter/albedo solution, the calculated STM values varied by +/-30 % and showed clear trends with wavelength and phase angle. Cybele can be considered as a nice textbook case for the thermophysical analysis of combined optical and thermal infrared observations.

The ISOGAL project is an infrared survey of specific regions sampling the Galactic Plane selected to provide information on Galactic structure,stellar populations,stellar mass-loss and the recent star formation history of the inner disk and Bulge of the Galaxy. ISOGAL combines 7 and 15 micron ISOCAM observations - with a resolution of 6'' at worst - with DENIS IJKs data to determine the nature of the sources and theinterstellar extinction. We have observed about 16 square degrees with a sensitivity approaching 10-20mJy, detecting ~10^5 sources,mostly AGB stars,red giants and young stars. The main features of the ISOGAL survey and the observations are summarized in this paper,together with a brief discussion of data processing and quality. The primary ISOGAL products are described briefly (a full description is given in Schuller et al. 2003, astro-ph/0304309): viz. the images and theISOGAL-DENIS five-wavelength point source catalogue. The main scientific results already derived or in progress are summarized. These include astrometrically calibrated 7 and 15um images,determining structures of resolved sources; identification and properties of interstellar dark clouds; quantification of the infrared extinction law and source dereddening; analysis of red giant and (especially) AGB stellar populations in the central Bulge,determining luminosity,presence of circumstellar dust and mass--loss rate,and source classification,supplemented in some cases by ISO/CVF spectroscopy; detection of young stellar objects of diverse types,especially in the inner Bulge with information about the present and recent star formation rate; identification of foreground sources with mid-IR excess. These results are the subject of about 25 refereed papers published or in preparation.

We present version 1.0 of the ISOGAL-DENIS Point Source Catalogue (PSC), containing more than 100,000 point sources detected at 7 and/or 15 micron in the ISOGAL survey of the inner Galaxy with the ISOCAM instrument on board the Infrared Space Observatory (ISO). These sources are cross-identified, wherever possible, with near-infrared (0.8-2.2 micron) data from the DENIS survey. The overall surface covered by the ISOGAL survey is about 16 square degrees, mostly (95%) distributed near the Galactic plane (|b| < 1 deg), where the source extraction can become confusion limited and perturbed by the high background emission. Therefore, special care has been taken aimed at limiting the photometric error to ~0.2 magnitude down to a sensitivity limit of typically 10 mJy. The present paper gives a complete description of the entries and the information which can be found in this catalogue, as well as a detailed discussion of the data processing and the quality checks which have been completed. The catalogue is available via the VizieR Service at the Centre de Donn\'ees Astronomiques de Strasbourg (CDS, http://vizier.u-strasbg.fr/viz-bin/VizieR/) and also via the server at the Institut d'Astrophysique de Paris (http://www-isogal.iap.fr/). A more complete version of this paper, including a detailed description of the data processing, is available in electronic form through the ADS service.

This article presents the study of the light-curves extracted from the MACHO database of a sample of stars observed by the ISO in the SMC. These stars belong to the ISO-Mini-Survey catalogue of the Magellanic Clouds (ISO-MCMS, Loup et al. in preparation). Most of them are in the AGB and supergiant phases. The dominant period and amplitude of pulsation have been derived and the stars have been classified as Mira or Semi-Regular pulsators. Furthermore, the cross-identification with near-infrared DENIS and 2MASS magnitudes available within the ISO-MCMS allowed us: (i) to investigate the properties of these stars in the combined near- and mid-infrared colour-magnitude diagrams, (ii) to derive the bolometric magnitude by integrating the spectral energy distribution and (iii) to estimate the mass-loss rate. The stars have been divided into carbon- (C-) and oxygen-rich (O-rich) using the (J-Ks, Ks) colour-magnitude diagram and their period and amplitude distributions have been compared. C-rich AGB stars have a sharp peak in their period distribution at about 250 days and have on average a larger amplitude than O-rich AGB stars. This effect, not previously detected from the study of similar stars in the LMC and in the Baade's window, might be closely related to the metallicity of the environment in which the stars have formed.

Near- and mid-IR survey data from DENIS and ISOGAL are used to investigate the structure and formation history of the inner 10 degree (1.4 kpc) of the Milky Way galaxy. Synthetic bolometric corrections and extinction coefficients in the near- and mid-IR are derived for stars of different spectral types, to allow the transformation of theoretical isochrones into observable colour-magnitude diagrams. The observed IR colour-magnitude diagrams are used to derive the extinction, metallicity and age for individual stars. The inner galaxy is dominated by an old population (> 7 Gyr). In addition, an intermediate-age population (200 Myr to 7 Gyr) is detected, which is consistent with the presence of a few hundred Asymptotic Giant Branch stars with heavy mass loss. Furthermore, young stars (< 200 Myr) are found across the inner Bulge. The metallicities of these stellar population components are discussed. These results can be interpreted in terms of an early epoch of intense star formation and chemical enrichment which shaped the bulk of the Bulge and nucleus, and a more continuous star formation history which gradually shaped the disk from the accretion of sub-solar metallicity gas from the halo. A possible increase in star formation about 200 Myr ago might have been triggered by a minor merger. Ever since the formation of the first stars, mechanisms have been at play that mix the populations from the nucleus, Bulge and disk. Luminosity functions across the inner galactic plane indicate the presence of an inclined (bar) structure at > 1 kpc from the galactic centre, near the inner Lindblad resonance. The innermost part of the Bulge, within about 1 kpc from the galactic centre, seems azimuthally symmetric.

We present the results of the first extensive mid-infrared (IR) imaging survey of the rho Ophiuchi embedded cluster, performed with the ISOCAM camera on board the ISO satellite. The main molecular cloud L1688, as well as L1689N and L1689S, have been completely surveyed for point sources at 6.7 and 14.3 micron. A total of 425 sources are detected including 16 Class I, 123 Class II, and 77 Class III young stellar objects (YSOs). Essentially all of the mid-IR sources coincide with near-IR sources, but a large proportion of them are recognized for the first time as YSOs. Our dual-wavelength survey allows us to identify essentially all the YSOs with IR excess in the embedded cluster down to Fnu ~ 10 - 15 mJy. It more than doubles the known population of Class II YSOs and represents the most complete census to date of newly formed stars in the rho Ophiuchi central region. The stellar luminosity function of the complete sample of Class II YSOs is derived with a good accuracy down to L= 0.03 Lsun. A modeling of this lumino- sity function, using available pre-main sequence tracks and plausible star for- mation histories, allows us to derive the mass distribution of the Class II YSOs which arguably reflects the IMF of the embedded cluster. We estimate that the IMF in rho Ophiuchi is well described by a two-component power law with a low- mass index of -0.35+/-0.25, a high-mass index of -1.7 (to be compared with the Salpeter value of -1.35), and a break occurring at M = 0.55+/-0.25 Msun. This IMF is flat with no evidence for a low-mass cutoff down to at least 0.06 Msun.

The obscured OH/IR star IRAS05298-6957 in the LMC was recently noticed to be member of the small double cluster HS 327 that also contains a carbon star (van Loon et al., 1998, A&A 329, 169). Hence they are coeval and have (nearly) the same progenitor mass, which can only be understood if Hot Bottom Burning (HBB) has prevented IRAS05298-6957 from being a carbon star. We present extensive visual and near-IR photometric data for >10^4 stars in and around HS 327, and spectroscopic data for some of the brightest AGB stars amongst these. Colour-magnitude diagrams are used to estimate the age for the cluster and its members, and luminosities are derived for the stars for which spectra have been obtained. The age for IRAS05298-6957 and the carbon star is estimated to be about 200 Myr. This corresponds to a Main-Sequence progenitor mass of about 4.0 M_sun --- the first direct measurement of the lower mass threshold for HBB. This agrees with stellar evolution models that, however, fail to reproduce the low luminosity of the carbon star.

By cross-correlating the results of two recent large-scale surveys, the general properties of a well defined sample of semi-regular variable stars have been determined. ISOGAL mid-infrared photometry and MACHO lightcurves are assembled for approximately 300 stars in the Baade's Windows of low extinction towards the Galactic bulge. These stars are mainly giants of late M spectral type, evolving along the asymptotic giant branch (AGB). They are found to possess a wide and continuous distribution of pulsation periods and to obey an approximate log~period -- bolometric magnitude relation or set of such relations. Approximate mass-loss rates in the range of 1e-8 to 5e-7 M_sun per year are derived from ISOGAL mid-infrared photometry and models of stellar spectra adjusted for the presence of optically-thin circumstellar silicate dust. Mass-loss rates depend on luminosity and pulsation period. Some stars lose mass as rapidly as short-period Miras but do not show Mira-like amplitudes. A period of 70 days or longer is a necessary but not a sufficient condition for mass loss to occur. For AGB stars in the mass-loss ranges that we observe, the functional dependence of mass-loss rate on temperature and luminosity is found to be in agreement with recent theoretical predictions. If we include our mass-loss rates with a sample of extreme mass-losing AGB stars in the Large Magellanic Cloud, we get the general result for AGB stars that mass-loss rate is proportional to luminosity^{2.7}, valid for AGB stars with 10^{-8} to 10^{-4} M_sun per year (Abridged).

In an attempt to determine the nature of the enigmatic cocoon stars in the Quintuplet Cluster, we have obtained mid-infrared imaging and spectrophotometry of the cluster, using the CAM and SWS instruments on ISO, using SpectroCam-10 on the Palomar 5m telescope, and NICMOS on HST. The spectra show smooth continua with various dust and ice absorption features. These features are all consistent with an interstellar origin, and there is no clear evidence for any circumstellar contribution to these features. We find no spectral line or feature that could elucidate the nature of these sources. Detailed modeling of the silicate absorption features shows that they are best reproduced by the mu Cep profile, which is typical of the interstellar medium, with tau(sil) \sim 2.9. The high spatial resolution mid-IR images show that three of the five cocoon stars have spatially extended and asymmetric envelopes, with diameters of \sim 20,000 AUs. A reddening law similar to that of Lutz (1999) but with silicate features based on the mu Cep profile and normalized to our value of tau(sil) is used to deredden the observed spectrophotometry. The dereddened energy distributions are characterised by temperatures of 750-925 K, somewhat cooler than determined from near IR data alone. Models of optically thin and geometrically thick dust shells, as used by Williams et al. (1987) for very dusty, late-type WC stars, reproduce the observed SEDs from 4 to 17 mic, and imply shell luminosities of log(L/L(sun)) \sim 4.5-4.9 for the brightest four components. An analysis of the various suggestions proposed to explain the nature of the cocoon stars reveals serious problems with all the hypotheses, and the nature of these sources remains an enigma.

The ISOGAL project is a survey of the stellar populations, structure, and recent star formation history of the inner disk and bulge of the Galaxy. ISOGAL combines 15 and 7micron ISOCAM observations with DENIS IJKs data to determine the nature of a source and the interstellar extinction. In this paper we report an ISOGAL study of a small field in the inner Galactic Bulge (l=0deg, b=1.0deg, area=0.035 sq. deg) as a prototype of the larger area ISOGAL survey of the inner Galaxy. The five wavelengths of ISOGAL+DENIS, together with the relatively low and constant extinction in front of this specific field, allow reliable determination of the nature of the sources. The primary scientific result of this paper is evidence that the most numerous class of ISOGAL 15micron sources are Red Giants in the Galactic bulge and central disk, with luminosities just above or close to the RGB tip and weak mass-loss rates. They form loose sequences in the magnitude-colour diagrams [15]/Ks-[15] and [15]/[7]-[15]. Their large excesses at 15micron with respect to 2micron and 7micron is due to circumstellar dust produced by mass-loss at low rates. These ISOGAL results are the first systematic evidence and study of dust emission at this early stage (''Intermediate'' AGB), before the onset of the large mass-loss phase. It is thus well established that efficient dust formation is already associated with such low mass-loss rates during this early phase.

We present ISO photometric and spectroscopic observations of a sample of 57 bright Asymptotic Giant Branch stars and red supergiants in the Large Magellanic Cloud, selected on the basis of IRAS colours indicative of high mass-loss rates. PHOT-P and PHOT-C photometry at 12, 25 and 60 $\mu$m and CAM photometry at 12 $\mu$m are used in combination with quasi-simultaneous ground-based near-IR photometry to construct colour-colour diagrams for all stars in our sample. PHOT-S and CAM-CVF spectra in the 3 to 14 $\mu$m region are presented for 23 stars. From the colour-colour diagrams and the spectra, we establish the chemical types of the dust around 49 stars in this sample. Many stars have carbon-rich dust. The most luminous carbon star in the Magellanic Clouds has also a (minor) oxygen-rich component. OH/IR stars have silicate absorption with emission wings. The unique dataset presented here allows a detailed study of a representative sample of thermal-pulsing AGB stars with well-determined luminosities.

The ISOGAL mid-infrared survey of areas close to the Galactic Plane aims to determine their stellar content and its possible bearing on the history of the Galaxy. The NGC6522 and Sgr I Baade's Windows of low obscuration towards the inner parts of the Bulge represent ideal places in which to calibrate and understand the ISOGAL colour-magnitude diagrams. The survey observations were made with the ISOCAM instrument of the ISO satellite. The filter bands chosen were LW2(~7 microns) and LW3(~15 microns). The results presented here show that most of the detected objects are late M-type giants on the AGB, with a cut-off for those earlier than M3-M4. The most luminous members of these two Bulge fields at 7 microns are the Mira variables. However, it is evident that they represent the end of a sequence of increasing 15 micron dust emission which commences with M giants of earlier sub-type. In observations of late-type giants the ISOCAM 15 micron band is mainly sensitive to the cool silicate or aluminate dust shells which overwhelm the photospheric emission. However, in ordinary M-giant stars, the 7 micron band is not strongly affected by dust emission and may be influenced instead by absorption. The nu2 band of water at 6.25 microns and the SiO fundamental at 7.9 microns are likely contribitors to this effect. A group of late M stars has been found which vary little or not at all but have infrared colours typical of well-developed dust shells. Their luminosities are similar to those of 200-300 day Miras but they have slightly redder [7]-[15] colours which form an extension of the ordinary M giant sequence. The Mira dust shells show a mid-infrared [7]-[15] colour-period relation. ca 700 days.

We describe ISO observations of the obscured Asymptotic Giant Branch (AGB) star IRAS04496-6958 in the Large Magellanic Cloud (LMC). This star has been classified as a carbon star. Our new ISOCAM CVF spectra show that it is the first carbon star with silicate dust known outside of the Milky Way. The existence of this object, and the fact that it is one of the highest luminosity AGB stars in the LMC, provide important information for theoretical models of AGB evolution and understanding the origin of silicate carbon stars.

We present initial results of a spectroscopic study of the Pistol and of the cocoon stars in the Quintuplet Cluster. From ISOCAM CVF 5--17 micron spectroscopy of the field of the Pistol Star, we have discovered a nearly spherical shell of hot dust surrounding this star, a probable LBV. This shell is most prominent at lambda > 12 micron, and its morphology clearly indicates that the shell is stellar ejecta. Emission line images show that most of the ionised material is along the northern border of this shell, and its morphology is very similar to that of the Pistol HII region (Yusef-Zadeh & Morris, 1987, AJ, 94, 1178). We thus confirm that the ionisation comes from very hot stars in the core of the Quintuplet Cluster. An SWS spectrum of the Pistol Nebula indicates a harder ionising radiation than could be provided by the Pistol Star, but which is consistent with ionisation from Wolf-Rayet stars in the Quintuplet Cluster. The CVF 5--17 micron spectra of the cocoon stars in the Quintuplet do not show any emission feature that could help elucidate their nature.