### 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 language | English |
---|---|

Article number | 044112 |

Journal | Journal of Chemical Physics |

Volume | 138 |

Issue number | 4 |

DOIs | |

Publication status | Published - Jan 28 2013 |

### Fingerprint

### ASJC Scopus subject areas

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

### Cite this

*Journal of Chemical Physics*,

*138*(4), [044112]. https://doi.org/10.1063/1.4776230

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

Research output: Contribution to journal › Article

*Journal of Chemical Physics*, vol. 138, no. 4, 044112. https://doi.org/10.1063/1.4776230

}

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 -