The eccentric massive binary HD152248 (also known as V1007 Sco), which hosts two O7.5 III-II(f) stars, is the most emblematic eclipsing O-star binary in the very young and rich open cluster NGC6231. Its properties render the system an interesting target for studying tidally induced apsidal motion. Measuring the rate of apsidal motion in such a binary system gives insight into the internal structure and evolutionary state of the stars composing it. A large set of optical spectra was used to reconstruct the spectra of the individual binary components and establish their radial velocities using a disentangling code. Radial velocities measured over seven decades were used to establish the rate of apsidal motion. We furthermore analysed the reconstructed spectra with the CMFGEN model atmosphere code to determine stellar and wind properties of the system. Optical photometry was analysed with the Nightfall binary star code. A complete photometric and radial velocity model was constructed in PHOEBE 2 to determine robust uncertainties. We find a rate of apsidal motion of $(1.843^{+0.064}_{-0.083})\deg$ yr$^{-1}$. The photometric data indicate an orbital inclination of $(67.6^{+0.2}_{-0.1})\deg$ and Roche-lobe filling factors of both stars of about 0.86. Absolute masses of $29.5^{+0.5}_{-0.4}$M$_\odot$ and mean stellar radii of $15.07^{+0.08}_{-0.12}$R$_\odot$ are derived for both stars. We infer an observational value for the internal structure constant of both stars of $0.0010\pm0.0001$. Our in-depth analysis of the massive binary HD152248 and the redetermination of its fundamental parameters can serve as a basis for the construction of stellar evolution models to determine theoretical rates of apsidal motion to be compared with the observational one. In addition, the system hosts two twin stars, which offers a unique opportunity to obtain direct insight into the internal structure of the stars.

Fundamental atomic transition parameters, such as oscillator strengths and wavelengths, play a key role in modelling and understanding the chemical composition of stars in the universe. Despite the significant work under way to produce these parameters for many ions, uncertainties in these parameters remain large and can limit the accuracy of chemical abundance determinations. The Belgian repository of fundamental atomic data and stellar spectra (BRASS) aims to provide a large systematic and homogeneous quality assessment of the atomic data available for quantitative spectroscopy. BRASS shall compare synthetic spectra against extremely high quality observed spectra, at a resolution of ~85000 and signal-noise ratios of ~1000, for around 20 bright BAFGK spectral type stars, in order to evaluate the atomic data available for over a thousand potentially useful spectral lines. A large-scale homogeneous selection of atomic lines is performed by synthesising theoretical spectra of literature atomic lines, for FGK-type stars including the Sun, resulting in a selection of 1091 theoretically deep and unblended lines, in the wavelength range 4200-6800~\AA, which may be suitable for quality assessment. Astrophysical log(gf) values are determined for the 1091 transitions using two commonly employed methods. The agreement of these log(gf) values are used to select well-behaving lines for quality assessment. 845 atomic lines were found to be suitable for quality assessment, of which 408 were found to be robust against systematic differences between analysis methods. Around 54% of the quality-assessed lines were found to have at least one literature log(gf) value in agreement with our derived values, though the remaining values can disagree by as much as 0.5 dex. Only 38% of FeI lines were found to have sufficiently accurate log(gf) values, increasing to ~70-75% for the remaining Fe-group lines.

Red supergiant stars are surrounded by a gaseous and dusty circumstellar environment created by their mass loss which spreads heavy elements into the interstellar medium. The structure and the dynamics of this envelope are crucial to understand the processes driving the red supergiant mass loss and the shaping of the pre-supernova ejecta. We have observed the emission from the CO $J = 2-1$ line from the red supergiant star $\mu$~Cep with the NOEMA interferometer. In the line the synthesized beam was $0.92 \times 0.72$~arcsec ($590 \times 462$~au at 641~pc). The continuum map shows only the unresolved contribution of the free-free emission of the star chromosphere. The continuum-subtracted channel maps reveal a very inhomogeneous and clumpy circumstellar environment. In particular, we detected a bright CO clump, as bright as the central source in the line, at 1.80~arcsec south-west from the star, in the blue channel maps. After a deprojection of the radial velocity assuming two different constant wind velocities, the observations were modelled using the 3D radiative transfer code \textsc{lime} to derive the characteristics of the different structures. We determine that the gaseous clumps observed around $\mu$~Cep are responsible for a mass loss rate of $(4.9 \pm 1.0) \times 10^{-7}~{\rm M}_\odot\,{\rm yr}^{-1}$, in addition to a spatially unresolved wind component with an estimated mass-loss rate of $2.0 \times 10^{-6}~{\rm M}_\odot\,{\rm yr}^{-1}$. Therefore, the clumps have a significant role in $\mu$~Cep's mass loss ($\ge 25 \%$). We cannot exclude that the unresolved central outflow may be made of smaller unresolved clumps.

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

Context: Post-Asymptotic Giant Branch (AGB) binaries are surrounded by stable dusty and gaseous disks similar to the ones around young stellar objects. Whereas significant effort is spent on modeling observations of disks around young stellar objects, the disks around post-AGB binaries receive significantly less attention, even though they pose significant constraints on theories of disk physics and binary evolution. Aims: We want to examine the structure of and phenomena at play in circumbinary disks around post-AGB stars. We continue the analysis of our near-infrared interferometric image of the inner rim of the circumbinary disk around IRAS08544-4431. We want to understand the physics governing this inner disk rim. Methods: We use a radiative transfer model of a dusty disk to reproduce simultaneously the photometry as well as the near-infrared interferometric dataset on IRAS08544-4431. The model assumes hydrostatic equilibrium and takes dust settling self-consistently into account. Results: The best-fit radiative transfer model shows excellent agreement with the spectral energy distribution up to mm wavelengths as well as with the PIONIER visibility data. It requires a rounded inner rim structure, starting at a radius of 8.25 au. However, the model does not fully reproduce the detected over-resolved flux nor the azimuthal flux distribution of the inner rim. While the asymmetric inner disk rim structure is likely to be the consequence of disk-binary interactions, the origin of the additional over-resolved flux remains unclear. Conclusions: As in young stellar objects, the disk inner rim of IRAS08544-4431 is ruled by dust sublimation physics. Additional observations are needed to understand the origin of the extended flux and the azimuthal perturbation at the inner rim of the disk.

Fundamental atomic parameters, such as oscillator strengths, play a key role in modelling and understanding the chemical composition of stars in the universe. Despite the significant work underway to produce these parameters for many astrophysically important ions, uncertainties in these parameters remain large and can propagate throughout the entire field of astronomy. The Belgian repository of fundamental atomic data and stellar spectra (BRASS) aims to provide the largest systematic and homogeneous quality assessment of atomic data to date in terms of wavelength, atomic and stellar parameter coverage. To prepare for it, we first compiled multiple literature occurrences of many individual atomic transitions, from several atomic databases of astrophysical interest, and assessed their agreement. Several atomic repositories were searched and their data retrieved and formatted in a consistent manner. Data entries from all repositories were cross-matched against our initial BRASS atomic line list to find multiple occurrences of the same transition. Where possible we used a non-parametric cross-match depending only on electronic configurations and total angular momentum values. We also checked for duplicate entries of the same physical transition, within each retrieved repository, using the non-parametric cross-match. We report the cross-matched transitions for each repository and compare their fundamental atomic parameters. We find differences in log(gf) values of up to 2 dex or more. We also find and report that ~2% of our line list and Vienna Atomic Line Database retrievals are composed of duplicate transitions. Finally we provide a number of examples of atomic spectral lines with different log(gf) values, and discuss the impact of these uncertain log(gf) values on quantitative spectroscopy. All cross-matched atomic data and duplicate transitions are available to download at brass.sdf.org.

The carriers of the diffuse interstellar bands (DIBs) are largely unidentified molecules ubiquitously present in the interstellar medium (ISM). After decades of study, two strong and possibly three weak near-infrared DIBs have recently been attributed to the C60+ fullerene based on observational and laboratory measurements. There is great promise for the identification of the over 400 other known DIBs, as this result could provide chemical hints towards other possible carriers. In an effort to systematically study the properties of the DIB carriers, we have initiated a new large-scale observational survey: the ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES). The main objective is to build on and extend existing DIB surveys to make a major step forward in characterising the physical and chemical conditions for a statistically significant sample of interstellar lines-of-sight, with the goal to reverse-engineer key molecular properties of the DIB carriers. EDIBLES is a filler Large Programme using the Ultraviolet and Visual Echelle Spectrograph at the Very Large Telescope at Paranal, Chile. It is designed to provide an observationally unbiased view of the presence and behaviour of the DIBs towards early-spectral-type stars whose lines-of-sight probe the diffuse-to-translucent ISM. Such a complete dataset will provide a deep census of the atomic and molecular content, physical conditions, chemical abundances and elemental depletion levels for each sightline. Achieving these goals requires a homogeneous set of high-quality data in terms of resolution (R ~ 70000 -- 100000), sensitivity (S/N up to 1000 per resolution element), and spectral coverage (305--1042 nm), as well as a large sample size (100+ sightlines). In this first paper the goals, objectives and methodology of the EDIBLES programme are described and an initial assessment of the data is provided.

Infrared observations of the dusty, massive Homunculus Nebula around the luminous blue variable $\eta$ Carinae are crucial to characterize the mass-loss history and help constrain the mechanisms leading to the Great Eruption. We present the 2.4 - 670 $\mu$m spectral energy distribution, constructed from legacy ISO observations and new spectroscopy obtained with the {\em{Herschel Space Observatory}}. Using radiative transfer modeling, we find that the two best-fit dust models yield compositions which are consistent with CNO-processed material, with iron, pyroxene and other metal-rich silicates, corundum, and magnesium-iron sulfide in common. Spherical corundum grains are supported by the good match to a narrow 20.2 $\mu$m feature. Our preferred model contains nitrides AlN and Si$_3$N$_4$ in low abundances. Dust masses range from 0.25 to 0.44 $M_\odot$ but $M_{\rm{tot}} \ge$ 45 $M_\odot$ in both cases due to an expected high Fe gas-to-dust ratio. The bulk of dust is within a 5$"$ $\times$ 7$"$ central region. An additional compact feature is detected at 390 $\mu$m. We obtain $L_{\rm{IR}}$ = 2.96 $\times$ 10$^6$ $L_\odot$, a 25\% decline from an average of mid-IR photometric levels observed in 1971-1977. This indicates a reduction in circumstellar extinction in conjunction with an increase in visual brightness, allowing 25-40\% of optical and UV radiation to escape from the central source. We also present an analysis of $^{12}$CO and $^{13}$CO $J = 5-4$ through $9-8$ lines, showing that the abundances are consistent with expectations for CNO-processed material. The [$^{12}$C~{\sc{ii}}] line is detected in absorption, which we suspect originates in foreground material at very low excitation temperatures.

Aims. The magnetic star HD148937 is the only Galactic Of?p star surrounded by a nebula. The structure of this nebula is particularly complex and is composed, from the center out outwards, of a close bipolar ejecta nebula (NGC6164/5), an ellipsoidal wind-blown shell, and a spherically symmetric Stromgren sphere. The exact formation process of this nebula and its precise relation to the star's evolution remain unknown. Methods. We analyzed infrared Spitzer IRS and far-infrared Herschel/PACS observations of the NGC6164/5 nebula. The Herschel imaging allowed us to constrain the global morphology of the nebula. We also combined the infrared spectra with optical spectra of the central star to constrain its evolutionary status. We used these data to derive the abundances in the ejected material. To relate this information to the evolutionary status of the star, we also determined the fundamental parameters of HD148937 using the CMFGEN atmosphere code. Results. The H$\alpha$ image displays a bipolar or "8"-shaped ionized nebula, whilst the infrared images show dust to be more concentrated around the central object. We determine nebular abundance ratios of N/O = 1.06 close to the star, and N/O = 1.54 in the bright lobe constituting NGC6164. Interestingly, the parts of the nebula located further from HD148937 appear more enriched in stellar material than the part located closer to the star. Evolutionary tracks suggest that these ejecta have occured $\sim$1.2-1.3 and $\sim$0.6 Myrs ago, respectively. In addition, we derive abundances of argon for the nebula compatible with the solar values and we find a depletion of neon and sulfur. The combined analyses of the known kinematics and of the new abundances of the nebula suggest either a helical morphology for the nebula, possibly linked to the magnetic geometry, or the occurrence of a binary merger.

Gaia is a space mission currently measuring the five astrometric parameters as well as spectrophotometry of at least 1 billion stars to G = 20.7 mag with unprecedented precision. The sixth parameter in phase space (radial velocity) is also measured thanks to medium-resolution spectroscopy being obtained for the 150 million brightest stars. During the commissioning phase, two fields, one around each ecliptic pole, have been repeatedly observed to assess and to improve the overall satellite performances as well as the associated reduction and analysis software. A ground-based photometric and spectroscopic survey was therefore initiated in 2007, and is still running in order to gather as much information as possible about the stars in these fields. This work is of particular interest to the validation of the Radial Velocity Spectrometer (RVS) outputs. The paper presents the radial velocity measurements performed for the Southern targets in the 12 - 17 R magnitude range on high- to mid-resolution spectra obtained with the GIRAFFE and UVES spectrographs.

Using the PACS and SPIRE spectrometers on-board the Herschel Space Observatory, we obtained spectra of two red supergiants (RSGs) in the Large Magellanic Cloud (LMC). Multiple rotational CO emission lines (J=6-5 to 15-14) and 15 H2O lines were detected from IRAS 05280-6910, and one CO line was detected from WOH G64. This is the first time CO rotational lines have been detected from evolved stars in the LMC. Their CO line intensities are as strong as those of the Galactic RSG, VY CMa. Modelling the CO lines and the spectral energy distribution results in an estimated mass-loss rate for IRAS 05280-6910 of 3x10^-4 Msun per yr. The model assumes a gas-to-dust ratio and a CO-to-H2 abundance ratio is estimated from the Galactic values scaled by the LMC metallicity ([Fe/H]~-0.3), i.e., that the CO-to-dust ratio is constant for Galactic and LMC metallicities within the uncertainties of the model. The key factor determining the CO line intensities and the mass-loss rate found to be the stellar luminosity.

Aims. We obtained Herschel/PACS imaging and spectroscopic observations of the nebula surrounding the candidate LBV HD168625. By combining these data with optical spectra of the star, we want to constrain the abundances in the nebula and in the star and compare them to trace back its evolution. Methods. We use CMFGEN to determine the fundamental parameters and the CNO abundances of the central star whilst the abundances of the nebula are derived from the emission lines present in the Herschel/PACS spectrum. Result. The FIR images show a nebula composed of an elliptical ring/torus of ejecta with a ESE-WNW axis and of a second perpendicular bipolar structure composed of empty caps/rings. We detect equatorial shells composed of dust and ionized material with different sizes when observed at different wavelengths, and bipolar caps more of less separated from the central star in H$\alpha$ and mid-IR images. This complex global structure seems to show 2 different inclinations: 40{\deg} for the equatorial torus and 60{\deg} for the bipolar ejections. From the Herschel/PACS spectrum, we determine nebular abundances of $\mathrm{N/H} = 4.1 \pm 0.8 \times~10^{-4}$ and $\mathrm{C/H} \sim 1.6_{-0.35}^{+1.16} \times~10^{-4}$, as well as a mass of ionized gas of $0.17\pm 0.04$M$_{\odot}$ and a neutral hydrogen mass of about $1.0\pm 0.3$M$_{\odot}$ which dominates. Analysis of the central star reveals Teff = 14000K, log g = 1.74 and log (L/L$_{\odot}$) = 5.58. We derive stellar abundances of about $\mathrm{N/H} = 5.0 \pm 1.5 \times~10^{-4}$, $\mathrm{C/H} = 1.4 \pm 0.5 \times~10^{-4}$ and $\mathrm{O/H} = 3.5 \pm 1.0\times~10^{-4}$. All these measurements taken together are compatible with the evolutionary tracks of a star with an initial mass between $28$ and $33$M$_{\odot}$ and with a crit. rotational rate between 0.3 and 0.4 that has lost its material during or just after the BSG phase.

We aim to constrain the temperature and velocity structures, and H2O abundances in the winds of a sample of M-type AGB stars. We further aim to determine the effect of H2O line cooling on the energy balance in the inner circumstellar envelope. We use two radiative-transfer codes to model molecular emission lines of CO and H2O towards four M-type AGB stars. We focus on spectrally resolved observations of CO and H2O from HIFI. The observations are complemented by ground-based CO observations, and spectrally unresolved CO and H2O observations with PAC. The observed line profiles constrain the velocity structure throughout the circumstellar envelopes (CSEs), while the CO intensities constrain the temperature structure in the CSEs. The H2O observations constrain the o-H2O and p-H2O abundances relative to H2. Finally, the radiative-transfer modelling allows to solve the energy balance in the CSE, in principle including also H2O line cooling. The fits to the line profiles only set moderate constraints on the velocity profile, indicating shallower acceleration profiles in the winds of M-type AGB stars than predicted by dynamical models, while the CO observations effectively constrain the temperature structure. Including H2O line cooling in the energy balance was only possible for the low-mass-loss-rate objects in the sample, and required an ad hoc adjustment of the dust velocity profile in order to counteract extreme cooling in the inner CSE. H2O line cooling was therefore excluded from the models. The constraints set on the temperature profile by the CO lines nevertheless allowed us to derive H2O abundances. The derived H2O abundances confirm previous estimates and are consistent with chemical models. However, the uncertainties in the derived abundances are relatively large, in particular for p-H2O, and consequently the derived o/p-H2O ratios are not well constrained.

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.

Infrared Herschel imaging and spectroscopic observations of the nebula M1-67 around the Wolf-Rayet star WR 124 have been obtained along with optical imaging observations. The infrared images reveal a clumpy dusty nebula that extends up to 1 pc. The comparison with the optical images shows that the ionized gas nebula coincides with the dust nebula, the dust and the gas being mixed together. A photodissociation region is revealed from the infrared spectroscopic analysis. The analysis of the infrared spectrum of the nebula, where forbidden emission lines of ionized elements were detected, showed that the nebula consists of mildly processed material with the calculated abundance number ratios being N/O = 1.0 +/- 0.5 and C/O = 0.46 +/- 0.27. Based on a radiative transfer model, the dust mass of the nebula was estimated to be 0.22 Msun with a population of large grains being necessary to reproduce the observations. The comparison of the mass-loss rate and the abundance ratios to theoretical models of stellar evolution led to the conclusion that the nebular ejection took place during a RSG/YSG evolutionary phase of a central star with an initial mass of 32 Msun.

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.

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.

Far-infrared Herschel PACS imaging and spectroscopic observations of the nebula around the luminous blue variable (LBV) star AG Car have been obtained along with optical imaging in the Halpha+[NII] filter. In the infrared light, the nebula appears as a clumpy ring shell that extends up to 1.2 pc with an inner radius of 0.4 pc. It coincides with the Halpha nebula, but extends further out. Dust modeling of the nebula was performed and indicates the presence of large grains. The dust mass is estimated to be ~ 0.2 Msun. The infrared spectrum of the nebula consists of forbidden emission lines over a dust continuum. Apart from ionized gas, these lines also indicate the existence of neutral gas in a photodissociation region that surrounds the ionized region. The abundance ratios point towards enrichment by processed material. The total mass of the nebula ejected from the central star amounts to ~ 15 Msun, assuming a dust-to-gas ratio typical of LBVs. The abundances and the mass-loss rate were used to constrain the evolutionary path of the central star and the epoch at which the nebula was ejected, with the help of available evolutionary models. This suggests an ejection during a cool LBV phase for a star of ~ 55 Msun with little rotation.

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.

S-type AGB stars have a C/O ratio which suggests that they are transition objects between oxygen-rich M-type stars and carbon-rich C-type stars. As such, their circumstellar compositions of gas and dust are thought to be sensitive to their precise C/O ratio, and it is therefore of particular interest to examine their circumstellar properties. We present new Herschel HIFI and PACS sub-millimetre and far-infrared line observations of several molecular species towards the S-type AGB star W Aql. We use these observations, which probe a wide range of gas temperatures, to constrain the circumstellar properties of W Aql, including mass-loss rate and molecular abundances. We used radiative transfer codes to model the circumstellar dust and molecular line emission to determine circumstellar properties and molecular abundances. We assumed a spherically symmetric envelope formed by a constant mass-loss rate driven by an accelerating wind. Our model includes fully integrated H2O line cooling as part of the solution of the energy balance. We detect circumstellar molecular lines from CO, H2O, SiO, HCN, and, for the first time in an S-type AGB star, NH3. The radiative transfer calculations result in an estimated mass-loss rate for W Aql of 4.0e-6 Msol yr-1 based on the 12CO lines. The estimated 12CO/13CO ratio is 29, which is in line with ratios previously derived for S-type AGB stars. We find an H2O abundance of 1.5e-5, which is intermediate to the abundances expected for M and C stars, and an ortho/para ratio for H2O that is consistent with formation at warm temperatures. We find an HCN abundance of 3e-6, and, although no CN lines are detected using HIFI, we are able to put some constraints on the abundance, 6e-6, and distribution of CN in W Aql's circumstellar envelope using ground-based data. We find an SiO abundance of 3e-6, and an NH3 abundance of 1.7e-5, confined to a small envelope.

The Helix Nebula (NGC 7293) is the closest planetary nebulae. Therefore, it is an ideal template for photochemical studies at small spatial scales in planetary nebulae. We aim to study the spatial distribution of the atomic and the molecular gas, and the structure of the photodissociation region along the western rims of the Helix Nebula as seen in the submillimeter range with Herschel. We use 5 SPIRE FTS pointing observations to make atomic and molecular spectral maps. We analyze the molecular gas by modeling the CO rotational lines using a non-local thermodynamic equilibrium (non-LTE) radiative transfer model. For the first time, we have detected extended OH+ emission in a planetary nebula. The spectra towards the Helix Nebula also show CO emission lines (from J= 4 to 8), [NII] at 1461 GHz from ionized gas, and [CI] (2-1), which together with the OH+ lines, trace extended CO photodissociation regions along the rims. The estimated OH+ column density is (1-10)x1e12 cm-2. The CH+ (1-0) line was not detected at the sensitivity of our observations. Non-LTE models of the CO excitation were used to constrain the average gas density (n(H2)=(1-5)x1e5 cm-3) and the gas temperature (Tk= 20-40 K). The SPIRE spectral-maps suggest that CO arises from dense and shielded clumps in the western rims of the Helix Nebula whereas OH+ and [CI] lines trace the diffuse gas and the UV and X-ray illuminated clumps surface where molecules reform after CO photodissociation. [NII] traces a more diffuse ionized gas component in the interclump medium.

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...

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.

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.

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.

We have obtained far-infrared Herschel PACS imaging and spectroscopic observations of the nebular environment of the luminous blue variable WRAY 15-751. These images clearly show that the main, dusty nebula is a shell of radius 0.5 pc and width 0.35 pc extending outside the H-alpha nebula. They also reveal a second, bigger and fainter dust nebula, observed for the first time. Both nebulae lie in an empty cavity, likely the remnant of the O-star wind bubble formed when the star was on the main sequence. The kinematic ages of the nebulae are about 20000 and 80000 years and each nebula contains about 0.05 Msun of dust. Modeling of the inner nebula indicates a Fe-rich dust. The far-infrared spectrum of the main nebula revealed forbidden emission lines coming from ionized and neutral gas. Our study shows that the main nebula consists of a shell of ionized gas surrounded by a thin photodissociation region illuminated by an "average" early-B star. The derived abundance ratios N/O=1.0+/-0.4 and C/O=0.4+/-0.2 indicate a mild N/O enrichment. We estimate that the inner shell contains 1.7+/-0.6 Msun of gas. Assuming a similar dust-to-gas ratio for the outer nebula, the total mass ejected by WRAY 15-751 amounts to 4+/-2 Msun. The measured abundances, masses and kinematic ages of the nebulae were used to constrain the evolution of the star and the epoch at which the nebulae were ejected. Our results point to an ejection of the nebulae during the RSG evolutionary phase of an ~ 40 Msun star. The presence of multiple shells around the star suggests that the mass-loss was not a continuous ejection but rather a series of episodes of extreme mass-loss. Our measurements are compatible with the recent evolutionary tracks computed for an 40 Msun star with little rotation. They support the O-BSG-RSG-YSG-LBV filiation and the idea that high-luminosity and low-luminosity LBVs follow different evolutionary paths.

We report on the Herschel Gould Belt survey (HGBS) of ChaII, focusing on the detection of Class I to III young stellar objects (YSOs). We aim at characterizing the circumstellar material around these YSOs and understanding which disk parameters are most likely constrained by the new HGBS data. We recovered 29 out of the 63 known YSOs in ChaII with a detection in at least one of the PACS/SPIRE pass-bands: 3 Class I YSOs (i.e.,100%), 1 Flat source (i.e., 50%), 21 Class II objects (i.e., 55%), 3 Class III objects (i.e, 16%) and the unclassified source IRAS 12522-7640. We explore PACS/SPIRE colors of this sample and present modeling of their SEDs using the RADMC-2D radiative transfer code. We find that YSO colors are confined in specific regions of PACS/SPIRE color-color diagrams. These color ranges are expected to be only marginally contaminated by extragalactic sources and field stars and, hence, provide a useful YSO selection tool when applied altogether. We were able to model the SED of 26 out of the 29 detected YSOs. We discuss the degeneracy/limitations of our SED fitting results and adopt the Bayesian method to estimate the probability of different values for the derived disk parameters. The Cha II YSOs present typical disk inner radii around 0.1 AU. The lower limit to Rc is typically around 50 AU. The lower limits to Mdisk are proportional to the stellar masses with a typical 0.3% ratio, i.e., in the range estimated in the literature for young Class II stars and brown dwarfs across a broad range of stellar masses. The estimated flaring angles, although very uncertain, point towards rather flat disks (1+phi less than 1.2), as found for low-mass M-type YSO samples in other star forming regions. Thus, our results support the idea that disk properties show a dependence on stellar properties.

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]

We investigate the effect of the "high-pass filter" data reduction technique on the Herschel PACS PSF and noise of the PACS maps at the 70, 100 and 160 um bands and in medium and fast scan speeds. This branch of the PACS Photometer pipeline is the most used for cosmological observations and for point-source observations.The calibration of the flux loss due to the median removal applied by the PACS pipeline (high-pass filter) is done via dedicated simulations obtained by "polluting" real PACS timelines with fake sources at different flux levels. The effect of the data reduction parameter settings on the final map noise is done by using selected observations of blank fields with high data redundancy. We show that the running median removal can cause significant flux losses at any flux level. We analyse the advantages and disadvantages of several masking strategies and suggest that a mask based on putting circular patches on prior positions is the best solution to reduce the amount of flux loss. We provide a calibration of the point-source flux loss for several masking strategies in a large range of data reduction parameters, and as a function of the source flux. We also show that, for stacking analysis, the impact of the high-pass filtering effect is to reduce significantly the clustering effect. The analysis of the global noise and noise components of the PACS maps shows that the dominant parameter in determining the final noise is the high-pass filter width. We also provide simple fitting functions to build the error map from the coverage map and to estimate the cross-correlation correction factor in a representative portion of the data reduction parameter space.

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.

Context: Observations of transitional disks give us an understanding of the formation of planets and planetary systems such as our own. But care must be taken in the identification of such sources: the higher spatial resolution of the Herschel Space Observatory provides a new view on the origin of the far-infrared and sub-millimeter excesses observed. Aims: We review the nature of previously known transitional disks in the Chamaeleon I star-forming region with Herschel data. Methods: We analyze Herschel PACS and SPIRE images of the young star T54 together with ancillary images. We also analyze its spectral energy distribution and indications from optical and mid-infrared spectroscopy. Results: We detect extended emission in the PACS 70 \mu m image ~6" off source at a position angle of 196{\deg} from T54. The emission detected at longer wavelength (PACS 100, 160, SPIRE 250 and 350 \mu m) is also offset from the position of the star. This suggests that the excess observed in the far-infrared part of the SED is not fully associated with T54. Conclusions: Herschel images show that the far-infrared excess seen in T54 is not due to a transitional disk but to extended emission south-west of the source. The object still shows point-like and now downscaled excess at mid-infrared wavelengths, but its origin cannot be constrained without higher spatial resolution data. However, different indications point towards an evolved disk or extended unresolved emission close to the source.

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.

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.

Aims. We analyze the surroundings of HD 97300, one of two intermediate-mass stars in the Chamaeleon I star-forming region. The star is known to be surrounded by a conspicuous ring of polycyclic aromatic hydrocarbons (PAHs). Methods. We present infrared images taken with Herschel and Spitzer using 11 different broad-band filters between 3.6 um and 500 um. We compare the morphology of the emission using cuts along different position angles. We construct spectral energy distributions, which we compare to different dust models, and calculate dust temperatures. We also derive opacity maps and analyze the density structure of the environment of HD 97300. Results. We find that HD 97300 has no infrared excess at or below 24 um, confirming its zero-age main-sequence nature. The morphology of the ring is very similar between 3.6 um and 24 um. The emission at these wavelengths is dominated by either PAH features or PAH continuum. At longer wavelengths, only the northwestern part of the ring is visible. A fit to the 100-500 um observations suggests that the emission is due to relatively warm (~26 K) dust. The temperature gradually decreases with increasing distance from the ring. We find a general decrease in the density from north to south, and an approximate 10% density increase in the northeastern part of the ring. Conclusions. Our results are consistent with the theory that the ring around HD 97300 is essentially a bubble blown into the surrounding interstellar matter and heated by the star.

The C-rich AGB star IRC+10216 undergoes strong mass loss, and quasi-periodic density enhancements in the circumstellar matter have been reported. CO is ubiquitous in the CSE, while CCH emission comes from a spatially confined shell. With the IRAM 30m telescope and Herschel/HIFI, we recently detected unexpectedly strong emission from the CCH N=4-3, 6-5, 7-6, 8-7, and 9-8 transitions, challenging the available chemical and physical models. We aim to constrain the physical properties of IRC+10216's CSE, including the effect of episodic mass loss on the observed emission. In particular, we aim to determine the excitation region and conditions of CCH and to reconcile these with interferometric maps of the N=1-0 transition. Via radiative-transfer modelling, we provide a physical description of the CSE, constrained by the SED and a sample of 20 high-resolution and 29 low-resolution CO lines. We further present detailed radiative-transfer analysis of CCH. Assuming a distance of 150pc, the SED is modelled with a stellar luminosity of 11300Lsun and a dust-mass-loss rate of 4.0\times10^{-8}Msun/yr. Based on the analysis of 20 high resolution CO observations, an average gas-mass-loss rate for the last 1000yrs of 1.5\times10^{-5}Msun/yr is derived. This gives a gas-to-dust-mass ratio of 375, typical for an AGB star. The gas kinetic temperature throughout the CSE is described by 3 powerlaws: it goes as r^{-0.58} for r<9R*, as r^{-0.40} for 9<=r<=65R*, and as r^{-1.20} for r>65R*. This model successfully describes all 49 CO lines. We show the effect of wind-density enhancements on the CCH-abundance profile, and the good agreement of the model with the CCH N=1-0 transition and with the lines observed with the 30m telescope and HIFI. We report on the importance of radiative pumping to the vibrationally excited levels of CCH and the significant effect this has on the excitation of all levels of the CCH-molecule.

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.

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.

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 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.

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.

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.

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

The properties of the early-type binary Cyg OB2 #5 have been debated for many years and spectroscopic and photometric investigations yielded conflicting results. We have attempted to constrain the physical properties of the binary by collecting new optical and X-ray observations. We find that the orbital period of the system slowly changes though we are unable to discriminate between several possible explanations of this trend. The best fit solution of the continuum light curve reveals a contact configuration with the secondary star being significantly brighter and hotter on its leading side facing the primary. The mean temperature of the secondary star turns out to be only slightly lower than that of the primary, whilst the bolometric luminosity ratio is found to be 3.1. The solution of the light curve yields a distance of 925 +/- 25 pc much lower than the usually assumed distance of the Cyg OB2 association. Whilst we confirm the existence of episodes of higher X-ray fluxes, the data reveal no phase-locked modulation with the 6.6 day period of the eclipsing binary nor any clear relation between the X-ray flux and the 6.7 yr radio cycle. The bright region of the secondary star is probably heated by energy transfer in a common envelope in this contact binary system as well as by the collision with the primary's wind. The existence of a common photosphere probably also explains the odd mass-luminosity relation of the stars in this system. Most of the X-ray, non-thermal radio, and possibly gamma-ray emission of Cyg OB2 #5 is likely to arise from the interaction of the combined wind of the eclipsing binary with at least one additional star of this multiple system.

We analyse observations, spanning 15 years, dedicated to the extreme emission-line object Hen3-209. Our photometric data indicate that the luminosity of the star undergoes marked variations with a peak-to-peak amplitude of 0.65mag. These variations are recurrent, with a period of 16.093+-0.005d. The spectrum of Hen3-209 is peculiar with many different lines (HI, HeI, FeII,...) showing P Cygni profiles. The line profiles are apparently changing in harmony with the photometry. The spectrum also contains [OIII] lines that display a saddle profile topped by three peaks, with a maximum separation of about 600km/s. Hen3-209 is most likely an evolved luminous object suffering from mass ejection events and maybe belonging to a binary system.

We summarize the results of a multiwavelength observing campaign on the massive eclipsing binary CPD-41deg7742, another remarkable object in the young open cluster NGC 6231. Our campaign relies on high resolution echelle spectroscopy, narrow-band optical photometry, and XMM-Newton X-ray observations. Combined with the spectroscopic analysis, the light curve analysis provides a direct measurement of the masses and sizes of the system components. However, the most outstanding results come from the XMM-Newton observations. Our 180 ks campaign towards NGC 6231, and CPD-41deg7742, provides an unprecedented phase coverage of such a close early-type binary. The EPIC-MOS light curves almost fully cover the 2.4 day period of the system and the brightness of the object is sufficient to yield a time resolution as tight as 1 ks. The X-ray flux presents clear variations along the orbit, that we interpret as the signature of an unusual wind interaction. We indeed expect that, in this O+B system, the dominant primary wind crashes into the secondary surface, leading to a wind-photosphere interaction. As a strong support to our interpretation, we provide a geometrical model that associates an extra X-ray emission to the secondary inner surface. Though quite simple, the present model matches the main features of the X-ray light curve.

In the first paper of this series, we presented a detailed high-resolution spectroscopic study of CPD - 41degr7742, deriving for the first time an orbital solution for both components of the system. In this second paper, we focus on the analysis of the optical light curve and on recent XMM-Newton X-ray observations. In the optical, the system presents two eclipses, yielding an inclination i ~ 77degr. Combining the constraints from the photometry with the results of our previous work, we derive the absolute parameters of the system. We confirm that the two components of CPD - 41degr7742 are main sequence stars (O9 V + B1-1.5 V) with masses (M_1 ~ 18 Msol and M_2 ~ 10 Msol) and respective radii (R_1 ~ 7.5 Rsol and R_2 ~ 5.4 Rsol) close to the typical values expected for such stars. We also report an unprecedented set of X-ray observations that almost uniformly cover the 2.44-day orbital cycle. The X-ray emission from CPD - 41degr7742 is well described by a two-temperature thermal plasma model with energies close to 0.6 and 1.0 keV, thus slightly harder than typical early-type emission. The X-ray light curve shows clear signs of variability. The emission level is higher when the primary is in front of the secondary. During the high emission state, the system shows a drop of its X-ray emission that almost exactly matches the optical eclipse. We interpret the main features of the X-ray light curve as the signature of a wind-photosphere interaction, in which the overwhelming primary O9 star wind crashes into the secondary surface. Alternatively the light curve could result from a wind-wind interaction zone located near the secondary star surface. As a support to our interpretation, we provide a phenomenological geometric model that qualitatively reproduces the observed modulations of the X-ray emission.

We present Gemini-N GMOS and CFHT MOS spectroscopy of Wolf-Rayet candidates in the Local Group dwarf galaxy IC10 that were previously identified by Massey et al. and Royer et al. From the present spectroscopic survey, the WC/WN ratio for IC10 remains unusually high, given its low metallicity, although none of the WC9 stars suspected from narrow-band imaging are confirmed. Our spectroscopy confirms 9 newly discovered Wolf-Rayet candidates from Royer et al., whilst spectral types of 14 Wolf-Rayet stars previously observed by Massey & Armandroff are refined here. In total, there are 26 spectroscopically confirmed Wolf-Rayet stars in IC~10. All but one of the fourteen WC stars are WC4-6 stars, the exception being # 10 from Massey et al., a broad-lined, apparently single WC7 star. There are a total of eleven WN stars, which are predominantly early WN3-4 stars, but include a rare WN10 star, #8 from Royer et al. #5 from Massey et al. is newly identified as a transition WN/C star. Consequently, the WC/WN ratio for IC10 is 14/11=1.3, unusually high for a metal-poor galaxy. Re-evaluating recent photometric data of Massey & Holmes, we suggest that the true WC/WN ratio may not be as low as 0.3. Finally, we present ground-based finding charts for all confirmed WR stars, plus HST/WFPC2 charts for twelve cases.

Based on imaging and spectroscopy obtained at the ESO NTT telescope and using an efficient image analysis algorithm, we study the core of the LMC OB association LH13, particularly the two compact stellar clusters Sk-6641 and HNT in the HII, region N11C. We resolve Sk-6641 into 15 components and for the first time the HNT cluster into 70 stars, and derive photometry for the members. Moreover, from medium resolution spectroscopy we determine the spectral types for sixteen stars in N11C. We compare the color-magnitude diagrams of the clusters with that of the field stars and discuss the cluster ages. With an age of ~100 Myr, the HNT cluster appears significantly older than the very young (< 5 Myr) Sk-6641 starburst. We suggest that most of the `field' O-stars in the core of N11C have actually been ejected from Sk-6641 through dynamical interactions in the compact cluster. The properties of the Sk-6641 and HNT clusters suggest that we are viewing different star formation regions lying at different distances along the same line of sight.

This paper addresses the question of the observations to be performed with the Optical Monitor (OM) of the X-ray Multi-Mirror Satellite (XMM) under several aspects. First, we discuss XMM-OM's photometric system and its colour transformations towards the standard UBV system. Second, we establish a set of procedures to determine the temperature and the amount of interstellar absorption affecting the observed stars. Last, we address the possibility of isolating quasars in multidimensional colour diagrams based on the XMM-OM filter set.

We present here tests of a five-filter photometric system aimed at WR classification. In addition to the well-known easy separation between the WN and WC spectral types, these tests indicate interesting potentialities in the discrimination of subgroups among the WN and the WC which look well related to the classical subtypes. The proposed combinations of filters (or derived ones) bear enough discriminating power to satisfy some evolutionary studies in crowded fields where spectroscopic follow-up is not possible.