/checkensemble

This software allows users to perform statistical test to determine if a given molecular simulation is consistent with the thermodynamic ensemble it is performed in.

Primary LanguagePythonGNU General Public License v2.0GPL-2.0

checkensemble

This software allows users to perform statistical test to determine if a given molecular simulation is consistent with the thermodynamic ensemble it is performed in.

Users should cite the JCTC paper: Shirts, "M. R. Simple Quantitative Tests to Validate Sampling from Thermodynamic Ensembles", J. Chem. Theory Comput., 2013, 9 (2), pp 909–926, http://dx.doi.org/10.1021/ct300688p

All examples can be run automatically by sourcing 'bash runexamples.sh inside the examples/ directory'

The main distribution consists of the following:

  • checkensemble/checkensemble.py
  • checkensemble/readmdfiles.py
  • checkensemble/timeseries.py
  • examples/analyze-md.py
  • examples/analyze-replica.py

and python code for toy models for validating the main tools:

  • examples/harmonic.py
  • examples/harmonic_pressure.py
  • examples/harmonic_mu.py

Until a proper setup.py is installed, add INSTALLDIR/checkensemble/checkensemble in the PYTHONPATH to run the scripts

Run 'analyze-md.py --help' for a listing and explanation of options of simulation output analysis.

Run 'analyze-replica.py --help' for a listing and explanation of options for temperature replica exchange analysis options

Dependencies: the code requires scipy and numpy. matplotlib is required to generate graphs.

For more data, a number of bins closer to 40 may be more appropriate, but the example data sets are smaller for size reasons, and thus require smaller bin numbers to get statistically meaningful results.

The timeseries.py code is the same as what is included as pymbar.py (at least as of 2012).