计算化学公社

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

[波函数分析求助] 福井函数和简缩福井函数计算结果矛盾求助

[复制链接 Copy URL]

6

帖子

0

威望

33

eV
积分
39

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
结构优化之后按照Multiwfn计算福井函数,给出等势面图,然后根据NPA电荷计算简缩福井函数,发现简缩福井函数大的地方等势面包裹并不是最充分的地方,按照理论上说不是应该等势面包裹越充分的地方简缩福井函数值越大吗?

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

2#
发表于 Post on 2020-1-21 00:16:55 | 只看该作者 Only view this author
计算简缩福井函数一定要用Hirshfeld电荷,不要用NPA。仔细看
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484http://bbs.keinsci.com/thread-13328-1-1.html
TCA上的一篇对比不同原子电荷预测反应位点、亲电/亲核性的文章
http://bbs.keinsci.com/thread-15512-1-1.html

那叫等值面图不叫等势面图,你计算的不是“势”

等值面看到的效果和等值面数值(isovalue)选取有关,而且要注意正负部分的贡献都得考虑,从图形上估计净值大小而不要只看正值部分

如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具体内容(仔细看http://bbs.keinsci.com/thread-9348-1-1.html),我已把你的不恰当标题“福井函数和简缩福井函数”改了,以后务必注意

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

帖子

0

威望

33

eV
积分
39

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2020-1-21 00:45:27 | 只看该作者 Only view this author
发帖我没注意看规则,我下次注意;另外,我换成hirshfeld电荷之后区别巨大,数据合理了很多,但是依旧不解决等值面和cFF的矛盾问题,我结构优化和电荷计算都是b3lyp/6-31+g(d,p)下做的,也并没有虚频,但是在输出等值面的时候无论怎么调整,正负值都看,依旧是cFF最大的地方没有等值面,而第二大的地方等值面明显,换成dual描述以后和Fukui的描述一致,都和cFF矛盾,这种情况理论上合理吗?

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

4#
发表于 Post on 2020-1-21 01:06:37 | 只看该作者 Only view this author
yynsz123 发表于 2020-1-21 00:45
发帖我没注意看规则,我下次注意;另外,我换成hirshfeld电荷之后区别巨大,数据合理了很多,但是依旧不解 ...

明显是计算过程有误。直接用http://sobereva.com/484文中说的功能,得到的CFF不可能和FF有这么大差别。说不清楚就把此功能涉及的三个gjf文件连同三个wfn文件一起压缩后上传。
北京科音自然科学研究中心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

3097

帖子

29

威望

1万

eV
积分
17098

Level 6 (一方通行)

5#
发表于 Post on 2020-1-21 14:56:23 | 只看该作者 Only view this author
yynsz123 发表于 2020-1-21 13:01
所有相关文件都放在里面了,输出是fch应该不影响吧,我看手册里面wfn和fch都有用,但是fch太大了,gjf完了 ...

没有波函数让别人帮你跑一遍么?用网盘传一下

6

帖子

0

威望

33

eV
积分
39

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2020-1-21 20:06:25 | 只看该作者 Only view this author
我反复跑了好几遍都不一样,等值面和cFF的差别不仅在于苯环上两个氯的间位C上,还有羧基C,无论是NPA还是hirshfeld电荷计算cFF都是一样存在这个问题

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

7#
发表于 Post on 2020-1-22 13:44:49 | 只看该作者 Only view this author
yynsz123 发表于 2020-1-21 13:01
所有相关文件都放在里面了,输出是fch应该不影响吧,我看手册里面wfn和fch都有用,但是fch太大了,gjf完了 ...

文件大可以压缩啊,论坛首页明确说了上传较大文本文件前必须压缩
北京科音自然科学研究中心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
积分
120109

管理员

公社社长

8#
发表于 Post on 2020-1-22 13:45:48 | 只看该作者 Only view this author
yynsz123 发表于 2020-1-21 20:06
我反复跑了好几遍都不一样,等值面和cFF的差别不仅在于苯环上两个氯的间位C上,还有羧基C,无论是NPA还是hi ...

不要再提NPA电荷,反复说了这不适合当前问题

直接用Multiwfn的现成的计算FF和CFF的功能,前面已经说了。严格照着下文操作
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484http://bbs.keinsci.com/thread-13328-1-1.html
如果想尽快解决问题,就不要做任何自行发挥,否则只会白费更多时间。

我不回复与文中操作不符的情况,尤其是自己写什么pop=hirshfeld乱七八糟的我博文或手册里都根本没有提到的做法
如果按照文中操作还有异常,把wfn文件打包上传。

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

帖子

0

威望

33

eV
积分
39

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2020-1-22 17:54:30 | 只看该作者 Only view this author
如帖子一步一步做的,还是等值面和F0对不上,而且结果和之前手动的差不多……

N.wfn

909.25 KB, 下载次数 Times of downloads: 3

N 1.wfn

1.77 MB, 下载次数 Times of downloads: 4

N-1.wfn

1.75 MB, 下载次数 Times of downloads: 3

CDFT_DCF-.txt

7.22 KB, 下载次数 Times of downloads: 5

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

10#
发表于 Post on 2020-1-22 19:07:20 | 只看该作者 Only view this author
你上传的N 1.wfn是什么文件?命名一定要清楚明白。用这个当文件名,Multiwfn的那个功能不可能能正常运行。

再次强调:
上传较大文本文件前必须压缩

另外你在6L贴的文件也莫名其妙,如果N状态就对应-1,怎么可能你标的N-1里的dcf_N 1.gjf里电荷成了-2?明明应当是0才对。
一定要注意N是指的电子数,不要在描述上引起任何混淆。


北京科音自然科学研究中心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
积分
120109

管理员

公社社长

11#
发表于 Post on 2020-1-22 19:28:50 | 只看该作者 Only view this author
倘若你之前传的N 1.wfn是N+1.wfn,结果完全没什么异常的

f+的0.001等值面,CFF最大数值完全对应于图中的情况,即H21附近具有很大的f+。那偏离原子中心的一大坨体现出电子在里德堡轨道上。(本身-1阴离子态结合电子能量就很弱,显然很难再结合一个电子,再结合的话必然容易分布在距离体系很远的地方。这体系-2电荷时HOMO轨道能量都是正值,因此实际上也根本不能稳定存在)



f- 0.005等值面。结果完全正常。虽然等值面图上N7看起来f-数值也较大,但负值部分很大程度抵消了正值部分,需要反复切换不同isovalue才能充分了解这一点。这也体现出直接看图不如直接看定量的数值便于考察。

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

帖子

0

威望

33

eV
积分
39

Level 2 能力者

12#
 楼主 Author| 发表于 Post on 2020-1-22 20:20:13 | 只看该作者 Only view this author
本帖最后由 yynsz123 于 2020-1-22 20:36 编辑

原来是这么回事,本来这个分子实际上是阴离子形态,所以我就使用了负离子作为最稳状态,论坛上载的时候+号显示不出来;另外,用中性分子且换成multiwfn算hirshfeld电荷之后结果就对上了,唯一对不上的原来是等值面选的问题。OK,全部搞定!祝老师新年快乐!

6

帖子

0

威望

33

eV
积分
39

Level 2 能力者

13#
 楼主 Author| 发表于 Post on 2020-1-23 08:22:57 | 只看该作者 Only view this author
追加两个问题:1,这种N+1完全不稳定的状态是否影响f+的理论意义?2. 羧基氧和苯环氢较大的f0不知道如何理解,其中苯环氢H21是和上一个问题联系在一起的,另外,羧基氧在讨论自由基反应,尤其是羟基自由基攻击的时候是被许多文献自动忽略的,这种忽略是合理的吗?

25

帖子

0

威望

405

eV
积分
430

Level 3 能力者

14#
发表于 Post on 2021-9-26 14:48:44 | 只看该作者 Only view this author
老师,您好,福井函数算出来的反应活性与化合物本身发生某一反应的反应能垒个别没有相关性,请问还需要考虑哪些因素?除了空间位阻,反应物的稳定性是否与福井函数的结果是否有重叠?

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

15#
发表于 Post on 2021-9-29 15:42:24 | 只看该作者 Only view this author
yinlina 发表于 2021-9-26 14:48
老师,您好,福井函数算出来的反应活性与化合物本身发生某一反应的反应能垒个别没有相关性,请问还需要考虑 ...

首先要明白,福井函数只适合分子内不同位点间对比,不同分子之间要用其它指标,诸如局部软度、ALIE、静电势、原子电荷等
而且即便是对同一份子,也不要总是指望能有很好的相关性,毕竟福井函数的思想很简单。如果能直接算出来势垒,发现福井函数(或其它指标)能解释差异就解释,不能解释就不要强行讨论。
北京科音自然科学研究中心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-14 13:35 , Processed in 0.167674 second(s), 22 queries , Gzip On.

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