SrTiO3中轴向EPR晶体场的起源：表面诱导的声子-弹性各向异性

《Physica B: Condensed Matter》：Surface-induced phonon-elastic anisotropy as the origin of axial EPR crystal fields in SrTiO 3

　　•提出了一种由边界引起的有限SrTiO3中轴向EPR场的机制。•将自由表面视为一个弹性边界值问题，而非受损的体材料。•该模型解释了SrTiO3板材中的(001)、(110)和(111)取向序列。•有限厚度效应和软晶格效应可解释立方型和轴向EPR分量。•该方法表明EPR可作为探测埋藏氧化物界面的局部探针。

引言
钛酸锶SrTiO3是一种高度极化的钙钛矿氧化物模型：在室温下为立方结构，对多种结构及极性不稳定性较为敏感，并且在Tc≈105 K附近存在著名的反铁畸变（AFD）相变[1]、[2]、[3]。在Tc以上，三个立方轴〈100〉在晶体学上是对等的。然而，电子顺磁共振（EPR）并非反映理想体材料的晶体学相位，而是反映顺磁杂质处的实际局部晶场。这一差异引发了一个纯粹的理论问题，与具体实验无关：真实的有限SrTiO3晶体的弹性及振动（声子）环境如何影响此类探针所检测到的晶场，尤其是边界会产生何种影响？

本文的核心观点是，自由表面并非体材料的缺陷版本，而是一个不同的弹性边界值问题[4]、[5]、[6]。即便是在理想无应力状态下，表面也会改变正常模式谱、局部声子态密度[7]以及应变相关性，且这种变化具有各向异性：表面法向方向与两个切向方向的驻波模式谱各不相同。当将这些效应投影到SrTiO3的立方轴上时，由边界引起的声子-弹性各向异性会生成一个局部四方对称性不变量，而这正是EPR自旋哈密顿量中轴向晶场项所敏感的对称性。从自由边界→驻波模式→正常应变与切向应变谱密度不等→Qz≠0→B20≠0，这一系列推导是该最小模型中边界条件所带来的对称性结果。本研究旨在将这一链条量化，解决一个能够体现这一现象的最简但明确的边界问题，并阐明其可验证的结论。

为检验这些预测，需要获取关于经过良好表征、具有不同晶体学取向、厚度、表面处理方式及形状的单晶SrTiO3样品的系统性EPR数据。这类数据是存在的。Gabbasov等人曾报道过在经过结构表征的定向单晶板中Fe3＋和Mn4+中心的EPR信号[8]、[9]；我们在此将其作为独立的、高质量的标准参照，用以验证理论并进一步拓展研究。相关实测结果构成一个连贯的整体，而非孤立的异常现象：(i) 在Tc以上，尽管理想体材料为立方结构，但EPR谱仍需包含轴向项B20；(ii) 两种不同的探针Fe3＋和Mn4+均显示出该效应，表明这是材料本身的特性而非某种离子的特性；(iii) 在(001)取向的板材中，法向[001]轴与平面内的轴不等价，而两个平面内的立方轴则依然对等；(iv) 在(110)取向的板材中，所选的四方对称轴为位于板材平面内的立方[001]轴，而非表面法向轴；(v) 在(111)取向的板材中，EPR对称性在可检测范围内仍保持立方对称性。同一组数据还表明，随着厚度减小，轴向分量会增强；较薄的样品会呈现双组分（立方型加轴向型）EPR谱；精细抛光会使轴向参数相对于粗磨后的数值有所降低；不对称的表面处理会扩大响应范围；此外，板材与棒材在EPR信号特征上也存在差异，包括符号方面。

这一机制具有更重要的意义，原因有二。首先，它能够做出超出现有数据的预测，包括不同取向之间的比值符号与大小，以及轴向参数和双组分转变的温度依赖性。其次，自由表面只是内部边界的极端情况（即真空覆盖层），同样的框架也可应用于埋藏界面，此时轴向响应由界面两侧的声阻抗差异决定。这表明EPR可作为探测氧化物异质结构界面的局部、对称性敏感型探针，可与能够解析结构与密度分布的X射线方法相互补充。本文的结构如下：首先建立理论框架，将EPR轴向参数视为由声子-弹性环境产生的四方对称性不变量的局部探针；接着介绍弹性格林函数、应变谱密度以及边界条件这些核心物理要素；随后明确求解最简半空间问题，展示自由边界如何产生非零的表面诱导四方对称性不变量，再将该不变量投影到SrTiO3的立方轴上，从而得出(001)、(110)和(111)表面对应的取向序列；之后分析有限厚度效应、软晶格放大效应，以及单组分和双组分EPR谱的出现机制；接着将理论结论与已有的SrTiO3板材实验观测结果进行对比；最后讨论可验证的预测、该机制在其他软晶格钙钛矿材料中的潜在应用，以及其在埋藏界面和氧化物异质结构研究中的前景。

章节节选

理论框架
在引入方程之前，先明确本文中使用的符号含义。x、y、z表示SrTiO3的立方晶体学轴，ex、ey、ez为相应的单位向量。向量n表示自由表面的外法向。上标“bulk”指的是无限大、具有平移对称性的晶体。上标“hs”则表示半空间几何结构，即晶体占据z>0区域，其自由表面位于z=0处；这只是一个标签，并非数值。

表面弹性问题及其后果
我们现在求解一个最简模型，以明确展现表面诱导的轴向不变量。考虑z>0半空间中的标量弹性场u(r,t)。该模型虽不能替代完整的矢量弹性理论，但能够捕捉关键的边界物理现象：自由表面会用驻波模式取代传播的体波。对于无应力边界，其标量形式的边界条件为?zu|z=0=0。正常模式可表示为uq,k(ρ,z)=2eiq?ρcos(kz)，其中k>0，ρ=(x,y)，q=(qx,qy)为平面内的坐标。

有限厚度与双组分谱
对于厚度为h的板材，假设两个大面会产生一个轴向源Q(z)=Q0f(z)+Q0f(h?z)，其中f(z)描述了表面修饰后的声子-弹性各向异性的渗透程度。在单组分情况下，观测到的EPR参数近似为B20ˉ(h)=1h∫0hB20(z)dz。若有效渗透长度为ξeff，且h?ξeff，则B20ˉ(h)∝ξeffh。而当h?2ξeff时，两个受表面控制的区域会重叠，此时轴向响应几乎可遍及整个体积。在双组分谱中，则会出现不同的情况。

与SrTiO3板材EPR测量的比较
现在将上述理论结论与已有文献中关于具有特定取向、经结构表征的单晶SrTiO3板材中Fe3＋和Mn4+中心的EPR数据进行了比较，这些数据反映了取向、厚度、表面处理方式及形状等因素对EPR信号的影响[8]。首先来看取向依赖性。无需考虑低温下的AFD相变，仅凭自由表面就能在受表面控制的区域内的中心点产生非零的Qzsurf，进而导致非零的B20，这正是上文提到的105K以上出现的轴向项。

可验证的预测与未来展望
上述比较所依据的是已有的数据。同样的理论框架还能做出超出该数据范围的预测，可用于验证或排除该机制。最直接的预测涉及(001)和(110)取向。根据ΔΦ<0[方程(28)]，最简模型预测这两种取向之间的轴向参数符号会发生反转，同时有B20,(110)/B20,(001)=?1/2[方程(40)]这一关系。从对称性层面来看，每种取向会选择不同的轴，而(111)取向则存在节点现象。

结论
基于EPR探测的是顺磁杂质处的局部晶场，而非理想体材料的晶体学相位这一事实，我们研究了在有限尺寸的SrTiO3中，边界如何改变此类探针所感知的声子-弹性环境。核心观点是，自由表面是一个独立的弹性边界值问题，而非体材料中的缺陷区域。在最简的半空间模型中，无应力边界会取代在体中传播的波。

作者贡献说明
A.G. Kiiamov：写作——审稿与编辑，写作——初稿撰写，结果验证，实验研究，形式分析，概念构建。R.F. Mamin：写作——审稿与编辑，研究指导，项目管理。

利益冲突声明
作者声明不存在任何可能影响本文研究的已知财务利益或个人关系。

致谢
本研究得到了俄罗斯科学院喀山科学中心FRC项目的政府资助支持。

A.G. Kiiamov|R.F. Mamin