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.     

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.     

Anharmonic analysis of the vibrational spectrum of ketene by density functional theory using second-order perturbative approach

The paper reports main results of a comprehensive study of the vibrational spectrum of ketene computed using second-order perturbation theory treatment based on quartic, cubic and semidiagonal quartic force constants. Two different models--a homogeneous model using the same density functionals and basis functions for the harmonic calculations and anharmonic corrections, and a hybrid model in which the two parts of the calculation are conducted using different density functionals and basis sets--have been employed in the present calculations. Different DFT and CCSD methods and DZ and TZ extended basis sets involving diffuse and polarization functions have been used to calculate optimized and vibrationally averaged geometrical parameters, the harmonic and anharmonic vibrational frequencies and the spectroscopic constants such as anharmonicity constants, rotational constants, rotation-vibration coupling constants, Nielsen's centrifugal distortion constants and Coriolis coupling constants. Homogeneous model is found to be superior to the hybrid model in several respects. Difficulties in the hybrid model may arise due to one of the following reasons: (a) the basic requirement that the geometry optimization and frequency calculations must be done at the same level of theory to have valid frequencies is not met in the hybrid model; (b) the molecular structure gets reoptimized at the low level for anharmonic corrections; (c) in addition, the perturbation could also diverge for the above reasons, particularly for the very low, very anharmonic terms where the harmonic approximation is not close enough to make the perturbation work.     

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.     

Approaching the basis set limit for transition metal compounds with highly polar bonds: a benchmark coupled-cluster study of the ScF3 and FeF3 molecular structures and spectra

The molecular equilibrium geometries, quadratic and cubic force constants, vibrational frequencies, and infrared intensities of scandium and iron trifluorides are determined ab initio with a sequence of atomic natural orbital basis sets using the CCSD(T) treatment of electron correlation. The largest basis set of spdf ghi quality contains 462 contracted Gaussian functions. Relativistic corrections are applied to compute the equilibrium geometries and vibrational frequencies. The cubic force constants are used to estimate vibrational corrections to the effective r(g) internuclear distances determined in the gas electron diffraction experiments. The computed molecular properties are extrapolated to the complete basis-set limit. The predicted values are compared to the available experimental data; uncertainties and inconsistencies in these data are then discussed.     

Anharmonic effects in ammonium nitrate and hydroxylammonium nitrate clusters

The covalent and ionic clusters of ammonium nitrate and hydroxyl ammonium nitrate are characterized using density functional theory and second-order vibrational perturbation theory. The most stable structures are covalent acid-base pairs for the monomers and ionic acid-base pairs for the dimers. The hydrogen-bonding distances are greater in the ionic dimers than in the covalent monomers, and the stretching frequencies are significantly different in the covalent and ionic clusters. The anharmonicity of the potential energy surfaces is found to influence the geometries, frequencies, and nuclear magnetic shielding constants for these systems. The inclusion of anharmonic effects significantly decreases many of the calculated vibrational frequencies in these clusters and improves the agreement of the calculated frequencies with the experimental data available for the isolated neutral species. The calculations of nuclear magnetic shielding constants for all nuclei in these clusters illustrate that quantitatively accurate predictions of nuclear magnetic shieldings for comparison to experimental data require the inclusion of anharmonic effects. These calculations of geometries, frequencies, and shielding constants provide insight into the significance of anharmonic effects in ionic materials and provide data that will be useful for the parametrization of molecular mechanical force fields for ionic liquids. Anharmonic effects will be particularly important for the study of proton transfer reactions in ionic materials.     

Hydrogen bonding in the urea dimers and adenine–thymine DNA base pair: anharmonic effects in the intermolecular H-bond and intramolecular H-stretching vibrations

The equilibrium structures, binding energies, vibrational harmonic frequencies, and the anharmonic corrections for two different (cyclic and asymmetric) urea dimers and for the adenine–thymine DNA base pair system have been studied using the second-order Møller–Plesset perturbation theory (MP2) method and different density functional theory (DFT) exchange–correlation (XC) functionals (BLYP, B3LYP, PBE, HCTH407, KMLYP, and BH and HLYP) with the D95V, D95V**, and D95V++** basis sets. The widely used a posteriori Boys–Bernardi or counterpoise correction scheme for basis set superposition error (BSSE) has been included in the calculations to take into account the BSSE effects during geometry optimization (on structure), on binding energies and on the different levels of approximation used for calculating the vibrational frequencies. The results obtained with the ab initio MP2 method are compared with those calculated with different DFT XC functionals; and finally the suitability of these DFT XC functionals to describe intermolecular hydrogen bonds as well as harmonic frequencies and the anharmonic corrections is assessed and discussed.     

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

CVRQD ab initio ground-state adiabatic potential energy surfaces for the water molecule

The high accuracy ab initio adiabatic potential energy surfaces (PESs) of the ground electronic state of the water molecule, determined originally by Polyansky et al. [Science 299, 539 (2003)] and called CVRQD, are extended and carefully characterized and analyzed. The CVRQD potential energy surfaces are obtained from extrapolation to the complete basis set of nearly full configuration interaction valence-only electronic structure computations, augmented by core, relativistic, quantum electrodynamics, and diagonal Born-Oppenheimer corrections. We also report ab initio calculations of several quantities characterizing the CVRQD PESs, including equilibrium and vibrationally averaged (0 K) structures, harmonic and anharmonic force fields, harmonic vibrational frequencies, vibrational fundamentals, and zero-point energies. They can be considered as the best ab initio estimates of these quantities available today. Results of first-principles computations on the rovibrational energy levels of several isotopologues of the water molecule are also presented, based on the CVRQD PESs and the use of variational nuclear motion calculations employing an exact kinetic energy operator given in orthogonal internal coordinates. The variational nuclear motion calculations also include a simplified treatment of nonadiabatic effects. This sophisticated procedure to compute rovibrational energy levels reproduces all the known rovibrational levels of the water isotopologues considered, H(2) (16)O, H(2) (17)O, H(2) (18)O, and D(2) (16)O, to better than 1 cm(-1) on average. Finally, prospects for further improvement of the ground-state adiabatic ab initio PESs of water are discussed.     

Exploring the effect of anharmonicity of molecular vibrations on thermodynamic properties

Thermodynamic properties of selected small and medium size molecules were calculated using harmonic and anharmonic vibrational frequencies. Harmonic vibrational frequencies were obtained by normal mode analysis, whereas anharmonic ones were calculated using the vibrational self-consistent field (VSCF) method. The calculated and available experimental thermodynamic data for zero point energy, enthalpy, entropy, and heat capacity are compared. It is found that the anharmonicity and coupling of molecular vibrations can play a significant role in predicting accurate thermodynamic quantities. Limitations of the current VSCF method for low frequency modes have been partially removed by following normal mode displacements in internal, rather than Cartesian, coordinates.     

BSSE‐corrected geometry and harmonic and anharmonic vibrational frequencies of formamide–water and formamide–formamide dimers

The basis set superposition error (BSSE) influence in the geometry structure, interaction energies, and intermolecular harmonic and anharmonic vibrational frequencies of cyclic formamide–formamide and formamide–water dimers have been studied using different basis sets (6‐31G, 6‐31G**, 6‐31++G**, D95V, D95V**, and D95V++**). The a posteriori “counterpoise” (CP) correction scheme has been compared with the a priori “chemical Hamiltonian approach” (CHA) both at the Hartree–Fock (HF) and second‐order Møller–Plesset many‐body perturbation (MP2) levels of theory. The effect of BSSE on geometrical parameters, interaction energies, and intermolecular harmonic vibrational frequencies are discussed and compared with the existing experimental data. As expected, the BSSE‐free CP and CHA interaction energies usually show less deep minima than those obtained from the uncorrected methods at both the HF and MP2 levels. Focusing on the correlated level, the amount of BSSE in the intermolecular interaction energies is much larger than that at the HF level, and this effect is also conserved in the values of the force constants and harmonic vibrational frequencies. All these results clearly indicate the importance of the proper BSSE‐free correlation treatment with the well‐defined basis functions. At the same time, the results show a good agreement between the a priori CHA and a posteriori CP correction scheme; this agreement is remarkable in the case of large and well‐balanced basis sets. The anharmonic frequency correction values also show an important BSSE dependence, especially for hydrogen bond stretching and for low frequencies belonging to the intermolecular normal modes. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Anharmonic rovibrational analysis for disilacyclopropenylidene (Si2CH2)

The global minimum on the Si(2)CH(2) electronic singlet potential energy surface has been theoretically predicted to be a peculiar hydrogen bridged (Si···H···Si) disilacyclopropenylidene structure (Si(2)CH(2)). An accurate quartic force field for Si(2)CH(2) has been determined employing ab initio coupled-cluster theory with single and double excitations and a perturbative treatment for triple excitations [CCSD(T)], in combination with the correlation consistent core-valence quadruple zeta (cc-pCVQZ) basis set. The vibration-rotation coupling constants, equilibrium and zero-point vibration corrected rotational constants, centrifugal distortion constants, and harmonic and fundamental vibrational frequencies for six isotopologues of Si(2)CH(2) are predicted using vibrational second-order perturbation theory (VPT2). The anharmonic corrections for the vibrational motions involving the H bridged bonds are found to be more than 5% with respect to the corresponding harmonic vibrational frequencies. In this light, an experimental detection and characterization of disilacyclopropenylidene (Si(2)CH(2)) is highly desired.     

A high‐accuracy theoretical study of the CHnP Systems n = 1–3

We have performed high‐level electronic structure computations on the most important species of the CH_nP systems n = 1–3 to characterize them and provide reliable information about the equilibrium and vibrationally averaged molecular structures, rotational constants, vibrational frequencies (harmonic and anharmonic), formation enthalpies, and vertical excitation energies. Those chemical systems are intermediates for several important reactions and also prototypical phosphorus‐carbon compounds; however, they are often elusive to experimental detection. The present results significantly complement their knowledge and can be used as an assessment of the experimental information when available. The explicitly correlated coupled‐cluster RCCSD(T)‐F12 method has been used for geometry optimizations and vibrational frequency calculations. Vibrational configuration interaction theory has been used to account for anharmonicity effects. Basis‐set limit extrapolations have been carried out to determine accurate thermochemical quantities. Electronic excited states have been calculated with coupled‐cluster approaches and also by means of the multireference configuration interaction method. © 2013 Wiley Periodicals, Inc.     

Molecular structures of the two most stable conformers of free glycine

The equilibrium molecular structures of the two lowest-energy conformers of glycine, Gly-Ip and Gly-IIn, have been characterized by high-level ab initio electronic structure computations, including all-electron cc-pVTZ CCSD(T) geometry optimizations and 6-31G* MP2 quartic force fields, the latter to account for anharmonic zero-point vibrational effects to isotopologic rotational constants. Based on experimentally measured vibrationally averaged effective rotational constant sets of several isotopologues and our ab initio data for structural constraints and zero-point vibrational shifts, least-squares structural refinements were performed to determine improved Born-Oppenheimer equilibrium (r(e)) structures of Gly-Ip and Gly-IIn. Without the ab initio constraints even the extensive set of empirical rotational constants available for 5 and 10 isotopologues of Gly-Ip and Gly-IIn, respectively, cannot satisfactorily fix their molecular structure. Excellent agreement between theory and experiment is found for the rotational constants of both conformers, the rms residual of the final fits being 7.8 and 51.6 kHz for Gly-Ip and Gly-IIn, respectively. High-level ab initio computations with focal point extrapolations determine the barrier to planarity separating Gly-IIp and Gly-IIn to be 20.5 +/- 5.0 cm(-1). The equilibrium torsion angle tau(NCCO) of Gly-IIn, characterizing the deviation of its heavy-atom framework from planarity, is (11 +/- 2) degrees. Nevertheless, in the ground vibrational state the effective structure of Gly-IIn has a plane of symmetry.     

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.     

Protonated carbonyl sulfide: prospects for the spectroscopic observation of the elusive HSCO+ isomer

Spurred by the apparent conflict between ab initio predictions and infrared spectroscopic evidence regarding the relative stability of isomers of protonated carbonyl sulfide, key stationary points on the isomerization surface of HOCS(+) have been examined via systematic extrapolations of ab initio energies. Electron correlation has been accounted for using second-order M?ller-Plesset perturbation theory and coupled cluster theory through triple excitations [CCSD, CCSD(T), and CCSDT] in conjunction with the correlation consistent hierarchy of basis sets, cc-pVXZ (X=D,T,Q,5,6). HSCO(+) is predicted to lie lower in energy than HOCS(+) by 4.86 kcal mol(-1), computed using the focal point extrapolation scheme of Allen and co-workers [J. Chem. Phys. 99, 4638 (1993)] with corrections for anharmonic zero-point vibrational energy, core correlation, non-Born-Oppenheimer, and scalar relativistic effects. A transition state has been located, constituting the barrier to isomerization of HSCO(+) to HOCS(+), lying 68.9 kcal mol(-1) higher in energy than HSCO(+). This is well above predicted exothermicity [DeltaH(r) (o)(0 K)=48.1 kcal mol(-1), cc-pVQZ CCSD(T)] for the reaction considered in the experiments (HSCO(+)+H(2)-->OCS+H(3) (+)). Though proton tunneling will lead to a lower effective barrier, this prediction is consistent with the lack of HSCO(+) in electrical discharges in H(2)OCS, since the relative populations of HOCS(+) and HSCO(+) will depend on the experimental details of the protonation route rather than the relative thermodynamic stability of the isomers. Anharmonic vibrational frequencies and vibrationally corrected rotational constants from cc-pVTZ CCSD(T) cubic and quartic force constants are provided, to aid in the spectroscopic observation of the energetically favorable but apparently elusive HSCO(+) isomer.     

Vibrational spectroscopy of protonated imidazole and its complexes with water molecules: ab initio anharmonic calculations and experiments

The results of anharmonic frequency calculations on neutral imidazole (C3N2H4, Im), protonated imidazole (ImH+), and its complexes with water (ImH+)(H2O)n, are presented and compared to gas phase infrared photodissociation spectroscopy (IRPD) data. Anharmonic frequencies are obtained via ab initio vibrational self-consistent field (VSCF) calculations taking into account pairwise interactions between the normal modes. The key results are: (1) Prediction of anharmonic vibrational frequencies on an MP2 ab initio potential energy surface show excellent agreement with experiment and outstanding improvement over the harmonic frequencies. For example, the ab initio calculated anharmonic frequency for (ImH+)(H2O)N2 exhibits an overall average percentage error of 0.6% from experiment. (2) Anharmonic vibrational frequencies calculated on a semiempirical potential energy surface fitted to ab initio harmonic data represents spectroscopy well, particularly for water complexes. As an example, anharmonic frequencies for (ImH+)H2O and (ImH+)(H2O)2 show an overall average deviation of 1.02% and 1.05% from experiment, respectively. This agreement between theory and experiment also supports the validity and use of the pairwise approximation used in the calculations. (3) Anharmonic coupling due to hydration effects is found to significantly reduce the vibrational frequencies for the NH stretch modes. The frequency of the NH stretch is observed to increase with the removal of a water molecule or replacement of water with N2. This result also indicates the ability of the VSCF method to predict accurate frequencies in a matrix environment. The calculation provides insights into the nature of anharmonic effects in the potential surface. Analysis of percentage anharmoncity in neutral Im and ImH+ shows a higher percentage anharmonicity in the NH and CH stretch modes of neutral Im. Also, we observe that anharmonicity in the NH stretch modes of ImH+ have some contribution from coupling effects, while that of neutral Im has no contribution whatsoever from mode-mode coupling. It is concluded that the incorporation of anharmonic effects in the calculation brings theory and experiment into much closer agreement for these systems.     

Ab initio vibrational 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.     

New vibrational self-consistent field program for large molecules

A recently developed, general computer program that performs vibrational self-consistent field (VSCF) calculations for large molecules is described. The program, which we refer to as VSCF―95, requires as its only input a force field in mass-scaled normal coordinates. Currently, it is limited to a maximum of 200 normal modes, and the force field is limited to coupling terms involving a maximum of six normal modes, with a maximum order of six in any normal mode. As output the program returns VSCF energies for specified quantum states. We illustrate the code with two new applications. The first is to HCO, for which we use a full sixth-order force field. The second is to a model of the fullerene, C₆₀, for which we have calculated a 75,731-term force field, which includes all anharmonic terms up to fifth order, and all two-mode coupling terms up to fourth order. © 1996 by John Wiley & Sons, Inc.     

A quartic force field coordinate substitution scheme using hyperbolic sine coordinates

Quartic force fields (QFF) are currently the most cost‐effective method for the approximation of potential energy surfaces for the calculation of anharmonic vibrational energies. It is known, although, that its performance can be less than satisfactory due to limitations related to slow convergence of the series. In this article, we present a coordinate substitution scheme using a combination of Morse and sinh coordinates, well adapted for its use with cartesian normal coordinates. We derive expressions for analytical integrals for use in VSCF and VCI calculations and show that the simultaneous substitution of symmetric and antisymmetric normal coordinates by Morse and sinh coordinates, respectively, significantly improves the vibrational transition frequencies for these modes in a well‐balanced fashion. The accuracy of this substitution scheme is demonstrated by comparing one and two‐dimensional sections of substituted and unsubstituted QFF with ab initio potential energy grids, as well as with vibrational energy calculations using as test cases two well‐studied benchmark molecules: water and formaldehyde. We conclude that the coordinate substitution scheme presented constitutes a very attractive alternative to simple QFFs in the context of cartesian normal coordinates.