版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

1maomao

金虫 (初入文坛)

应助: 2 (幼儿园)
金币: 1475.7
帖子: 27
在线: 22小时
虫号: 1475787
注册: 2011-11-04
性别: MM
专业: 理论和计算化学

[求助] 用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

回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:

新手请教：turbomole的基组如何选择？已经有5人回复
柔性扫描两个H2分子之间的势能曲线，该怎么写gauss配置文件？已经有7人回复
用MP2或者MP4进行单点能以及opt报错已经有4人回复
Pb2H4优化选择什么基组合适已经有10人回复
分子优化总是出错，请大家帮忙已经有9人回复
高斯计算分子间相互作用已经有11人回复
gaussian不同半经验优化得到的能量相差较大已经有4人回复
cc-pvdz-pp基集的使用已经有7人回复
急求Gaussian 03或者09软件，已经有4人回复
Gaussian 03优化结构已经有12人回复
Gaussian.03W-for.VC 已经有200人回复
利用DFT-D进行Gaussian优化，输入文件应该怎么设置呢已经有4人回复
电子轨道能量用什么方法计算比较好？已经有6人回复
this version requires full pseudpotential generation information 这是什么情况？已经有5人回复
农药增效剂有哪些，天然化合物和化学合成的都行已经有6人回复
全双工通信 full duplex在未来有前途么？可以作为研究方向么？已经有4人回复
gaussian 03和gaussian 09 优化分子的区别已经有5人回复
vasp优化结构中出现的问题已经有10人回复
求问一般优化小分子结构用什么方法和机组已经有10人回复
【讨论】用gaussian进行SCS-mp2计算已经有12人回复

1楼2012-07-08 17:15:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Siriusgoul

铜虫 (初入文坛)

应助: 10 (幼儿园)
金币: 266
散金: 200
帖子: 48
在线: 27.4小时
虫号: 1303115
注册: 2011-05-22
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

专家顾问 (职业作家)

地沟油冶炼专家

专家经验: +457
QC强帖: 12
应助: 86 (初中生)
金币: 12719.4
散金: 27307
红花: 300
帖子: 4663
在线: 925.3小时
虫号: 408372
注册: 2007-06-21
性别: GG
专业: 理论和计算化学
管辖: 计算模拟

★ ★ ★
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

金虫 (正式写手)

应助: 21 (小学生)
金币: 3200
散金: 35
红花: 8
帖子: 839
在线: 564.2小时
虫号: 428469
注册: 2007-08-04
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +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

金虫 (初入文坛)

应助: 2 (幼儿园)
金币: 1475.7
帖子: 27
在线: 22小时
虫号: 1475787
注册: 2011-11-04
性别: MM
专业: 理论和计算化学

送鲜花一朵

引用回帖:

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

1maomao

金虫 (初入文坛)

应助: 2 (幼儿园)
金币: 1475.7
帖子: 27
在线: 22小时
虫号: 1475787
注册: 2011-11-04
性别: MM
专业: 理论和计算化学

引用回帖:

4楼: Originally posted by nwwolfchj at 2012-07-09 09:22:52
SCF不能收敛，若你输入正确则可能存在不稳定态，用ZMAT坐标或者目前的CART坐标改变一下结构先跳到“阱”内试一下，若还不行，若还不行，换MP2到其它更多包含关联能方法尝试一下，还不行就有可能是排斥态。

但从 ...

帅哥，我刚学量化不久，不太会改格式，改了Gaussian 还是算不过去，麻烦您帮我看看哪里有问题。谢谢
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
S 1 1.00
   1.00000000          0.00000000
S 3 1.00
   40.03337600          49.98964900
   17.30057600          281.00655600
   8.85172000          61.41673900
P 4 1.00
   15.72014100          67.41623900
   15.20822200          134.80769600
   8.29418600          14.56654800
   7.75394900          28.96842200
D 4 1.00
   13.81775100          35.53875600
   13.58780500          53.33975900
   6.94763000          9.71646600
   6.96009900          14.97750000
F 4 1.00
   18.52295000          -20.17661800
   18.25103500          -26.08807700
   7.55790100          -0.22043400
   7.59740400          -0.22164600
****

赞一下

回复此楼

6楼2012-07-09 21:30:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

nwwolfchj

金虫 (正式写手)

应助: 21 (小学生)
金币: 3200
散金: 35
红花: 8
帖子: 839
在线: 564.2小时
虫号: 428469
注册: 2007-08-04
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
1maomao: 金币+10, ★★★★★最佳答案, 解决了我的燃眉之急 2012-07-11 11:29:14

引用回帖:

6楼: Originally posted by 1maomao at 2012-07-09 21:30:12
帅哥，我刚学量化不久，不太会改格式，改了Gaussian 还是算不过去，麻烦您帮我看看哪里有问题。谢谢
gjf文件

# mp2(full)/gen opt(maxcycle=400) freq scf(maxcycle=400)

Title Card Required

0 1
I ...

%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

赞一下

回复此楼

7楼2012-07-10 14:52:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

nwwolfchj

金虫 (正式写手)

应助: 21 (小学生)
金币: 3200
散金: 35
红花: 8
帖子: 839
在线: 564.2小时
虫号: 428469
注册: 2007-08-04
性别: GG
专业: 理论和计算化学

【答案】应助回帖

赞一下

回复此楼

8楼2012-07-10 14:54:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

nwwolfchj

金虫 (正式写手)

应助: 21 (小学生)
金币: 3200
散金: 35
红花: 8
帖子: 839
在线: 564.2小时
虫号: 428469
注册: 2007-08-04
性别: GG
专业: 理论和计算化学

最后结果

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

相关版块跳转我要订阅楼主 1maomao 的主题更新

返回列表