Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

SIMULATION OF THE INJECTION OF POSITRONS EMITTED BY RADIOISOTOPE SOURCE FOR INDUCTION CYCLIC ACCELERATOR

Laptev R.S. 1 Gogoleva T.S. 1 Rychkov M.M. 1 Shestak A.P. 1 Malikov E.L. 1 Lider A.M. 1
1 National Research Tomsk Polytechnic University
Induction cyclic accelerators of positrons (antibetatrons) are very promising technical basis as sources of monoenergetic particles with energies up to 10 MeV. Structurally, antibetatron is largely identical to the traditional betatron, except for particle injection system. The radioactive isotopes, emitting positrons, or positrons produced by the process of electron-positron pair production can be used as injectors for antibetatron. In this paper, the mathematical modeling of the injection of positrons emitted from radioisotope source 44Ti in the antibetatron accelerator chamber with azimuthally symmetric field and azimuthal variation of the field is done. As the accelerating system the electromagnetic system of MIB-6 betatron was selected. It is shown that the particles, making the betatron oscillations, after a number of turns again return to the source, and either, losing its energy at all and annihilate in the source material, or losing some of the energy and deposit on the chamber walls. The interactions of positrons with matter in these model are not considered. Application of the special radial ridge pole design in the antibetatron not allows to increase the efficiency of positrons capture during the acceleration. To ensure the sustainable capture of positrons the injection should be organized with the parameters close to the parameters of the electrons injected into the classical betatron.
Monte Carlo simulation
positron injector
betatron
antibetatron
1. Grafutin V.I., Prokop’ev E.P. Primenenie pozitronnoj annigiljacionnoj spektroskopii dlja izuchenija stroenija veshhestva // Uspehi fizicheskih nauk. 2002. t.172. no. 1. pp. 67–83.
2. Ivanilova T.S., Kashkovskij V.V. Modelirovanie zahvata jelektronov v uskorenie v betatrone s azimutal’noj variaciej polja //Nauchnaja sessija MIFI – 2009: Sbornik nauchnyh trudov – Moskva, 26-30 janvarja 2009. Moskva: MIFI, 2009. pp. 135–138
3. Ivanilova T.S., Kashkovskij V. V. Granica osvobozhdenija jelektronnogo puchka v betatrone s poljusami grebnevogo tipa // Izvestija TPU 2007 no. 2 Rezhim dostupa:: http://cyberleninka.ru/article/n/granitsa-osvobozhdeniya-elektronnogo-puchka-v-betatrone-s-polyusami-grebnevogo-tipa (data obrashhenija: 25.03.2013).
4. Ivanilova T.S. Modelirovanie zahvata i vyvoda jelektronnogo puchka betatrona s azimutal’noj variaciej polja: dis. kand. fiz.-mat. nauk. Tomsk, 2010. 116 p.
5. Kashkovskij V.V. Dinamika i parametry jelektronnyh puchkov betatronov. Tomsk: Izd-vo NTL. 2006. 196 p.
6. Krjoning H.M., Laptev R.S., Lider A.M. i dr. Issledovanie vozmozhnosti primenenija magnitnoj sistemy malogabaritnogo betatrona dlja indukcionnogo ciklicheskogo uskoritelja pozitronov // Kontrol’. Diagnostika. 2012. no. 13. pp. 192–196
7. Malogabaritnye impul’snye betatrony tipa MIB [Jelektronnyj resurs]. – Rezhim dostupa: http://www.epd.tpu.ru/science/mib1.htm (data obrashhenija: 25.03.13).
8. Malogabaritnye betatrony s azimutal’noj variaciej upravljajushhego magnitnogo polja [Jelektronnyj resurs] / V.L. Chahlov, A.A. Zvoncov. – Jelektronnye tekstovye dannye (1 fajl: 519 Kb) // Izvestija Tomskogo politehnicheskogo universiteta [Izvestija TPU]; Nerazrushajushhij kontrol’ i diagnostika: sbornik statej / Tomskij politehnicheskij universitet (TPU); pod red. V.L. Chahlova. Tomsk, 1998. [pp. 10–13]. – Zaglavie s titul’nogo lista. – Jelektronnaja versija pechatnoj publikacii.
9. Razrabotki laboratorii № 41: Kompaktnye uskoriteli jelektronov promyshlennogo i medicinskogo naznachenija [Jelektronnyj resurs]. – Rezhim dostupa: http://portal.tpu.ru/departments/institut/ink/Files/Lab41/RazLAb41.doc (data obrashhenija: 25.03.13).
10. Coleman P.G. Positron beams and their applications. // World Scientific Publishing. 2002. 322 p.
11. Krause-Rehberg R., Leipner H.S. Positron Annihilation in Semiconductors. // Berlin: Springer-Verlag. 1999. 397 p.

По своей природе позитроны обладают высокой чувствительностью к дефектам и могут использоваться для определения типа дефектов и их концентрации [1, 11]. Расширение сферы применения методов электрон-позитронной аннигиляции возможно при использовании пучков моноэнергетических позитронов с регулируемой энергией, что позволит проводить прецизионные исследования процесса аннигиляции на различных глубинах, то есть исследовать как поверхности и приповерхностные области, так и дефекты в объеме материалов [12]. Тем не менее, несмотря на потенциально широкую область использования пучков позитронов, их применение сильно ограничено из-за отсутствия доступной техники получения пучков позитронов с требуемыми параметрами. Использование индукционного циклического ускорителя позитронов – антибетатрона, один из наиболее перспективных способов получения моноэнергетических пучков позитронов.

Конструкционно антибетатрон во многом идентичен традиционному электронному бетатрону за исключением системы инжекции частиц. В качестве инжектора для антибетатрона возможно использовать радиоактивные изотопы, испускающие позитроны или позитроны, полученные в результате процесса рождения электрон-позитронных пар.

Целью данной работы было исследование возможности создания инжектора позитронов для индукционных циклических ускорителей на основе радиоактивных изотопов. Для оптимизации затрат на физическое моделирование предлагаемого инжектора позитронов целесообразно провести математическое моделирование процесса инжекции частиц от радиоизотопного источника для установления количественного соотношения между инжектированными и захваченными в ускорение позитронами.

В качестве ускоряющей системы выбран электромагнит серийного малогабаритного бетатрона МИБ-6 [7–9]. Проведено математическое моделирование процесса захвата позитронов для антибетатронов с азимутально симметричным полем и с азимутальной вариацией поля.

Параметры моделирования

Источник позитронов на основе изотопа 44Ti с активностью 0,91 МБк представляли в виде плоского диска диаметром 6 мм и толщиной 20 мкм, размещенного в вертикальной плоскости перпендикулярно орбитам ускоряемых частиц. Центр диска расположили на радиусе инжекции, где показатель спадания магнитного поля n = 1. Обычно на этом месте в традиционном бетатроне находится нить накала катода инжектора электронов.

Методом Монте-Карло осуществлялся выбор начальных параметров частиц от радиоизотопного источника: энергии из спектра источника, начальных координат r и z, начальных углов вылета частиц из источника по вертикали и горизонтали, времени влета, соответствующему энергии частиц на равновесной орбите. Начальные углы вылета частиц из источника в горизонтальной плоскости находились в интервале α = (0°,180°), а в вертикальной плоскости – в пределах γ = ±90° от горизонтальной плоскости. Напряженность магнитного поля соответствует движению позитрона с энергией 500 кэВ на равновесной орбите.

аpic_73.tifб
вpic_33.tifг

Рис. 1. Параметры частиц от радиоизотопного источника: а – начальные углы вылета частиц из источника в горизонтальной плоскости; б – начальные углы вылета частиц из источника в вертикальной плоскости; в – распределение частиц по координатам; г – распределение частиц по энергии

Электромагнитное поле антибетатрона считали постоянным. Движение одиночного позитрона в магнитном поле антибетатрона описывали с помощью систем дифференциальных уравнений (1) и (2):

Eqn117.wmf (1)

Eqn118.wmf (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):

Eqn119.wmf (3)

Eqn120.wmf (4)

Eqn121.wmf (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 инжектированных позитронов ни одна частица не попала в ускорение.

Рассмотрим параметры частиц, погибших на диске источника.

pic_34.tif

Рис. 2. Радиальное распределение вертикальной составляющей индукции магнитного поля, усредненное по азимуту

Таблица 2

Результаты моделирования процесса захвата позитронов от радиоизотопного источника магнитной системой МИБ-6 с азимутальной вариацией поля

Модель

МИБ-6 с азимутальной вариацией поля

Количество разыгранных историй

100000

Количество позитронов, сделавших 10000 оборотов

0

Количество позитронов, выпавших на стенку

99 802

Количество позитронов, погибших на источнике

198

Из рис. 3 видно, что произошло значительное перераспределение частиц по координатам. По сравнению с азимутально симметричным полем примерно на 70 % уменьшилось количество частиц погибших на источнике. Из рис. 4 и 5 видно, что в случае азимутально симметричного поля, значительная доля позитронов погибает на диске источника еще на первых оборотах (до 15 оборотов), причем в основном на источник попадают низкоэнергетические позитроны с энергией до 100 кэВ. Позитроны с энергиями, близкими к энергии на равновесной орбите (450–600 кэВ), совершают несколько оборотов и вновь попадают на диск источника.

аpic_35.tif б

Рис. 3. Распределение позитронов, погибших на диске источника по координатам: а – для азимутально симметричного поля; б – для поля с азимутальной вариацией

В случае же азимутальной вариации поля на диске источника погибают позитроны c энергиями 450–650 кэВ, что говорит о том, что позитроны были захвачены управляющим полем, однако, сделав несколько оборотов (до 50 оборотов), вновь попадают на диск источника.

а pic_36.tif б

Рис. 4. Энергетическое распределение позитронов, погибших на диске источника: а – для азимутально симметричного поля; б – для поля с азимутальной вариацией

аpic_37.tif б

Рис. 5. Распределение позитронов, погибших на диске источника, по числу полных оборотов: а – для азимутально симметричного поля; б – для поля с азимутальной вариацией

Заключение

Проведенное математическое моделирование процесса инжекции позитронов от радиоизотопного источника для антибетатронов с азимутально симметричным полем и с азимутальной вариацией поля позволяет сделать следующие выводы:

1. Радиоизотопный источник не удовлетворяет требованиям для осуществления устойчивого захвата частиц в ускорение, поскольку имеет непрерывный энергетический спектр и непрерывный угловой разброс.

2. Частицы, совершая бетатронные колебания, через некоторое число оборотов вновь возвращаются к источнику и, либо, потеряв всю энергию, погибают на диске источника, либо, потеряв часть энергии, оседают на стенках камеры. Процессы взаимодействия позитронов с веществом в данных моделях не учитывались.

3. Применение в антибетатроне специальной радиально-гребневой конструкции полюсов не позволит повысить эффективность захвата позитронов в ускорение.

4. Для обеспечения устойчивого захвата необходимо организовывать инжекцию позитронов с параметрами, близкими к параметрам инжекции электронов в классическом электронном бетатроне.

Таким образом, для повышения коэффициента захвата необходимо организовать импульсный ввод позитронов длительностью импульса 1–10 мкс и током не менее 109 частиц. Необходимо также иметь малый энергетический разброс (не более 1 %) и малый угол расходимости инжектируемого пучка (менее 10 градусов). Поэтому для создания ускорителя позитронов, способного ускорять 103–106 позитронов в секунду, целесообразна разработка специального импульсного инжектора позитронов с вышеуказанными параметрами.

Работа выполнена по государственному заданию «Наука», тема № 0.1146.2013.

Рецензенты:

Арефьев К.П., д.ф.-м.н., профессор, заведующий кафедрой высшей математики Физико-технического института, ФГБОУ ВПО «Национальный исследовательский Томский политехнический университет», г. Томск;

Потылицын А.П., д.ф.-м.н., профессор, заведующей кафедрой прикладной физики Физико-технического института, ФГБОУ ВПО «Национальный исследовательский Томский политехнический университет», г. Томск.

Работа поступила в редакцию 07.05.2013.