24小时热门版块排行榜    

CyRhmU.jpeg
查看: 921  |  回复: 10
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xuhu_11

金虫 (职业作家)

[交流] 【求助】金属配合物计算输出文件问题

前几天老板要求算下金属配合物的弱作用,于是初步选定了具有高自旋的四面体构型四氯合镍,其中2价Ni的3d轨道有2个孤电子.
     计算过程中出现: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

上述问题原因何在,如何解决,请虫友不吝赐教,谢谢!
回复此楼
千万别考博!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuhu_11

金虫 (职业作家)

解决了,加关键词nosymm scf=tight
看来不限制对称性经常容易让GS短路
千万别考博!
11楼2009-11-01 14:20:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

xuhu_11

金虫 (职业作家)

输入文件的附件没有贴上,直接写下面了
%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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuy069

铁杆木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
erylingjet(金币+1,VIP+0):感谢交流 10-31 19:19
用ADF计算吧
ADF中有个例子就是你这个化合物吧,或者四羰基合镍
用分片的思想
3楼2009-10-31 16:28:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)

b971/lanl2dz

Cl 0
6-31+G(d,p)
****

这两种方法不冲突吗?
4楼2009-10-31 19:19:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见