Density Profiles and Direct Correlation Functions from Density Functional Theory in Binary Hard-Sphere Crystals: Substitutional Solid and Interstitial Solid Solution

Dit artikel gebruikt klassieke dichtheidsfunctionaaltheorie om de evenwichtsdichtheidsprofielen en de inhomogene tweelichamsdirecte correlatiefuncties te bepalen voor twee typen binaire harde-sferen-kristallen, waarbij een fundamenteel verschil wordt aangetoond tussen de nauwelijks verstoorde vervangingsoplossingen en de sterk gedelokaliseerde interstitiële oplossingen.

De kern

Stel je voor dat je een enorme doos met balletjes hebt. Sommige balletjes zijn groot (zoals tennisballen) en sommige zijn klein (zoals knikkers). Als je deze doos schudt, gedragen ze zich als een vloeistof. Maar als je ze heel voorzichtig laat rusten, gaan ze vanzelf een strakke, kristallijne structuur vormen, net als de tegels op je vloer.

Deze wetenschappers hebben gekeken naar hoe deze balletjes zich gedragen in twee specifieke soorten kristallen, en hoe ze met elkaar "praten" (in de natuurkunde noemen we dit correlaties). Ze hebben dit niet gedaan door de balletjes fysiek te schudden, maar door een heel slim computerprogramma te gebruiken dat de wiskunde van de natuur nabootst.

Hier is wat ze hebben ontdekt, vertaald naar alledaags taal:

1. De Twee Soorten Kristallen

De onderzoekers keken naar twee manieren waarop de grote en kleine balletjes zich kunnen organiseren:

• Het "Verwisselbare" Kristal (Substitutional Crystal):
  Stel je een perfect geordend rijtje stoelen voor in een theater. Meestal zitten daar grote mensen. Maar soms zit er een klein persoon op een stoel die bedoeld was voor een grote persoon. Ze wisselen dus van plek. In dit geval gedragen de kleine balletjes zich heel veel als de grote: ze zitten strak op hun stoel en bewegen nauwelijks. Het is alsof je een rij tennisballen hebt, en af en toe vervang je er één door een knikker; de knikker zit net zo strak vast als de tennisbal.

• Het "Ingepakte" Kristal (Interstitial Solid Solution):
  Hier is de structuur anders. De stoelen (het rooster) zijn allemaal bezet door de grote tennisballen. Maar er zijn ook kleine ruimtes tussen de stoelen in (zoals de ruimte onder een stoel of in de gang). De kleine knikkers zitten hier niet op een vaste stoel, maar ze zwerven een beetje door die kleine ruimtes. Ze zijn niet vastgeklemd; ze zijn wat losser en kunnen makkelijker van de ene kleine ruimte naar de andere huppelen. Het is alsof de tennisballen een muur vormen en de knikkers als muisjes door de gaten in de muur rennen.

2. De "Lege Stoelen" (Vacatures)

In een perfect kristal zou elke stoel bezet moeten zijn. Maar in de echte wereld (en in hun berekeningen) zijn er altijd een paar lege stoelen.

• In het eerste type kristal zijn er heel weinig lege stoelen.
• In het tweede type kristal (waar de knikkers rondrennen) zijn er iets meer lege plekken voor de grote tennisballen.

De onderzoekers ontdekten iets fascinerends: hoe minder lege stoelen er zijn, hoe "angstiger" de balletjes worden om te bewegen. Als er bijna geen lege plekken zijn, is de kracht die ze op elkaar uitoefenen enorm groot. Het is alsof je in een volle trein staat: als er geen ruimte is om te bewegen, voelt elke kleine duw van een ander persoon als een enorme klap.

3. Het "Gesprek" tussen de Balletjes (Directe Correlatiefunctie)

Dit is het moeilijkste deel, maar we kunnen het vergelijken met een gesprek.
Stel je voor dat je een balletje op een bepaalde plek zet. Hoe beïnvloedt dat de kans dat er een ander balletje op een andere plek komt? Dit noemen ze de "directe correlatiefunctie".

• Bij het eerste type kristal: Omdat alles strak zit, is dit gesprek heel kort en krachtig. Als je een tennisbal verplaatst, merken de buren dat direct, maar het effect stopt snel.
• Bij het tweede type kristal: Hier is het gesprek heel anders.
  • De grote tennisballen praten nog steeds heel intens met elkaar (vanwege de lege stoelen).
  • Maar de kleine knikkers? Die praten heel anders. Omdat ze vrijer kunnen bewegen, is hun "gesprek" met de rest van het kristal veel minder strak. Ze gedragen zich bijna alsof ze in een vloeistof zitten, zelfs al zitten ze in een vast kristal. Ze kunnen makkelijk van de ene kant van de kamer naar de andere "huppelen" zonder dat de tennisballen er erg in hebben.

4. De "Tunnels" voor de Knikkers

De onderzoekers hebben ook gekeken naar hoe makkelijk de kleine knikkers kunnen bewegen. Ze ontdekten dat er tussen de vaste plekken van de knikkers (de octaëdrische gaten) en de andere plekken (de tetraëdrische gaten) een soort "tunnel" is.
Het kost de knikkers maar een klein beetje energie om door deze tunnel te gaan. Het is alsof je een heuvel moet beklimmen, maar de heuvel is zo laag dat je er zo overheen kunt springen. Dit betekent dat de kleine balletjes in dit kristal heel mobiel zijn; ze kunnen snel door het kristal "diffunderen" (verspreiden).

Samenvatting

Kortom:

1. Als je kleine en grote balletjes in een kristal stopt, kunnen ze op twee manieren zitten: of ze wisselen van plek (en gedragen zich allemaal strak), of de grote zitten vast en de kleine rennen vrij rond in de gaten.
2. In het kristal waar de kleine balletjes vrij rondrennen, gedragen ze zich bijna als een vloeistof, terwijl de grote balletjes nog steeds als een vast kristal zitten.
3. De "leegte" in het kristal (lege stoelen) bepaalt hoe sterk de balletjes op elkaar reageren. Hoe voller het is, hoe sterker de reactie.

Deze studie helpt ons beter te begrijpen hoe materialen zich gedragen op microscopisch niveau, wat belangrijk kan zijn voor het ontwikkelen van nieuwe materialen, zoals speciale verven, medicijnen of zelfs nieuwe soorten batterijen.

Technische samenvatting

Titel en Context

Titel: Dichtheidsprofielen en directe correlatiefuncties uit dichtheidsfunctionaalthorie in binaire harde-sferen-kristallen: Substitutionele vaste stof en interstitiële vaste oplossing.
Auteurs: Alessandro Simon en Martin Oettel (Universiteit Tübingen).
Kerngebied: Statistische mechanica, dichtheidsfunctionaalthorie (DFT), harde-sferen-systemen, kristallografie.

---

1. Het Probleem

Hard-sferen (HS) systemen zijn fundamenteel voor het begrijpen van de fasegedrag van colloïdale systemen en de structuur van vloeistoffen en vaste stoffen. Hoewel er uitgebreide analytische benaderingen (zoals de Percus-Yevick-oplossing) en numerieke simulaties bestaan voor vloeibare fasen, is de situatie voor kristallijne fasen complexer:

• Voor één-component kristallen is bekend dat de directe correlatiefunctie (DCF) sterk verschilt van die in de vloeibare fase en divergeert voor een ideaal (vacuüm-vrij) kristal.
• Er is echter weinig inzicht in de binaire (twee-componenten) harde-sferen kristallen, specifiek voor twee typen vaste oplossingen:
  1. Substitutionele kristallen (SC): Waar kleine sferen willekeurig grote sferen vervangen op roosterpunten.
  2. Interstitiële vaste oplossingen (ISS): Waar kleine sferen voornamelijk in de interstitiële holtes (octaëdrische gaten) van het rooster van de grote sferen zitten, zonder noodzakelijk aan een rooster gebonden te zijn.
• Er ontbreekt een volledig opgelost evenwichtsdichtheidsprofiel en de bijbehorende inhomogene twee-ligging directe correlatiefuncties ($c^{(2)}$) voor deze systemen, mede omdat de inhomogene Ornstein-Zernike-vergelijking computatieel zeer complex is om op te lossen via simulaties.

2. Methodologie

De auteurs gebruiken klassieke dichtheidsfunctionaalthorie (DFT) gebaseerd op de Fundamental Measure Theory (FMT), specifiek de White Bear II (WBII) functionaal.

• Model: Een binaire mengsel van harde sferen met diameters $\sigma_L$ (groot) en $\sigma_S = q\sigma_L$ (klein).
  • Substitutioneel (SC): $q=0.9$, samenstelling $x_S=0.3$.
  • Interstitieel (ISS): $q=0.3$, samenstelling $x_S=0.3$.
• Berekening:
  1. Het evenwichtsdichtheidsprofiel $\rho(r)$ wordt bepaald door het grand-potentiaal $\Omega[\rho]$ te minimaliseren.
  2. De minimalisatie gebeurt in twee stappen: eerst wordt het aantal deeltjes (en dus de vacuümconcentratie $n_{vac}$) vastgehouden om het dichtheidsprofiel te vinden, en vervolgens wordt geoptimaliseerd over $n_{vac}$.
  3. De ruimte wordt gediskretiseerd in een kubische eenheidscel ($128^3$ punten) met Gaussische pieken als startpunt.
  4. De directe correlatiefunctie $c^{(2)}_{ij}(r, r')$ wordt berekend via tweede-orde functionele differentiatie van het excess vrije-energie functioneel.
• Analyse: De auteurs analyseren de dichtheidsprofielen en de $c^{(2)}$-functies (een 6-dimensionale grootheid) langs kristallografische richtingen ([100], [110], [111]) en voor verschillende referentiepunten in de eenheidscel.

3. Belangrijkste Bijdragen en Resultaten

A. Dichtheidsprofielen en Vacuümconcentraties

• Substitutioneel Kristal (SC):
  • De dichtheidsprofielen voor zowel grote als kleine sferen lijken sterk op het één-component geval: smalle, bijna isotrope Gaussische pieken rond de fcc-roosterpunten.
  • De kleine sferen zijn iets minder gelokaliseerd door de grotere vrije ruimte rondom het roosterpunt.
  • De evenwichts vacuümconcentratie is zeer laag ($n_{vac} \approx 10^{-5}$), vergelijkbaar met het één-component kristal.
• Interstitiële Vaste Oplossing (ISS):
  • Grote sferen: Vormen een fcc-rooster met scherpe pieken, vergelijkbaar met het SC.
  • Kleine sferen: Vertonen een complexer, gedelokaliseerd gedrag. Ze bevinden zich niet strikt op roosterpunten, maar vullen de octaëdrische holtes en hebben een niet-verwaarloosbare dichtheid in de "bruggen" tussen deze holtes. Er is een secundaire piek bij de tetraëdrische holtes.
  • De vacuümconcentratie is hoger dan bij het SC ($n_{vac} \approx 3.7 \times 10^{-4}$), wat overeenkomt met simulaties.
  • De kleine sferen bezetten een klein percentage van de lege roosterpunten van de grote sferen, maar de meeste "leegte" betreft roosterpunten die niet bezet zijn door een groot deeltje.

B. Directe Correlatiefuncties (DCF)

De DCF is een matrix in de deeltjestypen ($LL, LS, SS$).

• Substitutioneel Kristal:
  • Alle componenten ($LL, LS, SS$) zijn zeer vergelijkbaar en lijken op het één-component geval.
  • De grootte is van de orde $O(1/n_{vac})$ (zeer groot in absolute waarde) met een kenmerkende dip bij $|r-r'| \approx \sigma/2$.
• Interstitiële Vaste Oplossing:
  • Groot-Groot ($LL$) component: Gedomineerd door de vacuümconcentratie. De grootte is $\sim 1/n_{vac}$. De functie kan worden geïnterpreteerd als een superpositie van kort-bereik, bijna isotrope basisfuncties gecentreerd op fcc-roosterpunten.
  • Klein-Klein ($SS$) en Groot-Klein ($LS$) componenten: Deze verschillen fundamenteel van het SC-geval. Ze tonen geen karakteristieke dip en zijn sterk anisotroop.
  • Vloeistof-achtig gedrag: Bij een kunstmatig hoge vacuümconcentratie ($n_{vac}=0.1$) blijken de correlaties die minstens één klein deeltje bevatten ($LS, SS$) perfect overeen te komen met de Percus-Yevick (PY) oplossing voor vloeistoffen. Dit suggereert dat de kleine sferen in de ISS-effectief in een vloeibare fase verkeren, ondanks de inhomogene dichtheidsverdeling. De grote sferen daarentegen vertonen kristalachtige correlaties.

C. Diffusiepaden

• De auteurs analyseren het energielandschap voor kleine sferen tussen octaëdrische en tetraëdrische holtes via de eerste-orde DCF ($c^{(1)}$).
• Er wordt een energiebarrière van ongeveer $2 k_B T$ gevonden voor het overschrijden van de octaëdrische naar de tetraëdrische holte.
• Dit impliceert dat kleine sferen in de ISS-fase zeer mobiel zijn en kunnen "hopen" tussen de interstitiële posities.

4. Significance en Conclusie

• Fundamenteel Inzicht: Het artikel levert de eerste volledig opgeloste, ruimtelijk gescheiden evenwichtsprofielen en DCF's voor binaire harde-sferen kristallen.
• Universeel Gedrag: Het bevestigt dat de divergentie van de DCF ($\sim 1/n_{vac}$) een universeel kenmerk is van kristallen met een roosterstructuur, maar dat de aard van de correlaties sterk afhangt van de lokale omgeving (roosterpunt vs. interstitiële holte).
• Toepassing: De resultaten zijn essentieel voor het berekenen van veralgemeende elastische constanten voor binaire kristallen, analoog aan eerdere werken voor één-component systemen.
• Fysisch beeld: Het biedt een geometrisch beeld van de 6-dimensionale DCF in kristallen, waarbij de correlaties worden verklaard door de waarschijnlijkheid van het invoegen van deeltjes in vacuümplaatsen en de superpositie van basisfuncties rond roosterpunten.

Samenvattend toont dit werk aan dat DFT met de White Bear II functionaal een krachtig instrument is om de complexe micro-structuur en correlaties in binaire kristallijne systemen te ontrafelen, waarbij een duidelijk onderscheid wordt gemaakt tussen het gedrag van substitutieel en interstitieel deeltjes.