Радиотехника и электроника, 2020, T. 65, № 7, стр. 725-728
Природа токов, индуцированных изменениями параметров образца. Емкостные и неемкостные токи
С. Г. Дмитриев *
Фрязинский филиал Института радиотехники и электроники им. В.А. Котельникова РАН
141190 Фрязино, Московской обл., пл. Введенского, 1, Российская Федерация
* E-mail: sgd@ms.ire.rssi.ru
Поступила в редакцию 08.07.2018
После доработки 08.07.2019
Принята к публикации 28.07.2019
Аннотация
В развитие идей теоремы Рамо рассмотрена природа токов во внешней цепи, возникающих при изменении параметров образца. Показано, что кроме токов, индуцированных движением зарядов в образце, и емкостных токов, возможны дополнительные. Приведены формулы и поясняющий пример, в котором при изменении параметров структуры токи во внешней цепи есть, хотя емкостные токи и индуцированные движением зарядов токи отсутствуют.
1. ТЕОРЕМА РАМО
Работа приборов в электронике определяется, как правило, движением зарядов в их активных частях. Поэтому формулы, описывающие связь между зарядами и токами в вакууме или в твердом теле (в полупроводниках и диэлектриках), с одной стороны, и токами, возникающими при этом во внешней цепи, – с другой, представляют особый интерес. Такие формулы, известные под названием теоремы Рамо (или Шокли–Рамо) [1, 2], рассматривалась в ряде работ для различных применений [1–10]. Впервые вклад в ток от двигающегося в вакууме одиночного заряда изучался в самом общем виде в работах [1, 2] в связи с дробовыми шумами в вакуумных электронных приборах. Эти результаты были обобщены и использованы для анализа работы вакуумных СВЧ-приборов во многих статьях и монографиях (см., например, [3–6] и цитированную там литературу). Применения к структурам с диэлектриками предложены в работах [7–9] для описания датчиков жесткого излучения [7–10]. Вопросы диагностики структур металл–диэлектрик–полупроводник (МДП), интегральных схем и анизотропных образцов рассматривались в работах [11–15].
Отметим, что при классическом рассмотрении движения отдельных электронов в рамках теоремы Рамо на электродах согласно теории индуцируются дробные заряды. Необходимы обобщения на одноэлектронный и вообще на квантовый случаи.
В наиболее общем виде для некоторой твердотельной структуры с N металлическими электродами и протекающими в ней токами теорема Рамо выражается следующим равенством [3, 4] (см. также [15]):
(1)
$\sum\limits_{k = 1}^N {\Phi _{k}^{{(1)}}{{I}_{k}}} = \iiint {({{{\vec {E}}}^{{(1)}}} \cdot {{{\vec {j}}}_{{\text{п}}}})dV},$(2)
${{\vec {j}}_{{\text{п}}}} = \vec {j} + {{\partial{ \vec {D}}} \mathord{\left/ {\vphantom {{\partial{ \vec {D}}} {\partial t}}} \right. \kern-0em} {\partial t}},$(4)
${{I}_{k}} = {{\partial {{Q}_{k}}} \mathord{\left/ {\vphantom {{\partial {{Q}_{k}}} {\partial t}}} \right. \kern-0em} {\partial t}} - {{i}_{k}},$Полезную формулу для тока на отдельный (α-й) электрод можно получить, выбирая вспомогательную задачу с $\Phi _{k}^{{(1)}} = 0$ при k ≠ α и $\Phi _{\alpha }^{{(1)}} = {{\Phi }_{0}} = 1\,\,{\text{В}}$ [9, 10]. Тогда из (1) следует :
(7)
${{\Phi }_{0}}{{I}_{\alpha }} = \iiint {({{{\vec {E}}}^{{(1\alpha )}}} \cdot {{{\vec {j}}}_{{\text{п}}}})dV},$Обычно теорему Рамо понимают в более узком смысле, удерживая в правых частях формул (1) и (7) только конвективные токи:
(8)
${{\Phi }_{0}}I_{\alpha }^{{({\text{к}})}} = \iiint {({{{\vec {E}}}^{{(1\alpha )}}} \cdot \vec {j})dV},$Для вывода обсуждаемых формул можно использовать интеграл
где интегрирование проводится по всему пространству, за исключением металлических электродов. Искомые формулы получаются отсюда с помощью формул векторного анализа и равенства ${\text{div}}{{\vec {j}}_{{\text{п}}}} = 0$ (предполагается также, что поверхностный интеграл по бесконечности равен нулю). Вспомогательная задача должна содержать те же электроды, что и в основной задаче, а поле в ней должно быть потенциальным (чтобы обеспечить постоянство потенциала на поверхностях электродов). Обычно в качестве вспомогательной рассматривается та же задача (в той же среде), но без зарядов и токов в образце. Временные изменения параметров вспомогательной задачи должны быть достаточно медленными, чтобы не нарушать потенциальность. В задачах диагностики, например, используются квазистационарные режимы. Поле же в исследуемой (основной) задаче может и не быть потенциальным.Однако токи во внешней цепи индуцируются не только конвективными токами, но и (см. формулы (1) и (2)) изменениями индукции, включая изменения спонтанной поляризации, тензора диэлектрической проницаемости, включая и изменения ориентаций осей тензора (см. формулы в [11, 13, 15]). Изменения поляризации имеют ту же природу, что и конвективные токи, а токи, связанные с изменениями второй части индукции (см. (3)), относят обычно к токам емкостной природы. В книгах по электронике можно найти много конкретных примеров на этот счет (см., например, [16–18]). Однако возникает вопрос, верно ли это в общем случае. Это рассматривается в следующих разделах.
2. ЕМКОСТНЫЕ ТОКИ
Рассмотрим сначала вопрос о точности. Вернемся к равенствам (4)–(6), описывающим сохранение заряда в k-м электроде. Сохранение заряда здесь обеспечивается равенством нулю потока полного тока ${{\vec {j}}_{{\text{п}}}}$ через поверхность электрода (что в свою очередь, связано с равенством ${\text{div}}{{\vec {j}}_{{\text{п}}}} = 0$ и, следовательно, с обращением в нуль объемного интеграла, в который преобразуется с помощью теоремы Остроградского–Гаусса поверхностный интеграл). Следовательно, соотношение (4) выполняется строго, если поверхность интегрирования охватывает весь электрод, кроме его контактов с подводящими заряд проводами (поверхность Sk), через которые и проходит в электрод ток Ik. Но тогда не выполняется уравнение (6) (хотя неточность, конечно, мала). Это можно подправить с помощью формального приема, проводя дополнительно в интегралах с ${{\partial{ \vec {D}}} \mathord{\left/ {\vphantom {{\partial{ \vec {D}}} {\partial t}}} \right. \kern-0em} {\partial t}}$ в области контактов две совпадающие геометрически поверхности (но с противоположными внешними нормалями), дополняющие поверхности Sk до полных. Тогда в левой части уравнения (1) появятся дополнительные слагаемые, но другого вида, так как поверхности проводов неэквипотенциальны. В пренебрежении неэквипотенциальностью электроды со сквозными токами (и непроницаемыми для зарядов поверхностями) дают, очевидно, нулевой вклад в левую часть (1).
Однако уравнение (7) для тока на отдельный электрод можно получить только для уединенного электрода. В первых работах [1, 2] и была предложена модель с разъединенными электродами, в рамках которой индуцированные заряды на электродах определялись при различных положениях неподвижного точечного заряда. Токи получались путем дифференцирования координаты заряда по времени. Вопрос же о влиянии подвода заряда не обсуждался.
В связи с этим отметим, что модель с изолированными электродами можно использовать только во вспомогательной задаче (не трогая основную). Левая часть уравнения (8) и в этом случае будет не точной (из-за проводов). Но влияние проводов обычно мало. В пренебрежении эффектами такого рода можно выделить емкостные токи, определив с этой целью и сами емкости отдельных электродов. Необходимо также, чтобы поле в основной задаче было потенциальным. Тогда
(10)
${{I}_{\alpha }} = \iiint {E_{i}^{{(\alpha )}}}\left( {{{j}_{i}} + \frac{{\partial {{\varepsilon }_{{ij}}}}}{{\partial t}}{{E}_{j}}} \right)dV + \sum\limits_{k = 1}^N {C_{k}^{\alpha }} \frac{{\partial {{\Phi }_{k}}}}{{\partial t}},$(11)
$E_{i}^{{(\alpha )}} = {{E_{i}^{{(1\alpha )}}} \mathord{\left/ {\vphantom {{E_{i}^{{(1\alpha )}}} {{{\Phi }_{0}}}}} \right. \kern-0em} {{{\Phi }_{0}}}},$(12)
$C_{k}^{\alpha } = {{Q_{k}^{{(1\alpha )}}} \mathord{\left/ {\vphantom {{Q_{k}^{{(1\alpha )}}} {{{\Phi }_{0}}}}} \right. \kern-0em} {{{\Phi }_{0}}}},$Отметим, что в уравнении (10) кроме емкостного слагаемого (последнее слагаемое в правой части (10)) и индуцированных движением зарядов в образце токов (первое слагаемое под интегралом) присутствует еще одно слагаемое с производной ${{\partial {{\varepsilon }_{{ij}}}} \mathord{\left/ {\vphantom {{\partial {{\varepsilon }_{{ij}}}} {\partial t}}} \right. \kern-0em} {\partial t}}$ (второе под интегралом). Это и не удивительно: ведь изменения емкости тоже индуцируют токи, а емкость конденсатора, например, зависит от диэлектрической проницаемости. Поэтому возникает вопрос – не связано ли обсуждаемое слагаемое с изменениями емкостей. Интересующая нас здесь формула (при тех же предположениях и с теми же замечаниями, что и выше) имеет вид
(13)
${{I}_{\alpha }} = \sum\limits_{k = 1}^N {\frac{\partial }{{\partial t}}} (С_{k}^{\alpha }{{\Phi }_{k}}) + \iiint {\left\{ {E_{i}^{{(\alpha )}}{{j}_{i}} - \frac{{\partial E_{i}^{{(\alpha )}}}}{{\partial t}}{{\varepsilon }_{{ji}}}{{E}_{j}}} \right\}}dV$(14б)
${{I}_{{\alpha 2}}} = \sum\limits_{k = 1}^N {\frac{\partial }{{\partial t}}} (С_{k}^{\alpha }{{\Phi }_{k}}),$обычно тоже относят к емкостным токам (см., например, [17]), хотя токи здесь связаны не только с изменениями потенциалов электродов (традиционный вклад), но и с изменениями емкостных коэффициентов, которые зависят в том числе и от диэлектрической проницаемости (явная зависимость токов от изменений тензора диэлектрической проницаемости представлена в [15]). Если емкостные коэффициенты зависят только от потенциалов электродов, то можно ввести более привычные понятия дифференциальных емкостей $С_{{k,{\text{д}}}}^{\alpha }$, определяемых следующим равенством:
(14в)
$C_{{k,{\text{д}}}}^{\alpha } = С_{k}^{\alpha } + \sum\limits_{j = 1}^N {{{\Phi }_{j}}} \frac{{\partial С_{j}^{\alpha }}}{{\partial {{\Phi }_{k}}}}.$Тогда (14б) принимает вид
(14г)
${{I}_{{\alpha 2}}} = \sum\limits_{k = 1}^N {С_{{k,{\text{д}}}}^{\alpha }} ({{\partial {{\Phi }_{k}}} \mathord{\left/ {\vphantom {{\partial {{\Phi }_{k}}} {\partial t}}} \right. \kern-0em} {\partial t}}).$Однако, как видно из (13), дополнительное слагаемое к емкостным и индуцированным движением зарядов токам (второе слагаемое в правой части под интегралом), хотя и другого вида, все равно остается:
(15)
${{I}_{{\alpha 3}}} = \iiint {\left( { - \frac{{\partial E_{i}^{{(\alpha )}}}}{{\partial t}}{{\varepsilon }_{{ji}}}{{E}_{j}}} \right)dV}.$В следующем разделе рассматривается пример, в котором возникают такие токи (т.е. токи неемкостного характера).
3. НЕЕМКОСТНЫЕ ТОКИ
Рассмотрим одномерную структуру с двумя смежными диэлектриками с относительными диэлектрическими проницаемостями ε1 и ε2 и с двумя полубесконечными металлами, контактирующими с этими диэлектриками. Толщина диэлектриков d/2. Ось Ox направим от первого диэлектрика ко второму, начало координат поместим на контакте первого диэлектрика с металлом. Если теперь поверхностная плотность заряда на границе второго электрода σ0 > 0 при x = d, то
(16)
${{\sigma }_{0}} = {{\varepsilon }_{0}}{{\varepsilon }_{2}}{{E}_{2}} = {{\varepsilon }_{0}}{{\varepsilon }_{1}}{{E}_{1}},$(17)
${{\varphi }_{0}} = {{{{\sigma }_{0}}} \mathord{\left/ {\vphantom {{{{\sigma }_{0}}} C}} \right. \kern-0em} C},$(18)
$C = \frac{{2{{\varepsilon }_{0}}{{\varepsilon }_{1}}{{\varepsilon }_{2}}}}{{d({{\varepsilon }_{1}} + {{\varepsilon }_{2}})}}$– емкость структуры на единицу площади (потенциал первого электрода равен φ(0) = 0). Добавим далее на границу между диэлектриками (x = d/2) отрицательный заряд с поверхностной плотностью –σ1 < 0. Если потенциалы электродов при этом неизменны, то добавленный заряд индуцирует на первом и втором электродах дополнительные положительные заряды σ11 и σ12 соответственно:
(19а)
${{\sigma }_{{11}}} = \frac{{{{\varepsilon }_{1}}{{\sigma }_{1}}}}{{{{\varepsilon }_{1}} + {{\varepsilon }_{2}}}} = \frac{{{{\sigma }_{1}}Cd}}{{2{{\varepsilon }_{0}}{{\varepsilon }_{2}}}},$(19б)
${{\sigma }_{{12}}} = \frac{{{{\varepsilon }_{2}}{{\sigma }_{1}}}}{{{{\varepsilon }_{1}} + {{\varepsilon }_{2}}}} = \frac{{{{\sigma }_{1}}Cd}}{{2{{\varepsilon }_{0}}{{\varepsilon }_{1}}}}.$При этом сохраняется электронейтральность:
Изменения диэлектрических проницаемостей приводят к изменению емкости построенной структуры, которое описывается следующей формулой:
(20)
$\frac{{\delta С}}{С} = \frac{{{{\varepsilon }_{1}}{{\varepsilon }_{2}}}}{{{{\varepsilon }_{1}} + {{\varepsilon }_{2}}}}\left( {\frac{{\delta {{\varepsilon }_{1}}}}{{\varepsilon _{1}^{2}}} + \frac{{\delta {{\varepsilon }_{2}}}}{{\varepsilon _{2}^{2}}}} \right).$Отсюда следует, что емкость остается неизменной, когда
(21)
$\frac{{\delta {{\varepsilon }_{1}}}}{{\varepsilon _{1}^{2}}} = - \frac{{\delta {{\varepsilon }_{2}}}}{{\varepsilon _{2}^{2}}}.$Если это условие выполняется, то изменение плотности заряда σ12 описывается выражением
(22)
$\frac{{\delta {{\sigma }_{{12}}}}}{{{{\sigma }_{{12}}}}} = - \frac{{\delta {{\varepsilon }_{1}}}}{{{{\varepsilon }_{1}}}}.$Напомним, что при постоянных емкости С и потенциале φ0 не меняется также и плотность заряда σ0. Тогда изменение полной плотности заряда ${{\Sigma }_{2}} = {{\sigma }_{0}} + {{\sigma }_{{12}}}$ на втором электроде сводится к $\delta {{\Sigma }_{2}} = \delta {{\sigma }_{{12}}}$ и может быть отлично от нуля, если выполняется условие (21) (при этом обсуждаемое изменение заряда на электроде обеспечивается токами из внешней цепи). Но если емкость и потенциал постоянны, то емкостных токов нет. Отсутствуют здесь и токи, обусловленные движением зарядов в образце (так как нет движения).
Таким образом, поясняющий пример с неемкостными токами построен: а именно описана структура с изменяющимися параметрами, в которой емкостные токи и токи, индуцированные движением зарядов в образце, отсутствуют, а ток во внешней цепи все же есть. Результаты этого раздела можно получить и более формальным образом с помощью формул общей теории из предыдущего раздела (см. формулы (13)–(15)).
ЗАКЛЮЧЕНИЕ
Проанализированы условия применимости теоремы Рамо и ее обобщений. Рассмотрена природа токов, возникающих при изменении параметров образца. Показано, что кроме емкостных токов и индуцированных движением зарядов токов возможны токи иной природы. Приведен поясняющий пример с токами неемкостного характера. Дополнительные токи могут сопутствовать, например, физическим процессам в образце, которые сопровождаются изменениями диэлектрической проницаемости (структурные изменения, фазовые переходы, химические реакции и т.д.). Их изучение полезно при диагностике процессов в электронных приборах.
Список литературы
Shockley W. // J. Appl. Phys. 1938. V. 9. № 10. P. 635.
Ramo S. // Proc. IRE. 1939. V. 27. № 9. P. 584.
Beck A.H.W. Thermionic Valves: Their Theory and Design. Cambridge: Cambridge University Press, 1953.
Jen C.K. // Proc. IRE. 1941. V. 29. P. 345.
Gabor D. // J. Inst. Electr. Engrs. 1944. V. 91. Pt. 3. № 15. P. 128.
Гвоздовер С., Лопухин В. // Изв. АН СССР. Сер. физ. 1946. Т. 10. № 1. С. 29.
Cavalleri G., Fabri G. et al. // Nucl. Instr. Meth. 1963. V. 21. P. 177.
Cavalleri G., Gatti E. // Nucl. Instr. Meth. 1971. V. 92. P. 137.
He Z. // Nucl. Instr. Meth. 2001. V. A463. № 1–2. P. 250.
Tavernier S. Experimental Techniques in Nuclear and Particle Physics. London: Springer, 2010.
Дмитриев С.Г. // ФТП. 2009. Т. 43. № 6. С. 854.
Дмитриев С.Г. // ФТП. 2011. Т. 45. № 2. С. 192.
Дмитpиeв C.Г. // PЭ. 2012. T. 57. № 11. C. 1229.
Дмитpиeв C.Г. // PЭ. 2013. T. 58. № 9. C. 983.
Дмитpиeв C.Г. // PЭ. 2018. T. 63. № 10. C. 1115.
Бонч-Бруевич В.Л., Калашников С.Г. Физика полупроводников. М.: Наука, 1990.
Зи С. Физика полупроводниковых приборов. М.: Мир, 1984.
Nicollian E.R., Brews J.R. MOS (Metal-Oxide-Semiconductor) Physics and Technology. N.Y.: J. Wiley & Sons, 1982.
Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред. М.: Наука, 1982.
Дмитpиeв C.Г. // PЭ. 2019. T. 64. № 9. C. 926.
Дополнительные материалы отсутствуют.
Инструменты
Радиотехника и электроника