Improved Ehrenfest Approach to Model Correlated Electron–Nuclear Dynamics

Abstract

Mixed quantum-classical mechanics descriptions are critical to modeling coupled electron-nuclear dynamics, i.e. non-adiabatic molecular dynamics, relevant to photochemical and photophysical processes. We argue that, for polyatomic molecules, such mixed dynamics can not be efficiently described in terms of a matrix gauge potential and develop the concept of a "length gauge" effective Hamiltonian, which helps clarifying certain aspects for the popular non-adiabatic computational approaches. In particular, within such an effective Hamiltonian formalism one readily derives the momentum rescaling boundary condition, used in the surface hopping algorithms. Furthermore, using the new formalism, we introduce a coupled Gaussian wavepacket parameterization of the nuclear wavefunction, which generalizes the Ehrenfest approach to account for electron-nuclei correlations. We test this new approach, Ehrenfest-Plus, on the standard set of model problems that probe electron-nuclear correlation in non-adiabatic transitions. The high accuracy of our approach, combined with mixed-quantum classical efficiency, opens a path for improved simulation of nonadiabatic molecular dynamics in realistic molecular systems. arXiv:1808.05263v2 [physics.chem-ph]

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

“…Thus, when substituting eq into the Schrödinger equation, eq , one should assume that ψ a are explicit functions of t , ψ a = ψ a ( r , R t ). Following such assumption, the substitution of eq into eq gives , with Note that in the expression for V ab the Hamiltonian, Ĥ el , contains R as an operator, while R t is a variable. Expansion of Ĥ el ( R ) in the vicinity of R t leads to an effective Hamiltonian, with In eq , parameters E a , f a , and κ̂ a are the a th expectation values of the electronic Hamiltonian, Ĥ el ( R t ), and its first and second derivatives, respectively.…”

A molecular perspective on Tully models for nonadiabatic dynamics

“…The diversity of strategies tested with some or all of Tully models is astonishing. It includes, among others, full multiple spawning, 15 semiclassical initial value representation, 23 symmetrical quasi-classical windowing combined with Meyer-Miller mapping Hamiltonian, 24 surface hopping Herman-Kluk semiclassical initial value representation method, 25 semiclassical Monte-Carlo, 26 dephasing representation of quantum fidelity, 27 counter-propagating wave methods trajectories, 28 Ehrenfest-Plus, 29 coupled-trajectory mixed quantum/classical dynamics, 30,31 quantum trajectory mean-field approach, 32 iterative linearized approach to nonadiabatic dynamics, 33 mean-field dynamics with stochastic decoherence, 34 nonadiabatic Bohmian dynamics, 35,36 mean-field molecular dynamics with surface hopping, 37 non-Hermitian surface hopping, 38 multi-state trajectory approach to nonadiabatic dynamics, 39…”

Experimental evidence of strong electromagnetic shower enhancement induced by high-energy photons in a thick oriented tungsten crystal

“…The enhancement of radiation and PP by highenergy e ± and photons in strong and sub-strong crystalline fields was first experimentally investigated in the 1980s. The studies on PP were firstly focused on lighter elements such as silicon and germanium [28], and only at a later stage on high-Z crystalline metals such as tungsten and iridium [29,30,31]. These studies were driven by the need to test the feasibility of an intense positron source [32] and to develop a compact photon converter, i.e., a device to separate the photon and neutral hadron beam components with minimal absorption or scattering of the latter [29,30], for the NA48 experiment at CERN; eventually, an iridium crystal with a thickness of 0.98 X 0 was employed [33].…”