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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 用Gaussian03或09，MP2(full)/aug-cc-pvdz如何优化HI？
51/1返回列表
查看: 2108  |  回复: 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的回帖

nwwolfchj

金虫 (正式写手)

【答案】应助回帖

%mem=1000MB
%nproc=2
# mp2(full)/genecp 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   9   1.00
      5.546500E+03           1.560000E-04     
      8.382140E+02           9.860000E-04     
      1.821870E+02           2.792000E-03     
      3.121230E+01          -4.325100E-02     
      1.953140E+01           2.341340E-01     
      8.240990E+00          -7.509430E-01     
      2.194550E+00           8.829680E-01     
      1.109110E+00           4.620610E-01     
      8.106100E-02           1.102000E-03     
S   9   1.00
      5.546500E+03          -7.300000E-05     
      8.382140E+02          -5.080000E-04     
      1.821870E+02          -1.158000E-03     
      3.121230E+01           1.219300E-02     
      1.953140E+01          -8.785400E-02     
      8.240990E+00           3.382000E-01     
      2.194550E+00          -5.765500E-01     
      1.109110E+00          -4.092980E-01     
      8.106100E-02           1.432310E-01     
S   1   1.00
      3.746410E-01           1.0000000        
S   1   1.00
      1.770800E-01           1.0000000        
S   1   1.00
      1.001000E-01           1.0000000        
S   1   1.00
      4.120000E-02           1.0000000        
P   9   1.00
      1.889880E+02           5.850000E-04     
      2.128680E+01           3.692300E-02     
      1.003960E+01          -2.353240E-01     
      3.451800E+00           3.414830E-01     
      1.974560E+00           5.347880E-01     
      1.024200E+00           2.651410E-01     
      4.494370E-01           2.578700E-02     
      1.866480E-01           5.220000E-04     
      7.348100E-02           6.060000E-04     
P   9   1.00
      1.889880E+02          -2.560000E-04     
      2.128680E+01          -1.168200E-02     
      1.003960E+01           8.319200E-02     
      3.451800E+00          -1.569700E-01     
      1.974560E+00          -2.245180E-01     
      1.024200E+00          -1.144510E-01     
      4.494370E-01           3.753560E-01     
      1.866480E-01           5.751360E-01     
      7.348100E-02           2.459170E-01     
P   1   1.00
      5.981000E-01           1.0000000        
P   1   1.00
      9.618000E-02           1.0000000        
P   1   1.00
      3.590000E-02           1.0000000        
D   7   1.00
      1.326620E+02           5.720000E-04     
      3.760540E+01           4.402000E-03     
      1.038910E+01          -4.092200E-02     
      6.490170E+00           9.966100E-02     
      3.454510E+00           3.226630E-01     
      1.844130E+00           4.003430E-01     
      9.624780E-01           2.683060E-01     
D   1   1.00
      4.728530E-01           1.0000000        
D   1   1.00
      1.932000E-01           1.0000000        
D   1   1.00
      7.820000E-02           1.0000000        
F   1   1.00
      4.064000E-01           1.0000000        
F   1   1.00
      1.848000E-01           1.0000000        
****

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
赞一下 回复此楼
8楼2012-07-10 14:54:07
已阅   回复此楼   关注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的回帖

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的回帖

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的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +15 木水思豆 2026-08-25 17/850 2026-08-25 23:24 by 实验小白ha
[基金申请] 今天务委会开完了，明天出结果吗 +17 angus9576 2026-08-25 21/1050 2026-08-25 22:39 by 悠悠美女
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗？20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭