No Arabic abstract
We investigate two approaches to derive the proper Floquet-based quantum-classical Liouville equation (F-QCLE) for laser-driven electron-nuclear dynamics. The first approach projects the operator form of the standard QCLE onto the diabatic Floquet basis, then transforms to the adiabatic representation. The second approach directly projects the QCLE onto the Floquet adiabatic basis. Both approaches yield a form which is similar to the usual QCLE with two modifications: 1. The electronic degrees of freedom are expanded to infinite dimension. 2. The nuclear motion follows Floquet quasi-energy surfaces. However, the second approach includes an additional cross derivative force due to the dual dependence on time and nuclear motion of the Floquet adiabatic states. Our analysis and numerical tests indicate that this cross derivative force is a factitious artifact, suggesting that one cannot safely exchange the order of Floquet state projection with adiabatic transformation. Our results are in accord with similar findings by Izmaylov et al., who found that transforming to the adiabatic representation must always be the last operation applied, though now we have extended this result to a time-dependent Hamiltonian. This paper and the proper derivation of the F-QCLE should lay the basis for further improvements of Floquet surface hopping.
Although the quantum classical Liouville equation (QCLE) arises by cutting off the exact equation of motion for a coupled nuclear-electronic system at order 1 (1 = $hbar^0$ ), we show that the QCLE does include Berrys phase effects and Berrys forces (which are proportional to a higher order, $hbar$ = $hbar^1$ ). Thus, the fundamental equation underlying mixed quantum-classical dynamics does not need a correction for Berrys phase effects and is valid for the case of complex Hamiltonians. Furthermore, we also show that, even though Tullys surface hopping model ignores Berrys phase, Berrys phase effects are included automatically within Ehrenfest dynamics. These findings should be of great importance if we seek to model coupled nuclear-electronic dynamics for systems with spin-orbit coupling, where the complex nature of the Hamiltonian is paramount.
Trajectory surface hopping (TSH) is one of the most widely used quantum-classical algorithms for nonadiabatic molecular dynamics. Despite its empirical effectiveness and popularity, a rigorous derivation of TSH as the classical limit of a combined quantum electron-nuclear dynamics is still missing. In this work we aim to elucidate the theoretical basis for the widely used hopping rules. Naturally, we concentrate thereby on the formal aspects of the TSH. Using a Gaussian wave packet limit, we derive the transition rates governing the hopping process at a simple avoided level crossing. In this derivation, which gives insight into the physics underlying the hopping process, some essential features of the standard TSH algorithm are retrieved, namely i) non-zero electronic transition rate (hopping probability) at avoided crossings; ii) rescaling of the nuclear velocities to conserve total energy; iii) electronic transition rates linear in the nonadiabatic coupling vectors. The well-known Landau-Zener model is then used for illustration.
We employ various quantum-mechanical approaches for studying the impact of electric fields on both nonretarded and retarded noncovalent interactions between atoms or molecules. To this end, we apply perturbative and non-perturbative methods within the frameworks of quantum mechanics (QM) as well as quantum electrodynamics (QED). In addition, to provide a transparent physical picture of the different types of resulting interactions, we employ a stochastic electrodynamic approach based on the zero-point fluctuating field. Atomic response properties are described via harmonic Drude oscillators - an efficient model system that permits an analytical solution and has been convincingly shown to yield accurate results when modeling non-retarded intermolecular interactions. The obtained intermolecular energy contributions are classified as field-induced (FI) electrostatics, FI polarization, and dispersion interactions. The interplay between these three types of interactions enables the manipulation of molecular dimer conformations by applying transversal or longitudinal electric fields along the intermolecular axis. Our framework combining four complementary theoretical approaches paves the way toward a systematic description and improved understanding of molecular interactions when molecules are subject to both external and vacuum fields.
Hydrogen diffusion on metals exhibits rich quantum behavior, which is not yet fully understood. Using simulations, we show that many hydrogen diffusion barriers can be categorized into those with parabolic-tops and those with broad-tops. With parabolic-top barriers, hydrogen diffusion evolves gradually from classical hopping to shallow tunneling to deep tunneling as the temperature decreases, and noticeable quantum effects persist at moderate temperatures. In contrast, with broad-top barriers quantum effects become important only at low temperatures and the classical to quantum transition is sharp, at which classical hopping and deep tunneling both occur. This coexistence indicates that more than one mechanism contributes to the quantum reaction rate. The conventional definition of the classical to quantum crossover temperature is invalid for the broad-tops, and we give a new definition. Extending this we propose a model to predict the transition temperature for broad-top diffusion, providing a general guide for theory and experiment.
We present a preliminary extension of the fewest switches surface hopping (FSSH) algorithm to the case of complex Hamiltonians as appropriate for modeling the dynamics of photoexcited molecules in magnetic fields. We make ansatze for the direction of momentum rescaling and we account for Berrys phase effects through magnetic forces as applicable in the adiabatic limit. Because Berrys phase is a nonlocal, topological characteristic of a set of entangled potential energy surfaces, we find that Tullys local FSSH algorithm can only partially capture the correct physics.