Приборы и техника эксперимента, 2019, № 5, стр. 100-105

МАГНИТНЫЙ СКАНЕР ДЛЯ ВЫСОКОЭНЕРГЕТИЧНЫХ ПУЧКОВ ТЯЖЕЛЫХ ИОНОВ

А. А. Голубев a, А. В. Канцырев a, В. А. Панюшкин a*, О. С. Сергеева a, В. С. Скачков a, В. В. Васильев a

a Институт теоретической и экспериментальной физики им. А.И. Алиханова Национального исследовательского центра “Курчатовский институт”
117218 Москва, ул. Б. Черемушкинская, 25, Россия

* E-mail: Vsevolod.Panyushkin@itep.ru

Поступила в редакцию 29.03.2019
После доработки 29.03.2019
Принята к публикации 13.04.2019

Полный текст (PDF)

Аннотация

Для решения широкого круга прикладных задач по материаловедению и радиобиологии выполнен расчет компактного магнитного сканера, предназначенного для однородного облучения неподвижных мишеней размером до 200 × 200 мм пучком тяжелых заряженных частиц (т.з.ч.) с отношением заряда к массе Z/A = 0.3–0.5. Сканер оптимизирован для работы с пучком тяжелых ионов с отношением Z/A = 0.3, энергией – 800 МэВ/а.е.м. и эмиттансом – 40 π · мм · мрад как в одиночном, так и в мультикадровом режиме облучения. Разработанный сканер состоит из двух компактных сканирующих диполей (один за другим) и может быть размещен вместе с другими ионно-оптическими элементами на ограниченном пространстве канала транспортировки пучков т.з.ч. к мишенной камере, где осуществляется развертка по двумерному полю облучения. Представлены также и обсуждены основные характеристики источника импульсного питания сканера при медленном выводе пучков ионов.

ВВЕДЕНИЕ

Пучки высокоэнергетичных тяжелых заряженных частиц (т.з.ч.) являются уникальным инструментом, позволяющим создавать высокие температуры, давления и радиационные повреждения в исследуемых материалах. Круг проводимых исследований определяется разнообразием физических задач, спектром облучаемых материалов и условиями проведения экспериментов [14]. В частности, пучки т.з.ч. используют для обработки различных материалов [2], для магниторезонансной диагностики и облучения злокачественных образований в организме человека [3, 4].

В данной работе описаны результаты разработки конструкции и режимов работы магнитного сканера – основного узла экспериментальной установки по проведению прикладных исследований по воздействию т.з.ч. с энергией до 0.8 ГэВ/а.е.м. на электронную компонентную базу космических и летательных аппаратов, спецтехнику и биологические объекты [1].

ОПИСАНИЕ УСТАНОВКИ

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

Таблица 1.

Параметры ионного пучка

Тип иона d+, C6+, Li3+ Ar, Kr, Xe, Au
Z/A 0.45–0.5 0.3–0.5
Параметры Твисса αx = 2.22, βх = 13.94; αy = –0.65, βy = 3.56
Энергетический разброс 0.2 %
Длительность выведенного пучка, с 2–30
Интенсивность пучка, ионов/(cм2 · с) 10–106
Поле облучения (макс.), мм2 Sscan=200 × 200
Неоднородность облучения ≤±15%
Энергия ионов, MэВ/а.е.м. 250 800
εxy = εx), π · мм · мрад 40 10

Наряду с пучками ускоренных протонов и легких ионов d+, C6+, Li3+ наиболее высокой зарядности, в подобных экспериментальных установках используют также пучки тяжелых ионов Ar, Kr, Xe, Au с массовым числом до А = 197 и отношением заряда к массе Z/A = 0.3–0.45. Время медленного вывода пучка интенсивностью до 106 ионов/(cм2 ⋅ с) на облучаемый объект составляет до 30 с. Неравномерность распределения флюенса на объекте с наибольшим размером LS × LS = 200 × 200 мм должна быть ≤±15%. В процессе облучения исследуемый объект остается неподвижным.

На рис. 1 показана магнитная жесткость Ri, соответствующая полностью или частично ободранным ионам, которая должна быть обеспечена магнитами сканера для двух значений энергии – T = 0.25 и 0.8 ГэВ/а.е.м. – ионов с различным отношением Z/A.

Рис. 1.

Магнитная жесткость для ряда ионов.

Расчет сканера выполнен для наибольшей жесткости, 16 Tл · м, соответствующей иону 197Au59+ с пониженной зарядностью Z/A = 0.3 и для пучков с наибольшим эмиттансом.

Ввиду обычного дефицита пространства на канале транспортировки длину LB сканера, включающую длины LM1 вертикально (VDM – Vertical Deflection Magnet) и LM2 горизонтально (HDM – Horizontal Deflection Magnet) отклоняющих магнитов, а также LD – участок дрейфа ионов до облучаемого объекта, минимизируют.

При фиксированной длине, например LB = 4 м, параметры VDM сканера: радиус R поворота иона, длина магнита LM, поперечное отклонение иона на выходе из него h = LBsinφ – (LS/2)cosφ – связаны следующим соотношением с магнитным полем B:

где φ – угол преломления траектории иона.

На рис. 2 показана зависимость этих параметров от угла φ. Предпочтительный режим возникает вблизи φ = 1.52°. Здесь поле составляет B ≈ 0.9 Тл, длина магнита LM ≈ 0.5 м, в то время как уширение апертуры вследствие смещения пучка в поперечной плоскости не превышает 2h ≈ 15 мм. Значительно более жесткий режим у HDM (табл. 2). Приведенные значения являются исходными для оптимизации сканера.

Рис. 2.

Зависимость магнитного поля и геометрических параметров VDM от угла отклонения.

Таблица 2.

Расчетные параметры магнитов сканера

Магнит B, Tл LM , м R, м h , мм φ, градус
VDM 0.9 0.5 19 7 1.53
HDM 0.9 0.6 19 10 1.86

РЕЖИМЫ ОБЛУЧЕНИЯ

Для повышения однородности облучения предпочтительным представляется режим, при котором образец облучается при движении пучка вдоль ряда горизонтальных строк, соединенных участками поворота пучка за пределами поля облучения. Такой режим наиболее эффективен в установках, где не предполагается использования каких-либо систем автоматизированного управления интенсивностью пучка т.з.ч., подобно тому, как это предложено в [3]. К тому же система управления возбуждением магнитов сканера значительно проще по сравнению с разверткой в полярных координатах [4].

Интенсивность пучка описывается двойным гауссовым распределением:

$f(x,y) = \frac{{{\text{exp}}\left[ {--~\frac{{{{{(x--{{x}_{0}})}}^{2}}}}{{2{\sigma }_{x}^{2}}}~--~\frac{{{{{(y--{{y}_{0}})}}^{2}}}}{{2{\sigma }_{y}^{2}}}} \right]}}{{2{\pi }{{{\sigma }}_{x}}{{{\sigma }}_{y}}}}{\text{\;}},$
где (x0, y0) – точка пересечения оси пучка с плоскостью облучения xy; σx, σy – горизонтальное и вертикальное r.m.s. отклонения от оси соответственно. Задача выбора способа и положения границ области перемещения пучка состоит в определении временной зависимости параметров x0, y0, при которых обеспечивается допуск на равномерность распределения флюенса по поверхности облучаемого объекта.

В случае бесконечной длины строки перемещения следа пучка интенсивностью IB со скоростью ${v}$ вдоль оси x в бесконечных пределах образуется одинокая строка облучения в виде гребня гауссова профиля. При пробеге параметром y0 ряда значений: 0, a, 2a, …, (n – 1)a, где a – дистанция между строками, вся площадь Sscan объекта будет покрыта n строками, и из (1) следует распределение флюенса:

${{{\Phi }}_{{ns}}} = \frac{{{{I}_{B}}}}{{v}}\mathop \sum \limits_{i = 0}^{n--1} \frac{1}{{{{{\sigma }}_{y}}\sqrt {2{\pi }} }}~\,{\text{exp}}\left[ {--\frac{{{{{(y--ia)}}^{2}}}}{{2{\sigma }_{y}^{2}}}} \right].$

Здесь параметр IB/v есть вклад каждой строки в линейную плотность распределения флюенса; результирующая линейная плотность равна nIB/v.

При небольшой величине расстояния a профиль функции Фns, согласно рис. 3a, имеет почти плоскую вершину. С возрастанием a плотность флюенса падает в промежутках между строками. При a = = 2.5σy неравномерность δfΣ = (Фns – Фns_max)/Фns_max достигает допуска величиной 15%, назначенного для этого источника неоднородности (рис. 3б).

Рис. 3.

Неравномерность облучения: а – при расстоянии между строками a = 1.5σy; б – при a = 2.5σy.

Если, например, выбрать σy = 10 мм, то максимальное расстояние между строками следует ограничить величиной a < 2.5σy = 25 мм, при этом число строк, покрывающих LS = 200 мм, должно быть n > 8.

В случае конечной длины строки интегрирование (1) в конечных пределах показывает, что уровню 85% на фронте строки соответствует координата, удаленная от координаты оси стартового положения на ∼σx внутрь строки. Поэтому длина LI, по которой перемещается ось пучка, должна быть не менее чем на ∼2σx больше длины LS облучаемого образца: LILS + 2σx, так что распределение флюенса в этом случае имеет вид

где

$~{{{\xi }}_{{{\text{min}}}}} = \frac{{x--{{l}_{S}}--{{{\sigma }}_{x}}}}{{\sqrt 2 {{{\sigma }}_{x}}}};\quad {{{\xi }}_{{{\text{max}}}}} = \frac{{x + {{l}_{S}} + {{{\sigma }}_{x}}}}{{\sqrt 2 {{{\sigma }}_{x}}}}.$

На рис. 4 приведен пример распределения флюенса при облучении объекта пучком с параметрами σx = σy = 5 мм и интенсивностью IB = 106 ионов/с. Число строк однократного сканирования n = 21. След пучка перемещается по поверхности Sscan со скоростью ${v}$ = 200 мм/с (см. табл. 1). В этом случае создается резерв времени около 0.5 с на каждую строку для перевода пучка на следующую строку, осуществляемого за пределами объекта облучения.

Рис. 4.

Распределение флюенса на мишени.

При рассмотрении непрерывного и мультикадрового облучения мишени предельное распределение флюенса возникает при фиксированном времени Tк одного полного кадра и безграничном возрастании числа строк, n → ∞. Тогда

где ${{{\eta }}_{{{\text{min}}}}} = \frac{{y - l - {{{\sigma }}_{y}}}}{{\sqrt 2 {{\sigma }_{y}}}};$ ${{\eta }_{{{\text{max}}}}} = \frac{{y + l + {{\sigma }_{y}}}}{{\sqrt 2 {{\sigma }_{y}}}}$.

Практически приблизиться к предельному распределению (4) возможно только за счет повышения скоростных и предельно допустимых электрических характеристик сканера и его источника питания. При этом часто [2] приходится значительно снижать величину поля и увеличивать длину сканера.

При многократном облучении площади Sscan влияние флуктуаций интенсивности пучка во времени уменьшается. Распределение полного флюенса в этом режиме по форме идентично (4): Фk = k Фпред, где k – число однотипных кадров.

ДИНАМИКА ПУЧКА В СКАНЕРЕ

На рис. 5 показан результат численного моделирования динамики пучка ионов 197Au59+ в режиме: LD = 4.0 м, эмиттанс 40 π · мм · мрад.

Рис. 5.

Транспортировка пучка к объекту облучения: момент сканирования на максимальные углы для HDM (a) и для VDM (б); фазовые портреты пучка на входе сканера (в) и его выходе (г).

Согласующий квадрупольный триплет (с.т.), предшествующий сканеру, оптимизируется по достижении удовлетворительного компромисса между размерами пучка в HDM и на объекте. Пара диполей сканера наиболее компактна. Их возбуждение соответствует полям до B1 = 1.31 Тл, B2 = 1.20 Тл и обеспечивает поворот пучка на максимальные углы в обеих плоскостях (рис. 5a, 5б). При этом кроссовер занимает положение вблизи плоскости объекта (рис. 5в, 5г). (Для удобства сравнения на диаграмме рис. 5г выходной фазовый портрет xx' отражен симметрично по x и x'.)

МАГНИТЫ СКАНЕРА И ТРИПЛЕТА

Разработка диполей сканера проводилась комплексно: параллельно с расчетами динамики ионного пучка в с.т. (рис. 5а, 5б) и процесса покрытия сканером площади Sscan. Эта процедура состояла в коррекции геометрических и магнитных параметров триплета и сканера для обеспечения надлежащего прохождения пучка через VDM и HDM: отклонении ионного пучка на необходимые углы при одновременном снижении габаритных размеров магнитов и минимизации размеров пучка на объекте. В результате оптимизации VDM и HDM получены наиболее компактные диполи. Конструкция обоих магнитов, изображенная на рис. 6а, однотипная, а их параметры приведены в табл. 3.

Рис. 6.

Оптимизированные диполь сканера (a) и квадруполь с.т. (б).

Таблица 3.

Параметры магнитов сканера

Параметры магнитов VDM HDM
Длина полюса, мм 280 320
Диаметр апертуры, мм 44 70
Bmax при максимальном возбуждении, Тл 1.5 1.5
Im, ток возбуждения (максимальный), кA 1 1
RD обмотки магнита, мОм 6 13
L магнита, мГ 4 13
τ магнита, с 0.8 1.1
Область 0.5%-ной однородности поля, мм 40 50

Диполи в канале расположены на минимальном расстоянии друг от друга: расстояние между их полюсами 210 мм, а между обмотками 10 мм. Согласующий квадрупольный триплет позволяет снизить диаметр пучка с эмиттансом εx = εy = = 40 π · мм · мрад в триплете до ≤70 мм; апертура квадруполей в с.т. должна быть ≥∅80 мм.

На рис. 6б изображен общий вид квадрупольного магнита с эффективной длиной 330 мм и апертурой (под пучок) 70 мм, основные параметры магнита приведены в табл. 4.

Таблица 4.

Параметры квадрупольного магнита

Максимальный градиент, Тл/м 22
Апертура и ее рабочая часть, мм Ø80/Ø70
Ток возбуждения (максимальный), А 500
Нелинейность поля в рабочей области <5 · 10–3
Число витков на полюс 31
Размеры шинки, мм 10 × 10 × Ø6
Плотность тока, А/мм2 7
Длина сердечника/габаритная, мм 280/380
Габариты (ширина/высота), мм 700 × 700

РЕЖИМ ПИТАНИЯ СКАНЕРА

Реактивная компонента L импеданса электромагнитов, использующихся в каналах транспортировки, часто имеет существенную величину, из-за чего постоянная времени τ этих магнитов обычно сравнима (или превышает) с длительностью Ts импульсов тока возбуждения, требующейся для формирования строки. Поэтому для линейного изменения тока в нагрузке с неизменной индуктивностью в диапазоне токов i(t) ∈ [–Im, Im] напряжение U(t) источника в течение промежутка времени Ts  должно быть линейной функцией времени:

$U(t) = {{R}_{D}}{{I}_{m}}\frac{{t + {\tau }--{{T}_{s}}{\text{/}}2}}{{{{T}_{s}}{\text{/}}2}}.$

На рис. 7 показаны диаграммы изменения тока возбуждения и напряжения на обмотке HDM такого сканера при Im = 1 кА, τ= Ts = 1 с, RD = 13 мОм, L = 13 мГ.

Рис. 7.

Изменение тока возбуждения и напряжения на обмотке HDM сканера.

По окончании перемещения следа пучка вдоль строки выполняется переход на следующую строку, когда напряжение U(t) на HDM поддерживается неизменным и равным напряжению на активной составляющей RD импеданса электромагнита. В течение этого времени выполняется перемещение положения следа пучка в вертикальном направлении с помощью VDM.

Для снижения базового расстояния LB в нашем сканере принят достаточно высокий уровень поля – до 1.5 Тл (номинальное значение составляет 1.3 Тл). Вследствие нелинейности намагничения магнитопровода диполей при линейном режиме питания сканера на краях строк возникнет замедление развертки пучка; скорость сканирования ${v}$ может снизиться вдвое, если применяются электротехнические стали с невысокой индукцией насыщения. В результате существенно повышается неоднородность облучения образца, обусловленная обратной зависимостью (3) плотности флюенса от скорости ${v}$. Для снижения влияния этого фактора предусматривается питание сканера от источника с регулируемым во времени выходным напряжением. Это позволяет увеличить скорость изменения тока возбуждения диполей в окрестностях точек поворота и выровнять ее со скоростью облучения центральной зоны исследуемого образца.

ВЫВОДЫ

Для решения широкого круга прикладных задач выполнен расчет компактного магнитного сканера и режимов его работы для однородного облучения пучком т.з.ч. исследуемых мишеней размером до 200 × 200 мм, с отношением заряда к массе Z/A = 0.3–0.5.

Список литературы

  1. Butenko A., Kozlov O., Syresin E., Trubnikov G., Bakhmutova A.V., Bogdanov A.V., Gavrilin R., Golubev A., Kantsyrev A.V., Liakin D.A., Markov N.V., Panyushkin V.A., Skachkov V., Skobljakov A., Visotski S.A. // Proc. of IPAC2017. Copenhagen, Denmark. accelconf.web.cern. ch/AccelConf/ipac2017/papers/tupva114.pdf

  2. Гикал Б.Н., Горбачев Е.В., Казаринов Н.Ю., Казача В.И., Казача Г.С., Лебедев Н.И., Макаров А.А., Мельников В.А., Миронов В.И., Рабцун С.В., Фатеев А.А. // Письма в ЭЧАЯ. 2005. Т. 2. № 3 (126). С. 97. http://www1.jinr. ru/Pepan_letters/panl_3_2005/17_gik.pdf

  3. Haberer Th., Becher W., Schardt D., Kraft G.// Nucl. Instrum. and Methods in Phys. Res. 1993. V. A330. P. 296.

  4. Kohno R., Hotta K., Dohmae T., Matsuzaki Y., Nishio T., Akimoto T., Tachikawa T., Asaba T., Inoue J., Ochi T., Yamada M., Miyanaga H. // International Journal of Particle Therapy. 2017. V. 3. № 4. P. 429. https://doi.org/10.14338/IJPT-16-00017.1

Дополнительные материалы отсутствуют.