Molecular second-quantized Hamiltonian: Electron correlation and non-adiabatic coupling treated on an equal footing.

We introduce a new theoretical and computational framework for treating molecular quantum mechanics without the Born-Oppenheimer approximation. The molecular wavefunction is represented in a tensor-product space of electronic and vibrational basis functions, with electronic basis chosen to reproduce the mean-field electronic structure at all geometries. We show how to transform the Hamiltonian to a fully second-quantized form with creation/annihilation operators for electronic and vibrational quantum particles, paving the way for polynomial-scaling approximations to the tensor-product space formalism. In addition, we make a proof-of-principle application of the new Ansatz to the vibronic spectrum of C2.

The Molpro quantum chemistry package.

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.

Comparing the accuracy of perturbative and variational calculations for predicting fundamental vibrational frequencies of dihalomethanes.

Three dihalogenated methane derivatives (CH2F2, CH2FCl, and CH2Cl2) were used as model systems to compare and assess the accuracy of two different approaches for predicting observed fundamental frequencies: canonical operator Van Vleck vibrational perturbation theory (CVPT) and vibrational configuration interaction (VCI). For convenience and consistency, both methods employ the Watson Hamiltonian in rectilinear normal coordinates, expanding the potential energy surface (PES) as a Taylor series about equilibrium and constructing the wavefunction from a harmonic oscillator product basis. At the highest levels of theory considered here, fourth-order CVPT and VCI in a harmonic oscillator basis with up to 10 quanta of vibrational excitation in conjunction with a 4-mode representation sextic force field (SFF-4MR) computed at MP2/cc-pVTZ with replacement CCSD(T)/aug-cc-pVQZ harmonic force constants, the agreement between computed fundamentals is closer to 0.3 cm-1 on average, with a maximum difference of 1.7 cm-1. The major remaining accuracy-limiting factors are the accuracy of the underlying electronic structure model, followed by the incompleteness of the PES expansion. Nonetheless, computed and experimental fundamentals agree to within 5 cm-1, with an average difference of 2 cm-1, confirming the utility and accuracy of both theoretical models. One exception to this rule is the formally IR-inactive but weakly allowed through Coriolis-coupling H-C-H out-of-plane twisting mode of dichloromethane, whose spectrum we therefore revisit and reassign. We also investigate convergence with respect to order of CVPT, VCI excitation level, and order of PES expansion, concluding that premature truncation substantially decreases accuracy, although VCI(6)/SFF-4MR results are still of acceptable accuracy, and some error cancellation is observed with CVPT2 using a quartic force field.

The PyPES library of high quality semi-global potential energy surfaces.

In this article, we present a Python-based library of high quality semi-global potential energy surfaces for 50 polyatomic molecules with up to six atoms. We anticipate that these surfaces will find widespread application in the testing of new potential energy surface construction algorithms and nuclear ro-vibrational structure theories. To this end, we provide the ability to generate the energy derivatives required for Taylor series expansions to sixth order about any point on the potential energy surface in a range of common coordinate systems, including curvilinear internal, Cartesian, and normal mode coordinates. The PyPES package, along with FORTRAN, C, MATLAB and Mathematica wrappers, is available at http://sourceforge.net/projects/pypes-lib.

Quadratic Corrections to Harmonic Vibrational Frequencies Outperform Linear Models.

Simulating accurate infrared spectra is a longstanding problem in computational quantum chemistry. Linearly scaling harmonic frequencies to better match experimental data is a popular way of approximating anharmonic effects while simultaneously attempting to account for deficiencies in ab initio method and/or basis set. As this approach is empirical, it is also nonvariational and unbounded, so it is important to separate and quantify errors as robustly as possible. Eliminating the confounding factor of methodological incompleteness enables us to explore the intrinsic accuracy of the scaling approach alone. We find that single-coefficient linear scaling methods systematically overcorrect low frequencies, while generally undercorrecting higher frequencies. A two-parameter polynomial model gives significantly better predictions without systematic bias in any spectral region, while a single-parameter quadratic scaling model is parameterized to minimize overcorrection errors while only slightly decreasing predictive power.