Polarizable Potentials For Metals: The Density Readjusting Embedded Atom Method (DR-EAM)
arXiv:1904.00263 · doi:10.1103/PhysRevB.99.094106
Abstract
In simulations of metallic interfaces, a critical aspect of metallic behavior is missing from the some of the most widely used classical molecular dynamics force fields. We present a modification of the embedded atom method (EAM) which allows for electronic polarization of the metal by treating the valence density around each atom as a fluctuating dynamical quantity. The densities are represented by a set of additional fluctuating variables (and their conjugate momenta) which are propagated along with the nuclear coordinates. This ``density readjusting EAM'' (DR-EAM) preserves nearly all of the useful qualities of traditional EAM, including bulk elastic properties and surface energies. However, it also allows valence electron density to migrate through the metal in response to external perturbations. We show that DR-EAM can successfully model polarization in response to external charges, capturing the image charge effect in atomistic simulations. DR-EAM also captures some of the behavior of metals in the presence of uniform electric fields, predicting surface charging and shielding internal to the metal. We further show that it predicts charge transfer between the constituent atoms in alloys, leading to novel predictions about unit cell geometries in layered L structures.
33 pages, 6 figures
References in corpus (3)
Cited by in corpus (4)
- Interatomic potentials: Achievements and challenges
- Heat Transfer in Gold Interfaces Capped with Thiolated Polyethylene Glycol: A Molecular Dynamics Study
- Embedding quantum statistical excitations in a classical force field
- A review of heat transport in solvated gold nanoparticles: Molecular dynamics modeling and experimental perspectives