Меняющийся во времени псевдопотенциал и его применение к описанию усредненного движения заряженных частиц. Ч. 2. Общая формула Текст научной статьи по специальности «Математика»

ISSN 0868-5886

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2011, том 21, № 3, c. 83-96

МАТЕМАТИЧЕСКИЕ МОДЕЛИ

УДК 537.534.7: 621.319.7 © А. С. Бердников

МЕНЯЮЩИЙСЯ ВО ВРЕМЕНИ ПСЕВДОПОТЕНЦИАЛ И ЕГО ПРИМЕНЕНИЕ К ОПИСАНИЮ УСРЕДНЕННОГО ДВИЖЕНИЯ ЗАРЯЖЕННЫХ ЧАСТИЦ. Ч. 2. ОБЩАЯ ФОРМУЛА

В этом цикле работ рассматривается новый метод управления движением заряженных частиц с помощью высокочастотных электрических полей, характеризующихся медленноменяющимися во времени эффективными потенциалами. Данная работа посвящена выводу общей формулы для медленноменяющегося псевдопотенциала. Показано, что классический эффективный потенциал является лишь главным членом ряда для эффективного гамильтониана, разложенного в двойной ряд: по малому параметру, характеризующему амплитуду высокочастотного электрического поля, и по обратным частотам высокочастотного электрического поля. Получена формула для эффективного потенциала (главного члена ряда для эффективного гамильтониана), позволяющего приближенно описывать усредненное движение заряженных частиц в высокочастотных электрических полях, характеризуемых "быстрым" и "медленным" временем характерного изменения поля.

Кл. сл.: масс-спектрометрические приборы, радиочастотные приборы, эффективный потенциал, гамильтонова динамика, гамильтоновы динамические системы, техники усреднения динамических гамильтоновых уравнений, высокочастотные электрические поля, радиочастотные электрические поля

ВВЕДЕНИЕ

В предыдущей работе этого цикла рассматривались теоретические основы нового способа управления движением заряженных частиц с помощью высокочастотных электрических полей. А именно: с помощью специальным образом организованных высокочастотных полей можно захватывать в локальные зоны устойчивости заряженные частицы с разными массами, энергиями и зарядами и синхронно перемещать указанные зоны захвата заряженных частиц в выбранной области пространства. Захваченные заряженные частицы перемещаются синхронно с перемещением областей захвата заряженных частиц, и тем самым осуществляется целенаправленное управление транспортировкой заряженных частиц. Характерной особенностью указанного способа является то, что в локальной области захвата объединяются как положительно, так и отрицательно заряженные частицы.

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

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

1. "МЕДЛЕННЫЕ" И "БЫСТРЫЕ" ФУНКЦИИ

Придадим конкретный математический смысл, что будет пониматься далее под "медленными" и "быстрыми" временами. Пусть имеются верхняя граница "медленных" частот 5 и нижняя граница "быстрых" частот А . Очевидным образом предполагается, что А»5. Функция h(t) называется

"медленной", если ее спектр равен нулю (или пренебрежимо мал) вне интервала частот о е (-5, +5) . Функция Н (t) называется "быстрой", если ее спектр равен нулю (или пренебрежимо мал) на интервале частот ое(-А, +А).

Т. е. функция Н (t) содержит только высокочас-

тотные компоненты и не содержит низкочастотных, а функция h (t) содержит только низкочастотные компоненты и не содержит высокочастотных. Интервалы частот для низкочастотных и для высокочастотных компонентов не пересекаются и достаточно далеко разнесены друг от друга.

Указанные ограничения на спектр функций влекут за собой справедливые "в среднем" неравенства:

^(г)/dt|У|h(t)|2 < 52 и (t)/dt| 2/\н (г )|2 >д2,

где усредненный квадрат модуля вещественного временного сигнала f (t) получается в результате

2 1 '+Т

предельного перехода ^(t)|2 ~ 11т — | f2 (т)dт.

ТТ г-т

В силу указанных соотношений для функций Н (г) и h (г) в выражении d (Н (г) h (г))/dt = = h (г) dH (г)/dt + Н (г) dh (г)Idt для полной производной по времени от произведения Н (г) h (г) вклад члена h (г) 6Н (г)/dt существенно превышает вклад члена Н (г) dh (г)/dt, поскольку выполняется соотношение

I 12 /| 12 I \2 !\ |2

|аН (г/|Н (г)| >> |dh(г/|h(г)| .

Этот факт может использоваться при процедуре усреднения уравнений движения для разделения членов уравнения, имеющих разные порядки малости.

Очевидно, что класс "медленных" функций образует линейное пространство (тождественный нуль — тоже "медленная" функция). Более того, произведение двух "медленных" функций также является "медленной функцией, но с более слабым показателем "медленности": после преобразования Фурье произведение функций превращается в свертку спектров, так что спектр произведения двух "медленных" функций будет равен нулю вне интервала частот юе(-25, +25). Соответственно спектр произведения N "медленных" функций будет равен нулю вне интервала частот ®е(-N5, + N5) .

В отличие от "медленных" функций, "быстрые" функции, как они определены выше, представляют собой более сомнительный математический объект. Класс "быстрых" функций не является линейным пространством (отсутствует нуль). Сумма двух "быстрых" функций может содержать в неявном виде "медленное" время:

cos(ю - 5)t + cos(a + 5)t = 2cosStcosat.

Произведение двух "быстрых" функций может быть как "быстрой" функцией (cos cot • sin cot = = (1/2)sin2at), так и содержать "медленную" компоненту (cos cot • cos cot = (1/2) - (1/2)cos2at).

Произведение "медленной" функции и "быстрой" функции является с формальной точки зрения "быстрой" функцией: в силу превращения произведения в свертку спектров, спектр произведения равен нулю на интервале частот

юе(-(А- 5),+(Д- 5)) , однако выделение "медленных" временных зависимостей на фоне "быстрых" как раз и является нашей задачей.

Поэтому в качестве "быстрых" функций мы будем рассматривать более узкий класс — набор гармоник cos akt и sin a>kt, к которому для полноты базиса добавлена "медленная" функция 1. Час-будут предполагаться "быстрыми"

тоты

c

и

"далеко

разнесенными"

c

будет

(Ы >Д )

(Vrn Ф n: |ют -юп\ >А). Набор частот предполагаться фиксированным, однако в процессе работы к нему по мере необходимости могут добавляться суммарные и разностные частоты вида ют ± юп.

^ т п

Предполагается, что расширенный набор частот по-прежнему удовлетворяет условию Vm Ф n : |ют - an | > А, однако при наличии несоизмеримых частот и при повторении указанного шага достаточное число раз это условие может начать нарушаться. Наиболее безопасным в этом смысле является набор кратных "быстрых" гармоник юк = ka , который изначально содержит в себе все возможные суммарные и разностные частоты и при этом удовлетворяет условию, что все частоты являются "далеко разнесенными" (Vm Ф п : |ют - юп | > a > А). Тем не менее, если добавление суммарных и разностных частот ограничивается несколькими шагами, вполне допустимо рассматривать наборы "быстрых" частот юк, которые не являются кратными.

В центре нашего внимания, однако, будет находиться отличный от этих двух классов третий класс функций — "быстро-медленные" функции времени, которые представимы в виде

A (t ) + Х( Ак (t) cos akt + Bk (t) sin akt), где Ak (t) и Bk (t) — медленные функции, а ak — фиксированный набор быстрых и далеко разнесенных частот. Такие функции образуют линейное пространство, а произведение этих функций принадлежит к подобному же классу, но с ослабленным параметром медленности и с расширенным набором час-

тот ak. Физически данный класс функций представляет собой функции времени, спектр которых отличен от нуля на далеко разнесенных и узких интервалах a е (ak - 8,ак + S) с фиксированными центрами ak и равен нулю (или пренебрежимо мал) для всех остальных частот. В следующей работе будет показано, что функции времени f (t) , у которых спектр обладает указанным свойством, действительно всегда могут быть приведены к виду

A0 (t) + X (Ak (t) cos Ckt + Bk (t) sin Ckt) .

2. МЕТОД УСРЕДНЕНИЯ УРАВНЕНИЙ ДВИЖЕНИЯ КРЫЛОВА—БОГОЛЮБОВА КАК СПОСОБ ПОЛУЧЕНИЯ ФОРМУЛ ДЛЯ ЭФФЕКТИВНОГО ПОТЕНЦИАЛА

Использование высокочастотных электрических полей для манипулирования заряженными частицами в масс-спектрометрии имеет давнюю традицию. Обзор основных классов подобных устройств и описание основных закономерностей их функционирования может быть найден, например, в [1-3].

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

Количественно данный эффект описывается с помощью теории эффективного потенциала, или псевдопотенциала, впервые введенного, по-видимому, П.Л. Капицей [4-6]. А именно, если частота со осцилляций электрического поля E (х, y, z, t) , происходящих по закону

E (х, y, z, t) = E0 (х, y, z) cos (at + p),

достаточно высока, а смещение заряженной частицы с массой m и зарядом e за один период изменения электрического поля мало, то движение заряженной частицы можно представить в виде "усредненного", или "медленного" движения, на которое наложено быстроосциллирующее движение, имеющее, впрочем, маленькую амплитуду. Уравнение для усредненного движения выглядит так, как будто движение происходит в электрическом поле с потенциалом U (х, y, z ) =

= e|E0 (х,y,z)| I(4ma2) . В более общем случае

высокочастотное электрическое поле представляется в виде ряда Фурье, как

E ( х, y, z, t) =

= X (E(k) (х, y, z) cos (kat) + E<k) (х,y, z)sin (kat)),

и тогда эффективный потенциал должен вычисляться по формуле

U ( х, y, z ) =

= eX (|E^k) (х,y, z)|2 + Ek) (х, y, z)|2)/(4ma2k2).

Если же высокочастотное электрическое поле имеет вид

E ( х, y, z, t) =

= X ( Ek) ( х, y, z ) cos akt + Ef) ( х, y, z ) sin akt),

где ak — произвольный набор частот, то эффективный потенциал, вычисляемый по формуле

U ( х, y, z ) =

= eX(|ECk) (х,y,z)|2 + E) (х,y,z)|2)/(4ma¡),

будет правильно описывать усредненное движение заряженных частиц только тогда, когда среди частот ak нет близких друг к другу. Подробности и обоснование теории эффективного потенциала можно найти в упомянутых выше ссылках [4-6], некоторые тонкие моменты указанной теории рассматриваются также в [7-13].

Для обобщения классической теории эффективного потенциала на случай высокочастотных электрических полей, характеризуемых быстрыми и медленными характерными временами изменения электрического поля, будет использоваться метод усреднения Крылова—Боголюбова [14-16] в соответствии с рекомендациями, указанными в [7]. Несмотря на то что получающиеся в результате формулы можно обосновать также более простым способом с помощью нестрогих и интуитивных соображений по аналогии с [4, 5], имеет смысл воспользоваться строгим и математически корректным подходом.1"1

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

Пусть дано некоторое устройство для управления движением заряженных частиц, обладающее набором электродов и создающее в объеме движения заряженных частиц высокочастотное электрическое поле при приложении к электродам устройства высокочастотных электрических напряжений. Будем считать, что сигналы, прикладываемые к электродам, относятся к специальному классу, характеризуемому медленными и быстрыми временами эволюции сигнала во времени, так что высокочастотное электрическое поле, возникающее в объеме устройства, может быть представлено, как

E ( x, y, z, t) = E0 ( x, y, z, t) +

+Z (Ekc) ( x, У, z, t) cos akt + Ek*) ( x, y, z, t) sin akt),

где компоненты высокочастотного электрического поля E0 (x, y, z, t), E(k} (x, y, z, t) и Ek*) (x, y, z, t) являются медленными функциями времени, а ak — быстрыми и далеко разнесенными частотами. (В следующей работе будет показано, как от электрических сигналов, характеризуемых медленными и быстрыми временами, можно перейти к такому представлению для электрического поля).

Опишем вкратце процедуру усреднения га-мильтоновых уравнений движения, используемую в данной работе. Для системы уравнений движения вводится каноническая замена гамильтоновых переменных

r = w + sR1 (w, v, t) + s2 R 2 (w, v, t) +...,

p = v + sP1 (w, v, t) + s2P2 (w, v, t) +...,

переводящая одни гамильтоновы уравнения в другие гамильтоновы уравнения (здесь (w, v) — новые координаты и импульсы, (r,p) — исходные координаты и импульсы). Функции Rk (w, v, t) и Pk (w, v, t) содержат зависимость от "быстрого" времени в виде быстроосциллирующих функций с нулевым средним значением, поэтому по своему физическому смыслу (w, v) — это усредненные функции (r,p , опреденные с помощью указанной замены переменных математически корректным способом. Функции Rk (w, v, t) и Pk (w, v, t) выбираются так, чтобы после использования этой подстановки результирующий гамильтониан не содержал "быстрого" времени (хотя, возможно, в нем останется зависимость от "медленного" времени). В результате применения указанного подхода динамические уравнения для усредненных переменных получаются в виде гамильтоновых

уравнений движения с медленно меняющимся во времени гамильтонианом, которым подчиняются усредненные координаты и импульсы (w, V) . Все инварианты и законы сохранения, справедливые для исходных уравнений, автоматически сохраняются для преобразованных уравнений, а сами усредненные уравнения оказываются автоматически представлены в гамильтоновой форме.

Возникает естественный вопрос: существуют ли требуемые нам замены переменных? Ответ на этот вопрос является положительным. Действительно, решение любой гамильтоновой системы может рассматриваться, как зависящее от времени преобразование начальных точек фазового пространства в конечные точки фазовой траектории. Указанная замена переменных будет канонической, т. е. обеспечивать преобразование любых исходных гамильтоновых уравнений в новые га-мильтоновы уравнения [17]. Преобразование га-мильтоновых координат, обратное к каноническому, тоже будет каноническим. Однако, если взять в качестве гамильтоновой системы, порождающей такое каноническое преобразование, саму систему уравнений движения заряженной частицы, то после преобразования в качестве системы уравнений должна получиться система уравнений с нулевыми правыми частями (начальные координаты траектории, выраженные через конечные точки траектории, являются первыми интегралами). Тем самым по крайней мере одна такая замена переменных существует, хотя ее нахождение будет эквивалентно решению исходной системы динамических уравнений в общем виде, и совсем не обязательно, что такая замена переменных будет оптимальной для наших целей.

Достаточно трудоемкие выкладки по нахождению приближенных выражений для требуемой нам замены переменных будут приведены в следующих разделах данной работы. В результате этих вычислений получим, что в первом приближении усредненное движение заряженной частицы представляет из себя движение в электрическом потенциале ио ( х, у, г, г) + и ( х, у, г, г), где

ио (х, у, г, г) — это медленный электрический потенциал, соответствующий низкочастотной компоненте электрического поля Е0 (х, у, г, г), а

и (х, у, г, г) — это медленный эффективный потенциал, который выражается через медленные амплитуды быстрых гармоник электрического поля, как

и ( х, у, г, г) =

= е!(|ЕС (х,у,г,г)|2 + Е) (х,у,г,г)|2)/(4шш1).

Данное выражение является прямым аналогом классического эффективного потенциала, обобщенного на случай быстро осциллирующего электрического поля, характеризуемого медленным и быстрым временами.

3. СИСТЕМА УРАВНЕНИЙ ДЛЯ УСРЕДНЯЮЩЕЙ ЗАМЕНЫ ПЕРЕМЕННЫХ

Рассмотрим движение заряженной частицы в быстро осциллирующем электрическом поле указанного вида. Из предположения, что напряженность высокочастотного электрического поля может быть записана в форме E (x, t ) =

= E0 (x, t) + (x,t) cos ((0/) + E(s) (x,t)sin((0/),

следует, что потенциал данного электрического

x,

поля имеет

(c)

вид U ( x, t ) = U0 ( x, t) + +S^Ujc) (x,t)cos(t) + U(f (x,t)sin((Ojt), где гармоники потенциала Uk (x, t) — медленные функции, (k — быстрые и далеко разнесенные друг от друга частоты, s — малый параметр, а x = ( х, y, z) — вектор, характеризующий положение частицы в пространстве.

Гамильтониан, описывающий движение заряженной частицы в таком электрическом поле, F(x,p,t) = F0(x,p,t) + sFa (x,p,t) + s2Fb (x,p,t) +...

новая система уравнений была гамильтоновой. Дифференциальные соотношения, которым должны в самом общем случае удовлетворять эти функции, приводятся, например, в [17]:

ЕРк5хк - G(^P,t)St =

- П .(Е РкЗХк - Н ( X, P, t ^ )- SF ( x, p, t),

где G(x,p,t) — новый гамильтониан, Н(X,P,t) — гамильтониан до замены переменных, F (x, p, t) — некоторая свободная функция, П — константа, а значок S в данном контексте обозначает варьирование функции по переменным x, p и t.

После замены переменных исходный гамильтониан приобретает вид

Н ( X, Г, t ) =

= Н0 (x,р,t) + еНа (x,р,t) + е2Нъ (x,р,t) + ...,

где е — малый параметр. Естественно искать решение в виде

G (x, р, t) = G0 (x, р, t) + еGа (x, р, t) + е2Gъ (x, р, t) + ..., П = По + еПа + е2 Пъ + ...,

имеет вид

H=¿ Е p )+eUо (хt)+

+se^(Uf (X,t)cos(Ojt + Ujs) (X,t)sinm}t),

где X — исходные пространственные координаты, P — соответствующие им импульсы, t — время, m — масса частицы, e — заряд частицы. Замена переменных записывается в форме:

Хк = Хк + еХа (x, р, t) + е2Хьк (x, р, t) +..., Рк = Рк + еРка (x, р, t) + е2Ркъ (x, р, t) + ..,

где x — усредненные координаты X, р — усредненные импульсы Р , а функции Х£ (x,р, t) ,

Хкъ (x,р, t), Ра (x,р, {), Ра (x,р, {), подлежат определению как функции независимых переменных x, р и t .

Указанная замена переменных не может быть произвольной, чтобы после замены переменных

2) Удобнее обозначить исходные переменные большими буквами, а усредненные переменные — маленькими буквами, поскольку вся основная работа будет происходить в пространстве усредненных переменных.

Подставив эти выражения в исходные дифференциальные соотношения и сгруппировав члены с разными степенями малого параметра при независимых дифференциалах Sxk, Spk и St, получим систему уравнений относительно частных производных д^ /дхк , /дрк , д^ /дt, дГа/дхк ,

дра/дРк , д , дхк , дръ/дРк , дt, стоящих в левой части. Эти вспомогательные уравнения приводятся далее в разделах 3.1, 3.2, 3.3. Разрешимость их требует отдельного анализа, который там и будет выполнен.

Отметим, однако, полезный прием, который будет многократно использоваться нами в дальнейшем. Пусть дана система уравнений в частных производных дF/дtk = Тк (t), где t = t2,...) —

независимые переменные, F (t) — неизвестная функция, Тк (t) — известные правые части, не зависящие от неизвестной функции F (t) . Тогда

система уравнений имеет решение, причем единственное с точностью до аддитивной константы, тогда и только тогда, когда Ут, п выполнены соотношения дТ/ дtn = дТ/ дtm .

т / п п / т

Указанные соотношения следуют из перестановочности частных производных д(дF|дtm )/дtn =

= d(dF/dt„)/dtm , и поэтому они очевидным образом являются необходимыми для разрешимости данной системы уравнений. Достаточность же этих условий следует, например, из общей теории совместности систем уравнений в частных производных [18-20].

Однако в данном случае достаточность можно легко доказать и без привлечения общей теории. Выберем одно из уравнений, например уравнение W/dtm = Tm (t) . Поскольку Vk dTm/dtk =STk/dtm ,

то после замены F(t) = F(t) -jTm (t)dtm в уравнениях dF/dtk = Tk (t) -j(öTm (t)/ötk)dtm, полученных для новой функции F, правые части

Tk = Tk -j(dTm/dtk )dtm не зависят от tm , для новых уравнений снова справедливы условия

ÖTk/ 5tn = dT„/ Ötk И, кроме TO^ dTkl 5tm = 0,

dFIdtm = 0. Последовательно исключая переменные t1, t2, ..., получим в конце концов одно уравнение с одной свободной переменной, которое имеет единственное решение с точностью до аддитивной константы.

Тем самым мы не только доказали разрешимость системы уравнений dF/dtk = Tk (t) при выполнении соотношений 6T /dt =6T /5t для

ml „ „Im ^

Vm,„ , но и нашли решение в явном виде.

3.1. Уравнения нулевого порядка для усредняющей замены переменных

Уравнения нулевого порядка имеют вид

5F0/ dpk = 0, 5Fjdxk =(D0-l)pk , dF0/dt = = G0 - D0H0. Из равенства смешанных производных d(dF0/dPk )/dxk = d(dF0/dxk )/dPk следует условие D0 = 1. Из равенства смешанных производных d(dF0/dPk)/dt = d(dF)/dt)/dp, и d(dF0/dxk)/dt = = d( dF0 / dt) / dxk получаются уравнения

dG()/ dPk = dH0l dPk , dG^ dxk = dH0/dxk , имеющие очевидное решение

G0 (x, p, t) = H0 (x, p, t) + const (t),

где произвольная аддитивная добавка const (t) никак не влияет на динамические уравнения, порождаемые гамильтонианом G (x,p, t) . Решение

F0 (x,p, t) = 0 получившейся системы уравнений,

не имеющее, впрочем, большого самостоятельного значения, завершает анализ уравнений нулевого порядка.

3.2. Уравнения первого порядка для усредняющей замены переменных

Уравнения первого порядка имеют вид:

дк/ дР„=Е Рк дХа/ дрп,

к

дFa/дхт = Р, + Dapm +ХРк дХЦ/дхт ,

к

дК/ дt = Ga - DaH0 - На +Х Рк дХка/дt.

к

Здесь учитывается условие D0 = 1, полученное ранее. Смешанные производные

д{дК1 дР„ )/дРт = д(дF0 /дРт )/дР„

дают условие дХ„а /дРт = дХат/дР„ . Смешанные производные

д^а/ дх„) дт = д(дFo/дх т )/дх „

дают условие дР„а /дхт = дР^/дх„ . Смешанные производные

д(д-^ /дхт )/дРт = д(дF0 /дРт )/^т

дают условие Da + дХ,/ дхст + дРта / дРт = 0, а смешанные производные

д(д-^ /^т )/дР„ = д(дF0 /дР„ )/^т —

условие дХ„а / дхт +дPma / ф„ = 0, т Ф „..

Из этих условий следует, что система уравнений

а^/ дРк = ХГ + Daxk| 2, д«,/дхк =-Рка -Dapк|2

относительно неизвестной функции (x, p, t) при заданных функциях Х, и Рк оказывается совместной и имеет единственное решение. С другой стороны, если замена переменных определена, как Х, =д«^дРк -Daxk¡2, Рка =-дБа1дхк -Dapk|2, где (x, p, ^ — произвольная функция, то автоматически выполнены нужные нам условия совместности для уравнений первого порядка.

С учетом этой подстановки оставшиеся неиспользованными смешанные производные

д^а/дt)/дРт =д(дР0/дРт)/дt и д^/д^/дхт =

= д(дF0/ дхт)/дt дают нам уравнения для функции Ga (x,I), 0:

д£а/ дРк = Da дН 0/дРк +дНа/дРк -д 2 дРк & , дСа/дхк = Da дН,/дхк + дН,,/дхк - д2дхкдt ,

очевидным решением которых является

аа (^ р, t) =

= ПаНо (X, р, t) + На (X, р, t) - д^а (^ р, t)/^ .

Правильным подбором функции Sа (x,р, t) и константы Па можно добиться, чтобы новый гамильтониан Оа (x,р, t) имел удобный для нас вид.

3.3. Уравнения второго порядка для усредняющей замены переменных

Уравнения второго порядка имеют вид

дРа/ дРп =

= ЕРк дХЦдРп +ЕРа дХка/дРп + ПаЕРк дХка1 дРп,

дРа/дХт = РЪт + ПаРт + ПъРт +

+Е Рк дХЦ дХт +Е Ра дХкадт + Па Е Рк дХка/дХт,

дFa|дt = аа -ПЬНо -ПаНа -Нъ +

+Е Рк дХЦ дt + Е Ра дХка / дt + Па Е Рк дХка / дt..

Смешанные производные д(дFa| Срп )/дРт = = д ( дF0 /дРт )/дРп дают условие

дХЦ дРт + Па дХа/дРт +Е(дР^/дРп )(Х/дРт ) =

к

= дХЦ дРп + Па дХт/дРп +Е(дРка/дРт )(дХка/дРп ).

к

Смешанные производные д(дFa| дхп )/дхт = = д ( дF0 /дхт Удхп дают условие

дРЫдхт + Па дР^дхт + Е^к^ )(дХа/дхп) =

к

= дРЫ дхп + Па дРа/дхп + Е(дРка/дхп )№/ дхт ).

к

Смешанные производные д(дFa| дхт )/дРт = = д(дF0 /дРт )/дхт дают условие

Пъ + ПадРта/дРт + Па дХ°т/дхт + д^/дРт +

+ дХЦдхт = Е(дРка/дхт )(дХа/дРт) -

к

-Е(дРка/дРт)(дХка/дхт) ,

к

а смешанные производные д(дРа/ дхт )/дРп = = д ( дF0 / дРп У дхт — условие

Па дРта /дРп + Па дХа /дхт + дРЫ дРп + дХЦ дхт =

= Е(дРка/ дхт )(дХка/дРп) -

к

-Е(дРка/дРп)(дХка/дхт) .

к

Можно проверить, что когда для функций Хъ и Ръ выполнены указанные здесь условия совместности, то система уравнений д£ъ/ Срп = Хъ - Мп, д£л/х = -Ръ + N относительно неизвестной

ъ т т т

функции Sъ (x, р, t) также оказывается совместной и имеет единственное решение, где функции Мп и выражаются через функции Ха и Рка по формулам

Мп=(12)Е(Ра дХа /дРп - дРка/дРп Ха) -

к

-ПъХпа -Пъхп/2,

и ^т = (1/2)Е(дРка/дхтХка - Рка дХаа/дхт )-к

- ПъРа -Пъ п /2.

ъ т ъ т

Точно также легко проверить, что если замена переменных определена как Хъ =д8ъ/Срп + Мп,

Рт =-д<^/дхт + ^т , где SЪ (X, P, t) — произвольная функция, а функции Мп и выражаются через функции Х а и Р а , как указано выше, то нужные нам условия совместности для уравнений второго порядка будут выполнены автоматически.

С учетом этой подстановки, оставшиеся неиспользованными смешанные производные

д(дРъ/дt ^дРт =д(дРъ/дРт )/дt и д(дРъ/а )/дхт =

= д(дРъ/ дхт дают нам для функции аъ (x, р, t) систему уравнений, для которых условия совместности оказываются всегда выполнены. Решением этих уравнений является функция

^ = ПЪН 0 + ПаНа + Нь-дsъ| а +

+(1/2)Е(ХадРка/дt - Ра дХа!^).

к

Правильным подбором функции Sъ (x,р, t) и константы Пъ можно добиться, чтобы новый гамильтониан аъ (x, р, t) имел удобный для нас вид.

3.4. Уравнения третьего и четвертого порядка для усредняющей замены переменных

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

Х°п (X,Р, 0 = д« (х,р, ^/дР„ +

+ (12)Е(Р/ дХа/дР„ -Х, дР//дР„) +

к

+ (12)К(Р/ дХЦдРп -ХЬ дР//дР„) +

к

+ (Da|2)Х(Рка дХка/дР„ -Ха дР//дР„) -

- DaX^l -DьXatt - Dcx„/2,

Рт (х,Р, 0 = -д«с (х,р, 0/дхт +

+ (12)Е(Ха дР/дхт -РЬ дХ1 /дхт) +

к

+ (12) Е( ХЬ дР//дхт - Рка дХЦ дхт) +

к

+ (Da|2)К(Хка дР//дхт -Р,а дХка/дхт) -

к

- D Рь - DhPa - D р /2,

ат Ьт с г т1 '

где « (х,р, ^ — произвольная вспомогательная функция.

Поправка Gc (х,р, ^ к гамильтониану определяется по формуле

Gc (х,р, 0 = DcH0 (х,р, ^ + DъHa (х,р, ^ + + DaHь (х,р, ^ + Нс (х,р, ^ - д«с (х,р, ^/дt +

+(12)Е(Ха дPkVдt - рь дХка/д^ +

к

+ (12)Е(ХЬ дPka|дt-Рка дХЦд) +

к

+ (Da|2)Х(Хка дР//дt -Рка дХка!дО .

к

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

Хап (х,р, 0 = два (х,р,^/ф„ +

+ (12)КРкс дХк/дР„ - Ха дРкс/дР„) +

к

+ (12)К(РкЬ дХЬ/дР„ -ХЬ дР//дР„) +

к

+ (12)КРка дХск/дР„ -Хс дР//дР„) +

к

+(Db|2)К(Р/ дХаа/дР„ - Ха дРка/дР„) +

к

+ (Da|2)К(РЬ дХка/дР„ -Хка дР//дР„) +

к

+ (Da|2)К(Рка дХЬ/дР„ - ХЬ дРка/дР„) -

- DaXn -DьX„Ь - ^Х„ь- Ddxn/2,

Рт (х,р, 0 = -д^ (х,р, 0/дхт +

+ (12)КХка дР//дхт - Ркс дХ,/дхт) +

к

+ (12)E(X¿ дР,*/^ - Pb X/dXm) +

к

+ (12)КXе дР//дт» -Р/ дХк/дх») +

к

+ (Db/2)К(Xa дР//дх» - Ра дХ1/дх») +

к

+ (Da/2)К(ха дР//дх» - РЬ дХк/дКт) +

к

+ ( DJ 2) К( Xb дР//дх» - Ра дХЦ дх») -

к

- D Рс - D^ - D Рь - D, p /2 ,

a rn b » с rn d г»! '

а поправка к гамильтониану вычисляется, как G¿ (x,p, í) = DdHo (x,p, í) + DcHa (x,p, í) + + DbHb (x,p, í) + DaHc (x,p, í) + Hd (x,p, í) -

- д^ (x,p,í)/dí +

+ (12) K(Xa дР// dí - Р/ dX^/dí) +

к

+ (12)K(Xb dFkVSí - рь dXjdí) +

к

+ (12)К(хс dP//dí -Pka dXl/dí) +

к

+(Db/2) К Xka dP//dí - pka ¿xa/Sí) +

к

+ (Da/2)К(Xka dP//dí -Pb dXkJdt) +

к

+ (DJ2)К(Xb дР//dí -Р/ дХЦdí) ,

к

где «S, (x,p, í) — произвольная вспомогательная функция.

4. МЕДЛЕННОМЕНЯЮЩИЙСЯ ЭФФЕКТИВНЫЙ ПОТЕНЦИАЛ

Наш исходный гамильтониан имеет вид

H = — (Р2 + Q2 + R2) + eU + 2»

+s • е^( V(k) cos wkí + Ж(k) sin wkí) ,

где X, 7, Z — пространственные координаты; P, Q,R — соответствующие им импульсы; í — время; U = U (X, Y, Z, í) — низкочастотный ("мед-ленноменяющийся") потенциал электрического поля; V(k) = V(k) (X,Y,Z,í) и Ж® = W(k)(X,Y,Z,í) —

"медленноменяющиеся" потенциалы, соответствующие гармоникам высокочастотного электрического поля; cok — (круговые) частоты соответствующих гармоник; rn — масса частицы; е — заряд частицы, а s — малый параметр. Сделаем для этого гамильтониана замену переменных:

X = х + ^.Х + s2 X, +..., Y = y + sY + s2Yb +...,

a b ' s a b '

Z = z + sZ + s2Zh + ..., P = p + sP + s2Ph + ... ,

a h ' .i a h '

Q = q + sQa + s2Qh + ... , Я = r + si?a + + ...,

где х, у, г — усредненные координаты Х, Y, Z, Р, q, г — усредненные импульсы Р, 2, ^ , а функции Ха (. ) , Хъ (. ) и т. д. подлежат определению как функции независимых переменных х, у, г , Р, q, г и t.

После подстановки новых (усредненных) переменных в старый гамильтониан он может быть представлен в виде ряда по малому параметру е , как Н = Н0 + еНа + е2 Нъ +..., где функции Н0 (...), На (...), Нъ (...), ... являются в конечном счете функциями новых переменных х, у, г , р, q, г и t. Выпишем в явном виде эти функции, поскольку они нам понадобятся в дальнейшем:

Н0 = 2т(р2 + + Г2) + (хУ*t) •

н. =-(pP. + qQ. + R) +

m

+ e (üxXa + U/„ + UzZ.) + + e^ ( V(k} cos cokt + Ж(k > sin cokt),

Hb = (P2 + Q2 + R + 2pPb + 2qQb + 2rR) + 2m

+ e (UxXb + UyX + UzZb) + +2 e (UxX2 + U^Y.2 + UZ +

+2их,ХаГа + 2UxzXaZa + ) +

+ e£[(v¿k + V(k Y. + V(k > Z.) cose + +(^(k+ W(kX + Ж®Z.)sinckt] .

Новый гамильтониан G разлагается в ряд по малому параметру s и может быть представлен, как

G(•••,t) = G0(•••,t) + sG.(•••,t) + s2Gb(•••,t) +...,

где функции G0 (...) , G. (...) , Gb (...) как функции новых переменных x, y, z , p, q, r и t нам, собственно, и надо определить. Из раздела 3.1 следует, что для элемента G0 (• • •, t) справедливо выражение

Go = Hо = (Р2 + q2 + r2) + eU(x,y,z,t) .

Нахождение остальных членов ряда для нового гамильтониана потребует существенно большей работы.

Из раздела 3.2 следует, что для первого порядка замена переменных Xa (•••,t), Ya (•••,t), Za (•••,t),

Pa (• • •, t) , Qa (• • •, t) , Ra (• • •, t) может быть записана с помощью вспомогательной функции Sa ( x, y, z, p, q, r, t) и константы Da, как

Xa =dSjdp - DaX/2, Ya =dSj dq - Day/2, Za = dSj5r - DaZ/2, Pa = - dSjdx - Dapl2, Qa =-CSjdy - Daq/2, Ra = -dSjdz - Dar/2.

При такой параметризации замены переменных гамильтониан равен

Ga (•, t ) = DaH 0 + Ha - dSj dt .

Из приведенного выше представления для функции Ha (•••,t) легко заметить, что для того чтобы

"убрать" из нового гамильтониана Ga (•••,t) все

быстроосциллирующие члены, нам надо использовать функции вида

Sa (-.., t ) =

= Sao) +Z( S^cos^t + S^sin^t),

где функции S^ (..., t), S« (..., t) , S« (• • •, t) зависят от координат x, y, z , импульсов p, q, r и "медленного" времени t .

Подставив это представление для функции Sa в формулы для вычисления функций Xa , Ya , Za , Pa, Qa, Ra, мы легко получим, что эти функции будут содержать только быстроосциллирующие члены, если выполнены соотношения

dSao)/dp = Dax/2, dSaoV dq = Day/2,

dSaoVdr = Daz/2, dSf/dx = -Dap/2, dSaoVdy = -Daql2 , dSao7dz = -Dar/2 .

Из условий совместности для этой системы уравнений в частных производных следует, что это возможно только при Da = 0 . В этом случае решением указанной системы уравнений является функция S^ = 0 (с точностью до аддитивной функции времени, которая на уравнения движения, выводимые из гамильтониана, никак не влияет).

Для того, чтобы в гамильтониане Ga (•••,t) отсутствовали "быстрые" функции, необходимо и достаточно, чтобы были выполнены уравнения

dS«/dt +

+ (p • dSVkay dx + q • dS^/dy + r • dS^ /dz)/m -

-dS«/ dp - ^ - SS«/ dq - dr

+

= eV(к),

дБ«/ дЛ +

+ (р • дБ« /дх + q • дБ« /ду + г • дБ« /дг)/га --еЦ. • дБ« /др - еиу • дБ« /дq - е^ • дБ« /дг -ю Б« = ^). ' '

Любое решение этих уравнений обеспечивает нам замену переменных с требуемыми свойствами, и, как легко понять, таких решений будет достаточно много. При этом сам гамильтониан Ga остается неизменным (Ga = 0), однако уравнения для компонентов гамильтониана более высоких порядков зависят от используемой на этом уровне замены переменных самым существенным образом.

Для нашей цели далеко не любое решение указанной системы дифференциальных уравнений в частных производных является допустимым и осмысленным. Из физических соображений следует, что полученные формулы должны оставаться регулярными при больших значениях частот высокочастотного электрического поля. В частности, следует ожидать, что при юк ^ да замена переменных, использованная для перехода от физических переменных к усредненным переменным, должна обращаться в тождество и что указанная замена переменных является регулярной в окрестности точки ак «да по меньшей мере. С учетом этого соображения удобно искать решение в виде ряда по обратным степеням круговых частот 1/юк . При таком представлении подходящие решения указанных дифференциальных соотношений оказываются строго единственными, причем для наших непосредственных целей достаточно ограничиться только главными членами разложения. Опуская ненужные технические подробности, получаем с точностью до квадратичных членов требуемый результат:

Бк «-е^(кV+

+ е (г/к) + ( рУ(к) + qV(k) + гУ(к)) /га) /юI Б« «+eV(k7юк +

+.

+ е (Щ(к) + (рЖ(к) + qW(k) + гЖ}к)) /га) /ю

+..

Из раздела 3.3 следует, что во втором порядке замена переменных может быть записана с помощью вспомогательной функции Бь (х, у, г, р, q, г, л) и константы Dь, как

Хь =дБь/др + Е(р) -DъXa -Dbx|2, Г =дБь/дq + Е« -DьYa -Dьу|2,

2ь = дБь/дг + Е(г) -DbZa -Dbz|2, Рь =-дБь/дх + Е(х) -DbPa -Dьp|2, Qь =-дБь/ду + Е(у) -DьQa -Dbq|2, Rь =-дБь/дг + Е(г) -DbRa -Dbr|2,

где

Е(р) = (12)(Ра дХа/др - Xa дРа/др + Qa д¥а/др -

-Уа дQa/др + Ra дZa|др - Za д^ /др) , Е(<г) = (12) (Ра дХа/дq - Xa дРа/дq + Qa дГа/дq -

-Уа дQa/дq + Ra дZa|дq - Za д^/дq), Е(Г) = (12) (Ра дХ а/ дг - Ха дР„ /дг + Qa дТ, / дг -

-Га дQa/дг + Ra дZ7 дг - Za дR7 дг) , Е(х) = (12)(Ха дРа/дх - Ра дХа/дх + Г дQa^дх -

^а дТа1 дх + Za дR7дх - Ra дZa /дх), Е(у) = (12)(Ха дРа/ду -Ра дХа/ду + Га дQa^ду -

^а дТа1 ду + Za дR7ду - Ra дZа /ду) , Е(г) = (12)(Ха дРа/дг -Ра дХа/дг + Га дQa/дг -- Qa дГа/дг + Zа дЯа1 дг - Ra дZa|дг),

а функции Ха, Га, Za, Ра, Qa, Ra соответствуют замене переменных, полученных при анализе уравнений для первого порядка. В случае такой параметризации для замены переменных новый гамильтониан равен

Оь (•••,Л) = DьЯ0 + DaHa + Нь + Е(л) -дБ^дЛ,

где

Е(Л) = (12) (Ха дРа/ дЛ - Ра дХа/ дЛ + Га дQa^ дЛ -

^а дГа/дЛ + Za дК^ дЛ - Ra дZa / дЛ) .

Поэтому для того, чтобы "убрать" из нового гамильтониана Оь (• • •, л) все быстроосциллирую-щие члены, нам надо использовать функции вида

Б (., Л) = ^ + Х( Б^есзЮкЛ +

к

+ Б« 81П юкЛ + Б2^ь 0С8 2ЮкЛ + б2^),ъ 8ш2ю^) +

+ ) С08 (юк + ю})Л + Б+^т (юк + ю,)^ +

к > ,

+ Б^ь)С08(юк -ю,)Л + Б-^т(юк -ю,)л] ,

где неизвестные функции зависят от координат х, у, г , импульсов р, q, г и "медленного" времени Л. Существенно, что в соответствии со сде-

ланным ранее предположением о том, что частоты (ок должны быть "далеко разнесенными", функции

cos(¡k -а>j) и sin(¡k -() будут "быстрыми".

Соответствующие выкладки требуют привлечения серьезных программ для символьных вычислений типа [21]. Легко получить, что S^'b = 0 и SWkb = 0 . Учитывая, что решения системы дифференциальных уравнений в частных производных для функций S%b, S2t)b, s+k/b), s+kj, s-kv',b), s-kj должны быть регулярными при сок,( ^<х>, можно вывести, что

- -e2 (V(k)W(k) + V(k)W(k) + V(kWk) )/(4mrn¡),

^ ^ («k»)!+(vk»)!»)2 -

-(Wf > )2-(W« > )2-(W,<k' )2),

S(k,j) - e2 (v(k)W(j) + V(k)W(j) + V(k)W(j) + W(k)V(j) + +Wy(k)V;j) + Wk)V0) )•(-¡j )7 (4m©k©2 ( +( ))

S(k^-j) - e2 (—V(k)V(j) -V(k)V(j) -V(k)V(j) + W(k)W(j)

+

д^/др = Пьх/2, дSЪ0)/дq = Пъу/2, д^/ дг = Вьг/2, дSЪ07 дх = - Пър/2, д^07 дУ = - 2, д^07 дг = - Пьг12,

откуда следуют условия Пъ = 0 и SЪ0) (•••,t) = 0. Из всех этих соотношений можно получить, что

G,

+Ж(к)цу(/) + цХЮцХЛ)•© -©) /(4т©©2(©к +©/)), S(k,j) «-е2 (у(к)у( /) + .) + У^У.) + Ц(к)Ц(.) + +цук)цС/) + ).© +©)У(4т©2©2(©к )),

S(к ^) да -е2 (ц(к )у (/) + ц(к )у (/) + ц(к (/) - у (к )Ц( /) -

-У^цШ -у(к)^^/)).(©к +©/)У(4т©2©2(©к -©/)).

Здесь сохранены члены, не более чем кубичные по степеням 1/©к .

Уравнения в частных производных для функции SЪ0) (•••,t), выражающие тот факт, что функции Хъ (•••,0, Yъ (•,0, ^ (•••,0, Ръ (•,t) , 2Ъ (•, t) , Rъ (•, t) содержат только быстроосцил-лирующие члены с нулевым средним значением, являются достаточно сложными, и, вообще говоря, для них не выполняются условия совместности в случае произвольных высокочастотных электрических полей. Однако, если в этих уравнениях ограничиться членами, не более чем квадратичными относительно обратных частот 1/©к , они принимают уже знакомый нам вид:

(•, t)« е2 Е^((Ух( к) )2 +(Уу(к) )2 +(Уг(к) )2 +

+ (Жк) )2 + (ж(к) )2 + (Ц(к) )2)/( 4т©2)

Здесь сохранены члены, не более чем квадратичные по степеням 1/©к .

5. ОКОНЧАТЕЛЬНЫЙ РЕЗУЛЬТАТ

Мы получили, что после введенной нами замены переменных исходная динамическая система с гамильтонианом Н (•, t) переходит в новую динамическую систему с гамильтонианом

а (•, t)« о0 (•, t)+еаа (•, t)+е2а (•, t)+•,

где функции 00 (•, t) , Оа (•, t), Оъ (•, t) вычисляются, как

°0 (•,t) = 2т(Р2 + г2 + ) + еи0 (х,У,г,t),

аа (•, t ) = 0,

аъ (•, t) - е2 Е^((Ух( к) )2 +(Уу(к) )2 +(Уг(к) )2 +

+ (Жк) )2 + (Ж(к) )2 + (Ц(к) )2)/( 4т©2)

Замена переменных целенаправленно организована так, чтобы новые переменные имели физический смысл исходных физических переменных, усредненных по быстрым осцилляциям. Величины

|Екс) (х, у, г, t)|2 = (дУ(к) (х, у, г, t)/дх)2 +

+

(dV(k) ( х, y, z, t )/cy ) + (CV(k) ( х, y, z, t )/dz )

и

|Eks) (х,y, z, t)|2 = (aW(k) (х, y, z, t)/cX)2 +

(a w (k) (х, y, z, t )/cy )2 + (a w (k) (х, y, z, t )/az)

+

имеют смысл квадратов модулей напряженности электрического поля (точнее, квадратов амплитуд напряженности высокочастотного электрического

2

2

поля). Функции U0 (x,y, z, t) , V(k) (x,y,z, t) и W(k) (x, y, z, t), использованные в этих формулах, берутся из представления потенциала U (x, y, z, t) высокочастотного электрического поля в виде разложения

U ( x, y, z, t) = U0 ( x, y, z, t) +

+£X(V(k) ( x, y, z, t) cos a>kt + W(k) ( x, y, z, t) sin a>kt)

по быстрым и далеко разнесенным вдоль оси частот гармоническим осцилляциям cos cokt и sin a>kt.

Тем самым решения динамических уравнений, полученных для усредненных переменных (координат и импульсов), будут такими, как будто заряженная частица с зарядом e и массой m двигается в электрическом поле с потенциалом U0 ( x, y, z, t) + s2U ( x, y, z, t), где U0 ( x, y, z, t) — это низкочастотная (квазистатическая) составляющая электрического поля, а эффективный потенциал U ( x, y, z, t) определяется выражением

U ( x, y, z, t) =

= el(| Ekc) ( x, y, z, t) |2 + E) ( x, y, z, t) |2) /( 4m®k2)

и также является медленной (низкочастотной и квазистатической) функцией времени. При этом указанный эффективный потенциал U (x, y, z, t)

уже не обязан удовлетворять уравнению Лапласа, что обеспечивает нам новые интересные возможности по управлению заряженными частицами.

Отметим, что "малость" малого параметра s при выводе соответствующих формул требует всего лишь, чтобы соответствующий формальный ряд был сходящимся, или чтобы по крайней мере следующий (опущенный) член разложения оказывался много меньше последнего сохраненного члена ряда — и, по-видимому, не является слишком обременительным условием. Более существенно, что при выводе этих формул использовались только главные члены разложения по степеням обратных частот 1/ сок , справедливых при условии сок ^ да .

Тот факт, что поправка Gb (• • •, t) зависит только

от усредненных координат и не зависит от усредненных импульсов (и, следовательно, может трактоваться как поправка к электрическому потенциалу), является следствием сохранения исключительно главных членов разложения и отбрасывания последующих членов ряда по степеням 1/сок . Сохранение отброшенных младших членов и переход тем самым от эффективного потенциала к эффективному гамильтониану может быть интересным направлением в плане уточнения и обоб-

щения классической теории эффективного потенциала.

Подведем промежуточные итоги. По определению, эффективный потенциал U (x, y, z, t) — это такая скалярная функция, высчитываемая по определенным правилам через имеющееся в системе высокочастотное поле, что усредненное движение заряженной частицы в данном высокочастотном поле описывается уравнениями движения заряженной частицы в псевдоэлектрическом поле U (x, y, z, t) с точностью до поправочных членов малого порядка (где вопрос о том, что именно считать членами малого порядка, является достаточно нечетко определенным и в значительной степени определяется априорными предположениями и вкусами исследователя). Математическое выражение для эффективного потенциала в случае использования класса высокочастотных электрических полей вида

E0 ( x, y, z, t) +

+ZE(kc) (x,y, z, t)cos (akt) + Eks) (x,y, z, t) sin (^t),

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

ЗАКЛЮЧЕНИЕ

В следующих работах данного цикла будет по-казанана универсальность представления временных сигналов, характеризуемых "медленными" и "быстрыми" временами, в форме

f (t) = Z An (t) cos КО + B (t) sin (ffl„t), где An (t) и Bn (t) — "медленные" функции, а a>n — "быстрые" и "далеко отстоящие" друг от друга частоты (и, как следствие, обоснована возможность представлений высокочастотных электрических полей, характеризуемых "медленными" и "быстрыми" временами, в виде

E ( x, y, z, t) = E0 ( x, y, z, t) +

+Z Enc) (x, y, z, t) cos Knt) + E(ns) (x, y, z, t) sin (^t),

n

где E0 (x, y, z, t), E(nc) (x, y, z, t), E(ns) (x, y, z, t) —

"медленные" функции, а соп — "быстрые" и "далеко разнесенные" друг от друга частоты). Будут рассмотрены технические способы создания требуемых электрических напряжений fk (t) и соответствующих им высокочастотных электрических полей в объеме устройства. Приведены примеры использования высокочастотных электрических полей указанного вида для управления движением заряженных частиц. Наконец, будут рассмотрены перспективные способы применения нового метода управления движением заряженных частиц для создания новых классов масс-спектрометрических приборов.

СПИСОК ЛИТЕРАТУРЫ

1. March R.E., Todd J.F.J. Quadrupole ion trap mass spectrometry, 2nd edition. Wiley-InterScience, 2005. 346 p.

2. Major F.J., Gheorghe V.N., Werth G. Charged particle traps. Springer, 2005. 354 p.

3. Werth G., Gheorghe V.N., Major F.J. Charged particle traps II. Springer, 2009. 275 p.

4. Ландау Л.Д., Лифшиц Е.М. Механика, сер. Теоретическая физика. М.: Физматлит, 2004. 220 c.

5. Заславский Г.М., Сагдеев Р.З. Введение в нелинейную физику: от маятника до турбулентности и хаоса. М.: Наука, 1988. 368 c.

6. Yavor M.I. Optics of charged particle analyzers, Ser. Advances of imaging and electron physics, V. 157. Elsevier, 2009.

7. Чирков А.Г. Асимптотическая теория взаимодействия заряженных частиц и квантовых систем с внешними электромагнитными полями. СПб.: Санкт-Петербургский государственный технический университет, 2001.

8. Гапонов В.А., Миллер М.А. О потенциальных ямах для заряженных частиц в высокочастотном поле // Журнал экспериментальной и технической физики. 1958. Т. 34, № 2. С. 242-243.

9. Миллер М.А. Движение заряженных частиц в высокочастотных электромагнитных полях // Известия вузов, сер. Радиофизика. 1958. Т. 1, № 3. С. 110123.

10. Бурштейн Э.Л., Соловьев Л.С. Гамильтониан усредненного движения // Доклады АН СССР. 1961. Т. 139, № 4. С. 855-858.

11. Литвак А.Г., Миллер, М.А., Шолохов Н.В. Уточнение усредненного уравнения движения заряженных частиц в поле стоячей электромагнитной волны // Известия вузов, сер. Радиофизика. 1962. Т. 5, № 6. С. 1160-1174.

12. Сивухин Д.В. Дрейфовая теория движения заряженной частицы в электромагнитных полях // Вопросы теории плазмы, вып. 1. М.: Госатомиздат, 1963. С. 7-97.

13. Морозов А.И., Соловьев Л.С. Движение заряженной частицы в электромагнитных полях // Вопросы теории плазмы, вып. 2. М.: Госатомиздат, 1963. С. 177-261.

14. Митропольский Ю.А. Метод усреднения в нелинейной механике. К.: Наукова Думка, 1971. 440 с.

15. Боголюбов Н.Н., Митропольский Ю.А. Асимптотические методы в теории нелинейных колебаний. М.: Наука, 1974. 503 с.

16. Моисеев Н.Н. Асимптотические методы нелинейной механики. М.: Наука, 1981. 400 с.

17. Гантмахер Ф.Р. Лекции по аналитической динамике. М.: Физматлит, 2002. 386 с.

18. Трикоми Ф. Лекции по уравнениям в частных производных. М.: ИЛ., 1957. 443 с.

19. Гюнтер Н.М. Интегрирование уравнений первого порядка в частных производных. М.-Л.: ОНТИ, 1934. 359 с.

20. Фиников С.П. Метод внешних форм Картана в дифференциальной геометрии. М.-Л.: ОГИЗ, 1948. 432 с.

21. Программа для символьных вычислений Mathema-йса. URL: (http://www.wolfram.com).

Институт аналитического приборостроения РАН, г. Санкт-Петербург

Контакты: Бердников Александр Сергеевич, [email protected]

Материал поступил в редакцию 20.04.2011.

TIME-DEPENDENT PSEUDOPOTENTIAL AND ITS APPLICATION FOR DESCRIPTION OF THE AVERAGED MOTION OF THE CHARGED PARTICLES.

PART 2. GENERAL EXPRESSION FOR TIME-DEPENDENT

PSEUDOPOTENTIALS

A. S. Berdnikov

Institute for Analytical Instrumentation of RAS, Saint Petersburg

The series of publications describes a new method to control the movement of the charged particles by high frequency electric fields. This new class of high frequency electric fields is characterized by pseudopotentials slowly evolving in time. The paper considers the deriving of the basic expression for slowly evolving pseudopotentials. It is shown that the classical pseudopotential is just the leading term of decomposition of the true pseu-do-Hamiltonian in a double series: with respect to a small parameter which characterizes the amplitude of the electric field, and with respect to the frequency of the electric field. The classical pseudopotential expression is generalized to include the case of high frequency electric fields which are characterized by "fast" and a "slow" characteristic time scales simultaneously.

Keywords: mass spectrometric devices, radio frequency devices, pseudopotentials, Hamiltonian dynamics, Hamiltonian dynamic systems, averaging techniques for Hamiltonian dynamics, high frequency electric fields, radio frequency electric fields