Operator splitting algorithm for isokinetic SLLOD molecular dynamics

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 per...

Full description

Saved in:

Bibliographic Details
Published in:	The Journal of chemical physics Vol. 122; no. 9; p. 094114
Main Authors:	Pan, Guoai, Ely, James F, McCabe, Clare, Isbister, Dennis J
Format:	Journal Article
Language:	English
Published:	United States 01.03.2005
ISSN:	0021-9606
Online Access:	Get more information
Tags:	Add Tag No Tags, Be the first to tag this record!

Description
Summary:	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.
Bibliography:	ObjectType-Article-1 SourceType-Scholarly Journals-1 ObjectType-Feature-2 content type line 23
ISSN:	0021-9606
DOI:	10.1063/1.1858861