BO approximation,

BO approximation vanishes even at low energies, hence the formulation of generalized BO equations become worth while considering. 

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... 

In this section, we prove that the non-adiabatic matiices have to be quantized ( similar to Bohr-Sommerfeld quantization of the angulai momentum) in order to yield a continous, uniquely defined, diabatic potential matrix W(i). In another way, the extended BO approximation will be applied only to those cases that fulfill these quantization rules. The ADT matrix A(s,so) transforms a given adiabatic potential matiix u(i) to a diabatic matiix W(s, so)... 

The ordinary BO approximate equations failed to predict the proper symmetry allowed transitions in the quasi-JT model whereas the extended BO equation either by including a vector potential in the system Hamiltonian or by multiplying a phase factor onto the basis set can reproduce the so-called exact results obtained by the two-surface diabatic calculation. Thus, the calculated hansition probabilities in the quasi-JT model using the extended BO equations clearly demonshate the GP effect. The multiplication of a phase factor with the adiabatic nuclear wave function is an approximate treatment when the position of the conical intersection does not coincide with the origin of the coordinate axis, as shown by the results of [60]. Moreover, even if the total energy of the system is far below the conical intersection point, transition probabilities in the JT model clearly indicate the importance of the extended BO equation and its necessity. 

Knowledge of the underlying nuclear dynamics is essential for the classification and description of photochemical processes. For the study of complicated systems, molecular dynamics (MD) simulations are an essential tool, providing information on the channels open for decay or relaxation, the relative populations of these channels, and the timescales of system evolution. Simulations are particularly important in cases where the Bom-Oppenheimer (BO) approximation breaks down, and a system is able to evolve non-adiabatically, that is, in more than one electronic state. 

Central to the description of this dynamics is the BO approximation. This separates the nuclear and electionic motion, and allows the system evolution to be described by a function of the nuclei, known as a wavepacket, moving over a PES provided by the (adiabatic) motion of the electrons. 

In this section, the basic theory of molecular dynamics is presented. Starting from the BO approximation to the nuclear Schrddinger equation, the picture of nuclear dynamics is that of an evolving wavepacket. As this picture may be unusual to readers used to thinking about nuclei as classical particles, a few prototypical examples are shown. 

Using the BO approximation, the Schrddinger equation describing the time evolution of the nuclear wave function, can be written... 

The adiabatic picture developed above, based on the BO approximation, is basic to our understanding of much of chemistry and molecular physics. For example, in spectroscopy the adiabatic picture is one of well-defined spectral bands, one for each electronic state. The smicture of each band is then due to the shape of the molecule and the nuclear motions allowed by the potential surface. This is in general what is seen in absorption and photoelectron spectroscopy. There are, however, occasions when the picture breaks down, and non-adiabatic effects must be included to give a faithful description of a molecular system [160-163]. 

In Section II, molecular dynamics within the BO approximation was introduced. As shown in Appendix A, the full nuclear Schrodinger equation is, however. 

The familiar BO approximation is obtained by ignoring the operators A completely. This results in the picture of the nuclei moving over the PES provided by the electrons, which are moving so as to instantaneously follow the nuclear motion. Another common level of approximation is to exclude the off-diagonal elements of this operator matrix. This is known as the Bom-Huang, or simply the adiabatic, approximation (see [250] for further details of the possible approximations and nomenclature associated with the nuclear Schrodinger equation). 

As described above in Appendix A, within the BO approximation the nuclear Schrddinger equation is... 

Within the Bom-Oppenheimer (BO) approximation, A) and B) may be written as the product of an electronic wave function, M)gj and a nuclear wave function M) . 

An alternative approximation scheme, also proposed by Bom and Oppenheimer [5-7], employed the straightforward perturbation method. To tell the difference between these two different BO approximation, we call the latter the crude BOA (CBOA). A main purpose of this chapter is to study the original BO approximation, which is often referred to as the crude BO approximation and to develop this approximation into a practical method for computing potential energy suifaces of molecules. 

Nevertheless, the examination of the applicability of the crude BO approximation can start now because we have worked out basic methods to compute the matrix elements. With the advances in the capacity of computers, the test of these methods can be done in lower and lower cost. In this work, we have obtained the formulas and shown their applications for the simple cases, but workers interested in using these matrix elements in their work would find that it is not difficult to extend our results to higher order derivatives of Coulomb interaction, or the cases of more-than-two-atom molecules. 

Until now we have implicitly assumed that our problem is formulated in a space-fixed coordinate system. However, electronic wave functions are naturally expressed in the system bound to the molecule otherwise they generally also depend on the rotational coordinate 4>. (This is not the case for E electronic states, for which the wave functions are invariant with respect to (j> ) The eigenfunctions of the electronic Hamiltonian, v / and v , computed in the framework of the BO approximation ( adiabatic electronic wave functions) for two electronic states into which a spatially degenerate state of linear molecule splits upon bending. 

Let us examine a special but more practical case where the total molecular Hamiltonian, H, can be separated to an electronic part, W,.(r,s Ro), as is the case in the usual BO approximation. Consequendy, the total molecular wave function fl(R, i,r,s) is given by the product of a nuclear wave function X uc(R, i) and an electronic wave function v / (r, s Ro). We may then talk separately about the permutational properties of the subsystem consisting of electrons, and the subsystemfs) formed of identical nuclei. Thus, the following commutative laws Pe,Hg =0 and =0 must be satisfied X =... 

---