计算化学公社

标题: 在molpro中，用多参考方法计算键解离过程势能曲线遇到问题 [打印本页]

作者
Author: Sswzar 时间: 2021-7-15 19:54
标题: 在molpro中，用多参考方法计算键解离过程势能曲线遇到问题
计算步骤就是先用gaussian做键解离的柔性扫描，把得到的曲线上的每一个点的结构拿到molpro里面用CASPT2(2e,2o)/aug-cc-pVDZ//CASSCF(2e,2o)/aug-cc-pVDZ做单点计算，但是计算结果好像还没有Gaussian里面的DFT合理，是不是我的操作步骤出了问题呢？求老师们帮忙看看吧

作者
Author: wzkchem5 时间: 2021-7-15 22:50
检查一下活性空间有没有选对，即使初始轨道是对的，CASSCF也有可能把活性轨道优化到其他地方去，尤其当你扫描的键长比较短的时候

作者
Author: Sswzar 时间: 2021-7-16 11:00

wzkchem5 发表于 2021-7-15 22:50
检查一下活性空间有没有选对，即使初始轨道是对的，CASSCF也有可能把活性轨道优化到其他地方去，尤其当你扫 ...

感谢您的回复，对于这个分子由于只涉及到碳碳键的断裂（一个四元环的开环），因此我直接选择了2e,2o，应该没什么问题吧
那么还想请问老师，怎么在molpro里面指定活化的轨道呢
这是我的输入
basis=AVDZ
{hf;wf,76,1,0;}
{casscf;occ,39;closed,37;}
{rs2,shift=0.2;occ,39;closed,37;orbital;maxit,500}
因为是新手，所以好多操作都不太明白，感谢老师的指导

作者
Author: wzkchem5 时间: 2021-7-16 14:54

Sswzar 发表于 2021-7-16 04:00
感谢您的回复，对于这个分子由于只涉及到碳碳键的断裂（一个四元环的开环），因此我直接选择了2e,2o，应 ...

把收敛的活性轨道的图作出来，看看是不是一个sigma轨道一个sigma*轨道，如果不是，通过旋转轨道让活性轨道是一个sigma轨道一个sigma*轨道，然后重新做CASSCF

作者
Author: Sswzar 时间: 2021-7-16 15:54
老师您好，我检查了轨道的图，活化的就是最后两个轨道吧，我觉得分别就是成键轨道和反键轨道，好像没有问题呢，这让我很疑惑

作者
Author: agoni 时间: 2024-8-30 21:36

Sswzar 发表于 2021-7-16 15:54
老师您好，我检查了轨道的图，活化的就是最后两个轨道吧，我觉得分别就是成键轨道和反键轨道，好像没有问题 ...

您好，请问这个问题您最后是怎么解决的呢？我也遇到了类似的问题

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