This is a Monte Carlo and Molecular Dynamics Simulation software used primarily for gas sorption in crystalline materials. It is a project that began as a re-write and expansion of Massively Parallel Monte Carlo (MPMC), another code developed and maintained by our laboratory, led by Brian Space at the University of South Florida, Dept. of Chemistry, Smart Metal-organic Materials Advanced Research and Technology Transfer (SMMARTT).
-
Download:
git clone https://github.com/khavernathy/mcmdor download .zip file -
Compile:
cd mcmd
cd src
bash compile.sh [ options ]
cd .. -
Run:
./mcmd mcmd.inp
cd mcmd
git pull
cd src
bash compile.sh [ options ]
Take a look at mcmd/src/compile.sh for different options in compilation (OS-specific, CUDA implementation, OpenMP, optimization on different HPC systems, etc.)
You can find details on available options, built-in potentials, etc. on the wiki page: https://github.com/khavernathy/mcmd/wiki
Douglas Franz: dfranz@mail.usf.edu University of South Florida Dept. of Chemistry
-> Monte Carlo simulation in NPT, NVT, NVE, and μVT ensembles.
-> Molecular Dynamics simulation in NVT, NVE, and μVT ensembles.
-> A crystal builder to create fully parameterized supercells from unit cells.
-> A fragment creator based around uniquely named atoms.
-> A LAMMPS input file exporter.
-> Trajectories and restart files in various formats.
-> Automatic radial distribution calculator
-> Hard-coded molecular models for easy input, including multi-sorbate support
-> Easy system basis parametrization via a, b, c, α, β, γ crystal values, or basis vectors
-> Quick routines for energy/force computation
-> Simulated annealing
-> Any periodic cell is supported for both MC and MD; non-periodic systems also supported.
-> Force-fields available are Lennard-Jones (12-6), Ewald electrostatics, and Thole-Applequist polarization.
-> Sample inputs are included. The program takes just one argument: the input file (which itself usually points to a file containing starting atoms).
The program outputs several quantities of interest:
->Uptake of sorbates in wt%, reduced wt%, cm^3/g, mmol/g and mg/g
->Excess adsorption ratio in mg/g
->Selectivities for multi-sorbate simulation
->Qst (heat of sorption) for sorbate
->Sorbate occupation about some site/atom (g(r))
->Diffusion coefficient and specific heat
->Trajectory and restart files to easily pickup a halted job and visualize simulation
->3D histogram data for visualization of sorbate occupation in a material (density visualization).
->Induced dipole strengths for polarization simulations
MCMD works out-of-the-box on
-> Linux (tested on Ubuntu 16.04)
-> Mac (tested on OS X El Capitan v10.11.6)
-> Windows (tested using Cygwin and Windows 7 with gcc 5.4.0 installed)
-> Raspberry Pi (3, using Raspian OS).
We recommend Visual Molecular Dynamics (VMD) for data visualization, but the output is compatible with most other software, e.g. Avogadro, Molden, Ovito, etc.
Below is a list of scientific publications/presentations that have been facilitated by this software.
-
Franz, D. M., Forrest, K. A., Pham, T., & Space, B. (2016). Accurate H2 Sorption Modeling in the rht-MOF NOTT-112 Using Explicit Polarization. Crystal Growth & Design, 16(10), 6024-6032.
-
Luna, I. B., Franz, D. M. & Space, B. (2017). Theoretical Investigation of Hydrogen Sorption in Metal-organic Framework NOTT-101. 2017 REU-RET Poster Symposium, University of South Florida, Tampa, FL.
-
Pham, T., Forrest, K. A., Franz, D. M., Guo, Z., Chen, B., & Space, B. (2017). Predictive models of gas sorption in a metal–organic framework with open-metal sites and small pore sizes. Physical Chemistry Chemical Physics, 19(28), 18587-18602.
-
Pham, T., Forrest, K. A., Franz, D. M., & Space, B. (2017). Experimental and theoretical investigations of the gas adsorption sites in rht-metal-organic frameworks. CrystEngComm, 19, 4646-4665.