Estimates of free energy from work measurements.
All works and free energies are in multiples of kT. (i.e. nats)
Estimate a free energy difference using the Bennett Acceptance Ratio method (BAR).
Three methods for estimating the error are provided. The original ‘BAR’ method [1], which underestimates the error if the posterior is not Gaussian (this happens when there is little overlap between the forward and negative reverse work distributions and therefore the uncertainties are large); the ‘MBAR’ [3] approach which is optimal in the large sample limit, but overestimates the error when the posteriors are not Gaussian, and ‘logistic’, which is the MBAR error, but with a correction for non-overlapping work distributions.
work_f – Measurements of work from forward protocol.
work_r – Measurements of work from reverse protocol.
weights_f – Optional weights for forward works
weights_r – Optional weights for reverse works
uncertainty_method – Method to calculate errors (“BAR” [1], “MBAR”[3], or “Logistic”[4])
Estimated free energy difference, and the estimated error
???, ???, ???, ???
Bayesian free energy estimate
work_f – Measurements of work from forward protocol.
work_r – Measurements of work from reverse protocol.
Posterior mean estimate of the free energy difference, and the estimated error
Calculate the free energy difference from a series of work measurements using the Jarzynski equality.
(Note that we can’t calculate an accurate error bounds for Jarzynski estimate using measurements from only a single protocol directions.)
work_f – Array of shape [N]
Delta free energy
TODO
Calculate the free energy difference from a series of work measurements under the assumption that the work distributions is Gaussian.
work_f – Array of shape [N]
Delta free energy
TODO
BAR for symmetric periodic protocols.
work_ab – Measurements of work from first half of protocol.
work_bc – Measurements of work from mirror image second half of protocol.
uncertainty_method – Method to calculate errors (BAR, MBAR, or Logistic)
Estimated free energy difference to the middle point of the protocol, and an estimated error
The bidirectional Minh-Chodera free energy estimate specialized to a symmetric protocol.
Delta F = (2/N) sum (e^W_ab + e^-W_bc)^-1)
work_ab – Measurements of work from first half of protocol.
work_bc – Measurements of work from mirror image second half of protocol.
Estimate of the free energy