第1章 偶极作用与距离核

第一部分

物理与实验基础

建立偶极相互作用、距离核、距离分布、自旋中心及实验前置表征的共同基础。

第一部分 · 第1章

第1章　电子—电子偶极作用与距离核

从磁偶极相互作用到脉冲偶极谱的距离编码

1.1 章节范围与符号约定

电子—电子距离测量利用磁偶极耦合对自旋间距和取向的依赖。对于空间位置明确、交换耦合可忽略的两个电子自旋，偶极频率按距离的三次方反比变化。CW线形分析直接观察偶极作用引起的谱线分裂或展宽；脉冲偶极谱在时域记录由同一相互作用产生的相位演化。

本章统一采用无量纲自旋算符 \(\hat{\mathbf S}_A\) 和 \(\hat{\mathbf S}_B\)。哈密顿量 \(\hat H\) 使用能量单位；\(\hat{\mathcal H}=\hat H/\hbar\) 使用角频率单位。普通频率记为 \(\nu=\omega/2\pi\)。连接两自旋中心的矢量为 \(\mathbf r=\mathbf r_B-\mathbf r_A\)，距离为 \(r=|\mathbf r|\)，单位矢量为 \(\mathbf e_r=\mathbf r/r\)。静磁场 \(\mathbf B_0\) 定义实验室坐标系的 \(z\) 轴，\(\theta\) 是 \(\mathbf r\) 与 \(\mathbf B_0\) 的夹角。

符号限制

不同教材对交换常数采用 \(J\)、\(-J\) 或 \(-2J\) 作为海森堡哈密顿量的系数。报告引用交换常数时同时写出哈密顿量定义，避免只比较数值和正负号。

1.2 经典磁偶极相互作用

磁矩 \(\boldsymbol\mu_A\) 在位置 \(\mathbf r\) 产生的磁场为

$$\mathbf B_A(\mathbf r)=\frac{\mu_0}{4\pi r^3}\left[3\mathbf e_r(\boldsymbol\mu_A\cdot\mathbf e_r)-\boldsymbol\mu_A\right].$$ (1.1)

第二个磁矩在该局域磁场中的相互作用能为 \(E=-\boldsymbol\mu_B\cdot\mathbf B_A\)，由此得到

$$E_{\mathrm{dd}}=\frac{\mu_0}{4\pi r^3}\left[\boldsymbol\mu_A\cdot\boldsymbol\mu_B-3(\boldsymbol\mu_A\cdot\mathbf e_r)(\boldsymbol\mu_B\cdot\mathbf e_r)\right].$$ (1.2)

式(1.2)包含两个特征：相互作用强度按 \(r^{-3}\) 衰减；能量同时依赖两个磁矩相对于自旋间矢量的方向。偶极作用因此同时携带距离信息和取向信息。

B₀ ∥ z A B r = r eᵣ θ μA μB

图1.1　两个磁矩、静磁场与自旋间矢量的几何定义。图中磁矩仅作示意；量子力学处理中以自旋算符和相应的g张量代替经典矢量。

1.3 量子哈密顿量与张量形式

1.3.1 各向同性g值的点偶极形式

电子磁矩算符可写为 \(\hat{\boldsymbol\mu}=-g\mu_B\hat{\mathbf S}\)。代入式(1.2)，两个点电子自旋之间的偶极哈密顿量为

$$\hat H_{\mathrm{dd}}=\frac{\mu_0}{4\pi}\frac{g_Ag_B\mu_B^2}{r^3}\left[\hat{\mathbf S}_A\cdot\hat{\mathbf S}_B-3(\hat{\mathbf S}_A\cdot\mathbf e_r)(\hat{\mathbf S}_B\cdot\mathbf e_r)\right].$$ (1.3)

定义能量尺度

$$d_0(r)=\frac{\mu_0}{4\pi}\frac{g_Ag_B\mu_B^2}{r^3},\qquad \omega_0(r)=\frac{d_0(r)}{\hbar},\qquad \nu_0(r)=\frac{d_0(r)}{h}.$$ (1.4)

式(1.4)分别对应能量、角频率和普通频率表示。后文的PDS频率以 \(\nu_0\) 表示。

1.3.2 各向异性g张量

当两个自旋具有各向异性g张量时，磁矩方向一般不与自旋角动量方向完全平行。偶极作用可写为双线性张量形式

$$\hat H_{\mathrm{dd}}=\hat{\mathbf S}_A^{\mathrm T}\mathbf D_{AB}\hat{\mathbf S}_B,$$ (1.5)
$$\mathbf D_{AB}=\frac{\mu_0\mu_B^2}{4\pi r^3}\mathbf g_A^{\mathrm T}\left(\mathbf 1-3\mathbf e_r\mathbf e_r^{\mathrm T}\right)\mathbf g_B.$$ (1.6)

式(1.6)显示，测得的偶极频率同时受分子取向、g张量主轴方向及自旋间矢量影响。对于g各向异性较小的有机自由基，常用有效g值近似；Cu(II)、低对称过渡金属中心和部分高自旋体系通常需要保留张量结构。

适用前提

点偶极模型把每个未配对电子的磁矩压缩到一个位置。该模型适合自旋密度空间尺度明显小于中心间距的情况。距离进入亚纳米范围或自旋密度在共轭骨架、金属簇和配体间显著离域时，应采用分布磁矩模型或电子结构计算得到的耦合张量。

1.4 高场近似、世俗项与赝世俗项

1.4.1 高场弱耦合条件

静磁场下的主要项是电子Zeeman哈密顿量

$$\hat{\mathcal H}_Z=\omega_A\hat S_{Az}+\omega_B\hat S_{Bz}.$$ (1.7)

当Zeeman能量显著大于偶极耦合时，与总磁量子数变化相关的快速振荡项在一阶处理中被平均。对于两个 \(S=1/2\) 自旋，保留到高场世俗近似后的偶极项可写为

$$\hat{\mathcal H}_{\mathrm{dd}}^{\mathrm{sec}}=d(\theta)\hat S_{Az}\hat S_{Bz}-\frac{d(\theta)}{4}\left(\hat S_A^+\hat S_B^-+\hat S_A^-\hat S_B^+\right),$$ (1.8)
$$d(\theta)=\omega_0(r)\left(1-3\cos^2\theta\right).$$ (1.9)

式(1.8)第一项是对角的 \(\hat S_{Az}\hat S_{Bz}\) 项。第二项交换两个自旋的磁量子数，常称赝世俗项或flip-flop项。Goldfarb与Stoll教材第20章采用 \(A\) 和 \(B\) 两类项描述这一结构。

1.4.2 Zeeman频率差与flip-flop项

定义 \(\Delta\omega_{AB}=|\omega_A-\omega_B|\)。当

$$\Delta\omega_{AB}\gg |d(\theta)|,$$ (1.10)

flip-flop过程偏离能量匹配条件，其一阶作用可忽略，得到PDS最常用的弱耦合哈密顿量

$$\hat{\mathcal H}_{\mathrm{dd}}^{\mathrm{weak}}=d(\theta)\hat S_{Az}\hat S_{Bz}.$$ (1.11)

当 \(\Delta\omega_{AB}\) 与 \(|d|\) 可比，式(1.11)不再充分。能级混合、赝世俗作用和跃迁概率变化会改变谱线分裂及偶极演化。短距离、窄线同类自由基、高自旋中心和频率接近的选定跃迁更容易进入该区域。

近似层级主要条件保留内容典型用途
完整张量哈密顿量一般取向、各向异性g、有限耦合完整 \(\mathbf D_{AB}\) 及Zeeman项数值模拟、金属中心、取向选择
高场世俗近似Zeeman作用远大于偶极作用\(S_{Az}S_{Bz}\)与flip-flop项单频率PDS、强耦合边界分析
弱耦合异核式近似\(|\Delta\omega_{AB}|\gg|d|\)仅 \(S_{Az}S_{Bz}\)常规DEER距离核
点偶极与各向同性g自旋局域、g各向异性较小\(r^{-3}(1-3\cos^2\theta)\)常规有机自由基距离分布
术语边界

“世俗”取决于所选零阶哈密顿量和旋转坐标系。某一文献中的赝世俗项在另一种强场或跃迁选择近似下可能被进一步删除。报告后续各方法章节会写明采用的零阶哈密顿量和频率分离条件。

1.5 角度依赖、魔角与Pake谱形

在式(1.11)的条件下，普通偶极频率为

$$\nu_{\mathrm{dd}}(r,\theta)=\nu_0(r)\left(1-3\cos^2\theta\right).$$ (1.12)
取向角度频率物理特征
自旋间矢量平行于磁场\(\theta=0^\circ\)\(-2\nu_0\)角度函数的一端极值
魔角\(\theta=\arccos(1/\sqrt3)=54.7356^\circ\)0一阶点偶极频率为零
自旋间矢量垂直于磁场\(\theta=90^\circ\)\(+\nu_0\)角度函数的另一端极值

各向同性冻结样品包含全部空间取向。对式(1.12)作粉末平均可得到Pake谱形，其特征位置对应平行和垂直取向。谱形不对称来自角度函数与球面取向权重共同作用。CW实验可直接观察由该分布形成的分裂或展宽；PDS实验先在时域记录振荡，再通过傅里叶变换得到相应的偶极频谱。

取向选择

标准Pake核要求被观测的分子取向近似均匀。g和A各向异性、有限脉冲带宽及特定磁场位置可能只激发部分取向。此时测得频率权重同时包含空间取向分布和激发选择函数，需要在具体方法章节中单独处理。

1.6 距离标度与时域距离核

1.6.1 偶极频率的数量级

对于 \(g_A\approx g_B\approx2.0023\)，点偶极常数可写成便于实验估算的形式

$$\nu_0(r)\approx\frac{52.04\ \mathrm{MHz\,nm^3}}{r^3}.$$ (1.13)
距离 \(r\) / nm\(\nu_0\) / MHz特征周期 \(1/\nu_0\)实验含义
152.0419.2 ns有限脉冲、交换及点偶极误差需要重点评估
26.5050.154 μsCW和短时间PDS可能同时适用
31.9270.519 μs常规PDS容易获得多个演化特征
40.8131.23 μs需要微秒量级相位记忆时间
50.4162.40 μs背景与时间窗开始明显影响分布宽度
80.1029.84 μs依赖长 \(T_m\)、氘代或多脉冲方案
100.052019.2 μs对标签相干性和信号稳定性要求很高

表中周期只反映 \(\theta=90^\circ\) 的 \(\nu_0\)，不等同于实验所需的最小演化时间。粉末平均、距离分布、背景衰减、噪声和相干衰减会共同决定可辨识距离。

1.6.2 单距离粉末平均核

对于各向同性取向、孤立自旋对和理想脉冲，单一距离 \(r\) 的归一化偶极核可写为

$$K(r,t)=\int_0^{\pi/2}\sin\theta\,\cos\left[2\pi\nu_0(r)(1-3\cos^2\theta)t\right]d\theta.$$ (1.14)

若样品具有归一化距离分布 \(P(r)\)，分子内偶极演化函数为

$$D(t)=\int_0^\infty P(r)K(r,t)dr,\qquad \int_0^\infty P(r)dr=1.$$ (1.15)

式(1.15)是多数双自旋PDS距离反演的基础。实际观测信号还包含调制深度、分子间背景、相干衰减、残留偶极路径和仪器响应。相应的前向模型在数据分析章节中展开。

距离分辨限制

有限时间窗相当于对 \(D(t)\) 截断，低频信息最先受影响。长距离对应较低的偶极频率，因此对可用演化时间和背景模型更敏感。距离分布的窄峰还会受到正则化、参数化模型及标签构象的影响。

1.7 交换耦合及其对距离测量的影响

交换作用来自电子不可分辨性、轨道重叠和电子间静电相互作用。常用各向同性交换哈密顿量之一为

$$\hat H_{\mathrm{ex}}=J\hat{\mathbf S}_A\cdot\hat{\mathbf S}_B.$$ (1.16)

式(1.16)的符号只适用于本章定义。其他文献使用 \(-J\hat{\mathbf S}_A\cdot\hat{\mathbf S}_B\) 或 \(-2J\hat{\mathbf S}_A\cdot\hat{\mathbf S}_B\)。交换作用通常随轨道重叠快速衰减，其大小和符号对桥连结构、共轭路径、构象与溶剂十分敏感，缺少可跨体系使用的单一距离标尺。

在高场世俗形式中，交换项可同时影响对角项和flip-flop项。Goldfarb与Stoll第20章在其约定下写为 \(a=d+2J\) 和 \(b_{\mathrm{ff}}=-d/2+2J\)。这说明实验频率可能同时包含具有角度依赖的偶极分量和近似各向同性的交换分量。

实验区域交换作用的表现对距离结论的影响建议处理
弱交换附加相位或小幅谱线位移可能造成距离峰偏移或展宽采用包含 \(J\) 的前向模型并与不同取向或序列比较
与偶极作用相当谱线与时域频率同时改变单一Pake核不再充分联合模拟 \(D\)、\(J\)、Zeeman差和有限脉冲
强交换形成总自旋多重态，跃迁结构重组独立 \(S=1/2\) 自旋对模型失效采用耦合基组和总自旋哈密顿量

1.8 自旋密度离域与有效偶极张量

点偶极公式默认两个磁矩位于明确的几何点。真实自由基和金属配合物的未配对自旋密度分布在原子轨道及配体上。两个归一化自旋密度 \(\rho_A(\mathbf r_A)\) 和 \(\rho_B(\mathbf r_B)\) 之间的有效偶极张量可概念性地写为

$$\mathbf D_{AB}^{\mathrm{eff}}=\frac{\mu_0\mu_B^2}{4\pi}\iint \rho_A(\mathbf r_A)\rho_B(\mathbf r_B)\frac{\mathbf g_A^{\mathrm T}[\mathbf 1-3\mathbf e\mathbf e^{\mathrm T}]\mathbf g_B}{|\mathbf r_B-\mathbf r_A|^3}\,d\mathbf r_A d\mathbf r_B,$$ (1.17)

其中 \(\mathbf e=(\mathbf r_B-\mathbf r_A)/|\mathbf r_B-\mathbf r_A|\)。式(1.17)说明实验得到的是两片自旋密度之间的磁相互作用。把该张量换算为单一 \(r_{\mathrm{eff}}\) 需要指定点中心和取向模型。

距离概念定义常见来源
电子自旋有效距离由测得偶极张量在指定模型下换算的距离PDS或CW线形
标签中心距离两个标签代表点之间的几何距离rotamer library、MD、结构模型
连接位点距离两个标记残基或配位位点之间的距离蛋白或核酸结构
骨架原子距离指定主链原子之间的距离晶体、cryo-EM、NMR或计算模型

这些距离具有不同定义。柔性连接臂、rotamer占比、局部静电环境和构象系综都会改变电子自旋距离分布。后续结构解释章节将把标签模型与分子构象模型分开处理。

1.9 高自旋、多自旋和分子间作用

1.9.1 高自旋中心

当一侧或两侧自旋满足 \(S>1/2\)，多个 \(m_S\leftrightarrow m_S\pm1\) 跃迁可能落入激发带宽。各跃迁具有不同的矩阵元、极化、弛豫速率和有效偶极频率。Gd(III)–Gd(III) RIDME中的基频及高次谐波属于多能级跃迁效应；Cu(II)和其他金属中心还会受到g、A各向异性与取向选择影响。

高自旋体系的距离核需要明确观测跃迁、伙伴自旋跃迁概率和混合时间条件。直接套用两个 \(S=1/2\) 自旋的式(1.15)可能产生系统性距离偏差。

1.9.2 多自旋体系

含三个及以上自旋的分子中，一个observer自旋可能同时与多个伙伴耦合。单泵浦、弱相关条件下，时域信号常可近似分解为多个自旋对因子的乘积。多次泵浦、较高极化或相关相干路径会引入额外偶极路径，简单相加多个距离分布不足以描述全部信号。

1.9.3 分子内信号与分子间背景

同一分子或复合物内、位于目标距离范围内的自旋构成分子内贡献。其他分子上的远距离自旋形成浓度相关的背景衰减。均匀三维分布常用指数背景描述；膜、表面、排除体积和非均匀分布可能需要伸展指数或显式空间模型。

$$V(t)=F_{\mathrm{intra}}(t)B_{\mathrm{inter}}(t)$$ (1.18)

式(1.18)是稀释自旋簇条件下常用的因子化近似。背景模型偏差会传递到低频偶极成分，对长距离和宽分布影响尤为明显。

1.10 偶极作用在各测距方法中的观测形式

方法偶极作用的观测形式本章公式的主要使用方式需另行处理的因素
CW偶极展宽频域分裂、展宽和线形变化完整或高场偶极哈密顿量交换作用、运动平均、本征线宽
DEER／PELDORpump翻转伙伴自旋后产生的observer相位调制弱耦合 \(S_{Az}S_{Bz}\) 核双频激发、取向选择、残留偶极路径
DQC双量子相干路径中的偶极相位演化高场偶极哈密顿量和密度算符传播强脉冲、相位瞬变、死时间
SIFTER／SIDRE单频率共同激发下的偶极重聚世俗项与flip-flop项的序列传播谱带覆盖、整形脉冲动态相位、过滤效应
RIDME伙伴自旋在混合期发生纵向弛豫跃迁后产生调制偶极频率与跃迁概率的组合核\(T_1\)、混合时间、高自旋谐波
光诱导PDS光生三重态或自由基出现前后的偶极相位差偶极哈密顿量加三重态ZFS非玻尔兹曼布居、三重态寿命和光激发效率

具体序列章节将从本章的相互作用哈密顿量出发，依次建立相干路径、相位累积区间、相位循环、实验参数、伪影和数据模型。相位表直接保存在相应方法章节中。

1.11 本章结论与适用边界

结论条件超出条件后的处理
偶极耦合按 \(r^{-3}\) 变化两个局域点磁矩使用自旋密度积分或电子结构计算
角度因子为 \(1-3\cos^2\theta\)各向同性g及高场量子化保留g张量和完整偶极张量
常规PDS可使用 \(S_{Az}S_{Bz}\) 项Zeeman频率差远大于偶极耦合保留flip-flop项并数值传播
Pake核描述均匀粉末取向全取向激发或选择效应可忽略加入取向选择函数和脉冲激发轮廓
单一双自旋核可反演距离分布孤立自旋对、背景可分离采用多路径、多自旋或谐波模型
交换作用可忽略轨道重叠很小且实验未显示交换特征联合拟合偶极与交换哈密顿量

本章建立的是后续实验章节共同使用的最低理论模型。具体实验中的可观测频率还会受到激发选择、相干路径、有限脉冲、弛豫和背景的共同影响。

返回第一部分目录