Editing Berendsen barostat
Jump to navigation
Jump to search
The edit can be undone. Please check the comparison below to verify that this is what you want to do, and then publish the changes below to finish undoing the edit.
Latest revision | Your text | ||
Line 1: | Line 1: | ||
{{Stub-general}} | |||
The '''Berendsen barostat''' <ref>[http://dx.doi.org/10.1063/1.448118 H. J. C. Berendsen, J. P. M. Postma, W. F. van Gunsteren, A. DiNola, and J. R. Haak "Molecular dynamics with coupling to an external bath", Journal of Chemical Physics '''81''' pp. 3684-3690 (1984)]</ref> is a method for controlling the [[pressure]] in a [[molecular dynamics]] simulation. | The '''Berendsen barostat''' <ref>[http://dx.doi.org/10.1063/1.448118 H. J. C. Berendsen, J. P. M. Postma, W. F. van Gunsteren, A. DiNola, and J. R. Haak "Molecular dynamics with coupling to an external bath", Journal of Chemical Physics '''81''' pp. 3684-3690 (1984)]</ref> is a method for controlling the [[pressure]] in a [[molecular dynamics]] simulation. | ||
The Berendsen [[barostats | barostat]] adds an extra term to to the equations of motion which effects the pressure change (Eq. 12): | The Berendsen [[barostats | barostat]] adds an extra term to to the equations of motion which effects the pressure change (Eq. 12): | ||
Line 6: | Line 7: | ||
where <math>P_0</math> is the reference pressure, i.e. the pressure of the external pressure "bath", and <math>P</math> is the instantaneous pressure. | where <math>P_0</math> is the reference pressure, i.e. the pressure of the external pressure "bath", and <math>P</math> is the instantaneous pressure. | ||
<math>\tau_P</math> is a time constant. | <math>\tau_P</math> is a time constant. | ||
Within this scheme the coordinates and the box sides are rescaled every | Within this scheme the coordinates and the box sides are rescaled every so-many steps. Assuming the system is isotropic and within a cubic box the scaling factor <math>\mu</math> is given by (Eq. 20): | ||
:<math> \mu = 1 - \frac{\kappa_T \Delta t}{3\tau_P} (P_0 -P)</math> | :<math> \mu = 1 - \frac{\kappa_T \Delta t}{3\tau_P} (P_0 -P)</math> |