理想共晶预测有多深？面向多组分混合物中理想低共熔溶剂（DES）预测的数值框架

《Journal of Chemical & Engineering Data》：How Deep is Your DES? A Numerical Framework for Ideal Eutectic Prediction in Multicomponent Mixtures

　　

研究人员在此提出一种简单且计算高效的流程，用于确定理想二元、三元及更高阶混合物的共晶温度，为评估低共熔溶剂（DES）行为提供定量参考。在DES文献中，三元体系常通过固定两组分摩尔比、仅沿固定组成线（等压线）改变第三组分而简化为拟二元体系。尽管这减少了实验工作量，却将液相线（liquidus）映射限制于单一组成路径，可能遗漏全局理想共晶。相比之下，该研究实现在假设拉乌尔定律（Raoult-law）理想性及纯组分、互不相溶固相的前提下，对所有组分求解经典理想液相线方程。利用纯组分熔化温度 Ti* 与摩尔熔化焓 ΔHfus,i，该方法可实现多维液相线曲面映射。根求解实现通过对解析对称混合物解、独立80位精度计算、MATLAB实现、置换检验、3521个压力测试体系及已发表基准进行验证。所有压力测试体系均收敛，已发表理想共晶温度复现误差在0.7 K以内。应用于文献报道的二元与三元混合物进一步表明，与全局理想共晶比较可将表观局部热偏差与非负理想组成惩罚相结合。该框架因而为多组分共晶筛选提供了透明且可扩展的理想平衡基线。

研究背景与问题提出

低共熔溶剂（Deep Eutectic Solvent，DES）作为可调液体介质在成本、毒性与生物降解性方面具备潜力，其热力学定义要求共晶温度低于对应理想液体参考值。然而DES文献中三元体系常被简化为拟二元等压线研究，仅沿固定组成路径映射液相线，无法保证识别全局理想共晶。同时，针对三元及更高阶体系的完整液相线映射与理想共晶求解在DES领域仍较稀缺，金属合金领域虽经CALPHAD发展，但DES所需的是一种透明、廉价、可复现的拉乌尔定律理想参考。为此，研究人员开展了面向二元至八元组分的理想共晶数值框架构建，以纯组分 T_i^* 与 ΔH_fus,i 为输入，求解全局理想共晶温度 T_e^id 与组成，并提供多维液相线场可视化。

主要关键技术方法

研究人员采用以下关键技术：基于化学势平衡的固液平衡方程，令活度系数 γ_i=1 导出理想液相线方程 ln x_i = ?(ΔH_fus,i/R)(1/T ? 1/T_i^*)；将共晶条件写为归一化方程 Σ exp[?(ΔH_fus,i/R)(1/T_e^id ? 1/T_i^*)] = 1，构造单调函数 G(T) 并用数学保证的根区间 [T_low, T_high] 包围唯一根；使用Python中 scipy.optimize.brentq（xtol=1e-12 K，rtol=4ε）与独立MATLAB fzero 交叉验证；三元采用重心晶格 n=80（3321点），四元四面体晶格 n=40（12341点）作图形离散，但共晶解独立于网格最小值；利用对称解析解、mpmath 80位二分、置换检验、3521随机与极值体系作数值鲁棒性验证。无生物样本队列。

研究结果

3.1 数值验证与基准测试

通过对称N组分解析解（x_i,e=1/N，T_e^id=[1/T^*+(R/ΔH_fus)ln N]^?1）校验，双精度下温度与组成误差为零。非对称三元6种置换均得 T_e^id=328.7065610691162 K，迭代数一致。3521个随机体系（2–8组分，T_i^*∈[200,900] K，ΔH_fus,i∈[2,60] kJ mol^?1）全部收敛，最大归一残差2.84e-14，与80位参考最大温差2.98e-13 K。外部基准Teixeira等奎宁体系复现误差≤0.7 K。Python与MATLAB最大温差3.13e-13 K。M5机二进制中位25.28 μs，八元组41.26 μs，t_μs=20.04+2.64N（R²=0.9987）。网格细化显示 n=80 三元RMSE 0.076 K，但网格最小节点比真根高0.775 K，证明共晶须由方程求根而非网格取最小。

3.2 理想多组分液相线图构建

以假想三元（T^*=400/450/500 K，ΔH_fus=5/7/9 kJ mol^?1）绘线框、着色面与等边三角等高投影三视图。二元共晶谷 T_ab=300 K、T_ac=322 K、T_bc=352 K 收敛至三元共晶 T_abc=277 K，证明新增组分经混合熵进一步压制熔点，且全局共晶不能由单一等压线获得。

讨论与结论翻译

讨论指出，文献报道转变温度 T_r 与理想液相线 T_m,χ 及全局理想共晶 T_e^id 之差（ΔT_χ=T_r?T_m,χ，ΔT_global=T_r?T_e^id）可将表观局部热偏差拆解为理想组成惩罚与非理想活度系数贡献；表观活度系数 γ_i,app 由单点反推，可能含多态、分解、动力学过冷等干扰，故不单独判正负偏差。结论为：该框架以纯组分 T_i^* 与 ΔH_fus,i 为唯一输入，在拉乌尔定律理想、固相纯且互不相溶近似下，可对任意2–8组分体系在<0.1 ms内求解全局理想共晶温度与组成，3521体系全收敛，跨语言一致至1e-13 K级，并为实验DES是否“深于理想”提供透明基线，远离理想最低温区的配方可快速剔除，需非理想处理者据此触发COSMO-RS等精细计算。论文发表于《Journal of Chemical & Engineering Data》。