Cis--Trans Rotational Isomerism of Seleno-, Thio-, and Formic Acids and Their Dimers: Chemical Kinetics under Interstellar Conditions

Technische Samenvatting: Cis–Trans Rotatie-isomerie van seleno-, thio- en mierenzuren en hun dimeren

Probleemstelling
De chemische kinetiek van rotatie-isomerie (cis–trans interconversie) in mierenzuur (HC(O)OH) en zijn zwavel- (HC(O)SH) en seleniumanalogen (HC(O)SeH) blijft een onderwerp van aanzienlijk debat binnen de astrochemie. Hoewel de trans-conformeer thermodynamisch stabiel is, is de cis-conformeer waargenomen in diverse interstellaire omgevingen, waaronder reusachtige moleculaire wolken en donkere moleculaire wolken zoals TMC-1. Echter, theoretische studies hebben tegenstrijdige resultaten opgeleverd met betrekking tot de snelheid van de cis-naar-trans isomerisatie. Sommige berekeningen suggereren dat de reactiesnelheidsconstanten te laag zijn ($k \le 4 \times 10^{-21} \text{ s}^{-1}$) om interconversie binnen de leeftijd van het universum mogelijk te maken, wat in strijd is met observationele gegevens. Daartegenover staan laboratoriumstudies op cryogene edelgasmatrices die een snelle isomerisatie aangeven via kwantumtunneling, met halfwaardetijden variërend van minuten tot uren, afhankelijk van de matrixomgeving. Er bestaat een kritische kloof in het begrip van de intrinsieke gasfase-kinetiek van deze monomeren en hun dimeren om als baseline te dienen voor het interpreteren van matrixeffecten en interstellaire chemie.

Methodologie
Deze studie maakt gebruik van hoogwaardige kwantumchemische berekeningen om de chemische kinetische reactiesnelheidsconstanten te bepalen voor de cis-naar-trans isomerisatie van monomerische en dimere mierenzuur-, thioformezuur- en selenoformezuur-varianten over een temperatuurbereik van 10 K tot 300 K.

• Elektronische Structuur: Geometrie-optimalisaties en frequentieberekeningen (harmonisch en anharmonisch) werden uitgevoerd met de dubbel-hybride functionaal B2PLYP met de Grimme's D3(BJ) dispersiecorrectie en de def2-TZVP basisset.
• Energieverfijning: Hogere-orde single-point energieën werden berekend met de DLPNO-CCSD(T) methode om nauwkeurigheid voor kinetische doeleinden te waarborgen. Voor dimere soorten werd B3LYP met D3(BJ) dispersie gebruikt voor robuustheid, hoewel Basis Set Superposition Error (BSSE) correcties om technische redenen werden genegeerd.
• Kinetische Theorie: Reactiesnelheidsconstanten werden berekend met behulp van Canonical Variational Transition State Theory (CVT), aangevuld met kwantummechanische tunneling-correcties. Specifiek maakte de studie gebruik van Small-Curvature Tunnelling (SCT), Quantised Reaction States (QRC) en het Interpolated Single-Point Energies (ISPE) algoritme (samengevat als CVT/SCT).
• Reactiecoördinaat: Intrinsic Reaction Coordinate (IRC) berekeningen werden gegenereerd op het B2PLYP–D3BJ/def2-TZVP niveau, gemoduleerd door de hoogwaardige CCSD(T) energieën.

Belangrijkste Resultaten

• Monomeer Kinetiek: De berekende reactiesnelheidsconstanten voor de cis-naar-trans isomerisatie van de monomeren laten zien dat kwantummechanische tunneling domineert bij lage temperaturen.

  • Bij 10 K is de berekende halfwaardetijd van cis-mierenzuur ongeveer 7,6 ms. Dit is aanzienlijk sneller dan eerdere gasfase theoretische voorspellingen (die levensduur van miljarden jaren suggereerden), maar langzamer dan sommige cryogene matrixexperimenten (die seconden tot uren rapporteren).
  • De reactie wordt gedreven door de barrièrebreedte in plaats van de barrièrehoogte. Ondanks vergelijkbare barrièrehoogtes voor HC(O)OH en HC(O)SeH, resulteert de smallere barrièrebreedte van de zuurstofsoort ($3,3 \text{ amu}^{1/2}\cdot\text{bohr}$) vergeleken met zwavel en selenium ($5,6$ en $6,3 \text{ amu}^{1/2}\cdot\text{bohr}$, respectievelijk) in een veel hogere tunneling-efficiëntie voor mierenzuur. Dit illustreert een "tweede-rij anomalie" waarbij zuurstof anders gedraagt dan de zwaardere Groep 16-elementen.
  • Thio- en selenoformezuren vertonen tragere isomerisatiesnelheden dan mierenzuur, waarbij tunneling het dominante mechanisme wordt onder respectievelijk 100 K en 75 K.

• Dimeer Kinetiek: De studie onderzocht ook drie dimere soorten (cyclische en acyclische vormen).

  • Tunneling is verantwoordelijk voor meer dan 99% van de reactiesnelheidsconstante bij temperaturen $\le 150$ K voor de dimeren.
  • De acyclische dimeer tc4 $\to$ tt3 isomeriseert ongeveer tien keer sneller dan de cyclische dimeer tc1 $\to$ tt2 bij 10 K. Dit komt overeen met experimentele argonmatrixgegevens waar tc4 een levensduur heeft van 8,0 minuten versus 24 minuten voor tc1.
  • De cc5 $\to$ tc3 transformatie vertoont een lossere overgangstoestand en intermediaire kinetiek.
  • Dimeer reactiesnelheidsconstanten zijn over het algemeen lager dan die van hun monomeren en de seleniumanalogen zijn ordes van grootte langzamer dan de zuurstofsoorten.

• Matrix versus Gasfase: De auteurs merken op dat hun gasfase-resultaten een bovengrens vormen voor de reactiesnelheid. Experimentele snelheden in cryogene matrices zijn over het algemeen trager (door factoren van $10^3$ tot $10^4$ voor monomeren) vanwege matrixstabilisatie en "solvatie"-effecten, die de effectieve barrière verhogen of de tunnelingwaarschijnlijkheid verminderen.

Betekenis en Claims
De primaire bijdrage van dit werk is het leveren van rigoureuze gasfase kinetische gegevens die dienen als gids voor het interpreteren van experimentele matrixstudies en astrochemische modellen.

1. Resolutie van Discrepanties: De auteurs claimen dat hun berekende reactiesnelheidsconstanten voor mierenzuur ($k \approx 10^2 \text{ s}^{-1}$ bij 10 K) zeer verschillend zijn (door factoren die de $10^{22}$ overstijgen in sommige vergelijkingen) van eerdere gasfase theoretische studies die een verwaarloosbare reactiviteit voorspelden. Zij schrijven dit toe aan de inclusie van geavanceerde tunneling-correcties (SCT en QRC), die essentieel zijn bij lage temperaturen.
2. Astrochemische Implicaties: Gezien de berekende halfwaardetijd van cis-mierenzuur van minder dan 10 ms bij 10 K in de gasfase, stellen de auteurs dat cis-mierenzuur niet detecteerbaar zou moeten zijn in het interstellaire medium (ISM), tenzij het continu wordt geregenereerd door een mechanisme dat sneller is dan de isomerisatie. Zij suggeren dat photoswitching (UV-excitatie) of specifieke protonering/deprotoneringscycli op ijskorrels noodzakelijk kunnen zijn om de observationele aanwezigheid te verklaren.
3. Voorspellende Waarde: De studie concludeert dat experimentele onderzoeken naar de isomerisatie van zwavel- en seleniumanalogen in standaard cryogene matrices waarschijnlijk geen productieve resultaten zullen opleveren vanwege de extreem trage voorspelde snelheden, terwijl mierenzuur een levensvatbare kandidaat blijft voor dergelijke studies.

Het artikel stelt geen nieuwe experimentele opstellingen voor, maar biedt een theoretisch kader om bestaande observaties te rationaliseren en suggereert dat de snelle gasfase-isomerisatie alternatieve vormings- of regeneratiepaden voor cis-conformeren in de ruimte vereist.