第3章 相干路径与相位循环

第一部分 · 第3章

第3章 脉冲、相干路径与相位循环基础

从旋转坐标系传播到离散相位投影

3.1 章节范围与符号约定

脉冲偶极谱通过一组微波脉冲建立、转换、储存并检测电子自旋相干。相位循环利用不同相干路径对脉冲相位的不同响应,将目标回波与FID、非目标回波、接收偏置和部分脉冲缺陷区分开。本章给出所有PDS序列共用的数学语言和验证规则。

旋转坐标系的横向轴记为 \(x,y\),静磁场方向记为 \(z\)。脉冲相位为 \(\phi\),翻转角为 \(\beta\),相干阶数为 \(p\)。密度算符记为 \(\hat\rho\),脉冲前后的相干阶数差记为 \(\Delta p=p_{\mathrm{after}}-p_{\mathrm{before}}\)。接收机复相信号为 \(S=I+iQ\)。

相位约定

主动旋转与被动坐标旋转、\(S^+=S_x+iS_y\) 的定义、傅里叶变换符号和接收机I/Q定义都会改变表中相位的正负号。任何完整相位表都应同时给出脉冲相位、接收相位、信号组合式和采用的约定。只抄录软件数组无法完成跨平台核对。

3.2 密度算符、自由演化与脉冲传播

3.2.1 自由演化

忽略弛豫时,密度算符满足Liouville–von Neumann方程

$$\frac{d\hat\rho}{dt}=-\frac{i}{\hbar}[\hat H,\hat\rho],\qquad \hat\rho(t)=\hat U(t)\hat\rho(0)\hat U^\dagger(t),\qquad \hat U(t)=\exp(-i\hat Ht/\hbar).$$ (3.1)

弛豫可通过主方程、Liouville空间超算符或经验衰减包络加入。相位循环本身只选择相位响应相符的信号;它不会恢复已经因 \(T_m\)、瞬时扩散或谱扩散丢失的相干。

3.2.2 矩形脉冲

对共振、理想硬脉冲,绕横向相位轴 \(\phi\) 旋转 \(\beta\) 的传播子为

$$\hat R(\beta,\phi)=\exp\{-i\beta[\hat S_x\cos\phi+\hat S_y\sin\phi]\}.$$ (3.2)

共振自旋的翻转角为 \(\beta=\omega_1t_p\)。对一般跃迁 \(m_S\leftrightarrow m'_S\),Rabi频率还包含跃迁矩阵元:

$$\omega_{\mathrm{nut}}=\omega_1\sqrt{S(S+1)-m_Sm'_S}.$$ (3.3)

离共振时,自旋绕由 \(\Delta\omega\) 和 \(\omega_1\) 合成的有效场旋转,\(\omega_{\mathrm{eff}}=(\Delta\omega^2+\omega_1^2)^{1/2}\)。有限脉冲因此同时改变翻转角、旋转轴和相位。

模型采用条件可解释内容遗漏因素
理想瞬时脉冲\(\omega_1\)远大于频率偏移和内相互作用相干阶数、回波时间、最低相位循环带宽、脉冲内演化、动态相位
有限矩形脉冲幅度和相位可视为分段常数激发轮廓、离共振响应、脉冲重叠谐振腔瞬态、放大器非线性
整形脉冲已知幅相波形和硬件响应宽带反转、扫频方向、动态相位未经测量的传输函数和B1不均匀性

3.3 相干阶数及其相位响应

相对于总电子自旋投影算符 \(\hat S_z\),相干阶数为 \(p\) 的算符分量 \(\hat O_p\) 满足

$$[\hat S_z,\hat O_p]=p\hat O_p.$$ (3.4)

绕 \(z\) 轴相移 \(\phi\) 后,该分量获得确定相位因子:

$$e^{-i\phi\hat S_z}\hat O_pe^{+i\phi\hat S_z}=e^{-ip\phi}\hat O_p.$$ (3.5)
相干阶数典型密度算符分量物理含义直接检测
\(p=0\)\(S_z\)、布居、纵向双自旋序布居或纵向储存不能由正交线圈直接检测
\(p=+1,-1\)\(S^+\)、\(S^-\)单量子横向相干形成复数FID或回波
\(p=+2,-2\)\(S_A^+S_B^+\)、\(S_A^-S_B^-\)双量子相干需转换回 \(p=\pm1\)
更高阶多自旋升降算符乘积多量子相干需经过路径选择和再转换

对一个脉冲,目标路径在该脉冲处发生 \(\Delta p\) 的相干阶数变化。改变脉冲相位 \(\Delta\phi\) 后,目标信号相位变化满足

$$\Delta\phi_{\mathrm{sig}}=-\Delta p\,\Delta\phi.$$ (3.6)

式(3.6)是教材第11.6.4节给出的相位循环基本规则。若采用共轭的复信号定义,等式整体符号相反,但不同路径之间的相对判别不变。

3.4 FID、初级回波与受激回波

3.4.1 单脉冲FID

平衡布居经 \(\pi/2\) 脉冲转换为可检测的单量子相干:

$$0\longrightarrow \pm1.$$ (3.7)

谱包的频率偏移不同,横向相位随时间散开。宽谱固体的FID可能落在谐振腔ringdown和接收死时间内,因此许多脉冲EPR实验采用回波检测。

3.4.2 Hahn回波

两脉冲序列 \((\pi/2)-\tau-(\pi)-\tau-\mathrm{echo}\) 的目标路径为

$$0\longrightarrow \pm1\longrightarrow \mp1.$$ (3.8)

\(\pi\) 脉冲把频率偏移导致的相位演化方向反转,谱包在第二个 \(\tau\) 后重聚。该重聚消除静态频率不均匀性;可逆或不可逆的相互作用能否重聚取决于相互作用哈密顿量在脉冲作用下的变换。

3.4.3 受激回波

三脉冲序列可将横向相干的一部分转换为 \(p=0\) 分量,在混合期储存,再由第三个脉冲读出。目标路径可概括为

$$0\longrightarrow \pm1\longrightarrow 0\longrightarrow \mp1.$$ (3.9)

同一组三个脉冲还会产生三个两脉冲回波、一个重聚的两脉冲回波和多个FID。延迟改变时,这些响应可能与目标受激回波相交,形成可重复的结构性伪影。

模块化含义

Hahn回波提供相干建立与检测模块;受激回波提供纵向储存模块。DEER、RIDME、DQC和SIFTER的长序列由这些模块及相互作用反转或相干转换步骤组合而成。序列分析应逐段标明相干阶数、有效哈密顿量和相位累积符号。

3.5 耦合自旋中的同相、反相与多量子项

对弱耦合的两个 \(S=1/2\) 自旋,采用角频率哈密顿量 \(\hat{\mathcal H}_{dd}=d\hat S_{Az}\hat S_{Bz}\)。当B自旋处于 \(m_B=\pm1/2\) 时,A自旋的条件频率偏移为 \(\pm d/2\)。A自旋单量子相干在演化中分解为同相与反相分量:

$$\hat S_{A+}(t)=\hat S_{A+}\cos\left(\frac{dt}{2}\right)-2i\hat S_{A+}\hat S_{Bz}\sin\left(\frac{dt}{2}\right).$$ (3.10)

式(3.10)的第二项是对B自旋状态呈相反相位的反相相干。伙伴自旋翻转、纵向弛豫跃迁或共同激发会改变两项的相对演化,由此形成DEER、RIDME及单频率PDS中的偶极调制。

算符类别示例相干阶数序列中的作用
同相单量子\(S_{Ax},S_{Ay}\)\(\pm1\)直接形成回波信号
反相单量子\(2S_{Ax}S_{Bz}\)\(\pm1\)记录伙伴自旋状态依赖的相位
零量子\(S_A^+S_B^-\)0flip-flop及相干交换
双量子\(S_A^+S_B^+\)+2DQC建立和过滤
系数约定

式(3.10)采用无量纲自旋算符及 \(\hat{\mathcal H}_{dd}=dS_{Az}S_{Bz}\)。若文献把哈密顿量写成 \(2dS_{Az}S_{Bz}\)、使用普通频率或把两条条件跃迁的间隔定义为偶极频率,三角函数自变量会出现2或 \(2\pi\) 的差异。比较序列推导时必须先统一哈密顿量单位。

3.6 相干转移路径与回波条件

含 \(n\) 个脉冲的路径可写为相干阶数序列 \(\{p_0,p_1,\ldots,p_n\}\),其中 \(p_0=0\)。第 \(k\) 个脉冲处的相干转移阶数为 \(\Delta p_k=p_k-p_{k-1}\)。若各自由演化区间长度为 \(t_j\),静态频率偏移产生的总相位为

$$\Phi_\Omega=\Omega\sum_j p_jt_j.$$ (3.11)

理想回波中心满足 \(\sum_jp_jt_j=0\)。偶极、超精细、交换和零场分裂项在各区间的符号未必与频率偏移相同,因此仍可在回波中心保留目标相位。

若第 \(k\) 个脉冲相位为 \(\phi_k\),该路径的脉冲相位标签为

$$\Phi_{\mathrm{path}}=-\sum_{k=1}^{n}\Delta p_k\phi_k.$$ (3.12)

接收相位 \(\phi_R\) 选取与目标路径同相的复数分量。实验相位表的每一行实质上都是式(3.12)的一次采样。

核查项目目标路径非目标路径
回波时间在采集窗内满足重聚条件可能随扫描变量穿过采集窗
脉冲相位标签与接收权重匹配经完整循环求和后抵消
偶极演化保留预定符号和时间长度可能形成额外偶极路径
弛豫权重在序列时标内可检测即使未完全抵消也可能很弱

3.7 相位循环的离散投影原理

3.7.1 单脉冲相位维度

令某一脉冲依次取 \(N\) 个等间隔相位 \(\phi_j=2\pi j/N\)。对目标转移阶数 \(\Delta p_t\),可用权重

$$w_j=\exp(+i\Delta p_t\phi_j)$$ (3.13)

组合复数瞬态:

$$S_{\mathrm{sel}}=\frac1N\sum_{j=0}^{N-1}w_jS(\phi_j).$$ (3.14)

具有转移阶数 \(\Delta p_u\) 的非目标路径保留系数为

$$\frac1N\sum_{j=0}^{N-1}e^{-i(\Delta p_u-\Delta p_t)2\pi j/N}=\begin{cases}1,&\Delta p_u\equiv\Delta p_t\pmod N\\0,&\text{其他情况。}\end{cases}$$ (3.15)

因此,\(N\) 步循环只能区分模 \(N\) 不同的相干转移。较高阶路径与目标路径同余时,需要增加相位步数、循环其他脉冲或利用时间位置、频率选择与梯度等条件分离。

3.7.2 多脉冲循环

多个相位维度的通用组合式为

$$S_{\mathrm{sel}}=\sum_{j=1}^{N_{\mathrm{cyc}}}w_jS(\phi_{1j},\phi_{2j},\ldots,\phi_{nj}),\qquad w_j=e^{i\phi_{R,j}}.$$ (3.16)

完整笛卡尔积循环的步数等于各相位维度步数的乘积。实际序列可利用对称性、预期很弱的路径或时间分离缩短循环;缩短后的版本属于实验条件相关的最低循环,应通过模拟和对照采集核查。

3.8 通用相位循环示例

3.8.1 二步相位交替

对一个使目标路径发生奇数阶转移的脉冲,取 \(+x\) 和 \(-x\) 两个相位并作差,可保留随该脉冲反号的信号,消除与该脉冲无关的静态偏置:

步脉冲相位接收权重组合
1\(+x\)+1\(S=S_{+x}-S_{-x}\)
2\(-x\)−1

该循环不能普遍区分所有回波,只适合删除对所循环脉冲不敏感的成分。其作用对象应由路径表确认。

3.8.2 三脉冲受激回波四步示例

Goldfarb与Stoll教材图11.13给出循环第二、第三脉冲的四步方案。采用教材的 \(x/-x\) 标记,组合关系为:

步\(\phi_1\)\(\phi_2\)\(\phi_3\)权重
1\(+x\)\(+x\)\(+x\)−1
2\(+x\)\(+x\)\(-x\)+1
3\(+x\)\(-x\)\(+x\)+1
4\(+x\)\(-x\)\(-x\)−1
$$S_{\mathrm{stim}}=-S_1+S_2+S_3-S_4.$$ (3.17)

该组合删除图示条件下的三个初级回波和重聚初级回波,保留受激回波。接收偏置、I/Q不平衡和仪器相位误差还可通过附加循环处理。该示例用于说明投影方法,不作为任何PDS序列的通用相位表。

完整数组的归属

DEER、DQC、SIFTER和RIDME的相位表包含序列特有的目标路径、多路径、扫描时间定义和伪回波。报告把这些数组放在相应方法章节中,并在表旁列出目标 \(\{p_j\}\)、接收相位及消除对象。光诱导方法暂不整理具体数组。

3.9 嵌套循环、时间平均与采集顺序

相位循环通常与核调制平均、正负时间采集、延迟步进、频率或磁场步进共同嵌套。总瞬态数可写为

$$N_{\mathrm{tr}}=N_tN_{\mathrm{pc}}N_{\tau}N_{\mathrm{field}}N_{\mathrm{avg}}N_{\mathrm{shot}}.$$ (3.18)
循环维度主要目的与相位循环的区别
脉冲相位选择相干路径、删除偏置和非目标回波按复相位权重线性组合
\(\tau\) 平均平均核调制或移动交叉回波改变物理演化时间,不等同于相位投影
正负时间轴校准零点、检验对称性或获取完整路径改变扫描脉冲的时间位置
频率/磁场循环测量选择性和取向依赖改变被激发的自旋包
扫描顺序随机化减轻慢漂移与时间轴的相关性改变采集次序,不改变目标路径

相位循环应尽可能作为最内层快速循环,使增益、相位和基线漂移在一个循环内近似恒定。若谱仪重编程时间很长,可采用分块采集,但每一块需要包含完整循环并记录时间戳。

3.10 相位相关伪影与诊断

现象可能来源诊断采集处理原则
固定时间尖峰或台阶非目标回波穿过采集窗、脉冲重叠改变固定延迟或扫描方向调整时间几何并使用针对性循环
奇偶点交替相位表顺序、触发或接收解调不同步保存每个相位步的未组合瞬态逐步核对相位与接收权重
虚部残留明显零阶相位错误、I/Q不平衡、脉冲相位瞬变比较复数回波及参考样品先校准仪器相位,再评估模型
改变泵浦功率后出现额外调制相干泵浦、部分激发或额外偶极路径功率、长度和频率阶梯实验使用有限脉冲/多路径模型
缩短循环后出现窄峰未消除路径与目标路径同余完整循环与缩短循环对照恢复必要的相位维度
整形脉冲扫频方向相关动态相位、脉冲内演化、硬件传输函数反转扫频方向和相位校正波形并进行全传播模拟

相位循环能消除相位标签不同的相干响应。与目标路径具有相同相位标签的背景、随机噪声、不可逆弛豫及某些脉冲内演化不会因增加循环步数而消失。

3.11 相位表的记录、验证与发布规则

记录字段最低内容核查方法
序列定义全部脉冲、频率通道、固定与扫描延迟与原论文图和时间定义逐项对应
相位约定\(+x,+y,-x,-y\) 与数值角度映射单脉冲或Hahn回波参考测量
目标路径每个区间的 \(p_j\) 和每个脉冲的 \(\Delta p_j\)式(3.11)和式(3.12)
相位表每一步的全部被循环脉冲相位离散投影或密度算符模拟
接收组合接收相位或复权重、加减公式从单步瞬态重算组合信号
消除对象FID、初级回波、交叉回波、偏置等删减循环对照和延迟位移实验
适用边界理想/有限脉冲、相位瞬变、整形脉冲条件功率、脉宽、偏置和扫频方向测试
报告采用的状态标记

已核实 表示相位表与原文、补充材料或教材原图一致,并完成路径检查;推导 表示由明确的路径和约定计算得到;待程序核对 表示来源只给出部分相位或厂商数组,尚缺接收组合或约定。

方法章节中的相位表将采用上述字段。若原论文只报告“采用若干步相位循环”而未给出数组,报告保留缺口并说明需要的补充材料,不凭经验补写成可直接执行的程序。

3.12 本章结论

结论实验含义
相干阶数描述密度算符对z轴相移的响应可用统一语言追踪不同PDS序列的相干转换
目标路径的脉冲相位由 \(-\sum\Delta p_k\phi_k\) 决定接收权重应与目标路径相位共轭匹配
等间隔相位循环是离散傅里叶投影最低步数取决于需区分的转移阶数及其模同余关系
延迟平均、正负时间采集和相位循环作用不同采集程序和报告中应分别记录
相位循环不修复弛豫、随机噪声和同标签背景需结合硬件校准、时间几何和前向模型诊断

返回第一部分目录