CC BY
69
ФИЗИКА НИЗКОРАЗМЕРНЫХ СТРУКТУР
УДК 544.225 ББК 22.368
ТЕМПЕРАТУРНАЯ ЗАВИСИМОСТЬ УДЕЛЬНОЙ ПРОВОДИМОСТИ БИСЛОЯ ГРАФЕНА 1
Г.С. Иванченко, Ю.В. Невзорова
В данной работе проведено теоретическое исследование транспортных свойств двойного графенового слоя с различной структурой расположения слоев друг относительно друга. Выявлена зонная структура такой системы и закономерности изменения удельной проводимости в зависимости от различных конфигураций графеновых слоев. Исследование проводилось методом молекулярных орбиталей в я-электронном приближении. Обнаружен металлический характер проводимости, а также ее слабая зависимость от выбранного направления в кристалле. Выявлена слабая зависимость удельной проводимости от взаимного расположения слоев графена.
Ключевые слова: бислой графена, электрическое поле, проводимость, нанолента, модель Хаббарда.
Введение
Проблема создания наноструктур с заданными свойствами и контролируемыми размерами входит в число важнейших проблем XXI века. Ее решение революционизирует элек-© тронику, материаловедение, механику, химию, медицину и биологию. В последние годы интерес к наноструктурам и нанокластерным системам значительно возрос с точки зрения pj их применения в нано- и оптоэлектронике. В 2004 г. удалось экспериментально получить 2 отдельные графеновые слои толщиной в один атом углерода [7-10]. Следовательно, воз-« никла необходимость изучения различных свойств и характеристик данного материала. о Графен интересен не только с точки зрения возможных приложений, но и с фундаменталь-8 ной точки зрения - вследствие своих уникальных электронных свойств. Вблизи уровня Ферми . электроны в графене обладают линейной дисперсией, а энергетическая щель между вален-^ тной зоной и зоной проводимости отсутствует. Модель зонной структуры графена послужи-§ ла стартовой площадкой для изучения свойств графита, но в многослойных стопках, в час-g тности в бислое, взаимодействие между слоями существенно искажает свойства графена S [ibid.]. И только после 2004 г., когда отдельные слои графена были успешно изолированы, ~ стало ясно, что на примере этой системы фактически можно изучать новый вид частиц -
безмассовые заряженные квазичастицы, нигде больше не существующие в природе. Такие свойства этих частиц, как двумерность, спинорная природа, нулевая масса и отсутствие щели в спектре, приводят к ряду электронных явлений, не имеющих аналогов ни в каких других физических системах [8].
На основе бислоя графена возможно создание различных наноэлектромеханических систем, например нановесы, нанодиоды, нанотранзисторы. Одной из важнейших характеристик таких приборов является электропроводность [2; 6].
1. Зонная структура и удельная проводимость
В настоящей работе рассматриваются четыре высокосимметричных случая расположения графеновых слоев друг относительно друга: без смещения (рис. 1а); со смещением в половину длины связи (б), так называемый графит Бернала; со смещением в длину связи (в) и 3/2 длины связи (г). Для расчета зонной структуры бислоя графена необходимо знать интеграл перекрывания для волновых функций 2р -орбиталей.
а)
б)
в) г)
Рис. 1. Различные варианты расположения слоев друг относительно друга: а) смещения по оси У нет; б) смещение по оси У на половину длины связи; в) смещение на длину связи;
г) смещение на 3/2 длины связи
Общий вид волновой функции для 2р -орбитали в сферической системе координат:
т2 -= 4К 12 /ехр
¡-/21»
os©, здесь X - порядковый номер углерода; ав = 0.5292А - пер-1т_
в
вый Боровский радиус; / - радиус-вектор. Зная волновую функцию, получим значения
ав
интеграла перекрывания для атомов углерода, принадлежащих разным слоям графена (см. таблицу).
Таблица
Значения интеграла перекрывания, d - расстояние между слоями графена, А -смещение одного слоя относительно другого
О 2 л/Зво / /2 а0
0,8~с1 7.29* Ю"4 4.59*10"" 1.87*10^ 1.22*ю"
с1 3.55*1 ОТ3 2.41*10 5 1.14*10"5 7.94*106
1. 2~с] 1.44*1 о6 1.05*10й 5.56*10-' 4.06*10-'
Для нахождения дисперсионных соотношений пользуемся методом МО ЛКАО, то есть волновую функцию записывают в виде детерминанта, построенного из одноэлектронных волновых функций. В кристаллах с трансляционной симметрией функции, зависящие только от координат электрона, называют блоховскими.
Для рассматриваемой системы базисные блоховские функции имеют следующий вид:
Ул (г ) = X ехр{'(к/Л + к/2с2 Ж (г " Ас " 12С2 )
УВ (7)=-ГТ X ехр{'(к/1С1 + к4С2 )}Рв (г - 11С1 - 12С2 )
•Л 1^12
УС X ехр{'(к11с1 + к12С2 )}Рс (г - - 11С1 -
(А
(г)=_1= X еХр{'(к/Л + к12С2 )}PD (г - Л2 - 11С1 - 12С2)
(1)
Здесь рЛ, рв, рс, рО - атомные орбитали р-типа. В таком базисе, используя приближение ближайших соседей, определяем дисперсионные уравнения четвертого порядка для законов дисперсии л-электронов в четырех высокосимметричных случаях, которые описывают электронные переходы в бислое графена. Метод получения дисперсионных соотношений наглядно описан в [2].
Переходы вдоль слоев и между ними взаимосвязаны. Зонная структура бислоя графена наглядно изображена на рисунке 2. Для расчетов использовались следующие значения параметров: химический потенциал: |а = |с = 0 эВ; кулоновская энергия: и = 10 эВ.
а) б)
Рис. 2. Зонная структура двойного графенового слоя:
а) по вертикальной оси отложена энергия, по горизонтальным - проекции волнового вектора кх, ку;
б) вид сверху
104 Г. С. Иванченко, Ю.В. Невзорова. Температурная зависимость удельной проводимости бислоя графена
На рисунке 2 видны «особые» точки энергетической поверхности, в которых энергия принимает экстремальные значения и имеет линейный характер (конусы). В конусах в приближении слабых возбуждений кристалла существуют безмассовые дираковские квазичастицы [8]. Поэтому двухслойный графен является перспективным материалом для создания современных устройств электроники.
Выражение для тензора удельной проводимости, полученное методом функций Грина [1], двухслойного графена в двухзонной модели Хаббарда [4; 5] имеет следующий вид [3]:
2
ааР =
г ж е
кТ V д хк,
ЕЕк (я К (к) Ад + (д К (к) Лс
+
+ К (д)К1 (к ) + Ка (д К(к),
Ад + Ас 2
Еп,
ch
+
дс дс
Е_ к кТ^
+1
sh
кТ
У У
(2)
где величины ха р - компоненты групповых скоростей электронов в зоне Бриллюена;
А - функции Грина операторов рождения и уничтожения электронов;
п - количество переходов между слоями;
Е - энергия электронов;
гас - полученные дисперсионные соотношения для переходов между слоями. Температурная зависимость электропроводности, рассчитанная по формуле (2), представлена на рисунке 3.
1,61,41,2 • 1,00,80,6-СО 0,40,2 0,0
ЕЕ
1 О
(О
Е та
\
V.
~1— 50
—I— 100
150
т,к
200
250
—I— 300
Рис. 3. Температурная зависимость электропроводности бислоя графена
и
Заключение
В заключении можно сделать следующие выводы: - Удельная проводимость двухслойного графена не зависит от выбранного направления в кристалле и носит металлический характер. Так как перпендикулярная составляющая волнового вектора к± имеет непрерывный набор значений, то при любом выделенном направлении существуют точки, при которых проводящие и валентные зоны пересекаются (рис. 2а). Но следует отметить, что в случае ограничения пространства для движения электронов, тем-
пературная зависимость проводимости может измениться и для некоторых направлений принять вид, характерный для полупроводниковых материалов.
-Предсказана слабая зависимость проводимости от взаимного расположения слоев графена. Это связано с тем, что между слоями присутствует лишь слабое Ван-дер-ваальсовое взаимодействие, что приводит к малому значению интеграла перекрывания, а следовательно, и энергии перехода между слоями. Однако включение внешних полей может существенно увеличить параметры системы, и межслоевые электронные переходы будут намного больше влиять на различные физические свойства данной системы.
ПРИМЕЧАНИЕ
1 Работа проведена в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы [проект № НК-16 (3)], а также поддержана грантом № 61-2010-а/ВолГУ
СПИСОК ЛИТЕРАТУРЫ
1. Бонч-Бруевич, В. Л. Метод функций Грина в статистической механике / В. Л. Бонч-Бруевич, С. В. Тяб-ликов. - М. : Гос. изд-во физ.-мат. лит., 1961. - 312 с.
2. Дьячков, П. Н. Углеродные нанотрубки. Строение, свойства, применение / П. Н. Дьячков. - М. : БИНОМ. Лаборатория знаний, 2006. - 290 с.
3. Иванченко, Г. С. Проводимость двухслойных углеродных нанотрубок в рамках модели Хаббарда / Г. С. Иванченко, Н. Г. Лебедев // Физика твердого тела. - 2007. - T. 49. - Вып. 1. - С. 183-189.
4. Изюмов, Ю. А. Магнетизм коллективизированных электронов / Ю. А. Изюмов, М. И. Кацнельсон, Ю. Н. Скрябин. - М. : Физматлит, 1994. - 368 с.
5. Изюмов, Ю. А. Статистическая механика магнитоупорядоченных систем / Ю. А. Изюмов, Ю. Н. Скрябин. - М. : Физматлит, 1987. - 264 с.
6. Харрис, П. Углеродные нанотрубы и родственные структуры. Новые материалы XXI века / П. Хар-рис. - М. : Техносфера, 2003. - 336 с.
7. Novoselov, K. S. Electric field effect in atomically thin carbon films / K. S. Novoselov [et al.] // Science. -2004. - V 306. - P. 666-669.
8. Novoselov, K. S. Two-dimensional gas of massless Dirac fermions in graphene / K. S. Novoselov [et al.] // Nature. - 2005. - V. 438. - P. 197-200.
9. Stankovich, S. Graphene-based composite materials / S. Stankovich [et al.] // Nature. - 2006. - V. 442. -P. 282-286.
10. Zhang, Y. Experimental observation of the quantum Hall effect and Berry's phase in graphene / Y. Zhang, J. W. Tan, H. L. Stormer & P. Kim // Nature. - 2005. - V 438. - P. 201-204.
TEMPERATURE DEPENDENCE OF CONDUCTIVITY OF THE DOUBLE WALLED GRAPHENE
G.S. Ivanchenko, Yu.V Nevzorova
In the paper, a theoretical study of the transport properties of double walled graphene with the different structure of the layers relative to each other has been carried out. The band structure of the system and the variation of conductivity depending on the different configurations of graphene layers have been revealed. The study was carried out with using the molecular orbitals in the framework of the л-electron approximation. The metallic conductivity and also its weak dependence on the direction in the crystal have been found. A weak dependence of the conductivity on the mutual arrangement of the graphene layers has been revealed.
Key words: graphene bilayer, electric field, conductivity, nanoribbon, Hubbard model.
