Commutator Matrix in Phase Space Mapping Models for Nonadiabatic Quantum Dynamics.
02 Aug 2021-Journal of Physical Chemistry A (American Chemical Society (ACS))-Vol. 125, Iss: 31, pp 6845-6863
TL;DR: In this article, a general phase space mapping Hamiltonian for nonadiabatic systems, which is reminiscent of the renowned Meyer-Miller mapping Hamiltonians, involves a commutator variable matrix rather than the conventional zero-point-energy parameter.
Abstract: We show that a novel, general phase space mapping Hamiltonian for nonadiabatic systems, which is reminiscent of the renowned Meyer-Miller mapping Hamiltonian, involves a commutator variable matrix rather than the conventional zero-point-energy parameter. In the exact mapping formulation on constraint space for phase space approaches for nonadiabatic dynamics, the general mapping Hamiltonian with commutator variables can be employed to generate approximate trajectory-based dynamics. Various benchmark model tests, which range from gas phase to condensed phase systems, suggest that the overall performance of the general mapping Hamiltonian is better than that of the conventional Meyer-Miller Hamiltonian.
Citations
More filters
TL;DR: In this paper, a unified framework for constructing the mapping Hamiltonian on phase space for coupled F-state systems where the renowned Meyer-Miller Hamiltonian model is a special case is presented.
Abstract: ConspectusNonadiabatic dynamical processes are one of the most important quantum mechanical phenomena in chemical, materials, biological, and environmental molecular systems, where the coupling between different electronic states is either inherent in the molecular structure or induced by the (intense) external field. The curse of dimensionality indicates the intractable exponential scaling of calculation effort with system size and restricts the implementation of "numerically exact" approaches for realistic large systems. The phase space formulation of quantum mechanics offers an important theoretical framework for constructing practical approximate trajectory-based methods for quantum dynamics. This Account reviews our recent progress in phase space mapping theory: a unified framework for constructing the mapping Hamiltonian on phase space for coupled F-state systems where the renowned Meyer-Miller Hamiltonian model is a special case, a general phase space formulation of quantum mechanics for nonadiabatic systems where the electronic degrees of freedom are mapped onto constraint space and the nuclear degrees of freedom are mapped onto infinite space, and an isomorphism between the mapping phase space approach for nonadiabatic systems and that for nonequilibrium electron transport processes. While the zero-point-energy parameter is conventionally assumed to be positive, we show that the constraint implied in the conventional Meyer-Miller mapping Hamiltonian requires that such a parameter can be negative as well and lies in (-1/F, +∞) for each electronic degree of freedom. More importantly, the zero-point-energy parameter should be interpreted as a special case of a commutator matrix in the comprehensive phase space mapping Hamiltonian for nonadiabatic systems. From the rigorous formulation of mapping phase space, we propose approximate but practical trajectory-based nonadiabatic dynamics methods. The applications to both gas phase and condensed phase problems include the spin-boson model for condensed phase dissipative two-state systems, the three-state photodissociation models, the seven-site model of the Fenna-Matthews-Olson monomer in photosynthesis of green sulfur bacteria, the strongly coupled molecular/atomic matter-optical cavity systems designed for controlling and manipulating chemical dynamical processes, and the Landauer model for a quantum dot state coupled with two electrodes. In these applications the overall performance of our phase space mapping dynamics approach is superior to two prevailing trajectory-based methods, Ehrenfest dynamics and fewest switches surface hopping.
16 citations
TL;DR: In this paper , the authors report progress on the phase space formulation of quantum mechanics with coordinatemomentum variables, focusing more on new theory of (weighted) constraint coordinate momentum phase space for discrete variable quantum systems.
Abstract: We report recent progress on the phase space formulation of quantum mechanics with coordinate‐momentum variables, focusing more on new theory of (weighted) constraint coordinate‐momentum phase space for discrete‐variable quantum systems. This leads to a general coordinate‐momentum phase space formulation of composite quantum systems, where conventional representations on infinite phase space are employed for continuous variables. It is convenient to utilize (weighted) constraint coordinate‐momentum phase space for representing the quantum state and describing nonclassical features. Various numerical tests demonstrate that new trajectory‐based quantum dynamics approaches derived from the (weighted) constraint phase space representation are useful and practical for describing dynamical processes of composite quantum systems in the gas phase as well as in the condensed phase.
7 citations
TL;DR: In this article , the generalized spin mapping representation for non-adiabatic dynamics is presented, where the Stratonovich-Weyl transform is used to map an operator in the Hilbert space to a continuous function on the SU(N) Lie group.
Abstract: We present the rigorous theoretical framework of the generalized spin mapping representation for non-adiabatic dynamics. Our work is based upon a new mapping formalism recently introduced by Runeson and Richardson [J. Chem. Phys. 152, 084110 (2020)], which uses the generators of the su(N) Lie algebra to represent N discrete electronic states, thus preserving the size of the original Hilbert space. Following this interesting idea, the Stratonovich-Weyl transform is used to map an operator in the Hilbert space to a continuous function on the SU(N) Lie group, i.e., a smooth manifold which is a phase space of continuous variables. We further use the Wigner representation to describe the nuclear degrees of freedom and derive an exact expression of the time-correlation function as well as the exact quantum Liouvillian for the non-adiabatic system. Making the linearization approximation, this exact Liouvillian is reduced to the Liouvillian of several recently proposed methods, and the performance of this linearized method is tested using non-adiabatic models. We envision that the theoretical work presented here provides a rigorous and unified framework to formally derive non-adiabatic quantum dynamics approaches with continuous variables and connects the previous methods in a clear and concise manner.
6 citations
TL;DR: A one-shot trajectory learning approach that allows us to directly make an ultrafast prediction of the entire trajectory of the reduced density matrix for a new set of such simulation parameters as temperature and reorganization energy.
Abstract: Nonadiabatic quantum dynamics is important for understanding light-harvesting processes, but its propagation with traditional methods can be rather expensive. Here we present a one-shot trajectory learning approach that allows us to directly make an ultrafast prediction of the entire trajectory of the reduced density matrix for a new set of such simulation parameters as temperature and reorganization energy. The whole 10-ps-long propagation takes 70 ms as we demonstrate on the comparatively large quantum system, the Fenna-Matthews-Olsen (FMO) complex. Our approach also significantly reduces time and memory requirements for training.
6 citations
TL;DR: The NQCDynamics.jl package as mentioned in this paper provides a framework for performing semiclassical and mixed quantum-classical dynamics in the condensed phase using the Julia programming language.
Abstract: Accurate and efficient methods to simulate nonadiabatic and quantum nuclear effects in high-dimensional and dissipative systems are crucial for the prediction of chemical dynamics in the condensed phase. To facilitate effective development, code sharing, and uptake of newly developed dynamics methods, it is important that software implementations can be easily accessed and built upon. Using the Julia programming language, we have developed the NQCDynamics.jl package, which provides a framework for established and emerging methods for performing semiclassical and mixed quantum-classical dynamics in the condensed phase. The code provides several interfaces to existing atomistic simulation frameworks, electronic structure codes, and machine learning representations. In addition to the existing methods, the package provides infrastructure for developing and deploying new dynamics methods, which we hope will benefit reproducibility and code sharing in the field of condensed phase quantum dynamics. Herein, we present our code design choices and the specific Julia programming features from which they benefit. We further demonstrate the capabilities of the package on two examples of chemical dynamics in the condensed phase: the population dynamics of the spin-boson model as described by a wide variety of semiclassical and mixed quantum-classical nonadiabatic methods and the reactive scattering of H2 on Ag(111) using the molecular dynamics with electronic friction method. Together, they exemplify the broad scope of the package to study effective model Hamiltonians and realistic atomistic systems.
5 citations
References
More filters
TL;DR: In this article, the Boltzmann formula for the probability of a configuration is given in classical theory by means of a probability function, and the result discussed is developed for the correction term.
Abstract: The probability of a configuration is given in classical theory by the Boltzmann formula $\mathrm{exp}[\ensuremath{-}\frac{V}{\mathrm{hT}}]$ where $V$ is the potential energy of this configuration. For high temperatures this of course also holds in quantum theory. For lower temperatures, however, a correction term has to be introduced, which can be developed into a power series of $h$. The formula is developed for this correction by means of a probability function and the result discussed.
6,791 citations
TL;DR: In this article, the photon statistics of arbitrary fields in fully quantum-mechanical terms are discussed, and a general method of representing the density operator for the field is discussed as well as a simple formulation of a superposition law for photon fields.
Abstract: Methods are developed for discussing the photon statistics of arbitrary fields in fully quantum-mechanical terms. In order to keep the classical limit of quantum electrodynamics plainly in view, extensive use is made of the coherent states of the field. These states, which reduce the field correlation functions to factorized forms, are shown to offer a convenient basis for the description of fields of all types. Although they are not orthogonal to one another, the coherent states form a complete set. It is shown that any quantum state of the field may be expanded in terms of them in a unique way. Expansions are also developed for arbitrary operators in terms of products of the coherent state vectors. These expansions are discussed as a general method of representing the density operator for the field. A particular form is exhibited for the density operator which makes it possible to carry out many quantum-mechanical calculations by methods resembling those of classical theory. This representation permits clear insights into the essential distinction between the quantum and classical descriptions of the field. It leads, in addition, to a simple formulation of a superposition law for photon fields. Detailed discussions are given of the incoherent fields which are generated by superposing the outputs of many stationary sources. These fields are all shown to have intimately related properties, some of which have been known for the particular case of blackbody radiation.
5,372 citations
TL;DR: In this article, a functional-integral approach to the dynamics of a two-state system coupled to a dissipative environment is presented, and an exact and general prescription for the reduction, under appropriate circumstances, of the problem of a system tunneling between two wells in the presence of dissipative environments to the spin-boson problem is given.
Abstract: This paper presents the results of a functional-integral approach to the dynamics of a two-state system coupled to a dissipative environment. It is primarily an extended account of results obtained over the last four years by the authors; while they try to provide some background for orientation, it is emphatically not intended as a comprehensive review of the literature on the subject. Its contents include (1) an exact and general prescription for the reduction, under appropriate circumstances, of the problem of a system tunneling between two wells in the presence of a dissipative environment to the "spin-boson" problem; (2) the derivation of an exact formula for the dynamics of the latter problem; (3) the demonstration that there exists a simple approximation to this exact formula which is controlled, in the sense that we can put explicit bounds on the errors incurred in it, and that for almost all regions of the parameter space these errors are either very small in the limit of interest to us (the "slow-tunneling" limit) or can themselves be evaluated with satisfactory accuracy; (4) use of these results to obtain quantitative expressions for the dynamics of the system as a function of the spectral density $J(\ensuremath{\omega})$ of its coupling to the environment. If $J(\ensuremath{\omega})$ behaves as ${\ensuremath{\omega}}^{s}$ for frequencies of the order of the tunneling frequency or smaller, the authors find for the "unbiased" case the following results: For $sl1$ the system is localized at zero temperature, and at finite $T$ relaxes incoherently at a rate proportional to $\mathrm{exp}\ensuremath{-}{(\frac{{T}_{0}}{T})}^{1\ensuremath{-}s}$. For $sg2$ it undergoes underdamped coherent oscillations for all relevant temperatures, while for $1lsl2$ there is a crossover from coherent oscillation to overdamped relaxation as $T$ increases. Exact expressions for the oscillation and/or relaxation rates are presented in all these cases. For the "ohmic" case, $s=1$, the qualitative nature of the behavior depends critically on the dimensionless coupling strength $\ensuremath{\alpha}$ as well as the temperature $T$: over most of the ($\ensuremath{\alpha}$,$T$) plane (including the whole region $\ensuremath{\alpha}g1$) the behavior is an incoherent relaxation at a rate proportional to ${T}^{2\ensuremath{\alpha}\ensuremath{-}1}$, but for low $T$ and $0l\ensuremath{\alpha}l\frac{1}{2}$ the authors predict a combination of damped coherent oscillation and incoherent background which appears to disagree with the results of all previous approximations. The case of finite bias is also discussed.
4,047 citations
TL;DR: In this article, a method for carrying out molecular dynamics simulations of processes that involve electronic transitions is proposed, where the time dependent electronic Schrodinger equation is solved self-consistently with the classical mechanical equations of motion of the atoms.
Abstract: A method is proposed for carrying out molecular dynamics simulations of processes that involve electronic transitions. The time dependent electronic Schrodinger equation is solved self‐consistently with the classical mechanical equations of motion of the atoms. At each integration time step a decision is made whether to switch electronic states, according to probabilistic ‘‘fewest switches’’ algorithm. If a switch occurs, the component of velocity in the direction of the nonadiabatic coupling vector is adjusted to conserve energy. The procedure allows electronic transitions to occur anywhere among any number of coupled states, governed by the quantum mechanical probabilities. The method is tested against accurate quantal calculations for three one‐dimensional, two‐state models, two of which have been specifically designed to challenge any such mixed classical–quantal dynamical theory. Although there are some discrepancies, initial indications are encouraging. The model should be applicable to a wide variety of gas‐phase and condensed‐phase phenomena occurring even down to thermal energies.
3,173 citations
2,236 citations