Kohn-Sham-Fock operator

The FH method assumes that the Fock/Kohn Sham (KS) operator can be written as... 

In (9.6) F k) is the matrix of the Hartree-Fock (HP) or Kohn-Sham (KS) operator. The former operator includes a nonlocal exchange part, depending on the density matrix p R,R ), whereas the latter operator involves the electron density p R) = p R, R), that is, it depends only on the diagonal elements of the density matrix, see Chapters 4 and 7. 

Since the Fock operator is a effective one-electron operator, equation (1-29) describes a system of N electrons which do not interact among themselves but experience an effective potential VHF. In other words, the Slater determinant is the exact wave function of N noninteracting particles moving in the field of the effective potential VHF.5 It will not take long before we will meet again the idea of non-interacting systems in the discussion of the Kohn-Sham approach to density functional theory. 

We may ask now, whether the same procedure may be applied to density-functional theory, just by replacing the Fock operator by the corresponding Kohn-Sham operator. To this end we have to look at the minimization of the total energy with respect to the density of a multi-determinantal wavefunction 4. We write the density as ... 

Note that the Kohn-Sham Hamiltonian hKS [Eq. (4.1)] is a local operator, uniquely determined by electron density15. This is the main difference with respect to the Hartree-Fock equations which contain a nonlocal operator, namely the exchange part of the potential operator. In addition, the KS equations incorporate the correlation effects through Vxc whereas they are lacking in the Hartree-Fock SCF scheme. Nevertheless, though the latter model cannot be considered a special case of the KS equations, there are some similarities between the Hartree-Fock and the Kohn-Sham methods, as both lead to a set of one-electron equations allowing to describe an n-electron system. 

Natural steric analysis57 allows quantitative evaluation of steric repulsion on the basis of this simple physical picture. Given the converged Fock (or Kohn-Sham) operator F, we can evaluate the average energy of each occupied NBO f2/NI 0) and the associated pre-orthogonal PNBO C/PNIi0j in the usual manner,... 

How do CMOs and LMOs differ The CMOs are symmetry-adapted eigenfunctions of the Fock (or Kohn-Sham) operator F, necessarily reflecting all the molecular point-group symmetries of F itself,26 whereas the LMOs often lack... 

NBO analysis can be used to quantify this phenomenon. Since tire NBOs do not diagonalize the Fock operator (or tire Kohn-Sham operator, if the analysis is carried out for DFT instead of HF), when the Fock matrix is formed in the NBO basis, off-diagonal elements will in general be non-zero. Second-order perturbation tlieory indicates that these off-diagonal elements between filled and empty NBOs can be interpreted as the stabilization energies... 

Use the K-S operator hKS and the basis functions (f>) to calculate Kohn-Sham matrix elements hrs (cf. Fock matrix elements Frs (Section 5.23.6),... 

In the above formula, HFvx([ J r) is the familiar nonlocal Fock exchange operator, but here built from the one-particle Kohn-Sham orbitals of instead of from the Hartree-Fock orbitals. 

The first derivative of the density matrix with respect to the magnetic induction (dPfiv/dBi) is obtained by solving the coupled-perturbed Hartree-Fock (or Kohn-Sham) equations to which the first derivative of the effective Fock (or Kohn-Sham) operator with respect to the magnetic induction contributes. Due to the use of GIAOs, specific corrections arising from the effective operator Hcnv describing the environment effects will appear. We refer to Ref. [28] for the PCM model and to Ref. [29] for the DPM within either a HF or DFT description of the solute molecule. 

A similar expression is obtained for the Kohn-Sham operator in density functional. One should note that because of the quadratic nature of these SCRF equations on the molecular orbital coefficients, the V2 factor in the total energy is not present in the Fock operator. 

For variational methods, such as Hartree-Fock (HF), multi-configurational self-consistent field (MCSCF), and Kohn-Sham density functional theory (KS-DFT), the initial values of the parameters are equal to zero and 0) thus corresponds to the reference state in the absence of the perturbation. The A operators are the non-redundant state-transfer or orbital-transfer operators, and carries no time-dependence (the sole time-dependence lies in the complex A parameters). Furthermore, the operator A (t)A is anti-Hermitian, and tlie exponential operator is thus explicitly unitary so that the norm of the reference state is preserved. Perturbation theory is invoked in order to solve for the time-dependence of the parameters, and we expand the parameters in orders of the perturbation... 

Because the convenience of the one-electron formalism is retained, DFT methods can easily take into account the scalar relativistic effects and spin-orbit effects, via either perturbation or variational methods. The retention of the one-electron picture provides a convenient means of analyzing the effects of relativity on specific orbitals of a molecule. Spin-unrestricted Hartree-Fock (UHF) calculations usually suffer from spin contamination, particularly in systems that have low-lying excited states (such as metal-containing systems). By contrast, in spin-unrestricted Kohn-Sham (UKS) DFT calculations the spin-contamination problem is generally less significant for many open-shell systems (39). For example, for transition metal methyl complexes, the deviation of the calculated UKS expectation values S (S = spin angular momentum operator) from the contamination-free theoretical values are all less than 5% (32). 

---