Reference Data
Download instructions
All tabulated data files are downloaded automatically by running the setup script from the repository root:
bash src/get_zalmoxis.sh
This invokes src/setup_zalmoxis.py, which downloads and extracts the required data into the data/ directory.
See Installation Step 4 of the installation guide for details.
Data inventory
The table below lists every data file used by Zalmoxis, organised by subdirectory.
Equation of state tables
| File | Location | Format | Source | Description |
|---|---|---|---|---|
eos_seager07_iron.txt |
data/EOS_Seager2007/ |
CSV (rho in g/cm^3, P in GPa) | Seager et al. (2007) | Fe epsilon pressure-density EOS |
eos_seager07_silicate.txt |
data/EOS_Seager2007/ |
CSV (rho in g/cm^3, P in GPa) | Seager et al. (2007) | MgSiO3 perovskite pressure-density EOS |
eos_seager07_water.txt |
data/EOS_Seager2007/ |
CSV (rho in g/cm^3, P in GPa) | Seager et al. (2007) | Water ice (phases VII, VIII, X) pressure-density EOS |
density_melt.dat |
data/EOS_WolfBower2018_1TPa/ |
TSV (P in Pa, T in K, rho in kg/m^3) | Wolf & Bower (2018) | MgSiO3 melt density EOS (P: 0--1 TPa, T: 0--16500 K). Out-of-bounds pressures are clamped to the table edge. |
density_solid.dat |
data/EOS_WolfBower2018_1TPa/ |
TSV (P in Pa, T in K, rho in kg/m^3) | Wolf & Bower (2018) | MgSiO3 solid (bridgmanite) density EOS, derived from Mosenfelder et al. (2009) (P: 0--1 TPa, T: 0--16500 K). Out-of-bounds pressures are clamped to the table edge. |
adiabat_temp_grad_melt.dat |
data/EOS_WolfBower2018_1TPa/ |
TSV (P in Pa, T in K, dT/dP in K/Pa) | Wolf & Bower (2018) | MgSiO3 melt adiabatic temperature gradient \((dT/dP)_S\) (P: 0--1 TPa, T: 0--16500 K). Used for "adiabatic" temperature mode. |
density_melt.dat |
data/EOS_RTPress_melt_100TPa/ |
TSV (P in Pa, T in K, rho in kg/m^3) | Extended RTpress melt table | MgSiO3 melt density EOS extended to 100 TPa (P: 1e3--1e14 Pa, T: 400--50000 K). Used by RTPress100TPa:MgSiO3. Solid phase uses the WolfBower2018 table above. |
adiabat_temp_grad_melt.dat |
data/EOS_RTPress_melt_100TPa/ |
TSV (P in Pa, T in K, dT/dP in K/Pa) | Extended RTpress melt table | MgSiO3 melt adiabatic temperature gradient \((dT/dP)_S\) extended to 100 TPa (P: 1e3--1e14 Pa, T: 400--50000 K). Used for "adiabatic" temperature mode with RTPress100TPa:MgSiO3. |
paleos_mgsio3_tables_pt_proteus_solid.dat |
data/EOS_PALEOS_MgSiO3/ |
Space-separated, 10 columns (P, T, rho, u, s, cp, cv, alpha, nabla_ad, phase_id) | PALEOS (Zenodo 18924171) | MgSiO3 solid-phase EOS with adiabatic gradient (P: 1 bar--100 TPa, T: 300 K--max(T_liquidus), 150 ppd log-uniform grid, 75% filled). |
paleos_mgsio3_tables_pt_proteus_liquid.dat |
data/EOS_PALEOS_MgSiO3/ |
Space-separated, 10 columns (same as solid) | PALEOS (Zenodo 18924171) | MgSiO3 liquid-phase EOS with adiabatic gradient (P: 1 bar--100 TPa, T: min(T_solidus)--100000 K, 150 ppd log-uniform grid, 53% filled). |
paleos_iron_eos_table_pt.dat |
data/EOS_PALEOS_iron/ |
Space-separated, 10 columns (P, T, rho, u, s, cp, cv, alpha, nabla_ad, phase) | PALEOS (Zenodo 19000316) | Unified Fe EOS, 5 stable phases (alpha-bcc, delta-bcc, gamma-fcc, epsilon-hcp, liquid). P: 1 bar--100 TPa, T: 300--100000 K. 75.5 MB. |
paleos_mgsio3_eos_table_pt.dat |
data/EOS_PALEOS_MgSiO3_unified/ |
Space-separated, 10 columns (same as iron) | PALEOS (Zenodo 19000316) | Unified MgSiO3 EOS, 6 stable phases (3 pyroxene, bridgmanite, postperovskite, liquid). P: 1 bar--100 TPa, T: 300--100000 K. 48.8 MB. |
paleos_water_eos_table_pt.dat |
data/EOS_PALEOS_H2O/ |
Space-separated, 10 columns (same as iron) | PALEOS (Zenodo 19000316) | Unified H2O EOS, 7 EOS (ice Ih--X, liquid, vapor, superionic). P: 1 bar--100 TPa, T: 100--100000 K. 147.3 MB. |
chabrier2021_H.dat |
data/EOS_Chabrier2021_HHe/ |
Space-separated, 10 columns (P, T, rho, u, s, cp, cv, alpha, nabla_ad, phase_id) | Chabrier et al. (2019), Chabrier & Debras (2021) | Pure H\(_2\) EOS (molecular, atomic, ionized). Grid: 121 x 441 (logT, logP), T: 100--\(10^8\) K, P: 1 Pa--\(10^{22}\) Pa. Registered as Chabrier:H. |
chabrier2021_HE.dat |
data/EOS_Chabrier2021_HHe/ |
Same as above | Same | Pure He EOS. Not registered (available for future use). |
chabrier2021_HHe_Y0275.dat |
data/EOS_Chabrier2021_HHe/ |
Same as above | Same | H/He mixture at \(Y = 0.275\). Not registered (available for future use). |
chabrier2021_HHe_Y0292.dat |
data/EOS_Chabrier2021_HHe/ |
Same as above | Same | H/He mixture at \(Y = 0.292\). Not registered (available for future use). |
chabrier2021_HHe_Y0297.dat |
data/EOS_Chabrier2021_HHe/ |
Same as above | Same | H/He mixture at \(Y = 0.297\). Not registered (available for future use). |
Melting curves
Tabulated (downloaded data files):
| File | Location | Format | Source | Description |
|---|---|---|---|---|
solidus.dat |
data/melting_curves_Monteux-600/ |
Space-separated (P in Pa, T in K) | Monteux et al. (2016) | MgSiO3 solidus curve (offset: solidus = liquidus - 600 K). Legacy; selected via rock_solidus = "Monteux600-solidus-tabulated". |
liquidus.dat |
data/melting_curves_Monteux-600/ |
Space-separated (P in Pa, T in K) | Monteux et al. (2016) | MgSiO3 liquidus curve. Legacy; selected via rock_liquidus = "Monteux600-liquidus-tabulated". |
Analytic (hardcoded, no data files needed):
| Identifier | Source | Description |
|---|---|---|
Monteux16-solidus |
Monteux et al. (2016) Eqs. 10/12 | Piecewise Simon-Glatzel solidus, valid to ~500 GPa. |
Monteux16-liquidus-A-chondritic |
Monteux et al. (2016) Eqs. 11/13 | A-chondritic liquidus, valid to ~500 GPa. |
Monteux16-liquidus-F-peridotitic |
Monteux et al. (2016) Eqs. 11/13 | F-peridotitic liquidus, valid to ~660 GPa. |
Stixrude14-solidus |
Stixrude (2014) Eqs. 1.9+1.10 | Simon-like power law with cryoscopic depression (\(x_0 = 0.79\)). Valid at all \(P\). |
Stixrude14-liquidus |
Stixrude (2014) Eq. 1.9 | Pure MgSiO\(_3\) Simon-like power law. Valid at all \(P\). |
PALEOS-liquidus |
Belonoshko et al. (2005) / Fei et al. (2021) | Piecewise Simon-Glatzel liquidus matching the unified PALEOS MgSiO\(_3\) table. Valid at all \(P\). |
See configuration guide for usage.
Mass-radius curves
| File | Location | Format | Source | Description |
|---|---|---|---|---|
massradiusEarthlikeRocky.txt |
data/mass_radius_curves/ |
Space-separated (M in M_Earth, R in R_Earth) | Zeng et al. (2019) | Earth-like rocky M-R curve (32.5% Fe + 67.5% MgSiO3) |
massradius_50percentH2O_300K_1mbar.txt |
data/mass_radius_curves/ |
Space-separated (M in M_Earth, R in R_Earth) | Zeng et al. (2019) | 50% H2O + 50% rocky M-R curve (300 K, 1 mbar surface) |
Radial profiles (validation benchmarks)
| File | Location | Format | Source | Description |
|---|---|---|---|---|
radiusdensitySeagerEarth.txt |
data/radial_profiles/ |
CSV (r, rho) | Seager et al. (2007) | Earth-like radial density profile |
radiusdensitySeagerEarthbymass.txt |
data/radial_profiles/ |
CSV (M, r, rho) | Seager et al. (2007) | Earth-like radial density profile (indexed by enclosed mass) |
radiusdensitySeagerwater.txt |
data/radial_profiles/ |
CSV (r, rho) | Seager et al. (2007) | Water-world radial density profile |
radiusdensitySeagerwaterbymass.txt |
data/radial_profiles/ |
CSV (M, r, rho) | Seager et al. (2007) | Water-world radial density profile (indexed by enclosed mass) |
radiusdensityEarthBoujibar.txt |
data/radial_profiles/ |
CSV (r, rho) | Boujibar et al. (2020) | Earth radial density profile |
radiuspressureEarthBoujibar.txt |
data/radial_profiles/ |
CSV (r, P) | Boujibar et al. (2020) | Earth radial pressure profile |
radiusdensityWagner.txt |
data/radial_profiles/ |
CSV (r, rho) | Wagner et al. (2011) | Earth radial density profile |
radiuspressureWagner.txt |
data/radial_profiles/ |
CSV (r, P) | Wagner et al. (2011) | Earth radial pressure profile |
radiusgravityWagner.txt |
data/radial_profiles/ |
CSV (r, g) | Wagner et al. (2011) | Earth radial gravity profile |
Seager EOS data
Zalmoxis uses the tabulated equation of state (EOS) data from Seager et al. (2007), providing pressure-density relations for iron (Fe epsilon), silicate (MgSiO3 perovskite), and water ice (ice VII). These tables are automatically downloaded during the setup process described above.
Analytic EOS
The Analytic:<material> EOS options (e.g. Analytic:iron, Analytic:MgSiO3, Analytic:H2O) use the analytic modified polytropic fits from Table 3 of Seager et al. (2007).
All fit parameters are hardcoded directly from the paper and do not require any data downloads.
This makes the analytic EOS suitable for quick setup, CI environments, or situations where the tabulated data files are not available.
Any EOS name beginning with Analytic: is self-contained.
The analytic EOS additionally provides three materials (MgFeSiO3, graphite, SiC) that are not available in the tabulated data.
Mass-radius relations
The following curves serve as reference benchmarks for validating the Zalmoxis internal structure model under both rocky and water-rich configurations:
-
Earth-like rocky exoplanets, with a 32.5% Fe + 67.5% MgSiO3 composition, taken from Zeng et al. (2019).
-
Water exoplanets, with a 50% Earth-like rocky core (32.5% Fe + 67.5% MgSiO3) and 50% water by mass at 300 K and surface pressure of 1 mbar.