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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Turbomole/ORCA » 【求助】用Orca计算一个片簇的单点能，不收敛。
51/1返回列表
查看: 4008  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ghcacj

荣誉版主 (著名写手)

[交流] 【求助】用Orca计算一个片簇的单点能，不收敛。

用Orca计算一个片簇的单点能，怎么算都不收敛，之前也算过其他的片簇，一般20次迭代就可以收敛了，这次计算的这个算了125次都不收敛，dE发生振荡。
用的是ubuntu 10.04+orca 2.8
计算输入文件如下：

! ECP{LANLDZ=Zn}
! UKS B3LYP/G 6-31G TightSCF Grid4 SP DFT
! PAL4
* xyz 0 1
C     14.7569      6.7818      3.1858
H     14.913       6.0981      3.7958
N     15.4572      7.9572      3.1056
C     14.7569     10.2092     13.8052
H     14.913      10.8929     13.1952
N     15.4572      9.0338     13.8854
C     13.8052      2.2341      6.7818
H     13.1952      2.0780      6.0981
N     13.8854      1.5338      7.9572
C     13.8052     14.7569     10.2092
H     13.1952     14.913      10.8929
N     13.8854     15.4572      9.0338
C     15.2773     11.6813      6.2614
H     14.5936     12.2913      6.4175
N     16.4527     11.6011      6.9617
C     15.2773      5.3097     10.7296
H     14.5936      4.6997     10.5735
N     16.4527      5.3899     10.0293
C     15.2773      6.2614     11.6813
H     14.5936      6.4175     12.2913
N     16.4527      6.9617     11.6011
H     15.6351      9.1972      9.7579
H     14.9045     10.0961     10.8131
H     15.2256     10.6636      9.3875
C     15.2773     10.7296      5.3097
H     14.5936     10.5735      4.6997
N     16.4527     10.0293      5.3899
H     15.6351      7.7938      7.2331
H     14.9045      6.8949      6.1779
H     15.2256      6.3274      7.6035
C     14.7569      3.1858      6.7818
H     14.913       3.7958      6.0981
H     16.2893      1.2624      9.8514
H     15.3904      2.3176     10.582
H     14.8229      0.8920     10.2609
N     15.4572      3.1056      7.9572
C     14.7569     13.8052     10.2092
H     14.913      13.1952     10.8929
H     16.2893     15.7286      7.1396
H     15.3904     14.6734      6.4090
H     14.8229     16.099      6.7301
N     15.4572     13.8854      9.0338
C     13.8052      6.7818      2.2341
H     13.1952      6.0981      2.0780
H     15.7286      9.8514      0.7017
H     14.6734     10.582      1.6006
H     16.099     10.2609      2.1681
N     13.8854      7.9572      1.5338
C     13.8052     10.2092     14.7569
H     13.1952     10.8929     14.913
H     15.7286      7.1396     16.2893
H     14.6734      6.4090     15.3904
H     16.099       6.7301     14.8229
N     13.8854      9.0338     15.4572
C     14.902       8.6297      2.0890
C     15.3955      9.9482      1.5955
C     14.902       8.3613     14.902
C     15.3955      7.0428     15.3955
C     14.902       2.0890      8.6297
C     15.3955      1.5955      9.9482
C     14.902      14.902       8.3613
C     15.3955     15.3955      7.0428
C     16.8568     10.5845     10.5845
C     15.5383     10.091      10.091
C     16.8568      6.4065      6.4065
C     15.5383      6.9000      6.9000
C     19.2251     10.2092      3.1858
H     19.069      10.8929      3.7958
N     18.5248      9.0338      3.1056
C     19.2251      6.7818     13.8052
H     19.069      6.0981     13.1952
N     18.5248      7.9572     13.8854
C     20.1768     14.7569      6.7818
H     20.7868     14.913      6.0981
N     20.0966     15.4572      7.9572
C     20.1768      2.2341     10.2092
H     20.7868      2.0780     10.8929
N     20.0966      1.5338      9.0338
C     18.7047     11.6813     10.7296
H     19.3884     12.2913     10.5735
N     17.5293     11.6011     10.0293
C     18.7047      5.3097      6.2614
H     19.3884      4.6997      6.4175
N     17.5293      5.3899      6.9617
H     18.3469      7.7938      9.7579
H     19.0775      6.8949     10.8131
H     18.7564      6.3274      9.3875
C     18.7047     10.7296     11.6813
H     19.3884     10.5735     12.2913
N     17.5293     10.0293     11.6011
H     18.3469      9.1972      7.2331
H     19.0775     10.0961      6.1779
H     18.7564     10.6636      7.6035
C     18.7047      6.2614      5.3097
H     19.3884      6.4175      4.6997
N     17.5293      6.9617      5.3899
C     19.2251      3.1858     10.2092
H     19.069       3.7958     10.8929
H     17.6927      1.2624      7.1396
H     18.5916      2.3176      6.4090
H     19.1591      0.8920      6.7301
N     18.5248      3.1056      9.0338
C     19.2251     13.8052      6.7818
H     19.069      13.1952      6.0981
H     17.6927     15.7286      9.8514
H     18.5916     14.6734     10.582
H     19.1591     16.099      10.2609
N     18.5248     13.8854      7.9572
C     20.1768      6.7818     14.7569
H     20.7868      6.0981     14.913
H     18.2534      9.8514     16.2893
H     19.3086      10.582     15.3904
H     17.883      10.2609     14.8229
N     20.0966      7.9572     15.4572
C     20.1768     10.2092      2.2341
H     20.7868     10.8929      2.0780
H     18.2534      7.1396      0.7017
H     19.3086      6.4090      1.6006
H     17.883       6.7301      2.1681
N     20.0966      9.0338      1.5338
C     19.08        8.3613      2.0890
C     18.5865      7.0428      1.5955
C     19.08        8.6297     14.902
C     18.5865      9.9482     15.3955
C     19.08       14.902       8.6297
C     18.5865     15.3955      9.9482
C     19.08        2.0890      8.3613
C     18.5865      1.5955      7.0428
C     17.1252     10.5845      6.4065
C     18.4437     10.091       6.9000
C     17.1252      6.4065     10.5845
C     18.4437      6.9000     10.091
Zn     16.991     12.7432      8.4955
Zn     16.991      4.2477      8.4955
Zn     16.991      8.4955     12.7432
Zn     16.991      8.4955      4.2477
H     13.3084     16.2196      8.7409
H     20.6736     16.2196      8.2501
H     13.3084      8.7409     16.2196
H     20.6736      8.7409      0.7714
H     13.3084      8.2501      0.7714
H     20.6736      0.7714      8.7409
H     20.6736      8.2501     16.2196
H     13.3084      0.7714      8.2501
*
片簇结构如下：


[ Last edited by ghcacj on 2011-3-2 at 16:05 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-03-02 16:01:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)
引用回帖:
Originally posted by recoli at 2011-03-04 00:33:15:
就你的结构来说，我认为电荷是4，闭壳层，所以输入文件写*xyz 4 1。

把C换成甲基是可以的，但是在N上添加H会影响整体的电荷。

不得不说，你的判断非常正确，认真分析了我要计算的结构之后，修改计算文件，真的很快就收敛了，13次循环，非常感谢。
赞一下 回复此楼
8楼2011-03-04 10:13:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

ghcacj

荣誉版主 (著名写手)
我想请教下各位，单点能计算不能收敛，有哪些常用的解决办法？比如Opt？还是调整片簇的切割方法？还是增大基组？
赞一下 回复此楼
2楼2011-03-02 16:39:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

recoli

金虫 (正式写手)
ghcacj(金币+3): 谢谢 2011-03-03 13:59:08
据我所知orca默认是开启RI近似的。但是6-31G这种基组没有辅助基组，所以应该关闭RI近似
%method
RI off
end
赞一下 回复此楼
3楼2011-03-02 17:29:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)
引用回帖:
Originally posted by recoli at 2011-03-02 17:29:36:
据我所知orca默认是开启RI近似的。但是6-31G这种基组没有辅助基组，所以应该关闭RI近似
%method
RI off
end

按照你的设置，基组也改为6-31G*，但是还是不行，有没有更好的建议？
赞一下 回复此楼
4楼2011-03-03 13:59:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 3/150 2026-11-09 21:05 by FChKX8meu1ds
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭