Stellar variability from pulsations and granulation presents a source of correlated noise that can impact the accuracy and precision of multi-band photometric transit observations of exoplanets. This can potentially cause biased measurements in the transmission or emission spectrum or underestimation of the final error bars on the spectrum. ARIEL is a future space telescope and instrument designed to perform a transit spectroscopic survey of a large sample of exoplanets. In this paper we perform simulations to assess the impact of stellar variability arising from pulsations and granulation on ARIEL observations of GJ 1214b and HD 209458b. We take into account the correlated nature of stellar noise, quantify it, and compare it to photon noise. In the range 1.95-7.8 \textmu m, stellar pulsation and granulation noise has insignificant impact compared to photon noise for both targets. In the visual range the contribution increases significantly but remains small in absolute terms and will have minimal impact on the transmission spectra of the targets studied. The impact of pulsation and granulation will be greatest for planets with low scale height atmospheres and long transit times around bright stars.

CubeSats are routinely used for low-cost photometry from space. Space-borne spectroscopy, however, is still the exclusive domain of much larger platforms. Key astrophysical questions in e.g. stellar physics and exoplanet research require uninterrupted spectral monitoring from space over weeks or months. Such monitoring of individual sources is unfortunately not affordable with these large platforms. With CUBESPEC we plan to offer the astronomical community a low-cost CubeSat solution for near-UV/optical/near-IR spectroscopy that enables this type of observations. CUBESPEC is a generic spectrograph that can be configured with minimal hardware changes to deliver both low resolution (R=100) with very large spectral coverage (200-1000nm), as well as high resolution (R=30,000) over a selected wavelength range. It is built around an off-axis Cassegrain telescope and a slit spectrograph with configurable dispersion elements. CUBESPEC will use a compact attitude determination and control system for coarse pointing of the entire spacecraft, supplemented with a fine-guidance system using a fast steering mirror to center the source on the spectrograph slit and to cancel out satellite jitter. An extremely compact optical design allows us to house this instrument in a 6U CubeSat with a volume of only 10x20x30cm$^{3}$, while preserving a maximized entrance pupil of ca. 9x19cm$^{2}$. In this contribution, we give an overview of the CUBESPEC project, discuss its most relevant science cases, and present the design of the instrument.

Measurements in the infrared wavelength domain allow us to assess directly the physical state and energy balance of cool matter in space, thus enabling the detailed study of the various processes that govern the formation and early evolution of stars and planetary systems in galaxies over cosmic time. Previous infrared missions, from IRAS to Herschel, have revealed a great deal about the obscured Universe, but sensitivity has been limited because up to now it has not been possible to fly a telescope that is both large and cold. SPICA is a mission concept aimed at taking the next step in mid- and far-infrared observational capability by combining a large and cold telescope with instruments employing state-of-the-art ultra-sensitive detectors. The mission concept foresees a 2.5-meter diameter telescope cooled to below 8 K. With cooling provided by mechanical coolers instead of depending on a limited cryogen supply, the mission lifetime can extend significantly beyond the required three years. SPICA offers instrumentation with spectral resolving powers ranging from R ~50 through 11000 in the 17-230 $\mu$m domain as well as R~28.000 spectroscopy between 12 and 18 $\mu$m. Additionally SPICA will provide efficient 30-37 $\mu$m broad band mapping, and polarimetric imaging in the 100-350 $\mu$m range. SPICA will provide unprecedented spectroscopic sensitivity of ~5 x $10^{-20}$ W/m$^2$ (5$\sigma$/1hr) - at least two orders of magnitude improvement over what has been attained to date. With this exceptional leap in performance, new domains in infrared astronomy will become accessible, allowing us, for example, to unravel definitively galaxy evolution and metal production over cosmic time, to study dust formation and evolution from very early epochs onwards, and to trace the formation history of planetary systems.

The Mid-Infrared Instrument (MIRI) Medium Resolution Spectrometer (MRS) is the only mid-IR Integral Field Spectrometer on board James Webb Space Telescope. The complexity of the MRS requires a very specialized pipeline, with some specific steps not present in other pipelines of JWST instruments, such as fringe corrections and wavelength offsets, with different algorithms for point source or extended source data. The MRS pipeline has also two different variants: the baseline pipeline, optimized for most foreseen science cases, and the optimal pipeline, where extra steps will be needed for specific science cases. This paper provides a comprehensive description of the MRS Calibration Pipeline from uncalibrated slope images to final scientific products, with brief descriptions of its algorithms, input and output data, and the accessory data and calibration data products necessary to run the pipeline.

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.

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

We describe the design and performance of the Medium Resolution Spectrometer (MRS) for the JWST-MIRI instrument. The MRS incorporates four coaxial spectral channels in a compact opto-mechanical layout that generates spectral images over fields of view up to 7.7 X 7.7 arcseconds in extent and at spectral resolving powers ranging from 1,300 to 3,700. Each channel includes an all-reflective integral field unit (IFU): an 'image slicer' that reformats the input field for presentation to a grating spectrometer. Two 1024 X 1024 focal plane arrays record the output spectral images with an instantaneous spectral coverage of approximately one third of the full wavelength range of each channel. The full 5 to 28.5 micron spectrum is then obtained by making three exposures using gratings and pass-band-determining filters that are selected using just two three-position mechanisms. The expected on-orbit optical performance is presented, based on testing of the MIRI Flight Model and including spectral and spatial coverage and resolution. The point spread function of the reconstructed images is shown to be diffraction limited and the optical transmission is shown to be consistent with the design expectations.

We present an estimate of the performance that will be achieved during on orbit operations of the JWST Mid Infrared Instrument, MIRI. The efficiency of the main imager and spectrometer systems in detecting photons from an astronomical target are presented, based on measurements at sub-system and instrument level testing, with the end-to-end transmission budget discussed in some detail. The brightest target fluxes that can be measured without saturating the detectors are provided. The sensitivity for long duration observations of faint sources is presented in terms of the target flux required to achieve a signal to noise ratio of 10 after a 10,000 second observation. The algorithms used in the sensitivity model are presented, including the understanding gained during testing of the MIRI Flight Model and flight-like detectors.

Photodissociation regions (PDRs) contain a large fraction of all of the interstellar matter in galaxies. Classical examples include the boundaries between ionized regions and molecular clouds in regions of massive star formation, marking the point where all of the photons energetic enough to ionize hydrogen have been absorbed. In this paper we determine the physical properties of the PDRs associated with the star forming regions IRAS 23133+6050 and S 106 and present them in the context of other Galactic PDRs associated with massive star forming regions. We employ Herschel PACS and SPIRE spectroscopic observations to construct a full 55-650 {\mu}m spectrum of each object from which we measure the PDR cooling lines, other fine- structure lines, CO lines and the total far-infrared flux. These measurements are then compared to standard PDR models. Subsequently detailed numerical PDR models are compared to these predictions, yielding additional insights into the dominant thermal processes in the PDRs and their structures. We find that the PDRs of each object are very similar, and can be characterized by a two-phase PDR model with a very dense, highly UV irradiated phase (n $\sim$ 10^6 cm^(-3), G$_0$ $\sim$ 10^5) interspersed within a lower density, weaker radiation field phase (n $\sim$ 10^4 cm^(-3), G$_0$ $\sim$ 10^4). We employed two different numerical models to investigate the data, firstly we used RADEX models to fit the peak of the $^{12}$CO ladder, which in conjunction with the properties derived yielded a temperature of around 300 K. Subsequent numerical modeling with a full PDR model revealed that the dense phase has a filling factor of around 0.6 in both objects. The shape of the $^{12}$CO ladder was consistent with these components with heating dominated by grain photoelectric heating. An extra excitation component for the highest J lines (J > 20) is required for S 106.

The Kepler space mission provided near-continuous and high-precision photometry of about 207,000 stars, which can be used for asteroseismology. However, for successful seismic modelling it is equally important to have accurate stellar physical parameters. Therefore, supplementary ground-based data are needed. We report the results of the analysis of high-resolution spectroscopic data of A- and F-type stars from the Kepler field, which were obtained with the HERMES spectrograph on the Mercator telescope. We determined spectral types, atmospheric parameters and chemical abundances for a sample of 117 stars. Hydrogen Balmer, Fe I, and Fe II lines were used to derive effective temperatures, surface gravities, and microturbulent velocities. We determined chemical abundances and projected rotational velocities using a spectrum synthesis technique. The atmospheric parameters obtained were compared with those from the Kepler Input Catalogue (KIC), confirming that the KIC effective temperatures are underestimated for A stars. Effective temperatures calculated by spectral energy distribution fitting are in good agreement with those determined from the spectral line analysis. The analysed sample comprises stars with approximately solar chemical abundances, as well as chemically peculiar stars of the Am, Ap, and Lambda Boo types. The distribution of the projected rotational velocity, Vsini, is typical for A and F stars and ranges from 8 to about 280 km/s, with a mean of 134 km/s.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Context: The 440 Myr old main-sequence A-star Fomalhaut is surrounded by an eccentric debris belt with sharp edges. This sort of a morphology is usually attributed to planetary perturbations, but the orbit of the only planetary candidate detected so far, Fomalhaut b, is too eccentric to efficiently shape the belt. Alternative models that could account for the morphology without invoking a planet are stellar encounters and gas-dust interactions. Aims: We aim to test the possibility of gas-dust interactions as the origin of the observed morphology by putting upper limits on the total gas content of the Fomalhaut belt. Methods: We derive upper limits on the CII 158 $\mu$m and OI 63 $\mu$m emission by using non-detections from the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel Space Observatory. Line fluxes are converted into total gas mass using the non-local thermodynamic equilibrium (non-LTE) code RADEX. We consider two different cases for the elemental abundances of the gas: solar abundances and abundances similar to those observed for the gas in the $\beta$ Pictoris debris disc. Results: The gas mass is shown to be below the millimetre dust mass by a factor of at least $\sim$3 (for solar abundances) respectively $\sim$300 (for $\beta$ Pic-like abundances). Conclusions: The lack of gas co-spatial with the dust implies that gas-dust interactions cannot efficiently shape the Fomalhaut debris belt. The morphology is therefore more likely due to a yet unseen planet (Fomalhaut c) or stellar encounters.

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

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 aimed to measure the H2O and dust production rates in C/2006 W3 (Christensen) with the Herschel Space Observatory at a heliocentric distance of ~ 5 AU. We have searched for emission in the H2O and NH3 ground-state rotational transitions at 557 GHz and 572 GHz, simultaneously, with HIFI onboard Herschel on UT 1.5 September 2010. Photometric observations of the dust coma in the 70 and 160 {\mu}m channels were acquired with the PACS instrument on UT 26.5 August 2010. A tentative 4-{\sigma} H2O line emission feature was found in the spectra obtained with the HIFI wide-band and high-resolution spectrometers, from which we derive a water production rate of $2.0(5) \times 10^{27}$ molec. s$^{-1}$. A 3-{\sigma} upper limit for the ammonia production rate of <$1.5 \times 10^{27}$ molec. s$^{-1}$ is obtained taking into account the contribution from all hyperfine components. The blueshift of the water line detected by HIFI suggests preferential emission from the subsolar point. However, it is also possible that water sublimation occurs in small ice-bearing grains that are emitted from an active region on the nucleus surface at a speed of ~ 0.2 km s$^{-1}$. The dust thermal emission was detected in the 70 and 160 {\mu}m filters, with a more extended emission in the blue channel. The dust production rates, obtained for a dust size distribution index that explains the fluxes at the photocenters of the PACS images, lie in the range 70-110 kg s$^{-1}$. Scaling the CO production rate measured post-perihelion at 3.20 and 3.32 AU, these values correspond to a dust-to-gas production rate ratio in the range 0.3-0.4. The dust production rates derived in August 2010 are roughly one order of magnitude lower than in September 2009, suggesting that the dust-to-gas production rate ratio remained approximately constant during the period when the activity became increasingly dominated by CO outgassing.

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.

Transit-spectroscopy of exoplanets is one of the key observational techniques to characterize the extrasolar planet and its atmosphere. The observational challenges of these measurements require dedicated instrumentation and only the space environment allows an undisturbed access to earth-like atmospheric features such as water or carbon-dioxide. Therefore, several exoplanet-specific space missions are currently being studied. One of them is EChO, the Exoplanet Characterization Observatory, which is part of ESA's Cosmic Vision 2015-2025 program, and which is one of four candidates for the M3 launch slot in 2024. In this paper we present the results of our assessment study of the EChO spectrometer, the only science instrument onboard this spacecraft. The instrument is a multi-channel all-reflective dispersive spectrometer, covering the wavelength range from 400 nm to 16 microns simultaneously with a moderately low spectral resolution. We illustrate how the key technical challenge of the EChO mission - the high photometric stability - influences the choice of spectrometer concept and drives fundamentally the instrument design. First performance evaluations underline the fitness of the elaborated design solution for the needs of the EChO mission.

Herschel-PACS measurements of the rotational R(0) and R(1) HD lines in the atmospheres of Uranus and Neptune are analyzed in order to derive a D/H ratio with improved precision for both planets. The derivation of the D/H ratio includes also previous measurements of the R(2) line by the Short Wavelength Spectrometer on board the Infrared Space Observatory (ISO). The available spectroscopic line information of the three rotational transitions is discussed and applied in the radiative transfer calculations. The best simultaneous fit of all three lines requires only a minor departure from the Spitzer temperature profile of Uranus and a departure limited to 2K from the Voyager temperature profile of Neptune (both around the tropopause). The resulting and remarkably similar D/H ratios for Uranus and Neptune are found to be (4.4$\pm$0.4)$\times10^{-5}$ and (4.1$\pm$0.4)$\times10^{-5}$ respectively. Although the deuterium enrichment in both atmospheres compared to the protosolar value is confirmed, it is found to be lower compared to previous analysis. Using the interior models of Podolak et al. (1995), Helled et al. (2011) and Nettelmann et al. (2013), and assuming that complete mixing of the atmosphere and interior occured during the planets history, we derive a D/H in protoplanetary ices between (5.75--7.0)$\times10^{-5}$ for Uranus and between (5.1--7.7)$\times10^{-5}$ for Neptune. Conversely, adopting a cometary D/H for the protoplanetary ices between (15-30)$\times10^{-5}$, we constrain the interior models of both planets to have an ice mass fraction of 14-32%, i.e. that the two planets are rock-dominated.

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 report a multisite photometric campaign for the Beta Cep stars V2052 Oph and V986 Oph. 670 hours of high-quality differential photoelectric Stromgren, Johnson and Geneva time-series photometry were obtained with eight telescopes on five continents during 182 nights. Frequency analyses of the V2052 Oph data enabled the detection of three pulsation frequencies, the first harmonic of the strongest signal, and the rotation frequency with its first harmonic. Pulsational mode identification from analysing the colour amplitude ratios confirms the dominant mode as being radial, whereas the other two oscillations are most likely l=4. Combining seismic constraints on the inclination of the rotation axis with published magnetic field analyses we conclude that the radial mode must be the fundamental. The rotational light modulation is in phase with published spectroscopic variability, and consistent with an oblique rotator for which both magnetic poles pass through the line of sight. The inclination of the rotation axis is 54o <i< 58o and the magnetic obliquity 58o <beta< 66o. The possibility that V2052 Oph has a magnetically confined wind is discussed. The photometric amplitudes of the single oscillation of V986 Oph are most consistent with an l=3 mode, but this identification is uncertain. Additional intrinsic, apparently temporally incoherent, light variations of V986 Oph are reported. Different interpretations thereof cannot be distinguished at this point, but this kind of variability appears to be present in many OB stars. The prospects of obtaining asteroseismic information for more rapidly rotating Beta Cep stars, which appear to prefer modes of higher l, are briefly discussed.

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.

Herschel images in six photometric bands show the thermal emission of the debris disk surrounding beta Pic. In the three PACS bands at 70 micron, 100 micron and 160 micron and in the 250 micron SPIRE band, the disk is well-resolved, and additional photometry is available in the SPIRE bands at 350 micron and 500 micron, where the disk is only marginally resolved. The SPIRE maps reveal a blob to the southwest of beta Pic, coinciding with submillimetre detection of excess emission in the disk. We investigated the nature of this blob. Our comparison of the colours, spectral energy distribution and size of the blob, the disk and the background sources shows that the blob is most likely a background source with a redshift between z =1.0 and z = 1.6.

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.

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.

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.

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.

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

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

In an effort to simultaneously study the gas and dust components of the disc surrounding the young Herbig Ae star HD 169142, we present far-IR observations obtained with the PACS instrument onboard the Herschel Space Observatory. This work is part of the Open Time Key Project GASPS, which is aimed at studying the evolution of protoplanetary discs. To constrain the gas properties in the outer disc, we observed the star at several key gas-lines, including [OI] 63.2 and 145.5 micron, [CII] 157.7 micron, CO 72.8 and 90.2 micron, and o-H2O 78.7 and 179.5 micron. We only detect the [OI] 63.2 micron line in our spectra, and derive upper limits for the other lines. We complement our data set with PACS photometry and 12/13CO data obtained with the Submillimeter Array. Furthermore, we derive accurate stellar parameters from optical spectra and UV to mm photometry. We model the dust continuum with the 3D radiative transfer code MCFOST and use this model as an input to analyse the gas lines with the thermo-chemical code ProDiMo. Our dataset is consistent with a simple model in which the gas and dust are well-mixed in a disc with a continuous structure between 20 and 200 AU, but this is not a unique solution. Our modelling effort allows us to constrain the gas-to-dust mass ratio as well as the relative abundance of the PAHs in the disc by simultaneously fitting the lines of several species that originate in different regions. Our results are inconsistent with a gas-poor disc with a large UV excess; a gas mass of 5.0 +/- 2.0 times 10^(-3) Msun is still present in this disc, in agreement with earlier CO observations.

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.

Context - Circumstellar discs are ubiquitous around young stars, but rapidly dissipate their gas and dust on timescales of a few Myr. The Herschel space observatory allows for the study of the warm disc atmosphere, using far-infrared spectroscopy to measure gas content and excitation conditions, and far-IR photometry to constrain the dust distribution. Aims - We aim to detect and characterize the gas content of circumstellar discs in four targets as part of the Herschel science demonstration phase. Methods - We carried out sensitive medium resolution spectroscopy and high sensitivity photometry at lambda ~60-190 micron using the Photodetector Array Camera and Spectrometer instrument on the Herschel space observatory. Results - We detect [OI] 63 micron emission from the young stars HD 169142, TW Hydrae, and RECX 15, but not HD 181327. No other lines, including [CII] 158 and [OI] 145, are significantly detected. All four stars are detected in photometry at 70 and 160 micron. Extensive models are presented in associated papers.

Planets are formed in disks around young stars. With an age of ~10 Myr, TW Hya is one of the nearest T Tauri stars that is still surrounded by a relatively massive disk. In addition a large number of molecules has been found in the TW Hya disk, making TW Hya the perfect test case in a large survey of disks with Herschel-PACS to directly study their gaseous component. We aim to constrain the gas and dust mass of the circumstellar disk around TW Hya. We observed the fine-structure lines of [OI] and [CII] as part of the Open-time large program GASPS. We complement this with continuum data and ground-based 12CO 3-2 and 13CO 3-2 observations. We simultaneously model the continuum and the line fluxes with the 3D Monte-Carlo code MCFOST and the thermo-chemical code ProDiMo to derive the gas and dust masses. We detect the [OI] line at 63 micron. The other lines that were observed, [OI] at 145 micron and [CII] at 157 micron, are not detected. No extended emission has been found. Preliminary modeling of the photometric and line data assuming [12CO]/[13CO]=69 suggests a dust mass for grains with radius < 1 mm of ~1.9 times 10^-4 Msun (total solid mass of 3 times 10^-3 Msun) and a gas mass of (0.5--5) times 10^-3 Msun. The gas-to-dust mass may be lower than the standard interstellar value of 100.

The Herschel GASPS Key Program is a survey of the gas phase of protoplanetary discs, targeting 240 objects which cover a large range of ages, spectral types, and disc properties. To interpret this large quantity of data and initiate self-consistent analyses of the gas and dust properties of protoplanetary discs, we have combined the capabilities of the radiative transfer code MCFOST with the gas thermal balance and chemistry code ProDiMo to compute a grid of 300 000 disc models (DENT). We present a comparison of the first Herschel/GASPS line and continuum data with the predictions from the DENT grid of models. Our objective is to test some of the main trends already identified in the DENT grid, as well as to define better empirical diagnostics to estimate the total gas mass of protoplanetary discs. Photospheric UV radiation appears to be the dominant gas-heating mechanism for Herbig stars, whereas UV excess and/or X-rays emission dominates for T Tauri stars. The DENT grid reveals the complexity in the analysis of far-IR lines and the difficulty to invert these observations into physical quantities. The combination of Herschel line observations with continuum data and/or with rotational lines in the (sub-)millimetre regime, in particular CO lines, is required for a detailed characterisation of the physical and chemical properties of circumstellar discs.

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.

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.

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.

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.

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.

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.

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.

Gamma Doradus stars are excellent targets for asteroseismology since the gravity modes present in these stars probe the deep stellar interiors. Mode identification will improve the knowledge of these stars considerably. A selected group of Gamma Doradus stars and some candidates were observed with the Mercator telescope to find and/or confirm the periodicities in the light variations and to derive reliable amplitude ratios in different pass bands. A frequency analysis was performed on all new data obtained in the Geneva photometric system. In order to have more reliable and accurate frequencies, the new data were combined with similar data from the literature and with Hipparcos observations. A set of frequencies that minimized the the residuals in a harmonic fit was searched for while allowing means and amplitudes to vary from one observation set to another. Frequencies and amplitudes in the photometric passbands of the Geneva system are given for 21 Gamma Doradus stars. We report the discovery of HD 74504 as a newly found Gamma Doradus star. Our study provides the first extensive multicolour database for the understanding of gravity modes in F-type stars.

The line profile variability and photometric variability of the O9.5 Vp star HD93521 are examined in order to establish the properties of the non-radial pulsations in this star. Fourier techniques are used to characterize the modulations of the He I 5876, 6678 and H-alpha lines in several spectroscopic time series and to search for variations in a photometric time series. Our spectroscopic data confirm the existence of two periods of 1.75 and 2.89 hr. The line profiles, especially those affected by emission wings, exhibit also modulations on longer time scales, but these are epoch-dependent and change from line to line. Unlike previous claims, we find no unambiguous signature of the rotational period in our data, nor of a third pulsation period (corresponding to a frequency of 2.66 day$^{-1}$). HD 93521 very likely exhibits non-radial pulsations with periods of 1.75 and 2.89 hr with $l \simeq 8 \pm 1$ and $l \simeq 4 \pm 1$ respectively. No significant signal is found in the first harmonics of these two periods. The 2.89 hr mode is seen at all epochs and in all lines investigated, while the visibility of the 1.75 hr mode is clearly epoch dependent. Whilst light variations are detected, their connection to these periodicities is not straightforward.

The fast classification of new variable stars is an important step in making them available for further research. Selection of science targets from large databases is much more efficient if they have been classified first. Defining the classes in terms of physical parameters is also important to get an unbiased statistical view on the variability mechanisms and the borders of instability strips. Our goal is twofold: provide an overview of the stellar variability classes that are presently known, in terms of some relevant stellar parameters; use the class descriptions obtained as the basis for an automated `supervised classification' of large databases. Such automated classification will compare and assign new objects to a set of pre-defined variability training classes. For every variability class, a literature search was performed to find as many well-known member stars as possible, or a considerable subset if too many were present. Next, we searched on-line and private databases for their light curves in the visible band and performed period analysis and harmonic fitting. The derived light curve parameters are used to describe the classes and define the training classifiers. We compared the performance of different classifiers in terms of percentage of correct identification, of confusion among classes and of computation time. We describe how well the classes can be separated using the proposed set of parameters and how future improvements can be made, based on new large databases such as the light curves to be assembled by the CoRoT and Kepler space missions.

We present results of a spectroscopic study of 37 southern (candidate) gamma Doradus stars based on echelle spectra. The observed spectra were cross-correlated with the standard template spectrum of an F0-type star for an easier detection of binary and intrinsic variations. We identified 15 objects as spectroscopic binaries, including 7 new ones, and another 3 objects are binary suspects. At least 12 objects show composite spectra. We could determine the orbital parameters for 9 binaries, of which 4 turn out to be ellipsoidal variables. For 6 binaries, we estimated the expected time-base of the orbital variations. Clear profile variations are observed for 17 objects, pointing towards stellar pulsation. For 8 of them, we have evidence that the main spectroscopic and photometric periods coincide. Our results, in combination with prior knowledge from the literature, lead to the classification of 10 objects as new bona-fide gamma Doradus stars, 1 object as new bona-fide delta Scuti star, and 8 objects as constant stars. Finally, we determined the projected rotational velocity with two independent methods. The resulting vsini values range from 3 to 135 km/s. For the bona-fide gamma Doradus stars, the majority has vsini below 60 km/s.

In this paper we study the extended atmosphere of the late-type supergiant Alpha Orionis. Infrared spectroscopy of red supergiants reveals strong molecular bands, some of which do not originate in the photosphere but in a cooler layer of molecular material above it. Lately, these layers have been spatially resolved by near and mid-IR interferometry. In this paper, we try to reconcile the IR interferometric and ISO-SWS spectroscopic results on Alpha Orionis with a thorough modelling of the photosphere, molecular layer(s) and dust shell. From the ISO and near-IR interferometric observations, we find that Alpha Orionis has only a very low density water layer close above the photosphere. However, mid-IR interferometric observations and a narrow-slit N-band spectrum suggest much larger extra-photospheric opacity close to the photosphere at those wavelengths, even when taking into account the detached dust shell. We argue that this cannot be due to the water layer, and that another source of mid-IR opacity must be present. We show that this opacity source is probably neither molecular nor chromospheric. Rather, we present amorphous alumina (Al2O3) as the best candidate and discuss this hypothesis in the framework of dust-condensation scenarios.

presented. The 2.38 -- 4.08 micron wavelength-range of band 1 of ISO-SWS (Short-Wavelength Spectrometers on board of the Infrared Space Observatory) in which many different molecules -- with their own dependence on each of the stellar parameters -- are absorbing, enables us to estimate the effective temperature, the gravity, the microturbulence, the metallicity, the CNO-abundances, the 12C/13C-ratio and the angular diameter from the ISO-SWS data. Using the Hipparcos' parallax, the radius, luminosity and gravity-inferred mass are derived. The stellar parameters obtained are in good agreement with other published values, though also some discrepancies with values deduced by other authors are noted. For a few stars (Delta Dra, Xi Dra, Alpha Tuc, H Sco and Alpha Cet) some parameters -- e.g. the CNO-abundances -- are derived for the first time. By examining the correspondence between different ISO-SWS observations of the same object and between the ISO-SWS data and the corresponding synthetic spectrum, it is shown that the relative accuracy of ISO-SWS in band 1 (2.38 -- 4.0 micron) is better than 2% for these high-flux sources. The high level of correspondence between observations and theoretical predictions, together with a confrontation of the estimated Teff(ISO) value with Teff-values derived from colours -- which demonstrates the consistency between V-K, BC(K), Teff and the angular diameter derived from optical or IR data -- proves that both the used MARCS models to derive the stellar quantities and the flux calibration of the ISO-SWS detectors have reached a high level of reliability.

Vega, Sirius, Beta Leo, Alpha Car and Alpha Cen A belong to a sample of twenty stellar sources used for the calibration of the detectors of the Short-Wavelength Spectrometer on board the Infrared Space Observatory (ISO-SWS). While general problems with the calibration and with the theoretical modelling of these stars are reported in Decin et al. (2002), each of these stars is discussed individually in this paper. As demonstrated in Decin et al. (2002), it is not possible to deduce the effective temperature, the gravity and the chemical composition from the ISO-SWS spectra of these stars. But since ISO-SWS is absolutely calibrated, the angular diameter of these stellar sources can be deduced from their ISO-SWS spectra, which consequently yields the stellar radius (R), the gravity-inferred mass (M) and the luminosity (L) for these stars. These deduced parameters are confronted with other published values and the goodness-of-fit between observed ISO-SWS data and the corresponding synthetic spectrum is discussed.

The fine calibration of the ISO-SWS detectors (Infrared Space Observatory - Short Wavelength Spectrometer) has proven to be a delicate problem. We therefore present a detailed spectroscopic study in the 2.38 -- 12 micron wavelength range of a sample of 16 A0 -- M2 stars used for the calibration of ISO-SWS. By investigating the discrepancies between the ISO-SWS data of these sources, the theoretical predictions of their spectra, the high-resolution FTS-KP (Kitt Peak) spectrum of Alpha Boo and the solar FTS-ATMOS (Atmospheric Trace Molecule Spectroscopy) spectrum, both calibration problems and problems in computing the theoretical models and the synthetic spectra are revealed. The underlying reasons for these problems are sought for and the impact on the further calibration of ISO-SWS and on the theoretical modelling is discussed extensively.

We present an atlas of spectra of O- and B-type stars, obtained with the Short Wavelength Spectrometer (SWS) during the Post-Helium program of the Infrared Space Observatory (ISO). This program is aimed at extending the Morgan & Keenan classification scheme into the near-infrared. Later type stars will be discussed in a seperate publication. The observations consist of 57 SWS Post-Helium spectra from 2.4 to 4.1 microns, supplemented with 10 spectra acquired during the nominal mission with a similar observational setting. For B-type stars, this sample provides ample spectral converage in terms of subtype and luminosity class. For O-type stars,the ISO sample is coarse and therefore is complemented with 8 UKIRT L'-band observations. In terms of the presence of diagnostic lines, the L'-band is likely the most promising of the near-infrared atmospheric windows for the study of the physical properties of B stars. Specifically, this wavelength interval contains the Brackett alpha, Pfund gamma, and other Pfund lines which are probes of spectral type, luminosity class and mass loss. Here, we present simple empirical methods based on the lines present in the 2.4 to 4.1 microns interval that allow the determination of: the spectral type of B dwarfs and giants to within two subtypes; the luminosity class of B stars to within two classes; the mass-loss rate of O stars and B supergiants to within 0.25 dex.

We present a detailed spectroscopic study of a sample of bright, mostly cool, stars observed with the Short-Wavelength Spectrometer (SWS) on board the Infrared Space Observatory (ISO), which enables the accurate determination of the stellar parameters of the cool giants, but also serves as a critical review of the ISO-SWS calibration.

A detailed spectroscopic study of the ISO-SWS data of the red giant Alpha Tau is presented, which enables not only the accurate determination of the stellar parameters of Alpha Tau, but also serves as a critical review of the ISO-SWS calibration. This study is situated in a broader context of an iterative process in which both accurate observations of stellar templates and cool star atmosphere models are involved to improve the ISO-SWS calibration process as well as the theoretical modelling of stellar atmospheres. Therefore a sample of cool stars, covering the whole A0 -- M8 spectral classification, has been observed in order to disentangle calibration problems and problems in generating the theoretical models and corresponding synthetic spectrum. By using stellar parameters found in the literature large discrepancies were seen between the ISO-SWS data and the generated synthetic spectrum of Alpha Tau. A study of the influence of various stellar parameters on the theoretical models and synthetic spectra, in conjunction with the Kolmogorov-Smirnov test to evaluate objectively the goodness-of-fit, enables us to pin down the stellar parameters with a high accuracy: Teff = 3850 +/- 70 K, log g = 1.50 +/- 0.15, M = 2.3 +/- 0.8 Msun, z = -0.15 +/- 0.20 dex, microturbulence = 1.7 +/- 0.3 km/s, 12C/13C= 10 +/- 1, abundance of C = 8.35 +/- 0.20 dex, abundance of N= 8.35 +/- 0.25 dex, abundance of O = 8.83 +/- 0.15 dex and the angular diameter is 20.77 +/- 0.83 mas. These atmospheric parameters were then compared with the results provided by other authors using other methods and/or spectra.