计算化学公社

标题: 计算resp电荷遇到GetVDW: no radius for atom [打印本页]

作者
Author: shaofeng 时间: 2018-10-24 16:09
标题: 计算resp电荷遇到GetVDW: no radius for atom
各位老师，
我在计算resp电荷的时候总是遇到GetVDW:  no radius for atom 的错误。我的输入文件如下，请指教哪里出了问题。

%nproc=14
%mem=30GB
%chk=H3AsO4.chk
# opt=tight hf/6-31g(d) nosymm pop=(mk,readradii) scf=qc iop(6/33=2,6/42=6,6/50=1)

Title

0 1
As                0.05640200 0.07479500 0.04530400
O                -0.07113100 0.18775100 1.74159400
O                1.74375700 -0.01676200 -0.18024800
O                -0.80927200 -1.07316900 -0.65003300
O                -0.28289300 1.67192900 -0.44506800
H                2.03738700 -0.90477700 -0.35653300
H                -0.65658000 1.71005800 -1.31957900
H                -0.89499300 -0.16104600 2.06633400

As 0.58

neutral.gesp

作者
Author: sobereva 时间: 2018-10-24 16:38
仔细看
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441

用Multiwfn计算RESP电荷是最完美、最省事、最灵活强大的，不需要考虑用其它程序。缺半径的时候程序还直接给你推荐半径，贴心得不得了

作者
Author: shaofeng 时间: 2018-10-24 21:17
我想用获得的gesp文件进一步得到力场文件，不知道Multiwfn有没有方便的途径能够实现。

作者
Author: sobereva 时间: 2018-10-25 09:05

shaofeng 发表于 2018-10-24 21:17
我想用获得的gesp文件进一步得到力场文件，不知道Multiwfn有没有方便的途径能够实现。

力场文件和gesp根本毫无直接关系
网上某些文章说要产生gesp（其实是多余的）的目的在于让antechamber调用resp模块产生RESP电荷
Multiwfn不是搞力场参数的，但可以给你RESP电荷
你搞小分子拓扑文件，不让antechamber产生RESP电荷就根本不需要什么gesp文件
“gesp文件进一步得到力场文件”这种描述是明显错误的

作者
Author: XiaoYan 时间: 2019-8-2 08:54
我也遇到了楼主相似的问题，试了很多遍，高斯手册里面提示原子半径默认是埃，末尾的0.58 是错误的

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