24小时热门版块排行榜    

查看: 1250  |  回复: 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的回帖

liujun2006

金虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
erylingjet(金币+1,VIP+0):感谢交流 11-1 08:46
L999步出错了,可能是优化进入局部最小而跳不出来。试一下用IOP控制,一般写IOP(1/8=?), ?为一个数字,具体作用为减小步长。你可以查看一下IOP(1/8)的解释。
6楼2009-10-31 20:24:41
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:39 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:03 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 00:51 by 9i7PCSZyMzot
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:15 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:03 by 9i7PCSZyMzot
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 4pRyVGCVerAA 2026-09-19 4/200 2026-09-20 00:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:51 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:39 by 9i7PCSZyMzot
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-18 7/350 2026-09-19 21:03 by 9i7PCSZyMzot
[教师之家] 现代”学阀”该如何界定 +3 iaeyuan 2026-09-19 3/150 2026-09-19 20:52 by 没昵称呀
[基金申请] 国社科系统bug了,是不是要放榜了? +7 kynobel 2026-09-16 9/450 2026-09-19 14:52 by 接好运婷
[考博] 申博发邮件 +8 Lmengk 2026-09-15 9/450 2026-09-19 09:27 by 中间体超市
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
信息提示
请填处理意见