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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Turbomole/ORCA » 【求助】用Orca计算一个片簇的单点能，不收敛。
51/1返回列表
查看: 3462  |  回复: 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(金币+5): 谢谢 2011-03-03 20:14:38
我看了一下你的结构，Zn应该都是Zn2+离子吧？中间作为桥的四个配体每个少一个H是-1价，那么总体电荷应该是+4？楼主认为呢？

电荷和多重度是很重要的，一定要设置正确，不然即使收敛了得到的也不是正确结果。
赞一下 回复此楼
5楼2011-03-03 17:50:42
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿085502，267分求调剂 +13 再忙也要吃饭啊 2026-04-08 14/700 2026-04-12 18:40 by laoshidan
[考研] 2本，初试303，0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 一志愿085802 323分求调剂 +9 drizzle_9 2026-04-12 10/500 2026-04-12 14:51 by drizzle_9
[考研] 一志愿华中农微生物，288分，三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[找工作] 山东高校教师考核超级无底线，员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 一志愿郑州大学 22408 305分求调剂 +5 安小满zzz 2026-04-08 5/250 2026-04-12 00:41 by 蓝云思雨
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 070300化学279求调剂 +19 哈哈哈^_^ 2026-04-08 20/1000 2026-04-11 20:43 by stoner78
[考研] 296求调剂 +14 汪！？！ 2026-04-08 15/750 2026-04-11 20:28 by dongdian1
[考研] 300分求调剂 （085501机械专硕，本科扬大） +8 xu@841019 2026-04-11 8/400 2026-04-11 10:46 by qingpingzhu
[考研] 材料复试求调剂 +20 xhhdjdjsjks 2026-04-09 20/1000 2026-04-10 10:25 by 孙小小12457
[考研] 297求调剂 +27 GENJIOW 2026-04-07 30/1500 2026-04-09 23:20 by wolf97
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
[考研] 一志愿华东理工085601材料工程303分求调剂 +15 a1708 2026-04-06 15/750 2026-04-08 16:23 by luoyongfeng
[考研] 313求调剂 +3 十六拾陆 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
[考研] 316求调剂 +4 15318418673 2026-04-07 4/200 2026-04-07 22:12 by hemengdong
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
信息提示
关闭
请填处理意见
关闭