Perspective on a challenge: predicting the photochemistry of cyclobutanone

Deze perspectiefartikel vat samen hoe een voorspellingsuitdaging voor de fotochemie van cyclobutanon, waarbij 15 theoriegroepen hun simulaties van tijd-opgeloste MeV-UED-signalen vergeleken met experimentele data, de kwalitatieve voorspellingskracht van niet-adiabatische moleculaire dynamica aantoont en het belang van zorgvuldige benchmarking van elektronische-structuurtheorie benadrukt.

De kern

De Grote Voorspellingswedstrijd: Hoe Cyclobutanon Zich Gedraagt onder een Lasmicroscoop

Stel je voor dat je een heel klein, kwetsbaar huisje hebt: cyclobutanon. Dit is een molecuul, een bouwsteen van de chemie. Nu, stel je voor dat je dit huisje plotseling raakt met een flits van blauw licht (een laserpuls). Wat gebeurt er dan? Verandert het huisje in een ander huisje? Breekt het uit elkaar? Of springt het in een vreemde dans?

In de wereld van de wetenschap willen we dit kunnen voorspellen voordat we het in het echt zien. Maar dit is lastig. Het is alsof je probeert te voorspellen hoe een dansgroep zich beweegt, terwijl je alleen een foto van het begin en een foto van het einde hebt, en je moet de hele dans in het midden uitvinden.

Dit artikel vertelt het verhaal van een enorme wetenschappelijke uitdaging die in 2023 werd georganiseerd.

De Wedstrijd: Een Voorspel-Feestje

De organisatoren (een groep van meer dan 70 wetenschappers) riepen de hele wereld van computerchemici op: "Kunnen jullie precies voorspellen wat er gebeurt met dit molecuul als we het raken met licht, en hoe het eruitziet als we het fotograferen met een supersnelle camera?"

Ze gaven 15 verschillende teams de opdracht om dit te simuleren. Elke groep gebruikte zijn eigen 'rekenmethode' en 'rekenmachine'. Het was als een kookwedstrijd waarbij iedereen een eigen recept had om hetzelfde gerecht te maken, maar niemand wist hoe het er echt uit zou zien.

De Drie Grote Uitdagingen (De 'Recepten')

Om de dans van het molecuul te voorspellen, moesten de teams drie moeilijke stappen nemen:

1. De Startpositie (De Foto voor de Dans)
Voordat het licht erop schijnt, moet je weten hoe het molecuul er precies uitziet. Is het perfect plat? Of staat het een beetje scheef?

• De Analogie: Stel je voor dat je een groep mensen vraagt om te dansen. Sommige teams dachten: "Ze staan allemaal perfect rechtop." Andere teams dachten: "Nee, ze wankelen een beetje en staan in een kromme houding."
• Het Resultaat: Het bleek dat deze kleine verschillen in de startpositie een groot effect hadden op hoe snel de dans begon.

2. De Rekenmachine (De Regels van de Dans)
Dit is het hart van de zaak. Hoe bereken je de krachten die het molecuul bewegen?

• De Simpele Regels (Single-Reference): Sommige teams gebruikten simpele regels. Alsof je zegt: "Het molecuul is een balletje dat rolt." Dit werkt goed als het balletje rustig is, maar faalt als het begint te springen en te veranderen.
• De Complexe Regels (Multireference): Andere teams gebruikten super-complexe regels. Alsof je zegt: "Het molecuul is een danser die tegelijkertijd in meerdere poses kan zijn." Dit is veel preciezer, maar ook veel moeilijker te berekenen.
• Het Grote Geheim: De teams die de 'complexe regels' gebruikten, voorspelden dat het molecuul snel uit elkaar viel. De teams met de 'simpele regels' dachten dat het molecuul langzaam bleef hangen. De echte experimenten (die later werden gedaan) bewezen dat de complexe regels het beste werkten.

3. De Camera (Het Filmen van de Dans)
De teams moesten niet alleen de dans voorspellen, maar ook voorspellen wat een supersnelle camera (MeV-UED) zou zien. Deze camera maakt 'moleculaire films' van atomen.

• De Analogie: Het is alsof je probeert te voorspellen hoe een danser eruitziet als je hem filmt met een wazige camera. Sommige teams vergeten dat de camera ook 'schaduwen' (elektronen) ziet, en niet alleen de mensen (atomen).
• Het Resultaat: De beste voorspellingen leken het meest op de echte film die later werd gemaakt.

Wat Leerden We? (De Lessen van de Wedstrijd)

Toen alle teams hun resultaten hadden ingeleverd en de echte experimenten werden gedaan, konden ze alles vergelijken. Hier zijn de belangrijkste lessen, vertaald naar alledaags taal:

• De Rekenmachine is Cruciaal: Als je de verkeerde regels gebruikt (te simpel), krijg je een verkeerd verhaal. Je moet de complexiteit van de moleculaire dans serieus nemen.
• Geen 'Zwarte Doos': Je kunt niet zomaar een computerprogramma aanklikken en hopen dat het goed werkt. Je moet precies weten welke instellingen je gebruikt. Twee teams met bijna dezelfde software kregen soms heel verschillende resultaten omdat ze één klein knopje anders hadden gedraaid.
• Samenwerking is Goud: Als één team dit had geprobeerd, hadden ze waarschijnlijk niet de juiste conclusie getrokken. Door 15 teams te laten concurreren en samen te werken, konden ze zien wat er goed ging en wat niet. Het was een enorme 'kalibratie-oefening' voor de hele wetenschap.
• We zijn bijna klaar: De vraag was: "Kunnen we dit nu al voorspellen?" Het antwoord is: Bijna. We kunnen het verhaal van de dans goed vertellen (kwalitatief), maar de exacte timing en snelheid (kwantitatief) moeten nog iets scherper.

Conclusie: De Dans Gaat Door

Dit artikel is een feest van samenwerking. Het laat zien dat wetenschappers uit de hele wereld (van Nederland tot China, van Amerika tot Zwitserland) samenwerken om de geheimen van het heelal te ontrafelen.

Ze hebben bewezen dat we met computers steeds beter kunnen voorspellen hoe moleculen reageren op licht. Dit is belangrijk voor de toekomst, bijvoorbeeld om nieuwe medicijnen te maken of betere zonnepanelen te ontwerpen. We zijn nog niet op het punt waar we elke dans perfect kunnen voorspellen, maar we hebben een enorme stap gezet in de goede richting.

Kortom: We hebben de dansvoorspelling bijna onder de knie, maar we moeten nog even oefenen met de juiste muziek (de rekenregels).

Technische samenvatting

Titel: Perspectief op een uitdaging: het voorspellen van de fotochemie van cyclobutanon

1. Het Probleem en de Context
Het veld van de niet-adiabatische moleculaire dynamica (NAMD) heeft de afgelopen decennia grote vooruitgang geboekt, maar de vraag of het veld nu voldoende volwassen is om voorspellend te werken (in plaats van alleen verklarend), bleef open. Om dit te testen, werd in 2023 een "voorspellingsuitdaging" georganiseerd voor de computationele fotochemie-gemeenschap.

• Doel: Voorspellen van de fotochemie van gasfase-cyclobutanon na excitatie met een 200 nm laserpuls, en het voorspellen van het bijbehorende tijdsopgeloste MeV-UED (Mega-electronvolt Ultrafast Electron Diffraction) signaal.
• Uitdaging: De experimentele resultaten waren op het moment van de uitdaging nog niet bekend (het experiment vond plaats eind januari 2024 bij SLAC).
• Deelnemers: 15 theoretische voorspellingen van meer dan 70 onderzoekers, gebruikmakend van een breed scala aan strategieën voor elektronische structuur en dynamica.

2. Methodologie
De bijdragen werden geanalyseerd op vier cruciale componenten die een NAMD-simulatie definiëren:

• A. Fotoexcitatie (Initiële Voorwaarden):

  • De meeste deelnemers gebruikten een harmonische Wigner-verdeling om de initiële atoomposities en -momenta te bemonsteren vanuit de grondtoestand.
  • Er waren variaties in hoe de ring-plooiing (ring-puckering) van cyclobutanon werd behandeld (soms als harmonisch, soms met kwantumcorrecties of anharmonische benaderingen).
  • Enkele groepen pasten een "energy windowing" toe om de bandbreedte van de laserpuls (200 nm) na te bootsen, terwijl anderen dit niet deden.

• B. Elektronische Structuur:

  • Er werd gebruik gemaakt van een breed scala aan methoden, variërend van enkel-referentie methoden (LR-TDDFT, EOM-CCSD, MP2) tot multiconfiguratie- en multireferentie methoden (SA-CASSCF, XMS-CASPT2, MRCI).
  • Een kritisch punt was het beschrijven van de overgang van de Rydberg-toestand (S2, 3s-karakter) naar een valentietoestand en de daaropvolgende ringopening.
  • Verschillende basissets (van cc-pVDZ tot aug-cc-pVDZ) en actieve ruimtes werden gebruikt.

• C. Niet-Adiabatische Dynamica:

  • De meest gebruikte methode was "Trajectory Surface Hopping" (FSSH), vaak met decoherentiecorrecties.
  • Andere methoden omvatten Ab Initio Multiple Spawning (AIMS), Ab Initio Multiple Cloning (AIMC), en DD-vMCG (Direct Dynamics variational Multiconfigurational Gaussian).
  • Twee groepen gebruikten MCTDH (Multiconfiguration Time-Dependent Hartree) met modelpotentiaals.
  • De simulatietijden varieerden sterk, van enkele honderden femtoseconden tot meerdere picoseconden.

• D. Berekening van Observabelen:

  • De theoretische signalen werden omgezet naar MeV-UED signalen ($\Delta I/I$ en $\Delta rPDF$) met behulp van het "Independent Atom Model" (IAM).
  • De resultaten werden vergeleken met twee experimentele datasets: Exp1 (SLAC, tot 2 ps) en Exp2 (Shanghai Jiao Tong University, tot 1.2 ps).

3. Belangrijkste Bijdragen en Resultaten

• Rol van de Elektronische Structuur: De keuze van de elektronische structuurmethode bleek de meest kritische factor.
  • SA-CASSCF: Voorspelde een zeer snelle verlaten van de Franck-Condon regio (binnen tientallen femtoseconden) en een snelle ringopening. Dit kwam niet overeen met de experimentele tijdschalen.
  • XMS-CASPT2: Voorspelde dat het golfpakketje enkele honderden femtoseconden in de S2(3s) Rydberg-minimum blijft, gevangen door een energiebarrière, voordat het adiabatisch evolueert naar een valentietoestand en ringopent. Dit kwam het beste overeen met de experimentele levensduur (~0.23 - 0.29 ps).
  • Enkel-referentie methoden (TDDFT/EOM-CCSD): Voorspelden vaak dat het molecuul te lang in de Rydberg-toestand bleef (picoseconden) of vertoonden inconsistenties bij het beschrijven van de barrière voor ringopening.
• Dynamica en Producten:
  • De consensus is dat intersysteemkruising (ISC) een verwaarloosbare rol speelt op de sub-picoseconde tot picoseconde tijdschalen.
  • De belangrijkste fotoproducten zijn CO + cyclopropaan/propaan (C3-kanaal) en ketene + ethyleen (C2-kanaal). De verhouding tussen deze kanalen varieerde tussen de voorspellingen, maar de algemene chemische uitkomst was consistent.
  • Er werd een barrière gevonden in de S2-toestand die de snelheid van de ringopening bepaalt. Een kleine fout in de berekening van deze barrière (enkele kJ/mol) leidt tot een grote fout in de voorspelde kinetiek (orde van grootte verschil).
• Observabelen:
  • De meeste voorspellingen misten het experimentele signaal bij lage impuls-overdracht ($s < 1.2 \text{ \AA}^{-1}$), dat wordt toegeschreven aan inelastische verstrooiing door het Rydberg-karakter. De meeste simulaties gebruikten alleen het elastische IAM-model.
  • De voorspellingen gebaseerd op XMS-CASPT2 toonden de beste overeenkomst met de experimentele tijdschalen voor de vorming van fotoproducten in het $\Delta rPDF$ signaal.

4. Significatie en Conclusies

• Kwalitatieve Voorspellingskracht: De uitdaging toonde aan dat NAMD kwalitatief voorspellend is. De meeste methoden voorspelden de juiste familie van fotoproducten en het gebrek aan intersysteemkruising.
• Kritieke Afhankelijkheid van Elektronische Structuur: De studie benadrukt dat de nauwkeurigheid van de dynamica direct afhankelijk is van de kwaliteit van de elektronische structuurtheorie, vooral in gebieden waar statische en dynamische correlatie beide belangrijk zijn (zoals bij het doorbreken van bindingen en het passeren van conische doorsnijdingen).
• Calibratie-oefening: Dit project diende als een "calibratie-oefening" voor de gemeenschap. Het onthulde dat kleine verschillen in implementatie (bijv. hoe niet-adiabatische koppelingen worden berekend of hoe decoherentie wordt behandeld) grote verschillen in resultaten kunnen veroorzaken.
• Toekomstperspectief:
  • Er is een dringende behoefte aan betere benchmarking van elektronische structuurmethoden buiten de Franck-Condon regio.
  • De gemeenschap moet meer samenwerken met ontwikkelaars van elektronische structuurcodes om methoden te testen op "pathologische" gevallen.
  • Toekomstige uitdagingen zouden moeten focussen op het isoleren van dynamische methoden (met vaste elektronische structuur) en het integreren van inelastische verstrooiingseffecten in de observabelen.

Samenvattend: Hoewel het veld nog niet volledig kwantitatief voorspellend is voor alle aspecten (vooral vanwege de gevoeligheid voor de elektronische structuur en de barrièrehoogten), heeft deze uitdaging bewezen dat NAMD een krachtig instrument is om complexe fotochemische mechanismen te ontrafelen, mits de juiste methodologische keuzes worden gemaakt en zorgvuldig worden gevalideerd.