计算化学公社

标题: 求助一个molpro的错误 [打印本页]

作者
Author: lastzealot 时间: 2019-10-29 23:31
标题: 求助一个molpro的错误
高手们好，请教一个casscf报错问题，屡次出现，不知如何是好

计算一个CrCl6_3-的Oh对称性离 (, 下载次数 Times of downloads: 5)

子的casscf（实际计算中为D2h），优化和hf过程都没问题，但是在multi的时候出现了问题。
在hf计算中显示，三个单电子分别占据在三个简并轨道B1g/B2g/B3g中，LUMO为一组二重简并轨道Ag/Ag。
那么在一个(3,5)活性空间中，不改变自旋方向情况下共有10种态，分别为一个B1g*B2g*B3g=Ag，三个Ag*B1g*B2g=B3g，三个Ag*B1g*B3g=B2g，三个Ag*B2g*B3g=B1g

在molpro中Ag,B1g,B2g,B3g对应的不可约表示分别为1,4,6,7
但是仍然出错，找不到原因。
所以来请高手帮帮忙。

输入文件如下
***,crcl6
gprint,orbital=30
symmetry,x,y,z
angstrom
geometry=zmat
r(i)=2.4702767791
basis=cc-pvdz
{hf
wf,spin=3,charge=-3
save,2100.2
}
{multi
occ,38,22,22,9,11,9,9,0
closed,36,22,22,8,11,8,8,0
frozen,36,22,22,8,11,8,8,0
start,2100.2

wf,129,4,3
state,3
wf,129,6,3
state,4
wf,129,7,3
state,3

save,2140.2
}

e_cas(3,5)_trip=energy

put,molden,cas_35_triplet.molden

输出文件结尾
Micro-Iterations solved by the two-step solver

ITER. MIC NCI NEG ENERGY(VAR) ENERGY(PROJ) ENERGY CHANGE GRAD(0) GRAD(ORB) GRAD(CI) STEP TIME

MORE STATES THAN CONFIGURATIONS SPECIFIED: 0 3

谢谢大家

作者
Author: lastzealot 时间: 2019-10-29 23:31
插入文件位置有点问题，失误，见谅

作者
Author: sobereva 时间: 2019-10-30 04:11

lastzealot 发表于 2019-10-29 23:31
插入文件位置有点问题，失误，见谅

发现帖子有问题就直接在原帖上重新编辑，不要在回帖里说明

作者
Author: qinjiu 时间: 2019-10-30 16:15
MORE STATES THAN CONFIGURATIONS SPECIFIED: 0 3

你报错信息写着实际没那么多个态，你要算3个state，只有0个state
自旋3是四重态，你看看是不是自旋写错了

作者
Author: lastzealot 时间: 2019-10-30 17:35

qinjiu 发表于 2019-10-30 16:15
MORE STATES THAN CONFIGURATIONS SPECIFIED: 0 3

你报错信息写着实际没那么多个态，你要算3个state ...

首先谢谢你的回答
报错信息中第一个数字0指的是实际的态的数量？

crcl6体系为d3，也就是下面三个简并轨道各有一个电子

现在还是晕的，我是按照群表里的直积表算的呀，不知道应该把用对称性设置为几

您能否帮看下，state后面应该怎么写？

谢谢

作者
Author: qinjiu 时间: 2019-10-30 20:30

lastzealot 发表于 2019-10-30 17:35
首先谢谢你的回答
报错信息中第一个数字0指的是实际的态的数量？

state 的意思是做态平均，给出同样自旋、对称性、电子数的态，按能量由低到高的顺序给出。
你的输入卡里面写 wf,129,4,3 ;state3是让程序算 129个电子，B1g对称性，四重态的三个分子态

程序报错，说满足这个条件的态是0

你要研究一下输入卡哪里写错了，同样电子组态直积出来的态，比如N原子2s2 2p3 这个组态，可以直积出
各种谱项，不同多重度都有可能。

此外，如果你要算的体系基态如果是开壳层的，你需要用open卡去指定HF，不然算出来的是闭壳层的激发态

作者
Author: lastzealot 时间: 2019-10-30 21:49

qinjiu 发表于 2019-10-30 20:30
state 的意思是做态平均，给出同样自旋、对称性、电子数的态，按能量由低到高的顺序给出。
你的输入卡里 ...

谢谢朋友
我在琢磨一下

作者
Author: Yin 时间: 2021-7-3 16:30
楼主你好，请问你的问题解决了吗？我也遇到了这个错误

作者
Author: wsmrt33 时间: 2024-9-26 14:08
可以试试把 frozen,36,22,22,8,11,8,8,0 去掉
根据 molpro 手册，closed 对 frozen 是包含的关系：
(, 下载次数 Times of downloads: 0)
因此你frozen,36,22,22,8,11,8,8,0 这行应该是多余的。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)