说明：本文华算科技主要介绍材料计算中为什么需要先建立晶胞，晶胞如何把无限晶体变成可计算的周期模型，以及晶胞参数、原子坐标、边界条件怎样影响能量、电子结构和后续判断。
晶胞先确定的到底是什么？
第一性原理计算面对的不是一个材料名称，而是一组明确的原子核位置、电子数和边界条件。晶胞的作用，是把宏观晶体压缩成一个可重复的基本计算单元。晶格矢量给出周期平移规则，原子坐标给出单元内部的基元，元素种类和价电子设置决定赝势与电子数，三者合起来才形成求解 Kohn-Sham 方程的对象。
如果只说“算 MoS2”或“算 TiO2”，计算程序并不知道材料采用哪种相、哪种空间群、哪一个晶面，也不知道原子之间的邻接关系从哪里开始重复。晶胞先把无限周期结构变成有限变量：三条晶格矢量、若干分数坐标、元素占位和对称信息。没有这些变量，哈密顿量中的外势项就没有清楚的空间排布。
分数坐标尤其能说明晶胞的意义。原子位置不是孤立写在笛卡尔空间里，而是绑定在 a、b、c 三个晶格方向上。晶胞拉伸、剪切或角度改变时，分数坐标相同的原子会对应新的真实键长和键角。晶胞不是画图外框，它直接参与定义近邻关系、配位环境、体积、应力和周期镜像距离。
图1. 周期晶体的几何描述符、晶胞参数和空间群信息对比。DOI：10.1038/s41598-024-59938-z
周期边界还带来一个容易忽略的事实：计算盒子外面不是空白，而是该晶胞在三个方向上的无限复制。对于体相晶体，重复方向决定长程周期的近似方式；对于二维材料和 slab，某些方向需要加真空层，把物理周期和数值边界区分开，后面的模型判断才有坐标基准。
总能量、DOS、声子谱和吸附能都沿着同一晶胞延拓来定义。总能量对应这个单元及其周期重复，DOS 来自这个周期势场下的电子态，声子谱来自该结构附近的力常数，吸附能绑定在某个表面单胞和覆盖度上。晶胞给出的不是几何准备，而是材料计算的数学边界。
为什么没有晶胞就没有可比较的能量？
能量比较必须指向同一类计算对象。形成能、相对能、吸附能、缺陷形成能看起来都是能量差，但每个能量项背后都有晶胞、原子数、体积、对称和参考态。晶胞把“同一组成”细化为“某一种具体结构”，随后能量才有可比较的基准。
二维材料里，同样的元素比例可能形成平面蜂窝结构、起伏结构、链状结构或含多元环结构；体相材料里，同一组成也可能有不同晶型。DFT 给出的低能构型不是由化学式单独决定，而是由晶胞内原子排列和周期重复方式共同决定。
图2. 二维材料结构搜索中由对称生成、晶体向量表征到 DFT 优化的闭环流程。DOI：10.1038/s41524-023-01122-4
结构搜索流程之所以要先生成候选晶胞，就是因为势能面上存在许多局部极小值。一个候选模型若晶胞太小，可能排不下真实低能配位；若对称约束过强，某些弛豫方向会被人为锁住；若初始键长异常，结构优化会从错误区域开始移动。能量优化是在既定晶胞自由度内寻找局部低能结构，它不会替代前面的模型定义。
图3. 多种二维化合物的低能结构搜索结果和单胞内原子数约束。DOI：10.1038/s41524-023-01122-4
能量数值必须说明“每什么单位”。eV/atom、eV/formula unit、单位面积能量和吸附前后能量差对应的归一化对象不同；缺陷形成能还会加入化学势和电荷态修正。晶胞给出原子数、面积、体积和周期镜像距离，归一化方式随之确定。
在 B-C-N 体系中，局部配位和多元环结构会改变形成能，不同组成点可以形成不同低能结构；sp2 配位、B-N 配对或 C/N 局部环境都会改变能量排序。这里参与比较的是原子排布后的周期模型，而不是几个元素符号本身。
图4. B-C-N 体系的形成能相图、候选结构和结构特征聚类。DOI：10.1038/s41524-023-01122-4
这就是材料计算里“先建模型”的根本原因。晶胞、坐标和约束条件先把问题限定清楚，能量、力、应力和电子态才有求解对象。若后面发现结果异常，常常不是某个输出数字突然失效，而是前面的晶胞、自由度或参考态没有对应目标材料状态。
晶胞怎样决定倒易空间、电子结构和图谱？
电子结构图谱同样从晶胞出发。实空间晶格矢量确定倒易晶格，倒易晶格再确定布里渊区、k 点采样和能带路径。能带图上的 Γ、M、K、X 等高对称点并不是通用坐标标签，而是当前晶胞和空间群在倒易空间中的特殊位置。晶胞换了，路径、能带折叠和能级简并关系都可能改变。
DOS 和能带也绑定在同一个周期势场里。晶胞内的配位数、键长、键角和局域结构单元会改变轨道重叠，轨道重叠再改变价带、导带、带隙和费米能级附近态密度。电子结构来自具体晶胞势场中的本征态集合，材料名称只提供化学背景，不能替代模型本身。
图5. B-C-N 体系带隙分布、配位特征和能带边态空间分布。DOI：10.1038/s41524-023-01122-4
B-C-N 体系中的N-C-C-N 链、B-N 多元环和不同配位 C 原子会改变带边态分布；相同元素组成下，若晶胞内连接方式变化，带边态的空间分布和 DOS 峰形也会改变。带隙大小要回到晶胞内的具体配位来解释。
二维材料和表面模型中的真空方向长度、面内晶格和 slab 厚度分别控制层间镜像距离、应变状态以及表面与内部层的分离程度。若晶胞法向太短，层间相互作用会进入本应隔开的真空区；若面内晶格被强行匹配，能带、功函数和吸附能可能带着额外应变。
图6. 多种低能 B-N 基结构的能带和 DOS 结果。DOI：10.1038/s41524-023-01122-4
能带和 DOS 图提供的是当前晶胞下的电子态分布。图中带隙、峰位、峰高和投影贡献可以支持电子结构判断，但它们不自动代表所有晶型、所有应变状态或所有表面终止。晶胞、路径、k 点密度和投影设置共同限定图谱含义，否则不同模型之间的图谱差异会被误读为材料本征差异。
建晶胞时最该守住哪些物理边界？
体相、二维材料、表面 slab 和分子盒子对应不同维度边界。周期方向和非周期方向必须与目标问题一致：体相材料需要三维周期；二维材料通常保留面内周期并加入足够真空；表面 slab 需要描述外露晶面、层数和真空区；分子大盒子只是用非相互作用镜像近似孤立环境。
结构优化的自由度决定结果对应哪一种材料状态：体相结构优化常允许晶胞体积和形状弛豫，固定压力条件下需要检查应力；表面 slab 往往固定部分底层原子，并且不随意释放真空方向晶格；二维异质结会面临晶格失配，过大的强制应变会改变能带和界面电荷。哪些晶格方向能动，哪些原子受约束，决定结构优化结束后的模型含义。
晶胞尺寸越大，能容纳的缺陷、吸附物、重构和长程扰动越丰富，但电子步、离子步、k 点和后处理成本也会上升。小晶胞计算快，却容易放大周期镜像、覆盖度或缺陷浓度；大晶胞更接近稀释极限，需要用收敛测试确认能量差已经稳定。
图7. 不同结构搜索方法在时间成本和最低能结构上的对比。DOI：10.1038/s41524-023-01122-4
参考状态需要和晶胞记录绑定：吸附模型要记录表面晶面、覆盖度和吸附位点；缺陷模型要记录超胞大小、化学势、电荷态和电荷修正；声子模型要记录力收敛、位移超胞和 q 点路径。同一个输出数值若离开晶胞记录，科学含义会迅速变窄甚至改变。
计算前先建晶胞，并不是软件流程里的形式步骤。它把材料的周期、局域配位、边界条件和可放松自由度写进模型；后面的总能、力、应力、能带、DOS、声子和吸附能，都只是这个模型在不同物理量上的投影。晶胞建得对，计算结果才有对应的材料对象；晶胞建得不对，再漂亮的曲线也只能描述那个错误模型。