Multiwfn对ORCA6.x的振子强度数据的读取可能存在bug

ABetaCarw

版本是windows Version 3.8(dev), update date: 2025-Jun-28
操作系统 Windows 10 专业版 22H2 19045.5854

在对6.0和6.1的ORCA SOC计算输出文件进行
18-15 Print major MO transitions in all excited states
分析单重态激发信息
分析中发现：
其关于激发能和波长可以正常读取，但是振子强度是错误的，比如水分子计算数据如下


分析结果是
#   1   8.6220 eV    143.80 nm   f=  0.00000   Spin multiplicity= 1:
   H -> L 98.8%
#   2  10.5210 eV    117.84 nm   f=  0.04603   Spin multiplicity= 1:
   H -> L+1 96.6%
#   3  10.8310 eV    114.47 nm   f=  0.00000   Spin multiplicity= 1:
   H-1 -> L 98.5%
#   4  12.8620 eV     96.40 nm   f=  0.00000   Spin multiplicity= 1:
   H-1 -> L+1 97.9%
#   5  13.8640 eV     89.43 nm   f=  0.00000   Spin multiplicity= 1:
   H-2 -> L 98.9%
#   6  16.3040 eV     76.05 nm   f=  0.12666   Spin multiplicity= 1:
   H-2 -> L+1 86.9%, H -> L+4 10.5%
#   7  16.5900 eV     74.73 nm   f=  0.00000   Spin multiplicity= 1:
   H -> L+2 98.1%
#   8  17.4640 eV     70.99 nm   f=  0.06499   Spin multiplicity= 1:
   H -> L+3 79.8%, H -> L+5 19.7%
#   9  18.7010 eV     66.30 nm   f=  0.00000   Spin multiplicity= 1:
   H -> L+4 59.5%, H-1 -> L+3 28.1%, H-1 -> L+5 8.3%
#  10  18.9740 eV     65.34 nm   f=  0.25420   Spin multiplicity= 1:
   H -> L+5 78.8%, H -> L+3 20.0%
#  11  19.3220 eV     64.17 nm   f=  0.00000   Spin multiplicity= 1:
   H-1 -> L+2 97.2%
#  12  19.8300 eV     62.52 nm   f=  0.10314   Spin multiplicity= 1:
   H-1 -> L+4 99.2%
#  13  20.7730 eV     59.69 nm   f=  0.00000   Spin multiplicity= 1:
   H-1 -> L+5 49.2%, H-1 -> L+3 45.2%
#  14  22.0030 eV     56.35 nm   f=  0.00000   Spin multiplicity= 1:
   H -> L+6 94.5%
#  15  22.4970 eV     55.11 nm   f=  0.00000   Spin multiplicity= 1:
   H-2 -> L+2 68.7%, H-1 -> L+3 9.7%, H-1 -> L+5 7.7%, H -> L+4 6.4%
#  16  23.0140 eV     53.87 nm   f=  0.00000   Spin multiplicity= 1:
   H -> L+7 99.3%
#  17  23.0520 eV     53.78 nm   f=  0.01424   Spin multiplicity= 1:
   H-2 -> L+3 76.0%, H-2 -> L+5 19.4%
#  18  23.2000 eV     53.44 nm   f=  0.00000   Spin multiplicity= 1:
   H-2 -> L+4 96.4%
#  19  23.2290 eV     53.37 nm   f=  0.00000   Spin multiplicity= 1:
   H-2 -> L+2 28.6%, H-1 -> L+5 27.4%, H -> L+4 12.6%, H-2 -> L+6 12.3%, H-1 -> L+3 11.5%
#  20  24.1880 eV     51.26 nm   f=  0.00000   Spin multiplicity= 1:
   H-1 -> L+6 82.4%, H-2 -> L+3 10.8%, H-2 -> L+5 5.1%




water.out (498.09 KB, 下载次数 Times of downloads: 2)

2025-7-8 02:00 上传 Uploaded

点击下载
Click to download

Oxygen

这个振子强度不是对的吗，问题在哪？

wal

Oxygen 发表于 2025-7-8 11:17
这个振子强度不是对的吗，问题在哪？

看起来波长列表正确筛掉了triplet的，但振子强度没有筛掉

sobereva

ORCA 6改变了输出格式所致，完全不是Multiwfn的bug，是ORCA开发者太任性
今日更新的Multiwfn兼容了ORCA 6

wzkchem5

sobereva 发表于 2025-7-10 11:25
ORCA 6改变了输出格式所致，完全不是Multiwfn的bug，是ORCA开发者太任性
今日更新的Multiwfn兼容了ORCA 6
...

ORCA团队比较建议的方法是通过property file或json文件读数据，可以保证将来的ORCA版本里格式不会变化

sobereva

wzkchem5 发表于 2025-7-10 19:54
ORCA团队比较建议的方法是通过property file或json文件读数据，可以保证将来的ORCA版本里格式不会变化

顺便咨询下，TDDFT的激发和去激发组态系数，有没有办法从property或json之类文件里，或者把某二进制文件以什么方式转成文本形式，然后读取到？ORCA TDDFT输出文件里不直接给出去激发组态系数

wzkchem5

sobereva 发表于 2025-7-10 22:38
顺便咨询下，TDDFT的激发和去激发组态系数，有没有办法从property或json之类文件里，或者把某二进制文件 ...

这个用json是可以的，json里有X+Y和X-Y的amplitudes。
用ORCA 6.1.0新出的python interface（https://www.faccts.de/docs/opi/1.0/docs/）也行，但缺点是不支持6.1.0以前的版本。
原则上*.cis文件（二进制文件）里也有amplitude信息，但orca官方还没公开*.cis文件的格式，也没有提供从*.cis文件转换为文本文件的方法

sobereva

wzkchem5 发表于 2025-7-11 09:30
这个用json是可以的，json里有X+Y和X-Y的amplitudes。
用ORCA 6.1.0新出的python interface（https://ww ...

谢谢，我有空研究下

sobereva

wzkchem5 发表于 2025-7-11 09:30
这个用json是可以的，json里有X+Y和X-Y的amplitudes。
用ORCA 6.1.0新出的python interface（https://ww ...

可惜当前版本有bug：https://orcaforum.kofo.mpg.de/viewtopic.php?f=8&t=13648，又不打算折腾python interface。Multiwfn支持ORCA的严格的TDDFT激发态分析只能等next release了，本来读取代码都写好了