| 查看: 3315 | 回复: 7 | |||
[交流]
【求助】用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人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有70人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
求助立方晶体In(OH)3的CIF卡片
已经有2人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有1人回复
-
阴离子交换膜电解二氧化碳还原
已经有0人回复
-
分子动力学模拟合作需求
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何处理自旋计算产生的能带
已经有3人回复
- 学位学历证书和成绩单的公证都是2010年12月20日的，现在签证能用吗？
已经有6人回复
- 请教：orca 中用CCSD(T)计算单点能出错--提示内存不足 该怎么办？
已经有10人回复
- 大家说说管理集群时的奇闻怪事吧
已经有22人回复
- 难道每次计算前都咬进行收敛性测试吗？
已经有14人回复
- 收敛精度
已经有7人回复
- 计算时候总是不能收敛如何解决？
已经有5人回复
- MOLCAS中采用CASPT2计算CO2二价阳离子1a2态不收敛怎么办
已经有4人回复
- vasp计算：优化能收敛，但同样的设置静态计算不收敛
已经有3人回复
- 【求助】Orca 2.8的一个计算赝势的例子报错
已经有12人回复
- 【ghcacj个人文集】单机环境下安装openmpi使用Orca并行运算的方法介绍
已经有12人回复
- 【求助】请教使用orca计算出的这个自旋密度图是怎么画的？
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
辽宁材料实验室高热通量框架复合材料方向招聘科研骨干2名（长期有效）
+5/230
-
几个高校工作如何选择?
+1/86
-
上海交通大学变革分子学中心申涛课题组2026秋季入学推荐-考核制博士招生（有机）
+1/81
-
招聘：中国科学院山西煤炭化学研究所
+1/74
-
[长期合作招募] 同济大学肖倩老师团队诚邀港澳学者学术交流
+1/57
-
香港理工大学招收电力系统优化及运筹学博士后
+1/55
-
坐标广州，诚征男友，大个子女生，非诚勿扰
+2/44
-
211双一流北京工业大学招计算机、AI、自动控制、电子信息等方向博士生（长期有效）
+1/40
-
[长期合作招募] 同济大学肖倩老师团队诚邀港澳学者学术交流
+1/34
-
西工大可塑性液体材料课题组招收软物质、流体力学、细胞工程等背景博士生
+1/33
-
【教授本人】南佛罗里达大学化学系刘文奇课题组 2026 Fall 招收有机/超分子方向博士生
+1/31
-
[长期合作招募] 同济大学肖倩老师团队诚邀港澳学者学术交流
+1/30
-
悉尼大学 AMME 机械工程 双ARC Future Fellows团队招收CSC博士生
+1/27
-
26博士申请-药物化学方向
+2/22
-
浙江大学信息光子材料与器件实验室诚聘博士后、科研助理
+1/15
-
悉尼麦考瑞大学Shujuan Huang 研究团队 博士生招生 CSC 奖学金 补齐全奖RMB55.5万/3年
+2/8
-
山东师范大学海外优青（校长团队）招聘硕士/博士/博后
+1/4
-
Plant communications难度如何？
+1/4
-
生殖医学与子代健康全国重点实验室华鹏课题组招收博士后及研究生(长期有效)
+1/2
-
【博士后/科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+1/2
2楼2011-03-02 16:39:11
3楼2011-03-02 17:29:36
4楼2011-03-03 13:59:41
5楼2011-03-03 17:50:42
|
麻烦帮我看看呢，量子化学方面我是入门选手。 我的结构是一个Zn作为节点，连接四个甲基咪唑，形成一种类似于沸石中Si-O西面体的结果，为了计算电荷，按照文献的方法，截取出这样的一个片段进行量化计算，因为是截取的片段，因此终端需要挂载一定的原子来补充，之前碰到的结果都是末端原本连接一个C原子，所以填上一个甲基，现在末端应该连接Zn，末端采用H作为补充不知道是否合适，但是看到文献里有这么做，没想到怎么计算都不收敛，按理说总体电荷应该为0，多重度为1，可是就是不行，正在尝试是否是开壳层导致的，有的人建议将同样的结构先用闭壳层计算，然后拿这个结果作为intial guess进行开壳层计算，但是结果还没出来。祈祷了。 谢谢你的关注。 [ Last edited by ghcacj on 2011-3-3 at 20:21 ] |
6楼2011-03-03 20:14:25
7楼2011-03-04 00:33:15
8楼2011-03-04 10:13:17