| 查看: 1148 | 回复: 10 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】金属配合物计算输出文件问题
|
|||
|
前几天老板要求算下金属配合物的弱作用，于是初步选定了具有高自旋的四面体构型四氯合镍，其中2价Ｎi的3d轨道有２个孤电子. 计算过程中出现：1）初始占据和空轨道不能识别，和2）最大距离和均方根距离收敛致零后，构型和能量不变，直到循环结束“Error termination request processed by link 9999”。 表达可能不够清楚，附件中贴出输入文件，输出文件中疑问如下： 1） Initial guess orbital symmetries: Alpha Orbitals: Occupied (A1) (T2) (T2) (T2) (T2) (T2) (T2) (A1) (?A) (?A) (?A) (T2) (?A) (?A) (?A) (T2) (?A) (?A) (?A) (T2) (?B) (T2) (T2) (T2) (A1) (?C) (?C) (?C) (?C) (?C) (?B) (?C) (?B) (?B) (?C) (?C) (?C) (?B) (?B) (?C) (?C) (?C) (?C) (?C) (?C) Virtual (T2) (T2) (T2) (A1) (T2) (T2) (T2) (T2) (T2) (T2) (T1) (T1) (T1) (?D) (?C) (?C) (?C) (?C) (?C) (?C) (?C) (?C) (?C) (?D) (?D) (?E) (?E) (?E) (?D) (?E) (?E) (?E) (?D) (?D) (?E) (?E) (?E) (?E) (?E) (?E) (?D) (?E) (?E) (?E) (?F) (?F) (?F) (?F) (?F) (E) (?F) (E) (?F) (?F) (?F) (?F) (?F) (?F) (T2) (T2) (T2) (?F) (?F) (?F) (A1) Beta Orbitals: Occupied (A1) (T2) (T2) (T2) (T2) (T2) (T2) (A1) (?A) (?A) (?A) (T2) (?A) (?A) (?A) (T2) (?A) (?A) (?A) (T2) (A1) (T2) (T2) (T2) (A1) (T2) (T2) (T2) (?B) (?C) (?B) (?B) (?C) (?C) (?C) (?C) (?C) (?B) (?C) (?C) (?C) (?B) (?C) Virtual (?C) (?C) (?C) (?C) (?C) (A1) (?C) (?C) (?C) (E) (E) (A1) (T1) (T1) (T1) (A1) (?C) (?C) (?C) (?C) (?C) (?C) (?C) (?C) (?C) (?B) (?B) (?C) (?C) (?C) (A1) (?C) (?C) (?C) (?B) (?B) (?C) (?C) (?C) (?C) (?C) (?C) (?B) (?C) (?C) (?C) (?D) (?D) (?D) (?D) (?D) (?D) (E) (E) (?D) (?D) (?D) (T2) (T2) (T2) (?B) (?B) (?B) (?C) (?C) (?C) (A1) 2） Item Value Threshold Converged? Maximum Force 0.010052 0.000450 NO RMS Force 0.005509 0.000300 NO Maximum Displacement 0.000000 0.001800 YES RMS Displacement 0.000000 0.001200 YES 上述问题原因何在，如何解决，请虫友不吝赐教，谢谢！ |
回复此楼» 猜你喜欢
中科院西安光机所——聚焦功能材料创新｜诚邀参与 2026 年度学术研讨会。
已经有1人回复
-
中科院西安光机所——促进学术交流，推动产学研融合，分享领域成果，加强校企合作。
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有52人回复
-
中科院西安光机所——促进学术交流，推动产学研融合，分享领域成果，加强校企合作。
已经有0人回复
-
分子动力学 Plumed教程求助
已经有0人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
申请2027年材料/化学类博士
已经有5人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有24人回复
|
输入文件的附件没有贴上，直接写下面了 %chk=ni.chk %mem=600MB # opt=maxcycle=200 freq b971/lanl2dz extrabasis test NiCl4 -2 3 Ni 0.00000000 0.00000000 0.00000000 Cl 0.00000000 0.00000000 2.14000000 Cl 2.01761157 0.00000000 -0.71333270 Cl -1.00880624 -1.74730262 -0.71333270 Cl -1.00880538 1.74730270 -0.71333371 Cl 0 6-31+G(d,p) **** |
2楼2009-10-31 16:17:51
wuy069
铁杆木虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 7921.9
- 散金: 15
- 红花: 10
- 帖子: 506
- 在线: 549.1小时
- 虫号: 760723
- 注册: 2009-04-30
- 性别: GG
- 专业: 理论和计算化学
3楼2009-10-31 16:28:21
erylingjet
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 5.372
- 金币: 9933.4
- 散金: 752
- 红花: 2
- 帖子: 2248
- 在线: 241.5小时
- 虫号: 415302
- 注册: 2007-06-28
- 性别: GG
- 专业: 计算化学
4楼2009-10-31 19:19:46
5楼2009-10-31 20:22:02
liujun2006
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 604.6
- 帖子: 289
- 在线: 84.8小时
- 虫号: 353658
- 注册: 2007-04-23
- 性别: GG
- 专业: physics
6楼2009-10-31 20:24:41
|
５楼说的是能量震荡，导致超过循环次数而不收敛，加入步长后再寻找的策略． 而我这个估计是进入了个死循环，因为计算很早就达到了 Item Value Threshold Converged? Maximum Force 0.010052 0.000450 NO RMS Force 0.005509 0.000300 NO Maximum Displacement 0.000000 0.001800 YES RMS Displacement 0.000000 0.001200 YES 而余下所有的循环所出现的都是一个结果（如下），并没有震荡， SCF Done: E(UB97-1+HF-B97-1) = -2010.00018366 A.U. after 1 cycles 我想是不是由于分子说明部分或路径哪块的说明不够确切，导致ＧＳ不识别初始轨道，才导致这种情况 具体怎么改，还是不知道－－头疼！ |
7楼2009-10-31 21:48:01
8楼2009-11-01 00:36:59
erylingjet
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 5.372
- 金币: 9933.4
- 散金: 752
- 红花: 2
- 帖子: 2248
- 在线: 241.5小时
- 虫号: 415302
- 注册: 2007-06-28
- 性别: GG
- 专业: 计算化学
9楼2009-11-01 08:45:45
10楼2009-11-01 08:46:26