计算化学公社

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

[波函数分析求助] 苯甲醚和N,N-二甲基苯胺的亲核性强弱与计算结果自相矛盾

[复制链接 Copy URL]

3

帖子

0

威望

131

eV
积分
134

Level 2 能力者

我基于http://sobereva.com/159中的方法,通过静电势和平均局部离子化能(ALIE)判断了苯甲醚和N,N-二甲基苯胺的亲电反应位点,计算结果表明二者都是邻对位取代基,这与试验结果一致。但横向对比能量大小时发现,结果与实验并不相符,如下图所示,苯甲醚和N,N-二甲基苯胺的静电势极小值为-18.386和-16.389 kcal/mol,平均局部离子化能极小值约为8.7和9.0 eV。计算表明苯甲醚的亲核性更强,但实际情况是二甲基苯胺的电子云密度应该比苯甲醚高,亲核性更强。
后续我又看了《亲电取代反应中活性位点预测方法的比较》(http://www.whxb.pku.edu.cn/CN/abstract/abstract28694.shtml)这篇文献以及其他人的问题解答,通过http://sobereva.com/484中描述的方法计算了两个化合物的福井函数(f-)与双描述符,以及局部硬度(福井函数乘以全局软度),这些结果也并无相关性,仅全局亲核指数可以用于判断亲核性大小。CDFT结果如附件所示,C1,C5是邻位,C2,C4是间位,C3是对位。
请问这一现象如何解释?

1.png (85.81 KB, 下载次数 Times of downloads: 0)

上:静电势表面和极小点,下:ALIE极小点

上:静电势表面和极小点,下:ALIE极小点

2.png (44.07 KB, 下载次数 Times of downloads: 0)

2.png

CDFT_苯甲醚.txt

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

CDFT_N,N-二甲基苯胺.txt

5.75 KB, 下载次数 Times of downloads: 0

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

6#
发表于 Post on 2025-12-13 04:21:46 | 只看该作者 Only view this author
Yijin 发表于 2025-12-12 16:17
谢谢您的回复,我根据2楼老师的回答,重新用无虚频的结构对比了不同化合物的静电势以及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

3

帖子

0

威望

131

eV
积分
134

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2025-12-12 16:17:28 | 只看该作者 Only view this author
本帖最后由 Yijin 于 2025-12-12 16:19 编辑
SharkYYX2025 发表于 2025-12-12 14:27
可以用这个试试,看看会不会比福井函数更符合你的期望
通过轨道权重福井函数和轨道权重双描述符预测亲核和 ...

谢谢您的回复,我根据2楼老师的回答,重新用无虚频的结构对比了不同化合物的静电势以及ALIE极小值点,都能很好的与实验匹配。condensed Fukui functions 和 condensed dual descriptors对于反应位点的判断相当准确,尤其是对于取代苯体系,但两者并不能在不同化合物之间进行对比。局部硬度也不能在不同化合物之间对比。此外,Condensed local nucleophilicity index同样可以判断不同化合物之间的亲核性强弱(我的实验中包括的其他结构是含硼化合物)。您提到的轨道权重福井函数和轨道权重描述符我后续尝试能否应用在我的体系中。

93

帖子

0

威望

453

eV
积分
546

Level 4 (黑子)

开心小鲨鱼 ORCA的饲养员/食物

4#
发表于 Post on 2025-12-12 14:27:11 | 只看该作者 Only view this author
可以用这个试试,看看会不会比福井函数更符合你的期望
通过轨道权重福井函数和轨道权重双描述符预测亲核和亲电反应位点 - 思想家公社的门口:量子化学·分子模拟·二次元  http://sobereva.com/533

3

帖子

0

威望

131

eV
积分
134

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-12-12 11:31:20 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-12-12 10:53
结构是如何构建、在什么级别计算的,是否几何优化收敛且振动分析确认无虚频?当前N,N-二甲基苯胺的氮的构型 ...

感谢您的回复,事实如您所说,N,N-二甲基苯胺确实结构有两个虚频,我改进一下初始结构再计算

1208

帖子

6

威望

2430

eV
积分
3758

Level 5 (御坂)

傻傻的木瓜

2#
发表于 Post on 2025-12-12 10:53:08 | 只看该作者 Only view this author
本帖最后由 Uus/pMeC6H4-/キ 于 2025-12-12 11:17 编辑

结构是如何构建、在什么级别计算的,是否几何优化收敛且振动分析确认无虚频?当前N,N-二甲基苯胺的氮的构型不对劲,两个甲基碳原子与氮所确定的平面不应当与苯环平面呈接近垂直的大角度。编辑:参考http://bbs.keinsci.com/thread-53539-1-1.html讲的苯胺的情形。
√546=23.36664289109

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

GMT+8, 2026-1-23 23:53 , Processed in 0.240927 second(s), 25 queries , Gzip On.

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