计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: rayinyin
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] ORCA的CASSCF,NEVPT2以及MRCI的表现到底如何?

[复制链接 Copy URL]

27

帖子

0

威望

299

eV
积分
326

Level 3 能力者

16#
发表于 Post on 2023-5-18 20:57:47 | 只看该作者 Only view this author
wzkchem5 发表于 2023-5-18 18:27
注意光做CASSCF不够,必须考虑动态相关,例如用NEVPT2。你就把NEVPT2等多参考态方法看作是比DFT更准的方 ...

好的,多谢老师解答

27

帖子

0

威望

299

eV
积分
326

Level 3 能力者

17#
发表于 Post on 2023-5-24 20:26:30 | 只看该作者 Only view this author
wzkchem5 发表于 2023-5-18 18:27
注意光做CASSCF不够,必须考虑动态相关,例如用NEVPT2。你就把NEVPT2等多参考态方法看作是比DFT更准的方 ...

老师您好,请问能推荐一下适合入门CASSCF的参考书吗?

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

18#
发表于 Post on 2023-5-24 20:32:55 | 只看该作者 Only view this author
SHANY 发表于 2023-5-24 13:26
老师您好,请问能推荐一下适合入门CASSCF的参考书吗?

orca论坛(见我签名档链接)的下载区有个CASSCF tutorial,可以看一下。不仅适用于orca里的CASSCF计算,对其他程序也有借鉴意义
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

27

帖子

0

威望

299

eV
积分
326

Level 3 能力者

19#
发表于 Post on 2023-5-24 21:55:47 | 只看该作者 Only view this author
wzkchem5 发表于 2023-5-24 20:32
orca论坛(见我签名档链接)的下载区有个CASSCF tutorial,可以看一下。不仅适用于orca里的CASSCF计算, ...

好的,我直接的任务就是参考那个tutorial做的,我想问一下有没有一些讲原理的参考书。
另外,我还想问两个问题,在做CASSCF计算时,roots应该怎样选择呢;
以及,如果我想确定系统的自选多重度,我需要在不同的多重度下先做结构优化,再计算CASSCF吗,还是可以看着结构差不多,直接算两次CASSCF就行呢?

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

20#
发表于 Post on 2023-5-24 22:24:05 | 只看该作者 Only view this author
本帖最后由 wzkchem5 于 2023-5-24 15:31 编辑
SHANY 发表于 2023-5-24 14:55
好的,我直接的任务就是参考那个tutorial做的,我想问一下有没有一些讲原理的参考书。
另外,我还想问两 ...

原理方面看orca手册detailed documentation底下的CASSCF章节就行
roots的选择规则三言两语不好概括完全,CASSCF tutorial里讲了一些。大体来说,先把roots设成你感兴趣的态的数目(例如你只感兴趣基态,就是1;你如果要算激发态,就设成足够涵盖你研究的光谱范围的值),再验证这个roots是否满足如下条件:如果N个态互为(近)简并态,其中一个被你选进来了,那么剩下N-1个必然也被你选进来了。如果不满足,那么增加roots直至这个条件刚好满足为止。当需要研究一系列体系(例如同一个单分子反应的反应物、过渡态、中间体、产物,或者你要比较若干个类似分子的光谱)时,roots一般而言需要统一,这样理论级别才能统一,所以选择roots时要兼顾你研究的所有体系,让你设置的roots对你研究的所有体系都同时适用。除此之外还有很多额外规则,尤其计算双分子反应,或者需要考虑SOC,或者可能有intruder states时更是如此,我就不一一列举了。
确定自旋多重度,严格的做法是分别做结构优化。当且仅当你不做结构优化就知道结构优化的影响不大时,可以只用一个结构做CASSCF(这时CASSCF有两种做法,分别算两个自旋多重度也可以,把两种自旋多重度放在一起做态平均、一次出结果也可以,什么情况下用哪种方法在CASSCF tutorial里应该讲了,这个要展开讲又是要罗列一大堆特殊情况,我就不讲了)。此外在做计算之前,应该充分利用化学常识对所有可能的自旋多重度做初步筛选,把根据常识即可知道肯定不是基态的自旋多重度先去掉。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

27

帖子

0

威望

299

eV
积分
326

Level 3 能力者

21#
发表于 Post on 2023-5-24 22:28:04 | 只看该作者 Only view this author
wzkchem5 发表于 2023-5-24 22:24
原理方面看orca手册detailed documentation底下的CASSCF章节就行
roots的选择规则三言两语不好概括完全 ...

好的谢谢老师

27

帖子

0

威望

299

eV
积分
326

Level 3 能力者

22#
发表于 Post on 2023-8-8 16:55:59 | 只看该作者 Only view this author
wzkchem5 发表于 2023-5-24 22:24
原理方面看orca手册detailed documentation底下的CASSCF章节就行
roots的选择规则三言两语不好概括完全 ...

老师您好,我还想请教一下nroots 的问题。请问如果我想算基态,我可不可以先设一个比较大的root,(比如说20,但是未完全收敛),然后发现前五个态能量相差不大(~0.04eV),第六个态和第五个态差大概0.4eV,我可不可以认为前五个态可能存在简并,之后就只算五个态呢。

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

23#
发表于 Post on 2023-8-8 17:50:12 | 只看该作者 Only view this author
SHANY 发表于 2023-8-8 09:55
老师您好,我还想请教一下nroots 的问题。请问如果我想算基态,我可不可以先设一个比较大的root,(比如说 ...

可以,与此同时也检查一下电子态的成分,看看这些态是不是确实近似为一个简并态的各个简并分量
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

27

帖子

0

威望

299

eV
积分
326

Level 3 能力者

24#
发表于 Post on 2023-8-8 19:46:50 | 只看该作者 Only view this author
wzkchem5 发表于 2023-8-8 17:50
可以,与此同时也检查一下电子态的成分,看看这些态是不是确实近似为一个简并态的各个简并分量

谢谢老师。请问怎么确定他是不是一个简并态的简并分量呢?

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

25#
发表于 Post on 2023-8-8 20:04:20 | 只看该作者 Only view this author
SHANY 发表于 2023-8-8 12:46
谢谢老师。请问怎么确定他是不是一个简并态的简并分量呢?

先从理论上推导,假如这个体系的对称性是精确的,那么最低的那几个电子态的波函数应该长什么样。如果其中两个态A和B因对称性而严格简并,那么你一旦发现你实际计算的一个电子态的成分和A相符,你就知道必须得把另外一个成分和B相符的电子态也考虑进来
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

27

帖子

0

威望

299

eV
积分
326

Level 3 能力者

26#
发表于 Post on 2023-8-8 20:36:05 | 只看该作者 Only view this author
wzkchem5 发表于 2023-8-8 20:04
先从理论上推导,假如这个体系的对称性是精确的,那么最低的那几个电子态的波函数应该长什么样。如果其中 ...

好的,谢谢老师

588

帖子

0

威望

3304

eV
积分
3892

Level 5 (御坂)

27#
发表于 Post on 2024-3-4 11:50:10 | 只看该作者 Only view this author
SHANY 发表于 2023-5-24 20:26
老师您好,请问能推荐一下适合入门CASSCF的参考书吗?

您好,请问CASSCF tutorial可以共享一下吗?现在网站进不去

3806

帖子

4

威望

7994

eV
积分
11880

Level 6 (一方通行)

MOKIT开发者

28#
发表于 Post on 2024-3-4 11:59:40 | 只看该作者 Only view this author
XXX程序的功能好/坏/慢,这类问题如今基本失去了意义,因为有了MOKIT后,可以充分利用每个量化程序最强的、特色的功能(比如在PySCF里做CASSCF计算,传CAS轨道给Molpro/OpenMolcas/ORCA做CASPT2计算),对于每个多组态/多参考方法都可实现世上最快的计算,完全没有必要把自己局限在一个程序中。
自动做多参考态计算的程序MOKIT

3806

帖子

4

威望

7994

eV
积分
11880

Level 6 (一方通行)

MOKIT开发者

29#
发表于 Post on 2024-3-4 12:04:34 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2024-3-4 12:06 编辑
xujian 发表于 2024-3-4 11:50
您好,请问CASSCF tutorial可以共享一下吗?现在网站进不去

CASSCF-Tutorial-4.2.pdf (2.15 MB, 下载次数 Times of downloads: 37)

geometries-4.2.zip (11.13 KB, 下载次数 Times of downloads: 13)
(官网上文件就是只更新到4.2。顺带提醒一下,不要再用RHF->MP2 NO作为CASSCF初始轨道这种落后和不是很合理的技术了)



自动做多参考态计算的程序MOKIT

103

帖子

0

威望

607

eV
积分
710

Level 4 (黑子)

30#
发表于 Post on 2024-8-28 04:34:02 | 只看该作者 Only view this author
wzkchem5 发表于 2023-5-17 20:23
你要算什么性质?UV-Vis,还是反应性,还是磁性,还是其他的什么性质?
光给定体系,是无法唯一确定应该 ...

请问王老师,在这篇文章的体系中你为什么认为组氨酸的氮原子的孤对电子需要被包含在活性轨道内而同为配体的谷氨酸羧基却不需要呢?是因为相比组氨酸羧基是弱场配体吗?

另外,如果某条轨道在优化中总是容易转出活性空间,是不是意味着其占据数过于接近2,本身就没有必要放入活性空间中?我在一个体系中使活性空间包含了三个过渡金属的全部d轨道,发现某些轨道总是倾向于转出活性空间。在更换优化方式后最终得到了活性空间包含全部d轨道的结果,果然容易转出的d轨道占据数都十分接近2,要么1.9999要么干脆就是2.000。对比转出了活性空间的结果,发现不限制旋转时能量更低,大约1 kcal/mol左右,这是否意味着如果不考虑UV-Vis等性质,单从计算能量的角度,相同活性空间下包含全部的d轨道并不是最好的选择?

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 03:43 , Processed in 0.229191 second(s), 23 queries , Gzip On.

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