为什么使用Multiwfn+VMD得到的总的跃迁偶极矩矢量与Gaussian log文件输出的结果不一致

chlei

Sob老师您好，我在使用Multiwfn+VMD绘制体系片段贡献的跃迁电偶极矩矢量（http://sobereva.com/396）时，通过在VMD界面输入“dip all”得到体系某个激发态总的电跃迁偶极矩矢量，但是发现一点问题：VMD中输出的结果与Gaussian计算log文件中“Ground to excited state transition electric dipole moments (Au):”相应激发态的电跃迁偶极矩矢量不一致，二者存在“很小”的差异。想请问：1）这种差异的来源是什么，是因为Gaussian计算时我所使用的关键词IOp(9/40=4)吗？2）应以哪个结果为准？非常感谢！

sobereva

不写IOp(9/40=4)差异只会更大，搞清楚写这个的目的是什么
想要差别更小写IOp(9/40=5)
直接拿具体数值说事，别就光说个“很小”，谁知道是多小

chlei

sobereva 发表于 2024-10-28 12:06
不写IOp(9/40=4)差异只会更大，搞清楚写这个的目的是什么
想要差别更小写IOp(9/40=5)
直接拿具体数值说事 ...

Sob老师，您好！谢谢回复。我的这个体系S0-S4激发log文件输出的电跃迁偶极矩矢量(x,y,z)分别为(0.1223, -0.0103, 1.1132)；而基于Multiwfn+VMD输出的结果为(0.12167099999999986, -0.014123000000000085, 1.1779890000000002)。

sobereva

chlei 发表于 2024-10-28 13:46
Sob老师，您好！谢谢回复。我的这个体系S0-S4激发log文件输出的电跃迁偶极矩矢量(x,y,z)分别为(0.1223, - ...

差异来源，组态系数输出的仍不够多（想进一步改进吻合程度用IOp(9/40=5)），或者其它你没提到的可能因素，诸如结构有轻微差异等等。没具体文件没法说