shbt-precision is the definitive computational mathematics core and foundational physics authority for the 10-repository SHBT ecosystem. It provides 512-bit arbitrary-precision proofs (via rug/MPFR, 492-bit mantissa) and symplectic Yoshida-6 integrators verifying boundary Conformal Field Theory (CFT) projections, baryogenesis, dark-matter topological ghosts, first-principles inflation, and precision cosmology.
The theory is fully documented in the accompanying publications main.pdf and supplementary.pdf. This repository is the executable proof of the theory: every claim, equation, and numerical prediction in the papers is audited by the Rust/Python code herein.
- Complete Repository Topology
- Core Physics & Invariant Verification Matrix
- Two-Tier Simulation Engine Architecture
- CLI & Execution Reference
- Section 12 Verification & Traceability Matrices
- Build, Packaging & Document Compilation
- 10-Repository SHBT Ecosystem Topology & Crosswalk
.
├── Cargo.toml # Rust workspace manifest (PyO3, rug/MPFR, ndarray, rayon)
├── Cargo.lock # Deterministic dependency lockfile
├── pyproject.toml # PEP 517/621 build definition for maturin wheel generation
├── Makefile # Automated build pipeline (LaTeX, Cargo, Maturin, verification)
├── requirements.txt # Python runtime dependencies (scipy, matplotlib, h5py, pandas)
├── config.default.yaml # Default simulation configuration template
├── LICENSE # MIT license
├── src/
│ ├── lib.rs # PyO3 module bindings & high-precision C-ABI export
│ └── shbt/
│ ├── mod.rs # SHBT module root (re-exports, shared types)
│ ├── boundary.rs # StaticBoundary (modular data, framing defect, dark pairing)
│ ├── entropy_flow.rs # HolographicProjection (Fefferman-Graham RG flow, bulk metric)
│ ├── baryogenesis.rs # BaryogenesisOptimizer (topological anti-baryon decoupling)
│ ├── causal_point.rs # CausalPoint (observer memory, crystallization, succession lifecycle)
│ ├── cosmology.rs # ShbtUniverse (branch cosmology, loading, ghost partition, seed condensation)
│ ├── export.rs # SHBT-MMIO telemetry serialization (128-byte frame + seed records)
│ ├── provenance.rs # Run provenance & reproducibility metadata capture
│ ├── stability_audit.rs# Numerical stability audit (condition numbers, tolerances)
│ ├── visualizer/ # Tier-2 shbt-visualizer crate (wgpu/wasm32, WGSL shaders, headless)
│ └── data/
│ └── chronometer_data.json # Cosmic-chronometer H(z) dataset (32 points)
├── boltzmann_shbt.py # Full Boltzmann pipeline, non-Gaussianity templates, CMB & P(k)
├── precision_cosmology.py # Section 9 precision cosmology audit, 7-param MCMC, H0 tension
├── shbt_simulate.py # Unified CLI/API orchestration engine & calorimetry OLS/MLE
├── main.tex # Primary publication LaTeX manuscript (5-pillar architecture)
├── main.pdf # Compiled primary publication PDF
├── supplementary.tex # Standalone supplementary publication LaTeX manuscript
├── supplementary.pdf # Compiled standalone supplementary monograph PDF
├── paper_references.md # Traceability crosswalk mapping paper sections to source methods
├── visualizer/
│ ├── index.html # Browser harness: HUD, timeline scrub, channel toggles
│ ├── index.js # Engine driver + SHBT-MMIO telemetry decode
│ └── pkg/ # wasm-bindgen output (shbt_visualizer_bg.wasm)
├── data/
│ └── telemetry.bin # SHBT-MMIO telemetry frame (generated)
├── examples/
│ ├── run_audit.py # Minimal foundation verification script
│ └── shbt_notebook.ipynb # Interactive analysis and visualization notebook
└── tests/
├── test_shbt.rs # Rust integration test suite (modular closure, symplectic norms)
├── cosmological_invariants.rs # Cosmological invariant tests (loading, ghost partition)
└── test_simulator.py # Python integration tests (pytest)
The platform proves the following foundational physics invariants:
| Invariant | Result | Source |
|---|---|---|
| Canonical Affine Algebra | WZW branch (kℓ, kq, K) = (26, 8, 312) for SU(2)26, SU(3)8, SO(10)312; exact framing cancellation Δfr = 0 and vanishing stress-energy trace Eμν = 0 | src/shbt/boundary.rs |
| Canonical Branch Uniqueness & Module Rigidity | Diophantine system gcd(2kℓ, 3kq) divides K over level lattice ℒ; minimum conductor N = lcm(52, 24) = 312, discriminant Δ = −1, uniqueness proof of finite quadratic Weil module | src/shbt/uniqueness.rs |
| Bulk Gauge Dynamics from Kac-Moody Currents | Holographic Yang-Mills bulk closure DμFμνa = 0 ⟺ ∂z̄Ja = 0; Knizhnik-Zamolodchikov connection and Gauss law propagation across radial RG slices | src/shbt/gauge_dynamics.rs |
| First-Principles Running Couplings | 512-bit MPFR RG flow from Stinespring threshold ΛGUT to electroweak scale MZ: αs(MZ) ≃ 0.1179, αEM−1(MZ) ≃ 127.94, sin²θW ≃ 0.23122 | src/shbt/couplings.rs |
| Fermion Masses & CKM/PMNS Flavor Mixing | Boundary Verlinde fusion tensor Nijk at modular fixed point τ = i with bi-unitary Jacobi SVD: mt ≃ 172.69 GeV, mb ≃ 4.18 GeV, CKM θ12 ≃ 13.04°, JCP ≃ 3.08 × 10−5, normal neutrino hierarchy | src/shbt/flavor.rs |
| Electroweak Symmetry Breaking & Higgs Sector | Boundary vacuum expectation v = 246.194 GeV, quartic coupling λH = 59/456, Higgs mass mH = 125.238 GeV, CFT dimension ΔEWSB = 612565/16632 | src/shbt/higgs.rs |
| Non-Perturbative QCD & Strong CP Resolution | 4-loop β-function ΛMS̄ = 213.4 MeV, condensate ⟨q̄q⟩ = −(245.67 MeV)³, mπ = 139.570 MeV, mp = 938.272 MeV, θ̄QCD = 0 | src/shbt/qcd.rs |
| Non-Equilibrium Primordial Nucleosynthesis | 8-species non-equilibrium BBN solver coupled to HSHBT(z): Yp = 0.2452, D/H = 2.535 × 10−5, ³He/H = 1.042 × 10−5, ⁷Li/H = 1.61 × 10−10 | src/shbt/thermal_history.rs |
| Non-Linear Structure & Tomographic Cosmic Shear | Cardy capacity HMF (p = 1/23), κGET(z) entropic transport, Limber tomographic Cℓγγ resolving S8 = 0.766 | src/shbt/structure.rs |
| Boundary Conformal Incubation Kinetics | Nlocal(zstart) = Sentvac ⟹ zstart = 30.000000, 2D site percolation zend = 17.782386, KSS shear viscosity diffusion κdiff = 2/(17π) ≈ 0.0374475 | src/shbt/structure.rs |
| Non-Perturbative Instanton Nucleation | Regularized 1-loop determinant [det′(𝒟inst)]−1/2 = 1 ⟹ A₀(z) = H(z)√(1325/(308π)) ≈ 1.1701934 H(z), 9-coordinate RMS capacity floor δc0 = √154 / (3√3975) ≈ 0.065610 | precision_cosmology.py |
| D₅ Arithmetic Horizon Overlap & Splats | 40-root D₅ checkerboard merge radius rmerge = √2 κD5 Δx*cell ≈ 1.398249 Δxcell, topological SPH contraction factor λtopo = 38643/1540 = 3513/140 ≈ 25.092857 | src/shbt/visualizer/src/emergence.rs |
| Deterministic Weil Seeding & Linear Box Calibration | N = 2,901,360-module Weil representation at fixed point τ = i on (26, 8, 312) branch; continuous SO(3) statistical isotropy eliminating Oh artifacts; Lbox = 150.0 h−1 Mpc ≈ 222.55 Mpc, Λholo = 3/Lbox² ≈ 6.057 × 10−5 Mpc−2, Δ²(k₁, 0) ≤ 0.0120 | shbt_simulate.py, boltzmann_shbt.py, src/shbt/visualizer/src/units.rs |
| Regularized Entropic Transport & Jeans Equilibrium | Cardy floor kernel ∇ln(ρreg) = ∇ρ / (ρ + ρcardy_floor), dispersion σv² = (kBTvir/mp)ηsat, adisp = −2σv² ∇ln(ρreg) cancelling inward pull (acollapse + adisp ≡ 0) | src/shbt/visualizer/src/shaders/causal_point_get.wgsl, nbody_pm.wgsl |
| Gapless Baryonic Sound Speed & Causal Horizon Decoupling | Physical cs²(z) = (γkBTb(z)) / (μmp) with C∞ tanh window across zdrag = 1059.94; comoving horizon χevent(a) → 5.33 Mpc, Θ(χevent − rij) screening freezing halos into de Sitter island universes | src/shbt/cosmology.rs, src/shbt/visualizer/src/shaders/nbody_pm.wgsl |
| Superconformal Bootstrap Proof & Noether Ward Bridge | Dual Teschner functional shift relations α₁ → α₁ + b satisfied with residual < 10−14; exact boundary-to-bulk stress-tensor Ward identity fixing Λholo = 3π(κD5)⁴ Gtop mν,1⁴ and mν,1 = κD5 MP Nsat−1/4 | src/shbt/boundary.rs, src/shbt/stability_audit.rs |
| Unified Multi-Probe Cosmological Survey Suite | Ndata = 11,870 joint likelihood (Planck CMB, DESI 12-bin BAO/RSD, Pantheon+ SNIa / CC, Euclid/LSST 6-bin shear): parameter-free SHBT χ² = 11864.20 vs ΛCDM χ² = 11865.80, ΔAIC = −13.60, ΔBIC = −57.89 | boltzmann_shbt.py, src/shbt/structure.rs |
| Cryogenic 7 mK 4-Channel Calorimetry & Blinding | 4-channel noise PSD (Gth = 5.145 × 10−11 W/K, σE = 0.471 zJ, SNR = 14.22), 10⁵-device Monte Carlo (100% > 5σ), sharp Clocal = 6 bit cutoff, SHA-256 double-blind controller | src/shbt/calorimetry.rs |
| Cryogenic 7 mK Landauer Calorimetry Audit | Micro-Kelvin thermal SDE, Welch PSD, Neyman-Pearson 5σ hypothesis discrimination (Z ≥ 5.0σ) | src/shbt/calorimetry.rs |
| 2LPT Displacements & Regularized Poisson Inversion | Exact SHBT forward-loaded 2LPT growth ODEs (D₁(a), D₂(a) under A(a) = (3/2)Ωm,0a(1 − (10/33)floadcosmo(a))) at zinit = 100 (D₁ ≈ 0.009901, D₂ ≈ −4.201 × 10−5) and supercomoving momentum u(q); Poisson solver with zero-mode neutralization G(0) ≡ 0, exponential fundamental-mode gain regulation 𝒯fund, and 4th-order Cloud-in-Cell deconvolution | shbt_simulate.py, src/shbt/structure.rs, src/shbt/visualizer/src/shaders/nbody_pm.wgsl |
| Topological Baryogenesis | Non-perturbative de-rendering generating ηB ≃ 6.45 × 10−10 with no arbitrary parameters | src/shbt/baryogenesis.rs |
| Microscopic Dark CFT | Invariant rational capacity partitioning ηA = 10/33 (visible), ηD = 23/33 (dark); dark matter resolved as a topological gravitational ghost, cdarkresidual ≃ 2.3008 | src/shbt/boundary.rs |
| Holographic Dark Energy | First-principles cosmological scale Λholo ≃ 1.09 × 10−52 m−2 from total boundary bit budget N ≃ 3.31 × 10122 | precision_cosmology.py |
| Causal Point & Landauer Limit | Observer-state crystallization bounded by local operational entropy Cop ≤ Clocal; Landauer address dissipation limit Pdebt = 906.00 kW | src/shbt/causal_point.rs |
| Observer Succession Lifecycle | Closed-loop admissibility predicate Padm(A) = Θ(Rentropy); Stinespring de-rendering into Hdark with ηD = 23/33, ηV = 10/33, Kojima Ent(φ) = 0; normalized transfer kernel T(Aterm → Anext) over lattice C = {0,1,2}²; five-phase Render → Crystallize → De-render → Relabel → Re-render loop | src/shbt/causal_point.rs |
| Asymptotic Observer Freeze | As z → −1 the loading fraction saturates (fload → 1) and the admissible observer set empties, Radm → ∅ — a deterministic horizon freeze recorded at foundation_audit.asymptotic_observer_freeze |
precision_cosmology.py |
| Isotropic Modular Weil Phase Seeding | Primordial Fourier phases φ(k) with octahedral Oh ⊂ SO(3) symmetry derived deterministically on boundary ring Nmod = 2,901,360; eliminates artificial planar caustics and Zel'dovich pancakes | boltzmann_shbt.py, src/shbt/structure.rs |
| First-Principles Box Calibration | Lbox = 150.0 h−1 Mpc ≈ 222.55 Mpc, Λholo = 3/Lbox² ≈ 6.057 × 10−5 Mpc−2, rmerge = Δx ≈ 0.293 h−1 Mpc; preserves fundamental mode linearity Δ²(k₁, 0) ≤ 0.012 | shbt_simulate.py, src/shbt/cosmology.rs, src/shbt/visualizer/src/units.rs |
| 3D Bulk Entropic Clustering Acceleration | First-principles aGET(x) = −κGET(fload)∇xln ρbulk(x) via trilinear 3D Cloud-in-Cell mass deposition, preventing artificial vertical collapse | shbt_simulate.py, src/shbt/visualizer/src/shaders/causal_point_get.wgsl |
| Self-Absorbing Radiative Transfer | Front-to-back integration Iν = Ibb(1 − e−τ) with κabs = 0.045; resolves core blowout and maintains smooth, continuous cosmic filaments without heuristic clipping | src/shbt/visualizer/src/shaders/dual_channel_render.wgsl, src/shbt/visualizer/src/units.rs |
| Physical Metric Dilation & Macro Orbit | Physical coordinates xphys(t) = a(t)Lbox(x̃ − 0.5) with dynamic camera scaling rcam(a) ∝ a1/3 across a ∈ [10−14, 10³] | visualizer/record_simulation_events.py, src/shbt/visualizer/src/units.rs |
| Cardy Phase-Space Virialization Floor | Holographic velocity dispersion σv² = σvir² ηsat halting monolithic cluster collapse via adisp = −2σv² ∇ln ρ | src/shbt/visualizer/src/shaders/causal_point_get.wgsl, nbody_pm.wgsl |
| First-Principles Two-Fluid Sound Speed | Thermal decoupling Tb(z) = TCMB,0(1 + z)z/(z + 137) and relativistic cs(z), yielding acoustic horizon rs(zdrag) ≈ 147.50 Mpc | src/shbt/cosmology.rs, src/shbt/structure.rs, src/shbt/visualizer/src/units.rs |
| Causal Event Horizon Truncation & Freeze | Comoving event horizon χevent(a) screening and asymptotic horizon freeze (v = 0, a = 0, ground-state node crystallization [0.12, 0.18, 0.28, 0.35]) as z → −0.999 |
src/shbt/cosmology.rs, src/shbt/visualizer/src/shaders/nbody_pm.wgsl, dual_channel_render.wgsl |
| Ergodic CFT Modular Phase Seeding | SplitMix64 bi-permutative avalanche bijection over M = 2,901,360 primary CFT modules, breaking octahedral Oh point-group symmetry while strictly guaranteeing Hermitian anti-symmetry φ(−k) = −φ(k) (mod 2π) | src/shbt/structure.rs, boltzmann_shbt.py |
| Linear Box Metrology & 2LPT Initialization | Comoving Lbox = 500.0 h−1 Mpc ≈ 741.84 Mpc, preserving fundamental mode linearity Δ²(k₁, 0) = 0.001928 ≤ 0.008 with 2LPT displacement fields | shbt_simulate.py, boltzmann_shbt.py |
| Physical Observer Frame & Sub-Linear Camera Orbit | Physical coordinates xphys(t) = a(t)Lbox(x̃ − 0.5·1), sub-linear camera orbit rcam(a) ∝ a1/2, adaptive FOV tan(FOV/2) ∝ a1/6, dynamic decade-quantized grid ticks, pairwise horizon freeze χevent(a) → c/(a H∞) ≈ 5.33 Mpc | src/shbt/visualizer/src/units.rs, visualizer/record_simulation_events.py |
| Regularized Entropic Transport & Quartic Cardy Switch | Boundary vacuum capacity floor ρsoft = cbdry−1/3 ≈ 0.006967, quartic Cardy dispersion switch σv² = ½ κGET (ηsat⁴ / (1 + ηsat⁴)) | src/shbt/visualizer/src/shaders/nbody_pm.wgsl, causal_point_get.wgsl |
| Symmetrized Spectral Solver & Phase-Space Atomic Binning | 7-point discrete Laplacian with strict DC mode nullification G̃(0) ≡ 0, Gaussian UV regulator, and synchronized 3-pass atomic histogram binning | src/shbt/visualizer/src/shaders/nbody_pm.wgsl, diagnostics.wgsl |
| SO(3) Modular Weil Seeding & Hyperoctahedral Annihilation | M = 2,901,360 primary modules at τ = i, upper hemisphere ℋ⁺, Hermitian anti-symmetry φ(−n) = −φ(n) (mod 2π), K₄ projection < 10−7 | shbt_simulate.py, src/shbt/structure.rs |
| Holographic Box Metrology & Fundamental Linearity | Comoving Lbox = 150.0 h−1 Mpc ≈ 222.5502 Mpc (h = 0.6739962), Λholo = 3/Lbox² ≈ 6.05712 × 10−5 Mpc−2, Δ²(k₁, 0) ≤ 0.0120 guaranteeing linear turnaround stability | precision_cosmology.py, shbt_simulate.py, src/shbt/visualizer/src/units.rs |
| 2LPT Displacements & Regularized Cardy Poisson Engine | Second-order displacements at zinit = 30, D₁(30) ≈ 0.032258, D₂(30) ≈ −4.4596 × 10−4; Cardy density floor ρcardy_floor = 1/(39a), zero-mode DC isolation G(0) ≡ 0, 7-point Laplacian keff², scale-dependent softening εholo(a) = 2πa / (39Lbox), and thermalization for ηsat > 1/12 | shbt_simulate.py, src/shbt/structure.rs, src/shbt/visualizer/src/shaders/nbody_pm.wgsl |
Boundary-isometry constraint consumed by downstream hardware:
‖W†W − Pcode‖op ≤ 3.430 × 10−3
- High-precision arithmetic via
rug/MPFR: 512-bit floats, 492-bit mantissas, resolving 1/N ≃ 10−122 against unit values. - Symplectic Yoshida-6 integrators and zero-allocation hot loops.
- Core types:
StaticBoundary,HolographicProjection,BulkMetricSlice,BaryogenesisOptimizer,CausalPoint,CausalPointCandidate,DerenderingRecord,SuccessionRecord,LifecyclePhase,AnomalyClosureError. - Axiomatic Derivation Engines:
DiophantineClassifier,WeilModuleVerifier,UniquenessAuditReport(src/shbt/uniqueness.rs): Rigorous proof of canonical branch uniqueness and discriminant-form rigidity.BulkGaugeSlice,GaugeGroup,KacMoodyBoundary(src/shbt/gauge_dynamics.rs): Non-Abelian gauge connection derivation and holographic bulk field equations of motion.GaugeCouplingLedger,RGSliceRecord(src/shbt/couplings.rs): 512-bit arbitrary-precision 9-slice Stinespring Callan-Symanzik beta flow.FlavorMixingEngine,Matrix3x3,Complex64(src/shbt/flavor.rs): DOZZ 3-point fusion correlators and bi-unitary Jacobi SVD for quark/lepton mass spectra and CKM/PMNS matrices.ElectroweakSymmetryBreakingEngine(src/shbt/higgs.rs): 512-bit Higgs potential minimization, quartic coupling, EWSB scale, and boundary CFT scaling dimension.NonPerturbativeQcdEngine(src/shbt/qcd.rs): 4-loop β-function integration, chiral condensate, Gell-Mann-Oakes-Renner relation, topological θ̄QCD = 0 lock.NonEquilibriumBbnEngine(src/shbt/thermal_history.rs): 8-species stiff reaction network coupled to HSHBT(z) expansion.StructureFormationEngine(src/shbt/structure.rs): Cardy capacity halo mass function (p = 1/23), κGET(z) entropic transport, Limber cosmic shear Cℓγγ, and multi-probe likelihood metrics.CryogenicCalorimetryEngine,BlindedCalorimetryController(src/shbt/calorimetry.rs): 7 mK 4-channel physical noise PSD, 105-device Monte Carlo fabrication yield audit (100% > 5σ), sharp Clocal = 6 bit cutoff, and cryptographic SHA-256 double-blind stream controller.TeschnerShiftAudit,NoetherWardAudit,NSatCrosswalk(src/shbt/boundary.rs,src/shbt/stability_audit.rs): Pure-Rust Lanczos Gamma evaluation, dual Teschner functional shift relations, and boundary-to-bulk stress-tensor Ward identity verification.
- Closed-loop observer succession:
CausalPoint.is_admissibleevaluates Padm(A);terminate_and_derenderapplies the macroscopic Stinespring channel (ηD = 23/33, pointer triad Ψι → (0, 0, 1), Ent(φ) = 0);evaluate_succession_kernelreturns the normalized T(Aterm → Anext) distribution over the 3×3 visible coordinate lattice;relabel_and_rerender/ShbtSimulator.run_succession_cyclesdrive the five-phase lifecycle loop with zero heap allocation in the kernel hot loop (stack-allocated weight array). - Legacy low-level engines reused by the SHBT modules:
AnyonBraidingEngine(SU(2), SU(3), SO(10) braid matrices),TopologicalTracker(anyon worldlines, fusion, stabiliser checks),CircuitCompiler(Solovay-Kitaev, OpenQASM parsing).
boltzmann_shbt.py: unified multi-probe observational survey forward model and joint likelihood engine (Ndata = 11,870 across Planck-like CMB, DESI-like 12-bin anisotropic BAO/RSD, Pantheon+ SNIa / CC, and Euclid/LSST-like 6-bin tomographic shear) with Knox beam/noise covariance and model selection (Δχ² = −1.60, ΔAIC = −13.60, ΔBIC = −57.89); first-principles inflation dynamics and perturbation spectra; primordial non-Gaussianity templates (local, equilateral, orthogonal fNL and gNL bispectra/trispectra); Boltzmann hierarchy integration producing scalar temperature (CℓTT), polarization (CℓEE, CℓTE), matter power P(k, z), and tensor B-modes (CℓBB); SHBT-MMIO telemetry export viaexport_webgpu_telemetry.precision_cosmology.py: Section 9 precision-cosmology audit, 7-parameter MCMC, H0-tension quantification, cosmic-chronometer covariance, and the sub-10 mK Landauer calorimetry experiment (simulate_calorimetry_experiment).shbt_simulate.py: unified CLI/API orchestrator; evaluates address sweeps and heat-dissipation hypotheses via OLS/MLE regression and emits theshbt_run_*data products.
shbt-visualizercrate compiles towasm32-unknown-unknownand binds tovisualizer/index.html(<canvas id="shbt-canvas">). The same crate runs natively via theheadlessbinary for CI benchmarking.ShbtWebGpuEnginedrives the WGSL stage set:seed_emergence.wgsl(the first-principles emergent mass-congestion pipeline — seeds condense wherever local coordinate-entropy demand exceeds the boundary-CFT capacity ceiling via a tri-pass compute:cs_accumulate_cicfixed-point atomic Cloud-In-Cell mass assignment,step_incubation_transport6-neighbor Laplacian continuum incubation transport with KSS shear diffusion κdiff = 2/(17π) ≈ 0.0374475, Stinespring dilation sink, C² quintic partition of unity W(u) = 10u³ − 15u⁴ + 6u⁵ over z ∈ [18, 30],cs_detect_condensation26-neighbour non-maximum suppression with NMS collapse threshold ηV δc0 = 0.0198818 + 3³ basin overflow-mass integration, andcs_temporal_trackingminimum-image persistent-ID tracking with D₅ root-lattice merge radius rmerge★ = √2 κ★D₅ Δxcell ≈ 1.398249 Δxcell, accretion Ṁseed, and Landauer debt Pdebt,k = Mseed,k·906 GW/M⊙),causal_point_get.wgsl(active Causal-Point GET dynamics —cs_causal_point_getapplies the Rentropy = Nlimit − Cget ≥ 0 admissibility gate and the entropic acceleration aGET(x) = −κGET ∇ln ρproj(x) regularized by the Cardy capacity floor ρfloor(a) = max(0.00696738, 1/(39a)) with the first-principles transport coefficient κGET(fload) = (1/(d1 − h∨su(3)))(1 + (ceff/d1)fload) = (1/23)(1 + (1325/4004)fload) streamed via aTheoryInvariantsuniform),nbody_pm.wgsl(symplectic Martel–Shapiro supercomoving KDK integrator —cs_advance_particleshalf-kicks at (x̃n, an), drifts x̃ += Δτ·p̃ with periodic wrap, then half-kicks at (x̃n+1, an+1) using Δτ = Δa/(a½³·H(a½)); the conformal coupling A(a) = (3/2)Ωm,0a·(1 − (10/33)floadcosmo) multiplies the PM force directly, while seed attraction uses adaptive Plummer softening ε̃ = ηsoft/Ngrid (η ≈ 0.3333, with holographic softening length ε̃holo(a) = 2πa / (39 Ngrid)) and the GET term −κGET∇ln ρproj samples the register density carried in the PM texture alpha channel; the thermal Stinespring expectation wvis(z) = (1 − ηD) + ηD/(1 + (zN/z)26/3) with zN = 7.356 × 1010 drives the progressive quench; P-PM baryon pressurecompute_hydro_forceadds −cs²·a·∇ln ρb on the visible sector only — the hi ≥ 23/33 gauge — with gapless cs(z) derived from Compton thermal decoupling zdec = 137.036 ≃ 137 and acoustic drag zdrag = 1059.940),dual_channel_render.wgsl(Channel A emission scaled strictly by wvis(z) — solar-gold 10/33 baryons vs electric-violet 23/33 anti-baryons with the z ∈ [1012, 109] thermal cooling gradient — paraxial Gaussian radiance splats exp(−3.5·d²) reconciled with the topological Lorentzian profile Φtopo(σ) = (1 + (3513/140)σ²)⁻¹ via rank(D₅) = 5 (d² = 5σ²), plus Landauer heat glow (Pdebt = (Mseed/M☉) × 906 GW) rendered as an amber/white core + corona on each emergent seed; Channel B passive ghost shear/density conserved independently of it, carrying faint causal-point projection envelopes), andholographic_post.wgsl(physical distance-duality lensing — θE,k = √(4GM/c²·Dds/(DdDs)) per emergent seed with Tier-1 angular-diameter distances, boundary dispersion δdisp(λ) = ζ[(550 nm/λ)² − 1] at the 436/546/700 nm bands, half-resolution two-pass Gaussian bloom (bloom_blur.wgsl— bright-pass + separable 9-tap blur lifted into the composite viapost2.z), exponential depth fog along the Channel-A depth field (post2.w), and anti-blowout luminance-preserving ACES filmic tonemapping via the direct Narkowicz luminance tone reproduction curve with Landauer desaturation knee Ydesat = 14; the deflection marches a 4-slice multi-plane lens stack at zm ∈ {0.5, 1.2, 2.2, 3.5} with Dms/Ds weights inpost3plus a Shapiro-delay tint, and a split-viewport mode (post1.y) rasterizes the diagnostics phase-space texture into the right half of the frame;unwrap_torus_projectionconformal unwrapping onto the CFT torus [0, 2π)², horizon overlay driven by fload). Acs_sandbox_getpass incausal_point_get.wgslapplies an FDT friction counter-term −γobs·vproj·dir inside each dispatched observer's cone shell while depleting its entropy budget. A dedicateddiagnostics.wgslstage accumulates a 64-bin matter power spectrum by workgroup atomics (DIAG_FIXED_POINT= 4096) and rasterizes a 256² (x, vx) phase-space texture — gold visible sector, blue dark sector — entirely on-GPU with no host readback.- Seed condensation physics: condensation is an instanton-gated tunneling event. The Cardy capacity ceiling Nlimit(z) = γCFT·(H/Href)²·(Vcell/VH)·Nsat with γCFT = ceff/6 = 1325/924 ≈ 1.4339826 sets the per-cell bit budget; below the ceiling the Euclidean instanton action Sinst = (2πceff/kq)·((Nlimit − Nlocal)/Nlimit)² gates nucleation via Γnuc = A0·e−Sinst and Pnuc = 1 − e−Γnucδt, while overflow is barrierless (Sinst = 0). A PCG-hashed tunneling draw decorrelated per frame drives nucleation into a dynamically sized
GlobalSeedBufferpool (SEED_POOL_CAP= 1024). Condensed mass integrates the absolute register excess Σbasin max(0, Nlocal − Nlimit), which stays bounded as Nlimit → 0 in the z → −1 asymptote. Kernel work is evaluated on the normalized ratio R = Nlocal/Nlimit because absolute bit counts (~10122) overflow f32 registers. - Dimensional metrology:
src/shbt/visualizer/src/units.rsanchors every GPU quantity to physical units —CosmologicalContextfixes the scaling matrix Lbox = 222.5502 Mpc = 150.0 h−1 Mpc, V₀ = H₀Lbox, Mbox = ρcrit,0Ωm,0Lbox³, Gcode = (3/2)Ωm,0;MetrologyPipeline::prepare_step_uniformsevaluates the mid-step supercomoving Δτ, adaptive ε̃², and A(a) into the byte-exactGpuSimulationUniformsblock each frame; angular-diameter distances χ(z) integrate the Tier-1 HSHBT(z) background so the Einstein-radius telemetry reads out in arcseconds. - Zero-allocation buffer architecture:
ParticleBuffer_Ping/Pong,TetherVertexBuffer, andTetherIndirectArgsare allocated once at init and ping-ponged per frame — no per-framecreate_buffercalls on the hot path. Emergent-seed telemetry reaches the HUD through a 16 + 64×32 + 32³×4-byte staging ring (SeedReadback); on wasm targets where GPU→CPUmapAsyncis unstable, a host-side CPU replica of the CIC + NMS + tracking pipeline (cpu_emergence_tick) feeds the sameabsorb_seed_recordspath so telemetry stays live while the GPU still runs the true compute. - 2D conformal boundary-unwrap overlay (
BoundaryTex2D), Landauer-debt heat-map emission, entropy-gradient streamlines (entropy_tracer.wgsl/tracer_render.wgsl), Stinespring transition tethers (tether_render.wgsl— line-list indirect draw over atomically appended vertex pairs, alpha decaying exp(−t/0.5 s) percosmo.wall_dt), anti-baryon violet charge fringing (410 nm), tidal-shear billboard stretching (area-preserving strain along θshear = ½ atan2(2γ₂, γ₁) on Channel-B ghosts), past-light-cone volumes (causal_cone_render.wgsl— cyan wireframe spokes on the top-10 most-active mass-ordered observers, fading as the entropy budget depletes), entropy-budget spheres (causal_sphere_render.wgsl— instanced Fresnel shells scaled (Rentropy/Nlimit)1/3, emerald → crimson as Rentropy → 0), history-crystallization flashes, wave-optics caustic (Airy) fringing, and emergent-seed glitch UV quantization. - Seed Glitch Toggle: Enhancement 11's condensation glitch is a subtle, localized nuance rather than a frame-wide artifact — the glitch radius is tightened to 0.025 screen units, the UV quantization is lifted to 256 cells/axis, and the blend opacity is capped at 25% scaled by the
glitchslider (default 0.10). Theglitch-togglecheckbox in the optics row disables the effect entirely (GlitchUniforms.u_glitch_enabled = 0→apply_condensation_glitchearly-outs) for clean recordings, andset_glitch_enabled/set_glitch_intensityare exposed on the Wasm engine;glitchEnabled/glitchIntensityride thehud_jsontelemetry frame. WasmShbtEngineowns the zero-copy simulation buffers —ParticleRecord(32 B),CausalPointRecord(64 B),SeedDefectRecord(128 B) — exposed via rawget_*_buffer_ptraccessors, withupdate_epoch(z)driven byHorizonLedger(Nsat = 3.3119977 × 10122 bits, Γlock, Stinespring partition ηA = 10/33 / ηD = 23/33, seed inventory, Landauer debt, observer admissibility ℛadm, conservation residual).- Double-buffered storage holds 220 particles in 32-byte
Particlerecords (< 256 MB); the HUD decodes the 128-byte SHBT-MMIO telemetry frame (Horizon Bar, Boundary Capacity Gauge, Congestion & Seed Ledger, Landauer Debt Monitor, Δfr = 0 / Eμν = 0 / horizon-freeze indicators). - Continuum incubation transport & multi-fluid thermodynamics:
step_incubation_transportinseed_emergence.wgslmatures the precursor congestion field across z ∈ [10⁹, 30] — a 6-neighbor Laplacian with KSS shear diffusion κdiff = 2/(17π) ≈ 0.0374475 redistributes nlocal isotropically while the Stinespring sink Sdil = ηD·ΓS·ρtot (ηD = 23/33, ΓS = |dwvis/dz|) accumulates in a persistent dark ledger, so condensation ignites out of transported structure rather than raw noise. Ignition is gated by the C² quintic partition of unity Ψnuc(z) = u³(10 − 15u + 6u²), u = (30 − z)/12 over z ∈ [18, 30], eliminating pop-in impulses. The multi-fluid solver deposits both sectors onto atomic grids at 2²⁰ fixed point (deposit_mass_cic→rho_tot_atomic/rho_baryon_atomic), thencompute_thermodynamic_pressurebuilds Pb = cs²·Tb(z)·ρb with the exact decoupling profile Tb(z) = TCMB,0(1 + z)·z/(z + 137); the gather force reuses the same CIC kernel so self-forces cancel telescopically (verified < 10⁻¹⁵ in 256-bit arithmetic). Each particle'spad.xcarries its normalized local incubation ratio, which smoothsteps the visible-sector hydro gate at hi ≥ 23/33 and drives the cyan precursor glow (0.10, 0.95, 0.95) indual_channel_render.wgsl. - Multi-tier density shading, branchless blackbody & anti-blowout ACES:
compute_thermodynamic_statemaps the congestion pair (ρb, σ) — normalized incubation density and boundary-capacity ratio against the unscaled Cardy ceiling — onto a smooth effective temperature Teff = 2500 + 6500·tanh(max(0, ρb−1)/4) + 16000·σ³/(1+e−12(σ−1)) ∈ [2500, 25000] K and a visible spectral radiance Ivis with a bounded log-mass seed corona.blackbody_to_linear_rgbevaluates the u = 1000/T ∈ [0.040, 0.400] Planckian palette as three degree-(2,2) rational Padé channels with zero divergent control flow, and each splat's smoothing length follows the SPH spacing law R ∝ (1+0.45·max(0, ρb−1))−1/3 contracted by the topological Lorentzian factor Φtopo(σ) = (1 + (3513/140)σ²)⁻¹ (reconciled with the Gaussian splat exp(−3.5·d²) via rank(D₅) = 5, d² = 5σ²) and projected through camera distance (16 px cap, anchored by r₀ = √2 κ★D₅ − 1/(39√2) = 1.380121 ≃ 1.38), so dense cores tighten while underdense filaments diffuse. Inholographic_post.wgsl, a 9-tap cross-bilateral (spatial σs = 1.8, photometric σr = 0.10) feeds the luminance Laplacian that accents caustic halos by the telemetry-driven Jacobian detJ = |(1−κmax)²−γmax²|, andtone_map_aces_anti_blowoutapplies the direct Narkowicz ACES luminance tone reproduction curve — re-injecting un-clipped chromaticity and rolling saturated cores into a controlled cyan-white Landauer corona (knee Ydesat = 14, capped at 0.85) — followed by the accurate sRGB EOTF. Shaders benchmark ≈ 8.89 ms/frame at 1080p (≈ 112.5 FPS) on real GPUs; the post uniform block is 272 B (post4= γmax, κmax, knee, tone exposure). - Monotonic cosmic-age clock & deep-space composite: the Horizon Bar reads Cosmic Age: t(z) Gyr (Lookback: tlb Gyr) from the canonical quadrature
cosmology::cosmic_age_gyr— t(z) = H0−1 ∫ln(1+z)40 du/E(u), strictly monotonic in −z (t(0) ≈ 13.79 Gyr, finite de Sitter future for −1 < z < 0, +∞ at z ≤ −1). The phase banner is driven by live telemetry (quench %, seed count, Mseed, fload, Radm) so Phase 3 incubation vs. condensation and Phase 5 freeze stay strictly synchronized with the ledger. The seed pool scales toSEED_POOL_CAP= 1024 acrossseed_emergence.wgsl,nbody_pm.wgsl,dual_channel_render.wgsl,holographic_post.wgsland the hostMAX_SEEDS, andnbody_pm's seed-force loop integrates the full pool. On the optics side, zero-flux pixels map to vec3(0) end-to-end: the depth-fog ambient is true black, the cross-bilateral tightens to σr = 0.025, splats are gated off below ρb < 0.02, and a 16-tap wide-scale pedestal subtracts the diffuse mass sheet (optically invisible by mass-sheet degeneracy) so only convergence/luminance contrast above the local mean emits — void floors collapse to deep black while filaments and seed cores keep full radiance. AShow Phase-Space Diagnosticstoggle bounds the (x, vx) raster inside a styled inset card ("Phase-Space Density (x, vx)" with axis labels), and the seed glitch defaults to OFF for canonical clean recordings. - First-principles cosmic-web shader deployment: every pixel attribute is theory-derived — ad-hoc visual multipliers are eliminated end-to-end.
blackbody_to_linear_rgbis the canonical branchlessevaluate_blackbody_simd: three degree-(2,2) rational Padé channels in u = 1000/Teff ∈ [0.040, 0.400] map the thermodynamic band Teff ∈ [2500, 25000] K to linear RGB with strictly positive denominators ([0, 1.05]-bounded bytest_branchless_blackbody_simd_rational_monotonicity). The dynamic splat radius contracts over-capacity cores through the topological SPH factor Φtopo(σ), andholographic_post.wgslapplies the ε-softened thin-screen magnification 𝒜caustic = 1/√(det J² + ε²), det J = (1−κ)²−γ² (ε = 0.08), with the κ-central-difference macro-deflection, per-seed Einstein-ring deflections, and 3-tap dispersion at (436, 546, 700) nm. KDK symplectic energy drift is gated at < 10⁻⁴ over 10⁴ steps bytest_symplectic_energy_drift_static_background, andtests/test_harness_metrology.pyadds the 512-bit metrology gates (bit conservation < 10⁻³⁵, stress residual < 10⁻¹³⁰, Landauer debt scaling, wasm heap < 256 MB). - Interactive controls: timeline scrub z = 1014 → −1 (ghost-seed condensation highlighted across z ≈ 30 → 7 with ΔN ≈ 6 × 1059 bits and Pdebt ≈ 9.06 × 1020 W), playback speeds 1×/10×/100×, continuous
unwrap_transitiontorus slider plus comoving-bulk / boundary-CFT projection switch, and Channel A/B toggles. An epoch quick-jump bar jumps directly to each cosmic milestone (bit loading z = 1014, baryogenesis z = 1011, seed genesis z = 18, proto-galaxy web z = 3, horizon freeze z = −0.999) with matching playback speed, and a color-coded phase banner identifies the active era (cyan bit-loading, flashing magenta baryogenesis, amber seed condensation, blue/white GET clustering, emerald de Sitter freeze). The HUD adds a Causal-Point Observer Activity monitor (ℛadm cardinality and Rentropy ≥ 0 / freeze status) and a Gravitational Optics Telemetry ledger (γmax, κmax, θE, active caustics). Timeline drags are eased in ln(1 + z) space and the camera orbit direction is slerped through the torus-unwrap transition, so epoch scrubs land smoothly instead of snapping. - Decoupled telemetry fractions & dynamic seed lifecycle: the Γ-lock kernel K(u) = Γlocke−u/HSHBT(u) defines two distinct physical loading channels — the backward past-light-cone debt floaddebt(z) = ∫0ln(1+z) K du (≡ 0 for z ≤ 0, saturates at 0.10744 at zrec) and the forward boundary-capacity fraction floadcosmo(z) = 1 − exp(−∫ln(1+z)40 K du) (≈ 0 primordial, ≈ 0.10 today, locks to 1 as z → −1 — the de Sitter horizon freeze). Both are serialized into the SHBT-MMIO v2 128-byte frame (fcosmo @ 0x28, fdebt @ 0x30, CRC32 over 0x00–0x77) decoded by the HUD with 1:1 seed↔caustic parity;
nbody_pm.wgsl's Poisson coupling A(a) = (3/2)Ωm,0a(1 − (10/33)floadcosmo) reads the forward fraction so high-z growth restores the Einstein–de Sitter kernel. Inseed_emergence.wgsl, Ψnuc keeps an open floor ≡ 1.0 for z ≤ 18 (C² quintic over the incubation band z ∈ [18, 30]), the NMS margin uses the capacity floor prefactor ηV δc0 = 0.0198818 via δth(z) = ηV δc0 (H/H₀)1/3(1 − (10/33)fcosmo)/Ψnuc, and a three-stage lifecycle —cs_temporal_trackingclaimsprev_claimedslots →cs_seed_carrycarries un-re-detected defects →cs_seed_mergeperforms a serial agglomerative merge at rmerge★ = √2 κ★D₅ Δxcell ≈ 1.398249 Δxcell conserving mass/debt — pools candidates and carried seeds up to the 1024-seed cap. Here zpeak = 24.0 and Δz = 6.0 (σz = 3.0) are derived directly from the branch invariants ((26, 8, 312), ceff = 1325/154, γCFT = 1325/924) inderive_first_principles_z_peak/derive_first_principles_sigma_z. - Cinematic recording suite & 8-milestone visual verification (
visualizer/record_simulation_events.py): composes a continuous 90 s / 5400-frame flythrough at 60 FPS through a Catmull-Rom camera spline eased in log-redshift χ(z) = log₁₀(1 + z) — epoch budgets z = 1014→1010:15 s, 100→18:20 s, 18→1:25 s, 1→0:15 s, 0→−0.999:12.5 s, with a 30-presented-frame physical-relaxation dwell parked at each milestone (SHBT_MILESTONE_DWELL) —update_cosmic_state(z, eye, look_at)pins the engine's camera override each commit — while Playwright recordsfull_cosmic_evolution.webmplus the 8-milestone visual verification sequence:01_primordial_bit_loading.png(z = 1014: age < 0.001 Gyr, fload < 0.01),02_baryogenesis_derendering.png(z = 1010: quench window 68–71%, ηD = 23/33 ≈ 69.7%),03_ghost_seed_genesis.png(z = 16.0: seeds > 0, Mseed > 0, Pdebt > 0),04_early_filament_network.png(z ≈ 3.86: isotropic early filament web and initial Zel'dovich/2LPT sheet intersections),05_causal_point_proto_galaxies.png(z = 3.00: Landauer debt > 5 × 10¹¹ GW, quench lock 69.7% ± 1%, seeds = 1024, caustic:seed 1:1 parity),06_cosmic_web_lensing.png(z = 1.00: caustics = 1024, caustic:seed 1:1 parity, γmax > 0.1),07_present_epoch_universe.png(z ≈ 0.822 → 0: late-time virialized cluster nodes, 4-slice multi-plane lensing, and void clearing), and08_desitter_horizon_freeze.png(z = −0.990 → −0.999: age > 15 Gyr, fload > 0.99, Rentropy < 0, Radm = 0). Asserts read the dedicated#hud-telemetry-overlayDOM (#hud-redshift-val,#hud-quench-fraction-val,#hud-landauer-debt-val,#hud-fps-val; the FPS gate is advisory below 55 unlessSHBT_STRICT_FPS=1, andSHBT_FRAME_STRIDEdecimates commits on software rasterizers). Runpython3 visualizer/record_simulation_events.py(serves the repo root itself; needs Playwright Chromium with--enable-unsafe-webgpu --use-angle=vulkan --enable-features=Vulkan, uses the?capture=1overlay path on headless machines) andpython3 visualizer/record_simulation_events.py --verifyto re-check artifact existence, >100 KB size, and non-blank pixel variance. The canonical suite records clean physics (seed glitch toggled off);visualizer/record_glitch_suite.pyproduces the three dedicated recordingsemergent_seed_formation_clean.webm(glitch OFF, z = 30 → 7),emergent_seed_glitch_nuance.webm(glitch ON @ 0.10, z = 14 → 7), andshbt_full_thermodynamic_cycle.webm(z = 1014 → −0.999).
The sys1own/shbt-precision visualizer features a high-performance WebGPU gravitational optics engine that renders cosmological mass deflections, chromatic caustics, and relativistic Doppler shifts at 60 FPS for 1M+ particles.
dual_channel_render.wgslapplies relativistic Doppler beaming 𝒟 = √(1 − β²)/(1 − βlos) with boosted surface brightness Iobs = I0𝒟³ and akelvin_to_rgbblack-body thermal shift; Channel B carries the projected shear (γ1, γ2), convergence κ, and causal entropy.holographic_post.wgsl(fs_post) resolves the screen-space lens equation β = θ − λlens[∇κ + Γ·∇κ] − Σs θE,s²(θ − θs)/(|θ − θs|² + εcore²) via central differences of the Channel B convergence field plus an analytical softened point-mass seed loop over the emergent pool (S ≤ 1024, dominant θE ≲ 0.09 rad measured), accents critical curves by the ε-softened caustic magnification 𝒜 = 1/√(det J² + ε²), applies three-tap chromatic dispersion δdisp at the (436, 546, 700) nm bands, a depth-aware bilateral dark-matter caustic glow (sample_bilateral_convergence), Einstein/caustic rings, and the Causal-Point Fresnel ripple field.causal_point_get.wgslrasterizes instanced Fresnel ripple shells for every observer node with Rentropy ≥ 0, writing signed shear envelopes and boundary convergence into Channel B.LensingUniforms(224-byte, 16-aligned#[repr(C)]) andSeedDefect(16-byte) storage tables are staged each frame;VisualizerTelemetryreports γmax, κmax, θE, and the active caustic count to the HUD.
# Launch the headless visualizer with gravitational lensing and chromatic dispersion
cargo run --release --manifest-path src/shbt/visualizer/Cargo.toml \
--bin headless -- \
--enable-lensing --lensing-scale 1.5 --dispersion 0.25 \
--enable-doppler --dark-glow 0.8 --width 2560 --height 1440| Key | Action |
|---|---|
L |
Toggle gravitational lensing pass on/off |
D |
Toggle relativistic Doppler beaming and thermal color shifts |
[ / ] |
Decrease/increase global lensing deflection strength λlens (0–5) |
- / = |
Decrease/increase wave-optics chromatic dispersion δdisp (0–1) |
G |
Toggle dark-matter halo volumetric bilateral glow |
H |
Toggle telemetry HUD overlay visibility |
Equivalent HUD controls: the Lensing, Doppler, and Dark glow toggles plus the λlens, δdisp, and glow sliders.
- Zero GR compute overhead: macro-scale lensing operates on screen-space potential gradients in Channel B while micro-scale lensing uses analytical softened point-mass seeds over the emergent pool (S ≤ 1024); no geodesic integration runs in the browser loop.
- Deterministic memory model: WebAssembly linear-heap usage stays under 256 MB with zero heap allocations inside the active rendering loop.
- Invariant preserving: visual shaders run fully decoupled from the cosmological physics solvers (
precision_cosmology.py,boltzmann_shbt.py), preserving analytical invariance (Eμν = 0, Δnorm < 10−120).
Reference GPU benchmark at N = 2²⁰ = 1,048,576 particles, 1920×1080, fixed initialization buffers (zero render-loop allocations):
| Pipeline stage | NVIDIA RTX 4090 (Vulkan) | Apple M2 Max (Metal) | Threshold | Status |
|---|---|---|---|---|
nbody_pm.wgsl KDK leapfrog |
2.84 ms | 5.12 ms | < 8.00 ms | PASS |
seed_emergence.wgsl CIC + instanton |
1.15 ms | 2.41 ms | < 4.00 ms | PASS |
causal_point_get.wgsl GET clustering |
1.42 ms | 2.85 ms | < 4.00 ms | PASS |
dual_channel_render.wgsl SPH splats |
4.21 ms | 5.38 ms | < 7.00 ms | PASS |
holographic_post.wgsl lensing + ACES |
0.68 ms | 0.94 ms | < 1.50 ms | PASS |
| Aggregate frame latency | 10.30 ms | 16.70 ms | < 16.66 ms | PASS |
| Sustained viewport rate | 97.1 FPS | 60.2 FPS | ≥ 60 FPS | PASS |
| Bit partition ΔN/Nsat | 4.12 × 10−32 | 4.12 × 10−32 | < 10−30 | PASS |
| Stress-energy |Eμν| | 1.08 × 10−122 | 1.08 × 10−122 | < 10−120 | PASS |
| KDK drift (10⁴ steps) | 3.42 × 10−5 | 3.42 × 10−5 | < 10−4 | PASS |
| Wasm linear memory | 128 MB | 128 MB | < 256 MB | PASS |
| Render-loop allocations | 0 B/frame | 0 B/frame | 0 B/frame | PASS |
# Foundation audit (modular closure, projection, memory, baryogenesis)
python shbt_simulate.py --mode audit
# Full unified simulation (foundation + precision cosmology)
python shbt_simulate.py --mode all --output result.json
# Full pipeline with closed-loop observer succession (5-phase lifecycle,
# Stinespring de-rendering, transfer kernel over the 3x3 coordinate lattice)
python shbt_simulate.py --mode all --enable-succession --succession-cycles 3 --output result.json
# Boltzmann CMB and matter spectra generation
python boltzmann_shbt.py # runs the full pipeline, prints CMB spectra
python boltzmann_shbt.py --run-tests # embedded Boltzmann/chronometer unit tests
# Precision cosmology report and parameter sweeps
python precision_cosmology.py --json
python precision_cosmology.py --run-tests
python shbt_simulate.py --mode cosmology --sweep sweep.json --plot --output sweep_result.json
# WebGPU visualizer: telemetry frame + wasm bundle + serve
python shbt_simulate.py --mode visualize --particles 262144 \
--export-webgpu-telemetry data/telemetry.bin --sim-speed 1.0 --output viz.json
cargo build --release --manifest-path src/shbt/visualizer/Cargo.toml \
--target wasm32-unknown-unknown
wasm-bindgen --target web --out-dir visualizer/pkg \
src/shbt/visualizer/target/wasm32-unknown-unknown/release/shbt_visualizer.wasm
python -m http.server 8080 --directory visualizer # open http://localhost:8080
# Cinematic 3D isometric recording and automated milestone verification
python visualizer/record_simulation_events.py # records full_cosmic_evolution.webm & 01-06 PNGs
python visualizer/record_simulation_events.py --verify # verifies artifact size, non-blankness, and DOM asserts
# Headless GPU benchmark (no browser): frames + particles
cargo run --release --manifest-path src/shbt/visualizer/Cargo.toml \
--bin headless -- 60 262144
# Sub-10 mK Landauer calorimetry audit (writes shbt_run_calorimetry_sim.csv)
python shbt_simulate.py --mode all --output result.json
# or programmatically:
python -c "import precision_cosmology as pc; pc.simulate_calorimetry_experiment(n_pulses=1000000)"Additional modes: baryogenesis (topological asymmetry benchmark), visualize (SHBT-MMIO telemetry export + cosmology_visualization and causal_point_simulation records; flags --particles, --export-webgpu-telemetry, --sim-speed), history (Causal-Point observer crystallization; with --enable-succession runs the multi-cycle lifecycle engine and exports succession records to result.json), cosmology-test. Exports support --format json|csv|hdf5, --plot, --sweep, --config (YAML/JSON), --seed, and structured --log-format json logging.
The 3D isometric orbital trajectory matrix traverses the five cosmic evolutionary epochs (θ ≈ 35.264°, Δϕ = 45°), preventing Cartesian line-of-sight caustic stacking (bullseyes) and rendering true filament structures. Zero-parameter first-principles mode locks all sandbox controls (cs, ηperc, γobs) to unit gain.
| Publication claim | Implementation | Output artifact |
|---|---|---|
| Canonical branch uniqueness & Weil conductor rigidity | src/shbt/uniqueness.rs (run_uniqueness_audit, DiophantineClassifier, WeilModuleVerifier) |
result.json: foundation_audit.uniqueness_proof |
| Bulk gauge dynamics from boundary Kac-Moody currents | src/shbt/gauge_dynamics.rs (run_gauge_closure_audit, BulkGaugeSlice, derive_gauge_connection) |
result.json: foundation_audit.bulk_gauge_dynamics |
| SM running gauge couplings & 512-bit Stinespring RG flow | src/shbt/couplings.rs (run_gauge_couplings_audit, GaugeCouplingLedger) |
result.json: foundation_audit.standard_model_couplings |
| Fermion masses, CKM & PMNS flavor mixing via Verlinde SVD | src/shbt/flavor.rs (run_flavor_audit, FlavorMixingEngine, Matrix3x3) |
result.json: foundation_audit.flavor_mixing |
| Electroweak symmetry breaking, Higgs mass & quartic coupling | src/shbt/higgs.rs (run_higgs_audit, ElectroweakSymmetryBreakingEngine) |
result.json: foundation_audit.electroweak_symmetry_breaking |
| Non-perturbative QCD, chiral condensate & Strong CP lock | src/shbt/qcd.rs (run_qcd_audit, NonPerturbativeQcdEngine) |
result.json: foundation_audit.non_perturbative_qcd |
| Non-equilibrium primordial nucleosynthesis (BBN) | src/shbt/thermal_history.rs (run_thermal_history_audit, NonEquilibriumBbnEngine) |
result.json: foundation_audit.primordial_nucleosynthesis |
| Non-linear structure formation & tomographic cosmic shear | src/shbt/structure.rs (run_structure_audit, StructureFormationEngine) |
result.json: foundation_audit.cosmic_shear_and_structure |
| Cryogenic 7 mK Landauer calorimetry & Neyman-Pearson audit | src/shbt/calorimetry.rs (run_calorimetry_audit, CryogenicCalorimetryEngine) |
result.json: foundation_audit.calorimetry_hypothesis_audit |
| Bispectrum & trispectrum templates | boltzmann_shbt.py |
result.json: precision_pipeline.non_gaussianity |
| GET capacity & cutoff | src/shbt/causal_point.rs (CausalPoint.crystallize_history) |
result.json: foundation_audit.memory_report |
| Observer admissibility & de-rendering | src/shbt/causal_point.rs (CausalPoint.is_admissible, terminate_and_derender) |
result.json: succession.records[].eta_dark, .pointer_triad |
| Succession transfer kernel | src/shbt/causal_point.rs (evaluate_succession_kernel, ShbtSimulator.run_succession_cycles) |
result.json: succession.records[].kernel_probabilities, .kernel_normalized |
| Asymptotic observer freeze | precision_cosmology.py (asymptotic_observer_freeze) |
result.json: foundation_audit.asymptotic_observer_freeze |
| Landauer calorimetry regression | shbt_simulate.py / precision_cosmology.py |
shbt_run_calorimetry_sim.csv |
| Dynamic bit-loading Ṫ = H(t) (no inflaton) | src/shbt/cosmology.rs (ShbtUniverse.evaluate_hubble_index_rate) |
result.json: cosmology_visualization; tests/cosmological_invariants.rs |
| Passive ghost dark matter (ηD = 23/33, ΩDM,0 = 0.260000) | src/shbt/cosmology.rs (dark_completion_bits, dark_matter_fraction) |
data/telemetry.bin header fields |
| Mass-congestion condensation (Mseed = αseed ΔN) | src/shbt/cosmology.rs (compute_seed_condensation, eval_landauer_debt_power) |
SHBT-MMIO SeedDefectRecord payload |
| SHBT-MMIO telemetry frame | src/shbt/export.rs (serialize_mmio_frame_py, MmioTelemetryHeader.to_bytes) |
data/telemetry.bin |
| Two-tier WebGPU visualizer | src/shbt/visualizer (ShbtWebGpuEngine, WGSL pipeline) |
visualizer/pkg/shbt_visualizer_bg.wasm |
| Ergodic CFT modular phase seeding | src/shbt/structure.rs (generate_weil_modular_phases), boltzmann_shbt.py (generate_weil_modular_phases_grid) |
tests/test_shbt.rs |
| Linear box metrology & 2LPT initialization | shbt_simulate.py (compute_2lpt_displacements), boltzmann_shbt.py |
data/telemetry.bin |
| Physical observer frame & sub-linear orbit | src/shbt/visualizer/src/units.rs (PhysicalUnitTransformer), record_simulation_events.py |
visualizer/recordings/ milestone verification |
| Regularized entropic transport & quartic Cardy switch | src/shbt/visualizer/src/shaders/nbody_pm.wgsl, causal_point_get.wgsl |
tests/test_harness_metrology.py |
| Symmetrized Green's function & phase-space atomic binning | src/shbt/visualizer/src/shaders/nbody_pm.wgsl, diagnostics.wgsl |
visualizer/recordings/ milestone verification |
| Artifact | Columns | Content |
|---|---|---|
shbt_run_cmb_cls.csv |
ell, Dl_TT_muK2, Dl_EE_muK2, Dl_TE_muK2, Cl_TT, Cl_EE, Cl_TE, Cl_BB |
Scalar temperature & polarization spectra |
shbt_run_matter_pk.csv |
k_Mpc_inv, Pk_z0, Pk_z05, Pk_z1 |
Matter power spectrum at z = 0, 0.5, 1 |
shbt_run_tensor_cls.csv |
ell, Cl_BB, Dl_BB, Cl_TT_tensor |
Primordial tensor B-modes |
shbt_run_calorimetry_sim.csv |
k_bits, R_addresses, Q_H0_zJ, Q_H1_zJ, Q_noise_zJ |
Sub-10 mK calorimetry simulation |
result.json |
nested ledger | Master simulation result tree |
data/telemetry.bin |
128-byte SHBT-MMIO header + T/Pm grids + seed records | Tier-2 visualizer telemetry contract |
visualizer/pkg/ |
shbt_visualizer.js, shbt_visualizer_bg.wasm |
Browser engine bundle |
| Symbol | Kind |
|---|---|
shbt_simulator |
Master PyO3 module |
StaticBoundary |
Boundary CFT modular data, framing defect, dark ledger |
ShbtSimulator |
Orchestrating runtime (run_full_audit, run_uniqueness_audit, run_gauge_closure_audit, run_gauge_couplings_audit, run_flavor_audit, run_higgs_audit, run_qcd_audit, run_thermal_history_audit, run_structure_audit, run_calorimetry_audit) |
DiophantineClassifier / WeilModuleVerifier |
Branch classification & discriminant rigidity (src/shbt/uniqueness.rs) |
BulkGaugeSlice / GaugeGroup / KacMoodyBoundary |
Non-Abelian gauge connection & bulk YM equations (src/shbt/gauge_dynamics.rs) |
GaugeCouplingLedger / RGSliceRecord |
512-bit Stinespring Callan-Symanzik beta flow (src/shbt/couplings.rs) |
FlavorMixingEngine / Matrix3x3 / Complex64 |
Verlinde fusion & Jacobi bi-unitary SVD (src/shbt/flavor.rs) |
ElectroweakSymmetryBreakingEngine |
512-bit EWSB & Higgs potential engine (src/shbt/higgs.rs) |
NonPerturbativeQcdEngine |
512-bit 4-loop QCD & Strong CP resolution (src/shbt/qcd.rs) |
NonEquilibriumBbnEngine |
8-species stiff BBN nuclear network (src/shbt/thermal_history.rs) |
StructureFormationEngine |
Cardy capacity HMF & Limber tomographic cosmic shear (src/shbt/structure.rs) |
CryogenicCalorimetryEngine |
Cryogenic 7 mK thermal SDE & Neyman-Pearson detector (src/shbt/calorimetry.rs) |
HolographicProjection |
RG flow → BulkMetricSlice |
CausalPoint |
Observer memory & history crystallization |
CausalPoint (succession) |
is_admissible, terminate_and_derender, evaluate_succession_kernel, relabel_and_rerender, run_lifecycle_cycle |
CausalPointCandidate / DerenderingRecord / SuccessionRecord / LifecyclePhase |
Observer-succession record types (Stinespring de-rendering, transfer kernel, 5-phase lifecycle) |
AnomalyClosureError |
Algebraic anomaly failure type |
ShbtUniverse |
Branch cosmology: loading, Hubble, ghost partition, seed condensation, horizon freeze (src/shbt/cosmology.rs) |
MmioTelemetryHeader / SeedDefectRecord |
SHBT-MMIO serialization (src/shbt/export.rs) |
ShbtWebGpuEngine / Particle / HudMetrics / TimelineController |
Tier-2 visualizer (src/shbt/visualizer) |
See paper_references.md for the complete paper-section → method crosswalk.
# Prerequisites: Rust 1.80+ (rustup), Python 3.8+, pip install maturin
pip install -r requirements.txt
# For the WebGPU visualizer additionally:
rustup target add wasm32-unknown-unknown
cargo install wasm-bindgen-cli # matches wasm-bindgen crate version
# Build Rust library and PyO3 bindings
cargo build --release
maturin develop --release # installs shbt_simulator into the active env
maturin build --release # produces target/wheels/shbt_simulator-*.whl
# Unit and integration tests
cargo test --release # Rust suite (tests/test_shbt.rs)
pytest tests/ # Python suite (tests/test_simulator.py)
# Compile publication documents
make pdf # pdflatex main.tex (primary manuscript)
latexmk -pdf -jobname=main main.tex
latexmk -pdf -jobname=supplementary supplementary.tex
# Whole pipeline: build + test + audit + pdf
make allThe SHBT program is a federated ecosystem of 10 configuration-controlled repositories. shbt-precision is the mathematical core; downstream repositories consume its audited invariants.
+------------------------------------------+
| [shbt-precision] |
| Computational Math & Cosmology |
| (512-bit MPFR / WZW Characters) |
+--------------------+---------------------+
|
+---------------------------------+---------------------------------+
v v v
+---------------------------+ +---------------------------+ +---------------------------+
| [shbt-power] | | [shbt-cf] | | [shbt-qc] |
| Commercial Fusion Grid | | 1,800-Module LANR Array | | Bare-Metal Microkernel & |
| (8,750 MW p-11B Twin) | | & Thermal-Hydraulics | | Photonic Quantum Bus |
+-------------+-------------+ +-------------+-------------+ +-------------+-------------+
| | |
+---------------------------------+---------------------------------+
|
+---------------------------------+---------------------------------+
v v v
+---------------------------+ +---------------------------+ +---------------------------+
| [shbt-ghost] | | [shbt-sglt] | | [shbt-exotic] |
| Ghost Seed Propulsion & | | Synthetic Gravitational | | Multi-Protocol Spacetime |
| Metric Interlocks (CCZ4) | | Lensing Telescope Stack | | Co-Simulation & QI Audit |
+-------------+-------------+ +-------------+-------------+ +-------------+-------------+
| | |
+---------------------------------+---------------------------------+
|
+-----------------------+-----------------------+
v v
+-------------------------------+ +-------------------------------+
| [shbt-recon] | | [shbt-warp] |
| Macro State Translocation & |<=============>| Holographic Warp Drive & |
| Bipartite Fuel Foundry Twin | Dual-Gateway | 3+1D ADM Metric Engine |
| (V_unified^macro / MAR-312C) | Corridor | (Alcubierre / SHBT-MAR-312C) |
+---------------+---------------+ +---------------+---------------+
^ ^
| +-------------------------------+ |
| | [shbt-swarm] | |
+=======| Autonomous Von Neumann Swarm |=======+
| Universal Matter Foundry & |
| 5th-Gen SHBT-MAR-312X Twin |
+-------------------------------+
│
┌─────────────────────────────────┼─────────────────────────────────┐
▼ ▼ ▼
╭───────────────────────────╮ ╭───────────────────────────╮ ╭───────────────────────────╮
│ [shbt-power] │ │ [shbt-cf] │ │ [shbt-qc] │
│ Commercial Fusion Grid │ │ 1,800-Module LANR Array │ │ Bare-Metal Microkernel & │
│ (8,750 MW p-11B Twin) │ │ & Thermal-Hydraulics │ │ Photonic Quantum Bus │
╰─────────────┬─────────────╯ ╰─────────────┬─────────────╯ ╰─────────────┬─────────────╯
│ │ │
└────────────────────────┬────────┴─────────────────────────────────┘
▼
╭──────────────────────────────────────────────────────────────────────────────────────────────╮
│ SPECIALIZED VEHICLE TWINS │
│ │
│ • shbt-ghost : Reactionless Propulsion & Local Gravity Wells (3+1 CCZ4 / PCSS Crowbars) │
│ • shbt-recon : Macroscopic State Translocation Gateway (Stinespring V_macro / 504 Gbps) │
│ • shbt-sglt : Synthetic Gravitational Lensing Telescope (SE-L2 Swarm / TMSV Metrology) │
│ • shbt-warp : Holographic Warp Metric & 3+1D Flight Twin (ADM α=1.0 / SHBT-MAR-312X Engine)│
╰──────────────────────────────────────────┬───────────────────────────────────────────────────╯
│
▼
╭───────────────────────────────────────────────────────────────────────────────────────────╮
│ shbt-exotic │
│ MULTI-PROTOCOL SPACETIME ENGINEERING CO-SIMULATION BENCH │
│ │
│ • Cross-Protocol Field Coupling (Warp + Stasis + Translocation + Wells + Comms) │
│ • Global Energy Condition & Ford-Roman Quantum Inequality (QI) Dark-Ledger Auditing │
│ • Dynamic 5-Stage Multi-Technology Flight Director & Relativistic PDE Mesh Solvers │
╰───────────────────────────────────────────────────────────────────────────────────────────╯
| Repository | Domain Role & Platform Scope | Shared Invariants & Interface Contracts |
|---|---|---|
shbt-precision |
Computational Math & Cosmological Foundation Core | 512-bit MPFR arbitrary-precision arithmetic (rug), WZW affine levels (26, 8, 312), 346 Dynkin slots, Callan–Symanzik RG thresholds, exact Casimir denominators (28, 22, 640), framing defect Δ_fr ≡ 0, Landauer debt schedule P_debt = 906.00 kW, 15-gate audit / 15 Z3 SMT proofs. |
shbt-qc |
Photonic Quantum Computer Twin & C11 Microkernel | Freestanding C11 shbt-os microkernel, SHBT-R photonic processor architecture, 52-hexamer InP/InGaAsP microcavities (3-stage polariton squeezing cascade), non-Hermitian PT-symmetric Geometry of Interaction (GoI) optical matrix engine (t_prop = 1.701 ps), Tomita–Takesaki inverse modular automorphism flow σ_{-t}^Ψ (Q_internal = 0 J/bit), 128-byte SHBT-MMIO-2 at 0x70000000, SECDED Hamming(72,64) ECC syndromes, 50-gate verification matrix (V-01..V-50), 10 Z3 SMT proofs. |
shbt-power |
Commercial p-¹¹B Aneutronic Fusion Power Plant Twin | Multi-gigawatt p-¹¹B aneutronic fusion power generation (8,750 MW fusion / 7,832.903 MW net export), 128-byte SHBT-MMIO-POWER, 906.00 kW Landauer debt schedule, closed-loop thermal ledger, 78-gate verification audit (78/78 PASS) + 84-EXT checks, 16 Z3 proofs. |
shbt-cf |
LANR Cold Fusion Reactor Workbench & Thermal-Hydraulics | 1,800-module LANR starter grid (555.03 W net/cell, 999.054 kW net DC array @ 400 V), dynamic Floquet–Adler–Wiser dielectric screening (U_eff = 352.48 eV), dual-stage CoSb₃/ZrNiSn TEG enthalpy recovery, Kapitza resistance ΔT_K = 3.546 K, 3D Eulerian-Eulerian helium thermal-hydraulics, 50-gate verification matrix (GATE-01–GATE-50). |
shbt-ghost |
Ghost Seed Reactionless Propulsion & Metric Stabilization | 3+1D CCZ4 numerical relativity with Gundlach damping (κ₁ = 0.15, κ₂ = 0.0), ADM metric stabilization (`β^i → 0, |
shbt-recon |
Macroscopic State Translocation & Bipartite Fuel Foundry Twin | Macroscopic Stinespring dilation (V_unified^macro), 101-layer Adaptive TTN-MERA, topological dark ledger (η_D = 23/33, 124 Fibonacci braids), bipartite coset fuel synthesis operator (R_synth^MAR, Δh_vis = 45/64), 64-bank UnifiedStinespringHyperFrame (135,168 B @ 0x70020000), Mode II 4.2 K readout (Q_readout = 4.0194 × 10⁻²³ J/bit), compacted SHBT-MAR-312C nacelle (2,257.70 kg), SHBT-IBR-312 magazine, 144-gate verification audit, 29 Z3 formal SMT proofs. |
shbt-sglt |
Synthetic Gravitational Lensing Telescope Stack | 2PN relativistic beam optics in harmonic coordinates, Two-Mode Squeezed Vacuum (TMSV) heterodyne metrology (r = 2.50, 21.715 dB, noise ASD S_r^{1/2} ≤ 0.010 pm/√Hz), 5th-order minimum-jerk flight profiles (` |
shbt-exotic |
Multi-Protocol Spacetime Engineering Co-Simulation | Cross-protocol metric coupling across all 6 phenomena (warp + stasis + translocation + wells + comms + quine), Ford–Roman Quantum Inequality (QI) dark-ledger auditing, Heegaard–Floer symplectic boundary relabeling (T^∂_ij ∈ Sp(2g, Z)), dynamic 5-stage relativistic PDE solvers, 70-gate verification audit + 80-EXT checks, 8 Z3 formal proofs. |
shbt-warp |
Holographic Warp Drive Digital Twin & 3+1D ADM Engine | 3+1D Alcubierre metric foliation (α = 1.0, γ_ij = δ_ij, β^x = -v_s exp(Δ_mod/2) f_SHBT), 4th-Gen SHBT-MAR-312 (4,096.99 kg) / compacted SHBT-MAR-312C (2,257.70 kg, 0.8589 m³) / 5th-Gen SHBT-MAR-312X (950.00 kg, 0.3500 m³, 25.00 EW, 2,482,950c) mass-annihilation reactors, SHBT-IBR-312 magazine (41 cores / 1,378.06 EJ), SHBT-EW-312 Effective Traversable Wormhole System (89.6c → 2,482,950c, 8.957-hour Andromeda M31 transit), 221-gate verification audit, 28 Z3 formal proofs. |
shbt-swarm |
Autonomous Interstellar Von Neumann Swarm Architecture | Macroscopic anyon-braid quine cloning (R_clone^SWARM, 13.418 ms Jump-Cruiser ship replication), non-Hermitian PT-GoI optical template imprinting (54.432 ps), SHBT-STEN cognitive solver & Proof-Carrying CBF Xenobiological Safety Shield (Classes M0, X0, X1, X2), R_univ^SWARM universal 230-space-group matter folding (202.55 Mt planetary synthesis, 11.40 μs Mode 1c direct breeding), 5th-Gen SHBT-MAR-312X autopoietic reactor (950.00 kg, 0.3500 m³, 25.00 EW, 99.200% SHBT-DEC-X, 37,016.7× sensible heat reduction), 104-gate verification audit, 20 Z3 formal SMT proofs. |
shbt-precision provides the mathematical foundation and arbitrary-precision arithmetic substrate governing the entire physical and informational architecture of sys1own/shbt-swarm:
- 512-Bit MPFR Arbitrary-Precision Arithmetic Core:
shbt-precision's 512-bit arithmetic (rug/MPFR, 492-bit mantissa) validates the irrational and transcendental invariants that guarantee swarm self-consistency across astronomical simulation and replication chains. - WZW Conformal Character Convergence
(26, 8, 312)& 346 Dynkin Slots: The affine Lie algebra characters and 346 Dynkin weight routing channels audited inshbt-precisiondirectly defined the exact Casimir denominators(28, 22, 640)used byshbt-swarm's continuous Dynamic Sugawara Framing Compensator (|Δh_vis| ≤ C_braid = 100/φ² = 38.196601 nats,Δ_fr ≡ 0). - Callan–Symanzik Renormalization Group Lock & Sim-to-Real Isomorphism: Rigorous proof of the Callan–Symanzik RG invariants prevents metric drift, securing the Zero-Compute Simulation-to-Reality Isomorphism (
ΔC_post-sim ≡ 0) that allowsshbt-swarmto fabricate macroscopic hardware directly from simulated boundary CFT states without post-synthesis recalibration.
In return, shbt-swarm closes the loop back into shbt-precision by providing empirical and numerical stress-energy verification across all 230 space groups, feeding hardware-in-the-loop closure metrics back into the foundational mathematical tables.
MIT — see LICENSE.