Chapitre 1

L'émergence de cônes de Dirac inclinés dans les ma-

tériaux à base de carbone

La recherche de nouveaux matériaux de Dirac au-delà du graphène est

en cours. En particulier, il y a eu récemment un regain d'intérêt pour une
classe de matériaux prédits pour héberger des cônes de Dirac inclinés dans la
structure de bande électronique. Cependant, malgré des eorts théoriques et
expérimentaux croissants, la recherche d'un matériau pratique bidimension-
nel à cônes de Dirac inclinés est toujours en cours. Dans ce chapitre, nous
discutons de l'émergence de cônes de Dirac inclinés dans deux nouveaux ma-
tériaux à base de carbone. Nous commençons par revoir la théorie de fond du
graphène, dévoilant l'existence de dispersions de bande linéaires isotropes ou
de cônes de Dirac au voisinage du niveau de Fermi dans la structure de bande
électronique. Nous discutons ensuite de l'émergence de cônes de Dirac incli-
nés et avec gap dans un réseau bipolaire à grille latérale en graphène. Dans
ce système, le gap de bande local et le paramètre d'inclinaison des cônes de
Dirac peuvent être ajustés après la fabrication du dispositif en contrôlant les
tensions de grille appliquées au super-réseau. Enn, nous démontrons l'exis-
tence de cônes de Dirac inclinés dans la structure de bande électronique d'un
réseau de nanotubes de carbone quasi-métalliques.

1.1 Théorie de fond : Les cônes de Dirac dans le graphène et

au-delà

Introduction
Dans cette section, nous calculons la structure de bande électronique du
graphène dans le formalisme de liaison forte. Nous trouvons que la structure
de bande n'a pas de gap électronique et est décrite par des dispersions de
bande coniques connues sous le nom de cônes de Dirac au voisinage du niveau
de Fermi (Section 1.1.2). Dans les matériaux bidimensionnels au-delà du
graphène, il est possible que les cônes de Dirac soient inclinés (Section 1.1.3)
ou aient un gap de bande (Section 1.1.4). Nous verrons que les cônes de Dirac
inclinés peuvent appartenir à une variété de classes ou de types (étiquetés I, II
et III), chacun possédant des dispersions de bande radicalement diérentes.

1
Cônes de Dirac dans la structure de bande électronique du gra-
phène
Structure de bande électronique du graphène
Le graphène est une monocouche bidimensionnelle d'atomes de carbone
arrangés dans une conguration en nid d'abeille (voir Fig. 1.(a)). La structure
en nid d'abeille du graphène est dénie par une maille

unitaire à deux √atomes
se répétant avec les vecteurs de réseau a1 = a 1, 0 et a2 = (a/2) 1, 3 où

a = 2.46 est la constante de réseau et les vecteurs sont dénis dans le plan
x − y . Les deux atomes de la maille unitaire ne peuvent pas être mappés
l'un sur l'autre et constituent donc deux sous-réseaux distincts, désormais
étiquetés A et B . Les positions
√ relatives de ces deux√atomes
dans la maille
unitaire sont δ A = (a/2 3) 0, −1 et δ B = (a/2 3) 0, 1 de sorte que

chaque
√ atome de carbone a trois plus proches voisins à une distance ann =
a/ 3 = 1.42.
Chaque atome de carbone dans le graphène a quatre électrons de valence ;
trois de ces électrons occupent des orbitales hybrides sp2 dans le plan qui
forment des liaisons fortes dictant la structure en nid d'abeille tandis que
l'électron restant occupe une orbitale pz hors plan. Les propriétés électro-
niques du graphène peuvent être susamment comprises en considérant un
seul électron pz délocalisé dans un paysage de potentiel périodique résultant
du réseau xe de noyaux chargés positivement. Cette approximation est jus-
tement connue sous le nom d'approximation à un électron. Nous commençons
par l'équation de Schrödinger indépendante du temps pour un électron :

Ĥ(r) |Ψ(r)⟩ = E |Ψ(r)⟩ , (1)

(1.1)
où Ĥ(r) est l'hamiltonien, |Ψ(r)⟩ est la fonction d'onde, E est l'énergie et
r est la coordonnée de position de l'électron. L'hamiltonien est déni comme

ℏ2 2
Ĥ(r) = − ∇ + V (r), (2)
2me
(1.2)
où me est la masse de l'électron et V (r) est le potentiel cristallin qui
a la périodicité du réseau en nid d'abeille V (r) = V (r + na1 + ma2 ) où n
et m sont des entiers. Suivant la procédure standard de liaison forte, nous
pouvons écrire la fonction d'onde comme une combinaison linéaire d'orbitales
atomiques [32], c'est-à-dire des sommes de Bloch de sous-réseau
1 X ik·(na1 +ma2 )
|ΦA (r, k)⟩ = √ e φ(r − na1 − ma2 − δ A ), (3)
N n,m

(1.3)
et

2
 1 X ik·(na1 +ma2 )
|ΦB (r, k)⟩ = √ e φ(r − na1 − ma2 − δ B ), (4)
N n,m
(1.4)
où N est le nombre de sites atomiques dans chaque sous-réseau, k est
le vecteur d'onde cristallin et φ(r − R) est une orbitale atomique située sur
le site de réseau R. En eet, la fonction d'onde |Ψ(r, k)⟩ sera répartie sur
les deux sous-réseaux et est donc écrite comme une combinaison linéaire des
deux sommes de Bloch de sous-réseau
 
cA (k)
|Ψ(r, k)⟩ = = cA (k) |ΦA (r, k)⟩ + cB (k) |ΦB (r, k)⟩ , (5)
cB (k)
(1.5)
où nous avons fourni deux notations équivalentes utilisant soit un vecteur
colonne, soit une notation de ket. On Rsuppose que les orbitales atomiques
sont hautement localisées de sorte que φ∗ (r − R)φ(r − R′ )d2 r = δR,R′ où
δi,j est le delta de Kronecker - c'est ce qu'on appelle le modèle de liaison forte
semi-empirique. Nous notons que les sommes de Bloch dans les Eqs.(1.3) o

Figure 1  (a) Le réseau en nid d'abeille du graphène est esquissé. La

maille unitaire contient deux atomes de carbone A (bleu) et B (gris) qui
sont situés aux positions δ A et δ B dans la maille unitaire de sorte que chaque
atome de carbone a trois plus proches voisins équidistants. La maille unitaire
du réseau est répétée à travers le réseau avec les vecteurs de réseau a1 et
a2 avec une constante de réseau a = |a1 | = |a2 | = 2.46. (b) La zone de
Brillouin de l'espace réciproque est esquissée en fonction du vecteur d'onde
k = (kx , ky ) normalisé par la constante de réseau a avec les vecteurs de réseau
réciproque b1 et b2 et les points de haute symétrie Γ, K et K′ . (c) Structure
de bande électronique (E ) du graphène tracée en unités d'électronvolts (eV)
- les états électroniques sont occupés (bleu) jusqu'au niveau de Fermi EF =
0 et inoccupés au-dessus (gris). La structure de bande près du niveau de
Fermi est décrite par des cônes de Dirac non équivalents aux points de haute
symétrie K et K′ .

et (1.4) démontrent que le vecteur d'onde n'est pas unique√- translater

√ le

vecteur d'onde par un vecteur de réseau réciproque b1 = 2π/ 3a 3, −1
√

ou b2 = 4π/ 3a (0, 1) n'aecte pas la fonction d'onde car ils satisfont la
relation suivante avec les vecteurs de réseau : ai · bj = 2πδi,j . Ces vecteurs de
réseau réciproque dénissent la zone de Brillouin du cristal (voir Fig. 1.1(b)).
Utilisant laP
relation de complétude, nous pouvons dénir la matrice iden-
tité 2 × 2 I = α=A,B |Φα (r, k)⟩ ⟨Φα (r, k)| qui peut être utilisée pour déve-
lopper l'hamiltonien Ĥ(r) dans la base des sommes de Bloch de sous-réseau,
donnant

3
 X
Ĥ(r) = |Φα (r, k)⟩ Hαβ (k) ⟨Φβ (r, k)| , (6)
α,β

(1.6)
avec

Hαβ (k) = ⟨Φα (r, k)| Ĥ(r) |Φβ (r, k)⟩ , (7)
(1.7)
où α, β = A, B . Combiner les Eqs. (1.5) et (1.6) nous permet de for-
muler l'Eq. (1.1) comme le problème aux valeurs propres H(k) |Ψ(k)⟩ =
E(k) |Ψ(k)⟩ où H(k) est l'hamiltonien de Bloch 2×2 et E(k) sont les valeurs
propres d'énergie. Écrire cette expression en utilisant la notation d'algèbre
linéaire donne
    
HAA (k) HAB (k) cA (k) cA (k)
= E(k) . (8)
HBA (k) HBB (k) cB (k) cB (k)
(1.8)
Pour calculer la structure de bande électronique du graphène, nous de-
vons déterminer les éléments de l'hamiltonien de Bloch. Combinant les Eqs.
(1.3), (1.4) et (1.7), nous pouvons écrire l'élément non diagonal de l'hamil-
tonien de Bloch

1 X X ik·[(n−n′ )a1 +(m−m′ )a2 ]

HAB (k) = e ×
N ′ ′ n,m
n ,m
Z
φ∗ (r − n′ a1 − m′ a2 − δ A )Ĥ(r)φ(r − na1 − ma2 − δ B )d2 r.
(9)

(1.9)
Utilisant l'intégrale de saut de liaison forte entre sites plus proches voisins

Z
φ∗ (r − R)Ĥ(r)φ(r − R′ )d2 r = t, quand |R − R′ | = ann , (10)

(1.10)
où t ≈ −3.0eV donne la solution HAB (k) = f (k) où
h √
i
f (k) = t 1 + 2e− 3iky a/2 cos kx2a . (11)

(1.11)
Répéter le même calcul pour l'autre élément non diagonal donne HBA (k) =
f ∗ (k). Pour déterminer les éléments diagonaux de l'hamiltonien de Bloch, la
même procédure

4
===== Page 5 =====
est réalisée - cette fois en utilisant l'énergie sur site
Z
φ∗ (r − R)Ĥ(r)φ(r − R′ )d2 r = E0 , quand R = R′ . (12)

(1.12)
La structure de bande électronique E(k) peut maintenant être trouvée
en résolvant le problème aux valeurs propres suivant
    
0 f (k) cA (k) cA (k)
= E(k) , (13)
f ∗ (k) 0 cB (k) cB (k)
(1.13)
où nous avons posé E0 = 0. L'équation (1.13) a des solutions non triviales
lorsque l'énergie et l'hamiltonien de Bloch satisfont det[H(k) − E(k)I] = 0
ce qui donne la structure de bande électronique

E± (k) = ± |f (k)| , (14)

(1.14)
où la structure de bande est à demi-remplie, plaçant le niveau de Fermi à
EF = 0 de sorte que l'indice de bande ± correspond à la bande de conduction
(+) ou de valence (−). La structure de bande électronique du graphène a
été esquissée dans la Fig. 1.1(c). Les fonctions d'onde correspondant à ces
valeurs propres prennent la forme
" #
1 ± |ff (k)|
(k)
|Ψ± (k)⟩ = √ . (15)
2 1
(1.15)
Il est intéressant de noter que la structure de bande électronique du
graphène dérivée dans cette section a été calculée pour la première fois en
1947 par Philip Russell Wallace [33]. Écrit plus d'un demi-siècle avant la
découverte du graphène, cet article avait initialement pour but de produire
une compréhension des propriétés électroniques du graphite en mettant de
côté les interactions intercouches.

Émergence des cônes de Dirac dans la structure de bande basse

énergie
Le graphène est souvent appelé un semi-conducteur à gap nul ou un semi-
métal en raison de sa structure de bande sans gap avec une densité d'états qui
s'annule au niveau de Fermi. Les excitations électroniques basse énergie dans
le graphène sont déterminées par des dispersions de bande coniques connues

5
sous le nom de cônes de Dirac aux deux points non équivalents de haute√sy-
métrie dans la zone de Brillouin : K = (4π/3a)(1, 0) et K′ = (2π/3a)(1, 3)
(voir l'encart de la Fig. 1.1(c)). Les électrons dans le graphène peuvent exis-
ter dans l'une de ces deux vallées non équivalentes qui sont dégénérées en
énergie mais largement séparées en vecteur d'onde cristallin. Eectuer un dé-
veloppement de Taylor de l'hamiltonien au voisinage des points K (s = +)
ou K′ (s = −) donne l'équation de Dirac pour les électrons dans le graphène

Hs (q) = ℏvF sqx σx + qy σy , (16)

(1.16)

===== Page 6 =====

où s = ± est un indice de vallée, q est
√ l'écart du vecteur d'onde cristallin
par rapport au point K ou K , vF = − 3at/2ℏ ≈ 10 ms est la vitesse de
′ 6 −1

Fermi et
   
0 1 0 −i
σx = , σy = , (17)
1 0 i 0
(1.17)
sont les matrices de Pauli. Les dispersions des cônes de Dirac peuvent
être trouvées en diagonalisant l'hamiltonien, ce qui donne

E± (q) = ±ℏvF q, (18)

(1.18) q
où q = qx2 + qy2 est le vecteur d'onde radial. En eet, les cônes de Dirac
coniques ont une dispersion de bande linéaire et n'ont donc pas de masse
eective. Les excitations électroniques basse énergie dans le graphène sont
régies par l'équation de Dirac - la même équation qui décrit les particules
relativistes de spin 1/2. C'est pour cette raison que les électrons dans le
graphène sont dits se comporter comme des fermions relativistes sans masse.
La fonction d'onde au voisinage des cônes de Dirac prend la forme d'un
spineur à deux composantes

1 ±se−isφq
 
∗
Ψ± (q) = √ , (19)
2 1
(1.19)
où φq = arg(qy /qx ) est l'angle du vecteur d'onde. Les deux composantes
d'un spineur correspondraient typiquement au spin électronique, cependant,
dans le cas du graphène, ces deux composantes décrivent la fonction d'onde à
travers les deux sommes de Bloch de sous-réseau (voir Eq. (1.5)). En raison de
sa similarité mathématique apparente avec le spin en mécanique quantique,

6
ce degré de liberté découlant des deux sous-réseaux du graphène est appelé
pseudospin.

===== Page 5 =====

est réalisée - cette fois en utilisant l'énergie sur site
Z
φ∗ (r − R)Ĥ(r)φ(r − R′ )d2 r = E0 , quand R = R′ . (20)

(1.12)
La structure de bande électronique E(k) peut maintenant être trouvée
en résolvant le problème aux valeurs propres suivant
    
0 f (k) cA (k) cA (k)
= E(k) , (21)
f ∗ (k) 0 cB (k) cB (k)
(1.13)
où nous avons posé E0 = 0. L'équation (1.13) a des solutions non triviales
lorsque l'énergie et l'hamiltonien de Bloch satisfont det[H(k) − E(k)I] = 0
ce qui donne la structure de bande électronique

E± (k) = ± |f (k)| , (22)

(1.14)
où la structure de bande est à demi-remplie, plaçant le niveau de Fermi à
EF = 0 de sorte que l'indice de bande ± correspond à la bande de conduction
(+) ou de valence (−). La structure de bande électronique du graphène a
été esquissée dans la Fig. 1.1(c). Les fonctions d'onde correspondant à ces
valeurs propres prennent la forme
" #
1 ± |ff (k)|
(k)
|Ψ± (k)⟩ = √ . (23)
2 1
(1.15)
Il est intéressant de noter que la structure de bande électronique du
graphène dérivée dans cette section a été calculée pour la première fois en
1947 par Philip Russell Wallace [33]. Écrit plus d'un demi-siècle avant la
découverte du graphène, cet article avait initialement pour but de produire
une compréhension des propriétés électroniques du graphite en mettant de
côté les interactions intercouches.

Émergence des cônes de Dirac dans la structure de bande basse

énergie
Le graphène est souvent appelé un semi-conducteur à gap nul ou un semi-
métal en raison de sa structure de bande sans gap avec une densité d'états qui

7
s'annule au niveau de Fermi. Les excitations électroniques basse énergie dans
le graphène sont déterminées par des dispersions de bande coniques connues
sous le nom de cônes de Dirac aux deux points non équivalents de haute√sy-
métrie dans la zone de Brillouin : K = (4π/3a)(1, 0) et K′ = (2π/3a)(1, 3)
(voir l'encart de la Fig. 1.1(c)). Les électrons dans le graphène peuvent exis-
ter dans l'une de ces deux vallées non équivalentes qui sont dégénérées en
énergie mais largement séparées en vecteur d'onde cristallin. Eectuer un dé-
veloppement de Taylor de l'hamiltonien au voisinage des points K (s = +)
ou K′ (s = −) donne l'équation de Dirac pour les électrons dans le graphène

Hs (q) = ℏvF sqx σx + qy σy , (24)

(1.16)

===== Page 6 =====

où s = ± est un indice de vallée, q est
√ l'écart du vecteur d'onde cristallin
par rapport au point K ou K , vF = − 3at/2ℏ ≈ 10 ms est la vitesse de
′ 6 −1

Fermi et
   
0 1 0 −i
σx = , σy = , (25)
1 0 i 0
(1.17)
sont les matrices de Pauli. Les dispersions des cônes de Dirac peuvent
être trouvées en diagonalisant l'hamiltonien, ce qui donne

E± (q) = ±ℏvF q, (26)

(1.18) q
où q = qx2 + qy2 est le vecteur d'onde radial. En eet, les cônes de Dirac
coniques ont une dispersion de bande linéaire et n'ont donc pas de masse
eective. Les excitations électroniques basse énergie dans le graphène sont
régies par l'équation de Dirac - la même équation qui décrit les particules
relativistes de spin 1/2. C'est pour cette raison que les électrons dans le
graphène sont dits se comporter comme des fermions relativistes sans masse.
La fonction d'onde au voisinage des cônes de Dirac prend la forme d'un
spineur à deux composantes

1 ±se−isφq
 
∗
Ψ± (q) = √ , (27)
2 1
(1.19)
où φq = arg(qy /qx ) est l'angle du vecteur d'onde. Les deux composantes
d'un spineur correspondraient typiquement au spin électronique, cependant,
dans le cas du graphène, ces deux composantes décrivent la fonction d'onde à

8
travers les deux sommes de Bloch de sous-réseau (voir Eq. (1.5)). En raison de
sa similarité mathématique apparente avec le spin en mécanique quantique,
ce degré de liberté découlant des deux sous-réseaux du graphène est appelé
pseudospin.

Cônes de Dirac inclinés

La recherche de nouveaux matériaux de Dirac bidimensionnels au-delà du
graphène a conduit à la prédiction de cônes de Dirac avec des dispersions de
bande inclinées. L'exemple peut-être le plus simple de ceci est le graphène
sous contrainte quinoïde, c'est-à-dire le graphène déformé le long de l'axe
cristallographique ŷ. En contraignant le réseau de graphène, en calculant la
structure de bande électronique et en développant au voisinage des dégéné-
rescences de bande (voir Réf. [12]), on peut trouver que les quasiparticules
électroniques obéissent à l'hamiltonien de Dirac de type graphène avec un
terme d'inclinaison ajouté. Cet hamiltonien décrit des cônes de Dirac inclinés
dans la structure de bande électronique, prenant la forme :

Hs (q) = ℏv sγηqx I + sηqx σx + qy σy , (28)

(1.20)
où γ est un paramètre d'inclinaison, η est un paramètre d'anisotropie et
tous les symboles restants sont dénis dans les sections précédentes. Alors que
la dispersion de bande d'un cône de Dirac dans le graphène est déterminée
par une seule vitesse (de Fermi) vF , les cônes de Dirac inclinés ont trois

===== Page 7 =====

vitesses associées parfois écrites comme : vy = v (vitesse le long de l'axe
du vecteur d'onde qy ), vx = ηv (vitesse le long de l'axe du vecteur d'onde
qx ), et vt = γηv (vitesse d'inclinaison). En diagonalisant l'hamiltonien dans
l'équation (1.20), nous pouvons obtenir les énergies propres
 q 
s
E± (q) = ℏv sγηqx ± η 2 qx2 + qy2 , (29)

(1.21)
et les vecteurs propres
" sηqx −iqy #
1 ±√ 2 2 2
Ψs± (q) = √ η qx +qy . (30)
2 1
(1.22)
Notez que, contrairement au cas non incliné, la dispersion de bande élec-
tronique d'un cône de Dirac incliné dépend de l'indice de vallée s.

9
 Bien que le graphène sous contrainte quinoïde présente un modèle simple
pour comprendre les cônes de Dirac inclinés, en pratique, le paramètre d'in-
clinaison prend de petites valeurs de l'ordre du paramètre de déformation (où
la déformation réalisable expérimentalement dans le graphène est ε ∼ 0.1)
[12]. Cependant, dans les matériaux bidimensionnels au-delà du graphène, ce
paramètre d'inclinaison peut prendre une gamme de valeurs. Le paramètre
d'inclinaison dénit trois types de cônes de Dirac inclinés : les cônes de Dirac
de type-I sous-critiquement inclinés (|γ| < 1), de type-III critique (|γ| = 1) et
de type-II super-critiquement inclinés (|γ| > 1), avec la dispersion de bande
de chacun tracée dans la Fig. 1.2. Alors que les cônes de Dirac de type-I
ont des contours isoénergie elliptiques fermés, les cônes de Dirac de type-III
et de type-II ont respectivement des contours isoénergie paraboliques et hy-
perboliques ouverts. Le paramètre d'anisotropie (η ) a pour eet de mettre à
l'échelle la dispersion de bande le long de l'axe du vecteur d'onde qx .
Il convient de souligner à ce stade que dans cette thèse, nous nous concen-
trons sur les cônes de Dirac inclinés dans les matériaux bidimensionnels et ne
considérons pas les cônes de Weyl tridimensionnels inclinés [34] ni les cônes
de Dirac bidimensionnels à la surface de matériaux tridimensionnels [35, 36].
En eet, la liste des matériaux bidimensionnels prédits pour héberger des
cônes de Dirac inclinés dans la structure de bande électronique ne cesse de
s'allonger avec

Figure 2  Dispersion de bande E(q) de cônes de Dirac inclinés sans gap

avec diérents paramètres d'inclinaison : pas d'inclinaison (|γ| = 0), type-
I (|γ| < 1), type-III (|γ| = 1) et type-II (|γ| > 1). Dans tous les tracés,
le paramètre d'anisotropie est égal à l'unité (η = 1). Le niveau de Fermi
(EF = 0) est représenté par un plan translucide avec les états électroniques
occupés et inoccupés ombrés en bleu et gris respectivement.

===== Page 8 =====

des cônes de Dirac de type-I prédits dans le borophène 8-P mmn [9, 10,
11], des cônes de Dirac de type-II dans certains ditellurures monocouches [15]
et bien d'autres [8, 16, 17, 18, 19, 20, 21, 22, 12, 13, 14, 16]. En eet, tout
au long de ce chapitre, nous avons l'intention d'ajouter à cette liste avec des
cônes de Dirac inclinés émergeant dans la structure de bande électronique de
super-réseaux de graphène latéraux (Section 1.2) et des réseaux de nanotubes
de carbone (Section 1.3).

Cônes de Dirac avec gap

La nature sans gap des cônes de Dirac dans le graphène est protégée
contre les petites perturbations par la présence de la symétrie de renverse-

10
ment du temps et de la symétrie d'inversion [37]. Dans le modèle de liaison
forte précédemment discuté pour le graphène, la symétrie de renversement
du temps correspond à la valeur réelle des paramètres de saut tandis que la
symétrie d'inversion correspond à l'invariance du nid d'abeille de graphène
sous la transformation r → −r autour du centre de la maille unitaire. Si l'une
de ces symétries est brisée, un gap peut s'ouvrir dans les cônes de Dirac. Par
exemple, nous pouvons postuler un matériau cristallin avec une structure en
nid d'abeille où les deux sous-réseaux A et B contiennent des atomes dié-
rents. Dans ce cas, les énergies sur site sur chaque sous-réseau ne sont plus
équivalentes, ce qui donne les éléments diagonaux de l'hamiltonien de Bloch
HAA = E0 + ∆/2 et HBB = E0 − ∆/2, ouvrant ainsi un gap de bande de ∆
dans la structure de bande électronique.
L'un des exemples les plus connus de matériau à cône de Dirac avec gap
est peut-être la monocouche de silicium en nid d'abeille bombée connue sous
le nom de silicène [38, 39, 40]. En théorie, le silicène vierge aurait une symé-
trie d'inversion intacte, hébergeant ainsi des cônes de Dirac sans gap dans
la structure de bande électronique. Cependant, la nature bombée des deux
sous-réseaux permet la rupture de la symétrie d'inversion en présence d'un
champ électrique externe - ouvrant ainsi un petit gap de bande de l'ordre de
centaines de milli-électronvolts [41, 42]. Que ce soit en permettant des géo-
métries cristallographiques, l'application de champs externes, ou l'ingénierie
de super-réseaux (voir Section 1.2), nous pouvons

Figure 3  Dispersion de bande E(q) de cônes de Dirac inclinés avec gap

avec diérents paramètres d'inclinaison : pas d'inclinaison (|γ| = 0), type-
I (|γ| < 1), type-III (|γ| = 1) et type-II (|γ| > 1). Dans tous les tracés,
le paramètre d'anisotropie est égal à l'unité (η = 1). Le niveau de Fermi
(EF = 0) est représenté par un plan translucide avec les états électroniques
occupés et inoccupés ombrés en bleu et gris respectivement.

===== Page 9 =====

considérer des cônes de Dirac qui sont simultanément avec gap et incli-
nés. Les quasiparticules électroniques dans de tels matériaux obéiraient à
l'hamiltonien suivant :

∆
Hs (q) = (31)


σz + ℏv sγηqx I + sηqx σx + qy σy ,
2
(1.23)
où la troisième matrice de Pauli est
 
1 0
σz = . (32)
0 −1

11
(1.24)
La dispersion de bande d'un cône de Dirac avec gap et incliné prend la
forme
" s  2 #
∆
s
E± (q) = ℏv sγηqx ± η 2 qx2 + qy2 + . (33)
2ℏv
(1.25)
En l'absence d'inclinaison du cône de Dirac (γ = 0), le terme ∆ décrit
simplement l'amplitude du gap de bande direct du cône de Dirac (voir Fig.
1.3). Les cônes de Dirac de type-I avec gap (|γ| < 1) ont un gap de bande
indirect avec les minima de bande de conduction (+p ) et les maxima de bande
de valence (−) prenant les valeurs ±λ∆/2 où λ = 1 − γ 2 . Pour le cas des
cônes de Dirac de type-III (|γ| = 1) et de type-II (|γ| > 1), la séparation de
bande ne correspond pas à un gap de bande. La présence du terme de sépa-
ration de bande change la forme de la fonction d'onde spineur par rapport
au cas sans gap
 
1  ±2(sηqx − iqy )
Ψs± (q) = q , (34)
N± ∓∆ ˜ + 4η 2 q 2 + 4q 2 + ∆ ˜2
x y

(1.26)
où ∆e = ∆/ℏv est la séparation de bande réduite et le facteur de norma-
lisation prend la forme
r  q 2
N± = 4η 2 qx2 + 4qy2 + ∓ ∆ ˜ + 4η 2 q 2 + 4q 2 + ∆
e2 . (35)
x y

(1.27)

1.1.5 Conclusions
Dans cette section, nous avons développé une théorie de bande eective
pour les matériaux bidimensionnels avec une structure de bande électronique
hébergeant des cônes de Dirac inclinés avec ou sans gap. Tout au long du reste
de ce chapitre, nous étudions deux systèmes : un super-réseau de graphène
électrostatique (Section 1.2) et un réseau planaire de nanotubes de carbone
(Section 1.3) qui hébergent des cônes de Dirac inclinés avec et sans gap dans
la structure de bande électronique. La théorie développée dans cette section
sera le point de départ à partir duquel nous étudierons l'absorption (Cha-
pitre 2) et l'optovalleytronique (Chapitre 3) des matériaux à cônes de Dirac
inclinés et le comportement de ces matériaux sous des champs magnétiques
et électriques externes (Chapitre 4).

12
===== Page 10 =====
1.2 Ingénierie des cônes de Dirac dans les super-réseaux bi-
polaires de graphène

Introduction
Actuellement, il existe une pénurie de systèmes électroniques pratiques
et accordables qui présentent des cônes de Dirac inclinés bidimensionnels.
Plusieurs matériaux théoriques avec des géométries de réseau spéciques
ont été prédits pour supporter des cônes de Dirac électroniques inclinés, le
plus célèbre étant la monocouche de bore 8-P mmn borophène [9, 10, 11].
Cependant, après la synthèse, les structures cristallines ne peuvent pas être
pratiquement modiées pour ajuster l'inclinaison ou modier le gap de ces
cônes. Plutôt que de placer de vrais atomes dans une conguration de réseau
particulière, nous proposons d'aborder le problème en utilisant des atomes
articiels ; à savoir, des états liés piégés à l'intérieur de puits et de barrières
de graphène organisés en un super-réseau latéral.
Contrairement aux systèmes non relativistes, les puits et barrières électro-
statiques dans le graphène supportent tous deux des états liés. Ces états liés
sont localisés autour du centre des potentiels de connement, un peu comme
les orbitales atomiques dans un cristal sont centrées sur leurs positions de
réseau. Les états connés d'un puits et d'une barrière se chevauchent, un
peu comme les orbitales atomiques adjacentes. Ce chevauchement peut être
caractérisé par le paramètre de saut dans le célèbre modèle de liaison forte.
Contrairement à un cristal réel, où le chevauchement entre orbitales adja-
centes est xe, le chevauchement entre les fonctions de puits et de barrière
peut être complètement contrôlé. Cela peut être réalisé en faisant varier la
hauteur et la profondeur des potentiels de connement via leurs tensions de
grille supérieures. Ainsi, la construction d'un super-réseau à partir de puits
et de barrières dans le graphène imite la structure de bande d'un réseau
atomique mais avec l'avantage d'une nouvell

===== Page 11 =====

tunabilité. Ainsi, faire progresser l'ingénierie de la structure de bande en
physique de la matière condensée dans la même direction que le contrôle
optique dans les métamatériaux sur mesure [43].
Dans ce qui suit, nous montrons qu'un super-réseau latéral composé de
paires de puits et de barrières se répétant, c'est-à-dire un réseau bipolaire
(voir Fig. 1.4(a)), héberge des cônes de Dirac inclinés dans la structure de
bande électronique (Section 1.2.2). En contrôlant les tensions de grille ap-
pliquées aux barrières et aux puits, les cônes peuvent être déplacés entre
le type-I, -II et -III ; ainsi le paramètre d'inclinaison peut maintenant être
considéré comme un degré de liberté supplémentaire (Section 1.2.3). Nous

13
démontrons également que ces cônes de Dirac peuvent devenir avec gap si le
réseau bipolaire manque de symétrie de réexion (Section 1.2.4) (voir Fig.
1.4). Ceci contraste fortement avec les super-réseaux précédemment étudiés
dans le graphène qui se concentraient sur des systèmes avec symétrie de ré-
exion, qui ne possèdent pas de cônes de Dirac avec gap. Bien que nous
démontrions l'existence de cônes de Dirac inclinés et avec gap accordables
dans la structure de bande électronique, nous soulignons que ces cônes sont
des satellites d'un cône central sans gap (Section 1.2.5). Contrairement aux
cônes de Dirac dans le graphène vierge, ces cônes centraux sont anisotropes
dans l'espace des impulsions, possédant des contours isoénergie elliptiques.

Modèle de liaison forte du super-réseau de graphène

Bien que les structures réalistes à grille supérieure dans le graphène gé-
nèrent des potentiels guides lisses [44, 45, 46, 47] (c'est-à-dire variant sur une
échelle de longueur beaucoup plus grande que la constante de réseau), elles
peuvent être modélisées comme des potentiels carrés avec une profondeur et
une largeur eectives. Cela parce que le nombre d'états liés dans un potentiel
est principalement déterminé par le produit de sa profondeur eective et de
sa largeur. Il convient de noter que dans ce qui suit, nous considérons des
potentiels carrés nets mais lisses, c'est-à-dire que nous négligeons la diusion
intervallée. Pour calculer la structure de bande d'un réseau bipolaire dans le
graphène, nous utiliserons un simple modèle de liaison forte aux plus proches
voisins. Considérons d'abord un atome articiel (c'est-à-dire un puits ou une
barrière carré) dans notre super-réseau. Les électrons dans le graphène sont
régis par l'équation de Dirac

(36)


ℏvF sq̂x σx + q̂y σy + U (y)I |Ψ(x, y)⟩ = E |Ψ(x, y)⟩ ,
(1.28)
où nous avons substitué les termes de vecteur d'onde qx et qy dans l'équa-
tion (1.16) par les opérateurs q̂x = −i∂x et q̂y = −i∂y respectivement et
le potentiel de puits ou de barrière unidimensionnel est U (y). Notez que
dans cette section, nous omettons l'indice de vallée (s) en exposant des opé-
rateurs et des fonctions d'onde par souci de concision. En présence d'un
potentiel unidimensionnel, nous pouvons écrire la fonction d'onde comme
|Ψ(x, y)⟩ = eiqx x |Ψ(y, qx )⟩ ce qui simplie le problème
h i
ĤG (qx ) + U (y)I |Ψ(y, qx )⟩ = E(qx ) |Ψ(y, qx )⟩ , (37)

(1.29)
où ĤG (qx ) = ℏvF sqx σx + q̂y σy est l'équation de Dirac unidimensionnelle

eective, qui agit sur la fonction d'onde à deux composantes |Ψ(y, qx )⟩. Tous
les paramètres restants ont été précédemment dénis dans le chapitre. Le
potentiel carré U (y) est déni

14
Figure 4  Schéma de deux réseaux planaires de guides d'ondes bipolaires
dans le graphène créés par des nanotubes de carbone polarisés en alternance.
Le potentiel électrostatique créé par les tensions de grille appliquées est re-
présenté ci-dessous. (a) Si les nanotubes de carbone sont espacés de manière
équidistante, le super-réseau possédera un plan de réexion autour du centre
de tout puits ou barrière. (b) Si ce n'est pas le cas, la symétrie de réexion
est brisée.

===== Page 12 =====

comme
(
U0 | y| ≤ W/2
U (y) = , (38)
0, ailleurs
(1.30)
où W est la largeur du potentiel et U0 = Uw < 0 pour un puits et
U0 = Ub > 0 pour une barrière. Les valeurs propres de l'équation (1.29) sont
discutées dans la référence [48].
Dans la gure 5, nous traçons les spectres d'énergie pour divers puits
et barrières carrés. Chaque potentiel a la même largeur et ne contient que
quelques modes à l'intérieur. Cela se produit lorsque le produit normalisé
de la hauteur du potentiel et de la largeur (|U0 |W/ℏvF ) est de l'ordre de
l'unité. En eet, des guides d'électrons lisses à quelques modes dans le gra-
phène peuvent être réalisés expérimentalement en utilisant des nanotubes de
carbone comme grilles supérieures [49].
Comme on peut le voir sur la gure 5, la variation des intensités de
potentiel du puits et de la barrière entraîne des vitesses de groupe diérentes
au point de croisement. Ainsi, le croisement formé par un puits et une barrière
isolés peut être commuté entre le type-I sous-critiquement incliné, le type-II
super-critiquement incliné ou le type-III critique en changeant simplement
l'intensité du potentiel du puits et de la barrière. Cependant, la structure de
bande est nettement

===== Page 13 =====

===== Page 14 =====
modiée lorsque le puits et la barrière sont rapprochés pour former un
guide d'ondes bipolaire. Dans ce cas, le chevauchement entre les états de bar-
rière et de puits conduit à un anti-croisement (pseudogap) apparaissant au
croisement de bande original (voir Fig. 1.6(a)) [50]. Pour les guides d'ondes
bipolaires créés par des grilles supérieures de nanotubes de carbone sur du
graphène, le pseudogap est de l'ordre de plusieurs térahertz (THz) [50 ; 51].
Cependant, un seul puits et une seule barrière ne constituent pas un dispositif

15
Figure 5  Le spectre énergétique des états connés dans un puits de ten-
sion appliquée U0 = −120meV (bleu) et trois barrières de forces 90meV (gris
plein), 60meV (gris pointillé) et 45meV (gris pointillé) dans le graphène -
dans chaque cas la largeur du puits/barrière est W = 15nm. Les dispersions
de bande sont esquissées en unités d'énergie E en milliélectronvolts (meV)
et de vecteur d'onde associé au mouvement le long des potentiels qx norma-
lisé par la largeur du puits/barrière. Les croisements de bande formés sont
respectivement de type-I, -III et -II. Ajuster la hauteur de la barrière et la
profondeur du puits change l'inclinaison des croisements de bande. Les ré-
gions grises contiennent des états continus à l'extérieur des potentiels guides.

macroscopique, et pour des applications THz réalistes, une question impor-

tante doit être répondue : Comment la structure de bande d'un seul guide
d'ondes bipolaire est-elle modiée lorsqu'il est placé dans un super-réseau ?
L'hamiltonien d'un super-réseau bipolaire planaire ĤSL (y, qx ) peut être
écrit comme

Ny
X
ĤSL (y, qx ) = ĤG (qx ) + [Ub (y − ny L) I + Uw (y − ny L − δ0 ) I] , (39)
ny =1

(1.31)
où Ub (y) et Uw (y) sont des barrières et des puits individuels (dénis par
l'équation (1.30)), centrés respectivement aux positions ny L et ny L + δ0 , où
L est la largeur de la maille unitaire du super-réseau et δ0 est la distance
entre le centre d'un puits et d'une barrière dans la maille unitaire. Étant
donné que la dispersion de tout potentiel guide réaliste est déterminée par le
produit de la profondeur et de la largeur du potentiel, nous xons désormais
la largeur de toutes les barrières et de tous les puits à la même valeur. Nous
notons que la coordonnée du puits doit satisfaire δ0 > W et que la largeur
totale de la maille unitaire doit être supérieure à la somme de la largeur du
puits et de la barrière L > 2W . De manière similaire à la division des niveaux
d'énergie atomique dans la formation d'un cristal, le rapprochement de Ny
guides d'ondes bipolaires entraîne la division de chaque niveau d'énergie du
puits et de la barrière en Ny sous-niveaux. Chaque sous-niveau correspond
à un eky quantié particulier. Dans la limite où Ny devient grand, e ky peut
être traité comme un vecteur d'onde de super-réseau continu sur un pied
d'égalité avec qx , le vecteur d'onde le long des potentiels guides.
Nous allons maintenant calculer les spectres d'énergie du réseau bipolaire
au voisinage des croisements de bande des puits et des barrières en utilisant
un simple modèle de liaison forte aux plus proches voisins. Les fonctions
d'onde du super-réseau peuvent être exprimées comme une combinaison li-
néaire de fonctions de Bloch individuelles de puits et de barrières :

16
Figure 6  La structure de bande et les schémas de (a) un guide d'ondes
bipolaire unique, (b) un réseau bipolaire avec symétrie de réexion, (c) un
réseau bipolaire sans symétrie de réexion. Le réseau bipolaire avec ou sans
symétrie de réexion possède un cône de Dirac sans gap au centre de la
structure de bande électronique. De plus, le réseau bipolaire avec symétrie
de réexion héberge des cônes de Dirac inclinés sans gap satellites tandis que
le réseau bipolaire sans symétrie de réexion ouvre un gap dans ces cônes de
Dirac. Les structures de bande sont tracées en termes d'énergie (E ) en unités
de milliélectronvolts (meV) et de vecteur d'onde le long des potentiels qx
normalisé par les largeurs de puits/barrière W = 15nm. Dans tous les cas, les
potentiels de barrière et de puits appliqués sont respectivement Ub = 90meV
et Uw = −120meV. Pour le réseau bipolaire, la largeur de la maille unitaire
est L = 90nm, la position du puits dans la maille unitaire est déterminée
par le paramètre δ0 = 45nm dans le panneau (b) et δ0 = 48nm dans le
panneau (c). Les zones grises dans les panneaux (b) et (c) correspondent
aux énergies et aux vecteurs d'onde qui supportent des solutions d'ondes
planes à travers tout le potentiel. La périodicité du super-réseau produit un
vecteur d'onde supplémentaire e ky ≤ π/L où L est la taille de la maille
unitaire. Les structures de bande des panneaux (b) et (c) ont été calculées
à l'aide d'un modèle de matrice de transfert. Ces panneaux achent des
projections orthographiques des structures de bande vues le long de l'axe e ky .
Les bords de bande sont représentés par des lignes pleines (e ky = 0) ou des
lignes pointillées (e
ky = ±π/L) avec des valeurs intermédiaires ombrées en
bleu.

E 1 X iny eky L
Θb (y, qx , e
ky ) = p e |Ψb (y − ny L, qx )⟩ , (40)
Ny ny

(1.32)
et
E 1 X iny eky L
Θw (y, qx , e
ky ) = p e |Ψw (y − ny L − δ0 , qx )⟩ , (41)
Ny ny

(1.33)
où eky ≤ π/L est le vecteur d'onde du super-réseau et les fonctions de
puits et de barrières |Ψw (y, qx )⟩ et |Ψb (y, qx )⟩ sont les solutions de l'équation
(1.29) avec les potentiels Uw (y) et Ub (y) respectivement pour un vecteur
d'onde donné le long des potentiels guides qx . Nous avons l'intention de
modéliser la structure de bande au voisinage des croisements de bande d'un
puits individuel

17
===== Page 15 =====
et barrière (voir Fig. 1.5). Ces croisements se produisent à des vecteurs
d'onde qx = Qx = se |Qx | où se = 1 ou −1. Par conséquent, nous restreignons
notre analyse à proximité des points de croisement (c'est-à-dire petit qex =
qx − Qx ) où la bande électronique d'un puits individuel
h i
Ew (qx ) = ⟨Ψw (y, qx )| ĤG (qx ) + Uw (y)I |Ψw (y, qx )⟩ , (42)
(1.34)
ou barrière
h i
Eb (qx ) = ⟨Ψb (y, qx )| ĤG (qx ) + Ub (y)I |Ψb (y, qx )⟩ , (43)
(1.35)
est approximée comme une dispersion linéaire Ew (e qx ) ≈ E0 + seℏvw qex et
qx ) ≈ E0 + seℏvb qex pour petit qex où vw > 0 et vb < 0 et E0 est une énergie
Eb (e
de décalage (voir Fig. 1.5).
Les valeurs propres du super-réseau en fonction de qx (vecteur d'onde
le long du potentiel guide) et e ky (le vecteur d'onde du super-réseau) sont
déterminées à partir de l'équation séculaire det H(qx , e ky ) = 0,
 
ky ) − E(qx , e
où les éléments de l'hamiltonien de Bloch 2 × 2 sont dénis par le produit
scalaire Hαβ (qx , e
ky ) = ⟨Θα (y, qx , e
ky )| ĤSL (y, qx ) |Θβ (y, qx , e ky )⟩, où α, β = w
ou b qui sont arrangés dans une matrice comme :
" #
H ww (q x , ky ) H wb (qx , k y )
(44)
e e
H(qx , e
ky ) = .
Hbw (qx , e ky ) Hbb (qx , e ky )
(1.36)
Pour une maille unitaire contenant un potentiel bipolaire, deux structures
de bande distinctes peuvent apparaître selon sa symétrie. Premièrement,
si le réseau bipolaire possède un plan de réexion, sa structure de bande
supportera des cônes de Dirac sans gap. Deuxièmement, l'absence d'un plan
de réexion ouvre le gap dans ces cônes de Dirac.

Émergence des cônes de Dirac inclinés

Initialement, nous considérons un réseau bipolaire idéalisé avec des puits
et des barrières séparés de manière équidistante (δ0 = L/2, voir Fig. 1.6(b)),
les cônes de Dirac se produisent à eky = 0 et qx = Qx = se |Qx |. Les éléments de
l'hamiltonien de Bloch peuvent être trouvés en utilisant les équations (1.31)-
(1.36). Nous concentrons ensuite notre analyse sur le point de croisement
qx = Qx (c'est-à-dire petit qex ) en utilisant les approximations linéaires pour
les dispersions de bande des puits individuels Ew (e qx ) et des barrières Eb (e
qx )
et l'approximation suivante pour les intégrales de saut aux plus proches
voisins :

18
 ⟨Ψb (y − ny L, Qx )| ĤSL (y, Qx ) |Ψw (y − ny L − L/2, Qx )⟩ = se
se
t, (45)

(1.37)
et

⟨Ψb (y − ny L, Qx )| ĤSL (y, Qx ) |Ψw [y − (ny − 1)L − L/2, Qx ]⟩ = −set. (46)

se

(1.38)

===== Page 16 =====

Ces paramètres de saut peuvent être ajustés pour correspondre à la struc-
ture de bande calculée par un modèle de matrice de transfert accompagnant
(dérivé dans l'annexe A.1) et ils peuvent être dérivés analytiquement pour le
cas de profondeur de puits et de hauteur de barrière égales Uw = −Ub (voir
annexes A.2 et A.3). L'hamiltonien de Bloch au voisinage du croisement de
bande, c'est-à-dire petit qex , prend la forme
" #
ℏv q s t (1 − e ie
ky L )
w x
(47)
e
H(e
qx , e
ky ) = se .
e
t(1 − e−iky L )
se ℏvb qex
e

(1.39)
Nous notons que l'énergie de décalage E0 a été omise par souci de conci-
sion. Le développement de Taylor autour de e
ky = qey = 0 et la transformation
unitaire H(e
e q) = UH(eq)U où
†

 
1 1 1
U= √ , (48)
2 1 −1
(1.40)
donne :

(49)


H(e
e q) = seℏv γηe
qx I + ηe
qx σx + se
qy σy ,
(1.41)
où qe = (e qx , qey ), v = −e
tL/ℏ est la vitesse de Fermi eective, γ =
(vw + vb )/(vw − vb ) et η = (vw − vb )/2v sont respectivement les paramètres
d'inclinaison et d'anisotropie. L'équation (1.41) est le bien connu hamiltonien
de cône de Dirac incliné discuté dans la section 1.1.3. Il convient de noter que
le pseudo-gap qui était présent dans le potentiel bipolaire individuel a été
complètement fermé. Ceci est un résultat de la symétrie de réexion présente
dans le réseau. Il convient également de noter que les conditions aux limites
des réseaux bipolaires nis et innis sont diérentes. À savoir, dans les ré-
seaux nis, la fonction d'onde doit décroître à l'extérieur des puits les plus

19
externes, tandis que dans le cas inni, le système est soumis aux conditions
aux limites de Born-von Karman. Par conséquent, dans les systèmes nis,
aucun mode guidé n'existe dans la région où |E| > ℏvF |qx | (régions grises de
la Fig. 1.6). Inversement, dans le cas inni, des modes guidés sont supportés
dans cette région.
La présence de cônes de Dirac inclinés est conrmée par un calcul de ma-
trice de transfert dérivé dans l'annexe A.1 et montré dans la Fig. 1.6(b). Nous
notons que ce modèle de liaison forte aux plus proches voisins correspond le
plus précisément au modèle de matrice de transfert lorsque les états liés dans
les puits et les barrières sont bien localisés dans une seule maille unitaire. À
l'extérieur du puits et de la barrière, les fonctions d'onde Ψw/b (y, qx ) au
vecteur d'onde de croisement qx = Qx avec énergie E0p sont proportionnelles
àe κ
e 0 y pour y < 0 ou e−e
κ 0 y pour y > 0 où κ
e0 = (1/ℏvF ) (ℏvF Qx )2 − E02 . Par
conséquent, le modèle de liaison forte aux plus proches voisins correspond
le plus ecacement aux calculs de matrice de transfert lorsque les fonctions
d'onde sont bien localisées dans une seule maille unitaire κ e0 L ≫ 1.
Contrairement aux matériaux de Dirac inclinés formés par des réseaux
cristallins, le paramètre d'inclinaison γ = (vw + vb )/(vw − vb ) apparaissant
dans l'équation (1.41) peut être ajusté en faisant varier les tensions de grille
appliquées du réseau bipolaire. Cela parce que faire varier les tensions de
puits et de barrière appliquées déplace le vecteur d'onde de croisement Qx ,
faisant ainsi varier les vitesses de groupe vb et vw au point de croisement
(comme démontré pour le cas de puits et barrières individuels*

===== Page 17 =====

dans la Fig. 1.5). Nous conrmons la prédiction d'un paramètre d'incli-
naison accordable γ en calculant la structure de bande d'un réseau bipolaire
à l'aide du modèle de matrice de transfert (voir Fig. 1.7). Comme on peut
le voir sur la Fig. 1.7, faire varier la hauteur de barrière Ub pour une pro-
fondeur de puits xe de Uw entraîne le passage des cônes du type-I (|γ| < 1)
au type-II (|γ| > 1) en passant par le cas de type-III à inclinaison critique
(|γ| = 1). Nous notons que dans le cas de type-III à inclinaison critique, une
branche de la dispersion de bande électronique apparaît quadratique plutôt
que linéaire. Ces croisements de bande sont parfois appelés points de Dirac
aux trois quarts qui ont des propriétés intéressantes telles que des niveaux de
Landau avec une énergie qui évolue comme la puissance quatre cinquièmes
de l'intensité du champ magnétique B et de l'indice de niveau de Landau n,
c'est-à-dire En ∝ (nB)4/5 [52, 53].
La nature relativiste des électrons dans les matériaux de Dirac inclinés
entraîne leurs excellentes propriétés électroniques. Cependant, pour les ap-
plications de dispositifs optoélectroniques, il serait souhaitable d'ouvrir un
gap accordable de l'extérieur dans le cône. Pour introduire le gap, nous de-
vons briser la symétrie de réexion du réseau bipolaire. Il convient de noter

20
que les études précédentes sur les super-réseaux de graphène se sont limitées
à des puits/barrières périodiques [54, 55, 56, 57], sinusoïdaux [58, 59], ou à
des potentiels pairs/impairs périodiques [60] qui ont tous une symétrie de
réexion, et n'ouvrent donc pas de gap dans les cônes de Dirac.

Figure 7  Projections orthographiques des structures de bande de trois

réseaux bipolaires vues le long du vecteur d'onde du super-réseau e ky . Les
bords de bande sont représentés par des lignes pleines (ky = 0) ou des lignes
e
pointillées (e
ky = ±π/L) avec des valeurs intermédiaires ombrées en bleu
où L est la taille de la maille unitaire du super-réseau. Les structures de
bande ont été calculées à l'aide d'une matrice de transfert et sont tracées en
termes d'énergie E en unités de milliélectronvolts (meV) et de vecteur d'onde
le long des potentiels guides qx normalisé par la largeur du puits/barrière
W = 15nm. Dans chaque tracé, la longueur de la maille unitaire est L =
90nm, l'emplacement du puits dans la maille unitaire est déterminé par δ0 =
45nm et la profondeur du puits est xée à Uw = −130meV. L'inclinaison
des croisements de bande varie du type-I avec une hauteur de barrière Ub =
90meV dans le panneau (a), en passant par le type-III à environ Ub = 55meV
dans le panneau (b) et jusqu'au type-II à Ub = 45meV dans le panneau
(c). Les zones grises correspondent aux énergies et aux vecteurs d'onde qui
supportent des solutions d'ondes planes à travers tout le potentiel.

===== Page 18 =====

Émergence des cônes de Dirac inclinés massifs
On peut envisager un réseau bipolaire créé en sandwichant une feuille
de graphène entre deux réseaux planaires de nanotubes, le réseau supérieur
étant polarisé à une polarité, le réseau inférieur à la polarité opposée. La
position relative de ces deux réseaux (paramétrée par δ0 ) sera xée après la
fabrication du dispositif. Dans un dispositif réaliste, il ne sera pas possible
d'aligner exactement les deux réseaux de telle manière que chaque tube soit
également séparé, en général les deux réseaux seront séparés par une distance
arbitraire (δ0 ̸= L/2). Suivant la même prescription que la section 1.2.3, mais
en utilisant les nouvelles intégrales de chevauchement :

⟨Ψb (y − ny L, Qx )| ĤSL (y, Qx ) |Ψw (y − ny L − δ0 , Qx )⟩ = s⃗se

t1 , (50)

(1.42)
et

⟨Ψb (y − ny L, Qx )| ĤSL (y, Qx ) |Ψw [y − (ny − 1)L − δ0 , Qx ]⟩ = s⃗ste2 , (51)

21
(1.43)
où e
t1 et e
t2 peuvent être ajustés à la structure de bande calculée par
un modèle de matrice de transfert ou peuvent être dérivés analytiquement
pour le cas de profondeur de puits et de hauteur de barrière égales (Uw =
−Ub ) comme on le voit dans les annexes A.2 et A.3. Notez que dans le cas
symétrique (δ0 = L/2), les paramètres de saut se réduisent à e t2 = −e t.
t1 = −e
L'hamiltonien de Bloch résultant est :
" #
ℏv q s( t + t eie
ky L )
w x 1 2
(52)
e e
H(eqx , e
ky ) = se ,
e
s(e t2 e−iky L )
t1 + e ℏvb qex
e

(1.44)
où nous avons encore une fois omis l'énergie de croisement E0 par souci
de concision. Après un développement de Taylor du second ordre autour de
ky = 0 et un changement de base via la matrice unitaire :
qey = e
!
1 e−iky l eiky l
e e
U= √ , (53)
2 e−iky l −eiky l
e e

(1.45)
avec
" p #
t2 + −e t1 e
t2
(54)
e
l= L,
2(e
t1 + e
t2 )
(1.46)
nous obtenons

∆
(55)


H(q̃)
e s σz + seℏv γηe
= se qx I + ηe
qx σx + se
qy σy ,
2
(1.47) p
où v = − −e t2 L/ℏ est la vitesse de Fermi eective. L'équation (1.47)
t1 e
a la forme de l'hamiltonien de cône de Dirac incliné avec un gap de bande
|∆| = 2 |e t2 |. Nous soulignons que le paramètre d'inclinaison γ est modié
t1 + e
par les tensions de grille appliquées comme discuté précédemment pour le
réseau bipolaire avec symétrie de réexion (voir Fig. 7). Il convient de no-
ter que les intégrales de chevauchement décroissent exponentiellement avec
l'augmentation de la taille de la maille unitaire L. Ceci est xé pendant le pro-
cessus de fabrication. Cependant, les intégrales de chevauchement peuvent
également

===== Page 19 =====

être contrôlées en faisant varier les tensions de grille appliquées des puits
et des barrières, ce qui signie que le gap de bande est néanmoins hautement
accordable après fabrication.

22
 Estimons maintenant la taille du gap de bande pour un dispositif réaliste.
Pour le cas de forces de puits et de barrière égales (|Uw | = |Ub |), le croisement
de bande se produit à énergie nulle (E0 = 0). Dans ce cas, les fonctions
d'onde de puits et de barrière et les paramètres de saut t̃1 et t̃2 peuvent être
trouvés analytiquement (voir annexes A.2 et A.3). Dans ce cas, le gap prend
l'expression analytique simple
 
|Qx ||L − 2δ0 | −|Qx |L/2
|∆| = ∆0 sinh e , (56)
2
(1.48)
où

4ℏ3 vF3 |Qx |Q

e2
∆0 = 2 e|Qx |W , (57)
Ub (1 + |Qx |W )
(1.49)
où |Q
qx | satisfait l'équation transcendante Q = −|Qx | tan(QW ) avec Q =
e e e
(1/ℏvF ) Ub2 − (ℏvF Qx )2 . Nous considérons les largeurs eectives de puits
et de barrière W = 10nm (qui sont trouvées expérimentalement dans la Réf.
[49]) avec une largeur de cellule de super-réseau L = 50nm où la barrière
est centrée au bord de la maille unitaire et le puits à δ0 = 27nm. Nous
considérons des intensités de potentiel réalistes, par exemple, Ub = −Uw =
200meV[49] ce qui donne un gap de 0,49THz. Diminuer la tension de 20% à
160meV augmente le gap à 1,03THz. Ces gaps sont en accord à 1% près avec
un calcul de matrice de transfert. Nos calculs indiquent qu'il sera possible
d'ajuster le gap de bande dans tout le régime térahertz 0,1-2THz.

Cônes de Dirac centraux sans gap

Bien que nous ayons discuté du rôle de la géométrie du super-réseau
dans l'ouverture de gaps de bande dans les cônes de Dirac, nous soulignons
que la structure de bande complète du réseau bipolaire reste sans gap. Ceci
est dû aux cônes de Dirac qui existent au centre de la structure de bande
électronique pour toutes les géométries de super-réseau [54, 55, 60] (voir Fig.
1.6(c)). À cet égard, les cônes de Dirac inclinés et avec gap précédemment
discutés sont des satellites d'un cône de Dirac central sans gap. Nos calculs
numériques indiquent que pourvu que le prol de tension du super-réseau
soit choisi de telle sorte que les cônes de Dirac satellites avec gap soient
sous-critiquement inclinés (type-I, |γ| < 1), les cônes de Dirac centraux et
satellites sont bien séparés dans l'espace des impulsions. Dans ce cas, le cône
de Dirac central (désigné par l'indice c) est bien décrit par l'hamiltonien

Hc (e
qc,x , qec,y ) = ℏ(svc,x qec,x σx + vc,y qec,y σy ), (58)
(1.50)

23
 où les vecteurs d'onde développés sont qec,x = qx et qec,y = qey et vc,x , vc,y ≤
vF sont des vitesses de groupe qui peuvent être déterminées à partir d'un
modèle de matrice de transfert. Nos calculs numériques indiquent que le
décalage d'énergie du cône de Dirac central est

===== Page 20 =====

égal au potentiel moyen du réseau bipolaire W (Ub +Uw )/L. Cette énergie
de décalage a été omise de l'équation (1.50) par souci de concision.

1.2.6 Conclusions
Dans cette section, nous avons proposé une méthode pour concevoir des
cônes de Dirac inclinés avec et sans gap dans un super-réseau de graphène
latéral. Contrairement aux monocouches atomiques cristallines, la structure
de bande électronique d'un super-réseau de graphène peut être modiée en
faisant varier le prol de tension appliqué - cela fournit un moyen pratique de
contrôler le paramètre d'inclinaison et le gap de bande des cônes de Dirac.
Nos calculs démontrent également que l'amplitude du gap de bande local
dans les cônes de Dirac avec gap peut être amenée à tomber dans le régime
térahertz. Une voie de travail future consistera à étudier les applications
térahertz facilitées par les transitions interbandes à travers les cônes de Dirac
avec gap et inclinés.
Cette capacité à modier le paramètre d'inclinaison des cônes de Dirac
dans le super-réseau bipolaire de graphène ouvre un nouveau degré de liberté
qui peut être exploité pour des applications de dispositifs. Par exemple, en
présence d'un champ magnétique orienté normalement au super-réseau, la
structure de bande électronique se quantie en niveaux de Landau discrets.
En faisant varier le paramètre d'inclinaison et le gap de bande des cônes de
Dirac satellites au moyen de la variation des tensions de grille électrostatiques
appliquées, il sera possible d'ajuster le spectre des niveaux de Landau. Cela
conduit à son tour à des spectres de niveaux de Landau accordables par grille
qui sont discutés au chapitre 4.
Bien que ce travail se soit limité à l'étude de super-réseaux latéraux
unidimensionnels dans le graphène, nous notons que les super-réseaux de
graphène bidimensionnels fournissent également une plateforme viable pour
réaliser des cônes de Dirac sur mesure [61]. En étudiant de nouvelles géomé-
tries de super-réseaux bidimensionnels, il pourrait également être possible
de réaliser des cônes de Dirac satellites inclinés et avec gap similaires à ceux
étudiés dans ce travail. Il est également noté que bien que dans cette sec-
tion nous ayons considéré des super-réseaux latéraux appliqués au graphène,
il pourrait également être possible de réaliser des cônes de Dirac avec gap
et inclinés dans des super-réseaux appliqués à d'autres systèmes bidimen-
sionnels tels que le graphène bicouche [62] ou éventuellement, des matériaux

24
bidimensionnels qui hébergent déjà des cônes de Dirac inclinés dans la struc-
ture de bande électronique, c'est-à-dire le borophène 8-P mmn [63].

1.3 Cônes de Dirac inclinés dans les réseaux planaires de na-

notubes de carbone

1.3.1 Introduction
Les nanotubes de carbone (CNT) sont connus comme l'allotrope uni-
dimensionnel du carbone en raison de leur rapport d'aspect considérable
déterminé par des diamètres de tube de l'ordre

===== Page 21 =====

de nanomètres et des longueurs pouvant dépasser des mètres [64, 65, 66].
D'un point de vue théorique, un nanotube de carbone peut être compris
comme une feuille de graphène enroulée (voir Fig. 8). Changer l'axe d'enrou-
lement d'un CNT entraîne des structures de tube radicalement diérentes
ou des chiralités, déterminant ainsi si le nanotube est un métal ou un semi-
conducteur [67, 68]. Depuis leur découverte, les techniques expérimentales
de production de CNT n'ont cessé de s'améliorer, avec des progrès incluant
la sélection de chiralité [69, 70, 71, 72] et la capacité à faire pousser des
tubes hautement alignés [73, 74, 75]. Tout comme l'utilisation du graphène
et d'autres monocouches cristallines dans les hétérostructures van der Waals
sur mesure [76, 77], les progrès dans la croissance et la sélection des CNT
pourraient un jour permettre aux CNT d'être utilisés comme constituants
dans de nouvelles structures composites maintenues ensemble par des forces
de van der Waals.
Dans ce travail, nous étudions un réseau van der Waals planaire de nano-
tubes de carbone à petit gap (quasi-métalliques) (Section 1.3.2). Nous notons
qu'il s'agit d'une structure bidimensionnelle planaire qui est signicative-
ment diérente des réseaux de nanotubes de carbone alignés verticalement
tridimensionnels plus communément connus [78, 79, 80]. Ce système a été
précédemment étudié à l'aide de la théorie de la fonctionnelle de la densité
(DFT) où il a été constaté que la structure de bande basse énergie est sans
gap et semble héberger des cônes super-critiquement inclinés au voisinage
du niveau de Fermi [20]. Pour conrmer la présence de ces croisements de
bande et prouver l'existence de cônes de Dirac de type-II, nous abordons
le problème théoriquement, en calculant la structure de bande électronique
basse énergie (Section 1.3.3) et en extrayant l'hamiltonien de cône de Dirac
incliné (Section 1.3.4).

25
Figure 8  Un nanotube de carbone peut être compris comme une feuille
de graphène enroulée. L'axe d'enroulement est déterminé par le vecteur chi-
ral Ch = (n, m) = na1 + ma2 qui est déni comme des multiples entiers
des vecteurs de réseau du graphène a1 et a2 . L'axe d'enroulement ou chiral
détermine si le nanotube de carbone sera zig-zag (n, 0), armchair (n, n) ou
chiral (n, m).

===== Page 22 =====

Théorie de fond : Nanotubes de carbone quasi-métalliques
Les nanotubes de carbone peuvent être modélisés comme une feuille de
graphène enroulée, c'est-à-dire sélectionner un atome à un point de réseau
R et l'amener à un autre site atomique situé à R + na1 + ma2 (voir Fig. 8).
En enroulant la feuille de graphène, nous avons construit un nanotube de
carbone dans lequel tout atome peut être mappé sur lui-même par le vecteur
chiral Ch = (n, m) = na1 + ma2 qui dénit la chiralité du tube. La chiralité
n'a pas seulement un eet structural sur un CNT, elle aecte également
radicalement sa structure de bande électronique. En enroulant une feuille de
graphène, nous appliquons une condition aux limites périodique au réseau
de graphène, ce qui entraîne une quantication du vecteur d'onde cristallin
Ch · kµ = 2πµ où µ est un entier. La dispersion de bande d'un CNT consiste
donc en des tranches unidimensionnelles quantiées de la structure de bande
électronique bidimensionnelle du graphène. Chaque chiralité de CNT a un
vecteur chiral distinct Ch entraînant des quantications de vecteur d'onde
et des structures de bande électronique diérentes. Par exemple, les CNT
armchair (n, n) sont métalliques tandis que les CNT zig-zag (n, 0) sont quasi-
métalliques si n est un multiple de trois et semi-conducteurs sinon.
Tout au long de cette section, nous nous concentrerons sur un sous-
ensemble de nanotubes de carbone quasi-métalliques (QMCNT) avec la chi-
ralité (n, 0) où n est un multiple de trois. Dans cette section, nous spécierons
une géométrie de CNT spécique (12, 0) (voir Fig. 9(a)) mais nous soulignons
que nos méthodes peuvent être utilisées pour tous les QMCNT. Pour la chi-
ralité (12, 0), le vecteur chiral prend
√ la forme Ch = 12a1 = 12ax̂, le vecteur
de réseau le long du tube est T = 3aŷ et il y a 48 atomes par maille unitaire
répartis entre deux sous-réseaux (A et B ) que nous énumérons comme jA ou
jB où j = 1, 2, ..., 24. Pour obtenir la structure de bande des nanotubes de
carbone, nous utilisons la méthode de repliement de zone [68] dans laquelle
nous prenons la structure de bande du graphène vierge et quantions le vec-
teur d'onde le long de l'axe chiral du nanotube Ch · kµ = 2πµ. Le vecteur
d'onde quantié résultant est
 µπ 
kµ = , ky , (59)
6a

26
(51)
où il y a 24 modes énumérés comme µ = 0, 1, ...23 et |ky | ≤ π/T où
T = |T|. La structure de bande d'un QMCNT (12, 0) peut être obtenue
en quantiant le vecteur d'onde dans la structure de bande du graphène
donnée dans l'équation (14) (voir Fig. 9(b)). En particulier, on peut voir
que les modes basse énergie dégénérés µ = 8 et 16 intersectent les cônes de
Dirac aux points K/K′ ce qui signie que selon les hypothèses actuelles, la
structure de bande est sans gap ou métallique (voir ligne pointillée dans la
Fig. 9(c)).
Dans le graphène vierge, tous les atomes de carbone se trouvent dans le
plan de la feuille de graphène et sont également séparés et ont donc le même
paramètre de saut t. En revanche, dans les QMCNT, certains atomes sont
séparés le long de la longueur du tube (voir atomes jA et jB ) tandis que
d'autres sont séparés azimutalement autour du tube (voir jA et (j ± 1)B ).
Il est clair que que ce soit par la structure incurvée des nanotubes (voir Fig.
9(d)) ou par les orientations diérentes des orbitales (voir Fig. 9(e)) que nous
devons attribuer deux

===== Page 23 =====

paramètres de saut diérents à ces cas. À partir de ce point, nous modi-
ons le modèle de graphène en utilisant deux paramètres de saut diérents :
t lorsque les atomes sont séparés le long de la longueur du tube et t′ au-
trement. Cela a pour eet de déplacer l'emplacement des cônes de Dirac
loin des points K/K′ ce qui ouvre à son tour un petit gap de bande dans
la structure de bande des QMCNT (voir ligne pleine dans la Fig. 1.9(c) et
aussi le panneau (f)). En substituant notre vecteur d'onde quantié donné
dans l'équation (1.51) dans l'hamiltonien de graphène donné dans l'équation
(1.13) avec les paramètres de saut ajustés, on obtient l'hamiltonien de Bloch
pour le mode µ
 
0 fµ (ky )
Hµ (ky ) = ∗ , (60)
fµ (ky ) 0
(1.52)

===== Page 23 =====

paramètres de saut diérents à ces cas. À partir de ce point, nous modi-
ons le modèle de graphène en utilisant deux paramètres de saut diérents :
t lorsque les atomes sont séparés le long de la longueur du tube et t′ au-
trement. Cela a pour eet de déplacer l'emplacement des cônes de Dirac
loin des points K/K′ ce qui ouvre à son tour un petit gap de bande dans
la structure de bande des QMCNT (voir ligne pleine dans la Fig. 1.9(c) et
aussi le panneau (f)). En substituant notre vecteur d'onde quantié donné

27
Figure 9  (a) Schéma d'un nanotube de carbone quasi-métallique (12, 0),
la maille unitaire (esquissée en couleur) contient 48 atomes de carbone éti-
quetés jA (bleu) et jB (noir) où j = 1, 2, ..., 24. (b) Dans l'approximation de
repliement de zone, la structure de bande d'un nanotube de carbone (12, 0)
peut être décrite en quantiant le vecteur d'onde du graphène (k → kµ où
µ = 0, 1, 2, ..., 23) - la structure de bande résultante consiste en 24 tranches
unidimensionnelles de la structure de bande bidimensionnelle du graphène.
La structure de bande basse énergie est décrite par les modes dégénérés µ = 8
et µ = 16 (bleu) qui intersectent les points K/K′ dans la zone de Brillouin
du graphène. (c) La structure de bande complète d'un nanotube de carbone
quasi-métallique (12, 0) est esquissée. Les paramètres de saut de liaison forte
(chevauchement des orbitales atomiques) des atomes séparés azimutalement
autour du nanotube (t′ ) sont plus grands que les paramètres de saut des
atomes séparés le long de la longueur du tube (t). Ceci est dû à (d) la réduc-
tion de la séparation des plus proches voisins atomiques due à la courbure
du nanotube ou (e) la rotation des orbitales atomiques. (f) Le paramètre de
saut modié t′ ̸= t déplace les cônes de Dirac loin des points K/K′ ce qui
signie que la structure de bande possède un mini-gap quasi-métallique.

dans l'équation (1.51) dans l'hamiltonien de graphène donné dans l'équation

(1.13) avec les paramètres de saut ajustés, on obtient l'hamiltonien de Bloch
pour le mode µ
 
0 fµ (ky )
Hµ (ky ) = ∗ , (61)
fµ (ky ) 0
(1.52)

===== Page 24 =====

avec
 µπ 
fµ (ky ) = t + 2t′ e−iky T /2 cos . (62)
12
(1.53)
La diagonalisation de l'hamiltonien donne la structure de bande pour le
mode µ comme Eµ,± (ky ) = ±|fµ (ky )| avec les états

±eiφµ (ky )
 
1
|Ψµ,± (ky )⟩ = √ , (63)
2 1
(1.54)
où eiφµ (ky ) = fµ (ky )/|fµ (ky )|. On peut voir à partir de la structure de
l'équation (1.53) que chaque mode µ = 1, 2, ..., 11 est dégénéré avec les modes
µ = 23, 22, ..., 13 tandis que les modes µ = 0 et µ = 12 sont non dégénérés.

28
Figure 10  (a) Schéma d'un nanotube de carbone quasi-métallique (12, 0),
la maille unitaire (esquissée en couleur) contient 48 atomes de carbone éti-
quetés jA (bleu) et jB (noir) où j = 1, 2, ..., 24. (b) Dans l'approximation de
repliement de zone, la structure de bande d'un nanotube de carbone (12, 0)
peut être décrite en quantiant le vecteur d'onde du graphène (k → kµ où
µ = 0, 1, 2, ..., 23) - la structure de bande résultante consiste en 24 tranches
unidimensionnelles de la structure de bande bidimensionnelle du graphène.
La structure de bande basse énergie est décrite par les modes dégénérés µ = 8
et µ = 16 (bleu) qui intersectent les points K/K′ dans la zone de Brillouin
du graphène. (c) La structure de bande complète d'un nanotube de carbone
quasi-métallique (12, 0) est esquissée. Les paramètres de saut de liaison forte
(chevauchement des orbitales atomiques) des atomes séparés azimutalement
autour du nanotube (t′ ) sont plus grands que les paramètres de saut des
atomes séparés le long de la longueur du tube (t). Ceci est dû à (d) la réduc-
tion de la séparation des plus proches voisins atomiques due à la courbure
du nanotube ou (e) la rotation des orbitales atomiques. (f) Le paramètre de
saut modié t′ ̸= t déplace les cônes de Dirac loin des points K/K′ ce qui
signie que la structure de bande possède un mini-gap quasi-métallique.

Les cas µ = 8 et 16 correspondent aux modes basse énergie qui décrivent

la structure de bande près du niveau de Fermi (voir ligne pleine dans la
Fig. 1.9(c)). En particulier, ces modes révèlent un mini-gap quasi-métallique
∆m = 2|t − t′ | à ky = 0. Comme il est de coutume avec les modèles de
liaison forte, nous ajustons les paramètres de saut à l'expérience : pour les
CNT, le paramètre de saut prend la valeur t ≈ −2.60eV et le mini-gap
est ∆m ≈ 3|t|a2nn /16R2 où R = |Ch |/2π est le rayon du tube [28]. Pour un
QMCNT (12, 0), le mini-gap est ∆m ≈ 45meV [28 ; 81 ; 82] que nous obtenons
en xant t′ = (1 + ϵ)t où ϵ ≈ 0.0086.
La fonction d'onde spineur donnée dans l'équation (1.54) est écrite dans
la base des sommes de Bloch du graphène |ΦA (k)⟩ et |ΦB (k)⟩ qui sont dénies
dans les équations (1.3) et (1.4) respectivement
1 
|Ψµ,± (ky )⟩ = √ ±eiφµ (ky ) |Φµ,A (ky )⟩ + |Φµ,B (ky )⟩ . (64)

2
(1.55)
Notez qu'à partir de maintenant, nous omettons la fonctionnalité de la
variable de position r par souci de concision. En mettant en ÷uvre l'approxi-
mation de repliement de zone, nous avons quantié le vecteur d'onde ce qui
signie que la forme des fonctions de Bloch a changé

29
 11
!
1 X imy ky T 1 X iµmx π/6 my
E
my
E
|Φµ,A (ky )⟩ = p e √ e φ(2mx +1)A
+ei(ky T /2+µπ/12) φ(2mx +2)A
,
Ny my 24 m =0 x
(65)
(1.56)
et

11
!
1 X imy ky T 1 X iµmx π/6 my
E
i(ky T /2+µπ/12) my
E
|Φµ,B (ky )⟩ = p e √ e φ(2mx +1)B
+e φ(2mx +2)B ,
Ny m y 24 m =0 x
(66)
(1.57)
où mx est un entier qui compte les sites atomiques dans la maille unitaire,
my est un entier qui compte sur les diérentes mailles unitaires et Ny = N/24
est le nombre de mailles unitaires.

===== Page 25 =====

Notez que nous étiquetons maintenant chaque orbitale atomique selon
son indice (j = 1, 2, ..., 24), son sous-réseau (A ou B ) et son indice de maille
m m
unitaire (my ) en utilisant la notation de ket |φjAy ⟩ ou |φjBy ⟩ (voir Fig. 1.9(a)).

Modèle de liaison forte eectif pour un réseau de nanotubes de

carbone (12, 0)
Nous construisons maintenant un modèle de liaison forte pour calculer la
structure de bande électronique d'un réseau van der Waals de nanotubes de
carbone quasi-métalliques. Développer des modèles de liaison forte simulta-
nés pour diverses géométries de QMCNT est une tâche dicile, donc, dans
ce travail, nous spécions un réseau de CNT (12, 0). Nous notons que les
mêmes méthodes peuvent être utilisées pour calculer la structure de bande
électronique de réseaux de CNT (6, 0) et (18, 0). Cependant, les réseaux de
nanotubes de carbone quasi-métalliques (9, 0) et (15, 0) nécessiteront un cal-
cul séparé en raison de leur géométrie contrastée qui manque de symétrie de
rotation de 180◦ autour de l'axe du tube.
Nous considérons un réseau planaire de QMCNT (12, 0) répétés le long
d'un axe z avec chaque centre de tube séparé par la distance L (voir Fig.
1.10). La largeur de la maille unitaire du réseau

===== Page 26 =====

est déterminée par les interactions de van der Waals ; des études de théorie
de la fonctionnelle de la densité estiment que la largeur de la maille unitaire

30
Figure 11  Schémas de face (a) et de dessus (b) d'un réseau planaire
de nanotubes de carbone quasi-métalliques (12, 0). Chaque maille unitaire
(ombragée) contient 24 atomes de carbone de chaque sous-réseau A (bleu)
et B (noir) étiquetés jA et jB où j = 1, 2, ..., 24 (voir Fig. 1.9(a)). La maille
unitaire a une largeur le long de l'axe du réseau (z ) L et une distance le long
de l'axe du tube (y ) T . Un modèle de liaison forte approximatif est construit
avec plusieurs paramètres de saut inter-tube : t1 (vert), t2 (vert clair), t3
(rouge) et t4 (rose).

est approximativement quatre fois le rayon du nanotube L ≈ 4R[20]. Nous

supposons que la structure s'auto-assemble en un réseau avec un nanotube de
carbone par maille unitaire. Dans ce calcul, nous supposons que les nanotubes
de carbone sont orientés dans la maille unitaire comme esquissé dans la Fig.
1.10. Nous notons qu'il pourrait être possible de considérer des réseaux de
nanotubes dans lesquels le tube est orienté à un angle diérent - déterminer
l'angle optimal dans lequel la structure s'auto-assemble est un sujet en dehors
de la portée de ce travail et pourrait être une future voie de travail pour des
études computationnelles.
Dans la section précédente, nous avons appris que la structure de bande
électronique d'un nanotube (12, 0) consiste en 24 modes avec la structure
de bande basse énergie régie par deux modes dégénérés µ = 8 et µ = 16.
Dans cette section, nous approximons la structure de bande électronique
du réseau de nanotubes au voisinage du niveau de Fermi en ne considérant
que ces deux modes basse énergie. Pour assurer la validité de notre approche,
nous comparerons nos résultats à un modèle de liaison forte numérique décrit
dans l'annexe B.
Suivant le traitement standard de liaison forte, nous construisons la fonc-
tion d'onde comme une combinaison linéaire de sommes de Bloch de réseau
1 X imz kz L mz E
|χµ,A (ky , kz )⟩ = √ e Φµ,A (ky ) , (67)
Nz m z

(1.58)
et
1 X imz kz L mz E
|χµ,B (ky , kz )⟩ = √ e Φµ,B (ky ) , (68)
Nz m z

(1.59) où mz est un entier qui compte sur les nanotubes dans le réseau,
Nz est le nombre de nanotubes dans le réseau et kz est le vecteur d'onde
cristallin du réseau qui peut prendre des valeurs |kz | ≤ π/L. Notez que les
sommes de Bloch de réseau |χµ,A (ky , kz )⟩ et |χµ,B (ky , kz )⟩ sont une combi-
naison linéaire des sommes de Bloch QMCNT |Φm µ,A (ky )⟩ et |Φµ,B (ky )⟩ sur le
z mz

site nanotube mz qui sont elles-mêmes des combinaisons linéaires d'orbitales

31
 m ,m m ,m
atomiques |φjAy z ⟩ et |φjBy z ⟩. Pour décrire la structure de bande électro-
nique basse énergie, nous ne considérons que les sommes de Bloch de réseau
correspondant aux modes µ = 8 et µ = 16.
L'hamiltonien pour un réseau planaire de nanotubes de carbone a deux
termes

Ĥ = Ĥ CNT + Ĥ Arr , (69)

(1.60)
le premier (Ĥ CNT ) encode les paramètres de saut entre sites atomiques
sur un seul nanotube tandis que le second (Ĥ Arr ) dénit les paramètres de
saut entre nanotubes voisins. Suivant la procédure standard de liaison forte,
ces deux hamiltoniens nous permettent d'écrire un hamiltonien de Bloch avec
la forme H(ky , kz ) = HCNT (ky , kz ) + HArr (ky , kz ) avec des éléments

H(µ′ ,α),(µ,β) (ky , kz ) = ⟨χµ′ ,α (ky , kz )| Ĥ CNT + Ĥ Arr |χµ,β (ky , kz )⟩ , (70)

(1.61)

===== Page 27 =====

où α, β = A, B et µ, µ′ = 8, 16 dénissant un hamiltonien de Bloch basse
énergie 4 × 4 avec des éléments arrangés dans la structure suivante :

 
H(8,A),(8,A) H(8,A),(8,B) H(8,A),(16,A) H(8,A),(16,B)
 H(8,B),(8,A) H(8,B),(8,B) H(8,B),(16,A) H(8,B),(16,B) 
H=
H(16,A),(8,A) H(16,A),(8,B) H(16,A),(16,A) H(16,A),(16,B)  ,
 (71)
H(16,B),(8,A) H(16,B),(8,B) H(16,B),(16,A) H(16,B),(16,B)

(1.62)
où la fonctionnalité (ky , kz ) a été omise par souci de concision. Le premier
terme (HCNT (ky , kz )) contient seulement les intégrales de saut à l'intérieur
d'un seul nanotube et se lit donc simplement comme les hamiltoniens de
Bloch des deux modes basse énergie pour un seul CNT (12, 0)

 
0 f8 (ky ) 0 0
f8∗ (ky )
HCNT (ky ) =
 
H8 (ky ) 02,2 0 0 0 
= ,
02,2 H16 (ky )  0 0 0 f16 (ky )
0 0 ∗ (k )
f16 0
y
(72)
(1.63)
où Hµ (ky ) et fµ (ky ) sont dénis dans les équations (1.52) et (1.53) res-
pectivement et 0n1 ,n2 est une matrice nulle n1 × n2 . Substituer ces identités
donne l'expression

32
 t − t′ e−iky T /2
 
0 0 0
t − t′ eiky T /2
HCNT (ky ) = 
0 0 0 
.
 0 0 0 ′ −ik
t−te y  T /2

0 0 t − t′ eiky T /2 0
(73)
(1.64)
Le second terme (HArr (ky , kz )) contient les paramètres de saut entre na-
notubes voisins dans le réseau. Il est plutôt arbitraire de choisir quels sites
voisins sont choisis - par exemple ajouter plus d'éléments de saut est équi-
valent à augmenter la force d'une quantité moindre de paramètres de saut.
C'est pour cette raison que nous choisissons un petit nombre de paramètres
de saut de réseau t1 à t4 (voir Fig. 1.10) qui se verront nalement attribuer
des valeurs pour ressembler aux résultats DFT mentionnés précédemment
[20]. Les paramètres de saut inter-tube esquissés dans la Fig. 1.10 sont dé-
nis à travers l'hamiltonien

===== Page 28 =====

écrit dans la base des orbitales atomiques :

Ĥ Arr =
X h m ,m my ,mz +1 m ,m my ,mz +1 m ,m my ,mz +1
t1 ⟨φ1Ay z |φ13A ⟩ + t1 ⟨φ1By z |φ13B ⟩ + t2 ⟨φ1Ay z |φ13B ⟩
my ,mz
m ,mz m ,mz +1 m ,mz m ,mz +1 m ,mz m ,mz +1
+ t2 ⟨φ1By y
|φ13A ⟩ + t3 ⟨φ2Ay y
|φ12A ⟩ + t3 ⟨φ2By y
|φ12B ⟩
my ,mz my ,mz +1 my ,mz my ,mz +1 m ,m my ,mz +1
+ t3 ⟨φ24A |φ14A ⟩ + t3 ⟨φ24B |φ14B ⟩ + t4 ⟨φ2Ay z |φ12B ⟩
i
my ,mz my ,mz +1 my ,mz my ,mz +1 my ,mz my ,mz +1
+t4 ⟨φ2B |φ12A ⟩ + t4 ⟨φ24A |φ14B ⟩ + t4 ⟨φ24B |φ14A ⟩
+ complexe conjugué. (74)

(1.65)
En utilisant ces paramètres de saut inter-tube, le second terme de l'ha-
miltonien de Bloch basse énergie peut être trouvé comme
 
(t1 − t3 ) (t2 − t4 ) (t1 + 2t3 ) (t2 + 2t4 )
HArr (kz ) = cos(kz L) 
 (t2 − t4 ) (t1 − t3 ) (t2 + 2t4 ) (t1 + 2t3 )
(t1 + 2t3 ) (t2 + 2t4 ) (t1 − t3 ) (t2 − t4 )  . (75)


(t2 + 2t4 ) (t1 + 2t3 ) (t2 − t4 ) (t1 − t3 )

(1.66)
L'expression totale de l'hamiltonien de Bloch basse énergie H(ky , kz ) =
H CNT (ky ) + HArr (kz ) est trouvée en sommant les expressions dans les équa-
tions (1.64) et (1.66). Pour le moment, cet hamiltonien de Bloch est une
matrice 4 × 4 entièrement occupée qui est dicile à analyser. Par une simple

33
transformation unitaire H(ky , kz ) = UH(ky , kz )U† avec
 
1 0 1 0
1 0 1 0 1
U= √  , (76)
2  i 0 −i 0 
0 i 0 −i
(1.67)
nous pouvons block-diagonaliser l'hamiltonien, donnant
 
H+ (ky , kz ) 02,2
H(ky , kz ) = , (77)
02,2 H− (ky , kz )
(1.68)
deux hamiltoniens block-diagonaux 2 × 2 que nous étiquetons + et − qui
prennent la forme
2t1 +t3 ′ e−iky T /2 + 2t2 +t4 cos(k L)
 
12 cos(k z L) t − t 12 z
H+ (ky , kz ) = .
t − t′ eiky T /2 + 2t212+t4 cos(kz L) 2t1 +t3
12 cos(kz L)
(78)
(1.69)

===== Page 29 =====

et

− t43 cos kz L t − t′ e−iky T /2 − t4

"

#
4 cos kz L
H − (ky , kz ) =
t4
t − t′ eiky T /2 − − t43



4 cos kz L cos kz L
(79)
respectivement.
La structure de bande basse énergie de ces deux hamiltoniens block-
diagonaux est tracée dans la Fig. 11. Nous notons que cette théorie basse
énergie décrit susamment les caractéristiques qualitatives de la partie basse
énergie d'une structure de bande complète obtenue par diagonalisation nu-
mérique d'un hamiltonien 48 × 48, Hnum (ky , kz ), dérivé dans l'annexe B et
tracé dans la Fig. 12. Nous notons que
les diérences entre les
modèles de
liaison forte basse énergie H (ky , kz ) et exact Hnum (ky , kz ) sont dues
aux interactions entre les modes basse énergie et haute énergie. Ces inter-
actions entraînent des changements quantitatifs subtils dans la structure de
bande basse énergie et pourraient être prises en compte via des techniques
perturbatives.
Les paramètres de saut inter-tube ont été choisis pour modéliser les ca-
ractéristiques qualitatives de la structure de bande électronique obtenue via
l'expérience DFT [20] : t1 = 0.25t, t2 = 0.35t1 , t3 = 0.65t1 et t4 = 0.35t3 .
Nous notons que les modèles de liaison forte présentés dans cette section
peuvent être ajustés pour correspondre plus étroitement aux résultats de
DFT.

34
Figure 12  Structure de bande énergétique basse énergie approximative
(E ) d'un réseau planaire de nanotubes de carbone quasi-métalliques (12, 0).
Les valeurs propres des hamiltoniens eectifs 2×2 H + (ky , kz ) et H − (ky , kz )
sont tracées en bleu et noir respectivement. Le panneau (a) esquisse la struc-
ture de bande électronique basse énergie le long de l'axe du vecteur d'onde du
réseau kz (pour ky = 0) tandis que le panneau (b) montre le vecteur d'onde
le long des nanotubes ky (pour kz = 0). La structure de bande héberge deux
paires de cônes de Dirac de type-II comme annoté. Dans les deux cas, les
vecteurs d'onde sont normalisés par la séparation des CNT L et le vecteur
unitaire des CNT T et l'énergie est tracée en unités d'électronvolts (eV).

===== Page 30 =====

expérience si nous réduisons articiellement le mini-gap ∆m d'un fac-
teur 210. Cette constatation peut être décrite par la possibilité d'une sous-
estimation du gap de bande dans certaines techniques DFT [83, 84]. Quelle
que soit la taille spécique du mini-gap, la structure de bande basse éner-
gie prédite par notre théorie conrme positivement les prédictions de quatre
croisements de bande super-critiquement inclinés dans la structure de bande
électronique. En utilisant cette théorie basse énergie, nous pouvons mainte-
nant développer les hamiltoniens au voisinage de ces croisements de bande
et prouver l'existence de cônes de Dirac de type-II.

Émergence des cônes de Dirac inclinés

Maintenant que nous avons divisé la structure de bande basse énergie en
deux hamiltoniens block-diagonaux 2 × 2 (H + (ky , kz ) et H − (ky , kz )), nous
pouvons analyser la structure de bande. Nous eectuons un développement
de Taylor sur l'hamiltonien H + (ky , kz ) autour des vecteurs d'onde ky = 0
et kz = s K+,z où
12 (t′ − t)
 
1
K+,z = arccos , (80)
L 2t2 + t4
et s = ±1 donne un hamiltonien de cône de Dirac incliné de la forme
h i
H + (q) = ℏv s γ η q z I + s η q z σx + q y σy . (81)

===== Page 31 =====

où q y = ky et q z = kz − sK+,z est l'écart du vecteur d'onde de chaque
point de Dirac, v = − ν2 T /2ℏ est la vitesse de Fermi eective, γ = (2t1 +
t3 )/(2t2 + t4 ) est le paramètre d'inclinaison et
s
(2t2 + t4 )L 144(t′ − t)2
η= 1 − , (82)
6t′ T (2t2 + t4 )2

35
Figure 13  Structure de bande basse énergie (E ) d'un réseau planaire de
nanotubes de carbone quasi-métalliques (12, 0) obtenue par diagonalisation
numérique de l'hamiltonien de Bloch complet du modèle de liaison forte
Hnum (ky , kz ). Le panneau (a) esquisse la structure de bande électronique
basse énergie le long de l'axe du vecteur d'onde du réseau kz (pour ky = 0)
tandis que le panneau (b) montre le vecteur d'onde le long des nanotubes
ky (pour kz = 0). Dans les deux cas, les vecteurs d'onde sont normalisés par
la séparation des CNT L et le vecteur unitaire des CNT T et l'énergie est
tracée en unités d'électronvolts (eV).

(1.73)
est le paramètre d'anisotropie (voir Fig. 1.11). Notez que l'énergie de dé-
calage dans l'équation (1.72) a été omise par souci de concision. Dans ce cas,
le paramètre s constitue un indice de vallée. En substituant les paramètres
de saut de réseau mentionnés précédemment t1 à t4 , nous voyons que le pa-
ramètre d'inclinaison est approximativement γ ≈ 2.86 ce qui signie que cet
hamiltonien décrit des cônes de Dirac inclinés de type-II (|γ| > 1). De plus,
comme le réseau planaire de nanotubes de carbone est une structure haute-
ment anisotrope (en raison de l'inadéquation de la taille des constantes de
réseau L et T ), nous obtenons un degré élevé d'anisotropie de la vitesse de
Fermi. Par exemple, si nous prenons la séparation inter-tube comme étant
quatre fois le rayon du tube L ≈ 4R = 18.8 et T = 4.26, alors le paramètre
d'anisotropie de la vitesse de Fermi est presque un ordre de grandeur plus
petit que l'unité (η ≈ 0.15). Nous notons que selon la forme de l'équation
(1.71), ces cônes de Dirac n'existent que si les paramètres de saut inter-tube
t2 et t4 et le mini-gap ∆m satisfont la condition 6∆m < |2t2 + t4 |. Le dé-
veloppement du second hamiltonien H− (ky , kz ) autour des vecteurs d'onde
ky = 0 et kz = sK−,z où

4(t − t′ )
 
1
K−,z = arccos , (83)
L t4
(1.74)
donne un second hamiltonien de cône de Dirac incliné qui décrit une
seconde paire de cônes de Dirac de type-II.

1.3.5 Conclusions
Dans cette section, nous avons développé un modèle analytique qui décrit
avec précision la structure de bande électronique basse énergie d'un réseau
planaire de nanotubes de carbone quasi-métalliques (12, 0). Curieusement,
un réseau de nanotubes de carbone avec gap (quasi-métalliques) forme une
structure de bande sans gap (métallique). Nos calculs démontrent que la
structure de bande contient deux paires de cônes de Dirac inclinés de type-II

36
au voisinage du niveau de Fermi - vériant les prédictions d'une étude DFT
précédente [20].
Bien que ce travail se soit concentré sur un réseau planaire de nanotubes
(12, 0), les mêmes méthodes peuvent être utilisées pour prouver l'existence
de cônes de Dirac inclinés dans la structure de bande électronique de réseaux
de nanotubes (6, 0) et (18, 0). Nous notons que nos résultats ne se généra-
lisent pas aux réseaux de CNT (9, 0) ou (15, 0), car ces nanotubes manquent
de symétrie de rotation de 180◦ et nécessitent donc un modèle analytique
séparé. Il a également été prédit que les réseaux de nanotubes (15, 0) pos-
séderont également une riche structure de bande électronique qui pourrait
potentiellement héberger des cônes de Dirac inclinés [20]. De plus, la même
étude a prédit que la formation d'un réseau à partir de nanotubes métal-
liques (15, 15) armchair peut entraîner l'ouverture d'un gap de bande. Ces
deux systèmes présentent des opportunités intéressantes

===== Page 32 =====

pour étudier théoriquement les réseaux planaires de nanotubes de car-
bone et constituent donc des voies possibles pour étendre ce travail. Il est
également noté que certains réseaux de nanotubes de carbone pourraient
subir une transition de Peierls et s'auto-assembler en un réseau avec deux
nanotubes par maille unitaire ; déterminer la structure relaxée d'un réseau de
nanotubes est probablement au-delà de la portée de la théorie et nécessiterait
donc des techniques computationnelles. La structure de bande électronique
d'un tel système pourrait être produite en combinant les méthodes discutées
dans cette section avec le modèle de Su-Schrieer-Heeger [85, 86, 87] (le mo-
dèle analytique utilisé pour produire la structure de bande électronique de
la molécule organique dimérisée polyacétylène).
Dans ce travail, nous avons montré que la structure de bande d'un ré-
seau planaire de nanotubes quasi-métalliques manque d'un gap de bande
électronique, limitant apparemment les applications térahertz. Cependant,
la présence de cônes de Dirac sans gap hautement anisotropes (par exemple,
anisotropie de la vitesse de Fermi η ≈ 0.15 pour L = 4R dans l'équation
(1.73)) permettra l'absorption sélective de photons térahertz polarisés le long
de l'axe du nanotube (l'absorption des matériaux à cônes de Dirac inclinés
est discutée au chapitre 2). La nature sans gap des cônes de Dirac de type-II
garantit qu'il n'y a pas de seuil de basse fréquence sur cette absorption, ce
qui signie que les réseaux de QMCNT pourraient trouver des applications
comme polariseurs térahertz à large bande en couche mince. Il est bien connu
que la taille du mini-gap dans un QMCNT peut être ajustée avec un champ
magnétique, ouvrant ainsi la possibilité de transitions optiques pour le rayon-
nement térahertz [88, 89, 90, 91, 92]. En appliquant un champ magnétique
dans le plan du réseau de QMCNT, il pourrait être possible de modier la
structure de bande basse énergie, fournissant une voie possible pour ajuster

37
la géométrie des cônes de Dirac inclinés et les propriétés optiques résultantes.

38