Rate Equation for the Transfer of Interstitials across Interfaces between Equilibrated Crystals

Technische Samenvatting: Snelheidsvergelijking voor de Overdracht van Interstitials over Interfaces tussen Gelijkwichtscristallen

Probleemstelling
Het artikel behandelt een discrepantie tussen experimentele waarnemingen en bestaande kinetische modellen met betrekking tot het laden en ontladen van interstitiële vaste oplossingen, specifiek metaalhydriden (bijv. Pd, PdPt). Experimenten wijzen uit dat de karakteristieke tijd ($T$) voor compositieveranderingen drastisch toeneemt (met ordes van grootte) wanneer het compositiebereik een faseovergang of de nabijheid van een kritisch punt van een mengbaarheidsgat omvat. De standaard Butler-Volmer (BV) vergelijking, die breed wordt gebruikt voor het analyseren van elektrochemisch gecontroleerde hydreringskinetiek, slaagt er niet in dit fenomeen te verklaren. Het BV-raamwerk gaat uit van overdracht tussen locaties met een vaste a priori bezetting en grondtoestandsenergie, waarbij de drijvende kracht uitsluitend wordt behandeld als het verschil in chemische potentialen ($\Delta \mu_S$). Echter, in nanoschaal metaalhydriden is interne diffusie orders van grootte sneller dan de interface-overdracht, wat betekent dat de bulkfasen quasi-instantaan een intern evenwicht bereiken. Bijgevolg vindt de interface-overdracht plaats tussen uitgebreide, ergodische ensembles van deeltjes, waarbij bezettingen en grondtoestandsenergieën fluctueren; een scenario dat niet wordt gevangen door standaard sprongkansmodellen.

Methodologie
De auteur leidt een nieuwe kinetische snelheidswet af op basis van statistische mechanica en het principe van gedetailleerd evenwicht, in plaats van fenomenologische sprongkansen. Het systeem wordt gemodelleerd als een composiet van twee interstitiële vaste oplossingfasen (A en B), gescheiden door een energiebarrière met een overgangstoestand ($T$).

1. Statistisch Raamwerk: De fasen worden behandeld als canonieke systemen bij een constante temperatuur ($T$) met geconserveerde aantallen deeltjes ($N_S$) en matrixlocaties ($N_M$). De bulkdiffusie wordt verondersteld snel genoeg te zijn om interne ergodiciteit binnen elke fase te handhaven.
2. Overgangstoestand Analyse: De waarschijnlijkheid van een "overdrachtconfiguratie" ($C_p$) — waarbij een deeltje zich in een geactiveerde staat bevindt op de bronlocatie en de doellocatie leeg is — wordt berekend door de sommatie over alle microtoestanden van het beperkte systeem.
3. Vrije Energie Decompositie: De vrije energieverandering ($\delta F$) die nodig is om de overdrachtconfiguratie te bereiken, wordt berekend via een vierstaps virtueel proces:
   • Het comprimeren van het deeltjesensemble in de bronfase om een locatie vrij te maken (entropische kosten).
   • Het comprimeren van het deeltjesensemble in de doelfase om een locatie vrij te maken (entropische kosten).
   • Het verwijderen van een deeltje uit de bronfase naar een referentietoestand (chemische potentiaal arbeid).
   • Het invoegen van het deeltje in de specifieke overgangstoestand (energiekosten).
4. Flux Afleiding: De netto fluxdichtheid ($J_B$) wordt afgeleid als het verschil tussen de voorwaartse en achterwaartse fluxen, gebruikmakend van de berekende waarschijnlijkheden van de overdrachtconfiguraties.

Belangrijkste Bijdragen en Resultaten
De primaire bijdrage is de afleiding van een nieuwe snelheidswet (Vergelijking 15) voor de interface-flux van ladingsneutrale deeltjes:
$$J_B = \nu_0 \rho_T (1 - x_A)(1 - x_B) \exp\left(-\frac{\epsilon_T}{k_B T}\right) \left( \exp\left(\frac{\mu_S^A}{k_B T}\right) - \exp\left(\frac{\mu_S^B}{k_B T}\right) \right)$$

Kernkenmerken van dit resultaat zijn:

• Expliciete Afhankelijkheid van Chemische Potentialen: In tegenstelling tot de BV-vergelijking, die afhangt van het verschil in chemische potentialen ($\Delta \mu_S$), hangt deze snelheidswet expliciet af van de individuele chemische potentialen ($\mu_S^A$ en $\mu_S^B$) van de twee fasen.
• Afhankelijkheid van Vacaturefractie: De flux wordt expliciet gemoduleerd door de vacaturefracties $(1-x_A)$ en $(1-x_B)$, wat de statistische waarschijnlijkheid reflecteert om een lege doellocatie te vinden in een ergodisch ensemble.
• Gedetailleerd Evenwicht: De vergelijking voldoet van nature aan het principe van gedetailleerd evenwicht, en wordt nul wanneer $\mu_S^A = \mu_S^B$.
• Verklaring van Kinetische Vertraging: Het artikel demonstreert dat nabij een kritisch punt van een mengbaarheidsgat of spinodale curve de deeltjessusceptibiliteit ($\chi_S = dx/d\mu_S$) divergeert. Aangezien de karakteristieke tijd $T$ omgekeerd evenredig is aan het product van de uitwisselingsfluxdichtheid en de susceptibiliteit ($T^{-1} \propto J_{eq}^\mu / \chi_S$), leidt de divergentie van $\chi_S$ tot een divergentie in $T$. Dit biedt een theoretische verklaring voor de experimenteel waargenomen drastische vertraging in laadtijden nabij kritische punten.
• Asymmetrie: De afgeleide wet voorspelt een asymmetrie in de karakteristieke tijd $T$ aan weerszijden van het kritische punt, een kenmerk dat experimenteel wordt waargenomen maar niet wordt gereproduceerd door de symmetrische BV-vergelijking.

Betekenis
Het artikel claimt dat de afgeleide snelheidswet een fysisch consistenter beschrijvend model biedt voor de interface-overdracht in systemen waar bulkfasen sneller in evenwicht komen dan de interface. Door de overdrachtslocaties te koppelen aan ergodische canonieke systemen in plaats van aan vaste locaties, houdt het model rekening met de fluctuerende grondtoestandsenergieën en bezettingen die inherent zijn aan interagerende vaste oplossingen. Het werk suggereert dat het falen van de Butler-Volmer vergelijking om kinetische anomalieën in metaalhydriden te verklaren, voortkomt uit het negeren van de expliciete afhankelijkheid van evenwichtscomposities en de statistische aard van de beschikbaarheid van vacatures. De resultaten bieden een theoretische basis voor het begrijpen van laad-/ontlaadkinetiek in energieopslagmaterialen, met name nabij faseovergangen, zonder dat daarvoor ad-hoc aanpassingen aan de overdrachtscoëfficiënt of de barrière-energie nodig zijn.