计算化学公社

标题: 求助Au2电子密度差图的分析 [打印本页]

作者
Author: lzeng 时间: 2025-6-10 09:18
标题: 求助Au2电子密度差图的分析
各位老师，本人最近在研究如下图所示的结构（Au2），通过LC-ωPBE-GD3BJ/def2qzvpp优化后，获取了体系（au2.fchk)和两个片段(au-1.fchk,au-2.fchk)的波函数文件，分析得到下图所示的电子密度差图，有没有老师可以指点一下，帮忙分析一下图中峰形成的原因。谢谢各位老师啦！

分析步骤：
输入au2.fchk
3
0
2
-,C:\Users\Administrator\Desktop\CDD\au2\x\def2qzvpp\au-1.fchk
-,C:\Users\Administrator\Desktop\CDD\au2\x\def2qzvpp\au-2.fchk
1
2
0,-1.24,1.81,4.71,-1.24,1.81
输出图片

作者
Author: sobereva 时间: 2025-6-10 12:29
我重算了一遍fch文件，结果没任何问题（G16 B.01，Multiwfn目前官网上最新版本）

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

相关文件
(, 下载次数 Times of downloads: 6)
绘制的是0,0,-4 Bohr到0,0,4 Bohr

作者
Author: lzeng 时间: 2025-6-10 15:24

sobereva 发表于 2025-6-10 12:29
我重算了一遍fch文件，结果没任何问题（G16 B.01，Multiwfn目前官网上最新版本）

感谢社长的回复！我还有一点疑惑，如果我的操作没有问题，那么不同点的就是坐标不一样了，我是让Au-Au键平行于x轴，而社长您给出的是平行于z轴，为什么坐标的变化影响这么大呢，我不太理解

作者
Author: sobereva 时间: 2025-6-10 21:41

lzeng 发表于 2025-6-10 15:24
感谢社长的回复！我还有一点疑惑，如果我的操作没有问题，那么不同点的就是坐标不一样了，我是让Au-Au键 ...

跟平行于哪个轴没直接关系

作者
Author: lzeng 时间: 2025-6-12 10:05
本帖最后由 lzeng 于 2025-6-12 10:06 编辑

sobereva 发表于 2025-6-10 21:41
跟平行于哪个轴没直接关系

好的，还想请问一下图中峰形成的原因是什么，麻烦社长了

作者
Author: sobereva 时间: 2025-6-13 02:39
我检查了一下，发现之前得到的Au的fch文件不对，在默认初猜下SCF没有收敛到正确波函数。应该以如下方式调换初猜轨道（以计算第一个Au为例），否则会收敛到4d9 5s2

%chk=C:\gtest\1.chk
#p lc-wpbe/def2qzvp guess=alter nosymm
Title Card Required
0 2
Au 0.00000000 0.00000000 1.24700009
5 10

复制代码

当前的波函数文件： (, 下载次数 Times of downloads: 4)

Multiwfn绘制的Au-Au键轴上的变形密度曲线图如下。可见这回正常了，在两个Au之间有电子密度增加区域，体现了共享电子作用
(, 下载次数 Times of downloads: 68)

作者
Author: lzeng 时间: 2025-6-13 15:20

sobereva 发表于 2025-6-13 02:39
我检查了一下，发现之前得到的Au的fch文件不对，在默认初猜下SCF没有收敛到正确波函数。应该以如下方式调换 ...

好的，感谢您

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