Vibrational structure theory: new vibrational wave function methods for calculation of anharmonic vibrational energies and vibrational contributions to molecular properties

A number of recently developed theoretical methods for the calculation of vibrational energies and wave functions are reviewed. Methods for constructing the appropriate quantum mechanical Hamilton operator are briefly described before reviewing a particular branch of theoretical methods for solving the nuclear Schr?dinger equation. The main focus is on wave function methods using the vibrational self-consistent field (VSCF) as starting point, and includes vibrational configuration interaction (VCI), vibrational M?ller-Plesset (VMP) theory, and vibrational coupled cluster (VCC) theory. The convergence of the different methods towards the full vibrational configuration interaction (FVCI) result is discussed. Finally, newly developed vibrational response methods for calculation of vibrational contributions to properties, energies, and transition probabilities are discussed.     

Response theory for vibrational wave functions

A formalism for deriving and implementing response functions for vibrational wave functions is described. The formalism utilizes a recently developed second-quantization formulation of many-mode dynamics to define nonredundant parameterizations for different types of approximate vibrational wave functions. The derived response functions cover the cases of an exact state, a vibrational self-consistent field state, and a vibrational configuration interaction state. Sample calculations are presented for the linear-response function and response excitation energies for a two-mode model system and for formaldehyde employing a quartic force field. The advantages and disadvantages of the response theoretical approach for describing excitation energies using different parameterizations are discussed.     

Vibrational absorption spectra from vibrational coupled cluster damped linear response functions calculated using an asymmetric Lanczos algorithm

We report the theory and implementation of vibrational coupled cluster (VCC) damped response functions. From the imaginary part of the damped VCC response function the absorption as function of frequency can be obtained, requiring formally the solution of the now complex VCC response equations. The absorption spectrum can in this formulation be seen as a matrix function of the characteristic VCC Jacobian response matrix. The asymmetric matrix version of the Lanczos method is used to generate a tridiagonal representation of the VCC response Jacobian. Solving the complex response equations in the relevant Lanczos space provides a method for calculating the VCC damped response functions and thereby subsequently the absorption spectra. The convergence behaviour of the algorithm is discussed theoretically and tested for different levels of completeness of the VCC expansion. Comparison is made with results from the recently reported [P. Seidler, M. B. Hansen, W. Gyo?rffy, D. Toffoli, and O. Christiansen, J. Chem. Phys. 132, 164105 (2010)] vibrational configuration interaction damped response function calculated using a symmetric Lanczos algorithm. Calculations of IR spectra of oxazole, cyclopropene, and uracil illustrate the usefulness of the new VCC based method.     

Vibrational coupled cluster theory

The theory and first implementation of a vibrational coupled cluster (VCC) method for calculations of the vibrational structure of molecules is presented. Different methods for introducing approximate VCC methods are discussed including truncation according to a maximum number of simultaneous mode excitations as well as an interaction space order concept is introduced. The theory is tested on calculation of anharmonic frequencies for a three-mode model system and a formaldehyde quartic force field. The VCC method is compared to vibrational self-consistent-field, vibrational M?ller-Plesset perturbation theory, and vibrational configuration interaction (VCI). A VCC calculation typically gives higher accuracy than a corresponding VCI calculation with the same number of parameters and the same formal operation count.     

Approximate inclusion of four-mode couplings in vibrational coupled-cluster theory

The vibrational coupled cluster (VCC) equations are analyzed in terms of vibrational Mo?ller-Plesset perturbation theory aiming specifically at the importance of four-mode couplings. Based on this analysis, new VCC methods are derived for the calculation of anharmonic vibrational energies and vibrational spectra using vibrational coupled cluster response theory. It is shown how the effect of four-mode coupling and excitations can be efficiently and accurately described using approximations for their inclusion. Two closely related approaches are suggested. The computational scaling of the so-called VCC[3pt4F] method is not higher than the fifth power in the number of vibrational degrees of freedom when up to four-mode coupling terms are present in the Hamiltonian and only fourth order when only up to three-mode couplings are present. With a further approximation, one obtains the VCC[3pt4] model which is shown to scale with at most the fourth power in the number of vibrational degrees of freedom for Hamiltonians with both three- and four-mode coupling levels, while sharing the most important characteristics with VCC[3pt4F]. Sample calculations reported for selected tetra-atomic molecules as well as the larger dioxirane and ethylene oxide molecules support that the new models are accurate and useful.     

Vibrational coupled cluster response theory: a general implementation

The calculation of vibrational contributions to molecular properties using vibrational coupled cluster (VCC) response theory is discussed. General expressions are given for expectation values, linear response functions, and transition moments. It is shown how these expressions can be evaluated for arbitrary levels of excitation in the wave function parameterization as well as for arbitrary coupling levels in the potential and property surfaces. The convergence of the method is assessed by benchmark calculations on formaldehyde. Furthermore, excitation energies and infrared intensities are calculated for the fundamental vibrations of furan using VCC limited to up to two-mode and up to three-mode excitations, VCC[2] and VCC[3], as well as VCC with full two-mode and approximate three-mode couplings, VCC[2pt3].     

Impact of nuclear quantum effects on the molecular structure of bihalides and the hydrogen fluoride dimer

The structural impact of nuclear quantum effects is investigated for a set of bihalides, [XHX](-), X = F, Cl, and Br, and the hydrogen fluoride dimer. Structures are calculated with the vibrational self-consistent-field (VSCF) method, the second-order vibrational perturbation theory method (VPT2), and the nuclear-electronic orbital (NEO) approach. In the VSCF and VPT2 methods, the vibrationally averaged geometries are calculated for the Born-Oppenheimer electronic potential energy surface. In the NEO approach, the hydrogen nuclei are treated quantum mechanically on the same level as the electrons, and mixed nuclear-electronic wave functions are calculated variationally with molecular orbital methods. Electron-electron and electron-proton dynamical correlation effects are included in the NEO approach using second-order perturbation theory (NEO-MP2). The nuclear quantum effects are found to alter the distances between the heavy atoms by 0.02-0.05 A for the systems studied. These effects are of similar magnitude as the electron correlation effects. For the bihalides, inclusion of the nuclear quantum effects with the NEO-MP2 or the VSCF method increases the X-X distance. The bihalide X-X distances are similar for both methods and are consistent with two-dimensional grid calculations and experimental values, thereby validating the use of the computationally efficient NEO-MP2 method for these types of systems. For the hydrogen fluoride dimer, inclusion of nuclear quantum effects decreases the F-F distance with the NEO-MP2 method and increases the F-F distance with the VSCF and VPT2 methods. The VPT2 F-F distances for the hydrogen fluoride dimer and the deuterated form are consistent with the experimentally determined values. The NEO-MP2 F-F distance is in excellent agreement with the distance obtained experimentally for a model that removes the large amplitude bending motions. The analysis of these calculations provides insight into the significance of electron-electron and electron-proton correlation, anharmonicity of the vibrational modes, and nonadiabatic effects for hydrogen-bonded systems.     

Anharmonic overtone and combination states of glycine and two model peptides examined by vibrational self-consistent field theory

In this paper, the application of the vibrational self-consistent field (VSCF) and correction-corrected VSCF methods for calculating anharmonic parameters, including transition frequency, transition intensity and dipole, and vibrational anharmonicity of 3N-6 normal modes for formamide, glycine, N-methylacetamide and their deuterated derivatives are explored mainly at the level of density functional theory. The computed fundamental anharmonic frequencies are found to be in reasonable agreement with experimental results. Diagonal anharmonicities of the second overtone states were examined for multiple normal modes, whose magnitudes were found to correlate well with those of the first overtone states in the three small molecules. The results show that the VSCF-based approach can be utilized to predict anharmonic parameters of higher vibrational states that are essential to understanding multi-pulse infrared nonlinear experiments of peptides.     

Calculating molecular vibrational spectra beyond the harmonic approximation

We present a new approach for calculating anharmonic corrections to vibrational frequency calculations. The vibrational wavefunction is modelled using translated Hermite functions thus allowing anharmonic effects to be incorporated directly into the wavefunction whilst still retaining the simplicity of the Hermite basis. We combine this new method with an optimised finite-difference grid for computing the necessary third and fourth nuclear derivatives of the energy. We compare our combined approach to existing anharmonic models—vibrational self-consistent field theory (VSCF), vibrational perturbation theory (VPT), and vibrational configuration interaction theory (VCI)—and find that it is more cost effective than these alternatives. This makes our method well-suited to computing anharmonic corrections for frequencies in medium-sized molecules. Contribution of the Mark S. Gordon 65th Birthday Festschrift Issue.     

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

Multiresolution potential energy surfaces for vibrational state calculations

A compact and robust many-mode expansion of potential energy surfaces (PES) is presented for anharmonic vibrations of polyatomic molecules, where the individual many-mode terms are approximated with various different resolutions, i.e., electronic structure methods, basis sets, and functional forms. As functional forms, the following three representations have been explored: numerical values on a grid, cubic spline interpolation, and a Taylor expansion. A useful index is proposed which rapidly identifies important many-mode terms that warrant a high resolution. Applications to water and formaldehyde demonstrate that the present scheme can increase the efficiency of the PES computation by a factor of up to 11 with the errors in anharmonic vibrational frequencies being no worse than ~ 10cm⁻¹.     

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

The VMFCI method: a flexible tool for solving the molecular vibration problem

The present article introduces a general variational scheme to find approximate solutions of the spectral problem for the molecular vibration Hamiltonian. It is called the "vibrational mean field configuration interaction" (VMFCI) method, and consists in performing vibrational configuration interactions (VCI) for selected modes in the mean field of the others. The same partition of modes can be iterated until self-consistency, generalizing the vibrational self-consistent field (VSCF) method. As in contracted-mode methods, a hierarchy of partitions can be built to ultimately contract all the modes together. So, the VMFCI method extends the traditional variational approaches and can be included in existing vibrational codes based on the latter approaches. The flexibility and efficiency of this new method are demonstrated on several molecules of atmospheric interest.     

Semiclassical many-mode floquet theory

It is shown that the single-mode Floquet formalism of Shirley can be extended to a generalized many-mode Floquet theory, yielding a practical non-perturbative technique for the semiclassical treatment of the interaction of a quantum system several monochromatic oscillating fields. The theory is illustrated by a detailed study of the population dynamics of a three-level system driven by two monochromatic radiation fields.     

Determination of molecular vibrational state energies using the ab initio semiclassical initial value representation: application to formaldehyde

We have demonstrated the use of ab initio molecular dynamics (AIMD) trajectories to compute the vibrational energy levels of molecular systems in the context of the semiclassical initial value representation (SC-IVR). A relatively low level of electronic structure theory (HF/3-21G) was used in this proof-of-principle study. Formaldehyde was used as a test case for the determination of accurate excited vibrational states. The AIMD-SC-IVR vibrational energies have been compared to those from curvilinear and rectilinear vibrational self-consistent field/vibrational configuration interaction with perturbation selected interactions-second-order perturbation theory (VSCF/VCIPSI-PT2) and correlation-corrected vibrational self-consistent field (cc-VSCF) methods. The survival amplitudes were obtained from selecting different reference wavefunctions using only a single set of molecular dynamics trajectories. We conclude that our approach is a further step in making the SC-IVR method a practical tool for first-principles quantum dynamics simulations.     

Size-extensive vibrational self-consistent field method

The vibrational self-consistent field (VSCF) method is a mean-field approach to solve the vibrational Schro?dinger equation and serves as a basis of vibrational perturbation and coupled-cluster methods. Together they account for anharmonic effects on vibrational transition frequencies and vibrationally averaged properties. This article reports the definition, programmable equations, and corresponding initial implementation of a diagrammatically size-extensive modification of VSCF, from which numerous terms with nonphysical size dependence in the original VSCF equations have been eliminated. When combined with a quartic force field (QFF), this compact and strictly size-extensive VSCF (XVSCF) method requires only quartic force constants of the ?(4)V/?Q(i)(2)?Q(j)(2) type, where V is the electronic energy and Q(i) is the ith normal coordinate. Consequently, the cost of a XVSCF calculation with a QFF increases only quadratically with the number of modes, while that of a VSCF calculation grows quartically. The effective (mean-field) potential of XVSCF felt by each mode is shown to be harmonic, making the XVSCF equations subject to a self-consistent analytical solution without matrix diagonalization or a basis-set expansion, which are necessary in VSCF. Even when the same set of force constants is used, XVSCF is nearly three orders of magnitude faster than VSCF implemented similarly. Yet, the results of XVSCF and VSCF are shown to approach each other as the molecular size is increased, implicating the inclusion of unnecessary, nonphysical terms in VSCF. The diagrams of the XVSCF energy expression and their evaluation rules are also proposed, underscoring their connected structures.     

Vibrational quasi-degenerate perturbation theory: applications to fermi resonance in CO2, H2CO, and C6H6

A quasi-degenerate perturbation method with vibrational self-consistent field (VSCF) reference wavefunction is developed. It simultaneously accounts for strong anharmonic mode-mode coupling among a few states (static correlation) by a configuration interaction theory and for weak coupling with a vast number of the other states (dynamic correlation) by a perturbation theory. A general formula is derived based on the van Vleck perturbation theory. An algorithm that selects a compact set of the most important VSCF configurations which contribute to the static correlation is proposed and a scheme to limit the number of configurations considered for dynamic correlation is also implemented. This method reproduces the vibrational frequencies of CO2 and H2CO that are subject to the strongest anharmonic mode-mode coupling within 10 cm(-1) of vibrational configuration interaction results in a computational expense reduced by a factor of one to two orders of magnitude. The method also reproduces the infrared absorption of C6H6 in the CH stretching (nu12) frequency region, in which combination tones nu13nu16 and nu2nu13nu18 appear on account of an intensity borrowing from nu12via the anharmonic coupling.     

Anharmonic vibrational spectroscopy and investigation of intramolecular mode couplings in adenine

Vibrational frequencies for the nucleobase adenine are calculated by the vibrational self-consistent field (VSCF) and correlation corrected vibrational self-consistent field (CC-VSCF) methods using Hartree-Fock (HF), density functional theory (DFT) and second order Møller-Plesset (MP2) theories. A large number of potential energy surface (PES) points were computed in the anharmonic calculations corresponding to each method. The quartic force field (QFF) approximation was used to generate the full grid of points for the VSCF solver. We have implemented our new procedure for computing the mode-mode coupling integrals in the 2-mode coupling representations of the quartic force field (2MR-QFF) for prediction of coupling magnitudes. Calculations were performed using the 6-31G(d,p) basis set. Comparison of the calculated ab initio anharmonic spectra with Ar matrix experimental data of adenine reported in the literature reveals that, the CC-VSCF (DFT) wavenumbers show the best agreement. The experimental geometric parameters of adenine are compared with the theoretically optimized molecular structural parameters. These are found to be in good agreement. Vibrational assignments are based on the calculated potential energy distribution (PED) values.