Journal Article

# 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
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.

Topics: Phase space (52%), ,
##### Citations
More

5 results found

Journal Article
Xinzijian Liu1, Linfeng Zhang, Jian Liu1Institutions (1)
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.

Topics: Phase space (64%), Quantum dynamics (62%), Equations of motion (55%) ... read more

7 Citations

Open accessJournal Article
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.

Topics: , Coherent states (54%), Spin-½ (54%) ... read more

2 Citations

Open accessJournal Article
Xin He1, Baihua Wu1, Zhihao Gong1, Jian Liu1Institutions (1)
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.

Topics: , Phase space (59%), Commutator (55%) ... read more

1 Citations

Open accessJournal Article
Jian Liu1, Xin He1, Baihua Wu1Institutions (1)
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.

Topics: Phase space (65%), Surface hopping (60%),  ... read more

Journal Article
Zhubin Hu1, Dominikus Brian1, Xiang Sun1Institutions (1)
Abstract: Model Hamiltonians constructed from quantum chemistry calculations and molecular dynamics simulations are widely used for simulating nonadiabatic dynamics in the condensed phase. The most popular two-state spin-boson model could be built by mapping the all-atom anharmonic Hamiltonian onto a two-level system bilinearly coupled to a harmonic bath using the energy gap time correlation function. However, for more than two states, there lacks a general strategy to construct multi-state harmonic (MSH) models since the energy gaps between different pairs of electronic states are not entirely independent and need to be considered consistently. In this paper, we extend the previously proposed approach for building three-state harmonic models for photoinduced charge transfer to the arbitrary number of electronic states with a globally shared bath and the system-bath couplings are scaled differently according to the reorganization energies between each pair of states. We demonstrate the MSH model construction for an organic photovoltaic carotenoid-porphyrin-C60 molecular triad dissolved in explicit tetrahydrofuran solvent. Nonadiabatic dynamics was simulated using mixed quantum-classical techniques, including the linearized semiclassical and symmetrical quasiclassical dynamics with the mapping Hamiltonians, mean-field Ehrenfest, and mixed quantum-classical Liouville dynamics in two-state, three-state, and four-state harmonic models of the triad system. The MSH models are shown to provide a general and flexible framework for simulating nonadiabatic dynamics in complex systems.

Topics: Anharmonicity (53%),
##### References
More

32 results found

Open accessJournal Article
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.

Topics: , Order (ring theory) (62%),  ... read more

3,752 Citations

Journal Article
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.

Topics: , Degrees of freedom (51%), Hartree (50%)

668 Citations

Journal Article
Yoshitaka Tanimura1, Ryogo Kubo1Institutions (1)
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.

Topics: Quantum system (55%), Redfield equation (52%), Time evolution (50%)

660 Citations

Journal ArticleDOI: 10.1063/1.437910
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.

Topics: , Equations of motion (53%), Schrödinger equation (51%) ... read more

623 Citations

Journal Article
Gerhard Stock1, Michael Thoss1Institutions (1)
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.