Diverse polymorphism in Ruddlesden-Popper chalcogenides

In deze studie wordt een machine-learned interatomair potentiaal gebruikt om de structurele diversiteit van Ruddlesden-Popper chalcogeniden te onderzoeken, waarbij nieuwe polymorfen worden voorspeld en unieke thermische gedragingen zoals negatieve thermische expansie en ongewone symmetriebreking worden geïdentificeerd.

De kern

De dans van de atomen: Een verhaal over mysterieuze kristallen en hun verborgen vormen

Stel je voor dat je een gigantische, driedimensionale legpuzzel hebt. De stukjes van deze puzzel zijn atomen, en ze vormen kristallen. In de wereld van de materialenwetenschap zijn er bepaalde kristallen die bekend staan om hun veelzijdigheid: de Ruddlesden-Popper (RP) chalcogeniden.

Voorheen wisten wetenschappers al veel over een soort kristallen die "oxide" heten (waar zuurstof in zit). Die konden ze als een Lego-set gebruiken om supergeleiders of speciale materialen te bouwen die niet uitzetten als ze warm worden. Maar er was een hele nieuwe, nog onontdekte versie van deze Lego-set: de chalcogeniden (waar zwavel in zit). Deze zijn net als de oude versie, maar dan met een andere "kleefstof" (zwavel in plaats van zuurstof), wat ze veelbelovend maakt voor zonnepanelen en sensoren. Het probleem? Niemand wist precies hoe ze zich gedroegen als ze warm werden of als je de puzzelstukjes anders stapelde.

De digitale voorspeller
Om dit mysterie op te lossen, hebben de auteurs van dit paper een slimme, digitale "profeet" gecreëerd. Ze hebben een machine learning-model (een soort super-slimme computer) getraind met duizenden berekeningen. Je kunt dit vergelijken met het trainen van een hond: je laat hem duizenden keren zien hoe atomen zich gedragen, zodat hij uiteindelijk zelf kan voorspellen wat er gebeurt zonder dat je elke keer de zware rekenkracht van een supercomputer nodig hebt.

Met deze digitale voorspeller hebben ze een enorme simulatie uitgevoerd. Ze keken naar een familie van kristallen, genaamd Ba(n+1)ZrnS(3n+1). De 'n' in deze naam is het belangrijkste: het vertelt je hoeveel lagen van het kristal er in één blok zitten.

• n=1: Een dunne laag (één blok).
• n=3: Drie lagen dik.
• n=6: Zes lagen dik.
• n=oneindig: Een heel groot blok (het oorspronkelijke kristal).

Wat hebben ze ontdekt? De verrassingen

1. De "Krimpende" Kristallen (Negatieve thermische expansie)
   Normaal gesproken gebeurt dit: als je een ballon opblaast of een metaal verwarmt, zet het uit. De atomen trillen harder en duwen elkaar uit elkaar.
   Maar bij de dunste kristallen (n=1) gebeurde er iets raars: ze krompen. Als je ze verwarmde, werden ze smaller in het midden.
   De analogie: Denk aan een elastiekje dat je verwarmt. In plaats van losser te worden, trekt het zich samen alsof het een koude douche krijgt. Dit is extreem zeldzaam en heel nuttig voor materialen die niet mogen vervormen bij temperatuurwisselingen.

2. De "Opeenvolgende" Dans (Symmetrie-breken)
   Meestal geldt: hoe warmer het wordt, hoe chaotischer en "vrijer" de atomen bewegen. Ze verliezen hun strakke structuur en gaan naar een hogere, simpele vorm (zoals een kubus).
   Maar bij deze kristallen gebeurde het tegenovergestelde. Bij het verwarmen werden de atomen strakker en complexer in hun dans. Ze braken hun eigen symmetrie op een manier die niemand had verwacht.
   De analogie: Stel je een groep dansers voor die in een perfecte cirkel staan. Normaal gesproken, als de muziek sneller gaat (warmte), gaan ze wild rondrennen en valt de cirkel uit elkaar. Maar hier, als de muziek sneller gaat, gaan ze ineens in een nog ingewikkelder, strakker formatie dansen. Het is alsof de hitte hen dwingt om zich nog meer te concentreren.

3. De "Rimpeling" en de "Til"
   Het geheim van deze gedragingen zit in twee bewegingen:

   • De Til (Tilting): De atoomblokken (oktaëders) kantelen, alsof ze op hun kop worden gezet.
   • De Rimpeling (Rumpling): De lagen van atomen gaan niet meer perfect vlak liggen, maar worden als een golfbeweging op en neer geduwd.
   • De interactie: De computer zag dat deze twee bewegingen met elkaar "vechten". Als de atomen te veel gaan kantelen, worden ze gedwongen om te rimpelen, en andersom. Bij de dunne lagen (n=1) gebeurt dit aan de randen van het kristal op een manier die de binnenste lagen beïnvloedt. Bij de dikkere lagen (n=3 en hoger) gedraagt het kristal zich meer als een groot blok, maar met een speciale "golf" aan de randen.

Waarom is dit belangrijk?
Dit onderzoek is als het vinden van de gebruiksaanwijzing voor een nieuwe, krachtige machine. Door te begrijpen hoe deze kristallen zich gedragen (waarom ze krimpen, waarom ze hun vorm veranderen), kunnen ingenieurs in de toekomst:

• Zonnepanelen maken die efficiënter zijn.
• Materialen ontwerpen die niet uitrekken of krimpen bij temperatuurveranderingen (cruciaal voor precisie-instrumenten).
• Nieuwe elektronica bouwen die sneller en stabieler werkt.

Kortom:
De auteurs hebben met een slimme computer ontdekt dat deze nieuwe kristalfamilie veel meer "personality" heeft dan gedacht. Ze zijn niet statisch; ze dansen, krimpen en veranderen van vorm op verrassende manieren als je ze verwarmt. Door deze dans te begrijpen, kunnen we in de toekomst betere technologieën bouwen. Het is alsof we net de taal hebben geleerd die deze atomen spreken, en ze vertellen ons nu hoe we ze kunnen gebruiken.

Technische samenvatting

Probleemstelling

Ruddlesden-Popper (RP) chalcogeniden (met de formule $Ba_{n+1}Zr_nS_{3n+1}$) vormen een opkomende klasse van gelaagde halfgeleiders met veelbelovende eigenschappen voor optoelektronica en thermoelektrica. Hoewel de structurele diversiteit van RP-oxiden de afgelopen vier decennia uitgebreid is bestudeerd om geavanceerde functionaliteiten te realiseren, blijft het structurele gedrag van RP-chalcogeniden slecht begrepen.
De belangrijkste uitdagingen zijn:

1. Het identificeren van toegankelijke polymorfen is moeilijk omdat de structurele verschillen vaak subtiel zijn en experimenteel of theoretisch moeilijk op te lossen.
2. Er is een kritieke kennislacune over de structurele evolutie en fase-overgangen van deze materialen bij variërende temperatuur en waarden van $n$ (het aantal perovskietlagen).
3. Bestaande experimentele data zijn beperkt en er is behoefte aan een nauwkeurige voorspelling van de grondtoestanden en thermische gedragingen.

Methodologie

De auteurs hebben een geavanceerde computationele aanpak ontwikkeld om deze lacune op te vullen:

• Machine-Learned Interatomic Potentiaal (NEP): Er is een hoog-nauwkeurige "Neuroevolution Potential" (NEP) ontwikkeld. Dit model is getraind op data afkomstig van Dichtheidsfunctionaaltheorie (DFT) berekeningen (HSE06 functional) voor 1375 unieke structuren, inclusief de $Ba_{n+1}Zr_nS_{3n+1}$ serie ($n=1$ tot $6$) en de parent perovskiet $BaZrS_3$. Het model bereikt een foutmarge van slechts 1,8 meV/atoom voor vormingsenergieën.
• Grootschalige Moleculaire Dynamica (MD): Met het getrainde NEP-model werden uitgebreide MD-simulaties uitgevoerd (NPT ensemble) bij temperaturen variërend van 0 K tot 1200 K. Dit maakte het mogelijk om de thermische evolutie van de kristalstructuren en fase-overgangen te bestuderen op een schaal die met pure DFT onhaalbaar zou zijn.
• Validatie: De resultaten werden gevalideerd door vergelijking met bestaande experimentele röntgendiffractie (XRD) data en door het analyseren van fononenspectra.

Belangrijkste Bijdragen en Resultaten

De studie onthult een rijk polymorf landschap en onthult nieuwe structurele fenomenen die afhankelijk zijn van de waarde van $n$:

1. Voorspelling van Grondtoestanden en Polymorfen:

   • Voor elke $n$-waarde (1 tot 6) zijn nieuwe grondtoestanden geïdentificeerd met specifieke octaëdrische kantelpatronen (tilt patterns).
   • Lage-$n$ fasen ($n=1-3$): Kenmerken zich door een enkel uit-fase kantelpatroon ($\phi$) langs één as.
   • Hoge-$n$ fasen ($n \geq 4$): Vertonen versterkte kanteling in beide slabs, met uit-fase kanteling langs de x- en y-as en in-fase kanteling ($\psi$) langs de z-as. Voor $n \geq 4$ convergeren de grondtoestanden naar die van de parent perovskiet ($Pnma$).

2. Temperatuurafhankelijke Fase-overgangen:

   • De materialen ondergaan complexe overgangen bij verhitting. Hoge-$n$ verbindingen vertonen twee overgangen: een eerste-orde overgang bij lage temperatuur en een tweede-orde overgang bij hoge temperatuur.
   • Bij $n=1$ ($Ba_2ZrS_4$) wordt negatieve thermische uitzetting in het vlak waargenomen tot 420 K, met een uitzettingscoëfficiënt van $\alpha \approx -3 \times 10^{-6} K^{-1}$.
   • Een opvallend fenomeen is de "oplopende symmetriebreking": Bij $Ba_2ZrS_4$ en $Ba_4Zr_3S_{10}$ leidt verhitting tot een verlaging van de symmetrie (naar $Cmca$ respectievelijk $P2_1/c$), wat ongebruikelijk is aangezien hogere symmetrie normaal gesproken bij hogere temperaturen wordt gestabiliseerd.

3. Laag-afhankelijke Kantelpatronen en Oppervlakte-overgangen:

   • Voor $n \geq 3$ worden nieuwe, laag-afhankelijke kantelpatronen ontdekt die nog niet eerder in anorganische RP-verbindingen zijn waargenomen.
   • In $Ba_4Zr_3S_{10}$ ($n=3$) treedt een overgang op waarbij de kanteling in de centrale laag van de slab verschilt van de buitenste lagen.
   • Bij hoge-$n$ materialen ($n \geq 4$) worden "oppervlakte-overgangen" waargenomen: de kanteling aan de perovskiet-rotsalt interface wordt onderdrukt terwijl de bulk-achtige lagen gekanteld blijven.

4. Mechanisme: Interplay tussen Kanteling en "Rumpling":

   • De auteurs identificeren de onderliggende oorzaak van dit gedrag als de competitie tussen octaëdrische rotaties en rumpling (een uit-plaatsing van de A-site kationen, hier Ba, in de rotsaltlagen).
   • Er is een inverse koppeling: wanneer de uit-plaatsing van de octaëders aan de interface afneemt, neemt de rumpling-amplitude toe. Dit mechanisme verklaart waarom bepaalde symmetrie-verlagende overgangen optreden en waarom de thermische uitzetting negatief kan zijn.

Significantie

Deze studie heeft een aanzienlijke impact op het veld van de materiaalkunde:

• Methodologische Doorbraak: Het succesvol toepassen van machine-learned potentiaalmodellen op RP-chalcogeniden bewijst dat complexe, temperatuur-afhankelijke structurele evoluties nauwkeurig kunnen worden voorspeld, wat een alternatief biedt voor kostbare en tijdrovende experimenten.
• Fundamenteel Begrip: Het onthult dat de interface tussen perovskiet- en rotsaltlagen een cruciale rol speelt in het bepalen van de structurele eigenschappen, wat een nieuw perspectief biedt op het ontwerp van gelaagde materialen.
• Toepassingsgericht: De ontdekking van negatieve thermische uitzetting en de mogelijkheid om eigenschappen te "tunen" via de $n$-waarde en interface-engineering (bijv. via selectieve doping) opent nieuwe wegen voor het ontwikkelen van materialen met specifieke thermische en functionele eigenschappen voor optoelektronische en thermoelektrische toepassingen.

Kortom, dit werk legt de basis voor een rationeel ontwerp van geavanceerde RP-chalcogeniden door het inzichtelijk maken van hun complexe polymorfe landschap en de onderliggende fysica van hun structurele stabiliteit.