Труды МАИ. Выпуск № 98
http://trudymai.ru/
УДК 533.9, 539.9, 539.17-539-19
Разработка кинетических моделей движущейся плазмы. Коэффициенты Эйнштейна для иона ксенона.
Кули-заде М.Е.*, Решетникова О.Ф.**, Скороход Е.П.***
Московский авиационный институт (национальный исследовательский университет), МАИ, Волоколамское шоссе, 4, Москва, A-80, ГСП-3, 125993, Россия
*e-mail: kuli_marina@mail. ru **e-mail: resholga@yandex.ru ** *e-mail: e.p.skorohod@mail. ru
Аннотация
Изучение потоков движущейся плазмы, как экспериментального, так и теоретического, а также фундаментальные вопросы плазмомеханики в настоящий момент относятся к быстро развивающимся областям науки, находящим применение при решении задач двигателестроения, ракетно-космической отрасли. В частности моделирование физико-химических процессов необходимо для задач многоуровневой кинетики, спектроскопии низкотемпературной плазмы, оптической диагностики плазменных потоков. Коэффициенты Эйнштейна спонтанных радиационных переходов являются основными слагаемыми многочисленных кинетических уравнений, используются для построения векторов состояний исследуемой плазмы при обработке спектров.
Приведены значения сил осцилляторов и вероятностей ~1000 радиационных переходов иона ксенона. Даны рекомендации по использованию приведённых
табличных данных Хе11 в соответствие со схемой учитываемых уровней и конфигураций.
Ключевые слова: плазма ксенона, моделирование радиационных процессов в кинетике, силы осцилляторов иона Хе11, коэффициенты Эйнштейна спонтанных радиационных переходов иона Хе11, вероятности фото-переходов иона Хе11, оптическая диагностика.
Введение
Настоящая работа является продолжением рассмотрения вопросов плазмомеханики [1- 4]. Разработка кинетических моделей [5-7, 22-25] по-прежнему актуальна, поскольку способствует развитию современной спектроскопии плазмы, становлению новых экспериментальных методов пассивной и активной оптической диагностики (например, метод лазерно-индуцированной флуоресценции) и др.
Изучение свойств низкотемпературной плазмы, в частности, вектора состояния [1- 4] движущейся ксеноновой плазмы, как экспериментального, так и теоретического, невозможно без знания вероятностей фото-переходов Лш
(коэффициентов Эйнштейна спонтанных радиационных переходов Лш), когда
вероятности спонтанного излучения (испускания), поглощения и вынужденного излучения характеризуются соответственно коэффициентами Эйнштейна Лш, В к1 ,
С3
В1к как glBlk = giBы = gJ Аы—— .
8жпу
Для расчёта коэффициентов Эйнштейна Лш необходимо определиться со схемой уровней, на базе которой формируются система кинетических уравнений и
релаксационная матрица [8].
Как отмечено в [9]: «Точное предсказание параметров спектральных атомных линий является одной из актуальных проблем современной атомной физики. В радиационной газовой динамике принципиально важно предсказание параметров как можно большего числа линий, в том числе достаточно слабых» (имеются в виду данные [10]).
В работах из списка, цитируемого в [1, 2], по изучению ксеноновой плазмы Холловского двигателя мало внимания уделяется вопросу: какие уровни и конфигурации следует учитывать в расчётах. Обычно в зарубежных публикациях рассматриваются усреднённые блоки, "полки", с использованием экспериментально полученных констант, отдельно пишется кинетика для атома, отдельно для иона (совместно не рассматриваются).
Следует заметить, что усреднение радиационных констант связано с определением < - фактора матричного элемента [11] (сложением моментов дипольных переходов). Предлагается более детально осуществить составление схемы уровней и выбор оптических констант для составления кинетических уравнений.
1. Схемы уровней ксенона
Общая схема уровней атома Хе1 и ионов ксенона Хе11 - Хе1У, принадлежащих основным электронным конфигурациям приведена в [3]. Атом
2 5 2
ксенона, относящийся к восьмой группе, имеет два предела ионизации: 5^ 5р Р3/2 , рис.1 и для смещенной (штрихованной) схемы уровней 5s2 5р5 2Р]/2, - рис.2.
я I
1(?смл
Пр /-
пс1
Л Г
Л А"
эВ
к ъ ь ь ъ ъ ъ ь ь ь ъ
'2 /2 12
J 1 2 0 1 1 2 2 3 0 1 1 2 2 3 3 4 1 2 2 3 3 4 4 5
&
О Т-
"- 12
4/
п=3
бр
5(1
п=2
бз
¡¿2$2рб3823р63с(04*4рб4Л05*5рб1*0
■П=4
I- и
- 10
- 9
- 8
Е, + и)4смл
10
Рис.1. Схема уровней атома XeI.
ил пр пс1 тр
/-^-л /-V Е>
К ъ 42 Ъ Ч2 Ъ т эВ
] 0 1 0 112 12 2 3
7в'
п=3 I
13 I*- п=4
бр'
5(1
- п=2
12
11
10
Рис.2. Смещённая схема уровней атома ксенона XeI.
Рис. 3. Схема уровней иона Хе11.
1<Г
см1
19 18 17 16 15 14 13 12 11
п (Л1
\
П 5
-\
2р 2В 2р ^ 2р 2В 2р
ъъ ъъъъъъ ъьъъъъъ
<&
бр'
5 а'
5s25p52P3
П 5
П ¿1
" 2И
>р" п
и Е
6р'
68
\
3/2
эВ
25 24 23 22 21 20 19 18 17 16 15 14 13
Рис. 4. Смещённые схемы уровней иона Хе11.
Схема уровней иона ксенона приведена на рис.3 и смещённые схемына
рис.4.
В рамках, так называемой, излучательно-столкновительной модели [4] рассматриваются: реакции спонтанных радиационных переходов, возбуждения электронным ударом, процессы ионизации и рекомбинации. Все бинарные взаимодействия в виде констант скоростей в слагаемых кинетических уравнений записываются согласно выбранной схеме уровней.
Рассмотрим пример составления схемы многоуровневой кинетической модели, которая включает уровни и (что принципиально) их группы, как для атома Хе1, так и для иона Хе11.
Схема атома Хе1 составляется из: 4 уровней (6, 7я конфигураций); 12 уровней для 6, 7 р конфигураций; 16 уровней для 5, 6 й конфигураций; 4 уровня для 6, 7 я' конфигураций; 8 уровней для 6, 7 р конфигураций; 8 уровней для 5, 6й' конфигураций, а также для конфигураций 8я, 8р, 7й, 8я', 8р', 7й' и водородоподобных уровней 5, 6 п (всего 60). Смещенные уровни можно учитывать выборочно, что уменьшит общее число до 41.
Для иона Хе11 включаются: - 10 уровней для 6, 7 я конфигураций; 26 уровней для 6, 7 р конфигураций; 34 уровня для 5, 6 й конфигураций; 4 уровня для 6, 7 я' конфигураций; 12 уровней для 6, 7 р' конфигураций; 14 уровней для 5, 6 й' конфигураций; 2 уровня для 6, 7 я'' конфигураций; 2 уровня для 6, 7 р'' конфигураций; 4 уровня для 5, 6 й'' конфигураций, а также для конфигураций 8 я, 8 р
, 7 ё, 8 s', 8 р', 7 ё' и водородоподобных уровней 6, 7 и (всего 108). Смещенные уровни с двумя штрихами можно не учитывать, что уменьшит общее число до 100.
Таким образом, размерность релаксационной матрицы - п = 168 х168 в первом случае и п = 141x141 - во втором.
Обратимся к экспериментальным исследованиям [1,2] плазмы Холловского двигателя. Регистрируемое число линий ксенона в спектре для атома составляет 95, а для первого иона - 149, т.е. вектор состояния, представлен 244 значениями. Для этих целей число кинетических уравнений составит 141 (второй вариант) для атома и первого иона и 12 уравнений для второго иона.
Чтобы решать теоретически задачу [1,2] необходимо определиться с многочисленными константами и решать задачу многоуровневой кинетики совместно с функцией распределения электронов по энергиям.
Сокращение числа уровней (рис.1.-рис.4.) с целью упрощения вычислительных алгоритмов уменьшает вкладываемую энергию в плазму, меняет фактор влияния согласно [8], которыми являются рекомбинационные процессы.
Одной из характеристик возбуждённых состояний атома являются вероятности радиационных переходов или коэффициенты Эйнштейна для спонтанного излучения.
где АЕ=Иа> разность энергий верхнего и нижнего уровней, Яу = 13.6 эВ, gl = 2J0 +1 -статистический вес нижнего уровня, gk = 2^ +1 - статистический вес верхнего
2. Силы осцилляторов и вероятности фото-переходов иона ХеП
(1)
уровня. Сила осциллятора в поглощении /1к перехода 1 - к в кулоновском приближении [11] равна
Л = в, £ к) . (2)
Здесь £, £ к - орбитальные моменты оптического электрона в начальном и конечном
состояниях, величина р(£ ^£к) указана в [11]. р(£ г,£к) находится как функция
эффективных главных квантовых чисел уровней п * и пк , которые указаны в табл.1-табл.4 электронного приложения, которое можно посмотреть на сайте:
http://tahir.sytes.net/!books/articles/Приложение к статье.pdf [26] (аналогичную
*
процедуру можно повторить и для эффективного главного числа п всей конфигурации).
В случае иона Хе11 угловой фактор О характеризуется ЬЯ - связью. Обозначим начальное состояние индексом «0», а конечное «1». Для перехода [ЯрЦр ] 1о Я 1о-10 ^ \5рЬр ] 11 Я Ц угловой фактор О равен
Оо1 = (2Ц + Ш +1){ЦЦ Я}2 & +1)(2^1 +1){ЦЦ Ц}2 (3)
где 10,11 - орбитальные моменты оптического электрона в начальном и конечном состояниях, Я - полный спин (рассматриваются переходы без изменения спина), Ц0, Ц - полные орбитальные моменты, Яр, Цр - полный спин и орбитальный момент родительского терма. О - фактор отражает угловую зависимость матричного элемента [11].
В приложении, табл.1.-табл.4, представлены результаты расчёта сил
осцилляторов в поглощении согласно (1), (2), а также расчёты вероятностей фотопереходов. В [2] проводится сравнение с данными других авторов [14-21].
Таблицы построены в зависимости от спектроскопических символов рассматриваемых переходов (не по длинам волн). Это удобно при формировании всех дипольных переходов с выбранного нижнего уровня.
В табл.1. - табл. 4 (см. приложение) приведены значения эффективных
* / * * \
главных квантовых чисел п (и* - нижнее, и* - верхнее), в соответствии с которыми определена энергия перехода кш = Ек-Е{, эВ. Указанная длина волны Л, выраженная в ангстремах (10-10 м), соответствует энергии 1 ^ к - го перехода (в поглощении) и может в последних знаках отличаться от данных [12]. Цифры, которые записаны курсивом, взяты из [12]. Цифры длин волн, обозначенные обычным шрифтом, - наш расчёт.
Обсуждение результатов
В табл. 1 - табл.4 (см. приложение) приведены длины волн для всех разрешенных дипольных переходов рис.1 - рис.4. В эксперименте спектральные линии с длинами волн более 1000 нм практически не регистрируются. Приведённые расчётные значения сил осцилляторов во всех случаях нужны для определения средних, а также оценочных значений вероятностей фото-переходов.
Согласно [11] при переходе с уровня на конфигурацию < - фактор определяется суммой, равной 1. В этом случае суммируем значения сил осцилляторов с конкретного 1 - го нижнего уровня на все возможные в таблицах к -состояния и затем просто делим на их число.
Так, например, рассмотрим в табл.1 (см. приложение) номера переходов: 5; 13; 21;29 с уровня 7s P1/2 на всю конфигурацию 8p. Сила осциллятора равна 0.071.
Заключение
Важным этапом развития плазмомеханики является организация банка оптических констант. Представленные многочисленные радиационные константы иона ксенона необходимы при разработке кинетических моделей.
Детальное рассмотрение радиационных констант предваряет обсуждение схемы уровней атома и иона ксенона. Предложенная схема уровней иона ксенона и указанные в таблицах значения эффективного главного квантового числа позволяют рассчитывать в кулоновском приближении переходы, отсутствующие в справочном издании [12].
Сформулированное понятие «вектор состояния» [1, 2] способствует с единых позиций изучать распределения заселённостей атома и иона ксенона с использованием радиационных констант.
Библиографический список
1. Криворучко Д.Д., Скрылев А.В., Скороход Е.П. Определение концентраций возбужденных состояний и вероятностей радиационных переходов Хе1 плазмы Холловских двигателей // Труды МАИ. 2017. № 92. URL: http://trudymai.ru/published.php?ID=76859
2. Криворучко Д.Д., Кули-заде М.Е., Скороход Е.П., Скрылев А.В. Вероятности фото-переходов иона XeII и распределения возбуждённых состояний в
низкотемпературной плазме холловского двигателя // Труды МАИ. 2017. № 94. URL: http://trudymai.ru/published.php?ID=80962
3. Кули-заде М.Е., Скороход Е.П. Разработка кинетических моделей движущейся плазмы. Константы радиационных d-p переходов иона ксенона // Труды МАИ. 2017. № 95. URL: http: //trudymai.ru/published. php?ID=84423
4. Гаврилова А.Ю., Киселёв А.Г., Скороход Е.П. Диаграммы метаравновесных состояний тяжёлых инертных газов // Теплофизика высоких температур. 2014. Т. 52. № 2. С. 174 - 185.
5. Лосев С.А., Сергиевская А.Л. Каталог моделей физико-химических процессов: наполнение интернет - реализации // Материалы VII Международной конференции по неравновесным процессам в соплах и струях, Алушта, 2008. С. 286 - 288.
6. Черный Г.Г., Лосев С.А. Физико-химические процессы в газовой динамике.- М.: Изд-во МГУ, 1995, Т.1. - 343 с.
7. Черный Г.Г., Лосев С.А. Физико-химические процессы в газовой динамике. - М.: Научно-издательский центр механики, 2002. - Т. 2. - 368 с
8. Скороход Е.П. Спектроскопические методы исследования физико-химических и тепловых процессов в плазменных устройствах. Дисс. ...докт. техн. наук. - М., МАИ, 2002. - 41 с.
9. Суржиков С.Т. Сравнительный анализ роли спектральных линий атомов и ионов в радиационном нагреве поверхностей четырех типов спускаемых космических аппаратов // Теплофизика высоких температур. 2016. Т. 54. № 2. С. 249 - 266.
10. Wiese W.L., Smith M.W., Glennon B.M. Atomic Transition Probabilities, New York, Iov. Print. Opt, 1966, 153 p.
11. Собельман И.И. Введение в теорию атомных спектров. - М.: Наука, 1977. - 319 с.
12. Стриганов А.Р., Свентицкий Н.С. Таблицы спектральных линий нейтральных и ионизированных атомов. - М.: Атомиздат, 1966. - 899 с.
13. Гаврилова А.Ю., Скороход Е.П. Сечения и константы скоростей плазмохимичекских реакций инертных газов. - М.: Изд-во МАИ, 2011. -192 с.
14. Tharwat M. El-Sherbini. Transition probabilities and radiative lifetimes for singly ionized xenon // Journal of Physics B: Atomic and Molecular Physics, 1976, vol. 9, no. 10, pp. 1665.
15. Wiese W.L., Martin G.A. Wavelengths and Transition Probabilities for Atoms and Atomic Ions, Part II: Transition Probabilities, United States National Bureau of Standards NSRDS-NBS 68, 1980, 148 p.
16. Michael W. Winter, Christoph Eichhorn, Monika Auweter-Kurtz and Thomas Pfrommer, Status on Plasma Diagnostic Measurements on a RIT-10 Ion Thruster IEPC-2007-173 // 30th International Electric Propulsion Conference, Stuttgart, Germany, September 17-20, 2007. URL: http://erps.spacegrant.org/uploads/images/images/iepc articledownload 1988-2007/2007index/IEPC-2007- 173.pdf
17. Gigosos M., Mar S., Perez C., de la Rosa I. Experimental Stark widths and shifts and transition probabilities of several XeII lines // Physical Review, 1994, vol. E 49, no. 2, pp. 1575.
18. Manzella D.H. Stationary Plasma Thruster Plume Emissions. Presented as IEPC 93097 // Proceeding of the 23rd International Electric Propulsion Conference. September 1993.
19. Miller M.H., Roig R.A. Transition Probabilities of Xe I and Xe II // Physical Review A, 1973, vol. 8, pp. 480 - 486.
20. Konjevic N., Lesage A., Fuhr J.R. and Wiese W.L. Experimental Stark Widths and Shifts for Spectral Lines of Neutral and Ionized Atoms (A Critical Review of Selected Data for the Period 1989 Through 2000) // Journal of Physical and Chemical Reference, 2002, vol. 31, no. 3, URL: https://www.nist.gov/sites/default/files/documents/srd/ipcrd622.pdf
21. Di Rocco H.O., Iriarte D.I., Pomarico J.A. Lifetimes and transition probabilities of XeII: Experimental measurements and theoretical calculations // European Physical Journal D, 2000, vol.10 (1), pp. 19 - 26.
22. Шахатов В.А. Уровневые полуэмпирические столкновительно-излучательные модели в оптической диагностике неравновесных газовых разрядов: Автореферат дисс.....докт. ф-м. наук. - М.: МГУ, 2016. - 50 с.
23. Шибкова Л.В. Физические процессы в движущейся плазме многокомпонентных инертных и химически активных смесей: Автореферат дисс. докт. ф-м. наук. - М.: Объединенный институт высоких температур РАН, 2007. - 43 с.
24. Дьячков Л.Г. Развитие квазиклассической теории радиационных свойств низкотемпературной атомарной плазмы: Автореферат дисс. ... докт. ф-м. наук. - М.: Объединенный институт высоких температур РАН, 2001. - 36 с.
25. Гидаспов В.Ю. Численное моделирование стационарных волн горения и детонации в смеси частиц бора с воздухом // Труды МАИ. 2016. № 91. URL: http: //trudymai. ru/published. php?ID=75562
26. Приложение к статье «Разработка кинетических моделей движущейся плазмы. Коэффициенты Эйнштейна для иона ксенона». URL: http://tahir.sytes.net/!books/articles/Приложение к статье.pdf