计算化学公社

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

[波函数分析交流] Bond Polarity Index(BPI)怎么算呢

[复制链接 Copy URL]

27

帖子

0

威望

199

eV
积分
226

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
看到文献算[img=0,1]file:///C:\Users\junzhu\AppData\Roaming\Tencent\Users\1135261755\QQ\WinTemp\RichOle\28PI2PK%%LES5%TGD5P[P3J.png[/img] ,这个EI在高斯里该如何计算呢~附EI理论背景一张 ,谢谢大家啦

6万

帖子

99

威望

5万

eV
积分
120100

管理员

公社社长

2#
发表于 Post on 2015-11-10 14:10:00 | 只看该作者 Only view this author
最好把文献直接传上来
北京科音自然科学研究中心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

27

帖子

0

威望

199

eV
积分
226

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2015-11-10 18:27:04 | 只看该作者 Only view this author
sobereva 发表于 2015-11-10 14:10
最好把文献直接传上来

好哒!

27

帖子

0

威望

199

eV
积分
226

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2015-11-10 18:29:56 | 只看该作者 Only view this author
附件是关于计算BPI的理论方法

j100377a037.pdf

791.74 KB, 下载次数 Times of downloads: 13

27

帖子

0

威望

199

eV
积分
226

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2015-11-10 19:32:08 | 只看该作者 Only view this author

Sob老师,您是要写进Multiwfn嘛?

6万

帖子

99

威望

5万

eV
积分
120100

管理员

公社社长

6#
发表于 Post on 2015-11-11 04:54:12 | 只看该作者 Only view this author
文献里那公式纯粹是把问题搞复杂了,让人看起来糊涂,其实原理非常简单。BPI只要你得到了EI就能马上得到,这里就说EI。Multiwfn极为灵活强大,算这个都不需要专门写新的代码。

f_ij那个量其实就是i分子轨道中j原子轨道的成份而已,也不知道作者怎么想的非要写成那个怪模样。用Multiwfn直接就得到了轨道成份,见《谈谈轨道成份的计算方法》(http://sobereva.com/131),手册4.8节有大量例子。文中用的计算轨道成份方法其实就是Mulliken方法,实际上你用Multiwfn支持的其它计算原子所占成分的方法也可以,诸如Hirshfeld。

EI_A的分子项就是:A原子在每个价层分子轨道中的轨道成份乘上这个分子轨道的占据数和轨道能,然后对所有价层分子轨道加和。用Multiwfn算出一批轨道中某原子的轨道成份后,拷到excel里,然后用主功能6的选项3可以显示所有分子轨道的能量,也拷到excel里,excel里编辑个简单公式让它们相乘再乘上轨道占据数后一加和立马就得到了。
EI_A的分母项就是:A原子在所有价层轨道中的布居数之和。在Multiwfn主功能6里选34就能把内层轨道占据数清零,然后你进主功能7用与获得轨道成份相同的方法做布居分析就行了(比如Hirshfeld算的轨道成份你就用Hirshfeld布居分析)。

EI_A的分子项和分母项一除,物理意义就是A原子平均每个价电子的能量。其实这和Multiwfn支持的平均局部离子化能(ALIE)极其相似,只不过ALIE是考虑所有分子轨道,并且考虑的是三维空间中每个点的电子能量。

BPI_AB反映的就是A、B两个原子的平均价电子能量的差异,由此试图衡量键的极性。但其实衡量键极性的方法有的是,没必要非得用这个定义。诸如NLMO键级就能反映键的极性,可以参看J. Phys. Chem. A, 117, 3100−3108 (2013)。
北京科音自然科学研究中心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
积分
120100

管理员

公社社长

7#
发表于 Post on 2015-11-11 06:22:45 | 只看该作者 Only view this author
花了也就半小时,写了50行代码,已经在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

27

帖子

0

威望

199

eV
积分
226

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2015-11-11 09:01:43 | 只看该作者 Only view this author
sobereva 发表于 2015-11-11 04:54
文献里那公式纯粹是把问题搞复杂了,让人看起来糊涂,其实原理非常简单。BPI只要你得到了EI就能马上得到, ...

爱死您了!!!

27

帖子

0

威望

199

eV
积分
226

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2015-11-11 09:42:41 | 只看该作者 Only view this author
sob老师,其实我是在找group electronegativity 的时候追踪到这篇文献的,附件的第二篇文献讲的就一种group electronegativity算法,文献中也提到了其他很多算法,您看这个group electronegativity 合不合适在Multiwfn中实现呢

j100180a032.pdf

920.39 KB, 下载次数 Times of downloads: 9

BPI: Application to group electronegativity

ja00311a008.pdf

563.59 KB, 下载次数 Times of downloads: 11

Calculation of group electronegativity

6万

帖子

99

威望

5万

eV
积分
120100

管理员

公社社长

10#
发表于 Post on 2015-11-11 17:52:51 | 只看该作者 Only view this author
郝玉蕾 发表于 2015-11-11 09:42
sob老师,其实我是在找group electronegativity 的时候追踪到这篇文献的,附件的第二篇文献讲的就一种group ...

这个概念我不打算implement进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
积分
120100

管理员

公社社长

11#
发表于 Post on 2015-11-14 00:37:48 | 只看该作者 Only view this author
已经在主页上更新了3.3.8(dev)版,可以算EI和BPI了,见手册3.200.12节的介绍和4.200.12节的例子,可以重复出原文的结果。基团电负性后来看了看就是算自由基的链接原子的EI而已,也是直接能算的。
北京科音自然科学研究中心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

27

帖子

0

威望

199

eV
积分
226

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2016-5-19 10:43:26 | 只看该作者 Only view this author
sobereva 发表于 2015-11-14 00:37
已经在主页上更新了3.3.8(dev)版,可以算EI和BPI了,见手册3.200.12节的介绍和4.200.12节的例子,可以重复 ...

真抱歉,好久没上来了,真的谢谢您!马上就硕士毕业了,虽然不再继续做科研,但是真的真的特别感谢sob老师,感谢思想家公社群,感谢计算化学公社!祝福您身体健康,科研顺心,衷心祝愿Mutiwfn越来越强大,我们课题组一直在用哦,也发了不少文章

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 21:09 , Processed in 0.191340 second(s), 24 queries , Gzip On.

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