MateriApps LIVE! 上で、ASE を用いてLAMMPSを実行する
Last Update:2023/07/13
LAMMPSでNaCl分子に応力を加える計算を行う。
- MateriApps LIVE! にASEをインストールする
pip install ase
- LAMMPSrun のページを参考に下記の内容でexamle.pyを作成する
from ase import Atoms, Atom
from ase.calculators.lammpsrun import LAMMPS
a = [6.5, 6.5, 7.7]
d = 2.3608
NaCl = Atoms([Atom('Na', [0, 0, 0]),
Atom('Cl', [0, 0, d])],
cell=a, pbc=True)
calc = LAMMPS()
NaCl.calc = calc
print(NaCl.get_stress())
- ASEを実行する
export ASE_LAMMPSRUN_COMMAND="mpirun lmp"
python3 example.py
実行結果
[0. 0. 0.0004212 0. 0. 0. ]