多元线性回归创新应用于二价金属离子与阴离子表面活性剂相互作用分析

《Analytical Science Advances》：Innovative Application of Multiple Linear Regression for Analysis of Interactions of Divalent Metal Ions With Anionic Surfactant

本综述创新性地应用多元线性回归（MLR）方法，深入分析了十四种二价金属离子（Me2+）与十二烷基硫酸盐（DS）三聚体头基的相互作用。研究通过质谱（MS）数据构建模型，揭示了金属离子结合强度与其水结构性质（如熵贡献Sstr、ΔGHB）间的线性关联，为理解表面活性剂在溶液中的行为及其在环境修复、材料科学等领域的应用提供了新的化学计量学视角。

　　

引言

离子表面活性剂，特别是阴离子表面活性剂，因其在食品、工业、农业、材料和药理学等领域的广泛应用而成为科学研究中最常使用的表面活性剂。其在溶液中的性质和聚集行为会受到无机离子的显著影响。多价金属阳离子与带负电的表面活性剂头基之间存在强烈的静电相互作用，这种作用远强于一价阳离子，因此能更显著地改变阴离子表面活性剂的性质。选择合适的含多价离子的表面活性剂体系在工业、农业和材料科学中具有重要意义。近年来，多价阳离子与阴离子表面活性剂之间的相互作用在环境保护方面也受到关注，特别是在检测和去除废水中的多价金属（如有毒重金属镉、铅、汞、镍、锌、铜）方面显示出应用潜力。

在这些体系中，水是普遍存在的介质。溶质和离子会影响水的结构，基于此效应，溶质被分为结构制造者（cosmotropic）或结构破坏者（chaotropic）。本研究所涉及的金属离子均属于水结构制造者。深入了解多价金属离子与阴离子表面活性剂的相互作用，对于相关应用至关重要。本研究旨在通过创新的化学计量学方法，分析十四种二价金属离子与十二烷基硫酸盐（DS）三聚体的相互作用。

材料与方法

试剂与化学品

研究所用电解质包括：BeCl₂、MgCl₂、CaCl₂、SrCl₂、Ba(NO₃)₂、MnCl₂、FeCl₂、CoCl₂、NiCl₂、CuCl₂、ZnCl₂、CdCl₂、Hg(NO₃)₂、PbCl₂，以及表面活性剂十二烷基硫酸钠（SDS）。所有化学品均未经进一步纯化直接使用。所研究的水溶液中含有浓度均为1 × 10^?4M的表面活性剂和给定盐。

实验方法——串联质谱（MS/MS）测量

实验数据通过串联质谱测量获得。将含有特定电解质的水性表面活性剂溶液注入配备电喷雾离子源的Q-ToF质谱仪。样品在300°C下去除水分，随后复合离子[Me(DS)₃]^?发生裂解。在负离子模式下，碰撞能量（E_C）从0 eV增加到59 eV（步长为1 eV）采集MS/MS谱图。从每个碰撞能量下的质谱图中读取前体离子（即[Me(DS)₃]^?）和碎片离子（即DS^?）的峰强度（I），并根据公式Y = I_precur/ (I_precur+ I_fragm)计算分子离子存活产率（Y）。Y值既是[Me(DS)₃]^?离子稳定性的量度，也反映了Me²⁺与DS三聚体头基的结合强度。

数据分析

数据分析包括数据预处理和多元线性回归（MLR）分析。

首先进行数据预处理，包括平滑强度-碰撞能量曲线以去除噪声，并根据平滑后的曲线计算特定碰撞能量下的离子峰强度和Y值。随后对数据进行标准化处理。

本研究应用多元线性回归（MLR）分析来估算对于给定的Me²⁺阳离子，前体离子和碎片离子对离子存活产率的贡献。所使用的MLR方程为：Y = b₀+ b_precur* I_precur+ b_fragm* I_fragm。其中，b_precur和b_fragm参数分别代表前体离子和碎片离子对Y值的贡献份额。分析所用数据集为Y与E_C呈线性下降关系的中间碰撞能量范围内的I_precur、I_fragm和Y值。同时计算了决定系数R²以评估模型拟合优度。

结果与讨论

MLR分析结果（见表1与图1）显示，对于所有研究的金属离子，b_precur为正值，b_fragm为负值，其绝对值均小于1，且│b_precur│ + │b_fragm│之和近似等于1。这与理论预期相符，即b_i值代表各离子对Y的贡献份额，且Y值随前体离子峰强度增加而增加，随碎片离子峰强度增加而减少。所有模型的R²值均高于0.95（大部分高于0.99），证明MLR模型能很好地描述实验数据。

研究的核心目标是发现Me²⁺离子的水结构性质与质谱数据之间的线性关联。如图2所示，研究人员绘制了金属离子的水结构对熵（S_str）和热容（C_p,str）的贡献、其表面张力增量（k_i）及其ΔG_HB值与MLR得到的b_precur值之间的关系图。ΔG_HB被定义为离子从重水转移到轻水时，水分子参与的平均氢键数量的变化，是离子水结构效应的定量指标，ΔG_HB正值越大的离子是更好的水结构制造者。

分析发现：

1. 1.
   就S_str和ΔG_HB与b_precur的关系而言（图2a, d），数据点分布在两条明确的回归线上，从而将所研究的结构制造阳离子分为了两个不同的组。一条回归线（r²= 0.91）贯穿Mg、Cd、Zn、Mn、Ca、Sr、Ba和Pb离子，另一条（r²= 0.85）则贯穿Ni、Fe、Co、Cu和Hg离子。这表明S_str和ΔG_HB在离子分类上具有相似的物理化学意义。此外，在这两个离子组中，b_precur值均随Me²⁺离子半径的增大而增加。
2. 2.
   就C_p,str与b_precur的关系而言（图2b），10种金属离子的数据点形成一条线（金属组1），而Be²⁺、Co²⁺、Cu²⁺和Hg²⁺的数据点则明显偏离该线（金属组2）。
3. 3.
   此项分类与已有文献部分吻合。例如，已有研究将Ca、Mg、Mn分为一组，将Ni²⁺和Co²⁺分为另一组。本研究中，Ca、Mg、Mn同属一组（图2a, d），而Ni和Co同属另一组，结果一致。

结论

本研究成功应用多元线性回归（MLR）分析了14种二价金属离子与DS三聚体头基的相互作用，重点分析了作为复合离子稳定性及金属结合强度量度的分子离子存活产率（Y）。MLR模型参数b_precur和b_fragm有效量化了前体离子和碎片离子对Y的贡献，高决定系数（R²> 0.95）证明了模型对实验数据的良好描述。

创新性地应用化学计量学分析后，成功发现了Me²⁺离子的水结构性质（特别是熵贡献S_str和氢键自由能变化ΔG_HB）与质谱结果之间的线性相关关系，并将这些水结构制造阳离子分为两个明确的组别，此分类与文献报道部分一致。

本研究提出的MLR创新应用，为质谱数据提供了更深入、更合理的物理解释。这种数学分析方法有助于进一步理解二价金属离子与阴离子表面活性剂的相互作用，从而推动表面活性剂在环境修复、材料科学等多个领域的应用。