Adiabatic state representation trajectory surface hopping

Quantum chemical methods, exemplified by CASSCF and other MCSCF methods, have now evolved to an extent where it is possible to routinely treat accurately the excited electronic states of molecules containing a number of atoms. Mixed nuclear dynamics, such as swarm of trajectory based surface hopping or Ehrenfest dynamics, or the Gaussian wavepacket based multiple spawning method, use an approximate representation of the nuclear wavepacket based on classical trajectories. They are thus able to use the infoiination from quantum chemistry calculations required for the propagation of the nuclei in the form of forces. These methods seem able to reproduce, at least qualitatively, the dynamics of non-adiabatic systems. Test calculations have now been run using duect dynamics, and these show that even a small number of trajectories is able to produce useful mechanistic infomiation about the photochemistry of a system. In some cases it is even possible to extract some quantitative information. 

The surface-hopping trajectories obtained in the adiabatic representation of the QCLE contain nonadiabatic transitions between potential surfaces including both single adiabatic potential surfaces and the mean of two adiabatic surfaces. This picture is qualitatively different from surface-hopping schemes [2,56] which make the ansatz that classical coordinates follow some trajectory, R(t), while the quantum subsystem wave function, expanded in the adiabatic basis, is evolved according to the time dependent Schrodinger equation. The potential surfaces that the classical trajectories evolve along correspond to one of the adiabatic surfaces used in the expansion of the subsystem wavefunction, while the subsystem evolution is carried out coherently and may develop into linear combinations of these states. In such schemes, the environment does not experience the force associated with the true quantum state of the subsystem and decoherence by the environment is not automatically taken into account. Nonetheless, these methods have provided com-... 

Several algorithms have been constructed to simulate the solution of the QCLE. The simulation methods usually utilize particular representations of the quantum subsystem. Surface-hopping schemes that make use of the adiabatic basis have been constructed density matrix evolution has been carried out in the diabatic basis using trajectory-based methods, some of which make use of a mapping representation of the diabatic states.A representation of the dynamics in the force basis has been implemented to simulate the dynamics using the multithreads algorithm. ... 

---