说明：本文华算科技主要介绍晶格常数的定义、它在第一性原理计算中的来源，以及晶胞参数变化怎样影响结构模型、倒易空间和后续材料性质分析。
晶格常数到底定义了什么结构尺寸？
晶格常数又称晶胞参数，指周期晶体重复单元的几何尺寸。三维晶胞由三条晶格矢量描述，矢量长度写作 a、b、c，矢量夹角写作 α、β、γ。六个参数确定晶胞形状与体积，原子坐标在这一几何基底中以分数坐标或笛卡尔坐标记录。
立方晶体的三个边长相等、三个夹角均为 90°，一个 a 值即可描述晶胞尺度。四方晶体由 a 和 c 给出，六方晶体常用 a、c 与 γ=120° 记录。低对称晶体中，a、b、c 以及三个夹角可各不相同；缺少 b、c 与夹角时，晶胞几何无法还原。
图1. 单层 M X 化合物的晶体结构与原胞示意，包含晶格常数、原胞选择和内部原子坐标三类信息。DOI：10.1038/s41598-019-49890-8
把晶格常数当作最近邻原子距离，容易在含多原子基元或低对称晶体中误读结构。少数简单结构中，两者近似对应。对于金刚石、钙钛矿、层状材料或含多原子基元的晶体，最近邻键长由晶格参数和内部原子坐标共同决定。晶格常数定义平移周期，并非某一根化学键的长度；键长、键角、配位数来自晶胞内原子位置的几何计算。
同一周期结构在计算文件中可采用原胞、常规胞或超胞三种写法，边长数值必须按晶胞类型对齐。原胞体积最小，常规胞用于表达对称性，超胞用于缺陷、掺杂、表面或声子计算。晶格常数必须和所选晶胞写法绑定，把原胞 a 值和常规胞 a 值混合比较，容易把同一周期结构误读成不同材料。
晶格常数在第一性原理中从哪里来？
第一性原理计算中的初始晶格常数多来自实验晶体结构、数据库条目、文献模型或人工构造的周期模型。初始值只是结构优化的起点；晶胞自由度开放时，程序根据总能量、原子受力和应力张量寻找局部平衡结构。此时得到的 平衡晶格常数对应给定泛函、赝势、磁态、压力和约束条件下的几何最小点。
体相材料多通过晶胞优化或状态方程拟合获得平衡体积。不同体积下的总能量曲线在最低点附近给出平衡体积，随后换算成晶格常数；曲线曲率与体模量有关。总能量最低点给出当前理论模型对平衡结构的预测，并非实验值复刻。
图2. Fe 的 DFT 总能量随体积变化曲线，最低点对应平衡体积和晶格参数。DOI：10.1002/qua.24675
晶胞弛豫的充分程度由原子受力和应力张量共同约束，最后一步总能量变化很小只反映能量迭代趋于平缓，晶胞几何是否平衡仍由应力收敛给出。原子力约束晶胞内部原子位置，残余应力反映晶胞边长或夹角仍承受拉伸、压缩或剪切趋势。在固定晶胞条件下只放松原子坐标，结果属于固定晶格下的局部结构；开放晶胞自由度时，a、b、c 和夹角随应力释放而改变。
平衡体积对泛函、赝势、色散修正和自旋设定敏感，同一材料换用交换相关泛函或赝势后，优化得到的 a、b、c 常出现系统偏移。LDA 倾向于较小体积，部分 GGA 给出偏大的晶格常数，PBEsol、SCAN 或 vdW 修正会改变误差方向。晶格常数带有方法依赖，来自结构模型与电子结构近似的耦合结果。
固定晶胞、只放松原子位置、各向同性改变体积、允许晶胞形状完全弛豫，这几种设置对应不同结构含义。层状材料的面内 a 值和层间 c 值宜分开评估；表面 slab 的真空层属于模型留白，不应并入材料本体的 晶胞自由度。
图3. 不同赝势数据相对全电子结果的晶格常数误差，对应离子势处理带来的平衡结构偏差。DOI：10.1016/j.commatsci.2013.08.053
晶格常数变化会牵动哪些计算对象？
晶格常数在优化、应变加载或外压条件下改变时，实空间几何直接体现在晶胞体积、密度和局域配位上。体相晶体的晶胞体积、密度、键长、键角和配位多面体随几何约束变化；二维材料的面内 a 值改变对应拉伸或压缩状态，层状材料的 c 值和层间距影响层间弱相互作用。实空间尺度改变后，电子云重叠、轨道杂化和力常数随结构重新分布。
晶格常数进入倒易基矢定义后，能带、声子和 k 点采样使用同一套倒易空间尺度。简单晶格中，实空间周期增大，对应倒易基矢长度减小；实空间周期缩短，对应倒易空间范围增大。k 点密度、能带路径长度、声子 q 点采样和超胞折叠都由倒易基矢标定。相同的 KPOINTS 网格用于不同大小晶胞时，实际倒易空间分辨率不同。
图4. 广义 Monkhorst-Pack 网格在倒易空间中的 k 点分布，k 点间距随晶胞尺度设定。DOI：10.1016/j.commatsci.2020.110100
弹性、声子、能带和表面计算都把晶格常数当作几何基准。弹性常数来自小应变下的能量或应力响应，声子频率来自平衡结构附近的力常数，能带色散由晶格周期下的 Bloch 态给出，表面 slab 的真空层厚度和覆盖度常用晶胞面积表示。晶胞面积、体积和倒易基矢改变后，后续图谱的坐标基准随之变化。
图5. 石墨烯、MoS2 和磷烯中的缺陷与吸附超胞模型，对应超胞面积、缺陷间距和覆盖度定义。DOI：10.1038/s41598-018-27632-6
在缺陷和掺杂模型中，晶格常数决定超胞体积、缺陷之间的周期镜像距离和名义浓度。一个 2×2×2 超胞与 4×4×4 超胞来自同一原胞，缺陷浓度、弹性相互作用和电荷修正仍然不同；表面模型中的 晶胞面积决定覆盖度定义中的面积基准，并影响吸附能和表面自由能的数值参照。
晶格常数变化与具体结构变量相连。晶格常数增大时，部分键长增加，部分体系主要改变空隙体积或层间距；带隙、磁矩、吸附能和扩散能垒的变化方向并不统一。性质分析应比较具体键长、配位、态密度、声子模式或反应路径，把 a 值变化直接写成性能变化，会掩盖真正发生改变的结构对象。
比较晶格常数时哪些条件必须同场出现？
晶格常数是结构比较中最常被引用的数值，温度、压力、相结构和晶胞写法需保持一致。实验晶格常数多来自有限温 XRD、电子衍射或中子衍射，静态 DFT 优化多对应 0 K 附近的势能面局部极小点。温度和压力改变时，热膨胀、相变和应力状态会改变测得或算得的晶胞参数。
泛函或赝势的比较要按目标性质分开处理。某个泛函给出的晶格常数接近实验，只能表明它在该结构尺度上误差较小；带隙、形成能、缺陷电荷态或磁序能量差对应不同误差来源。目标性质不同，结构误差、电子自相互作用误差和色散相互作用误差应分别评估。
图6. 不同密度泛函给出的晶格常数和体模量误差分布，呈现结构尺度和力学响应中的方法误差。DOI：10.1002/qua.24675
磁性体系、强关联氧化物、层状材料和含缺陷超胞必须标明额外条件。磁态和电荷态改变局域键长，U 值影响 d/f 电子局域程度，范德华修正显著影响层间距，缺陷或掺杂浓度则由超胞尺寸决定。缺少这些条件时，同一个 a 值可能对应不同磁态、电荷态、层间距或缺陷浓度。
图7. α-S8 体积随温度变化的 XRD 与 QHA 结果，对应实验晶格参数与 0 K 优化结构的温度参照差异。DOI：10.48550/arXiv.1609.07980
晶格常数与空间群、原子分数坐标、优化约束、泛函、赝势、色散修正、磁态和温度参照一起给出，结构状态才明确。a、b、c 和夹角对应具体周期模型、平衡条件和目标性质；能带、声子、吸附或缺陷结论依托相应结构状态和计算证据。
多个材料或多个构型并列比较时，a、b、c 的变化需与空间群、原子位移和体积同步读取。晶格常数变化很小而内部原子坐标明显移动时，局部键长和配位多面体仍可改变；晶胞体积变化明显时，电子结构、声子频率和弹性响应随几何基准重排。结构基准确认后，能带、DOS、声子或吸附结果才具备明确的几何参照。