基于原子尺度位移定量测量的极性构型熵研究

一、 作者、机构与发表信息

本研究由宾夕法尼亚州立大学材料科学与工程系的Debangshu Mukherjee、Leixin Miao、Venkatraman Gopalan和Nasem Alem，列日大学（比利时）的Sergei Prokhorenko和Eric Bousquet，以及阿肯色大学的Sergei Prokhorenko合作完成。该研究成果以题为“Atomic-scale measurement of polar entropy”的论文形式，于2019年9月3日发表在物理学领域知名期刊《Physical Review B》（第100卷，第104102期）上。

二、 学术背景与研究目的

本研究属于凝聚态物理与材料科学的交叉领域，具体聚焦于铁电材料的基础物理性质研究。熵是热力学中的一个基本物理量，它量化了系统可能存在的微观状态数量，系统的稳定性由其总熵的大小决定。然而，传统的熵变测量（如量热法）主要是宏观尺度的，直接在原子尺度上定量测量熵的变化极为罕见。铁电材料具有自发且可翻转的电极化，其根源在于晶格畸变导致的离子位移。长期以来，典型的单轴铁电体（如铌酸锂，LiNbO₃）被视为伊辛（Ising）型铁电体，即其极化矢量被限制在晶体学允许的单一轴向（如〈0001〉方向）上，在畴壁处极化方向发生180°反转。然而，近期理论研究表明，在其他铁电体（如钛酸铅）中，极化矢量可能存在偏离理想伊辛方向的涨落（如布洛赫型和奈尔型分量）。这种局域极化序参量的可变性（即“极性无序”）会引入构型熵（configurational entropy），可能对材料性能产生重要影响。尽管理论上有所预测，但在LiNbO₃这类经典“硬性”铁电体中，实验上直接观测并定量这种原子尺度的极性无序及其对应的熵，仍是一个巨大的挑战。

本研究的主要目的是：利用具有亚埃分辨率的球差校正扫描透射电子显微镜（STEM），在原子尺度上直接测量LiNbO₃单晶中铌（Nb）和氧（O）原子柱相对于其顺电位置的位移，从而量化由极性位移变化引起的构型熵（文中称为“极性熵”）。研究旨在探索在体单畴区和畴壁附近，极化矢量是否真的被严格限制在伊辛方向，以及这种极性无序的程度如何，并通过结合第一性原理计算和平均场有效哈密顿量方法，阐明极性熵在稳定局域对称破缺中的作用。

三、 详细研究流程

本研究是一个结合了先进实验表征、精密数据分析和理论计算的系统性工作，主要流程可分为以下几个关键步骤：

1. 样品制备与表征：

   • 研究对象： 采用商业购买的周期性极化铌酸锂（LiNbO₃）单晶，其空间群为R3c。研究关注的是晶体中的180°铁电畴壁以及远离畴壁的体单畴区域。
   • 样品制备： 使用聚焦离子束（FIB）技术制备电子透明的透射电镜样品。为确保原子尺度成像（尤其是轻元素氧）的质量并最小化样品表面损伤，最终抛光步骤采用了低电压（0.5 kV和2 kV）离子束。通过电子能量损失谱（EELS）测量非弹性平均自由程，确定样品厚度约为20-25纳米。
   • 初步表征： 首先使用常规透射电镜（CTEM）在稍离焦条件下成像，利用衍射衬度定位并识别畴壁。选区电子衍射（SAED）确认观察方向为[11̄00]晶带轴，此方向使极化位移（沿〈0001〉方向）位于成像平面内，便于测量。

2. 原子尺度位移测量：

   • 实验设备与方法： 在球差校正的FEI Titan³透射电子显微镜中进行扫描透射电子显微镜（STEM）成像。同时采集环形暗场像（ADF-STEM）和明场像（BF-STEM）。本研究主要利用BF-STEM图像进行位移定量，因为该技术可以同时精确定位重原子（Nb）和轻原子（O）的柱位置，且对样品倾斜和离焦变化相对不敏感。
   • 数据采集： 在畴壁区域和远离畴壁（约100纳米）的体单畴区域分别获取高分辨BF-STEM图像。成像后，使用专门的算法对图像进行扫描漂移校正。
   • 位移计算流程： a. 原子位置精修： 使用自行开发的MATLAB脚本，首先定位图像中原子柱的强度峰值，然后通过拟合二维高斯函数对原子强度分布进行亚像素精度的细化，使位移测量精度达到约2皮米（pm）。 b. 建立平均晶胞： 将整幅图像中所有测量到的原子位置分配到对应的LiNbO₃晶胞位置，然后对所有晶胞进行平均，得到一个“平均晶胞”作为参考基准。 c. 计算位移矢量： 对于每个独立的晶胞，分别计算其Nb原子和O原子的质心相对于平均晶胞中心的位移矢量。 d. 提取极化位移： 将每个晶胞的Nb位移矢量减去O位移矢量，得到该晶胞的Nb-O极性位移矢量。 e. 分解位移分量： 将每个极性位移矢量分解为沿〈0001〉方向的“伊辛（Ising）分量”和沿〈112̄0〉方向的“奈尔（Néel）分量”。前者对应经典理论预期的极化方向，后者则代表垂直于理想极化方向的横向涨落。

3. 第一性原理与有效哈密顿量理论计算：

   • 计算目的： 从理论上理解实验中观察到的极性不稳定性的起源，并预测材料固有的位移能量学。
   • 方法： 采用基于密度泛函理论（DFT）的ABINIT软件包，计算顺电相LiNbO₃的声子谱。识别出不稳定的声子模式，特别是驱动铁电性的A₂u模式（对应伊辛位移）和简并的Eu模式（对应奈尔和布洛赫位移）。
   • 能量景观构建： 通过计算原子沿不同模式位移时的总能量变化，构建了与Eu模式相关的能量景观。结果显示，在抑制了伊辛位移的〈0001〉方向后，系统能量在垂直于该方向的平面内呈现著名的“戈德斯通草帽势（Goldstone sombrero potential）”形状，表明存在稳定的非伊辛位移（奈尔/布洛赫分量），其径向幅度在约8 pm时达到能量最小值。
   • 平均场有效哈密顿量模拟： 基于DFT计算结果，构建了包含模式相互作用和应变耦合的有效哈密顿量模型，并在平均场近似下进行模拟，用于计算位移的概率分布和理论熵值。

4. 极性熵的量化分析：

   • 概率分布计算： 将实验测量得到的所有晶胞的伊辛和奈尔位移幅度（由于STEM是二维投影，无法直接测量沿电子束方向的布洛赫位移）按0.1 pm的精度分箱，统计每个位移配置出现的频率，从而得到位移的联合概率分布函数（PDF）。
   • 熵计算公式： 采用吉布斯-玻尔兹曼公式 S = -k_B Σ ρ log(ρ) 计算极性构型熵。其中，k_B 是玻尔兹曼常数，ρ 是观测到特定位移配置的概率，求和遍及所有可能的位移状态。
   • 仪器效应去卷积： 认识到实验测量的位移概率分布包含了材料本征无序、电子显微镜不稳定性和电子束-样品相互作用引入的额外涨落。为了估算并扣除非本征贡献，研究团队假设仪器效应可用一个二维高斯函数（点扩散函数，PSF）来描述。他们采用理查森-露西（Richardson-Lucy）迭代反卷积算法，从实测概率分布中反演出本征的概率分布。通过寻找使反卷积后熵值降低最小的最佳高斯函数宽度（σ），来保守地估计仪器效应的贡献（体畴区σ约为3.⁶⁄₈.7 pm，畴壁附近σ约为11.7 pm）。最终用于计算熵的是反卷积后的概率分布。

四、 主要研究结果

1. 畴壁与体畴区的极性位移特征：

   • 畴壁结构： 实验直接观测到LiNbO₃中180°畴壁的原子结构。畴壁处伊辛位移方向在1-2个晶胞内发生反转，Nb和O原子协同运动，产生的净伊辛位移约为55 pm。畴壁并非完全平直，存在弯曲和扭折。
   • 非伊辛位移的普遍存在： 最关键的发现是，不仅在畴壁处，甚至在远离畴壁的体单畴区域，都观测到了显著的奈尔分量位移（即极化矢量偏离〈0001〉方向）。这与经典伊辛模型的预期相悖。
   • 位移幅度的空间变化： 定量分析表明，奈尔位移的绝对值在畴壁附近显著增加（可达约10 pm），并伴有更大的离散度；而在体畴区（距离畴壁约100 nm），奈尔位移减小至5 pm以下，但并未完全消失。伊辛位移在体畴区内基本保持恒定，但在畴壁处也显示出更大的变异性。

2. 理论计算对位移能量学的解释：

   • DFT计算确认了顺电相LiNbO₃中存在不稳定的Eu声子模式，其对应的原子位移方向垂直于铁电极化主轴（〈0001〉方向）。这意味着，即使在体畴区，由于Eu模式的不稳定性，局域偶极子方向偏离理想的伊辛轴在能量上也是允许的，从而为实验中观测到的非伊辛涨落提供了理论依据。
   • 计算得到的Eu模式能量景观表明，在抑制伊辛分量后，系统在垂直于〈0001〉的平面内具有“草帽”形的势能面，其能量最小值对应的径向位移约为8 pm，与实验在畴壁处观测到的奈尔位移幅度（~10 pm）接近。

3. 极性熵的定量结果：

   • 体畴区的熵： 基于反卷积后的位移概率分布，计算得到体畴区的极性熵。实验测得伊辛分量贡献的熵约为0.40 meV/K，奈尔分量贡献的熵约为0.31 meV/K。理论有效哈密顿量模拟预测的对应值分别为0.44 meV/K和0.40 meV/K，两者吻合度在10%以内，强烈支持了所观测到的极性无序是材料本征性质。
   • 畴壁附近的熵： 在畴壁附近（±10 nm范围内），实验测得的极性熵比体畴区高出约28%。熵的增加源于畴壁处位移概率分布更加弥散（双峰分布，源于两个畴的叠加）和位移幅度的更大涨落。
   • 熵的物理意义： 这些非零的熵值明确证实，即使在看似均匀的LiNbO₃单晶单畴区内，也存在着可观的、由局域极化矢量涨落引起的构型熵。畴壁作为缺陷，进一步增加了系统的无序度，从而具有更高的极性熵。

五、 研究结论与意义

本研究的核心结论是：在经典的、被认为是理想伊辛型的铁电体LiNbO₃中，存在本征的、原子尺度的极性无序（非伊辛涨落），并首次在实验上直接测量并量化了由此产生的极性构型熵。 畴壁作为拓扑缺陷，会进一步增强这种无序，导致局部熵的增加。

这项工作的科学价值和应用意义在于： * 基础科学层面： 挑战了对传统单轴铁电体极化行为的简化认知（纯伊辛模型），揭示了即使在高温下表现出高稳定性的铁电体中，有限温度下也存在显著的局域对称破缺和序参量涨落。它将熵这一核心热力学概念的研究从宏观尺度推进到了原子尺度，建立了通过原子位移计量学直接量化构型熵的新范式。 * 方法论层面： 发展了一套结合亚埃分辨率STEM位移测量、概率统计分析、仪器效应反卷积和第一性原理计算的研究方法，为在更广泛的材料体系（不限于铁电体）中研究原子尺度的无序和熵提供了强大的工具。使电子显微镜不仅能用于成像，还能成为进行原子分辨率热力学量化的仪器。 * 材料性能关联： 极性无序（熵）是影响铁电、压电、电热等性能的关键因素。本研究揭示的这种本征无序可能存在于更多被认为“均匀”的铁电材料中，为理解和设计具有优异功能特性的材料（如高性能压电、电热材料）提供了新的微观视角和调控思路。

六、 研究亮点

1. 开创性的直接测量： 首次在原子尺度上实验测量并定量了铁电材料中的极性构型熵，实现了从宏观热力学测量到微观状态统计的跨越。
2. 颠覆性的实验发现： 在公认的“硬性”伊辛铁电体LiNbO₃的体单畴区内，明确观测到非伊辛型极化涨落（奈尔分量），修正了对其极化行为的传统理解。
3. 严谨的多方法交叉验证： 研究将顶尖的STEM实验技术、精密的位移计量学、严谨的概率统计分析与第一性原理计算紧密结合，实验与理论结果相互印证，结论坚实可靠。
4. 创新的数据分析流程： 开发了完整的从原子位置提取到熵计算的数据处理流程，并创新性地采用反卷积方法处理仪器效应，力求提取材料的本征熵信息。
5. 方法的普适性： 所建立的方法学框架具有通用性，可推广应用于研究其他涉及原子位移、有序-无序相变、局域对称性破缺的凝聚态体系。

七、 其他有价值的内容

研究在附录中提供了大量补充信息，进一步支撑了主要结论的可靠性，包括： * 样品厚度与制备细节： 详细说明了FIB制备参数和EELS测厚结果，论证了样品质量足以进行氧原子尺度的精确测量。 * 多个畴壁与体畴区域的补充数据： 展示了除正文主要区域外的另外四个畴壁和两个体畴区域的位移图，证明了所观察到的非伊辛位移和熵增现象在样品中具有普遍性和可重复性。 * 静电势能估算： 基于玻恩有效电荷和测量的位移，估算了由极性涨落引起的电荷积累和静电势能。结果显示，在300 K时，熵项（-TΔS）对自由能的贡献（约-200 meV量级）远大于静电排斥能的增加（<0.5 meV），从能量上解释了为何静电相互作用不会抑制这种极性涨落。 * 模拟验证： 通过多片层模拟验证了样品倾斜对测量结果的影响很小，排除了实验假象的可能性。

这项研究是一项里程碑式的工作，它通过精妙的实验设计和深入的理论分析，在原子尺度上揭示了铁电材料中隐藏的复杂性，并将“熵”这一基本概念与实空间的原子排列直接联系起来，为理解和调控复杂功能材料提供了全新的视角和强有力的工具。