6.17 Complete Active Space Methods

6.17.2 Theory

The electronic energy is an exact functional of the 1-RDM and 2-RDM

E=12pqrsΓpqrsgpqrs+pqDpqhpq, (6.48)

Given the 1- and 2-PDMs, the generalized Fock matrices may be generated for this MCSCF. The derivation and further details are neatly described by Helgaker, Jorgensen, and Olsen 365, but the key results are summarized here. In the following, m,n,p,q, are general indices, i,j,k, are inactive indices, t,u,v,w, are active indices, and a,b,c, are virtual indices.

The generalized Fock matrix is defined as:

Fmn=qDpqhpq+qrsΓmqrsgnqrs (6.49)

where hpq are the 1-electron integrals and gnqrs are the 2-electron integrals and all indices run over all orbital classes (inactive, active, and virtual). This, generally non-symmetric, matrix can be simplified by taking advantage of the fact that the form of the density matrices when some indices are inactive or virtual are much simpler than when the indices are active. When the first index of the generalized Fock matrix is inactive and the second is general:

Fin=2(FniI+FniA) (6.50)

where the inactive and active Fock matrices are

FmnI =hmn+i(2gmnii-gmiin) (6.51)
FmnA =vwDvw(gmnvw-gmwvn) (6.52)

In other words, the inactive Fock matrix is the Fock matrix formed from using only the inactive density and the active Fock matrix is sum of J and K matrices built from the active space 1-PDM. When the first index is active, and the second index is general, we have

Ftn=uFnuIDvu+Qtn (6.53)

where the auxiliary Q matrix is

Qtm=u,v,wΓtuvwgmuvw (6.54)

and finally, if the first index is virtual then Fan=0. This formulation of the generalized Fock matrix is quite useful because it only requires density matrices with all indices active and two-electron integrals in the MO basis with three indices active and one general index, greatly reducing the storage and computational cost of the MO transformation.

The orbital gradient is then given by

EΔpq=2(Fpq-Fqp) (6.55)