The vibrational auto-adjusting perturbation theory

In this work a new method to calculate anharmonic vibrational ground and excited state energies is proposed. The method relies on the auto-adjusting perturbation theory (APT) which has been successfully used to diagonalize square matrices. We use as zeroth order correction the self-consistent vibrational energies, and with the APT approach we calculate the vibrational anharmonic correlation correction to any desired order. In this paper we present the methodology and apply it to a model system and formaldehyde. Vibrational APT approach shows a robust convergent behavior even for the states where the standard (Rayleigh-Schrödinger) vibrational Møller-Plesset perturbation theory is clearly divergent.     

Iterative natural orbitals for configuration interaction using perturbation theory

Approximate natural orbitals are determined iteratively from CI expansions constructed using first-order perturbation theory in order to investigate the possibility of eliminating the complete transformation of MO integrals on each iteration. Results on LiH and H₂O are compared with fully variationally determined NO's to assess questions of convergence.     

Extended Koopmans' theorem at the second-order perturbation theory

The extended Koopmans' theorem (EKT), when combined with the second-order Møller−Plesset (MP2) perturbation theory through the relaxed density matrix approach [J. Cioslowski, P. Piskorz, and G. Liu, J. Chem. Phys. 1997, 107, 6,804], provides a straightforward way to calculate the ionization potentials (IPs) as an one electron quantity. However, such an EKT-MP2 method often suffers from the negative occupation problem, failing to provide the complete IP spectra for a system of interest. Here a small positive number scheme is proposed to cure this problem so as to remove the associated unphysical results. In order to obtain an in-depth physical interpretation of the EKT-MP2 method, we introduce a Koopmans-type quantity, named KT-MP2, based on which the respective contribution from the relaxation and the correlation parts in the EKT-MP2 results are recognized. Furthermore, the close relationship between the EKT-MP2 method and the derivative approach of the MP2 energy with respect to the orbital occupation numbers [N. Q. Su and X. Xu, J. Chem. Theory Comput. 2015, 11, 4,677] is revealed. When these MP2-based methods are applied to a set of atoms and molecules, new insights are gained on the role played by the relaxation and the correlation effects in the electron ionization processes.     

Near-degeneracy corrections for second-order perturbation theory: comparison of two approaches

 We compare two approximate perturbation schemes which were developed recently to deal with the (quasi)degeneracy problem in many-body perturbation theory. We conclude that although the two methods were introduced on quite different theoretical grounds, their performances are quite similar, and present an improvement over traditional perturbation theory. Both methods are cheap in computation time, but cannot compete in accuracy with more sophisticated schemes such as complete-active-space perturbation theory or dressed particle theories. Received: 1 August 2000 / Accepted: 2 August 2000 / Published online: 19 January 2001     

Dual-basis second-order Moller-Plesset perturbation theory: A reduced-cost reference for correlation calculations

The resolution-of-the-identity (RI) approximation has placed the onus of the cost of a second-order Moller-Plesset (MP2) calculation on the underlying self-consistent field (SCF) calculation for many moderately sized molecules. A dual-basis approach to the SCF calculation, based on previous methods demonstrated for density functional theory, is combined with RI-MP2 calculations, and small basis subsets for cc-pVTZ, cc-pVQZ, and 6-311++G(3df,3pd) are presented. These subsets provide time savings of greater than 90%, with negligible errors in absolute and relative energies, compared to the associated full-basis counterpart. The method is tested with a series of rotational barriers, relative conformational energies of alanine tetrapeptides, as well as the full G3/99 molecular set. RI-MP2 calculations on alanine octapeptides (40 heavy atoms, 3460 basis functions), using cc-pVQZ, are presented. Results improve upon previous methods that diagonalize the virtual space separately.     

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-Schro?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.     

Analytic geometrical derivatives of second-order molecular properties from perturbation theory

Assuming the Born-Oppenheimer approximation for molecular wavefunctions satisfies the Hellmann-Feynman theorem, Rayleigh-Schrödinger perturbation theory is employed to develop an analytic formula for derivatives of expectation values and second-order properties with respect to nuclear coordinates.     

One-electron properties of molecules calculated using second-order multireference perturbation theory

A new form of second-order multireference perturbation theory coupled with finite-field perturbation theory is applied to evaluate some one-electron molecular properties. Several possible definitions of the zeroth-order Hamiltonian are considered and results tested against bench-mark full CI calculations. © 1995 John Wiley & Sons, Inc.     

Dissociation energies within selected configuration interaction and perturbation theory

Selected configuration interaction (CI) calculations and second-order perturbational theory are used to truncate systematically multireference single and double excitation CI (MRCI) expansions in the calculation of the bond dissociation energies of several systems like the single-bonded LiF molecule or the multiple-bonded N₂, NO and O₂ diatomic systems. The method is extended to compute the CH bond dissociation energy ofethene C₂H₄. It is shown how the proposed scheme (perturbation-selected MRCI (MRCI-PS)) is able to reproduce the accuracy of complete MRCI expansions with only a small number of configurations variationally evaluated.     

Application of perturbation theory in large configuration interaction calculations

Rayleigh-Schrödinger perturbation theory has been applied through fifth order in the energy, to the problem of estimating the roots of the secular equation in large configuration interaction calculations. The NO₂⁺, O₃ and H₂O molecules are used as test cases, with accuracy as good as 0.01 eV, with appropriate choice of zero order problem.     

Time-dependent perturbation theory for vibrational energy relaxation and dephasing in peptides and proteins

Without invoking the Markov approximation, we derive formulas for vibrational energy relaxation (VER) and dephasing for an anharmonic system oscillator using a time-dependent perturbation theory. The system-bath Hamiltonian contains more than the third order coupling terms since we take a normal mode picture as a zeroth order approximation. When we invoke the Markov approximation, our theory reduces to the Maradudin-Fein formula which is used to describe the VER properties of glass and proteins. When the system anharmonicity and the renormalization effect due to the environment vanishes, our formulas reduce to those derived by and Mikami and Okazaki [J. Chem. Phys. 121, 10052 (2004)] invoking the path-integral influence functional method with the second order cumulant expansion. We apply our formulas to VER of the amide I mode of a small amino-acid like molecule, N-methylacetamide, in heavy water.     

Explicitly correlated second-order perturbation theory using density fitting and local approximations

Three major obstacles in electronic structure theory are the steep scalings of computer time with respect to system size and basis size and the slow convergence of correlation energies in orbital basis sets. Three solutions to these are, respectively, local methods, density fitting, and explicit correlation; in this work, we combine all three to produce a low-order scaling method that can achieve accurate MP2 energies for large systems. The errors introduced by the local approximations into the R12 treatment are analyzed for 16 chemical reactions involving 21 molecules. Weak pair approximations, as well as local resolution of the identity approximations, are tested for molecules with up to 49 atoms, over 100 correlated electrons, and over 1000 basis functions.     

An application of perturbation theory ideas in configuration interaction calculations

Configuration interaction calculations of electronic wave functions for atoms and molecules have generally been limited to relatively small basis sets because of the exponential increase in the number of configurations as basis functions are added. While higher than quadruply excited configurations are of negligible importance in CI wave functions, it is shown that the effect of triple and quadruple excitation configurations can be substantially included even when the matrix elements between such configurations are neglected, leaving only their diagonal elements and the elements connecting them with the single and double excitations. This approximation is seen to be formally practically equivalent to a first-order perturbation expression for the wave function (second-order for the energy) based on an optimum linear combination of the zero, single, and double excitation configurations as the zero-order function. If suitable procedures are used, the amount of computational effort involved in such a calculation is roughly proportional to the fourth power of the number of basis functions employed, thus preventing the CI stage of the calculation from increasing in magnitude much faster than the stages involving the calculation and manipulation of the elementary integrals.     

The shifted scheme in the general-model-space diagrammatic perturbation theory

The shifted scheme of many-body perturbation theory is applied to open-shell states within the framework of the general-model-space theory. Rules for shifting the denominators of folded diagrams. which appear in open-shell perturbation expansions, are given. The finite-order energies in the shifted scheme obtained in two equivalent representations may differ. This happens, for instance, in the case of triplet states. For ³Σ_u⁺ states of the He₂, differences up to 0.07 mhartree have been found in third order. A similar phenomenon is the size inconsistency of the shifted scheme observed by Silver in the ground state of He₂. A possible advantage of the shifted scheme is its faster convergence for excited states.     

Quasi-degenerate second-order perturbation theory for occupation restricted multiple active space self-consistent field reference functions

A multi-configuration quasi-degenerate second-order perturbation method based on the occupation restricted multiple active space (ORMAS-PT/ORMAS) reference wavefunction is presented. ORMAS gives one the ability to approximate a complete active space self-consistent field (CASSCF) wavefunction using only a subset of the configurations from the CASSCF space. The essential idea behind ORMAS-PT is to use the multi-reference M?ller-Plesset formalism to correct the ORMAS reference energy. A computational scheme employing direct CI methodology is presented. Several tests are presented to demonstrate the performance of the ORMAS-PT method.     

Configuration selection as a route towards efficient vibrational configuration interaction calculations

A configuration selective vibrational configuration interaction (CI) approach is presented that efficiently reduces the variational space and thus leads to significant speedups in comparison to standard vibrational CI implementations. Deviations with respect to reference calculations are well below the accuracy of the underlying electronic structure calculations for the potential and hence are essentially negligible. Parallel implementations of the presented configuration selective vibrational CI approaches lead to further significant time savings. Benchmark calculations based on potential energy surfaces of coupled-cluster quality are presented for the fundamental modes of cis- and trans-difluoroethylene. The size-consistency error within the vibrational configuration interaction calculations of the difluoroethylene dimer has been studied in dependence on the excitation level.