KCdP中压力诱导的拓扑Dirac半金属相

 

PHYS. STAT. SOL. RRL 20, e202500398 (2026)

KCdP中压力诱导的拓扑Dirac半金属相

Pressure-Induced Topological Dirac Semimetallic Phase in KCdP

导读：Dirac半金属因其无质量Dirac费米子和拓扑保护的高迁移率表面态，在量子器件和自旋电子学领域具有重要应用前景。本文通过第一性原理计算，系统研究了KCdP在负三轴压力下的电子结构演化。核心发现：3%负压力下，无SOC时KCdP转变为三重简并Fermion半金属，含SOC后转变为Dirac半金属（vF=1.425x10^5 m/s）；5%压力下出现四个Dirac点。C6v对称性保护Dirac节点，P的s-和px+y-轨道之间的能带反转驱动拓扑相变。HSE06验证拓扑特征定性不变，声子计算确认全压力范围动力学稳定。

图片

一、前言背景

Dirac半金属：拓扑材料的重要分支

Dirac半金属（DSM）是一类具有四重简并零维节点和线性色散关系的拓扑半金属，其低能激发可描述为无质量Dirac费米子。与Weyl半金属（二重简并节点）和节线半金属（一维节点线/环）不同，DSM同时保留时间反演和反演对称性，展示巨抗磁性、量子磁阻、高载流子迁移率和Fermi弧表面态等独特性质。

DSM的形成有三种途径：(1) 非点式对称性在Brillouin区高对称点强制四重简并（如BiO2, BiZnSiO4）；(2) 拓扑相变临界点处的偶然Dirac节点（需要极端参数微调，实验实现困难）；(3) 能带反转诱导的偶然能带交叉DSM，由适当轴向旋转对称性保护（如Na3Bi, Cd3As2）。

本文研究的KCdP属于第三类：通过负三轴压力诱导能带反转，C6v对称性保护Dirac节点，在有限参数范围内保持拓扑稳健性。

压力诱导拓扑相变：从半导体到Dirac半金属

核心问题：外压如何驱动KCdP从普通半导体转变为拓扑Dirac半金属？SOC在这一转变中扮演什么角色？

方法体系：VASP（GGA-PBE + HSE06，电子结构） -> Phonopy（DFPT，动力学稳定性） -> Wannier90（紧束缚模型） -> WannierTools（表面态、Fermi面）。

关键发现：(1) 3%负三轴压力下，无SOC时出现三重简并Fermion（triple point semimetal），含SOC时转变为Dirac半金属；(2) 5%压力下出现双三重简并点（无SOC）->四个Dirac点（含SOC）；(3) HSE06确认拓扑特征定性不变，仅临界压力值从3%移至10%。

figure

KCdP压力诱导拓扑Dirac半金属研究流程。晶体结构（P63/mmc） -> VASP DFT（PBE+HSE06） -> 负三轴压力调控（0%->3%->5%->7%->10%） -> C6v对称性分析（无SOC: 三重简并Fermion; 有SOC: 四重简并Dirac点） -> 关键结果：Dirac锥（vF=1.425e5 m/s）、多Dirac点、表面态（001/010/100）、Fermi弧。Phonopy验证全压力范围动力学稳定性。

二、研究方法

DFT计算设置：VASP + HSE06 + Phonopy

VASP计算：PAW赝势，GGA-PBE泛函，ENCUT=500 eV，Gamma-centered 12x12x6 k点网格。结构弛豫使用共轭梯度法，力收敛标准0.005 eV/A。能带计算包含SOC。

HSE06验证：杂化泛函HSE06用于验证GGA-PBE结果的定性正确性。HSE06给出更大的带隙（常压896 meV vs PBE 220 meV），但拓扑相变的定性特征不变。临界负压力从PBE的3%移动到HSE06的10%。

Phonopy：3x3x2超胞，DFPT方法，在常压、3%和10%负三轴压力下计算声子谱，确认无虚频，验证动力学稳定性。

Wannier90 + WannierTools：紧束缚模型与拓扑分析

Wannier90：从VASP能带构建最大局域化Wannier函数（MLWFs），生成紧束缚模型。WannierTools用于计算拓扑表面态和Fermi面。

表面态计算：(001)面沿M-Gamma-K、(010)面沿Y-Gamma-Y、(100)面沿A-Gamma-A。三个表面均在Gamma点处展示Dirac点，Fermi面上出现Fermi弧特征。

Fermi速度：沿Dirac锥方向提取dE/dk斜率，得到vF = 1.425 x 10^5 m/s。该值与其他三维DSM相当，但小于石墨烯（约10^6 m/s），符合三维Dirac材料的预期。

formula

C6v对称性分析：C6z与sigma_v反对易保证沿kz轴的二重简并，非简并Gamma1导带与二重简并Gamma6价带交叉形成三重简并Fermion（无SOC）

formula

SOC下的Kramers定理：Theta^2=-1沿kz轴强制执行二重简并，C6v仅允许二维不可约表示，两套二重简并带交叉形成四重简并Dirac节点

formula

Fermi速度：从Dirac锥附近线性色散的斜率提取，vF=1.425x10^5 m/s

三、核心结果

figure

图 1：不同类型拓扑半金属的Fermi面特征示意图。(a) Dirac半金属在动量空间中具有离散的点状Fermi面。(b) Weyl半金属中Dirac点分裂为相反手性的Weyl节点对，表面Brillouin区出现开放的Fermi弧。(c) 节线半金属中能带交叉沿连续线或闭合环发生，产生扩展的环状Fermi面和鼓面表面态。

图片

图 2：KCdP的晶体结构和Brillouin区。(a) 晶体结构。(b) 不同角度视图。(c) 六方晶系的体相和表面Brillouin区[(100),(010),(001)]。KCdP属于P63/mmc空间群，a=4.44 A, c=10.18 A。

KCdP的晶体结构与电子性质

KCdP为六方晶系，空间群P63/mmc（No. 194）。优化晶格参数a=4.44 A, c=10.18 A。常压下为直接带隙半导体，GGA-PBE带隙220 meV，HSE06带隙896 meV。

常压能带（含SOC）：展示半导体特征，导带在Gamma点处呈锥状结构，暗示在合适的扰动下可能出现拓扑相变。无虚频声子谱确认常压动力学稳定性。

figure

图 3：KCdP的含SOC电子能带结构。(a) 常压KCdP--展示半导体特征，导带在Gamma点处呈锥状结构。(b) 3%负三轴压力下KCdP--导带和价带在Fermi能级处接触，形成拓扑Dirac点。

3%负三轴压力：从半导体到Dirac半金属

随着负三轴压力的施加，晶格常数增大（a=4.57 A, c=10.48 A），导带和价带逐渐靠近。在3%负压力下，能带在Fermi能级处交叠，实现拓扑相变。

无SOC：非简并Gamma1导带与二重简并Gamma6价带交叉，形成三重简并Fermion（triple point semimetal）。含SOC：Kramers定理使所有能带变为二重简并，三重简并点转变为四重简并Dirac点。

Dirac点在Fermi能级处出现，线性色散关系在所有方向成立，确认为Type-I DSM。3D能带色散图直观展示Dirac锥的线性E-k关系。

图片

图 4：KCdP的声子能带结构。(a) 常压KCdP--无虚频，确认动力学稳定。(b) 3%负三轴压力下KCdP--无虚频，确认压力诱导相在动力学上稳定。

声子稳定性：从常压到压力态

常压和3%负压力下的声子谱均无虚频，确认KCdP在所有考虑条件下动力学稳定。这一结果至关重要：许多理论预测的拓扑材料在实际中因结构不稳定而无法合成，声子验证是排除"虚假设"的关键步骤。

进一步地，10%负压力下的声子谱（HSE06临界压力点）同样无虚频，确认即使在大压力下KCdP仍保持动力学稳定。这表明KCdP的压力诱导DSM相在实验上可实现。

figure

图 5：KCdP的能带反转和Dirac半金属特征。(a) 含SOC的3D能带反转体相能带结构，展示两个对称性保护的Dirac点沿Gamma-A线出现。(b) 能带反转诱导的Dirac半金属：P的s-轨道和px+y-轨道之间的能带反转。

能带反转机制：s-轨道与px+y-轨道的角色

能带反转分析揭示了KCdP中DSM相的微观起源：P原子的s-轨道和px+y-轨道在负压力下发生能带反转。这与Na3Bi和Cd3As2中的能带反转机制类似（两者均为Type-I DSM）。

反转的物理图像：负压力增大晶格常数，减小价带和导带之间的晶体场劈裂，导致s-轨道（导带特征）和px+y-轨道（价带特征）的能量顺序反转。反转后，具有相反宇称的轨道在Fermi能级处交叉，形成拓扑非平庸的Dirac节点。

C6v对称性保护：Dirac节点沿Gamma-A线（C6轴）受到C6v小群对称性保护。两套二重简并带属于C6v的不同不可约表示，禁止杂化导致带隙打开，确保Dirac节点的稳健性。

figure

图 6：KCdP的表面态和Fermi面。(a) (001)面沿M-Gamma-K的表面态。(b) (010)面沿Y-Gamma-Y的表面态。(c) (100)面沿A-Gamma-A的表面态。(d-f) 对应表面的Fermi面（E=Ef），展示Fermi弧特征。

表面态与Fermi弧：拓扑Dirac半金属的确证

(001)、(010)和(100)三个表面的表面态计算均在Gamma点处展示Dirac点，Fermi面上出现Fermi弧特征。表面态的存在是拓扑Dirac半金属的标志性特征，来源于体-边对应原理。

Fermi面图（d-f）展示E=Ef处的二维谱函数A(k)，在Gamma点周围出现明显的亮环，对应于体Dirac锥在表面上的投影。Fermi弧连接Dirac点的投影，是ARPES实验可观测的拓扑特征。

Type-I DSM分类：Dirac锥无倾斜，Lorentz不变性保持，这与Na3Bi和Cd3As2同类。Type-I DSM在所有方向展示负磁阻，是实验验证的重要特征。

图片

图 7：GGA-PBE与HSE06能带对比。(a,b) GGA-PBE常压和3%负压力。(c,d) HSE06常压和10%负压力。拓扑相变的定性特征在两种泛函下一致，仅临界压力值不同。

HSE06验证：泛函依赖性分析

GGA-PBE vs HSE06的核心差异：(1) 带隙大小：常压下PBE=220 meV, HSE06=896 meV；(2) 临界压力：PBE=3%, HSE06=10%。但拓扑相变的定性特征完全一致--能带排序、轨道特征和Fermi能级附近的对称性均不变。

这一结果具有重要的方法论意义：虽然GGA-PBE系统性地低估带隙，但对于拓扑相变的定性预测（能带是否反转、Dirac点是否存在），PBE和HSE06给出相同结论。因此，GGA-PBE用于初步高通量筛选拓扑材料是可靠的，但定量预测（临界压力、带隙）需要HSE06校核。

作者的策略值得借鉴：用GGA-PBE进行主计算（成本低），用HSE06进行关键验证（确保定性正确），这在高通量拓扑材料筛选中是高效的。

figure

图 8：10%负三轴压力下KCdP的声子能带结构。无虚频，确认HSE06临界压力点（10%）的动力学稳定性。

DFT Tips

【DFT Tip 1】负压力DFT计算：概念与实现

负压力（Negative Pressure）在DFT中通过扩大晶格常数实现。本文使用"负三轴压力"（negative triaxial pressure），即沿三个晶轴方向均匀扩大晶格常数。

VASP实现：在POSCAR中手动修改晶格常数（如a=4.44->4.57 A），然后弛豫原子内部坐标（保持晶格常数固定）。注意：这是"施加应变"而非真正的"负压力"--DFT计算的应力张量分量可能不为零。

常见陷阱：负压力不同于静水压（hydrostatic pressure）。静水压下所有应力分量相等，而均匀晶格膨胀可能产生各向异性应力。对于六方晶系，a和c的膨胀比例可能不同，需要分别优化。

【DFT Tip 2】SOC对拓扑分类的定性影响：三重简并->Dirac点

本文展示了SOC如何改变拓扑相的类型：无SOC时三重简并Fermion，含SOC后转变为四重简并Dirac点。这一转变来源于SOC对能带简并度的重组（Kramers定理强制二重简并）。

DFT实践：必须同时计算无SOC和含SOC能带，对比分析。INCAR设置：LSORBIT=.TRUE.开启SOC。注意：SOC为无SOC的2-4倍，且收敛更困难。

常见错误：仅在无SOC下计算就声称发现了Dirac半金属。由于SOC打开能隙或改变简并度，无SOC的Dirac点可能在SOC下变为有能隙的拓扑绝缘体或Weyl半金属。

【DFT Tip 3】HSE06与GGA-PBE在拓扑材料中的选择

HSE06的带隙通常比GGA-PBE大0.5-1.5 eV，对于半导体-半金属转变，这将显著改变临界压力/应变/掺杂浓度。但拓扑相变的定性特征（能带是否反转）通常不变。

建议策略：(1) 用GGA-PBE进行大范围参数扫描（压力、应变、掺杂），确定拓扑相变的大致位置；(2) 在关键参数点上用HSE06验证；(3) 如果HSE06确认定性特征不变，则GGA-PBE的相图可作为定性参考。

注意：HSE06为PBE的10-20倍，k点网格需要对应减小。对于金属体系，HSE06的收敛比PBE更困难。

【DFT Tip 4】声子动力学稳定性：拓扑材料DFT的"必修课"

许多理论预测的拓扑材料在实空间结构不稳定（虚频），因此声子计算是拓扑材料DFT论文的标准要求。本文在常压、3%和10%三个压力点均计算了声子谱，确认无虚频。

DFT设置：Phonopy + VASP（DFPT或有限位移法），超胞大小至少3x3x2（如本文），力收敛标准<0.001 eV/A。对于极性材料，还需考虑LO-TO劈裂和非解析项修正。

常见陷阱：超胞不够大导致虚假虚频；有限位移法中位移幅度过大或过小；DFPT中未包含极性修正。建议：至少用两种方法（如DFPT和有限位移法）交叉验证声子谱。

【DFT Tip 5】Wannier90拟合：从VASP能带到紧束缚模型的关键步骤

KCdP的Wannier化涉及K-s、Cd-d、P-p轨道。投影选择：(1) K: s轨道；(2) Cd: d轨道（5个）；(3) P: s和p轨道（4个）。共约10-20个Wannier轨道。

关键步骤：(1) 解缠窗口（disentanglement）必须覆盖Fermi能级附近的所有目标能带；(2) 冻结窗口（frozen window）应包含Dirac点附近的能带；(3) 检查Wannier内插能带与DFT能带的一致性；(4) 确认Wannier函数的对称性（s轨道应为球对称，p轨道应具有方向性）。

常见陷阱：投影轨道选择不当导致Wannier拟合失败或虚假的拓扑性质。建议：先用VASPKIT提取轨道投影，确认目标能带的主要轨道成分，再选择投影。

【DFT Tip 6】Fermi速度的计算：从Dirac锥斜率到实验可测量

Fermi速度vF = (1/hbar) * |dE/dk|。在DFT中，从能带数据沿Dirac锥方向线性拟合提取斜率。对于KCdP，vF=1.425x10^5 m/s。

注意：(1) vF是方向依赖的，Dirac锥在不同方向上的斜率可能不同（各向异性Dirac锥）；(2) vF的DFT计算值受泛函选择影响--PBE可能低估vF；(3) vF是实验可测量（ARPES、磁阻振荡、量子振荡），DFT与实验的直接比较是验证Dirac半金属分类的重要方式。

常见陷阱：仅取单点斜率而非线性拟合区域，导致vF误差大。建议：在Dirac点附近|E-E_D|<0.1 eV范围内进行线性拟合。

【DFT Tip 7】表面态计算的数值收敛性

WannierTools通过迭代Green函数方法计算半无限体系的表面态。关键参数：(1) 表面终止面的选择（本文计算了001、010、100三个面）；(2) k点密度（表面Brillouin区至少100x100）；(3) 能量展宽（约0.005-0.01 eV）。

表面态的计算对Wannier函数的质量敏感。如果Wannier拟合质量差（spread大），表面态可能出现虚假的能隙或表面态。建议：在计算表面态之前，先确认Wannier内插能带与DFT能带的一致性。

注意：DFT计算的是单粒子表面态，不包含多体效应（如表面重构、电荷转移、自旋-轨道极化）。实验ARPES可能观察到与DFT表面态不同的特征。

【DFT Tip 8】拓扑材料DFT论文的"完整证据链"

本文为拓扑材料DFT论文提供了一个标准模板：(1) 晶体结构和电子结构（常压）；(2) 声子稳定性验证；(3) 参数调控下的能带演化（压力/应变）；(4) 对称性分析（群论解释能带简并和交叉）；(5) 泛函验证（HSE06 vs GGA-PBE）；(6) 表面态和Fermi面计算；(7) Fermi速度等可实验观测量。

缺失的环节（本文未包含但值得补充）：(1) Z2拓扑不变量或陈数计算；(2) 非平庸拓扑的输运特征（如反常Hall效应）；(3) 与已知实验数据的比较（KCdP尚无实验数据，但可比较同类材料）。

建议：DFT拓扑材料论文至少包含上述(1)-(6)项。(7)和额外的拓扑不变量/输运计算是加分项。

知识扩展

【知识扩展 1】拓扑半金属的分类：Dirac、Weyl与节线半金属

【理论解释】拓扑半金属根据能带交叉的维度分类：(1) Dirac半金属--四重简并零维节点，同时具有TRS和IS，两个相反手性Weyl点重叠；(2) Weyl半金属--二重简并零维节点，TRS或IS破缺，手性相反的Weyl点分离；(3) 节线半金属--一维节线或环，需额外对称性保护。

【Dirac半金属的三种类型】Type-I：无倾斜Dirac锥，Lorentz不变性保持，负磁阻（如Na3Bi, Cd3As2, KCdP）；Type-II：强倾斜Dirac锥，Lorentz不变性破缺，电子和空穴口袋共存（如PtSe2, VAl3）；Type-III（临界DSM）：倾斜但具有线状Fermi面，拓扑不变量N2=1（如Zn2In2S5）。

【经典参考】Young et al., PRL 108, 140405 (2012)--3D DSM理论预测；Liu et al., Science 343, 864 (2014)--Na3Bi实验发现；Liu et al., Nat. Mater. 13, 677 (2014)--Cd3As2实验发现；Armitage et al., RMP 90, 015001 (2018)--拓扑半金属综述。

【迁移能力】压力诱导拓扑相变的方法可推广到其他窄带隙半导体，特别是具有六方/四方对称性和重元素的体系。

【知识扩展 2】压力诱导拓扑相变：从等静压到应变工程

【理论解释】压力通过改变原子间距、轨道重叠和晶体场劈裂，驱动能带反转和拓扑相变。当价带和导带在压力下交换能量顺序时，能带反转发生，通常伴随拓扑不变量的改变。

【方法比较】静水压：通过DAC实验实现，各向同性压缩，理论计算简单但实验条件苛刻。单轴应变：沿特定方向施加，可选择性调控特定能带，实验相对容易。负压力/膨胀：晶格常数扩大，在DFT中易于实现但实验上难以直接施加，通常通过晶格失配外延或化学掺杂间接实现。

【实验验证】SnTe在应变下被ARPES确认为拓扑晶体绝缘体（Fragkos et al., PRM 3, 104201）。LaAs在压力下被预测为拓扑半金属（Khalid et al., PRB 98, 220102）。NaBiBa在压力下被预测为拓扑绝缘体（Sun et al., PRB 93, 205303）。

【经典参考】Bansil et al., RMP 88, 021004 (2016)--拓扑材料综述；2016年Nobel Prize专题--拓扑相变理论。

科研经验

【科研经验 1】拓扑材料DFT论文中的"对称性分析"如何写

问题：许多DFT论文的对称性分析部分只是引用群论表格，缺乏物理洞察。如何让对称性分析既有理论深度又易于理解？

原因：群论是抽象的，但拓扑物理本质上是对称性物理。如果对称性分析不能解释"为什么能带交叉是受保护的"和"为什么SOC改变简并度"，审稿人会觉得分析流于形式。

解决方案：本文的对称性分析值得借鉴：(1) 从具体对称操作出发（C6z、sigma_v），而非抽象群论；(2) 展示反对易关系如何强制简并（C6z与sigma_v反对易->双二重简并）；(3) 展示SOC如何改变简并度（Kramers定理Theat^2=-1）；(4) 解释为什么两套二重简并带交叉形成Dirac点（不同IR间无杂化）。

建议：对称性分析的核心价值在于解释"为什么"而非"是什么"。在写对称性分析时，始终问自己：读者能否从这段分析中理解拓扑保护的微观机制？

【科研经验 2】压力诱导拓扑相变：从DFT预测到实验验证的鸿沟

问题：DFT预测的压力诱导拓扑相变，在实验中如何验证？负压力在实验中如何实现？

原因：(1) 负压力（晶格膨胀）在实验中无法直接施加，需要间接方法；(2) ARPES需要高质量单晶表面，压力下测量困难；(3) 输运测量（如SdH振荡、负磁阻）需要高迁移率样品。

解决方案：(1) 替代实现：晶格失配外延（如SiC衬底上外延生长）、化学掺杂（大原子替代小原子）、插层膨胀；(2) 实验验证：首先确认材料可合成（本文的声子验证是第一步），然后测量常压下的输运和ARPES，再通过施加外压逼近临界点；(3) 间接证据：压力下的电阻率异常、Hall效应变化、量子振荡。

建议：DFT论文中应明确讨论实验实现的可行路径，而不仅仅停留在理论预测。审稿人通常会问"这个预测如何被实验验证"。

 

Gupta, Singh, Sen, Patra & Singh, VNIT Nagpur | phys. stat. sol. RRL 20, e202500398 (2026) | Dirac半金属 压力诱导拓扑相变 KCdP VASP HSE06