基团贡献法预测咪唑类离子液体的环境足迹：LCA-GC耦合模型构建及应用

本文创新性地将生命周期评估（LCA）与线性基团贡献（GC）模型结合，首次构建了可用于预测咪唑类离子液体（ILs）生产过程全球变暖潜势（GWP）、人类毒性（HTP）及生态毒性等环境影响的预测模型。该研究为在ILs设计早期阶段快速评估其环境可持续性性能提供了概念验证，使研究者能在功能属性筛选的同时兼顾环境影响，推动了可持续化学设计，并支持可持续发展目标（SDGs）9、12和13的实现。

　　

简介

离子液体（ILs）是一种熔点通常低于100°C的有机盐，凭借其可调的物理化学性质（如低蒸汽压、高热稳定性、良好的导电性和宽电化学窗口），为可持续创新提供了无限可能。它们被视为常规溶剂、催化剂、电解质以及润滑油、涂料等配方产品中添加剂的潜在“绿色”替代品，在生物质转化、气体捕获与储存、能源存储与转换等诸多领域具有应用前景，并已出现工业化应用案例，如BASF的BASIL?技术和Chevron/Honeywell UOP的ISOALKY?技术。尽管目前约有500种离子液体实现不同程度的商业化生产，且存在庞大的阳离子-阴离子组合可能，但其生命周期环境可持续性评估研究仍不充分，主要受限于数据稀缺和系统复杂性导致传统LCA研究耗时费力。

现有的预测方法，如定量构效关系（QSPR/QSAR）、分子模拟、分子基状态方程和基团贡献（GC）模型，多用于预测离子液体的物理化学性质或固有毒性，而忽略了其生产过程中涉及的前体及其他生命周期阶段的上游环境负担。同时，离子液体“绿色溶剂”的印象主要源于其可忽略的挥发性，但对其合成、可回收性和处置等其他可持续性指标的评估同样重要。研究表明，离子液体的合成通常涉及多个步骤和众多前体，与常规溶剂相比可能产生更高的环境影响。然而，鉴于离子液体组合的庞大规模，量化数百种离子液体的环境影响是一项艰巨的任务，需要专业知识和大量资源，使得传统的LCA方法在其设计早期阶段难以实施。

方法学

本研究采用的方法学包含四个主要步骤：第一步是离子液体样本选择。研究选取了咪唑类阳离子作为研究对象，因其是被研究最广泛、数据最易获得的离子液体家族。样本包含30种不同的咪唑基离子液体，考虑了阳离子烃基侧链（如烷基、芳基、多侧链和链长变化）和三种常用阴离子（氯离子Cl^-、四氟硼酸根[BF₄]^-和六氟磷酸根[PF₆]^-）的变体，以确保化学基团在不同分子间重复出现，促进GC模型参数的可迁移性。

第二步是合成路线选择与主要合成步骤识别。研究假设所有咪唑基离子液体通过烷基化和复分解反应合成，并将合成路线分解为四个可区分的步骤：咪唑环形成、甲基化、烷基化和阴离子复分解。其中，咪唑环形成通过乙二醛、甲醛（或乙醛，用于合成1-丁基-2,3-二甲基咪唑鎓盐）和氨的Debus–Radziszewski反应实现；甲基化步骤是咪唑环与氯甲烷反应生成1-甲基咪唑；烷基化步骤是1-甲基咪唑与不同的烷基氯反应生成氯化咪唑鎓盐；最后，氯化咪唑鎓盐通过与NaBF₄或LiPF₆进行阴离子复分解反应，分别得到四氟硼酸盐或六氟磷酸盐离子液体。净化与干燥步骤因方法多样且缺乏细分能耗数据而未纳入本研究的系统边界。

第三步是离子液体的环境影响评估。采用生命周期评估（LCA）方法，依据ISO 14040和14044标准，评估了30种咪唑基离子液体的环境影响。系统边界为“摇篮到大门”，功能单位设定为“每摩尔产品”以在分子水平评估累积贡献，并补充“每千克产品”以方便文献比较。生命周期清单数据基于上述合成步骤的化学计量计算和能量消耗估算获得，背景数据来自Ecoinvent数据库v3.8。使用CML-IA 2001: August 2016方法评估了全球变暖潜势（GWP）、人类毒性潜势（HTP）、淡水水生生态毒性潜势（FAETP）、海水水生生态毒性潜势（MAETP）和陆地生态毒性潜势（TETP）这几类环境影响。

第四步是基团贡献分析与模型验证。基于LCA估算的环境影响数据，采用线性基团贡献模型，将环境影响与每个合成步骤中引入的功能化学基团相关联。模型形式为环境影响等于各化学基团贡献值乘以其在分子中出现次数的总和。共定义了10个不同的化学基团。从30种离子液体中随机选取24种用于模型训练和验证，其余6种作为测试集。为减少过拟合，采用了k折交叉验证和自助抽样法进行重采样。模型拟合时，咪唑基团的贡献是单独计算的，等于环形成和甲基化步骤贡献的平均值加上线性模型的截距，不参与优化。模型预测准确性最终通过测试集进行评估。

结果与讨论

离子液体的环境影响评估

全球变暖潜势

对于咪唑基氯化物离子液体，其全球变暖潜势估计在每摩尔0.73至1.07千克二氧化碳当量之间。咪唑环形成步骤是主要热点，贡献了总影响的44%至63%，这主要归因于乙二醛的使用。甲基化步骤的贡献约为每摩尔0.19千克二氧化碳当量，占总影响的17%至37%。除了[Hmim]Cl（仅一个取代基），烷基化步骤是第二大热点，贡献了21%至39%的总影响。侧链的变化显著影响环境影响，例如，具有最长烷基链（[C₁₀mim]Cl）或芳基（[Phmim]Cl, [Bzmim]Cl）的离子液体显示出最高的全球变暖潜势。取代基数量对全球变暖潜势的影响似乎很小，而侧链中是否存在双键影响也不大。

当将氯离子替换为[BF₄]^-和[PF₆]^-时，所得离子液体的全球变暖潜势分别增加了约2.7倍和4.8倍。这主要是因为复分解步骤所需的钠盐或锂盐（NaBF₄, LiPF₆）本身具有很高的碳足迹。阴离子的选择对离子液体的总体全球变暖潜势有巨大影响。

人类毒性影响

所有咪唑基氯化物离子液体的人类毒性潜势估计约为每摩尔6千克DCB当量。咪唑环形成步骤占主导地位，贡献了超过98%的总影响，这主要是由于乙二醛生产过程中释放的环氧乙烷排放物毒性极高。甲基化和烷基化步骤的贡献可忽略不计。对于咪唑基四氟硼酸盐，环形成步骤仍是主要热点（69%），而对于咪唑基六氟磷酸盐，使用LiPF₆盐贡献了约52%的影响，成为主要热点。与全球变暖潜势类似，[PF₆]^-阴离子由于其前体盐（LiPF₆）生产过程中释放的氢氟酸排放，导致人类毒性潜势比相应的氯化物对应物高出约2倍。

生态毒性影响

对于咪唑基氯化物，淡水水生生态毒性潜势和海水水生生态毒性潜势分别估计在每摩尔0.2-0.3千克DCB当量和0.4-0.6吨DCB当量。环形成步骤再次占主导地位，分别贡献了54-95%和50-87%的总影响。烷基化步骤贡献显著（淡水水生生态毒性潜势 29-53%），贡献大小取决于烷基链长度。甲基化步骤由于副产盐酸的系统信用，贡献可忽略甚至为负。陆地生态毒性潜势结果在每摩尔3.85至6.85克DCB当量之间，环形成步骤贡献52-74%，烷基化步骤贡献28-51%。

对于含[BF₄]^-和[PF₆]^-的离子液体，阴离子复分解反应因其前体盐的生态毒性而成为所有生态毒性影响类别的绝对主导因素。与相应的氯化物相比，四氟硼酸盐离子液体的陆地生态毒性潜势、淡水水生生态毒性潜势和海水水生生态毒性潜势分别高出约3倍、5倍和19倍；而六氟磷酸盐离子液体则分别高出约7倍、9倍和116倍。再次表明[PF₆]^-阴离子在所有环境影响类别中均导致更高的负担。

与其它离子液体及常规有机溶剂的比较

将本研究中[Hmim]Cl、[C₄mim]Cl和[C₄mim][BF₄]的全球变暖潜势结果（以每千克计）与文献报道的其他离子液体以及一组常规有机溶剂进行了比较。比较显示，不同研究、不同合成路线（特别是卤化物路线与无卤路线）得出的结果存在显著差异。总体而言，本研究的全球变暖潜势估值处于文献报道的类似离子液体的范围之内或偏高。更重要的是，本研究中评估的离子液体的全球变暖潜势显著高于所比较的常规溶剂（如己烷、甲醇）的平均水平，凸显了在评估离子液体“绿色”声称时考虑其生产阶段环境影响的重要性。

基团贡献预测模型结果与讨论

研究成功开发了用于预测咪唑基离子液体全球变暖潜势、人类毒性潜势和淡水水生生态毒性潜势的线性基团贡献模型。模型训练和验证显示了良好的性能。对于全球变暖潜势，模型预测与LCA计算值高度一致，测试集的平均绝对百分比误差较低，表明模型具有良好的预测能力。化学基团贡献分析揭示了不同基团对环境影响的定量贡献。例如，在全局变暖潜势模型中，甲基、亚甲基等饱和烷基链基团的贡献为正值且随链长增加，而芳基（C_b-(C)C_b-(N)）和乙烯基（C_v-(C)）等基团也显示出特定的贡献值。模型成功捕捉了阴离子的巨大影响，[BF₄]^-和[PF₆]^-基团的贡献值远高于其他基团，与LCA结果一致。这些模型参数使得仅根据离子液体的分子结构即可快速估算其生产的环境影响。

研究意义与展望

本研究首次将基团贡献法应用于预测离子液体的环境影响，扩展了离子液体性能预测工具包，将全球变暖等环境可持续性指标纳入其中。这为研究人员和制造商在离子液体设计早期阶段，同步筛选功能属性和环境性能提供了实用工具，是迈向其在净零技术图景中商业化应用的重要一步。该方法有助于将定量环境指标嵌入可访问的预测框架，推动可持续创新。

然而，本研究也存在一些局限性，例如样本仅限于咪唑类阳离子和三种阴离子，合成路线假设较为简化，且未包括净化和干燥步骤。未来工作建议扩大样本范围以涵盖其他阳离子头基（如铵、磷），比较替代合成路线（如无卤路线），并纳入净化步骤以扩展分析。尽管存在这些局限，本研究提出的基团贡献-生命周期评估耦合框架为高效、高通量评估离子液体的环境足迹奠定了坚实基础，有助于实现真正可持续的化学品设计。