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

回复此楼

» 猜你喜欢

北京师范大学胡德平课题组以及袁萌伟课题组考研招生已经有0人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有223人回复
祈福论文小修已经有21人回复
香港理工大学计算材料与物理课题组招聘博后已经有15人回复
香港理工大学计算材料与物理课题组招聘博后已经有15人回复
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生已经有0人回复
英国利物浦大学屠昕和雍雪课题组招收CSC博士生已经有26人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）已经有0人回复
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人已经有41人回复

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

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

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

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

查看全部 9 个回答

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