PMASS = [Real]
Default: PMASS = 1000
When running NpT simulations with a Langevin thermostat (MDALGO=3), using the method of Parrinello and Rahman, a fictitious mass (in amu) for the lattice degrees-of-freedom has to be assigned using the PMASS-tag. The friction coefficient for lattice degrees-of-freedom have to be specified (in ps-1) by means of the LANGEVIN_GAMMA_L-tag.
The friction coefficients γ for the atomic degrees-of-freedom are specified using the LANGEVIN_GAMMA-tag.
The optimal setting for PMASS depends very much on the particular system at hand and can be considered as a compromise between two opposing factors: too large values lead to very slow variation of lattice degrees of freedom (and hence the sampling becomes inefficient) while too small value can lead to too large geometric changes in a MD step and hence may cause numerical problems. We strongly recommend to make careful tests with various settings before performing the production run.
Related Tags and Sections
- ↑ M. P. Allen and D. J. Tildesley, Computer simulation of liquids, Oxford university press: New York, 1991.
- ↑ M. Parrinello and A. Rahman, Phys. Rev. Lett. 45, 1196 (1980).
- ↑ M. Parrinello and A. Rahman, J. Appl. Phys. 52, 7182 (1981).