计算化学公社

标题: PWPB95 不支持计算原子受力情况 [打印本页]

作者
Author: neon 时间: 2023-3-9 08:20
标题: PWPB95 不支持计算原子受力情况
大家好，我计算的体系是含有一个Zn原子的有机体系，在计算单点能的同时希望也能计算原子受力情况，关键词输入PWPB95 Extrapolate(3/4,def2) EnGrad ，报错信息如下：WARNING: Gradient is not yet implemented for the PWPB95, PW6B95 and DSD-PBEB95 funtionals Please choose numerical gradient (! NumGrad)

再次计算时，关键词输入PWPB95 Extrapolate(3/4,def2) NumGrad EnGrad, 报错信息如下:Error (ORCA): Basis set extrapolation not yet implemented in Hessian.

使用其他方法，如 r2scan-3c，可以得到.engrad 文件
想问一下有没有其他方式可以得到PWPB95/CBS 条件下的原子受力情况？

作者
Author: wjc404 时间: 2023-3-9 09:38
可以试试这个链接里面的方法能不能套用到ORCA上http://bbs.keinsci.com/thread-29330-1-1.html。

作者
Author: wzkchem5 时间: 2023-3-9 16:07
计算triple zeta、quadruple zeta的梯度，然后手动套用外推的公式就行了。程序之所以不支持自动外推梯度，是因为外推梯度的应用场景不多，一般来说triple zeta的梯度足够好了。
另外经过检查代码，这个报错信息里的Hessian是写错了（毕竟这个计算不需要Hessian），是想写numerical gradient。这个我已经改过来了

作者
Author: neon 时间: 2023-3-9 20:48
好的，谢谢大家的回复

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