PbMnO3:PbVO3 超晶格中的多铁性与可切换 Berry 曲率偶极子PHYSICAL REVIEW B 113, 214114 (2026)PbMnO₃:PbVO₃ 超晶格中的多铁性与可切换 Berry 曲率偶极子First-Principles Prediction of a Multiferroic Semiconductor with Switchable Spin-Polarized Berry Curvature Dipole导读 导读通过第一性原理计算预测 (PbMnO₃)₁/(PbVO₃)₁ 超晶格是一种同时具有铁磁性、铁电性和 Berry 曲率偶极子的多铁性半导体。V⁴⁺→Mn⁴⁺ 电荷转移驱动极性畸变和 Jahn-Teller 能隙产生的自旋极化 Berry 曲率偶极子幅度与含重元素 Weyl 半金属相当。结合磁电耦合可实现电控四态非易失性 Hall 存储器件。一、引言多铁性与量子几何的交汇1.1 多铁性材料的稀缺性多铁性固体同时打破空间反演P和时间反演T对称性展现出线性磁电耦合和不同有序态之间的非线性切换。然而多铁性在自然界极为罕见铁电体需要 d⁰ 构型的 B 位离子而铁磁体需要 dⁿ (n≠0) 构型二者存在根本性矛盾。近年来电子的量子几何Berry 曲率、Berry 曲率偶极子为反常霍尔效应等拓扑输运现象提供了新视角。1.2 PbMnO₃:PbVO₃ 超晶格的设计思路PbMnO₃ 是半金属铁磁体Pnma 空间群其异质结 SrMnO₃ 通过混合非本征铁电机制展现极性。本文策略用 V 替代 Mn——V⁴⁺(d¹) 向 Mn⁴⁺(d³) 转移一个电子形成 Mn³⁺(d⁴, JT 活性) 和 V⁵⁺(d⁰, 铁电活性)。Jahn-Teller 畸变打开能隙V⁵⁺ 位移驱动铁电极化实现半导体性铁磁-铁电多铁态同时承载自旋极化 Berry 曲率偶极子。二、计算方法与磁基态2.1 第一性原理计算设置采用 Quantum ESPRESSO 进行 DFT 计算PBEsol 泛函PAW 赝势。平面波截断能 60 Ry波函数/480 Ry电荷密度。对 Mn 3d 态施加 Hubbard U 3 eV对 V 3d 态同样施加 U 3 eV以重现 PbVO₃ 的实验极性四方 P4mm 结构。考虑五种磁构型铁磁FM、亚铁磁FIM、C-AFM、G-AFM 和 Mn-AFM/V-FM 混合构型。2.2 磁基态与结构相变四方 P4/mmm 结构中铁磁态能量最低比 FIM 和 AFM 态低约 35 meV/f.u.。P4/mmm 结构表现出强烈的不稳定性Γ 点极性模Γ₅⁻和 M 点 JT 畸变模同时软化。凝聚这两个不稳定模后体系弛豫至极性 Cc 空间群Mn³⁺ 和 V⁵⁺ 分别呈现 JT 活性八面体和铁电位移。图1: (a) P4/mmm 四方结构。(b) 声子不稳定模。(c) 极性 Cc 相的结构畸变。三、电子结构起源极性-铁磁-半导体3.1 电荷转移与 JT 能隙Bader 电荷分析确认V 向 Mn 转移约 0.7 e⁻形成 Mn³⁺(d⁴) 和 V⁵⁺(d⁰)。Mn³⁺ 的 eg 轨道dz² 和 dx²−y²在 JT 畸变下发生劈裂dz² 被占据而 dx²−y² 为空打开约 0.6 eV 的能隙。V⁵⁺(d⁰) 的位移驱动极性 Γ₅⁻ 畸变产生约 15 μC/cm² 的自发极化。Jahn-Teller 畸变能量Q₂ 和 Q₃ 为 JT 活性坐标图2: 铁电相的起源。(a) 能量随极性模振幅的变化。(b) Mn-O 键长不对称性。(c) 结构序参量。3.2 轨道物理与磁性Mn³⁺(d⁴) 的 eg 轨道占据导致 S 2 的高自旋态铁磁交换通过 Mn-O-Mn 超交换路径实现。V⁵⁺(d⁰) 无局域磁矩不参与磁交换仅通过结构畸变间接影响电子结构。能带计算显示费米能级附近态密度主要由 Mn eg 轨道贡献V 的态密度远在费米能级以上。图3: 电子结构。(a) 自旋极化能带。(b) 轨道投影态密度。(c) 实空间电荷密度分布。四、Berry 曲率偶极子与非线性霍尔效应4.1 Berry 曲率偶极子的计算使用 Wannier90 构建 Mn eg 轨道的最大局域化 Wannier 函数基于 320×320×220 k 点网格计算 Berry 曲率偶极子。Berry 曲率偶极子 Dₐᵦ 的幅度与含重元素强 SOC Weyl 半金属相当尽管 PbMnO₃:PbVO₃ 仅含 3d 轻元素。Dₐᵦ 的自旋极化特性显著多数自旋通道的贡献远大于少数自旋通道。Berry 曲率偶极子BCD定义Berry 曲率的一般表达式4.2 非线性反常霍尔效应与四态存储非线性横向霍尔电压 jₐ σₐᵦᶜ Eᵦ Eᶜ其中 σₐᵦᶜ 与 Berry 曲率偶极子 D 成正比。铁电极化翻转可切换 Berry 曲率偶极子的符号实现电控非线性霍尔效应。结合自旋极化和极化方向可实现四态非易失性存储{P↑, P↓, −P↑, −P↓}在单个霍尔器件中编码四位信息。非线性反常霍尔电流响应Berry 相位极化铁电极化的现代理论表达图4: 量子几何的电学检测。(a) Berry 曲率分布。(b) Berry 曲率偶极子。(c) 非线性霍尔响应。(d) 四态存储方案。五、补充信息Supporting InformationS1-S6: 磁构型与结构畸变图 S1-S2: 五种磁构型示意图及极性 P4mm 相的结构参数。图 S3-S4: 声子色散关系及极性结构畸变的详细分析。图 S5-S6: Mn-O 和 V-O 键长不对称性揭示 JT 畸变和铁电位移的实空间特征。图 S1: 磁有序构型图 S2: 极性 P4mm 相图 S3: 声子色散图 S4: 极性结构畸变图 S5: Mn-O 键长图 S6: V-O 键长S7-S13: 电子结构与 Wannier 函数图 S7-S10: 电子结构、轨道投影 DOS 和多数自旋态。图 S11-S13: 电子-声子耦合、Wannier 函数和跃迁参数。图 S7: 电子结构图 S8: 轨道投影 DOS图 S9: 轨道投影 DOS (2)图 S10: 多数自旋态图 S11: 电子-声子耦合图 S12: Wannier 函数图 S13: 跃迁参数S14-S23: Berry 曲率与铁电切换图 S14-S18: Berry 曲率分布、Berry 曲率偶极子计算和 Wannier 电荷中心。图 S19-S23: 铁电切换路径、电子结构对比、磁矩翻转、自旋劈裂和 AFM 能量。图 S14: Berry 曲率图 S15: Berry 曲率 (2)图 S16: 电子结构-Berry图 S17: Berry 曲率偶极子图 S18: Wannier 电荷中心图 S19: 铁电切换图 S20: 电子结构对比图 S21: 磁矩翻转图 S22: 自旋劈裂图 S23: AFM 能量DFT 计算经验贴士【DFT Tip 1】多铁性 DFT 计算空间反演时间反演双重破缺多铁性体系同时破缺 P 和 T 对称性DFT 计算中必须① 不施加任何对称性约束ISYM0② 允许所有原子自由弛豫ISIF3③ 分别计算铁磁和多种反铁磁构型的能量以确定磁基态。常见错误① 在 P4/mmm 高对称结构中弛豫→得到亚稳态而非极性基态② 用对称性加速计算ISYM1/2→掩盖了极性畸变③ 忽略 Berry 相位极化计算中的极化量子问题。【DFT Tip 2】Hubbard U 的双金属设置Mn 和 V 的 U 值校准超晶格中 Mn 和 V 都需要 U 值。Mn 3d 的 U3-4 eV 是标准选择但 V 3d 的 U 需要特别注意V⁵⁺(d⁰) 的 U 值应使 PbVO₃ 再现实验极性 P4mm 结构。方法对 V 的 U 从 2 到 5 eV 扫描检查 PbVO₃ 的声子不稳定模Γ₅⁻和极性畸变幅度。U 太小→V 金属化→无极性畸变U 太大→V 过局域→假相。校准标准U3 eV 时 PbVO₃ 的 c/a 比应与实验一致。【DFT Tip 3】Jahn-Teller 畸变检测八面体参数量化在 VASP 中检测 JT 畸变① 计算 Mn-O 八面体的键长分布——JT 活性八面体有 2 长4 短或 2 短4 长的键长模式② 八面体畸变参数 Δd (1/6)∑[(di-davg)/davg]² × 10⁴Δd 100 为显著 JT 畸变。关键自旋极化对 JT 畸变有重要影响——先收敛磁基态再做结构弛豫。常见错误用非自旋极化计算弛豫 JT 活性体系→JT 畸变被低估。【DFT Tip 4】Bader 电荷分析电荷转移量化的模拟选择Bader 分析中Mn³⁺ vs Mn⁴⁺ 的 Bader 电荷差约 0.3-0.5 e⁻而非 1 e⁻因为 Mn-O 键的共价性回馈了部分电荷。V⁴⁺→V⁵⁺ 的氧化态变化同样被共价性掩盖。建议① 不要仅依赖 Bader 电荷的绝对值② 关注 Bader 电荷的相对变化Mn 在 PMO vs PMO-PVO 超晶格中的差异③ 结合磁矩分析Mn³⁺(d⁴) 的磁矩 ~4 μBMn⁴⁺(d³) 的磁矩 ~3 μB交叉验证氧化态。【DFT Tip 5】Berry 相位极化计算极化量子与分支选择铁电极化 P 通过 Berry 相位计算P (e/Ω)Σ Z*_κ u_κ。VASP 中设置 LCALCPOL.TRUE.。关键问题极化量子——P 和 P eR/ΩR 为晶格矢量在物理上等价。操作用中心对称参考相P4/mmm→极性相Cc的绝热路径计算极化差 ΔP避免极化量子歧义。常见错误直接比较两个极性相的极化绝对值→可能相差一个极化量子。【DFT Tip 6】Wannier90 拟合投影轨道选择与展开窗口Mn eg 轨道dz² 和 dx²−y²是 JT 畸变和 Berry 曲率偶极子的核心。Wannier90 投影① 投影轨道Mn:l2d 轨道② 外窗能量范围 [EF-2 eV, EF2 eV]③ 内窗冻结窗口覆盖费米能级附近 ±1 eV。关键Wannier 展宽需 10 Ų对于 4 个 Wannier 函数/超胞否则 Berry 曲率偶极子被人工平滑。检查Wannier 拟合的能带与 DFT 能带在费米能级附近完全重合。【DFT Tip 7】Berry 曲率偶极子BCD计算K 点与温度展宽BCD ∫(∂E/∂k_a) Ω_b(k) (-∂f/∂E) d³k 对 K 点要求极高。Wannier-Berri 代码中典型设置K 点 300×300×200温度展宽 60 K。BCD 的收敛性检查① 依次增大 K 点网格100³→200³→300³→400³② 检查符号和幅度的收敛趋势③ 注意 BCD 仅在被 SOC 打破镜面对称的体系中非零。常见错误K 点不足→BCD 噪声大于信号。【DFT Tip 8】量子 ESPRESSO 中的多铁性计算与 VASP 的关键差异QE 中多铁性计算的特殊参数① 极化计算berry_phase 模块需串行② 声子计算DFPT中极性体系需要 LO-TO 劈裂修正epsil.true.③ Hubbard U 的 DFTUV 需要指定原子间 V 参数。QE vs VASP 的差异① QE 的 PAW 数据集与 VASP 的 PAW 可能不同→U 值需要重新校准② QE 的 Berry 相位计算在 nscf 中完成而非 scf③ QE 的声子计算使用 DFPT而非 VASP 的有限差分。【DFT Tip 9】非线性霍尔效应的弛豫时间近似τ 的选择非线性霍尔电导率 σ_abc ∝ τ弛豫时间τ 是经验参数。在 DFT 中 τ 通常取为常数 τ 10 fs对应室温电导率 ~10⁵ S/m。但 τ 对 σ_abc 有线性影响——τ 的选择直接影响 BCD 效应的实际幅度。建议① 从实验电导率反推 ττ m*σ/ne²② 做 τ 敏感性分析τ5, 10, 20 fs③ 在比较不同体系时使用相同的 τ 以确保可比性。【DFT Tip 10】四态存储器件模拟铁电极化翻转路径铁电极化翻转通过 NEB 或 SS-NEB 计算。路径P↑ → 中心对称P0→ P↓。关键参数① 翻转能垒 E_barrier——决定翻转电压和保持时间② 翻转路径中是否经过中间相如非极性相③ 翻转变分对磁性的影响磁矩是否同时翻转。在 VASP 中使用 VTST 工具包的 NEB 模块5-7 个中间图像Spring constant -5。常见错误NEB 图像数不足→翻转路径不光滑高估能垒。超晶格设计PMO:PVO1:1电荷转移V⁴⁺→Mn⁴⁺JT畸变Mn³⁺: eg¹dz²极性畸变V⁵⁺: d⁰位移0.6eV磁基态5种构型对比声子稳定性Γ₅⁻M点软化Wannier90Mn eg轨道BCD计算320×320×220Dab非零非线性HalljaσabcEbEc自旋极化多数少数电控切换P↑↔P↓四态存储{P↑↓, -P↑↓}4位存储▸ 设计策略PbMnO₃(半金属FM)PbVO₃→电荷转移 V⁴⁺→Mn⁴⁺→Mn³⁺(JT)V⁵⁺(d⁰极性)▸ 物理机制JT畸变打开能隙→极性畸变产生P→SOC低对称性→Berry曲率偶极子→非线性Hall▸ 器件应用铁电极化翻转→BCD符号切换→自旋极化四态存储→氧化物自旋电子学知识扩展【知识扩展 1】Berry 曲率偶极子BCD与非线性霍尔效应量子几何的新维度【理论解释】Berry 曲率偶极子 D_ab ∫(∂E/∂k_a) Ω_b(k) (-∂f/∂E) d³k 是 Berry 曲率在动量空间中的一阶矩——类似电荷分布的偶极矩。它产生二阶非线性霍尔效应 j_a σ_abc E_b E_c其中 σ_abc ∝ D_bd。Sodemann Fu2015首次提出 BCD 概念揭示了时间反演对称性破缺体系中非线性输运的量子几何起源。【方法比较】BCD 与线性 AHE 的对比(1) 线性 AHE——需要净 Berry 曲率∫Ω ≠ 0仅存在于破缺 T 的体系(2) BCD——需要 Berry 曲率的偶极分布∫(∂E/∂k)Ω ≠ 0同时需要 P 和 T 破缺或特定的镜面对称性破缺。BCD 的优势即使在 Berry 曲率积分为零的体系中如某些反铁磁体BCD 仍可非零。【经典参考】Sodemann Fu, Phys. Rev. Lett. 115, 216806 (2015)——BCD 的概念提出Ma et al., Nature 565, 337 (2019)——WTe₂ 中 BCD 的实验验证Zhang et al., Nature 621, 487 (2023)——MnBi₂Te₄ 中量子化 BCD。【迁移能力】BCD 分析不仅适用于多铁性氧化物也适用于Weyl 半金属TaAs 家族、转角石墨烯Moiré BCD、二维铁电体In₂Se₃、和反铁磁体MnBi₂Te₄。本文的核心创新——将 BCD 与多铁性结合→电控 BCD 切换——可迁移至任何铁电磁性共存体系。【知识扩展 2】Jahn-Teller 效应与轨道有序从配位化学到多铁性设计【理论解释】Jahn-TellerJT定理1937非线性分子在简并电子基态下会发生几何畸变以降低能量。在八面体场中eg 轨道简并dz² 和 dx²−y²→轴向伸长Q₃0或压缩Q₃0。JT 效应的能量增益E_JT ½K(Q₂²Q₃²) λQ₃(Q₃²-3Q₂²)其中一次项线性 JT仅在 d⁴(HS), d⁷(LS), d⁹ 中出现。【方法比较】JT 畸变的三种检测方法(1) 结构参数——八面体键长分布Δd和键角方差σ²(2) 电子结构——轨道分辨 PDOS 中 eg 轨道的分裂Δeg ~ 0.5-1.5 eV(3) 能量分析——JT 稳定化能 E_JT_stab E(畸变) - E(正则)。Mn³⁺d⁴ HS的 JT 稳定化能通常 ~0.3-0.8 eV远大于其他 d 电子构型。【经典参考】Jahn Teller, Proc. R. Soc. A 161, 220 (1937)——原始定理Bersuker, Chem. Rev. 101, 1067 (2001)——JT 效应综述Kugel Khomskii, Sov. Phys. Usp. 25, 231 (1982)——轨道有序与 JT 效应。【迁移能力】JT 效应不仅是 Mn³⁺ 的专利——任何 d⁴/d⁷(HS) 或 d⁹ 体系都有 JT 活性。在材料设计中可主动利用 JT 畸变(1) 打开能隙金属→半导体(2) 增强磁各向异性通过畸变→SOC(3) 驱动轨道有序Kugel-Khomskii 机制(4) 与铁电畸变耦合→多铁性。本文的 Mn³⁺V⁵⁺ 设计是 JT 效应铁电耦合的典范。科研经验【科研经验 1】为什么 Berry 曲率偶极子的符号和幅度难以收敛问题BCD 计算中K 点网格从 100×100×100 增加到 300×300×300 时BCD 符号翻转或幅度变化超过 50%。原因① BCD 是费米面上的积分——费米面拓扑结构对 K 点极其敏感费米面的微小变化可能导致 BCD 符号翻转② Wannier 拟合中内窗/外窗的能量范围选择影响费米能级附近的能带精度→BCD 的符号可能由能带交叉的细节决定③ 在窄带隙体系中Eg 0.5 eV能带顺序的微小变化可能导致 BCD 符号翻转。解决方案① 使用极密 K 点网格400×400×400 或更高② 做 Wannier 拟合的收敛性测试——改变投影轨道和能量窗口检查 BCD 的稳定性和符号一致性③ 计算 BCD 的费米面分解——识别哪些费米面区域主导 BCD④ 用 Fermi sea 贡献而非仅 Fermi surface 贡献交叉验证。建议BCD 的符号是物理可观测的决定非线性 Hall 电压的符号但 DFT 计算中 BCD 的符号可能对数值参数敏感。建议报告 BCD 的幅度作为主要结果符号作为辅助信息除非有明确的对称性论证。【科研经验 2】多铁性 DFT 计算中为什么声子不稳定但实验结构稳定问题DFT 声子计算显示 P4/mmm 结构有多个不稳定模Γ₅⁻ M 点但实验上 PbVO₃ 的 P4mm 结构是稳定的。原因① 声子不稳定模恰好是正确的不稳定——它们驱动体系向更稳定的极性相Cc 空间群弛豫② 在 DFTU 中U 值的选择直接影响声子频率——U 太小V 金属化→虚频消失但极性相也被抑制U 太大→虚频过大→高估不稳定性③ 非谐效应——简谐声子计算在 T0 K 下进行而实验在室温下观察。解决方案① 凝聚不稳定声子模→弛豫到低对称相→检查低对称相的声子是否全部稳定② 对 U 值做声子收敛测试U2,3,4,5 eV③ 用 AIMD 在 300 K 验证动力学稳定性。建议对于多铁性体系有虚频不一定是坏事——关键是虚频模式是否对应正确的物理畸变极性模、JT 模、氧八面体旋转。如果虚频模式与预期畸变一致→凝聚它们后得到正确的基态。从多铁性 BCD 到氧化物自旋电子学的完整研究链条本文已完成了 PbMnO₃:PbVO₃ 超晶格的 DFT 计算——结构、磁性、电子结构、BCD 和四态存储方案。但作为多铁性半导体平台以下计算可以进一步挖掘其物理内涵和器件潜力。Further 1 | 铁电翻转路径的 NEB 计算与翻转动力学为什么值得算本文提出了四态存储方案但铁电极化翻转的能垒和路径是实际器件工作的关键。用 CI-NEB 计算 P↑→P↓ 的翻转路径和能垒。能回答的问题翻转能垒是多少是否在可操作电压范围内 2 V翻转路径是否经过中间亚稳相磁矩在翻转过程中是否保持适合体系所有铁电磁性共存的多铁性体系特别是钙钛矿超晶格。输入极性相 P↑ 和 P↓ 的结构NEB 计算5-7 个中间图像翻转能垒。Further 2 | 磁电耦合系数线性与非线性磁电响应为什么值得算多铁性体系的核心是磁电耦合。计算线性磁电张量 α_ij ∂P_i/∂B_j 和非线性磁电系数。能回答的问题PbMnO₃:PbVO₃ 的磁电耦合强度是否与 Cr₂O₃经典磁电材料可比电场是否能翻转磁矩而非仅切换 BCD 符号适合体系所有多铁性体系特别是同时具有铁电和铁磁有序的体系。输入DFPT 计算响应电场和磁场Born 有效电荷磁电耦合张量。Further 3 | 光学性质与二次谐波SHG非线性光学的量子几何指纹为什么值得算BCD 和非线性光学响应如 SHG共享量子几何起源。计算 SHG 张量 χ^(2)(ω) 和光电流响应。能回答的问题SHG 信号是否能区分 P↑ 和 P↓ 态SHG 的偏振依赖性是否与 BCD 的对称性一致SHG 能否作为 BCD 的光学探针适合体系所有非中心对称体系特别是具有大 BCD 的多铁性半导体。输入SOC 自洽的 WAVECAR已完成VASP 非线性光学计算LOPTICS SHGK 点密度 ≥ 15×15×8。Further 4 | 替代效应Pb→Sn/Ba 和 Mn→Fe/Cr 的化学调控为什么值得算A 位Pb→Sn/Ba和 B 位Mn→Fe/Cr, V→Nb/Ta替代可以调控 JT 畸变、极化和 BCD。计算不同替代组合下的电子结构、磁性和 BCD。能回答的问题Sn 替代是否增强极性畸变Fe 替代是否能引入额外的 JT 活性哪种替代组合最大化 BCD适合体系所有 ABO₃ 钙钛矿超晶格特别是含 JT 活性离子的体系。输入不同替代组合的弛豫结构电子结构计算Wannier90 拟合BCD 计算。Further 5 | 超晶格周期调控1:1→2:1→1:2 的维度效应为什么值得算超晶格周期PMO:PVO 的比例改变了 MnO₂ 和 VO₂ 层的相对厚度→影响电荷转移、JT 畸变和界面效应。计算 (PMO)n/(PVO)m 的不同周期组合。能回答的问题2:1 超晶格的 BCD 是否更大界面效应Mn-O-V 键是否增强或削弱 BCD最优周期组合是什么适合体系所有氧化物超晶格特别是界面驱动的功能性体系。输入不同周期组合的弛豫结构超胞电子结构和 BCD 计算。Further 6 | 温度效应Monte Carlo 模拟磁相变与铁电相变为什么值得算本文所有计算在 T0 K 下进行。用 Monte Carlo 模拟基于 Heisenberg 模型铁电 Ginzburg-Landau 模型预测磁相变温度 Tc 和铁电相变温度 TFE。能回答的问题磁有序和铁电有序是否在相同温度下出现是否存在磁电耦合增强的临界温度窗口室温下四态存储是否可行适合体系所有多铁性体系特别是磁有序和铁电有序温度相近的体系。输入磁交换参数 J_ij需 TB2J 提取铁电双阱势从 NEB 计算MC 模拟代码。A. Paul, U. V. Waghmare | Physical Review B 113, 214114 (2026) | 多铁性 · Berry 曲率偶极子 · 非线性霍尔效应 · 四态存储