Pressure-Temperature Phase Diagram and $\lambda$-Transition in Liquid Sulfur

Mithilfe von Molekulardynamik-Simulationen liefert diese Arbeit eine mikroskopische Erklärung für den $\lambda$-Übergang in Schwefel und zeigt auf, wie die thermisch induzierte Polymerisation in Abhängigkeit von Druck und Temperatur verläuft.

Das Wesentliche

Das Geheimnis des „tanzenden“ Schwefels: Eine Geschichte von Ringen und Ketten

Stellen Sie sich vor, Sie beobachten eine riesige, perfekt organisierte Tanzveranstaltung in einer Turnhalle. Das ist unser kristalliner Schwefel bei niedrigen Temperaturen. Alle Tänzer (die Schwefelatome) sind in kleinen, geschlossenen Gruppen von acht Personen unterwegs. Sie bilden perfekte, geschlossene Kreise – wie kleine, unbewegliche Ringe, die ordentlich in Reihen aufgestellt sind. Alles ist ruhig, geordnet und sehr vorhersehbar.

1. Der große Umbruch (Der Schmelzvorgang)

Wenn wir die Temperatur erhöhen, wird die Musik lauter und schneller. Die Tänzer werden unruhig. Irgendwann wird es so heiß, dass die Ordnung zusammenbricht: Die Reihen lösen sich auf, die Ringe fangen an zu rollen. Das ist das Schmelzen. Jetzt haben wir eine „molekulare Flüssigkeit“ – eine Menge kleiner, gelber Ringe, die wild durcheinanderrollen.

2. Die „Lambda-Transition“: Wenn aus Kreisen Schlangen werden

Jetzt kommt der spannendste Teil, den die Forscher die „$\lambda$-Transition“ nennen. Stellen Sie sich vor, die Musik wird plötzlich extrem intensiv und fast schon aggressiv.

Anstatt dass die Tänzer einfach nur in ihren Achtertruppen weiterrollen, passiert etwas Seltsames: Einige der Ringe „platzen“ auf. Aber sie bleiben nicht als Einzelteile liegen. Stattdessend fangen sie an, sich an die Hände der Nachbarn zu greifen. Aus den kleinen, geschlossenen Kreisen werden plötzlich riesige, endlose Schlangenlinien (Polymerketten).

Das ist wie bei einem riesigen „Stopptanz“, bei dem plötzlich alle Teilnehmer versuchen, sich zu einer einzigen, gigantischen Schlange zu verbinden, die durch die ganze Halle zieht. Die Flüssigkeit verändert sich schlagartig: Sie wird zäher, dunkler und verhält sich völlig anders.

3. Was haben die Forscher neu entdeckt? (Die „Zündkerzen“)

Früher dachte man, die Ringe würden einfach so zufällig aufbrechen. Die Forscher haben aber mit einem extrem leistungsstarken „Super-Mikroskop“ (einer KI-gestützten Computersimulation) etwas Neues gesehen:

• Die „unperfekten“ Ringe als Zündkerzen: Es sind nicht die perfekten 8er-Ringe, die das Chaos auslösen. Es sind die „unperfekten“ Ringe (z. B. mit 7 oder 9 Atomen). Diese kleinen Gruppen sind instabil – sie sind wie kleine, nervöse Tänzer, die zuerst die Ordnung stören und so die Kettenbildung überhaupt erst „anzünden“.
• Das Chaos im Festkörper: Wenn man den Schwefel mit extrem hohem Druck presst, passiert etwas Erstaunliches. Die Kettenbildung beginnt nicht erst, wenn der Schwefel flüssig ist, sondern sie fängt schon innerhalb des festen Kristalls an! Es ist, als würden die Tänzer schon anfangen, sich zu Schlangen zu verbinden, während sie noch stramm in ihren Reihen stehen. Erst wenn die Ketten groß genug sind, bricht die ganze Struktur zusammen und wird flüssig.

Zusammenfassend: Was bedeutet das?

Die Forscher haben eine Art „Landkarte“ (das Phasendiagramm) erstellt, die zeigt, wie Schwefel unter verschiedenen Temperaturen und Drücken reagiert. Sie haben verstanden, dass Schwefel nicht einfach nur von „fest“ zu „flüssig“ wechselt, sondern eine faszinierende Verwandlung von „kleinen Kreisen“ zu „endlosen Ketten“ durchmacht.

Die Metapher zum Mitnehmen:
Schwefel ist wie eine Party, die von einem geordneten Kreisspiel (Kristall) über ein wildes Rollen von Reifen (molekulare Flüssigkeit) bis hin zu einem riesigen, zähen Schlangen Tanz (Polymer-Flüssigkeit) führt – und der Druck bestimmt, ob die Party erst im Sitzen oder erst beim Tanzen losgeht!

Technische Zusammenfassung

Zusammenfassung: Druck-Temperatur-Phasendiagramm und $\lambda$-Übergang in flüssigem Schwefel

Problemstellung

Schwefel weist ein komplexes Phasendiagramm auf, das durch eine enorme molekulare Vielfalt geprägt ist – von zyklischen Ringen ($S_8, S_6$) bis hin zu linearen Polymerketten. Ein zentrales Phänomen ist der sogenannte $\lambda$-Übergang in der flüssigen Phase bei Umgebungsdruck: Bei etwa 432 K findet eine temperaturinduzierte Polymerisation statt, bei der sich $S_8$-Ringe öffnen und zu langen Ketten verbinden. Dieser Übergang ist durch drastische Änderungen der physikalischen Eigenschaften (Viskosität, Wärmekapazität, Farbe) gekennzeichnet.

Bisherige Herausforderungen waren:

1. Modellierung: Die Simulation dieses Prozesses auf atomarer Ebene ist extrem rechenintensiv, da die Zeitskalen der Polymerisation und die benötigten Systemgrößen die Kapazitäten klassischer Dichtefunktionaltheorie (DFT)-Molekulardynamik (MD) übersteigen.
2. Mechanismus: Es war unklar, ob die Polymerisation durch die Öffnung von $S_8$-Ringen (Tobolsky-Eisenberg-Modell) oder durch andere Spezies (wie nicht-oktazyklische Ringe) eingeleitet wird.
3. Hochdruckregime: Die Frage, wie der Übergang verläuft, wenn die Schmelzkurve und die Polymerisationslinie bei moderaten Drücken (ca. 0,1 GPa) konvergieren.

Methodik

Die Autoren nutzen eine hochmoderne computergestützte Methode, um die Lücke zwischen Genauigkeit und Effizienz zu schließen:

• Machine-Learned Interatomic Potential (MLIP): Es wurde ein lokales, äquivariantes ML-Potenzial (das Allegro-Modell) entwickelt. Dieses wurde auf einem Datensatz trainiert, der mittels Aktivem Lernen (Query-by-Committee-Strategie) aus DFT-Rechnungen (GGA-PBE-Niveau) gewonnen wurde.
• Molekulardynamik (MD): Mit dem ML-Potenzial wurden großskalige MD-Simulationen durchgeführt. Dabei wurden Temperaturrampen (T-Ramps) bei verschiedenen Drücken (von 1 atm bis 1,0 GPa) verwendet, um sowohl das Schmelzen als auch den $\lambda$-Übergang spontan zu beobachten.
• Strukturelle Analyse: Zur Charakterisierung der Phasen wurden topologische Analysen (Koordinationszahlen), die Bestimmung der Ringgrößen ($S_\pi$) und die SOAP-Ähnlichkeitsanalyse (Smooth Overlap of Atomic Positions) eingesetzt.

Wesentliche Ergebnisse

1. Mikroskopischer Mechanismus des $\lambda$-Übergangs: Die Simulationen bestätigen, dass nicht-$S_8$-Ringe (insbesondere $S_7$ und andere Ringe der Größe 6–20) als reaktive Zentren fungieren. Diese Spezies sind am Beginn der Polymerisation überrepräsentiert und dienen als "Keime" für die Kettenbildung. Dies stellt das klassische Modell der direkten $S_8$-Ringöffnung teilweise in Frage.
2. Reproduktion experimenteller Signaturen: Das Modell reproduziert die scharfe Spitze der Wärmekapazität ($C_s$) und die starke Abhängigkeit der Übergangstemperatur $T_\lambda$ von der Heizrate. Durch Extrapolation der Heizraten wurde ein Wert von $T_\lambda \approx 500\text{ K}$ ermittelt, was nahe am experimentellen Wert liegt.
3. Phasendiagramm und kritischer Punkt: Die Autoren rekonstruierten das Phasendiagramm und zeigten, dass die Polymerisationstemperatur mit steigendem Druck moderat sinkt, bis sie mit der Schmelzkurve verschmilzt.
4. Polymerisation im Festkörper (Hochdruckregime): Bei Drücken $\geq 0,5\text{ GPa}$ beobachteten die Autoren, dass die Polymerisation bereits im kristallinen Zustand (vor dem Schmelzen) beginnt. Nicht-$S_8$-Ringe und Polymerketten bilden sich innerhalb des Gitters und passen sich dessen Geometrie an (topochemische Polymerisation). Das Schmelzen erfolgt dann direkt aus diesen polymerisierten Anteilen.

Bedeutung der Arbeit

Diese Arbeit liefert das erste atomare Bild des gesamten Schwefel-Phasendiagramms, einschließlich des komplexen Übergangs zwischen molekularer und polymerer Phase. Durch den Einsatz von Machine-Learning-Potenzialen wurde gezeigt, dass es möglich ist, chemische Bindungsbruch- und -bildungsprozesse in komplexen Systemen mit der Genauigkeit der Quantenchemie, aber der Effizienz klassischer Mechanik zu simulieren. Die Erkenntnisse über die Rolle von Nicht-$S_8$-Ringen und die Polymerisation im Festkörper bieten neue theoretische Grundlagen für das Verständnis von Schwefel und ähnlichen Elementen unter extremen Bedingungen.