A Chemical Space Perspective on Diastereomeric Barriers in Alkylperoxy-to-Hydroperoxyalkyl Isomerization

Dit onderzoek benadrukt het cruciale belang van stereochemie in de lage-temperatuur auto-oxidatie van koolwaterstoffen door aan te tonen dat diastereomere barrières in de isomerisatie van alkylperoxy- naar hydroperoxyalkyl-radicalen variëren van degeneraat tot meer dan 60 kcal/mol, wat aantoont dat constitutioneel samengevoegde moleculaire representaties kinetisch relevante reactiekanalen systematisch kunnen missen.

De kern

Stel je voor dat je een gigantische, chaotische fabriek hebt waar brandstof (zoals benzine) wordt omgezet in energie. In deze fabriek werken miljoenen kleine, onzichtbare robotjes (moleculen) die constant van baan wisselen. Soms botsen ze, soms splitsen ze, en soms veranderen ze van vorm.

Dit artikel van Raghunathan Ramakrishnan gaat over een heel specifiek, maar cruciaal moment in deze fabriek: het moment waarop een molecuul een zuurstofatoom oppikt en vervolgens zijn eigen atomen herschikt om een nieuwe vorm aan te nemen. Dit proces is essentieel voor hoe auto's branden en hoe de lucht in onze atmosfeer veroudert.

Hier is de kern van het verhaal, vertaald naar alledaagse taal met een paar creatieve vergelijkingen:

1. Het Probleem: De "Spiegelbeeld"-Illusie

Stel je voor dat je een knoop in een touw moet losmaken. Er zijn twee manieren om dat te doen: je kunt het touw naar links draaien of naar rechts. Als je touw perfect symmetrisch is, maakt het niet uit welke kant je kiest; het resultaat is hetzelfde.

Maar wat als je touw een knoop heeft die eruitziet als een hand? Een linkerhand en een rechterhand zijn spiegelbeelden, maar ze zijn niet uitwisselbaar. Je kunt je linkerhand niet in een rechterhandschoen steken.

In de chemie noemen we dit stereochemie. De meeste computersimulaties van brandstofverbranding kijken alleen naar de "kabels" (welke atomen aan welke vastzitten), maar negeren de "handen" (de 3D-oriëntatie). Ze denken: "Oh, dit is dezelfde reactie, dus het gaat even snel."

De auteur zegt: "Nee, dat is niet zo!"

2. De Ontdekking: Twee Wegen, Twee Kosten

De auteur en zijn team hebben een enorme database gemaakt (ze noemen het SEARS, wat klinkt als een supermarkt voor chemische reacties). Ze hebben gekeken naar bijna 5.000 verschillende moleculen en hun reactiepaden.

Ze ontdekten iets verrassends:
Soms heeft een molecuul twee mogelijke routes om van A naar B te gaan. Deze routes lijken op het eerste gezicht hetzelfde, maar omdat de moleculen in 3D-ruimte anders liggen, zijn het eigenlijk twee verschillende "deuren".

• Situatie A (De gelijke deuren): Soms zijn beide deuren even groot en even makkelijk te openen. Dan telt je computer twee keer mee, en de reactie gaat twee keer zo snel.
• Situatie B (De gesloten deur): Soms is één deur een smalle, krampachtige spleet (door drukte in de moleculaire "kamer") en de andere een brede, open poort. De moleculen die proberen door de smalle spleet te gaan, stoten zich de kop. Die route is zo moeilijk (soms wel 60 keer moeilijker!) dat hij in de praktijk niet bestaat.

De metafoor:
Stel je voor dat je een drukke treinreis maakt.

• De oude manier: De computer zegt: "Er is één trein naar Amsterdam."
• De nieuwe manier: De auteur zegt: "Er zijn twee treinen. De ene is een snelle, comfortabele hogesnelheidstrein. De andere is een oude, stoffige fiets die door een modderig veld moet. Als je alleen naar de 'trein' kijkt, mis je dat de fiets bijna niemand haalt."

3. Waarom is dit belangrijk?

Als je een computermodel maakt om te voorspellen hoe een motor werkt of hoe de lucht vervuilt, en je negeert deze "handen" (de stereochemie), dan maak je twee fouten:

1. Je denkt dat een reactie te snel gaat (als je de moeilijke route meetelt die in werkelijkheid niet gebeurt).
2. Je denkt dat een reactie te langzaam gaat (als je twee gemakkelijke routes als één ziet, terwijl ze samen de snelheid verdubbelen).

De studie toont aan dat voor veel moleculen, vooral die in ringen of met veel vertakkingen, deze "twee deuren" er echt toe doen. Het maakt niet uit of je linksom of rechtsom draait; de energie die je nodig hebt om die draai te maken, kan enorm verschillen.

4. De Oplossing: De "SEARS"-Supermarkt

De auteur heeft een nieuwe methode ontwikkeld om deze moleculen niet als platte tekeningen te bekijken, maar als 3D-figuurtjes. Ze hebben een dataset gemaakt waar elke mogelijke "hand" en elke mogelijke "deur" apart wordt geteld.

Ze hebben bewezen dat:

• Soms de twee routes bijna even goed zijn (dan verdubbelt de snelheid).
• Soms één route zo slecht is dat hij vergeten kan worden.
• De grootte van dit verschil afhangt van hoe "vol" de moleculaire kamer is (sterische druk). Als het er druk is, wordt de ene route een onbegaanbare doolhof.

Conclusie in één zin

Dit onderzoek zegt tegen de chemici en ingenieurs: "Stop met het tellen van moleculen alsof ze platte papieren zijn; ze zijn 3D-acteurs, en soms is hun linkerhand een superkracht, terwijl hun rechterhand een onoverkomelijke muur is."

Door dit in de gaten te houden, kunnen we betere motoren bouwen, schoner branden en de atmosfeer beter begrijpen. Het is een stap van "plat denken" naar "diep, 3D-denken" in de wereld van chemie.

Technische samenvatting

Probleemstelling

Bij de auto-oxidatie van koolwaterstoffen bij lage temperaturen spelen radicaal-intermediaten een cruciale rol. Een van de bepalende stappen is de isomerisatie van alkylperoxy-radicalen ($ROO^\bullet$) naar hydroperoxyalkyl-radicalen ($^\bullet QOOH$) via intramoleculaire waterstofoverdracht. De kinetiek van deze reacties hangt sterk af van de moleculaire structuur.

Het huidige probleem in de verbrandingschemie en geautomatiseerde mechanismegeneratie is dat stereochemische informatie vaak wordt genegeerd of "ingestort" (collapsed) naar constitutionele connectiviteit. Hoewel symmetrie-equivalente paden correct worden behandeld via degeneratiefactoren, worden diastereomere paden vaak ten onrechte als één effectief pad behandeld. Wanneer een $ROO^\bullet$-intermediair chirale centra bevat, kan de waterstofoverdracht via twee verschillende diastereomere overgangstoestanden (TS) plaatsvinden die leiden tot hetzelfde $^\bullet QOOH$-product, maar die verschillen in de driedimensionale rangschikking van atomen. Als deze verschillen niet expliciet worden behandeld, kunnen kinetisch relevante reactiekanalen worden gemist, wat leidt tot onnauwkeurige snelheidsconstanten en voorspellingen voor auto-oxidatie en verbranding.

Methodologie

De auteurs hebben een nieuwe dataset ontwikkeld, genaamd SEARS (Stereochemically Expanded Autooxidation Reaction Space), om dit probleem aan te pakken. De methodologie omvat de volgende stappen:

1. Dataset Selectie en Enumeratie:

   • Startpunt was de bigQM7$\omega$-dataset (12.880 verbindingen).
   • Er werden 498 niet-aromatische koolwaterstoffen (C1–C7) geselecteerd, exclusief aromaten en sterk gespannen gefuseerde ringen.
   • Er werden stereochemisch onderscheiden isomeren gegenereerd (geometrische isomeren en diastereomeren), resulterend in 590 unieke parent-structuren.
   • Van deze parent-structuren werden systematisch $R^\bullet$, $ROO^\bullet$ en $^\bullet QOOH$-radicalen gegenereerd.

2. Generatie van Overgangstoestanden (TS):

   • Voor de $ROO^\bullet \rightarrow ^\bullet QOOH$-isomerisatie werden cyclische overgangstoestanden (ringgrootte 4 tot 8 atomen) geconstrueerd.
   • Een unieke aanpak werd gebruikt waarbij een tijdelijk stikstofatoom ([N]) in de SMILES-structuur fungeerde als brug tussen het zuurstofatoom van de peroxygroep en het koolstofatoom dat het waterstofdonor is. Dit dwong de cyclische topologie af tijdens de constructie.
   • Na generatie van de 3D-geometrieën werd het stikstofatoom vervangen door waterstof voor de kwantumchemische berekeningen.
   • Dit resulteerde in 5.590 kandidaat-TS-structuren, waarvan 3.808 succesvol als 3D-geometrieën konden worden gegenereerd.

3. Kwantumchemische Berekeningen (DFT):

   • Optimalisatie: Geometrie-optimalisaties en frequentie-analyses werden uitgevoerd met behulp van het $\omega$B97M-D4 functioneel en het def2-SV(P) basisstel (via ORCA).
   • Validatie: Intrinsic Reaction Coordinate (IRC) berekeningen werden uitgevoerd om te verifiëren dat elke gevonden TS daadwerkelijk de beoogde $ROO^\bullet$ en $^\bullet QOOH$ minima verbindt.
   • Energieberekening: Eindelijke elektronische energieën werden berekend op het $\omega$B97M-D4/def2-TZVP niveau.
   • Filtering: Structuren die tijdens de IRC naar een ongewenst kanaal (zoals O-O bindingsbreuk) evolueerden, werden verwijderd.

Belangrijkste Bijdragen

• SEARS Dataset: Een grote, openbare dataset met 5.356 gespecificeerde soorten, inclusief 2.324 gevalideerde diastereomere overgangstoestanden die corresponderen met 1.162 unieke $ROO^\bullet \rightarrow ^\bullet QOOH$ reacties.
• Expliciete Stereochemische Behandeling: Het eerste grootschalige bewijs dat diastereomere paden systematisch worden onderscheiden in plaats van te worden samengevoegd, wat essentieel is voor nauwkeurige kinetische modellering.
• Validatie van "Ephemeral Diastereomers": Het bevestigen dat deze stereochemische verschillen zich uitsluitend manifesteren in de overgangstoestand en niet in het isolabele intermediair, maar toch een grote impact hebben op de reactiesnelheid.

Resultaten

De analyse van de dataset leverde de volgende inzichten op:

1. Energieverschillen (Diastereomere Splijting):

   • De energieverschillen tussen de twee diastereomere TS-paden ($\Delta\Delta E_{barrier}$) variëren enorm: van degeneraat (0 kcal/mol) tot meer dan 60 kcal/mol.
   • Hoewel veel paren klein zijn gesplitst (32% heeft een verschil $\le$ 1 kcal/mol), is er een significante staart van grote verschillen.
   • Grote verschillen treden vaak op wanneer de peroxygroep gebonden is aan een $sp^3$-gehybridiseerd koolstofatoom en de waterstof wordt onttrokken aan een secundair koolstofcentrum, vooral in cyclische systemen waar ring-beperkingen en sterische hindering de ene diastereomeer sterk destabiliseren.

2. Invloed op Kinetiek:

   • Wanneer de barrières vergelijkbaar zijn (nabij-degeneraat), dragen beide paden bij aan de totale snelheid, wat de totale snelheid kan verdubbelen vergeleken met het nemen van slechts één pad.
   • Bij grote energieverschillen (> 5-10 kcal/mol) wordt het pad met de hogere barrière kinetisch verwaarloosbaar.
   • Het negeren van deze splitsing leidt tot systematische fouten in de schatting van reactiesnelheden.

3. Structuur-Activiteit Relaties:

   • Ringgrootte: 6-ring TS's hebben over het algemeen lagere barrières dan 5-ring TS's, maar de diastereomere splijting hangt meer af van lokale stereochemie en spanning dan alleen van de ringgrootte.
   • Hybridisatie: Systemen met een peroxygroep op een $sp^2$-koolstof tonen over het algemeen kleinere splijtingen dan die op $sp^3$-koolstof.

Significantie en Conclusie

De studie concludeert dat constitutioneel samengevoegde moleculaire representaties systematisch kinetisch relevante reactiekanalen missen. Voor de ontwikkeling van nauwkeurige verbrandingsmodellen en geautomatiseerde mechanismegeneratie is het essentieel om stereochemie expliciet te behandelen, zelfs voor kortstondige overgangstoestanden.

• Voor Verbrandingsmodellering: Het SEARS-dataset biedt een fundament voor het ontwikkelen van stereochemie-bewuste snelheidsregels en het verbeteren van voorspellende modellen voor auto-oxidatie.
• Voor Geautomatiseerde Workflows: Bestaande workflows die alleen constitutionele connectiviteit volgen, moeten worden uitgebreid om diastereomere paden te onderscheiden.
• Toekomstige Toepassingen: De dataset dient als een waardevol benchmark voor machine-learning modellen die niet alleen gemiddelde barrières moeten leren, maar ook de complexe patronen van stereochemisch geïnduceerde splitsing en nabije degeneratie.

Kortom, stereochemie is geen louter labelkwestie, maar een fundamentele bron van barrièrediversiteit die de reactiekinetiek van koolwaterstoffen diepgaand beïnvloedt.