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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Turbomole/ORCA » 【求助】用Orca计算一个片簇的单点能，不收敛。
51/1返回列表
查看: 3407  |  回复: 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的回帖

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的回帖
查看全部 8 个回答

ghcacj

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

recoli

金虫 (正式写手)
ghcacj(金币+5): 谢谢 2011-03-03 20:14:38
我看了一下你的结构，Zn应该都是Zn2+离子吧？中间作为桥的四个配体每个少一个H是-1价，那么总体电荷应该是+4？楼主认为呢？

电荷和多重度是很重要的，一定要设置正确，不然即使收敛了得到的也不是正确结果。
赞一下 回复此楼
5楼2011-03-03 17:50:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工科 320总分 求调剂 +11 梨花珞晚风 2026-03-17 11/550 2026-03-22 17:42 by luoyongfeng
[考研] 一志愿中南化学（0703）总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考研] 384求调剂 +3 子系博 2026-03-22 4/200 2026-03-22 11:04 by 搏击518
[考研] 0856材料专硕353求调剂 +4 NIFFFfff 2026-03-20 4/200 2026-03-22 09:49 by 2026paper
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕，总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 085410人工智能专硕317求调剂（0854都可以） +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 一志愿福大288有机化学，求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
信息提示
关闭
请填处理意见
关闭