Secular matrix

Later in this book, we shall need to find the roots of the secular matrix... 

This equation is a quadratic and has two roots. For quantum mechanical reasons, we are interested only in the lower root. By inspection, x = 0 leads to a large number on the left of Eq. (1-10). Letting x = leads to a smaller number on the left of Eq. (1-10), but it is still greater than zero. Evidently, increasing a approaches a solution of Eq. (1-10), that is, a value of a for which both sides are equal. By systematically increasing a beyond 1, we will approach one of the roots of the secular matrix. Negative values of x cause the left side of Eq. (1-10) to increase without limit hence the root we are approaching must be the lower root. 

If we divide each element of the secular matrix by p and perform the substitution... 

Write the secular matrix, compute the eigenvalues, and draw the energy level diagram for fulvcnc. 

One restriction imposed by Huckel theory that is rather easy to release is that of zero overlap for nearest-neighbor interactions. One can retain a — as the diagonal elements in the secular matrix and replace p by p — EjS as nearest-neighbor elements where S is the overlap integral. Now,... 

Hoffman s extended Huckel theory, EHT (Hoffman, 1963), includes all bonding orbitals in the secular matrix rather than just all n bonding orbitals. This inclusion increases the complexity of the calculations so that they are not practical without a computer. The basis set is a linear combination that includes only valence orbitals... 

The orbitals used for methane, for example, are four Is Slater orbitals of hydrogen and one 2s and three 2p Slater orbitals of carbon, leading to an 8 x 8 secular matrix. Slater orbitals are systematic approximations to atomic orbitals that are widely used in computer applications. We will investigate Slater orbitals in more detail in later chapters. 

We fill the secular matrix H with elements Hy over the entire set of valence orbitals. The diagonal elements are... 

There are two functions, so we shall obtain two eigenvalues. The ground-state energy will be the lower of the two. The full secular matrix is... 

The various basis sets used in a calculation of the H and S integrals for a system are attempts to obtain a basis set that is as close as possible to a complete set but to stay within practical limits set by the speed and memory of contemporary computers. One immediately notices that the enterprise is directly dependent on the capabilities of available computers, which have become more powerful over the past several decades. The size and complexity of basis sets in common use have increased accordingly. Whatever basis set we choose, however, we are attempting to strike a balance. If the basis set is too small, it is inaeeurate if it is too large, it exceeds the capabilities of our computer. Whether our basis set is large or small, if we attempt to calculate all the H and S integrals in the secular matrix without any infusion of empirical information, the procedure is described as ab initio. 

The canonical molecular orbitals of any molecule can by obtained by computer calculations. All MO methods involve the diagonalization of a secular matrix. It can be said that by moving from AOs to FOs to BOs basis sets one proceeds through the various stages of this diagonalization process, as the number of non-zero off-diagonal overlap matrix elements decreases. 

The local mode Hamiltonian (5.16) includes only the operator C[2 that is, interactions of the Casimir type between bonds 1 and 2. One may wish, in some cases, to include also interactions of this type between bonds 1 and 3, C13, and 2 and 3, C23. These can be included by diagonalizing the secular matrix obtained by evaluating the matrix elements of C13 and C23 in the basis (5.4). These matrix elements are given by (5.15a) and (5.15b). 

This operator is used only for the construction of the d matrix. The gradient terms in the secular matrix are, of course, computed exactly. 

In fact the relative coefficients within a set of symmetrically equivalent atoms, such as Cl, C3, C5, and Cg in para-benzosemiquinone, can be determined by group theory alone. The appropriate set of symmetrized orbitals, also listed in Table 1, can be obtained by use of character tables and procedures described in Refs. 3 to 6. In matrix notation, the symmetrized combinations are <1> = Up, where <1) and p are column vectors and U is the transformation matrix giving the relations shown in part (c) of Table 1. The transformation U and its transpose U can be used to simplify the solution of the secular matrix X since the matrix multiplication UXU gives the block diagonal form shown in part (6) of Table 1. 

To calculate numerically the quantum dynamics of the various cations in time-dependent domain, we shall use the multiconfiguration time-dependent Hartree method (MCTDH) [79-82, 113, 114]. This method for propagating multidimensional wave packets is one of the most powerful techniques currently available. For an overview of the capabilities and applications of the MCTDH method we refer to a recent book [114]. Additional insight into the vibronic dynamics can be achieved by performing time-independent calculations. To this end Lanczos algorithm [115,116] is a very suitable algorithm for our purposes because of the structural sparsity of the Hamiltonian secular matrix and the matrix-vector multiplication routine is very efficient to implement [6]. 

Application of the variational method leads the secular matrix equation... 

One sees from the functional dependence of the energy, eq (44), on trigonometric functions that the energies of the bands are of the standard form. Because of the one-dimensional nature of the problem and the 2x2 dimensions of the secular matrix, it is possible simply to include account of overlap automatically without specifically using methods of orthonormalization, such as Lowdin s method [42-44]. 

In constrast to the matrix H - E0 in (1.19), the secular matrix has a complicated, nonlinear energy dependence, and the one-electron energies E must be found individually by tracing the roots of the determinant of as a 10... 

In 1973 Andersen and Woolley [1.25] extended the LCMTO method to molecular calculations. At the end of their paper they introduced that choice of MTO tail, i.e. proportional to J = 9i/j/9E, which in a natural fashion ensured orthogonality to the core states and at the same time led to an accurate and elegant formulation of linear methods. The resulting, technique was immediately developed in a paper by Andersen [1.26] which, in a condensed form, contains most of what one need know about the simple concepts of linear band theory. Thus, we find here the KKR equation within the atomic-sphere approximation at this stage is called ASM the LCMTO secular matrix, latter called the LMTO matrix the energy-independent structure constants and the canonical bands and the Laurent expansion of the logarithmic-derivative function and the corresponding potential parameters. 

In Chap.5 we derive the LCMTO equations in a form not restricted to the atomic-sphere approximation, and use the , technique introduced in Chap.3 to turn these equations into the linear muffin-tin orbital method. Here we also give a description of the partial waves and the muffin-tin orbitals for a single muffin-tin sphere, define the energy-independent muffin-tin orbitals and present the LMTO secular matrix in the form used in the actual programming, Sect.9.3. 

In Chap.6 the atomic-sphere approximation is introduced and discussed, canonical structure constants are presented, and it is shown that the LMTO-ASA and KKR-ASA equations are mathematically equivalent in the sense that the KKR-ASA matrix is a factor of the LMTO-ASA secular matrix. In addition, we treat muffin-tin orbitals in the ASA, project out the i character of the eigenvectors, derive expressions for the spherically averaged electron density, and develop a correction to the ASA. 

---