next up previous contents index
Next: Full CI Algorithm Up: Restricted Active Space CI Previous: Restricted Active Space CI

Olsen's Full CI tex2html_wrap_inline4376 Equations

  We now turn our attention to the general formulation of the full CI problem in terms of alpha and beta strings. Later on, we will consider how to modify our results for RAS CI's. We begin by describing Olsen's expressions for tex2html_wrap_inline4376, which is the action of the Hamiltonian tex2html_wrap_inline3269 on the CI vector (or matrix) tex2html_wrap_inline4274. First we express tex2html_wrap_inline3269 in second-quantized  form (see section 5).
 equation1857
where tex2html_wrap_inline4386 is the shift operator
 equation1879
Again, tex2html_wrap_inline4376 is given by
 equation1890
As we will proceed to show, tex2html_wrap_inline4376 can be split up into three terms: one involving only beta components of the linear group generators (tex2html_wrap_inline4392), one involving only alpha components of the generators (tex2html_wrap_inline4394), and one involving mixtures of the two (tex2html_wrap_inline4396). Inserting equation (6.16) into equation (6.18) yields


 eqnarray1906

Now expanding the shift operators according to equation (6.17), we write tex2html_wrap_inline4376 as a sum of three terms
equation1936
where
eqnarray1946
and
eqnarray1978
and
equation2010

Obviously, these expressions can be simplified further. First, we observe that the expression for tex2html_wrap_inline4392 contains no tex2html_wrap_inline4116 operators. Therefore, tex2html_wrap_inline4404 unless tex2html_wrap_inline4406. Using this fact, and integrating out the tex2html_wrap_inline4116 part, we obtain
equation2042

Taking out the Kronecker delta term and rearranging, we have
eqnarray2071
Now the second term can be combined with the first to give equation (9a) of reference [17].
 eqnarray2114

Similarly, tex2html_wrap_inline4394 can be simplified to
 eqnarray2148

For efficient implementation, it is convenient to precompute the quantities
equation2182
Finally, we simplify tex2html_wrap_inline4396. It may be rewritten as
eqnarray2189
Since we sum over all ijkl, we can permute i and j with k and l. We can also swap tex2html_wrap_inline4424 and tex2html_wrap_inline4426 as can easily be verified. This yields equation (9c) from reference [17].
eqnarray2234

Thus we have written the action of the Hamiltonian on the current CI vector in terms of alpha and beta strings and alpha and beta shift operators. The product tex2html_wrap_inline4376 is written as a sum of three terms: the first (tex2html_wrap_inline4392) involves only beta shift operators, the second (tex2html_wrap_inline4394) involves only alpha shift operators, and the third (tex2html_wrap_inline4396) involves both alpha and beta shift operators. Note that, except for the factor tex2html_wrap_inline4436, tex2html_wrap_inline4392 is independent of tex2html_wrap_inline4173 so that the algorithm for computing tex2html_wrap_inline4392 is vectorizable (see below). Analogous results hold for tex2html_wrap_inline4394. This situation does not obtain in the computation of tex2html_wrap_inline4396, however, which is the rate-limiting step.


next up previous contents index
Next: Full CI Algorithm Up: Restricted Active Space CI Previous: Restricted Active Space CI

CCQC WWW repository
Thu Aug 14 12:57:45 EDT 1997