エネルギー最小化with reaxFF

何もかも手探りで面白い。lammps配布サイトにあるドキュメントを読みながら一つ一つ理解していこう。まずはエネルギー最小化ができたのでメモ。これが正解かどうかは保証しない。

# minimize graphene
# .....
units real #単位系を実単位系に

atom_style charge
read_data data.graphene #座標データを読み込む

pair_style reax/c lmp_control　#力場の指定
pair_coeff * * ffield.reax.graphene 2 1 3　#力場の指定

neighbor 2 bin #隣接原子の計算に関する設定
neigh_modify every 10 delay 0 check no

fix 1 all nve
fix 2 all qeq/reax 1 0.0 10.0 1e-6 param.qeq

#Run minimization
reset_timestep 0
thermo 10
thermo_style custom step pe lx ly lz press pxx pyy pzz c_eatoms evdwl ecoul　#書き出したいデータの指定

min_style cg　#最小化アルゴリズムを共益勾配法にする
minimize 1e-25 1e-25 50000 100000　#最適化実行

#最適化後の座標書き出し
dump 1 all atom 1 final.dump
run 0
undump 1

ユニットセルの最小化もしたいのだけれど、
reaxFFでは応力計算の部分が実装されてないので、
現状では無理みたいだ。こんなエラーがでる。

WARNING: Pair reax/c cannot yet compute per-atom energy or stress

知ってる方がいたら助けて下さい。