A Theoretical Framework for the Coupling of Macroscale-Nanoscale Mechanochemical Phenomena in Condensed Matter

Technische Samenvatting: Een Theoretisch Kader voor de Koppeling van Macro- naar Nano-schaal Mechanochemische Fenomenen in Gecondenseerde Materie

Probleemstelling
Het vakgebied van de covalente mechanochemie is getransformeerd van fundamentele wetenschap naar technische toepassingen, maar het ontbreekt aan een robuust theoretisch kader voor het voorspellen van reactiekinetiek in gecondenseerde materie. Bestaande analytische modellen, zoals het Bell-model en het getitleerde Potentiële Energie Oppervlak (PES)-model, berusten op de aanname dat de toegepaste kracht ($F$) exact langs de reactiecoördinaat ($\xi$) werkt. In realistische technische scenario's waarbij gecondenseerde materie betrokken is, induceren macroscopische rekken complexe, niet-lineaire en meerdimensionale moleculaire deformaties (combinaties van buiging, torsie en rek). In deze scenario's is de toegepaste kracht zelden uitgelijnd met de intrinsieke reactiecoördinaat, waardoor standaard perturbatiemodellen ontoereikend zijn voor het voorspellen van veranderingen in de activeringsbarrière. De primaire kenniskloof is de afwezigheid van een theoretisch kader dat in staat is om macroscopische krachten te transduceren naar covalente bindingen wanneer de deformatie geen zuivere trekkracht langs het reactiepad is.

Methodologie
De auteurs ontwikkelen een niet-perturbatief theoretisch kader waarin de verandering in de activeringsbarrière ($\Delta E^\ddagger$) wordt uitgedrukt als een functie van de stijging in moleculaire rekenergie ($\Delta U_{sys}$), in plaats van direct als een functie van de toegepaste kracht. Deze benadering herformuleert het probleem met behulp van het activatie-rekmodel (activation-strain model), waarbij de barrière wordt gewijzigd door de rek die nodig is om de overgangstoestand te bereiken.

De kern van de methodologie bestaat uit het uitbreiden van de afgeleide van de activeringsbarrière met betrekking tot de systeemrekenergie ($\frac{\partial \Delta E^\ddagger_o}{\partial \Delta U_{sys}}$) met behulp van de kettingregel naar twee afzonderlijke vormen:

1. Kracht-gebaseerde expansie: $\frac{\partial \Delta E^\ddagger_o}{\partial F_\xi} \frac{\partial F_\xi}{\partial \Delta U_{sys}}$. Deze vorm is toepasbaar wanneer collectieve variabelen expliciete krachtanalyse langs het reactiepad mogelijk maken.
2. Coördinaat-gebaseerde expansie: $\frac{\partial \Delta E^\ddagger_o}{\partial \xi} \frac{\partial \xi}{\partial \Delta U_{sys}}$. Deze vorm wordt gepresenteerd als praktischer voor simulatie-workflows en hoogwaardige kwantumchemische methoden. Het berust op de observatie dat $\frac{\partial \Delta E^\ddagger_o}{\partial \xi}$ vaak onafhankelijk is van de specifieke deformatiemodus ($\epsilon$), waardoor het eenmaal gesampled kan worden, terwijl $\frac{\partial \xi}{\partial \Delta U_{sys}}$ kan worden geëxtraheerd uit steered molecular dynamics (MD) of geometrische optimalisaties onder beperkingen.

Het kader werd gevalideerd met behulp van klassieke moleculaire dynamica-simulaties van het mechanochromische polymeer spiropyrane (SP). Vijf verschillende deformatiepaden werden gesimuleerd met het ReaxFF-krachtveld binnen LAMMPS:

• Trekken langs de ruggengraat (baseline).
• "Centrale Twist" (rotatie rond het spiro-atoom).
• "Secundaire Ring Buckle" (out-of-plane buckling van de rechter ring).
• "Centrale Buckle" (in-plane vouwen van de twee helften).
• Twee "Penta" ring-gebaseerde buckle-deformaties.

De CoGEF-methode (Constrained Geometries to Simulate External Force) werd gebruikt om bindingsbreuk te induceren, terwijl Many-Bodied Steered Molecular Dynamics (MBsMD) de aanvullende niet-reactieve deformaties induceerde.

Belangrijkste Resultaten

• Ontkoppeling van Rek en Reactie: De simulaties toonden aan dat verschillende deformatiepaden resulteren in zeer verschillende veranderingen aan de activeringsbarrière voor dezelfde stijging in interne energie ($\Delta U_{sys}$).
• Validatie van de Kettingregel: De studie bevestigde dat $\frac{\partial \Delta E^\ddagger_o}{\partial \xi}$ (de gevoeligheid van de barrière voor de reactiecoördinaat) grotendeels onafhankelijk is van de specifieke deformatiemodus. Data van alle vijf de deformatiepaden voldeden aan één enkele lineaire trend wanneer ze tegen $\xi$ werden uitgezet, wat de aanname valideert dat deze term benaderd kan worden door een eenvoudige lineaire functie.
• Projectie-effecten: De term $\frac{\partial \xi}{\partial \Delta U_{sys}}$ vangt effectief hoe een specifieke macroscopische deformatie projecteert op de reactiecoördinaat.
  • De "Secondary Ring Buckle" vertoonde minimale projectie op de C-O binding rekcoördinaat, wat resulteerde in een verwaarloosbare reductie van de barrière.
  • "Central Twist" en "Central Buckle" vertoonden sterke projectie, waardoor de barrière werd verlaagd. "Central Twist" vertoonde super-lineaire groei (hoge efficiëntie), terwijl "Central Buckle" een sub-lineaire respons vertoonde (dalende meeropbrengst) vanwege structurele relaxatiemechanismen die de bindingsrek onderdrukken.
• Computationele Efficiëntie: Door het probleem te scheiden in een reactiecoördinaat-gevoeligheidsterm en een rek-mappingterm, vermindert het kader de computationele kosten. Het elimineert de noodzaak om voor elke unieke deformatiemodus en -grootte een dure volledige reactiepad-sampling uit te voeren.

Betekenis en Claims
Het artikel claimt een algemeen, parametriseerbaar theoretisch kader te bieden dat de kloof overbrugt tussen macroscopische rek en nano-schaal mechanochemische reactiviteit. De primaire betekenis ligt in:

1. Generaliseerbaarheid: In tegen tegenstelling tot eerdere modellen, gaat dit kader er niet vanuit dat de toegepaste kracht uitgelijnd is met de reactiecoördinaat, waardoor het geschikt is voor complexe, niet-lineaire deformaties in gecondenseerde materie.
2. Toegankelijkheid: Door de barrièreverandering uit te drukken in termen van algemene observeerbare grootheden (rekenergie en reactiecoördinaatverplaatsing), maakt de methode de voorspelling van mechanochemische effecten mogelijk met eenvoudige structuuroptimalisatieberekeningen. Dit vergemakkelijkt het gebruik van hoogwaardige kwantumchemische methoden die voorheen te computationeel duur waren voor dergelijke multiscale problemen.
3. Interpreteerbaarheid: De benadering biedt een analytische, interpreteerbare manier om te begrijpen hoe specifieke intra-moleculaire rektoestanden de chemische reactiviteit beïnvloeden, waarmee verder wordt gekeken dan "black box"-simulaties.

De auteurs merken op dat hoewel de methode breed bruikbaar is, toekomstig werk vereist is om statistische variantie, fluctuaties bij eindige temperatuur en het modelleren van distributies van intra-moleculaire rektoestanden binnen representatieve volumeelementen in grootschalige systemen van gecondenseerde materie aan te pakken.