We present a modified approach for simulating electronically nonadiabatic dynamics based on the Nakajima-Zwanzig generalized quantum master equation (GQME). The modified approach utilizes the fact that the Nakajima-Zwanzig formalism does not require casting the overall Hamiltonian in system-bath form, which is arguably neither natural nor convenient in the case of the Hamiltonian that governs nonadiabatic dynamics. Within the modified approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density operator is fully captured by a memory kernel super-operator. A methodology for calculating the memory kernel from projection-free inputs is developed. Simulating the electronic dynamics via the modified approach, with a memory kernel obtained using exact or approximate methods, can be more cost effective and/or lead to more accurate results than direct application of those methods. The modified approach is compared to previously proposed GQME-based approaches, and its robustness and accuracy are demonstrated on a benchmark spin-boson model with a memory kernel which is calculated within the Ehrenfest method.
We describe a general-purpose framework for formulating the dynamics of any subset of electronic reduced density matrix elements in terms of a formally exact generalized quantum master equation (GQME). Within this framework, the effect of coupling to the nuclear degrees of freedom, as well as to any projected-out electronic reduced density matrix elements, is captured by a memory kernel and an inhomogeneous term, whose dimensionalities are dictated by the number of electronic reduced density matrix elements included in the subset of interest. We show that the memory kernel and inhomogeneous term within such GQMEs can be calculated from projection-free inputs of the same dimensionality, which can be cast in terms of the corresponding subsets of overall system two-time correlation functions. The applicability and feasibility of such reduced-dimensionality GQMEs is demonstrated on the two-state spin-boson benchmark model. To this end, we compare and contrast the following four types of GQMEs: (1) a full density matrix GQME, (2) a single-population scalar GQME, (3) a populations-only GQME, and (4) a subset GQME for any combination of populations and coherences. Using a method based on the mapping Hamiltonian approach and linearized semiclassical approximation to calculate the projection-free inputs, we find that while single-population GQMEs and subset GQMEs containing only one population are less accurate, they can still produce reasonable results and that the accuracy of the results obtained via the populations-only GQME and a subset GQME containing both populations is comparable to that obtained via the full density matrix GQMEs.
In this work, we investigate the ability of different quasiclassical mapping Hamiltonian methods to simulate the dynamics of electronic transitions through conical intersections. The analysis is carried out within the framework of the linear vibronic coupling (LVC) model. The methods compared are the Ehrenfest method, the symmetrical quasiclassical method, and several variations of the linearized semiclassical (LSC) method, including ones that are based on the recently introduced modified representation of the identity operator. The accuracy of the various methods is tested by comparing their predictions to quantum-mechanically exact results obtained via the multiconfiguration time-dependent Hartree (MCTDH) method. The LVC model is found to be a nontrivial benchmark model that can differentiate between different approximate methods based on their accuracy better than previously used benchmark models. In the three systems studied, two of the LSC methods are found to provide the most accurate description of electronic transitions through conical intersections.
The generalized quantum master equation (GQME) provides a powerful framework for simulating electronically nonadiabatic molecular dynamics. Within this framework, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density matrix is fully captured by a memory kernel superoperator. In this paper, we consider two different procedures for calculating the memory kernel of the GQME from projection-free inputs obtained via the combination of the mapping Hamiltonian (MH) approach and the linearized semiclassical (LSC) approximation. The accuracy and feasibility of the two procedures are demonstrated on the spin-boson model. We find that although simulating the electronic dynamics by direct application of the two LSC-based procedures leads to qualitatively different results that become increasingly less accurate with increasing time, restricting their use to calculating the memory kernel leads to an accurate description of the electronic dynamics. Comparison with a previously proposed procedure for calculating the memory kernel via the Ehrenfest method reveals that MH/LSC methods produce memory kernels that are better behaved at long times and lead to more accurate electronic dynamics.
The generalized quantum master equation (GQME) provides a general and formally exact framework for simulating the reduced dynamics of open quantum systems. The recently introduced modified approach to the GQME (M-GQME) corresponds to a specific implementation of the GQME that is geared toward simulating the dynamics of the electronic reduced density matrix in systems governed by an excitonic Hamiltonian. Such a Hamiltonian, which is often used for describing energy and charge transfer dynamics in complex molecular systems, is given in terms of diabatic electronic states that are coupled to each other and correspond to different nuclear Hamiltonians. Within the M-GQME approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic density matrix is fully captured by a memory kernel superoperator, which can be obtained from short-lived (compared to the time scale of energy/charge transfer) projection-free inputs. In this paper, we test the ability of the M-GQME to predict the energy transfer dynamics within a seven-state benchmark model of the Fenna–Matthews–Olson (FMO) complex, with the short-lived projection-free inputs obtained via the Ehrenfest method. The M-GQME with Ehrenfest-based inputs is shown to yield accurate results across a wide parameter range. It is also found to dramatically outperform the direct application of the Ehrenfest method and to provide better-behaved convergence with respect to memory time in comparison to an alternative implementation of the GQME approach previously applied to the same FMO model.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.