Суперкомпьютерное исследование движения ионов в ловушках Кингдона и Пеннинга с полным учетом кулоновского взаимодействия тема автореферата и диссертации по физике, 01.04.17 ВАК РФ

Рюмин, Павел Александрович АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2011 ГОД ЗАЩИТЫ
   
01.04.17 КОД ВАК РФ
Диссертация по физике на тему «Суперкомпьютерное исследование движения ионов в ловушках Кингдона и Пеннинга с полным учетом кулоновского взаимодействия»
 
Автореферат диссертации на тему "Суперкомпьютерное исследование движения ионов в ловушках Кингдона и Пеннинга с полным учетом кулоновского взаимодействия"

Российская академия наук Институт энергетических проблем химической физики

На правах рукописи

005002719 <2._____________1

РЮМИН ПАВЕЛ АЛЕКСАНДРОВИЧ

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

01.04.17 - химическая физика, горение и взрыв, физика экстремальных состояний вещества

автореферат диссертации на соискание ученой степени кандидата физико-математических наук

1 7 НОЯ 2011

Москва, 2011г.

005002719

Работа выполнена в Учреждении Российской академии наук Институте энергетических проблем химической физики РАН

Научный руководитель:

доктор физико-математических наук, профессор Николаев Евгений Николаевич

Официальные оппоненты: доктор физико-математических наук Разников Валерий Владиславович

доктор физико-математических наук Трахтенберг Леонид Израилевич

Ведущая организация:

Московский физико-технический институт(государственный университет)

Защита состоится «07» декабря 2011г. в 11 час. 00 мин. на заседании диссертационного совета Д 002.112.01 при Институте энергетических проблем химической физики Российской академии наук по адресу: 119334, г. Москва, Ленинский проспект, д. 38, корп.2, ИНЭП ХФ РАН.

С диссертацией можно ознакомиться в библиотеке Института химической физики им. H.H. Семенова Российской академии наук.

Автореферат разослан «3 » ноября 2011г.

Ученый секретарь Г\

диссертационного совета Д 002.112.01 \ \

кандидат физико-математических наук Ларичев М.Н.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ.

Введение. Актуальность проблемы.

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

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

масс по известным соотношениям, связывающим частоту с отношением массы к заряду ионов. Большое количество ионов в ионных ловушках типа Пеннинга (ИЦР ПФ) и Кингдона (Орбитрэп) приводит к таким явлениям как коалесценция пиков в масс-спектрах и разрушение ионных облаков из-за их взаимодействия, что не позволяет достигать высокую разрешающую способность одновременно с достижением большого динамического диапазона. Показано, что при наличии в ловушке ионов с разным отношением массы к заряду кулоновское взаимодействие между образуемыми ими облаками приводит к частотным сдвигам, которые ограничивают точность измерения масс. Даже в том случае, когда в ловушке изолируются только ионы с одним отношением массы к заряду, при значительном радиусе возбуждения взаимодействие со стенками электродов так же приводит к сдвигу циклотронной частоты [1]. Подобные эффекты наблюдаются также и в ловушке Кингдона[2].

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

До недавнего времени моделирование одновременного движения большого числа ионов было невозможно из-за недостаточной производительности компьютеров. Первые работы по моделированию были посвящены движению отдельных ионов в различных электрических полях. Появление коммерческой программы 81МЮЫ [3] значительно упростило понимание различных аспектов движения единичных ионов в ионных ловушках произвольной геометрии. Эта программа способна

рассчитывать траектории отдельных частиц в различных измерительных ячейках ИЦР с довольно высокой точностью. Однако программа не предоставляет достаточно адекватных возможностей учитывать кулоновское взаимодействие между ионами.

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

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

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

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

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

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

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

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

Личный вклад автора. Автор внес основной вклад в разработку алгоритма по учету взаимодействия ионов с наведенными зарядами-изображениями на электродах ионных ловушек произвольной геометрии и его программную реализацию, а так же алгоритма, который учитывает симметрию задачи. Исследование движения ионных ансамблей в ловушке Кингдона и возможности аксиального детектирования в ловушке Пеннинга с динамической гармонизацией выполнено совместно с И.А. Бодциным (ИНЭП ХФ РАН).

Защищаемые положения.

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

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

Обнаружены эффекты коллективного взаимодействия ионов в ионных ансамблях в ловушке Кингдона.

Исследована возможность аксиального детектирования в ловушке Пеннинга с динамической гармонизацией как потенциального масс-спектрометрического метода.

Апробация работы. Результаты работы докладывались на следующих российских и международных конференциях: 58-ая Ежегодная конференция американского масс-спектрометрического общества «Масс-спектрометрия и смежные темы», Солт Лейк Сити, США, 23-27 мая 2010; 4-я Всероссийская конференция «Фундаментальные вопросы масс-спектрометрии, и ее аналитические применения», Звенигород, Россия, 10 -14 октября 2010; 57-ая Ежегодная конференция американского масс-спектрометрического общества «Масс-спектрометрия и смежные темы»,

Филадельфия, США, июнь 2009, 56-ая Ежегодная конференция американского масс-спектрометрического общества «Масс-спектрометрия и смежные темы», Денвер, США, июнь 2008, 55-ая Ежегодная конференция американского масс-спектрометрического общества «Масс-спектрометрия и смежные темы», Индианаполис, США, июнь 2007, 8th European FTMS conference, Москва, август 2007 .

Публикации. Основное содержание диссертации опубликовано в работах, список которых приводится в конце автореферата. Работы, результаты которых изложены в диссертации, выполнены при поддержке грантов РФФИ.

Структура и объем диссертации. Работа изложена на 93 страницах, содержит 32 рисунка и 3 таблицы. Диссертация состоит из пяти глав, заключения и списка цитируемой литературы. ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ.

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

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

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

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

Обзор существующих методов компьютерного моделирования.

Наиболее простым подходом к решению задачи о движении ионов с учетом их взаимодействия является метод РР (Particle-Particle, частица частица). Ион-ионное взаимодействие учитывается путем расчета электрического поля, создаваемого каждым ионом в точках, где находятся

все остальные ионы. При этом взаимодействие с зарядами, наводимыми на электродах, вообще не учитывается, то есть метод целесообразно использовать, только если ионы не подходят близко к стенкам электродов. Другой недостаток этого метода - быстрый рост объемов вычислений с ростом числа частиц (0(N2), N - число частиц). Обычно метод используется, если количество частиц составляет несколько тысяч; для расчета движения десятков тысяч частиц с помощью метода РР используется распараллеливание вычислений на суперкомпьютерах. Метод, позволяющий избежать этой проблемы и работать с большим количеством частиц - это метод PIC (Particle-in-Cell - частица в ячейке).

Метод PIC известен с 1955 г. и широко применяется в физике плазмы. Идея метода состоит в том, чтобы не рассчитывать ион-ионные взаимодействия по одному (как в методе РР), а интерполировать заряды в узлы вычислительной сетки. Тогда количество зарядов получается равным числу узлов сетки. Затем на этой сетке решается уравнение Пуассона . Для этого могут быть использованы различные алгоритмы, но самым быстрым при использовании параллельных вычислений является метод быстрого преобразования Фурье (БПФ). БПФ является прямым методом, то есть после выполнения определённого числа шагов, он дает результат, точность которого зависит от мелкости сетки разбиения, в отличие от итерационных методов, где каждый шаг алгоритма уменьшает ошибку в определенное число раз. Количество вычислительных операций в методе PIC пропорционально (N + Ng log Ng), где N - количество частиц, a Ng -число узлов сетки.

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

Обзор предыдущих работ

Одной из первых попыток моделировать кулоновское взаимодействие ионов в ловушке была работа группы Николаева [4]. Моделировалось поведение системы, включающей до 1024 ионов на ЭВМ с 1024 процессорами (25 MHz SPARC, 16 MB RAM). Поле ловушки вычислялось в узлах сетки и интерполировалось на частицы; кулоновские взаимодействия рассчитывались по методу РР, взаимодействие с зарядами-изображениями не учитывалось; число шагов по времени было ограничено 50000. Для расчета большего количества частиц использовалось приближение квази-ионов (замену группы ионов одним многозарядным ионом).

Следующим шагом в области моделирования систем большого числа частиц, заключенных в ловушку ИЦР ПФ масс-спектрометра, была работа по двумерному моделированию, выполненная Митчеллом и Смитом [5]. Для своего двумерного кода они использовали модель PIC, до этого не применявшуюся в масс-спектрометрии. Возможности программы были продемонстрированы при исследовании вопросов, связанных с равновесием ионных облаков в отсутствии возбуждения их циклотронного движения. Первое реалистичное трехмерное моделирование ИЦР ПФ масс-спектрометра было выполнено в работе Митчелла [6], в которой он представил новый трехмерный код PIC3D как расширение двумерного кода. Численные эксперименты включали до 500000 частиц, выполнялось 100000 итераций; расчеты проводились на рабочей станции на базе процессора 500 MHz Dec/Alpha и занимали от одного до четырех дней.

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

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

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

частицы, находящиеся в ш

окружении электродов, Ж^ЛШш^Щ{

индуцируют заряды на Ш

поверхности этих электродов,

таким образом ЧТО потенциал Рис' 1а(слевя)- Ловушка Киигдоиа (ОгЫП-ар) в разрезе. Рис.

* * 1б(справа) Поверхности, образующие электроды, помещаются в

соответсвующуга прямоугольную вычислительную сетку, и ОКаЗЫВаеТСЯ ПОСТОЯННЫМ ВДОЛЬ электроды заменяются набором точек. На рисунке заменен

только внутренний электрод.

электродов, что требует

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

Рассмотрим метод емкостной матрицы на примере ловушки Кингдона (рис. 1а). Сначала нужно поместить поверхности электродов в прямоугольную вычислительную сетку (на этой сетке будет решаться уравнения Пуассона методом БПФ) и заменить электроды емкостями (рис. 16).

Каждый элемент поверхности (емкость) моделируется двумя близко лежащими точками - в одной будет вычисляться потенциал , в другую помещаются заряды. Для того, чтобы найти емкостную матрицу, нужно разместить в соответствующую точку каждого элемента поверхности единичный заряд и вычислить методом БПФ потенциалы, наведенные этим зарядом на каждую емкость и повторить эту процедуру для каждой емкости. Значения потенциалов в соответствующих точках образуют столбец матрицы А (рис. 2).

(1)

Ж (р 1"

42 Ч>2

Как видно, матрица А есть матрица обратная матрице С:

<Р\

1-1

А

4г 42

А

-1

<Р 2

<71 42

с

(рх <Р2

(2)

(3)

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

Полный алгоритм расчета движения ионов на каждом временном шаге выглядит так:

1) Интерполяция зарядов ионов в узлы вычислительной сетки;

2) Решение уравнения Пуассона методом БПФ;

3) Экстраполяция значений потенциалов из узлов сетки в точки емкостей (алгоритм идентичен алгоритму интерполяции поля в точки, где находятся

ионы);

4) Определение столбца зарядов емкостей путем умножения матрицы С на столбец потенциалов, представляющий собой разность изначально заданного потенциала электрода и потенциала, полученного в пункте 3;

5) Интерполяция зарядов ионов и зарядов емкостей на узлы вычислительной сетки;

6) Решение уравнения Пуассона методом FFT с учетом как ионов, так и зарядов емкостей;

7) Вычисление напряженности электрического поля в узлах вычислительной сетки;

8) Экстраполяция напряженности поля на ионы;

9) Интегрирование уравнений движения;

Помимо электрического поля, .

создаваемого самими ионами, в масс-

спектрометрических устройствах надо

учитывать поле электродов: на электроды

может подаваться переменное или постоянное

напряжение, которое создает дополнительное

электрическое поле в ловушке. При

вычислении суммарного электрического поля Рис. 2 На рисунке схематически

изображено приближение электродов МОЖНО использовать два разных подхода. 2-мя наборами точек и вычисление

матрицы А. Черные сплошные точки Один из НИХ был описан выше: значения — места, где измеряется потенциал,

незакрашенные точки — места куда

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

' (1' ' Ol '

In =

_ 0

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

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

Для проверки описанного алгоритма было проведено моделирование движения ионов в кубической ячейке ИЦР, и его результат был сравнен с моделированием методом PIC.

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

приводятся характерные времена вычислений на каждом шаге, которые показывают, что с определенного момента процедура решения системы линейных уравнений(СЛАУ) является лимитирующей.

Количество емкостей, аппроксимирующих электроды 294 726 1536 3174 5766

Решение уравнения Пуассона методом FFT 0,04с 0,04с 0,04с 0,04с 0,04с

Решение СЛАУ 0,005с 0,015с 0,045с 0,3с 0,9с

Остальные процедуры 0,2с 0,2с 0,2с 0,2с 0,2с

Суммарное время 0,25с 0,26с 0,29с 0,54с 1,14с

Табл. 1 Характерные времена вычисления различных процедур, выполняемых на каждом шаге, измеренные в секундах. Вычисления производились на 1 процессоре 1,1 Ггц серверной системы IBM eServer pSeries 690

В третьей главе предлагается улучшенный емкостной метод,

И

1« 4* 4 • 3»

"•R' I \

4* 1« > 4«

который позволяет использовать симметрию задачи.

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

учетом этой симметрии можно существенно 2< ^

сократить вычислительную сложность задачи. Рис. з Преораювания

индексов, которые образуют Основная идея СОСТОИТ В ТОМ, чтобы абелеву группу С4

преобразовать существующую СЛАУ, решить полученную систему со

значительно меньшими вычислительными затратами и преобразовать

обратно полученное решение. Впервые идея использовать симметрию в

задачах такого рода была предложена Алговером [8].

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

Матрица Л называется эквивариантной под действием группы С,

если

Ац = Л1я ;я, для всех ¿,У £ / и д в С, (4)

где I - множество индексов, а С - группа симметрии.

При выборе схемы дискретизации, которая «_ учитывает симметрию матрица А эквивариантна под ' •

действием группы, описывающей эту симметрию.

требуется место только под определенные столбцы, рис. 4 Набор

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

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

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

Отсюда следует, что для хранения матрицы А

независимых индексов.

Где х, - элемент столбца, С - группа |С|, порядок группы, Л множество неприводимых представлений, с![( - степень представления. Для каждого 1 6 I и g 6 С выполняются уравнения:

Где trace(X|(R)) - след матрицы, в случае абелевой группы равный самому элементу Xj(R), Xj - элемент столбца, G - группа, |G| - порядок группы, Л - множество неприводимых представлений, dR - степень представления.

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

где (С ( - порядок группы, Л - множество неприводимых

представлений, с!к - степень представления, I - множество индексов, 5 -

множество независимых индексов.

Применяя (10), мы можем привести исходную систему(1) к блочно-диагональному виду:

/ЛЙ1 - 0 \ ^ фг

! ! х : = : , (9)

\ о ••■ АПп) Чп Фп

Где АКк - матрица размером 5x5. Такая система уравнений

решается значительно быстрее исходной. Матрицы АЛ|( содержат только

действительные числа, тогда как преобразованные столбцы потенциалов и

зарядов могут быть мнимыми, так что для каждой матрицы Айк, требуется

решить 2 системы уравнений.

х(Я,ё0 = ^Жги),

Тогда обратное преобразование может быть определено как:

(6)

(7)

(8)

Общий алгоритм решения СЛАУ выглядит так:

• Выполнить обобщенное дискретное преобразование Фурье от матрицы

А из (1) воспользовавшись (8). Этот шаг может быть выполнен на предварительной стадии, а его результат сохранен в память компьютера.

• Выполнить обобщенное дискретное преобразование Фурье от столбца потенциалов на каждом шаге, воспользовавшись (5).

• Решить 2п систем уравнений (одна для действительной части, другая для мнимой), где п - количество неприводимых представлений (в случае абелевой группы равное порядку группы)

• Сложить действительную и мнимую части решения.

• Используя (6) восстановить все значения столбца ц

• Воспользовавшись (7), посчитать обратное преобразование столбца ц

Следует отметить, что данный алгоритм не накладывает никаких

ограничений на столбец потенциалов из (1).

Количество точек Прямое решение СЛАУ Используя симметрию задачи

2500 0.016с 0.004с

6400 0.17с 0,016с

10000 0,38с 0,028с

Табл. 2 Времена решения СЛАУ в случае использования симметрии задачи и в случае прямого вычисления.

Количество точек Прямое решение СЛАУ Используя симметрию задачи

2500 146 Мб 1.2 Мб

6400 952Мб 10Мб

10000 2.3 Гб 20 Мб

Табл. 3 Использование оперативной памяти при решении СЛАУ в случае использования симметрии

задачи и в случае прямого вычисления.

Этот алгоритм запрограммирован на языке Фортран. Сравнительные времена вычисления операции на каждом шаге и используемая память приведены в табл.2 и табл. 3 соответственно.

17

В данном примере использовалась схема дискретизации, применимая при аппроксимации цилиндрической ячейки ИЦР. Аналогичная схема дискретизации может быть использована при аппроксимации ячейки Пеннинга с динамической гармонизацией, с заданными другим образом граничными условиями. В случае 2500 точек использовалось 25 независимых индексов, в случаях 6400 и 10000 80 и 100 соответственно.

В четвертой главе представлены результаты моделирования движения различного количества ионов в ловушке Кингдона.

Ловушка Кингдона моделировалась 1600 емкостями. Вычисления проводились в кубе со стороной 40,6 мм. Форма электродов задавалась уравнением:

([ - номер электрода, 1 соответствует внешнему электроду, 2 -внутреннему), Ш = 15 мм, Ш = 6 мм, К] = 20,3 мм. Поле электродов задавалось аналитически формулой:

к=1,412-107 В/м. Напряжение между электродами 2 кВ. Для исследования взаимодействия ионных облаков близких масс в ловушку запускались ионы двух масс: 100 Да и 100,5 Да. Начальная скорость ионов была подобрана таким образом, чтобы они двигались по окружности (в координатах ху) - 4,65-104 м/с (соответствует энергии 1,12 кэВ). Был использован временной шаг 3-Ю"8 с. Для получения временного сигнала вычислялась разность суммы зарядов на емкостях одной половины ловушки и суммы зарядов противоположной половины. Рассчитанное время движения 4 мс.Для исследования влияния эффектов объемного заряда на движение ионных облаков было рассчитано движение ионных облаков, состоящих из различного количества ионов. При 1000 ионов

(10)

(И)

каждой массы движение облаков, временной сигнал и Фурье-спектр выглядят совершенно невозмущенными ион-ионными взаимодействиями (рис. 5). Для моделирования большего количества ионов нами была применена модель квази-ионов, то есть 10б ионов с массой 100 и зарядом 1 моделировались 103 ионами с массой 105 и

Рнс. 5 Ионные облака(сверху проекция XV, снизу \/. 1, временном сигнал и частотный спектр 2000 ионов.

зарядом 10 . Адекватность модели была проверена путем сравнения динамики движения ионных облаков в случае 103 квази-ионов с массой 105 зарядом 103 и в случае 104 квази-ионов с массой 104 зарядом 102; расчет показал, что движение облаков в этих случаях идентично. На рис. 6 представлены результаты расчета с квазиионами массой 105 и зарядом 103, моделирующими 106 ионов каждой массы. При таком количестве ионов эффекты объемного заряда играют существенную роль. Ионное облако разделяется на три части: две - это кольца с равномерным

распределением ИОНОВ, а третья Рис. 6 Ионные облака(сверху проекция

XV, снизу X/.). временной сигнал и представляет собой небольшой СГуСТОК. частотный спектр 2 млн. ионов.

Этот сгусток состоит из частиц обеих масс и в течение эксперимента

разделяется на два сгустка, намного медленнее, чем основная масса частиц

разделяется на кольца. Сгусток двигается с частотой отличной от частоты

движения колец, образованных ионами обеих масс и на Фурье-спектре

дает третий пик в дополнение к пикам, соответствующим ионам двух

измеряемых масс.

I

л

1 .

' ™

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

Для того, чтобы измерять частоту аксиального движения нужно, чтобы эта частота не зависела от амплитуды. Известно поле, в котором частота не зависит от амплитуды - это гиперболическое поле :

Ф(г, г) = ^у(2г2 — г2), (12)

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

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

Болдиным и Николаевым[10.11] был предложен метод, в котором создается не настоящее гармоническое поле в ячейке, а поле, которое становится гармоническим после усреднения по циклотронному движению ионов. Такое поле получается при сегментировании боковых

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

ЗаЧйМЯегтьШ эпекград

Запирающее напряжение

Электрод сферпческоЛ

Рис, 7 Основные принципы гармонизации, сегментирование боковых электродов, торцевые электроды сферической формы

На торцевые электроды и вогнутые электроды подается запирающий потенциал V, а центральные выпуклые электроды заземляются. Такая

конфигурация делает возможным

ЯЧеИКе Рис. 8 Масс-,

спектр и сигнал

аксиальное детектирование

полученный в численном

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

динамической гармонизацией.

УСЛОВИЯ усреднения. Так как частота Разрешающая способность 104,

время измерения сигнала 0.88 с

аксиальных колебаний не зависит от

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

Возможность использования такого типа измерительной ячейки в масс-спектрометрах ИЦР ПФ была подтверждена численными экспериментами. Была продемонстрирована разрешающая способность

104. Спектр и сигнал изображены на рис. 8. В численном эксперименте моделировалось движение ионов массой 1000 и 1000.1 Да. Магнитное поле в эксперименте было равным 5Т. Длина ячейки 150мм, внутренний диаметр равен 56мм. На запирающие электроды подавался потенциал равный 1кВ.

Рис. 9 Схемы измерительной ячейки при различных количествах секторов.

В предложенном масс-спектрометре ионы удерживаются в аксиальном направлении гармонизированным запирающим потенциалом, а в радиальном направлении магнитным полем. Напряженность магнитного поля должна быть достаточной для усреднения запирающего электрического поля и удержания ионов в радиальном направлении. Для удержания ионов с большим m/z требуется большее магнитное поле при одинаковом запирающем потенциале. При заданном магнитном поле и запирающем электрическом потенциале в обычной ячейке Пеннинга ионы с отношением m/z больше определенного не удерживаются. Эта зависимость выражается формулой [12]:

Г/2) крИт. = ^ (13)

где q - заряд, а - длина ячейки, В0 - величина магнитного поля, V -запирающее напряжение, а а - эмпирический коэффициент, зависящий от геометрии ячейки.

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

Ячейку Пеннинга с динамической гармонизацией можно

ИЗГОТОВЛЯТЬ С раЗНЫМ р,)с. 10 Зависимость критической массы от запирающего

потенциала, при разном количестве секторов. Верхняя количеством боковых СеКТОрОВ, штриховая линия - теоретический предел удержания ионов,

далее предел удержания ионов в ловушках с 32,16,8,4

Два варианта показаны на рис. 9. секторамии соответственно.

Для ловушек с различным количеством секторов была промоделирована зависимость критической массы от запирающего напряжения(рис. 10). Магнитное поле в численном эксперименте было равно 5Т, длина ячейки 150мм. В заключении приводятся основные выводы:

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

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

3) Исследовано влияние объемного заряда на движение ионов в ловушке Кингдона.

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

5) Выявлены особенности движения больших ионных ансамблей в ловушке Кингдона.

Цитируемая литература

1. Wong, R. L.; Amster, I. J., Int. J. Mass Spectrom. 2007,265, 99-105.

2. Makarov, A.; Denisov, E.; Lange, О.; Horning, S., J. Am. Soc. Mass Spectrom. 2006, 17,977-982.

3. Dahl DA. SIMION 3D version 7.0 User's manual, 1995.

4. Nikolaev EN, Miluchihin NV, Inoue M. Int. J. Mass Spectrom. Ion Processes 1995; 148: 145.

5. Mitchell, D. W.; Smith, R. D. Int. J. Mass Spectrom. Ion Processes 1997, 165/166, 271-297.

6. Mitchell, D. W. J. Am. Soc. Mass Spectrom., 1999,10,136-152.

7. Nikolaev E.N., Heeren R.M.A., Popov A.M., Pozdneev A.V., Chingin K., , Rapid Communications in Mass Spectrometry 21 (2007) 3527-3546.

8. E. L. Allgower, K. Böhmer, К. Georg and R. Miranda, Exploiting symmetry in boundary element methods, SIAM J. Numer. Anal., 29 (1992), pp. 534552.

9. Van Dyck, R.S., Schwinberg, P.B.; Dehmelt, H.G. Source: Physical Review Letters, 14 Feb. 1977, p 310-14

lO.Ivan A. Boldin, Eugene N. Nikolaev. Rapid Communications in Mass Spectrometry, Vol. 25, p. 122-126, 2011.

11 .Eugene N. Nikolaev, Ivan A. Boldin, Roland Jertz and Gökhan Baykut. J. Am. Soc. Mass Spectrom. Published online April 2011, DOI: 10.1007/sl3361-011-0125-9

12.Marshall, A. G.; Hendrickson, С. L.; Jackson, G. S. Fourier transform ion cyclotron resonance mass spectrometry: a primer. Mass Spectrom. Rev., 1998, 17, 1-35.

Публикации:

13.Рюмин П. А., Болдин И. А., Автономов Д. М., Николаев Е. Н.Метод емкостей в моделировании движения ионных ансамблей в ионных ловушках и системах транспорта ионов с электродами произвольной формы. Труды МФТИ, Том 3, No. 3(11), стр. 17-21,2011.

Н.Рюмин П. А., Болдин И. А., Николаев Е. Н. Аксиальное детектирование в ловушке Пеннинга. Четвертая Всероссийская конференция «Фундаментальные вопросы масс-спектрометрии и ее аналитические применения». Звенигород, Россия, 10 -14 октября, 2010.

15.Pavel Ryumin; Ivan Boldin ; Eugene Nikolaev. New possibilities of capacity method improvement for ion clouds dynamics simulations in traps with arbitrary geometry electrodes. 58th Amer. Soc. Mass Spectrom. Annual

Conf. on Mass Spectrometry & Allied Topics, Salt Lake City, UT, May 23 -27, 2010.

16.Pavel Ryumin, Eugene Nikolaev. New method of ion motion simulation for arbitrary electrodes geometry with total account for space and image charge interactions. 59th Amer. Soc. Mass Speetrom. Annual Conf. on Mass Spectrometry & Allied Topics, Denver, Colorado, USA, 2011.

17.Pavel Ryumin, Eugene Nikolaev, Alexander Pozdneev, Ivan Boldin, Dmitry Avtonomov. The comparison of 3 different approaches of realistic ion motion modeling in cubic ICR cell ("Particle-particle" approach, "Particle-in-cell" approach, "Capacity method" approach). 8th European FTMS Conference, Moscow, Russia, August 28 - September 1, 2007.

18.Eugene Nikolaev; Ivan Boldin; Pavel Ryumin; Igor Popov. Characterization of axial motion frequency analyses in harmonized cylindrical Penning traps as potential mass spectrometry method. 58th Amer. Soc. Mass Speetrom. Annual Conf. on Mass Spectrometry & Allied Topics, Salt Lake City, UT, May 23-27,2010.

19.Eugene Nikolaev, Ivan Boldin, Pavel Ryumin, Gleb Vladimirov, Ron M.a. Heeren, Alexander Pozdneev, Dmitry Avtonomov. New approaches to supercomputer modeling of fields and ion cloud dynamics with total account for ion-ion and image charge interactions. 56th American Society for Mass Spectrometry Conference. Denver, Colorado, USA, June 1-5, 2008.

20.Дмитрий Автономов, Иван Болдин, Павел Рюмин, Александр Позднеев, Евгений Николаев. Метод ёмкостей в моделировании движения ионных ансамблей в ионных ловушках и системах транспорта ионов с электродами произвольной формы. 3-ий съезд ВМСО, 2-ая всеросийская конференция с международным участием «масс-спектрометрия и ее прикладные проблемы», 03-07 сентября, 2007 г., г. Москва.

21.Павел Рюмин, Иван Болдин, Дмитрий Автономов, Александр Позднеев, Евгений Николаев. Сравнительный анализ результатов моделирования движения ионных ансамблей в ИЦР ловушке кубической геометрии методом Р1С(частица в ячейке) и методом емкостей. 3-ий съезд ВМСО, 2-ая всеросийская конференция с международным участием «масе-спектрометрия и ее прикладные проблемы», 03-07 сентября, 2007 г., г. Москва.

22.Dmitry Avtonomov, Pavel Ryumin, Ivan Boldin, Alexander Pozdneev, Eugene Nikolaev, Capacitance matrix algorithm implementation for supercomputer simulations of ion dynamics in Kingdon trap. 8th European FTMS Conference, Moscow, Russia, August 28 - September 1, 2007.

23.Eugene N Nikolaev, Ron Heeren, Alexander Pozdneev, Gleb Vladimirov, Ivan Boldin, Pavel Ryumin, Dmitry Avtonomov The new possibilities in ion clouds dynamic simulation using supercomputers. Application to FTICR,

Kingdon trap and accumulation quadrupole devices. 55th American Society for Mass Spectrometry Conference, Indianapolis, IN, USA, June 3-7, 2007.

24.Болдин И.А., Рюмин П.А., Автономов Д.М., Позднеев A.B., Николаев E.H. Метод ёмкостей в моделировании движения ионных ансамблей в ионных ловушках и системах транспорта ионов с электродами произвольной формы. Труды 50-ой научной конференции МФТИ, Москва - Долгопрудный, 2007.

Подписано в печать:

03.11.2011

Заказ №6195 Тираж - 50 экз. Печать трафаретная. Типография «11-й ФОРМАТ» ИНН 7726330900 115230, Москва, Варшавское ш., 36 (499) 788-78-56 \vww.autoreferat. ги

 
Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Рюмин, Павел Александрович

Введение.

Глава 1. Движение заряженных частиц в ионных ловушках Пеннинга и Кингдона. Влияние объемного заряда. Методы моделирования движения ионов в масс-спектрометрических устройствах.

Глава 2. Моделирование движения ионов в масс-спектрометрических устройствах с электродами произвольной формы. Учет взаимодействия ионного облака со стенками электродов.

Глава 3. Симметрия задачи. Улучшенный алгоритм, позволяющий уменьшить вычислительную сложность, используя симметрию задачи.

Глава 4. Моделирование влияния объемного заряда на движение ионных облаков различных размеров в ловушке Кингдона.

Глава 5. Аксиальное детектирование в ловушке Пеннинга с динамической гармонизацией.

Выводы.

 
Введение диссертация по физике, на тему "Суперкомпьютерное исследование движения ионов в ловушках Кингдона и Пеннинга с полным учетом кулоновского взаимодействия"

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

Наилучшими аналитическими характеристиками среди всех видов масс-спектрометров, такими как разрешающая способность и точность определения масс, обладают масс-спектрометры с преобразованием Фурье. К таким масс-спектрометрам относятся приборы, основанные на принципе ионного циклотронного резонанса (ИЦР) и на принципе орбитальной ловушки Кингдона (Орбитрэп). В масс-спектрометрах такого типа ионы запираются в определенной области пространства электромагнитным или электрическим полем, в котором совершают периодические движения. Для, измерения отношения массы к заряду регистрируется частота этих движений по наведенному на электродах ловушки заряду изображения. После регистрации сигнал подвергается частотному анализу, чаще всего для этого используют преобразование Фурье. Разрешающая способность при таком методе измерения зависит от длительности измеряемого сигнала. В общем случае лимитирующими факторами для получения - длительного сигнала являются недостаточный вакуум, а также неидеальности электрического и магнитного поля. В экспериментах с большим количеством частиц на движение ионов также начинает влиять кулоновское взаимодействие и взаимодействие с зарядом, наводимым ионами на стенках электродов. Получаемые в результате частотного анализа спектры частот преобразуются 3 в спектры масс по известному соотношению, связывающему частоту с отношением массы к заряду ионов. Большое количество ионов в ионных ловушках, например, типа Пеннинга (ИЦР ПФ) приводит к таким явлениям как коалесценция пиков в масс-спектрах и разрушение ионных облаков из-за их взаимодействия, что не позволяет достигать высокую разрешающую способность одновременно с достижением большого динамического диапазона. Показано, что при наличии в ловушке ионов с, разным отношением массы к заряду кулоновское взаимодействие между образуемыми ими облаками приводит к частотным сдвигам, которые ограничивают точность измерения масс. Даже в том случае, когда в ловушке изолируются только ионы с одним отношением массы к заряду, при значительном радиусе возбуждения взаимодействие со стенками электродов так же приводит к сдвигу циклотронной частоты [1-5]. Подобные эффекты наблюдаются также и в ловушке Кингдона [6].

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

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

Появление коммерческой программы БГМКЖ [7-8] значительно упростило понимание различных аспектов движения единичных ионов в ионных ловушках произвольной геометрии. Эта программа способна рассчитывать траектории отдельных частиц в различных измерительных ячейках ИЦР с довольно высокой точностью. Однако программа не предоставляет достаточно адекватных возможностей учитывать кулоновское взаимодействие между ионами.

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

В физике плазмы был разработан подход «частица в ячейке»(Р1С), в котором время вычислений зависит линейно от количества частиц. С использованием этого метода несколькими группами было промоделировано движение ионных ансамблей в кубической ловушке Пеннинга,' использовавшейся в масс-спектрометрии ИЦР. Основным недостатком метода «частица в ячейке» является отсутствие возможности моделирования движения ионов в ловушках с непрямоугольной геометрией электродов, как в большинстве реальных устройств.

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

Цель работы

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

Научная новизна работы

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

Практическая значимость работы

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

Структура диссертации

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

Во второй главе излагается метод емкостной матрицы, а так же его ограничения. Также дается описание алгоритма поверхностных зарядов, аналогичного емкостному методу, который позволяет более точно вычислять электрическое поле от электродов ловушки. Приводятся рассчитанные методом поверхностных зарядов поля. Оценивается точность вычисления полей, приводятся результаты сравнения полей рассчитанных методом поверхностных зарядов и коммерческой программой БГШО^Т. Приводятся результаты сравнения разработанного подхода для моделирования движения ионов в масс-спектрометрических устройствах с методом частица в ячейке.

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

В четвертой главе представлены результаты исследования влияния кулоновского взаимодействия на движение ионов в ловушке Кингдона.

В пятой главеописываются возможности аксиального детектирования в ячейке Пеннинга с динамической гармонизацией. Исследуется зависимость критической массы от количества сегментов ячейки с динамической гармонизацией.

 
Заключение диссертации по теме "Химическая физика, в том числе физика горения и взрыва"

Заключение. Выводы.

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

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

3) Исследовано влияние объемного заряда на движение ионов в ловушке Кингдона.

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

5) Выявлены особенности движения больших ионных ансамблей в ловушке Кингдона.

 
Список источников диссертации и автореферата по физике, кандидата физико-математических наук, Рюмин, Павел Александрович, Москва

1. Цитируемая литература

2. Wineland, D.; Dehmelt, Н. Line Shifts and Widths of Axial, Cyclotron and G-2 Resonances in Tailored, Stored Electron (Ion) Cloud. Int. J. Mass Spectrom. Ion Phys. 1975,16, 338-342.

3. Wong, R. L.; Amster, I. J. Experimental Evidence for Space-charge Effects Between Ions of the Same Mass-to-Charge in Fourier-Transform Ion Cyclotron Resonance Mass Spectrometry. Int. J. Mass Spectrom. 2007, 265, 99-105.

4. X. Xiang; P. B. Grosshans; A. G. Marshall. Image charge-induced ion cyclotron orbital frequency shift for orthorhombic and cylindrical FT-ICR ion traps. Int. J. Mass Spectrom. Ion Processes, 1993, 125, 33.

5. Makarov, A.; Denisov, E.; Lange, O.; Horning, S. Dynamic Range of Mass Accuracy in LTQ Orbitrap Hybrid Mass Spectrometer. Journal of the American Society for Mass Spectrometry 2006, 17, 977-982.

6. Dahl, D. A. Simion 3D Version 6.0 User's Manual, Idaho Falls, ID: Idaho National Engineering Laboratory Chemical Materials and Processes Department, Lockheed Idaho Technologies Company, 1995.

7. Dahl, D. A. SIMION for the personal computer in reflection, J. Mass Spectrom., 2000, 200, 3-25.

8. Penning, F.M., Philips techn. Rev. 2 1937, 201.

9. O.Marshall, A. G.; Hendrickson, C. L.; Jackson, G. S. Fourier transform ion cyclotron resonance mass spectrometry: a primer. Mass Spectrom. Rev., 1998, 17, 1-35.

10. L. S. Brown; G. Gabrielse, Geonium theory: Physics of a single electron or ion in a Penning trap, Reviews of Modern Physics, 1986, 58, 233

11. E.O. Lawrence; M.S. Livingston, Phys. Rev., 1932, 40, 19.

12. Hippie H. A., Sommer H., Thomas J. A., A Precise Method of Determining the Faraday by Magnetic Resonance. Phys. Rev., 1949, 76, 1877-1878.

13. Sommer H, Hippie H A and Thomas J A. The Measurement of e/M by Cyclotron Resonance. Phys. Rev. 1951, 82, 697.

14. M.B. Comisarow; A.G. Marshall. Fourier Transform Ion Cyclotron Resonance Spectroscopy. Chem. Phys. Lett., 1974, 25, 282.

15. Comisarow MB; Marshall AG. Frequency-sweep Fourier transform ion cyclotron resonance spectroscopy. Chem. Phys. Lett. 1974; 26: 489.

16. Huang, J.; Tiedemann, P. W.; Land, D. P.; Mclver, R T.; Hemminger, J. C. Dynamics of ion coupling in an FTMS ion trap and resulting effects on mass spectra, including isotope ratios. Int. J. Mass Spectrom. Ion Processes 1994, 134, 11.

17. Naito Y.; Inoue M. Mass Spectrometry Using Fourier Transform . J. Mass Spectrom. Soc. Jpn., 1994, 42, 291.

18. R. S. Van Dyck Jr.; D. J. Wineland; P. A. Ekstrom; H. G. Dehmelt. High mass resolution with a new variable anharmonicity Penning trap. Applied Physics Letters, 1976, 28, 446.

19. Van Dyck, R. S.; Jr., Schwinberg, P. B. Preliminary Proton/Electron Mass Ratio using a Compensated Quadring Penning Trap. Phys. Rev. Lett. 1981, 47, 395.

20. Yin, W. W.; Wang, M.; Marshall, A. G.; Ledford, E. B. Experimental evaluation of a hyperbolic ion trap for Fourier transform ion cyclotron resonance mass spectrometry. J. Am. Soc. Mass Spectrom., 1992, 3, 188.

21. Precise measurements of axial, magnetron, cyclotron, and spin-cyclotron-beat frequencies on an isolated 1 meV electron R.S. Van Dyck Jr.; P.B. Schwinberg; H.G. Dehmelt, Physical Review Letters, 1977, 14 Feb,310-14

22. L. Schweikhard; M. Blundschling; R. Jertz; H.-J. Kluge, Fourier transform mass spectrometry without ion cyclotron resonance: direct observation of the trapping frequency of trapped ions Int. J. Mass Spectrom. Ion Processes, 1989, 89, R7-12

23. Kingdon, A method for the Neutralization of Electron Space Charge by Positive Ionization at very low gas pressures, Phys. Rev. 1923,21, 408418

24. Sekioka, T.; Terasama, M.; Awaya, Ion storage in Kingdon trap,Y. Rad. Effects Defects Solids 1991,117, 253-259.

25. Richard H. Perry; R. Graham Cooks; Robert J. Noll, Orbitrap mass spectrometry: instrumentation, ion motion and applications, Mass Spectrometry Reviews, 2008, 27, 661- 699

26. Gillig, K. J.; Bluhm, B. K.; Russell, D. H. Ion motion in a Fourier transform ion cyclotron resonance wire ion guide cell, Int. J. Mass Spectrom. Ion Processes 1996,157/158, 129-147

27. Knight, R. D., Storage of ions from laser-produced plasmas, Appl. Phys. Lett. 1981, 38, 221-222.

28. Mitchell D. W.; Smith R. D. Cyclotron motion of two Coulombically interacting ion clouds with implications to Fourier-transform ion cyclotron resonance mass spectrometry. Phys. Rev. E 1995, 52, 4366.

29. X. Xiang; P. B. Grosshans; A. G. Marshall. Image charge-induced ion cyclotron orbital frequency shift for orthorhombic and cylindrical FT-ICR ion traps. Int. J. Mass Spectrom. Ion Processes, 1993, 125, 33.

30. E. N. Nikolaev; R. M. A. Heeren; A. M. Popov; A. V. Pozdneev; K. S. Chingin. 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, 21, 3527.

31. Mitchell D. W., Smith R. D. Two dimensional many particle simulation of trapped ions. Int. J. Mass Spectrom. Ion Processes 1997, 165/166, 271.

32. Cornell, E. A; Boyce, K. R; Fygenson, D. L.; Pritchard, D. E. Two ions in a Penning trap: Implications for precision mass spectroscopy. Phys. Rev. A 1992, 45, 3049.

33. Huang, J.; Tiedemann, P. W.; Land, D. P.; Mclver, R T.; Hemminger, J. C. Dynamics of ion coupling in an FTMS ion trap and resulting effects on mass spectra, including isotope ratios. Int. J. Mass Spectrom. Ion Processes 1994, 134, 11.

34. Peurrung, A. J.; Kouzes, R. T. Long-term coherence of the cyclotron mode in a trapped ion cloud. Phys. Rev. E1994, 49, 4362.

35. Ландау Jl. Д., Лифшиц Е. М. Теоретическая физика. Т. II. Теория поля. 8-е изд. М.: ФИЗМАТЛИТ, 2001.

36. Pasa-Tolic L.; Huang Y.; Guan S.; Kim H. S., Marhsall A. G. Ultrahighresolution matrix-assisted laser desorption/ionization Fourier transform ion cyclotron resonance mass spectra of peptides. J. Mass Spectrom. 1995, 30, 825.

37. Naito Y.; Inoue M. Mass Spectrometry Using Fourier Transform . J. Mass Spectrom. Soc. Jpn. 1994, 42, 291.

38. Naito Y., Inoue M. Collective motion of ions in an ion trap for Fourier transform ion cyclotron resonance mass spectrometry. Int. J. Mass Spectrom. Ion Processes 1996; 157/158: 85.

39. A. Boldin; E. N. Nikolaev. Theory of peak coalescence in Fourier transform ion cyclotron resonance mass spectrometry. Rapid Commun. Mass Spectrom. 2009, 23, 3213.

40. Jeffries, J. В.; Barlow, S. E.; Dunn, G. H., Theory of space-charge shift of ion cyclotron resonance frequencies, Int. J. Mass Spectrom. Ion Proc.1983, 54, 169-187.

41. Ledford, E. В.; Rempel, D. E.; Gross, M. L., Space-charge effects in Fourier transform mass spectrometry mass calibration, Anal. Chem.1984, 56, 2744-2748.

42. Francl, T. J.; Sherman, M. G.; Hunter, R. L.; Locke, M. J.; Bowers, W. D.; Mclver, R.T., Experimental determination of the effects of spacecharge on ion-cyclotron resonance frequencies, Int. J. Mass Spectrom. IonProc. 1983, 54, 189-199.

43. Easterling, M. L.; Mize, T. H.; Amster, I. J., Routine Part-per-Million Mass Accuracy for High-Mass Ions: Space-Charge Effects in MALDI FT-ICR, ¿wa/. Chem. 1999, 71, 624-632.

44. Masselon, C.; Tolmachev, A. V.; Anderson, G. A.; Harkewicz, R.; Smith, R. D., Mass measurement errors caused by "local" frequency perturbations in FTICR mass spectrometry, J. Am. Soc. Mass Spectrom. 2002, 13,99-106.

45. Muddiman, D. C.; Oberg, A. L., Statistical evaluation of internal and external mass calibration laws utilized in Fourier transform ion cyclotron resonance mass spectrometry, Anal. Chem. 2005, 77, 2406-2414.

46. Wong, R. L; Amster, I. J., Sub part-per-million mass accuracy by using stepwiseexternal calibration in Fourier transform ion cyclotron resonance mass spectrometry, J. Am. Soc. Mass Spectrom. 2006, 17, 1681-1691.

47. Miluchihin N. V., Miura K., Inoue M. Application of a parallel computer to simulation of ion trajectories in an ion cyclotron resonance spectrometer. Rapid Commun. Mass Spectrom. 1993, 7, 966.

48. Birdsall CK, Langdon AB. Plasma Physics via Computer Simulation. McGraw-Hill, New York, 1985.

49. Hockney RW, Eastwood JW. Computer Simulation Using Particles. Adam Hilger: New York, 1988.

50. Birdsall CK. Particle-in-cell charged-particle simulations, plus Monte Carlo collisions with neutral atoms, PIC-MCC. IEEE. Transactions on Plasma Science 1991, 19, 65.

51. Cooley, J. W.; Tukey, J. W. An algorithm for the machine computation of the complex Fourier series. Math. Computation, 1965, 19, 297 301.

52. Frigo, M.; Johnson, S. J. "The design and implementation of FFTW3." Proc. IEEE, 2005, 93, 216-231.

53. E. L. Allgower; K. Böhmer; K. Georg ; R. Miranda, Exploiting symmetry in boundary element methods, SIAMJ. Numer. Anal., 1992, 29, pp. 534552.

54. Eugene L, Allgower; Kurt Georg, Exploiting Symmetry in BEM, doi 10.1.1.47.6679,1996

55. K. Ahlander ; H. Munthe-Kaas, Applications of the generalized Fourier transform in numerical linear algebra, BIT Numerical Mathematics , 2005 45: 819-850

56. H Z Munthe-Kaas; On group Fourier analysis and symmetry preserving discretizations ofPDEs, J. Phys. A: Math. Gen. 2006, 39, 5563-5584

57. Alexander Makarov, Eduard Denisov, and Oliver Lange, Performance Evaluation of a High-field Orbitrap Mass Analyzer, Journal of The American Society for Mass Spectrometry, 2009 20, 8

58. H. F. Weinberger, A first course in partial differential equations with complex variables and transform method, Dover Publications,199571.0ртега Дж. Введение в параллельные и векторные методы решения линейных систем: Пер. с англ.-М.: Мир, 1991

59. Серр Ж.-П. Линейные представления конечных групп, М.: Мир, 1970