PALEOS MgSiO₃ EoS Lookup Table for PROTEUS

Some points are missing where the EoS solver failed to converge; when reconstructing the full grid for interpolation, these entries should be filled with NaN. The nonconverging regions should not affect the physically relevant domain (beyond liquidus for the solid table, vapor/supercritical phase for the liquid table). Columns Each row contains 10 columns, all in SI units: P — Pressure [Pa] T — Temperature [K] rho — Density [kg/m³] u — Specific internal energy [J/kg] s — Specific entropy [J/(kg·K)] cp — Isobaric heat capacity [J/(kg·K)] cv — Isochoric heat capacity [J/(kg·K)] alpha — Thermal expansion coefficient [K⁻¹] nabla_ad — Adiabatic gradient (d ln T / d ln P)_S [dimensionless] phase — Phase identifier Phase selection For the solid table, phase selection at each (P, T) point follows the full PALEOS MgSiO₃ phase diagram.

At temperatures above the PALEOS melting curve (which can occur since the solid table extends to the PROTEUS liquidus), the code falls back to the solid phase stable just below the melting curve and evaluates that EoS at the requested temperature. The liquid table uses Wolf+Bower18 everywhere. Generating the tables The notebook to regenerate these tables is available in the PALEOS repository at utils/mgsio3_tables_proteus.ipynb.

It contains a single clearly marked User Input cell where all parameters (file paths, P/T ranges, resolution, output filenames) can be adjusted.

• *
• Version log (1.1.0): covered all PROTEUS solidi/liquidi, implemented consistent U₀/S₀ per-phase shifts to ensure dS/dT > 0 everywhere and Gibbs-energy continuity at triple points, and generated high-resolution (600 ppd) tables. Mirror of Zenodo deposit 10.5281/zenodo.19680050. Zenodo is the primary source.