计算化学公社

标题: ORCA S1 T1 分子轨道计算问题 [打印本页]

作者
Author: 是aweia 时间: 2022-4-25 10:39
标题: ORCA S1 T1 分子轨道计算问题
本帖最后由是aweia 于 2022-4-25 10:45 编辑

各位老师好我想计算分子在S1 T1态下的分子轨道图我计算T1态时用如下输入
%tddft
nroots 5
IROOTMULT TRIPLET
* XYZFILE 0 1 an.xyz
S1态计算时把
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态计算我该如何写输入文件？

作者
Author: ionexchangeC 时间: 2022-4-25 12:48
本帖最后由 ionexchangeC 于 2022-4-25 12:50 编辑

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

作者
Author: 是aweia 时间: 2022-4-25 21:54

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

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

作者
Author: wzkchem5 时间: 2022-4-25 22:12

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

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

作者
Author: 是aweia 时间: 2022-4-25 22:16

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

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

作者
Author: zjxitcc 时间: 2022-4-25 22:20
本帖最后由 zjxitcc 于 2022-4-25 22:22 编辑

“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态的正则轨道”。

作者
Author: sobereva 时间: 2022-4-26 08:00
(, 下载次数 Times of downloads: 14)

作者
Author: 是aweia 时间: 2022-4-26 14:38

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

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

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