24小时热门版块排行榜

返回列表

xinzifish

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1141.5
帖子: 38
在线: 62.4小时
虫号: 536494
注册: 2008-03-31
专业: 理论和计算化学

[求助] Molpro 优化问题

我用Molpro 在MP2下优化碘苯这个分子，基组（H,C=cc-pVTZ, I=cc-pVTZ-PP), 具体输入文件如下：

***, optimize C2v-ground state PhI-S0 with cc-pVTZ-pp basis
memory, 500, m
geomtyp=xyz
geom={
12
C6H5I
C                0.00000000 0.00000000 -0.54455300
C                0.00000000 1.21269600 -1.23228600
C                0.00000000 1.20457900 -2.62540200
C                0.00000000 0.00000000 -3.32520200
C                0.00000000 -1.20457900 -2.62540200
C                0.00000000 -1.21269600 -1.23228600
I                0.00000000 0.00000000 1.53997500
H                0.00000000 2.14760700 -0.69036300
H                0.00000000 2.14437300 -3.16047300
H                0.00000000 0.00000000 -4.40619400
H                0.00000000 -2.14437300 -3.16047300
H                0.00000000 -2.14760700 -0.69036300
}

Basis={
ecp,I,28,4,3
1
2,1.000000,0.000000
3
2,40.033376,49.989649
2,17.300576,281.006556
2,8.851720,61.416739
4
2,15.720141,67.416239
2,15.208222,134.807696
2,8.294186,14.566548
2,7.753949,28.968422
4
2,13.817751,35.538756
2,13.587805,53.339759
2,6.947630,9.716466
2,6.960099,14.977500
4
2,18.522950,-20.176618
2,18.251035,-26.088077
2,7.557901,-0.220434
2,7.597404,-0.221646
4
2,15.720141,-134.832478
2,15.208222,134.807696
2,8.294186,-29.133096
2,7.753949,28.968422
4
2,13.817751,-35.538756
2,13.587805,35.559839
2,6.947630,-9.716466
2,6.960099,9.985000
4
2,18.522950,13.451079
2,18.251035,-13.044039
2,7.557901,0.146956
2,7.597404,-0.110823
! IODINE basis set
s,I,5.546500E+03,8.382140E+02,1.821870E+02,3.121230E+01,1.953140E+01,8.240990E+00,2.194550E+00,1.109110E+00,3.746410E-01,1.770800E-01,8.106100E-02,1.001000E-01;
c,1.11,1.560000E-04,9.860000E-04,2.792000E-03,-4.325100E-02,2.341340E-01,-7.509430E-01,8.829680E-01,4.620610E-01,2.228600E-02,-4.353000E-03,1.102000E-03;
c,1.11,-7.300000E-05,-5.080000E-04,-1.158000E-03,1.219300E-02,-8.785400E-02,3.382000E-01,-5.765500E-01,-4.092980E-01,5.674590E-01,6.124890E-01,1.432310E-01;
c,9.9,1.000000E+00;
c,10.10,1.000000E+00;
c,12.12,1.000000E+00;
p,I,1.889880E+02,2.128680E+01,1.003960E+01,3.451800E+00,1.974560E+00,1.024200E+00,4.494370E-01,1.866480E-01,7.348100E-02,5.981000E-01,9.618000E-02;
c,1.9,5.850000E-04,3.692300E-02,-2.353240E-01,3.414830E-01,5.347880E-01,2.651410E-01,2.578700E-02,5.220000E-04,6.060000E-04;
c,1.9,-2.560000E-04,-1.168200E-02,8.319200E-02,-1.569700E-01,-2.245180E-01,-1.144510E-01,3.753560E-01,5.751360E-01,2.459170E-01;
c,10.10,1.000000E+00;
c,11.11,1.000000E+00;

d,I,1.326620E+02,3.760540E+01,1.038910E+01,6.490170E+00,3.454510E+00,1.844130E+00,9.624780E-01,4.728530E-01,1.932000E-01;
c,1.9,5.720000E-04,4.402000E-03,-4.092200E-02,9.966100E-02,3.226630E-01,4.003430E-01,2.683060E-01,8.484700E-02,7.632000E-03;
c,8.8,1.000000E+00;
c,9.9,1.000000E+00;

f,I,4.064000E-01;
c,1.1,1.000000E+00;
H=cc-pVTZ
C=cc-pVTZ
}

wf, 94, 1, 0
hf
mp2
optg
frequencies

它可以计算hf, mp2的单点能量，可是到optg优化这步总是出现无故断掉，查看LOG文件也没给出任何出错信息，我也不知道是什么原因出错。后来换了其它方法优化例如B3LYP,也是一到OPTG这步就死掉，屏幕上就会出现GLOBAL ERROR fehler on processor 0 等错误，查看log文件，给出错误信息也就是：ERROR IN QEXIT: ACTIVE=TWOEL  CALLED WITH，不知道是什么原因造成的。即使只做HF优化，它也会说：Geometry written to block  1 of record 700
Norm of gradient contribution is huge!
Probably due to wrong coordinates.
可是即使换成Z-MATRIX的格式来做也是一样。而我在Gaussian09里却可以正常完成这个优化作业，不知道是我的MOLPRO里输入有问题，还是其它撒问题，请各位赐教！

回复此楼

» 猜你喜欢

英国利物浦大学屠昕和雍雪课题组招收CSC博士生已经有27人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有233人回复
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建已经有0人回复
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人已经有45人回复
英国贝尔法斯特女王大学招收 2026 年 CSC/BBSRC phd 机器学习催化剂设计方向已经有0人回复
法国里昂大学招生CSC 资助博士生已经有0人回复
法国里昂大学接收CSC 资助访问学者已经有0人回复
北京工业大学材料科学与工程学院《院士团队》招收工程专业博士研究生（计算优先）已经有1人回复
帮导师招2026CSC博士（巴塞罗那自治大学UAB-CSC博士项目）已经有3人回复

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

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

1楼 2013-09-24 21:50:25

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

chembetsey

木虫 (小有名气)

应助: 125 (高中生)
金币: 3820.1
散金: 5
红花: 9
帖子: 262
在线: 281.8小时
虫号: 1781694
注册: 2012-04-27
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
gmy1990: 金币+2 2013-09-30 22:07:59
xinzifish: 金币+10, ★有帮助, 谢谢回复，我后来也发现电子数目不对，改成了66，可是在G09和MOLPRO12算的能量还是有所差别，虽然用的是同样的构型和基组 2013-10-01 21:28:27

出现这种情况是因为你对I原子使用了cc-pVTZ-PP这个赝势基组，导致“总电子数”不再是94了。

fix的方法是不写“wf, 94, 1, 0”，让程序自己算。你会发现程序给出的电子数是66。(为什么是66呢？怎么少了94-66=28 呢？如果你了解赝势基组的意思的话，大概能猜到）

PS1：我用的Molpro里面已经自带cc-pVTZ-PP这个基组的信息了，所以不用写那么一大堆，直接I=cc-pVTZ-PP就好。
PS2：建议你把molpro算出的HF能量跟Gaussian的比较一下。如果你用WF,94,1,0选项算的话，结果应该跟高斯的差别很大（上百个hartree）
PS3：记得发金币给我！

赞一下(1人)

回复此楼

2楼2013-09-30 19:26:05

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

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

返回列表