О возможности измерения работы выхода электрона труднолетучих материалов методом высокотемпературной масс-спектрометрии
- Авторы: Борщевский А.Я.1
-
Учреждения:
- Московский государственный университет имени М. В. Ломоносова
- Выпуск: Том 98, № 9 (2024)
- Страницы: 69-81
- Раздел: 100-ЛЕТИЮ ЛАБОРАТОРИИ ХИМИЧЕСКОЙ ТЕРМОДИНАМИКИ МГУ
- Статья получена: 07.03.2025
- Статья одобрена: 07.03.2025
- Статья опубликована: 30.12.2024
- URL: https://journal-vniispk.ru/0044-4537/article/view/282760
- DOI: https://doi.org/10.31857/S0044453724090101
- EDN: https://elibrary.ru/ONVGTB
- ID: 282760
Цитировать
Полный текст
Аннотация
Путем численных экспериментов изучена возможность измерения работы выхода металлов и солевых систем через измерения парциальных давлений ионной составляющей газовой фазы в рамках эффузионного метода Кнудсена в сочетании с масс-спектрометрическим анализом продуктов испарения. Проанализированы принципиальные источники возникновения ошибочных данных. Рассмотрены различные определения работы выхода электрона и связь между ними.
Полный текст
ВВЕДЕНИЕ
Определение работы выхода (РВ) электрона металлов, твердых неорганических соединений и систем на их основе – важная задача физической химии и техники, поскольку эта величина является фундаментальной характеристикой поверхности твердого тела, определяющей эмиссионную способность материала.
Понятие РВ применимо те только к электрону, но и любой заряженной частице, способной переходить из вещества через его поверхность в вакуум или иную среду. В данной работе речь идет об удалении электрона в вакуум. РВ известна для большинства химических элементов и многих соединений –оксидов, боридов, карбидов, нитридов, силицидов [1]. Нас будут интересовать исследования, выполненные в рамках эффузионного метода Кнудсена в сочетании с масс-спектрометрической регистрацией нейтральных и заряженных составляющих газовой фазы, получившего название высокотемпературная масс-спектрометрия (ВТМС) [2]. Основная цель –выяснить, в какой степени можно считать результаты измерений РВ, величинами, объективно характеризующими исследуемые вещества.
ЭКСПЕРИМЕНТАЛЬНАЯ ЧАСТЬ
Эффузионный метод
Применительно к ионам, присутствующим в насыщенном паре конденсированных веществ, метод предполагает вытягивание ионов из эффузионной ячейки (ЭЯ) электростатическим полем через небольшое отверстие, ускорение их до энергии (2–4 кэВ), необходимой для работы масс-анализатора, и измерение ионных токов в приемнике ионов. Принципиальная схема ионного источника (ИИ) и важнейшие его параметры показаны на рис. 1. Парциальные давления заряженных составляющих пара обычно на 6–10 порядков ниже давлений нейтральных составляющих.
Рис. 1. Схема ионного источника; ВЭ – вытягивающий электрод, ФЭ – фокусирующий электрод, УЭ – ускоряющий электрод, МА – масс-анализатор.
Ввиду того, что мы рассматриваем принципиальные моменты, будем считать, что ИИ работает идеально: ионы истекают только из эффузионного отверстия (ЭО) и светосила ионной оптики одинакова в режимах вытягивания ионов обоих знаков. В частности, одинакова чувствительность прибора к положительным и отрицательным ионам.
Различные определения РВ электрона
Разобраться в методах получения РВ невозможно без анализа различных определений данной величины. Общим для них является факт, что заряженная частица обладает некоторой положительной энергией связи с веществом. Это очевидно уже из того, что точечный заряд притягивается ко всякому нейтральному телу.
В ВТМС рабочие температуры T ~ 1000 K, так что любые вещества обладают некоторой проводимостью, и поэтому электрический потенциал одинаков на поверхности образца. То же относится к потенциалу в объеме. Отсюда следует, что с феноменологической точки зрения нет необходимости различать проводники, полупроводники и диэлектрики. Это уместно делать для интерпретации РВ в рамках зонной теории кристаллов.
Наиболее общее определение РВ дано в [3]: работа, которую необходимо совершить над частицей для ее удаления обратимым образом. Будучи термодинамической величиной, РВ зависит от температуры, давления, условий, в которых происходит удаление, – при постоянстве заряда тела или его потенциала. В изобарно-изотермических условиях минимальная работа, затрачиваемая при переводе 1 моль зарядов из тела в вакуум, равна увеличению мольной энергии Гиббса (∆Gm > 0). Речь идет о квазистатическом переходе электронов в состояние идеального газа под действием воображаемого силового поля [4]. Этот переход следует отличать от эмиссии электронов, когда электронный газ является компонентом насыщенного пара конденсированного тела, находится в равновесии с последним без всякого поля, и ∆Gm = 0.
Выход заряженной частицы из вещества сопряжен с преодолением потенциального барьера, возникающего благодаря существованию на поверхности двойного электрического слоя (ДЭС). Он характеризуется поверхностной плотностью дипольного момента, которая обеспечивает скачок потенциала по обе стороны поверхности. Структура ДЭС весьма чувствительна к внешним загрязнениям и зависит от кристаллографического направления. Ясно поэтому, что РВ определяется состоянием поверхности проводника. По современным представлениям РВ можно условно разделить на две части – работу по преодолению ДЭС и поляризационную, обусловленную притяжением уходящего электрона к телу [5].
Термоэлектрическая РВ (ϕTE), определяемая по температурной зависимости тока эмиссии электронов из нагретого металла с использованием известного теоретического выражения для плотности тока эмиссии (формулы Ричардсона – Дэшмана) [6]:
(1)
Такой эксперимент аналогичен процедуре определения энтальпии сублимации вещества по температурной зависимости упругости его пара.
Формула (1) приведена в современном виде, отличающимся от [6] множителем 2 перед экспонентой. Эта поправка возникает благодаря существованию у электрона спина – обстоятельству, которое не могло быть учтено Дэшманом в силу исторических причин.
В работах [7–13], выполненных эффузионным методом Кнудсена, которые нас здесь интересуют, авторы понимают РВ как энтальпию сублимации заряженной частицы (∆sHi°). Так, в [9] РВ электрона кристаллических галогенидов щелочных металлов (MX) определена опосредованно, с помощью специального термохимического цикла, в который входят измеренные температурные зависимости токов ионов M2X+ или MX2–, присутствующих в насыщенном паре соединений. Во-первых, правильнее было бы говорить о парциальной мольной энтальпии сублимации. Во-вторых, в реальных условиях ячейка с образцом соли, очевидно, не изолирована, а поддерживается под фиксированным потенциалом, поэтому потеря заряда при “сублимации” ионов определенного знака должна компенсироваться его притоком из резервуара (электротехнической земли). Только тогда имеет место аналогия с сублимацией. Нельзя не отметить, что ϕTE в (1) согласно основополагающей работе [6] относится к абсолютному нулю температуры, в отличие от ричардсоновской РВ ϕRICH, фигурирующей в формуле для j, выведенной Ричардсоном ранее. Связь между этими величинами дается равенством ϕTE = ϕRICH – (3/2) kBT/e.
Насколько нам известно, ранее определение РВ через ∆sHi°в литературе отсутствовало. Обычно РВ определяют как энергию связи электрона в кристалле или энергию, связанную с термодинамически равновесным переходом электрона из кристалла во внешнюю среду [14]. Тогда разность между энергией и энтальпией в расчете на 1 моль есть PV = RT, так что эту поправку легко учесть. При T ≈ 1000K она составляет примерно 8.3 кДж моль–1 ≈ 0.1 эВ.
Еще одно определение РВ электрона [14] – “минимальная энергия, необходимая для перемещения его с поверхности Ферми в твердом теле в точку пространства, где поле отсутствует”. Оно относится и к металлам, и полупроводникам, и диэлектрикам. Имеется в виду перевод электрона с уровня, отвечающего энергии Ферми εF, в свободное состояние в вакууме с нулевой кинетической энергией, причем переход должен быть адиабатическим, а не вертикальным. На самом деле не имеет значения, с какого именно уровня энергетической зоны удален электрон (образовавшаяся вакансия заполняется электроном с более высоких уровней), и только тогда данное определение имеет термодинамический смысл. Однако полностью придать ему этот смысл невозможно, так как εF есть химический потенциал электронного газа в кристалле при T = 0, т. е. не зависит от температуры. Другие смыслы РВ видны из рис. 2. Величину χ называют сродством к электрону [14], она могла бы служить еще одним определением РВ – для диэлектриков с достаточно широкой запрещенной зоной.
Рис. 2. К определениям РВ электрона проводников (a), полупроводников или диэлектриков (б); wFE – фотоэлектронная РВ; ∆εg – ширина запрещенной зоны; χ – энергия сродства к электрону вещества.
Таким образом, понятие РВ неоднозначно. Ради справедливости следует сказать, что различия в значениях РВ, связанные с различием определений, скорее всего, не превышают погрешностей измерения РВ известными способами, если отбросить фотоэлектронную РВ (ϕFE) как величину, не имеющую отношения к рассматриваемому вопросу.
Формула (1) выведена с использованием распределения Ферми–Дирака, и в ней ϕTE = ε0 – εF (ε0 – дно потенциальной ямы, в которой находится электрон в твердом теле), так что определения 2 и 4 следует считать идентичными. По-видимому, термоэлектрическая РВ – единственная энергетическая характеристика связи электрона с кристаллом, доступная измерению и имеющая термодинамический смысл. Далее мы будем под РВ подразумевать именно эту величину и обозначать ее посредством ϕ.
РВ связана с парциальным давлением электронов в насыщенном паре внутри замкнутой полости в веществе соотношением [15]
(2)
Давления заряженных составляющих в замкнутой полости
Обычно ЭЯ изготавливается из тугоплавких и инертных по отношению к образцу металлов (Ni, Mo, W, Pt) или соединений того же свойства. Многие авторы справедливо считают [11, 16], что в условиях ВТМС внутренняя поверхность ячейки покрыта адсорбционным слоем толщиной в несколько диаметров молекул загруженного в нее образца. Поэтому ионный состав пара определяется эмиссионными свойствами изучаемого вещества.
Возьмем в качестве вещества чистую кристаллическую соль MX(cr), где M – щелочной металл, X – галоген. Основной компонент пара – молекулы MX. Нас будут интересовать также продукты гомолитической (M, X) и гетеролитической (M+, X–) диссоциации этих молекул и, в связи с соотношением (2), электроны (e–). Эксперименты показывают [17], что помимо вышеуказанных ионов в соизмеримых количествах имеются ионы M2X+ и MX2–, учет которых для демонстрационных целей не принципиален.
Свойства газовой фазы, которая представляет собой слабо ионизированный газ, существенно зависят от соотношения длины Дебая
(3)
и характерного линейного размера L полости, содержащей газ. В формуле (3) ni – плотность зарядов сорта i (см–3), zi = ±1 – зарядовые числа, P – суммарное давление заряженных составляющих (дин см–2), e = 4.8×10–10 ед. СГСЭ. Важнейшие параметры плазмы при 1000 K приведены в табл. 1. Условие идеальности плазмы – ND >> 1.
Таблица 1. Характеристики ионизированного газа при 1000 K (zi = ±1)
n = Σni±, см–3 | P = ΣPi±, атм | rD, cм | ND | r, см | λ, см |
7.34×1013 | 10–5 | 2.55×10–5 | 5.1 | 1.48×10–5 | 9.15×10–4 |
7.34×1012 | 10–6 | 8.06×10–5 | 16.1 | 3.19×10–5 | 7.03×10–3 |
7.34×1011 | 10–7 | 2.55×10–4 | 51 | 6.88×10–5 | 5.71×10–2 |
7.34×1010 | 10–8 | 8.06×10–4 | 161 | 1.48×10–4 | 0.48 |
7.34×109 | 10–9 | 2.55×10–3 | 510 | 3.19×10–4 | 4.15 |
7.34×108 | 10–10 | 8.06×10–3 | 1610 | 6.88×10–4 | 36.5 |
7.34×107 | 10–11 | 2.55×10–2 | 5.10×103 | 1.48×10–3 | 326 |
7.34×106 | 10–12 | 8.06×10–2 | 1.61×104 | 3.19×10–3 | 2.94×103 |
7.34×105 | 10–13 | 2.55×10–1 | 5.10×104 | 6.88×10–3 | 2.68×104 |
7.34×104 | 10–14 | 0.81 | 1.61×105 | 1.48×10–2 | 2.45×105 |
7.34×103 | 10–15 | 2.55 | 5.10×105 | 3.19×10–2 | 2.28×106 |
7.34×102 | 10–16 | 8.06 | 1.61×106 | 6.88×10–2 | 2.12×107 |
7.34×101 | 10–17 | 25.5 | 5.10×106 | 1.48×10–1 | 1.98×108 |
Примечание. ; ND– среднее число заряженных частиц в объеме сферы радиуса rD; r – среднее расстояние между ионами, λ – длина свободного пробега ионов [27].
При высоком давлении rD L (в опытах L ≈ 10 мм) и внутренняя область газа нейтральна, а избыточный заряд сосредоточен в пристеночном слое толщиной ~rD. Если же rD p L (низкое давление) то нейтральность не гарантирована. Поскольку результат расчета давлений заранее не известен, его надо считать численным экспериментом, который покажет его применимость.
Низкие давления. Между интересующими нас составляющими пара возможны следующие независимые реакции:
{1} MX(cr) = MX(g), K1
{2} MX(g) = M(g) + X(g), K2
{3} M(g) = M+(g) + e–(g), K3
{4} X–(g) = X(g) + e–(g), K4
где K1, K2, K3 и K4 – константы равновесия. Значения констант для NaF и CsI для численного моделирования взяты из справочника [18] и приведены в табл. 2. Эти данные соответствуют высшему классу точности I.
Таблица 2. Константы равновесия (атм) реакций с участием галогенидов щелочных металлов
Соединение | T, K | –lgK1 | –lgK2 | –lgK3 | –lgK4 |
NaF | 800 | 9.9813 | 26.3844 | 31.5999 | 19.6445 |
900 | 8.0090 | 22.8621 | 27.8748 | 17.1299 | |
1000 | 6.4437 | 20.0417 | 24.8827 | 15.1059 | |
CsI | 800 | 5.0974 | 17.6074 | 23.7557 | 17.6224 |
900 | 3.7538 | 15.1581 | 20.9022 | 15.3500 | |
1000 | 2.8374 | 13.1968 | 18.6073 | 13.5200 |
Соли MX сублимируют конгруэнтно, поэтому парциальные давления компонентов связаны условием равенства элементных составов обеих фаз в соотношении M/X = 1/1. В результате получаем недоопределенную систему из четырех уравнений с пятью неизвестными парциальными давлениями:
P(M)P(X) = K2P(MX), (4)
P(M+)P(e–) = K3P(M), (5)
P(X)P(e–) = K4P(X–), (6)
P(M) + P(M+) = P(X) + P(X–). (7)
В обычной практике соединения MX справедливо считают однокомпонентными, однако в рассматриваемом случае, когда необходимо учитывать электроны, положение меняется. Давления P(MX) и P(e–) в соответствии с равновесием {1} и формулой (2) независимо определяются одной лишь температурой, и можно произвольно задать либо ϕ, либо P(e–). Рассматриваемая двухфазная система в реальных условиях является открытой, так как обмен заряженными частицами между газовой и твердой фазами (например, при изменении температуры) должен компенсироваться электронами из внешнего резервуара (электротехнической земли). Результаты расчета для NaF и CsI приведены в строках (a) табл. 3, 4 при каждом заданном P(e–) и отвечающей ему РВ.
Таблица 3. Парциальные давления (атм) в паре NaF(cr) при 1000 K; P(NaF) = 3.60×10–7 атм
P(e–) | ϕ, эВ | P(Na) | P(F) | P(Na+) | P(F–) | rD, мм | |
10–28 | 5.82 | (a) | 1.58×10–15 | 2.07×10–12 | 2.07×10–12 | 2.64×10–25 | 0.56 |
(b) | 6.00×10–14 | 5.45×10–14 | 7.86×10–11 | 6.96×10–27 | 0.09 | ||
(c) | 3.03×10–16 | 1.08×10–11 | 3.96×10–13 | 1.38×10–24 | 1.28 | ||
10–22 | 4.63 | (a) | 5.71×10–14 | 5.72×10–14 | 7.49×10–17 | 7.30×10–21 | 93.1 |
(b) | 6.00×10–14 | 5.45×10–14 | 7.86×10–17 | 6.96×10–21 | 90.9 | ||
(c) | 5.88×10–14 | 5.56×10–14 | 7.71×10–17 | 7.09×10–21 | 91.8 | ||
10–20 | 4.23 | (a) | 5.72×10–14 | 5.72×10–14 | 7.49×10–19 | 7.30×10–19 | 660 |
(b) | 6.00×10–14 | 5.45×10–14 | 7.86×10–19 | 6.96×10–19 | 660 | ||
(c) | 6.00×10–14 | 5.45×10–14 | 7.86×10–19 | 6.96×10–19 | 660 | ||
10–18 | 3.83 | (a) | 5.72×10–14 | 5.72×10–14 | 7.50×10–21 | 7.29×10–17 | 93.7 |
(b) | 6.11×10–14 | 5.72×10–14 | 8.01×10–21 | 6.83×10–17 | 96.8 | ||
(c) | 6.00×10–14 | 5.45×10–14 | 7.86×10–21 | 6.96×10–17 | 95.9 | ||
10–14 | 3.04 | (a) | 2.12×10–13 | 1.54×10–14 | 2.78×10–24 | 1.97×10–13 | 1.77 |
(b) | 1.18×10–12 | 2.78×10–15 | 1.54×10–23 | 3.55×10–14 | 3.78 | ||
(c) | 6.00×10–14 | 5.45×10–14 | 7.86×10–25 | 6.96×10–13 | 0.96 |
Обозначения: (a) – замкнутая полость, (b) – вытягивание отрицательных ионов, (c) – вытягивание положительных ионов.
Таблица 4. Парциальные давления (атм) в паре CsI(cr) при 800 K, P(CsI) = 7.99×10–6 атм
P(e–) | ϕ, эВ | P(Cs) | P(I) | P(Cs+) | P(I–) | rD, мм | |
10–28 | 4.62 | (a) | 3.35×10–14 | 5.89×10–10 | 5.88×10–10 | 2.47×10–20 | 2.66×10–2 |
(b) | 4.49×10–12 | 4.39×10–12 | 7.89×10–8 | 1.84×10–22 | 2.29×10–3 | ||
(c) | 6.19×10–15 | 3.19×10–9 | 1.09×10–10 | 1.34×10–19 | 6.17×10–2 | ||
10–22 | 3.67 | (a) | 4.40×10–12 | 4.48×10–12 | 7.73×10–14 | 1.88×10–16 | 2.32 |
(b) | 4.50×10–12 | 4.39×10–12 | 7.89×10–14 | 1.84×10–16 | 2.29 | ||
(c) | 3.64×10–12 | 5.43×10–12 | 6.38×10–14 | 2.27×10–16 | 2.55 | ||
10–20 | 3.35 | (a) | 4.45×10–12 | 4.43×10–12 | 7.81×10–16 | 1.86×10–14 | 4.73 |
(b) | 4.77×10–12 | 4.14×10–12 | 8.37×10–16 | 1.73×10–14 | 4.90 | ||
(c) | 4.48×10–12 | 4.40×10–12 | 7.87×10–16 | 1.85×10–14 | 4.74 | ||
10–18 | 3.03 | (a) | 5.29×10–12 | 3.73×10–12 | 9.29×10–18 | 1.56×10–12 | 0.52 |
(b) | 1.66×10–11 | 1.19×10–12 | 2.91×10–17 | 5.00×10–13 | 0.91 | ||
(c) | 4.49×10–12 | 4.39×10–12 | 7.89×10–18 | 1.84×10–12 | 0.48 | ||
10–14 | 2.39 | (a) | 2.88×10–10 | 6.86×10–14 | 5.05×10–20 | 2.88×10–10 | 0.04 |
(b) | 1.59×10–9 | 1.24×10–14 | 2.80×10–19 | 5.19×10–11 | 0.09 | ||
(c) | 4.49×10–12 | 4.39×10–12 | 7.89×10–22 | 1.84×10–8 | 0.005 |
Обозначения: (a) – замкнутая полость, (b) – вытягивание отрицательных ионов, (c) – вытягивание положительных ионов.
Высокие давления. При rD L к уравнениям (4)–(7) следует добавить условие нейтральности газовой фазы:
P(M+) = P(X–) + P(e–). (8)
Тогда система уравнений полностью определена при каждой температуре. Рассчитанные парциальные давления и некоторые другие параметры газовой фазы приведены в табл. 5, 6.
Таблица 5. Сравнение результатов расчета давлений (атм) продуктов диссоциации NaF в условиях электронейтральности пара
P(Na) | P(F) | P(Na+) | P(F–) | |
T = 800 K, PNaF = 1.04×10–10, Pe = 2.35×10–26, ϕ = 4.24 эВ, rD = 5500 см | ||||
(a) | 6.56×10–19 | 6.56×10–19 | 7.03×10–25 | 6.79×10–25 |
(b) | 6.64×10–19 | 6.49×10–19 | 7.11×10–25 | 6.72×10–25 |
(c) | 6.64×10–19 | 6.49×10–19 | 7.11×10–25 | 6.72×10–25 |
T = 900 K, PNaF = 9.79×10–9, Pe = 3.11×10–23, ϕ = 4.24 эВ, rD = 1300 см | ||||
(a) | 3.67×10–16 | 3.67×10–16 | 1.57×10–21 | 1.54×10–21 |
(b) | 3.71×10–16 | 3.63×10–16 | 1.59×10–21 | 1.52×10–21 |
(c) | 3.71×10–16 | 3.63×10–16 | 1.59×10–21 | 1.52×10–21 |
T = 1000 K, PNaF = 3.60×10–7, Pe = 1.01×10–20, ϕ = 4.23 эВ, rD = 66 см | ||||
(a) | 5.72×10–14 | 5.72×10–14 | 7.44×10–19 | 7.34×10–19 |
(b) | 5.79×10–14 | 5.65×10–14 | 7.53×10–19 | 7.26×10–19 |
(c) | 5.78×10–14 | 5.65×10–14 | 7.53×10–19 | 7.26×10–19 |
Обозначения: (a) – замкнутая полость, (b) – вытягивание отрицательных ионов, (c) – вытягивание положительных ионов.
Таблица 6. Сравнение результатов расчета давлений (атм) продуктов диссоциации CsI в условиях электронейтральности пара
P(Cs) | P(I) | P(Cs+) | P(I–) | |
T = 800 K, PCsI = 7.99×10–6, Pe = 2.05×10–21, ϕ = 3.46 эВ, rD = 7.38 мм | ||||
(a) | 4.44×10–12 | 4.44×10–12 | 3.81×10–15 | 3.81×10–15 |
(b) | 4.55×10–12 | 4.34×10–12 | 3.90×10–15 | 3.72×10–15 |
(c) | 4.44×10–12 | 4.45×10–12 | 3.81×10–15 | 3.81×10–15 |
T = 900 K, PCsI = 1.76×10–4, Pe = 7.48×10–19, ϕ = 3.45 эВ, rD = 0.67 мм | ||||
(a) | 3.50×10–10 | 3.50×10–10 | 5.86×10–13 | 5.86×10–13 |
(b) | 3.63×10–10 | 3.38×10–10 | 6.08×10–13 | 5.65×10–13 |
(c) | 3.45×10–10 | 3.55×10–10 | 5.79×10–13 | 5.94×10–13 |
T = 1000 K, PCsI = 1.45×10–3, Pe = 8.64×10–17, ϕ = 3.45 эВ, rD = 0.11 мм | ||||
(a) | 9.61×10–9 | 9.61×10–9 | 2.75×10–11 | 2.75×10–11 |
(b) | 1.01×10–8 | 9.12×10–9 | 2.90×10–11 | 2.61×10–11 |
(c) | 9.33×10–9 | 9.90×10–9 | 2.67×10–11 | 2.83×10–11 |
Обозначения: (a) – замкнутая полость, (b) – вытягивание отрицательных ионов, (c) – вытягивание положительных ионов.
Давления заряженных составляющих в условиях вытягивания ионов
Низкие давления. В условиях отбора ионов из ЭЯ с помощью электрического поля их истечение из ЭО нельзя считать эффузией, так как на выходе из отверстия заряды имеют энергию в несколько эВ, что примерно на два порядка больше тепловой энергии, отвечающей температурам ~1000 K. Этот вопрос качественно затрагивался в [8, 19, 20]. Количественная картина (для потенциала U) проникновения однородного электростатического поля (напряженностью E) в ЭЯ, рассчитанная нами по формулам для бесконечной проводящей плоскости с круглым отверстием [3], показана на рис. 3 и в табл. 7. За исходные параметры взяты величины, типичные для экспериментов.
Рис. 3. Сечения эквипотенциальных поверхностей электрического поля внутри ЭЯ; a = 0.5 мм, E = = 31.45 В мм–1. Потенциалы указаны в мВ.
Таблица 7. Зависимость потенциала внутри ЭЯ с круглым отверстием диаметром 0.5 мм от глубины проникновения поля
r, мм | U, мВ | |
θ = 0 | θ = 45° | |
0 | 5000 | 5000 |
0.5 | 1073 | 1123 |
1 | 364 | 301 |
2 | 100 | 74.3 |
3 | 45.5 | 32.9 |
4 | 25.8 | 18.5 |
5 | 16.6 | 11.8 |
6 | 11.5 | 8.19 |
7 | 8.48 | 6.02 |
8 | 6.50 | 4.61 |
9 | 5.14 | 3.64 |
10 | 4.16 | 2.95 |
Формулы хорошо описывают реальное поле в области, не слишком близко прилегающей к стенкам камеры, и при условии rD p L, когда экранированием поля пространственным зарядом можно пренебречь. В целом поле внутри ЭЯ носит характер дипольного, и U зависит от расстояния (r) от центра ЭО по закону обратного квадрата. Эквипотенциаль kBT/e = 86 мВ, отвечающая средней тепловой энергии при 1000 K, расположена в заштрихованной области II. Обозначим ее площадь через SU.
Несмотря на вышесказанное, описание истечения ионов из ЭЯ в терминах эффузии возможно благодаря тому, что поле в ячейке быстро убывает по мере удаления от ЭО, и объем ячейки можно достаточно резко разграничить на три части. В первой из них, прилегающей к стенкам (I на рис. 3), энергия направленного движения ионов много меньше тепловой энергии, газ можно считать невозмущенным полем и распределение ионов по скоростям соответствует максвелловскому. В области II это соотношение обратно. Между ними расположена весьма непротяженная “переходная” область III. Работу ЭЯ можно рассматривать как эффузию ионов из I в II через переходную зону. Тогда выражение для потока заряженных частиц будет отличаться от молекулярного эффузионного потока [2] лишь заменой Sэфф на SU, играющую роль эффективной площади отбора ионов [19]:
, (9)
где mi – масса заряда сорта i. Величина SU в отличие от Sэфф может меняться в ходе эксперимента в зависимости от вытягивающего напряжения U0.
При вытягивании ионов следует заменить уравнение (7) для давлений уравнениями для потоков. С учетом (9) имеем для случаев отбора отрицательных и положительных ионов, соответственно:
Для вычислений принята оценочная величина α = 30. Результаты расчетов представлены в табл. 3–6 в строках (b) и (c).
Высокие давления. При rD < L в плазме невозможно макроскопическое пространственное разделение разноименных зарядов. При любой интенсивности эмиссии ионов и электронов ДЭС претерпевает такое изменение потенциала, которое способствует прохождению в газовую фазу зарядов одного знака и препятствует зарядам другого знака, обеспечивая нейтральность газа во внутренней области ЭЯ. Из табл. 1 видно, при каких давлениях и концентрациях ионов начинает выполняться указанное условие.
При наличии отверстия в полости поток ионов через него определяется не только эмиссионной способностью стенок. Если имеется внешнее вытягивающее поле, то оно экранируется собственным электрическим полем плазмы, и уменьшение напряженности поля, проникающего в ЭЯ, происходит быстрее, чем в пустой ячейке. Для демонстрации характерных параметров, описывающих проникновение поля в плазму, на рис. 4а (кривая Ia) представлены результаты расчета для равновесной плазмы в однородном поле, направленном вдоль оси x, по формуле из [21]
,
где U0 – потенциал поля на границе плазмы. Для сравнения показана кривая IIa, отвечающая простому экспоненциальному закону затухания поля U = U0exp(–x/xD). В плоскости x = 0 находится непроницаемая для частиц перегородка; среднее направленное движение ионов в плазме отсутствует. За начало отсчета потенциала U принято его значение на бесконечной глубине. Величина U выражена в единицах kBT/e, x – в единицах дебаевской длины xD. На рис. 4б показано, на какую глубину проникает в плазму эквипотенциаль kBT/e при заданном потенциале внешнего поля U0 на границе (кривая Ib). В отличие от экспоненциального закона (кривая IIb) эта глубина не может быть больше 1.4xD и практически не зависит от внешнего поля при U0 > 10.
Рис. 4. Глубина проникновения электростатического поля в плазму в закрытой полости (a) и в условиях вытягивания ионов (б).
Процесс вытягивания плазмы в вакуум (случай ЭЯ) довольно хорошо описан в литературе [22] при анализе работы электрических зондов, плазменных источников ионов, и т. п. Непосредственным источником ионов является граница плазмы – поверхность, на которой поле обращается в нуль. В пространстве между границей плазмы и вытягивающим электродом происходит ее разрушение, и имеется направленный поток заряженных частиц. Плотность потока описывается выражением, не зависящим от вытягивающего напряжения U0:
,
Величиной, зависящей от U0, является толщина слоя разрушения плазмы xd. Связь между ni±, U0 и x в рассматриваемом случае дается уравнением [22]
. (10)
Учитывая, что , используя (3) и подставляя в (10), получим
. (11)
Из (11) видно, что в отличие от случая с непроницаемой перегородкой поле может проникать в плазму на сколь угодно большую глубину. Под U0 следует понимать потенциал в центре ЭО относительно ЭЯ, и он имеет иной смысл, чем вытягивающее напряжение в ВТМС. Согласно табл. 7 значение U0 близко к 5 В, что отвечает 58 в единицах kBT/e при 1000 K.
Таким образом, в ЭЯ роль SU играет граница плазмы, площадь которой, в отличие от случая малых давлений, зависит не только от U0, но и от концентрации ионов. Это означает, что коэффициент чувствительности экспериментальной установки (κ) изменяется с давлением ионов (например, при изменении температуры опыта), то есть может нарушаться линейная связь между ионными токами Ii и парциальными давлениями Pi ионов, лежащая в основе всех эффузионных экспериментов [20]:
В силу особенностей конфигурации ЭЯ отбор ионов из различных ее областей имеет разный характер. Вблизи краев ЭО плазма полностью разрушена, и отбор ионов происходит как с открытой поверхности, а из удаленных от ЭО областей – как из плазмы. Отметим, что это давно обнаружено при изучении электронной эмиссии полых катодов [23].
ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ
В настоящей работе расчеты состава газовой фазы проведены без каких-либо предположений за исключением химического равновесия. В замкнутой полости этот вопрос не возникает. В ЭЯ одно из условий равновесия состоит в выполнении неравенства λ >> dэфф (см. табл. 1) для нейтральных составляющих пара. Для заряженных составляющих dэфф следует заменить характерным линейным размером области III (рис. 3), а в случае высоких давлений – границы плазмы. Свойства этой области III определены таким образом, что за ее внешними пределами максвелловское распределение ионов по скоростям не нарушается. Отметим, что принятая в расчетах величина α = SU/Sэфф = 30 является верхней оценкой; линейный размер области III вряд ли превышает dэфф более чем в 5 раз. Если форму эквипотенциали или границы плазмы считать сферой, то α ≈ 25. Таким образом, химическое равновесие не нарушается по указанным причинам.
Другая причина, которая может препятствовать достижению равновесия в полости с ЭО – низкая ионная проводимость полимолекулярного адсорбционного слоя образца (NaF и CsI в нашем случае), поверхность которого эмитирует заряженные частицы. В работе [10] убедительно показано, что этот фактор не оказывает заметного влияния.
Обращает на себя внимание, что вытягивание ионов оказывает сильное влияние на отношения давлений P(M)/P(X) и P(M+)/P(X–) только тогда, когда сами эти отношения либо очень велики, либо очень малы. Это имеет место при слишком низких или слишком высоких давлениях электронов. Например, для NaF при 1000 K и P(e–) = 10–28 атм. давление Na+ выше, чем F– на 13 порядков (см. табл. 3). При этом нарушается нейтральность пара, хотя радиус rD мал по сравнению с размером полости, что должно было ее обеспечивать. При “высоком” P(e–) = 10–14 атм. картина прямо противоположная. Зато для промежуточного значения P(e–) = 10–20 атм. указанные выше отношения близки к единице, и нейтральность соблюдается, несмотря на большой дебаевский радиус (66 см). Неудивительно поэтому, что расчет всех давлений сразу (в том числе электронов), предполагающий наложение условия нейтральности, дает тот же самый результат (см. табл. 5). Получается, что это требование соблюдается автоматически, и конгруэнтная сублимация NaF имеет место не только в обычном смысле, но и означает практически одинаковую интенсивность эмиссии ионов разного знака. Аналогичное поведение наблюдается и в случае CsI (см. табл. 4, 6).
Другой важный результат состоит в том, что согласно табл. 5, 6 для всех случаев РВ оказывается одинаковой с точностью до 0.01 эВ при различных температурах и давлениях P(e–), отличающихся на три порядка. Строго говоря, РВ любого эмиттера зависит от температуры и обычно ее аппроксимируют линейной зависимостью ϕ(T) = ϕ(0) + aT. По данным [24] средний температурный коэффициент различных эмиттеров лежит в пределах a = ±(10–3–10–6) эВ×K–1. Для NaF и CsI коэффициент a, разумеется, не известен, но при ∆T = 100–200 K эффект мог бы быть заметен. Покажем, однако, что его отсутствие является следствием связи (2) между давлением электронов, температурой и РВ. Перепишем (2) в виде
(12)
Выражение Φ можно рассматривать как функцию трех аргументов, один из которых (ϕ) сам является функцией остальных двух. Вычислим так называемую полную частную производную DΦ/DT по температуре:
Но второй член в (12) есть производная ∂ϕ/∂T, так что DΦ/DT = 0. Несмотря на это, можно утверждать, что вычисленная РВ, если и отличается от истинного значения, то не сильно.
Проблема измерения РВ электрона эффузионным методом состоит в возможном искажении этой величины по сравнению со значением, получаемом путем измерения тока эмиссии в условиях вакуума по методике Ричардсона–Дэшмана. Для этого имеются две фундаментальные причины: 1) Наличие плотности заряда вблизи внутренней поверхности ЭЯ (если rD ‘ L), что приводит к появлению дополнительного энергетического барьера, не имеющего отношения к свойствам ДЭС вещества как такового; 2) Адсорбция продуктов диссоциации молекул (M, X, M+, X–), меняющая структуру ДЭС за счет отличия P(M)/P(X) и P(M+)/P(X–) от единицы в объеме газовой фазы. Если в первом случае можно оценить возникающий потенциал, то во втором такая возможность не просматривается. Влияние обоих факторов будет минимальным, если измерения проводятся в нейтральном паре, что, очевидно, не может быть гарантированно во всех опытах. Рассмотрим с этой точки зрения результаты вычислений, проделанных в настоящей работе.
Из табл. 3 видно, что неравенство rD p 10 мм хорошо соблюдается в диапазоне P(e–) = 10–22–10–18 атм. В остальных случаях rD d 10 мм, и давления ионов разных знаков различаются на много порядков, то есть отсутствует нейтральность пара. При низком давлении электронов (10–28 атм.) общее давление ионов определяется катионами Na+, а при высоком (10–14 атм) – анионами F–. Весьма значительно, хотя и в меньшей степени, различаются давления Na и F. В промежуточном случае (10–20 атм) давления Na+ и F– близки, хотя rD (66 см) сильно превышает размер полости. Очевидно, именно давление электронов, близкое данному, нивелирует дисбаланс между энергией ионизации IE(Na) (5.14 эВ) и сродством к электрону атома EA(F) (3.45 эВ). Соответственно, близкими оказываются давления продуктов гомолитической диссоциации NaF по реакции {2}. Поскольку избыточный заряд в объеме полости практически отсутствует, можно надеяться, что полученная РВ (4.23 эВ) близка в величине, определяемой внутренними свойствами соли. Более точное давление электронов 1.01×10–20 атм (табл. 5), полученное принудительным наложением условия (8), совпадает с приведенным выше, так что избыточный заряд отсутствует несмотря на большой радиус Дебая. Неудивительно поэтому, что указанное выше значение ϕ(NaF) совпадает с полусуммой [IE(Na) + EA(F)]/2 = (5.14 + 3.45)/2 = 4.23 эВ – факт, отмеченный в [13], хотя там почему-то приведено значение 4.27 эВ.
Таким образом, в случае NaF (табл. 5) давления заряженных составляющих оказываются столь низкими, что rD превышает размер ЭЯ на 2–4 порядка. Замечательно, что нейтральность пара при определенном Pe имеет место без требования (8). Это дает основания надеяться, что достаточно точные измерения парциальных давлений эффузионным методом привели бы к величине ϕ образца, определяемой самим кристаллом.
Похожая картина получается для CsI, но при более низкой температуре 800 K (см. табл. 6). При 1000 K и даже 900 K рассмотрение не имеет смысла, так как упругость пара данного вещества выходит за пределы применимости эффузионного метода. Так, P(CsI, 900 K) = 1.76×10–4 атм. Отличие от NaF состоит в том, что даже при самой низкой из взятых температур (800 K) rD ~ 10 мм, т. е. реализуется промежуточный случай и наложение условия (8) не вполне правомерно.
Ясно, что в таких простых случаях, как бинарные соединения MX, для нахождения РВ никаких особенных масс-спектрометрических измерений делать не надо, так как фактическая нейтральность пара приведет к значению ϕ = [IE(M) + EA(X)]/2. Однако при более богатом ионном и молекулярном составе пара на достоверные результаты можно рассчитывать лишь тогда, когда опыт позволяет убедиться в той степени нейтральности газовой фазы, которая обеспечивает отсутствие некомпенсированного заряда вблизи поверхности образца.
Что касается фактора адсорбции, то оценка искажений РВ за счет него не поддается количественному анализу. Полученные в настоящей работе парциальные давления это лишь теоретический расчет, по которому можно определить тенденцию возможных изменений ϕ. Так, из табл. 5, 6 следует, что вытягивание ионов в условиях, близких к нейтральности пара, и ϕ = const, слабо влияет на давления. Но даже весьма незначительное нарушение нейтральности может привести к изменению эффективной РВ. В [7] отмечается, что вытягивание отрицательных ионов способствует увеличению давлений ионов указанного знака, а не уменьшению, как показывает расчет равновесия. Это означает (в нашем примере), что в реальности поверхность откликается на избыток M и M+ понижением РВ. Вытягивание положительных ионов дает обратный эффект. В ВТМС возможны только последовательные измерения ионных токов разноименных зарядов, поэтому результаты будут относиться к поверхностям с разной структурой ДЭС. В [8] проведен разбор влияния давления щелочного металла на давления ионов на примере системы NaF–AlF3 с помощью модельных представлений. Рассматривая ионный кристалл как полупроводник с широкой запрещенной зоной, слабо допированный Na, и применяя закон Генри, автор установил, что давление ионов AlF4–не зависит от давления Na. Роль щелочного металла сводится только к понижению РВ поверхности изучаемой солевой системы. Тем самым он подтвердил, что получить значение ϕ для чистого кристалла в данном случае невозможно.
Для пущей наглядности возможностей эффузионного метода стоит рассмотреть некоторые другие результаты, имеющиеся в литературе. В работах [12, 13] использован любопытный способ расчета парциального давления электронов и других заряженных составляющих в насыщенном паре вольфрама и некоторых галогенидов щелочных металлов, основанный на принципе экстремума общего давления при конгруэнтной сублимации. Результаты для этого второго случая представлены в табл. 8. В первом случае металл рассматривается как бинарная система W+–e–, в которой имеется “соединение” W с узкой областью гомогенности. Парциальные давления Pi определяются из уравнений, получающихся приравниванием к нулю производной ∂P/∂Pi от полного давления P по Pi.
Таблица 8. Рассчитанные парциальные давления (атм) составляющих насыщенного пара некоторых галогенидов щелочных металлов по литературным данным
MX | T, K | P(M) = P(X) | P(M+) | P(e–) | ϕ0,*эВ | Литература |
NaF | 1100 | 3.5×10–12 | 1.1×10–16 | 1.2×10–18 | 4.27 | [13] |
CsI | 800 | 4.4×10–12 | 6.8×10–15 | 3.0×10–21 | 3.48 | [13] |
* ϕ0 = [IE(M) + EA(X)]/2 – полусумма энергии (IE) ионизации и сродства к электрону (EA)
Хорошо известно, что указанный принцип применим к любым двухфазным системам с любым числом компонентов. Он утверждает равенство составов сосуществующих фаз при экстремуме температуры или давления. В случае вольфрама, обозначая мольную долю электронов в металле через x, а в паре – через y, должно быть x = y. Предложенный в [12, 13] метод представляется эквивалентным использованию обобщенного уравнения Ван-дер-Ваальса, которое не требует громоздких выкладок и, главное, дополнительных предположений для получения численных результатов. В самом деле, запишем для данного случая это уравнение в переменных газовой фазы [25], положив в нем T = const (m указывает на молярные значения объема и энергии Гиббса):
(13)
Учитывая, что Pe = yP, Vm(g) = RT/P, Gm(g) = (1 – y)µ(g)(W+) + yµ(g)(e–), можно переписать (13) в виде
. (14)
Производная (14) равна нулю, когда x = Pe/P, но это и означает, что составы газовой и твердой фаз совпадают, поскольку конгруэнтная сублимация вольфрама по утверждению авторов [12] предполагает равенство x(W+, cr.) = x(e–, cr.). Тогда без всяких дополнительных предположений парциальные давления заряженных составляющих можно вычислить таким же простым способом, как это сделано в настоящей работе для солей MX, записав соответствующие равновесия и решив простейшую систему уравнений:
{5} W(cr) = W(g), P(W) = K5 = 2.477×10–15, (15)
{6} W(g) = W+(g) + e–(g), P(W+)P(e–)/P(W) = K6 = 2.477×10–15, (16)
{7} W–(g) = W(g) + e–(g), P(W)P(e–)/P(W–) = K7 = 0.8624, (17)
P(W+) = P(e–). (18)
Константы равновесия реакций {5} и {6} рассчитаны нами по данным [18] для T = 2000 K по формуле
,
где Φ° – приведенный термодинамический потенциал.
В систему включена реакция {7}, так как в [12] учтено присутствие в паре аниона W–. Данные о термодинамических функциях W–(g) в [18] отсутствуют. При практически одинаковых массах W и W– отличие Φ°(W–, g) от Φ°(W, g) сводится к учету статистических весов основного состояния (6S0 и 5D0, соответственно). Эта поправка равна Rln(6/5). Сродство к электрону EA(W) = 0.816 эВ = 78.73 кДж моль–1 взято из [26] как наиболее надежное. Собранные таким образом данные для реакций {5}–{7} приведены в табл. 9.
Таблица 9. Термодинамические данные реакций в газовой фазе вольфрама
Реакция | ∆rΦ°, Дж×моль–1×K–1 | lgK | K |
{5}W(cr) = W(g) | 145.42 | –14.606 | 2.48×10–15 |
{6}W(g) = W+(g) +e–(g) | 41.47 | –17.944 | 1.14×10–18 |
{7}W–(g) = W(g) +e–(g) | 38.13 | –0.062 | 0.866 |
Решая систему (15)–(18), находим P(W+) = = P(e–) = 5.31×10–17 атм., P(W–) = 1.52×10–31 атм. (см. табл. 10, случай 1) В этой же таблице представлены результаты расчетов при других предположениях (случаи 2 и 3), заменяющих собой условие (18). Близкие значения давлений получены в [12], но в приближении K5/K7 << 1, которое понадобилось авторам после применения уравнения ∂P/∂Pe = 0. Давлению электронов отвечает РВ ϕ(W) = 7.3 эВ (в работе [12] почему-то 8.6 эВ), а общее давление заряженных частиц дает rD = 320 мм. В любом случае ϕ резко отличается от значения 4.55 эВ, рекомендованного во многих справочных изданиях, и полученного измерениями тока эмиссии, то есть методом Ричардсона–Дэшмана. Объяснение авторов [12] следующее: “рассчитанная величина относится к области, отстоящей от поверхности металла на расстояние, превышающее rD”. Иначе говоря, к области пространства, где должна соблюдаться электронейтральность, а в значение ϕ автоматически включена потенциальная энергия, связанная с преодолением разности потенциалов, создаваемой областью с отличной от нуля общей плотностью заряда. Утверждение совершенно верное, хотя в наш расчет условие электронейтральности не закладывалось, а ее фактическое наличие связано с ничтожным давлением P(W–). Но rD на полтора порядка превышает размер ЭЯ, значит, измерения парциальных давлений эффузионным методом привели бы к совершенно неверному результату для РВ. Стоит заметить, что плотность заряда в ЭЯ пренебрежимо мала, а значительная разность потенциалов возникает за счет гигантской протяженности “приповерхностного” слоя.
Таблица 10. Парциальные давления (атм) в паре вольфрама при 2000 K
P(W) | P(W+) | P(W–) | P(e–) | ϕ, эВ | rD, мм | Лит. |
Случай 1: P(W+) = P(e–) | ||||||
2.47×10–15 | 5.31×10–17 | 1.52×10–31 | 5.31×10–17 | 7.29 | 156 | |
2.5×10–15 | 6.8×10–17 | 1.1×10–27 | 6.8×10–17 | 8.6 | 320* | [12] |
Случай 2: P(W+) = P(W–) | ||||||
2.47×10–15 | 2.84×10–24 | 2.84×10–24 | 9.92×10–10 | 4.40 | 0.05 | |
2.8×10–22 | 2.8×10–22 | 1.7×10–11 | 4.40 | [12] | ||
Случай 3: P(W+) = P(W–) + P(e–) | ||||||
2.47×10–15 | 5.31×10–17 | 1.52×10–31 | 5.31×10–17 | 7.29 | 156 |
* Рассчитано неверно, без учета давлений всех заряженных составляющих.
ЗАКЛЮЧЕНИЕ
Результаты измерения РВ электрона эффузионным методом могут быть достоверными только тогда, когда эксперимент позволяет убедиться в электронейтральности газовой фазы. Это условие необходимое, но не достаточное, так как извлечение ионов разного знака электрическим полем может значительно влиять на РВ в противоположных направлениях. Но даже в самых благоприятных случаях полученные данные вряд ли заслуживают внесения в реестр справочных данных, поскольку не являются свойствами соединений как таковых.
Работа выполнена при поддержке проекта “Химическая термодинамика и теоретическое материаловедение (№ 121031300039-1). Автор выражает благодарность к. х. н. М. И. Никитину за полезные обсуждения.
Об авторах
А. Я. Борщевский
Московский государственный университет имени М. В. Ломоносова
Автор, ответственный за переписку.
Email: andrey.borschevsky@gmail.com
химический факультет
Россия, МоскваСписок литературы
- Фоменко В.С. Эмиссионные свойства материалов. Справочник. Киев: Наукова думка, 1981.
- Сидоров Л.Н., Коробов М.В., Журавлева Л.В. Масс-спектральные термодинамические исследования. М.: Изд-во МГУ, 1985.
- Ландау Л.Д, Лифшиц Е.М. Электродинамика сплошных сред. М.: Физматлит, 2005.
- Леонтович М.А. Введение в термодинамику. Статистическая физика. М.: Наука, 1983.
- Владимиров В.В. Физическая электроника. Часть 1. Термоэлектронная эмиссия. М.: Физматлит, 2005.
- Dushman S. // Phys. Rev. 1923. V. 21. № 6. P. 623.
- Sidorov L.N. // Int. J. Mass Spectrom. Ion Phys. 1981. V. 38. № 1. P. 49.
- Никитин М.И. // Теплофизика высоких температур. 1986. Т. 24. № 2. С. 247.
- Кудин Л.С., Дунаев А.М., Бутман М.Ф. и др. // Журн.физ.химии. 2011. Т. 85. № 2. С. 313.
- Никитин М.И., Збежнева С.Г. // Теплофизика высоких температур. 2014. Т. 52. № 2. С. 229.
- Никитин М.И., Збежнева С.Г. // Там же. 2014. Т. 52. № 4. С. 556.
- Никитин М.И., Чилингаров Н.С., Осина Е.Л. и др. // Там же. 2018. Т. 56. № 4. С. 1.
- Осина Е.Л., Никитин М.И., Скокан Е.В. и др. // Там же. 2023. Т. 61. № 1. С. 30.
- Зарубский А. А, Чернова А.П. Структура и электронные свойства твердых тел. СПб., 2007.
- Зандберг Э.Я., Ионов Н.И. Поверхностная ионизация. М.: Наука, 1969.
- Бутман М.Ф. Масс-спектрометрическое исследование ионной компоненты в парах гидроксидов щелочных металлов методом ионно-молекулярных равновесий. Дисс. … канд. хим. наук. Иваново: ИХТИ. 1985.
- Кудин Л.С., Гусаров А.В., Горохов Л.Н. // Теплофизика высоких температур. 1973. Т. 11. № 1. С. 59.
- Термодинамические свойства индивидуальных веществ: Справочное издание в 4-х т. / Под ред. В.П. Глушко. М.: Наука, 1978–1984.
- Sidorov L.N., Sorokin I.D., Nikitin M.I. et al. // Int. J. Mass Spectrom. Ion Phys. 1981. V. 39. № 1. P. 311.
- Sidorov L.N., Zhuravleva I.D., Sorokin I.D. // Mass Spectrom. Rev. 1986. V. 5. No 1. P. 73.
- Смирнов Б.М. Физика слабоионизованного газа. М.: Наука, 1978.
- Габович М.Д. Физика и техника плазменных источников ионов. М.: Атомиздат, 1972.
- Poole K.M. // J. Appl. Phys. 1955. V. 26. No 9. P. 1176.
- Kawano H., Page F.M. // Int. J. Mass Spectrom. Ion. Phys. 1983. V. 50. P. 1.
- Морачевский А.Г., Смирнова Н.А., Пиотровская Е.М. и др. Термодинамика равновесия жидкость–пар. Л.: Химия, 1989.
- Смирнов Б.М. Физика атома и иона. М.: Энергоатомиздат, 1986.
- Арцимович Л.А., Лукьянов С.Ю. Движение заряженных частиц в электрических и магнитных полях. М.: Наука, 1972.
Дополнительные файлы
