УДК 533.95
ЭФФЕКТ ГРУППИРОВАНИЯ ЭЛЕКТРОНОВ ПРИ УСКОРЕНИИ В КИЛЬВАТЕРНОЙ ПЛАЗМЕННОЙ ВОЛНЕ
Н. Е. Андреев, С. В. Кузнецов
Исследован процесс группирования электронов в кильватерной волне. Получены аналитические формулы, описывающие сжатие электронного сгустка, и выявлены условия реализации эффекта сильного группирования электронов. Результаты численного моделирования ускорения электронов подтверждают основные изученные закономерности.
Изучение процесса ускорения электронов в кильватерной волне пространственного заряда приобретает все больший интерес в связи с недавними успешными экспериментами по генерации сверхсильных ускоряющих полей (~ 100 ГэВ/м), создаваемых в плазме мощными короткими лазерными импульсами [1-5]. Для многих приложений (см., например, [6] и цитируемую там литературу) важной характеристикой ускоренных электронов является размер их сгустка (т.е. длительность электронного импульса). Поскольку в процессе ускорения в неоднородном поле кильватерной волны может проявляться эффект группирования электронов в пространстве, а также в распределении по энергиям [7, 8], лазерное ускорение в кильватерной волне открывает интересную возможность создания сверхкоротких сгустков релятивистских электронов [9, 10].
В настоящем сообщении исследуются основные закономерности группирования электронного сгустка конечной длины и предложена схема многократного сжатия электронного пучка и формирования последовательности сфазированных плотных коротких сгустков электронов, когда начальная длина пучка превышает длину плазменной волны.
Рассматривая процесс продольного сжатия сгустка электронов в направлении их ускорения, т.е. в направлении фазовой скорости волны Урь (вдоль оси 02), ограничимся сначала одномерным случаем, когда потенциал волны ф(£) зависит только от
сопутствующей переменной £ = кр(г - Ур^), где кр = шр/с, и>р = ^Атге2пе/гпс - ленг-мюровская частота электронов плазмы, и фазовая скорость волны близка к скорости света (Урн ~ с), как это имеет место для кильватерного поля, возбуждаемого лазерным импульсом в разреженной плазме (и>0 >> шр, ш0 - частота лазерного излучения).
Физический механизм сжатия пучка обусловлен продольной неоднородностью силы 1'\ = \е\дф/дг, когда на ускоряемые электроны в передней части сгустка действует сила, меньшая, чем на электроны в задней части. Поэтому наиболее сильно эффект группирования будет проявляться для электронов, проходящих область наибольшего отрицательного градиента силы (где д2ф/дг2 < 0), т.е. окрестность максимума потенциала ф((). Для моноэнергетического однородного сгустка электронов, инжектированного в окрестность максимума потенциала [£то — Д£0,£т], гДе Ф{£т) = Фтах, Д£о = крЬьо начальная безразмерная длина сгустка (Д£о < 1),в результате интегрирования уравнений движения электронов в потенциальном поле волны [7] находим изменение его длины Ьь по мере ускорения (движения относительно волны по координате £):
Ьь(0 _ Е1п,р'(0 + ИтМО
'60
ЕШп3 2 тсЕ(0
д2ф((т)
дС
крЬьо = аг (О + а2(ОкрЬь о- (1)
Здесь Е^Пу - начальная энергия инжектированных электронов и Е(£) - текущая энергия в лабораторной системе координат того электрона, который в начальный момент времени был инжектирован в точку максимума потенциала £т кильватерной волны, р'(£) его импульс в системе координат, связанной с волной, р'1П] - начальный импульс электрона в системе координат волны, 7 = 1/\/1 — Р2 — шр/и>о - гамма-фактор, определяемый фазовой скоростью волны /3 = Урь/с, /(£) - интеграл вдоль траектории движения
ВД = / + |в|7(^)-^т))]2/(шС2)2-1}
-3/2
где Е[П1 - начальная энергия электронов в системе координат волны. Предполагается, что энергия инжекции меньше резонансной Е1ти- < 7тс2 (т.е. электроны в системе координат волны движутся в области £ < £т).
Важным следствием формулы (1) является резкое уменьшение длины сгустка в окрестности точки поворота движения электронов относительно волны (их захвата), когда р'(£) и, следовательно, аг(^) стремятся к нулю. При этом минимальная длина сгустка определяется только его начальной длиной и формой потенциала волны в
точках максимума и захвата электронов:
Ьь = = —--дтг)т ■ (2)
В частности, для потенциала линейной гармонической волны вида ф(£) = — <^0со8(£) с максимумом в точке £т = тг максимальное сжатие сгустка согласно (2) достигается при захвате электронов в точке £(г = 7г/2, когда его минимальная длина
и = крЬ2ьо/2. (3)
Отметим, что в однородной плазме, как и в широком (по сравнению с к~г) плазменном канале, трехмерная структура линейной кильватерной волны такова [11, 12], что точка поворота = 7г/2 соответствует началу области фазы волны £ £ [тг/2, 37г/2], в которой радиальные силы кильватерной волны являются фокусирующими (что необходимо для устойчивого ускорения электронов в реальной трехмерной геометрии).
Соотношение (3) показывает, что для эффективного сжатия начальная длина элек тронного сгустка не должна превышать скиновой длины А;"1, причем степень сжатия (Ььо/Ьь) увеличивается с уменьшением начальной длины сгустка и не зависит непосредственно от амплитуды и фазовой скорости кильватерной волны (при выполнении, разумеется, условия захвата инжектированных электронов).
Для того, чтобы электроны с начальной энергией Е^ (инжектированные в окрест ности максимума потенциала фтах) были захвачены в заданной точке Т кильватерной волны, перепад потенциала должен удовлетворять условию:
|е|[<^тох - <К6г)] = - [(1 - 7~2)(Е1; - т^4)]1/2 - тс2/7. (4)
При этом для устойчивого сжатия электронного сгустка (с учетом его поперечного движения) точку захвата следует выбирать так, чтобы вся область движения электронов относительно волны в продольном направлении £ Е £т] находилась в фокусирующей фазе.
На рис. 1 в соответствии с формулой (4) изображена зависимость амплитуды потенциала линейной гармонической волны ф0 от энергии инжекции электронов при = 7г/2 для различных значений 7 = 15; 50; 100. Из графика видно, что с ростом фазовой скорости волны и уменьшением энергии инжектируемого пучка для осуществления его захвата необходима кильватерная волна с большей амплитудой. Однако для реально существующих в настоящее время инжекторов с Егп] и Юте2 амплитуда кильватерной волны, достаточная для захвата инжектируемого пучка, оказывается относительно
0.10
0.08
Г1
О
0.06
с
^ 0.04
0.02
0.0
0.0
0
10 20
30
40
0.0 0.05 0.10 0.15 0.20 0.25
Рис. 1. Зависимость амплитуды потенциала линейной кильватерной волны ф0, необходимой для захвата электронов на границе фокусирующей фазы при - 7г/2, от энергии инжекции электронов Е1П] для значений 7 = 15; 50; 100.
Рис. 2. Изменение длины сгустка по мере увеличения длины ускорения в кильватерной волм с параметрами 7 = 50, |е|ф0/тс2 = 0.01 при его инжекции в максимум потенциала с энергией Е1П] = 20 тс1; сплошная линия - аналитическая формула (1), точечные линии - трехмерно1 численное моделирование: квадраты - Ль/Яр = 0.067, кружки - 726/Лр = 0.033.
небольшой. Например, для пучка с энергией Е^ = 20 тс2 при 7 = 50 достаточно использовать кильватерную волну с амплитудой \е\ф0/тс2 > 0.01.
На рис. 2 показано, как изменяется в процессе захвата длина сгустка электронов (при начальной длине = ОЛА:"1) в кильватерной волне с 7 = 50 и \е\ф01тс2 = 0.01 по мере увеличения нормированной длины ускорения Lacc|Lvh, где Ьасс длина, прой денная электронным сгустком в лабораторной системе координат, Ьр^ — 1
максимальная длина ускорения (соответствующая прохождению электроном полного перепада потенциала от фт{п до фтах)- Сплошная линия показывает аналитический результат. вычисленный по формуле (1), в которой координата £ посредством решения уравнений движения [7] выражена через Ьасс/ЬР1г, а точечные кривые демонстрируют результаты трехмерного численного моделирования динамики электронов [8. 13] для двух значений относительной ширины электронного сгустка Яь/Я-р = 0.033 (кружки) и 0.067 (прямоугольники) при крЯр > 1, где Яр - характерный размер поперечного распре-
деления потенциала кильватерной волны (которое в данных расчетах предполагалось гауссовым: ф{£,г) = — ф0 соз(£) ехр[—(г/Др)2]). Для электронного сгустка, первоначально достаточно короткого (крЬьо <<1, что является условием применимости формулы (1)), рис. 2 показывает хорошее согласие одномерной аналитической теории (1) с результатами трехмерного численного моделирования при уменьшении относительной ширины пучка электронов во всем диапазоне изменения длины ускорения вплоть до захвата сгустка (Ьасс/Ь^ ~ 0.23).
Рис. 3. Коэффициент сжатия электронного сгустка в зависимости от его начальной длины; сплошная линия - аналитическая формула (3), отдельные точки - результаты численного моделирования при длине ускорения Ьасс/Ьрн = 0.234, отвечающей захвату сгустка при
= 1.67 яз 7г/2, в кильватерной волне с теми же параметрами, что и на рис. 2, и при той же энергии инжекции электронов.
Рис. 4. Распределение в фазовом пространстве электронов пучка (с начальной длиной, превышающей длину плазменной волны и Е= 20 тс2) в окрестности точки захвата = 1.67 для тех же параметров кильватерной волны, что и на рис. 2 (7 = 50, \е\ф0/тс2 = 0.01 при 1>асс/ЬрН = 0.234. Пунктиром изображена плотность электронов согласно формуле (5), заштрихованная область соответствует электронам, "вырезаемы" границей фокусирующей фазы волны.
Следует отметить, что область применимости простых асимптотических формул (2), (3) для минимальной длины электронного сгустка, достигаемой в процессе его сжа-
тия в момент захвата, оказывается шире ограничения крЬьо < 1, предполагавшегося при получении формул (1) и (2). Действительно, как показывает рис. 3, на котором представлена зависимость коэффициента сжатия ас = Ьь/Ььо от первоначальной длины электронного сгустка, выражение для ас — крЬьо/2 (которое следует из формулы (3) и изображено на рисунке сплошной линией) хорошо описывает результаты одномерного численного моделирования в широком диапазоне начальных длин сгустка вплоть до значений крЬьо ~ 7г/2 для Ьасс1Ьри = 0.234 (что отвечает захвату при = 1.67 » 7г/2 электрона, первоначально находившегося в точке максимума потенциала = 7г). Для меньших длин ускорения (т.е. до захвата электронов) выражение для коэффициента сжатия, полученное с помощью формулы (1), количественно хорошо описывает результаты моделирования при крЬьо <<1.
Рассмотренные выше характеристики процесса сжатия короткого (крЬьо < 1) сгустка электронов являются интегральными. Формула (2) позволяет определить также локальную плотность электронного сгустка в окрестности точки его захвата £(г (когда достигается максимальное сжатие):
где по начальная (предполагающаяся однородной) погонная плотность инжектируемого электронного пучка и 012(^т) определено соотношением (2) (для линейной волны «2(6т — 7Г/2) — 0.5). Распределение (5), имеющее интегрируемую особенность в точке захвата электронов, инжектированных в область максимума потенциала £т, описывает эффект пространственного группирования первоначально однородного электронного пучка. Этот эффект может быть использован для формирования короткого сгустка электронов с плотностью, превышающей начальную в инжектируемом пучке, и в случае когда длина инжектированного электронного пучка больше длины плазменной волны.
На рисунке 4 изображена часть фазового пространства (£,Рг) изначально однородного моноэнергетического электронного пучка (Е1ги- = 20 тс2) с длиной, превышающей длину плазменной волны, в окрестности точки захвата £гг = 1.67 для тех же параметров кильватерной волны, что и на рис. 3 (7 = 50, \е\ф0/тс2 = 0.01), при £асс/ЬРл = 0.234. Кружками показаны электроны, которые в момент инжекции занимали область фазы волны £ € [тт/2,7г], квадратами - область £ Е [7г,37г/2], т.е. все отмеченные электроны в начальный момент времени находились в фокусирующей фазе линейной кильватерной волны. Темные кружки и квадраты соответствуют электронам, которые находятся в фокусирующей фазе и в момент захвата. Пунктирная линия на рис. 4 демонстрирует
п
(О=0.5поМ6т)(6Т-ОГ1/2,
(5)
профиль погонной плотности электронов сгустка п(£)/п0 согласно формуле (5) с уче том двух ветвей распределения электронов в фазовом пространстве (что приводит к удвоению численного множителя в распределении (5)).
Электроны, попавшие в дефокусирующую фазу волны, выбрасываются из приосевой области сильного ускоряющего поля и будут потеряны для дальнейшего эффективного ускорения. Поэтому границей фокусирующей фазы осуществляется "вырезание" нам более сгруппированной части электронного пучка после его сжатия (при г < 7г/2). Таким образом, оба эти эффекта - сжатие и "вырезание" - позволяют выделить короткий сгусток электронов повышенной плотности из первоначально протяженного и менее плотного электронного пучка (см. заштрихованную область на рис. 4). Среднюю погонную плотность получаемого сгустка можно оценить с помощью распределения (5):
где крЬь(^Т) = — 7г/2 - длина сгустка, "вырезаемая" границей фокусирующей фазы (коэффициент 2 учитывает обе ветви распределения электронов в фазовом пространстве). Эта формула показывает, что при приближении точки захвата пучка к границе фокусирующей фазы средняя плотность сжатого короткого сгустка электронов возрастает. При этом, однако, полное число электронов в сгустке уменьшается как корень из конечной длины сгустка:
где Лр = 2тгк~1п0 - число электронов исходного пучка, приходящееся на длину плазменной волны. В случае, когда длина исходного инжектированного пучка электронов превышает длину кильватерной волны, рассмотренный механизм сжатия будет приводить к формированию нескольких коротких плотных сгустков электронов на границе каждой фокусирующей фазы волны (в окрестности £ = 7г/2 + 27гп, п — 0,1,2,...). При этом образованный цуг коротких сгустков будет сфазирован с кильватерной плазменной волной, что позволяет использовать все электроны цуга для дальнейшего эффективного ускорения в кильватерной волне с той же фазовой скоростью.
Для приведенного на рис. 4 примера средняя плотность сжатого сгустка электронов согласно (6) превышает начальную в ~ 8.9 раза. При этом его длина составляет ~ 0.016 от плазменной длины волны (крЬь(&г = 1-67) ~ 0.1), а полное число электронов
(6)
(7)
в сгустке согласно (7) ~ 14.2% от числа электронов инжектированного пучка, приходящихся на длину волны. Эти аналитически полученные значения хорошо согласуются с результатами численного моделирования, изображенными на рис. 4: Л/уЛ'р ~ 14.8%, < п > /гг0 — 9.3. При возбуждении кильватерной плазменной волны импульсом лазера с длиной волны 1 мкм указанные параметры соответствуют длине сжатого сгустка ~ 0.8мкм, длине плазменной волны Лр = 50мкм и длине группировочной секции
La.cc ~ 3 СМ.
Если рассмотренный выше способ сжатия и "вырезания" не позволяет сразу сформировать электронный сгусток нужной плотности (или нужного полного заряда) и достаточно малой длины, то можно применить его повторное сжатие в кильватерной волне с большим 7-фактором.
Исследованные процессы группирования первоначально протяженных однородных электронных пучков в кильватерной волне позволяют как эффективно сжать одиночной сгусток электронов, так и сформировать последовательность коротких плотных сгустков, что необходимо для эффективного моноэнергетического ускорения электронов в неоднородном поле кильватерной волны.
Работа выполнена при частичной поддержке РФФИ, грант N 98-02-16263.
ЛИТЕРАТУРА
[1] Nakajima К. et al. Phys. Rev. Lett., 74, 4428 (1995).
[2] M о d e n a A. et al. Nature, 337, 606 (1995); Clayton C. et al. 7th Workshop on Advanced Accelerator Concepts, Lake Tahoe, US, 12-18 October, 1996.
[3] U m s t a d t e r D. et al. Science, 273, 472 (1996).
[4] E s a г e у E. et al. IEEE Trans. Plasma Sei., 24, 252 (1996).
[5] A h д p e e в H. E., Горбунов Л. M. УФН, 168, (1998).
[6] О g a t a A. et al. IEEE Trans. Plasma Sei., 24, 252 (1996).
[7] Андреев H. E., Горбунов Л. M., Кузнецов С. В. Физика плазмы, 21, 1099 (1995).
[8] Andreev N. Е., Gorbunov L. М., Kuznetsov S. V. IEEE Trans. On Plasma Sei., 24, 448 (1996).
[9] С 1 а у t о n С. E., S е г a f i n i L. IEEE Trans. On Plasma Sei., 24, 440 (1996); К a t s о u 1 e a s Т. C., et al. IEEE Trans. On Plasma Sei., 24, 443 (1996).
[10] Ogata A., Nakajima К., Andreev N. J. Nuclear Materials, 248, 392 (1997).
[11] Горбунов JI. M., Кирсанов В. И. ЖЭТФ, 93, 509 (1987); Труды ФИАН, 219, 3 (1992).
[12] Andreev N. Е. et al. Phys. Plasmas, 4, 1145 (1997).
[13] Andreev N. E. et al. Proceedings of LASERS'97, New-Orleans, LU, December 1997.
Поступила в редакцию 30 октября 1998 г.