Skinner, Pudritz, & Cloutier (2026) Mass-Radius Relations

Mass-Radius Curves from Skinner, Pudritz, & Cloutier (2026, second revision submitted to MNRAS) and planet profiles for every model planet as well as a validation sample of solar system objects.Formatting:See https://github.com/Bennett-Skinner/SkinnerPudritzCloutier2026-MR-curves to download non-zipped files (but full interior structure profiles are not provided).File Structure:Each zipped folder contains information for a specified zero albedo equilibrium temperature, with one compositional variable varied (Core mass fraction for CMF, water mass fraction for WMF, and H/He mass fraction for HHE). Once unzipping a folder, you will find many files with names formatted NAME_WX_Y_Z.txt, where NAME is either "many_radii", "inward_planet_profiles", or "outward_planet_profiles", W is the value of the varied compositional parameter, X is the name of the varied compositional parameter, Y is the value of the parameters for the rest of the composition (generally Earthlike), and Z is the zero albedo equilibrium temperature. If Z is omitted, it has Earth's temperature."many_radii" give the input parameters, masses, and radii of planets; "inward_planet_profiles" give inward profiles of the planetary interior, and "outward_planet_profiles" give outward profiles of the planetary interior (see below for column-by-column descriptions). The full planet profile is the concatenation of the inward and outward profiles, although note that each thread includes the fitting mass and thus plots may look strange if one of the two datapoints at the same mass is not removed.The earth-like isocomposition curve is found under WMF curves as a WMF of 0. In addition, one zipped folder contains all the validation profiles, with the formatting following NAME_OBJECTNAME_INFO, where NAME is the same as above, OBJECTNAME is the name of the object, and INFO is a string describing any specific attributes of the modelled planet, corresponding to discussion in Skinner, Pudritz, & Cloutier (2026).Info (many_radii):The columns are as follows:mass - Mass of the planet (Earth radii)envelope_frac - Envelope mass fraction of the planetwater_frac - Water mass fraction of the planetmantle_frac - Mantle mass fraction of the planetcore_frac - Core mass fraction of the planetH_frac_envelope - H fraction of the envelope, by massHe_frac_envelope - He fraction of the envelope, by massMgO_frac_mantle - MgO fraction of the mantle, by numberSiO2_frac_mantle - SiO2 fraction of the mantle, by numberFeO_frac_mantle - FeO fraction of the mantle, by numberFe_frac_solid_core - Atomic Fe fraction of the solid core, by numberFeS_frac_solid_core - FeS fraction of the solid core, by numberFeO_frac_solid_core - FeO fraction of the solid core, by numberFe_frac_liquid_core - Atomic Fe fraction of the liquid core, by numberFeS_frac_liquid_core - FeS fraction of the liquid core, by numberFeO_frac_liquid_core - FeO fraction of the liquid core, by numberT_eq - Zero-albedo equilibrium temoerature (if there was no internal luminosity) (K)albedo - Bond Albedorotation_period - Rotation period (days)age - Age (Gyr)atmospheric_Z - Metallicity in atmosphere used to calulate opacity (log solar)100pa_radius - Radius at a pressure of 100 Pa (Earth radii)2000pa_radius - Radius at a pressure of 2000 Pa (Earth radii)transit_radius - Radius at an optical depth of 2/3 for a transiting ray of light (Earth radii) (THIS IS PROBABLY THE RADIUS YOU WANT, but read section 2.8)twothirdsopt_radius - Radius at an optical depth of 2/3 for a ray of light travelling directly to the center of the planet (Earth radii)condensed_radius - Layer at which planet is made of condensed matter (Earth radii)diff_v_2000pa - (transit_radius-2000pa_radius)/transit_radius\*100diff_v_100pa - (transit_radius-100pa_radius)/transit_radius\*100inner_pressure - Pressure at center of model planet (GPa)inner_temperature - Temperature at center of model planet (K)Info (inward_planet_profiles & outward_planet_profiles):The columns are as follows:m[Frac_of_total] - Mass as a fraction of total planetary radiusr[R_e] - Radius in Earth radiiP[GPa] - Pressure in GigapascalsT[K] - Temperature in KelvinOptical_Depth[unitless] - Optical depth, as measured along a line directly going to the planetary coredensity[kg/m3] - Density in kg/m3ad_grad[unitless] - Unitless adiabatic gradient, defined as adgrad = P/T\*dT/dPopacity[m2/kg] - Opacity in m2/kgopacity_density[kg/m3] - The density from Freedman+2014 (if not in atmosphere, just the density). Note this is only present for inward planet profiles, not outward planet profiles.MgSiO3_melting[K] - Melting temeprature of mantleFe_melting[K] - Melting temperature of corephase[see_doc] - Phase of material, see below for guide

Codes for phases:H/He (-11):All H/He: -11 (Everything done in one table, phase transitions not exposed to planetary solver)

Water (-10 to 5) (From AQUA):Ice-Ih: -1Ice-II: -2Ice-III: -3Ice-V: -5Ice-VI: -6Ice-VII: -7Ice-X: -10        Vapour: 3Liquid: 4Supercritical/superionic: 5

Mantle (6 to 9):Perple_X: 6Post-Perovskite and High-Pressure Phases (Mg2SiO4 and MgSi2O5) : 7Dissociated: 8Liquid: 9

Iron (10 to 11):Solid: 10Liquid: 11Changelog for v1.0:Initial version. This was published alongside the second version of the arXiv posting.Changelog for v1.1.:No changes impact the mass-radius curves themselves, but several minor points crucial for interpretation in the other columns have been changed. This is the version that accompanies the accepted paper.- Core mass fractions column fixed (previously did not subtract off H/He mass fraction)- Phase of iron now always displayed consistently (previously smoothing prescription led to the phase of iron listed in the last column not always being consistent with whether Iron was above its melting temperature)- For condensed planets, 100 Pascal radius = 2000 Pascal radius = 2/3 Transit radius = 2/3 Inwards radius now strictly enforced (previously, the columns containing the percentage differences were low but nonzero values due to internal details of the solver, even though the listed radii were identical)- Saturn's listed atmospheric metallicity fixed to value used (was previously the default 50)- New Earth-like quick access folder (completely redundant with the 0% water mass fraction varied lines but allows for much faster finding of most desirable mass-radius curve)- Junk files intended for internal use removed- File names changed to remove spaces for better compatibility with most OSPlease direct all correspondance to [e-mail] and cite Skinner, Pudritz, & Cloutier (2026, in press.) if using this work. Mass-radius curves with any of the listed free parameters varied can be generated upon request. Feedback is welcome.