EBSD 实操：花样形成物理
Kikuchi 明带从哪里来：二束动力学与 Bloch 波逐步推导
从两条耦合振幅方程走到消光振荡，再用原子位置的电子密度理解通道效应、明暗带和真实 EBSD 强度。
先记住这三点
- 单次与多次散射
- 二束耦合与消光距离
- 波场、热漫散射和互易
1. Bragg 定律确定位置，但没有给出整幅强度
一条Kikuchi带可先用几何构造理解：晶体内的散射源向许多方向供给电子，满足Bragg条件的方向构成一对圆锥；圆锥与屏幕相交形成带边，晶面经过点源的延伸面给出带中心线。这解释“带在哪里”和“大致多宽”，却不能解释为何带内常是亮的、边缘常是暗的，以及为什么不同材料有不同精细结构。
原因是电子与晶体相互作用较强。已经散射的波还能再次散射回直射方向，或者散射到另一反射方向。能量在不同波束间交换，不再是各原子一次散射强度的简单相加。下面用只有直射束0和一个衍射束g的理想模型建立直觉：忽略吸收、非弹性损失和其他强反射，固定样品厚度方向z。这些前提必须随公式一起保留。
2. 两束振幅为什么构成一对耦合方程
把波写成两项：ψ=C₀(z)e^(2πik₀·r)+C_g(z)e^(2πik_g·r)。这里沿用|k|=1/λ的不含2π波矢，故相位指数显式带2π。晶体势的周期性可把相差倒易矢量的平面波互相耦合。采用一组常见相位约定，振幅满足 dC₀/dz=i(π/ξ_g)C_g，dC_g/dz=i(π/ξ_g)C₀+2πisC_g。ξ_g称消光距离，单位长度；s是激发误差，单位长度⁻¹。因此方程两侧均为振幅/长度。整体相位和s符号随约定可变，但本节最终强度只含s²。
初始条件取C₀(0)=1、C_g(0)=0。若恰好满足Bragg条件s=0，对第一式再求一次导数并代入第二式，得到 C₀″=−(π/ξ_g)²C₀。与简谐振子方程相同，解为 C₀=cos(πz/ξ_g)，C_g=i sin(πz/ξ_g)。于是I₀=cos²(πz/ξ_g)，I_g=sin²(πz/ξ_g)，相加恒为1，正是无吸收的强度守恒。
3. 不在精确 Bragg 条件时，自己求本征值
将方程写为 C′=iπMC，M=[[0,1/ξ_g],[1/ξ_g,2s]]。求 det(M−μI)=0，得到 μ²−2sμ−1/ξ_g²=0，两根 μ±=s±√(s²+1/ξ_g²)。共同的s只给整体相位；两根间距决定振幅交换频率。用初始条件解两个本征模式的系数，得到 I_g(t)=sin²[πt√(s²+1/ξ_g²)]/[1+(sξ_g)²]。
这个式子同时给出两种效应：分母让离Bragg越远时最大可转移强度变小；正弦项让强度随厚度振荡。令w=sξ_g是无量纲偏离参数，可写成 I_g=sin²[(πt/ξ_g)√(1+w²)]/(1+w²)。读文献时要先查s还是w，不能把一个0.1 nm⁻¹直接填进无量纲w的位置。
4. 手算消光振荡，知道“厚一些更亮”何时失效
取ξ_g=100 nm、s=0。t=25 nm时I_g=sin²(π/4)=0.5；t=50 nm时I_g=1；t=100 nm时I_g=0。薄膜变厚可能使某反射先变亮再变暗，因此不能直接用一条反射的亮度比较两个位置的相含量。真实多束、吸收和厚度分布会改变上述完全回零行为，但不会恢复简单的“亮度只看原子数”规则。
再取s=0.01 nm⁻¹，w=1，衍射束最大强度不超过1/2；t=50 nm时I_g=½sin²(π√2/2)≈0.317。这与同厚度、s=0时的1相差很大，说明局部取向、厚度和结构因子可同时改变亮度。在πt√(s²+ξ_g⁻²)≪1的薄极限下可用sin x≈x，式子回到I_g≈(πt/ξ_g)²；强偏离Bragg时仅有t≪ξ_g还不足以保证这个近似。此时单次散射近似才逐渐合理。
5. 从多束 Fourier 展开到 Bloch 波
本段统一改用带2π的角波矢K和G，定义K=2πk₀、G=2πg，二者单位仍为长度⁻¹。把周期势写作V(r)=Σ_GV_Ge^(iG·r)，波函数写作ψ=Σ_GC_Ge^[i(K+G)·r]。代入稳态Schrödinger方程 [−ℏ²∇²/(2m)+V]ψ=Eψ，再比较每一个e^[i(K+G)·r]的系数，得到 [ℏ²|K+G|²/(2m)−E]C_G+Σ_HV_(G−H)C_H=0。
这一步的关键是乘法产生了G−H差项，所以每一束与其他束相耦合。截取有限反射集后，它成为矩阵本征问题，特征向量给出Bloch模式中各波束权重。高能电子计算还要采用相对论修正、吸收和边界匹配；这些不是把m任意改个数字就全部解决。课程给出方程如何形成，实际多束数值计算应交给经过验证的衍射程序，并检查反射截断、能量和深度积分的收敛。
6. 为什么原子之间的波场更容易通道传播
两束波相加后的概率密度包含交叉项：|ψ|²=|C₀|²+|C_g|²+2Re[C₀C_g*e^(−iG·r)]。交叉项使电子密度在晶格周期内有高低起伏。若高密度主要落在原子列附近，电子更容易发生热漫散射；若主要落在原子列之间，散射损失较少，波能更深入传播。这是通道效应的波动解释，并不是把电子想成沿一条空管直线滑行。
对某个出射方向u，EBSD强度可概念性写成 I(u)=∫w(r,E)S(r,E)|ψ_u(r,E)|²drdE。S表示散射源，w包含深度、能量和逸出权重。互易原理允许用反向入射到该u方向所形成的晶体波场计算出射概率。ECP记录改变入射方向时的背散射总量，EBSD记录许多出射方向；两者相关，但不是同一个测量操作。
7. 这套推导怎样帮助读真实花样
实际应用先看晶面几何是否匹配，再问强度差来自哪里。热漫散射提供非相干电子源，多束干涉决定精细带形，探测器点扩散与能量分布进一步模糊细节。因此弱带、强背景、变形造成的带展宽不能互相替代。动态背景应先用低通滤波估计缓慢变化的背景，再相减或相除；高通结果才是去背景后的细节图。
Winkelmann等的动力学模拟与教材铝花样比较说明，模拟精细强度可用于检验晶体结构和伪对称候选；书中同时指出实验的过剩/亏损不对称还受非均匀照明影响。实际教程应保留原花样、同一几何的模拟花样以及残差图：若残差集中在所有带的一侧，先排查照明和几何；若集中在结构敏感的特定反射，才进一步比较结构候选。不要把动力学模拟拟合出的一个最优解当成无误差的结构证明。
来源与延伸阅读
以下链接用于核对概念、标准状态或原始文献。本页为中文教学整理，不复制受版权保护的标准正文。