HIGH-RESOLUTION IR ABSORPTION SPECTROSCOPY OF POLYCYCLIC AROMATIC HYDROCARBONS: THE REALM OF ANHARMONICITY

Vibrational spectroscopy of interstellar PAHs: 1-cyanonaphthalene and 2-cyanonaphthalene

In this work, we present a comprehensive experimental and theoretical study of the vibrational spectra of PAH molecules recently detected in the interstellar medium: 1-cyanonaphthalene and 2-cyanonaphthalene. The room temperature IR spectra of 1- and 2-cyanonaphthalene in the region 100-3100 cm-1 and their vibrational Raman spectra in the region 35-3100 cm-1 are reported here for the first time. A detailed spectral analysis is carried out using quantum chemical calculations employing the DFT methodology. Anharmonic corrections using the VPT2 method yield excellent agreement with the experimental spectra. A re-investigation of the vibrational spectrum of the parent molecule: naphthalene validates the experimental and theoretical methods used. A consistent set of assignments is reported for the fundamental bands of 1- and 2-cyanonapththalene. The experimental and theoretical data presented here would be useful inputs for modelling the role of cyanonaphthalene in astrophysical processes.

Performance of Effective Harmonic Oscillator Approach for the Calculations of Vibrational Transition Energies of Large Molecules

The accuracy and performance of the effective harmonic oscillator approximation for the description of anharmonic vibrational structure calculations are tested for large molecular systems and compared with experimental values along with vibrational self-consistent field and second-order perturbation theories. The effective harmonic oscillator approach is an effective single-particle approximation where the variational parameters are the centroids and widths of the multidimensional Gaussian product functions posited as the vibrational wave functions. A comprehensive calculation for 849 transitions that include the fundamentals, two and three quanta overtone transitions, and several combination bands of three polyaromatic hydrocarbons and one DNA nucleobase with a total of 231 normal modes are assessed. A comparison of EHO results with the experimental values is done for the polyaromatic hydrocarbons, and a close agreement is found between the two results. It also offers anharmonic eigenstates and eigenfunctions that are nearly identical with vibrational self-consistent field theory. An extensive analysis on the resultant wave functions of the excited states is performed. The overall root-mean-square deviation (RMSD) between these two methods for 849 transitions understudy is only about 8.3 cm-1, suggesting the effective harmonic oscillator as a viable alternative for the reliable calculations of transition energies of large molecular systems.

Picking up Good Vibrations through Quartic Force Fields and Vibrational Perturbation Theory

Quartic force fields (QFFs) define sparse potential energy surfaces (compared to semiglobal surfaces) that are the cheapest and easiest means of computing anharmonic vibrational frequencies, especially when utilized with second-order vibrational perturbation theory (VPT2). However, flat and shallow potential surfaces are exceedingly difficult for QFFs to treat through a combination of numerical noise in the often numerically computed derivatives and in competing energy factors in the composite energies often utilized to provide high-level spectroscopic predictions. While some of these issues can be alleviated with analytic derivatives, hybrid QFFs, and intelligent choices in coordinate systems, the best practice is for predicting good molecular vibrations via QFFs is to understand what they cannot do, and this manuscript documents such cases where QFFs may fail.

Assigning the CH stretch overtone spectrum of benzene and naphthalene with extension to anthracene and tetracene using 2- and 3-quanta anharmonic quantum chemical computations

The CH stretch overtone region (5750-6300 cm-1) of benzene and naphthalene is assigned herein using anharmonic quantum chemical computations, and the trend of how this extends to larger polycyclic aromatic hydrocarbons (PAHs) is established. The assignment of all experimental bands to specific vibrational states is performed for the first time. Resonance polyads and the inclusion of 3-quanta vibrational states are both needed to compute accurate vibrational frequencies with the proper density-of-states to match the experimental band shape. Hundreds of 3-quanta states produce the observed band structure in naphthalene, anthracene, and tetracene, and this number is expected to increase drastically for larger PAHs. The width and shape of the main peak are consistent from naphthalene to anthracene, necessitating further exploration of this trend to confirm whether it is representative of all PAHs in the CH stretch overtone region. Understanding observations of PAH sources in the 1-3 μm region from the NIRSpec instrument aboard JWST requires new computational data, and this study provides a benchmark and foundation for their computation.

The infrared absorption spectrum of phenylacetylene and its deuterated isotopologue in the mid- to far-IR

Anharmonicity strongly influences the absorption and emission spectra of polycyclic aromatic hydrocarbon (PAH) molecules. Here, IR-UV ion-dip spectroscopy experiments together with detailed anharmonic computations reveal the presence of fundamental, overtone, as well as 2- and 3-quanta combination band transitions in the far- and mid-infrared absorption spectra of phenylacetylene and its singly deuterated isotopologue. Strong absorption features in the 400-900 cm-1 range originate from CH(D) in-plane and out-of-plane wags and bends, as well as bending motions including the C≡C and CH bonds of the acetylene substituent and the aromatic ring. For phenylacetylene, every absorption feature is assigned either directly or indirectly to a single or multiple vibrational mode(s). The measured spectrum is dense, broad, and structureless in many regions but well characterized by computations. Upon deuteration, large isotopic shifts are observed. At frequencies above 1500 cm-1 for d1-phenylacetylene, a one-to-one match is seen when comparing computations and experiments with all features assigned to combination bands and overtones. The C≡C stretch observed in phenylacetylene is not observed in d1-phenylacetylene due to a computed 40-fold drop in intensity. Overall, a careful treatment of anharmonicity that includes 2- and 3-quanta modes is found to be crucial to understand the rich details of the infrared spectrum of phenylacetylene. Based on these results, it can be expected that such an all-inclusive anharmonic treatment will also be key for unraveling the infrared spectra of PAHs in general.

Quantum Chemistry and Astrochemistry: A Match Made in the Heavens

Quantum chemistry can uniquely answer astrochemical questions that no other technique can provide. Computations can be parallelized, automated, and left to run continuously providing exceptional molecular throughput that cannot be done through experimentation. Additionally, the granularity of the individual computations that are required of potential energy surfaces, reaction mechanism pathways, or other quantum chemically derived observables produces a unique mosaic that make up the larger whole. These pieces can be dissected for their individual contributions or evaluated in an ad hoc fashion for each of their roles in generating the larger whole. No other scientific approach is capable of reporting such fine-grained insights. Quantum chemistry also works from a bottom-up approach in providing properties directly from the desired molecule instead of a top-down perspective as required of experiment where molecules have to be linked to observed phenomena. Furthermore, modern quantum chemistry is well within the range of "chemical accuracy" and is approaching "spectroscopic accuracy." As such, the seemingly difficult questions asked by astrochemistry that would not be asked initially for any other application require quantum chemical reference data. While the results of quantum chemical computations are needed to interpret astrochemical observation, modeling, or laboratory experimentation, such hard questions, regardless of the original need to answer them, produce unique solutions. While questions in astrochemistry often require novel developments in and implementations of quantum chemistry as outlined herein, the applications of these solutions will stretch beyond astrochemistry and may yet impact fields much closer to Earth.

An approximation to the vibrational coupled-cluster method for CH-stretching of large molecules: application to naphthalene and anthracene

We propose an approximation to the vibrational coupled-cluster method (VCCM) to describe the CH-stretching region of the vibrational spectrum of large molecules. The vibrational modes of a molecule are divided into two sets: the target set and the bath set. The target set includes the CH stretches and the modes that are strongly coupled with the CH stretches and/or involve strong Fermi resonances with a CH stretch fundamental. The rest of the modes are in the bath set. First, the effective harmonic oscillator (EHO) approximation is invoked for the whole system to obtain the zeroth-order frequencies and modified potentials. The effects of interaction between the bath set and the target sets are included in the modified potential from the EHO calculation. The VCCM equations are constructed with the modified potential from the EHO calculations and for the target set only. The transition energies and intensities are calculated using such a truncated VCCM approximation. The proposed method is applied to calculate the IR spectra of naphthalene and anthracene. The results with three different criteria for selecting the modes in the target set are compared with the experimental IR spectra.

Anharmonicity and the IR Emission Spectrum of Neutral Interstellar PAH Molecules

The characteristics of the CH stretching and out-of-plane bending modes in polycyclic aromatic hydrocarbon molecules are investigated using anharmonic density functional theory (DFT) coupled to a vibrational second-order perturbation treatment taking resonance effects into account. The results are used to calculate the infrared emission spectrum of vibrationally excited species in the collision-less environment of interstellar space. This model follows the energy cascade as the molecules relax after the absorption of a UV photon in order to calculate the detailed profiles of the infrared bands. The results are validated against elegant laboratory spectra of polycyclic aromatic hydrocarbon absorption and emission spectra obtained in molecular beams. The factors which influence the peak position, spectral detail, and relative strength of the CH stretching and out-of-plane bending modes are investigated, and detailed profiles for these modes are derived. These are compared to observations of astronomical objects in space, and the implications for our understanding of the characteristics of the molecular inventory of space are assessed.

Modeling the infrared cascade spectra of small PAHs: the 11.2 μm band

The profile of the 11.2 μm feature of the infrared (IR) cascade emission spectra of polycyclic aromatic hydrocarbon (PAH) molecules is investigated using a vibrational anharmonic method. Several factors are found to affect the profile including: the energy of the initially absorbed ultraviolet (UV) photon, the density of vibrational states, the anharmonic nature of the vibrational modes, the relative intensities of the vibrational modes, the rotational temperature of the molecule, and blending with nearby features. Each of these factors is explored independently and influence either the red or blue wing of the 11.2 μm feature. The majority impact solely the red wing, with the only factor altering the blue wing being the rotational temperature.

Predicting Infrared Spectra with Message Passing Neural Networks

Infrared (IR) spectroscopy remains an important tool for chemical characterization and identification. Chemprop-IR has been developed as a software package for the prediction of IR spectra through the use of machine learning. This work serves the dual purpose of providing a trained general-purpose model for the prediction of IR spectra with ease and providing the Chemprop-IR software framework for the training of new models. In Chemprop-IR, molecules are encoded using a directed message passing neural network, allowing for molecule latent representations to be learned and optimized for the task of spectral predictions. Model training incorporates spectra metrics and normalization techniques that offer better performance with spectral predictions than standard practice in regression models. The model makes use of pretraining using quantum chemistry calculations and ensembling of multiple submodels to improve generalizability and performance. The spectral predictions that result are of high quality, showing capability to capture the extreme diversity of spectral forms over chemical space and represent complex peak structures.

Understanding of the C-H stretch region of infra-red spectroscopy: an analysis of the final state wavefunctions

The nature of the wavefunctions associated with the final states in the CH stretching region of several medium sized molecules is analysed. The number of optically bright transitions is much larger than the number of CH oscillators present in the molecule, and they are spread over a range of about 300 cm^-1. Several of them are clustered together within about 5 cm^-1 with near equal intensities. The final states of all these transitions are superpositions of multiple zeroth order states. In almost all of such superpositions, no single zeroth order state has more than 50% weight. Several multiquantum states, with three to four quanta of excitation dominate the final states, with the CH chromophore contributing only a small weightage. Thus the band structure of the CH stretch region is due to several optically bright transitions whose final states are superpositions of low frequency multiquantum states with the CH chromophore contributing only a small weight to make them spectroscopically active.

A Spectroscopic View on Cosmic PAH Emission

ConspectusPolycyclic aromatic hydrocarbon molecules (PAHs) are ubiquitously present at high abundances in the Universe. They are detected through their infrared (IR) fluorescence UV pumped by nearby massive stars. Hence, their infrared emission is used to determine the star formation rate in galaxies, one of the key indicators for understanding the evolution of galaxies. Together with fullerenes, PAHs are the largest molecules found in space. They significantly partake in a variety of physical and chemical processes in space, influencing star and planet formation as well as galaxy evolution.Since the IR features from PAHs originate from chemical bonds involving only nearest neighbor atoms, they have only a weak dependence on the size and structure of the molecule, and it is therefore not possible to identify the individual PAH molecules that make up the cosmic PAH family. This strongly hampers the interpretation of their astronomical fingerprints. Despite the lack of identification, constraints can be set on the characteristics of the cosmic PAH family thanks to a joint effort of astronomers, physicists, and chemists.This Account presents the spectroscopic properties of the cosmic PAH emission as well as the intrinsic spectroscopic properties of PAHs and astronomical modeling of the PAH evolution required for the interpretation of the cosmic PAH characteristics. We discuss the observed spectral signatures tracing PAH properties such as charge, size, and structure and highlight the related challenges. We discuss the recent success of anharmonic calculations of PAH infrared absorption and emission spectra and outline the path forward. Finally, we illustrate the importance of models on PAH processing for the interpretation of the astronomical data in terms of the charge balance and PAH destruction.Throughout this Account, we emphasize that huge progress is on the horizon on the astronomical front. Indeed, the world is eagerly awaiting the launch of the James Webb Space Telescope (JWST). With its incredible improvement in spatial resolution, combined with its complete spectral coverage of the PAH infrared emission bands at medium spectral resolution and superb sensitivity, the JWST will revolutionize PAH research. Previous observations could only present spectra averaged over regions with vastly different properties, thus greatly confusing their interpretation. The amazing spatial resolution of JWST will disentangle these different regions. This will allow us to quantify precisely how PAHs are modified by the physical conditions of their host environment and thus trace how PAHs evolve across space. However, this will only be achieved when the necessary (and still missing) fundamental properties of PAHs, outlined in this Account, are known. We strongly encourage you to join this effort.

Anharmonic Infrared Spectra of Thermally Excited Pyrene (C16H10): A Combined View of DFT-Based GVPT2 with AnharmonicCaOs, and Approximate DFT Molecular dynamics with DemonNano

The study of the Aromatic Infrared Bands (AIBs) in astronomical environments has opened interesting spectroscopic questions on the effect of anharmonicity on the infrared (IR) spectrum of hot polycyclic aromatic hydrocarbons (PAHs) and related species in isolated conditions. The forthcoming James Webb Space Telescope will unveil unprecedented spatial and spectral details in the AIB spectrum; significant advancement is thus necessary now to model the infrared emission of PAHs, their presumed carriers, with enough detail to exploit the information content of the AIBs. This requires including anharmonicity in such models, and to do so systematically for all species included, requiring a difficult compromise between accuracy and efficiency. We performed a benchmark study to compare the performances of two methods in calculating anharmonic spectra, comparing them to available experimental data. One is a full quantum method, AnharmoniCaOs, relying on an ab initio potential, and the other relies on Molecular Dynamics simulations using a Density Functional based Tight Binding potential. The first one is found to be very accurate and detailed, but it becomes computationally very expensive for increasing temperature; the second is faster and can be used for arbitrarily high temperatures, but is less accurate. Still, its results can be used to model the evolution with temperature of isolated bands. We propose a new recipe to model anharmonic AIB emission using minimal assumptions on the general behaviour of band positions and widths with temperature, which can be defined by a small number of empirical parameters. Modelling accuracy will depend critically on these empirical parameters, allowing for an incremental improvement in model results, as better estimates become gradually available.

How to VPT2: Accurate and Intuitive Simulations of CH Stretching Infrared Spectra Using VPT2+K with Large Effective Hamiltonian Resonance Treatments

This article primarily discusses the utility of vibrational perturbation theory for the prediction of X-H stretching vibrations with particular focus on the specific variant, second-order vibrational perturbation theory with resonances (VPT2+K). It is written as a tutorial, reprinting most important formulas and providing numerous simple examples. It discusses the philosophy and practical considerations behind vibrational simulations with VPT2+K, including but not limited to computational method selection, cost-saving approximations, approaches to evaluating intensity, resonance identification, and effective Hamiltonian structure. Particular attention is given to resonance treatments, beginning with simple Fermi dyads and gradually progressing to arbitrarily large polyads that describe both Fermi and Darling-Dennison resonances. VPT2+K combined with large effective Hamiltonians is shown to be a reliable framework for modeling the complicated CH stretching spectra of alkenes. An error is also corrected in the published analytic formula for the VPT2 transition moment between the vibrational ground state and triply excited states.

Dependence upon charge of the vibrational spectra of small Polycyclic Aromatic Hydrocarbon clusters: the example of pyrene

Infrared spectra are computed for neutral and cationic clusters of Polycyclic Aromatic Hydrocarbon clusters, namely( C 16 H 10 ) n = 1 , 4 ( 0 / + ) , using the Density Functional based Tight Binding scheme combined with a Configuration Interaction scheme (DFTB-CI) in the double harmonic approximation. Cross-comparison is carried out with DFT and simple DFTB. Similarly to the monomer cation, the IR spectra of cluster cations are characterized by a depletion of the intensity of the CH stretch modes around 3000 cm-1, with a weak revival for n = 3 and 4. The in-plane CCC modes in the region 1400-2000 cm-1 are enhanced while the CH bending modes in the range 700-1000 cm-1 are significantly weakened with respect to the monomer cation, in particular for n = 2. Finally, soft modes corresponding to diedral fluctuations of the monomers within the central stack of the ion structure, possibly mixed with monomer folding, are also observed in the region 70-120 cm-1.

Anharmonic Frequencies and Spectroscopic Constants of OAlOH and AlOH: Strong Bonding but Unhindered Motion

The astrophysical buildup of premineral nanocrystals from atoms to the smallest network-covalent solids will require observations of various small molecules containing the most common elements in minerals including aluminum and oxygen. The present work utilizes high-level quantum chemical quartic force field (QFF) approaches to produce anharmonic vibrational frequencies and spectroscopic constants for such species. The computed B_eff for the astrochemically known AlOH molecule at 15780.5 MHz is a mere 40 MHz above the experimental value implying that the B_eff for OAlOH at 5580.9 MHz is similarly accurate. The additional 7.31 D dipole moment in OAlOH implies that this molecule is a viable target for interstellar observation. Unlike the other anharmonic vibrational frequencies reported in this work, the Al-O-H bending frequencies in both AlOH and OAlOH are poorly described in the present QFF results. However, this failing actually highlights the fact that these bends are exceptionally floppy yet with counterintuitive exceedingly strong bonding. The Al-O bond energies are 128.2 and 107.2 kcal/mol, respective of AlOH and OAlOH, while the barriers to linearity are meager 16.6 and 380.7 cm^-1 (0.1 and 1.1 kcal/mol).

Gas-Phase Infrared Spectroscopy of Neutral Peptides: Insights from the Far-IR and THz Domain

Gas-phase, double resonance IR spectroscopy has proven to be an excellent approach to obtain structural information on peptides ranging from single amino acids to large peptides and peptide clusters. In this review, we discuss the state-of-the-art of infrared action spectroscopy of peptides in the far-IR and THz regime. An introduction to the field of far-IR spectroscopy is given, thereby highlighting the opportunities that are provided for gas-phase research on neutral peptides. Current experimental methods, including spectroscopic schemes, have been reviewed. Structural information from the experimental far-IR spectra can be obtained with the help of suitable theoretical approaches such as dynamical DFT techniques and the recently developed Graph Theory. The aim of this review is to underline how the synergy between far-IR spectroscopy and theory can provide an unprecedented picture of the structure of neutral biomolecules in the gas phase. The far-IR signatures of the discussed studies are summarized in a far-IR map, in order to gain insight into the origin of the far-IR localized and delocalized motions present in peptides and where they can be found in the electromagnetic spectrum.

Anharmonic vibrational spectroscopy of polycyclic aromatic hydrocarbons (PAHs)

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 infrared transitions including bands becoming active due to resonances as well as overtone, 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₁₆H₁₀) and coronene (C₂₄H₁₂), 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 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 toward calculating moderately high-temperature spectra of PAHs and other similarly rigid molecules using Monte Carlo sampling.

The anharmonic quartic force field infrared spectra of hydrogenated and methylated PAHs

Reproducing and explaining the complex infrared spectra of functionalized polycyclic aromatic hydrocarbons through proper treatment of Fermi resonances.