mindspore/model_zoo/research/hpc/molecular_dynamics

• Contents
• • Description
  • Model Architecture
  • Dataset
  • Environment Requirements
  • Script Description
  • • Script and Sample Code
    • Training Process
    • Evaluation Process
    • Result
  • ModelZoo Homepage

Contents

• Description
• Model Architecture
• Dataset
• Environment Requirements
• Script Description
  • Script and Sample Code
  • Training Process
  • Evaluation Process
  • Result
• ModelZoo Homepage

Description

Molecular Dynamics (MD) is playing an increasingly important role in the research of biology, pharmacy, chemistry, and materials science. The architecture is based on DeePMD, which using an NN scheme for MD simulations, which overcomes the limitations associated to auxiliary quantities like the symmetry functions or the Coulomb matrix. Each environment contains a number of atoms, whose local coordinates are arranged in a symmetry preserving way following the prescription of the Deep Potential method. According to the atomic position, atomic types and box tensor to construct energy, force and virial.

Thanks a lot for DeePMD team's help.

[1] Paper: L Zhang, J Han, H Wang, R Car, W E. Deep potential molecular dynamics: a scalable model with the accuracy of quantum mechanics. Physical review letters 120 (14), 143001 (2018).

[2] Paper: H Wang, L Zhang, J Han, W E. DeePMD-kit: A deep learning package for many-body potential energy representation and molecular dynamics. Computer Physics Communications 228, 178-184 (2018).

Model Architecture

The overall network architecture of MD simulation is show below.

Link

Dataset

Dataset used: deepmodeling/deepmd-kit/examples/water/data

The data is generated by Quantum Espresso and the input of Quantum Espresso is setted manually.

The directory structure of the data is as follows:

└─data
    ├─type.raw
    ├─set.000
    │   ├──box.npy
    │   ├──coord.npy
    │   ├──energy.npy
    │   └──force.npy
    ├─set.001
    ├─set.002
    └─set.003

In deepmodeling/deepmd-kit/source:

• Use train/DataSystem.py to get coord and atype.
• Use function compute_input_stats in train/DataSystem.py to get avg and std.
• Use op/descrpt_se_a.cc to get nlist.
• Save coord, atype, avg, std and nlist as Npz file for infer.

Environment Requirements

• Hardware (Ascend)
• Framework
  • MindSpore
• For more information, please check the resources below:
  • MindSpore Tutorials
  • MindSpore Python API

Script Description

Script and Sample Code

├── md
    ├── README.md               # descriptions about MD
    ├── script
    │   ├── eval.sh             # evaluation script
    ├── src
    │   ├── descriptor.py       # descriptor function
    │   ├── virial.py           # calculating virial function
    │   └── network.py          # MD simulation architecture
    └── eval.py                 # evaluation interface

Training Process

To Be Done

Evaluation Process

After installing MindSpore via the official website, you can start evaluation as follows:

python eval.py --dataset_path [DATASET_PATH] --checkpoint_path [CHECKPOINT_PATH]

checkpoint can be trained by using DeePMD-kit, and convert into the ckpt of MindSpore.

Result

推理的结果如下：

atom_ener: -94.38766   -94.294426  -94.39194   -94.70758   -94.51311   -94.457954 ...
force: 1.64911175 -1.09822524  0.46055657 -1.34915102 -0.33827361 -0.97184098 ...
virial: -11.736662   -4.286214    2.8852937  -4.286209  -10.408775   -5.6738234 ...

ModelZoo Homepage

Please check the official homepage.