用MOKIT程序生成的输入文件跑Molpro报错

sxj-111 · 发表于 Post on 2024-5-13 16:17:44

Molpro报错如下：
? Error ? Cannot assign closed-shell space. Please used closed directive
? The problem occurs in find_scfcls

ts12_s0001.inp (13.52 KB, 下载次数 Times of downloads: 5)

ts12_s0001.out (26.79 KB, 下载次数 Times of downloads: 5)

zjxitcc · 发表于 Post on 2024-5-13 17:08:04

能否把你的gjf文件展示一下？不然不确定你想要做什么计算。

sxj-111 · 发表于 Post on 2024-5-13 17:48:31

zjxitcc 发表于 2024-5-13 17:08
能否把你的gjf文件展示一下？不然不确定你想要做什么计算。

老师好，是不是用Gaussian算ROHF的gjf文件？上面的inp文件是Molpro的输入文件，我想用Molpro算这个结构的单点能，结果不收敛，就采取了fch2com的方法。

zjxitcc · 发表于 Post on 2024-5-13 19:31:34

sxj-111 发表于 2024-5-13 17:48
老师好，是不是用Gaussian算ROHF的gjf文件？上面的inp文件是Molpro的输入文件，我想用Molpro算这个结构的 ...

哦，我以为你是用MOKIT的autosr自动做的单参考计算。看你的意思，没用到autosr，是自己做了一个ROHF计算然后用fch2com传轨道给Molpro，然后在输入文件末尾加关键词，也可以的。生成的com文件本身没什么问题，是不同版本的Molpro对输入信息要求不同，Molpro 2015要求这么写
{HF;start,2100.2;closed,33}
就是加上双占据轨道的数目，对于这个体系而言是33个。但更高版本的Molpro就没有这种要求。

你使用的是非常旧的MOKIT版本，不建议继续使用。早在v1.2.6rc16的版本中就会在这里加上closed,33。建议及时更新到最新版（目前是v1.2.6rc30）。

sxj-111 · 发表于 Post on 2024-5-15 20:26:01

zjxitcc 发表于 2024-5-13 19:31
哦，我以为你是用MOKIT的autosr自动做的单参考计算。看你的意思，没用到autosr，是自己做了一个ROHF计算 ...

感谢老师的解答，加上之后可以进行计算了。另外，如何查看双占据轨道的数目呀

，再次谢过老师了。

zjxitcc · 发表于 Post on 2024-5-15 20:40:45

sxj-111 发表于 2024-5-15 20:26
感谢老师的解答，加上之后可以进行计算了。另外，如何查看双占据轨道的数目呀，再次谢过老师了 ...

以ROHF为例，假设体系自旋多重度为3，则有2个单占据轨道，也即有2个单电子；再假设体系有50个电子，由于2个单占据轨道已经占用了2个单电子，还剩48个电子，那就是24个双占据轨道。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[辅助/分析程序] 用MOKIT程序生成的输入文件跑Molpro报错