第一性原理模拟研究铁电钙钛矿中局域畸变对X射线吸收的影响
主要作者与发表信息
本研究由加利福尼亚大学圣迭戈分校（University of California San Diego）纳米与化学工程系的Pedram Abbasi、David P. Fenning和Tod A. Pascal共同完成，其中Tod A. Pascal为通讯作者。该论文于2022年10月15日投稿，2023年2月10日被接收，并于2023年2月13日发表于《Nanoscale》期刊第15卷第5193至5200页。
学术背景与研究目的
铁电材料在传感器、非易失性存储器等电子与能量转换器件中有着广泛应用，其中很大一部分属于具有ABO₃通式的钙钛矿（perovskite）家族。以钛酸钡（BaTiO₃）为代表的钛酸盐铁电体中，铁电性源于Ti⁴⁺阳离子从其中心对称位置向两个亚稳态构型的位移，即正向和负向极化。这种极化翻转（polarization switching）效应源于钛-氧（Ti–O）键的杂化作用，深刻理解其电子结构对于阐明极化翻转效应至关重要。
X射线吸收光谱（X-ray absorption spectroscopy, XAS）是一种能够以高灵敏度探测选定元素电子结构的独特工具，对理解局域结构和短程相互作用导致的铁电极化效应具有重要意义。其中，氧K边（O–K edge, 1s→2p激发）和钛L边（Ti–L edge, 2p→3d激发）是最受关注的谱学特征。然而，由于实验方法在所需分辨率上的困难，以及对结果进行解释的难度，这些表面的纳米级形貌及其对外加电场的响应尚未得到适当阐明。特别是钛L边涉及自旋-轨道耦合和强多重态效应等更复杂的物理过程，使得该边的极化翻转研究更具挑战性。
在计算模拟领域，模拟XAS的理论方法可分为单粒子方法和多体方法两类。本研究旨在回答两个同等重要的问题：其一，利用实验和基于第一性原理的模拟XAS，在O–K边和Ti–L边研究极化翻转对BaTiO₃薄膜表面电子态的调控作用；其二，对两种引入多体物理的XAS计算方法进行基准测试（benchmark），一种是基于多体行列式Delta自洽场（many-body ΔSCF determinant）的方法，另一种是基于GW/Bethe–Salpeter方程（GW/BSE）的多体微扰理论方法。
研究方法与详细流程
本研究的实验数据来源于同步辐射X射线扫描隧道显微镜（Synchrotron X-ray Scanning Tunneling Microscopy, SX-STM）采集的外延生长（001）BaTiO₃薄膜的表面敏感Ti–L边和O–K边XAS谱。通过近场模式下光激发电子向STM针尖的隧穿效应，该方法能够以高灵敏度捕获表面极化翻转效应。
计算模拟方面，研究构建了具有向上和向下极化方向的BaTiO₃平板模型（polarized slab models）。为充分阐明局域极化对电子结构的影响，研究设置了两种不同程度的钛原子偏离中心对称位置的局域畸变，记作P₁和P₂，分别对应约18 μC cm⁻²和32 μC cm⁻²的极化强度。研究者系统地开展了两类计算方法的对比研究。
第一类方法是基于多体行列式的ΔSCF方法（mb-ΔSCF）。该方法属于约束占据的Delta自洽场密度泛函理论框架，其核心是将XAS末态处理为有效的电子-空穴对（electron-hole pairs）的叠加。研究中系统考察了价带空轨道数量（n×m，其中m为基态占据轨道数）和激发阶数（fᵃ，即模拟行列式谱时考虑的电子-空穴对数目）对模拟谱的收敛性影响。在O–K边计算中，研究者发现当n=3时谱图在吸收边以上20 eV范围内达到收敛，且二阶激发（f²）即可实现收敛。作为参考指标，研究还计算了基于初态近似的投影态密度（projected density of states, PDOS）谱，以及单粒子末态谱，用于与多体行列式谱进行对比。在Ti–L边模拟中，同样采用收敛参数（f², n=3）对不同极化方向和程度的平板模型进行计算。
第二类方法基于多体微扰理论的GW/BSE方法。该方法通过求解Bethe–Salpeter方程，包含准粒子能量以及核心空穴与激发光电子之间的相互作用，能够显式地考虑核心空穴（core-hole）和多重态（multiplet）效应。研究利用GW/BSE方法计算了极化平板模型的O–K边和Ti–L边XAS谱，并将结果与行列式方法及实验数据进行系统性对比。
主要研究结果
实验数据显示，Ti–L边XAS的特征为两对峰，分别位于456.4 eV和458.53 eV（属于L₃态的t₂g和e_g轨道），以及461.93 eV和464.06 eV（属于L₂态的t₂g和e_g轨道）。O–K边的t₂g和e_g精细结构特征峰分别位于529.65 eV和532.1 eV。整体而言，将极化方向从向下翻转向向上时，t₂g和e_g态的峰强度比（记作P⁺/P⁻）出现了可重复且显著的增加，其中以Ti–L₂边最为敏感。
在多体行列式ΔSCF方法的模拟结果中，O–K边谱的对比分析具有决定性意义：相比多体行列式谱，PDOS谱出现了约7 eV的红移（red shift），且完全缺失了关键的t₂g和e_g特征；单粒子XAS方法则未能捕获主要的低能t₂g精细结构，也无法再现实验的t₂g/e_g比值。而多体行列式方法显著改善了主t₂g峰的强度，与实验数据更为吻合。在极化效应方面，O–K边谱显示当极化从向下（P⁻）翻转向向上（P⁺）时，t₂g和e_g峰强度均增加，且P⁺/P⁻比值随畸变程度从P₁到P₂的增强而进一步增大，证实峰调制确实源于局域极化。然而，该方法在Ti–L边模拟中表现出显著局限性：仅能再现位于456.3 eV和458.4 eV的L₃态峰对，完全无法捕获L₂峰，其根本原因在于缺乏对自旋-轨道耦合（spin–orbit coupling）的描述，无法正确反映晶体场效应和d态的多重态分裂。此外，L₃峰的t₂g/e_g比值约为1.9，与实验值0.55相比出现反转。
GW/BSE方法的结果则展现出显著改进。在O–K边，低能t₂g峰强度相比行列式方法明显增强，更接近实验数据，且极化翻转效应更加显著。P₁极化下t₂g和e_g特征峰分别位于529.2 eV和532.2 eV，而P₂的更高畸变程度导致e_g峰蓝移约1 eV，高能峰线形也因多重散射效应而发生调制。最为关键的是Ti–L边模拟：BSE方法成功捕获了全部四个锐利峰——456.3 eV和458.5 eV归属于L₃的t₂g和e_g态，461.8 eV和463.9 eV归属于L₂的t₂g和e_g态，且L₂和L₃峰的t₂g/e_g比值均与实验更为吻合。极化从P₁增至P₂对Ti–L谱t₂g峰的影响最为显著，与实验发现一致。BSE方法通过捕获自旋-轨道耦合效应，显著改善了物理描述的完备性。
研究结论与价值
本研究通过实验和两种先进计算方法，系统探究了BaTiO₃铁电薄膜中晶格畸变和极化翻转对Ti–L边和O–K边X射线吸收的影响。多体行列式计算方法虽较单粒子方法有所改进，在O–K边实现了更好的t₂g/e_g比值并再现了向上极化表面的更高激发强度，但因缺乏自旋-轨道耦合而无法再现Ti–L边的L₂,₃分离。相比之下，多体BSE方法同时具备了自旋-轨道耦合和更优越的电子-空穴屏蔽与光激发动力学描述，在O–K边精细结构和Ti–L边L₂,₃分离及t₂g/e_g比值方面均与实验数据高度吻合。研究表明，对于包含强自旋耦合效应的系统（如4d和5d过渡金属），BSE方法是准确捕获物理本质的关键。
科学价值方面，该研究证明了模拟XAS是捕获复杂材料纳米级结构的强有力方法，并明确了多体微扰方法在描述核心空穴和多重态效应方面的重要性。应用价值方面，研究成果为利用第一性原理光谱方法作为独立和互补探针，研究复杂体系的结构和电子性质提供了方法学基础。
研究亮点
本研究的亮点体现在三个方面。其一，通过系统性基准测试，首次明确指出多体行列式ΔSCF方法虽在O–K边表现可接受，但在Ti–L边因自旋-轨道耦合缺失而存在原理性缺陷，而GW/BSE方法在此类问题上的显著优越性得到了实验验证。其二，成功建立了局域极化畸变程度（P₁与P₂）与XAS峰强度调制之间的定量关联，为通过谱学特征反推局域结构提供了依据。其三，同步辐射X射线扫描隧道显微镜（SX-STM）的实验数据与第一性原理模拟的紧密结合，代表了从纳米尺度理解表面极化和电子结构调控的有效策略。研究还指出了未来在4d和5d过渡金属等强自旋耦合体系中进一步应用的广阔前景。