Schrodinger equation body-fixed Hamiltonian

This work introduced the concept of a vibronic R-matrix, defined on a hypersurface in the joint coordinate space of electrons and intemuclear coordinates. In considering the vibronic problem, it is assumed that a matrix representation of the Schrodinger equation for N+1 electrons has been partitioned to produce an equivalent set of multichannel one-electron equations coupled by a matrix array of nonlocal optical potential operators [270], In the body-fixed reference frame, partial wave functions in the separate channels have the form p(q xN)YL(0, radial channel orbital function i/(q r) and antisymmetrized in the electronic coordinates. Here 0 is a fixed-nuclei A-electron target state or pseudostate and Y] is a spherical harmonic function. Both and i r are parametric functions of the intemuclear coordinate q. It is assumed that the target states 0 for each value of q diagonalize the A-electron Hamiltonian matrix and are orthonormal. 

Equation (1-13) or its body-fixed equivalent is of little use for Van der Waals complexes, as it discriminates one nuclear coordinate, e.g. y = 1. Specific mathematical forms of Hamiltonians describing the nuclear motions in Van der Waals dimers have been developed (7). This point will be discussed in more details in Section 12.4. Here we only want to stress that whatever the mathematical form of the Hamiltonian is used to solve the problem of nuclear motions, the results will be the same, if the Schrodinger equation is solved exactly. However, in weakly bound complexes there is a hierarchy of motions due to the strong intramolecular forces which determine the internal vibrations of the molecules, and to much weaker intermolecular forces which determine their relative translations and rotations. This hierarchy allows to make a separation between the intramolecular vibrations with high frequencies and the intermolecular modes with much lower frequencies. Such a separation of the fast intramolecular vibrations and slow rotation-vibration-tunneling motions can be performed if a suitable form of the Hamiltonian for the nuclear motions in Van der Waals molecules is used. 

For a complete treatment of a laser-driven molecule, one must solve the many-body, multidimensional time-dependent Schrodinger equation (TDSE). This represents a tremendous task and direct wavepacket simulations of nuclear and electronic motions under an intense laser pulse is presently restricted to a few bodies (at most three or four) and/or to a model of low dimensionality [27]. For a more general treatment, an approximate separation of variables between electrons (fast subsystem) and nuclei (slow subsystem) is customarily made, in the spirit of the BO approximation. To lay out the ideas underlying this approximation as adapted to field-driven molecular dynamics, we will consider from now on a molecule consisting of Nn nuclei (labeled a, p,...) and Ne electrons (labeled /, j,...), with position vectors Ro, and r respectively, defined in the center of mass (rotating) body-fixed coordinate system, in a classical field E(f) of the form Eof t) cos cot). The full semiclassical length gauge Hamiltonian is written, for a system of electrons and nuclei, as [4]... 

---