计算化学公社

 找回密码 Forget password
 注册 Register
Views: 18492|回复 Reply: 17
打印 Print 上一主题 Last thread 下一主题 Next thread

[波函数分析求助] 做电子密度差分图时考虑的跃迁组成

[复制链接 Copy URL]

294

帖子

2

威望

3012

eV
积分
3346

Level 5 (御坂)

请问在做反映基态向激发态跃迁过程的电子密度差分图时,是用激发态的总电子密度减去基态的吗?还是只考虑跃迁组成不低于2%的?
在进行TD-DFT计算时,每个电子态跃迁都可能对应一系列轨道跃迁组合,每个轨道跃迁组合的跃迁系数平方乘以二,就是跃迁组成。
但输出文件(log)文件中只会列出跃迁组成不低于2%的,所以做电子密度差分图时,是否不考虑跃迁组成小于2%的?

6万

帖子

99

威望

5万

eV
积分
120073

管理员

公社社长

2#
发表于 Post on 2015-1-7 22:46:48 | 只看该作者 Only view this author
如果你是通过Multiwfn和.wfn文件,令激发态的密度减去基态的密度,那么可以认为所有轨道跃迁全都考虑了。(实际情况更复杂些,默认情况下Gaussian产生的激发态密度是弛豫的密度,因此还有额外的处理,但这就是题外话了)
如果你是通过Multiwfn的electron-hole分析功能来绘制密度差图,那么考虑多少轨道跃迁,就看Gaussian输出文件里记录了多少轨道跃迁。默认只输出系数绝对值大于0.1的,此时密度差图不准确,需要用比如IOp(9/40=4)输出更多的系数结果才准确。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

294

帖子

2

威望

3012

eV
积分
3346

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2015-1-7 22:55:26 | 只看该作者 Only view this author
就是前者,用Multiwfn得到cub文件,然后用GV画图。顺便问一下,用GaussView做这种图时蓝色是电子密度减小,紫色是电子密度增加吗?

6万

帖子

99

威望

5万

eV
积分
120073

管理员

公社社长

4#
发表于 Post on 2015-1-7 23:19:10 | 只看该作者 Only view this author
你用Multiwfn产生电子密度cube文件,放到gv里显示等值面,看到的颜色就是正值区域的颜色。然后密度差的等值面图中哪些区域是正是负就清楚了。
正的就是激发后电子密度增加区域,负的就是减小区域
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

294

帖子

2

威望

3012

eV
积分
3346

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2015-2-6 22:38:40 | 只看该作者 Only view this author
sobereva 发表于 2015-1-7 23:19
你用Multiwfn产生电子密度cube文件,放到gv里显示等值面,看到的颜色就是正值区域的颜色。然后密度差的等值 ...

试过之后发现,正值的电子密度图就是紫色的。在这里写下,希望对有类似问题的后来者有帮助。

6万

帖子

99

威望

5万

eV
积分
120073

管理员

公社社长

6#
发表于 Post on 2015-2-7 08:17:37 | 只看该作者 Only view this author
preference-colors-surface colors里也可以设定对于不同情况的等值面的颜色,左侧的颜色就是正值,右侧是负值。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

144

帖子

0

威望

2291

eV
积分
2435

Level 5 (御坂)

7#
发表于 Post on 2015-3-7 06:25:32 | 只看该作者 Only view this author
Multiwfn激发态跃迁过程的电子密度差分图 有相关发表的文献吗?若有,能给个链接吗?

294

帖子

2

威望

3012

eV
积分
3346

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2015-4-23 13:34:53 | 只看该作者 Only view this author
sobereva 发表于 2015-1-7 22:46
如果你是通过Multiwfn和.wfn文件,令激发态的密度减去基态的密度,那么可以认为所有轨道跃迁全都考虑了。( ...

用IOp(9/40=3)是否可以?
之前没注意看,一直写成IOp(9/40=3)。是否需要重算?

6万

帖子

99

威望

5万

eV
积分
120073

管理员

公社社长

9#
发表于 Post on 2015-4-23 14:54:48 | 只看该作者 Only view this author
虎王 发表于 2015-4-23 13:34
用IOp(9/40=3)是否可以?
之前没注意看,一直写成IOp(9/40=3)。是否需要重算?

对于图形描述一般不会有什么问题,看不出差异来。

如果你计算electron和hole格点数据后屏幕显示的积分值偏离1比较明显,说明设成3有点小,如果偏离1不明显就无所谓。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

5万

eV
积分
120073

管理员

公社社长

10#
发表于 Post on 2015-4-23 14:59:20 | 只看该作者 Only view this author
ZHANGZY 发表于 2015-3-7 06:25
Multiwfn激发态跃迁过程的电子密度差分图 有相关发表的文献吗?若有,能给个链接吗?

这样的文章太多了,目前引用Multiwfn的600篇文章里大抵有二三十篇都是这样做的。由于文章太多,具体是哪些文章这么做我一时没法完全找出来。本版有个引用multiwfn的451篇文章的合集可以下载,可以用acrobat的目录里全部pdf搜索功能在里面搜。

手头能直接找到的是这两篇:
J Mol Model (2014) 20:2309
DOI: 10.1039/C4CP03355H
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

294

帖子

2

威望

3012

eV
积分
3346

Level 5 (御坂)

11#
 楼主 Author| 发表于 Post on 2015-4-23 15:05:45 | 只看该作者 Only view this author
sobereva 发表于 2015-4-23 14:54
对于图形描述一般不会有什么问题,看不出差异来。

如果你计算electron和hole格点数据后屏幕显示的积分 ...

您指的是看以下数据:
There are   20 transition modes, analyze which one?  e.g. 2
1
Multiplicity of this excited state is   3
Excitation energy is   3.9452000 eV
There are    3179 orbital pairs in this transition mode
The sum of square of excitation coefficients:  0.508537
The negative of the sum of square of de-excitation coefficients: -0.009812
The sum of above two values  0.498725

6万

帖子

99

威望

5万

eV
积分
120073

管理员

公社社长

12#
发表于 Post on 2015-4-23 15:39:24 | 只看该作者 Only view this author
虎王 发表于 2015-4-23 15:05
您指的是看以下数据:
There are   20 transition modes, analyze which one?  e.g. 2
1

不是这个,是指之后选0计算格点数据之后显示的。
但你看你提贴的这个也可以,最后一行偏离0.5太大说明离归一化条件偏得较远。当前偏离比较小,可以忽略。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

294

帖子

2

威望

3012

eV
积分
3346

Level 5 (御坂)

13#
 楼主 Author| 发表于 Post on 2015-4-23 20:40:57 | 只看该作者 Only view this author
sobereva 发表于 2015-4-23 15:39
不是这个,是指之后选0计算格点数据之后显示的。
但你看你提贴的这个也可以,最后一行偏离0.5太大说明离 ...

如果The sum of above two values 超过0.498,是否就可以接受?

6万

帖子

99

威望

5万

eV
积分
120073

管理员

公社社长

14#
发表于 Post on 2015-4-23 21:04:20 | 只看该作者 Only view this author
虎王 发表于 2015-4-23 20:40
如果The sum of above two values 超过0.498,是否就可以接受?

可以
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

294

帖子

2

威望

3012

eV
积分
3346

Level 5 (御坂)

15#
 楼主 Author| 发表于 Post on 2015-4-26 11:08:04 | 只看该作者 Only view this author

那加了IOP(9/40=4)之后,有的激发态会出现 The sum of above two values  0.500002,也就是超过了50%,这也正常吧。

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-13 00:14 , Processed in 0.160785 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list