计算化学公社

标题: Multiwfn对ORCA6.x的振子强度数据的读取可能存在bug [打印本页]

作者
Author: ABetaCarw 时间: 2025-7-8 02:01
标题: Multiwfn对ORCA6.x的振子强度数据的读取可能存在bug
版本是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
分析单重态激发信息
分析中发现：
其关于激发能和波长可以正常读取，但是振子强度是错误的，比如水分子计算数据如下
(, 下载次数 Times of downloads: 259)

分析结果是
# 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%

(, 下载次数 Times of downloads: 2)

作者
Author: Oxygen 时间: 2025-7-8 11:17
这个振子强度不是对的吗，问题在哪？

作者
Author: wal 时间: 2025-7-8 13:03

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

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

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

作者
Author: wzkchem5 时间: 2025-7-10 19:54

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

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

作者
Author: sobereva 时间: 2025-7-10 22:38

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

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

作者
Author: wzkchem5 时间: 2025-7-11 09:30

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文件转换为文本文件的方法

作者
Author: sobereva 时间: 2025-7-11 22:32

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

谢谢，我有空研究下

作者
Author: sobereva 时间: 2025-7-17 02:08

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了，本来读取代码都写好了

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