msuruzhon / alchemlyb

the simple alchemistry library

alchemlyb: the simple alchemistry library

alchemlyb makes alchemical free energy calculations easier to do by leveraging the full power and flexibility of the PyData stack. It includes:

1. Parsers for extracting raw data from output files of common molecular dynamics engines such as GROMACS, AMBER, NAMD and other simulation codes.
2. Subsamplers for obtaining uncorrelated samples from timeseries data (including extracting independent, equilibrated samples [Chodera2016] as implemented in the pymbar package).
3. Estimators for obtaining free energies directly from this data, using best-practices approaches for multistate Bennett acceptance ratio (MBAR) [Shirts2008] and BAR (from pymbar) and thermodynamic integration (TI).

Installation

Install via pip from PyPi (alchemlyb)

pip install alchemlyb

or as a conda package from the conda-forge (alchemlyb) channel

conda install -c conda-forge alchemlyb

Update with pip

pip install --update alchemlyb

or with conda run

conda update -c conda-forge alchemlyb

to get the latest released version.

Getting involved

Contributions of all kinds are very welcome.

If you have questions or want to discuss alchemlyb please post in the alchemlyb Discussions.

If you have bug reports or feature requests then please get in touch with us through the Issue Tracker.

We also welcome code contributions: have a look at our Developer Guide. Open an issue with the proposed fix or change in the Issue Tracker and submit a pull request against the alchemistry/alchemlyb GitHub repository.

References

[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.

[Chodera2016]

Chodera, J.D. (2016). A Simple Method for Automated Equilibration Detection in Molecular Simulations. Journal of Chemical Theory and Computation 12, 1799–1805.