计算化学公社

标题: 求助过渡金属体系电离能计算偏离实验值比较多且得不到dyson轨道 [打印本页]

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

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

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

作者
Author: 天真有邪 时间: 2025-4-17 20:59

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方法计算过暂未收敛
辛苦老师指正~

作者
Author: sobereva 时间: 2025-4-18 05:39

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

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

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

作者
Author: 天真有邪 时间: 2025-4-24 15:30

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

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

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