计算化学公社

标题: 用MOKIT程序生成的输入文件跑Molpro报错 [打印本页]

作者
Author: sxj-111 时间: 2024-5-13 16:17
标题: 用MOKIT程序生成的输入文件跑Molpro报错
Molpro报错如下：
? Error ? Cannot assign closed-shell space. Please used closed directive
? The problem occurs in find_scfcls
(, 下载次数 Times of downloads: 5)

(, 下载次数 Times of downloads: 5)

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

作者
Author: sxj-111 时间: 2024-5-13 17:48

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

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

作者
Author: zjxitcc 时间: 2024-5-13 19:31

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）。

作者
Author: sxj-111 时间: 2024-5-15 20:26

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

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

，再次谢过老师了。

作者
Author: zjxitcc 时间: 2024-5-15 20:40

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

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

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