计算化学公社

标题: 用orca计算开壳单线态和闭壳单线态的单点能量 [打印本页]

作者
Author: 米宝舞 时间: 2018-8-5 21:36
标题: 用orca计算开壳单线态和闭壳单线态的单点能量
本帖最后由米宝舞于 2018-8-5 22:12 编辑

请问各位老师，用orca计算开壳单线态和闭壳单线态的单点能量时，输入文件可以像如下这样写吗：
开壳单线态：
#
# 1a_21u_b3lypoptPWPB95sp
#

!UKS RI-PWPB95 D3BJ def2-tzvp def2/JK def2-tzvp/C RIJK TIGHTSCF Grid5 FINALGrid6
%PAL NPROCS 6
end%scf
Flipspin 0
FinalMs 0.0
%maxcore 1500
%cpcm smd true
solvent "acetonitrile"
end

* xyz 2 1
C                2.29090200 1.77490300 -1.48680300
C                3.75512200 2.02545800 -1.82179200
......
*
闭壳单线态：
#
# 1a_21r_b3lypoptPWPB95sp
#

!RKS RI-PWPB95 D3BJ def2-tzvp def2/JK def2-tzvp/C RIJK TIGHTSCF Grid5 FINALGrid6
%PAL NPROCS 6
end
%scf
Flipspin 0
FinalMs 0.0
%maxcore 1500
%cpcm smd true
solvent "acetonitrile"
end

* xyz 2 1
C                2.29090200 1.77490300 -1.48680300
C                3.75512200 2.02545800 -1.82179200
......
*
我用orca计算三重态的输入文件为：
#
# 1a_b3lypoptPWPB95sp
#

!UKS RI-PWPB95 D3BJ def2-tzvp def2/JK def2-tzvp/C RIJK TIGHTSCF Grid5 FINALGrid6
%PAL NPROCS 6
end
%maxcore 1500
%cpcm smd true
solvent "acetonitrile"
end

* xyz 2 3
C                2.29090200 1.77490300 -1.48680300
C                3.75512200 2.02545800 -1.82179200
......
*

作者
Author: sobereva 时间: 2018-8-6 06:00
算开壳层单重态spinflip后面显然不能是0，否则达不到令初猜波函数处于对称破缺态的目的。再仔细理解下ORCA手册里的例子

作者
Author: 米宝舞 时间: 2018-8-6 08:26

sobereva 发表于 2018-8-6 06:00
算开壳层单重态spinflip后面显然不能是0，否则达不到令初猜波函数处于对称破缺态的目的。再仔细理解下ORCA ...

收到谢谢sob老师的解答我再回去看一下orca的手册

作者
Author: 米宝舞 时间: 2018-8-9 21:55

sobereva 发表于 2018-8-6 06:00
算开壳层单重态spinflip后面显然不能是0，否则达不到令初猜波函数处于对称破缺态的目的。再仔细理解下ORCA ...

老师，我根据看的手册改了一下，您看开壳单线态和闭壳单线态这样输入可以吗
开壳单线态：
!UKS RI-PWPB95 D3BJ def2-tzvp def2/JK def2-tzvp/C RIJK TIGHTSCF Grid5 FINALGrid6
%PAL NPROCS 6
end
%scf
Flipspin 1
FinalMs 0.0
%maxcore 1500
%cpcm smd true
solvent "acetonitrile"
end

闭壳单线态：
!RKS RI-PWPB95 D3BJ def2-tzvp def2/JK def2-tzvp/C RIJK TIGHTSCF Grid5 FINALGrid6
%PAL NPROCS 6
end
%scf
Flipspin 0
FinalMs 0.0
%maxcore 1500
%cpcm smd true
solvent "acetonitrile"
end

作者
Author: sobereva 时间: 2018-8-10 05:03

米宝舞发表于 2018-8-9 21:55
老师，我根据看的手册改了一下，您看开壳单线态和闭壳单线态这样输入可以吗
开壳单线态：
!UKS RI-PWPB ...

flipspin后面设的是自旋要翻转的原子序号，该设什么看你具体体系
不确信的话就算算试试，之后用Multiwfn绘制自旋密度分布，看看是否已经成了双自由基状态

作者
Author: 米宝舞 时间: 2018-8-10 07:13

sobereva 发表于 2018-8-10 05:03
flipspin后面设的是自旋要翻转的原子序号，该设什么看你具体体系
不确信的话就算算试试，之后用Multiwfn ...

收到谢谢老师的解答

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