Sampling the isothermal-isobaric ensemble by Langevin dynamics
- Laboratory of Computational Physics, Huayuan Road 6, Beijing 100088 (China)
- Institute of Applied Physics and Computational Mathematics, Fenghao East Road 2, Beijing 100094 (China)
We present a new method of conducting fully flexible-cell molecular dynamics simulation in isothermal-isobaric ensemble based on Langevin equations of motion. The stochastic coupling to all particle and cell degrees of freedoms is introduced in a correct way, in the sense that the stationary configurational distribution is proved to be consistent with that of the isothermal-isobaric ensemble. In order to apply the proposed method in computer simulations, a second order symmetric numerical integration scheme is developed by Trotter’s splitting of the single-step propagator. Moreover, a practical guide of choosing working parameters is suggested for user specified thermo- and baro-coupling time scales. The method and software implementation are carefully validated by a numerical example.
- OSTI ID:
- 22657847
- Journal Information:
- Journal of Chemical Physics, Vol. 144, Issue 12; Other Information: (c) 2016 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Isobaric-isothermal Monte Carlo simulations from first principles: Application to liquid water at ambient conditions
Generalized rate-code model for neuron ensembles with finite populations