Automated calculation of fundamental frequencies: Application to AlH3 using the coupled-cluster singles-and-doubles with perturbative triples method

Perturb-Then-Diagonalize Vibrational Engine Exploiting Curvilinear Internal Coordinates

The present paper is devoted to the implementation and validation of a second-order perturbative approach to anharmonic vibrations, followed by variational treatment of strong couplings (GVPT2) based on curvilinear internal coordinates. The main difference with respect to the customary Cartesian-based formulation is that the kinetic energy operator is no longer diagonal, and has to be expanded as well, leading to additional terms which have to be taken into proper account. It is, however, possible to recast all the equations as well-defined generalizations of the corresponding Cartesian-based counterparts, thus achieving a remarkable simplification of the new implementation. Particular attention is paid to the treatment of Fermi resonances with significant number of test cases analyzed fully, validating the new implementation. The results obtained in this work confirm that curvilinear coordinates strongly reduce the strength of inter-mode couplings compared to their Cartesian counterparts. This increases the reliability of low-order perturbative treatments for semi-rigid molecules and paves the way toward the reliable representation of more flexible molecules where small- and large-amplitude motions can be safely decoupled and treated at different levels of theory.

Matrix Product State Formulation of the MCTDH Theory in Local Mode Representations for Anharmonic Potentials

The matrix product state formulation of the multiconfiguration time-dependent Hartree theory, MPS-MCTDH, reported previously [Kurashige, J. Chem. Phys. 2018, 19, 194114] is extended to realistic anharmonic potentials with n-mode representations beyond the linear vibronic coupling model. For realistic vibrational potentials, the local mode representation should give a more compact representation of the potentials, i.e., lowering the dimensionality of the entanglements, than the normal coordinates, and the MPS-MCTDH formulation should work more efficiently and maintain the accuracy with a small bond dimension of the MPS ansatz. In fact, it was confirmed that the use of the local coordinates made the interaction matrices diagonal dominant and the number of terms in the n-body expansion of the potentials was significantly reduced. The method was applied to the IR spectrum of the CH₂O molecule, the zero-point energies, and the vibrational energy redistribution dynamics of polyenes C_2nH_2n+2. The results showed that the efficiency of the MPS-MCTDH method is significantly accelerated by the use of local coordinates even if the long-range interactions are included in the potential.

General Perturb-Then-Diagonalize Model for the Vibrational Frequencies and Intensities of Molecules Belonging to Abelian and Non-Abelian Symmetry Groups

In this paper, we show that the standard second-order vibrational perturbation theory (VPT2) for Abelian groups can be used also for non-Abelian groups without employing specific equations for two- or threefold degenerate vibrations but rather handling in the proper way all the degeneracy issues and deriving the peculiar spectroscopic signatures of non-Abelian groups (e.g., l -doubling) by a posteriori transformations of the eigenfunctions. Comparison with the results of previous conventional implementations shows a perfect agreement for the vibrational energies of linear and symmetric tops, thus paving the route to the transparent extension of the equations already available for asymmetric tops to the energies of spherical tops and the infrared and Raman intensities of molecules belonging to non-Abelian symmetry groups. The whole procedure has been implemented in our general engine for vibro-rotational computations beyond the rigid rotor/harmonic oscillator model and has been validated on a number of test cases.

Multiresolution continuous wavelet transform for studying coupled solute-solvent vibrations via ab initio molecular dynamics

Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.

Finite Field Method for Nonlinear Optical Property Prediction Using Rational Function Approximants

The finite field (FF) method is a quick, easy-to-implement tool for the prediction of nonlinear optical properties. Here, we present and explore a novel variant of the FF method, which uses a rational function to fit a molecule's energy with respect to an electric field. Similarly to previous FF methods, factors crucial for the method's accuracy were tuned. These factors include the number of terms in the function, the distribution of fields used to construct the approximation, and the initial field in the approximation. It was found that the approximant form that best fits the energy has four numerator terms and three denominator terms. To determine a reasonable field distribution, the common ratio of a geometric progression was optimized to √2. Finally, an algorithm for determining a good initial field guess was devised. The optimized FF method was used to compute the polarizability and second hyperpolarizability for a set of 121 molecules and the first hyperpolarizability for a set of 91 molecules. The results from this were compared to a previous polynomial-based FF method. It was found that using a rational function gives higher errors compared to the polynomial model. However, unlike the polynomial model, no subsequent refinement steps were needed to obtain usable results. An overall comparison of the behavior of the two methods also shows that the rational function is less sensitive to the chosen initial field, making it a good choice for new quantum chemistry codes.

Rovibrational and Temperature Effects in Theoretical Studies of NMR Parameters

The demand for high precision calculations of NMR shieldings (or their related values, chemical shifts δ) and spin–spin coupling constants facilitating and supporting detailed interpretations of NMR spectra increases hand in hand with the development of computational techniques and hardware resources. Highly sophisticated calculations including even relativistic effects are nowadays possible for these properties. However, NMR parameters depend not only on molecular structure and environment but also on molecular flexibility and temperature and the apparent success of theoretical predictions for molecular equilibrium geometries creates a demand for zero-point vibrational and temperature corrections. In this chapter we describe briefly the theory behind rovibrational corrections and review then some important contributions to this field.

CC/DFT Route toward Accurate Structures and Spectroscopic Features for Observed and Elusive Conformers of Flexible Molecules: Pyruvic Acid as a Case Study

The structures and relative stabilities as well as the rotational and vibrational spectra of the three low-energy conformers of pyruvic acid (PA) have been characterized using a state-of-the-art quantum-mechanical approach designed for flexible molecules. By making use of the available experimental rotational constants for several isotopologues of the most stable PA conformer, Tc-PA, the semiexperimental equilibrium structure has been derived. The latter provides a reference for the pure theoretical determination of the equilibrium geometries for all conformers, thus confirming for these structures an accuracy of 0.001 Å and 0.1 deg for bond lengths and angles, respectively. Highly accurate relative energies of all conformers (Tc-, Tt-, and Ct-PA) and of the transition states connecting them are provided along with the thermodynamic properties at low and high temperatures, thus leading to conformational enthalpies accurate to 1 kJ mol(-1). Concerning microwave spectroscopy, rotational constants accurate to about 20 MHz are provided for the Tt- and Ct-PA conformers, together with the computed centrifugal-distortion constants and dipole moments required to simulate their rotational spectra. For Ct-PA, vibrational frequencies in the mid-infrared region accurate to 10 cm(-1) are reported along with theoretical estimates for the transitions in the near-infrared range, and the corresponding infrared spectrum including fundamental transitions, overtones, and combination bands has been simulated. In addition to the new data described above, theoretical results for the Tc- and Tt-PA conformers are compared with all available experimental data to further confirm the accuracy of the hybrid coupled-cluster/density functional theory (CC/DFT) protocol applied in the present study. Finally, we discuss in detail the accuracy of computational models fully based on double-hybrid DFT functionals (mainly at the B2PLYP/aug-cc-pVTZ level) that avoid the use of very expensive CC calculations.

Generalized Vibrational Perturbation Theory for Rotovibrational Energies of Linear, Symmetric and Asymmetric Tops: Theory, Approximations, and Automated Approaches to Deal with Medium-to-Large Molecular Systems

Models going beyond the rigid-rotor and the harmonic oscillator levels are mandatory for providing accurate theoretical predictions for several spectroscopic properties. Different strategies have been devised for this purpose. Among them, the treatment by perturbation theory of the molecular Hamiltonian after its expansion in power series of products of vibrational and rotational operators, also referred to as vibrational perturbation theory (VPT), is particularly appealing for its computational efficiency to treat medium-to-large systems. Moreover, generalized (GVPT) strategies combining the use of perturbative and variational formalisms can be adopted to further improve the accuracy of the results, with the first approach used for weakly coupled terms, and the second one to handle tightly coupled ones. In this context, the GVPT formulation for asymmetric, symmetric, and linear tops is revisited and fully generalized to both minima and first-order saddle points of the molecular potential energy surface. The computational strategies and approximations that can be adopted in dealing with GVPT computations are pointed out, with a particular attention devoted to the treatment of symmetry and degeneracies. A number of tests and applications are discussed, to show the possibilities of the developments, as regards both the variety of treatable systems and eligible methods.

New developments of a multifrequency virtual spectrometer: stereo-electronic, dynamical, and environmental effects on chiroptical spectra

Computational spectroscopy has recently evolved from a field reserved for specialists toward a general tool allowing interpretations and analyses of experimental results. However, the current practice of providing tables of transitions for rigid geometries, possibly tuned by phenomenological broadening, is by far too naive. In order to improve this situation in the last few years we have been developing a general, robust, and user-friendly virtual spectrometer (VS) able to complement experimental studies for complex systems in condensed phases. The VS is based on flexible graphical pre- and postprocessing tools interfaced with general number-crunching software. This last tool is rooted in several electronic structure methodologies (DFT, TD-DFT, post-Hartree-Fock), powerful discrete/continuum models for describing environmental effects, and general vibrational and vibronic models. These last topics are the main focus of this work, which sketches our latest developments related to effective inclusion of anharmonic contributions, together with time-independent and/or time-dependent descriptions of vibronic transitions including Franck-Condon, Herzberg-Teller, and Duschinsky effects. Some test cases are described in some detail with the aim of showing the role of different effects in ruling vibrational (VCD) and electronic (ECD, CPL) chiral spectroscopies.

Dispersion corrected DFT approaches for anharmonic vibrational frequency calculations: nucleobases and their dimers

Computational spectroscopy techniques have become in the last few years an effective means to analyze and assign infrared (IR) spectra of molecular systems of increasing dimensions and in different environments. However, transition from compilation of harmonic data to fully anharmonic simulations of spectra is still underway. The most promising results for large systems have been obtained, in our opinion, by perturbative vibrational approaches based on potential energy surfaces computed by hybrid (especially B3LYP) density functionals and medium size (e.g. SNSD) basis sets. In this framework, we are actively developing a comprehensive and robust computational protocol aimed at quantitative reproduction of the spectra of nucleic acid base complexes and their adsorption on solid supports (organic/inorganic). In this contribution we report the essential results of the first step devoted to isolated monomers and dimers. It is well known that in order to model the vibrational spectra of weakly bound molecular complexes dispersion interactions should be taken into proper account. In this work we have chosen two popular and inexpensive approaches to model dispersion interactions, namely the semi-empirical dispersion correction (D3) and pseudopotential based (DCP) methodologies both in conjunction with the B3LYP functional. These have been used for simulating fully anharmonic IR spectra of nucleobases and their dimers through generalized second order vibrational perturbation theory (GVPT2). We have studied, in particular, isolated adenine, hypoxanthine, uracil, thymine and cytosine, the hydrogen-bonded and stacked adenine and uracil dimers, and the stacked adenine-naphthalene heterodimer. Anharmonic frequencies are compared with standard B3LYP results and experimental findings, while the computed interaction energies and structures of complexes are compared to the best available theoretical estimates.

Fully anharmonic IR and Raman spectra of medium-size molecular systems: accuracy and interpretation

Computation of full infrared (IR) and Raman spectra (including absolute intensities and transition energies) for medium- and large-sized molecular systems beyond the harmonic approximation is one of the most interesting challenges of contemporary computational chemistry. Contrary to common beliefs, low-order perturbation theory is able to deliver results of high accuracy (actually often better than those issuing from current direct dynamics approaches) provided that anharmonic resonances are properly managed. This perspective sketches the recent developments in our research group toward the development of a robust and user-friendly virtual spectrometer rooted in second-order vibrational perturbation theory (VPT2) and usable also by non-specialists essentially as a black-box procedure. Several examples are explicitly worked out in order to illustrate the features of our computational tool together with the most important ongoing developments.

The Dalton quantum chemistry program system

Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, Møller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

Optimized Structure and Vibrational Properties by Error Affected Potential Energy Surfaces

The precise theoretical determination of the geometrical parameters of molecules at the minima of their potential energy surface and of the corresponding vibrational properties are of fundamental importance for the interpretation of vibrational spectroscopy experiments. Quantum Monte Carlo techniques are correlated electronic structure methods promising for large molecules, which are intrinsically affected by stochastic errors on both energy and force calculations, making the mentioned calculations more challenging with respect to other more traditional quantum chemistry tools. To circumvent this drawback in the present work, we formulate the general problem of evaluating the molecular equilibrium structures, the harmonic frequencies, and the anharmonic coefficients of an error affected potential energy surface. The proposed approach, based on a multidimensional fitting procedure, is illustrated together with a critical evaluation of systematic and statistical errors. We observe that the use of forces instead of energies in the fitting procedure reduces the statistical uncertainty of the vibrational parameters by 1 order of magnitude. Preliminary results based on variational Monte Carlo calculations on the water molecule demonstrate the possibility to evaluate geometrical parameters and harmonic and anharmonic coefficients at this level of theory with an affordable computational cost and a small stochastic uncertainty (<0.07% for geometries and <0.7% for vibrational properties).

A second-order perturbation theory route to vibrational averages and transition properties of molecules: general formulation and application to infrared and vibrational circular dichroism spectroscopies

A general formulation to compute anharmonic vibrational averages and transition properties at the second-order of perturbation theory is derived from the Rayleigh-Schrödinger development. This approach is intended to be applicable to any property expanded as a Taylor series up to the third order with respect to normal coordinates or their associated momenta. The equations are straightforward to implement and can be easily adapted to various properties, as illustrated for the case of electric and magnetic dipole moments. From those, infrared and vibrational circular dichroism spectra can be readily obtained. This fully automatic procedure has been applied to several chiral molecules of small-to-medium sizes and compared to the standard double harmonic approximation and to experimental data.

Environmental effects in computational spectroscopy: accuracy and interpretation

Spectroscopic techniques are valuable tools for understanding the structure and dynamics of complex systems, such as biomolecules or nanomaterials. Most of the current research is devoted to the development of new experimental techniques for improving the intrinsic resolution of different spectra. However, the subtle interplay of several different effects acting at different length and time scales still makes the interpretation and analysis of such spectra a very difficult task. In this respect, computational spectroscopy is becoming a needful and versatile tool for the assignment and interpretation of experimental spectra. It is in fact possible nowadays to model with relatively high accuracy the physical-chemical properties of complex molecules in different environments, and to link spectroscopic evidence directly to the structural and dynamical properties of optically or magnetically active solvated probes. In this Review, significant steps toward the simulation of entire spectra in condensed phases are presented together with some basic aspects of computational spectroscopy, which highlight how intramolecular and intermolecular degrees of freedom influence several spectroscopic parameters.

Harmonic and Anharmonic Vibrational Frequency Calculations with the Double-Hybrid B2PLYP Method: Analytic Second Derivatives and Benchmark Studies

This work aims to provide reliable benchmark data on the accuracy of harmonic and anharmonic vibrational frequencies computed with the B2PLYP double-hybrid density functional method. The exchange-correlation contributions required for the B2PLYP analytical second derivatives are presented here, which allow for the effective calculation of harmonic frequency as well as cubic and semidiagonal quartic force fields. The latter, in turn, are necessary to compute the anharmonic vibrational frequencies with the perturbative approach (VPT2). The quality of harmonic vibrational frequencies computed in conjunction with basis sets of double- to quadruple-ζ quality has been checked against reference data from the F38 benchmark set. Then, for an additional set of small closed- and open-shell systems, both harmonic frequencies and anharmonic contributions computed at the B2PLYP/N07D and the B2PLYP/aug-cc-pVTZ levels have been compared to their CCSD(T) counterparts. Moreover, for selected medium-size molecules (furan, pyrrole, thiophene, uracil, anisole, phenol, and pyridine), anharmonic frequencies have been compared to well established experimental results. Such benchmark studies have shown that the B2PLYP/N07D model provides good quality harmonic frequencies and describes correctly anharmonic contributions, the latter being of similar accuracy to their B3LYP/N07D counterparts, but obtained at significantly larger computational cost. Additionally, increased accuracy can be obtained by adopting hybrid models where the B2PLYP/N07D anharmonic contributions are combined with harmonic frequencies computed with more accurate quantum mechanical (QM) approaches or by B2PLYP with larger basis sets. This work confirmed also that most of the recently developed density functionals are significantly less suited for vibrational computations, while the B2PLYP method can be recommended for spectroscopic studies where a good accuracy of vibrational properties is required.

Efficient configuration selection scheme for vibrational second-order perturbation theory

A fast algorithm of vibrational second-order Moller-Plesset perturbation theory is proposed, enabling a substantial reduction in the number of vibrational self-consistent-field (VSCF) configurations that need to be summed in the calculations. Important configurations are identified a priori by assuming that a reference VSCF wave function is approximated well by harmonic oscillator wave functions and that fifth- and higher-order anharmonicities are negligible. The proposed scheme has reduced the number of VSCF configurations by more than 100 times for formaldehyde, ethylene, and furazan with an error in computed frequencies being not more than a few cm(-1).

Anharmonic vibrational properties by a fully automated second-order perturbative approach

This paper describes the implementation of a fully automated code for the building of anharmonic force constants and their use in a second-order perturbative evaluation of vibrorotational parameters. Next, a number of test applications are discussed, which show the strengths and limits of various computational levels.

Automatic generation of force fields and property surfaces for use in variational vibrational calculations of anharmonic vibrational energies and zero-point vibrational averaged properties

An automatic and general procedure for the calculation of geometrical derivatives of the energy and general property surfaces for molecular systems is developed and implemented. General expressions for an n-mode representation are derived, where the n-mode representation includes only the couplings between n or less degrees of freedom. The general expressions are specialized to derivative force fields and property surfaces, and a scheme for calculation of the numerical derivatives is implemented. The implementation is interfaced to electronic structure programs and may be used for both ground and excited electronic states. The implementation is done in the context of a vibrational structure program and can be used in combination with vibrational self-consistent field (VSCF), vibrational configuration interaction (VCI), vibrational Moller-Plesset, and vibrational coupled cluster calculations of anharmonic wave functions and calculation of vibrational averaged properties at the VSCF and VCI levels. Sample calculations are presented for fundamental vibrational energies and vibrationally averaged dipole moments and frequency dependent polarizabilities and hyperpolarizabilities of water and formaldehyde.

Franck-Condon factors based on anharmonic vibrational wave functions of polyatomic molecules

Franck-Condon (FC) integrals of polyatomic molecules are computed on the basis of vibrational self-consistent-field (VSCF) or configuration-interaction (VCI) calculations capable of including vibrational anharmonicity to any desired extent (within certain molecular size limits). The anharmonic vibrational wave functions of the initial and final states are expanded unambiguously by harmonic oscillator basis functions of normal coordinates of the respective electronic states. The anharmonic FC integrals are then obtained as linear combinations of harmonic counterparts, which can, in turn, be evaluated by established techniques taking account of the Duschinsky rotations, geometry displacements, and frequency changes. Alternatively, anharmonic wave functions of both states are expanded by basis functions of just one electronic state, permitting the FC integral to be evaluated directly by the Gauss-Hermite quadrature used in the VSCF and VCI steps [Bowman et al., Mol. Phys. 104, 33 (2006)]. These methods in conjunction with the VCI and coupled-cluster with singles, doubles, and perturbative triples [CCSD(T)] method have predicted the peak positions and intensities of the vibrational manifold in the X 2B1 photoelectron band of H2O with quantitative accuracy. It has revealed that two weakly visible peaks are the result of intensity borrowing from nearby states through anharmonic couplings, an effect explained qualitatively by VSCF and quantitatively by VCI, but not by the harmonic approximation. The X 2B2 photoelectron band of H2CO is less accurately reproduced by this method, likely because of the inability of CCSD(T)/cc-pVTZ to describe the potential energy surface of open-shell H2CO+ with the same high accuracy as in H2O+.

Highly accurate potential-energy and dipole moment surfaces for vibrational state calculations of methane

Full-dimensional ab initio potential-energy surface (PES) and dipole moment surface are constructed for a methane molecule at the CCSD(T)/cc-pVTZ and MP2/cc-pVTZ levels of theory, respectively, by the modified Shepard interpolation method based on the fourth-order Taylor expansion [MSI(4th)]. The reference points for the interpolation have been set in the coupling region of CH symmetric and antisymmetric stretching modes so as to reproduce the vibrational energy levels related to CH stretching vibrations. The vibrational configuration-interaction calculations have been performed to obtain the energy levels and the absorption intensities up to 9000 cm(-1) with the use of MSI(4th)-PES. The calculated fundamental frequencies and low-lying vibrational energy levels show that MSI(4th) is superior to the widely employed quartic force field, giving a better agreement with the experimental values. The absorption bands of overtones as well as combination bands, which are caused by purely anharmonic effects, have been obtained up to 9000 cm(-1). Strongly coupled states with visible intensity have been found in the 6500-9000 cm(-1) region where the experimental data are still lacking.

Prediction of accurate anharmonic experimental vibrational frequencies for water clusters, (H2O)n, n=2-5

Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2-5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.

Accurate vibrational spectra and magnetic properties of organic free radicals: The case of H2CN

We present the structural, magnetic, and vibrational properties of H2CN computed using a second-order perturbative approach in which equilibrium values and harmonic frequencies evaluated at the coupled-cluster level are combined with anharmonic and vibrational averaging contributions obtained by hybrid Hartree-Fock/Kohn-Sham methods. Our computations lead to remarkably accurate results and suggest some revision of the experimental vibrational assignments.

Internal rotation in peroxynitrous acid (ONOOH)

Using higher levels of wave-function-based electronic structure theory than previously applied, as well as density functional theory (B-LYP and B3-LYP functionals), all theoretical models conclude that three ONOOH conformers are stationary point minima, in disagreement with some of the previous studies that we survey. In order of increasing energy, these are the cis-cis, cis-perp, and trans-perp conformers. Basis sets including diffuse functions seem to be needed to obtain a qualitatively correct representation of the internal rotation potential energy surface at higher levels of theory. Internal rotation about the peroxide bond involving the cis-cis, cis-gauche transition structure (TS), cis-perp, and cis-trans TS conformers is studied in detail. To help ascertain the relative stability of the cis-perp conformer, multireference configuration interaction energy calculations are carried out, and rule of thumb estimates of multireference character in the ground-state wave functions of the ONOOH conformers are considered. CCSD(T)/aug-cc-pVTZ physical properties (geometries, rotational constants, electric dipole moments, harmonic vibrational frequencies, and infrared intensities) are compared with the analogous experimental data wherever possible, and also with density functional theory. Where such experimental data are nonexistent, the CCSD(T) and B3-LYP results are useful representations. For example, the electric dipole moment |mu(e)| of the cis-cis conformer is predicted to be 0.97+/-0.03 D. CCSD(T) energies, extrapolated to the aug-cc-pVNZ limit, are employed in isodesmic reaction schemes to derive zero Kelvin heats of formation and bond dissociation energies of the ONOOH stationary point minima. In agreement with recent gas-phase experiments, the peroxide bond dissociation energies of the cis-cis and trans-perp conformers are calculated as 19.3+/-0.4 and 16.0+/-0.4 kcalmol, respectively. The lowest energy cis-cis conformer is less stable than nitric acid by 28.1+/-0.4 kcalmol at 0 K.

Coupled-cluster connected quadruples and quintuples corrections to the harmonic vibrational frequencies and equilibrium bond distances of HF, N2, F2, and CO

Using the coupled-cluster method, we have examined the contributions from the connected quadruple and quintuple virtual excitations to the harmonic vibrational frequencies and equilibrium bond distances of HF, N(2), F(2), and CO. Whereas the largest quadruples contributions are -18.8 cm(-1) to the harmonic frequency of N(2) and 0.43 pm to the bond distance of F(2), the largest quintuples contributions are -3.9 cm(-1) to the harmonic frequency of N(2) and 0.03 pm to the bond distances of N(2) and F(2). As we improve the description by going from the coupled-cluster singles-and-doubles model with a perturbative triples correction [CCSD(T)] to the coupled-cluster singles-doubles-triples-and-quadruples model, the mean and maximum absolute errors in the calculated frequencies relative to experiment are reduced from 11.3 and 15.8 cm(-1), respectively, to 3.2 and 4.7 cm(-1); for the bond distances, the mean and maximum absolute errors are reduced from 0.16 and 0.47 pm, respectively, to 0.04 and 0.13 pm. The calculations presented here confirm previous observations that, to some extent, the relatively small errors in the CCSD(T) equilibrium bond distances and harmonic frequencies arise from a cancellation of errors in the approximate (perturbative) treatment of the connected triples and the neglect of higher-order connected excitations. Further inclusion of quintuples contributions, relativistic corrections, and adiabatic corrections reduces the mean and maximum absolute errors to 1.1 and 2.3 cm(-1), respectively, for the harmonic frequencies and to 0.02 and 0.05 pm, respectively, for the bond distances.

Ab initiovibrational state calculations with a quartic force field: Applications to H2CO, C2H4, CH3OH, CH3CCH, and C6H6

For polyatomic molecules, n-mode coupling representations of the quartic force field (nMR-QFF) are presented, which include terms up to n normal coordinate couplings in a fourth-order polynomial potential energy function. The computational scheme to evaluate third-and fourth-order derivatives by finite differentiations of the energy is fully described. The code to generate the nMR-QFF has been implemented into GAMESS program package and interfaced with the vibrational self-consistent field (VSCF) and correlation corrected VSCF (cc-VSCF) methods. As a demonstration, fundamental frequencies have been calculated by the cc-VSCF method based on 2MR-QFF for formaldehyde, ethylene, methanol, propyne, and benzene. The applications show that 2MR-QFF is a highly accurate potential energy function, with errors of 1.0-1.9% relative to the experimental value in fundamental frequencies. This approach will help quantitative evaluations of vibrational energies of a general molecule with a reasonable computational cost.