Inleiding tot 2D-materialen en warmtegeleiding

Tweedimensionale (2D) materialen, zoals grafeen, hexagonale boordioxyde (h-BN), en transitiemetaaldichhalcogeniden (TMD's) zoals molybdeendisulfide (MoS2), hebben een intense onderzoeksinteresse aangetrokken vanwege hun uitzonderlijke elektronische, optische en mechanische eigenschappen. Hun atomaire dikte . Vaak leidt slechts een of enkele atomen . . . tot quantuminvloeden die aanleiding geven tot nieuwe fysische fenomenen. Onder deze, thermische geleidbaarheid speelt een cruciale rol bij het bepalen van de prestaties en betrouwbaarheid van apparaten gebouwd uit deze materialen. Bijvoorbeeld, grafene vertoont een zeer hoge intrinsieke thermische geleidbaarheid (rond 3000.

De rol van de dichtheid functionele theorie

Stichtingen van DFT

Density Functional Theory is een quantum-mechanical modeling methode die het veel-elektron probleem in kaart brengt op een systeem van niet-interactieve elektronen die in een effectief potentieel bewegen. De kern van DFT is de Kohn .Sham benadering, die de totale energie van een systeem uitdrukt als een functioneel van de elektronendichtheid. Door het oplossen van de Kohn .Sham vergelijkingen zelf-consistent, onderzoekers verkrijgen de elektronische grondtoestand, waaronder de elektronendichtheid, totale energie, en krachten op atomen. Moderne DFT implementaties gebruiken uitwisseling . Correlatie functionele functies zoals de lokale dichtheid benadering (LDA) of de algemene gradiënt benadering (GGA, bv. PBE of PBEsol). Voor 2D materialen, van der Waal correcties (bijv. DFT‐D3) of vdW‐DF) zijn vaak noodzakelijk om de interlayer interacties nauwkeurig te beschrijven bij het modelleren van enkele lagen. DFT codes die veel voor dit doel worden gebruikt OJeasture: ESP[FLT], [FLT] en CFTP

Phononberekeningen met DFT

De temperatuurafhankelijke thermische geleidbaarheid van kristallen wordt gedomineerd door roostertrillingen of fononen. Om de fononeigenschappen van de eerste principes te berekenen, moet men de krachtconstanten van het kristal berekenen. Twee belangrijke methoden worden gebruikt: de eindige-verdrijvingsmethode en dichtheid functionele perturbatietheorie (DFPT)[. Bij de eindige-verdrijvingsaanpak worden kleine atomaire verschuivingen geïntroduceerd in een supercel en de resulterende krachten worden berekend uit DFT om harmonische en anharmonische IFC's te extraheren. DFPT, aan de andere kant, combineert de dynamische matrix door lineaire responstheorie, waardoor de behoefte aan grote supercellen voor lange-golflengte fonons wordt vermeden. Beide methoden vereisen zorgvuldige convergentie met betrekking tot k‐puntsteek, vlak-golfafsneiging en supercelgrootte. De fonondisperatieverhoudingen worden verkregen door diagonaliserende matrix, en de groepvelocities (v g) en fononononen worden afgeleid van de harmonischen van de thermische klanken van de derde geleiding

Methode voor Thermische geleidbaarheid Simulatie

Structuuroptimalisatie

De eerste stap in een thermische geleidbaarheidsstudie op basis van DFT is het verkrijgen van de volledig ontspannen atomaire structuur van het 2D materiaal. Voor monolagen wordt een vacuümgebied van minstens 15 Å loodrecht op het vlak toegevoegd om ongewenste interacties tussen periodieke beelden te elimineren. De parameters en atomaire posities in het vlak worden geoptimaliseerd totdat de resterende krachten onder een drempel liggen (gewoonlijk 10−4 eV/Å) en de spanningen zijn verwaarloosbaar. De keuze van de wisselwerking functioneel en van der Waals correctie kan significant invloed hebben op de ontspannen roosterconstanten en, bijgevolg, de fononfrequenties. Bijvoorbeeld, de PBE functionele neiging om de continue constanten van grafeen licht te overschatten, terwijl PBEsol een betere overeenkomst oplevert met experiment voor sommige TMDs. Een zorgvuldige convergentietest op het energieafsnijden en k‐puntrooster (bijv., 24×24×1 voor een eenheid grafeen) is essentieel om structurele stabiliteit en nauwkeurige krachten te garanderen.

Phonon Dispersion and Force Constanten

Zodra de evenwichtsstructuur is bepaald, worden harmonische fononberekeningen uitgevoerd. Met behulp van DFPT wordt de dynamische matrix berekend op een grof q-puntsrooster (bijv. 6×6×1) en vervolgens Fourier geïnterpoleerd op een fijnere mesh voor fonon dispersie plots. Akoestische somregels moeten worden afgedwongen om imaginaire modi (negatieve frequenties) die structurele instabiliteit zou aangeven te elimineren. Voor 2D materialen, de ZA (out-of-plane akoestische) tak heeft vaak een kwadratische dispersie in de buurt van het Γ punt, een handtekening van de verminderde dimensionaliteit. Voor anharmonische berekeningen, een supercel (bijv., 4×4×1 voor grafeen) wordt geconstrueerd, en een set van ontheemde configuraties wordt gegenereerd om derde-orde IFCs te berekenen via eindige verschillen. De cutoff radius voor anharmonische interacties wordt meestal gekozen om alle neighbors binnen een paar angstroms te omvatten; een onvoldoende cutoff kan leiden tot onnauwkeurige scatteringssnelheden, vooral voor materialen met lange-range interacties zoals grafeen.

Oplossen van de Boltzmann-vervoersvergelijking

De warmtegeleiding van het rooster κ lat wordt verkregen door het oplossen van de Boltzmann Transport Equation (BTE) voor de fononen. In de ontspanningstijd-apparatatie (RTA), κ = (1/N q V) Σ {q,ν} C {q,ν} v {q,ν}2 τ {q,ν}, waar C de modus-specifieke warmtecapaciteit is, v de groepssnelheid, en de levensduur van de fonon van alle verstrooiende mechanismen (anharmonisch, isotop, grens) te bepalen. Echter, de RTA negeert de correctie van de fononverdelingsfunctie van het evenwicht, die significant kan zijn voor materialen met sterke normale verstrooiing (bijv., grafeen). Een iteratieve oplossing van de BTE, zoals geïmplementeerd in codes zoals ShengBTE (die gebruik maakt van de volledige fononverdeling en derde-orde IFC's als input).

Praktische overwegingen

Verschillende numerieke parameters beïnvloeden de nauwkeurigheid van de berekende thermische geleidbaarheid:

  • Supercell size: Voor harmonische berekeningen is een voldoende grote supercel (of een dichte DFPT q‐mesh) nodig om verpakkingsfouten te voorkomen. Voor anharmonische krachtconstanten moet de supercel groter zijn dan de ruimtelijke omvang van de interatomaire interacties die in overweging worden genomen.
  • k-puntbemonstering: Een dicht k-puntmaas (bv. 12×12×1 voor een 2D-supercel) zorgt voor convergerende krachten en totale energie.
  • Spreukenschema: Gaussian of Methfessel
  • Isotopenverstrooiing: Natuurlijke isotopische mengsels (bv. 98,9% 12C en 1,1% 13C in koolstof) verminderen thermische geleidbaarheid met 10
  • Temperatuurafhankelijkheid: De levensduur van de fonon en de warmtecapaciteit zijn temperatuurafhankelijk, zodat κ lat wordt berekend bij meerdere temperaturen (bv. 100 K, 300 K, 500 K) om de volledige thermische geleidbaarheidscurve te verkrijgen.

Uitdagingen in DFT-gebaseerde thermische geleidbaarheid Simulaties

De servers zijn een combinatie van de printplaat en de printplaat.Ondanks het succes van de printplaat, kunnen de printplaten en de printplaat ook een printplaat vormen.De printplaat is een combinatie van de printplaat en de printplaat.De printplaat is een printplaat die de printplaat bevat.De printplaat is een printplaat die de printplaat bevat.De printplaat is een printplaat die de printplaat van de printplaat bevat.De printplaat is een printplaat die de printplaat bevat.De printplaat is een printplaat die de printplaat van de printplaat bevat.De printplaat is ook een printplaat die de printplaat bevat.De printplaat is een printplaat die de printplaat bevat.

Toepassingen en toekomstige aanwijzingen

Thermisch beheer in elektronica

De extreem hoge thermische geleidbaarheid van grafeen maakt het een natuurlijke kandidaat voor warmtespreiders in hoog vermogen elektronica. DFT simulaties helpen het gebruik van grafeen te optimaliseren door te voorspellen hoe substraatkoppeling, aantal lagen en korrelgrootte invloed hebben op warmtedissipatie. Zo tonen simulaties aan dat het schorsen van grafeen het substraat-geïnduceerde fononverstrooiing elimineert, waardoor de intrinsieke geleidbaarheid wordt hersteld, terwijl het op SiO2 plaatst, de geleidbaarheid vermindert tot 500 .1000 W/m·K door het verstrooien van de interface. Soortgelijke studies over hexagonale boriumferroceen en transition metaaldichhalcogeniden begeleiden het ontwerp van alle-2D heterostructuren waar thermisch beheer cruciaal is.

Thermo-elektrische materialen

Tweedimensionale materialen met een lage thermische geleidbaarheid zijn aantrekkelijk voor thermo-elektrische energieconversie, die een hoge Seebeck-coëfficiënt, hoge elektrische geleidbaarheid en lage thermische geleidbaarheid vereist. DFT-gebaseerde thermische geleidbaarheid screening heeft veelbelovende 2D-thermo-elektrische apparaten zoals SnSe (die een extreem lage κ lat van ~0,4 W/m·K in de richting van het vliegtuig), Bi2Te3 en weinig-laag WSe2 geïdentificeerd. Door de berekening van de fonon dispersie en verstrooiing rates, kunnen onderzoekers de mechanismen identificeren die leiden tot lage κ (bijvoorbeeld zachte akoestische modi, ratelende atomen, of complexe unitcellen) en chemische vervangingen voorstellen om thermische geleidbaarheid verder te verminderen zonder vernederende elektronische transport. De voorspellingen leiden experimentele synthese, versnellen de zoektocht naar hoog-efficiënte thermo-elektrische materialen.

Machine learning en high-throughput Screening

De hoge rekenkosten van BTE-berekeningen op basis van DFT hebben geleid tot de ontwikkeling van interatomaire mogelijkheden voor machinelearning (ML) die DFT-krachten getrouw kunnen reproduceren terwijl ze orden van grootte zijn. Modellen zoals het Behler.Parrinello-neurale netwerk, het Gaussiaanse approaching potential (GAP) en het diep potentieel (DP)[] zijn opgeleid op DFT-gegevens voor 2D-materialen en gebruikt om de eigenschappen van het fonontransport met bijna-DFT-nauwkeurigheid te berekenen. Deze ML-potentieelen maken thermische geleidbaarheidsberekeningen mogelijk voor grotere supercellen (bijv. 10 000 atomen) en langere tijdschalen, waardoor de studie van graangrenzen, nanotubes en verstoorde legeringen mogelijk is. Gekoppeld met hoge-doorlaatsframes zoals het Materials Project en AFLOWW, onderzoeken onderzoekers duizenden hypothetische 2D-verbindingen voor extreme thermische transporteigenschappen, die kunnen worden geïdentificeerd als borofene, silicene, silicene.

Conclusie

Door thermische geleidbaarheid in 2D-materialen te simuleren met de dichtheidsfunctionele theorie, kunnen onderzoekers een gedetailleerd, atomistisch begrip bieden van warmtetransportmechanismen die de prestaties van nanoschaalapparatuur ondersteunen. Door DFT-afgeleide krachtconstanten te combineren met de Boltzmann Transport Equation, kunnen onderzoekers intrinsieke thermische geleidbaarheidswaarden voorspellen, de invloed van defecten en belasting onderzoeken en de engineering van materialen voor thermisch beheer of thermo-elektrische conversie begeleiden. Ondanks de rekenuitdagingen, zullen vooruitgang in algoritmen, machine learning en high-performance computing de nauwkeurigheid en toegankelijkheid van deze simulaties gestaag verbeteren. Aangezien de bibliotheek van bekende 2D-materialen blijft uitbreiden, zullen op DFT-gebaseerde thermische geleidbaarheidsssimulaties een onmisbaar instrument blijven voor zowel wetenschappelijke ontdekking als technologische innovatie. De inzichten die uit deze simulaties worden verkregen, zijn al informatie over het ontwerp van de volgende generatie elektronica, energieconverters en flexibele sensoren, waarbij atomic-schaal begrip wordt omgezet in real-world toepassingen.