ORCA5.0如何通过设置OrbWin来控制电子激发行为

wangxubin

大家好，我想请教一下关于ORCA5.0中OrbWin参数的使用。
具体以计算乙烷S1与T1态的自旋轨道耦合相互作用为例，我的目的是想限制电子的激发仅在HOMO->LUMO之间，通过查阅手册，发现应使用OrbWin参数进行设置，但并未理解手册中的示例，希望可以有大佬解释一下，麻烦大家了。
以下是我初步制作的输入文件：
! B3LYP/G TZVP pal4
%tddft
orbwin = #此处不知该如何设置参数
mode stddft
nroots 5
dosoc true
tda false
end
* xyz 0 1
C          1.00610        0.06479       -0.07614
C          2.51817        0.06479       -0.07614
H          0.62183        0.70746       -0.87383
H          0.62183        0.43427        0.87927
H          0.62183       -0.94736       -0.23387
H          2.90244        1.07694        0.08158
H          2.90245       -0.30469       -1.03155
H          2.90245       -0.57788        0.72154
*

wzkchem5

首先确认一下你要做的计算到底是限制只有HOMO->LUMO的激发，还是限制只有以HOMO->LUMO的激发为主要成分的激发态。一般极少有研究会涉及“限制电子的激发仅在HOMO->LUMO之间”这种计算，因为计算结果不对应任何实验可观测量。所以一定要再三确认你确实要做这个，再讨论怎么做的问题

wangxubin

wzkchem5 发表于 2021-12-18 01:04
首先确认一下你要做的计算到底是限制只有HOMO->LUMO的激发，还是限制只有以HOMO->LUMO的激发为主要成分的激 ...

你好，我确定是“仅有HOMO->LUMO的激发”(与之前的工作对应)，但不知道如何设置输入文件。

wzkchem5

wangxubin 发表于 2021-12-18 08:11
你好，我确定是“仅有HOMO->LUMO的激发”(与之前的工作对应)，但不知道如何设置输入文件。

%tddft
OrbWin[0] = a0,b0,c0,d0
OrbWin[1] = a1,b1,c1,d1
end

比如这个表示允许alpha电子从第a0个轨道（含）到第b0个轨道（含）之中的任何一个占据轨道，跃迁到第c0个轨道（含）到第d0个轨道（含）之中的任何一个空轨道；与此同时，允许beta电子从第a1个轨道（含）到第b1个轨道（含）之中的任何一个占据轨道，跃迁到第c1个轨道（含）到第d1个轨道（含）之中的任何一个空轨道。如果b0或b1写成-1，代表HOMO；如果c0或c1写成-1，代表LUMO（此时程序自动查询HOMO和LUMO的实际编号，用实际编号替换这个-1）
如果是RKS计算，不能写OrbWin[1]这行，OrbWin[1]会自动设成和OrbWin[0]一致。另外注意轨道编号是从0开始的，不是从1开始的。
这个确实属于说明书没写清楚，而且说明书有的地方少写了那个等号，如果不写等号会报错。我这边改一下说明书，从orca 5.1.0开始应该就清楚了

wangxubin

wzkchem5 发表于 2021-12-19 00:39
%tddft
OrbWin[0] = a0,b0,c0,d0
OrbWin[1] = a1,b1,c1,d1

明白了，谢谢老师