24小时热门版块排行榜    

查看: 2850  |  回复: 11

mengsk

金虫 (著名写手)

[交流] H和烷基链的给电子能力 已有2人参与

H和烷基链的给电子能力:
H, methyl, ethyl, hexyl, isopropyl,问了一下搞有机化学的人,都说给电子能力依次增强。
但最近我用gaussian B3LYP/6-31G(d,p) 算了一下NBO,得到的给电子能力顺序为:
hexyl=isopropyl<ethyl<methyl<H
顺序正好颠倒过来,请高手解释一下,谢谢啦!

如果计算方法有问题,请给个可以正确模拟该规律的方法和基组。
回复此楼
奋斗中的小木虫
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

zhou2009

版主 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
mengsk: 金币+5, 非常感谢!非常详细,我也想做一下电子密度差试一下,请问怎么个计算步骤?请不吝赐教 2014-07-20 17:50:23
gmy1990: 金币+5 2014-07-28 04:29:09
判断H和烷基链的给电子能力,实际上是比较它们的电负性,电负性小则给电子能力大。这在许多年前量子化学计算还高不可攀时,有大量人工计算基团电负性的工作,比较它们的电负性。现在量子化学计算普及了,算一个分子,电子分布一目了然,再也不去考虑电负性大小和给电子能力的大小了。
但这毕竟是一个基本问题,值得算一下。
可以比较H与-CH3谁的电负性大,用hf/6-311++g(d,p)算了一下CH4,并作了它与CH3、H的电子密度差。CH4中的H上竟增加了电子,表明H的电负性大于基团-CH3!
-CH3不是SP^3杂化、电负性微弱大于H吗?实际上基团-CH3的三个CH键电子对另一个H有排斥,以至于基团-CH3电负性比H的电负性还稍小。基团-CH2CH3因为排斥更大电负性更小。于是H, methyl, ethyl, hexyl, isopropyl依次电负性减小,即给电子能力依次增强。当然,烷基链再增长,排斥作用就变化不大了。
Mulliken atomic charges(净电荷)与NBO电荷,在上述问题上人为地在空间不合理划分电子归属,是有缺点不靠谱的。电子密度差是较可信的。
注:在作电子密度差时,如果等值面取默认值0.0004,看见的只是H3C-H的负值外壳,因为H3C-H比自由的H总体上电子是向内收缩的,外围是减少的。等值面取0.004时才能看见H上正的增值和C上电子的减少。
2楼2014-07-19 20:29:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
换了不同的取代基,原子和坐标都发生了变化,现在是包含了不同的取代基的分子,要重新计算、重新作密度差,各自独立。但可以将密度差相互比较。
4楼2014-07-21 14:30:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
gmy1990: 金币+2 2014-07-28 04:29:16
算一个F原子,还是要加IOP(6/7=3),算单个原子都要加。

算F原子,不可加nosymm来定位。因为F有P电子,加了nosymm,即不做对称性分析,F的P电子取向将不在坐标xyz的方向上,对进一步做电子密度差会有问题。比如算HF分子的密度差,如果HF的坐标在z轴上,我们要求F的单电子必须在z轴上。
不加nosymm,F会自动放在坐标原点,单电子在z轴上。
算HF分子时,F要放在坐标原点,H按优化的键长安排在z轴上。这时要加nosymm,将HF固定住,不进行对称性分析,否则坐标就变了。
10楼2014-07-23 12:30:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

mengsk

金虫 (著名写手)

引用回帖:
2楼: Originally posted by zhou2009 at 2014-07-19 20:29:29
判断H和烷基链的给电子能力,实际上是比较它们的电负性,电负性小则给电子能力大。这在许多年前量子化学计算还高不可攀时,有大量人工计算基团电负性的工作,比较它们的电负性。现在量子化学计算普及了,算一个分子 ...

我做过分子不同激发态的电子密度差,不同激发态时,分子中各个原子核的位置不变,只是电子密度发生变化;然而,换了不同的取代基,原子和坐标都发生变化,这是还是用同样的方法做电子密度差吗?两个分子的cube文件,相减,是这样吗?还是有别的数学处理?谢谢!
奋斗中的小木虫
3楼2014-07-21 11:22:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengsk

金虫 (著名写手)

引用回帖:
4楼: Originally posted by zhou2009 at 2014-07-21 14:30:57
换了不同的取代基,原子和坐标都发生了变化,现在是包含了不同的取代基的分子,要重新计算、重新作密度差,各自独立。但可以将密度差相互比较。

比如我要比较这两个分子:
1. backbone-H
2. backbone-CH3

按您的说法应该比较1和2的密度差。
我的疑问是,如何做出1的密度差,和2的密度差?
就拿1来说吧,是将1的哪两个电子密度进行相减呢?
谢谢!
奋斗中的小木虫
5楼2014-07-21 15:26:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
gmy1990: 金币+2 2014-07-28 04:29:25
就拿1来说吧,是将1的 backbone-H的电子密度,分别减去 backbone和H的电子密度,得到电子密度差。
这种电子密度差在multiwfn中是用计算时分别得到的各个wfn文件相减完成的。

计算时要注意:
backbone和H分别取 backbone-H中的坐标,不能优化,并且加关键词nosymm。
backbone和H分别成为有单电子的了,电荷虽然还是0,但多重度为2。
算单个原子H,还要加IOP(6/7=3)。

» 本帖已获得的红花(最新10朵)

6楼2014-07-21 17:52:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanan2015

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
请教各位大侠,一些有机分子,比如丙酮酸分子和乳酸分子的电子亲和能力可以计算出来吗?谢谢!
7楼2014-07-21 20:45:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengsk

金虫 (著名写手)

送红花一朵
引用回帖:
6楼: Originally posted by zhou2009 at 2014-07-21 17:52:09
就拿1来说吧,是将1的 backbone-H的电子密度,分别减去 backbone和H的电子密度,得到电子密度差。
这种电子密度差在multiwfn中是用计算时分别得到的各个wfn文件相减完成的。

计算时要注意:
backbone和H分别取 ...

中国良心啊,国人都这么搞科研的话还能落后吗!
又回过头来看了自己的帖子,版主的耐心和渊博,不仅解决了我目前遇到的实际问题,更拓宽了我的见识,谢谢版主!
奋斗中的小木虫
8楼2014-07-22 14:08:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengsk

金虫 (著名写手)

送红花一朵
引用回帖:
6楼: Originally posted by zhou2009 at 2014-07-21 17:52:09
就拿1来说吧,是将1的 backbone-H的电子密度,分别减去 backbone和H的电子密度,得到电子密度差。
这种电子密度差在multiwfn中是用计算时分别得到的各个wfn文件相减完成的。

计算时要注意:
backbone和H分别取 ...

按照你的方法可以做出电子密度差图。
请问IOP(6/7=3)是什么指令,为什么在gaussian IOP中找不到啊?
还有,我想算一个F原子的,应该使用哪个指令?
奋斗中的小木虫
9楼2014-07-23 11:06:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mengsk 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 8/400 2026-09-01 12:58 by ssyjh
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见