Радиотехника и электроника, 2020, T. 65, № 7, стр. 725-728

Природа токов, индуцированных изменениями параметров образца. Емкостные и неемкостные токи

С. Г. Дмитриев *

Фрязинский филиал Института радиотехники и электроники им. В.А. Котельникова РАН
141190 Фрязино, Московской обл., пл. Введенского, 1, Российская Федерация

* E-mail: sgd@ms.ire.rssi.ru

Поступила в редакцию 08.07.2018
После доработки 08.07.2019
Принята к публикации 28.07.2019

Полный текст (PDF)

Аннотация

В развитие идей теоремы Рамо рассмотрена природа токов во внешней цепи, возникающих при изменении параметров образца. Показано, что кроме токов, индуцированных движением зарядов в образце, и емкостных токов, возможны дополнительные. Приведены формулы и поясняющий пример, в котором при изменении параметров структуры токи во внешней цепи есть, хотя емкостные токи и индуцированные движением зарядов токи отсутствуют.

1. ТЕОРЕМА РАМО

Работа приборов в электронике определяется, как правило, движением зарядов в их активных частях. Поэтому формулы, описывающие связь между зарядами и токами в вакууме или в твердом теле (в полупроводниках и диэлектриках), с одной стороны, и токами, возникающими при этом во внешней цепи, – с другой, представляют особый интерес. Такие формулы, известные под названием теоремы Рамо (или Шокли–Рамо) [1, 2], рассматривалась в ряде работ для различных применений [110]. Впервые вклад в ток от двигающегося в вакууме одиночного заряда изучался в самом общем виде в работах [1, 2] в связи с дробовыми шумами в вакуумных электронных приборах. Эти результаты были обобщены и использованы для анализа работы вакуумных СВЧ-приборов во многих статьях и монографиях (см., например, [36] и цитированную там литературу). Применения к структурам с диэлектриками предложены в работах [79] для описания датчиков жесткого излучения [710]. Вопросы диагностики структур металл–диэлектрик–полупроводник (МДП), интегральных схем и анизотропных образцов рассматривались в работах [1115].

Отметим, что при классическом рассмотрении движения отдельных электронов в рамках теоремы Рамо на электродах согласно теории индуцируются дробные заряды. Необходимы обобщения на одноэлектронный и вообще на квантовый случаи.

В наиболее общем виде для некоторой твердотельной структуры с N металлическими электродами и протекающими в ней токами теорема Рамо выражается следующим равенством [3, 4] (см. также [15]):

(1)
$\sum\limits_{k = 1}^N {\Phi _{k}^{{(1)}}{{I}_{k}}} = \iiint {({{{\vec {E}}}^{{(1)}}} \cdot {{{\vec {j}}}_{{\text{п}}}})dV},$
где ${{\vec {j}}_{{\text{п}}}}(t,\vec {r})$ – плотность полного тока, определяемая известной формулой:
(2)
${{\vec {j}}_{{\text{п}}}} = \vec {j} + {{\partial{ \vec {D}}} \mathord{\left/ {\vphantom {{\partial{ \vec {D}}} {\partial t}}} \right. \kern-0em} {\partial t}},$
$\vec {j}(t,\vec {r})$ – плотность конвективного тока, $\vec {D}(t,\vec {r})$ – электрическая индукция, которая в отсутствие спонтанной поляризации равна
(3)
${{D}_{i}} = {{\varepsilon }_{{ij}}}{{E}_{j}},$
${{\varepsilon }_{{ij}}}(t,\vec {r})$ – тензор диэлектрической проницаемости (по одинаковым тензорным индексам здесь и далее предполагается суммирование), ${{E}_{j}}(t,\vec {r})$ – электрическое поле, ${{\vec {E}}^{{(1)}}} = - {\text{grad}}{{\varphi }^{{(1)}}}$, ${{\varphi }^{{(1)}}}(t,\vec {r})$ – соответственно электрическое поле и его потенциал, $\Phi _{k}^{{(1)}}(t)$ – потенциал k-го электрода (k = 1, 2, …, N) из некоторой вспомогательной задачи (которую мы уточним ниже), Ik – ток, втекающий в k‑й электрод из внешней цепи, определяемый равенством
(4)
${{I}_{k}} = {{\partial {{Q}_{k}}} \mathord{\left/ {\vphantom {{\partial {{Q}_{k}}} {\partial t}}} \right. \kern-0em} {\partial t}} - {{i}_{k}},$
где ik – ток, втекающий из рассматриваемой области в k-й электрод:
(5)
${{i}_{k}} = - \iint\limits_{{{S}_{k}}} {(\vec {j} \cdot \vec {n})}dV,$
Qk – заряд k-го электрода, равный
(6)
${{Q}_{k}} = \iint\limits_{{{S}_{k}}} {(\vec {D} \cdot \vec {n})}dV,$
Sk – поверхность k-го электрода, $\vec {n}$ – внешняя нормаль к ней.

Полезную формулу для тока на отдельный (α-й) электрод можно получить, выбирая вспомогательную задачу с $\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},$
где ${{\vec {E}}^{{(1\alpha )}}}$ – поле во вспомогательной задаче в рассматриваемом случае.

Обычно теорему Рамо понимают в более узком смысле, удерживая в правых частях формул (1) и (7) только конвективные токи:

(8)
${{\Phi }_{0}}I_{\alpha }^{{({\text{к}})}} = \iiint {({{{\vec {E}}}^{{(1\alpha )}}} \cdot \vec {j})dV},$
($I_{\alpha }^{{({\text{к}})}}$ – вклад во внешний ток на α-й электрод от конвективных токов в образце), или даже, как это делалось в первых работах, рассматривая движение только одного точечного заряда [1, 2].

Для вывода обсуждаемых формул можно использовать интеграл

(9)
${{J}_{1}} = - \iiint {{\text{div}}({{\varphi }^{{(1)}}}{{{\vec {j}}}_{{\text{п}}}})}dV,$
где интегрирование проводится по всему пространству, за исключением металлических электродов. Искомые формулы получаются отсюда с помощью формул векторного анализа и равенства ${\text{div}}{{\vec {j}}_{{\text{п}}}} = 0$ (предполагается также, что поверхностный интеграл по бесконечности равен нулю). Вспомогательная задача должна содержать те же электроды, что и в основной задаче, а поле в ней должно быть потенциальным (чтобы обеспечить постоянство потенциала на поверхностях электродов). Обычно в качестве вспомогательной рассматривается та же задача (в той же среде), но без зарядов и токов в образце. Временные изменения параметров вспомогательной задачи должны быть достаточно медленными, чтобы не нарушать потенциальность. В задачах диагностики, например, используются квазистационарные режимы. Поле же в исследуемой (основной) задаче может и не быть потенциальным.

Однако токи во внешней цепи индуцируются не только конвективными токами, но и (см. формулы (1) и (2)) изменениями индукции, включая изменения спонтанной поляризации, тензора диэлектрической проницаемости, включая и изменения ориентаций осей тензора (см. формулы в [11, 13, 15]). Изменения поляризации имеют ту же природу, что и конвективные токи, а токи, связанные с изменениями второй части индукции (см. (3)), относят обычно к токам емкостной природы. В книгах по электронике можно найти много конкретных примеров на этот счет (см., например, [1618]). Однако возникает вопрос, верно ли это в общем случае. Это рассматривается в следующих разделах.

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}}}},$
где Фk – потенциалы электродов в основной задаче, $E_{i}^{{(\alpha )}}$ – вспомогательное нормированное поле, $Q_{k}^{{(1\alpha )}}$ – заряд k-го электрода в вспомогательной задаче в рассматриваемом случае, $C_{k}^{\alpha }$ – коэффициенты емкости (при α = k) и коэффициенты электростатической индукции (при α ≠ k) (см., например, [19]). Мы предполагаем здесь также, что заряды и поляризация в образце во вспомогательной задаче отсутствуют, а тензор диэлектрической проницаемости симметричен (εij = εji). Более общее рассмотрение этого вопроса представлено в [15].

Отметим, что в уравнении (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$
(подробнее см. в [20]). Здесь слагаемое
(14а)
${{I}_{{\alpha 1}}} = \iiint {(E_{i}^{{(\alpha )}}{{j}_{i}})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}},$
где ε0 – диэлектрическая постоянная вакуума, а E1 > 0 и E2 > 0 – амплитуды электрических полей в первом и втором диэлектрике соответственно (поверхностная плотность заряда на границе с первым электродом равна –σ0). Тогда (как видно из (16)) потенциал второго электрода φ(d) = φ0 > 0 относительно первого электрода равен
(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}}}}.$

При этом сохраняется электронейтральность:

(19в)
${{\sigma }_{{11}}} + {{\sigma }_{{12}}} = {{\sigma }_{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)).

ЗАКЛЮЧЕНИЕ

Проанализированы условия применимости теоремы Рамо и ее обобщений. Рассмотрена природа токов, возникающих при изменении параметров образца. Показано, что кроме емкостных токов и индуцированных движением зарядов токов возможны токи иной природы. Приведен поясняющий пример с токами неемкостного характера. Дополнительные токи могут сопутствовать, например, физическим процессам в образце, которые сопровождаются изменениями диэлектрической проницаемости (структурные изменения, фазовые переходы, химические реакции и т.д.). Их изучение полезно при диагностике процессов в электронных приборах.

Список литературы

  1. Shockley W. // J. Appl. Phys. 1938. V. 9. № 10. P. 635.

  2. Ramo S. // Proc. IRE. 1939. V. 27. № 9. P. 584.

  3. Beck A.H.W. Thermionic Valves: Their Theory and Design. Cambridge: Cambridge University Press, 1953.

  4. Jen C.K. // Proc. IRE. 1941. V. 29. P. 345.

  5. Gabor D. // J. Inst. Electr. Engrs. 1944. V. 91. Pt. 3. № 15. P. 128.

  6. Гвоздовер С., Лопухин В. // Изв. АН СССР. Сер. физ. 1946. Т. 10. № 1. С. 29.

  7. Cavalleri G., Fabri G. et al. // Nucl. Instr. Meth. 1963. V. 21. P. 177.

  8. Cavalleri G., Gatti E. // Nucl. Instr. Meth. 1971. V. 92. P. 137.

  9. He Z. // Nucl. Instr. Meth. 2001. V. A463. № 1–2. P. 250.

  10. Tavernier S. Experimental Techniques in Nuclear and Particle Physics. London: Springer, 2010.

  11. Дмитриев С.Г. // ФТП. 2009. Т. 43. № 6. С. 854.

  12. Дмитриев С.Г. // ФТП. 2011. Т. 45. № 2. С. 192.

  13. Дмитpиeв C.Г. // PЭ. 2012. T. 57. № 11. C. 1229.

  14. Дмитpиeв C.Г. // PЭ. 2013. T. 58. № 9. C. 983.

  15. Дмитpиeв C.Г. // PЭ. 2018. T. 63. № 10. C. 1115.

  16. Бонч-Бруевич В.Л., Калашников С.Г. Физика полупроводников. М.: Наука, 1990.

  17. Зи С. Физика полупроводниковых приборов. М.: Мир, 1984.

  18. Nicollian E.R., Brews J.R. MOS (Metal-Oxide-Semiconductor) Physics and Technology. N.Y.: J. Wiley & Sons, 1982.

  19. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред. М.: Наука, 1982.

  20. Дмитpиeв C.Г. // PЭ. 2019. T. 64. № 9. C. 926.

Дополнительные материалы отсутствуют.