• 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.
  • The predicted orbital period histogram of an sdB population is bimodal with a peak at short (< 10 days) and long (> 250 days) periods. Observationally, there are many short-period sdB systems known, but only very few long-period sdB binaries are identified. As these predictions are based on poorly understood binary interaction processes, it is of prime importance to confront the predictions to observational data. In this contribution we aim to determine the absolute dimensions of the long-period sdB+MS binary system PG1104+243. High-resolution spectroscopy time-series were obtained with HERMES at the Mercator telescope at La Palma, and analyzed to obtain radial velocities of both components. Photometry from the literature was used to construct the spectral energy distribution (SED) of the binary. Atmosphere models were used to fit this SED and determine the surface gravity and temperature of both components. The gravitational redshift provided an independent confirmation of the surface gravity of the sdB component. An orbital period of 753 +- 3 d and a mass ratio of q = 0.637 +- 0.015 were found from the RV-curves. The sdB component has an effective temperature of Teff = 33500 +- 1200 K and a surface gravity of logg = 5.84 +- 0.08 dex, while the cool companion is found to be a G-type star with Teff = 5930 +- 160 K and logg = 4.29 +- 0.05 dex. Assuming a canonical mass of Msdb = 0.47 Msun, the MS component has a mass of 0.74 +- 0.07 Msun, and its Teff corresponds to what is expected for a terminal age main-sequence star with sub-solar metalicity. PG1104+243 is the first long-period sdB binary in which accurate physical parameters of both components could be determined, and the first sdB binary in which the gravitational redshift is measured. Furthermore, PG1104+243 is the first sdB+MS system that shows consistent evidence for being formed through stable Roche-lobe overflow.
  • 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.
  • We present a light curve analysis and radial velocity study of KOI-74, an eclipsing A star + white dwarf binary with a 5.2 day orbit. Aside from new spectroscopy covering the orbit of the system, we used 212 days of publicly available Kepler observations and present the first complete light curve fitting to these data, modelling the eclipses and transits, ellipsoidal modulation, reflection, and Doppler beaming. Markov Chain Monte Carlo simulations are used to determine the system parameters and uncertainty estimates. Our results are in agreement with earlier studies, except that we find an inclination of 87.0 \pm 0.4\degree, which is significantly lower than the previously published value. We find that the mass ratio derived from the radial velocity amplitude (q=0.104 \pm 0.004) disagrees with that derived from the ellipsoidal modulation (q=0.052 \pm 0.004} assuming corotation). This was found before, but with our smaller inclination, the discrepancy is even larger than previously reported. Accounting for the rapid rotation of the A-star is found to increase the discrepancy even further by lowering the mass ratio to q=0.047 \pm 0.004. These results indicate that one has to be extremely careful in using the amplitude of an ellipsoidal modulation signal in a close binary to determine the mass ratio, when a proof of corotation is not firmly established. The radial velocities that can be inferred from the detected Doppler beaming in the light curve are found to be in agreement with our spectroscopic radial velocity determination. We also report the first measurement of R{\o}mer delay in a light curve of a compact binary. This delay amounts to -56 \pm 17 s and is consistent with the mass ratio derived from the radial velocity amplitude. The firm establishment of this mass ratio at q=0.104 \pm 0.004 leaves little doubt that the companion of KOI-74 is a low mass white dwarf.
  • Recent studies have shown that stellar chromospheric activity, and its effect on convective energy transport in the envelope, is most likely the cause of significant radius and temperature discrepancies between theoretical evolution models and observations. We aim to determine absolute dimensions and abundances for the solar-type detached eclipsing binary EF Aqr, and to perform a detailed comparison with results from recent stellar evolutionary models. uvby-beta standard photometry was obtained with the Stromgren Automatic Telescope. The broadening function formalism was applied on spectra observed with HERMES at the Mercator telescope in La Palma, to obtain radial velocity curves. Masses and radii with a precision of 0.6% and 1.0% respectively have been established for both components of EF Aqr. The active 0.956 M_sol secondary shows star spots and strong Ca II H and K emission lines. The 1.224 M_sol primary shows signs of activity as well, but at a lower level. An [Fe/H] abundance of 0.00+-0.10 is derived with similar abundances for Si, Ca, Sc, Ti, V, Cr, Co, and Ni. Solar calibrated evolutionary models such as Yonsei-Yale, Victoria-Regina and BaSTI isochrones and evolutionary tracks are unable to reproduce EF Aqr, especially for the secondary, which is 9% larger and 400 K cooler than predicted. Models adopting significantly lower mixing length parameters l/H_p remove these discrepancies, as seen in other solar type binaries. For the observed metallicity, Granada models with a mixing length of l/H_p=1.30 (primary) and 1.05 (secondary) reproduce both components at a common age of 1.5+-0.6 Gyr. Observations of EF Aqr suggests that magnetic activity, and its effect on envelope convection, is likely to be the cause of discrepancies in both radius and temperature, which can be removed by adjusting the mixing length parameter of the models downwards.
  • S-type AGB stars are thought to be in the transitional phase between M-type and C-type AGB stars. Because of their peculiar chemical composition, one may expect a strong influence of the stellar C/O ratio on the molecular chemistry and the mineralogy of the circumstellar dust. In this paper, we present a large sample of 87 intrinsic galactic S-type AGB stars, observed at infrared wavelengths with the Spitzer Space Telescope, and supplemented with ground-based optical data. On the one hand, we derive the stellar parameters from the optical spectroscopy and photometry, using a grid of model atmospheres. On the other, we decompose the infrared spectra to quantify the flux-contributions from the different dust species. Finally, we compare the independently determined stellar parameters and dust properties. For the stars without significant dust emission, we detect a strict relation between the presence of SiS absorption in the Spitzer spectra and the C/O ratio of the stellar atmosphere. These absorption bands can thus be used as an additional diagnostic for the C/O ratio. For stars with significant dust emission, we define three groups, based on the relative contribution of certain dust species to the infrared flux. We find a strong link between group-membership and C/O ratio. We show that these groups can be explained by assuming that the dust-condensation can be cut short before silicates are produced, while the remaining free atoms and molecules can then form the observed magnesium sulfides or the carriers of the unidentified 13 and 20 micron features. Finally, we present the detection of emission features attributed to molecules and dust characteristic to C-type stars, such as molecular SiS, hydrocarbons and magnesium sulfide grains. We show that we often detect magnesium sulfides together with molecular SiS and we propose that it is formed by a reaction of SiS molecules with Mg.
  • 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.
  • SuWt2 is a planetary nebula (PN) consisting of a bright ionized thin ring seen nearly edge-on. It has a bright (V=12) central star, too cool to ionize the PN, which we discovered to be an eclipsing binary. A spectrum from IUE did not reveal a UV source. We present extensive ground-based photometry and spectroscopy of the central binary collected over the ensuing two decades, resulting in the determination that the orbital period of the eclipsing pair is 4.9 d, and consists of two nearly identical A1 V stars, each of mass ~2.7 M_sun. The physical parameters of the A stars, combined with evolutionary tracks, show that both are in the short-lived "blue-hook" evolutionary phase that occurs between the main sequence and the Hertzsprung gap, and that the age of the system is about 520 Myr. One puzzle is that the stars' rotational velocities are different from each other, and considerably slower than synchronous with the orbital period. It is possible that the center-of-mass velocity of the eclipsing pair is varying with time, suggesting that there is an unseen third orbiting body in the system. We propose a scenario in which the system began as a hierarchical triple, consisting of a ~2.9 M_sun star orbiting the close pair of A stars. Upon reaching the AGB stage, the primary engulfed the pair into a common envelope, leading to a rapid contraction of the orbit and catastrophic ejection of the envelope into the orbital plane. In this picture, the exposed core of the initial primary is now a white dwarf of ~0.7 M_sun, orbiting the eclipsing pair, which has already cooled below the detectability possible by IUE at our derived distance of 2.3 kpc and a reddening of E(B-V)=0.40. The SuWt2 system may be destined to perish as a Type Ia supernova. (Abridged)
  • 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 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.
  • 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.
  • The aim of this work is to make available unpublished non-Fe+ emission line fluxes from optical spectra of the symbiotic nova RR Tel which were taken in 2000, and to compare them with fluxes of the same lines from spectra taken in 1996. After leaving out blends and misidentifications, as well as the unreliable far-red and violet lines, we present the log (F2000/F1996) flux ratios for identified non-Fe+ lines. Mean values of log (F2000/F1996) for different ionization potential ranges of the ions producing the lines are shown separately for the permitted and forbidden lines. All means show fading, which is larger in the lowest range of ionization potential. Provisional interpretations are suggested. We also measured the values of FWHM in 2000; the previously known decrease with time of FWHM of lines due to the same ion has continued.
  • Powerful radio galaxies and radio-loud quasars at high redshifts are frequently associated with extended emission-line regions (EELRs). Here we investigate the [O II] EELR around the quasar 3C 196 at z=0.871 using integral field spectroscopy. We also detect extended [Ne II] emission at a distance of about 30 kpc from the core. The emission is aligned with the radio hot spots and shows a redshifted and a blueshifted component with a velocity difference of ~800 km s^-1. The alignment effect and large velocities support the hypothesis that the EELR is caused by a jet-cloud interaction, which is furthermore indicated by the presence of a pronounced bend in the radio emission at the location of the radio hot spots. We also report observations of two other systems which do not show as clear indications of interactions. We find a weaker alignment of an [O II] EELR from the z=0.927 quasar 3C 336, while no EELR is found around the core-dominated quasar OI 363 at z=0.63.
  • Infrared and optical photometric and spectroscopic observations of the symbiotic nova RR Tel are used to study the effects and properties of dust in symbiotic binaries containing a cool Mira component, as well as showing "obscuration events" of increased absorption, which are typical for such Miras. A set of photometric observations of the symbiotic nova RR Tel in different wavelength bands - visual from 1949 to 2002 and near-infrared (JHKL) from 1975 to 2002 - are presented. The variability due to the normal Mira pulsation was removed from the JHKL data, which were then compared with the AAVSO visual light curve. The changes of the Fe II emission line fluxes during the 1996-2000 obscuration episode were studied in the optical spectra taken with the Anglo-Australian telescope. We discuss the three periods during which the Mira component was heavily obscured by dust as observed in the different wavelength bands. A change in the correlations of J with other infrared magnitudes was observed with the colour becoming redder after JD2446000. Generally, J-K was comparable, while K-L was larger than typical values for single Miras. A distance estimate of 2.5 kpc, based on the IR data, is given. A larger flux decrease for the permitted than for the forbidden Fe II lines, during the obscuration episode studied, has been found. There is no evidence for other correlations with line properties, in particular with wavelength, which suggests obscuration due to separate optically thick clouds in the outer layers.