| 查看: 3403 | 回复: 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 ] |
回复此楼» 猜你喜欢
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有142人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
分子动力学模拟合作需求
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何处理自旋计算产生的能带
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
祈福中标
+1/233
-
电力全国重点实验室双一流A类长江学者团队招2026年全日制博士1名（4月3日截止）
+2/96
-
双一流南京医科大学招计算机、AI、统计、生物信息等方向26年9月入学博士
+1/80
-
“双一流”高校化工类博士“申请-考核”招生 江苏南京
+5/75
-
江苏科技大学能源材料化学课题组郭兴梅教授招收博士研究生1-2名
+1/44
-
湘潭大学“过程强化与绿色化工”创新团队补招2026年秋入学博士生
+1/35
-
上海交通大学沈道智副教授招收2026年联培博士/硕士研究生（微纳器件与制造方向）
+1/33
-
福建师范大学闽江学者陈巍教授团队招收2026年化学、化工或材料专业的学硕6名
+1/31
-
巴黎西岱大学（Université Paris Cité）Fayna MAMMERI 导师寻找CSC公派博士生
+1/22
-
中山大学院士团队微纳器件、脑机接口方向博士招生、博后招聘
+1/18
-
西南某211招27届博士
+1/12
-
南昌一本市政工程接受调剂
+1/12
-
哈尔滨工业大学(深圳)国家级青年人才机械工程、热能与能源等专业26年入学博士招生
+1/11
-
江西赣州-赣南师范大学-学校推荐
+1/7
-
福州大学国家级人才汤育欣教授课题组2026年招收”申请考核制“储能电池方向博士研究生
+1/6
-
【第三轮招生】澳科大诚招2026年秋季药剂学硕士研究生（3月31日下午18:00截止）
+1/5
-
上海交大招收26年级秋季入学博士生，CET-6过，生物医用材料方向（有高分子合成基础）
+1/5
-
总分265计算机求跨调环境
+1/3
-
2026考研的同学，分数不够？想调剂上岸的同学，来看
+1/2
-
西安交通大学生命学院智能化诊疗装备研究中心2026年博士生招生
+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
