Negative Zero-Point-Energy Parameter in the Meyer-Miller Mapping Model for Nonadiabatic Dynamics
05 Mar 2021-Journal of Physical Chemistry Letters (American Chemical Society (ACS))-Vol. 12, Iss: 10, pp 2496-2501
TL;DR: In this article, it was shown that a negative zero-point energy (ZPE) parameter could lead to a reasonably good performance in describing dynamic behaviors in typical spin-boson models for condensed-phase two-state systems, even at challenging zero temperature.
Abstract: The celebrated Meyer-Miller mapping model has been a useful approach for generating practical trajectory-based nonadiabatic dynamics methods. It is generally assumed that the zero-point-energy (ZPE) parameter is positive. The constraint implied in the conventional Meyer-Miller mapping Hamiltonian for an F-electronic-state system actually requires γ∈(-1/F, ∞) for the ZPE parameter for each electronic degree of freedom. Both negative and positive values are possible for such a parameter. We first establish a rigorous formulation to construct exact mapping models in the Cartesian phase space when the constraint is applied. When nuclear dynamics is approximated by the linearized semiclassical initial value representation, a negative ZPE parameter could lead to reasonably good performance in describing dynamic behaviors in typical spin-boson models for condensed-phase two-state systems, even at challenging zero temperature.
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: Proof-of-concept applications to realistic molecules demonstrate that machine learning phase space dynamics approaches are possible as well as competent in producing reasonably accurate results with a modest computation effort.
Abstract: Derived from phase space expressions of the quantum Liouville theorem, equilibrium continuity dynamics is a category of trajectory-based phase space dynamics methods, which satisfies the two critical fundamental criteria: conservation of the quantum Boltzmann distribution for the thermal equilibrium system and being exact for any thermal correlation functions (even of nonlinear operators) in the classical and harmonic limits. The effective force and effective mass matrix are important elements in the equations of motion of equilibrium continuity dynamics, where only the zeroth term of an exact series expansion of the phase space propagator is involved. We introduce a machine learning approach for fitting these elements in quantum phase space, leading to a much more efficient integration of the equations of motion. Proof-of-concept applications to realistic molecules demonstrate that machine learning phase space dynamics approaches are possible as well as competent in producing reasonably accurate results with a modest computation effort.
13 citations
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.
12 citations
TL;DR: In this article, a spin-mapping-based nonadiabatic ring polymer molecular dynamics (NRPMD) method was proposed, which is based on the spin mapping formalism.
Abstract: We present a new non-adiabatic ring polymer molecular dynamics (NRPMD) method based on the spin mapping formalism, which we refer to as the spin mapping NRPMD (SM-NRPMD) approach. We derive the path-integral partition function expression using the spin coherent state basis for the electronic states and the ring polymer formalism for the nuclear degrees of freedom. This partition function provides an efficient sampling of the quantum statistics. Using the basic properties of the Stratonovich–Weyl transformation, we further justify a Hamiltonian that we propose for the dynamical propagation of the coupled spin mapping variables and the nuclear ring polymer. The accuracy of the SM-NRPMD method is numerically demonstrated by computing the nuclear position and population auto-correlation functions of non-adiabatic model systems. The results obtained using the SM-NRPMD method agree very well with the numerically exact results. The main advantage of using the spin mapping variables over the harmonic oscillator mapping variables is numerically demonstrated, where the former provides nearly time-independent expectation values of physical observables for systems under thermal equilibrium. We also explicitly demonstrate that SM-NRPMD provides invariant dynamics upon various ways of partitioning the state-dependent and state-independent potentials.
7 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
References
More filters
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 paper, a test system is assumed to interact with a heat bath consisting of harmonic oscillators or an equivalent bath with a proper frequency spectrum producing a Gaussian-Markoffian random perturbation, and the effect of reaction of the test system to the bath is considered in the high temperature approximation.
Abstract: A test system is assumed to interact with a heat bath consisting of harmonic oscillators or an equivalent bath with a proper frequency spectrum producing a Gaussian-Markoffian random perturbation. The effect of reaction of the test system to the bath is considered in the high temperature approximation. Elimination of the bath using the influence functional method of Feynman and Vernon yields a continuous fraction expression for the reduced density matrix of the test system. The result affords a basis to clarify the relationship between the stochastic and the dynamical approaches to treat the problem of partial destruction of quantum coherence of a system interacting with its environment.
787 citations
TL;DR: In this paper, a multilayer formulation of the multiconfiguration time-dependent Hartree (MCTDH) theory is presented, where the single-particle (SP) functions in the original MCTDH method are further expressed employing a timedependent multi-figurational expansion, and the Dirac-Frenkel variational principle is applied to optimally determine the equations of motion.
Abstract: A multilayer (ML) formulation of the multiconfiguration time-dependent Hartree (MCTDH) theory is presented. In this new approach, the single-particle (SP) functions in the original MCTDH method are further expressed employing a time-dependent multiconfigurational expansion. The Dirac–Frenkel variational principle is then applied to optimally determine the equations of motion. Following this strategy, the SP groups are built in several layers, where each top layer SP can contain many more Cartesian degrees of freedom than in the previous formulation of the MCTDH method. As a result, the ML-MCTDH method has the capability of treating substantially more physical degrees of freedom than the original MCTDH method, and thus significantly enhances the ability of carrying out quantum dynamical simulations for complex molecular systems. The efficiency of the new formulation is demonstrated by converged quantum dynamical simulations for systems with a few hundred to a thousand degrees of freedom.
768 citations
TL;DR: In this paper, a formal exact classical analog for a finite dimensional (in Hilbert space) quantum mechanical system is defined for the electronic degrees of freedom in a molecular collision system, and the combination of this with the usual classical description of the heavy particle (i.e., nuclear) motion provides a completely classical model for electronic and heavy particle degree of freedom.
Abstract: It is shown how a formally exact classical analog can be defined for a finite dimensional (in Hilbert space) quantum mechanical system. This approach is then used to obtain a classical model for the electronic degrees of freedom in a molecular collision system, and the combination of this with the usual classical description of the heavy particle (i.e., nuclear) motion provides a completely classical model for the electronic and heavy particle degrees of freedom. The resulting equations of motion are shown to be equivalent to describing the electronic degrees of freedom by the time‐dependent Schrodinger equation, the time dependence arising from the classical motion of the nuclei, the trajectory of which is determined by the quantum mechanical average (i.e., Ehrenfest) force on the nuclei. Quantizing the system via classical S‐matrix theory is shown to provide a dynamically consistent description of nonadiabatic collision processes; i.e., different electronic transitions have different heavy particle trajectories and, for example, the total energy of the electronic and heavy particle degrees of freedom is conserved. Application of this classical model for the electronic degrees of freedom (plus classical S‐matrix theory) to the two‐state model problem shows that the approach provides a good description of the electronic dynamics.
688 citations
TL;DR: In this article, a semiclassical approach is presented that allows us to extend the usual Van Vleck-Gutzwiller formulation to the description of nonadiabatic quantum dynamics on coupled potential energy surfaces.
Abstract: A semiclassical approach is presented that allows us to extend the usual Van Vleck--Gutzwiller formulation to the description of nonadiabatic quantum dynamics on coupled potential-energy surfaces. Based on Schwinger's theory of angular momentum, the formulation employs an exact mapping of the discrete quantum variables onto continuous degrees of freedom. The resulting dynamical problem is evaluated through a semiclassical initial-value representation of the time-dependent propagator. As a first application we have performed semiclassical simulations for a spin-boson model, which reproduce the exact quantum-mechanical results quite accurately.
521 citations