计算化学公社

标题: orca能否直接计算黄昆因子 [打印本页]

作者
Author: 好兄弟 时间: 2024-6-25 11:33
标题: orca能否直接计算黄昆因子
尊敬的各位老师：用gaussian计算激发态t2能量总比t1低，现在用orca计算出的t1 t2是正常的。现在我想用orca的计算结果分析HR，.out文件里输出频率和displacement项，推测这是可行的，请问有什么方法可以用orca直接计算黄昆因子吗？

作者
Author: cokie 时间: 2024-6-25 15:52
可以尝试FCclasses，载入.hess生成的state，做个emi即可得到HuangRhys的输出。

作者
Author: wzkchem5 时间: 2024-6-25 16:31
把ESD模块的printlevel设成3（%esd printlevel 3 end）就可以打印Huang-Rhys因子，以及重组能的分解结果。这个功能手册里没写，下一版会加上

作者
Author: 陈AG 时间: 2025-8-28 09:12

wzkchem5 发表于 2024-6-25 16:31
把ESD模块的printlevel设成3（%esd printlevel 3 end）就可以打印Huang-Rhys因子，以及重组能的分解结果。 ...

老师，请问打印的Huang-Rhys因子是搜索out文件的什么关键词呢，就是Huang-Rhys吗

作者
Author: wal 时间: 2025-8-28 11:28

陈AG 发表于 2025-8-28 09:12
老师，请问打印的Huang-Rhys因子是搜索out文件的什么关键词呢，就是Huang-Rhys吗

我翻了翻我以前的输出应该是这样的

The Huang-Rhys factor per mode is:

mode HR-factor Freq
---- --------- ----

   1 0.116 18.06
   2 0.083 20.47
   3 0.121 46.69
   4 0.014 49.16

作者
Author: 陈AG 时间: 2025-8-28 11:35

wal 发表于 2025-8-28 11:28
我翻了翻我以前的输出应该是这样的

大佬，你这是有算磷光光谱吗，我现在快一个月了，还在Computing transition dipole derivatives             ...

            <<Displacing mode 1 (+)>>
            <<Displacing mode 2 (+)>>这些地方，是算完这个才会有HR因子吗，如果你结构不保密的话，可不可以发个out文件我看看，找找HR因子是在哪个位置

作者
Author: wal 时间: 2025-8-28 11:41

陈AG 发表于 2025-8-28 11:35
大佬，你这是有算磷光光谱吗，我现在快一个月了，还在Computing transition dipole derivatives ...

没算磷光，我这个是IC
按理说这个之前就打印了

作者
Author: 陈AG 时间: 2025-8-28 11:43

wal 发表于 2025-8-28 11:41
没算磷光，我这个是IC
按理说这个之前就打印了

! PBE0 DKH2 DKH-def2-TZVP(-f) SARC/J RIJCOSX tightSCF CPCM(toluene) ESD(PHOSP) miniprint
%maxcore  10000
%pal nprocs 10 end
%basis
NewGTO Pt "SARC-DKH-TZVP" end
end
%scf
ConvForced true
end
%TDDFT
NROOTS 5
DOSOC  TRUE
TDA FALSE
IROOT 1
printlevel 3

END
%ESD
GSHESSIAN "/data/3.AG-Pt-B/AG-Pt-B013/orca/8-27/S0-OPT-FREQ.hess"
TSHESSIAN "/data/3.AG-Pt-B/AG-Pt-B013/orca/8-27/T1-OPT-FREQ.hess"
DOHT TRUE
DELE 22120.51631
COORDSYS CARTESIAN
LINES    VOIGT
LINEW    75
INLINEW 200
END
* XYZFILE 0 1 /data/3.AG-Pt-B/AG-Pt-B013/orca/8-27/S0-OPT-FREQ.xyz，输入是这个，我搜HR没有，huang也没有

作者
Author: wal 时间: 2025-8-28 11:53

陈AG 发表于 2025-8-28 11:43
! PBE0 DKH2 DKH-def2-TZVP(-f) SARC/J RIJCOSX tightSCF CPCM(toluene) ESD(PHOSP) miniprint
%maxcor ...

你可以换个小点的模型试一下另外我的PRINTLEVEL是4 而且是写在ESD而非TDDFT里的

作者
Author: 陈AG 时间: 2025-8-29 11:51

wal 发表于 2025-8-28 11:53
你可以换个小点的模型试一下另外我的PRINTLEVEL是4 而且是写在ESD而非TDDFT里的

大佬，你知道为什么频率有299个mode，但是HR和重组能只有294个吗
! PBE0 DKH2 DKH-def2-TZVP(-f) SARC/J RIJCOSX tightSCF CPCM(toluene) ESD(PHOSP) miniprint
%maxcore  5000
%pal nprocs 20 end
%basis
NewGTO Pt "SARC-DKH-TZVP" end
end
%scf
ConvForced true
end
%TDDFT
NROOTS 5
DOSOC  TRUE
TDA FALSE
IROOT 1
END
%ESD
GSHESSIAN "/data/3.AG-Pt-B/AG-Pt-B013/orca/8-27/S0-OPT-FREQ.hess"
TSHESSIAN "/data/3.AG-Pt-B/AG-Pt-B013/orca/8-27/T1-OPT-FREQ.hess"
DOHT TRUE
DELE 22120.51631
COORDSYS CARTESIAN
LINES    VOIGT
LINEW    75
INLINEW 200
printlevel 4
END
* XYZFILE 0 1 /data/3.AG-Pt-B/AG-Pt-B013/orca/8-27/S0-OPT-FREQ.xyz

作者
Author: wal 时间: 2025-8-29 11:56

陈AG 发表于 2025-8-29 11:51
大佬，你知道为什么频率有299个mode，但是HR和重组能只有294个吗
! PBE0 DKH2 DKH-def2-TZVP(-f) SARC/J ...

ORCA的频率输出应该是包含平动/转动的，这些不参与分解重组能计算

作者
Author: 陈AG 时间: 2025-8-29 12:14

wal 发表于 2025-8-29 11:56
ORCA的频率输出应该是包含平动/转动的，这些不参与分解重组能计算

那是不是前面那几个频率不参与呢，因为前面那几个频率为0，任何后面那294个才跟HR和重组能一一对应

作者
Author: wal 时间: 2025-8-29 12:17

陈AG 发表于 2025-8-29 12:14
那是不是前面那几个频率不参与呢，因为前面那几个频率为0，任何后面那294个才跟HR和重组能一一对应

对的

作者
Author: sxkl 时间: 2025-9-29 17:35
大佬有没有完整输入文件参考一下呀？想计算TADF材料的HR因子

作者
Author: wal 时间: 2025-9-29 17:39

sxkl 发表于 2025-9-29 17:35
大佬有没有完整输入文件参考一下呀？想计算TADF材料的HR因子

问之前不翻翻帖子么楼上有人贴了输入可自行修改

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