Известия РАН. Серия физическая, 2019, T. 83, № 2, стр. 211-213

Компенсация изменения фокусировки ондулятора при перестройке длины волны излучения Новосибирского ЛСЭ

А. С. Матвеев 1 2*, О. А. Шевченко 1, Н. А. Винокуров 1

1 Федеральное государственное бюджетное учреждение науки Институт ядерной физики им. Г.И. Будкеpа Сибирского отделения Российской академии наук
Новосибирск, Россия

2 Федеральное государственное автономное образовательное учреждение высшего образования “Новосибирский национальный исследовательский государственный университет”
Новосибирск, Россия

* E-mail: matveev.a.s@ya.ru

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

Аннотация

Перестройка длины волны излучения лазера на свободных электронах (ЛСЭ) приводит к изменению жесткости фокусировки ондуляторов. В работе рассмотрены два способа компенсации этого эффекта: сохранение транспортной матрицы и согласование параметров Твисса. Проверка этих методов была проведена на первом лазере новосибирского ЛСЭ. В результате удалось сохранить прохождение пучка при перестройке длины волны излучения в диапазоне 120–130 мкм.

ВВЕДЕНИЕ

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

$\lambda = {{\lambda }_{w}}\frac{{1 + {{{{K}^{2}}} \mathord{\left/ {\vphantom {{{{K}^{2}}} 2}} \right. \kern-0em} 2}}}{{2{{\gamma }^{2}}}},$
где $\lambda $ – длина волны излучения, ${{\lambda }_{w}}$ – длина периода ондулятора, $K$ – параметр ондулятора, $\gamma $ – лоренц-фактор. Имеются два способа перестроить длину волны излучения: изменить энергию электронного пучка, либо изменить параметр ондулятора K. Второй способ более предпочтителен, так как требует согласования магнитно-оптической системы локально на участке, содержащем ондуляторы.

Установка Новосибирский ЛСЭ включает в себя три лазера, способных генерировать излучение в трех различных диапазонах длин волн [1]. Данная работа проводилась на первом ЛСЭ, генерирующем ТГц-излучение в диапазоне 90–220 мкм. Дорожка этого ЛСЭ содержит два ондулятора, разделенных пустым промежутком, и 8 квадрупольных линз, по 4 с каждой стороны от ондуляторов (см. рис. 1).

Рис. 1.

Схема дорожки ондуляторов первого ЛСЭ.

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

МЕТОД СОХРАНЕНИЯ ТРАНСПОРТНОЙ МАТРИЦЫ

В первом подходе для сохранения транспортной матрицы достаточно 6 квадрупольных линз, так как отсутствует связь поперечных степеней свободы, и, значит, для каждой из них транспортная матрица 2 × 2 унимодулярная [2]. Таким образом, имеется следующая математическая постановка задачи:

$M\left( {{{g}_{1}}, \ldots ,{{g}_{6}},I} \right) = {{M}_{0}},$
где ${{g}_{k}}$ – градиент магнитного поля квадрупольной линзы k, $I$ – ток ондулятора, ${{M}_{0}}$ – исходная транспортная матрица. Мы имеем систему из 6 независимых уравнений для элементов матрицы M. Однако их выбор в общем случае не произволен: если имеется нулевой элемент, то диагональный к нему не может быть найден из условия унимодулярности. Поэтому в качестве независимых выбираются элементы той диагонали (главной или побочной), произведение элементов которой минимально по абсолютной величине, и один элемент с другой диагонали. В этом случае вычислительные ошибки в определении независимых элементов не приведут к существенной ошибке в оставшемся элементе, вычисляемом из условия унимодулярности.

Взяв полную производную по току ондулятора, получим систему уравнений:

$\frac{{\partial {{m}_{i}}}}{{\partial {{g}_{1}}}}\frac{{\partial {{g}_{1}}}}{{\partial I}} + \ldots + \frac{{\partial {{m}_{i}}}}{{\partial {{g}_{6}}}}\frac{{\partial {{g}_{6}}}}{{\partial I}} = - \frac{{\partial {{m}_{i}}}}{{\partial I}},\,\,\,i = 1, \ldots 6.$

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

$\frac{{\partial {{g}_{i}}}}{{\partial I}} = f({{g}_{1}}, \ldots ,{{g}_{6}},I),\,\,\,\,i = 1, \ldots 6.$

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

$\Delta = \sqrt {\sum\limits_{i,j} {\left( {M - {{M}_{0}}} \right)_{{ij}}^{2}} } .$
Рис. 2.

Зависимость токов квадрупольных линз от тока ондуляторов (а), ошибка в сохранении транспортной матрицы (б).

Данный подход обеспечивает точное решение, однако на практике не для всех начальных условий оно существует.

МЕТОД СОГЛАСОВАНИЯ ПАРАМЕТРОВ ТВИССА

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

Оптимальные параметры Твисса в центре пустого промежутка между ондуляторами были определены следующим образом. В виду симметрии альфа- и гамма-функции равны нулю, а бета-функция находится из условия минимума ее максимального значения в ондуляторах. Полученная зависимость оптимальной бета-функции в центре между ондуляторами от их тока представлена на рис. 3.

Рис. 3.

Оптимальное значение бета-функции в центре между ондуляторами в зависимости от их тока.

По аналогии с предыдущим методом получаем систему дифференциальных уравнений:

$\left( {\frac{{\partial {{M}_{{St}}}}}{{\partial {{g}_{k}}}}\frac{{\partial {{g}_{k}}}}{{\partial I}} + \frac{{\partial {{M}_{{St}}}}}{{\partial I}}} \right)\left( {\begin{array}{*{20}{c}} {{{\beta }_{0}}} \\ {{{\alpha }_{0}}} \\ {{{\gamma }_{0}}} \end{array}} \right) = \frac{\partial }{{\partial I}}{{\left( {\begin{array}{*{20}{c}} {\beta {\text{*}}} \\ {\alpha {\text{*}}} \\ {\gamma {\text{*}}} \end{array}} \right)}_{I}},$
где $\left( {{{\alpha }_{0}},{{\beta }_{0}},{{\gamma }_{0}}} \right)$ и $\left( {\alpha {\text{*}},\beta *,\gamma {\text{*}}} \right)$ – параметры Твисса на входе (выходе) и в центре пустого промежутка между ондуляторами соответственно, а ${{M}_{{St}}}$ – матрица Штеффена [2]. Таким образом, решение данной системы для левой и правой частей дорожки ондуляторов даст правило перестройки токов 8 квадрупольных линз, при котором параметры Твисса на входе и выходе будут согласованы с оптимальными параметрами в ондуляторах.

РЕЗУЛЬТАТЫ

На рис. 4 представлены результаты эксперимента для энергии пучка 12.7 МэВ, в котором перестраивалась длина волны излучения ЛСЭ: с компенсацией изменения фокусирующей силы ондуляторов (сплошная линия), без изменения токов квадрупольных линз (штриховая линия). Ток ондуляторов варьировался в пределах 1300–1500 А, что соответствует изменению длины волны излучения в диапазоне 120–130 мкм. Наблюдалось изменение тока пучка в поглотителе, который определяет прохождение пучка в канале. Компенсацией фокусирующей силы ондуляторов удалось сохранить прохождение почти на начальном уровне, для сравнения при ее отсутствии было потеряно более четверти пучка (см. рис. 4а). В ходе эксперимента сохранение транспортной матрицы контролировалось диагностикой матрицы откликов [4]. Для ее построения использовалось 4 корректора (2 по каждой поперечной оси) и 2 датчика положения пучка (ДПП). Ошибка в сохранении матрицы откликов $\Delta $ в ходе перестройки длины волны излучения показана на рис. 4б.

Рис. 4.

Эксперимент: сохранение транспортной матрицы: а – измеренный сигнал, пропорциональный току пучка в поглотителе; б – параметр, характеризующий изменение измеренной матрицы отклика. Штриховая линия соответствует отсутствию компенсации фокусировки.

Метод согласования параметров Твисса с оптимальными в ондуляторах также был экспериментально проверен, однако сохранить прохождение пучка при перестройке тока ондуляторов не удалось (см. рис. 5). Причиной данного результата может являться некорректность предположения об оптимальных параметрах Твисса для рассмотренного режима оптической системы.

Рис. 5.

Эксперимент: согласование параметров Твисса.

ЗАКЛЮЧЕНИЕ

В итоге удалось компенсировать изменение фокусирующей силы ондуляторов и сохранить прохождение пучка при перестройке длины волны излучения в диапазоне 120–130 мкм путем сохранения транспортной матрицы. Метод согласования параметров Твисса пучка требует дальнейшего развития диагностики параметров пучка в канале дорожки ондуляторов.

Работа выполнена с использованием оборудования ЦКП “СЦСТИ” на базе НЛСЭ ИЯФ СО РАН, поддержанного Минобрнауки России (уникальный идентификатор проекта RFMEFI62117X0012).

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

  1. Shevchenko O.A. et al. // Proc. RUPAC 2012. St. Petersburg. Russia. 2012. P. 136.

  2. Штеффен К. Оптика пучков высокой энергии. М.: Мир, 1989. 23 с.

  3. Калиткин Н.Н. Численные методы. М.: Наука, 1978. 237 с.

  4. Safranek J. // Proc. of the Joint US–CERN–Japan–Russia School on Particle Accelerators. Switzerland, 1998. 298 p.

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