From f5099b7c16f813c073eb8dddb885429ce50e0799 Mon Sep 17 00:00:00 2001 From: lasergyro <31989178+lasergyro@users.noreply.github.com> Date: Tue, 18 Feb 2025 17:15:01 +0100 Subject: [PATCH] Update fix_press_langevin.rst to include full equation 13 Completest equation 13 to match the [paper](https://pubs.aip.org/jcp/article/141/19/194108/152571/Constant-pressure-and-temperature-discrete-time) and the [implementation](https://github.com/lammps/lammps/blob/36e739b734d2c0b12f6ad56da2b3a59a2b35d7b5/src/fix_press_langevin.cpp#L485-L487). --- doc/src/fix_press_langevin.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/doc/src/fix_press_langevin.rst b/doc/src/fix_press_langevin.rst index 8438d72192..02437bd731 100644 --- a/doc/src/fix_press_langevin.rst +++ b/doc/src/fix_press_langevin.rst @@ -54,7 +54,7 @@ the Langevin equation such as: f_P = & \frac{N k_B T_{target}}{V} + \frac{1}{V d}\sum_{i=1}^{N} \vec r_i \cdot \vec f_i - P_{target} \\ Q\ddot{L} + \alpha{}\dot{L} = & f_P + \beta(t)\\ - L^{n+1} = & L^{n} + bdt\dot{L}^{n} \frac{bdt^{2}}{2Q} \\ + L^{n+1} = & L^{n} + bdt\dot{L}^{n} + \frac{bdt^{2}}{2Q} f^{n}_{P} + \frac{bdt}{2Q} \beta^{n+1} \\ \dot{L}^{n+1} = & \alpha\dot{L}^{n} + \frac{dt}{2Q}\left(a f^{n}_{P} + f^{n+1}_{P}\right) + \frac{b}{Q}\beta^{n+1} \\ a = & \frac{1-\frac{\alpha{}dt}{2Q}}{1+\frac{\alpha{}dt}{2Q}} \\ b = & \frac{1}{1+\frac{\alpha{}dt}{2Q}} \\