NOTE: GDPy
is the name of this repository and the corresponding python package name is gdpx
.
NOTE: gdpx
is under active development and has not been released. The APIs are frequently changed and we cannot ensure any
backward compatibility.
gdpx
is a pure python package. Since it does not include any codes that actually perform calculations and training, for example, VASP
and DEEPMD
, you should install them by yourselves.
$ python -m pip install git+https://github.com/hsulab/GDPy.git
$ conda install gdpx -c conda-forge
Documentation: https://gdpx.readthedocs.io (Changed from gdpyx)
GDPy stands for Generating Deep Potential with Python (GDPy/GDP¥), including a set of tools and Python modules to automate the structure exploration and the training for machine learning interatomic potentials (MLIPs).
It mainly focuses on the applications in heterogeneous catalysis. The target systems are metal oxides, supported clusters, and solid-liquid interfaces.
We do not implement any MLIP but offers a unified interface to access. Certain MLIP could not be utilised before corresponding required packages are
installed correctly.The calculations are performed by ase calculators
using either python built-in codes (PyTorch, TensorFlow)
or File-IO based external codes (e.g. lammps).
Supported MLIPs:
MLIPs | Representation | Regressor | Implemented Backend |
---|---|---|---|
eann | (Rescursive) Embedded Atom Descriptor | NN/PyTorch | ASE/Python, ASE/LAMMPS |
deepmd | Deep Potential Descriptors | NN/Tensorflow | ASE/Python, ASE/LAMMPS |
lasp | Atom-Centered Symmetry Functions | NN/LASP | ASE/LASP |
nequip | E(3)-Equivalent Message Passing | NN/PyTorch | ASE/Python, ASE/LAMMPS |
NOTE: We use a modified eann package to train and utilise.
NOTE: Allegro is supported as well through the nequip manager.
Other Potentials: Some potentials besides MLIPs are supported. Force fields or semi-empirical potentials are used for pre-sampling to build an initial dataset. Ab-initio methods are used to label structures with target properties (e.g. total energy, forces, and stresses).
Name. | Description | Backend | Notes |
---|---|---|---|
reax | Reactive Force Field | LAMMPS | |
xtb | Tight Binding | xtb | Under development |
VASP | Plane-Wave Density Functional Theory | VASP | |
CP2K | Density Functional Theory | CP2K |
We take advantage of codes in well-established packages (ASE and LAMMPS) to perform basic minimisation and dynamics. Meanwhile, we have implemented several complicated alogirthms in GDPy itself. | Name | Current Algorithm | Backend |
---|---|---|---|
Molecular Dynamics (md) | Brute-Force/Biased Dynamics | ASE, LAMMPS | |
Evolutionary Global Optimisation (evo) | Genetic Algorithm | ASE/GDPy | |
Basin Hopping | Monte Carlo like Global Optimisation | GDPy | |
Adsorbate Configuration (ads) | Adsorbate Configuration Graph Search | GDPy | |
Reaction Event Exploration (rxn) | Artificial Force Induced Reaction (AFIR) | GDPy | |
Grand Cononical Monte Carlo (gcmc) | Monte Carlo with Variable Composition | GDPy |
There are two kinds of workflows according to the way they couple the expedition and the training. Offline workflow as the major category separates the expedition and the training, which collects structures from several expeditions and then trains the MLIP with the collective dataset. This process is highly parallelised and is usually aimed at a general dataset. Online workflow, a really popular one, adopts an on-the-fly strategy to build a dataset during the expedition, where a new MLIP is trained to continue exploration once new candidates are selected (sometimes only one structure every time!). Thus, it is mostly used to train an MLIP for a particular system.
Type | Supported Expedition |
---|---|
Offline | md, evo, ads, rxn |
Online | md |
Jiayan Xu (jxu15@qub.ac.uk)
under the supervision of Prof. P. Hu at Queen's University Belfast.
GDPy project is under the GPL-3.0 license.