求助分子电离前后结构变化的问题

pet

wzkchem5 发表于 2021-1-11 17:14
电子密度差分只需要电子密度，原则上不需要任何轨道。如果你用来做密度差分的软件没有办法直接读取你的量 ...

十分感谢老师的耐心回复，我优化用的M06-2X，单点能用的CCSD(T)/cc-pVTZ，如果单点能没有加density pop=NO（或NOAB)，他是直接算的HF轨道吗(如果是的话我上面就写错了，右边的不是自然轨道是HF轨道) ， 这个结果也可以直接拿来做电子密度差吗？

而且我CCSD(T)/cc-pVTZ加density pop=NOAB 总是报错， 显示：Post-SCF densities or gradients only with Real MP2, MP3, MP4SDQ, CI, CCD, and QCI.   
老师您的意思是加了NBO关键词的就不能用还是说可以不用呢？因为我发现加了跟没加的算出来不一样，加了的好像更好用来解释我的问题（垂直电离过程C2-C4比C2-C3掉的电子更多）。

snljty

pet 发表于 2021-1-11 17:45
十分感谢老师的耐心回复，我优化用的M06-2X，单点能用的CCSD(T)/cc-pVTZ，如果单点能没有加density pop=N ...

Gaussian没有CCSD(T)的一阶解析导数，但是计算后HF密度必须计算一阶导数。非要用CCSD(T)级别，换PSI4。看卢老师博文《在Multiwfn中分析比CCSD更高级别波函数的方法》 http://bbs.keinsci.com/thread-7396-1-1.html

wzkchem5

pet 发表于 2021-1-11 17:45
十分感谢老师的耐心回复，我优化用的M06-2X，单点能用的CCSD(T)/cc-pVTZ，如果单点能没有加density pop=N ...

不加pop=NO就是HF轨道。HF轨道可以用来做电子密度差，但是只有HF的质量，而HF比DFT还差，所以你的CCSD(T)相当于白跑了，还不如跑个M06-2X/cc-pVTZ。
这个报错是说高斯不支持CCSD(T)水平下的密度。按我的理解，如果只是定性说明问题，M06-2X足够了，甚至基组也不需要cc-pVTZ，类似6-311G(d,p)级别的基组就足够好了。
至于用了NBO结果会变化，我之前没有画过加了NBO的输出文件的密度，所以不是太清楚发生了什么，怀疑是计算密度的时候错把NBO按NO的公式去算密度了。如果是这样的话，那么不能用，不知道有没有人能科普一下这种情况下读到的密度到底是什么。不过总之不写NBO肯定是更保险的。