计算化学公社

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

[Multiwfn使用咨询] fukui函数测亲电性问题

[复制链接 Copy URL]

22

帖子

0

威望

35

eV
积分
57

Level 2 能力者

sob老师,我想再请教下关于fukui函数测算亲电性,亲核性的具体。我已经可以测算出N,N+1,N-1分子对应的每个原子的Hirshfeld电荷,但是我不清楚您给的图中f-,f+底下对应的每个原子的Hirshfeld电荷是怎么来的(下图红色标出的),我按照浮亏函数直接先测算亲电性,让后返回操作界面后直接输入 7 // Population analysis,1 // Hirshfeld population,但是得出的结果显示的原子电荷值跟单独的N是一样的,请问老师到底应该怎么操作啊

求助.png (112.65 KB, 下载次数 Times of downloads: 65)

求助.png

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

2#
发表于 Post on 2016-5-20 23:04:44 | 只看该作者 Only view this author
看手册4.7.3节的例子
北京科音自然科学研究中心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

22

帖子

0

威望

35

eV
积分
57

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2016-5-29 17:05:02 | 只看该作者 Only view this author
sobereva 发表于 2016-5-20 23:04
看手册4.7.3节的例子

又打搅了sob老师,我看了,但是没有写出怎么得到对应的具体原子N-1,N+1状态下的Hirshfeld电荷,我也翻了下4.5.4节,也没有发现额,求老师可以具体一点么,我卡在这块了

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

4#
发表于 Post on 2016-5-29 18:26:50 | 只看该作者 Only view this author
miya 发表于 2016-5-29 17:05
又打搅了sob老师,我看了,但是没有写出怎么得到对应的具体原子N-1,N+1状态下的Hirshfeld电荷,我也翻了 ...


你载入N-1和N+1状态的.wfn文件,按照手册4.7节的例子就算出电荷来了,超级容易
会算中性的不可能不会算离子态的,操作过程完全一样,仅输入文件不一样而已
北京科音自然科学研究中心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

22

帖子

0

威望

35

eV
积分
57

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2016-5-30 15:22:31 | 只看该作者 Only view this author
sobereva 发表于 2016-5-29 18:26
你载入N-1和N+1状态的.wfn文件,按照手册4.7节的例子就算出电荷来了,超级容易
会算中性的不可能不会 ...

谢谢sob老师,搞定啦

78

帖子

0

威望

814

eV
积分
892

Level 4 (黑子)

6#
发表于 Post on 2018-4-24 12:06:54 | 只看该作者 Only view this author
借楼主问的Fukui函数预测反应亲电性的问题再请教下Sob老师。博文http://sobereva.com/402这段有些迷糊。。先在Gaussian中在B3LYP/6-31G*下优化中性状态的晕苯,将得到的chk文件转换为fch后命名为N.fch。然后基于优化的中性的晕苯的结构,将电荷和自旋多重度改为1 2(即+1阳离子态),算个单点任务(切勿做优化任务),将得到的chk转换并命名为N-1.fch
单点计算N-1的时候基组可以与输入文件为N时的一致么,只是不优化。比如输入文件为N时用的#p m062X/6-311+G(d,p) opt freq   N-1时直接写 #p m062x/6-311+G(d,p),还是得算单点能提高一个级别才可行呢?谢谢老师
自己决定自己的DNA

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

7#
发表于 Post on 2018-4-24 14:19:20 | 只看该作者 Only view this author
QuantumHan 发表于 2018-4-24 12:06
借楼主问的Fukui函数预测反应亲电性的问题再请教下Sob老师。博文http://sobereva.com/402这段有些迷糊。。 ...


需要一致
北京科音自然科学研究中心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

78

帖子

0

威望

814

eV
积分
892

Level 4 (黑子)

8#
发表于 Post on 2018-4-24 15:16:37 | 只看该作者 Only view this author

老师不好意思,我没表述清楚。我是想问,优化用的#p m062X/6-311+G(d,p) opt freq . 将优化后的N和N-1两个文件保存成两个单点输入文件,这两个计算单点能否还是仍用#p m062X/6-311+G(d,p)呢,还是需要提升一个级别比如用双杂化泛函?因为在算势能面能垒的时候用的后HF方法,而后HF方法计算又显得比较昂贵,只是想得到这个N-1的.chk文件作分析
自己决定自己的DNA

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

9#
发表于 Post on 2018-4-24 17:09:58 | 只看该作者 Only view this author
QuantumHan 发表于 2018-4-24 15:16
老师不好意思,我没表述清楚。我是想问,优化用的#p m062X/6-311+G(d,p) opt freq . 将优化后的N和N-1两 ...


#p m062X/6-311+G(d,p) opt freq优化N,直接就得到了N.fch
然后再用N.fch保存成新的.gjf文件,设成N-1状态的单点任务,还是用m062X/6-311+G(d,p),算完了就有了N-1.fch。通过这俩fch求差得到福井函数。
没必要用后HF。电子密度相关的量对计算级别敏感度远远低于能量。
北京科音自然科学研究中心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

78

帖子

0

威望

814

eV
积分
892

Level 4 (黑子)

10#
发表于 Post on 2018-4-27 20:44:44 | 只看该作者 Only view this author
sobereva 发表于 2018-4-24 17:09
#p m062X/6-311+G(d,p) opt freq优化N,直接就得到了N.fch
然后再用N.fch保存成新的.gjf文件,设成N-1 ...

谢谢老师的指导。
自己决定自己的DNA

188

帖子

0

威望

922

eV
积分
1110

Level 4 (黑子)

11#
发表于 Post on 2019-12-21 14:32:26 | 只看该作者 Only view this author

老师为了更准确我可以用b3lyp/genecp优化出来的结构,然后把N,N-1状态下都在b2plypd3/def2tzvp 基组下计算单点能来考查亲电位点吗?

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

12#
发表于 Post on 2019-12-22 00:36:09 | 只看该作者 Only view this author
liuyjhx 发表于 2019-12-21 14:32
老师为了更准确我可以用b3lyp/genecp优化出来的结构,然后把N,N-1状态下都在b2plypd3/def2tzvp 基组下计 ...

没意义
算电子能量更准的不代表算的波函数也更准,何况波函数分析对波函数质量要求并不高,用普通泛函足够了。而且本来Gaussian算双杂化就很慢,产生双杂化的密度更耗时
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-17 01:43 , Processed in 0.179281 second(s), 23 queries , Gzip On.

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