Aim. We investigate the role of PAHs as a sink for deuterium in the interstellar medium and study UV photolysis as a potential process in the variations of the deuterium fractionation in the ISM. Methods. The UV photo-induced fragmentation of various isotopologues of D-enriched, protonated anthracene and phenanthrene ions was recorded in a FTICR mass spectrometer. IRMPD spectroscopy using FELIX provided the IR spectra that were compared to DFT vibrational spectra; reaction barriers and rates were also calculated and related to the product abundances. Results. The mass spectra for both UV and IRMPD photolysis show the loss of H from [D-C$_{14}$H$_{10}$]$^+$, whereas [H-C$_{14}$D$_{10}$]$^+$ shows a strong preference for D loss. Calculations reveal facile 1,2-H and -D shift reactions, with barriers lower than the energy supplied by the photo-excitation process. Together with confirmation of the ground-state structures via the IR spectra, we determined that the photolytic processes in the 2 PAHs are largely governed by scrambling where the H and the D atoms relocate between different peripheral C atoms. The $\sim$0.1 eV difference in zero-point energy between C-H and C-D bonds ultimately leads to faster H scrambling than D scrambling, and increased H atom loss compared to D. Conclusion. Scrambling is common in PAH cations under UV radiation. Upon photoexcitation of deuterium-enriched PAHs, the scrambling results in a higher probability for the aliphatic D atom to migrate to an aromatic site, protecting it from elimination. This could lead to increased deuteration as a PAH moves towards more exposed interstellar environments. Also, large, compact PAHs with an aliphatic C-HD group on solo sites might be responsible for the majority of aliphatic C-D stretching bands seen in astronomical spectra.

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

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

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

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

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

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

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

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

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

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

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

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

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

The 11.2 $\mu$m emission band belongs to the family of the `Unidentified' Infrared (UIR) emission bands seen in many astronomical environments. In this work we present a theoretical interpretation of the band characteristics and profile variation for a number of astrophysical sources in which the carriers are subject to a range of physical conditions. The results of Density Functional Theory (DFT) calculations for the solo out-of-plane (OOP) vibrational bending modes of large polycyclic aromatic hydrocarbon (PAH) molecules are used as input for a detailed emission model which includes the temperature and mass dependence of PAH band wavelength, and a PAH mass distribution that varies with object. Comparison of the model with astronomical spectra indicates that the 11.2 $\mu$m band asymmetry and profile variation can be explained principally in terms of the mass distribution of neutral PAHs with a small contribution from anharmonic effects.

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

This work aims to disentangle the morphological, kinematic, and chemical components of the nebula M1-67 to shed light on its process of formation around the central Wolf-Rayet (WR) star WR124. We have carried out integral field spectroscopy observations over two regions of M1-67, covering most of the nebula in the optical range. Maps of electron density, line ratios, and radial velocity were created to perform a detailed analysis of the two-dimensional structure. We studied the physical and chemical properties by means of integrated spectra selected over the whole nebula. Photoionization models were performed to confirm the empirical chemical results theoretically. In addition, we analysed infrared spectroscopic data and the MIPS 24micron image of M1-67 from Spitzer. We find that the ionized gas of M1-67 is condensed in knots aligned in a preferred axis along the NE-SW direction, like a bipolar structure. Both electron density and radial velocity decrease in this direction when moving away from the central star. From the derived electron temperature, Te~8200 K, we have estimated chemical abundances, obtaining that nitrogen appears strongly enriched and oxygen depleted. From the last two results, we infer that this bipolarity is the consequence of an ejection of an evolved stage of WR124 with material processed in the CNO cycle. The infrared study has revealed that the bipolar axis is composed of ionized gas with a low ionization degree that is well mixed with warm dust and of a spherical bubble surrounding the ejection at 24micron. Taking the evolution of a 60 Mo star and the temporal scale of the bipolar ejection into account, we propose that the observed gas was ejected during an eruption in the luminous blue variable. The star has entered the WR phase recently without apparent signs of interaction between WR-winds and interstellar material.

The spatial distribution of 3.3 \mum PAH and associated emission in the 3.3" x 6.0" inner region of the Red Rectangle nebula has been determined using the UIST imager-spectrometer at the United Kingdom Infrared Telescope (UKIRT). Interpretation of the 3.3 \mum feature as comprising two spectroscopic components centred at 3.30 \mum and 3.28 \mum, as put forward by Song et al. (2003, 2007), is supported by these data which reveal that they have different spatial distributions. It is deduced that there are two classes of 3.3 \mum band carrier with a peak wavelength separation of ~0.02 \mum. From comparison of the 3.3 \mum observations with laboratory and theoretical spectra for a range of PAH molecules it is proposed that the 3.28 \mum and 3.30 \mum components arise from 'bay' and 'non-bay' hydrogen sites, respectively, on the periphery of small neutral PAHs. Observational data are also obtained for L-band continuum emission and for the Pfund \epsilon hydrogen recombination line.