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

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

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

回复此楼

» 猜你喜欢

298分 070300求调剂已经有15人回复
湖南大学张世国团队招26年第二批博士—离子液体、高分子与凝胶材料已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有179人回复
中国电机工程学报的投稿系统已经有0人回复
湖南农业大学能源材料创新团队招收2025化学专业研究生（学硕）调剂已经有0人回复
2026年西华大学材料学院-材料近净成形与表面技术研究团队-招收调剂研究生已经有9人回复
调剂信息已经有0人回复
河北农大生物与医药专业招收调剂，名额充足已经有0人回复
齐鲁工业大学国家重点实验室-新能源、锂离子电池、材料、化学方向已经有0人回复
杰青/长江团队招收调剂已经有0人回复
中南林业科技大学罗勇锋教授团队2026年招收材料类调剂研究生已经有0人回复

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

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

新手请教：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

专家顾问 (职业作家)

地沟油冶炼专家

专家经验: +458
QC强帖: 12
应助: 88 (初中生)
金币: 13045.7
散金: 27307
红花: 300
帖子: 4673
在线: 927.1小时
虫号: 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 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 求调剂，一志愿厦门大学，生物与医药，总分272，本科211 +3	Electron1cc 2026-04-01	4/200	2026-04-05 20:24 by lys0704
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4	cs0106 2026-04-04	4/200	2026-04-05 18:46 by imissbao
[考研] 304求调剂（085602，过四级，一志愿985） +15	化工人999 2026-04-04	15/750	2026-04-05 16:25 by 我是电风扇r
[考研] 367求调剂 +3	谢28 2026-03-30	3/150	2026-04-05 13:27 by huangmoli
[考研] 电子信息调剂交叉学科有推荐吗 +6	jhtfeybgj 2026-04-01	9/450	2026-04-05 11:13 by 猪会飞
[考研] 302求调剂一志愿华中师范大学 +8	小江小江江江 2026-04-02	8/400	2026-04-04 19:50 by 蓝云思雨
[考研] 考研调剂 +4	美丽的youth_ 2026-04-04	5/250	2026-04-04 17:16 by imissbao
[考研] 319求调剂 +4	星星不眨眼喽 2026-04-03	4/200	2026-04-04 16:25 by 中飞院空管学院�
[考研] 考研调剂 +3	15615482637 2026-04-03	3/150	2026-04-03 22:50 by ms629
[考研] 求调剂机会 +5	意染ivy 2026-04-03	5/250	2026-04-03 15:13 by qoooooo614
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +9	哇呼哼呼哼 2026-04-02	9/450	2026-04-03 12:05 by 1753564080
[考研] 325分化学调剂 +5	15771691647 2026-04-02	5/250	2026-04-03 09:58 by ChemPharm
[考研] 专硕 351 086100 也是考的材科基本科也是材料 +8	202451007219 2026-04-02	8/400	2026-04-03 09:50 by 蓝云思雨
[考研] 考研调剂 +3	李木子0120 2026-04-02	5/250	2026-04-02 21:45 by dongzh2009
[考研] 一志愿武汉理工0856，初试334 +3	26考研材料 2026-04-02	3/150	2026-04-02 21:22 by dongzh2009
[考研] 348求调剂 +11	zzzzyk123 2026-04-01	11/550	2026-04-02 16:52 by Wang200018
[考研] 能源动力调剂 +3	不破不立0 2026-04-02	3/150	2026-04-02 12:46 by ffffjjjj
[考研] 材料专硕322分 +11	哈哈哈吼吼吼哈 2026-04-01	11/550	2026-04-02 10:52 by lnilvy
[考研] 302求调剂一志愿北航070300，本科郑大化学 +8	圣日耳曼条 2026-04-01	11/550	2026-04-02 07:40 by chemdavid
[考研] 085601英二数二求调剂总分325 +4	余航航 2026-03-31	4/200	2026-03-31 17:38 by 唐沐儿