🔬 physics

Improved Correlations for the Static Dielectric Constant of Hydrogen

This paper presents improved correlations for the static dielectric constant of para- and ortho-hydrogen over a temperature range of 1 to 1000 K by recalculating low-density contributions with modern ab initio methods and developing new functional representations based on an extensive review of experimental data, thereby significantly enhancing accuracy for cryogenic and liquid states compared to previous models.

Original authors: Guinevere M. Sellner, Liam D. Tenardi, Paul L. Stanwix, Eric F. May, Allan H. Harvey

Published 2026-08-11

📖 3 min read☕ Coffee break read

CC BY 4.0

Original authors: Guinevere M. Sellner, Liam D. Tenardi, Paul L. Stanwix, Eric F. May, Allan H. Harvey

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to listen to a whisper in a crowded room. To hear that whisper clearly, you need a microphone that doesn't distort the sound, and you need to know exactly how the air in the room behaves when it gets cold or hot. In the world of science, this "microphone" is often an electric field, and the "air" is a gas or liquid. Scientists use a property called the static dielectric constant to measure how a material reacts to an electric field. Think of it as a "squishiness" factor: when you push on a material with electricity, how much does it squish or stretch? This number is crucial for building sensors that can detect tiny leaks in hydrogen tanks or measure how much fuel is left in a rocket.

However, hydrogen is a tricky customer. It comes in two "personalities" or spin-isomers: para-hydrogen and ortho-hydrogen. They are chemically identical but behave differently when they spin, much like how a left-handed glove fits a left hand but feels weird on a right hand. For decades, scientists have used a map to predict how these hydrogen personalities react to electricity, but that map was drawn with a straight ruler. It worked okay for warm, gassy hydrogen, but when things got super cold—like the freezing temperatures needed to turn hydrogen into a liquid for rockets—the straight-line map started to fail, leading to confusing errors in our sensors.

This paper is about redrawing that map with a flexible, curvy ruler that actually fits the terrain. The authors, a team of researchers from Australia and the United States, realized that the old rules didn't work for the extreme cold of liquid hydrogen, especially for the "ortho" version which is common in normal hydrogen but behaves strangely at low temperatures. They went back to the drawing board, using modern super-computer calculations to figure out exactly how the hydrogen molecules vibrate and spin at different temperatures. They found that the "squishiness" of hydrogen isn't a straight line; it curves and twists, especially when it's freezing.

By combining these new computer calculations with a massive review of old and new experiments, they created a brand-new set of formulas, which they jokingly named the SITH function (after their initials). This new model is much better at predicting how hydrogen behaves from a chilly 1 Kelvin up to a scorching 1000 Kelvin. The paper shows that their new map fits the experimental data much better than the old one, particularly for liquid hydrogen and for mixtures of the two hydrogen types.

However, the authors are careful not to claim they have solved every mystery. They point out that while their new model is a huge improvement, there are still some gaps in the data, especially for the "ortho" hydrogen at high densities. Because they couldn't measure pure ortho-hydrogen directly in many experiments, they had to estimate its behavior by mixing it with normal hydrogen in their calculations. The paper suggests that while their new correlations are the best we have right now and are ready to help build better sensors, future scientists will need to take more measurements in the middle-temperature ranges to make the map even more perfect. For now, though, this new "curvy ruler" gives us a much clearer way to listen to the whispers of hydrogen in our future energy systems.

Technical Summary: Improved Correlations for the Static Dielectric Constant of Hydrogen

Problem Statement
Accurate prediction of the static dielectric constant of hydrogen is critical for fundamental thermophysical calculations and emerging dielectric-based sensing applications, particularly in hydrogen liquefaction, storage, and processing systems. Previous reference correlations, specifically the 2005 model by Harvey and Lemmon [6], exhibit significant deficiencies at cryogenic conditions. These limitations include:

A lack of specific correlations for normal hydrogen (n-H₂) and ortho-hydrogen (o-H₂), restricting application to para-hydrogen (p-H₂) and equilibrium hydrogen (e-H₂).
An oversimplified low-density calculation that considers only the ground vibrational state, neglecting excited states that become populated at higher temperatures (e.g., 500 K).
A linear temperature dependence for the low-density contribution that fails to capture the non-linear behavior observed at cryogenic temperatures.
Insufficient use of available experimental data, particularly for liquid states and high-pressure compressed gases, leading to poor agreement with recent measurements of spin-isomer dielectric contrasts.

Methodology
The authors developed improved reference correlations for the static dielectric constant of p-H₂ and o-H₂ over a temperature range of 1 K to 1000 K. The methodology involved three primary stages:

Recalculation of Low-Density Limits ( $A_\varepsilon$ ):
The first dielectric virial coefficient, $A_\varepsilon(T)$ , was recalculated using modern ab initio rovibrational energy levels (Komasa et al. [9]) and state-dependent molecular polarizabilities. Unlike previous models, this approach included:
- Nonadiabatic, relativistic, and quantum electrodynamics (QED) effects in energy levels.
- Population of the first excited vibrational state ( $\nu=1$ ) and higher rotational states, which are non-negligible at temperatures up to 1000 K.
- Adjustments to polarizability values to account for the breakdown of the Born–Oppenheimer approximation for light nuclei.
- Distinct partition functions for p-H₂ (even rotational states) and o-H₂ (odd rotational states) to account for nuclear spin statistics.
Development of a New Functional Form (SITH Function):
To represent the temperature dependence of $A_\varepsilon(T)$ accurately across the full range, a modified Padé approximant, termed the "SITH" function, was developed. This functional form captures the strong curvature and asymptotic behavior at cryogenic temperatures that linear approximations miss.
Experimental Database Review and High-Density Fitting:
An extensive review of historical and modern experimental data for p-H₂ and n-H₂ was conducted. This involved:
- Converting historical temperature scales (NBS-55, IPTS-68) to the International Temperature Scale of 1990 (ITS-90).
- Re-evaluating density calculations using the reference equation of state by Leachman et al. [21].
- Correcting for calibration issues in historical capacitance measurements.
- Fitting the higher-order virial coefficients ( $B_\varepsilon$ and $C_\varepsilon^+$ ) to the experimental data. For o-H₂, where direct high-density measurements are unavailable, correlations were derived indirectly by fitting n-H₂ data using a linear polarization mixing rule with the p-H₂ correlation.

Key Contributions

The SITH Function: A new, continuous functional representation for the first dielectric virial coefficient ( $A_\varepsilon$ ) that accurately models the temperature dependence of both p-H₂ and o-H₂ from 1 K to 1000 K, replacing the inadequate linear forms of previous models.
Spin-Isomer Specific Correlations: The development of distinct reference correlations for p-H₂ and o-H₂, addressing the previous inability to model n-H₂ and o-H₂ accurately at cryogenic temperatures.
Expanded Experimental Database: The integration of a broader range of experimental data, including recent high-precision measurements of dielectric contrast between spin isomers and historical liquid hydrogen data, with rigorous uncertainty analysis and temperature scale conversions.
Improved High-Density Representation: Revised coefficients for the second ( $B_\varepsilon$ ) and effective third ( $C_\varepsilon^+$ ) dielectric virial coefficients that better reproduce cryogenic liquid and supercritical data.

Results

Low-Density Accuracy: The SITH function reproduces ab initio calculations within their estimated uncertainties. It predicts a low-density polarizability difference between n-H₂ and p-H₂ of 0.13% at 20.4 K, consistent with experimental reports of ~0.15%.
High-Density Fit Quality:
- For p-H₂, the new correlation reduced the average absolute deviation (AAD) in dielectric constant from 0.0098% (Harvey and Lemmon) to 0.0060%. The root-mean-square of normalized residuals ( $RMS_z$ ) improved from 0.52 to 0.45.
- For o-H₂ (derived via n-H₂), the AAD improved from 0.0389% to 0.0215%, with $RMS_z$ decreasing from 2.11 to 1.03.
Independent Validation: The new correlations show substantially improved agreement with independent literature measurements not used in the fitting process, including:
- Recent measurements of the dielectric contrast between n-H₂ and p-H₂ by Sellner et al. [5].
- Dielectric constant measurements of non-equilibrium o-p mixtures by Milenko and Sibileva [38], where the total residual scatter was approximately halved compared to the Harvey and Lemmon model.

Significance and Limitations
The authors state that these correlations provide the most comprehensive representation currently available for the static dielectric constant of hydrogen spin isomers, significantly improving predictive capability for cryogenic applications such as liquefaction, storage, and sensing. The work specifically addresses the deficiencies of the 2005 model in the liquid and supercritical regions.

However, the paper maintains a modest stance regarding the limitations of the current work:

Weak Constraints on Higher-Order Terms: The temperature dependences of the higher-density terms ( $B_\varepsilon$ and $C_\varepsilon^+$ ) remain only weakly constrained by the available experimental database, particularly for intermediate temperatures and densities.
Empirical Nature: The higher-density terms rely on empirical forms originally developed for gaseous systems, and it remains unclear if these fully capture liquid-region behavior.
Indirect o-H₂ Determination: The o-H₂ correlation is derived indirectly from n-H₂ data due to a lack of direct pure o-H₂ measurements, introducing reliance on mixing assumptions.
Need for Future Data: The authors note that future improvements require new experimental measurements at moderate densities and intermediate temperatures, as well as theoretical ab initio calculations of the second dielectric virial coefficient ( $B_\varepsilon$ ) to replace empirical treatments.

More like this