从费米能级到理查逊直线：热电子发射逸出功的两次拟合与避坑指南

简介这份资源是华南理工大学国际校区物理实验课程的金属逸出功实验报告面向修读大学物理实验的理工科学生尤其适合需要参考规范报告结构、理解查逊直线法测钨逸出功原理与数据处理的学习者。压缩包内为1个docx文档约1.11MB内容涵盖实验目的、仪器说明、费米-狄拉克分布与理查逊-热西曼公式推导、零场电流测量及肖脱基效应修正并附有完整的实验步骤与数据处理思路可直接对照实验仪器操作与报告撰写。已有572人学习下载适合作为实验预习、报告模板参考与考前复盘的实用材料。1. 金属逸出功为什么一根烧红的钨丝能“漏”出电子把一根钨丝通电加热到一千多摄氏度它周围的真空里会凭空多出一股电流——没有外加电源推动电子自己就“漏”出来了。这个现象叫热电子发射而决定它强弱的那个材料常数就是逸出功。很多人第一次做这个实验时以为难点在接线和读数真正让人翻车的是后面那两步拟合先要从一组阳极电压里外推出“零场电流”再把五个温度点压成一条直线去抠斜率。这两步背后各藏着一个物理机制一个是空间电荷一个是肖脱基效应任何一个理解偏了算出来的逸出功就会偏出好几电子伏。这份实验报告把理查逊直线法的完整链路——从费米-狄拉克分布到两次取对数拟合——都摊开了值得逐层拆开看它为什么能 work。2. 从费米能级到理查逊公式逸出功到底在量什么2.1 费米分布给出的能量账本金属里的传导电子不是随便什么能量都能占它们按费米-狄拉克分布排布。绝对零度时电子从最低能级一路填到费米能级 EF 为止再往上全是空的。这时候能量最高的那批电子能量恰好等于 EF。但“最高”不等于“能跑出去”——金属表面对电子来说是一道势垒势垒顶端在真空能级 Evac 处两者之差就是逸出功W Evac - EF 单位 eV数值上等于逸出电位 φ这个式子是整份报告的物理地基。它说明逸出功不是某个电子的能量而是一个“门槛高度”处于费米能级的电子至少还要再吸收 W 这么多能量才能越过势垒跑到真空里。温度升高时分布曲线在 EF 上方拖出一条指数衰减的尾巴能量大于势垒的电子数从“几乎为零”变成“可观测”这就是热电子发射的来源。我一般会这么理解逸出功量的是材料把电子“关”在里面的强度数值越大同样温度下漏出来的电子越少。2.2 理查逊-热西曼公式的结构把上面这套统计物理的结论落到宏观可测量上就得到理查逊-热西曼公式I A·S·T²·exp(-eφ / (k·T))其中 A 是与阴极表面化学纯度有关的系数S 是阴极有效发射面积T 是阴极绝对温度k 是玻尔兹曼常数eφ 就是逸出功。这个公式最值得盯住的是那个 T² 因子和指数项的组合温度既通过 T² 线性放大发射又通过 exp(-eφ/kT) 指数级放大发射后者才是主导。也就是说逸出功 eφ 出现在指数的分子上它对电流的影响是“决定性”的——eφ 稍微估错一点反推出来的温度或电流就会差出一大截。这正是为什么实验必须用拟合而不是单点计算单点测量把 T² 和指数项混在一起根本分不开。2.3 取两次对数把曲线压成直线直接拿 I 对 T 做图是一条陡峭的指数曲线肉眼和最小二乘都不好处理。报告里的做法是两边除以 T² 再取对数lg(I / T²) lg(A·S) - (eφ / (2.303·k)) · (1/T)这一步是整份报告的方法论核心。取对数之后lg(I/T²) 与 1/T 变成线性关系斜率里只含 eφ截距里只含 A·S。于是“测逸出功”被转化成一个纯粹的直线拟合问题横轴是 1/T纵轴是 lg(I/T²)斜率 k_fit 一旦拟合出来逸出功就能直接反解。报告里给出的换算式是 Ф 5039·kk 取斜率绝对值这个 5039 就是把 2.303·k 的倒数换算成以 eV 为单位后的常数。理解了这个推导后面所有操作——为什么要测多个温度、为什么要外推零场电流——就都有了目的一切是为了让这条直线成立且干净。3. 零场电流的测量空间电荷与肖脱基效应怎么互相打架3.1 空间电荷为什么让电流测不准热电子从阴极飞向阳极的过程中会在两极之间堆积成一团空间电荷。这团负电荷产生的电场方向与外加加速场相反会阻碍后续电子飞向阳极结果就是你测到的阳极电流 Is 并不等于阴极真正发射的电流而是被空间电荷“压”下来的值。温度越高、发射越强空间电荷越厚压制越明显。如果直接拿这个被压制的电流去拟合等于把两个效应混在一起斜率必然失真。所以实验的第一步不是测 I而是想办法把空间电荷的影响剥掉——这就是为什么要在阳极和阴极之间加一个加速场 Ua。3.2 肖脱基效应加速场带来的“副作用”加了加速场 Ua 之后电子确实能被迅速拉向阳极空间电荷被压下去了。但 Ua 同时会降低阴极表面的势垒高度让电子更容易逸出这就是肖脱基效应。它意味着你测到的电流 Ia 是“在加速场作用下的发射电流”比真正的零场电流 I 要大。两个效应方向相反空间电荷让电流偏小肖脱基效应让电流偏大。实验设计的巧妙之处在于它不去单独消除某一个而是利用肖脱基效应与 Ua 的确定关系把零场电流外推出来。3.3 用 Ua^1/2 外推到 Ua0对共轴圆柱形电极忽略接触电位差阴极表面加速场 Ea 与阳极电压 Ua 的关系可以写成 Ea ∝ Ua^1/2 / (r1·ln(r2/r1))其中 r1、r2 分别是阴极和阳极半径。把这个关系代入肖脱基效应的对数形式得到lg(Ia) lg(I) C · Ua^1/2C 是一个与电极尺寸和温度有关的常数。这个式子的含义非常直接在固定温度 T 下lg(Ia) 对 Ua^1/2 作图是一条直线而这条直线在 Ua0 处的截距就是 lg(I)——也就是零场电流的对数。报告里每个灯丝电流下测了 16V 到 121V 共七组阳极电压正是为了画出这七条直线、各自外推出截距。这一步是整个实验里最容易踩坑的地方如果 Ua 的取值范围太窄直线斜率估不准截距外推就会飘如果某个点读数抖动整条线的截距都会被带偏。3.4 温度从哪来灯丝电流与 T 的标定外推出零场电流之后还需要每个电流对应的阴极温度 T。报告用的是标准二极管的经验标定式T 900 1430 · If If 单位 AT 单位 K这个线性关系是厂家针对特定管型标定出来的不是普适公式。它意味着温度不是直接测的而是通过灯丝电流 If 换算的。这里有一个容易被忽略的细节不同电子管即使灯丝电流相同灯丝温度也不完全一致因为生产工艺有差异。报告在思考题里专门提醒了这一点并建议用多个电子管取平均来减小误差。换句话说T 的误差会直接进入 1/T 横轴而 1/T 在高温端被压缩得很厉害温度标定的小偏差在拟合斜率上会被放大。4. 两次拟合的源码级拆解从原始表格到 4.1627V4.1 第一次拟合每个温度外推零场电流报告给出的原始数据是一张 If × Ua 的电流表五个灯丝电流0.600A 到 0.700A步长 0.025A每个温度下七个阳极电压。第一次拟合的任务是对每一行做 lg(Ia) 对 Ua^1/2 的线性回归取截距。用 Python 复现这个过程的骨架如下import numpy as np # Ua 的七个取值单位 V Ua np.array([16, 25, 36, 49, 64, 81, 100]) sqrt_Ua np.sqrt(Ua) # 报告表格中的一行示例If 0.600A 时的阳极电流mA Ia_0600 np.array([0.050, 0.053, 0.055, 0.058, 0.061, 0.062, 0.062]) # 对 lg(Ia) 与 sqrt(Ua) 做一次线性拟合 coeff np.polyfit(sqrt_Ua, np.log10(Ia_0600), 1) slope, intercept coeff[0], coeff[1] # 截距即 lg(I)反对数得到零场电流 ImA I_zero 10 ** intercept print(fslope{slope:.4f}, lgI{intercept:.4f}, I{I_zero:.4f} mA)这段代码里np.polyfit返回的是一次多项式系数intercept对应 Ua0 处的 lg(I)。报告表格里已经给出了五个温度下的 lgIs 值-1.3586 到 -0.3945说明第一次拟合的截距是逐行算出来的。这里的关键参数是sqrt_Ua——必须对 Ua 开方后再拟合直接拿 Ua 做横轴会得到一条弯曲线截距完全不对。另一个细节是电流单位表格里是 mA取对数时单位不影响斜率但会影响截距的数值反解 I 时要和后续公式的单位保持一致。4.2 第二次拟合斜率里抠出逸出功拿到五个温度下的 lgIs 之后第二次拟合是 lg(Is/T²) 对 1/T 做线性回归。报告表格里已经算好了 T、1/T 和 lg(Is/T²) 三列直接拟合即可# 五个灯丝电流对应的温度K由 T 900 1430*If 算出 T np.array([1758, 1793.75, 1829.5, 1865.25, 1901]) # 第一次拟合得到的 lgIsmA lgIs np.array([-1.3586, -1.0834, -0.8414, -0.6001, -0.3945]) # 纵轴lg(Is/T^2)横轴1/T y lgIs - 2 * np.log10(T) x 1.0 / T coeff2 np.polyfit(x, y, 1) k_fit coeff2[0] # 斜率负值 phi 5039 * abs(k_fit) # 逸出电位单位 V print(fslope{k_fit:.4f}, phi{phi:.4f} V)注意y lgIs - 2*np.log10(T)这一行它等价于 lg(Is) - lg(T²) lg(Is/T²)是把 T² 因子从电流里除掉。这一步不能省否则拟合出来的斜率混入了 T² 的贡献反解出的逸出功会系统性偏大。报告最终得到 Ф 4.1627V对应逸出功 4.103eV因为以 eV 为单位的逸出功数值上等于以 V 为单位的逸出电位。这个结果与钨的理论值相比相对偏差报告写的是 0.0831也就是 8.31%——报告正文里“83.1%”的写法是笔误从数值看应是 8.31%。4.3 数据链路一览把两次拟合串起来整条数据流是这样的阶段输入操作输出原始测量If、Ua、Ia逐点记录5×7 电流表第一次拟合lg(Ia) vs Ua^1/2线性回归取截距5 个 lgIs温度换算IfT 900 1430·If5 个 T第二次拟合lg(Is/T²) vs 1/T线性回归取斜率斜率 k_fit反解k_fitФ 5039·k_fit逸出电位与逸出功这张表说明一件事逸出功不是“测”出来的是“拟合”出来的。任何一个环节的单位、开方、取对数出错最终结果都会偏。报告把中间量 lgIs 和 lg(Is/T²) 都列在表里其实是在给读者留一条可复核的链路——你可以拿原始电流自己重算一遍看能不能复现 4.1627V。5. 避坑与排查这份实验最容易翻车的五个地方5.1 现象拟合出的逸出功偏离理论值好几伏原因横轴用了 Ua 而不是 Ua^1/2。肖脱基效应的推导里加速场 Ea 正比于 Ua^1/2所以外推零场电流时必须对 Ua 开方。直接拿 Ua 做横轴lg(Ia) 与 Ua 不是直线强行拟合得到的截距没有物理意义后续所有温度点的 lgIs 都跟着错。解决画图前先对阳极电压列做一次开方确认 lg(Ia) 对 sqrt(Ua) 的散点大致落在一条直线上再拟合。如果散点明显弯曲先检查是不是某个电压点读数异常。5.2 现象第二次拟合的直线斜率不稳定换一组数据就变原因1/T 在高温端被压缩。五个温度从 1758K 到 1901K1/T 只从 0.0005688 变到 0.0005260跨度不到 8%。横轴跨度这么小纵轴的微小读数误差在斜率上会被放大。这是热电子发射实验的固有难点不是操作失误。解决尽量保证每个温度点的电流读数精度尤其是高温端有条件时增加温度点数量或者用多个电子管的数据合并拟合用样本量换斜率稳定性。5.3 现象零场电流外推出来是负值或明显偏大原因Ua 的取值范围选得不对。如果最低电压太高直线段太短截距外推距离远误差大如果某个高压点已经进入非肖脱基区比如空间电荷完全消失后电流饱和该点会把斜率拉偏。解决先画散点看线性区确认七个点大致共线再拟合。报告选的 16V 到 121V 是覆盖了线性段的换设备时要重新确认这个区间。5.4 现象温度换算出来的 T 和实际灯丝温度对不上原因T 900 1430·If 是针对特定管型的经验标定式不是普适物理定律。不同批次、不同型号的电子管同样的 If 对应的 T 不同。报告思考题里明确提醒了“不同电子管相同灯丝电流灯丝温度不相同”。解决不要跨管型套用这个式子如果换了电子管要么用厂家给的标定曲线要么用多个管子取平均来压低单管偏差。5.5 现象结果偏差 8% 左右怀疑自己算错了原因8% 量级的偏差在这个实验里是正常范围主要来自预热时间不足、温度标定误差和读数精度。报告结论里把“预热时间可能过短”列为误差来源之一是有道理的——阴极温度没稳定费米分布没达到稳态发射电流就不可靠。解决严格预热二十分钟以上等灯丝电流和阳极电流读数都稳定再记录数据处理时保留中间量方便定位是第一次拟合还是第二次拟合引入的偏差。提示实验过程中机箱长期处于高压状态温度较高数据采集结束后要及时降压或关闭仪器接线和拔线时不要触碰金属部分。6. 设计取舍为什么这个实验用“两次外推”而不是直接测回过头看这个实验最值得琢磨的设计取舍是它明明可以直接测某一温度下的发射电流却非要绕两道弯——先用 Ua^1/2 外推零场电流再用 1/T 拟合斜率。这不是为了把实验做复杂而是因为直接测量根本分离不出逸出功。单点测量里电流同时受温度、空间电荷、肖脱基效应三重影响你无法判断测到的值里有多少是材料属性、有多少是测量条件。两次外推的本质是把“测量条件”当作可控变量逐步剥离第一次外推固定温度、消掉加速场的影响得到纯粹的零场发射第二次拟合固定材料、消掉温度的影响只留下斜率里的 eφ。每一步只让一个变量变化这是实验物理里最朴素也最可靠的方法论。这套方法也有它的代价。两次拟合意味着误差会累积第一次外推的截距误差会原封不动进入第二次拟合的纵轴而 1/T 横轴跨度又小斜率对误差极其敏感。所以 8% 左右的偏差不是失败而是这个方法在当前测量精度下的合理预期。想进一步压低偏差方向不是“算得更细”而是“测得更多”——增加温度点、多管取平均、延长预热让温度真正稳态。我做完这类多级拟合的实验后养成了一个习惯先把每一级的中间量单独存一份画一遍散点确认线性再往下走。因为一旦最后结果不对你唯一能定位问题的线索就是这些中间量。希望这份拆解能帮你在下一次面对理查逊直线时知道每一步为什么要这么走。本文还有配套的精品资源点击获取