A Self-Consistent Radiative-Hydrodynamic and Magnetocentrifugal Computational Framework for the Nuclear Impact Hypothesis: Global Sensitivity Analysis, Ejection Viability Mapping, and Falsifiability Tests

The Nuclear Impact Hypothesis (NIH) posits that a hypervelocity compact impactor (v0 ∼ 1500 km/s, m0 ∼ 10²⁴–10²⁶ g) striking a quiescent proto-Sun could trigger thermonuclear ignition while asymmetrically ejecting nascent planetary embryos through magnetocentrifugal thrust. This paper develops a self-consistent computational implementation of that hypothesis, in which (i) the impact velocity is treated in its correct dynamical regime (β0 = v0/c ≈ 5 × 10⁻³, γ0 − 1 ≈ 1.25 × 10⁻⁵), with the governing equations retained in general relativistic form solely so that the same formalism extends to genuinely relativistic progenitors; (ii) the impactor's mass, bulk density, and cross-section are enforced to be mutually consistent, r_imp = (3m0/4πρ0)^(1/3), rather than treated as independent free inputs; (iii) the magnetocentrifugal thrust term is given a toroidal radial dependence consistent with rotational field winding, Bφ(r) ∝ 1/r, together with a finite physical ceiling, removing the 1/ρ divergence that an unbounded thrust term would otherwise produce as the ambient density falls toward the photosphere; and (iv) all reported uncertainties come from executed variance-based sensitivity analysis and Monte Carlo forward propagation rather than from qualitative robustness claims.Under this formulation, the ingress sub-model gives a penetration depth δ = 0.067 R⊙ for the fiducial parameter set, and δ = 0.105 ± 0.091 R⊙ (90% interval [0.039, 0.186] R⊙) under parameter uncertainty. A Sobol-type global sensitivity analysis (Saltelli sampling, N = 128) shows δ is controlled predominantly by the electromagnetic charging-law parameters (S^T_B0 = 0.30, S^T_q = 0.33, S^T_βref = 0.57) rather than by the impact velocity itself (S^T_v0 ≈ 0). For the egress (ejection) sub-model, a 9 × 9 scan over the rotation rate ω and poloidal field B0 shows that escape from r0 = 0.1 R⊙ occurs in only ≈ 3.7% of the physically plausible grid, confined to a narrow corridor along which the surviving terminal velocities are v∞ ∼ 4.0–4.7 × 10³ km/s, substantially higher and more tightly localized in parameter space than a generic magnetocentrifugal-ejection scenario would suggest.All results reported here are produced by, and reproducible from, the Python code listed in the appendix; no numerical result in the text is asserted without a corresponding code path that generates it. The paper is explicitly framed as a hypothesis-testing exercise: its purpose is to expose the model's true sensitivities and failure modes so that the viability of the mechanism can be judged against future data, rather than to argue for its correctness.