计算化学公社

标题: 请教如何从orca输出文件.engrad中读取原子受力 [打印本页]

作者
Author: yuchanpei 时间: 2022-9-20 20:22
标题: 请教如何从orca输出文件.engrad中读取原子受力
本帖最后由 yuchanpei 于 2022-9-20 20:25 编辑

刚开始用ORCA计算体系中原子受力，使用的关键词是：! r2SCAN-3c engrad，了解到受力信息是输出在.engrad文件中，但是不知道文件中的信息与原子受力之间的对应关系。
请问，如果我想得到每个原子的受力信息，该怎么从输出文件中进行提取呢？烦请大家能够帮忙解答一下，万分感谢。

作者
Author: RAL 时间: 2022-9-20 21:56
engrad文件有3×N行对应体系梯度，N为体系原子数，顺序是g1x, g1y, g1z, g2x, g2y, g2z, ..., gNx, gNy, gNz

作者
Author: yuchanpei 时间: 2022-9-21 09:00

RAL 发表于 2022-9-20 21:56
engrad文件有3×N行对应体系梯度，N为体系原子数，顺序是g1x, g1y, g1z, g2x, g2y, g2z, ..., gNx, gNy, gN ...

谢谢回复。
想确认一下，当我需要使用原子受力的时候，用的应该是负梯度对么？万分感谢。

作者
Author: RAL 时间: 2022-9-21 09:43

yuchanpei 发表于 2022-9-21 09:00
谢谢回复。
想确认一下，当我需要使用原子受力的时候，用的应该是负梯度对么？万分感谢。

是的

作者
Author: yuchanpei 时间: 2022-9-21 09:56

RAL 发表于 2022-9-21 09:43
是的

收到，感谢回复，谢谢指导。

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