По своей природе позитроны обладают высокой чувствительностью к дефектам и могут использоваться для определения типа дефектов и их концентрации [1, 11]. Расширение сферы применения методов электрон-позитронной аннигиляции возможно при использовании пучков моноэнергетических позитронов с регулируемой энергией, что позволит проводить прецизионные исследования процесса аннигиляции на различных глубинах, то есть исследовать как поверхности и приповерхностные области, так и дефекты в объеме материалов [12]. Тем не менее, несмотря на потенциально широкую область использования пучков позитронов, их применение сильно ограничено из-за отсутствия доступной техники получения пучков позитронов с требуемыми параметрами. Использование индукционного циклического ускорителя позитронов – антибетатрона, один из наиболее перспективных способов получения моноэнергетических пучков позитронов.
Конструкционно антибетатрон во многом идентичен традиционному электронному бетатрону за исключением системы инжекции частиц. В качестве инжектора для антибетатрона возможно использовать радиоактивные изотопы, испускающие позитроны или позитроны, полученные в результате процесса рождения электрон-позитронных пар.
Целью данной работы было исследование возможности создания инжектора позитронов для индукционных циклических ускорителей на основе радиоактивных изотопов. Для оптимизации затрат на физическое моделирование предлагаемого инжектора позитронов целесообразно провести математическое моделирование процесса инжекции частиц от радиоизотопного источника для установления количественного соотношения между инжектированными и захваченными в ускорение позитронами.
В качестве ускоряющей системы выбран электромагнит серийного малогабаритного бетатрона МИБ-6 [7–9]. Проведено математическое моделирование процесса захвата позитронов для антибетатронов с азимутально симметричным полем и с азимутальной вариацией поля.
Параметры моделирования
Источник позитронов на основе изотопа 44Ti с активностью 0,91 МБк представляли в виде плоского диска диаметром 6 мм и толщиной 20 мкм, размещенного в вертикальной плоскости перпендикулярно орбитам ускоряемых частиц. Центр диска расположили на радиусе инжекции, где показатель спадания магнитного поля n = 1. Обычно на этом месте в традиционном бетатроне находится нить накала катода инжектора электронов.
Методом Монте-Карло осуществлялся выбор начальных параметров частиц от радиоизотопного источника: энергии из спектра источника, начальных координат r и z, начальных углов вылета частиц из источника по вертикали и горизонтали, времени влета, соответствующему энергии частиц на равновесной орбите. Начальные углы вылета частиц из источника в горизонтальной плоскости находились в интервале α = (0°,180°), а в вертикальной плоскости – в пределах γ = ±90° от горизонтальной плоскости. Напряженность магнитного поля соответствует движению позитрона с энергией 500 кэВ на равновесной орбите.
![pic_73.tif](/fs/i/2013/6-4/pic_73_fmt.jpg)
![pic_33.tif](/fs/i/2013/6-4/pic_33_fmt.jpg)
Рис. 1. Параметры частиц от радиоизотопного источника: а – начальные углы вылета частиц из источника в горизонтальной плоскости; б – начальные углы вылета частиц из источника в вертикальной плоскости; в – распределение частиц по координатам; г – распределение частиц по энергии
Электромагнитное поле антибетатрона считали постоянным. Движение одиночного позитрона в магнитном поле антибетатрона описывали с помощью систем дифференциальных уравнений (1) и (2):
(1)
(2)
Здесь E – кинетическая энергия позитрона; E0 – кинетическая энергия, соответствующая движению частицы по равновесной орбите; c – скорость света в вакууме; α – угол между направлением проекции скорости электрона на плоскость равновесной орбиты Vr и направлением радиуса – вектора r; γ – угол между направлением полной скорости Vr и направлением ее проекции на плоскость равновесной орбиты Vr [4–5].
Моделирование процесса инжекции и захвата позитронов в антибетатроне с азимутально-симметричным полем
Конфигурация магнитного поля в рабочем зазоре антибетатрона с азимутально симметричным полем подробно описана в работе [6]. Результаты математического моделирования процесса захвата позитронов от радиоизотопного источника магнитной системой МИБ-6 с азимутально симметричным полем представлены в табл. 1.
Таблица 1
Результаты моделирования процесса захвата позитронов от радиоизотопного источника магнитной системой МИБ-6 с азимутально симметричным полем
Модель |
МИБ-6 с азимутально-симметричным полем |
Количество разыгранных историй |
100000 |
Количество позитронов, сделавших 10000 оборотов |
0 |
Количество позитронов, выпавших на стенку |
99 654 |
Количество позитронов, погибших на источнике |
346 |
Анализ результатов показал, что инжекция позитронов от радиоизотопного источника в ускорительную камеру антибетатрона носит случайный характер. Частицы влетают в поле антибетатрона со случайной величиной энергии, случайными значениями координат и углов вылета из источника и случайным значением времени влета, часто с энергией, не соответствующей энергии частиц на равновесной орбите. Дальнейшее движение позитронов в магнитном поле антибетатрона носит сложный характер: двигаясь по разнообразным орбитам, частицы одновременно совершают вертикальные, радиальные и продольные колебания относительно положения равновесной орбиты. Амплитуда колебаний при этом определяется начальными условиями вылета частиц из источника. Большая часть позитронов теряется на стенках камеры или гибнет на источнике еще на первых оборотах после инжекции в камеру.
Таким образом, для осуществления захвата позитронов от радиоизотопного источника в антибетатроне необходимо использовать дополнительные методы или устройства, улучшающие технику ввода позитронов в ускорительную камеру. Так как большая часть позитронов выпадает на стенку ускорительной камеры, один из возможных вариантов повышения эффективности захвата позитронов – это увеличение линейных размеров самой камеры. Другой способ ‒ это увеличение фокусирующих сил, действующих на пучок. Использование азимутальной вариации позволяет существенно увеличить фокусирующие силы поля бетатрона за счет применения радиально-гребневой конструкции.
Моделирование процесса инжекции и захвата позитронов в антибетатроне с азимутальной вариацией поля
Расчеты проводились для антибетатрона на базе магнитной системы бетатрона типа МИБ-6 с азимутальной вариацией поля [2, 4, 8].
Для моделирования процесса инжекции и захвата позитронов в ускорение в антибетатроне с азимутальной вариацией поля была определена оптимальная величина вариации поля по азимуту, создаваемой гребневой конструкцией ускорителя, которая составила d(r) = 0,10.
Для описания конфигурации магнитного поля в рабочем зазоре антибетатрона с азимутальной вариацией использовали теоретическую модель поля, согласно которой выражения для компонент магнитной индукции имеют вид (3)–(5):
(3)
(4)
(5)
где r, φ, z – радиальная, азимутальная и вертикальная координаты соответственно; r0 и B0 – радиус равновесной орбиты антибетатрона и величина индукции на этом радиусе соответственно; δ и ξ – определяют относительный вклад в суммарное поле ускорителя магнитных полей, обеспечивающих периодичную по азимуту структуру поля, определяемых r и z – компонентами векторного потенциала A; a1, a2, a3, b1, b2, b3 – амплитуды; k1, k2, k3 – параметры разделения; J0(k1r), N0(k1r), J1(k1r), N1(k1r) – цилиндрические функции Бесселя и Неймана нулевого и первого порядка, а Jν(k3r), Nν(k3r), Jν+1(k2r), Nν+1(k2r) – функции Бесселя и Неймана порядка ν и (ν + 1), n0 – усредненный по азимуту показатель спадания поля на равновесной орбите [2–5].
Параметры разделения и амплитуды в выражениях (1–3) определяются из краевых условий, которым должно удовлетворять магнитное поле бетатрона на равновесной орбите. Выражения (1–3) позволяют достаточно точно описать распределение магнитного поля как в рабочем зазоре ускорителя, так и за его пределами.
На рис. 4 приведено экспериментально снятое распределение вертикальной составляющей индукции магнитного поля по радиусу Bz(r, z = 0), усредненное по азимуту в медианной плоскости ускорителя.
При моделировании источник располагали на азимуте φ = 0° и радиусе инжекции rинж, а электромагнитное поле антибетатрона считали постоянным.
Результаты математического моделирования процесса захвата позитронов от радиоизотопного источника магнитной системой МИБ-6 с азимутальной вариацией поля представлены в табл. 2.
Использование радиально-гребневой конструкции поля не позволило увеличить коэффициент захвата позитронов, из 105 инжектированных позитронов ни одна частица не попала в ускорение.
Рассмотрим параметры частиц, погибших на диске источника.
Рис. 2. Радиальное распределение вертикальной составляющей индукции магнитного поля, усредненное по азимуту
Таблица 2
Результаты моделирования процесса захвата позитронов от радиоизотопного источника магнитной системой МИБ-6 с азимутальной вариацией поля
Модель |
МИБ-6 с азимутальной вариацией поля |
Количество разыгранных историй |
100000 |
Количество позитронов, сделавших 10000 оборотов |
0 |
Количество позитронов, выпавших на стенку |
99 802 |
Количество позитронов, погибших на источнике |
198 |
Из рис. 3 видно, что произошло значительное перераспределение частиц по координатам. По сравнению с азимутально симметричным полем примерно на 70 % уменьшилось количество частиц погибших на источнике. Из рис. 4 и 5 видно, что в случае азимутально симметричного поля, значительная доля позитронов погибает на диске источника еще на первых оборотах (до 15 оборотов), причем в основном на источник попадают низкоэнергетические позитроны с энергией до 100 кэВ. Позитроны с энергиями, близкими к энергии на равновесной орбите (450–600 кэВ), совершают несколько оборотов и вновь попадают на диск источника.
а б
Рис. 3. Распределение позитронов, погибших на диске источника по координатам: а – для азимутально симметричного поля; б – для поля с азимутальной вариацией
В случае же азимутальной вариации поля на диске источника погибают позитроны c энергиями 450–650 кэВ, что говорит о том, что позитроны были захвачены управляющим полем, однако, сделав несколько оборотов (до 50 оборотов), вновь попадают на диск источника.
а б
Рис. 4. Энергетическое распределение позитронов, погибших на диске источника: а – для азимутально симметричного поля; б – для поля с азимутальной вариацией
а б
Рис. 5. Распределение позитронов, погибших на диске источника, по числу полных оборотов: а – для азимутально симметричного поля; б – для поля с азимутальной вариацией
Заключение
Проведенное математическое моделирование процесса инжекции позитронов от радиоизотопного источника для антибетатронов с азимутально симметричным полем и с азимутальной вариацией поля позволяет сделать следующие выводы:
1. Радиоизотопный источник не удовлетворяет требованиям для осуществления устойчивого захвата частиц в ускорение, поскольку имеет непрерывный энергетический спектр и непрерывный угловой разброс.
2. Частицы, совершая бетатронные колебания, через некоторое число оборотов вновь возвращаются к источнику и, либо, потеряв всю энергию, погибают на диске источника, либо, потеряв часть энергии, оседают на стенках камеры. Процессы взаимодействия позитронов с веществом в данных моделях не учитывались.
3. Применение в антибетатроне специальной радиально-гребневой конструкции полюсов не позволит повысить эффективность захвата позитронов в ускорение.
4. Для обеспечения устойчивого захвата необходимо организовывать инжекцию позитронов с параметрами, близкими к параметрам инжекции электронов в классическом электронном бетатроне.
Таким образом, для повышения коэффициента захвата необходимо организовать импульсный ввод позитронов длительностью импульса 1–10 мкс и током не менее 109 частиц. Необходимо также иметь малый энергетический разброс (не более 1 %) и малый угол расходимости инжектируемого пучка (менее 10 градусов). Поэтому для создания ускорителя позитронов, способного ускорять 103–106 позитронов в секунду, целесообразна разработка специального импульсного инжектора позитронов с вышеуказанными параметрами.
Работа выполнена по государственному заданию «Наука», тема № 0.1146.2013.
Рецензенты:
Арефьев К.П., д.ф.-м.н., профессор, заведующий кафедрой высшей математики Физико-технического института, ФГБОУ ВПО «Национальный исследовательский Томский политехнический университет», г. Томск;
Потылицын А.П., д.ф.-м.н., профессор, заведующей кафедрой прикладной физики Физико-технического института, ФГБОУ ВПО «Национальный исследовательский Томский политехнический университет», г. Томск.
Работа поступила в редакцию 07.05.2013.