D+A-组成的波函数查看自旋密度图

赵云跳槽 · 发表于 Post on 2017-5-8 20:09:32

使用Qchem计算D+和A-的束缚能，使用片段局域方法，但可惜的是产生的波函数文件只是最后一个片段的波函数文件
没办法，计算两次，分别得到了D+的波函数，和A-的波函数文件

现在使用Multiwfn将其组合成一个文件，使用波函数运算符号“+”，然后做自旋密度图
Multiwfn的依次操作为
D+.fchk
5
0  // 运算
1
+，A-.fchk
2 //格点选择
5  // 自旋密度图
2  // 导出cube文件，放到VMD中作图

原理上，由于两个分子都带有电荷，那么得到的自旋密度图应该呈现alpha电子在D上，beta电子在A上
但现在显示的自旋密度图就只有alpha电子，没有beta电子
这个有什么好的办法？

左边是单独展示的效果，右边是电子密度相加后的效果，实际上我想把右边的那个水分子变成蓝色

赵云跳槽 · 发表于 Post on 2017-5-10 00:02:11

谢谢两位的讨论和回答

sobereva · 发表于 Post on 2017-5-9 12:47:25

赵云跳槽发表于 2017-5-9 11:07
D+的波函数加上负的A-的波函数，这样就能够看到轨道的不同自旋方向了

应当是 D+的波函数对应的自旋密度减去 D-波函数对应的自旋密度
不是波函数直接的加减

量化程序输出的波函数里，都是令alph电子a>=beta电子，所以当前必须做减法，才能令D-体系beta电子和alpha电子调换，使其beta电子>=alpha电子

冰释之川 · 发表于 Post on 2017-5-9 11:30:52

赵云跳槽发表于 2017-5-9 11:07
D+的波函数加上负的A-的波函数，这样就能够看到轨道的不同自旋方向了

我怎么感觉 “负的A-的波函数”代表在A-的波函数里alpha自旋与beta自旋人为发生反转从而体现beta>alpha，也许默认的自旋都是alpha>beta导致做的自旋密度清一色的只看到alpha？

赵云跳槽 · 发表于 Post on 2017-5-9 11:07:34

alonewolfyang 发表于 2017-5-9 11:03
减号是神么意思

D+的波函数加上负的A-的波函数，这样就能够看到轨道的不同自旋方向了

alonewolfyang · 发表于 Post on 2017-5-9 11:03:14

减号是神么意思

赵云跳槽 · 发表于 Post on 2017-5-8 22:26:37

sobereva 发表于 2017-5-8 21:10
把运算符从+改成-就行了

还真可以，太有喜感了

sobereva · 发表于 Post on 2017-5-8 21:10:24

把运算符从+改成-就行了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Multiwfn使用咨询] D+A-组成的波函数查看自旋密度图