簡(jiǎn)介這份資源面向凝聚態(tài)物理實(shí)驗(yàn)方向的研究人員與具備MATLAB基礎(chǔ)的高年級(jí)學(xué)生用于處理霍爾效應(yīng)測(cè)量中的雙載流子擬合問(wèn)題。當(dāng)半導(dǎo)體等材料中電子與空穴同時(shí)貢獻(xiàn)電導(dǎo)時(shí)單載流子模型往往無(wú)法準(zhǔn)確解釋霍爾電壓與電阻隨磁場(chǎng)的變化需要將兩類(lèi)載流子的密度與遷移率分別建模求解。壓縮包內(nèi)共1個(gè)文件為MATLAB腳本.m整體約3KB體積輕量便于直接嵌入現(xiàn)有分析流程。腳本預(yù)期覆蓋實(shí)驗(yàn)數(shù)據(jù)導(dǎo)入、異常值清洗與預(yù)處理、基于非線(xiàn)性最小二乘的雙載流子擬合、霍爾電阻隨磁場(chǎng)變化及電子空穴貢獻(xiàn)的可視化以及載流子濃度與遷移率等關(guān)鍵參數(shù)的輸出。已有273人學(xué)習(xí)適合需要快速搭建雙載流子分析鏈路、對(duì)照理論模型驗(yàn)證實(shí)驗(yàn)結(jié)果的科研場(chǎng)景。1. 雙載流子霍爾擬合從一張 Rxx 曲線(xiàn)里拆出兩種載流子的濃度和遷移率你拿到一張霍爾測(cè)試的 Rxx-B 曲線(xiàn)發(fā)現(xiàn)它不按套路出牌——低場(chǎng)斜率是正的高場(chǎng)斜率卻翻成負(fù)的或者曲線(xiàn)明顯往一邊鼓。單載流子模型擬合出來(lái)的載流子濃度和遷移率換一個(gè)磁場(chǎng)區(qū)間就變一個(gè)值根本沒(méi)法用。這時(shí)候大概率不是設(shè)備壞了而是樣品里同時(shí)存在兩種載流子電子和空穴都在導(dǎo)電各自有不同的濃度和遷移率霍爾信號(hào)是它們疊加的結(jié)果。hall_twocarrier_rxx 這個(gè)方向要解決的就是這件事從一條 Rxx或 Hall 電阻隨磁場(chǎng)變化的曲線(xiàn)里用雙載流子模型把兩套載流子參數(shù)拆出來(lái)。它適合做半導(dǎo)體薄膜、二維材料、拓?fù)浣^緣體、有機(jī)導(dǎo)體這類(lèi)多載流子共存體系的同學(xué)尤其是那些單載流子擬合殘差大、變溫?cái)?shù)據(jù)對(duì)不上的情況。載流子濃度這個(gè)熱詞在這里不是單一數(shù)字而是一組隨磁場(chǎng)和溫度變化的擬合參數(shù)理解這一點(diǎn)是后面所有操作的前提。2. 雙載流子模型的物理圖像與公式推導(dǎo)為什么 Rxx 會(huì)拐彎2.1 兩種載流子各自怎么貢獻(xiàn)電導(dǎo)先回到最基本的 Drude 模型。單一載流子體系里電導(dǎo)率張量的兩個(gè)分量是σ_xx n e μ / (1 (μ B)^2) σ_xy n e μ^2 B / (1 (μ B)^2)其中 n 是載流子濃度μ 是遷移率B 是磁感應(yīng)強(qiáng)度?；魻栯娮杪?ρ_xy 和橫向電阻率 ρ_xx 通過(guò)張量求逆得到ρ_xx σ_xx / (σ_xx^2 σ_xy^2) ρ_xy σ_xy / (σ_xx^2 σ_xy^2)當(dāng)體系里同時(shí)有兩種載流子時(shí)總電導(dǎo)是各自電導(dǎo)的簡(jiǎn)單相加σ_xx_total σ_xx1 σ_xx2 σ_xy_total σ_xy1 σ_xy2注意這里有個(gè)容易翻車(chē)的地方電阻率不能直接相加必須先把各自的電導(dǎo)率算出來(lái)求和之后再求逆得到總電阻率。很多人第一次寫(xiě)代碼就是在這里把 ρ_xx 直接相加結(jié)果擬合出來(lái)的曲線(xiàn)形狀完全不對(duì)。雙載流子模型的核心參數(shù)一共四個(gè)n1、μ1、n2、μ2。載流子類(lèi)型電子還是空穴通過(guò) μ 的符號(hào)或者電荷符號(hào)來(lái)區(qū)分。通常約定電子遷移率為負(fù)、空穴為正或者反過(guò)來(lái)只要在代碼里保持一致就行。2.2 從 Rxx 曲線(xiàn)形狀反推載流子類(lèi)型Rxx 隨磁場(chǎng)的變化形狀本身就攜帶信息。如果兩種載流子的遷移率差異很大Rxx 會(huì)在高場(chǎng)出現(xiàn)明顯的飽和或者非單調(diào)變化。具體來(lái)說(shuō)當(dāng)兩種載流子濃度相當(dāng)、遷移率接近時(shí)Rxx 在低場(chǎng)緩慢上升高場(chǎng)趨于飽和曲線(xiàn)看起來(lái)像被壓扁的拋物線(xiàn)。當(dāng)一種載流子濃度遠(yuǎn)高于另一種但遷移率低很多時(shí)低場(chǎng)由高遷移率載流子主導(dǎo)Rxx 上升快高場(chǎng)低遷移率載流子貢獻(xiàn)變大Rxx 趨于平緩。如果兩種載流子類(lèi)型相反電子和空穴霍爾電阻率 ρ_xy 會(huì)出現(xiàn)符號(hào)反轉(zhuǎn)對(duì)應(yīng)的 Rxx 也會(huì)在反轉(zhuǎn)場(chǎng)附近出現(xiàn)一個(gè)特征性的鼓包。這些形狀特征就是擬合時(shí)判斷初值好不好的依據(jù)。我一般會(huì)先目測(cè)曲線(xiàn)拐點(diǎn)大概在哪個(gè)磁場(chǎng)然后估算兩種載流子的遷移率量級(jí)再讓擬合去細(xì)化。2.3 擬合殘差為什么比單載流子小單載流子模型只有兩個(gè)自由參數(shù)n 和 μ雙載流子有四個(gè)。參數(shù)多了擬合能力自然強(qiáng)但這不是隨便加參數(shù)的借口。判斷是否真的需要雙載流子模型看兩個(gè)指標(biāo)第一單載流子擬合的殘差是否隨磁場(chǎng)系統(tǒng)性偏大。如果殘差在低場(chǎng)和高場(chǎng)符號(hào)相反說(shuō)明模型缺了一項(xiàng)。第二擬合出來(lái)的遷移率是否隨擬合區(qū)間變化。單載流子模型擬合 0-1T 和 0-9T 得到的 μ 差了兩倍以上基本可以確定是多載流子體系。提示不要因?yàn)殡p載流子擬合殘差更小就直接采用。參數(shù)多了容易過(guò)擬合尤其是數(shù)據(jù)點(diǎn)少、噪聲大的時(shí)候。先用單載流子試殘差有系統(tǒng)性結(jié)構(gòu)再上雙載流子。3. 用 Python 實(shí)現(xiàn)雙載流子擬合從數(shù)據(jù)讀入到參數(shù)輸出3.1 數(shù)據(jù)準(zhǔn)備與預(yù)處理假設(shè)你從霍爾測(cè)試系統(tǒng)導(dǎo)出的數(shù)據(jù)是兩列磁場(chǎng) B單位 T和 Rxx單位 Ω 或 Ω/sq。常見(jiàn)格式是 CSV 或 Excel。先做幾件事去掉零場(chǎng)附近的異常點(diǎn)有些設(shè)備在 B0 附近有剩磁或偏移檢查磁場(chǎng)是否對(duì)稱(chēng)如果只測(cè)了正場(chǎng)最好補(bǔ)測(cè)負(fù)場(chǎng)或者至少確認(rèn)曲線(xiàn)在正負(fù)場(chǎng)下對(duì)稱(chēng)。import numpy as np import pandas as pd from scipy.optimize import curve_fit import matplotlib.pyplot as plt # 讀入數(shù)據(jù)假設(shè)文件有兩列B_T 和 Rxx_ohm df pd.read_csv(hall_rxx_data.csv) B df[B_T].values Rxx_exp df[Rxx_ohm].values # 去掉零場(chǎng)附近可能受剩磁影響的點(diǎn)可選 mask np.abs(B) 0.05 B B[mask] Rxx_exp Rxx_exp[mask] # 按磁場(chǎng)絕對(duì)值排序方便后續(xù)畫(huà)圖 idx np.argsort(B) B B[idx] Rxx_exp Rxx_exp[idx] print(f數(shù)據(jù)點(diǎn)數(shù): {len(B)}) print(f磁場(chǎng)范圍: {B.min():.2f} T 到 {B.max():.2f} T) print(fRxx 范圍: {Rxx_exp.min():.2f} 到 {Rxx_exp.max():.2f} Ω)這段代碼做了三件事讀數(shù)據(jù)、去零場(chǎng)附近點(diǎn)、排序。去零場(chǎng)附近點(diǎn)是因?yàn)楹芏嗷魻栂到y(tǒng)在 B0 附近有剩磁Rxx 會(huì)出現(xiàn)一個(gè)不合理的凹陷或尖峰留著會(huì)帶偏擬合。閾值 0.05 T 不是固定的如果你的最高場(chǎng)只有 1 T這個(gè)閾值要相應(yīng)減小。3.2 雙載流子電導(dǎo)率模型函數(shù)這是整個(gè)擬合的核心。輸入磁場(chǎng) B 和四個(gè)參數(shù)輸出模型預(yù)測(cè)的 Rxx。def two_carrier_conductivity(B, n1, mu1, n2, mu2): 雙載流子電導(dǎo)率模型。 n1, n2: 載流子濃度單位 cm^-3 mu1, mu2: 遷移率單位 cm^2/(V·s)電子為負(fù)、空穴為正 B: 磁場(chǎng)單位 T 返回: Rxx 模型值單位 Ω假設(shè)樣品為二維用方塊電阻 e 1.602e-19 # 元電荷C # 把濃度從 cm^-3 轉(zhuǎn)成 m^-3 n1_si n1 * 1e6 n2_si n2 * 1e6 # 遷移率從 cm^2/(V·s) 轉(zhuǎn)成 m^2/(V·s) mu1_si mu1 * 1e-4 mu2_si mu2 * 1e-4 # 各自電導(dǎo)率 sigma_xx1 n1_si * e * np.abs(mu1_si) / (1 (mu1_si * B)**2) sigma_xy1 n1_si * e * mu1_si**2 * B / (1 (mu1_si * B)**2) sigma_xx2 n2_si * e * np.abs(mu2_si) / (1 (mu2_si * B)**2) sigma_xy2 n2_si * e * mu2_si**2 * B / (1 (mu2_si * B)**2) # 總電導(dǎo) sigma_xx sigma_xx1 sigma_xx2 sigma_xy sigma_xy1 sigma_xy2 # 求逆得到電阻率 rho_xx sigma_xx / (sigma_xx**2 sigma_xy**2) return rho_xx參數(shù)說(shuō)明n1、n2 用 cm^-3 是文獻(xiàn)里最常見(jiàn)的單位代碼內(nèi)部轉(zhuǎn)成 SI 單位計(jì)算。μ 的符號(hào)區(qū)分載流子類(lèi)型電子為負(fù)、空穴為正。注意 σ_xx 里用的是 |μ|因?yàn)殡妼?dǎo)率是正值符號(hào)只影響 σ_xy 的方向。如果你的樣品是三維體材料Rxx 和 ρ_xx 之間還差一個(gè)厚度因子需要把厚度乘進(jìn)去。3.3 初值估計(jì)與擬合調(diào)用雙載流子擬合對(duì)初值敏感直接扔隨機(jī)初值大概率不收斂。我一般用兩步法先目測(cè)曲線(xiàn)估一個(gè)大概范圍再用網(wǎng)格搜索找最優(yōu)初值。# 第一步根據(jù)曲線(xiàn)形狀估初值 # 假設(shè)低場(chǎng)斜率對(duì)應(yīng)高遷移率載流子高場(chǎng)飽和對(duì)應(yīng)低遷移率載流子 # 這里給一組典型初值實(shí)際要根據(jù)你的數(shù)據(jù)調(diào)整 p0 [1e18, 1000, 1e19, 100] # n1, mu1, n2, mu2 # 參數(shù)邊界濃度 1e15 到 1e22遷移率 1 到 10000 bounds ([1e15, -10000, 1e15, -10000], [1e22, 10000, 1e22, 10000]) # 第二步擬合 popt, pcov curve_fit(two_carrier_conductivity, B, Rxx_exp, p0p0, boundsbounds, maxfev20000) n1_fit, mu1_fit, n2_fit, mu2_fit popt perr np.sqrt(np.diag(pcov)) print(f載流子1: n {n1_fit:.3e} cm^-3, mu {mu1_fit:.1f} cm^2/(V·s)) print(f載流子2: n {n2_fit:.3e} cm^-3, mu {mu2_fit:.1f} cm^2/(V·s)) print(f擬合誤差: n1 ±{perr[0]:.2e}, mu1 ±{perr[1]:.1f}, fn2 ±{perr[2]:.2e}, mu2 ±{perr[3]:.1f}) # 畫(huà)圖對(duì)比 B_fit np.linspace(B.min(), B.max(), 500) Rxx_fit two_carrier_conductivity(B_fit, *popt) plt.figure(figsize(8, 5)) plt.plot(B, Rxx_exp, o, markersize3, label實(shí)驗(yàn)數(shù)據(jù)) plt.plot(B_fit, Rxx_fit, -, linewidth2, label雙載流子擬合) plt.xlabel(磁場(chǎng) B (T)) plt.ylabel(Rxx (Ω)) plt.legend() plt.tight_layout() plt.savefig(two_carrier_fit.png, dpi150) plt.show()這段代碼的關(guān)鍵在初值和邊界。p0 里的四個(gè)數(shù)不是隨便寫(xiě)的n1 和 mu1 對(duì)應(yīng)低場(chǎng)主導(dǎo)的載流子n2 和 mu2 對(duì)應(yīng)高場(chǎng)主導(dǎo)的載流子。如果你不確定哪個(gè)是哪個(gè)可以交換 p0 里兩組參數(shù)的位置再擬合一次看哪個(gè)殘差小。bounds 給了一個(gè)很寬的物理合理范圍防止擬合跑到負(fù)濃度或者離譜的遷移率上去。3.4 擬合結(jié)果的可信度檢查擬合完不能只看殘差平方和要做幾項(xiàng)檢查第一把擬合曲線(xiàn)和實(shí)驗(yàn)數(shù)據(jù)畫(huà)在一起目測(cè)有沒(méi)有系統(tǒng)性偏離。如果低場(chǎng)擬合好、高場(chǎng)差說(shuō)明高場(chǎng)載流子參數(shù)沒(méi)抓準(zhǔn)。第二看參數(shù)的誤差棒。如果某個(gè)參數(shù)的誤差比參數(shù)本身還大說(shuō)明這個(gè)參數(shù)沒(méi)被數(shù)據(jù)約束住可能是過(guò)擬合。第三換一個(gè)磁場(chǎng)區(qū)間重新擬合看參數(shù)是否穩(wěn)定。比如用 0-6T 和 0-9T 分別擬合如果 n1 變了三倍說(shuō)明數(shù)據(jù)質(zhì)量不夠或者模型不對(duì)。第四檢查擬合出來(lái)的遷移率是否和文獻(xiàn)里同類(lèi)材料的量級(jí)一致。如果擬合出 10000 cm^2/(V·s) 的遷移率但你的材料是普通多晶薄膜大概率是過(guò)擬合了。4. 避坑與排查雙載流子擬合里最容易翻車(chē)的五個(gè)地方4.1 現(xiàn)象擬合不收斂報(bào)錯(cuò) “Optimal parameters not found”原因初值離真實(shí)值太遠(yuǎn)或者邊界設(shè)得太窄把真值排除在外了。雙載流子模型的參數(shù)空間是四維的而且 σ_xx 和 σ_xy 對(duì)參數(shù)的敏感度不一樣梯度下降很容易卡在局部極小值。解決先用單載流子模型擬合一遍把得到的 n 和 μ 作為其中一個(gè)載流子的初值另一個(gè)載流子的初值設(shè)成濃度差一個(gè)量級(jí)、遷移率差三到五倍。如果還不行用 scipy.optimize.differential_evolution 做全局搜索雖然慢但不容易卡住。4.2 現(xiàn)象擬合出來(lái)的兩種載流子濃度幾乎一樣遷移率也幾乎一樣原因數(shù)據(jù)里其實(shí)只有一種載流子雙載流子模型把一組參數(shù)拆成了兩組相近的相當(dāng)于用四個(gè)參數(shù)擬合兩個(gè)參數(shù)就能描述的數(shù)據(jù)。這是典型的過(guò)擬合。解決先做單載流子擬合看殘差有沒(méi)有系統(tǒng)性結(jié)構(gòu)。如果沒(méi)有就別用雙載流子。如果有檢查是不是測(cè)量噪聲導(dǎo)致的假結(jié)構(gòu)?？梢阅靡唤M已知單載流子的樣品數(shù)據(jù)跑一下雙載流子擬合看看它會(huì)不會(huì)也拆出兩組相近參數(shù)——如果會(huì)說(shuō)明你的擬合流程需要加正則化或者參數(shù)約束。4.3 現(xiàn)象Rxx 曲線(xiàn)在低場(chǎng)有個(gè)尖峰擬合怎么都對(duì)不上原因低場(chǎng)尖峰通常不是載流子輸運(yùn)的本征行為而是測(cè)量偽影。常見(jiàn)來(lái)源包括樣品接觸不良導(dǎo)致的非歐姆接觸、零場(chǎng)附近的剩磁、電流過(guò)大導(dǎo)致的自加熱、或者樣品本身在低場(chǎng)有磁阻效應(yīng)比如弱局域化。解決先檢查測(cè)量條件。降低測(cè)試電流看尖峰是否變小。如果尖峰在正負(fù)場(chǎng)下不對(duì)稱(chēng)基本可以確定是接觸問(wèn)題或剩磁。把低場(chǎng)數(shù)據(jù)通常 |B| 0.1T去掉再擬合不要硬擬合。4.4 現(xiàn)象擬合參數(shù)隨溫度劇烈跳動(dòng)沒(méi)有物理趨勢(shì)原因雙載流子擬合的參數(shù)相關(guān)性很強(qiáng)n1 和 n2 之間、μ1 和 μ2 之間存在 trade-off。數(shù)據(jù)信噪比不夠時(shí)擬合會(huì)在多組參數(shù)之間跳。解決把不同溫度的數(shù)據(jù)放在一起看趨勢(shì)。如果 n1 隨溫度升高而增加、n2 隨溫度升高而減少但總濃度守恒那可能是擬合把兩種載流子的角色互換了。可以在擬合時(shí)加約束比如固定 n1 n2 總濃度如果總濃度能從其他測(cè)量獨(dú)立得到或者固定一種載流子的遷移率隨溫度的變化規(guī)律。4.5 現(xiàn)象擬合殘差很小但參數(shù)物理上不合理比如遷移率超過(guò) 10^5原因參數(shù)邊界設(shè)得太寬擬合跑到了非物理區(qū)域?；蛘吣Ｐ屠锷倏紤]了串聯(lián)電阻、接觸電阻等寄生效應(yīng)。解決收緊邊界。常見(jiàn)半導(dǎo)體材料的遷移率室溫下很少超過(guò) 2000 cm^2/(V·s)二維材料可能到 10^4 量級(jí)但 10^5 基本不可能。把邊界上限設(shè)成 5×10^4 再擬合。如果擬合結(jié)果貼著邊界說(shuō)明數(shù)據(jù)里沒(méi)有足夠信息約束這個(gè)參數(shù)需要補(bǔ)測(cè)更高場(chǎng)的數(shù)據(jù)或者換測(cè)量方式。5. 進(jìn)階技巧用變溫 Rxx 數(shù)據(jù)交叉驗(yàn)證雙載流子參數(shù)單溫度下的雙載流子擬合只能給你一組參數(shù)但物理上是否自洽要靠變溫?cái)?shù)據(jù)來(lái)驗(yàn)證。我一般會(huì)測(cè) 5-10 個(gè)溫度點(diǎn)每個(gè)溫度下擬合出 n1、μ1、n2、μ2然后看它們隨溫度的變化趨勢(shì)是否符合預(yù)期。一個(gè)典型的驗(yàn)證方法是畫(huà) Arrhenius 圖把 n1 和 n2 取對(duì)數(shù)后對(duì) 1/T 作圖。如果兩種載流子都來(lái)自熱激發(fā)log(n) vs 1/T 應(yīng)該是直線(xiàn)斜率給出激活能。如果 n1 的激活能和 n2 差很多說(shuō)明它們來(lái)自不同的能級(jí)或不同的導(dǎo)電通道雙載流子模型是合理的。如果 n1 和 n2 的激活能幾乎一樣那可能它們其實(shí)是同一種載流子的不同分布模型需要簡(jiǎn)化。另一個(gè)技巧是用 Hall 電阻率 ρ_xy 的符號(hào)反轉(zhuǎn)場(chǎng)來(lái)約束參數(shù)。雙載流子體系中如果電子和空穴共存ρ_xy 會(huì)在某個(gè)磁場(chǎng)下過(guò)零。這個(gè)過(guò)零場(chǎng) B0 和兩種載流子的濃度、遷移率滿(mǎn)足一個(gè)近似關(guān)系B0 ≈ (n2 μ2^2 - n1 μ1^2) / (e (n2 μ2 - n1 μ1) ... )具體形式取決于符號(hào)約定但核心思路是過(guò)零場(chǎng)的位置對(duì)參數(shù)組合很敏感可以把它作為一個(gè)額外的約束條件加進(jìn)擬合。我通常的做法是先用 Rxx 擬合出一組參數(shù)然后檢查這組參數(shù)預(yù)測(cè)的 ρ_xy 過(guò)零場(chǎng)是否和實(shí)驗(yàn)測(cè)到的過(guò)零場(chǎng)一致。如果不一致說(shuō)明 Rxx 擬合的參數(shù)有問(wèn)題需要把 ρ_xy 數(shù)據(jù)也納入聯(lián)合擬合。聯(lián)合擬合的代碼結(jié)構(gòu)是在目標(biāo)函數(shù)里同時(shí)計(jì)算 Rxx 和 ρ_xy 的殘差然后加權(quán)求和。權(quán)重取各自測(cè)量誤差的倒數(shù)平方。這樣得到的參數(shù)比單獨(dú)擬合 Rxx 更可靠尤其是當(dāng) Rxx 對(duì)某些參數(shù)不敏感時(shí)ρ_xy 可以提供額外約束。最后說(shuō)一個(gè)我自己的習(xí)慣每次擬合完把參數(shù)、誤差、殘差、擬合區(qū)間、溫度都記在一個(gè)表里存成 CSV。做變溫系列的時(shí)候把表畫(huà)成圖一眼就能看出哪個(gè)溫度點(diǎn)的擬合出了問(wèn)題。這個(gè)習(xí)慣幫我省了很多回頭重測(cè)的時(shí)間。雙載流子擬合不是一錘子買(mǎi)賣(mài)數(shù)據(jù)質(zhì)量、初值選擇、約束條件、交叉驗(yàn)證每一步都影響最終結(jié)果的可信度。希望幫到你。本文還有配套的精品資源點(diǎn)擊獲取