DeePMD-kit Manual
Table of contents
About DeePMD-kit
DeePMD-kit is a package written in Python/C++, designed to minimize the effort required to build deep learning based model of interatomic potential energy and force field and to perform molecular dynamics (MD). This brings new hopes to addressing the accuracy-versus-efficiency dilemma in molecular simulations. Applications of DeePMD-kit span from finite molecules to extended systems and from metallic systems to chemically bonded systems.
For more information, check the documentation.
Highlights in DeePMD-kit v2.0
- Model compression. Accelerate the efficiency of model inference for 4-15 times.
- New descriptors. Including
se_e2_r
andse_e3
. - Hybridization of descriptors. Hybrid descriptor constructed from concatenation of several descriptors.
- Atom type embedding. Enable atom type embedding to decline training complexity and refine performance.
- Training and inference the dipole (vector) and polarizability (matrix).
- Split of training and validation dataset.
- Optimized training on GPUs.
Highlighted features
- interfaced with TensorFlow, one of the most popular deep learning frameworks, making the training process highly automatic and efficient, in addition Tensorboard can be used to visualize training procedure.
- interfaced with high-performance classical MD and quantum (path-integral) MD packages, i.e., LAMMPS and i-PI, respectively.
- implements the Deep Potential series models, which have been successfully applied to finite and extended systems including organic molecules, metals, semiconductors, and insulators, etc.
- implements MPI and GPU supports, makes it highly efficient for high performance parallel and distributed computing.
- highly modularized, easy to adapt to different descriptors for deep learning based potential energy models.
License and credits
The project DeePMD-kit is licensed under GNU LGPLv3.0.
If you use this code in any future publications, please cite this using
Han Wang, Linfeng Zhang, Jiequn Han, and Weinan E. "DeePMD-kit: A deep learning package for many-body potential energy representation and molecular dynamics." Computer Physics Communications 228 (2018): 178-184.
Deep Potential in a nutshell
The goal of Deep Potential is to employ deep learning techniques and realize an inter-atomic potential energy model that is general, accurate, computationally efficient and scalable. The key component is to respect the extensive and symmetry-invariant properties of a potential energy model by assigning a local reference frame and a local environment to each atom. Each environment contains a finite number of atoms, whose local coordinates are arranged in a symmetry preserving way. These local coordinates are then transformed, through a sub-network, to a so-called atomic energy. Summing up all the atomic energies gives the potential energy of the system.
The initial proof of concept is in the Deep Potential paper, which employed an approach that was devised to train the neural network model with the potential energy only. With typical ab initio molecular dynamics (AIMD) datasets this is insufficient to reproduce the trajectories. The Deep Potential Molecular Dynamics (DeePMD) model overcomes this limitation. In addition, the learning process in DeePMD improves significantly over the Deep Potential method thanks to the introduction of a flexible family of loss functions. The NN potential constructed in this way reproduces accurately the AIMD trajectories, both classical and quantum (path integral), in extended and finite systems, at a cost that scales linearly with system size and is always several orders of magnitude lower than that of equivalent AIMD simulations.
Although being highly efficient, the original Deep Potential model satisfies the extensive and symmetry-invariant properties of a potential energy model at the price of introducing discontinuities in the model. This has negligible influence on a trajectory from canonical sampling but might not be sufficient for calculations of dynamical and mechanical properties. These points motivated us to develop the Deep Potential-Smooth Edition (DeepPot-SE) model, which replaces the non-smooth local frame with a smooth and adaptive embedding network. DeepPot-SE shows great ability in modeling many kinds of systems that are of interests in the fields of physics, chemistry, biology, and materials science.
In addition to building up potential energy models, DeePMD-kit can also be used to build up coarse-grained models. In these models, the quantity that we want to parameterize is the free energy, or the coarse-grained potential, of the coarse-grained particles. See the DeePCG paper for more details.
Download and install
Please follow our GitHub webpage to download the latest released version and development version.
DeePMD-kit offers multiple installation methods. It is recommend using easily methods like offline packages, conda and docker.
One may manually install DeePMD-kit by following the instuctions on installing the Python interface and installing the C++ interface. The C++ interface is necessary when using DeePMD-kit with LAMMPS, i-PI or GROMACS.
Use DeePMD-kit
A quick-start on using DeePMD-kit can be found as follows:
A full document on options in the training input script is available.
Advanced
- Installation
- Data
- Model
- Training
- Freeze and Compress
- Test
- Inference
- Integrate with third-party packages
Code structure
The code is organized as follows:
-
data/raw
: tools manipulating the raw data files. -
examples
: examples. -
deepmd
: DeePMD-kit python modules. -
source/api_cc
: source code of DeePMD-kit C++ API. -
source/ipi
: source code of i-PI client. -
source/lib
: source code of DeePMD-kit library. -
source/lmp
: source code of Lammps module. -
source/gmx
: source code of Gromacs plugin. -
source/op
: tensorflow op implementation. working with library.
Troubleshooting
- Model compatibility
- Installation
- The temperature undulates violently during early stages of MD
- MD: cannot run LAMMPS after installing a new version of DeePMD-kit
- Do we need to set rcut < half boxsize?
- How to set sel?
- How to control the number of nodes used by a job?
- How to tune Fitting/embedding-net size?
Contributing
See DeePMD-kit Contributing Guide to become a contributor!