Simulating the influence of stoichiometry on the spectral emissivity of Mo$_x$Si$_y$ thin films

In dit werk wordt met behulp van dichtheidsfunctionale perturbatietheorie gesimuleerd hoe de stoichiometrie en kristalstructuur de spectrale emissiviteit van Mo$_x$Si$_y$-dunne films beïnvloeden, waarbij wordt aangetoond dat hexagonale MoSi$_2$ een lagere emissiviteit heeft dan de tetragonale fase en dat defecten de emissiviteit aanzienlijk kunnen verhogen.

De kern

De Kern: Een Warmte-Debat tussen Molybdeen en Silicium

Stel je voor dat je een heel dun laagje (slechts 20 nanometer dik, dat is 50.000 keer dunner dan een haar) van een speciaal metaal op een voorwerp spuit. Dit metaal is een mix van Molybdeen (Mo) en Silicium (Si).

De onderzoekers wilden weten: Hoe goed kan dit laagje warmte afgeven als straling?
In de wereld van warmte noemen we dit emissiviteit.

• Een emissiviteit van 1 is als een perfecte 'zwarte doos' die alle warmte afgeeft (zoals een gloeiende kolenhaard).
• Een emissiviteit van 0 is als een spiegel die de warmte terugkaatst en niets kwijtraakt.

Deze laagjes zijn belangrijk voor dingen die het heel heet moeten houden, zoals motoren of ruimteschepen. Als ze te heet worden, smelten ze. Als ze warmte goed kunnen afgeven (stralen), blijven ze koeler.

Het Grote Geheim: Het is niet alleen de 'receptuur'

De onderzoekers dachten eerst: "Als we meer Molybdeen toevoegen, wordt het laagje waarschijnlijk een betere warmte-afgever."
Maar hun computer-simulaties (die heel precies de atomen nabootsen) toonden iets verrassends: Het hangt niet alleen af van hoeveel Molybdeen erin zit.

Het is alsof je twee cakes bakt met exact hetzelfde ingrediëntenmengsel, maar je bakt ze op een heel andere manier. De ene cake is zacht en vochtig, de andere is hard en droog. Ze hebben dezelfde ingrediënten, maar een heel ander resultaat.

Bij deze metalen laagjes maakt de kristalstructuur (hoe de atomen precies in elkaar zitten) het verschil:

1. De vierkante structuur (Tetragonaal): Deze atomen zitten netjes in een vierkant patroon. Dit gedraagt zich als een goede warmte-afgever.
2. De zeshoekige structuur (Hexagonaal): Hier zitten de atomen in een zeshoekig patroon. Dit gedraagt zich als een slechte warmte-afgever (het houdt de warmte vast).

De Drie Manieren waarop Warmte 'Lekt'

De onderzoekers keken naar drie manieren waarop deze laagjes warmte kunnen uitstralen, en zagen hoe belangrijk ze zijn:

1. De trillende atomen (Ionen): Stel je voor dat de atomen als kleine balletjes aan veren trillen. Bij kamertemperatuur helpen deze trillingen een beetje, maar zodra het heet wordt (900 graden), stoppen ze met helpen. Ze zijn als een kleine rups die in de winter warmte geeft, maar in de zomer verdwijnt.
2. Elektronen die springen (Inter-band): Elektronen die van de ene plek naar de andere springen. Dit is de belangrijkste bron van warmte-afgifte bij de meeste materialen.
3. Elektronen die zwieren (Intra-band/Drude): Dit is als een zwerm bijen die in een pot ronddraait. Hoe sneller ze draaien, hoe meer warmte ze afgeven. Dit is heel belangrijk voor metalen.

De verrassing: Bij het materiaal dat het best warmte afgeeft (de vierkante structuur), is het juist deze 'zwerm bijen' die de warmte wegdraait. Bij het materiaal dat het slechtst afgeeft (de zeshoekige structuur), is er een kleine 'barrière' (een bandgap) die de bijen tegenhoudt. Ze kunnen niet vrij bewegen, dus ze stralen minder warmte uit.

De Magie van de Dikte

De dikte van het laagje is cruciaal.

• Is het laagje te dik? Dan werkt het alsof het een blok metaal is.
• Is het te dun? Dan werkt het niet goed.
• Het gouden punt: De onderzoekers ontdekten dat bij een dikte van 5 tot 10 nanometer de warmte-afgifte het hoogst is.

Dit is als een akoestische kamer. Als je te veel geluid hebt, klinkt het rommelig. Als je te weinig hebt, klinkt het dood. Maar op het juiste moment (de juiste dikte) kaatst het geluid (of in dit geval, de warmte) perfect heen en weer binnen het laagje, waardoor het veel efficiënter wordt afgegeven.

De 'Gekke' Defecten: Een Gebroken Spiegel

Tot slot keken ze naar wat er gebeurt als het kristal niet perfect is. Stel je voor dat je een perfect glazen raam hebt. Dat reflecteert heel goed. Maar als je er een paar krassen in maakt (defecten), verandert het gedrag.

In hun simulaties maakten ze 'krassen' in het kristalrooster (door atomen te verplaatsen of extra atomen toe te voegen).
Het resultaat? De warmte-afgifte steeg enorm!
De 'krassen' creëerden nieuwe plekken waar elektronen zich konden bevinden, waardoor ze makkelijker warmte konden uitstralen.

Conclusie voor de praktijk:
Als je een heel heet apparaat wilt maken dat niet smelt, wil je misschien niet de perfecte kristallen laagjes maken. Je wilt juist een beetje 'rommel' of onvolkomenheden in het materiaal, en je wilt dat het laagje heel dun is (rond de 5-10 nm). Dan stralen ze de warmte het beste af en blijven je machines koel.

Samenvattend in één zin:

De onderzoekers ontdekten dat je de warmte-afgifte van deze speciale metalen laagjes niet alleen kunt regelen door de ingrediënten te mengen, maar vooral door te spelen met de dikte van het laagje en door kleine onvolkomenheden in de structuur te creëren, net zoals je met een gebroken spiegel meer licht kunt verspreiden dan met een perfect gladde.

Technische samenvatting

Titel: Simulatie van de invloed van stoichiometrie op de spectrale emissiviteit van MoxSiy-dunne films

Auteurs: Zahra Golsanamlou et al. (Universiteit Twente en ASML Research)

---

1. Probleemstelling

Materialen die bestand zijn tegen hoge temperaturen zijn van toenemend belang voor toepassingen zoals turbinebladen, verbrandingskamers en EUV-lithografie (pellicles). Een kritieke uitdaging bij het gebruik van dunne films (ca. 20 nm) in deze omgeving is het thermisch beheer. Om oververhitting te voorkomen, moet er een thermisch evenwicht worden gevonden tussen verwarming en koeling. Passieve radiatieve koeling, bepaald door de emissiviteit van het materiaal, speelt hierbij een cruciale rol.

Hoewel experimenten hebben aangetoond dat de emissiviteit van molybdenum-silicium (MoxSiy) films sterk varieert afhankelijk van de kristalfase (bijv. hexagonaal vs. tetragonaal MoSi2) en de stoichiometrie, ontbreekt er een fundamenteel inzicht in hoe de atomaire structuur en samenstelling de diëlektrische eigenschappen en daarmee de emissiviteit beïnvloeden. Bestaande theoretische studies zijn beperkt en hebben vaak geen rekening gehouden met zowel elektronische als ionische bijdragen, noch met de specifieke effecten van dunne films (zoals meervoudige reflecties).

2. Methodologie

De auteurs hebben een eerste-principes (ab initio) modelleringsschema ontwikkeld om de emissiviteit van dunne MoxSiy-films te simuleren. De aanpak bestaat uit twee hoofdniveaus:

• Bulk-kristal berekeningen (DFT & DFPT):

  • Methode: Dichtetheorie-functie (DFT) met de Projector Augmented Wave (PAW) methode en het PBEsol-functie voor exchange-correlatie.
  • Systemen: Er zijn acht verschillende kristalfasen onderzocht (o.a. MoSi2, Mo5Si3, Mo3Si, MoSi3) met variërende stoichiometrieën ($0.25 \le x/(x+y) \le 0.75$).
  • Diëlektrische functie: De totale diëlektrische functie $\varepsilon(\omega)$ wordt berekend als de som van drie componenten:
    1. Inter-band elektronische bijdrage: Berekend met Dichtetheorie-stoornistheorie (DFPT).
    2. Intra-band elektronische bijdrage: Gemodelleerd met het Drude-model (voor metalen) met een plasmafrequentie ($\omega_p$) en dempingsfactor ($\gamma$).
    3. Ionische bijdrage: Berekend via harmonische fononmodi en Born-effectieve ladingstensors.
  • Defectstudies: Voor tetragonaal MoSi2 werden supercel-berekeningen uitgevoerd met defecten (Mo-overschot en Mo/Si-uitwisseling) om het effect van imperfecte kristalstructuren te testen.

• Dunne film simulatie (Transfer Matrix Method - TMM):

  • De berekende bulk-diëlektrische functies worden gebruikt als input voor de Transfer Matrix Method.
  • Dit model houdt rekening met meervoudige reflecties binnen de film, wat essentieel is omdat de optische doordringingsdiepte vaak groter is dan de filmdikte (ca. 20 nm).
  • De emissiviteit wordt afgeleid uit de absorptie ($A = 1 - T - R$) volgens de wet van Kirchhoff, waarbij rekening wordt gehouden met temperatuur (tot 3000 K) en invalshoek.

3. Belangrijkste Resultaten

• Metalliciteit en Bandstructuur:

  • De meeste onderzochte MoxSiy-fasen zijn metallisch met een hoge dichtheid van toestanden (DOS) bij het Fermi-niveau, behalve hexagonaal MoSi2, dat een kleine bandkloof (0,2 eV) heeft en zich als halfgeleider gedraagt.
  • Er is geen eenvoudige lineaire correlatie gevonden tussen de emissiviteit en het Mo-gehalte. De kristalstructuur en de specifieke elektronische bandstructuur zijn bepalender.

• Invloed van de Kristalfase (MoSi2):

  • Tetragonaal MoSi2 vertoont een aanzienlijk hogere emissiviteit dan hexagonaal MoSi2.
  • Reden: Tetragonaal MoSi2 heeft een gesloten bandkloof en een hoge DOS bij het Fermi-niveau, wat leidt tot sterke intra-band (Drude) en inter-band bijdragen. Hexagonaal MoSi2 heeft een bandkloof en een lage DOS, wat de emissiviteit onderdrukt.
  • De berekende emissiviteit voor een 20 nm tetragonale film bij 900 K is 0,37, wat goed overeenkomt met experimentele waarden (~0,33). Voor hexagonaal MoSi2 is de berekende waarde (0,06) lager dan experimenten (0,28), wat wordt toegeschreven aan defecten in de experimentele films.

• Invloed van Filmdikte:

  • Voor metallische films neemt de emissiviteit toe naarmate de dikte afneemt tot een kritieke waarde, door meervoudige reflecties.
  • De maximale emissiviteit wordt voorspeld bij een dikte van 5–10 nm. Dikkere films (bijv. 20-40 nm) hebben een lagere emissiviteit, terwijl zeer dunne films (<5 nm) weer afnemen door kwantumopsluitingseffecten (hoewel dit in de simulatie benaderd wordt).

• Invloed van Defecten:

  • De introductie van defecten (zoals atomaire uitwisseling of extra atomen) in tetragonaal MoSi2 leidt tot een drastische toename van de emissiviteit.
  • Defecten creëren lokale toestanden dicht bij het Fermi-niveau, wat de inter-band overgangen in het infrarood (IR) versterkt en de plasmafrequentie verlaagt. Dit suggereert dat films met een hoge emissiviteit vaak defecte of deels amorfe structuren hebben in plaats van perfecte kristallen.

• Temperatuurafhankelijkheid:

  • Bij kamertemperatuur is de bijdrage van ionische trillingen verwaarloosbaar.
  • Bij hogere temperaturen (900 K) domineren de elektronische bijdragen. De intra-band (Drude) bijdrage is cruciaal bij lagere temperaturen, terwijl inter-band overgangen bij hogere temperaturen belangrijker worden.

4. Bijdragen en Significance

• Fundamenteel Inzicht: Het werk koppelt voor het eerst de atomaire stoichiometrie en kristalfase van MoxSiy direct aan de spectrale emissiviteit van dunne films, inclusief de bijdrage van ionische trillingen en meervoudige reflecties.
• Validatie van Experimenten: De simulaties bevestigen de experimentele observatie dat de emissiviteit sterk afhankelijk is van de kristalfase (tetragonaal > hexagonaal) en verklaren de waargenomen waarden kwalitatief goed.
• Ontwerprichting voor Materialen:
  • Voor toepassingen die hoge emissiviteit vereisen (bijv. radiatieve koeling), is het niet noodzakelijk om perfecte kristallen te groeien; juist gecontroleerde defecten of een minder geordende structuur kunnen de emissiviteit maximaliseren.
  • De optimale filmdikte voor maximale emissiviteit ligt in het bereik van 5 tot 10 nm.
• Methodologische Vooruitgang: Het gekoppelde schema (DFT/DFPT + TMM) biedt een robuust kader om de optische eigenschappen van andere dunne filmmaterialen te voorspellen zonder afhankelijk te zijn van empirische fittingparameters.

Conclusie

De studie concludeert dat de emissiviteit van MoxSiy-dunne films niet triviaal te voorspellen is op basis van de chemische samenstelling alleen, maar sterk wordt bepaald door de kristalstructuur en de aanwezigheid van defecten. Tetragonaal MoSi2 met een dikte van 5-10 nm en een zekere mate van defecten biedt de beste prestaties voor hoge emissiviteit bij hoge temperaturen. Deze inzichten zijn waardevol voor het ontwikkelen van thermisch beschermende coatings en materialen voor extreme omstandigheden.