简介这份资料面向凝聚态物理与材料计算方向的学习者聚焦狄拉克半金属中狄拉克电导率的理论分析与数值计算。包内共2个文件含1个MATLAB脚本与1篇PDF文献压缩包约1.12MB。脚本可用于模拟电导率其中对虚部采用阶跃函数近似处理便于简化复数频率下的吸收与散射计算PDF则围绕三维狄拉克半金属薄膜的介电响应与新颖电磁模式展开涉及狄拉克锥、狄拉克点及边界效应等关键概念。已有710人学习下载适合希望从能带结构出发理解零质量电子输运、并借助脚本复现电导率曲线的读者参考也可作为薄膜电磁模式研究的入门材料。1. 狄拉克电导率到底在算什么从半金属到狄拉克半金属的一条主线第一次接触“狄拉克电导率”这个词很多人会以为它只是石墨烯的专属名词。其实只要体系里存在线性色散的低能激发也就是能量和动量近似成正比、能带在一点或一条线上简并那套用狄拉克方程描述的电导率公式就能派上用场。狄拉克半金属正是这样一类材料它的导带和价带在布里渊区里以四重简并的狄拉克点接触低能载流子表现得像无质量的相对论粒子。于是问题就变成——这种“无质量”的载流子在电场里到底怎么输运电导率随频率、温度、化学势怎么变。搞懂这条主线你才能判断一个材料值不值得做器件、该测哪个频段、拟合时该锁哪些参数。这篇笔记就按“概念立住 → 公式落地 → 参数怎么取 → 坑在哪”的顺序把狄拉克电导率从半金属一路讲到狄拉克半金属的实操路径。2. 从半金属到狄拉克半金属能带、态密度与电导率的物理起点2.1 半金属和普通金属的差别在哪普通金属的费米面穿过能带态密度在费米能级处是个有限值载流子有有效质量电导率的 Drude 形式大家都很熟。半金属不一样它的导带底和价带顶在能量上轻微重叠费米能级处的态密度趋近于零但不完全为零载流子浓度比普通金属低两三个数量级。这个“低载流子浓度 小费米面”的组合直接导致半金属的电导率对温度、杂质和缺陷特别敏感。你如果在实验里看到电阻率随温度先降后升、低温残余电阻比很高大概率就是半金属体系里载流子被冻析或者缺陷散射主导了。狄拉克半金属是半金属里更特殊的一支。它的导带和价带不是简单重叠而是在动量空间某些高对称点或高对称线上形成四重简并的狄拉克点。低能激发满足 E ±v_F |p|v_F 是费米速度p 是相对狄拉克点的动量。这个线性色散让态密度在狄拉克点附近按 |E|² 变化三维或 |E|二维和普通抛物带金属的 √E 完全不同。态密度的形状一变电导率对化学势、频率和温度的依赖关系就全变了。2.2 为什么狄拉克半金属的电导率值得单独算普通金属的光电导率在低频下基本是平的由 Drude 权重决定。狄拉克半金属不一样它的带间跃迁在低频就有贡献而且带内和带间两部分会随频率、温度、化学势互相竞争。更麻烦的是三维狄拉克半金属的带间电导率在低频极限下趋近于一个常数这个常数只和费米速度、普朗克常数有关和载流子浓度无关。这意味着你没法靠调掺杂把低频电导率压下去器件的暗电流、热噪声都会受这个常数限制。做红外探测器或者高频器件的人如果忽略这一点设计出来的响应度可能和预期差一个量级。从建模角度狄拉克半金属的电导率通常写成带内 Drude 项加带间跃迁项。带内项描述费米面附近的载流子被电场加速带间项描述电子从价带被光子激发到导带。在有限温度下带间项还要乘一个费米-狄拉克分布差把热激发和泡利阻塞都算进去。这套公式在石墨烯里已经验证得很充分搬到三维狄拉克半金属时主要变化是态密度从线性变成二次积分维度从二变三。2.3 一个最小可用的电导率表达式实操里最常用的形式是 Kubo 公式在清洁极限下的结果。对三维狄拉克半金属带内电导率可以写成import numpy as np def dirac_intraband_sigma(omega, mu, temp, vf, tau): 三维狄拉克半金属带内电导率简化 Drude 形式 omega: 角频率 (rad/s) mu: 化学势 (eV) temp: 温度 (K) vf: 费米速度 (m/s) tau: 弛豫时间 (s) hbar 1.0545718e-34 kB 8.617333e-5 # eV/K e 1.6021766e-19 # 有限温度下的载流子浓度近似 if mu 0: n (mu**3) / (3 * np.pi**2 * (hbar * vf)**3) else: n 0.0 # Drude 权重 D e**2 * vf * n / (hbar * 3) # 三维狄拉克的 Drude 权重系数 sigma D * tau / (1 - 1j * omega * tau) return sigma这段代码里mu是化学势实验上通过栅压或者掺杂调节vf是费米速度不同材料差别很大常见范围在 10^5 到 10^6 m/stau是弛豫时间和杂质浓度、温度有关。注意这里用的是简化 Drude 形式只适合低频、带内主导的情况。如果你要算中远红外响应必须把带间项加回来否则会低估高频电导率。提示三维狄拉克半金属的 Drude 权重和载流子浓度 n 成正比而 n 又和 mu³ 成正比。所以化学势稍微动一点低频电导率就会明显变化。做栅控实验时栅压范围要算够。3. 用 Kubo 公式算狄拉克电导率从带间跃迁到有限温度修正3.1 带间电导率的积分长什么样带间项来自价带到导带的垂直跃迁。在清洁极限下三维狄拉克半金属的带间电导率实部可以写成def dirac_interband_sigma_real(omega, mu, temp, vf): 三维狄拉克半金属带间电导率实部清洁极限零温近似 omega: 角频率 (rad/s) mu: 化学势 (eV) temp: 温度 (K) vf: 费米速度 (m/s) hbar 1.0545718e-34 e 1.6021766e-19 kB 8.617333e-5 # eV/K # 零温下的带间电导率omega 2*mu 时才有贡献 if omega * hbar 2 * abs(mu): return 0.0 # 三维狄拉克带间电导率趋近常数 sigma_inter e**2 * omega / (12 * np.pi * hbar * vf) return sigma_inter这个表达式最反直觉的地方是当光子能量超过两倍化学势时带间电导率随频率线性增长但除以频率后得到的“有效电导率”趋近常数。换句话说三维狄拉克半金属在低频带间吸收是平的不像半导体那样有吸收边。这个平带吸收在实验上表现为反射率在低频段有一个宽平台拟合时如果只放 Drude 项残差会很大。3.2 有限温度怎么加进去零温公式在室温下不够用因为热激发会让费米面变模糊。标准做法是在带间项里乘一个费米-狄拉克分布差def fermi_diff(E, mu, temp): 费米-狄拉克分布差 f(E-mu) - f(Emu) kB 8.617333e-5 # eV/K if temp 1e-6: return 1.0 if abs(E) abs(mu) else 0.0 f1 1.0 / (np.exp((E - mu) / (kB * temp)) 1) f2 1.0 / (np.exp((E mu) / (kB * temp)) 1) return f1 - f2然后在带间电导率的积分里把fermi_diff乘进去再对能量积分。实操中为了省时间很多人直接用零温公式乘一个温度修正因子但温度超过 100 K 以后误差会超过 20%。如果你的实验温度在液氮以上建议老老实实做数值积分。3.3 参数怎么取费米速度、化学势、弛豫时间费米速度 v_F 是狄拉克半金属最核心的参数。常见做法是用角分辨光电子能谱测出狄拉克锥的斜率或者用红外光谱拟合 Drude 权重反推。不同材料的 v_F 差别很大比如 Cd₃As₂ 的 v_F 大约在 10^6 m/s 量级而 Na₃Bi 稍低一些。如果你手头没有实验值可以从第一性原理能带计算里提取但要注意计算出的 v_F 通常比实验值高 10% 到 30%因为交换关联势和实验有偏差。化学势 mu 在输运实验里通常由载流子浓度反推。三维狄拉克半金属的载流子浓度和 mu 的关系是 n mu³ / (3π²ħ³v_F³)。如果你知道霍尔浓度可以直接解出 mu。注意这个关系只在 mu 远大于温度展宽时成立如果 mu 接近零载流子进入本征激发区公式要换成热激发形式。弛豫时间 tau 是最难确定的参数。它和杂质浓度、温度、样品质量都有关。常见做法是从直流电导率反推sigma_DC D * tau其中 D 是 Drude 权重。如果你测了直流电阻率和霍尔浓度就能算出 tau。但要注意直流电导率里可能还有带间项的贡献特别是当 mu 比较小的时候。稳妥的做法是同时拟合直流和低频交流电导率让 tau 和 mu 一起优化。注意三维狄拉克半金属的带间电导率在低频趋近常数这个常数和 tau 无关。所以如果你只拟合低频反射率tau 和 mu 可能会简并拟合结果不唯一。最好把中红外频段的数据也放进去。4. 狄拉克半金属电导率的避坑与排查五个血泪教训4.1 现象拟合出的费米速度和文献差一倍原因把带间项的贡献误算进了 Drude 权重。三维狄拉克半金属的带间电导率在低频是平的如果你用纯 Drude 模型拟合拟合算法会把带间平台当成一个超宽的 Drude 峰导致 Drude 权重被高估反推的 v_F 和 mu 全偏。解决拟合时强制加一个频率无关的带间项或者把拟合频段限制在带间贡献可以忽略的低频区。判断标准是如果拟合残差在低频段有系统性的正偏差说明带间项没算进去。4.2 现象低温下电导率随温度升高反而下降原因半金属里载流子冻析和缺陷散射的竞争。低温下热激发载流子少本征载流子浓度低电导率应该下降。但如果样品里有大量缺陷缺陷散射的迁移率随温度升高而下降可能把载流子浓度的增加抵消掉甚至反超。解决先测霍尔浓度随温度的变化确认载流子浓度是不是真的在增加。如果浓度在增加但电导率在下降说明迁移率主导这时候要重点分析缺陷散射机制而不是硬套狄拉克电导率公式。4.3 现象化学势调不动栅压加了几十伏电导率没反应原因三维狄拉克半金属的载流子浓度通常比二维材料高两三个数量级栅极电容根本调不动费米面。如果你用二维材料的栅控经验来套三维狄拉克半金属会发现栅压效率极低。解决要么用离子液体栅极提高电容要么直接换掺杂方式比如元素替代或者表面吸附。如果坚持用固态栅极先算一下需要多少栅压才能把 mu 移动 10 meV再判断可行性。4.4 现象带间电导率公式算出来的值比实验大一个量级原因零温公式没有考虑温度展宽和杂质展宽。零温下带间跃迁在 omega 2mu 处有一个尖角但实际样品里这个尖角会被温度展宽和杂质散射抹平峰值降低、宽度增加。如果你直接拿零温公式和实验比高频段会高估。解决在带间电导率里加一个洛伦兹展宽因子展宽宽度取 max(kB*T, ħ/tau)。这个修正对低频影响不大但对中红外频段的拟合至关重要。4.5 现象不同文献报道的同一材料电导率差好几倍原因样品质量、费米速度定义、化学势提取方法都不一样。有的文献用第一性原理的 v_F有的用实验拟合值有的把 mu 定义在狄拉克点有的定义在费米面。这些定义差异会导致电导率数值差很多。解决比较不同文献时先统一 v_F 和 mu 的定义。最稳妥的做法是直接比较无量纲量比如 sigma / (e²/ħ) 和 omega / mu 的关系这样能消掉大部分参数定义带来的差异。5. 进阶用法用光电导率验证狄拉克半金属的线性色散5.1 为什么光电导率比直流电导率更能说明问题直流电导率只能告诉你载流子浓度和迁移率的乘积没法区分抛物带和线性色散。光电导率不一样它的频率依赖直接反映态密度和跃迁矩阵元。三维狄拉克半金属的带间光电导率在低频趋近常数而普通半导体的带间光电导率在吸收边以上按 (ħω - Eg)^(1/2) 增长。这个定性差别是判断一个材料是不是狄拉克半金属的关键证据。实操中你可以用傅里叶变换红外光谱测反射率然后通过 Kramers-Kronig 变换提取光电导率。如果光电导率实部在低频有一个平台而且平台高度和温度、化学势的关系符合狄拉克公式那基本可以确认线性色散。5.2 一个可复现的验证流程第一步测变温反射率温度范围覆盖 10 K 到 300 K。第二步用金镜或者铝镜做参考确保反射率绝对值准确。第三步做 Kramers-Kronig 变换注意高频外推要用 x 射线散射数据或者原子散射因子。第四步把提取的光电导率和狄拉克模型对比重点看三个特征低频平台、平台高度随温度的变化、带间吸收边随化学势的移动。def drude_dirac_fit(omega, sigma_exp, mu, temp, vf, tau, sigma_inter_const): 拟合光电导率Drude 带间常数项 omega: 角频率数组 sigma_exp: 实验光电导率实部 mu, temp, vf, tau: 拟合参数 sigma_inter_const: 带间平台高度 sigma_drude dirac_intraband_sigma(omega, mu, temp, vf, tau) sigma_total sigma_drude.real sigma_inter_const residual sigma_exp - sigma_total return residual拟合时sigma_inter_const不要设成自由参数先用零温公式算一个初值再让它在 ±20% 范围内浮动。tau和mu的初值从直流输运数据里取。如果拟合残差在低频段还有结构说明模型里少了东西可能是声子吸收或者多带效应。5.3 我自己的习惯先看量纲再调参数做了这么多年拟合我最大的教训是不要一上来就扔给算法。先把实验数据的量纲和公式对一遍确认单位制统一。狄拉克电导率公式里 ħ、e、v_F 的单位特别容易混有人用 eV 有人用 J有人用 cm/s 有人用 m/s差一个量级结果就全错了。我一般会先手算一个特征频率下的电导率和实验值比个数量级确认没问题再让算法细调。这个习惯帮我省了很多后悔药。希望帮到你。本文还有配套的精品资源点击获取 狄拉克电导率从半金属到狄拉克半金属的实操指南
简介这份资料面向凝聚态物理与材料计算方向的学习者聚焦狄拉克半金属中狄拉克电导率的理论分析与数值计算。包内共2个文件含1个MATLAB脚本与1篇PDF文献压缩包约1.12MB。脚本可用于模拟电导率其中对虚部采用阶跃函数近似处理便于简化复数频率下的吸收与散射计算PDF则围绕三维狄拉克半金属薄膜的介电响应与新颖电磁模式展开涉及狄拉克锥、狄拉克点及边界效应等关键概念。已有710人学习下载适合希望从能带结构出发理解零质量电子输运、并借助脚本复现电导率曲线的读者参考也可作为薄膜电磁模式研究的入门材料。1. 狄拉克电导率到底在算什么从半金属到狄拉克半金属的一条主线第一次接触“狄拉克电导率”这个词很多人会以为它只是石墨烯的专属名词。其实只要体系里存在线性色散的低能激发也就是能量和动量近似成正比、能带在一点或一条线上简并那套用狄拉克方程描述的电导率公式就能派上用场。狄拉克半金属正是这样一类材料它的导带和价带在布里渊区里以四重简并的狄拉克点接触低能载流子表现得像无质量的相对论粒子。于是问题就变成——这种“无质量”的载流子在电场里到底怎么输运电导率随频率、温度、化学势怎么变。搞懂这条主线你才能判断一个材料值不值得做器件、该测哪个频段、拟合时该锁哪些参数。这篇笔记就按“概念立住 → 公式落地 → 参数怎么取 → 坑在哪”的顺序把狄拉克电导率从半金属一路讲到狄拉克半金属的实操路径。2. 从半金属到狄拉克半金属能带、态密度与电导率的物理起点2.1 半金属和普通金属的差别在哪普通金属的费米面穿过能带态密度在费米能级处是个有限值载流子有有效质量电导率的 Drude 形式大家都很熟。半金属不一样它的导带底和价带顶在能量上轻微重叠费米能级处的态密度趋近于零但不完全为零载流子浓度比普通金属低两三个数量级。这个“低载流子浓度 小费米面”的组合直接导致半金属的电导率对温度、杂质和缺陷特别敏感。你如果在实验里看到电阻率随温度先降后升、低温残余电阻比很高大概率就是半金属体系里载流子被冻析或者缺陷散射主导了。狄拉克半金属是半金属里更特殊的一支。它的导带和价带不是简单重叠而是在动量空间某些高对称点或高对称线上形成四重简并的狄拉克点。低能激发满足 E ±v_F |p|v_F 是费米速度p 是相对狄拉克点的动量。这个线性色散让态密度在狄拉克点附近按 |E|² 变化三维或 |E|二维和普通抛物带金属的 √E 完全不同。态密度的形状一变电导率对化学势、频率和温度的依赖关系就全变了。2.2 为什么狄拉克半金属的电导率值得单独算普通金属的光电导率在低频下基本是平的由 Drude 权重决定。狄拉克半金属不一样它的带间跃迁在低频就有贡献而且带内和带间两部分会随频率、温度、化学势互相竞争。更麻烦的是三维狄拉克半金属的带间电导率在低频极限下趋近于一个常数这个常数只和费米速度、普朗克常数有关和载流子浓度无关。这意味着你没法靠调掺杂把低频电导率压下去器件的暗电流、热噪声都会受这个常数限制。做红外探测器或者高频器件的人如果忽略这一点设计出来的响应度可能和预期差一个量级。从建模角度狄拉克半金属的电导率通常写成带内 Drude 项加带间跃迁项。带内项描述费米面附近的载流子被电场加速带间项描述电子从价带被光子激发到导带。在有限温度下带间项还要乘一个费米-狄拉克分布差把热激发和泡利阻塞都算进去。这套公式在石墨烯里已经验证得很充分搬到三维狄拉克半金属时主要变化是态密度从线性变成二次积分维度从二变三。2.3 一个最小可用的电导率表达式实操里最常用的形式是 Kubo 公式在清洁极限下的结果。对三维狄拉克半金属带内电导率可以写成import numpy as np def dirac_intraband_sigma(omega, mu, temp, vf, tau): 三维狄拉克半金属带内电导率简化 Drude 形式 omega: 角频率 (rad/s) mu: 化学势 (eV) temp: 温度 (K) vf: 费米速度 (m/s) tau: 弛豫时间 (s) hbar 1.0545718e-34 kB 8.617333e-5 # eV/K e 1.6021766e-19 # 有限温度下的载流子浓度近似 if mu 0: n (mu**3) / (3 * np.pi**2 * (hbar * vf)**3) else: n 0.0 # Drude 权重 D e**2 * vf * n / (hbar * 3) # 三维狄拉克的 Drude 权重系数 sigma D * tau / (1 - 1j * omega * tau) return sigma这段代码里mu是化学势实验上通过栅压或者掺杂调节vf是费米速度不同材料差别很大常见范围在 10^5 到 10^6 m/stau是弛豫时间和杂质浓度、温度有关。注意这里用的是简化 Drude 形式只适合低频、带内主导的情况。如果你要算中远红外响应必须把带间项加回来否则会低估高频电导率。提示三维狄拉克半金属的 Drude 权重和载流子浓度 n 成正比而 n 又和 mu³ 成正比。所以化学势稍微动一点低频电导率就会明显变化。做栅控实验时栅压范围要算够。3. 用 Kubo 公式算狄拉克电导率从带间跃迁到有限温度修正3.1 带间电导率的积分长什么样带间项来自价带到导带的垂直跃迁。在清洁极限下三维狄拉克半金属的带间电导率实部可以写成def dirac_interband_sigma_real(omega, mu, temp, vf): 三维狄拉克半金属带间电导率实部清洁极限零温近似 omega: 角频率 (rad/s) mu: 化学势 (eV) temp: 温度 (K) vf: 费米速度 (m/s) hbar 1.0545718e-34 e 1.6021766e-19 kB 8.617333e-5 # eV/K # 零温下的带间电导率omega 2*mu 时才有贡献 if omega * hbar 2 * abs(mu): return 0.0 # 三维狄拉克带间电导率趋近常数 sigma_inter e**2 * omega / (12 * np.pi * hbar * vf) return sigma_inter这个表达式最反直觉的地方是当光子能量超过两倍化学势时带间电导率随频率线性增长但除以频率后得到的“有效电导率”趋近常数。换句话说三维狄拉克半金属在低频带间吸收是平的不像半导体那样有吸收边。这个平带吸收在实验上表现为反射率在低频段有一个宽平台拟合时如果只放 Drude 项残差会很大。3.2 有限温度怎么加进去零温公式在室温下不够用因为热激发会让费米面变模糊。标准做法是在带间项里乘一个费米-狄拉克分布差def fermi_diff(E, mu, temp): 费米-狄拉克分布差 f(E-mu) - f(Emu) kB 8.617333e-5 # eV/K if temp 1e-6: return 1.0 if abs(E) abs(mu) else 0.0 f1 1.0 / (np.exp((E - mu) / (kB * temp)) 1) f2 1.0 / (np.exp((E mu) / (kB * temp)) 1) return f1 - f2然后在带间电导率的积分里把fermi_diff乘进去再对能量积分。实操中为了省时间很多人直接用零温公式乘一个温度修正因子但温度超过 100 K 以后误差会超过 20%。如果你的实验温度在液氮以上建议老老实实做数值积分。3.3 参数怎么取费米速度、化学势、弛豫时间费米速度 v_F 是狄拉克半金属最核心的参数。常见做法是用角分辨光电子能谱测出狄拉克锥的斜率或者用红外光谱拟合 Drude 权重反推。不同材料的 v_F 差别很大比如 Cd₃As₂ 的 v_F 大约在 10^6 m/s 量级而 Na₃Bi 稍低一些。如果你手头没有实验值可以从第一性原理能带计算里提取但要注意计算出的 v_F 通常比实验值高 10% 到 30%因为交换关联势和实验有偏差。化学势 mu 在输运实验里通常由载流子浓度反推。三维狄拉克半金属的载流子浓度和 mu 的关系是 n mu³ / (3π²ħ³v_F³)。如果你知道霍尔浓度可以直接解出 mu。注意这个关系只在 mu 远大于温度展宽时成立如果 mu 接近零载流子进入本征激发区公式要换成热激发形式。弛豫时间 tau 是最难确定的参数。它和杂质浓度、温度、样品质量都有关。常见做法是从直流电导率反推sigma_DC D * tau其中 D 是 Drude 权重。如果你测了直流电阻率和霍尔浓度就能算出 tau。但要注意直流电导率里可能还有带间项的贡献特别是当 mu 比较小的时候。稳妥的做法是同时拟合直流和低频交流电导率让 tau 和 mu 一起优化。注意三维狄拉克半金属的带间电导率在低频趋近常数这个常数和 tau 无关。所以如果你只拟合低频反射率tau 和 mu 可能会简并拟合结果不唯一。最好把中红外频段的数据也放进去。4. 狄拉克半金属电导率的避坑与排查五个血泪教训4.1 现象拟合出的费米速度和文献差一倍原因把带间项的贡献误算进了 Drude 权重。三维狄拉克半金属的带间电导率在低频是平的如果你用纯 Drude 模型拟合拟合算法会把带间平台当成一个超宽的 Drude 峰导致 Drude 权重被高估反推的 v_F 和 mu 全偏。解决拟合时强制加一个频率无关的带间项或者把拟合频段限制在带间贡献可以忽略的低频区。判断标准是如果拟合残差在低频段有系统性的正偏差说明带间项没算进去。4.2 现象低温下电导率随温度升高反而下降原因半金属里载流子冻析和缺陷散射的竞争。低温下热激发载流子少本征载流子浓度低电导率应该下降。但如果样品里有大量缺陷缺陷散射的迁移率随温度升高而下降可能把载流子浓度的增加抵消掉甚至反超。解决先测霍尔浓度随温度的变化确认载流子浓度是不是真的在增加。如果浓度在增加但电导率在下降说明迁移率主导这时候要重点分析缺陷散射机制而不是硬套狄拉克电导率公式。4.3 现象化学势调不动栅压加了几十伏电导率没反应原因三维狄拉克半金属的载流子浓度通常比二维材料高两三个数量级栅极电容根本调不动费米面。如果你用二维材料的栅控经验来套三维狄拉克半金属会发现栅压效率极低。解决要么用离子液体栅极提高电容要么直接换掺杂方式比如元素替代或者表面吸附。如果坚持用固态栅极先算一下需要多少栅压才能把 mu 移动 10 meV再判断可行性。4.4 现象带间电导率公式算出来的值比实验大一个量级原因零温公式没有考虑温度展宽和杂质展宽。零温下带间跃迁在 omega 2mu 处有一个尖角但实际样品里这个尖角会被温度展宽和杂质散射抹平峰值降低、宽度增加。如果你直接拿零温公式和实验比高频段会高估。解决在带间电导率里加一个洛伦兹展宽因子展宽宽度取 max(kB*T, ħ/tau)。这个修正对低频影响不大但对中红外频段的拟合至关重要。4.5 现象不同文献报道的同一材料电导率差好几倍原因样品质量、费米速度定义、化学势提取方法都不一样。有的文献用第一性原理的 v_F有的用实验拟合值有的把 mu 定义在狄拉克点有的定义在费米面。这些定义差异会导致电导率数值差很多。解决比较不同文献时先统一 v_F 和 mu 的定义。最稳妥的做法是直接比较无量纲量比如 sigma / (e²/ħ) 和 omega / mu 的关系这样能消掉大部分参数定义带来的差异。5. 进阶用法用光电导率验证狄拉克半金属的线性色散5.1 为什么光电导率比直流电导率更能说明问题直流电导率只能告诉你载流子浓度和迁移率的乘积没法区分抛物带和线性色散。光电导率不一样它的频率依赖直接反映态密度和跃迁矩阵元。三维狄拉克半金属的带间光电导率在低频趋近常数而普通半导体的带间光电导率在吸收边以上按 (ħω - Eg)^(1/2) 增长。这个定性差别是判断一个材料是不是狄拉克半金属的关键证据。实操中你可以用傅里叶变换红外光谱测反射率然后通过 Kramers-Kronig 变换提取光电导率。如果光电导率实部在低频有一个平台而且平台高度和温度、化学势的关系符合狄拉克公式那基本可以确认线性色散。5.2 一个可复现的验证流程第一步测变温反射率温度范围覆盖 10 K 到 300 K。第二步用金镜或者铝镜做参考确保反射率绝对值准确。第三步做 Kramers-Kronig 变换注意高频外推要用 x 射线散射数据或者原子散射因子。第四步把提取的光电导率和狄拉克模型对比重点看三个特征低频平台、平台高度随温度的变化、带间吸收边随化学势的移动。def drude_dirac_fit(omega, sigma_exp, mu, temp, vf, tau, sigma_inter_const): 拟合光电导率Drude 带间常数项 omega: 角频率数组 sigma_exp: 实验光电导率实部 mu, temp, vf, tau: 拟合参数 sigma_inter_const: 带间平台高度 sigma_drude dirac_intraband_sigma(omega, mu, temp, vf, tau) sigma_total sigma_drude.real sigma_inter_const residual sigma_exp - sigma_total return residual拟合时sigma_inter_const不要设成自由参数先用零温公式算一个初值再让它在 ±20% 范围内浮动。tau和mu的初值从直流输运数据里取。如果拟合残差在低频段还有结构说明模型里少了东西可能是声子吸收或者多带效应。5.3 我自己的习惯先看量纲再调参数做了这么多年拟合我最大的教训是不要一上来就扔给算法。先把实验数据的量纲和公式对一遍确认单位制统一。狄拉克电导率公式里 ħ、e、v_F 的单位特别容易混有人用 eV 有人用 J有人用 cm/s 有人用 m/s差一个量级结果就全错了。我一般会先手算一个特征频率下的电导率和实验值比个数量级确认没问题再让算法细调。这个习惯帮我省了很多后悔药。希望帮到你。本文还有配套的精品资源点击获取