计算化学公社

标题: D+A-组成的波函数查看自旋密度图 [打印本页]

作者
Author: 赵云跳槽 时间: 2017-5-8 20:09
标题: D+A-组成的波函数查看自旋密度图
使用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电子
这个有什么好的办法？
(, 下载次数 Times of downloads: 34) (, 下载次数 Times of downloads: 30)
左边是单独展示的效果，右边是电子密度相加后的效果，实际上我想把右边的那个水分子变成蓝色

作者
Author: sobereva 时间: 2017-5-8 21:10
把运算符从+改成-就行了

作者
Author: 赵云跳槽 时间: 2017-5-8 22:26

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

还真可以，太有喜感了

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

作者
Author: alonewolfyang 时间: 2017-5-9 11:03
减号是神么意思

作者
Author: 赵云跳槽 时间: 2017-5-9 11:07

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

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

作者
Author: 冰释之川 时间: 2017-5-9 11:30

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

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

作者
Author: sobereva 时间: 2017-5-9 12:47

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

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

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

作者
Author: 赵云跳槽 时间: 2017-5-10 00:02
谢谢两位的讨论和回答

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