Matrix orthonormality

U(qJ is referred to as an adiabatic-to-diabatic transformation (ADT) matrix. Its mathematical sbucture is discussed in detail in Section in.C. If the electronic wave functions in the adiabatic and diabatic representations are chosen to be real, as is normally the case, U(q ) is orthogonal and therefore has n n — l)/2 independent elements (or degrees of freedom). This transformation mabix U(qO can be chosen so as to yield a diabatic electronic basis set with desired properties, which can then be used to derive the diabatic nuclear motion Schrodinger equation. By using Eqs. (27) and (28) and the orthonormality of the diabatic and adiabatic electronic basis sets, we can relate the adiabatic and diabatic nuclear wave functions through the same n-dimensional unitary transformation matrix U(qx) according to... 

If the functions Oj are orthonormal, then the overlap matrix S reduces to the unit matrix and the above generalized eigenvalue problem reduces to the more familiar form ... 

For the hermitian matrix in review exereise 3b show that the pair of degenerate eigenvalues ean be made to have orthonormal eigenfunetions. 

After forming the overlap matrix, the new orthonormal funetions are defined as follows ... 

In the unlikely event that none of the basis functions overlap, then S is a unit matrix. We usually require the LCAO orbitals . to be orthonormal... 

Roothaan actually solved the problem by allowing the LCAO coefficients to vary, subject to the LCAO orbitals remaining orthonormal. He showed that the LCAO coefficients are given from the following matrix eigenvalue equation ... 

We then allow Ri and R2 to vary, subject to orthonormality, just as in the closed-shell case. Just as in the closed-shell case, Roothaan (1960) showed how to write a Hamiltonian matrix whose eigenvectors give the columns U] and U2 above. 

The bra n denotes a complex conjugate wave function with quantum number n standing to the of the operator, while the ket m), denotes a wave function with quantum number m standing to the right of the operator, and the combined bracket denotes that the whole expression should be integrated over all coordinates. Such a bracket is often referred to as a matrix element. The orthonormality condition eq. (3.5) can then be written as. 

Quadralically Convergent or Second-Order SCF. As mentioned in Section 3.6, the variational procedure can be formulated in terms of an exponential transformation of the MOs, with the (independent) variational parameters contained in an X matrix. Note that the X variables are preferred over the MO coefficients in eq. (3.48) for optimization, since the latter are not independent (the MOs must be orthonormal). The exponential may be written as a series expansion, and the energy expanded in terms of the X variables describing the occupied-virtual mixing of the orbitals. 

Finally, since the MOs are orthonormal, the derivatives of the coefficients may be replaced by derivatives of the overlap matrix. 

The variational problem is to minimize the energy of a single Slater determinant by choosing suitable values for the MO coefficients, under the constraint that the MOs remain orthonormal. With cj) being an MO written as a linear combination of the basis functions (atomic orbitals) /, this leads to a set of secular equations, F being the Fock matrix, S the overlap matrix and C containing the MO coefficients (Section 3.5). 

In order to try to approach the HF scheme as much as possible, we will now introduce the basic orthonormal set fc which has maximum occupation numbers. Let U be the unitary matrix which brings the hermitean matrix (ylk) to diagonal form ... 

This procedure would generate the density amplitudes for each n, and the density operator would follow as a sum over all the states initially populated. This does not however assure that the terms in the density operator will be orthonormal, which can complicate the calculation of expectation values. Orthonormality can be imposed during calculations by working with a basis set of N states collected in the Nxl row matrix (f) which includes states evolved from the initially populated states and other states chosen to describe the amplitudes over time, all forming an orthonormal set. Then in a matrix notation, (f) = (f)T (t), where the coefficients T form IxN column matrices, with ones or zeros as their elements at the initial time. They are chosen so that the square NxN matrix T(f) = [T (f)] is unitary, to satisfy orthonormality over time. Replacing the trial functions in the TDVP one obtains coupled differential equations in time for the coefficient matrices. 

After solving this equation for the density matrix, its diagonalization provides the matrices of coefficients and the weights needed to reconstruct orthonormal density amplitudes. The density operator follows from r(f) = (<))r(f)( (f) ... 

The matrix of a self-adjoint operator in any orthonormal basis is a symmetric matrix. 

Since Doo and are constructed with the same set of orthonormal spinorbitals, the two first matrix elements can easily rewritten, according to the Slater s rules [13], as ... 

The calculation of the cross matrix elements (6) is somewhat more difficult, because the Slater Determinants involved in them are constructed with two sets of non-orthonormal spinorbitals. This calculation, however, may be greatly simplified, if the two sets are assumed to be corresponding, that is, if they fulfill the following condition [14] ... 

So far we have considered an orthonormal basis set x In actnal calcnlations, employing non orthogonal sets of Gaussian fnnctions with overlap matrix... 

The columns of the orthonormal matrix Vp are linear combinations of reaction invariants . In fact, the only invariants for the batch reaction being analyzed can be stoichiometric coefficients. Hence the matrix Vp may be interpreted as containing the stoichiometric information (Waller and Makila (1981)) and its rank Nr can be considered to be equal to the number of independent... 

A matrix is called orthonormal if additionally we obtain that ... 

The 3x3 matrix U shown below is both row- and column-orthonormal ... 

---