企业官网建设流程全解析

FeSe超薄膜在CaF2衬底上的电子结构DFT研究

MATER. RES. EXPRESS 13, 126002 (2026)

FeSe超薄膜在CaF2衬底上的电子结构DFT研究

DFT-Based Analysis of Substrate Effects on FeSe Ultrathin Films on CaF2

导读 导读：FeSe是最简单的铁基超导体，其超导性质强烈依赖于衬底选择。与FeSe/SrTiO3（Tc~65-100 K）不同，FeSe/CaF2的Tc仅约15 K，但CaF2的化学惰性和宽禁带使其成为研究FeSe本征超导性的理想平台。本文通过Quantum ESPRESSO DFT计算，揭示了CaF2极性表面重构产生的界面偶极子如何在标准NM DFT中消除Gamma空穴口袋--无需Hubbard U或电子掺杂，与FeSe/STO需要强关联修正形成鲜明对比。随FeSe层数增加（1L->8L），空穴口袋逐渐恢复，暗示从incipient-band配对到传统多带s+/-配对的维度跨越。

一、前言背景

铁基超导体FeSe：从体材料到单层极限

FeSe（硒化铁）是铁基超导体家族中最简单的成员，具有PbO型四方晶体结构。体材料FeSe的超导转变温度Tc仅约8 K，但通过化学掺杂、外加压力或外延应变可将Tc显著提升至约37 K。其超导配对机制被认为是非传统的，由反铁磁自旋涨落介导，配对对称性为s+/-波。

然而，当FeSe被制备为超薄膜时，其超导性质表现出强烈的衬底依赖性。在SrTiO3（STO）衬底上，单层FeSe的Tc可达65-100 K以上，这一发现引发了铁基超导领域的研究热潮。而在双层石墨烯上，超导性被抑制，单层极限下Tc低于2.2 K。

CaF2（氟化钙）是近年来新兴的FeSe衬底材料。与STO不同，CaF2是宽禁带离子绝缘体，具有化学惰性和极性表面，与FeSe的晶格失配极小。实验上FeSe/CaF2的Tc约15 K，提供了研究FeSe本征超导性的理想平台，避免了STO体系中极端的界面声子增强效应。

核心科学问题：衬底如何改变FeSe的费米面拓扑？

核心问题：CaF2衬底对FeSe超薄膜电子结构的影响机制是什么？与STO衬底有何本质不同？为什么FeSe/CaF2的Tc（~15 K）远低于FeSe/STO（~65 K）？

方法体系：Quantum ESPRESSO（DFT, PBE+US赝势, ENCUT=49 Ry wavefunction/442 Ry charge, 3x3x1 k点, 9层CaF2 slab + 1-8层FeSe）-> Bader电荷分析（界面电荷转移）-> 非磁/棋盘AFM磁构型对比。

关键发现：(1) CaF2衬底通过表面重构产生的界面偶极子，在标准NM DFT中即可消除Gamma点空穴口袋--不需要Hubbard U或电子掺杂；(2) 与FeSe/STO需要强关联修正才能消除空穴口袋形成鲜明对比；(3) 随FeSe层数增加（2L->8L），空穴口袋逐渐恢复，7L/8L时重新穿过费米面。

FeSe/CaF2衬底效应与电子结构研究流程。CaF2(100)极性表面重构（50%F转移）-> FeSe/CaF2异质结构建（a=3.86A，晶格失配极小）-> Quantum ESPRESSO DFT计算（PBE+US, NM/AFM）-> DOS+PDOS+能带+Bader电荷分析。核心发现：纯静电效应在NM DFT中即可消除Gamma空穴口袋。层数依赖：1L纯电子口袋 -> 7L-8L空穴口袋恢复。

二、研究方法

Quantum ESPRESSO计算设置：异质结建模

计算平台：Quantum ESPRESSO（QE），而非VASP。这是本文的一个重要方法学特征。QE使用平面波基组+赝势方法，与VASP功能等价但参数设置不同。注意：QE中截断能以Ry为单位（1 Ry = 13.606 eV），49 Ry ~ 667 eV，电荷密度截断442 Ry ~ 6014 eV。

CaF2衬底建模：9层CaF2(100) slab，晶格常数a=5.46 A。CaF2(100)表面是本征极性表面（Tasker type III），由交替的Ca2+和F-平面组成，产生宏观偶极子。通过将顶面50%的F原子转移到底面进行表面重构，补偿表面偶极子，确保静电稳定性。

FeSe/CaF2异质结：将1-8层FeSe置于CaF2 slab顶部。面内晶格常数a=3.86 A（a_CaF2/sqrt(2)），与FeSe体材料晶格常数失配极小。真空层厚度至少等于异质结总厚度。所有原子在结构优化中完全弛豫。

NM vs 棋盘AFM：磁构型的影响

非磁（NM）计算：关闭自旋极化，直接计算电子结构。本文的关键发现在于：NM DFT中Gamma点空穴口袋已经消失，说明纯静电效应（而非磁关联效应）是主要驱动力。

棋盘AFM（checkerboard AFM）：Fe原子在面内形成最近邻反平行排列。这是FeSe体系中常用的有序近似，用于模拟顺磁态的自旋涨落效应。AFM有序进一步抑制价带，引入自旋相关的能带劈裂。

与FeSe/STO对比：STO体系中，NM DFT仍然保留Gamma空穴口袋，需要Hubbard U或AFM+电子掺杂才能消除。CaF2体系中NM DFT即可消除空穴口袋，说明CaF2的界面静电效应更强。

晶格失配定义：delta为CaF2与FeSe晶格常数之差除以FeSe晶格常数，CaF2(100)的a/sqrt(2)=3.86A与FeSe的a匹配极好

Born-Huang力学稳定性判据：C11>0且C11>|C12|，确保2D材料弹性稳定

三、核心结果

图 1：(a) 单层FeSe/CaF2异质结的晶体结构（1L FeSe置于9层CaF2衬底上）。(b) 棋盘反铁磁（AFM）自旋构型。(c) FeSe/CaF2体系的二维布里渊区和高对称点。

图 2：非磁单层FeSe/CaF2的态密度（DOS）。(a) 原子投影DOS：CaF2衬底在费米面附近无电子态贡献，Fe态主导。(b) Fe原子轨道分辨DOS：Fe-dxz, dyz, dxy轨道在费米面附近占主导，与铁基超导体的已知特征一致。

DOS分析：Fe-d轨道与衬底惰性

CaF2衬底在费米面附近完全没有电子态，这是其作为"理想衬底"的关键优势。与STO不同，STO的O-2p轨道可能与Fe-3d轨道杂化，CaF2的F-2p轨道能量极低（~-20 eV以下），完全不会与FeSe的能带混合。

Fe-d轨道中dxz、dyz和dxy在费米面附近占主导，与铁基超导体的标准图像一致。这些轨道形成的费米面口袋（电子口袋在M点，空穴口袋在Gamma点）的嵌套决定了自旋涨落和超导配对强度。

轨道分辨PDOS是理解铁基超导体的第一窗口。在分析任何铁基超导体时，PDOS是必须首先检查的量，它可以揭示：(1) 哪些轨道贡献费米面态密度；(2) 是否存在轨道选择的Mott物理；(3) 衬底是否引入了不需要的界面态。

图 3：单层FeSe的电子能带结构对比。(a) 1L FeSe/CaF2 NM态：Gamma点空穴口袋消失，价带顶低于费米面。(b) 1L FeSe/CaF2 AFM态：价带进一步被抑制，出现自旋相关劈裂。(c) Free-standing 1L FeSe（使用与CaF2相同的面内晶格常数）：展示清晰的Gamma空穴口袋，作为参考。

Gamma空穴口袋消失：纯静电效应 vs 强关联效应

这是本文最核心的发现。在NM DFT计算中，FeSe/CaF2的Gamma点价带顶已经低于费米面，空穴口袋消失。这与FeSe/STO形成鲜明对比：STO体系中NM DFT仍然保留空穴口袋，需要Hubbard U修正或AFM有序+电子掺杂才能消除。

物理机制：CaF2(100)极性表面重构后，50%F原子转移在界面处产生残余电荷不对称性，形成界面偶极子。该偶极子修改了宏观静电势，导致FeSe的能带整体发生刚性下移（downward band shift），Gamma点价带移出费米面。

Bader电荷分析确证：界面净电荷转移仅约0.01e/unit cell，极其微小。这说明空穴口袋的消失不是由电荷转移驱动的，而是由界面偶极子产生的静电势偏移驱动的。这是极性界面物理学的一个优美案例。

图 4：FeSe/CaF2界面的电荷密度重分布。(a) NM态。(b) AFM态。黄色等值面表示电荷积累，青色表示电荷耗尽。Bader电荷分析显示净电荷转移仅约0.009-0.011e/unit cell。AFM态中Fe原子周围出现更显著的电荷重分布。

界面电荷重分布：偶极子形成的微观证据

电荷密度差图展示了界面偶极子的微观起源：Se原子下方出现电荷积累（黄色），相邻区域出现电荷耗尽（青色），形成空间电荷分离。这正是界面偶极子的微观特征。

AFM态中Fe原子周围的电荷重分布比NM态更显著，反映了磁有序对局域电子环境的影响。但净电荷转移仍然极小，说明磁有序主要影响电荷的空间分布而非总转移量。

这种电荷重分布模式与STO界面完全不同。STO界面通常涉及O空位、Ti3+形成等复杂化学过程，而CaF2界面是纯粹的物理静电效应。这使CaF2成为研究FeSe本征超导性的更干净平台。

图 5：不同层数FeSe/CaF2（NM态）的电子能带结构。(a) 2L: 空穴口袋仍缺失。(b) 3L: 价带开始上移。(c) 4L: 价带接近费米面但未穿过。(d) 6L: 价带非常接近费米面。(e) 7L: 空穴口袋重新出现。(f) 8L: 空穴口袋清晰，体材料特征恢复。

层数依赖的维度跨越：从incipient band到多带s+/-

1L-4L：衬底效应主导。Gamma空穴口袋缺失，费米面仅由M点电子口袋组成。这一费米面拓扑结构暗示超导配对可能由"incipient band"机制驱动：配对发生在电子口袋之间，而非传统的电子-空穴口袋嵌套。

6L-8L：衬底效应衰减。7L和8L时Gamma空穴口袋重新穿过费米面，标志着从准二维向体材料行为的维度跨越。此时传统多带s+/-配对机制可能恢复。

超导配对机制的推测：CaF2上Tc~15 K的温和增强，可能正是因为缺乏STO体系中的强界面声子耦合和关联驱动重整化。超薄FeSe/CaF2的电子-口袋-only配对效率低于电子-空穴口袋嵌套，导致Tc低于STO体系。

DFT Tips

【DFT Tip 1】Quantum ESPRESSO vs VASP：截断能单位与设置差异

QE使用Ry作为能量单位，VASP使用eV。1 Ry = 13.606 eV。本文WFC截断49 Ry ~ 667 eV，电荷密度截断442 Ry ~ 6014 eV。QE中波函数和电荷密度使用不同截断能（ecutwfc和ecutrho），而VASP中ENCUT同时控制两者（电荷密度默认为4xENCUT）。

QE赝势选择：本文使用US（超软）赝势，而非PAW。US赝势的截断能通常低于PAW，但需要更高的电荷密度截断。对于FeSe体系，Fe的3d电子需要特别注意赝势的精度，建议使用建议的截断能或进行收敛测试。

常见陷阱：从VASP文献直接迁移参数到QE时，务必注意单位转换。VASP的ENCUT=500 eV对应QE的ecutwfc~37 Ry。使用QE的SSSP效率赝势库（PBE efficiency）可降低截断能要求。

【DFT Tip 2】极性表面建模：Tasker分类与偶极子补偿

CaF2(100)是Tasker type III极性表面，由交替的带电平面（Ca2+和F-）组成，在未重构的slab中会产生发散的宏观偶极子。DFT计算中必须进行表面重构以消除偶极子。

常见补偿方案：(1) 本文的50%F转移（顶面->底面）；(2) 表面羟基化/氢钝化（对于氧化物表面更常见）；(3) 分数氢原子（伪氢钝化，用于半导体表面）。

验证偶极子补偿是否成功：(1) 检查slab中心区域的静电势是否平坦（bulk-like）；(2) 计算slab两侧真空能级是否对齐；(3) 检查总能量是否随slab厚度线性收敛。如果静电势在slab内部仍然倾斜，说明偶极子未完全补偿。

【DFT Tip 3】异质结计算中的晶格失配与应变处理

FeSe/CaF2的晶格失配极小（a_CaF2/sqrt(2)=3.86A vs FeSe a~3.77A），~2.4%。这种小失配使得可以忽略应变效应，直接使用衬底晶格常数约束FeSe。

当晶格失配>5%时，必须考虑以下选项：(1) 构建超胞匹配（如4x4 FeSe on 3x3 substrate）；(2) 使用应变FeSe（计算不同晶格常数下的能量，取最小值）；(3) 使用Moire超胞（计算成本急剧增加）。

本文的简化处理（将free-standing FeSe也用CaF2的a=3.86A计算）是合理的，因为失配很小，且实验上FeSe薄膜确实被衬底应变。但需要注意：free-standing参考态也使用相同晶格常数，排除了应变作为混淆变量的可能性。

【DFT Tip 4】棋盘AFM作为顺磁态的有序近似

铁基超导体中，棋盘AFM（checkerboard）是常用的长程磁有序近似，用于模拟顺磁态中动态自旋涨落的平均效应。这种方法可以捕获自旋涨落对能带结构的主要修正，但只是近似。

棋盘AFM的局限性：(1) 它引入了人工的长程磁有序，可能高估某些能带劈裂；(2) 它无法捕获自旋涨落的动态（频率依赖）特性；(3) 对于接近磁量子临界点的体系，有序近似可能定性不正确。

更准确的方法：DFT+DMFT（动力学平均场）可以同时处理局域关联和巡游磁性，但。对于FeSe/CaF2这类以定性分析为目标的工作，棋盘AFM是合理的折中。

【DFT Tip 5】Bader电荷分析在异质结界面的正确使用

本文的Bader分析揭示净电荷转移仅~0.01e/unit cell，是核心证据之一。但需要注意：Bader电荷是拓扑划分，绝对值有方法依赖性，只有相对变化趋势有意义。

QE中Bader分析流程：(1) 自洽计算生成电荷密度文件；(2) 使用pp.x后处理生成cube格式电荷密度；(3) 使用Bader分析程序（如Henkelman组的代码）进行拓扑划分。

常见陷阱：对于异质结，如果只做整体Bader分析，无法区分界面区域的电荷重分布。更精细的方法是将电荷密度差（delta_rho = rho_FeSe/CaF2 - rho_FeSe - rho_CaF2）可视化，直观展示界面的电荷积累/耗尽区域。

【DFT Tip 6】铁基超导体能带分析的关键高对称k点

铁基超导体的费米面拓扑分析集中在两个高对称点：(1) Gamma点（布里渊区中心），通常出现空穴口袋（dxz/dyz轨道）；(2) M点（布里渊区角），通常出现电子口袋（dxy/dxz/dyz轨道）。

嵌套（nesting）分析：Gamma空穴口袋和M电子口袋的面积/形状匹配程度决定了自旋涨落的强度。嵌套越完美，Tc越高。在FeSe/CaF2中，1L时Gamma口袋消失，嵌套被破坏，暗示配对机制从传统的电子-空穴嵌套转变为incipient-band配对。

QE中能带计算路径：Gamma(0,0,0)->X(0.5,0,0)->M(0.5,0.5,0)->Gamma(0,0,0)。对于2D正方晶格，这是标准的高对称路径。

【DFT Tip 7】Hubbard U在铁基超导体中的争议

铁基超导体中是否需要Hubbard U是一个长期争议的问题。支持方：U可以纠正Fe-3d电子的自相互作用误差，改善磁矩和能带宽度。反对方：Fe-3d电子在铁基超导体中相对巡游，U可能过分局域化，破坏费米面拓扑。

本文的价值在于：在NM DFT中已经观察到Gamma空穴口袋消失，不依赖U修正。这避免了U值选择的主观性，使结论更为稳健。

建议：对于铁基超导体，应先做NM DFT获取基准结果，观察关键物理特征是否已经出现。如果NM DFT中已经出现目标特征，则不依赖U的结论更可信。如果必须使用U，建议在U+/-1 eV范围内做敏感性测试。

【DFT Tip 8】slab计算中的真空层收敛性

异质结slab计算中，真空层厚度必须足够大以消除周期性镜像的虚假相互作用。本文要求真空层至少等于异质结总厚度，这对于9层CaF2+1-8层FeSe的体系，真空层~20-30 A。

收敛检查：(1) 比较不同真空层厚度下的总能量和功函数；(2) 检查静电势在真空区域是否平坦；(3) 确保费米面附近的能带不因真空层厚度变化而显著移动。

对于极性表面slab，偶极子校正（QE中的assume_isolated="esm"或dipole correction）可能需要。但本文通过表面重构已经补偿了偶极子，额外的偶极子校正可能不是必需的。

知识扩展

【知识扩展 1】铁基超导体的配对机制：s+/-波与自旋涨落

【理论解释】铁基超导体的超导配对被认为由反铁磁自旋涨落介导。在体材料FeSe中，Gamma空穴口袋和M电子口袋之间的嵌套产生强自旋涨落，导致s+/-波配对对称性：Gamma口袋和M口袋上的超导能隙符号相反。

【incipient band配对】当空穴口袋下沉到费米面以下时（如1L FeSe/CaF2），它不再参与费米面态的配对。但Bang等人提出，即使能带不穿过费米面（incipient band），只要其能量足够接近费米面，仍可通过虚过程参与配对。

【方法比较】s+/-波（电子-空穴嵌套）vs incipient band配对（电子-电子口袋）：前者配对强度更强但需要Gamma口袋存在，后者配对强度较弱但更鲁棒。FeSe/CaF2的温和Tc~15 K可能正是incipient band配对效率较低的体现。

【经典参考】Graser et al., NJP 11, 025016 (2009) -- 铁基超导体配对对称性；Bang, NJP 18, 113054 (2016) -- incipient band配对；Lee et al., Nature 515, 245 (2014) -- FeSe/STO界面声子增强。

【迁移能力】incipient band配对的概念可推广至其他超薄超导体、界面超导体系以及掺杂Mott绝缘体。

【知识扩展 2】极性界面物理学：从LaAlO3/SrTiO3到FeSe/CaF2

【理论解释】极性界面（polar interface）是指由带电原子平面组成的两种材料之间的界面。当极性材料与非极性材料接触时，界面处会出现"极性灾难"（polar catastrophe）：静电势发散，需要通过电子重构或原子重构来补偿。

【LaAlO3/SrTiO3】最著名的极性界面体系。LaAlO3由(LaO)+和(AlO2)-交替平面组成，SrTiO3由中性(SrO)0和(TiO2)0平面组成。界面处出现二维电子气（2DEG），具有超导性和磁性。

【FeSe/CaF2与LAO/STO的对比】LAO/STO中电子重构是主要补偿机制（产生2DEG），而FeSe/CaF2中原子重构是主要补偿机制（50%F转移消除偶极子）。CaF2的宽禁带限制了电子重构，因此界面电荷转移极小。

【经典参考】Ohtomo & Hwang, Nature 427, 423 (2004) -- LAO/STO界面2DEG；Noguera, JPCM 12, R367 (2000) -- 极性氧化物表面综述；Tasker, J. Phys. C 12, 4977 (1979) -- 极性表面分类。

【迁移能力】极性界面物理学适用于所有离子型异质结体系，包括钙钛矿氧化物界面、氟化物/氧化物界面、氮化物/半导体界面等。

科研经验

【科研经验 1】铁基超导体DFT计算中的"Gamma空穴口袋"争议

问题：不同文献报道的FeSe单层Gamma空穴口袋存在与否的结论相互矛盾，应该如何理解？

原因：(1) 不同衬底（STO vs CaF2 vs free-standing）的界面效应不同；(2) 是否使用Hubbard U；(3) 晶格常数（应变）的影响；(4) 磁构型的选择（NM vs AFM vs SDW）。这些因素都会显著改变Gamma点价带相对于费米面的位置。

解决方案：(1) 系统比较不同衬底、不同U值、不同磁构型下的能带结构，找出稳健的物理趋势；(2) 明确说明结论的适用范围（"在NM DFT+CaF2衬底下，Gamma空穴口袋消失"）；(3) 利用Bader电荷和静电势分析揭示背后的物理机制，使结论超越参数选择。

建议：不要仅凭一组参数的结果就下结论。铁基超导体的能带结构对计算参数极其敏感，需要在参数空间中做系统扫描才能得出可靠结论。本文通过对比NM/AFM、不同衬底、不同层数，建立了一个相对完整的物理图像。

【科研经验 2】Quantum ESPRESSO与VASP的互操作性

问题：课题组使用VASP，但本文使用QE，如何理解和复现QE的结果？

原因：QE和VASP虽然都是DFT平面波代码，但赝势库、截断能单位、输入文件格式完全不同。直接迁移参数可能导致错误结果。

解决方案：(1) 截断能转换：QE的49 Ry ~ VASP的667 eV，如果自己的VASP计算使用500 eV，可能需要稍微提高；(2) 赝势选择：QE的US赝势与VASP的PAW赝势不同，建议先用简单体系（如体FeSe）做基准对比；(3) k点设置：3x3x1对于slab计算偏小，复现时建议使用更密的k点（如6x6x1）做收敛测试。

建议：不要盲目信任文献中的参数。每个课题组应建立自己的收敛测试标准。对于铁基超导体这类参数敏感的体系，复现实验时建议从收敛测试开始。

【科研经验 3】DFT预测超导配对机制的局限性

问题：DFT本身不包含超导性（超导是BEC/BCS凝聚现象），如何从DFT能带结构推断超导配对机制？

原因：DFT计算的是正常态的电子结构。超导配对机制需要从正常态性质（费米面拓扑、态密度、自旋涨落、EPC）间接推断。本文的"incipient band配对"和"s+/-配对"都是基于DFT能带结构的推测，不是DFT的直接输出。

解决方案：(1) 明确说明"配对机制讨论仅为推测，非定量预测"；(2) 如果需要更直接地研究超导性，需要结合以下方法：Migdal-Eliashberg理论（EPC超导）、RPA自旋涨落计算（非传统超导）、BdG方程（实空间超导态）；(3) 理想情况下，DFT提供正常态参数，后续模型给出超导预测。

建议：在DFT论文中讨论超导机制时，一定要区分"DFT的直接结果"和"基于DFT结果的推测"。本文在第3.3节明确标注"speculative rather than quantitative prediction"，是值得学习的严谨做法。

如果是我，我还会继续算

【继续算 1】SOC效应：自旋轨道耦合对能带的影响

为什么值得算：本文未包含SOC，但已知FeSe中SOC能级在20-40 meV量级。SOC可以引起Gamma和M点附近的能带劈裂，虽不改变费米面拓扑，但可能影响自旋涨落和配对对称性。

能回答的问题：SOC是否改变Gamma空穴口袋消失的结论？SOC引起的能带劈裂是否影响incipient band的有效能量距离？

适合体系：所有含重元素的铁基超导体。输入：QE中设置lspinorb=.true.和noncolin=.true.

【继续算 2】EPC计算：界面声子模式的作用

为什么值得算：FeSe/STO的高Tc被归因于界面声子模式（F-K模式）的强前向散射EPC。FeSe/CaF2中是否存在类似的界面声子模式？其EPC强度如何？

能回答的问题：CaF2的声子模式是否与FeSe电子耦合？界面EPC的强度与STO相比如何？这能否解释Tc差异？

适合体系：所有异质结超导体系。输入：QE的PHonon模块+EPW（需大超胞+密集k/q网格）。

【继续算 3】Hubbard U的敏感性分析

为什么值得算：本文使用NM DFT（无U）获得主要结论，但铁基超导体中U=2-4 eV是常见选择。系统研究U对Gamma空穴口袋的影响可以验证结论的鲁棒性。

能回答的问题：多大幅度的U值会使Gamma空穴口袋重新出现？CaF2衬底效应是否在U存在时仍然主导？

适合体系：所有铁基超导体。输入：QE中设置Hubbard U（lda_plus_u=.true.）。

【继续算 4】费米面嵌套函数与自旋磁化率

为什么值得算：超导配对强度直接与自旋磁化率chi(q)相关。计算Lindhard响应函数可以量化费米面嵌套的强度，给出超导Tc的定性趋势。

能回答的问题：不同层数FeSe/CaF2的自旋磁化率峰值在什么q矢量？嵌套强度如何随层数演化？incipient band配对对应的自旋涨落能谱是什么？

适合体系：所有超导体。输入：Wannier90紧束缚模型 + Lindhard函数计算。

【继续算 5】载流子掺杂效应：静电栅压调控费米面

为什么值得算：通过静电栅压注入电子或空穴，可以移动费米能级，调控Gamma空穴口袋相对于费米面的位置。这是实验上调控FeSe超导性的标准手段。

能回答的问题：电子掺杂能否使Gamma空穴口袋远离费米面，增强incipient band配对？空穴掺杂能否恢复Gamma口袋，切换回s+/-配对？

适合体系：所有超薄超导体。输入：QE中通过tot_charge标签模拟掺杂，或使用刚性能带模型。

【继续算 6】替代衬底比较：MgO、BaF2、LiF

为什么值得算：系统比较不同离子绝缘体衬底（MgO、BaF2、LiF）对FeSe电子结构的影响，建立衬底-界面偶极子-费米面拓扑的定量关系。

能回答的问题：界面偶极子强度与衬底的Madleung势有何关系？是否存在最优衬底使Tc最大化？

适合体系：所有异质结体系。输入：为每种衬底构建slab模型+异质结（需多个独立计算）。

Jacosalem, Domato, Ambolode | Mater. Res. Express 13, 126002 (2026) | FeSe CaF2衬底 Quantum ESPRESSO 界面偶极子