close

Вход

Забыли?

вход по аккаунту

?

Вычисление распределений локальных токов и полей в слоистых средах в условиях квантового эффекта Холла.

код для вставкиСкачать
Физика, математика, техника, технология
УДК 537.7
ББК В. 22
М. Ю. Малакеева
г. Чита, Россия
В. Е. Архинчеев
г. Улан-Удэ, Россия
Вычисление распределений локальных токов и полей
в слоистых средах в условиях квантового эффекта Холла
Рассмотрено протекание тока в условиях квантового эффекта Холла, исследован вопрос о формировании значения эффективной холловской проводимости в средах с границами. Установлено, что в этом случае эффективная холловская проводимость принимает
ненулевое значение. Показано, что оно обусловлено протеканием холловского тока через
конечное число сингулярных точек (в нашей модели это углы на стыке фаз).
Методом конформных преобразований получены локальные распределения токов и
полей в слоистых средах со смещением полос в шахматном порядке.
Ключевые слова: квантовый эффект Холла (КЭХ), эффективная проводимость, конформное преобразование, условие Коши-Римана, слоистые среды.
M. Yu. Malakeeva
Chita, Russia
V. Ye. Arkhincheev
Ulan-Ude, Russia
Calculation of Local Current and Field Division in the Stratified Medium under
Quantum Hall Effect Conditions
The current percolation has been considered in the medium with boundaries under
Quantum Hall Effect conditions. It has been shown that in that case the effective Hall
conductivity possessed a non-zero value due to percolation of the Hall current through the
finite number of singular points (in our model these are corners at the phase joints).
Within the conformal transformation method, local distributions of current and field were
found in layered media, in which regions are shifted in a staggered order.
Keywords: Quantum Hall Effect (QHE), effective conductivity, conformal mapping,
Cauchy-Riemann conditions, layered media.
1. Введение
При исследовании протекания тока в средах с границами одной из основных проблем является проблема нахождения эффективной проводимости. В двумерном случае благодаря дуальной
симметрии уравнений постоянного тока были получены точные результаты для эффективной проводимости. К ним относятся теорема Келлера [13, с. 548] и общий подход, развитый в работах Дыхне
[5, с. 111]. Было показано, что эффективная проводимость двухфазной среды с проводимостями фаз
и на пороге протекания (при равных концентрациях фаз) равна:
σt (0) =
√
σ1 σ2 .
(1)
При произвольных концентрациях фаз было установлено дуальное соотношение для эффективной проводимости:
σt (ǫ) σe (−ǫ) = σ1 σ2 ,
108
(2)
© Малакеева М. Ю., Архинчеев В. Е., 2011
Ученые записки ЗабГГПУ
где ǫ = x − xc – отклонение от порога протекания xc = 1/2. В работах [6; 7] методом конформных преобразований были найдены также локальные распределения токов и полей в двоякопериодических структурах. Было также показано, что формула (2) имеет смысл стационарной точки
преобразований Дыхне [1, с. 951].
В общем случае были установлены нижняя и верхняя границы для эффективной проводимости:
<
1 −1
> ≤ σe ≤< σ > .
σ
(3)
Здесь < α > означает среднее значение величины α . Для установления граничных значений
для эффективной проводимости было рассмотрено выражение для джоулева тепла [4, с. 265]:
Z 1
~ E
~ = σe E~2 .
~j, ē dV = J,
Q=
(4)
V
Если подставить в подынтегральное выражение формулы (4) вместо локального электрического тока его среднее значение, то получим нижнюю границу в формуле (3). И наоборот, если вместо
локального электрического поля подставить среднее значение для поля, то получим верхнюю границу в формуле (3). Уточнения полученных границ связаны с дальнейшими предположениями о
симметрии задач ( подробнее см [10]).
Однако вышеизложенные рассуждения не применимы для исследования вопроса о локальных
распределениях токов и полей в неоднородной двумерной среде в условиях квантового эффекта
Холла: σxx = 0, σxy = const. Как известно в этом случае электрический ток всегда перпендикулярен
электрическому полю:
~j = σxe [~n~e] .
(5)
Здесь ~n – единичный вектор, направленный перпендикулярно двумерной плоскости по магнитному полю. Поэтому джоулево тепло при протекании тока в режиме КЭХ всегда равно нулю.
Q = 0.
(6)
Т. е. холловские фазы являются бездиссипативными. Следовательно, вопрос о значении холловской проводимости в условиях КЭХ не может быть исследован описанным выше способом, поскольку необходимо учитывать особенности протекания тока в условиях квантового эффекта Холла.
Особенности протекания тока в режиме квантового эффекта Холла исследовались во многих
работах как экспериментально, так и теоретически, причем исследования проводились в том числе
и с помощью макроскопических уравнений электростатики на основе метода Дыхне [2; 12], см также
[14; 11].
Важно отметить еще одну характерную особенность протекания тока в режиме квантового
эффекта Холла. Дело в том, что из обычных граничных условий непрерывности тангенциальных
компонент электрического поля e1t = e2t и нормальных компонент электрического тока j1n = j2n
при протекании тока в режиме КЭХ вытекают новые граничные условия:
j1n = j2n = 0.
(7)
Другими словами, холловский ток не может пересекать границу раздела фаз и всегда обтекает
неоднородности за исключением конечного числа сингулярных точек. Поэтому на первый взгляд
кажется, что холловская проводимость в средах с границами раздела фаз всегда должна обращаться в ноль из-за невозможности протекания через границы фаз. Как будет показано в работе,
проводимость отлична от нуля и соответственно возникает вопрос о механизме протекания тока и
формировании холловской проводимости.
Целью настоящей работы является исследование протекания тока в условиях КЭХ, вычисление
локальных распределений токов и полей методом конформных преобразований.
2. Методы теории функции комплексного переменного
в двумерных задачах проводимости
Чтобы понять, как локальные поля и токи распределены в слоистых средах с межфазными
границами в условиях квантового эффекта Холла, рассмотрим модель, состоящую из слоистых
сред со смещением слоев относительно друг друга в начале координат (рис. 1).
109
Физика, математика, техника, технология
Рис. 1. Слоистые среды со смещением слоев в начале координат
Для нахождения локальных распределений токов (полей) в слоистых средах со смещением
полос в шахматном порядке в начале координат будут использованы методы теории функции комплексного переменного. Для этой цели переформулируем проблему проводимости в терминах теории
функции комплексных переменных [7; 8; 9].
Перейдем в плоскость комплексного переменного z = x + iy и введем комплексные выражения
для двумерных электрических токов и электрических полей:
(k)
j (k) (z) = jx(k) − ijy(k) , e(k) (z) = e(k)
x − iey
jx =
j + ~j
j − ~j
; jy =
.
2
2i
(8)
Здесь k = 1, 2. Такое представление величин возможно для функций, удовлетворяющих условиям Коши-Римана. Легко проверить выполнение этих условий для введенных выше функций электрических токов и полей. Действительно, одно из условий Коши-Римана дает уравнение неразрывности тока:
(k)
(k)
∂jy
∂jx
=−
,
(9)
∂x
∂y
а второе условие следует из уравнения rot~e = 0, равносильного для сред с кусочно-непрерывной
проводимостью уравнению rotj̄ = 0 :
(k)
(k)
∂jx
∂jy
=
.
∂y
∂x
(10)
Аналогично уравнения rot~e = 0 div~e = 0, второе из которых следует из уравнения div~j = 0,
приводят к условиям Коши-Римана для компонент электрического поля:
(k)
(k)
(k)
(k)
∂ex
∂ey ∂ex
∂ey
=
;
=−
.
∂y
∂x
∂x
∂y
Вид этих условий объясняет наличие отрицательного знака перед мнимыми частями функций
j(z) и e(z). Закон Ома в магнитном поле также может быть представлен в комплексной форме:
j (z) =
σk
E (z)
1 − iβk
(11)
В последнем соотношении учитывается эффект Холла. Степень влияния магнитного поля на
распределение электрических токов в сплошной среде определяется величиной параметра Холла β,
равного отношению циклотронной частоты ω к частоте столкновений электронов ν: β = ω/ν. Параметр Холла βk так же, как и электропроводность σk , изменяясь дискретно, принимает различные
значения (k=1, 2) в смежных ячейках.
110
Ученые записки ЗабГГПУ
3. Конформное отображение плоскости Z на плоскость ξ
Вследствие имеющейся геометрической симметрии во всех одноцветных ячейках картина поля
двоякопериодически повторяется, поэтому достаточно рассмотреть элементарную ячейку, состоящую из двух прилегающих полубесконечных полос с различными проводимостями. Пусть на рис. 1
это будут ячейки abcd и abc′ d′ (bc = bc′ = L) Ниже мы построим конформное отображение внутренних областей прилегающих полуполос с проводимостями и на нижнюю и верхнюю полуплоскости.
В нашем случае такое отображение осуществляется функцией:
πz ξ = cosh
.
(12)
L
Здесь cosh (x) – гиперболический косинус, L – ширина полосы.
После отображения полуполоса abcd перейдет в нижнюю, а полуполоса abc′ d′ – в верхнюю
полуплоскость с разрезом (−∞, 1) на оси Oξ. Соответствие точек отображения показано на рис. 2,
а вычисления координат отображенных точек приведено в приложении 1.
Рис. 2. Соответствие точек при конформом отображении
4. Решения граничной задачи в матричном виде
Граничные условия также сохраняются и на плоскости ζ = ξ + iη для электрического тока
(k)
(k)
j (k) (ζ) = jξ (ξ, η) − ijη (ξ, η). На различных частях оси Oξ они имеют различный вид:
j1 − j¯1 = 0, j2 − j¯2 = 0,
(13)
j1 + j¯1 = 0, j2 + j¯2 = 0,
(14)
на L1 , включающее в себя лучи [ba) , [cd) , [c′ d′ ). Черта над функцией указывает на комплексное
сопряжение.
′
на L2 , включающее в себя лучи [bc] , [bc ].
Таким образом, исходные начальные граничные условия сведены к граничной задаче Римана
для кусочно – непрерывной функции j (ζ). Коэффициенты граничной задачи терпят скачки на
оси Oξ в точках b, c, c′ . Как следует из граничных условий в задаче имеются четыре неизвестные
функции: jk , j¯k , k = 1, 2. В соответствие с [7, с. 164] представим функции j2 , j¯1 в виде функций от
их комплексно сопряженных величин j1 , j¯2 . Введем кусочно-аналитическую вектор-функцию:
Φ1 (ζ)
j (ζ)
Φ (ξ) =
= ¯2
(15)
j1 (ζ)
Φ2 (ζ)
+ − j (ξ)
j (ξ)
+
−
На оси ξ она принимает предельные значения Φ (ξ) =
; Φ (ξ) =
.
j − (ξ)
j + (ξ)
Введенная функция удовлетворяет условию симметрии:
0
1
Φ̄ ζ̄ =
Φ (ζ) .
(16)
1
0
Соответственно, на оси Oξ эта функция принимает значения:
j (ξ)
Φ+ (ξ) = ¯2
j1 (ξ)
(17)
111
Физика, математика, техника, технология
и
Φ− (ξ) =
j1 (ξ)
.
j¯2 (ξ)
(18)
Таким образом, вышеуказанные граничные условия формулируют векторно-матричную краевую задачу Римана. Кратко она записывается в виде:
Φ+ (ξ) = Gi Φ− = (ξ) ,
(19)
σ2 /σ1
µ (ξ) Q
где коэффициент Gi есть матрица Gi = P1
.
−µ (ξ) Q
1
В исследуемом нами случае протекания тока в режиме КЭХ эти матрицы имеют простой вид.
На интервале L1
0 1
G1 =
,
(20)
1 0
а на интервале L2
G2 =
0 −1
.
−1 0
(21)
Векторно-матричная краевая задача Римана (19) в общем виде не имеет решения, поэтому
для решения необходимо провести соответствующую диагонализацию матриц граничных условий
и привести их к стандартным скалярным краевым задачам. Для этого введем новую кусочнонепрерывную аналитическую функцию:
(
Φ (ξ) ,
если Im (ζ) > 0,
Ψ (ξ) =
(22)
M Φ (ξ) , если Im (ζ) < 0.
Здесь М-матрица, совпадающая с матрицей Gi на отрезках −∞ < ξ < −1 и 0 < ξ < 1 оси Oξ.
В результате преобразований для этой функции получим следующие граничные условия:
Ψ+ (ξ) = Ψ− (ξ) на интервале L1 .
(23)
Это означает, что на этом интервале функция аналитически продолжима из одной плоскости
в другую. На остальной части действительной оси из интервала L2 получим:
Ψ+ (ξ) = G2 M −1 Ψ− (ξ).
(24)
Дальнейшее упрощение задачи можно получить редукцией ее к задаче с матрицей диагонального типа. С этой целью введем еще одну функцию равенством: Φ(ξ) = SF (ξ).
Здесь матрица S подбирается так, чтобы диагонализовать матрицу границ условий на участке
L2 . В результате получим систему из двух независимых граничных задач Римана.
Ψ1 (ξ) = λ1 F1 (ξ),
Ψ2 (ξ) = λ2 F2 (ξ).
(25)
Согласно алгоритмудиагонализации
матриц вычислим собственные значения, которые определяются из уравнения: λE − M −1 G2 = 0.
Здесь E – единичная матрица.
В исследованном нами случае КЭХ собственные значения имеют простой вид и равны: λ1,2 = −1.
Таким образом, исходная краевая задача Римана в матричном виде сведена к паре скалярных
задач Римана в следующем виде:
Ψ1 (ξ) = λ1 F1 (ξ) ,
Ψ2 (ξ) = λ2 F2 (ξ).
(26)
Решение краевых задач (26) строится в классе функций, допускающих интегрируемые особенности в узловых точках a, b, c(c′ ) и d(d′ ). Используется теория скалярной краевой задачи Римана
[17, 18]. В нашем случае, решения краевых скалярных задач Римана имеют следующий вид:
F (ξ) = C1 X(ζ),
112
F2 (ξ) = C2 X −1 (ζ),
X(ξ) =
(ζ − 1) (ζ + 1)
ζ
1/2
.
(27)
Ученые записки ЗабГГПУ
Напомним, что переменная ζ связана с исходной переменной z соотношением (12). Далее возвращаемся к исходной вектор функции. Определим фазы коэффициентов C1 , C2 . Согласно условиям
симметрии коэффициенты удовлетворяют соотношениям:
π
C1
=
exp
i
,
2
C¯1
π
C2
=
exp
−i
.
2
C¯2
(28)
Таким образом, получаем выражения для локальных токов:
π
π
j1 (ζ) = |C1 | exp i
X (ζ) = C1 X(ζ), j2 (ζ) = |C2 | exp −i
X −1 (ζ) = CX −1 (ζ).
4
4
(29)
Модули коэффициентов C1 , C2 определяются внешними токами (полями). А именно,
< j1 > + < j2 >= Jx − iJy = C1 < X(ζ) > +C2 < X −1 (ζ) >
(2)
−1
< e1 > + < e2 >= Ex − iEy = ρ(1)
(ζ) > .
xy C1 < X(ζ) > +ρxy C2 < X
(30)
Введем тензор эффективной холловской проводимости:
e
0
σxy
Jx
Ex
=
.
e
Jy
σyx
0
Ey
Таким образом, согласно формулам (28)–(30), компоненты эффективного тензора холловской
проводимости равны:
e
σxy
=
IRe(C1 + C2 )
(1)
Im(ρxy C1
+
(2)
ρxy C2 )
e
, σyx
=
IIm(C1 + C2 )
(1)
(2)
Re(ρxy C1 + ρxy C2 )
.
(31)
Здесь I =< X(ζ) >=< X −1 (ζ) > .
Во-первых, как можно видеть из формул (31), компоненты тензора холловской проводимости
e
e
определяются через константы |C1 |, |C2 | и не равны друг другу в общем случае: σxy
6= σyx
.
Во-вторых, значения константы |C1 |, |C2 | определяются условиями на бесконечности. В качестве примера рассмотрим два вида условий на бесконечности.
1) Пусть задано среднее значение тока на бесконечности.
J1 = const, Jy = 0
В этом случае из формулы получим:
(|C1 | + |C2 |) 2cos (π/4) I = Jx , (|C1 | − |C2 |) 2sin (π/4) I = Jy = 0.
Поэтому следующее соотношение получается: |C1 | = |C2 |. В случае равенства этих констант
среднему геометрическому значению
q
(1) (2)
|C1 | = |C2 | = σxy σxy .
(32)
получим выражение для компоненты тензора:
(1) (2)
⊥
σxy
=< ρ−1
xy >=
2σxy σxy
(1)
(2)
σxy + σxy
.
(33)
Это значение соответствует решению с постоянным электрическим током.
2) Рассмотрим второе условие, когда задано электрическое поле на бесконечности. В этом случае получим:
(2)
(1)
(2)
ρ(1)
|C
|
+
ρ
|C
|
2sin
(π/4)
I
=
E
,
ρ
|C
|
−
ρ
|C
|
2cos (π/4) I = Ey = 0.
1
2
x
1
2
xy
xy
xy
xy
Поэтому следующее соотношение получается:
(2)
ρ(1)
xy |C1 | = ρxy |C2 |
113
Физика, математика, техника, технология
В случае равенства этих констант следующим значениям:
|C1 | =
1
(1)
ρxy
, |C2 | =
1
(34)
(2)
ρxy
получим выражение для компоненты тензора:
⊥
σxy
=< σxy >=
1 (1)
(2)
σxy + σxy
.
2
(35)
Это значение соответствует решению с постоянным электрическим полем.
5. Заключение
В настоящей работе исследован вопрос о формировании ненулевого значения эффективной
холловской проводимости в условиях квантового эффекта Холла. Показано, что оно обусловлено протеканием холловского тока через конечное число сингулярных точек. Согласно результатам
вычислений, задание условий на бесконечности определяет характер протекания тока через сингулярные точки – в нашей модели это углы на стыке фаз. Существует две возможности протекания
краевых холловских токов. Первая их них – протекание краевого тока последовательно из первой во
вторую фазу через угловые контакты. Этот результат соответствует решению с постоянным током
на бесконечности (33). В такой же системе возможно, однако и другое решение, когда задано среднее поле на бесконечности (35). В этом случае краевые холловские токи текут отдельно – каждый
по своей фазе с проникновением через углы полуполос, расположенных в шахматном порядке.
Полученные результаты могут быть использованы для описания постоянной Холла в режиме
сильных магнитных полей.
Приложение
Соответствие точек при конформом отображении
Найдем координаты отображение точек в комплексной плоскости ξ = ξ + iη.
π
ez L
−z π
∞
−∞
L
+e
πz
=
= e +e
= ∞;
L
2
2
При z = α = (∞; 0), z = ∞, ζα = cosh
π
При z = b = (0; 0), z = 0, ζb =
π
z
−z
L
e L +e
2
При z = c = (0; L), z = Li, ζc =
=
e0 +e−0
2
Li π
−Li π
L
e L +e
2
При z = d = (∞; L), z = ∞ + Li, ζd =
При z = c′ = (0; −L), z = −Li, ζc′ =
=
e
zπ
−z π
L
e L +e
2
−Li π
Li π
L +e
L
e
2
При z = d′ = (∞; −L), z = ∞ − Li, ζd′ =
= 1;
iπ
+e−iπ
2
=
=
zπ
−z π
L
e L +e
2
= cosπ = −1;
∞ iπ
e
e
−iπ
+e−∞ eiπ
2
iπ
e
+e
2
=
∞ −iπ
e
e
= e∞ cosπ = −∞;
= cosπ = −1;
+e∞ eiπ
2
= e∞ cosπ = −∞;
Список литературы
1. Архинчеев В. Е. О неподвижных точках, инвариантах преобразований Дыхне и устойчивости решений задач эффективной проводимости // Письма в ЖЭТФ, 1998. Т. 67. Н7–8.
С. 951–958.
2. Архинчеев В . Е ., Батыев Э . Г . // Sol. St. Commun, 1989. Vol. 12. P. 1059–1060.
3. Балагуров Б. Я. Соотношепия взаимности в двумерной теории протекания // Журнал
экспериментальной и теоретической физики, 1981. Т. 81. Выи. 2(8). С. 665–671.
4. Дыхне А. М. О вычислении кинетических коэффициентов сред со случайными неоднородностями // Журнал экспериментальной и теоретической физики, 1967. Т.52. № 1. С. 264–266.
5. Дыхне А. М. Проводимость двумерной двухфазной системы // Журнал экспериментальной и теоретической физики, 1970. Т. 59. № 1. С. 110–115.
6. Емец Ю. П. Преобразование симметрии двухмерной двухкомпонентной электропроводной системы // ЖЭТФ, 1989. Т. 96. С. 701–711.
7. Емец Ю. П. Электрические характеристики композиционных материалов с регулярной
структурой. Киев: Наукова думка, 1986. 191 с.
8. Лаврентьев M. A., Шабат Б. В. Методы теории функций комплексного переменного. М.,
1973. 716 с.
114
Ученые записки ЗабГГПУ
9. Маркушевич А. И. Краткий курс теории аналитических функций. Наука, М., 1978. 416 с.
10. Фокин А. Г. Макроскопическая проводимость случайно-неоднородных сред. Методы
расчета // Успехи физических наук, 1996. 166 б. В. 10. 1069–1093.
11. Cabo Montes de Osa A., D. Martinez-Pedrera, Phys, Rev. B, 2003. V. 67, 245310.
12. Dykhne A. M., Ruzin I. M., Phys. Review B, V. 50, P. 2369 (1994).
13. Keller J. B. , J.Math. Phys. 1964. V. 5, P. 548.
14. Kivelson S., Dung-Hai Lee, Shou-Cheng Zhang, Phys. Rev. B, 1992. V. 46, P. 2223.
Рукопись поступила в редакцию 28 апреля 2011 г.
115
Документ
Категория
Без категории
Просмотров
8
Размер файла
460 Кб
Теги
локальные, условия, поле, вычисления, слоистых, среда, квантовое, холла, распределение, эффекты, токов
1/--страниц
Пожаловаться на содержимое документа