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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 逼近收敛，但是就是不收敛，求助高手
北京石油化工学院2026年研究生招生接收调剂公告
141/2返回列表12下一页
查看: 2316  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 332求调剂 +14 小小孟... 2026-04-05 15/750 2026-04-05 22:27 by hemengdong
[考研] 348求调剂 +3 车厘子zzz 2026-04-05 3/150 2026-04-05 20:30 by 啵啵啵0119
[考研] 308求调剂 +4 maverick^_^ 2026-04-03 4/200 2026-04-05 19:08 by 蓝云思雨
[考研] 调剂 +3 好好读书。 2026-04-02 3/150 2026-04-05 13:02 by arrow8852
[考研] 298分 070300求调剂 +15 zwen03 2026-04-02 15/750 2026-04-05 12:52 by Hdyxbekcb
[考研] 调剂求助 +10 想换手机不想解� 2026-04-02 13/650 2026-04-05 09:41 by sam3303
[考研] 材料调剂 +18 一样YWY 2026-04-02 19/950 2026-04-04 22:14 by hemengdong
[考研] 求调剂 +4 晟功? 2026-04-03 4/200 2026-04-04 21:58 by hemengdong
[考研] 315求调剂 +13 小羊小羊_ 2026-04-02 14/700 2026-04-04 20:30 by 蓝云思雨
[考研] 调剂 +4 是可乐不是可乐 2026-04-04 4/200 2026-04-04 19:41 by 唐沐儿
[考研] 一志愿武理材料工程302调剂环化或化工 +19 Doleres 2026-03-31 20/1000 2026-04-04 16:44 by 啊俊！
[考研] 319求调剂 +4 星星不眨眼喽 2026-04-03 4/200 2026-04-04 16:25 by 中飞院空管学院�
[考研] 283求调剂 +3 jiouuu 2026-04-03 4/200 2026-04-03 13:28 by jiouuu
[考研] 313求调剂 +3 ～微微凉～ 2026-04-03 3/150 2026-04-03 11:25 by 啵啵啵0119
[考研] 330求调剂 +3 白神呜呼呼 2026-04-02 3/150 2026-04-03 10:15 by 蓝云思雨
[考研] 326求调剂 +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[考研] 一志愿安徽大学计算机科学与技术学硕，331分求调剂 +5 蒋昌鹏qtj 2026-04-01 5/250 2026-04-02 08:10 by fxue1114
[考研] 求调剂 +4 图鉴212 2026-03-30 5/250 2026-04-01 15:32 by 图鉴212
[考研] 326求调剂 +4 崽崽仔 2026-03-31 4/200 2026-04-01 09:58 by 我的船我的海
[考研] 085601英二数二求调剂 总分325 +4 余航航 2026-03-31 4/200 2026-03-31 17:38 by 唐沐儿
信息提示
关闭
请填处理意见
关闭