Bennett Acceptance Ratio

From AlchemistryWiki
Jump to navigation Jump to search



The Bennett Acceptance Ratio (BAR) is one of the earliest free energy methods which draws on data from multiple states to estimate the free energy difference; This method has significantly improved results over EXP. Both EXP and TI require the ensemble average from a single state to estimate free energies. Although TI needs the derivatives at state [math]\displaystyle{ k }[/math], it does not require the configurations from any neighboring state; BAR however, requires configuration information from two states to estimate free energy differences.

BAR works under the principal that at the same configuration, [math]\displaystyle{ \vec{q} }[/math], at two separate states, [math]\displaystyle{ i }[/math] and [math]\displaystyle{ j }[/math], there is a pathway connecting the two potentials and a difference of [math]\displaystyle{ \Delta U_{ij}(\vec{q}) }[/math]. Because the states are in the same configuration, there is a exact relationship between the distributions of potential energy differences [math]\displaystyle{ \Delta U_{ij}(\vec{q}) }[/math] of states sampled from [math]\displaystyle{ i }[/math] and [math]\displaystyle{ \Delta U_{ji}(\vec{q}) }[/math] the distribution of potential energy differences sampled from the state [math]\displaystyle{ j }[/math].[1]. Since it's an exact function of distributions, statistics can be applied to find the optimal way to use the information between two states, improving the free energy estimate. This is where Bennett started his derivation and, since he was the first to derive this, the method was named after him.[2]

Derivation

This derivation starts from a slightly modified version of the core free energy equation. Taking the properties of expectation values, we can write the free energy difference as

[math]\displaystyle{ \Delta A_{ij} = -k_{B} T \ln \frac{Q_j}{Q_i} = k_{B}T \ln \frac{\left\langle\alpha(\vec{q}) \exp[-\beta U_{i}(\vec{q})]\right\rangle_j}{\left\langle\alpha(\vec{q}) \exp[-\beta U_{j}(\vec{q})]\right\rangle_i} }[/math]

which is true for any [math]\displaystyle{ \alpha(\vec{q})\gt 0 }[/math] for all [math]\displaystyle{ \vec{q} }[/math]. This was Bennett's start point and he then used variational calculus to select the value pf [math]\displaystyle{ \alpha(\vec{q}) }[/math] minimizing the variance of the free energy. The end result is an implicit function of the free energy, and the number of samples at each state, [math]\displaystyle{ n_i }[/math] and [math]\displaystyle{ n_j }[/math], and is

[math]\displaystyle{ \sum_{i=1}^{n_i} \frac{1}{1 + \exp(\ln(n_i/n_j) + \beta \Delta U_{ij} - \beta \Delta A))} - \sum_{j=1}^{n_j} \frac{1}{1 + \exp(\ln(n_j/n_i) - \beta \Delta U_{ji} + \beta \Delta A))} = 0 }[/math]

which must be solved numerically. This is the full BAR equation and its full derivation can be found in Bennett's paper.[2] To get the equation above from Bennett's paper, you will have to take the exponential of both sides of his self-consistent Equation 12, then a bit of algebra to get the end result. This equation can also be derived from a maximum likelihood approach.[3] There are a number of other equivalent ways to express this implicit equation.

Comparison with Other Methods

BAR has been shown to be superior to EXP in every way (save maybe simplicity). BAR is not only better from a practical and theoretical aspect,[4][5] and it is shown to converge to EXP in assuming all samples come from one state.[2][3] Even so, less phase overlap is needed to run BAR effectively.

Comparing TI and BAR is not such a simple thing since each approach requires very different information. Based on experience, BAR will do better than TI on average, although more details are needed to make a better comparison. One such detail is that BAR can give in fewer intermediate states the same statistical precision as TI. However, if the integrand is very smooth for TI, it will perform at the same level as BAR (examples of this are changes, not removing, non- and bonded parameters).[3][6] One advantage that could be attributed to BAR over TI is that you do not need to calculate [math]\displaystyle{ \frac{dU}{d\lambda} }[/math] and so you do not have to modify your code to do so.

A variant of BAR which takes in data from more than two states has been developed, it is coincidentally called the Multistate Bennett Acceptance Ratio or MBAR for short.

Estimating Accurate Free Energies with BAR

BAR has a very straightforward method to calculate the free energy difference between two states. It still suffers from the same general problems as the other methods do with regards to variance along the pathway, what parameters are changing, sequence of change etc. and must still be treated with care. BAR can estimate free energies between many states, however, it must do this a pair at a time. As such, BAR requires information collected from [math]\displaystyle{ k-1 }[/math], [math]\displaystyle{ k }[/math], and [math]\displaystyle{ k+1 }[/math] state for each configuration stored.

BAR require an iterative solution but all the samples will be correlated since results are taken from adjacent states. Unlike in TI, however, the results are not so clean to get an error estimate. More indirect methods, such as Bootstrap Sampling are required to get variances. Furthermore, subsampling is needed since each state can have different correlation times.

Download BAR

Although you could write your own BAR code for analysis, there is a free, python-based code available for download at https://simtk.org/home/pymbar. Although this site is mostly for MBAR, there is a BAR implementation bundled with it along with examples.

References

  1. Crooks, G. E. (2000) Path-ensemble averages in systems driven far from equilibrium. Phys. Rev. E 61, 2361–2366. - Find at Cite-U-Like
  2. 2.0 2.1 2.2 Bennett, C. H. (1976) Efficient Estimation of Free Energy differences from Monte Carlo Data. J. Comput. Phys. 22, 245–268. - Find at Cite-U-Like
  3. 3.0 3.1 3.2 Shirts, M. R., Bair, E., Hooker, G., and Pande, V. S. (2003) Equilibrium free energies from nonequilibrium measurements using maximum-likelihood methods. Phys. Rev. Lett 91, 140601. - Find at Cite-U-Like
  4. Shirts, M. R., and Pande, V. S. (2005) Comparison of efficiency and bias of free energies computed by exponential averaging, the Bennett acceptance ratio, and thermodynamic integration. J. Chem. Phys. 122, 144107. - Find at Cite-U-Like
  5. Lu, N. D., Singh, J. K., and Kofke, D. A. (2003) Appropriate methods to combine forward and reverse free-energy perturbation averages. J. Chem. Phys. 118, 2977–2984. - Find at Cite-U-Like
  6. Ytreberg, F. M., Swendsen, R. H., and Zuckerman, D. M. (2006) Comparison of free energy methods for molecular systems. J. Chem. Phys. 125, 184114. - Find at Cite-U-Like