# Skyrmion

Als Skyrmion (nach Tony Skyrme) wird in der theoretischen Physik ein Modell topologisch stabiler Solitonen-Wirbel in Feldern bezeichnet. Diese Wirbel verhalten sich wie Teilchen bzw. Quasiteilchen endlicher Masse.
Skyrmionen verwendete man als Modell ab 1958 bei Versuchen, die bis dahin rätselhafte Starke Wechselwirkung zu erklären. Ihr unterliegen insbesondere Protonen, Neutronen und Pionen. Tony Skyrme wollte die Starke Wechselwirkung damit erklären, dass Protonen und Neutronen Wirbel in Pionenfeldern wären. Man nannte die stabilen Wirbel "Skyrmionen". Um 1965 wurde klar, dass Protonen, Neutronen und Pionen aus Quarks bestehen. Damit bedurfte man in der Kernphysik keiner Skyrmionen als Erklärungsmodell mehr.
Ab den 1980er Jahren übernahm man den Modellbegriff in der Festkörperphysik, und er wurde auch in der Teilchenphysik mit Arbeiten von Edward Witten und verschiedenen Bag-Modellen für Hadronen populär (siehe auch Kenneth A. Johnson). In der Festkörperphysik wurde er u. a. beim Quanten-Hall-Effekt in zweidimensionalen Elektronengasen diskutiert. Derzeit untersucht man Skyrmionen auch an Oberflächen und Grenzflächen magnetischer Systeme.
Anfang 2009 konnte an der TU München von Sebastian Mühlbauer, Christian Pfleiderer, Peter Böni, dem Theoretiker Achim Rosch (Universität zu Köln) und anderen erstmals ein Skyrmionengitter in einem magnetischen Festkörper (Mangansilizium bei −245 °C und in einem Magnetfeld von 0,2 Tesla) direkt nachgewiesen werden. Eine im September 2010 eingereichte und im Juli 2011 veröffentlichte Publikation einer Forschergruppe der Universitäten Kiel und Hamburg sowie des Forschungszentrums Jülich beschreibt den ersten Nachweis von Skyrmionen ohne externes Magnetfeld. 2013 gelang es an der Universität Hamburg, Skyrmionen gezielt auf Oberflächen zu erzeugen und zu löschen. Da man stabile Skyrmionen auch bei Zimmertemperatur nachwies, erscheint ihr Einsatz in schnellen Informationsspeichern künftig möglich. Hierbei unterscheidet man in Kristallen Néel- und Bloch-Skyrmionen sowie Anti-Skyrmionen als Mischung aus Néel- und Bloch-Zuständen. 2019 gelang die dreidimensionale Auflösung der magnetischen Struktur von Skyrmionen, wobei die rund 100 nm großen Skyrmionen in Vielfachschichten von Ta/CoFeB/MgO untersucht wurden. Dabei stellte sich die Dipol-Dipol-Wechselwirkung zusammen mit der Wechselwirkung mit dem äußeren magnetischen Feld als besonders wichtig für die Stabilisierung heraus.

## Simulation in FDTD (2D)

Ein 2D-Skyrmion ergibt sich aus der Rotation des Solitons der Sinus-Gordon-Gleichung um den Mittelpunkt. Es kann in einem FDTD-Simulationsmedium nachgebildet werden, das in 2D einem Rechenkästchengitter ähnelt. In jedem Kästchen befindet sich auf der 2D-Fläche ein 3D-Vektorpfeil mit der Einheitslänge 1. Die dritte Koordinatenachse ermöglicht es, dass an einem Ort ohne Teilchen der Pfeil in eine der beiden neutralen Richtungen zeigt (z. B. Down). Oft wird das Skyrmion als ein auf die Fläche abgewickelter Igel gezeigt. Dabei kann man den Igel auch so realisieren, dass die Pfeile kästchenweise angeordnet sind. Im Schaubild sind meist alle Pfeile, die nach Down zeigen, weggelassen.
FDTD bedeutet, dass z. B. eine kleine Welt mit 360 × 240 Pixeln entsteht. Die wenigen in der Simulation verwendeten Formeln sind in dieser Welt Weltformeln. Der Autor der Simulation kann durch die Initialisierung des Mediums (Hauptspeicher) festlegen, was sich zur Laufzeit ereignen wird. Das Besondere der Sinus-Gordon-Gleichung ist, dass sie als Weltformel ruhende Partikel ermöglicht.
In den gezeigten Videos der Simulation wird ein Double-Buffer-Algorithmus verwendet. In Takt 1 erfolgt das Bilden von Differenzen der in den Pixeln gespeicherten Werte (1. Ableitung). Danach bildet ein im Double-Buffer gewöhnlich nicht vorkommender Schritt den Mittelwert von einem Pixel und seinen acht Nachbarpixeln. Als zweiter Schritt wird die Sinus-Gordon-Gleichung sowie ein Wellenalgorithmus angewendet, welche sich überlagern. Die Rechnung kommt ohne expliziten Aufruf der sin-Funktion aus, sie wird ja durch Simulation nachgebildet. Einfach gesagt ist jedes Pixel nur von seiner engen Umgebung aus dem Vorgängerbild abhängig. Es wird nicht abhängig von beliebiger Zeit gerechnet. Einschränkend ergibt sich die Notwendigkeit, wegen Schrittabweichungen durch Näherung die Einheitslänge der 3D-Vektoren bei jedem Schritt zu normalisieren. Dazu ist die Quadratwurzel notwendig. Auch enthält die Simulation einen Dämpfungsterm. Ränder sind bis jetzt nicht terminiert (siehe Tastkopf Terminierung). Durch die Nachbarschafts-Differenzbildung entsteht ein ungerades Gitter. Es überlappen sich Pixel aus Schritt eins und Schritt zwei. Zur Rechenzeitoptimierung zeigt die x-Achse 45° von links oben nach rechts unten. Die y-Achse zeigt 45° von links unten nach rechts oben. Von der Simulation wird im Zeitraffer nur jedes zweite Bild in das Video gerendert.
Die Variablen
Schritt 1 (für jedes Pixel in jedem Bild des Videos)
Schritt 2 (für jedes Pixel in jedem Bild des Videos)