Физика графена

Материал из Википедии — свободной энциклопедии
Перейти к навигации Перейти к поиску
Графен
См. также: Портал:Физика

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

Кристаллическая решётка[править | править код]

GrapheneLattice.png GrapheneBZ.png
Рис. 1. Слева. Изображение гексагональной решётки графена. Жёлтым цветом показана элементарная ячейка. На красной окружности расположены ближайшие к центральному атомы решётки с векторами: . Красные и синие круги соответствуют различным подрешёткам кристалла. — базис решётки. Этот базис нельзя выбрать однозначным образом, поэтому в литературе встречаются несколько вариантов[1][2]. Справа. Первая зона Бриллюэна выделена жёлтым с векторами обратной решётки Показаны точки высокой симметрии Γ, Μ, Κ, Κ'[1].

Кристаллическая решётка графена (см. рис. 1) состоит из правильных шестиугольников и её можно представить как аналог пчелиных сот, что эквивалентно двумерной гексагональной решётке с атомами углерода расположенными в узлах кристалла. В элементарной ячейке кристалла находятся два типа атомов, обозначенные A и B. Каждый из этих атомов при сдвиге на вектора трансляций (любой вектор вида , где m и n — любые целые числа) образует треугольную подрешётку из эквивалентных ему атомов, то есть свойства кристалла независимы от точек наблюдения, расположенных в эквивалентных узлах кристалла. На рисунке 1 представлены две подрешётки атомов, закрашенные разными цветами: синим и красным. Например, красный атом четырёхвалентного углерода связан ковалентно с тремя соседними синими атомами углерода, расположенными в плоскости, поэтому угол между связями составляет 120°, а четвёртый электрон делокализован по всему кристаллу. Такая конфигурация 2s и двух 2p атомных орбиталей называется sp²-гибридизацией. Четвёртый электрон занимает |2pz> состояние, — эта орбиталь ориентирована перпендикулярно плоскости графена. Именно эти электроны отвечают за уникальные электронные свойства графена и формируют π-зону.

Расстояние между ближайшими атомами углерода в шестиугольниках, обозначенное a, составляет 0,142 нм. Это расстояние занимает промежуточное место между двойной связью (длина С=С 0.135 нм) и одинарной (длина С-С 0.147 нм)[3]. Постоянную решётки (a0) можно получить из геометрических соображений: она равна , то есть 0,246 нм. Площадь элементарной ячейки составляет 0,051 нм² и концентрация атомов 3,9×1015 см−2. Если определить за начало координат точку, соответствующую узлу кристаллической решётки (подрешётка A), из которой начинаются векторы элементарных трансляций с длиной векторов, равной a0 и ввести двумерную декартову систему координат в плоскости графена с осью ординат, направленной вверх, и осью абсцисс, направленной по вдоль вектора , то тогда координаты базисных векторов, запишутся в виде[1]:

(1.1)

а соответствующие им векторы обратной решётки[4]:

(1.2)

В декартовых координатах положение ближайших к узлу подрешётки A (все атомы которой на рисунке 3 показаны красным) в начале координат атомов из подрешётки B (показаны соответственно синим цветом) задаётся в виде:

(1.3)

Для гексагональной решётки известно, что её обратная решётка тоже будет гексагональной. Динамика электронов в кристалле определяется первой зоной Бриллюэна, которая является шестиугольником. Можно выделить несколько точек высокой симметрии в этой зоне, а именно Γ — в центре зоны Бриллюэна, и несколько точек на краях зоны Μ — седловая точка или особенность ван Хова (англ. Van Hove singularity), Κ, Κ' — дираковские точки с координатами

Экспериментально кристаллическую структуру графена наблюдали с помощью просвечивающего электронного микроскопа. Наблюдения с атомарным разрешением показали высокое качество плёнок графена полученных методом механического расщепления[5]. Альтернативный способ визуализации с использованием сканирующей туннельной микроскопии позволил исследовать не только кристаллическую структуру, но и электронный спектр графена. Используя атомно-силовой микроскоп, возможно получить изображение графена в прямом пространстве, а в ультра-высоком вакууме дифракция медленных электронов позволяет получать информацию о качестве кристалла в обратном пространстве при росте графена во время термического разложения карбида кремния[6].

Зонная структура[править | править код]

Рис. 2. Зонная структура графена.

Зонная структура графена впервые рассчитана в статье[2] в приближении сильно связанных электронов. На внешней оболочке атома углерода находится 4 электрона, три из которых образуют связи с соседними атомами в решётке при перекрывании sp²-гибридизированных орбиталей, а оставшийся электрон находится в состоянии |2pz> (именно это состояние отвечает в графите за образование межплоскостных связей). В приближении сильно связанных электронов полная волновая функция всех электронов кристалла записывается в виде суммы волновых функций электронов из разных подрешёток с учётом только ближайших соседей.

где коэффициент λ — некий неизвестный (вариационный) параметр, который определяется из минимума энергии. Волновые функции электронов подрешёток имеют вид[7]:

где C отвечает за нормировку полной волновой функции,  — двумерный волновой вектор, вектор трансляций, пробегающий по всем элементарным ячейкам кристалла,  — вектора, направленные на два атома из подрешёток A и B в элементарной ячейке.

В приближении сильно связанных электронов интеграл перекрытия между соседними атомами (), то есть сила взаимодействия, быстро спадает на межатомных расстояниях и следующие за ближайшими атомы можно не учитывать. Другими словами — взаимодействие волновой функции центрального атома с волновыми функциями атомов, расположенных на красной окружности (см. рис. 3), вносит основной вклад в формирование зонной структуры графена, изображённой на рисунке 2.

где

Индексы c и v относятся к π*-зоне (зоне проводимости) и π-зоне (валентной зоне). Ноль энергии выбран в центре зоны для нелегированного графена. Уровень Ферми разделяет полностью заполненную электронами валентную зону с отрицательными энергиями от полностью свободной зоны проводимости с положительными энергиями при нуле температуры. Сама точка с энергией ноль называется дираковской точкой или точкой электронейтральности. Уровень Ферми пересекает особые точки зонной диаграммы и , где валентная и проводимости зоны соприкасаются. Это происходит из-за того, что количество 2pz электронов в кристалле равно половине доступных состояний с учётом вырождения по спину. Около этих точек зонная диаграмма графена приобретает вид конусов. Благодаря такому виду закона дисперсии квазичастицы в графене при низких энергиях подчиняются уравнению Дирака, а не уравнению Шрёдингера. Так как и расположены на краю зоны Бриллюэна, волновой вектор имеет амплитуду сравнимую с вектором обратной решётки. Несмотря на это в низкоэнергетическом приближении вблизи дираковских точек можно разложить полный волновой вектор на два, а именно , где малый вектор это отклонение полного волнового вектора от дираковской точки . Дираковские точки и образуют две независимые долины, движение в которых разнонаправленно. Наличие двух долин приводит к дополнительному двукратному вырождению спектра. Если пренебречь процессами перехода электронов между долинами, то рассеяние квазичастиц происходит только вблизи уровня Ферми и наличие второй долины просто добавляет коэффициент 2 к току, поэтому часто вторая долина игнорируется в расчётах. Следует заметить, что данное приближение теряет смысл при снятии долинного вырождения.

Зонная структура двумерного графена является первым шагом к расчёту зонной структуры трёхмерного кристалла графита[8]. Накладывая периодические граничные условия по выбранному направлению можно получить закон дисперсии для одномерных нанотрубок. Вводя дополнительно пентагоны вместо шестиугольников, получают дискретный спект нульмерных фуллеренов.

Фотоэлектронная спектроскопия с угловым разрешением — прямой способ измерить зонную структуру материала, что и было сделано для графена выращенного на карбиде кремния[9]. Было продемонстрировано согласие теоретических предсказаний, наличие линейного спектра и измеренных свойств материала.

Уравнение Дирака[править | править код]

Рис. 3. Изолинии постоянной энергии (формула (2.4)). Жирный чёрный шестиугольник — первая зона Бриллюэна. Показаны также красные окружности на краях первой зоны Бриллюэна, где закон дисперсии носителей линеен. и обозначают две долины в -пространстве с неэквивалентными волновыми векторами

Из уравнения (2.4) следует, что вблизи точек соприкосновения валентной зоны и зоны проводимости ( и ) закон дисперсии для носителей (электронов) в графене представляется в виде:

где  — скорость Ферми (экспериментальное значение[10] =106 м/с, то есть в 300 раз меньше скорости света в вакууме и формально электроны нерелятивистские по определению специальной теории относительности),  — модуль волнового вектора в двумерном пространстве с компонентами (κx, κy) отсчитанного от K или K' точек Дирака,  — редуцированная постоянная Планка. Здесь следует отметить, что такого рода спектром обладает фотон, поэтому говорят, что квазичастицы (электроны и дырки, энергия для последних выражается формулой ) в графене обладают нулевой эффективной массой. Скорость Ферми играет роль «эффективной» скорости света. Хотя Филип Уоллес (англ. P. R. Wallace) в 1947 году первым вывел закон дисперсии для графена[8] написали уравнение Дирака для носителей тока другие исследователи в 1984[11][12]. Здесь также следует обратить внимание на тот факт, что появление линейного закона дисперсии при рассмотрении гексагональной решётки не является уникальной особенностью для данного типа кристаллической структуры, а может появляться и при существенном искажении решётки вплоть до квадратной решётки[13][14]. Наличие одинаковых атомов в двух подрешётках графена делает конический спектр защищённым благодаря симметрии: возмущения, инвариантные относительно одновременного действия инверсии времени и пространственной инверсии, не могут привести к образованию щели в спектре, но если симметрия между подрешётками нарушена, то в этом случае отсутствие пространственной инверсии приведёт к возникновению запрещённой зоны[15].

Уравнение Дирака, линейный закон дисперсии и наличие двух долин следует напрямую из уравнения Шрёдингера для графена, зонной структуры при малых энергиях электронов. Новые квазичастицы возникающие при этом предельном переходе описываются двумерным уравнением Дирака для безмассовых частиц (состоит из четырёх дифференциальных уравнений первого порядка), причём спин электрона не учтённый в уравнении Шрёдингера не входит и в уравнение Дирака. Но это уравнение обладает аналогичной характеристикой, называемой псевдоспином, связанной физически с наличием двух подрешёток (рис. 1) в кристаллической структуре графена. В качестве античастиц в отличие от трёхмерного уравнения Дирака в графене появляются дырки, хотя в основном уравнении их не было. Область применимости этого приближения описывается условием .

Обычно спин электрона не принимают во внимание (когда отсутствуют сильные магнитные поля или пренебрегают спин-орбитальным взаимодействием), и гамильтониан уравнения Дирака записывается в виде[16]:

где  — вектор-строка, состоящий из матриц Паули. Этот гамильтониан описывает свободные квазичастицы в графене и чтобы добавить в него потенциал, нужно провести формальный переход от точного уравнения Шрёдингера с потенциалом к низкоэнергетичному приближению. Для слабых (по сравнению с ) и медленно меняющихся на расстоянии a такой переход сделать легко, а для дефектов кристаллической структуры таких как границы кристалла и точечные дефекты нужно исходить из точного уравнения для нахождения правильного вида уравнения Дирака. Применимость уравнения Дирака можно расширить если разложить точный гамильтониан кристалла не до первого порядка малости (соответствует уравнению Дирака), а до второго порядка по , что приведёт к значительному усложнению задачи, но позволит учесть треугольную деформацию конического закона дисперсии (рис. 3.), это приближение применяется для исследования слабой локализации в графене и оптике. Для кулоновского потенциала возникают определённые сложности связанные с расходимостью потенциала при малых расстояниях если примесь находится вблизи от решётки. Уравнение Дирака не применимо для исследования оптических свойств когда энергия квантов сравнима с .

Спиноры[править | править код]

Волновая функция для гамильтониана имеет вид столбца[16]:

где индексы соответствуют подрешёткам кристалла в прямом пространстве: A и B, а также долинам в обратном пространстве: и . Гамильтониан для долины можно записать кратко

Этот двумерный гамильтониан аналог уравнения Дирака для безмассовых частиц, за исключением скорости света, в качестве которой выступает скорость Ферми. Из трёхмерного уравнения Дирака следует существование ферми-частиц, то есть частиц обладающих полуцелым спином. В графене из формально подобного уравнения следует существование характеристики названной псевдоспин, которая имеет отношение только к распределению электронной плотности между подрешётками кристалла. Таким образом состояние псевдоспин вверх означает подрешётку A, а псевдоспин вниз — подрешётку B. Для двух долин в k-пространстве вводят характеристику изоспин, и электроны конечно имеют внутреннюю степень свободы: спин (не отражённую в этом гамильтониане для графена).

Решения для свободных частиц для долин и имеют различный вид для положительной энергии (электроны) и отрицательной энергии (дырки):

Здесь  — полярный угол волнового вектора.

Полный гамильтониан можно представить в более симметричной форме:

где единичная матрица τ0 действует на индексы долин. Тогда спинор имеет вид[16]:

Киральность[править | править код]

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

где для электронов киральность положительна, а для дырок — отрицательна. Матрицы Паули здесь не имеют отношения к спину электрона, а отражают вклад двух подрешёток в формирование двухкомпонентной волновой функции частицы. Матрицы Паули являются операторами псевдоспина по аналогии со спином электрона. Так как оператор киральности коммутирует с гамильтонианом, то киральность сохраняется, что в графене приводит к такому явлению, как парадокс Клейна. В квантовой механике с этим явлением связано нетривиальное поведение коэффициента прохождения релятивистской частицей потенциальных барьеров, высота которых больше, чем удвоенная энергия покоя частицы. Частица более легко преодолевает более высокий барьер. В графене можно построить в задаче о преодолении потенциального барьера при нормальном падении не происходит отражения[17].

Фаза Берри[править | править код]

Дираковская точка[править | править код]

Плотность состояний и концентрация[править | править код]

Линейный закон дисперсии приводит к линейной зависимости плотности состояний от энергии, в отличие от обычных двумерных систем с параболическим законом дисперсии, где плотность состояний не зависит от энергии. Плотность состояний в графене задаётся стандартным способом

где выражение под интегралом и есть искомая плотность состояний (на единицу площади)[18]:

где и  — спиновое и долинное вырождение соответственно, а модуль энергии появляется, чтобы описать электроны и дырки одной формулой. Отсюда видно, что при нулевой энергии плотность состояний равна нулю, то есть отсутствуют носители (при нулевой температуре).

Концентрация электронов задаётся интегралом по энергии

где  — уровень Ферми. Если температура мала по сравнению с уровнем Ферми, то можно ограничиться случаем вырожденного электронного газа

Примечания[править | править код]

  1. 1 2 3 Katsnelson, 2012, с. 6.
  2. 1 2 Wallace P. R. The Band Theory of Graphite // Phys. Rev.. — 1947. — Т. 71. — С. 622—634. — doi:10.1103/PhysRev.71.622.
  3. Fuchs J., Goerbig M. O. Introduction to the Physical Properties of Graphene (англ.) : Lecture Notes. — 2008.
  4. Katsnelson, 2012, с. 7.
  5. Castro Neto et. al., 2009, с. 132.
  6. Andrei, 2012.
  7. Shung K. W. Dielectric function and plasmon structure of stage-1 intercalated graphite (англ.) // Phys. Rev. B. — 1986. — Vol. 34. — P. 979—993. — doi:10.1103/PhysRevB.34.979.
  8. 1 2 Wallace P. R. The Band Theory of Graphite (англ.) // Phys. Rev. 71, 622–634 (1947). — 1947. — Vol. 71. — P. 622—634. — doi:10.1103/PhysRev.71.622.
  9. Castro Neto et. al., 2009, с. 120.
  10. Novoselov et. al. Nature, 2005.
  11. DiVincenzo D. P., Mele E. J. Self-consistent effective-mass theory for intralayer screening in graphite intercalation compounds (англ.) // Phys. Rev. B. — 1984. — Vol. 29. — P. 1685—1694. — doi:10.1103/PhysRevB.29.1685.
  12. Semenoff G. W. Condensed-Matter Simulation of a Three-Dimensional Anomaly (англ.) // Phys. Rev. Lett. 53, 2449–2452 (1984). — 1984. — Vol. 53. — P. 2449—2452. — doi:10.1103/PhysRevLett.53.2449.
  13. Hatsugai Y., Fukui T., Aoki H. Topological analysis of the quantum Hall effect in graphene: Dirac-Fermi transition across van Hove singularities and edge versus bulk quantum numbers (англ.) // Phys. Rev. B. — 2006. — Vol. 74. — P. 205414. — doi:10.1103/PhysRevB.74.205414. — arXiv:cond-mat/0607669.
  14. Gusynin, 2007.
  15. Katsnelson, 2012, с. 13—14.
  16. 1 2 3 Katsnelson, 2012, с. 8—11.
  17. Katsnelson, 2012, с. 77—90.
  18. Ando T. Screening Effect and Impurity Scattering in Monolayer Graphene J. Phys. Soc. Jpn. 75, 074716 (2006) doi:10.1143/JPSJ.75.074716

Литература[править | править код]