Interstellar carbonaceous dust is mainly formed in the innermost regions of circumstellar envelopes around carbon-rich asymptotic giant branch (AGB) stars. In these highly chemically stratified regions, atomic and diatomic carbon, along with acetylene are the most abundant species after H$_2$ and CO. In a previous study, we addressed the chemistry of carbon (C and C$_2$) with H$_2$ showing that acetylene and aliphatic species form efficiently in the dust formation region of carbon-rich AGBs whereas aromatics do not. Still, acetylene is known to be a key ingredient in the formation of linear polyacetylenic chains, benzene and polycyclic aromatic hydrocarbons (PAHs), as shown by previous experiments. However, these experiments have not considered the chemistry of carbon (C and C$_2$) with C$_2$H$_2$. In this work, by employing a sufficient amount of acetylene, we investigate its gas-phase interaction with atomic and diatomic carbon. We show that the chemistry involved produces linear polyacetylenic chains, benzene and other PAHs, which are observed with high abundances in the early evolutionary phase of planetary nebulae. More importantly, we have found a non-negligible amount of pure and hydrogenated carbon clusters as well as aromatics with aliphatic substitutions, both being a direct consequence of the addition of atomic carbon. The incorporation of alkyl substituents into aromatics can be rationalized by a mechanism involving hydrogen abstraction followed by methyl addition. All the species detected in gas phase are incorporated into the nanometric sized dust analogues, which consist of a complex mixture of sp, sp$^2$ and sp$^3$ hydrocarbons with amorphous morphology.

As a part of interstellar dust, polycyclic aromatic hydrocarbons (PAHs) are processed by the interaction with vacuum ultraviolet (VUV) photons that are emitted by hot young stars. This interaction leads to the emission of the well-known aromatic infrared bands but also of electrons, which can significantly contribute to the heating of the interstellar gas.Our aim is to investigate the impact of molecular size on the photoionization properties of cationic PAHs.Trapped PAH cations of sizes between 30 and 48 carbon atoms were submitted to VUV photons in the range of 9 to 20 eV from the DESIRS beamline at the synchrotron SOLEIL. All resulting photoproducts including dications and fragment cations were mass-analyzed and recorded as a function of photon energy.Photoionization is found to be predominant over dissociation at all energies, which differs from an earlier study on smaller PAHs. The photoionization branching ratio reaches 0.98 at 20 eV for the largest studied PAH. The photoionization threshold is observed to be between 9.1 and 10.2 eV, in agreement with the evolution of the ionization potential with size. Ionization cross sections were indirectly obtained and photoionization yields extracted from their ratio with theoretical photoabsorption cross sections, which were calculated using time-dependent density functional theory. An analytical function was derived to calculate this yield for a given molecular size.Large PAH cations could be efficiently ionized in H I regions and provide a contribution to the heating of the gas by photoelectric effect. Our work provides recipes to be used in astronomical models to quantify these points.

Polycyclic aromatic hydrocarbons (PAHs) are key species in astrophysical environments in which vacuum ultraviolet (VUV) photons are present, such as star-forming regions. The interaction with these VUV photons governs the physical and chemical evolution of PAHs. Models show that only large species can survive. However, the actual molecular properties of large PAHs are poorly characterized and the ones included in models are only an extrapolation of the properties of small and medium-sized species. We discuss here experiments performed on trapped ions including some at the SOLEIL VUV beam line DESIRS. We focus on the case of the large dicoronylene cation, C48H20+ , and compare its behavior under VUV processing with that of smaller species. We suggest that C2H2 is not a relevant channel in the fragmentation of large PAHs. Ionization is found to largely dominate fragmentation. In addition, we report evidence for a hydrogen dissociation channel through excited electronic states. Although this channel is minor, it is already effective below 13.6 eV and can significantly influence the stability of astro-PAHs. We emphasize that the competition between ionization and dissociation in large PAHs should be further evaluated for their use in astrophysical models.

The JWST will deliver large data sets of high-quality spectral data over the 0.6-28 $\mu$m range. It will combine sensitivity, spectral and spatial resolution. Specific tools are required to provide efficient scientific analysis of such large data sets. Our aim is to illustrate the potential of unsupervised learning methods to get insights into chemical variations in the populations that carry the aromatic infrared bands (AIBs), more specifically PAH species and carbonaceous very small grains (VSGs). We present a method based on linear fitting and blind signal separation (BSS) for extracting representative spectra for a spectral data set. The method is fast and robust, which ensures its applicability to JWST spectral cubes. We tested this method on a sample of ISO-SWS data, which resemble most closely the JWST spectra in terms of spectral resolution and coverage. Four representative spectra were extracted. Their main characteristics appear consistent with previous studies with populations dominated by cationic PAHs, neutral PAHs, evaporating VSGs, and large ionized PAHs, known as the PAH$^x$ population. In addition, the 3 $\mu$m range, which is considered here for the first time in a BSS method, reveals the presence of aliphatics connected to neutral PAHs. Each representative spectrum is found to carry second-order spectral signatures (e.g. small bands), which are connected with the underlying chemical diversity of populations. However, the precise attribution of theses signatures remains limited by the combined small size and heterogeneity of the sample of astronomical spectra available in this study. The upcoming JWST data will allow us to overcome this limitation. The large data sets of hyperspectral images provided by JWST analysed with the proposed method, which is fast and robust, will open promising perspectives for our understanding of the chemical evolution of the AIB carriers.

While powerful techniques exist to accurately account for anharmonicity in vibrational molecular spectroscopy, they are computationally very expensive and cannot be routinely employed for large species and/or at non- zero vibrational temperatures. Motivated by the study of Polycyclic Aromatic Hydrocarbon (PAH) emission in space, we developed a new code, which takes into account all modes and can describe all IR transitions including bands becoming active due to resonances as well as overtones, combination and difference bands. In this article, we describe the methodology that was implemented and discuss how the main difficulties were overcome, so as to keep the problem tractable. Benchmarking with high-level calculations was performed on a small molecule. We carried out specific convergence tests on two prototypical PAHs, pyrene (C$_{16}$H$_{10}$) and coronene (C$_{24}$H$_{12}$), aiming at optimising tunable parameters to achieve both acceptable accuracy and computational costs for this class of molecules. We then report the results obtained at 0 K for pyrene and coronene, comparing the calculated spectra with available experimental data. The theoretical band positions were found to be significantly improved compared to harmonic Density Functional Theory (DFT) calculations. The band intensities are in reasonable agreement with experiments, the main limitation being the accuracy of the underlying calculations of the quartic force field. This is a first step towards calculating moderately high-temperature spectra of PAHs and other similarly rigid molecules using Monte Carlo sampling.

In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines. We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J lines that trace warm molecular gas at PDR interfaces. We observed the 12CO and 13CO ladders in two prototypical PDRs, the Orion Bar and NGC 7023 NW using the instruments onboard Herschel. We also considered line emission from key species in the gas cooling of PDRs (C+, O, H2) and other tracers of PDR edges such as OH and CH+. All the intensities are collected from Herschel observations, the literature and the Spitzer archive and are analyzed using the Meudon PDR code. A grid of models was run to explore the parameter space of only two parameters: thermal gas pressure and a global scaling factor that corrects for approximations in the assumed geometry. We conclude that the emission in the high-J CO lines, which were observed up to Jup=23 in the Orion Bar (Jup=19 in NGC7023), can only originate from small structures of typical thickness of a few 1e-3 pc and at high thermal pressures (Pth~1e8 K cm-3). Compiling data from the literature, we found that the gas thermal pressure increases with the intensity of the UV radiation field given by G0, following a trend in line with recent simulations of the photoevaporation of illuminated edges of molecular clouds. This relation can help rationalising the analysis of high-J CO emission in massive star formation and provides an observational constraint for models that study stellar feedback on molecular clouds.

We present a proof of concept on the coupling of radio astronomical receivers and spectrometers with chemical reactorsand the performances of the resulting setup for spectroscopy and chemical simulations in laboratory astrophysics. Several experiments including cold plasma generation and UV photochemistry were performed in a 40\,cm long gas cell placed in the beam path of the Aries 40\,m radio telescope receivers operating in the 41-49 GHz frequency range interfaced with fast Fourier transform spectrometers providing 2 GHz bandwidth and 38 kHz resolution. The impedance matching of the cell windows has been studied using different materials. The choice of the material and its thickness was critical to obtain a sensitivity identical to that of standard radio astronomical observations. Spectroscopic signals arising from very low partial pressures of CH3OH, CH3CH2OH, HCOOH, OCS,CS, SO2 (<1E-03 mbar) were detected in a few seconds. Fast data acquisition was achieved allowing for kinetic measurements in fragmentation experiments using electron impact or UV irradiation. Time evolution of chemical reactions involving OCS, O2 and CS2 was also observed demonstrating that reactive species, such as CS, can be maintained with high abundance in the gas phase during these experiments.

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

Far-UV photons (E<13.6 eV) from hot massive stars regulate, or at least influence, the heating, ionization, and chemistry of most of the neutral interstellar medium. Investigating the interaction between FUV radiation and interstellar matter thus plays an important role in astrochemistry. We have used ALMA to mosaic a small field of the Orion Bar where the critical transition from atomic to molecular gas takes place. These observations provide an unprecedented sharp view of this transition layer (~1" resolution or ~414 AU). The resulting images (so far in the rotational emission of CO, HCO+, H13CO+, SO+, SO, and reactive ions SH+ and HOC+) show the small-scale structure in gas density and temperature, and the steep abundance gradients. The images reveal a pattern of high-density substructures, photo-ablative gas flows and instabilities at the edge of the molecular cloud. We have also used the IRAM 30m telescope to carry out a line-survey of the illuminated edge of the Bar in the mm domain. Our observations reveal the presence of complex organic molecules that were not expected in such a harsh environment. In particular, we have reported the first detection of the unstable cis conformer of formic acid (HCOOH) in the ISM. The energy barrier to internal rotation (the conversion from trans to cis) is approximately 4827 cm-1. Hence, this detection is surprising. The low inferred trans-to-cis abundance ratio of 2.8+/-1.0 supports a photoswitching mechanism: a given conformer absorbs a FUV stellar photon that radiatively excites the molecule to electronic states above the interconversion barrier. Subsequent fluorescent decay leaves the molecule in a different conformer form. This mechanism, which we have specifically studied with ab initio quantum calculations, was not considered so far in astrochemistry although it can affect the structure of a variety of molecules in PDRs. (abridged)

Emission of fullerenes in their infrared vibrational bands has been detected in space near hot stars. The proposed attribution of the diffuse interstellar bands at 9577 and 9632 \AA\ to electronic transitions of the buckminsterfullerene cation (i.e. C$_{60}^+$ ) was recently supported by new laboratory data, confirming the presence of this species in the diffuse interstellar medium (ISM). In this letter, we present the detection, also in the diffuse ISM, of the 17.4 and 18.9 $\mu$m emission bands commonly attributed to vibrational bands of neutral C$_{60}$. According to classical models that compute the charge state of large molecules in space, C$_{60}$ is expected to be mostly neutral in the diffuse ISM. This is in agreement with the abundances of diffuse C$_{60}$ we derive here from observations.

We developed a new analytical experimental setup called AROMA (Astrochemistry Research of Organics with Molecular Analyzer) that combines laser desorption/ionization techniques with ion trap mass spectrometry. We report here on the ability of the apparatus to detect aromatic species in complex materials of astrophysical interests and characterize their structures. A limit of detection of 100 femto-grams has been achieved using pure polycyclic aromatic hydrocarbon (PAH) samples, which corresponds to 2x10^8 molecules in the case of coronene (C24H12). We detected the PAH distribution in the Murchison meteorite, which is made of a complex mixture of extraterrestrial organic compounds. In addition, collision induced dissociation experiments were performed on selected species detected in Murchison, which led to the first firm identification of pyrene and its methylated derivatives in this sample.

Interstellar polycyclic aromatic hydrocarbons (PAHs) are strongly affected by the absorption of vacuum ultraviolet (VUV) photons in the interstellar medium (ISM), yet the branching ratio between ionization and fragmentation is poorly studied. This is crucial for the stability and charge state of PAHs in the ISM in different environments, affecting in turn the chemistry, the energy balance, and the contribution of PAHs to the extinction and emission curves. We studied the interaction of PAH cations with VUV photons in the 7-20 eV range from the synchrotron SOLEIL beamline, DESIRS. We recorded by action spectroscopy the relative intensities of photo-fragmentation and photo-ionization for a set of eight PAH cations ranging in size from 14 to 24 carbon atoms, with different structures. At photon energies below ~13.6 eV fragmentation dominates for the smaller species, while for larger species ionization is immediately competitive after the second ionization potential (IP). At higher photon energies, all species behave similarly, the ionization yield gradually increases, leveling off between 0.8 and 0.9 at ~18 eV. Among isomers, PAH structure appears to mainly affect the fragmentation cross section, but not the ionization cross section. We also measured the second IP for all species and the third IP for two of them, all are in good agreement with theoretical ones confirming that PAH cations can be further ionized in the diffuse ISM. Determining actual PAH dication abundances in the ISM will require detailed modeling. Our measured photo-ionization yields for several PAH cations provide a necessary ingredient for such models.

Linear carbon chains are common in various types of astronomical molecular sources. Possible formation mechanisms involve both bottom-up and top-down routes. We have carried out a combined observational and modeling study of the formation of carbon chains in the C-star envelope IRC+10216, where the polymerization of acetylene and hydrogen cyanide induced by ultraviolet photons can drive the formation of linear carbon chains of increasing length. We have used ALMA to map the emission of 3 mm rotational lines of the hydrocarbon radicals C2H, C4H, and C6H, and the CN-containing species CN, C3N, HC3N, and HC5N with an angular resolution of 1". The spatial distribution of all these species is a hollow, 5-10" wide, spherical shell located at a radius of 10-20" from the star, with no appreciable emission close to the star. Our observations resolve the broad shell of carbon chains into thinner sub-shells which are 1-2" wide and not fully concentric, indicating that the mass loss process has been discontinuous and not fully isotropic. The radial distributions of the species mapped reveal subtle differences: while the hydrocarbon radicals have very similar radial distributions, the CN-containing species show more diverse distributions, with HC3N appearing earlier in the expansion and the radical CN extending later than the rest of the species. The observed morphology can be rationalized by a chemical model in which the growth of polyynes is mainly produced by rapid gas-phase chemical reactions of C2H and C4H radicals with unsaturated hydrocarbons, while cyanopolyynes are mainly formed from polyynes in gas-phase reactions with CN and C3N radicals.

Protoplanetary disks undergo substantial mass-loss by photoevaporation, a mechanism which is crucial to their dynamical evolution. However, the processes regulating the gas energetics have not been well constrained by observations so far. We aim at studying the processes involved in disk photoevaporation when it is driven by far-UV photons. We present a unique Herschel survey and new ALMA observations of four externally-illuminated photoevaporating disks (a.k.a. proplyds). For the analysis of these data, we developed a 1D model of the photodissociation region (PDR) of a proplyd, based on the Meudon PDR code and computed the far infrared line emission. We successfully reproduce most of the observations and derive key physical parameters, i.e. densities at the disk surface of about $10^{6}$ cm$^{-3}$ and local gas temperatures of about 1000 K. Our modelling suggests that all studied disks are found in a transitional regime resulting from the interplay between several heating and cooling processes that we identify. These differ from those dominating in classical PDRs, i.e. grain photo-electric effect and cooling by [OI] and [CII] FIR lines. This energetic regime is associated to an equilibrium dynamical point of the photoevaporation flow: the mass-loss rate is self-regulated to set the envelope column density at a value that maintains the temperature at the disk surface around 1000 K. From our best-fit models, we estimate mass-loss rates - of the order of $10^{-7}$ $\mathrm{M}_\odot$/yr - that are in agreement with earlier spectroscopic observation of ionised gas tracers. This holds only if we assume an evaporation flow launched from the disk surface at sound speed (supercritical regime). We have identified the energetic regime regulating FUV-photoevaporation in proplyds. This regime could be implemented into models of the dynamical evolution of protoplanetary disks.

Photon Dominated Regions (PDRs) are interfaces between the mainly ionized and mainly molecular material around young massive stars. Analysis of the physical and chemical structure of such regions traces the impact of far-ultraviolet radiation of young massive stars on their environment. We present results on the physical and chemical structure of the prototypical high UV-illumination edge-on Orion Bar PDR from an unbiased spectral line survey with a wide spectral coverage. A spectral scan from 480-1250 GHz and 1410-1910 GHz at 1.1 MHz resolution was obtained by the HIFI instrument onboard the Herschel Space Observatory. For molecules with multiple transitions we used rotational diagrams to obtain excitation temperatures and column densities. For species with a single detected transition we used an optically thin LTE approximation. In case of species with available collisional rates, we also performed a non-LTE analysis to obtain kinetic temperatures, H2 volume densities, and column densities. About 120 lines corresponding to 29 molecules (including isotopologues) have been detected in the Herschel/HIFI line survey, including 11 transitions of CO, 7 transitions of 13CO, 6 transitions of C18O, 10 transitions of H2CO, and 6 transitions of H2O. Most species trace kinetic temperatures in the range between 100 and 150 K and H2 volume densities in the range between 10^5 and 10^6 cm^-3. The species with temperatures and / or densities outside of this range include the H2CO transitions tracing a very high temperature (315 K) and density (1.4x10^6 cm^-3) component and SO corresponding to the lowest temperature (56 K) measured as a part of this line survey. The observed lines/species reveal a range of physical conditions (gas density /temperature) involving structures at high density / high pressure, obsoleting the traditional 'clump / interclump' picture of the Orion Bar.

The abundance of CH+ and OH and excitation are predicted to be enhanced by the presence of vibrationally excited H2 or hot gas (~500-1000 K) in PDRs with high incident FUV radiation field. The excitation may also originate in dense gas (>10^5 cm-3) followed by nonreactive collisions. Previous observations suggest that the CH+ and OH correlate with dense and warm gas, and formation pumping contributes to CH+ excitation. We examine the spatial distribution of the CH+ and OH emission in the Orion Bar to establish their physical origin and main formation and excitation mechanisms. We present spatially sampled maps of the CH+ J=3-2 transition at 119.8 {\mu}m and the OH {\Lambda}-doublet at 84 {\mu}m in the Orion Bar over an area of 110"x110" with Herschel (PACS). We compare the spatial distribution of these molecules with those of their chemical precursors, C+, O and H2, and tracers of warm and dense gas. We assess the spatial variation of CH+ J=2-1 velocity-resolved line profile observed with Herschel (HIFI). The OH and CH+ lines correlate well with the high-J CO emission and delineate the warm and dense molecular region. While similar, the differences in the CH+ and OH morphologies indicate that CH+ formation and excitation are related to the observed vibrationally excited H2. This indicates that formation pumping contributes to the excitation of CH+. Interestingly, the peak of the rotationally excited OH 84 {\mu}m emission coincides with a bright young object, proplyd 244-440, which shows that OH can be an excellent tracer of UV-irradiated dense gas. The spatial distribution of CH+ and OH revealed in our maps is consistent with previous modeling studies. Both formation pumping and nonreactive collisions in a UV-irradiated dense gas are important CH+ J=3-2 excitation processes. The excitation of the OH {\Lambda}-doublet at 84 {\mu}m is mainly sensitive to the temperature and density.

The Orion Bar is the archetypal edge-on molecular cloud surface illuminated by strong ultraviolet radiation from nearby massive stars. Owing to the close distance to Orion (about 1,350 light-year), the effects of stellar feedback on the parental cloud can be studied in detail. Visible-light observations of the Bar(1) show that the transition between the hot ionised gas and the warm neutral atomic gas (the ionisation front) is spatially well separated from the transition from atomic to molecular gas (the dissociation front): about 15 arcseconds or 6,200 astronomical units (one astronomical unit is the Earth-Sun distance). Static equilibrium models(2,3) used to interpret previous far-infrared and radio observations of the neutral gas in the Bar(4,5,6) (typically at 10-20 arcsecond resolution) predict an inhomogeneous cloud structure consisting of dense clumps embedded in a lower density extended gas component. Here we report one-arcsecond-resolution millimetre-wave images that allow us to resolve the molecular cloud surface. In contrast to stationary model predictions(7,8,9), there is no appreciable offset between the peak of the H2 vibrational emission (delineating the H/H2 transition) and the edge of the observed CO and HCO+ emission. This implies that the H/H2 and C+/C/CO transition zones are very close. These observations reveal a fragmented ridge of high-density substructures, photoablative gas flows and instabilities at the molecular cloud surface. The results suggest that the cloud edge has been compressed by a high-pressure wave that currently moves into the molecular cloud. The images demonstrate that dynamical and nonequilibrium effects are important for the cloud evolution.

[Abridged] Fullerenes have been recently detected in various circumstellar and interstellar environments, raising the question of their formation pathway. It has been proposed that they can form by the photo-chemical processing of large polycyclic aromatic hydrocarbons (PAHs). Following our previous work on the evolution of PAHs in the NGC 7023 reflection nebula, we evaluate, using photochemical modeling, the possibility that the PAH C$_{66}$H$_{20}$ (i.e. circumovalene) can lead to the formation of C$_{60}$ upon irradiation by ultraviolet photons. The chemical pathway involves full dehydrogenation, folding into a floppy closed cage and shrinking of the cage by loss of C$_2$ units until it reaches the symmetric C$_{60}$ molecule. At 10" from the illuminating star and with realistic molecular parameters, the model predicts that 100% of C$_{66}$H$_{20}$ is converted into C$_{60}$ in $\sim$ 10$^5$ years, a timescale comparable to the age of the nebula. Shrinking appears to be the kinetically limiting step of the whole process. Hence, PAHs larger than C$_{66}$H$_{20}$ are unlikely to contribute significantly to the formation of C$_{60}$, while PAHs containing between 60 and 66 C atoms should contribute to the formation of C$_{60}$ with shorter timescales, and PAHs containing less than 60 C atoms will be destroyed. Assuming a classical size distribution for the PAH precursors, our model predicts absolute abundances of C$_{60}$ are up to several $10^{-4}$ of the elemental carbon, i.e. less than a percent of the typical interstellar PAH abundance, which is consistent with observational studies. According to our model, once formed, C$_{60}$ can survive much longer than other fullerenes because of the remarkable stability of the C$_{60}$ molecule at high internal energies.Hence, a natural consequence is that C$_{60}$ is more abundant than other fullerenes in highly irradiated environments.

A gas-phase optical spectrum of a large polycyclic aromatic hydrocarbon (PAH) cation - C78H26 +- in the 410-610 nm range is presented. This large all-benzenoid PAH should be large enough to be stable with respect to photodissociation in the harsh conditions prevailing in the interstellar medium (ISM). The spectrum is obtained via multi-photon dissociation (MPD) spectroscopy of cationic C78H26 stored in the Fourier Transform Ion Cyclotron Resonance (FT-ICR) cell using the radiation from a mid-band optical parametric oscillator (OPO) laser. The experimental spectrum shows two main absorption peaks at 431 nm and 516 nm, in good agreement with a theoretical spectrum computed via time-dependent density functional theory (TD-DFT). DFT calculations indicate that the equilibrium geometry, with the absolute minimum energy, is of lowered, nonplanar C2 symmetry instead of the more symmetric planar D2h symmetry that is usually the minimum for similar PAHs of smaller size. This kind of slightly broken symmetry could produce some of the fine structure observed in some diffuse interstellar bands (DIBs). It can also favor the folding of C78H26 + fragments and ultimately theformation of fullerenes. This study opens up the possibility to identify the most promising candidates for DIBs amongst large cationic PAHs.

Evolved stars are primary sources for the formation of polycyclic aromatic hydrocarbons (PAHs) and dust grains. Their circumstellar chemistry is usually designated as either oxygen-rich or carbon-rich, although dual-dust chemistry objects, whose infrared spectra reveal both silicate- and carbon-dust features, are also known. The exact origin and nature of this dual-dust chemistry is not yet understood. Spitzer-IRS mid-infrared spectroscopic imaging of the nearby, oxygen-rich planetary nebula NGC6720 reveals the presence of the 11.3 micron aromatic (PAH) emission band. It is attributed to emission from neutral PAHs, since no band is observed in the 7 to 8 micron range. The spatial distribution of PAHs is found to closely follow that of the warm clumpy molecular hydrogen emission. Emission from both neutral PAHs and warm H2 is likely to arise from photo-dissociation regions associated with dense knots that are located within the main ring. The presence of PAHs together with the previously derived high abundance of free carbon (relative to CO) suggest that the local conditions in an oxygen-rich environment can also become conducive to in-situ formation of large carbonaceous molecules, such as PAHs, via a bottom-up chemical pathway. In this scenario, the same stellar source can enrich the interstellar medium with both oxygen-rich dust and large carbonaceous molecules.

In 1985, "During experiments aimed at understanding the mechanisms by which long-chain carbon molecules are formed in interstellar space and circumstellar shells", Harry Kroto and his collaborators serendipitously discovered a new form of carbon: fullerenes. The most emblematic fullerene (i.e. C$_{60}$ "buckminsterfullerene"), contains exactly 60 carbon atoms organized in a cage-like structure similar to a soccer ball. Since their discovery impacted the field of nanotechnologies, Kroto and colleagues received the Nobel prize in 1996. The cage-like structure, common to all fullerene molecules, gives them unique properties, in particular an extraordinary stability. For this reason and since they were discovered in experiments aimed to reproduce conditions in space, fullerenes were sought after by astronomers for over two decades, and it is only recently that they have been firmly identified by spectroscopy, in evolved stars and in the interstellar medium. This identification offers the opportunity to study the molecular physics of fullerenes in the unique physical conditions provided by space, and to make the link with other large carbonaceous molecules thought to be present in space : polycyclic aromatic hydrocarbons.

The near- and mid-IR spectrum of many astronomical objects is dominated by emission bands due to UV-excited polycyclic aromatic hydrocarbons (PAH) and evaporating very small grains (eVSG). Previous studies with the ISO, Spitzer and AKARI space telescopes have shown that the spectral variations of these features are directly related to the local physical conditions that induce a photo-chemical evolution of the band carriers. Because of the limited sensitivity and spatial resolution, these studies have focused mainly on galactic star-forming regions. We discuss how the advent of JWST will allow to extend these studies to previously unresolved sources such as near-by galaxies, and how the analysis of the infrared signatures of PAHs and eVSGs can be used to determine their physical conditions and chemical composition.

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.

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.

In photon-dominated regions (PDRs), UV photons from nearby stars lead to the evaporation of very small grains (VSGs) and the production of gas-phase polycyclic aromatic hydrocarbons (PAHs). Our goal is to achieve better insight into the composition and evolution of evaporating very small grains (eVSGs) and PAHs through analyzing the infrared (IR) aliphatic and aromatic emission bands. We combined spectro-imagery in the near- and mid-IR to study the spatial evolution of the emission bands in the prototypical PDR NGC 7023. We used near-IR spectra obtained with AKARI to trace the evolution of the 3.3$\mu$m and 3.4$\mu$m bands, which are associated with aromatic and aliphatic C-H bonds on PAHs. The spectral fitting involves an additional broad feature centred at 3.45$\mu$m. Mid-IR observations obtained with Spitzer are used to discriminate the signatures of eVSGs, neutral and cationic PAHs. We correlated the spatial evolution of all these bands with the intensity of the UV field to explore the processing of their carriers. The intensity of the 3.45$\mu$m plateau shows an excellent correlation with that of the 3.3$\mu$m aromatic band (correlation coefficient R = 0.95), indicating that the plateau is dominated by the emission from aromatic bonds. The ratio of the 3.4$\mu$m and 3.3$\mu$m band intensity ($I_{3.4}/I_{3.3}$) decreases by a factor of 4 at the PDR interface from the more UV-shielded to the more exposed layers. The transition region between the aliphatic and aromatic material is found to correspond spatially with the transition zone between neutral PAHs and eVSGs. We conclude that the photo-processing of eVSGs leads to the production of PAHs with attached aliphatic sidegroups that are revealed by the 3.4$\mu$m emission band. Our analysis provides evidence for the presence of very small grains of mixed aromatic and aliphatic composition in PDRs.

We study the spatial distribution and chemistry of small hydrocarbons in the Orion Bar PDR. We used the IRAM-30m telescope to carry out a millimetre line survey towards the Orion Bar edge, complemented with ~2'x2' maps of the C2H and c-C3H2 emission. We analyse the excitation of the detected hydrocarbons and constrain the physical conditions of the emitting regions with non-LTE radiative transfer models. We compare the inferred column densities with updated gas-phase photochemical models including 13CCH and C13CH isotopomer fractionation. ~40% of the lines in the survey arise from hydrocarbons (C2H, C4H, c-C3H2, c-C3H, C13CH, 13CCH, l-C3H and l-H2C3). We detect new lines from l-C3H+ and improve its rotational spectroscopic constants. Anions or deuterated hydrocarbons are not detected: [C2D]/[C2H]<0.2%, [C2H-]/[C2H]<0.007% and [C4H-]/[C4H]<0.05%. Our gas-phase models can reasonably match the observed column densities of most hydrocarbons (within factors <3). Since the observed spatial distribution of the C2H and c-C3H2 emission is similar but does not follow the PAH emission, we conclude that, in high UV-flux PDRs, photodestruction of PAHs is not a necessary requirement to explain the observed abundances of the smallest hydrocarbons. Instead, gas-phase endothermic reactions (or with barriers) between C+, radicals and H2 enhance the formation of simple hydrocarbons. Observations and models suggest that the [C2H]/[c-C3H2] ratio (~32 at the PDR edge) decreases with the UV field attenuation. The observed low cyclic-to-linear C3H column density ratio (<3) is consistent with a high electron abundance (Xe) PDR environment. In fact, the poorly constrained Xe gradient influences much of the hydrocarbon chemistry in the more UV-shielded gas. We propose that reactions of C2H isotopologues with 13C+ and H atoms can explain the observed [C13CH]/[13CCH]=1.4(0.1) fractionation level.

The determination of the physical conditions in molecular clouds is a key step towards our understanding of their formation and evolution of associated star formation. We investigate the density, temperature, and column density of both dust and gas in the photodissociation regions (PDRs) located at the interface between the atomic and cold molecular gas of the NGC 7023 reflection nebula. We study how young stars affect the gas and dust in their environment. Our approach combining both dust and gas delivers strong constraints on the physical conditions of the PDRs. We find dense and warm molecular gas of high column density in the PDRs.

Clusters of polycyclic aromatic hydrocarbons (PAHs) have been proposed as candidates for evaporating very small grains, which are thought to be precursors of free-flying PAHs. Evaporation rates have been calculated so far only for species containing up to a few 100 C atoms, whereas interstellar PAH clusters could contain up to ~1000 C atoms. We present a method that generalises the calculation of the statistical evaporation rate of large PAH clusters and provides rates for species containing up to ~1000 C-atoms. The evaporation of non-rotating neutral homo-molecular PAH clusters containing up to 12 molecules from a family of highly symmetric compact PAHs is studied. Statistical calculations were performed and completed with molecular dynamics simulations at high internal energies to provide absolute values for the evaporation rate and distributions of kinetic energy released. The calculations used explicit atom-atom Lennard-Jones potentials in the rigid molecule approximation. A new method is proposed to take both inter- and intra-molecular vibrations into account. Without any parameter adjustment, the calculated evaporation rates agree well with available experimental data. We find that the non-rotation assumption has a limited impact on the evaporation rates. The photostability of PAH clusters increases dramatically with the size of molecules in the clusters, and to a lesser extent with the number of molecules in the clusters. For values of the UV radiation field that are typical of the regions where evaporating very small grains are observed, the smallest clusters in this study (~50 C-atoms) are found to be quickly photo-evaporated, whereas the largest clusters (~1000 C-atoms) are photostable. Our results support the idea that large PAH clusters are good candidates for evaporating very small grains.

Context. Monoceros R2 (Mon R2), at a distance of 830 pc, is the only ultra-compact H ii region (UC H ii) where its associated photon-dominated region (PDR) can be resolved with the Herschel Space Observatory. Aims. Our aim is to investigate observationally the kinematical patterns in the interface regions (i.e., the transition from atomic to molecular gas) associated with Mon R2. Methods. We used the HIFI instrument onboard Herschel to observe the line profiles of the reactive ions CH+, OH+ and H2O+ toward different positions in Mon R2. We derive the column density of these molecules and compare them with gas-phase chemistry models. Results. The reactive ion CH+ is detected both in emission (at central and red-shifted velocities) and in absorption (at blue-shifted velocities). OH+ is detected in absorption at both blue- and red-shifted velocities, with similar column densities. H2O+ is not detected at any of the positions, down to a rms of 40 mK toward the molecular peak. At this position, we find that the OH+ absorption originates in a mainly atomic medium, and therefore is associated with the most exposed layers of the PDR. These results are consistent with the predictions from photo-chemical models. The line profiles are consistent with the atomic gas being entrained in the ionized gas flow along the walls of the cavity of the H ii region. Based on this evidence, we are able to propose a new geometrical model for this region. Conclusions. The kinematical patterns of the OH+ and CH+ absorption indicate the existence of a layer of mainly atomic gas for which we have derived, for the first time, some physical parameters and its dynamics.

We report the detection in IRC+10216 of lines of HNC $J$=3-2 pertaining to 9 excited vibrational states with energies up to $\sim$5300 K. The spectrum, observed with ALMA, also shows a surprising large number of narrow, unidentified lines that arise in the vicinity of the star. The HNC data are interpreted through a 1D--spherical non--local radiative transfer model, coupled to a chemical model that includes chemistry at thermochemical equilibrium for the innermost regions and reaction kinetics for the external envelope. Although unresolved by the present early ALMA data, the radius inferred for the emitting region is $\sim$0.06'' (i.e., $\simeq$ 3 stellar radii), similar to the size of the dusty clumps reported by IR studies of the innermost region ($r <$ 0.3''). The derived abundance of HNC relative to H$_2$ is $10^{-8} <$ $\chi$(HNC) $< 10^{-6}$, and drops quickly where the gas density decreases and the gas chemistry is dominated by reaction kinetics. Merging HNC data with that of molecular species present throughout the inner envelope, such as vibrationally excited HCN, SiS, CS, or SiO, should allow us to characterize the physical and chemical conditions in the dust formation zone.

The far-IR range is a critical wavelength range to characterize the physical and chemical processes that transform the interstellar material into stars and planets. Objects in the earliest phases of stellar and planet evolution release most of their energy at these long wavelengths. In this contribution we briefly summarise some of the most relevant scientific advances achieved by the Herschel Space Observatory in the field. We also anticipate those that will be made possible by the large increase in sensitivity of SPICA cooled telescope. It is concluded that only through sensitive far-IR observations much beyond Herschel capabilities we will be able to constrain the mass, the energy budget and the water content of hundreds of protostars and planet-forming disks.

In the early 90s, C60+ was proposed as the carrier of two diffuse interstellar bands (DIBs) at 957.7 and 963.2 nm, but a firm identification still awaits gas-phase spectroscopic data. Neutral C60, on the other hand, was recently detected through its infrared emission bands in the interstellar medium and evolved stars. In this contribution, we present the detection of C60+ through its infrared vibrational bands in the NGC 7023 nebula, based on spectroscopic observations with the Spitzer space telescope, quantum chemistry calculation, and laboratory data from the literature. This detection supports the idea that C60+ could be a DIB carrier, and provides robust evidence that fullerenes exist in the gas-phase in the interstellar medium. Modeling efforts to design specific observations, combined with new gas-phase data, will be essential to confirm this proposal. A definitive attribution of the 957.7 and 963.2 nm DIBs to C60+ would represent a significant step forward in the field.

We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact HII region Mon R2. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate their chemistry. We present an observational study of CH, CCH and c-C$_3$H$_2$ in Mon R2 combining data obtained with the IRAM 30m telescope and Herschel. We determine the column densities of these species, and compare their spatial distributions with that of polycyclic aromatic hydrocarbon (PAH). We compare the observational results with different chemical models to explore the relative importance of gas-phase, grain-surface and time-dependent chemistry in these environments. The emission of the small hydrocarbons show different patterns. The CCH emission is extended while CH and c-C$_3$H$_2$ are concentrated towards the more illuminated layers of the PDR. The ratio of the column densities of c-C$_3$H$_2$ and CCH shows spatial variations up to a factor of a few, increasing from $N(c-C$_3$H$_2$)/N(CCH)\approx0.004$ in the envelope to a maximum of $\sim0.015-0.029$ towards the 8$\mu$m emission peak. Comparing these results with other galactic PDRs, we find that the abundance of CCH is quite constant over a wide range of G$_0$, whereas the abundance of c-C$_3$H$_2$ is higher in low-UV PDRs. In Mon R2, the gas-phase steady-state chemistry can account relatively well for the abundances of CH and CCH in the most exposed layers of the PDR, but falls short by a factor of 10 to reproduce c-C$_3$H$_2$. In the molecular envelope, time-dependent effects and grain surface chemistry play a dominant role in determining the hydrocarbons abundances. Our study shows that CCH and c-C$_3$H$_2$ present a complex chemistry in which UV photons, grain-surface chemistry and time dependent effects contribute to determine their abundances.

We observationally investigate the relation between the photoelectric heating efficiency in PDRs and the charge of PAHs, which are considered to play a key role in photoelectric heating. Using PACS onboard Herschel, we observed six PDRs spanning a wide range of FUV radiation fields (G_0=100-10^5). To measure the photoelectric heating efficiency, we obtained the intensities of the main cooling lines, i.e., the [OI]63um, 145um, and [CII]158um, as well as the FIR continuum intensity. We used Spitzer/IRS spectroscopic mapping observations to investigate the MIR PAH features in the same regions. We decomposed the MIR PAH emission into that of neutral (PAH^0) and positively ionized (PAH^+) species to derive the fraction of the positively charged PAHs, and compare it to the photoelectric heating efficiency. The heating efficiency traced by ([OI]63um+[OI]145um+[CII]158um) / TIR, ranges between 0.1% and 0.9% in different sources, and the fraction of PAH^+ relative to (PAH^0 + PAH^+) spans from 0(+11)% to 87(+/-10)%. All positions with a high PAH^+ fraction show a low heating efficiency, and all positions with a high heating efficiency have a low PAH^+ fraction, supporting the scenario in which a positive grain charge results in a decreased heating efficiency. Theoretical estimates of the photoelectric heating efficiency show a stronger dependence on the charging parameter gamma=G_0 T^{1/2}/n_e than the observed efficiency reported in this study, and the discrepancy is significant at low gamma. The photoelectric heating efficiency on PAHs, traced by ([OI]63um+[OI]145um+[CII]158um) / (PAH+[OI]63um+[OI]145um+[CII]158um), shows a much better match between the observations and the theoretical estimates. The good agreement of the photoelectric heating efficiency on PAHs with a theoretical model indicates the dominant contribution of PAHs to the photoelectric heating. (abridged for arXiv)

Various studies have emphasized variations of the charge state and composition of the interstellar polycyclic aromatic hydrocarbon (PAH) population in photodissociation regions (PDRs). We aim to model the spatial evolution of the charge and hydrogenation states of PAHs in PDRs. We focus on the specific case of the north-west (NW) PDR of NGC 7023 and also discuss the case of the diffuse interstellar medium (ISM). The physical conditions in NGC 7023 NW are modelled using a state-of-the-art PDR code. We then use a new PAH chemical evolution model that includes recent experimental data on PAHs and describes multiphoton events. We consider a family of compact PAHs bearing up to 96 carbon atoms. The calculated ionization ratio is in good agreement with observations in NGC 7023 NW. Within the PDR, PAHs evolve into three major populations: medium-sized PAHs (50<Nc<90) are normally hydrogenated, larger PAHs (Nc>90) can be superhydrogenated, and smaller species (Nc<50) are fully dehydrogenated. In the cavity, where the fullerene C60 was recently detected, all the studied PAHs are found to be quickly fully dehydrogenated. PAH chemical evolution exhibits a complex non-linear behaviour as a function of the UV radiation field because of multiphoton events. Steady state for hydrogenation is reached on timescales ranging from less than a year for small PAHs, up to 10000 years for large PAHs at Av=1. We identified critical reactions that need more studies. Our new model allows us to rationalize the observational constraints without any fitting parameter. PAHs smaller than 50 carbon atoms are not expected to survive in the NGC 7023 NW PDR. A similar conclusion is obtained for the diffuse ISM. Carbon clusters turn out to be end products of PAH photodissociation, and the evolution of these clusters needs to be investigated further to evaluate their impact on the chemical and physical evolution of PDRs.

The abundances of interstellar CH+ and SH+ are not well understood as their most likely formation channels are highly endothermic. Using data from Herschel, we study the formation of CH+ and SH+ in a typical high UV-illumination photon-dominated region (PDR), the Orion Bar. Herschel/HIFI provides velocity-resolved data of CH+ 1-0 and 2-1 and three hyperfine transitions of SH+. Herschel/PACS provides information on the excitation and spatial distribution of CH+ (up to J=6-5). The widths of the CH+ 2-1 and 1-0 transitions are of ~5 km/s, significantly broader than the typical width of dense gas tracers in the Orion Bar (2-3 km/s) and are comparable to the width of tracers of the interclump medium such as C+ and HF. The detected SH+ transitions are narrower compared to CH+ and have line widths of 3 km/s, indicating that SH+ emission mainly originates in denser condensations. Non-LTE radiative transfer models show that electron collisions affect the excitation of CH+ and SH+, and that reactive collisions need to be taken into account to calculate the excitation of CH+. Comparison to PDR models shows that CH+ and SH+ are tracers of the warm surface region (AV<1.5) of the PDR with temperatures between 500-1000 K. We have also detected the 5-4 transition of CF+ (FWHM=1.9 km/s) with an intensity that is consistent with previous observations of lower-J CF+ transitions toward the Orion Bar. A comparison to PDR models indicate that the internal vibrational energy of H2 can explain the formation of CH+ for typical physical conditions in the Orion Bar near the ionization front. H2 vibrational excitation is the most likely explanation of SH+ formation as well. The abundance ratios of CH+ and SH+ trace the destruction paths of these ions, and through that, indirectly, the ratios of H, H2 and electron abundances as a function of depth into the cloud.

Buckminsterfullerene (C60) was recently detected through its infrared emission bands in the interstellar medium (ISM), including in the proximity of massive stars, where physical conditions could favor the formation of the cationic form, C60+. In addition, C60+ was proposed as the carrier of two diffuse interstellar bands in the near-IR, although a firm identification still awaits for gas-phase spectroscopic data. We examined in details the Spitzer IRS spectra of the NGC 7023 reflection nebula, at a position close (7.5") to the illuminating B star HD 200775, and found four previously unreported bands at 6.4, 7.1, 8.2 and 10.5 \mu m in addition to the classical bands attributed to Polycylic Aromatic Hydrocarbons (PAHs) and neutral C60. These 4 bands are observed only in this region of the nebula, while C60 emission is still present slightly further away from the star, and PAH emission even further away. Based on this observation, on theoretical calculations we perform, and on laboratory studies, we attribute these bands to C60+. The detection of C60+ confirms the idea that large carbon molecules exist in the gas-phase in these environments. In addition, the relative variation of the C60, and C60+, band intensities constitutes a potentially powerful probe of the physical conditions in highly UV-irradiated regions.

Large spatial-spectral surveys are more and more common in astronomy. This calls for the need of new methods to analyze such mega- to giga-pixel data-cubes. In this paper we present a method to decompose such observations into a limited and comprehensive set of components. The original data can then be interpreted in terms of linear combinations of these components. The method uses non-negative matrix factorization (NMF) to extract latent spectral end-members in the data. The number of needed end-members is estimated based on the level of noise in the data. A Monte-Carlo scheme is adopted to estimate the optimal end-members, and their standard deviations. Finally, the maps of linear coefficients are reconstructed using non-negative least squares. We apply this method to a set of hyperspectral data of the NGC 7023 nebula, obtained recently with the HIFI instrument onboard the Herschel space observatory, and provide a first interpretation of the results in terms of 3-dimensional dynamical structure of the region.

(Abridged) We have observed velocity resolved spectra of four ro-vibrational far-infrared transitions of C3 between the vibrational ground state and the low-energy nu2 bending mode at frequencies between 1654--1897 GHz using HIFI on board Herschel, in DR21(OH), a high mass star forming region. Several transitions of CCH and c-C3H2 have also been observed with HIFI and the IRAM 30m telescope. A gas and grain warm-up model was used to identify the primary C3 forming reactions in DR21(OH). We have detected C3 in absorption in four far-infrared transitions, P(4), P(10), Q(2) and Q(4). The continuum sources MM1 and MM2 in DR21(OH) though spatially unresolved, are sufficiently separated in velocity to be identified in the C3 spectra. All C3 transitions are detected from the embedded source MM2 and the surrounding envelope, whereas only Q(4) & P(4) are detected toward the hot core MM1. The abundance of C3 in the envelope and MM2 is \sim6x10^{-10} and \sim3x10^{-9} respectively. For CCH and c-C3H2 we only detect emission from the envelope and MM1. The observed CCH, C3, and c-C3H2 abundances are most consistent with a chemical model with n(H2)\sim5x10^{6} cm^-3 post-warm-up dust temperature, T_max =30 K and a time of \sim0.7-3 Myr. Post warm-up gas phase chemistry of CH4 released from the grain at t\sim 0.2 Myr and lasting for 1 Myr can explain the observed C3 abundance in the envelope of DR21(OH) and no mechanism involving photodestruction of PAH molecules is required. The chemistry in the envelope is similar to the warm carbon chain chemistry (WCCC) found in lukewarm corinos. The observed lower C3 abundance in MM1 as compared to MM2 and the envelope could be indicative of destruction of C3 in the more evolved MM1. The timescale for the chemistry derived for the envelope is consistent with the dynamical timescale of 2 Myr derived for DR21(OH) in other studies.

Context. Mon R2 is the only ultracompact HII region (UCHII) where the associated photon-dominated region (PDR) can be resolved with Herschel. Due to its brightness and proximity, it is the best source to investigate the chemistry and physics of highly UV-irradiated PDRs. Aims. Our goal is to estimate the abundance of H2O and NH3 in this region and investigate their origin. Methods. We present new observations obtained with HIFI and the IRAM-30m telescope. Using a large velocity gradient approach, we model the line intensities and derive an average abundance of H2O and NH3 across the region. Finally, we model the line profiles with a non-local radiative transfer model and compare these results with the abundance predicted by the Meudon PDR code. Results. The variations of the line profiles and intensities indicate complex geometrical and kinematical patterns. The H2O lines present a strong absorption at the ambient velocity and emission in high velocity wings towards the HII region. The spatial distribution of the o-H2^18O line shows that the its emission arises in the PDR surrounding the HII region. By modeling the o-H2^18O emission we derive a mean abundance of o-H2O of ~10^-8 relative to H2. The ortho-H2O abundance is however larger, ~1x10^-7, in the high velocity wings. Possible explanations for this larger abundance include an expanding hot PDR and/or an outflow. Ammonia seems to be present only in the envelope with an average abundance of ~2x10^-9 relative to H2. Conclusions. The Meudon PDR code can account for the measured water abundance in the high velocity gas as long as we assume that it originates from a <1 mag hot expanding layer of the PDR, i.e. that the outflow has only a minor contribution to this emission. To explain the abundances in the rest of the cloud the molecular freeze out and grain surface chemistry would need to be included.

Context. In photo-dissociation regions (PDRs), Polycyclic Aromatic Hydrocarbons (PAHs) could be produced by evaporation of Very Small Grains (VSGs) by the impinging UV radiation field from a nearby star. Aims. We investigate quantitatively the transition zone between evaporating Very Small Grains (eVSGs) and PAHs in several PDRs. Methods. We study the relative contribution of PAHs and eVSGs to the mid-IR emission in a wide range of excitation conditions. We fit the observed mid-IR emission of PDRs by using a set of template band emission spectra of PAHs, eVSGs and gas lines. The fitting tool PAHTAT (PAH Toulouse Astronomical Templates) is made available to the community as an IDL routine. From the results of the fit, we derive the fraction of carbon f_eVSG locked in eVSGs and compare it to the intensity of the local UV radiation field. Results. We show a clear decrease of f_eVSG with increasing intensity of the local UV radiation field, which supports the scenario of photo-destruction of eVSGs. Conversely, this dependence can be used to quantify the intensity of the UV radiation field for different PDRs, including non resolved ones. Conclusions. PAHTAT can be used to trace the intensity of the local UV radiation field in regions where eVSGs evaporate, which correspond to relatively dense (nH = [100, 10^5 ] cm-3) and UV irradiated PDRs (G0 = [100, 5x10^4]) where H2 emits in rotational lines.

Ultracompact (UC) HII regions constitute one of the earliest phases in the formation of a massive star and are characterized by extreme physical conditions (Go>10^5 Habing field and n>10^6 cm^-3). The UC HII Mon R2 is the closest one and therefore an excellent target to study the chemistry in these complex regions. We carried out a 3mm and 1mm spectral survey using the IRAM 30-m telescope towards three positions that represent different physical environments in Mon R2: (i) the ionization front (IF) at (0",0"); two peaks in the molecular cloud (ii) MP1 at the offset (+15",-15") and (iii) MP2 at the farther offset (0",40"). In addition, we carried out extensive modeling to explain the chemical differences between the three observed regions. We detected more than thirty different species. We detected SO+ and C4H suggesting that UV radiation plays an important role in the molecular chemistry of this region. We detected the typical PDR molecules CN, HCN, HCO, C2H, and c-C3H2. While the IF and the MP1 have a chemistry similar to that found in high UV field and dense PDRs like the Orion Bar, the MP2 is more similar to lower UV/density PDRs like the Horsehead nebula. We also detected complex molecules that are not usually found in PDRs (CH3CN, H2CO, HC3N, CH3OH and CH3C2H). Sulfur compounds CS, HCS+, C2S, H2CS, SO and SO2 and the deuterated species DCN and C2D were also identified. [DCN]/[HCN]=0.03 and [C2D]/[C2H]=0.05, are among the highest in warm regions. Our results show that the high UV/dense PDRs present a different chemistry from that of the low UV case. Abundance ratios like [CO+]/[HCO+] or [HCO]/[HCO+] are good diagnostics to differentiate between them. In Mon R2 we have the two classes of PDRs, a high UV PDR towards the IF and the adjacent molecular bar and a low-UV PDR which extends towards the north-west following the border of the cloud.

We present spatially resolved Herschel/PACS observations of the Orion Bar. We have characterise the emission of the far-infrared fine-structure lines of [CII] (158um), [OI] (63 and 145um), and [NII] (122um) that trace the gas local conditions. The observed distribution and variation of the lines are discussed in relation to the underlying geometry and linked to the energetics associated with the Trapezium stars. These observations enable us to map the spatial distribution of these fine-structure lines with a spatial resolution between 4" and 11" and covering a total square area of about 120"x105". The spatial profile of the emission lines are modelled using the radiative transfer code Cloudy. We find that the spatial distribution of the [CII] line coincides with that of the [OI] lines. The [NII] line peaks closer to the ionising star than the other three lines, but with a small region of overlap. We can distinguish several knots of enhanced emission within the Bar indicating the presence of an inhomogenous and structured medium. The emission profiles cannot be reproduced by a single photo-dissociation region, clearly indicating that, besides the Bar, there is a significant contribution from additional photo-dissociation region(s) over the area studied. The combination of both the [NII] and [OI] 145um lines can be used to estimate the [CII] emission and distinguish between its ionised or neutral origin. We have calculated how much [CII] emission comes from the neutral and ionised region, and find that at least 82% originates from the photo-dissocciation region. Together, the [CII] 158um and [OI] 63 and 145um lines account for 90% of the power emitted by the main cooling lines in the Bar (including CO, H2, etc...), with [OI] 63um alone accounting for 72% of the total.

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.

As part of a far-infrared (FIR) spectral scan with Herschel/PACS, we present the first detection of the hydroxyl radical (OH) towards the Orion Bar photodissociation region (PDR). Five OH rotational Lambda-doublets involving energy levels out to E_u/k~511 K have been detected (at ~65, ~79, ~84, ~119 and ~163um). The total intensity of the OH lines is I(OH)~5x10^-4 erg s^-1 cm^-2 sr^-1. The observed emission of rotationally excited OH lines is extended and correlates well with the high-J CO and CH^+ J=3-2 line emission (but apparently not with water vapour), pointing towards a common origin. Nonlocal, non-LTE radiative transfer models including excitation by the ambient FIR radiation field suggest that OH arises in a small filling factor component of warm (Tk~160-220 K) and dense (n_H~10^{6-7} cm^-3) gas with source-averaged OH column densities of ~10^15 cm^-2. High density and temperature photochemical models predict such enhanced OH columns at low depths (A_V<1) and small spatial scales (~10^15 cm), where OH formation is driven by gas-phase endothermic reactions of atomic oxygen with molecular hydrogen. We interpret the extended OH emission as coming from unresolved structures exposed to far-ultraviolet (FUV) radiation near the Bar edge (photoevaporating clumps or filaments) and not from the lower density "interclump" medium. Photodissociation leads to OH/H2O abundance ratios (>1) much higher than those expected in equally warm regions without enhanced FUV radiation fields.

We have examined methanol emission from Orion KL with of the {\em Herschel}/HIFI instrument, and detected two methanol bands centered at 524 GHz and 1061 GHz. The 524 GHz methanol band (observed in HIFI band 1a) is dominated by the isolated $\Delta$J$=$0, K$=-4\rightarrow$-3, v$_t$$=$0 Q branch, and includes 25 E-type and 2 A-type transitions. The 1061 GHz methanol band (observed in HIFI band 4b) is dominated by the $\Delta$J$=$0, K$=7\rightarrow$6, v$_t$$=$0 Q branch transitions which are mostly blended. We have used the isolated E-type v$_t$$=$0 methanol transitions to explore the physical conditions in the molecular gas. With HIFI's high velocity resolution, the methanol emission contributed by different spatial components along the line of sight toward Orion KL (hot core, low velocity flow, and compact ridge) can be distinguished and studied separately. The isolated transitions detected in these bands cover a broad energy range (upper state energy ranging from 80 K to 900 K), which provides a unique probe of the thermal structure in each spatial component. The observations further show that the compact ridge is externally heated. These observations demonstrate the power of methanol lines as probes of the physical conditions in warm regions in close proximity to young stars.

We present spectroscopic observations obtained with the infrared Spitzer Space Telescope, which provide insight into the H$_2$ physics and gas energetics in photodissociation Regions (PDRs) of low to moderate far-ultraviolet (FUV) fields and densities. We analyze data on six well known Galactic PDRs (L1721, California, N7023E, Horsehead, rho Oph, N2023N), sampling a poorly explored range of excitation conditions ($\chi \sim 5-10^3$), relevant to the bulk of molecular clouds in galaxies. Spitzer observations of H$_2$ rotational lines are complemented with H$_2$ data, including ro-vibrational line measurements, obtained with ground-based telescopes and ISO, to constrain the relative contributions of ultraviolet pumping and collisions to the H$_2$ excitation. The data analysis is supported by model calculations with the Meudon PDR code. The observed column densities of rotationally excited H$_2$ are observed to be much higher than PDR model predictions. In the lowest excitation PDRs, the discrepancy between the model and the data is about one order of magnitude for rotational levels $J \ge $3. We discuss whether an enhancement in the H$_2$ formation rate or a local increase in photoelectric heating, as proposed for brighter PDRs in former ISO studies, may improve the data-model comparison. We find that an enhancement in the H$_2$ formation rates reduces the discrepancy, but the models still fall short of the data. This large disagreement suggests that our understanding of the formation and excitation of H$_2$ and/or of PDRs energetics is still incomplete. We discuss several explanations, which could be further tested using the Herschel Space Telescope

We report observations of the reactive molecular ions OH$^+$, H$_2$O$^+$, and H$_3$O$^+$ towards Orion KL with Herschel/HIFI. All three $N=1-0$ fine-structure transitions of OH$^+$ at 909, 971, and 1033GHz and both fine-structure components of the doublet {\it ortho}-H$_2$O$^+$ $1_{11}-0_{00}$ transition at 1115 and 1139GHz were detected; an upper limit was obtained for H$_3$O$^+$. OH$^+$ and H$_2$O$^+$ are observed purely in absorption, showing a narrow component at the source velocity of 9 kms$^{-1}$, and a broad blueshifted absorption similar to that reported recently for HF and {\it para}-H$_{2}^{18}$O, and attributed to the low velocity outflow of Orion KL. We estimate column densities of OH$^+$ and H$_2$O$^+$ for the 9 km s$^{-1}$ component of $9 \pm 3 \times 10^{12}$cm$^{-2}$ and $7 \pm 2 \times 10^{12}$cm$^{-2}$, and those in the outflow of $1.9 \pm 0.7 \times 10^{13}$cm$^{-2}$ and $1.0 \pm 0.3 \times 10^{13}$cm$^{-2}$. Upper limits of $2.4\times 10^{12}$cm$^{-2}$ and $8.7\times 10^{12}$cm$^{-2}$ were derived for the column densities of {\it ortho} and {\it para}-H$_3$O$^+$ from transitions near 985 and 1657GHz. The column densities of the three ions are up to an order of magnitude lower than those obtained from recent observations of W31C and W49N. The comparatively low column densities may be explained by a higher gas density despite the assumption of a very high ionization rate.

The fullerene C$_{60}$ has four infrared-active vibrational transitions at 7.0, 8.5, 17.4 and 18.9 $\mu$m. We have previously observed emission features at 17.4 and 18.9 $\mu$m in the reflection nebula NGC 7023 and demonstrated spatial correlations suggestive of a common origin. We now confirm our earlier identification of these features with C$_{60}$ by detecting a third emission feature at 7.04 $\pm$ 0.05 $\mu$m in NGC 7023. We also report the detection of these three C$_{60}$ features in the reflection nebula NGC 2023. Our spectroscopic mapping of NGC 7023 shows that the 18.9 $\mu$m C$_{60}$ feature peaks on the central star and that the 16.4 $\mu$m emission feature due to polycyclic aromatic hydrocarbons peaks between the star and a nearby photodissociation front. The observed features in NGC 7023 are consistent with emission from UV-excited gas-phase C$_{60}$. We find that 0.1-0.6% of interstellar carbon is in C$_{60}$; this abundance is consistent with those from previous upper limits and possible fullerene detections in the interstellar medium. This is the first firm detection of neutral C$_{60}$ in the interstellar medium.

We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species.

We report the detection of more than 48 velocity-resolved ground rotational state transitions of H2(16)O, H2(18)O, and H2(17)O - most for the first time - in both emission and absorption toward Orion KL using Herschel/HIFI. We show that a simple fit, constrained to match the known emission and absorption components along the line of sight, is in excellent agreement with the spectral profiles of all the water lines. Using the measured H2(18)O line fluxes, which are less affected by line opacity than their H2(16)O counterparts, and an escape probability method, the column densities of H2(18)O associated with each emission component are derived. We infer total water abundances of 7.4E-5, 1.0E-5, and 1.6E-5 for the plateau, hot core, and extended warm gas, respectively. In the case of the plateau, this value is consistent with previous measures of the Orion-KL water abundance as well as those of other molecular outflows. In the case of the hot core and extended warm gas, these values are somewhat higher than water abundances derived for other quiescent clouds, suggesting that these regions are likely experiencing enhanced water-ice sublimation from (and reduced freeze-out onto) grain surfaces due to the warmer dust in these sources.

To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular cores, whose asymmetries trace infall and expansion motions. The high-mass star forming region SgrB2(M) was observed with Herschel/HIFI (HEXOS key project) in various lines of HCN and its isotopologues, complemented by APEX data. The observations are compared to spherically symmetric, centrally heated models with density power-law gradient and different velocity fields (infall or infall+expansion), using the radiative transfer code RATRAN. The HCN line profiles are asymmetric, with the emission peak shifting from blue to red with increasing J and decreasing line opacity (HCN to H$^{13}$CN). This is most evident in the HCN 12--11 line at 1062 GHz. These line shapes are reproduced by a model whose velocity field changes from infall in the outer part to expansion in the inner part. The qualitative reproduction of the HCN lines suggests that infall dominates in the colder, outer regions, but expansion dominates in the warmer, inner regions. We are thus witnessing the onset of feedback in massive star formation, starting to reverse the infall and finally disrupting the whole molecular cloud. To obtain our result, the THz lines uniquely covered by HIFI were critically important.

Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes -- grain surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr~B2(M) presented here allow the determination of the HDO abundance throughout the envelope, which has not been possible before with ground-based observations only. The abundance structure has been modeled with the spherical Monte Carlo radiative transfer code RATRAN, which also takes radiative pumping by continuum emission from dust into account. The modeling reveals that the abundance of HDO rises steeply with temperature from a low abundance ($2.5\times 10^{-11}$) in the outer envelope at temperatures below 100~K through a medium abundance ($1.5\times 10^{-9}$) in the inner envelope/outer core, at temperatures between 100 and 200~K, and finally a high abundance ($3.5\times 10^{-9}$) at temperatures above 200~K in the hot core.

Using Herschel's HIFI instrument we have observed [C II] along a cut through S140 and high-J transitions of CO and HCO+ at two positions on the cut, corresponding to the externally irradiated ionization front and the embedded massive star forming core IRS1. The HIFI data were combined with available ground-based observations and modeled using the KOSMA-tau model for photon dominated regions. Here we derive the physical conditions in S140 and in particular the origin of [C II] emission around IRS1. We identify three distinct regions of [C II] emission from the cut, one close to the embedded source IRS1, one associated with the ionization front and one further into the cloud. The line emission can be understood in terms of a clumpy model of photon-dominated regions. At the position of IRS1, we identify at least two distinct components contributing to the [C II] emission, one of them a small, hot component, which can possibly be identified with the irradiated outflow walls. This is consistent with the fact that the [C II] peak at IRS1 coincides with shocked H2 emission at the edges of the outflow cavity. We note that previously available observations of IRS1 can be well reproduced by a single-component KOSMA-tau model. Thus it is HIFI's unprecedented spatial and spectral resolution, as well as its sensitivity which has allowed us to uncover an additional hot gas component in the S140 region.

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.

The HIFI instrument on board the Herschel Space Observatory has been used to observe interstellar nitrogen hydrides along the sight-line towards G10.6-0.4 in order to improve our understanding of the interstellar chemistry of nitrogen. We report observations of absorption in NH N=1-0, J=2-1 and ortho-NH2 1_1,1-0_0,0. We also observed ortho-NH3 1_0-0_0, and 2_0-1_0, para-NH3 2_1-1_1, and searched unsuccessfully for NH+. All detections show emission and absorption associated directly with the hot-core source itself as well as absorption by foreground material over a wide range of velocities. All spectra show similar, non-saturated, absorption features, which we attribute to diffuse molecular gas. Total column densities over the velocity range 11-54 km/s are estimated. The similar profiles suggest fairly uniform abundances relative to hydrogen, approximately 6*10^-9, 3*10^-9, and 3*10^-9 for NH, NH2, and NH3, respectively. These abundances are discussed with reference to models of gas-phase and surface chemistry.

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 discuss the detection of absorption by interstellar hydrogen fluoride (HF) along the sight line to the submillimeter continuum sources W49N and W51. We have used Herschel's HIFI instrument in dual beam switch mode to observe the 1232.4762 GHz J = 1 - 0 HF transition in the upper sideband of the band 5a receiver. We detected foreground absorption by HF toward both sources over a wide range of velocities. Optically thin absorption components were detected on both sight lines, allowing us to measure - as opposed to obtain a lower limit on - the column density of HF for the first time. As in previous observations of HF toward the source G10.6-0.4, the derived HF column density is typically comparable to that of water vapor, even though the elemental abundance of oxygen is greater than that of fluorine by four orders of magnitude. We used the rather uncertain N(CH)-N(H2) relationship derived previously toward diffuse molecular clouds to infer the molecular hydrogen column density in the clouds exhibiting HF absorption. Within the uncertainties, we find that the abundance of HF with respect to H2 is consistent with the theoretical prediction that HF is the main reservoir of gas-phase fluorine for these clouds. Thus, hydrogen fluoride has the potential to become an excellent tracer of molecular hydrogen, and provides a sensitive probe of clouds of small H2 column density. Indeed, the observations of hydrogen fluoride reported here reveal the presence of a low column density diffuse molecular cloud along the W51 sight line, at an LSR velocity of ~ 24kms-1, that had not been identified in molecular absorption line studies prior to the launch of Herschel.

We report a detection of the fundamental rotational transition of hydrogen fluoride in absorption towards Orion KL using Herschel/HIFI. After the removal of contaminating features associated with common molecules ("weeds"), the HF spectrum shows a P-Cygni profile, with weak redshifted emission and strong blue-shifted absorption, associated with the low-velocity molecular outflow. We derive an estimate of 2.9 x 10^13 cm^-2 for the HF column density responsible for the broad absorption component. Using our best estimate of the H2 column density within the low-velocity molecular outflow, we obtain a lower limit of ~1.6 x 10^-10 for the HF abundance relative to hydrogen nuclei, corresponding to 0.6% of the solar abundance of fluorine. This value is close to that inferred from previous ISO observations of HF J=2--1 absorption towards Sgr B2, but is in sharp contrast to the lower limit of 6 x 10^-9 derived by Neufeld et al. (2010) for cold, foreground clouds on the line of sight towards G10.6-0.4.

We have observed CH absorption lines ($J=3/2, N=1 \leftarrow J=1/2, N=1$) against the continuum source Sgr~B2(M) using the \textit{Herschel}/HIFI instrument. With the high spectral resolution and wide velocity coverage provided by HIFI, 31 CH absorption features with different radial velocities and line widths are detected and identified. The narrower line width and lower column density clouds show `spiral arm' cloud characteristics, while the absorption component with the broadest line width and highest column density corresponds to the gas from the Sgr~B2 envelope. The observations show that each `spiral arm' harbors multiple velocity components, indicating that the clouds are not uniform and that they have internal structure. This line-of-sight through almost the entire Galaxy offers unique possibilities to study the basic chemistry of simple molecules in diffuse clouds, as a variety of different cloud classes are sampled simultaneously. We find that the linear relationship between CH and H$_2$ column densities found at lower $A_V$ by UV observations does not continue into the range of higher visual extinction. There, the curve flattens, which probably means that CH is depleted in the denser cores of these clouds.

Mon R2, at a distance of 830 pc, is the only ultracompact HII region (UC HII) where the photon-dominated region (PDR) between the ionized gas and the molecular cloud can be resolved with Herschel. HIFI observations of the abundant compounds 13CO, C18O, o-H2-18O, HCO+, CS, CH, and NH have been used to derive the physical and chemical conditions in the PDR, in particular the water abundance. The 13CO, C18O, o-H2-18O, HCO+ and CS observations are well described assuming that the emission is coming from a dense (n=5E6 cm-3, N(H2)>1E22 cm-2) layer of molecular gas around the UC HII. Based on our o-H2-18O observations, we estimate an o-H2O abundance of ~2E-8. This is the average ortho-water abundance in the PDR. Additional H2-18O and/or water lines are required to derive the water abundance profile. A lower density envelope (n~1E5 cm-3, N(H2)=2-5E22 cm-2) is responsible for the absorption in the NH 1_1-0_2 line. The emission of the CH ground state triplet is coming from both regions with a complex and self-absorbed profile in the main component. The radiative transfer modeling shows that the 13CO and HCO+ line profiles are consistent with an expansion of the molecular gas with a velocity law, v_e =0.5 x (r/Rout)^{-1} km/s, although the expansion velocity is poorly constrained by the observations presented here.

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 first detection of the ground-state rotational transition of the methylidyne cation CH+ towards the massive star-forming region DR21 with the HIFI instrument onboard the Herschel satellite. The line profile exhibits a broad emission line, in addition to two deep and broad absorption features associated with the DR21 molecular ridge and foreground gas. These observations allow us to determine a CH+ J=1-0 line frequency of 835137 +/- 3 MHz, in good agreement with a recent experimental determination. We estimate the CH+ column density to be a few 1e13 cm^-2 in the gas seen in emission, and > 1e14 cm^-2 in the components responsible for the absorption, which is indicative of a high line of sight average abundance [CH+]/[H] > 1.2x10^-8. We show that the CH+ column densities agree well with the predictions of state-of-the-art C-shock models in dense UV-illuminated gas for the emission line, and with those of turbulent dissipation models in diffuse gas for the absorption lines.

We present Herschel/HIFI observations of the fundamental rotational transitions of ortho- and para-H$_2^{16}$O and H$_2^{18}$O in absorption towards Sagittarius~B2(M) and W31C. The ortho/para ratio in water in the foreground clouds on the line of sight towards these bright continuum sources is generally consistent with the statistical high-temperature ratio of 3, within the observational uncertainties. However, somewhat unexpectedly, we derive a low ortho/para ratio of $2.35 \pm 0.35$, corresponding to a spin temperature of $\sim$27~K, towards Sagittarius~B2(M) at velocities of the expanding molecular ring. Water molecules in this region appear to have formed with, or relaxed to, an ortho/para ratio close to the value corresponding to the local temperature of the gas and dust.

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.

H2O+ has been observed in its ortho- and para- states toward the massive star forming core Sgr B2(M), located close to the Galactic center. The observations show absorption in all spiral arm clouds between the Sun and Sgr B2. The average o/p ratio of H2O+ in most velocity intervals is 4.8, which corresponds to a nuclear spin temperature of 21 K. The relationship of this spin temperature to the formation temperature and current physical temperature of the gas hosting H2O+ is discussed, but no firm conclusion is reached. In the velocity interval 0 to 60 km/s, an ortho/para ratio of below unity is found, but if this is due to an artifact of contamination by other species or real is not clear.

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.

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 a power spectrum analysis of the Herschel-SPIRE observations of the Polaris flare, a high Galactic latitude cirrus cloud midway between the diffuse and molecular phases. The SPIRE images of the Polaris flare reveal for the first time the structure of the diffuse interstellar medium down to 0.01 parsec over a 10 square degrees region. These exceptional observations highlight the highly filamentary and clumpy structure of the interstellar medium even in diffuse regions of the map. The power spectrum analysis shows that the structure of the interstellar medium is well described by a single power law with an exponent of -2.7 +- 0.1 at all scales from 30" to 8 degrees. That the power spectrum slope of the dust emission is constant down to the SPIRE angular resolution is an indication that the inertial range of turbulence extends down to the 0.01 pc scale. The power spectrum analysis also allows the identification of a Poissonian component at sub-arcminute scales in agreement with predictions of the cosmic infrared background level at SPIRE wavelengths. Finally, the comparison of the SPIRE and IRAS 100 micron data of the Polaris flare clearly assesses the capability of SPIRE in maping diffuse emission over large areas.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

We report the detection of strong absorption by interstellar hydrogen fluoride along the sight-line to the submillimeter continuum source G10.6-0.4 (W31C). We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1232.4763 GHz J=1-0 HF transition in the upper sideband of the Band 5a receiver. The resultant spectrum shows weak HF emission from G10.6-0.4 at LSR velocities in the range -10 to -3 km/s, accompanied by strong absorption by foreground material at LSR velocities in the range 15 to 50 km/s. The spectrum is similar to that of the 1113.3430 GHz 1(11)-0(00) transition of para-water, although at some frequencies the HF (hydrogen fluoride) optical depth clearly exceeds that of para-H2O. The optically-thick HF absorption that we have observed places a conservative lower limit of 1.6E+14 cm-2 on the HF column density along the sight-line to G10.6-0.4. Our lower limit on the HF abundance, 6E-9 relative to hydrogen nuclei, implies that hydrogen fluoride accounts for between ~ 30 and 100% of the fluorine nuclei in the gas phase along this sight-line. This observation corroborates theoretical predictions that - because the unique thermochemistry of fluorine permits the exothermic reaction of F atoms with molecular hydrogen - HF will be the dominant reservoir of interstellar fluorine under a wide range of conditions.

We present far-infrared spectra and maps of the DR21 molecular cloud core between 196 and 671 microns, using the Herschel-SPIRE spectrometer. Nineteen molecular lines originating from CO, 13CO, HCO+ and H2O, plus lines of [N II] and [CI] were recorded, including several transitions not previously detected. The CO lines are excited in warm gas with Tkin ~ 125 K and nH2 ~ 7 x 10^4 cm-3, CO column density N(CO) ~ 3.5 x 10^18 cm^-2 and a filling factor of ~ 12%, and appear to trace gas associated with an outflow. The rotational temperature analysis incorporating observations from ground-based telescopes reveals an additional lower excitation CO compoment which has a temperature ~ 78 K and N(CO) ~ 4.5 x 10^21 cm^-2. Astronomy & Astrophysics HERSCHEL special Issue, in press.

Context: Sh2-104 is a Galactic H ii region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observational confirmation of the collect-and-collapse model of triggered star-formation. Aims: We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In addition, we investigate the relationship between the dust emissivity index {\beta} and the dust temperature, T_dust. Methods: Using Herschel PACS and SPIRE images at 100, 160, 250, 350 and 500 {\mu}m we determine T_dust and {\beta} throughout Sh2-104, fitting the spectral energy distributions (SEDs) obtained from aperture photometry. With the SPIRE Fourier transform spectrometer (FTS) we obtained spectra at different positions in the Sh2-104 region. We detect J-ladders of CO and 13CO, with which we derive the gas temperature and column density. We also detect proxies of ionizing flux as the [NII] 3P1-3P0 and [CI] 3P2-3P1 transitions. Results: We find an average value of {\beta} ~ 1.5 throughout Sh2-104, as well as a T dust difference between the photodissociation region (PDR, ~ 25 K) and the interior (~ 40 K) of the bubble. We recover the anti-correlation between {\beta} and dust temperature reported numerous times in the literature. The relative isotopologue abundances of CO appear to be enhanced above the standard ISM values, but the obtained value is very preliminary and is still affected by large uncertainties.

The electronic spectra of gas-phase cationic polycyclic aromatic hydrocarbons (PAHs), trapped in the Fourier Transform Ion Cyclotron Resonance cell of the PIRENEA experiment, have been measured by multiphoton dissociation spectroscopy in the 430-480 nm spectral range using the radiation of a mid-band optical parametric oscillator laser. We present here the spectra recorded for different species of increasing size, namely the pyrene cation (C16H10+), the 1-methylpyrene cation (CH3-C16H9+), the coronene cation (C24H12+), and its dehydrogenated derivative C24H10+. The experimental results are interpreted with the help of time-dependent density functional theory calculations and analysed using spectral information on the same species obtained from matrix isolation spectroscopy data. A kinetic Monte Carlo code has also been used, in the case of pyrene and coronene cations, to estimate the absorption cross-sections of the measured electronic transitions. Gas-phase spectra of highly reactive species such as dehydrogenated PAH cations are reported for the first time.

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 molecular gas in the DR21 massive star formation region is known to be affected by the strong UV field from the central star cluster and by a fast outflow creating a bright shock. The relative contribution of both heating mechanisms is the matter of a long debate. By better sampling the excitation ladder of various tracers we provide a quantitative distinction between the different heating mechanisms. HIFI observations of mid-J transitions of CO and HCO+ isotopes allow us to bridge the gap in excitation energies between observations from the ground, characterizing the cooler gas, and existing ISO LWS spectra, constraining the properties of the hot gas. Comparing the detailed line profiles allows to identify the physical structure of the different components. In spite of the known shock-excitation of H2 and the clearly visible strong outflow, we find that the emission of all lines up to > 2 THz can be explained by purely radiative heating of the material. However, the new Herschel/HIFI observations reveal two types of excitation conditions. We find hot and dense clumps close to the central cluster, probably dynamically affected by the outflow, and a more widespread distribution of cooler, but nevertheless dense, molecular clumps.

We use the SPIRE Fourier-Transform Spectrometer (FTS) on-board the ESA Herschel Space Telescope to analyse the submillimetre spectrum of the Ultra-compact HII region G29.96-0.02. Spectral lines from species including 13CO, CO, [CI], and [NII] are detected. A sparse map of the [NII] emission shows at least one other HII region neighbouring the clump containing the UCHII. The FTS spectra are combined with ISO SWS and LWS spectra and fluxes from the literature to present a detailed spectrum of the source spanning three orders of magnitude in wavelength. The quality of the spectrum longwards of 100 {\mu}m allows us to fit a single temperature greybody with temperature 80.3\pm0.6K and dust emissivity index 1.73\pm0.02, an accuracy rarely obtained with previous instruments. We estimate a mass of 1500 Msol for the clump containing the HII region. The clump's bolometeric luminosity of 4 x 10^6 Lsol is comparable to, or slightly greater than, the known O-star powering the UCHII region.

Aims. To follow the species chemistry arising in diverse sources of the Galaxy with Herschel. Methods. SPIRE FTS sparse sampled maps of the Orion bar & compact HII regions G29.96-0.02 and G32.80+0.19 have been analyzed. Results. Beyond the wealth of atomic and molecular lines detected in the high-resolution spectra obtained with the FTS of SPIRE in the Orion Bar, one emission line is found to lie at the position of the fundamental rotational transition of CH+ as measured precisely in the laboratory (Pearson & Drouion 2006). This coincidence suggests that it is the first detection of the fundamental rotational transition of CH+. This claim is strengthened by the observation of the lambda doublet transitions arising from its relative, CH, which are also observed in the same spectrum. The broad spectral coverage of the SPIRE FTS allows for the simultaneous measurement of these closely related chemically species, under the same observing conditions. The importance of these lines are discussed and a comparison with results obtained from models of the Photon Dominated Region (PDR) of Orion are presented. The CH+ line also appears in absorption in the spectra of the two galactic compact HII regions G29.96-0.02 and G32.80+0.19, which is likely due to the presence of CH+ in the the Cold Neutral Medium of the galactic plane. These detections will shed light on the formation processes and on the existence of CH+, which are still outstanding questions in astrophysics.

Context. RCW 120 is a well-studied, nearby Galactic HII region with ongoing star formation in its surroundings. Previous work has shown that it displays a bubble morphology at mid-infrared wavelengths and has a massive layer of collected neutral material seen at sub-mm wavelengths. Given the well-defined photo-dissociation region (PDR) boundary and collected layer, it is an excellent laboratory to study the "collect and collapse" process of triggered star formation. Using Herschel Space Observatory data at 100, 160, 250, 350, and 500 micron, in combination with Spitzer and APEX-LABOCA data, we can for the first time map the entire spectral energy distribution of an HII region at high angular resolution. Aims. We seek a better understanding of RCW120 and its local environment by analysing its dust temperature distribution. Additionally, we wish to understand how the dust emissivity index, beta, is related to the dust temperature. Methods. We determine dust temperatures in selected regions of the RCW 120 field by fitting their spectral energy distribution (SED), derived using aperture photometry. Additionally, we fit the SED extracted from a grid of positions to create a temperature map. Results. We find a gradient in dust temperature, ranging from >30 K in the interior of RCW 120, to ~20K for the material collected in the PDR, to ~10K toward local infrared dark clouds and cold filaments. Our results suggest that RCW 120 is in the process of destroying the PDR delineating its bubble morphology. The leaked radiation from its interior may influence the creation of the next generation of stars. We find support for an anti-correlation between the fitted temperature and beta, in rough agreement with what has been found previously. The extended wavelength coverage of the Herschel data greatly increases the reliability of this result.

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.

Mid-infrared (IR) observations of polycyclic aromatic hydrocarbons (PAHs) and molecular hydrogen emission are a potentially powerful tool to derive physical properties of dense environments irradiated by intense UV fields. We present new, spatially resolved, \emph{Spitzer} mid-IR spectroscopy of the high UV-field and dense photodissocation region (PDR) around Monoceros R2, the closest ultracompact \hII region, revealing the spatial structure of ionized gas, PAHs and H$_2$ emissions. Using a PDR model and PAH emission feature fitting algorithm, we build a comprehensive picture of the physical conditions prevailing in the region. We show that the combination of the measurement of PAH ionization fraction and of the ratio between the H$_2$ 0-0 S(3) and S(2) line intensities, respectively at 9.7 and 12.3 $\mu$m, allows to derive the fundamental parameters driving the PDR: temperature, density and UV radiation field when they fall in the ranges $T = 250-1500 $K, $n_H=10^4-10^6$cm$^{-3}$, $G_0=10^3-10^5$ respectively. These mid-IR spectral tracers thus provide a tool to probe the similar but unresolved UV-illuminated surface of protoplanetary disks or the nuclei of starburst galaxies.

Polycyclic Aromatic Hydrocarbons (PAHs) are widely accepted as the carriers of the Aromatic Infrared Bands (AIBs), but an unambiguous identification of any specific interstellar PAH is still missing. For polar PAHs, pure rotational transitions can be used as fingerprints for identification. Combining dedicated experiments, detailed simulations and observations, we explored the mm domain to search for specific rotational transitions of corannulene (C20H10). We performed high-resolution spectroscopic measurements and a simulation of the emission spectrum of UV-excited C20H10 in the environment of the Red Rectangle, calculating its synthetic rotational spectrum. Based on these results, we conducted a first observational campaign at the IRAM 30m telescope towards this source to search for several high-J rotational transitions of (C20H10). The laboratory detection of the J = 112 <- 111 transition of corannulene showed that no centrifugal splitting is present up to this line. Observations with the IRAM 30m telescope towards the Red Rectangle do not show any corannulene emission at any of the observed frequencies, down to a rms noise level of Tmb = 8 mK for the J =135 -> 134 transition at 137.615 GHz. Comparing the noise level with the synthetic spectrum, we are able to estimate an upper limit to the fraction of carbon locked in corannulene of about 1.0x10(-5) relative to the total abundance of carbon in PAHs. The sensitivity achieved shows that radio spectroscopy can be a powerful tool to search for polar PAHs. We compare this upper limit with models for the PAH size distribution, emphasising that small PAHs are much less abundant than predicted. We show that this cannot be explained by destruction but is more likely related to the chemistry of their formation in the environment of the Red Rectangle.

Polycyclic Aromatic Hydrocarbons (PAHs) are considered as a major constituent of interstellar dust. They have been proposed as the carriers of the Aromatic Infrared Bands (AIBs) observed in emission in the mid-IR. They likely have a significant contribution to various features of the extinction curve such as the 220 nm bump,the far-UV rise and the diffuse interstellar bands. Emission bands are also expected in the far-IR, which are better fingerprints of molecular identity than the AIBs. They will be searched for with the Herschel Space Observatory. Rotational emission is also expected in the mm range for those molecules which carry significant dipole moments. Despite spectroscopic studies in the laboratory, no individual PAH species could be identified. This emphasises the need for an investigation on where interstellar PAHs come from and how they evolve due to environmental conditions: ionisation and dissociation upon UV irradiation, interactions with electrons, gas and dust. There is also evidence for PAH species to contribute to the depletion of heavy atoms from the gas phase, in particular Si and Fe. This paper illustrates how laboratory work can be inspired from observations. In particular there is a need for understanding the chemical properties of PAHs and PAH-related species, including very small grains, in physical conditions that mimic those found in interstellar space. This motivates a joint effort between astrophysicists, physicists and chemists. Such interdisciplinary studies are currently performed, taking advantage of the PIRENEA set-up, a cold ion trap dedicated to astrochemistry.

We completed a search for circumstellar disks around Herbig Be stars using the NRAO Very Large Array (VLA) and the IRAM Plateau de Bure (PdB) interferometers. We present our new VLA and PdBI data for the three objects MWC 297, Z CMa, and LKHa 215. We constructed the SED from near-IR to centimeter wavelengths by adding our millimeter and centimeter data to the available data at other wavelengths, mainly Spitzer images. The entire SED was fitted using a disk+envelope model. In addition, we compiled all the disk millimeter observations in the literature and completed a statistical analysis of all the data. We show that the disk mass is usually only a small percentage (less than 10%) of the mass of the entire envelope in HBe stars. For the disks, there are large source-to-source variations. Two disks in our sample, R Mon and Z CMa, have similar sizes and masses to those found in T Tauri and Herbig Ae stars. The disks around MWC 1080 and MWC 297 are, however, smaller (rout<100 AU). We did not detect the disks towards MWC 137 and LkHa 215 at millimeter wavelengths, which limits the mass and the size of the possible circumstellar disks. A comparison between our data and previous results for T Tauri and Herbig Ae stars indicates that although massive disks (0.1 Msun) are found in young objects (10^4 yr), the masses of the disks around Herbig Be stars are usually 5-10 times lower than those around lower mass stars. We propose that disk photoevaporation is responsible for this behavior. In Herbig Be stars, the UV radiation disperses the gas in the outer disk on a timescale of a few 10^5 yr. Once the outer part of the disk has vanished, the entire gaseous disk is photoevaporated on a very short timescale (10^5 yr) and only a small, dusty disk consisting of large grains remains.

In this Paper we analyze the mid-infrared (mid-IR) emission of very small dust particles in a sample of 12 protoplanetary disks to see how they are connected to interstellar dust particles and to investigate the possibility that their emission can be used as a probe of the physical conditions and evolution of the disk. We define a basis made of three mid-IR template spectra PAH$^0$, PAH$^+$ and VSGs that were derived from the analysis of reflection nebulae, and an additional PAH$^x$ spectrum that was introduced by Joblin et al. (2008) for the analysis of the spectra of planetary nebulae. From the optimization of the fit of 12 star+disk spectra, using a linear combination of the 4 template spectra, we found that an additional small grain component with a broad feature at 8.3 $\mu$m is needed. We find that the fraction of VSG emission in disks decreases with increasing stellar temperature. VSGs appear to be destroyed by UV photons at the surface of disks, thus releasing free PAH molecules, which are eventually ionized as it is observed in photodissociation regions. On the opposite, we observe that the fraction of PAH$^x$ increases with increasing star temperature except in the case of B stars where they are absent. We argue that this is compatible with the identification of PAH$^x$ as large ionized PAHs, most likely emitting in regions of the disk that are close to the star. Finally, we provide a UV-dependant scheme to explain the evolution of PAHs and VSGs in protoplanetary disks. We show that A stars modify the size spectrum of PAHs and VSGs in favor of large PAHs while B stars destroy even the largest PAHs up to large radii in the disk. These results allow us to put new constrains on the properties of two sources: IRS 48 and "Gomez's Hamburger" which are poorly characterized.

We investigate the presence of silicon atoms adsorbed on the surface of interstellar polycyclic aromatic hydrocarbons (PAHs) to form SiPAH pi-complexes. We use quantum chemistry calculations to obtain structural, thermodynamic and mid-IR properties of neutral and cationic SiPAH complexes. The binding energy was found to be at least 1.5 eV for [SiPAH]+ complexes whereas it is roughly 0.5 eV for their neutral counterparts. From the spectral analysis of the calculated IR spectra, we found that the coordination of silicon to PAH+ does not strongly affect the intensities of the PAH+ spectra, but systematically introduces blueshifts of the C-C in-plane and the C-H out-of-plane bands. The thermodynamic data calculated for [SiPAH]+ complexes show that these species are stable and can be easily formed by radiative association of Si+ and PAH species that are known to be abundant in photodissociation regions. Their mid-IR fingerprints show features induced by the coordination of silicon that could account for (i) the blueshifted position of the 6.2 micron AIB and (ii) the presence of satellite bands observed on the blue side of the 6.2 and 11.2 micron AIBs. From such an assignment, we can deduce that typically 1% of the cosmic silicon appears to be attached to PAHs.

It has been shown that the diversity of the aromatic emission features can be rationalized into different classes of objects, in which differences between circumstellar and interstellar matter are emphasised. We probe the links between the mid-IR emitters observed in planetary nebulae (PNe) and their counterparts in the interstellar medium in order to probe a scenario in which the latter have been formed in the circumstellar environment of evolved stars. The mid-IR (6-14 um) emission spectra of PNe and compact HII regions were analysed on the basis of previous work on photodissociation regions (PDRs). Galactic, Large Magellanic Cloud (LMC), and Small Magellanic Cloud (SMC) objects were considered in our sample.We show that the mid-IR emission of PNe can be decomposed as the sum of six components. Some components made of polycyclic aromatic hydrocarbon (PAH) and very small grain (VSG) populations are similar to those observed in PDRs. Others are fitted in an evolutionary scenario involving the destruction of the aliphatic component observed in the post-AGB stage, as well as strong processing of PAHs in the extreme conditions of PNe that leads to a population of very large ionized PAHs. This species called PAH^x are proposed as the carriers of a characteristic band at 7.90 um. This band can be used as part of diagnostics that identify PNe in nearby galaxies and is also observed in galactic compact HII regions. These results support the formation of the aromatic very small dust particles in the envelopes of evolved stars, in the Milky Way, as well as in the LMC and SMC, and their subsequent survival in the interstellar medium.

Recent calculations have shown that the UV bump at about 217.5 nm in the extinction curve can be explained by a complex mixture of PAHs in several charge states. Other studies proposed that the carriers are a restricted population made of neutral and singly-ionised dehydrogenated coronene molecules (C24Hn, n less than 3), in line with models of the hydrogenation state of interstellar PAHs predicting that medium-sized species are highly dehydrogenated. To assess the observational consequences of the latter hypothesis we have undertaken a systematic study of the electronic spectra of dehydrogenated PAHs. We use our first results to see whether such spectra show strong general trends upon dehydrogenation. We used state-of-the-art techniques in the framework of the density functional theory (DFT) to obtain the electronic ground-state geometries, and of the time- dependent DFT to evaluate the electronic excited-state properties. We computed the absorption cross-section of the species C24Hn (n=12,10,8,6,4,2,0) in their neutral and cationic charge-states. Similar calculations were performed for other PAHs and their fullydehydrogenated counterparts. pi electron energies are always found to be strongly affected by dehydrogenation. In all cases we examined, progressive dehydrogenation translates into a correspondingly progressive blue shift of the main electronic transitions. In particular, the pi-pi* collective resonance becomes broader and bluer with dehydrogenation. Its calculated energy position is therefore predicted to fall in the gap between the UV bump and the far-UV rise of the extinction curve. Since this effect appears to be systematic, it poses a tight observational limit on the column density of strongly dehydrogenated medium-sized PAHs.

Aims: We explore the relation between charge state of polycyclic aromatic hydrocarbons (PAHs) and extinction curve morphology. Methods: We fit extinction curves with a dust model including core-mantle spherical particles of mixed chemical composition (silicate core, $sp^2$ and $sp^3$ carbonaceous layers), and an additional molecular component. We use exact methods to calculate the extinction due to classical particles and accurate computed absorption spectra of PAHs in different charge states, for the contribution due to the molecular component, along a sample of five rather different lines of sight. Results: A combination of classical dust particles and mixtures of real PAHs satisfactorily matches the observed interstellar extinction curves. Variations of the spectral properties of PAHs in different charge states produce changes consistent with the varying relative strengths of the bump and non-linear far-UV rise.

We present the detection and characterization of a peculiar low-mass protostar (IRAS 22129+7000) located ~0.4 pc from Ced 201 Photodissociation Region (PDR) and ~0.2 pc from the HH450 jet. The cold circumstellar envelope surrounding the object has been mapped through its 1.2 mm dust continuum emission with IRAM-30m/MAMBO. The deeply embedded protostar is clearly detected with Spitzer/MIPS (70 um), IRS (20-35 um) and IRAC (4.5, 5.8, and 8 um) but also in the K_s band (2.15 um). Given the large "near- and mid-IR excess" in its spectral energy distribution, but large submillimeter-to-bolometric luminosity ratio (~2%), IRAS 22129+7000 must be a transition Class 0/I source and/or a multiple stellar system. Targeted observations of several molecular lines from CO, 13CO, C18O, HCO+ and DCO+ have been obtained. The presence of a collimated molecular outflow mapped with the CSO telescope in the CO J=3-2 line suggests that the protostar/disk system is still accreting material from its natal envelope. Indeed, optically thick line profiles from high density tracers such as HCO+ J=1-0 show a red-shifted-absorption asymmetry reminiscent of inward motions. We construct a preliminary physical model of the circumstellar envelope (including radial density and temperature gradients, velocity field and turbulence) that reproduces the observed line profiles and estimates the ionization fraction. The presence of both mechanical and (non-ionizing) FUV-radiative input makes the region an interesting case to study triggered star formation.

Extended Red Emission (ERE) was recently attributed to the photo-luminescence of either doubly ionized Polycyclic Aromatic Hydrocarbons (PAH$^{++}$), or charged PAH dimers. We analysed the visible and mid-infrared (mid-IR) dust emission in the North-West and South photo-dissociation regions of the reflection nebula NGC 7023.Using a blind signal separation method, we extracted the map of ERE from images obtained with the Hubble Space Telescope, and at the Canada France Hawaii Telescope. We compared the extracted ERE image to the distribution maps of the mid-IR emission of Very Small Grains (VSGs), neutral and ionized PAHs (PAH$^0$ and PAH$^+$) obtained with the Spitzer Space Telescope and the Infrared Space Observatory. ERE is dominant in transition regions where VSGs are being photo-evaporated to form free PAH molecules, and is not observed in regions dominated by PAH$^+$. Its carrier makes a minor contribution to the mid-IR emission spectrum. These results suggest that the ERE carrier is a transition species formed during the destruction of VSGs. Singly ionized PAH dimers appear as good candidates but PAH$^{++}$ molecules seem to be excluded.

The use of Blind Signal Separation methods (ICA and other approaches) for the analysis of astrophysical data remains quite unexplored. In this paper, we present a new approach for analyzing the infrared emission spectra of interstellar dust, obtained with NASA's Spitzer Space Telescope, using FastICA and Non-negative Matrix Factorization (NMF). Using these two methods, we were able to unveil the source spectra of three different types of carbonaceous nanoparticles present in interstellar space. These spectra can then constitute a basis for the interpretation of the mid-infrared emission spectra of interstellar dust in the Milky Way and nearby galaxies. We also show how to use these extracted spectra to derive the spatial distribution of these nanoparticles.

We study the evolution of the Aromatic Infrared Bands (AIBs) emitters across the illuminated edge of the Horsehead nebula and especially their survival and properties in the HII region. We present spectral mapping observations taken with the Infrared Spectrograph (IRS) at wavelengths 5.2-38 microns. A strong AIB at 11.3 microns is detected in the HII region, relative to the other AIBs at 6.2, 7.7 and 8.6 microns. The intensity of this band appears to be correlated with the intensity of the [NeII] at 12.8 microns and of Halpha, which shows that the emitters of the 11.3 microns band are located in the ionised gas. The survival of PAHs in the HII region could be due to the moderate intensity of the radiation field (G0 about 100) and the lack of photons with energy above about 25eV. The enhancement of the intensity of the 11.3 microns band in the HII region, relative to the other AIBs can be explained by the presence of neutral PAHs. Our observations highlight a transition region between ionised and neutral PAHs observed with ideal conditions in our Galaxy. A scenario where PAHs can survive in HII regions and be significantly neutral could explain the detection of a prominent 11.3 microns band in other Spitzer observations.

This work was conducted as part of the SPECPDR program, dedicated to the study of very small particles and astrochemistry, in Photo-Dissociation Regions (PDRs). We present the analysis of the mid-IR spectro-imagery observations of Ced 201, NCG 7023 East and North-West and rho-Ophiuchi West filament. Using the data from all four modules of the InfraRed Spectrograph onboard the Spitzer Space Telescope, we produced a spectral cube ranging from 5 to 35 um for each one of the observed PDRs. The resulting cubes were analysed using Blind Signal Separation methods (NMF and FastICA). For Ced 201, rho-Ophiuchi West filament and NGC 7023 East, we find that two signals can be extracted from the original data cubes, which are 5 to 35 um spectra. The main features of the first spectrum are a strong continuum emission at long wavelengths, and a broad 7.8 um band. On the contrary, the second spectrum exhibits the classical Aromatic Infrared Bands (AIBs) and no continuum. The reconstructed spatial distribution maps show that the latter spectrum is mainly present at the cloud surface, close to the star whereas the first one is located slightly deeper inside the PDR. The study of the spectral energy distribution of Ced 201 up to 100 um suggests that, in cool PDRs, the 5-25 um continuum is carried by Very Small Grains (VSGs). The AIB spectra in the observed objects can be interpreted as the contribution of neutral and positively-charged Polycyclic Aromatic Hydrocarbons (PAHs). We extracted the 5 to 25 um emission spectrum of VSGs in cool PDRs, these grains being most likely carbonaceous. We show that the variations of the mid-IR (5-35 um) spectra of PDRs can be explained by the photo-chemical processing of these VSGs and PAHs, VSGs being the progenitors of free PAH.

We present an on-line database of computed molecular properties for a large sample of polycyclic aromatic hydrocarbons (PAHs) in four charge states: -1, 0, +1, and +2. At present our database includes 40 molecules ranging in size from naphthalene and azulene (C10H8) up to circumovalene (C66H20). We performed our calculations in the framework of the density functional theory (DFT) and the time-dependent DFT to obtain the most relevant molecular parameters needed for astrophysical applications. For each molecule in the sample, our database presents in a uniform way the energetic, rotational, vibrational, and electronic properties. It is freely accessible on the web at http://astrochemistry.ca.astro.it/database/ and http://www.cesr.fr/~joblin/database/.

Aims. We present a systematic theoretical study on 40 polycyclic aromatic hydrocarbons dications (PAHs++) containing up to 66 carbon atoms. Methods. We performed our calculations using well established quantum-chemical techniques in the framework of the density functional theory (DFT) to obtain the electronic ground-state properties, and of the time-dependent DFT (TD-DFT) to evaluate the excited\textendash state properties. Results. For each PAH++ considered we computed the absolute visible-UV photo-absorption cross-section up to about 30 eV. We also evaluated their vibrational properties and compared them to those of the corresponding neutral and singly-ionised species. We estimated the adiabatic and vertical second ionisation energy DeltaI through total energy differences. Conclusions. The DeltaI values obtained fall in the energy range 8-13 eV, confirming that PAHs could reach the doubly-ionised state in HI regions. The total integrated IR absorption cross-sections show a marked increase upon ionization, being on the average about two and five times larger for PAHs++ than for PAHs+ and PAHs, respectively. The visible-UV photo-absorption cross-sections for the 0, +1 and +2 charge-states show comparable features but PAHs++ are found to absorb slightly less than their parent neutrals and singly ionized species between ~7 and ~12 eV. Combining these pieces of information we found that PAHs++ should actually be more stable against photodissociation than PAHs and PAHs+,_if_ dissociation tresholds are nearly unchanged by ionization.

The competition between the formation and destruction of coronene clusters under interstellar conditions is investigated theoretically. The unimolecular nucleation of neutral clusters is simulated with an atomic model combining an explicit classical force field and a quantum tight-binding approach. Evaporation rates are calculated in the framework of the phase space theory and are inserted in an infrared emission model and compared with the growth rate constants. It is found that, in interstellar conditions, most collisions lead to cluster growth. The time evolution of small clusters (containing up to 312 carbon atoms) was specifically investigated under the physical conditions of the northern photodissociation region of NGC 7023. These clusters are found to be thermally photoevaporated much faster than they are reformed, thus providing an interpretation for the lowest limit of the interstellar cluster size distribution inferred from observations. The effects of ionizing the clusters and density heterogeneities are also considered. Based on our results, the possibility that PAH clusters could be formed in PDRs is critically discussed.

Context. In the framework of the interstellar PAH hypothesis, far-IR skeletal bands are expected to be a fingerprint of single species in this class. Aims. A detailed model of the photophysics of interstellar PAHs is required for such single-molecule identification of their far-IR features in the presently available Infrared Space Observatory data and in those of the forthcoming Herschel Space Observatory mission. Methods. We modelled the detailed photophysics of a vast sample of species in different radiation fields, using a compendium of Monte-Carlo techniques and quantum-chemical calculations. This enabled us to validate the use of purely theoretical data and assess the expected accuracy and reliability of the resulting synthetic far-IR emission spectra. Results. We produce positions and intensities of the expected far-IR features which ought to be emitted by each species in the sample in the considered radiation fields. A composite emission spectrum for our sample is computed for one of the most favourable sources for detection, namely the Red Rectangle nebula. The resulting spectrum is compared with the estimated dust emission in the same source, to assess the dependence of detectability on key molecular parameters. Conclusions. Identifying specific PAHs from their far-IR features is going to be a difficult feat in general, still it may well be possible under favourable conditions.

Context. In the framework of the interstellar polycyclic aromatic hydrocarbons (PAHs) hypothesis, far-IR skeletal bands are expected to be a fingerprint of single species in this class. Aims. We address the question of detectability of low energy PAH vibrational bands, with respect to spectral contrast and intensity ratio with ``classical'' Aromatic Infrared Bands (AIBs). Methods. We extend our extablished Monte-Carlo model of the photophysics of specific PAHs in astronomical environments, to include rotational and anharmonic band structure. The required molecular parameters were calculated in the framework of the Density Functional Theory. Results. We calculate the detailed spectral profiles of three low-energy vibrational bands of neutral naphthalene, and four low-energy vibrational bands of neutral anthracene. They are used to establish detectability constraints based on intensity ratios with ``classical'' AIBs. A general procedure is suggested to select promising diagnostics, and tested on available Infrared Space Observatory data for the Red Rectangle nebula. Conclusions. The search for single, specific PAHs in the far-IR is a challenging, but promising task, especially in view of the forthcoming launch of the Herschel Space Observatory.

Following the tentative identification of the blue luminescence in the Red Rectangle by Vijh et al. (2005), we compute absolute fluxes for the vibrational IR emission and phosphorescence bands of three small polycyclic aromatic hydrocarbons. The calculated IR spectra are compared with available ISO observations. A subset of the emission bands are predicted to be observable using presently available facilities, and can be used for an immediate, independent, discriminating test on their alleged presence in this well-known astronomical object.