J. Phys. Soc. Jpn. Downloaded from journals.jps.jp by Univ of Southern California on 11/07/19
Journal of the Physical Society of Japan 84, 094002 (2015) http://dx.doi.org/10.7566/JPSJ.84.094002
Analysis of the Trajectory Surface Hopping Method from the Markov State Model Perspective
Alexey V. Akimov1, Dhara Trivedi2, Linjun Wang1, and Oleg V. Prezhdo1+
1Department of Chemistry, University of Southern California, Los Angeles, CA 90089, U.S.A. 2Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627, U.S.A. (Received May 20, 2015; accepted July 22, 2015; published online August 18, 2015) We analyze the applicability of the seminal fewest switches surface hopping (FSSH) method of Tully to modeling quantum transitions between electronic states that are not coupled directly, in the processes such as Auger recombination. We address the known deficiency of the method to describe such transitions by introducing an alternative definition for the surface hopping probabilities, as derived from the Markov state model perspective. We show that the resulting transition probabilities simplify to the quantum state populations derived from the time-dependent Schrödinger equation, reducing to the rapidly switching surface hopping approach of Tully and Preston. The resulting surface hopping scheme is simple and appeals to the fundamentals of quantum mechanics. The computational approach is similar to the FSSH method of Tully, yet it leads to a notably different performance. We demonstrate that the method is particularly accurate when applied to superexchange modeling. We further show improved accuracy of the method, when applied to one of the standard test problems. Finally, we adapt the derived scheme to atomistic simulation, combine it with the time-domain density functional theory, and show that it provides the Auger energy transfer timescales which are in good agreement with experiment, significantly improving upon other considered techniques.
estimated in FSSH theory because the transitions through 1. Introduction such intermediate states are inhibited by the probabilities of Theory of molecular dynamics with quantum transitions is overcoming the high-energy barrier. Simulations of non- an actively developing field of computational chemistry.1–14) adiabatic processes using FSSH are usually performed in the The interest in this subject is constantly stimulated by the adiabatic basis, in which FSSH shows reasonable accuracy. need for accurate and computationally tractable simulation However, the use of diabatic states (e.g., donor and acceptor methodologies for modeling processes involving coupled states) may be advantageous, both in terms of computational electron–nuclear dynamics in contemporary materials. In efforts and for interpretation of the problem physics. Diabatic particular, such techniques are required for studying photo- states are used regularly in phenomenological models of induced charge transfer in photovoltaic and photocatalytic charge and energy transfer because they have well-defined materials,15–20) long-distance energy transfer in biological physical meaning. Superexchange effects are traditionally photosynthetic complexes, or photo- and electro-stimulated described in a diabatic representation. An adiabatic picture mechanical response in nanoscale systems.20,21) can lead to spurious effects arising from artificial delocaliza- A great body of methods for modeling molecular dynamics tion of an adiabatic state between donor and acceptor sites with quantum transitions has been developed over the last that are very distant from each other. Finally, the adiabatic several decades. We refer the reader to specialized reviews representation leads to numerical problems with calculation discussing some of these methods.22–27) Among the variety of non-adiabatic coupling.45,46) Under the conditions when of techniques, Tully’s fewest switches surface hopping diabatic states are strongly localized and many pairs of states (FSSH)1) remains one of the most utilized methods, due to are decoupled, the accuracy of the original FSSH method is its conceptual and practical simplicity, good accuracy, and often unacceptable when diabatic states are used, motivating high computational efficiency. The FSSH method is not a development of the alternative formulations.28,47) panacea, and there are examples where it breaks down. In Recently, the problem of superexchange has been particular, the superexchange problem28) requires a more addressed by Wang et al.28) who proposed a global flux advanced treatment. surface hopping (GFSH) algorithm. In it, hopping between The need for consideration of superexchange effects is a pair of states is determined by a relative change of state motivated in applications by excitonic, many-particle proc- population during the infinitesimal time slice. The method esses, such as Auger recombination and energy transfer,29–35) avoids computing the hopping probabilities based on one- singlet fission36–39) or Raman scattering. The superexchange particle properties (such as non-adiabatic couplings between mechanism is also important for tunneling and conduc- one-particle states), thus solving the superexchange problem tivity in molecular electronics.40–44) In all these processes, phenomenologically. In the alternative second-quantized transitions between excitonic states can involve more than surface hopping (SQUASH) approach,47) the superexchange one simultaneous single-particle transition. A direct consid- problem was solved by extending the quantum dynamics to eration of the dynamics of such excitations is prohibited by the space of second-quantized states representing coupled the Slater rules, and only sequential mechanisms are viable. trajectories. The method introduces many-particle states, The problem of superexchange arises when there is at least allowing a single transition between a pair of such states to one pair of states that are not coupled to each other directly. capture a many-body transition. The states can be either diabatic or adiabatic. The transition Although the developed techniques provide notable between these states is mediated by an intermediate high- improvement of the accuracy in modeling the superexchange lying state. The overall transition rates are strongly under- problem, they have certain limitations. Specifically, the 094002-1 ©2015 The Physical Society of Japan J. Phys. Soc. Jpn. Downloaded from journals.jps.jp by Univ of Southern California on 11/07/19
J. Phys. Soc. Jpn. 84, 094002 (2015) A. V. Akimov et al.
GFSH approach relies on an ad hoc scheme for splitting total MSSH method to the classical path approximation (CPA), probability to leave a state i into probabilities of doing this which neglects electron back-reaction on nuclei, and apply via different channels i ! j. As it is pointed out by Wang it to compute Auger-type energy transfer dynamics in et al.,28) there are, in principle, many possibilities of such atomistic quantum dots (QDs). We conclude that the method splitting. Different splitting schemes define different hopping can be useful for a number of systems and processes, where probabilities, which may affect the results of surface hopping the standard FSSH may have problems. We also discuss simulations. Thus, the question of uniqueness, physical limitations of the proposed approach. meaning and rigor of the definition of surface hopping probabilities remains open, presenting a conceptual chal- 2. Theory lenge. On the other hand, SQUASH47) introduces coupled 2.1 Analysis of the FSSH trajectories that may exchange energy. The number of N- One of the basic elements of Tully’s FSSH method is the particle states describing N coupled trajectories increases definition of hopping probabilities. The method is derived the dimensionality of the problem very rapidly, making from the following assumption of hopping out of a single the computations costly. The method also relies on the state. Assume that the number of trajectories in the current conceptually challenging prescription of the energy exchange state i at time t is NiðtÞ. If at the later time t þ dt this number in this ensemble-based approach to quantum dynamics. decreases to Niðt þ dtÞ, with Niðt þ dtÞ