请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 442|回复: 5

[波函数分析求助] 不同泛函的波函数信息

[复制链接]

177

帖子

0

威望

1817

eV
积分
1994

Level 5 (御坂)

发表于 2018-1-10 17:13:13 | 显示全部楼层 |阅读模式
本帖最后由 赵云跳槽 于 2018-1-10 17:15 编辑

背景:使用DFT或MP2的fchk(实际使用SCF density或MP2 density,以前文献都是使用MP2的)作为波函数文件导入到DMA(分布式多极分析)中,得到原子的电荷,偶极,四极等信息,通过计算后发现,DFT或MP2的fchk对这些多极信息的影响不敏感,从而计算得到的分子间相互作用变化不大

解决:虽然不敏感,但还是得解释下原因
          对此,本人的解决方案是通过波函数的对应性分析来说明
          选择了两个泛函M11和B3LYP(与金标准比较的一好一坏的两泛函)
          (选取比较的是HOMO-4 ~ LUMO+4)
QQ截图20180110170526.png

问题:1. 我使用的是SCF density或MP2 density,但分析的时候是通过自然轨道或正则轨道,是否合适?
           2. 前面说过 DFT或MP2的fchk对结果不敏感,通过我图片中的数据,能否说明 “在波函数层次上就比较小,使用这些波函数产生的结果自然也小” 这样的结论?

1万

帖子

25

威望

1万

eV
积分
33013

管理员

公社社长

发表于 2018-1-10 17:36:13 | 显示全部楼层
1 不理解你的意思
2 如果你是指DFT和MP2的DMA结果相差不大,这容易理解,因为二者都考虑了电子相关作用。在我的比较原子电荷的文章中也体现出DFT和MP2结果差异不大,但是和HF差异非常大。不建议从轨道角度上讨论,本身不好分析,也不容易直接说明问题。(对不同级别的密度做密度差图,或者再取模做个积分之类,更容易考察差异程度)
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

177

帖子

0

威望

1817

eV
积分
1994

Level 5 (御坂)

 楼主| 发表于 2018-1-10 19:27:27 | 显示全部楼层
sobereva 发表于 2018-1-10 17:36
1 不理解你的意思
2 如果你是指DFT和MP2的DMA结果相差不大,这容易理解,因为二者都考虑了电子相关作用。 ...

我再问下,如果做密度差图,对于MP2方法来说,它使用的是fchk中的SCF density还是MP2 density ?

1万

帖子

25

威望

1万

eV
积分
33013

管理员

公社社长

发表于 2018-1-10 19:32:24 | 显示全部楼层
赵云跳槽 发表于 2018-1-10 19:27
我再问下,如果做密度差图,对于MP2方法来说,它使用的是fchk中的SCF density还是MP2 density ?

如果你是指Multiwfn,载入高斯直接产生的MP2任务的fch,用的是HF密度,因为此时fch里的轨道是HF的
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

177

帖子

0

威望

1817

eV
积分
1994

Level 5 (御坂)

 楼主| 发表于 2018-1-10 21:36:44 | 显示全部楼层
sobereva 发表于 2018-1-10 19:32
如果你是指Multiwfn,载入高斯直接产生的MP2任务的fch,用的是HF密度,因为此时fch里的轨道是HF的

有点迷惑,我加了density关键词,算出来有MP2 density,我使用的就是这块的密度,怎么就扯到HF密度上了

1万

帖子

25

威望

1万

eV
积分
33013

管理员

公社社长

发表于 2018-1-11 05:45:19 | 显示全部楼层
赵云跳槽 发表于 2018-1-10 21:36
有点迷惑,我加了density关键词,算出来有MP2 density,我使用的就是这块的密度,怎么就扯到HF密度上了


详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379
这里详细说了,直接载入输入文件,Multiwfn读取的只是轨道,而不是密度矩阵。轨道对应什么波函数,Multiwfn就分析的什么波函数

如果要利用fch里的不同的密度矩阵,需要按本文做法把密度矩阵专门转化成对应的轨道
在Multiwfn中基于fch产生自然轨道的方法与激发态波函数、自旋自然轨道分析实例
http://sobereva.com/403
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-8-18 16:35 , Processed in 0.197460 second(s), 26 queries .

快速回复 返回顶部 返回列表