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 Information）
S1-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 计算中必须：① 不施加任何对称性约束（ISYM=0），② 允许所有原子自由弛豫（ISIF=3），③ 分别计算铁磁和多种反铁磁构型的能量以确定磁基态。
常见错误：① 在 P4/mmm 高对称结构中弛豫→得到亚稳态而非极性基态，② 用对称性加速计算（ISYM=1/2）→掩盖了极性畸变，③ 忽略 Berry 相位极化计算中的"极化量子"问题。
【DFT Tip 2】Hubbard U 的双金属设置：Mn 和 V 的 U 值校准
超晶格中 Mn 和 V 都需要 U 值。Mn 3d 的 U=3-4 eV 是标准选择，但 V 3d 的 U 需要特别注意：V⁵⁺(d⁰) 的 U 值应使 PbVO₃ 再现实验极性 P4mm 结构。
方法：对 V 的 U 从 2 到 5 eV 扫描，检查 PbVO₃ 的声子不稳定模（Γ₅⁻）和极性畸变幅度。U 太小→V 金属化→无极性畸变，U 太大→V 过局域→假相。校准标准：U=3 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 μB，Mn⁴⁺(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:l=2（d 轨道），② 外窗能量范围 [EF-2 eV, EF+2 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 的 DFT+U+V 需要指定原子间 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↑ → 中心对称（P=0）→ P↓。关键参数：① 翻转能垒 E_barrier——决定翻转电压和保持时间，② 翻转路径中是否经过中间相（如非极性相），③ 翻转变分对磁性的影响（磁矩是否同时翻转）。
在 VASP 中：使用 VTST 工具包的 NEB 模块，5-7 个中间图像，Spring constant = -5。常见错误：NEB 图像数不足→翻转路径不光滑，高估能垒。
超晶格设计PMO:PVO=1: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 & Fu（2015）首次提出 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-Teller（JT）定理（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 空间群）弛豫，② 在 DFT+U 中，U 值的选择直接影响声子频率——U 太小（V 金属化）→虚频消失但极性相也被抑制，U 太大→虚频过大→高估不稳定性，③ 非谐效应——简谐声子计算在 T=0 K 下进行，而实验在室温下观察。
解决方案：① 凝聚不稳定声子模→弛豫到低对称相→检查低对称相的声子是否全部稳定，② 对 U 值做声子收敛测试（U=2,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 + SHG），K 点密度 ≥ 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 模拟磁相变与铁电相变
为什么值得算：本文所有计算在 T=0 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 曲率偶极子 · 非线性霍尔效应 · 四态存储
