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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 逼近收敛，但是就是不收敛，求助高手
51/1返回列表
查看: 2590  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

skysky112211

木虫 (小有名气)

NaNa

[求助] 逼近收敛，但是就是不收敛，求助高手

我有一个分子，计算的话，只有一个指标还差最后一点就收敛了，但是就是不能收敛，请问是怎么回事！
Item                                   Value     Threshold  Converged?
Maximum Force            0.000017     0.000015     NO
RMS     Force                0.000007     0.000010     YES
Maximum Displacement     0.000003     0.000060     YES
RMS     Displacement     0.000001     0.000040     YES

坐标

Standard orientation:                        
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
      1         27           0       -1.066625    0.164072    0.000000
      2         27           0        0.820412   -0.768407    0.000000
      3          6           0       -1.326215    1.925930    0.000000
      4          8           0       -1.411484    3.079329    0.000000
      5          6           0       -0.542787   -0.678856    1.591775
      6          8           0       -0.542787   -1.109606    2.676258
      7          6           0        2.380903    0.070865    0.000000
      8          8           0        3.351191    0.700203    0.000000
      9          6           0       -0.542787   -0.678856   -1.591775
     10          8           0       -0.542787   -1.109606   -2.676258

输入文件
# opt=(noraman,tight,cartesian,maxcycle=5000) freq b3lyp/gen symm
=loose SCF=(maxcycle=5000,qc)
基组
Co 0
S 5 1.00
       270991.000000000     0.000310000
        39734.800000000     0.002420000
         9057.460000000     0.012380000
         2598.210000000     0.048490000
          868.200000000     0.146720000
S 1 1.00
          323.431000000     1.000000000
S 1 1.00
          130.860000000     1.000000000
S 1 1.00
           56.121900000     1.000000000
S 1 1.00
           18.921900000     1.000000000
S 1 1.00
            7.952380000     1.000000000
S 1 1.00
            2.197540000     1.000000000
S 1 1.00
            0.846713000     1.000000000
S 1 1.00
            0.122266000     1.000000000
S 1 1.00
            0.044172000     1.000000000
P 4 1.00
         1636.210000000     0.002960000
          390.903000000     0.023360000
          127.884000000     0.103430000
           49.241300000     0.279540000
P 1 1.00
           20.751200000     1.000000000
P 1 1.00
            9.203680000     1.000000000
P 1 1.00
            3.817790000     1.000000000
P 1 1.00
            1.587620000     1.000000000
P 1 1.00
            0.624660000     1.000000000
P 1 1.00
            0.245800000     1.000000000
P 1 1.00
            0.096700000     1.000000000
D 4 1.00
           44.977400000     0.026270000
           12.569000000     0.141820000
            4.244220000     0.359420000
            1.443300000     0.466130000
D 1 1.00
            0.449965000     1.000000000
D 1 1.00
            0.121900000     1.000000000
****
C 0
S 6 1.00
   4232.610000000    0.002029000
    634.882000000    0.015535000
    146.097000000    0.075411000
     42.497400000    0.257121000
     14.189200000    0.596555000
      1.966600000    0.242517000
S 1 1.00
      5.147700000    1.000000000
S 1 1.00
      0.496200000    1.000000000
S 1 1.00
      0.153300000    1.000000000
P 4 1.00
     18.155700000    0.018534000
      3.986400000    0.115442000
      1.142900000    0.386206000
      0.359400000    0.640089000
P 1 1.00
      0.114600000    1.000000000
D 1  1.00
       0.75000000    1.00000000
****
O 0
S 6 1.00
   7816.540000000    0.002031000
   1175.820000000    0.015436000
    273.188000000    0.073771000
     81.169600000    0.247606000
     27.183600000    0.611832000
      3.413600000    0.241205000
S 1 1.00
      9.532200000    1.000000000
S 1 1.00
      0.939800000    1.000000000
S 1 1.00
      0.284600000    1.000000000
P 4 1.00
     35.183200000    0.019580000
      7.904000000    0.124189000
      2.305100000    0.394727000
      0.717100000    0.627375000
P 1 1.00
      0.213700000    1.000000000
D 1  1.00
       0.85000000        1.00000000
****

使用 nolinear 情况类似
曾经出现过一次这样的情况，反复算了20+次最后才收敛，求解，谢谢
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
NaNa
1楼 2013-07-25 22:09:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

亍寞

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
skysky112211: 金币+3, ★★★★★最佳答案, 谢谢 2013-07-28 14:29:21
opt里面又不是ts，为啥要noraman。。。而且一般来说不需要tight吧，默认就行，窃以为做个频率，没虚频的话就拿来用吧，能量还是几何都差不了多少。。。。。
赞一下 回复此楼
4楼2013-07-26 09:13:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

t13340033021

金虫 (正式写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
skysky112211: 金币+1, ★有帮助, 谢谢 2013-07-28 14:29:59
分子不复杂  看看有没有振荡
赞一下 回复此楼
甲烷无敌
2楼2013-07-25 23:22:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kekexili_08

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
gkf高: 金币+1 2013-07-26 07:47:59
gkf高: 金币+1 2013-07-26 07:48:28
skysky112211: 金币+1, ★★★很有帮助, 谢谢 2013-07-28 14:29:54
就算按照严格的标准， 也可以认为早就收敛了。 差那么点，何必计较呢。
赞一下 回复此楼
3楼2013-07-26 06:13:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skysky112211

木虫 (小有名气)

NaNa
引用回帖:
2楼: Originally posted by t13340033021 at 2013-07-25 23:22:35
分子不复杂  看看有没有振荡

没有振荡，反而几乎是一条直线
赞一下 回复此楼
NaNa
5楼2013-07-26 10:56:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +23 医学老男孩 2026-08-13 49/2450 2026-08-14 20:02 by zhaifei
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] filecode +6 cratir 2026-08-14 10/500 2026-08-14 18:29 by 笑叹辞穷
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +5 工厂打螺丝 2026-08-08 5/250 2026-08-14 11:03 by 玖戈弋
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
信息提示
关闭
请填处理意见
关闭