| 查看: 1372 | 回复: 1 | ||
[求助]
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里输入有问题，还是其它撒问题，请各位赐教！ |
回复此楼» 猜你喜欢
如何掌握一门外语？
已经有25人回复
-
国家级人才团体课题组招收2026届博士
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有60人回复
-
气相色谱点不着火
已经有4人回复
-
化工安全有关的复习提纲
已经有0人回复
-
国家级人才课题组招收2026年入学博士
已经有2人回复
-
散金
已经有87人回复
-
浙江师范大学国家杰青杨启华教授团队招收2026年博士研究生
已经有34人回复
-
南通大学电气与自动化学院2026年优秀博士招聘公告
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 高精度量子化学计算软件Molpro 入门及理论精讲课程
已经有3人回复
- ANSYS单位问题以及参数化优化问题
已经有9人回复
- 固定键长优化问题
已经有7人回复
- 关于晶体优化问题
已经有6人回复
- 为什么meta-GGA一般不建议用来做结构优化？
已经有11人回复
- 优青应该是一场混战吧
已经有32人回复
- Chem3D中优化的问题
已经有6人回复
- 全国第二届理论与计算化学软件Molpro Molcas培训班
已经有19人回复
- molpro 结构优化CCSD不收敛，求指教
已经有4人回复
- vasp 结构优化问题
已经有14人回复
- MOLPRO优化计算为什么如此神奇？急啊！
已经有8人回复
- 结构优化中问题
已经有6人回复
- VASP优化结构遇到的问题。
已经有9人回复
- 【讨论】molpro, molcas和gaussian的优缺点
已经有3人回复
- 【求助】请教molpro计算中出现的一个问题
已经有7人回复
- 【求助】molpro or molcas??
已经有7人回复
- 【求助】Domol的优化问题
已经有10人回复
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
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：记得发金币给我！ |
2楼2013-09-30 19:26:05