群发资讯网

晶格概述:晶胞、基矢表述及实空间、倒易空间与电子结构的关联

说明:本文华算科技主要介绍晶格的基本含义、晶胞和基矢在材料计算中的写法,以及晶格怎样连接实空间结构、倒易空间和电子结构图

说明:本文华算科技主要介绍晶格的基本含义、晶胞和基矢在材料计算中的写法,以及晶格怎样连接实空间结构、倒易空间和电子结构图谱。

很多人第一次接触晶体结构时,会把晶格直接当成原子排布图;在材料计算里,这种混用很容易让晶胞、坐标、空间群、k 点路径和结构优化结果挤在一起。晶格先描述周期平移规则,原子种类和原子位置再把这种规则填成具体材料。把这个层次分开,能更自然地理解 CIF 文件、VASP 结构模型、二维材料超胞、缺陷超胞和能带图的共同基础。

这里讨论的“晶格”不是某个软件里的按钮名称,也不是论文图中的装饰外框;它指向一个更基础的问题:晶体为什么能用有限信息表示无限延展的固体。第一性原理计算无法真的放入无限多个原子,只能用周期边界条件把一个有限晶胞复制到整个空间。晶格正是这套复制关系的几何语言。

一、晶格到底是什么:先把基础概念分清

晶格最容易被误读的地方,是把它和真实的原子球棍图混在一起;在晶体学和第一性原理计算里,晶格是由平移矢量生成的周期点阵。若用三个基矢 a1、a2、a3 描述三维晶体,那么任意一个等价晶格点都能写成 R = n1a1 + n2a2 + n3a3,其中 n1、n2、n3 为整数。这个式子没有指定原子是什么,也没有指定某个原子在晶胞里的位置,它只说明“往哪些方向平移多少距离后,环境重新等价”。

图1. 二维周期结构中的晶格矢量、边长和夹角示意,适合用来区分平移规则与图案本身。DOI:10.1038/s41597-025-06150-x

完整晶体结构还应在晶格点上放置基元。基元可以是一个原子,也可以是多个原子、分子片段或更复杂的局域结构。晶体结构 = 晶格 + 基元,这句话背后有一个很实际的含义:同一种平移框架可以承载不同化学组成,不同基元也能放进相似的几何周期中。只看晶格参数 a、b、c、α、β、γ,往往不能判断材料的成键方式和电子结构。

晶胞则是从无限晶格中截取出来的有限重复单元。原胞含有最少的晶格点,常规胞更重视晶体对称性和读图习惯,超胞把多个原胞拼在一起以容纳缺陷、吸附物、磁序或长周期畸变。晶胞是晶格的一种书写方式,同一个周期结构常常可以换用不同晶胞表示;只要平移关系、基元和坐标变换保持一致,描述的物理对象仍然相同。

二维材料、表面模型和体相晶体在这个层面并没有断开。二维晶格只在平面内周期重复,第三个方向通常加入真空层以削弱周期镜像之间的相互作用;表面模型在两个方向保留平移,垂直方向变成 slab 和真空层组合。周期方向的个数决定模型边界,这也是同一套晶格语言能够覆盖体相、薄膜、界面和吸附体系的原因。

图2. 金刚石结构、Wyckoff 操作和表格化结构表达的对应关系,用于说明“晶格 + 基元 + 对称操作”共同定义晶体。DOI:10.1038/s41524-025-01931-9

理解晶格时还应区分“平均周期”和“局域结构”。实际材料中可能存在热振动、缺陷、应变、位错和短程有序,实验晶体结构常给出平均意义上的周期排布;DFT 计算中输入的晶格和原子坐标,是对研究对象作出的模型化选择。晶格给出周期骨架,局域配位、键长、键角、电荷分布和磁序才决定许多具体性质。

二、晶胞、基矢和空间群怎样写进 DFT 模型

第一性原理计算读入晶体结构时,通常看到的是一个三维晶胞矩阵和一组原子坐标;晶胞矩阵的三列就是 a1、a2、a3,它们共同给出晶格长度、夹角和体积。分数坐标依附于晶胞基矢,例如 (0.25, 0.25, 0.25) 并不直接等于某个固定的笛卡尔距离,它表示沿三个基矢分别走四分之一。晶胞改变后,分数坐标的几何含义随之改变。

计算程序把这个有限晶胞沿基矢方向无限复制,电子密度、波函数和势场在平移后满足周期关系。体相材料能用一个晶胞描述宏观晶体,根源就在这里。周期边界条件把有限结构扩展为无限固体,而晶格决定了镜像之间的相对位置。对于缺陷或吸附体系,超胞越大,缺陷镜像或吸附物镜像之间的距离越远,模型更接近孤立扰动。

图3. 基于晶胞尺寸、空间群和 Wyckoff 位置的晶体结构表达,说明计算模型中的晶格并不只是三维外框。DOI:10.1038/s41524-025-01940-8

空间群把平移、旋转、反演、镜面和螺旋轴等对称操作组合起来,告诉我们哪些原子位置可以互相生成。Wyckoff 位置则把对称等价位置写成更压缩的形式。晶格参数和空间群共同限制原子自由度,这使得结构数据库、原型识别和材料生成模型能够用较少信息表达大量晶体。对 DFT 而言,对称性还能减少不可约 k 点和等价原子数量。

图4. 晶体原型识别中,晶格参数、空间群和内部自由度共同参与结构分类。DOI:10.1038/s41524-020-00483-4

原胞和常规胞之间的差别,在材料计算中经常带来表面上的混乱。原胞更小,计算代价低;常规胞更贴近晶体学分类,晶轴方向和对称元素更直观。单胞选择会改变输入文件外观,却不必然改变材料本身。真正要保持一致的是晶格变换、原子坐标变换、k 点采样和性质归一化方式。

图5. 原胞与常规胞之间的几何投影关系,说明同一周期结构可以有不同晶胞写法。DOI:10.1038/s41524-023-01148-8

结构优化时,晶格也不是固定背景。若允许晶胞形状和体积变化,计算会根据总能、应力张量和外加压力调整 a、b、c 与夹角;若只优化原子位置,晶格框架保持不变,体系只能在给定晶胞内寻找较低能量构型。是否放开晶格自由度,会影响应变材料、相变材料、二维材料和多铁材料的结构结果。

材料数据库中的结构条目常包含晶格常数、空间群、原子坐标、占位和磁性信息。进入 DFT 之前,这些信息还会经过标准化、去重、对称识别和初始磁矩设置。晶格是计算模型的起始几何信息,但真正的计算对象还包含电子数、赝势、交换关联泛函、磁序和边界条件。单独拿晶格参数讨论全部性质,很容易把模型层次混在一起。

三、晶格为什么会连到倒易空间和能带图

晶格定义在实空间,电子能带却常画在倒易空间,这两个空间之间通过傅里叶关系相连。只要实空间具有周期平移,电子波函数就可以按 Bloch 形式写成周期函数与相因子的乘积。倒易晶格描述周期波矢的重复规则,它由实空间晶格基矢决定。实空间晶胞越大,倒易空间尺度越小;超胞计算中能带折叠就是这个关系的直接结果。

布里渊区是倒易晶格中的基本区域,能带图横轴上的 Γ、X、L、K 等点,来自晶体对称性和倒易晶格几何。不同晶系、不同 Bravais 格子、不同标准化单胞会给出不同高对称点标签和路径。能带路径不是普通坐标轴,它是一条穿过布里渊区高对称区域的采样路线,用来观察电子能级沿代表性方向的色散。

图6. 不同晶格和对称类别对应的布里渊区路径选择,说明能带横轴来自倒易空间。DOI:10.1038/s41524-020-00383-7

平面波 DFT 中,波函数展开依赖倒易矢 G,k 点网格负责采样第一布里渊区。晶格形状改变时,倒易格矢、k 点密度和高对称路径都会随之改变。同样的 k 点数字不等于同样的采样精度,因为采样间距要和倒易空间尺度一起看。对于长轴晶胞、二维 slab 或大超胞,k 点设置与晶格长度之间的关系尤其明显。

能带图常被读成“材料导电或绝缘的图”,但它其实承载了更多晶格信息。带隙是电子占据与未占据能级之间的能量间隔,色散斜率反映有效质量,简并和劈裂常与对称性、SOC、磁序或畸变有关。晶格对称性会约束能带简并,晶格畸变则可能打开带隙、移动能谷、改变轨道杂化强度。

图7. 高对称路径变化对能带图呈现方式的影响,体现晶格、倒易空间和电子能级之间的联系。DOI:10.1038/s41524-020-00383-7

当晶格与能带放在同一幅物理图像中时,周期势场把原子基元和电子能级连接起来;原子基元在晶格中重复,形成周期性的电子势,电子在这个势场中运动,允许能级变成随 k 变化的能带。实空间周期决定倒易空间结构,倒易空间采样又决定我们怎样观察电子态。声子谱、光学跃迁、费米面和电输运计算,也都建立在这套实空间与倒易空间的互相转换上。

对于非晶、液体、纳米团簇或强无序体系,严格晶格对称性会减弱甚至消失,普通能带路径的物理意义随之下降。研究者有时改用态密度、局域态分析、展宽谱函数或大尺度统计描述。晶格适用于长程周期清晰的体系,一旦研究对象转向强无序或局域缺陷,图谱解释就要补充局域结构信息。

四、在材料计算中使用晶格时要看哪些边界

晶格给了材料计算一个清晰的几何框架,但它并不能替代化学和物理分析。相同或相近的晶格常数下,不同元素会带来不同价电子数、电负性、轨道能级和成键方式;相同化学式下,不同晶型也可能拥有完全不同的配位环境和电子态。晶格相似不等于性质相同,结构判断还应纳入基元、对称性、局域配位和电子占据。

结构松弛后的晶格参数常被用来检查计算模型是否合理,但它只检验平衡几何的一部分。泛函选择、色散修正、磁序、U 值、压力条件和温度效应都会影响晶胞体积与键长。单一晶格指标不能推出全部性质,尤其不能直接替代能带、态密度、声子、弹性常数、吸附能或反应自由能等专门分析。

图8. 材料数据库中的三维晶体结构与计算流程关系,说明晶格信息经过结构松弛后才成为可分析的计算对象。DOI:10.1038/s41597-021-01022-6

在缺陷、掺杂、表面和界面模型中,晶格不再只是原始体相单胞的重复,而是被改写成带有镜像关系的超胞框架。缺陷超胞通过扩大周期单元来降低镜像相互作用,表面 slab 通过真空层隔开上下表面,异质结模型要处理晶格失配、应变分配和界面取向。超胞晶格本身就是物理假设,它决定了周期镜像、缺陷浓度、界面周期和应变状态。

晶格也不是永远刚性的背景。热膨胀、相变、软模、铁电畸变、Jahn-Teller 畸变和电化学界面重构都会改变晶体的平均结构或局域结构。静态 DFT 常从零温近似结构出发,分子动力学、声子自由能或准谐近似才会把温度引入结构讨论。晶格描述的是模型中的周期骨架,研究热运动和动态重构时还应引入时间尺度和统计信息。

在阅读材料计算结果时,较好的习惯是先识别三个层次:晶格基矢决定周期边界,基元和坐标决定局域配位,电子结构计算给出能量、电荷和轨道响应。晶格负责把有限模型连接到无限晶体,它是很多图谱的共同背景,却不是所有结论的唯一来源。把这三个层次放稳,晶胞参数、空间群、k 点路径、能带图和缺陷超胞就会变成同一套语言中的不同部分。