Anharmonic Vibrational Frequencies Research Articles

On the Formation and Detectability of H2CNCN and Its Progenitors

New highly exothermic formation pathways incorporating both thermodynamic and kinetic control for the newly astronomically detected H2CNCN molecule are paired with extremely accurate quantum chemical rovibrational spectroscopic computations. The reactions between astronomically known CH2CN/CH2CCH + HNCN follow effectively identical pathways and proceed through stable intermediates and over deeply submerged transition states to form H2CNCN and HCN/HCCH coproducts. Similarly, the reaction between CH2CN and NCN− can also form H2CNCN, although this pathway first requires the initial formation of NCN−, which is currently undetected in space, via HNCN + CN−. This two-step mechanism uses the highly abundant CN− as the catalyst. Incredibly accurate quantum chemical spectroscopic data are reported for all reactants and products of these reactions, with errors between experimental values and the computations herein on the order of 0.1% or less. Anharmonic vibrational frequencies and intensities are also reported in order to guide experimental and observational searches for these molecules that have mostly been detected in the radio but may now be detectable via JWST.

Anharmonic Vibrational Spectroscopy of Germanium-Containing Clusters, GexC4-x and GexSi4-x (x = 0-4), for Interstellar Detection.

An extensive, high-level theoretical study on tetra-atomic germanium carbide/silicide clusters is presented. Accurate harmonic and anharmonic vibrational frequencies and rotational constants are calculated at the CCSD(T)-F12a(b)/cc-pVT(Q)Z-F12 levels of theory. With growing capabilities to discern more of the chemical composition of the interstellar medium (ISM), an accurate database of reference material is required. The presence of carbon is ubiquitous in the ISM, and silicon is known to be present in interstellar dust grains; however, germanium-containing molecules remain elusive. To begin understanding the presence and role of germanium in the ISM, we present this study of the vibrational and rotational spectroscopic properties of various germanium-containing molecules to aid in their potential identification in the ISM with modern observational tools such as the James Webb Space Telescope. Structures studied herein include rhomboidal (r-), diamond (d-), and trapezoidal (t-) tetra-atomic molecules of the form GexC4-x and GexSi4-x, where x = 0-4. The most promising structure for detection is r-Ge2C2 via the ν4 mode with a frequency of 802.7 cm-1 (12.5 μm) and an intensity of 307.2 km mol-1. Other molecules that are potentially detectable, i.e., through vibrational modes or rotational transitions, include r-Ge3C, r-GeSi3, d-GeC3, r-GeC3, and t-Ge2C2.

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.

Molecular Hydrogen Acts as a Hydrogen Bond Proton Acceptor: From Protonated Betaine Tagging to the Weakest Hydrogen Bond.

In an attempt to gain further insights into the intermolecular interactions implied by Rizzo's group's cautionary tale related to molecular tagging in infrared multiple photon dissociation (IRMPD) spectroscopy with molecular messengers [Masson, A. . J. Chem. Phys. 2015, 143, 104313], in the present study, we provide an in-depth analysis of the noncovalent interaction between the molecular hydrogen and protonated betaine molecule in the gas phase. We aim to shed some new light on the fundamental issues concerning the wide diapason of hydrogen-bonding-type intermolecular interactions, with a wide variety of proton acceptors. We demonstrate that in the course of tagging the protonated betaine with molecular hydrogen from the OH group side, it is the σ bond of molecular hydrogen that plays the role of hydrogen-bonding proton acceptor. The tagging thus induces a small yet significant red shift of the protonated betaine O-H stretching mode. We investigate the performance of a wide range of density functional theory (DFT) functionals for the calculation of anharmonic vibrational frequency shifts of the studied system, which are essential for the correct interpretation of the experimental IRMPD data. For an accurate prediction of the OH stretching frequency shifts, specifically designed functionals such as Handy's group HCTH/407 should be applied. The empirical dispersion correction enhances the systematic overestimation of the anharmonic frequency shift, characteristic of the most widely used DFT functionals. Combining the full-wave function approach with the charge field perturbation and natural bond orbital (NBO) deletion analyses, we demonstrate that the frequency shift in the OH-tagged structure is governed by the σHH → σ*OH intermolecular charge transfer. This interaction stabilizes the OH-tagged dimer as well, in contrast to the dipole-quadrupole electrostatic interaction energy term. Topological analysis of the electron density reveals the presence of an intermolecular bond critical point with a positive value of the density Laplacian very close to the lower limit for hydrogen bonds. NCI analyses demonstrate that the OH···H2 interaction is weaker than the intramolecular CH···O one within the protonated betaine molecule, with the through of reduced density gradient appearing at less negative sign(λ2)·ρ values. Analyzing the O-H stretching vibrational potential with the second-generation absolutely localized molecular orbitals energy decomposition analysis (ALMO-EDA 2) revealed that in the case of betaineH(+) tagged from the OH group side, the permanent electrostatics (ΔEelec), polarization (ΔEpol), and charge-transfer (ΔEct) contributions to the total intermolecular interaction energy contribute favorably to the weak hydrogen bond formation and to the red shift of the fundamental O-H stretching frequency, the ΔEct contribution being the most significant in the last context. The Pauli repulsion term, on the other hand, favors an O-H stretching frequency blue shift as a consequence of the vibrational confinement effects.

Performance of EOM-CCSD(T)(a)*-Based Quartic Force Fields in Computing Fundamental, Anharmonic Vibrational Frequencies of Molecular Electronically Excited States with Application to the Ã1A″ State of :CCH2 (Vinylidene).

Highly accurate anharmonic vibrational frequencies of electronically excited states are not as easily computed as their ground electronic state counterparts, but recently developed approximate triple excited state methods may be changing that. One emerging excited state method is equation of motion coupled cluster theory at the singles and doubles level with perturbative triples computed via the (a)* formalism, EOMEE-CCSD(T)(a)*. One of the most employed means for the ready computation of vibrational anharmonic frequencies for ground electronic states is second-order vibrational perturbation theory (VPT2), a theory based on quartic force fields (QFFs),fourth-order Taylor series expansions of the potential portion of the internuclear Watson Hamiltonian. The combination of these two is herein benchmarked for its performance for use as a means of computing rovibrational spectra of electronically excited states. Specifically, the EOMEE-CCSD(T)(a)* approach employing a complete basis set extrapolation along with core electron inclusion and relativity (the so-called "CcCR" approach) defining the QFF produces anharmonic fundamental vibrational frequencies within 2.83%, on the average, of reported gas-phase experimentally assigned values for the test set including the states of HCF, HCCl, HSiF, HNO, and HPO. However, some states have exceptional accuracy in the fundamentals, most notably for ν2 of HCCl in which the CcCR QFF value is within 1.8 cm-1 at 927.9 cm-1 (or 0.2%) of the experiment. Additionally, this approach produces rotational constants to, on the absolute average, within 0.41% of available experimental data, showcasing notable accuracy in the computation of rovibronic spectral data. Furthermore, utilizing a hybrid approach composed of harmonic CcCR force constants along with a set of simple EOMEE-CCSD(T)(a)*/aug-cc-pVQZ QFF cubic and quartic force constants is faster than using pure CcCR and better represents those modes that suffer from numerical instability in the anharmonic portion of the QFF, implying that this so-called "CcCR + QZ" QFF approach may be the best for future applications. Finally, complete, rovibrational spectral data are provided for :CCH2, a molecule of potential astrochemical interest, in order to aid in its potential future experimental rovibronic characterization.

Diborane anharmonic vibrational frequencies and Intensities: Experiment and theory

The vibrational dynamics of diborane have been extensively studied both theoretically and experimentally ever since the bridge structure of diborane was established in the 1950s. Numerous infrared and several Raman spectroscopic studies have followed in the ensuing years at ever increasing levels of spectral resolution. In parallel, ab initio computations of the underlying potential energy surface have progressed as well as the methods to calculate the anharmonic vibration dynamics beyond the double harmonic approximation. Nevertheless, even 70 years after the bridge structure of diborane was established, there are still significant discrepancies between experiment and theory for the fundamental vibrational frequencies of diborane. In this work we use parahydrogen (pH2) matrix isolation infrared spectroscopy to characterize six fundamental vibrations of B2H6 and B2D6 and compare them with results from configuration-selective vibrational configuration interaction theory. The calculated frequencies and intensities are in very good agreement with the pH2 matrix isolation spectra, even several combination bands are well reproduced. We believe that the reason discrepancies have existed for so long is related to the large amount of anharmonicity that is associated with the bridge BH stretching modes. However, the calculated frequencies and intensities reported here for the vibrational modes of all three boron isotopologues of B2H6 and B2D6 are within ± 2.00 cm−1 and ± 1.44 cm−1, respectively, of the experimental frequencies and therefore a refined vibrational assignment of diborane has been achieved.

Revisiting the vibrational spectrum of formic acid anhydride

The potential energy surface of formic acid anhydride has been investigated at highly accurate quantum chemical methods. The rotation of CHO group in global minimum conformer of formic acid anhydride can lead two different local minimum conformers. Both local minimum conformers are around three kcal/mol higher in energy than the global minimum conformer. One conformer is planar and another has two rotamers which rapidly interconvert to each others via planar transition state. In some earlier studies, this planar transition state has incorrectly assigned to be a local minimum structure. We calculated anharmonic vibrational frequencies for the global and local minimum energy conformers and compared the theoretical wavenumbers with experimental gas phase and argon matrix measurements. Our results suggest that some of the experimentally detected peaks are overtone and combination bands. In previous studies with harmonic calculations, those peaks have assigned to be fundamental bands with zero intensities. We confirmed that the higher energy conformer produced in argon matrix by ultra violet induced rotamerization of global minimum conformer belongs to the planar conformer.

DFT + F12 QFFs for Cost-Effective Rovibrational Spectral Data Predictions of Ground and Excited Electronic States.

The quest for faster computation of anharmonic vibrational frequencies of both ground and excited electronic states has led to combining coupled cluster theory harmonic force constants with density functional theory cubic and quartic force constants for defining a quartic force field (QFF) utilized in conjunction with vibrational perturbation theory at second order (VPT2). This work shows that explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples levels [CCSD(T)-F12] provides accurate anharmonic vibrational frequencies and rotational constants when conjoined with any of B3LYP, CAM-B3LYP, BHandHLYP, PBE0, and ωB97XD for roughly one-quarter of the computational time of the CCSD(T)-F12 QFF alone for our test set. As the number of atoms in the molecule increases, however, the anharmonic terms become a greater portion of the QFF, and the cost comparison improves with HOCO+ and formic acid, requiring less than 15 and 10% of the time, respectively. In electronically excited states, PBE0 produces more consistently accurate results. Additionally, as the size of the molecule and, in turn, QFF increase, the cost savings for utilizing such a hybrid approach for both ground- and excited-state computations grows. As such, these methods are promising for predicting accurate rovibrational spectral properties for electronically excited states. In cases where well-behaved potentials for a small selection of targeted excited states are needed, such an approach should reduce the computational cost compared to that of methods requiring semiglobal potential surfaces or variational treatments of the rovibronic Hamiltonian. Such applications include spectral characterization of comets, exoplanets, or any situation in which gas phase molecules are being excited by UV-vis radiation.

Exploring the Tl H potential energy surface: A comparative analysis with group 13 systems and experiment.

Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet Tl H potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug-cc-pwCVQZ-PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero-point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 M H (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist.

Lowest electronic states of neutral and ionic LiN

AbstractWe have investigated the potential energy curves (PECs) of the LiN heteronuclear diatomic molecule, including its ionic species LiN+ and LiN−, using explicitly correlated multi‐reference configuration interaction (MRCI‐F12) calculations in conjunction with the correlation consistent quintuple‐𝜁 basis set. The effect of core–valence correlation, scalar relativistic effects, and the size of the basis sets has been investigated. A comprehensive set of spectroscopic constants determined based on the above‐mentioned calculations are also reported for the lowest electronic states and all systems, including dissociation energies, harmonic and anharmonic vibrational frequencies, and rotational constants. Additional parameters, such as the dipole moments, equilibrium spin‐orbit constants, excitation energies, and rovibrational energy levels, are also documented. We found that the three triplet states of LiN, namely, X 3∑−, A 3Π, and 2 3∑−, exhibit substantial potential wells in the PEC diagrams, while the quintet states are repulsive in nature. The ground state of the anion also shows a deep potential well in the vicinity of its equilibrium geometry. In contrast, the ground and excited states of the cation are very loosely bound. Charge transfer properties of each of these states are also analyzed to obtain an in‐depth understanding of the interatomic interactions. We found that the core–valence correlation has a substantial effect on the calculated spectroscopic constants.

Taming semi-empirical methods for PAHs and vibrational spectra

Semi-empirical methods offer a cost-effective means of computing explicit, anharmonic vibrational frequencies for large molecules, such as polycyclic aromatic hydrocarbons (PAHs), but their default parameters produce insufficiently accurate results for comparison to experiment, especially in the hydride stretching region where the NIRSpec instrument on JWST is most effective. This work delivers several reparameterized variants of the PM6 semi-empirical method trained to reproduce the experimental vibrational frequencies of 5 small hydrocarbon molecules. When benchmarked on the experimental frequencies of benzene and naphthalene, these reparameterized methods match the empirical values to within 38.7 cm−1 on average. Combining these values with the default PM6 frequencies below 3000 cm−1 brings the average deviation below 22 cm−1 for naphthalene, comparing favorably with the existing state of the art in B3LYP/4-31G, for a two order of magnitude decrease in the computational cost. As such, the present work offers a promising means of extending the computation of explicit, anharmonic vibrational frequencies to PAHs larger than those that can be examined anharmonically via B3LYP. Such large and accurate data sets are necessary to disentangle the unidentified spectral features observed around myriad astronomical bodies and the influx of observational data from JWST.

Ab initio calculations on the [formula omitted] → [formula omitted] photoelectron spectrum of thioformaldehyde negative ion

The equilibrium structures, multidimensional potential energy surfaces, and anharmonic vibrational frequencies of thioformaldehyde (H2CS) and its 2B1 radical anion were obtained at the (U)CCSD(T)-F12a/cc-pVTZ-F12 level; and the photoelectron spectrum of H2CS−X∼2B1→ H2CS X∼1A1 was simulated by calculating the Franck-Condon factors (FCFs). The additional electrons cause a significant change in the bond length of the C = S bond of H2CS, which affects the potential energy surface and the C = S bond stretching mode, while the interaction between electrons turns the originally planar structure of H2CS into a slightly out-of-plane umbrella structure, resulting in an unusual change in the vibrational frequency of the out-of-plane bending mode. Finally, a time-independent eigenstate-free Raman wave function approach (RWF) was used to calculate the photoelectron spectrum of H2CS and to point out the changes in the band under the influence of different temperatures. These data provides a theoretical basis for interstellar observations and experimental studies.

The performance of hybrid and F12∗/F12c explicitly correlated coupled cluster methods for use in anharmonic vibrational frequency computations

AbstractA hybrid quartic force field approach produces the same accuracies as non‐hybrid methods but for less than one quarter of the computational time. This method utilizes explicitly correlated coupled cluster theory at the singles and doubles level inclusive of perturbative triples (CCSD(T)‐F12b) in conjunction with a triple‐ basis set, core electron correlation, and scalar relativity for the harmonic terms and CCSD(T)‐F12b with a valence double‐ basis set for the cubic and quartic terms. There is no sacrifice in the prediction of fundamental anharmonic vibrational frequencies or vibrationally‐averaged rotational constants as compared to experiment, but the time saved is notable. Other hybrid methods are examined involving different sizes of basis sets and composite terms included or excluded. Not one is more accurate; only one is faster. F12 (also called F12c) is tested as well, but it has an increase in computational time for no increase in accuracy. As such, this work reports a hybrid and composite approach (F12‐TcCR+DZ) in the computation of rovibrational spectral data which can be applied to the observation of novel molecules in the gas phase in the laboratory and potentially even in astrophysical environments.

Molecular Structures and Vibrational Spectra of trans- and cis-Polyacetylene and Their Oligoenes Revisited Using Density Functional Theory Calculations.

Polyacetylene, the most representative synthetic conducting polymer, has attracted much attention because it exhibits high conductivity upon doping. In this paper, molecular structures, electronic excitation energies, and Raman and infrared spectra were calculated using density functional theory for trans- and cis-oligoenes with various chain lengths up to the number of C═C bonds (n) of 100 and trans- and cis-polyacetylenes under one-dimensional periodic boundary condition. The harmonic vibrational frequencies obtained at the B3LYP/6-311G(d,p) level were scaled by the scaling factors determined with respect to the anharmonic vibrational frequencies using the B2PLYP method, in which the coefficients of the functional were optimized for trans-oligoenes. The calculated infrared and Raman frequencies reproduce reasonably well the observed frequencies for trans- and cis-polyacetylene. Based on the chain-length dependence of the calculated Raman spectra of trans-oligoenes, we proposed the possibility of longer conjugated trans-segments observed in the resonance Raman spectra of trans-polyacetylene excited at longer wavelengths of 647.1 and 1064 nm. We also elucidated the origin of the excitation-wavelength dependence of the resonance Raman spectra of trans-polyacetylene and the structure of isomerization intermediates from cis-form to trans-form. In addition, the previous assignments of Raman and infrared spectra of trans- and cis-polyacetylene were reexamined in the present study based on the chain-length dependence of the spectra.

F12+EOM Quartic Force Fields for Rovibrational Predictions of Electronically Excited States.

Quartic force fields (QFFs) constructed using a sum of ground-state CCSD(T)-F12b energies with EOM-CCSD excitation energies are proposed for computation of spectroscopic properties of electronically excited states. This is dubbed the F12+EOM approach and is shown to provide similar accuracy to previous methodologies at lower computational cost. Using explicitly correlated F12 approaches instead of canonical CCSD(T), as in the corresponding (T)+EOM approach, allows for 70-fold improvement in computational time. The mean percent difference between the two methods for anharmonic vibrational frequencies is only 0.10%. A similar approach is also developed herein which accounts for core correlation and scalar relativistic effects, named F12cCR+EOM. The F12+EOM and F12cCR+EOM approaches both match to within 2.5% mean absolute error of experimental fundamental frequencies. These new methods should help in clarifying astronomical spectra by assigning features to vibronic and vibrational transitions of small astromolecules when such data are not available experimentally.

Positioning of grid points for spanning potential energy surfaces-How much effort is really needed?

The positions of grid points for representing a multidimensional potential energy surface (PES) have a non-negligible impact on its accuracy and the associated computational effort for its generation. Six different positioning schemes were studied for PESs represented by n-mode expansions as needed for the accurate calculation of anharmonic vibrational frequencies by means of vibrational configuration interaction theory. A static approach, which has successfully been used in many applications, and five adaptive schemes based on Gaussian process regression have been investigated with respect to the number of necessary grid points and the accuracy of the fundamental modes for a small set of test molecules. A comparison with a related, more sophisticated, and consistent approach by Christiansen et al. is provided. The impact of the positions of the abinitio grid points is discussed for multilevel PESs, for which the computational effort of the individual electronic structure calculations decreases for increasing orders of the n-mode expansion. As a result of that, the ultimate goal is not the maximal reduction of grid points but rather the computational cost, which is not directly related.

Structure, spectroscopic signatures, and formation of hydroxy-azirine: a potential interstellar prebiotic molecule

ABSTRACT Hydroxy-azirine (C2H3NO) is a -OH derivative of azirine (C2H3N), molecule that has been the subject of several unfruitful searches in space. Hydroxy-azirine is an isomer of the detected prebiotic species methyl isocyanate, CH3NCO, and glycolonitrile, HOCH2CN, as well as the yet undetected imine acetaldehyde, NHCHCHO. However, the lack of preliminary spectroscopic data on hydroxy-azirine has prevented its astronomical search. The aim of this study is to provide high-level theoretical spectroscopic signatures of the most stable hydroxy-azirine isomers to enable their eventual interstellar search. A total of 12 isomers have been characterized for hydroxy-azirine and their isomerization processes have been analysed at the CCSD(T)-F12/cc-pVTZ-F12 level. The most stable structures are 3-hydroxy-2H-azirine (I) and 2-hydroxy-2H-azirine (II) in their syn- and anti-configurations, which are suggested as the most relevant candidates for laboratory and interstellar detection. To ease their identification by means of rotational spectroscopy, we report a set of the required spectroscopic parameters using state-of-the-art composite and coupled-cluster approaches. For astronomical purposes, we provide a complete line list for I-syn and I-antihydroxy-azirine up to 50 GHz, which takes the hyperfine structure into account, and will be essential to hunt for these interstellar candidates experimentally. In addition, anharmonic vibrational frequencies and intensities are reported to predict a trustworthy vibrational spectra and to estimate the vibrational partition function. Finally, we analyse the possibility of formation of hydroxy-azirine from the reaction of azirine with the hydroxyl radical in the gas-phase and on the surface of ices, finding for the latter a feasible formation route under interstellar conditions.

Elucidating NO x Surface Chemistry at the Anatase (101) Surface in TiO2 Nanoparticles.

Understanding NO x chemistry at titania nanoparticle surfaces is important for photocatalytic environmental remediation processes. We focus on this problem and put forward an experimental-computational approach based on vibrational spectroscopy grounds. Temperature-dependent IR experiments of NO x adsorption on shape-engineered nanoparticle (101) anatase surfaces are paired with power spectra obtained from Born-Oppenheimer trajectories. Then, the harmonic versus anharmonic vibrational frequencies of several adsorption scenarios are directly compared with the IR experiments. We conclude that molecules are adsorbed mainly by the N-end side and both the intermolecular interactions between adsorbed molecules and (NO)2 dimer formation are responsible for the main NO adsorption spectroscopic features. We also investigate the spectroscopy and the mechanism of formation on defective anatase surfaces of the long-lived greenhouse gas N2O.

The spectral features and detectability of small, cyclic silicon carbide clusters

Rovibrational spectral data for several tetra-atomic silicon carbide clusters (TASCCs) are computed in this work using a CCSD(T)-F12b/cc-pCVTZ-F12 quartic force field. Accurate theoretical spectroscopic data may facilitate the observation of TASCCs in the interstellar medium which may lead to a more complete understanding of how the smallest silicon carbide (SiC) solids are formed. Such processes are essential for understanding SiC dust grain formation. Due to SiC dust prevalence in the interstellar medium, this may also shed light on subsequent planetary formation. Rhomboidal Si2C2 is shown here to have a notably intense (247 km mol−1) anharmonic vibrational frequency at 988.1 cm−1 (10.1 μm) for ν2, falling into one of the spectral emission features typically associated with unknown infrared bands of various astronomical regions. Notable intensities are also present for several of the computed anharmonic vibrational frequencies including the cyclic forms of C4, SiC3, Si3C, and Si4. These features in the 6–10 μm range are natural targets for infrared observation with the James Webb Space Telescope (JWST)’s MIRI instrument. Additionally, t-Si2C2, d-Si3C, and r-SiC3 each possess dipole moments of greater than 2.0 D making them interesting targets for radioastronomical searches especially since d-SiC3 is already known in astrophysical media.

Reparameterized semi-empirical methods for computing anharmonic vibrational frequencies of multiply-bonded hydrocarbons

Abstract Reparameterized semi-empirical methods can reproduce gas-phase experimental vibrational frequencies to within 24 cm−1 or better for a 100-fold decrease in computational cost in the anharmonic fundamental vibrational frequencies. To achieve such accuracy and efficiency, the default parameters in the PM6 semi-empirical model are herein optimized to reproduce the experimental and high-level theoretical vibrational spectra of three small hydrocarbon molecules, C2H2, c-C3H2, and C2H4, with the hope that these same parameters will be applicable to large polycyclic aromatic hydrocarbons (PAHs). This massive cost reduction allows for the computation of explicit anharmonic frequencies and the inclusion of resonance corrections that have been shown to be essential for accurate predictions of anharmonic frequencies. Such accurate predictions are necessary to help to disentangle the heretofore unidentified infrared spectral features observed around diverse astronomical bodies and hypothesized to be caused by PAHs, especially with the upcoming influx of observational data from the James Webb Space Telescope. The optimized PM6 parameters presented herein represent a substantial step in this direction with those obtained for ethylene (C2H4) yielding a 37% reduction in the mean absolute error of the fundamental frequencies compared to the default PM6 parameters.