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

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:The Journal of chemical physics Jg. 122; H. 9; S. 094114
Hauptverfasser: Pan, Guoai, Ely, James F, McCabe, Clare, Isbister, Dennis J
Format: Journal Article
Sprache:Englisch
Veröffentlicht: United States 01.03.2005
ISSN:0021-9606
Online-Zugang:Weitere Angaben
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung: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.
Bibliographie:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 23
ISSN:0021-9606
DOI:10.1063/1.1858861