通过其他过程得到的结构,例如压缩、拉伸,或一定温度压强下弛豫完成等等,得到所需的结构,这里只演示基于该结构如何计算应力。应力原则上跟温度、压强没有直接关系,仅仅由结构(当然力场也会影响结果)决定。温度、压强等是通过影响结构,而产生了影响应力的效果。因此Task设置为Single Point即可,不需要设置温度压强。
本文以ReaxFF作为范例,其他模块例如机器学习势、DFTB、基于BAND模块则属于AIMD,除了Main的参数设置不同,其他设置完全一致。
SCM → Output,在底部搜索框搜索“Stress Tensor”:
Stress Tensor (hartree/bohr^3) Index 1 2 3 1 0.0000828213 0.0000031971 0.0000009257 2 0.0000031971 0.0000926634 0.0000007136 3 0.0000009257 0.0000007136 0.0001368110
可以看到3×3应力张量矩阵,单位为Hartree/Bohr3,1 Hartree/Borh3 = 2.94315 x 1013 Pa,可见应力的原子单位是一个很大的量纲。