第4章 · 4.4

4.4 弛豫、线宽与饱和

Relaxation, Linewidth and Saturation
本页范围与学习目标
纵向磁化恢复、横向相干衰减和静态场不均匀分别影响信号的强度与线形。本节从 Bloch 方程推导指数恢复与衰减,说明 Lorentz 线宽公式的条件,并用两个信号的完整重复周期计算差异饱和造成的积分偏差。
先读图解和核心解释;需要细节时再展开。
先看图 · 阅读入口
打开本节交互图

恢复、衰减与谱线宽度分开理解

重复太快会使不同信号恢复不充分,从而改变积分比例。

恢复、衰减与谱线宽度分开理解 比较 1 纵向恢复 沿静磁场恢复 影响重复采集响应 比较 2 横向衰减 相干响应随时间减弱 与线形相关 比较 3 场不均匀 增加表观失相 需要区分本征弛豫 恢复、衰减与谱线宽度分开理解 比较 1 纵向恢复 沿静磁场恢复 影响重复采集响应 比较 2 横向衰减 相干响应随时间减弱 与线形相关 比较 3 场不均匀 增加表观失相 需要区分本征弛豫
图 4.4-1 本页概念关系图。三个并列框分别表示不同的对象或判断,按各自条件使用。

读图时保留这项条件:线宽还受场不均匀、交换及处理影响,不能全部归入同一弛豫时间。

两个原创指数模型恢复和衰减,是两个不同方向的过程教学设定:两种特征时间都为 1 s纵向:逐渐恢复1005 s横向:逐渐衰减1005 s纵轴分别按对应初值或平衡值归一化
图 4.4-2 两个原创指数模型。纵向从零恢复到平衡值,横向从初值衰减;各自归一化,不表示两个实验信号在同一绝对尺度上。实际弛豫时间、初始条件与场不均匀还需分别核查。
动手看一看 · 教学模型

同时看纵向恢复和横向衰减

改变两种特征时间,再移动时间指针,对照各自归一化响应。

同时看纵向恢复和横向衰减两种响应各自归一化纵向恢复横向剩余
恢复与衰减
恢复与衰减归一响应00.250.50.75102.557.510时间 / s
当前纵向恢复63.21%
当前横向剩余1.832%
本征模型 FWHM0.6366 Hz
2 s
0.5 s
2 s
同时看纵向恢复和横向衰减两种响应各自归一化纵向恢复横向剩余

图 4.4-3 独立指数模型;线宽换算只对应均匀场中的本征横向衰减,未加入交换或场不均匀。

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\)表示氢核旋磁比,实验室坐标系中的现象学形式为:

展开说明与推导
\[\frac{\mathrm d\mathbf M}{\mathrm dt}=\gamma\mathbf M\times\mathbf B-\frac{M_x\hat{\mathbf x}+M_y\hat{\mathbf y}}{T_2}-\frac{M_z-M_0}{T_1}\hat{\mathbf z}.\]
(4.4-1)

磁化可用安培每米表示,磁场用特斯拉表示,\(\gamma\)常以弧度每秒每特斯拉表示。叉乘项只描述转动,后两项分别使横向分量衰减、纵向偏差消失。Bloch(1946)提出了这类带有限弛豫时间的宏观描述;它在本页作为模型使用,尚未计算每种分子相互作用给出的速率。

射频关闭后,在绕 z 轴旋转的坐标系中取有符号偏置 \(\Omega\),并约定横向相位随时间变化为 \(-\Omega t\)。忽略交换、交叉弛豫与耦合调制,\(T_1\)、\(T_2\)、\(M_0\)在考察期间为常量,分量方程成为:

\[\frac{\mathrm dM_z}{\mathrm dt}=-\frac{M_z-M_0}{T_1},\qquad \frac{\mathrm dM_x}{\mathrm dt}=\Omega M_y-\frac{M_x}{T_2},\qquad \frac{\mathrm dM_y}{\mathrm dt}=-\Omega M_x-\frac{M_y}{T_2}.\]
(4.4-2)

先令 \(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_+\)。因此:

\[M_z(t)=M_0+[M_z(0)-M_0]e^{-t/T_1},\qquad M_+(t)=M_+(0)e^{-t/T_2}e^{-\mathrm i\Omega t}.\]
(4.4-3)

若理想 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\)。

以各自弛豫时间归一化的恢复与衰减 左图从零开始的纵向磁化在一个T1后恢复到平衡值的百分之63.2;右图横向磁化在一个T2后衰减到初始值的百分之36.8。两图使用不同的归一化时间轴,不能读成T1等于T2。 纵向恢复:初值为零横向衰减:初值归一为 1 Mz / M0|Mxy| / |Mxy(0)| 1010 63.2% 恢复36.8% 残余 012345 012345 时间 t / T1时间 t / T2 示意:每幅图按各自的时间常数缩放,两条时间轴不可直接比较秒数。
图 4.4-4 单指数模型的两种时间过程。曲线由本页公式绘制;纵向磁化的恢复对象是平衡值,横向磁化的衰减对象是相干分量。教学示意不表示同一样品的实测参数。
习题 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 包络的复数形式可写成:

\[F(t)=e^{-t/T_2}G(t),\qquad G(t)=\int_{-\infty}^{\infty}p(\delta\nu)e^{-\mathrm i2\pi\delta\nu t}\,\mathrm d\delta\nu.\]
(4.4-4)

\(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}\)衰减;时间不变的频率偏置能够抵消,随机的相干损失不能被同一操作全部恢复。

静态偏置可重聚,横向弛豫仍然损失幅度 四个时刻分别为90度脉冲后、自由演化一段时间、180度脉冲后,以及再演化同一段时间。静态不同偏置的矢量先散开,180度脉冲使相位镜像,之后重聚;最后矢量缩短表示T2弛豫不能被回波恢复。 90° 后演化 τ理想 180° 后再演化 τ 相位一致不同偏置累积相位差相位关于 x 轴镜像相位重聚 t = 0τ−τ+2τ 示意:忽略耦合、扩散与交换;末端长度按 T2 衰减,重聚不等于恢复全部幅度。
图 4.4-5 静态场不均匀与自旋回波的区别。中间脉冲改变相位的位置,偏置频率保持原值;两段等长演化抵消静态偏置造成的相位差。箭头为示意,不表示单个质子的经典轨迹。
回波测量的边界

本图忽略扩散、交换、耦合演化和脉冲误差。真实回波衰减可能含这些贡献,因而回波拟合值也需核查模型;不能把任意回波序列拟合出的一个时间直接称为纯粹 \(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 核变换得到:

展开说明与推导
\[\widetilde S(\nu)=\int_0^\infty S(t)e^{\mathrm i2\pi\nu t}\,\mathrm dt=\frac{S_0}{R-\mathrm i2\pi(\nu-\nu_0)},\qquad A(\nu)=\frac{S_0R}{R^2+[2\pi(\nu-\nu_0)]^2}.\]
(4.4-5)

实部 \(A(\nu)\)是 Lorentz 线形(Lorentzian line shape),峰高为 \(S_0/R\)。半高条件要求分母变为 \(2R^2\),故两侧偏移各为 \(R/(2\pi)\),合起来得到:

\[\Delta\nu_{1/2}=\frac{R}{\pi}=\frac{1}{\pi T}\quad(\mathrm{Hz}).\]
(4.4-6)

教学设定 \(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=1-e^{-TR/T_{1,i}},\qquad I_i=I_i^{(0)}f_i,\qquad r_{\mathrm{obs}}=\frac{I_A}{I_B}=r_0\frac{f_A}{f_B},\quad r_0=\frac{I_A^{(0)}}{I_B^{(0)}}.\]
(4.4-7)

\(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)
重复激发使纵向磁化未充分恢复而造成的信号抑制。

参考文献与延伸阅读

  1. Bloch, F. 1946. Nuclear Induction. Physical Review 70, 460. 原始论文 DOI。宏观进动与有限弛豫时间描述的来源。
  2. Hahn, E. L. 1950. Spin Echoes. Physical Review 80, 580. 原始论文 DOI。静态频率分布、回波与扩散影响;链接不等同于授权转载全文。
  3. Keeler, James. 2002. Understanding NMR Spectroscopy, Chapter 8: Relaxation; Chapter 9: Phase Cycling and Gradient Pulses, §9.2.4 Lineshapes. University of Cambridge. 弛豫讲义;线形推导。本文重新绘图和推导,不转载讲义插图。
  4. University of Oxford NMR Facility. 2025. Measuring relaxation times — Inversion recovery for T₁. 官方讲义。重复周期与反转恢复测量;平台操作需结合本地规程。
  5. University of Oxford NMR Facility. 2017. Quantitative NMR Spectroscopy. 官方讲义。充分弛豫、采集时间与定量采集及处理条件。