Orbitals analytic

The electron—nucleus potential is not quite trivial. For larger atoms the radial functions that are large at the nucleus are affected by the finite charge distribution of the nucleus. It is sufficient to use the potential for a uniform charge distribution of radius R 

The extension of the matrix solution of section 4.3 for one-electron bound states to the Hartree—Fock problem has many advantages. It results in radial orbitals specified as linear combinations of analytic functions, usually normalised Slater-type orbitals (4.38). This is a very convenient form for the computation of potential matrix elements in reaction theory. The method has been described by Roothaan (1960) for a closed-shell or single-open-shell structure. 

We illustrate the method by applying it to the simplest closed-shell form (5.29) of the Hartree—Fock problem. The radial orbitals u (r) are expressed as a linear combination of basis radial orbitals /i/(r). 

We define a matrix notation, omitting explicit reference to the coordinate 

In order to take into account the two-electron interaction terms we define a supermatrix notation that we illustrate for The components of a supervector t are the matrix elements t ij in dictionary order. Supervectors are transformed by supermatrices. We need the supermatrix V, defined by 

---

Our main reason for introducing the Slater atomic orbitals, however, is that linear combinations of them have often been used to approximate the SCF Hartree-Fock numerical orbitals. If x represents a set of analytical orbitals, such as the STOs of equation (6.56), then we may expand 0 in the Slater determinant (6.36) in terms of x,... 

It is of course possible to solve the Dirac—Fock problem with a linear combination of analytic orbitals. However, owing to the rapid variation of the orbitals near the nucleus it requires an awkwardly-large basis. If an analytic representation is convenient for a reaction calculation it may be obtained by a least-squares fit to a numerical orbital. 

Table 5.1 illustrates the frozen-core approximation for the case of sodium using a simple Slater (4.38) basis in the analytic-orbital representation. The core (Is 2s 2p ) is first calculated by Hartree—Fock for the state characterised by the 3s one-electron orbital, which we call the 3s state. The frozen-core calculation for the 3p state uses the same core orbitals and solves the 3p one-electron problem in the nonlocal potential (5.27) of the core. Comparison with the core and 3p orbitals from a 3p Hartree—Fock calculation illustrates the approximation. The overwhelming component of the 3p orbital agrees to almost five significant figures. 

The theory and computation of state-specific wavefuncfions for discrete and continuous atomic spectra are discussed in Section 5. These are constructed in terms of numerical as well as of analytic orbitals, each set... 

As a result, we advocate the use of plane-wave or blip [1] basis sets because they are easy to manipulate numerically and increase in number to improve accuracy. They are exponential type atomic orbitals with complex exponents. They may be used to represent the analytical orbitals over a numerical grid to arbitrary accuracy. 

The gradient of the PES (force) can in principle be calculated by finite difference methods. This is, however, extremely inefficient, requiring many evaluations of the wave function. Gradient methods in quantum chemistiy are fortunately now very advanced, and analytic gradients are available for a wide variety of ab initio methods [123-127]. Note that if the wave function depends on a set of parameters X], for example, the expansion coefficients of the basis functions used to build the orbitals in molecular orbital (MO) theory. 

Each of these factors can be viewed as combinations of CSFs with the same Cj and Cyj coefficients as in F but with the spin-orbital involving basis functions that have been differentiated with respect to displacement of center-a. It turns out that such derivatives of Gaussian basis orbitals can be carried out analytically (giving rise to new Gaussians with one higher and one lower 1-quantum number). 

Recently Thiel and Voityuk have constructed a workable NDDO model which also includes d-orbitals for use in connection with MNDO, called MNDO/d. With reference to the above description for MNDO/AM1/PM3, it is clear that there are immediately three new parameters Cd, Ud and (dd (eqs. (3.82) and (3.83)). Of the 12 new one-centre two-electron integrals only one (Gjd) is taken as a freely varied parameter. The other 11 are calculated analytically based on pseudo-orbital exponents, which are assigned so that the analytical formulas regenerate Gss, Gpp and Gdd. 

Many complexes of metals with organic ligands absorb in the visible part of the spectrum and are important in quantitative analysis. The colours arise from (i) d- d transitions within the metal ion (these usually produce absorptions of low intensity) and (ii) n->n and n n transitions within the ligand. Another type of transition referred to as charge-transfer may also be operative in which an electron is transferred between an orbital in the ligand and an unfilled orbital of the metal or vice versa. These give rise to more intense absorption bands which are of analytical importance. 

The Physiome Project should be viewed as both a vision and a route. It has been portrayed as consisting of two parts (Bassingthwaighte et al. 1998) (i) the databasing of biological information (the Mechanic s touch ), and (ii) the development of descriptive and, ultimately, analytical models of biological function (the Orbiter s view ). These are by no means sequential stages of the development. 

The analytical determination of the derivative dEtotldrir of the total energy Etot with respect to population n, of the r-th molecular orbital is a very complicated task in the case of methods like the BMV one for three reasons (a), those methods assume that the atomic orbital (AO) basis is non-orthogonal (b), they involve nonlinear expressions in the AO populations (c) the latter may have to be determined as Mulliken or Lbwdin population, if they must have a physical significance [6]. The rest of this paper is devoted to the presentation of that derivation on a scheme having the essential features of the BMV scheme, but simplified to keep control of the relation between the symbols introduced and their physical significance. Before devoting ourselves to that derivation, however, we with to mention the reason why the MO occupation should be treated in certain problems as a continuous variable. 

The method presented here allows, starting with trial gaussian functions, a partial analytical treatment which we have used to improve the LCAO-GTO orbitals (trial functions) essentially obtained from all ab initio quantum chemistry programs. As in r-representation, trial functions (t>i( Hp) (Eq. 21) are conveniently expressed as linear combinations of m functions Xi(P) themselves written as linear combinations of Gt gaussian functions (LCAO-GTO approximation) gta(P). 

---