计算化学公社

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

[波函数分析求助] 有关缩减福井函数(CFF)的计算公式疑问

[复制链接 Copy URL]

6

帖子

0

威望

13

eV
积分
19

Level 1 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 djs1209 于 2016-8-1 23:26 编辑

我看到Multiwfn 3.3.8版本的Manual 4.7.3节有关于缩减福井函数(CFF)的详细计算公式,本人初次涉及,因而在这方面有一些疑问。为什么使用手册和sober大神的诸多论文如(“亲电取代反应中活性位点预测方法的比较”、“预测亲核反应位点方法的比较”)里的CFF公式都是这样的;

而大多数文章里面都是符号相反,如4.7.3节引用的example里的J. Phys. Chem. A, 106, 3885中的计算公式为 ,符号完全是相反的。我现在计算一个分子时(hirshfeld电荷),如果按照Manual 4.7.3节CFF的计算公式,得到的部分原子f0值都是负值(但是绝对值很大),按照判断原则这些原子显然不是活性位点,这些原子却恰恰是实际反应的活性位点。如果按照其他文献中的计算方法符号相反,正好得到正值,而且数值也是最大的,恰恰呢符合实际情况的活性位点,以致现在很困惑。

6

帖子

0

威望

13

eV
积分
19

Level 1 能力者

11#
 楼主 Author| 发表于 Post on 2016-8-2 19:51:26 | 只看该作者 Only view this author
sobereva 发表于 2016-8-2 16:04
结果不会相差多少,也没有哪种泛函更合适之说

另外,condensed要翻译成简缩而不是缩减

好的,非常感谢sober老师的解答

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

10#
发表于 Post on 2016-8-2 16:04:37 | 只看该作者 Only view this author
djs1209 发表于 2016-8-2 00:18
再请教老师您一个问题,以前计算电荷时,有审稿人提出这样的问题:1.        The authors selected B3LYP. This  ...


结果不会相差多少,也没有哪种泛函更合适之说

另外,condensed要翻译成简缩而不是缩减
北京科音自然科学研究中心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

管理员

公社社长

9#
发表于 Post on 2016-8-2 16:02:23 | 只看该作者 Only view this author
djs1209 发表于 2016-8-2 00:10
谢谢老师的回答,这下我明白了,N指的是电子数。还有一个问题:就是在优化一个分子时,自动断裂成2个分子 ...

当然会影响,结构不同原子电荷自然不同
北京科音自然科学研究中心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

威望

13

eV
积分
19

Level 1 能力者

8#
 楼主 Author| 发表于 Post on 2016-8-2 00:18:33 | 只看该作者 Only view this author
sobereva 发表于 2016-8-2 00:02
根据你的文件名,体系电荷设定你都弄错了
手册和文献中N都是指电子数,SMX_N-1比中性少一个电子,体系理应 ...

再请教老师您一个问题,以前计算电荷时,有审稿人提出这样的问题:1.        The authors selected B3LYP. This functional tends to overestimate bond lengths and activation barriers. For at least one case the authors should redo calculations using of one of the plethora of modern functionals (e.g. one of the Truhlar's functionals) and compare the results.  建议增加Truhlar函数方法(如m06)进行校核。我想问老师的是如果计算CFF时,使用B3LYP函数和M06方法,哪种函数更加适合?

6

帖子

0

威望

13

eV
积分
19

Level 1 能力者

7#
 楼主 Author| 发表于 Post on 2016-8-2 00:10:03 | 只看该作者 Only view this author
sobereva 发表于 2016-8-2 00:02
根据你的文件名,体系电荷设定你都弄错了
手册和文献中N都是指电子数,SMX_N-1比中性少一个电子,体系理应 ...

谢谢老师的回答,这下我明白了,N指的是电子数。还有一个问题:就是在优化一个分子时,自动断裂成2个分子,会不会影响电荷的结果,这个时的电荷进行CFF计算还能不能作为依据了。

6

帖子

0

威望

13

eV
积分
19

Level 1 能力者

6#
 楼主 Author| 发表于 Post on 2016-8-2 00:06:09 | 只看该作者 Only view this author
sobereva 发表于 2016-8-1 23:59
你拿什么电荷算的简缩福井函数?

PS:上传较大文本文件(>2MB)前请压缩以节省论坛空间。

谢谢老师。现在这个log文件里是NBO的,hirshfeld电荷我用Multiwin也生成算了,结果都是f0为负值。我的分子在计算N-1时,会自动在7S处断开成2个分子,会不会是这种原因导致了7S处的电荷在N-1时始终小于N+1时,因而得到的f0值始终是负值?

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

5#
发表于 Post on 2016-8-2 00:02:01 | 只看该作者 Only view this author
根据你的文件名,体系电荷设定你都弄错了
手册和文献中N都是指电子数,SMX_N-1比中性少一个电子,体系理应带+1净电荷,但你设的净电荷居然是-1
你先把Multiwfn手册里计算简缩福井函数的例子完整重复一遍,以免操作步骤犯错
北京科音自然科学研究中心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

管理员

公社社长

4#
发表于 Post on 2016-8-1 23:59:14 | 只看该作者 Only view this author
djs1209 发表于 2016-8-1 23:56
非常感谢sober老师的讲解,但是为什么我计算出来的分子三种(N、N+1和N-1)所得到的f0值(以7S为例)为什 ...


你拿什么电荷算的简缩福井函数?

PS:上传较大文本文件(>2MB)前请压缩以节省论坛空间。
北京科音自然科学研究中心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

威望

13

eV
积分
19

Level 1 能力者

3#
 楼主 Author| 发表于 Post on 2016-8-1 23:56:44 | 只看该作者 Only view this author
sobereva 发表于 2016-8-1 23:41
Multiwfn手册里的q是指原子电荷,你那个文献里的Q是指布居数,当然是反着的
Multiwfn手册和我发的文章里的 ...

非常感谢sober老师的讲解,但是为什么我计算出来的分子三种(N、N+1和N-1)所得到的f0值(以7S为例)为什么是负值呢,正常这个分子里面S肯定是最大的,也是最强的活性位点。麻烦老师帮忙解惑,非常感谢。 SMX N 1 NEW.LOG (2.39 MB, 下载次数 Times of downloads: 2)
SMX N-1.LOG (4.08 MB, 下载次数 Times of downloads: 2)

SMX.LOG (3.12 MB, 下载次数 Times of downloads: 1)








6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

2#
发表于 Post on 2016-8-1 23:41:28 | 只看该作者 Only view this author
Multiwfn手册里的q是指原子电荷,你那个文献里的Q是指布居数,当然是反着的
Multiwfn手册和我发的文章里的公式都是绝对正确的。f怎么算手册里有例子。
北京科音自然科学研究中心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

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

GMT+8, 2025-8-14 04:25 , Processed in 0.202915 second(s), 24 queries , Gzip On.

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