应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 计算出错
111/2返回列表12下一页
查看: 1606  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

w824194485

木虫 (正式写手)

[求助] 计算出错 已有1人参与

输入文件
%NOSAVE
%mem=500MB
# OPT=gdiis freq CBS-QB3

Title Card Required

0 2
N              
S                  1            B1
O                  2            B2    1            A1

   B1             1.50866118
   B2             1.45832917
   A1           124.00909224
  输出文件
-----------------------------------------------------------
#N Geom=AllCheck Guess=TCheck SCRF=Check CCSD(T)/6-31+G(d')
-----------------------------------------------------------
1/18=20,19=11,29=7,38=1,40=1/1;
2/12=2,40=1/2;
3/5=11,6=6,7=11,11=9,14=-4,16=1,25=1,30=1,70=2,116=-2/1,2,3;
4/5=101/1;
5/5=2/2;
8/6=4,9=120000,10=1/1,4;
9/5=7,14=2/13;
6/7=2,8=2,9=2,10=2/1;
99/5=1,9=1/99;
-------------------
Title Card Required
-------------------
Redundant internal coordinates taken from checkpoint file:
/data/tmp/Gau-29637.chk
Charge =  0 Multiplicity = 2
N,0,-0.0063063572,0.,-0.0065648995
S,0,0.003093632,0.,1.5051797991
O,0,1.2120929758,0.,2.3343866301
Recover connectivity data from disk.
                          Input orientation:                          
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
      1          7           0       -0.006306    0.000000   -0.006565
      2         16           0        0.003094    0.000000    1.505180
      3          8           0        1.212093    0.000000    2.334387
---------------------------------------------------------------------
                    Distance matrix (angstroms):
                    1          2          3
     1  N    0.000000
     2  S    1.511774   0.000000
     3  O    2.639044   1.466037   0.000000
Stoichiometry    NOS(2)
Framework group  CS[SG(NOS)]
Deg. of freedom     3
Full point group                 CS      NOp   2
Largest Abelian subgroup         CS      NOp   2
Largest concise Abelian subgroup C1      NOp   1
                         Standard orientation:                        
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
      1          7           0        1.401806   -0.231005    0.000000
      2         16           0        0.000000    0.335034    0.000000
      3          8           0       -1.226580   -0.467938    0.000000
---------------------------------------------------------------------
Rotational constants (GHZ):     68.7554633      9.7059073      8.5052574
Standard basis: 6-31+(d') (6D, 7F)
There are    45 symmetry adapted basis functions of A'  symmetry.
There are    16 symmetry adapted basis functions of A"  symmetry.
Integral buffers will be    131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
    61 basis functions,   120 primitive gaussians,    61 cartesian basis functions
    16 alpha electrons       15 beta electrons
       nuclear repulsion energy        96.6357984684 Hartrees.
NAtoms=    3 NActive=    3 NUniq=    3 SFac= 7.50D-01 NAtFMM=   80 NAOKFM=F Big=F
One-electron integrals computed using PRISM.
NBasis=    61 RedAO= T  NBF=    45    16
NBsUse=    61 1.00D-06 NBFU=    45    16
Initial guess read from the checkpoint file:  /data/tmp/Gau-29637.chk
B after Tr=     0.000000    0.000000    0.000000
         Rot=    1.000000    0.000000    0.000000    0.000000 Ang=   0.00 deg.
Initial guess orbital symmetries:
Alpha Orbitals:
       Occupied  (A') (A') (A') (A') (A') (A") (A') (A') (A') (A')
                 (A') (A") (A') (A') (A") (A')
       Virtual   (A") (A') (A') (A") (A') (A') (A') (A") (A') (A")
                 (A') (A') (A") (A') (A') (A') (A") (A') (A') (A')
                 (A") (A') (A') (A') (A") (A') (A') (A") (A") (A')
                 (A") (A') (A') (A") (A') (A") (A') (A') (A') (A')
                 (A') (A') (A') (A') (A')
Beta  Orbitals:
       Occupied  (A') (A') (A') (A') (A') (A") (A') (A') (A') (A')
                 (A") (A') (A') (A") (A')
       Virtual   (A') (A") (A') (A') (A") (A') (A') (A') (A") (A')
                 (A") (A') (A') (A") (A') (A') (A") (A') (A') (A')
                 (A') (A") (A') (A') (A') (A") (A') (A") (A') (A")
                 (A') (A") (A') (A') (A") (A') (A") (A') (A') (A')
                 (A') (A') (A') (A') (A') (A')
Initial guess <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7644 S= 0.5072
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Keep R1 and R2 ints in memory in canonical form, NReq=4427933.
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UHF) =  -526.685148166     A.U. after  129 cycles
             Convg  =    0.3454D-05             -V/T =  2.0011
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.9197 S= 0.5815
<L.S>= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation     0.9197,   after     0.7571
Convergence failure -- run terminated.
Error termination via Lnk1e in /opt/g09/l502.exe at Tue Oct  6 13:37:29 2015.
回复此楼

» 猜你喜欢
我就是我
1楼 2015-10-06 14:45:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

小范范1989

木虫 (著名写手)
http://sobereva.com/61
赞一下(1人) 回复此楼
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
2楼2015-10-06 15:01:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)
引用回帖:
3楼: Originally posted by w824194485 at 2015-10-06 15:14:52
小白一个，求直接指点怎么弄...

你看看帖子就知道了，就是你的计算没有收敛，
这不是你采用了opt=gdiis的方式来调节收敛了，
不收敛的问题分为两个：几何优化不收敛和scf不收敛。
我的建议是：（个人见解，不一定对哈）
1：opt=(maxstep=10,notrust) 这个是因为默认补偿是0.03bohr，我们调节为0.01，减小步长，看看能否成功。
2：可以调节对称性： symm=loose
3：我曾经成功的方式：先用小的机组，再用大的机组。
4： freq和opt是否可以分开计算？
5：使用maxcyc来调节步数，基本不管用。

以上是个人见解，不一定对，关于两个收敛的调节方式有很多帖子，楼主可以多看看帖子，都是从不会，到慢慢熟悉的，我的感觉是这样哈
赞一下(1人) 回复此楼
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
7楼2015-10-06 17:33:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

w824194485

木虫 (正式写手)
引用回帖:
2楼: Originally posted by 小范范1989 at 2015-10-06 15:01:59
http://sobereva.com/61

小白一个，求直接指点怎么弄
赞一下 回复此楼
我就是我
3楼2015-10-06 15:14:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongma2008

木虫 (著名写手)

论文？苦！主。

【答案】应助回帖

感谢参与，应助指数 +1
SCF Done:  E(UHF) =  -526.685148166     A.U. after  129 cycles
系统默认128步，增大循环个数即可。maxcyc 参数
赞一下 回复此楼
4楼2015-10-06 15:26:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w824194485

木虫 (正式写手)
引用回帖:
4楼: Originally posted by yongma2008 at 2015-10-06 15:26:04
SCF Done:  E(UHF) =  -526.685148166     A.U. after  129 cycles
系统默认128步，增大循环个数即可。maxcyc 参数

不行哦
回复此楼
我就是我
5楼2015-10-06 15:43:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w824194485

木虫 (正式写手)
大神呢，求解啊
回复此楼
我就是我
6楼2015-10-06 16:15:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w824194485

木虫 (正式写手)
引用回帖:
7楼: Originally posted by 小范范1989 at 2015-10-06 17:33:17
你看看帖子就知道了，就是你的计算没有收敛，
这不是你采用了opt=gdiis的方式来调节收敛了，
不收敛的问题分为两个：几何优化不收敛和scf不收敛。
我的建议是：（个人见解，不一定对哈）
1：opt=(maxstep=10,n ...

1.2.4都试了还是不行，只算OPT不算频率是同样的错误，这个结构是我先用B3LYP6-311+g(3df)算过的，你说先用小基组，但是我不知道比CBS-QB3小的基组是什么
赞一下 回复此楼
我就是我
8楼2015-10-06 18:34:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w824194485

木虫 (正式写手)
引用回帖:
2楼: Originally posted by 小范范1989 at 2015-10-06 15:01:59
http://sobereva.com/61

用帖子的方法解决了，谢谢
赞一下 回复此楼
我就是我
9楼2015-10-06 19:04:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)
引用回帖:
9楼: Originally posted by w824194485 at 2015-10-06 19:04:23
用帖子的方法解决了，谢谢...

请问楼主是采用哪条解决的？谢谢指点
赞一下 回复此楼
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
10楼2015-10-06 19:48:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 w824194485 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +3 gouxfjh 2026-08-28 3/150 2026-08-29 01:01 by jurkat.1640
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
关闭
请填处理意见
关闭