🔬 materials science

The crystallographic quaternions and their product law

本文通过将四元数乘法法则推广到直接在晶体基底内运行，并利用一个源自度量张量的“交叉张量”，引入了晶体学四元数，从而消除了在计算晶体学旋转复合时对坐标变换的需求。

原作者： Cyril Cayron

发布于 2026-07-21

📖 1 分钟阅读☕ 轻松阅读

原作者： Cyril Cayron

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一下，你正试图描述一个玩具陀螺如何旋转，或者一名舞者如何进行旋转跳跃。在物理学和工程学的世界里，我们有一个特殊的数学工具叫做“四元数”（quaternion）。你可以把四元数想象成一本极其高效的旋转指令手册。与其使用一组臃肿且复杂的数字（比如矩阵）来告诉计算机如何旋转物体，四元数使用一种紧凑的代码，计算速度更快，且不易出错。通常，这些指令是基于一个完美的方格网格构建的世界——就像坐标纸一样，每一条线都与另一条线呈完美的 90 度角。这被称为“笛卡尔基底”（Cartesian basis），也是大多数数学运算的标准游乐场。

然而，现实世界的晶体并不总是由完美的正方形构成的。许多晶体（例如金属或矿物中的晶体）其原子排列呈现出奇特的倾斜形状。它们的内部网格可能是一个拉伸的长方形、一个倾斜的盒子，或者是一个六边形。如果你尝试用标准的“方格网格”旋转指令来处理这些倾斜的晶体，你就必须不断地在倾斜的晶体世界和虚构的正方形世界之间进行数学转换，才能得到答案。这就像是用方形的模具烤蛋糕，但你却不得不先在圆碗里测量配料，然后再转换回方块单位，接着再测量一次。这样做虽然可行，但既混乱、缓慢，又容易出错。这就是从事电子背散射衍射（EBSD）研究的科学家们一直面临的难题——EBSD 是一种通过拍摄晶体表面图像来观察其取向的技术。他们需要精确知道这些倾斜的晶体是如何旋转的，但标准的数学工具感觉像是说错了语言。

于是，来自洛桑联邦理工学院的研究员西里尔·凯伦（Cyril Cayron）决定停止这种“翻译”，直接开始使用晶体的母语进行交流。在他的论文中，他提出了对四元数系统的一次巧妙升级。他主张我们不需要为了做数学题而强行将晶体塞进一个正方形的盒子里。相反，我们可以推广四元数的规则，使其能够完美地运行在晶体实际存在的那些倾斜、形状怪异的网格上。

这种新方法的“秘密武器”是凯伦引入的一对数学“魔杖”。第一个工具是科学家们已经熟知的度规张量（metric tensor），你可以把它想象成一把知道如何在倾斜表面上测量距离和角度的尺子。第二个是凯伦发明的新工具——交叉张量（cross tensor）。如果说度规张量是尺子，那么交叉张量就是那个知道如何在倾斜网格上处理“叉积”（一种寻找两个方向的垂直方向的方法）的量角器。凯伦展示了这两个工具实际上是亲戚：交叉张量只是度规张量逆矩阵的一个特定的、经过计算的版本。

通过结合这些工具，凯伦推导出了一个新的四元数乘法公式。这个新公式允许科学家直接在晶体自身的倾斜坐标系内计算两种旋转是如何结合的（比如先旋转陀螺，然后再倾斜它）。不再需要为了回到正方形网格而反复翻译了！论文表明，对于非立方晶体，这种直接方法不仅在数学上是严谨的，而且在计算机程序中也可能更加高效。

该论文还解决了一个常见的误区：即认为你必须使用笛卡尔（正方形）基底来测量晶体中的角度或距离。凯伦认为这种观念是不必要的。虽然旧有的使用“罗德里格斯向量”（Rodrigues vectors，另一种旋转代码）的方法非常紧凑，但它有一个致命缺陷：如果你尝试旋转正好 180 度，它会变成无穷大，从而导致计算机程序崩溃。四元数避开了这个无穷大的问题，但此前却一直受困于需要向正方形网格进行转换。

凯伦通过计算来验证他的“晶体学四元数”是否真的更快。他发现，对于单次旋转计算，新方法大约需要 46 个数学步骤（浮点运算），而旧方法（先转换到正方形网格再转回）则需要大约 43 个步骤。乍看之下，旧方法似乎略快一些。然而，论文指出这种比较是不公平的，因为它忽略了隐藏的成本。如果算上在开始时将数据转换到正方形网格，以及在结束时转换出来的额外步骤，旧方法实际上做了更多的功。此外，新方法使得识别对称性（例如当晶体在旋转后看起来依然一样时）变得更加容易，因为数字在晶体自身的语言中保持得简单且纯净。

简而言之，这篇论文建议我们可以停止强迫晶体进入正方形的盒子。通过引入“交叉张量”并更新四元数的乘法规则，我们可以直接在晶体天然的倾斜网格上进行旋转数学运算。这使得分析电子背散射衍射中的晶体取向的过程变得更加清晰、优雅，并且在处理构成我们物质世界的复杂非立方晶体时，效率可能更高。作者总结道，这种方法可以取代目前用于此类分析的计算机程序中的标准方法，为理解微观世界的旋转提供一种更直接、更稳健的方式。

技术摘要：晶体学四元数及其乘法法则

问题陈述
单位四元数是科学领域中用于编码旋转的标准工具，因其与矩阵或罗德里格斯矢量（Rodrigues vector）方法相比具有更高的计算效率和稳定性。然而，标准的四元数公式本质上是与笛卡尔基底绑定的。在晶体学中，特别是对于非立方晶系，晶体基底是非笛卡尔的。电子背散射衍射（EBSD）等领域的现有实践需要使用“结构张量”在晶体基底与任意附加的笛卡尔基底之间进行转换，以执行旋转组合。这种在坐标系之间来回转换的过程被认为效率低下，且依赖于并非晶格本身基本属性的约定（即结构张量）。本文旨在解决直接在晶体基底内进行旋转组合而无需经过中间笛卡尔转换的需求。

方法论
作者提出了一种广义的四元数乘法法则，以适应非笛卡尔晶体基底。该方法涉及三个主要的理论步骤：

重构四元数定义： 旋转四元数相对于晶体基底 ( $\mathcal{B}_c$ ) 被重新定义。向量 $\mathbf{u}$ 不再由标准的欧几里得范数归一化，而是通过度量张量 $\mathcal{M}$ 进行归一化，使得 $\|\tilde{\mathbf{u}}\| = \sqrt{\tilde{\mathbf{u}}^t \mathcal{M} \tilde{\mathbf{u}}} = 1$ 。
引入叉积张量： 虽然度量张量 $\mathcal{M}$ 在计算非笛卡尔基底下的标量积时已广为人知（ $\mathbf{u} \cdot \mathbf{v} = \mathbf{u}^t \mathcal{M} \mathbf{v}$ ），但文献中尚不存在用于叉积的等效张量。作者引入了“叉积张量”（ $\mathcal{X}$ ），其定义为由晶体基底向量的叉积（ $\mathbf{b} \times \mathbf{c}, \mathbf{c} \times \mathbf{a}, \mathbf{a} \times \mathbf{b}$ ）所构成的矩阵。
推导广义乘法法则： 作者证明了叉积张量与度量张量的逆成比例，并受单位晶胞体积 ( $V$ ) 的缩放，即 $\mathcal{X} = \sqrt{\det(\mathcal{M})} \mathcal{M}^{-1}$ 。这种关系使得标准的四元数组合公式（依赖于点积和叉积）可以利用 $\mathcal{M}$ 和 $\mathcal{X}$ 进行重写，从而实现在整个晶体基底内进行四元数乘法。

核心贡献

叉积张量： 主要的理论贡献是叉积张量 $\mathcal{X}$ 的定义与推导。该张量允许使用与向量相同的坐标表示来计算非笛卡尔基底下的叉积，类似于度量张量处理标量积的方式。
晶体学四元数乘法法则： 本文提供了两个晶体学四元数之积（ $\mathbf{q}_3 = \mathbf{q}_1 \mathbf{q}_2$ ）的显式公式。标量部分使用度量张量进行点积，向量部分使用叉积张量进行叉积。
罗德里格斯矢量的推广： 本文将罗德里格斯矢量的组合法则扩展到晶体学语境中，提供了一个同样依赖于 $\mathcal{M}$ 和 $\mathcal{X}$ 的公式。
向量旋转公式： 给出了一个在晶体基底中绕轴 $\tilde{\mathbf{u}}$ 旋转角度 $\alpha$ 的向量 $\mathbf{v}$ 的广义公式，其中用叉积张量操作取代了标准叉积。

结果与计算分析
文中包含了计算成本分析（以浮点运算次数 FLOPs 衡量），将新的晶体学四元数方法与传统的笛卡尔四元数方法进行了对比：

直接乘法成本： 计算两个晶体学四元数的乘积需要 46 FLOPs（28 次乘法，18 次加法），高于笛卡尔四元数所需的 28 FLOPs。
系统性效率： 作者认为，公平的比较必须包含传统方法的开销。在晶体基底与笛卡尔基底之间转换旋转轴需要结构张量（增加约 15 FLOPs）。此外，初始获取 EBSD 数据（将菊池花样转换为欧拉角）通常涉及与结构张量的矩阵乘法（增加约 45 FLOPs）。
净优势： 当考虑完整的流程时，晶体学四元数方法完全避免了这种“来回”转换。所有计算都保留在晶体基底内，从而可能降低整个 EBSD 处理流水线的总计算负荷。
对称性处理： 文中指出，使用晶体学四元数处理对称操作更加优雅且计算更简便，因为晶体基底中的对称轴通常具有整数坐标（0, ±1），而在笛卡尔表示中则涉及无理数。

意义与主张
本文声称，引入叉积张量和广义乘法法则消除了在晶体学旋转计算中使用任意笛卡尔基底和结构张座的必要性。作者建议，晶体学四元数可以有利地取代笛卡尔四元数，并应用于 EBSD 软件中。通过直接在晶体基底内工作，该方法与衍射数据的本质（这些数据通常通过球面索引直接从晶体基底中获得）相一致，并消除了不必要的坐标转换。本文认为，这种方法不仅在理论上更具根本性（依赖于度量张量而非依赖于约定的结构张量），而且对于非立方晶系取向图的完整处理过程而言，在效率上也具有潜在优势。

技术摘要：晶体学四元数及其乘法法则

类似论文