这篇研究发表于《自然·材料》（*Nature Materials*）期刊的2025年9月卷，由来自美国康奈尔大学、斯坦福大学和SLAC国家加速器实验室、阿肯色大学等多个机构的Harikrishnan K.P.、Ruijuan Xu、David A. Muller等学者共同完成。该工作报道了一项关于无铅铁电材料铌酸钠（NaNbO₃）薄膜的突破性研究。
学术背景 本研究的核心科学领域是凝聚态物理和材料科学，具体聚焦于复杂氧化物材料的铁电性（ferroelectricity）微观机制。铌酸钠是一种备受关注的无铅铁电/反铁电候选材料，具有丰富的相图和多个相互竞争、微妙平衡的基态。在体材料中，铌酸钠在室温下通常呈现反铁电态，但通过外延应变（epitaxial strain）可以稳定其宏观铁电态。然而，这种应变稳定铁电相的微观结构一直不明确，特别是决定其极化（polarization）的核心原子位移细节。
理解这一问题至关重要，因为传统上对钙钛矿铁电体的研究常以阳离子亚晶格（如B位阳离子相对于A位阳离子的位移）的偏移作为极化的替代指标。这种方法对于阳离子位移明显的典型位移型铁电体（displacive ferroelectrics）是有效的。但该假设并不具有普适性，可能忽视了由阴离子位移主导的机制。由于氧原子等轻元素散射能力弱，在电子显微镜等表征技术中难以被精确观测，这使得准确解析阴离子的贡献充满挑战。因此，本研究旨在利用最先进的显微成像技术，直接、精准地揭示应变稳定铌酸钠铁电相的原子级三维结构，特别是氧原子的位置，以阐明其铁电性的微观起源。
详细研究流程 本研究是一项综合性实验与理论研究，主要流程包括样品制备、多尺度结构表征、先进显微成像、理论计算和数据分析。
1. 样品制备与初步表征： 研究团队首先使用脉冲激光沉积（pulsed laser deposition）方法，在（110）取向的DyScO₃衬底上外延生长了58纳米厚的铌酸钠薄膜。DyScO₃衬底提供了双轴拉伸应变。通过X射线衍射（X-ray diffraction, XRD）和三维倒易空间映射（3D reciprocal-space mapping）证实了薄膜具有高质量的结晶性和与衬底相干的拉伸应变。利用压电力显微镜（piezoresponse force microscopy, PFM）对薄膜进行了宏观铁电畴表征，观察到了清晰的面内极化畴结构，平均畴尺寸约为50纳米，这证实了薄膜在室温下确实具有宏观铁电性。
2. 电子显微学表征与多切片电子叠层衍射成像技术（Multislice Electron Ptychography, MEP）： 这是本研究的核心技术环节。为了在原子尺度解析所有原子（特别是轻的氧和钠原子）的位置，研究采用了扫描透射电子显微镜（scanning transmission electron microscopy, STEM）下的多切片电子叠层衍射成像技术。这是一种革命性的成像方法。
技术原理与创新性： 传统STEM成像模式，如环形暗场像（annular dark field, ADF）、环形明场像（annular bright field, ABF）或积分微分相位衬度像（integrated differential phase contrast, iDPC），在观测轻原子或存在多重散射（multiple scattering）的较厚样品时存在局限。ADF对重原子（如铌）敏感但对轻原子几乎无信号；ABF和iDPC虽能显示轻原子，但分辨率有限，且对样品厚度、衍射条件、像差和表面结构极为敏感，易产生系统误差和假象，难以进行皮米级精度的定量分析。 多切片电子叠层衍射成像技术（MEP）从根本上克服了这些限制。其实验设置如文中图1a所示：在STEM模式下，用会聚的电子束在样品上逐点扫描，在每个扫描点，使用像素化探测器（本研究中采用EMPAD探测器）记录完整的二维衍射图案，从而获得一个四维数据集（两个实空间扫描坐标，两个倒空间衍射坐标）。通过迭代算法求解电子散射的逆问题，并利用多切片法（multislice method）在正向模型中考虑电子束传播过程中的多重散射和通道效应，最终能够重建出样品电势的三维分布。该方法不仅突破了衍射极限，实现了亚埃级的横向空间分辨率，还能沿电子束方向进行“深度切片”（depth sectioning），恢复三维结构信息。
样品处理与数据采集： 用于STEM分析的样品通过聚焦离子束（focused ion beam）技术制备成横截面薄片。MEP数据在300 kV的球差校正STEM上采集，使用30 mrad的半会聚角，步长为0.42 Å。每个扫描点的衍射图样被记录下来，用于后续的重建。
数据处理与重建： 相位恢复和三维电势重建使用最大似然多切片叠层衍射算法（maximum-likelihood multislice ptychography algorithm），在fold-slice软件包中实现。该算法采用基于梯度的方法进行位置校正，并使用多个探针模式（probe modes）来考虑电子束的部分相干性。研究人员还利用贝叶斯优化模型来估计重建的最佳参数（如会聚角、离焦量、样品厚度）。最终的重建结果以一系列沿电子束方向的“切片”形式呈现，每个切片对应样品不同深度的电势投影。
3. 原子位置测量与极化计算： 从重建的三维电势图中，使用二维高斯拟合（2D Gaussian fitting）精确定位所有原子柱（包括Na、Nb、O）的位置。对于每个以铌原子为中心的赝立方晶胞，计算周围四个氧原子的质心以及四个钠原子的质心。微观偶极矩由两个矢量决定：1）氧原子质心到铌原子的矢量；2）氧原子质心到钠原子质心的矢量。通过将形式电荷或基于密度泛函理论（DFT）计算得到的玻恩有效电荷（Born effective charges）分配给各原子，可以估算每个晶胞的偶极矩，进而计算出极化强度。
4. 密度泛函理论（DFT）计算与声子谱分析： 为了支持实验观察并从理论上理解其机制，研究团队进行了系统的第一性原理计算。他们模拟了在不同外延应变下铌酸钠薄膜的结构，搜索与实验结构相匹配的相，并通过分解晶体结构来分析其对应的不可约表示（irreducible representations）。此外，还计算了立方相铌酸钠的声子色散关系（phonon dispersion），以理解其结构不稳定性。
5. 界面结构与八面体旋转（octahedral rotation）耦合研究： 利用MEP的深度切片和高分辨率能力，研究团队直接成像了NaNbO₃薄膜与DyScO₃衬底之间的界面原子结构。通过分析界面附近氧原子柱的椭圆度（ellipticity）和八面体旋转角的变化，绘制了衬底八面体旋转模式向薄膜中“模板化”（templating）的过渡行为，确定了其耦合长度尺度。
主要研究结果 1. MEP技术优势的直接证明： 研究首先直观比较了MEP与常规STEM成像模式（ADF, ABF, iDPC）。结果显示（图1b-e），只有MEP能够以无与伦比的分辨率和清晰度同时分辨出Na、Nb、O所有原子柱，甚至可以直接解读Nb-O键长和键角的细微变化。而ABF像显示被抑制的八面体旋转幅度，iDPC像中的氧原子柱则出现沿旋转方向的模糊拉长，这些均是由于对样品表面松弛结构敏感而产生的假象。
2. 三维深度切片揭示表面松弛与纳米尺度相分离： MEP的深度切片能力使研究者能够将样品的铁电体内部结构与表面松弛结构分离开来。如图2所示，在样品的入口和出口表面（即薄片的两个新表面），氧八面体旋转（OOR）的幅度被抑制，原子位置发生弛豫。如果只关注对表面敏感的传统STEM图像，就会错误地表征内部体结构。更重要的是，MEP在样品内部深度方向发现了纳米尺度区域，其中氧八面体旋转的方向发生反转，而极化方向保持不变（图2d-f）。这直接证明氧八面体旋转与极化这两个序参量（order parameters）并非线性耦合。
3. 极化源自阴离子（氧）位移的核心发现： 这是本研究最关键的成果。通过精确测量所有原子位置并排除表面效应，研究人员发现（图3）： * 铌酸钠薄膜中存在大的氧八面体旋转（约7.3°）以及钠、铌阳离子的反极性位移（anti-polar displacements）。这两种模式传统上被认为会抑制铁电极化。 * 然而，拉伸应变诱导了面内极化。这种极化并非源于阳离子亚晶格的大幅偏移。相反，钠和铌阳离子亚晶格几乎保持完美的中心对称性，阳离子之间的净位移极小（约2-6 pm），低于大多数原子尺度散射技术的检测阈值，且与更大尺度的、保持中心对称的反极性位移卷积在一起。 * 极化主要来源于氧原子相对于铌亚晶格的、近乎刚性的面内整体偏移。 如图3c-d所示，如果将纯旋转的氧八面体覆盖在实验图像上，会发现氧原子的实测位置存在系统性偏移。正是这种阴离子（氧）相对于阳离子（铌）质心的位移，打破了反演对称性，产生了偶极矩。极化矢量图（图3f）显示主要为面内极化，与PFM结果一致。而仅基于阳离子（Nb相对于Na质心）位移绘制的矢量图（图3g）则比真实极化小约五倍且有序性差，符号也相反。这彻底颠覆了传统电子显微学中“极化与阳离子间位移成正比”的假设。
4. 单斜相结构与随机八面体旋转模式的直接可视化： 通过MEP成像，研究者确定所研究的铌酸钠相具有面内同相（in-phase）的八面体旋转模式（a⁺b⁺c⁻或a⁺b⁺c⁺，对应单斜Pc或Cm相）。DFT计算表明，在DyScO₃对应的应变下，多个能量极其接近的相（能量差<14 meV/化学式）可能共存，解释了实验观察到的结构复杂性。此外，研究在另一样品区域直接观测到了一种沿面内方向为反相位（anti-phase）旋转的结构。沿[110]方向观察时，氧原子呈现“之”字形排列，并且沿生长方向，相邻晶胞间的八面体倾斜模式（同相或反相）呈现随机序列（图4g），即随机八面体旋转图案。DFT计算的声子谱显示，立方相在布里渊区M点（同相倾斜）和R点（反相倾斜）之间存在一条异常平缓的声子支（图4h），表明两种倾斜模式的能量成本相近，从而导致了实验观测到的随机性。这是首次在实空间直接可视化这种随机倾斜模式及其空间范围。
5. 界面八面体旋转耦合的长度尺度： 通过MEP成像NaNbO₃/DyScO₃界面（图5），研究人员绘制了衬底反相位倾斜模式向薄膜同相位倾斜模式过渡的过程。量化分析表明，氧原子柱椭圆度和八面体旋转角在距离界面约10个晶胞（~4 nm）的尺度上达到渐近值。这为“八面体旋转工程”（octahedral tilt engineering）设定了厚度上限：当薄膜厚度超过约4 nm时，其八面体将恢复“体材料”构型，衬底的模板效应失效。
结论与意义 本研究得出核心结论：在拉伸应变稳定的铌酸钠铁电相中，宏观铁电性源于一种非常规的微观结构，其极化主要由阴离子（氧）相对于阳离子（铌）亚晶格的位移所主导，而阳离子亚晶格本身仅呈现微小的反极性位移并近乎保持中心对称。这一发现揭示了铁电性的一种新机制。
其科学价值在于： 1. 挑战传统认知： 打破了通过阳离子位移推断极化的传统范式，证明了一类铁电体（其极化由阴离子位移主导，且与八面体旋转非线性耦合）的存在，拓宽了对铁电性微观起源的理解。 2. 方法学突破： 充分展示了多切片电子叠层衍射成像技术（MEP）在解析复杂氧化物材料，特别是精确测定轻原子三维位置方面的强大能力。该方法能够克服多重散射，分离表面与体效应，实现皮米精度的定量分析，是未来研究材料结构-性能关系的强有力工具。 3. 对材料设计的启示： 提示在寻找新型（尤其是无铅）铁电材料时，不应局限于阳离子位移明显的体系，那些阳离子亚晶格看似对称但可能隐藏着阴离子主导极化的材料值得重新审视。研究确定的界面耦合长度尺度也为通过超薄层设计调控材料性能提供了具体指导。
研究亮点 1. 革命性的发现： 首次在原子尺度直接实验证实了铁电性可以由阴离子位移主导，而阳离子亚晶格保持准中心对称。 2. 顶尖的技术应用： 将最先进的多切片电子叠层衍射成像技术（MEP）应用于复杂铁电氧化物体系，充分发挥其三维、高精度、对轻元素敏感的优势，解决了长期存在的表征难题。 3. 全面的机制阐释： 结合高分辨率成像、深度切片分析、第一性原理计算和声子谱分析，不仅揭示了主导机制，还解释了随机相共存、序参量解耦等复杂现象背后的物理根源（如平坦声子支）。 4. 重要的附加成果： 直接可视化了氧化物中罕见的随机八面体旋转图案，并定量给出了界面结构耦合的关键长度尺度，兼具基础科学意义和器件应用参考价值。