thermoextrap: Thermodynamic Extrapolation/Interpolation Library

This repository contains code used and described in references ^{1 2}.

Overview

If you find this code useful in producing published works, please provide an appropriate citation. Note that the second citation is focused on adding features that make use of GPR models based on derivative information produced by the core code base. For now, the GPR code, along with more information, may be found under here. In a future release, we expect this to be fully integrated into the code base rather than a standalone module.

Code included here can be used to perform thermodynamic extrapolation and interpolation of observables calculated from molecular simulations. This allows for more efficient use of simulation data for calculating how observables change with simulation conditions, including temperature, density, pressure, chemical potential, or force field parameters. Users are highly encourage to work through the Jupyter Notebooks presenting examples for a variety of different observable functional forms. We only guarantee that this code is functional for the test cases we present here or for which it has previously been applied Additionally, the code may be in continuous development at any time. Use at your own risk and always check to make sure the produced results make sense. If bugs are found, please report them. If specific features would be helpful just let us know and we will be happy to work with you to come up with a solution.

Features

• Fast calculation of derivatives

Status

This package is actively used by the author. Please feel free to create a pull request for wanted features and suggestions!

Quick start

Use one of the following to install thermoextrap:

conda install -c conda-forge thermoextrap

or

pip install thermoextrap

Additional dependencies

To utilize the full potential of thermoextrap, additional dependencies are needed. This can be done via pip by using:

pip install thermoextrap[all]

If using conda, then you'll have to manually install some dependencies. For example, you can run:

conda install bottleneck dask "pymbar>=4.0"

At this time, it is recommended to install the Gaussian Process Regression (GPR) dependencies via pip, as the conda-forge recipes are slightly out of date:

pip install tensorflow tensorflow-probability "gpflow>=2.6.0"

Example usage

import thermoextrap

Documentation

See the documentation for a look at thermoextrap in action.

To have a look at using thermoextrap with Gaussian process regression, look in the gpr and gpr_active_learning directories.

License

This is free software. See LICENSE.

Related work

This package extensively uses the cmomy package to handle central comoments.

Contact

Questions may be addressed to Bill Krekelberg at william.krekelberg@nist.gov or Jacob Monroe at jacob.monroe@uark.edu.

Credits

This package was created using Cookiecutter with the usnistgov/cookiecutter-nist-python template.

Footnotes

1. Extrapolation and Interpolation Strategies for Efficiently Estimating Structural Observables as a Function of Temperature and Density ↩

2. Leveraging Uncertainty Estimates and Derivative Information in Gaussian Process Regression for Expedited Data Collection in Molecular Simulations. In preparation. ↩