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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 用Gaussian03或09，MP2(full)/aug-cc-pvdz如何优化HI？
51/1返回列表
查看: 1689  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

1maomao

金虫 (初入文坛)

[求助] 用Gaussian03或09，MP2(full)/aug-cc-pvdz如何优化HI？

用Gaussian优化含I的分子总是出错，请各位高人帮帮忙，看看怎么能优化出来
gjf文件
# mp2(full)/gen opt(maxcycle=400) freq scf(maxcycle=400)

Title Card Required

0 1
I                 -0.94029852    2.25373131    0.00000000
H                 -2.57029852    2.25373131    0.00000000

H     0
aug-cc-pvdz
****
I     0
S   6   1.00
      2.449790E+03           4.190000E-04     
      3.598080E+02           2.240000E-03     
      1.440580E+01           3.972230E-01     
      9.076320E+00          -9.322490E-01     
      2.088100E+00           9.371380E-01     
      1.034980E+00           3.920860E-01     
S   6   1.00
      2.449790E+03           1.750000E-04     
      3.598080E+02           1.057000E-03     
      1.440580E+01           1.690000E-01     
      9.076320E+00          -4.217930E-01     
      2.088100E+00           6.388640E-01     
      1.034980E+00           3.201150E-01     
S   1   1.00
      3.162840E-01           1.0000000        
S   1   1.00
      1.217190E-01           1.0000000        
S   1   1.00
      4.200000E-02           1.0000000        
P   5   1.00
      1.953010E+01           5.893400E-02     
      1.108820E+01          -2.309300E-01     
      2.715630E+00           6.648010E-01     
      1.204300E+00           4.506730E-01     
      3.399450E-01           2.898000E-02     
P   5   1.00
      1.953010E+01          -1.883600E-02     
      1.108820E+01           8.000600E-02     
      2.715630E+00          -3.066520E-01     
      1.204300E+00          -1.475940E-01     
      3.399450E-01           6.075060E-01     
P   1   1.00
      1.108810E-01           1.0000000        
P   1   1.00
      3.380000E-02           1.0000000        
D   5   1.00
      4.547650E+01           4.266000E-03     
      1.319280E+01          -1.362500E-02     
      4.227410E+00           3.097560E-01     
      1.942800E+00           5.097720E-01     
      8.397710E-01           2.974610E-01     
D   1   1.00
      3.000000E-01           1.0000000        
D   1   1.00
      1.191000E-01           1.0000000        
****

out出错信息
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RHF) =  -1631.30068077     A.U. after  401 cycles
             Convg  =    0.4417D+00             -V/T =  9.3532
             S**2   =   0.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in e:\G03W\l502.exe at Sun Jul 08 16:59:17 2012.
Job cpu time:  0 days  0 hours  0 minutes 55.0 seconds.
File lengths (MBytes):  RWF=     11 Int=      0 D2E=      0 Chk=      4 Scr=      1
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-07-08 17:15:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家
★ ★ ★
gmy1990: 金币+3 2012-07-08 21:55:56
这样的基组显然不是全电子基组

1，s,p收缩函数太少
2, 最陡的s函数收缩因子太小
3, I原子全电子计算的总能量应该在-6000 Hartree
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

1maomao
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2012-07-08 21:32:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

Siriusgoul

铜虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+2 2012-07-08 21:55:47
1maomao: 金币+5, ★有帮助 2012-07-09 21:52:06
1maomao: 金币+2, ★有帮助 2012-07-11 08:10:45
你这个是自洽场没有收敛的问题，可以在本论坛搜一搜解决方法。
但刚看过类似的问题在如下的帖子里
http://muchong.com/bbs/viewthread.php?tid=4696851&fpage=1
先可以在scf里加xqc等关键字试试，然后是再看看对碘原子的赝势处理对不对，一般碘这样的重原子容易出问题
赞一下 回复此楼
2楼2012-07-08 20:55:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nwwolfchj

金虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+3 2012-07-09 09:47:53
SCF不能收敛，若你输入正确则可能存在不稳定态，用ZMAT坐标或者目前的CART坐标改变一下结构先跳到“阱”内试一下，若还不行，若还不行 ，换MP2到其它更多包含关联能方法尝试一下，还不行就有可能是排斥态。

但从你的输入看，似乎有问题，I元素你用的为aug-cc-pVDZ的ECP基组，你只输入了价电子基组，对核电子的ECP没有输入，这样不会收敛的，即使你做了上述的尝试，把这段按照格式加在后边应该就没问题了，注意空格和****

I     0
I-ECP     4     28
g-ul potential
  1
2      1.00000000             0.00000000      
s-ul potential
  3
2     40.03337600            49.98964900      
2     17.30057600           281.00655600      
2      8.85172000            61.41673900      
p-ul potential
  4
2     15.72014100            67.41623900      
2     15.20822200           134.80769600      
2      8.29418600            14.56654800      
2      7.75394900            28.96842200      
d-ul potential
  4
2     13.81775100            35.53875600      
2     13.58780500            53.33975900      
2      6.94763000             9.71646600      
2      6.96009900            14.97750000      
f-ul potential
  4
2     18.52295000           -20.17661800      
2     18.25103500           -26.08807700      
2      7.55790100            -0.22043400      
2      7.59740400            -0.22164600
赞一下 回复此楼
4楼2012-07-09 09:22:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1maomao

金虫 (初入文坛)
送鲜花一朵
引用回帖:
3楼: Originally posted by beefly at 2012-07-08 21:32:47
这样的基组显然不是全电子基组

1，s,p收缩函数太少
2, 最陡的s函数收缩因子太小
3, I原子全电子计算的总能量应该在-6000 Hartree

改变收缩因子对计算结果影响太大，谢谢指导。
赞一下 回复此楼
5楼2012-07-09 19:41:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭