No Arabic abstract
We have obtained new STIS/HST spectra to search for structure in the ultraviolet interstellar extinction curve, with particular emphasis on a search for absorption features produced by polycyclic aromatic hydrocarbons (PAHs). The presence of these molecules in the interstellar medium has been postulated to explain the infrared emission features seen in the 3-13 $mu$m spectra of numerous sources. UV spectra are uniquely capable of identifying specific PAH molecules. We obtained high S/N UV spectra of stars which are significantly more reddened than those observed in previous studies. These data put limits on the role of small (30-50 carbon atoms) PAHs in UV extinction and call for further observations to probe the role of larger PAHs. PAHs are of importance because of their ubiquity and high abundance inferred from the infrared data and also because they may link the molecular and dust phases of the interstellar medium. A presence or absence of ultraviolet absorption bands due to PAHs could be a definitive test of this hypothesis. We should be able to detect a 20 AA wide feature down to a 3$sigma$ limit of $sim$0.02 A$_V$. No such absorption features are seen other than the well-known 2175 AA bump.
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.
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)
We present a new method to accurately describe the ionization fraction and the size distribution of polycyclic aromatic hydrocarbons (PAHs) within astrophysical sources. To this purpose, we have computed the mid-infrared emission spectra of 308 PAH molecules of varying sizes, symmetries, and compactness, generated in a range of radiation fields. We show that the intensity ratio of the solo CH out-of-plane bending mode in PAH cations and anions (referred to as the 11.0 $mu$m band, falling in the 11.0-11.3 $mu$m region for cations and anions) to their 3.3 $mu$m emission, scales with PAH size, similarly to the scaling of the 11.2/3.3 ratio with the number of carbon atoms (N$_{mathrm{C}}$) for neutral molecules. Among the different PAH emission bands, it is the 3.3 $mu$m band intensity which has the strongest correlation with N$_{mathrm{C}}$, and drives the reported PAH intensity ratio correlations with N$_{mathrm{C}}$ for both neutral and ionized PAHs. The 6.2/7.7 intensity ratio, previously adopted to track PAH size, shows no evident scaling with N$_{mathrm{C}}$ in our large sample. We define a new diagnostic grid space to probe PAH charge and size, using the (11.2+11.0)/7.7 and (11.2+11.0)/3.3 PAH intensity ratios respectively. We demonstrate the application of the (11.2+11.0)/7.7 - (11.2+11.0)/3.3 diagnostic grid for galaxies M82 and NGC 253, for the planetary nebula NGC 7027, and the reflection nebulae NGC 2023 and NGC 7023. Finally, we provide quantitative relations for PAH size determination depending on the ionization fraction of the PAHs and the radiation field they are exposed to.
The mid-infrared spectra of large PAHs ranging from C54H18 to C130H28 are determined computationally using Density Functional Theory. Trends in the band positions and intensities as a function of PAH size, charge and geometry are discussed. Regarding the 3.3, 6.3 and 11.2 micron bands similar conclusions hold as with small PAHs. This does not hold for the other features. The larger PAH cations and anions produce bands at 7.8 micron and, as PAH sizes increases, a band near 8.5 micron becomes prominent and shifts slightly to the red. In addition, the average anion peak falls slightly to the red of the average cation peak. The similarity in behavior of the 7.8 and 8.6 micron bands with the astronomical observations suggests that they arise from large, cationic and anionic PAHs, with the specific peak position and profile reflecting the PAH cation to anion concentration ratio and relative intensities of PAH size. Hence, the broad astronomical 7.7 micron band is produced by a mixture of small and large PAH cations and anions, with small and large PAHs contributing more to the 7.6 and 7.8 micron component respectively. For the CH out-of-plane vibrations, the duo hydrogens couple with the solo vibrations and produce bands that fall at wavelengths slightly different than their counterparts in smaller PAHs. As a consequence, previously deduced PAH structures are altered in favor of more compact and symmetric forms. In addition, the overlap between the duo and trio bands may reproduce the blue-shaded 12.8 micron profile.
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.