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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » Molpro 优化问题
21/1返回列表
查看: 1430  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xinzifish

金虫 (初入文坛)

[求助] 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里输入有问题，还是其它撒问题，请各位赐教！
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-09-24 21:50:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chembetsey

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 评审有感 +16 popular289 2026-05-18 27/1350 2026-05-24 17:34 by hhs666
[教师之家] 论文撤稿了 +4 bjvtcliu 2026-05-24 7/350 2026-05-24 17:29 by bjvtcliu
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 hvkbtfonbv 2026-05-23 4/200 2026-05-24 17:21 by 75ui6h7z2t
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:10 by 75ui6h7z2t
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:01 by 75ui6h7z2t
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 5/250 2026-05-24 16:21 by hhx1yx9evi
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 4/200 2026-05-24 16:16 by hhx1yx9evi
[基金申请] 河北省自然科学基金 +6 Peterchao 2026-05-18 9/450 2026-05-24 16:02 by 130067131
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 pmo95bazuy 2026-05-23 8/400 2026-05-24 15:56 by 1uy1ht2y9r
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！ +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[基金申请] 青B发送上会通知了吗 +5 chemBioBro 2026-05-22 7/350 2026-05-23 12:35 by zhuifengzhy
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[论文投稿] 投稿求助，期刊 +4 希冀，有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀，有书读
[基金申请] 面上本子正文33页，违规吗？会被低分嘛？ +14 1234567wang 2026-05-17 16/800 2026-05-21 17:58 by 脆脆的饼干
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[考博] 如果工作了想读博，可以边工作边读全日制嘛？ 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
信息提示
关闭
请填处理意见
关闭