Баллистическая модель движения электронов над потенциальным барьером. Преамбула. (

advertisement
Баллистическая модель движения электронов над потенциальным барьером.
Станислав Ордин (stas_ordin@mail.ru)
Преамбула.
Как уже не раз писал в своих статьях и комментариях (например, комментарий к
выступлению Жореса Алферова на пленарном заседании Государственной Думы), меня
интересуют принципиальные моменты и в самом «теле» науки и в её деятельности –
взаимодействии с окружающей средой. И так как принципиально новые моменты в море
современных «научных» публикаций встречал всё реже и реже, я лично отказался от
потери времени в виде одного рабочего дня в неделю, проводимого в молодости в БАНе.
Но мой коллега и соавтор Андрей Иванович Шалых, даже перевалив за 80 лет, не изменил
свои правила жизни и периодически, поработав в БАНе, с воодушевлением вручал мне
очередную ссылку на свежую публикацию по вопросам, сопряженным с моими
исследованиями. И каждый раз, после ознакомления со ссылками из области
термоэлектричества, я видел одну и ту же ошибку, которую отразил в своей рецензии для
«Review of Scientific Instruments» на манускрипт «Examining thermal transport through a
frequency-domain representation of time-domain thermoreflectance data»
Author: Kimberlee Collins, Alex Maznev, John Cuffe, Keith Nelson, and Gang Chen
[Manuscript #A142252]:
I.e. my assessment of this work, in principle, the positive. Work is performed, in general, at the
modern high level, "but" without special claims for fundamental nature, and without the exit
claims for a framework used in modern physics of representations.
Noted "but", is defined not so much by the level of this work, how many level of the majority of
modern scientific publications.
The used basic macroscopic representations fragmentary, are not fully considered in the analysis.
In particular, in the interpretation used by authors of the manuscript communication of heat
conductivity with thermoEMF is not considered:
1. S.V. Ordin, Thermoelectric Waves in Anisotropic Cristal of Higher Manganese Silicide, Proc.
XV the Int. Conf. on Thermoelectrics, p.212-214, (ICT'95), St. - Petersburg, Russia, 1995.
Besides, basic macroscopic representations in the majority of works are transferred on nano scale without the bases, sufficient on that, and it is often simply wrong.
Many my scientific publications earlier devoted to thermoelectric effects on micro and nano level, for example were devoted to this problem:
2. Ordin S.V., Zjuzin A. Yu., Ivanov Yu. and Yamaguchi S., Nano-structured materials for
thermoelectric devices, NATO Workshop "Advanced Materials and Technologies
for Micro/Nano-Devices, Sensors and Actuators"
June 29 - July 2, 2009, St Petersburg, Russia.
3. Ordin S.V., Zjuzin A. Yu., Ivanov Yu. and Yamaguchi S., Nano-structured materials for
thermoelectric devices, ITC2010, Shanghai, China, July, 2010, 10 pp.
4. S.V. Ordin, W.N. Wang, "Thermoelectric Effects on Micro and Nano Level.", J. Advances in
Energy Research, Volume 9, 2011, 30 pp.
To this problem it is devoted much and my popular scientific articles, for example:
1. A nano or new mentality - http://www.rusnor.org/pubs/articles/7556.htm
2. Thermodynamics and nano - http://www.rusnor.org/pubs/articles/7565.htm
3. Size effects and NANO - http://www.rusnor.org/pubs/articles/7724.htm
4. NANO and new Dimensional - http://www.rusnor.org/pubs/articles/8751.htm
This "but", in my opinion, does not belittle in any way value and usefulness of this manuscript to
which, I will repeat I give appreciation. Simply to both authors and edition I suggest not to be
afraid to see things more widely. It concerns not only this magazine. I glanced, in connection
with preparation of the review of the site MRS "Symposium BB: Nanoscale Heat Transport--From Fundamentals to Devices" and there the same, listed above by me.
Я уже и раньше отмечал, что принципиальная ошибка заключается в том, что
ультрасовременные эксперименты с нано-образцами и прецизионные их расчёты строятся
на базе макроскопических представлений. Этот вывод я сделал после анализа ряда
базовых моделей, о чём неоднократно докладывал и на семинарах, и на конференциях. Но
после ухода из жизни корифея термоэлектричества Лазаря Соломоновича Стильбанса,
термоэлектрическое сообщество фактически лишилось головы и стало невосприимчиво ко
всякого рода новым веяниям, затрагивающим базовые модели. Но эти базовые модели
были и разработаны для макроскопики, и справедливы лишь для макроскопики. Чтобы
выйти на уровень корректного описания термоэлектрических процессов на микро/наномасштабах потребовался скрупулёзный перерасчёт базовых коэффициентов и
закономерностей данного масштаба. Начало такого расчёта представлено в данной работе.
Линейная фенология.
Макроскопические процессы описываются в рамках линейной фенологии, связывающей
термодинамические потоки Qi с обобщёнными термодинамическими силами F j градиентами потенциалов  j . Для электрического E , температурного T и
концентрационного N поля в канонической форме обобщенные силы имеют следующий
вид:



T
N
FE   E , FT   2   T , FN   2   N
T
T
T
T
T
При этом в общей фенологии описываются как прямые эффекты, типа
электропроводность, теплопроводность и диффузия, являющиеся коэффициентами
линейной связи потока с градиентом сопряжённого потенциала: Qi  Lii  Fi ( i  j ), так и
перекрёстные эффекты, являющиеся коэффициентами линейной связи потока с
градиентами других потенциалов: Qi  Lij  Fj ( i  j ).
Коэффициенты линейной связи Lij называются кинетическими коэффициентами и
соответствуют коэффициентам трения различных солитонов: электронных, фононных,
концентрационных и т. д. и т.п. Макроскопические термоэлектрические процессы, по
определению перекрёстные, исторически описываются с помощью двух коэффициентов:
Зеебека и Пельтье:  и  , которые, как показал Онзагер – симметрийно вырождены.
В рамках квантово-механических представлений разработаны модели (содержащие много
предположений), которые позволили рассчитать (весьма грубо) макроскопические
кинетические коэффициенты. На масштабах меньше длины свободного пробега эти
модели, естественно, в принципе, и теоретически не работают, и не описывают результаты
экспериментов. В данной работе, оставаясь в рамках разработанных для
макроскопических кристаллов кантово-механических моделей, были проведены расчёты,
учитывающие баллистические эффекты, естественно, за рамками приближения времени
релаксации, и получена общая связь потоков и потенциалов для прямых эффектов.
Электропроводность.
Для расчёта электропроводности электронов над симметричным прямоугольным
потенциальным барьером, длина которого меньше длины свободного пробега электрона,
используем закон дисперсии свободных электронов. Использованная прямоугольная
форма барьера в дальнейшем, конечно, потребует уточнения, т.к. уже на микро-масштабе
необходимо учитывать, что реальные барьеры, как показали расчёты на базе потенциалов
атомов, далеки от идеального. Но прямоугольный барьер был использован во многих
ранее проведённых расчётах и, поэтому, качественные отличия баллистического тока от
диффузного он позволяет наглядно продемонстрировать. При приложении к барьеру
тянущего электроны электрического поля справа, его правый край и дно электронной
зоны опускаются на величину, равную разнице потенциалов на барьере (рис.2.1а).
Рис.2.1. Электронные переходы над барьером под действием электрического поля при
рассмотрении в рамках: а - стандартного подхода, б – в рамках модифицированного,
учитывающего полярность закона дисперсии, подхода.
При этом, для расчётов потока электронов над барьером обычно используют формулу
Ричардсона-Дешмана. Её квантово-механический расчёт основан на предположении
выхода в область над барьером всех электронов имеющих положительную компоненту
скорости в направлении границы барьера. При этом берётся интеграл от половины объёма
электронной зоны, изображённой на границе барьера. В одномерном QR1D и трёхмерном
случаях QR3D поток Ричардсона без учёта высоты барьера  равен
 Vx 2
b
b
2 ,
0
2
2
2
   Vx Vy Vz 
b 2
b
   e
Vx dVz dVy dVx 
2
0 
QR1D  e
QR3 D
Vx dVx 
QR3 D QR1D  b
,
где b 
2kT
m
А с учётом высоты барьера  поток Ричардсона равен
1D
R
Q
 Vx 2
b
 e
0
 
b 
3D
e kT Vx dVx   e kT , QR     e
2
0 



 Vx 2 Vy 2 Vz 2

b
e


kT
b2  kT
Vx dVz dVy dVx 
e
2
Полученный ричардсоновский ток насыщения J R  e  QR , естественно, не зависит от
приложенного электрического поля E . Для получения зависимости тока от поля берут
разницу токов насыщения Ричардсона с двух обкладок барьера с учётом высоты
потенциального барьера:


 eU
*
J R (U )  J R 1  e kBT   J R 1  e  eU  eU *  J R


Таким образом, получают линейную зависимость тока при напряжениях eU  kBT , что
соответствует независящей от приложенного напряжения электропроводности:
J (U ) e  J R
 R


U
kT

1D
e2

m
и

3D
e2 2kT


m m
При мысленном изначальном, в отсутствии электрического поля, эксперименте
Ричардсона бессмысленно рассматривать электроны с отрицательными компонентами
скорости относительно направления в сторону границы барьера. Но при изначальном
учёте в рассмотрении электрического поля нет никаких разумных оснований их не
учитывать. Более того, как будет показано ниже, для получения выражений для тока над
барьером, работающих как при малых перепадах потенциала на нано- и микро-барьерах,
так и при больших, это необходимо.
Воспользуемся методом, аналогичным расчёту поглощения при оптических переходах.
Но, в отличии от обычно рассматриваемого в оптических переходах изотропного
воздействия света, для кинетических явлений на микро/нано уровне требуется учесть
полярность (относительно внешнего электрического поля) кинетической энергии
mVx2
электронов  C*  Sign[Vx ] 
. Тогда законы дисперсии электронов на границах барьеров
2
трансформируются как показано на рис. 2.1б, а электронные переходы как для
положительных скоростей, так и для отрицательных скоростей отображаются
горизонтальными линиями. При этом, как показано на рисунке, из модифицированного
закона сохранения энергии:
mVx2 [ x]
mVx2 [0]
 C* [ x]   C* [0]  eEx
Sign[Vx [ x]] 
 Sign[Vx [0]] 
 eEx

2
2
следует, что изменение (приращение) скорости электронов сильно зависит от их
начальной скорости V0 . Их конечную скорость над барьером Vx получим, решив
уравнение для модифицированного закона сохранения энергии:
2eEx 2kT Ex

e
 b  eU * .
Vx  Sign Sign V0 V02  A Abs Sign V0 V02  A , где A 
m
m
kT
Рис.3. Зависимости конечной скорости электронов Vx (а) и её приращения Vx (б) от их
начальной тепловой скорости V0 . Разность потенциалов меняется в 2 раза.
Так как при нулевой разнице потенциалов интеграл по всем начальным скоростям (на
рис.3а – прямая линия) равен нулю, то общее выражение для суммарного тока электронов
можно получить, проинтегрировав приращения скорости (рис.3б).
С учетом больцмановского распределения электронов по энергиям для одномерного
случая
Sign Sign V V  beU    Abs Sign V V
Получим Q  eU   Q    eU  , где

1D
Q
 e
Vx 2
b
2
x
*
x
x
x

1D
 1D
*
1D
R
1D
2


 beU *   Vx dVx
*


   3 
 1
*
*
 2 HypergeometricU   , 0, eU   MeijerG  {},   , 0,1 ,{} , eU  
 2


   2 


*


eU
*
*
2
 e 2   BesselI 0, eU   BesselI 1, eU   eU *



 2 


2







Зависимость полученного полного баллистического потока от приведённого напряжения
на барьере показана на рис.4. Как видно из рисунка, этот поток существенно нелинейный
при малых напряжениях и превышает ричардсоновский.
QD
QR
6
5
4
3
2
1
eU
0.0
0.5
1.0
1.5
2.0
2.5
3.0
Рис.4. Зависимости разностного ричардсоновского потока (нижняя кривая) и полного
потока (верхняя кривая) в приведённых к потоку Ричардсона единицах от приведённой
разности потенциалов на барьере.
Качественно данное превышение суммарного микроскопического тока над
ричардсоновским (полученным с использованием как принято при макроскопическом
подходе средних значений) связано с тем, что при eU * 1 основной вклад даёт слой
низкоэнергетических электронов, который имеет и максимальную концентрацию, и
максимальное приращение скорости, а при eU *  1 - то, что насыщения тока не
происходит.
Q
QR
20.0
15.0
10.0
7.0
5.0
10 5
10 4
eU
0.001
0.01
0.1
1
10
Рис.5. Зависимость отношения суммарного тока к току Ричардсона от эффективного
напряжения на барьере.
Отношение полученных потоков в двойном логарифмическом масштабе (рис.5)
*
демонстрирует, что при eU  kT имеется минимум и неограниченное возрастание этого
отношения при напряжениях больше/меньше kT .
Полученное отношение микроскопического тока и ричардсоновского сохраняется и для
3х-мерного случая,
  
J
3D
x
3

 m   ( 2 kT ( vx 02 vy 02 vz 02 ) kT )
 e    vx  vx 0 , x  2   e
dvz 0dvy 0dvx 0
h
 

m
Проведя интегрирование, получим аналитическое выражение для зависимости тока от
приложенного напряжения (рис.10), которое является обобщённой и исправленной
формулой для разностного тока Ричардсона:
J x3 D  J R  3xD eU * 
где:
1
2
1
2
3
2
3xD eU *    ( HypergeometricU [ , 0, eU * ]  MeijerG[{{},{ }},{{0,1},{}}, eU * ]
1  eU
 e 2
2
*
eU *
eU *
]  BesselI [1,
])eU * )
2
2

3

eU *

1
1 2,0
1
eU *
eU *  

3D
*
*
*
*
2
 x eU    U ( , 0, eU )  G (eU 2 ) 
 e eU  I 0 (
)  I1 (
) 
2
2 1,2
2
2
2  


0,1


 ( BesselI [0,
- hypergeometric function, G - Meijer Gamma function, I - modified Bessel function.
При равенстве температур электродов, учитывая изменение вероятности перехода
электронов за счет их ускорения можно получить более строгое выражение  3D
x .
U
1
2
 3xD eU *    HypergeometricU [ , 0, eU * ] 

3
(4MeijerG[{{},{ }},{{0,1},{}}, eU * ]
8
2
*
eU
3
1  eU *
eU *
*
e MeijerG[{{},{ }},{{0,1},{}}, 2eU ])  e  ( 1  2e 2 ( BesselI [0,
]
2
4
2
eU *
1
 BesselI [1,
]))eU *  ( BesselK [0, eU * ]  BesselK [1, eU * ])eU *
2
4
eU *
Рис.10. Зависимость тока симметричного барьера от напряжения при равенстве
температур на границах барьера: 1 - разностный ток Ричардсона без учета
динамической добавки, 2 – с учетом динамической добавки, 3 – с учетом изменения
вероятности перехода , 4 – корневая зависимость.
Обобщение на прямые эффекты.
Как известно, диффузность макроскопических кинетических коэффициентов проявляется
в их взаимозависимости. В частности, существует соотношение Эйнштейна,
устанавливающее связь между подвижностью носителей тока u и коэффициентом их
диффузии De : u De  e kT и дающее взаимно-однозначную связь между независящими
от электрического поля электропроводностью и электронной теплопроводность.
Полученное выражение для микроскопического баллистического тока даёт нелинейную,
но взаимно-однозначную связь тока и градиента электрического потенциала. Полученное
выражение позволяет получить общую прямую связь между любыми градиентами
потенциала и соответствующими ему термодинамическими потоками, если градиент
электрического потенциала на границах барьера  E заменить в выражении
 E
eEx FE x
eU * 


kT
kT
kT 2
на градиенты других потенциалов, а поток взять соответствующий выбранному
потенциалу. Для качественного понимания потока концентрации электронов эта замена не
вызывает трудностей, т.к. он собственно и получен Q1D  J e 
  
Q1D  eU *   Q1D   N2  , а больцмановское распределение концентрации носителей
 kT 
наглядно демонстрирует разницу концентраций на границах барьеров при приложении
электрического поля на всех уровнях энергии. Но разница концентраций на границах
барьера при всех уровнях энергии может быть создана и за счёт разницы температуры на
этих границах, и аналогично ричардсоновской формуле от напряжения на барьере


 eU
*
J R (U )  J R 1  e kBT   J R 1  e  eU  eU *  J R


можно получить ричардсоноскую линейную зависимость от перепада температуры,
  
которой соответствует полученная нелинейная зависимость Q1D  eU *   Q1D   T2  .
 kT 
Справедливость данного соотношения для малых перепадов температуры качественно
также прозрачна.
Но в практическом плане проще воспользоваться известным алгоритмом для разностного
тока Ричардсона, заменив его на обобщённую формулу Ричардсона.
Перекрёстные баллистические эффекты.
Полученная обобщённая формула Ричардсона позволяет рассчитать, в первом
приближении, локальные термо-ЭДС, возникающие на симметричном барьере. Для этого
достаточно решить уравнение
 e  (U ) 
Q1D 
0
 k (T  T ) 


 1

eU
 2 HypergeometricU   , 0,




 2 k (T  T ) 


   3 


eU
0
 MeijerG  {},   , 0,1 ,{} ,



2m
k (T  T ) 
   2 


*
  eU2






eU
eU

e

BesselI
0,

BesselI
1,
eU



 2k (T  T ) 
 2k (T  T )  







Мощности имеющегося в моём распоряжении компьютера оказалось не достаточно для
аналитического решения этого уравнения и для одномерного случая. Попытки же
численного решения этого уравнения при задании пробных параметров выявили дребезг
 e  (U ) 
величины Q1D 
 при прохождении нуля по траектории U (T ) .
 k (T  T ) 
Возможные причины этого дребезга – строгость полученного аналитического выражения
для тока или строгость вычисления его значений при прохождении нуля – разделить не
получилось по той же причине: недостаточная мощность компьютера.
Download