PbTe的局域极性关联不靠非谐性：机器学习势模拟显示，去掉高阶力常数后关联反而保留
核心概要
立方PbTe在平均结构上保持中心对称，却存在有限尺度的局域极性关联，过去常把这种关联归因于晶格非谐性。用基于密度泛函数据训练的机器学习原子间势做十万原子级分子动力学，再把势映射成可解释的原子间力常数模型，可以逐项开关高阶项。结果显示，去掉三阶和四阶非谐力常数后纳米尺度极性关联依然存在，而把二阶谐波相互作用截断到第三近邻，关联完全消失；低温下关联范围超过50 Å。
深度剖析
PbTe中有限尺度的局域极性关联由长程谐波相互作用支撑，而不是由通常与极性纳米区形成联系在一起的高阶晶格非谐性驱动。 此前的解释把Pb–Te键长分布的非高斯展宽和光学声子峰劈裂当作晶格非谐性的直接后果，并据此推断局域偶极的关联来源。 在200 K下比较包含与不包含三阶、四阶力常数的力常数分子动力学轨迹，排除非谐项后纳米尺度极性关联图样仍然存在，且10 K与200 K之间没有定性差别。
实现方式是把机器学习原子间势的轨迹翻译成按近邻壳层截断的原子间力常数模型，从而能单独保留或删除某一阶、某一距离的相互作用。 机器学习势能给出大尺度、高精度的分子动力学，但本身不直接说明是哪些相互作用造成了观察到的关联。 只有把二阶谐波相互作用保留到第14近邻（无位移对称构型下12.9 Å）时，力常数轨迹才重现机器学习势轨迹中的局域关联行为；截断到第8近邻（9.7 Å）时关联范围明显缩小，截断到第3近邻（5.6 Å）时各方向关联消失。
在倒空间中，这段长程谐波作用表现为横光学声子支在Γ点附近的高曲率，即Γ点频率低、随波矢迅速硬化。 仅凭Γ点横光学模的频率高低不足以判断关联范围，因为把短程力常数整体缩放也能把该频率调到同一数值。 把第3近邻以内的力常数整体乘以0.05，可复现14近邻情形下0.5 THz的横光学模频率，但得到的轨迹中局域极性没有可分辨的空间关联。
启示与展望
结论来自对立方PbTe的机器学习势分子动力学与力常数分子动力学模拟，温度覆盖10 K至300 K，并与同结构NaCl的密度泛函力常数和电势分布作比较；作者把所捕获的关联范围视为下界。
作者指出，当前机器学习势能的最大相互作用范围受训练晶胞尺寸和截断半径等因素限制，即使消息传递层数足够也是如此；他们希望后续实验与理论研究能分辨低温局域对称破缺及相关材料响应。