CC BY
69
Естественные науки
ИЗВЕСТИЯ
ТОМСКОГО ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ ПОЛИТЕХНИЧЕСКОГО ИНСТИТУТА им. С.М. КИРОВА
Том 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 cab culated and the critical anp-les for channelling for sodium
and potassium halo-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]. Для таких частиц главные кристаллографические направления в отличие от частиц, кана-
гт-1 м /- м
лирующих через кристаллы, оказываются закрытыми. Такое явление получило название блокировки или "эффекта теней" [13].
Отметим далее, что все исследования эффектов каналирования частиц к настоящему времени проведены в узком кругу простых веществ, таких как: Si, Ge, А^ С^ обладающих различной кристаллической структурой и при различных энергиях частиц. Поэтому результаты работ трудно систематизировать и невозможно установить связь между изменением свойств кристаллической среды и ее влиянием на движение частиц. С этой точки зрения целесообразно исследовать эффекты каналирования в монокристаллических химических соединениях, обладающих одинаковой кристаллической решеткой, но с различными ее параметрами. Исследование химических соединений тем более интересно, что можно проследить влияние кристаллической структуры в ряду веществ при условии постоянства одной из компонент соединения при изменении других.
Исходя из этих условий, нами в качестве объектов исследования выбраны монокристаллы галогенидов щелочных металлов. №0, NaBr, КС1, Ог, Ы. Все эти соединения относятся к типу ионных кристаллов и обладают простой кубической решеткой.
Некоторые характеристики данных соединений представлены в табл. 1. Наиболее открытыми направлениями в простой кубической решетке являются направления [100] и [110].
Таблица 1
Соединение
Характе- х ристики х, кристаллов ч. N301 ЫаВг КС! КВг К1
а- 10-е см гк-10™8 см ГЛ гь 2,81 0,98 0,54 2,98 0,98 0,5 3,14 1,33 0,73 3,29 1,33 0,68 3,53 1,33 0,61
Ч тпя. °к вд °К 0,657 (077 306 0,76 1028 243 0,55 1049 221 0,57 1013 ¡85 0,61 996 155
Здесь
а—постоянная решетки [14], га—радиус аииоиа в соединении [14], '"к—радиус катиона и соединении [14], д—коэффициент компактности структуры [15] Т„я.-^температура плавления [16],
—характеристическая температура Дебая [17].
На рис. 1 показана структура канала [100] в монокристалле и схематично показано распределение электронной плотности в плоскости (100).
Можно заметить, что относительно свободное от электронов пространство несимметрично относительно оси канала вследствие различных порядковых номеров ионов, расположенных в узлах решетки.
ектории протонов, Н — радиус канала, в — расстояние между соседними ионами в ряду, X — длина волны колебания траектории протонов, 12' и 12" — порядковые номера ионов
Механизм каналирования
Анализ движения частиц по каналам решетки в классическом приближении, справедливом для протонов [18], показывает, что когда частица влетает в кристалл под достаточно малым углом к оси канала, она в течение относительно длительного времени будет двигаться по осциллирующей траектории с амплитудой в некотором среднем потенциальном поле У(р), обусловленном всеми атомами в ряду или плоскости. При этом удобно рассматривать только поперечное движение частицы с энергией Е± =ЕТ2, где Е — энергия частицы и Т — угол между направлением оси канала и направлением пучка частиц.
Потенциалы атомных рядов и атомных плоскостей легко получить, зная потенциалы взаимодействия частицы и(г) с отдельными атомами.
Известно, что при рассеянии протонов, с энергиями в интервале от 1 кэв до 10 Мэв на малые углы особенно важным становится учет экранирования ядер атомными электронами.
Бор [19] для учета экранирования ввел потенциал взаимодействия между частицей и атомом вида:
Здесь и далее применяем следующие обозначения: Т1е — заряд протона, Т2е — заряд рассеивающего атома, г — расстояние от центра атома до налетающей частицы, аБ — радиус экранирования.
где а0 — боровский радиус атома водорода.
В большинстве случаев более точным, особенно для тяжелых атомов и на малых расстояниях, является потенциал Томаса-Ферми, записанный в виде [20]:
Здесь х I — I — функция экранирования Томаса-Ферми для отдельного атома; аТФ — радиус экранирова-
V атФ )
ния Томаса-Ферми.
—1
В дальнейшем так как Т1<Т2, то аТФ = 0,468 3. Для значений г, сравнимых с радиусом экранирования, функцию х I — I — можно аппроксимировать выражением
V атФ )
и потенциал запишется в виде
(5)
(6)
Для больших значений г, а именно для 0<г<6аТФ более точным является потенциал, полученный Мольер [21]
0,1 ёхр[ ——\+ 0,55 ехр/ —] + 0,35 ехр «тф i \ атф }
1,2
> «тф /!
(7)
Нами были получены выражения для потенциалов атомных рядов и атомных плоскостей в монокристаллах двухкомпонентных химических соединений, составленных из ионов с различными порядковыми номерами 7г и 7'г.
Для атомных рядов потенциалы Бора (1), Мольер (7) и Нильсена (6) имеют вид:
(9)
(10)
Здесь К(1) - функция Макдональда [22], р - расстояние по радиусу от центра атомного ряда.
Для атомных плоскостей для потенциалов Бора, Мольер и Нильсена получены следующие выражения:
У(р)ла. =
2;аБехр [--) -\-Z\as ехр (— р/'яв) | ■
ту, . < 0,1 / ■
а I 6 V «тф
X
6 V «№/ 1.2 \ ЙТФ
+ ¿1аТФ X
\ ■ 0^5
) 0,3 ехр( <1тфI
_ом
(11) (12)
(13)
— lгZ е1
v(fj)пл = —— [г>ТФ + 2;аТФ].
2а2
Рассчитанные зависимости потенциалов F(гmiI1) для атомных рядов [100] и атомных плоскостей (100) в КВг представлены на рис. 2 и 3.
Сравним потенциалы рядов и плоскостей между собой. На рисунках видно, что потенциалы рядов и плоскостей малы до некоторых значении р = яТФ, а затем наблюдается их резкое возрастание. Следует отметить, что потенциалы рядов уменьшаются быстрее с расстоянием, чем потенциалы плоскостей; оба потенциала становятся равными на расстоянии ~ й, отсюда эффекты ряда должны быть более чувствительными к нарушениям решетки. Движение в межплоскостном пространстве, кроме того, является двухмерным. Отметим, что потенциалы рядов зависят только от заряда на единицу ряда —- и не зависят от энер-
й
гии. Таким образом, каналы, расположенные вдоль главных кристаллографических направлений, имеют не только большие размеры поперечного сечения, но и являются более эффективными с точки зрения влияния на характер движения частиц.
В рядах рассматриваемых соединений величины потенциалов возрастают с увеличением порядковых номеров.
Приближение непрерывных рядов справедливо на расстояниях р, больших, чем радиус экранирования атомов.
Для того, чтобы это приближение было справедливым, необходимо потребовать, чтобы рассеяние частицы рядом определялось одновременно многими атомами в ряду. Это условие, очевидно, лучше выполняется при больших энергиях частиц. Можно показать, что каналирование становится возможным, если угол влета частицы в канал ¥ меньше некоторого предельного угла ¥мах [23]
(14)
где b =
E
-диаметр столкновения в лабораторной системе координат.
Величину предельного угла каналирования можно найти, если известно минимальное расстояние подхода частицы к ряду, при котором непрерывное описание рядов и плоскостей еще выполняется. Согласно Линдхарду [23] для устойчивого каналирования рш,п=аТф. С другой стороны, если траектория частицы лежит в плоскости атомного ряда, минимальное расстояние можно найти из решения уравнения
V(p) =
Отсюда
Е
(15)
(16)
Если теперь найти значение потенциала F(rmin), то можно получить значение максимального угла каналирования при данной энергии частиц.
Частицы, движущиеся по каналам решетки в среднем поле V(r) с поперечной энергией E¥2, подходя к рядам или плоскостям на расстояние min, встречают как бы потенциальный барьер. При высоких скоростях частиц потенциальные барьеры по порядку величины равны:
4Z,Z261
для рядов
: ZyZ2- 20 эв
(17)
при й~5а0;
для плоскостей ~ '■ ■ N • ¿„атч. л- 1 м, (18)
где йп - расстояние между плоскостями, Шп - среднее число атомов на единицу площади плоскости.
гн
2.2
£
•а
I
3
20
/8
(6
Уо'Як
ЮМ А,
{М А
2-Ю
4 5 6 7 Ъ 911)'г
У рад
г/о4
Рис. 4. Амплитуда колебаний траектории протона в канале [100] КВг в зави симости от угла влета протона в канал при различных значениях энергии
протонов
Потенциальные барьеры для плоскостей ниже Ь ~ 42\ раз барьеров рядов и отсюда предельные углы
1
каналирования для плоскостей должны быть ниже в ~ 222 раза предельных углов для рядов.
Если использовать выражение (16) для вычисления амплитуды колебания траектории частицы в канале у как функции угла ¥ при различных энергиях частиц Е, то можно получить кривые, представленные на рис. 4 и 5. Расчет выполнен для протонов с энергиями 1 Мэв, 5 Мэв и 10 Мэв для каналов [100] и [100] в КВг и для потенциала Бора.
Из рисунков видно, что для различных энергий кривые расположены близко друг к другу и, следовательно, амплитуды у слабо зависят от угла ¥. Сближение кривых по мере приближения амплитуды к значению радиуса канала Як и их насыщение обусловлено резким возрастанием потенциалов, начиная со значений
Принимая это значение для р™, можно найти предельные углы каналирования протонов в монокристаллах. В случае КВг для протонов с энергией 6,72 Мэв углы для каналирования вдоль [100] направления и между плоскостями (100) равны
Время одного колебания траектории частицы приблизительно равно [18]
Ау (21)
и, следовательно, длина волны траектории Я равна
(22)
Соотношение между у и ¥, представленное на рис. 4 и 5, говорит о том, что Я уменьшается с ростом ¥, что и следовало ожидать. Из рис. 4 и 5 можно найти, что для протонов с энергией 6,72 Мэв для монокристалла КВг
Соотношение между амплитудой у и углом ¥ в ряду соединений позволяет сделать вывод о том, что при увеличении порядковых номеров компонент величина длины волны незначительно возрастает.
Теперь качественно рассмотрим поведение пучка частиц, проходящих через решетку кристаллов. При углах падения частиц на кристалл относительно главных направлений или плоскостей, меньших предельного угла ¥<¥Мах, частицы захватываются в каналы решетки. Подход их к атомным рядам на расстояние, меньшее радиуса экранирования, затруднен, и угловое распределение таких частиц будет отличным от углового распределения частиц в аморфной среде. Такой пучок получал название центрального пучка частиц. Если ¥>¥мах, то частицы, проходя через решетку, испытывают обычное рассеяние на отдельных атомах и в большинстве случаев их движение аналогично движению частиц в аморфной среде. Такой пучок называют нормальным или хаотическим. Интересно отметить, что эти два пучка не связаны между собой
и переход частиц из одного пучка в другой в совершенном кристалле невозможен. Действительно, чтобы частица могла из центрального пучка перейти в хаотический, ей необходимо испытать очень сильное рассеяние на отдельном атоме, а для этого ей необходимо подойти к атому на рас стояние р<аТФ, что для центральных частиц невозможно. По этой же причине затруднен и переход частиц из нормального пучка в центральный. Область углов является промежуточной и рассмотрение поведения таких частиц
чрезвычайно затруднено вследствие наложения многих эффектов.
До сих пор мы рассматривали модель жесткой решетки с атомами, неподвижно закрепленными в своих узлах. Оценим влияние тепловых колебаний атомов решением на величины потенциалов.
Рассмотрим гармонические колебания решетки с классической точки зрения. Такое рассмотрение дает хорошее приближение для температур больших температур Дебая 0д. Полагаем, что смещения атомов подчиняются распределению Гаусса со среднеквадратичным смещением
йз = 2 ехр •
V 2-м
_
(24)
При таком распределении смещении атомов, потенциал вида где А и а - постоянные величины, преобразуется в потенциал вида [24]
(25)
где sl - среднеквадратичное смещение в направлении перпендикулярном траектории частицы в канале. Согласно динамической теории кристаллических решеток для отклонения атома от своего положения равновесия можно записать [25].
' Т при Т вд
Г1 ^ г.п.г 'I и (26)
■И
Это выражение справедливо только для простых кубических решеток, где выполняется условие ^ =1 я2-
Чтобы оценить влияние тепловых колебаний на величину У(р), рассмотрим изменение потенциала в монокристалле КС1 с изменением температуры Т. Для оценки используем потенциал Мольер (14) для атомной плоскости (100). Результаты расчета приведем в табл. 2.
Таблица 2
Потенциал Мольер для атомных плоскостей (100) в монокристалле КС! с учетом тепловых колебании ионов
р/Ртф V (р)ш. эв I7' <Р) 1л с учетам тсллоимх колебаний
г=о°к 7=30:)^ 7 =700°К 7=Ю49°К
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
Можно заметить, что тепловые колебания атомов приводят к незначительному росту потенциала, и этот рост вызывает соответственное и еще более слабое увеличение угла ^м^. С другой стороны, тепловые колебания атомов приводят к уменьшению размеров канала Лк на величину средней амплитуды колебания Б1. Теперь для устойчивого каналирования расстояние ближнего подхода к стенке канала р^ будет определяться выражением [26].
Рш1п ^ , (27)
где с - малая постоянная величина. Из простых соображений следует ожидать, что с будет являться функцией длины каналированной траектории протона ¡, которая приблизительно равна толщине монокристалла. Зависимость с от / означает, что ^шах будет уменьшаться с увеличением ¡, и, следовательно, влияние температуры для толстых образцов более существенно, чем для тонких. Можно показать, что зависимость ^шах от температуры дается выражением
(28)
Таким образом, повышение температуры образца приведет к снижению интенсивности каналирован-ного пучка протонов при данной толщине образца.
Потери энергии протонов в монокристаллах соединений
У тяжелых заряженных частиц при прохождении их через вещество энергия теряется в двух процессах: в неупругих столкновениях с атомными электронами (электронное торможение) и в упругих столкновениях с ядрами атомов (ядерное торможение). При высоких энергиях частиц электронное торможение преобладает и энергия частиц тратится на возбуждение или вырывание электронов в атомах. В хаотической системе ядерное торможение преобладает над электронным при 5<10, где 5 - безразмерная величина, характеризующая энергию, введенная Томасом-Ферми.
Е
ад,»2 м^мг
(29)
где М\ - масса налетающей частицы (в нашем случае - протона), М2 - масса ядра атома среды.
Так как нас интересует прохождение быстрых протонов через вещество, где ядерное торможение незначительно и составляет лишь ~ 0,1% от общих потерь энергии, то остановимся на ионизационных потерях энергии частиц.
При рассмотрении потерь энергии удобно выделить потери энергии на близкие столкновения с большой передачей импульса электронам и далекие - с малой передачей импульса электронам. Такое деление обусловлено тем, что в близких столкновениях необходимо учитывать квантовомеханические явления, в то время как в далеких столкновениях они не играют роли в определении потерь энергии. Такое разделение удобно еще и тем, что каналированные частицы движутся в областях, где электронная плотность значительно ниже средней и, следовательно, роль близких столкновений в каналах невелика, в то время как далекие столкновения мало чувствительны к структуре вещества.
При больших энергиях и малых атомных номерах 2Х частиц для расчета средних потерь энергии справедлива формула Бете-Блоха, имеющая вид [23]
йЕ
= ■ N ■ г., = 4 ■ /V • в.
(30)
йх даоа
Здесь N - число атомов в единице объема, - поперечное сечение электронного торможения для
2щу2 —
аморфного вещества, В - коэффициент торможения. В = Ъ21п —-— I (г)- средний потенциал ионизации
атома с атомным номером г.
Для сложных соединений величину I находим по правилу Брэгга [27].
2 1п ¡1
1П/ :
2 N¿1
(31)
где N1 число атомов в единице объема с атомным номером I, средний потенциал ионизации атомов элемента.
I айяица 3
Массовая тормозная способность соединений
Соединение 1 ,9в Число электронов Н СМ" 1 ¡йЕ \ мэн р [¡¿х) 2 см?
ЫаС! 189 6,248 42,9
№Вг 316 8,647 34,¡85
КС! 219 5,790 40,12
КВг 345 7, £80 34,32
К[ 479 8,179 29,86
Нами вычислены потери энергии протонов с энергией 6,72 Мэв во всех исследуемых веществах с учетом поправок на эффект связи ¿электронов (табл. 3). Поправками на эффект плотности среды и на эффект связи Х-электронов в атомах при данной энергии протонов можно пренебречь [20]. Бете [28] при выводе формулы (31) для случая электронного газа установил, что при высоких скоростях частиц потери энергии на электронные столкновения асимптотически складываются из равных долей от далеких и близких столкновений. Недавно Линдхард [29] подтвердил это правило для случая атомных систем. Исходя из этих замечаний, удобно разделить потери энергии на далекие и близкие [23].
¿Е , , „ — И -
ах
-(Х)р(г)],
(32)
где р(г) - плотность электронов в некоторой точке пространства, через которое проходит частица. Второй член в выражении определяет потери энергии в близких столкновениях с электронами. Коэффициент а<1/2 и при увеличении энергии приближается к значению 1/2. Это выражение показывает, что для быстрых частиц, проходящих по каналам решетки, потери энергии должны уменьшаться в два раза, поскольку вероятность близких столкновений для центрального пучка незначительна. Потери почти полностью будут определяться далекими столкновениями частиц с электронами среды. Так как далекие столкновения не чувствительны к точному расположению траектории протона в канале, то изменение температуры Т не должно заметно влиять на потери энергии протонов в монокристаллах.
Можно ввести эффективный атомный номер ^.Е), который определяет число электронов на атом, находящихся на расстояниях больших р™ от ряда, когда
Тогда поперечное сечение электронного торможения можно переписать в виде [23]
I (1-а)
г,
(33)
(Е±)
Можно показать, что отношение —-—- дается выражением
¿2 №) Л " Рщ1а V' (Рш1п)
2-,
21,7-, еа
(34)
и приблизительно равно
Р1ТГ!
если потенциал изменяется как р-
Потери энергии являются статистическим процессом и поэтому формула (31) определяет лишь средние потери энергии на единицу пути. Для тяжелых заряженных частиц при прохождении ими достаточных толщин функция распределения потерь энергии около среднего значения имеет формулу распределения Гаусса. Флуктуации потерь энергии на единице пробега частицы обусловлены только близкими столкновениями и, следовательно, мы вправе заменить среднюю плотность электронов на р(г).
Тогда
с{х.
или, используя эффективный атомный номер,
= 4 т^е'Аг)
(35)
(36)
Таким образом, полуширина аномального пика должна быть уже, чем нормального.
Отметим, что при по физическим соображениям должны наблюдаться потери энергии выше,
чем у нормального пучка частиц вследствие блокировки на атомных рядах или плоскостях. Если угловое распределение каналированных частиц должно иметь максимум при ¥=0 и затем быстро уменьшаться с ростом угла, то угловое распределение блокированных частиц должно иметь максимум при углах Такое поведение частиц наблюдали в монокристалле 81 [4].
Экспериментальные результаты
Эксперименты по исследованию потерь энергии и каналированию протонов в монокристаллах галоге-нидов щелочных металлов нами проведены на экспериментальной установке, схема которой приведена на рис. 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 мг/см2 и КВг толщиной 41,2 мг/см2, представлено на рис. 7 и 8. В правой части рисунков приведен энергетический спектр протонов, падающих на образцы.
Образцы поворачивали в горизонтальной плоскости вокруг оси [001] так, что пучок протонов оставался параллельным плоскостям (010). Кривые 2, 3, 4, 5 представляют энергетические спектры протонов при углах поворота 0; 0,6; 1,2 и 6 градусов соответственно. Все спектры нормированы к одному максимальному значению. Незначительный сдвиг спектров в сторону меньших значений энергии и их уширение обусловлены увеличением толщины образца при его повороте.
Таблица 4
Потери энергии протонов в монокристаллах соединений
'оединетзие Угол поворота н градусах Толщина образца мг/см1 1Ете„р. кэв эксп- кэв кэв ^нор*.
0 17,3 742 760±40 6б0±40 0,28
МаС1 0,6 —^17,3 -742 760 £40 710±40 0,14
1,2 —17,3 —742 760±4О 740 ±40 0,04
К —17,31 -743 770±40 77О±40 -0,0
0 41,2 1410 1420 ±60 ■1250±60 0,36
КВг 0,6 -41,2 ~ 1410 !420±60 Ш0±60 0,266
1,2 —41 ,2 -1410 1420±60 1380±60 0,11
6 -41,24 ~1411 1440±60 1440^60 ~0,в
Из рисунков видно, что протоны, падающие на монокристаллы в направлении, достаточно далеком от направления [100], имеют энергетическое распределение, ожидаемое для поликристаллических и аморфных веществ. Когда направление [100] монокристаллов совпадает с направлением пучка, спектр искажается, и на высокоэнергетичной части спектра появляется "хвост"; указывающий на то, что некоторая до-
ля частиц теряет меньшую энергию, чем частицы в нормальном спектре. При увеличении угла поворота "хвост" постепенно исчезает. Число частиц в центральном пучке можно определить простым вычитанием нормального спектра из нормированного.
В табл. 4 приведены средние потери энергии протонов в №С1 и КВг, вычисленные по формуле Бете-Блоха и найденные экспериментально, а также значение коэффициента анизотропии, представляющего собой отношение числа частиц в аномальной части спектра к общему числу частиц в спектре.
Можно отметить, что большая часть протонов находится в нормальном пучке и ее энергетические потери довольно хорошо согласуются с расчетными. Средние потери энергии пучка протонов, найденные из энергетических спектров, оказываются несколько ниже за счет той доли частиц, которые каналируют через кристаллы. Доля частиц в каналированном пучке при нулевом угле равна 0,36 для КВг и 0,28 для №С1. При увеличении угла поворота кристалла доля каналированных частиц уменьшается и средние потери протонов все более приближаются к расчетным.
На рис. 9 представлена зависимость коэффициента анизотропии от угла влета протонов в канал [100] для №С1 и КВг. Для обоих кристаллов ход кривых аналогичен.
Угловые полуширины кривых для №С1 и КВг равны ~0,68° и 0,8° соответственно. Согласно теоретическим оценкам (20) полуширины кривых должны составлять величины ~ 0,25° для КВг и ~ 0,02° для
NaCl. Тот факт, что "хвосты" исчезают значительно медленнее, чем предсказывает теория, подтверждает вывод о том, что при высоких энергиях протонов ориентационные эффекты обусловлены главным образом атомными плоскостями, а не атомными рядами монокристаллов.
Ниже обсудим некоторые факторы, которые могут повлиять на изменение формы спектров. Поскольку потери энергии протонов существенно определяются распределением электронной плотности в пространстве между плоскостями (010) или в направлении [100], мы вправе полагать, что изменения в распределении электронной плотности в ряду соединений от NaCl до Nal и от КС1 до KI приведет и к изменению формы энергетического спектра. Для наглядности обратимся к рис. 1. Мы уже отмечали, что пространство, свободное от электронов, в ионных кристаллах несимметрично относительно оси [100], что обусловлено различием в порядковых номерах ионов в узлах решетки. Несимметрия распределения электронной плотности возрастает в рядах галогенидов натрия и калия с ростом порядковых номеров галоидов. По мере движения частицы вдоль направления [100] картина распределения электронной плотности будет поворачиваться на угол 90° при переходе от одного атомного слоя к другому за счет перемены ионов в узлах решетки. Становится очевидным, что все пространство, заключенное между атомными рядами [100], удобно разбить на три части.
Частица, траекторию которой можно заключить в цилиндр с осью вдоль [100] направления и с радиусом основания р2 движется в области, где теряет энергию преимущественно на далекие столкновения с атомными электронами. Потери энергии у таких частиц составляют приближенно половину потерь энергии нормального пучка частиц. При больших углах траектория частицы может оказаться в области, ограниченной цилиндром с радиусом р2, где кроме далеких столкновений в потерях энергий начинают играть роль и близкие столкновения с электронами, принадлежащими отрицательно заряженным ионам с большими порядковыми номерами. Потери энергии у таких частиц несколько больше, чем потери энергии у первой группы частиц, но меньше, чем у частиц в нормальном пучке. И, наконец, при углах влета в кристалл мы имеем дело с нормальным пучком частиц. Соотношение между этими тремя областями в межионном пространстве и определяет форму энергетического спектра протонов. Большая несимметричность в распределении электронной плотности в монокристалле NaCl, чем в монокристалле КВг, и объясняет тот факт, что аномальная часть спектра в случае КВг выражена отчетливее в случае NaCl.
Заключение
1. Рассмотрено влияние кристаллической структуры соединений галогенидов щелочных металлов, обладающих простой кубической структурой на прохождение быстрых протонов.
2. Для соединений NaCl, NaBr, Nal, KC1, KBr, KI найдены потенциалы атомных рядов и атомных плоскостей. С увеличением порядкового номера галоида в ряду потенциалы возрастают.
3. Найдены предельные углы каналирования быстрых протонов в ионных кристаллах вдоль главных направлений [100] и [110] и между плоскостями (100) и (110). В ряду соединений с увеличением порядкового номера галоида предельный угол возрастает. Захват протонов в каналы решетки наиболее эффективен в соединениях с большими порядковыми номерами.
4. Экспериментально исследованы потери энергии протонов с первичной энергией 6,72 Мэв в рядах соединений. При совпадении оси [100] монокристаллов с направлением пучка протонов часть протонов каналирует через монокристалл, теряя вдвое меньшую энергию по сравнению с неканалированным пучком.
5. Доля каналированных протонов возрастает в ряду соединений с увеличением порядкового номера галоида, что подтверждает теоретические оценки.
6. Сделан вывод об определяющем влиянии на каналирование быстрых протонов атомных плоскостей кристаллов.
7. Каналирование в рядах галогенидов натрия и калия наиболее выражено в монокристаллах NaCl и KC1 и уменьшается в рядах соединений при увеличении порядкового номера галоида с изменением соотношения между радиусами ионов в соединении.
ЛИТЕРАТУРА
1. Р. Н. Гарбер, А. И. Федоренко. УФН, т. 83, № 3, 385 (1964).
2. С. Erginsoy, H. Е. Wegner, W. М. Gibson. Phys. Rev. Letters, 13, № 17, 530 (1964).
3. G. Dearnaley. IEEE, Trans. Nucl. Sci., 11, 249 (1964).
4. W. М. Gibson, С. Erginsoy, H. E. Wegner, В. R. Appletton. Phys. Rev. Letters, 15, № 8, 357 (1965).
5. S. Datz, T. S. Noggle, C. D. Moak. Phys. Rev. Letters, 15, 254 (1964).
6. W. Brandt, I. М. Khan, D. L. Potter, R. D. Wcrty, H. P. Smith. Phys. Rev. Letters, 14, №2, 42 (1965).
7. E. Bogh, I. A. Davies, H. L. Nielson. Phys. Letters, 12, № 2, 129 (1964).
8. М. W. Thompson. Phys. Rev. Letters, 13, 256 (1964).
9. O. S. Oen, T. S. Noggle. Phys. Rev. Letters, 16, № 10, 395 (1966).
10. W. М. Gibson. IEEE, Trans. Nucl. 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. Дж. Кэй, Т. Лэби. Таблицы физических и химических постоянных. Физматгиз, М., 1962.
17. А. А. Воробьев. Физические свойства ионных кристаллических диэлектриков, Томск, 1960.
18. М. W. Thompson, R. S. Nielson. Phil. Mag., № 8 94, 1677 (1963).
19. H. Бор. Прохождение атомных частиц через вещество. Техиздат, Москва, 1957.
20. С. В. Стародубцев, А. М. Романов. Прохождение заряженных частиц через вещество. Изд. АН Уз. ССР, Ташкент, 1962.
21. G. Moliere. Z. Naturforsch, 2а, 133 (1947).
22. Б. И. Сегал, К. А. Семендяев. Пятизначные математические таблицы, ГИМФЛ, М., 1962.
23. I. Lindhard. Phys. Letters, 12, № 2, 126 (1964); Mat. Fys. Medd. Dan. Vid. Selsk., 34, № 14, (1965).
24. I. B. Sanders, I. М. Flint. Physica, 30, 129 (1964).
25. Дж. Лейбфрид. Макроскопическая теория механических и тепловых свойств кристаллов. Физматгиз, М., 1963.
26. В. W. Farmery, R. S. Nelson, R. Sizman, М. W. Thompson. Nucl. Instr. and Meth., 38, 231, (1965).
27. А. А. Воробьев, Б. А. Кононов. Прохождение электронов через вещество, Томск, 1966.
28. W. Bethe. Z. Naturforsch, 4a, 88, (1949); 5a, 8, (1950).
29. L. Lindhard, A. Winther. Mat.-Fys. Medd. Dan. Vid. Selsk., 34, №4, (1964).
