Multireference scheme

The development of a general multiroot multireference scheme for treating electron correlation effects usually starts from consideration of the Schrodinger equation for a number (d) of target states. 

The methodology presented here expands the recent CASPT2 approach to more flexible zeroth-order variational spaces for a multireference perturbation, either in the Moller-Plesset scheme or in Epsein-Nesbet approach [70-72]. Furthermore, it allows for the use of a wide set of possible correlated orbitals. These two last points were discussed elsewhere [34]. 

Hamiltonian proposed by Muller and Plesset gives rise to a very successful and efficient method to treat electron correlation effects in systems that can be described by a single reference wave function. However, for a multireference wave function the Moller-Plesset division can no longer be made and an alternative choice of B(0> is needed. One such scheme is the Complete Active Space See-ond-Order Perturbation Theory (CASPT2) developed by Anderson and Roos [3, 4], We will briefly resume the most important definitions of the theory one is referred to the original articles for a more extensive description of the method. The reference wave function is a CASSCF wave function that accounts for the largest part of the non-dynamical electron correlation. The zeroth-order Hamiltonian is defined as follows and reduces to the Moller-Plesset operator in the limit of zero active orbitals ... 

An alternative contraction scheme which has received more attention is internally contracted multireference CISD (usually denoted simply CMRCI), which was first discussed by Meyer99 and Siegbahn.100 This method applies the single and double excitation operators to a single multiconfigurational reference wavefunction as a whole, including the reference coefficients. Thus, if the reference wavefunction is... 

Recent developments include exact [12-14, 44, 90, 91] and approximate [14, 90, 92-94] iterative schemes to determine Hg, the intermediate Hamiltonian method [21, 24, 95], the use of incomplete model spaces [43, 44] and some multireference open-shell coupled-cluster (CC) formalisms [16-20, 96, 97]. Only some eigenvalues of the intermediate Hamiltonian H, are also eigenvalues of H. The corresponding model eigenvectors of H, are related to their true counterparts as in Bloch s theory. Provided effective operators a are restricted to act solely between these model eigenvectors, the possible a definitions from Bloch s formalism (see Section VI.A) can be used. 

In many multireference perturbation theories, the energy spectrum computed with the approximated Hamiltonian is ill conditioned. Functions outside the reference space become artificially degenerate with reference wavefunction. The phenomenon is known in the literature as intruder state, and it results in unphysical energy corrections and spurious bumps along a potential energy surface. Several schemes, such as ad hoc level shift parameter, have been proposed [97-100], yet... 

---