请教如何从orca输出文件.engrad中读取原子受力

yuchanpei · 发表于 Post on 2022-9-20 20:22:56

本帖最后由 yuchanpei 于 2022-9-20 20:25 编辑

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

RAL · 发表于 Post on 2022-9-20 21:56:15

engrad文件有3×N行对应体系梯度，N为体系原子数，顺序是g1x, g1y, g1z, g2x, g2y, g2z, ..., gNx, gNy, gNz

yuchanpei · 发表于 Post on 2022-9-21 09:00:13

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

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

RAL · 发表于 Post on 2022-9-21 09:43:26

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

是的

yuchanpei · 发表于 Post on 2022-9-21 09:56:59

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

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

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

[ORCA] 请教如何从orca输出文件.engrad中读取原子受力