^Эбэорыша^аг
1ЬНЫХ]
проблем нашейной динамики
Изв.вузов «ПНД», т. 2, № 6,1994
УДК 621.385.633
КОНКУРЕНЦИЯ И КООПЕРАЦИЯ МОД В ЛАЗЕРАХ НА СВОБОДНЫХ ЭЛЕКТРОНАХ
Н.С.Гинзбург, М.И.Петелин
Рассмотрены особенности конкуренции и кооперации мод в ЛСЭ генераторах, работающих в режимах как однородного, так и неоднородцрго уширения линии активного вещества. С использованием аналитических и численных методов показано, что в режиме однородного уширения линии по мере превышения тока пучка над стартовым значением проходится следующая последовательность бифуркаций: режим одномодовой одночасготной генерации, режим периодической автомодуляции (самосинхронизации мод), режим стохастической автомодуляции. Эффективность энергообмена в последнем режиме может превышать КПД режима одномодовой генерации за счет реализации механизма стохастического торможения, когда электроны, провзаимодействовав с одной модой и отдав ей часть своей энергии, оказываются в резонансе с другой модой и т.д. КПД при этом растет пропорционально ширине спектра излучения. Показано, что при использовании некогерентной накачки или стохастических ондуляторов аналогичный механизм торможения и повышения КПД может быть реализован при монохроматической сигнальной волне.
Влсдснис
Данный обзор посвящен динамике генераторов когерентного электромагнитного излучения, основанных на взаимодействии электромагнитных волн с потоками электронов, обладающих релятивистской поступательной скоростью. Благодаря доплеровскому смещению частоты последнее обстоятельство делает устройства такого рода пригодными для получения весьма высоких частот - вплоть до диапазона световых волн. Это и послужило основанием для того, чтобы назвать эти генераторы лазерами на свободных электронах (ЛСЭ) [1].
В ЛСЭ могут использоваться различные механизмы индуцированного излучения электронов: ондуляторное, циклотронное, черенковское излучение, вынужденное рассеяние волн. Нередко понятие ЛСЭ сводят к приборам, основанным только на ондуляторном излучении [1,2]. Однако в данном обзоре, за исключением раздела 5, последнее ограничение не подразумевается. Дело в том, что в довольно широкой области параметров теория всех ЛСЭ может быть построена универсальным способом. Для стационарных режимов это было сделано в [3-6]. Данный обзор посвящен главным образом нестационарным процессам, таким как нелинейная конкуренция мод, параметрическая неустойчивость рабочей моды, периодическая и стохастическая автомодуляция, взаимная синхронизация мод, некогерентная накачка и т.д.
Следует отметить, что указанные выше эффекты присущи весьма широкому классу распределенных автоколебательных систем, исследуемых и используемых в гидрофизике, физике плазмы, вакуумной электронике и в традиционных (квантовых) лазерах [7-18]. Некоторая специфика ЛСЭ ( по сравнению с квантовыми генераторами) состоит, очевидно, в том, что здесь активная среда образована релятивистскими частицами и поэтому обладает сильной пространственной дисперсией, масштаб которой соизмерим с длиной генератора.
Для описания многочастотных процессов в ЛСЭ может быть использован как модовый, так и пространственно-временной подход. %
Очевидно, что модовый подход полезен для описания сильного взаимодействия небольшого числа мод или произвольного числа слабосвязанных мод [19-21]. Напомним, что связь мод называется слабой [10, 17], если каждая мода черпает энергию на свое возбуждение из своей части (в пространстве скоростей или координат) электронного пучка. В этом случае моды сосуществуют почти независимо в одном генераторе. Связь называется сильной, когда электронный поток не может быть разделен на части, принадлежащие различным модам. Аналогично квантовым генераторам сильная связь мод реализуется в режиме однородного уширения линии активного вещества, а слабая связь мод возникает в режиме неоднородного уширения.
Сильная связь мод, имеющих сильно неэквидистантный спектр собственных частот, может быть сведена к совокупности бинарных взаимодействий и описана в рамках модового подхода. Однако, в случае сильной связи большого числа эквидистантных мод, обладающих почти одинаковой добротностью (случай типичный для квазиоптических резонаторов), модовый подход не столь привлекателен. Дело в том, что если возбуждены две произвольные моды, то возникает нелинейная комбинационная компонента электронного тока, которая возбуждает третью моду и т.д. В результате может возникать каскадное перераспределение энергии по спектру, поэтому трудно предсказать, какие моды будут участвовать во взаимодействии. В таких условиях более целесообразным становится пространственно-временной подход.
В традиционном виде пространственно-временной подход состоит в компьютерном моделировании взаимодействия электронов с электромагнитной волной на каждом проходе волны по резонатору [22-27]. Однако, чем ближе коэффициенты отражения от зеркал к единице и чем меньше огибающая изменяется на одном проходе, тем большее число проходов занимает переходный процесс
1
N =—-.
В этом случае удобно перейти [4,28,29] от дискретной переменной ■ омера прохода) к непрерывной (медленному времени) подобно тому, как в электродинамике описывается экспоненциальное затухание поля в холодном резонаторе.
В разделах 1-3 настоящего обзора в рамках пространственно-временного подхода анализируется сильная связь мод в условиях однородного уширения линии активного вещества. Исследован как режим стационарной, так и периодической импульсной инжекции электронного потока. В случае стационарной инжекции в рамках клистронной аппроксимации пространства взаимодействия аналитически исследована устойчивость режима одномодовой генерации и описан механизм возникновения автомодуляции.
В разделе 4 рассмотрена слабая связь мод в условиях неоднородного уширения линии активного вещества. Здесь для описания эффектов насыщения и определения стационарного распределения энергии излучения по спектру мод со слабо коррелированными фазами модовый подход дополнен квазилинейным приближением
121].
Характерной особенностью ЛСЭ как с однородным, так и неоднородным уширением линии является то, что при больших превышениях над порогом генерации, когда число возбуждающихся мод очень велико, реализуется механизм стохастического торможения электронов. В этом режиме электрон, провзаимодействовав с одной модой и отдав ей часть энергии, оказывается в
синхронизме с другой модой и т.д. В результате электронный КПД растет > пропорционально ширине спектра излучения и может существенно превосходить КПД режима одномодовой одночастотной генерации. Таким образом КПД оказывается тем выше, чем ниже степень когерентности излучения. Однако в ЛСЭ, основанных на вынужденном ондуляторном излучении и- вынужденном рассеянии волн, механизм стохастического торможения может быть реализован и при когерентной монохроматической сигнальной волне, если использовать стохастическую накачку (стохастические ондуляторы) [30]. Подобный метод повышения КПД анализируется в разделе 5 настоящего обзора.
1. Вывод нестационарных уравнений ЛСЭ-генераторов
Пусть генерируемое излучение заключено между двумя зеркалами с коэффициентами отражения 2 (рис.1). В таком резонаторе поле может быть представлено в виде двух встречных квазиоптических пучков
А=11е{уо[Ат+ (2,г)ехр;((Оо/- Л+ + Ат~ (2,/)ехр/(соо/+/го2)]},
где й)0 - несущая частота, /г0=Л(юо). Условия синхронизма в рассматриваемом классе генераторов могут быть Рис. 1. Схема ЛСЭ с двухзеркальным квазиоп-Записаны в виде тическим резонатором. Я\ 2 - коэффициенты
отражения от зеркал, I - расстояние между Ю-/Шц =£},., (1) зеркалами.
где - частота осцилляции электронов. Для ЛСЭ, основанных на вынужденном рассеянии волн и вынужденном ондуляторном излучении, частота осцилляций ^х^Ф+й^н, гДе Юь /г, - частота и волновое число поля накачки1, и условие (1) эквивалентно условию синхронизма между частицами и виртуальной комбинационной волной
=г-'с=юс//гс, (2)
где ас=(й - ее»}, /гс=/г+/?; - частота и волновое число комбинационной волны.
Взаимодействие волнового пучка на т-ом проходе по резонатору с электронным потоком описывается с учетом дисперсионного расплывания следующими уравнениями
Эат+
дат
дг
(-
v,
8Г
дt
+ I
Э2а„
2 Эг*
=/-
Юр12
4ю0Ц
зе/„
ВТ
д 1 д
-— + — — )2в=к2цх1т(а,„+ ехр/в),
О! О?
(3)
(4)
е12=0=е0 е [0,2тс], (-
Э 1 Э ч Ът т
где ь^йЫйт - групповая скорость электромагнитной волны, - пара-
метр дисперсионного расплывания, ае - параметр связи, пропорциональный
1 Здесь и далее ондуляторное поле рассматривается как волна накачки, имеющая нулевую
частоту.
осцилляторной скорости частиц, р. - параметр инерционной группировки, ат= =еАт/(тес2у0), к=а/с, 0=соо/- 1г0г - \о.±с1г - фаза электрона в поле синхронной волны, 5=(ю0- А0Оц " П±о)/г>|| " начальная расстройка синхронизма на несущей частоте, юр=(47геро/ме)1/2 - плазменная частота, сор1=(йр/уо1/2, р0 - невозмущенная плотность
2к
электронного пучка, /и=1/я|ехр(-/8)с/0о -амплитуда первой гармоники ВЧ-тока.
При записи (2)-(4) полагали, что изменения энергии Е=тес2у электронов малы, а также пренебрегали ближним кулоновским взаимодействием частиц. В таких предположениях уравнения движения частиц в различных модификациях ЛСЭ могут быть с точностью до величин параметров связи и группировки (см. [3-5]) приведены к универсальному виду (4).
Пучок Ат~, распространяющийся в отрицательном направлении оси г, не взаимодействует с электронами; его эволюция описывается уравнением
Эа ,„"
дг ¡V Э?
1 да,„- V Э2а„,"
+ I----=0.
Э/2
(5)
Далее уравнения (2), (5) следует проинтегрировать вдоль характеристик + Принимая во внимание граничные условия на зеркалах
ат (1,0 =Я2а т+(1,0, ат+1+(0,0=Я1ат- (0,0,
(6)
где I - расстояние между зеркалами, получим, что значения амплитуды волны на двух соседних проходах по резонатору связаны соотношением
. Им , д2а,
ат+1+- ат+ = - (1-ЗДК+ + —- (й^г
ЭС2
Э2а„" , юп.2 [ . , ч + ^ \kIJZ. (7)
ЭС2 4ю0авг о
При выводе (7) мы считали, что дисперсионное расплывание волнового пучка на одном проходе мало
г Э2ат± , Э2а,„± ] -г-— с/г в/
0 эс2
эр
Допустим теперь, что модули коэффициентов отражения близки к единице /?1>2 = 1. В таких условиях изменения амплитуды огибающей сигнала на одном проходе по резонатору мало: 1(ат+1+- ат+)/ат+1«1, что позволяет перейти от дискретной переменной т к медленному времени т, выбрав за единицу отсчета времени период обращения волны по резонатору Тк =21/В результате из уравнения в конечных разностях (7) получаем дифференциальное уравнение [4,28, 29]
/
Эа ~ЭГ
соо Ы)„г Э2а /ах, .2 а+г-
26
2 ЭС2 8юо/ о
1ае/щ^г,
(8)
где а=а(х,0, 1 -/?1/?2)) - добротность резонатора. Уравнение (8)
следует дополнить начальными и граничными условиями
а1т=о=ао(С),
а(т,0=а(т,С+Га).
Условие цикличности (10) позволяет разложить решение в ряд Фурье
(9) (10)
ос(т,С) =2 а,(т)ехр(- /Д9С), (11)
q= -°о
где Д?=(2п/Тк)ц. Амплитуду каждой из гармоник можно рассматривать как амплитуду моды резонатора, возбуждаемой электронным потоком.
Уравнения движения электронов (3) в переменных г, 1=1- г!х)р принимают
вид
Э 1 1 э
(— + (---) — )20 = £2цае1т(аехр/9). (12)
дг vw 61;
При этом малая величина изменения амплитуды огибающей сигнала на одном проходе позволяет не учитывать в уравнениях движения (12) явной зависимости амплитуды от продольной координаты. Уравнения (8), (12) совместно с начальными и граничными условиями (9), (11) образуют самосогласованную систему уравнений, описывающую эволюцию поля в ЛСЭ-генераторах.
В режиме одномодовой генерации, представляя решение (9) и (12) в виде
а=а9(т)ехр(-/А,С), + Д^,
приведем эти уравнения к форме
ах 2() 2 8со0/ о
= А:2цае1ш(а9схргтЭ), (13)
г , Эв
£1г=0 =Л0е [0,2л], — 12_0 =5 - Аг
Уравнения (13) совпадают с хорошо известными уравнениями, описывающими установление колебаний в одномодовых резонансных СВЧ-генераторах [31].
2. Режим стационарной инжекции электронов
В этом разделе анализируются многочастотные процессы в ЛСЭ-генераторах, на вход которых поступает стационарный электронный поток. В п.2.1 в рамках клистронной аппроксимации пространства взаимодействия генератора аналитически исследована устойчивость рабочей моды и найдены основные закономерности перехода к многочастотным автомодуляционным режимам генерации. В п.2.2 эти результаты подтверждены численным моделированием в рамках более универсальной и реалистичной модели ЛСЭ, в которой предполагается, что взаимодействие электронов с излучением непрерывно на всей длине генератора.
2.1. Анализ устойчивости режима одномодовой генерации
Для определения причин и условий возникновения автомодуляции полезен анализ режима стационарной одномодовой генерации на устойчивость. Аналитически такой анализ может быть легко проведен в рамках клистронной аппроксимации пространства взаимодействия, то есть в предположении, что взаимодействие электронов с волной происходит в двух коротких зазорах, разделенных областью свободного дрейфа частиц [32,33]. На практике такая
конфигурация может быть, например, реализована в убитронах с ондулятором, состоящим из двух коротких секций, которые помещены в начале и в конце пространства взаимодействия. Клистронной аппроксимации соответствует параметр связи в виде двух дельта-функций
ас=ае1б(2) + 1Ат). (14)
Будем далее также считать, что величина параметра связи в выходной секции значительно превосходит соответствующую величину в модулирующей секции ге1«яе2. При параметре связи, заданном соотношением (14), уравнения движения электронов легко интегрируются и система уравнений (8), (12) редуцируется к одному уравнению с запаздывающим аргументом [34]
да" 1 дГ Э2а~ ».а*,,,
—— +-а"-----= -¿XI— /х(1а1) ¡---г- и. (15)
Эг 20 2 Э£~2 1а1 и 4 >
где Х=(цж1ге2(йрх2)/(16с2) - параметр, пропорциональный току пучка, - функция Бесселя,
1/аег)"2г>, Ьйг=к16т.
При записи уравнения (15) мы выбрали в качестве несущей частоту, на которой расстройка синхронизма 8 равна нулю. КПД ЛСЭ в -рамках клистронной модели определяется соотношениями
1 } оГ*
Т1 =к3^х1х21Лтцп, Т1П = — дафа{-—, ]. (16)
ТКо !аI
Для определения характеристик стационарного режима генерации представим решение уравнения (15) в виде а*=а0ехрг(^2т*- Д?£"). Тогда для амплитуды стационарных колебаний и для электронной перестройки частоты имеем
1/2(2 = -X", (17а)
' £2+(иД,2)/2=-х', (176)
где
Х=Х'- а'-Хехр/ФЛЫ/оо (18)
- комплексная восприимчивость, вносимая электронным потоком в резонатор, Ф=А9/(1г - угол пролета электронов относительно рабочей моды. В режиме малого сигнала из соотношения (17а) получаем стартовые условия
<2Х8ШФ>1. (19)
Будем называть далее область углов пролета (область расстроек синхронизма), где условие (19) выполнено, областью самовозбуждения. В случае малых превышений тока пучка над минимальным пороговым значением <2Х=1.0 из (17а) получаем следующее соотношение между амплитудой моды и ее положением в зоне самовозбуждения (рис.2)
Оо2?=4[(АФ*)2- (ДФ)2], (20)
где ДФ* - ширина зоны самовозбуждения, а АФ - сдвиг частоты рабочей моды
i i м ГН
lililí
гттгтп
г
Рис. 2. Активная часть электронной восприимчивости %" как функция угла пролета рабочей моды Ф и процессы распада при умеренном (я) и сильном (б) превышении тока пучка над стартовым значением; 2ЛФ* - полоса самовозбуждения, Ф5 - полоса устойчивой одномодовой генерации (заштрихована), ДФ - сдвиг по углу пролета частоты рабочей моды от центра полосы самовозбуждения Ф = к/2, ¥ - сдвиг по углу пролета частоты моды-сателлита от частоты рабочей моды
относительно центра (Ф=я/2) указанной зоны. Нормированный КПД в одномодовом режиме генерации равен
rin=sin<í>ao/i(ao).
Максимальный КПД rjnmax=1.25sinO достигается при (X<2)crisnví>=2.25 (рис. 3).
Линеаризуя уравнение (15) вблизи стационарного состояния (17), то есть представляя а=(оо+е)ехр/(£2Г- Д9С"), где leí« 1, для малых возмущений получим
Эе ~дг
■ + iílz +
2Q-- 0'\ - "Э£г)2е= (21)
где р=Х71(а0)/(2а0), 5= -^'(а0)а0их(с1о). Заметим, что в условиях достижения максимального КПД параметр 5 = 1 (см.рис.З). Первый член в правой части уравнения (21) описывает изменение усиления затравочной волны (моды) в присутствии сильной волны. Второй член в отличие от первого существен лишь тогда, когда спектр мод резонатора близок к эквидистантному. Этот член ответствен за четырехквантовый параметрический процесс, состоящий в том, что два кванта интенсивной волны распадаются на два кванта волн сателлитов
2со = (о+ + ю.,
(22)
частоты которых <в± = а>±Д(о равно отстоят от частоты основной волны со. Соответственно, представляя решение (21)в виде
Рис. 3. Нормализованный КПД Т1П и параметр 5 для рабочей моды, расположенной в центре полосы самовозбуждения, как функции параметра ОХ., характеризующего превышение над порогом генерации; (£ХХ)сг1 - порог возникновения автомодуляции при слабой дисперсии, (£?Х)сг2 - тот же порог при сильной дисперсии
е=е+(г)ехр(-/Д"С~) + е. (г)ехр(/Д~^*),
(23)
где величина Л~ =Дсо/соо(с/уц - с/ь&) характеризует отстройку частот сателлитов
от частоты основной моды, дня амплитуд волн сателлитов приходим к уравнению связанных волн
с/е. е+ го*
+ 20~+ /Пе+ + (Л?+л~)2е+= ~[р ехР'(ф+*ЖМе+-(1+*)еЛ
(24)
с/е е /и"
где Ч?=А~Ь&. Решение (24) очевидно следует искать в виде е+,е.*~ехр(Гт"). В результате для инкремента распадной неустойчивости ИеГ получим следующее дисперсионное уравнение
[Г+(20)" 1 + /£2 + /(м-/2)(Д9+Д~)2 + /р( 1 -фхр/(Ф+Ч0]*
*[Г+(2<2)-1 - /£2 - ;(\)72)(Д?- Д~)2 - /р(1- 5)ехр(-/(Ф - ^))]=/?2(1+5)2ехр(2№). (25)
Проанализируем вначале случай пренебрежимо малой дисперсии г>"<< 1. Используя метод Д-разбиения, найдем, что волны-сателлиты нарастают (ЫеГ>0), если вьшолнен критерий
#
(1+£)(1+.у - 2gsУ■ > 8т2Ф(2-£- ^)(3 + 4^), ^совУ. (26)
В случае малых превышений над порогом генерации, когда значение амплитуды рабочей моды определяется соотношением (20), полученный критерий приобретает вид
Зср2- \|/2/2 > 1, (27)
где ф=ДФ/ДФ* и у =ХР/ДФ*. Согласно (27), рабочая мода неустойчива, когда ее частота достаточно сильно сдвинута от центра полосы самовозбуждения: <р>3"1/2 (рис.4, а). Моды, расположенные вблизи центра полосы, устойчивы.
Центральная мода теряет устойчивость при значительных превышениях тока над стартовым значением. Так для рабочей моды, расположенной строго в центре полосы (Ф=я/2), критерий неустойчивости имеет вид
«хю¥<-1. (28)
При ^ >1 эта мода распадается на два сателлита, расположенных в полосе поглощения Ч* « ± я (рис.4, б). Подчеркнем, что центральная мода теряет устойчивость точно в тот момент, коща ток пучка начинает превышать значение: <2Хсг1=2.25, при котором достигается максимум КПД в режиме одночастотной стационарной генерации.
В случае сильной дисперсии \Г» 1 и сильно неэквидистантном спектре мод резонатора взаимодействие между волнами - сателлитами пренебрежимо мало. При этом система уравнений (24) распадается на два независимых уравнения, и инкремент нарастания возмущений определяется соотношением
ИеГ=/?( 1 - ¿)яп(Ф+ч9 -1 /20. (29)
При малых превышениях тока пучка над стартовым значением, как и в случае слабой дисперсии моды, расположенные вблизи центра полосы самовозбуждения, устойчивы, а неустойчивость развивается только для мод
-1.0
2.0 1.0
"оо
Л 1 -1.0
-2.0
Рис. 4. Границы областей стабильности одномодового режима генерации: а - при малом превышении тока пучка над стартовым значением и произвольном положении рабочей моды в полосе самовозбуждения; б - при сильном превышении тока пучка над стартовым значением, когда частота рабочей моды лежит в центре полосы самовозбуждения. Кривые а - границы области стабильности при слабой ( *>«1) дисперсии, кривые в - те" же границы при сильной (1)»1) дисперсии
достаточно сильно сдвинутых от центра: <р>2"1/2. Нарастают моды с углами пролета, лежащими, согласно (29), в интервале (см.рис. 4, а)
Ч1<У<Щ, где \)/1>2=-ф± (2ср2-1)1Д. (30)
Для моды в центре полосы критерий неустойчивости приобретает вид (ср.с (28))
(5 - 1)С08ЧЧ-2
(31)
В этом случае боковые сателлиты 0¥=±л) возбуждаются при 5>3 (£Ж>4.15 ), то есть при больших превышениях над порогом генерации, чем при слабой дисперсии (см.рис.4, б).
Отсутствие устойчивых одночастотных решений означает, что многочастотные автомодуляционные режимы генерации должны реализоваться в случае больших превышений тока пучка над стартовым значением (0Х>2.25 в случае слабой и £?Х>4.15 в случае сильной дисперсии). Характерный период автомодуляции выходного излучения может быть определен по отстройке частоты волн сателлитов и основной, питающей их волны
Тэм
(32)
Очевидно, автомодуляция возникает, когда точка, соответствующая стационарному режиму генерации лежит на падающем участке зависимости высокочастотного тока пучка от,амплитуды поля: ¿1/щ1/с/ао=//(оо)<0. Таким образом, одна из причин возникновения автомодуляции - перегруппировка электронов. Другая причина - влияние значения амплитуды волны на одном участке огибающей на соседние участки [11,34]. Такое влияние обусловлено проскальзыванием волны относительно электронов из-за отличия групповой скорости волны от поступательной скорости частиц. Характерный временной масштаб этого влияния Тсот~1йх(1!ь\\ -1/авг) .
Рис. 5. Пространственно-временная диаграмма, поясняющая механизм автомодуляции в ЛСЭ -генераторах
Механизм автомодуляции может быть проиллюстрирован с помощью рис.5, на котором на плоскости показаны характеристики волны £=const и электронов £=z(l/t>n -1 /vgr) + const. Предположим, что амплитуда поля превысила стационарный уровень а0 на характеристике £ =0. Если с1\1а\!с1щ< 0, тогда амплитуда высокочастотного тока уменьшится в точке В. Соответственно, поле на характеристике £,=ТС0[ уменьшится, приводя к увеличению амплитуда ВЧ тока в точке С и возрастанию амплитуды поля на характеристике С=2Тсот. Далее процесс повторяется с периодом, определяемым соотношением (32).
. 2.2. Численное моделирование многочас-Готных процессов
Описанная выше модуляционная неустойчивость рабочей моды, расположенной в центре полосы усиления, должна приводить к реализации автомодуляционных режимов генерации. Исследуем здесь конкуренцию мод и развитие автомодуляционных процессов в рамках модели непрерывного взаимодействия электронного пучка с электромагнитной волной на всей длине генератора ае=сош1. Пренебрегая дисперсией (и=0) и уменьшая число независимых параметров путем замены переменных
х'=(соо/2<2)т, г'=кгО, С'=ю0С(с/г>ц - с/ь^'Ю, £=(со0!с)Ю, а=аецосО"2,
8'=(с/то)8£>-1, ТК'=щТц(с/г>1' - с/ьРУЮ,
о=[ж2цшр12С/(4(1-/г1й2)ш02гд]1/з,
опуская штрихи, приведем систему уравнений (8), (12) к виду
ь
да/дт + аы\ (33)
о
(д/ЬХ + Э/Э£)20 =1ш(аехрг'9), (34)
01^ - е0е [о,2я], (э/эг+ э/эс)е12=0 =6,
КПД генератора в новых переменных задается соотношениями
£) 1 2п
Л = -гг- Л~> л ~ = к 1 яВДЭ/Эг+Э/Эф- 5].
1 1) 2тс7К о о
2 Тсс
С
~уЪ
у
/
/
/
/
/
г У
~7
/
у
у
L dr г
Характер режима генерации определяется двумя параметрами: приведенной длиной генератора Ь и временем обращения волны по резонатору Тк. Начальная расстройка синхронизма на несущей частоте существенного влияния на режим генерации не оказывает, поэтому далее положим 5 = 0.
В использованных переменных стартовое условие имеет вид
1=£3(- х"),
Х"(Ф) =
е{ (1-созФ)
£Й>
ф2
Минимальная стартовая длина генератора Ьт1п=1.97 достигается при угле пролета Ф=А^=0.8я. В режиме одномодовой генерации максимум приведенного КПД реализуется при ¿=3.2, то есть примерно при троекратном превышении тока пучка над стартовьм значением. Величина 2п/Гк определяет межмодовое расстояние. Пример расположения мод внутри полосы активного вещества показан на рис. 6, б. В отсутствие дисперсии спектр мод эквидистантен.
При численном моделировании уравнений (33), (34) применялась алгоритмическая схема, подробно описанная в [35]. Для задания начальных условий и интерпретации результатов использовалось разложение Фурье (И), позволяющее определить модовый состав излучения в заданный момент времени. Как показывают расчеты, по мере увеличения тока пучка (параметра V) происходит последовательная смена режимов генерации.
тг
\ая\ 1.0
0.5
0.0
✓ тг
1
1 1 / X / / / /
1 - / / _
1 / и/ / ___4
/7 / 2 _________ 1 - 1 ""
1.0
0.5
0.0
4(0)1 ХУ~\ 0.5 \ I 1 I 1 \ 1ПШ
-6-4 -2 0 -0.5 V 2 4 4 /
4.0
8.0
12.0
б
Рис. 6. Установление режима одночастотной одпомодовой генерации в ЛСЭ с непрерывной инжекцией электронного пучка: а - эволюция амплитуд мод и приведенного КПД; б - начальные амплитуды мод и их расположение в полосе самовозбуждения (Ь = 2.8,Тц = 17.0)
1. Автоколебания возбуждаются при Ь> 1.97, и в интервале 1.97 < Ь< 3.2 устанавливается стационарная генерация на одной из мод, расположенной вблизи центра полосы усиления (см. рис.6). Этот результат справедлив при любом числе мод (М=ТК/Ь), попадающих в полосу усиления. Моды, расположенные на краю полосы усиления, в соответствии с приведенным в п.2.1 анализом неустойчивы. На рис.7 показан процесс распада такой моды, в результате которого рождается и выходит на стационарный уровень генерации мода, расположенная вблизи центра полосы. В целом картина конкуренции мод в этой области параметров аналогична имеющей место при несинхронном энергетическом взаимодействии мод [10,17,20], когда, черпая энергию от общего источника, моды взаимно подавляют друг друга.
Следует однако заметить, что с увеличением числа мод в полосе и соответственно с уменьшением расстояния между модами, время установления режима одночастотной одномодовой генерации возрастает (см. также [36,37]). Это обусловлено тем, что по мере приближения частоты паразитной моды к частоте
УД 0.5 к м м м\ - I I I м и
-6-4-2 0 -0.5
Рис. 7. а - распад моды q=l, расположенной на краю полосы усиления, и установление генерации на моде <7 = 3, находящейся в центре полосы; б - начальная амплитуда моды <р1 соответствует стационарному значению в рамках одномодовой теории, «затравочные» возмущения заданы для мод д = 0 и ц = 2 (I = 2.8, Гк = 17.0)
рабочей моды уменьшается декремент затухания паразитной моды. В частности, на рис.6 моды с номерами <?=2;4, расположенные вблизи рабочей моды, затухают медленнее, чем более отдаленные по частоте моды = 1;5.
2. При умеренном превышении параметра Ь над величиной 3.2 между модами в полосе усиления сохраняются отношения конкуренции. Однако победившая в процессе конкуренции мода рождает ряд сателлитов, расположенных в зоне поглощения (рис.8). Разность между номерами гармоник, амплитуды которых отличны от нуля, совпадает с числом мод М в полосе усиления. В установившемся режиме амплитуды всех гармоник постоянны, а их фазы жестко связаны. Поэтому такой режим может быть назван режимом самосинхронизации мод (ср.с [15,16]). Огибающая сигнала в этом режиме представляет собой серию импульсов постоянной формы (их число равно М). Важно подчеркнуть, что при заданных значениях параметров Ь и £> существует несколько режимов периодической автомодуляции, которые можно идентифицировать по номеру основной моды, лежащей в полосе усиления.
Величина области превышений над порогом Ь, в которой реализуется рассмотренный выше режим, убывает с ростом числа мод в полосе усиления и при М > 4-г5 становится исчезающе малой.
агяя
1.0 |в<Н 0.5
Л"
А,
-2 0к
6 8 <7
4.0 8.0 £ 12.0
а б
Рис. 8. Режим периодической автомодуляции (самосинхронизации мод): а - стационарная форма огибающей волны; б - спектр огибающей (£ = 4.0, Гк = 19.0)
\а,\
0.25
м
ш
11
8 12 16 Ц
24.0 32.0
а
Рис. 9. Режим стохастической автомодуляции: а - эволюция амплитуд первых четырех мод и приведенного КПД; б - спектр огибающей при т = 30.0 (Ь = 4.8, Тк = 19.0)
3. Дальнейшее увеличение превышения над порогом Ь или увеличение числа мод М в полосе усиления приводит к тому, что амплитуды мод начинают меняться во времени сначала периодически, а затем (по мере роста Ь) и стохастически (рис. 9, а). Кроме того, происходит рождение (если даже в начальный момент времени их амплитуды были равны нулю) всех прочих волн, расположенных в полосе усиления вместе с их сателлитами. Соответственно эволюция огибающей приобретает сложный стохастический характер. В пределе очень больших превышений над порогом генерируемое излучение распределено в спектральном интервале, существенно превышающем полосу линейного усиления (рис. 9, б).
Интегральный КПД при этом может быть значительно выше, чем в режиме одномодовой генерации, поскольку электрон, провзаимодействовав с одной модой и отдав ей часть своей энергии, оказывается в синхронизме и попадает в область захвата следующей моды, которой соответствует меньшая фазовая скорость комбинационной волны и т.д. В результате процесс торможения электронов приобретает стохастический характер. Возникающая при этом диффузия частиц в пространстве скоростей аналогична рассмотренной в разделе 4. Подчеркнем, что описанный выше механизм роста КПД экспериментально наблюдался в [38], где в ЛСЭ инфракрасного диапазона по мере увеличения превышения над порогом был зарегистрирован значительный рост КПД, сопровождавшийся соответствующим уширением спектра излучения.
3. Режим периодической импульсной инжекции электронов
Развитый в разделах 1,2 подход к описанию многочастотных процессов в генераторах с высокодобротными резонаторами может быть использован и для описания ЛСЭ, запитываемых последовательностью коротких импульсов, длительность которых Тс существенно меньше периода их следования Тг В таких генераторах период обращения волны по резонатору Тц подбирается примерно равным периоду следования импульсов тока 7^. На спектральном языке это означает, что частота воздействующей на систему внешней силы (в данном случае частота модуляции тока) совпадает с межмодовым расстоянием, то есть выполнено необходимое условие вынужденной синхронизации мод [15].
В результате синхронизации в резонаторе образуется электромагнитный импульс с длительностью, близкой к длительности импульса тока. Оба импульса (электровный и электромагнитный) вместе бегут по резонатору, лишь немного смещаясь друг относительно друга из-за отличия групповой скорости световой волны от поступательной скорости электронов. Достигнув правого зеркала, электронный импульс покидает резонатор, а электромагнитный отражается и возвращается назад, достигая левого зеркала как раз в момент прихода следующего импульса тока (рис. 10). Таким образом, с точки зрения
ч.
1=0
1= 7Ь
¡'¡к 1 4 Л-
7к
ЗТя
г»
электромагнитного импульса электронный пучок практически эквивалентен сплошному. Описанный выше режим генерации был впервые предложен и реализован в эксперименте [1] и затем использовался в большинстве ЛСЭ оптического и инфракрасного диапазонов (см.обзор [2]), в которых источниками питания служили линейные ускорители, накопительные кольца и т.п.
Процесс взаимодействия электронного и электромагнитного импульсов может быть описан с помощью системы уравнений (33)-(34), в которой, принимая во внимание расстройку между периодом следования импульсов тока и периодом обращения волны по резонатору, задаваемую параметром
а=20£(с/у„ - -ТЛ)/ТК,
удобно перейти к новой независимой переменной ох. В новых
переменных уравнения принимают вид
Рис. 10. Диаграмма, поясняющая синхронизацию светового и электронных импульсов в ЛСЭ, запитываемом периодической последовательностью электронных импульсов; А/ - временной сдвиг, обусловленный различием периодов следования импульсов тока т, и времепи обращения светового импульса по резонатору
Рис. 11. Режим периодической импульсной инжекции электронов: характеристики уравнения для волны на плоскости х, £ и характеристики уравнений движения частиц на плоскости Ъ, с,, а также стационарная форма электромагнитного импульса Волна и частицы взаимодействуют на интервале 0< % < Гс' + Ь
да/дх - да/д$ + а = /| (35)
о
(Э/Эг + Э/Э£)28 = 1т(аехр/0). (36)
Здесь введена функция С(^), характеризующая форму токового импульса. В проводимых ниже расчетах предполагалось, что импульс имеет прямоугольную форму и задан на интервале 0<^<ГС', где 7'с'=ЮоХс(с/г)ц -- с/ЬрУЮ. Это позволяет записать граничное условие к уравнению (35) в виде2
Важно подчеркнуть [22,3], что при периодической инжекции электронов генерация возникает только тогда, когда период следования импульсов тока несколько превышает период обращения волны о>0 (предполагается Ьц<г)р ). При этом, как видно из рис.11, на котором приведены характеристики уравнений (35) и (36), возмущения переносятся электронами и волной в противоположных направлениях, то есть реализуется механизм обратной связи, аналогичный имеющему место в ЛОВ (ср. с [34]).
2 Поскольку длина электромагнитного импульса существенно короче расстояния между зеркалами, условие цикличности (10) можно считать выполненным автоматически.
На рис. 1.2 на плоскости параметров Ь и Те' показаны зона генерации импульсов стационарной формы и зона, в которой имеет место автомодуляция формы импульса. Заметим, что при больших длительностях импульса тока Тс' стартовое значение параметра примерно совпадает со значением, характеризующим порог одномодовой генерации в условиях непрерывной инжекции = 2.0). Рис. 13
иллюстрирует пространственно-
временную эволюцию формы импульса в различных режимах генерации, возникающих по мере превышения над порогом. Видно, что при больших превышениях имеет место
стохастическая автомодуляция формы импульса.
Рис. 12. Границы области самовозбуждения (сплошная линия) и области автомодуляции
(пунктир) о = 2.0
на плоскости параметров ь, тс
Рис. 13. Эволюция электромагнитного импульса в ЛСЭ в режиме импульсной инжекции: а -установление стационарной формы светового импульса (¿=2.4, тс' = 4.0, 0=2.0); б - периодическая автомодуляция формы импульса, (¿=5.1, Гс'=4.0, 0=2.0); в - стохастическая автомодуляция формы импульса (ь = 8.0, тс' = 4.8, о =0.8)
4. Многомодовая генерация в ЛСЭ с неоднородно-уширенной линией активного вещества
Выше были рассмотрены многочастотные процессы в ЛСЭ, работающих в режиме однородного уширения линии активного вещества, который реализуется, когда разброс параметров электронного потока достаточно мал. Исследуем теперь многомодовую генерацию в ЛСЭ с неоднородно-уширенной линией активного вещества, в которых дисперсия параметров играет существенную роль.
Обобщение уравнений (33)-(34) на случай пучка, имеющего конечный разброс параметров, в принципе, элементарно: пучок нужно разбить на отдельные фракции, отличающиеся расстройкой синхронизма 5, а правую часть уравнения для амплитуды поля (33) дополнить интегралом по 5 с соответствующей весовой функцией, определяемой функцией распределения электронов на входе в резонатор. Однако по мере увеличения дисперсии энергий электронов и, как следствие этого, перехода к кинетическому режиму взаимодействия, число электронных фракций будет возрастать, что в конечном итоге приводит к необходимости использования приближенных методов.
Заметим, что в кинетическом режиме взаимодействия источником энергии для каждой моды является своя электронная фракция, резонансная только этой моде, и нелинейная связь мод слаба. Поэтому можно считать, что возбуждение каждой из мод происходит независимо от остальных и описывается линейным инкрементом, а ограничение усиления возникает из-за деформации усредненной
функции распределения. Подобный подход в физике плазмы принято называть квазилинейным (см., например, [12-14]). Модифицируем здесь указанный подход применительно к возбуждению высокодобротных резонаторов размытыми по энергиям РЭП [21], рассматривая для определенности ЛСЭ, основанные на вынужденном рассеянии волн (вынужденном ондуляторном излучении).
Пусть поле накачки представляет собой линейно поляризованную волну: А|=Ке(у0Л1схрг(сй|?+/г|2)}. Поле, возбуждаемое электронным потоком в высокодобротном квазиоптическом резонаторе, можно записать в виде суммы мод
А, =Ке[2ХуИ^п(/г92)ехр(/йу)], (37)
<7=1
частоты которых образуют эквидистантный спектр (оч=щ/1, а фазы в пренебрежении нелинейной связью мод случайны.
Усредненное движение релятивистских электронов в полях сигнала и накачки определяется пондеромоторной силой
е2 Э ~
=- -т-— V- ИЯе^^ехр/'((0С9/- 1гсг1г)], (38)
2тес2у ог <г=1
представляющей собой ансамбль комбинационных волн с частотами (осч=и>ч- со, и волновыми числами АС9=А9+ /г,. Это движение может быть описано с помощью кинетического уравнения для одномерной функции распределения Д/,2,рг) по продольным координатам и импульсам
Э/ Э/ Э/
эг+1,'э;+/^га (39)
Амплитуды мод могут быть найдены из уравнений возбуждения резонаторов
(40)
л йч~тч о
2к
где;Ш¥=1/п|/^ехр(-/(со¥Г- кчг))с1(й^ - амплитуда Фурье-гармоники тока,}у = -р0/ру/фг,
о
у>=(е/с)Ке[А1ехрг'(ю,Г+/г^)] - осцилляторная скорость, приобретаемая электроном в поле накачки.
Определим сначала инкременты возбуждения мод. С этой целью линеаризуем кинетическое уравнение (39), представив функцию распределения в виде
/=/0 + £ Ке[/(?ехр/(шс^- ксдг)], (41)
где/о - начальная функция распределения,
еу^ »
Подставляя (41), (42) в правую часть уравнения (40) и производя интегрирование по рг с учетом правил Ландау для обхода полюсов: ьг=ис<=(ося/ксч , приведем это уравнение к форме
где Г - инкремент кинетического режима неустойчивости
к юр2 еЛ\ Э/о
В линейном приближении инкремент Г не зависит от продольной координаты 2. Однако с учетом обратного воздействия волн на частицы, приводящего к деформации функции распределения, величина становится функцией г. В таком квазилинейном приближении под величиной/0 в (44) следует понимать усредненную величину, эволюция которой описывается диффузионным уравнением
Э/о Э «■
^-= " <Рг (45)
02 дрг 9=1
где угловые скобки означают операцию усреднения. В уравнении (45) производная Э/о/Эг исключена, поскольку в высокодобротном резонаторе с коэффициентами
отражения зеркал К12 ~ 1 характерное время изменения амплитуд мод Т0=Оя/а<!-1-Л1Л2) существенно превышает время пролета электронов через резонатор
Т=11ьг.
Вычисляя с помощью (38), (42) выражение в скобках, получим
1 ~ 1е*К2\Ач\2Щ2 Э/0
8 <,=1 (1гчхк- соС9)Е2 Эр,
Перейдем далее от суммирования к интегрированию по Н¥ принимая во внимание, что в интервале волновых чисел А/г содержится Л<7=(//тс)Л/г мод резонатора. Это позволяет привести уравнение (45) к виду
- Э/о Э (46)
Ъг Эр2 Эр.
где коэффициент диффузии
пе*ИМ>\2
£>= , ' Ж,
ф
8(с- ог)Е2
пропорционален спектральной плотности электромагнитного поля N^=(//71)^,12 в резонансной точке спектра ю=(сО| + г;гЛ,)/(с- .
Уравнения (43), (46) составляют самосогласованную систему квазилинейных уравнений многомодового генератора с высокодобротным резонатором. Эти уравнения должны быть дополнены начальными и граничными условиями
/ои=/о0, ид,=0=1^ о,
где /п° - невозмущенная функция распределения электронов на входе в резонатор, М'е,0 - начальный уровень электромагнитных шумов.
Рассмотрим далее некоторые характеристики стационарного режима генерации. Полагая в уравнениях (43), (46) Э/Э/=0 и считая пучок ультрарелятивистским т»1, после перехода к безразмерным переменным приведем эти уравнения к виду
1
1=П/(Г1— К, (47)
о
dF Э dF
К Эу
Эу
Здесь n=(7t/16)(mp2^2|Ai|27,0)/(me2c4Q1), =Ю; + h(Vz - частота осцилляций электронов в поле накачки; F-mecf - функция распределения электронов по энергиям, удовлетворяющая условию нормировки: ¡Fdy= 1; {¡=z/Z - нормированная длина; W(aA=(8e4\Ai\2Q.i2l2)/(m^cw)Wa - нормированная спектральная плотность электромагнитной энергии; со= 2y2i2j_.
Пренебрегая деформацией функции распределения, получим из уравнения (47) стартовое условие
dFo
1=П(7" 1 ——)1т=(ш/(2п1))1д (49)
Эу
для моды, частоте со которой соответствует резонансная энергия частиц: у= =(co/(2Q1))1/2. Очевидно, существует минимальное значение параметра П, при котором стартовое условие выполнено только для одной моды с максимальным инкрементом. Резонансная энергия частиц для этой моды лежит вблизи точки перегиба функции F°(y).
При n>nmin стационарная спектральная плотность энергии излучения отлична от нуля в конечном интервале значений (рис. 14). Чтобы определить этот интервал и соответствующий интервал резонансных энергий электронов, обозначим его нижнюю границу уъ верхнюю - у2 и ?ведем функцию
1
Ф(у)=| Felt,.
о
Эта функция характеризует число частиц с энергией у на всей длине пространства взаимодействия. Вне интервала У1<У<У2 коэффициент диффузии равен нулю, и функция распределения не деформируется. Следовательно, дня таких у: Ф(у)=р°(у). Внутри упомянутого интервала функция распределения находится из соотношения (47)
Ф=(П" !у" 2)/2 + const. (50)
Требуя непрерывности функции Ф(у) в точках у12 и принимая во внимание
вытекающее из условия сохранения числа частиц условие нормировки
¡фс1у=1, получим следующие уравнения для определения величин
НгО + П"1(у22/2-у,2/2) = F°( у2),
72 Ь
J[F°(Yi) + П"1 СУ2/ 2 -у,2/2)]с/у = J/70 (у) с/у. (51)
Решение уравнений (51) для начальной гауссовой функции распределения
1 /у- ус)2
Рис. 14. Начальная функция распределения F°(у) F°(у) =-. , ехрр- —--——]
и возмущенная функция распределения Ф(у) при Vя/ v°eYcJ
П/ус = 0.2. Штриховая линия - равновесная
спектральная плотность энергии электро- представлено на рис. 15. С увеличением магнитного поля Wm (у) (качественно)
0.5 Virc t5 ^
генерируемое (усиливаемое) излучение будет являться монохроматическим. Механизм повышения КПД обусловлен тем, что расширение спектра поля накачки даже при монохроматической сигнальной волне приводит к обогащению спектра фазовых скоростей, синхронных с электронами комбинационных волн. В результате, провзаимодействовав с одной спектральной составляющей ансамбля комбинационных волн и отдав ей часть энергии, электроны оказываются в синхронизме с другой составляющей и т.д. В таких условиях КПД оказывается пропорциональным ширине спектра поля накачки.
Для описания ЛСЭ с некогерентным (случайным) полем накачки можно, как в разделе 4, использовать квазилинейный подход, полагая, что усиление сигнальной волны определяется инкрементом кинетического режима взаимодействия, а эффекты насыщения обусловлены деформацией усредненной функции распределения электронов.
Исследуем здесь усилительный вариант ЛСЭ. Пусть монохроматическая сигнальная волна задана вектор-потенциалом
А8=Ке[уоДв(г)ехр/(юГ - 1п)\. (54)
Предполагая, что поле накачки состоит из ряда дискретных линий с некоррелированными между собой фазами
А, =Ке[£уД?ехр/(со^ + (55)
7=1
будем считать, что расстояние между отдельными гармониками 8И\ много меньше полной ширины пакета А/г,.
Усредненное движение релятивистских электронов в полях сигнала (54) и накачки (55) определяется пондеромоторной силой (38), где ссц=со - со1(?, /гс?=/г+/г1<?, и может быть описано с помощью кинетического уравнения (39) для одномерной функции распределения №,г,рг). Производя далее усреднение по комбинационным волнам, приходим к диффузионному уравнению
(56)
где коэффициент диффузии
Здесь /5=1А$|2со2/(8яс) - интенсивность сигнальной волны, 1=\Л1\2(л2/(8ксЬЬ,) спектральная интенсивность поля накачки. В предположении о фиксированной амплитуде поля накачки усиление сигнальной волны описывается уравнением
(II, к2е2<йр2
с12 тсс$ А о ю,2 др.,
~ (57)
г
В ультрарелятивистском пределе система уравнений (56)-(57) после перехода к безразмерным переменным приобретает вид
^ Э Л^о(у),
зг Зу Эу Й 2 (¡Й ! Эу
параметра П ширина интервала у, где спектральная плотность энергии излучения отлична от нуля, возрастает. Заметим также, что нижняя граница этого интервала у1 движется быстрее, чем верхняя у2.
Согласно соотношениям (47), (48), внутри интервала уе[уьу2] спектральная плотность энергии излучения определяется разностью функций распределения на входе и выходе резонатора и задается интегральным соотношением
1 7
ЮЛу) = — 1(ПС= 1 (52)
Пу5
Рис. 15. Зона резонансных энергий электронов (заштрихована), которой соответствует отличная от нуля спектральная плотность энергии электромагнитного поля, как функция параметра П/ус, характеризующего превышение над порогом генерации
Как следует из (52), спектральная
плотность сначала возрастает по у (и соответственно по ш), достигает максимума в области /7(^=1) = а затем спадает до нуля при у=у2-
Определить вид функции распределения на выходе из резонатора можно в предельном случае очень больших превышений над порогом П»1. В этом случае образующаяся в результате диффузии функция распределения должна иметь вид плато и может быть приближенно задана в виде
^/(ъ-кО-
Вычисляя изменение потока кинетической энергии частиц
К=(р0/е)}/и1.с3у/7й?у,
(53)
найдем, что интегральный КПД ЛСЭ т) =1-К12=(/К12=0 при у1«у2 стремится к 50%. Разумеется, излучение в этом случае распределено в очень широком спектральном интервале: Дш/и-Г. При конечных превышениях спектр излучения уже, но одновременно меньше и интегральный КПД.
5. Использование режима стохастического торможения электронов в поле некогерентной накачки для увеличения эффективности ЛСЭ
Выше в разделах 2 и 4 было показано, что реализация режима стохастического торможения электронов позволяет увеличить эффективность энергоотдачи электронного пучка. Однако это достигается за счет расширения спектра и ухудшения когерентных свойств генерируемого излучения. В ЛСЭ, основанных на вынужденном рассеянии волн и вынужденном ондуляторном излучении3, существует возможность реализации режима стохастического торможения и повышения КПД при использовании некогерентной накачки или стохастических ондуляторов. При электродинамическом разрежении спектра мод в генераторах или подаче входного сигнала от внешнего источника в усилителях
3 В ЛСЭ черенковского типа с этой же целью могут использоваться нерегулярные замедляющие системы.
где 0~=(4Р5Р11г2у2)/к, - функции распределения электронов по энергиям,
2=(со/с)г, Р&=Ь\21Р* - мощность сигнальной волны в единицах Р,=те2с5/е2,
проходящая через площадь, равную X2, Р, - спектральная мощность
накачки в той же нормировке.
Для проведения численного моделирования начальная функция распределения электронов по энергиям задавалась следующим выражением:
=2/аесо82[(«-1 )я/ае]
в интервале [- ое/2, oJ'2\ в окрестности энергии центральной фракции ус, где ое= =Ду/ус - ширина разброса электронов по энергиям, и=у/ус- нормированная энергетическая переменная.
Спектральное распределение поля накачки аппроксимировалось гауссовой формой линии /(=/; '"'^(/г,), где I,- интегральная интенсивность накачки,
описывающую
^¡И/^ЛАОехр(-(А, - ЛЕЖАЛО2).
В используемых переменных нормированную функцию, спектральное распределение накачки, следует переписать в виде
где <^=Л/г|//г, - относительная ширина спектра накачки, а величина щ=у/ус (у;=/?/(4/гк))1/2) характеризует центр спектра поля накачки в пространстве энергий.
С учетом закона сохранения энергии в системе пучок - сигнальная волна уравнения (58) удобно привести к виду, содержащему минимальное число независимых параметров
Э ^
Э и
Э э^
[(л+^МФ*^—],
Эн
(59)
где Ус Х21^г)/(4пс2Р,), />;о=(2Р8Осо2)/(л:юр2 ус), ц=1-\иРс1и - электронный
КПД. о
Результаты численного моделирования уравнений (59), представленные на рис. 16, 17, показывают, что при увеличении ширины спектра накачки сть эффективность трансформации энергии возрастает и при относительной ширине спектра накачки 40% максимальный КПД достигает 15%. Вместе с тем увеличение ширины спектра накачки приводит к падению инкрементов и, соответственно, увеличивается длина пространства взаимодействия, на которой
л,%
ю.о -
л,%
10.0
10.0 20.0 30.0 ,40.0 'С,
Рис. 16. Зависимости КПД от длины области усиления при различных ширинах спектра стохастической накачки О; и начальных дисперсий частиц по энергиям: а - ое =15%, б - сте = 30%
реализуется максимальный КПД. Важно подчеркнуть, что рассматриваемый метод повышения КПД в широких пределах некритичен к качеству электронных пучков. Как видно из сравнения рис.16, а и б максимальный КПД практически не зависит от параметра ае, характеризующего дисперсию параметров электронного потока. Действительно, насыщение усиления обусловлено образованием плато на функции распределения электронов (см.рис. 17), ширина которого и, следовательно, максимальный КПД определяются шириной спектра накачки, а также расстоянием (в пространстве энергий) между центром линии накачки и центром функции распределения частиц. Интенсивность накачки и плотность электронного пучка оказывают влияние в основном на
использовании в качестве источников накачки СВЧ-генераторов или оптических лазеров управление шириной спектра накачки представляет собой достаточно трудную в техническом отношении задачу. По-видимому, значительно проще расширение спектра накачки может быть обеспечено в убитронах с нерегулярными (стохастическими) ондуляторами, в которых расстояние между магнитами есть случайная величина. Рассмотренный здесь метод повышения КПД можно рассматривать как альтернативу режиму синхронного торможения [39, 40], когда период ондулятора меняется таким расчетом, чтобы электроны длительное время оставались захваченными комбинационной волной. Захват значительной доли электронного пучка такой волной возможен, если энергетический разброс электронов достаточно мал в масштабе амплитуды этой волны. В убитроне со стохастическим ондулятором разброс параметров пучка может быть существенно выше, что позволяет считать такой метод повышения КПД перспективньЛ! при использовании в качестве источников питания ЛСЭ сильноточных релятивистских электронных пучков.
Библиографический список
1. Deacon D.A.G., Elias L.R., Madey J.M.J., et.al I/ First operation of free-electron laser // Phys.Rev.Lett. 1977. Vol.38, № 16. P.892.
2. Roberson C.W., Sprangle P. A review of free-electron lasers // Phys. Fluids. 1989. Vol.1, №1.P.3.
3. Bratman V.L., Ginzburg N.S., Petelin M.I. Common properties of free electron lasers // Optics Commun. 1979. Vol.30, № 3. P.409.
4. Братман В.Л., Гинзбург H.C., Петелин М.И. Лазеры на свободных электронах: перспективы продвижения классических электронных генераторов в коротковолновые диапазоны // Изв. АН СССР. Сер. физ. 1980. Т.44, № 8. С.1593.
5. Братман В.Л., Гинзбург Н.С., Ковалев Н.Ф., Нусиновин Г.С., Петелин М.И. Общие свойства коротковолновых приборор с длительной инерционной группировкой электронов // Релятивистская высокочастотная электроника. Горький: ИПФ АН СССР, 1979. Вып.1. С.249.
6. Gover A., Sprangle P. A unified theory of magnetic brems strahlung, electro-
Рис. 17. Эволюция функции распределения электронов при относительной ширине начального энергетического разброса частиц ое=15% и относительной ширине спектра накачки а; = 20%: 1 - С=0.0; 2- £=10.0; 3 -С= 12.0; 4 - £ = 18.0; 5 - £=20.0
инкремент и длину области рассеяния. Отметим в заключение, что при
static bremsstrahlung, Compton-Raman scattering, and Chcrenkov-Smith-Purcell free-electron laser // IEEE J. Quant. Electr. 1981. Vol.QE-17, № 8. P.1196.
7. Ruell D., Takens F. On the nature of turbulence// Comm. Math. Phys. 1971. Vol.20.P.167.
8. Feigenbaum M.J. Quantitative universality for a class of nonlinear transformation // J. Stat. Phys. 1978. Vol.19, № 1. P.25.
9. Рабинович М.И. Стохастические колебания и турбулентность // УФН. 1978. Т.13, № 1. С.123.
10. Рабинович М.И., Трубецков Д.И. Введение в теорию колебаний и волн. М.: Наука, 1984.
11. Кузнецов С. П. Сложная динамика генераторов с запаздывающей обратной связью // Изв. вузов.Сер. Радиофизика. 1982. Т.25, № 12. С.1410.
12. Кадомцев Б.В. Коллективные явления в плазме.М.: Наука, 1976.
13. Цытовпч В.II. Нелинейные эффекты в плазме. М.: Наука, 1973.
14. Шапиро В.Д., Шевченко В.И. Взаимодействие волна-частица в неравновесных средах // Изв. вузов. Сер. Радиофизика. 1976. Т.19, № 5-6. С.767.
15. Ярив А. Квантовая электроника. М.: Сов. Радио, 1980.
16. Файн В.М., Ханин Я.И. Квантовая радиофизика. М.: Сов. Радио, 1965.
17. Lamb W.E. Theory of optical masers // Phys.Rev.A. 1964. Vol.134. P.1429.
18. Беспалов В.И., Кубарев В.А., Пасманик Г.А. Вынужденное рэлеевское рассеяние света // Изв. вузов. Сер. Радиофизика. 1970. Т.13, № 10. С.1433.
19. Datolli G., Reniery A. The free-electron laser single-particale multimode classical theory // Nouvo Cimento. 1981. Vol.61 B, № 2. P.153.
20. Нусинович Г.С. Взаимодействие мод в лазерах на свободных электронах // Письма в ЖТФ. 1980. Т.6, № 14. С.848.
21. Ginzburg N.S., Shapiro М.А. Quasi-linear theory of multimode FELs with an inhomogeneous frequency broadening // Optics Commun. 1982. Vol.40, № 3. P.215.
22. Al-Abawi II., HopfF.A., Moore G.T., Scully M.O. Coherent transients in the free-electron laser: laser lethargy and coherence brightening // Optics Commun. 1979. Vol.30, № 2. P.235.
23. HopfF.A., Kufer T.G., Moore G.T., Scully M.O. Free-electron laser from a laser physics perspective // Free-Electron Generators of Coherent Radiation. Physics of Quantum Electronics. Addison-Wesley, 1980. Vol.7. P.31.
24. Colson W.B. Optical pulse evolution in the Stanford free-electron laser and in tapered wiggler // Free-Electron Generators of Coherent Radiation. Physics of Quantum Electronics. Addison-Wesley, 1982. Vol.9. P.457.
25. Colson W.B., Richardson J.L. Multimode theory of the free-electron laser oscillators // Phys.Rev.Lett. 1983. Vol.50, № 14. P.1059.
26. Colson W.B., Freedman R.A. Synchratron instability for long pulses in free-electron lasers // Optics Commun.1983.Vol.36, № 3.P.209.
27. Mani S.A., Korff D.A., Blimnell J. Resonator mode structure // Free-Electron Generators of Coherent Radiation. Physics of Quantum Electronics. Addison-Wesley, 1982. Vol.9. P.557.
28. Bogomolov Ya.L., Bratman V.L., Ginzburg N.S., Petelin M.I., Yunakovsky A.D. Nonstationary generation in free-electron lasers // Optics Commun. 1981. Vol.36, № 3. P.209.
29. Ginzburg N.S., Petelin M.I. Multifrequency generation in free-electron lasers with quasi-optical resonators // Int. J. of Electronics. 1985. Vol.59, № 3. P.291.
30. Гинзбург H.C. Об использовании в лазерах на свободных электронах некогерентной накачки // Письма в ЖТФ. 1984. Т.10, № 10. С.584.
31. Вайнштейн JI.A. Общая теория резонансных электронных автогенераторов // Электроника больших мощностей. М.: Наука, 1969. Вып.6. С.84.
32. Петелин М.И. Электронная селекция мод в гиротронах // Гиротрон. Горький: ИПФ АН СССР, 1981. С.77.
33. Ginzburg N.S., Petelin M.I., Sergeev A.S. On the mechanism of selfmodulation of radiation in free-electron lasers // Optics Commun. 1985. Vol.55, № 4. P.283.
34. Гинзбург Н.С., Кузнецов С. П. Периодическая и стохастическая автомодуляция излучения в электронных генераторах с распределенным взаимодействием // Релятивистская высокочастотная электроника. Горький: ИПФ АН СССР, 1981. Вып.2. С.101.
35. Bogomolov Ya.L., Yunakovsky A.D. Numerical simulation of nonstationary processes in free-electron lasers // J. Сотр. Phys. 1985. Vol.58, № 1. P.80.
36. Levush В., Antonsen T.M. Mode competition and suppression in free-electron laser oscillators // Phys. Fluids. 1989. Vol.1, № 5. P.1097.
37. Levush В., Antonsen T.M. Mode competition and control in free-electron laser oscillators // Phys.Rev.Lett. 1989. Vol.62, № 13. P.1448.
38. Chatс P. FEL theory and experiment in Bruyeres-le-Chatel // 16-th Int. FEL Conference. Stanford, 1994. Abstr.P.90.
39. Kroll N.M., Morton Ph.,L., Rosenbluth M.N. Enhanced energy extraction in free-electron lasers by means of adiabatic decrease of resonant energy // Free-Electron Generators of Coherent Radiation. Physics of Quantum Electronics. Addison-Wesley, 1980. Vol.7. P.113.
40. Sprangle P., Tang C.M., Manheumer W.M. Nonlinear theory of free-electron laser and efficiency enhancement // Phys. Rev. A. 1980. Vol.21, № 1. P.302.
Иститут прикладной физики РАН, Поступила в редакцию 2.11.94 г.
Нижний Новгород
MODE COMPETITION AND COOPERATION IN FREE-ELECTRON LASERS
N.S. Ginzburg, M.I. Petelin
Mode interaction in free-electron lasers operating under homogeneous and in-homogeneous broading of amplification line is reviewed. Analytical and numerical methods demonstrate that in the case of homogeneous broading with increase of electron current the following stages take place successively: regime a) single mode oscillations, b) periodical selfmodulation ( mode phase-locking), c) chaotic selfmodulation. It is important that in the last regime electron efficiency may essentially exceed efficiency of the single-mode operation regime. The enhancement is caused by the effect of stochastic particles" deceleration in the multi-frequency chaotic radiated field, when an electron, having lost some energy to one mode, is trapped in synchronism by another mode, which corresponds to a lower phase velocity. It is shown that in the case of the monochromatic signal wave similar mechanism of deceleration and efficiency enhancement may be realized using incoherent pump wave or stochastic undulators.
Гинзбург Наум Самуилович, 1952 г. рождения, ведущий научный сотрудник Института прикладной физики РАН, доктор физ.-мат. наук. Область интересов - лазеры на свободных электронах, мазеры на циклотронном резонансе, релятивистские гиротроны, многочастотные процессы в электронных СВЧ-приборах с распределенным взаимодействием, эффекты канализации и сверхизлучения, плазменно-пучковые неустойчивости. Лауреат премии им. Ленинского комсомола (1980). Автор более 120 статей и обзоров по указанной проблематике
Петелин Михаил Иванович, 1937 г. рождения, заведующий отделом релятивистской высокочастотной электроники Института прикладной физики РАН, профессор Нижегородского государственного университета. Область интересов - электроника больших мощностей и электродинамика. Автор фундаментальных работ по теории индуцированного излучения ансамблей неизохронных электронов-осцилляторов, приведших к созданию нового класса приборов СВЧ-мазеров на циклотронном резонансе. Лауреат государственной премии СССР (1967). Автор более 100 статей и обзоров, а также около 20 авторских свидетельств на изобретение.