粗粒化原子密度场与配位数在晶体缺陷描述中的关系与区别

《RSC Advances》：On the relationship and distinction between atomic density and coordination number in describing crystal defectsOpen Access

　　

晶体缺陷常通过断键模型和最近邻配位数来描述。本研究强调，晶体缺陷导致的原子间距内在异质性（这在断键模型中常被忽略）对缺陷性质有决定性影响。研究人员分析了大量由分子静力学弛豫的BCC-Fe晶界的过剩自由体积、能量和熵，并证明作为一个系统性粗粒化场变量，原子密度场与晶界性质的相关性远优于配位数。研究表明，原子密度场能够同时捕捉键缺失和间距变化，从而涵盖并统一了结构和体积信息。这些结果确立了粗粒化原子密度场作为描述晶界热力学的原则性替代变量，激励在界面能量学、过剩性质、偏析和相行为的预测理论中对断键规则进行补充，并为从原子尺度到介观尺度的模型提供了直接桥梁，其潜在应用可扩展至其他微结构缺陷。

一、研究背景与科学问题

晶体缺陷的传统描述依赖于断键模型（broken-bond model），该模型将界面视为二维晶体学缺陷，其中原子配位环境发生重排，并通过配位数（coordination number, CN）的变化来衡量界面形成的能量代价。配位数作为局部结构描述符，广泛应用于晶界能量、溶质偏析（segregation）等理论模型中。然而，配位数方法本质上存在局限：它通过离散的邻居计数来分类原子环境，却无法捕捉原子间距（interatomic spacing）、局域堆积（local packing）或体积畸变（volumetric distortion）的变化。这一局限性促使研究者探索粗粒化场描述符，如原子密度（atomic density），它既能保留邻居多重性的信息，又能反映空间排列特征。为此，本研究旨在揭示配位数与粗粒化原子密度在描述晶体缺陷时的关系与区别。

二、研究方法概述

研究人员基于一套包含408个不同BCC-Fe晶界（含倾斜、扭转和混合类型）的数据集，这些晶界通过分子静力学（molecular statics）生成并弛豫。研究采用粗粒化原子密度场（coarse-grained atomic density field）作为结构描述符，将其定义为三维高斯分布函数卷积原子位置得到的连续场变量，并与传统的第一近邻配位数进行对比。通过计算晶界平面处的平均原子密度ρ^GB和平均配位数Z₁^GB，研究人员系统评估了这两个描述符与晶界过剩自由体积（excess free volume per unit area ΔV）、晶界能量（grain boundary energy γ^GB）以及两种构型熵（configurational entropy）指标——基于原子密度场的S_ρ^GB和基于壳层描述符的和熵（sum entropy）之间的相关性。

三、主要研究结果

代表性晶界示例（Fig. 1）： 对于Σ3、Σ5和Σ9晶界，原子密度场在界面处呈现出清晰的耗尽区，从而稳健地定位晶界区域；而配位数Z₁仅记录最近邻连接性发生变化的原子环境。尽管两个描述符均指示了界面的大致位置，但原子密度场为所有示例提供了一致的晶界区域定义，而配位数的变化则更具个案特异性。

密度与配位的相关性（Fig. 2）： 密度亏损（1 - ρ^GB）与配位亏损（8 - Z₁^GB）/8之间存在松散的整体相关性，表明配位损失确实贡献于原子密度亏损，但数据点散布显著。相似的密度亏损可能对应截然不同的配位亏损，证实ρ^GB与Z₁^GB是相关但不等价的结构描述符。

过剩自由体积与密度及配位的关系（Fig. 3）： 过剩自由体积ΔV与密度亏损（1 - ρ^GB）呈近乎完美的线性关系，单一晶界平面量ρ^GB即可可靠预测积分过剩体积。相反，ΔV与配位亏损之间的关系非常分散，表明配位损失单独不足以决定界面的体积状态。

晶界能量与密度及配位的关系（Fig. 4）： 晶界能量γ^GB与密度亏损和配位亏损均呈现总体趋势，但基于密度的表征显示出更清晰、更明显的家族依赖趋势，而基于配位的表征则相当松散和高散。

构型熵与密度及配位的关系（Fig. 5）： 和熵随密度亏损和配位亏损的增加整体上升，反映了随密度耗尽而增长的拓扑距离相关无序。值得注意的是，配位亏损引起的熵贡献远小于原子间距扰动引起的熵贡献（约5至10倍），证实当间距波动占主导时，配位数本身会低估构型无序。

密度熵与和熵的等价性（Fig. 6）： 基于原子密度场的熵S_ρ^GB与基于壳层描述符的和熵沿y=x线紧密聚集，斜率和截距均接近理想值，验证了密度熵与和熵的数值等价性，从而确认原子密度场紧凑地捕获了由壳层描述符编码的联合拓扑-距离相关无序。

四、讨论与结论

本研究的核心发现是配位损失与密度耗尽之间存在松散而非唯一的关系。这符合物理预期：最近邻多重性的改变直接影响局域原子密度，而原子密度场额外包含了原子间距变化的信息，后者在配位数描述中完全缺失。通过比较两种描述符与独立确定的晶界过剩性质，研究明确证明了原子密度场比配位数更能忠实地反映晶界结构。具体而言，原子密度场与过剩自由体积呈近乎完美的线性关系，在晶界能量表征中展现出更清晰的家族依赖趋势，并且其熵度量与包含多壳层信息的和熵高度一致。这些结果确立了粗粒化原子密度场作为描述晶界热力学的原则性替代变量，激励在界面能量学、过剩性质、偏析和相行为的预测理论中补充断键规则，并为从原子尺度到介观尺度的模型提供了直接桥梁，其潜在应用可扩展至其他微结构缺陷。