A Stochastic Algorithm for the Isobaric-Isothermal Ensemble with Ewald Summations for All Long Range Forces

Michele Di Pierro*, Ron Elber, Benedict Leimkuhler

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

Abstract / Description of output

We present an algorithm termed COMPEL (COnstant Molecular Pressure with Ewald sum for Long range forces) to conduct simulations in the NPT ensemble. The algorithm combines novel features recently proposed in the literature to obtain a highly efficient and accurate numerical integrator. COMPEL exploits the concepts of molecular pressure, rapid stochastic relaxation to equilibrium, exact calculation of the contribution to the pressure of long-range nonbonded forces with Ewald summation, and the use of Trotter expansion to generate a robust, highly stable, symmetric, and accurate algorithm. Explicit implementation in the MOIL program and illustrative numerical examples are discussed.

Original languageEnglish
Pages (from-to)5624-5637
Number of pages14
JournalJournal of Chemical Theory and Computation
Volume11
Issue number12
Early online date25 Nov 2015
DOIs
Publication statusPublished - 8 Dec 2015

Fingerprint

Dive into the research topics of 'A Stochastic Algorithm for the Isobaric-Isothermal Ensemble with Ewald Summations for All Long Range Forces'. Together they form a unique fingerprint.

Cite this