УДК 538.9, 530.145
Магический радиус квантовой точки двухслойного АА графена Сбойчаков А.О. * ^ Рожков А.В.1, Рахманов А.Л.1
1 Федеральное государственное бюджетное учреждение науки Институт теоретической и прикладной
электродинамики Российской академии наук
Статья поступила в редакцию 01.09.2022 Одобрена после рецензирования 19.10.2022 Принята к публикации 31.10.2022
Аннотация
В статье исследуются электронные свойства квантовой точки, изготовленной из двухслойного AA графена круглой формы. Наблюдается дискретный набор радиусов точки, при которых электронные состояния с низкой энергией вырождены по квантовому числу, характеризующему четность волновой функции относительно перестановки слоев. По аналогии с "магическими углами" в подкрученном двухслойном графене мы называем эти радиусы "магическими". Такая особенность уникальна для структур AA и связана со специфической слоевой симметрией двухслойного AA графена: четность самого высокого занятого уровня энергии изменяется с симметричного по индексу слоя на антисимметричный, когда радиус точки AA равен ее магическому значению. Исследуется аналогия в электронной структуре между подкрученным двухслойным графеном при магическом угле поворота и квантовой точкой AA графена с магическим радиусом. Эта аналогия может быть полезна для теоретического описания электронных свойств подкрученного двухслойного графена.
Ключевые слова: двухслойный графен, электронная структура, квантовая точка
Введение
В последние годы мы стали свидетелями большого успеха в изготовлении двумерных (2Б) ван-дер-ваальсовых гетероструктур. Современный уровень технологий позволяет обеспечить исключительную степень контроля над геометрией и составом искусственно созданных 2Б-систем. Эти достижения открывают новые направления исследований и ставят оригинальные вопросы как перед теоретиками, так и перед экспериментаторами. Руководствуясь этими факторами, в данной статье мы исследуем интересный мезоскопический объект - квантовая точка, вырезанная из двухслойного АА графена (АА-БЬС). Двухслойный АА графен представляет собой один из двух типов высокосимметричной упаковки двухслойного графена (другой тип упаковки называется АБ или фаза Бернала) [1]. Хотя АА-БЬС менее изучен экспериментально, существуют работы [2] - [5], рассказывающие об изготовлении образцов АА-БЬС.
В работе исследуются электронные свойства квантовых точек круглой формы, изготовленных из двухслойного АА графена(АА^Б). Показывается, что существует дискретный набор радиусов точек, для которых электронные состояния точки являются вырожденными. По аналогии с "магическими углами" в теории подкрученного двухслойного графена (1БЬС) мы называем эти радиусы "магическими".
Вырождение при магическом радиусе происходит из-за исчезновения щели между самым низким пустым и самым высоким занятым одноэлектронными состояниями. Такое свойство является эксклюзивным только для АА^Б. Это связано со специфическим топологическим свойством киральных электронов в АА-БЬС: симметрия волновой функции самого высокого занятого уровня меняется с симметричного относительно перестановки слоев на антисимметричный, когда радиус АА^Б равен его магическому значению.
Идентификация последовательности магических радиусов является основным результатом этой работы. Вырождение основного состояния квантовой точки позволяет стабилизировать экзотические упорядоченные состояния, обусловленные взаимодействием. Кроме того, можно ожидать, что наше исследование будет полезно для дальнейшего теоретического описания электронных свойств 1БЬС с магическими значениями угла подкрутки. Мы утверждаем, что существует простая связь между магическими углами поворота 1БЬС и магическими радиусами АА^Б. Это соответствие основано
*Автор, ответственный за переписку: Артём Олегович Сбойчаков, [email protected]
6
©2022 ИТПЭ РАН
на том факте, что образец 1БЬС можно рассматривать как структуру, состоящую из чередующихся областей с АА и АБ упаковками, при этом электроны проводимости находятся в основном в областях с АА упаковкой [1].
Гамильтониан квантовой точки двухслойного АА графена
При исследовании электронных свойств АА^Б мы будем работать в низкоэнергетическом континуальном приближении. В однослойном графене электроны на низких энергиях описываются безмассовым гамильтонианом Дирака, где компоненты спинора фа соответствуют волновым функциям электрона в подрешетке а = а, Ь. Так, гамильтониан для электронов с импульсами в окрестности дираковской точки К (долина £=+1) имеет вид
^вьа = _-Л]р ^ (1)
где V - матрицы Паули, а Ур - скорость Ферми. В окрестности точки К' (долина £ = - 1) имеем
Нзь° = _ ьР а* V (2)
В двухслойном АА графене имеется два графеновых слоя, при этом атомы в разных слоях, но в одной подрешетке располагаются один над другим. Волновая функция электрона в каждой долине описывается четырехкомпонентным спинором вида
( Ф^" \
Фц1ъ Ф&а V ф£2Ь )
(3)
Предполагая, что перескоки электронов между слоями осуществляются только между ближайшими атомами углерода в разных слоях (амплитуда перескока Ь о), гамильтониан двухслойного АА графена в континуальном приближении можно записать в виде
hp
( 0 -дх -idy 1/г 0 ^
—дх +^ду 0 0 V/
Vi 0 0 — дх -£ду
\ 0 1/ — дх +^ду 0 J
(4)
В этом уравнении используются безразмерные единицы энергии и длины. Энергия измеряется в единицах t о, а длина в единицах I = vp ft о. Подставляя характерные значения амплитуд перескоков электронов внутри слоя и между слоями, t = 2.7 эВ, tо = 0.35 эВ, а также имея в виду, что vp = 3tао/2, получим: I « 12 ао, где ао - расстояние между ближайшими атомами углерода в графене.
Волновая функция квантовой точки двухслойного АА графена. Граничные условия и уровни энергии
Волновое уравнение Н^ Ф^ = ЕФ^ есть система четырех дифференциальных уравнений первого порядка. Для круговой квантовой точки удобнее всего работать в полярных координатах х = г cos в, у = г sin в . В общем виде решение есть
Фцы = [CpJn ((1 - Е) г) + BMJM ((1 + Е) г)] Ф51Ь = -гe^+W C Jm+I ((1 - Е)г) - ((1 + Е) г)]
Ф52а = [CpJn ((1 -Е)г) - ((1 + Е) г)] ' (5)
фе2Ь = -гe^+W C JM+I ((1 - Е) г) + B^+i ((1 + Е) г)]
где ц - орбитальный момент частицы, J^(z) - функция Бесселя, а Cм и Вц - комплексные коэффициенты, которые находятся из граничных условий и условия нормировки. В случае, когда Вц = 0, волновая функция симметрична относительно перестановки слоев, а если Cм = 0, то она антисимметрична.
Радиальная компонента тока в каждой долине, = ах cos в + £ау sin в, должна быть равна нулю на границе точки, г = R. Следуя работе [6], сформулируем граничные условия для волновой функции (5). Мы интересуемся низкоэнергетическим спектром. Поэтому для достаточно гладкой границы междолинным рассеянием можно пренебречь. При этих предположениях граничные условия можно записать в виде
М5 (в) — I
(R, в) =0
(6)
для любого угла в. Матрица М^ (в) равна
М5 (в)
mu
—i m\2¿
—mise'
\ mi^e
i^e -if
m-f)
i mi2e —mii í(í if
г mise'
if
mi4 ei(f -^в) \
—mi4e i mise
— mi4 mss —i ms4e
i(f+&)
—г mise'
if
г ms4e
(7)
— mss
/
где фаза f произвольна, а вещественные параметры maß связаны между собой соотношениями
ii
+ 'm^2 + mis + mi4 = 1
2 4 2 4
mss + ms4 + m2s + m24 = 1
ms (mn + mss) — m4 (mi2 + ms4) = 0 ms (mi2 — ms4) + mi4 (mn — mss) = 0
(8)
Матрица (7) с условиями (8) есть общий вид унитарной и эрмитовой матрицы 4х4, которая антиком-мутирует с радиальной компонентой оператора тока . Последнее условие гарантирует, что ток на границе точки обнуляется для каждого слоя графена.
Рассмотрим здесь простейший случай, когда недиагональные элементы матрицы М^ (в), перепутывающие слои, равны нулю, т.е. т\3 = т\4 = 0. При этом будем считать, что тц = mss = совф, а т\2 = тз4 = Бтф, где -я < ф < я - феноменологический параметр, определяемый детальной структурой границы точки. В принципе, ф может быть функцией в, но в статье такая возможность не рассматривается.
Уравнение (6) формирует систему четырех линейных уравнений для двух независимых кон-
М5 (в) — I
= 0, в
(9)
стант Вм и С^. Однако из условий, наложенных на матрицу М^, следует, что det
результате чего ранг системы редуцируется до двух. Это позволяет однозначно найти энергетический спектр квантовой точки. Оказывается, что уровни энергии удовлетворяют одному из двух следующих уравнений:
" ((1 -E)R) - ctg2 JM+i ((1 -E) R) = 0 ^ Bß = 0, _ Jß ((1 + E)R) + ctg2J+ ((1 + E) R) = 0 ^ СM = 0.
Если уровень энергии удовлетворяет верхнему уравнению, то волновая функция оказывается симметричной относительно перестановки слоев. Если же уровень энергии удовлетворяет второму из уравнений (9), то волновая функция, соответствующая этому уровню, антисимметрична относительно перестановки слоев. Таким образом, можно ввести квантовое число с = ±1, характеризующее четность волновой функции (если с = +1, то волновая функция симметрична относительно перестановки слоев, а если с = -1, то волновая функция антисимметрична). Помимо этого, существуют еще два квантовых числа: орбитальный момент ^ и главное квантовое число п(-ж < п < ж) . Уровень энергии, характеризующийся квантовыми числами п, ц и с есть
(R) = c+ *Jf, где zß;n - n-тый корень характеристического уравнения
(z) + ctg7¡J^+i (z) =0.
(10)
(11)
Будем считать, что все упорядочены в порядке возрастания главного квантового числа п, при этом z^■n есть первый положительный корень уравнения (11).
Магические радиусы
Рассмотрим зависимость низкоэнергетических уровней энергии как функций радиуса квантовой точки R. Для произвольного ф, вообще говоря, спектр не показывает электрон-дырочной симметрии. Электрон-дырочная симметрия, однако, существует при ф = 0, я и при ф = я/2. При ф = я имеем на границе ф^а в) = 0, а при ф = 0, ф^ъ в) = 0 . Если же ф = я/2, то на границе имеет место равенство 1ф^га $)| = 1Фе,гъ На рисунке 1 представлена зависимость первых нескольких близких
к нулю уровней энергии как функций радиуса точки R, рассчитанных при ф = я/2.
2
Поскольку в системе есть электрон-дырочная симметрия, то в недопированном образце все уровни энергии меньше нуля заняты электронами, а все уровни энергии выше нуля пусты (каждый уровень энергии четырехкратно вырожден, дважды по спину и дважды по долине). При малых радиусах точки заполненными уровнями энергии являются все уровни с с = + 1. При увеличении радиуса точки происходит пересечение уровня с уровнем Е—_;0. Это происходит при Я = Я_ ~ 1.4 (в единицах /).
Будем называть этот радиус первым магическим радиусом. При дальнейшем увеличении Я наступает момент, когда Е_+__ = Е—_0. Это происходит при Я = Я2 ~ 2.5 (второй магический радиус). При еще большем увеличении Я пересекаются некоторые другие уровни энергии и т.д. Таким образом, имеется бесконечный набор магических радиусов Щ. При прохождении Я через Щ симметрия волновой функции ближайшего к нулю снизу уровня энергии меняется с симметричного на антисимметричный. Отметим, что магические радиусы квантовой точки, сделанной из двухслойного АА графена, существуют, вообще говоря, при любом ф.
Рисунок 1 - Первые несколько близких к нулю уровней энергии Е£;п как функции радиуса квантовой точки Я,
рассчитанные при ф = п/2
Магические радиусы и магические углы
Таким образом, мы показали, что для квантовой точки, сделанной из двухслойного АА графена, существует бесконечный набор магических радиусов Щ, при которых основное состояние системы является вырожденным. Это обстоятельство чрезвычайно существенно, поскольку при наличии вырождения даже сколь угодно слабое взаимодействие электронов будет это вырождение снимать (по крайней мере частично), формируя некое нетривиальное основное состояние. Анализ основного состояния при учете электрон-электронного взаимодействия выходит за рамки данной статьи и будет произведен отдельно. Здесь же мы хотим указать на связь между магическими радиусами квантовой точки АА графена и магическими углами подкрученного двухслойного графена.
Чтобы установить связь между АА^Б и 1БЪС, рассмотрим следующие обстоятельства. Подкрученный двухслойный графен при малых углах подкрутки можно рассматривать как периодическую структуру областей с АА, АБ и ВА упаковками [1]. Известно, что низкоэнергетические электронные состояния локализованы в областях с АА упаковкой. Следовательно, можно предвидеть, что соответствующие волновые функции могут быть описаны формализмом, разработанным в данной статье. Гибридизация между ячейками муара расширяет локализованные в АА области одноэлектронные состояния на восемь когерентных зон (по одной зоне на долину, на спин и на четность слоев). Эти зоны отделены от других электронных состояний щелью порядка 0.1эВ. Это устанавливает желаемое отображение между спектром "магической" АА^Б и низкоэнергетической зонной структурой tBLG при магическом угле подкрутки. Когда система выходит из "магического состояния" (либо магического
радиуса, либо магического угла поворота), расположение зон резко меняется: вырождение исчезает, а ширина зоны увеличивается на порядок. Этот факт также имеет простое объяснение в рамках предлагаемой схемы. Действительно, можно ожидать, что ширина зоны будет наименьшей, когда вырождение зон будет максимальным. Это происходит, когда размер области AA, Raa, близок к магическому радиусу R\. Поскольку Raa пропорционально l/в, последнее условие может быть использовано для связи магического угла и магического радиуса области AA. Характерный размер ячейки муара можно оценить как L « аол/3/в, и для первого магического угла (в « l0) имеем L « l00ao. Что касается Raa, то он, естественно, меньше, чем L. Величину Raa можно оценить следующим образом. Известно [1], что при малых углах подкрутки сверхъячейку подкрученного двухслойного графена можно рассматривать как состоящую из областей с AA, AB и BA упаковками равной площади. Площадь ячейки сверхрешетки равна Ssc = \/3L2/2. Таким образом, имеем оценку: Raa ~ \/Ssc/3ii « 0.3L « 30ao. Согласно нашим результатам, в зависимости от граничных условий для первого магического радиуса мы имеем l7ao К R\ К 23ao, что согласуется с нашей оценкой для Raa подкрученного двухслойного графена при первом магическом угле.
Заключение
Итак, в статье изучены электронные свойства круглой квантовой точки, изготовленной из двухслойного АА графена. Наблюдается дискретный набор "магических" радиусов, при которых основное состояние вырожденопо четности волновой функции относительно перестановки слоев. Существует аналогия между "магическими углами" подкрученного двухслойного графена и "магическими радиусами" двухслойной графеновой квантовой точки типа AA. Существование "магических радиусов" уникально для структур AA и связано со специфической слоевой симметрией двухслойного AA графена. Аналогия между "магическими углами" и "магическими радиусами" может быть полезна для описания электронных свойств подкрученного двухслойного графена.
Список литературы
[1] A. V. Rozhkov, A. O. Sboychakov, A. L. Rakhmanov, and F. Nori, Electronic properties of graphene-based bilayer systems // Phys. Reports.-201б.- № 648.P.1.
[2] H.-V. Roy, C. Kallinger, and K. Sattler, Study of singleand multiple foldings of graphitic sheets // Surf. Sci. - 1998. -№ 407.P. 1.
[3] J.-K. Lee, S.-C. Lee, J.-P. Ahn, S.-C. Kim, J. I. B. Wilson,and P. John, The growth of AA graphite on (111) diamond // J. Chem. Phys. - 2008. -№129. P. 234Т09.
[4] Z. Liu, K. Suenaga, P. J. F. Harris, and S. Iijima, Openand Closed Edges of Graphene Layers // Phys. Rev. Lett. - 2009. - №102.P. 015501.
[5] J. Borysiuk, J. Soltys, and J. Piechota, Stacking sequencedependence of graphene layers on SiC (0001) -Experimental and theoretical investigation //J. Appl. Phys. - 2011. - № 109.P. 093523.
[6] A. R. Akhmerov and C. W. J. Beenakker, Boundary conditionsfor Dirac fermions on a terminated honeycomb lattice // Phys. Rev. B. - 2008. - № ТТ. P. 085423.
Magic radius of AA bilayer graphene quantum dot SboychakovA.O.*, RozhkovA.V., RakhmanovA.A.
Institute for Theoretical and Applied Electromagnetics of RAS, Moscow, Russia
CoBpeMeHHaa э.пектроlцннамнка, № 1 (1), 0KTa6pb-H0a6pb 2022
Abstract
The paper investigates the electronic properties of a quantum dot made of AA bilayer graphene of round shape. A discrete set of point radii is observed at which low-energy electronic states are degenerate with respect to the parity of the wave function with respect to the permutation of layers. By analogy with the "magic angles" in twisted bilayer graphene, we call these radii "magic". This feature is unique for AA structures and is associated with the specific layer symmetry of the AA bilayer graphene: the parity of the highest occupied energy level changes from symmetric in the layer index to antisymmetric when the radius of the AA dot is equal to its magic value. The analogy in the electronic structure between twisted bilayer graphene with a magic twist angle and the quantum dot of AA bilayer graphene with a magic radius is investigated. This analogy can be useful for the theoretical description of the electronic properties of twisted bilayer graphene.
Key words: Bilayergraphene, electronic structure, quantum dot