Extending numerical simulations in SIMPSON: Electron paramagnetic resonance, dynamic nuclear polarisation, propagator splitting, pulse transients, and quadrupolar cross terms

Dit paper introduceert een nieuwe C++-versie van de SIMPSON-software met geavanceerde functies voor het simuleren en interpreteren van hybride NMR/EPR-experimenten, waaronder dynamische nucleaire polarisatie, propagatorsplitsing en kwadrupolaire kruistermen.

De kern

SIMPSON 6.0: De Super-Update voor het Spioneren op Atomen

Stel je voor dat atomen en elektronen een enorme, ingewikkelde danszaal zijn. Wetenschappers proberen deze dans te begrijpen door er met speciale microfoons (magnetische resonantie) naar te luisteren en met flitslichten (pulsen) de dansers te sturen. De software SIMPSON is de beroemde "rekenmachine" die wetenschappers gebruiken om te voorspellen hoe deze dansers zich zullen gedragen voordat ze het echt in het lab doen.

Deze nieuwe paper introduceert SIMPSON versie 6.0. Het is alsof ze de rekenmachine niet alleen hebben opgefrist, maar er een compleet nieuw brein, snellere benen en een nieuwe taal voor hebben ingebouwd. Hier is wat er nieuw is, vertaald naar alledaagse taal:

1. Van C naar C++: De Motor wordt Vervangen

Vroeger was SIMPSON geschreven in een wat ouderwetse programmeertaal (C). Nu is het volledig herschreven in C++.

• De analogie: Denk aan een oude, betrouwbare Volkswagen Kever. Die doet het prima, maar als je er een Formule-1-auto van wilt maken, moet je de motor, het chassis en de elektronica vervangen. C++ is die nieuwe Formule-1-motor. Het maakt de software sneller, flexibeler en makkelijker voor andere programmeurs om er aan te sleutelen (net als Lego-blokken die makkelijker te koppelen zijn).

2. Elektronen en Kernen: Een Nieuwe Danspartner

Tot nu toe kon SIMPSON vooral goed kijken naar atoomkernen (NMR). Maar nu kunnen ze ook elektronen (EPR) en de interactie tussen beide (DNP) simuleren.

• De analogie: Stel je voor dat je vroeger alleen kon kijken naar hoe een zware olifant (de atoomkern) dansde. Nu kun je ook kijken naar een snelle, wilde muis (het elektron) die op de olifant springt. Soms helpen ze elkaar: de muis geeft een duwtje aan de olifant, waardoor de olifant veel harder danst. Dit helpt wetenschappers om heel zwakke signalen (zoals van een heel klein monster) toch helder te horen.

3. Propagator Splitting: De Kunst van het Opdelen

Een groot probleem bij het simuleren van complexe dansen is dat de berekeningen enorm langzaam zijn. De software moet elke seconde van de dans in miljoenen stukjes opbreken.

• De analogie: Stel je voor dat je een gigantische pizza moet eten. Als je probeert de hele pizza in één hap te nemen, ben je er nooit klaar mee. De nieuwe methode, Propagator Splitting, zegt: "Laten we de pizza eerst in grote plakken snijden, en die plakken pas in hapjes."
  • Het berekent eerst de simpele bewegingen apart, en dan pas de complexe interacties.
  • Resultaat: De berekeningen gaan veel sneller, zonder dat de pizza (de nauwkeurigheid) minder lekker wordt.

4. Pulse Transients: De "Trage" Flits

Wanneer je een flitslicht (een puls) aanzet, duurt het even voordat het echt fel brandt. Het licht flikkert eerst even. In de echte wereld is dit een probleem, maar in de oude software werd dit vaak genegeerd.

• De analogie: Het is alsof je een auto start. Je draait de sleutel, maar de motor hoest even voordat hij goed loopt. De nieuwe SIMPSON houdt hier rekening mee. Hij zegt: "Oké, de flits is niet perfect, hij heeft een 'hik'."
  • De software kan nu zelfs de pulse corrigeren. Het berekent precies hoe de flits eruit moet zien om die "hik" te compenseren, zodat de atomen toch perfect dansen, zelfs als de apparatuur niet perfect is.

5. Optimal Control: De Dansmeester

Soms willen wetenschappers een heel specifieke dansbeweging maken die niet in de standaardhandleiding staat. Ze gebruiken "Optimal Control" om een perfecte dans te ontwerpen.

• De analogie: Stel je voor dat je een dansmeester bent die een choreografie bedenkt voor een groep dansers die allemaal net iets anders reageren op muziek (sommige zijn traag, sommige snel). De nieuwe software helpt de dansmeester om een choreografie te vinden die voor iedereen perfect werkt, zelfs als de muziek soms haperend is of de dansers moe zijn.

6. Vierkante Kernen: De Moeilijke Dansers

Sommige atoomkernen (zoals die van zuurstof of stikstof) zijn "vierkant" in plaats van rond. Ze zijn veel moeilijker te simuleren.

• De analogie: De meeste atomen zijn als balletdansers: rond en soepel. Deze "vierkante" atomen zijn als dansers met zware blokken aan hun voeten. De nieuwe SIMPSON kan nu berekenen hoe die blokken de dans beïnvloeden, zelfs als ze heel snel draaien. Dit is cruciaal voor het bestuderen van nieuwe batterijen en materialen.

Waarom is dit belangrijk?

Vroeger moesten wetenschappers vaak wachten tot hun computer klaar was met rekenen, of ze moesten het doen met een "ruwe schatting". Met deze update:

1. Het gaat sneller: Berekeningen die uren duurden, gaan nu in minuten.
2. Het is accurater: Ze kunnen nu kijken naar de kleine details (zoals de "hik" van de flits) die eerder werden gemist.
3. Het is breder: Ze kunnen nu onderzoek doen naar nieuwe gebieden zoals quantumcomputers en superkrachtige batterijen, omdat ze zowel elektronen als kernen in één simulatie kunnen laten dansen.

Kortom: SIMPSON 6.0 is de nieuwe, superkrachtige simulator die wetenschappers helpt om de dans van de atomen niet alleen te voorspellen, maar ook te perfectioneren, zelfs als de apparatuur niet perfect is.

Technische samenvatting

Titel: Uitbreiding van numerieke simulaties in SIMPSON: EPR, DNP, propagator-splitsing, puls-transiënten en kruistermen van quadrupolaire kernen.

Auteurs: David L. Goodwin et al. (Aarhus University en Charles University in Praag)
Datum: Februari 2026 (voorgesteld)

---

1. Het Probleem

Numerieke simulaties zijn onmisbaar voor het ontwikkelen, begrijpen en interpreteren van geavanceerde experimenten in magnetische resonantie (NMR, EPR en DNP). Hoewel er al diverse softwarepakketten bestaan (zoals Spinach, EasySpin, SPINEVOLUTION), is de software SIMPSON, een standaard voor vaste-stof NMR, sinds de laatste grote update (versie 4.0) meer dan een decennium geleden niet significant uitgebreid.

De huidige uitdagingen zijn:

• Complexiteit van experimenten: De opkomst van Dynamic Nuclear Polarisation (DNP) en Pulsed EPR vereist simulaties die zowel elektronen- als kernspinsystemen koppelen, wat rekenkundig zeer intensief is.
• Rekenkundige bottlenecks: Het berekenen van de tijdsontwikkeling van de dichtheidsmatrix (via matrix-exponentiële functies) wordt prohibitief duur naarmate de dimensie van de Hilbert-ruimte toeneemt (bijv. bij grote spin-systemen of optimale controle).
• Experimentele realiteit: Simulaties moeten rekening houden met niet-ideale pulsen, zoals puls-transiënten (vervorming door de resonator) en ruimtelijke inhomogeniteit van het RF/MW-veld.
• Hoge-orde effecten: Voor quadrupolaire kernen (spin > 1/2) zijn tweede-orde kruistermen (bijv. tussen quadrupolaire koppeling en CSA of dipolaire koppeling) essentieel voor nauwkeurige spectraalanalyse.
• Onderhoud en extensie: De oude code was geschreven in C, wat community-bijdragen en verdere ontwikkeling bemoeilijkte.

2. Methodologie

De auteurs introduceren SIMPSON versie 6.0, een fundamentele herziening van de software met de volgende methodologische pijlers:

• Programmeertaal en Architectuur:

  • Migratie van C naar C++ met object-georiënteerde principes. Dit verbetert de snelheid, leesbaarheid en maakt de code toegankelijker voor community-contributions.
  • De besturingsscripttaal blijft TCL om compatibiliteit met bestaande inputbestanden te behouden, maar de kern is herschreven.
  • Beschikbaarheid via Docker, NMRbox en GitLab voor broncode-gebaseerde compilatie.

• Nieuwe Fysieke Modellen:

  • EPR en DNP: Implementatie van elektronenspin-dynamica, inclusief Zeeman-interactie, hyperfijne koppeling, elektron-elektron dipolaire koppeling en Heisenberg-uitwisseling.
  • Frame-selectie: Introductie van de DNPframe methode, waarbij elektronen in een roterend frame worden berekend en kernen in het laboratoriumframe.

• Propagator-Splitsing (Propagator Splitting):

  • Om de kosten van matrix-exponentiële functies te verminderen, wordt de totale Hamiltoniaan opgesplitst in een deel $A$ (controle-pulsen en offset, analytisch oplosbaar) en een deel $B$ (interacties).
  • Toepassing van Lie-Trotter en Strang-splitsing (en hogere-orde varianten tot orde 6) om tijdspropagatie te benaderen zonder volledige diagonalisatie. Dit versnelt berekeningen aanzienlijk, vooral bij optimale controle.

• Puls-Transiënten en Veldinhomogeniteit:

  • Implementatie van een vervormingsoperator ($\phi_{nm}$) die de convolutie van een ideale puls met een impulsresponsfunctie ($h(t)$) simuleert.
  • Voor optimale controle wordt de gradiënt van de fideliteit via de kettingregel teruggetransformeerd van het gedetailleerde (vervormde) tijdsnet naar het ruwe (rechthoekige) inputnet.
  • Ondersteuning voor ruimtelijke RF-veldverdelingen (bijv. door MAS-rotatie geïnduceerde modulaties).

• Quadrupolaire Kruistermen:

  • Toevoeging van Hamiltonian-termen voor tweede-orde kruistermen tussen quadrupolaire koppeling en chemische verschuiving (CSA) of dipolaire koppeling. Dit is cruciaal voor 2D STMAS-experimenten en residuële dipolaire splitsing.

3. Belangrijkste Bijdragen

1. SIMPSON v6.0 Release: Een volledig herschreven C++-versie die EPR, DNP en geavanceerde NMR in één pakket integreert.
2. Efficiëntie door Splitsing: De introductie van propagator-splitsing ($prop\_split$) reduceert de rekentijd voor complexe spin-systemen en optimale controle met een orde van grootte, zonder significant verlies aan nauwkeurigheid (afhankelijk van de gekozen splitsingsorde).
3. Realistische Pulsmodellering: Mogelijkheid om pulsen te optimaliseren en te simuleren rekening houdend met instrumentele transiënten en veldinhomogeniteit, wat essentieel is voor de overdracht van simulatie naar experiment.
4. Geavanceerde Kwantificering: Ondersteuning voor tweede-orde kruistermen bij quadrupolaire kernen, wat de interpretatie van complexe vaste-stof spectra mogelijk maakt.
5. Community en Visualisatie: Integratie met EasyNMR (web-based workflow) en SimView (Python GUI), en verbeterde documentatie voor community-ontwikkeling via GitLab.

4. Resultaten

De auteurs demonstreren de nieuwe functionaliteit via diverse voorbeelden:

• ESEEM en DNP Simulaties: Succesvolle simulatie van Electron Spin Echo Envelope Modulation (ESEEM) en diverse pulsed DNP-sequenties (NOVEL, BEAM, PLATO, cRW-OPT1). De software kan breedbandige excitatie en polarisatietransfer efficiënt modelleren.
• Optimale Controle (Optimal Control):
  • Een voorbeeld van een "composite inversion pulse" die RF-inhomogeniteit compenseert, geoptimaliseerd met de simplex-methode.
  • SORDOR-pulsen: Optimalisatie van fase-only pulsen met hoge fideliteit, waarbij de gebruikte propagator-splitsing (orde 2 en 3) de rekentijd drastisch verkortte vergeleken met exacte diagonalisatie (diag).
  • Resultaten tonen aan dat splitsing van orde 2 of 3 vaak voldoende is voor ensemble-optimatie en tot 50-80% tijdswinst oplevert.
• Transiënt-compensatie: Simulaties tonen aan dat optimalisatie die rekening houdt met puls-transiënten leidt tot robuustere pulsen die beter presteren in de praktijk dan die geoptimaliseerd voor ideale pulsen.
• Quadrupolaire Effecten: 2D STMAS-spectra tonen de splitsing van pieken veroorzaakt door kruistermen (quadrupool-CSA en quadrupool-dipool), wat overeenkomt met theoretische verwachtingen. De gebruikte splitsingsmethode bleek ook hier aanzienlijk sneller dan de standaard diagonalisatie, vooral bij grotere matrixgroottes (bijv. 16x16).

5. Significantie

Deze publicatie markeert een mijlpaal in de numerieke simulatie van magnetische resonantie:

• Interdisciplinaire Integratie: Het overbrugt de kloof tussen NMR, EPR en DNP, wat essentieel is voor het groeiende veld van hybride kwantumsensoren en materialenonderzoek.
• Rekenkundige Efficiëntie: De implementatie van propagator-splitsing maakt het mogelijk om grotere en complexere spin-systemen te simuleren binnen redelijke tijdsbestekken, wat eerder onmogelijk was met traditionele methoden.
• Experimentele Relevantie: Door rekening te houden met transiënten en veldinhomogeniteit, sluit de software beter aan bij de realiteit van moderne experimenten, wat de interpretatie van data en het ontwerp van nieuwe pulsen verbetert.
• Toekomstbestendigheid: De overstap naar C++ en de focus op community-bijdragen zorgen voor de levensvatbaarheid en verdere ontwikkeling van SIMPSON als open-source standaard in de wetenschappelijke gemeenschap.

Kortom, SIMPSON v6.0 biedt een krachtig, snel en veelzijdig platform voor de volgende generatie magnetische resonantie-experimenten, met name in de context van dynamische nucleaire polarisatie en geavanceerde kwantumsimulaties.