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.

Potential energy surfaces and dynamic properties via ab initio composite and density functional approaches

Vibrational spectroscopy enables critical insight into the structural and dynamic properties of molecules. Presently, the majority of theoretical approaches to spectroscopy employ wavefunction-based ab initio or density functional methods that rely on the harmonic approximation. This approximation breaks down for large molecules with strongly anharmonic bonds or for molecules with large internuclear separations. An alternative to these methods involves generating molecular anharmonic potential energy surfaces (potentials) and using them to extrapolate the vibrational frequencies. This study examines the efficacy of density functional theory (DFT) and the correlation consistent Composite Approach (ccCA) in generating anharmonic frequencies from potentials of small main group molecules. Vibrational self-consistent field Theory (VSCF) and post-VSCF methods were used to calculate the fundamental frequencies of these molecules from their potentials. Functional choice, basis set selection, and mode-coupling are also examined as factors in influencing accuracy. The absolute deviations for the calculated frequencies using potentials at the ccCA level of theory were lower than the potentials at the DFT level. With DFT resulting in bending modes that are better described than those of ccCA, a multilevel DFT:ccCA approach where DFT potentials are used for single vibrational mode potentials and ccCA is used for vibrational mode-mode couplings can be utilized for larger polyatomic systems. The frequencies obtained with this multilevel approach using VCIPSI-PT2 were closer to experimental frequencies than the scaled harmonic frequencies, indicating the success of utilizing post-VSCF methods to generate more accurate representations of computed infrared spectra.

The Near-Sightedness of Many-Body Interactions in Anharmonic Vibrational Couplings

Couplings between vibrational motions are driven by electronic interactions, and these couplings carry special significance in vibrational energy transfer, multidimensional spectroscopy experiments, and simulations of vibrational spectra. In this investigation, the many-body contributions to these couplings are analyzed computationally in the context of clathrate-like alkali metal cation hydrates, including Cs+(H2O)20, Rb+(H2O)20, and K+(H2O)20, using both analytic and quantum-chemistry potential energy surfaces. Although the harmonic spectra and one-dimensional anharmonic spectra depend strongly on these many-body interactions, the mode-pair couplings were, perhaps surprisingly, found to be dominated by one-body effects, even in cases of couplings to low-frequency modes that involved the motion of multiple water molecules. The origin of this effect was traced mainly to geometric distortion within water monomers and cancellation of many-body effects in differential couplings, and the effect was also shown to be agnostic to the identity of the ion. These outcomes provide new understanding of vibrational couplings and suggest the possibility of improved computational methods for the simulation of infrared and Raman spectra.

Pitfalls in the n-mode representation of vibrational potentials

Simulations of anharmonic vibrational motion rely on computationally expedient representations of the governing potential energy surface. The n-mode representation (n-MR)-effectively a many-body expansion in the space of molecular vibrations-is a general and efficient approach that is often used for this purpose in vibrational self-consistent field (VSCF) calculations and correlated analogues thereof. In the present analysis, a lack of convergence in many VSCF calculations is shown to originate from negative and unbound potentials at truncated orders of the n-MR expansion. For cases of strong anharmonic coupling between modes, the n-MR can both dip below the true global minimum of the potential surface and lead to effective single-mode potentials in VSCF that do not correspond to bound vibrational problems, even for bound total potentials. The present analysis serves mainly as a pathology report of this issue. Furthermore, this insight into the origin of VSCF non-convergence provides a simple, albeit ad hoc, route to correct the problem by "painting in" the full representation of groups of modes that exhibit these negative potentials at little additional computational cost. Somewhat surprisingly, this approach also reasonably approximates the results of the next-higher n-MR order and identifies groups of modes with particularly strong coupling. The method is shown to identify and correct problematic triples of modes-and restore SCF convergence-in two-mode representations of challenging test systems, including the water dimer and trimer, as well as protonated tropine.

Accuracy of Different Electronic Basis Set Families for Anharmonic Molecular Vibrations: A Comprehensive Benchmark Study

In this work, the accuracy and convergence of different electronic basis set families for the computation of anharmonic molecular vibrational spectroscopic calculations are benchmarked. A series of 39 different basis sets from different families following their hierarchy are assessed on VSCF and VSCF-PT2 algorithms with commonly used MP2 and DFT based B3LYP-D potentials for a set of molecular systems. Such an effort has been validated in a previous work ( J. Phys. Chem. A 2020, 124, 9203-9221) with split-valence basis sets for fundamentals and intensities. Here, fundamental transitions, vibrationally excited states, and intensities are compared with the experimental data to estimate the accuracy for a series of Jensen, Dunning, Calendar, Karlsruhe, and Sapporo basis set families. The convergence of basis sets are also compared with the large ANO basis set. Comprehensive statistical error analysis in terms of accuracy and precision was carried out to assess the performance of each basis set. It is observed that the improvement for the calculated harmonic and anharmonic values from the smaller basis sets to the medium (i.e., triple-ξ) is considerable. Beyond this, from medium to large basis sets, the convergence is slow and mostly posits nearly converged values. Basis sets with and without diffuse functions offer characteristically different accuracies and convergence patterns. Finally, recommendations are given on the choice of basis set chosen as black-box which can balance between accuracy and computational time, estimation of the errors, and their selections especially for large molecules.

Accelerating Anharmonic Spectroscopy Simulations via Local-Mode, Multilevel Methods

Ab initio computer simulations of anharmonic vibrational spectra provide nuanced insight into the vibrational behavior of molecules and complexes. The computational bottleneck in such simulations, particularly for ab initio potentials, is often the generation of mode-coupling potentials. Focusing specifically on two-mode couplings in this analysis, the combination of a local-mode representation and multilevel methods is demonstrated to be particularly symbiotic. In this approach, a low-level quantum chemistry method is employed to predict the pairwise couplings that should be included at the target level of theory in vibrational self-consistent field (and similar) calculations. Pairs that are excluded by this approach are "recycled" at the low level of theory. Furthermore, because this low-level pre-screening will eventually become the computational bottleneck for sufficiently large chemical systems, distance-based truncation is applied to these low-level predictions without substantive loss of accuracy. This combination is demonstrated to yield sub-wavenumber fidelity with reference vibrational transitions when including only a small fraction of target-level couplings; the overhead of predicting these couplings, particularly when employing distance-based, local-mode cutoffs, is a trivial added cost. This combined approach is assessed on a series of test cases, including ethylene, hexatriene, and the alanine dipeptide. Vibrational self-consistent field (VSCF) spectra were obtained with an RI-MP2/cc-pVTZ potential for the dipeptide, at approximately a 5-fold reduction in computational cost. Considerable optimism for increased accelerations for larger systems and higher-order couplings is also justified, based on this investigation.

Performance of Vibrational Self-Consistent Field Theory for Accurate Potential Energy Surfaces: Fundamentals, Excited States, and Intensities

The performance of vibrational structure calculations beyond harmonic approximation in the framework of the vibrational self-consistent field method with second-order perturbation corrections (VSCF-PT2) is investigated in conjunction with very accurate potential energy surfaces (PESs) given by various coupled-cluster electronic structure theories. The quality of anharmonic calculations depends on the accuracy of the underlying multidimensional PES obtained from its functional form, which is given by the level of electronic structure theory. Two such highest levels of typical coupled-cluster electronic structure methods, CCSD and the ″gold standard″ CCSD(T), along with their variants such as CCD, CR-CCL (completely renormalized CR-CC(2,3) approach), and CCSD(TQ) are tested for the construction of accurate anharmonic potentials without any fitting or ad hoc scaling and using cc-pVTZ basis sets. The accuracy of VSCF-PT2 theory in comparison to experimental values is tested for a series of 16 molecules with 135 fundamental bands, 64 overtones, and combination bands and also for 39 intensities. It is found that CCD and CCSD bind the potential tighter than CCSD(T) and the computed VSCF-PT2 transitions are more blue-shifted showing higher deviation from the experiment. In general, VSCF-PT2 results computed at the CCSD(T) potential offer a good cost/accuracy ratio, with the mean absolute deviation and the mean absolute percentage error with the experiment being ∼16 cm^-1 and 1.38, respectively, for fundamentals. Additionally, while the CR-CCL and CCSD(TQ) methods offer similar levels of accuracies as compared to CCSD(T), the former offers a better accuracy/cost ratio than the latter and is a suitable alternative to CCSD(T).

Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case

The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)^MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm^-1, has been obtained by combining CCSD(T)/(CBS+CV)^MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.

Advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase

The purpose of this review is to demonstrate advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase. Molecular spectroscopy, particularly vibrational spectroscopy and electronic spectroscopy, has been used extensively for a wide range of areas of chemical sciences and materials science as well as nano- and biosciences because it provides valuable information about structure, functions, and reactions of molecules. In the meantime, quantum chemical approaches play crucial roles in the spectral analysis. They also yield important knowledge about molecular and electronic structures as well as electronic transitions. The combination of spectroscopic approaches and quantum chemical calculations is a powerful tool for science, in general. Thus, our article, which treats various spectroscopy and quantum chemical approaches, should have strong implications in the wider scientific community. This review covers a wide area of molecular spectroscopy from far-ultraviolet (FUV, 120-200 nm) to far-infrared (FIR, 400-10 cm-1)/terahertz and Raman spectroscopy. As quantum chemical approaches, we introduce several anharmonic approaches such as vibrational self-consistent field (VSCF) and the combination of periodic harmonic calculations with anharmonic corrections based on finite models, grid-based techniques like the Numerov approach, the Cartesian coordinate tensor transfer (CCT) method, Symmetry-Adapted Cluster Configuration-Interaction (SAC-CI), and the ZINDO (Semi-empirical calculations at Zerner's Intermediate Neglect of Differential Overlap). One can use anharmonic approaches and grid-based approaches for both infrared (IR) and near-infrared (NIR) spectroscopy, while CCT methods are employed for Raman, Raman optical activity (ROA), FIR/terahertz and low-frequency Raman spectroscopy. Therefore, this review overviews cross relations between molecular spectroscopy and quantum chemical approaches, and provides various kinds of close-reality advanced spectral simulation for condensed phases.

Theoretical Simulation of Near-Infrared Spectrum of Piperine: Insight into Band Origins and the Features of Regression Models

We investigated the near-infrared spectrum of piperine using quantum mechanical calculations. We evaluated two efficient approaches, DVPT2//PM6 and DVPT2//ONIOM [PM6:B3LYP/6-311++G(2df, 2pd)] that yielded a simulated spectrum with varying accuracy versus computing time factor. We performed vibrational assignments and unveiled complex nature of the near-infrared spectrum of piperine, resulting from a high level of band convolution. The most meaningful contribution to the near-infrared absorption of piperine results from binary combination bands. With the available detailed near-infrared assignment of piperine, we interpreted the properties of partial least square regression models constructed in our earlier study to describe the piperine content in black pepper samples. Two models were compared with spectral data sets obtained with a benchtop and a miniaturized spectrometer. The two spectrometers implement distinct technology which leads to a profound instrumental difference and discrepancy in the predictive performance when analyzing piperine content. We concluded that the sensitivity of the two instruments to certain types of piperine vibrations is different and that the benchtop spectrometer unveiled higher selectivity. Such difference in obtaining chemical information from a sample can be one of the reasons why the benchtop spectrometer performs better in analyzing the piperine content of black pepper. This evidenced direct correspondence between the features critical for applied near-infrared spectroscopic routine and the underlying vibrational properties of the analyzed constituent in a complex sample.

My Trajectory in Molecular Reaction Dynamics and Spectroscopy

This is the story of a career in theoretical chemistry during a time of dramatic changes in the field due to phenomenal growth in the availability of computational power. It is likewise the story of the highly gifted graduate students and postdoctoral fellows that I was fortunate to mentor throughout my career. It includes reminiscences of the great mentors that I had and of the exciting collaborations with both experimentalists and theorists on which I built much of my research. This is an account of the developments of exciting scientific disciplines in which I was involved: vibrational spectroscopy, molecular reaction mechanisms and dynamics, e.g., in atmospheric chemistry, and the prediction of new, exotic molecules, in particular noble gas molecules. From my very first project to my current work, my career in science has brought me the excitement and fascination of research. What a wonderful pursuit!

Dual Basis Approach for Ab Initio Anharmonic Calculations of Vibrational Spectroscopy: Application to Microsolvated Biomolecules

A dual electronic basis set approach is introduced for more efficient but accurate calculations of the anharmonic vibrational spectra in the framework of the vibrational self-consistent field (VSCF) theory. In this approach, an accurate basis set is used to compute the vibrational spectra at the harmonic level. The results are used to scale the potential surface from a more modest but much more efficient basis set. The scaling is such that at the harmonic level the new, scaled potential agrees with one of the accurate basis sets. The approach is tested in the application of the microsolvated, protected amino acid Ac-Phe-OMe, using the scaled anharmonic hybrid potential in the VSCF and VSCF-PT2 algorithms. The hybrid potential method yields results that are in good accord with the experiment and very close to those obtained in calculations with the high-level, very costly potential from the large basis set. At the same time, the hybrid potential calculations are considerably less expensive. The results of the hybrid calculations are much more accurate than those computed from the potential surface corresponding to the modest basis set. The results are very encouraging for using the hybrid potential method for inexpensive yet sufficiently accurate anharmonic calculations for the spectra of large biomolecules.

Gas phase dynamics, conformational transitions and spectroscopy of charged saccharides: the oxocarbenium ion, protonated anhydrogalactose and protonated methyl galactopyranoside

Protonated intermediates are postulated to be involved in the rate determining step of many sugar reactions. This paper presents a study of protonated sugar species, isolated in the gas phase, using a combination of infrared multiple photon dissociation (IRMPD) spectroscopy, classical ab initio molecular dynamics (AIMD) and quantum mechanical vibrational self-consistent field (VSCF) calculations. It provides a likely identification of the reactive intermediate oxocarbenium ion structure in a d-galactosyl system as well as the saccharide pyrolysis product anhydrogalactose (that suggests oxocarbenium ion stabilization), along with the spectrum of the protonated parent species: methyl d-galactopyranoside-H+. Its vibrational fingerprint indicates intramolecular proton sharing. Classical AIMD simulations for galactosyl oxocarbenium ions, conducted in the temperature range ∼300-350 K (using B3LYP potentials on-the-fly) reveal efficient transitions on the picosecond timescale. Multiple conformers are likely to exist under the experimental conditions and along with static VSCF calculations, they have facilitated the identification of the individual structural motifs of the galactosyl oxocarbenium ion and protonated anhydrogalactose ion conformers that contribute to the observed experimental spectra. These results demonstrate the power of experimental IRMPD spectroscopy combined with dynamics simulations and with computational spectroscopy at the anharmonic level to unravel conformer structures of protonated saccharides, and to provide information on their lifetimes.

Efficiently Calculating Anharmonic Frequencies of Molecular Vibration by Molecular Dynamics Trajectory Analysis

Two efficient methods, the Eckart frame algorithm and the multiorder derivative algorithm, for vibrational frequency calculation directly based on the raw data of atomic trajectory from the state-of-the-art first-principles molecular dynamics simulation are presented. The Eckart frame approach is robust to retrieve the full set of anharmonic fundamental frequencies of any molecule from the atomic trajectory for a sufficiently long molecular dynamics simulation at a temperature close to 0 K. In addition to the fundamental vibrational frequencies, the multiorder derivative approach is universal for the calculations of vibrational frequencies based on the molecular dynamics result in a wide range of temperatures. The accuracy, efficiency, and applicability of these two methods are demonstrated through several successful examples in calculating the anharmonic fundamental vibrational frequencies of methane, ethylene, water, and cyclobutadiene.

Intrinsic structure of pentapeptide Leu-enkephalin: geometry optimization and validation by comparison of VSCF-PT2 calculations with cold ion spectroscopy

The intrinsic structure of an opioid peptide [Ala2, Leu5]-leucine enkephalin (ALE) has been investigated using first-principles based vibrational self-consistent field (VSCF) theory and cold ion spectroscopy. IR-UV double resonance spectroscopy revealed the presence of only one highly abundant conformer of the singly protonated ALE, isolated and cryogenically cooled in the gas phase. High-level quantum mechanical calculations of electronic structures in conjunction with a systematic conformational search allowed for finding a few low-energy candidate structures. In order to identify the observed structure, we computed vibrational spectra of the candidate structures and employed the theory at the semi-empirically scaled harmonic level and at the first-principles based anharmonic VSCF levels. The best match between the calculated "anharmonic" and the measured spectra appeared, indeed, for the most stable candidate. An average of two spectra calculated with different quantum mechanical potentials is proposed for the best match with experiment. The match thus validates the calculated intrinsic structure of ALE and demonstrates the predictive power of first-principles theory for solving structures of such large molecules.

Anharmonic vibrational spectra from double incremental potential energy and dipole surfaces

Using incremental approaches, size limitations for property surface generations are pushed significantly, enabling accurate large molecule anharmonic vibrational spectra calculations.

A Decapeptide Hydrated by Two Waters: Conformers Determined by Theory and Validated by Cold Ion Spectroscopy

The intrinsic structures of biomolecules in the gas phase may not reflect their native solution geometries. Microsolvation of the molecules bridges the two environments, enabling a tracking of molecular structural changes upon hydration at the atomistic level. We employ density functional calculations to compute a large pool of structures and vibrational spectra for a gas-phase complex, in which a doubly protonated decapeptide, gramicidin S, is solvated by two water molecules. Though most vibrations of this large complex are treated in a harmonic approximation, the water molecules and the vibrations of the host ion coupled to them are locally described by a quantum mechanical vibrational self-consistent field theory with second-order perturbation correction (VSCF-PT2). Guided and validated by the available cold ion spectroscopy data, the computational analysis identifies structures of the three experimentally observed conformers of the complex. They, mainly, differ by the hydration sites, of which the one at the Orn side chain is the most important for reshaping the peptide toward its native structure. The study demonstrates the ability of a quantum chemistry approach that intelligently combines the semiempirical and ab initio computations to disentangle a complex interplay of intra- and intermolecular hydrogen bonds in large molecular systems.

Anharmonic Theoretical Vibrational Spectroscopy of Polypeptides

Because of the size of polypeptides and proteins, the quantum-chemical prediction of their vibrational spectra presents an exceptionally challenging task. Here, we address one of these challenges, namely, the inclusion of anharmonicities. By performing the expansion of the potential energy surface in localized-mode coordinates instead of the normal-mode coordinates, it becomes possible to calculate anharmonic vibrational spectra of polypeptides efficiently and reliably. We apply this approach to calculate the infrared, Raman, and Raman optical activity spectra of helical alanine polypeptides consisting of up to 20 amino acids. We find that while anharmonicities do not alter the band shapes, simple scaling procedures cannot account for the different shifts found for the individual bands. This closes an important gap in theoretical vibrational spectroscopy by making it possible to quantify the anharmonic contributions and opens the door to a first-principles calculation of multidimensional vibrational spectra.