CC BY
11ISSN 0868-5886 НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2019, том 29, № 1, c. 143-148
МАТЕМАТИЧЕСКИЕ МЕТОДЫ И МОДЕЛИРОВАНИЕ В ПРИБОРОСТРОЕНИИ
УДК 537.534.3:621.384.8 (075.8) © К. В. Соловьев, М. В. Виноградова, 2019
УТОЧНЕНИЕ УСЛОВИЙ УСТОЙЧИВОСТИ ДВИЖЕНИЯ
ИОНОВ В ЭЛЕКТРОСТАТИЧЕСКОЙ ЛОВУШКЕ, ИНТЕГРИРУЕМОЙ В ЭЛЛИПТИЧЕСКИХ КООРДИНАТАХ
Для электростатических полей идеально фокусирующих ловушек, реализующих разделяемое в эллиптической системе координат движение заряженной частицы, уточнены условия, при которых обеспечивается финитное по всем координатам движение частицы. Для нулевых и близких к ним значений параметра поля ^ найдены аналитические точные и приближенные соотношения, связывающие начальные данные и параметры поля, гарантирующие финитность движения.
Кл. сл.: идеальная пространственно-временная фокусировка, масс-спектрометрия, ионная ловушка
ВВЕДЕНИЕ
ИССЛЕДУЕМАЯ СИСТЕМА
Впервые предложенный Ю.К. Голиковым и в 1985 г. запатентованный [1] принцип идеальной пространственно-временной фокусировки (ИПВФ) ионов в квадратичном по одной из координат электростатическом поле является в настоящее время теоретической основой синтеза обширного класса перспективных масс-спектрометрических приборов [2]. Интерес к полям данного типа не ослабевает, в том числе благодаря появившимся в последнее время новым технологическим возможностям, делающим потенциально доступным изготовление с хорошей точностью систем с криволинейными трехмерными электродами [3].
Среди множества реализующих принцип ИПВФ электростатических потенциалов вида
(( x, y, z) = z2 + f (x, y)
(1)
Потенциал исследуемой системы в декартовых координатах имеет вид [6]:
2 2
x, ч 2 x2 + y ф( x, y, z) = z---+
2
+ Re
í
=( мы (í )+£+ 1S )
где
g = x +i y,
, (2)
(3)
выделяется класс полей, обеспечивающих движение, интегрируемое методом разделения переменных в уравнении Гамильтона—Якоби [4]. Свойства ионных ловушек, создаваемых на базе таких полей, часто могут быть исследованы аналитически, что позволяет получить прозрачные связи между параметрами задачи.
Настоящая статья продолжает изучение интегрируемых электростатических ловушек, предпринятое ранее в работах [4-8]. В данной работе делается попытка уточнить условия устойчивости движения заряженной частицы в интегрируемой ловушке с разделением переменных в эллиптических координатах, ранее рассматривавшейся в [6, 7].
Re — реальная часть комплексной величины, ц , s и S — вещественные константы.
В эллиптических координатах а = u + i v, связанных с декартовыми (3) соотношением д = cha , или x = ch u cos v , y = sh u sin v, потенциал (2) приобретает вид
ф(и, v, z) = z2 -
sh2 u + cos2 v
2
+
1 (pu + e)sh2u + (pv+5)sin2v ,.)
+ _ i 2 -2 ' ( )
2 sh u + sin v
Движение иона в поле (4) — синусоидальное по координате z, по u и v оно определяется уравнениями
2 2 + U (u) = C, + V (v) = -C, (5)
где pu, pv — и, v импульсы, а U(и), V(v) — соответствующие эффективные потенциалы:
U(и) = -—sh4 и - I — + E |sh2 и + 2 12 )
+ —(ци +£^2и,
V(v) = —sin4 v - I — + E |sin2 v + 2 12 )
+ —(uv+5)sin2v. 2
(6)
(7)
Величины С, Е в (5)-(7) — константы разделения переменных в уравнении Гамильтона—Якоби, связанные с начальными условиями движения частицы.
Возможность удержания иона в объеме ловушки, гарантированная по г квадратичностью соответствующей составляющей потенциала, определяется фактом одновременного наличия ям потенциалов (6), (7). При этом дно по крайней мере одной из этих ям должно находиться в области отрицательных значений, а берега другой обеспечивать удержание иона в области значений положительных. Эти требования вытекают из специфики уравнений (5), в правые части которых входит с противоположным знаком одна и та же константа С.
В работе [6] третьи слагаемые эффективных потенциалов (6), (7) рассмотрены как возмущения структуры ям, определяемой первыми двумя слагаемыми (6), (7). Этот подход позволил подтвердить удерживающие способности обеих ям, выявив условие Е <-1/2 в качестве необходимого для наличия ямы по и. В данной работе мы хотим продвинуться дальше, уточнив условия финитно-сти движения заряженных частиц в поле ловушки для случая 8 Ф 0 .
УТОЧНЕНИЕ УСЛОВИЙ УСТОЙЧИВОСТИ
Потенциалы (6), (7) существенно упрощаются при / = 0, вновь открывая возможности для использования аналитических подходов.
Начнем с рассмотрения движения по V. Преобразуем (7) с учетом /и = 0:
V (v) = —cos 4v + E cos 2v + 16 2
5 -o —(r — + — sin2v—I E + —
2 21 8
личие ямы при любых значениях E, 5 . Очевидно, при 5 = 0 значение V(0) = 0 для любых E, причем для E <-1/2 это значение является минимумом V(v) — дном v-ямы (поскольку V"(0) = -2E -1; при E > -1/ 2 это уже не так, имеем локальный максимум). Наибольшее значение Vup = V/2) = -E
(для 5 = 0, E < 0) позволяет определить глубину ямы при 5 = 0. В случае 5 Ф 0 определение параметров ямы требует решения уравнения
V '(v) = -1 sin4v - E sin 2v + 5 cos2v = 0, (9)
которое подстановкой q = exp 2iv (в предположении q Ф 0) приводится к разрешимому в радикалах алгебраическому уравнению четвертой степени с комплексными коэффициентами
q4 + 4(E - i 5) q3 - 4(E + i5) q -1 = 0.
(10)
Решая (10) для заданных Е, 8 и выделяя корни qi, для которых V'(— 1п qi /2) = 0, получаем возможность построить поверхности Vшf(E,8), К,ир(Е,8), определяющие глубину v-ямы и диапазон возможных значений эффективного потенциала V (рис. 1).
-1.0
Периодичность функции V(v) гарантирует на-
(8) Рис. 1. Параметры УМ(Е,8) (дно), 1 [щ1(/<.д) (тах уровень удержания в финитном состоянии) потенциальной ямы
Рассмотрим далее функцию (6), определяющую движение иона по и. Поскольку и (и) ^-<х> при и ^ +<» и и(0) = 0, (6) имеет по крайней мере один экстремум (максимум, если экстремум единственный). Необходимые для получения ямы три экстремума (два максимума и минимум между ними) при / = 0, £ = 0 обеспечиваются положительностью коэффициента при sh2 и в (6), т.е. выполнением условия Е < -1/ 2 [6].
С учетом / = 0 удобно переписать (6) в виде
1 E
U (и) =--ch4u--ch 2u +
16 2
+ — sh 2u + — | E + -|.
2 21 8 J
(11)
Положение экстремумов определяется нулями производной
U'(u) = - 1sh4u -Esh2u + sch2u . (12)
При s = 0 для E >-1/2 имеем только максимум Umax = U(0) = 0, для E < -1 / 2 — минимум Umin =
= U (0) = 0
и
два
максимума
U_ =
= и1агсМ-2£)j = 1 ^Е2 + Е +1 j (рис. 2).
Исследование случая £ Ф 0 требует решения уравнения и '(и) = 0, трансформирующегося в результате подстановки q = ехр2и в (12) (с учетом q Ф 0) в уравнение
д4 + 4(Е-£)q3 -4(Е + £)q-1 = 0. (13)
Рис. 2. Эволюция конфигурации потенциальной ямы (6) для случая /и = 0
В отличие от периодической функции (8), здесь наличие ямы при £ Ф 0 для Е <-1/2 не гарантировано. Яма будет существовать для тех комбинаций Е, £ , при которых (13) обладает по крайней мере тремя вещественными положительными корнями, обеспечивая три вещественных решения уравнения и '(и) = 0. Поскольку произведение всех корней уравнения (13) равно минус единице, этот вариант возможен, (13) должно иметь три положительных и один отрицательный корень. Область финитности движения по и в координатах (Е, £ ) может быть построена селекцией "правильных" вариантов решений (13). По факту эта область может быть задана неравенством — условием наличия только вещественных корней у резольвенты (13) [9, (1.8-15)] — кубического уравнения у3 +16(£2 - Е2 +1/4)у - 64Е£ = 0. Данное неравенство имеет вид в [9] — условие Q < 0 , где Q вычисляется согласно [9, п. 1.8-3, (1.8-7)]):
4—2 -E2 +1/4)3 + 27(Es)2 <0,
(14)
выделяя совместно с неравенством E < -1/ 2 область допустимых значений E, s (рис. 3).
Вычисление значений эффективного потенциала U в точках экстремумов позволяет построить
граниЦЫ изменения U: Umin <U < Uminmax =
= min(Umax1,Umax2), в рамках которых реализуется финитное u-движение (рис. 4).
Далее рассмотрим случай малых p Ф 0. Оставаясь в рамках аналитического подхода, можем использовать лишь аппроксимации u « sh 2u / 2, v « sin 2v / 2 , причем приближения вводим непосредственно в выражения для производных (что
Рис. 3. Область значений (£, Е), обеспечивающих устойчивость движения иона по координате и для значений /и = 0 (точно) и / = ±0.1 (приближенно, согласно (16))
1" И
Е +i 8 - и / 2 , „ - 4---q -1 = 0,
1 - И
а (13) — уравнением
(15)
4 лE-е-и/2 з
q + 4-----q3 -
1 - И
„Е+е-и/2 , „
-4--— q -1 = 0.
1 - И
(16)
рационныи процесс
U = 1 Arch (-2Е + и (2ui-1 cth 2ui1 +1)), i = 1, 2,...; u0 = ^Arch(-2E),
Рис. 4. Параметры Umin(E,e) (дно),
U_(E,e) = min(Umx1(E,е), Ц^ВД) (max уровень удержания в финитном состоянии) потенциальной ямы (6) для случая и = 0
дает более точный результат, чем аппроксимация и и V в исходных потенциалах (6) и (7)). В этом случае (10) заменяем уравнением
4 Е -15 - и/ 2 3 q + 4---q -
Приближения достаточно грубы, однако применение большего количества членов разложения приводит к уравнениям порядка выше четвертого, решение которых в радикалах становится невозможным. Впрочем, до значений примерно | и | < 1/4 результат данной аппроксимации вполне может служить в качестве оценки либо точки старта уточняющего итерационного процесса.
Для меньших единицы значений и при е = 0 можно приближенно определить положение максимумов (6). Используя уравнение и'(и)|Е=0 = 0 или эквивалентное ему при условии и Ф 0 уравнение Л 2и = -2Е + и (2и сЬ 2и +1), организуем ите-
полагая, что под знаком арккосинуса стоит не меньшая единицы величина. Уже формула для и2 дает приемлемое приближение для правого максимума, здесь не приводим ее ввиду громоздкости.
Поскольку условия финитности по и и V связаны общей константой разделения С, входящей в (5), есть возможность построить проекцию зоны устойчивости в пространство ( е, 5, С) при фиксированном значении Е. Действительно, полагая С > ^ из (5) имеем ит1П <С < Vnf <
< -С < Vsup или -< С < -Vinf. Следовательно, для каждого е, 5 можно получить тах(ит1П, - ^р) <
< С < т1п(итахт1„, - определяя тем самым границы изменения С. При ¡и = 0, е = 0, 5 = 0 дно обеих ям расположено на нулевом уровне, следовательно, С = 0 и возможны лишь колебания ионов на оси г с нулевым фазовым объемом. Неравенство нулю любого из параметров ¡и, е, 5 устраняет это ограничение.
Заметим, что условие и Ф 0 нарушает периодичность (7) и снимает ограничение по V на допустимые значения С ввиду линейного (с коэффициентом и / 2) нарастания амплитуды члена 8т 2у в (7). Для движения по и условие и Ф 0 приводит к сдвигу области устойчивости по Е, при е = 0 критерий существования и-ямы меняется на Е < и-1/2 (рис. 3).
ЗАКЛЮЧЕНИЕ
Таким образом, в работе определены новые и уточнены найденные ранее в [6] условия финит-ности движения ионов в интегрируемой эллиптической электростатической ловушке. Для полей с нулевым и близким к нулю значением параметра и в (2), (4) удалось получить результаты в аналитическом виде, дающем ясные взаимосвязи между характеристиками поля и допустимыми начальными данными. Применение новых подстановок в ходе решения задачи позволило упростить алгоритм ее решения. Результаты статьи могут быть использованы при синтезе новых масс-спектро-метрических систем на базе электростатических ловушек с идеальной пространственно-временной фокусировкой.
Данная работа частично выполнена в рамках государственного задания № 075-00780-19-00 для Института аналитического приборостроения РАН.
СПИСОК ЛИТЕРАТУРЫ
1. Галль Л.Н., Голиков Ю.К., Александров М.Л. и др. Времяпролетный масс-спектрометр. Авторское свидетельство СССР № 1247973 от 01.04.1986. Приоритет от 16.01.1985.
2. Hu Q., Noll R., Li H., Makarov A., Hardman M., Cooks G.R. The Orbitrap: a new mass spectrometer // J. Mass Spectrom. 2005. Vol. 40. P. 430-443. DOI: 10.1002/jms.856
3. Nikolaev E., Sudakov M., Vladimirov G., et al. Multi-electrode harmonized kingdon traps // J. Am. Soc. Mass Spectrom. 2018. Vol 29, no. 11. P. 2173-2181. DOI: 10.1007/s13361-018-2032-9
4. Голиков Ю.К., Краснова Н.К., Соловьев К.В., Никитина Д.В. Интегрируемые ионные ловушки // Прикладная физика. 2006. № 5. C. 51-57.
5. Голиков Ю.К., Соловьев К.В. Электростатические ионные ловушки с разделением переменных в параболических координатах // Письма в ЖТФ. 2010. Т. 36, № 7. С. 82-88.
6. Голиков Ю.К., Соловьев К.В. Критерий поперечной устойчивости в ионных ловушках с интегрируемым в эллиптических координатах движением // Письма в ЖТФ. 2011. Т. 37, № 22. С. 43-49.
7. Соловьев К.В., Виноградова М.В. Об уточнении критериев устойчивости движения ионов в электростатических ловушках с интегрируемым движением // Тезисы VIII съезда ВМСО 9-13 окт. 2017. М., 2017. С. 110. URL: http://www. spsl.nsc.ru/FullText/konfe/VMSO _2017.pdf
8. Соловьев К.В., Виноградова М.В. Условия финитности движения иона в электростатической ловушке с разделением переменных в параболических координатах // Письма в ЖТФ. 2018. Т. 44, № 14. С. 34-41. DOI: 10.21883/PJTF.2018.14.46342.17215
9. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. М.: Наука, 1984. 833 с.
Санкт-Петербургский политехнический университет Петра Великого, Санкт-Петербург
(Соловьев К.В., Виноградова М.В.)
Институт аналитического приборостроения РАН, Санкт-Петербург (Соловьев К.В.)
Контакты: Соловьев Константин Вячеславович, k-solovyev@mail. ru
Материал поступил в редакцию 18.09.2018
ISSN 0868-5886
NAUCHNOE PRIBOROSTROENIE, 2019, Vol. 29, No. 1, pp. 143-148
UPDATING THE ION MOTION STABILITY CONDITIONS FOR ELECTROSTATIC TRAP INTEGRABLE IN ELLIPTIC COORDINATES
K. V. Solovyev12, M. V. Vinogradova1
1 Peter the Great St. Petersburg Polytechnic University, Saint-Petersburg, Russia 2Institute for Analytical Instrumentation of RAS, Saint-Petersburg, Russia
For ideal focusing traps electrostatic fields having charged particle motion separated in elliptic coordinates the ion motion stability conditions are updated. Analytical exact and approximate expressions are found for fields with zero and close to zero field parameter p values. These expressions connect initial data and field parameters to guarantee the ion motion finiteness.
Keywords: ideal space-time focusing, mass spectrometry, ion trap
REFERENСES
1. Gall' L.N., Golikov Yu.K., Aleksandrov M.L. et al. Vre-myaproletnyj mass-spektrometr. Copyright certificate of the USSR no. 1247973. [Time-of-flight mass spectrometer]. Prioritet 16.01.1985. (In Russ.).
2. Hu Q., Noll R., Li H., Makarov A., Hardman M., Cooks G.R. The Orbitrap: a new mass spectrometer. J.
Mass Spectrom., 2005, vol. 40, pp. 430-443. DOI: 10.1002/jms.856
3. Nikolaev E., Sudakov M., Vladimirov G., et al. Multi-electrode harmonized kingdon traps. J. Am. Soc. Mass Spectrom, 2018, vol. 29, no. 11, pp. 2173-2181. DOI: 10.1007/s13361-018-2032-9
4. Golikov Yu.K., Krasnova N.K., Solovyev K.V., Nikitina D.V. [Integrating electrostatic traps]. Prikladnaya fizika [Applied physics], 2006, no. 5, pp. 51-57. (In Russ.).
5. Golikov Yu.K., Solovyev K.V. [Electrostatic ionic traps with a separation of variables in parabolic coordinates]. Pis'ma v ZhTF [Applied Physics Letters], 2010, vol. 36, no. 7, pp. 82-88. (In Russ.).
6. Golikov Yu.K., Solovyev K.V. [Criterion of cross stabili-
ty in ionic traps with the driving integrated in elliptical axials]. Pis'ma v ZhTF [Applied Physics Letters], 2011, vol. 37, no. 22, pp. 43-49. (In Russ.). Solovyev K.V., Vinogradova M.V. [About specification of stability criterions of driving of ions in electrostatic traps with integrable driving]. Tezisy VIII s'ezda VMSO [Theses of the VIII congress of VMSO], 9-13 October 2017, Moscow, 2017, pp. 110. URL: http://www. spsl.nsc.ru/FullText/konfe/VMS0_2017.pdf (In Russ.).
Solovyev K.V., Vinogradova M.V. [Conditions of a finit-nost of driving of an ion in an electrostatic trap with a separation of variables in parabolic coordinates]. Pis'ma v ZhTF [Applied Physics Letters], 2018, vol. 44, no. 14, pp. 34-41. DOI: 10.21883/PJTF.2018.14.46342.17215 (In Russ.).
Korn G., Korn T. Mathematical Handbook for Scientists and Engineers: Definitions, Theorems, and Formulas for Reference and Review. General Publishing Company, 2000. 1151 p. (Russ. ed.: Korn G., Korn T. Spravochnik po matematike dlya nauchnyh rabotnikov i inzhenerov. Moscow, Nauka Publ., 1984. 833 p.). (In Russ.).
Contacts: Solovyev Konstantin Vyacheslavovitch, k-solovyev@mail. ru
Article received by editing board 18.09.2018
