Non-adiabatic coupling matrix

Some final comments on the relevance of non-adiabatic coupling matrix elements to the nature of the vector potential a are in order. The above analysis of the implications of the Aharonov coupling scheme for the single-surface nuclear dynamics shows that the off-diagonal operator A provides nonzero contiibutions only via the term (n A n). There are therefore no necessary contributions to a from the non-adiabatic coupling. However, as discussed earlier, in Section IV [see Eqs. (34)-(36)] in the context of the x e Jahn-Teller model, the phase choice t / = —4>/2 coupled with the identity... 

A. The Quantization of the Non-Adiabatic Coupling Matrix Along a Closed Path... 

Hence, in order to contract extended BO approximated equations for an N-state coupled BO system that takes into account the non-adiabatic coupling terms, we have to solve N uncoupled differential equations, all related to the electronic ground state but with different eigenvalues of the non-adiabatic coupling matrix. These uncoupled equations can yield meaningful physical... 

We concentrate on an adiabatic tri-state model in order to derive the quantization conditions to be fulfilled by the eigenvalues of the non-adiabatic coupling matrix and finally present the extended BO equation. The starting point is the 3x3 non-adiabatic coupling matrix,... 

The adiabatic coupled SE for the above 3x3 non-adiabatic coupling matrix are... 

The matrix of vectors F is thus the defining quantity, and is called the non-adiabatic coupling matrix. It gives the strength (and direction) of the coupling between the nuclear functions associated with the adiabatic electronic states. 

The elements of the matrix G can be written in terms of F, which is called the non-adiabatic coupling matrix. For a particular coordinate, a, and dropping the subscript for clarity,... 

One of the main outcomes of the analysis so far is that the topological matrix D, presented in Eq. (38), is identical to an adiabatic-to-diabatic transformation matrix calculated at the end point of a closed contour. From Eq. (38), it is noticed that D does not depend on any particular point along the contour but on the contour itself. Since the integration is carried out over the non-adiabatic coupling matrix, x, and since D has to be a diagonal matrix with numbers of norm 1 for any contour in configuration space, these two facts impose severe restrictions on the non-adiabatic coupling terms. 

In this section, we intend to show that for a certain type of models the above imposed restrictions become the ordinary well-known Bohr-Sommerfeld quantization conditions [82]. For this purpose, we consider the following non-adiabatic coupling matrix x ... 

The non-adiabatic coupling matrix t will be defined in a way similar to that in the Section V.A [see Eq. (51)], namely, as a product between a vector-function t(i) and a constant antisymmetric matrix g written in the form... 

We intend to show that an adiabatic-to-diabatic transformation matrix based on the non-adiabatic coupling matrix can be used not only for reaching the diabatic fi amework but also as a mean to determine the minimum size of a sub-Hilbert space, namely, the minimal M value that still guarantees a valid diabatization. 

Equations (169) and (171), together with Eqs. (170), fomi the basic equations that enable the calculation of the non-adiabatic coupling matrix. As is noticed, this set of equations creates a hierarchy of approximations starting with the assumption that the cross-products on the right-hand side of Eq. (171) have small values because at any point in configuration space at least one of the multipliers in the product is small [115]. 

Figure 12, Results for the C2H molecule as calculated along a circle surrounding Che 2 A -3 A conical intersection, The conical intersection is located on the C2v line at a distance of 1,813 A from the CC axis, where ri (=CC distance) 1.2515 A. The center of the circle is located at the point of the conical intersection and defined in terms of a radius < . Shown are the non-adiabatic coupling matrix elements tcp((p ) and the adiabatic-to-diabatic transformation angles y((p i ) as calculated for (ii) and (b) where q = 0.2 A (c) and (d) where q = 0.3 A (e) and (/) where q = 0.4 A. Also shown are the positions of the two close-by (3,4) conical intersections (designated as X34).

In Section IV, we introduced the topological matrix D [see Eq. (38)] and showed that for a sub-Hilbert space this matrix is diagonal with (-1-1) and (—1) terms a feature that was defined as quantization of the non-adiabatic coupling matrix. If the present three-state system forms a sub-Hilbert space the resulting D matrix has to be a diagonal matrix as just mentioned. From Eq. (38) it is noticed that the D matrix is calculated along contours, F, that surround conical intersections. Our task in this section is to calculate the D matrix and we do this, again, for circular contours. 

---