Моделирование движения ионов в орбитальной ионной ловушке с учетом объемного заряда и неидеальности удерживающего электрического поля Текст научной статьи по специальности «Физика»

УДК 537.291

Г.Н. Владимиров1,2, А.В. Харченко3, Р. Хеерен3, Е.Н. Николаев1,2

1 Институт биохимической физики им. Н.М. Эмануэля РАН

2 Институт энергетических проблем химической физики РАН

3 Institute for Atomic and Molecular Physics, Netherlands

Моделирование движения ионов в орбитальной ионной ловушке с учетом

объемного заряда и неидеальности удерживающего электрического поля

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

Ключевые слова: масс-спектрометрия, орбитальная ионная ловушка, метод «частица в ячейке», объемный заряд.

Введение. Использование масс-спектрометрии для анализа биологических объектов предъявляет высокие требования к разрешению (более 106) и точности измерения масс (выше ppm). Таким высоким требованиям долгое время удовлетворяла только масс-спектрометрия ионного циклотронного резонанса с преобразованием Фурье (ИЦР-ФП) [1]. Однако, обладая рекордным разрешением и точностью измерения масс, этот метод не лишен некоторых недостатков, среди которых самым главным является необходимость использования крио-магнитов с относительно высокими (3 Тесла и выше) магнитными полями, что приводит к высоким эксплуатационным расходам при обслуживании магнитов. Поэтому большой интерес представляет другой тип фу-рье-масс-спектрометров — орбитальная ионная ловушка (orbitrap, орбитрэп), изобретенная относительно недавно [2]. Орбитальная ионная ловушка приближается по своим характеристикам к ИЦР-ФП приборам, использующим ячейки без динамической гармонизации [3]. Приборы, использующие ионные ловушки типа «орбитрэп», получили широкое распространение благодаря работам Макарова, который реализовал орбитальный подход к удержанию ионов, предложенный ранее Кингдоном [4].

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

Устройство ионной орбитальной ловушки. Работа масс-спектрометра на основе орбитальной ловушки основана на измерении наведенных токов при движении ионов в электростатическом поле, создаваемом потенциалом вида [5]:

где к, Кт, С — константы, причем Кт — характеристический радиус, к — кривизна поля;

ООО

г= хО + у2; х, у, г — декартовы координаты.

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

(1)

где К\, Я 2 — параметры, характеризующие геометрию орбитальной ловушки.

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

[5]: _______

ш = л/ (я/т)к (3)

— это частота продольных колебаний. Эта частота зависит от т/г, измерение именно этой частоты используется в масс-анализаторе орбитальной ионной ловушки:

“г = - 2

— это частота радиальных колебаний, где К — радиус круговой орбиты:

(4)

— это частота вращения вокруг центрального электрода.

Таким образом, чтобы задать геометрию орбитрэпа, необходимо выбрать коэффициенты к и Кт. Из литературы [7] известно, что для коммерчески доступной орбитальной ловушки К1 = 6 мм, Ко = 15 мм, Кт ~ \[2Ко; напряжение, прикладываемое к центральному электроду и = 3500 В. Зная эти величины можно вычислить к и частоту продольных колебаний. С учетом того, что внешний электрод заземлен, а на внутренний электрод подается напряжение, имеем: ^(К2) = 0^(К1) = и, тогда выразив 2 через К1 и К2, используя (2), имеем

к =____________2и____________ (5)

Кт 1п(§) - 0,5(К2 - К2)' ( )

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

относительно центрального электрода У^, которую можно определить из (4): т/О^‘р = к(КЛ-К2).

Из работы [7] известно также, что система инжекции настроена так, чтобы конечный радиус вращения относительно центрального электрода составлял величину в диапазоне К = 9-10 мм.

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

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

Начальные координаты — равномерное распределение ионов внутри сферы радиуса г = 0,12 мм (с координатами центра сферы К = 9 мм; 2 = К1 =6 мм) [7], всем ионам сообщалась скорость по направлению вращения ионов вокруг центрального электрода. Если скорость вращения вокруг центрального электрода в начальный момент времени будет отличаться от величины, даваемой выражением (4), то вращение вокруг этого электрода будет осуществляться не по круговым, а по эллиптическим траекториям. Поэтому скорость подбиралась так, чтобы эллиптичность траектории относительно центрального электрода была минимальна (рис. 1). На рис. 1 изображена схема численного эксперимента: модель движения ионов, которая была выбрана за основу для расчетов; область моделирования; проекции (рис. 1а). На рис. 1 также изображена эволюция ионных облаков в «орбитрэпе» для разных начальных скоростей вращения (рис. 1б, в), видно, что для начальной энергии вращения 1300 эВ эллиптичность траекторий значительно выше, чем в случае начальной скорости вращения 1600 эВ.

Частотные сдвиги. Одна из причин ограничения точности измерения масс при использовании орбитальной ловушки в качестве масс-спектрометра — это вызванные объемным зарядом отклонение частоты колебаний в потенциальной аксиальной яме от частоты колебаний в пределе одиночного иона. Поэтому повышение точности измерения масс требует понимания всех процессов, приводящих к таким систематическим отклонениям, для устранения их с помощью калибровки. В работе [8] экспериментально исследованы частотные сдвиги в орбитальной ионной ловушке, вызванные объемным зарядом с определением количества ионов по порядку величины при помощи automatic gain control (AGC). Исследование частотных сдвигов с помощью моделирования позволяет более детально, чем в эксперименте, изучить механизм частотных сдвигов, предоставляя точную информацию о количестве ионов.

Нами определены частотные сдвиги для 104; 105, 106; 2 106; 3 106; 4 106; 5 106; 6 106 зарядов в облаке в орбитальной ловушке. На рис. 2 изображены частотные сдвиги для инжекции при Up = 1600 эВ.

Рис. 2. Частотные сдвиги в орбитальной ионной ловушке: Я\ =6 мм, К2 = 15 мм, Ят ~ у/2я2 , и = 3500 В для 10 000; 105; 106; 2106; 3106; 4106; 5106; 6106 зарядов в облаке при т/г = 500 Да, 2 = 10 е

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

1600 эВ (в); т/г = 500 Да, 2 = 1 е

Использовалось программное обеспечение для моделирования трехмерного движения заряженных частиц в электро-магнитных полях, разработанное ранее в [9], основанное на методе «частица в ячейке» [10]. Уравнение Пуассона при вычислении электрического поля, создаваемого ионными облаками, решалось на трехмерной кубической сетке 64 х 64 х 64. Область моделирования имела размеры куба со стороной 3,81 см. Частотные сдвиги определялись для времени детектирования порядка 0.02 с, что соответствовало 100 000 шагам интегрирования движения. Величина шага подбиралась так, чтобы на один оборот вокруг центрального электрода приходилось порядка 100 шагов.

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

7^ + 772 + '''7П + т + + ''' + , где 51, 42, ..., 5П — виртуальные заряды; Р1, ро,

..., рп — заряды моделируемых частиц; гц — расстояния между рассматриваемым виртуальным зарядом и другими виртуальными зарядами или зарядом моделируемых частиц; ^ — это потенциал одного из двух электродов: р1 или р2. Кроме того, известно, что для системы виртуальных зарядов их сумма равна нулю, а сумма потенциалов детектирующих электродов, на которых они расположены, также равна нулю. В результате мы имеем систему п + 2 уравнений для п + 2 неизвестных, из которой находим величины зарядов, наводимых на каждом из электродов. Система уравнений решалась методом псевдоинверсии. Зная на каждом шаге величины зарядов, находящихся на детектирующих электродах, мы можем вычислить ток, наводимый ионами.

Для точного определения частот, соответствующих детектируемому сигналу по короткому куску сигнала, использовался метод фильтрации-диагонализации (ФДМ) [11]. В отличие от преобразования Фурье, в методе ФДМ используются некоторые предположения о природе порождающего сигнал процесса, например, предполагается, что отсчеты сгенерированы линейной ди-

намической системой. Создание этого метода было мотивировано задачей представления затухающих сигналов рядом затухающих по амплитуде синусоид. Гармонический затухающий сигнал, возникающий, например, при импульсном воздействии на колебательный контур, представляется в виде суммы затухающих синусоид, характеризуемых частотой, начальной фазой и добротностью, обратно пропорциональной коэффициенту затухания. Метод ФДМ представляет собой способ вычисления коэффициентов в разложении сигнала в ряд экспоненциально затухающих синусоид. Вычисление фазы, частоты и декремента затухания для членов этого ряда осуществляется с помощью анализа матрицы автокорреляции сигнала, строки и столбцы которой состоят из значений функции автокорреляции для разного положения сдвига сигнала (функция автокорреляции — это сумма произведений величин сигнала в разные моменты времени для разной величины сдвига). Такую матрицу можно привести к диагональному виду и найти ее собственные числа, зная которые можно определить фазы, частоты и декремент затухания [11]. Так как задача вычисления коэффициентов разложения в большом окне частот требует работы с большим базисом, то это порождает большую матрицу автокорреляции. Однако для такой матрицы вычисление собственных чисел оказывается вычислительно громоздким, поэтому используется поиск решения по небольшим частотным «окнам», а общий спектр «склеивается» из таких частичных решений, что уменьшает вычислительную сложность задачи, так как вычислительная сложность растет нелинейно с ростом размера матрицы. В случае более точного описания движения отдельного иона в орбитальной ионной ловушке необходимо также учитывать влияние объемного заряда ионного облака на рассматриваемый ион. Если предположить, что ионное облако дает добавку в потенциал вида р = -С1р5г2, где р — плотность распределения заряженных частиц, С1 — константа, характеризующая размер и форму ионного облака; тогда, добавив это слагаемое в уравнение (1) и проведя замену обозначений, мы можем привести систему к старому виду. В результате для измеряемой частоты, которая описывается выражением (3), получим: “ = 2п/ = д/(д/т)к* = д/(д/т)к - 2С15р. Тогда для зависимости т/г от измеряемой частоты имеем: т/г = толтгг—2с1Рд , что соответствует калибровочному закону т/г к С2//2 - NС'3//4,

( ’ 1+ (2Л/)2

где N — число зарядов в ионном облаке, С\ — константы.

Из результатов моделирования (рис. 2) можно сделать заключение, что частотные сдвиги аппроксимируются зависимостью вида f = 277789 - 2,88N 10-6 Гц. Это показывает, что сдвиги частоты можно аппроксимировать выражением вида f = С4 - NC5, что соответствует калибровочному закону т/г к р - .

Аппроксимация, основанная на экспериментальных данных [8], дает т/г = С§//2-\^Сд//2.

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

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

В реальном приборе длина электродов ограничена по 2, и поэтому создаваемый потенциал будет отличаться от потенциала, даваемого выражением (1), на некоторую величину. Будем аппроксимировать отклонение поля орбитрэпа от идеального с помощью потенциала электростатического поля в пространстве, не содержащем зарядов. Аппроксимировать электростатический потенциал в области, свободной от зарядов, можно только с помощью функций, являющихся решениями уравнения Лапласа. Электростатическое поле, выраженное с помощью функций электростатического потенциала, описывается благодаря уравнениям Максвелла. В области с нулевым суммарным зарядом уравнения Максвелла сводятся к уравнению Лапласа для электростатического потенциала [12]. Решения уравнения Лапласа для систем, обладающих симметрией

вращения (то есть в нашем случае), удобно выражать в виде ряда по сферическим гармоникам вида [12]:

то I

р(г,вр) = ЕЕ [А,тГ1 + Вт г-((+1) ]У1т (М, (6)

1=0 т=-1

где У[т — сферические гармоники, Л[т, Вт — константы.

АА=0.01 ;В=0; ю3зарядов БА=0.03;В=0; 103 зарядов В А=0.05 В=0; 103 зарядов -----------------1 " >-

{ детектируемый

" сигнал

детектируемый сигнал

t=0.04

Ь (О)-^ (©

й (©)

Х,¥''-'

Г А=0;В=0.001; 103 зарядовД А=0;В=0.001; 105 зарядов Е А=В=0; 105 зарядов

детектируемый детектируемый

сигнал сигнал

1=0.04 1=0.04

сек сек

(о

ХЛ''----^

е о

Рис. 3. Детектируемый сигнал и форма ионных облаков для разных моментов времени. Поведение ионных облаков при разной величине вклада слагаемых Ак (2г3 - 3гг2) и (8г4 - 24г2г2 + 3г4) в суммарное поле для орбитальной ионной ловушки. Измерение в течение 0.05 с при Л = 0,01 (а); Л = 0,03 (б); Л = 0,05 (в);

В = 0,001 (г, д); Л = В = 0 (е). Используется 1000 частиц в облаке (а, б, в, г), используется 100 000 частиц в облаке (д, е), т/г = 500 Да, 2 = 1 е

Заметим, что | (г2 - ^) — слагаемое потенциала идеального орбитрэпа (1) является членом ряда (6) Л20г2У2о. Для простоты будем аппроксимировать рассматриваемую неидеальность поля членами Л30г3У30 и Л40г4У40 ряда (6).

Рассмотрим динамику движения ионов в орбитальной ловушке с потенциалом вида

2

\ к К2 1п / \ + Лк (2г3 _ 3гг2) + ,

ь2

Р(г,г) = 2 (г2 - у) + 2Кт111

Кт

+ Лк (2г3 - 3гг2) + ^ (8г4 - 24г2г2 + 3г4) + С,

Ь Т 2

^2Ко, и = 3500 В; Л, В

где к — это коэффициент (5) для К1 =6 мм, К2 = 15 мм, Я„

безразмерные коэффициенты, Ь — длина моделируемой области.

Можно видеть, что наличие в потенциале орбитальной ловушки дополнительных слагаемых вида Вк (8г4 - 24г2г2 + 3г4) и ^ (2г3 - 3гг2) приводит к разрушению ионных облаков (рис. 3а,

б, в, г). Видно также, что чем больше вклад этих слагаемых в суммарное электростатическое поле, тем быстрее происходит расфазировка ионов в ионном облаке вдоль оси 2 и тем сильнее вызванное этим ограничение на разрешающую способность.

При увеличении числа частиц в ионном облаке с 1000 до 100 000 (m/z 500, Z = 1) при наличии в электростатическом потенциале орбитальной ловушки слагаемого jj (8z4—24z2r2+3r4) происходит ускорение расфазировки ионов в ионных облаках вдоль оси Z (рис. 3 г, д). Также видно, что при отсутствии дополнительных слагаемых в потенциале ионное облако, состоящее из 100 000 частиц (m/z 500, Z = 1) устойчиво (рис. 3е).

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

m/z ~ f 2 — f 3 .

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

Работа была Выполнена при поддержке РФФИ (Гранты 10-04-13306-РТ_оми, 09-04-00725-а, 09-03-92500-ИК_а), CRDF (Грант RUC1-2941-M0-09), Министерства образования и науки Российской Федерации (Федеральная целевая программа «Научные и научно-педагогические кадры, инновационной России на 2009-2013 годы», ГК 14.740.11.0755, 16.740.11.0369), Российской академии наук (Программа Президиума РАН «Фундаментальные науки - медицине»).

Литература

1. Marshall A.G., Hendrickson C.L., Jackson G.S. Fourier transform ion cyclotron resonance mass spectrometry: A primer // Mass Spectrom. Rev. — 1998. — V. 17. — P. 1-35.

2. Makarov A., Denisov E., Lange O., Horning S. Dynamic Range of Mass Accuracy in LTQ Orbitrap Hybrid Mass Spectrometer // J. Am. Soc. Mass. Spectrom. — 2006. — V. 17. — P. 977-982.

3. Boldin I.A., Nikolaev E.N. Fourier transform ion cyclotron resonance cell with dynamic harmonization of the electric field in the whole volume by shaping of the excitation and detection electrode assembly // Rapid Commun. Mass Spectrom. — 2011. — V. 25. — P. 122-126.

4. Kingdon K.H. A Method for the Neutralization of Electron Space Charge by Positive Ionization at Very Low Gas Pressures // Phys. Rev. — 1923. — V. 21. — P. 408-418.

5. Makarov A. Electrostatic Axially Harmonic Orbital Trapping: A High-Performance Technique of Mass Analysis // Anal. Chem. — 2000. — V. 72. — P. 1156-1162.

6. Lewis R.R. Motion of ions in the Kingdon trap // J. Appl. Phys. — 1982. — V. 53(6). — P. 3975--3980.

7. Makarov A, Denisov E., Lange O. Performance Evaluation of a High-field Orbitrap Mass Analyzer // J. Am. Soc. Mass. Spectrom. — 2009. — V. 20. — P. 1391-1396.

8. Gorshkov M.V., Good D.M., Lyutvinskiy Y., Yang H., Zubarev R.A. Calibration function for the Orbitrap FTMS accounting for the space charge effect // Journal J. Am. Soc. Mass. Spectrom. — 2010. — V. 21(11). — P. 1846-1851.

9. Nikolaev E.N., Heeren R.M.A., Popov A.M., Pozdneev A.V., Chingin K.S. Realistic modeling of ion cloud motion in a Fourier transform ion cyclotron resonance cell by use of a particle-in-cell approach // Rapid Commun. Mass Spectrom. — 2007. — V. 21. — P. 3527-3546.

10. Harlow F.H. The particle-in-cell method for numerical solution of problems in fluid dynamics

// Methods in Computational Physics. — 1964. — P. 319-343.

11. Mandelshtam V.A., Taylor H.S. Harmonic inversion of time signals and its applications // J.

Chem. Phys. — 1997. — V. 107. — P. 6756-6759.

12. Jackson J.D. Classical Electrodynamics. — 3rd Ed. John Wiley & Sons, 1998.

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