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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » Molpro 优化问题
北京石油化工学院2026年研究生招生接收调剂公告
21/1返回列表
查看: 1398  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +5 辰沐5211314 2026-03-26 5/250 2026-03-28 01:48 by 飞行日记西
[考研] 化学调剂 +4 爱吃番茄的旭 2026-03-24 5/250 2026-03-27 17:50 by kiokin
[考研] 材料与化工085600，总分304，本科有两篇sci参与，求调剂 +10 幸运的酱酱 2026-03-22 12/600 2026-03-27 16:08 by muchong357
[考研] 085600材料与化工调剂 +10 A-哆啦Z梦 2026-03-23 16/800 2026-03-27 15:13 by caszguilin
[考研] 308求调剂 +7 墨墨漠 2026-03-25 7/350 2026-03-27 14:47 by 狂炫麦当当
[考研] 0703一志愿9，初试成绩：338，四六级已过，有科研经历，求调剂！ +3 Zuhui0306 2026-03-25 3/150 2026-03-27 14:09 by shangxh
[考研] 085601 材料工程 313分 求调剂 +5 Ong3 2026-03-27 5/250 2026-03-27 12:24 by goldfish51
[考研] 085600材料与化工306 +10 z1z2z3879 2026-03-21 11/550 2026-03-27 11:31 by wangjy2002
[考研] 324求调剂 +8 hanamiko 2026-03-26 10/500 2026-03-27 08:06 by hypershenger
[考研] 329求调剂 +7 钮恩雪 2026-03-25 7/350 2026-03-27 04:28 by wxiongid
[考研] 求调剂 一志愿 本科 北科大 化学 343 +6 13831862839 2026-03-24 7/350 2026-03-26 22:57 by 不吃魚的貓
[考研] 081700 调剂 267分 +11 迷人的哈哈 2026-03-23 11/550 2026-03-26 15:41 by zzll406
[考研] 一志愿天津大学339材料与化工求调剂 +3 江往卖鱼 2026-03-26 3/150 2026-03-26 09:42 by 王小欠i
[考研] 332求调剂 +6 032500 2026-03-25 6/300 2026-03-25 22:45 by 418490947
[考研] 求调剂 +3 QiMing7 2026-03-25 3/150 2026-03-25 21:13 by 给你你注意休息
[考研] 302求调剂 +4 锦衣卫藤椒 2026-03-25 4/200 2026-03-25 16:29 by 功夫疯狂
[考研] 生物学学硕求调剂 +7 小羊睡着了? 2026-03-23 10/500 2026-03-25 02:24 by 清风拂扬。 m
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 一志愿深大，0703化学，总分302，求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
信息提示
关闭
请填处理意见
关闭