VASP可以算解析梯度或者数值梯度，在输出文件找到各个原子的梯度吗？

Jennie · 发表于 Post on 2020-6-24 10:50:51

求助！
VASP可以算解析梯度或者数值梯度，在输出文件找到各个原子的梯度吗？

卡开发发 · 发表于 Post on 2020-6-24 16:41:38

OUTCAR里面有记录。

Jennie · 发表于 Post on 2020-6-24 17:58:42

卡开发发发表于 2020-6-24 16:41
OUTCAR里面有记录。

是那个TOTAL-FORCE吗？那个好像不是解析梯度

卡开发发 · 发表于 Post on 2020-6-24 19:50:30

Jennie 发表于 2020-6-24 17:58
是那个TOTAL-FORCE吗？那个好像不是解析梯度

这个就是对核坐标导数，是个3*Natom的矩阵。

Jennie · 发表于 Post on 2020-6-24 19:58:21

卡开发发发表于 2020-6-24 19:50
这个就是对核坐标导数，是个3*Natom的矩阵。

比方Ni111晶面吸附HCl,读H、Cl、Ni每个原子的受力，TOTAL-FORCE的结果Ni好像是正确的，但是HCl的梯度值不太对啊

Jennie · 发表于 Post on 2020-6-24 20:04:12

Jennie 发表于 2020-6-24 19:58
比方Ni111晶面吸附HCl,读H、Cl、Ni每个原子的受力，TOTAL-FORCE的结果Ni好像是正确的，但是HCl的梯度值不 ...

是只能读晶体原子的梯度吗？不能读吸附的气相原子的

卡开发发 · 发表于 Post on 2020-6-24 20:06:37

Jennie 发表于 2020-6-24 20:04
是只能读晶体原子的梯度吗？不能读吸附的气相原子的

没有啊，全部晶格内的原子都可以，但是哪些是吸附的或者不是吸附的序号你得想办法自己索引后提取。比较快的话可以考虑对z方向的坐标设置个阈值。

Jennie · 发表于 Post on 2020-6-24 20:18:17

卡开发发发表于 2020-6-24 20:06
没有啊，全部晶格内的原子都可以，但是哪些是吸附的或者不是吸附的序号你得想办法自己索引后提取。比较快 ...

好吧，我算的吸附原子的梯度不太对，我再想想，谢谢你了！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] VASP可以算解析梯度或者数值梯度，在输出文件找到各个原子的梯度吗？