Can Imaginary Instantaneous Normal Mode Frequencies Predict Barriers to Self-Diffusion?
J. Daniel Gezelter, Eran Rabani, and B. J. Berne
Department of Chemistry and Center for Biomolecular Simulation,
Columbia University, 3000 Broadway, New York, New York 10027
We discuss whether or not local information on the potential energy
surface embodied by the distribution of unstable instantaneous normal
modes can be used to predict the hopping rates and barrier heights for
Zwanzig's model of self-diffusion [R. Zwanzig, J. Chem. Phys.
79, 4507 (1983)]. Results from a set of simulations done at
multiple temperatures and densities are presented. These simulations
show that the theories which predict barrier heights from the
distribution of imaginary frequencies are questionable. This
discrepancy is due to the presence of imaginary frequency
instantaneous normal modes which persist into the solid phase. Model
systems are used to show that imaginary frequency instantaneous normal
modes (and even those at the top of the barrier along that mode) are
not necessarily indicators of diffusive barrier crossing as used in
Zwanzig's model. These false barriers are shown to be the cause of
all of the imaginary frequency zero-force modes in the solid as well
as many of the imaginary frequency modes in the high-density
supercooled liquid. We therefore dispute their utility as predictors
of barrier heights or hopping rates in related liquid systems. We
also show that attempts to separate the modes that are truly diffusive
from those with false barriers using a frequency cutoff or local
information on the potential energy surface are not successful at
removing all of the non-barrier modes.