g09报错Unable to orthonormalize MBS in NAO

ipoy

老师好：
萌新求教！软件g09，超算运行一个半小时才报错，何解？

报错：
Unable to orthonormalize MBS in NAO.
Error termination via Lnk1e in /public/home/mahw/app/g09/l607.exe at Wed Dec  6 17:13:36 2017.
Job cpu time:  0 days  1 hours 28 minutes 31.0 seconds.
----------------------------------------------------------------------------------------------------------------------------------
输入文件：
%chk=DPESIOPRO3_BULI_631mshsgD.chk
#p opt freq=noraman b3lyp/6-31+g(d) nosymm geom=connectivity
pop=(full,nboread)

DPESIOPRO3_BULI_631mshsgD

【分子坐标】

C H Li O Si 0
6-31+g(d)
****

$nbo
bndidx
plot file=DPESIOPRO3_BULI_631mshsgD
$end
-------------------------------------------------------------------------------------------------------------------------------------------------------------------------
请老师帮忙看看！

sobereva

noraman是完全多余的
你既然都要用自定义基组，干嘛还在#里指定基组？
基组去掉弥散再试。如果去了还不行，用更新版本NBO

ipoy

回复错误，见楼下。

ipoy

sobereva 发表于 2017-12-6 19:05
noraman是完全多余的
你既然都要用自定义基组，干嘛还在#里指定基组？
基组去掉弥散再试。如果去了还不行 ...

这代码都是照葫芦画瓢自己瞎写的，接触计算机不多，本人计算机二级水平(很菜很菜)。
看到sob大神回复，我明白noraman和下面自定义的基组是重复的。
sob大神几个问题：
第一：如果改进想要删除norman代码，noraman 代码是写到这吗？noraman b3lyp/6-31+g(d)，我要是删除是删除到这个位置吗？
第二：我要是用noraman，是不是可以把自定义基组全部删除？
C H Li O Si 0
6-31+g(d)
****
这三行完全删除？
第三：我用 b3lyp/6-31g(d)计算过，程序计算没问题，能有效的反馈结果，但电荷分布的结果跟我们预测的不一样，所以我想把弥散加上去再试一试，因为我这里面有若键结合，所以，sob大神是否能推荐个基组。整个体系150个原子一下，有C, H, Li, O, Si, Na.
第四：我用的是学校超算，不知道学校的g09是否能更新。

再次感谢sob大神无私指导。谢谢！！

sobereva

ipoy 发表于 2017-12-6 20:34
这代码都是照葫芦画瓢自己瞎写的，接触计算机不多，本人计算机二级水平(很菜很菜)。
看到sob大神回复， ...

noraman和自定义基组完全无关，看此文
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

Mulliken等原子电荷计算方法，加了弥散，结果完全无意义，看此文：http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml
而对于ADCH、NPA等不怕弥散函数的方法，加了弥散函数对结果基本没影响，加了白加，白费时间，还不容易SCF收敛，还容易造成一些数值问题

ipoy

sobereva 发表于 2017-12-6 20:36
noraman和自定义基组完全无关，看此文
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/3 ...

万分感谢！！！我学习一下。