Warning: This library is young. It is not API stable. It is a nucleation point. By all means use and help improve it, but note that it will change with time.
alchemlyb is an attempt to make alchemical free energy calculations easier to do by leveraging the full power and flexibility of the PyData stack. It includes:
- Parsers for extracting raw data from output files of common molecular dynamics engines such as GROMACS [Abraham2015].
- Subsamplers for obtaining uncorrelated samples from timeseries data.
- Estimators for obtaining free energies directly from this data, using best-practices approaches for multistate Bennett acceptance ratio (MBAR) [Shirts2008] and thermodynamic integration (TI).
In particular, it uses internally the excellent pymbar library for performing MBAR and extracting independent, equilibrated samples [Chodera2016].
- Abraham2015
Abraham, M.J., Murtola, T., Schulz, R., Páll, S., Smith, J.C., Hess, B., and Lindahl, E. (2015). GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 1–2, 19–25.
- Chodera2016
Chodera, J.D. (2016). A Simple Method for Automated Equilibration Detection in Molecular Simulations. Journal of Chemical Theory and Computation 12, 1799–1805.
- Shirts2008
Shirts, M.R., and Chodera, J.D. (2008). Statistically optimal analysis of samples from multiple equilibrium states. The Journal of Chemical Physics 129, 124105.