CI wavefunction norm in Gaussian



    Dear fellow netters,
          I  am running CI(SD) calculations using G90, and  I
     am   trying  to  extract  the  coefficient  of  the   HF
     configuration. The manual says (page 165) that "Norm(A1)
     - 1 gives a measure of the correlation correction to the
     wavefunction." But what exactly is the  coefficient  for
     the  HF configuration? My guess is that it is  given  by
     1/Norm(A1), but I am not sure. Any comments?
          Thank you,
                    Gernot Frenking
                    frenking { *at * } dmrhrz11.bitnet
 Your guess is correct.  This is one of the most commonly asked questions,
 and the G92 manual has the equation for Norm(A) in an attempt to make the
 relationship clearer.  The wavefunction is stored in intermediate
 normalization.  That is:
 Psi(CISD) = Psi(0) + Sum(ia) Tia Psi(i-->a) + Sum (ijab) Tijab
 Psi(ij-->ab)
 Where Psi(0) is the Hartree-Fock determinant and has a coefficient of 1
 (that's what intermediate normalization means).  Norm(A) is the factor
 by which to divide the wavefunction as given above to fully normalize it.
 Thus
 Norm(A) = Sqrt (1 + Sum(ia) Tia*Tia + Sum(ijab) Tijab*Tijab)
 The coefficient of the Hartree-Fock determinant in the fully normalized
 wavefunction is then 1/Norm(A), the coefficient of singly excited
 determinant Psi(i-->a) is Tia/Norm(A), etc.
 Another pitfall to be careful of is that the sum over the doubly-excited
 determinants is really restricted to the unique determinants.  It is
 really
 Sum(i<j,a<b) Psi(ij-->ab ijab all alpha spin) +
 Sum(all ijab) Psi(ij-->ab i,a alpha j,b beta) +
 Sum(i<j,a<b) Psi(ij-->ab ijab all beta spin)
 For closed shell, the all beta coefficients are equal to the all alpha
 coefficients and if you request that the coefficients be printed (using
 IOP(9/28=N) to get coefficients larger than 10^-N or IOP(9/28=-1) to get all
 coefficients) only the alpha-alpha and alpha-beta coefficients will be
 printed.
 Mike Frisch
 -------
 ---
 Administrivia: This message is automatically appended by the mail exploder.
 CHEMISTRY { *at * } ccl.net --- everybody; CHEMISTRY-REQUEST { *at * } ccl.net
 --- coordinator only
 OSCPOST { *at * } ccl.net : send something from chemistry; FTP: www.ccl.net
 ---