Skip to content

tomswinburne/tammber

Repository files navigation

TAMMBER

Temperature Accelerated Markov Model construction with Bayesian Estimation of Rates

TAMMBER is a heavily modified variant of the ParSplice code (see below).

Designed for massively parallel deployment, TAMMBER optimally manages thousands of workers performing accelerated molecular dynamics (TAD) and minimum energy path routines (NEB). This sampling data is constantly collated and processed using Bayesian techniques to yeild uncertainty-controlled kMC/Markov models of complex atomistic target systems, with minimal end-user involvement. Recent developements allow for the autonomous construction and convergence of arbitrarily complex diffusion tensors.

To be efficient at the peta- or exa-scale requires optimal worker management. TAMMBER achieves this by treating as-yet-unseen configuration space as an absorbing sink, then calculating (given the current information) where additional sampling effort will maximally increase the time-to-absorbtion through discovery or reduction of uncertainty.

TAMMBER efficiently scales to 100,000+ cores, though can be profitibly employed on a few hundred. See the publications below for more detail.

Questions / bugs ? Raise a github issue or email thomas "dot" swinburne "at" cnrs "dot" fr

User Guide

  • If you are using TAMMBER for diffusion studies please read the "Cluster Definitions For Diffusion" section!
  • It is highly recommended to run the tammber-md test routine first!

- Output Analysis : Run online with Binder: Binder

  • If the Binder server is slow the notebook can be viewed here, or run the notebook locally after downloading this repository.

Coming Soon

  • Updated analysis scripts with custom point groups and disconnectivity graphs
  • More detailed analysis of NEB simulations "on-the-fly" (specifically for multiple minima)

Publications:

[1] Thomas D Swinburne and Danny Perez, Self-optimized construction of transition rate matrices from accelerated atomistic simulations with Bayesian uncertainty quantification, Physical Review Materials 2018, preprint

Original publication detailing the optimal worker management and Bayesian uncertainty quantification

[2] Thomas D Swinburne and Danny Perez, Automated Calculation of Defect Transport Tensors, NPJ Computational Materials, 2020. article

Incorporation of crystal symmetries to efficiently and autonomously calculate diffusion tensors. See the Getting Started section


ParSplice code

The ParSplice code implements the Parallel Trajectory Splicing algorithm described in [3]. This method is part of the Accelerated Molecular Dynamics family of techniques developed in Los Alamos National Laboratory over the last 16 years. These methods aim at generating high-quality trajectories of ensembles of atoms in materials. ParSplice uses multiple independent replicas of the system in order to parallelize the generation of such trajectories in the time domain, enabling simulations of systems of modest size over very long timescales. ParSplice includes capabilities to store configurations of the system, to generate and distribute tasks across a large number of processors, and to harvest the results of these tasks to generate long trajectories. ParSplice is a management layer that orchestrate large number of calculations, but it does not perform the actual molecular dynamics itself; this is done by external molecular dynamics engines.

[3] Danny Perez, Ekin D Cubuk, Amos Waterland, Efthimios Kaxiras, Arthur F Voter, Long-time dynamics through parallel trajectory splicing, Journal of chemical theory and computation 12, 18 (2015)

About

TAMMBER fork of ParSplice - Accelerated, massively parallel construction of Markov/kMC Models from MD

Resources

License

Stars

Watchers

Forks

Packages

No packages published