| 查看: 3326 | 回复: 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 ] |
回复此楼» 猜你喜欢
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
求助火焰封管的时候管子炸了
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有131人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有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/172
-
北京理工大学 珠海校区全职院士招数名博士生--申请考核制-半导体、光学、微电子
+2/118
-
征婚
+1/87
-
深圳大学信息功能聚合物电介质方向“申请-考核制”博士生招生
+1/78
-
大连理工大学智能系统实验室优秀硕博研究生招生
+1/76
-
[长期合作招募] 同济大学肖倩老师团队诚邀港澳学者学术交流
+1/53
-
[长期合作招募] 同济大学肖倩老师团队诚邀港澳学者学术交流
+1/34
-
中国农业科学院招聘化学，化工，材料相关专业科研助理 （北京）
+2/34
-
上海交通大学 碳纳米管、阵列、器件应用博士后招聘
+1/32
-
北欧岗位制博士招生（锂离子电池正极材料和粘结剂）
+1/30
-
中国科学院深圳先进技术研究院合成生物学研究所陈禹课题组招聘博士后
+2/26
-
学生申博
+1/21
-
2026申博自荐 本硕双一流学科 纳米药物递送方向 一篇一区TOP 两个国家奖学金
+1/17
-
华南理工大学宋波教授招收2026年博士生（二氧化碳有机、高分子和光电生物催化还原）
+1/8
-
华南理工大学宋波教授招收2026年博士生（二氧化碳转化方向优先）
+1/6
-
复旦大学2026年博士研究生，“药学 + 合成生物学”
+1/5
-
山东师范大学海外优青（校长团队）招聘硕士/博士/博后
+1/2
-
河北工业大学层状材料与器件团队诚聘二维材料与原位电子显微学方向青年教师与博士后
+1/2
-
换工作
+1/2
-
深圳理工大学招聘AI4Science与分子模拟方向科研人员
+1/1
|
麻烦帮我看看呢，量子化学方面我是入门选手。 我的结构是一个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
2楼2011-03-02 16:39:11
3楼2011-03-02 17:29:36
4楼2011-03-03 13:59:41