Модель нелинейного дрейфа ионов в спектрометрии приращения ионной подвижности тема автореферата и диссертации по физике, 01.04.01 ВАК РФ

ИНСТИТУТ ЭНЕРГЕТИЧЕСКИХ ПРОБЛЕМ ХИМИЧЕСКОИ ФИЗИКИ РАН

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

ООЗОБ4ЭЭЭ

Шибков Сергей Викторович

МОДЕЛЬ НЕЛИНЕЙНОГО ДРЕЙФА ИОНОВ В СПЕКТРОМЕТРИИ ПРИРАЩЕНИЯ ИОННОЙ ПОДВИЖНОСТИ

01.04.01 - приборы и методы экспериментальной физики

Автореферат

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

Москва - 2007

003054999

Работа выполнена в Институте криптографии, связи и информатики Академии ФСБ России

Научный руководитель: доктор физико-математических наук,

С. Д. Бенеславский

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

В. В. Зеленов

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

Ведущая организация: Институт аналитического приборостроения РАН

Защита состоится « 17» января 2007 г. в 12 часов 00 минут на заседании диссертационного совета Д 002.112.01 при Институте энергетических проблем химической физики РАН по адресу: 119334, Москва, В-334, Ленинский проспект, д. 38, корп. 2.

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

2006 г.

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

диссертационного совета Д 002.112.01, кандидат химических наук

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность темы

Впервые метод разделения ионов по виду зависимости их подвижности от напряженности электрического поля при атмосферном давлении был реализован Горшковым М.П. [1]. Экспериментальная установка представляла собой спектрометр приращения ионной подвижности с плоской дрейф-камерой, в котором источником ионизирующего излучения служил радиоактивный изотоп 63Ni, а в качестве системы регистрации ионов использовался электрометрический усилитель. Дальнейшее развитие данный способ разделения ионов получил в работе [2], где впервые был предложен метод идентификации ионов третичных аминов, определение которых является сложной задачей из-за высокой реакционной и адсорбционной способности. Для ионизации использовался /^-источник 3Н, а также молибденовый термоэмиттер для поверхностной ионизации. В работах [3] и [4] с помощью присоединенного масс-спектрометра был исследован состав газовой смеси на выходе дрейф-камеры спектрометра приращения ионной подвижности (СПИП-спектрометра), что позволило рассмотреть ионные процессы, сопровождающие процесс дрейфа ионов.

Для снижения относительного числа ионов, рекомбинирующих на обкладках плоской дрейф-камеры, были разработаны спектрометры с цилиндрической дрейф-камерой. Впервые прибор с такой конфигурацией электродов был создан Буряковым И.А., Крыловым Е.В., Солдатовым В.П. [5], а затем Карнаганом Б.Л. и Тарассовым A.C. [6]. В работе [7], а также независимо [8] было экспериментально исследовано явление фокусировки ионов, которое позволяло существенно снизить потери ионов на обкладках.

В дальнейшем на основе СПИП-спектрометров в России и за рубежом был создан ряд портативных приборов для обнаружения следовых количеств различных веществ: для экологического мониторинга [9], [10], для решения

поисковых задач [11], [12], для экспрессного и лабораторного медицинского анализа [13], [14].

В работах [15], [16] исследовались чувствительность, разрешающая способность СПИП-спектрометров, зависимость спектрометрических данных от расхода и температуры газа-носителя. Одновременно расширялся спектр веществ, которые обнаруживались спектрометрами приращения ионной подвижности [17]. Появлялись тандемные устройства [18], [19], [20], [21], в составе которых СПИП-спектрометры использовались в качестве устройства концентрации пробы, ионного фильтра на входе масс-спектрометров различных типов, а также как детекторы на выходе скоростной хроматографической установки или источника ионизации в виде электроспрея. Исследовались различные источники ионов для СПИП-спектрометров: электроспрей [22], коронный разряд, поверхностная ионизация [23], лазерная ионизация [24], высокочастотный емкостной разряд [25].

Рассматривались такие физические явления, как захват ионов в дрейф-камере спектрометров приращения ионной подвижности [26], влияние влажности на данные спектрометрии [27]. В работах [28] и [29] обсуждались способы уменьшения диффузионных потерь, проводился сравнительный анализ спектрометров с плоской и цилиндрической дрейф-камерой. Некоторые вопросы теории СПИП были рассмотрены в работах [30], [31], [32], [33], [34].

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

Цель работы

Целью работы явилось создание модели нелинейного дрейфа ионов в СПИП-спектрометрах, имеющих плоскую и цилиндрическую дрейф-камеры и

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

Научная новизна

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

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

3. Для описания ансамбля ионов в дрейф-камере СПИП-спектрометра применены такие понятия теории динамических систем как отображение Пуанкаре, фазовый портрет, предельный цикл, фазовое многообразие, аттрактор, положение равновесия, бифуркация рождения пары предельных циклов и связанная с ней перестройка фазового портрета динамической системы.

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

1. В работе проведено теоретическое рассмотрение процесса разделения ионов в сильных пространственно однородных и неоднородных электрических полях при атмосферном давлении. Рассчитаны основные величины, наблюдаемые в экспериментальных установках,

использующих данный принцип разделения ионов - СПИП-спектрометрах с различной геометрией дрейф-камеры.

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

Основные положения, выносимые на защиту

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

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

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

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

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

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

Апробация работы

Результаты диссертационной работы докладывались и обсуждались на семинаре Секции Прикладных проблем при Президиуме РАН (Москва, 2002 г.), на XXIX Звенигородской конференции по физике плазмы и управляемому термоядерному синтезу (Звенигород, 2002 г.), на Межведомственных конференциях по проблемам технического обеспечения безопасности (Москва, 2002 г. и 2004 г.), на Международной конференции по масс-спектрометрии (1МБС-16, Эдинбург, 2003 г.), на Втором международном семинаре-школе «Масс-спектрометрия в химической физике, биофизике и экологии» (Звенигород, 2004 г.), на Всероссийской конференции «Масс-спектрометрия и ее прикладные проблемы», проводившейся Всероссийским масс-спектрометрическим обществом (Москва, 2005 г.).

Публикации

Основные результаты работы изложены в 4 публикациях, две из которых опубликованы в отечественном реферируемом журнале и две в международном реферируемом журнале. Список публикаций приведен в конце автореферата. Работа была выполнена в Институте криптографии, связи и информатики Академии ФСБ России.

Объем и структура диссертации

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

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В разделе 1.1 первой главы формулируется модель динамики ионов в СПИП-спектрометрах с плоской дрейф-камерой.

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

*(£) = *0(1+«(£)),

(1)

где к о — подвижность в слабых полях, се(Е) - переменная составляющая подвижности. Разделение ионов в спектрометрах приращения ионной

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

Процесс разделения ионов в СПИП-спектрометре (Рис. 1) происходит, например [11], следующим образом: проба воздуха, содержащая примеси анализируемых веществ, попадает в камеру ионизации 1; под действием излучения /^-источника 2, молекулы примеси ионизируются; в области газового

5ашр!е _/ ^"у

гз; О О 3

у* С

Т"

--

А Лл / Л В

Ш 5 и£) = и0АЪ

т,«тс

\ис 6 ___ 1

0 Тс

Рис. 1 Схема СПИП-спектрометра.

затвора и системы транспортировки 3, ионы, двигаясь под действием продольного электрического поля навстречу запирающему газовому потоку, попадают из ионизационной в дрейф-камеру 4; к коаксиальным цилиндрическим обкладкам дрейф-камеры, приложено большое по амплитуде высокочастотное разделяющее напряжение иБ и малое медленно меняющееся компенсирующее напряжение ис, генерируемые источниками 5 и 6 соответственно; ионы осциллируют в радиальном и увлекаются газом-носителем в аксиальном направлении; ионы, для которых выполнено условие отбора, сформулированное ниже, достигают выхода дрейф-камеры и коллектора ионов 8; для каждого значения С/с, ток прошедших ионов регистрируется электрометром 7, на выходе которого возникает ионограмма

7(С/с); газ-носитель внутреннего контура, разбавляющий малый входной поток, подвергается очистке с помощью фильтра 9.

Анализ условий дрейфа ионов в типичных СПИП-спектрометрах показывает, что исходными физическими положениями модели могут быть приняты следующие: 1) движение ионов в газовой среде происходит в сильном электрическом поле; 2) газ, содержащий ионы, может рассматриваться в пренебрежении эффектами сжимаемости, то есть в гидродинамическом приближении; 3) течение газа имеет ламинарный характер с однородным вдоль всей дрейф-камеры профилем скоростей; 4) рассматриваются ионные токи, при которых влиянием пространственного заряда на электрическое поле внутри дрейф-камеры и движение ионов можно пренебречь; 5) не учитывается тепловая диффузия ионов. В рамках этих предположений, описание поведения ансамбля ионов, даваемое уравнением конвекции-диффузии, посредством метода характеристик сводится к задаче движения отдельного иона в вязкой среде. Для решения этой задачи применяется стробоскопический метод усреднения быстрых осцилляций. С помощью построенной модели рассчитывается вид ионного пика на выходе спектрометра с плоской дрейф-камерой. В пренебрежении диффузией и влиянием объемного заряда пик имеет треугольную форму.

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

При анализе газовой смеси, содержащей примесь определенного вещества, местоположение максимума ионного тока этого вещества 1/с, зависит от амплитуды приложенного разделяющего напряжения С/5. Проведя измерения для набора значений ¡У5, получим набор ионограмм, аналогичный изображенному на Рис. 2(а). Используя один из алгоритмов интерполяции, строим график зависимости ис(и$) (Рис. 2(6)), проходящий через

экспериментально полученные точки С/од)- Записав условие отбора ионов для плоской дрейф-камеры - радиальное смещение за период С/у равно нулю, -получим обобщенное уравнение, связывающее зависимость компенсирующего напряжения, соответствующего положению ионного пика, от амплитуды разделяющего напряжения ис(и$) и а(Е)

МС{Е)

(2)

где (см. Рис. 1, зависимость £/$(/)) М - участки постоянства функции Х0> определяющей временную зависимость разделяющего напряжения; -значения напряженности на участках <Г2у - участки взаимно однозначного

А М Иф

1,„.

л

/ ист \

ист Чт

(Л'<х Л

и с (А) Л

С/с, Г

-I -2

-3

500~\1000 1500

10

20

30

Рнс. 2 Ионограммы, зависимость 11с{Ц5), аналитическая зависимость а(Е) и последовательные итерации для уравнения Вольтерра (указан номер итерации).

соответствия переменных ? и £', на которых существуют обратные функции (=<р)(Е'). Решение ЫЕ) уравнения ( 2) дает искомую переменную составляющую подвижности ионов.

В работе показано, что если функция_/(?) обладает участками постоянного значения, аналогичными ДГ| и на Рис. 1, то уравнение ( 2) сводится к интегральному уравнению Вольтерра II рода. Напротив, широко используемая в спектрометрах зависимость Дг) в виде суперпозиции двух косинусов, не принадлежит к вышеуказанному классу функций, в силу чего не существует

однозначной связи ис (С/?) и а(Е) и, следовательно, данный вид разделяющего напряжения не может использоваться для восстановления зависимости подвижности от напряженности поля.

Далее в работе был проведен численный эксперимент по восстановлению зависимости подвижности от напряженности поля с помощью уравнения ( 2). Задавшись аналитическим видом а(Е), изображенным сплошной кривой на Рис. 2(е), и временной зависимостью /(7), относящейся к указанному выше типу функций, аналогичной изображенной на Рис. 1, но с отфильтрованными высокочастотными осцилляциями, с помощью выражений для ионного пика из раздела 1.1 Главы 1, получим серию пиков (Рис. 2(а)), затем график ис(и5) (Рис. 2(6)). Используя (С/у) в правой части ( 2) и упомянутую выше временную зависимость_Д/) для нахождения (рЕ') в левой части, получим из ( 2) уравнение Вольтерра II рода, результаты решения которого методом последовательных приближений приведены на Рис. 2(в). Благодаря свойству сжимаемости интегрального оператора Вольтерра II рода, любое начальное приближение демонстрирует сходимость к искомой функции.

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

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

зависимости подвижности от напряженности поля является важной общефизической и прикладной задачей.

Во второй главе формулируется модель динамики ионов в СПИП-спектрометрах с цилиндрической дрейф-камерой.

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

± = *(£(г,0)Е(гД

т

(3)

где г - радиальная координата иона, к(Е) определено согласно ( 1), Е{г,?) -радиальная компонента вектора напряженности электрического поля в дрейф-камере. Уравнение ( 3) описывает движение отдельного иона в вязкой среде. Это нелинейное дифференциальное уравнение является неавтономной динамической системой, для интегрирования которой не существует общих методов. Ионы, находясь в дрейф-камере, под действием разделяющего поля совершают несколько тысяч осцилляций. Поэтому, несмотря на малость нелинейной составляющей подвижности а(Е), интегрирование уравнения ( 3) приводит к относительной погрешности, сравнимой с единицей. Эта проблема является разновидностью известной в физике задачи учета влияния малых возмущений на больших временах, которая, по высказыванию А. Пуанкаре, является основной задачей динамики. Для устранения этой трудности в диссертации анализ системы (3), был проведен с помощью стробоскопического метода усреднения быстрых осцилляций.

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

Исходная динамическая система ( 3) обладает цилиндрическим фазовым многообразием (см. Рис. 3), где угол <р = (2я//Т)шос1(271) - фаза разделяющего напряжения, а координата х - положение частицы в зазоре дрейф-камеры.

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

Центральным понятием стробоскопического метода усреднения и основной характеристикой динамической системы является функция смещения ионов Е(х0,иа1) = х(х0,иС0)-х,) за период разделяющего напряжения Т, которая

определяется с помощью интегрирования уравнения ( 3) на отрезке [0;7] (здесь

использованы безразмерные величины: начальная координата частицы Х(, и известная из теории динамических систем функция последования или отображение Пуанкаре х(х0,исо), значение которого определяет координату в конце периода как функцию начальной координаты). Вид функции смещения при некотором значении С/со изображен на Рис. 4(е). Интегрирование на отрезке [0;7], например численное, может быть проведено с гораздо более высокой точностью, чем в течение всего времени пребывания ионов в дрейф-камере. Приведенное выше определение функции смещения может быть использовано в качестве разностной аппроксимации дифференциального уравнения

Рис. 3 Цилиндрическое фазовое многообразие

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

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

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

4 и, и2 и,-

Рис. 4 Построение ионного пика, функция смешения, усреднение.

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

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

Траектории ионов в приближении постоянной подвижности k0, то есть в случае а(Е) = 0, будут периодическими при Uc = 0. Назовем множество траекторий x(t), по которым ионы могут двигаться, не сталкиваясь с обкладками дрейф-камеры, областью выживания. Она ограничена, очевидно, двумя периодическими траекториями xSD\(t) и xSca(t), касающимися внутренней и внешней обкладки дрейф-камеры. Подчеркнем, что область выживания есть область в фазовом пространстве соответствующей динамической системы (см. Рис. 3). Назовем окном выживания отрезок радиальных координат [xSm< xsm]> заключенный между двумя траекториями, ограничивающими область выживания, в начале каждого периода разделяющего поля (ср = 0).

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

Если у - корень функции F, то есть F(y',UCo) = 0, то, у(<%) =у* = const является решением уравнения ( 4), называемым положением равновесия. Система ( 3) при этом имеет периодическое решение x\t) = x*{t+T), называемое предельным циклом.

Введем параметр т](иСо) = оЕ(у*,исаУ8у, который назовем интенсивностью фокусировки, знак которой определяет устойчивость положения равновесия, а величина характеризует поведение траекторий в окрестности предельного цикла (прямая на Рис. 4(е)).

Назовем траекторией предельного цикла в пространстве параметров динамической системы зависимость координаты положения равновесия усредненной системы от компенсирующего напряжения у\исо) (Рис. 4(6)), которая определяется как неявная функция уравнением Е(у'(исо),исо) ~ 0.

Рассмотрим процедуру построения ионного пика на конкретном примере. Примем зависимость от времени разделяющего напряжения кусочно-постоянной Щг) = ¡У0./(г), где функция/(г) = {1, 0 < т< Т\, -Т\1{\-Т\), Т\ < т< 1}, т= 1/Т - безразмерное время, Г1 - момент переключения полярности напряжения. Зависимость подвижности от величины напряженности электрического поля предположим полиномиальной простейшего вида к(Е) = ко(\ + 0(Е/Е\)2), Р — малый безразмерный коэффициент, характеризующий относительную величину отклонения подвижности от постоянного значения к0 в пределе слабого поля, Е\ - величина напряженности поля, при которой подвижность существенно отличается от кй. Функция смещения, как показано в работе автора [2], имеет вид

Р(у,исо) = сисо+/31в 1п

V

1 + ^

У

(5)

где /?, =/3{Е(/Е,)\ £0 = £Мг21п(г2/г,)), с=2къЕ0Т/г2 и С/со= ис/и0, в= (1-2г,)/ (1-Г1)2, г ь г2 - радиус внутренней и внешней обкладки дрейф-камеры, у = (г/г2)2 -безразмерная координата иона. Используем функцию смещения ( 5) в уравнении ( 4). Линеаризуя уравнение ( 4) в окрестности положения равновесия у{исо), а затем интегрируя полученное линейное уравнение, получим траектории усредненной системы

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

У*Фсо)= Т~~Т, \ ' =

Ехр

сУс о рхв

-1

2Д0

(7)

Интегральные кривые ( 6) линеаризованного уравнения ( 4) с правой частью ( 5) в случае устойчивого положения равновесия (т](исо) < 0) изображены на Рис. 4(а) и (е).

На Рис. 4(а) положение равновесия у (£/3) находится внутри окна выживания. Траектории стремятся к / и все частицы, попавшие в окно выживания, дают вклад в ионный ток на выходе. Как видно из ( 6) и Рис. 4(а) ширина Ь ионного шнура ограниченного жирными линиями экспоненциально убывает: ¿(£) = Ь0ехр{т!(исо)^}, где Ь0 - ширина ионного шнура на входе дрейф-камеры и г/(иСо)< 0. Спектрометр работает при этом в режиме фокусировки. Очевидно, что при любом значении компенсирующего напряжения и со из интервала [и^и2] ток на выходе спектрометра остается постоянным (Рис. 4(д)), так как все частицы, попавшие в окно выживания на входе дрейф-камеры, достигают выхода. Таким образом, режим фокусировки характеризуется наличием платообразного участка ионного пика.

На Рис. 4(в) положение равновесия >■ (С/4) находится вне окна выживания. Траектории по-прежнему стремятся к у", но часть частиц достигают границы окна выживания и погибают на обкладке дрейф-камеры в результате рекомбинации. В результате при удалении положения равновесия от окна выживания все большая доля частиц не достигает выхода дрейф-камеры. Величина ионного тока убывает по сравнению со значением, соответствующим

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

Аналогичным образом может быть получен вид пика в случае положительной интенсивности фокусировки (rj(UCo) > 0), то есть режима дефокусировки. При этом из ( 6) видно, что траектории ионов экспоненциально удаляются от неустойчивого положения равновесия. Вследствие этого ток на выходе спектрометра уменьшается с увеличением времени нахождения ионов в дрейф-камере.

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

Интересно отметить, что, например, в случае 0<ri<0.5 для отрицательных ионов (с<0) явление фокусировки (устойчивый предельный цикл) наблюдается только при сублинейном < 0) характере подвижности, пик на ионограмме при этом находится в области положительных компенсирующих напряжений ( Uco>0)-

Как видно из ( 7), в случае с < 0, 0< ri <0.5, р t <0 при UCx> —» °° х0 —> - стх и T](Uco)<0, таким образом при - сг, е [х\, х2 ] будет происходить «зависание» устойчивого предельного цикла в окне выживания. При этом на эксперименте должно наблюдаться практически равномерное засвечивание всего дрейф-спектра.

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

В разделе 3.1. третьей главы проводится исследование динамической системы, описывающей поведения ансамбля ионов, обладающих немонотонной зависимостью подвижности от напряженности поля. В этом случае функция смещения F(y, Ucо) может обладать немонотонной зависимостью от начальной координаты у (см. Рис. 5(а)).

Рассмотрим функцию смещения F(y,Uçx>) как функцию одной переменной F(y), зависящую от Uco как от параметра. Уменьшая значение параметра С/со мы наблюдаем, к примеру, следующие характерные случаи: С/со > Ub — положения равновесия отсутствуют, вид функции смещения F(y) представлен на Рис. 5(а), кривая 2; UCo < Uh - существует два корня функции смещения F(y) (Рис. 5(а), кривая 1), соответствующие неустойчивому yt * и устойчивому^* положению равновесия динамической системы. При прохождении компенсирующего напряжения Uco значения U/, происходит явление называемое бифуркацией рождения пары предельных циклов. Траектория предельного цикла в пространстве параметров y'{Ucо) в отличие от случая монотонной зависимости к(Е) (соответствующая траектория изображена пунктирной линией на Рис. 5(в)) становится двузначной функцией (сплошная кривая на Рис. 5(е)). Такая функция смещения может быть аппроксимирована в окрестности максимума (точка Ъ) параболой (линия, изображенная точками на Рис. 5(е)), что приводит к уравнению ( 4) следующего вида

dy

-— = c{U-Ul) + a(y-ysm)(y - у„2 )

(8)

где а = {\l2){d2F{y0,Uai)ldy2), уо = (Кяп+Кякг)^, с определено аналогично выражению ( 5). Интегральные кривые уравнения ( 8) изображены для случая 1 на Рис. 5(6), а для случая 2 на Рис. 5(г). Одновременное существование двух

предельных циклов (случай 1) приводит к фокусировке ионных траекторий в окрестности устойчивого и дефокусировке в окрестности неустойчивого предельного цикла. В отсутствии предельных циклов (случай 2) все траектории имеют смещение одного знака. Качественная перестройка интегральных кривых динамической системы, происходящая при прохождении £/с через значение (Д, является следствием бифуркации рождения пары предельных циклов. Траектории частиц, достигающих выхода дрейф-камеры и дающих вклад в ионный ток, ограничены жирными линиями на Рис. 5(6,г).

Далее, действуя аналогично Главе 2, легко получить зависимость ионного тока /(С/со) при наличии бифуркации рождения пары предельных циклов.

и, = и,

Рис. 5 Функции смещения, траектории и ионный пик при осуществлении бифуркации рождения пары предельных циклов.

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

1

о 2

2

и5,кУ

о

В разделе 3.2. модель нелинейного дрейфа ионов, построенная в Главе 2, обобщена на случай больших значений переменной составляющей зависимости подвижности от напряженности поля, которыми характеризуются ионы легких веществ. С помощью теории возмущений получено выражение для функции смещения, в виде ряда по малому параметру {На), которое является противоположным предельным случаем по отношению к случаю малых а, рассмотренному в Главе 2. Построен ионный пик, имеющий вид, сходный с изображенным на Рис. 4(<Э).

В разделе 3.3. Главы 3 проведено сравнение выводов теории ц У спектрометрии

приращения ионной подвижности с

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

Ионограммы, получаемые в

Рис. 6 Пики бензола. 1 - расчет, 2 — эксперимент. экспериментах со

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

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

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

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

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

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

Публикации.

[1] Елистратов A.A., Шибков C.B. Анализ метода спектрометрии нелинейного дрейфа ионов для газодетекторов с плоской геометрией разделяющей камеры. Письма в ЖТФ. 2003. Т. 29. № 2.С. 88-94.

[2] Елистратов A.A., Шибков C.B. Модель метода спектрометрии нелинейного дрейфа ионов для газоанализаторов с гщлиндрической геометрией дрейф-камеры. Письма в ЖТФ. 2004, Т. 30. № 5. С. 23-29.

[3] A.A. Elistratov, S.V. Shibkov, E.N. Nikolaev Determination of the non-constant component of the ion mobility using the spectrometer of ion mobility increment. European Journal of Mass Spectrometry. Volume 12, Issue 3 (2006), p.143-151.

[4] A.A. Elistratov, S.V. Shibkov, E.N. Nikolaev Analysis of non-linear ion drift in spectrometers of ion mobility increment with cylindrical drift chamber. European Journal of Mass Spectrometry. Volume 12, Issue 3 (2006), p.153-160.

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

1 Горшков М.П. А. с. 966583 СССР // Б.И. 1982. № 38. ( М.Р. Gorshkov. Patent of USSR, # 966583 (1982).)

2 Буряков И.А., Крылов Е.В., Макась А.Л., Назаров Э.Г., Первухин В.В., Расулев У.Х. Разделение ионов по подвижности в сильных переменных электрических полях. Письма в ЖТФ. 1991. Т. 17. № 12. С. 60-65.

3 Буряков И.А., Крылов Е.В., Макась АЛ., Назаров Э.Г., Первухин В.В., Расулев У.Х. Дрейф-спектрометр для контроля следовых количеств аминов в атмосфере воздуха. ЖАХ. 1993. Т. 48. № 1.С. 156-165.

4 I.A. Buryakov, E.V. Krylov, E.G. Nazarov, U.K. Rasulev. A New Method of Separation of Multi-Atomic Ions by Mobility at Atmospheric Pressure Using a High-Frequency Amplitude-Asymmetric Strong Electric Field. Int. J. Mass Spectrom. Ion Processes. 128 (1993) 143-148.

5 И.А. Буряков, E.B. Крылов, В.П. Солдатов. Патент РФ № 1485808, МКИ G 01 N 27/62. Способ анализа микропримесей веществ в газах / Заяв. 30.03.87. 0пуб.08.02.89.

6 B.L. Carnahan, A.S. Tarassov. US Patent #5420424, 1995 (Ion Mobility Spectrometer).

7 R. Guevremont, R. W. Purves. Atmospheric Pressure Ion Focusing in a High-Field Asymmetric Waveform Ion Mobility Spectrometer. Rev. Sci. Instrum. 70 (1999) 1370-1383.

8 A. Kudryavtsev, A. Makas. Ion Focusing in a Ion Mobility Increment Spectrometer (IMIS) with Non-Uniform Electric Fields: Fundamental Considerations. International Journal of Ion Mobility Spectrometry. 4(2) (2001) 117-120.

9 G.A. Eiceman, A. Tarassov, P.A. Funk, S.E. Hughs, E.G.Nazarov, R.A. Miller. Discrimination of Combustion Fuel Sources Using Gas Chromatography- Planar Field Asymmetric Waveform Ion Mobility Spectrometry. J. Sep. Sci. 26 (2003) 585-593.

10 W. Gabryelski, F. Wu and K.L. Froese. Comparison of High-Field Asymmetric Waveform Ion Mobility Spectrometry with GC Methods in Analysis of Haloacetic Acids in Drinking Water. Anal. Chem. 75 (2003) 2478-2486.

11 Буряков И. А., Коломиец Ю.Н., Jlynny В.Б. Обнаружение паров взрывчатых веществ в воздухе с помощью спектрометра нелинейности дрейфа ионов. ЖАХ. 2001. Т. 56. № 4. С. 381-385.

12 I.A. Buryakov. Express Analysis of Explosives, Chemical Warfare Agents and Drugs with Multicapillary Column Gas Chromatography and Ion Mobility Increment Spectrometry. J.

Chromatogr. B, 800 (2004) 75-82.

13 R.W. Purves, D.A. Barnett, R. Guevremont. Separation of Protein Conformers Using Electrospray-High Field Asymmetric Waveform Ion Mobility Spectrometry-Mass Spectrometry. Int. Jour, of Mass Spectrom. 197 (2000) 163-177.

14 R. Guevremont, D.A. Barnett, R.W. Purves, J. Vandermey. Analysis of a Tryptic Digest of Pig Hemoglobin Using ESI-FAIMS-MS. Anal. Chem. 72 (2000)4577-4584.

15 I.A. Buryakov. Ion Current Amplitude and Resolution of Ion Mobility Increment Spectrometer (IMIS). International Journal of Ion Mobility Spectrometry. 4(2) (2001) 112-116.

16 E.G.Nazarov, R.A. Miller, G.A. Eiceman, E. Krylov, B. Tadjikov. Effect of the Electric Field Strength, Drift Gas Flow Rate, and Temperature on RF IMS Response. International Journal of Ion Mobility Spectrometry. 4(2) (2001) 43-46.

17 M. McCooeye, L. Ding, G.J. Gardner, C.A. Fraser, J. Lam, R. Sturgeon, Z. Mester. Separation and Quantitation of the Stereoisomers of Ephedra Alkaloids in Natural Health Products Using

Flow Injection-Electrospray Ionization-High Field Asymmetric Waveform Ion Mobility Spectrometry-Mass Spectrometry. Anal. Chem. 75 (2003) 2538-2542.

18 R. Handy, D.A. Barnett, R.W. Purves, G. Horlick, R. Guevremont. Determination of nanomolar levels of perchlorate in water by ESI-FAIMS-MS. J. Anal. At. Spectrom. 15 (2000) 907-911.

19 G.A. Eiceraan, E.G.Nazarov, R.A. Miller. A Micro-Machined Ion Mobility Spectrometer-Mass Spectrometer. International Journal of Ion Mobility Spectrometry. 3(1) (2001) 15-27.

20 R.A. Miller, G.A. Eiceman, E.G.Nazarov, A. Zapata, E. Krylov, B. Tadjikov. A Micromachined Radio Frequency Ion Mobility Spectrometer as a Gas Chromatograph Detector. International Journal of Ion Mobility Spectrometry. 4(2) (2001) 58-61.

21 I. A. Buryakov, Yu.N. Kolomiets, V.B. Louppou. Ion Non-linear Drift Spectrometer (INLDS) - a Selective Detector for High-Speed Gas Chromatography. International Journal of Ion Mobility Spectrometry. 4(1) (2001) 13-15.

22 R.W. Purves, R. Guevremont. Electrospray Ionization High-Field Asymmetric Waveform Ion Mobility Spectrometry-Mass Spectrometry. Anal. Chem. 71 (1999) 2346-2357.

23 O.A. Банных, К.Б. Поварова, В.И. Капустин. Новый подход к поверхностной ионизации и дрейф-спектроскопии органических молекул. ЖТФ. 72(12) (2002) 88-93.

24 Графский Е.Б., Клочков Д.В., Котковский Г.Е., Мочкин B.C., Першенков B.C., Тананина Е.С. Чистяков А.А. О возможности регистрации ультрамалых количеств нитросоединений методом многоступенчатой лазерной ионизации в сочетании со спектрометрией ионной подвижности. Инженерная физика. 4 (2000) 49-54.

25 Крылов Е.В. Газоразрядный ионизатор в аргоновом ионизационном детекторе. ЖТФ. 70(11) (2000) 126-132.

26 R. Guevremont, R.W. Purves, D.A. Barnett, L. Ding. Ion Trapping at Atmospheric Pressure (760 Torr) and Room Temperature with a High-Field Asymmetric Waveform Ion Mobility Spectrometer. Int. Jour, of Mass Spectrom. 193 (1999)45-56.

27 N. Krylova, E. Krylov, G.A. Eiceman, J. A. Stone. Effect of Moisture on the Field Dependence of Mobility for Gas-Phase Ions of Organophosphorus Compounds at Atmospheric Pressure with Field Asymmetric Ion Mobility Spectrometry. J. Phys. Chem. A. 107 (2003) 3648-3654.

28 Крылов Е.В. Способ уменьшения диффузионных потерь в дрейф-спектрометре. ЖТФ. 69(1) (1999) 124-127.

29 E.V. Krylov. Comparison of the Planar and Coaxial Field Asymmetrical Waveform Ion Mobility Spectrometer (FAIMS). Int. Jour, of Mass Spectrom. 225 (2003) 39-51.

30 Буряков И.А. Определение кинетических коэффициентов переноса ионов в воздухе как функций напряженности электрического поля и температуры. ЖТФ. 74(8) (2004) 15-20.

31 Буряков И. А. Решение уравнения непрерывности для ионов в газе при их движении в знакопеременном периодическом несимметричном по полярности электрическом поле. ПЖТФ. 32(2) (2006) 39-44.

32 G.E. Spangler. Relationships for Ion Dispersion in Ion Mobility Spectrometry. International Journal of Ion Mobility Spectrometry. 4(1) (2001) 71-74.

33 G.E. Spangler, R.A. Miller. Application of Mobility Theory to the Interpretation of Data Generated by Linear and RF Excited Ion Mobility Spectrometers. International Journal of Mass Spectrometry. 214 (2002) 95-104.

34 E.N. Nikolaev, A.A. Vedenov, I.A. Vedenova. The Theory of FAIMS in Coaxial Cylinders Configuration. Proceedings of the 52nd ASMS Conference on Mass Spectrometry and Allied Topics, Nashville, TN, USA, May 2004.

к исполнению 06/12/2006 Заказ №984

Исполнено 07/12/2006 Тираж: 100 экз.

Отпечатано в ООО «Стиль» г. Москва

ВВЕДЕНИЕ. ОБЗОР МЕТОДОВ РАЗДЕЛЕНИЯ ИОНОВ ПРИ АТМОСФЕРНОМ ДАВЛЕНИИ.

ГЛАВА 1. ДРЕЙФ ИОНОВ В ДРЕЙФ-КАМЕРАХ С ПЛОСКОЙ ГЕОМЕТРИЕЙ. МЕТОД ОПРЕДЕЛЕНИЯ ЗАВИСИМОСТИ ОТ НАПРЯЖЕННОСТИ ПОЛЯ НЕЛИНЕЙНОЙ СОСТАВЛЯЮЩЕЙ ПОДВИЖНОСТИ ИОНОВ.

Введение.

1.1. Модель спектрометра приращения ионной подвижности с плоской геометрией дрейф-камеры.

1.1.1. Постановка задачи. Область применимости модели.

1.1.2. Основное уравнение.

1.1.3. Основные понятия модели нелинейного дрейфа ионов.

1.1.4. Ионограмма спектрометра приращения ионной подвижности с плоской геометрией дрейф-камеры.

1.2. Определение зависимости подвижности от напряженности поля.

1.2.1. Постановка задачи.

1.2.2. Обобщенное уравнение для нелинейной составляющей подвижности.

1.2.3. Вывод интегрального уравнения для модельных и реальных временных зависимостей разделяющего напряжения.

1.2.3.1. Кусочно-постоянное разделяющее напряжение.

1.2.3.2. Трапециевидное напряжение.

1.2.3.3. Вырожденное трапециевидное напряжение.

1.2.3.4. Асимметричный меандр с экспоненциальными фронтами.

1.2.3.5. Суперпозиция двух косинусов.

1.2.4. Численное исследование.

2.1. Постановка задачи. Область применимости модели.47

2.2. Основное уравнение.51

2.3. Основные понятия теории метода спектрометрии приращения ионной подвижности. .53

2.3.1. Структура ионного шнура.53

2.3.2. Область и окно выживания.54

2.3.3. Цилиндрическое фазовое многообразие.54

2.3.4. Усреднение быстрых осцилляций.55

2.3.5. Функция смещения. Уравнение в медленном времени.56

2.3.6. Условие отбора. Предельный цикл и его устойчивость. Положение равновесия усредненной системы.56

2.3.7. Интенсивность фокусировки. Фокусировка и дефокусировка.57

2.3.8. Траектория предельного цикла в пространстве параметров динамической системы. .57 2

2.4. Модельная функция смещения.58

2.4.1. Вид модельной функции смещения.58

2.4.2. Ионограмма в режиме фокусировки.58

2.4.3. Ионограмма в режиме дефокусировки.59

2.5. Модель спектрометра для веществ с монотонной зависимостью нелинейной составляющей подвижности от напряженности поля.60

2.5.1 Общий метод получения вида пика на ионограмме.60

2.5.2. Общий вид ионограммы в режиме фокусировки.61

2.5.3. Общий вид ионограммы в режиме дефокусировки.65

2.5.4. Реальная функция смещения для цилиндрической дрейф-камеры.67

2.5.5. Ионограмма в режиме фокусировки в случае реальной функции смещения.68

2.5.6. Ионограмма в режиме дефокусировки в случае реальной функции смещения.69

Заключение.69

Рисунки к главе 2.71

ГЛАВА 3. МОДЕЛЬ СПЕКТРОМЕТРА ДЛЯ ВЕЩЕСТВ С НЕМОНОТОННОЙ ЗАВИСИМОСТЬЮ НЕЛИНЕЙНОЙ СОСТАВЛЯЮЩЕЙ ПОДВИЖНОСТИ ОТ НАПРЯЖЕННОСТИ ПОЛЯ. СЛУЧАЙ БОЛЬШИХ НЕЛИНЕЙНЫХ СОСТАВЛЯЮЩИХ ПОДВИЖНОСТИ.79

Введение.79

3.1. Модель спектрометра для веществ с немонотонной зависимостью нелинейной составляющей подвижности от напряженности поля.79

3.1.1. Введение.79

3.1.2 Реальная функция смещения в случае немонотонной зависимости подвижности от напряженности поля.80

3.1.3. Аппроксимация реальной функции смещения.82

3.1.4. Модельная немонотонная функция смещения.83

3.1.5. Траектории усредненной системы.84

3.1.6. Ионограмма в случае бифуркации рождения пары предельных циклов.86

3.2. Случай больших нелинейных составляющих подвижности.88

3.2.1. Введение.88

3.2.2. Анализ динамической системы. Функция смещения в случае больших подвижностей. .88

3.3. Сравнение с результатами экспериментальных исследований.90

3.3.1. Вид пика на ионограмме спектрометров приращения ионной подвижности с плоской геометрией дрейф-камеры.90

3.3.2. Вид пика на ионограмме и зависимость UdUs) для спектрометров приращения ионной подвижности с цилиндрической дрейф-камерой.91

Заключение.93

Рисунки к главе 3.94

ЗАКЛЮЧЕНИЕ.105

ОСНОВНЫЕ ПОЛОЖЕНИЯ, ВЫНОСИМЫЕ НА ЗАЩИТУ.107

ПРИЛОЖЕНИЕ 1.108

ПРИЛОЖЕНИЕ 2.113

ПУБЛИКАЦИИ.115

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

Введение. Обзор методов разделения ионов при атмосферном давлении.

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

Впервые метод разделения ионов в сильном переменном электрическом поле при атмосферном давлении был реализован Горшковым М.П. [1]. Экспериментальная установка представляла собой спектрометр приращения ионной подвижности с плоской геометрией дрейф-камеры. В качестве источника электронов, ионизирующих нейтральные молекулы, использовался радиоактивный источник 63Ni, а в качестве системы регистрации ионов -электрометрический усилитель. Дальнейшее развитие данный способ разделения ионов получил в работе Расулева У.Х., Бурякова И.А., Крылова Е.В., Макася А.Л., Назарова Э.Г., Первухина В.В. [2], где предложен метод регистрации ионов третичных аминов, образующихся при ионизации с помощью /^-источника, а также с помощью методов поверхностной ионизации, а также работах Бурякова И.А., Крылова Е.В., Макася A.JL, Назарова Э.Г., Первухина В.В., Расулева У.Х. [3] и [41], где были исследованы ионные процессы, сопровождающие процесс дрейфа ионов, и с помощью присоединенного масс-спектрометра был исследован состав газовой смеси на выходе дрейф-камеры спектрометра приращения ионной подвижности.

Плоская геометрия дрейф-камеры, использовавшаяся первоначально в экспериментальных установках, обладала следующим существенным недостатком -пространственно однородное электрическое поле внутри камеры не оказывало на ионы фокусирующего действия, что приводило к большим потерям ионов вследствие рекомбинации на электродах дрейф-камеры. Данная трудность была преодолена с помощью использования цилиндрической геометрии дрейф-камеры. Впервые прибор с такой конфигурацией электродов был создан Буряковым И.А., Крыловым Е.В., Солдатовым В.П. [4], а затем Карнаганом Б.Л. и Тарассовым А.С. [5]. В работе Говремонта Р., Пурвса Р.В. [42], а также независимо Кудрявцева А. и Макася A.J1. [6] было исследовано явление фокусировки ионов, которое позволяло существенно снизить потери ионов на обкладках дрейф-камеры спектрометра приращения ионной подвижности.

Исследования нового метода анализа газовой смеси при атмосферном давлении имели несколько основных направлений.

1. На основе спектрометров приращения ионной подвижности с плоской и цилиндрической геометрией дрейф-камеры в России и за рубежом был создан целый ряд портативных приборов для обнаружения следовых количеств различных веществ: для экологического мониторинга Назаровым Э.Г., Тарассовым А.С., Миллером Р.А., Айсманом Г.А., Фанком П.А., Хьюгсом С.Е. [7], а также Габриэльским В., By Ф., Фросом К. [8], для решения поисковых задач Буряковым И.А., Коломийцем Ю.Н. и Луппу В,Б. [9], [33], для экспрессного и лабораторного медицинского анализа Говремонтом Р., Пурвсом Р.В., Барнеттом Р., Вандермеем Дж. [10], [11].

2. Исследовались возможности и приборные характеристики созданных спектрометров приращения ионной подвижности: в работе Бурякова И.А. [12], в статье Назарова Э.Г., Миллера Р.А., Айсмана Г.А., Крылова Е., Таджикова Б. [13].

3. Расширялся спектр веществ, которые обнаруживались спектрометрами приращения ионной подвижности: в работах Маккой М., Динга Л., Гарднера Г.Дж., Фрейзер К.А., Лэма Дж., Стургеона Р.Е., Местера 3., Эллса Б., Говремонта Р., Пурвса Р.В., Барнетта Р. [14], [35].

4. Появлялись возможные тандемные устройства, в составе которых спектрометры приращения ионной подвижности использовались либо как предварительный фильтр или устройство концентрации пробы перед вводом в масс-спектрометры различных типов, либо как детекторы на выходе скоростной хроматографической установки или источника ионизации в виде электроспрея. На этот счет сведения содержатся в работах Ханди Р., Говремонта Р., Пурвса Р.В., Барнетта Р., Хорлика Г. [15], Назарова Э.Г., Миллера Р.А., Айсмана Г.А., Крылова Е., Таджикова Б., Запаты А. [16] [17], Бурякова И.А., Коломийца Ю.Н. и ЛуппуВ.Б. [18].

5. Исследовались новые источники ионов, применяемые для спектрометров приращения ионной подвижности, такие как электроспрей в работе Говремонта Р., Пурвса Р.В. [43], коронный разряд, поверхностная ионизация в работе Капустина В.И., Банных О.А., Поваровой К.Б. [19], лазерная ионизация в работе Графского Е.Б., Клочкова Д.В., Котковского Г.Е., Мочкина B.C., Першенкова B.C., Тананиной Е.С. Чистякова А.А. [20], высокочастотный емкостной разряд в работе Крылова Е.В. [21].

6. Рассматривались такие физические явления, как захват ионов в дрейф-камере спектрометров приращения ионной подвижности в работе Говремонта Р., Пурвса Р.В., Барнетта Р., Динга J1. [22], а также влияние влажности на данные спектрометрии в работе Айсмана Г.А., Крылова Е., Крыловой Н., Стоуна Дж.А. [23].

7. Наконец, рассматривались пути развития спектрометрии приращения ионной подвижности в работах Крылова Е.В. [24] и [25]. Рассматривались также вопросы теории спектрометрии приращения ионной подвижности в работах Бурякова И.А. [26] и [27], в работах Спенглера Г.Е. [28] и [39], а также в работе Николаева Е.Н., Веденова А.А., Веденовой И.А. [47].

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

Цель работы

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

Научная новизна

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

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

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

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

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

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

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

Апробация работы

Результаты диссертационной работы докладывались и обсуждались на семинаре Секции Прикладных проблем при Президиуме РАН (Москва, 2002 г.), на XXIX Звенигородской конференции по физике плазмы и управляемому термоядерному синтезу (Звенигород, 2002 г.), на Межведомственных конференциях по проблемам технического обеспечения безопасности (Москва, 2002 г. и 2004 г.), на Международной конференции по масс-спектрометрии (IMSC-16, Эдинбург, 2003 г.), на Втором международном семинаре-школе «Масс-спектрометрия в химической физике, биофизике и экологии» (Звенигород, 2004 г.), на Всероссийской конференции «Масс-спектрометрия и ее прикладные проблемы», проводившейся Всероссийским масс-спектрометрическим обществом (Москва, 2005 г.).

Публикации

Основные результаты работы изложены в 4 публикациях, две из которых в отечественном реферируемом журнале и две в международном реферируемом журнале. Список публикаций приведен в конце диссертации. Работа была выполнена в Институте криптографии, связи и информатики Академии ФСБ России.

Объем и структура диссертации

Диссертация состоит из введения, трех глав, заключения, двух приложений и списка цитируемой литературы. Объем диссертации составляет 119 страниц, включая 61 рисунок. Список литературы содержит 48 наименований.

Заключение.

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

Рис. 61 Алгоритм работы спектрометра приращения ионной подвижности.

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

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

Есть данные (фирма «СИБЭЛ»), что вновь создаваемые приборы будут иметь датчики физических величин, характеризующих условия дрейфа (р, ср, t). Это позволяет предполагать, что производители используют алгоритм сравнения полученной ионограммы с ионограммой, содержащейся в базе данных эталонных ионограмм. Недостатком такого подхода является приборная зависимость базы данных ионных пиков, и, что еще хуже, слишком большое количество параметров (р, (р, t, Uc, Us, Q, d), которые необходимо задать для правильного выбора пика из базы данных.

Предложенная модель позволяет закладывать в алгоритм сравнения не базу данных по пикам 7(С/с), а базу данных по зависимостям подвижности от напряженности поля к(Е). Подвижность - физическая, приборно-независимая величина, зависящая, судя по всему лишь от (р, ср, t), что резко сокращает объем базы данных и делает ее универсальной, пригодной для использования в любом приборе.

Один из используемых путей сокращения числа параметров (группа А.А. Чистякова из Московского инженерно-физического института, группа И.А. Бурякова из Конструкторско-технологического института геофизического и экологического приборостроения СОРАН, г. Новосибирск) - стабилизация свойств газа-носителя (р, (р, t) путем использования внутреннего контура чистого воздуха и газового затвора. Этот путь, к сожалению, ведет к ухудшению массо-габаритных, энергетических и ценовых показателей прибора, к необходимости регулярного технического обслуживания и пополнения расходных материалов, что существенно сужает круг потенциальных потребителей портативного оборудования, создаваемого на таком принципе.

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

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

Основные положения, выносимые на защиту.

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

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

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

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

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

1. Горшков М.П. А. с. 966583 СССР // Б.И. 1982. № 38. (М.Р. Gorshkov. Patent of USSR, # 966583 (1982).)

2. Буряков И.А, Крылов Е.В, Макась A.J1, Назаров Э.Г, Первухин В.В, Расулев У.Х. Разделение ионов по подвижности в сильных переменных электрических полях. Письма в ЖТФ. 1991. Т. 17. № 12. С. 60-65.

3. Буряков И.А, Крылов Е.В., Макась A.JI, Назаров Э.Г, Первухин В.В, Расулев УХ. Дрейф-спектрометр для контроля следовых количеств аминов в атмосфере воздуха. ЖАХ. 1993. Т. 48. № 1.С. 156-165.

4. И.А. Буряков, Е.В. Крылов, В.П. Солдатов. Патент РФ № 1485808, МКИ G 01 N 27/62. Способ анализа микропримесей веществ в газах / Заяв. 30.03.87. 0пуб.08.02.89.

5. B.L. Carnahan, A.S. Tarassov. US Patent #5420424, 1995 (Ion Mobility Spectrometer).

6. A. Kudryavtsev, A. Makas. Ion Focusing in a Ion Mobility Increment Spectrometer (IMIS) with Non-Uniform Electric Fields: Fundamental Considerations. International Journal of Ion Mobility Spectrometry. 4(2) (2001) 117-120.

7. G.A. Eiceman, A. Tarassov, P.A. Funk, S.E. Hughs, E.G.Nazarov, R.A. Miller. Discrimination of Combustion Fuel Sources Using Gas Chromatography- Planar Field Asymmetric Waveform Ion Mobility Spectrometry. J. Sep. Sci. 26 (2003) 585-593.

8. W. Gabryelski, F. Wu and K.L. Froese. Comparison of High-Field Asymmetric Waveform Ion Mobility Spectrometry with GC Methods in Analysis of Haloacetic Acids in Drinking Water. Anal. Chem. 75 (2003) 2478-2486.

9. Буряков И.А, Коломиец Ю.Н, Jlynny В.Б. Обнаружение паров взрывчатых веществ в воздухе с помощью спектрометра нелинейности дрейфа ионов. ЖАХ. 2001. Т. 56. № 4. С. 381-385.

10. R.W. Purves, D.A. Barnett, R. Guevremont. Separation of Protein Conformers Using Electrospray-High Field Asymmetric Waveform Ion Mobility Spectrometry-Mass Spectrometry. Int. Jour, of Mass Spectrom. 197 (2000) 163-177.

11. R. Guevremont, D.A. Barnett, R.W. Purves, J. Vandermey. Analysis of a Tryptic Digest of Pig Hemoglobin Using ES1-FAIMS-MS. Anal. Chem. 72 (2000) 4577-4584.

12. I. A. Buryakov. Ion Current Amplitude and Resolution of Ion Mobility Increment Spectrometer (IMIS). International Journal of Ion Mobility Spectrometry. 4(2) (2001) 112-116.

13. E.G.Nazarov, R.A. Miller, G.A. Eiceman, E. Krylov, B. Tadjikov. Effect of the Electric Field Strength, Drift Gas Flow Rate, and Temperature on RF IMS Response. International Journal of Ion Mobility Spectrometry. 4(2) (2001) 43-46.

14. R. Handy, D.A. Barnett, R. W. Purves, G. Horlick, R. Guevremont. Determination of nanomolar levels of perchlorate in water by ESI-FAIMS-MS. J. Anal. At. Spectrom. 15 (2000) 907-911.

15. G.A. Eiceman, E.G.Nazarov, R.A. Miller. A Micro-Machined Ion Mobility Spectrometer-Mass Spectrometer. International Journal of Ion Mobility Spectrometry. 3(1) (2001) 15-27.

16. R.A. Miller, G.A. Eiceman, E.G.Nazarov, A. Zapata, E. Krylov, B. Tadjikov. A Micromachined Radio Frequency Ion Mobility Spectrometer as a Gas Chromatograph Detector. International Journal of Ion Mobility Spectrometry. 4(2) (2001) 58-61.

17. I.A. Buryakov, Yu.N. Kolomiets, V.B. Louppou. Ion Non-linear Drift Spectrometer (INLDS) a Selective Detector for High-Speed Gas Chromatography. International Journal of Ion Mobility Spectrometry. 4(1) (2001) 13-15.

18. O.A. Банных, К.Б. Поварова, В.И. Капустин. Новый подход к поверхностной ионизации и дрейф-спектроскопии органических молекул. ЖТФ. 72(12) (2002) 88-93.

19. Крылов Е.В. Газоразрядный ионизатор в аргоновом ионизационном детекторе. ЖТФ. 70(11) (2000) 126-132.

20. R. Guevremont, R.W. Purves, D.A. Barnett, L. Ding. Ion Trapping at Atmospheric Pressure (760 Torr) and Room Temperature with a High-Field Asymmetric Waveform Ion Mobility Spectrometer. Int. Jour, of Mass Spectrom. 193 (1999) 45-56.

21. Крылов Е.В. Способ уменьшения диффузионных потерь в дрейф-спектрометре. ЖТФ. 69(1) (1999) 124-127.

22. E.V. Krylov. Comparison of the Planar and Coaxial Field Asymmetrical Waveform Ion Mobility Spectrometer (FAIMS). Int. Jour, of Mass Spectrom. 225 (2003) 39-51.

23. Буряков И.А. Определение кинетических коэффициентов переноса ионов в воздухе как функций напряженности электрического поля и температуры. ЖТФ. 74(8) (2004) 15-20.

24. Буряков И. А. Решение уравнения непрерывности для ионов в газе при их движении в знакопеременном периодическом несимметричном по полярности электрическом поле. ПЖТФ. 32(2) (2006) 39-44.

25. G.E. Spangler. Relationships for Ion Dispersion in Ion Mobility Spectrometry. International Journal of Ion Mobility Spectrometry. 4(1) (2001) 71-74.

26. И. Мак-Даниель, Э. Мэзон. Подвижность и диффузия ионов в газах. М.: Мир, 1976 г., 424 с. (Mason Е.А., McDaniel E.W. Transport Properties of Ions in Gases. John Willey & Sons: New York, 1988. 560 p.)

27. Ландау Л. Д., Лифшиц Е. М. Гидродинамика. Теоретическая физика, т. VI. М.: Наука, 1986.

28. Андронов А.А., Леонтович Е.А., Гордон И.И. Майер А.Г. Теория бифуркаций динамических систем на плоскости. М.: Наука. 1967. 488 с.

29. Федоренко Р.П. Введение в вычислительную физику. М.: МФТИ, 1994.

30. I.A. Buryakov. Express Analysis of Explosives, Chemical Warfare Agents and Drugs with Multicapillary Column Gas Chromatography and Ion Mobility Increment Spectrometry. J. Chromatogr. B, 800 (2004) 75-82.

31. D.A. Barnett, B. Ells, R. Guevremont, R.W. Purves, L.A. Viehland. Evaluation of Carrier Gases for Use in High-Field Asymmetric Waveform Ion Mobility Spectrometry. J. Am. Soc. Mass Spectrom. 11 (2000) 1125-1133.

32. B. Ells, D.A. Barnett, R.W. Purves, R. Guevremont. Detection of Nine Chlorinated and Brominated Haloacetic Acids at Part-per-Trillion Levels Using ESI-FAIMS-MS. Anal. Chem. 72 (2000)4555-4559.

33. Буряков И.А. Экспериментальное определение зависимости коэффициентов подвижности ионов в газе от напряженности электрического поля. ЖТФ. 2002. Т. 72. № 11. С. 109-113.

34. L.A. Viehland, R. Guevremont, R.W. Purves, D.A. Barnett. Comparison of High-Field Ion Mobility Obtained from Drift Tubes and a FAIMS Apparatus. International Journal of Mass Spectrometry. 197 (2000) 123-130.

35. G.E. Spangler, R.A. Miller. Application of Mobility Theory to the Interpretation of Data Generated by Linear and RF Excited Ion Mobility Spectrometers. International Journal of Mass Spectrometry. 214 (2002) 95-104.

36. R. Guevremont, R.W. Purves. Atmospheric Pressure Ion Focusing in a High-Field Asymmetric Waveform Ion Mobility Spectrometer. Rev. Sci. Instrum. 70 (1999) 1370-1383.

37. R.W. Purves, R. Guevremont. Electrospray Ionization High-Field Asymmetric Waveform Ion Mobility Spectrometry-Mass Spectrometry. Anal. Chem. 71 (1999) 2346-2357.

38. Колмогоров А.Н, Фомин С.В. Элементы теории функций и функционального анализа. М.:Наука, 1989, 624с.

39. Андронов А.А, Витт А.А, Хайкин С.Э. Теория колебаний. М.: ГИФМЛ. 1959. 916 с.

40. Андронов А.А, Леонтович Е.А, Гордон И.И, Майер А.Г. Качественная теория динамических систем второго порядка. М.: Наука. 1966. 568 с.

41. E.N. Nikolaev, А.А. Vedenov, I.A. Vedenova. The Theory of FAIMS in Coaxial Cylinders Configuration. Proceedings of the 52nd ASMS Conference on Mass Spectrometry and Allied Topics, Nashville, TN, USA, May 2004.

42. В.И. Арнольд. Геометрические методы в теории обыкновенных дифференциальных уравнений. М.: Наука, 1999. 400 с.