Many-body perturbation theory computational components

If effects from electron correlation on parity violating potentials shall be accounted for in a four-component framework, the situation becomes more complicated than in the Dirac Hartree-Fock case. This is related to the fact, that in four-component many body perturbation theory (MBPT) or in a four-component coupled cluster (CC) scheme the reduced density matrices on the respective computational level are required in order to determine the parity violating potentials. Since these densities were not available in analytic form, Thyssen, Laerdahl and Schwerdtfeger [153] used a finite field approach to compute parity violating potentials in a four-component framework on a correlated level. This amounts to adding the parity violating operator with different scaling factors A to the... 

In the earlier sections of this chapter we reviewed the many-electron formulation of the symmetry-adapted perturbation theory of two-body interactions. As we saw, all physically important contributions to the potential could be identified and computed separately. We follow the same program for the three-body forces and discuss a triple perturbation theory for interactions in trimers. We show how the pure three-body effects can be separated out and give working equations for the components in terms of molecular integrals and linear and quadratic response functions. These formulas have a clear, partly classical, partly quantum mechanical interpretation. The exchange terms are also classified for the explicit orbital formulas we refer to Ref. (302). 

---