龙哥觉得，搞第一性原理计算的人都知道，能量对原子位移的泰勒展开引出两个漂亮的不变性：平移不变性和旋转不变性。平移不变性就是整个系统平移一下能量不变，对应的就是声学求和规则（Acoustic Sum Rule, ASR），这个通常在DFT代码里已经修好了。但旋转不变性——系统转一转能量也不变——却常常被忽略。为啥？因为修起来麻烦，而且以前大家觉得影响不大，直到研究低维材料时，柔性声学模式的二次色散怎么算都不对，这事儿才被拎到台面上。
从理论到实践：倒易空间中的旋转不变性修正
传统的修正方案是在实空间里对原子间力常数（Interatomic Force Constants, IFCs）加约束，比如Quantum Espresso里用的最小二乘法。但实空间方法有个致命伤：修正会不均衡地施加到最远处的IFCs上，导致傅里叶变换后在倒易空间出现非解析行为，算出来的声子谱反而更糟。图8：用简化模型说明实空间施加旋转不变性的挑战。(a)初始IFCs和修正后IFCs；(b)修正量在实空间随距离线性增长，导致倒易空间非解析。本论文剑走偏锋，直接在倒易空间里动手。核心思想是：旋转不变性条件可以写成对倒易空间IFC矩阵及其波矢导数的约束。具体来说，论文构造了一个包含区中心IFCs、一阶导数和二阶导数（均对波矢）的向量x（见公式截图），然后构建条件矩阵A，使得当x满足平移和旋转不变性时，Ax=0。公式(11)：变量向量x的定义，包含区中心IFCs的实部、一阶导数虚部、二阶导数实部。其中κ,κ′为原子索引，α,β,γ,δ为笛卡尔方向。公式(6)：三维情况下旋转不变性在倒易空间的条件：一阶导数叉乘的ϵ求和为零。ɛβγδ是Levi-Civita张量。论文提出了两种修正方案。第一种是Moore-Penrose伪逆法：计算条件矩阵A的伪逆A⁺，然后对初始x施加最小二乘修正 Δ = -A⁺Ax，使得修正后向量满足所有条件。这种方案通用但计算量稍大。公式(12)：伪逆修正公式。第二种是直接修正法，更接近传统ASR修正的“简单粗暴”风格。论文将一阶导数分解成对称和反对称部分（见下面公式），旋转不变性只约束反对称部分。通过引入轴向量形式，反对称部分修正可以直接施加在对角线上，就像修正ASR一样简单。公式(14)：一阶导数分解为对称和反对称部分。公式(18)：轴向量的旋转不变性条件——简单求和为零。公式(19)：直接修正：轴向量的对角元素减去所有非对角元素的和。两种方法修正完区中心及其导数后，需要通过对称化操作将修正传播到整个布里渊区，生成新的IFCs矩阵（见下面公式）。公式(22)：修正后的IFC矩阵表达式。Δ(i)是不同阶的修正项。龙哥看到这个思路，忍不住竖大拇指：倒易空间只处理原胞内原子，不受IFC截断半径影响，计算复杂度与原子数立方成正比，比实空间方法高效得多。
高端的食材往往只需要最朴素的烹饪方式……但是高精度的声子计算却难逃长程静电的“魔爪”。对于半导体和绝缘体，偶极子-偶极子（dipole-dipole）相互作用导致IFC在实空间长程衰减，必须用解析公式处理。常规做法是把IFC分成短程（短程IFC）和长程（静电部分）两部分。公式(30)：偶极-偶极贡献的解析表达式，包含Born有效电荷Z*和介电常数张量ε∞。论文发现一个让人头疼的事实：即使DFT计算的短程部分通过伪逆法完美修正了旋转不变性，但长程静电部分本身就不是旋转不变的！因为静电贡献依赖于介电常数张量，当材料具有各向异性介电张量时（绝大多数低维材料都是），旋转系统不旋转介电张量，能量就会变。更具体地说，传统偶极-偶极公式（公式30）在介电张量各向异性时，即便在q=0处也会产生非零的一阶矩，破坏旋转不变性。图1：各向同性（左）和各向异性（右）介电环境中点电荷静电势等势面。各向异性情况下，旋转系统需要同时旋转介电张量才能保持能量不变。论文在磷烯（phosphorene）上定量验证了这一现象。TABLE I对比了IFC实空间矩和长波方法计算的导数，发现二阶导数条件（即公式9）破坏严重：表I：磷烯中一阶和二阶IFC导数与实空间矩的对比。一阶条件基本满足，二阶条件偏差很大（如(x,x)列：0.1946 vs 0.0671 GPa），说明长程静电贡献导致二阶导数不满足旋转不变性。由于这个原因，论文目前只对IFC的一阶导数施加旋转不变性修正，二阶导数暂时保留原样。这虽然能恢复柔性模式的二次色散，但二阶部分的残留误差对垂直于平面的模式（如ZA）仍有一定影响。
[1] B. Van Troeye, X. Gonze, G. Pourtois, "Reciprocal-space implementation of rotational invariance for first-principles phonons and application to low-dimensional materials", arXiv:2606.16320 (2026).[2] ABINIT software package, https://www.abinit.org[3] P. Giannozzi et al., "QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials", J. Phys.: Condens. Matter 21, 395502 (2009).[4] X. Gonze et al., "Interatomic force constants from first principles: the case of the acoustic sum rule and rotational invariance", Phys. Rev. B 89, 045205 (2014).
