Title: Classical and nonclassical effects in surface hopping methodology for simulating coupled electronic-nuclear dynamics
In this paper, we analyze the detailed quantum-classical behavior of two alternative approaches to simulating molecular dynamics with electronic transitions: the popular fewest switches surface hopping (FSSH) method, introduced by Tully in 1990 [Tully, \textit{J.~Chem.~Phys.}, 1990, \textbf{93}, 1061] and our recently developed quantum trajectory surface hopping (QTSH) method [Martens, \textit{J.~Phys.~Chem.~A}, 2019 \textbf{123}, 1110]. Both approaches employ an independent ensemble of trajectories that undergo stochastic transitions between electronic surfaces. The methods differ in their treatment of energy conservation, with FSSH imposing conservation of the classical kinetic plus potential energy by rescaling of the classical momenta when a surface hop occurs, while QTSH incorporates quantum forces throughout the dynamics which lead naturally to the conservation of the full quantum-classical energy. We investigate the population transfer and energy budget of the surface hopping methods for several simple model systems and compare with exact quantum result. In addition, the detailed dynamics of the trajectory ensembles in phase space are compared with the quantum evolution in the Wigner representation. Conclusions are drawn. more »« less
Tao, Zhen; Bian, Xuezhi; Wu, Yanze; Rawlinson, Jonathan; Littlejohn, Robert G; Subotnik, Joseph E
(, The Journal of Chemical Physics)
We show that standard Ehrenfest dynamics does not conserve linear and angular momentum when using a basis of truncated adiabatic states. However, we also show that previously proposed effective Ehrenfest equations of motion [M. Amano and K. Takatsuka, “Quantum fluctuation of electronic wave-packet dynamics coupled with classical nuclear motions,” J. Chem. Phys. 122, 084113 (2005) and V. Krishna, “Path integral formulation for quantum nonadiabatic dynamics and the mixed quantum classical limit,” J. Chem. Phys. 126, 134107 (2007)] involving the non-Abelian Berry force do maintain momentum conservation. As a numerical example, we investigate the Kramers doublet of the methoxy radical using generalized Hartree–Fock with spin–orbit coupling and confirm that angular momentum is conserved with the proper equations of motion. Our work makes clear some of the limitations of the Born–Oppenheimer approximation when using ab initio electronic structure theory to treat systems with unpaired electronic spin degrees of freedom, and we demonstrate that Ehrenfest dynamics can offer much improved, qualitatively correct results.
Villaseco_Arribas, Evaristo; Maitra, Neepa T; Agostini, Federica
(, The Journal of Chemical Physics)
Advances in coherent light sources and development of pump–probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron–nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. This work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.
Griffin, Hadassah B; Kryjevski, Andrei B; Kilina, Svetlana V; Kilin, Dmitri S
(, ECS Meeting Abstracts)
The work provides computational arguments in support of excitonic approach for the treatment of the photo-induced processes in semiconductor quantum dots. The non-radiative relaxation, non-radiative recombination, and photo-luminescence quantum yield are computed for a range of atomistic models of semiconductor quantum dots (QDs) in the quantum confinement regime. The excitonic (EX) approach is compared to independent orbital approximation (IOA) approach. Both approaches address dissipation of the electronic energy from electronic degrees of freedom to thermal vibrations of the lattice. The difference of two approaches appears in treatment of energies of electronic states and in a way how the electron-phonon interaction is taken into account. IOA approach uses energies of Kohn-Sham orbitals and on the fly non-adiabatic couplings. [1-3] EX approach uses Bethe-Salpeter equation (BSE) for energies.[4-6] The excitonic wavefunctions from BSE is used to construct a linear transformation matrix that transforms IOA-based non-adiabatic couplings into an excitonic basis. Both approaches are compared in application to untrasmall 1 nm diameter Si QD.Results include an evidence that hot excitons relax sooner in the excitonic picture than in the IOA picture. The observed effect is rationalized via smaller subgaps and different available relaxation pathways in the excitonic picture. The most surprising result is found for the simulated emission spectrum. The spectum in the excitonic picture demonstrates intensity in several 5 orders of magnitude higher than in the IOA picture. This observation is related to formation of a bright exciton in the lowest excitation of the ultra-small Si QDs. Obtained evidence favors excitonic approach and promises a reliable interpretation and prediction of time-dependent observables in a range of semiconductor quantum dots of different composition, sizes, and surface environment.[7] Most intriguing results are expected for QDs representing interface between PbSe and CdSe. [8] Support of National Science foundation via NSF CHE-2004197 is gratefully acknowledged. [1] D. S. Kilin and D. A. Micha, “Relaxation of photoexcited electrons at a nanostructured Si(111) surface”, J. Phys. Chem. Lett. 1, 1073-1077 (2010). [2] D. J. Vogel and D. S. Kilin, "First-Principles Treatment of Photoluminescence in Semiconductors" J. Phys. Chem. C 119, 50, 27954–27964 (2015). [3] D. J. Vogel, A. B. Kryjevski, T. M. Inerbaev, and D. S. Kilin, "Photoinduced Single- and Multiple-Electron Dynamics Processes Enhanced by Quantum Confinement in Lead Halide Perovskite Quantum Dots", J. Phys. Chem. Lett. 8, 13, 3032–3039(2017). [4] A. B. Kryjevski and Dmitri Kilin, "Enhanced multiple exciton generation in amorphous silicon nanowires and films", Molec. Phys. 114, 365-379 (2016). [5] M. Rohlfing and S. G. Louie, "Electron-Hole Excitations in Semiconductors and Insulators", Phys. Rev. Lett. 81, 2312-2315 (1998). [6] T. Sander, G. Kresse, "Macroscopic dielectric function within time-dependent density functional theory—Real time evolution versus the Casida approach", J. Chem. Phys. 146, 064110 (2017). [7] S. V. Kilina, P. K. Tamukong, and D. S. Kilin, "Surface Chemistry of Semiconducting Quantum Dots: Theoretical Perspectives", Acc. Chem. Res. 49, 10, 2127–2135 (2016). [8] H. B. Griffin, A. B. Kryjevski, and Dmitri S. Kilin, "Ab initio calculations of through-space and through-bond charge-transfer properties of interacting Janus-like PbSe and CdSe quantum dot heterostructures", Molec. Phys., e2273415 (2023).
Abstract In this work, we report the development and assessment of the nonadiabatic molecular dynamics approach with the electronic structure calculations based on the linearly scaling subsystem density functional method. The approach is implemented in an open-source embedded Quantum Espresso/Libra software specially designed for nonadiabatic dynamics simulations in extended systems. As proof of the applicability of this method to large condensed-matter systems, we examine the dynamics of nonradiative relaxation of excess excitation energy in pentacene crystals with the simulation supercells containing more than 600 atoms. We find that increased structural disorder observed in larger supercell models induces larger nonadiabatic couplings of electronic states and accelerates the relaxation dynamics of excited states. We conduct a comparative analysis of several quantum-classical trajectory surface hopping schemes, including two new methods proposed in this work (revised decoherence-induced surface hopping and instantaneous decoherence at frustrated hops). Most of the tested schemes suggest fast energy relaxation occurring with the timescales in the 0.7–2.0 ps range, but they significantly overestimate the ground state recovery rates. Only the modified simplified decay of mixing approach yields a notably slower relaxation timescales of 8–14 ps, with a significantly inhibited ground state recovery.
Hu, Deping; Chng, Benjamin_X K; Ying, Wenxiang; Huo, Pengfei
(, The Journal of Chemical Physics)
We benchmark the accuracy of various trajectory-based non-adiabatic methods in simulating the polariton relaxation dynamics under the collective coupling regime. The Holstein–Tavis–Cummings Hamiltonian is used to describe the hybrid light–matter system of N molecules coupled to a single cavity mode. We apply various recently developed trajectory-based methods to simulate the population relaxation dynamics by initially exciting the upper polariton state and benchmark the results against populations computed from exact quantum dynamical propagation using the hierarchical equations of motion approach. In these benchmarks, we have systematically varied the number of molecules N, light–matter detunings, and the light–matter coupling strengths. Our results demonstrate that the symmetrical quasi-classical method with γ correction and spin-mapping linearized semi-classical approaches yield more accurate polariton population dynamics than traditional mixed quantum-classical methods, such as the Ehrenfest and surface hopping techniques.
Martens, Craig C. Classical and nonclassical effects in surface hopping methodology for simulating coupled electronic-nuclear dynamics. Retrieved from https://par.nsf.gov/biblio/10106649. Faraday Discussions . Web. doi:10.1039/C9FD00042A.
Martens, Craig C. Classical and nonclassical effects in surface hopping methodology for simulating coupled electronic-nuclear dynamics. Faraday Discussions, (). Retrieved from https://par.nsf.gov/biblio/10106649. https://doi.org/10.1039/C9FD00042A
@article{osti_10106649,
place = {Country unknown/Code not available},
title = {Classical and nonclassical effects in surface hopping methodology for simulating coupled electronic-nuclear dynamics},
url = {https://par.nsf.gov/biblio/10106649},
DOI = {10.1039/C9FD00042A},
abstractNote = {In this paper, we analyze the detailed quantum-classical behavior of two alternative approaches to simulating molecular dynamics with electronic transitions: the popular fewest switches surface hopping (FSSH) method, introduced by Tully in 1990 [Tully, \textit{J.~Chem.~Phys.}, 1990, \textbf{93}, 1061] and our recently developed quantum trajectory surface hopping (QTSH) method [Martens, \textit{J.~Phys.~Chem.~A}, 2019 \textbf{123}, 1110]. Both approaches employ an independent ensemble of trajectories that undergo stochastic transitions between electronic surfaces. The methods differ in their treatment of energy conservation, with FSSH imposing conservation of the classical kinetic plus potential energy by rescaling of the classical momenta when a surface hop occurs, while QTSH incorporates quantum forces throughout the dynamics which lead naturally to the conservation of the full quantum-classical energy. We investigate the population transfer and energy budget of the surface hopping methods for several simple model systems and compare with exact quantum result. In addition, the detailed dynamics of the trajectory ensembles in phase space are compared with the quantum evolution in the Wigner representation. Conclusions are drawn.},
journal = {Faraday Discussions},
author = {Martens, Craig C.},
}
Warning: Leaving National Science Foundation Website
You are now leaving the National Science Foundation website to go to a non-government website.
Website:
NSF takes no responsibility for and exercises no control over the views expressed or the accuracy of
the information contained on this site. Also be aware that NSF's privacy policy does not apply to this site.
An official website of the United States government
