УДК 537.538
А.С. МУСТАФАЕВ, д-р физ.-мат. наук, профессор, [email protected] Санкт-Петербургский государственный горный институт (технический университет)
A.S. MUSTAFAEV, Dr. inрhys. and math., professor, [email protected] Saint Petersburg State Mining Institute (Technical University)
МЕТОД ДИАГНОСТИКИ ЭМИССИОННЫХ СВОЙСТВ ПОВЕРХНОСТИ ТЕРМОЭМИССИОННЫХ КАТОДОВ
Представлен метод диагностики эмиссионных свойств поверхности катодов с помощью поперечного магнитного поля. В качестве одномерной модели, удобной для проведения теоретических и экспериментальных исследований, рассматривается кнудсеновский Cs-Ba-диод с поверхностной ионизацией (КДПИ).
Ключевые слова: поверхностная ионизация, работа выхода катода, коэффициент отражения электронов, кнудсеновская плазма, термоэлектронная эмиссия.
DIAGNOSTIC METHOD FOR MEASUREMENTS OF EMISSION PARAMETERS OF THERMOEMISSION CATHODES
This article deals with the diagnostic method of emission parameters of thermo emission cathodes by the transverse magnetic field. The knudsen Cs-Ba-diode with the surface ionization are investigated.
Key words: surface ionization, work function of cathodes, coefficients of reflection of heat electrons, Knudsen plasma, thermo electronic emission.
Идеальный кнудсеновский Cs-Ba-диод с поверхностной ионизацией - одномерная модель приэлектродной плазмы. КДПИ является удобной моделью изучения сильнонеравновесной плазмы с поверхностной ионизацией. В кнудсеновском режиме плазма синтезируется на поверхности катода. Ионы образуются за счет поверхностной ионизации атомов плазмообразующего компонента (цезия), электроны - за счет термоэлектронной эмиссии, при этом работа выхода катода регулируется изменением давления паров бария так, что длина свободного пробега электронов определяется только давлением паров цезия. Анод, поглощающий заряженные частицы, располагается параллельно катоду на расстоянии, меньшем длины свободного пробега электронов. Предполагается, что катод и анод однородны по работе выхода и обладают ну-
левым коэффициентом отражения от поверхности. Эта достаточно простая одномерная модель КДПИ допускает проведение теоретических и экспериментальных исследований в широком диапазоне изменения степени компенсации у. Здесь у = п1+/пв+ - отношение концентраций ионов и электронов, рожденных на поверхности катода.
Бинарный Cs-Ba-наполнитель КДПИ решает одну из важных проблем катодной электроники, снимая ограничения эмиссионной способности при низких температурах катода и давлении паров цезия.
Силы, действующие на заряженные частицы в плазме, обусловлены самосогласованным электрическим полем. Независимое изменение параметра Кнудсена и степени компенсации удобно для проведения исследования особенностей распределения
_ 83
Санкт-Петербург. 2010
потенциала, условий возникновения двойных слоев в приэлектродной плазме и их влияния на кинетику процессов токопере-носа. Варьируя давление паров бария можно в широких пределах изменять концентрации заряженных частиц, напряженность электрического поля в приэлектродном слое, направленную скорость электронов в плазме.
Однако малые пространственные размеры приэлектродных областей, высокие концентрации заряженных частиц, функции распределения которых, как правило, анизотропны, сильные электрические поля и многообразие протекающих процессов делают практически невозможным применение традиционных методов диагностики плазмы КДПИ.
Для исследования кнудсеновской анизотропной плазмы, синтезированной на поверхности катода, в работе применяется поперечное магнитное поле, с помощью которого обычно измеряются удельный заряд, масса и скорость монокинетичных пучков заряженных частиц.
\
1 2 3 4 S
Рис. 1. Магнитная характеристика КДПИ ^(5) в перекомпенсированном режиме для различных значений у (см. таблицу)
84 _
Принципы магнитной диагностики плазмы с поверхностной ионизацией. В
КДПИ без магнитного поля в межэлектродном промежутке (за исключением узких приэлектродных областей) потенциал в плазме постоянен и электроны движутся прямолинейно. На прикатодном скачке потенциала электроны приобретают нормальную к поверхности катода составляющую скорости (разброс по скоростям определяется температурой катода). В поперечном магнитном поле электроны движутся по криволинейным траекториям с ларморовским радиусом кривизны. С ростом напряженности магнитного поля ларморовский радиус становится меньше межэлектродного зазора и часть электронов, не достигнув анода, возвращается на катод. Это приводит к уменьшению анодного тока электронов. Исследование зависимости относительного ослабления электронного тока насыщения в магнитном поле (магнитной характеристики) положено в основу метода диагностики. Пространственное перераспределение концентрации электронов под действием магнитного поля приводит, в свою очередь, к изменению распределения потенциала в межэлектродном зазоре и последующему изменению траекторий движения электронов. Таким образом, поведение электронов описывается самосогласованной системой стационарных кинетических уравнений Власова с учетом внешнего магнитного поля и уравнением Пуассона. Теоретический анализ поведения электронов в магнитном поле [2-3] показал, что для определения доли электронов, возвращенных магнитным полем, достаточно знать потенциал только некоторых характерных точек в межэлектродном промежутке. Выполнена классификация распределений потенциала КДПИ в присутствии магнитного поля. С учетом этого обстоятельства задача сведена к решению трансцендентных уравнений. В результате самосогласованного решения этой системы построены магнитные характеристики ^(S) при различных значениях степени компенсации у (рис.1). Здесь
S = Jí/JQ -
ISSN 0135-3500. Записки Горного института. Т.187
Условные обозначения к рис. 1
Кривая у Кривая у Кривая у Кривая у
1 1 6 20 11 400 16 5 • 103
2 3 7 40 12 500 17 104
3 4 8 50 13 103 18 2 • 104
4 5 9 100 14 2 • 103 19 4 • 104
5 10 10 200 15 4 • 103 20 5 • 104
отношение электронного тока насыщения в магнитном поле к току в его отсутствие;
^ = еНй^^ - безразмерный параметр л/2 ткТ ё
магнитного поля; Н - напряженность магнитного поля; d - межэлектродный зазор; Тё - температура катода; е и т - заряд и масса электрона соответственно; к - постоянная Больцмана.
С увеличением давления паров бария снижается работа выхода катода и КДПИ переходит из перекомпенсированного (у > 1) в недокомпенсированный (у < 1) режим. В режиме у < 1 распределение потенциала существенно зависит от степени заполнения потенциальной ямы ионами. Теоретический анализ проведен для двух предельных условий: нулевого и полного заполнения ямы ионами. Показано, что при нулевом заполнении потенциальной ямы в недокомпенсированном режиме формируется двойной слой (виртуальный катод), перепад потенциала на внешней части которого растет с уменьшением у. Это приводит соответственно к немонотонному изменению электронного тока и направленной скорости электронов в плазме с изменением у. В результате самосогласованных расчетов [3] построены магнитные характеристики ) в недокомпенсированном режиме при различных значениях степени компенсации у и степени заполнения потенциальной ямы ионами (рис.2).
Сравнение теоретических и экспериментальных магнитных характеристик может дать ответ на вопрос о реальной степени заполнения потенциальной ямы ионами в эксперименте.
Детальные экспериментальные исследования КДПИ в магнитном поле выявили расхождения реальных магнитных характеристик с идеальными. Причины отклонения реальных магнитных характеристик, связанные с конечностью ширины приэлек-тродных областей, неоднородностью катода по работе выхода, с существованием поверхностного коэффициента отражения, были учтены в теории [1, 2] и, как следствие, привели к увеличению точности метода диагностики. Это позволило определять степень компенсации, ионный ток с поверхности катода, прикатодный скачок потенциала, потенциал и концентрацию плазмы, направленную скорость электронов в плазме, ко-
Ъ
Рис.2. Зависимость ) при нулевом заполнении потенциальной ямы ионами. Штриховая кривая соответствует полному заполнению ямы ионами
эффициенты отражения и неоднородность катода по работе выхода.
Таким образом, с помощью поперечного магнитного поля может быть выполнена диагностика анизотропной плазмы с поверхностной ионизацией.
Учет влияния реальных свойств поверхности катода. Основными особенностями идеального КДПИ являются отсутствие поверхностного коэффициента отражения и однородность электродов по работе выхода. Проведено сопоставление ионного тока, определенного при сравнении экспериментальных данных с рассчитанными магнитными характеристиками идеального КДПИ, с ионным током, вычисленным по формуле Саха - Ленгмюра. Показано, что отношение этих двух величин растет по мере увеличения степени компенсации, что является проявлением аномального эффекта Шоттки. Таким образом, реальные свойства катода настолько сильно влияют на зависимость тока насыщения от напряженности магнитного поля, что уместно ставить вопрос о восстановлении свойств катода по магнитным характеристикам. Для этого необходимо учесть в расчетах такие свойства катода, как поверхностный коэффициент отражения и неоднородность по работе выхода. Последнее особенно существенно в глубоко перекомпенсированном режиме, когда внутреннее поле пятен скомпенсировано сильным электрическим полем, возникающим в прикатодной области. Следует заметить, что при больших потенциалах анода реальные свойства анода не должны сказываться на магнитных характеристиках, так как отраженные от поверхности анода электроны рассеиваются под различными углами к поверхности анода и не могут преодолеть задерживающий при-анодный барьер. Таким образом, практически все электроны, попавшие в при-анодную область, попадают на анод.
Рассмотрим, как электроны, отраженные от поверхности катода, влияют на вид магнитных характеристик. Пусть при заданном значении 5 на анод приходит доля £ электронов. При этом из первичных элек-
тронов на катод возвращается доля 1 - £. Если г - средний коэффициент отражения от поверхности катода, то в результате отражения с катода выходит доля (1 - £) г вторичных электронов. Поскольку возвращающиеся на катод электроны имеют распределение по скоростям, близкое к мак-свелловскому с температурой катода, то при любом механизме отражения вторичные электроны будут иметь максвеллов-ское распределение. В результате воздействия магнитного поля на катод вернется доля (1 - £) всех вторичных электронов. Продолжая рассуждения дальше и суммируя вклад от всех этих электронов в проходящий ток, получим с учетом коэффициента отражения
£r = £[1 + (1-£)r + (1-£)2 r2 + ...] =
1-(1-%)r
. (1)
Отсюда видно, что при £, близких к единице, £г = £, а при £ ^ 0 £г = £/(1 - г ). Следовательно, при увеличении напряженности магнитного поля должно возрастать отклонение реальных зависимостей £г(5) от идеальных £(5). По этому отклонению можно измерить величину коэффициента отражения г . Из (1) найдем
r =
\ -Г
(1
(2)
r J
Таким образом, реальные магнитные характеристики оказываются более пологими, чем идеальные. Без учета этого обстоятельства направленная скорость электронов в плазме будет несколько завышенной.
Рассмотрим теперь, как на магнитные характеристики влияет неоднородность работы выхода катода. Пусть диод находится в сильно перекомпенсированном режиме, а электрическое поле у поверхности катода таково, что все пятна раскрыты. Прежде всего, неоднородность катода приводит к зависимости направленной скорости электронов в плазме от того, с какого участка катода электроны поступают. Следовательно, воздействие магнитного поля на эти группы будет различ-
86 _
ISSN 0135-3500. Записки Горного института. Т.187
ным. Однако то обстоятельство, что основной вклад в ток дают пятна с малой работой выхода, позволяет оценивать высоту скачка по эффективной средней работе выхода катода.
Значительно существеннее влияние неоднородности проявляется в том, что после поворота магнитным полем электроны, поступающие с участка поверхности с работой выхода ф, попадают на участок катода с некоторой работой выхода ф'. Если ф' < ф, то независимо от начальной скорости электрона он достигает поверхности катода. В противном случае (ф' > ф) электроны, обладающие малой начальной скоростью, встречают на своем пути к катоду задерживающий потенциальный барьер, который они не могут преодолеть. Отразившиеся от этого барьера электроны снова поступают в межэлектродное пространство. Таким образом, неоднородность по работе выхода катода ведет к увеличению эффективного коэффициента отражения, что, в свою очередь, влияет на магнитные характеристики. Для определения эффективного коэффициента отражения по магнитным характеристикам можно пользоваться формулой (2), заменив в ней г на г эф.
Представляет интерес связать этот эффективный коэффициент отражения с функцией распределения пятен по поверхности или с контрастностью катода. Введем следующие обозначения: g(ф^ф -вероятность встретить на поверхности катода работу выхода в интервале ф^ (ф + dф), г 0 - поверхностный коэффициент отражения.
В качестве примера рассмотрим g(ф) в виде распределения Гаусса около некоторого среднего значения ф0, т.е.
£ (ф) = С (р)е
2 2 ,-Р 2(ф-ф0)
(3)
Здесь р = 1/Аф (Аф характеризует контрастность катода); С(Р) - нормировочная константа (функция £(ф) нормирована на единицу).
Плотность тока с катода
ад
] = А(Т) { £(ф)е-аЧф.
Здесь а = а/(кТ); А(Т) - предэкспоненци-альный множитель в формуле Ричардсона. Интегрирование в формуле (4) может быть распространено до -ад, так как £ (ф) быстро
убывает, если ф стремится к нулю. Подставляя (3) в (4) и выполняя интегрирование, получим
] = А(Т )е
-а| ф0 —
4р2
Тогда эффективная работа выхода катода по электронному току
фэф = ф0
а
4р2
Рассмотрим электроны, выходящие с участка поверхности катода с работой выхода ф. В присутствии поперечного магнитного поля часть электронов возвращается к поверхности катода. Рассмотрим те из них, которые подошли к участку поверхности с работой выхода ф'. Считаем, что возвращенные магнитным полем электроны имеют полу-максвелловское распределение по скоростям. В случае больших значений напряженности магнитного поля S, когда почти все электроны возвращаются к катоду, это действительно так. При не слишком больших значениях S функция распределения возвращенных электронов не совпадает с максвелловской. Однако и в этом случае она представляет собой часть максвелловского распределения [2]. Существенно, что вырезанными оказываются области с большой энергией электронов. Поэтому поправки в расчете эффективного коэффициента отражения, связанные с отклонением функции распределения от мак-свелловской, должны быть небольшими.
Рассчитаем эффективный коэффициент отражения для данной группы электронов. Если ф' < ф, то
г ^ ф') = г0.
(5)
а
Если ф' > ф, то потенциальный барьер у катода преодолевают не все электроны и поверхности катода достигает доля электронного тока, описываемого выражением
J exP
л/ехрКф-ф')/m]
f
\
mv 2kT
2 Л
v, dv.
J exP
mv 2kT
2
ехр[-а(ф'- ф)],
vdvz
z z
J
и, следовательно,
r(ф, ф') =1 - (1 - ro)e-
а(ф'-ф)
(6)
(7)
При выводе зависимостей (6) и (7) предполагалось, что потенциальный барьер вблизи катода параллелен его поверхности.
Вероятность вылета электронов с участка поверхности ф', ф' + dф' есть g (ф)dфх х g(ф>')dty'. Интегрируя эффективный коэффициент отражения для данной группы электронов, описанной формулами (5) и (7), с учетом их вклада в полный ток и вероятности встречи такой группы, получим выражение для эффективного коэффициента отражения поверхности неоднородного катода:
ад ад
| | £(ф)£(ф'Кафг(ф, ф>ф?ф'
г _ -ад -ад_
уб ад •
| £ (ф)е-афdф
-ад
Его можно преобразовать, не используя конкретного вида £(ф) . Тогда
г.6 = 2r0 -1 +
2(1 - r0)W
(8)
где
W = J g (ф)е-а(^ф| g (фУф';
_ -аф эф
z = e
Для гауссова распределения
ВД ВД 2 2
J g(ф')¿ф' = C(Р)J e(ф'-фо) ¿ф' =
ф Ф
_1 - erf [р(ф-фо)]
и
W = z -22
и(Р) _С(Р) |eгf[p(ф-ф0)]e-aфe-p (ф-фо) dф. (9)
—ад
Введем функцию
ад 2 2
и*(0 _С(Р) |eгf(t(ф-ф0))e-aфe-p (ф-фо) dф,
дифференцируя ее по I и заменяя переменную интегрирования лучим
х = 12/(р2 +12), по-
а
(
U (р) = - — С (P)exp 2р
аф 0 +
а
2 Л
4р2
1/2 e 4P2
4Х
dx = e
_ -афуб
eгf
а
8р2
(10)
Из уравнений (8)-(10) следует, что
г _ г0 + (1 - г^г^у /422), у _ — _ — ^. уб 04 ^ Л ' 2Р кТ 2
Видно, что если у = 0, то г эф = г 0; если у ^ ад, то г эф ^ 1. На рис.3 построены зависимости г эф(у) для нескольких значений г 0.
Таким образом, величина г эф определяется двумя параметрами г 0 и у и определение г эф по магнитным характеристикам в случае сильного электрического поля у катода не позволяет однозначно восстановить оба эти параметра.
Представляет интерес построить зависимости, аналогичные рис.3, а, для случая, когда электрическое поле у катода равно нулю. Соответствующая ситуация реализуется в КДПИ при переходе из перекомпенсированного в недокомпенсированный режим. В этом случае на некотором расстоянии от катода устанавливается потенциал, равный средней работе выхода катода ф0. При этом участки катода с работой выхода, меньшей ф0, оказываются прикрытыми локальным полем пятен и дают такой же вклад в электронный ток, как и участки с работой выхода ф0.
ISSN 0135-3500. Записки Горного института. Т.187
ад
ад
X
2
а
--—х
X
z
-ад
2
Рис.3. Зависимость гэф(у) для ряда значений Г0 в случае сильного электрического поля у катода (а)
и при нулевом электрическом поле (б)
Таким образом, в случае гауссова распределения работ выхода на поверхности катода и равенства нулю внешнего электрического поля для £ (ф) имеем
g (Ф) = - 8(ф-ф о) +
+ | 0(Ф-Фс)с (ß)e-Р2(Ф-Ф0)2;
0( х) =
0, х < 0,
(11)
1, х > 0.
Плотность электронного тока в этом случае
j = A(T)Z, Z = е-аф01 Ex(y2);
<■> 2 Ex(y2) = 1 + ey erfc(y).
Подставляя (11) в (12), получим
Гэф = ro +(1 - ro)F(y),
(12)
где
■/V2
F (y) =
J Ex(y 2 - х2)dx
2 о_
л/Ж 1 + Ex( y2)
Отсюда видно, что Г(0) = 0, следовательно, г^ | ^=0 = г0. В другом предельном
случае (у ^ ад) можно использовать асимптотическое представление функции Ех(у2),
а именно Ех(у2)| у^ад я
ди
1/( ул/Ж) . Тогда
F (y)
y—
о y/4l
2 J
V/3
2
= = — arcsin
-2 ж
(1/V2) = 2.
Следовательно.
— — + —
yo | y—2 2
Зависимости y) для ряда значений Г) представлены на рис.3, б.
Измеряя с помощью поперечного магнитного поля Гэф для двух предельных случаев (нулевого и сильного электрических полей у катода) можно определить основные параметры у и /), характеризующие поверхность катода.
ЛИТЕРАТУРА
1. Мустафаев А.С. Методы диагностики анизотропной плазмы в термоэмиссионной плазменной энергетике / Автореф. дис. ... д-ра физ.-мат. наук / ОНТИ ФЭИ. Обнинск, 2003.
2.МустафаевА.С. Интегральные методы решения кинетических уравнений в кнудсеновском режиме ТЭП / А.С.Мустафаев, В.И.Бабанин, А.Я.Эндер и др.; ФТИ АН СССР. Л., 1971.
1
Г
3. Мустафаев А. С. Влияние поперечного магнитного поля в недокомпенсированном режиме ТЭП / А.С.Мустафаев, В.И.Бабанин, А.Я.Эндер и др.; ФТИ АН СССР. Л., 1971.
REFERENCES
1. Mustafaev A.S.. Diagnostics methods of anisotropic plasmas in thermionic plasma energetics: Referat
... Doctor of Phys.-math. Sci. ONTI FEI. Obninsk, 2003.
2. Mustafaev A.S. Integral methods resolving of kinetic equations / Mustafaev A.S., Babanin V.I., Ender A.Ya.; FTI AN USSR. Leningrad, 1971.
3. Mustafaev A.S. Influence of transverse magnetic fields in undercompensation regime of TEC / Mustafaev A.S., Babanin V.I., Ender A.Ya.; FTI AN USSR. Leningrad, 1971.
90 _
ISSN 0135-3500. Записки Горного института. Т.187