计算化学公社

标题: orca做casscf计算时报错 [打印本页]

作者
Author: 澄海霸王龙 时间: 2023-8-6 00:25
标题: orca做casscf计算时报错
各位老师好！
最近学习了一下如何在orca上进行casscf的计算，我的体系是一个Fe-bpy体系的催化剂分子，以下是我的输入文件以及报错信息，请各位老师批评指导！

! def2-TZVP def2/JK Conv Patom RIJK
%maxcore 15000
%pal nprocs 10 end
* xyz 1 4
Fe             -0.00001500 1.08294000 0.00000800
N                1.90201300 1.00290800 -0.89363900
N                0.95185800 -0.68719400 1.05611000
...
O                -0.67967100 2.02672300 -1.13886400
*

%casscf
nel 3
norb 5
nroots 1
mult 2,6
end

作者
Author: Satoru 时间: 2023-8-6 03:59
3个电子为什么能产生6重态？

有个语法我不太确定，nroots一行可能需要写"nroots 1,1"

作者
Author: zjxitcc 时间: 2023-8-6 09:10
可以看此贴中我在9L的回复，自动CASSCF计算，也可以调用ORCA。http://bbs.keinsci.com/thread-38828-1-1.html

作者
Author: 澄海霸王龙 时间: 2023-8-6 09:51

Satoru 发表于 2023-8-6 03:59
3个电子为什么能产生6重态？

有个语法我不太确定，nroots一行可能需要写"nroots 1,1"

谢谢老师，因为氧上也有两个自由基，我尝试一下

作者
Author: 澄海霸王龙 时间: 2023-8-6 09:53

zjxitcc 发表于 2023-8-6 09:10
可以看此贴中我在9L的回复，自动CASSCF计算，也可以调用ORCA。http://bbs.keinsci.com/thread-38828-1-1.ht ...

谢谢老师，只是该楼中的人报错信息跟我不一样，都可以用MoKIT解决吗

作者
Author: zjxitcc 时间: 2023-8-6 09:58

澄海霸王龙发表于 2023-8-6 09:53
谢谢老师，只是该楼中的人报错信息跟我不一样，都可以用MoKIT解决吗

MOKIT有自己的严谨计算流程，调用ORCA一般不会报错，无需解决报错信息。

作者
Author: hebrewsnabla 时间: 2023-8-6 10:10
活性空间内只有3个电子的话，六重态是不可能的。

作者
Author: 澄海霸王龙 时间: 2023-8-6 17:10

hebrewsnabla 发表于 2023-8-6 10:10
活性空间内只有3个电子的话，六重态是不可能的。

老师您好，我这个是4重态，活性空间是3个电子也是报的同样的错误

作者
Author: 澄海霸王龙 时间: 2023-8-6 17:11

zjxitcc 发表于 2023-8-6 09:58
MOKIT有自己的严谨计算流程，调用ORCA一般不会报错，无需解决报错信息。

好的好的，谢谢老师

作者
Author: hebrewsnabla 时间: 2023-8-6 17:26

澄海霸王龙发表于 2023-8-6 17:10
老师您好，我这个是4重态，活性空间是3个电子也是报的同样的错误

你写的mult 2,6 代表casscf计算二重态和六重态，根本没有计算四重态。

作者
Author: 澄海霸王龙 时间: 2023-8-7 11:19

hebrewsnabla 发表于 2023-8-6 17:26
你写的mult 2,6 代表casscf计算二重态和六重态，根本没有计算四重态。

原来是这样，谢谢老师，我将mult修改为2,4重新提交了，谢谢老师

作者
Author: 澄海霸王龙 时间: 2023-8-7 16:30

zjxitcc 发表于 2023-8-6 09:58
MOKIT有自己的严谨计算流程，调用ORCA一般不会报错，无需解决报错信息。

老师，我按照您的教程用Mokit进行计算，但是出现了这个报错，以下是我的报错信息以及输入文件，请老师帮我看看哪里有问题

作者
Author: hebrewsnabla 时间: 2023-8-7 16:35

澄海霸王龙发表于 2023-8-7 16:30
老师，我按照您的教程用Mokit进行计算，但是出现了这个报错，以下是我的报错信息以及输入文件，请老师帮 ...

你在用高斯运行automr的输入文件？不是这样的，需要用automr运行。

作者
Author: 澄海霸王龙 时间: 2023-8-7 16:52

hebrewsnabla 发表于 2023-8-7 16:35
你在用高斯运行automr的输入文件？不是这样的，需要用automr运行。

啊？我是按照那个老师给的那个教程写的输入文件，是哪里错了吗，请老师指教

作者
Author: hebrewsnabla 时间: 2023-8-7 16:58
本帖最后由 hebrewsnabla 于 2023-8-7 17:00 编辑

澄海霸王龙发表于 2023-8-7 16:52
啊？我是按照那个老师给的那个教程写的输入文件，是哪里错了吗，请老师指教

那个帖子写了如何运行automr的输入文件。

automr Ir.gjf >Ir.out 2>&1 &

复制代码

当然，你得先安装MOKIT及其依赖。

作者
Author: 澄海霸王龙 时间: 2023-8-7 17:27

hebrewsnabla 发表于 2023-8-7 16:58
那个帖子写了如何运行automr的输入文件。

是的是的，我忘记把那一行加入脚本了，谢谢老师指正

作者
Author: zjxitcc 时间: 2023-8-7 19:02

澄海霸王龙发表于 2023-8-7 16:52
啊？我是按照那个老师给的那个教程写的输入文件，是哪里错了吗，请老师指教

MOKIT的automr用的也是高斯gjf文件，这是为了方便用户没有学习成本，而且GaussView也可以直接打开。实际运行时是用automr，而非g16。

作者
Author: 澄海霸王龙 时间: 2023-8-7 19:27
本帖最后由澄海霸王龙于 2023-8-7 20:04 编辑

zjxitcc 发表于 2023-8-7 19:02
MOKIT的automr用的也是高斯gjf文件，这是为了方便用户没有学习成本，而且GaussView也可以直接打开。实际 ...

了解了老师。谢谢老师指点。我还有一个疑问，我计算后在生成casscf自然轨道时脚本报了一个错误是说EOF while scanning triple-quoted string literal，但是我找不到不成对的双引号，请问如何解决呢

作者
Author: zjxitcc 时间: 2023-8-7 20:01

澄海霸王龙发表于 2023-8-7 19:27
了解了老师。谢谢老师指点。我还有一个疑问，我计算后报了一个错误是说EOF while scanning triple-quoted ...

最好截图多展示几行报错信息

作者
Author: 澄海霸王龙 时间: 2023-8-7 20:06

zjxitcc 发表于 2023-8-7 20:01
最好截图多展示几行报错信息

老师，我已经修改了帖子，麻烦您看一下

作者
Author: zjxitcc 时间: 2023-8-7 21:50
本帖最后由 zjxitcc 于 2023-8-7 21:51 编辑

澄海霸王龙发表于 2023-8-7 20:06
老师，我已经修改了帖子，麻烦您看一下

gjf文件内容有误。def2TZVP对于Fe和O都是全电子基组，没有赝势，所以你应该写gen，而非genecp；坐标底下只能给基组信息，不要写赝势信息（你对Fe和O试图定义def2赝势，这是错误的）。

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