heliobench
Solar System Navigation Physics Boundary Benchmark v2.0
太阳系全域导航物理边界基准计算系统
A production-grade Rust library providing authoritative physical boundary calculations across the full extent of the solar system — from the solar photosphere to the theoretical limit of the Oort Cloud — built on three authoritative standards: CODATA 2018, IAU 2015 Resolution B3, and JPL DE440 ephemeris.
What This Library Does
Most solar system models treat physical boundaries as lookup tables or hardcoded constants. This library derives every boundary from first principles, using the actual physical equations that govern each zone. The termination shock distance is not a number — it is the solution to a pressure balance equation. The Alfvén surface is not a range — it is the unique root where solar wind velocity equals local Alfvén velocity, solved by bisection under PSP-calibrated coronal parameters. Ice line positions are not empirical fits — they emerge from radiation equilibrium between stellar flux and sublimation thermodynamics.
Seven-Zone Architecture
Zone 01 — Solar Atmosphere & Alfvén Surface (1–25 R☉)
Saito et al. (1977) multi-component coronal magnetic field replacing the naive 1/r² monopole approximation. Leblanc et al. (1998) two-component electron density model (A/r² + B/r⁴ + C/r⁶) replacing constant-speed wind assumption. Parker spiral magnetic field with source surface boundary at 2.5 R☉. PSP-calibrated Alfvén surface: solar minimum 12 R☉, ascending phase 19 R☉, solar maximum 15.8 R☉, consistent with Kasper et al. (2021) first-crossing at 18.8 R☉. Physical coronal heating: Poynting-flux magnetic reconnection (Sweet-Parker / Petschek), Alfvén wave dissipation from amplitude and damping scale, heating-radiation equilibrium corona temperature 1.3–2.2 × 10⁶ K.
Zone 02 — Inner Planetary Belt (0.05–2.7 AU)
Vis-viva equation, Keplerian orbital mechanics with GM_SUN = 1.32712440018 × 10²⁰ m³s⁻² (JPL DE440). Full planet database: Mercury, Venus, Earth, Mars with surface pressure, magnetic field, obliquity. Synodic periods, sphere of influence, Hill radius. Eccentricity vector in correct 3D cross-product formulation.
Zone 03 — Main Asteroid Belt (2.7–3.6 AU)
Kirkwood resonance positions with corrected period ratio definition P_body/P_Jupiter. Tisserand parameter with full eccentricity term: T = a_J/a + 2√(a/a_J × (1−e²)) cos(i). Jupiter resonances: 4:1 at 2.065 AU, 3:1 at 2.501 AU, 5:2 at 2.825 AU, 2:1 at 3.278 AU, 3:2 Hilda at 3.971 AU.
Zone 04 — Giant Planet Region (3.6–30 AU)
Chapman-Ferraro magnetopause standoff with distance-scaled solar wind density. Jupiter magnetosphere cross-validated against Juno (Rutala 2025, JGR): median 71±24 R_J, compressed 51±23 R_J, Joy (2002) dual-mode distribution 63±4 / 92±6 R_J. Saturn cross-validated against Cassini (Achilleos 2008, Jackman 2019): typical 18–25 R_S. Laplace resonance verification for Io/Europa/Ganymede. Neptune-Pluto 3:2 resonance. Giant planet Hill radii.
Zone 05 — Kuiper Belt & Scattered Disk (30–1500 AU)
Neptune resonance positions with corrected period ratios: Plutinos 3:2 at 39.4 AU, 2:1 at 47.7 AU, 4:1 at 75.8 AU. TNO dynamical classification: Plutino, Cubewano, Scattered Disk, Detached, Extreme TNO.
Zone 06 — Oort Cloud (2000–200000 AU)
Solar tidal radius from Galactic mass density (Eilers et al. 2019, ρ = 0.10 M☉/pc³): r_tidal = (M☉/4πρ)^(1/3) = 191,143 AU = 3.02 light-years. Dynamical stability boundary at 1/3 Hill radius: 63,714 AU = 1.01 light-years. Density sensitivity analysis: ρ = 0.08–0.30 M☉/pc³ → 132,531–205,903 AU. Escape velocities from 1000 to 50,000 AU.
Zone 07 — Heliosphere (80–200 AU)
LISM parameters from Voyager in-situ measurements: ion density 0.06 cm⁻³, neutral hydrogen 0.18 cm⁻³, magnetic field 3.2 μG, ISM velocity 26.3 km/s, temperature 8000 K. Three-component ISM pressure: thermal + magnetic (P = B²/2μ₀) + ram. Termination shock via pressure-balance bisection: 94.1 AU (vs Voyager 1: 94.0 AU, Voyager 2: 83.7 AU). Heliopause via total heliosheath pressure bisection: 126.6 AU (vs Voyager 1: 121.6 AU, Voyager 2: 119.0 AU). IMF-derived heliosheath magnetic field from shock compression ratio.
Navigation & Communication
Relativistic Doppler shift with classical/relativistic comparison at 1 km/s (5.56 ppt correction) and 100 km/s (55.6 ppb correction). Free-space path loss (FSPL) in dB for X-band (8.4 GHz): 274.4 dB at 1 AU, 304.0 dB at 30 AU. Light-time delays for all eight planets. Stellar habitability comparison across G2V, K0V, K5V, M0V, M2V, M5V spectral types.
Standards & Validation
| Standard | Application |
|---|---|
| CODATA 2018 | G, c, k_B, μ₀, m_p, m_e, e |
| IAU 2015 Resolution B3 | M☉, R☉, L☉, T_eff = 5778 K |
| JPL DE440 | GM_SUN, planetary masses |
| Saito et al. (1977) | Coronal magnetic field |
| Leblanc et al. (1998) | Coronal electron density |
| Kasper et al. (2021) | Alfvén surface PSP calibration |
| Eilers et al. (2019) | Galactic mass density |
| Voyager 1 & 2 in-situ | Heliosphere boundary validation |
| Juno / Cassini | Giant planet magnetosphere validation |
65 unit tests, all passing.
Part of the Soulweaver Embodied AI Infrastructure Stack
This library is one component of a broader Rust-native infrastructure stack for embodied intelligence systems, published under the Embodied AI Memory Infrastructure Stack dataset on IEEE DataPort.