如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南

说明:本文华算科技主要介绍 HOMOLUMO 用于溶剂分子初筛时的计算含义、可比条件、氧化还原量换算,以及盐、局域构型、电极表面和配方性质对候选排序的修正。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南

HOMO与LUMO是什么?

HOMO 是基态分子中能量最高的占据轨道,LUMO 是能量最低的未占据轨道。前者标记电子从哪个占据态离开,后者标记外来电子优先进入哪个空态。筛选表里的数值通常来自优化构型的 Kohn–Sham 轨道本征值,单位为 eV,并以真空能级或计算程序采用的统一零点为参照。轨道等值面还会显示电子增减集中在哪个官能团,能量相近但空间分布不同的两个分子,后续断键位置可能并不相同。

在同一计算口径下,HOMO 上移通常对应较低的失电子能量,LUMO 下移通常对应较强的受电子倾向。主溶剂的筛选目标偏向较低 HOMO 和较高 LUMO;成膜添加剂则要求电子转移发生在主体溶剂之前,并生成预期界面产物。主溶剂、共溶剂和成膜添加剂应分表排序,同一组轨道数值用于不同分子角色会得到相反的取舍。

可比表至少要固定交换相关泛函、基组、色散处理、分子电荷、自旋多重度和溶剂模型。柔性分子还需保留多个低能构象,因为羰基、醚氧或砜基的转动会改变轨道离域范围。若一部分数值来自气相 B3LYP,另一部分来自隐式溶剂中的 M06-2X,把两列直接排序得到的差别可能主要来自方法。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
图1. 离子液体候选由分子池依次经过相态、离子电导率和 HOMO/LUMO 电化学窗口筛选。DOI:10.1038/s41467-023-38493-7

图 1 的 2220 种离子液体依次接受室温相态、离子电导率和电化学窗口筛选,最终保留 49 种组合。HOMO/LUMO 推得的窗口位于流程末端,用于删除明显越过电位范围的分子;黏度、盐溶解度和离子迁移由独立变量约束。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
图2. 离子液体的电导率—黏度关系、HOMO/LUMO 计算窗口与实测窗口偏差及聚类结果。DOI:10.1038/s41467-023-38493-7

图 2c、d 中,按阳离子和阴离子类型统计的计算窗口与数据库实测值平均绝对误差约为 1.0–1.1 V;图 2b 的电化学窗口与电导率未呈单调关系。HOMO/LUMO 适合划出明显淘汰区和待复核区,小于方法误差的名次差不宜写成稳定性优劣。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南

怎样由轨道排序估算氧化还原量?

计算成本最低的换算采用 Koopmans 型近似:电离势 IP 约为 −εHOMO,电子亲和能 EA 约为 −εLUMO。它忽略电子移走或加入后的轨道弛豫,对常见近似泛函还会混入自相互作用误差。这两个式子适用于同系列分子的相对预筛;相对于 Li+/Li 或 Na+/Na 的电极电位需要用带电态总能量和参比反应计算。

ΔSCF 以中性态与带电态的总能量差替代单个轨道本征值。保持中性构型不动得到垂直 IP 或垂直 EA,允许带电态弛豫得到绝热量;溶液氧化还原电位的换算包含溶剂化自由能、热修正、标准态修正和参比电极。垂直量接近电子转移发生瞬间,绝热量描述核坐标完成响应后的热力学差值,两者对应不同时间尺度。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
图3. 两类阴离子—溶剂对在 200 个构型上的垂直电离势分布及初末态能量分布。DOI:10.1038/s41467-019-11317-3

图 3 对 TFSI—PC 和 PF6—DME 各抽取 200 个分子动力学构型,垂直 IP 的展宽接近 1 eV。电解液在有限温度下会持续访问不同取向和距离,某个最低能优化结构给出的单点值无法覆盖这种波动。若研究氧化起始,分布低能侧的少量易氧化构型可能比均值更早贡献电流;候选总排序应报告均值、离散范围和采样数。

受控分子系列仍能体现轨道初筛的实际价值。图 4 在相同 EC/DEC 隐式介质和相同理论方法下比较 DME、呋喃、苯甲醚与 1,2-二氟苯。较高 LUMO 的呋喃和苯甲醚在 Li|Cu 循环后保留比例高于 LUMO 更低的 1,2-二氟苯,分子余量给出了轨道排序之外的化学证据。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
图4. DME、呋喃、苯甲醚与 1,2-二氟苯的 LUMO 能量、循环后余量及还原分解示意。DOI:10.1038/s41467-022-32192-5

这类相关性要求分子角色、介质和界面环境接近。若用同一阈值处理碳酸酯、离子液体、腈类和含氟稀释剂,方法误差、溶剂化差别和反应路径差别会共同影响轨道排序。同类分子可用 HOMO/LUMO 删除远离目标区间的对象,临界候选改用 ΔSCF 或完整反应自由能;名次处于不确定区间内时保留并列。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南

盐和局域构型怎样改变单分子排名?

盐溶解后,靠近电极的电子转移对象可能是自由溶剂、阴离子、接触离子对或更大的聚集体。此时应检查复合物的 HOMO/LUMO 位于哪个组分,并计算复合体系的电离或还原能。单分子表中的最高 HOMO 不一定属于实际先失电子的物种,配位和组分间静电稳定能会重新安排电荷落点。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
图5. 四种阴离子与 DMSO、DME、PC、ACN 构成分子对后的垂直电离势,以及孤立组分参照。DOI:10.1038/s41467-019-11317-3

图 5 中,TDI 较易失电子,四种溶剂对的 IP 主要随阴离子变化;换成 BF4 或 PF6 后,多数组合由溶剂失电子,分子对 IP 比孤立溶剂低。该研究用“阴离子 IP”与“溶剂 IP 减去约 2.8 eV”中的较小值拟合这组体系。2.8 eV 是该数据集中的电荷转移稳定项,不能移作所有盐和溶剂的固定修正。

TFSI—DMSO 对靠近时,溶剂失电子后形成的 AS+ 电荷转移态得到静电稳定;两者逐渐分开,稳定项按距离衰减,体系在约 7.2 Å 附近切换为阴离子失电子。一个孤立溶剂的 HOMO 数值无法包含这条距离依赖关系。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
图6. TFSI—DMSO 间距改变时的垂直电离势与氧化对象切换。DOI:10.1038/s41467-019-11317-3

构建复合物时不能只挑一个手工摆放姿态。分子动力学或构象搜索应覆盖自由溶剂、溶剂分离离子对、接触离子对和聚集体,再从各类构型抽样做量子化学计算。带电碎片的半局域泛函可能出现电子过度离域,看到 HOMO 跨多个远距离分子铺开时,需要用长程校正泛函、不同泛函对照或高精度波函数方法检查。还原侧也要检查阳离子配位后的空轨道位置,因为 Li+、Na+ 与羰基或醚氧配位会改变受电子态和断键反应。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南

哪些候选值得进入配方验证?

进入配方验证的候选应满足电子稳定区间和液体电解质的物性约束。主溶剂需要兼顾盐溶解、黏度、介电响应、供体能力和温区;共溶剂还要检查混溶性、是否进入第一溶剂化壳层,以及是否稀释局部离子网络;成膜分子则要追踪优先断裂的化学键、产物组成与电子绝缘能力。轨道阈值只占候选矩阵中的两列,其余列决定分子能否作为配方组分工作。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
图7. 非溶剂化共溶剂的锂离子溶剂化、混溶性、Lewis 碱度及循环库仑效率关系。DOI:10.1038/s41467-022-32192-5

图 7 比较了 14 种共溶剂。Kamlet–Taft Lewis 碱度 β 与(CIP+AGG)/SSIP 比值、Li|Cu 平均库仑效率的相关系数分别为 −0.887 和 −0.847,而气相 Li+—溶剂结合能对应的相关系数只有 0.233 和 0.273。这里的差异来自考察对象:单分子结合能描述一个局域配位,β 和归一化电子跃迁能 ETN包含混合环境中的相互作用与混溶信息。LUMO 再加入这两个量,才形成针对非溶剂化共溶剂的筛选坐标。

高电位正极上的脱氢、亲核进攻和表面催化过程涉及吸附态与过渡结构,低电位负极上的开环、脱氟和聚合则受表面电子态、缺陷和局域电场影响。临界候选应采用溶剂—盐—表面模型,比较反应自由能和必要的能垒,并用线性扫描伏安、浮充电流或产物分析核对。

如何通过HOMO与LUMO筛选溶剂分子?从计算到实际体系修正全指南
图8. 碳酸酯与磺酸酯电解液的 SEI 来源、NCM811 表面氧化稳定性及双描述符分区。DOI:10.1038/s41467-023-37999-4

图 8b 直接计算 EC、DMC、TM 和 TTMS 在充电态 NCM811 表面的氧化反应能;图 8c 又加入反应指数 RI 与氢转移能,区分溶剂成膜、溶剂/FSI混合成膜和阴离子成膜区域。TM、TTMS 在该体系中的位置由氧化反应与界面产物共同控制,数值分区只适用于论文采用的盐浓度、溶剂组合和正极模型。

候选进入实验配方前,可保留三档记录:同口径 HOMO/LUMO 明显越界者删除;轨道排序稳定、ΔSCF 与复合物结果方向一致者进入优先组;构象波动、盐配位或表面反应造成名次翻转者进入条件组。条件组要写明盐种、浓度、电极、温度和参比电位,避免把某一配方中的稳定区间移到另一体系。实验测量应对应计算目标:耐氧化看高电位电流与气体、耐还原看分子余量和 SEI 产物,传输问题另查电导率、黏度与迁移数。

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
(0)
上一篇 2小时前
下一篇 2025年7月28日 下午5:27

相关推荐