Martyna-Tuckerman-Tobias-Klein barostat
Jump to navigation
Jump to search
|
This article is a 'stub' page, it has no, or next to no, content. It is here at the moment to help form part of the structure of SklogWiki. If you add sufficient material to this article then please remove the {{Stub-general}} template from this page. |
Martyna-Tuckerman-Tobias-Klein barostat [1] [2] has the following equations of motion (Eq.13):
![{\displaystyle {\dot {\mathbf {p} }}_{i}={\mathbf {F} }_{i}-{\frac {{\overline {\mathbf {p} }}_{g}}{W_{g}}}{\mathbf {p} }_{i}-\left({\frac {1}{N_{f}}}\right){\frac {\mathrm {Tr} [{\overline {\mathbf {p} }}_{g}]}{W_{g}}}-{\frac {p_{\xi }}{Q}}{\mathbf {p} }_{i}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/047bebf5a66e401b154a0b8a23ac0d2cad4d422b)
![{\displaystyle {\dot {\overline {\mathbf {p} }}}_{g}=V\left({\overline {\mathbf {p} }}_{\mathrm {int} }-{\overline {\mathbf {I} }}P_{\mathrm {ext} }\right)+\left[{\frac {1}{N_{f}}}\sum _{i=1}^{N}{\frac {{\mathbf {p} }_{i}^{2}}{m_{i}}}\right]{\overline {\mathbf {I} }}-{\frac {p_{\xi }}{Q}}{\overline {\mathbf {p} }}_{g}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/44de86d15a6d2135673a4ca48cc3847978088631)
![{\displaystyle {\dot {p}}_{\xi }=\sum _{i=1}^{N}{\frac {{\mathbf {p} }_{i}^{2}}{m_{i}}}+{\frac {1}{W_{g}}}\mathrm {Tr} \left[{\overline {\mathbf {p} }}_{g}^{t}{\overline {\mathbf {p} }}_{g}\right]-(N_{f}+d^{2})kT}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7f5f62291bbc47e174e19150dbc279935a0ebad6)
References[edit]
- ↑ Glenn J. Martyna, Douglas J. Tobias, and Michael L. Klein "Constant pressure molecular dynamics algorithms", Journal of Chemical Physics 101 pp. 4177-4189 (1994)
- ↑ G. J. Martyna, M. E. Tuckerman, D. J. Tobias and M. L. Klein "Explicit reversible integrators for extended systems dynamics", Molecular Physics 87 pp. 1117-1157 (1996)