本研究由 Giulia Venditti（意大利SPIN-CNR研究所）、Maria Eleonora Temperini（罗马大学物理系）、Paolo Barone（SPIN-CNR研究所）、Jose Lorenzana（ISC-CNR研究所）以及 Maria N. Gastiasoro（ISC-CNR研究所/西班牙DIPC）等学者合作完成。该成果于2023年1月25日发表在 Journal of Physics: Materials 期刊上，标题为《Anisotropic Rashba coupling to polar modes in KTaO₃》。
近年来，在KTaO₃（KTO）界面及离子门控掺杂的表面发现了一类新型超导体，其二维电子气的超导临界温度（Tc）比同类体系SrTiO₃（STO）高出一个数量级，且强烈依赖于晶向。这一发现激发了学界对其配对机制的巨大兴趣。此前有研究提出，电子通过与软横光学声子（TO）模式，即铁电（FE）模式的线性耦合来介导库珀对，或许能解释Tc的强方向依赖性。
这种耦合机制被称为Rashba耦合，它源于自旋-轨道耦合（SOC）辅助下电子与极性模式的相互作用。在STO中，该机制已被证实能提供足够的BCS耦合常数并解释Tc随载流子浓度的穹顶状变化。基于这些成果，本研究旨在深入探索KTO中三个SOC能带与极性模式之间的线性耦合，超越先前研究的粗略估算，为理解KTO基异质结中的超导机制提供更精确的微观理论依据。
1. DFT电子结构参数化与模型构建 研究团队首先利用第一性原理（ ab initio ）密度泛函理论（DFT）计算了块体KTO的电子能带结构（图1a），并使用包含原子SOC的紧束缚（tight-binding）模型来拟合其低能色散部分。该模型纳入了Ta原子5d轨道的t₂g子空间（轨道$\mu = x, y, z$），并考虑了最近邻和次近邻的跃迁项。通过拟合参数，包括跃迁能量$t_1$至$t_4$和有效原子SOC参数$\xi$，模型在低能区域与DFT计算结果高度吻合（图1b虚线）。在布里渊区中心，SOC将t₂g六重简并态分裂为能量较低的$J=3⁄2$四重态（公式6,7）和能量高出400 meV的$J=1⁄2$双重态（公式8）。这一模型是后续分析电子与极性模式耦合的基础。
2. 极性软模的对称化坐标描述 为研究极性声子，研究引入了适用于立方钙钛矿的对称化坐标，将任意极性畸变分解为三个属于$T_{1u}$不可约表示的独立模式：Slater模式（$\bar{s}_1$，即Ta原子与氧笼的相对振动）、钾原子与刚性TaO₃笼的反相振动（$\bar{s}_2$）以及氧笼的形变（$\bar{s}_3$）。这些模式的振幅分别由$u_1$、$u_2$、$u_3$描述（公式9-11），其归一化保证了质量加权位移中心不变。
3. 通过“冻结声子”计算提取线性Rashba耦合 这是研究的核心实验部分。团队采用了“冻结声子”方法，即在DFT计算中人为引入一个静态的、波矢q=0的极性畸变$u_i$（沿[001]方向），然后计算并重构由此导致的电子能带分裂。 首先，根据对称性，一个极性畸变与SOC电子的线性耦合在动量空间中表现为一种类似Rashba的效应，可导致一个与动量k线性相关的能带分裂$\delta e_n(k)$。对于传统的各向同性Rashba模型，该分裂在垂直于极性轴的平面内具有轴对称性（公式14, 15）。然而，对三个$\bar{s}_i$模式的DFT冻结声子计算结果（图2a-c）显示：只有最高能带（$n=3$, $J=1⁄2$双重态）符合各向同性Rashba描述；而较低的两个能带（$n=1, 2$, $J=3⁄2$四重态）的能带分裂表现出强烈的各向异性，其沿[110]和[100]方向的线性斜率明显不同。 为解释这一现象，研究者在微观耦合矩阵中加入了更高阶的角动量修正项——$l=3$的立方谐波$k_3(\hat{k})$（公式16, 17）。通过用新的表达式（包含$\tau^{si}{n,1}$和$\tau^{si}{n,3}$两个参量）拟合DFT数据，成功量化了各能带与各极性模式的各向同性及立方谐波修正耦合强度（表1）。
4. 基于微观玩具模型的机理探索 为阐明各向异性的起源，研究构建了一个微观玩具模型。该模型在原有的紧束缚哈密顿量$H0+H{SOC}$基础上，增加了一项由极性位移诱导的、奇宇称且自旋守恒的轨道间跃迁项$H^{mic}_u$（公式21, 22），其强度由一个参数$\partial t / \partial u$控制。 通过对该模型进行简并微扰理论分析，发现其核心机制在于：非相互作用的电子跃迁项$H_0$会以依赖于动量方向$\hat{k}$的方式打破$J=3⁄2$多态在$\Gamma$点的简并。当将$H^{mic}_u$投影到这个与$\hat{k}$相关的微扰基上时，自然导致了$J=3⁄2$能带的Rashba耦合呈现出角度依赖性，而$J=1⁄2$能带因缺乏额外简并度则保持各向同性。这一结果成功再现了DFT计算中各向异性出现与否的定性特征（图3e），并给出了沿[100]和[110]方向分裂的解析表达式（公式24-26）。
5. 电子-声子耦合哈密顿量量化与超导配对强度估算 最后，研究对极性声子进行了量子化，推导出电子-声子耦合函数$g^\alpha_n(k,q)$（公式29）。通过假设软TO模是纯粹的某种$\bar{s}i$模式，并结合实验观测到的低温声子频率$\omega{TO} = 2.5$ meV，估算出其特征零點運動長度和各模式的电子-声子矩阵元（表1, 图3右侧坐标）。结果显示，Slater模式（$\bar{s}_1$）的耦合最强（约100 meV）。在此基础上，通过取各向同性极限并对s波超导配对进行BCS理论估算，发现配对强度与STO体系极为相似，表明掺杂的块体KTO是探索超导电性的一个有前景的方向。
本研究系统地揭示了立方KTaO₃中电子与极性模式耦合的高度各向异性特征，其核心结论与价值如下： 1. 发现了各向异性Rashba耦合：首次通过DFT计算证实，KTO中$J=3⁄2$能带的电子与软FE模式的耦合远非传统的各向同性Rashba模型所能描述，必须引入高阶（$l=3$）立方谐波修正，并定量提取了耦合参数。 2. 阐明了各向异性的微观起源：提出的微观玩具模型清晰地表明，各向异性源于动量依赖的电子跃迁项对$J=3⁄2$能带简并度的打破。这一发现具有普适性，暗示在任何具有高简并度SOC多重态的材料（如立方STO、半赫斯勒材料）中均可能存在类似的耦合各向异性。 3. 揭示了KTO的强大配对潜力：估算表明，KTO中的Rashba配对相互作用强度与已被广泛研究的STO相当，这为在掺杂的块体KTO中寻找超导电性提供了强有力的理论支持。 4. 为理解异质结超导提供了新视角：研究指出的块体材料中固有的各向异性Rashba耦合，为理解KTO异质结中Tc对晶体取向的极端敏感性提供了关键线索。界面处的额外晶体场效应将与此各向异性耦合共同作用，这是未来研究的重要方向。 5. 方法学贡献：本工作建立了一套将DFT冻结声子计算、对称性分析、紧束缚模型和微扰理论相结合的系统性方法，用于精确提取和解析材料的电子-声子相互作用，对相关领域的研究具有重要的方法学参考价值。