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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 计算出错
111/2返回列表12下一页
查看: 1410  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 321求调剂 +4 wasdssaa 2026-03-26 4/200 2026-03-26 16:25 by 不吃魚的貓
[考研] 考研调剂 +7 呼呼？~+123456 2026-03-24 7/350 2026-03-26 15:24 by 小小麦片
[考研] 求调剂 +3 芦lty 2026-03-25 4/200 2026-03-25 23:25 by 芦lty
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研经历 +6 cccchenso 2026-03-23 6/300 2026-03-25 22:25 by 544594351
[考研] 机械学硕总分317求调剂！！！！ +4 Acaciad 2026-03-25 4/200 2026-03-25 19:59 by hanserlol
[考研] 26考研-291分-厦门大学（085601）-柔性电子学院材料工程专业求调剂 +3 min3 2026-03-24 4/200 2026-03-25 18:22 by xcjcqu
[考研] 0854AI CV方向招收调剂 +4 章小鱼567 2026-03-23 4/200 2026-03-25 17:04 by CoderLoser
[考研] 333求调剂 +6 87639 2026-03-21 11/550 2026-03-25 16:17 by 87639
[考研] 一志愿吉林大学材料与化工303分求调剂 +4 为学666 2026-03-24 4/200 2026-03-25 11:27 by BruceLiu320
[考研] 求调剂323材料与化工 +4 1124361 2026-03-24 4/200 2026-03-25 11:19 by shulmg
[考研] 289材料与化工（085600）B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 070300化学求调剂 +9 苑豆豆 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
信息提示
关闭
请填处理意见
关闭