Departmental Papers (MSE)

Document Type

Journal Article

Date of this Version



We describe an atomistic method for computing the viscosity of highly viscous liquids based on activated state kinetics. A basin-filling algorithm allowing the system to climb out of deep energy minima through a series of activation and relaxation is proposed and first benchmarked on the problem of adatom diffusion on a metal surface. It is then used to generate transition state pathway trajectories in the potential energy landscape of a binary Lennard-Jones system. Analysis of a sampled trajectory shows the system moves from one deep minimum to another by a process that involves high activation energy and the crossing of many local minima and saddle points. To use the trajectory data to compute the viscosity we derive a Markov Network model within the Green-Kubo formalism and show that it is capable of producing the temperature dependence in the low-viscosity regime described by molecular dynamics simulation, and in the high-viscosity regime (10(2)-10(12) Pa s) shown by experiments on fragile glass-forming liquids. We also derive a mean-field-like description involving a coarse-grained temperature-dependent activation barrier, and show it can account qualitatively for the fragile behavior. From the standpoint of molecular studies of transport phenomena this work provides access to long relaxation time processes beyond the reach of current molecular dynamics capabilities. In a companion paper we report a similar study of silica, a representative strong liquid. A comparison of the two systems gives insight into the fundamental difference between strong and fragile temperature variations.


Copyright 2009 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics.

Reprinted from:
Computing the viscosity of supercooled liquids Akihiro Kushima, Xi Lin, Ju Li, Jacob Eapen, John C. Mauro, Xiaofeng Qian, Phong Diep, and Sidney Yip, J. Chem. Phys. 130, 224504 (2009), DOI:10.1063/1.3139006



Markov processes, molecular dynamics method, potential energy surfaces, silicon compounds, supercooling, viscosity, FINDING SADDLE-POINTS, GLASS-FORMING LIQUIDS, TRANSITION, SURFACES, DIFFUSION, MODELS, STATES



Date Posted: 06 July 2009

This document has been peer reviewed.