Operators Kohn-Sham

Properties can be computed by finding the expectation value of the property operator with the natural orbitals weighted by the occupation number of each orbital. This is a much faster way to compute properties than trying to use the expectation value of a multiple-determinant wave function. Natural orbitals are not equivalent to HF or Kohn-Sham orbitals, although the same symmetry properties are present. 

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. 

The term in square brackets defines the Kohn-Sham one-electron operator and equation (7-1) can be written more compactly as... 

This step is similar to what we have done in equation (7-7) where we obtained the matrix representation of the Kohn-Sham operator. If we insert expression (7-14) for the charge density in terms of the LCAO functions and make use of the density matrix P defined in equation (7-15), we arrive at... 

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... 

Various reasons have been advanced for the relative accuracy of spin-polarized Kohn-Sham calculations based on local (spin) density approximations for E c- However, two very favourable aspects of this procedure are clearly operative. First, the Kohn-Sham orbitals control the physical class of density functions which are allowed (in contrast, for example, to simpler Thomas-Fermi theories). Second, local density approximations for are mild-mannered,... 

Comparison of the Kohn-Sham and Skyrme functionals leads to a natural question why these two functionals exploit, for the time-dependent problem, so different sets of basic densities and currents If the Kohn-Sham functional is content with one density, the Skyrme forces operate with a diverse set of densities and currents, both T-even and T-odd. Then, should we consider T-odd densities as genuine for the description of dynamics of finite many-body systems or they are a pequliarity of nuclear forces This question is very nontrivial and still poorly studied. We present below some comments which, at least partly, clarify this point. 

All electron calculations were carried out with the DFT program suite Turbomole (152,153). The clusters were treated as open-shell systems in the unrestricted Kohn-Sham framework. For the calculations we used the Becke-Perdew exchange-correlation functional dubbed BP86 (154,155) and the hybrid B3LYP functional (156,157). For BP86 we invoked the resolution-of-the-iden-tity (RI) approximation as implemented in Turbomole. For all atoms included in our models we employed Ahlrichs valence triple-C TZVP basis set with polarization functions on all atoms (158). If not noted otherwise, initial guess orbitals were obtained by extended Hiickel theory. Local spin analyses were performed with our local Turbomole version, where either Lowdin (131) or Mulliken (132) pseudo-projection operators were employed. Broken-symmetry determinants were obtained with our restrained optimization tool (136). Pictures of molecular structures were created with Pymol (159). 

In the Kohn Sham equations (A.116) [324, 325], the core Hamiltonian operator h( 1) has the same definition as in HF theory (equation A.6), as does the Coulomb operator, 7(1), although the latter is usually expressed as... 

Note that because the E of Eq. (8.14) that we are minimizing is exact, the orbitals / must provide die exact density (i.e., the minimum must correspond to reality). Further note that it is these orbitals that form the Slatcr-dctcrminantal eigenfunction for the separable noninteracting Hamiltonian defined as the sum of the Kohn-Sham operators in Eq. (8.18), i.e.. 

It is important to emphasize that nearly all applications of DFT to molecular systems are undertaken within the context of the Kohn-Sham SCF approach. The motivation for this choice is that it permits the kinetic energy to be computed as the expectation value of the kinetic-energy operator over the KS single determinant, avoiding the tricky issue of... 

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... 

A DFT-based third order perturbation theory approach includes the FC term by FPT. Based on the perturbed nonrelativistic Kohn-Sham orbitals spin polarized by the FC operator, a sum over states treatment (SOS-DFPT) calculates the spin orbit corrections (35-37). This approach, in contrast to that of Nakatsuji et al., includes both electron correlation and local origins in the calculations of spin orbit effects on chemical shifts. In contrast to these approaches that employed the finite perturbation method the SO corrections to NMR properties can be calculated analytically from... 

Kohn-Sham orbitals (18)), Vn is the external, nuclear potential, and p is the electronic momentum operator. Hence, the first integral represents the kinetic and potential energy of a model system with the same density but without electron-electron interaction. The second term is the classical Coulomb interaction of the electron density with itself. Exc> the exchange-correlation (XC) energy, and ENR are functionals of the density. The exact functional form for Exc is unknown it is defined through equation 1 (79), and some suitable approximation has to be chosen in any practical application of... 

QR Method. The first relativistic method is the so-called quasi-relativistic (QR) method. It has been developed by Snijders, Ziegler and co-workers (13). In this approach, a Pauli Hamiltonian is included into the self-consistent solution of the Kohn-Sham equations of DFT. The Pauli operator is in a DFT framework given by... 

---