分子间作用力变强、变弱的原因?电子结构与分子构型核心解析

说明:本文华算科技主要介绍分子间作用增强或减弱的电子与结构原因,以及溶剂、温度和浓度对分子缔合的影响。

分子间相互作用能怎样定义?

两个分子靠近时,各自的原子核和电子发生相互作用,分子所受的力随距离和取向改变。科研中讨论“结合变强”,通常指某种能量差变得更负:相互作用能描述给定构型下两个片段的净作用,缔合自由能描述特定温度和介质中的结合平衡。力的单位与能量的单位不同,势能最低构型处的净力为零。

分子间作用力变强、变弱的原因?电子结构与分子构型核心解析
图1. 含分子内氢键的双酚及参照单酚,分别与双环叔胺形成分子间氢键的缔合平衡。DOI:10.1021/jacs.2c08120

对由 A、B 组成的复合物,冻结单体内部几何后,相互作用能写为 Eint=EAB−EA−EB。其中单体保持在复合物内的构型和指定电子态。复合物与单体采用相同计算方法,能量相减才有统一参照。酚羟基与双环叔胺的氮原子形成 O—H···N 接触,就是可以明确指定两个片段的分子缔合实例。

图中 K 与 K′分别描述两组溶液缔合平衡,其差别包含周围溶剂的贡献。静态计算得到的 Eint 没有包含这些贡献,也没有包含分子结合时平移、转动和构象分布的变化。相同温度、溶剂和标准态下的自由能差适合比较溶液结合倾向;气相电子能差回答的是所选分子构型的能量变化。

自由分子结合时,原有键长和构象可能偏离单体最低能结构。相对于各自优化的单体,缔合电子能满足 ΔEassoc=Eint+Edef,其中 Edef 是分子形变所需的能量。Br 与 CF3OH、CF3SH 接近时,O—H 或 S—H 键伸长;相互作用的能量降低与单体形变的能量升高共同形成缔合能曲线的最低点。

分子间作用力变强、变弱的原因?电子结构与分子构型核心解析
图2. Br 与 CF3OH、CF3SH 复合物的距离扫描,包括完整量子化学计算、点电荷模型及分子形变能。DOI:10.1007/s00894-022-05189-6

这两组曲线中的 Nuc 项标记形变能;PC 曲线来自用点电荷替代 Br 的简化模型,其静电和极化项均可产生吸引。形变代价与交换排斥具有不同定义,不能把橙色曲线改称排斥能。柔性分子即使形成更有利的局部接触,扭转骨架所需的能量也可能抵消一部分结合收益。

DFT 的近似交换关联泛函、色散处理和基组会改变上述能量差。局域基组中,复合物能借用另一片段的基函数,产生基组重叠误差;对照校正与基组收敛可以检验这部分偏差。讨论取代基效应时,统一方法下的构型与能量对照才能把化学结构变化和计算近似造成的差别区分开。

哪些分子结构变化会增强或削弱净作用?
局部电荷分布与极化响应

氢键、卤键和芳环接触通常同时包含静电、诱导、色散及交换排斥。静电来自未相互作用片段的电荷分布;诱导来自对方电场引起的电子密度响应;色散来自电子涨落关联。SAPT 可按这些物理贡献分析净能量,一个吸引分量增加时,总相互作用仍可能变弱,因为排斥或其他分量也在改变。

吸电子取代基常使 O—H 给体的氢端更缺电子,有利于吸引电子富集受体;受体若因取代而失去电子供给能力,氢键则可能减弱。对芳环,杂原子所处位置和相对取向会改变局部静电作用。苯与吡啶形成 T 形接触时,相同的氮原子可以增强某些位移处的吸引,也可以削弱另一些位移处的吸引。

分子间作用力变强、变弱的原因?电子结构与分子构型核心解析
图3. 苯与不同取向吡啶形成 T 形二聚体的相互作用能,横轴为水平位移。DOI:10.1021/jacs.5c12769

在 T 形接触中,氮位点的局部偶极与苯电场同向时产生额外静电吸引,反向时则削弱结合,图内各水平位移处的垂直间距分别经过优化。用分子总偶极矩大小排序会遗漏这些局部几何差别。胺质子化后,原来的受体位点和片段净电荷随之改变,应按新的带电分子计算其接触,而非继续沿用中性分子的氢键排序。

接触距离、位阻与分子尺寸

在较远距离处,分子靠近常增强吸引;进入短程重叠区后,交换排斥迅速增加,分子形变所需的能量也可能上升。距离缩短并不保证净结合增强。大取代基遮挡受体位点时,分子可能无法达到有利的氢键排列;芳环侧向滑移则可能减少不利重叠,在保留接触的条件下降低总能量。

分子间作用力变强、变弱的原因?电子结构与分子构型核心解析
图4. 平行苯二聚体的二维相互作用势能面,以及苯与不同杂环沿各自最低能量路径的能量曲线。DOI:10.1021/jacs.5c12769

平行苯二聚体的二维势能面以水平位移 x 和垂直间距 z 为变量,最低能量路径在两侧形成能谷。沿 x 改变时,最优 z 也随之变化,下方曲线并非固定层间距的扫描。换成对苯醌后,局部静电作用可以使正对排列成为最低点;“芳环都应错位排列”的判断只适用于部分分子与模型。

增加共轭骨架或可极化原子,往往会增加电子云响应和色散吸引,但分子尺寸变大也会增加近程接触与排斥。只有新增部分能靠近另一分子,较大的可接触表面才可能产生额外稳定化。相对分子质量没有直接决定色散强度;极化率、激发能谱和接触构型才参与色散能的计算。组成相同的紧凑支链与伸展骨架,因接触方式不同可以具有不同结合能。

相邻氢键的协同与竞争

一个羟基接受氢键后,局部电子密度重排可增强它向下一个受体供氢的能力,产生正协同。酚链与胺的缔合实验中,新增分子内氢键会改变终端分子间氢键的强度;新增接触不必直接连到胺,仍可经邻近羟基的极化影响末端结合。

分子间作用力变强、变弱的原因?电子结构与分子构型核心解析
图5. 不同长度酚链与三类胺形成 1∶1 复合物的缔合常数对数,N 为链中的酚单元数。DOI:10.1021/jacs.2c08120

酚单元从一个增至两个时,双环叔胺对应的缔合常数大幅增加,继续延长后的增幅减小,三个与四个酚单元的结果接近。该系列体现出随链长衰减的协同增强。若分子内氢键占用了原本用于分子间结合的给体或受体,外部缔合又可能减弱;预组织、位点占用和极化需要按具体构型分别计算。

溶剂、温度和浓度为何改变分子缔合?
溶剂对给体与受体的竞争

在溶液中,两个分子结合会替换部分溶质—溶剂接触,并使溶剂重新排列。溶质间氢键的能量收益要与去溶剂化代价比较:强受氢键溶剂可与目标受体竞争给体,强供氢键溶剂可竞争受体。单个分子对在真空中吸引很强,进入溶液后仍可能只有较小的缔合常数。

分子间作用力变强、变弱的原因?电子结构与分子构型核心解析
图6. 氢键给体与受体在溶剂中的缔合,比较游离状态和结合状态的四类接触。DOI:10.1039/d1sc03392a

溶剂竞争模型用 α、β 描述给体和受体的经验氢键参数,并比较游离态与结合态的四类接触。介电常数无法单独描述位点竞争,显式溶剂分子能进一步呈现配位与氢键占位;连续介质模型主要描述平均介质响应。对疏水表面,溶剂释放与水分子重新排列也会贡献缔合自由能,不能把水中聚集的全部增益记作分子对色散能。

温度、浓度与结合比例

升温改变热运动和构象分布,结合自由能按 ΔG°=ΔH°−TΔS° 变化。对于焓变为负、温区内焓变近似不变的缔合,升温通常降低平衡常数;吸热且由熵增驱动的聚集可以随升温增强。固定几何下的电子能计算不包含这种统计平均,温度效应需由热力学修正或溶液构型采样描述。

稀溶液中,用 Kc=[AB]/([A][B]) 表示 1∶1 缔合常数。保持 Kc 不变而增加 B 的浓度,也能提高 A 的结合比例。双酚—双环叔胺滴定中,核磁共振峰随加入的胺增多而移动,多个质子的变化可用同一个结合平衡拟合。

分子间作用力变强、变弱的原因?电子结构与分子构型核心解析
图7. 正辛烷中双酚与双环叔胺的复合物结构、核磁共振滴定谱及 1∶1 结合等温线拟合。DOI:10.1021/jacs.2c08120

该实验在正辛烷、298 K 下进行,分子间 NOE 约束结合位置,化学位移随浓度趋于饱和。快交换条件下,峰位是游离态与结合态占比加权的平均值,随滴定移动无需假定每个复合物内的氢键逐渐增强。浓度较高时若出现额外聚集、溶剂组成变化或活度偏离,单一 1∶1 模型得到的表观常数才可能随拟合区间改变。

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
(0)
上一篇 1天前
下一篇 1天前

相关推荐