Eigenvalue equation matrix form

Now that we have established the equivalence between an eigenvalue equation (operator form or matrix form) and a stationary-value condition (for functional variation or linear parameter variation) we turn to the general problem of finding the stationary values of the energy functional when all types of parameter, linear and non-linear, are admitted. A very convenient machinery for this purpose has been developed by Moccia (1974), whose approach we now adopt. 

Remember that ai is the representation of g(x) in the fi basis. So the operator eigenvalue equation is equivalent to the matrix eigenvalue problem if the functions fi form a complete set. 

In order to find the normal modes of vibration, I am going to write the above equations in matrix form, and then find the eigenvalues and eigenvectors of a certain matrix. In matrix form, we write... 

We want an eigenvalue equation because (cf. Section 4.3.4) we hope to be able to use the matrix form of a series of such equations to invoke matrix diagonalization to get eigenvalues and eigenvectors. Equation (5.35) is not quite an eigenvalue equation, because it is not of the form operation on function = k x function, but rather operation on function = sum of (k x functions). However, by transforming the molecular orbitals to a new set the equation can be put in eigenvalue form (with a caveat, as we shall see). Equation 5.35 represents a system of equations... 

The Hartree-Fock equations (5.47) (in matrix form Eqs. 5.44 and 5.46) are pseudoeigenvalue equations asserting that the Fock operator F acts on a wavefunction i//, to generate an energy value ,-, times i/q. Pseudoeigenvalue because, as stated above, in a true eigenvalue equation the operator is not dependent on the function on which it acts in the Hartree-Fock equations F depends on i// because (Eq. 5.36) the operator contains J and K, which in turn depend (Eqs. 5.29 and 5.30) on i//. Each of the equations in the set (5.47) is for a single electron ( electron 1 is indicated, but any ordinal number could be used), so the Hartree-Fock operator F is a one-electron operator, and each spatial molecular orbital i// is a one-electron function (of the coordinates of the electron). Two electrons can be placed in a spatial orbital because the, full description of each of these electrons requires a spin function 7 or jl (Section 5.2.3.1) and each electron moves in a different spin orbital. The result is that the two electrons in the spatial orbital i// do not have all four quantum numbers the same (for an atomic Is orbital, for example, one electron has quantum numbers n= 1, / = 0, m = 0 and s = 1/2, while the other has n= l,l = 0,m = 0 and s = —1/2), and so the Pauli exclusion principle is not violated. 

CCSD(T) method. The question then naturally arises as to how these methods can be extended to excited states. For the iterative methods, the extension is straightforward by analyzing the correspondence between terms in the CC equations and in H, one can define an H matrix for these methods, even though it is not exactly of the form of a similarity-transformed Hamiltonian. If one follows the linear-response approach, one arrives at the same matrix in the linear response theory, one starts from the CC equations, rather than the CC wave function, and no CC wave function is assumed. This matrix also arises in the equations for derivatives of CC amplitudes. In linear response theory, this matrix is sometimes called the Jacobian [19]. The upshot is that excited states for methods such as CCSDT-1, CCSDT-2, CCSDT-3, and CC3 can be obtained by solving eigenvalue equations in a manner similar to those for methods such as CCSD and CCSDT. 

The matrix eigenvalues are analyzed for the stationary state in respect to the concentrations of aU independent intermediates. In the thermody namic representation of kinetic equations, independent variables that describe the system evolution in time are not the concentrations but ther modynamic rushes A of the intermediates. Hence, the analysis of the sta bihty criterion in terms of thermodynamic variables needs an inspection of eigenvalues of matrix M p = However, the specific form of the... 

For every square matrix A of order m there is an eigenvalue, equation of the form Ax he. (4-4.1)... 

The matrix U is obtained by assuming that the matrix A satisfies three eigenvalue equations of the form... 

Thus, if A solves an eigenvalue equation of the form (E2.8), then a coordinate transformation can be found which will yield a diagonalized form of the matrix. The matrix A is itself a tensor (see Section E5.2) which can be constructed by taking the product of a column and a row vector... 

Substitutions of Eqs. (78-85) in Eq. (73) lead to three-term recurrence relations of the expansion coefficients contained in Eqs. (32-39) in Ref. [6] and not reproduced here for the sake of space. The three-term recurrence relations can be cast onto the form of matrix eigenvalue equations with eigenvalues li and the eigenvectors, just like in the case of Section 2.3. 

These results indicate that the double-commutator Liouvillian matrix has rather different properties compared to the ordinary matrix (2.64), and that the eigenvalues in the form of energy differences have a different form. It should also be observed that the solutions to the linear equation system (2.83) may have a rather different character.13... 

To apply these equations, we need the wavefuncdons m> in order to get the dipole moment transition elements and the frequencies spectral series, where only the ground state need be near-exact. This is done by diagonalizing the Hamiltonian matrix formed from a large number of basis functions (which implicitly include the interelectronic coordinate and thus electron correlation). We do this for each symmetry state that is involved. All the ensuing eigenvalues and eigenvectors are then used in the sum-over-states expressions. For helium we require S, P, and D states and for H2 (or D2) E, II, and A states. 

We subdivide the system into two parts, one built from metal nd orbitals (d) and another composed of valence metal ( + l)s and ( + l)p and ligand functions (v). Then the eigenvalue problem (Equation 7) can be represented in the form of the pseudo-eigenvalue Equation 8, completely restricted to the d-subspace, with the explicit form of given by Equation 9. is a A(j x Nj matrix. Thus, for a d system, for instance, Nj = 45. (H j) are rectangular... 

Written in matrix-form, the new eigenvalue equation has a simple structure, with all r-dependent parts in the off-diagonal of the matrix ... 

The derived EOM equations (6) [(7)] are linear matrix eigenvalue equations for the exact excitation energies (ionization potentials and electron affinities). The eigenvectors give an satisfying conditions (12) and (13) (with o 0 =0 for IPs and EAs). The most general Oj is expressed in terms of the basis operators X >[Pg.12]

The exchange matrix, K, is just the rate, k, times the unit matrix. In block form, the full matrix for two sites is given in the eigenvalue equation, (B2.4.38). 

The method is of algebraic order 6. It has been proved that for a = 3355 the phase-lag is of order 8. The interval of periodicity is equal to 0,6(231 9 -"" - = (0. 7.26).This method is rather more complicated in its matrix form. The eigenvalue equation now takes the form... 

Here Fdimer and Sdimer have the same definitions as before, the unitary matrix U is formed from the solutions Uj of the generalized eigenvalue equation... 

---