Single-reference function perturbation

Hubac and his co-workers222"231 have explored the use of Brillouin-Wigner perturbation theory in solving the coupled cluster equations. For the case of a single reference function, this approach is entirely equivalent to other formulations of the coupled cluster equations. However, for the multireference case, the Brillouin-Wigner coupled cluster theory shows some promise in that it appears to alleviate the intruder state problem. No doubt perturbative analysis will help to gain a deeper understanding of this approach. 

We turn now to the Brillouin-Wigner perturbation theory for a system described in zero order by a multi-reference function. The multi-reference formalism closely parallels that given in the previous section for the case of a single-reference function. Let us begin by defining a reference space V. Let... 

Brillouin-Wigner perturbation theory can be developed for both the single reference function case and the multireference function case using a common formalism. This contrasts with the situation for Rayleigh-Schrodinger perturbation theory. We shall, therefore, consider the single reference and multireference formalisms together. 

We have included the parameter A in eq. (1.13) which is set equal to unity in order to recover the perturbed problem. Equation (1.13) is the basic formula of the Brillouin-Wigner perturbation theory for a single-reference function. 

We are now in a position to obtain perturbation expansions by expanding the inverse operator in the effective Hamiltonian, the wave operator and the reaction operator. We begin, as we did in our discussion of the partitioning technique, by considering the case of a single-reference function and then turn our attention to the multi-reference function case. 

In the M0ller-Plesset formalism, a single-reference function is employed and the partition of the Hamiltonian into a reference or zero-order operator and a perturbation uses the Hartree-Fock model to define the reference. Third-order theory (mp3) and fourth-order theory (mp4) are computationally tractable. 

Abstract The Brillouin-Wigner many-body problem in atomic and molecular physics and in quantum chemistry is described. The use of coupled cluster expansions, configuration interaction and perturbation series is considered both for the single-reference function case and for those cases requiring the use of a multi-reference formalism. 

By following procedures similar to those employed in the case of a single-reference function, the exact wave function f o, for a = 1,2,. .., d, in the Brillouin-Wigner perturbation theory can be written as the expansion... 

Just as single reference Cl can be extended to MRCI, it is also possible to use perturbation methods with a multi-detenninant reference wave function. Formulating MR-MBPT methods, however, is not straightforward. The main problem here is similar to that of ROMP methods, the choice of the unperturbed Hamilton operator. Several different choices are possible, which will give different answers when the tlieory is carried out only to low order. Nevertheless, there are now several different implementations of MP2 type expansions based on a CASSCF reference, denoted CASMP2 or CASPT2. Experience of their performance is still somewhat limited. 

If we except the Density Functional Theory and Coupled Clusters treatments (see, for example, reference [1] and references therein), the Configuration Interaction (Cl) and the Many-Body-Perturbation-Theory (MBPT) [2] approaches are the most widely-used methods to deal with the correlation problem in computational chemistry. The MBPT approach based on an HF-SCF (Hartree-Fock Self-Consistent Field) single reference taking RHF (Restricted Hartree-Fock) [3] or UHF (Unrestricted Hartree-Fock ) orbitals [4-6] has been particularly developed, at various order of perturbation n, leading to the widespread MPw or UMPw treatments when a Moller-Plesset (MP) partition of the electronic Hamiltonian is considered [7]. The implementation of such methods in various codes and the large distribution of some of them as black boxes make the MPn theories a common way for the non-specialist to tentatively include, with more or less relevancy, correlation effects in the calculations. 

The tautomerism of furoxan (l,2,5-oxadiazole-2-oxide) has been investigated by different computational methods comprising modern density functions as well as single-reference and multi-reference ab initio methods. The ring-opening process to 1,2-dinitrosoethylene is the most critical step of the reaction and cannot be treated reliably by low-level computations (Scheme 2). The existence of cis-cis-trans- 1,2-dinitrosoethylene as a stable intermediate is advocated by perturbational methods, but high-level coupled-cluster calculations identify this as an artifact <2001JA7326>. 

The second step of the calculation involves the treatment of dynamic correlation effects, which can be approached by many-body perturbation theory (62) or configuration interaction (63). Multireference coupled-cluster techniques have been developed (64—66) but they are computationally far more demanding and still not established as standard methods. At this point, we will only focus on configuration interaction approaches. What is done in these approaches is to regard the entire zeroth-order wavefunc-tion Tj) or its constituent parts double excitations relative to these reference functions. This produces a set of excited CSFs ( Q) that are used as expansion space for the configuration interaction (Cl) procedure. The resulting wavefunction may be written as... 

CCSD is similarly sensitive to multireference character, aldiough it is less obvious that this should be so based on the formalism presented above. However, inclusion of triples in the CCSD wave function is usually very effective in correcting for a single-reference treatment of a weakly to moderately multireference problem. Of course, die most common way to include die triples is by perturbation theory, i.e., CCSD(T), and as noted above, this level too can be unstable if singles amplitudes are large. In such an instance, BD(T) calculations, which eliminate die singles amplitudes, can be efficacious. 

---