The Orion Molecular Cloud is the nearest massive-star forming region. Massive stars have profound effects on their environment due to their strong radiation fields and stellar winds. Velocity-resolved observations of the [CII] $158\,\mu\mathrm{m}$ fine-structure line allow us to study the kinematics of UV-illuminated gas. Here, we present a square-degree-sized map of [CII] emission from the Orion Nebula complex obtained by the upGREAT instrument onboard SOFIA, covering the entire Orion Nebula (M42) plus M43 and the nebulae NGC 1973, 1975, and 1977. We compare the stellar characteristics of these three regions with the kinematics of the expanding bubbles surrounding them. The bubble blown by the O7V star $\theta^1$ Ori C in the Orion Nebula expands rapidly, at $13\,\mathrm{km\,s^{-1}}$. Simple analytical models reproduce the characteristics of the hot interior gas and the neutral shell of this wind-blown bubble and give us an estimate of the expansion time of $0.2\,\mathrm{Myr}$. M43 with the B0.5V star NU Ori also exhibits an expanding bubble structure, with an expansion velocity of $6\,\mathrm{km\,s^{-1}}$. Comparison with analytical models for the pressure-driven expansion of H\,{\sc ii} regions gives an age estimate of $0.02\,\mathrm{Myr}$. The bubble surrounding NGC 1973, 1975, and 1977 with the central B1V star 42 Orionis expands at $1.5\,\mathrm{km\,s^{-1}}$, likely due to the over-pressurized ionized gas as in the case of M43. We derive an age of $0.4\,\mathrm{Myr}$ for this structure. We conclude that the bubble of the Orion Nebula is driven by the mechanical energy input by the strong stellar wind from $\theta^1$ Ori C, while the bubbles associated with M43 and NGC 1977 are caused by the thermal expansion of the gas ionized by their central later-type massive stars.

Strong winds and ultraviolet (UV) radiation from O-type stars disrupt and ionize their molecular core birthplaces, sweeping up material into parsec-size shells. Owing to dissociation by starlight, the thinnest shells are expected to host low molecular abundances and therefore little star formation. Here, we expand previous maps taken with the IRAM 30m telescope and present square-degree 12CO and 13CO (J=2-1) maps of the wind-driven "Veil bubble'' that surrounds the Trapezium cluster and its natal Orion molecular core (OMC). Although widespread and extended CO emission is largely absent from the Veil, we show that several CO "globules'' exist and are embedded in the [CII]158um-bright shell that confines the bubble. This includes the first detection of quiescent CO at negative LSR velocities in Orion. Given the harsh UV irradiation conditions in this translucent material, the detection of CO globules is surprising. These globules are small (R=7,100 AU), not massive (M=0.3M_Sun), and are moderately dense: n_ H=4x10^4 cm^-3 (median values). They are confined by the external pressure of the shell, P_ext/k~10^7 cm^-3 K, and are likely magnetically supported. They are either transient objects formed by instabilities or have detached from pre-existing molecular structures, sculpted by the passing shock associated with the expanding shell and by UV radiation from the Trapezium. Some represent the first stages in the formation of small pillars, others of isolated small globules. Although their masses do not suggest they will form stars, one globule matches the position of a known YSO. The lack of extended CO in the "Veil shell'' demonstrates that feedback from massive stars expels, agitates, and reprocesses most of the disrupted molecular cloud gas, thereby limiting the star-formation rate in the region. The presence of globules is a result of this feedback.

We have observed the massive protostar AFGL 2136 IRS 1 in multiple wavelength windows in the near-to-mid-infrared at high ($\sim3$ km s$^{-1}$) spectral resolution using VLT+CRIRES, SOFIA+EXES, and Gemini North+TEXES. There is an abundance of H$_2$O absorption lines from the $\nu_1$ and $\nu_3$ vibrational bands at 2.7 $\mu$m, from the $\nu_2$ vibrational band at 6.1 $\mu$m, and from pure rotational transitions near 10-13 $\mu$m. Analysis of state-specific column densities derived from the resolved absorption features reveals that an isothermal absorbing slab model is incapable of explaining the relative depths of different absorption features. In particular, the strongest absorption features are much weaker than expected, indicating optical depth effects resulting from the absorbing gas being well-mixed with the warm dust that serves as the "background" continuum source at all observed wavelengths. The velocity at which the strongest H$_2$O absorption occurs coincides with the velocity centroid along the minor axis of the compact disk in Keplerian rotation recently observed in H$_2$O emission with ALMA. We postulate that the warm regions of this dust disk dominate the continuum emission at near-to-mid infrared wavelengths, and that H$_2$O and several other molecules observed in absorption are probing this disk. Absorption line profiles are not symmetric, possibly indicating that the warm dust in the disk that produces the infrared continuum has a non-uniform distribution similar to the substructure observed in 1.3 mm continuum emission.

We report on a common fragment ion formed during the electron-ionization-induced fragmentation of three different three-ring polycyclic aromatic hydrocarbons (PAHs), fluorene (C$_{13}$H$_{10}$), 9,10-dihydrophenanthrene (C$_{14}$H$_{12}$), and 9,10-dihydroanthracene (C$_{14}$H$_{12}$). The infrared spectra of the mass-isolated product ions with $m/z=165$ were obtained in a Fourier transform ion cyclotron resonance mass spectrometer whose cell was placed inside the optical cavity of an infrared free-electron laser, thus providing the high photon fluence required for efficient infrared multiple-photon dissociation. The infrared spectra of the $m/z=165$ species generated from the three different precursors were found to be similar, suggesting the formation of a single C$_{13}$H$_{9}^+$ isomer. Theoretical calculations using density functional theory (DFT) revealed the fragment's identity as the closed-shell fluorenyl cation. Decomposition pathways from each parent precursor to the fluorenyl ion are proposed on the basis of DFT calculations. The identification of a single fragmentation product from three different PAHs supports the notion of the existence of common decomposition pathways of PAHs in general and can aid in understanding the fragmentation chemistry of astronomical PAH species.

Besides buckminsterfullerene (C60), other fullerenes and their derivatives may also reside in space. In this work, we study the formation and photo-dissociation processes of astronomically relevant fullerene/anthracene (C14H10) cluster cations in the gas phase. Experiments are carried out using a quadrupole ion trap (QIT) in combination with time-of-flight (TOF) mass spectrometry. The results show that fullerene (C60, and C70)/anthracene (i.e., [(C14H10)nC60]+ and [(C14H10)nC70]+), fullerene (C56 and C58)/anthracene (i.e., [(C14H10)nC56]+ and [(C14H10)nC58]+) and fullerene (C66 and C68)/anthracene (i.e., [(C14H10)nC66]+ and [(C14H10)nC68]+) cluster cations, are formed in the gas phase through an ion-molecule reaction pathway. With irradiation, all the fullerene/anthracene cluster cations dissociate into mono$-$anthracene and fullerene species without dehydrogenation. The structure of newly formed fullerene/anthracene cluster cations and the bonding energy for these reaction pathways are investigated with quantum chemistry calculations. Our results provide a growth route towards large fullerene derivatives in a bottom-up process and insight in their photo-evolution behavior in the ISM, and clearly, when conditions are favorable, fullerene/PAH clusters can form efficiently. In addition, these clusters (from 80 to 154 atoms or ~ 2 nm in size) offer a good model for understanding the physical-chemical processes involved in the formation and evolution of carbon dust grains in space, and provide candidates of interest for the DIBs that could motivate spectroscopic studies.

We study the impact of active galactic nuclei (AGN) on the CO Spectral Line Energy Distribution (SLED) of high-$z$ galaxies. In particular, we want to assess if the CO SLED can be used as a probe of AGN activity. To this purpose, we develop a semi-analytical model that takes into account the radiative transfer and the clumpy structure of giant molecular clouds where the CO lines are excited, their distribution in the galaxy disk, and the torus obscuration of the AGN radiation. We study the joint effect on the CO SLED excitation of (i) the X-ray luminosity from the AGN ($L_{X}$), (ii) the size of the molecular disk, (iii) the inclination angle ($\Omega$) of the torus with respect to the molecular disk, and (iv) the GMC mean density. We also discuss the possibility of an enhanced Cosmic Ray Ionization Rate (CRIR). We find that the X-ray Dominated Region (XDR) generated by the AGN in every case enhances the CO SLED for $J>5$, with increasing excitation of high-$J$ CO lines for increasing X-ray luminosity. Because high-$z$ galaxies are compact, the XDR region typically encloses the whole disk, thus its effect can be more important with respect to lower redshift objects. The impact of the XDR can be disentangled from an enhanced CRIR either if $L_X>10^{44} \rm \,erg\, s^{-1}$, or if $\Omega \geq 60^{\circ}$. We finally provide predictions on the CO(7-6)/[CII] and CO(17-16)/[CII] ratios as a function of $L_X$, which can be relevant for ALMA follow up of galaxies and quasars previously detected in [CII].

(abridged) Radio recombination lines (RRLs) at frequencies $\nu$ < 250 MHz trace the cold, diffuse phase of the ISM. Next generation low frequency interferometers, such as LOFAR, MWA and the future SKA, with unprecedented sensitivity, resolution, and large fractional bandwidths, are enabling the exploration of the extragalactic RRL universe. We observed the radio quasar 3C 190 (z~1.2) with the LOFAR HBA. In reducing this data for spectroscopic analysis, we have placed special emphasis on bandpass calibration. We devised cross-correlation techniques to significantly identify the presence of RRLs in a low frequency spectrum. We demonstrate the utility of this method by applying it to existing low-frequency spectra of Cassiopeia A and M 82, and to the new observations of 3C 190. RRLs have been detected in the foreground of 3C 190 at z = 1.12355 (assuming a carbon origin), owing to the first detection of RRLs outside of the local universe (first reported in Emig et al. 2019). Towards the Galactic supernova remnant Cas A, we uncover three new detections: (1) C$\epsilon$-transitions ($\Delta$n = 5) for the first time at low radio frequencies, (2) H$\alpha$-transitions at 64 MHz with a FWHM of 3.1 km/s, the most narrow and one of the lowest frequency detections of hydrogen to date, and (3) C$\alpha$ at v$_{LSR}$ = 0 km/s in the frequency range 55-78 MHz for the first time. Additionally we recover C$\alpha$, C$\beta$, C$\gamma$, and C$\delta$ from the -47 km/s and -38 km/s components. In the nearby starburst galaxy, M 82, we do not find a significant feature. Our current searches for RRLs in LOFAR observations are limited to narrow (< 100 km/s) features, owing to the relatively small number of channels available for continuum estimation. Future strategies making use of larger contiguous frequency coverage would aid calibration to deeper sensitivities and broader features.

The hydrogen fluoride (HF) molecule is seen in absorption in the interstellar medium (ISM) along many lines of sight. Surprisingly, it is observed in emission toward the Orion Bar, which is an interface between the ionized region around the Orion Trapezium stars and the Orion molecular cloud. We aim to understand the origin of HF emission in the Orion Bar by comparing its spatial distribution with other tracers. We examine three mechanisms to explain the HF emission: thermal excitation, radiative dust pumping, and chemical pumping. We used a Herschel/HIFI strip map of the HF J = 1-0 line, covering 0.5' by 1.5' that is oriented perpendicular to the Orion Bar. We used the RADEX non-local thermodynamic equilibrium (non-LTE) code to construct the HF column density map. We use the Meudon PDR code to explain the morphology of HF. The bulk of the HF emission at 10 km s^-1 emerges from the CO-dark molecular gas that separates the ionization front from the molecular gas that is deeper in the Orion Bar. The excitation of HF is caused mainly by collisions with H2 at a density of 10^5 cm^-3 together with a small contribution of electrons in the interclump gas of the Orion Bar. Infrared pumping and chemical pumping are not important. We conclude that the HF J = 1-0 line traces CO-dark molecular gas. Similarly, bright photodissociation regions associated with massive star formation may be responsible for the HF emission observed toward active galactic nuclei.

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.

[abridged]CO rovibrational lines are efficient probes of warm molecular gas and can give unique insights into the inner 10 AU of proto-planetary disks. Recent studies have found a relation between the ratio of lines originating from the second and first vibrationally excited state, denoted as $v2/v1$, and the emitting radius of CO. In disks around Herbig Ae stars the vibrational excitation is low when CO lines come from close to the star, and high when lines only probe gas at large radii (more than 5 AU). We aim to find explanations for the observed trends between CO vibrational ratio, emitting radii, and NIR excess, and identify their implications in terms of the physical and chemical structure of inner disks around Herbig stars. Slab models and full disk thermo chemical models are calculated. Simulated observations from the models are directly compared to the data. Broad CO lines with low vibrational ratios are best explained by a warm (400-1300 K) inner disk surface with gas-to-dust ratios below 1000; no CO is detected within/at the inner dust rim, due to dissociation at high temperatures. In contrast, explaining the narrow lines with high vibrational ratios requires an inner cavity of a least 5 AU in both dust and gas, followed by a cool (100-300 K) molecular gas reservoir with gas-to-dust ratios greater than 10000 at the cavity wall. In all cases the CO gas must be close to thermalization with the dust. The high gas-to-dust ratios needed to explain high $v2/v1$ in narrow CO lines for a subset of group I disks can naturally be interpreted as due to the dust traps that have been proposed to explain millimeter dust cavities. The broad lines seen in most group II objects indicate a very flat disk in addition to inner disk substructures within 10 AU that can be related to the substructures recently observed with ALMA.

The rich molecular structures of polycyclic aromatic hydrocarbons -- essentially planar flakes of fused benzene rings -- and their fullerene cousins are revealed through their vibrational and electronic spectra.

Aims. In this work we determine the effects of anharmonicity on the mid-infrared spectra of the linear polycyclic aromatic hydrocarbons (PAHs) naphthalene, anthracene, tetracene and pentacene recorded using the free electron laser FELIX. Methods. Comparison of experimental spectra obtained under supersonic jet conditions with theoretically predicted spectra was used to show if anharmonicity explicitly needs to be taken into account. Results. Anharmonic spectra obtained using second-order vibrational perturbation theory agree on average within 0.5% of the experimental frequencies. Importantly, they confirm the presence of combination bands with appreciable intensity in the 5-6 $\mu$m region. These combination bands contain a significant fraction of the IR absorption, which scales linearly with the size of the PAH. Detection and assignment of the combination bands are a preliminary indication of the accuracy of far-IR modes in our anharmonic theoretical spectra. Detailed analysis of the periphery-sensitive CH out-of-plane band of naphthalene reveals that there is still room for improvement of the VPT2 approach. In addition, the implications of our findings for the analysis of the aromatic infrared bands are discussed.

We present a study on the prospects of observing carbon, sulfur, and other lower abundance elements (namely Al, Ca, Ti and Ni) present in the interstellar medium using future X-ray instruments. We focus in particular on the detection and characterization of interstellar dust along the lines of sight. We compare the simulated data with different sets of dust aggregates, either obtained from past literature or measured by us using the SOLEIL-LUCIA synchrotron beamline. Extinction by interstellar grains induces modulations of a given photolelectric edge, which can be in principle traced back to the chemistry of the absorbing grains. We simulated data of instruments with characteristics of resolution and sensitivity of the current Athena, XRISM and Arcus concepts. In the relatively near future, the depletion and abundances of the elements under study will be determined with confidence. In the case of carbon and sulfur, the characterization of the chemistry of the absorbing dust will be also determined, depending on the dominant compound. For aluminum and calcium, despite the large depletion in the interstellar medium and the prominent dust absorption, in many cases the edge feature may not be changing significantly with the change of chemistry in the Al$^-$ or Ca$^-$ bearing compounds. The exinction signature of large grains may be detected and modeled, allowing a test on different grain size distributions for these elements. The low cosmic abundance of Ti and Ni will not allow us a detailed study of the edge features.

Dense, fast-moving ejecta knots in supernova remnants are prime sites for molecule and dust formation. We present SOFIA far-IR spectrometer FIFI-LS observations of CO-rich knots in Cas A which cover a ~1 square arc minute area of the northern shell, in the [O III] 52 and 88 micron and [O I] 63 micron lines. The FIFI-LS spectra reveal that the line profiles of [O III] and [O I] are similar to those of the Herschel PACS [O III] and CO lines. We find that the [O III] maps show very different morphology than the [O I] map. The [O III] maps reveal diffuse, large-scale structures and the ratio of the two [O III] lines imply the presence of gas with a range of density 500 - 10,000 per cm^3 within the mapped region. In contrast, the [O I] map shows bright emission associated with the dense CO-rich knots. The 63 micron [O I] line traces cooled, dense post-shocked gas of ejecta. We find that IR-dominated [O III] emission is from post-shocked gas based on its morphology, high column density, and velocity profile. We describe multi-phase ejecta knots, a lifetime of clumps, and survival of dust in the young supernova remnants.

We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $\mu$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a velocity resolution of 0.17 km s$^{-1}$. The observations cover a 0.25$^\circ$ by 0.28$^\circ$ area with central position {\it l} = 297.34$^\circ$, {\it b} = -0.60$^\circ$. The kinematics show that bright [C\,{\sc ii}] structures are spatially and spectrally correlated with the surfaces of CO clouds, tracing the photodissociation region and ionization front of each molecular cloud. Along 7 lines of sight that traverse Tr 14 into the dark ridge to the southwest, we find that the [C\,{\sc ii}] luminosity from the HII region is 3.7 times that from the PDR. In same los we find in the PDRs an average ratio of 1:4.1:5.6 for the mass in atomic gas:dark-CO gas: molecular gas traced by CO. Comparing multiple gas tracers including HI 21cm, [C\,{\sc ii}], CO, and radio recombination lines, we find that the HII regions of the Carina Nebula Complex are well-described as HII regions with one-side freely expanding towards us, consistent with the champagne model of ionized gas evolution. The dispersal of the GMC in this region is dominated by EUV photoevaporation; the dispersal timescale is 20-30 Myr.

The composition and properties of interstellar silicate dust are not well understood. In X-rays, interstellar dust can be studied in detail by making use of the fine structure features in the Si K-edge. The features in the Si K-edge offer a range of possibilities to study silicon-bearing dust, such as investigating the crystallinity, abundance, and the chemical composition along a given line of sight. We present newly acquired laboratory measurements of the silicon K-edge of several silicate-compounds that complement our measurements from our earlier pilot study. The resulting dust extinction profiles serve as templates for the interstellar extinction that we observe. The extinction profiles were used to model the interstellar dust in the dense environments of the Galaxy. The laboratory measurements, taken at the Soleil synchrotron facility in Paris, were adapted for astrophysical data analysis and implemented in the SPEX spectral fitting program. The models were used to fit the spectra of nine low-mass X-ray binaries located in the Galactic center neighborhood in order to determine the dust properties along those lines of sight. Most lines of sight can be fit well by amorphous olivine. We also established upper limits on the amount of crystalline material that the modeling allows. We obtained values of the total silicon abundance, silicon dust abundance, and depletion along each of the sightlines. We find a possible gradient of $0.06\pm0.02$ dex/kpc for the total silicon abundance versus the Galactocentric distance. We do not find a relation between the depletion and the extinction along the line of sight.

Context. The combined use of carbon radio recombination lines (CRRLs) and the 158 $\mu$m-[CII] line is a powerful tool for the study of the energetics and physical conditions (e.g., temperature and density) of photodissociation regions (PDRs). However, there are few observational studies that exploit this synergy. Aims. Here we explore the relation between CRRLs and the 158 $\mu$m-[CII] line in light of new observations and models. Methods. We present new and existing observations of CRRLs in the frequency range 0.15--230 GHz with ALMA, VLA, the GBT, Effelsberg 100m, and LOFAR towards Orion~A (M42). We complement these observations with SOFIA observations of the 158 $\mu$m-[CII] line. We studied two PDRs: the foreground atomic gas, known as the Veil, and the dense PDR between the HII region and the background molecular cloud. Results. In the Veil we are able to determine the gas temperature and electron density, which we use to measure the ionization parameter and the photoelectric heating efficiency. In the dense PDR, we are able to identify a layered PDR structure at the surface of the molecular cloud to the south of the Trapezium cluster. There we find that the radio lines trace the colder portion of the ionized carbon layer, the C$^{+}$/C/CO interface. By modeling the emission of the $158$~$\mu$m-[CII] line and CRRLs as arising from a PDR we derive a thermal pressure $>5\times10^{7}$ K cm$^{-3}$ and a radiation field $G_{0}\approx10^{5}$ close to the Trapezium. Conclusions. This work provides additional observational support for the use of CRRLs and the 158 $\mu$m-[CII] line as complementary tools to study dense and diffuse PDRs, and highlights the usefulness of CRRLs as probes of the C$^{+}$/C/CO interface.

A significant fraction of the molecular gas in star-forming regions is irradiated by stellar UV photons. In these environments, the electron density (n_e) plays a critical role in the gas dynamics, chemistry, and collisional excitation of certain molecules. We determine n_e in the prototypical strongly irradiated photodissociation region (PDR), the Orion Bar, from the detection of new millimeter-wave carbon recombination lines (mmCRLs) and existing far-IR [13CII] hyperfine line observations. We detect 12 mmCRLs (including alpha, beta, and gamma transitions) observed with the IRAM 30m telescope, at ~25'' angular resolution, toward the H/H2 dissociation front (DF) of the Bar. We also present a mmCRL emission cut across the PDR. These lines trace the C+/C/CO gas transition layer. As the much lower frequency carbon radio recombination lines, mmCRLs arise from neutral PDR gas and not from ionized gas in the adjacent HII region. This is readily seen from their narrow line profiles (dv=2.6+/-0.4 km/s) and line peak LSR velocities (v_LSR=+10.7+/-0.2 km/s). Optically thin [13CII] hyperfine lines and molecular lines - emitted close to the DF by trace species such as reactive ions CO+ and HOC+ - show the same line profiles. We use non-LTE excitation models of [13CII] and mmCRLs and derive n_e = 60-100 cm^-3 and T_e = 500-600 K toward the DF. The inferred electron densities are high, up to an order of magnitude higher than previously thought. They provide a lower limit to the gas thermal pressure at the PDR edge without using molecular tracers. We obtain P_th > (2-4)x10^8 cm^-3 K assuming that the electron abundance is equal or lower than the gas-phase elemental abundance of carbon. Such elevated thermal pressures leave little room for magnetic pressure support and agree with a scenario in which the PDR photoevaporates.

Massive stars inject mechanical and radiative energy into the surrounding environment, which stirs it up, heats the gas, produces cloud and intercloud phases in the interstellar medium, and disrupts molecular clouds (the birth sites of new stars). Stellar winds, supernova explosions and ionization by ultraviolet photons control the lifetimes of molecular clouds. Theoretical studies predict that momentum injection by radiation should dominate that by stellar winds, but this has been difficult to assess observationally. Velocity-resolved large-scale images in the fine-structure line of ionized carbon ([C II]) provide an observational diagnostic for the radiative energy input and the dynamics of the interstellar medium around massive stars. Here we report observations of a one-square-degree region (about 7 parsecs in diameter) of Orion molecular core -- the region nearest to Earth that exhibits massive-star formation -- at a resolution of 16 arcseconds (0.03 parsecs) in the [C II] line at 1.9 terahertz (158 micrometres). The results reveal that the stellar wind originating from the massive star ${\theta}^{1}$ Orionis C has swept up the surrounding material to create a bubble roughly four parsecs in diameter with a 2,600-solar-mass shell, which is expanding at 13 kilometres per second. This finding demonstrates that the mechanical energy from the stellar wind is converted very efficiently into kinetic energy of the shell and causes more disruption of the Orion molecular core 1 than do photo-ionization and evaporation or future supernova explosions.

The fullerene C$_{60}$, one of the largest molecules identified in the interstellar medium (ISM), has been proposed to form top-down through the photo-chemical processing of large (more than 60 C-atoms) polycyclic aromatic hydrocarbon (PAH) molecules. In this article, we focus on the opposite process, investigating the possibility that fullerenes form from small PAHs, in which bowl-forming plays a central role. We combine laboratory experiments and quantum chemical calculations to study the formation of larger PAHs from charged fluorene clusters. The experiments show that with visible laser irradiation, the fluorene dimer cation - [C$_{13}$H$_{9}$$-$C$_{13}$H$_{9}$]$^+$ - and the fluorene trimer cation - [C$_{13}$H$_{9}$$-$C$_{13}$H$_{8}$$-$C$_{13}$H$_{9}$]$^+$ - undergo photo-dehydrogenation and photo-isomerization resulting in bowl structured aromatic cluster-ions, C$_{26}$H$_{12}$$^+$ and C$_{39}$H$_{20}$$^+$, respectively. To study the details of this chemical process, we employ quantum chemistry that allows us to determine the structures of the newly formed cluster-ions, to calculate the hydrogen loss dissociation energies, and to derive the underlying reaction pathways. These results demonstrate that smaller PAH clusters (with less than 60 C-atoms) can convert to larger bowled geometries that might act as building blocks for fullerenes, as the bowl-forming mechanism greatly facilitates the conversion from dehydrogenated PAHs to cages. Moreover, the bowl-forming induces a permanent dipole moment that - in principle - allows to search for such species using radio astronomy.

Recombination lines involving high principal quantum numbers populate the radio spectrum in large numbers. Low-frequency (<1 GHz) observations of radio recombination lines (RRLs) primarily from carbon and hydrogen offer a new, if not unique, way to probe cold, largely atomic gas and warm, ionised gas in other galaxies. Furthermore, RRLs can be used to determine the physical state of the emitting regions, such as temperature and density. These properties make RRLs, potentially, a powerful tool of extragalactic ISM physics. At low radio frequencies, its conceivable to detect RRLs out to cosmological distances when illuminated by a strong radio continuum. However, they are extremely faint (tau ~ 1e-3 -- 1e-4) and have so far eluded detection outside of the local universe. With LOFAR observations of the radio quasar 3C 190 (z=1.1946), we aim to demonstrate that the ISM can be explored out to great distances through low-frequency RRLs. We report the detection of RRLs in the frequency range 112--163 MHz in the spectrum of 3C 190. Stacking 13 a-transitions with principal quantum numbers n=266-301, a peak 6sigma feature of optical depth, tau(peak) = (1.0 +\- 0.2) x 1e-3 and FWHM = 31.2 +/- 8.3 km/s was found at z=1.124. This corresponds to a velocity offset of -9965 km/s with respect to the systemic redshift of 3C 190. We consider three interpretations of the origin of the RRL emission: an intervening dwarf-like galaxy, an AGN-driven outflow, and the inter-galactic medium. We argue that the RRLs most likely originate in a dwarf-like galaxy (M ~ 1e9 Msun) along the line of sight, although we cannot rule out an AGN-driven outflow. We do find the RRLs to be inconsistent with an inter-galactic medium origin. With this detection, we have opened up a new way to study the physical properties of cool, diffuse gas out to cosmological distances.

Sulfur has been observed to be severely depleted in dense clouds leading to uncertainty in the molecules that contain it and the chemistry behind their evolution. Here, we aim to shed light on the sulfur chemistry in young stellar objects (YSOs) by using high-resolution infrared spectroscopy of absorption by the $\nu_3$ rovibrational band of SO$_2$ obtained with the Echelon-Cross-Echelle Spectrograph on the Stratospheric Observatory for Infrared Astronomy. Using local thermodynamic equilibrium models we derive physical parameters for the SO$_2$ gas in the massive YSO MonR2 IRS3. This yields a SO$_2$/$\mathrm{H}$ abundance lower limit of $5.6\pm0.5\times10^{-7}$, or $>\!4\%$ of the cosmic sulfur budget, and an intrinsic line width (Doppler parameter) of $b<3.20\;\mathrm{km\;s}^{-1}$. The small line widths and high temperature ($T_\mathrm{ex}=234\pm15\;\mathrm{K}$) locate the gas in a relatively quiescent region near the YSO, presumably in the hot core where ices have evaporated. This sublimation unlocks a volatile sulfur reservoir (e.g., sulfur allotropes as detected abundantly in comet 67P/Churyumov--Gerasimenko), which is followed by SO$_2$ formation by warm, dense gas-phase chemistry. The narrowness of the lines makes formation of SO$_2$ from sulfur sputtered off grains in shocks less likely toward MonR2 IRS3.

We have performed a 5-8 $\mu$m spectral line survey of the hot molecular core associated with the massive protostar AFGL 2591, using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We have supplemented these data with a ground based study in the atmospheric M band around 4.5 $\mu$m using the iSHELL instrument on the Infrared Telescope Facility (IRTF), and the full N band window from 8-13 $\mu$m using the Texas Echelon Cross Echelle Spectrograph (TEXES) on the IRTF. Here we present the first detection of ro-vibrational transitions of CS in this source. The absorption lines are centred on average around -10 kms$^{-1}$ and the line widths of CS compare well with the hot component of $^{13}$CO (around 10 kms$^{-1}$). Temperatures for CS, hot $^{13}$CO and $^{12}$CO v=1-2 agree well and are around 700 K. We derive a CS abundance of 8$\times$10$^{-3}$ and 2$\times$10$^{-6}$ with respect to CO and H$_2$ respectively. This enhanced CS abundance with respect to the surrounding cloud (1$\times$10$^{-8}$) may reflect sublimation of H$_2$S ice followed by gas-phase reactions to form CS. Transitions are in LTE and we derive a density of $>$10$^7$ cm$^{-3}$, which corresponds to an absorbing region of $<$0.04$''$. EXES observations of CS are likely to probe deeply into the hot core, to the base of the outflow. Submillimeter and infrared observations trace different components of the hot core as revealed by the difference in systemic velocities, line widths and temperatures, as well as the CS abundance.

The infrared (IR) emission of polycyclic aromatic hydrocarbons (PAHs) permeates our universe; astronomers have detected the IR signatures of PAHs around many interstellar objects. The IR emission of interstellar PAHs differs from their emission as seen under conditions on Earth, as they emit through a collisionless cascade down through their excited vibrational states from high internal energies. The difficulty in reproducing interstellar conditions in the laboratory results in a reliance on theoretical techniques. However, the size and complexity of PAHs requires careful consideration when producing the theoretical spectra. In this work we outline the theoretical methods necessary to lead to a fully theoretical IR cascade spectra of PAHs including: an anharmonic second order vibrational perturbation theory (VPT2) treatment; the inclusion of Fermi resonances through polyads; and the calculation of anharmonic temperature band shifts and broadenings (including resonances) through a Wang--Landau approach. We also suggest a simplified scheme to calculate vibrational emission spectra that retains the essential characteristics of the full IR cascade treatment and can directly transform low temperature absorption spectra in IR cascade spectra. Additionally we show that past astronomical models were in error in assuming a 15 cm$^{-1}$ correction was needed to account for anharmonic emission effects.

We present mid-infrared (MIR) spectra of HII regions within star-forming galaxies M83 and M33. Their emission features are compared with Galactic and extragalactic HII regions, HII-type galaxies, starburst galaxies, and Seyfert/LINER type galaxies. Our main results are as follows: (i) the M33 and M83 HII regions lie in between Seyfert/LINER galaxies and HII-type galaxies in the 7.7/11.3 - 6.2/11.3 plane, while the different sub-samples exhibiting different 7.7/6.2 ratios; (ii) Using the NASA Ames PAH IR Spectroscopic database, we demonstrate that the 6.2/7.7 ratio does not effectively track PAH size, but the 11.3/3.3 PAH ratio does; (iii) variations on the 17 $\mu$m PAH band depends on object type; however, there is no dependence on metallicity for both extragalactic HII regions and galaxies; (iv) the PAH/VSG intensity ratio decreases with the hardness of the radiation field and galactocentric radius (Rg), yet the ionization alone cannot account for the variation seen in all of our sources; (v) the relative strength of PAH features does not change significantly with increasing radiation hardness, as measured through the [NeIII]/[NeII] ratio and the ionization index; (vi) We present PAH SFR calibrations based on the tight correlation between the 6.2, 7.7, and 11.3 $\mu$m PAH luminosities with the 24 $\mu$m luminosity and the combination of the 24 $\mu$m and H$\alpha$ luminosity; (vii) Based on the total luminosity from PAH and FIR emission, we argue that extragalactic HII regions are more suitable templates in modeling and interpreting the large scale properties of galaxies compared to Galactic HII regions.

Aims. Infrared (IR) spectroscopy is a powerful tool to study molecules in space. A key issue in such analyses is understanding the effect that temperature and anharmonicity have on different vibrational bands, and thus interpreting the IR spectra for molecules under various conditions. Methods. We combined second order vibrational perturbation theory and the Wang-Landau random walk technique to produce accurate IR spectra of highly excited PAHs. We fully incorporated anharmonic effects, such as resonances, overtones, combination bands, and temperature effects. Results. The results are validated against experimental results for the pyrene molecule (C16H10). In terms of positions, widths, and relative intensities of the vibrational bands, our calculated spectra are in excellent agreement with gas-phase experimental data.

In this work, we study the photodissociation processes of small PAH clusters (e.g., pyrene clusters). The experiments are carried out using a quadrupole ion trap in combination with time-of-flight (QIT-TOF) mass spectrometry. The results show that pyrene clusters are converted into larger PAHs under the influence of a strong radiation field. Specifically, pyrene dimer cations (e.g., [C$_{16}$H$_{10}$$-$C$_{16}$H$_{9}$]$^+$ or C$_{32}$H$_{19}$$^+$), will photo-dehydrogenate and photo-isomerize to fully aromatic cations (PAHs) (e.g., C$_{32}$H$_{16}$$^+$) with laser irradiation. The structure of new formed PAHs and the dissociation energy for these reaction pathways are investigated with quantum chemical calculations. These studies provide a novel efficient evolution routes for the formation of large PAHs in the interstellar medium (ISM) in a bottom-up process that will counteract the top-down conversion of large PAHs into rings and chains, and provide a reservoir of large PAHs that can be converted into C$_{60}$ and other fullerenes and large carbon cages.

We present the first infrared (IR) gas phase spectrum of a large and astronomically relevant PAH cation (C$_{42}$H$_{18}$$^+$, HBC$^+$) and its di-cation (C$_{42}$H$_{18}$$^{2+}$, HBC$^{2+}$). The spectra are recorded via infrared multi-photon dissociation (IRMPD) spectroscopy of ions stored in a quadrupole ion trap, using the intense infrared radiation of a free electron laser in the 530$-$1800 cm$^{-1}$ (5.6$-$18.9 $\mu$m) range. HBC$^{+}$ shows main intense absorption peaks at 762 (13.12), 1060 (9.43), 1192 (8.39), 1280 (7.81), 1379 (7.25) and 1530 (6.54) cm$^{-1}$($\mu$m), in good agreement with DFT calculations that after scaling to take the anharmonicities effect into account. HBC$^{2+}$ has its main absorption peaks at 660 (15.15), 766 (13.05), 1054 (9.49), 1176 (8.50), 1290 (7.75), 1370 (7.30) and 1530 (6.54) cm$^{-1}$($\mu$m). Given the similarity in the cationic and di-cationic spectra, we have not identified an obvious diagnostic signature to the presence of multiply charged PAHs in space. While experimental issues associated with the IRMPD technique preclude a detailed comparison with interstellar spectra, we do note that the strong bands of HBC$^+$ and HBC$^{2+}$ at $\sim$ 6.5, 7.7, 8.4 and 13.1 $\mu$m coincide with prominent aromatic infrared bands (AIBs). HBC has only trio CH groups and the out-of-plane CH bending mode of both HBC cations is measured at 13.1 $\mu$m, squarely in the range predicted by theory and previously found in studies of small (substituted) PAHs. This study therefore supports the use of AIBs observed in the 11$-$14 $\mu$m range as a diagnostic tool for the edge topology of large PAHs in space.

The first gas-phase infrared spectra of two isolated astronomically relevant and large PAH cations - diindenoperylene (DIP) and dicoronylene (DC) - in the 530$-$1800 cm$^{-1}$ (18.9$-$5.6 $\mu$m) range - are presented. Vibrational band positions are determined for comparison to the aromatic infrared bands (AIBs). The spectra are obtained via infrared multiphoton dissociation (IRMPD) spectroscopy of ions stored in a quadrupole ion trap (QIT) using the intense and tunable radiation of the free electron laser for infrared experiments (FELIX). DIP$^{+}$ shows its main absorption peaks at 737 (13.57), 800 (12.50), 1001 (9.99), 1070 (9.35), 1115 (8.97), 1152 (8.68), 1278 (7.83), 1420 (7.04) and 1550 (6.45) cm$^{-1}$($\mu$m), in good agreement with DFT calculations that are uniformly scaled to take anharmonicities into account. DC$^+$ has its main absorption peaks at 853 (11.72), 876 (11.42), 1032 (9.69), 1168 (8.56), 1300 (7.69), 1427 (7.01) and 1566 (6.39) cm$^{-1}$($\mu$m), that also agree well with the scaled DFT results presented here. The DIP$^+$ and DC$^+$ spectra are compared with the prominent infrared features observed towards NGC 7023. This results both in matches and clear deviations. Moreover, in the 11.0$-$14.0 $\mu$m region, specific bands can be linked to CH out-of-plane (oop) bending modes of different CH edge structures in large PAHs. The molecular origin of these findings and their astronomical relevance are discussed.

The physical and chemical conditions in photodissociation regions (PDRs) are largely determined by the influence of far ultraviolet radiation. Far-UV photons can efficiently dissociate molecular hydrogen, a process that must be balanced at the HI/H2 interface of the PDR. Given that reactions involving hydrogen atoms in the gas phase are highly inefficient under interstellar conditions, H2 formation models mostly rely on catalytic reactions on the surface of dust grains. Additionally, molecular hydrogen formation in polycyclic aromatic hydrocarbons (PAHs) through the Eley-Rideal mechanism has been considered as well, although it has been found to have low efficiency in PDR fronts. In a previous work, we have described the possibility of efficient H2 release from medium to large sized PAHs upon photodissociation, with the exact branching between H-/H2-loss reactions being molecule dependent. Here we investigate the astrophysical relevance of this process, by using a model for the photofragmentation of PAHs under interstellar conditions. We focus on three PAHs cations (coronene, ovalene and circumcoronene), which represent three possibilities in the branching of atomic and molecular hydrogen losses. We find that, for ovalene (H2-loss dominated) the rate coefficient for H2 formation reaches values of the same order as H2 formation in dust grains. This result suggests that this hitherto disregarded mechanism can account, at least partly, for the high level of molecular hydrogen formation in dense PDRs.

PAHs constitute a major component of the interstellar medium carbon budget, locking up to 10--20% of the elemental carbon. Sequential fragmentation induced by energetic photons leads to the formation of new species, including fullerenes. However, the exact chemical routes involved in this process remain largely unexplored. In this work, we focus on the first photofragmentation steps, which involve the dehydrogenation of these molecules. For this, we consider a multidisciplinary approach, taking into account the results from experiments, DFT calculations, and modeling using dedicated Monte-Carlo simulations. By considering the simplest isomerization pathways --- i.e., hydrogen roaming along the edges of the molecule --- we are able to characterize the most likely photodissociation pathways for the molecules studied here. These comprise nine PAHs with clearly different structural properties. The formation of aliphatic-like side groups is found to be critical in the first fragmentation step and, furthermore, sets the balance of the competition between H- and H2-loss. We show that the presence of trio hydrogens, especially in combination with bay regions in small PAHs plays an important part in the experimentally established variations in the odd-to-even H-atom loss ratios. In addition, we find that, as PAH size increases, H2 formation becomes dominant, and sequential hydrogen loss only plays a marginal role. We also find disagreements between experiments and calculations for large, solo containing PAHs, which need to be accounted for. In order to match theoretical and experimental results, we have modified the energy barriers and restricted the H-hopping to tertiary atoms. The formation of H2 in large PAHs upon irradiation appears to be the dominant fragmentation channel, suggesting an efficient formation path for molecular hydrogen in PDRs.

The bipolar nebula Menzel 3 (Mz 3) was observed as part of the \textit{Herschel} Planetary Nebula Survey (\textit{HerPlaNS}), which used the PACS and SPIRE instruments aboard the \textit{Herschel Space Observatory} to study a sample of planetary nebulae (PNe). In this paper, one of the series describing \textit{HerPlaNS} results, we report the detection of H I recombination lines (HRLs) in the spectrum of Mz 3. Inspection of the spectrum reveals the presence of 12 HRLs in the 55 to 680 $\mu$m range covered by the PACS and SPIRE instruments (H11$\alpha$ to H21$\alpha$ and H14$\beta$). The presence of HRLs in this range is unusual for PNe and has not been reported in Mz 3 before. Our analysis indicates that the HRLs we observed are enhanced by laser effect occurring in the core of Mz 3. Our arguments for this are: (i) the available Mz 3 optical to submillimetre HRL $\alpha$ line intensity ratios are not well reproduced by the spontaneous emission of optically thin ionized gas, as would be typical for nebular gas in PNe; (ii) the compact core of Mz 3 is responsible for a large fraction of the Herschel HRLs emission; (iii) the line intensity ratios for Mz 3 are very similar to those in the core emission of the well known star MWC 349A, where laser effect is responsible for the enhancement of HRLs in the Herschel wavelength range; (iv) the physical characteristics relevant to cause laser effect in the core of MWC 349A are very similar to those in the core of Mz 3.

Aims. We aim to elucidate the spectral changes in the 3 micron region that result from chemical changes in the molecular periphery of polycyclic aromatic hydrocarbons (PAHs) with extra hydrogens (H-PAHs) and methyl groups (Me-PAHs). Methods. Advanced laser spectroscopic techniques combined with mass spectrometry were applied on supersonically cooled 1,2,3,4-tetrahydronaphthalene, 9,10-dihydroanthracene, 9,10-dihydrophenathrene, 1,2,3,6,7,8-hexahydropyrene, 9-methylanthracene, and 9,10-dimethylanthracene, allowing us to record mass-selective and conformationally selective absorption spectra of the aromatic, aliphatic, and alkyl CH-stretches in the 3.175-3.636 micron region with laser-limited resolution. We compared the experimental absorption spectra with standard harmonic calculations and with second-order vibrational perturbation theory anharmonic calculations that use the SPECTRO program for treating resonances. Results. We show that anharmonicity plays an important if not dominant role, affecting not only aromatic, but also aliphatic and alkyl CH-stretch vibrations. The experimental high-resolution data lead to the conclusion that the variation in Me- and H-PAHs composition might well account for the observed variations in the 3 micron emission spectra of carbon-rich and star-forming regions. Our laboratory studies also suggest that heavily hydrogenated PAHs form a significant fraction of the carriers of IR emission in regions in which an anomalously strong 3 micron plateau is observed.

Quantitative understanding of the interstellar medium requires knowledge of its physical conditions. Low frequency carbon radio recombination lines (CRRLs) trace cold interstellar gas, and can be used to determine its physical conditions (e.g., electron temperature and density). In this work we present spatially resolved observations of the low frequency ($\leq390$ MHz) CRRLs centered around C$268\alpha$, C$357\alpha$, C$494\alpha$ and C$539\alpha$ towards Cassiopeia A on scales of $\leq1.2$ pc. We compare the spatial distribution of CRRLs with other ISM tracers. This comparison reveals a spatial offset between the peak of the CRRLs and other tracers, which is very characteristic for photodissociation regions and that we take as evidence for CRRLs being preferentially detected from the surfaces of molecular clouds. Using the CRRLs we constrain the gas electron temperature and density. These constraints on the gas conditions suggest variations of less than a factor of two in pressure over $\sim1$ pc scales, and an average hydrogen density of $200$-$470$ cm$^{-3}$. From the electron temperature and density maps we also constrain the ionized carbon emission measure, column density and path length. Based on these, the hydrogen column density is larger than $10^{22}$ cm$^{-2}$, with a peak of $\sim4\times10^{22}$ cm$^{-2}$ towards the South of Cassiopeia A. Towards the southern peak the line of sight length is $\sim40$ pc over a $\sim2$ pc wide structure, which implies that the gas is a thin surface layer on a large (molecular) cloud that is only partially intersected by Cassiopeia A. These observations highlight the utility of CRRLs as tracers of low density extended HI and CO-dark gas halo's around molecular clouds.

The efficiency of radial transport of icy solid material from outer disk to the inner disk is currently unconstrained. Efficient radial transport of icy dust grains could significantly alter the composition of the gas in the inner disk. Our aim is to model the gaseous CO$_2$ abundance in the inner disk and use this to probe the efficiency of icy dust transport in a viscous disk. Features in the simulated CO$_2$ spectra are investigated for their dust flux tracing potential. We have developed a 1D viscous disk model that includes gas and grain motions as well as dust growth, sublimation and freeze-out and a parametrisation of the CO$_2$ chemistry. The thermo-chemical code DALI was used to model the mid-infrared spectrum of CO$_2$, as can be observed with JWST-MIRI. CO$_2$ ice sublimating at the iceline increases the gaseous CO$_2$ abundance to levels equal to the CO$_2$ ice abundance of $\sim 10^{-5}$, which is three orders of magnitude more than the gaseous CO$_2$ abundances of $\sim 10^{-8}$ observed by Spitzer. Grain growth and radial drift further increase the gaseous CO$_2$ abundance. A CO$_2$ destruction rate of at least $10^{-11}$ s$^{-1}$ is needed to reconcile model prediction with observations. This rate is at least two orders of magnitude higher than the fastest known chemical destruction rate. A range of potential physical mechanisms to explain the low observed CO$_2$ abundances are discussed. Transport processes in disks can have profound effects on the abundances of species in the inner disk. The discrepancy between our model and observations either suggests frequent shocks in the inner 10 AU that destroy CO$_2$, or that the abundant midplane CO$_2$ is hidden from our view by an optically thick column of low abundance CO$_2$ in to the disk surface XDR/PDR. Other molecules, such as CH$_4$ or NH$_3$, can give further handles on the rate of mass transport.

Observations towards L1630 in the Orion B molecular cloud, comprising the iconic Horsehead Nebula, allow us to study the interplay between stellar radiation and a molecular cloud under relatively benign conditions, that is, intermediate densities and an intermediate UV radiation field. Contrary to the well-studied Orion Molecular Cloud 1 (OMC1), which hosts much harsher conditions, L1630 has little star formation. We aim to relate the [CII] fine-structure line emission to the physical conditions predominant in L1630 and compare it to studies of OMC1. The [CII] $158\,\mu\mathrm{m}$ emission from an area of $12' \times 17'$ in L1630 was observed using the upGREAT instrument onboard SOFIA. Of the [CII] emission from the mapped area 95%, $13\,L_{\odot}$, originates from the molecular cloud; the adjacent HII region contributes only 5%, that is, $1\,L_{\odot}$. From comparison with other data (CO (1-0)-line emission, far-infrared (FIR) continuum studies, emission from polycyclic aromatic hydrocarbons (PAHs)), we infer a gas density of the molecular cloud of $n_{\mathrm{H}}\sim 3\cdot 10^3\,\mathrm{cm^{-3}}$, with surface layers, including the Horsehead Nebula, having a density of up to $n_{\mathrm{H}}\sim 4\cdot 10^4\,\mathrm{cm^{-3}}$. The temperature of the surface gas is $T\sim 100\,\mathrm{K}$. The average [CII] cooling efficiency within the molecular cloud is $1.3\cdot 10^{-2}$. The fraction of the mass of the molecular cloud within the studied area that is traced by [CII] is only $8\%$. Our PDR models are able to reproduce the FIR-[CII] correlations and also the CO (1-0)-[CII] correlations. Finally, we compare our results on the heating efficiency of the gas with theoretical studies of photoelectric heating by PAHs, clusters of PAHs, and very small grains, and find the heating efficiency to be lower than theoretically predicted, a continuation of the trend set by other observations.

During its 2012 decline the R Coronae Borealis star (RCB) V854 Cen was spectroscopically monitored with X-shooter on the ESO Very Large Telescope. The obscured optical and near-infrared spectrum exhibits many narrow and several broad emission features, as previously observed. The envelope is spatially resolved along the slit and allows for a detailed study of the circumstellar material. In this Letter we report on the properties of a number of unidentified emission features (UFs), including the detection of a new one at 8692 {\AA}. These UFs have been observed in the Red Rectangle, but their chemical and physical nature is still a mystery. The previously known UFs behave similarly in the Red Rectangle and V854 Cen, but are not detected in six other observed RCBs. Possibly the presence of some hydrogen is required for the formation of their carrier(s). The $\lambda 8692$ UF is present in all RCBs. Its carrier is likely of a carbonaceous molecular nature, presumably different from that of the other UFs.

We use the Spitzer SAGE survey of the Magellanic Clouds to evaluate the relationship between the 8-mic PAH emission, 24-mic hot dust emission, and HII region radiative transfer. We confirm that in the higher-metallicity Large Magellanic Cloud, PAH destruction is sensitive to optically thin conditions in the nebular Lyman continuum: objects identified as optically thin candidates based on nebular ionization structure show 6 times lower median 8-mic surface brightness (0.18 mJy arcsec^-2) than their optically thick counterparts (1.2 mJy arcsec^-2). The 24-mic surface brightness also shows a factor of 3 offset between the two classes of objects (0.13 vs 0.44 mJy arcsec^-2, respectively), which is driven by the association between the very small dust grains and higher density gas found at higher nebular optical depths. In contrast, PAH and dust formation in the low-metallicity Small Magellanic Cloud is strongly inhibited such that we find no variation in either 8-mic or 24-mic emission between our optically thick and thin samples. This is attributable to extremely low PAH and dust production together with high, corrosive UV photon fluxes in this low-metallicity environment. The dust mass surface densities and gas-to-dust ratios determined from dust maps using Herschel HERITAGE survey data support this interpretation.

We present a study of carbon radio recombination lines towards Cassiopeia A using LOFAR observations in the frequency range 10-33 MHz. Individual carbon $\alpha$ lines are detected in absorption against the continuum at frequencies as low as 16 MHz. Stacking several C$\alpha$ lines we obtain detections in the 11-16 MHz range. These are the highest signal-to-noise measurements at these frequencies. The peak optical depth of the C$\alpha$ lines changes considerably over the 11-33 MHz range with the peak optical depth decreasing from 4$\times10^{-3}$ at 33 MHz to 2$\times10^{-3}$ at 11 MHz, while the line width increases from 20 km s$^{-1}$ to 150 km s$^{-1}$. The combined change in peak optical depth and line width results in a roughly constant integrated optical depth. We interpret this as carbon atoms close to local thermodynamic equilibrium. In this work we focus on how the 11-33 MHz carbon radio recombination lines can be used to determine the gas physical conditions. We find that the ratio of the carbon radio recombination lines to that of the 158 $\mu$m [CII] fine-structure line is a good thermometer, while the ratio between low frequency carbon radio recombination lines provides a good barometer. By combining the temperature and pressure constraints with those derived from the line width we are able to constrain the gas properties (electron temperature and density) and radiation field intensity. Given the 1$\sigma$ uncertainties in our measurements these are; $T_{e}\approx68$-$98$ K, $n_{e}\approx0.02$-$0.035$ cm$^{-3}$ and $T_{r,100}\approx1500$-$1650$ K. Despite challenging RFI and ionospheric conditions, our work demonstrates that observations of carbon radio recombination lines in the 10-33 MHz range can provide insight into the gas conditions.

We present 5-20 micron spectral maps of the reflection nebula NGC2023 obtained with the Infrared Spectrograph SL and SH modes on board the Spitzer Space Telescope which reveal emission from polycyclic aromatic hydrocarbons (PAHs), C60, and H2 superposed on a dust continuum. We show that several PAH emission bands correlate with each other and exhibit distinct spatial distributions revealing a spatial sequence with distance from the illuminating star. We explore the distinct morphology of the 6.2, 7.7 and 8.6 micron PAH bands and find that at least two spatially distinct components contribute to the 7--9 micron PAH emission in NGC2023. We report that the PAH features behave independently of the underlying plateaus. We present spectra of compact oval PAHs ranging in size from C_66 to C_210, determined computationally using density functional theory, and investigate trends in the band positions and relative intensities as a function of PAH size, charge and geometry. Based on the NASA Ames PAH database, we discuss the 7--9 micron components in terms of band assignments and relative intensities. We assign the plateau emission to very small grains with possible contributions from PAH clusters and identify components in the 7--9 micron emission that likely originates in these structures. Based on the assignments and the observed spatial sequence, we discuss the photochemical evolution of the interstellar PAH family as they are more and more exposed to the radiation field of the central star in the evaporative flows associated with the PDRs in NGC2023.

We study the absorption and scattering of X-ray radiation by interstellar dust particles, which allows us to access the physical and chemical properties of dust. The interstellar dust composition is not well understood, especially on the densest sight lines of the Galactic Plane. X-rays provide a powerful tool in this study. We present newly acquired laboratory measurements of silicate compounds taken at the Soleil synchrotron facility in Paris using the Lucia beamline. The dust absorption profiles resulting from this campaign were used in this pilot study to model the absorption by interstellar dust along the line of sight of the low-mass X-ray binary (LMXB) GX 5-1. The measured laboratory cross-sections were adapted for astrophysical data analysis and the resulting extinction profiles of the Si K-edge were implemented in the SPEX spectral fitting program. We derive the properties of the interstellar dust along the line of sight by fitting the Si K-edge seen in absorption in the spectrum of GX 5-1. We measured the hydrogen column density towards GX 5-1 to be $3.40\pm0.1\times10^{22} \rm cm^{-2}$. The best fit of the silicon edge in the spectrum of GX 5-1 is obtained by a mixture of olivine and pyroxene. In this study, our modeling is limited to Si absorption by silicates with different Mg:Fe ratios. We obtained an abundance of silicon in dust of $4.0\pm0.3\times10^{-5}$ per H atom and a lower limit for total abundance, considering both gas and dust, of $>4.4\times10^{-5}$ per H atom, which leads to a gas to dust ratio of >0.22. Furthermore, an enhanced scattering feature in the Si K-edge may suggest the presence of large particles along the line of sight.

The 2175 \AA\ UV extinction feature was discovered in the mid-1960s, yet its physical origin remains poorly understood. One suggestion is absorption by Polycyclic Aromatic Hydrocarbons (PAH) molecules, which is supported by theoretical molecular structure computations and by laboratory experiments. PAHs are positively detected by their 3.3, 6.2, 7.7 8.6, 11.3 & 12.7 $\mu$m IR emission bands, which are specified by their modes of vibration. A definitive empirical link between the 2175 \AA\ UV extinction, and the PAH IR emission bands, however, is still missing. We present a new sample of hot stars that have both 2175 \AA\ absorption and PAH IR emission. We find significant shifts of the central wavelength of the UV absorption feature, up to 2350 \AA, but predominantly in stars that also have IR PAH emission. These UV shifts depend on stellar temperature in a fashion that is similar to the shifts of the 6.2 and 7.7$\mu$m PAH IR bands, namely the features are increasingly more red-shifted as the stellar temperature decreases, but only below $\sim 15$ kK. Above 15 kK both UV and IR features retain their nominal values. Moreover, we find a suggestive correlation between the UV and IR shifts. We hypothesize that these similar dependences of both the UV and IR features on stellar temperature hint to a common origin of the two in PAH molecules and may establish the missing link between the UV and IR observations. We further suggest that the shifts depend on molecular size, and that the critical temperature of $\sim 15$ kK above which no shifts are observed is related to the onset of UV driven hot-star winds and their associated shocks.

We use the Low Frequency Array to perform a systematic high spectral resolution investigation of the low-frequency 33-78 MHz spectrum along the line of sight to Cassiopeia A. We complement this with a 304-386 MHz Westerbork Synthesis Radio telescope observation. In this first paper we focus on the carbon radio recombination lines. We detect Cn$\alpha$ lines at -47 and -38 km s$^{-1}$ in absorption for quantum numbers n=438-584 and in emission for n=257-278 with high signal to noise. These lines are associated with cold clouds in the Perseus spiral arm component. Hn$\alpha$ lines are detected in emission for n=257-278. In addition, we also detect Cn$\alpha$ lines at 0 km s$^{-1}$ associated with the Orion arm. We analyze the optical depth of these transitions and their line width. Our models show that the carbon line components in the Perseus arm are best fit with an electron temperature 85 K and an electron density 0.04 cm$^{-3}$ and can be constrained to within 15\%. The electron pressure is constrained to within 20\%. We argue that much of these carbon radio recombination lines arise in the CO-dark surface layers of molecular clouds where most of the carbon is ionized but hydrogen has made the transition from atomic to molecular. The hydrogen lines are clearly associated with the carbon line emitting clouds, but the low-frequency upperlimits indicate that they likely do not trace the same gas. Combining the hydrogen and carbon results we arrive at a firm lower limit to the cosmic ray ionization rate of 2.5$\times$10$^{-18}$ s$^{-1}$, but the actual value is likely much larger.

In this work we report on high-resolution IR absorption studies that provide a detailed view on how the peripheral structure of irregular polycyclic aromatic hydrocarbons (PAHs) affects the shape and position of their 3 {\mu}m absorption band. To this purpose we present mass-selected, high-resolution absorption spectra of cold and isolated phenanthrene, pyrene, benz[a]antracene, chrysene, triphenylene, and perylene molecules in the 2950-3150 cm-1 range. The experimental spectra are compared with standard harmonic calculations, and anharmonic calculations using a modified version of the SPECTRO program that incorporates a Fermi resonance treatment utilizing intensity redistribution. We show that the 3 {\mu}m region is dominated by the effects of anharmonicity, resulting in many more bands than would have been expected in a purely harmonic approximation. Importantly, we find that anharmonic spectra as calculated by SPECTRO are in good agreement with the experimental spectra. Together with previously reported high-resolution spectra of linear acenes, the present spectra provide us with an extensive dataset of spectra of PAHs with a varying number of aromatic rings, with geometries that range from open to highly-condensed structures, and featuring CH groups in all possible edge configurations. We discuss the astrophysical implications of the comparison of these spectra on the interpretation of the appearance of the aromatic infrared 3 {\mu}m band, and on features such as the two-component emission character of this band and the 3 {\mu}m emission plateau.

In the second paper of the series, we have modeled low frequency carbon radio recombination lines (CRRL) from the interstellar medium. Anticipating the LOw Frequency ARray (LOFAR) survey of Galactic CRRLs, we focus our study on the physical conditions of the diffuse cold neutral medium (CNM). We have used the improved departure coefficients computed in the first paper of the series to calculate line-to-continuum ratios. The results show that the line width and integrated optical depths of CRRL are sensitive probes of the electron density, gas temperature, and the emission measure of the cloud. Furthermore, the ratio of CRRL to the [CII] at 158 $\mu$m line is a strong function of the temperature and density of diffuse clouds. Guided by our calculations, we analyze CRRL observations and illustrate their use with data from the literature.

In the first paper of this series, we study the level population problem of recombining carbon ions. We focus our study on high quantum numbers anticipating observations of Carbon Radio Recombination Lines to be carried out by the LOw Frequency ARray (LOFAR). We solve the level population equation including angular momentum levels with updated collision rates up to high principal quantum numbers. We derive departure coefficients by solving the level population equation in the hydrogenic approximation and including low temperature dielectronic recombination effects. Our results in the hydrogenic approximation agree well with those of previous works. When comparing our results including dielectronic recombination we find differences which we ascribe to updates in the atomic physics (e.g., collision rates) and to the approximate solution method of the statistical equilibrium equations adopted in previous studies. A comparison with observations is discussed in an accompanying article, as radiative transfer effects need to be considered.

Massive stars are key players in the evolution of galaxies, yet their formation pathway remains unclear. In this work, we use data from several galaxy-wide surveys to build an unbiased dataset of ~700 massive young stellar objects (MYSOs), ~200 giant molecular clouds (GMCs), and ~100 young (<10 Myr) optical stellar clusters (SCs) in the Large Magellanic Cloud. We employ this data to quantitatively study the location and clustering of massive star formation and its relation to the internal structure of GMCs. We reveal that massive stars do not typically form at the highest column densities nor centers of their parent GMCs at the ~6 pc resolution of our observations. Massive star formation clusters over multiple generations and on size scales much smaller than the size of the parent GMC. We find that massive star formation is significantly boosted in clouds near SCs. Yet, whether a cloud is associated with a SC does not depend on either the cloud's mass or global surface density. These results reveal a connection between different generations of massive stars on timescales up to 10 Myr. We compare our work with Galactic studies and discuss our findings in terms of GMC collapse, triggered star formation, and a potential dichotomy between low- and high-mass star formation.

We present mid-infrared photometry of the Orion Bar obtained with FORCAST aboard SOFIA at 6.4, 6.6, 7.7, 19.7, 31.5 and 37.1 \um. By complementing this observations with archival FORCAST and \emph{Herschel}/PACS images we are able to construct a complete infrared spectral energy distribution of the Huygens region in the Orion nebula By comparing the infrared images with gas tracers, we find that PACS maps trace the molecular cloud, while the FORCAST data trace the photodissociation region (PDR) and HII region. Analysis of the energetics of the region reveal that the PDR extends for 0.28~pc along the line-of-sight and that the Bar is inclined at an angle of $4\degr$. The infrared and submillimeter images reveal that the Orion Bar represents a swept up shell with a thickness of 0.1~pc. The mass of the shell implies a shock velocity of $\simeq 3$ km/s and an age of $\simeq 10^5$ yr for the HII region. Our analysis shows that the UV and infrared dust opacities in the HII region and the PDR are a factor 5 to 10 lower than in the diffuse interstellar medium. In the ionized gas, Ly$\alpha$ photons are a major source of dust heating at distances larger than $\simeq0.06$~pc from \toc. Dust temperatures can be explained if the size of the grains is between 0.1 to 1~\um. We derive the photo-electric heating efficiency of the atomic gas in the Orion Bar. The results are in good qualitative agreement with models and The quantitative differences indicate a decreased PAH abundance in this region.

The Orion star-forming region is the nearest active high-mass star-forming region and has created a large superbubble, the Orion-Eridanus superbubble. Recent work by Ochsendorf et al. has extended the accepted boundary of the superbubble. We fit Kompaneets models of superbubbles expanding in exponential atmospheres to the new, larger shape of the Orion-Eridanus superbubble. We find that this larger morphology of the superbubble is consistent with the evolution of the superbubble being primarily controlled by expansion into the exponential Galactic disk ISM if the superbubble is oriented with the Eridanus side farther from the Sun than the Orion side. Unlike previous Kompaneets model fits that required abnormally small scale heights for the Galactic disk (<40 pc), we find morphologically consistent models with scale heights of 80 pc, similar to that expected for the Galactic disk.

IRAS 19312+1950 is a peculiar object that has eluded firm characterization since its discovery, with combined maser properties similar to an evolved star and a young stellar object (YSO). To help determine its true nature, we obtained infrared spectra of IRAS 19312+1950 in the range 5-550 $\mu$m using the Herschel and Spitzer space observatories. The Herschel PACS maps exhibit a compact, slightly asymmetric continuum source at 170 $\mu$m, indicative of a large, dusty circumstellar envelope. The far-IR CO emission line spectrum reveals two gas temperature components: $\approx0.22M_{\odot}$ of material at $280\pm18$ K, and $\approx1.6M_{\odot}$ of material at $157\pm3$ K. The OI 63 $\mu$m line is detected on-source but no significant emission from atomic ions was found. The HIFI observations display shocked, high-velocity gas with outflow speeds up to 90 km s$^{-1}$ along the line of sight. From Spitzer spectroscopy, we identify ice absorption bands due to H$_2$O at 5.8 $\mu$m and CO$_2$ at 15 $\mu$m. The spectral energy distribution is consistent with a massive, luminous ($\sim2\times10^4L_{\odot}$) central source surrounded by a dense, warm circumstellar disk and envelope of total mass $\sim500$-$700M_{\odot}$, with large bipolar outflow cavities. The combination of distinctive far-IR spectral features suggest that IRAS 19312+1950 should be classified as an accreting high-mass YSO rather than an evolved star. In light of this reclassification, IRAS 19312+1950 becomes only the 5th high-mass protostar known to exhibit SiO maser activity, and demonstrates that 18 cm OH maser line ratios may not be reliable observational discriminators between evolved stars and YSOs.

NGC 7023 is a well-studied reflection nebula, which shows strong emission from polycyclic aromatic hydrocarbon (PAH) molecules in the form of aromatic infrared bands (AIBs). The spectral variations of the AIBs in this region are connected to the chemical evolution of the PAH molecules which, in turn, depends on the local physical conditions. We use the capabilities of SOFIA to observe a 3.2' x 3.4' region of NGC 7023 at wavelengths that we observe with high spatial resolution (2.7") at 3.3 and 11.2 um. We compare the SOFIA images with existing images of the PAH emission at 8.0 um (Spitzer), emission from evaporating very small grains (eVSG) extracted from Spitzer-IRS spectral cubes, the ERE (HST and CFHT), and H_2 (2.12 um). We create maps of the 11.2/3.3 um ratio to probe the morphology of the PAH size distribution and the 8.0/11.2 um ratio to probe the PAH ionization. We make use of an emission model and of vibrational spectra from the NASA Ames PAHdb to translate the 11.2/3.3 um ratio to PAH sizes. The 11.2/3.3 um map shows the smallest PAH concentrate on the PDR surface (H_2 and extended red emission) in the NW and South PDR. We estimated that PAHs in the NW PDR bear, on average, a number of carbon atoms (N_c) of ~70 in the PDR cavity and ~50 at the PDR surface. In the entire nebula, the results reveal a factor of 2 variation in the size of the PAH. We relate these size variations to several models for the evolution of the PAH families when they traverse from the molecular cloud to the PDR. The PAH size map enables us to follow the photochemical evolution of PAHs in NGC 7023. Small PAHs result from the photo-evaporation of VSGs as they reach the PDR surface. Inside the PDR cavity, the PAH abundance drops as the smallest PAH are broken down. The average PAH size increases in the cavity where only the largest species survive or are converted into C_60 by photochemical processing.

We present a sample of resolved galactic HII regions and photodissociation regions (PDRs) observed with the Spitzer infrared spectrograph (IRS) in spectral mapping mode between the wavelengths of 5--15 $\mu$m. For each object we have spectral maps at a spatial resolution of $\sim$4" in which we have measured all of the mid-infrared emission and absorption features. These include the PAH emission bands, primarily at 6.2, 7.7, 8.6, 11.2 and 12.7 $\mu$m, as well as the spectral emission lines of neon and sulfur and the absorption band caused by silicate dust at around 9.8 $\mu$m. In this work we describe the data in detail, including the data reduction and measurement strategies, and subsequently present the PAH emission band intensity correlations for each of the objects and the sample as a whole. We find that there are distinct differences between the sources in the sample, with two main groups, the first comprising the HII regions and the second the reflection nebulae (RNe). Three sources, the reflection nebula NGC~7023, the Horsehead nebula PDR (an interface between the HII region IC~434 and the Orion B molecular cloud) and M 17, resist this categorization, with the Horsehead PDR points mimicking the RNe and the NGC~7023 fluxes displaying unique bifurcated appearance in our correlation plots. These discrepancies seem to be due to the very low radiation field experienced by the Horsehead PDR and the very clean separation between the PDR environment and a diffuse environment in the NGC~7023 observations.

The amount of deuterium locked up in polycyclic aromatic hydrocarbons (PAHs) has to date been an uncertain value. We present a near-infrared (NIR) spectroscopic survey of HII regions in the Milky Way, Large Magellanic Cloud (LMC), and Small Magellanic Cloud (SMC) obtained with AKARI, which aims to search for features indicative of deuterated PAHs (PAD or Dn-PAH) to better constrain the D/H ratio of PAHs. Fifty-three HII regions were observed in the NIR (2.5-5 {\mu}m), using the Infrared Camera (IRC) on board the AKARI satellite. Through comparison of the observed spectra with a theoretical model of deuterated PAH vibrational modes, the aromatic and (a)symmetric aliphatic C-D stretch modes were identified. We see emission features between 4.4-4.8 {\mu}m, which could be unambiguously attributed to deuterated PAHs in only six of the observed sources, all of which are located in the Milky Way. In all cases, the aromatic C-D stretching feature is weaker than the aliphatic C-D stretching feature, and, in the case of M17b, this feature is not observed at all. Based on the weak or absent PAD features in most of the observed spectra, it is suggested that the mechanism for PAH deuteration in the ISM is uncommon.

The journey from dust particle to planetesimal involves physical processes acting on scales ranging from micrometers (the sticking and restructuring of aggregates) to hundreds of astronomical units (the size of the turbulent protoplanetary nebula). Considering these processes simultaneously is essential when studying planetesimal formation. We develop a novel, global, semi-analytical model for the evolution of the mass-dominating dust particles in a turbulent protoplanetary disk that takes into account the evolution of the dust surface density while preserving the essential characteristics of the porous coagulation process. This panoptic model is used to study the growth from submicron to planetesimal sizes in disks around Sun-like stars. For highly porous ices, unaffected by collisional fragmentation and erosion, rapid growth to planetesimal sizes is possible in a zone stretching out to ${\sim}10\mathrm{~AU}$ for massive disks. When porous coagulation is limited by erosive collisions, the formation of planetesimals through direct coagulation is not possible, but the creation of a large population of aggregates with Stokes numbers close to unity might trigger the streaming instability (SI). However, we find that reaching conditions necessary for SI is difficult and limited to dust-rich disks, (very) cold disks, or disks with weak turbulence. Behind the snow-line, porosity-driven aggregation of icy grains results in rapid (${\sim}10^{4}\mathrm{~yr}$) formation of planetesimals. If erosive collisions prevent this, SI might be triggered for specific disk conditions. The numerical approach introduced in this work is ideally suited for studying planetesimal formation and pebble delivery simultaneously and will help build a coherent picture of the start of the planet formation process.

We report on an experimental and theoretical investigation of the importance of anharmonicity in the 3 micron CH stretching region of Polycyclic Aromatic Hydrocarbon (PAH) molecules. We present mass-resolved, high-resolution spectra of the gas-phase cold (~4K) linear PAH molecules naphthalene, anthracene, and tetracene. The measured IR spectra show a surprisingly high number of strong vibrational bands. For naphthalene, the observed bands are well separated and limited by the rotational contour, revealing the band symmetries. Comparisons are made to the harmonic and anharmonic approaches of the widely used Gaussian software. We also present calculated spectra of these acenes using the computational program SPECTRO, providing anharmonic predictions enhanced with a Fermi-resonance treatment that utilises intensity redistribution. We demonstrate that the anharmonicity of the investigated acenes is strong, dominated by Fermi resonances between the fundamental and double combination modes, with triple combination bands as possible candidates to resolve remaining discrepancies. The anharmonic spectra as calculated with SPECTRO lead to predictions of the main modes that fall within 0.5% of the experimental frequencies. The implications for the Aromatic Infrared Bands, specifically the 3 micron band are discussed.

The formation of dust gaps in protoplanetary disks is one of the most important signposts of disk evolution and possibly the formation of planets. We aim to characterize the 'flaring' disk structure around the Herbig Ae/Be stars HD 100453 and HD 34282. Their spectral energy distributions (SEDs) show an emission excess between 15-40{\mu}m, but very weak (HD 100453) and no (HD 34282) signs of the 10 and 20 {\mu}m amorphous silicate features. We investigate whether this implies the presence of large dust gaps. In this work, spatially resolved mid-infrared Q-band images taken with Gemini North/MICHELLE are investigated. We perform radiative transfer modeling and examine the radial distribution of dust. We simultaneously fit the Q-band images and SEDs of HD 100453 and HD 34282. Our solutions require that the inner-halos and outer-disks are likely separated by large dust gaps that are depleted wih respect to the outer disk by a factor of 1000 or more. The inner edges of the outer disks of HD 100453 and HD 34282 have temperatures of about $160 \pm 10$ K and $60 \pm 5$ K respectively. Because of the high surface brightnesses of these walls, they dominate the emission in the Q-band. Their radii are constrained at 20+2 AU and 92+31 AU, respectively. We conclude that, HD 100453 and HD 34282 likely have disk dust gaps and the upper limit on the dust mass in each gap is estimated to be about $10^{-7}$M$_{\odot}$. We find that the locations and sizes of disk dust gaps are connected to the SED, as traced by the mid-infrared flux ratio F30/F13.5. We propose a new classification scheme for the Meeus groups (Meeus et al. 2001) based on the F30/F13.5 ratio. The absence of amorphous silicate features in the observed SEDs is caused by the depletion of small (smaller than 1 {\mu}m) silicate dust at temperatures above 160 K, which could be related to the presence of a dust gap in that region of the disk.

Planets are thought to form in the gas and dust disks around young stars. In particular, it has been proposed that giant planets can form through the gravitational instability of massive extended disks around intermediate-mass stars. However, we still lack direct observations to constrain this mechanism. We have spatially resolved the 8.6 and 11.2 $\mu$m emission of a massive protoplanetary disk seen edge on around an A star, Gomez's Hamburger (GoHam), using VISIR at the Very Large Telescope. A compact region situated at a projected distance of $350\pm50$ AU south of the central star is found to have a reduced emission.This asymmetry is fully consistent with the presence of a cold density structure, or clump, identified in earlier CO observations, and we derive physical characteristics consistent with those observations: a mass of 0.8-11.4 Jupiter masses (for a dust-to-gas mass ratio of 0.01), a radius of about 10$^2$ astronomical units, and a local density of about $10^{7}$ cm$^{-3}$. Based on this evidence, we argue that this clump, which we call GoHam b, is a promising candidate for a young protoplanet formed by gravitational instability that might be representative of the precursors of massive planets observed around A stars, such as HR 8799 or Beta pictoris. More detailed studies at high angular resolution are needed to better constrain the physical properties of this object to confirm this proposal.

The Orion-Eridanus superbubble is the prototypical superbubble due to its proximity and evolutionary state. Here, we provide a synthesis of recent observational data from WISE and Planck with archival data, allowing to draw a new and more complete picture on the history and evolution of the Orion-Eridanus region. We discuss the general morphological structures and observational characteristics of the superbubble, and derive quantitative properties of the gas- and dust inside Barnard's Loop. We reveal that Barnard's Loop is a complete bubble structure which, together with the lambda Ori region and other smaller-scale bubbles, expands within the Orion-Eridanus superbubble. We argue that the Orion-Eridanus superbubble is larger and more complex than previously thought, and that it can be viewed as a series of nested shells, superimposed along the line of sight. During the lifetime of the superbubble, HII region champagne flows and thermal evaporation of embedded clouds continuously mass-load the superbubble interior, while winds or supernovae from the Orion OB association rejuvenate the superbubble by sweeping up the material from the interior cavities in an episodic fashion, possibly triggering the formation of new stars that form shells of their own. The steady supply of material into the superbubble cavity implies that dust processing from interior supernova remnants is more efficient than previously thought. The cycle of mass-loading, interior cleansing, and star formation repeats until the molecular reservoir is depleted or the clouds have been disrupted. While the nested shells come and go, the superbubble remains for tens of millions of years.

Interstellar polycyclic aromatic hydrocarbons (PAHs) are expected to be strongly processed by vacuum ultraviolet photons. Here, we report experimental studies on the ionization and fragmentation of coronene (C24H12), ovalene (C32H14) and hexa-peri-hexabenzocoronene (HBC; C42H18) cations by exposure to synchrotron radiation in the range of 8--40 eV. The results show that for small PAH cations such as coronene, fragmentation (H-loss) is more important than ionization. However, as the size increases, ionization becomes more and more important and for the HBC cation, ionization dominates. These results are discussed and it is concluded that, for large PAHs, fragmentation only becomes important when the photon energy has reached the highest ionization potential accessible. This implies that PAHs are even more photo-stable than previously thought. The implications of this experimental study for the photo-chemical evolution of PAHs in the interstellar medium are briefly discussed.

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.

Recent studies of dust in the interstellar medium have challenged the capabilities and validity of current dust models, indicating that the properties of dust evolve as it transits between different phases of the interstellar medium. We conduct a multi-wavelength study of the dust emission from the ionized gas of the IC 434 emission nebula, and combine this with modeling, from large scales that provide insight into the history of the IC 434/L1630 region, to small scales that allow us to infer quantitative properties of the dust content inside the H II region. The dust enters the H II region through momentum transfer with a champagne flow of ionized gas, set up by a chance encounter between the L1630 molecular cloud and the star cluster of $\sigma$ Ori. We observe two clearly separated dust populations inside the ionized gas, that show different observational properties, as well as contrasting optical properties. Population A is colder ($\sim$ 25 K) than predicted by widely-used dust models, its temperature is insensitive to an increase of the impinging radiation field, is momentum-coupled to the gas, and efficiently absorbs radiation pressure to form a dust wave at 1.0 pc ahead of $\sigma$ Ori AB. Population B is characterized by a constant [20/30] flux ratio throughout the HII region, heats up to $\sim$ 75 K close to the star, and is less efficient in absorbing radiation pressure, forming a dust wave at 0.1 pc from the star. We conclude that the dust inside IC 434 is bimodal. The characteristics of population A are remarkable and can not be explained by current dust models. Population B are grains that match the classical description of spherical, compact dust. Our results confirm recent work that stress the importance of variations in the dust properties between different regions of the interstellar medium.

The Square Kilometre Array (SKA) will transform our understanding of the role of the cold, atomic gas in galaxy evolution. The interstellar medium (ISM) is the repository of stellar ejecta and the birthsite of new stars and, hence, a key factor in the evolution of galaxies over cosmic time. Cold, diffuse, atomic clouds are a key component of the ISM, but so far this phase has been difficult to study, because its main tracer, the HI 21 cm line, does not constrain the basic physical information of the gas (e.g., temperature, density) well. The SKA opens up the opportunity to study this component of the ISM through a complementary tracer in the form of low-frequency (<350 MHz) carbon radio recombination lines (CRRL). These CRRLs provide a sensitive probe of the physical conditions in cold, diffuse clouds. The superb sensitivity, large field of view, frequency resolution and coverage of the SKA allows for efficient surveys of the sky, that will revolutionize the field of low-frequency recombination line studies. By observing these lines with the SKA we will be able determine the thermal balance, chemical enrichment, and ionization rate of the cold, atomic medium from degree-scales down to scales corresponding to individual clouds and filaments in our Galaxy, the Magellanic Clouds and beyond. Furthermore, being sensitive only to the cold, atomic gas, observations of low-frequency CRRLs with the SKA will aid in disentangling the warm and cold constituents of the HI 21 cm emission.

The coagulation of microscopic dust into planetesimals is the first step towards planet formation. The size and shape of the growing aggregates determine the efficiency of this early growth. It has been proposed that fluffy ice aggregates can grow very efficiently, suffering less from the bouncing and radial drift barriers. While the collision velocity between icy aggregates of similar size is thought to stay below the fragmentation threshold, they may nonetheless lose mass from collisions with much smaller projectiles. We investigate the effect of these erosive collisions on the ability of porous ice aggregates to cross the radial drift barrier. We develop a Monte Carlo code that calculates the evolution of the growing aggregates, while resolving the entire mass distribution at all times. The aggregate's porosity is treated independently of its mass, and is determined by collisions, gas compaction, and eventually self-gravity compaction. For erosion threshold velocities of 20-40 m/s, high-velocity collisions with small projectiles prevent the largest aggregates from growing when they start to drift. In these cases, our local simulations result in a steady-state distribution, with the majority of the dust mass in particles with Stokes numbers close to unity. Only for the highest erosion threshold considered (60 m/s), do porous aggregates manage to cross the radial drift barrier in the inner 10 AU of MMSN-like disks. Erosive collisions are more effective in limiting the growth than fragmentary collisions between similar-size particles. Conceivably, erosion limits the growth before the radial drift barrier, although the robustness of this statement depends on (uncertain) material properties of icy aggregates. If erosion inhibits planetesimal formation through direct sticking, the sea of ${\sim}10^9$ g, highly porous particles appears well-suited for triggering streaming instability.

Interstellar molecules are thought to build up in the shielded environment of molecular clouds or in the envelope of evolved stars. This follows many sequential reaction steps of atoms and simple molecules in the gas phase and/or on (icy) grain surfaces. However, these chemical routes are highly inefficient for larger species in the tenuous environment of space as many steps are involved and, indeed, models fail to explain the observed high abundances. This is definitely the case for the C$_{60}$ fullerene, recently identified as one of the most complex molecules in the interstellar medium. Observations have shown that, in some PDRs, its abundance increases close to strong UV-sources. In this letter we report laboratory findings in which C$_{60}$ formation can be explained by characterizing the photochemical evolution of large PAHs. Sequential H losses lead to fully dehydrogenated PAHs and subsequent losses of C$_{2}$ units convert graphene into cages. Our results present for the first time experimental evidence that PAHs in excess of 60 C-atoms efficiently photo-isomerize to Buckminsterfullerene, C$_{60}$. These laboratory studies also attest to the importance of top-down synthesis routes for chemical complexity in space.

The freezeout of gas-phase species onto cold dust grains can drastically alter the chemistry and the heating-cooling balance of protostellar material. In contrast to well-known species such as carbon monoxide (CO), the freezeout of various carriers of elements with abundances $<10^{-5}$ has not yet been well studied. Our aim here is to study the depletion of chlorine in the protostellar core, OMC-2 FIR 4. We observed transitions of HCl and H2Cl+ towards OMC-2 FIR 4 using the Herschel Space Observatory and Caltech Submillimeter Observatory facilities. Our analysis makes use of state of the art chlorine gas-grain chemical models and newly calculated HCl-H$_{2}$ hyperfine collisional excitation rate coefficients. A narrow emission component in the HCl lines traces the extended envelope, and a broad one traces a more compact central region. The gas-phase HCl abundance in FIR 4 is 9e-11, a factor of only 0.001 that of volatile elemental chlorine. The H2Cl+ lines are detected in absorption and trace a tenuous foreground cloud, where we find no depletion of volatile chlorine. Gas-phase HCl is the tip of the chlorine iceberg in protostellar cores. Using a gas-grain chemical model, we show that the hydrogenation of atomic chlorine on grain surfaces in the dark cloud stage sequesters at least 90% of the volatile chlorine into HCl ice, where it remains in the protostellar stage. About 10% of chlorine is in gaseous atomic form. Gas-phase HCl is a minor, but diagnostically key reservoir, with an abundance of <1e-10 in most of the protostellar core. We find the 35Cl/37Cl ratio in OMC-2 FIR 4 to be 3.2\pm0.1, consistent with the solar system value.

Carbon radio recombination lines (RRLs) at low frequencies (<=500 MHz) trace the cold, diffuse phase of the interstellar medium, which is otherwise difficult to observe. We present the detection of carbon RRLs in absorption in M82 with LOFAR in the frequency range of 48-64 MHz. This is the first extragalactic detection of RRLs from a species other than hydrogen, and below 1 GHz. Since the carbon RRLs are not detected individually, we cross-correlated the observed spectrum with a template spectrum of carbon RRLs to determine a radial velocity of 219 +- 9 km/s . Using this radial velocity, we stack 22 carbon-{\alpha} transitions from quantum levels n = 468-508 to achieve an 8.5 sigma detection. The absorption line profile exhibits a narrow feature with peak optical depth of 0.003 and FWHM of 31 km/s. Closer inspection suggests that the narrow feature is superimposed on a broad, shallow component. The total line profile appears to be correlated with the 21 cm H I line profile reconstructed from H I absorption in the direction of supernova remnants in the nucleus. The narrow width and centroid velocity of the feature suggests that it is associated with the nuclear starburst region. It is therefore likely that the carbon RRLs are associated with cold atomic gas in the direction of the nucleus of M82.

Recent studies have confirmed the presence of buckminsterfullerene (C$_{60}$) in different interstellar and circumstellar environments. However, several aspects regarding C$_{60}$ in space are not well understood yet, such as the formation and excitation processes, and the connection between C$_{60}$ and other carbonaceous compounds in the interstellar medium, in particular polycyclic aromatic hydrocarbons (PAHs). In this paper we study several photodissociation regions (PDRs) where C$_{60}$ and PAHs are detected and the local physical conditions are reasonably well constrained, to provide observational insights into these questions. C$_{60}$ is found to emit in PDRs where the dust is cool ($T_d = 20-40$ K) and even in PDRs with cool stars. These results exclude the possibility for C$_{60}$ to be locked in grains at thermal equilibrium in these environments. We observe that PAH and C$_{60}$ emission are spatially uncorrelated and that C$_{60}$ is present in PDRs where the physical conditions (in terms of radiation field and hydrogen density) allow for full dehydrogenation of PAHs, with the exception of Ced 201. We also find trends indicative of an increase in C$_{60}$ abundance within individual PDRs, but these trends are not universal. These results support models where the dehydrogenation of carbonaceous species is the first step towards C$_{60}$ formation. However, this is not the only parameter involved and C$_{60}$ formation is likely affected by shocks and PDR age.

Context. The level of porosity of interstellar ices - largely comprised of amorphous solid water (ASW) - contains clues on the trapping capacity of other volatile species and determines the surface accessibility that is needed for solid state reactions to take place. Aims. Our goal is to simulate the growth of amorphous water ice at low temperature (10 K) and to characterize the evolution of the porosity (and the specific surface area) as a function of temperature (from 10 to 120 K). Methods. Kinetic Monte Carlo simulations are used to mimic the formation and the thermal evolution of pores in amorphous water ice. We follow the accretion of gas-phase water molecules as well as their migration on surfaces with different grid sizes, both at the top growing layer and within the bulk. Results. We show that the porosity characteristics change substantially in water ice as the temperature increases. The total surface of the pores decreases strongly while the total volume decreases only slightly for higher temperatures. This will decrease the overall reaction efficiency, but in parallel, small pores connect and merge, allowing trapped molecules to meet and react within the pores network, providing a pathway to increase the reaction efficiency. We introduce pore coalescence as a new solid state process that may boost the solid state formation of new molecules in space and has not been considered so far.

In this Letter we report the results of density functional theory calculations on medium-sized neutral Polycyclic Aromatic Hydrocarbon (PAH) molecules with armchair edges. These PAH molecules possess strong C-H stretching and bending modes around 3 {\mu}m and in the fingerprint region (10-15 {\mu}m), and also strong ring deformation modes around 12.7 {\mu}m. Perusal of the entries in the NASA Ames PAHs Database shows that ring deformation modes of PAHs are common - although generally weak. We then propose that armchair PAHs with NC >65 are responsible for the 12.7 {\mu}m Aromatic Infrared Band in HII regions and discuss astrophysical implications in the context of the PAH life-cycle.

We present an analysis of mid-and far-infrared (IR) spectrum and spectral energy distribution (SED) of the LBV R71 in the LMC.This work aims to understand the overall contribution of high-mass LBVs to the total dust-mass budget of the interstellar medium (ISM) of the LMC and compare this with the contribution from low-mass asymptotic giant branch (AGB) stars. As a case study, we analyze the SED of R71. We compiled all the available photometric and spectroscopic observational fluxes from various telescopes for a wide wavelength range (0.36 -- 250\,$\mu$m). We determined the dust composition from the spectroscopic data, and derived the ejected dust mass, dust mass-loss rate, and other dust shell properties by modeling the SED of R71. We noted nine spectral features in the dust shell of R71 by analyzing Spitzer spectroscopic data. Among these, we identified three new crystalline silicate features. We computed our model spectrum by using 3D radiative transfer code MCMax. Our model calculation shows that dust is dominated by amorphous silicates, with some crystalline silicates, metallic iron, and a very tiny amount of polycyclic aromatic hydrocarbon (PAH) molecules. The presence of both silicates and PAHs indicates that the dust has a mixed chemistry. We derived a dust mass of 0.01 M$_\odot$, from which we arrive at a total ejected mass of $\approx$ 5 M$_\odot$. This implies a time-averaged dust mass-loss rate of 2.5$\times$10$^{-6}$ M$_\odot$\,yr$^{-1}$ with an explosion about 4000 years ago. We assume that the other five confirmed dusty LBVs in the LMC loose mass at a similar rate, and estimate the total contribution to the mass budget of the LMC to be $\approx$ 10$^{-5}$ M$_\odot$\,yr$^{-1}$, which is comparable to the contribution by all the AGB stars in the LMC. Based on our analysis on R71, we speculate that LBVs as a class may be an important dust source in the ISM of the LMC.

We model multi-wavelength broadband UBVIJHKs and Spitzer IRAC and MIPS photometry and IRS spectra from the SAGE and SAGE-Spec observing programs of two oxygen-rich asymptotic giant branch (O-rich AGB) stars in the Large Magellanic Cloud (LMC) using radiative transfer models of dust shells around stars. We chose a star from each of the bright and faint O-rich AGB populations found by earlier studies of the SAGE sample in order to derive a baseline set of dust properties to be used in the construction of an extensive grid of radiative transfer models of the O-rich AGB stars found in the SAGE surveys. From the bright O-rich AGB population we chose HV 5715, and from the faint O-rich AGB population we chose SSTISAGE1C J052206.92-715017.6 (SSTSAGE052206). We found the complex indices of refraction of oxygen-deficient silicates from Ossenkopf et al (1992) and a "KMH"-like grain size distribution with gamma of -3.5, a_min of 0.01 microns, and a_0 of 0.1 microns to be reasonable dust properties for these models. There is a slight indication that the dust around the faint O-rich AGB may be more silica-rich than that around the bright O-rich AGB. Simple models of gas emission suggest a relatively extended gas envelope for the faint O-rich AGB star modeled, consistent with the relatively large dust shell inner radius for the same model. The models compute the dust mass-loss rates for the two stars to be 2.0x10^-9 M_Sun/yr and 2.3x10^-9 M_Sun/yr, respectively. When a dust-to-gas mass ratio of 0.002 is assumed for SSTSAGE052206 and HV 5715, the dust mass-loss rates imply total mass-loss rates of 1.0x10^-6 M_Sun/yr and 1.2x10^-6 M_Sun/yr, respectively. These properties of the dust shells and stars, as inferred from our models of the two stars, are found to be consistent with properties observed or assumed by detailed studies of other O-rich AGB stars in the LMC and elsewhere.

Comparison of observations of radio recombination lines in the interstellar medium with theoretical models can be used to constrain electron temperature and density of the gas. An important component of the models is spontaneous transition rates between bound levels. Calculating these rates relies on accurate bound-bound oscillator strengths, which can be cast in terms of the Gaunt factor. The Gaunt factor contains terminating hypergeometric functions that cannot be calculated with sufficient accuracy for high quantum levels ($n \gtrsim 50$) by standard machine-precision methods. Methods to overcome the accuracy problem have been developed, which include asymptotic expansions and recursion relations. These methods, used in astrophysical models to calculate oscillator strengths, can introduce errors, sometimes up to as much as $\sim 8$ per cent. Detections of radio recombination lines with the new Low Frequency Array (LOFAR) has prompted an examination of theoretical models of the interstellar medium. We revisit the calculation of the Gaunt factor, employing modern arbitrary-precision computational methods to tabulate the Gaunt factor for transitions up to quantum level $n=2000$, sufficient to model low frequency Carbon radio recombination lines. The calculations provide a relative error of $\sim3\times 10^{-4}$ when compared to more detailed calculations including relativistic corrections. Our values for the Gaunt factor are provided for download in a tabular format to be used for a wide range of applications.

Observations with the Spitzer Space Telescope and the WISE satellite have revealed a prominent arc-like structure at 50" ($\simeq0.1$ pc) from the O9.5V/B0.5V system $\sigma$ Ori AB. We attribute this dust structure to the interaction of radiation pressure from the star with dust carried along by the IC 434 photo-evaporative flow of ionized gas from the dark cloud L1630. We have developed a quantitative model for the interaction of a dusty ionized flow with nearby (massive) stars where radiation pressure stalls dust, piling it up at an appreciable distance (> 0.1 pc), and force it to flow around the star. The model demonstrates that for the conditions in IC 434, the gas will decouple from the dust and will keep its original flow lines. We argue that this dust structure is the first example of a dust wave created by a massive star moving through the interstellar medium. Dust waves (and bow waves) stratify dust grains according to their radiation pressure opacity, which reflects the size distribution and composition of the grain material. Comparison of our model with observations implies that dust-gas coupling through Coulomb interaction is less important than previously thought, challenging our understanding of grain dynamics in hot, ionized regions of space. We describe the difference between dust (and bow) waves and classical bow shocks. We conclude that dust waves and bow waves should be common around stars showing the weak-wind phenomenon and that these structures are best observed at mid-IR to FIR wavelengths. In particular, dust waves and bow waves are most efficiently formed around weak-wind stars moving through a high density medium. Moreover, they provide a unique opportunity to study the direct interaction between a (massive) star and its immediate surroundings.

Massive stars drive the evolution of the interstellar medium through their radiative and mechanical energy input. After their birth, they form bubbles of hot gas surrounded by a dense shell. Traditionally, the formation of bubbles is explained through the input of a powerful stellar wind, even though direct evidence supporting this scenario is lacking. Here we explore the possibility that interstellar bubbles seen by the Spitzer- and Herschel space telescopes, blown by stars with log(L/L_sun) < 5.2, form and expand due to the thermal pressure accompanying ionization of the surrounding gas. We show that density gradients in the natal cloud or a puncture in the swept up shell lead to an ionized gas flow through the bubble into the general interstellar medium, which is traced by a dust wave near the star, demonstrating the importance of radiation pressure during this phase. Dust waves provide a natural explanation for the presence of dust inside H II bubbles, offer a novel method to study dust in H II regions and provide direct evidence that bubbles are relieving their pressure into the ISM through a champagne flow, acting as a probe of the radiative interaction of a massive star with its surroundings. We create a parameter space connecting the ambient density, the ionizing source luminosity, and the position of the dust wave, while using the well-studied H II bubbles RCW 120 and RCW 82 as benchmarks of our model. Finally, we briefly examine the implications of our study for the environments of super star clusters formed in UltraLuminous InfraRed Galaxies (ULIRGs), merging galaxies and the early Universe, which occur in very luminous and dense environments and where radiation pressure is expected to dominate the dynamical evolution.

The spectra of protoplanetary disks contain mid- and far- infrared emission features produced by forsterite dust grains. The spectral features contain information about the forsterite temperature, chemical composition and grain size. We aim to characterize how the 23 and 69 micron features can be used to constrain the physical locations of forsterite in disks. We check for consistency between two independent forsterite temperature measurements: the 23/69 feature strength ratio and the shape of the 69 micron band. We performed radiative transfer modeling to study the effect of disk properties to the forsterite spectral features. Temperature-dependent forsterite opacities were considered in self-consistent models to compute forsterite emission from protoplanetary disks. Modelling grids are presented to study the effects of grain size, disk gaps, radial mixing and optical depth to the forsterite features. Independent temperature estimates derived from the 23/69 feature strength ratio and the 69 micron band shape are most inconsistent for HD141569 and Oph IRS 48. A case study of the disk of HD141569 shows two solutions to fit the forsterite spectrum. A model with T ~ 40 K, iron-rich (~0-1 % Fe) and 1 micron forsterite grains, and a model with warmer (T ~ 100 K), iron-free, and larger (10 micron) grains. We find that for disks with low upper limits of the 69 micron feature (most notably in flat, self-shadowed disks), the forsterite must be hot, and thus close to the star. We find no correlation between disk gaps and the presence or absence of forsterite features. We argue that the 69 micron feature of the evolved transitional disks HD141569 and Oph IRS 48 is most likely a tracer of larger (i.e. ~10 micron) forsterite grains.

We report the first detections of OH$^+$ emission in planetary nebulae (PNe). As part of an imaging and spectroscopy survey of 11 PNe in the far-IR using the PACS and SPIRE instruments aboard the Herschel Space Observatory, we performed a line survey in these PNe over the entire spectral range between 51 and 672$\mu$m to look for new detections. OH$^+$ rotational emission lines at 152.99, 290.20, 308.48, and 329.77$\mu$m were detected in the spectra of three planetary nebulae: NGC 6445, NGC 6720, and NGC 6781. Excitation temperatures and column densities derived from these lines are in the range of 27 to 47 K and 2$\times$10$^{10}$ to 4 $\times$10$^{11}$ cm$^{-2}$, respectively. In PNe, the OH+ rotational line emission appears to be produced in the photodissociation region (PDR) in these objects. The emission of OH+ is observed only in PNe with hot central stars (T$_{eff}$ > 100000 K), suggesting that high-energy photons may play a role in the OH+ formation and its line excitation in these objects, as it seems to be the case for ultraluminous galaxies.

This is the first of a series of investigations into far-IR characteristics of 11 planetary nebulae (PNs) under the Herschel Space Observatory Open Time 1 program, Herschel Planetary Nebula Survey (HerPlaNS). Using the HerPlaNS data set, we look into the PN energetics and variations of the physical conditions within the target nebulae. In the present work, we provide an overview of the survey, data acquisition and processing, and resulting data products. We perform (1) PACS/SPIRE broadband imaging to determine the spatial distribution of the cold dust component in the target PNs and (2) PACS/SPIRE spectral-energy-distribution (SED) and line spectroscopy to determine the spatial distribution of the gas component in the target PNs. For the case of NGC 6781, the broadband maps confirm the nearly pole-on barrel structure of the amorphous carbon-richdust shell and the surrounding halo having temperatures of 26-40 K. The PACS/SPIRE multi-position spectra show spatial variations of far-IR lines that reflect the physical stratification of the nebula. We demonstrate that spatially-resolved far-IR line diagnostics yield the (T_e, n_e) profiles, from which distributions of ionized, atomic, and molecular gases can be determined. Direct comparison of the dust and gas column mass maps constrained by the HerPlaNS data allows to construct an empirical gas-to-dust mass ratio map, which shows a range of ratios with the median of 195+-110. The present analysis yields estimates of the total mass of the shell to be 0.86 M_sun, consisting of 0.54 M_sun of ionized gas, 0.12 M_sun of atomic gas, 0.2 M_sun of molecular gas, and 4 x 10^-3 M_sun of dust grains. These estimates also suggest that the central star of about 1.5 M_sun initial mass is terminating its PN evolution onto the white dwarf cooling track.

Planet-forming disks of gas and dust around young stars contain polycyclic aromatic hydrocarbons (PAHs). We aim to characterize how the charge state of PAHs can be used as a probe of flows of gas through protoplanetary gaps. In this context, our goal is to understand the PAH spectra of four transitional disks. In addition, we want to explain the observed correlation between PAH ionization (traced by the 6.2/11.3 feature ratio) and the disk mass (traced by the 1.3 mm luminosity). We implement a model to calculate the charge state of PAHs in the radiative transfer code MCMax. The emission spectra and ionization balance are calculated. A benchmark modeling grid is presented that shows how PAH ionization and luminosity behave as a function of star and disk properties. The PAH ionization is most sensitive to ultraviolet (UV) radiation and the electron density. In optically thick disks, where the UV field is low and the electron density is high, PAHs are predominantly neutral. Ionized PAHs trace low-density optically thin disk regions where the UV field is high and the electron density is low. Such regions are characteristic of gas flows through the gaps of transitional disks. We demonstrate that fitting the PAH spectra of four transitional disks requires a contribution of ionized PAHs in gas flows through the gap. The PAH spectra of transitional disks can be understood as superpositions of neutral and ionized PAHs. For HD97048, neutral PAHs in the optically thick disk dominate the spectrum. In the cases of HD169142, HD135344B and Oph IRS 48, small amounts of ionized PAHs located in the gas flows through the gap are strong contributors to the total PAH luminosity. The observed trend between the disk mass and PAH ionization may imply that lower-mass disks have larger gaps. Ionized PAHs in gas flows through these gaps contribute strongly to their spectra.

We present the first detection of carbon radio recombination line absorption along the line of sight to Cygnus A. The observations were carried out with the LOw Frequency ARray in the 33 to 57 MHz range. These low frequency radio observations provide us with a new line of sight to study the diffuse, neutral gas in our Galaxy. To our knowledge this is the first time that foreground Milky Way recombination line absorption has been observed against a bright extragalactic background source. By stacking 48 carbon $\alpha$ lines in the observed frequency range we detect carbon absorption with a signal-to-noise ratio of about 5. The average carbon absorption has a peak optical depth of 2$\times$10$^{-4}$, a line width of 10 km s$^{-1}$ and a velocity of +4 km s$^{-1}$ with respect to the local standard of rest. The associated gas is found to have an electron temperature $T_{e}\sim$ 110 K and density $n_{e}\sim$ 0.06 cm$^{-3}$. These properties imply that the observed carbon $\alpha$ absorption likely arises in the cold neutral medium of the Orion arm of the Milky Way. Hydrogen and helium lines were not detected to a 3$\sigma$ peak optical depth limit of 1.5$\times$10$^{-4}$ for a 4 km s$^{-1}$ channel width. Radio recombination lines associated with Cygnus A itself were also searched for, but are not detected. We set a 3$\sigma$ upper limit of 1.5$\times$10$^{-4}$ for the peak optical depth of these lines for a 4 km s$^{-1}$ channel width.

Aims. This article aims to provide an alternative method of measuring the porosity of multi-phase composite ices from their refractive indices and of characterising how the abundance of a premixed contaminant (e.g., CO2) affects the porosity of water-rich ice mixtures during omni-directional deposition. Methods. We combine optical laser interference and extended effective medium approximations (EMAs) to measure the porosity of three astrophysically relevant ice mixtures: H2O:CO2=10:1, 4:1, and 2:1. Infrared spectroscopy is used as a benchmarking test of this new laboratory-based method. Results. By independently monitoring the O-H dangling modes of the different water-rich ice mixtures, we confirm the porosities predicted by the extended EMAs. We also demonstrate that CO2 premixed with water in the gas phase does not significantly affect the ice morphology during omni-directional deposition, as long as the physical conditions favourable to segregation are not reached. We propose a mechanism in which CO2 molecules diffuse on the surface of the growing ice sample prior to being incorporated into the bulk and then fill the pores partly or completely, depending on the relative abundance and the growth temperature.

We report a Herschel detection of high-J rotational CO lines from a dense knot in the supernova remnant Cas A. Based on a combined analysis of these rotational lines, and previously observed ro-vibrational CO lines, we find the gas to be warm (two components at 400 and 2000 K) and dense (1e6-7 cm-3), with a CO column density of 5e17 cm-2. This, along with the broad line widths (400 kms-1), suggests that the CO emission originates in the post-shock region of the reverse shock. As the passage of the reverse shock dissociates any existing molecules, the CO has most likely reformed in the last few years, in the post-shock gas. The CO cooling time is comparable to the CO formation time, so possible heating sources (UV photons from the shock front, X-rays, electron conduction) to maintain the large column density of warm CO are discussed.

Infrared photometry and spectroscopy (1-25 um) of background stars reddened by the Lupus molecular cloud complex are used to determine the properties of the grains and the composition of the ices before they are incorporated into circumstellar envelopes and disks. H2O ices form at extinctions of Ak=0.25+/-0.07 mag (Av=2.1+/-0.6). Such a low ice formation threshold is consistent with the absence of nearby hot stars. Overall, the Lupus clouds are in an early chemical phase. The abundance of H2O ice (2.3+/-0.1*10^-5 relative to Nh) is typical for quiescent regions, but lower by a factor of 3-4 compared to dense envelopes of YSOs. The low solid CH3OH abundance (<3-8% relative to H2O) indicates a low gas phase H/CO ratio, which is consistent with the observed incomplete CO freeze out. Furthermore it is found that the grains in Lupus experienced growth by coagulation. The mid-infrared (>5 um) continuum extinction relative to Ak increases as a function of Ak. Most Lupus lines of sight are well fitted with empirically derived extinction curves corresponding to Rv~ 3.5 (Ak=0.71) and Rv~5.0 (Ak=1.47). For lines of sight with Ak>1.0 mag, the tau9.7/Ak ratio is a factor of 2 lower compared to the diffuse medium. Below 1.0 mag, values scatter between the dense and diffuse medium ratios. The absence of a gradual transition between diffuse and dense medium-type dust indicates that local conditions matter in the process that sets the tau9.7/Ak ratio. This process is likely related to grain growth by coagulation, as traced by the A7.4/Ak continuum extinction ratio, but not to ice mantle formation. Conversely, grains acquire ice mantles before the process of coagulation starts.

Broadband spectral surveys of protostars offer a rich view of the physical, chemical and dynamical structure and evolution of star-forming regions. The Herschel Space Observatory opened up the terahertz regime to such surveys, giving access to the fundamental transitions of many hydrides and to the high-energy transitions of many other species. A comparative analysis of the chemical inventories and physical processes and properties of protostars of various masses and evolutionary states is the goal of the Herschel CHEmical Surveys of Star forming regions (CHESS) key program. This paper focusses on the intermediate-mass protostar, OMC-2 FIR 4. We obtained a spectrum of OMC-2 FIR 4 in the 480 to 1902 GHz range with the HIFI spectrometer onboard Herschel and carried out the reduction, line identification, and a broad analysis of the line profile components, excitation, and cooling. We detect 719 spectral lines from 40 species and isotopologs. The line flux is dominated by CO, H2O, and CH3OH. The line profiles are complex and vary with species and upper level energy, but clearly contain signatures from quiescent gas, a broad component likely due to an outflow, and a foreground cloud. We find abundant evidence for warm, dense gas, as well as for an outflow in the field of view. Line flux represents 2% of the 7 L_Sol luminosity detected with HIFI in the 480 to 1250 GHz range. Of the total line flux, 60% is from CO, 13% from H2O and 9% from CH3OH. A comparison with similar HIFI spectra of other sources is set to provide much new insight into star formation regions, a case in point being a difference of two orders of magnitude in the relative contribution of sulphur oxides to the line cooling of Orion KL and OMC-2 FIR 4.

The spatial variations in PAH band intensities are normally attributed to the physical conditions of the emitting PAHs, however in recent years it has been suggested that such variations are caused mainly by extinction. To resolve this question, we have obtained near-infrared (NIR), mid-infrared (MIR) and radio observations of the compact H ii region IRAS 12063-6259. We use these data to construct multiple independent extinction maps and also to measure the main PAH features (6.2, 7.7, 8.6 and 11.2 {\mu}m) in the MIR. Three extinction maps are derived: the first using the NIR hydrogen lines and case B recombination theory; the second combining the NIR data with radio data; and the third making use of the Spitzer/IRS MIR observations to measure the 9.8 {\mu}m silicate absorption feature using the Spoon method and PAHFIT (as the depth of this feature can be related to overall extinction). The silicate absorption over the bright, southern component of IRAS 12063-6259 is almost absent while the other methods find significant extinction. While such breakdowns of the relationship between the NIR extinction and the 9.8 {\mu}m absorption have been observed in molecular clouds, they have never been observed for HII regions. We then compare the PAH intensity variations in the Spitzer/IRS data after dereddening to those found in the original data. It was found that in most cases, the PAH band intensity variations persist even after dereddening, implying that extinction is not the main cause of the PAH band intensity variations.

*Context The evolution of young massive protoplanetary disks toward planetary systems is expected to include the formation of gaps and the depletion of dust and gas. *Aims A special group of flaring disks around Herbig Ae/Be stars do not show prominent silicate emission features. We focus our attention on four key Herbig Ae/Be stars to understand the structural properties responsible for the absence of silicate feature emission. *Methods We investigate Q- and N-band images taken with Subaru/COMICS, Gemini South/T-ReCS and VLT/VISIR. Our radiative transfer modeling solutions require a separation of inner- and outer- disks by a large gap. From this we characterize the radial density structure of dust and PAHs in the disk. *Results The inner edge of the outer disk has a high surface brightness and a typical temperature between ~100-150 K and therefore dominates the emission in the Q-band. We derive radii of the inner edge of the outer disk of 34, 23, 30 and 63 AU for HD97048, HD169142, HD135344B and Oph IRS 48 respectively. For HD97048 this is the first detection of a disk gap. The continuum emission in the N-band is not due to emission in the wings of PAHs. This continuum emission can be due to VSGs or to thermal emission from the inner disk. We find that PAH emission is not always dominated by PAHs on the surface of the outer disk. *Conclusions. The absence of silicate emission features is due to the presence of large gaps in the critical temperature regime. Many, if not all Herbig disks with Spectral Energy Distribution (SED) classification `group I' are disks with large gaps and can be characterized as (pre-) transitional. An evolutionary path from the observed group I to the observed group II sources seems no longer likely. Instead, both might derive from a common ancestor.

Optical and infrared emission lines from HII regions are an important diagnostic used to study galaxies, but interpretation of these lines requires significant modeling of both the internal structure and dynamical evolution of the emitting regions. Most of the models in common use today assume that HII region dynamics are dominated by the expansion of stellar wind bubbles, and have neglected the contribution of radiation pressure to the dynamics, and in some cases also to the internal structure. However, recent observations of nearby galaxies suggest that neither assumption is justified, motivating us to revisit the question of how HII region line emission depends on the physics of winds and radiation pressure. In a companion paper we construct models of single HII regions including and excluding radiation pressure and winds, and in this paper we describe a population synthesis code that uses these models to simulate galactic collections of HII regions with varying physical parameters. We show that the choice of physical parameters has significant effects on galactic emission line ratios, and that in some cases the line ratios can exceed previously claimed theoretical limits. Our results suggest that the recently-reported offset in line ratio values between high-redshift star-forming galaxies and those in the local universe may be partially explained by the presence of large numbers of radiation pressured-dominated HII regions within them.

The emission line ratios [OIII]5007/H-beta and [NII]6584/H-alpha have been adopted as an empirical way to distinguish between the fundamentally different mechanisms of ionization in emission-line galaxies. However, detailed interpretation of these diagnostics requires calculations of the internal structure of the emitting HII regions, and these calculations depend on the assumptions one makes about the relative importance of radiation pressure and stellar winds. In this paper we construct a grid of quasi-static HII region models to explore how choices about these parameters alter HII regions' emission line ratios. We find that, when radiation pressure is included in our models, HII regions reach a saturation point beyond which further increases in the luminosity of the driving stars does not produce any further increase in effective ionization parameter, and thus does not yield any further alteration in an HII region's line ratio. We also show that, if stellar winds are assumed to be strong, the maximum possible ionization parameter is quite low. As a result of this effect, it is inconsistent to simultaneously assume that HII regions are wind-blown bubbles and that they have high ionization parameters; some popular HII region models suffer from this inconsistency. Our work in this paper provides a foundation for a companion paper in which we embed the model grids we compute here within a population synthesis code that enables us to compute the integrated line emission from galactic populations of HII regions.

Proplyds are photodissociation region (PDR)-like cometary cocoons around young stars which are thought to originate through photo-evaporation of the central protoplanetary disk by external UV radiation from the nearby OB stars. This letter presents spatially resolved mid-infrared imaging and spectroscopy of the proplyd HST10 obtained with the VLT/VISIR instrument. These observations allow us to detect Polycyclic Aromatic Hydrocarbons (PAH) emission in the proplyd photodissociation region and to study the general properties of PAHs in proplyds for the first time. We find that PAHs in HST10 are mostly neutral and at least 50 times less abundant than typical values found for the diffuse ISM or the nearby Orion Bar. With such a low PAH abundance, photoelectric heating is significantly reduced. If this low abundance pertains also to the original disk material, gas heating rates could be too low to efficiently drive photoevaporation unless other processes can be identified. Alternatively, the model behind the formation of proplyds as evaporating disks may have to be revised.

Context. The interaction between stellar winds and the interstellar medium (ISM) can create complex bow shocks. The photometers on board the Herschel Space Observatory are ideally suited to studying the morphologies of these bow shocks. Aims. We aim to study the circumstellar environment and wind-ISM interaction of the nearest red supergiant, Betelgeuse. Methods. Herschel PACS images at 70, 100, and 160 micron and SPIRE images at 250, 350, and 500 micron were obtained by scanning the region around Betelgeuse. These data were complemented with ultraviolet GALEX data, near-infrared WISE data, and radio 21 cm GALFA-HI data. The observational properties of the bow shock structure were deduced from the data and compared with hydrodynamical simulations. Results. The infrared Herschel images of the environment around Betelgeuse are spectacular, showing the occurrence of multiple arcs at 6-7 arcmin from the central target and the presence of a linear bar at 9 arcmin. Remarkably, no large-scale instabilities are seen in the outer arcs and linear bar. The dust temperature in the outer arcs varies between 40 and 140 K, with the linear bar having the same colour temperature as the arcs. The inner envelope shows clear evidence of a non-homogeneous clumpy structure (beyond 15 arcsec), probably related to the giant convection cells of the outer atmosphere. The non-homogeneous distribution of the material even persists until the collision with the ISM. A strong variation in brightness of the inner clumps at a radius of 2 arcmin suggests a drastic change in mean gas and dust density some 32 000 yr ago. Using hydrodynamical simulations, we try to explain the observed morphology of the bow shock around Betelgeuse. Conclusions: [abbreviated]

Over the last few years, the chemistry of molecules other than CO in the planet-forming zones of disks is starting to be explored with Spitzer and high-resolution ground-based data. However, these studies have focused only on a few simple molecules. The aim of this study is to put observational constraints on the presence of more complex organic and sulfur-bearing molecules predicted to be abundant in chemical models of disks and to simulate high resolution spectra in view of future missions. High S/N Spitzer spectra at 10-30 micron of the near edge-on disks IRS46 and GVTau are used to search for mid-infrared absorption bands of various molecules. These disks are good laboratories because absorption studies do not suffer from low line/continuum ratios that plague emission data. Simple LTE slab models are used to infer column densities (or upper limits) and excitation temperatures. Bands of HCN, C2H2 and CO2 are clearly detected toward both sources. The HCN and C2H2 absorption arises in warm gas with Tex of 400-700 K, whereas the CO2 absorption originates in cooler gas of app. 250 K (as in Lahuis 2006). No other absorption features are detected. Limits of those molecules are determined and compared with disk models. The inferred ratios wrt. C2H2 and HCN are roughly consistent with models of the chemistry in high-T gas. Models of UV irradiated disk surfaces generally agree better than pure X-ray models. The limit on NH3/HCN implies that evaporation of NH3 containing ices is only a minor contributor. The inferred ratios also compare well with those found in comets, suggesting that part of the cometary material may derive from warm inner disk gas. High resolution simulations show that future instruments on JWST, ELTs, SOFIA and SPICA can probe up to an order of magnitude lower ratios and put important new constraints on the models, especially if pushed to high S/N ratios.

The disk around the Herbig Ae star HD\,169142 was imaged and resolved at 18.8 and 24.5\,$\mu$m using Subaru/COMICS. We interpret the observations using a 2D radiative transfer model and find evidence for the presence of a large gap. The MIR images trace dust that emits at the onset of the strong rise in the spectral energy distribution (SED) at 20\,$\mu$m, therefore are very sensitive to the location and characteristics of the inner wall of the outer disk and its dust. We determine the location of the wall to be 23$^{+3}_{-5}$\,AU from the star. An extra component of hot dust must exist close to the star. We find that a hydrostatic optically thick inner disk does not produce enough flux in the NIR and an optically thin geometrically thick component is our solution to fit the SED. Considering the recent findings of gaps and holes in a number of Herbig Ae/Be group I disks, we suggest that such disk structures may be common in group I sources. Classification as group I should be considered a support for classification as a transitional disk, though improved imaging surveys are needed to support this speculation.

We present 75"x75" size maps of M82 at 6.4 micron, 6.6 micron, 7.7 micron, 31.5 micron, and 37.1 micron with a resolution of ~4" that we have obtained with the mid-IR camera FORCAST on SOFIA. We find strong emission from the inner 60" (~1kpc) along the major axis, with the main peak 5" west-southwest of the nucleus and a secondary peak 4" east-northeast of the nucleus. The detailed morphology of the emission differs among the bands, which is likely due to different dust components dominating the continuum emission at short mid-IR wavelengths and long mid-IR wavelengths. We include Spitzer-IRS and Herschel/PACS 70 micron data to fit spectral energy distribution templates at both emission peaks. The best fitting templates have extinctions of A_V = 18 and A_V = 9 toward the main and secondary emission peak and we estimated a color temperature of 68 K at both peaks from the 31 micron and 37 micron measurement. At the emission peaks the estimated dust masses are on the order of 10^{4} M_sun.

Buckminsterfullerene (C60) was recently confirmed to be the largest molecule identified in space. However, it remains unclear how, and where this molecule is formed. It is generally believed that C60 is formed from the build up of small carbonaceous compounds, in the hot and dense envelopes of evolved stars. Analyzing infrared observations, obtained by Spitzer and Herschel, we found that C60 is efficiently formed in the tenuous and cold environment of an interstellar cloud illuminated by strong ultraviolet (UV) radiation fields. This implies that another formation pathway, efficient at low densities, must exist. Based on recent laboratory and theoretical studies, we argue that Polycyclic Aromatic Hydrocarbons are converted into graphene, and subsequently C60, under UV irradiation from massive stars. This shows that alternative - top-down - routes are key to understanding the organic inventory in space.

We present 15-20 um spectral maps towards the reflection nebula NGC2023 obtained with the Infrared Spectrograph in short-wavelength, high-resolution mode on board the Spitzer Space Telescope. These spectra reveal emission from PAHs, C60, and H2 superposed on a dust continuum. These emission components exhibit distinct spatial distributions: with increasing distance from the illuminating star, we observe the PAH emission followed by the dust continuum emission and the H2 emission. The C60 emission is located closest to the illuminating star in the south while in the north, it seems to be associated with the H/H2 transition. Emission from PAHs and PAH-related species produce features at 15.8, 16.4, 17.4, and 17.8 um and the 15-18 um plateau. These different PAH features show distinct spatial distributions. The 15.8 um band and 15-18 um plateau correlate with the 11.2 um PAH band and with each other, and are attributed to large, neutral PAHs. Conversely, the 16.4 um feature correlates with the 12.7 um PAH band, suggesting that both arise from species that are favored by the same conditions that favor PAH cations. The PAH contribution to the 17.4 um band is displaced towards the illuminating star with respect to the 11.2 and 12.7 um emission and is assigned to doubly ionized PAHs and/or a subset of cationic PAHs. The spatial distribution of the 17.8 um band suggests it arises from both neutral and cationic PAHs. In contrast to their intensities, the profiles of the PAH bands and the 15-18 um plateau do not vary spatially. Consequently, we conclude that the carrier of the 15-18 um plateau is distinct from that of the PAH bands.

We present results of mid-infrared spectroscopic mapping observations of six star-forming regions in the Small Magellanic Cloud from the Spitzer Spectroscopic Survey of the SMC (S4MC). We detect the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) in all of the mapped regions, greatly increasing the range of environments where PAHs have been spectroscopically detected in the SMC. We investigate the variations of the mid-IR bands in each region and compare our results to studies of the PAH bands in the SINGS sample and in a sample of low-metallicity starburst galaxies. PAH emission in the SMC is characterized by low ratios of the 6-9 micron features relative to the 11.3 micron feature and weak 8.6 and 17.0 micron features. Interpreting these band ratios in the light of laboratory and theoretical studies, we find that PAHs in the SMC tend to be smaller and less ionized than those in higher metallicity galaxies. Based on studies of PAH destruction, we argue that a size distribution shifted towards smaller PAHs cannot be the result of processing in the interstellar medium, but instead reflects differences in the formation of PAHs at low metallicity. Finally, we discuss the implications of our observations for our understanding of the PAH life-cycle in low-metallicity galaxies---namely that the observed deficit of PAHs may be a consequence of PAHs forming with smaller average sizes and therefore being more susceptible to destruction under typical interstellar medium conditions.

Aims. In this paper we study the Spitzer and TIMMI2 infrared spectra of post-AGB disc sources, both in the Galaxy and the LMC. Using the observed infrared spectra we determine the mineralogy and dust parameters of the discs, and look for possible differences between the Galactic and extragalactic sources. Methods. Modelling the full spectral range observed allows us to determine the dust species present in the disc and different physical parameters such as grain sizes, dust abundance ratios, and the dust and continuum temperatures. Results. We find that all the discs are dominated by emission features of crystalline and amorphous silicate dust. Only a few sample sources show features due to CO2 gas or carbonaceous molecules such as PAHs and C60 fullerenes. Our analysis shows that dust grain processing in these discs is strong, resulting in large average grain sizes and a very high crystallinity fraction. However, we do not find any correlations between the derived dust parameters and properties of the central source. There also does not seem to be a noticeable difference between the mineralogy of the Galactic and LMC sources. Even though the observed spectra are very similar to those of protoplanetary discs around young stars, showing similar mineralogy and strong grain processing, we do find evidence for differences in the physical and chemical processes of the dust processing.

The Small Magellanic Cloud (SMC) provides a unique laboratory for the study of the lifecycle of dust given its low metallicity (~1/5 solar) and relative proximity (~60 kpc). This motivated the SAGE-SMC (Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud) Spitzer Legacy program with the specific goals of studying the amount and type of dust in the present interstellar medium, the sources of dust in the winds of evolved stars, and how much dust is consumed in star formation. This program mapped the full SMC (30 sq. deg.) including the Body, Wing, and Tail in 7 bands from 3.6 to 160 micron using the IRAC and MIPS instruments on the Spitzer Space Telescope. The data were reduced, mosaicked, and the point sources measured using customized routines specific for large surveys. We have made the resulting mosaics and point source catalogs available to the community. The infrared colors of the SMC are compared to those of other nearby galaxies and the 8 micron/24 micron ratio is somewhat lower and the 70 micron/160 micron ratio is somewhat higher than the average. The global infrared spectral energy distribution shows that the SMC has ~3X lower aromatic emission/PAH (polycyclic aromatic hydrocarbon) abundances compared to most nearby galaxies. Infrared color-magnitude diagrams are given illustrating the distribution of different asymptotic giant branch stars and the locations of young stellar objects. Finally, the average spectral energy distribution (SED) of HII/star formation regions is compared to the equivalent Large Magellanic Cloud average HII/star formation region SED. These preliminary results are expanded in detail in companion papers.

The aromatic infrared bands (AIBs) observed in the mid infrared spectrum are attributed to Polycyclic Aromatic Hydrocarbons (PAHs). We observe the NGC 7023-North West (NW) PDR in the mid-infrared (10 - 19.5 micron) using the Infrared Spectrometer (IRS), on board Spitzer. Clear variations are observed in the spectra, most notably the ratio of the 11.0 to 11.2 micron bands, the peak position of the 11.2 and 12.0 micron bands, and the degree of asymmetry of the 11.2 micron band. The observed variations appear to change as a function of position within the PDR. We aim to explain these variations by a change in the abundances of the emitting components of the PDR. A Blind Signal Separation (BSS) method, i.e. a Non-Negative Matrix Factorization algorithm is applied to separate the observed spectrum into components. Using the NASA Ames PAH IR Spectroscopic Database, these extracted signals are fit. The observed signals alone were also fit using the database and these components are compared to the BSS components. Three component signals were extracted from the observation using BSS. We attribute the three signals to ionized PAHs, neutral PAHs, and Very Small Grains (VSGs). The fit of the BSS extracted spectra with the PAH database further confirms the attribution to ionized and neutral PAHs and provides confidence in both methods for producing reliable results. The 11.0 micron feature is attributed to PAH cations while the 11.2 micron band is attributed to neutral PAHs. The VSG signal shows a characteristically asymmetric broad feature at 11.3 micron with an extended red wing. By combining the NASA Ames PAH IR Spectroscopic Database fit with the BSS method, the independent results of each method can be confirmed and some limitations of each method are overcome.

The dust size distribution in molecular clouds can be strongly affected by ice-mantle formation and (subsequent) grain coagulation. Following previous work where the dust size distribution has been calculated from a state-of-the art collision model for dust aggregates that involves both coagulation and fragmentation (Paper I), the corresponding opacities are presented in this study. The opacities are calculated by applying the effective medium theory assuming that the dust aggregates are a mix of 0.1{\mu}m silicate and graphite grains and vacuum. In particular, we explore how the coagulation affects the near-IR opacities and the opacity in the 9.7{\mu}m silicate feature. We find that as dust aggregates grow to {\mu}m-sizes both the near-IR color excess and the opacity in the 9.7 {\mu}m feature increases. Despite their coagulation, porous aggregates help to prolong the presence of the 9.7{\mu}m feature. We find that the ratio between the opacity in the silicate feature and the near-IR color excess becomes lower with respect to the ISM, in accordance with many observations of dark clouds. However, this trend is primarily a result of ice mantle formation and the mixed material composition of the aggregates, rather than being driven by coagulation. With stronger growth, when most of the dust mass resides in particles of size 10{\mu}m or larger, both the near-IR color excess and the 9.7{\mu}m silicate feature significantly diminish. Observations at additional wavelengths, in particular in the sub-mm range, are essential to provide quantitative constraints on the dust size distribution within dense cores. Our results indicate that the sub-mm index {\beta} will increase appreciably, if aggregates grow to ~100{\mu}m in size.

Mid-infrared astronomy (operating at wavelengths ranging from 2 to 25 $\mu$m) has progressed significantly in the last decades, thanks to the improvement of detector techniques and the growing diameter of telescopes. Space observatories benefit from the absence of atmospheric absorption, allowing to reach the very high sensitivities needed to perform 3D hyperspectral observations, but telescopes are limited in diameter ($< 1$ meter) and therefore provide observations at low angular resolution (typically a few seconds of arc). On the other hand, ground-based facilities suffer from strong atmospheric absorption but use large telescopes (above 8m diameter) to perform sub-arcsecond angular resolution imaging through selected windows in the mid-infrared range. In this Paper, we present a method based on Lee and Seung's Non-negative Matrix Factorization (NMF) to merge data from space and ground based mid-infrared (mid-IR) telescopes in order to combine the best sensitivity, spectral coverage and angular resolution. We prove the efficiency of this technique when applied to real mid-IR astronomical data. We suggest that this method can be applied to any combination of low and high spatial resolution positive hyperspectral datasets, as long as the spectral variety of the data allows decomposition into components using NMF.

We present a galaxy (SAGE1CJ053634.78-722658.5) at a redshift of 0.14 of which the IR is entirely dominated by emission associated with the AGN. We present the 5-37 um Spitzer/IRS spectrum and broad wavelength SED of SAGE1CJ053634, an IR point-source detected by Spitzer/SAGE (Meixner et al 2006). The source was observed in the SAGE-Spec program (Kemper et al., 2010) and was included to determine the nature of sources with deviant IR colours. The spectrum shows a redshifted (z=0.14+-0.005) silicate emission feature with an exceptionally high feature-to-continuum ratio and weak polycyclic aromatic hydrocarbon (PAH) bands. We compare the source with models of emission from dusty tori around AGNs from Nenkova et al. (2008). We present a diagnostic diagram that will help to identify similar sources based on Spitzer/MIPS and Herschel/PACS photometry. The SED of SAGE1CJ053634 is peculiar because it lacks far-IR emission and a clear stellar counterpart. We find that the SED and the IR spectrum can be understood as emission originating from the inner ~10 pc around an accreting black hole. There is no need to invoke emission from the host galaxy, either from the stars or from the interstellar medium, although a possible early-type host galaxy cannot be excluded based on the SED analysis. The hot dust around the accretion disk gives rise to a continuum, which peaks at 4 um, whereas the strong silicate features may arise from optically thin emission of dusty clouds within ~10 pc around the black hole. The weak PAH emission does not appear to be linked to star formation, as star formation templates strongly over-predict the measured far-IR flux levels. The SED of SAGE1CJ053634 is rare in the local universe but may be more common in the more distant universe. The conspicuous absence of host-galaxy IR emission places limits on the far-IR emission arising from the dusty torus alone.

While unbiased surveys observable from ground-based telescopes have previously been obtained towards several high mass protostars, very little exists on low mass protostars. To fill up this gap, we carried out a complete spectral survey of the bands at 3, 2, 1 and 0.8 mm towards the solar type protostar IRAS16293-2422. The observations covered about 200\,GHz and were obtained with the IRAM-30m and JCMT-15m telescopes. Particular attention was devoted to the inter-calibration of the obtained spectra with previous observations. All the lines detected with more than 3 sigma and free from obvious blending effects were fitted with Gaussians to estimate their basic kinematic properties. More than 4000 lines were detected (with sigma \geq 3) and identified, yielding a line density of approximatively 20 lines per GHz, comparable to previous surveys in massive hot cores. The vast majority (~2/3) of the lines are weak and due to complex organic molecules. The analysis of the profiles of more than 1000 lines belonging 70 species firmly establishes the presence of two distinct velocity components, associated with the two objects, A and B, forming the IRAS16293-2422 binary system. In the source A, the line widths of several species increase with the upper level energy of the transition, a behavior compatible with gas infalling towards a ~1 Mo object. The source B, which does not show this effect, might have a much lower central mass of ~0.1 Mo. The difference in the rest velocities of both objects is consistent with the hypothesis that the source B rotates around the source A. This spectral survey, although obtained with single-dish telescope with a low spatial resolution, allows to separate the emission from 2 different components, thanks to the large number of lines detected. The data of the survey are public and can be retrieved on the web site http://www-laog.obs.ujf-grenoble.fr/heberges/timasss.

This paper presents spectra in the 2 to 20 micron range of quiescent cloud material located in the IC 5146 cloud complex. The spectra were obtained with NASA's Infrared Telescope Facility (IRTF) SpeX instrument and the Spitzer Space Telescope's Infrared Spectrometer. We use these spectra to investigate dust and ice absorption features in pristine regions of the cloud that are unaltered by embedded stars. We find that the H2O-ice threshold extinction is 4.03+/-0.05 mag. Once foreground extinction is taken into account, however, the threshold drops to 3.2 mag, equivalent to that found for the Taurus dark cloud, generally assumed to be the touchstone quiescent cloud against which all other dense cloud and embedded young stellar object observations are compared. Substructure in the trough of the silicate band for two sources is attributed to CH3OH and NH3 in the ices, present at the ~2% and ~5% levels, respectively, relative to H2O-ice. The correlation of the silicate feature with the E(J-K) color excess is found to follow a much shallower slope relative to lines of sight that probe diffuse clouds, supporting the previous results by Chiar et al. (2007).

The recent findings of gas giant planets around young A-type stars suggest that disks surrounding Herbig Ae/Be stars will develop planetary systems. An interesting case is HD142527, for which previous observations revealed a complex circumstellar environment and an unusually high ratio of infrared to stellar luminosity. Its properties differ considerably from other Herbig Ae/Be stars. This suggests that the disk surrounding HD142527 is in an uncommon evolutionary stage. We aim for a better understanding of the geometry and evolutionary status of the circumstellar material around the Herbig Ae/Be star HD142527. We map the composition and spatial distribution of the dust around HD142527. We analyze SEST and ATCA millimeter data, VISIR N and Q-band imaging and spectroscopy. We gather additional relevant data from the literature. We use the radiative transfer code MCMax to construct a model of the geometry and density structure of the circumstellar matter, which fits all of the observables satisfactorily. We find that the disk of HD142527 has three geometrically distinct components separated by a disk gap running from 30 to 130 AU. There is a geometrically flat inner disk running from 0.3 AU up to 30 AU; an optically thin halo-like component of dust in the inner disk regions; and a massive self-shadowed outer disk running from 130 AU up to 200 AU. We derived a total dust mass in small grains of 10^-3 M0 and a vertical height of the inner wall of the outer disk of h = 60 AU. ....

The relation between ices in the envelopes and disks surrounding YSOs and those in the quiescent interstellar medium is investigated. For a sample of 31 stars behind isolated dense cores, ground-based and Spitzer spectra and photometry in the 1-25 um wavelength range are combined. The baseline for the broad and overlapping ice features is modeled, using calculated spectra of giants, H2O ice and silicates. The adopted extinction curve is derived empirically. Its high resolution allows for the separation of continuum and feature extinction. The extinction between 13-25 um is ~50% relative to that at 2.2 um. The strengths of the 6.0 and 6.85 um absorption bands are in line with those of YSOs. Thus, their carriers, which, besides H2O and CH3OH, may include NH4+, HCOOH, H2CO and NH3, are readily formed in the dense core phase, before stars form. The 3.53 um C-H stretching mode of solid CH3OH was discovered. The CH3OH/H2O abundance ratios of 5-12% are larger than upper limits in the Taurus molecular cloud. The initial ice composition, before star formation occurs, therefore depends on the environment. Signs of thermal and energetic processing that were found toward some YSOs are absent in the ices toward background stars. Finally, the peak optical depth of the 9.7 um band of silicates relative to the continuum extinction at 2.2 um is significantly shallower than in the diffuse interstellar medium. This extends the results of Chiar et al. (2007) to a larger sample and higher extinctions.

Context: PAHs are thought to be a ubiquitous and important dust component of the interstellar medium. However, the effects of their immersion in a hot (post-shock) gas have never before been fully investigated. Aims: We study the effects of energetic ion and electron collisions on PAHs in the hot post-shock gas behind interstellar shock waves. Methods: We calculate the ion-PAH and electron-PAH nuclear and electronic interactions, above the carbon atom loss threshold, in H II regions and in the hot post-shock gas, for temperatures ranging from 10^3 to 10^8 K. Results: PAH destruction is dominated by He collisions at low temperatures (T < 3x10^4 K), and by electron collisions at higher temperatures. Smaller PAHs are destroyed faster for T < 10^6 K, but the destruction rates are roughly the same for all PAHs at higher temperatures. The PAH lifetime in a tenuous hot gas (n_H ~ 0.01 cm^-3, T ~ 10^7 K), typical of the coronal gas in galactic outflows, is found to be about thousand years, orders of magnitude shorter than the typical lifetime of such objects. Conclusions: In a hot gas, PAHs are principally destroyed by electron collisions and not by the absorption of X-ray photons from the hot gas. The resulting erosion of PAHs occurs via C_2 loss from the periphery of the molecule, thus preserving the aromatic structure. The observation of PAH emission from a million degree, or more, gas is only possible if the emitting PAHs are ablated from dense, entrained clumps that have not yet been exposed to the full effect of the hot gas.

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

Context: Cosmic rays are present in almost all phases of the ISM. PAHs and cosmic rays represent an abundant and ubiquitous component of the interstellar medium. However, the interaction between them has never before been fully investigated. Aims: To study the effects of cosmic ray ion (H, He, CNO and Fe-Co-Ni) and electron bombardment of PAHs in galactic and extragalactic environments. Methods: We calculate the nuclear and electronic interactions for collisions between PAHs and cosmic ray ions and electrons with energies between 5 MeV/nucleon and 10 GeV, above the threshold for carbon atom loss, in normal galaxies, starburst galaxies and cooling flow galaxy clusters. Results: The timescale for PAH destruction by cosmic ray ions depends on the electronic excitation energy Eo and on the amount of energy available for dissociation. Small PAHs are destroyed faster, with He and the CNO group being the more effective projectiles. For electron collisions, the lifetime is independent of the PAH size and varies with the threshold energy To. Conclusions: Cosmic rays process the PAHs in diffuse clouds, where the destruction due to interstellar shocks is less efficient. In the hot gas filling galactic halos, outflows of starburst galaxies and intra-cluster medium, PAH destruction is dominated by collisions with thermal ions and electrons, but this mechanism is ineffective if the molecules are in denser cloudlets and isolated from the hot gas. Cosmic rays can access the denser clouds and together with X-rays will set the lifetime of those protected PAHs. This limits the use of PAHs as a`dye' for tracing the presence of cold entrained material.

Context: PAHs appear to be an ubiquitous interstellar dust component but the effects of shocks waves upon them have never been fully investigated. Aims: To study the effects of energetic (~0.01-1 keV) ion (H, He and C) and electron collisions on PAHs in interstellar shock waves.Methods: We calculate the ion-PAH and electron-PAH nuclear and electronic interactions, above the threshold for carbon atom loss from a PAH, in 50-200 km/s shock waves in the warm intercloud medium. Results: Interstellar PAHs (Nc = 50) do not survive in shocks with velocities greater than 100 km/s and larger PAHs (Nc = 200) are destroyed for shocks with velocities greater/equal to 125 km/s. For shocks in the ~75 - 100 km/s range, where destruction is not complete, the PAH structure is likely to be severely denatured by the loss of an important fraction (20-40%) of the carbon atoms. We derive typical PAH lifetimes of the order of a few x10^8 yr for the Galaxy. These results are robust and independent of the uncertainties in some key parameters that have yet to be well-determined experimentally. Conclusions: The observation of PAH emission in shock regions implies that that emission either arises outside the shocked region or that those regions entrain denser clumps that, unless they are completely ablated and eroded in the shocked gas, allow dust and PAHs to survive in extreme environments.

We present the comparison of the three most important ice constituents (water, CO and CO2) in the envelopes of massive Young Stellar Objects (YSOs), in environments of different metallicities: the Galaxy, the Large Magellanic Cloud (LMC) and, for the first time, the Small Magellanic Cloud (SMC). We present observations of water, CO and CO2 ice in 4 SMC and 3 LMC YSOs (obtained with Spitzer-IRS and VLT/ISAAC). While water and CO2 ice are detected in all Magellanic YSOs, CO ice is not detected in the SMC objects. Both CO and CO2 ice abundances are enhanced in the LMC when compared to high-luminosity Galactic YSOs. Based on the fact that both species appear to be enhanced in a consistent way, this effect is unlikely to be the result of enhanced CO2 production in hotter YSO envelopes as previously thought. Instead we propose that this results from a reduced water column density in the envelopes of LMC YSOs, a direct consequence of both the stronger UV radiation field and the reduced dust-to-gas ratio at lower metallicity. In the SMC the environmental conditions are harsher, and we observe a reduction in CO2 column density. Furthermore, the low gas-phase CO density and higher dust temperature in YSO envelopes in the SMC seem to inhibit CO freeze-out. The scenario we propose can be tested with further observations.

We present the classification of 197 point sources observed with the Infrared Spectrograph in the SAGE-Spec Legacy program on the Spitzer Space Telescope. We introduce a decision-tree method of object classification based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, which is used to classify the SAGE-Spec sample of point sources. The decision tree has a broad application to mid-infrared spectroscopic surveys, where supporting photometry and variability information are available. We use these classifications to make deductions about the stellar populations of the Large Magellanic Cloud and the success of photometric classification methods. We find 90 asymptotic giant branch (AGB) stars, 29 young stellar objects, 23 post-AGB objects, 19 red supergiants, eight stellar photospheres, seven background galaxies, seven planetary nebulae, two HII regions and 12 other objects, seven of which remain unclassified.

We present Herschel observations of the water molecule in the massive star-forming region W3 IRS5. The o-H17O 110-101, p-H18O 111-000, p-H2O 22 202-111, p-H2O 111-000, o-H2O 221-212, and o-H2O 212-101 lines, covering a frequency range from 552 up to 1669 GHz, have been detected at high spectral resolution with HIFI. The water lines in W3 IRS5 show well-defined high-velocity wings that indicate a clear contribution by outflows. Moreover, the systematically blue-shifted absorption in the H2O lines suggests expansion, presumably driven by the outflow. No infall signatures are detected. The p-H2O 111-000 and o-H2O 212-101 lines show absorption from the cold material (T ~ 10 K) in which the high-mass protostellar envelope is embedded. One-dimensional radiative transfer models are used to estimate water abundances and to further study the kinematics of the region. We show that the emission in the rare isotopologues comes directly from the inner parts of the envelope (T > 100 K) where water ices in the dust mantles evaporate and the gas-phase abundance increases. The resulting jump in the water abundance (with a constant inner abundance of 10^{-4}) is needed to reproduce the o-H17O 110-101 and p-H18O 111-000 spectra in our models. We estimate water abundances of 10^{-8} to 10^{-9} in the outer parts of the envelope (T < 100 K). The possibility of two protostellar objects contributing to the emission is discussed.

We present a 2Dust model for the dust shell around a LMC long-period variable (LPV) previously studied as part of the OGLE survey. OGLE LMC LPV 28579 (SAGE J051306.40-690946.3) is a carbon-rich asymptotic giant branch (AGB) star for which we have photometry and spectra from the Spitzer SAGE and SAGE-Spec programs along with UBVIJHK_s photometry. By modeling this source, we obtain a baseline set of dust properties to be used in the construction of a grid of models for carbon stars. We reproduce its spectral energy distribution using a mixture of AmC and SiC (15% by mass). The grain sizes are distributed according to the KMH model. The best-fit model has an optical depth of 0.28 for the shell at the peak of the SiC feature, with R_in~1430 R_sun or 4.4 R_star. The temperature at this inner radius is 1310 K. Assuming an expansion velocity of 10 km s^-1, we obtain a dust mass-loss rate of 2.5x10^-9 M_sun yr-1. We calculate a 15% variation in this rate by testing the fit sensitivity against variation in input parameters. We also present a simple model for the molecular gas in the extended atmosphere that could give rise to the 13.7 \mu m feature seen in the spectrum. We find that a combination of CO and C_2H_2 gas at an excitation temperature of about 1000 K and column densities of 3x10^21 cm^-2 and 10^19 cm^-2 respectively are able to reproduce the observations. Given that the excitation temperature is close to T_dust(R_in), most of the molecular contribution probably arises from the inner shell region. The luminosity corresponding to the first epoch of SAGE observations is 6580 L_sun. For an effective temperature of about 3000 K, this implies a stellar mass of 1.5-2 M_sun and an age of 1-2.5 Gyr. For a gas:dust ratio of 200, we obtain a gas mass-loss rate of 5.0x10^-7 M_sun yr^-1, consistent with the gas mass-loss rates estimated from the period, color and 8 \mu m flux of the source.

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars, NGC 1333 IRAS 2A, IRAS 4A, and IRAS 4B, obtained as part of the Water In Star-forming regions with Herschel (WISH) key program. The spectrally-resolved HIFI data are complemented by ground-based observations of lower-J CO and isotopologue lines. The 12CO 10-9 profiles are dominated by broad (FWHM 25-30 km s^-1) emission. Radiative transfer models are used to constrain the temperature of this shocked gas to 100-200 K. Several CO and 13CO line profiles also reveal a medium-broad component (FWHM 5-10 km s^-1), seen prominently in H2O lines. Column densities for both components are presented, providing a reference for determining abundances of other molecules in the same gas. The narrow C18O 9-8 lines probe the warmer part of the quiescent envelope. Their intensities require a jump in the CO abundance at an evaporation temperature around 25 K, thus providing new direct evidence for a CO ice evaporation zone around low-mass protostars.

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

We present the first results of a high-spectral-resolution survey of the carbon-rich evolved star IRC+10216 that was carried out with the HIFI spectrometer onboard Herschel. This survey covers all HIFI bands, with a spectral range from 488 to 1901GHz. In this letter we focus on the band-1b spectrum, in a spectral range 554.5-636.5GHz, where we identified 130 spectral features with intensities above 0.03 K and a signal-to-noise ratio >5. Detected lines arise from HCN, SiO, SiS, CS, CO, metal-bearing species and, surprisingly, silicon dicarbide (SiC2). We identified 55 SiC2 transitions involving energy levels between 300 and 900 K. By analysing these rotational lines, we conclude that SiC2 is produced in the inner dust formation zone, with an abundance of ~2x10^-7 relative to molecular hydrogen. These SiC2 lines have been observed for the first time in space and have been used to derive an SiC2 rotational temperature of ~204 K and a source-averaged column density of ~6.4x10^15 cm^-2. Furthermore, the high quality of the HIFI data set was used to improve the spectroscopic rotational constants of SiC2.

In this paper we present mid-infrared spectra of a comprehensive set of Herbig Ae/Be stars observed with the Spitzer Space Telescope. The signal-to-noise ratio of these spectra is very high, ranging between about a hundred and several hundreds. During the analysis of these data we tested the validity of standard protoplanetary dust models and studied grain growth and crystal formation. On the basis of the analyzed spectra, the major constituents of protoplanetary dust around Herbig Ae/Be stars are amorphous silicates with olivine and pyroxene stoichiometry, crystalline forsterite and enstatite and silica. No other solid state features, indicating other abundant dust species, are present in the Spitzer spectra. Deviations of the synthetic spectra from the observations are most likely related to grain shape effects and uncertainties in the iron content of the dust grains. Our analysis revealed that larger grains are more abundant in the disk atmosphere of flatter disks than in that of flared disks, indicating that grain growth and sedimentation decrease the disk flaring. We did not find, however, correlations between the value of crystallinity and any of the investigated system parameters. Our analysis shows that enstatite is more concentrated toward the warm inner disk than forsterite, in contrast to predictions of equilibrium condensation models. None of the three crystal formation mechanisms proposed so far can alone explain all our findings. It is very likely that all three play at least some role in the formation of crystalline silicates.

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines is significantly diminished in these bands compared to lower frequency Herschel/HIFI bands.

'Water In Star-forming regions with Herschel' (WISH) is a key programme dedicated to studying the role of water and related species during the star-formation process and constraining the physical and chemical properties of young stellar objects. The Heterodyne Instrument for the Far-Infrared (HIFI) on the Herschel Space Observatory observed three deeply embedded protostars in the low-mass star-forming region NGC1333 in several H2-16O, H2-18O, and CO transitions. Line profiles are resolved for five H16O transitions in each source, revealing them to be surprisingly complex. The line profiles are decomposed into broad (>20 km/s), medium-broad (~5-10 km/s), and narrow (<5 km/s) components. The H2-18O emission is only detected in broad 1_10-1_01 lines (>20 km/s), indicating that its physical origin is the same as for the broad H2-16O component. In one of the sources, IRAS4A, an inverse P Cygni profile is observed, a clear sign of infall in the envelope. From the line profiles alone, it is clear that the bulk of emission arises from shocks, both on small (<1000 AU) and large scales along the outflow cavity walls (~10 000 AU). The H2O line profiles are compared to CO line profiles to constrain the H2O abundance as a function of velocity within these shocked regions. The H2O/CO abundance ratios are measured to be in the range of ~0.1-1, corresponding to H2O abundances of ~10-5-10-4 with respect to H2. Approximately 5-10% of the gas is hot enough for all oxygen to be driven into water in warm post-shock gas, mostly at high velocities.

We performed a sensitive search for the ground-state emission lines of ortho- and para-water vapor in the DM Tau protoplanetary disk using the Herschel/HIFI instrument. No strong lines are detected down to 3sigma levels in 0.5 km/s channels of 4.2 mK for the 1_{10}--1_{01} line and 12.6 mK for the 1_{11}--0_{00} line. We report a very tentative detection, however, of the 1_{10}--1_{01} line in the Wide Band Spectrometer, with a strength of T_{mb}=2.7 mK, a width of 5.6 km/s and an integrated intensity of 16.0 mK km/s. The latter constitutes a 6sigma detection. Regardless of the reality of this tentative detection, model calculations indicate that our sensitive limits on the line strengths preclude efficient desorption of water in the UV illuminated regions of the disk. We hypothesize that more than 95-99% of the water ice is locked up in coagulated grains that have settled to the midplane.

We present initial results from the Herschel GT key program: Herschel observations of EXtra-Ordinary Sources (HEXOS) and outline the promise and potential of spectral surveys with Herschel/HIFI. The HIFI instrument offers unprecedented sensitivity, as well as continuous spectral coverage across the gaps imposed by the atmosphere, opening up a largely unexplored wavelength regime to high-resolution spectroscopy. We show the spectrum of Orion KL between 480 and 560 GHz and from 1.06 to 1.115 THz. From these data, we confirm that HIFI separately measures the dust continuum and spectrally resolves emission lines in Orion KL. Based on this capability we demonstrate that the line contribution to the broad-band continuum in this molecule-rich source is ~20-40% below 1 THz and declines to a few percent at higher frequencies. We also tentatively identify multiple transitions of HD18O in the spectra. The first detection of this rare isotopologue in the interstellar medium suggests that HDO emission is optically thick in the Orion hot core with HDO/H2O ~ 0.02. We discuss the implications of this detection for the water D/H ratio in hot cores.

We performed Herschel/HIFI observations of several CO lines in the far-infrared/sub-mm in the protoplanetary nebula CRL618. The high spectral resolution provided by HIFI allows measurement of the line profiles. Since the dynamics and structure of the nebula is well known from mm-wave interferometric maps, it is possible to identify the contributions of the different nebular components (fast bipolar outflows, double shells, compact slow shell) to the line profiles. The observation of these relatively high-energy transitions allows an accurate study of the excitation conditions in these components, particularly in the warm ones, which cannot be properly studied from the low-energy lines. The 12CO J=16-15, 10-9, and 6-5 lines are easily detected in this source. 13CO J=10-9 and 6-5 are also detected. Wide profiles showing spectacular line wings have been found, particularly in 12CO 16-15. Other lines observed simultaneously with CO are also shown. Our analysis of the CO high-J transitions, when compared with the existing models, confirms the very low expansion velocity of the central, dense component, which probably indicates that the shells ejected during the last AGB phases were driven by radiation pressure under a regime of maximum transfer of momentum. No contribution of the diffuse halo found from mm-wave data is identified in our spectra, because of its low temperature. We find that the fast bipolar outflow is quite hot, much hotter than previously estimated; for instance, gas flowing at 100 km/s must have a temperature higher than ~ 200 K. Probably, this very fast outflow, with a kinematic age < 100 yr, has been accelerated by a shock and has not yet cooled down. The double empty shell found from mm-wave mapping must also be relatively hot, in agreement with the previous estimate.

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

We report the detection of absorption by interstellar hydroxyl cations and water cations, along the sight-line to the bright continuum source W49N. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 972 GHz N = 1 - 0 transition of OH+ and the 1115 GHz 1(11) - 0(00) transition of ortho-H2O+. The resultant spectra show absorption by ortho-H2O+, and strong absorption by OH+, in foreground material at velocities in the range 0 to 70 km/s with respect to the local standard of rest. The inferred OH+/H2O+ abundance ratio ranges from ~ 3 to ~ 15, implying that the observed OH+ arises in clouds of small molecular fraction, in the 2 - 8% range. This conclusion is confirmed by the distribution of OH+ and H2O+ in Doppler velocity space, which is similar to that of atomic hydrogen, as observed by means of 21 cm absorption measurements, and dissimilar from that typical of other molecular tracers. The observed OH+/H abundance ratio of a few E-8 suggests a cosmic ray ionization rate for atomic hydrogen of (0.6 - 2.4) E-16 s-1, in good agreement with estimates inferred previously for diffuse clouds in the Galactic disk from observations of interstellar H3+ and other species.

We present high resolution HIFI spectroscopy of the nucleus of the archetypical starburst galaxy M82. Six 12CO lines, 2 13CO lines and 4 fine-structure lines are detected. Besides showing the effects of the overall velocity structure of the nuclear region, the line profiles also indicate the presence of multiple components with different optical depths, temperatures and densities in the observing beam. The data have been interpreted using a grid of PDR models. It is found that the majority of the molecular gas is in low density (n=10^3.5 cm^-3) clouds, with column densities of N_H=10^21.5 cm^-2 and a relatively low UV radiation field (GO = 10^2). The remaining gas is predominantly found in clouds with higher densities (n=10^5 cm^-3) and radiation fields (GO = 10^2.75), but somewhat lower column densities (N_H=10^21.2 cm^-2). The highest J CO lines are dominated by a small (1% relative surface filling) component, with an even higher density (n=10^6 cm^-3) and UV field (GO = 10^3.25). These results show the strength of multi-component modeling for the interpretation of the integrated properties of galaxies.

We present spectrally resolved observations of triatomic carbon (C3) in several ro-vibrational transitions between the vibrational ground state and the low-energy nu2 bending mode at frequencies between 1654-1897 GHz along the sight-lines to the submillimeter continuum sources W31C and W49N, using Herschel's HIFI instrument. We detect C3 in absorption arising from the warm envelope surrounding the hot core, as indicated by the velocity peak position and shape of the line profile. The sensitivity does not allow to detect C3 absorption due to diffuse foreground clouds. From the column densities of the rotational levels in the vibrational ground state probed by the absorption we derive a rotation temperature (T_rot) of ~50--70 K, which is a good measure of the kinetic temperature of the absorbing gas, as radiative transitions within the vibrational ground state are forbidden. It is also in good agreement with the dust temperatures for W31C and W49N. Applying the partition function correction based on the derived T_rot, we get column densities N(C3) ~7-9x10^{14} cm^{-2} and abundance x(C3)~10^{-8} with respect to H2. For W31C, using a radiative transfer model including far-infrared pumping by the dust continuum and a temperature gradient within the source along the line of sight we find that a model with x(C3)=10^{-8}, T_kin=30-50 K, N(C3)=1.5 10^{15} cm^{-2} fits the observations reasonably well and provides parameters in very good agreement with the simple excitation analysis.

Recently, we have discovered an error in our Monte-Carlo spectral fitting routine, more specifically where the errors on the fluxes were rescaled to get a reduced chi2 of 1. The rescaled errors were too big, resulting in too wide a range of good fits in our 100 step Monte-Carlo routine. This problem affects Figs. 7-9 and Tables A.1, A.2 in Gielen et al. (2008), Table 3 in Gielen et al. (2009a), and Table 4 in Gielen et al. (2009b). We corrected for this error and present the new values and errors in the tables below. The new values and errors nearly all fall within the old error range. Our best chi2 values and overall former scientific results are not affected. With these new errors some possible new trends in the dust parameters might be observed. These will be discussed in an upcoming paper where we extend the sample presented in Gielen et al. (2008) with newly obtained SPITZER-IRS data.

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

The chemistry of astronomical hydrocarbons, responsible for the well-known infrared emission features detected in a wide variety of targets, remains enigmatic. Here we focus on the group of young intermediate-mass Herbig Ae stars. We have analyzed the aliphatic and polycyclic aromatic hydrocarbon (PAH) emission features in the infrared spectra of a sample of 53 Herbig Ae stars, obtained with the Infrared Spectrograph aboard the Spitzer Space Telescope. We confirm that the PAH-to-stellar luminosity ratio is higher in targets with a flared dust disk. However, a few sources with a flattened dust disk still show relatively strong PAH emission. Since PAH molecules trace the gas disk, this indicates that gas disks may still be flared, while the dust disk has settled due to grain growth. There are indications that the strength of the 11.3 um feature also depends on dust disk structure, with flattened disks being less bright in this feature. We confirm that the CC bond features at 6.2 and 7.8 um shift to redder wavelengths with decreasing stellar effective temperature. Moreover, we show that this redshift is accompanied by a relative increase of aliphatic CH emission and a decrease of the aromatic 8.6 um CH feature strength. Cool stars in our sample are surrounded by hydrocarbons with a high aliphatic/aromatic CH ratio and a low aromatic CH/CC ratio, and vice versa for the hot stars. We conclude that, while the overall hydrocarbon emission strength depends on the dust disk's geometry, the relative differences seen in the IR emission features in disks around Herbig Ae stars are mainly due to chemical differences of the hydrocarbon molecules induced by the stellar UV field. Strong UV flux reduces the aliphatic component and emphasizes the spectral signature of the aromatic molecules in the IR spectra.

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory's PACS and SPIRE photometry cameras, we imaged a 2x8 square degree strip through the LMC, at a position angle of ~22.5 degrees as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 {\mu}m and SPIRE 250, 350 and 500 {\mu}m. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 {\mu}m emission is in excess of the models by 6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_PAH) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J=1-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65(+15,-18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287(+25,-42) that is more consistent with a GDR inferred from extinction.

We report the detection of absorption lines by the reactive ions OH+, H2O+ and H3O+ along the line of sight to the submillimeter continuum source G10.6$-$0.4 (W31C). We used the Herschel HIFI instrument in dual beam switch mode to observe the ground state rotational transitions of OH+ at 971 GHz, H2O+ at 1115 and 607 GHz, and H3O+ at 984 GHz. The resultant spectra show deep absorption over a broad velocity range that originates in the interstellar matter along the line of sight to G10.6$-$0.4 as well as in the molecular gas directly associated with that source. The OH+ spectrum reaches saturation over most velocities corresponding to the foreground gas, while the opacity of the H2O+ lines remains lower than 1 in the same velocity range, and the H3O+ line shows only weak absorption. For LSR velocities between 7 and 50 kms$^{-1}$ we estimate total column densities of $N$(OH+) $> 2.5 \times 10^{14}$ cm$^{-2}$, $N$(H2O+) $\sim 6 \times 10^{13}$ cm$^{-2}$ and $N$(H3O+) $\sim 4.0 \times 10^{13}$ cm$^{-2}$. These detections confirm the role of O$^+$ and OH$^+$ in initiating the oxygen chemistry in diffuse molecular gas and strengthen our understanding of the gas phase production of water. The high ratio of the OH+ by the H2O+ column density implies that these species predominantly trace low-density gas with a small fraction of hydrogen in molecular form.

We present clumps of dust emission from Herschel observations of the Large Magellanic Cloud (LMC) and their physical and statistical properties. We catalog cloud features seen in the dust emission from Herschel observations of the LMC, the Magellanic type irregular galaxy closest to the Milky Way, and compare these features with HI catalogs from the ATCA+Parkes HI survey. Using an automated cloud-finding algorithm, we identify clouds and clumps of dust emission and examine the cumulative mass distribution of the detected dust clouds. The mass of cold dust is determined from physical parameters that we derive by performing spectral energy distribution fits to 250, 350, and 500 micronm emission from SPIRE observations using DUSTY and GRASIL radiative transfer calculation with dust grain size distributions for graphite/silicate in low-metallicity extragalactic environments. The dust cloud mass spectrum follows a power law distribution with an exponent of gamma=-1.8 for clumps larger than 400 solar mass and is similar to the HI mass distribution. This is expected from the theory of ISM structure in the vicinity of star formation.

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158\mu\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (>100 micron). Only with SPIRE on Herschel has it been possible to get high spatial resolution at 200 to 500 micron that is beyond the peak (~160 micron) of dust emission in most galaxies. We investigate the differences in the fitted dust temperatures and masses determined using only <200 micron data and then also including >200 micron data (new SPIRE observations) to determine how important having >200 micron data is for deriving these dust properties. We fit the 100 to 350 micron observations of the Large Magellanic Cloud (LMC) point-by-point with a model that consists of a single temperature and fixed emissivity law. The data used are existing observations at 100 and 160 micron (from IRAS and Spitzer) and new SPIRE observations of 1/4 of the LMC observed for the HERITAGE Key Project as part of the Herschel Science Demonstration phase. The dust temperatures and masses computed using only 100 and 160 micron data can differ by up to 10% and 36%, respectively, from those that also include the SPIRE 250 & 350 micron data. We find that an emissivity law proportional to lambda^-1.5 minimizes the 100-350 micron fractional residuals. We find that the emission at 500 micron is ~10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the fractional 500 micron excess is weakly anti-correlated with MIPS 24 micron flux and the total gas surface density. This argues against a flux calibration error as the origin of the 500 micron excess. Our results do not allow us to distinguish between a systematic variation in the wavelength dependent emissivity law or a population of very cold dust only detectable at lambda > 500 micron for the origin of the 500 micron excess.

Context. Based on the far infrared excess the Herbig class of stars is divided into a group with flaring circumstellar disks (group I) and a group with flat circumstellar disks (group II). Dust sedimentation is generally proposed as an evolution mechanism to transform flaring disks into flat disks. Theory predicts that during this process the disks preserve their gas content, however observations of group II Herbig Ae stars demonstrate a lack of gas. Aims. We map the spatial distribution of the gas and dust around the group II Herbig Ae star HD 95881. Methods. We analyze optical photometry, Q-band imaging, infrared spectroscopy, and K and N-band interferometric spectroscopy. We use a Monte Carlo radiative transfer code to create a model for the density and temperature structure which quite accurately reproduces all the observables. Results. We derive a consistent picture in which the disk consists of a thick puffed up inner rim and an outer region which has a flaring gas surface and is relatively void of 'visible' dust grains. Conclusions. HD 95881 is in a transition phase from a gas rich flaring disk to a gas poor self-shadowed disk.

The SAGE-Spec Spitzer Legacy program is a spectroscopic follow-up to the SAGE-LMC photometric survey of the Large Magellanic Cloud carried out with the Spitzer Space Telescope. We present an overview of SAGE-Spec and some of its first results. The SAGE-Spec program aims to study the life cycle of gas and dust in the Large Magellanic Cloud, and to provide information essential to the classification of the point sources observed in the earlier SAGE-LMC photometric survey. We acquired 224.6 hours of observations using the InfraRed Spectrograph and the SED mode of the Multiband Imaging Photometer for Spitzer. The SAGE-Spec data, along with archival Spitzer spectroscopy of objects in the Large Magellanic Cloud, are reduced and delivered to the community. We discuss the observing strategy, the specific data reduction pipelines applied and the dissemination of data products to the scientific community. Initial science results include the first detection of an extragalactic "21 um" feature towards an evolved star and elucidation of the nature of disks around RV Tauri stars in the Large Magellanic Cloud. Towards some young stars, ice features are observed in absorption. We also serendipitously observed a background quasar, at a redshift of z~0.14, which appears to be host-less.

The aim of this work is to link the broad lambda 5450 diffuse interstellar band (DIB) to a laboratory spectrum recorded through an expanding acetylene plasma. Cavity ring-down direct absorption spectra and astronomical observations of HD 183143 with the HERMES spectrograph on the Mercator Telescope in La Palma and the McKellar spectrograph on the DAO 1.2 m Telescope are compared. In the 543-547 nm region a broad band is measured with a band maximum at 545 nm and FWHM of 1.03(0.1) nm coinciding with a well-known diffuse interstellar band at lambda 5450 with FWHM of 0.953 nm. A coincidence is found between the laboratory and the two independent observational studies obtained at higher spectral resolution. This result is important, as a match between a laboratory spectrum and a - potentially lifetime broadened - DIB is found. A series of additional experiments has been performed in order to unambiguously identify the laboratory carrier of this band. This has not been possible. The laboratory results, however, restrict the carrier to a molecular transient, consisting of carbon and hydrogen.

We present spectroscopic observations of a sample of 15 embedded young stellar objects (YSOs) in the Large Magellanic Cloud (LMC). These observations were obtained with the Spitzer Infrared Spectrograph (IRS) as part of the SAGE-Spec Legacy program. We analyze the two prominent ice bands in the IRS spectral range: the bending mode of CO_2 ice at 15.2 micron and the ice band between 5 and 7 micron that includes contributions from the bending mode of water ice at 6 micron amongst other ice species. The 5-7 micron band is difficult to identify in our LMC sample due to the conspicuous presence of PAH emission superimposed onto the ice spectra. We identify water ice in the spectra of two sources; the spectrum of one of those sources also exhibits the 6.8 micron ice feature attributed to ammonium and methanol. We model the CO_2 band in detail, using the combination of laboratory ice profiles available in the literature. We find that a significant fraction (> 50%) of CO_2 ice is locked in a water-rich component, consistent with what is observed for Galactic sources. The majority of the sources in the LMC also require a pure-CO_2 contribution to the ice profile, evidence of thermal processing. There is a suggestion that CO_2 production might be enhanced in the LMC, but the size of the available sample precludes firmer conclusions. We place our results in the context of the star formation environment in the LMC.

Methanol and its precursor formaldehyde are among the most studied organic molecules in the interstellar medium and are abundant in the gaseous and solid phases. We recently developed a model to simulate CO hydrogenation via H atoms on interstellar ice surfaces, the most important interstellar route to H2CO and CH3OH, under laboratory conditions. We extend this model to simulate the formation of both organic species under interstellar conditions, including freeze-out from the gas and hydrogenation on surfaces. Our aim is to compare calculated abundance ratios with observed values and with the results of prior models. Simulations under different conditions, including density and temperature, have been performed. We find that H2CO and CH3OH form efficiently in cold dense cores or the cold outer envelopes of young stellar objects. The grain mantle is found to have a layered structure with CH3OH on top. The species CO and H2CO are found to exist predominantly in the lower layers of ice mantles where they are not available for hydrogenation at late times.

Aims: By studying the photospheric abundances of 4 RV Tauri stars in the LMC, we test whether the depletion pattern of refractory elements, seen in similar Galactic sources, is also common for extragalactic sources. Since this depletion process probably only occurs through interaction with a stable disc, we investigate the circumstellar environment of these sources. Methods: A detailed photospheric abundance study was performed using high-resolution UVES optical spectra. To study the circumstellar environment we use photometric data to construct the spectral energy distributions of the stars, and determine the geometry of the circumstellar environment, whereas low-resolution Spitzer-IRS infrared spectra are used to trace its mineralogy. Results: Our results show that, also in the LMC, the photospheres of RV Tauri stars are commonly affected by the depletion process, although it can differ significantly in strength from source to source. From our detailed disc modelling and mineralogy study, we find that this process, as in the Galaxy, appears closely related to the presence of a stable Keplerian disc. The newly studied extragalactic objects have similar observational characteristics as Galactic post-AGB binaries surrounded by a dusty disc, and are therefore also believed to be part of a binary system. One source shows a very small infrared excess, atypical for a disc source, but still has evidence for depletion. We speculate this could point to the presence of a very evolved disc, similar to debris discs seen around young stellar objects.

Thanks to high sensitivity and angular resolution and broad spectral coverage, SPICA will offer a unique opportunity to better characterize the nature of polycyclic aromatic hydrocarbons (PAHs) and very small grains (VSGs), to better use them as probes of astrophysical environments. The angular resolution will enable to probe the chemical frontiers in the evolution process from VSGs to neutral PAHs, to ionized PAHs and to "Grand-PAHs" in photodissotiation regions and HII regions, as a function of G$_0$/n (UV radiation field / density). High sensitivity will favor the detection of the far-IR skeletal emission bands of PAHs, which provide specific fingerprints and could lead to the identification of individual PAHs. This overall characterization will allow to use PAH and VSG populations as tracers of physical conditions in spatially resolved protoplanetary disks and nearby galaxies (using mid-IR instruments), and in high redshift galaxies (using the far-IR instrument), thanks to the broad spectral coverage SPICA provides. Based on our previous experience with ISO and Spitzer we discuss how these goals can be reached.

[abridged] We present 52-93 micron spectra obtained with Spitzer in the MIPS-SED mode, of a representative sample of luminous compact far-IR sources in the LMC. These include carbon stars, OH/IR AGB stars, post-AGB objects and PNe, RCrB-type star HV2671, OH/IR red supergiants WOHG064 and IRAS05280-6910, B[e] stars IRAS04530-6916, R66 and R126, Wolf-Rayet star Brey3a, Luminous Blue Variable R71, supernova remnant N49, a large number of young stellar objects, compact HII regions and molecular cores, and a background galaxy (z~0.175). We use the spectra to constrain the presence and temperature of cold dust and the excitation conditions and shocks within the neutral and ionized gas, in the circumstellar environments and interfaces with the surrounding ISM. Evolved stars, including LBV R71, lack cold dust except in some cases where we argue that this is swept-up ISM. This leads to an estimate of the duration of the prolific dust-producing phase ("superwind") of several thousand years for both RSGs and massive AGB stars, with a similar fractional mass loss experienced despite the different masses. We tentatively detect line emission from neutral oxygen in the extreme RSG WOHG064, with implications for the wind driving. In N49, the shock between the supernova ejecta and ISM is revealed by its strong [OI] 63-micron emission and possibly water vapour; we estimate that 0.2 Msun of ISM dust was swept up. Some of the compact HII regions display pronounced [OIII] 88-micron emission. The efficiency of photo-electric heating in the interfaces of ionized gas and molecular clouds is estimated at 0.1-0.3%. We confirm earlier indications of a low nitrogen content in the LMC. Evidence for solid state emission features is found in both young and evolved object; some of the YSOs are found to contain crystalline water ice.

An updated Science Vision for the SOFIA project is presented, including an overview of the characteristics and capabilities of the observatory and first generation instruments. A primary focus is placed on four science themes: 'The Formation of Stars and Planets', 'The Interstellar Medium of the Milky Way', 'Galaxies and the Galactic Center' and 'Planetary Science'.

We present IRS spectra and revised MIPS photometry for the 18 members of the Eta Chamaeleontis cluster. Aged 8 Myr, the Eta Cha cluster is one of the few nearby regions within the 5-10 Myr age range, during which the disk fraction decreases dramatically and giant planet formation must come to an end. For the 15 low-mass members, we measure a disk fraction ~50%, high for their 8 Myr age, and 4 of the 8 disks lack near-IR excesses, consistent with the empirical definition of "transition'' disks. Most of the disks are comparable to geometrically flat disks. The comparison with regions of different ages suggests that at least some of the "transition" disks may represent the normal type of disk around low-mass stars. Therefore, their flattened structure and inner holes may be related to other factors (initial masses of the disk and the star, environment, binarity), rather than to pure time evolution. We analyze the silicate dust in the disk atmosphere, finding moderate crystalline fractions (~10-30%) and typical grain sizes ~1-3 micron, without any characteristic trend in the composition. These results are common to other regions of different ages, suggesting that the initial grain processing occurs very early in the disk lifetime (<1 Myr). Large grain sizes in the disk atmosphere cannot be used as a proxy for age, but are likely related to higher disk turbulence. The dust mineralogy varies between the 8-12micron and the 20-30 micron features, suggesting high temperature dust processing and little radial mixing. Finally, the analysis of IR and optical data on the B9 star Eta Cha reveals that it is probably surrounded by a young debris disk with a large inner hole, instead of being a classical Be star.

In dense molecular clouds collisions between dust grains alter the ISM-dust size distribution. We study this process by inserting the results from detailed numerical simulations of two colliding dust aggregates into a coagulation model that computes the dust size distribution with time. All collisional outcomes -- sticking, fragmentation (shattering, breakage, and erosion) -- are included and the effects on the internal structure of the aggregates are also tabulated. The dust aggregate evolution model is applied to an homogeneous and static cloud of temperature 10 K and gas densities between 10^3 and 10^7 cm^-3. The coagulation is followed locally on timescales of ~10^7 yr. We find that the growth can be divided into two stages: a growth dominated phase and a fragmentation dominated phase. Initially, the mass distribution is relatively narrow and shifts to larger sizes with time. At a certain point, dependent on the material properties of the grains as well as on the gas density, collision velocities will become sufficiently energetic to fragment particles, halting the growth and replenishing particles of lower mass. Eventually, a steady state is reached, where the mass distribution is characterized by a mass spectrum of approximately equal amount of mass per logarithmic size bin. The amount of growth that is achieved depends on the cloud's lifetime. If clouds exist on free-fall timescales the effects of coagulation on the dust size distribution are very minor. On the other hand, if clouds have long-term support mechanisms, the impact of coagulation is important, resulting in a significant decrease of the opacity on timescales longer than the initial collision timescale between big grains.

We present a Spitzer Space Telescope spectroscopic study of a sample of 25 planetary nebulae in the Magellanic Clouds. The low-resolution modules are used to analyze the dust features present in the infrared spectra. This study complements a previous work by the same authors where the same sample was analyzed in terms of neon and sulfur abundances. Over half of the objects (14) show emission of polycyclic aromatic hydrocarbons, typical of carbon-rich dust environments. We compare the hydrocarbon emission in our objects to those of Galactic HII regions and planetary nebulae, and LMC/SMC HII regions. Amorphous silicates are seen in just two objects, enforcing the now well-known-fact that oxygen-rich dust is less common at low metallicities. Besides these common features, some planetary nebulae show very unusual dust. Nine objects show a strong silicon carbide feature at 11um and twelve of them show magnesium sulfide emission starting at 25um. The high percentage of spectra with silicon carbide in the Magellanic Clouds is not common. Two objects show a broad band which may be attributed to hydrogenated amorphous carbon and weak low-excitation atomic lines. It is likely that these nebulae are very young. The spectra of the remaining eight nebulae are dominated by the emission of fine-structure lines with a weak continuum due to thermal emission of dust, although in a few cases the S/N in the spectra is low, and weak dust features may not have been detected.

This white paper is submitted to the Astronomy and Astrophysics 2010 Decadal Survey (Astro2010)1 Committee on the State of the Profession to emphasize the potential of the Stratospheric Observatory for Infrared Astronomy (SOFIA) to contribute to the training of instrumentalists and observers, and to related technology developments. This potential goes beyond the primary mission of SOFIA, which is to carry out unique, high priority astronomical research. SOFIA is a Boeing 747SP aircraft with a 2.5 meter telescope. It will enable astronomical observations anywhere, any time, and at most wavelengths between 0.3 microns and 1.6 mm not accessible from ground-based observatories. These attributes, accruing from the mobility and flight altitude of SOFIA, guarantee a wealth of scientific return. Its instrument teams (nine in the first generation) and guest investigators will do suborbital astronomy in a shirt-sleeve environment. The project will invest $10M per year in science instrument development over a lifetime of 20 years. This, frequent flight opportunities, and operation that enables rapid changes of science instruments and hands-on in-flight access to the instruments, assure a unique and extensive potential - both for training young instrumentalists and for encouraging and deploying nascent technologies. Novel instruments covering optical, infrared, and submillimeter bands can be developed for and tested on SOFIA by their developers (including apprentices) for their own observations and for those of guest observers, to validate technologies and maximize observational effectiveness.

We present a mid-infrared spectroscopic data cube of the central part of 30 Doradus, observed with Spitzer's IRS and MIPS/SED mode. Aromatic dust emission features and emission lines from molecular and atomic hydrogen are detected but not particularly strong. The dominant spectral features are emission lines from moderately ionized species of argon, neon, and sulphur, which are used to determine the physical conditions in the ionized gas. The ionized gas excitation shows strong variations on parsec scales, some of which can plausibly be associated with individual hot stars. We fit the ionic line strengths with photoionization and shock models, and find that photoionization dominates in the region. The ionization parameter U traces the rim of the central bubble, as well as highlighting isolated sources of ionization, and at least one quiescent clump. The hardness of the ionizing radiation field T_rad reveals several "hot spots" that are either the result of individual very hot stars or trace the propagation of the diffuse ionizing field through the surrounding neutral cloud. Consistent with other measurements of giant molecular hydrogen regions, log(U) ranges between -3 and -0.75, and T_rad between 30000 and 85000K.

(abridged) We study photoelectric heating throughout the Large Magellanic Cloud. We quantify the importance of the [CII] cooling line and the photoelectric heating process of various environments in the LMC and investigate which parameters control the extent of photoelectric heating. We use the BICE [CII] map and the Spitzer/SAGE infrared maps. We examine the spatial variations in the efficiency of photoelectric heating: photoelectric heating rate over power absorbed by grains. We correlate the photoelectric heating efficiency and the emission from various dust constituents and study the variations as a function of H\alpha emission, dust temperatures, and the total infrared luminosity. From this we estimate radiation field, gas temperature, and electron density. We find systematic variations in photoelectric efficiency. The highest efficiencies are found in the diffuse medium, while the lowest coincide with bright star-forming regions (~1.4 times lower). The [CII] line emission constitutes 1.32% of the far infrared luminosity across the whole of the LMC. We find correlations between the [CII] emission and ratios of the mid infrared and far infrared bands, which comprise various dust constituents. The correlations are interpreted in light of the spatial variations of the dust abundance and by the local environmental conditions that affect the dust emission properties. As a function of the total infrared surface brightness, S_{TIR}, the [CII] surface brightness can be described as: S_{[CII]}=1.25 S_{TIR}^{0.69} [10^{-3} erg s^{-1} cm^{-2} sr^{-1}]. The [CII] emission is well-correlation with the 8 micrometer emission, suggesting that the polycyclic aromatic hydrocarbons play a dominant role in the photoelectric heating process.

Circumstellar disks provide the material reservoir for the growth of young stars and for planet formation. We combine a high-level radiative transfer program with a thermal-chemical model of a typical T Tauri star disk to investigate the diagnostic potential of the far-infrared lines of water for probing disk structure. We discuss the observability of pure rotational H2O lines with the Herschel Space Observatory, specifically the residual gas where water is mainly frozen out. We find that measuring both the line profile of the ground 110-101 ortho-H2O transition and the ratio of this line to the 312-303 and 221-212 line can provide information on the gas phase water between 5-100 AU, but not on the snow line which is expected to occur at smaller radii.

Astronomical mid-IR spectra show two minor PAH features at 5.25 and 5.7 $\mu$m (1905 and 1754 cm$^{\rm - 1}$) that hitherto have been little studied, but contain information about the astronomical PAH population that complements that of the major emission bands. Here we report a study involving both laboratory and theoretical analysis of the fundamentals of PAH spectroscopy that produce features in this region and use these to analyze the astronomical spectra. The ISO SWS spectra of fifteen objects showing these PAH features were considered for this study, of which four have sufficient S/N between 5 and 6 $\mu$m to allow for an in-depth analysis. All four astronomical spectra show similar peak positions and profiles. The 5.25 $\mu$m feature is peaked and asymmetric, while the 5.7 $\mu$m feature is broader and flatter. Detailed analysis of the laboratory spectra and quantum chemical calculations show that the astronomical 5.25 and 5.7 $\mu$m bands are a blend of combination, difference and overtone bands primarily involving CH stretching and CH in-plane and CH out-of-plane bending fundamental vibrations. The experimental and computational spectra show that, of all the hydrogen adjacency classes possible on PAHs, solo and duo hydrogens consistently produce prominent bands at the observed positions whereas quartet hydrogens do not. In all, this a study supports the picture that astronomical PAHs are large with compact, regular structures. From the coupling with primarily strong CH out-of-plane bending modes one might surmise that the 5.25 and 5.7 $\mu$m bands track the neutral PAH population. However, theory suggests the role of charge in these astronomical bands might also be important.

To show the importance of high-spatial resolution observations of HII regions when compared with observations obtained with larger apertures such as ISO, we present mid-infrared spectra of two Magellanic Cloud HII regions, N88A and N160A. We obtained mid-infrared (8-13 um), long-slit spectra with TIMMI2 on the ESO 3.6m telescope. These are combined with archival spectra obtained with the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope, and are compared with the low-spatial resolution ISO-SWS data. An inventory of the spectra in terms of atomic fine-structure lines and molecular bands is presented. Concerning N88A, an isolated HII region with no adjacent infrared sources, the observations indicate that the line fluxes observed by ISO-SWS and Spitzer-IRS come exclusively from the compact HII region of about 3 arcsec in diameter. This is not the case for N160A, which has a more complex morphology. We have spectroscopically isolated for the first time the individual contributions of the three components of N160A, two high-excitation blobs, A1 and A2, and the young stellar object N160A-IR. In addition, extended [SIV] emission is observed with TIMMI2 and is most likely associated with the central star cluster located between A1 and A2. We show the value of these high-spatial resolution data in determining source characteristics, such as the degree of ionization of each high-excitation blob or the bolometric luminosity of the YSO. This luminosity is about one order of magnitude lower than previously estimated. For each high-excitation blob, we also determine the electron density and the elemental abundances of Ne, S, and Ar.

Infrared Dark Clouds (IRDCs) represent the earliest observed stages of clustered star formation, characterized by large column densities of cold and dense molecular material observed in silhouette against a bright background of mid-IR emission. Up to now, IRDCs were predominantly known toward the inner Galaxy where background infrared emission levels are high. We present Spitzer observations with the Infrared Camera Array toward object G111.80+0.58 (G111) in the outer Galactic Plane, located at a distance of ~3 kpc from us and ~10 kpc from the Galactic center. Earlier results show that G111 is a massive, cold molecular clump very similar to IRDCs. The mid-IR Spitzer observations unambiguously detect object G111 in absorption. We have identified for the first time an IRDC in the outer Galaxy, which confirms the suggestion that cluster-forming clumps are present throughout the Galactic Plane. However, against a low mid-IR back ground such as the outer Galaxy it takes some effort to find them.

Herbig Ae/Be stars are a class of young pre-main sequence stellar objects of intermediate mass and are known to have varying amounts of natal cloud material still present in their direct vicinity. We characterise the IR emission bands, due to fluorescence by PAH molecules, in the spectra of Herbig Ae/Be stars and link observed variations to spatial aspects of the mid-IR emission. We analysed two PAH dominated spectra from a sample of 15 Herbig Ae/Be stars observed with Spitzer and derive profiles of the major PAH bands. The shape and the measured band characteristics show pronounced variations between the two Spitzer spectra. Those variations parallel those found between three ISO spectra of other, well-studied, Herbig Ae/Be stars. The derived profiles are compared to those from a broad sample of sources. The Spitzer and ISO spectra exhibit characteristics commonly interpreted respectively as interstellar matter-like (ISM), non-ISM-like, or a combination of the two. We argue that the PAH emission detected from the sources exhibiting a combination of ISM-like and non-ISM-like characteristics indicates the presence of two dissimilar, spatially separated, PAH families. As the shape of the individual PAH band profiles reflects the composition of the PAH molecules involved, this demonstrates that PAHs in subsequent, evolutionary linked stages of star formation are different from those in the general ISM, implying active chemistry. None of the detected PAH emission can be associated with the (unresolved) disk and is thus associated with the circumstellar cloud. This implies that chemical changes may already occur in the natal cloud and not necessarily in the disk.

We use observations of the CI, CII, HI, and H_2 column densities along lines of sight in the Galactic plane to determine the formation rate of H_2 on grains and to determine chemical reaction rates with Polycyclic Aromatic Hydrocarbons. Photodissociation region models are used to find the best fit parameters to the observed columns. We find the H_2 formation rate on grains has a low rate (R ~ 1 x 10^(-17) cm^(3) s^(-1)) along lines of sight with low column density (A_V < 0.25) and low molecular fraction (f_(H_2) < 10^(-4)). At higher column densities (0.25 < A_V <2.13), we find a rate of R ~ 3.5x10^(-17) cm^(3) s^(-1). The lower rate at low column densities could be the result of grain processing by interstellar shocks which may deplete the grain surface area or process the sites of H +H formation, thereby inhibiting H_2 production. Alternatively, the formation rate may be normal, and the low molecular fraction may be the result of lines of sight which graze larger clouds. Such lines of sight would have a reduced H_2 self-shielding compared to the line-of-sight column. We find the reaction C^+ +PAH^- --> C + PAH^0 is best fit with a rate 2.4 x 10^(-7) \Phi_PAH T_2^(-0.5) cm^(3) s^(-1) with T_2= T/100 K and the reaction C^+ + PAH^0 --> C + PAH^+ is best fit with a rate 8.8x 10^(-9)\Phi_PAH cm^(3) s^(-1). In high column density gas we find \Phi_PAH ~ 0.4. In low column density gas, \Phi_PAH is less well constrained with \Phi_PAH ~ 0.2 - 0.4.

We present a mechanism for chondrules to stick together by means of compaction of a porous dust rim they sweep up as they move through the dusty nebula gas. It is shown that dust aggregates formed out of micron-sized grains stick to chondrules, forming a porous dust rim. When chondrules collide, this dust can be compacted by means of rolling motions within the porous dust layer. This mechanism dissipates the collisional energy, compacting the rim and allowing chondrules to stick. The structure of the obtained chondrule-dust agglomerates (referred to as compounds) then consists of three phases: chondrules, porous dust, and dust that has been compacted by collisions. Subsequently, these compounds accrete their own dust and collide with other compounds. The evolution of the compound size distribution and the relative importance of the phases is calculated by a Monte Carlo code. Growth ends, and a simulation is terminated when all the dust in the compounds has been compacted. Numerous runs are performed, reflecting the uncertainty in the physical conditions at the chondrule formation time. It is found that compounds can grow by 1-2 orders of magnitudes in radius, upto dm-sizes when turbulence levels are low. However, relative velocities associated with radial drift form a barrier for further growth. Earlier findings that the dust sweep-up by chondrules is proportional to their sizes are confirmed. We contrast two scenarios regarding how this dust evolved further towards the densely packed rims seen in chondrites.

We present the results of a systematic study of mid-IR spectra of Galactic regions, Magellanic HII regions, and galaxies of various types (dwarf, spiral, starburst), observed by the satellites ISO and Spitzer. We study the relative variations of the 6.2, 7.7, 8.6 and 11.3 micron features inside spatially resolved objects (such as M82, M51, 30 Doradus, M17 and the Orion Bar), as well as among 90 integrated spectra of 50 objects. Our main results are that the 6.2, 7.7 and 8.6 micron bands are essentially tied together, while the ratios between these bands and the 11.3 micron band varies by one order of magnitude. This implies that the properties of the PAHs are remarkably universal throughout our sample, and that the relative variations of the band ratios are mainly controled by the fraction of ionized PAHs. In particular, we show that we can rule out both the modification of the PAH size distribution, and the mid-infrared extinction, as an explanation of these variations. Using a few well-studied Galactic regions (including the spectral image of the Orion Bar), we give an empirical relation between the I(6.2)/I(11.3) ratio and the ionization/recombination ratio G0/ne.Tgas^0.5, therefore providing a useful quantitative diagnostic tool of the physical conditions in the regions where the PAH emission originates. Finally, we discuss the physical interpretation of the I(6.2)/I(11.3) ratio, on galactic size scales.

We have examined the relationship between the optical depth of the 9.7 micron silicate absorption feature (tau_9.7) and the near-infrared color excess, E(J-Ks) in the Serpens, Taurus, IC 5146, Chameleon I, Barnard 59, and Barnard 68 dense clouds/cores. Our data set, based largely on Spitzer IRS spectra, spans E(J-Ks)=0.3 to 10 mag (corresponding to visual extinction between about 2 and 60 mag.). All lines of sight show the 9.7 micron silicate feature. Unlike in the diffuse ISM where a tight linear correlation between the 9.7 micron silicate feature optical depth and the extinction (Av) is observed, we find that the silicate feature in dense clouds does not show a monotonic increase with extinction. Thus, in dense clouds, tau_9.7 is not a good measure of total dust column density. With few exceptions, the measured tau_9.7 values fall well below the diffuse ISM correlation line for E(J-Ks) > 2 mag (Av >12 mag). Grain growth via coagulation is a likely cause of this effect.

High resolution spectra of the Spitzer Space Telescope show vibration-rotation absorption bands of gaseous C2H2, HCN, and CO2 molecules toward a sample of deeply obscured (U)LIRG nuclei. The observed bands reveal the presence of dense (n>~ 10^7 cm^-3), warm (T_ex = 200-700 K) molecular gas with high column densities of these molecules ranging from a few 10^15 - 10^17 cm^-2. Abundances relative to H2, inferred from the silicate optical depth, range from ~10^-7 to 10^-6 and show no correlation with temperature. Theoretical studies show that the high abundances of both C2H2 and HCN exclude a X-ray dominated region (XDR) associated with the toroid surrounding an AGN as the origin of this dense warm molecular gas. Galactic massive protostars in the so-called Hot Core phase have similar physical characteristics with comparable high abundances of C2H2, HCN, and CO2 in the hot phase. However, the abundances of C2H2 and HCN and the C2H2/CO2 and HCN/CO2 ratios are much higher toward the (U)LIRGs in the cooler (T_ex <= 400 K) phase. We suggest that the warm dense molecular gas revealed by the mid-IR absorption lines is associated with a phase of deeply embedded star formation where the extreme pressures and densities of the nuclear starburst environment have inhibited the expansion of HII regions and the global disruption of the star forming molecular cloud cores, and `trapped' the star formation process in an `extended' Hot Core phase.

We present results of the analysis of cometary X-ray spectra with an extended version of our charge exchange emission model (Bodewits et al. 2006). We have applied this model to the sample of 8 comets thus far observed with the Chandra X-ray observatory and ACIS spectrometer in the 300-1000 eV range. The surveyed comets are C/1999 S4 (LINEAR), C/1999 T1 (McNaught-Hartley), C/2000 WM1 (LINEAR), 153P/2002 (Ikeya-Zhang), 2P/2003 (Encke), C/2001 Q4 (NEAT), 9P/2005 (Tempel 1) and 73P/2006-B (Schwassmann-Wachmann 3) and the observations include a broad variety of comets, solar wind environments and observational conditions. The interaction model is based on state selective, velocity dependent charge exchange cross sections and is used to explore how cometary X-ray emission depend on cometary, observational and solar wind characteristics. It is further demonstrated that cometary X-ray spectra mainly reflect the state of the local solar wind. The current sample of Chandra observations was fit using the constrains of the charge exchange model, and relative solar wind abundances were derived from the X-ray spectra. Our analysis showed that spectral differences can be ascribed to different solar wind states, as such identifying comets interacting with (I) fast, cold wind, (II), slow, warm wind and (III) disturbed, fast, hot winds associated with interplanetary coronal mass ejections. We furthermore predict the existence of a fourth spectral class, associated with the cool, fast high latitude wind.

Context: Existing SWAS observations and future HIFI/Herschel data require a clear sense of the information content of water emission and absorption lines. Aims: Investigate wether the ground-state transition of ortho-H2O, 1_(10)-1_(01), at 557GHz can be used to measure the column density throughout an interstellar cloud. Methods: We make use of a multi-zone escape probability code suitable for the treatment of molecular line emission. Results: For low abundances, i.e., X(H2O)=<10^(-9), the intensity of the 1_(10)-1_(01) transition scales with the total column of H2. However, this relationship breaks down with increasing abundance, i.e., optical depth, due to line trapping and -- for T_(dust)>=25K, X(H2O)=<10^(-8) and n~10^4 cm^(-3) -- absorption of the dust continuum. Conclusions: An observed decline in intensity per column, expected if H2O is a surface tracer, does not necessarily mean that the water is absent in the gas phase at large columns, but can be caused by line trapping and subsequent collisional de-excitation. To determine the amount of water vapour in the interstellar medium, multiple line measurements of optically thin transitions are needed to disentangle radiative transfer and local excitation effects.

Stars form in the densest, coldest, most quiescent regions of molecular clouds. Molecules provide the only probes which can reveal the dynamics, physics, chemistry and evolution of these regions, but our understanding of the molecular inventory of sources and how this is related to their physical state and evolution is rudimentary and incomplete. The Spectral Legacy Survey (SLS) is one of seven surveys recently approved by the JCMT Board. Starting in 2007, the SLS will produce a spectral imaging survey of the content and distribution of all the molecules detected in the 345 GHz atmospheric window (between 332 GHz and 373 GHz) towards a sample of 5 sources. Our intended targets are: a low mass core (NGC1333 IRAS4), 3 high mass cores spanning a range of star forming environments and evolutionary states (W49, AFGL2591, and IRAS20126), and a PDR (the Orion Bar). The SLS will use the unique spectral imaging capabilities of HARP-B/ACSIS to study the molecular inventory and the physical structure of these objects, which span different evolutionary stages and physical environments, to probe their evolution during the star formation process. As its name suggests, the SLS will provide a lasting data legacy from the JCMT that is intended to benefit the entire astronomical community. As such, the entire data set (including calibrated spectral datacubes, maps of molecular emission, line identifications, and calculations of the gas temperature and column density) will be publicly available.

The formation of planets is directly linked to the evolution of the circumstellar (CS) disk from which they are born. The dissipation timescales of CS disks are, therefore, of direct astrophysical importance in evaluating the time available for planet formation. We employ Spitzer Space Telescope spectra to complete the CS disk census for the late-type members of the ~8 Myr-old eta Chamaeleontis star cluster. Of the 15 K- and M-type members, eight show excess emission. We find that the presence of a CS disk is anti-correlated with binarity, with all but one disk associated with single stars. With nine single stars in total, about 80% retain a CS disk. Of the six known or suspected close binaries the only CS disk is associated with the primary of RECX 9. No circumbinary disks have been detected. We also find that stars with disks are slow rotators with surface values of specific angular momentum j = 2-15 j_sun. All high specific angular momentum systems with j = 20-30 j_sun are confined to the primary stars of binaries. This provides novel empirical evidence for rotational disk locking and again demonstrates the much shorter disk lifetimes in close binary systems compared to single star systems. We estimate the characteristic mean disk dissipation timescale to be ~5 Myr and ~9 Myr for the binary and single star systems, respectively.

Context: Sticking of colliding dust particles through van der Waals forces is the first stage in the grain growth process in protoplanetary disks, eventually leading to the formation of comets, asteroids and planets. A key aspect of the collisional evolution is the coupling between dust and gas motions, which depends on the internal structure (porosity) of aggregates. Aims: To quantify the importance of the internal structure on the collisional evolution of particles, and to create a new coagulation model to investigate the difference between porous and compact coagulation in the context of a turbulent protoplanetary disk. Methods: We have developed simple prescriptions for the collisional evolution of porosity of grain-aggregates in grain-grain collisions. Three regimes can then be distinguished: `hit-and-stick' at low velocities, with an increase in porosity; compaction at intermediate velocities, with a decrease of porosity; and fragmentation at high velocities. (..) Results: (..) We can discern three different stages in the particle growth process (..) We find that when compared to standard, compact models of coagulation, porous growth delays the onset of settling, because the surface area-to-mass ratio is higher, a consequence of the build-up of porosity during the initial stages. As a result, particles grow orders of magnitudes larger in mass before they rain-out to the mid-plane. Depending on the turbulent viscosity and on the position in the nebula, aggregates can grow to (porous) sizes of ~ 10 cm in a few thousand years. We also find that collisional energies are higher than in the limited PCA/CCA fractal models, thereby allowing aggregates to restructure. It is concluded that the microphysics of collisions plays a key role in the growth process.

Color-magnitude diagrams (CMDs) are presented for the Spitzer SAGE (Surveying the Agents of a Galaxy's Evolution) survey of the Large Magellanic Cloud (LMC). IRAC and MIPS 24 um epoch one data are presented. These data represent the deepest, widest mid-infrared CMDs of their kind ever produced in the LMC. Combined with the 2MASS survey, the diagrams are used to delineate the evolved stellar populations in the Large Magellanic Cloud as well as Galactic foreground and extragalactic background populations. Some 32000 evolved stars brighter than the tip of the red giant branch are identified. Of these, approximately 17500 are classified as oxygen-rich, 7000 carbon-rich, and another 1200 as ``extreme'' asymptotic giant branch (AGB) stars. Brighter members of the latter group have been called ``obscured'' AGB stars in the literature owing to their dusty circumstellar envelopes. A large number (1200) of luminous oxygen--rich AGB stars/M supergiants are also identified. Finally, there is strong evidence from the 24 um MIPS channel that previously unexplored, lower luminosity oxygen-rich AGB stars contribute significantly to the mass loss budget of the LMC (1200 such sources are identified).

We are performing a uniform and unbiased, ~7x7 degrees imaging survey of the Large Magellanic Cloud (LMC), using the IRAC and MIPS instruments on board the Spitzer Space Telescope in order to survey the agents of a galaxy's evolution (SAGE), the interstellar medium (ISM) and stars in the LMC. The detection of diffuse ISM with column densities >1.2x10^21 H cm^-2 permits detailed studies of dust processes in the ISM. SAGE's point source sensitivity enables a complete census of newly formed stars with masses >3 solar masses that will determine the current star formation rate in the LMC. SAGE's detection of evolved stars with mass loss rates >1x10^-8 solar masses per year will quantify the rate at which evolved stars inject mass into the ISM of the LMC. The observing strategy includes two epochs in 2005, separated by three months, that both mitigate instrumental artifacts and constrain source variability. The SAGE data are non-proprietary. The data processing includes IRAC and MIPS pipelines and a database for mining the point source catalogs, which will be released to the community in support of Spitzer proposal cycles 4 and 5. We present initial results on the epoch 1 data with a special focus on the N79 and N83 region. The SAGE epoch 1 point source catalog has ~4 million sources. The point source counts are highest for the IRAC 3.6 microns band and decrease dramatically towards longer wavelengths consistent with the fact that stars dominate the point source catalogs and that the dusty objects, e.g. young stellar objects and dusty evolved stars that detected at the longer wavelengths, are rare in comparison. We outline a strategy for identifying foreground MW stars, that may comprise as much as 18% of the source list, and background galaxies, that may comprise ~12% of the source list.

Context: Mass loss plays a dominant role in the evolution of low mass stars while they are on the Asymptotic Giant Branch (AGB). The gas and dust ejected during this phase are a major source in the mass budget of the interstellar medium. Recent studies have pointed towards the importance of variations in the mass-loss history of such objects. Aims: By modelling the full line profile of low excitation CO lines emitted in the circumstellar envelope, we can study the mass-loss history of AGB stars. Methods: We have developed a non-LTE radiative transfer code, which calculates the velocity structure and gas kinetic temperature of the envelope in a self-consistent way. The resulting structure of the envelope provides the input for the molecular line radiative calculations which are evaluated in the comoving frame. The code allows for the implementation of modulations in the mass-loss rate. This code has been benchmarked against other radiative transfer codes and is shown to perform well and efficiently. Results: We illustrate the effects of varying mass-loss rates in case of a superwind phase. The model is applied to the well-studied case of VY CMa. We show that both the observed integrated line strengths as the spectral structure present in the observed line profiles, unambiguously demonstrate that this source underwent a phase of high mass loss (~ 3.2E-4 Msun/yr) some 1000 yr ago. This phase took place for some 100 yr, and was preceded by a low mass-loss phase (~ 1E-6 Msun/yr) taking some 800 yr. The current mass-loss rate is estimated to be in the order of 8E-5 Msun/yr. Conclusions: In this paper, we demonstrate that both the relative strength of the CO rotational line profiles and the (non)-occurrence of spectral structure in the profile offer strong diagnostics to pinpoint the mass-loss history.

We present late-time optical and mid-infrared observations of the Type-II supernova 2003gd in NGC 628. Mid-infrared excesses consistent with cooling dust in the ejecta are observed 499-678 days after outburst, and are accompanied by increasing optical extinction and growing asymmetries in the emission-line profiles. Radiative-transfer models show that up to 0.02 solar masses of dust has formed within the ejecta, beginning as early as 250 days after outburst. These observations show that dust formation in supernova ejecta can be efficient and that massive-star supernovae can be major dust producers throughout the history of the Universe.

In order to show the importance of high spatial resolution observations of extra-galactic sources when compared to observations obtained with larger apertures such as ISO, we present N-band spectra (8-13 um) of some locations in three starburst galaxies. In particular, the two galactic nuclei of the spiral galaxy NGC3256, the compact IR supernebula in the dwarf galaxy IIZw40 and the two brightest IR knots in the central starburst of the WR galaxy He2-10. The spectra have been obtained with TIMMI2 on the ESO 3.6m telescope. An inventory of the spectra in terms of atomic fine-structure lines and molecular bands is presented. We show the great value of these high spatial resolution data at constraining properties such as the extinction in the mid-IR, metallicity or stellar content (age, IMF, etc.). Regarding this, we have constrained the stellar content of the IR compact knot in IIZw40 by using the mid-IR fine-structure lines and setting restrictions on the nebular geometry. Considering the PAH bands, we have constructed a new mid-/far-IR diagnostic diagram based on the 11.2 um PAH and continuum, accessible to ground-based observations. Finally, we find a dependence between the presence or non-presence of PAHs and the hardness of the radiation field as measured by the [SIV]/[NeII] ratio. In particular, sources with PAH emission have in general a [SIV]/[NeII] ratio < 0.35. We investigate possible origins for this relation and conclude that it does not necessarily imply PAH destruction, but could also be explained by the PAH-dust competition for FUV photons.

Context : Despite the low cosmic abundance of deuterium (D/H ~ 1e-5), large degrees of deuterium fractionation in molecules are observed in star forming regions with enhancements that can reach 13 orders of magnitude, which current models have difficulties to account for. Aims : Multi-isotopologue observations are a very powerful constraint for chemical models. The aim of our observations is to understand the processes forming the observed large abundances of methanol and formaldehyde in low-mass protostellar envelopes (gas-phase processes ? chemistry on the grain surfaces ?) and better constrain the chemical models. Methods : Using the IRAM 30m single-dish telescope, we observed deuterated formaldehyde (HDCO and D2CO) and methanol (CH2DOH, CH3OD, and CHD2OH) towards a sample of seven low-mass class 0 protostars. Using population diagrams, we then derive the fractionation ratios of these species (abundance ratio between the deuterated molecule and its main isotopologue) and compare them to the predictions of grain chemistry models. Results : These protostars show a similar level of deuteration as in IRAS16293-2422, where doubly-deuterated methanol -- and even triply-deuterated methanol -- were first detected. Our observations point to the formation of methanol on the grain surfaces, while formaldehyde formation cannot be fully pined down. While none of the scenarii can be excluded (gas-phase or grain chemistry formation), they both seem to require abstraction reactions to reproduce the observed fractionations.

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.

We have analyzed the 9.7 and ``18'' micron interstellar silicate absorption features along the line of sight toward four heavily extincted galactic WC-type Wolf-Rayet (WR) stars. We construct two interstellar extinction curves from 1.25 to 25 micron using near-IR extinction measurements from the literature along with the silicate profiles of WR 98a (representing the local ISM) and GCS 3 (representing the Galactic Center). We have investigated the mineralogy of the interstellar silicates by comparing extinction profiles for amorphous silicates with olivine and pyroxene stochiometry to the 9.7 and ``18'' micron absorption features in the WR 98a spectrum. In this analysis, we have considered solid and porous spheres and a continuous distribution of ellipsoids. While it is not possible to simultaneously provide a perfect match to both profiles, we find the best match requires a mixture of these two types of compounds. We also consider iron oxides, aluminosilicates and silicate carbide (SiC) as grain components. Iron oxides cannot be accommodated in the observed spectrum, while the amount of Si in SiC is limited to <4%. Finally, we discuss the cosmic elemental abundance constraints on the silicate mineralogy, grain shape and porosity.

Far infrared fine structure line data from the ISO archive have been extracted for several hundred YSOs and their outflows, including molecular (CO) outflows, optical jets and Herbig-Haro (HH) objects. Given the importance of these lines to astrophysics, their excitation and transfer ought to be investigated in detail and, at this stage, the reliability of the diagnostic power of the fine structure transitions of O I and C II has been examined. Several issues, such as the extremely small intensity ratios of the oxygen 63 micron to 145 micron lines, are still awaiting an explanation. It is demonstrated that, in interstellar cloud conditions, the 145 micron line is prone to masing, but that this effect is likely an insufficient cause of the line ratio anomaly observed from cold dark clouds. Very optically thick emission could in principle also account for this, but would need similar, prohibitively high column densities and must therefore be abondoned as a viable explanation. One is left with [O I] 63 micron self absorption by cold and tenuous foreground gas, as has been advocated for distant luminous sources. Recent observations with the submillimeter observatory Odin support this scenario also in the case of nearby dark molecular clouds. On the basis of this large statistical material we are led to conclude that in star forming regions, the [O I] and [C II] lines generally have only limited diagnostic value.

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.

We present observations of methanol lines in a sample of Class 0 low mass protostars. Using a 1-D radiative transfer model, we derive the abundances in the envelopes. In two sources of the sample, the observations can only be reproduced by the model if the methanol abundance is enhanced by about two order of magnitude in the inner hot region of the envelope. Two other sources show similar jumps, although at a lower confidence level. The observations for the other three sources are well reproduced with a constant abundance, but the presence of a jump cannot be ruled out. The observed methanol abundances in the warm gas around low mass protostars are orders of magnitude higher than gas phase chemistry models predict. Hence, in agreement with other evidences, this suggest that the high methanol abundance reflects recent evaporation of ices due to the heating by the newly formed star. The observed abundance ratios of CH3 OH, H2 CO, and CO are in good agreement with grain surface chemistry models. However, the absolute abundances are more difficult to reproduce and may point towards the presence of multiple ice components in these regions.

Dust emission from the Type II supernova SN 2002hh in NGC 6946 has been detected at mid-infrared wavelengths by the Spitzer Space Telescope from 590 to 758 days after outburst and confirmed by higher angular resolution Gemini-N mid-IR observations. The day-600 5.8-24 micron emission can be fit by a 290-K blackbody having a luminosity of 1.6E+7 L_Sol. The minimum emitting radius of 1.1E+17 cm is too large for the emitting dust to have been formed in the supernova ejecta. Using radiative transfer models and realistic dust grain parameters, fits to the observed flux distribution could be obtained with an optically thick dust shell having a mass of 0.10-0.15 M_Sol, corresponding to a total dust+gas mass in excess of 10 M_Sol, suggesting a massive M supergiant or luminous blue variable precursor to this self-obscured object.

We report the spectroscopic detection of silicate emission at 10 and 18 microns in five PG quasars, the first detection of these two features in galaxies outside the Local Group. This finding is consistent with the unification model for Active Galactic Nuclei (AGNs), which predicts that an AGN torus seen pole-on should show a silicate emission feature in the mid-infrared. The strengths of the detected silicate emission features range from 0.12 to 1.25 times the continuum at 10 microns and from 0.20 to 0.79 times the continuum at 18 microns. The silicate grain temperatures inferred from the ratio of 18-to-10 micron silicate features under the assumption of optically thin emission range from 140 to 220K.

An overview is given of ISO results on regions of high excitation ISM and gas, i.e. HII regions, the Galactic Centre and Supernovae Remnants. IR emission due to fine-structure lines, molecular hydrogen, silicates, polycyclic aromatic hydrocarbons and dust are summarised, their diagnostic capabilities illustrated and their implications highlighted.

We reconsider H2 formation on grain surfaces. We develop a rate equation model which takes into account the presence of both physisorbed and chemisorbed sites on the surface, including quantum mechanical tunnelling and thermal diffusion. In this study, we took into consideration the uncertainties on the characteristics of graphitic surfaces. We calculate the H2 formation efficiency with the Langmuir Hinshelwood and Eley Rideal mechanisms, and discuss the importance of these mechanisms for a wide range of grain and gas temperatures. We also develop a Monte Carlo simulation to calculate the H2 formation efficiency and compare the results to our rate equation models. Our results are the following: (1) Depending on the barrier against chemisorption, we predict the efficiency of H2 formation for a wide range of grain and gas temperatures. (2) The Eley-Rideal mechanism has an impact on the H2 formation efficiency at high grain and gas temperatures. (3) The fact that we consider chemisorption in our model makes the rate equation and Monte Carlo approaches equivalent.

The results from a single-dish molecular line survey of a set of 18 deeply embedded young stellar objects are summarized. More than 40 lines from 16 different species were observed with the JCMT, Onsala, IRAM 30m and SEST telescopes. The multi-transition data are analyzed using a temperature and density structure derived from models of the dust continuum emission. For the outer envelope (>300 AU), the data indicate a `drop' abundance profile for many species, with normal abundances in the outer- and innermost regions and highly depleted abundances in an intermediate zone. This zone is bounded at the outer than the lifetime of the core, and at the inner edge by the evaporation temperature of the species involved. In the innermost envelope (<300 AU), all ices evaporate resulting in jumps in the abundances of complex organic molecules such as CH3OH. A key project for Herschel will be to survey gas-phase water in these objects, whose abundance shows extreme variations with temperature. ALMA wil be able to directly image the chemical variations throughout the envelope and zoom in on the inner hot core and protoplanetary disks on scales on tens of AU.

We present IRAM 30m and JCMT observations of HDO lines towards the solar-type protostar IRAS 16293-2422. Five HDO transitions have been detected on-source, and two were unfruitfully searched for towards a bright spot of the outflow of IRAS 16293-2422. We interpret the data by means of the Ceccarelli, Hollenbach and Tielens (1996) model, and derive the HDO abundance in the warm inner and cold outer parts of the envelope. The emission is well explained by a jump model, with an inner abundance of 1e-7 and an outer abundance lower than 1e-9 (3 sigma). This result is in favor of HDO enhancement due to ice evaporation from the grains in theinner envelope. The deuteration ratio HDO/H2O is found to be f_in=3% and f_out < 0.2% (3 sigma) in the inner and outer envelope respectively and therefore, the fractionation also undergoes a jump in the inner part of the envelope. These results are consistent with the formation of water in the gas phase during the cold prestellar core phase and storage of the molecules on the grains, but do not explain why observations of H2O ices consistently derive a H2O ice abundance of several 1e-5 to 1e-4, some two orders of magnitude larger than the gas phase abundance of water in the hot core around IRAS 16293-2422.

A solid-state feature was detected at around 2175 cm-1 towards 30 embedded young stellar objects in spectra obtained using the ESO VLT-ISAAC. We present results from laboratory studies of CO adsorbed at the surface of Zeolite wafers, where absorption bands were detected at 2177 and 2168 cm-1 (corresponding to CO chemisorbed at the Zeolite surface), and 2130 cm-1 (corresponding to CO physisorbed at the Zeolite surface), providing an excellent match to the observational data. We propose that the main carrier of the 2175-band is CO chemisorbed at bare surfaces of dust grains in the interstellar medium. This result provides the first direct evidence that gas-surface interactions do not have to result in the formation of ice mantles on interstellar dust. The strength of the 2175-band is estimated to be ~ 4 x 10-19 cm molecule-1. The abundance of CO adsorbed at bare grain surfaces ranges from 0.06 to 0.16 relative to H2O ice, which is, at most, half of the abundance (relative to H2O ice) of CO residing in H2O-dominated ice environments. These findings imply that interstellar grains have a large (catalytically-active) surface area, providing a refuge for interstellar species. Consequently the potential exists for heterogeneous chemistry to occur involving CO molecules in unique surface chemistry pathways not currently considered in gas-grain models of the interstellar medium.

Complex organic molecules have previously been discovered in solar type protostars, raising the questions of where and how they form in the envelope. Possible formation mechanisms include grain mantle evaporation, interaction of the outflow with its surroundings or the impact of UV/X-rays inside the cavities. In this Letter we present the first interferometric observations of two complex molecules, CH3CN and HCOOCH3, towards the solar type protostar IRAS16293-2422. The images show that the emission originates from two compact regions centered on the two components of the binary system. We discuss how these results favor the grain mantle evaporation scenario and we investigate the implications of these observations for the chemical composition and physical and dynamical state of the two components.

We report the detection of complex molecules (HCOOCH_3, HCOOH and CH_3CN), signposts of a "hot core" like region, toward the low mass, Class 0 source NGC1333-IRAS4A. This is the second low mass protostar where such complex molecules have been searched for and reported, the other source being IRAS16293-2422. It is therefore likely that compact (few tens of AUs) regions of dense and warm gas, where the chemistry is dominated by the evaporation of grain mantles, and where complex molecules are found, are common in low mass Class 0 sources.Given that the chemical formation timescale is much shorter than the gas hot core crossing time, it is not clear whether the reported complex molecules are formed on the grain surfaces (first generation molecules) or in the warm gas by reactions involving the evaporated mantle constituents (second generation molecules). We do not find evidence for large differences in the molecular abundances, normalized to the formaldehyde abundance, between the two solar type protostars, suggesting perhaps a common origin.

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.

IR emission bands at 3.3, 6.2, 7.7, 8.6 and 11.3 um are generally attributed to IR fluorescence from (mainly) FUV pumped PAHs. As such, they trace the FUV stellar flux and are a measure of star formation. We examined the IR spectral characteristics of Galactic star forming regions, normal and starburst galaxies, AGNs and ULIRGs. The goal is to analyze if PAH bands are a good qualitative and/or quantitative tracer of star formation and hence the application of PAH bands as a diagnostic in order to identify the dominant processes contributing to the IR emission from Seyfert's and ULIRGs. We develop a MIR/FIR diagnostic and compare it to known diagnostics, with these also applied to the Galactic sample. This diagnostic is based on the FIR normalized 6.2 um PAH flux and the FIR normalized 6.2 um continuum flux. The Galactic sources form a sequence spanning a range of 3 orders of magnitude, from embedded compact HII regions to exposed PDRs and the (D)ISM. The variation in the 6.2 um PAH/continuum ratio is relative small. Normal and starburst galaxies ressemble exposed PDRs. While Seyfert-2's coincide with the starburst trend, Seyfert-1's are displaced by at least a factor 10 in 6.2 um continuum flux. ULIRGs show a diverse spectral appearance (AGN hot dust continuum, starburst-like or strong dust obscuration in the nucleus). ULIRGs also seems to have more prominent FIR emission than either starburst galaxies or AGNs. We discuss the observed variation in the Galactic sample in view of the evolutionary state and the PAH/dust abundance and the use of PAHs as quantitative tracers of star formation activity. We find that PAHs may be better suited as a tracer of B stars, which dominate the Galactic stellar energy budget, than as a tracer of massive star formation (O stars).

We present spectra of the 3.3 $\mu$m and 11.2 $\mu$m PAH features of a large number of (extra-) galactic sources, obtained with ISO-SWS. Clear variations are present in the profiles of these features. The sources are classified independently based on the 3.3 and 11.2 $\mu$m feature profiles and peak positions. Correlations between these classes and those based on the 6--9 $\mu$m features (Peeters et al. 2002) are found. Also, these classifications depend on the type of object. The observed pronounced contrast in the spectral variations for the CH modes (3.3 and 11.2 $\mu$m bands) versus the CC modes (6.2, 7.7 and 8.6 $\mu$m bands) is striking : the peak wavelengths of the features attributed to CC modes vary by $\sim$15--80 cm$^{-1}$, while for the CH modes the variations are $\sim$4--6.5 cm$^{-1}$. We summarize existing laboratory data and theoretical calculations of PAH molecules and complexes. In contrast to the 6.2 and 7.7 $\mu$m components which are attributed to PAH cations, the 3.3 $\mu$m feature appears to originate in neutral and/or negatively charged PAHs. We attribute the variations in peak position and profile of these features to the composition of the PAH family. The variations in FWHM of the 3.3 $\mu$m feature remains an enigma while those of the 11.2 $\mu$m can be explained by anharmonicity and molecular structure. The possible origin of the observed contrast in profile variations between the CH modes and the CC modes is highlighted.

We have studied the dust along the line-of-sight towards the Galactic Center using Short Wavelength Spectrometer (SWS) data obtained with the Infrared Space Observatory (ISO). We focussed on the wavelength region from 8-13 micron which is dominated by the strong silicate absorption feature. Using the absorption profiles observed towards Galactic Center Sources (GCS) 3 and 4, which are C-rich Wolf-Rayet Stars, as reference objects, we are able to disentangle the interstellar silicate absorption and the silicate emission intrinsic to the source, toward Sgr A* and derive a very accurate profile for the intrinsic 9.7 micron band. The interstellar absorption band is smooth and featureless and is well reproduced using a mixture of 15.1% amorphous pyroxene and 84.9% of amorphous olivine by mass, all in spherical sub-micron-sized grains. There is no direct evidence for substructure due to interstellar crystalline silicates. We are able to determine an upper limit to the degree of crystallinity of silicates in the diffuse interstellar medium (ISM), and conclude that the crystalline fraction of the interstellar silicates is 0.2% (+/- 0.2%) by mass. This is much lower than the degree of crystallinity observed in silicates in the circumstellar environment of evolved stars, the main contributors of dust to the ISM. There are two possible explanations for this discrepancy. First, an amorphization process occurs in the ISM on a time scale significantly shorter than the destruction time scale, possibly caused by particle bombardment by heavyweight ions. Second, we consider the possibility that the crystalline silicates in stellar ejecta are diluted by an additional source of amorphous silicates, in particular supernovae.

We present a systematic study of the effect of metallicity on the stellar spectral energy distribution (SED) of O main sequence (dwarf) stars, focussing on the hydrogen and helium ionizing continua, and on the optical and near-IR lines used for spectral classification. The spectra are based on non-LTE line blanketed atmosphere models with stellar winds calculated using the CMFGEN code of Hillier & Miller (1998). We draw the following conclusions. First, we find that the total number of Lyman photons emitted is almost independent of line blanketing effects and metallicity for a given effective temperature. This is because the flux that is blocked by the forest of metal lines at wavelengths shorter than 600 Angstrom is redistributed mainly within the Lyman continuum. Second, the spectral type, as defined by the ratio of the equivalent widths of HeI 4471 Angstrom and HeII 4542 Angstrom, is shown to depend noticeably on the microturbulent velocity in the atmosphere, on metallicity and, within the luminosity class of dwarfs, on gravity. Third, we confirm the decrease in the effective temperature for a given spectral type due to the inclusion of line blanketing recently found by e.g. Martins et al. (2002). Finally, we find that the SED below ~450 Angstrom is highly dependent on metallicity. This is reflected in the behaviour of nebular fine-structure line ratios such as [NeIII]/[NeII] 15.5/12.8 and [ArIII]/[ArII] 9.0/7.0 micron. This dependence complicates the use of these nebular ratios as diagnostic tools for the effective temperature determination of the ionizing stars in HII regions and for age dating of starburst regions in galaxies.

Infrared spectroscopy is a powerful tool to probe the inventory of solid state and molecular species in circumstellar ejecta. Here we analyse the infrared spectrum of the Asymptotic Giant Branch star W Hya, obtained by the Short and Long Wavelength Spectrometers on board of the Infrared Satellite Observatory. These spectra show evidence for the presence of amorphous silicates, aluminum oxide, and magnesium-iron oxide grains. We have modelled the spectral energy distribution using laboratory measured optical properties of these compounds and derive a total dust mass loss rate of 3E-10 Msol/yr. We find no satisfactory fit to the 13 micron dust emission feature and the identification of its carrier is still an open issue. We have also modelled the molecular absorption bands due to H2O, OH, CO, CO2, SiO, and SO2 and estimated the excitation temperatures for different bands which range from 300 to 3000K. It is clear that different molecules giving rise to these absorption bands originate from different gas layers. We present and analyse high resolution Fabry-Perot spectra of the three CO2 bands in the 15 micron region. In these data, the bands are resolved into individual Q-lines in emission, which allows the direct determination of the excitation temperature and column density of the emitting gas. This reveals the presence of a warm (about 450K) extended layer of CO2, somewhere between the photosphere and the dust formation zone. The gas in this layer is cooler than the 1000K CO2 gas responsible for the low-resolution absorption bands at 4.25 and 15 micron. The rotational and vibrational excitation temperatures derived from the individual Q-branch lines of CO2 are different (450K and 150K, respectively) so that the CO2 level population is not in LTE.

The Infrared Space Observatory (ISO) has provided the first complete mid-IR spectra for a wide range of objects. Almost all of these spectra are dominated by the well-known infrared emission features at 3.3, 6.2, 7.7, and 11.2 $\mu$m, the so-called Unidentified Infra-Red (UIR) features. Besides the major features, there is an array of minor features and broad plateaux stretching from 3 to 20 $\mu$m which reveal subtle details of conditions in the emission zones and properties of the carriers. Generally attributed to the vibrational relaxation of UV-pumped Polycyclic Aromatic Hydrocarbon molecules (PAHs) containing some 50--100 C-atoms, these UIR spectra are a treasure trove of information. The ISO spectra have, for the first time, allowed a systematic analysis of the spectral characteristics of the UIR features in a wide variety of environments. The peak positions, profiles, and relative strengths of the major features vary from source to source and spatially within sources. These specific profiles are not unique to certain object types but can occur within each individual source. Here, we review ISO and recent ground-based observations and assess some of their implications.

We present a survey of the formaldehyde emission in a sample of eight Class 0 protostars obtained with the IRAM and JCMT telescopes. The data have been analyzed with three different methods with increasing level of sophistication. We first analyze the observed emission in the LTE approximation, and derive rotational temperatures between 11 and 40 K, and column densities between 1 and 20 x 10^13 cm^-2. Second, we use a LVG code and derive larger kinetic temperatures, between 30 and 90 K, consistent with subthermally populated levels and densities from 1 to 6 x 10^5 cm^-3. The column densities from the LVG modeling are within a factor of 10 with respect to those derived in the LTE approximation. Finally, we analyze the observations based upon detailed models for the envelopes surrounding the protostars, using temperature and density profiles previously derived from continuum observations. We approximate the formaldehyde abundance across the envelope with a jump function, the jump occurring when the dust temperature reaches 100 K, the evaporation temperature of the grain mantles. The observed formaldehyde emission is well reproduced only if there is a jump, more than two orders of magnitude, in four sources. In the remaining four sources the data are consistent with a formaldehyde abundance jump, but the evidence is more marginal (~2 sigma). The inferred inner H2CO abundance varies between 1 x 10^-8 and 6 x 10^-6. We discuss the implications of these jumps for our understanding of the origin and evolution of ices in low mass star forming regions. Finally, we give predictions for the submillimeter H2CO lines, which are particularly sensitive to the abundance jumps.

We report the first detection of triply-deuterated methanol, with 12 observed transitions, towards the low-mass protostar IRAS 16293-2422, as well as multifrequency observations of 13CH3OH, used to derive the column density of the main isotopomer CH3OH. The derived fractionation ratio [CD3OH]/[CH3OH] averaged on a 10'' beam is 1.4%. Together with previous CH2DOH and CHD2OH observations, the present CD3OH observations are consistent with a formation of methanol on grain surfaces, if the atomic D/H ratio is 0.1 to 0.3 in the accreting gas. Such a high atomic ratio can be reached in the frame of gas-phase chemical models including all deuterated isotopomers of H3+.

We have re-analyzed the 6-12 micron ISO spectrum of the ultra-luminous infrared galaxy Arp220 with the conclusion that it is not consistent with that of a scaled up version of a typical starburst. Instead, both template fitting with spectra of the galaxies NGC4418 and M83 and with dust models suggest that it is best represented by combinations of a typical starburst component, exhibiting PAH emission features, and a heavily absorbed dust continuum which contributes ~40% of the 6-12 micron flux and likely dominates the luminosity. Of particular significance relative to previous studies of Arp220 is the fact that the emission feature at 7.7 micron comprises both PAH emission and a broader component resulting from ice and silicate absorption against a heavily absorbed continuum. Extinction to the PAH emitting source, however, appears to be relatively low. We tentatively associate the PAH emitting and heavily dust/ice absorbed components with the diffuse emission region and the two compact nuclei respectively identified by Soifer et al. (2002) in their higher spatial resolution 10 micron study. Both the similarity of the absorbed continuum with that of the embedded Galactic protostars and results of the dust models imply that the embedded source(s) in Arp220 could be powered by, albeit extremely dense, starburst activity. Due to the high extinction, it is not possible with the available data to exclude that AGN(s) also contribute some or all of the observed luminosity. In this case, however, the upper limit measured for its hard X-ray emission would require Arp220 to be the most highly obscured AGN known.

We present ground-based 2.1 to 4.2 microns observations of four low-mass protostars. We searched for the 4.1 microns OD stretch band, characteristic of solid HDO in grain mantles. We did not detect solid HDO in any of the four sources, but we derive 3-sigma upper limits from 0.5% to 2% for the HDO/H2O ratio depending on the source. These ratios provide strong constraints to solid-state deuteration models when compared to deuterium fractionation values observed in the gas phase. We discuss various scenarios that could lead to such a low water deuteration compared to the high formaldehyde and methanol deuteration observed in the gas phase.

Sun like stars are born from the collapse of fragment of molecular clouds. During the first embedded phase (the so called class 0/1), the protostar is surrounded by a collapsing envelope, whose physical and chemical structure would set up the initial conditions of the proto-stellar disk which may eventually form planets. In this contribution, we show how spectroscopic far infrared to millimeter observations can be a powerful diagnostic to derive the chemical and physical structure of the innermost regions of envelopes collapsing around solar type protostars, where disk are believe to form. Our claim is based on the simultaneous modeling of the dynamics, chemistry and radiative transfer in a protostellar envelope. Here we present a study of a class 0 protostar NGC1333-IRAS4, obtained using ISO-LWS observations of water far-infrared lines. We show the derived physical and chemical structure of the protostar, namely the density and temperature profiles, as well as the water abundance across the envelope. Besides, the mass of the central protostar and its accretion rate can be constrained by these observations and modeling, yielding an alternative method to measure these two poorly know parameters.

Sun like stars are born from the collapse of fragment of molecular clouds. During the first embedded phase (the so called class 0/1), the protostar is surrounded by a collapsing envelope, whose physical and chemical structure would set up the initial conditions of the proto-stellar disk which may eventually form planets. In this contribution, we show how spectroscopic far infrared to millimeter observations can be a powerful diagnostic to derive the chemical and physical structure of the innermost regions of envelopes collapsing around solar type protostars, where disk are believe to form. Our claim is based on the simultaneous modeling of the dynamics, chemistry and radiative transfer in a protostellar envelope. Here we present a study of a class 0 protostar NGC1333-IRAS4, obtained using ISO-LWS observations of water far-infrared lines. We show the derived physical and chemical structure of the protostar, namely the density and temperature profiles, as well as the water abundance across the envelope. Besides, the mass of the central protostar and its accretion rate can be constrained by these observations and modeling, yielding an alternative method to measure these two poorly know parameters.

The elemental abundances of ten planetary nebulae, derived with high accuracy including ISO and IUE spectra, are analysed with the aid of synthetic evolutionary models for the TP-AGB phase. Model prescriptions are varied until we achieve the simultaneous reproduction of all elemental features, which allows placing important constraints on the characteristic masses and nucleosynthetic processes experienced by the stellar progenitors. First of all, it is possible to separate the sample into two groups of PNe, one indicating the occurrence of only the third dredge-up during the TP-AGB phase, and the other showing also the chemical signature of hot-bottom burning. The former group is reproduced by stellar models with variable molecular opacities (see Marigo 2002), adopting initial solar metallicity, and typical efficiency of the third dredge-up 0.3-0.4. The latter group of PNe, with extremely high He content 0.15<=He/H<=0.20 and marked oxygen deficiency, is consistent with original sub-solar metallicity (i.e. LMC composition). Moreover, we are able to explain quantitatively both the N/H-He/H correlation and the N/H-C/H anti-correlation, thus solving the discrepancy pointed out long ago by Becker & Iben (1980). This is obtained only under the hypothesis that intermediate-mass TP-AGB progenitors (M > 4.5-5.0 Mo) with LMC composition have suffered a number of very efficient, carbon-poor, dredge-up events. Finally, the neon abundances of the He-rich PNe can be recovered by invoking a significant production of 22Ne during thermal pulses, which would imply a reduced role of the 22Ne(alpha,n)25Mg reaction as neutron source to the s-process nucleosynthesis in these stars.

Medium resolution (lambda/Delta lambda = 5000-10000) VLT-ISAAC M-band spectra are presented of 39 young stellar objects in nearby low-mass star forming clouds showing the 4.67 micron stretching vibration mode of solid CO. By taking advantage of the unprecedentedly large sample, high S/N ratio and high spectral resolution, similarities in the ice profiles from source to source are identified. It is found that excellent fits to all the spectra can be obtained using a phenomenological decomposition of the CO stretching vibration profile at 4.67 micron into 3 components. All observed interstellar CO profiles can thus be uniquely described by a model depending on only 3 linear fit parameters, indicating that a maximum of 3 specific molecular environments of solid CO exist under astrophysical conditions. A simple physical model of the CO ice is presented, which shows that the 2139.9 cm-1 component is indistinguishable from pure CO ice. It is concluded, that in the majority of the observed lines of sight, 60-90% of the CO is in a nearly pure form. In the same model the 2143.7 cm-1 component can possibly be explained by the longitudinal optical (LO) component of the vibrational transition in pure crystalline CO ice which appears when the background source is linearly polarised. Furthermore the absorption band from solid 13CO at 2092 cm-1 is detected towards IRS 51 in the Rho Ophiuchi cloud complex and an isotopic ratio of 12CO/13CO=68+/-10 is derived. It is shown that all the observed solid 12CO profiles, along with the solid 13CO profile, are consistent with grains with an irregularly shaped CO ice mantle simulated by a Continuous Distribution of Ellipsoids (CDE), but inconsistent with the commonly used models of spherical grains in the Rayleigh limit. [Abstract abridged]

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.

A high quality, medium-resolution K-band spectrum has been obtained of the ultracompact HII region G29.96-0.02 with the Very Large Telescope (VLT). The slit was positioned along the symmetry axis of the cometary shaped nebula. Besides the spectrum of the embedded ionizing O star, the long-slit observation revealed the rich emission-line spectrum produced by the ionized nebula with sub-arcsec spatial resolution. The nebular spectrum includes Br$\gamma$, several helium emission lines and a molecular hydrogen line. A detailed analysis is presented of the variation in strength, velocity and width of the nebular emission lines along the slit. The results are consistent with previous observations, but the much better spatial resolution allows a critical evaluation of models explaining the cometary shape of the nebula. Our observations support neither the wind bow shock model nor the champagne flow model. The measured line ratios of the nebular hydrogen and helium lines are compared to predictions from case B recombination-line theory. The results indicate an electron temperature between 6400 and 7500 K, in good agreement with other determinations and the Galactocentric distance of 4.6 kpc. The He$^{+}$/H$^{+}$ ratio is practically constant over the slit; we argue that He is singly ionized throughout the nebula. We review the various observational constraints on the effective temperature of the ionizing star and show that these are in agreement with its K-band spectral type of O5-O6V.

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 a detailed study of the circumstellar envelope of the post-asymptotic giant branch ``21 micron object'' HD 56126. We build a detailed dust radiative transfer model of the circumstellar envelope in order to derive the dust composition and mass, and the mass-loss history of the star. To model the emission of the dust we use amorphous carbon, hydrogenated amorphous carbon, magnesium sulfide and titanium carbide. We present a detailed parametrisation of the optical properties of hydrogenated amorphous carbon as a function of H/C content. The mid-infrared imaging and spectroscopy is best reproduced by a single dust shell from 1.2 to 2.6 arcsec radius around the central star. This shell originates from a short period during which the mass-loss rate exceeded 10^(-4) M_sun/yr. We find that the strength of the ``21'' micron feature poses a problem for the TiC identification. The low abundance of Ti requires very high absorption cross-sections in the ultraviolet and visible wavelength range to explain the strength of the feature. Other nano-crystalline metal carbides should be considered as well. We find that hydrogenated amorphous carbon in radiative equilibrium with the local radiation field does not reach a high enough temperature to explain the strength of the 3.3-3.4 and 6-9 micron hydrocarbon features relative to the 11-17 micron hydrocarbon features. We propose that the carriers of these hydrocarbon features are not in radiative equilibrium but are transiently heated to high temperature. We find that 2 per cent of the dust mass is required to explain the strength of the ``30'' micron feature, which fits well within the measured atmospheric abundance of Mg and S. This further strengthens the MgS identification of the ``30'' micron feature.

Much of the interstellar medium in disk galaxies is in the form of neutral atomic hydrogen, H I. This gas can be in thermal equilibrium at relatively low temperatures, T < 300 K (the cold neutral medium, or CNM) or at temperatures somewhat less than 10^4 K (the warm neutral medium, or WNM). These two phases can coexist over a narrow range of pressures, P_min < P < P_max. We determine P_min and P_max in the plane of the Galaxy as a function of Galactocentric radius R using recent determinations of the gas heating rate and the gas phase abundances of interstellar gas. We provide an analytic approximation for P_min as a function of metallicity, far-ultraviolet radiation field, and the ionization rate of atomic hydrogen. Over most of the disk of the Galaxy, the H I must be in two phases: the weight of the H I in the gravitational potential of the Galaxy is large enough to generate thermal pressures exceeding P_min, so that turbulent pressure fluctuations can produce cold gas that is thermally stable; and the mean density of the H I is too low for the gas to be all CNM. Our models predict the presence of CNM gas to R = 16-18 kpc, somewhat farther than previous estimates. We also examine the potential impact of turbulent heating on our results and provide expressions for the heating rate as a function of Galactic radius.

We report ISO-LWS far infrared observations of CO, water and oxygen lines towards the protobinary system IRAS4 in the NGC1333 cloud. We detected several water, OH, CO rotational lines, and two [OI] and [CII] fine structure lines. Given the relatively poor spectral and spatial resolution of these observations, assessing the origin of the observed emission is not straightforward. In this paper, we focus on the water line emission and explore the hypothesis that it originates in the envelopes that surround the two protostars, IRAS4 A and B, thanks to an accurate model. The model reproduces quite well the observed water line fluxes, predicting a density profile, mass accretion rate, central mass, and water abundance profile in agreement with previous works. We hence conclude that the emission from the envelopes is a viable explanation for the observed water emission, although we cannot totally rule out the alternative that the observed water emission originates in the outflow.

We report the first detection of doubly-deuterated methanol (CHD2OH), as well as firm detections of the two singly-deuterated isotopomers of methanol (CH2DOH and CH3OD), towards the solar-type protostar IRAS16293-2422. From the present multifrequency observations, we derive the following abundance ratios: [CHD2OH]/[CH3OH] = 0.2 +/- 0.1, [CH2DOH]/[CH3OH] = 0.9 +/- 0.3, [CH3OD]/[CH3OH] = 0.04 +/- 0.02. The total abundance of the deuterated forms of methanol is greater than that of its normal hydrogenated counterpart in the circumstellar material of IRAS16293-2422, a circumstance not previously encountered. Formaldehyde, which is thought to be the chemical precursor of methanol, possesses a much lower fraction of deuterated isotopomers (~ 20%) with respect to the main isotopic form in IRAS16293-2422. The observed fractionation of methanol and formaldehyde provides a severe challenge to both gas-phase and grain-surface models of deuteration. Two examples of the latter model are roughly in agreement with our observations of CHD2OH and CH2DOH if the accreting gas has a large (0.2-0.3) atomic D/H ratio. However, no gas-phase model predicts such a high atomic D/H ratio, and hence some key ingredient seems to be missing.

We have developed a model for molecular hydrogen formation under astrophysically relevant conditions. This model takes fully into account the presence of both physisorbed and chemisorbed sites on the surface, allows quantum mechanical diffusion as well as thermal hopping for absorbed H-atoms, and has been benchmarked versus recent laboratory experiments on H2 formation on silicate surfaces. The results show that H2 formation on grain surface is efficient in the interstellar medium up to some 300K. At low temperatures (<100K), H2 formation is governed by the reaction of a physisorbed H with a chemisorbed H. At higher temperatures, H2 formation proceeds through reaction between two chemisorbed H atoms. We present simple analytical expressions for H2 formation which can be adopted to a wide variety of surfaces once their surfaces characteristics have been determined experimentally.

The high resolution (R=25,000) infrared M band spectrum of the massive protostar NGC 7538 : IRS9 shows a narrow absorption feature at 4.779 um (2092.3 cm-1) which we attribute to the vibrational stretching mode of the 13CO isotope in pure CO icy grain mantles. This is the first detection of solid 13CO in the interstellar medium. The 13CO band is a factor of 2.3 narrower than the apolar component of the 12CO band. With this in mind, we discuss the mechanisms that broaden solid state absorption bands. It is shown that ellipsoidally shaped pure CO grains fit the bands of both isotopes at the same time. Reasonable fits are also obtained by CO embedded in N2-rich ices and thermally processed O2-rich ices. In addition, we report new insights into the the nature and evolution of interstellar CO ices by comparing the high quality multi-component solid 12CO spectrum of NGC 7538 : IRS9 with that of the previously studied low mass source L1489 IRS. The narrow absorption by apolar CO ices is present in both spectra, but much stronger in NGC 7538. It is superposed on a smooth broad absorption feature well fitted by a combination of CO2 and H2O-rich laboratory CO ices. The abundances of the latter two ices relative to the total H2O ice column are the same in both sources. We thus suggest that thermal processing manifests itself as evaporation of apolar ices only, and not the formation of CO2 or polar ices. Finally, the decomposition of the 12CO band is used to derive the 12CO/13CO abundance ratio in apolar ices. A ratio of 71+/-15 is deduced, in good agreement with gas phase CO studies (~77) and the solid 12CO2/13CO2 ratio of 80+/-11 found in the same line of sight. The implications for the chemical path along which CO2 is formed are discussed.

Observations of infrared fine-structure lines provide direct information on the metallicity and ionization structure of HII regions and indirectly on the hardness of the radiation field ionizing these nebulae. We have analyzed a sample of Galactic and Magellanic Cloud HII regions observed by the Infrared Space Observatory (ISO) to examine the interplay between stellar content, metallicity and the ionization structure of HII regions. The observed [SIV]10.5/[SIII]18.7 mum and [NeIII]15.5/[NeII]12.8 mum line ratios are shown to be highly correlated over more than two orders of magnitude. We have compared the observed line ratios to the results of photoionization models using different stellar energy distributions. The derived characteristics of the ionizing star depend critically on the adopted stellar model as well as the (stellar) metallicity. We have compared the stellar effective temperatures derived from these model studies for a few well-studied HII regions with published direct spectroscopic determinations of the spectral type of the ionizing stars. This comparison supports our interpretation that stellar and nebular metallicity influences the observed infrared ionic line ratios. We can explain the observed increase in degree of ionization, as traced by the [SIV]\[SIII] and [NeIII]\[NeII] line ratios, by the hardening of the radiation field due to the decrease of metallicity. The implications of our results for the determination of the ages of starbursts in starburst galaxies are assessed.

Here, we present 6-9 um spectra of a sample of RNe, HII regions, YSOs, evolved stars and galaxies that show strong unidentified infrared bands, obtained with the ISO-SWS. The IR emission features show pronounced variations. 1) The 6.2 um feature shifts and clearly shows profile variations. 2) The 7.7 um complex is comprised of at least two subpeaks peaking at 7.6 and one longwards of 7.7 um. In some cases the main peak can apparently shift up to 8 um. Two sources do not exhibit a 7.7 um complex but instead show a broad emission feature at 8.22 um. 3) The 8.6 um feature has a symmetric profile in all sources and some sources exhibit this band at slightly longer wavelengths. For the 6.2, 7.7 and 8.6 um features, the sources have been classified independently based on their profile and peak position. The classes derived for these features are directly linked with each other. Furthermore, the observed 6--9 um spectrum depends on the type of object. We attribute the observed 6.2 um profile and peak position to the combined effect of a PAH family and anharmonicity with pure PAHs representing the 6.3 um component and substituted/complexed PAHs representing the 6.2 um component. The 7.6 um component is well reproduced by both pure and substituted/complexed PAHs but the 7.8 um component remains an enigma. In addition, the exact identification of the 8.22 um feature remains unknown. The observed variations in the characteristics of the IR emission bands are linked to the local physical conditions. Possible formation and evolution processes that may influence the interstellar PAH class are highlighted.

We present 2-45 micron spectra of a large sample of carbon-rich evolved stars in order to study the ``30'' micron feature. We find the ``30'' micron feature in sources in a wide range of sources: low mass loss carbon stars, extreme carbon-stars, post-AGB objects and planetary nebulae. We extract the profiles from the sources by using a simple systematic approach to model the continuum. We find large variations in the wavelength and width of the extracted profiles of the ``30'' micron feature. We modelled the whole range of profiles in a simple way by using magnesium sulfide (MgS) dust grains with a MgS grain temperature different from the continuum temperature. The systematic change in peak positions can be explained by cooling of MgS grains as the star evolves off the AGB. In several sources we find that a residual emission excess at ~26 micron can also be fitted using MgS grains but with a different grains shape distribution. The profiles of the ``30'' micron feature in planetary nebulae are narrower than our simple MgS model predicts. We discuss the possible reasons for this difference. We find a sample of warm carbon-stars with very cold MgS grains. We discuss possible causes for this phenomenon. We find no evidence for rapid destruction of MgS during the planetary nebula phase and conclude that the MgS may survive to be incorporated in the ISM.

Mid infrared spectra provide a powerful probe of the conditions in dusty galactic nuclei. They variously contain emission features associated with star forming regions and absorptions by circumnuclear silicate dust plus ices in cold molecular cloud material. Here we report the detection of 6-8um water ice absorption in 18 galaxies observed by ISO. While the mid-IR spectra of some of these galaxies show a strong resemblance to the heavily absorbed spectrum of NGC 4418, other galaxies in this sample also show weak to strong PAH emission. The 18 ice galaxies are part of a sample of 103 galaxies with good S/N mid-IR ISO spectra. Based on our sample we find that ice is present in most of the ULIRGs, whereas it is weak or absent in the large majority of Seyferts and starburst galaxies. This result is consistent with the presence of larger quantities of molecular material in ULIRGs as opposed to other galaxy types. Like NGC 4418, several of our ice galaxy spectra show a maximum near 8um that is not or only partly due to PAH emission. While this affects only a small part of the galaxy population studied by ISO, it stresses the need for high S/N data and refined diagnostic methods, to properly discriminate spectra dominated by PAH emission and spectra dominated by heavy obscuration. The spectral variation from PAH emission to absorbed continuum emission near 8um shows strong similarities with Galactic star forming clouds. This leads us to believe that our classification of ice galaxy spectra might reflect an evolutionary sequence from strongly obscured beginnings of star formation (and AGN activity) to a less enshrouded stage of advanced star formation (and AGN activity), as the PAHs get stronger and the broad 8um feature weakens.

We have carried out a quantitative trend analysis of the crystalline silicates observed in the ISO spectra of a sample of 14 stars with different evolutionary backgrounds. We have modeled the spectra using a simple dust radiative transfer model and have correlated the results with other known parameters. We confirm the abundance difference of the crystalline silicates in disk and in outflow sources, as found by Molster et al. (1999, Nature 401, 563). We found some indication that the enstatite over forsterite abundance ratio differs, it is slightly higher in the outflow sources with respect to the disk sources. It is clear that more data is required to fully test this hypothesis. We show that the 69.0 micron feature, attributed to forsterite, may be a very suitable temperature indicator. We found that the enstatite is more abundant than forsterite in almost all sources. The temperature of the enstatite grains is about equal to that of the forsterite grains in the disk sources but slightly lower in the outflow sources. Crystalline silicates are on average colder than amorphous silicates. This may be due to the difference in Fe content of both materials. Finally we find an indication that the ratio of ortho to clino enstatite, which is about 1:1 in disk sources, shifts towards ortho enstatite in the high luminosity (outflow) sources.

This is the first paper in a series of three where we present the first comprehensive inventory of solid state emission bands observed in a sample of 17 oxygen-rich circumstellar dust shells surrounding evolved stars. The data were taken with the Short and Long Wavelength Spectrographs on board of the Infrared Space Observatory (ISO) and cover the 2.4 to 195 micron wavelength range. The spectra show the presence of broad 10 and 18 micron bands that can be attributed to amorphous silicates. In addition, at least 49 narrow bands are found whose position and width indicate they can be attributed to crystalline silicates. Almost all of these bands were not known before ISO. We have measured the peak positions, widths and strengths of the individual, continuum subtracted bands. Based on these measurements, we were able to order the spectra in sequence of decreasing crystalline silicate band strength. We found that the strength of the emission bands correlates with the geometry of the circumstellar shell, as derived from direct imaging or inferred from the shape of the spectral energy distribution. This naturally divides the sample into objects that show a disk-like geometry (strong crystalline silicate bands), and objects whose dust shell is characteristic of an outflow (weak crystalline silicate bands). All stars with the 33.6 micron forsterite band stronger than 20 percent over continuum are disk sources. We define spectral regions (called complexes) where a concentration of emission bands is evident, at 10, 18, 23, 28, 33, 40 and 60 micron. We derive average shapes for these complexes and compare these to the individual band shapes of the programme stars.

A (sub-)millimeter line and continuum study of the class I protostar Elias 29 in the Rho Ophiuchi molecular cloud is presented, whose goals are to understand the nature of this source, and to locate the ices that are abundantly present along this line of sight. Within 15-60" beams, several different components contribute to the line emission. Two different foreground clouds are detected, an envelope/disk system and a dense ridge of HCO+ rich material. The latter two components are spatially separated in millimeter interferometer maps. We analyze the envelope/disk system by using inside-out collapse and flared disk models. The disk is in a relatively face-on orientation (<60 degrees), which explains many of the remarkable observational features of Elias 29, such as its flat SED, its brightness in the near infrared, the extended components found in speckle interferometry observations, and its high velocity molecular outflow. It cannot account for the ices seen along the line of sight, however. A small fraction of the ices is present in a (remnant) envelope of mass 0.12-0.33 Msun, but most of the ices (~70%) are present in cool (T<40 K) quiescent foreground clouds. This explains the observed absence of thermally processed ices (crystallized H2O) toward Elias 29. Nevertheless, the temperatures could be sufficiently high to account for the low abundance of apolar (CO, N2, O2) ices. This work shows that it is crucial to obtain spectrally and spatially resolved information from single-dish and interferometric molecular gas observations in order to determine the nature of protostars and to interpret infrared ISO satellite observations of ices and silicates along a pencil beam.

This is the second paper in a series of three in which we present an exhaustive inventory of the 49 solid state emission bands observed in a sample of 17 oxygen-rich dust shells surrounding evolved stars. Most of these emission bands are concentrated in well defined spectral regions (called complexes). We define 7 of these complexes; the 10, 18, 23, 28, 33, 40 and 60 micron complex. We derive average properties of the individual bands. Comparison with laboratory data suggests that both olivines (Mg(2x)Fe(2-2x)SiO(4)) and pyroxenes (Mg(x)Fe(1-x)SiO(3)) are present, with x close to 1, i.e. the minerals are very Mg-rich and Fe-poor. This composition is similar to that seen in disks surrounding young stars and in the solar system comet Hale-Bopp. A significant fraction of the emission bands cannot be identified with either olivines or pyroxenes. Possible other materials that may be the carriers of these unidentified bands are briefly discussed. There is a natural division into objects that show a disk-like geometry (strong crystalline silicate bands), and objects whose dust shell is characteristic of an outflow (weak crystalline silicate bands). In particular, stars with the 33.5 micron olivine band stronger than about 20 percent over continuum are invariably disk sources. Likewise, the 60 micron region is dominated by crystalline silicates in the disk sources, while it is dominated by crystalline H(2)O ice in the outflow sources. We show that the disk and outflow sources have significant differences in the shape of the emission bands. This difference must be related to the composition or grain shapes of the dust particles. The incredible richness of the crystalline silicate spectra observed by ISO allows detailed studies of the mineralogy of these dust shells, and is the origin and history of the dust.

We present a fit to the spectral energy distribution of OH 127.8+0.0, a typical asymptotic giant branch star with an optically thick circumstellar dust shell. The fit to the dust spectrum is achieved using non-spherical grains consisting of metallic iron, amorphous and crystalline silicates and water ice. Previous similar attempts have not resulted in a satisfactory fit to the observed spectral energy distributions, mainly because of an apparent lack of opacity in the 3--8 micron region of the spectrum. Non-spherical metallic iron grains provide an identification for the missing source of opacity in the near-infrared. Using the derived dust composition, we have calculated spectra for a range of mass-loss rates in order to perform a consistency check by comparison with other evolved stars. The L-[12 micron] colours of these models correctly predict the mass-loss rate of a sample of AGB stars, strengthening our conclusion that the metallic iron grains dominate the near-infrared flux. We discuss a formation mechanism for non-spherical metallic iron grains.

Based on the ISO spectral catalogue of compact HII regions by Peeters et al. (2001), we present a first analysis of the hydrogen recombination and atomic fine-structure lines originated in the ionized gas. The sample consists of 34 HII regions located at galactocentric distances between Rgal = 0 and 15 kpc. The SWS HI recombination lines between 2 and 8 mum are used to estimate the extinction law at these wavelengths for 14 HII regions. An extinction in the K band between 0 and $\sim$ 3 mag. has been derived. The fine-structure lines of N, O, Ne, S and Ar are detected in most of the sources. Most of these elements are observed in two different ionization stages probing a range in ionization potential up to 41 eV. The ISO data, by itself or combined with radio data taken from the literature, is used to derive the elemental abundances relative to hydrogen. The present data thus allow us to describe for each source its elemental abundance, its state of ionization and to constrain the properties of the ionizing star(s).

The first results from a VLT-ISAAC program on L- and M-band infrared spectroscopy of deeply-embedded young stellar objects are presented. The advent of 8-m class telescopes allows high S/N spectra of low-luminosity sources to be obtained. In our first observing run, low- and medium-resolution spectra have been measured toward a dozen objects, mostly in the Vela and Chamaeleon molecular clouds. The spectra show strong absorption of H2O and CO ice, as well as weak features at `3.47' and 4.62 mu. No significant solid CH3OH feature at 3.54 mu is found, indicating that the CH3OH/H2O ice abundance is lower than toward some massive protostars. Various evolutionary diagnostics are investigated for a set of sources in Vela.

We report the discovery of the "21" micron and "30" micron features in the planetary nebulae around the hydrogen-deficient stars HD 826 and HD 158269. The carriers of these features are known to be produced in outflows around carbon-rich stars. This discovery demonstrates that the bulk of the dust in these nebulae has been produced during a carbon-rich phase before the atmospheres of these stars became hydrogen poor. This is the first time that the "21" micron feature has been detected in any planetary nebula. It shows that once formed its carrier can survive the formation of the nebula and the exposure to the UV radiation of the hot central star. This means that the carrier of "21" micron feature is not transient: the absence of the feature sets limits on the production of its carrier.

We present 10-15 micron spectra of a sample of H II regions, YSOs and evolved stars that show strong unidentified infrared emission features, obtained with the ISO/SWS spectrograph on-board ISO. These spectra reveal a plethora of emission features with bands at 11.0, 11.2, 12.0, 12.7, 13.5 and 14.2 micron. These features are observed to vary considerably in relative strength to each-other from source to source. In particular, the 10-15 micron spectra of the evolved stars are dominated by the 11.2 micron band while for H II regions the 12.7 is typically as strong as the 11.2 micron band. Analysing the ISO data we find a good correlation between the 11.2 micron band and the 3.3 micron band, and between the 12.7 micron and the 6.2 micron band. There is also a correlation between the ratio of the UIR bands to the total dust emission and the 12.7 over 11.2 micron ratio. Bands in the 10-15 micron spectral region are due to CH out-of-plane (OOP) bending modes of polycyclic aromatic hydrocarbons (PAHs). We summarise existing laboratory data and theoretical quantum chemical calculations of these modes for neutral and cationic PAHs. Due to mode coupling, the exact peak position of these bands depends on the number of adjacent CH groups and hence the observed interstellar 10-15 micron spectra can be used to determine the molecular structure of the interstellar PAHs emitting in the different regions. We conclude that evolved stars predominantly inject compact ~100-200 C-atom PAHs into the ISM where they are subsequently processed, resulting in more open and uneven PAH structures.

Infrared Space Observatory (ISO) observations have shown that O-rich Asymptotic Giant Branch (AGB) stars exhibit crystalline silicate features in their spectra only if their mass-loss rate is higher than a certain threshold value. Usually, this is interpreted as evidence that crystalline silicates are not present in the dust shells of low mass-loss rate objects. In this study, radiative transfer calculations have been performed to search for an alternative explanation to the lack of crystalline silicate features in the spectrum of low mass-loss rate AGB stars. It is shown that due to a temperature difference between amorphous and crystalline silicates it is possible to include up to 40% of crystalline silicate material in the circumstellar dust shell, without the spectra showing the characteristic spectral features. Since this implies that low mass-loss rate AGB stars might also form crystalline silicates and deposit them into the Interstellar Medium (ISM), the described observational selection effect may put the process of dust formation around AGB stars and the composition of the predominantly amorphous dust in the Interstellar Medium in a different light. Our model calculations result in a diagnostic tool to determine the crystallinity of an AGB star with a known mass-loss rate.

We report the detection of absorption features in the 6-8 micron region superimposed on a featureless mid-infrared continuum in NGC 4418. For several of these features this is the first detection in an external galaxy. We compare the absorption spectrum of NGC 4418 to that of embedded massive protostars and the Galactic centre, and attribute the absorption features to ice grains and to hydrogenated amorphous carbon grains. From the depth of the ice features, the powerful central source responsible for the mid-infrared emission must be deeply enshrouded. Since this emission is warm and originates in a compact region, an AGN must be hiding in the nucleus of NGC 4418.

The VLA and the BIMA array were used to obtain high resolution (10-20") observations of C+, traced by the C91a recombination line at 8.6 GHz, and the dense molecular gas, traced by HCN and HCO+(1-0), of the photon dominated region (PDR) associated with the reflection nebula NGC 2023. Using the VLA, continuum emission is detected at 8.6 GHz from a faint HII region associated with HD 37903. The C91a emission originates from a 0.4 pc long filament, extending from the east to the south of the exciting star HD 37903. Within the filament three C91a clumps can be distinguished, each associated with filamentary vibrationally excited H2 emission in the direction toward HD 37903. The HCO+ emission has a clumpy appearance superimposed on a more extended component. C91a is, in general, closer to the exciting star than HCO+ emission as expected from PDR models. The morphologies of HCO+ and HCN are quite similar. Based on the C91a linewidth towards one of the clumps a limit of 170 K on the kinetic temperature in the ionized carbon layer can be derived. This value is consistent with PDR models with H2 densities of about 10^5/cm^3. However, this result suggests surprisingly low limits on the turbulence in the PDR. We detected a compact 3 mm continuum source in the PDR, which appears to be a cold "core" of density 10^7/cm^3, 0.03 parsec diameter, and 6 solar masses. We conclude that it may have formed within the PDR. In an appendix, observations of the C91a recombination line toward five additional PDRs using the Effelsberg 100m telescope are described.

We have obtained the full 1-200 um spectrum of the low luminosity (36 Lsun) Class I protostar Elias 29 in the Rho Ophiuchi molecular cloud. It provides a unique opportunity to study the origin and evolution of interstellar ice and the interrelationship of interstellar ice and hot core gases around low mass protostars. We see abundant hot CO and H2O gas, as well as the absorption bands of CO, CO2, H2O and ``6.85 um'' ices. We compare the abundances and physical conditions of the gas and ices toward Elias 29 with the conditions around several well studied luminous, high mass protostars. The high gas temperature and gas/solid ratios resemble those of relatively evolved high mass objects (e.g. GL 2591). However, none of the ice band profiles shows evidence for significant thermal processing, and in this respect Elias 29 resembles the least evolved luminous protostars, such as NGC 7538 : IRS9. Thus we conclude that the heating of the envelope of the low mass object Elias 29 is qualitatively different from that of high mass protostars. This is possibly related to a different density gradient of the envelope or shielding of the ices in a circumstellar disk. This result is important for our understanding of the evolution of interstellar ices, and their relation to cometary ices.

We present the results of ISO SWS and LWS grating scans towards the embedded Young Stellar Objects (YSOs) S106 IR and Cep A East. Emission from the pure rotational lines of H2 and the infrared fine structure lines of [C II], [O I], [S I], [Si II] and [Fe II], as well as absorption bands due to H2O, CO and CO2 ice were detected toward Cep A. In S106 we detected emission lines of H2, CO, H I, and a large number of ionized species including Fe, O, N, C, Si, S, Ne and Ar. S106 also shows many of the infrared PAH bands in emission. Excitation temperatures and molecular hydrogen masses were derived from the low-lying pure rotational levels of H2 and are 500 and 730 K and 8 and 3 x 10^{-3} solar masses for S106 and Cep A, respectively. Since both objects are expected to have several solar masses of H2 in their environment, we conclude that in both cases the bulk of the H2 is cooler than a few hundred Kelvins. Excitation temperatures and line ratios were compared with those predicted by theoretical models for PDRs and dissociative and non-dissociative shocks. The [S I] 25.2 micron/[Si II] 34.8 micron ratio is a particularly useful shock versus PDR discriminant and we conclude that S106 IR is dominated by PDR emission while Cep A East has a large shock component. From an analysis of the ionic lines in S106 we conclude that the central star must have a temperature around 37,000 K, corresponding to a spectral type of O8. From its luminosity it is concluded that the driving source of Cep A must also be a massive early-type star. The absence of strong high-ionization ionic lines in its ISO spectrum shows that Cep A has not yet created a significant H II region and must be younger than S106, illustrating the process of the clearing of the surroundings of a massive young star.

We have observed the central region of the nearby starburst galaxy NGC4945 with the mid-infrared spectrometers SWS and PHT-S aboard ISO. We do not find any evidence for the existence of the powerful AGN, inferred from hard X-ray observations. The upper limits on our AGN tracers [NeV] 14.32um & 24.32um and [NeVI] 7.65um imply an A(V)>160 towards the NLR, assuming the NLR to be of equal strength as in the Circinus galaxy. Other possibilities are discussed. The starburst excitation indicators [NeIII]15.56um/[NeII]12.81um and L(bol)/L(lyc) suggest that the starburst in the central region is at least 5x10^6 yrs old, and that it accounts for at least half of the nuclear bolometric luminosity. The starburst might well power the entire bolometric luminosity, but the available constraints are also consistent with an up to 50% contribution of the embedded AGN. With PHT-S, at a resolution of approximately 90, we detect strong absorption features of water ice, and, for the first time in an external galaxy, of CO_2 and CO. The same PHT-S spectrum also reveals strong emission from the family of PAH features. Finally, we have observed and detected several pure rotational and ro-vibrational H_2 lines, two of which, the (0-0) S(0) & S(1) lines, allow us to determine the excitation temperature (160K) and warm H_2 mass (2.4x10^7 Msun). The low excitation temperature shows Orion-like shocks not to be representative for the entire emission of the central region of the galaxy and fairly normal PDRs to be perhaps more typical.

We report results from a Hubble Space Telescope (HST) and Near-Infrared Camera and Multiobject Spectrometer (NICMOS) program to study the distribution of hot neutral (molecular hydrogen) and ionized circumstellar material in the young planetary nebulae NGC 7027. HST/NICMOS provided very high spatial resolution imaging in line and continuum emission, and the stability and large dynamic range needed for investigating detailed structures in the circumstellar material. We present dramatic new images of NGC 7027 that have led to a new understanding of the structure in this important planetary nebula. The central star is clearly revealed, providing near-infrared fluxes that are used to directly determine the stellar temperature very accurately (T = 198,000 K). It is found that the photodissociation layer as revealed by near-infrared molecular hydrogen emission is very thin, and is biconical in shape. The interface region is structured and filamentary, suggesting the existence of hydrodynamic instabilities. We discuss evidence for the presence of one or more highly collimated, off-axis jets that might be present in NGC 7027. The evolutionary future of NGC 7027 is discussed.

We discuss the composition of dust and ice along the line of sight to the Galactic Center (GC) based on analysis of mid-infrared spectra (2.4-13 micron) from the Short Wavelength Spectrometer on the Infrared Space Observatory (ISO). We have analyzed dust absorption features arising in the molecular cloud material and the diffuse interstellar medium along the lines of sight toward Sagittarius A* and the Quintuplet sources, GCS3 and GCS4. It is evident from the depth of the 3.0 micron H2O and the 4.27 micron CO2 ice features that there is more molecular cloud material along the line of sight toward Sgr A* than GCS3 and 4. In fact, Sgr A* has a rich infrared ice spectrum with evidence for the presence of solid CH4, NH3, and possibly HCOOH. Hydrocarbon dust in the diffuse interstellar medium along the line of sight to the GC is characterized by absorption features centered at 3.4, 6.85, and 7.25 micron. Ground-based studies have identified the 3.4 micron feature with aliphatic hydrocarbons, and ISO has given us the first meaningful observations of the corresponding modes at longer wavelengths. The integrated strengths of these three features suggest that hydrogenated amorphous carbon is their carrier. We attribute an absorption feature centered at 3.28 micron in the GCS3 spectrum to the C-H stretch in aromatic hydrocarbons. This feature is not detected, and its C-C stretch counterpart appears to be weaker, in the SgrA* spectrum. One of the key questions which now arises is whether aromatics are a widespread component of the diffuse interstellar medium, analogous to aliphatic hydrocarbons.

Images of Hubble 12 were obtained with the HST/NICMOS instrument in the F110W, F164N, and F166N filters of NIC1, and the F160W, F187N, F190N, F212N, F215N filters of NIC2. The images show the structure of the inner ``torus'' and lobes much clearer than previous ground-based images. In particular, the [Fe II] lobes are clearly resolved and shown to be distinct from the molecular hydrogen emission structures. Apparent changes in the inner ionized and neutral bipolar shell implies periodic mass loss or changes in stellar wind shape or direction. The molecular hydrogen in the ``eye'' is radiatively excited and shows a complex morphology that suggests that several mass ejection events are responsible for producing this structure. The position angle of the molecular hydrogen and [Fe II] lobes differ, indicating a possible precession of the ejection axis.

We present observations of the stretching mode of 13CO2 ice along 13 lines of sight in the Galaxy, using the Short Wavelength Spectrometer on board of the Infrared Space Observatory. Remarkable variations are seen in the absorption band profile in the different lines of sight. The main feature is attributed to 13CO2 mixed with polar molecules such as H2O, and CH3OH. Several high-mass protostars show an additional narrow substructure at 4.382 um, which we attribute to a polar CO2 ice, that experienced heating. This heating effect is sustained by a good correlation of the strength of the substructure with dust and CO gas temperatures along the line of sight, and anti-correlation with ice abundances. Thus, our main conclusion is that interstellar CO2 ices around luminous protostars are subjected to, and altered by, thermal processing and that it may reflect the evolutionary stage of the nearby protostar. In contrast, the ices around low mass protostars and in a quiescent cloud in our sample do not show signs of thermal processing. Furthermore, we determine for the first time the Galactic 12C/13C ratio from the solid state. The 12CO2/13CO2 ratio for the local ISM (69+/-15), as well as the dependence on Galacto-centric radius, are in good agreement with gas phase (C18O, H2CO) studies. For the few individual objects for which gas phase values are available, the 12C/13C ratios derived from CO2 tend to be higher compared to CO studies (albeit with 2.5 sigma significance only). We discuss the implications of this possible difference for the chemical origin of interstellar CO2.

We have studied molecular hydrogen emission in a sample of 21 YSOs using spectra obtained with the Infrared Space Observatory (ISO). H2 emission was detected in 12 sources and can be explained as arising in either a shock, caused by the interaction of an outflow from an embedded YSO with the surrounding molecular cloud, or in a PDR surrounding an exposed young early-type star. The distinction between these two mechanisms can not always be made from the pure rotational H2 lines alone. Other tracers, such as PAH emission or [S I] 25.25 micron emission, are needed to identify the H2 heating mechanism. No deviations from a 3:1 ortho/para ratio of H2 were found. Both shocks and PDRs show a warm and a hot component in H2, which we explain as thermal emission from warm molecular gas (warm component), or UV-pumped infrared fluorescence in the case of PDRs and the re-formation of H2 for shocks (hot component).

We present the results of ISO SWS and LWS grating scans towards the three brightest members of the BD+40 4124 group in the infrared: BD+40 4124 (B2Ve), LkHalpha 224 (A7e) and the embedded source LkHalpha 225. Emission from the pure rotational lines of H_2, from ro-vibrational transitions of CO, from PAHs, from H I recombination lines and from the infrared fine structure lines of [Fe II], [Si II], [S I], [O I], [O III] and [C II] was detected. These emission lines arise in the combination of a low-density (approx 10^2 cm^-3) H II region with a clumpy PDR in the case of BD+40 4124. The lower transitions of the infrared H I lines observed in BD+40 4124 are optically thick; most likely they arise in either a dense wind or a circumstellar disk. This same region is also responsible for the optical H I lines and the radio continuum emission. In the lines of sight towards LkHalpha 224 and LkHalpha 225, the observed emission lines arise in a non-dissociative shock produced by a slow (approx 20 km/s) outflow arising from LkHalpha 225. Toward LkHalpha 225 we also observe a dissociative shock, presumably located closer to the outflow source than the non-dissociative shock. In the line of sight towards LkHalpha 225 we observed absorption features due to solid water ice and amorphous silicates, and due to gas-phase H_2O, CO and CO_2. No solid CO_2 was detected towards LkHalpha 225, making this the first line of sight where the bulk of the CO_2 is in the gas-phase.

Infrared spectra of a number of evolved O-rich stars have been obtained with the Short- and Long- Wavelength spectrometers on board the Infrared Space Observatory. The very broad wavelength coverage (2.4-197 micron) obtained by combining observations made with the two spectrometers includes practically all of the flux emitted by the sources, and allows us to determine the emission and absorption features of the dense circumstellar dust shells. Agreement between the fluxes obtained by the two instruments is generally very good; the largest discrepancies are probably due to source variability. Our sample of oxygen-rich AGB stars exhibits a wealth of spectral features due to crystalline silicates and crystalline water ice in emission and absorption. In this study a qualitative overview of all features due to crystalline silicates and water ice in these high mass loss rate objects is presented. It seems that there is a certain onset value for the mass loss rate above which these features appear in the spectrum. Moreover, crystalline silicate emission features have been detected for the first time at wavelengths where the amorphous silicates are still in absorption, implying different spatial distributions for the two materials. A spherically symmetric and an axi-symmetric geometry are proposed.

In this paper we discuss the propagation of dust through the interstellar medium (ISM), and describe the destructive effects of stellar winds, jets, and supernova shock waves on interstellar dust. We review the probability that grains formed in stellar outflows or supernovae survive processing in and propagation through the ISM, and incorporate themselves relatively unprocessed into meteoritic bodies in the solar system. We show that very large (radii >= 5 micron) and very small grains (radii <= 100 Angstrom) with sizes similar to the pre-solar SiC and diamond grains extracted from meteorites, can survive the passage through 100\kms shock waves relatively unscathed. High velocity (>= 250 km/s) shocks destroy dust efficiently. However, a small (~10%) fraction of the stardust never encountered such fast shocks before incorporation into the solar system. All grains should therefore retain traces of their passage through interstellar shocks during their propagation through the ISM. The grain surfaces should show evidence of processing due to sputtering and pitting due to small grain cratering collisions on the micron-sized grains. This conclusion seems to be in conflict with the evidence from the large grains recovered from meteorites which seem to show little interstellar processing.

We reexamine the spectroscopic underpinnings of recent suggestions that [O I] and [Fe II] lines from the Orion H II region are produced in gas where the iron-carrying grains have been destroyed and the electron density is surprisingly high. Our new observations show that previous detections of [O I] 5577 were dominated by telluric emission. Our limits are consistent with a moderate density (10^4 cm^{-3}) photoionized gas. We show that a previously proposed model of the Orion H II region reproduces the observed [O I] and [Fe II] spectrum. These lines are fully consistent with formation in a dusty region of moderate density.