Een nauwkeurige DFT-1/2-benadering voor ondiepe defecttoestanden: efficiënte berekening van donorbindingsenergieën in silicium

Waarom kleine aanpassingen in silicium ertoe doen

Elke computerchip en zonnecel hangt af van zorgvuldig toegevoegde onzuiverheden, of “dopanten”, die bepalen hoe gemakkelijk elektriciteit stroomt. Voor sommige geavanceerde technologieën – van ultrazuinige transistoren tot donor-gebaseerde qubits – moeten we precies weten hoe sterk een extra elektron aan een dopantatoom in een halfgeleiderkristal gebonden is. Dit artikel introduceert een snellere en praktischer manier om die bindingsenergie met hoge nauwkeurigheid te berekenen, vooral voor dopanten in silicium, de ruggengraat van moderne elektronica.

Atomen die extra elektronen doneren

In zuiver silicium delen atomen elektronen in nette, herhalende bindingen, en het materiaal geleidt slecht bij kamertemperatuur. Voeg een klein beetje van een element uit groep V toe, zoals fosfor, arseen, antimoon of bismut, en elke dopant brengt één extra elektron mee. Dat extra elektron beweegt zich niet vrij; in plaats daarvan zit het in een waterstofachtige wolk die losjes aan de dopant en het omringende silicium gebonden is. De sterkte van deze binding – de donorbindingsenergie – bepaalt hoe gemakkelijk het elektron kan worden vrijgemaakt om stroom te voeren of deel te nemen aan kwantumoperaties. Het meten van deze energieën in het laboratorium is goed onderbouwd, maar ze betrouwbaar voorspellen vanuit eerst\-principesberekeningen blijkt moeilijk en kostbaar.

Waarom standaardberekeningen tekortschieten

Computermodellen gebaseerd op dichtheidsfunctionaaltheorie (DFT) zijn het werkpaard van materiaalkunde, maar ze hebben de neiging de lokalisatie van elektronen te onderschatten en de randen van de energiebanden in halfgeleiders verkeerd te plaatsen. Voor ondiepe donors, waarvan de elektronenwolken zich over vele tientallen atomen uitstrekken, betekent dit dat DFT meestal bindingsenergieën voorspelt die veel te klein zijn. Geavanceerdere methoden, zoals hybride functionalen en GW-berekeningen, kunnen deze problemen verhelpen maar tegen een enorme rekenprijs, vooral wanneer grote simulatiecellen met duizenden atomen nodig zijn om de uitgebreide donorstaat vast te leggen. Eerdere ‘tandem’-benaderingen moesten verschillende theoretische niveaus in verschillende celgroottes combineren en de resultaten aan elkaar kleven, wat de werkwijze complex en systeemafhankelijk maakte.

Een eenvoudige correctie met groot effect

De auteurs bouwen voort op een techniek die DFT-1/2 heet, waarbij een benaderde zelfenergiecorrectie direct in standaard DFT wordt ingebracht. In praktische termen passen zij het effectieve potentiaal van geselecteerde atomen licht aan door conceptueel een half elektron uit bepaalde atomaire orbitalen te verwijderen. Eerst passen ze deze correctie toe op bulk-silicium zodat de berekende bandopening veel beter met het experiment overeenkomt en een betrouwbaar referentiepunt voor de geleidingsband oplevert. Vervolgens onderzoeken ze het elektronische karakter van de donorstaat en vinden dat deze, voor alle groep‑V-dopanten, wordt gedomineerd door het s-orbitaal van het dopantatoom. Ze passen dan een op maat gemaakte half-elektron-correctie toe op dat orbitaal en stellen één afkastradius fijn af om de scheiding tussen het donorlevel en de dichtstbijzijnde lege geleidende toestand te maximaliseren. Belangrijk is dat deze geoptimaliseerde correctie geldig blijft wanneer de simulatiebox wordt vergroot, zodat ze hergebruikt kan worden in supercellen met tot duizenden atomen.

Hoe goed de methode werkt

Met deze tweestapscorrectie – eerst voor het gastmateriaal silicium, daarna voor de donor zelf – levert de methode donorbindingsenergieën die nauw aansluiten bij experimentele waarden. Voor arseen in silicium wijkt de voorspelde energie slechts 0,3 millielectronvolt af van het experiment, in wezen perfecte overeenstemming en vergelijkbaar met veel duurdere hybride berekeningen. Voor antimoon en fosfor zijn de fouten ongeveer 5 en 8 millielectronvolt, respectievelijk, een grote verbetering ten opzichte van ongecorrigeerde DFT. Voor bismut, een zeer zware dopant, includeren de auteurs ook spin–baankoppeling, een relativistisch effect dat energieniveaus licht vervormt. Dit brengt de berekende bindingsenergie binnen ongeveer 5 millielectronvolt van het experiment en belicht een fysica die eerdere, zwaardere methoden verwaarloosden. Om aan te tonen dat de benadering niet beperkt is tot silicium, passen ze dezelfde werkwijze met succes toe op een waterstofdonor in zinkoxide, waarbij opnieuw de gemeten bindingsenergieën binnen enkele millielectronvolt worden gereproduceerd.

Een praktisch hulpmiddel voor het ontwerpen van toekomstige chips

Voor niet‑specialisten is de kernboodschap dat de auteurs een recept hebben ontwikkeld dat de lage kosten en eenvoud van standaard DFT-berekeningen behoudt, terwijl het de nauwkeurigheid van veel zwaardere methoden bereikt. Door zowel de algemene bandstructuur van het gastmateriaal als de lokale omgeving rond de dopant systematisch te corrigeren, levert hun DFT-1/2-protocol betrouwbare donorbindingsenergieën in zeer grote simulatiecellen. Dit maakt het tot een krachtig en algemeen hulpmiddel voor het bestuderen van dopanten die dagelijkse elektronica en opkomende kwantumapparaten aansturen, en helpt ingenieurs materialen te ontwerpen waarin individuele onzuiverheden precies zo werken als bedoeld.

Bronvermelding: Claes, J., Partoens, B., Lamoen, D. et al. An accurate DFT-1/2 approach for shallow defect states: efficient calculation of donor binding energies in silicon. npj Comput Mater 12, 153 (2026). https://doi.org/10.1038/s41524-026-02003-2

Trefwoorden: ondiepe donors, silicium-dopanten, dichtheidsfunctionaaltheorie, kwantummaterialen, halfgeleiderdefecten