| 查看: 3324 | 回复: 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 ] |
回复此楼» 猜你喜欢
香港科技大学（广州）诚招电催化方向博士生（2026秋入学）
已经有0人回复
-
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有289人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
求助Cu2+1O的CIF文件(PDF: 05-0667)
已经有1人回复
-
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
-
KAUST（阿卜杜拉国王科技大学）MXene 器件方向博士后招聘
已经有0人回复
-
沙特阿拉伯阿卜杜拉国王科技大学（KAUST）电池方向博士后招聘
已经有0人回复
-
分子动力学模拟合作需求
已经有0人回复
-
福州大学新能源材料与工程研究院招收2026年入学博士
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
浙江师范大学招收大气遥感、光学或人工智能方向的2026年物理学博士生1名
+2/274
-
专任教师招聘
+1/172
-
北京理工大学 珠海校区全职院士招数名博士生--申请考核制-半导体、光学、微电子
+2/130
-
保障电池制造超纯环境：威格科技电池生产全线解决方案
+1/86
-
几个高校工作如何选择?
+1/76
-
大连理工大学智能系统实验室优秀硕博研究生招生
+1/76
-
中国科协青年托举（针对32岁青年人才的）是不是不搞了，25年都没有通知
+1/68
-
辽宁材料实验室框架复合材料课题组招收联合培养研究生（长期有效）
+2/68
-
华中科技大学周英教授招博后
+1/30
-
美国密苏里大学堪萨斯城分校（UMKC）生物材料诚聘全奖博士
+1/29
-
[长期合作招募] 同济大学肖倩老师团队诚邀港澳学者学术交流
+1/25
-
专科起点双非硕士，申博求建议
+1/20
-
中国地质大学（武汉）分析地球化学团队招收博士生1名、硕士生3名
+2/12
-
中山大学农业与生物技术学院周潇峰课题组诚聘微生物/植物病理学方向科研助理
+1/11
-
复旦大学2026年博士研究生，“药学 + 合成生物学”
+1/5
-
国防科大超声速燃烧课题组招收博士后
+1/5
-
0854电子信息330求调剂
+1/4
-
南开大学齐迹课题组诚聘分子生物学、免疫学、有机分子合成相关方向的博士后、博士硕士
+1/3
-
河北工业大学层状材料与器件团队诚聘二维材料与原位电子显微学方向青年教师与博士后
+1/2
-
多伦多城市大学计算机视觉方向博士后
+1/1
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