计算化学公社

标题: 使用gaussian 09 E.01和antechamber拟合RESP电荷失败求助 [打印本页]

作者
Author: 乐神666 时间: 2019-1-21 16:21
标题: 使用gaussian 09 E.01和antechamber拟合RESP电荷失败求助
使用高斯E.01版本，方法/基组、和关键字选择的是： #HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) iop(6/50=1) opt
高斯总是这个报错，虽然是l602,的错误但却不是缺少原子半径的问题，我曾尝试在我的输入文件（gjf）的最后加上**.gesp的文件名字，但是也是不行。

Free-field routines initialized.
FFRead: read line ""
Blank file name read.
Error termination via Lnk1e in /home/ionlab/software/g09/l602.exe

我用同样的关键词同样的体系，使用高斯B.01版本，就能跑完，但是在我最后用 antechamber想要写出RESP的时候，它或会卡在一步，

Info:the number of the path atoms exceeds MAXPATHATOMNUM(7890000) foratom[71],extend the size and reallocate the memory automatically。

然后也是无法写出RESP。

希望哪位大佬给解答一下，不胜感激。

作者
Author: sobereva 时间: 2019-1-21 16:23
改用Multiwfn计算RESP电荷是最佳选择，比antechamber方便好用百倍，看文中例子两分钟就学会
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）

而且antechamber默认产生的那些Gaussian输入文件里的关键词在现在来看都是坑爹的，谁这年头还用HF/6-31G*优化？而且这年都RESP电荷都应当用DFT结合隐式溶剂模型的波函数来算。还用HF/6-31G*优化和其波函数产生RESP电荷，这年头极大几率会被审稿人狠批。
antechamber计算RESP电荷这点已经完全过时了。

而且什么gesp文件，纯粹也都是多余的

PS：注意求助帖应当在帖子标题明确体现出“提问”、“求助”，给你改了

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