The anharmonic quartic force field infrared spectra of three polycyclic aromatic hydrocarbons: Naphthalene, anthracene, and tetracene

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.

Size distribution of polycyclic aromatic hydrocarbons in space: an old new light on the 11.2/3.3 μm intensity ratio

The intensity ratio of the 11.2/3.3 μm emission bands is considered to be a reliable tracer of the size distribution of polycyclic aromatic hydrocarbons (PAHs) in the interstellar medium (ISM). This paper describes the validation of the calculated intrinsic infrared (IR) spectra of PAHs that underlie the interpretation of the observed ratio. The comparison of harmonic calculations from the NASA Ames PAH IR spectroscopic database to gas-phase experimental absorption IR spectra reveals a consistent underestimation of the 11.2/3.3 μm intensity ratio by 34%. IR spectra based on higher level anharmonic calculations, on the other hand, are in very good agreement with the experiments. While there are indications that the 11.2/3.3 μm ratio increases systematically for PAHs in the relevant size range when using a larger basis set, it is unfortunately not yet possible to reliably calculate anharmonic spectra for large PAHs. Based on these considerations, we have adjusted the intrinsic ratio of these modes and incorporated this in an interstellar PAH emission model. This corrected model implies that typical PAH sizes in reflection nebulae such as NGC 7023 - previously inferred to be in the range of 50 to 70 carbon atoms per PAH are actually in the range of 40 to 55 carbon atoms. The higher limit of this range is close to the size of the C60 fullerene (also detected in reflection nebulae), which would be in line with the hypothesis that, under appropriate conditions, large PAHs are converted into the more stable fullerenes in the ISM.

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.

Rotamers of Methanediol: Composite Ab Initio Predictions of Structures, Frequencies, and Rovibrational Constants

Geminal diols are known to be important intermediates in atmospheric ozonolysis and the aerosol cycle. Recently, the simplest member of this class, methanediol, was interrogated in the gas phase with infrared spectroscopy. To aid in future spectroscopic investigations of methanediol, including in the interstellar medium, we report fundamental frequencies and rovibrational constants for the two rotamers of this molecule using ab initio composite methods along with vibrational perturbation theory. Sensitivity of the predictions to the level of theory and the treatment of anharmonic resonances are carefully assessed. The OH stretching harmonic frequencies of both rotamers are particularly sensitive to the level of theory. The CH stretches of the C_s rotamer are sensitive to the treatment of anharmonic resonances with VPT2-based effective Hamiltonian models. Equilibrium bond distances and harmonic frequencies are converged conservatively to within 0.0005 Å and 3 cm^-1, respectively. The effect of tunneling on the rotational constants is investigated with a 2D variational calculation, based on a relaxed hydroxyl torsional potential energy surface. Tunneling is found to be negligible in the lower energy C₂ rotamer but should modify the rotational constants of the C_s rotamer on the order of MHz, giving rise to rotational line splittings of the same order. The rovibrational constants of the C_s rotamer are dominated by hydroxyl torsional effects, and here we see evidence for the breakdown of vibrational perturbation theory.

Aromatic and Acetylenic C-H or C-D Stretching Bands Anharmonicity Detection of Phenylacetylene by UV Laser-Induced Vibrational Emission

The anharmonic infrared (IR) emission spectra of phenylacetylene C₆H₅CCH and an isotopologue C₆H₅CCD induced by 193 nm UV excitation have been investigated in the gas phase. The study has been operated with a homemade IR spectrometer enabling to record time- and wavelength-resolved spectra between 2.5 and 4.5 μm, emitted all along the collisional cooling. The analysis is supported by a kinetic Monte Carlo simulation in the vibrational harmonic approximation. For both species, the anharmonic shifts of the acetylenic C-H or C-D stretching modes and the aromatic C-H stretching modes are studied for band positions and bandwidths in terms of the internal energy. For C₆H₅CCD, the internal energy dependence of the emission intensity band ratio is investigated and rationalized. This work demonstrates the potential of time-resolved IR emission spectroscopy to explore anharmonicity of astrophysically relevant molecules.

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.

F12-TZ-cCR: A Methodology for Faster and Still Highly Accurate Quartic Force Fields

The F12-TZ-cCR quartic force field (QFF) methodology, defined here as CCSD(T)-F12b/cc-pCVTZ-F12 with further corrections for relativity, is introduced as a cheaper and even more accurate alternative to more costly composite QFF methods like those containing complete basis set extrapolations within canonical coupled cluster theory. F12-TZ-cCR QFFs produce B₀ and C₀ vibrationally averaged principal rotational constants within 7.5 MHz of gas-phase experimental values for tetraatomic and larger molecules, offering higher accuracy in these constants than the previous composite methods. In addition, F12-TZ-cCR offers an order of magnitude decrease in the computational cost of highly accurate QFF methodologies accompanying this increase in accuracy. An additional order of magnitude in cost reduction is achieved in the F12-DZ-cCR method, while also matching the accuracy of the traditional composite method's B₀ and C₀ constants. Finally, F12-DZ and F12-TZ are benchmarked on the same test set, revealing that both methods can provide anharmonic vibrational frequencies that are comparable in accuracy to all three of the more expensive methodologies, although their rotational constants lag behind. Hence, the present work demonstrates that highly accurate theoretical rovibrational spectral data can be obtained for a fraction of the cost of conventional QFF methodologies, extending the applicability of QFFs to larger molecules.

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.

Infrared Spectroscopy of Chemically Diverse Carbon Clusters: A Data-Driven Approach

Carbon clusters exhibit a broad diversity of topologies and shapes, encompassing fullerene-like cages, graphene-like flakes, and more disordered pretzel-like and branched structures. Here, we examine computationally their infrared spectra in relation with these structures from a statistical perspective. Individual spectra for broad samples of isomers were determined by means of the self-consistent charge density functional-based tight-binding method, and an interpolation scheme is designed to reproduce the spectral features by regression on a much smaller subset of the sample. This interpolation proceeds by encoding the structures using appropriate descriptors and selecting them through principal component analysis, Gaussian regression or inverse distance weighting providing the nonlinear weighting functions. Metric learning is employed to reduce the global error on a preselected testing set. The interpolated spectra satisfactorily reproduce the specific spectral features and their dependence on the size and shape, enabling quantitative prediction away from the testing set. Finally, the classification of structures within the four proposed families is critically discussed through a statistical analysis of the sample based on iterative label spreading.

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.

Optical Cycling Functionalization of Arenes

Closed, laser-induced optical transitions ("optical cycling transitions") of molecules can be used for state preparation and measurement in quantum information science and quantum sensing. Increasingly complex molecular species supporting optical cycling can provide new capabilities for quantum science, and it is not clear if there is a limit on their size or complexity. We explore Ca-O-L molecular constructs to support the optical cycling center, Ca, with ligands, L, being arenes. We find that L can be as large as coronene (i.e., CaOC₂₄H₁₁) without losing the diagonality of the Franck-Condon factor (FCF). Furthermore, L can be substituted with electron-withdrawing groups to improve the FCF. Larger L, beyond ∼7 rings, can disrupt the diagonality of the FCF by closing the HOMO-LUMO ligand electronic state gap and reordering with the local states on the cycling center. Overall, we find that optical cycling can be retained for arenes, and we offer a principle for their design.

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.

Non long-range corrected density functionals incorrectly describe the intensity of the CH stretching band in polycyclic aromatic hydrocarbons

We present a comprehensive study of the most relevant numerical aspects influencing frequencies and intensities in the infrared spectrum of isolated polycyclic aromatic hydrocarbons (PAHs) regarding the overestimate of the IR CH-stretching bands. We use naphthalene as benchmark and show the validity of our results to different members of the PAH family. Our analysis relies on widely employed density functional theory methods and second-order vibrational perturbational theory for the computation of vibrational eigenstates. We have focused on the elucidation of the origin of the systematic overestimate of the intensities in the CH-stretching region. To rule out nonfundamental numerical errors, we have initially considered the influence of the electronic basis set and various other parameters on the different stages of the vibrational analysis. In a second stage, we have benchmarked the results of different density functional theory functionals with respect to the aforementioned overestimate taken as the ratio between the most prominent features of the spectrum, the CH-bending and the CH-stretching bands. Our results unambiguously indicate that the long-range correction plays a major role in this spurious numerical issue. More specifically, this phenomenon is due to an incorrect description of the charge distribution (and hence dipole) within the symmetrically relevant CH bonds. Long-range correction specifically remedies this issue. It improves the description of the intensities in the stretching region while at the same time it does not perturb significantly the rest of the spectrum. With respect to the frequencies, we have observed an overall improvement when compared to noncorrected functionals.

The unsolved issue with out-of-plane bending frequencies for CC multiply bonded systems

More than 30 years ago two groups independently identified a problem in the calculation of the out-of-plane bending (OPB) vibrational frequencies for the ethylene molecule using correlated electronic structure methods. Several studies have been done in the meantime to try and understand and resolve this issue. In so doing this problem has been found to be far more insidious than previously realized for acetylene-like and benzene-like molecules, which can become non-linear and non-planar, respectively. The one common feature that all molecules with this problem have is that they contain CC multiple bonds, and so this has been called the "CC multiple bond OPB frequency issue" or "the CC OPB problem." Various explanations for this problem have been advanced such as basis set superposition error, basis set incompleteness error, linear dependences in the basis set, proper balancing of the basis set between saturation and inclusion of higher angular momentum functions, etc. and possible solutions have arisen from these suggestions. All of these proposed solutions, however, amount to one main point connecting them all: modifying the one-particle basis set in some way. None of the explanations that have been advanced, however, really fit all of the data for all of the molecules where this problem has been identified, and importantly, none of these diagnostic tests have been applied to similar molecules where this issue does not appear. In this review, the studies over the last 30 plus years are discussed and relevant data from each of these is compared and contrasted. It is hoped that by collecting and analyzing the data from these studies a path forward to understanding and resolving this issue will become evident.

Correlated Triplet Pair Formation Activated by Geometry Relaxation in Directly Linked Tetracene Dimer (5,5'-Bitetracene)

Singlet fission (SF) materials have the potential to overcome the traditional external quantum efficiency limits of organic light-emitting diodes (OLEDs). In this study, we theoretically designed an intramolecular SF molecule, 5,5'-bitetracene (55BT), in which two tetracene units were directly connected through a C-C bond. Using quantum chemical calculation and the Fermi golden rule, we show that 55BT undergoes efficient SF induced by geometry relaxation in a locally excited singlet state, ¹(S₀S₁). Compared with another high-performing SF system, the tetracene dimer in the crystalline state, 55BT has advantages when used in doped systems owing to covalent bonding of the two tetracene units. This feature makes 55BT a promising candidate triplet sensitizer for near-infrared OLEDs.

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.

Rovibrational Quantum Chemical Treatment of Inorganic and Organometallic Astrochemicals

ConspectusOur two groups have both independently and collaboratively been pushing quantum-chemical techniques to produce highly accurate predictions of anharmonic vibrational frequencies and spectroscopic constants for molecules containing atoms outside of the typical upper p block. Methodologies employ composite approaches, relying on various levels of coupled cluster theory-most often at the singles, doubles, and perturbative triples level-and quartic force field constructions of the potential portion of the intramolecular Watson Hamiltonian. Such methods are known to perform well for organic species, and we have extended this to molecules containing atoms outside of this realm.One notable atom that has received much attention in this application is magnesium. Mg is the second-most-abundant element in the Earth's mantle, and while molecules containing this element are among the confirmed astrochemicals, its further atomic abundance in the galaxy implies that many more molecules (both purely inorganic and organometallic) containing element 12 exist in astrophysical regions in chemical sizes between those of atoms and dust-sized nanocrystals. Our approach discussed herein is producing quality benchmarks and predicting novel data for magnesium-bearing molecules.The story is similar for Al and Si, which are also notably abundant in both rocky bodies and the universe at large. While Na, Sc, and Cu may not be as abundant as Mg, Al, and Si, molecules containing Na and transition metals have also previously been reported to be detected beyond the Earth. Consequently, the need to produce spectral reference data for molecules containing such atoms is growing. While several experimental groups (including, notably, the groups in Arizona, Boston, and France/Spain) have clearly led the way in detection of inorganic/organometallic molecules in space, computational support and even rational design can provide novel avenues for the detection of molecules containing atoms not typically studied in most laboratories. The application of quantum chemistry to other elements beyond carbon and its cronies at the top right of the periodic table promises a better understanding of the observable universe. It will also provide novel and fundamental chemical insights pushing the "central science" into new molecular territory.

Structural investigation of doubly-dehydrogenated pyrene cations

The vibrationally resolved spectra of the pyrene cation and doubly-dehydrogenated pyrene cation (C₁₆H₁₀˙⁺; Py⁺ and C₁₆H₈˙⁺; ddPy⁺) are presented. Infrared predissociation spectroscopy is employed to measure the vibrational spectrum of both species using a cryogenically cooled 22-pole ion trap. The spectrum of Py⁺ allows a detailed comparison with harmonic and anharmonic density functional theory (DFT) calculated normal mode frequencies. The spectrum of ddPy⁺ is dominated by absorption features from two isomers (4,5-ddPy⁺ and 1,2-ddPy⁺) with, at most, minor contributions from other isomers. These findings can be extended to explore the release of hydrogen from interstellar PAH species. Our results suggest that this process favours the loss of adjacent hydrogen atoms.

Efficient Analytic Second Derivative of Electrostatic Embedding QM/MM Energy: Normal Mode Analysis of Plant Cryptochrome

Analytic second derivatives of electrostatic embedding (EE) quantum mechanics/molecular mechanics (QM/MM) energy are important for performing vibrational analysis and simulating vibrational spectra of quantum systems interacting with an environment represented as a classical electrostatic potential. The main bottleneck of EE-QM/MM second derivatives is the solution of coupled perturbed equations for each MM atom perturbation. Here, we exploit the Q-vector method [J. Chem. Phys., 2019, 151, 041102] to workaround this bottleneck. We derive the full analytic second derivative of the EE-QM/MM energy, which allows us to compute QM, MM, and QM-MM Hessian blocks in an efficient and easy to implement manner. To show the capabilities of our method, we compute the normal modes for the full Arabidopsis thaliana plant cryptochrome. We show that the flavin adenine dinucleotide vibrations (QM subsystem) strongly mix with protein modes. We compute approximate vibronic couplings for the lowest bright transition, from which we extract spectral densities and the homogeneous broadening of FAD absorption spectrum in protein using vibrationally resolved electronic spectrum simulations.

The Production and Potential Detection of Hexamethylenetetramine-Methanol in Space

Numerous laboratory studies of astrophysical ice analogues have shown that their exposure to ionizing radiation leads to the production of large numbers of new, more complex compounds, many of which are of astrobiological interest. We show here that the irradiation of astrophysical ice analogues containing H₂O, CH₃OH, CO, and NH₃ yields quantities of hexamethylenetetramine-methanol (hereafter HMT-methanol; C₇N₄H₁₄O) that are easily detectible in the resulting organic residues. This molecule differs from simple HMT, which is known to be abundant in similar ice photolysis residues, by the replacement of a peripheral H atom with a CH₂OH group. As with HMT, HMT-methanol is likely to be an amino acid precursor. HMT has tetrahedral (T_d) symmetry, whereas HMT-methanol has C₁ symmetry. We report the computed expected infrared spectra for HMT and HMT-methanol obtained using ab initio quantum chemistry methods and show that there is a good match between the observed and computed spectra for regular HMT. Since HMT-methanol lacks the high symmetry of HMT, it produces rotational transitions that could be observed at longer wavelengths, although establishing the exact positions of these transitions may be challenging. It is likely that HMT-methanol represents an abundant member of a larger family of functionalized HMT molecules that may be present in cold astrophysical environments.

Experimental Approach to the Study of Anharmonicity in the Infrared Spectrum of Pyrene from 14 to 723 K

Quantifying the effect of anharmonicity on the infrared spectrum of large molecules such as polycyclic aromatic hydrocarbons (PAHs) at high temperatures is the focus of a number of theoretical and experimental studies, many of them motivated by astrophysical applications. We recorded the IR spectrum of pyrene C16H10 microcrystals embedded in KBr pellets over a wide range of temperatures (14-723 K) and studied the evolution of band positions, widths, and integrated intensities with temperature. We identified jumps for some of the spectral characteristics of some bands in the 423-473 K range. These were attributed to a change of phase from crystal to molten in condensed pyrene, which appears to affect more strongly bands involving large CH motions. Empirical anharmonicity factors that quantify the linear evolution of band positions and widths with temperature for values larger than ∼150-250 K, depending on the band, were retrieved from both phases and averaged to provide recommended values for these anharmonicity factors. The derived values were found to be consistent with available gas phase data. We conclude about the relevance of the methodology to produce data that can be compared with calculated anharmonic IR spectra and provide input for models that simulate the IR emission of astro-PAHs.

Tracking down the origin of peculiar vibrational spectra of aromatic self-assembled thiolate monolayers

Several studies have previously observed surprisingly low frequencies for the C-H stretching modes of self-assembled monolayers (SAMs) prepared from aromatic thiols. The reason for this property has so far remained elusive. Therefore, we report a novel study of the vibrational spectra of SAMs prepared on Au from two different aromatic thiols, namely, 4'-nitro-1,1'-biphenyl-4-thiol (NBPT) and 4-aminothiophenol (ATP). The SAMs were prepared by vapor deposition (VD) in ultrahigh vacuum (UHV) as well as by the solution method (SM) and their quality was controlled by X-ray photoelectron spectroscopy (XPS). In addition, amino terminated SAMs were also obtained by electron irradiation and by chemical reduction of NBPT SAMs. Beside infrared reflection absorption spectroscopy (IRRAS), we have employed high resolution electron energy loss spectroscopy (HREELS), by which VD SAMs can be studied in situ, i.e. without exposing them to air. Hence, we can exclude possible contributions of solvent molecules to the vibrational spectra. Nonetheless, HREELS in fact reveals the same large red shift of the C-H stretching modes in the SAMs as also observed in ex situ IRRAS experiments. In contrast, HREELS for physisorbed ATP and ATP in a KBr pellet measured by transmission infrared spectroscopy exhibit the expected aromatic bands. Using a computational approach, we can exclude molecular packing effects as origin of this shift. Therefore, we propose chemical changes in the aromatic rings during SAM formation as an alternative explanation for the observed frequency shift. As another striking effect, the N-H stretching vibrational modes of the amino-terminated SAMs are extremely weak in both IRRAS and HREELS despite the fact that XPS confirms the presence of amino groups. A very weak signal is observed only in the case of an electron irradiated NBPT SAM. In contrast, an energy loss ascribed to the N-H stretching vibrations is clearly observed in HREELS of ATP physisorbed on an ATP SAM and on graphite as well as in the transmission infrared spectrum of ATP in KBr. The extremely low intensity of these vibrations in the SAM is traced back to the inherently low transition dipole moment for the excitation of N-H stretching modes in free N-H groups. Furthermore, the calculations suggest that the much stronger signals of N-H stretching modes involved in hydrogen-bonding with adjacent amino groups are suppressed because these vibrations are oriented parallel to the surface.

The Astrochemistry Implications of Quantum Chemical Normal Modes Vibrational Analysis

Understanding the molecular vibrations underlying each of the unknown infrared emission (UIE) bands (such as those found at 3.3, 3.4, 3.5, 6.2, 6.9, 7.7, 11.3, 15.8, 16.4, 18.9 μ m) observed in or towards astronomical objects is a vital link to uncover the molecular identity of their carriers. This is usually done by customary classifications of normal-mode frequencies such as stretching, deformation, rocking, wagging, skeletal mode, etc. A large literature on this subject exists and since 1952 ambiguities in classifications of normal modes via this empirical approach were pointed out by Morino and Kuchitsu New ways of interpretation and analyzing vibrational spectra were sought within the theoretical framework of quantum chemistry. Many of these methods cannot easily be applied to the large, complex molecular systems which are one of the key research interests of astrochemistry. In considering this demand, a simple and new method of analyzing and classifying the normal mode vibrational motions of molecular systems was introduced. This approach is a fully quantitative method of analysis of normal-mode displacement vector matrices and classification of the characteristic frequencies (fundamentals) underlying the observed IR bands. Outcomes of applying such an approach show some overlap with customary empirical classifications, usually at short wavelengths. It provides a quantitative breakdown of a complex vibration (at longer wavelengths) into the contributed fragments such as their aromatic or aliphatic components. In addition, in molecular systems outside the classical models of chemical bonds and structures where the empirical approach cannot be applied, this quantitative method enables an interpretation of vibrational motion(s) underlying the IR bands. As a result, further modifications in the structures (modeling) and the generation of the IR spectra (simulating) of the UIE carriers, initiated by proposing a PAH model, can be implemented in an efficient way. Here fresh results on the vibrational origin of the spectacular UIE bands based on astrochemistry molecular models, explored through the lens of the quantitative method applied to thousands of different vibrational motion matrices are discussed. These results are important in the context of protoplanetary nebulae and planetary nebulae where various molecular species have been uncovered despite their harsh environments.

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.

Direct Evidence of the Benzylium and Tropylium Cations as the Two Long-Lived Isomers of C7 H7. Chemphyschem 2018;19:3182-3185. [PMID: 30238585 PMCID: PMC6420061 DOI: 10.1002/cphc.201800744]

Disentangling the isomeric structure of C7 H7 + is a longstanding experimental issue. We report here the full mid-infrared vibrational spectrum of C7 H7 + tagged with Ne obtained with infrared-predissociation spectroscopy at 10 K. Saturation depletion measurements were used to assign the contribution of benzylium and tropylium isomers and demonstrate that no other isomer is involved. Recorded spectral features compare well with density functional theory calculations. This opens perspectives for a better understanding and control of the formation paths leading to either tropylium or benzylium ions.

Photoelectron spectroscopy of anthracene and fluoranthene radical anions

We report the slow electron velocity map imaging spectroscopy of cryogenically cooled anthracene and fluoranthene radical anions, two similarly sized polycyclic aromatic hydrocarbon molecules. The results allow us to examine the lowest energy singlet and triplet states in the neutral molecules on equal footing from the anionic ground state. The analysis of the experimental spectra is aided by harmonic calculations and Franck-Condon simulations, which generally show good agreement with experimental values and spectra. The electron affinity of fluoranthene is measured to be 0.757(2) eV, which is larger than that of anthracene at 0.532(3) eV. The lowest energy triplet state in anthracene is observed at 1.872(3) eV above the singlet ground state, while that of fluoranthene is observed at 2.321(2) eV above its singlet ground state. Comparisons of experimental and calculated spectra show that in addition to the Franck-Condon active modes, there is a clear presence of vibrational modes that gain intensity via vibronic coupling in both the singlet and triplet states in both molecules. In addition, the triplet state generally exhibits increased vibronic coupling compared to the singlet state, with the fluoranthene triplet state exhibiting evidence of distortion from C_2v symmetry.

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.