Atomistic substrate relaxation effects in the band gaps of graphene on hexagonal boron nitride

Dit onderzoek toont aan dat atomaire substraatrelaxatie een cruciale invloed heeft op de grootte en het teken van de bandgaten in graphene op hexagonaal boor-nitride, waarbij de relaxatie de primaire bandgap bij nul rotatiehoek aanzienlijk vergroot en een kleine, persistente bandgap over een breed hoekbereik mogelijk maakt.

De kern

Titel: De Dans van de Atomen: Waarom een Vloer van Stikstof en Boor de Eigenschappen van Grafen Verandert

Stel je voor dat je twee zeer dunne, bijna onzichtbare vellen papier op elkaar legt. Het bovenste vel is grafen (een net van koolstofatomen, zo sterk als staal maar dun als een haar). Het onderste vel is hexagonaal boor-nitride (h-BN), een soort "plaatje" van boor- en stikstofatomen dat als een perfect vlakke ondergrond dient.

Wanneer je deze twee vellen op elkaar legt, gebeurt er iets fascinerends: ze vormen een moirépatroon. Denk hierbij aan twee truien met een fijn patroon die je over elkaar trekt; je ziet dan een nieuw, groter patroon ontstaan dat eruitziet als een ruitjespatroon. In de wereld van atomen is dit patroon cruciaal voor hoe elektriciteit zich door het materiaal beweegt.

Het Grote Geheim: De "Vloer" is niet Stevig

In het verleden dachten wetenschappers dat de ondergrond (de h-BN) als een stijf, onbewegelijk betonnen vloertje fungeerde. Maar dit nieuwe onderzoek laat zien dat die "vloer" eigenlijk meer lijkt op een zacht matras.

Wanneer je de bovenste laag (grafen) een klein beetje draait ten opzichte van de onderlaag, gaan de atomen niet stilzitten. Ze glijden, duwen en trekken aan elkaar, net als mensen in een drukke menigte die proberen een comfortabele plek te vinden. Dit noemen we atomaire relaxatie.

Wat hebben de onderzoekers ontdekt?

1. De "Gaten" in de Elektronenstroom:
   Grafen heeft normaal gesproken geen "gat" in zijn energiestructuur (elektronen kunnen er vrij doorheen). Maar door het moirépatroon ontstaan er kleine gaten (bandgaps) waar elektronen niet kunnen komen.

   • De ontdekking: Als je de ondergrond als een stijf betonnen blok behandelt, zijn deze gaten heel klein. Maar als je rekening houdt met het feit dat de ondergrond meebuigt en meedraait (zoals een matras), worden deze gaten veel groter. Het is alsof je een deur dichtslaat: als de vloer meebeweegt, sluit de deur veel beter.

2. Het "Gouden Hoekje" van 0,6 graden:
   De onderzoekers draaiden de bovenste laag in verschillende hoeken. Ze ontdekten dat er een heel specifiek hoekje is (ongeveer 0,6 graden) waar de atomen het meest tevreden zijn.

   • De analogie: Stel je voor dat je twee tandwielen met een iets ander aantal tanden op elkaar zet. Op de meeste hoeken passen ze niet perfect. Maar op precies 0,6 graden vallen de "tanden" van het moirépatroon en de atomen van het grafen zo mooi in elkaar dat het systeem extra stabiel wordt. Op dit punt is de "deur" (het gat) het grootst.

3. De Rol van de Ondergrond:
   Als je de ondergrond volledig stijf houdt (alsof het aan de grond is gelijmd), verdwijnt dit speciale effect bijna. Maar als je toestaat dat de ondergrond meebeweegt, blijft er zelfs bij grote draaihoeken een klein, maar belangrijk gat over. Dit is cruciaal voor het bouwen van nieuwe, snelle elektronische apparaten.

Waarom is dit belangrijk voor ons?

Dit onderzoek is als het vinden van de perfecte instelling voor een radio. Als je de "knop" (de hoek en de flexibiliteit van de ondergrond) niet goed begrijpt, krijg je ruis of geen geluid.

• Voor de toekomst: Om superkrachtige, energiezuinige computers en sensoren te maken met grafen, moeten we precies weten hoe de atomen zich gedragen.
• De les: Je kunt niet zomaar aannemen dat een ondergrond stijf is. Soms moet je de "zachte matras" meenemen in je berekeningen, anders mis je de belangrijkste details.

Kort samengevat:
Deze wetenschappers hebben laten zien dat de manier waarop atomen in een laagje grafen en een laagje boor-nitride met elkaar "danssen" (zich verplaatsen en aanpassen), bepaalt of het materiaal elektriciteit goed of slecht geleidt. Door de "zachte" kant van de ondergrond mee te nemen in hun berekeningen, kunnen ze nu veel nauwkeuriger voorspellen hoe deze materialen zich gedragen, wat een enorme stap is voor de technologie van de toekomst.

Technische samenvatting

Titel: Atomaire substraat-relaxatie-effecten in de bandgaten van grafreen op hexagonaal boor-nitride (G/h-BN)

Auteurs: Jiaqi An, Nicolas Leconte, Srivani Javvaji, Youngju Park, en Jeil Jung (Universiteit van Seoel).

---

1. Het Probleem

Grafreen op hexagonaal boor-nitride (G/h-BN) vormt een prototypische van der Waals-heterostructuur waarbij de roosterverschil (lattice mismatch) tussen de twee materialen een moiré-potentiaal induceert. Dit leidt tot het openen van elektronische bandgaten bij het Dirac-punt (het primaire gat) en bij hogere energieën (het secundaire gat).

Er bestaat echter een aanzienlijke discrepantie tussen verschillende benaderingen voor het voorspellen van de grootte van deze gaten:

• Experimenten met tunnel-spectroscopie meten primaire gaten van ongeveer 14 meV.
• Continuümmodellen (gebaseerd op DFT-geïnspireerde benaderingen) voorspellen vaak lagere waarden rond 8 meV.
• Atomaire simulaties met volledige rooster-relaxatie voorspellen veel hogere waarden, variërend van 20 tot 35 meV.

De kernvraag is hoe de behandeling van substraat-relaxatie (de vervorming van het h-BN-ondersteunende laag) de grootte en het gedrag van deze bandgaten beïnvloedt, en waarom eerdere atomaire studies mogelijk onnauwkeurigheden vertonen door het substraat als te star te modelleren.

2. Methodologie

De auteurs hanteren een geavanceerde, hybride benadering die atomaire precisie combineert met efficiënte elektronische structuurberekeningen:

• Atomaire Relaxatie (MD): Ze gebruiken de LAMMPS-pakket met krachtvelden (DRIP voor interlaag-interacties, ExTeP en REBO2 voor intra-laag interacties) om de energie te minimaliseren. Ze onderzoeken vier specifieke configuraties om het effect van het substraat te isoleren:
  1. Volledig gerelaxeerd (Suspended): Beide lagen (grafreen en h-BN) zijn vrij om te relaxeren.
  2. Rigide contactlaag: De interface-laag van h-BN is star, maar een daarboven gelegen "remote" laag mag relaxeren.
  3. Rigide remote laag: Een extra h-BN-laag (op 180° gedraaid) wordt star gehouden om de relaxatie van de interface-laag te beperken.
  4. Volledig rigide: Geen enkele laag mag relaxeren.
• Elektronische Structuur (HTC Model): Ze gebruiken een Hybride Tight-Binding (HTC) model. Dit combineert:
  • Korte-afstand intra-laag termen gebaseerd op DFT-Wannier-functies (F2G2-model).
  • Inter-laag termen gebaseerd op een verfijnd twee-centrum (TC) model, gecorrigeerd met een exponentiële prefactor om tunneling bij verschillende inter-laag-afstanden nauwkeurig te beschrijven.
• Berekening van het Massaterm: Om de primaire bandgap efficiënt te berekenen over een groot bereik van draaihoeken, gebruiken ze de Truncated Atomic Plane Wave (TAPW) methode. Dit projecteert het volledige atomaire Hamiltoniaan op een gereduceerde subspace rond de Dirac-punten, waardoor de grootte van de matrix wordt verkleind zonder nauwkeurigheid te verliezen. De gemiddelde massaterm ($\bar{m}$) wordt direct afgeleid uit deze projectie.
• Supercellen: Ze construeren commensurate moiré-supercellen met vier gehele getallen ($p, q, p', q'$) om een breed scala aan draaihoeken ($\theta$) te bestrijken, variërend van 0° tot 60°.

3. Belangrijkste Bijdragen en Resultaten

A. Invloed van Substraat-Relaxatie op de Primaire Gap

De studie toont aan dat de behandeling van het substraat cruciaal is voor de voorspelde gap-grootte:

• Volledig gerelaxeerd (Suspended): De primaire gap is het grootst, beginnend bij ~30 meV bij 0° draaihoek.
• Rigide remote laag (Substraat-effect): Wanneer de remote h-BN-laag star wordt gehouden (wat een realistischere benadering van een substraat is), daalt de gap naar ~9 meV.
• Volledig rigide: In een volledig star systeem daalt de gap verder naar ~3 meV.
• Conclusie: De aanwezigheid van een substraat (die de relaxatie van de interface-laag beperkt) vermindert de primaire gap met een factor 2 tot 4 ten opzichte van een vrij hangend systeem. Dit verklaart de discrepantie tussen atomaire simulaties (vaak suspended) en experimenten.

B. Hoekafhankelijkheid en het "0.6°"-Minimum

• Bij een 0° draaihoek is de gap maximaal in gerelaxeerde systemen.
• Er wordt een niet-monotoon gedrag waargenomen: de primaire gap toont een lichte lokale piek (verhoging) rond 0.6°.
• Dit coincideert met een energetische stabilisatie waarbij de moiré-roostervectoren zich uitlijnen met de grafreen-roostervectoren. Deze uitlijning leidt tot een vergrote BA-stacking regio en een lokale versterking van de rooster-reconstructie.
• Bij grotere hoeken (>1°) neemt de primaire gap geleidelijk af tot waarden rond ~1 meV, maar sluit hij niet volledig. Dit staat in contrast met continuümmodellen die een volledige sluiting voorspellen bij grote hoeken.

C. Secundaire Gap

• De secundaire gap (bij hogere energieën) is minder gevoelig voor de specifieke relaxatie-scheme.
• Deze gap daalt van ongeveer 12 meV bij kleine hoeken naar 0 meV bij hoeken groter dan ~1°, ongeacht of het substraat star of flexibel is.

D. Rol van Lokale Rotatie en Stacking

• De auteurs kwantificeren de lokale rotatiehoek ($\theta_R$) en tonen aan dat bij kleine hoeken (rond 0.6°) de lokale rooster-reconstructie sterker is dan bij perfect uitgelijnde 0°.
• De grootte van de primaire gap correleert sterk met het oppervlak van de BA-stacking regio (de energetisch meest stabiele configuratie).
• Bij hoeken groter dan 20° wordt de gap gevoelig voor de keuze van het rotatiecentrum (AA, AB of BA stacking), wat leidt tot een "even-odd" onderscheid in de bandstructuur.

E. Massaterm en Sublattice-voorkeur

• De primaire gap wordt voor meer dan 70% bepaald door de gemiddelde massaterm ($\bar{m}$), die voortkomt uit het potentiaalverschil tussen de twee grafreen-sublattices veroorzaakt door het h-BN-substraat.
• Voor hoeken van 0° tot 30° behoudt de golf functie zijn voorkeur voor dezelfde koolstof-sublattice (geen tekenwisseling). Boven 30° wisselt de voorkeur van sublattice als gevolg van de trigonale symmetrie en de ongelijkwaardigheid van B en N atomen.

4. Betekenis en Conclusie

De studie biedt een cruciale link tussen atomaire simulaties, continuümmodellen en experimentele data:

1. Kwaliteit van Modellen: Het benadrukt dat het negeren van substraat-relaxatie (of het verkeerd modelleren ervan) leidt tot aanzienlijke fouten in de voorspelling van elektronische eigenschappen van moiré-materialen. Een "rigide remote layer" benadering lijkt een betere benadering voor experimentele setups dan een volledig suspended model.
2. Stabilisatie bij 0.6°: De ontdekking van een energetisch minimum en een lokale piek in de gap bij ~0.6° suggereert dat het mogelijk is om een niet-nul draaihoek te stabiliseren in heterobilayers, mits de vervorming van het interface-substraat onderdrukt kan worden.
3. Robuustheid: De primaire gap blijft zelfs bij grote hoeken (tot 30°) bestaan (ongeveer 1 meV), wat in strijd is met sommige theoretische voorspellingen van volledige sluiting, en wijst op de blijvende invloed van de substraat-potentiaal.

Samenvattend biedt dit werk een kwantitatief inzicht in hoe atomaire relaxatie en substraat-interacties de elektronische bandgaten van G/h-BN heterostructuren moduleren, en lost het de discrepantie op tussen eerdere theoretische voorspellingen en experimentele waarnemingen.