Ajustement de la fonction de travail du graphène par la taille, la fonctionnalisation, les défauts et le dopage : une étude par principes premiers

Résumé

Nous utilisons la théorie de la fonctionnelle de la densité (DFT) aux principes premiers pour concevoir et évaluer systématiquement la fonction de travail (WF) du graphène pour applications en dispositifs électroniques. En explorant des structures de graphène pur, hydroxyle-modifié, défectueux et dopé, nous révélons que la WF dépend fortement de la largeur du ruban, de la fonctionnalisation de surface, de la position des défauts et du type de dopant. Le graphène zigzag pur présente une WF plus élevée que le graphène armchair, et l'augmentation de la largeur du ruban abaisse constamment la WF. Les groupes hydroxyle augmentent la WF, avec un effet amplifié par le nombre de groupes et leur regroupement. Les défauts réduisent généralement la WF, les lacunes centrales ayant l'impact le plus prononcé. Le dopage substitutionnel par B ou Al (type p) élève significativement la WF, tandis que N ou P (type n) la diminue fortement ; le dopage Si a un effet modeste et stable. Nous calculons aussi les énergies de formation pour évaluer la faisabilité des dopants, montrant que B et N s'incorporent aisément, alors que Al est difficile. Ces insights fournissent une base théorique robuste pour adapter les propriétés électroniques du graphène en ingénierie de dispositifs.

Contexte

Les propriétés électriques, mécaniques et chimiques exceptionnelles du graphène l'ont propulsé dans des applications diverses — capteurs, transistors à effet de champ (FET), électrodes photovoltaïques, diodes Schottky, tubes sous vide et jonctions LED — surpassant souvent les matériaux conventionnels. Cependant, la performance des dispositifs repose sur la fonction de travail (WF) du graphène, qui gouverne l'injection de charge, la hauteur de barrière Schottky et l'alignement de bande global. Les valeurs expérimentales typiques de WF varient de 4,2 à 4,8 eV. De nombreuses stratégies ont été rapportées pour ajuster la WF, incluant le dopage chimique, la fonctionnalisation de surface, l'ingénierie de défauts, la contrainte, le gating électrostatique et l'irradiation laser. Pourtant, les études comparatives à travers la chiralité, la taille et la position des défauts du graphène restent limitées, et la faisabilité thermodynamique des dopants est souvent négligée. Ce travail comble ces lacunes par une analyse exhaustive aux principes premiers.

Méthodes

Tous les calculs ont été effectués avec le code CASTEP utilisant la théorie de la fonctionnelle de la densité (DFT). Nous avons employé l'approximation de gradient généralisé (GGA) pour l'échange-corrélation, choisie pour sa gestion supérieure des densités électroniques inhomogènes. L'espacement de vide a été fixé à 15 Å pour éliminer les interactions entre couches, et des pseudopotentiels ultrasoft avec une coupure à 340 eV ont été utilisés. Un maillage k de Monkhorst-Pack 9×9×1 et un lissage Methfessel-Paxton (0,05 eV) ont assuré la convergence. Les fonctions de travail ont été obtenues à partir de la différence entre le niveau de vide et l'énergie de Fermi (WF = V₀ – E_f).

Résultats et discussion

Effet de la largeur du ruban et de la chiralité sur la WF

Nous avons modélisé des nanorubans de graphène (GNR) zigzag et armchair de 1 à 7 cellules unitaires en largeur. La WF diminue monotone avec l'augmentation de la largeur pour les deux chiralités ; les rubans plus larges approchent la WF du graphène massif (~4,48 eV). Les GNR zigzag affichent constamment une WF plus élevée que les homologues armchair, attribuable à leurs symétries cristallines distinctes. Les analyses de gap de bande (tableau 1) confirment que des largeurs plus grandes réduisent les gaps, soutenant les tendances de WF.

Influence de l'hydroxylation, des défauts et de leurs positions

La fonctionnalisation hydroxyle augmente la WF, un seul OH la portant à 4,50 eV et plusieurs groupes OH atteignant 5,10 eV. Des arrangements OH symétriques et groupés produisent la WF la plus élevée, tandis que des groupes dispersés la réduisent légèrement (4,66–4,83 eV). Les défauts, modélisés comme lacunes de carbone uniques, abaissent généralement la WF. Les lacunes centrales donnent la WF la plus basse (4,34 eV), alors que les lacunes décentrées augmentent légèrement la WF par rapport aux GNR purs (4,36 eV). Ces résultats s'alignent avec des observations expérimentales antérieures de réduction de WF induite par défauts.

Effets de dopage par B, N, Al, Si et P

Les concentrations de dopage ont varié de 2,4 % à 14,6 % (1–6 atomes dopants dans une supercellule 4×4). Les dopants type p B et Al augmentent la WF, B atteignant 5,15 eV à 14,6 % de dopage. Les dopants type n N et P diminuent la WF dramatiquement, descendant à 3,23 eV pour P à la concentration la plus élevée. Le dopage Si montre une WF modeste et stable. Les calculs d'énergie de formation révèlent que les dopants B et N sont thermodynamiquement favorables, tandis que Al, Si et P présentent des énergies de formation élevées, indiquant une moindre probabilité d'incorporation.

Conclusions

Notre étude systématique aux principes premiers démontre que la fonction de travail du graphène peut être finement ajustée par la largeur du ruban, la fonctionnalisation hydroxyle, l'ingénierie de défauts et le dopage sélectif. Ces résultats fournissent des directives claires pour optimiser les propriétés électroniques du graphène dans les architectures de dispositifs de prochaine génération.