ISSN:
0192-8651
Keywords:
constraints
;
molecular dynamics
;
Langevin dynamics
;
SHAKE
;
Chemistry
;
Theoretical, Physical and Computational Chemistry
Source:
Wiley InterScience Backfile Collection 1832-2000
Topics:
Chemistry and Pharmacology
,
Computer Science
Notes:
In this article, we present a new LINear Constraint Solver (LINCS) for molecular simulations with bond constraints. The algorithm is inherently stable, as the constraints themselves are reset instead of derivatives of the constraints, thereby eliminating drift. Although the derivation of the algorithm is presented in terms of matrices, no matrix matrix multiplications are needed and only the nonzero matrix elements have to be stored, making the method useful for very large molecules. At the same accuracy, the LINCS algorithm is three to four times faster than the SHAKE algorithm. Parallelization of the algorithm is straightforward. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1463-1472, 1997
Additional Material:
2 Ill.
Type of Medium:
Electronic Resource
Permalink