Orbital energy functions

The reference state T is a Slater determinant constructed as a normalized antisymmetrized product of N orthonormal spin-indexed orbital functions orbital orthonormality constraint (i j) = StJ, imposed by introducing a matrix of Lagrange multipliers. The variational condition is... 

Since the exact ground-state electronic wave function and density can only be approximated for most A-electron systems, a variational theory is needed for the practical case exemplified by an orbital functional theory. As shown in Section 5.1, any rule T 4> defines an orbital functional theory that in principle is exact for ground states. The reference state for any A-electron wave function T determines an orbital energy functional E = Eq + Ec,in which E0 = T + Eh + Ex + V is a sum of explicit orbital functionals, and If is aresidual correlation energy functional. In practice, the combination of exchange and correlation energy is approximated by an orbital functional Exc. 

In each SCF iteration, the Hartree-Fock orbital-energy function (10.7.98) may be minimized using any of the standard iterative methods of uneonstrained optimization. Moreover, since the optimization of Eg(X) requires no re-evaluation of the Fock matrix, each step is relatively inexpensive. 

SmarglassI E and Madden P A 1994 Orbital-free kinetic-energy functionals for first-principles molecular dynamics Phys. Rev. B 49 5220-6... 

Wang Y A, Govind N and Carter E A 1998 Orbital-free kinetic energy functionals for the nearly-free electron gas Phys. Rev. B 58 13 465... 

VV e now wish to establish the general functional form of possible wavefunctions for the two electrons in this pseudo helium atom. We will do so by considering first the spatial part of the u a efunction. We will show how to derive functional forms for the wavefunction in which the i change of electrons is independent of the electron labels and does not affect the electron density. The simplest approach is to assume that each wavefunction for the helium atom is the product of the individual one-electron solutions. As we have just seen, this implies that the total energy is equal to the sum of the one-electron orbital energies, which is not correct as ii ignores electron-electron repulsion. Nevertheless, it is a useful illustrative model. The wavefunction of the lowest energy state then has each of the two electrons in a Is orbital ... 

In Ecjuation (3.47) we have written the external potential in the form appropriate to the interaction with M nuclei. , are the orbital energies and Vxc is known as the exchange-correlation functional, related to the exchange-correlation energy by ... 

One of the advantages of this method is that it breaks the many-electron Schrodinger equation into many simpler one-electron equations. Each one-electron equation is solved to yield a single-electron wave function, called an orbital, and an energy, called an orbital energy. The orbital describes the behavior of an electron in the net field of all the other electrons. 

In this formulation, the electron density is expressed as a linear combination of basis functions similar in mathematical form to HF orbitals. A determinant is then formed from these functions, called Kohn-Sham orbitals. It is the electron density from this determinant of orbitals that is used to compute the energy. This procedure is necessary because Fermion systems can only have electron densities that arise from an antisymmetric wave function. There has been some debate over the interpretation of Kohn-Sham orbitals. It is certain that they are not mathematically equivalent to either HF orbitals or natural orbitals from correlated calculations. However, Kohn-Sham orbitals do describe the behavior of electrons in a molecule, just as the other orbitals mentioned do. DFT orbital eigenvalues do not match the energies obtained from photoelectron spectroscopy experiments as well as HF orbital energies do. The questions still being debated are how to assign similarities and how to physically interpret the differences. 

Semiempirical programs often use the half-electron approximation for radical calculations. The half-electron method is a mathematical technique for treating a singly occupied orbital in an RHF calculation. This results in consistent total energy at the expense of having an approximate wave function and orbital energies. Since a single-determinant calculation is used, there is no spin contamination. 

Once you have calculated an ab initio or a semi-empirical wave function via a single point calculation, geometry optimization, molecular dynamics or vibrations, you can plot the electrostatic potential surrounding the molecule, the total electronic density, the spin density, one or more molecular orbitals /i, and the electron densities of individual orbitals You can examine orbital energies and select orbitals for plotting from an orbital energy level diagram. 

Multiple solutions /j and 8j are possible for this last equation. The wave functions for individual electrons, /j, are called molecular orbitals, and the energy, 8j, of an electron in orbital /j is called the orbital energy. 

For a quantum mechanical calculation, the single point calculation leads to a wave function for the molecular system and considerably more information than just the energy and gradient are available. In principle, any expectation value might be computed. You can get plots of the individual orbitals, the total (or spin) electron density and the electrostatic field around the molecule. You can see the orbital energies in the status line when you plot an orbital. Finally, the log file contains additional information including the dipole moment of the molecule. The level of detail may be controlled by the PrintLevel entry in the chem.ini file. 

Here, c is a column vector of LCAO coefficients and e is called the orbital energy. If we start with n basis functions, then there are exactly n different c s (and e s) and the m lowest-energy solutions of the eigenvalue problem correspond to the doubly occupied HF orbitals. The remaining n — m solutions are called the virtual orbitals. They are unoccupied. 

So far there have not been any restrictions on the MOs used to build the determinantal trial wave function. The Slater determinant has been written in terms of spinorbitals, eq. (3.20), being products of a spatial orbital times a spin function (a or /3). If there are no restrictions on the form of the spatial orbitals, the trial function is an Unrestricted Hartree-Fock (UHF) wave function. The term Different Orbitals for Different Spins (DODS) is also sometimes used. If the interest is in systems with an even number of electrons and a singlet type of wave function (a closed shell system), the restriction that each spatial orbital should have two electrons, one with a and one with /3 spin, is normally made. Such wave functions are known as Restricted Hartree-Fock (RHF). Open-shell systems may also be described by restricted type wave functions, where the spatial part of the doubly occupied orbitals is forced to be the same this is known as Restricted Open-shell Hartree-Fock (ROHF). For open-shell species a UHF treatment leads to well-defined orbital energies, which may be interpreted as ionization potentials. Section 3.4. For an ROHF wave function it is not possible to chose a unitary transformation which makes the matrix of Lagrange multipliers in eq. (3.40) diagonal, and orbital energies from an ROHF wave function are consequently not uniquely defined, and cannot be equated to ionization potentials by a Koopman type argument. 

The foundation for the use of DFT methods in computational chemistry was the introduction of orbitals by Kohn and Sham. 5 The main problem in Thomas-Fermi models is that the kinetic energy is represented poorly. The basic idea in the Kohn and Sham (KS) formalism is splitting the kinetic energy functional into two parts, one of which can be calculated exactly, and a small correction term. 

Kolos, W., J. Chem. Phys. 27, 591, Excitation energies of C2H4. The correlation between electrons with opposite spins is estimated by multiplying the usual orbital wave functions by the inter-electronic distance. 

The goodness of the PP representation can be checked by comparing the all-electron and PP orbital energies and relative stability of atomic states. The comparison is shown in Table 4, and is seen to be very satisfying. For a balanced treatment, also the carbon and oxygen atoms were treated by a PP, as described in previous work5.3d functions were not introduced in the sulphur basis set, mainly because they were not deemed necessary for the illustrative purposes of this chapter. Also, the derivation of a PP representation for polarization functions is not a straightforward matter. The next section is devoted to the discussion of this point. 

Figure 6.14a shows the sp and d bands of a transition metal (e.g. Pt), i.e. the density of states (DOS) as a function of electron energy E. It also shows the outer orbital energy levels of a gaseous CO molecule. Orbitals 4a, l7t and 5cr are occupied, as indicated by the arrows, orbital 27c is empty. The geometry of these molecular orbitals is shown in Figure 6.14b. 

---