恢复、衰减与谱线宽度分开理解
重复太快会使不同信号恢复不充分,从而改变积分比例。
读图时保留这项条件:线宽还受场不均匀、交换及处理影响,不能全部归入同一弛豫时间。
同时看纵向恢复和横向衰减
改变两种特征时间,再移动时间指针,对照各自归一化响应。
4.4.1 纵向恢复与横向衰减各测什么
上一节用理想脉冲把平衡磁化转入横向平面。本节补上脉冲之后的时间过程:信号能持续多久,下一次激发又能得到多少磁化。先修是自由感应衰减的概念与旋转坐标系和翻转角。读完本节,应能区分三个弛豫参数,从指数衰减推得线宽,并估算重复激发给两个信号的积分比带来的偏差。
展开说明与推导
弛豫(relaxation)描述偏离平衡的自旋体系趋向平衡的过程。静磁场取沿 z 轴,热平衡时纵向磁化为 \(M_0\),横向净磁化为零。纵向弛豫时间(longitudinal relaxation time) \(T_1\)刻画 \(M_z\) 返回 \(M_0\) 的速度,常称自旋—晶格弛豫时间。这里的晶格是环境的历史名称;溶液中的分子运动也能参与能量交换,并不要求样品是固体。
横向弛豫时间(transverse relaxation time) \(T_2\)刻画横向相干磁化的衰减,常称自旋—自旋弛豫时间。横向分量按相位相加;各部分相位逐渐分散时,净分量可减小,而纵向布居尚未恢复到平衡。纵向恢复与横向衰减可以同时发生,不能把横向信号的消失直接解释为全部磁化已经返回 z 轴。Keeler(2002)的弛豫讲义以能级布居与横向相干分别描述这两件事。
本页的 \(T_1\)、\(T_2\)是单一组分、固定条件下的指数时间常数,单位均为秒;速率 \(R_1=1/T_1\)、\(R_2=1/T_2\)的单位为每秒。两者描述不同分量,不普遍相等。表观横向弛豫时间(effective transverse relaxation time) \(T_2^*\)由 FID 的表观包络得到,还包含场不均匀等去相位贡献,必须另行标记。
展开说明与推导
同一分子不同氢环境也可能有不同 \(T_1\)。弛豫速率取决于局部相互作用与分子运动,溶剂、温度和顺磁物质等条件都可能改变它。需要定量时,应在目标样品条件下测量或验证相关信号的恢复,不能把某个官能团的经验时间当作固定常数。本节不给出供所有样品套用的弛豫参数表。
4.4.2 Bloch 方程与两种指数解
Bloch 方程(Bloch equations)把宏观磁化的进动与返回平衡的速率合在一起。用矢量 \(\mathbf M=(M_x,M_y,M_z)\)表示磁化,\(\mathbf B\)表示磁场,\(\gamma\)表示氢核旋磁比,实验室坐标系中的现象学形式为:
展开说明与推导
磁化可用安培每米表示,磁场用特斯拉表示,\(\gamma\)常以弧度每秒每特斯拉表示。叉乘项只描述转动,后两项分别使横向分量衰减、纵向偏差消失。Bloch(1946)提出了这类带有限弛豫时间的宏观描述;它在本页作为模型使用,尚未计算每种分子相互作用给出的速率。
射频关闭后,在绕 z 轴旋转的坐标系中取有符号偏置 \(\Omega\),并约定横向相位随时间变化为 \(-\Omega t\)。忽略交换、交叉弛豫与耦合调制,\(T_1\)、\(T_2\)、\(M_0\)在考察期间为常量,分量方程成为:
先令 \(u=M_z-M_0\)。方程给出 \(\mathrm du/u=-\mathrm dt/T_1\),积分后 \(u(t)=u(0)e^{-t/T_1}\)。再令 \(M_+=M_x+\mathrm iM_y\),两个横向方程合为 \(\mathrm dM_+/\mathrm dt=-(1/T_2+\mathrm i\Omega)M_+\)。因此:
若理想 90° 脉冲把纵向分量变为零,恢复曲线为 \(M_z/M_0=1-e^{-t/T_1}\),一个 \(T_1\)后约恢复 63.2%;横向包络在一个 \(T_2\)后剩下约 36.8%。时间常数给的是指数变化的尺度,不能读成某时刻突然完成恢复。理想 180° 反转则给 \(M_z(0)=-M_0\),从而 \(M_z/M_0=1-2e^{-t/T_1}\),过零时间为 \(T_1\ln 2\)。改变初值会改变曲线,方程中的恢复速率仍是同一个 \(1/T_1\)。
习题 4.4-1 · 从初值计算恢复
设 \(T_1=2.00\,\mathrm s\)。分别从 \(M_z(0)=0\)和 \(M_z(0)=-M_0\)出发,计算两秒后的 \(M_z/M_0\),并求反转曲线的过零时间。
参考解答两秒等于一个 \(T_1\),所以从零恢复为 \(1-e^{-1}=0.63212\);从负平衡值恢复为 \(1-2e^{-1}=0.26424\)。过零条件为 \(1-2e^{-t/T_1}=0\),解得 \(t=2.00\ln 2=1.38629\,\mathrm s\)。两种初值使用同一个时间常数,却给出不同磁化,不能只报时间而遗漏制备方式。
4.4.3 场不均匀怎样进入 T₂*
实际线圈同时接收样品多个位置的信号。即使这些位置含有同一种氢,略微不同的静磁场也使它们有不同偏置。脉冲后各位置起初相位相近,经过时间 \(t\)便积累不同相位;线圈接收的是复数分量之和。这个总和减小所需的时间,可以短于各位置本身的横向弛豫时间。匀场改善的是空间磁场分布,不能据此认为分子的固有 \(T_2\)必然改变。
展开说明与推导
把共同的化学共振频率提出,设场不均匀给出的额外频率分布为 \(p(\delta\nu)\),以赫兹为自变量,并归一到积分为一。忽略位置间扩散且各位置有同一 \(T_2\),归一化 FID 包络的复数形式可写成:
\(G(t)\)是静态偏置分布带来的去相位因子。只有它也能写成单指数 \(e^{-t/T_{\mathrm{inh}}}\)时,乘积才给出 \(1/T_2^*=1/T_2+1/T_{\mathrm{inh}}\)。这是速率相加的特例。若偏置呈高斯(Gaussian)分布,\(G(t)\)含与 \(t^2\)有关的指数项;其他分布还可能产生振荡。此时单个 \(T_2^*\)可能只是指定拟合范围内的表观参数,应同时给出包络模型,不能强迫所有 FID 都服从同一指数。
自旋回波(spin echo)可以检验静态去相位与不可重聚衰减的差别。在理想 \(90^\circ-\tau-180^\circ-\tau\)序列中,各位置在第一段积累的相位差经 180° 脉冲镜像,第二段仍按原偏置进动,在 \(2\tau\)时重聚。Hahn(1950)的原始论文讨论了这种回波及弛豫测量。回波幅度仍受 \(e^{-2\tau/T_2}\)衰减;时间不变的频率偏置能够抵消,随机的相干损失不能被同一操作全部恢复。
本图忽略扩散、交换、耦合演化和脉冲误差。真实回波衰减可能含这些贡献,因而回波拟合值也需核查模型;不能把任意回波序列拟合出的一个时间直接称为纯粹 \(T_2\)。FID 较短、回波较长支持静态去相位的存在,但不能单独确定分子结构。
4.4.4 从指数 FID 推得 Lorentz 线宽
半高全宽(full width at half maximum, FWHM)是谱线达到峰高一半时两侧频率的间隔。本节只推导独立、正确相位的吸收型单线。设其复数 FID 为 \(S(t)=S_0e^{-Rt}e^{-\mathrm i2\pi\nu_0t}\),\(R=1/T\),\(T\)可为符合单指数模型的 \(T_2\)或 \(T_2^*\)。假定从零采集到无穷、无额外窗函数,以正号 Fourier 核变换得到:
展开说明与推导
实部 \(A(\nu)\)是 Lorentz 线形(Lorentzian line shape),峰高为 \(S_0/R\)。半高条件要求分母变为 \(2R^2\),故两侧偏移各为 \(R/(2\pi)\),合起来得到:
教学设定 \(T_2=0.400\,\mathrm s\)时,固有单线宽度约为 0.796 Hz;若额外指数去相位使 \(T_2^*=0.100\,\mathrm s\),该 Lorentz 模型的观测宽度约为 3.183 Hz。不能把后者倒算出的 0.100 s 当作分子的 \(T_2\)。非指数包络对应的非 Lorentz 谱线没有这一通用换算,非对称场分布甚至可使线形明显偏斜。
线宽还受到有限采集、数字采样、窗函数与相位等处理条件影响。一个多重峰的总跨度包含耦合裂分,不能代入单条谱线的宽度公式。对于式(4.4-5)的理想全线积分,面积为 \(S_0/2\),不随 \(R\)改变;宽峰低一些并不意味着氢数少。真实积分窗有限,展宽后更多尾部可能落在窗外,或者与相邻信号重叠,面积便会偏离这一理想结果。相关处理条件将在4.5至4.7中展开。
习题 4.4-2 · 线宽能够说明什么
某独立吸收线的 FWHM 为 2.00 Hz。若已经验证单指数 FID、足够采集和无额外展宽,求对应时间常数。若另一条线呈高斯形,可否按同式求 \(T_2^*\)?若该信号是跨度 10 Hz 的三重峰,可否用跨度求 \(T_2\)?
参考解答Lorentz 条件下 \(T=1/(\pi\times2.00)=0.15915\,\mathrm s\)。尚未排除场不均匀时应称表观 \(T_2^*\),不直接称 \(T_2\)。高斯线形对应不同包络,不能套用单指数换算;需写出其线形或时域模型。三重峰总跨度主要包括裂分位置的差异,不能当成每条谱线的 FWHM。应独立量取或拟合单线,并检查重叠。
4.4.5 重复时间、饱和与两信号积分比
重复时间(repetition time, TR)是相邻两次对应激发之间的完整周期。简单单脉冲实验中,它通常包括弛豫延迟 d1、采集时间 AQ 以及其他序列时间,不能只读参数 d1。Oxford NMR Facility(2025)的 T₁ 讲义明确提醒重复速率通常由 d1 与 AQ 之和决定。采集横向信号的同时,纵向磁化也在恢复;这段时间仍属于两次激发的间隔。
展开说明与推导
设每次激发均为理想 90°,前一周期留下的横向磁化在下次脉冲前已消失或被破坏,脉冲长度相对弛豫时间可忽略。每次脉冲后 \(M_z\)归零,下一次脉冲前便恢复到 \(M_0[1-e^{-TR/T_1}]\)。在相同激发与线性接收条件下,完整积分与该磁化成比例。重复过快导致恢复不足,称为饱和(saturation);不同信号恢复程度不同时会出现差异饱和。
\(f_i\)为无量纲恢复因子,\(I_i^{(0)}\)为充分恢复时的积分,相同线性接收带来的公共比例因子在比值中消去。这里比较的是信号积分;若要比较不同分子的物质的量,还必须除以各信号代表的氢数。公式假定没有选择性激发差异、重叠、基线误差或磁化转移;增加扫描次数只会把同一种偏差重复累加,不能使 \(f_A/f_B\)自动回到一。
教学设定真积分比 \(r_0=3.000\),\(T_{1,A}=1.000\,\mathrm s\),\(T_{1,B}=4.000\,\mathrm s\),完整 \(TR=2.000\,\mathrm s\)。先算指数比:\(TR/T_{1,A}=2\),\(TR/T_{1,B}=0.5\);再算 \(f_A=0.864665\)、\(f_B=0.393469\)。因此 \(I_A:I_B=(3\times0.864665):(1\times0.393469)\),观测比为 6.592621,相对真比值偏高约 119.754%。长 \(T_1\)的 B 信号受抑制更多。把 \(TR\)改为 20 s,A 近于完全恢复,B 恢复 0.993262,观测比仍约为 3.020351,偏高约 0.678%。这些数值来自模型计算,未代表仪器测量。
展开说明与推导
最长相关 \(T_1\)的五倍是常用设置起点,但 \(1-e^{-5}=0.993262\),即约 99.326% 恢复,并非精确的完全恢复。若只要求某信号的剩余恢复缺口不超过 \(\varepsilon\),模型给出 \(TR\geq -T_1\ln\varepsilon\);例如缺口不超过 0.1% 要求约 \(6.908T_1\)。这个单信号要求不能自动充当整个定量方法的误差保证,积分比也仍取决于两个恢复因子的比。
教学工具:恢复因子与观测积分比
输入 A、B 的 T₁(秒)、完整 TR(秒),以及充分恢复时的真积分比 A/B(无量纲)。本页采用理想 90°、残余横向磁化已消失或被破坏、单指数恢复、相同线性接收的稳态模型;TR 包括采集等整个周期。可用范围:T₁ 为 0.001–10000 s,TR 为 0.001–100000 s,真比值为 0.000001–1000000。超出范围或空白时不计算。
恢复因子与比值均无量纲;显示值经过舍入,接近 100% 时不表示数学上完全恢复。工具只显示当前输入的理论预测,不作为真实 qNMR 自动校正或精度保证。
本工具适用于已进入重复周期的理想 90° 单脉冲模型。首次由完全平衡开始的扫描可以具有不同强度。更复杂翻转角下纵向磁化未必归零,残余横向磁化、交换与重叠也会改变结果;本页不以计算器自动校正真实 qNMR。Oxford 的 qNMR 讲义(2017)还要求核查激发、匀场、信噪与处理条件,相关定量验证见7.1至7.3。
习题 4.4-3 · 检查积分偏差和完整周期
两信号真比为 \(3:1\),\(T_{1,A}=1\,\mathrm s\),\(T_{1,B}=4\,\mathrm s\)。简单实验设 d1=1 s、AQ=1 s,其他时间忽略。按本页模型求观测比。若两信号 \(T_1\)都改为 2 s、其余不变,比值又是多少?解释为什么准确比值不表示完全恢复。
参考解答完整 \(TR=d1+AQ=2\,\mathrm s\),恢复因子为 \(1-e^{-2}=0.864665\)和 \(1-e^{-0.5}=0.393469\),观测比为 \(3\times0.864665/0.393469=6.592621\)。误用 d1 会给另一个错误预测。若 \(T_1\)均为 2 s,则两恢复因子均为 \(1-e^{-1}=0.632121\),在比值中抵消,观测比仍为 3。每个信号强度均只有充分恢复时的约 63.2%,所以相同程度的饱和能保留比值,却降低单次信号。
4.4.6 把弛豫条件写进证据和实验决策
从结构解析角度,弱峰与少氢不能直接画等号。若一个候选需要某组氢而观察不到清楚信号,应先检查它是否被差异饱和压低,是否因展宽而混入基线,或被另一信号覆盖。增加 TR 后积分比趋于稳定支持恢复不足曾影响旧结果;只有在相位、基线、区域与样品条件可比时,这个比较才有解释力。相同温度下延长周期还应注意样品是否随时间反应或降解。
展开说明与推导
线宽异常时先看影响范围。多个无关信号都呈相似偏斜,通常优先检查匀场、样品均匀性及处理;局部信号展宽则需要考虑交换或该位点的相干衰减,但单凭外观不能确定机制。理想静态去相位与 \(T_1\)恢复是不同问题:延长 TR 能改变饱和,却通常不能修复磁场分布造成的线形。匀场改善峰形,也不能保证不同信号已经充分纵向恢复。
需要估计 \(T_1\)时,可按反转恢复原理改变等待时间,用末端 90° 脉冲把当时的 \(M_z\)转成可观察信号,再拟合恢复曲线;需保留负信号的相位和符号。Oxford 的 2025 年讲义给出实验与处理流程,但实际参数应由本地训练和样品条件决定。这里说明测量原理,并未在本页运行仪器或软件。拟合残差若随时间系统变化,应先检查单指数模型和脉冲条件,而非只抄一个拟合常数。
记录一份用于比较积分的谱图时,至少保留脉冲角、完整 TR 的组成、目标信号 \(T_1\)的来源及温度,并保存原始 FID、采集长度和加窗设置。已知线形时注明 FWHM 的测量对象是独立谱线,已知积分受重叠限制时注明区域。本页的恢复计算能够指出需要延长周期的方向;可靠含量测量还要在第 7 章加入内标选择、信号分离和不确定度记录。
关键术语
- 纵向弛豫时间 (T₁)
- 纵向磁化偏离平衡值的单指数恢复时间常数。
- 横向弛豫时间 (T₂)
- 均匀条件下横向相干磁化的衰减时间常数。
- 表观横向时间 (T₂*)
- 包含场不均匀等去相位贡献的表观 FID 时间尺度;非指数时须注明定义。
- Bloch 方程 (Bloch equations)
- 用进动与现象学弛豫项描述宏观磁化演化的方程。
- 去相位 (dephasing)
- 各信号贡献的相位分散使它们的矢量和减小。
- 自旋回波 (spin echo)
- 脉冲使可重聚的相位差在指定时刻抵消而形成的信号。
- 半高全宽 (FWHM)
- 同一条谱线峰高一半处两侧频率之差,通常以 Hz 表示。
- Lorentz 线形 (Lorentzian line shape)
- 单指数 FID 的正确相位吸收型 Fourier 线形。
- 重复时间 (TR)
- 相邻对应激发之间的完整周期,包含采集等序列时间。
- 饱和 (saturation)
- 重复激发使纵向磁化未充分恢复而造成的信号抑制。
参考文献与延伸阅读
- Bloch, F. 1946. Nuclear Induction. Physical Review 70, 460. 原始论文 DOI。宏观进动与有限弛豫时间描述的来源。
- Hahn, E. L. 1950. Spin Echoes. Physical Review 80, 580. 原始论文 DOI。静态频率分布、回波与扩散影响;链接不等同于授权转载全文。
- Keeler, James. 2002. Understanding NMR Spectroscopy, Chapter 8: Relaxation; Chapter 9: Phase Cycling and Gradient Pulses, §9.2.4 Lineshapes. University of Cambridge. 弛豫讲义;线形推导。本文重新绘图和推导,不转载讲义插图。
- University of Oxford NMR Facility. 2025. Measuring relaxation times — Inversion recovery for T₁. 官方讲义。重复周期与反转恢复测量;平台操作需结合本地规程。
- University of Oxford NMR Facility. 2017. Quantitative NMR Spectroscopy. 官方讲义。充分弛豫、采集时间与定量采集及处理条件。