NequIP is an open-source code for building E(3)-equivariant interatomic potentials.
PLEASE NOTE: the NequIP code is under active development and is still in beta versions 0.x.x. In general changes to the patch version (the third number) indicate backward compatible beta releases, but please be aware that file formats and APIs may change. Bug reports are also welcomed in the GitHub issues!
NequIP requires:
- Python >= 3.6
- PyTorch >= 1.8, <=1.10.*. PyTorch can be installed following the instructions from their documentation. Note that neither
torchvision
nortorchaudio
, included in the default install command, are needed for NequIP. NequIP is also not currently compatible with PyTorch 1.10; PyTorch 1.9 can be specified withpytorch==1.9
in the install command.
To install:
- We use Weights&Biases to keep track of experiments. This is not a strict requirement โ you can use our package without it โ but it may make your life easier. If you want to use it, create an account here and install the Python package:
pip install wandb
- Install NequIP
git clone https://github.com/mir-group/nequip.git
cd nequip
pip install .
The easiest way to check if your installation is working is to train a toy model:
$ nequip-train configs/minimal.yaml
If you suspect something is wrong, encounter errors, or just want to confirm that everything is in working order, you can also run the unit tests:
pip install pytest
pytest tests/unit/
To run the full tests, including a set of longer/more intensive integration tests, run:
pytest tests/
Note: the integration tests have hung in the past on certain systems that have GPUs. If this happens to you, please report it along with information on your software environment in the Issues!
The best way to learn how to use NequIP is through the Colab Tutorial. This will run entirely on Google Hardware, you will not need to install anything, but can instead simply run it in your browser.
! PLEASE NOTE: the first few training epochs/calls to a NequIP model can be painfully slow. This is expected behaviour as the profile-guided optimization of TorchScript models takes a number of calls to warm up before optimizing the model. This occurs regardless of whether the entire model is compiled because many core components from e3nn are compiled and optimized through TorchScript.
To train a network, you run nequip-train
with a YAML config file that describes your data set, model hyperparameters, and training options.
$ nequip-train configs/example.yaml
A number of example configuration files are provided:
configs/minimal.yaml
: A minimal example of training a toy model on force data.configs/minimal_eng.yaml
: The same, but for a toy model that predicts and trains on only energy labels.configs/example.yaml
: Training a more realistic model on forces and energies. Start here for real models.configs/full.yaml
: A complete configuration file containing all available options along with documenting comments.
Training runs can be restarted using nequip-restart
; training that starts fresh or restarts depending on the existance of the working directory can be launched using nequip-requeue
. All nequip-*
commands accept the --help
option to show their call signatures and options.
The nequip-evaluate
command can be used to evaluate a trained model on a specified dataset, optionally computing error metrics or writing the results to an XYZ file for further processing.
The simplest command is:
$ nequip-evaluate --train-dir /path/to/training/session/
which will evaluate the original training error metrics over any part of the original dataset not used in the training or validation sets.
For more details on this command, please run nequip-evaluate --help
.
The nequip-deploy
command is used to deploy the result of a training session into a model that can be stored and used for inference.
It compiles a NequIP model trained in Python to TorchScript.
The result is an optimized model file that has no dependency on the nequip
Python library, or even on Python itself:
nequip-deploy build path/to/training/session/ where/to/put/deployed_model.pth
For more details on this command, please run nequip-deploy --help
.
An ASE calculator is also provided in nequip.dynamics
.
NequIP is integrated with the popular Molecular Dynamics code LAMMPS which allows for MD simulations over large time- and length-scales and gives users access to the full suite of LAMMPS features.
The interface is implemented as pair_style nequip
. Using it requires two simple steps:
- Deploy a trained NequIP model, as discussed above.
nequip-deploy build path/to/training/session/ path/to/deployed.pth
The result is an optimized model file that has no Python dependency and can be used by standalone C++ programs such as LAMMPS.
- Change the LAMMPS input file to the nequip
pair_style
and point it to the deployed NequIP model:
pair_style nequip
pair_coeff * * deployed.pth <NequIP type for LAMMPS type 1> <NequIP type for LAMMPS type 2> ...
For installation instructions, please see the pair_nequip
repository.
The theory behind NequIP is described in our preprint (1). NequIP's backend builds on e3nn, a general framework for building E(3)-equivariant neural networks (2). If you use this repository in your work, please consider citing NequIP (1) and e3nn (3):
NequIP is being developed by:
- Simon Batzner
- Albert Musaelian
- Lixin Sun
- Anders Johansson
- Mario Geiger
- Tess Smidt
under the guidance of Boris Kozinsky at Harvard.
If you have questions, please don't hesitate to reach out at batzner[at]g[dot]harvard[dot]edu.
If you find a bug or have a proposal for a feature, please post it in the Issues. If you have a question, topic, or issue that isn't obviously one of those, try our GitHub Disucssions.