An Efficient Integrator Scheme for Sampling the (Quantum) Isobaric–Isothermal Ensemble in (Path Integral) Molecular Dynamics Simulations
Uloženo v:
| Název: | An Efficient Integrator Scheme for Sampling the (Quantum) Isobaric–Isothermal Ensemble in (Path Integral) Molecular Dynamics Simulations |
|---|---|
| Autoři: | Weihao Liang, Sihan Wang, Cong Wang, Weizhou Wang, Xinchen She, Chongbin Wang, Jiushu Shao, Jian Liu |
| Zdroj: | Journal of Chemical Theory and Computation. 21:6394-6409 |
| Publication Status: | Preprint |
| Informace o vydavateli: | American Chemical Society (ACS), 2025. |
| Rok vydání: | 2025 |
| Témata: | Chemical Physics (physics.chem-ph), Statistical Mechanics (cond-mat.stat-mech), Biological Physics (physics.bio-ph), Physics - Chemical Physics, Classical Physics (physics.class-ph), FOS: Physical sciences, Physics - Biological Physics, Physics - Classical Physics, Computational Physics (physics.comp-ph), Physics - Computational Physics, Condensed Matter - Statistical Mechanics |
| Popis: | Because most chemical or biological experiments are performed under conditions of controlled pressure and temperature, it is important to simulate the isobaric-isothermal ensemble at the atomic level to reveal the microscopic mechanism. By extending our configuration sampling protocol for the canonical ensemble, we propose a unified middle scheme to sample the coordinate (configuration) and volume distribution and thereby are able to accurately simulate either classical or quantum isobaric-isothermal processes. Various barostats and thermostats can be employed in the unified middle scheme for simulating real molecular systems with or without holonomic constraints. In particular, we demonstrate the recommended middle scheme by employing the Martyna-Tuckerman-Tobias-Klein barostat and stochastic cell-rescaling barostat, with the Langevin thermostat, in molecular simulation packages (DL_POLY, Amber, Gromacs, etc.). Benchmark numerical tests show that, without additional numerical effort, the middle scheme is competent in increasing the time interval by a factor of 5~10 to achieve the same accuracy of converged results for most thermodynamic properties in (path integral) molecular dynamics simulations. |
| Druh dokumentu: | Article |
| Jazyk: | English |
| ISSN: | 1549-9626 1549-9618 |
| DOI: | 10.1021/acs.jctc.5c00573 |
| DOI: | 10.48550/arxiv.2504.08342 |
| Přístupová URL adresa: | http://arxiv.org/abs/2504.08342 |
| Rights: | STM Policy #29 CC BY |
| Přístupové číslo: | edsair.doi.dedup.....9d19d2139a3e8596f0311d7494d8d2f9 |
| Databáze: | OpenAIRE |
Nájsť tento článok vo Web of Science | Abstrakt: | Because most chemical or biological experiments are performed under conditions of controlled pressure and temperature, it is important to simulate the isobaric-isothermal ensemble at the atomic level to reveal the microscopic mechanism. By extending our configuration sampling protocol for the canonical ensemble, we propose a unified middle scheme to sample the coordinate (configuration) and volume distribution and thereby are able to accurately simulate either classical or quantum isobaric-isothermal processes. Various barostats and thermostats can be employed in the unified middle scheme for simulating real molecular systems with or without holonomic constraints. In particular, we demonstrate the recommended middle scheme by employing the Martyna-Tuckerman-Tobias-Klein barostat and stochastic cell-rescaling barostat, with the Langevin thermostat, in molecular simulation packages (DL_POLY, Amber, Gromacs, etc.). Benchmark numerical tests show that, without additional numerical effort, the middle scheme is competent in increasing the time interval by a factor of 5~10 to achieve the same accuracy of converged results for most thermodynamic properties in (path integral) molecular dynamics simulations. |
|---|---|
| ISSN: | 15499626 15499618 |
| DOI: | 10.1021/acs.jctc.5c00573 |