Ein Zwei-Spin-System ist ein nüchterner Prüfstein für NMR-Simulationen: vier Zustände, eine 4×4-Dichtematrix und Relaxationsraten, die von Temperatur, Feldstärke und Umgebung abhängen. Wer hier reproduzierbare Werte braucht, kommt mit einer einzelnen Formel nicht weit. Genau an dieser Stelle wird der Rechenweg entscheidend.
Das Problem: Kleine Systeme, große Fehlerquellen
Bei zwei Spins wirkt die Aufgabe überschaubar. Doch schon die Unterscheidung zwischen T1- und T2-Relaxation verlangt, dass Populations- und Kohärenzterme getrennt behandelt werden. Ein vertauschtes Vorzeichen in der Relaxationsmatrix fällt in einer breiten Linie kaum auf, verschiebt aber das gesamte Relaxationsprofil. Für Lehre und Methodenvergleich ist das ein handfester Fehler.
Ein typischer Praxisvergleich zeigt die Folge: Zwei Arbeitsgruppen simulieren dasselbe Zwei-Spin-System, nutzen aber unterschiedliche Reihenfolgen für Operatoren. Die Kurven ähneln sich, die Relaxationszeiten weichen jedoch um mehrere Prozent ab. Ohne dokumentierte Zwischenschritte lässt sich nicht klären, ob die Abweichung aus der Physik oder aus der Konvention stammt.
Die Ursache: Dichtematrizen koppeln viele Beiträge
Die Dichtematrix beschreibt nicht nur Besetzungen, sondern auch Phasenbeziehungen zwischen Zuständen. Im Zwei-Spin-Fall ergeben sich 16 Matrixelemente, von denen mehrere über Relaxationskanäle miteinander verknüpft sind. Wer alle Terme von Hand sortiert, verliert leicht den Überblick. Die Folge sind uneinheitliche Konventionen, doppelt gezählte Raten oder eine nicht-physikalische Spur. Genau darunter leidet die Vergleichbarkeit von Simulationsergebnissen.
Erschwerend kommt hinzu, dass Relaxation kein einzelner Term ist. Populationen und Kohärenzen unterliegen unterschiedlichen Dämpfungen; Kreuzrelaxation koppelt selbst im kleinen System mehrere Matrixelemente. Wird das nicht getrennt erfasst, entsteht eine Matrix, die formal rechnet, aber physikalisch nicht mehr sauber interpretierbar ist.
Die Lösung: spinformula als strukturierte Rechenhilfe
Hier setzt spinformula an: Der Ansatz zwingt dazu, Zustände, Operatoren und Relaxationsbeiträge in einer festen Reihenfolge zu erfassen. Statt Formeln frei zu kombinieren, wird der Rechenweg in prüfbaren Abschnitten aufgebaut. Das macht die Simulation nicht automatisch richtig, aber es macht Fehler sichtbar. Besonders im Minimalbeispiel mit zwei Spins lässt sich jeder Schritt nachvollziehen.
Für den Einstieg reicht ein festes Schema: erst Basis, dann Operatoren, dann Relaxationskanäle, zuletzt Kontrollen. Dieses Schema ist bewusst schmal gehalten. Es ersetzt keine vollständige Relaxationstheorie, verhindert aber, dass grundlegende Konventionen unbemerkt durcheinanderlaufen.
Schritt für Schritt zum Relaxationsprofil
- Zustandsbasis festlegen: Die vier Zustände |↑↑⟩, |↑↓⟩, |↓↑⟩ und |↓↓⟩ werden eindeutig geordnet. Ohne diese Festlegung ist jede spätere Matrix interpretationslos.
- Operatoren zuordnen: Für jeden Spin werden Leiteroperatoren und z-Komponenten eingetragen. Die Reihenfolge entscheidet über Vorzeichen und Kopplungsterme.
- Relaxationsbeiträge eintragen: T1-Prozesse ändern Besetzungen, T2-Prozesse dämpfen Kohärenzen. Beide Kanäle werden getrennt notiert, bevor sie in die Matrix übergehen.
- Spur und Symmetrie prüfen: Die Dichtematrix muss Spur eins behalten, ihre Hermitizität darf nicht verloren gehen. Diese Kontrollen sind kein Formalismus, sondern ein früher Filter gegen Rechenfehler.
Was der Minimaltest zeigt
Das Zwei-Spin-Beispiel liefert kein vollständiges Modell realer Proben, aber einen belastbaren Maßstab. Wer die vier Schritte sauber durchläuft, erkennt, ob die Relaxationsmatrix konsistent aufgebaut ist. Im Vergleich zweier Parametersätze bleibt dann nachvollziehbar, welche Änderung wirklich aus der Physik kommt und welche aus der Rechenstruktur. Für die Lehre ist das wertvoller als eine glatte Kurve ohne Herleitung.