计算化学公社

标题: 重叠密度与跃迁偶极矩密度图怎么做 [打印本页]

作者
Author: 叫我易小星 时间: 2016-10-24 13:10
标题: 重叠密度与跃迁偶极矩密度图怎么做
(, 下载次数 Times of downloads: 41)
有人做过这种图吗？看着应该是在G-view中得到的，但这种重叠密度及跃迁偶极矩密度不知道是怎么设的，求高手指点

作者
Author: 小范范1989 时间: 2016-10-24 15:10
我记得是sob老师手册 4.18章节。
理论是3.21章节。
之前我作图的时候做过，记不太清出来，你查查是哪个章节。手机回复，没看手册。

作者
Author: sobereva 时间: 2016-10-24 15:37
Multiwfn手册4.18.1节有例子
可以直接在Multiwfn里显示，也可以让Multiwfn导出相应的cube文件后在gview、VMD等第三方程序里绘制等值面

作者
Author: 叫我易小星 时间: 2016-10-24 16:13

小范范1989 发表于 2016-10-24 15:10
我记得是sob老师手册 4.18章节。
理论是3.21章节。
之前我作图的时候做过，记不太清出来，你查查是哪个章 ...

嗯嗯，谢谢

作者
Author: 叫我易小星 时间: 2016-10-24 16:14

sobereva 发表于 2016-10-24 15:37
Multiwfn手册4.18.1节有例子
可以直接在Multiwfn里显示，也可以让Multiwfn导出相应的cube文件后在gview、V ...

嗯，非常感谢!

作者
Author: 冰释之川 时间: 2016-10-24 19:48
本帖最后由冰释之川于 2016-10-24 19:52 编辑

顺便插一句，为什么 overlap density 会有两种颜色（正负）？我看 4.18.1例子里就一个颜色，代表重叠区域啊？是不是因为定义是min(hole,electron)，有可能会取到 hole部分？

作者
Author: 叫我易小星 时间: 2016-10-24 19:59

冰释之川发表于 2016-10-24 19:48
顺便插一句，为什么 overlap density 会有两种颜色（正负）？我看 4.18.1例子里就一个颜色，代表重叠区 ...

(, 下载次数 Times of downloads: 44)
这是手册上写的

作者
Author: 冰释之川 时间: 2016-10-24 20:41

叫我易小星发表于 2016-10-24 19:59
这是手册上写的

那是针对电子密度差图的，我好奇重叠密度是怎么回事，重叠密度里的黄/蓝代表什么？也是电子密度增加、减小？

作者
Author: 叫我易小星 时间: 2016-10-24 20:56

冰释之川发表于 2016-10-24 20:41
那是针对电子密度差图的，我好奇重叠密度是怎么回事，重叠密度里的黄/蓝代表什么？也是电子密度增加、 ...

确实是，我也不理解怎么会有两种颜色，文章中说的是代表正负

作者
Author: 叫我易小星 时间: 2016-10-24 21:13

sobereva 发表于 2016-10-24 15:37
Multiwfn手册4.18.1节有例子
可以直接在Multiwfn里显示，也可以让Multiwfn导出相应的cube文件后在gview、V ...

Sob老师您好，我想问下在gview里做这种图要怎么操作啊？还有就是教程做出来的重叠密度图只有一种颜色，与FIG.4中的怎么不太一样啊？

作者
Author: sobereva 时间: 2016-10-25 07:03

冰释之川发表于 2016-10-24 20:41
那是针对电子密度差图的，我好奇重叠密度是怎么回事，重叠密度里的黄/蓝代表什么？也是电子密度增加、 ...

“重叠密度”这种语焉不详，不给出明确定义的概念不要用它

作者
Author: sobereva 时间: 2016-10-25 07:06

叫我易小星发表于 2016-10-24 21:13
Sob老师您好，我想问下在gview里做这种图要怎么操作啊？还有就是教程做出来的重叠密度图只有一种颜色，与 ...

gview没法直接作出这种图，必须用Multiwfn做。至多只能用gview载入multiwfn产生的对应的cube文件再显示等值面，gview才没那能力去计算相应的数据。gview载入cube文件后在results-surface/contour界面里设定等值面的显示。

文中根本没写清楚怎么定义所谓的重叠密度的。八成他指的重叠密度就是跃迁密度。手册4.18.1节的例子给了跃迁密度图（You will see left graph below, which exhibits the transition density in real space representation那段话下面作图），明显有正负。
我特别反感文献乱用名词还不写清楚定义，造成读者confusing。

作者
Author: 叫我易小星 时间: 2016-10-25 08:34
本帖最后由叫我易小星于 2016-10-25 08:41 编辑

sobereva 发表于 2016-10-25 07:06
gview没法直接作出这种图，必须用Multiwfn做。至多只能用gview载入multiwfn产生的对应的cube文件再显示 ...

嗯，谢谢老师，还想问下分析跃迁偶极矩密度时X,Y ,Z分量怎么判断自己要选取哪个？

作者
Author: sobereva 时间: 2016-10-25 09:26

叫我易小星发表于 2016-10-25 08:34
嗯，谢谢老师，还想问下分析跃迁偶极矩密度时X,Y ,Z分量怎么判断自己要选取哪个？

up to you。
可以都考察，或者比如只有一个大的，就考察那个大的。

作者
Author: 叫我易小星 时间: 2016-10-25 12:40

sobereva 发表于 2016-10-25 09:26
up to you。
可以都考察，或者比如只有一个大的，就考察那个大的。

OK，再次感谢！

作者
Author: liuxm 时间: 2016-10-25 14:31
sob老师，我想接着这个话题问一句，在计算跃迁偶极矩的时候，加上关键词density=rhoci，td激发的时候会给出各个态的Ground to excited state transition electric dipole moments (Au)； Ground to excited state transition velocity dipole moments（Au）; Ground to excited state transition magnetic dipole moments (Au)还有最后给出的 Dipole moment (field-independent basis, Debye):
X= -0.0222 Y= -0.0027 Z= 0.0029 Tot= 0.0225

那么哪一个偶极矩的量才能表示基态s0->激发态s1的跃迁偶极矩呢？另外如果不加那个关键词density=rhoci，只有最后的这一项 Dipole moment (field-independent basis, Debye)不同外前三者都相同，这是为什么？这个关键词限制了什么么？谢谢sob老师

作者
Author: liuxm 时间: 2016-10-25 14:37
@sobereva

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