A Dark Material Model of the Universe
A substrate formed of microscopic orbital systems fills the gaps in physics
AI disclosure: I used Claude (Anthropic PBC) for help with content, graphics, and with equations.
About This Project
This is a compilation of the work of many great scientists. I’m a 61 year old Engineering Fellow at Posit PBC - a code first data science company that supports open source. I’m pretty good at software, but a novice when it comes to Physics and this is an open source, solo side project.
In late June, 2025 after seeing one of the first Vera Rubin photos on a hi-res monitor, I intuited that dark matter had to be dark material orbital systems, they had to be in the atom and in space. In my mind’s eye, I saw why it has to be this way. The universe is not weird, and there’s an explanation for cellular dynamics, the gulf stream, magnetic fields.
I eventually found these five amazing theories that when combined paved the last mile:
- Pilot wave hydrodynamics by Bush, Oza, et al.1
- Volovik’s theories from “The Universe in a Helium Droplet”2
- Simeonov’s mapping of the Madelung equations as a fluid form of the Schrödinger equation3
- Khoury’s papers on Dark Matter Superfluidity4
- Larichev & Reznik, Saffman, Aftalion, Sonin, Blatter, Fetter for the math behind modons and vortices5
From this foundation I built the substrate framework and in dialog with Claude, refined it till it first bridged particle physics with general relativity, then accurately modeled perplexing dark energy measurements.
Why you should keep reading
Find perhaps the most intuitive interpretation of mass, quantum potential, spin statistics, quarks, gluons, super conductivity and the nature of the universe.
These results connect through a single chain. One measures the lattice cell size in the substrate using the number of dc1 particles per unit volume with the model for gravity. The other from the mass of the electron with the angle that measures how electromagnetism interacts with the weak nuclear force (the Weinberg angle).
Using the same angle, the substrate model predicts the location of the moraine crust from the last Big Bubble and that matches the latest observations of dark energy, the latest DESI 2 results and resolves a discrepancy in today’s best model predicts of the clumpiness of matter from what we measure.
\sin^2\theta_W \;\xrightarrow{\text{C8}}\; \alpha_{mf} \;\xrightarrow{\text{bridge}}\; \rho_\text{DM} \;\xrightarrow{a_0}\; \text{galactic dynamics} \;\xrightarrow{f(z)}\; \text{dark energy} \;\xrightarrow{S_8}\; \text{structure formation}
Seven domains — electroweak symmetry breaking, quantum mechanics, general relativity, cosmology, galactic dynamics, dark energy evolution, and structure formation — connected through one superfluid, with one geometric backbone (f = 4\pi/(K\sqrt{2})) and one dynamical thread (\alpha_{mf}). Zero adjustable parameters, off by 1-3%. A mathematical model that explains recent data from DESI 2, S_8 tension, tidal dwarf galaxies, the bullet cluster, and Maisie’s galaxy’s unusual redshift.
One mind-bending result is the size of the lattice cell - \xi \approx 100\;\mum | - about the width of a human hair. It is the wake distance of a photon/modon in the dc1 sea, creating a distorted, confusing view of the inside from the outside. This happens because of the stiffness of the lattice layers due to the close-packing of dc1 vortices inside the lattice, pinned by larger orbital systems. The exact mechanics inside the lattice - the mass or spacing dag - the lattice spacing d is <= 7\mu m. A open problem is how to convert the 2D modon equations to 3D in this substrate. Essentially how do we add mass into the bridge equation.
What changes
| Standard physics | Substrate framework | |
|---|---|---|
| Spacetime | Warped geometry (no medium) | Acoustic geometry of a superfluid |
| Time dilation | Fundamental, unexplained | Pressure-dependent clock rate in the medium |
| Gravity | Curvature of spacetime | Ebbing current through counter-spinning vortex boundaries |
| Dark matter | Unknown particle (not yet detected) | The substrate itself (n_1 m_1 = \rho_{DM}) |
| Galaxy rotation curves | Dark matter halos (not yet detected) or modified gravity (MOND, ad hoc) | Boundary parity → quadratic CPR → MOND field equation; a_0 = c\sqrt{G\rho_\text{DM}} |
| Tidal dwarf galaxies | Should be dark-matter-free (no halos) → purely baryonic dynamics | Substrate is the medium, not a halo — TDGs immersed like everything else → same MOND, same a_0, zero new parameters |
| Dark energy | Fine-tuned to 10^{-122} | Zero at equilibrium; observed \Lambda from residual disequilibrium; DESI-confirmed C = 1 |
| S_8 tension | Unexplained 2–3\sigma gap between Planck and weak lensing | Moraine Crust disruption from \mathcal{B}^{-1}: \eta_\text{crust} = 2\alpha_{mf}^2 → S_8 = 0.7788 (zero new parameters) |
| The Big Bang | Singular origin, no location, expanding into itself | A local bubble nucleation in a substrate that boils — one of many, not unique |
| Quantum vacuum | Abstract field with 10^{122}\times too much energy | Superfluid at measured \rho_{DM} |
| Wave–particle duality | Complementarity principle | Vibrating particle in a responsive medium |
| Wave function collapse | Measurement problem (unresolved) | Orbital system reorganization at boundary |
| Spin | “Intrinsic” (no classical analog) | Counter-rotating boundary layer angular momentum |
| Quantization | Postulated (Born rule) | Standing-wave boundary matching (geometric) |
| Entanglement | Nonlocal correlations | Topologically protected vortex channel |
| Mathematical frameworks | Two (QFT + GR), incompatible | One (fluid dynamics at all scales) |
Framework Thesis
If dark matter is really a dark material, composed of at least two particles — nicknamed dc1 (“dark carbon”) and dag (“dark silver”) — much smaller than an electron and of sufficiently different mass, they would have formed a superfluid substrate during the creation event, with enormous rotational energy. By itself dc1 forms vortices; dag pins the lattice that organizes them at the ~100 μm scale. Electrons are persistent dc1 vortex storms — topologically protected circulation patterns, not objects with a center. Atomic nuclei are vortex, orbital system complexes, stabilized by counter-spinning dc1 vortex layers at every turbulent boundary.
There is no mechanism that could have removed this material from the atom, or from space.
How light travels through the substrate, and crossing boundaries without losing energy:
Energy transfer follows the boundary equations for modons — dipole vortex streams that travel long distances against the flow in ocean currents, bracketed by the dag organized lattice. These “dark modons” are vortex/orbital system pair complexes absorbed and emitted between systems, conveying packets of energy with large wakes of interference due to the speed, stiffness, and energetic nature of the superfluid. They are transmitted freely in (and propelled by) the energy of the substrate, crossing boundary layers frictionlessly by flip-flopping spin direction, producing the lowest-energy transmission: Lorentz Invariance.
The substrate framework asks: how do co-rotating systems in low-dissipation environments create and maintain boundary layers, and what is the energy budget of those boundaries? In all analogous macro-scale systems — binary star collisions, vortex streets, Jupiter’s atmosphere, modons, pilot wave pairs — counter-rotating structures form spontaneously at boundaries between co-rotating regions. They self-organize into the lowest-energy configuration consistent with the boundary conditions, and they transport excess energy away as propagating vortex pairs (modons). The math is the same across scales: Euler equations with a vorticity source term at boundaries.
This same mechanism operating at the Planck scale in a dc1 superfluid, organized by the dag lattice reproduces quantum mechanics à la Madelung/Bohm6, where the quantum potential is the reaction force of the counter-rotating boundary layer. Quantized states emerge from wave interference in the dc1/dag medium, not from imposed quantum rules. Quantization is a geometric consequence of oscillatory solutions enclosed by decaying solutions, joined at a boundary.
Spacetime is the acoustic geometry of the dc1/dag substrate. Curvature is the gradient of the current that ebbs through the boundary, then falls in the laminar stream in the substrate by the energetic system inbetween. And Einstein’s equations are understood at a deeper level. The substrate reframes General Relativitity as the low-energy effective theory of quasiparticle propagation.7 Einstein’s equations are the self-consistency condition for the substrate’s response to organized energy. And the features of our universe that ΛCDM takes as given — \Lambda, dark matter, flatness, inflation — emerge from the material properties of the substrate.
The equations show how effective mass disappears, and how a dipole vortex spinning at 0.776c will flow through an superfluid substrate at a constant c.
the substrate forms 2D lattices of orbital systems and counter-spinning vortices
The math behind the model makes a dimensional leap, and has lots of interesting open problems and refinements.
The dc1/dag substrate is dark matter. The “missing mass” in galaxy rotation curves, cluster dynamics, and CMB anisotropies is the substrate itself.
- Collisionless: counter-rotating boundary layers do not interact with electromagnetic modons
- Pressureless on galactic scales: bulk flow is coherent, P_\text{eff} \approx 0 for structure formation
- Self-gravitating: vortices and orbital system rotational energy generates ebbing currents
What it predicts
In cosmology:
| Prediction | Expression | Value | Observation | Domain |
|---|---|---|---|---|
| MOND acceleration a_0 | c\sqrt{G\rho_\text{DM}} | 1.16 \times 10^{-10} m/s² | (1.20 \pm 0.02) \times 10^{-10} m/s² (~3%)8 | Galactic dynamics |
| Flat rotation curves | Derived from boundary parity | Observed universally | — | |
| Baryonic Tully-Fisher M_b \propto v^4 | Derived | M_b \propto v^{3.98 \pm 0.06} | — | |
| a_0/(cH_0) ratio | \sqrt{3\Omega_\text{DM}/(8\pi)} = 0.179 | 0.179 | 0.15% | |
| All Dark Energy is transient (C = 1) | Volovik self-tuning | Equilibrium DE = 0 | DESI DR2 best fit: C = 1.09 | Dark energy evolution |
| TDG rotation curves | Same a_0, same RAR — substrate is medium, not halo | Observed: TDGs follow same relation as normal galaxies | — | |
| Growth suppression S_8 | \eta_\text{crust} = 2\alpha_{mf}^2 | S_8 = 0.7788 | Weak lensing: 0.76–0.79 | Structure formation |
The S_8 result is striking: the Weinberg angle — measured at particle colliders — determines how much the previous cycle’s remnant boundary suppressed galaxy formation in this cycle. The disruption efficiency \eta_\text{crust} = 2\alpha_{mf}^2 = 2(\sin^2\theta_W/(1 - \sin^2\theta_W))^2 = 0.181 arises from a two-step mutual friction coupling, with zero new parameters. This resolves the 2–3\sigma tension between Planck CMB (S_8 = 0.832) and weak lensing surveys.
From one measured input (\sin^2\theta_W = 0.2312) and zero new free parameters:
| Quantity | Predicted | Measured | Discrepancy |
|---|---|---|---|
| Fine structure constant \alpha | 1/135.1 | 1/137.036 | +1.45% |
| Anomalous magnetic moment (g-2)/2 | 0.001178 | 0.001160 | +1.6% |
| Core-boundary asymmetry \eta | 0.03432 | 0.03406 | +0.8% |
Free parameters: The Standard Model + ΛCDM requires ~25 free parameters. The substrate framework has 5 truly unconstrained parameters (M_d, n_d, \delta, B, z_s), of which three (M_d, n_d, \delta) appear in no current prediction and two (B, z_s) describe the previous cycle’s crust — inherently cycle-dependent. All quantitative results above use only measured constants as inputs.
See the bridge equation that ties it all together across seven domains, and the constraint summary.
This model has not been independently validated and has some open problems.
The key identifications:
| Standard Physics | Substrate Framework |
|---|---|
| Quantum vacuum | dc1 vortices organized by the dag-pinned lattice |
| Photon | Modon (counter-rotating vortex dipole) |
| Electron | One effective quantum (\sim 10^9 dc1, mass m_\text{eff} = 1.7 MeV/c^2) a large dc1 vortex circulating at r_\text{eff} = 150 fm with \hbar angular momentum, enclosed in a perturbation envelope with radius \xi \approx 100\;\mum |
| Quantum potential Q | Reaction force from counter-rotating layer |
| Planck’s constant \hbar | 2m \cdot D (diffusion constant of counter-rotating layer) |
| Speed of light c | \hbar/(m_1 \cdot \xi) — ratio of Planck’s constant to dc1 mass times coherence length (Volovik quasiparticle speed in BEC regime) |
| Gravitational constant G | Parametrizes dc1 leak current through boundaries |
| Wave function \psi | \sqrt{\rho} \cdot \exp(iS/\hbar) — amplitude + phase of co-rotating layer |
| Spin | Angular momentum of effective quantum about its axis |
| Chirality | Direction of spin relative to direction of motion |
| Higgs field | Local chirality state of the dc1/dag substrate |
| Measurement | Interaction that couples particle vortex/orbital system to detector field |
| Wave function collapse | Orbital system reorganization upon boundary interaction |
Fermions, including the electron, proton, and quark, are polarized orbital systems or dc1 vortices — with an unbalanced topological charge and an odd number of counter-rotating boundary layers separating their internal co-rotating flow from the external substrate. They repel each other according to the Pauli exclusion principle because two same-state fermions would create an irreconcilable boundary conflict.
Bosons, including the photon, W, Z and Higgs, are balanced opposite-spinning pairs or larger systems that have formed an orbital system complex. They move through the fluid with zero forward momentum — massless — due to the counter-balancing energetic effects and even boundary parity.
Mass is rotational energy - the orbital kinetic energy of substrate particles spinning in organized systems. An electron’s 0.511 MeV is entirely accounted for by one effective quantum (~8.3 × 10⁸ condensed dc1 particles) orbiting at 0.776c.
The speed of light is a substrate property: c = \hbar/(m_1 \cdot \xi), the ratio of Planck’s constant to the dc1 mass times the coherence length. It is the maximum speed at which organized disturbances propagate through the superfluid — set by the medium, not by geometry.
Planck’s constant is the minimal action of a counter-rotating boundary layer: \hbar = 2m \cdot D, where D is the diffusion constant of the counter-spinning eddies. Quantization is not imposed — it emerges because only discrete standing-wave patterns survive the boundary matching between co-rotating interior and decaying exterior.
Photons are modons — counter-rotating vortex dipoles ejected when an orbital system reorganizes across boundary layers. They form from the electron’s coherence dress, ejected as a compact vortex dipole that propagates through the ξ-scale perturbation envelope, which provides the boundary-matching domain for quantization and speed c.
In standard QM, the wavefunction \psi is a probability amplitude with no agreed-upon physical meaning. In the substrate framework, \psi = \sqrt{\rho}\,\exp(iS/\hbar) decomposes into two measurable quantities — \rho is the co-rotating substrate density, and S is the phase of the pilot wave flow. The quantum potential Q emerges from the interaction between the co-rotating flow and the counter-rotating eddies at its boundaries. Nothing is mysterious. Nothing requires interpretation. It’s fluid dynamics.
And this results in these phenomena:
- Free-flowing, Lorentz Invariant energy transport; gravitational lensing based on pressure/energy changes in the substrate; a way to model subtle observed CMB effects without warping space-time or invoking nonlocality — instead creating hidden pathways of wave energy in an active dynamic substrate that reacts like pilot-wave hydrodynamics.
- Gravity acts like an “ebbing” or tidal force applied to boundaries, not individual particles. It applies to the contained mass of the orbital system enclosed by that boundary. The theory predicts gravity’s weakness as a consequence of boundary layer efficiency. And once the boundary layer has been penetrated, it accelerates through the boundary to the next layer.
- QFT is incredibly accurate but combines two dynamical layers into one effective description, producing terms (Q, \mathbf{A}) whose physical origin is obscured because the counter-rotating layer’s degrees of freedom have been integrated out. This framework proposes the microscopic content that would make QFT a complete theory — analogous to how kinetic theory provides the microscopic content behind thermodynamics.
- Tunneling occurs because the counter-rotating layer at the turning point isn’t a perfect wall. It’s a dynamic, fluctuating boundary whose eddies occasionally create momentary gaps.
- Entanglement occurs because a singlet’s topology guarantees signal fidelity through a topologically protected vortex channel — the half-integer winding protects information in transit the same way it protects a qubit in a topological quantum computer. The substrate predicts a drop-off in Bell correlations beyond long distances.
- The spin-statistics connection — fermion/boson distinction, Pauli exclusion, 720° rotation — emerges as a topological consequence of how many counter-rotating boundary layers separate a particle’s internal flow from the external substrate.
- The Higgs mechanism is the local chirality ordering of the dc1/dag substrate, with spontaneous symmetry breaking arising from same-chirality clustering.
- The Aharonov-Bohm effect shows that the substrate field polarity is an almost immeasurable difference in a field that makes a big difference in electron spin.
- The Big Bang as classically conceived is a singular origin with no place to be — spacetime expanding everywhere all at once. The substrate framework’s bang is a local bubble popping in a universe that boils. It pops wherever the pressure has built up enough to nucleate a pocket of the other phase. There is no contradiction between this and what we observe. There is, instead, a different relationship between the observer and the cosmos: we are not surveying the aftermath of a unique creation. We are inside one of the substrate’s normal relaxations, looking outward at the wall that made us, and asking — reasonably, but mistakenly — where it came from. It came from the substrate boiling. That is the whole story. The rest is hydrodynamics.
The DESI satellite may have seen the evidence. The dark energy evolution it measured — phantom behavior at high redshift, a crossing near z \approx 0.5, the return toward w = -1 today — is not what a cosmological constant does. It is what a bubble that passed through the crust of the previous cycle does. The crust parameters (B = 1.88, z_s = 0.63) are the first quantitative fingerprint of whatever came before us. See A Universe That Boils and Dark Energy and the Crust.
The framework now has a consistent story from the Planck scale (substrate particles) through nuclear physics (quark confinement), atomic physics (hydrogen orbitals), condensed matter (conductivity, superconductivity), particle physics (electroweak symmetry breaking, Higgs mechanism, weak force chirality), galactic dynamics (flat rotation curves, the MOND acceleration scale, Tully-Fisher, tidal dwarf galaxies — all from boundary parity with zero new parameters), dark energy evolution (DESI-confirmed C = 1 self-tuning plus crust encounter), and the growth of cosmic structure (S_8 = 0.7788 from the Weinberg angle). Seven domains. One superfluid. Zero adjustable parameters.
Footnotes
Bush, J.W.M. & Oza, A.U., “Hydrodynamic Quantum Analogs,” Annual Review of Fluid Mechanics 52, 2020. A vibrating particle in a responsive medium reproduces quantized orbits, tunneling, and interference — the template for the electron. See also Dagan & Bush, “Hydrodynamic quantum field theory: the free particle,” Comptes Rendus Mécanique, 2020. [R1, R1b]↩︎
Volovik, G.E., The Universe in a Helium Droplet, Oxford University Press, 2003. The single most important source: emergent speed of light (Ch. 7), two-fluid model (Ch. 4–5), gauge fields from vortex cores (Ch. 22–25), and cosmological constant self-tuning (Ch. 29–30). [R2]↩︎
Simeonov, L., “Quantum mechanics as a two-fluid stochastic theory,” arXiv:2509.02868, 2025. Shows that the osmotic velocity of a counter-rotating fluid derives from the HVBK mutual friction force, formally bridging superfluid hydrodynamics to quantum mechanics. [R3]↩︎
Khoury, J., “Dark Matter Superfluidity,” arXiv:1507.01860, 2015; Berezhiani, L. & Khoury, J., “Theory of Dark Matter Superfluidity,” Phys. Rev. D 92, 103510, 2015. Dark matter as a superfluid on galactic scales, with a MOND-like phonon-mediated force and CDM-to-MOND transition at the Landau critical velocity. [R4, R26]↩︎
Larichev, V.D. & Reznik, G.M., “Two-dimensional solitary Rossby waves,” Doklady Akademii Nauk SSSR 231, 1976 (modon boundary matching, K = j_{11}^2+1); Saffman, P.G., Vortex Dynamics, Cambridge, 1992 (dipole dynamics, lattice stability); Aftalion, A. et al., Phys. Rev. A 71, 023611, 2005 (vortex lattice energy functional); Fetter, A.L., Rev. Mod. Phys. 81, 647, 2009 (GP healing length, 1/\sqrt{2} factor). [R5, R6, R7, R8]↩︎
Nelson, E., “Derivation of the Schrödinger Equation from Newtonian Mechanics,” Phys. Rev. 150, 1079, 1966; Bohm, D. & Vigier, J.-P., Phys. Rev. 96, 208, 1954. The formal bridge from superfluid hydrodynamics to QM is completed by Simeonov [R3]. [R18, R19]↩︎
Barceló, C., Liberati, S. & Visser, M., “Analogue gravity,” Living Reviews in Relativity 8, 12, 2005. Any barotropic, irrotational, inviscid fluid produces an acoustic metric formally identical to curved Lorentzian spacetime. The substrate satisfies this. [R13]↩︎
McGaugh, S.S., Lelli, F. & Schombert, J.M., “Radial Acceleration Relation in Rotationally Supported Galaxies,” PRL 117, 201101, 2016. Measured from 2693 data points across 153 galaxies. [R24]↩︎
DESI Collaboration, DR2 combined analysis, 2025. CPL parameterization: w_0 = -0.42 \pm 0.21, w_a = -1.75 \pm 0.58, with ΛCDM excluded at > 2\sigma.↩︎