正格子、倒格子与布里渊区的几何关系及其在第一性原理计算中的应用

说明：本文华算科技主要介绍正格子、倒格子与布里渊区的几何关系，以及它们在第一性原理计算中的不同用途。

正格子、倒格子和布里渊区是什么

晶体的周期平移可写成 R = n₁a₁ + n₂a₂ + n₃a₃，其中 n_i 为整数，a_i 为三条平移基矢。所有 R 构成正格子，也称实空间格子或直接格子。它记录周期平移的位置集合；原子种类和原胞内坐标属于基元，二者共同给出具体晶体结构。

倒格子的基矢 b_j 由 a_i 唯一确定，并满足 a_i·b_j = 2πδ_ij。离散矢量 G = hb₁ + kb₂ + lb₃ 构成倒格子，其单位通常为 Å⁻¹。任意 G 与正格子平移相乘都有 e^iG·R = 1，这使倒格子适合描述周期函数的傅里叶分量、衍射条件和周期势中的平面波。

布里渊区取自倒易空间。以某个倒格点为中心作 Wigner–Seitz 构造：连接它与周围倒格点，在连线中点作垂直平分面，由这些平面围出的区域就是第一布里渊区。第一布里渊区是倒格子的一个原胞，它是有限区域；倒格子则是向整个倒易空间延伸的离散点阵。

Bloch 波矢 k 可在倒易空间连续取值，G 是离散倒格矢。由于 k 与 k + G 在正格子平移下具有相同相位标签，能带可折回第一布里渊区。Γ、X、M、K 等符号标记区内的高对称点；它们依赖晶格类型和采用的原胞约定，不是任意材料都共享同一组坐标。

图2中的灰色矩形围出第一布里渊区，蓝色 Γ–X–M–Y 区域利用点群对称性缩减了独立波矢范围。不可约布里渊区只保留对称不重复部分，它与第一布里渊区的几何大小不同，程序通过对称操作和权重恢复全区积分。

正格子怎样确定倒格子的尺度和方向？

三维原胞体积 Ω = a₁·(a₂ × a₃)。采用含 2π 的定义时，b₁ = 2π(a₂ × a₃)/Ω，其余两条循环置换。这个叉乘关系给出两点：每条倒易基矢垂直于另外两条实空间基矢张成的平面，其长度受对应晶面间距控制。

正格子原胞与倒易原胞的体积满足 ΩΩ^* = (2π)³。某一实空间周期增大时，与其共轭方向相关的倒易尺度减小；剪切或改变夹角时，倒易基矢的方向和夹角随之重排。倒易空间保留完整的度量信息，不只是一幅把长度取倒数的示意图。

蜂窝状原子排布常让人把六边形当作正格子的最小平移单元。石墨烯的原胞由两条平移基矢和 A、B 两个碳原子组成，两原子基元产生蜂窝结构；相应倒格子为三角点阵，围绕 Γ 点的 Wigner–Seitz 原胞呈六边形。K 点处的 π 与 π^* 能带相接，反映的是该材料的电子结构，不是六方布里渊区必然具有的能谱。

采用原胞时，第一布里渊区具有最大的周期等价区域；扩胞引入更长的实空间平移，倒易原胞随之缩小，原有能带折叠到较小区域。能带条数增加与布里渊区缩小是一组几何伴随结果，缺陷、磁序和界面模型中的超胞计算经常出现这种情况。

硅烷与石墨烯都可采用六方二维原胞，原胞内的元素、起伏构型和轨道组成却不同。两者可使用相同名称的 Γ、K、M 点，高对称点名称来自晶格对称性；带隙、带边曲率和轨道杂化则由周期势及基元决定。

布里渊区为何只是倒格子的一个原胞？

倒格子可用任意一组原始倒易基矢围成平行多面体原胞，也可用最近邻垂直平分面围成 Wigner–Seitz 原胞。后者就是第一布里渊区。两种原胞体积相同，外形和对称信息的显露方式不同。布里渊区的区面是倒格点间的等距面，在近自由电子图像中，这些面也是满足 Bragg 散射条件的区域。

第一布里渊区之外仍有第二、第三布里渊区，它们按离 Γ 点的几何关系划分倒易空间。周期哈密顿量允许把区外波矢减去某个 G 映射回第一区，约化区方案把独立 Bloch 标签集中在第一布里渊区。扩展区方案保留连续延伸的波矢，两种画法描述同一周期体系。

高对称路径、均匀网格和不可约区域属于三类计算对象。Γ–K–M–Γ 路径是一维折线，用于读取指定方向的色散；Monkhorst–Pack 网格是一组带权离散点，用于近似全区积分；不可约区域由空间群或点群把对称等价的 k 点归并。一条路径无法代替布里渊区积分，均匀网格也不会自动给出常用高对称线上的连续能带。

图5左侧的六边形与菱形都是倒易空间中的原胞选择，绿色区域为不可约楔形，红点代表 5 × 5 网格经对称归并后的采样位置。右侧另一组特殊点采用不同积分构造。采样点是数值求积节点，其数量、权重和偏移方式会影响总能、力与电荷密度的收敛。

低对称结构的不可约区域往往占据更大比例，Γ 点居中的网格并非适合所有采样密度和晶格对称性。金属费米面附近的占据变化尖锐，通常要求更密的 k 点或合适的展宽；大尺寸绝缘体超胞的布里渊区较小，所需网格常可减少。网格密度应以倒易空间间距和目标物理量的收敛为依据，单独比较各方向点数容易忽略晶胞尺寸差异。

三种几何对象怎样进入DFT计算？

结构优化、缺陷形成能、表面吸附和界面匹配以正格子原胞或超胞承载原子坐标与周期性条件。计算程序读取晶格矩阵后生成倒易基矢，平面波基组写成 e^i(k+G)·r。截断能限定可用的 k + G 平面波集合，它控制基组规模；k 点网格控制布里渊区求积，二者对应不同的数值误差来源。

自洽计算在离散 k 点上求解 Kohn–Sham 方程，并按权重累加电子密度、总能和力。能带计算通常固定已经收敛的电荷密度，再沿选定高对称线求本征值。自洽网格负责全区积分，能带路径负责方向分辨，两套 k 点列表承担的任务不同。

图6左图把硅烷能带限制在 K–Γ–M–K 折线上，右侧三个能量面覆盖整个第一布里渊区。折线适合标认高对称方向的带边和简并，二维分布则能检验带极值是否偏离该路径，并给出不同方向的曲率。有效质量来自带边附近的二阶曲率，只读取一条路径时，结论对应那条路径经过的方向。

声子计算使用波矢 q 标记原子位移扰动的周期性。q 仍位于倒易空间并可折回第一布里渊区，但它标记晶格振动，不承担电子 Bloch 态的标签。电子 k 路径与声子 q 路径可共享同一组高对称点，纵轴分别对应电子本征能量和振动频率。

构建超胞后，实空间复制倍数增大，倒易原胞与第一布里渊区缩小，电子能带和声子支发生折叠。缺陷或界面破坏原胞平移对称后，折叠态还会混合。比较原胞与超胞能谱时需建立波矢映射，带展开可把超胞本征态的谱权重投影回原胞布里渊区，区分几何折叠与真实杂化。

正格子决定原子周期性条件和实空间模型，倒格子给出傅里叶展开、衍射矢量与平面波基组，第一布里渊区承载不重复的 Bloch 波矢及全区积分。三者由同一组平移对称性联结，却分别对应点阵、点阵的互易表示和互易点阵原胞。