Degenerate perturbation theory corrections for the vibrational self-consistent field approximation: Method and applications

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.

The first HyDRA challenge for computational vibrational spectroscopy

Vibrational spectroscopy in supersonic jet expansions is a powerful tool to assess molecular aggregates in close to ideal conditions for the benchmarking of quantum chemical approaches. The low temperatures achieved as well as the absence of environment effects allow for a direct comparison between computed and experimental spectra. This provides potential benchmarking data which can be revisited to hone different computational techniques, and it allows for the critical analysis of procedures under the setting of a blind challenge. In the latter case, the final result is unknown to modellers, providing an unbiased testing opportunity for quantum chemical models. In this work, we present the spectroscopic and computational results for the first HyDRA blind challenge. The latter deals with the prediction of water donor stretching vibrations in monohydrates of organic molecules. This edition features a test set of 10 systems. Experimental water donor OH vibrational wavenumbers for the vacuum-isolated monohydrates of formaldehyde, tetrahydrofuran, pyridine, tetrahydrothiophene, trifluoroethanol, methyl lactate, dimethylimidazolidinone, cyclooctanone, trifluoroacetophenone and 1-phenylcyclohexane-cis-1,2-diol are provided. The results of the challenge show promising predictive properties in both purely quantum mechanical approaches as well as regression and other machine learning strategies.

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.

Amino Acids Compete with Ammonia in Sulfuric Acid-Based Atmospheric Aerosol Prenucleation: The Case of Glycine and Serine

We present a computational investigation of the sulfuric acid, glycine, serine, ammonia, and water system to understand if this system can form prenucleation clusters, which are precursors to larger aerosols in the atmosphere. We have performed a comprehensive configurational search of all possible clusters in this system, starting with the four different monomers and zero to five waters. Accurate Gibbs free energies of formation have been calculated with the DLPNO-CCSD(T)/complete basis set (CBS) method on ωb97xd/6-31++G** geometries. For the dry dimers of sulfuric acid, the weakest base, serine, is found to form the most stable complex, which is a consequence of the strong di-ionic complex formed between the bisulfate ion and the protonated serine cation. For the dry dimers without sulfuric acid, the glycine-serine complex is more stable than the glycine-ammonia or serine-ammonia complexes, stemming from the detailed structure and not related to base strength. For the larger complexes, sulfuric acid deprotonates and the proton is shifted to glycine, serine, or ammonia. The two amino acids and ammonia are almost interchangeable and there is no easy way to predict which molecule will be protonated without the calculated results. Assuming reasonable starting concentrations and a closed system of sulfuric acid, glycine, serine, ammonia, and five waters, we predict the concentrations of all possible complexes at two temperatures spanning the troposphere. The most negative ΔG° values are a function of the detailed molecular interactions of these clusters. These details are more important than the base strength of ammonia, glycine, and serine.

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).

Investigation of the Failure of the MP2 Method to Describe the Out-of-Plane Bending Motions of Carbon-Carbon Double-Bonded Molecules: The Role of Atomic Orbitals

We present a detailed study on the role of atomic orbitals in the unphysical description of the out-of-plane bending (OOPB) vibrations with the MP2 method. The anharmonicities of the OOPB vibrations are found to be unphysically large for several basis sets with the MP2 method. We find that the inclusion of additional valence s and p orbitals to the basis set leads to the lowering of π* orbitals' energies along with the generation of several spurious low-energy virtual molecular orbitals (VMOs). Such erroneous VMOs create a disbalance between the σ and π correlations for the planar structure of the molecules. Since the OOPB vibrations distort the π* orbitals, the disbalance in the σ and π correlations leads to errors in the derivatives of the correlation energy with respect to the displacement along with OOPB modes. The diffuse functions in the basis set enhance the unphysical anharmonicities of the OOPB modes since these functions are of s- and p-type atomic orbitals. The polarization atomic orbitals of higher angular momentum improve the description of the π* orbitals and thereby reduce the disbalance between the σ and π correlations in the MP2 calculations. We find that the unphysical frequencies for the OOPB modes are significantly eliminated with the d and f orbitals on the C atoms in the basis set.

Reduced Two-Electron Interactions in Anharmonic Molecular Vibrational Calculations Involving Localized Normal Coordinates

Spatially localized vibrational normal mode coordinates are shown to reduce the importance of calculating the full set of two-electron terms in the molecular electronic Schrödinger equation. Electron correlation and dispersion interactions become less significant in (E,E)-1,3,5,7-octatetraene vibrational self-consistent field calculations when displacing remote atoms along multiple coordinates. Electron correlation interactions between spatially remote modes are also found to be less important compared to their corresponding uncorrelated interaction terms. Attenuation of the Coulomb operator indicates that the two-electron terms between remote electrons become less important for accurately describing the strongly contributing mode-coupling terms between sets of localized vibrational modes.

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.

Probing vibrational coupling via a grid-based quantum approach-an efficient strategy for accurate calculations of localized normal modes in solid-state systems

In this work an approach to investigate the properties of strongly localized vibrational modes of functional groups in bulk material and on solid-state surfaces is presented. The associated normal mode vectors are approximated solely on the basis of structural information and obtained via diagonalization of a reduced Hessian. The grid-based Numerov procedure in one and two dimensions is then applied to an adequate scan of the respective potential surface yielding the associated vibrational wave functions and energy eigenvalues. This not only provides a detailed description of anharmonic effects but also an accurate inclusion of the coupling between the investigated vibrational states on a quantum mechanical level. All results obtained for the constructed normal modes are benchmarked against their analytical counterparts obtained from the diagonalization of the total Hessian of the entire system. Three increasingly complex systems treated at quantum chemical level of theory have been considered, namely the symmetric and asymmetric stretch vibrations of an isolated water molecule, hydroxyl groups bound to the surface of GeO2 (001), α-quartz(001) and Rutil (001) as well as crystalline Li2 NH serving as an example for a bulk material. While the data obtained for the individual systems verify the applicability of the proposed methodology, comparison to experimental data demonstrates the accuracy of this methodology despite the restriction to limit this methodology to a few selected vibrational modes. The possibility to investigate vibrational phenomena of localized normal modes without the requirement of executing costly harmonic frequency calculations of the entire system enables the application of this method to cases in which the determination of normal modes is prohibitively expensive or not available for a particular level of theory. © 2018 Wiley Periodicals, Inc.

Uncertainty quantification for quantum chemical models of complex reaction networks

For the quantitative understanding of complex chemical reaction mechanisms, it is, in general, necessary to accurately determine the corresponding free energy surface and to solve the resulting continuous-time reaction rate equations for a continuous state space. For a general (complex) reaction network, it is computationally hard to fulfill these two requirements. However, it is possible to approximately address these challenges in a physically consistent way. On the one hand, it may be sufficient to consider approximate free energies if a reliable uncertainty measure can be provided. On the other hand, a highly resolved time evolution may not be necessary to still determine quantitative fluxes in a reaction network if one is interested in specific time scales. In this paper, we present discrete-time kinetic simulations in discrete state space taking free energy uncertainties into account. The method builds upon thermo-chemical data obtained from electronic structure calculations in a condensed-phase model. Our kinetic approach supports the analysis of general reaction networks spanning multiple time scales, which is here demonstrated for the example of the formose reaction. An important application of our approach is the detection of regions in a reaction network which require further investigation, given the uncertainties introduced by both approximate electronic structure methods and kinetic models. Such cases can then be studied in greater detail with more sophisticated first-principles calculations and kinetic simulations.

Reduced Basis Set Dependence in Anharmonic Frequency Calculations Involving Localized Coordinates

Localized normal coordinates are known to be effective in speeding up anharmonic frequency calculations by reducing the complexity of the nuclear Hamiltonian and wave function. Displacing atoms in localized coordinates can also cause relatively small changes in the electronic structure, which can be exploited for further computational efficiency improvements during ab initio electronic structure calculations of the potential energy surface by reducing the electronic basis set dependence. Three different schemes for reducing the basis set dependence have been investigated in this work. These include combining localized coordinate schemes with general mixed basis sets, distance based force-field reductions, and using coordinate specific basis sets. The importance of accurately describing electronic interactions is found to diminish both for multicoordinate terms involving the displacement of remote atoms and when describing the interactions between more remote atoms within specific coordinates.

Intermediate vibrational coordinate localization with harmonic coupling constraints

Optimized normal coordinates can significantly improve the speed and accuracy of vibrational frequency calculations. However, over-localization can occur when using unconstrained spatial localization techniques. The unintuitive mixtures of stretching and bending coordinates that result can make interpreting spectra more difficult and also cause artificial increases in mode-coupling during anharmonic calculations. Combining spatial localization with a constraint on the coupling between modes can be used to generate coordinates with properties in-between the normal and fully localized schemes. These modes preserve the diagonal nature of the mass-weighted Hessian matrix to within a specified tolerance and are found to prevent contamination between the stretching and bending vibrations of the molecules studied without a priori classification of the different types of vibration present. Relaxing the constraint can also be used to identify which normal modes form specific groups of localized modes. The new coordinates are found to center on more spatially delocalized functional groups than their fully localized counterparts and can be used to tune the degree of vibrational correlation energy during anharmonic calculations.

Proton relays in anomalous carbocations dictate spectroscopy, stability, and mechanisms: case studies on C2H5+ and C3H3+

We present a detailed exposition of “Grothuss-like” proton shuttles in C₂H₅⁺ and C₃H₃⁺ that result in anomalous structural and spectral behavior.

Exploring the relationship between vibrational mode locality and coupling using constrained optimization

The use of alternate coordinate systems as a means to improve the efficiency and accuracy of anharmonic vibrational structure analysis has seen renewed interest in recent years. While normal modes (which diagonalize the mass-weighted Hessian matrix) are a typical choice, the delocalized nature of this basis makes it less optimal when anharmonicity is in play. When a set of modes is not designed to treat anharmonicity, anharmonic effects will contribute to inter-mode coupling in an uncontrolled fashion. These effects can be mitigated by introducing locality, but this comes at its own cost of potentially large second-order coupling terms. Herein, a method is described which partially localizes vibrations to connect the fully delocalized and fully localized limits. This allows a balance between the treatment of harmonic and anharmonic coupling, which minimizes the error that arises from neglected coupling terms. Partially localized modes are investigated for a range of model systems including a tetramer of hydrogen fluoride, water dimer, ethene, diphenylethane, and stilbene. Generally, partial localization reaches ∼75% of maximal locality while introducing less than ∼30% of the harmonic coupling of the fully localized system. Furthermore, partial localization produces mode pairs that are spatially separated and thus weakly coupled to one another. It is likely that this property can be exploited in the creation of model Hamiltonians that omit the coupling parameters of the distant (and therefore uncoupled) pairs.

Examining the impact of harmonic correlation on vibrational frequencies calculated in localized coordinates

Carefully choosing a set of optimized coordinates for performing vibrational frequency calculations can significantly reduce the anharmonic correlation energy from the self-consistent field treatment of molecular vibrations. However, moving away from normal coordinates also introduces an additional source of correlation energy arising from mode-coupling at the harmonic level. The impact of this new component of the vibrational energy is examined for a range of molecules, and a method is proposed for correcting the resulting self-consistent field frequencies by adding the full coupling energy from connected pairs of harmonic and pseudoharmonic modes, termed vibrational self-consistent field (harmonic correlation). This approach is found to lift the vibrational degeneracies arising from coordinate optimization and provides better agreement with experimental and benchmark frequencies than uncorrected vibrational self-consistent field theory without relying on traditional correlated methods.

Optimizing Vibrational Coordinates To Modulate Intermode Coupling

The choice of coordinate system strongly affects the convergence properties of vibrational structure computations. Two methods for efficient generation of improved vibrational coordinates are presented and justified by analysis of a model anharmonic two-mode Hessian and numerical computations on polyatomic molecules. To produce optimal coordinates, metrics which quantify off-diagonal couplings over a grid of Hessian matrices are minimized through unitary rotations of the vibrational basis. The first proposed metric minimizes the total squared off-diagonal coupling, and the second minimizes the total squared change in off-diagonal coupling. In this procedure certain anharmonic modes tend to localize, for example X-H stretches. The proposed methods do not rely on prior fitting of the potential energy, vibrational structure computations, or localization metrics, so they are unique from previous vibrational coordinate generation algorithms and are generally applicable to polyatomic molecules. Fitting the potential to the approximate n-mode representation in the optimized bases for all-trans polyenes shows that off-diagonal anharmonic couplings are substantially reduced by the new choices of coordinate system. Convergence of vibrational energies is examined in detail for ethylene, and it is shown that coupling-optimized modes converge in vibrational configuration interaction computations to within 1 cm(-1) using only 3-mode couplings, where normal modes require 4-mode couplings for convergence. Comparison of the vibrational configuration interaction convergence with respect to excitation level for the two proposed metrics shows that minimization of the total off-diagonal coupling is most effective for low-cost vibrational structure computations.

Scaling Factors and Uncertainties for ab Initio Anharmonic Vibrational Frequencies

To predict the vibrational spectra of molecules, ab initio calculations are often used to compute harmonic frequencies, which are usually scaled by empirical factors as an approximate correction for errors in the force constants and for anharmonic effects. Anharmonic computations of fundamental frequencies are becoming increasingly popular. We report scaling factors, along with their associated uncertainties, for anharmonic (second-order perturbation theory) predictions from HF, MP2, and B3LYP calculations using the 6-31G(d) and 6-31+G(d,p) basis sets. Different scaling factors are appropriate for low- and high-frequency vibrations. The method of analysis is based upon the Guide to the Expression of Uncertainty in Measurement, published by the International Organization for Standardization (ISO). The data used are from the Computational Chemistry Comparison and Benchmark Database (CCCBDB), maintained by the National Institute of Standards and Technology, which includes more than 3939 independent vibrations for 358 molecules.

Efficient anharmonic vibrational spectroscopy for large molecules using local-mode coordinates

This article presents a general computational approach for efficient simulations of anharmonic vibrational spectra in chemical systems. An automated local-mode vibrational approach is presented, which borrows techniques from localized molecular orbitals in electronic structure theory. This approach generates spatially localized vibrational modes, in contrast to the delocalization exhibited by canonical normal modes. The method is rigorously tested across a series of chemical systems, ranging from small molecules to large water clusters and a protonated dipeptide. It is interfaced with exact, grid-based approaches, as well as vibrational self-consistent field methods. Most significantly, this new set of reference coordinates exhibits a well-behaved spatial decay of mode couplings, which allows for a systematic, a priori truncation of mode couplings and increased computational efficiency. Convergence can typically be reached by including modes within only about 4 Å. The local nature of this truncation suggests particular promise for the ab initio simulation of anharmonic vibrational motion in large systems, where connection to experimental spectra is currently most challenging.

Largely reduced grid densities in a vibrational self-consistent field treatment do not significantly impact the resultingwavenumbers

Especially for larger molecules relevant to life sciences, vibrational self-consistent field (VSCF) calculations can become unmanageably demanding even when only first and second order potential coupling terms are considered. This paper investigates to what extent the grid density of the VSCF’s underlying potential energy surface can be reduced without sacrificing accuracy of the resulting wavenumbers. Including single-mode and pair contributions, a reduction to eight points per mode did not introduce a significant deviation but improved the computational efficiency by a factor of four. A mean unsigned deviation of 1.3% from the experiment could be maintained for the fifteen molecules under investigation and the approach was found to be applicable to rigid, semi-rigid and soft vibrational problems likewise. Deprotonated phosphoserine, stabilized by two intramolecular hydrogen bonds, was investigated as an exemplary application.

The impact of highly correlated potential energy surfaces on the anharmonically corrected IR spectrum of acetonitrile

This paper discusses the quality and feasibility of highly correlated ab initio techniques in a vibrational self-consistent field (VSCF) approach using acetonitrile as a model system. The topical renormalized coupled-cluster technique exploiting the similarity-transformed Hamiltonian's left eigenstates (i.e. CR-CC(2,3)) is investigated alongside the well-known Hartree-Fock (HF), Møller-Plesset second-order perturbation theory (MP2) and coupled cluster (CCSD(T)) methods. The inclusion of mode triple interactions is discussed and it is found that the use of an effective core potential (ECP) serves as a viable compromise during the highly demanding task of computing such contributions, thus enabling a grid-based evaluation of three mode interaction terms with coupled cluster techniques also for larger molecules. In this context, a previously proposed reduced coupling scheme [1] is investigated, confirming the applicability of this technique to a system exhibiting a rather complex electronic structure. A combination of Ahlrichs' triple-ζ valence polarized (TZVP) basis set with Dunning's set of core-valence correlation functions is found to deliver results in good agreement with experiment while being computationally very feasible. Since CH3CN exhibits four degenerate vibrational degrees of freedom, it serves as an ideal model system for critically assessing the qualities of the degenerate second-order perturbation theory corrected (DPT2) VSCF technique. Besides fundamental vibrations, a thorough investigation of overtone transitions and combination bands is conducted by means of comparing the results to both available and newly recorded experimental data.

Approximate first-principles anharmonic calculations of polyatomic spectra using MP2 and B3LYP potentials: comparisons with experiment

Anharmonic vibrational spectroscopy calculations using MP2 and B3LYP computed potential surfaces are carried out for a series of molecules, and frequencies and intensities are compared with those from experiment. The vibrational self-consistent field with second-order perturbation correction (VSCF-PT2) is used in computing the spectra. The test calculations have been performed for the molecules HNO3, C2H4, C2H4O, H2SO4, CH3COOH, glycine, and alanine. Both MP2 and B3LYP give results in good accord with experimental frequencies, though, on the whole, MP2 gives very slightly better agreement. A statistical analysis of deviations in frequencies from experiment is carried out that gives interesting insights. The most probable percentage deviation from experimental frequencies is about -2% (to the red of the experiment) for B3LYP and +2% (to the blue of the experiment) for MP2. There is a higher probability for relatively large percentage deviations when B3LYP is used. The calculated intensities are also found to be in good accord with experiment, but the percentage deviations are much larger than those for frequencies. The results show that both MP2 and B3LYP potentials, used in VSCF-PT2 calculations, account well for anharmonic effects in the spectroscopy of molecules of the types considered.