求助过渡金属体系电离能计算偏离实验值比较多且得不到dyson轨道

天真有邪

过渡金属氧化物体系 用newtonx 算了dyson轨道 任务正常结束 但是电离能计算偏离实验值比较多而且轨道系数是0 就还是得不到dyson轨道 找不到问题 请问各位老师们 有没有遇到过这种问题该如何解决？是不是DFT的方法不适用这个体系 需要改用多参考的方法计算dyson轨道？求助该如何计算

sobereva

帖子标题要写成完整的句子，别就随便摆三个关键词 过渡金属体系 多参考计算 dyson轨道。这次给你改了，以后注意！认真看http://bbs.keinsci.com/thread-9348-1-1.html

说清楚算Dyson轨道的目的。
算电离能误差大有的是可能原因，诸如泛函不合适、基组不合适、自旋多重度不对或波函数不稳定、构型/构象不对、和实验缺乏可比性，等等。完全不提体系和计算细节，以及理论和实验数据具体值，根本没人能确切回答。认真看下文
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）

天真有邪

sobereva 发表于 2025-4-17 16:06
帖子标题要写成完整的句子，别就随便摆三个关键词 过渡金属体系 多参考计算 dyson轨道。这次给你改了，以后 ...

不好意思sob老师谢谢您的提示，我一定改正。
对于我的计算体系是过渡金属氧化物体系PtO阴离子，在激光作用下失去电子成为中性分子。使用高斯B3LYP/aug-cc-pvtz(-pp)结构优化且波函数稳定检测过是稳定的 后续基于优化后的中性阴离子结构计算dyson轨道 使用newtonx计算，电荷自旋多重度设置应该没问题 是阴离子基态二重态开壳层，中性三重态。使用DFT方法，基组不变，以下附输入文件设置：
&dat
nact = 2
numat = 2       
npoints = 500
file_geom = geom
iprog = 4
file_nmodes = freq.log
anh_f = 1
temp = 13
ics_flg = y
chk_e = 0
prog = 6.5       
iseed = 0
nis_mol = 1
nfs_cat = 21
ip_type = 2
do_type = 2
eq_geom = y
lvprt = 1
/
文件正常结束但是第一电离能偏离实验值较大且第一个电离dyson轨道系数为0；
后续我还在采用qchem软件 eom-ip-ccsd方法计算过 暂未收敛
辛苦老师指正~

sobereva

天真有邪 发表于 2025-4-17 20:59
不好意思sob老师谢谢您的提示，我一定改正。
对于我的计算体系是过渡金属氧化物体系PtO阴离子，在激光作 ...

算PtO这么小的体系的第一IP，没有任何理由用DFT来算，相对来说太糙了
高精度计算可以用MRCI，黑箱方法里至少也该用个CCSD(T)或IP-EOM-CCSD（Q-Chem如果不成功可以尝试用ORCA或NWChem的）

我仍然不知道你算Dyson轨道的实际目的何在，说不定其实都可以不用算，有其它方法能讨论。诸如对电离过程中电子密度的改变感兴趣，直接用Multiwfn对前后两个态求密度差就能考察

天真有邪

sobereva 发表于 2025-4-18 05:39
算PtO这么小的体系的第一IP，没有任何理由用DFT来算，相对来说太糙了
高精度计算可以用MRCI，黑箱方法里 ...

实际上老师我们想看阴离子在激光作用下电离电子的轨道 所以想计算dyson轨道直观考察 电离能附加信息看一下结果是否准确。电荷密度差也有计算 主要是想直观的考察所以要看dyson轨道
另外本体系静态相关较强 需要使用多参考方法来算 所以想请教老师如何使用多参考方法计算dyson轨道。