Third and fourth density and acoustic virial coefficients of neon from first-principles calculations

In dieser Studie wurden die dritten und vierten Dichte- und akustischen Virialkoeffizienten von Neon im Temperaturbereich von 10 bis 5000 K mithilfe von Pfadintegral-Monte-Carlo-Simulationen und hochpräzisen ab-initio-Potentialen erstmals aus ersten Prinzipien bestimmt, wobei die resultierenden Unsicherheiten deutlich unter denen fast aller experimenteller Daten liegen.

Das Wesentliche

Die Geschichte von den drei Neon-Partnern

Stellen Sie sich Neon als eine riesige Menge winziger, unsichtbarer Bälle vor, die in einem Raum schweben. Diese Bälle sind die Atome. Wenn man Neon als Gas verwendet (zum Beispiel in Thermometern, die die Temperatur messen), ist es wichtig zu wissen, wie sich diese Bälle gegenseitig beeinflussen.

In der Wissenschaft gibt es eine einfache Regel: Je näher die Bälle kommen, desto mehr stoßen sie sich ab oder ziehen sich an.

1. Das alte Spiel: Nur zu zweit (Paar-Wechselwirkung)

Bisher haben Wissenschaftler meistens nur betrachtet, wie sich zwei Neon-Bälle verhalten, wenn sie sich nähern. Das ist wie ein Tanz zu zweit. Wenn zwei Partner tanzen, wissen wir genau, wie sie sich bewegen. Das war schon sehr gut erforscht (durch eine vorherige Studie von Hellmann et al.).

Aber in der echten Welt gibt es selten nur zwei Bälle. Oft sind es drei, vier oder mehr, die gleichzeitig im Raum sind. Und hier wird es kompliziert.

2. Das neue Spiel: Die Dreier- und Vierer-Gruppen

Die Autoren dieser neuen Studie haben sich gefragt: „Was passiert, wenn drei Neon-Bälle gleichzeitig in einem Dreieck stehen? Oder sogar vier in einer Pyramide?"

Stellen Sie sich vor, drei Freunde (die Neon-Atome) stehen in einem Kreis.

• Wenn nur zwei von ihnen sich unterhalten, ist das einfach.
• Aber wenn der dritte dazukommt, verändert sich das Gespräch. Der dritte Freund beeinflusst, wie die ersten beiden miteinander reden. Das nennt man nicht-additive Wechselwirkung. Es ist nicht einfach die Summe von „Freund A + Freund B" und „Freund B + Freund C". Der ganze Kreis hat eine eigene Dynamik.

Die Wissenschaftler haben mit super-leistungsfähigen Computern (die wie riesige Rechenmaschinen funktionieren) genau berechnet, wie diese Dreier-Gruppen und Vierer-Gruppen sich verhalten. Sie haben dafür die modernsten mathematischen Werkzeuge benutzt, die es gibt, um die kleinsten Details der atomaren Kräfte zu verstehen.

3. Warum ist das so wichtig? (Die Thermometer-Frage)

Warum interessiert das jemanden, der kein Physiker ist?

Neon wird in der Metrologie (der Wissenschaft vom Messen) verwendet, um Temperatur und Druck extrem genau zu messen. Man nutzt dafür sogenannte „akustische Gasthermometer". Das sind Geräte, die messen, wie schnell sich Schallwellen durch Neon-Gas bewegen.

• Das Problem: Wenn man die Temperatur messen will, muss man genau wissen, wie sich das Gas verhält. Wenn man die „Dreier-Gruppen" ignoriert, ist die Berechnung wie eine Landkarte, auf der die Berge fehlen. Die Vorhersage ist dann nicht präzise genug für die modernsten Messgeräte der Welt.
• Die Lösung: Mit diesen neuen, extrem genauen Berechnungen können Wissenschaftler nun die Eigenschaften von Neon so genau vorhersagen, dass sie sogar genauer sind als die besten Messgeräte, die man heute bauen kann.

4. Der Vergleich mit anderen Gasen

• Helium: Das ist der „Superstar". Man kann sein Verhalten so genau berechnen, dass es fast perfekt ist. Aber Helium ist sehr leicht und empfindlich.
• Argon: Das ist schwerer und robuster, aber man kann sein Verhalten nicht so genau berechnen wie bei Helium, weil es zu viele Elektronen (die kleinen Teilchen im Atom) hat, die das Rechnen schwer machen.
• Neon: Es liegt genau in der Mitte. Es ist schwerer als Helium (also weniger empfindlich für Verunreinigungen), aber man kann es fast so genau berechnen wie Helium. Das macht es zum perfekten Kandidaten für die Zukunft der Messung.

5. Was haben die Forscher konkret getan?

1. Sie haben die „Dreier-Regeln" erfunden: Sie haben berechnet, wie drei Neon-Atome zusammenwirken. Das war wie das Lösen eines riesigen 3D-Puzzles mit tausenden von Teilen.
2. Sie haben die „Vierer-Regeln" geprüft: Sie haben auch geschaut, was bei vier Atomen passiert. Das Ergebnis? Das ist so winzig klein, dass es für die meisten Zwecke fast keine Rolle spielt, aber sie haben es trotzdem berechnet, um sicherzugehen.
3. Sie haben Unsicherheiten berechnet: Nicht nur den Wert selbst, sondern auch: „Wie sicher sind wir uns?" Die Antwort ist: Wir sind uns so sicher, dass unsere Berechnung besser ist als fast alle bisherigen Experimente.

Das Fazit in einem Satz

Diese Studie hat die „Spielregeln" für Neon-Atome so genau berechnet, dass wir nun bessere Thermometer bauen können, die die Temperatur der Welt mit einer Präzision messen, die früher unmöglich schien – alles dank einer genauen Analyse davon, wie drei oder vier Neon-Atome zusammen tanzen.

Kurz gesagt: Die Forscher haben das Verhalten von Neon-Atomen in Gruppen so perfekt verstanden, dass wir jetzt noch genauere Messgeräte für Temperatur und Druck bauen können.

Technische Zusammenfassung

Titel: Dritte und vierte Dichte- und akustische Virialkoeffizienten von Neon aus Erstprinzipien-Rechnungen

Autoren: Robert Hellmann und Giovanni Garberoglio
Institutionen: Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg (Deutschland) und FBK-ECT* (Italien)

---

1. Problemstellung und Motivation

Edelgase, insbesondere Neon, sind entscheidende Modellsubstanzen für die Untersuchung thermophysikalischer Eigenschaften und werden zunehmend in der Gas-Metrologie (Temperatur- und Druckmessung) als Arbeitsgase eingesetzt. Während die Eigenschaften von Helium bereits mit extrem hoher Genauigkeit aus quantenchemischen ab initio-Rechnungen vorhergesagt werden können, ist dies für Neon aufgrund der höheren Elektronenzahl und der damit verbundenen Komplexität der Elektronenkorrelationen schwieriger.

Bisherige Arbeiten (z. B. Hellmann et al., 2021) lieferten hochpräzise Paarpotentiale und den zweiten Virialkoeffizienten ($B$). Für Anwendungen in der Metrologie bei höheren Dichten und Temperaturen ist jedoch die Kenntnis höherer Virialkoeffizienten (dritte und vierte Ordnung) unerlässlich. Diese hängen nicht nur von Paarkräften, sondern auch von nicht-additiven Mehrkörperwechselwirkungen (Drei- und Vier-Körper-Effekte) ab. Bisherige experimentelle Daten für diese Koeffizienten weisen oft große Unsicherheiten auf oder sind inkonsistent.

Ziel der Arbeit: Die Bestimmung der dritten ($C$) und vierten ($D$) Dichte-Virialkoeffizienten sowie der dritten ($\gamma_a$) und vierten ($\delta_a$) akustischen Virialkoeffizienten für Neon im Temperaturbereich von 10 bis 5000 K mit bisher unerreichter Genauigkeit aus ersten Prinzipien.

---

2. Methodik

Die Studie kombiniert hochpräzise quantenchemische Berechnungen mit fortgeschrittenen statistisch-mechanischen Simulationsmethoden.

A. Potentiale und Wechselwirkungen

1. Paarpotential ($V_{12}$): Es wurde das bereits von Hellmann et al. (2021) entwickelte ab initio-Paarpotential wiederverwendet. Dieses basiert auf Supermolekül-Rechnungen mit Basis-Sätzen bis zur Oktuple-Zeta-Qualität und Theorien bis hin zu CCSDTQ(P) (Coupled Cluster mit Single-, Double-, Triple-, Quadruple- und perturbativen Pentuple-Anregungen). Es berücksichtigt relativistische Effekte, Verzögerungseffekte (Retardation) und post-Born-Oppenheimer-Korrekturen.
2. Nicht-additives Drei-Körper-Potential ($\Delta V_{123}$): Ein neues Potential wurde entwickelt, basierend auf 2550 Konfigurationen von Neon-Tripletts. Die Berechnungen umfassten CCSD(T), CCSDT und CCSDT(Q) mit Basis-Sätzen bis zur Sextuple-Zeta-Qualität. Die Ergebnisse wurden auf den vollständigen Basissatz-Limit (CBS) extrapoliert. Eine analytische Anpassungsfunktion (inspiriert von Axilrod-Teller-Muto für den Fernbereich) wurde erstellt.
3. Nicht-additives Vier-Körper-Potential ($\Delta V_{1234}$): Da dieser Beitrag als sehr klein erwartet wurde, wurde er stark vereinfacht behandelt. Es wurden nur reguläre Tetraeder-Konfigurationen mit CCSD(T) berechnet. Eine modifizierte Bade-Potential-Funktion wurde angepasst, um den Kurz- und Fernbereich korrekt abzubilden.

B. Berechnung der Virialkoeffizienten

• Methode: Path-Integral-Monte-Carlo (PIMC). Dies ermöglicht eine vollquantenmechanische Behandlung der Teilchen, was bei tiefen Temperaturen entscheidend ist.
• Dichte-Virialkoeffizienten ($C, D$): Diese wurden direkt durch PIMC-Simulationen berechnet. Dabei wurden die ring-polymerischen Darstellungen der Quantenteilchen verwendet, um die mehrdimensionalen Integrale über die Konfigurationsräume zu lösen.
• Akustische Virialkoeffizienten ($\gamma_a, \delta_a$):
  • Der dritte akustische Koeffizient ($\gamma_a$) wurde direkt via PIMC berechnet.
  • Der vierte akustische Koeffizient ($\delta_a$) wurde indirekt über thermodynamische Beziehungen aus den Dichte-Virialkoeffizienten und deren Temperaturableitungen abgeleitet, da eine direkte PIMC-Berechnung zu rechenintensiv wäre.
• Isotopenmischung: Berechnungen wurden für alle natürlichen Isotopenkombinationen ($^{20}$Ne, $^{21}$Ne, $^{22}$Ne) sowie für die gemittelte molare Masse durchgeführt.

C. Unsicherheitsanalyse

Ein wesentlicher Aspekt ist die rigorose Fortpflanzung der Unsicherheiten der Potentiale in die Virialkoeffizienten. Anstatt nur das Potential zu verschieben, wurde eine funktionale Variation verwendet, um die Unsicherheit als Integral über den Betrag der Ableitung des Koeffizienten nach dem Potential zu bestimmen. Dies liefert eine realistischere Obergrenze für die Unsicherheit, insbesondere bei tiefen Temperaturen.

---

3. Wichtige Beiträge und Ergebnisse

• Präzision: Die berechneten Unsicherheiten für die Virialkoeffizienten sind deutlich kleiner als die der meisten verfügbaren experimentellen Daten.
• Einfluss der Mehrkörpereffekte:
  • Der Beitrag der nicht-additiven Drei-Körper-Wechselwirkungen ist signifikant und übersteigt den Quanteneffekt bei weitem.
  • Der Vier-Körper-Beitrag ist extrem klein (ca. zwei Größenordnungen kleiner als der Drei-Körper-Beitrag) und vernachlässigbar für die Gesamtunsicherheit.
• Vergleich mit Experimenten:
  • Dritter Virialkoeffizient ($C$): Die neuen PIMC-Werte stimmen hervorragend mit re-analysierten Daten von Gaiser und Fellmuth (dielektrische Gasthermometrie) sowie mit Daten von von Preetzmann und Span überein. Bisherige theoretische Studien (z. B. Bich et al., Wiebke et al.) wichen aufgrund einer weniger genauen Behandlung der Drei-Körper-Effekte stärker von den experimentellen Daten ab.
  • Vierter Virialkoeffizient ($D$): Es gibt kaum verlässliche experimentelle Daten. Die neuen Werte stimmen mit früheren theoretischen Berechnungen von Wiebke et al. überein, zeigen aber eine viel kleinere Unsicherheit.
  • Akustische Koeffizienten ($\gamma_a, \delta_a$): Der Vergleich mit hochpräzisen Schallgeschwindigkeitsmessungen von Dietl et al. zeigt eine gute Übereinstimmung innerhalb der erweiterten experimentellen Unsicherheit ($k=2$).
• Analytische Anpassungen: Für alle Koeffizienten wurden analytische Funktionen bereitgestellt, die eine einfache Interpolation im gesamten Temperaturbereich (10–5000 K) ermöglichen.

---

4. Bedeutung und Ausblick

Diese Arbeit stellt einen Meilenstein in der theoretischen Beschreibung von Neon dar.

• Metrologie: Die Ergebnisse ermöglichen die Nutzung von Neon als alternatives Arbeitsgas in der akustischen Gasthermometrie und bei Druckstandards, insbesondere dort, wo Helium aufgrund von Verunreinigungen oder Isotopenvariationen Nachteile hat.
• Validierung: Die hohe Genauigkeit der Berechnungen dient als Referenz, um experimentelle Messungen zu validieren und bestehende Datenbestände zu überprüfen.
• Methodischer Fortschritt: Die rigorose Unsicherheitsfortpflanzung und die Kombination aus hochrangigen CCSDTQ(P)-Berechnungen mit PIMC setzen einen neuen Standard für die Berechnung thermophysikalischer Eigenschaften von Fluiden aus ersten Prinzipien.

Zusammenfassend liefert das Paper die derzeit genauesten Vorhersagen für die thermodynamischen und akustischen Eigenschaften von Neon, die für die nächste Generation von Primärmessstandards in der Temperatur- und Druckmetrologie essenziell sind.