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...
Uloženo v:
| Vydáno v: | The Journal of chemical physics Ročník 122; číslo 9; s. 094114 |
|---|---|
| Hlavní autoři: | , , , |
| Médium: | Journal Article |
| Jazyk: | angličtina |
| Vydáno: |
United States
01.03.2005
|
| ISSN: | 0021-9606 |
| On-line přístup: | Zjistit podrobnosti o přístupu |
| Tagy: |
Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
|
| Shrnutí: | 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. |
|---|---|
| Bibliografie: | ObjectType-Article-1 SourceType-Scholarly Journals-1 ObjectType-Feature-2 content type line 23 |
| ISSN: | 0021-9606 |
| DOI: | 10.1063/1.1858861 |