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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » Molpro 优化问题
21/1返回列表
查看: 1394  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290分材料工程085601求调剂 数二英一 +3 llx0610 2026-03-02 3/150 2026-03-02 14:15 by yc258
[考研] 化工专硕348，一志愿985求调剂 +6 弗格个 2026-02-28 9/450 2026-03-02 14:09 by liyongv
[考研] 材料学硕318求调剂 +9 February_Feb 2026-03-01 9/450 2026-03-02 13:31 by njzyff
[考研] 272求调剂 +7 材紫有化 2026-02-28 7/350 2026-03-02 12:48 by 无际的草原
[考研] 085600 英一数二272求调剂 5+4 vida_a 2026-03-01 10/500 2026-03-02 12:33 by vida_a
[考研] 材料工程274求调剂 +4 Lilithan 2026-03-01 4/200 2026-03-02 12:06 by yuchj
[基金申请] 此成果不能导入原因：元数据必填信息不完整，可 进行补充。 +4 Kittylucky 2026-03-02 5/250 2026-03-02 11:07 by jurkat.1640
[考研] 284求调剂 +10 天下熯 2026-02-28 11/550 2026-03-02 11:03 by 无际的草原
[考研] 欢迎采矿、地质、岩土、计算机、人工智能等专业的同学报考 +5 pin8023 2026-02-28 7/350 2026-03-02 10:33 by ZY，先生
[考研] 材料类求调剂 +11 wana_kiko 2026-02-28 14/700 2026-03-02 08:46 by 聪明的大松鼠
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考研] 材料复试调剂 +4 学材料的点 2026-03-01 5/250 2026-03-02 08:26 by houyaoxu
[考研] 0857调剂 +4 一ll半 2026-02-28 5/250 2026-03-02 02:33 by 908055542
[考研] 290求调剂 +9 材料专硕调剂； 2026-02-28 11/550 2026-03-01 17:21 by sunny81
[考研] 285求调剂 +8 满头大汗的学生 2026-02-28 8/400 2026-03-01 16:47 by caszguilin
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[考研] 295复试调剂 +3 简木ChuFront 2026-03-01 3/150 2026-03-01 14:27 by zzxw520th
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[硕博家园] 2025届双非化工硕士毕业，申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
信息提示
关闭
请填处理意见
关闭