skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Electron-phonon interaction within classical molecular dynamics

; ; ; ; ;
Publication Date:
Research Org.:
Energy Frontier Research Centers (EFRC) (United States). Energy Dissipation to Defect Evolution (EDDE)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22)
OSTI Identifier:
DOE Contract Number:
Resource Type:
Journal Article
Resource Relation:
Journal Name: Physical Review B; Journal Volume: 94; Journal Issue: 2; Related Information: EDDE partners with Oak Ridge National Laboratory (lead); Lawrence Livermore National Laboratory; University of Michigan; University of Tennessee; University of Wisconsin; University of Wyoming; Virginia Tech
Country of Publication:
United States
phonons, thermal conductivity, nuclear (including radiation effects), defects, spin dynamics, materials and chemistry by design, synthesis (novel materials)

Citation Formats

Tamm, A., Samolyuk, G., Correa, A. A., Klintenberg, M., Aabloo, A., and Caro, A.. Electron-phonon interaction within classical molecular dynamics. United States: N. p., 2016. Web. doi:10.1103/PhysRevB.94.024305.
Tamm, A., Samolyuk, G., Correa, A. A., Klintenberg, M., Aabloo, A., & Caro, A.. Electron-phonon interaction within classical molecular dynamics. United States. doi:10.1103/PhysRevB.94.024305.
Tamm, A., Samolyuk, G., Correa, A. A., Klintenberg, M., Aabloo, A., and Caro, A.. 2016. "Electron-phonon interaction within classical molecular dynamics". United States. doi:10.1103/PhysRevB.94.024305.
title = {Electron-phonon interaction within classical molecular dynamics},
author = {Tamm, A. and Samolyuk, G. and Correa, A. A. and Klintenberg, M. and Aabloo, A. and Caro, A.},
abstractNote = {},
doi = {10.1103/PhysRevB.94.024305},
journal = {Physical Review B},
number = 2,
volume = 94,
place = {United States},
year = 2016,
month = 7
  • Here, we present a model for nonadiabatic classical molecular dynamics simulations that captures with high accuracy the wave-vector q dependence of the phonon lifetimes, in agreement with quantum mechanics calculations. It is based on a local view of the e-ph interaction where individual atom dynamics couples to electrons via a damping term that is obtained as the low-velocity limit of the stopping power of a moving ion in a host. The model is parameter free, as its components are derived from ab initio-type calculations, is readily extended to the case of alloys, and is adequate for large-scale molecular dynamics computermore » simulations. We also show how this model removes some oversimplifications of the traditional ionic damped dynamics commonly used to describe situations beyond the Born-Oppenheimer approximation.« less
  • Adiabatic mixed quantum/classical (MQC) molecular dynamics (MD) simulations were used to generate snapshots of the hydrated electron in liquid water at 300 K. Water cluster anions that include two complete solvation shells centered on the hydrated electron were extracted from the MQC MD simulations and embedded in a roughly 18 Angstrom x 18 Angstrom x 18 Angstrom matrix of fractional point charges designed to represent the rest of the solvent. Density functional theory (DFT) with the Becke-Lee-Yang-Parr functional and single-excitation configuration interaction (CIS) methods were then applied to these embedded clusters. The salient feature of these hybrid DFT(CIS)/MQC MD calculationsmore » is significant transfer ({approx}18%) of the excess electron's charge density into the 2p orbitals of oxygen atoms in OH groups forming the solvation cavity. We used the results of these calculations to examine the structure of the singly occupied and the lower unoccupied molecular orbitals, the density of states, the absorption spectra in the visible and ultraviolet, the hyperfine coupling (hfcc) tensors, and the infrared (IR) and Raman spectra of these embedded water cluster anions. The calculated hfcc tensors were used to compute electron paramagnetic resonance (EPR) and electron spin echo envelope modulation (ESEEM) spectra for the hydrated electron that compared favorably to the experimental spectra of trapped electrons in alkaline ice. The calculated vibrational spectra of the hydrated electron are consistent with the red-shifted bending and stretching frequencies observed in resonance Raman experiments. In addition to reproducing the visible/near IR absorption spectrum, the hybrid DFT model also accounts for the hydrated electron's 190-nm absorption band in the ultraviolet. Thus, our study suggests that to explain several important experimentally observed properties of the hydrated electron, many-electron effects must be accounted for: one-electron models that do not allow for mixing of the excess electron density with the frontier orbitals of the first-shell solvent molecules cannot explain the observed magnetic, vibrational, and electronic properties of this species. Despite the need for multielectron effects to explain these important properties, the ensemble-averaged radial wavefunctions and energetics of the highest occupied and three lowest unoccupied orbitals of the hydrated electrons in our hybrid model are close to the s- and p-like states obtained in one-electron models. Thus, one-electron models can provide a remarkably good approximation to the multielectron picture of the hydrated electron for many applications; indeed, the two approaches appear to be complementary.« less
  • Phonon dispersion and phonon density of states of graphene are calculated using classical molecular dynamics simulations. In this method, the dynamical matrix is constructed based on linear response theory by computing the displacement of atoms during the simulations. The computed phonon dispersions show excellent agreement with experiments. The simulations are done in both NVT and NPT ensembles at 300 K and found that the LO/TO modes are getting hardened at the Γ point. The NPT ensemble simulations capture the anharmonicity of the crystal accurately and the hardening of LO/TO modes is more pronounced. We also found that at 300 Kmore » the C-C bond length reduces below the equilibrium value and the ZA bending mode frequency becomes imaginary close to Γ along K-Γ direction, which indicates instability of the flat 2D graphene sheets.« less
  • Electron-phonon interaction is considered within the framework of the fluctuating valence of Cu atoms. Anderson's lattice Hamiltonian is suitably modified to take this into account. Using Green's function technique the possible quasiparticle excitations are determined. The quantity 2{Delta}{sub k}(0)/k{sub B}T{sub c} is calculated for T{sub c} = 40 K. The calculated values are in good agreement with the experimental results.