УДК 621.373.8
УПОРЯДОЧЕННЫЕ СОСТОЯНИЯ В МОДЕЛИ АА ГРАФЕНА С ГРУППОЙ СИММЕТРИЙ SU(4)
А.В. Рожков * ^ A.O. Сбойчаков1
1 Федеральное государственное бюджетное учреждение науки Институт теоретической и прикладной электродинамики Российской академии наук, Москва, Россия
Статья поступила в редакцию 15.05.2024 Одобрена после рецензирования 20.05.2024 Принята к публикации 21.05.2024
Аннотация
В работе мы обсуждаем классификацию упорядоченных состояний для модели двухслойного АА графена, основанную на матричной группе Яи(4). Теория со столь широкой группой симметрий возникает из-за объединения спинового и т.н. долинного вырождений в широкий мультиплет, типичный для графена и других двумерных материалов на его основе. Мы описываем построение эффективной модели, демонстрирующей Ви(4)-инвариантность, и намечаем этапы исследования полученного гамильтониана на основе метода среднего поля. Упорядоченные состояния такой системы описываются единым Яи(4)-инвариантным уравнением самосогласования. Мы явно сконструируем три различных класса решений выведенного уравнения самосогласования, а также обобщим имеющуюся в литературе топологическую схему классификации упорядоченных фаз на случай двухслойного АА графена.
Ключевые слова: двухслойный АА графен, упорядоченное состояние, группа Яи(4)
EDN DOKXRI
doi:10.24412/2949-0553-2024-210-18-22
1. Введение
Данная статья посвящена теоретическому исследованию электронной жидкости двухслойного АА-графена (АА-Г). Конкретный вопрос, который мы хотели бы рассмотреть, - это проблема классификации низкотемпературных несверхпроводящих упорядоченных состояний АА-Г. Как мы увидим, изучаемая система обладает значительным набором таких состояний, и в такой ситуации задача классификации стоит весьма остро. Кроме этого, исследование относительно простого АА-Г может помочь нам распространить предлагаемый подход на другие типы порядков (в частности, сверхпроводящие) и другие типы систем, таких как подкрученные двухслойки и двухслойный АБ графен. В данной статье мы изложим лишь основные технические идеи подхода. Дополнительные подробности можно почерпнуть в работе [1].
2. Описание системы и одночастичной модели
Мы исследуем возможные упорядоченные состояния двухслойного графена с АА упаковкой. Кристаллическая структура этого материала изображена на Рис. 1. Гамильтониан сильной связи для такой кристаллической структуры имеет следующий вид
Й° = - Е (^Аа <1Ва +Ь.С.) - ¿0 £ (^а, ¿п2а<х + Ь.С.) . (1)
(тп)гст паи
В данном выражении й П 1аа. и ¿п1аа — операторы рождения и уничтожения электрона со спином а в слое I = 1, 2 на подрешетке а = А, В в элементарной ячейке п. Запись (шп) обозначает, что узел подрешетки А в элементарной ячейке ш и узел подрешетки В в элементарной ячейке п образуют пару ближайших соседей. Амплитуда Ь = 2, 7 эВ (¿о = 0, 35 эВ) в ур. (1) описывает перескок на ближайшие атомы внутри (между) слоями.
*Автор, ответственный за переписку: Александр Владимирович Рожков, [email protected]
18
©2024 ИТПЭ РАН
Рисунок 1 — Слева — кристаллическая структура АА-Г. Видны два слоя графена, уложенные один на другой так, что над каждым атомом нижнего слоя находится один атом верхнего слоя. Элементарная ячейка состоит из четырех атомов, обозначенных А1, В1, А2 и В2. Амплитуда внутрислоевого перескока обозначена t, а t0 обозначает амплитуду межслоевого перескока. Справа - зона Бриллюэна (голубой шестиугольник) и элементарная ячейка Вигнера-Зейца (ромб) для АА-Г. Показаны элементарные вектора обратной решетки Ь^г. Символы К^г отмечают точки Дирака (углы зоны Бриллюэна). Долины вокруг каждой точки Дирака заштрихованы.
Известно, что данный тип двухслойного графена обладает так называемым идеальным нестингом: в каждой долине имеется один электронный лист поверхности Ферми (ПФ), а также почти точно совпадающий с ним дырочный лист ПФ. Чтобы доказать это, необходимо диагонализовать гамильтониан (1). В результате можно получить следующую форму гамильтониана:
Н
6 к 7к
7к
В этой формуле зонные энергии ек5) равны
^ кх) = - о — А /к |,
е к2) = - о + А /к|, е к3) =+¿0 -А /к|, е к4 =+¿0 + А /к|
_(3) к
(4)
(2)
(3)
где функция /к = 1 + 2ехр(3гкха0/2)еоэ (л/3куа0/2) задается через расстояние между ближайшими
атомами углерода а0. Опираясь на хорошо известные свойства функции /, несложно установить, что
(2 3)
зоны б к ' ) доходят до уровня Ферми ер = 0, при этом возникшая поверхность Ферми (ПФ) задается уравнением г>р^| = ¿0, где q = Кх;2 — к, а скорость Ферми равна г>р = 3а0£/2. Здесь важно отметить, что уравнение на ПФ (1) не зависит от номера зоны , а также (п) не зависит от номера точки Дирака. Обстоятельство (1) прямо указывает на нестинг ПФ в нашей модели АА-Г. Что касается (И), данное свойство является проявлением т.н. долинного вырождения, постоянно возникающего в графене и материалах на его основе. Формально долинное вырождение может быть выражено в виде формулы
ек8) « ( —1)>Р¿0),
(4)
верной для зон = 2, 3 вблизи обоих точек Дирака.
3. 8и(4)-инвариантная модель
Нестинг делает систему восприимчивой к диэлектрическим неустойчивостям разного вида. Из-за восьмикратного вырождения ПФ (двукратное вырождение по спину, по зоне и по точке Дирака) количество упорядоченных состояний может быть велико. В такой ситуации возникает вопрос о возможности классификации упорядоченных фаз. Мы предлагаем схему такой классификации на основе группы симметрий Яи(4). Оказывается, что при некоторых предположениях группа симметрий гамильтониана включает в себя матричную группу Яи(4), действующую в четырехмерном спин-долинном пространстве. Подобную модель мы сконструируем, упрощая гамильтониан (2) так, чтобы новый (эффективный) гамильтониан учитывал долинное вырождение как свою явную симметрию. С этой целью рассмотрим гамильтониан вида:
Н0 « £ (—1)>р|к| — ¿0) , (5)
в=2,3,|к|<К0 &
где мы ввели долинный индекс £, принимающий значения £ = К1 или £ = К2. Геометрическое построение, вводящее долины в обратном пространстве, представлено справа на Рис. 1. Таким образом, при £ = К1 оператор 7^рождает электрон в зоне в со спином а и полным импульсом К1 + к.
Аналогично при £ = К2 оператор рождает электрон в зоне в со спином а и полным импульсом
К 2 + к.
Для ясности и компактности обозначений удобно ввести мульти-индекс т, объединяющий индексы £ и а. Это позволяет упростить ур. (5), записав его таким образом:
Н,
о - Е
s=2,3,|k|<K0
(-!)■> F|k|- io)E
(6)
где суммирование по^иа символически заменено на суммирование по т.
Такая запись позволяет явно продемонстрировать Би(4)-инвариантность рассматриваемой модели. С этой целью возьмем произвольную матрицу и, принадлежащую группе Би(4): и € Би(4). Обозначим ее матричные элементы как итт' и введем следующее преобразование Боголюбова:
7ks
'7ks
(7)
Легко проверить, что такое преобразование оставляет гамильтониан (6) неизменным. Мы доказали, что группа Би(4) входит в группу симметрий эффективной одночастичной модели АА-Г.
Теперь необходимо изучить вопрос об инвариантности гамильтониана взаимодействия относительно преобразований (7). В самом общем виде взаимодействие можно выразить так:
Н
1
2N~C
Е
k 11' аа'
^kL' Pkl аР-kl 'а'
(8)
Вектор k - передаваемый импульс, функция V^aa, - фурье-образ потенциальной энергии (R), описывающей взаимодействие электрона в слое I на подрешетке а с другим электроном в слое I' на подрешетке а'. Оператор pkiа - это компонента Фурье одноузлового оператора плотности частиц. Анализ показывает, что оператор взаимодействия H;nt самого общего вида не является инвариантным относительно преобразований (7). Однако, если допустить, что взаимодействие VЦ'(R) является дальнодействующим, то можно показать, что неинвариантные члены малы. Инвариантная часть взаимодействия может быть записана как сумма трех слагаемых:
Н
eff
Я" eff | йeff | йeff dir + Hex + Hu >
(9)
где прямое взаимодействие —обменное взаимодействие —^х и ишЫарр -и" задаются следующими выражениями:
1
й eff Hdir
16 Nc
Е
qq'kmm'
v+(k)
1 + e
Я 0q-0q + k)
1 + ег( Фq' k) 7q+k2m + ^q 3m 7q+k3m)
X ^7q'2m'7q'-k2m' + 7q'3m'7q'-k3m'^ '
Н
eff
16 NC
E
qq' kmm'
V-(k)
1 _ g»(0q-0q+k) 1 _ gi( 0q'-0q'-k)
X ^7q2m7q+k3m7q'3m/7q'-k2m' + 7q3m7q+k2m7q'2m'7q'-k3m')
Н
eff u
16 NC
E
qq' kmm'
V-(k)
1 _ ei(<Pq-<Pq+k) 1 _ еК0q'-0q'-k)
X ^7q2m7q+k3m7q'2m/7q'-k3m' + lq3m7q+k2mTq'3m'7q'-k2m')
(10)
(11)
(12)
где У±(к) = ^2аа, (ук^а' ± УкОа) /4. Суммирование по парным мульти-индексам в уравнениях (10-12) гарантирует Яи(4)-инвариантность этих выражений. Доказательство этого утверждения полностью аналогично доказательству инвариантности ур. (6).
4. Метод среднего поля
Для исследования упорядоченных фаз в предложенной Яи(4)-инвариантной модели АА-Г мы будем опираться на приближение среднего поля при нулевой температуре Т = 0. Для этого нужно
m m
1
1
ввести параметр порядка. Как следствие широкой группы симметрий матричной природы, удобно считать, что параметр порядка Q - это матрица 4x4, преобразующаяся под действием Би(4). Параметр порядка связывается с «аномальными» средними следующим образом:
Q = ^ £(*+(<>> + У-(ЭР>). (13)
с р
В этой формуле - взаимодействие, усредненное по ПФ, а «аномальные» средние (Эр> являются математическим ожиданием (средним значением), вычисленным по отношению к упорядоченному основному состоянию, для операторной матрицы ©ч, заданной своими матричными элементами ©чтт/ = 7<^3т7<2т/. Это позволяет вывести следующую форму среднеполевого гамильтониана взаимодействия:
Я
МР
£ Тт^в« + в<э) . (14)
Учитывая эту формулу, а также выражения (6) и (13), мы выводим уравнение самосогласования на параметр порядка Q:
(У+ - l^2)Q^(QtQ) = F+Q -У^, (15)
в котором функция к задана интегралом:
fe(QtQ) = W ; =, (16)
2
так что ^(е) к еД^/З^2) - это плотность состояний для графена, вычисленная на одну спиновую проекцию и на одну долину. Интегрирование выполняется от нижнего до верхнего края зоны.
Уравнение самосогласования может быть несколько упрощено, если воспользоваться сингулярным разложением матрицы параметр порядка: Q = иОУ''' и Qt = УОи^ где матрицы и иУ унитарны, а О диагональна, с неотрицательными элементами ^ > 0, г = 1,..., 4. Тогда ур. (15) может быть переписано в следующем виде:
- (у2 - У2)к(О2) О = ЖО1У, (17)
где ^ = и^ € и(4).
5. Решения уравнения самосогласования
Наша следующая задача - решить полученное уравнение самосогласования. Каждое найденное решение соответствует низкотемпературной упорядоченной фазе АА-Г. К сожалению, нам не удалось доказать, что мы знаем все возможные решения уравнения самосогласования. Тем не менее, мы смогли явно сконструировать три типа решений, каждый такой тип обладает своим физическим смыслом.
Во-первых, можно искать эрмитовы решения ур. (15). Как показано в работах [1,2], эрмитовы решения могут быть сгруппированы в три топологически неэквивалентных класса. Первый класс описывает волну зарядовой плотности. Второй класс включает в себя волну спиновой плотности, волну долинной плотности и волну спин-долинной плотности. Упорядоченные фазы, образующие класс номер три, - это комбинации всех вышеупомянутых типов волн плотности.
Во-вторых, матрица Q может быть антиэрмитовой. В такой ситуации упорядоченные фазы характеризуются спонтанными межслоевыми кольцевыми токами. В зависимости от топологического класса решения эти токи могут переносить зарядовые, спиновые, долинные и спин-долинные кванты.
И в-третьих, матрица параметра порядка может быть ни эрмитовой, ни антиэрмитовой, но диагонализуемой, причем такой, что ее собственные значения либо чисто действительные, либо чисто мнимые. Упорядоченные состояния, соотносящиеся с такими Q, будут одновременно демонстрировать токовые и зарядовые упорядочения.
6. Заключение
Благодаря особенностям сотовой решетки, одноэлектронная дисперсия в графене и системах на его основе характеризуется дополнительным квантовым числом, долинным индексом £. Хотя во многих
о
отношениях индекс долины отличается от спинового квантового числа, для графена и материалов на его основе мы можем сформулировать теорию, которая включает эти два квантовых числа на равных. В нашей статье мы описали построение такой теории для АА-Г. Конечно же, Яи(4)-симметричная теория для материала на основе графена не может являться точной моделью, детально описывающей электронные свойства материала. Тем не менее это полезный теоретический инструмент. Исследование Яи(4)-симметричных моделей следует рассматривать как физически мотивированный подход, позволяющий разработать классификационную схему упорядоченных состояний в материалах на основе графена. Предложенный подход сводит воедино как уже хорошо известные «зарядово-упорядоченные», так и менее изученные «токо-упорядоченные» состояния, а также обнаруживает ранее, по-видимому, не обсуждавшиеся гибридные виды состояний со спонтанно нарушенной симметрией, объединяющие в себе как «зарядовые», так и «токовые» параметры порядка. Представляется перспективным поиск новых типов решений уравнения самосогласования (или доказательства, что таковых не существует), а также расширение формализма на сверхпроводящие параметры порядка и на другие материалы.
Список литературы
[1] Rozhkov A. V., Sboychakov A. O., Rakhmanov A. L. Ordering in the SU(4)-symmetric model of AA bilayer graphene // Phys. Rev. B. — 2023. — Nov. — Vol. 108. — P. 205153. — Access mode: https://link.aps.org/doi/10.1103/PhysRevB.108.205153.
[2] Nandkishore R., Levitov L. Quantum anomalous Hall state in bilayer graphene // Phys. Rev. B.— 2010. —Sep. —Vol. 82, no. 11. —P. 115124._
ORDERED STATES IN THE SU(4)-SYMMETRIC MODEL OF AA BILAYER GRAPHENE
Rozhkov A.V.1*, Sboychakov A.O.1
1 Institute for Theoretical and Applied Electromagnetics of RAS, Moscow, Russia
* [email protected] Abstract
In this paper we discuss the classification of ordered states for model of bilayer AA graphene based on the SU(4) matrix group. Such a wide group of symmetries arises due to the unification of spin and the so-called valley degeneracy into a broad multiplet, typical for graphene and other two-dimensional materials based on it. We describe the construction of an effective theory demonstrating the SU(4) invariance, and discuss the mean field investigation of the resulting Hamiltonian. Ordered states of such a system are represented by a unifying SU(4)-invariant self-consistency equation. We will explicitly construct three different classes of solutions for the derived self-consistency equation. Using the SU(4) invariance, we generalize a previously published topological classification scheme for non-superconducting ordered states in bilayer AA graphene.
Key words: bilayer AA graphene, ordered state, SU(4) group