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

24小时热门版块排行榜    

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

金虫 (正式写手)
最后结果

Predicted change in Energy=-4.623138D-13
Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
--------------------------                            --------------------------
! Name  Definition              Value          Derivative Info.                !
--------------------------------------------------------------------------------
! R1    R(1,2)                  1.6102         -DE/DX =    0.0                 !
--------------------------------------------------------------------------------
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

1|1|UNPC-CHJ-PC|Freq|RMP2-Full|GenECP|H1I1|CHJ|10-Jul-2012|0||#N Geom=
AllCheck Guess=TCheck SCRF=Check GenChk RMP2-Full/GenECP Freq||Title C
ard Required||0,1|I,-0.9501980107,2.25373131,0.|H,-2.5603990293,2.2537
3131,0.||Version=IA32W-G09RevA.02|State=1-SG|HF=-295.2438808|MP2=-295.
4938741|RMSD=3.419e-010|RMSF=2.582e-007|ZeroPoint=0.0054477|Thermal=0.
0078083|Dipole=-0.1931939,0.,0.|DipoleDeriv=-0.0256766,0.,0.,0.,-0.063
4913,0.,0.,0.,-0.0634913,0.0256766,0.,0.,0.,0.0634913,0.,0.,0.,0.06349
13|Polar=37.0356885,0.,34.7230088,0.,0.,34.7230088|PG=C*V [C*(H1I1)]|N
Imag=0||0.21637266,0.,0.00000015,0.,0.,0.00000015,-0.21637266,0.,0.,0.
21637266,0.,-0.00000015,0.,0.,0.00000015,0.,0.,-0.00000015,0.,0.,0.000
00015||0.00000045,0.,0.,-0.00000045,0.,0.|||@


TO SEE A WORLD IN A GRAIN OF SAND
AND A HEAVEN IN A WILD FLOWER
HOLD INFINITY IN THE PALM OF YOUR HAND
AND ETERNITY IN AN HOUR

                       -- WILLIAM BLAKE
Job cpu time:  0 days  0 hours  1 minutes 28.0 seconds.
File lengths (MBytes):  RWF=    286 Int=      0 D2E=      0 Chk=      1 Scr=      1
Normal termination of Gaussian 09 at Tue Jul 10 02:53:34 2012.
赞一下 回复此楼
9楼2012-07-10 14:55:02
 已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭