24小时热门版块排行榜

返回列表

lozai

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 13.5
帖子: 41
在线: 49.2小时
虫号: 1999719
注册: 2012-09-14
性别: GG
专业: 半导体材料

[求助] 求助计算以中性分子优化结构下带一个单位正电荷能量

各位虫友，小弟现在遇到一个比较棘手的问题，就是计算以中性分子优化结构为基础求带一个单位正电荷能量，计算结果总是出错，带一个单位负电荷的能量可以计算出来，但是正电荷始终算不出来（结构都已优化好，且结构不能改变）初始使用的基组是B3LYP/6-31G(d)，计算不出后后又尝试了B3LYP/6-311+g(d,p),下面是输入/输出文件，还请大家指点一二。四种比较相近的分子，都是这种情况，其他能量都可以计算出来，但是中性分子优化结构为基础带一个单位正电荷能量始终搞不定。
# b3lyp/6-311+g(d,p)

[No Title]

1 2
C                0.63276700 1.12586400 -0.11620400
N                -0.70828100 1.18890400 -0.06786600
C                -1.31161900 -0.00047200 0.00014400
N                -0.70740100 -1.18939600 0.06824300
C                0.63359200 -1.12536800 0.11677800
N                1.35092400 0.00051200 0.00035200
N                -2.70889500 -0.00100200 0.00006400
……
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -1829.37779363    A.U. after  129 cycles
         Convg  = 0.7260D-03          -V/T =  2.0045
= 0.0000 = 0.0000 = 0.5000 = 0.7568 S= 0.5034
= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation    0.7568, after    0.7500
DiagDN has N=  1134 LTot=  3436 but NE2=    3 cannot use DSYEVD.
Convergence failure -- run terminated.
Error termination via Lnk1e in /opt/g09/l502.exe at Tue Dec  4 07:28:15 2012.
Job cpu time:  4 days 15 hours 33 minutes 43.5 seconds.
File lengths (MBytes):  RWF= 1169 Int=    0 D2E=    0 Chk=    30 Scr=    1

回复此楼