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.

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.

The PyVCI program package is a general purpose open-source code for simulating accurate molecular spectra, based upon force field expansions of the potential energy surface in normal mode coordinates. It includes harmonic normal coordinate analysis and vibrational configuration interaction (VCI) algorithms, implemented primarily in Python for accessibility but with time-consuming routines written in C. Coriolis coupling terms may be optionally included in the vibrational Hamiltonian. Non-negligible VCI matrix elements are stored in sparse matrix format to alleviate the diagonalization problem. CPU and memory requirements may be further controlled by algorithmic choices and/or numerical screening procedures, and recommended values are established by benchmarking using a test set of 44 molecules for which accurate analytical potential energy surfaces are available.Force fields in normal mode coordinates are obtained from the PyPES library of high quality analytical potential energy surfaces (to 6th order) or by numerical differentiation of analytic second derivatives generated using the GAMESS quantum chemical program package (to 4th order).

This work describes the benchmarking of a VCI algorithm that combines the favourable computational scaling of VPT2 with the algorithmic robustness of VCI, in which VCI basis states are selected according to the magnitude of their contribution to the VPT2 energy, for the ground state and fundamental excited states. Particularly novel aspects of this work include: expanding the potential to 6th order in normal mode coordinates, using a double-iterative procedure in which configuration selection and VCI wavefunction updates are performed iteratively (micro-iterations) over a range of screening threshold values (macro-iterations), and characterisation of computational resource requirements as a function of molecular size. Computational costs may be further reduced by a priori truncation of the VCI wavefunction according to maximum extent of mode coupling, along with discarding negligible force constants and VCI matrix elements, and formulating the wavefunction in a harmonic oscillator product basis to enable efficient evaluation of VCI matrix elements. Combining these strategies, we define a series of screening procedures that scale as O(N 6 mode ) − O(N 9 mode ) in run time and O(N 6 mode ) − O(N 7 mode ) in memory, depending on the desired level of accuracy. Our open-source code is freely available for download from http://www.sourceforge.net/projects/pyvci-vpt2.

In this paper, we outline a general, scalable, and black-box approach for calculating high-order strongly coupled force fields in rectilinear normal mode coordinates, based upon constructing low order expansions in curvilinear coordinates with naturally limited mode-mode coupling, and then transforming between coordinate sets analytically. The optimal balance between accuracy and efficiency is achieved by transforming from 3 mode representation quartic force fields in curvilinear normal mode coordinates to 4 mode representation sextic force fields in rectilinear normal modes. Using this reduced mode-representation strategy introduces an error of only 1 cm(-1) in fundamental frequencies, on average, across a sizable test set of molecules. We demonstrate that if it is feasible to generate an initial semi-quartic force field in curvilinear normal mode coordinates from ab initio data, then the subsequent coordinate transformation procedure will be relatively fast with modest memory demands. This procedure facilitates solving the nuclear vibrational problem, as all required integrals can be evaluated analytically. Our coordinate transformation code is implemented within the extensible PyPES library program package, at http://sourceforge.net/projects/pypes-lib-ext/.

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 non-variational 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.

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 C 2 .

Three dihalogenated methane derivatives (CHF, CHFCl, and CHCl) 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 on average, with a maximum difference of 1.7 cm. 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, with an average difference of 2 cm, 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.