目录

在一定应力下的结构优化

注意显示坐标轴(View → Axis),隐藏晶格矢量(View → Periodic → Show Lattice Vector),同时显示Cell(View → Periodic → Show Unit Cell),这样便于用户清楚地知道应力的施加方向。

本文以ReaxFF作为范例,其他模块例如机器学习势、DFTB、基于BAND模块则属于AIMD,除了Main的参数设置不同,其他设置完全一致。

参数设置

勾选优化晶格常数按钮,并设置应力,共6个数值,分别为xx、yy、zz、yz、xz、xy,本例中第三个值不为0,表示应力施加在z方向。

其中应力的单位1 a.u. = 1 Hartree/Bohr3 = 2.94315 x 1013 Pa,负值为拉伸,正值为压缩。 保存并提交作业。

如果用户不想换算,也可以在保存作业后生成的脚本里面直接修改,例如:

EngineAddons
    ExternalStress
       StressTensorVoigt 0 0 10 0 0 0 [GPa]
    End
End

修改之后,在AMSJobs窗口直接提交即可

结果查看

SCM → Movie → View → Periodic → Show Unit Cell,显示Cell的形状,隐藏晶格矢量(View → Periodic → Show Lattice Vector),显示坐标轴(View → Axis),便于我们观察在应力下,Cell的变化:

对比第一帧与最后一帧,可以明显看到Cell在z方向被拉长了。File → Update Geometry in Input,将该结构更新到AMSinput窗口中,在Model → Lattice可以看到新的晶格常数(当然在Update之前,应该看看旧的晶格常数,作为对比)。