Improved Correlations for the Static Dielectric Constant of Hydrogen

Technische Samenvatting: Verbeterde Correlaties voor de Statische Diëlektrische Constante van Waterstof

Probleemstelling
Nauwkeurige voorspelling van de statische diëlektrische constante van waterstof is cruciaal voor fundamentele thermodynamische berekeningen en opgedifferentieerde sensortoepassingen op basis van diëlektrica, met name in systemen voor de liquefactie, opslag en verwerking van waterstof. Eerdere referentiecorrelaties, specifiek het model van Harvey en Lemmon uit 2005 [6], vertonen significante tekortkomingen bij cryogene condities. Deze beperkingen omvatten:

• Een gebrek aan specifieke correlaties voor normale waterstof (n-H₂) en ortho-waterstof (o-H₂), wat de toepassing beperkt tot para-waterstof (p-H₂) en evenwichtswaterstof (e-H₂).
• Een te versimpelde berekening bij lage dichtheid die slechts de grondtoestand van de vibratie beschouwt, waarbij de geëxciteerde toestanden die bij hogere temperaturen (bijv. 500 K) worden bezet, worden genegeerd.
• Een lineaire temperatuurafhankelijkheid voor de bijdrage bij lage dichtheid die het niet-lineaire gedrag dat bij cryogene temperaturen wordt waargenomen, niet kan vatten.
• Onvoldoende gebruik van beschikbare experimentele gegevens, met name voor vloeistoffen en samengeperste gassen onder hoge druk, wat leidt tot een slechte overeenstemming met recente metingen van de diëlektrische contrasten van spin-isomeren.

Methodologie
De auteurs hebben verbeterde referentiecorrelaties ontwikkeld voor de statische diëlektrische constante van p-H₂ en o-H₂ over een temperatuurbereik van 1 K tot 1000 K. De methodologie bestond uit drie primaire fasen:

1. Herberekening van de Limieten bij Lage Dichtheid ($A_\varepsilon$):
   De eerste diëlektrische virialcoëfficiënt, $A_\varepsilon(T)$, werd herberekend met behulp van moderne ab initio rotatie-vibratie energieniveaus (Komasa et al. [9]) en toestandsafhankelijke moleculaire polariseerbaarheden. In tegenstelling tot eerdere modellen omvatte deze aanpak:

   • Niet-adiabatische, relativistische en kwantumelektrodynamische (QED) effecten in energieniveaus.
   • De bezetting van de eerste geëxciteerde vibratietoestand ($\nu=1$) en hogere rotationele toestanden, die niet verwaarloosbaar zijn bij temperaturen tot 1000 K.
   • Aanpassingen aan polariseerbaarheidswaarden om rekening te houden met de breuk van de Born–Oppenheimer-benadering voor lichte kernen.
   • Distincte partitiefuncties voor p-H₂ (even rotationele toestanden) en o-H₂ (oneven rotationele toestanden) om rekening te houden met de kernspinstatistiek.

2. Ontwikkeling van een Nieuwe Functionele Vorm (SITH-functie):
   Om de temperatuurafhankelijkheid van $A_\varepsilon(T)$ nauwkeurig weer te geven over het volledige bereik, werd een gemodificeerde Padé-benadering, de "SITH"-functie genoemd, ontwikkeld. Deze functionele vorm vangt de sterke kromming en asymptotische gedragingen bij cryogene temperaturen die lineaire benaderingen missen.

3. Review van de Experimentele Database en Fitting bij Hoge Dichtheid:
   Er is een uitgebreide review uitgevoerd van historische en moderne experimentele gegevens voor p-H₂ en n-H₂, waarbij het volgende is gebeurd:

   • Conversie van historische temperatuurschalen (NBS-55, IPTS-68) naar de Internationale Temperatuurschaal van 1990 (ITS-90).
   • Herwaardering van dichtheidsberekeningen met behulp van de referentie-toestandsvergelijking door Leachman et al. [21].
   • Correctie voor kalibratieproblemen in historische capaciteitsmetingen.
   • Het fitten van de hogere-orde virialcoëfficiënten ($B_\varepsilon$ en $C_\varepsilon^+$) aan de experimentele gegevens. Voor o-H₂, waar directe metingen bij hoge dichtheid niet beschikbaar zijn, werden correlaties indirect afgeleid door n-H₂ gegevens te fitten met behulp van een lineaire polarisatie-mengregel met de p-H₂ correlatie.

Belangrijkste Bijdragen

• De SITH-functie: Een nieuwe, continue functionele representatie voor de eerste diëlektrische virialcoëfficiënt ($A_\varepsilon$) die de temperatuurafhankelijkheid van zowel p-H₂ als o-H₂ van 1 K tot 1000 K nauwkeurig modelleert, ter vervanging van de inadequate lineaire vormen van eerdere modellen.
• Spin-Isomeer Specifieke Correlaties: De ontwikkeling van onderscheidende referentiecorrelaties voor p-H₂ en o-H₂, wat de eerdere onmogelijkheid om n-H₂ en o-H₂ nauwkeurig te modelleren bij cryogene temperaturen aanpakt.
• Uitgebreide Experimentele Database: De integratie van een breder scala aan experimentele gegevens, inclus\nibil de recente hoogprecisie metingen van het diëlektrische contrast tussen spin-isomeren en historische vloeibare waterstofgegevens, met een strikte onzekerheidsanalyse en temperatuurschaalconversies.
• Verbeterde Representatie bij Hoge Dichtheid: Herziene coëfficiënten voor de tweede ($B_\varepsilon$) en effectieve derde ($C_\varepsilon^+$) diëlektrische virialcoëfficiënt die de cryogene vloeistof- en superkritische data beter reproduceren.

Resultaten

• Nauwkeurigheid bij Lage Dichtheid: De SITH-functie reproduceert ab initio berekeningen binnen hun geschatte onzekerheden. Het voorspelt een verschil in polariseerbaarheid bij lage dichtheid tussen n-H₂ en p-H₂ van 0,13% bij 20,4 K, wat consistent is met experimentele rapporten van ~0,15%.
• Kwaliteit van de High-Density Fit:
  • Voor p-H₂ verminderde de nieuwe correlatie de gemiddelde absolute afwijking (AAD) in de diëlektrische constante van 0,0098% (Harvey en Lemmon) naar 0,0060%. De wortel van de gemiddelde gekwadrateerde residuen van genormaliseerde residuen ($RMS_z$) verbeterde van 0,52 naar 0,45.
  • Voor o-H₂ (afgeleid via n-H₂), verbeterde de AAD van 0,0389% naar 0,0215%, waarbij $RMS_z$ daalde van 2,11 naar 1,03.
• Onafhankelijke Validatie: De nieuwe correlaties vertonen een aanzienlijk verbeterde overeenstemming met onafhankelijke literatuurmetingen die niet voor de fitting zijn gebruikt, waaronder:
  • Recente metingen van het diëlektrische contrast tussen n-H₂ en p-H₂ door Sellner et al. [5].
  • Metingen van de diëlektrische constante van niet-evenwichtso-p mengsels door Milenko en Sibileva [38], waarbij de totale residu-verstrooiing ongeveer gehalveerd was vergeleken met het Harvey en Lemmon-model.

Betekenis en Beperkingen
De auteurs stellen dat deze correlaties de meest uitgebreide representatie bieden die momenteel beschikbaar is voor de statische diëlektrische constante van waterstof spin-isomeren, wat de voorspellende kracht voor cryogene toepassingen zoals liquefactie, opslag en sensing aanzienlijk verbetert. Het werk adresseert specifiek de tekortkomingen van het model uit 2005 in de vloeistof- en superkritische regio's.

Het artikel behoudt echter een bescheiden houding met betrekking tot de beperkingen van het huidige werk:

• Zwakke Beperkingen op Hogere-Orde Termen: De temperatuurafhankelijkheden van de termen met hogere dichtheid ($B_\varepsilon$ en $C_\varepsilon^+$) blijven slechts zwak begrensd door de beschikbare experimentele database, met name voor intermediaire temperaturen en dichtheden.
• Empirisch Karakter: De termen voor hogere dichtheid vertrouwen op empirische vormen die oorspronkelijk voor gasvormige systemen zijn ontwikkeld, en het blijft onduidelijk of deze volledig de gedragingen in de vloeistofregio kunnen vatten.
• Indirecte o-H₂ Bepaling: De o-H₂ correlatie is indirect afgeleid van n-H₂ gegevens vanwege een gebrek aan directe zuivere o-H₂ metingen, wat een afhankelijkheid introduceert van mengingsaannames.
• Behoefte aan Toekomstige Data: De auteurs merken op dat toekomstige verbeteringen nieuwe experimentele metingen vereisen bij matige dichtheden en intermediaire temperaturen, evenals theoretische ab initio berekeningen van de tweede diëlektrische virialcoëfficiënt ($B_\varepsilon$) om empirische behandelingen te vervangen.