Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

Simulation of environmental effects on coherent quantum dynamics in many-body systems

Jeanne M. Riga, Craig C. Martens

Published Mar 29, 2004• 28 citations
Publisher

Abstract

In this paper we describe an application of the trajectory-based semiclassical Liouville method for modeling coherent molecular dynamics on multiple electronic surfaces to the treatment of the evolution and decay of quantum electronic coherence in many-body systems. We consider a model representing the coherent evolution of quantum wave packets on two excited electronic surfaces of a diatomic molecule in the gas phase and in rare gas solvent environments, ranging from small clusters to a cryogenic solid. For the gas phase system, the semiclassical trajectory method is shown to reproduce the evolution of the electronic-nuclear coherence nearly quantitatively. The dynamics of decoherence are then investigated for the solvated systems using the semiclassical approach. It is found that, although solvation in general leads to more rapid and extensive loss of quantum coherence, the details of the coupled system-bath dynamics are important, and in some cases the environment can preserve or even enhance quantum coherence beyond that seen in the isolated system.

DOI
10.1063/1.1651472

Related papers

Simulation of vibrational dephasing of I2 in solid Kr using the semiclassical Liouville method

Jeanne M. Riga, Erick Fredj, Craig C. Martens

In this paper, we present simulations of the decay of quantum coherence between vibrational states of I(2) in its ground (X) electronic state embedded in a cryogenic Kr matrix. We employ a numerical method based on the semiclassical limit of the quantum Liouville equation, which allows the simulation of the evolution and decay of quantum vibrational coherence using classical trajectories and ensemble averaging. The vibrational level-dependent interaction of the I(2)(X) oscillator with the rare-gas environment is modeled using a recently developed method for constructing state-dependent many-body potentials for quantum vibrations in a many-body classical environment [J. M. Riga, E. Fredj, and C. C. Martens, J. Chem. Phys. 122, 174107 (2005)]. The vibrational dephasing rates gamma(0n) for coherences prepared between the ground vibrational state mid R:0 and excited vibrational state mid R:n are calculated as a function of n and lattice temperature T. Excellent agreement with recent experiments performed by Karavitis et al. [Phys. Chem. Chem. Phys. 7, 791 (2005)] is obtained.

The Journal of Chemical Physics•Feb 9, 2006•23 citations

Theory and simulation of the loss of coherence in thermal and nonequilibrium environments

Craig C. Martens

In this paper, we describe recent work on modelling the loss of quantum coherence in condensed phase systems. We review numerical methodology based on a semiclassical limit of the quantum Liouville equation for simulating quantum coherent processes using classical-like molecular dynamics simulation and ensemble averaging, and apply the approach to simulating vibrational dephasing of I 2 in cryogenic rare gas matrices and the quantum vibrations of OH stretches of HOD in D 2 O. We then describe a simple analytic and numerical model that highlights novel behaviour that can be exhibited by quantum coherent processes in the presence of a nonequilibrium environment.

Journal of Physics B Atomic Molecular and Optical Physics•Jul 30, 2012•12 citations

Classical mechanics and the spreading of localized wave packets in condensed phase molecular systems

Michael F. Herman, D. F. Coker

The relationship between the diverging of classical trajectories in chaotic many-body systems, the spreading of quantum wave packets, and the validity and use of classical molecular dynamics is explored. This analysis, which is based on the semiclassical description of wave function propagation in terms of a weighted integration over a traveling fixed width coherent state basis, suggests that the exponential divergence of nearby classical trajectories in chaotic many-body systems should result in the rapid delocalization of an initially localized quantum wave packet describing the state of the system. Thus the justification for the use of classical molecular dynamics procedures for these supposedly classical systems cannot be based on the picture of the system wave function remaining localized as its center follows a nearly classical trajectory. The quantum evolution of the system density, on the other hand, requires two propagators, and each of these propagators is represented as an integration over trajectories in the semiclassical picture. The interference between the contributions from these two integrations over classical trajectories focuses the analysis on the most important points in this trajectory pair space, which are shown to occur when both trajectories in the pair are the same. Given reasonable assumptions for the initial density for a system that is expected to be well described by classical molecular dynamics, and given an appropriate choice for the width of the coherent state basis which is employed in the semiclassical description, it is shown that the semiclassical expressions for time dependent observables and correlation functions reduce the purely classical expressions, despite the fact that an initially localized wave packet would rapidly delocalize for the same system.

The Journal of Chemical Physics•Aug 1, 1999•43 citations

Non-Born-Oppenheimer Liouville-von Neumann Dynamics. Evolution of a Subsystem Controlled by Linear and Population-Driven Decay of Mixing with Decoherent and Coherent Switching

Chaoyuan Zhu, Ahren W. Jasper, Donald G. Truhlar

Electronic energy flow in an isolated molecular system involves coupling between the electronic and nuclear subsystems, and the coupled system evolves to a statistical mixture of pure states. In semiclassical theories, nuclear motion is treated using classical mechanics, and electronic motion is treated as an open quantal system coupled to a "bath" of nuclear coordinates. We have previously shown how this can be simulated by a time-dependent Schrödinger equation with coherent switching and decay of mixing, where the decay of mixing terms model the dissipative effect of the environment on the electronic subdynamics (i.e., on the reduced dynamics of the electronic subsystem). In the present paper we reformulate the problem as a Liouville-von Neumann equation of motion (i.e., we propagate the reduced density matrix of the electronic subsystem), and we introduce the assumption of first-order linear decay. We specifically examine the cases of equal relaxation times for both longitudinal (i.e., population) decay and transverse decay (i.e., dephasing) and of longitudinal relaxation only, yielding the linear decay of mixing (LDM) and the population-driven decay of mixing (PDDM) schemes, respectively. Because we do not generally know the basis in which coherence decays, that is, the pointer basis, we judge the semiclassical methods in part by their ability to give good results in both the adiabatic and diabatic bases. The accuracy in the prediction of physical observables is shown to be robust not only with respect to basis but also with respect to the way in which demixing is incorporated into the master equation for the density matrix. The success of the PDDM scheme is particularly interesting because it incorporates the least amount of decoherence (i.e., the PDDM scheme is the most similar of the methods discussed to the fully coherent semiclassical Ehrenfest method). For both the new and previous decay of mixing schemes, four kinds of decoherent state switching algorithms are analyzed and compared to one another: natural switching (NS), self-consistent switching (SCS), coherent switching (CS), and globally coherent switching (GCS). The CS formulations are examples of a non-Markovian method, in which the system retains some memory of its history, whereas the GCS, SCS, and NS schemes are Markovian (time local). These methods are tested against accurate quantum mechanical results using 17 multidimensional atom-diatom test cases. The test cases include avoided crossings, conical interactions, and systems with noncrossing diabatic potential energy surfaces. The CS switching algorithm, in which the state populations are controlled by a coherent stochastic algorithm for each complete passage through a strong interaction region, but successive strong-interaction regions are not mutually coherent, is shown to be the most accurate of the switching algorithms tested for the LDM and PDDM methods as well as for the previous decay of mixing methods, which are reformulated here as Liouville-von Neumann equations with nonlinear decay of mixing (NLDM). We also demonstrate that one variant of the PDDM method with CS performs almost equally well in the adiabatic and diabatic representations, which is a difficult objective for semiclassical methods. Thus decay of mixing methods provides powerful mixed quantum-classical methods for modeling non-Born-Oppenheimer polyatomic dynamics including photochemistry, charge-transfer, and other electronically nonadiabatic processes.

Journal of Chemical Theory and Computation•Jun 14, 2005•141 citations

Lessons on electronic decoherence in molecules from exact modeling

Wenxiang Hu, Bing Gu, Ignacio Franco

Electronic decoherence processes in molecules and materials are usually thought and modeled via schemes for the system-bath evolution in which the bath is treated either implicitly or approximately. Here we present computations of the electronic decoherence dynamics of a model many-body molecular system described by the Su-Schrieffer-Heeger Hamiltonian with Hubbard electron-electron interactions using an exact method in which both electronic and nuclear degrees of freedom are taken into account explicitly and fully quantum mechanically. To represent the electron-nuclear Hamiltonian in matrix form and propagate the dynamics, the computations employ the Jordan-Wigner transformation for the fermionic creation/annihilation operators and the discrete variable representation for the nuclear operators. The simulations offer a standard for electronic decoherence that can be used to test approximations. They also provide a useful platform to answer fundamental questions about electronic decoherence that cannot be addressed through approximate or implicit schemes. Specifically, through simulations, we isolate basic mechanisms for electronic coherence loss and demonstrate that electronic decoherence is possible even for one-dimensional nuclear bath. Furthermore, we show that (i) decreasing the mass of the bath generally leads to faster electronic decoherence; (ii) electron-electron interactions strongly affect the electronic decoherence when the electron-nuclear dynamics is not pure-dephasing; (iii) classical bath models with initial conditions sampled from the Wigner distribution accurately capture the short-time electronic decoherence dynamics; (iv) model separable initial superpositions often used to understand decoherence after photoexcitation are only relevant in experiments that employ delta-like laser pulses to initiate the dynamics. These insights can be employed to interpret and properly model coherence phenomena in molecules.

arXiv:1801.05846•The Journal of Chemical Physics•Apr 3, 2018•35 citations

Coherent quantum processes in thermal and nonequilibrium environments

Craig C. Martens

Abstract In this article, we describe recent work on investigating the role of the environment in influencing coherent quantum dynamics. We review numerical methodology based on a semiclassical limit of the quantum Liouville equation for simulating quantum coherent processes using classical‐like molecular dynamics simulation and ensemble averaging and apply the approach to simulate vibrational dephasing of I2 in cryogenic rare gas matrices and the quantum vibrations of OH stretches of HOD in D2O. We then describe a simple analytic and numerical model that highlights novel behavior that can be exhibited by quantum coherent processes in the presence of an environment that is not at thermal equilibrium. © 2012 Wiley Periodicals, Inc.

International Journal of Quantum Chemistry•Jun 13, 2012•1 citation
In this paper we describe an application of the trajectory-based semiclassical Liouville method for modeling coherent molecular dynamics on multiple electronic surfaces to the treatment of the evolution and decay of quantum electronic coherence in many-body systems. We consider a model representing the coherent evolution of quantum wave packets on two excited electronic surfaces of a diatomic molecule in the gas phase and in rare gas solvent environments, ranging from small clusters to a cryogenic solid. For the gas phase system, the semiclassical trajectory method is shown to reproduce the evolution of the electronic-nuclear coherence nearly quantitatively. The dynamics of decoherence are then investigated for the solvated systems using the semiclassical approach. It is found that, although solvation in general leads to more rapid and extensive loss of quantum coherence, the details of the coupled system-bath dynamics are important, and in some cases the environment can preserve or even enhance quantum coherence beyond that seen in the isolated system.
In this paper, we present simulations of the decay of quantum coherence between vibrational states of I(2) in its ground (X) electronic state embedded in a cryogenic Kr matrix. We employ a numerical method based on the semiclassical limit of the quantum Liouville equation, which allows the simulation of the evolution and decay of quantum vibrational coherence using classical trajectories and ensemble averaging. The vibrational level-dependent interaction of the I(2)(X) oscillator with the rare-gas environment is modeled using a recently developed method for constructing state-dependent many-body potentials for quantum vibrations in a many-body classical environment [J. M. Riga, E. Fredj, and C. C. Martens, J. Chem. Phys. 122, 174107 (2005)]. The vibrational dephasing rates gamma(0n) for coherences prepared between the ground vibrational state mid R:0 and excited vibrational state mid R:n are calculated as a function of n and lattice temperature T. Excellent agreement with recent experiments performed by Karavitis et al. [Phys. Chem. Chem. Phys. 7, 791 (2005)] is obtained.
In this paper, we describe recent work on modelling the loss of quantum coherence in condensed phase systems. We review numerical methodology based on a semiclassical limit of the quantum Liouville equation for simulating quantum coherent processes using classical-like molecular dynamics simulation and ensemble averaging, and apply the approach to simulating vibrational dephasing of I 2 in cryogenic rare gas matrices and the quantum vibrations of OH stretches of HOD in D 2 O. We then describe a simple analytic and numerical model that highlights novel behaviour that can be exhibited by quantum coherent processes in the presence of a nonequilibrium environment.
The relationship between the diverging of classical trajectories in chaotic many-body systems, the spreading of quantum wave packets, and the validity and use of classical molecular dynamics is explored. This analysis, which is based on the semiclassical description of wave function propagation in terms of a weighted integration over a traveling fixed width coherent state basis, suggests that the exponential divergence of nearby classical trajectories in chaotic many-body systems should result in the rapid delocalization of an initially localized quantum wave packet describing the state of the system. Thus the justification for the use of classical molecular dynamics procedures for these supposedly classical systems cannot be based on the picture of the system wave function remaining localized as its center follows a nearly classical trajectory. The quantum evolution of the system density, on the other hand, requires two propagators, and each of these propagators is represented as an integration over trajectories in the semiclassical picture. The interference between the contributions from these two integrations over classical trajectories focuses the analysis on the most important points in this trajectory pair space, which are shown to occur when both trajectories in the pair are the same. Given reasonable assumptions for the initial density for a system that is expected to be well described by classical molecular dynamics, and given an appropriate choice for the width of the coherent state basis which is employed in the semiclassical description, it is shown that the semiclassical expressions for time dependent observables and correlation functions reduce the purely classical expressions, despite the fact that an initially localized wave packet would rapidly delocalize for the same system.
Electronic energy flow in an isolated molecular system involves coupling between the electronic and nuclear subsystems, and the coupled system evolves to a statistical mixture of pure states. In semiclassical theories, nuclear motion is treated using classical mechanics, and electronic motion is treated as an open quantal system coupled to a "bath" of nuclear coordinates. We have previously shown how this can be simulated by a time-dependent Schrödinger equation with coherent switching and decay of mixing, where the decay of mixing terms model the dissipative effect of the environment on the electronic subdynamics (i.e., on the reduced dynamics of the electronic subsystem). In the present paper we reformulate the problem as a Liouville-von Neumann equation of motion (i.e., we propagate the reduced density matrix of the electronic subsystem), and we introduce the assumption of first-order linear decay. We specifically examine the cases of equal relaxation times for both longitudinal (i.e., population) decay and transverse decay (i.e., dephasing) and of longitudinal relaxation only, yielding the linear decay of mixing (LDM) and the population-driven decay of mixing (PDDM) schemes, respectively. Because we do not generally know the basis in which coherence decays, that is, the pointer basis, we judge the semiclassical methods in part by their ability to give good results in both the adiabatic and diabatic bases. The accuracy in the prediction of physical observables is shown to be robust not only with respect to basis but also with respect to the way in which demixing is incorporated into the master equation for the density matrix. The success of the PDDM scheme is particularly interesting because it incorporates the least amount of decoherence (i.e., the PDDM scheme is the most similar of the methods discussed to the fully coherent semiclassical Ehrenfest method). For both the new and previous decay of mixing schemes, four kinds of decoherent state switching algorithms are analyzed and compared to one another: natural switching (NS), self-consistent switching (SCS), coherent switching (CS), and globally coherent switching (GCS). The CS formulations are examples of a non-Markovian method, in which the system retains some memory of its history, whereas the GCS, SCS, and NS schemes are Markovian (time local). These methods are tested against accurate quantum mechanical results using 17 multidimensional atom-diatom test cases. The test cases include avoided crossings, conical interactions, and systems with noncrossing diabatic potential energy surfaces. The CS switching algorithm, in which the state populations are controlled by a coherent stochastic algorithm for each complete passage through a strong interaction region, but successive strong-interaction regions are not mutually coherent, is shown to be the most accurate of the switching algorithms tested for the LDM and PDDM methods as well as for the previous decay of mixing methods, which are reformulated here as Liouville-von Neumann equations with nonlinear decay of mixing (NLDM). We also demonstrate that one variant of the PDDM method with CS performs almost equally well in the adiabatic and diabatic representations, which is a difficult objective for semiclassical methods. Thus decay of mixing methods provides powerful mixed quantum-classical methods for modeling non-Born-Oppenheimer polyatomic dynamics including photochemistry, charge-transfer, and other electronically nonadiabatic processes.
Electronic decoherence processes in molecules and materials are usually thought and modeled via schemes for the system-bath evolution in which the bath is treated either implicitly or approximately. Here we present computations of the electronic decoherence dynamics of a model many-body molecular system described by the Su-Schrieffer-Heeger Hamiltonian with Hubbard electron-electron interactions using an exact method in which both electronic and nuclear degrees of freedom are taken into account explicitly and fully quantum mechanically. To represent the electron-nuclear Hamiltonian in matrix form and propagate the dynamics, the computations employ the Jordan-Wigner transformation for the fermionic creation/annihilation operators and the discrete variable representation for the nuclear operators. The simulations offer a standard for electronic decoherence that can be used to test approximations. They also provide a useful platform to answer fundamental questions about electronic decoherence that cannot be addressed through approximate or implicit schemes. Specifically, through simulations, we isolate basic mechanisms for electronic coherence loss and demonstrate that electronic decoherence is possible even for one-dimensional nuclear bath. Furthermore, we show that (i) decreasing the mass of the bath generally leads to faster electronic decoherence; (ii) electron-electron interactions strongly affect the electronic decoherence when the electron-nuclear dynamics is not pure-dephasing; (iii) classical bath models with initial conditions sampled from the Wigner distribution accurately capture the short-time electronic decoherence dynamics; (iv) model separable initial superpositions often used to understand decoherence after photoexcitation are only relevant in experiments that employ delta-like laser pulses to initiate the dynamics. These insights can be employed to interpret and properly model coherence phenomena in molecules.
Abstract In this article, we describe recent work on investigating the role of the environment in influencing coherent quantum dynamics. We review numerical methodology based on a semiclassical limit of the quantum Liouville equation for simulating quantum coherent processes using classical‐like molecular dynamics simulation and ensemble averaging and apply the approach to simulate vibrational dephasing of I2 in cryogenic rare gas matrices and the quantum vibrations of OH stretches of HOD in D2O. We then describe a simple analytic and numerical model that highlights novel behavior that can be exhibited by quantum coherent processes in the presence of an environment that is not at thermal equilibrium. © 2012 Wiley Periodicals, Inc.