24小时热门版块排行榜    

查看: 3573  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

piaoxue001

木虫 (正式写手)

[交流] 【求助】高斯计算过渡态不收敛,请问如何解决? 已有6人参与

高斯计算过渡态设定maxcyc=800,出现link9999错误
我查看了一下能量收敛的情况,发现能量在optimization step number=200后趋于平缓,但是能量依然有下降趋势,到300左右报link9999错误结束,我不明白的是我设置的800,为什么到300就结束了,尽管我查阅资料,有人说200基本上就可以了。但看我的趋势,应该在300以上。另外能量未出现振荡,构型未偏离预计很大,前两个收敛标准已经达到。
请教高手指点一下!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

piaoxue001

木虫 (正式写手)

非常感谢

引用回帖:
Originally posted by jiewei at 2010-08-18 14:02:26:
很简单, 你看看输出坐标,是不是有-179.9以下的角或者二面角? 如果有,那把输入结构调整一下; 还有一个可能是你这个构型是亚稳态(你预计构型它就是一个亚稳态), 然后你的优化步长太小,它跳不出去,所以报错, 办法就是 ...

我查看了一下结构坐标
C                  1              B1
C                  2              B2    1              A1
C                  3              B3    2              A2    1              D1
C                  4              B4    3              A3    2              D2
C                  5              B5    4              A4    3              D3
H                  1              B6    2              A5    3              D4
H                  2              B7    1              A6    6              D5
H                  4              B8    3              A7    2              D6
H                  5              B9    4              A8    3              D7
H                  6             B10    5              A9    4              D8
C                  3             B11    2             A10    1              D9
H                 12             B12    3             A11    2             D10
H                 12             B13    3             A12    2             D11
Pd                12             B14    3             A13    2             D12
P                 15             B15   12             A14    3             D13
H                 16             B16   15             A15   12             D14
H                 16             B17   15             A16   12             D15
H                 16             B18   15             A17   12             D16
N                 15             B19   12             A18    3             D17
C                 20             B20   15             A19   12             D18
C                 20             B21   15             A20   12             D19
C                 21             B22   20             A21   15             D20
O                 21             B23   20             A22   15             D21
C                 22             B24   20             A23   15             D22
O                 22             B25   20             A24   15             D23
H                 23             B26   21             A25   20             D24
H                 23             B27   21             A26   20             D25
H                 25             B28   22             A27   20             D26
H                 25             B29   22             A28   20             D27
Br                15             B30   12             A29    3             D28
H                  4             B31    3             A30    2             D29
C                  4             B32    3             A31    2             D30
C                 33             B33    4             A32    3             D31
H                 33             B34    4             A33    3             D32
H                 34             B35   33             A34    4             D33
C                 34             B36   33             A35    4             D34
H                 37             B37   34             A36   33             D35
H                 37             B38   34             A37   33             D36
Sn                37             B39   34             A38   33             D37
C                 40             B40   37             A39   34             D38
H                 41             B41   40             A40   37             D39
H                 41             B42   40             A41   37             D40
H                 41             B43   40             A42   37             D41
C                 40             B44   37             A43   34             D42
H                 45             B45   40             A44   37             D43
H                 45             B46   40             A45   37             D44
H                 45             B47   40             A46   37             D45
C                 40             B48   37             A47   34             D46
H                 49             B49   40             A48   37             D47
H                 49             B50   40             A49   37             D48
H                 49             B51   40             A50   37             D49

   B1             1.39122455
   B2             1.41350259
   B3             1.41335341
   B4             1.39219717
   B5             1.39894083
   B6             1.08637467
   B7             1.08731246
   B8             1.08786146
   B9             1.08640513
   B10            1.08650274
   B11            1.45879801
   B12            1.09359886
   B13            1.08613443
   B14            2.13850356
   B15            2.37061717
   B16            1.41991627
   B17            1.41865635
   B18            1.42363184
   B19            2.16544705
   B20            1.38307158
   B21            1.34996533
   B22            1.54732165
   B23            1.21635040
   B24            1.52282849
   B25            1.24454931
   B26            1.09449281
   B27            1.09458810
   B28            1.09485850
   B29            1.09521374
   B30            2.58593546
   B31            4.02981891
   B32            3.84770751
   B33            1.34143270
   B34            1.08876746
   B35            1.09254371
   B36            1.49821885
   B37            1.09999507
   B38            1.10069955
   B39            2.10784808
   B40            2.08871486
   B41            1.09929576
   B42            1.09926763
   B43            1.09947263
   B44            2.08376834
   B45            1.09908327
   B46            1.09939817
   B47            1.09920896
   B48            2.08010944
   B49            1.09932878
   B50            1.09907994
   B51            1.09934648
   A1           120.61395031
   A2           118.68469811
   A3           120.39649223
   A4           120.23659543
   A5           119.88643000
   A6           120.24431737
   A7           119.19252022
   A8           119.73807165
   A9           120.02378131
   A10          120.62682148
   A11          115.54078604
   A12          115.39438141
   A13          114.83577002
   A14           99.33750922
   A15          119.47275684
   A16          114.40111541
   A17          118.19436539
   A18           85.69175075
   A19          146.90152732
   A20           99.23651705
   A21          107.65583440
   A22          126.31412802
   A23          111.37519101
   A24          120.21916599
   A25          108.71873850
   A26          108.95307532
   A27          109.56096205
   A28          109.92894906
   A29          144.04807010
   A30          133.32337737
   A31          118.35032597
   A32           92.06457721
   A33           87.15548802
   A34          117.58023049
   A35          126.80950799
   A36          109.78775566
   A37          109.56722851
   A38          114.89610134
   A39          108.18912521
   A40          114.17933236
   A41          114.07988253
   A42          114.32752839
   A43          109.43121527
   A44          113.62686422
   A45          114.36534258
   A46          114.25697485
   A47          109.57470314
   A48          114.59512459
   A49          114.04994755
   A50          113.56515823
   D1            -0.28650311
   D2             0.12204403
   D3             0.02233741
   D4          -179.50023444
   D5          -179.76580633
   D6           179.09296519
   D7           179.46132663
   D8           179.70311340
   D9           178.46793401
   D10         -156.14200549
   D11          -20.14922024
   D12          100.13039483
   D13          -17.19505409
   D14           87.53883784
   D15         -153.57153588
   D16          -34.97533972
   D17         -123.73390869
   D18          100.96846522
   D19          -90.93846392
   D20          166.88759735
   D21          -13.45206325
   D22         -172.76545623
   D23            7.12895548
   D24         -121.29333609
   D25          122.29630860
   D26          121.20101761
   D27         -121.44413902
   D28           76.10572769
   D29          177.20762037
   D30         -177.16710682
   D31          -43.27840833
   D32           78.44492964
   D33          -93.16072218
   D34           90.48342815
   D35         -144.57043355
   D36          -29.73292330
   D37           89.76983800
   D38           57.58776868
   D39         -178.97560458
   D40          -59.20034253
   D41           61.02218138
   D42          176.74680962
   D43           53.95111138
   D44          174.18871411
   D45          -65.37768679
   D46          -61.72298361
   D47          175.80255414
   D48          -63.64145914
   D49           55.62488734
请问我如何调整?
另外,我的优化步长为maxstep=5。请问加大循环步数有可能收敛吗?
3楼2010-08-19 08:43:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

jiewei

专家顾问 (著名写手)

小才


小木虫(金币+0.5):给个红包,谢谢回帖交流
很简单, 你看看输出坐标,是不是有-179.9以下的角或者二面角? 如果有,那把输入结构调整一下; 还有一个可能是你这个构型是亚稳态(你预计构型它就是一个亚稳态), 然后你的优化步长太小,它跳不出去,所以报错, 办法就是你还是调整输入构型,并且加大优化步长去优化
我想做个科学家,可我又喜欢写诗;我想做一个诗人,可我却读的理工科..
2楼2010-08-18 14:02:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiewei

专家顾问 (著名写手)

小才


小木虫(金币+0.5):给个红包,谢谢回帖交流
输出坐标,指的是你出错前一步的坐标,你这个是? 看起来不像哈
你比如说4          -179.50023444
   D5          -179.76580633
   D6           179.09296519
   D7           179.46132663
   D8           179.70311340
这几个角,经过优化后,可能就出现-180.01或者180.01 这个角度是高斯不认的,它就报错.其实-180.01不就是0.1么?
我想做个科学家,可我又喜欢写诗;我想做一个诗人,可我却读的理工科..
4楼2010-08-19 10:37:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

piaoxue001

木虫 (正式写手)

引用回帖:
Originally posted by jiewei at 2010-08-19 10:37:41:
输出坐标,指的是你出错前一步的坐标,你这个是? 看起来不像哈
你比如说4          -179.50023444
   D5          -179.76580633
   D6           179.09296519
   D7           179.46132663
   D8         ...

你好 这就是我link 9999错误最后的坐标。没有你指出的那种情况
5楼2010-08-23 20:16:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见