ORCA S1 T1 分子轨道计算问题

是aweia

1. 各位老师好 我想计算分子在S1 T1态下的分子轨道图 我计算T1态时用如下输入
   
2. %tddft
   
3. nroots 5
   
4. IROOTMULT TRIPLET
   
5. * XYZFILE 0 1 an.xyz
   
6. 
   
7. S1态计算时把
   
8. IROOTMULT TRIPLET   改为 IROOT 1

复制代码

将产生的gbw文件转为molden文件后用multiwfn打开绘图发现 S1和T1态的分子轨道图一样，
1、我猜测是否是我的S1 T1计算的到的gbw文件中储存的都是”参考态轨道“，请问我的输入情况下，参考态就是基态么？那么如何指定参考态呢？如果我想把S1、T1态下的分子轨道保存进gbw文件而不是参考态轨道，我应该加如什么参数么？




2、另外我希望得到alpha和beta各自的分子轨道图，因为我使用的是闭壳层，是否只有开壳层方法才能得到分离的alpha和beta轨道
我直接将

• * XYZFILE 0 1 an.xyz 改为
• * XYZFILE 0 3 an.xyz
  

并删除TDDFT部分，确实可以得到alpha和beta分离的轨道那么请问如果对于S1态计算 我该如何写输入文件？

ionexchangeC

计算T1的分子轨道直接把自旋多重度改成3就行
计算S1的分子轨道比较麻烦，deltaSCF可以做到但必须体系适合（相关：http://sobereva.com/265 乱谈激发态的计算方法）（ORCA大概做不到）

是aweia

ionexchangeC 发表于 2022-4-25 12:48
计算T1的分子轨道直接把自旋多重度改成3就行
计算S1的分子轨道比较麻烦，deltaSCF可以做到但必须体系适合 ...

谢谢您的回复 那如果是Gaussian呢 可以做到么

wzkchem5

是aweia 发表于 2022-4-25 14:54
谢谢您的回复 那如果是Gaussian呢 可以做到么

能做到的概率应该不会比orca大。
注意对于有的激发态（尤其是比较高的激发态），deltaSCF从数学上就不可能收敛到你要的激发态。对于deltaSCF能收敛到的激发态，又分为两种情况，一种是靠普通SCF程序，只要准备合适的初猜就能收敛的，这种高斯、orca都能做；另一种是要靠MOM等方法才能收敛的，这种高斯、orca都做不到，需要用BDF、Q-Chem等软件。
关于MOM以及相关的一些方法，参见doi: 10.1021/ct500066t以及这篇文章引用的文献

是aweia

wzkchem5 发表于 2022-4-25 22:12
能做到的概率应该不会比orca大。
注意对于有的激发态（尤其是比较高的激发态），deltaSCF从数学上就不可 ...

好的 非常感谢您的回复和指导

zjxitcc

“S1、T1态下的分子轨道”是模糊表述，你没说要哪种类型的轨道，那TDDFT的S1态的自然轨道是很好获得的（自然轨道没有轨道能量，有轨道占据数；对应地，没有HOMO、LUMO，有HONO、LUNO），高斯和ORCA都行。注意：
（1）自旋多重度写3，用基态UHF/UDFT做计算，这样是算T1的经济性做法（也是最常用的），这得到的是三重态UHF/UDFT正则轨道，有对应的轨道能量。跟TDDFT计算是两码事。
（2）自旋多重度写1，外加TDDFT计算S1或T1，你可以再加关键词获得S1或T1态的自然轨道。但不存在“S1或T1态的正则轨道”。

sobereva

是aweia

zjxitcc 发表于 2022-4-25 22:20
“S1、T1态下的分子轨道”是模糊表述，你没说要哪种类型的轨道，那TDDFT的S1态的自然轨道是很好获得的（自 ...

谢谢您的指导  对于您说的第二个情况    我查询了一下手册 还是搞不清楚我应该使用哪个关键词来获得S1、T1的自然轨道  想请教您一下是哪个关键词 非常感谢