分离混合物先检查模型能否区分组分
拟合得到一组比例后,还要检验是否存在同样好的其他分配。
读图时保留这项条件:相同或十分相近的参考谱会限制可辨识性,低残差不能独自解决歧义。
让比例变化与可辨识性同时可见
改变峰间距和已知构造比例;在完全重合时比较不同配比。
7.5.1 混合谱叠加需要先定义响应
混合样品把多个组分的信号放在同一个谱轴上。若各组分相互作用可以忽略,测量响应保持线性,且混合后各自的位移、耦合和线形与所用参照相容,可以把谱写成若干参照谱的加权叠加。权重与组分量之间的关系,仍取决于参照的归一方式和定量条件。一个拟合系数不能凭其名字自动变成摩尔分数。
展开说明与推导
设在离散谱点或完整积分区域 \(k\) 上,观测为 \(y_k\)。组分 \(j\) 的参照响应为 \(s_{kj}\),按明确量尺度定义的系数为 \(c_j\),背景为 \(b_k\),未被模型解释的变化为 \(\varepsilon_k\)。一个基本线性模型为:
\(m\) 是纳入的候选组分数,\(k\) 表示所用观测。若参照为每单位浓度的响应,\(c_j\) 可对应浓度;若所有参照先被任意归一到相同总面积,则系数主要描述对应信号贡献,还需氢数及响应换算。背景与残差也不同,不能把每一项未解释信号都事先定义为噪声。
本节使用小型教学区域模型,让读者先理解可辨识性。Dashti 等(2017)的自旋系统建模工作和 GISSMO 可以提供结构化谱参数;BMRB 的小分子实验条目提供条件明确的参照资料。采用参照时须保持溶剂、温度、酸碱状态与场强相容,混合物中的状态变化还可能需要新的模型。
在指定数据、模型和误差范围内,不同组分量能否得到可区分的预测,称为当前问题的可辨识性。软件找到一个解,不能证明其他解都已被排除;能检测混合信号,也不能保证各个贡献可分别确定。
展开说明与推导
可辨识性有结构与实际两层。无噪声模型中存在多组完全相同预测,属于模型本身的信息不足;理论上唯一但受噪声、重叠和参照误差影响很大,则是实际估计不稳定。增加扫描次数主要降低部分噪声,对两组完全相同预测无能为力。需要先辨认限制来自哪里,才知道应提高信噪比还是增加新的区分观测。
组分数也不能由峰数直接读出。一个分子可有很多氢环境和多重峰,多种分子又可能共享重叠区域;快速交换还可能把同一物质的状态平均。分析前区分化合物种类与当前可见化学状态。如果混合后发生结合或酸碱变化,纯组分参照可能不再适用,需要明确状态模型和共同物料关系,而非任意增加一列未知谱来改善拟合。
7.5.2 用三个区域检验两个组分
教学设定:在相容定量条件下,A 每单位量在区域一、二、三的响应为 \((1,0,1)\),B 为 \((0,1,1)\)。区域一只报告 A,区域二只报告 B,区域三同时报告二者。观测为 \((0.70,0.30,1.00)\)。前两行直接给 \(c_A=0.70,c_B=0.30\),第三行预测 \(c_A+c_B=1.00\),与观测相容。
| 区域 | A 单位响应 | B 单位响应 | 混合观测 |
|---|---|---|---|
| 一 | 1 | 0 | 0.70 |
| 二 | 0 | 1 | 0.30 |
| 三 | 1 | 1 | 1.00 |
展开说明与推导
第三行在这个模型中是前两行之和,所以并没有提供另一种区分 A、B 的方向。它仍有检查模型和量关系的价值:若第三区域观测为 1.20,模型会留下 0.20 的偏差。原因可能包括未知组分、背景、积分不完整或响应不相容。不能直接把这项偏差分到 A、B 上,再同时宣称前两区域的直接测量仍被精确满足。
若采用全部谱点,很多点属于同一条谱线或共同峰形,相邻点之间可能相关。采样点数很大,并不等于拥有同样多条独立化学约束。增加数字点密度可以改善表示与拟合,但不能把一个完全重叠区域变成多个独有峰。区分能力取决于组分响应的差异,误差分析还要考虑数据的相关性。
实际参照可能包含不同数量的可见氢,因此使用同一归一方式前要检查其含义。若 A 的参照总面积代表六氢、B 代表十二氢,把二者总面积各归一到一后获得相同信号份额,并不表示相同摩尔量。应保存参照的量尺度和报告氢数,把拟合输出换算为所需比例,而非在最终一步才猜测系数的单位。
用固定积分窗口时,峰移动可使一部分面积越过边界,看起来像组分量变化。可先核对完整峰组和参照,再选择适当对齐、较宽区间或具有物理依据的峰模型。对齐也不能无限灵活:若把每个局部区域独立拉到参照位置,可能抹掉真实化学状态变化,甚至使错误组分得到吻合。处理选择需要保留原始坐标和依据。
习题 7.5-1 · 总量区域的作用
沿用表7.5-1,前两区域不变,第三区域为1.20。能否把系数改成0.80与0.40,使第三行精确吻合?
参考解答这样会把前两行分别预测成0.80与0.40,均偏离原来观测。应整体比较误差并检查背景、遗漏组分与响应条件。第三行是模型相容性检查,不能单独决定系数后忽略另外两行。若有噪声,可作加权拟合,但必须保留各行残差。
7.5.3 完全相同的响应允许多组组成
再设 A、B 在所有可用区域的单位响应均为 \((1,2,0)\),观测也是 \((1,2,0)\)。第一行给 \(c_A+c_B=1\),第二行只把相同等式乘以二,第三行没有响应。于是 \((0.20,0.80)\)、\((0.50,0.50)\)、\((0.90,0.10)\) 都给完全相同观测。非负约束把范围限制到零至一,仍不产生唯一组成。
展开说明与推导
把这种响应排列成表或矩阵时,每列表示一个组分,每行表示一个观测。独立列的数量称为列秩。在刚才的模型中两列相同,只有一个独立响应方向,不能分别确定两个系数。读者无需先掌握完整矩阵代数,也可直接通过三行等式检查:所有信息都只测量了二者之和。
在本节完全相同响应模型中,设总量为一。拖动 A 的百分比,B 自动补足;显示的三个区域响应由模型计算。
展开说明与推导
算法若输出其中一组比例,可能来自初始值、优化规则或额外先验。先验可以合理使用,例如已知组分不存在或某个量已由独立方法测得,但需说明它来自哪里。不能把程序选择的一组可行解当成谱图独自证明的事实。强制系数之和为一也只是增加一个量定义或约束,未必能消除同方向响应的歧义。
增加另一条只报告 A 的观测后,模型就可能得到新的区分方向。若该观测支持 A 为0.30,原来总量为一便给 B为0.70。新的信息可以来自独有一维峰、合适的异核相关、选择性实验或其他测量。若 A、B 仅为在非手性条件下不可区分的对映体,则通常需要手性依据,单纯增加常规非手性谱的分辨率不足以解决。
习题 7.5-2 · 相同响应再加约束
两组分响应完全相同,总量为一。某独立测量只给 A 在0.20至0.40之间。能否得到唯一组成?
参考解答只能将范围限制为 \(c_A\in[0.20,0.40]\),\(c_B=1-c_A\in[0.60,0.80]\),仍有多组相同谱预测。应报告这个共同约束范围。独立测量若给可靠精确值及不确定度,才可进一步收窄,而其不确定度须传播到B。
7.5.4 近似相同的响应会放大误差
实际问题更多是组分谱相似而非完全相同。教学 A、B 在两区域的单位响应分别为 \((1,0.9)\) 与 \((1,1.1)\),观测 \((1,1)\) 给 \(c_A+c_B=1\) 及 \(0.9c_A+1.1c_B=1\),解为各0.50。第二区域若变为1.02而第一不变,就得到 \(c_B=(1.02-0.9)/0.2=0.60\),A为0.40。
展开说明与推导
第二区域仅改变0.02,系数却各改变0.10;当前区分来自两个响应之间只有0.2的差。当参照误差或基线变化与这项差异相当时,比例估计会不稳定。这种问题常称为病态或条件不佳。拟合残差仍可能很小,因为多个比例变化能互相补偿;小残差不单独证明每个系数准确。
可以通过已知组成、重复制样、合理改变处理和参照条件,检查估计范围。参考谱若允许各组分独立位移、线宽、强度及背景同时自由变化,模型会更灵活,也可能使组分之间相互替代。应限制参数到有依据的范围,保留残差与参数变化,并检查新增自由度是否得到独立数据支持。
还可保留没有参与拟合的报告区域检验预测。若一个系数集合只解释被选择的两个重叠区,却不能复现同组分的独有峰,它的身份或量尺度仍有问题。留出区域要具有可比响应并排除额外污染,不能随意挑最有利结果。已知比例样品和未用于调参数的样品,可帮助判断方法是否只贴合了当前资料。
库中没有的真实组分可能被错误分配到已有组分,或进入背景与未解释信号。先检查结构化残差和候选范围。允许噪声项并不意味着所有被划入该项的信号都来自随机噪声。
展开说明与推导
Domżał 等的 Magnetstein 研究使用已知或预期组分参照谱,采用最优传输相关方法处理位移、线形和分辨率差异。它提供了另一种量估计思路,输入仍需要混合谱与组分库;没有纳入的结构不会自动变成完整指认。方法对谱差异的容忍与组分区分能力需要共同检验,容忍范围过大也可能使相似候选难以分开。
比较方法时,应使用同样明确的任务与数据划分:已知组分量估计、有限库检索和未知组分结构解析具有不同目标。对纯¹H输入有效的结果,不应扩写成所有复杂混合物都能从头解决。研究中的数据和代码可用于独立复核,只有实际运行相应流程并留下输出后,才能报告当前环境复现表现。
7.5.5 用新的观测改善区分能力
当限制来自相似响应时,优先寻找对组分变化具有不同敏感性的观测。改变溶剂或酸碱条件可能分开位移,但也可能改变样品状态;二维异核谱可利用碳分散,选择性方法可减少某些重叠,扩散实验可提供运动维度。任何增强区分的方法都应说明适用域,并检查新增资料与原有编号和量模型是否相容。
展开说明与推导
独有峰也不是永久属性。某区域在当前候选库中只含 A,真实样品的未知物可能仍贡献相同位置。通过多个独立报告核、加标或相关网络检查,可以提高归属可信度;仍然无法排除的贡献则进入结果范围。目标是说明数据实际能够确定哪些量,避免用复杂拟合掩盖未解决的组成歧义。
保存参照谱条件、归一方式、候选库、处理参数和逐区域残差。说明系数对应信号份额、浓度还是摩尔比例。多个组成都在误差范围内相容时,报告范围或候选,列出可提供新区分方向的观测。
展开说明与推导
下一节讨论 DOSY 与扩散。它可以为某些混合物增加区分维度,仍然受重叠、交换、运动模型和拟合条件限制。先理解本节的组成歧义,才能判断一个漂亮的二维扩散图究竟增加了多少可用信息。
习题 7.5-3 · 检查敏感性
沿用近似相同响应模型,第二区域为0.98时两系数是多少?这组变化对报告精度有什么提示?
参考解答\(c_B=(0.98-0.9)/0.2=0.40\),A为0.60。第二区域正负0.02变化使系数在0.40与0.60之间改变,说明当前区分较敏感。需要参照和测量误差依据来确定比例范围,不能只报软件显示的多位小数。
关键术语
- 谱叠加 (spectral superposition)
- 在适用线性条件下组合各组分响应。
- 参照谱 (reference spectrum)
- 具有明确条件与量尺度的组分响应资料。
- 可辨识性 (identifiability)
- 模型和数据能否区分不同组成参数。
- 列秩 (column rank)
- 响应排列中独立组分方向的数量。
- 病态问题 (ill-conditioned problem)
- 输入小变化引起较大参数变化的问题。
- 残差 (residual)
- 观测与模型预测之间未解释的差。
- 非负约束 (nonnegative constraint)
- 要求组分量系数不小于零的限制。
- 先验信息 (prior information)
- 当前谱以外、具有明确来源的附加知识。
参考文献与延伸阅读
- Dashti H, Westler WM, Tonelli M, Wedell JR, Markley JL, Eghbalnia HR. 2017. Spin System Modeling of Nuclear Magnetic Resonance Spectra for Applications in Metabolomics and Small Molecule Screening. Analytical Chemistry, 89: 12201–12208. 自旋模型与跨场强响应。
- Domżał B, Nawrocka EK, Gołowicz D, Ciach MA, Miasojedow B, Kazimierczuk K, Gambin A. 2023 在线发表. Magnetstein: An Open-Source Tool for Quantitative NMR Mixture Analysis Robust to Low Resolution, Distorted Lineshapes, and Peak Shifts. Analytical Chemistry. 混合物定量原始研究。
- BMRB. Metabolomics. 条件明确的小分子实验资料。
- Claridge TDW. 2016. High-Resolution NMR Techniques in Organic Chemistry, 3rd edition. Elsevier. 混合物、相关与选择性实验。
- Odell B, Claridge T; University of Oxford NMR Facility. NMR Handouts. 实验选择与处理背景。