Random close packing fraction of bidisperse discs: Theoretical derivation and exact bounds

Dit artikel presenteert een theoretische afleiding en exacte bovengrenzen voor de maximale dichtheid van willekeurige, ongeordende pakkingen van bidisperse schijven, waarbij de dichtheid wordt uitgedrukt als een functie van de grootteverhouding en concentratie.

De kern

De Dichtste Stapel: Een Reis door de Wiskunde van Willekeurige Schijven

Stel je voor dat je een grote doos hebt vol met ronde koekjes. Je wilt ze zo strak mogelijk in de doos proppen, maar er is één belangrijke regel: je mag ze niet in een perfect, saai patroon leggen (zoals een ruitjespatroon). Ze moeten eruitzien alsof ze willekeurig zijn neergegooid, maar toch zo dicht mogelijk bij elkaar zitten.

Dit klinkt als een simpel puzzeltje, maar voor wiskundigen en natuurkundigen is het al decennia lang een enorme uitdaging. Hoe dicht kun je ze eigenlijk stapelen zonder dat ze "vastlopen" in een geordend kristal?

In dit nieuwe artikel legt de auteur, Raphael Blumenfeld, uit hoe hij dit probleem oplost voor twee soorten koekjes: kleine en grote. Hij gebruikt een slimme wiskundige truc om de absolute grens te vinden van hoe vol zo'n doos kan worden.

Hier is de uitleg, vertaald naar alledaags taal:

1. Het Probleem: De "Willekeurige" Stapel

Als je alleen maar één soort koekje hebt (allemaal even groot), proberen ze van nature om in een perfect hexagonaal patroon te vallen (zoals een honingraat). Dat is de dichtste manier, maar het is niet "willekeurig" meer. Om een willekeurige, rommelige stapel te krijgen, gebruiken mensen vaak een mix van kleine en grote koekjes. De grote vullen de gaten, de kleine vullen de rest.

Maar hoe weet je of je stapel echt willekeurig is? En wat is de theoretische limiet? Als je te veel probeert te persen, krijg je per ongeluk toch een geordend patroon. De auteur wil weten: Wat is de maximale dichtheid die je kunt bereiken zonder dat het een kristal wordt?

2. De Oplossing: Kijken naar de "Gaten" in plaats van de Koekjes

In plaats van naar de koekjes zelf te kijken, kijkt de auteur naar de gaten tussen de koekjes.

• Als je de middelpunten van drie aan elkaar liggende koekjes met elkaar verbindt, krijg je een driehoek.
• Als je vier koekjes verbindt, krijg je een vierkant.
• Deze driehoeken en vierkanten noemt hij "cellen".

De auteur introduceert een concept dat hij de Cell Order Distribution noemt. Denk hierbij aan een telling: hoeveel driehoekjes (3-cellen) en hoeveel vierkantjes (4-cellen) heb je in je stapel?

• Driehoekjes zijn erg compact. Als je alleen maar driehoekjes hebt, is je stapel heel dicht.
• Vierkantjes zijn wat minder compact.

De wiskunde laat zien: hoe meer driehoekjes je hebt, hoe voller je doos wordt.

3. De "Orde-Detectie" Regel

Het grootste gevaar is dat je stapel te geordend wordt. Als je te veel identieke driehoekjes naast elkaar hebt (bijvoorbeeld allemaal driehoekjes van drie kleine koekjes), begint je stapel te kristalliseren. Het wordt dan geen willekeurige stapel meer, maar een kristal.

De auteur bedacht een slimme alarmbel:

"Als een driehoekje gemiddeld meer dan één 'tweeling' (een identiek buur-driehoekje) heeft, is het gevaar te groot dat je een kristal vormt."

Hij gebruikt deze regel om een veilige zone te definiëren. Binnen deze zone is het garandeerd willekeurig. Buiten deze zone riskeer je dat je per ongeluk een kristal maakt.

4. De Resultaten: De Uitersten

De auteur berekent nu twee uitersten voor de dichtheid van de stapel:

• De Uiterste Bovenkant (Het Droomscenario):
  Stel dat je een stapel hebt die alleen uit driehoekjes bestaat, maar die net zo willekeurig is dat het alarmbel niet afgaat. Dit is de theoretisch hoogste dichtheid die mogelijk is. Het is als het "gouden punt" waar je zo dicht mogelijk bij elkaar zit zonder saai te worden.

  • Voorbeeld: Voor een bepaalde verhouding tussen grote en kleine koekjes, kun je tot wel 94% van de ruimte vullen.

• De Uiterste Onderkant (Het Zekere Scenario):
  Dit is de dichtheid die je altijd kunt bereiken, zelfs als je niet perfect bent in het vermijden van kristallen. Het is de "veilige ondergrens".

5. Wat betekent dit voor de praktijk?

De auteur laat zien dat er een "sweet spot" is. Als je een mengsel maakt van grote en kleine koekjes, moet je de verhouding precies goed kiezen.

• Als je te veel kleine koekjes neemt, krijg je te veel identieke driehoekjes en kristalliseer je.
• Als je te weinig kleine koekjes neemt, blijven er te veel grote gaten over.

Hij geeft een formule die precies aangeeft hoeveel kleine koekjes je nodig hebt voor elke grootte van de grote koekjes, om de dichtste mogelijke willekeurige stapel te krijgen.

Samenvattend: De Metafoor van de Tuin

Stel je voor dat je een tuin wilt volplanten met bloemen.

• De monodisperse situatie (alleen grote bloemen) is alsof je alleen rozen plant. Ze vallen snel in een perfect rijtje (kristal).
• De bidisperse situatie (grote en kleine bloemen) is alsof je rozen en viooltjes mixt.
• De wiskunde van de auteur is als een tuinarchitect die precies zegt: "Als je 30% viooltjes en 70% rozen gebruikt, en je plant ze zo dat er geen lange rijen identieke bloemen naast elkaar staan, kun je de tuin 94% vullen met bloemen."

Als je meer viooltjes plant, krijg je een rij viooltjes (kristal). Als je er minder plant, blijven er gaten over.

Conclusie:
Deze paper geeft voor het eerst een exact wiskundig antwoord op de vraag: "Hoe vol kan een willekeurige stapel van twee soorten schijven worden?" Het lost het probleem op door niet naar de chaos te kijken, maar naar de geometrie van de gaten tussen de schijven, en biedt een veiligheidsnet om te voorkomen dat de chaos per ongeluk in orde verandert.

Technische samenvatting

Titel: Random close packing fractie van bidisperse schijven: Theoretische afleiding en exacte grenzen

1. Het Probleem

Het voorspellen van de dichtste mogelijke willekeurige (ongeordende) pakking, bekend als Random Close Packing (RCP), van monodisperse (gelijkgrootte) schijven in 2D en bollen in 3D, is een langdurig onopgelost probleem in de natuurkunde en wiskunde.

• Uitdaging: Experimenten en simulaties geven vaak waarden voor de pakkingfractie ($\phi_{RCP}$) in een smal bereik (bijv. 0,82–0,86 voor 2D schijven), maar er is geen theoretisch bewijs dat dit het absolute maximum is.
• Beperkingen van bestaande methoden:
  • Het parameter-ruimte van mogelijke pakkingprocessen is oneindig groot, waardoor trial-and-error benaderingen ondoenlijk zijn.
  • Het garanderen van "wanorde" is moeilijk; gelijkgrootte deeltjes kristalliseren vaak tot geordende structuren (zoals hexagonale roosters) die dichter zijn dan wanordelijke toestanden.
  • Bestaande wanorde-criteria zijn vaak a-posteriori (ze testen wanorde nadat het systeem is gepakt), wat leidt tot inefficiënte zoektochten.
• Specifiek doel: Het artikel richt zich op bidisperse schijven (twee verschillende groottes). Dit wordt vaak gebruikt in experimenten om kristallisatie te minimaliseren. Het doel is om een theoretische bovengrens af te leiden voor $\phi_{RCP}$ als functie van de grootteverhouding ($D$) en de concentratie van de kleine schijven ($p$).

2. Methodologie

De auteur gebruikt een theoretisch raamwerk gebaseerd op de Cell Order Distribution (COD) om de oneindige parameter-ruimte te omzeilen.

• Cell Order Distribution (COD):

  • Door de centra van contacterende schijven te verbinden, ontstaat een grafiek. De kleinste gaten (voids) die door deze lijnen worden omsloten, worden "cellen" genoemd.
  • De "orde" ($k$) van een cel is het aantal randen die deze cel omringen. De COD ($Q_k$) is de verdeling van het fractie van cellen met orde $k$.
  • Voordeel: De COD correleert direct met de dichtheid en fungeert als een parameter die alle mogelijke pakkingen beschrijft, ongeacht het specifieke proces dat ze heeft gegenereerd.

• A-priori Wanorde-criterium:

  • Om kristallisatie te voorkomen zonder het systeem eerst te genereren, stelt de auteur een criterium op: een cel van identieke schijven mag gemiddeld slechts één identieke buur hebben.
  • Dit beperkt de kans op het vormen van grote clusters van identieke 3-cellen (die leiden tot kristallen).
  • Dit criterium definieert een toegestaan bereik voor de concentratie $p$ voor elke grootteverhouding $D$, binnen welke wanorde gegarandeerd is.

• Afleiding van de Pakkingfractie:

  • De totale pakkingfractie $\phi$ wordt uitgedrukt als een functie van de gemiddelde celoppervlakken en de oppervlakten die door de schijven worden ingenomen binnen deze cellen.
  • Voor een pakking bestaande uit alleen 3-cellen (de dichtst mogelijke configuratie) wordt een exacte formule afgeleid.
  • De auteur introduceert een hulpvariabele $u$ (een effectieve concentratie binnen de 3-cellen) om de onderlinge afhankelijkheid van aangrenzende cellen te modelleren.

3. Belangrijkste Bijdragen

1. Theoretisch kader voor bidisperse pakkingen: De eerste analytische afleiding van de maximale willekeurige pakkingfractie voor bidisperse schijven in 2D, waarbij $D$ en $p$ expliciet worden meegenomen.
2. Garantie voor wanorde: Een nieuw, a-priori criterium dat het bereik van $p$ definieert waarbinnen een pakking wiskundig gegarandeerd wanordelijk is, zonder kristallisatie.
3. Exacte grenzen: Afleiding van exacte boven- en ondergrenzen voor $\phi_{RCP}(p, D)$.
4. Validatie van de COD-methode: Demonstratie dat de COD-methode robuust is en niet afhankelijk van het specifieke pakkingproces.

4. Resultaten

• Bovengrens ($\phi_3$): De theoretisch hoogst mogelijke pakkingfractie wordt bereikt door een pakking die uitsluitend bestaat uit 3-cellen. Hoewel dergelijke volledig wanordelijke pakkingen topologisch misschien niet mogelijk zijn, vormen ze een exacte bovengrens.
  • De maximale waarde wordt bereikt bij een specifieke concentratie $p_{max}$ die afhangt van $D$.
  • Voor $D=1$ (monodisperse) wordt de waarde van het trigonale kristal ($\pi/2\sqrt{3} \approx 0.907$) hersteld als bovengrens.
• Ondergrens ($\phi_{low}$): Een exacte ondergrens wordt berekend door aan te nemen dat de pakking bestaat uit 3- en 4-cellen, waarbij het aandeel 3-cellen ($Q_3$) zo laag is dat wanorde nog steeds gegarandeerd is (volgens het gestelde criterium).
• Gedrag van $\phi_{RCP}$:
  • De berekende $\phi_{RCP}(p, D)$ valt binnen de afgeleide boven- en ondergrenzen.
  • Voor een groot deel van het bereik van $p$ (waar wanorde gegarandeerd is) coincideert de berekende $\phi_{RCP}$ met de bovengrens $\phi_3$. Dit betekent dat in dit bereik de dichtst mogelijke wanordelijke toestand topologisch toegankelijk is.
  • Invloed van $D$: Er is een kritieke grootteverhouding ($D \approx 3$). Voor $D \lesssim 3$ en wanneer de twee schijftypen een gelijk oppervlak innemen (een veelgemaakte keuze in experimenten), neemt het risico op kristallisatie toe en wordt het bereik van veilige $p$-waarden smaller.
• Numerieke waarden: De paper levert numerieke tabellen en grafieken voor $1 < D < D_{max}$ (waar $D_{max} = \sqrt{3}/(2-\sqrt{3})$ zodat er geen "rattlers" in 3-cellen passen).

5. Betekenis en Toepassingen

• Richting voor Experimenten: De paper biedt een leidraad voor het ontwerpen van experimenten en simulaties. Door de concentratie $p$ te kiezen binnen het bereik waar wanorde gegarandeerd is en dicht bij $p_{max}$, kunnen onderzoekers proberen dichter bij de theoretische limiet te komen.
• Detectie van Kristallisatie: Als een experimentele pakking een fractie oplevert die hoger is dan de hier voorspelde $\phi_{RCP}$, betekent dit dat er uitgebreide kristallijne domeinen aanwezig zijn, zonder dat men de structuur in detail hoeft te analyseren.
• Verbetering van Bestaande Literatuur: De afgeleide grenzen verbeteren de bestaande schattingen in de literatuur en bieden een wiskundig onderbouwde basis in plaats van puur empirische waarden.
• Uitbreidbaarheid: De methode is in principe uitbreidbaar naar tridisperse of polydisperse systemen, hoewel de wiskundige complexiteit exponentieel toeneemt met het aantal deeltjessizes.

Conclusie:
Blumenfeld presenteert een doorbraak in het theoretisch begrijpen van willekeurige pakkingen door de complexiteit van het proces te reduceren tot een statistische verdeling van celordes. De afgeleide formules geven voor het eerst exacte, proces-onafhankelijke grenzen voor de dichtheid van bidisperse schijven in 2D, wat essentieel is voor het optimaliseren van materialen en het begrijpen van faseovergangen in granulaire systemen.