请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 385|回复: 5

[ORCA] 用orca计算开壳单线态和闭壳单线态的单点能量

[复制链接]

11

帖子

0

威望

131

eV
积分
142

Level 2 能力者

发表于 2018-8-5 21:36:41 | 显示全部楼层 |阅读模式
本帖最后由 米宝舞 于 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
......
*




要计算单点的中间体.zip

1.25 MB, 下载次数: 1

1万

帖子

25

威望

1万

eV
积分
34842

管理员

公社社长

发表于 2018-8-6 06:00:46 | 显示全部楼层
算开壳层单重态spinflip后面显然不能是0,否则达不到令初猜波函数处于对称破缺态的目的。再仔细理解下ORCA手册里的例子
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

11

帖子

0

威望

131

eV
积分
142

Level 2 能力者

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

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

11

帖子

0

威望

131

eV
积分
142

Level 2 能力者

 楼主| 发表于 2018-8-9 21:55:39 | 显示全部楼层
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

1万

帖子

25

威望

1万

eV
积分
34842

管理员

公社社长

发表于 2018-8-10 05:03:14 | 显示全部楼层
米宝舞 发表于 2018-8-9 21:55
老师,我根据看的手册改了一下,您看开壳单线态和闭壳单线态这样输入可以吗
开壳单线态:
!UKS RI-PWPB ...

flipspin后面设的是自旋要翻转的原子序号,该设什么看你具体体系
不确信的话就算算试试,之后用Multiwfn绘制自旋密度分布,看看是否已经成了双自由基状态
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

11

帖子

0

威望

131

eV
积分
142

Level 2 能力者

 楼主| 发表于 2018-8-10 07:13:50 | 显示全部楼层
sobereva 发表于 2018-8-10 05:03
flipspin后面设的是自旋要翻转的原子序号,该设什么看你具体体系
不确信的话就算算试试,之后用Multiwfn ...

收到 谢谢老师的解答
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-10-18 05:41 , Processed in 0.196740 second(s), 27 queries .

快速回复 返回顶部 返回列表