Local-Basis-Function Equation of State for Ice VII–X to 450 GPa at 300 K

Abstract

Helmholtz energy of ice VII−X is determined in a pressure regime extending to 450 GPa at 300 K using local-basis-functions in the form of b-splines. The new representation for the equation of state is embedded in a physics-based inverse theory framework of parameter estimation. Selected pressures as a function of volume from 14 prior experimental studies and two theoretical studies constrain the behavior of Helmholtz energy. Separately measured bulk moduli, not used to construct the representation, are accurately replicated below about 20 GPa and above 60 GPa. In the intermediate range of pressure, the experimentally determined moduli are larger and have greater scatter than values predicted using the Helmholtz representation. Although systematic error in the determination of elastic moduli is possible and likely, the alternative hypothesis is a slow relaxation time associated with changes in proton mobility or the ice VII to X transition. A correlation is observed between anomalies in the pressure derivative of the predicted bulk modulus and previously suggested higher-order phase transitions. Improved determinations of elastic properties at high pressure would allow refinement of the current equation of state. More generally, the current method of data assimilation is broadly applicable to other materials in high-pressure studies and for investigations of planetary interiors.

Keywords

• equation of state
• helmholtz energy
• phase transition
• ice vii
• ice x
• nacl
• exoplanets
• icy/ocean worlds
• local-basis-function
• b-spline
• tikhonov inverse