CC BY
7ISSN 0868-5886
НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2015, том 25, № 2, c. 69-90
МАТЕМАТИЧЕСКИЕ МЕТОДЫ -
И МОДЕЛИРОВАНИЕ В ПРИБОРОСТРОЕНИИ
УДК 681.787.7
© А. С. Бердников, Н. К. Краснова
ДОСТАТОЧНЫЙ КРИТЕРИЙ УСТОЙЧИВОСТИ И КОМПАКТНОСТИ ПЛОСКИХ ИОННЫХ ПУЧКОВ В ТРЕХМЕРНЫХ ЭЛЕКТРИЧЕСКИХ И МАГНИТНЫХ ПОЛЯХ С ПЛОСКОСТЬЮ СИММЕТРИИ
Рассматривается проблема устойчивости пучков заряженных частиц, которая возникает при отклонении движения частиц от плоскости симметрии г = 0 электрического и/или магнитного поля. Для электрических полей показано, что если для электростатического потенциала и(х, у), заданного в плоскости симметрии г = 0, выполнено условие ихх + иуу < -е2 < 0, то это будет достаточным условием, чтобы при малом отклонении положительно заряженных частиц от плоскости симметрии г = 0 пучок, двигающийся в данном электрическом поле, оставался в ограниченной окрестности плоскости симметрии. Дополнительно возможно построить конструктивные оценки, насколько далеко от плоскости симметрии может удаляться траектория при заданном возмущении начальных условий. Указанный достаточный критерий устойчивости и компактности пучков может быть обобщен для магнитных полей с плоскостью симметрии и для совмещенных электростатических и магнитостатических полей с общей плоскостью симметрии.
Кл. сл.: электронная оптика, компактные пучки заряженных частиц, устойчивость движения, синтез оптимальных электронно-оптических систем
ПОСТАНОВКА ЗАДАЧИ
Существуют эффективные алгоритмы синтеза электростатических и магнитостатических полей, обеспечивающих желаемое двумерное движение заряженных частиц в плоскости симметрии электрического поля, магнитного поля или совмещенных электрических и магнитных полей с общей плоскостью симметрии [1-10]. Однако, после того как обеспечено желаемое движение пучка заряженных частиц в плоскости симметрии, возникает важный вопрос: что происходит с траекториями заряженных частиц при их незначительном отклонении от плоскости симметрии? Насколько поведение реального пучка с конечной толщиной соответствует идеализированному поведению бесконечно тонкого пучка, которое было рассчитано теоретически?
Рассмотрим типовую схему использования метода обратной задачи по созданию электронно-или ионно-оптической системы, решающей ту или иную практическую проблему. Имеется априорно заданное желаемое поведение пучка заряженных частиц, требуется сконструировать электрические и магнитные поля (а также порождающие их электроды, магниты и катушки с током), которые обеспечат именно такое движение заряженных частиц. При использовании методов [1-10] безупречное поведение бесконечно тонкого пучка гарантируется самой процедурой синтеза распреде-
ления электрического и/или магнитного полей для плоскости симметрии г = 0 . По известным данным о синтезированных в плоскости симметрии электрическом и магнитном полях и с учетом требования симметричности поля можно однозначным образом восстановить поле во всем пространстве или по крайней мере в некоторой окрестности плоскости симметрии. Как результат на выходе можно получить конфигурацию электродов, полюсных наконечников и катушек с током, нужных для создания искомого трехмерного поля.
Однако устойчивость пучка к малым отклонениям от плоскости симметрии будет реализовы-ваться далеко не для всех синтезированных симметричных электромагнитных полей. Задачей данной работы является получение достаточного критерия устойчивости и компактности электронных и ионных пучков в трехмерных электрических и магнитных полях с плоскостью симметрии. Этот критерий удается сформулировать в таком виде, что для его проверки достаточно знания о распределении электрического и магнитного полей в плоскости симметрии (а также знания базовых траекторий в плоскости симметрии); восстанавливать электрическое и магнитное поля во всем трехмерном пространстве для проверки устойчивости при выходе за пределы плоскости симметрии не требуется. В качестве примера полезность рассматриваемого критерия будет продемонстрирована на частном случае анализа работы двух
отклоняющих систем в виде электростатических зеркал, которые характеризуются электрическими потенциалами, однородными по Эйлеру с нулевым показателем однородности (теория таких отклоняющих систем подробно излагается в [3, 10-12]).
УСТОЙЧИВОСТЬ и компактность ПЛОСКОГО ПУЧКА В ЭЛЕКТРОСТАТИЧЕСКИХ ПОЛЯХ С ПЛОСКОСТЬЮ СИММЕТРИИ
Ограничимся для начала частным случаем движения заряженных частиц в чисто электростатическом поле с плоскостью симметрии. Для того чтобы гарантированно обеспечивать устойчивость (локализацию в окрестности плоскости симметрии) пучка заряженных частиц при малых отклонениях от плоскости симметрии г = 0, полезна следующая теорема.
Теорема. Имеется трехмерный электрический потенциал и (х, у, г) с плоскостью симметрии г = 0 : и ( х, у, - г ) = и (х, у, + г) и рассматривается
движение заряженных частиц в окрестности плоскости г = 0. Чтобы при малых отклонениях начальных условий заряженной частицы от плоскости OXY движение оставалось глобально локализованным в окрестности плоскости г = 0, достаточно выполнения на множестве точек плоскости ОХУ, заполненного лежащими в плоскости симметрии траекториями, условия д2и/дг2 > е2 > 0 в случае положительно заряженных частиц и условия д2и/дг2 <-е2 < 0 в случае отрицательно заряженных частиц.
Из этой теоремы, доказательство которой будет приведено далее, имеются полезные следствия. Как правило, эффектами объемного заряда можно пренебречь, поэтому электрический потенциал и (х, у, г) должен удовлетворять уравнению Лапласа ихх + иуу + игг = 0 (где нижние индексы для
и обозначают частные производные по соответствующим переменным). Поэтому если известна функция и0 ( х, у ) = и (х, у,0 ) — распределение
электрического потенциала для плоскости симметрии ОХУ, но не известен сам электрический потенциал и (х, у, г) , то проверка сводится к условию Ух, у: и°. + и^у < -е2 < 0 в случае положительно заряженных частиц и к условию Ух, у : и1 + и°уу > +е2 > 0 в случае отрицательно заряженных частиц. Если же объемным зарядом р( х, у, г) пренебречь нельзя и электрический потенциал и (х, у, г) удовлетворяет уравнению Пу-
ассона ихх + иуу + игг = -4жр (предполагается отсутствие диэлектрической проницаемости среды), проверка сводится к условию Ух, у:
то
и 0
+ иу + 4жр0 <-е2 <0
заряженных частиц
и 0
+ и° + 4жр0 >+е2 > 0
в случае положительно
и к условию Ух, у:
в случае отрицательно
заряженных частиц (где р0(х,у) = р(х,у,0) —
распределение плотности объемного заряда в плоскости симметрии).
Перед тем как приступить к доказательству теоремы, рассмотрим физический смысл сформулированного условия устойчивости пучка. Для простоты ограничимся положительно заряженными частицами. В силу симметрии электростатического потенциала выполнено условие иг| = 0.
Поэтому, во-первых, частица, стартующая в плоскости симметрии г = 0, продолжает движение в плоскости симметрии, а во-вторых, разложение потенциала при малых отклонениях от плоскости ОХУ имеет вид:
и ( х, у, г )« и ( х, у,0 ) + 2 г 2Ц^ ( х, у,0) +
+ ¿г4игггг (х,у,0) + -,
где разложение содержит только четные степени г в силу симметричности потенциала относительно плоскости ОХУ . Это означает, что в случае положительного значения коэффициента и гг(х, у,0) при малых отклонениях г независимо от знака возникает электрическая сила еЕ2 = -е (ди/дг) « -еги г (х, у, 0) , возвращающая заряженную частицу в сторону к плоскости ОХУ . В случае же отрицательного значения коэффициента игг(х,у,0) при малых отклонениях г возникает электрическая сила, еще более уводящая частицу от плоскости ОХУ .
Само по себе наличие возвращающей силы в любой точке траектории, вообще говоря, не означает устойчивости движения, т. к. вполне возможны случаи, когда
а) возвращающей силы недостаточно, чтобы вернуть частицу к плоскости ОХУ ;
б) возвращающая сила слишком сильна, чтобы локализовать частицу в окрестности плоскости ОХУ , так что частица пересекает плоскость и уходит от нее вдаль с другой стороны плоскости симметрии;
в) происходит резонансная раскачка колебаний частицы относительно плоскости симметрии, когда за счет возвращающей силы, разнонаправленной с двух сторон от плоскости симметрии, части-
ца многократно пересекает плоскость симметрии со все более возрастающей амплитудой колебаний.
Точно так же само по себе отрицательное значение и(х,у,0) не означает запрета на локали-зованность движения частицы в окрестности плоскости ОХУ, поскольку на некотором расстоянии уводящая от плоскости сила может смениться возвращающей силой (соответствующий пример будет рассмотрен ниже). Наконец, знак величины игг (х, у,0) может меняться вдоль траектории частицы, так что зоны локальной устойчивости сменяются зонами локальной неустойчивости и ответ на вопрос, каким будет совокупный эффект, требует отдельного скрупулезного исследования.
Для строгого доказательства нашей теоремы воспользуемся теоремой Чаплыгина о дифференциальных неравенствах [13]:
Теорема Чаплыгина. Пусть имеется обыкновенное дифференциальное уравнение п -го порядка f(п)- F (*, f, f' ,..., f (п-1)) = 0 для функции / ^ )
с начальными условиями / (t0 ) = /0, /' (^ ) = /0', •••, /(п-1)(^) = /0(п1). Пусть также имеется функция g ^) , удовлетворяющая начальным условиям g{Ч) = /о, g'(^) = /0, •, ¿п-1)(^) = /о(п-1) и неравенству g(п) -F(t,g,g',...,g|(п-1))> 0 на интервале t e[t0, Т] при t > t0 (где Т > t0). Пусть, кроме того, функции / ^) и g (t) дифференцируемы и непрерывны вплоть до п -го порядка, а функция F (t, /, /' ,..., /(п-1)) и ее частные производные по
аргументам /, /' ,..., / ^п-1) ограничены для интересующих нас диапазонов значений /,/',...,/(п-1) и независимой переменной t. Тогда существует такое значение Т0 > t0, Т0 < Т, что для всех t Т0 ], t > t0, выполнено неравенство g ^) > / (t) . Аналогично, если есть функция h (t) , удовлетворяющая начальным условиям h ) = /0,
А '(^о ) = /0, •, h(n-1)(tо ) = /о(п-1) h(n) -F(^А,А',...,А^п-1))<0, то существует такое значение Т0 > t0, Т0 < Т , что для всех t е[^, Т0 ], t > t0, выполнено неравенство А (t) < / ^) .
Теорема Чаплыгина может использоваться для оценки (мажорирования) решений дифференци-
ального уравнения. А именно, пусть /^) — решение уравнения /^ = F (t, /, /',.., /(п-1)), g (t) — решение уравнения gl-n) = G g, g',..., gl(n-1)) = = 0, а А ^) — решение уравнения А( п) = = Н (t, А, А' ,..., А^п-1)) с одними и теми же начальными условиями. Если для интересующей нас области значений функций и независимого переменного выполнено неравенство Нt0,tn-1 )<
<F(^t0,t1,.,tn-1 )<G(t,t0,t1,.,tn-1), то из теоремы Чаплыгина следует, что существует такое Т0 > t0 (возможно, равное бесконечности), что А (V)< / (V)< g (t) на рассматриваемом интервале
t е[^, Т0 ], t > Искусство состоит в подборе таких функций G и Н , чтобы, во-первых, соответствующие уравнения имели простые решения, во-вторых, чтобы функции G и Н достаточно тесно "зажимали" между собой функцию F, в-третьих, чтобы критическое значение времени Т0, за пределами которого теорема Чаплыгина перестает быть верной, было расположено вне интересующего нас интервала времени.
Примечание. Хотя теорема Чаплыгина и постулирует существование предельного значения Т0, выше которого рассматриваемые неравенства перестают быть справедливыми, но ее доказательство включает в себя вполне конструктивные средства для вычисления предельного значения Т0 по заданному дифференциальному уравнению. Однако для нашей задачи это ограничение несущественно — как правило, предельное значение Т0 много больше времени Т движения заряженной частицы вдоль рассматриваемой траектории, поэтому в дальнейшем мы не будем делать разницы между Т0 и Т . Подробности см. в [13].
Приступим к доказательству нашей теоремы. Пусть х ^), у (t), г ^) — решение системы уравнений движения для положительно заряженной частицы с массой т и зарядом е
и неравенству
тх
" = -еих (х,у,г), ту" = -е^ (х,у,г), тг " = -еиг ( х, у, г)
с начальными условиями, лежащими в окрестности плоскости ОХУ . Зафиксируем функции х ^), у ^) и рассмотрим отдельно уравнение для
функции г (t), которое при заданных х ^), у ^) становится самодостаточным:
mz" = -eUz (х(t),y (t),z) .
(1) выполняются следующие оценки: Uz (х,y, z)
>
По условию теоремы для любых х, у , лежащих в области движения частицы, выполняется требование игг (х,у,0)> е2 > 0. Выберем диапазон значений переменной г е[0, к] достаточно малым, чтобы была справедлива оценка
У Г е^, Т ], У г е[0, к]:
и2 (х (X), у (X), г)/г >е?2 > 0, (2)
либо в более сильной форме
Ух, у е ОХУ, У г е [0, к]:
и2 (х,у,г)/г >е?2 > 0. (2а)
Это возможно, поскольку при г е [0, к] и к ^ 0 величина
иг (х, у, г) = г • игг (х, у,0) + +1 г3 • игг (х,у,0) + - = г ив (х,у,0) + О(к3)
оказывается сколь угодно близкой к величине гиг(х,у,0) . Поэтому величина иг(х,у,г)/г
будет регулярной функцией в окрестности точки = 0 , а поскольку в соответствии с заданным условием величина игг(х,у,0)> е2 > 0 положитель-
iin1(Uz ( х, У, z)/z) пРи
min
e[0,h], ie[i0,T]
min
min U (х, y,z )/z
ze[0,h], х, yeOXYX
min (U (х, y,z)/z)= liL > 0.
ze[0,h], х, yeOXYV zV " > П
(3а)
> z • s,.
при г > 0, иг (х, у, г)< г •ек при г < 0, иг (х, у, г)< г • при г > 0, иг (х, у, г)> г • при г < 0.
Пусть г0, г'0 — начальные условия для уравнения (1). Рассмотрим в качестве функций, используемых в теореме Чаплыгина, функции и (X) и ^ (X), имеющие те же начальные условия г0, г'0 и удовлетворяющие уравнениям
ти" + eF (и) и = 0, mw" + eG w = 0, (4)
где положительные кусочно-постоянные функции F (г ) и G ( г ) определены как
sl при z > 0, G (z) = _Jl1 пРи z > 0, (5)
F ( z )=, 2
I Lh при z < 0;
на, то величина тт \иг(х,у,г)/г) при достаточ-
ге[0,к]
но малых к тоже обязана быть положительной. Очевидно, что существуют разнообразные сочетания к иг?, обеспечивающие выполнение оценки (2) или (2а), однако обсуждение оптимального подбора к и е отложим до конца рассуждений.
В результате можно определить верхние и нижние грани величины иг (х, у, г)/г, каждая из которых будет положительной при выборе значения к достаточно малым:
(и(х(X),у(X),г)/г) = еЛ2 > 0,
(3)
[0,Ч, Ц,,^ ( х ( X ), у ^ ), г Vг ) = ^ > 0,
(иг(х,у,г)/г) = е2 >0,
В таком случае У г е [-к, +к] (в силу симметричности потенциала и соответственно антисимметричности компоненты напряженности поля иг )
sh при z < 0.
Поскольку Vz e [-h, +h] Uz (х,y, z)> F(z)• z, Uz (х, y, z)< G (z)• z, то при сравнении уравнений (4) с уравнением (1) в соответствии с теоремой Чаплыгина получаем оценку w (t) < z (t) < u (t) ,
т. е. |z(t)| < maxЦu (t)|,|w (t)|j, если, конечно, u(t)e[-h, +h] и w(t)e[-h, +h]. Однако функции u (t) и w (t), являющиеся решениями уравнений (4) с кусочно-постоянными функциями (5), — это кусочно-синусоидальные функции, глобально ограниченные по амплитуде пропорционально начальным условиям (ниже в статье приведены их графические образы). Поэтому всегда можно так масштабировать начальные условия z0, z" (подобрать их достаточно маленькими), чтобы требования u(t)e[-h, +h] и w(t)e[-h, +h] были выполнены. На этом доказательство можно считать законченным: показано, что при выполнении условий теоремы и при достаточно малых возмущениях начальных условий пучка в перпендикулярном направлении траектории будут оставаться глобально локализованными в окрестности плоскости симметрии (в пределах z(t)e[-h, +h]) в течение
всего времени движения заряженных частиц. Сам же параметр h можно выбрать сколь угодно близким к нулю, так что при выборе возмущений начальных условий пучка в перпендикулярном направлении достаточно малыми пучок можно сделать сколь угодно тонким.
Примечание. Для применения теоремы Чаплыгина также требуется, чтобы функции u (t) и w (t) были дважды непрерывно-дифференцируемы, что
эквивалентно требованию, чтобы правые части как исходного дифференциального уравнения, так и дифференциальных уравнений, используемых для сравнения решений, были непрерывно-дифференцируемы. Для решений уравнений (4), представляющих собой кусочно-синусоидальные функции, и для кусочно-разрывных функций (5), используемых для правых частей уравнений (4), эти условия не выполнены. Однако этот момент не является принципиальным: функции (5) можно локально изменить в окрестности y = 0 таким образом, чтобы а) они стали непрерывными, б) сохранились неравенства Vz e[-h, +h]:
Uz (x, y, z)> F (z)• z, Uz (x, y, z)< G (z)• z, в) по крайней мере, на фиксированном интервале времени было обеспечено сколь угодно малое отклонение решений модифицированных уравнений (4) от кусочно-синусоидальных функций, соответствующих точным уравнениям (4).
Рассмотрим теперь более подробно поведение решений уравнений (4), которые ограничивают в вертикальном направлении траектории возмущенного пучка. Пусть решение u (t) рассматривается на интервале времени, где оно удовлетворяет условию u (t) > 0. Для этого начальные условия z0, z'0 должны удовлетворять определенным требованиям: либо z0 > 0, а z' любое, либо z0 = 0, а z'0 > 0 . Если z0 > 0, то можно оттрасировать частицу назад по времени и найти момент времени, в котором будет выполнено условие z0 = 0 (в силу
колебательного характера функции u (t) такой
момент времени обязательно существует). Поэтому без ограничения общности можно считать, что мажорирующее решение u (t) стартует с начальными условиями t = t0, z0 = 0 и z'0 = z'0 > 0 . В таком случае u (t) = (i'/ю+) sin (©+(t -10)) , где
a>+ = s^e/m , вплоть до пересечения оси u (t) = 0 и изменения знака в момент времени tj = t0 + njm+ . При этом в момент t = tj изменения
знака решения скорость u'(t) становится равной -z'0, т. е. просто меняет знак по сравнению с начальной скоростью. Отрицательная полуволна решения конструируется аналогичным образом
с заменой частоты а>+ на частоту со_ = ¿ihyJe/m , а результирующая формула u (t ) = _( z' / ю_)х х sin (ю_ (t - tj)) будет справедлива вплоть до следующего изменения знака в момент t2 = tj + л/ю_ .
После этого снова последует положительная полуволна u (t) = (z'0lю+)sin(©+(t -12)) и так далее. Решение является периодическим с периодом T = л/®++ n¡a_ . Максимум положительной полуволны равен, как легко видеть, величине +(z'0lю+) , а минимум отрицательной полуволны равен величине -(z'0¡ю-) . В случае отрицательной начальной скорости отрицательная и положительная полуволны меняются местами. По той же схеме рассматривается поведение решения w (t).
Вывод: решения z (t) для уравнений (1) ограничены полосой ±| z¿| ja+ в окрестности плоскости OXY , заполняемой кусочно-синусоидальными решениями уравнений (4). Тем самым возмущенное движение оказывается ограниченным узкой окрестностью плоскости OXY , причем, как легко понять, ширина полосы стремится к нулю, когда возмущение начальных условий траектории z0, z'0 стремится к нулю. Подобные пучки, которые глобально локализованы в окрестности плоскости симметрии, а их толщина стремится к нулю при стремлении к нулю перпендикулярного возмущения начальных условий, будем называть плоскими.
Чтобы завершить анализ, следует установить ограничения на возмущения начальных условий z0, z', при которых мажорирующие решения u (t)
и w (t ) оказываются ограниченны полосой z e[-h, +h], в пределах которой только и имеют силу сделанные оценки, а также соотнести начальные условия z0, z' с шириной полосы, в пределах которой лежит возмущенная траектория. Без ограничения общности можно считать, что либо z0 > 0, а z' любое, либо z0 = 0, а z' > 0 (в силу симметричности электрического потенциала при изменении знака начальных условий z0, z' траектория симметрично отражается относительно плоскости OXY ). При таких условиях начальная полуволна мажорирующей функции u (t) , описываемой уравнениями (4), (5), является ее максимальной полуволной и совпадает с синусоидой, порожденной дифференциальным уравнением mu" + esh (u ) u = 0 и записанной в виде аналитического решения как
u (t) = z0cos (®+ (t -10)) + (z'/ ®+) sin (®+ (t -10 )) = = >/z02 +(z'/®+)2 cos(®+ (t-10)-arctg((y'J)/У0)),
где = .
Максимум +tJz2 + (z'/®+)2 этой синусоиды
и есть глобальный максимум верхней мажорирующей функции и (X) . Соответственно минимум
г2 + (г",/о)+ )2 этой синусоиды — это глобальный минимум функции w (X), ограничивающей решение г (X) снизу. В итоге неравенство
г02 +(г"/®+)2 < г(X)<+>/г02 +(г0/ю+)2 задает
тавляли заряженную частицу в пределах полосы
2 , (_< /„ \2 ^ 7.2
-к < г (X)<+к ,
записывается как г,
+ (г,/®+)2 <к2.
В более наглядной форме "энергетического неравенства" ограничение на допустимые начальные условия выглядит как
т(г")2/2 + е(е2г02/2)<е(е2к2/2) .
полосу, в пределах которой оказывается траектория с возмущенными начальными условиями г0, г',. Требование, чтобы начальные условия ос-
Графики решений и (X) и w (X) показаны на рис. 1.
Рис. 1. Графики функций )
(а) и w(x) (б) при различных
начальных условиях у 0, у 0
при фиксированной энергии возмущения
Е = т(г0)2/2 + е(екгЦ2) ,
описывающей отклонение начальных условий от плоскости ОХУ
Из "энергетического неравенства", в частности, становится понятно, как именно следует выбирать компромисс между величиной h и параметром sh для неравенства (2): надо максимизировать площадь эмиттанса допустимого возмущения z0, z0,
описываемого неравенством m(z0)2/2 + e(s2z2/2)< < e(s2h2/2) • Задача поиска оптимальных параметров h и sh не вполне тривиальна, поскольку, приближая sh как можно ближе к нулю, мы увеличиваем h и тем самым увеличиваем диаметр фазового эллипса, равный его большей полуоси. Но, с другой стороны, приближение sh как можно ближе к нулю уменьшает до нуля величину , делая эллипс узким и сводя его площадь к нулю. Окончательная же оценка зоны устойчивости плоского пучка заряженных частиц получается объединением всех возможных эллипсов
m (z0 )2 /2 + e Sz2 /2) < e (s^h2/2) при всех возможных оценочных параметрах h , для которых величина s^ = min (Uz (x, y, z)/z) в неравенствах
(3), (3а) остается положительной.
Примечание. Аналогичным образом можно было бы построить приближенные мажорирующие (снизу) решения для случая Uzz (x, y,0)<-s2 < 0. Однако результат оценки не
будет работоспособным, поскольку соответствующие решения — экспоненциально растущие и быстро выходят за границы полосы z е +h],
в пределах которой сделанные оценки только и остаются справедливыми. При более детальном рассмотрении оказывается, что хотя условие Uzz (х,у,0)> 82 > 0 является достаточным для устойчивости плоского пучка, но обратное ему условие Uzz (х, у,0)<-е2 < 0 отнюдь не является достаточным для неустойчивости плоского пучка. В качестве примера на рис. 2 показаны осцилляции заряженной частицы в потенциальной яме
с профилем и (z ) = и о (-Az2 + Bz4), которые являются глобально ограниченными по амплитуде, несмотря на условие и А < 0. Следует, однако,
\г-0
заметить, что локальная ограниченность плоского пучка для этого случая несколько отличается от рассмотренного ранее случая. А именно, асимптотическая толщина пучка, отнюдь, не стремится к нулю при стремлении к нулю начального размера пучка, а стремится к вполне конечной величине, так что тонкий пучок на входе системы быстро превращается в тонкое кольцо конечного радиуса, а не в тонкую ленту. Аналогом этого эффекта является фиксированная ограниченность амплитуды при устойчивости в виде предельного цикла как противовес сколь угодно малой амплитуде при устойчивости в виде фокуса [14-20].
л
0.6 и\у)
0.4 Ео /
\ 0.2
-1.о\ -0.5 \ -0.2 0.5 11.0
КО
7
Рис. 2. Осцилляции заряженной частицы в потенциальной яме и(у) = ио (-Ау2 + БуА).
а — профиль потенциальной ямы и уровень полной энергии, соответствующий начальным условиям; б — нелинейные строго периодические осцилляции у(), глобально ограниченные по амплитуде
б
а
ЭНЕРГЕТИЧЕСКОЕ ОБОСНОВАНИЕ УСТОЙЧИВОСТИ И КОМПАКТНОСТИ ПЛОСКОГО ПУЧКА
Если пренебречь математической строгостью, результаты предыдущего раздела могут быть обоснованы проще.
Пусть отклонения г (X) от плоскости ОХУ имеют величину порядка Л, где Л — малый параметр. После подстановки г (Л, X ) = Лг1 (Л, X) (где г1 (Л,X) — функция порядка О(1) при Л^ 0) в силу симметрии электрического потенциала и, как следствие, в силу условия иг| = 0, закон сохранения энергии можно записать как
т (х"2 + у "2 + г'2)
Е = —^-'- + еи ( х, у, г )*
т ( х "2 + у "2 + Л2 г/2)
2
+ еи ( х, у,0) +
1
+ 2 Л2 &иа ( х, у,0 ) + • а уравнения движения как
(6)
тх
" = -еих (х,у,г)',
*-еих(х,у,0)-1Л2г^ • е^(х,у,0)--, ту" = -еиу (х,у,г)*
* -еиу (х,у, 0) -1 Л2г^ • еи^ (х,у, 0)-
(7)
тг
" = т Лг[ = -еиг (х, у, г) *
* -Лг • еи^ (х,у,0) - 6Л3г3 • еиш (х,у,0) -
Сохраняя в (7) члены не выше Л , а в (6) члены не выше Л2 (т. е. первые ненулевые поправки), получаем соотношения
тх
" = -еих (х,у,0), ту" = -еиу (х,у,0),
р" = -р^ еи^ (х, у,0 ),
(8)
Е =
т ( х "2 + у "2)
+ еи ( х, у,0)
+
+ Л2
тг" 1 2 ТТ / пЧ
+ 2 г'еигг (^ y,0)
+
О (Л4).
Однако из первых двух уравнений (8) следует, т ( х"2 + у "2)
—--- + еи ( х, у,0) — первый инте-
что Е0 =
грал усеченных уравнений движения, поэтому Л(Е - Е0 ) = АЖ =
тг1 2
1 г^ ( х, у,0) + О (Л2)
- + - г 2
I = const.
(9)
Известно, что Ух,у: иг (х,у,0)> е2 > 0, поэтому из условия (9) следует, что при ненулевых значениях г1, г" и достаточно малых Л величина АЕ всегда положительна. В случае, когда известна априорная оценка 0 < АЖ < АЖ0, допустимые значения г1 и г^, удовлетворяющие условию (9), образуют компактный эллипс, зависящий, впрочем, от текущих значений х(X),у (X) . По результатам решения системы (8) г1 (X) и г^ (X) лежат в пределах 4-мерной полосы, составленной из фазовых эллипсов (9) для переменных г1, г" с центрами г1 = 0, г" = 0, которые размещаются в каждой точке х,у плоскости ОХУ . "Толщина" этой
полосы определяется величиной АЕ = Л2 • А Ж, которая представляет собой малую энергетическую поправку порядка Л2, описывающую возмущение начальных условий движения в направлении, перпендикулярном плоскости ОХУ .
Возвращаясь к переменным г(X), г '(X) , из равенства (9) следует, что
|г(X)|<^/2АЕ/(еиг (х(X),у (X),0)) <д/2АЕ/ее2 ,
т. е. представляет собой ограниченную величину порядка Л при условии, что АЕ представляет собой величину порядка ЛЛ. Этот результат как раз и означает, что пучок заряженных частиц будет оставаться локализованным в окрестности плоскости ОХУ, если возмущение начальных условий, выводящее пучок за пределы плоскости ОХУ, достаточно мало.
Уточним энергетическую оценку размера пучка с возмущенными начальными условиями. Пусть х(X),у (X), г (X) = 0 — решение уравнений движения с начальными условиями х0, х'0, у0, у'0, г0 = = 0, г'0 = 0, лежащее целиком в плоскости симметрии г = 0 . Добавим возмущение начальных условий в поперечном направлении: г0 = Лг0, г'0 = Лг<0, где Л — малый параметр. Представим решение трехмерной системы уравне-
2
ний движения в виде ряда по малому параметру Л как
х (Л, t) = х (t) + Лх1 (t) + Л2 х2 (t) + - • •, у (Л, t) = у ) + Лу1 (t) + Л2 у 2 (t) + -, 2(Л,t) = Лгх (t) + Л2г2 ^) + —.
В силу симметричности электрического потенциала и ( х, у, 2) в разложении для х (Л, t) и у (Л, t) будут присутствовать только четные степени Л, а в разложении для 2(Л,t) будут присутствовать только нечетные степени Л. Ограничиваясь членами, не выше чем квадратичными, для закона сохранения энергии получаем соотношения:
Е = Е0 + Л Е2 + Л Е4 +— = сопй,
т ( х'2 + у '2)
Ео -+ еи ( х, у,0 ),
Е2 = тх'х'2 + ту'у'2 + еих (х,у,0) • х2 + т2' 2 1
+ еиу (х,у,0) • у2 + + -еиг1 (х,у,0)^ 22
т ( х22 + у22 ) 1 Еа -+ - еих (х, у,0) • х2 +
2
2
1
+и (х y,0)• х2 у 2 + 2 еиуу (х y,0)• у2,
(Выражение для Е4 является не вполне точным, поскольку должно содержать также и коэффициенты х4 (t), у4 (t), 23 (t) разложения решения в ряд
по малому параметру Л , однако ранее мы договорились этими поправками пренебречь).
Выражение для Е2 описывает компактный эллипс для 21, , если значения х2, у2 меняются в ограниченном диапазоне или если значения х2, у2 пренебрежимо малы на фоне значений 21 (второе условие, конечно же, является более сильным, чем первое). Выражение для Е4 описывает компактный четырехмерный эллипс для х2, у2, х2, у2 и тем самым обеспечивает глобальную ограниченность для х2, у2, х2, у2, если квадратичная форма и^х^ + 2ихух2 у2 + иявляется
положительно определенной. Это выполнено тогда и только тогда, когда их > 0, иу > 0,
и и - и2 > 0.
хх уу ху
Однако эти условия находятся в очевидном противоречии с условием и^ + и <- е2 < 0
(в частности, поскольку и^ <-е2 - иуу, то
и и - и2 < -е2и - е 2и2 - и2 < 0). Соответ-
хх уу ху уу уу ху '
ственно использованное ранее для доказательства глобальной ограниченности 21, предположение, что значения х2, у2 меняются в глобально ограниченном диапазоне, оказывается под вопросом. Это не означает, что траектория будет неустойчивой, это лишь означает, что мы не можем вывести устойчивость траектории, исходя только из энергетических соображений.
Кроме того, это наводит нас на мысль, что одновременная глобальная ограниченность возмущения траектории как внутри плоскости, так и перпендикулярно плоскости невозможна. Этот вывод будет почти очевидным, если мы вспомним, что электрический потенциал удовлетворяет уравнению Лапласа и тем самым не может иметь минимумов и максимумов, но лишь седловые точки. Опять-таки, отсутствие глобальной ограниченности траектории сразу во всех направлениях по отношению к отклонению пучка от плоскости симметрии не означает, что невозможно добиться фокусировки пучка по двум направлениям сразу: ничто не мешает возмущенным траекториям сперва расходиться, а потом снова сходиться в точку. Принципиальным вопросом, однако, является то, что, по-видимому, возмущенный пучок не может быть глобально узким вдоль базовой траектории по двум направлениям одновременно.
Очевидный недостаток всех этих выкладок в том, что в качестве исходного предположения выбирается именно тот тезис ("малость" 2 ^) при
всех t), который требуется доказать. Достоинством же подхода, изложенного в предыдущем разделе, является не только математическая строгость выкладок, но и возможность конструктивно оценить полосу устойчивости движения в окрестности плоскости OXY , а также получить в явном виде зависимость толщины плоского пучка от энергии возмущения при отклонении частицы от плоскости ОХТ .
ОЦЕНКА ЛОКАЛЬНОЙ ТОЛЩИНЫ ПУЧКА
При разложении решения в ряд по малому параметру Л и представлении его в виде
х (t) = Лх1 ^ ) + Л2 х2 ^) + —, у (t ) = Лу1 ^) + Л2 у 2 (t) + -, 2 ^) = Л21 (t) + Л222 ^) + —
для функции 21 (t) получается самодостаточное дифференциальное уравнение тг" =
= -еигг (х(t),у(у),0)г,. Поэтому изменение во
ч гаг,'2 1 т 2 времени энергии Е ^ ) =—-—+ — еи^г, для возмущения, перпендикулярного плоскости симмет-
1 2
рии, описывается уравнением Е = — ег, х
х(ихггх ' (t) + иуггу'(t)) . Считая, что г, ^) меняется в ограниченных пределах, так что у0 < т <t: - Zmюí ^)< г, (т)< 2тах (t), можно, интегрируя функцию времени f' ^ ) = ихггх' (t) + +иуггу' ^) вдоль конкретной траектории, определить максимально достижимое значение энергии Е ^) в каждый момент времени и тем самым максимальную толщину возмущенного пучка в каждой точке базовой траектории.
Задача упрощается, если учесть, что для получения максимальной оценки значения энергии
Е ^) интегрировать функцию ихггх' ^) + иуггу' (t)
следует лишь по тем участкам траектории, где эта функция положительна, и что функция
У ) = ихггх ' (t) + иуг?у' ^ ) — ЭТ° произвОДная по времени от функции У (t) = иш (х (t), у (у)) . Поэтому интересующая нас вариация 5У ^) — это сумма приращений величины игг (х(t), у ^)) вдоль тех участков траектории, где функция У (t) = иш (х (t), у (у)) возрастает.
Тем самым в каждой точке траектории мож-
тг,' 2 1 ТТ 2
1 + "еи^ <
2 2
но построить эллипс
< Е20 +—е •5/ ^) 2ттах (t) для допустимых значений
г,, , описывающий возможную толщину сопровождающего пучка в конкретной точке конкрет-
тг '2 1
ной траектории. Здесь Е° = ^ + —еи^г^ — это
начальное значение величины
Е2 (t),
описываю-
самосогласованной оценки 2тах (t) представляется
достаточно сложной задачей в рамках используемого здесь "энергетического" формализма, особенно если функция иг (х(t), у ^)) не является монотонной вдоль траектории и, как следствие, локальная амплитуда колебаний функции г, ^) также не является монотонной функцией.
Если считать, что иг (х ^), у ^)) мало меняется на протяжении одной полуволны решения уравнения тг' = -е иг (х (у), у ^) ,0) г,, то легко
понять, что при локальном увеличении иг амплитуда полуволны падает, а при локальном уменьшении игг амплитуда полуволны растет. В частном случае, если игг (х ^), у (t)) монотонно падает
вдоль траектории, то 5У ^) = 0 и значение 2тах ^) становится несущественным. Тем самым локальная оценка толщины пучка в заданной точке траектории вычисляется как гтах = ^2Е°/еигг , где
Е\ — это начальное значение энергии Е2 ^) для
отклонения траектории от плоскости симметрии. Как и следовало ожидать, эта функция монотонно растет вдоль траектории в силу монотонного убывания иг (хотя с учетом того, что Е2 (t) — монотонно убывающая функция, данная оценка слегка завышена).
Ожидается, что для монотонно растущей функции У ^) = иг (х(у), у ^)) амплитуда колебаний г (t) будет монотонно падать по ходу траектории
(более-менее строгое обоснование этого факта см. далее), поэтому максимальное значение амплитуды колебаний (t) = Е20/\ 1 еи"г | достигается
щей энергию отклонения траектории от плоскости симметрии; 2тах (t) — оценка максимального по модулю значения г, (t) на интервале времени от точки старта до текущей точки; а для сокращения записи введено новое обозначение и°г = иг (х0,у0,0). В свою очередь требуется оценка для величины 2тах ^) максимального отклонения траектории от плоскости симметрии между начальной и текущей точкой. Получение же
2
в начальной точке, вариация энергии Е2 ^) равна 5/ (t) = игг (х (t), у (у))- и]2, а оценка для толщины пучка в заданной точке траектории сводится к неравенству гтах < ^2Е° /еи° , соответствующему оценке амплитуды колебаний функции г, ^)
в начальной точке. Данная оценка является основательно завышенной и не соответствует нашим ожиданиям, что толщина пучка должна падать вдоль длины траектории при монотонно растущей функции иг. Тем не менее эта весьма грубая оценка правильно описывает глобальную ограниченность пучка и стремление толщины пучка к нулю при стремлении к нулю энергетического возмущения, связанного с отклонением начальных условий от плоскости симметрии.
Более содержательную информацию о поведении толщины пучка, т. е. об амплитуде колебаний функции 21 ^) вдоль рассматриваемой траектории, можно получить, аккуратно проанализировав параксиальные уравнения траекторий.*'1 Один из способов это сделать — использовать подстановку Крылова—Боголюбова [21-24].
При применении этой процедуры уравнение
т^" = -еи2 (х^),у (t),0)21 для параксиального решения 21 (t) удобно записать в виде системы
tf (t) = ~ f (t)Sln2 7 + Я> ), A (t ) = " 2¡k f (t) sin (2yt + 2tf);
tf2 (t) = - f (t) tf (t) Sin (2yt + 2p),
7
A (t) = -- f (t)tf(t)cos(27t + 2tfo); 7
Г' = -z, Í72 +
- (72 + ¿f (t)),
(10)
где Л — малый параметр, у2 > 0 — оценочное значение величины (е/т(х(t),у ^),0) вдоль конкретной базовой траектории; а медленно меняющаяся во времени величина Лf ^) описывает
изменение (е/т)и22 (х^),у(t),0) вдоль выбранной базовой траектории по отношению к ее оценочному значению у2 > 0 . После этого решение для (10) ищется в форме 21 ^) = ехрА(t)• smр(t), Г (t) = уехрА(t)• соэ^^), где А(t) и р^) — две новые функции, однозначным образом связанные с прежними функциями 21 (t) и г (t) . После элементарных выкладок для функций А({) и р^) получаются уравнения
A' (t ) = -— f (t) sin2tf(t),
27
tf (t ) = 7 + ¿f (t) sin2 tf( t).
(11)
Остается представить функции А^) и р^) в виде рядов по малому параметру Л таким образом, чтобы нулевое приближение соответствовало точному решению системы ^ = г1, г1 ' = -у221, получающейся при Л=0:
А^ ) = А0 + ЛА1 (t) + Л2 А2 (t) + —, р^) = р0 + уt + Лр ^) + Л2р2 (t) + —.
В результате получается цепочка уравнений
) В частности, это позволяет обосновать предположение, что с локальным увеличением локальная амплитуда колебаний падает, а с локальным уменьшением и — растет.
или с учетом того, что р(t0 ) = 0, A1 (t0 ) = 0, tf2 (t0) = 0, A2 (t0) = 0, ..., цепочка рекуррентных формул
1 1
tf (t) = - j f (t)sin2 (7t + tf)dt,
7 t
1 '
A (t) = -—j f (t) sin (27t + 2tf0)dt;
t0
1 t
tf (t) = - j f (t) tf (t) sin (27t + 2tf0) dt, (12)
7 t0
1t
A (t) = — j f (t) tf (t)cos (27t + 2tf0)dt; 7 i
откуда сперва можно получить решения для поправочных членов A (t) и tf (t), потом для поправочных членов A2 (t) и р2 (t), и т. д.
В случае, когда за период T = n¡7 изменение функции f (t) мало (функция является "медленной" по сравнению с частотой колебаний), при интегрировании выражений (12) функцию f (t)
можно вынести за знак интеграла. Более аккуратные результаты можно получить с помощью формулы интегрирования по частям [25]:
jf (г)F(г)dr = f (t)F (t)-jf '(t)F (t)dr =
t0 *O
= f (t)F (t) - f '(t)F2 (t) + jf "(г)F2 (t)dr =
*0
= f (t)F1 (t)-f"(t)F2 (t) +
+ f "(t) F3 (t )-jf"(r) F3 (r) dr =..., (13)
Z1= r1
где Рк (t)=ГРк-1 (т)dт , Р0 ^) = Р(t) — повторные интегралы от функции Р (), обращающиеся в ноль при t = tQ. Если У() — "медленная" функция по сравнению с рассматриваемым масштабом времени, ее производными начиная с некоторого номера можно пренебречь и тем самым из формул (13) можно корректным способом вывести удобные приближенные выражения для решений системы (12).
Точку отсчета ^ = 0 удобно сместить так, что начальная фаза ср0 = 0. Тогда формулы (12) в сочетании с (13) дают приближенные решения
/ \ t г/ sm2уt
р (О— У (0 1 -
2у'
-—У' ()
4у
2^
1 -
1 - соэ2^
2,2
<+—У (t)
2у 1
(
1 -
2у I
sm2уt 2yt
А1 ^ ) = -^ У (t) sin2 ^ +
2У
(14)
4у
^ У' (t)
1 -
sin2уt 2уt
+ • •• ?
2у
-У()sin2 уt,
£(0
V У V
V
V II
Рис. 3. Осцилляции решений линейного дифференциального уравнения ^) = -^)(1 + Ы) с нормированными начальными условиями £ (0) = 0, £' (0) = 1 при (а) ы = +0.01 и (б) ы = -0.01
Соотношения (14), в частности, подтверждают тот факт, что с адиабатическим увеличением значения и22 (т. е. при У(7)> 0) локальная амплитуда колебаний ехр А (t) падает, а с уменьшением (т. е. при У ^)< 0) — растет. Локальная частота с(t) = (р(t)-Р0)/1, как легко видеть, ведет себя
прямо противоположным образом. Для иллюстрации на рис. 3 приведены графики решений нормированного уравнения (t) = -£(t)(1 + Ы) при положительных и отрицательных значениях параметра ы , которые для этого частного случая могут быть выражены в аналитической форме через специальные функции.
Альтернативный вариант — использование метода ВКБ (Вентцеля—Крамерса—Бриллюэна) [26], который обычно применяется для решения уравнения Шредингера. А именно для уравнения
(10), записанного в форме г' = -г, (у2 + ХУ^)) ,
при условии положительности значения
у2 + ХУ (t) решение ищется в виде г, (t) = ехр^у(t)) . После этой подстановки уравнение приобретает вид -(у '2) + ^ " +
+ (у2 + ХУ(t)) = 0. Предполагая, что у" мало по
сравнению с остальными членами (т. к. при Х « 0 у точного решения имеем у"«0), получаем
в первом приближении у' = у2 + () ,
у = ±j^У+Хf(tjdt. Два разных знака соответствуют двум независимым приближенным решениям, и при переходе к исходному линейному дифференциальному уравнению для г, ^) эти два решения должны быть скомбинированы в одно решение общего вида в виде линейной суперпозиции двух независимых решений. Однако, как легко видеть, это решение слишком грубое и не отражает всех особенностей ожидаемого поведения функции г, () : в частности, при медленном изме-
а
2
+
+
1
нении у2 + Лf (t) вдоль невозмущенной траектории меняется локальная частота осцилляций, но не их амплитуда.
Условие применимости полученного приближения (малость у ") выглядит как
у
л/ )
<< у2 + Л/(t). Это
позволяет
У+/)
найти более точное приближение, а именно, пред. . 1 лг(г)
положив, что № « ±-
уравнение для
2 у/)
функции у( t) дает уточненное решение в виде
(у ')2 - у2 + Л/^)± 1
1
Л/' (t)
№ «±^7/) + \
^у2 + Л/ ^)
Л/'(t) 4 у2 + Л/^) ,
у ( t)« + ^1п (у2 + Л/ ^))
(фактически здесь использован еще один шаг итерационного уточнения решения, когда от предыдущей итерации берется значение у " и с помощью квадратур вычисляется новое значение у(t)). Возвращаясь к функции 21 (t) и комбинируя независимые решения, соответствующие разным знакам перед корнем, получаем приближенное решение
21 (t)',
1
(у2+Л/ (t ))14
Са ехр (+i{д/у2 + Л/(х^)
+
С
+СЬ ехр(-i{^у2 + Л/^)dt)
((0 + у2 + Л/ (t )dt)
ПРИМЕР. ЭЛЕКТРОСТАТИЧЕСКИЕ ЗЕРКАЛА С ПОТЕНЦИАЛАМИ ДОНКИНА
В качестве примера рассмотрим системы с электростатическими полями, синтезированными на основе формулы Донкина [10-12]. А именно рассмотрим электростатические потенциалы и (х, у, 2) , которые представляют собой
функцию, однородную по Эйлеру с показателем однородности, равным нулю:
У Л,Ух, у, 2: и (Лх,Лу,Л2) = и (х,у, 2) .
Подобные потенциалы обладают специфическими полезными электронно-оптическими свойствами [10-12, 27-30].
В силу того, что потенциалы и ( х, у, 2 ) должны удовлетворять уравнению Лапаласа
иж + иу + и2 = 0 , можно показать, что любой
такой потенциал можно выразить с помощью формулы ([31, пример 23 к главе 18]):
и ( х, у, 2) = (
= и0 ¥
у
х + ^/х2 + у2 + 22
х + д/х2 + у2 + 22
(15)
(у2+Л/ ^))
которое адекватно описывает изменение амплитуды и частоты локальных колебаний функции 21 ^)
при медленном изменении величины у2 + Л/ ^) вдоль базовой траектории.
Более детальное рассмотрение данного вопроса выходит за рамки задачи, поставленной в этой статье.
где ¥ (р, q) — надлежащим образом выбранная гармоническая функция, т. е. функция, удовлетворяющая двумерному уравнению Лапласа ¥рр + = 0 . Легко проверить, что при произвольной гармонической функции ¥ (р, q) электрический потенциал (15) будет удовлетворять трехмерному уравнению Лапаласа и при этом будет представлять собой функцию, однородную по Эйлеру с показателем однородности, равным нулю. Для того чтобы электрический потенциал (15) был симметричен относительно плоскости 2 = 0, необходимо и достаточно, чтобы функция ¥ (р, q) была четной по аргументу q . Для того, чтобы электрический потенциал (15) был антисимметричен относительно плоскости у = 0 и чтобы соответственно плоскость у = 0 могла использоваться как входной (заземленный) электрод электростатического зеркала, необходимо и достаточно, чтобы функция ¥ (р, q) была нечетной по аргументу р .
Простейший электрический потенциал вида (15) получается при выборе ¥ (р, q) = р (как легко видеть, эта функция одновременно и четна по аргументу q , и нечетна по аргументу р). В [11, 12] в качестве одного из примеров электронно-оптических систем на основе потенциалов, основанных на формуле Донкина, рассматриваются
2
Рис. 4. Эквипотенциальные линии и эквипотенциальные поверхности для электрического потенциала (16).
а — эквипотенциальные линии в сечении ОУ2 при х = 1; б — эквипотенциальные поверхности (полезадающие электроды) для значений потенциала и = 0.0 и и = 0.9
электростатические зеркала, описываемые потенциалом
, иУ 2 2 при У * 0, , ,
х + у х + У + z (16)
U (х, y, z ) =
Эквипотенциальные линии и эквипотенциальные поверхности для электрического потенциала (16) в сечениях х = 1 и z = 0 показаны на рис. 4. (Трехмерные эквипотенциальные поверхности для потенциалов вида (15) представляют собой конические поверхности с центром в начале координат и проходящие через эквипотенциальные линии в фиксированном сечении х = const). На рис. 5 показан профиль электрического потенциала (16) при фиксированных значениях х = 1, z = 0 по мере углубления во внутреннюю область зеркала (т. е. при возрастании координаты y), из которого следует, что при U0 > 0 на положительно заряженные частицы электрическое поле (16) действует, как зеркало.
На рис. 6 показаны траектории ионов в плоскости симметрии OXY этого электростатического зеркала. Как и следует из общей теории [10-12],
Рис. 5. График распределения электрического потенциала (16) вдоль прямой линии х = 1, 2 = 0
0
при y < 0.
Рис. 6. Три параллельные моноэнергетические пучка, двигающиеся в плоскости ОХУ и отклоняемые зеркалом на основе электрического поля с потенциалом (16) на разные углы в зависимости от энергии
входной моноэнергетический пучок параллельных траекторий преобразуется в выходной пучок параллельных траекторий, где угол поворота и точка выхода пучка из зеркала зависит от энергии пучка, а степень сжатия пучка определяется законом подобия траекторий, характерным для потенциалов Донкина.
Однако для потенциала (16) величина
(ихх + иуу )|.
и0 у
>/х2 + у2 (х + >/х2 + У 2)
является строго положительной в области у > 0 и, как результат, следует ожидать неустойчивость траекторий при отклонении частиц от плоскости симметрии. Этот факт иллюстрируется на рис. 7, где показаны траектории с одинаковыми начальными координатами и углами в плоскости ОХТ, смещенные в начальный момент вдоль оси OZ с сохранением прежнего направления скорости, которые монотонно расходятся.
Этот недостаток можно исправить, если в формуле (15) использовать более сложные гармонические функции ¥ (р, д), обеспечивающие отрицательный знак для значения (и„ + и„,) . Так, выбе-
V хх =0
рем в качестве прародителя гармоническую функцию
¥ (p, д )=1
(
(
агС£
(
\\
1 + д
+ aгctg
1 - д.
р.
Такой выбор соответствует электростатическому зеркалу с электрическим потенциалом
и ( х, у, 2 ) =
(
и
aгctg — -х
2 у
Л
х + л[х2 + у2 + 22
при у > 0,
при у < 0.
(17)
Рис. 7. Траектории заряженных частиц в электрическом поле с потенциалом (16), которые стартуют с тремя разными энергиями и со смещениями по вертикали вдоль оси OZ .
а — трехмерные траектории; б — проекция на плоскость ОХТ ; в — проекция на плоскость OXZ ; г — проекция на плоскость ОТХ
б
2
0
в
г
Легко проверить, что для симметричного по 2 электрического потенциала (17) значение
будет строго отрицательным при и0 > 0 и у > 0, т. е. что достаточный критерий устойчивости плоского пучка для электрического потенциала (17) выполняется.
На рис. 8 для двумерных сечений х = 1 и 2 = 0 показаны эквипотенциальные линии электрического потенциала (17). На рис. 9 показан профиль этого потенциала вдоль линии х = 1, 2 = 0 при возрастании координаты у, из которого следует, что на положительно заряженные частицы при и0 > 0 это электрическое поле будет действовать как зеркало. На рис. 10 показан стартующий в плоскости ОХУ и отражаемый зеркалом параллельный пучок траекторий с тремя разными энергиями, целиком лежащий в плоскости ОХУ , для которого начальные координаты смещены друг от друга вдоль оси ОХ. Из рисунка видно, что, как и должно
Ц» + )|
2и0 у
4х2 + у2 (х + х2 + у2 )
1.0
0.5
0.0
-0.5
0,
0.0 0.5
1.0
2.0 У
Рис. 8. Эквипотенциальные линии и эквипотенциальные поверхности для электрического потенциала, порожденного функцией (17).
а — эквипотенциальные линии в сечении ОУХ при х = 1; б — эквипотенциальные поверхности (полеза-дающие электроды), соответствующие значениям потенциала и = 0 и и = 0.3
Рис. 9. График распределения электрического потенциала, порожденного функцией (17), вдоль прямой линии х = 1, 2 = 0
Рис. 10. Три параллельные моноэнергетические пучка, двигающиеся в плоскости ОХУ и отклоняемые зеркалом на основе электрического поля с потенциалом (17) на разные углы в зависимости от энергии
а
быть в соответствии с теорией, входной моноэнергетический параллельный пучок, лежащий в плоскости симметрии ОХТ , преобразуется в выходной параллельный пучок, лежащий в плоскости симметрии ОХТ , причем угол поворота пучка зависит от его энергии. На рис. 11 показаны траектории с тремя разными энергиями и фиксированными начальными условиями в плоскости ОХТ , у которых начальные координаты смещены по оси ОХ . Видно, что по сравнению с рис. 7, в, где траектории последовательно расходятся в вертикальном направлении, на рис. 11, в, имеется промежуточная фокусировка траекторий в вертикальном направлении (положение кроссовера для траекторий с разными энергиями, естественно, разное).
УСТОЙЧИВОСТЬ И КОМПАКТНОСТЬ
ПЛОСКОГО ПУЧКА В МАГНИТОСТАТИЧЕСКИХ ПОЛЯХ С ПЛОСКОСТЬЮ СИММЕТРИИ
Рассматриваемая теорема может быть расширена для движения заряженных частиц в симметричном магнитном поле, а также в совмещенных симметричных электрических и магнитных полях.
Пусть индукция магнитного поля В(х,у,2) = -УФ(х,у,2) задана с помощью скалярного магнитного потенциала Ф(х,у,2) [1], которого плоскость ОХТ (т. е. 2 = 0) является
2
в
2
г
Рис. 11. Траектории заряженных частиц в электрическом поле с потенциалом (17), которые стартуют с тремя разными энергиями и со смещениями по вертикали вдоль оси ОХ .
а — трехмерные траектории; б — проекция на плоскость ОХТ ; в — проекция на плоскость ОХХ ; г — проекция на плоскость ОТХ
плоскостью антисимметрии: ф( x, y,-z)=-0( x, y, z).
Движение заряженной частицы в электромагнитном поле [32] без учета релятивистских или квантовых эффектов описывается уравнениями движения Ньютона mr = eE(r) + (e/c)[r x B(r)\, или в скалярных обозначениях
=Л
dx dt
dp e e г i
— =—EX ( ^ y, z) +—I qBz ( x, y, z)- rBY ( ^ y, z)\, dt m mcL J
dy
dt
=q,
dq e e
— =—ey (x y, z) +—|_-PBz ( x y, z)+ rBx ( x y, z)\
(18)
dt m mc
dz = dt " r, dr e
— = - Ez ( x, ^ z) +—L PBY (x, y, z)- qBx (x y, z)\,
dt m
mc
удовлетворяющий
условию
dz dt dr
■ = r,
e
— « z — dt m
fô 2U ô2U - + -
ôx ôy
1
+ -c
( ôBZ
P-
ôy
■q-
ôB
z\ >
ôx
//
(19)
поскольку с учетом того, что электрический скалярный потенциал и (х, у, 2) разлагается в окрестности точки 2 = 0 по четным степеням 2 , магнитный скалярный потенциал Ф( х, у, 2) разлагается в окрестности точки 2 = 0 по нечетным степеням 2 и 8Ф(х,у,2)/& = -ВХ (х,у,2) , выполнены приближенные соотношения:
, ди (х, у, 2) 82и ( х, у,0) Е = -- ~ " ■
ôz
■ - z •-
ôz2
2
= z •
ô 2U ( x, y,0) ô 2U ( x, y,0)
ôx
ôy 2
Bx = -
ôФ(x, y, z) ô2Ф(x, y,0 )
ôx
- z-
ôxôz
где х, у, 2 — декартовы координаты частицы; х ^ ), у (t), 2 (t) — траектория частицы; р (t), д (t), г ^ ) — компоненты скорости частицы; т — масса частицы; е — заряд частицы; t — время; Е ( х, у, 2, t ) = ( ЕХ, ЕТ, ЕХ ) — напряженность электрического поля в точке х, у, 2 в момент времени t; В (х, у, 2, t) = (ВХ, ВТ, ВХ ) — индукция магнитного поля в точке х, у, 2 в момент времени t;
ЕХ =-8и/8х ; ЕТ = -8и/8у ; ЕХ =-8и/8г ; ВХ = -8Ф/8х ; ВТ = -8Ф/8у ; ВХ = -8Ф/8г ; и ( х, у, 2 ) — скалярный электрический потенциал,
= z
ôBZ (x, y,0)
ôx
By =-
ôФ(x, y, z) ô^(x, y,0 )
- z-
= z
8у 8у82
8ВХ (х, у,0) 8у
Как и для чисто электрического поля, согласно уравнениям (19), при малых отклонениях частицы от плоскости симметрии возникает сила, пропорциональная 2 , направление которой зависит от знака выражения I ( х, у )
симметрии
f
и(х,у,-2) = и(х,у,2) ; Ф(х,у,2) — скалярный магнитный потенциал, удовлетворяющий условию антисимметрии Ф(х, у, -2) = -Ф(х, у, 2) .
Как и раньше, заряженные частицы, стартовавшие в плоскости ОХТ , будут оставаться в ней, поскольку в плоскости ОХТ в силу антисимметричности скалярного магнитного потенциала Ф выполнены условия ВХ (х,у,0) = 0 , ВТ (х,у,0) = 0, а в силу симметричности скалярного электрического потенциала и — условие ЕХ (х, у,0) = 0 .
При малых отклонениях от плоскости ОХТ уравнения движения заряженной частицы вдоль координаты 2 принимают вид:
I =
+-c
ô2U ( x, y, 0) ô2U ( x, y,0)
ôx2
ôy
p(x, y )
ôBZ ( x, y, 0) ôy
- q ( x y)
ôBZ (x, y,0)
ôx
//
,(20)
где р ( х, у ) и д ( х, у) — скорости р ^) и д (t) , "приписанные" к точкам х ^), у (t) базовых траекторий, плотно заполняющих соответствующую подобласть плоскости ОХТ , так что каждой рабочей точке х, у соответствует некоторое допустимое множество скоростей. Поэтому, будет ли пучок, описываемый уравнениями (19), плоским целиком, определяется знаком величины (20).
z=0
1
z=0
А именно, для того чтобы пучок был плоским, т. е. локализованным в окрестности плоскости ОХУ, для положительно заряженных частиц достаточно выполнения условия
Ух,у,р(х,у),q(х,у): I <-е2 < 0, а для отрицательно заряженных частиц — условия
Ух,у,р(х,у),q(х,у): I > +е2 > 0 .
Строгое доказательство достаточности этого критерия совершенно аналогично электрическому случаю.
Как и в чисто электрическом случае, величина (20) зависит от данных электромагнитного поля, заданных исключительно в плоскости симметрии, — распределения электрического потенциала и0 (х,у) = и(х,у,0) и распределения нормальной компоненты индукции магнитного поля В^ (х,у) = В2 (х,у,0) . В отличие от чисто электрического поля, в выражение (20) входят еще и скорости заряженных частиц в соответствующих точках базовых траекторий, поэтому для одних семейств траекторий движение может оказаться плоским, а для других — нет.
КВАЗИСТАТИЧЕСКИЕ ЭЛЕКТРИЧЕСКИЕ И МАГНИТНЫЕ ПОЛЯ С ПЛОСКОСТЬЮ СИММЕТРИИ
Как легко понять, все эти выкладки остаются практически без изменений, если рассматриваются меняющиеся во времени симметричные электрические потенциалы и (х, у, г, t) и антисимметричные магнитные потенциалы Ф(х, у, 2, t), удовлетворяющие условию квазистатичности, т. е. когда напряжения на электродах и токи в катушках должны меняться достаточно медленно, чтобы в уравнениях для электромагнитного поля можно было пренебречь максвелловскими эффектами. Физический смысл "медленности" (а именно за время распространения электромагнитного возмущения — со скоростью света (!) — между электродами и катушками напряжения и токи для электродов и катушек должны меняться незначительно) показывает, что все используемые на практике электронно-оптические и ионно-оптические приборы, включая радиочастотные квадруполи и аналогичные им устройства, вполне могут рассматриваться в качестве квазистатических.
Посвящение. Лицензионная чистота
Авторы посвящают эту работу памяти Юрия Константиновича Голикова, создателя и бессменного руководителя аналитической школы оптики заряженных
частиц, к которой оба автора имеют честь и счастье принадлежать (в [33] — интернет-адрес списка трудов Ю.К. Голикова).
При выполнении моделирований траекторий и создании рисунков для данной публикации использовалась лицензионно-чистая программа Wolfram Mathematica версии 9.0 [34], приобретенная одним из авторов, и свободно распространяемая программа для редактирования растровых изображений Paint.NET версии 3.5 [35].
СПИСОК ЛИТЕРАТУРЫ
1. Голиков Ю.К., Уткин К.Г., Чепарухин В.В. Расчет элементов электростатических электронно-оптических систем. Учебное пособие. Л.: Изд-во ЛПИ им. М.И. Калинина, 1984. 80 с.
2. Голиков Ю.К., Соловьев К.В. Электростатические ионные ловушки. Санкт-Петербург: Изд-во Политехнического университета, 2008. 153 с.
3. Голиков Ю.К., Краснова Н.К. Теория синтеза электростатических энергоанализаторов. Санкт-Петербург: Изд-во Политехнического университета, 2010. 409 с.
4. Голиков Ю.К. Энергоанализирующие свойства электростатических полей с плоскостью симметрии. Дис. ... канд. физ.-мат. наук по спец. 01.04.04. Л.: ЛПИ им. М.И. Калинина, 1977. 105 с.
5. Чепарухин В.В. Некоторые вопросы теории синтеза электростатических дисперсионных энергоанализаторов. Дис. ... канд. физ.-мат. наук по спец. 01.04.04. Л.: ЛПИ им. М.И. Калинина, 1979. 132 с.
6. Голиков Ю.К. Определение электростатических полей по заданным характеристикам движения заряженных и дипольных частиц. Дис. ... д-ра физ.-мат. наук по спец. 01.04.04. Л.: ЛПИ им. М.И. Калинина, 1985. 254 с.
7. Спивак-Лавров И.Ф. Корпускулярная оптика статических ионно-оптических систем со средней плоскостью. Дис. ... д-ра физ.-мат. наук по спец. 01.04.04. Актюбинск, 1999. 232 с.
8. Соловьев К.В. Псевдооднородные электромагнитные структуры в инструментальной электронной оптике. Дис. ... канд. физ.-мат. наук по спец. 01.04.04. Санкт-Петербург: Санкт-Петербургский политехнический университет, 2004. 185 с.
9. Саулебеков А.О. Электронно-оптические свойства некоторых осесимметричных электростатических полей и разработка на их основе энергоанализаторов для исследования поверхности твердого тела. Дис. ... д-ра физ.-мат. наук по спец. 01.04.04. Алма-ты, 2007. 216 с.
10. Краснова Н.К. Теория и синтез диспергирующих и фокусирующих электронно-оптических сред. Дис. ... д-ра физ.-мат. наук по спец. 01.04.04. Санкт-Петербург, Санкт-Петербургский политехнический университет, 2014. 259 с.
11. Габдуллин П.Г., Голиков Ю.К., Краснова Н.К., Давыдов С.Н. Применение формулы Донкина в теории энергоанализаторов I // Журнал технической физики. 2000. Т. 70, № 2. С. 91-94.
12. Габдуллин П.Г., Голиков Ю.К., Краснова Н.К., Да-
выдов С.Н. Применение формулы Донкина в теории энергоанализаторов II // Журнал технической физики. 2000. Т. 70, № 3. С. 44-47.
13. Чаплыгин С.А. Новый метод приближенного интегрирования дифференциальных уравнений. М.—Л.: ГИТТЛ, 1950. 103 с.
14. Андронов А.А., Витт А.А., Хайкин С.Э. Теория колебаний. М., 1959. 916 с.
15. Андронов А.А., Леонтович Е.А., Гордон И.И., Май-ер А.Г. Качественная теория динамических систем второго порядка. М.: Наука, 1967. 568 с.
16. Андронов А.А., Леонтович Е.А., Гордон И.И., Май-ер А.Г. Теория бифуркаций динамических систем на плоскости. М.: Наука, 1967. 487 с.
17. Блакьер О. Анализ нелинейных систем. М.: Мир, 1969. 400 с.
18. Арнольд В.И. Геометрические методы в теории обыкновенных дифференциальных уравнений. Ижевск: РХД, 1999. 400 с.
19. Арнольд В.И., Афраймович В.С., Ильяшенко Ю.С., Шильников Л.П. Теория бифуркаций // Сб. "Динамические системы", серия "Итоги науки и техники. Современные проблемы математики. Фундаментальные направления". Т. 5. М.: ВИНИТИ, 1986. С. 5-218.
20. Марсден Дж., Мак-Кракен М. Бифуркация рождения цикла и ее приложения. М.: Мир, 1980. 368 с.
21. Крылов Н.М., Боголюбов Н.Н. Новые методы нелинейной механики. М.—Л.: Гостехтеориздат, 1934. 243 с.
22. Крылов Н.М., Боголюбов Н.Н. Введение в нелинейную механику. Киев: Изд-во АН УССР, 1937. 365 с.
23. Боголюбов Н.Н., Митропольский Ю.А. Асимптотические методы в теории нелинейных колебаний. М.: Наука, 1974. 503 с.
24. Моисеев Н.Н. Асимптотические методы нелинейной механики. М.: Наука, 1981. 379 с.
25. Бердников А.С. Управление транспортировкой заряженных частиц высокочастотными электрическими полями с квазидискретным спектром. Дис. ... д-ра физ.-мат. наук по спец. 01.04.01. Санкт-Петербург: Институт аналитического приборостроения РАН, 2013. 294 с.
26. Мэтьюз Дж., Уокер Р. Математические методы физики. М.: Атомиздат, 1972. 397 с.
27. Голиков Ю.К., Чепарухин В.В. Задача Коши для однородных гармонических потенциалов нулевой кратности // Сб. "Новые методы расчета ЭОС". Наука, 1983. С. 166-168.
28. Голиков Ю.К., Краснова Н.К. Электрические поля, однородные по Эйлеру, для электронной спектрографии // Журнал технической физики. 2011. Т. 81, № 2. С. 9-15.
29. Голиков Ю.К., Краснова Н.К., Абраменок О.А. Электрические спектрографы потоков заряженных частиц с потенциалами эйлерова типа // Прикладная физика. 2011. № 5. С. 69-73.
30. Голиков Ю.К., Краснова Н.К. Аналитические структуры электрических обобщенно-однородных спектрографических сред // Научное приборостроение. 2014. Т. 24, № 1. С. 50-58.
31. Уитеккер Э.Т., Ватсон Дж. Курс современного анализа. Часть 2: Трансцендентные функции. М.: ГИФМЛ, 1963. 516 с.
32. Арцимович Л.А., Лукьянов С.Ю. Движение заряженных частиц в электрических и магнитных полях. 2-е изд. М.: Наука, 1978. 224 с.
33. Список трудов проф. Ю.К. Голикова. URL: (http://iairas.ru/labs/mass_spectr_golikov_mem.pdi).
34. Программа для символьных вычислений Wolfram Mathematica. URL: (http://www.wolfram.com).
35. Свободно распространяемая программа для редактирования изображений Paint.NET. URL: (http://www.getpaint.net).
Институт аналитического приборостроения РАН,
Санкт-Петербург (Бердников А.С.)
Санкт-Петербургский государственный политехнический университет (Краснова Н.К.)
Контакты: Бердников Александр Сергеевич, [email protected]
Материал поступил в редакцию: 2.02.2015
ISSN 0868-5886
NAUCHNOE PRIBOROSTROENIE, 2015, Vol. 25, No. 2, pp. 69-90
SUFFICIENT CRITERIA FOR STABILITY AND NARROWNESS OF THE BAND-SHAPED ION BEAMS IN 3D ELECTRIC AND MAGNETIC FIELDS WITH PLANE OF SYMMETRY
1 2 A. S. Berdnikov , N. K. Krasnova
1 Institute for Analytical Instrumentation of RAS, Saint-Petersburg, Russia 2Saint-Petersburg State Polytechnic University, Russia
The problem of stability of the charged particle beams is considered when they are shifted outside the plane of symmetry z = 0 of the electric and/or magnetic field. It is shown that if the electrostatic potential U(x, y) defined along the plane of symmetry z = 0 satisfies the relation Uxx + Uyy < -s2 < 0 then the trajectories of charged particles moving in this electric field are located in the vicinity of the plane of symmetry z = 0 providing the perturbations of the initial conditions of the particles are small enough. Moreover, it is possible to make constructive evaluations how far the trajectory can go outside the plane of symmetry z = 0 for specified perturbations of the initial conditions. The said sufficient criterion for the stability and the compact size of the beam is generalized for magnetostatic fields with a plane of symmetry and for combined electrostatic and magnetostatic fields possessing the same plane of symmetry.
Keywords: charged particle optics, compact beams of charged particles, stability of the motion, optimal synthesis of electron-optical systems
REFERENCES
1. Golikov Yu.K., Utkin K.G., Cheparukhin V.V. [Calculation of the elements of electrostatic electron optical systems. Textbook]. L., Leningrad Polytechnic Institute Publ., 1984. 80 p. (In Russ.).
2. Golikov Yu.K., Solov'ev K.V. [Electrostatic ion traps]. St.Petersburg, St.Petersburg Polytechnic University Publ., 2008. 153 p. (In Russ.).
3. Golikov Yu.K., Krasnova N.K. [The theory of electrostatic energy analyzers synthesis]. St.Petersburg, St.Petersburg Polytechnic University Publ., 2010. 409 p. (In Russ.).
4. Golikov Yu.K. [Energy analyzing properties of electrostatic fields with plane of symmetry. Doct. diss.]. L., Leningrad Polytechnic Institute, 1977. 105 p. (In Russ.).
5. Cheparukhin V.V. [Some problems in the theory of synthesis of electrostatic dispersive energy analyzers. Doct. diss.]. L., Leningrad Polytechnic Institute, 1979. 132 p. (In Russ.).
6. Golikov Yu.K. [Finding of the electrostatic fields using the specified characteristics of the motion of charged and dipole particles. Dr. phys. and math. sci. diss.]. L., Leningrad Polytechnic Institute, 1985. 254 p. (In Russ.).
7. Spivak-Lavrov I.F. [Charged particle optics of static ion-optical systems with the median plane. Dr. phys. and math. sci. diss.]. Aktobe, 1999. 232 p. (In Russ.).
8. Solov'ev K.V. [Pseudo uniform electro magnetic structures for electron optical instruments. Doct. diss.]. St.Petersburg, St.Petersburg Polytechnic University, 2004. 185 p. (In Russ.).
9. Saulebekov A.O. [Electron optical properties of some
axially symmetrical electrostatic fields and its applications to the development of energy analyzers to study the solid surfaces. Dr. phys. and math. sci. diss.]. Almaty, 2007. 216 p. (In Russ.).
10. Krasnova N.K. [Theory and synthesis of dispersing and focusing electron optical media. Dr. phys. and math. sci. diss.]. St.Petersburg, St.Petersburg Polytechnic University, 2014. 259 p. (In Russ.).
11. Gabdullin P.G., Golikov Yu.K., Krasnova N.K., Davydov S.N. [Application of Donkin's formula in the theory of energy analyzers: Part I]. Zhurnal technicheskoy fiziki [The Russian Journal of Applied Physics], 2000, vol. 70, no. 2, pp. 91-94. (In Russ.).
12. Gabdullin P.G., Golikov Yu.K., Krasnova N.K., Davydov S.N. [Application of Donkin's formula in the theory of energy analyzers: Part II]. Zhurnal technicheskoy fiziki [The Russian Journal of Applied Physics], 2000, vol. 70, no. 3, pp. 44-47. (In Russ.).
13. Chaplygin S.A. [The new method for approximate integration of differential equations]. M., L., GITTL Publ., 1950. 103 p. (In Russ.).
14. Andronov A.A., Vitt A.A., Khaikin S.E. Theory of oscillators. Pergamon Press, Oxford, 1966. 482 p. (Russ. ed.: Andronov A.A., Vitt A.A., Khaikin S.E. Teoriya kolebaniy. M., 1959. 916 p.).
15. Andronov A.A., Leontovich E.A., Gordon I.I., Maier A.G. Qualitative theory of second-order dynamic systems. Halsted Press, New York, 1973. (Russ. ed.: Andronov A.A., Leontovich E.A., Gordon I.I., Maier A.G. Kachestvennaya teoriya dinamicheskich sistem vtorogo poryadka. M., Nauka Publ., 1967. 568 p.).
16. Andronov A.A., Leontovich E.A., Gordon I.I., Maier A.G. Theory of bifurcations of dynamic systems on a plane. Halsted Press, New York, 1973. (Russ. ed.:
Andronov A.A., Leontovich E.A., Gordon I.I., Maier A.G. Teoriya bifurkaziy dinamicheskich sistem naploskosti. M., Nauka Publ., 1967. 487 p.).
17. Blaquiere A. Nonlinear system analysis. London, Academic Press, 1966. (Russ. ed.: Blaquiere A. Analiz nelineynych system. M., Mir Publ., 1969. 400 p.).
18. Arnold V.I. Geometrical methods in the theory of ordinary differential equations. Springer-Verlag, Berlin, 1988. (Russ. ed.: Arnol'd V.I. Geometricheskie metody v teorii obyknovennych differenzial'nych uravneniy. Izhevsk, RChD, 1999. 400 p.).
19. Arnol'd V.I., Afraymovich V.S., Il'yashenko Yu.S., Shil'nikov L.P. [Theory of bifurcations]. Sb. "Dinamicheskie sistemy", seriya "Itogi nauki i techniki. Sovremennye problemy matematiki. Fundamental'nye napravleniya" [Collection of articles "Dynamic systems"], vol. 5, M., VINITI Publ., 1986. pp. 5-218. (In Russ.).
20. Marsden J.E., McCracken M. The hopf bifurcation and its applications. Ser. "Applied Mathematical Sciences", vol. 19, Springer-Verlag , New York, 1976. 409 p. (Russ. ed.: Marsden Dzh., Mak-Kraken M. Bifurkaziya rozhdeniya zikla i eeprilozheniya. M., Mir Publ., 1980. 368 p.).
21. Krylov N.M., Bogolyubov N.N. Novye metody nelinejnoj mekhaniki [New methods in nonlinear mechanics]. M.-L., Gostechizdat Publ., 1934. 243 p. (In Russ.).
22. Krylov N.M., Bogolyubov N.N. Introduction to nonlinear mechanics. Princeton University Press, Princeton, 1947. (Russ. ed.: Krylov N.M., Bogolyubov N.N. Vvedenie v nelineynuyu mechaniku. Kiev, AN USSR Publ., 1937. 365 p.).
23. Bogolyubov N.N., Mitropolsky Yu.A. Asymptotic methods in the theory of nonlinear oscillations. Gordon and Breach/Hindustan Publ., New York, Delhi, 1961. (Russ. ed.: Bogolyubov N.N., Mitropolsky Yu.A. Asimptoticheskie metody v teorii nelineynych ko-lebaniy. M., Nauka Publ., 1974. 503 p.).
24. Moiseev N.N. [Asymptotic methods of nonlinear mechanics]. M., Nauka Publ., 1981. 379 p. (In Russ.).
25. Berdnikov A.S. [Control of charged particle transportation by means of high frequency electric fields with quasi-discrete spectrum. Dr. phys. and math. sci.
Contacts: Berdnikov Aleksandr Sergeevich, [email protected]
diss.]. St.Petersburg, Institute for analytical instrumentation, 2013. 294 p. (In Russ.).
26. Mathews J., Walker R.L. Mathematical methods of physics. W.A. Benjamin Inc., New-York—Amsterdam, 1964. (Russ. ed.: Met'yuz Dzh., Uoker R. Matematicheskie metody fiziki. M., Atomizdat Publ., 1972. 397 p.).
27. Golikov Yu.K., Cheparukhin V.V. [Cauchy's task for uniform harmonious potentials of zero frequency rate]. Sb. "Novye metody rascheta EOS" [New methods of calculation of electron-optical systems], Nauka Publ., 1983. pp. 166-168. (In Russ.).
28. Golikov Y.K., Krasnova N.K. [Application of electric fields uniform in the Euler sense in electron spectrography]. Zhurnal technicheskoy fiziki [The Russian Journal of Applied Physics], 2011, vol. 81, no. 2, pp. 9-15. (In Russ.).
29. Golikov Yu.K., Krasnova N.K., Abramyonok O.A. [Electric spectrographs for the flows of charged particles using the Euler' type potentials]. Prikladnaya fizika [Applied physics], 2011, no. 5, pp. 69-73. (In Russ.).
30. Golikov Yu.K., Krasnova N.K. [Analytical frameworks for electric generalized-homogeneous spectrographic media]. Nauchnoe priborostrienie [Science Instrumentation], 2014, vol. 24, no. 1, pp. 50-58. (In Russ.).
31. Whitakker E.T., Watson G.A. Course of modern analysis. Cambridge, Cambridge University Press, 1927. (Russ. ed.: Uitekker E.T., Vatson Dzh. Kurs sovremennogo analiza. Chast' 2: Transzendentnye funkzii. M., GIFVL Publ., 1963. 516 p.).
32. Artzimovich L.A., Lukyanov S.Yu. [Motion of charged particles in electric and magnetic fields. 2nd edition]. M., Nauka Publ., 1978. 224 p. (In Russ.).
33. List of works of the prof. Yu. K. Golikov. URL: (http://iairas.ru/labs/mass_spectr_golikov_mem.pdf).
34. The software for symbolic computations Wolfram Mathematica. URL: (http://www.wolfram.com).
35. The free editor for bitmap images Paint.NET. URL: (http://www.getpaint.net).
Article received in edition: 2.02.2015
