Phonon-Mediated Thermal Transport in Nanocrystalline Silicon Using Machine-Learning Interatomic Potentials

Technische Samenvatting: Fonoon-gemedieerd thermisch transport in nanocristallijn silicium met behulp van machine-learning interatomaire potentialen

Probleemstelling
Thermisch beheer in nanoschaalapparaten, met name die betrokken zijn bij nanocristallijn silicium, wordt fundamenteel beperkt door fonoonschattering bij korrelgrenzen (GB's) en interfaces. Hoewel klassieke interatomaire potentialen (zoals Stillinger–Weber en Tersoff) grootschalige simulaties mogelijk hebben gemaakt, schieten hun beperkte functionele vormen vaak tekort in het nauwkeurig vastleggen van vibratie-eigenschappen, anharmoniciteit en interface-wanorde. Dit leidt tot significante onnauwkeurigheden in het voorspellen van fonoondispersies, levensduur en thermische grensweerstand (TBR), vooral in structureel complexe omgevingen waar de lokale atomaire coördinatie afwijkt van de ideale bulk-symmetrie. Er is een kritieke behoefte aan interatomaire modellen die een nauwkeurigheid bieden die nabij die van first-principles, terwijl ze de computationele efficiëntie behouden die vereist is voor grootschalige, niet-evenwichtssimulaties.

Methodologie
De auteurs hebben een verenigd machine-learning (ML) framework ontwikkeld dat Gaussian Approximation Potentials (GAP) en equivariant neural-network-gebaseerde MACE-modellen integreert met roosterdynamische berekeningen en non-equilibrium molecular dynamics (NEMD).

1. Modeltraining en Validatie:

   • Een verenigde dataset van energieën en krachten werd gegenereerd uit bulk silicium configuraties (diamantkubische structuur) met behulp van Density Functional Theory (DFT) referentiedata.
   • Vier potentialen werden gebenchmarkt: een hertraind GAP-model, een MACE-model, en de empirische Stillinger–Weber (SW) en Tersoff potentialen.
   • Harmonische en anharmonische krachtconstanten werden berekend met behulp van de finite-displacement methode (Phonopy/Phono3py) op 3×3×3 supercells. Deze werden gebruikt om fonoondispersies, dichtheden van toestanden (DOS), levensduur en roosterthermische geleidbaarheid ($\kappa$) te berekenen.

2. Korrelgrens Modellering (Roosterdynamica):

   • Symmetrische kantelkorrelgrenzen (symmetric tilt grain boundaries) werden geconstrueerd met misoriëntiehoeken van 10°, 20°, 30° en 40°.
   • Het systeem werd gepartitioneerd in "nabij-GB" en "GB-kern" regio's om ruimtelijk opgeloste vibratie-eigenschappen te analyseren.
   • Een MACE-potentiaal, specifiek getraind op GB-configuraties, werd ingezet om fonoonlevensduur, gemiddelde vrije paden (MFP's) en groepsnelheden binnen de gedisorderde interface-regio's te berekenen.

3. Interfaciale Transport (NEMD):

   • ML-NEMD simulaties werden uitgevoerd in LAMMPS met de QUIP/QUIPY interface.
   • Bicristallijne systemen met gecontroleerde interface-ruwheid (sinusoïdale verplaatsingen met amplitudes $A = 1, 2, 3$ Å) werden onderworpen aan een temperatuurgradiënt ($T_{hot}=400$ K, $T_{cold}=300$ K).
   • De thermische grensweerstand (TBR) werd geëxtraheerd uit de steady-state temperatuurdiscontinuïteit ($\Delta T_j$) en de warmteflux ($J_x$).

Belangrijkste Resultaten

• Bulk Validatie: De MLIP's (GAP en MACE) reproduceerden accuraat experimentele fonoondispersies en DOS, waarbij de akoestische en optische takken correct werden gescheiden. In tegen tegenstelling hiertoe overschatten de SW en Tersoff potentialen systematisch de frequenties van de optische modi en onderschatten zij de anharmonische verstrooiing. Bij consequent waren de MLIP's in staat om fonoonlevensduur en roosterthermische geleidbaarheid ($\kappa$) in uitstekende overeenstemming met experimentele data en first-principles berekeningen te voorspellen, terwijl de klassieke potentialen $\kappa$ significant overschatten door kunstmatig lange fonoonlevensduur.
• Korrelgrens Vibratie-eigenschappen: In nanocristallijn silicium vertoonde de GB-kern significante spectrale modificaties, inclus\ de piekverbreding en spectrale gewichtsherverdeling, met name in laagfrequente akoestische modi. De fonoonlevensduur in de GB-kern werd drastisch gereduceerd tot de picoseconde-schaal over alle frequenties, wat wijst op sterke verstrooiing. Hoewel misoriëntiehoeken de mate van wanorde beïnvloedden, werd de primaire verstrooiingsmechanisme geïdentificeerd als lokale atomaire wanorde in plaats van de specifieke kantelhoek.
• Thermische Grensweerstand (TBR): NEMD-simulaties toonden aan dat TBR zeer gevoelig is voor zowel de interface-ruwheid als de keuze van de interatomaire potentiaal.
  • Afhankelijkheid van Ruwheid: Het vergroten van de ruwheidsamplitude ($A$) verhoogde systematisch de TBR voor alle potentialen.
  • Afhankelijkheid van Potentiaal: Het GAP-model voorspelde aanzienlijk hogere TBR-waarden (tot 1,66 m² K GW⁻¹ voor $A=3$ Å) vergeleken met SW (1,08) en Tersoff (1,17). Dit verschil komt voort uit het feit dat MLIP's sterkere anharmonische verstrooiing vastleggen, wat leidt tot grotere temperatuurspringen ($\Delta T_j$) en lagere effectieve warmtefluxen.
  • Experimentele Consistentie: De door het GAP-model voorspelde TBR-waarden voor ruwe interfaces komen goed overeen met experimentele bereiken (0,3–2,3 m² K GW⁻¹), terwijl klassieke potentialen de weerstand in sterk gedisorderde morfologieën de neiging hebben te onderschatten.

Betekenis en Claims
Het artikel claimt een verenigd, intern consistent framework te hebben gevestigd voor het modelleren van thermisch transport dat de kloof tussen nauwkeurigheid en schaal overbrugt. Door aan te tonen dat MLIP's een fysisch realistischerere beschrijving bieden van zowel harmonische als anharmonische roosterdynamica, betogen de auteurs dat voorspellende modellering van interfaciale warmteoverdracht in laagdimensionale materialen vereist dat men verder gaat dan empirische potentialen. De studie benadrukt dat de keuze van de interatomaire potentiaal de voorspellingen van de thermische grensweerstand aanzienlijk verandert, met name voor complexe, ruwe interfaces. Het werk valideert dat ML-afgeleide krachtconstanten betrouwbaar de onderdrukking van thermische geleidbaarheid in nanocristallijn silicium kunnen vastleggen, wat een robuust instrument biedt voor de thermische engineering van de volgende generatie nanomechanische en elektronische apparaten.