Faster strain fluctuation methods through partial volume updates
J. Chem. Phys. 130, 194706 (2009); doi:10.1063/1.3122383
Published 19 May 2009
You are not logged in to this journal. Log in
Elastic systems that are spatially heterogeneous in their mechanical response pose special challenges for molecular simulations. Standard methods for sampling thermal fluctuations of a system's size and shape proceed through a series of homogeneous deformations, whose magnitudes can be severely restricted by its stiffest parts. Here we present a Monte Carlo algorithm designed to circumvent this difficulty, which can be prohibitive in many systems of modern interest. By deforming randomly selected subvolumes alone, it naturally distributes the amplitude of spontaneous elastic fluctuations according to intrinsic heterogeneity. We describe in detail implementations of such “slice moves” that are consistent with detailed balance. Their practical application is illustrated for crystals of two-dimensional hard disks and random networks of cross-linked polymers.
©2009 American Institute of Physics
| History: | Received 13 October 2008; accepted 31 March 2009; published 19 May 2009 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/130/194706/1 |
REFERENCES (22)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- M. Parrinello and A. Rahman, Phys. Rev. Lett. 45, 1196 (1980).
- J. Ray and A. Rahman, J. Chem. Phys. 80, 4423 (1984).
- J. Ray and A. Rahman, J. Chem. Phys. 82, 4243 (1985).
- R. Najafabadi and S. Yip,
Scr. Metall. 17, 1199 (1983) . - J. Lill and J. Broughton, Phys. Rev. B 49, 11619 (1994).
- K. W. Wojciechowski and A. C. Brańka,
Phys. Lett. A 134, 314 (1989) . - A. C. Brańka and K. W. Wojciechowski,
Mol. Phys. 78, 1513 (1993) . - M. L. Gardel, J. H. Shin, F. C. MacKintosh, L. Mahadevan, P. Matsudaira, and D. A. Weitz,
Science 304, 1301 (2004) . - J. Wilhelm and E. Frey, Phys. Rev. Lett. 91, 108103 (2003).
- D. A. Head, F. C. MacKintosh, and A. J. Levine, Phys. Rev. E 68, 025101 (2003).
- For NPT simulations we are not restricted to choosing a rectangular box, but can choose any unit cell that is scaled by a vector li; the rectangular case was chosen here for notational convenience.
- D. Frenkel and B. Smit, Understanding Molecular Simulation, 2nd ed. (Academic, London, 2002).
- Note that the ensemble of box deformations at constant applied isotropic pressure differs from that at fixed thermodynamic tension. Computing elastic constants from strain fluctuations at constant pressure thus requires a careful accounting of contributions to corresponding compliances from external forces (Refs. 2,6).
- M. A. Bates and D. Frenkel, Phys. Rev. E 61, 5223 (2000).
- S. Sengupta, P. Nielaba, and K. Binder, Phys. Rev. E 61, 6294 (2000).
- K. W. Wojciechowski, K. V. Tretiakov, A. C. Brańka, and M. Kowalik, J. Chem. Phys. 119, 939 (2003).
- D. R. Nelson, Defects and Geometry in Condensed Matter (Cambridge University Press, Cambridge, 2002).
- P. M. Chaikin and T. C. Lubensky, Principles of Condensed Matter Physics (Cambridge University Press, Cambridge, 1995).
- B. J. Alder and T. E. Wainwright,
Phys. Rev. 127, 359 (1962) . - J. Marko and E. Siggia,
Macromolecules 28, 8759 (1995) . - S. Whitelam and P. Geissler, J. Chem. Phys. 127, 154101 (2007).
- R. Swendsen and J. Wang, Phys. Rev. Lett. 58, 86 (1987).








