计算化学公社

标题: 关于自由基复合反应体系进行统一活性空间 [打印本页]

作者
Author: CCLI 时间: 2025-7-4 13:10
标题: 关于自由基复合反应体系进行统一活性空间
各位老师好！我在使用多参考态方法做自由基复合反应体系时，使用了MOKIT分别对TS、反应物和产物进行casscf/cc-pVQZ的计算，目前还没有进行更为精准的基组计算(def2-TZVP)。没指定活性空间时，MOKIT分别对TS、反应物和反应物选择了（2,2）、(2,2)和（10,10）作为活性空间，由于本人在后hf计算经验严重不足，以下几点是我想问的：
1.轨道顺序选择，活性空间统一需要以过渡态为标准，我查看轨道发现TS与反应物相同，但两者等值面是相反的，即TS的HONO对应反应物的LUNO，而产物完全对应不上，对于前者需要交换轨道吗？对于后者需要继续搜寻轨道吗？(同TS和反应物相同的轨道)
2.反应物和TS用MOKIT计算没有问题，但产物的结果出现了一些提示和部分文件夹没有自动删除，并且前者我都是用xxx_gvb4_s.fch查看轨道的，而后者是通过xxx_rhf_gvb83_CASSCF_NO.fch查看轨道(因为没有生成xxx_gvb4_s.fch)，对于没有生成xxx_gvb4_s.fch，是不是没有形成明显sigma键，pi键，因此才不产生？
3.如果是统一到(2,2)是不是不可信？我听说，(2,2)是10多年的计算水平了，现在起码是(6,6),如果我选择了(2,2),审稿人问我，我该如何去回答？
4.关于后续的NEVPT2的计算，我知道有三种方法：DLPNO-NEVPT2、FIC-NEVPT2和SC_NEVPT2，我体系的原子数为66，用orca在wb97x-d3/def2-TZVP水平下计算，ts和反应物为开壳层单重态，产物为闭壳层单重态，根据上面的信息，我应该选择哪种方法去计算，我们这边最好的机子有两台：36核物理核心，125GB内存以及16核物理核心，250GB内存。
以下是我处理好的的轨道图和异常的结果信息，以及我的计算文件通过网盘分享的文件：CASSCF-E
链接: https://pan.baidu.com/s/160IbAIk3in_iySX8wiUROw 提取码: ciri

作者
Author: CCLI 时间: 2025-7-4 13:13
这里漏了产物的轨道排序图

作者
Author: hebrewsnabla 时间: 2025-7-4 20:54
2. 肯定会生成*_gvb*_s.fch 的。只是不一定是gvb4。
3. 没有这种说法。该用多大就用多大，并不是越大越好，并不是越大代表计算水平越高。

粗看一下的话，似乎(2,2)是足够的。产物可以去gvb*_s.fch里面找到对应过渡态的那一对轨道。

作者
Author: CCLI 时间: 2025-7-4 21:19
本帖最后由 CCLI 于 2025-7-4 21:22 编辑

hebrewsnabla 发表于 2025-7-4 20:54
2. 肯定会生成*_gvb*_s.fch 的。只是不一定是gvb4。
3. 没有这种说法。该用多大就用多大，并不是越大越好 ...

确实，都有*_gvb*_s.fch，请问这些前面的命名是不是与开壳层、闭壳层体系有关，见下图红框
[attach]115702[/attach] (, 下载次数 Times of downloads: 24)

作者
Author: hebrewsnabla 时间: 2025-7-4 21:27

CCLI 发表于 2025-7-4 21:19
确实，都有*_gvb*_s.fch，请问这些前面的命名是不是与开壳层、闭壳层体系有关，见下图红框

是的。

作者
Author: CCLI 时间: 2025-7-4 22:32

hebrewsnabla 发表于 2025-7-4 21:27
是的。

好的，谢谢老师

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