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

24小时热门版块排行榜    

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

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的回帖

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的回帖

skysky112211

木虫 (小有名气)

NaNa
引用回帖:
4楼: Originally posted by 亍寞 at 2013-07-26 09:13:40
opt里面又不是ts，为啥要noraman。。。而且一般来说不需要tight吧，默认就行，窃以为做个频率，没虚频的话就拿来用吧，能量还是几何都差不了多少。。。。。

据说noraman可以节省一些时间，然后要求要tight，还没有加上UltraFine呢
赞一下 回复此楼
NaNa
6楼2013-07-26 10:57:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skysky112211

木虫 (小有名气)

NaNa
引用回帖:
3楼: Originally posted by kekexili_08 at 2013-07-26 06:13:05
就算按照严格的标准， 也可以认为早就收敛了。 差那么点，何必计较呢。

感觉总是不好的
回复此楼
NaNa
7楼2013-07-26 10:58:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

亍寞

金虫 (小有名气)
引用回帖:
6楼: Originally posted by skysky112211 at 2013-07-26 10:57:53
据说noraman可以节省一些时间，然后要求要tight，还没有加上UltraFine呢...

noraman节省的是你有类似opt=calcfc或者freq然后算hessian的时候，不计算光谱那部分。。。。。于是会快。。。嘛。。反正也没出错，无所谓了。。。要还用tight，那就要么方法基组换一个，然后用那个的构型再来这个方法基组opt，或者就在这附近微调吧。。。
赞一下 回复此楼
8楼2013-07-26 11:39:24
已阅   回复此楼   关注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的回帖

t13340033021

金虫 (正式写手)
那可能是势能面平缓
赞一下 回复此楼
甲烷无敌
10楼2013-07-28 14:40:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 skysky112211 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +14 爱看书的可乐 2026-08-26 16/800 2026-08-26 17:44 by skyjudy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 9/450 2026-08-26 15:33 by process2012
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国际合作可查了，中了面上 +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
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +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
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭