docs/reference-manual/expanded-ensemble.rst

   1 Expanded Ensemble
   2 -----------------
   3
   4 In an expanded ensemble simulation \ :ref:`68 <refLyubartsev1992>`, both the
   5 coordinates and the thermodynamic ensemble are treated as configuration
   6 variables that can be sampled over. The probability of any given state
   7 can be written as:
   8
   9 .. math:: P(\vec{x},k) \propto \exp\left(-\beta_k U_k + g_k\right),
  10
  11 where :math:`\beta_k = \frac{1}{k_B T_k}` is the :math:`\beta`
  12 corresponding to the :math:`k`\ th thermodynamic state, and :math:`g_k`
  13 is a user-specified weight factor corresponding to the :math:`k`\ th
  14 state. This space is therefore a *mixed*, *generalized*, or *expanded*
  15 ensemble which samples from multiple thermodynamic ensembles
  16 simultaneously. :math:`g_k` is chosen to give a specific weighting of
  17 each subensemble in the expanded ensemble, and can either be fixed, or
  18 determined by an iterative procedure. The set of :math:`g_k` is
  19 frequently chosen to give each thermodynamic ensemble equal probability,
  20 in which case :math:`g_k` is equal to the free energy in non-dimensional
  21 units, but they can be set to arbitrary values as desired. Several
  22 different algorithms can be used to equilibrate these weights, described
  23 in the mdp option listings.
  24
  25 In |Gromacs|, this space is sampled by alternating sampling in the
  26 :math:`k` and :math:`\vec{x}` directions. Sampling in the
  27 :math:`\vec{x}` direction is done by standard molecular dynamics
  28 sampling; sampling between the different thermodynamics states is done
  29 by Monte Carlo, with several different Monte Carlo moves supported. The
  30 :math:`k` states can be defined by different temperatures, or choices of
  31 the free energy :math:`\lambda` variable, or both. Expanded ensemble
  32 simulations thus represent a serialization of the replica exchange
  33 formalism, allowing a single simulation to explore many thermodynamic
  34 states.