Nonadiabatic couplings

Dardi P S and Dahler J S 1993 A model for nonadiabatic coupling in the photodissociation of l2-solvent complexes J. Chem. Phys. 98 363-72... 

The full dynamical treatment of electrons and nuclei together in a laboratory system of coordinates is computationally intensive and difficult. However, the availability of multiprocessor computers and detailed attention to the development of efficient software, such as ENDyne, which can be maintained and debugged continually when new features are added, make END a viable alternative among methods for the study of molecular processes. Eurthemiore, when the application of END is compared to the total effort of accurate determination of relevant potential energy surfaces and nonadiabatic coupling terms, faithful analytical fitting and interpolation of the common pointwise representation of surfaces and coupling terms, and the solution of the coupled dynamical equations in a suitable internal coordinates, the computational effort of END is competitive. 

Ignoring all nonadiabatic couplings to higher electronic states, the nuclear motion in a two-state elechonic manifold is described explicitly in the adiabatic representation by... 

By following Section II.B, we shall be more specific about what is meant by strong and weak interactions. It turns out that such a criterion can be assumed, based on whether two consecutive states do, or do not, form a conical intersection or a parabolical intersection (it is important to mention that only consecutive states can form these intersections). The two types of intersections are characterized by the fact that the nonadiabatic coupling terms, at the points of the intersection, become infinite (these points can be considered as the black holes in molecular systems and it is mainly through these black holes that electronic states interact with each other.). Based on what was said so far we suggest breaking up complete Hilbert space of size A into L sub-Hilbert spaces of varying sizes Np,P = 1,..., L where... 

Figure 6, Equi-nonadiabatic coupling lines for the terms ri2(r,v) and r23(r,y) as calculated for the C2H molecule for a fixed C—C distauce, that is, rcc — 1-35 A. (a) Equi-non-adiabaric coupling term lines for the ri2(r,y). (b) Equi-non-adiabatic coupling term lines for X23(r,y). The Cartesian coordinates (x,y) are related to (q, 9) as follows x — qcosO] y —, sin9, where q and 9 are measured with respect to the midpoint between the two carbons.

The extended Bom-Oppenheimer approximation based on the nonadiabatic coupling terms was discussed on several occasions [23,25,26,55,56,133,134] and is also presented here by Adhikari and Billing (see Chapter 3). 

In Chapter VI, Ohm and Deumens present their electron nuclear dynamics (END) time-dependent, nonadiabatic, theoretical, and computational approach to the study of molecular processes. This approach stresses the analysis of such processes in terms of dynamical, time-evolving states rather than stationary molecular states. Thus, rovibrational and scattering states are reduced to less prominent roles as is the case in most modem wavepacket treatments of molecular reaction dynamics. Unlike most theoretical methods, END also relegates electronic stationary states, potential energy surfaces, adiabatic and diabatic descriptions, and nonadiabatic coupling terms to the background in favor of a dynamic, time-evolving description of all electrons. 

For example, the ZN theory, which overcomes all the defects of the Landau-Zener-Stueckelberg theory, can be incorporated into various simulation methods in order to clarify the mechanisms of dynamics in realistic molecular systems. Since the nonadiabatic coupling is a vector and thus we can always determine the relevant one-dimensional (ID) direction of the transition in multidimensional space, the 1D ZN theory can be usefully utilized. Furthermore, the comprehension of reaction mechanisms can be deepened, since the formulas are given in simple analytical expressions. Since it is not feasible to treat realistic large systems fully quantum mechanically, it would be appropriate to incorporate the ZN theory into some kind of semiclassical methods. The promising semiclassical methods are (1) the initial value... 

Figure 6. Initial rovibrational state specified reaction probabilities. Solid line exact quantum mechanical numerical solution. Solid line with solid square generalized TSH with use of the nonadiabatic coupling vector. Solid line with open circle generalized TSH with use of Hessian. Sur= 1(2) means the ground (excited) potential energy surface. Taken from Ref. [51].

The aim of this work is to obtain the four lowest E curves and wavefunctions of BH at the same level of accuracy and to bring out the interplay of ionic, Rydberg and valence states at energies and internuclear distances which were not previously investigated. We have therefore made use of a method, already put forward by us [16,17] to determine at once quasi-diabatic and adiabatic states, potential energy cnrves and approximate nonadiabatic couplings. We have analogously determined the first three E+ states, of which only the lowest had been theoretically studied... 

The C matrix, the columns ofwhich, Cj(, are the eigenvectors of H, is normally not too different from the matrix defined above. However, the QDPT treatment, applied either to an adiabatic or to a diabatic zeroth-order basis, is necessary in order to prevent serious artefacts, especially in the case of avoided crossings [27]. The preliminary diabatisation makes it easier to interpolate the matrix elements of the hamiltonian and of other operators as functions of the nuclear coordinates and to calculate the nonadiabatic coupling matrix elements ... 

Fig. 8 shows the nonadiabatic coupling functions between states, which are completely dominated by the double crossing feature. The coefficient mixing term is a very good approximation of the total matrix elements, thus confirming the considerations already put forward with regard to the singlets. 

---