{
 "tillotson_granite_params_cgs": {
  "la": 0.5,
  "lb": 1.3,
  "e_0": 160000000000.0,
  "rho_0": 2.7,
  "A": 180000000000.0,
  "B": 180000000000.0,
  "e_s": 35000000000.0,
  "e_s_prime": 180000000000.0,
  "alpha": 5.0,
  "beta": 5.0,
  "c_v": 7900000.0
 },
 "source": "Microphysics/EOS/tillotson/_parameters (Reinhardt & Stadel 2017 Table A1)",
 "energy_convention_finding": "actual_eos_init sets EOSData::mine = eos_e_0 = 1.6e11 erg/g. Initialising quiescent granite at e = eos_e_0 gives P(rho_0, e_0) = 49.7 GPa -- rock at rest under 50 GPa, wrong. Initialising with e measured from ZERO gives P(rho_0,0) = 0.00 GPa and reproduces the published granite Hugoniot (U_s = 3.68 + 1.24 u_p km/s) to a few percent. Use the e-from-zero convention; if the code floors e at eos_e_0, that floor must be lowered (eos.eos_e_0 runtime param) or the target will be pre-pressurised.",
 "hugoniot_e_from_zero": {
  "5": {
   "v_imp_kms": 5.0,
   "u_p_kms": 2.5,
   "U_s_kms": 6.365488262612243,
   "eta": 1.6467488270965631,
   "P_GPa": 42.967045772632645,
   "e_erg_g": 31250000000.0
  },
  "10": {
   "v_imp_kms": 10.0,
   "u_p_kms": 5.0,
   "U_s_kms": 10.225185356184513,
   "eta": 1.9569038530053324,
   "P_GPa": 138.04000230849093,
   "e_erg_g": 125000000000.0
  },
  "15": {
   "v_imp_kms": 15.0,
   "u_p_kms": 7.5,
   "U_s_kms": 13.95152406367227,
   "eta": 2.162516007997485,
   "P_GPa": 282.5183622893635,
   "e_erg_g": 281250000000.0
  },
  "18": {
   "v_imp_kms": 18.0,
   "u_p_kms": 9.0,
   "U_s_kms": 16.106872543455857,
   "eta": 2.266379823891364,
   "P_GPa": 391.39700280597737,
   "e_erg_g": 405000000000.0
  },
  "20": {
   "v_imp_kms": 20.0,
   "u_p_kms": 10.0,
   "U_s_kms": 17.51376468954798,
   "eta": 2.3308907602483875,
   "P_GPa": 472.87164661779553,
   "e_erg_g": 500000000000.0
  },
  "25": {
   "v_imp_kms": 25.0,
   "u_p_kms": 12.5,
   "U_s_kms": 20.94421819492513,
   "eta": 2.4803028192133096,
   "P_GPa": 706.8673640787232,
   "e_erg_g": 781250000000.0
  },
  "30": {
   "v_imp_kms": 30.0,
   "u_p_kms": 15.0,
   "U_s_kms": 24.278655012156534,
   "eta": 2.616613612678517,
   "P_GPa": 983.2855279923397,
   "e_erg_g": 1125000000000.0
  }
 },
 "hugoniot_e0_offset_WRONG": {
  "5": {
   "P_GPa": 85.10779557709778,
   "U_s_kms": 12.608562307718191
  },
  "10": {
   "P_GPa": 176.6492968717221,
   "U_s_kms": 13.085133101609044
  },
  "15": {
   "P_GPa": 316.68800546348535,
   "U_s_kms": 15.638913850048658
  },
  "18": {
   "P_GPa": 422.9889758763839,
   "U_s_kms": 17.406953739768884
  },
  "20": {
   "P_GPa": 502.8874276258031,
   "U_s_kms": 18.625460282437157
  },
  "25": {
   "P_GPa": 733.5352487743978,
   "U_s_kms": 21.734377741463636
  },
  "30": {
   "P_GPa": 1007.4019511452453,
   "U_s_kms": 24.874122250499884
  }
 },
 "impact": {
  "a_cm": 5000.0,
  "D_m": 100.0,
  "rho_i": 2.7,
  "rho_t": 2.7,
  "v_imp_cms": 1800000.0,
  "g_moon": 162.0,
  "m_i_g": 1414000000000.0,
  "KE_erg": 2.29e+24,
  "KE_Mt": 54.7
 },
 "crater_scaling_SUPERSEDED": {
  "K1": 0.095,
  "mu": 0.55,
  "nu": 0.4,
  "Y_dyncm2": 100000000.0,
  "pi2": 2.5e-07,
  "pi3": 1.143e-05,
  "regime": "strength-dominated",
  "pi3_over_pi2": 45.7,
  "V_tr_cm3": 457000000000000.0,
  "D_tr_km": 1.2,
  "D_final_km": 1.51,
  "source": "Holsapple 1993 pi-group scaling, competent rock",
  "why_wrong": "Two errors, both found by the Phase-0 theory track and independently reproduced here: (1) the K2=0.257 coefficient was omitted from the strength term of Holsapple (1993) eq.18, inflating it by 0.257^-1.275 = 5.65x and flipping the regime to 'strength'; (2) the paraboloid volume->diameter conversion dropped a factor sqrt(2), biasing D by 2^(1/6)=1.12. Also note pi3/pi2 = 45.7 is arithmetically correct but is NOT the regime test, because the strength term enters eq.18 raised to (2+mu)/2 = 1.275 and carries K2."
 },
 "crater_scaling": {
  "K1": 0.095,
  "K2": 0.257,
  "mu": 0.55,
  "nu": 0.33,
  "Y_dyncm2": 100000000.0,
  "gravity_term": 2.5e-07,
  "strength_term": 8.846e-08,
  "strength_over_gravity": 0.354,
  "regime": "gravity-dominated (marginal: strength supplies ~1/3 of resistance)",
  "V_tr_cm3": 764700000000000.0,
  "D_tr_km": 1.766,
  "d_tr_km": 0.62,
  "D_final_km": 2.21,
  "D_final_range_km": [
   1.52,
   2.35
  ],
  "Y_range_dyncm2": [
   10000000.0,
   1000000000.0
  ],
  "simple_to_complex_km": [
   15,
   20
  ],
  "crater_type": "simple",
  "source": "Holsapple 1993 eq.18 pi-group scaling; Hard Rock coefficients (K1=0.095, K2=0.257, mu=0.55, nu=0.33, Y=10 MPa) READ DIRECTLY from Singer et al. 2020 (arXiv:2008.10132) Table S1 p.41, fetched in-session as artifact singer2020.pdf. Paraboloid V=pi D^3/(16 sqrt2) with d/D=1/(2 sqrt2); D_final=1.25 D_tr. Note Table S1 tabulates Hard Rock with rho=3200 kg/m^3; rho_t=2.7 was used here for consistency with the Tillotson granite (-3.5% on D_final)."
 },
 "decay_law": {
  "P0_GPa": 391.4,
  "r_ic_over_a": 1.442,
  "n_central": 2.0,
  "n_bracket": [
   1.44,
   3.0
  ],
  "n_P97_at_18kms": 1.436,
  "n_Marcus2009_at_18kms": 1.375,
  "fit_window_r_over_a": [
   3,
   20
  ],
  "note": "n=1.44 with hemispherical isobars implies 130 impactor masses of melt = 108% of KE, so n=2 adopted as the shallowest exponent that closes the energy budget."
 },
 "melt_vapor": {
  "P_incipient_melt_GPa": 46,
  "P_complete_melt_GPa": 56,
  "P_incipient_vapor_GPa": 119,
  "P_complete_vapor_GPa": 579,
  "firmness": "46/56 GPa firm (P97 via Silber 2021); 119/579 GPa derived, order-of-magnitude",
  "predicted_melt_impactor_masses": 28,
  "predicted_melt_range": [
   10,
   30
  ],
  "predicted_vapor_impactor_masses": 9,
  "predicted_vapor_range": [
   5,
   10
  ],
  "complete_vaporisation": "none -- peak 391 GPa < 579 GPa threshold",
  "tillotson_trap": "EOS region boundary e_s maps to 47 GPa on this Hugoniot -- that is the MELT pressure, not vaporisation. Classify by the pressure thresholds, not EOS region flags."
 },
 "citation_provenance": {
  "audit_date": "2026-07-29 (superseded by a Phase-2 retrieval; corrected)",
  "why": "A reviewer flagged that several Author-Year citations in the Phase-0 theory artifacts were not traceable to any retrieval in the transcript. Audited: only Pierazzo et al. 1997 (doi 10.1006/icar.1997.5713) and arXiv:1603.08682 were actually FETCHED by the theory track. The rest were written from model knowledge. Verified afterwards via arXiv/Crossref:",
  "VERIFIED_real_and_on_topic": {
   "Pierazzo et al. 1997, Icarus, 'A Reevaluation of Impact Melt Production": "doi 10.1006/icar.1997.5713 -- fetched in-session AND independently confirmed",
   "Singer et al. 2020, arXiv:2008.10132": "FETCHED AND READ in-session by the Phase-2 track (artifact singer2020.pdf). Table S1 p.41 'Hard Rock' row confirmed by direct inspection: K1=0.095, K2=0.257, mu=0.55, nu=0.33, Y=10 MPa. This is now a retrieved source, not model knowledge.",
   "Kraus et al. 2012, JGR 117": "doi 10.1029/2012JE004082, 'Shock vaporization of silica and the thermodynamics of planetary impact events' -- real and directly on topic; not fetched"
  },
  "MISATTRIBUTED": {
   "'Chandnani et al. 2019' (doi 10.1029/2019JE006227)": "That DOI resolves to Davies et al. (2020), 'Silicate Melting and Vaporization During Rocky Planet Formation', JGR Planets -- a real paper but NOT the cited author/year. Treat the simple-to-complex transition corroboration as UNVERIFIED."
  },
  "UNVERIFIED_no_identifier_checked": [
   "Holsapple 1993, Annu. Rev. Earth Planet. Sci. 21, 333-373 (eq.18) -- the scaling law itself is standard and the arithmetic was verified numerically in-session, but the equation number and the coefficient values were not read from the paper",
   "Pike 1977 (via Lunar Sourcebook ch.4) -- simple-to-complex transition 15-20 km",
   "Silber et al. 2017 -- mean transition ~19 km",
   "Schmidt & Housen 1987; Holsapple & Schmidt 1982 -- underlying experiments"
  ],
  "status_of_the_numbers": "The PHYSICS numbers are unaffected: the Tillotson Hugoniot, the impedance-match peak state and the crater-scaling arithmetic were all computed and cross-checked numerically in-session (two independent implementations agreeing to <1e-3 %). What is uncertain is the literature PROVENANCE of the adopted coefficients and thresholds, not the algebra applied to them. Anyone reusing K1/K2/mu/nu or the melt/vapour thresholds should verify them against the primary source.",
  "CORRECTION": "My earlier audit was STALE. It was performed after Phase 0 and before Phase 2 finished. The Phase-2 analysis track subsequently FETCHED Singer et al. 2020 (saved as artifact singer2020.pdf, 88948987-00e9-4d1d-a476-6948d8d3ccca, 47 pp) and read the coefficient set directly from Table S1. I have now confirmed page 41 by eye: the 'Hard Rock' row reads K1=0.095, K2=0.257, mu=0.55, nu=0.33, Y=10 MPa -- exactly the values used. So the crater-scaling coefficients ARE sourced from a retrieved primary reference, and my statement that 'only Pierazzo 1997 was fetched' was wrong.",
  "target_density_nuance": "Singer Table S1 pairs the Hard Rock coefficient set with a target density of 3200 kg/m^3, whereas this study used rho_t = 2.7 g/cm^3 to stay self-consistent with the simulation's Tillotson granite (rho_0 = 2.7). Recomputed sensitivity: D_final = 2.21 km at rho_t=2.7 vs 2.13 km at 3.2, i.e. -3.5% -- well inside the 1.52-2.35 km spread from the order-of-magnitude uncertainty in Y, and 2.7 is the defensible choice. Recorded for transparency."
 }
}