验证快速电子的动量与动能的相对论关系

实验四 验证快速电子的动量与动能的相对论关系【实验目的】1、 验证通过对快速电子的动量及动能的同时测定验证动量和动能之间的相对论关系；2、 了解β磁谱仪测量原理、闪烁记数器的使用方法及一些实验数据处理的思想方法。

【实验原理】经典力学总结了低速物理的运动规律，它反映了牛顿的绝对时空观：认为时间和空间是两个独立的观念，彼此之间没有联系；同一物体在不同惯性参照系中观察到的运动学量(如坐标、速度)可通过伽利略变换而互相联系。

这就是力学相对性原理：一切力学规律在伽利略变换下是不变的。

19世纪末至20世纪初，人们试图将伽利略变换和力学相对性原理推广到电磁学和光学时遇到了困难；实验证明对高速运动的物体伽利略变换是不正确的，实验还证明在所有惯性参照系中光在真空中的传播速度为同一常数。

在此基础上，爱因斯坦于1905年提出了狭义相对论；并据此导出从一个惯性系到另一惯性系的变换方程即“洛伦兹变换”。

洛伦兹变换下，静止质量为m 0，速度为v 的物体，狭义相对论定义的动量p 为：p m v mv =−=012β (1) 式中m m v c =−=012/,/ββ。

相对论的能量E 为：E mc =2 (2)这就是著名的质能关系。

mc 2是运动物体的总能量，当物体静止时v=0，物体的能量为E 0=m 0c 2称为静止能量；两者之差为物体的动能E k ，即E mc m c m c k =−=−−222200111()β (3)当β« 1时，式（3）可展开为 E m c v c m c m v p m k =++−≈=00022222201121212()L (4) 即得经典力学中的动量—能量关系。

由式(1)和(2)可得：E c p E 22202−= (5)这就是狭义相对论的动量与能量关系。

而动能与动量的关系为：E E E c p m c m c k =−=+−02242020 (6)这就是我们要验证的狭义相对论的动量与动能的关系。

对高速电子其关系如图所示，图中pc 用MeV 作单位，电子的m 0c 2=0.511MeV 。

式(4)可化为：E p c m c p c k ==×1220511222220.以利于计算。

【实验装置及方法】实验装置主要由以下部分组成：①真空、非真空半圆聚焦β磁谱仪；②β放射源90Sr—90Y(强度≈1毫居里)，定标用γ放射源137Cs和60Co(强度≈2微居里)；③200μmAl窗NaI(Tl)闪烁探头；④数据处理计算软件；⑤高压电源、放大器、多道脉冲幅度分析器。

β源射出的高速β粒子经准直后垂直射入一均匀磁场中(BV⊥)，粒子因受到与运动方向垂直的洛伦兹力的作用而作圆周运动。

如果不考虑其在空气中的能量损失(一般情况下为小量)，则粒子具有恒定的动量数值而仅仅是方向不断变化。

粒子作圆周运动的方程为：dp dtev B=−×(7)e为电子电荷，v为粒子速度，B为磁场强度。

由式(1)可知p=mv，对某一确定的动量数值P，其运动速率为一常数，所以质量m是不变的，故dp dt m dv dt =, 且dv dt v R=2所以 p eBR = (8)式中R 为β粒子轨道的半径，为源与探测器间距的一半。

在磁场外距β源X 处放置一个β能量探测器来接收从该处出射的β粒子，则这些粒子的能量(即动能)即可由探测器直接测出，而粒子的动量值即为：p eBR eB X ==Δ/2。

由于β源38903990Sr Y −(0～2.27MeV)射出的β粒子具有连续的能量分布(0～2.27MeV)，因此探测器在不同位置(不同ΔX)就可测得一系列不同的能量与对应的动量值。

这样就可以用实验方法确定测量范围内动能与动量的对应关系，进而验证相对论给出的这一关系的理论公式的正确性。

【实验内容】1． 测量快速电子的动量。

2． 测量快速电子的动能。

3． 验证快速电子的动量与动能之间的关系符合相对论效应。

【实验步骤】1． 检查仪器线路连接是否正确，然后开启高压电源，开始工作；2． 打开60Co γ定标源的盖子，移动闪烁探测器使其狭缝对准60Co 源的出射孔并开始记数测量；3． 调整加到闪烁探测器上的高压和放大数值，使测得的60Co 的1.33MeV 峰位道数在一个比较合理的位置（建议：在多道脉冲分析器总道数的50%～70%之间，这样既可以保证测量高能β粒子(1.8~1.9MeV)时不越出量程范围，又充分利用多道分析器的有效探测范围）；4． 选择好高压和放大数值后，稳定10～20分钟；5． 正式开始对NaI(Tl)闪烁探测器进行能量定标，首先测量60Co 的γ能谱，等1.33MeV 光电峰的峰顶记数达到1000以上后（尽量减少统计涨落带来的误差），对能谱进行数据分析，记录下1.17和1.33MeV 两个光电峰在多道能谱分析器上对应的道数CH 3、CH 4；6． 移开探测器，关上60Co γ定标源的盖子，然后打开137Cs γ定标源的盖子并移动闪烁探测器使其狭缝对准137Cs 源的出射孔并开始记数测量，等0.661MeV 光电峰的峰顶记数达到1000后对能谱进行数据分析，记录下0.184MeV 反散射峰和0.661 MeV 光电峰在多道能谱分析器上对应的道数CH 1、CH 2；7． 关上137Cs γ定标源，打开机械泵抽真空（机械泵正常运转2～3分钟即可停止工作）；8． 盖上有机玻璃罩，打开β源的盖子开始测量快速电子的动量和动能，探测器与β源的距离ΔX 最近要大于9cm 、最远要小于24cm ，保证获得动能范围0.4～1.8MeV 的电子；9． 选定探测器位置后开始逐个测量单能电子能峰，记下峰位道数CH 和相应的位置坐标X ；10．全部数据测量完毕后关闭β源及仪器电源，进行数据处理和计算。

【数据处理】1.真空状态下P 与ΔX 的关系的合理表述由于工艺水平的限制，磁场的非均匀性(尤其是边缘部分)无法避免，直接用p eBR eB X ==Δ/2来求动量将产生一定的系统误差；因此需要采取更为合理的方式来表述P 与ΔX 的关系。

设粒子的真实径迹为aob ，位移ds 与Y 轴的夹角为θ，如上图所示；则ds 在X 轴上的投影为sin θ⋅ds 。

显然有：Δx ds ds =⋅≈⋅∫∫sin sin θθπθ001 ()Q θπ1≈ ( 9) 又因为ds R d =⋅θ以及eB P R /=，(其中R 、B 分别为ds 处的曲率半径和磁场强度)，则有：ΔX P eB d P e Bd =⋅⋅=⋅∫∫sin sin θθθθππ00 (Q 真空中P 为定值) (10) 所以有： P e X B d Be X =⋅=∫ΔΔ/sin θθπ120 (1120B B d =⋅∫sin θθπ) (11) 把1B 改写成： ∫∫⋅⋅=ππθθθθ00sin /sin 1d d B B， 则物理含义更为明显：即B /1为粒子在整个路径上的磁场强度的倒数以各自所处位置处的位移与Y 轴夹角的正弦为权重的加权平均值。

显然，B 相当于均匀磁场下公式p eBR eB X ==Δ/2中的磁场强度B ；即只要求出B ，就能更为确切地表述P 与ΔX 的关系，进而准确地确定粒子的动量值。

实际计算操作中还需要把求积分进一步简化为求级数和；即可把画在磁场分布图上直径为ΔX 的半圆弧作N 等分(间距取10毫米左右为宜)，依此读出第i 段位移所在处的磁场强度B i ，再注意到：θπi N i =−()1以及Δθπi N =，则最后求和可以得到： 11212121101B B d N N i B N N i B i N i i i N ≈′′⋅′≈−=−==∑∫∑sin sin[()]/sin[()]/θθπππππ (12)所以：P Ne xN i B ii N =−=∑Δππsin[()]/11 (13) 2.β粒子动能的测量β粒子与物质相互作用是一个很复杂的问题，如何对其损失的能量进行必要的修正十分重要。

①β粒子在Al 膜中的能量损失修正在计算β粒子动能时还需要对粒子穿过Al 膜(220μm ：200μm 为NaI(Tl)晶体的铝膜密封层厚度，20μm 为反射层的铝膜厚度)时的动能予以修正，计算方法如下。

设β−粒子在Al 膜中穿越Δx 的动能损失为ΔE ，则：ΔΔE d E d x x =ρρ (14) 其中dE dx ρ(dE dx ρ<0)是Al 对β−粒子的能量吸收系数，(ρ是Al 的密度)，dE dx ρ是关于E 的函数，不同E 情况下dE dx ρ的取值可以通过计算得到。

可设dE dx K E ρρ=()，则ΔE=K(E)Δx ；取Δx →0，则β−粒子穿过整个Al 膜的能量损失为：∫+=−d x x dx E K E E )(12 (15)； 即∫+−=dx xdx E K E E )(21 (16) 其中d 为薄膜的厚度，E 2为出射后的动能，E 1为入射前的动能。

由于实验探测到的是经Al 膜衰减后的动能，所以经公式(4─9)可计算出修正后的动能(即入射前的动能)。

下表列出了根据本计算程序求出的入射动能E 1和出射动能E 2之间的对应关系：②β粒子在有机塑料薄膜中的能量损失修正此外，实验表明封装真空室的有机塑料薄膜对β存在一定的能量吸收，尤其对小于0.4MeV的β粒子吸收近0.02MeV 。

由于塑料薄膜的厚度及物质组分难以测量，可采用实验的方法进行修正。

实验测量了不同能量下入射动能E k 和出射动能E 0(单位均为MeV)的关系，采用分段插值的方法进行计算。

具体数据见下表： E k (MeV) 0.382 0.581 0.777 0.973 1.173 1.367 1.567 1.752 E 0(MeV) 0.365 0.571 0.770 0.966 1.166 1.360 1.557 1.7473．数据处理的计算方法和步骤（举例说明）：设对探测器进行能量定标（操作步骤中的第5、6步）的数据如下：能量(MeV)0.184 0.661 1.17 1.33 道数(CH) 48 152 262 296实验测得当探测器位于21cm 时的单能电子能峰道数为204，求该点所得β粒子的动能、动量及误差，已知β源位置坐标为6cm 、该点的有效磁场强度为620高斯(Gs)。

1） 根据能量定标数据求定标曲线已知；48,184.011==CH MeV E ；152,661.022==CH MeV E ；262,17.133==CH MeV E ；296,33.144==CH MeV E ；根据最小二乘原理用线性拟合的方法求能量E 和道数CH 之间的关系：E a b C H =+×可以推导，其中：])([12∑∑∑∑⋅⋅−⋅Δ=ii i i i i i i i E CH CH E CH a ])([1∑∑∑⋅−⋅Δ=ii i i i i i E CH E CH n b Δ=−∑∑n C H C H i i i i22()代入上述公式计算可得：CH E ⋅+−=0046.0038613.02）求β粒子动能对于X=21cm 处的β粒子：① 将其道数204代入求得的定标曲线，得动能E 2=0.8998MeV ，注意：此为β粒子穿过总计220μm 厚铝膜后的出射动能，需要进行能量修正；② 在前面所给出的穿过铝膜前后的入射动能E 1和出射动能E 2之间的对应关系数据表中取E 2=0.8998MeV 前后两点作线形插值，求出对应于出射动能E 2=0.8998MeV 的入射动能E 1=0.9486MeVE 1(MeV) E 2(MeV)0.937 0.8500.988 0.900③ 上一步求得的E 1为β粒子穿过封装真空室的有机塑料薄膜后的出射动能E 0，需要再次进行能量修正求出之前的入射动能E k ，同上面一步，取E 0=0.9486MeV 前后两点作线形插值，求出对应于出射动能E 0=0.9486MeV 的入射动能E k =0.9556MeV ； E k (MeV) 0.777 0.973E 0(MeV) 0.770 0.966E k =.09556MeV 才是最后求得的β粒子动能。