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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 用Gaussian03或09，MP2(full)/aug-cc-pvdz如何优化HI？
51/1返回列表
查看: 2141  |  回复: 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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 mibUvS8DDCwf 2026-09-15 3/150 2026-09-16 23:28 by Tql5LQhh5rLK
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-16 22:40 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-16 20:50 by MHf40yA0Uhyb
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +8 QUjhNVAcOSff 2026-09-13 11/550 2026-09-16 18:37 by rRvxEj3zorCy
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:15 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:03 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 16:51 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-16 06:48 by LH5NkK5Ud4bh
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:36 by LH5NkK5Ud4bh
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-15 14:49 by VMDqgnlKimtZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭