Skip to main content

Analysis of lnPi results from TMMC simulation

Project description

Repo Docs PyPI license PyPI version Conda (channel only) Code style: black

tmmc-lnpy

Overview

A package to analyze $\ln \Pi(N)$ data from Transition Matrix Monte Carlo simulation. The main output from TMMC simulations, $\ln \Pi(N)$, provides a means to calculate a host of thermodynamic properties. Moreover, if $\ln \Pi(N)$ is calculated at a specific chemical potential, it can be reweighted to provide thermodynamic information at a different chemical potential

Features

tmmc-lnpy provides a wide array of routines to analyze $\ln \Pi(N)$. These include:

  • Reweighting to arbitrary chemical potential
  • Segmenting $\ln \Pi(N)$ (to identify unique phases)
  • Containers for interacting with several values of $\ln \Pi(N)$ in a vectorized way.
  • Calculating thermodynamic properties from these containers
  • Calculating limits of stability, and phase equilibrium

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

pip install tmmc-lnpy

or

conda install -c conda-forge tmmc-lnpy

Example usage

Note that the distribution name tmmc-lnpy is different than the import name lnpy due to name clashing on pypi.

>>> import numpy as np
>>> import lnpy
>>> import lnpy.examples

>>> ref = lnpy.examples.load_example_lnpimasked("lj_sub")

>>> phase_creator = lnpy.PhaseCreator(nmax=1, ref=ref)
>>> build_phases = phase_creator.build_phases_mu([None])
>>> collection = lnpy.lnPiCollection.from_builder(
...     lnzs=np.linspace(-10, 3, 5), build_phases=build_phases
... )


# Collections are like pandas.Series
>>> collection
<class lnPiCollection>
lnz_0   phase
-10.00  0        [-10.0]
-6.75   0        [-6.75]
-3.50   0         [-3.5]
-0.25   0        [-0.25]
 3.00   0          [3.0]
dtype: object


# Access xarray backend for Grand Canonical properties with `xge` accessor
>>> collection.xge.betaOmega()
<xarray.DataArray 'betaOmega' (lnz_0: 5, phase: 1)> Size: 40B
array([[-2.3245e-02],
       [-6.0370e-01],
       [-1.8552e+02],
       [-1.5447e+03],
       [-2.9580e+03]])
Coordinates:
  * lnz_0    (lnz_0) float64 40B -10.0 -6.75 -3.5 -0.25 3.0
  * phase    (phase) int64 8B 0
    beta     float64 8B 1.372
    volume   float64 8B 512.0
Attributes:
    dims_n:         ['n_0']
    dims_lnz:       ['lnz_0']
    dims_comp:      ['component']
    dims_state:     ['lnz_0', 'beta', 'volume']
    dims_rec:       ['sample']
    standard_name:  grand_potential
    long_name:      $\beta \Omega(\mu,V,T)$

Documentation

See the documentation for a look at tmmc-lnpy in action.

License

This is free software. See LICENSE.

Related work

This package is used for with thermoextrap to analyze thermodynamically extrapolated macro state probability distributions.

Contact

The author can be reached at wpk@nist.gov.

Credits

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

Changelog

Changelog for lnpy

Unreleased

See the fragment files in changelog.d

v0.7.0 — 2024-03-28

Added

  • Added submodule lnpy.combine to combine $\ln\Pi$ from multiple simulations.

v0.6.0 — 2023-08-24

Added

  • Added type hints to most all code. Passing mypy (with strict) and pyright (non-strict).
  • Clean up doc strings in several places.
  • Added nbval testing.
  • Ran linters across all code and notebooks.

v0.5.0 — 2023-07-06

Added

Full set of changes: v0.4.0...0.5.0

v0.4.0 — 2023-05-12

Added

  • Package now available on conda-forge Full set of changes: v0.3.0...0.4.0

Changed

  • Changed examples.load_example_maskddata to examples.load_example_lnpimasked for consistency with other method names.

v0.3.0 — 2023-05-02

Added

  • Added support for python3.11

  • Moved _docstrings -> docstrings to make available

  • Moved from local docfiller to module_utilities.docfiller

  • Moved from local cached module to module-utilities.cached

  • Add support for python3.11

Changed

  • Update package layout
  • New linters via pre-commit
  • Development env now handled by tox

Full set of changes: v0.2.2...0.3.0

v0.2.2 - 2023-04-05

Full set of changes: v0.2.1...v0.2.2

v0.2.1 - 2023-04-04

Full set of changes: v0.2.0...v0.2.1

v0.2.0 - 2023-04-04

Full set of changes: v0.1.5...v0.2.0

v0.1.5 - 2022-09-28

Full set of changes: v0.1.4...v0.1.5

v0.1.4 - 2022-09-26

Full set of changes: v0.1.3...v0.1.4

v0.1.3 - 2022-09-15

Full set of changes: v0.1.2...v0.1.3

v0.1.2 - 2022-09-14

Full set of changes: v0.1.1...v0.1.2

v0.1.1 - 2022-09-13

Full set of changes: v0.1.0...v0.1.1

v0.1.0 - 2022-09-13

Full set of changes: v0.0.1...v0.1.0

v0.0.1 - 2022-09-13

This software was developed by employees of the National Institute of Standards and Technology (NIST), an agency of the Federal Government. Pursuant to title 17 United States Code Section 105, works of NIST employees are not subject to copyright protection in the United States and are considered to be in the public domain. Permission to freely use, copy, modify, and distribute this software and its documentation without fee is hereby granted, provided that this notice and disclaimer of warranty appears in all copies.

THE SOFTWARE IS PROVIDED 'AS IS' WITHOUT ANY WARRANTY OF ANY KIND, EITHER EXPRESSED, IMPLIED, OR STATUTORY, INCLUDING, BUT NOT LIMITED TO, ANY WARRANTY THAT THE SOFTWARE WILL CONFORM TO SPECIFICATIONS, ANY IMPLIED WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE, AND FREEDOM FROM INFRINGEMENT, AND ANY WARRANTY THAT THE DOCUMENTATION WILL CONFORM TO THE SOFTWARE, OR ANY WARRANTY THAT THE SOFTWARE WILL BE ERROR FREE. IN NO EVENT SHALL NIST BE LIABLE FOR ANY DAMAGES, INCLUDING, BUT NOT LIMITED TO, DIRECT, INDIRECT, SPECIAL OR CONSEQUENTIAL DAMAGES, ARISING OUT OF, RESULTING FROM, OR IN ANY WAY CONNECTED WITH THIS SOFTWARE, WHETHER OR NOT BASED UPON WARRANTY, CONTRACT, TORT, OR OTHERWISE, WHETHER OR NOT INJURY WAS SUSTAINED BY PERSONS OR PROPERTY OR OTHERWISE, AND WHETHER OR NOT LOSS WAS SUSTAINED FROM, OR AROSE OUT OF THE RESULTS OF, OR USE OF, THE SOFTWARE OR SERVICES PROVIDED HEREUNDER.

Distributions of NIST software should also include copyright and licensing statements of any third-party software that are legally bundled with the code in compliance with the conditions of those licenses.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

tmmc_lnpy-0.7.0.tar.gz (876.3 kB view details)

Uploaded Source

Built Distribution

If you're not sure about the file name format, learn more about wheel file names.

tmmc_lnpy-0.7.0-py3-none-any.whl (802.5 kB view details)

Uploaded Python 3

File details

Details for the file tmmc_lnpy-0.7.0.tar.gz.

File metadata

  • Download URL: tmmc_lnpy-0.7.0.tar.gz
  • Upload date:
  • Size: 876.3 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/5.0.0 CPython/3.11.8

File hashes

Hashes for tmmc_lnpy-0.7.0.tar.gz
Algorithm Hash digest
SHA256 aca7eaffcd9c77f36db5b440dec2159f8b80cfd38ef4dd01f0117a42351fa68a
MD5 07250e38c855b128362f7e698c01368a
BLAKE2b-256 52941018be82c316b4145b9c714882a9ed8d56cdcdbbdf700b8f13502de73017

See more details on using hashes here.

File details

Details for the file tmmc_lnpy-0.7.0-py3-none-any.whl.

File metadata

  • Download URL: tmmc_lnpy-0.7.0-py3-none-any.whl
  • Upload date:
  • Size: 802.5 kB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/5.0.0 CPython/3.11.8

File hashes

Hashes for tmmc_lnpy-0.7.0-py3-none-any.whl
Algorithm Hash digest
SHA256 7d3844ddd42043361b3eb32303e1c7fba3ab9960aa8715a7e5a1d77ee5bcbd3e
MD5 a75b2252158ff7dd566dc108e10e8154
BLAKE2b-256 8cb3130f6ae0fce17d5e3b717bd240766a6853860ba9005867fb4be41e0cd9db

See more details on using hashes here.

Supported by

AWS Cloud computing and Security Sponsor Datadog Monitoring Depot Continuous Integration Fastly CDN Google Download Analytics Pingdom Monitoring Sentry Error logging StatusPage Status page