![](/user_photo/_userpic.png)
Елютин Експерименталная физика.Лабораторный практикум 2011
.pdf![](/html/65386/144/html_spNKoIRX_y.1NDB/htmlconvd-8wNt7n31x1.jpg)
Схема установки изображена на рис. 34.8. Излучение источника 1 – галогенной лампы накаливания, проходит через систему нейтральных фильтров и пленку, укрепленных в держателях, и через входную щель попадает в монохроматор 2. Из входной щели монохроматора излучение направляется в блок приемника 3, основой которого служит фотоэлемент. Выходной сигнал с фотоэлемента после усиления через кабель и разъем подается на гнездо блока питания 4. Источник излучения 1 также подключается к блоку питания 4. Регистрация интенсивности излучения в относительных единицах осуществляется на шкале цифрового вольтметра 5. Вход цифрового вольтметра соединен с гнездами блока питания 4. Необходимые длины волн устанавливаются рукояткой 2а монохроматора и отсчитываются по цифровому механическому счетчику с точностью 0,2 нм.
Рис. 34.8
Чувствительность установки регулируется диафрагмой 3а и ступенчатым переключателем 3б. Рукоятка 3в устанавливает «нуль» вольтметра. На лицевой панели блока питания 4 имеется тумблер 4а – «сеть» и сигнальная лампа, тумблер питания фотоэлемента 4б «накал-недокал», обеспечивающий работу источника излучения в
31
![](/html/65386/144/html_spNKoIRX_y.1NDB/htmlconvd-8wNt7n32x1.jpg)
двух режимах, и гнезда для подключения вольтметра 4в. Рабочий предел вольтметра – 20 В.
Снятие спектров пропускания осуществляется в режиме «недокал» источника излучения 1 в области 450-630 нм, которая является областью наибольшей чувствительности установки. В этом диапазоне длин волн используется входная щель монохроматора шириной 0,25 мм, а выходная щель шириной 0,05 мм. Во избежание перегрузки вольтметра, световой поток ослабляется с помощью двух шестикратных нейтральных фильтров. В этом случае положение ступенчатого переключателя чувствительности 3б – «3».
До начала измерений необходимо установить «0» на табло вольтметра с помощью рукоятки 3в, предварительно перекрыв доступ света в монохроматор. Следует помнить, что с вольтметра производится отсчет установившегося значения напряжения (для чего требуется выждать 10 – 15 секунд).
Измерения производят вначале без пленки с тем, чтобы снять аппаратную функцию спектральной чувствительности фотоэлемента U0 ( ). При перегрузке вольтметра из-за большой величины све-
тового потока необходимо ввести один, а при необходимости и два светофильтра. Коэффициент ослабления света k указан на самом светофильтре. Светофильтры не убираются до конца измерения.
После этого производят измерение спектра пленки As2S3 толщиной 1,8 мкм U ( ) в том же диапазоне длин волн. На коротковолновом участке диапазона показания вольтметра могут стать малыми, что потребует снятия фильтра, с тем, чтобы восстановить значащие цифры в разрядах дисплея вольтметра.
Коэффициент пропускания вычисляется по формуле |
|
D( ) U ( ) U0 ( ). |
(34.6) |
Если в ходе выполнения работы снимался фильтр, то необходимо вводить поправочный коэффициент, т.е. использовать формулу
D( ) kU ( ) U0 ( ). |
(34.7) |
Порядок проведения эксперимента
1. Включите блок питания монохроматора, вольтметр и источник света. Дайте приборам прогреться не менее 15 минут.
32
![](/html/65386/144/html_spNKoIRX_y.1NDB/htmlconvd-8wNt7n33x1.jpg)
2.Установите на счетчике монохроматора длину волны 600 нм. Произведите установку «0» вольтметра, перекрыв предварительно доступ света в монохроматор.
3.Измерьте без пленки аппаратную функцию U0 ( ) в диапазоне
длин волн 630-450 нм через каждые 10 нм. Результат занесите в таблицу 34.2.
4. Произведите с пленкой измерения спектральной функции U ( ) в том же диапазоне длин волн. Отметьте измерения, где был снят светофильтр. Результат занесите в таблицу 34.2.
Таблица 34.2
, , нм |
U0 ( ) |
U ( ) |
D( ) |
|
|
|
|
630 |
|
|
|
|
|
|
|
620 |
|
|
|
|
|
|
|
|
|
|
|
460 |
|
|
|
|
|
|
|
450 |
|
|
|
|
|
|
|
5.Вычислите D( ) по формулам (34.6) и (34.7).
6.Постройте график D( ) и по графику определите ГР и .
7.Вычислите n по формуле (34.5) и Eg формуле (34.2). Срав-
ните полученный результат с табличным.
Контрольные вопросы
1.Как образуются энергетические зоны в кристалле?
2.Объясните «прозрачность» веществ с точки зрения зонной теории.
3.Объясните физический смысл ГР.
4.Что происходит с запрещенной зоной кристаллического полупроводника при появлении дефектов в его структуре?
33
Литература
1. Яворский Б.М., Пинский А.А. Основы физики в 2-х томах.
Т2. М.: Физматлит. 2003.
2.Савельев И.В. Курс общей физики. Т 3. М.: Наука. 1989. С. 197-210.
3.Физические величины: Справочник. / Под ред. И. С. Григорьева, Е.З. Мейлихова. М.: Энергоатомиздат. 1991.
РАБОТА 35
КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ТРАЕКТОРИЙ ЗАРЯЖЕННЫХ ЧАСТИЦ В ОДНОРОДНЫХ ПОЛЯХ
Цель – исследование траектории движения заряженных частиц в однородных электрическом и магнитном полях с помощью компьютерного моделирования.
Введение
В повседневной научной практике часто встречаются задачи, когда имеется возможность записать исходные уравнения и известны все начальные и граничные условия, но решение которых затруднительно, либо совсем невозможно осуществить аналитическим образом. В этих случаях применяют компьютерное моделирование. В данной работе решаются довольно сложные системы уравнений, решить которые учащимся на данном этапе затруднительно, но с помощью специальных алгоритмов и программы моделирования это становится возможным.
Движение частицы в однородном электрическом поле с учетом сопротивления среды
При моделировании движения положительно заряженного тела (частицы) в однородном электрическом поле (для определенности поле направлено вертикально вниз) сила сопротивления среды, в которой движется частица, была выбрана пропорциональной квадрату скорости частицы. Расчеты показали сильное влияние сопротивления на форму траектории. Модель позволяет определить время движения частицы в заданной области пространства, мгновен-
34
![](/html/65386/144/html_spNKoIRX_y.1NDB/htmlconvd-8wNt7n35x1.jpg)
ную скорость, геометрические параметры траектории, величину диссипации энергии при движении в среде.
Решение этой задачи без учета сил сопротивления хорошо известно. Так, расстояние от точки влета заряженной частицы в область однородного поля до вторичного пересечения той же эквипотенциальной плоскости задается выражением
L |
v2 sin2 |
, |
|
|
||
0 |
|
|
|
|
||
qE m |
|
|||||
|
|
|
|
|||
с максимальным значением |
|
|||||
L |
|
|
v2 |
|
, |
(35.1) |
|
0 |
|
||||
|
|
|
||||
МАКС |
|
qE m |
|
|||
|
|
|
|
|||
где v0 |
– скорость влета частицы, |
– угол влета по отношению к |
эквипотенциальной поверхности. Максимальное смещение частицы против направления поля составит
H v02 sin2 . 2 qEm
Время движения
T |
|
|
2v0 sin |
, |
||||
|
|
|
|
|
||||
ПОЛЕТА |
|
|
|
qE m |
|
|||
|
|
|
|
|
|
|||
время движения при максимальной дальности: |
||||||||
|
|
|
|
|
|
v0 |
. |
|
T |
|
|
|
|
2 |
(35.2) |
||
|
|
|
|
|
||||
МАКС |
|
|
qE m |
|
||||
|
|
|
|
Траектория движения частицы без учета сил сопротивления выглядит следующим образом:
y |
y |
0 |
x tg |
x2 (qE m) |
. |
(35.3) |
||
2v2 cos2 |
|
|||||||
|
|
|
|
|
||||
|
|
|
|
0 |
|
|
|
Уравнение второго закона Ньютона для движения с учетом силы сопротивления задается векторным равенством
ma qE F |
. |
(35.4) |
СOПР |
|
|
Предполагается, что сила сопротивления FСОПР направлена по касательной к траектории против движения тела и может быть выбрана либо пропорциональной вектору скорости тела с обратным знаком, т.е. в виде
35
![](/html/65386/144/html_spNKoIRX_y.1NDB/htmlconvd-8wNt7n36x1.jpg)
FСОПР rЛИНv,
либо пропорциональной квадрату модуля скорости тела в виде
|
|
v |
|
|
F |
r |
v2. |
||
|
||||
СОПР |
КВ v |
В зависимости от выбора, вектор силы сопротивления имеет следующие проекции:
FСОПР rЛИНvx , rЛИНvy ,
или
FСОПР rКВvx v2x v2y , rКВvy
vx2 vy2 , ,
где rЛИН, КВ – соответствующий коэффициент сопротивления среды,
в которой движется тело.
Уравнение (35.4) в проекциях на оси координат приводит к системе дифференциальных уравнений по времени для координат x и y
d 2 x |
|
|
|
r |
|
|
|
|
||
|
|
|
|
ЛИН |
|
v |
x |
|
|
|
dt2 |
|
|
|
|
||||||
|
|
|
m |
|
|
|||||
d 2 y |
|
|
|
r |
|
|
|
qE |
||
|
|
|
|
ЛИН |
v |
y |
|
|
||
dt2 |
|
|
m |
|||||||
|
|
|
m |
|
|
|
|
r |
|
|
|
|
|
|
|
|
vx2 vy2 |
|
|
|
|||||
|
|
КВ |
vx |
|
|
, |
|
|
|
|
|
|
|||||||
|
|
m |
|
|
|
|
(35.5) |
||
|
|
r |
|
|
|
|
qE |
||
|
vx2 v2y |
|
|
|
|||||
|
|
КВ |
vy |
|
|
|
. |
||
|
|
||||||||
|
|
m |
|
|
|
|
m |
Кинетическая энергия частицы равна
K |
mv2 |
(35.6) |
||
|
. |
|||
2 |
||||
|
|
|
||
Потенциальная энергия заряженной частицы: |
|
|||
qE( y y0 ). |
(35.7) |
Для удобства дальнейшего рассмотрения перейдем к безразмерным переменным:
Y |
( y y0 ) |
, X |
(x x0 ) |
, |
t |
, |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
LМАКС |
TМАКС |
|
|
||||||||||||||||||||
|
LМАКС |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
V |
|
|
|
v |
|
|
,V |
|
|
|
|
|
|
|
vy |
|
|
|
|
|
|
|
|
(35.8) |
|||||
2 |
x |
|
y |
|
2 |
, |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
x |
|
|
|
v0 |
|
|
|
|
|
|
|
|
|
|
v0 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
v |
|
|
|
|
|
|
|
|
|
v2 |
|
|
|
H |
1 |
|
||
b |
r |
|
|
|
|
|
2 |
0 |
, b |
r |
, h |
|
sin2 . |
||||||||||||||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
||||||||||||||||
|
|
|
|
qE |
|
|
L |
|
2 |
||||||||||||||||||||
ЛИН |
|
ЛИН |
|
|
|
|
|
КВ |
|
КВ qE |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
МАКС |
|
|
36
![](/html/65386/144/html_spNKoIRX_y.1NDB/htmlconvd-8wNt7n37x1.jpg)
В обозначениях (35.8) уравнение (35.3) для движения без учета силы сопротивления среды примет следующий вид:
Y Y X tg |
X 2 |
|
. |
(35.9) |
2cos2 |
|
|||
0 |
|
|
Уравнения Ньютона (35.5) при учете силы сопротивления среды можно записать в безразмерных переменных (35.8) следующим образом:
dVx |
|
|
|
|
|
|
|
|||
|
|
2 |
2 |
|
|
|||||
|
|
|
bЛИНVx |
bКВVx Vx |
Vy |
, |
||||
d |
||||||||||
dVy |
|
2 bКВVy |
|
|
, |
|||||
2 bЛИНVy |
Vx2 Vy2 |
|||||||||
|
|
|||||||||
d |
|
|
|
|
|
|
|
|
||
dX |
Vx , |
|
|
|
|
(35.10) |
||||
|
|
|
|
|
||||||
d |
|
|
|
|
|
|
|
|
dY
d
Vy .
Выражения для потенциальной и кинетической энергии частицы в обозначениях (35.8) будут выглядеть как
Wp |
|
|
2Y , |
|
|
(35.11) |
||
mv2 |
2 |
|
|
|||||
|
0 |
|
|
Vx2 |
Vy2 |
|
|
|
W |
|
K |
|
|
. |
(35.12) |
||
mv2 |
2 |
|
2 |
|||||
k |
|
|
|
|
||||
|
0 |
|
|
|
|
|
|
Потери энергии на работу против силы сопротивления среды составят
W 1 2Y |
Vx2 Vy2 |
. |
(35.13) |
|
|||
2 |
|
|
Систему дифференциальных уравнений (35.10) представим в виде уравнений в конечных разностях, когда дифференциал функции заменяется разностью ее последовательных значений в узлах (n 1) и (n) расчетной сетки, разделенных временным шагом h. Таким образом, система уравнений в конечных разностях для движения заряженной частицы в однородном электрическом поле будет иметь следующий вид:
37
![](/html/65386/144/html_spNKoIRX_y.1NDB/htmlconvd-8wNt7n38x1.jpg)
|
Vx n 1 |
Vx n |
|
|
|
V n |
|
b V n |
|
|
|
|
|
|
|
||
|
|
b |
|
V n 2 V n 2 |
, |
|
|
||||||||||
h |
|
|
|||||||||||||||
|
|
ЛИН x |
|
КВ x |
|
x |
y |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Vy n 1 |
Vy n |
|
|
|
|
|
|
n |
|
|
|
|
|||||
|
|
|
n |
2 b V |
V n 2 V |
n 2 |
, |
||||||||||
|
|
|
|
2 bЛИНVy |
|
КВ |
y |
x |
y |
|
|
||||||
|
h |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
X n 1 |
X n |
|
Vx n , |
Y n 1 |
Y n |
Vy n , |
|
n 1 n h, |
|
(35.14) |
||||||
|
|
h |
|
|
|
h |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n 2 n 2
W n 1 2Y n Vx Vy . 2
Система уравнений (35.14) решается численно. Для этого первые пять уравнений системы (35.14) разрешаются относительно величин с индексом (n 1). Получившиеся алгебраические уравнения используются как рекуррентные соотношения для получения значений переменных для всех n, т.е. получения численного решения в виде массива чисел. В качестве начальных значений искомых функций задавались следующие величины:
X (0) 0, Y (0) 0, (0) 0,Vx(0) |
2 cos ,Vy(0) |
2 sin . |
Точность такой численной процедуры (схемы Эйлера) составляет o(h).
Сравнение численного решения с табуляцией классической формулы (35.9) показывает влияние силы сопротивления среды. Заметим, что координата X ( ) в уравнении (35.10) получается как решение соответствующего дифференциального уравнения. Эта величина не может быть заменена известным выражением
X ( ) 2 cos из-за наличия ускорения вдоль оси X , как результата действия силы сопротивления. Результат одного из расчетов (b 1, 45 ) представлен на рисунке 35.1. График 1 – без учета сопротивления среды, график 2 – с учетом сопротивления среды, график 3 – потери энергии на сопротивлении среды.
Моделирование пузырьковой камеры
Пузырьковой камерой называется прибор для регистрации следов (треков) заряженных частиц, как правило, в магнитном поле. Действие этого прибора основано на вскипании перегретой жидко-
38
![](/html/65386/144/html_spNKoIRX_y.1NDB/htmlconvd-8wNt7n39x1.jpg)
сти вблизи траектории частицы. Перегрев жидкости осуществляется быстрым понижением давления в камере.
Рис. 35.1
Прохождение заряженной частицы через перегретую жидкость приводит к образованию вдоль траектории частицы «зародышевых» центров кипения. За время порядка 0,5-3 мс пузырьки, образующиеся на зародышах, достигают размеров 50-300 мкм и могут быть сфотографированы при освещении их импульсным источником света.
Во время движения в пузырьковой камере на частицу действуют сила Лоренца q[v B] и сила сопротивления среды, выбранная в на-
шей задаче в виде r v v2 (здесь v – вектор скорости частицы,
КВ v
q – ее заряд, B – индукция магнитного поля, rКВ – коэффициент сопротивления среды). Следовательно, второй закон Ньютона для частицы можно записать в виде
|
|
v |
|
|
|
ma |
q[v B] r |
v2. |
(35.15) |
||
|
|||||
|
КВ v |
|
Будем рассматривать случай, когда частица влетает в камеру со скоростью v0 перпендикулярно к направлению вектора индукции
магнитного поля B. В этом случае движение частицы будет происходить в плоскости, перпендикулярной к векторуB. Движение частицы удобно рассматривать в системе координат, у которой ось Z
39
![](/html/65386/144/html_spNKoIRX_y.1NDB/htmlconvd-8wNt7n40x1.jpg)
совпадает с направлением вектора B, а ось Y составляет с начальным направлением вектора скорости v0 угол (рис. 35.2).
Рис. 35.2
Запишем уравнение (35.13) в проекции на оси такой системы координат:
|
d 2 x |
|
|
dy |
|
|
|
|
|
|
|
|
|
|||
m |
|
qB |
r v |
|
|
|
v2 v2 , |
|
|
|||||||
dt2 |
|
|
|
|
|
|
|
|||||||||
|
|
|
dt |
КВ |
x |
|
x |
y |
|
|
||||||
|
d 2 y |
|
|
dx |
|
|
|
|
|
|
|
|
|
|||
m |
|
qB |
r |
|
v |
|
|
v2 |
v2 |
, |
(35.16) |
|||||
dt2 |
|
dt |
|
|
||||||||||||
|
|
|
|
КВ |
|
y |
|
x |
|
y |
|
|
||||
m |
d 2 z |
|
0. |
|
|
|
|
|
|
|
|
|
|
|||
dt2 |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Начальные условия для задачи (35.16) выбраны следующими: |
||||||||||||||||
x(0) x0 , y(0) y0 , vx |
v0 sin , vy |
v0 cos . |
Перейдя к безразмерным переменным в первых двух уравнениях системы (35.16), можно получить следующую систему:
|
dX |
|
|
dVx |
|
2 Vy bКВVx |
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
||||
|
|
|
|
2 |
|
2 |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
Vx , |
|
|
Vx |
Vy |
|
|
|
|
|
|
|
|
||||||||||
|
d |
d |
|
|
|
(35.17) |
|
||||||||||||||||||
|
dY |
|
|
dVy |
|
2 Vx bКВVy |
|
|
|
|
|
|
|
, |
|
|
|
|
|||||||
|
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|||||||
|
|
Vy , |
|
|
|
|
Vx |
Vy |
|
|
|
|
|
|
|
|
|||||||||
|
d |
|
|
d |
|
|
|
|
|
|
|
vx, y |
|
||||||||||||
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
r v2 |
|
r v2 |
|
|
|
||||
где введены переменные |
|
|
, b |
|
|
|
КВ 0 |
|
КВ 0 |
, V |
x,y |
|
|
|
, |
||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
T |
|
|
|
КВ |
|
2 mv0 T |
|
qBv0 |
|
v0 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
40