EBSD 实操:错取向统计

错取向直方图该和谁比较:随机、非关联与邻接分布

从 Haar 概率推到随机转角,再构建保留织构的配对基线,分清角度峰来自总体织构、空间邻接还是孪晶。

阅读设置

先记住这三点

  1. 随机取向的错取向角分布不平坦;对称约化会再次改变分布。
  2. 非关联配对保留所选权重下的织构,去掉指定的邻接信息。
  3. 真正讨论晶界偏好,需要同时看角、轴、邻接以及统计分母。
01

先固定旋转定义

设 RA、RB 都把各自晶体正交坐标变换到试样坐标。把晶粒 A 中同一物理向量改写到晶粒 B 坐标,用 MAB=RBᵀRA。相同相的对称等价形式为 SBᵀMABSA,其中 SA、SB 是各自允许的晶体正旋转对称操作。对所有等价形式取最小转角,才得到常用的 disorientation。

任何旋转 M 的特征值包含 1 和一对 exp(±iθ),所以 trM=1+2cosθ,得到 θ=arccos[(trM−1)/2]。这个角不包含转轴信息。两个边界即使同为 60°,轴不同,晶体学关系也可能完全不同。

02

从四元数体积元推导随机转角

先暂时不考虑晶体对称。旋转四元数写 q=(cos(θ/2),a sin(θ/2)),其中 a 为单位轴,0≤θ≤π,并只取代表 q 与 −q 中一个的半个三维单位球面。令 r=θ/2,球面体积元随 sin²r dr dΩa 变化;轴方向积分给 4π,半球总体积为 π²。

于是 pθ(θ)dθ=(4π/π²)sin²(θ/2)(dθ/2),即 pθ=(2/π)sin²(θ/2)=(1−cosθ)/π。积分得 F(θ)=(θ−sinθ)/π,角度必须用弧度。即使完全随机,接近零的小转角也比大转角少,这与均匀画一条水平直方图不同。

03

为什么立方晶体需要 Mackenzie 型约化

上式对每个真实旋转只保留轴角的 θ≤π,但尚未利用晶体对称。立方晶体有 24 个正旋转对称操作;每个随机旋转须映到离单位旋转最近的对称等价区域。这个区域的边界切掉大转角方向,因此不能沿用未约化的 pθ。

严格定义可写 Fcub(t)=∫SO(3)1[minS∈Gc angle(MS)≤t]dμ(M),dμ 是归一化 Haar 测度。小角区域中,围绕 24 个群元素的半径 t 球尚未重叠,因此当 0≤t≤π/4 时 Fcub(t)=24(t−sin t)/π。超过 π/4 球开始重叠,需要减去重叠体积,形成分段的 Mackenzie 分布;不能把前一段公式外推到全部角度。

04

手算随机低角边界的几何背景

取 t=10°=0.174533 rad,则 t−sin t≈0.0008847。按立方随机取向基线,Fcub≈24×0.0008847/π≈0.00676,即约 0.676%。这是随机独立取向、同为立方且按规定对称约化时的小角概率。

这个计算不意味着真实材料的低角晶界必定只有 0.676%。变形、再结晶、孪晶和空间聚集都会改变邻接关系;强织构即使完全随机安排晶粒位置,也会提高某些取向对出现的概率。此时需要下一层基线。

05

非关联分布是取向分布的配对卷积

令 f(R) 是相同权重定义下的归一化 ODF。独立取两个取向 R1、R2,令 M=R2ᵀR1,即 R1=R2M。利用 Haar 测度对乘法不变,M 的密度为 h(M)=∫f(R)f(RM)dμ(R),再按相应晶体对称约化。

若 f=1,h也为1,约化后的角分布回到随机基线;若 f 有峰,h 会反映峰内配对及不同峰之间的关系。这正是“非关联”与“随机织构”不同的原因:非关联指取向配对方式独立,并未要求 ODF 均匀。

06

有限列表怎样配对才可解释

有 N 个晶粒取向及权重 wi 时,选择所有 i≠j 的有向对。目标角区间 B 的非关联分数为 Σi≠jwiwj1[θij∈B]/[(Σwi)²−Σwi²]。排除自配对是为了避免每个晶粒与自己产生人为零角峰;若允许重复独立抽样,则必须说明这种自配对处理。

所有成对结果不是 N(N−1) 个独立实验,因为许多取向对共享同一个晶粒。不能把成对数量直接当二项式样本量计算极窄误差条。按独立视场或空间块讨论不确定性,通常更符合测量来源。

07

邻接分布还包含边界权重

设真实晶界边集合为 E,每条边的长度为 lij,则 Fadj(B)=Σ(i,j)∈E lij1[θij∈B]/Σ(i,j)∈E lij。若按界面条数统计,把 lij 换成 1;若按三维面积统计,应换成真实界面面积。把一条边切成很多段后按段数计数,会使结果依赖软件分段方式。

直接将长度加权邻接分布与晶粒等权非关联分布比较,可能混入大小和连接度的影响。更严格的空间零假设可固定同一邻接图及边长度,在图上打乱取向标签,重复计算邻接分布;若大小与取向已有强关联,还需在合理尺寸组内置换。零假设应写清保留哪些组织特征、打乱哪些关系。

08

六个晶粒的课堂配对

设 6 个同大小晶粒,其中 3 个为取向 A,3 个为取向 B,A、B 的对称约化差为 60°/〈111〉。排除自配对的 15 个无序对中,同类共有 C(3,2)+C(3,2)=6 对,异类有 9 对;非关联分布因此有 40% 的 0° 对和 60% 的孪晶关系对。

若排成一条 A–B–A–B–A–B 的链,5 条邻接边全是异类,邻接关系峰达到 100%;若排成 A–A–A–B–B–B,只有 1/5=20% 为异类。取向列表完全相同,空间排列改变了边界统计。真实晶粒不会靠标签简单拼成这种理想链,这个例子只隔离“配对方式”的数学作用。

09

为什么还要看轴和低角误差

从旋转矩阵提取单位轴可用 a=(M32−M23,M13−M31,M21−M12)/(2sinθ),远离 0 和 π 时成立。θ 很小时分母很小,相同矩阵噪声会被放大,轴方向的不确定性约随噪声/θ增长;这解释了低角轴图常比角图更不稳定。

因此一个 60° 峰不能独立认定 Σ3,必须在规定轴坐标与对称处理中检查〈111〉附近;一个 1° 左右的轴聚集也不能忽略测量精度。按角区间分别画轴分布,能看出同一角峰里是否混合了不同晶体学关系。

10

Randle 黄铜案例应当怎样读

Randle 图 7.8 对退火 α 黄铜同时给出实际邻界分布、非关联配对分布和轴分布。实际邻界在较大角度处有显著峰,轴分析及 CSL 分类将其与退火孪晶联系;非关联分布则更接近随机基线。这里的论据来自三个视角相互配合,并非看到一个柱子高就直接命名。

实操阅读顺序是:先确定相和对称群,再确认统计边界与权重,叠加保留织构的基线,检查峰的轴组成,然后回地图看这些边界的位置、形貌和连通。只有到这一步,才有条件讨论某种边界形成或存留的偏好。

来源与延伸阅读

以下链接用于核对概念、标准状态或原始文献。本页为中文教学整理,不复制受版权保护的标准正文。