Polaron formation

Ehrenfest dynamics vs. exact results

Guangqi Li, Bijan Movaghar, Abraham Nitzan, Mark A Ratner

Research output: Contribution to journalArticle

16 Citations (Scopus)

Abstract

We use a one-dimensional tight binding model with an impurity site characterized by electron-vibration coupling, to describe electron transfer and localization at zero temperature, aiming to examine the process of polaron formation in this system. In particular we focus on comparing a semiclassical approach that describes nuclear motion in this many vibronic-states system on the Ehrenfest dynamics level to a numerically exact fully quantum calculation based on the Bonca-Trugman method [J. Bonča and S. A. Trugman, Phys. Rev. Lett. 75, 2566 (1995)]10.1103/PhysRevLett.75.2566. In both approaches, thermal relaxation in the nuclear subspace is implemented in equivalent approximate ways: In the Ehrenfest calculation the uncoupled (to the electronic subsystem) motion of the classical (harmonic) oscillator is simply damped as would be implied by coupling to a Markovian zero temperature bath. In the quantum calculation, thermal relaxation is implemented by augmenting the Liouville equation for the oscillator density matrix with kinetic terms that account for the same relaxation. In both cases we calculate the probability to trap the electron by forming a polaron and the probability that it escapes to infinity. Comparing these calculations, we find that while both result in similar long time yields for these processes, the Ehrenfest-dynamics based calculation fails to account for the correct time scale for the polaron formation. This failure results, as usual, from the fact that at the early stage of polaron formation the classical nuclear dynamics takes place on an unphysical average potential surface that reflects the distributed electronic population in the system, while the quantum calculation accounts fully for correlations between the electronic and vibrational subsystems.

Original languageEnglish
Article number044112
JournalJournal of Chemical Physics
Volume138
Issue number4
DOIs
Publication statusPublished - Jan 28 2013

Fingerprint

Gene Conversion
Electrons
Liouville equation
electronics
Liouville equations
infinity
harmonic oscillators
escape
baths
electron transfer
electrons
oscillators
traps
Impurities
impurities
vibration
Temperature
Kinetics
temperature
kinetics

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

Polaron formation : Ehrenfest dynamics vs. exact results. / Li, Guangqi; Movaghar, Bijan; Nitzan, Abraham; Ratner, Mark A.

In: Journal of Chemical Physics, Vol. 138, No. 4, 044112, 28.01.2013.

Research output: Contribution to journalArticle

Li, Guangqi ; Movaghar, Bijan ; Nitzan, Abraham ; Ratner, Mark A. / Polaron formation : Ehrenfest dynamics vs. exact results. In: Journal of Chemical Physics. 2013 ; Vol. 138, No. 4.
@article{df93a6f99feb4379a93bd75d041b9d77,
title = "Polaron formation: Ehrenfest dynamics vs. exact results",
abstract = "We use a one-dimensional tight binding model with an impurity site characterized by electron-vibration coupling, to describe electron transfer and localization at zero temperature, aiming to examine the process of polaron formation in this system. In particular we focus on comparing a semiclassical approach that describes nuclear motion in this many vibronic-states system on the Ehrenfest dynamics level to a numerically exact fully quantum calculation based on the Bonca-Trugman method [J. Bonča and S. A. Trugman, Phys. Rev. Lett. 75, 2566 (1995)]10.1103/PhysRevLett.75.2566. In both approaches, thermal relaxation in the nuclear subspace is implemented in equivalent approximate ways: In the Ehrenfest calculation the uncoupled (to the electronic subsystem) motion of the classical (harmonic) oscillator is simply damped as would be implied by coupling to a Markovian zero temperature bath. In the quantum calculation, thermal relaxation is implemented by augmenting the Liouville equation for the oscillator density matrix with kinetic terms that account for the same relaxation. In both cases we calculate the probability to trap the electron by forming a polaron and the probability that it escapes to infinity. Comparing these calculations, we find that while both result in similar long time yields for these processes, the Ehrenfest-dynamics based calculation fails to account for the correct time scale for the polaron formation. This failure results, as usual, from the fact that at the early stage of polaron formation the classical nuclear dynamics takes place on an unphysical average potential surface that reflects the distributed electronic population in the system, while the quantum calculation accounts fully for correlations between the electronic and vibrational subsystems.",
author = "Guangqi Li and Bijan Movaghar and Abraham Nitzan and Ratner, {Mark A}",
year = "2013",
month = "1",
day = "28",
doi = "10.1063/1.4776230",
language = "English",
volume = "138",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "4",

}

TY - JOUR

T1 - Polaron formation

T2 - Ehrenfest dynamics vs. exact results

AU - Li, Guangqi

AU - Movaghar, Bijan

AU - Nitzan, Abraham

AU - Ratner, Mark A

PY - 2013/1/28

Y1 - 2013/1/28

N2 - We use a one-dimensional tight binding model with an impurity site characterized by electron-vibration coupling, to describe electron transfer and localization at zero temperature, aiming to examine the process of polaron formation in this system. In particular we focus on comparing a semiclassical approach that describes nuclear motion in this many vibronic-states system on the Ehrenfest dynamics level to a numerically exact fully quantum calculation based on the Bonca-Trugman method [J. Bonča and S. A. Trugman, Phys. Rev. Lett. 75, 2566 (1995)]10.1103/PhysRevLett.75.2566. In both approaches, thermal relaxation in the nuclear subspace is implemented in equivalent approximate ways: In the Ehrenfest calculation the uncoupled (to the electronic subsystem) motion of the classical (harmonic) oscillator is simply damped as would be implied by coupling to a Markovian zero temperature bath. In the quantum calculation, thermal relaxation is implemented by augmenting the Liouville equation for the oscillator density matrix with kinetic terms that account for the same relaxation. In both cases we calculate the probability to trap the electron by forming a polaron and the probability that it escapes to infinity. Comparing these calculations, we find that while both result in similar long time yields for these processes, the Ehrenfest-dynamics based calculation fails to account for the correct time scale for the polaron formation. This failure results, as usual, from the fact that at the early stage of polaron formation the classical nuclear dynamics takes place on an unphysical average potential surface that reflects the distributed electronic population in the system, while the quantum calculation accounts fully for correlations between the electronic and vibrational subsystems.

AB - We use a one-dimensional tight binding model with an impurity site characterized by electron-vibration coupling, to describe electron transfer and localization at zero temperature, aiming to examine the process of polaron formation in this system. In particular we focus on comparing a semiclassical approach that describes nuclear motion in this many vibronic-states system on the Ehrenfest dynamics level to a numerically exact fully quantum calculation based on the Bonca-Trugman method [J. Bonča and S. A. Trugman, Phys. Rev. Lett. 75, 2566 (1995)]10.1103/PhysRevLett.75.2566. In both approaches, thermal relaxation in the nuclear subspace is implemented in equivalent approximate ways: In the Ehrenfest calculation the uncoupled (to the electronic subsystem) motion of the classical (harmonic) oscillator is simply damped as would be implied by coupling to a Markovian zero temperature bath. In the quantum calculation, thermal relaxation is implemented by augmenting the Liouville equation for the oscillator density matrix with kinetic terms that account for the same relaxation. In both cases we calculate the probability to trap the electron by forming a polaron and the probability that it escapes to infinity. Comparing these calculations, we find that while both result in similar long time yields for these processes, the Ehrenfest-dynamics based calculation fails to account for the correct time scale for the polaron formation. This failure results, as usual, from the fact that at the early stage of polaron formation the classical nuclear dynamics takes place on an unphysical average potential surface that reflects the distributed electronic population in the system, while the quantum calculation accounts fully for correlations between the electronic and vibrational subsystems.

UR - http://www.scopus.com/inward/record.url?scp=84873589568&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84873589568&partnerID=8YFLogxK

U2 - 10.1063/1.4776230

DO - 10.1063/1.4776230

M3 - Article

VL - 138

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 4

M1 - 044112

ER -