计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2868|回复 Reply: 13
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] 新手求助,orca利用多参考计算能量

[复制链接 Copy URL]

51

帖子

0

威望

657

eV
积分
708

Level 4 (黑子)

已毕业

本帖最后由 lm1997 于 2022-6-18 17:22 编辑

各位老师们好,小生最近学习使用orca用多参考态计算能量(某烷加氧过程),状况频出,特来请教。
我用casscf对反应分子进行了柔性扫描,原想利用mrci计算下能量一直遇到不收敛的情况,就改换成了NEVPT2,但是遇到计算很快的问题大概几十秒就完成了计算,我导估计每个点大概需要几个小时的计算,所以我觉得可能是我输入文件的错误,故向各位老师请教一下,希望各位老师批评指正。输入文件如下:




18.1.out

398.74 KB, 下载次数 Times of downloads: 8

1.3.out

430.04 KB, 下载次数 Times of downloads: 2

51

帖子

0

威望

657

eV
积分
708

Level 4 (黑子)

已毕业

14#
 楼主 Author| 发表于 Post on 2022-6-18 21:04:30 | 只看该作者 Only view this author
ionexchangeC 发表于 2022-6-18 17:38
为什么mrci计算时不用ORCA5?ORCA5的CASSCF能收敛。
另外,做NEVPT2计算,基组至少3zeta

是用的 orca5.01算的 大概我有三个结构只有一个结构能算出来,其他的算到中间就卡住了,需要算dz和tz之后外推到完全基组

218

帖子

2

威望

2039

eV
积分
2297

Level 5 (御坂)

13#
发表于 Post on 2022-6-18 17:38:18 | 只看该作者 Only view this author
lm1997 发表于 2022-6-18 17:22
输出文件 已经上传 您看一下

为什么mrci计算时不用ORCA5?ORCA5的CASSCF能收敛。
另外,做NEVPT2计算,基组至少3zeta

51

帖子

0

威望

657

eV
积分
708

Level 4 (黑子)

已毕业

12#
 楼主 Author| 发表于 Post on 2022-6-18 17:22:29 | 只看该作者 Only view this author
ionexchangeC 发表于 2022-6-18 16:35
mrci的不收敛是指casscf步不收敛还是后面的mrci步不收敛?

输出文件 已经上传 您看一下

218

帖子

2

威望

2039

eV
积分
2297

Level 5 (御坂)

11#
发表于 Post on 2022-6-18 16:35:57 | 只看该作者 Only view this author
mrci的不收敛是指casscf步不收敛还是后面的mrci步不收敛?

4291

帖子

4

威望

9568

eV
积分
13939

Level 6 (一方通行)

MOKIT开发者

10#
发表于 Post on 2022-6-18 16:30:28 | 只看该作者 Only view this author
计算时间没问题。记得检查CASSCF活性轨道(算前、算后)合理性,不要当成HF/MP2那样黑箱式地算下去。
自动做多参考态计算的程序MOKIT

51

帖子

0

威望

657

eV
积分
708

Level 4 (黑子)

已毕业

9#
 楼主 Author| 发表于 Post on 2022-6-18 16:21:24 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-18 16:10
所有计算步骤都正常结束了,结果确实是可以用的。
不过既然SC-NEVPT2做得这么快的话,建议用更可靠的FIC ...

好的老师,我去学习一下,十分感谢各位老师的帮助,解决小生心中的疑惑。

1万

帖子

0

威望

9902

eV
积分
22161

Level 6 (一方通行)

8#
发表于 Post on 2022-6-18 16:10:47 | 只看该作者 Only view this author

所有计算步骤都正常结束了,结果确实是可以用的。
不过既然SC-NEVPT2做得这么快的话,建议用更可靠的FIC-NEVPT2,SC-NEVPT2不满足轨道酉不变性,做势能面可能出现势能面不连续的问题
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

51

帖子

0

威望

657

eV
积分
708

Level 4 (黑子)

已毕业

7#
 楼主 Author| 发表于 Post on 2022-6-18 15:59:31 | 只看该作者 Only view this author
hebrewsnabla 发表于 2022-6-18 15:56
那很正常,就是这么快

第一次做这个,我觉得这种方法起码都要几十分钟起步的,竟然是秒级别的,超出我的认知,谢谢老师

51

帖子

0

威望

657

eV
积分
708

Level 4 (黑子)

已毕业

6#
 楼主 Author| 发表于 Post on 2022-6-18 15:58:25 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-18 15:55
NEVPT2本来就是远快于MRCI的。如果是比CASSCF还快,那才是有问题。
如果还不确定,贴完整输出文件上来

老师已经上传

928

帖子

1

威望

8273

eV
积分
9221

Level 6 (一方通行)

5#
发表于 Post on 2022-6-18 15:56:54 | 只看该作者 Only view this author
那很正常,就是这么快

1万

帖子

0

威望

9902

eV
积分
22161

Level 6 (一方通行)

4#
发表于 Post on 2022-6-18 15:55:35 | 只看该作者 Only view this author
NEVPT2本来就是远快于MRCI的。如果是比CASSCF还快,那才是有问题。
如果还不确定,贴完整输出文件上来
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

51

帖子

0

威望

657

eV
积分
708

Level 4 (黑子)

已毕业

3#
 楼主 Author| 发表于 Post on 2022-6-18 15:52:24 | 只看该作者 Only view this author
zjxitcc 发表于 2022-6-18 15:51
体系多少原子、多少基函数

C6H3O3

4291

帖子

4

威望

9568

eV
积分
13939

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2022-6-18 15:51:06 | 只看该作者 Only view this author
体系多少原子、多少基函数
自动做多参考态计算的程序MOKIT

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-28 15:43 , Processed in 0.187191 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list