计算化学公社

标题: Li原子算频率分析报错，度娘无解 [打印本页]

作者
Author: Aristotler 时间: 2016-10-12 17:17
标题: Li原子算频率分析报错，度娘无解
本帖最后由 Aristotler 于 2016-10-12 17:19 编辑

各位老师，同学：如下任务出错，恳求各位帮忙~

gjf文件：
%chk=Li.chk
# freq b3lyp/6-311+g(d,p) test
Li freq
0 2
C-Bq                -1.10420210 0.94209391 -0.03364862
C-Bq                0.29092994 0.93644804 -0.04044060
C-Bq                0.99358445 2.14029676 0.00695796
C-Bq                0.30107251 3.35059926 0.05998261
C-Bq                -1.09372712 3.35618600 0.06629111
C-Bq                -1.79644188 2.15184103 0.01986883
H-Bq                -1.65799338 -0.00716609 -0.07057423
H-Bq                0.83638686 -0.01749265 -0.08186174
H-Bq                0.85530588 4.29966718 0.09728941
H-Bq                -1.63979731 4.30988215 0.10806715
O-Bq                2.42355987 2.13457915 0.00082444
O-Bq                -3.22641328 2.15787552 0.02660656
Li                      4.24565990 2.12729367 -0.00699100

log文件结尾内容：
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Search for a local minimum.
Step number 1 out of a maximum of 2
All quantities printed in internal units (Hartrees-Bohrs-Radians)
Second derivative matrix not updated -- analytic derivatives used.
NTrRo=3 is not 5 or 6.
Error termination via Lnk1e in /vol-th/home/chenjun/gaussian09/g09/l103.exe at Wed Oct 12 16:59:09 2016.
Job cpu time:  0 days  0 hours 38 minutes 13.7 seconds.
File lengths (MBytes):  RWF=    68 Int=    0 D2E=    0 Chk=    5 Scr=    1
附件为全部log文件内容

作者
Author: liyuanhe211 时间: 2016-10-12 17:26
本帖最后由 liyuanhe211 于 2016-10-12 17:29 编辑

不明白你放那么多赝原子只算Li的一个核俩电子想体现什么问题？

P.S. 珍爱生命，远离百度。

作者
Author: Aristotler 时间: 2016-10-12 18:05

liyuanhe211 发表于 2016-10-12 17:26
不明白你放那么多赝原子只算Li的一个核俩电子想体现什么问题？

P.S. 珍爱生命，远离百度。

老师，我看到一个老师的高斯讲义中写到：如果想准确算AB的结合能需要考虑BSSE矫正（不然结合能会变大），解决办法是在A,B各自加对方的Bq，所以尝试一下。附：老师ps部分亮了，哈哈哈

作者
Author: liyuanhe211 时间: 2016-10-12 19:38

Aristotler 发表于 2016-10-12 18:05
老师，我看到一个老师的高斯讲义中写到：如果想准确算AB的结合能需要考虑BSSE矫正（不然结合能会变大）， ...

1. Li离子与你的体系的相互作用应该是较强的，不需要做BSSE。
2. BSSE一般指是对电子能的校正，在单点上做就可以了（虽然对势能面也有影响，见第四条）
3. Gaussian自己带着做counterpoise的功能，不要自己用Bq原子实现
4. Gaussian 用counterpoise时没有解析梯度和频率，不要用它做freq或优化，也没有任何必要
5. 真要带着counterpoise做频率也不是你这个做法。。。

作者
Author: sobereva 时间: 2016-10-12 22:52

Aristotler 发表于 2016-10-12 18:05
老师，我看到一个老师的高斯讲义中写到：如果想准确算AB的结合能需要考虑BSSE矫正（不然结合能会变大）， ...

那个讲义纯粹坑你，仔细看4L写的。
切勿随意相信网上来路不明的，尤其是中文的计算化学资料，绝大多数都有错

搞科研的必须用google，用百度不仅几乎得不到真正有用的东西，还净搜出坑害自己的东西，除非还是本科科研阶段。

作者
Author: Aristotler 时间: 2016-10-13 09:19

sobereva 发表于 2016-10-12 22:52
那个讲义纯粹坑你，仔细看4L写的。
切勿随意相信网上来路不明的，尤其是中文的计算化学资料，绝大多数 ...

谢谢sobereva老师的谆谆教诲，感恩~

作者
Author: Aristotler 时间: 2016-10-13 09:21

liyuanhe211 发表于 2016-10-12 19:38
1. Li离子与你的体系的相互作用应该是较强的，不需要做BSSE。
2. BSSE一般指是对电子能的校正，在单点上 ...

十分感谢老师的回复~~

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