We report the first spectroscopic observations of a background star seen through the region between the ionization front and the dissociation front of the nebula IC 63. This photodissociation region (PDR) exhibits intense extended red emission (ERE) attributed to fluorescence by large molecules/ions. We detected strong diffuse interstellar bands (DIB) in the stellar spectrum, including an exceptionally strong and broad DIB at $\lambda$4428. The detection of strong DIBs in association with ERE could be consistent with the suggestion that the carriers of DIBs and ERE are identical. The likely ERE process is recurrent fluorescence, enabled by inverse internal conversions from highly excited vibrational levels of the ground state to low-lying electronic states with subsequent transitions to ground. This provides a path to rapid radiative cooling for molecules/molecular ions, greatly enhancing their ability to survive in a strongly irradiated environment. The ratio of the equivalent widths (EW) of DIBs $\lambda$5797 and $\lambda$5780 in IC 63 is the same as that observed in the low-density interstellar medium with UV interstellar radiation fields (ISRF) weaker by at least two orders of magnitude. This falsifies suggestions that the ratio of these two DIBs can serve as a measure of the UV strength of the ISRF. Observations of the nebular spectrum of the PDR of IC 63 at locations immediately adjacent to where DIBs were detected failed to reveal any presence of sharp emission features seen in the spectrum of the Red Rectangle nebula. This casts doubts upon proposals that the carriers of these features are the same as those of DIBs seen at slightly shorter wavelengths.

(Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada.

The interstellar medium mediates galactic evolution as the reservoir of material for future star formation and the repository of energy and matter output by stellar processes. Canadians have played leading roles in ISM science for decades. The Canadian Galactic Plane Survey identified a wealth of small-scale structure in H I emission as well as self-absorption and in the structure of polarized emission. These observations demonstrated that no phase of the ISM, including the transition from atomic gas to star formation, can be understood in isolation. Canadians have also played leading roles in the characterization of dust with Planck and balloon-borne telescopes. Canadians have also used pulsar scintillometry to measure structure in the ISM at the smallest scales, below 1 AU. The 2020s offer many opportunities for ISM science in Canada. A major but cost-effective upgrade to the Synthesis Telescope with broadband (400-1800 MHz) single-pixel feeds would enable broadband polarimetry as well as wide-area, arcminute surveys of radio recombination lines. The next generation of balloon-borne telescopes will investigate magnetic fields and dust properties. Large single dishes, particularly the Green Bank Telescope, remain essential for our understanding of the diffuse structure which characterizes the ISM. Very long baseline interferometry capability enables parallax measurements of pulsars and masers and for further progress in scintillometry. ISM astronomers will continue to participate in cosmological experiments including CHIME and CHORD. Protecting quiet radio frequency interference environments will be ever more important as broadband observations are ever more central to ISM science. Computational capability is essential both for numerical work and for handling the observational data.

We present a detailed analysis of deep VLT/X-Shooter observations of the planetary nebula Tc 1. We calculate gas temperature, density, extinction, and abundances for several species from the empirical analysis of the total line fluxes. In addition, a spatially resolved analysis of the most intense lines provides the distribution of such quantities across the nebula. The new data reveal that several lines exhibit a double peak spectral profile consistent with the blue- and red-shifted components of an expanding spherical shell. The study of such components allowed us to construct for the first time a three-dimensional morphological model, which reveals that Tc 1 is a slightly elongated spheroid with an equatorial density enhancement seen almost pole on. A few bright lines present extended wings (with velocities up to a few hundred km/s), but the mechanism producing them is not clear. We constructed photoionization models for the main shell of Tc 1. The models predict the central star temperature and luminosity, as well as the nebular density and abundances similar to previous studies. Our models indicate that Tc 1 is located at a distance of approximately 2 kpc. We report the first detection of the [Kr III] 6825 A emission line, from which we determine the Krypton abundance. Our model indicates that the main shell of Tc 1 is matter bounded; leaking H ionizing photons may explain the ionization of its faint AGB-remnant halo.

We present the highest resolution single-dish submillimetre observations of the detached shell source U Antliae to date. The observations were obtained at $450~\micron$ and $850~\micron$ with SCUBA-2 instrument on the James Clerk Maxwell Telescope as part of the Nearby Evolved Stars Survey. The emission at $850~\micron$ peaks at $40\arcsec$ with hints of a second peak seen at $\sim 20\arcsec$. The emission can be traced out to a radius of $56\arcsec$ at a $3\sigma$ level. The outer peak observed at $850~\micron$ aligns well with the peak observed at Herschel/PACS wavelengths. With the help of spectral energy distribution fitting and radiative transfer calculations of multiple-shell models for the circumstellar envelope, we explore the various shell structures and the variation of grain sizes along the in the circumstellar envelope. We determine a total shell dust mass of $(2.0 \pm 0.3) \times 10^{-5}$ M$_{\odot}$ and established that the thermal pulse which gave rise to the detached shell occurred 3500 $\pm$ 500 years ago.

Recent advances in laboratory spectroscopy lead to the claim of ionized Buckminsterfullerene (C60+) as the carrier of two diffuse interstellar bands (DIBs) in the near-infrared. However, irrefutable identification of interstellar C60+ requires a match between the wavelengths and the expected strengths of all absorption features detectable in the laboratory and in space. Here we present Hubble Space Telescope (HST) spectra of the region covering the C60+ 9348, 9365, 9428 and 9577 {\AA} absorption bands toward seven heavily-reddened stars. We focus in particular on searching for the weaker laboratory C60+ bands, the very presence of which has been a matter for recent debate. Using the novel STIS-scanning technique to obtain ultra-high signal-to-noise spectra without contamination from telluric absorption that afflicted previous ground-based observations, we obtained reliable detections of the (weak) 9365, 9428 {\AA} and (strong) 9577 {\AA} C60+ bands. The band wavelengths and strength ratios are sufficiently similar to those determined in the latest laboratory experiments that we consider this the first robust identification of the 9428 {\AA} band, and a conclusive confirmation of interstellar C60+.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

We report a computational study of the stability and infrared (IR) vibrational spectra of neutral and singly ionised fullerene cages containing between 44 and 70 carbon atoms. The stability is characterised in terms of the standard enthalpy of formation per CC bond, the HOMO-LUMO gap, and the energy required to eliminate a C$_2$ fragment. We compare the simulated IR spectra of these fullerene species to the observed emission spectra of several planetary nebulae (Tc 1, SMP SMC 16, and SMP LMC 56) where strong C$_{60}$ emission has been detected. Although we could not conclusively identify fullerenes other than C$_{60}$ and C$_{70}$, our results point to the possible presence of smaller (44, 50, and 56-atom) cages in those astronomical objects. Observational confirmation of our prediction should become possible when the James Webb Space Telescope comes online.

We report new spectroscopic observations of the magnetic Of?p star HD 148937 obtained since 2015 that differ qualitatively from its extensive historical record of weak, periodic spectral variations. This remarkable behaviour represents clear evidence for an unprecedented change in the character of variability of the star. In this paper we describe the new spectral properties and compare them to the previous line profiles. Based on measurements of the radial velocities of the C iii/N iii} emission lines near 4640 \AA\ and the C iv absorption lines near 5800 \AA, we infer that HD 148937 is likely a high-mass, double-lined spectroscopic binary. Combining the spectroscopic orbit with an archival interferometric measurement of the apparent separation of the equal-brightness components, we tentatively conclude that HD 148937 consists of two O-type stars with masses of approximately 34 and $49~M_\odot$, orbiting in an eccentric ($e=0.75$), long-period ($P_{\rm orb}\sim 26$ y) orbit. We discuss the potential relationship of the binary system to the peculiar properties of HD 148937, and propose future observations to refine the orbital and stellar properties.

We report cosmic ray ionization rates towards ten reddened stars studied within the framework of the EDIBLES (ESO Diffuse Interstellar Bands Large Exploration Survey) program, using the VLT-UVES. For each sightline, between 2 and 10 individual rotational lines of OH$^+$ have been detected in its (0,0) and (1,0) $A^3\Pi-X^3\Sigma^-$ electronic band system. This allows constraining of OH$^+$ column densities towards different objects. Results are also presented for 28 additional sightlines for which only one or rather weak signals are found. An analysis of these data makes it possible to derive the primary cosmic ray ionization rate $\zeta_p$ in the targeted diffuse interstellar clouds. For the ten selected targets, we obtain a range of values for $\zeta_p$ equal to $(3.9-16.4) \times 10^{-16}~\mathrm{s}^{-1}$. These values are higher than the numbers derived in previous detections of interstellar OH$^+$ in the far-infrared / sub-millimeter-wave regions and in other near-ultraviolet studies. This difference is a result of using new OH$^+$ oscillator strength values and a more complete picture of all relevant OH$^+$ formation and destruction routes (including the effect of proton recombinations on PAHs), and the relatively high $N$(OH$^+$) seen toward those ten targets.

About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides in this intermediate region. We put forward that the next frontier for AMO science is to explore the AMO complexity that describes most of the Cosmos.

We present JCMT SCUBA-2 $450\mu$m and $850\mu$m observations of 14 Asymptotic Giant Branch (AGB) stars (9 O--rich, 4 C-rich and 1 S--type) and one Red Supergiant (RSG) in the Solar Neighbourhood. We combine these observations with \emph{Herschel}/PACS observations at $70\mu$m and $160\mu$m and obtain azimuthally-averaged surface-brightness profiles and their PSF subtracted residuals. The extent of the SCUBA-2 850 $\mu$m emission ranges from 0.01 to 0.16 pc with an average of $\sim40\%$ of the total flux being emitted from the extended component. By fitting a modified black-body to the four-point SED at each point along the radial profile we derive the temperature ($T$), spectral index of dust emissivity ($\beta$) and dust column density ($\Sigma$) as a function of radius. For all the sources, the density profile deviates significantly from what is expected for a constant mass-loss rate, showing that all the sources have undergone variations in mass-loss during this evolutionary phase. In combination with results from CO line emission, we determined the dust-to-gas mass ratio for all the sources in our sample. We find that, when sources are grouped according to their chemistry, the resulting average dust-to-gas ratios are consistent with the respective canonical values. However we see a range of values with significant scatter which indicate the importance of including spatial information when deriving these numbers.

The so-called C2-DIBs are a class of very weak bands that fall in the blue part of the optical spectrum and are associated with high column densities of the C2 molecule. DIB profile structures constrain potential molecular carriers, but their measurement requires high S/N and spectra and the use of sightlines without Doppler splitting, as typical for a single-cloud situation. Spectra from the ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES) conducted at the VLT (ESO/Paranal) were explored to identify single-cloud and high C2 column sightlines, extract the corresponding C2-DIBs and study their strengths and profiles, and to investigate in detail any sub-structures. The target selection was made based on profile-fitting of the Sodium doublets and the detection of C2 lines. The C2 (2-0) Phillips system was fitted using a physical model of the host cloud. C2 column densities, temperatures as well as gas densities were derived for each sightline. 18 known C2-DIBs and 8 strong non-C2 DIBs were extracted towards 8 targets, comprising 7 single-cloud and one multi-cloud line-of-sights. Correlational studies revealed a tight association of the former group with the C2 columns, whereas the non-C2 DIBs are primarily correlated with reddening. We report three new weak diffuse band candidates. We show for the first time that at least 14 C2-DIBs exhibit spectral sub-structures which are consistent with unresolved rotational branches of molecular carriers. The variability of their peak separations among the bands for a given sightline implies that their carriers are different molecules with quite different sizes. We also illustrate how profiles of the same DIB vary among targets and as a function of physical parameters and provide tables defining the sub-structures to be compared with future models and experimental results.

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

Gas phase spectroscopic laboratory experiments for the buckminsterfullerene cation C60+ resulted in accurate rest wavelengths for five C60+ transitions that have been compared with diffuse interstellar bands (DIBs) in the near infra-red. Detecting these in astronomical spectra is difficult due to the strong contamination of ground-based spectra by atmospheric water vapor, to the presence of weak and shallow stellar lines and/or blending with other weak DIBs. The detection of the two strong bands has been claimed by several teams, and the three additional and weaker bands have been detected in a few sources. Certain recent papers have argued against the identification of C60+ based on spectral analyses claiming (i) a large variation in the ratio between the equivalent widths of the 9632 and 9577\AA\: bands, (ii) a large redshift of the 9632\AA\: band for the Orion star HD 37022, and (iii) the non-detection of the weaker 9428\AA~DIB. Here we address these three points. (i) We show that the model stellar line correction for the 9632\AA~DIB overestimates the difference between the strengths of the lines in giant and dwarf star spectra, casting doubts on the conclusions about the ratio variability. (ii) Using high quality stellar spectra from the ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES), recorded with the ESO/Paranal Ultraviolet Echelle Spectrograph (UVES) in about the same atmospheric conditions, we find no wavelength shift in the 9632\AA\ band towards HD 37022. (iii) Using EDIBLES spectra and data from the Echelle SpectroPolarimetric Device for the Observation of Stars (ESPaDOnS) at CFHT we show that the presence of a weak 9428\AA\ band cannot be ruled out, even in the same observations that a previous study claimed it was not present.

It is well known that the dust properties of the diffuse interstellar medium exhibit variations towards different sight-lines on a large scale. We have investigated the variability of the dust characteristics on a small scale, and from cloud-to-cloud. We use low-resolution spectro-polarimetric data obtained in the context of the Large Interstellar Polarisation Survey (LIPS) towards 59 sight-lines in the Southern Hemisphere, and we fit these data using a dust model composed of silicate and carbon particles with sizes from the molecular to the sub-micrometre domain. Large (> 6 nm) silicates of prolate shape account for the observed polarisation. For 32 sight-lines we complement our data set with UVES archive high-resolution spectra, which enable us to establish the presence of single-cloud or multiple-clouds towards individual sight-lines. We find that the majority of these 35 sight-lines intersect two or more clouds, while eight of them are dominated by a single absorbing cloud. We confirm several correlations between extinction and parameters of the Serkowski law with dust parameters, but we also find previously undetected correlations between these parameters that are valid only in single-cloud sight-lines. We find that interstellar polarisation from multiple-clouds is smaller than from single-cloud sight-lines, showing that the presence of a second or more clouds depolarises the incoming radiation. We find large variations of the dust characteristics from cloud-to-cloud. However, when we average a sufficiently large number of clouds in single-cloud or multiple-cloud sight-lines, we always retrieve similar mean dust parameters. The typical dust abundances of the single-cloud cases are [C]/[H] = 92 ppm and [Si]/[H] = 20 ppm.

Polarimetric studies of light transmitted through interstellar clouds may give constraints on the properties of the interstellar dust grains. Traditionally, broadband linear polarisation (BBLP) measurements have been considered an important diagnostic tool for the study of the interstellar dust, while comparatively less attention has been paid to spectropolarimetric measurements. However, spectropolarimetry may offer stronger constraints than BBLP, for example by revealing narrowband features, and by allowing us to distinguish the contribution of dust from the contribution of interstellar gas. Therefore, we have decided to carry out a Large Interstellar Polarisation Survey (LIPS) using spectropolarimetric facilities in both hemispheres. Here we present the results obtained in the Southern Hemisphere with the FORS2 instrument of the ESO Very Large Telescope. Our spectra cover the wavelength range 380--950\,nm at a spectral resolving power of about 880. We have produced a publicly available catalogue of 127 linear polarisation spectra of 101 targets. We also provide the Serkowski-curve parameters, as well as the wavelength gradient of the polarisation position angle for the interstellar polarisation along 76 different lines of sight. In agreement with previous literature, we found that the best-fit parameters of the Serkowski-curve are not independent of each other. However, the relationships that we obtained are not always consistent with what was found in previous studies.

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

Due to recent advances in laboratory spectroscopy, the first optical detection of a very large molecule has been claimed in the diffuse interstellar medium (ISM): C60+ (ionized Buckminsterfullerene). Confirming the presence of this molecule would have significant implications regarding the carbon budget and chemical complexity of the ISM. Here we present results from a new method for ultra-high signal-to-noise (S/N) spectroscopy of background stars in the near infrared (at wavelengths 0.9-1 micron), using the Hubble Space Telescope Imaging Spectrograph (STIS) in a previously untested `STIS scan' mode. The use of HST provides the crucial benefit of eliminating the need for error-prone telluric correction methods in the part of the spectrum where the C60+ bands lie, and terrestrial water vapor contamination is severe. Our STIS spectrum of the heavily-reddened B0 star BD63\,1964 reaches an unprecedented S/N for this instrument ($\sim600-800$), allowing the detection of the diffuse interstellar band (DIB) at 9577 \AA\ attributed to C60+ as well as new DIBs in the near-IR. Unfortunately, the presence of overlapping stellar lines, and the unexpected weakness of the C60+ bands in this sightline, prevents conclusive detection of the weaker C60+ bands. A probable correlation between the 9577 \AA\ DIB strength and interstellar radiation field is identified, which suggests that more strongly-irradiated interstellar sightlines will provide the optimal targets for future C60+ searches.

We present a principal component analysis of 23 line of sight parameters (including the strengths of 16 diffuse interstellar bands, DIBs) for a well-chosen sample of single-cloud sightlines representing a broad range of environmental conditions. Our analysis indicates that the majority ($\sim$93\%) of the variations in the measurements can be captured by only four parameters The main driver (i.e., the first principal component) is the amount of DIB-producing material in the line of sight, a quantity that is extremely well traced by the equivalent width of the $\lambda$5797 DIB. The second principal component is the amount of UV radiation, which correlates well with the $\lambda$5797/$\lambda$5780 DIB strength ratio. The remaining two principal components are more difficult to interpret, but are likely related to the properties of dust in the line of sight (e.g., the gas-to-dust ratio). With our PCA results, the DIBs can then be used to estimate these line of sight parameters.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We performed a detailed spectroscopic analysis of the fullerene C60-containing planetary nebula (PN) Lin49 in the Small Magellanic Cloud using XSHOOTER at the ESO VLT and the Spitzer/IRS instruments. We derived nebular abundances for nine elements. We used TLUSTY to derive photospheric parameters for the central star. Lin49 is C-rich and metal-deficient PN (Z~0.0006). The nebular abundances are in good agreement with Asymptotic Giant Branch nucleosynthesis models for stars with initial mass 1.25 Msun and metallicity Z = 0.001. Using the TLUSTY synthetic spectrum of the central star to define the heating and ionising source, we constructed the photoionisation model with CLOUDY that matches the observed spectral energy distribution (SED) and the line fluxes in the UV to far-IR wavelength ranges simultaneously. We could not fit the ~1-5 um SED using a model with 0.005-0.1 um-sized graphite grains and a constant hydrogen density shell owing to the prominent near-IR excess, while at other wavelengths the model fits the observed values reasonably well. We argue that the near-IR excess might indicate either (1) the presence of very small particles in the form of small carbon clusters, small graphite sheets, or fullerene precursors, or (2) the presence of a high-density structure surrounding the central star. We found that SMC C60 PNe show a near-IR excess component to lesser or greater degree. This suggests that these C60 PNe might maintain a structure nearby their central star.

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.

We present the results of a computational study to investigate the infrared spectroscopic properties of a large number of polycyclic aromatic hydrocarbon (PAH) molecules and their fully dehydrogenated counterparts. We constructed a database of fully optimized geometries for PAHs that is complete for eight or fewer fused benzene rings, thus containing 1550 PAHs and 805 fully dehydrogenated aromatics. A large fraction of the species in our database have clearly non-planar or curved geometries. For each species, we determined the frequencies and intensities of their normal modes using density functional theory calculations. Whereas most PAH spectra are fairly similar, the spectra of fully dehydrogenated aromatics are much more diverse. Nevertheless, these fully dehydrogenated species show characteristic emission features at 5.2$\mu$m, 5.5$\mu$m and 10.6$\mu$m; at longer wavelengths, there is a forest of emission features in the 16--30$\mu$m range that appears as a structured continuum, but with a clear peak centered around 19$\mu$m. We searched for these features in Spitzer-IRS spectra of various positions in the reflection nebula NGC 7023. We find a weak emission feature at 10.68$\mu$m in all positions except that closest to the central star. We also find evidence for a weak 19$\mu$m feature at all positions that is not likely due to C$_{60}$. We interpret these features as tentative evidence for the presence of a small population of fully dehydrogenated PAHs, and discuss our results in the framework of PAH photolysis and the formation of fullerenes.

We present a comprehensive and sensitive unbiased survey of interstellar features in the near-UV range (3050-3700 {\AA}). We combined a large number of VLT/UVES archival observations of a sample of highly reddened early type stars -- typical diffuse interstellar band (DIB) targets -- and unreddened standards. We stacked the individual observations to obtain a reddened "superspectrum" in the interstellar rest frame with a signal-to-noise (S/N) ratio exceeding 1500. We compared this to the analogous geocentric and stellar rest frame superspectra as well as to an unreddened superspectrum to find interstellar absorption features. We find 30 known features (11 atomic and 19 molecular) and tentatively detect up to 7 new interstellar absorption lines of unknown origin. Our survey is sensitive to narrow and weak features; telluric residuals preclude us from detecting broader features. For each sightline, we measured fundamental parameters (radial velocities, line widths, and equivalent widths) of the detected interstellar features. We also revisit upper limits for the column densities of small, neutral polycyclic aromatic hydrocarbon (PAH) molecules that have strong transitions in this wavelength range.

HR 4049 is a peculiar evolved binary which is surrounded by a circumbinary disk. Mid-infrared observations show that the disk is rich in molecular gas and radially extended. To study the properties of this disk, we re-analyzed a set of near-infrared observations at high spectral resolution obtained with Gemini-Phoenix. These data cover absorption lines originating from the first overtone of CO and from \hoh in the 2.3 $\mu$m region as well as more complex emission-absorption profiles from \hoh and the fundamental mode of CO near 4.6 $\mu$m. By using an excitation diagram and from modeling the spectrum, we find that most of the CO overtone and \hoh absorption originates from hot gas ($T_{\rm ex} \approx 1000$~K) with high column densities, consistent with the mid-infrared data. The strong emission in the wavelength range of the CO fundamental furthermore suggests that there is a significant quantity of gas in the inner cavity of the disk. In addition, there is a much colder component in the line of sight to the disk. A detailed analysis of the overtone line profiles reveals variations in the line widths which are consistent with a radially extended disk in Keplerian rotation with hotter gas closer to the central star. We estimate the mass of the primary to be $\sim 0.34$~M$_{\odot}$ and discuss the implications for its evolutionary status.

This paper presents a spectral survey of diffuse interstellar bands (DIBs) in the NIR range, from 0.9 to 2.5 micron. The observations were designed to detect new DIBs, confirm previously proposed NIR DIBs, and characterise their behaviour with respect to known line-of-sight properties (including the optical DIBs present in our spectra). X-shooter at the VLT was used to obtained medium-resolution spectra of eight known DIB targets and one telluric reference star. In addition to the known 9577, 9632, 10780, 11797, and 13175 Angstroms NIR DIBs, we confirm 9 out of the 13 NIR DIBs that were presented by Geballe and co-workers in 2011. Furthermore, we report 12 new NIR DIB candidates. The strengths of the strongest NIR DIBs show a general correlation with reddening, E(B-V), but with a large scatter. Several NIR DIBs are more strongly correlated with the 5780 Angstroms DIB strength than with E(B-V); this is especially the case for the 15268 Angstroms DIB. The NIR DIBs are strong: the summed equivalent widths of the five strongest NIR DIBs represent a few percent of the total equivalent width of the entire average DIB spectrum (per unit reddening). The NIR DIBs towards the translucent cloud HD 147889 are all weak with respect to the general trend. No direct match was found between observed NIR DIBs and laboratory matrix-isolation spectroscopic data of PAHs. The strong correlation between the 5780-15268 DIB pair implies that (Nf)_5780 / (Nf)_15268 = 14. However, the reduced strength of the 15268 Angstroms DIB in HD 147889 rules out a common carrier for these two DIBs. Since the ionisation fraction for small PAHs in this translucent cloud is known to be low compared to diffuse clouds, the weakness of the 15268 Angstrom DIB suggests that an ionised species could be the carrier of this NIR DIB. (abridged).

Fullerenes are a particularly stable class of carbon molecules in the shape of a hollow sphere or ellipsoid that might be formed in the outflows of carbon stars. Once injected into the interstellar medium (ISM), these stable species survive and are thus likely to be widespread in the Galaxy where they contribute to interstellar extinction, heating processes, and complex chemical reactions. In recent years, the fullerene species C60 (and to a lesser extent C70) have been detected in a wide variety of circumstellar and interstellar environments showing that when conditions are favourable, fullerenes are formed efficiently. Fullerenes are the first and only large aromatics firmly identified in space. The detection of fullerenes is thus crucial to provide clues as to the key chemical pathways leading to the formation of large complex organic molecules in space, and offers a great diagnostic tool to describe the environment in which they reside. Since fullerenes share many physical properties with PAHs, understanding how fullerenes form, evolve and respond to their physical environment will yield important insights into one of the largest reservoirs of organic material in space. In spite of all these detections, many questions remain about precisely which members of the fullerene family are present in space, how they form and evolve, and what their excitation mechanism is. We present here an overview of what we know from astronomical observations of fullerenes in these different environments, and discuss current thinking about the excitation process. We highlight the various formation mechanisms that have been proposed, discuss the physical conditions conducive to the formation and/or detection of fullerenes in carbon stars, and their possible connection to PAHs, HACs and other dust features.

Infrared spectra of carbon-rich objects which have evolved off the asymptotic giant branch reveal a range of dust properties, including fullerenes, polycyclic aromatic hydrocarbons (PAHs), aliphatic hydrocarbons, and several unidentified features, including the 21 um emission feature. To test for the presence of fullerenes, we used the position and width of the feature at 18.7-18.9 um and examined other features at 17.4 and 6-9 um. This method adds three new fullerene sources to the known sample, but it also calls into question three previous identifications. We confirm that the strong 11 um features seen in some sources arise primarily from SiC, which may exist as a coating around carbonaceous cores and result from photo-processing. Spectra showing the 21 um feature usually show the newly defined Class D PAH profile at 7-9 um. These spectra exhibit unusual PAH profiles at 11-14 um, with weak contributions at 12.7 um, which we define as Class D1, or show features shifted to ~11.4, 12.4, and 13.2 um, which we define as Class D2. Alkyne hydrocarbons match the 15.8 um feature associated with 21 um emission. Sources showing fullerene emission but no PAHs have blue colors in the optical, suggesting a clear line of sight to the central source. Spectra with 21 um features and Class D2 PAH emission also show photometric evidence for a relatively clear line of sight to the central source. The multiple associations of the 21 um feature to aliphatic hydrocarbons suggest that the carrier is related to this material in some way.

We searched the Spitzer Space Telescope data archive for Galactic planetary nebulae (PNe), that show the characteristic 17.4 and um features due to C60, also known as buckminsterfullerene. Out of 338 objects with Spitzer/IRS data, we found eleven C60-containing PNe, six of which (Hen2-68, IC2501, K3-62, M1-6, M1-9, and SaSt2-3) are new detections. We analyzed the spectra, along with ancillary data, using the photo-ionization code CLOUDY to establish the atomic line fluxes, and determine the properties of the radiation field, as set by the effective temperature of the central star. In addition, we measured the infrared spectral features due to dust grains. We find that the Polycyclic Aromatic Hydrocarbon (PAH) profile over 6-9 um in these C60-bearing carbon-rich PNe is of the more chemically-processed class A. The intensity ratio of 3.3 um to 11.3 um PAH indicates that the number of C-atoms per PAH in C60-containing PNe is small compared to that in non-C60 PNe. The Spitzer spectra also show broad dust features around 11 and 30 um. Analysis of the 30-um feature shows that it is strongly correlated with the continuum, and we propose that a single carbon-based carrier is responsible for both the continuum and the feature. The strength of the 11-um feature is correlated to the temperature of the dust, suggesting that it is at least partially due to a solid-state carrier. The chemical abundances of C60-containing PNe can be explained by AGB nucleosynthesis models for initially 1.5-2.5 Msun stars with Z=0.004. We plotted the locations of C60-containing PNe on a face-on map of the Milky Way and we found that most of these PNe are outside the solar circle, consistent with low metallicity values. Their metallicity suggests that the progenitors are an older population.

We present a detailed analysis of the mid-infrared spectrum of the peculiar evolved object HR 4049. The full Spitzer-IRS high-resolution spectrum shows a wealth of emission with prominent features from CO$_2$ and H$_2$O and possible contributions from HCN and OH. We model the molecular emission and find that it originates from a massive ($M \gtrsim 8 \times 10^{-3}$ M$_{\odot}$), warm ($T_{\rm ex} \approx 500$ K) and radially extended gas disk that is optically thick at infrared wavelengths. We also report less enrichment in $^{17}$O and $^{18}$O than previously found and a comparison of the Spitzer observations to earlier data obtained by ISO-SWS reveals that the CO$_2$ flux has more than doubled in 10 years time, indicating active and ongoing chemical evolution in the circumbinary disk. If the gas originates from interaction between the stellar wind and the dust, this suggests that the dust could be oxygen-rich in nature. The molecular gas plays a crucial role in the thermal properties of the circumbinary disk by allowing visible light to heat the dust and then trapping the infrared photons emitted by the dust. This results in higher temperatures and a more homogeneous temperature structure in the disk.

Fullerenes have recently been identified in space and they may play a significant role in the gas and dust budget of various astrophysical objects including planetary nebulae (PNe), reflection nebulae (RNe) and H II regions. The tenuous nature of the gas in these environments precludes the formation of fullerene materials following known vaporization or combustion synthesis routes even on astronomical timescales. We have studied the processing of hydrogenated amorphous carbon (a-C:H or HAC) nano-particles and their specific derivative structures, which we name "arophatics", in the circumstellar environments of young, carbon-rich PNe. We find that UV-irradiation of such particles can result in the formation of fullerenes, consistent with the known physical conditions in PNe and with available timescales.

In recent years, the fullerene species C60 (and to a lesser extent also C70) has been reported in the mid-IR spectra of various astronomical objects. Cosmic fullerenes form in the circumstellar material of evolved stars, and survive in the interstellar medium (ISM). It is not entirely clear how they form or what their excitation mechanism is.

We compare and analyze the Spitzer mid-infrared spectrum of three fullerene-rich planetary nebulae in the Milky Way and the Magellanic Clouds; Tc1, SMP SMC16, and SMP LMC56. The three planetary nebulae share many spectroscopic similarities. The strongest circumstellar emission bands correspond to the infrared active vibrational modes of the fullerene species C60 and little or no emission is present from Polycyclic Aromatic Hydrocarbons (PAHs). The strength of the fullerene bands in the three planetary nebulae is very similar, while the ratio of the [NeIII]15.5um/[NeII]12.8um fine structure lines, an indicator of the strength of the radiation field, is markedly different. This raises questions about their excitation mechanism and we compare the fullerene emission to fluorescent and thermal models. In addition, the spectra show other interesting and common features, most notably in the 6-9um region, where a broad plateau with substructure dominates the emission. These features have previously been associated with mixtures of aromatic/aliphatic hydrocarbon solids. We hypothesize on the origin of this band, which is likely related to the fullerene formation mechanism, and compare it with modeled Hydrogenated Amorphous Carbon that present emission in this region.

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

We report the detection of carbon monoxide (CO) emission from the young supernova remnant Cassiopeia A (Cas A) at wavelengths corresponding to the fundamental vibrational mode at 4.65 micron. We obtained AKARI Infrared Camera spectra towards 4 positions which unambiguously reveal the broad characteristic CO ro-vibrational band profile. The observed positions include unshocked ejecta at the center, indicating that CO molecules form in the ejecta at an early phase. We extracted a dozen spectra across Cas A along the long 1 arcmin slits, and compared these to simple CO emission models in Local Thermodynamic Equilibrium to obtain first-order estimates of the excitation temperatures and CO masses involved. Our observations suggest that significant amounts of carbon may have been locked up in CO since the explosion 330 years ago. Surprisingly, CO has not been efficiently destroyed by reactions with ionized He or the energetic electrons created by the decay of the radiative nuclei. Our CO detection thus implies that less carbon is available to form carbonaceous dust in supernovae than is currently thought and that molecular gas could lock up a significant amount of heavy elements in supernova ejecta.

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.

Carbon-rich evolved stars from the asymptotic giant branch to the planetary nebula phase are characterized by a rich and complex carbon chemistry in their circumstellar envelopes. A peculiar object is the preplanetary nebula SMP LMC 11, whose Spitzer-IRS spectrum shows remarkable and diverse molecular absorption bands. To study how the molecular composition in this object compares to our current understanding of circumstellar carbon chemistry, we modeled this molecular absorption. We find high abundances for a number of molecules, perhaps most notably benzene. We also confirm the presence of propyne (CH3C2H) in this spectrum. Of all the cyanopolyynes, only HC3N is evident; we can detect at best a marginal presence of HCN. From comparisons to various chemical models, we can conclude that SMP LMC 11 must have an unusual circumstellar environment (a torus rather than an outflow).

We recently identified several emission bands in the Spitzer-IRS spectrum of the unusual planetary nebula Tc 1 with the infrared active vibrational modes of the neutral fullerene species C60 and C70. Since then, the fullerene bands have been detected in a variety of sources representing circumstellar and interstellar environments. Abundance estimates suggest that C60 represents ~0.1%-1.5% of the available carbon in those sources. The observed relative band intensities in various sources are not fully compatible with single-photon heating and fluorescent cooling, and are better reproduced by a thermal distribution at least in some sources. The observational data suggests that fullerenes form in the circumstellar environments of evolved stars, and survive in the interstellar medium. Precisely how they form is still a matter of debate.

Context. The circumstellar environment of evolved stars is generally rich in molecular gas and dust. Typically, the entire environment is either oxygen-rich or carbon-rich, depending on the evolution of the central star. Aims. In this paper we discuss three evolved disc sources with evidence of atypical emission lines in their infrared spectra. The stars were taken from a larger sample of post-AGB binaries for which we have Spitzer infrared spectra, characterised by the presence of a stable oxygen-rich circumbinary disc. Our previous studies have shown that the infrared spectra of post-AGB disc sources are dominated by silicate dust emission, often with an extremely high crystallinity fraction. However, the three sources described here are selected because they show a peculiar molecular chemistry. Methods. Using Spitzer infrared spectroscopy, we study in detail the peculiar mineralogy of the three sample stars. Using the observed emission features, we identify the different observed dust, molecular and gas species. Results. The infrared spectra show emission features due to various oxygen-rich dust components, as well as CO2 gas. All three sources show the strong infrared bands generally ascribed to polycyclic aromatic hydrocarbons. Furthermore, two sample sources show C60 fullerene bands. Conclusions. Even though the majority of post-AGB disc sources are dominated by silicate dust in their circumstellar environment, we do find evidence that, for some sources at least, additional processing must occur to explain the presence of large carbonaceous molecules. There is evidence that some of these sources are still oxygen-rich, which makes the detection of these molecules even more surprising.

Context: Mass-loss occurring in red supergiants (RSGs) is a major contributor to the enrichment of the interstellar medium in dust and molecules. The physical mechanism of this mass loss is however relatively poorly known. Betelgeuse is the nearest RSG, and as such a prime object for high angular resolution observations of its surface (by interferometry) and close circumstellar environment. Aims: The goal of our program is to understand how the material expelled from Betelgeuse is transported from its surface to the interstellar medium, and how it evolves chemically in this process. Methods: We obtained diffraction-limited images of Betelgeuse and a PSF calibrator (Aldebaran) in six filters in the N band (7.76 to 12.81 mic) and two filters in the Q band (17.65 and 19.50 mic), using the VLT/VISIR instrument. Results: Our images show a bright, extended and complex circumstellar envelope at all wavelengths. It is particularly prominent longwards of 9-10 mic, pointing at the presence of O-rich dust, such as silicates or alumina. A partial circular shell is observed between 0.5 and 1.0" from the star, and could correspond to the inner radius of the dust envelope. Several knots and filamentary structures are identified in the nebula, one of the knots, located at a distance of 0.9" west of the star, being particularly bright and compact. Conclusions: The circumstellar envelope around Betelgeuse extends at least up to several tens of stellar radii. Its relatively high degree of clumpiness indicates an inhomogeneous spatial distribution of the material lost by the star. Its extension corresponds to an important intermediate scale, where most of the dust is probably formed, between the hot and compact gaseous envelope observed previously in the near infrared and the interstellar medium.

The origin of the narrow optical emission bands seen toward the Red Rectangle is not yet understood. In this paper we investigate further the proposal that these are due to luminescence of large carbonaceous molecules. Polarised signals of several percent could be expected from certain asymmetric molecular rotators. The ESPaDOnS echelle spectrograph mounted at the CFHT was used to obtain high-resolution optical spectropolarimetric data of the Red Rectangle nebular emission. The RRBs at 5800, 5850, and 6615 Angstrom are detected in spectra of the nebular emission 7" and 13" North-East from the central star. The 5826 and 6635 Angstrom RRB are detected only at the position nearest to the central star. For both positions the Stokes Q and U spectra show no unambiguous polarisation signal in any of the RRBs. We derive an upper limit of 0.02% line polarisation for these RRBs. A tentative feature with peak polarisation of 0.05% is seen for the 5800 RRB at 7" offset. However, the Null spectra suggest that this may be an instrumental artifact. The lack of a clear polarisation signal for the five detected RRBs implies that, if the emission is due to luminescence of complex organics, these gas-phase molecular carriers are likely to have a high degree of symmetry, as they do not exhibit a Q-branch in their rotational profile, although this may be modified by statistical effects.

Studying the composition of dust in the interstellar medium (ISM) is crucial in understanding the cycle of dust in our galaxy. The mid-infrared spectral signature of amorphous silicates, the most abundant dust species in the ISM, is studied in different lines-of-sight through the Galactic plane, thus probing different conditions in the ISM. We have analysed 10 spectra from the Spitzer archive, of which 6 lines-of-sight probe diffuse interstellar medium material and 4 probe molecular cloud material. The 9.7 um silicate absorption features in 7 of these spectra were studied in terms of their shape and strength. In addition, the shape of the 18 um silicate absorption features in 4 of the diffuse sightline spectra were analysed. The 9.7 um silicate absorption bands in the diffuse sightlines show a strikingly similar band shape. This is also the case for all but one of the 18 um silicate absorption bands observed in diffuse lines-of-sight. The 9.7 um bands in the 4 molecular sightlines show small variations in shape. These modest variations in the band shape are inconsistent with the interpretation of the large variations in {\tau}_9.7/E(J-K) between diffuse and molecular sightlines in terms of silicate grain growth. Instead, we suggest that the large changes in {\tau}_9.7 / E(J-K) must be due to changes in E(J-K).

Submillimetre continuum radiation allows us to probe cold objects, particularly the earliest, dusty phases of star formation, high-redshift galaxies and circumstellar disks. The submillimetre window gives a unique view of the physical and dynamical conditions in the neutral and molecular interstellar medium. In the next decade a combination of wide-field surveys with single-dish telescopes and targeted follow-up with ALMA and other facilities should enable rapid progress in answering questions about the origins of planetary systems, stars and galaxies.

Context: Betelgeuse is one the largest stars in the sky in terms of angular diameter. Structures on the stellar photosphere have been detected in the visible and near-infrared as well as a compact molecular environment called the MOLsphere. Mid-infrared observations have revealed the nature of some of the molecules in the MOLsphere, some being the precursor of dust. Aims: Betelgeuse is an excellent candidate to understand the process of mass loss in red supergiants. Using diffraction-limited adaptive optics (AO) in the near-infrared, we probe the photosphere and close environment of Betelgeuse to study the wavelength dependence of its extension, and to search for asymmetries. Methods: We obtained AO images with the VLT/NACO instrument, taking advantage of the "cube" mode of the CONICA camera to record separately a large number of short-exposure frames. This allowed us to adopt a "lucky imaging" approach for the data reduction, and obtain diffraction-limited images over the spectral range 1.04-2.17 $\mu$m in 10 narrow-band filters. Results: In all filters, the photosphere of Betelgeuse appears partly resolved. We identify an asymmetric envelope around the star, with in particular a relatively bright "plume" extending in the southwestern quadrant up to a radius of approximately six times the photosphere. The CN molecule provides an excellent match to the 1.09 mic bandhead in absorption in front of the stellar photosphere, but the emission spectrum of the plume is more difficult to interpret. Conclusions: Our AO images show that the envelope surrounding Betelgeuse has a complex and irregular structure. We propose that the southwestern plume is linked either to the presence of a convective hot spot on the photosphere, or to the rotation of the star.

Context: Red super-giant (RSG) stars exhibit significant mass loss through a slow and dense wind. They are often considered to be the more massive counter parts of Asymptotic Giant Branch (AGB) stars. While the AGB mass loss is linked to their strong pulsations, the RSG are often only weakly variable. Aim: To study the conditions at the base of the wind, by determining the dust composition in a sample of RSG. The dust composition is thought to be sensitive to the density, temperature and acceleration at the base of the wind. Method: We compile a sample of 27 RSG infrared spectra (ISO-SWS) and supplement these with photometric measurements to obtain the full spectral energy distribution (SED). These data are modelled using a dust radiative transfer code. The results are scrutinised for correlations. Results: We find (1) strong correlations between dust composition, mass-loss rate and stellar luminosity, roughly in agreement with the theoretical dust condensation sequence, (2) the need for a continuous (near-)IR dust opacity and tentatively propose amorphous carbon, and (3) significant differences with AGB star winds: presence of PAHs, absence of 'the' 13 micron band, and a lack of strong water bands. Conclusions: Dust condensation in RSG is found to experience a similar freeze-out process as in AGB stars. Together with the positive effect of the stellar luminosity on the mass-loss rate, this suggests that radiation pressure on dust grains is an important ingredient in the driving mechanism. Still, differences with AGB stars are manifold and thus the winds of RSG deserve separate studies.

The study of the atmosphere of red supergiant stars in general and of Betelgeuse (alpha Orionis) in particular is of prime importance to understand dust formation and how mass is lost to the interstellar medium in evolved massive stars. A molecular shell, the MOLsphere (Tsuji, 2000a), in the atmosphere of Betelgeuse has been proposed to account for the near- and mid-infrared spectroscopic observations of Betelgeuse. The goal is to further test this hypothesis and to identify some of the molecules in this MOLsphere. We report on measurements taken with the mid-infrared two-telescope beam combiner of the VLTI, MIDI, operated between 7.5 and 13.5 $\mu$m. The data are compared to a simple geometric model of a photosphere surrounded by a warm absorbing and emitting shell. Physical characteristics of the shell are derived: size, temperature and optical depth. The chemical constituents are determined with an analysis consistent with available infrared spectra and interferometric data. We are able to account for the measured optical depth of the shell in the N band, the ISO-SWS spectrum and K and L band interferometric data with a shell whose inner and outer radii are given by the above range and with the following species: H2O, SiO and Al2O3. These results confirm the MOLsphere model. We bring evidence for more constituents and for the presence of species participating in the formation of dust grains in the atmosphere of the star, i.e. well below the distance at which the dust shell is detected. We believe these results bring key elements to the understanding of mass loss in Betelgeuse and red supergiants in general and bring support to the dust-driven scenario.

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

We present the first mid-infrared spectra of SMP LMC 11 in the Large Magellanic Cloud. While this object resembles a planetary nebula in the optical, its infrared properties are more similar to an object in transition from the asymptotic giant branch to the planetary nebula phase. A warm dust continuum dominates the infrared spectrum. The peak emission corresponds to a mean dust temperature of 330 K. The spectrum shows overlapping molecular absorption bands from 12 to 17 um corresponding to acetylene and polyacetylenic chains and benzene. This is the first detection of C4H2, C6H2, C6H6 and other molecules in an extragalactic object. The infrared spectrum of SMP LMC 11 is similar in many ways to that of the pre-planetary nebula AFGL 618. The IRS spectrum shows little evidence of nitrogen-based molecules which are commonly seen in Galactic AGB stars. Polycyclic aromatic hydrocarbons are also absent from the spectrum. The detection of the [NeII] 12.8 um line in the infrared and other forbidden emission lines in the optical indicates that an ionized region is present.

We report on our first results from a mid-infrared spectroscopic study of ISM features in a sample of deeply obscured ULIRG nuclei using the InfraRed Spectrograph (IRS) on the Spitzer Space Telescope. The spectra are extremely rich and complex, revealing absorption features of both amorphous and crystalline silicates, aliphatic hydrocarbons, water ice and gas phase bands of hot CO and warm C_2H_2, HCN and CO_2. PAH emission bands were found to be generally weak and in some cases absent. The features are probing a dense and warm environment in which crystalline silicates and water ice are able to survive but volatile ices, commonly detected in Galactic dense molecular clouds, cannot. If powered largely by star formation, the stellar density and conditions of the gas and dust have to be extreme not to give rise to the commonly detected emission features associated with starburst.

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.

Silicates are an important component of interstellar dust and the structure of these grains -- amorphous versus crystalline -- is sensitive to the local physical conditions. We have studied the infrared spectra of a sample of ultra-luminous infrared galaxies. Here, we report the discovery of weak, narrow absorption features at 11, 16, 19, 23, and 28 microns, characteristic of crystalline silicates, superimposed on the broad absorption bands at 10 and 18 microns due to amorphous silicates in a subset of this sample. These features betray the presence of forsterite (Mg_2SiO_4), the magnesium-rich end member of the olivines. Previously, crystalline silicates have only been observed in circumstellar environments. The derived fraction of forsterite to amorphous silicates is typically 0.1 in these ULIRGs. This is much larger than the upper limit for this ratio in the interstellar medium of the Milky Way, 0.01. These results suggest that the timescale for injection of crystalline silicates into the ISM is short in a merger-driven starburst environment (e.g., as compared to the total time to dissipate the gas), pointing towards massive stars as a prominent source of crystalline silicates. Furthermore, amorphization due to cosmic rays, which is thought to be of prime importance for the local ISM, lags in vigorous starburst environments.

Analysis of high spectral resolution observations of the $\lambda$6614 diffuse interstellar band (DIB) line profile show systematic variations in the positions of the peaks in the substructure of the profile. These variations -- shown here for the first time -- can be understood most naturally in the framework of rotational contours of large molecules, where the variations are caused by changes in the rotational excitation temperature. We show that the rotational excitation temperature for the DIB carrier is likely significantly lower than the gas kinetic temperature -- indicating that for this particular DIB carrier angular momentum buildup is not very efficient.

We report the detection of strong absorption and weak emission features in the 4--27 micron Spitzer-IRS spectrum of the distant ultraluminous infrared galaxy (ULIRG) IRAS F00183--7111 (z=0.327). The absorption features of CO2 and CO gas, water ice, hydrocarbons and silicates are indicative of a strongly obscured (A[9.6]>=5.4; A[V]>=90) and complex line of sight through both hot diffuse ISM and shielded cold molecular clouds towards the nuclear power source. From the profile of the 4.67 micron CO fundamental vibration mode we deduce that the absorbing gas is dense (n~10^6 cm^-3) and warm (720 K) and has a CO column density of ~10^19.5 cm^-2, equivalent to N[H]~10^23.5 cm^-2. The high temperature and density, as well as the small infered size (<0.03pc), locates this absorbing gas close to the power source of this region. Weak emission features of molecular hydrogen, PAHs and Ne+, likely associated with star formation, are detected against the 9.7 micron silicate feature, indicating an origin away from the absorbing region. Based on the 11.2 micron PAH flux, we estimate the star formation component to be responsible for up to 30% of the IR luminosity of the system. While our mid-infrared spectrum shows no tell-tale signs of AGN activity, the similarities to the mid-infrared spectra of deeply obscured sources (e.g. NGC4418) and AGN hot dust (e.g. NGC1068), as well as evidence from other wavelength regions, suggest that the power source hiding behind the optically thick dust screen may well be a buried AGN.

We present submillimeter observations of rotational transitions of carbon monoxide from J = 2 -> 1 up to 7 -> 6 for a sample of Asymptotic Giant Branch stars and red supergiants. It is the first time that the high transitions J = 6 -> 5 and 7 -> 6 are included in such a study. With line radiative transfer calculations, we aim to determine the mass-loss history of these stars by fitting the CO line intensities. We find that the observed line intensities of the high transitions, including the J = 4 -> 3 transition, are significantly lower than the predicted values. We conclude that the physical structure of the outflow of Asymptotic Giant Branch stars is more complex than previously thought. In order to understand the observed line intensities and profiles, a physical structure with a variable mass-loss rate and/or a gradient in stochastic gas velocity is required. A case study of the AGB star WX Psc is performed. We find that the CO line strengths may be explained by variations in mass-loss on time scales similar to those observed in the separated arc-like structures observed around post-AGB stars. In addition, a gradient in the stochastic velocity may play a role. Until this has been sorted out fully, any mass loss determinations based upon single CO lines will remain suspect.

The composition of ice grains provides an important tool for the study of the molecular environment of star forming regions. Using ISAAC at the VLT to obtain spectra around 4.65 microns we have detected for the first time `XCN' and CO ice in an extragalactic environment: the nuclear region of the nearby dusty starburst/AGN galaxy NGC4945. The profile of the solid CO band reveals the importance of thermal processing of the ice while the prominence of the XCN band attests to the importance of energetic processing of the ice by FUV radiation and/or energetic particles. In analogy to the processing of ices by embedded protostars in our Galaxy, we attribute the processing of the ices in the center of NGC4945 to ongoing massive star formation. Our M-band spectrum also shows strong HI Pfund-beta and H2 0-0 S(9) line emission and gas phase CO absorption lines. The HI, H2, PAH, gas phase CO and the ices seem to be embedded in a rotating molecular disk which is undergoing vigorous star formation. Recently, strong OCN- absorption has been detected in the spectrum of the Galactic center star GC:IRS19. The most likely environment for the OCN- absorption is the strongly UV-exposed GC molecular ring. The presence of processed ice in the center of NGC4945 and our Galactic center leads us to believe that processed ice may be a common characteristic of dense molecular material in star forming galactic nuclei.

We present the Spectral Energy Distribution of HR 4049 based on literature data and new continuum measurements at 850 micron. The SED shows variable absorption in the UV, and a large IR excess, both caused by circumstellar dust. The shape of the IR excess from 1 micron all the way down to 850 micron can be nearly perfectly fitted with a single blackbody function at T ~ 1150 K or alternatively with a sum of blackbodies in a narrow temperature range. The energy emitted in this IR continuum radiation is about one-third of the stellar luminosity. We show that this blackbody radiation must be due to the presence of a circumbinary disk with a large height. This disk must also be gas-rich, in agreement with the observations of molecular bands in the ISO-SWS spectrum. We present two possible scenario's for explaining the shape and the intensity of the IR excess. The first scenario involves large grains (a >1 mm) that each radiate like a blackbody. The second scenario argues that the blackbody radiation is due to a very optically thick circumbinary disk. We investigate if such a disk would indeed produce blackbody radiation by presenting results from radiative transfer calculations. We further quantify the properties of such a disk and its stability in the framework of (hydro)dynamics, grain settling, radiation pressure and grain drift. The virtues and shortcomings of both models for the origin of the IR blackbody are discussed by contrasting them with other observations and assessing them in the framework of (binary) (post-)AGB evolution.

The time variation in the water-vapour bands in oxygen-rich Mira variables has been investigated using multi-epoch ISO/SWS spectra of four Mira variables in the 2.5-4.0 micron region. All four stars show H2O bands in absorption around minimum in the visual light curve. At maximum, H2O emission features appear in the ~3.5-4.0 micronm region, while the features at shorter wavelengths remain in absorption. These H2O bands in the 2.5-4.0 micron region originate from the extended atmosphere. The analysis has been carried out with a disk shape, slab geometry model. The observed H2O bands are reproduced by two layers; a `hot' layer with an excitation temperature of 2000 K and a `cool' layer with an excitation temperature of 1000-1400 K. The radii of the `hot' layer (R_hot) are ~1 R_* at visual minimum and 2 R_* at maximum, where R_* is a radius of background source of the model. The time variation of R_hot/R_* from 1 to 2 is attributed to the actual variation in the radius of the H2O layer. A high H2O density shell occurs near the surface of the star around minimum, and moves out with the stellar pulsation. This shell gradually fades away after maximum, and a new high H2O density shell is formed in the inner region again at the next minimum. Due to large optical depth of H2O, the near-infrared variability is dominated by the H2O layer, and the L'-band flux correlates with the area of the H2O shell. The infrared molecular bands trace the structure of the extended atmosphere and impose appreciable effects on near-infrared light curve of Mira variables.

We analyze ISO archive data of the nearby bright emission-line star 51 Oph, previously classified as a proto-planetary system similar to beta Pic. The infrared spectrum reveals the presence of gas-phase emission bands of hot (approx 850 K) CO, CO_2, H_2O and NO. In addition to this, partially crystalline silicate dust is present. The solid-state bands and the energy distribution are indicative of dust that has formed recently, rather than of debris dust. The presence of hot molecular gas and the composition of the circumstellar dust are highly unusual for a young star and are reminiscent of what is found around evolved (AGB) stars, although we exclude the possibility of 51 Oph belonging to this group. We suggest several explanations for the nature of 51 Oph, including a recent episode of mass loss from a Be star, and the recent destruction of a planet-sized body around a young star.