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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 逼近收敛，但是就是不收敛，求助高手
51/1返回列表
查看: 2579  |  回复: 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的回帖

kekexili_08

木虫 (正式写手)
引用回帖:
7楼: Originally posted by skysky112211 at 2013-07-26 10:58:11
感觉总是不好的...

最严格的的判据是所有四项全是0，可是谁又真的全设置成零呢？ 浮点数数值计算，达到一定精度就可以了。
忘掉这个感觉就好了。
赞一下 回复此楼
9楼2013-07-26 15:57:42
已阅   回复此楼   关注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的回帖

亍寞

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +12 医学老男孩 2026-08-13 32/1600 2026-08-14 08:19 by 医学老男孩
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 30/1500 2026-08-14 07:55 by 医学老男孩
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
信息提示
关闭
请填处理意见
关闭