Classically exact surface diffusion constants at arbitrary temperature

Abstract

An expression is presented for computing the classical diffusion constant of a point defect (e.g., an adatom) in an infinite lattice of binding sites at arbitrary temperature. The transition state theory diffusion constant is simply multiplied by a dynamical correction factor that is computed from short-time classical trajectories initiated at the site boundaries. The time scale limitations of direct molecular dynamics are thus avoided in the low- and middle-temperature regimes. The expression results from taking the time derivative of the particle mean-square displacement in the lattice-discretized coordinate system. Applications are presented for surface diffusion on fcc(100) and fcc(111) Lennard-Jones crystal faces.