科学报告：钙钛矿氧化物中铁电性的起源
主要作者及发表信息
本项开创性研究由华盛顿卡内基研究所（Carnegie Institution of Washington）地球物理实验室（Geophysical Laboratory）的Ronald E. Cohen独立完成。该研究成果以题为“Origin of ferroelectricity in perovskite oxides”（钙钛矿氧化物中铁电性的起源）的论文，于1992年7月9日发表在《Nature》杂志第358卷第136至138页。这一研究在当时为理解铁电材料的微观机制提供了全新的视角。
研究背景与目的
铁电材料（ferroelectric materials）是一类具有可翻转宏观极化特性的重要功能材料，在电容器、存储器、传感器等电子器件领域有着广泛的应用前景。技术应用上最重要的铁电体大多是具有钙钛矿结构（perovskite structure）的氧化物。然而，长期以来，这些材料中铁电行为的起源一直不明确，一个关键的未解之谜是：为什么化学组成相似但不同的钙钛矿材料会表现出截然不同的铁电行为——例如，钛酸钡（BaTiO₃）和钛酸铅（PbTiO₃）虽然具有相似的晶胞体积等凝聚特性，其相变行为却大相径庭：BaTiO₃在降温过程中依次经历立方相到四方相（393 K）、四方相到正交相（278 K）、正交相到菱方相（183 K）的三次铁电相变，而PbTiO₃仅在766 K发生一次立方相到四方相的相变，且其四方相的铁电畸变更大，c/a应变高达6%（BaTiO₃仅约1%）。此外，学术界对于这些相变应被描述为位移型（displacive）还是有序-无序型（order-disorder）也存在争议。早先的研究提出了“八位点”模型（eight-site model），认为高温立方相中原子沿八个[111]方向随机畸变，四方相则沿四个对角方向有序排列。然而，该模型是否对两种材料均适用尚未得到验证。
本研究的目的正是通过高精度的第一性原理计算，揭示BaTiO₃和PbTiO₃中铁电性的微观起源，并解释二者铁电相变行为差异的本质原因。研究假设共价杂化（covalency）和电子结构的差异可能是决定铁电行为的关键因素。
研究方法与技术流程
为精确描述铁电相变中长程库仑力（long-range Coulomb forces，有利于铁电态）和短程排斥力（short-range repulsions，有利于非极性立方结构）之间的微妙平衡，Cohen采用了全电子、全势、线性化缀加平面波方法（all-electron, full-potential, linearized augmented plane-wave method，LAPW方法），并基于局域密度近似（local density approximation, LDA）处理电子间的交换关联作用。这一计算方案不做任何关于键合、离子性、电荷密度或势函数形状的人为假设，能够精确处理电荷畸变和共价效应，并达到了高度收敛。该方法的计算量极为庞大，在当时的高性能计算条件下是一项挑战性工作，研究使用了国家超级计算应用中心的Cray 2超级计算机完成。
研究以实验观测的畸变结构为指引，重点考察了菱方（rhombohedral）和四方（tetragonal）两种布里渊区中心（zone-centre）畸变模式。具体计算流程包括：
第一步，以实验结构为初始构型，计算BaTiO₃和PbTiO₃在有无四方c/a应变条件下，能量随软模（soft mode）畸变幅度变化的势能面。计算涵盖压缩体积和实验体积两种条件。
第二步，针对PbTiO₃的实验四方相结构，采用“冻结声子”（frozen phonon）方法计算动力学矩阵（dynamical matrix）。具体操作是对各个原子施加微小位移，计算由此产生的力，进而构建力常数矩阵，再通过对角化求得谐波近似下的本征模频率（eigenmode frequencies），并与实验拉曼光谱数据进行系统比对。
第三步，计算并分析两种材料的电子态密度（density of states, DOS），包括总态密度和投影到各原子轨道的分态密度（partial density of states）。为将Ba和Pb的差异纯粹地分离出来，Cohen还使用PbTiO₃的实验结构参数，仅将Pb原子替换为Ba原子，重新进行了态密度计算。
第四步，对PbTiO₃在不同铁电畸变幅度下的价电荷密度（valence charge density）进行三维可视化分析，以直观展示共价键合和离子极化的空间分布特征。
第五步，进行两项控制性验证计算：一是采用假定离子为球形且无共价性的“势诱导呼吸”（potential induced breathing）模型进行计算；二是在LAPW计算中将Ti的3d轨道能量参数人为提高以限制其变分自由度，观察铁电不稳定性是否消失。这两项验证旨在确认Ti-O杂化的必要性。
主要研究结果
能量-畸变关系的计算结果（论文图1）揭示了若干关键发现：
在不施加应变时，无论是BaTiO₃还是PbTiO₃，能量最低的结构均为菱方相。但两者的势阱深度存在巨大差异：PbTiO₃的势阱远深于BaTiO₃。当加入四方c/a应变后，情况发生根本性改变。在PbTiO₃中，四方应变使四方相的能量显著降低，从而稳定了四方相，使其相较菱方相更为稳定——带应变的实验四方结构比立方结构能量低逾1,800 K。而在BaTiO₃中，四方应变对能量的影响要小得多，菱方相仍为能量最低相，该能量差仅约150 K。这一结果完美解释了PbTiO₃具有单一四方铁电相而BaTiO₃最终趋向菱方相的实验事实。此外，两种材料均表现出强烈的体积效应：在压缩条件下，势阱深度减小，软模势能面的多重阱结构逐渐消失，这与实验上钛酸盐在高压下失去铁电性的现象一致。
声子计算结果显示，对于PbTiO₃计算得到的A1和E模频率与实验拉曼数据吻合良好，验证了计算方法的可靠性。最低频率模的稳定存在表明PbTiO₃在四方相不存在进一步相变的不稳定性；而BaTiO₃在四方相则仍保持动力学不稳定，这与它继续发生相变的实验事实相符。
态密度分析（论文图2）揭示了铁电性的电子起源：Ti 3d态与O 2p态之间存在强烈的杂化（hybridization），且这种杂化在铁电畸变下进一步增强。铁电畸变导致能带特征值之和降低，从而降低了体系总能量。更具启示性的是PbTiO₃和BaTiO₃的对比分析：在PbTiO₃中，Pb 6s态与O 2p态也存在显著的杂化，而Ba的5p态则与价带完全无杂化，呈现纯粹的离子性特征。由于Pb-O的共价键合作用以及Pb²⁺的离子半径小于Ba²⁺，PbTiO₃产生了更大的四方应变，进而缩短了部分Pb-O距离。同时，Ti-O之间的排斥作用阻止了体积过度收缩以稳定立方相。Pb-O杂化还对Ti-O相互作用产生了间接影响——即使采用相同的原子位移，PbTiO₃中Ti 3d轨道的本征值也低于BaTiO₃。
电荷密度分布图（论文图3）直观展示了PbTiO₃中随着铁电畸变增加，Pb离子周围出现的偶极电子密度分布，而BaTiO₃中的Ba则不发生此类极化。对自洽电荷密度的拟合给出了O的电荷为-1.63、Ti的电荷为+2.89，表明Ti并非完全电离的Ti⁴⁺离子，这与传统离子模型的观点形成鲜明对比。
两项验证性计算均证实，一旦Ti-O杂化被抑制，铁电不稳定性即告消失，立方相变得最为稳定。这确凿地证明了Ti 3d态与O 2p态的杂化是产生铁电性的必要条件。
研究结论与意义
基于上述结果，本研究提出了关于钙钛矿氧化物铁电性起源的统一性理论框架：
第一，铁电性得以发生的根本机制在于B位阳离子（b cation，如Ti⁴⁺, Nb⁵⁺, Zr⁴⁺等）与氧离子的d轨道杂化。这种共价键合削弱了短程B-O排斥力，使得长程库仑力能够驱动铁电畸变。大多数铁电氧化物钙钛矿的B位阳离子均具有空的d态，这为d-p杂化提供了前提条件。
第二，A位阳离子（a cation）通过两种途径调控基态结构和相变性质：（1）增大铁弹应变（ferroelastic strain），使其与铁电畸变耦合；（2）与氧形成共价杂化，间接改变B-O相互作用。
第三，短程力（有利于立方相）和长程力（有利于铁电态）之间的微妙平衡，使铁电相变对缺陷（改变短程相互作用）和载流子（屏蔽长程电场）高度敏感。
本研究不仅解开了长期困扰学界的BaTiO₃与PbTiO₃相变行为差异之谜，还建立了一个具有普适性的理论框架，对理解其他钙钛矿氧化物乃至相关材料（如铜氧化物超导体）的晶格不稳定性具有重要的借鉴意义。这一工作代表了将第一性原理计算方法成功应用于复杂氧化物功能材料研究的早期典范，有力地推动了从“唯象描述”向“微观机理理解”的转变。
研究亮点
本研究的突出亮点体现在以下几个方面：第一，其核心发现在于明确揭示了Ti 3d与O 2p之间的共价杂化是钙钛矿铁电性的必要条件，彻底修正了此前以离子模型为基础的理解范式。第二，方法学上，研究采用了全电子全势LAPW方法这一在二十世纪九十年代初期极为先进的电子结构计算手段，其计算精度和对共价效应的描述能力远超同时代常规方法。第三，研究设计了巧妙的对比分析策略——通过将BaTiO₃置于PbTiO₃的结构中进行计算，实现了对A位阳离子化学效应的隔离分析。第四，研究将宏观相变行为（相变温度、相序、应变）与微观电子结构（态密度、轨道杂化、电荷密度）建立了定量的因果联系，提供了完整的物理图像。第五，论文最后明确指出，具有底层非谐波势能面的高温相——其中原子几乎不处于理想平衡位置——并非铁电体独有的特征，也是钙钛矿旋转不稳定性和铜氧化物超导体倾斜不稳定性的典型表现，这一洞察将该研究的理论价值扩展至更广阔的凝聚态物理领域。