Отклонение от линейного протекания тока в упорядоченной двумерной среде Текст научной статьи по специальности «Физика»

УДК 537.311.3

С.В. Хорьков

ОТКЛОНЕНИЕ ОТ ЛИНЕИНОГО ПРОТЕКАНИЯ ТОКА В УПОРЯДОЧЕННОЙ ДВУМЕРНОЙ среде

Нижегородский государственный технический университет им. Р.Е. Алексеева

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

Ключевые слова: неоднородные проводящие структуры, протекание, проводимость, нелинейная проводимость.

Проводящие среды интенсивно изучаются в связи с широким практическим применением. Ввиду особой роли дуальной симметрии [1] оказываются выделенными двумерные системы. Неоднократно отмечалась глубокая аналогия между двумерными случайными системами и периодическим решетками [1-4], позволяющие отнести их к одному классу универсальности в критической области, т.е. в близи перехода металл-диэлектрик. Двумерные системы привлекательны и с точки зрения экспериментальной реализации, т.к. методы современной нанотехнологии позволяют создавать плёнки с заданной текстурой.

В данной работе рассматривается плёнка с чередующимися в шахматном порядке ячейками. Пусть проводимости чёрной ячейки равна а1, а белой соответственно - о 2. Фрагмент такой плёнки представлен на рис. 1. Данная среда находится на пороге протекания, поскольку концентрации компонент материалов равны. Как показано в работе [1], эффективная линейная проводимость рассматриваемой среды зависит от проводимостей ячеек по закону oeff = л/о,о2 . Отметим, что такое же выражение для проводимости справед

ливо и для двумерной среды со случайным, но одинаковым по концентрации разбросом компонент [1]. Аналогично ведут себя и квадратичные корреляторы по-

Рис. 1

лей

2)

и токов [J

12

где скобки (...}}

обозначают усреднение по площади плёнки локального значения квадратов поля е и тока )2 соответственно. Это и позволяет отнести данные среды к одному классу универсальности.

В ряде экспериментальных работ [5-8] был обнаружен нелинейный отклик неупорядоченных плёнок, находящихся в критической области в зависимости от концентрации хо-

© С.В. Хорьков, 2016.

рошо проводящей компоненты. В работе [7] было показано, что нелинейная эффективная проводимость х^ неоднородной среды определяется коррелятором (с

Х* =(хе 4)/(ё4), (1)

где х - локальная нелинейная проводимость компонент, а с -локальное электрическое поле эквивалентной среды, полученное в линейном приближении. Выражение (1) справедливо только в случае слабой нелинейности, когда в локальном законе Ома добавляется следующее

слагаемое в разложении тока ] по полю с , т.е.

] = ъе + хе3. (2)

Эффективная линейная и нелинейная проводимости многокомпонентной среды входят в выражение для усреднённого тока J :

J = ^ Е + XeJfE3, (3)

где J = ^ Е = е} . Нелинейные отклик среды будет существенным, если в выражении (3) слагаемые станут одного порядка. Определим критическое поле нелинейности выражением

Ec=/% • (4)

а критический ток

Jc = G effEc . (5)

Задачей данной работы является определение нелинейной эффективной проводимости %eff среды по выражению (1), фрагмент которой, представлен на рис.1. Как отмечалось выше

для этого необходимо вычислить распределение электрического поля e в среде в линейном приближении. Такая задача была решена в работе [3]. Точное решение было получено с привлечением методов теории краевых задач аналитических функций. В силу двумерности среды удобно воспользоваться комплексным представлением векторов поля e(z) = ex - iey и тока j(z) = jx - ijy (здесь z = x + iy), которые являются аналитическими функциями внутри

квадратов (исключая границы). Для комплексных векторов закон Ома сохраняет прежний вид: j(z ) = o(z )e(z). Приведем выражение для электрического поля в ячейках с проводимостя-ми Gj и g2 соответственно

e1 (z) = C(cos 0 + sin 0)exp [in/ 2 (l/2 - y)]X (z) + + C(cos0 - sin 0)exp [in/2 (y - l/2)]X- (z), (6)

С

e2( z) = с (cos 0 + sin 0)exp [in/ 2(l/2 + y)]X (z) + Vh

C

+ C (cos 0 - sin 0)exp [in/ 2 (-y - 1/2)]X - (z), (7)

Vh

где

-,2y

X (z) =

cn

(Kz/L, k )

(8)

ап {Кг/Ь, к ]ёп(Кг!Ь, к )

5п(.),сп(.)^п(.) - эллиптические функции Якоби; К - полный эллиптический интеграл с модулем к (для квадрата к = 1/л/2, К « 1.8541 ), Ь - длина стороны квадратной ячейки; 0 - угол

.3/2

; I = -

между средним полем Е и осью " ох"; С = Ел/1 + Н21;

' 2К cos(лу)

Г (.) - гамма-функция. Параметр у связан с отношением проводимостей ячеек Н = —^ урав-

а,

Г ,3+1) Г (3 -Г

4 2 ) I 4 2 ,

нением

(9)

Проанализируем распределение локального электрического поля в среде. Первое и второе слагаемые в (6) и (7) представляют собой частные решения, которые соответствуют направлениям внешнего поля по диагонали " ac" и "bd" соответственно (см. рис. 1). Функция X(г) вблизи вершин "a" и "c" имеет интегрируемые особенности

х (*) «4у-'

(10)

а в смежных углах 'Ъ" и "d":

X(*)к *2у. (11)

Понятно, что для функции X _1( *), соответствующей направлению внешнего поля по диагонали "bd" , особыми являются углы и 'Ъ". Таким образом, при произвольном направлении внешнего поля общее решение имеет особенности во всех углах структуры. На рис. 2 и рис. 3 иллюстрируется распределение локального электрического поля в смежных ячейках при направлении внешнего тока вдоль оси "ox". Особенности типа (10) в окрестности особого угла будут играть важную роль в формировании нелинейного режима протекания тока.

Рис. 2

Рис. 3

Качественное объяснение асимптотик, состоит в следующем. На расстояниях вблизи угла, много меньших чем размер ячейки, в рассматриваемой системе отсутствуют характерные пространственные масштабы. Поэтому решение уравнения Лапласа в этой области масштабов имеет степенной вид (10). Показатель степени зависит от распределения зарядов на границах областей с разной проводимостью и, тем самым, от параметра Н .

Приступим к непосредственному вычислению нелинейной эффективной проводимости рассматриваемой среды. Для простоты будем полагать, что средний ток J направлен по

диагонали «ас» квадратных ячеек (см. рис. 1). Как показано в работе [3], в этом случае локальные комплексные поля в соседних ячейках связаны соотношением

е2( г ) = еД г ), (12)

где черта сверху означает комплексное сопряжение. Запишем абсолютную величину электрического поля в квадрате с проводимостью в виде

^ 2 ) = IХ ( 2 ^ Е

г)| X(г)|Е , (13)

где

3/2

I (у) = *

Г 3+Х) Г Г 3

I 4 2, ( 4 2 ,

-1

, (14)

2К С08(лу )

Г ( л) - гамма-функция. Используя (12), (14) и (1) получим выражение для нелинейной проводимости, удобное для анализа:

х>+1)2 (|X(г)'). (15)

Как уже отмечалось, функция Х (г) имеет особенность и при некоторых значениях параметра у интеграл по г в (г)|4^ может расходиться. Можно показать, что величина

нелинейной проводимости увеличивается с уменьшением к, и при стремлении к к некоторому пороговому значению кс ведет себя как

С

х е, = щ- <16>

где

С =^2^ +х2){кс + 1) I = I{у ) (17)

С 41с4К24ЙС , с {Ус). (17)

Из выражения (9) нетрудно получить кс = (^2 -1)2 . При значениях к < кс интеграл для хег расходится. Поясним природу расходимости. Из общего выражения для х, нетрудно заметить, что аномальный рост нелинейной проводимости связан с ростом корреляторов /|е(г)|4\ , так как все остальные сомножители в формуле конечны. Основной вклад в

интеграл (|е( г)|у для х дают области вблизи особых углов с асимптотическим поведе-

1,2

4а,

—2 у

нием решения пропорционального г . Поэтому для оценки корpелятоpа ^и к ^ кс достаточно pассмотреть малые окрестности этих точек, а точнее окрестность одного особого угла в элементарной ячейке. При любом направлении внешнего поля решение в этой области имеет вид (10). На основе этих рассуждений получим

(|е(г)|4)12 г|Му <х |г"8у+^г. (18)

Интеграл (18) расходится при у> 1/4. Из соотношения (9) нетрудно получить соответствующий интервал по к : к < кс. Для определения степени расходимости по (к — кс), разложим величину у в показателе подынтегральной функции (18) в ряд по малому параметру (к — кс) и получим

(К .J r - * * .

С

При стремлении Н ^ Нс расходится лишь нелинейная проводимость, поэтому из формул (4), (5) для Ес и 3 с видно, что зависимость критического поля и критического тока нелинейности вблизи Нс будет полностью определяться поведением %, то есть

Ес, Зс к(Н - Нс У2. (19)

Выводы

Хорошо известно, что в однородных проводящих средах электрический ток пропорционален приложенному напряжению. Такая зависимость имеет место в слабых электрических полях. В сильных полях зависимость тока от напряжения может стать нелинейной. Сильным полем называется такое, в результате которого за время свободного пробега (время релаксации) носители тока приобретают скорость по порядку величины сравнимую со скоростью хаотического движения. Создание сильных полей в однородных проводниках практически невозможно из-за большого тепловыделения и вследствие этого их плавления.

В неоднородных проводящих средах нелинейные эффекты наблюдаются в сравнительно слабых полях [5]. Одно из возможных объяснений состоит в том, что в среде имеются области, в которых локальное электрическое поле может сильно превысить среднее поле в данном образце. Вследствие этого электронный газ сильно разогревается. Это приводит к зависимости локальной проводимости от температуры.

В данной работе изучается нелинейный отклик неоднородной двумерной среды со структурой шахматной доски. Получено выражение для эффективной нелинейной проводимости среды. Было показано, что ввиду особенностей локального тока и поля вблизи углов структуры линейная зависимость между средним током и средним полем невозможна при конечном значении отношения линейных проводимостей компонент ячеек. Найдены зависимости критического поля и тока нелинейности от параметра Н. Отметим, что в работе [ 10] изучалась нелинейная эффективная проводимость бинарной среды со случайным, но равным по концентрации распределением её компонент. Результаты, полученные там, существенно отличаются от результатов данной работы.

Библиографический список

1. Дыхне, A.M. Проводимость двумерной двухфазной системы // ЖЭТФ. - 1970. - T. 59. - C. 110-115.

2. Балагуров, Б.Я. О проводимости двумерных систем с макроскопическими неоднородностя-ми // ЖЭТФ. - 1980. - Т. 79. - С. 1561-1572.

3. Емец, Ю.П. Преобразование симметрии двухмерной двухкомпонентной электропроводной системы // ЖЭТФ. - 1989. - Т. 96. - С. 701-711.

4. Емец, Ю.П. Электрические характеристики композиционных материалов с регулярной структурой / Ю.П. Емец. - Киев: Наукова думка, 1986.

5. Gefen, Y. Nonlinear behavior near the percolation metal-insulator transition / Y. Gefen, W.H. Shih, R.B. Laibowitz, J.M. Viggiano // Phys. Rev. Lett. - 1986. - V. 57. - P. 3097-3100.

6. Chakrabarty, R.K. Nonlinear I — V characteristics near the percolation threshold / R.K.Chakrabarty, K.K. Bardhan, A. Basu // Phys. Rev. - 1991. - V. B44. - P. 6773-6779.

7. Stroud, D. Nonlinear susceptibilities of granular matter / D. Stroud, P.M. Hui // Phys.Rev-1988. - V. B37. - P. 8719-8724.

8. Aharony, A. Crossover from linear to nonlinear resistance near percolation // Phys. Rev. Lett. -1987. - V. 58. - P. 2726-2727.

9. Жигальский, Г.П. Шум вида 1/f и нелинейные эффекты в тонких металлических пленках // УФН. - 1997. - Т. 167. No.6. - C. 623-648.

10.Сатанин, A.M. Нелинейная проводимость неупорядоченной среды на пороге протекания A.M. Сатанин, С.В. Хорьков, А.Ю. Угольников // Письма в ЖЭТФ. - 1995. - Т. 62. - С. 301-304.

Дата поступления в редакцию 23.04.2016

S.V. Khorkov

THE DEVIATION FROM THE LINEAR CURRENT FLOW IN AN ORDERED TWO-DIMENSIONAL ENVIRONMENT

Nizhny Novgorod state technical university n.a. R.E. Alexeyev

We investigate the impact of specific areas of periodic two-dimensional environment on its nonlinear effective conductivity. It is shown that the distribution of local currents and fields substantially depends on the ratio of the linear conductivities of the components. This dependence leads to the divergence of the correlator to the fourth power of the field and, hence, to divergence of the response function associated with the correlator. Defined value of the ratio of the linear conductivities of the component at which the effective nonlinear conductivity of the studied environment detects a peculiarity in the behavior.

Kew words: inhomogeneous conducting structures, percolation, conductivity, nonlinear conductivity.