Hamiltonian operator free energy calculations

In Fenske and Hall s parameter-free SCF calculations (80-84), the He1t 1-electron operator is substituted by a model 1-electron operator that has a kinetic energy and potential energy term for each atomic center in the complex. This approach assumes that the electron density may be assigned to appropriate centers. The partitioning of electron density is done through Mulliken population analyses (163) until self-consistency is obtained. The Hamiltonian elements are evaluated numerically, and the energies of the MO s depend only on the choice of basis functions and the intemuclear distance. 

The Kohn-Sham-Dirac equation (27) introduced in the last section is the basis of most relativistic electronic structure calculations in solid state theory. There are certain aspects which make the numerical solution of this four-component equation more involved than its non-relativistic coimterpart The Hamiltonian of the Kohn-Sham-Dirac equation is, unlike its Schrodinger equivalent and unlike the field-theoretical Hamiltonian (7) with the properly chosen normal order, not bounded below. In the limit of free, non-interacting particles the solutions of the Kohn-Sham-Dirac equation are plane waves with energies e(k) = cVk -I- c, where positive energies correspond to electrons and states with negative energy can be interpreted as positrons. For numerical procedures, which preferably use variational techniques to find electronic solutions, this property of the Dirac operator causes a severe problem, which can be circumvented by certain techniques like the application of a squared Dirac operator or a projection onto the properly chosen electronic states according to their above definition after Eq. (19). 

All these basis sets are essentially optimized for the calculation of electronic energies and are therefore able to represent the operators included in the field-free electronic Hamiltonian reasonably well. However, in the calculation of molecular electromagnetic properties it is necessary also to represent other operators such as the electric dipole operator, the electronic angular momentum operator, the Fermi-contact operator and more. Most of these basis sets are a priori not optimized for this and have to be extended. 

---