Mean-field trajectory method results

Figure 3. Time-dependent simulations of the nonadiabatic photoisomerization dynamics exhibited by Model III, comparing results of the mean-field-trajectory method (dashed lines), the surface-hopping approach (thin lines), and exact quanmm 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 Pdsit) that the sytem remains in the initially prepared cis conformation.

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. 

Figure 19. Time-dependent (a) diabatic and (b) adiabatic electronic excited-state populations and (c) vibrational mean positions as obtained for Model 1. Shown are results of the mean-field trajectory method (dotted lines), the quasi-classical mapping approach (thin full lines), and exact quantum calculations (thick full lines).

PESs, and (iv) obey the principle of microreversibility. Section 3 describes in some detail the mean-field trajectory method and also discusses its connection to time-dependent self-consistent-field schemes. The surface-hopping method is considered in Sec. 4, which discusses various motivations of the ansatz as well as several variants of the implementation. Section 5 gives a brief account on the quantum-classical Liouville description and considers the possibility of an exact stochastic realization of its equation of motion. The mapping formalism, its relation to other formulations, and its quasi-classical implementation is introduced in Sec. 6. Section 7 is concerned with the semiclassical description of nonadiabatic quantum mechanics. Section 8 summarizes our results and concludes with some general remarks. 

Finally, we consider Model III, which describes an ultrafast photoin-duced isomerization process. Figure 13 shows quantum mechanical results as well as results of the ZPE-corrected mapping approach for three different observables the adiabatic and diabatic population of the excited state, and Pcis> the probability that the system remains in the initially prepared cis-conformation. A ZPE correction of 7 = 0.5 has been used in the mapping calculation, based on the criterion to reproduce the quantum mechanical long-time limit of the adiabatic population. It is seen that the ZPE-corrected mapping approach represents an improvement compared to the mean-field trajectory method (cf. Fig. 3), in particular for the adiabatic population. The influence of the ZPE correction on the dynamics of the observables is, however, not as large as in the two other models. 

Finally, we consider Model V by describing two examples of outer-sphere electron-transfer in solution. Figures 7 and 8 display results for the diabatic electronic population for Models Va and Vb, respectively. Similar to the mean-field trajectory calculations, for Model Va the SH results are in excellent agreement with the quantum calculations, while for Model Vb the SH method only is able to describe the short-time dynamics. As for the three-mode Model IVb discussed above, the SH calculations in particular predict an incorrect long-time limit for the diabatic population. The origin of this problem will be discussed in more detail in Section VI in the context of the mapping formulation. 

To summarize, it has been found that the SH method is able to at least qualitatively describe the complex photoinduced electronic and vibrational relaxation dynamics exhibited by the model problems under consideration. The overall quality of SH calculations is typically somewhat better than the quality of the mean-field trajectory results. In particular, this holds in the case of several curve crossings (see Fig. 2) as well as when the dynamics and the observables of interest are essentially of adiabatic nature— for example, for the calculation of the adiabatic population dynamics associated with a conical intersection (see Figs. 3 and 12). Furthermore, we have briefly discussed various consistency problems of a simple quasi-classical SH description. It has been shown that binned electronic population probabilities and no momentum adjustment for classically forbidden transitions help us to improve this matter. There have been numerous suggestions to further improve the hopping algorithm [70-74] however, the performance of all these variants seems to depend largely on the problem under consideration. 

---