|
|
笔者在用ORCA计算一个过渡态为自旋极化单重态的反应的IRC,
过渡态的<S**2>在1左右,波函数稳定,
笔者摹仿这篇帖子过渡态是自旋极化单重态时找TS和跑IRC流程分享(以坛友的实际问题为例) - 量子化学 (Quantum Chemistry) - 计算化学公社在Gaussian中的操作,
用ORCA计算过渡态为自旋极化单重态的反应的IRC,读取了之前优化自旋极化单重态的过渡态的波函数,
请问在跑IRC时需要特意两个方向分开来跑吗?
因为笔者想尝试用过渡态与反应物/产物之间的几个结构进行波函数分析,查看反应过程中自旋密度的变化,
由于产物与反应物都是闭壳层单重态,笔者担心算完一个方向后,算另一个方向就会以计算闭壳层为前提。
关键词如下:
! UKS B3LYP D4 def2-TZVP def2/J RIJCOSX IRC tightSCF SlowConv noautostart miniprint moread
%base "...IRC"
%moinp "...TS.gbw"
%maxcore 20000
%pal nprocs 20 end
%irc
MaxIter 50
Direction forward
InitHess read
Hess_Filename "...TS.hess"
Init_Displ DE
Scale_Init_Displ 0.05
DE_Init_Displ 2.0
end
|
|