Fock-space coupled cluster method equations

The relativistic coupled cluster method starts from the four-component solutions of the Drrac-Fock or Dirac-Fock-Breit equations, and correlates them by the coupled-cluster approach. The Fock-space coupled-cluster method yields atomic transition energies in good agreement (usually better than 0.1 eV) with known experimental values. This is demonstrated here by the electron affinities of group-13 atoms. Properties of superheavy atoms which are not known experimentally can be predicted. Here we show that the rare gas eka-radon (element 118) will have a positive electron affinity. One-, two-, and four-components methods are described and applied to several states of CdH and its ions. Methods for calculating properties other than energy are discussed, and the electric field gradients of Cl, Br, and I, required to extract nuclear quadrupoles from experimental data, are calculated. 

Equations for the Fock space coupled cluster method, including all single, double, and triple excitations (FSCCSDT) for ionization potentials [(0,1) sector], are presented in both operator and spin orbital form. Two approximations to the full FSCCSDT equations are described, one being the simplest perturbative inclusion of triple excitation effects, FSCCSD+T(3), and a second that indirectly incorporates certain higher-order effects, FSCCSD+T (3). 

The basic relativistic equations are described in Sec. 2, and the Fock-space coupled cluster method is discussed in Sec. 3. The recently developed intermediate Hamiltonian approach is described and illustrated by several... 

L. Meissner and R. J. Bartlett, J. Chem. Phys., 94, 6670 (1991). Transformation of the Hamiltonian in Excitation Energy Calculations Comparison Between Fock-Space Multireference Coupled-Cluster and Equation-of-Motion Coupled-Cluster Methods. 

Kiylov, A I. (2008). Equation-of-Motion Coupled-Cluster Methods for Open-Shell and Electronically Excited Species The fJitchhiker s Guide to Fock Space, Annu. Rev. Phys. Chem., 59,433-462. 

In this section we will introduce some wavefunction-based methods to calculate photoabsorption spectra. The Hartree-Fock method itself is a wavefunction-based approach to solve the static Schrodinger equation. For excited states one has to account for time-dependent phenomena as in the density-based approaches. Therefore, we will start with a short review of time-dependent Hartree-Fock. Several more advanced methods are available as well, e.g. configuration interaction (Cl), multireference configuration interaction (MRCI), multireference Moller-Plesset (MRMP), or complete active space self-consistent field (CASSCF), to name only a few. Also flavours of the coupled-cluster approach (equations-of-motion CC and linear-response CQ are used to calculate excited states. However, all these methods are applicable only to fairly small molecules due to their high computational costs. These approaches are therefore discussed only in a more phenomenological way here, and many post-Hartree-Fock methods are explicitly not included. 

One of the original approximate methods is the wavefunction-theory-based Hartree-Fock (HF) method [40]. The HF method is a single determinant method that does not include any correlation interactions between the electrons, and as such has limited accuracy [41, 42]. Higher level wavefunction-based methods such as coupled cluster [43 5], configuration interaction [40,46,47], and complete active space [48-50] methods include multiple determinants to incorporate some of the electron-electron correlation. Methods based on perturbation theory, such as second order Mpller-Plesset perturbation theory [51], go beyond the HF method by perturbatively adding electron correlation. These correlated wavefunction-based methods have well-defined ways in which they approach the exact solution to the Schrodinger equation and thus have the potential to be extremely accurate, but this accuracy comes at a price [52]. 

Applications to atoms are in most cases based on the publicly available programs using finite difference methods for integration in the solution of the (multi-configurational) Dirac-Hartree-Fock equations. The problem of introducing electron correlation in this framework is most successfully accomplished by employing complete active space (CAS) and restricted active space (RAS) techniques (see Ref. 84 for a recent application with further references to the literature) or coupled-cluster techniques. ... 

The preceding step to both MP2 and coupled-cluster calculations is to solve the Hartree-Fock equations. The standard approach is, of course, to solve the equations in a basis set expansion (Roothaan-Hall method), using atom-centered basis functions. This set of basis functions is used to expand the molecular orbitals and we will call it orbital basis set (OBS). It spans the computational (finite) orbital space. Occupied spin orbitals will be denoted (pi and virtual (unoccupied) spin orbitals pa- In order to address the terms that miss in a finite OBS expansion, the set of virtual spin orbitals in a formally complete space is introduced, pa- If we exclude from this space all those orbitals which can be represented by the OBS, we obtain the complementary space, with orbitals denoted cp i. The subdivision of the orbital space and the index conventions are summarized in the left part of Fig. 2. 

---