| 查看: 1058 | 回复: 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 上述问题原因何在，如何解决，请虫友不吝赐教，谢谢！ |
回复此楼» 猜你喜欢
26年储能、电池方向博士申请
已经有2人回复
-
大连工业大学杰青/长江团队-生物质材料-储能电池方向招收2026级博士生
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有172人回复
-
欢迎硕博咨询
已经有26人回复
-
求审搞机会
已经有4人回复
-
昆明理工大学冶能院离子液体冶金课题组招收博士生
已经有2人回复
-
湖南大学刘巧玲课题组2026年博士研究生招生信息
已经有28人回复
-
锂电池正极材料前驱体法制备正极材料
已经有2人回复
-
深圳理工大学-湖南大学招收电化学方向联培博士生1名
已经有2人回复
-
南方科技大学表面催化反应动力学研究课题组招收2026届博士生
已经有0人回复
|
输入文件的附件没有贴上，直接写下面了 %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