Том 170
ИЗВЕСТИЯ
ТОМСКОГО ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ ПОЛИТЕХНИЧЕСКОГО ИНСТИТУТА им. С. М. КИРОВА
1969
ВЛИЯНИЕ КРИСТАЛЛИЧЕСКОЙ СТРУКТУРЫ ВЕЩЕСТВА НА ПОТЕРИ ЭНЕРГИИ БЫСТРЫМИ ПРОТОНАМИ В МОНОКРИСТАЛЛАХ ГАЛОГЕНИДОВ ЩЕЛОЧНЫХ МЕТАЛЛОВ
Э. Т. ШИПАТОВ, Б. А. КОНОНОВ
В работе рассмотрен вопрос о каналировании быстрых протонов в ионных кристаллах, обладающих простой кубической структурой. Рассмотрены случаи прохождения протонов в направлении [100] и между плоскостями (100) монокристаллов. Рассчитаны потенциалы для атомных плоскостей и атомных рядок кристаллов и определены предельные углы каналирования для соединений га • догенидов натрия и калия.
Исследованы энергетические спектры протонов, прошедших монокристаллические мишени различных соединений, и найдена корреляция между формой энергетического спектра и соотношением радиусов ионов в соединении, между долей каналирования частиц и порядковым номером галоида в рядах.
The channelling of fast protons in the simple ionic crystals of the cubical structure is considered in this paper.
The penetration of protons along the [100] directions and between the ( 100> planes of monocrystals is considered.
The Potentials of the atomic planes and the atomic rows of crystals are calculated and the critical angles for channelling for sodium and potassium haîo-genides combinations are defined.
The energy spectra of protons having passed through the monocrystals of various compositions are examined and the correlation between the shape of the energy spectrum and the ratio of radii in composition, between the share of particles channelling and the ordinal number of halide in rews is discovered.
В теории прохождения частиц через вещество основным приближением является замена действительной среды, через которую проходит частица, средой изотропной, гомогенной, составленной из отдельных атомов, беспорядочно расположенных по всему объему. Процессы рассеяния частиц на отдельных атомах независимы друг от друга и характеризуются дифференциальным поперечным сечением рассеяния. Торможение частиц в такой аморфной среде не зависит от направления движения частиц и распределение потерь энергии и угловое распределение частиц зависят только от числа атомов в единице объема или от толщины слоя, пройденного частицей.
Однако в большинстве своем, реальные твердые тела обладают вполне определенной кристаллической структурой, то есть являются представителями другого вида среды — среды гетерогенной, анизотропной и упорядоченной. Развитие радиационной физики в настоящее время остро ставит вопрос о роли кристаллической структуры во всех процессах взаимодействия частиц с веществом.
Первые исследования процессов взаимодействия частиц с монокрис-таялическими твердыми телами были связаны с, катодным распылением
веществ и привели к открытию Венером в 1954 году так называемых фокусирующих столкновений атомов [1].
Выполненные в течение последних лет теоретические и экспериментальные исследования прохождения быстрых заряженных частиц через моиокристаллические вещества показали, что кристаллическая структура оказывает сильное влияние и на характер движения частиц. При совпадении направлений пучка частиц и какой-либо главной кристаллографической оси монокристалла наблюдаются аномально низкие потери энергии как у протонов [2 4], так и у тяжелых ионов [5], уменьшаются выходы атомных [6] и ядерных [7ч-8] процессов, значительно снижается скорость образования радиационных дефектов [9], а также наблюдаются аномалии в рассеянии частиц на монокристаллических мишенях [10]. Все эти явления анизотропии связываются с движением частиц в относительно свободных пространствах решетки — «каналах», образованных совокупностями соседних, параллельных друг другу атомных рядов или атомных плоскостей, лежащих в определенных кристаллографических направлениях.
Захват частиц в каналы осуществляется последовательными скользящими столкновениями частиц с атомами рядов или плоскостей, образующих стенки каналов. Подобное явление получило название «кана-лирования» частиц в кристаллических решетках.
Следует отметить и другую группу явлений, связанных с влиянием кристаллической структуры на движение частиц, образующихся внутри самих кристаллов в ядерных реакциях [11] или рассеянных на большие углы [12]. Для таких частиц главные кристаллографические направления в отличие от частиц, каналирующих через кристаллы, оказываются за крытыми. Такое явление получило название «блокировки» или «эффек та теней.» [13].
Таблица 1
Соединение
№С1 ЫаВг КС1 КВг К1
Характе-ристики кристаллов N.
см 2,81 2,98 3,14 3,29 3,53
/V Ю-8 см 0,98 0,98 1,33 1,33 1,33
гл Г'А 0,54 0,5 0,73 0,68 0,61
Я 0,657 0,76 0,55 0,57 0,61
Тпл. °К 1077 1028 1049 1043 996
ед *к 306 243 221 185 155
Здесь
а—постоянная решетки [14], га—радиус аниона в соединении [14], гк—радиус катиона в соединении [14],
<7—коэффициент компактности структуры [15] Тпл.^-температура плавления [16], вд—характеристическая температура Дебая [17].
Отметим далее, что все исследования эффектов каналирования частиц к настоящему времени проведены в узком кругу простых веществ.
таких как: Si, Ge, Au, Си, обладающих различной кристаллической структурой и при различных энергиях частиц. Поэтому результаты работ трудно систематизировать и невозможно установить связь между изменением свойств кристаллической среды и ее влиянием на движение частиц. С этой точки зрения целесообразно исследовать эффекты каналиро-вания в монокристаллических химических соединениях, обладающих одинаковой кристаллической решеткой, но с различными ее параметрами. Исследование химических соединений тем более интересно, что можно проследить влияние кристаллической структуры в ряду веществ при условии постоянства одной из компонент соединения при изменении других.
Исходя из этих условий, нами в качестве объектов исследования выбраны монокристаллы галогенидов щелочных металлов. NaCl, NaBr, KCl, KBr, KI. Все эти соединения относятся к типу ионных кристаллов и обладают простой кубической решеткой.
Некоторые характеристики данных соединений представлены и табл. 1.
Наиболее открытыми направлениями в простой кубической решетке являются направления [100] и [110].
На рис. 1 показана структура канала [100] в монокристалле NaCl и схематично показано распределение электронной плотности в плоскости (100).
Рис. 1. Структура канала [100] в простой кубической решетке.
Ч*" — угол влета протона в канал, У о—амплитуда колебания траектории протонов, I? — радиус канала, с! — расстояние между соседними ионами в ряду, Я — длина волны колебания траектории протонов, и — порядковые номера ионов
Можно заметить, что относительно свободное от электронов пространство несимметрично относительно оси канала вследствие различных порядковых номеров ионов, расположенных в узлах решетки.
• Механизм каналирования
Анализ движения частиц по каналам решетки в классическом приближении, справедливом для протонов [18], показывает, что когда частица влетает в кристалл под достаточно малым углом к оси канала, она
55
в течение относительно длительного времени будет двигаться по осциллирующей траектории с амплитудой Уо в некотором среднем потенциальном поле ^(р), обусловленном всеми атомами в ряду или плоскости. При этом удобно рассматривать только поперечное движение частицы с энергией Е± =ЕХ¥2, где Е—энергия частицы и 4я — угол между направлением оси канала и направлением пучка частиц.
Потенциалы атомных рядов и атомных плоскостей легко получить, зная потенциалы взаимодействия частицы и (г) с отдельными атомами.
Известно, что при рассеянии протонов, с энергиями в интервале от от 1 кэв до 10 Мэв на малые углы особенно важным становится учет экранирования ядер атомными электронами.
Бор [19] для учета экранирования ввел потенциал взаимодействия между частицей и атомом вида:
= (I)
г \ ав]
Здесь и далее применяем следующие обозначения: 1хе — заряд протона, 12е — заряд рассеивающего атома,
г — расстояние от центра атома до налетающей частицы, ав — радиус экранирования.
Яв = а0 (2? + (2)
где а0 = 0,529-Ю-8 см —боровский радиус атома водорода.
В большинстве случаев более точным, особенно для тяжелых атомов и на малых расстояниях, является потенциал Томаса-Ферми, записанный в виде [20]
V (3)
Г \атф/
Здесь (—Функция экранирования Томаса-Ферми для отдель-
\ аТФ/ ного атома.
аТФ — радиус экранирования Томаса-Ферми.
аТФ = 0,8853 а0 (¿1 + Zf)-i. (4)
В дальнейшем так как то аТФ = 0,468Для значений г,
г
ЯТФ
сравнимых с радиусом экранирования, функцию х(—-) можно ап-
проксимировать выражением
г \ атФ
атФ / 2г и потенциал запишется в виде
(5)
и{г)=2*г>е1аТФ (6)
2 г
Для больших значений г, а именно для 0<г<6йтф более точным является потенциал, полученный Мольер [21]
0,1 еХр(- —) + 0,55 ехр(- 0,35 ехр ('
\ йтф / \ атФ / \ аТф / .
(7)
и(г) =
ZJZ2£2
Нами были получены выражения для потенциалов атомных рядор и атомных плоскостей в монокристаллах двухкомпонентных химических соединений, составленных из ионов с различными порядковыми номерами и Z¡.
Для атомных рядов потенциалы Бора (1), Мольер (7) и Нильсена (6) имеют вид:
+ 0,35/0,
\атФ
V(p) =
F(P) = 4-z:
Z,t?5 a
a 0Ж
Z.
аБ 0,1 K0 6p
Z:k0 л + Z:K0 -l.
aB
(H)
/6p_ \«тф
f 0,55A"0 (
Лтф
-5Л-„ ('
йтф
0.55Л;,
0,3/
Отф
0.35/0,
/азр
* (¿то
V(P)
7zZxe2 4 р-а
[Z^a-хф -г 22Лтф| .
9) (10)
Здесь
К0(г) — функция Макдональда [22], р — расстояние по радиусу от центра атомного ряда. Для атомных плоскостей для потенциалов Бора, Мольер и Ннль сена получены следующие выражения:
V{ р),
-■Z,e2
Z'3aB exp
X
0,55 1,2
0,1
— exp 6
л2
"Й(р)пл KZxe2 а2
ехр(- аТФ /
6Р х \ .0,55 ) 1,2
ЙТФ ,
Лб
' 0,1 ¿2атФ — exp 6
Z",aB exp (• 6P \
«Б)
атФ
0.35 0,3
exp r.Z-ув2
2a2
exp
0,3p
^ТФ
(И)
(12)
ЙТФ
, 0,35
'--exp
+ Zm-тф X 0,3p
/ 1 0,3
\Z\asx<b + 7ЛаТФ\.
<1тф
(13)
Рассчитанные зависимости потенциалов V($) для атомных рядов |100j и атомных плоскостей (100) в КВг представлены на рис. 2 и 3.
Сравним потенциалы рядов и плоскостей между собой. На рисунках видно, что потенциалы рядов и плоскостей малы до некоторых значений р ^ йтф , а затем наблюдается их резкое возрастание. Следует отметить, что потенциалы рядов уменьшаются быстрее с расстоянием, чем потенциалы поскостей; оба потенциала становятся равными на рас-d
стоянии-— 2 .отсюда эффекты ряда должны быть более чувствительны
ми к нарушениям решетки. Движение в межплоскостном пространстве, кроме того, является двухмерным. Отметим, что потенциалы рядов за-
Z,e
висят только от заряда на единицу ряда и не зависят от
энергии. Таким образом, каналы, расположенные вдоль главных кристаллографических направений, имеют не только большие размеры попе-
речного сечения, но и являются более эффективными с точки зрении влияния на характер движения частиц.
В рядах рассматриваемых соединений величины потенциалов возрастают с увеличением порядковых номеров.
Рис. 2. Средние потенциалы У(р) для рядов [100] в КВг.
Здссь н на следующем рисунке приняты следующие обозначения: 1 —потенциал Бора, 2 — потенциал Нельсона 3 — потенциал Мольер
Приближение непрерывных рядов справедливо на расстояниях р, больших, чем радиус экранирования атомов.
Для того чтобы это приближение было справедливым, необходимо потребовать, чтобы рассеяние частицы рядом определялось одновременно многими атомами в ряду. Это условие, очевидно, лучше выполняется при больших энергиях частиц. Можно показать, что каналирование становится возможным, если угол влета частицы в канал Ч1* меньше некоторого предельного угла *Р мах [23]
(Н)
£ е2
г де Ь - ——--диаметр столкновения в лабораторной системе коор-
Е
динат.
Величину предельного угла каналирования можно найти, если известно минимальное расстояние подхода частицы к ряду, при котором непрерывное описание рядов и плоскостей еще выполняется. Согласно Линдхарду [23] для устойчивого каналирования р т'т — йтф- С другой стороны, если траектория частицы лежит в плоскости атомного ряда, минимальное расстояние можно найти из решения уравнения
Рис. 3. Средние потенциалы V (р) пл для плоскостей (100)
в КВг
Отсюда
У(о) - ЕЧГК
Т
Е
(15)
(16)
Если теперь найти значение потенциала V (рт1П)> то можно получить значение максимального угла каналирования при данной энергии частиц.
Частицы, движущиеся по каналам решетки в среднем поле У(р) с поперечной энергией ЕЧ72, подходя к рядам или плоскостям на расстояние рт}п, встречают как бы потенциальный барьер. При высоких
скоростях частиц потенциальные барьеры по порядку величины равны:
для рядов
Л
2^2-20 эв
(18)
при й — 5а0;
для плоскостей ^^¿^.¿е^М-ипатф—^^р эв,
где йа — расстояние между плоскостями,
— среднее число атомов на единицу площади плоскости. Потенциальные барьеры для плоскостей ниже Ь ~ 41% раз барьеров рядов и отсюда предельные углы каналирования для плоскостей должны быть ниже в — 2Zí/6 раза предельных углов для рядов.
Если использовать выражение (16) для вычисления амплитуды, колебания траектории частицы в канале у ка.< функции угла У при различных энергиях частиц то можно получить кривые, представленные на рис. 4 и 5. Расчет выполнен для протонов с энергиями 1 Мэв, 5Мэв и ЮМэв для каналов [100] и (100] в КВг и для потенциала Бора.
24
22
\ го
/.8
16
10 М А, эб
1М А
210
4 5 6 7 Ъ 9 ю
У рад
210
Рис. 4. Амплитуда колебаний траектории протона в канале [100] КВг в зависимости от угла влета протона в кянял при различных значениях энергии
протонов
Из рисунков видно, что для различных энергий кривые расположены близко друг к другу и, следовательно, амплитуды у слабо зависят от угла ЧГ. Сближение кривых по мере приближения амплитуды к значению радиуса канала /?к и их насыщение обусловлено резким возрастанием потенциалов, начиная со значений
Ртт = Як — У0 ^ аТф. (19)
Принимая это значение для рт|п, можно найти предельные углы каналирования протонов в монокристаллах. В случае КВг для про -
тонов с энергией 6,72 Мэв углы для каналирования вадль ¡1001 направления и между плоскостями (100) равны
^[100] = 0,49°,
,1;тах(100) = 0,17°. (20)
Время одного колебания траектории частицы приблизительно равно [18]
4 у
S> /2
510 7 8 9 Ю
V рад
710
Рис. 5. Амплитуда колебаний траектории протона между плоскостями (100) КВг в зависимости от угла влета при различных значениях энергии протонов
и. следовательно, длина волны траектории ^ равна
л^. (22)
W
Соотношение между у и представленное h-í рис. 4 и 5, говорит о том, что л уменьшается с ростом ХЪ\ что и следовало ожидать. Из рис. 4 и 5 можно найти, что для протонов с энергией (5.72 Мэв для монокристалла КВг
ХП(Ю] ^ 2-10см,
(23)
/.(loo) — 3,9-10"5 см.
Соотношение между амплитудой У и углом "Ф" в ряду соединений позволяет сделать вывод о том, что при увеличении порядковых номеров компонент величина длины волны незначительно возрастает.
Теперь качественно рассмотрим поведение пучка частиц, проходящих через решетку кристаллов. При углах падения частиц на кристалл относительно главных направлений или плоскостей, меньших предель-
ного угла ^Р^Ч'ма* , частицы захватываются в каналы решетки. Подход их к атомным рядам на расстояние, меньшее радиуса экранирования, затруднен, и угловое распределение таких частиц будет отличным от углового распределения частиц в аморфной среде. Такой пучок получил название центрального пучка частиц. Если Ч1* ма,, то частицы, проходя через решетку, испытывают обычное рассеяние на отдельных атомах и в большинстве случаев их движение аналогично движению частиц в аморфной среде. Такой пучок называют нормальным или хаотическим. Интересно отметить, что эти два пучка не связаны между собой и переход частиц из одного пучка в другой в совершенном кристалле невозможен. Действительно, чтобы частица могла из центрального пучка перейти в хаотический, ей необходимо испытать очень сильное рассеяние на отдельном атоме, а для этого ей необходимо подойти к атому на рас стояние р<ат Ф, что для центральных частиц невозможно. По этой же причине затруднен и переход частиц из нормального пучка в центральный. Область углов 4я мах является промежуточной и рассмотрение поведения таких частиц чрезвычайно затруднено вследствие наложения многих эффектов.
До сих пор мы рассматривали модель жесткой решетки с атомами, неподвижно закрепленными в своих узлах. Оценим влияние тепловых колебаний атомов решением на величины потенциалов.
Рассмотрим гармонические колебания решетки с классической точки зрения. Такое рассмотрение дает хорошее приближение для температур больших температур Деоая 11олагаем, что смещения атомов подчиняются распределению Гаусса со среднеквадратичным смещением в2.
<24>
При таком распределении смещений атомов, потенциал вида
Т/(о) = Л ехр(— р/й),
где А и а - постоянные величины, преобразуется в потенциал вида [24]
Г(?)=У(р)ехр(^), (25)
где в'х — среднеквадратичное смещение в направлении перпендикулярном траектории частицы в канале. Согласно динамической теории кристаллических решеток для отклонения атома от своего положения равновесия можно записать [25]
(Т при Т > ед
о ^
10-2
^ при Т < вд. (26)
4
Это выражение справедливо только для простых кубических решеток, где выполняется условие = - $2.
3
Чтобы оценить влияние тепловых колебаний на величину- У(р), рассмотрим изменение потенциала в монокристалле КС1 с изменением температуры 7\ Для оценки используем потенциал Мольер (14) для атомной плоскости (100). Результаты расчета приведем в табл. 2.
Можно заметить, что тепловые колебания атомов приводят к незначительному росту потенциала, и этот рост вызывает соответственное и еще более слабое увеличение угла Ч'мах. С другой стороны.
тепловые колебания атомов приводят к уменьшению размеров канала Як на величину средней амплитуды колебания 5х. Теперь для устойчивого каналирования расстояние ближнего подхода к стенке канала рт$п будет определяться выражением [26]..
Ршт>С5±, (27)
где с — малая постоянная величина. Из простых соображений следует ожидать, что с будет являться функцией длины каналированной траек
Таблица с2
Потенциал Мольер для атомных плоскостей (100) в монокристалле КС1 с учетом тепловых колебаний ионов
р/Отф У( р)пл эв У (Р) 1л с учетом тепловых колебаний
г=зоэ°к Г=700°К Г=1049°К
1,0 29,7 30.6 31,6 31,2 36,8
2,0 20,2 20,9 21,41 23,2 25,0
3,0 14,4 14,9 15,3 16,51 17,7
4,0 10,52 10,7 11,2 12,1 13,1
5,0 7,73 7,96 8,21 8,9 9,6
гории протона которая приблизительно равна толщине монокристалла Зависимость с от t означает, что тах будет уменьшаться с увеличением и, следовательно, влияние температуры для толстых образцов более существенно, чем для тонких. Можно показать, что зависимость ^ тах от температуры дается выражением
тр (28)
х шах 1
Таким образом, повышение температуры образца приведет к снижению интенсивности каналированного пучка протонов при данной толщине образца
Потери энергии протонов в монокристаллах соединений
У тяжелых заряженных частиц при прохождении их через вещество энергия теряется в двух процессах: в неупругих столкновениях с атомными электронами (электронное торможение) и в упругих столкновениях с ядрами атомов (ядерное торможение). При высоких энергиях частиц электронное торможение преобладает и энергия частиц тратится на возбуждение или вырывание электронов в атомах. В хаотической системе ядерное торможение преобладает над электронным при г<С Ю, где е—безразмерная величина,, характеризующая энергию, введенная То масом-Ферми,
.. М> , (29)
ад*2 м{+мг
где
Мх — масса налетающей частицы (в нашем случае—протона),
М2 — масса ядра атома среды.
Так как нас интересует прохождение быстрых протонов через вещество, где ядерное торможение незначительно и составляеят лишь ~ 0,1% от
63
общих потерь энергии, то остановимся на ионизационных потерях энергии частиц.
При рассмотрении потерь энергии удобно выделить потери энергии на близкие столкновения с большой передачей импульса электронам и далекие — с малой передачей импульса электронам. Такое деление обусловлено тем, что в близких столкновениях необходимо учитывать кван-говомеханические явления, в то время как в далеких столкновениях они не играют роли в определении потерь энергии. Такое разделение удобно еще и тем, что каналированные частицы движутся в областях, где электронная плотность значительно ниже средней и, следовательно, роль близких столкновений в каналах невелика, в то время как далекие столкновения мало чувствительны к структуре вещества.
При больших энергиях и малых атомных номерах частиц для расчета средних потерь энергии справедлива формула Бете-Блоха, имеющая вид [23]
— = 5ЭЛ. • N • г,, = -1--N • В.
с1х эл ^ тъ2
(30)
Здесь
N—число атомов в единице объема,
5ЭЛ--поперечное сечение электронного торможения для аморфного вещества,
—
В— коэффициент торможения. В ^ 221п—— I (г) средний потенциал ионизации атома с атомным номером_г.
Для сложных соединений величину I находим по правилу Брэгга [27].
2 N¿11п Iг
1П / = ' 2 N#1 ' (31)
где N4 число атомов в единице объема с атомным номером И Н средний потенциал ионизации атомов элемента.
Нами вычислены потери энергии протонов с энергией 6,72 Мэв во всех исследуемых веществах с учетом поправок на эффект связи /(-электронов (табл. 3). Поправками на эффект плотности среды и на
Т а б л и ц а 3
Массовая тормозная способность соединений
Соединение I эв Число электронов в см3 иг«- юэз 1 [йЕ\ мэв р \ах) 2 см2
N301 189 6,248 42,9
№Вг 316 8,647 34,185
КС1 219 5,790 40,12
КВг 345 7,Г80 34,32
К1 479 8,179 29,86
эффект связи Ь — электронов в атомах при данной энергии протонов можно пренебречь [20]. Бете [28] при выводе формулы (31) для случая электронного газа установил, что при высоких скоростях частиц потери энергии на электронные столкновения асимптотически складываются из
равных долей от далеких и близких столкновений. Недавно Линдхард [29] подтвердил это правило для случая атомных систем. Исходя из этих замечаний, удобно разделить потери энергии на далекие и близкие [23].
~ (г) = 5ЭЛ [аЛГ • + (1 - а) Р(г)], (32)
ах
где р (г) —плотность электронов в некоторой точке пространства, через которое проходит частица. Второй член в выражении определяет потери
1 '
энергии в близких столкновениях с электронами. Коэффициента^ —
и при увеличении энергии приближается к значению 1/2. Это выражение показывает, что для быстрых частиц, проходящих по каналам решетки, потери энергии должны уменьшаться в два раза, поскольку вероятность близких столкновений для центрального пучка незначительна. Потери почти полностью будут определяться далекими столкновениями частиц с электронами среды. Так как далекие столкновения не чувствительны к точному расположению траектории протона в канале, то изменение температуры Т не должно заметно влиять на потери энергии протонов в монокристаллах.
Можно ввести эффективный атомный номер (£_|_), который определяет число электронов на атом, находящихся на расстояниях больших ртт от рЯДЭ, когдв
^(ртш) = Е±ш
Тогда поперечное сечение электронного торможения можно переписать в виде [23]
(Ел_) =--- Г< 1 - а) + а МЫ. V ¿"2
Можно показать, что отношение дается выражением
V" (ргп!п)
г,
(34)
и приблизительно равно —5-, если потенциал изменяется как р
Ех¥>
шах
• Потери энергии являются статистическим процессом и поэтому формула (31) определяет лишь средние потери энергии на единицу пути. Для тяжелых заряженных частиц при прохождении ими достаточных толщин функция распределения потерь энергии около среднего значения имеет формулу распределения Гаусса. Флуктуации потерь энергии на единице пробега частицы обусловлены только близкими столкновениями и, следовательно, мы вправе заменить среднюю плотность электронов Z2УV на р(г)
Тогда
¿(А Е)2
1
Лх ■ эл
или, используя эффективный атомный номер,
¿(Д£)а
г\ех?(г) (35)
Ах
= 4ъ11е*Е2(Е±). (36)
эл
Таким образом, полуширина аномального пика должна быть уже, чем нормального.
Заказ 7504
Отметим, что при 4х ^ Чт тах по физическим соображениям, должны наблюдаться потери энергии выше, чем у нормального пучка частиц вследствие блокировки на атомных рядах или плоскостях. Если угловое распределение каналированных частиц должно иметь максимум при ¥ = 0 и затем быстро уменьшаться с ростом угла, то угловое распределение блокированных частиц должно иметь максимум при углах ^щах• Такое поведение частиц наблюдали в монокристалле Б! [4].
Экспериментальные результаты
Эксперименты по исследованию потерь энеррш и каналированию протонов в монокристаллах галогенидов щелочных металлов нами проведены на экспериментальной установке, схема которой приведена на
* ' Рис. 6. Схема экспериментальной установки
рис. 6. Протоны 1 из циклотрона Томского политехнического института проходят в камеру рассеяния, где падают на тонкую платиновую мишень 2 толщиной около 3 мг/см2, установленную в центре камеры. Пучок протонов, рассеянный под углом 30° от мишени, коллимировался тонкими танталовыми диафрагмами и проходил в камеру гониометра, где помещался исследуемый монокристаллический образец 3. Диаметр пятна на образце 2 мм, энергия пучка протонов 6,72 Мэв, и угловая расходимость пучка меньше+ 0,1°.
Тонкие монокристаллические образны различных соединений с толщинами от 15мг/см2 до 50 мг/см2, предварительно ориентированные осью [100] по направлению пучка, укреплены на поворотном устройстве гониометра, имеющем регулировку в горизонтальной и вертикальной плоскостях с точностью установки углов+ 0,05 и+ 0,1 градусу соответственно. Энергию протонов после прохождения образцов определяли сцинтцлля-дионным спектрометром 4. Энергетическое разрешение спектрометра для протонов с энергией 6,72 Мэв составляет 2|8%.
Энергетическое распределение протонов, прошедших монокристаллические образцы ЫаС1 толщиной 17,3 мг/см3 и КВт толщиной 41,2 мг/см2, представлено на рис, 7 и 8. В правой части рисунков приведен энергетический спектр протонов, падающих на образцы.,
Образцы поворачивали в горизонтальной плоскости вокруг оси [001] гак, что пучок протонов оставался параллельным плоскостям (010). Кривые 2, 3, 4, 5 представляют энергетические спектры протонов при углах поворота 0; 0,6; 1,2 и 6 градусов соответственно. Все спектры нормированы к одному максимальному значению. Незначительный сдвиг спектров в сторону меньших значений энергии и их уширение обусловлены увеличением толщины образца при его повороте.
¿ООО
(Ш
тол
«с-Ч»
3
500
Энергия протонов Мз$
Рис. 7. Энергетическое распределение протонов с энергией 6,72 Мэв после прохождения монокристалла КаС1 под различными углами относительна направления [100] 0° — 2, 0,6°— 3, 1,2° — 4, 6° — 5
Из рисунков видно, что протоны, падающие на монокристаллы в направлении, достаточно далеком от направления [100], имеют энергетическое распределение, ожидаемое для поликристаллических и аморфных веществ. Когда направление [100] монокристаллов совпадает с направлением пучка, спектр искажается и* на высокоэнергетич-ной части спектра появляется «хвост»; указывающий на то, что некоторая доля частиц теряет меньшую энергию, чем частицы в нормальном спектре. При увеличении угла поворота «хвост» постепенно исчезает. Число частиц в центральном пучке можно определить простым вычитанием нормального спектра из нормированного.
В табл. 4 приведены средние потери энергии протонов в КаС1 и КВг, вычисленные по формуле Бете-Блоха и найденные экспериментально, а также значение коэффициента анизотропии, представляющего собой отношение числа частиц в аномальной части спектра к общему числу частиц в спектре.
Можно отметить, что большая часть протонов находится в нормаль* ном пучке и ее энергетические потери довольно хорошо согласуются с расчетными. Средние потери энергии пучка протонов, найденные из энергетических спектров, оказываются несколько ниже за счет той доли частиц, которые каналируют через кристаллы. Доля частиц в канали-
5*. 67—
рованном пучке при нулевом угле равна 0,36 для КВг и 0,28 для №С1. При увеличении угла поворота кристалла доля каналированных частиц уменьшается и средние потери протонов все более приближаются к расчетным.
60 * 6.72
Энергия протонов мэб
Рис. 8. Энергетическое распределение, протонов с энергией 6,72 Мэв после прохождения монокристалла КВг под различными углами относительно направления [100] 0°—2, 0,6°—3, 1,2°—4, 6°—5
На рис. 9 представлена зависимость коэффициента анизотропии от угла влета протонов в канал [100] для ЫаС1 и КВг. Для обоих кристал-
Таблица 4
Потери энергии протонов в монокристаллах соединений
Соединение Угол поворота в градусах Толщина образца мг\см* Л£теор. кэв ^норм. эксп. кэв ^■^ЭКСП. кэв N норм.
0 17,3 742 760+40 650+40 0,28
N301 0,6 -17,3 —742 760 ±40 710+40 0,14
1,2 -17,3 -742 760+40 740+40 0,04
6 — 17,31 —743 770+40 770+40 -0,0
0 41.2 1410 1420+60 1250+60 0,36
КВг 0,6 -41,2 - 1410 1420+60 1320+60 0,266
1,2 —41,2 -1410 1420+60 1380+60 0,11
6 -41,24 -1411 1440+60 1440+60 -0,6
лов ход кривых аналогичен. Угловые полуширины кривых для ЫаС1 и КВг равны ~ 0,68° и 0,8° соответственно. Согласно теоретическим оценкам (20) полуширины кривых должны составлять величины ~ 0,25°
для КВг и ~ 0,0 2° для №С1. Тот факт, что «хвосты» исчезают значительно медленнее, чем предсказывает теория, подтверждает вывод о том, что при высоких энергиях протонов ориентационные эффекты обусловлены главным образом атомными плоскостями, а не атомными рядами монокристаллов.
У^град
Рис. 9. Зависимость доли каналированных протонов в монокристаллах KBr (1) и NaCl (2) от угла влета протонов в канал [100]
Ниже обсудим некоторые факторы, которые могут повлиять на изменение формы спектров. Поскольку потери энергии протонов существенно определяются распределением электронной плотности в пространстве между плоскостями (010) или в направлении [100], мы вправе полагать, что изменения в распределении электронной плотности в ряду соединений от NaO до Nal и от КО до К1 приведет и к изменению формы энергетического спектра. Для наглядности обратимся к рис. 1. Мы уже отмечали, что пространство, свободное от электронов, в ионных кристаллах несимметрично относительно оси [100], что обусловлено различием в порядковых номерах ионов в узлах решетки. Несимметрия распределения электронной плотности возрастает в рядах галогенидов натрия и калия с ростом порядковых номеров галоидов. По мере движения частицы вдоль направления [100] картина распределения электронной плотности будет поворачиваться на угол 90° при переходе от одного атомного слоя к другому за счет перемены ионов в узлах решетки. Ста-
новится очевидным, что все пространство, заключенное между атомными рядами [100], удобно разбить на три части.
Частица, траекторию которой можно заключить в цилиндр с осью вдоль [100] направления и с радиусом основания рь движется в области, где теряет энергию преимущественно на далекие столкновения с атомными электронами. Потери энергии у таких частиц составляют приближенно половину потерь энергии нормального пучка частиц. При больших углах траектория частицы может оказаться в области, ограниченной цилиндром с радиусом рг, тде кроме далеких стсЛкновений в потерях энерщй начинают играть роль и близкие столкновения с электронами, принадлежащими отрицательно заряженным ионам с большими порядковыми номерами. Потери энергии у таких частиц несколько больше, чем потери энергии у первой группы частиц, но меньше, чем у частиц в нормальном пучке. И, наконец, при углах влета в фисталл мы имеем дело с нормальным пучком частиц. Соотношение между этими тремя областями в межионном пространстве и определяет форму энергетического спектра протонов. Большая несимметричность в распределении электронной плотности в монокристалле NaCl, чем в монокристалле КВг, и объясняет тот факт, что аномальная часть спектра в случае КВг выражена отчетливее в случае NaCl.
Заключение
1. Рассмотрено влияние кристаллической структуры соединений галогенидов щелочных металлов, обладающих простой кубической структурой на прохождение быстрых протонов.
2. Для соединений NaCl, NaBr, Nal, KCl, КВг, KI найдены потенциалы атомных рядов и атомных плоскостей, С увеличением порядкового номера галоида в ряду потенциалы возрастают.
3. Найдены предельные углы каналирования быстрых протонов в ионных кристаллах вдоль главных направлений [100] и [110] и между плоскостями (100)' и (110). В ряду соедйнений с увеличением порядкового номера галоида предельный угол возрастает. Захват протонов в каналы решетки наиболее эффективен в соединениях с большими порядковыми номерами.
4. Экспериментально исследованы потери энергии протонов с первичной энергией 6,72 Мэв в рядах соединений. При совпадении оси [100] монокристаллов с направлением пучка протонов часть протонов каналирует через монокристалл, теряя вдвое меньшую энергию по сравнению с неканалированным пучком.
5. Доля каналированных протонов возрастает в ряду соединений с увеличением порядкового номера галоида, что подтверждаеттеоретические -оценки.
• 6; Сдел ан выйод об' Определяющем вл иянии на каналированнё^бы-стрых протонов атоМй&х плоскостей кристаллов.
7. Каналирование в рядах галогенидов натрия и калия найб&лее выражено в монокристаллах NaCl и KCl и> уменьшается-в -рядйх^ соединений При увеличении порядкового номера галоида1 с изменений соотношения : между- радиусами ионов, в соединении.
литература
гзрбер, а! и.-Федоренко. уфн, т, 83, Кг 3, 385 (1964).
.2/G ТЕr.giihsqy/vH.VЕ. íWé g n е г/ W.M. Gibson.' Phys. Rev. Liters,: 13,
47,J 53(У (Í964).
' З.ЭДУ0 é а т ri á:l e у. ШЕЕ: Tr¡ms:'NúcI. Sei.; 11, 249 (1964).
4. W. M. Gibson, С. E r gin soy, H. E Wegner, B. R. Appletton. Phys. Rev. Letters, 15, № 8, 357 (1965).
5. S. Datz, T. S. Noggie, C. D. Mo a k. Phys. Rev. Letters, 15, 254 (1964).
6. W. Brandt, 1. M. Khan, D. L. Potter, R. D. Worly H. P. Smith. Phys. Rev. Letters, 14, №2, 42'(1965).
7. E. Bogh, I. A. Da vies, H. L. Nielson. Phys. Letters, 12, № 2, 129 (1964).
S. M. W. Thompson. Phys. Rev. Letters, 13, 256 (1964).
9. O. S. О en, T. S. Noggie. Phys. Rev. Letters, 16, № 10, 395 (1966).
10. W. M. Gibson. IEEE, Trans. Nuci. Sci., 13, № 3, 162 (1966).
11. B. Domeij, K. Bijoqvist. Phys. Letters., 14, № 2, 127 (1965).
12. E. Bogh, E. Uggerhoj. Phys. Letters, 17, № 2, 116 (1965).
13. А. Ф. Ту л и нов. УФН, т. 87, вып. 4, 585, 1965.
14. Ч. Киттель. Введение в физику твердого тела. Изд-во «Мир», Москва, 1963.
15. Г. С. Жданов. Физика твердого тела. Изд. МГУ, М., 1961. '
16. Д ж. К эй, Т. JI э б,и. Таблицы физических и химических постоянных. Физмат-гиз, М., 1962.
17. А. А. Воробьев. Физические свойства ионных кристаллических диэлектриков, Томск, 1960.
18. М. W. Thompson, R. S. Nielson. Phil. Mag., № 8 94, 1^77 (1963).
19. H. Б о p. .Прохождение атомных частиц через вещество. Техиздат, Москва, 1957.
20. С. В. Стародубцев, А. М. Романов. Прохождение заряженных частиц через вещество. Изд. АН Уз. ССР, Ташкент, 1962.
21. G. Moliere. Z .Naturforsch, 2а, 133 (1947).
22. Б. И. С е г а л, К. А. С е м е н д я е в. Пятизначные математические таблицы, ГИМФЛ, М., 1962.
23 I. Lin d hard. Phys. Letters, 12, № 2, 126 (1964).
Mat. Fys! Medd. Dan. Vid. Selsk., 34, № -14, (1965).
24. i; B, Sanders, I. M. Flint. Physica, 30, 129 (1964).
25. Дж. Лейбфрид. Макроскопическая теория механических и тепловых свойств кристаллов. Физматгиз, М., 1963.
26. В. W. Farmery, R. S. Nelson, R. Sizman, M. W. Thompson. Nuci Instr. and Meth., 38, 231, (1965). !
27. К. А. Воробьев, Б. А. Кононов. Прохождение электронов через вещество, Томск, 1&66.
28. W. Bethe. Z. Naturforsch, 4а, 88, (1949); 5а, 8, (1950).
'29. I. Lindhard, A. Winther.Mat.-Fys. Medd. Dan. Vid. Selsk., 34, Ns 4, (1964).