Operator splitting algorithm for isokinetic SLLOD molecular dynamics

Guoai Pan*, James F. Ely, Clare McCabe, Dennis J. Isbister

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

27 Citations (Scopus)


We apply an operator splitting method to develop a simulation algorithm that has complete analytical solutions for the Gaussian thermostated SLLOD equations of motion [D. J. Evans and G. P. Morriss, Phys. Rev. A 30, 1528 (1984)] for a system under shear. This leads to a homogeneous algorithm for performing both equilibrium and nonequilibrium isokinetic molecular dynamics simulation. The resulting algorithm is computationally efficient. In particular, larger integration time steps can be used compared to simulations with regular Gaussian thermostated SLLOD equations of motion. The utility and accuracy of the algorithm are demonstrated through application to the Weeks-Chandler-Anderson fluid. Although strict conservation of the kinetic energy suppresses thermal fluctuations in the system, this algorithm does not allow simulations at lower shear rates than those normally afforded by older nonequilibrium molecular dynamics simulations.

Original languageEnglish
Article number094114
JournalJournal of Chemical Physics
Issue number9
Publication statusPublished - 1 Mar 2005

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry


Dive into the research topics of 'Operator splitting algorithm for isokinetic SLLOD molecular dynamics'. Together they form a unique fingerprint.

Cite this