orca能否直接计算黄昆因子

好兄弟

尊敬的各位老师：用gaussian计算激发态t2能量总比t1低，现在用orca计算出的t1 t2是正常的。现在我想用orca的计算结果分析HR，.out文件里输出频率和displacement项，推测这是可行的，请问有什么方法可以用orca直接计算黄昆因子吗？

cokie

可以尝试FCclasses，载入.hess生成的state，做个emi即可得到HuangRhys的输出。

wzkchem5

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

陈AG

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

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

wal

陈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

陈AG

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

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

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

wal

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

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

陈AG

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也没有

wal

陈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里的

陈AG

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

wal

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

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

陈AG

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

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

wal

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

对的

sxkl

大佬有没有完整输入文件参考一下呀？想计算TADF材料的HR因子

wal

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

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