基于Matlab的夫兰克-赫兹实验的计算机模拟.docxVIP

PAGE

1-

基于Matlab的夫兰克-赫兹实验的计算机模拟

一、引言

(1)夫兰克-赫兹实验是量子力学的重要实验之一，该实验由德国物理学家詹姆斯·夫兰克和戈特哈德·赫兹在1914年进行。实验的主要目的是为了验证原子能级量子化的理论，即原子的能量状态是离散的，而不是连续的。实验过程中，通过向汞蒸气中发射电子束，测量电子能量损失，从而确定汞原子的能级。这一实验为量子力学的发展奠定了坚实的基础，并为后续的研究提供了宝贵的实验依据。

(2)在夫兰克-赫兹实验中，电子束在穿过汞原子时，可能会与汞原子发生碰撞。根据量子力学的原理，当电子的能量等于汞原子的某个能级差时，电子与汞原子会发生弹性或非弹性碰撞，导致电子的能量降低。如果电子的能量低于能级差，则电子不会与汞原子发生作用，电子能量损失为零。而当电子能量恰好等于能级差时，会发生非弹性碰撞，此时电子将损失与其能量相等的能量给汞原子，汞原子从基态跃迁到激发态。

(3)夫兰克-赫兹实验的具体过程如下：将汞蒸气充满一个低压的真空管，然后在真空管两端加上直流高压电源。在真空管的一端，产生电子束，当电子束穿过汞蒸气时，会受到汞原子的阻碍。通过改变电子束的加速电压，可以调节电子的能量。当电子能量恰好等于汞原子的能级差时，电子与汞原子发生非弹性碰撞，电子的能量损失可以通过测量电子束的能量变化来确定。实验结果显示，当加速电压在4.9伏特左右时，电子束的能量损失最为明显，这对应着汞原子从基态跃迁到第一激发态的能量。这一结果与当时的理论预测相吻合，从而证实了量子力学的能级量子化假设。

二、夫兰克-赫兹实验原理及模型建立

(1)夫兰克-赫兹实验的原理基于能量量子化的概念，即原子的能量状态是离散的。在实验中，电子通过加速电压获得动能，当电子能量达到特定值时，与汞原子发生非弹性碰撞，导致汞原子从基态跃迁到激发态。这一跃迁需要特定的能量，称为激发能，对于汞原子来说，这个能量约为4.9电子伏特。实验通过测量电子束的能量损失来验证这一跃迁的发生。

(2)实验模型建立的核心是计算电子与汞原子碰撞时的能量转移。根据量子力学原理，碰撞过程中电子的动能部分或全部转化为汞原子的激发能。模型中假设碰撞是完全非弹性的，即电子失去的能量等于汞原子的激发能。通过实验测量不同电压下电子束的电流，可以构建能量损失与电子能量的关系曲线，从而验证能量量子化的假设。

(3)在实验模型中，电子的能量由电子加速电压决定，电子的动能与电压成正比。通过调节加速电压，可以改变电子的能量，从而实现与汞原子能级相匹配的碰撞。实验数据表明，在电压为4.9伏特时，电子束的电流出现显著下降，这对应于汞原子吸收了足够的能量发生激发跃迁。此外，实验还展示了随着电压的进一步增加，电流不再增加，表明所有的汞原子都已经被激发，验证了能量量子化的理论预测。

三、Matlab编程实现

(1)在Matlab中实现夫兰克-赫兹实验的计算机模拟，首先需要建立电子能量与加速电压之间的关系。根据经典电动力学，电子在电场中的加速电压V下，其动能可以表示为\(E_k=eV\)，其中\(e\)是电子电荷。在模拟中，我们通过改变电压V来计算不同能量水平的电子与汞原子碰撞时的行为。

为了模拟这一过程，我们首先定义了一个电子能量数组，从0电子伏特开始，以0.1电子伏特的步长增加，直到一个预定的上限，例如10电子伏特。接着，我们使用一个循环来模拟每个能量水平的电子与汞原子的碰撞。在循环中，我们计算电子在碰撞后可能获得的能量，以及汞原子可能吸收的能量。如果电子的能量等于汞原子的激发能（4.9电子伏特），则模拟该电子与汞原子发生了激发跃迁。

例如，假设我们使用了一个包含1000个能量点的数组，那么循环中的代码可能如下所示：

```matlab

e=0:0.1:10;%电子能量数组，单位为电子伏特

current=zeros(1,length(e));%存储不同能量下的电流

excitation_energy=4.9;%汞原子的激发能，单位为电子伏特

fori=1:length(e)

%计算电子在碰撞后的能量

energy_after_collision=e(i)-excitation_energy;

%如果能量非负，则表示发生了激发跃迁

ifenergy_after_collision=0

current(i)=current(i)+1;

end

end

```

(2)在Matlab中，我们还可以通过绘制能量与电流的关系图来直观地展示模拟结果。使用绘图函数`plot`，我们可以将电子能量作为横坐标，对应的电流作为纵坐标，生成一个曲线图。例如，以下代码展示了如何生成并显示这个曲线图：

```matlab

plot(e,current);

xlabel(电子能量(eV));

ylabel(电流)