Jacobi transformation

THE JACOBI TRANSFORMATION, DIAGONALIZATION OF A SYMMETRIC MATRIX AND CANONICAL ORTHOGONALIZATION... 

Calculate the rotation angle for the Jacobi transformation of equation 3.25, with the conditional formula... 

For the jsto-ng) 2s and Is basis sets of Table 1.6, equation 3.15 under the Jacobi transformation becomes the diagonal matrix of cells H 10 to I 11 in Figure 3.15 and 3.16,... 

Secondly, the canonical orthonormalization procedure to diagonalize the overlap matrix and then the application of the Jacobi transformation to diagonalize the Fock matrix in the eigenfunctions of the overlap matrix, returns two eigenvalues, the values —0.50000 and —0.12352 Hartrees, in canonical B 18 and B 19. This is the important elementary point that we can make two linear combinations of two functions and so there are two possible eigenvalues to be calculated. These eigenvalues, of course, are present in the calculation set out in the other worksheet, based on the Schmidt procedure. The Is... 

Figure 4.20a Application of the Jacobi transformation and the canonical orthonormalization procedure to the calculation of the Is orbital energy in hydrogen using the Dunning (4s 2s) Gaussian basis (46). Note that the calculation returns a second set of coefficients in cells H 22 and H 23, which are simply a product of the calculation procedure and give rise to a primitive approximation to the 2s orbital (22).

Exercise 4.15. Application of the Jacobi transformation and canonical orthonormalization to the calculation of the Is orbital energy in hydrogen using the Dunning (4s 2s) basis. 

Partial derivatives (dfcan be transformed to (df/Byi)yjjiyi with the help of functional determinants (Jacobi transformation) if the functions jc, = x, (yy ) are known. For practical use all partial derivatives of energy functions U, H, A, and G and of entropy S are reduced to functions of the tabulated material properties a (thermal expansivity coefficient), p (isothermal compresibility coefficient), and Cp (heat capacity at constant pressure) ... 

Proof of Proposition I.2.I It sufflces to prove the assertion for a domain D which has the form of an infinitesimal rectangle. It is well known that the change of the rectangle area under a shift is measured by the determinant of the Jacobi transformation matrix. Consequently, it suffices to calculate this determinant for an infinitesimal transformation... 

This matrix is positive semi-definite and so it is expected to have three positive eigenvalues. The solution of this cubic equation is not a trivial matter. Instead of estimating the roots of the characteristic equation in this case, it is preferable to use a more general method and estimate the eigenvalues of matrix A, defined by equation (8), by use of the Jacobi transformation of the matrix. The Jacobi transformation (or orthogonalisation) is an iterative application of rotations to a matrix until all its off-diagonal values are zero at machine precision. 

By the introduction of the (x, y) coordinate system, one has reduced the problem to the motion of a particle of mass (i in a two-dimensional rectilinear space (x, y). Thus, the problem of the collision between an atom and a diatomic molecule in a collinear geometry has been converted into a problem of a single particle on the potential energy surface expressed in terms of the coordinates x and y rather than the coordinates rAB and rBc The coordinates x and y which transform the kinetic energy to diagonal form in such way that the kinetic energy contains only one (effective) mass are referred to as mass scaled Jacobi coordinates. 

Although diagrams like Fig. 6.1 are especially convenient to illustrate the qualitative features of TST and VTST, the solution of the equations of motion in (rAB,rBc) coordinates is complicated due to cross terms coupling the motions of the different species. It is for that reason we introduced mass scaled Jacobi coordinates in order to simplify the equations of motion. So, one now asks what does the potential function for reaction between A and BC look like in these new mass scaled Jacobi coordinates. To illustrate we construct a graph with axes designated rAB and rBc within the (x,y) coordinate system. In the x,y space lines of constant y are parallel to the x axis while lines of constant x are parallel to the y axis. The rAB and rBc axes are constructed in similar fashion. Lines of constant rBc are parallel to the rAB axis while lines of constant rAB are parallel are parallel to the rBc axis. From the above transformation, Equations 6.10 to 6.13... 

These two systems are examples from non-linear physics, where the equations can be solved in terms of elliptic functions and elliptic integrals. The reader who is not familiar with these functions, which do not arise in the same way as the previously mentioned special functions, is referred to the excellent book by Whittaker and Watson [6]. In that book, the reader will see that there are two flavours of elliptic functions, the Weierstrass and Jacobi representations, three kinds of elliptic integrals, and six kinds of pseudo-periodic functions, the Weierstrass zeta and sigma functions and the four kinds of Jacobi theta functions. Of historical interest for theoretical chemists is the fact that Jacobi s imaginary transformation of the theta functions is the same as the Ewald transformation of crystal physics [7]. 

If reactant coordinates are to be used for the propagation, then this function is propagated forward in time and analyzed [151] after it has reached the product region at large values of r. If state-to-state reactive probabilities and cross sections are required, the initial wavepacket must be transformed to the product Jacobi coordinates and the propagation must be performed in these coordinates. The general form of the transformation from reactant to product Jacobi coordinates is... 

In order to transform to the body-fixed representation, we will need to relate the angular functions Wj (R,r) to angular functions defined relative to the body-fixed axes [L., J,K,M,p)QjK ), where J,K,M,p) are the parity-adapted total angular momentum eigenfunctions of Eq. (4.5) and x(0) normalized associated Legendre polynomials of the body-fixed Jacobi angle]. 

The development just given illustrates the fact that the topology of the vacuum determines the nature of the gauge transformation, field tensor, and field equations, as inferred in Section (I). The covariant derivative plays a central role in each case for example, the homogeneous field equation of 0(3) electrodynamics is a Jacobi identity made up of covariant derivatives in an internal 0(3) symmetry gauge group. The equivalent of the Jacobi identity in general relativity is the Bianchi identity. 

The absence of an E(3) field does not affect Lorentz symmetry, because in free space, the field equations of both 0(3) electrodynamics are Lorentz-invariant, so their solutions are also Lorentz-invariant. This conclusion follows from the Jacobi identity (30), which is an identity for all group symmetries. The right-hand side is zero, and so the left-hand side is zero and invariant under the general Lorentz transformation [6], consisting of boosts, rotations, and space-time translations. It follows that the B<3) field in free space Lorentz-invariant, and also that the definition (38) is invariant. The E(3) field is zero and is also invariant thus, B(3) is the same for all observers and E(3) is zero for all observers. 

The eigenvalues of A can be find by solving the characteristic equation of (1.61). It is much more efficient to look for similarity transformations that will translate A into the diagonal form with the eigenvalues in the diagonal. The Jacobi method involves a sequence of orthonormal similarity transformations, 12,... such that A(<+1 = TTkAkTk. The matrix Tk differs from the identity... 

In the Jacobi method, a series of similarity transformations is carried out. It is easily proven that similar matrices have the same eigenvalues. Let A = P, BP. The eigenvalues of A satisfy the secular equation (2.38) ... 

Usually on the first iteration of an SCF calculation W is computed by the Schmidt orthogonalization method but thereafter W is chosen to be the C matrix from the previous iteration. This produces an F matrix which is nearly diagonal so the Jacobi method becomes quite efficient after the first iteration. Further, in the Jacobi method, F is diagonalized by an iterative sequence of simple plane-rotation transformations... 

The Jacobi method works by successively transforming the matrices A and v in place , in a way which ensures that the non-diagonality measure... 

In the 0(3) invariant theory of the Aharonov-Bohm effect, the holonomy consists of parallel transport using 0(3) covariant derivatives and the internal gauge space is a physical space of three dimensions represented in the basis ((1),(2),(3)). Therefore, a rotation in the internal gauge space is a physical rotation, and causes a gauge transformation. The core of the 0(3) invariant explanation of the Aharonov-Bohm effect is that the Jacobi identity of covariant derivatives [46]... 

This coupling potential is smooth everywhere, which allows numerical calculations with high precision. There is no nonadiabatic coupling since the basis functions [0< )( 2C) are independent of p in each sector. The solution I Wf/o, 2C) is connected smoothly, in principle, from sector to sector by a unitary frame transformation from the /th set of channels to the (/ + l)st set [97-99]. The coordinate system is transformed from the hyperspherical to the Jacobi coordinates at some large p, beyond which the conventional close-coupling equations are employed for determining the asymptotic form of the wavefunction appropriate for the scattering boundary condition [100]. 

In order to simplify the evaluation of overlap integrals between bound and continuum wavefunctions, it is advisable (although not necessary) to describe both wavefunctions by the same set of coordinates. Usually, the calculation of continuum, i.e., scattering, states causes far more problems than the calculation of bound states and therefore it is beneficial to use Jacobi coordinates for both nuclear wavefunctions. If bound and continuum wavefunctions are described by different coordinate sets, the evaluation of multi-dimensional overlap integrals requires complicated coordinate transformations (Freed and Band 1977) which unnecessarily obscure the underlying dynamics. 

---