CC BY
10ISSN 0868-5886 НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2013, том 23, № 3, c. 69-75
= ТЕОРЕТИЧЕСКИЕ ИССЛЕДОВАНИЯ =
УДК 621.384.8
© Н. Д. Семкин, Д. В. Родин, И. В. Пияков
МЕТОД РАСЧЕТА ПОТЕНЦИАЛОВ ПОЛЕЗАДАЮЩИХ ЭЛЕМЕНТОВ НЕЛИНЕЙНОГО ИОННОГО ЗЕРКАЛА
Описан метод расчета неоднородного электростатического поля отражателя масс-рефлектрона. Метод позволяет рассчитать распределение осевого потенциала в отражателе, теоретически исключающее дисперсию суммарного времени пролета для осевых ахроматичных ионных пакетов в заданном диапазоне начальной энергии ионов. Методом численного моделирования проведен анализ влияния погрешностей физического синтеза поля в близком к реальному отражателе на точность фокусирования. Предложена методика расчета отражателя с неоднородным полем. Приведено описание алгоритма оптимизации потенциалов полезадаю-щих элементов. Даны рекомендации по уменьшению времени выполнения алгоритма. Описана конструкция, используемая для проверки алгоритма оптимизации, даны результаты моделирования процессов фокусировки для ионов различных масс и различных начальных энергетических разбросов.
Кл. сл.: времяпролетный масс-спектрометр, масс-анализатор, масс-рефлектрон, нелинейное ионное зеркало, генетический алгоритм оптимизации
ВВЕДЕНИЕ
Использование времяпролетных масс-спектрометров позволяет минимизировать массогабарит-ные и энергетические параметры аппаратуры при сохранении ее высоких функциональных характеристик, что особенно важно при проведении экспериментов по исследованию космического пространства, изучению химического состава газопылевой компоненты собственной внешней атмосферы космических аппаратов.
Применение ионных зеркал в конструкции прибора позволяет улучшить характеристики масс-спектрометра за счет увеличения пути ионов на бесполевых участках, а также реализации принципа пространственно-временной фокусировки ионных пакетов в плоскости детектора. Использование этих принципов обеспечивает разрешение R ~ 500 при начальном энергетическом разбросе ионов до 15 % от ускоряющего напряжения [1-3]. Однако с ростом энергетического разброса увеличивается также ошибка компенсации времени пролета, что ограничивает возможности применения данного типа анализаторов, поскольку исследование проб с высоким начальным энергетическим разбросом (например, полученных в результате высокоскоростного взаимодействия микрометеорита и мишени) требует использования высоких ускоряющих напряжений и быстродействующей аппаратуры обработки ионных спектров.
Метод расчета нелинейного ионного зеркала, предложенный в [4], позволяет обеспечить незави-
симость суммарного времени пролета ионов от их начальной энергии. Фокусировка осуществляется за счет того, что ионам с большей энергией соответствуют более длинные траектории внутри зеркала, таким образом компенсируется меньшее время пролета бесполевого участка.
Авторами был разработан метод аналитического расчета нелинейного электростатического зеркала, лишенный вышеприведенных недостатков [5]. Данный метод позволяет получить форму потенциала (1), осуществляющую точную компенсацию времени пролета для ионов с различной начальной энергией:
= - \to^4p> - Lo
П I
arcsin,
Р
1 + р
- Li
^/р + (1 + р) arcsin
Р
п
1 + р 2
Р
(1)
где — время пролета частицы с нулевой начальной энергией; а — скорость частицы после ускоряющего промежутка; ф — нормированный потенциал; L0 — суммарная длина бесполевых участков; L1 — суммарная длина участков ускорения и торможения.
Однако физический синтез поля с заданным распределением осевого потенциала путем установки потенциалов полезадающих элементов по приведенному выше выражению невозможен. Это вызвано тем, что из-за просадок поля, вызванных
z
различными факторами, осевое распределение отличается от расчетного.
АЛГОРИТМ ОПТИМИЗАЦИИ ПОТЕНЦИАЛОВ ПОЛЕЗАДАЮЩИХ ЭЛЕМЕНТОВ
Для решения задачи восстановления расчетного распределения осевого потенциала авторами была разработана программа расчета потенциалов поле-задающих элементов с использованием генетического алгоритма оптимизации: на каждой итерации работы алгоритма генерируется набор потенциалов для полезадающих элементов. Для каждого члена набора вычисляется поле на расчетной сетке, и осевое распределение потенциала сравнивается с целевым. Четыре члена набора с минимальным отклонением осевого распределения от расчетного участвуют в генерации набора для следующей итерации: их потенциалы "скрещиваются" и вносятся псевдослучайные "мутации" (поэтому алгоритм и получил название "генетический"). После нахождения минимума отклонения осевого распределения алгоритм переходит в режим работы с целевой функцией максимизации разрешения. Для этого в программе был реализован модуль трассировки ионов и расчета разрешения. Окончанием работы алгоритма служит установившийся режим по максимизации разрешения, т. е. от итерации к итерации увеличения разрешения не происходит.
Электростатическое поле внутри отражателя было рассчитано экстраполяционным методом Либмана, который был предложен Франкелом [6] и Янгом [7]. Данный метод является наиболее универсальным, поскольку обеспечивает достаточно быструю сходимость при сохранении устойчивости решения:
An+1 _ An I
An,k - An,k +
+ 4 ( aCu + <k+i + ACU + A"nk-1 + hW - Ak)
(2)
или сходимости, определяющим скорость сходимости. Значение должно находиться в пределах 1^2. При а = 1 уравнение (2) превращается в уравнение для обычного метода Либмана, а при а > 2 процесс решения становится неустойчивым.
Значение а можно выбрать способом, предложенным Янгом [8]:
f
>2
1 - ^ + Л
\p q
\
(3)
где аь — оптимальное значение коэффициента сходимости; р — количество узлов сетки по вертикали; q — количество узлов сетки по горизонтали.
После расчета матрицы потенциалов требуется перейти к напряженности электрического поля. В соответствии с уравнением Максвелла, можно перейти от частных производных к приращениям (при условии их достаточной малости)
E «-
AV
Ах '
E»-*V-,
у Ау
(4)
Здесь а является коэффициентом релаксации,
где V — потенциал электрического поля.
Поскольку задача вычисления поля для сеток с большим количеством элементов является требовательной к вычислительным ресурсам, а метод оптимизации подразумевает большое количество итераций таких вычислений, код метода вычисления поля был портирован на язык CUDA, что позволило существенно ускорить решение задачи. Для большего ускорения возможно использование промежуточных расчетов на более грубых сетках с последующей интерполяцией результатов. Сравнительное время выполнения одной итерации вычисления для кода CPU, CUDA и CUDA c промежуточной интерполяцией приведено в таблице. Рост ускорения с ростом размерности сетки при переходе на CUDA объясняется уменьшением влияния накладных расходов при пересылке блоков данных между host-памятью и device-памятью.
а
b
Сравнительное время выполнения одной итерации вычисления (с)
Способ вычислений Размерность сетки
64X64 128X128 192X192 256X256
CPUX4, 3.4 ГГц 3.56 39.11 3-60 + 1.51 10-60 + 5.51
CUDA GF560 1.18 3.66 11.61 26.58
Interpolate 0.99 2.9 7.5 16.1
МЕТОД РАСЧЕТА ПОТЕНЦИАЛОВ...
71
В
MCP
Ш
L L
I I
1о1
125
¿12
-I & а-
Ф1 Л Ф1
L Li i I 11 11 iL
i
i ■ п i
■н
l23
ф1
-0 11
I
34
I I I I I I I I
фп
+0
п
- ■> Z
Рис. 1. Схема исследуемого масс-спектрометра
Расчет траекторий производился с использованием известных уравнений движения заряженных частиц в электрических и магнитных полях. Численный расчет траектории производился с шагом по времени, обеспечивающим, по крайней мере, 100 шагов траектории на каждую ячейку расчетной сетки. Поле в точках местонахождения ионов рассчитывалось методом бикубической интерполяции.
Воспользовавшись уравнениями (4) и переходя от дифференциала к приращению времени, получим следующие уравнения, на основе которых будет составляться алгоритм:
x(t + At) = x(t) +1 vx (t) + At I At.
y(t + At) = y (t) + Vv (t) + —^ At
m
qEy
(5)
m
At,
1 N
T . = — У T ,
ср< \T /—i J '
Nj-1 J
(6)
где Т — время прихода у-го иона.
Среднеквадратическое отклонение определялось согласно формуле
GT =
1 N _
N у (т - т„ )2
j=1
(7)
Разрешение по сигме:
T - T
г\ _ ср.+1 срг
Rg' = 2G
где GT — среднеквадратическое отклонение.
(8)
где ух, уу — соответственно горизонтальная и вертикальная составляющие скорости частицы в текущий момент времени.
Величина пакетов выбиралась, исходя из условий гладкости полученного максвелловского распределения по энергии, и составила 100 000 шт. на каждую массу. Далее для Мi массы вычислялось среднее время пролета
РЕЗУЛЬТАТ РАБОТЫ АЛГОРИТМА
Рассмотрим работу программы на примере конструкции масс-спектрометра, приведенной на рис. 1. Здесь р — потенциал выталкивающего поля; рп — максимальный эквивалентный потенциал, соответствующий начальному энергетическому разбросу ионов; l12 — участок бесполевого пространства; l0i — участок ускорения ионов; l23 — линейный участок электронно-оптической фокусирующей системы; l34 — нелинейный участок фокусирующей системы; l25 — участок бесполевого пространства при движении ионов от отражателя к приемнику (MCP).
Для расчета были выбраны следующие геометрические размеры конструкции: суммарная длина бесполевого участка L = l12 + l25, равная 12.5 см, ширина тормозящего промежутка l23 = 2.5 см, ширина ускоряющего промежутка l01 условно равна нулю.
АЕ/Е
Rma- 1500 7
Рис. 2. Промежуточные результаты процесса оптимизации.
По левой оси ординат — средняя относительная ошибка восстановления расчетного потенциала; по правой оси ординат — разрешение прибора на 1 а.е.м.; по оси абсцисс — номер итерации; кривая 1 — зависим ость разрешения от номера итерации; кривая 2 — зависимость средней относительной ошибки от номера итерации
4600 Е, В/м
4500 4400
X, м
Рис. 3. Зависимость относительной ошибки от координаты внутри зеркала и целевое распределение поля внутри зеркала.
Кривая 1 — целевая зависимость поля от координаты; кривая 2 — синтезированное осевое распределение внутри отражателя; кривая 3 — зависимость относительной ошибки восстановления от координаты при оптимизации по целевой функции минимизации ошибки восстановления; кривая 4 — зависимость относительной ошибки восстановления от координаты при оптимизации по целевой функции максимизации разрешения
Р ис. 4. Распределение поля внутри нелинейного и онного зеркала в зависимости от радиальной и продольной координат
Рис. 5. Распределение относительной ошибки внутри нелинейного ионного зеркала в зависимости от радиальной и продольной координат
Результатом работы данной программы является набор потенциалов для полезадающих элементов, при использовании которых отклонение распределения осевого потенциала от расчетного является минимальным. Промежуточные результаты процесса оптимизации приведены на рис. 2.
На рис. 2 хорошо виден момент перехода программы в режим работы с целевой функцией максимизации разрешения: он приходится примерно на 190-ю итерацию. Небольшое ухудшение средней относительной ошибки восстановления рас-
четного распределения потенциала при одновременном улучшении разрешения может быть объяснено тем, что при работе программы по целевой функции уменьшения средней относительной ошибки не учитывается распределение этой ошибки по координате внутри зеркала. При переходе программы в режим работы по максимизации разрешения ошибка восстановления может быть видоизменена к распределению, меньше влияющему на время пролета иона в зеркале, при этом суммарная ошибка может возрасти, но разрешение повыситься.
МЕТОД РАСЧЕТА ПОТЕНЦИАЛОВ.
73
На рис. 3 приведено целевое и полученное в результате работы алгоритма осевые распределения электростатического поля, а также ошибки восстановления поля для конечных моментов работы алгоритма по целевой функции минимизации ошибки восстановления осевого распределения потенциала и по максимизации разрешения прибора.
Можно видеть, что наше предположение о распределении ошибки восстановления формы потенциала подтвердилось, средняя величина ошибки возросла, поскольку после оптимизации по целевой функции распределение ошибки восстановления было знакопеременным, а после оптимизации по разрешению ошибка имеет постоянный знак, при этом ее средняя величина увеличивается незначительно. Как будет показано далее, такое
изменение распределения ошибки приводит к возрастанию разрешающей способности.
Используя полученные значения потенциалов полезадающих колец, было рассчитано распределение поля в зеркале и рассчитана ошибка восстановления, приведенные на рис. 4 и 5.
Из рис. 4 хорошо видно, что распределение осевой составляющей поля в объеме зеркала являются достаточно гладким и для околоосевых траекторий ионов влияние неточностей восстановления сведено к минимуму. Ошибка восстановления поля на оси прибора не превышает 3 % вблизи сеток, ограничивающих зеркало, и 0.5 % в остальной части отражателя. На периферии рефлектора относительная ошибка восстановления не превышает 10 %.
2000 -| 1800 -1600 -1400 -1200 -1000 -800 -600 -400 -200 -0 -
R
Рис. 6. Разрешение в зависимости от энергии ионов. 1 — нелинейное зеркало с оптимизацией по разрешению; 2 — нелинейное зеркало с оптимизацией по минимальной ошибке; 3 — линейное зеркало
80 100 ДЕ, эВ
10000
1000
100
10
0.1
1
2
Рис. 7. Разрешение в зависи-
3 мости от массы.
1 — нелинейное зеркало, 100 % разброс; 2 — нелиней-
4 ное зеркало, 50 % разброс; 3 — линейное зеркало, 50 % разброс; 4 — линейное зерка-
20
40
60
М,
100 а.е.м.
Зависимости разрешения от энергетического разброса для линейного зеркала и нелинейного зеркала с оптимизацией по целевой функции минимизации ошибки и по разрешению приведены на рис. 6. Как видно из приведенных графиков, для результирующего распределения с минимальной ошибкой восстановления максимальное разрешение достигается при энергетическом разбросе ионов, не превышающем 10 % от ускоряющего напряжения, затем разрешение резко снижается. Для распределения, оптимизированного по разрешению, фокусировка ионов по энергиям превосходит результаты первого способа оптимизации более чем в 2 раза на всем диапазоне энергий, за исключением диапазона малых разбросов, где результаты первого способа превосходят результаты второго.
Зависимости разрешения от массы для 50 и 100 % начального энергетического разброса для традиционного рефлектрона и масс-спектрометра с нелинейным зеркалом приведены на рис. 7. Можно видеть, что для нелинейного зеркала зависимость разрешения от энергетического разброса минимальна и превосходит 10 единиц даже для масс порядка 100 а.е.м.
ЗАКЛЮЧЕНИЕ
СПИСОК ЛИТЕРАТУРЫ
1. Каратаев В.И., Мамырин Б.А., Шмикк Д.В. Новый принцип фокусировки ионных пакетов во времяпролeтных масс-спектрометрах // ЖТФ. 1971. Т. 41, вып. 7. С. 1498-1501.
2. Каратаев В.И., Мамырин Б.А., Шмикк Д.В., Загу-лин В.А. Новый безмагнитный времяпролeтный масс-спектрометр с высокой разрешающей способностью // ЖТФ. 1973. Т. 64, вып. 1. С. 82-89.
3. Шмикк Д.В., Дубенский Б.Н. Отражатель масс-рефлектрона // ЖТФ. 1984. Т. 58, вып. 1. С. 912916.
4. Глащенко В.П., Сёмкин Н.Д. Расширение энергетического диапазона фокусируемых ионов во времяпролeтном масс-спектрометре // ЖТФ. 1985. Т. 55, вып. 5. С. 904-907.
5. Сёмкин Н.Д., Пияков И.В., Родин Д.В., Помельни-ков Р.А. Аналитический метод расчета распределения электростатического поля отражателя время-пролетного масс-спектрометра // ЖТФ. 2012. Т. 82, вып.10. С. 79-84.
6. Frankel S.P. Convergence rates of iterative treatments of partial differential equations // Math. Tables and Other Aids to Computation. 1950. N 4. Р. 65-74.
7. Southwell R. V. Relaxation methods in theoretical physics. University Press, Oxford, 1946. 260 р.
8. Шмикк Д.В., Дубенский Б.Н. Отражатель масс-рефлектрона // ЖТФ. 1984. Т. 54, вып. 5. С. 913916.
Как видно из приведенных графиков, предложенный метод превосходит классический линей- Самарский государственный аэрокосмический уни-ный рефлектрон в 10-50 раз в зависимости от на- верситет им. акад. С.П. Королева (Национальный чального энергетического разброса ионов. Данный исследовательский университет) метод аналитического расчета нелинейного распределения потенциала электростатического ионного зеркала в совокупности с применением алгоритма генетической оптимизации потенциалов
полезадающих элементов позволяет существенно Контакты: Родин Дмитрий Владимирович, расширить диапазон разрешаемых масс, а также rodin.SSaU@gmail.COm расширить диапазон допустимых начальных энергий ионов вплоть до 100 %. Материал поступил в редакцию 5.04.2013
THE METHOD OF CALCULATION OF POTENTIALS FOR THE FIELD DETERMINATION ELEMENTS OF NON-LINEAR ION MIRROR
N. D. Semkin, D. V. Rodin, I. V. Piyakov
Samara State Aerospace University named after academician S.P. Korolyov (National Research University)
A method for calculating of the non-uniform electrostatic field of the mass reflection reflector is described. The method allows to calculate the axial potential distribution in the reflector theoretically eliminating dispersion of the total time of flight for axial achromatic ion packets in a given range of the initial energy of the ions.
METQA PACTETA nQTEHUHAHQB.
75
By the method of the numerical simulation the analysis of the errors influence for the field physical synthesis close to the real reflector on the focusing accuracy is carried out. The calculation procedure of the reflector with an inhomogeneous field is proposed. The description of the potentials optimization algorithm for field-determination elements is given. The recommendations to reduce execution time of the algorithm are given. The design used for verification of the optimization algorithm is described; the results of simulation of the focusing procedures for ions of different masses and different initial energy spreads are given.
Keywords: time of flight mass spectrometer, mass analyzer, mass-reflectron, non-linear ion mirror, genetic optimization algorithm
HAYHHQE nPHEQPQCTPQEHHE, 2013, tom 23, № 3
