计算化学公社

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

[新手求助] 关于电子密度的问题

[复制链接 Copy URL]

8

帖子

0

威望

195

eV
积分
203

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 zengxingye 于 2016-8-18 16:32 编辑

各位老师:有问题请教,在一篇文献中看到使用PM3计算苯并噻吩中硫的电子密度为5.739,请问他这是用什么方法得到电子密度,因为我使用Mulliken和NPA电荷去分析时,分别才0.245和0.421.

ef000096i.pdf

250.58 KB, 下载次数 Times of downloads: 25

BT-PM3.log

108.05 KB, 下载次数 Times of downloads: 6

BT.log

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

99

帖子

0

威望

1092

eV
积分
1191

Level 4 (黑子)

2#
发表于 Post on 2016-8-18 16:55:32 | 只看该作者 Only view this author
两个概念不同,第一个是电子密度,指的是空间中某一点上电子密度,第二个是原子电荷,对一个原子周围的电子密度积分得到

评分 Rate

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

查看全部评分 View all ratings

3097

帖子

29

威望

1万

eV
积分
17098

Level 6 (一方通行)

3#
发表于 Post on 2016-8-18 21:56:44 | 只看该作者 Only view this author
电子密度是三维函数,用某个值表示肯定要说明是什么值,如某点的值、某区域的积分等等,而你的全文中都未说明它到底考察的是什么。

唯一沾边的描述是“electron densities on the sulfur atom”,或许可以指S原子核处的电子密度,但一方面核处的电子密度对化学反应几乎没有意义(实际上电子密度本身对分析化学反应就用处很小,有机化学里那些按照“电子密度高低”分析反应性根本就是在意淫,实际做个电子密度分布图就知道了),另一方面核处的电子密度一般是很大的数,不会有“5.739”这样,应该是1E+4、1E+5的样子。

另外文中希望讨论的问题是被氧化的难易程度,这在热力学上直接推测产物、计算氧化电势就可以了。不知道产物可以用电离势IP作为一个近似,如果嫌IP算起来稍麻烦,用HOMO做一个很粗糙的近似也比这文章用什么电子密度瞎扯淡好。

评分 Rate

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

查看全部评分 View all ratings

6万

帖子

99

威望

5万

eV
积分
120157

管理员

公社社长

4#
发表于 Post on 2016-8-18 22:02:09 | 只看该作者 Only view this author
作者明显就把密度和布居数给彻底搞混了
估计他说的是硫的电子布居数。(由于半经验方法忽略了内核电子,所以布居数是接近于其价电子数6的值)

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
liyuanhe211 + 2 ( ⊙ o ⊙ )原来是内层省了

查看全部评分 View all ratings

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

8

帖子

0

威望

195

eV
积分
203

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2016-8-19 09:58:24 | 只看该作者 Only view this author
liyuanhe211 发表于 2016-8-18 21:56
电子密度是三维函数,用某个值表示肯定要说明是什么值,如某点的值、某区域的积分等等,而你的全文中都未说 ...

谢谢老师的耐心指导,您看得非常仔细,我们确实是想分析氧化程度的难易,我等下在后面再说明一下。

8

帖子

0

威望

195

eV
积分
203

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2016-8-19 10:21:36 | 只看该作者 Only view this author
本帖最后由 zengxingye 于 2016-8-19 10:24 编辑

谢谢各位老师的指导,明白了很多,因为我们主要是做实验的,对一些概念没有弄明白。原来这个文献(Otsuki-2000)确实没有说清楚怎么来的电子密度,后来又有一个日本人 Shiraishi(文献 :Shiraishi-2002)中说明更清楚点,他是这么描述的:The net electron density is the sum of all the occupied orbitals on the sulfur atom.  意思是在硫原子上所有占据轨道净电子密度之和(当然我也不知道怎么算的,老师们知道的话,请解释一下),有点像老师们说的价电子之和了。

补充一下,这个日本人Otsuki(文献:Otsuki-2000)是最早提出这样个观点来解释这个噻吩化合物氧化难易的,后面很多人都根据他这个观点来解释实验现象(文献:Zhang-2014,Xun-2016),但是也有不少文献发现这个是相反的(文献:zhu-2015),但还是提到硫原子电子密度这个问题。我们在实验中也发现反应与这个电子密度好像关系不大,只是根据Otsuki-2000解释的人多了,搞得大家好像都得提一下不可,结果仔细来研究时,发现根本得不到他列出来的数据,才提出贴子主题这个问题。


总结的结果是做实验得和做理论的人多沟通才行,很容易张冠李戴,后面研究这个问题时再向各位老师请教,还是计算IP或者HOMO可能比较好点,免得不要有这么多误会。

Otsuki-2000.pdf

255.63 KB, 下载次数 Times of downloads: 2

Shiraishi-2002.PDF

429.63 KB, 下载次数 Times of downloads: 1

Zhang-2014.pdf

2.94 MB, 下载次数 Times of downloads: 1

Xun-2016.pdf

3.65 MB, 下载次数 Times of downloads: 1

zhu-2015.pdf

2.62 MB, 下载次数 Times of downloads: 1

3097

帖子

29

威望

1万

eV
积分
17098

Level 6 (一方通行)

7#
发表于 Post on 2016-8-19 10:38:47 | 只看该作者 Only view this author
zengxingye 发表于 2016-8-19 10:21
谢谢各位老师的指导,明白了很多,因为我们主要是做实验的,对一些概念没有弄明白。原来这个文献(Otsuki-20 ...

文中的这个参数很容易算,就是Mulliken布居,它就是总电子数减去Mulliken电荷。例如你算出来S的Mulliken电荷0.245,就是16-0.245=15.755;不过由于Sob提到的,半经验里没有内层电子,所以文中程序显示出来的是减去10个电子之后的5.755。

(Again,这参数虽然容易算但是几乎没用,从原理上也不可能与氧化电位的趋势吻合的好)

评分 Rate

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

查看全部评分 View all ratings

8

帖子

0

威望

195

eV
积分
203

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2016-8-20 11:51:34 | 只看该作者 Only view this author
liyuanhe211 发表于 2016-8-19 10:38
文中的这个参数很容易算,就是Mulliken布居,它就是总电子数减去Mulliken电荷。例如你算出来S的Mulliken ...

非常感谢!

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 06:41 , Processed in 0.199860 second(s), 29 queries , Gzip On.

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