应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】金属配合物计算输出文件问题
111/2返回列表12下一页
查看: 1229  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

xuhu_11

金虫 (职业作家)

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

前几天老板要求算下金属配合物的弱作用，于是初步选定了具有高自旋的四面体构型四氯合镍，其中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

上述问题原因何在，如何解决，请虫友不吝赐教，谢谢！
回复此楼

» 猜你喜欢
千万别考博！
1楼 2009-10-31 16:12:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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中有个例子就是你这个化合物吧，或者四羰基合镍
用分片的思想
赞一下(11人) 回复此楼
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的回帖

xuhu_11

金虫 (职业作家)
ADF -!?
应该和所有的软件关系不大吧
谢谢2楼,我查查看


to 3 楼
用了额外基组关键词 extrabasis,理解为额外对Cl原子l用6-31+(d,p)基组
赞一下 回复此楼
千万别考博！
5楼2009-10-31 20:22:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liujun2006

金虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+1,VIP+0):感谢交流 11-1 08:46
L999步出错了，可能是优化进入局部最小而跳不出来。试一下用IOP控制，一般写IOP(1/8=?)， ？为一个数字，具体作用为减小步长。你可以查看一下IOP(1/8)的解释。
赞一下(11人) 回复此楼
6楼2009-10-31 20:24:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuhu_11

金虫 (职业作家)
５楼说的是能量震荡，导致超过循环次数而不收敛，加入步长后再寻找的策略．
而我这个估计是进入了个死循环，因为计算很早就达到了
     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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gonemrc

捐助贵宾 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
偶菜鸟，不过前天在版上问过这个问题。在后面读入赝势有没有用的？
赞一下(1人) 回复此楼
8楼2009-11-01 00:36:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)
引用回帖:
Originally posted by xuhu_11 at 2009-10-31 20:22:


to 3 楼
用了额外基组关键词 extrabasis,理解为额外对Cl原子l用6-31+(d,p)基组

呵呵,不好意思,没看到.
赞一下 回复此楼
9楼2009-11-01 08:45:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuhu_11

金虫 (职业作家)
不用写　lanl2dz会自动读入基组信息和赝势
当然你觉得不妥，也可以用gencep重新定义基组和赝势
赞一下 回复此楼
千万别考博！
10楼2009-11-01 08:46:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xuhu_11 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭