Construction of basis functions

Because H and its eigenvalues are invariant when a symmetry operator T acts on the physical system, T j s is also an eigenfunction of II with the same eigenvalue E, and therefore it is a linear combination of the (j s,  

The problem is this how can we generate j ]p, p=l,. ..,/,-, the set of lj orthonormal functions which form a basis for the yth IR of the group of the Schrodinger equation We start with any arbitrary function o defined in the space in which the set of function operators T operate. Then 

---

The utility of such lists is emphasised in a further example, the construction of basis functions transforming with Tiu symmetry for the two orbit fullerene, Cso. The 80-vertex structure is formed by combination of the O20 orbit and an Oeo orbit of In point symmetry to realize the three-valent fullerene. The permutation character over the full set of 80 vertices is the direct sum... 

The situation is quite similar in chemistry. Due to decades of experience with Hartree-Fock and Cl calculations much is known about the construction of basis functions that are suitable for molecules. Almost all of this continues to hold in DFT — a fact that has greatly contributed to the recent popularity of DFT in chemistry. Chemical basis functions are classified with respect to their behaviour as a function of the radial coordinate into Slater type orbitals (STOs), which decay exponentially far from the origin, and Gaussian type orbitals (GTOs), which have a gaussian behaviour. STOs more closely resemble the true behaviour of atomic wave functions [in particular the cusp condition of Eq. (19)], but GTOs are easier to handle numerically because the product of two GTOs located at different atoms is another GTO located in between, whereas the product of two STOs is not an STO. The so-called contracted basis functions , in which STO basis functions are reexpanded in... 

It is now necessary to construct linear combinations of the d-orbitals which transform according to the representations eg and t2g. Generally speaking, the construction of basis functions may be quite tedious, apart from a number of simple cases where it may be done practically by inspection (as for example, the one-dimensional representations). Basis sets for the common situations are tabulated in various places e. g., Koster et al. (34), Ballhausen (2), Griffith (21). It will be sufficient for our purpose to give the basis sets for eg and t2g and to demonstrate that they satisfy the necessary requirements. Since our discussion is confined to systems of d-electrons, all states will be of even parity or g-states. To simplify the notation we shall henceforth suppress the parity index, unless specifically needed. 

The phrase symmetry adapted basis functions refers to those linear combinations of basis functions (on several atoms) that transform like the particular irreducible representation of the appropriate point group. Molecular symmetry is used at various points in these calculations twenty years ago I would have had to write several chapters on molecular symmetry, point groups, constructing symmetry-adapted combinations of basis functions, factoring a Hamiltonian matrix using symmetry and related topics. The point is that twenty... 

To construct the Fock matrix, eq. (3.51), integrals over all pairs of basis functions and the one-electron operator h are needed. For M basis functions there are of the order of of such one-electron integrals. These one-integrals are also known as core integrals, they describe the interaction of an electron with the whole frame of bare nuclei. The second part of the Fock matrix involves integrals over four basis functions and the g two-electron operator. There are of the order of of these two-electron integrals. In conventional HF methods the two-electron integrals are calculated and saved before the... 

The conclusion is that for every particular set of basis functions and given data, there exists an appropriate size of G that can approximate both accurately and smoothly this data set. A decisive advantage would be if there existed a set of basis functions, which could probably represent any data set or function with minimal complexity (as measured by the number of basis functions for given accuracy). It is, however, straightforward to construct different examples that acquire minimal representations with respect to different types of basis functions. Each basis function for itself is the most obvious positive example. A Gaussian (or discrete points... 

The simplest approximation corresponds to a single-determinant wave function. The best possible approximation of this type is the Hartree-Fock (HF) molecular-orbital determinant. The HF wavefunction is constructed from the minimal number of occupied MOs (i.e., NI2 for an V-eleclron closed-shell system), each approximated as a variational linear combination of the chosen set of basis functions (vide infra). 

From this wave function, one sees how even in the early beginning of molecular quantum mechanics, atomic orbitals were used to construct molecular wave functions. This explains why one of the first AIM definitions relied on atomic orbitals. Nowadays, molecular ab initio calculations are usually carried out using basis sets consisting of basis functions that mimic atomic orbitals. Expanding the electron density in the set of natural orbitals and introducing the basis function expansion leads to [15]... 

Now, if a suitable space of basis functions is used (a space of basis functions that is closed under the symmetry operations of the group), we can construct a set of representations (each one consisting of 48 matrices) for this space that is particularly useful for our purposes. It is especially relevant that the matrices of each one of these representations can be made equivalent to matrices of lower dimensions. 

This set of representations is usually known as a representation of the group. Obviously, if we choose anotiier space of basis functions., anotiier representation of the group can be constructed, and so an infinite number of representations is possible for a given symmetry group. 

As an example, to construct the character table for the Oh symmetry group we could apply the symmetry operations of the ABg center over a particularly suitable set of basis functions the orbital wavefunctions s, p, d,... of atom (ion) A. These orbitals are real functions (linear combinations of the imaginary atomic functions) and the electron density probability can be spatially represented. In such a way, it is easy to understand the effect of symmetry transformations over these atomic functions. 

Projection operators are a technique for constructing linear combinations of basis functions that transform according to irreducible representations of a group. Projection operators can be used to form molecular orbitals from a basis set of atomic orbitals, or to form normal modes of vibration from a basis of displacement vectors. With projection operators we can revisit a number of topics considered previously but which can now be treated in a uniform way. 

Gaussian. A function of the form x y z" exp (ar ) where 1, m, n are integers (0,1,2. . .) and a is a constant. Used in the construction of Basis Sets for Hartree-Fock, Density Functional, MP2 and other Correlated Models. 

The number of basis functions (defined by the chosen basis sets) used to construct the molecular orbitals also strongly affects the effort/accuracy ratio. The use of minimal basis sets yielded wrong results (56), whereas reasonable agreement with experiment is obtained when double zeta plus polarization basis sets are applied. Correlated methods require larger basis sets to include as much electron correlation as possible. This implies that in addition to the increased computational demand of such methods, a further increase of the computational cost results due to the requirement of using larger basis sets. 

The basis set is the set of madiematical functions from which the wave function is constructed. As detailed in Chapter 4, each MO in HF theory is expressed as a linear combination of basis functions, the coefficients for which are determined from the iterative solution of the HF SCF equations (as flow-charted in Figure 4.3). The full HF wave function is expressed as a Slater determinant formed from the individual occupied MOs. In the abstract, the HF limit is achieved by use of an infinite basis set, which necessarily permits an optimal description of the electron probability density. In practice, however, one cannot make use of an infinite basis set. Thus, much work has gone into identifying mathematical functions that allow wave functions to approach the HF limit arbitrarily closely in as efficient a manner as possible. 

Note, however, that since we now work with only the trace of the matrix, we have no information about off-diagonal elements of the irrep matrices and hence no way to construct shift operators. The business of establishing symmetry-adapted functions therefore involves somewhat more triad and error than the approach detailed above. Character projection necessarily yields a function that transforms according to the desired irrep (or zero, of course), but application of character projection to different functions will be required to obtain a set of basis functions for a degenerate irrep, and the resulting basis functions need not be symmetry adapted for the full symmetry species (irrep and row) obtained above. 

A technique for direct computations of the eigenvalues Er —zT/2 of H(6 = 0) with the outgoing-wave boundary condition is reviewed in detail in a chapter in Part I of this two-volume special issue of Advances in Quantum Chemistry on Unstable States in the Continuous Spectra [27]. Determination of the wavefunction of Eq. (2) with a real eigenvalue EQ using a judiciously chosen real, square-integrable basis set, followed by diagonalization of a complex Hamiltonian matrix for the whole Hq + H constructed in terms of basis functions of complex-rotated coordinates, is shown to be quite useful. 

A means to extrapolate the attractive part of Udip(p) is to remove from the adiabatic potentials the contributions from the arrangement e+ + He+. This may be effected by constructing the basis functions for the adiabatic hyperangular states in Eq. (93) in terms only of the Jacobi coordinates for the arrangement Ps + He. The two lowest potential curves (henceforth denoted by V2(p) and V3(p)) obtained in this way are included in Figure 4.16 as red broken curves. The curve approaching Ps(n = 3), i.e., V3(p), can be... 

In 1951 Roothaan and Hall independently pointed out [26] that these problems can be solved by representing MO s as linear combinations of basis functions (just as in the simple Hiickel method, in Chapter 4, the % MO s are constructed from atomic p orbitals). Roothaan s paper was more general and more detailed than Hall s, which was oriented to semiempirical calculations and alkanes, and the method is sometimes called the Roothaan method. For a basis-function expansion of MO s we write... 

---