第9章 · 9.4

9.4 深入脉冲实验理论

Density Operators and Coherence Transfer
本页范围与学习目标
用可逐步核验的理想二自旋计算建立密度矩阵、乘积算符与相干转移的联系,说明相位循环和梯度如何选择实验通路,以及这些理论对二维谱解释的条件。
先读图解和核心解释;需要细节时再展开。
先看图 · 阅读入口

脉冲序列沿状态变化逐步检查

复杂序列可以拆成状态、演化和检测三步,逐步核对符号与条件。

脉冲序列沿状态变化逐步检查 步骤 1 规定初始状态 核与算符的约定 保持坐标和相位一致 步骤 2 逐步演化 脉冲、耦合与位移 分别改变相应状态 步骤 3 选择并检测 相位循环与梯度 筛选可观察通路 脉冲序列沿状态变化逐步检查 步骤 1 规定初始状态 核与算符的约定 保持坐标和相位一致 步骤 2 逐步演化 脉冲、耦合与位移 分别改变相应状态 步骤 3 选择并检测 相位循环与梯度 筛选可观察通路
图 9.4-1 本页概念关系图。箭头表示本页采用的检查顺序。

读图时保留这项条件:理想算符通路不是可直接输入仪器的完整程序,还需实际校准与控制。

9.4.1 从一个矢量到自旋体系的状态

一维脉冲的矢量图能解释磁化如何转入横向平面,但二维谱还需要追踪两个核之间的相关状态。一个核的箭头不能同时记录耦合伙伴处于不同状态时产生的相反相位。本节用密度算符和乘积算符描述这一过程,并说明实验为什么能选择特定相关通路。

展开说明与推导

密度矩阵描述样品中许多分子的统计状态。它在所选基底中的对角元素是布居,非对角元素记录相干;完整矩阵必须厄米、半正定且迹为一。相干的具体矩阵元素依赖基底,因此应说明采用什么能级或自旋基底,不能把任意非对角数都直接称为可接收的信号。

高温液体核磁常把状态写成单位矩阵部分加一个很小的偏差。单位矩阵经理想旋转仍不变,对无迹磁化算符的期望值也为零;分析脉冲时可以只追踪偏差,但这不意味着样品已经完全极化。教材中的 Iz 是一条初始偏差分量,并非归一化后的全部密度矩阵。

(9.4-1)\[\rho(t)=U(t)\rho(0)U^\dagger(t),\qquad U(t)=\exp[-i\mathcal H t],\qquad \langle A\rangle=\operatorname{Tr}[\rho A].\]
展开说明与推导

这里使用无量纲自旋算符,ℋ 为能量哈密顿量除以 ℏ,单位是 rad·s⁻¹;t 用秒,指数为无量纲。自由演化、射频脉冲和梯度都改变传播算符。公式先忽略弛豫与交换,这些过程需要额外的动力学描述,不能仅靠一个幺正矩阵模拟全部损失。

两核体系可按双方的高低能自旋状态排列为四个基底状态,每个算符因此是四行四列。单核算符作用在自己的核上,对另一核保留单位矩阵;双核相关项才同时作用于两者。计算时固定这组排列,脉冲、状态和接收算符使用同一排列,才能比较矩阵系数。

还可用矩阵内积判断状态包含多少目标项:将目标算符与当前状态相乘取迹,再除以目标算符自身的平方迹。正交基中的不同项互不贡献,这使传播结果能与乘积算符推导逐项比较。只比较整份矩阵的某一个元素,容易漏掉其他残留项。

定义 · 状态与可观测量

状态描述体系储存的全部信息;接收线圈只检测其中与横向磁化相关的部分。实验可以先把信息存入暂时不可直接接收的相关状态,再用后续脉冲转为可检测磁化。

9.4.2 乘积算符给相关状态命名

对于两个自旋二分之一的核 I 和 S,六个单核算符 Ix、Iy、Iz、Sx、Sy、Sz 加上九个双核乘积和单位矩阵,可以展开任意四阶密度矩阵。实际推导只追踪当前存在的几项,不必每一步都把十六个矩阵元素写开。

展开说明与推导

Ix 表示 I 多重峰的同相横向磁化。2IxSz 表示 I 磁化与 S 的纵向状态相关:S 两种状态对应的 I 跃迁有相反符号,通常称反相磁化。它仍可在频率分辨谱中表现为正负分裂线,只是简单积分可能相互抵消。

同时含两个横向因子的项需要进一步分解。例如 IxSx 可展开为零量子和双量子分量。只看到两个横向字母就给整个实算符指定一个相干级数会出错。采用升降算符 I±=Ix±iIy 后,各纯相干分量的级数更容易识别。

弱耦合异核模型的耦合部分是 2πJIzSz,J 用 Hz。由此 Ix 随时间变为同相与反相的组合。πJt 出现在以下正弦余弦中;不能因哈密顿量里有 2π,就把这里的角度再乘二。不同自旋算符的归一化会影响书写系数,需始终使用同一约定。

(9.4-2)\[I_x\longrightarrow I_x\cos(\pi Jt)+2I_yS_z\sin(\pi Jt).\]
展开说明与推导

当 t=1/(2J) 时,同相项为零,得到纯反相项。这是理想模型的转化时间,不是任意多自旋样品的最佳延迟。若一个 I 同时耦合多个 S,演化会产生更多乘积项;若同核频差很小,则还应保留强耦合哈密顿量,简单二项规则不再充分。

9.4.3 一个可以逐步核验的转移例子

取 J=140 Hz 的理想异核二自旋体系,从 Iz 分量开始。规定正旋转为 U=exp(−iθIα),忽略化学位移演化、弛豫、有限脉宽和其他耦合。先施加 I 核绕 y 轴的九十度脉冲,得到 Ix;再耦合演化 3.571 ms,得到 2IySz。

展开说明与推导

接着对 I 施加 x 轴九十度,对 S 施加 y 轴九十度。前者把 Iy 转为 Iz,后者把 Sz 转为 Sx,于是相关项变成 2IzSx。再经历相同的耦合时间,得到 Sy。原先属于 I 的偏差贡献已经转为 S 的可接收横向磁化。

本书用四阶矩阵实际传播上述五个状态,并与预期乘积算符逐个比较。系数均为一,最大矩阵元素差小于 2×10⁻¹⁶,属于浮点舍入量级。这核验了所用符号和理想转移链;没有采集仪器数据,也没有验证实际样品的转移效率。

若真实初始状态还含有 S 的热平衡磁化,应把这一分量单独经过同样序列,然后与 I 来源的结果相加。它可能形成不希望的直接信号,需要实验选择抑制。线性传播允许分开计算各项,但最终接收包含所有未被选择消除的贡献。

图 9.4-2 五个实际矩阵核验的状态 IzIx2IySz2IzSxSy I: 90°y耦合 3.571 msI: 90°xS: 90°y耦合 3.571 msJ = 140 Hz;理想瞬时脉冲与耦合演化
图 9.4-2 教学转移链采用本节的主动旋转约定。它展示 INEPT 类转移的算符逻辑,并非可直接输入仪器的完整脉冲程序。
展开说明与推导

实际实验通常用回波把化学位移演化重聚焦。在两核同时施加理想一百八十度脉冲时,两者的纵向符号都翻转,位移项反号,而乘积 IzSz 保持不变,因此可以保留耦合演化。若只反转一个核,耦合项也反号,结果不同。

仪器序列常把总耦合演化分成回波两侧的延迟,因此图上一个延迟的名称未必等于本例的总时间。读方法部分时应加上各段持续时间,并确认有限脉冲期间是否计入有效演化。把单侧延迟与总延迟混同,会使本应转为反相的时刻落在另一状态。

在 HSQC 中,磁化经过异核单量子相干的间接演化,再转回质子检测;在 HMBC 中,较长转移延迟偏重较小的长程耦合。算符分析能解释这些选择怎样发生,也提醒我们:转移效率随 J、脉冲和弛豫变化,相关峰强度不等于所连质子数。

方法 · 读脉冲序列

先写出初始偏差项、每段作用的哈密顿量、脉冲旋转轴及接收算符,再逐段传播。核对希望储存和检测的项是否出现,最后加入相干选择、弛豫与仪器条件。

9.4.4 相干级数与相位循环

相干级数 p 根据共同 z 轴旋转时的响应定义。令 Fz 为所讨论同种核纵向算符之和,纯分量在角度 φ 下乘 exp(−ipφ)。纵向磁化的 p 为零,I+ 为正一,I− 为负一,I+S+ 为正二;异核梯度处理还需分别记录各核的贡献。

(9.4-3)\[e^{-i\phi F_z}\rho_p e^{i\phi F_z}=e^{-ip\phi}\rho_p.\]
展开说明与推导

相位循环改变某些脉冲相位并配合接收相位,使目标通路在多次扫描相加时同相,其他通路抵消。必须先找到能让目标和干扰获得不同相位的步骤;若两条通路对所循环脉冲完全相同,增加扫描数也无法将它们分离。

考察一个抽象四步相位滤波:目标分量 p=1,φ 依次取零、π/2、π、3π/2,接收端补偿目标相位。各分量留下的平均系数是 exp[−i(p−1)φ] 的四项平均。p=1 得一;p=0、−1、2 得零;p=5 也得一,因为相位差每步都是整周。

图 9.4-3 离散相位滤波的选择与混叠 p = −101251101补偿目标 p = 1;保留系数的模每四级重复:仍须考虑所有实际可达通路
图 9.4-3 本节实际计算的抽象相位滤波。五个级数的结果显示有限步数具有周期性选择;是否存在 p=5 取决于实际自旋体系。
展开说明与推导

这个平均不是一份具体 HSQC 的完整相位表。真实相位循环按脉冲处相干级数的变化设计,还可能同时处理轴峰、脉冲误差和正交检测。部分完成一个规定循环会破坏相消条件;只从最终谱看不见干扰,也不能证明所有非目标通路已被抑制。

二维数据的间接维还需要区分正负频率,常用规定的相位递增或成对采集实现。相干选择表与间接维正交检测都涉及相位,却回答不同问题。修改相位循环时若同时破坏了频率符号编码,二维峰可能镜像或产生不正确的色散分量,不能靠简单转置坐标修复。

9.4.5 梯度把通路写成空间相位

z 方向梯度使共振频率随位置变化。对同核纯相干分量,一段理想梯度产生与 pγz∫Gdt 成比例的空间相位。样品不同位置的贡献相加时可能互相抵消;后续梯度若补偿目标通路的累计相位,该通路重聚焦,而不满足条件的通路仍被分散。

(9.4-4)\[\Phi(z)=-z\sum_k p_k\gamma\!\int_kG(t)\,dt,\qquad \sum_kp_k\!\int_kG(t)\,dt=0\quad\text{(同核重聚焦条件)}.\]
展开说明与推导

γ 采用 rad·s⁻¹·T⁻¹,梯度用 T·m⁻¹,空间位置用米。异核通路需把各核的 pγ 分别求和,不能仅按整数级数比设梯度比。有限脉宽、梯度形状、涡流和仪器规定会影响实际实现,以上等式只是静止样品的理想相位条件。

例如同核通路在两段相等面积梯度中依次具有 p=1 与 p=−1,空间相位可以抵消;若两段都是 p=1,则同向等面积梯度不能重聚焦它。改变脉冲造成的相干级数,与改变梯度方向是两种不同操作,读序列时应同时记账。

梯度还会对运动敏感。即使静止位置上的相位完全补偿,分子扩散或对流使其经历不同位置,也可能造成信号损失。这与第7章 DOSY 的扩散编码相关,但相干选择实验中的损失不能直接当成经过校准的扩散系数;需要相应序列和拟合模型。

9.4.6 理想通路与实验交叉峰之间

理想传播回答某种相干能否由脉冲链产生,实验还要回答它能否在噪声、重叠和损失下被识别。有限带宽使不同位移的翻转角不完全一致,脉冲误差会留下不期望的项;转移延迟过长时横向弛豫可能削弱目标,交换也可能改变相关状态。

展开说明与推导

缺失交叉峰因此有多种原因。应先看该核是否处于激发带宽内,耦合是否与延迟匹配,峰是否因相位或处理相消,再判断结构证据。反之,残留的非目标通路、折叠或重叠也可能造成貌似相关的峰,不能仅凭一幅二维图就指定完整连接。

当一个异常峰只在某些相位步出现时,保留逐步扫描或适当控制数据有助于判断它的来源。若只保存最终相加后的谱,原本可以区分的响应已经混在一起。对需要比较转移效率的研究,宜同时记录目标与控制条件的原始数据,而不是只保存裁剪后的正峰区域。

脉冲理论的实际用途包括设置合理延迟、识别可观察项及规划补充实验。报告应保留序列名称和版本、相位与梯度选择、延迟、核种及处理方式。若进一步使用 Spinach 等工具,需明确哈密顿量、初始状态、弛豫模型和接收算符,计算输出才与所问实验相对应。

误判 · 理论允许不等于实测出现

一个非零转移系数说明理想模型存在该路径;观察到的峰仍需实验识别,未观察到的峰也需评估检测和转移条件。算符计算不能替代数据质量检查。

9.4.7 习题与解答

习题9.4-1 耦合时间

按式9.4-2,J=8 Hz 时,纯反相时刻与同相项第一次回到反号的时刻是多少?

解答纯反相在 t=1/(2J)=62.5 ms,此时正弦为一;在 t=1/J=125 ms,余弦为负一、正弦为零,得到 −Ix。实际长延迟需另外考虑弛豫。

习题9.4-2 四步滤波

图 9.4-3 的滤波保留 p=−3 吗?

解答保留,因为 p−1=−4,四个离散相位的因子都为一。有限循环按模四选择;只有排除其他可达级数或增加选择条件,才能把目标通路与这类混叠分量分开。

习题9.4-3 梯度面积

同核两段通路分别 p=2 和 p=−1,第一段面积为 A,第二段应取多少面积才能满足理想重聚焦?

解答2A−A₂=0,故第二段为同方向的 2A。异核时须带入相应旋磁比;扩散、有限脉冲与梯度误差还需另评估。

9.4.8 术语

密度矩阵
描述自旋体系统计状态的矩阵。
偏差密度算符
从单位矩阵背景分离出的状态偏差。
传播算符
将初始状态演化到后续状态的算符。
乘积算符
以单核算符乘积构成的状态展开基。
反相磁化
与耦合伙伴纵向状态相关的横向磁化。
相干级数
纯相干分量对共同纵向旋转的相位响应整数。
相干转移通路
脉冲实验中相干类型随时间变化的序列。
相位循环
改变脉冲及接收相位以选择信号贡献的多次扫描组合。
梯度选择
利用空间相位与重聚焦选择通路的方法。

9.4.9 参考文献

  1. Sørensen O W, Eich G W, Levitt M H, Bodenhausen G, Ernst R R. Product operator formalism for the description of NMR pulse experiments. Prog Nucl Magn Reson Spectrosc, 16:163–192(原文印刷年1983). 出版标识;原文首页与理论。
  2. Keeler J. Understanding NMR spectroscopy, Chapter 6: Product Operators. 2004版本. 作者讲义。
  3. Keeler J. Understanding NMR spectroscopy, Chapter 9: Coherence Selection: Phase Cycling and Gradient Pulses. 2004课程版本. 作者讲义。
  4. Morris G A, Freeman R. Enhancement of nuclear magnetic resonance signals by polarization transfer. J Am Chem Soc, 1979, 101:760–762. INEPT 原始研究。
  5. Levitt M H. Spin Dynamics: Basics of Nuclear Magnetic Resonance, 2nd ed. Wiley, 2008. 出版预览。