Mean-field trajectory method nonadiabatic quantum dynamics

In the MQC mean-field trajectory scheme introduced above, all nuclear DoF are treated classically while a quantum mechanical description is retained only for the electronic DoF. This separation is used in most implementations of the mean-field trajectory method for electronically nonadiabatic dynamics. Another possibility to separate classical and quantum DoF is to include (in addition to the electronic DoF) some of the nuclear degrees of freedom (e.g., high frequency modes) into the quantum part of the calculation. This way, typically, an improved approximation of the overall dynamics can be obtained—albeit at a higher numerical cost. This idea is the basis of the recently proposed self-consistent hybrid method [201, 202], where the separation between classical and quantum DoF is systematically varied to improve the result for the overall quantum dynamics. For systems in the condensed phase with many nuclear DoF and a relatively smooth distribution of the electronic-vibrational coupling strength (e.g.. Model V), the separation between classical and quanmm can, in fact, be optimized to obtain numerically converged results for the overall quantum dynamics [202, 203]. 

The mapping approach outlined above has been designed to furnish a well-defined classical limit of nonadiabatic quantum dynamics. The formalism applies in the same way at the quantum-mechanical, semiclassical (see Section VIII), and quasiclassical level, respectively. Most important, no additional assumptions but the standard semiclassical and quasi-classical approximations are needed to get from one level to another. Most of the established mixed quantum-classical methods such as the mean-field-trajectory method or the surface-hopping approach do invoke additional assumptions. The comparison of the mapping approach to these formulations may therefore (i) provide insight into the nature of these additional approximation and (ii) indicate whether the conceptual virtues of the mapping approach may be expected to result in practical advantages. 

Fig. 3. Time-dependent simulations of the nonadiabatic photoisomerization dynamics exhibited by Model 111, comparing results of the mean-field-trajectory method (dashed lines), the siu-face-hopping approach (thin hues), and exact quantum calculations (full lines). Shown are the population probabilities of the initially prepared (a) adiabatic and (b) diabatic electronic state, respectively, as well as (c) the probability Pcis(f) that the system remains in the initially prepared cis conformation.

---