计算化学公社

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

作者
Author: Jennie 时间: 2020-6-24 10:50
标题: VASP可以算解析梯度或者数值梯度，在输出文件找到各个原子的梯度吗？
求助！
VASP可以算解析梯度或者数值梯度，在输出文件找到各个原子的梯度吗？

作者
Author: 卡开发发 时间: 2020-6-24 16:41
OUTCAR里面有记录。

作者
Author: Jennie 时间: 2020-6-24 17:58

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

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

作者
Author: 卡开发发 时间: 2020-6-24 19:50

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

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

作者
Author: Jennie 时间: 2020-6-24 19:58

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

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

作者
Author: Jennie 时间: 2020-6-24 20:04

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

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

作者
Author: 卡开发发 时间: 2020-6-24 20:06

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

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

作者
Author: Jennie 时间: 2020-6-24 20:18

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)