Controlling Pressure

Overview

Teaching: 10 min
Exercises: 0 min
Questions
  • Why do we want to control the pressure of our MD-simulations?

  • What pressure control algorithms are commonly used?

  • What are the strengths and weaknesses of these common barostats?

Objectives
  • Identify the appropriate method for controlling pressure

Introduction

The role of pressure control algorithms is to keep pressure in the simulation system constant or to apply an external stress to the simulated system.

$ \qquad {P}=\frac{NK_{B}T}{V}+\frac{1}{3V}\langle\sum{r_{ij}F_{ij}}\rangle$

Pressure Control Algorithms

  1. Weak coupling methods
  2. Extended system methods
  3. Stochastic methods
  4. Monte-Carlo methods

1. Weak coupling methods

Berendsen pressure bath coupling.

Downsides:

The time constant for pressure bath coupling is the main parameter of the Berendsen thermostat. The pressure of the system is corrected such that the deviation exponentially decays with a lifetime defined by this constant.

Reference: Molecular dynamics with coupling to an external bath

2. Extended system methods

Parrinello-Rahman barostat

Downsides:

Nosé-Hoover barostat

References: [Hoover, 1986], [Martyna, 1994].

MTTK (Martyna-Tuckerman-Tobias-Klein) barostat.

3. Stochastic methods

Langevin piston pressure control.

Reference: Constant pressure molecular dynamics simulation: The Langevin piston method

MTTK and Langevin barostats produce identical ensembles Langevin barostat oscillates less then MTTK and converges faster due to stochastic collisions and damping.

Comparison of Barostats

Reprinted with permission from Rogge et al. 2015, A Comparison of Barostats for the Mechanical Characterization of Metal−Organic Frameworks, J Chem Theory Comput. 2015;11: 5583-97. doi:10.1021/acs.jctc.5b00748. Copyright 2015 American Chemical Society.

Stochastic Cell Rescaling

Reference: [Bernetti and Bussi (2020)]

4. Monte-Carlo pressure control.

References:

  1. Molecular dynamics simulations of water and biomolecules with a Monte Carlo constant pressure algorithm
  2. Constant pressure hybrid Molecular Dynamics–Monte Carlo simulations

Pitfalls


Selecting barostats in molecular dynamics packages

Thermostat\MD package GROMACS NAMD AMBER
Berendsen pcoupl = Berendsen BerendsenPressure on barostat = 1
Stoch. cell rescaling pcoupl = C-rescale    
Langevin   LangevinPiston on  
Monte-Carlo     barostat = 2
Parrinello-Rahman pcoupl = Parrinello-Rahman    
MTTK pcoupl = MTTK    

Key Points

  • Each barostat or thermostat technique has its own limitations and it is your responsibility to choose the most appropriate method or their combination for the problem of interest.