A derivation of the frozen-orbital unrestricted open-shell and restricted closed-shell second-order perturbation theory analytic gradient expressions

A detailed derivation of the frozen-orbital second-order perturbation theory (MP2) analytic gradient in the spin-orbital basis is presented. The summation ranges and modification of the MP2 gradient terms that result from the frozen-orbital approximation are clearly identified. The frozen-orbital analytic gradients for unrestricted MP2 and closed-shell MP2 are determined from the spin-orbital derivation. A discussion of useful implementation procedures is included. Timings from full and frozen-orbital MP2 gradient calculations on the molecule silicocene (the silicon analog of the sandwich compound ferrocene) are also presented.     

Orbital-invariant formulation and second-order gradient evaluation in Møller-Plesset perturbation theory

Based on the Hylleraas functional form, the second and third orders of Møller-Plesset perturbation theory are reformulated in terms of arbitrary (e.g., localized) internal orbitals, and atomic orbitals in the virtual space. The results are strictly equivalent to the canonical formulation if no further approximations are introduced. The new formalism permits the extension of the local correlation method to Møller-Plesset theory. It also facilitates the treatment of weak pairs at a lower (e.g., second order) level of theory in CI and coupled cluster methods. Based on our formalism, an MP2 gradient algorithm is outlined which does not require the storage of derivative integrals, integrals with three external MO indices, and, using the method of Handy and Schaefer, the repeated solution of the coupled-perturbed SCF equations.     

A constant denominator perturbation theory for molecular energy

A constant denominator perturbation theory is developed based on a zeroth order Hamiltonian characterized by degenerate subsets of orbitals. Such a formulation allows for a decoupling of the numerators of the perturbation sequence, allowing for much more rapid evaluation of the resultant sums. For example, the full fourth order theory can be evaluated as an N ⁶ step rather than N ⁷, where N is proportional to the basis set.Although the theory is general, a constant denominator is chosen for this study as the difference between the average occupied and average virtual orbital energies scaled so that the first order wavefunction yields the lowest possible variational bound. The third order correction then appears naturally as a scaled Langhoff-Davidson correction. The full fourth order with this partitioning is developed. Results are presented within the localized bond model utilizing both the Pariser-Parr-Pople and CNDO/2 model Hamiltonians. The second order theory presents a useful bound, usually containing a good deal of the basis set correlation. In all cases examined the fourth order theory shows remarkable stability, even in those cases in which the Nesbet-Epstein partitioning seems poorly convergent, and the Moller-Plesset theory uncertain.     

Comparison of CEPA and CP-MET methods

The known CEPA variants CEPA (v) withv = 0,1,2,3 and two new ones withv = 4, 5 are compared both formally and for various numerical examples with CP-MET. The main conclusions are: 1. In those situations where both CP-MET and the CEPA variants are justified (i.e. for “good” closed shell states) the correlation energies obtained with the 7 different schemes differ very little (by something like ±2%), with CEPA (1) closest to CP-MET (difference usually a fraction of 1%) and CEPA (4) nearly as close; this is rather insensitive to whether one uses canonical or localized orbitals. Even CEPA (3) is not too far from CP-MET, which confirms an earlier suggestion of Kelly. 2. In those cases where one of the 7 schemes fails (e.g. due to near degeneracy as in covalent molecules at large internuclear distances) the other 6 usually fail as well, though CEPA (0) is then somewhat poorer than the other schemes. Then no longer CEPA (1) but rather CEPA (3) is closest to CP-MET and then all schemes converge much better in a localized representation. 3. CEPA (2) usually leads to best agreement with experiment since it simulates to some extent triple substitutions. In none of the studied examples does CP-MET show a significant superiority as compared to the other schemes. Possible improvements to extend the domain of applicability of these methods are discussed.     

On the Hellmann–Feynman theorem and the corrections to the energy in the Rayleigh–Schrödinger perturbation theory

In this work we present an alternative method, based on the Hellmann–Feynman theorem, to generate energy corrections within the standard Rayleigh–Schrödinger perturbation theory. As a result, compact expressions for the corrections to the energy at different orders are obtained. We also review a method that allows us to calculate the corrections to the wave function needed in the energy calculations. Finally, our results are compared with those obtained by other authors by a different technique.     

Many-body perturbation theory for open-shell systems. Expansion through fourth-order

All of the diagrams which arise in the many-body perturbation theory of open-shell systems using a restricted Hartree-Fock reference function are given through fourth-order in the energy. New effects which arise in fourthorder are discussed.S.E.R.C. Advanced Fellow     

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     

Multi‐state multi‐reference Møller–Plesset second‐order perturbation theory for molecular calculations

This work presents multi‐state multi‐reference Møller–Plesset second‐order perturbation theory as a variant of multi‐reference perturbation theory to treat electron correlation in molecules. An effective Hamiltonian is constructed from the first‐order wave operator to treat several strongly interacting electronic states simultaneously. The wave operator is obtained by solving the generalized Bloch equation within the first‐order interaction space using a multi‐partitioning of the Hamiltonian based on multi‐reference Møller–Plesset second‐order perturbation theory. The corresponding zeroth‐order Hamiltonians are nondiagonal. To reduce the computational effort that arises from the nondiagonal generalized Fock operator, a selection procedure is used that divides the configurations of the first‐order interaction space into two sets based on the strength of the interaction with the reference space. In the weaker interacting set, only the projected diagonal part of the zeroth‐order Hamiltonian is taken into account. The justification of the approach is demonstrated in two examples: the mixing of valence Rydberg states in ethylene, and the avoided crossing of neutral and ionic potential curves in LiF. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

A practical and efficient method to calculate AIM localization and delocalization indices at post-HF levels of theory

A practical and efficient method is proposed for calculating localization and delocalization indices at post-Hartree-Fock levels, and the method is tested at the CISD/6-311G++(2d, 2p) level for a large set of molecules. Our method, which utilizes wave functions written in the natural molecular orbital format and obtained with GAUSSIAN 94 or GAUSSIAN 98, convincingly extends the concepts established at the HF level.     

An overview of coupled cluster theory and its applications in physics

Summary What has since become known as the normal coupled cluster method (NCCM) was invented about thirty years ago to calculate ground-state energies of closed-shell atomic nuclei. Coupled cluster (CC) techniques have since been developed to calculate excited states, energies of open-shell systems, density matrices and hence other properties, sum rules, and the sub-sum-rules that follow from imbedding linear response theory within the NCCM. Further extensions deal both with systems at nonzero temperature and with general dynamical behaviour. More recently, a new version of CC theory, the so-called extended coupled cluster method (ECCM) has been introduced. It has the potential to describe such global phenomena as phase transitions, spontaneous symmetry breaking, states of topological excitation, and nonequilibrium behaviour. CC techniques are now widely recognized as providing one of the most universally applicable, most powerful, and most accurate of all microscopicab initio methods in quantum many-body theory. The number of successful applications within physics is now impressively large. In most such cases the numerical results are either the best or among the best available. A typical case is the electron gas, where the CC results for the correlation energy agree over the entire metallic density range to within less than 1 millihartree (or <1%) with the essentially exact Green's function Monte Carlo results. The role of CC theory within modern quantum many-body theory is first surveyed, by a comparison with other techniques. Its full range of applications in physics is then reviewed. These include problems in nuclear physics, both for finite nuclei and infinite nuclear matter; the electron gas; various integrable and nonintegrable models; various relativistic quantum field theories; and quantum spin chain and lattice models. Particular applications of the ECCM include the quantum hydrodynamics of a zero-temperature, strongly-interacting condensed Bose fluid; a charged impurity in a polarizable medium (e.g., positron annihilation in metals); and various anharmonic oscillator and spin systems.     

Comparison of the Hartree-Fock,Møller-Plesset,and Hartree-Fock-Slater method with respect to electrostatic properties of small molecules

Summary The results of various quantum chemical calculations, the Hartree-Fock (HF) method, the Møller-Plesset perturbation theory (MP2), and the Hartree-Fock-Slater (HFS) method are compared. Atomic charges, dipole moments, topological properties of the electron density distribution and polarizabilities, and first hyperpolarizabilities are calculated. Atomic charges obtained with the HFS method are found to be very close to those calculated with the MP2 method, from which we conclude that the HFS method describes to some extent electron correlation effects. Performing an MP2 calculation after an HF calculation improves the molecular dipole moments considerably, yielding values close to the experimental ones. HFS calculations are computationally less demanding than MP2 and yield comparable results for the electron density distributions, dipole moments and polarizabilities.