Two Counter-Rotating Fluids → Quantum Potential
Starting Point: Simeonov’s Framework
Simeonov showed that two coupled fluids reproduce the quantum potential exactly. The substrate framework provides the physical content behind this mathematical result: Fluid 1 is the co-rotating vortex flow of dc1 (the “particle”), and Fluid 2 is the counter-rotating boundary eddies (the source of quantum behavior).
Fluid 1 (co-rotating layer): density \rho_1, velocity \mathbf{v}_1 — carries energy and momentum.
Fluid 2 (counter-rotating layer): density \rho_2, velocity \mathbf{v}_2 — forms as boundary eddies between co-rotating regions, responding to gradients in \rho_1.
Fluid 1 obeys the classical Euler equation with a reaction force from Fluid 2:
\frac{\partial \mathbf{v}_1}{\partial t} + (\mathbf{v}_1 \cdot \nabla)\mathbf{v}_1 = -\frac{1}{\rho_1}\nabla P - \nabla U + \mathbf{F}_\text{reaction}
Fluid 2 diffuses in response to density gradients of Fluid 1:
\mathbf{v}_2 = -D \cdot \nabla(\ln \rho_1) \quad\text{[osmotic velocity]}
where D = \hbar/(2m) is the diffusion constant. The reaction force from Fluid 2 on Fluid 1 is the quantum potential. With \rho = \rho_1 and R = \sqrt{\rho}:
Q = -\frac{\hbar^2}{2m} \cdot \frac{\nabla^2 R}{R}
The quantum potential is not imposed — it emerges from the counter-rotating layer’s response to density curvature. Three cases illustrate the physics:
Near a density maximum (center of an orbital, peak of |\psi|^2): R is large and \nabla^2 R < 0 (concave down), so Q > 0. The quantum potential adds to the effective potential energy, creating a repulsive “quantum pressure” that prevents collapse. This is why electrons don’t spiral into the nucleus. In substrate terms: at the center of a co-rotating region (high \rho_1), the counter-rotating eddies are compressed and their back-pressure pushes outward.
Near a density minimum (node of a wavefunction): R is small and \nabla^2 R > 0 (concave up), so Q is large and negative. But the quantum force is -\nabla Q, not Q itself. Near a node, Q has a sharp negative dip whose gradient points away from the node on both sides — the quantum force repels particles from nodes, maintaining the zero. In substrate terms: at a boundary between co-rotating regions (low \rho_1), the counter-rotating layer is strongest, and the steep gradients push co-rotating flow away from the boundary.
In a uniform region (\rho = \text{constant}): \nabla^2 R = 0, so Q = 0. No quantum effects where there are no boundaries — exactly what the substrate picture predicts.
The Fourth-Order Structure
The quantum force has a distinctive mathematical signature:
\mathbf{F}_\text{reaction} = -\nabla Q = \frac{\hbar^2}{2m} \cdot \nabla\!\left(\frac{\nabla^2 R}{R}\right)
This is a fourth-order spatial derivative of the density — the counter-rotating layer responds to the curvature of the curvature of the co-rotating density. In Simeonov’s framework, this sensitivity emerges naturally: the osmotic velocity \mathbf{v}_2 = -D \cdot \nabla(\ln \rho_1) generates \nabla^2 R / R terms when its divergence and gradient are taken.
The same structure must emerge from the HVBK mutual friction formalism. Starting from the mutual friction force and taking its divergence in steady state should yield:
\nabla \cdot \mathbf{F}_{ns} \propto \nabla^2\!\left(\frac{\nabla^2 R}{R}\right)
which upon integration gives Q. This is the connection point: Simeonov’s abstract “two fluids” become the HVBK co-rotating and counter-rotating components, and the quantum potential becomes the mutual friction reaction force.
A First-Order Equation: the Quaternion Packaging
The route above reaches the Schrödinger equation through the Madelung pair — continuity for \rho_1 and the Euler equation for \mathbf v_1 with Q added — and the framework takes its relativistic spectrum, E^2=\mu^2+c^2p^2, from Volovik (Emergent Speed of Light). Both are second order in space. What the paper had not written is the first-order equation underneath them: the Dirac-type equation that the two-fluid velocity field itself obeys. Danielewski and Sapa’s quaternion quantum mechanics ([R179]) supplies the packaging. Their medium is different — an ideal elastic solid of Planck masses at the Planck length, with the wavefunction a rescaled deformation potential — but the algebraic move carries over, and it lands on the two degrees of freedom this chapter already has: the breath and the circulation.
The wavefunction is the two velocities. Write \Psi=R\,e^{iS/\hbar} as in Madelung. Fluid 1 moves at \mathbf v_1=\nabla S/m and fluid 2 at the osmotic velocity \mathbf v_2=-(\hbar/m)\nabla\ln R (the D\,\nabla\ln\rho_1 of the opening section, with \rho_1=R^2). The gradient of the complex logarithm of \Psi is exactly the pair:
\mathbf w \;\equiv\; -\frac{i\hbar}{m}\,\nabla\ln\Psi \;=\; \mathbf v_1 + i\,\mathbf v_2 .
\Psi is not a probability amplitude with a fluid reading attached afterwards. Its log-gradient is the two-fluid velocity, and the imaginary unit does one job: it keeps the bulk velocity and the boundary-layer velocity from being summed as a single vector. With this substitution the Schrödinger equation becomes an equation for \mathbf w alone,
\frac{\partial \mathbf w}{\partial t} + (\mathbf w\cdot\nabla)\mathbf w \;=\; \frac{i\hbar}{2m}\,\nabla^2\mathbf w \;-\; \frac{1}{m}\nabla U ,
a complex Burgers equation: the Euler equation for the two-fluid velocity with an imaginary viscosity \hbar/2m — the same D=\hbar/2m as the osmotic law. The quantum potential has dissolved: the real part of (\mathbf w\cdot\nabla)\mathbf w together with the real part of the viscous term is exactly -\nabla Q/m. This is first order in time, but still second order in space, and its i is bookkeeping rather than geometry.
The quaternion derivative of a velocity field is (breath, vorticity). Hamilton’s units i,j,k (i^2=j^2=k^2=ijk=-1) let a scalar and a vector share one object, q=s+\mathbf v. Danielewski and Sapa’s Cauchy–Riemann operator is the quaternion gradient \partial = i\,\partial_x + j\,\partial_y + k\,\partial_z. Acting on a full q=s+\mathbf v it gives
\partial q \;=\; -\nabla\!\cdot\mathbf v \;+\; \nabla s \;+\; \nabla\times\mathbf v , \qquad \partial\partial = -\nabla^2 ,
and on a pure velocity field, \partial\mathbf v = -\nabla\!\cdot\mathbf v+\nabla\times\mathbf v. Read against a superfluid, the two pieces are the two degrees of freedom of the cell. The scalar part -\nabla\!\cdot\mathbf v is the breath: by continuity, \partial_t\ln\rho=-\nabla\!\cdot\mathbf v, the compression rate of the anti-phase Compton breath. The vector part \nabla\times\mathbf v is the circulation, carried in this framework by the counter-rotating boundary layer. One operator, one derivative, and the split into compression and twist that Danielewski and Sapa impose by Helmholtz decomposition falls out of the multiplication table. \partial\partial=-\nabla^2 says \partial is the square root of the Laplacian: the operator a first-order equation needs.
The quaternion norm is the energy. For the stiff substrate equation of state P=\rho c^2, set s=c\,\delta\rho/\rho_0. Then
\tfrac12\rho_0\,|q|^2 \;=\; \tfrac12\rho_0|\mathbf v|^2 + \tfrac12\frac{c^2}{\rho_0}\,\delta\rho^2 ,
kinetic plus compressional energy. The quaternion norm |q|^2=s^2+|\mathbf v|^2 is the acoustic energy density. Hurwitz’s theorem says \mathbb R,\mathbb C,\mathbb H,\mathbb O are the only real algebras in which such a norm is multiplicative — the only ways to give a multi-component field a real energy that behaves under products. A field with one scalar and three vector components has exactly one choice, \mathbb H.
The first-order acoustic equation. With s=c\,\delta\rho/\rho_0, the linearized continuity and Euler equations of the bulk collapse to a single quaternion equation:
\frac1c\,\frac{\partial \bar q}{\partial t} \;=\; \partial q \;-\; \nabla\times\mathbf v , \qquad \bar q = s-\mathbf v .
The scalar part is continuity and the vector part is Euler (checked symbolically). For fluid 1 alone — irrotational except at its vortex cores — the last term vanishes and the equation closes: \partial_t\bar q/c=\partial q, a massless first-order equation with one signal speed c, the structure of the Weyl equation. Applying \partial again returns the wave equation, so this is the square root of the acoustic wave equation, not a new dynamics. The conjugate on the left is the helicity: q and \bar q differ by the sign of the vector part, which is the sense of rotation.
The term that spoils closure is the vorticity, and the vorticity is the boundary layer. Fluid 2’s circulation enters the bulk’s first-order equation in exactly one place — the slot \partial has for \nabla\times\mathbf v — and nowhere else. The two-fluid coupling is not added by hand; the algebra has room for exactly it.
The Dirac reading. Assemble the two fluids as a pair of counter-rotating components, q_+ (co-rotating bulk) and q_- (counter-rotating boundary layer). By the equation above each carries its own signal at c; the framework already has this fact from the other side — the Compton breath’s internal disturbance travels at c (Special Relativity), which is what standard physics writes as the Dirac velocity operator having eigenvalues \pm c. Couple the two conservatively at a rate \omega and the pair, in Weyl form, is
i\hbar\,\partial_t\,\psi_\pm \;=\; \mp\, i\hbar c\,\boldsymbol\sigma\!\cdot\!\nabla\,\psi_\pm \;+\; \hbar\omega\,\psi_\mp ,
with \boldsymbol\sigma\!\cdot\!\nabla the Pauli-matrix form of \partial (the Pauli matrices are, up to a factor of i, Hamilton’s units). Squaring gives E^2=(\hbar\omega)^2+c^2p^2 — Volovik’s spectrum with \hbar\omega=\mu=mc^2. So the first-order equation underneath the spectrum the framework imports reads, on the substrate: two counter-rotating fluids, each signalling at c, trading at the Compton frequency.
- The Dirac chirality label \pm is literally the rotation sense, co- and counter-rotating. Chirality means handedness.
- The velocity operator’s \pm c is the rim speed and the breath’s signal speed.
- The mass term \hbar\omega=mc^2 is the exchange between the two components at the Compton frequency \omega_C=mc^2/\hbar. It is the anti-phase breath of Mass as Rotational Energy — energy passing from the contracted all-rotation phase to the expanded all-boundary phase and back, the zitterbewegung at 2\omega_C. A stable particle needs this term Hermitian, so it is the reactive (B') channel of mutual friction, not the dissipative one, that supplies the exchange.
- Spin \tfrac12 appears as the quaternion double cover: a rotation of the medium by the unit quaternion u turns the vector part by u\,\mathbf v\,\bar u, so u and -u are the same rotation and the state needs 720° to return. This is the algebraic face of the odd boundary parity argument of Spin-Statistics: an odd number of sign flips between core and background.
What is derived and what is identified. Three statements here are exact: the operator identity \partial\mathbf v=(-\nabla\!\cdot\mathbf v,\ \nabla\times\mathbf v), the norm |q|^2 as the acoustic energy, and the first-order acoustic equation with the vorticity as its only coupling term. The Dirac step is an identification, not yet a derivation. The coupling rate \omega is set to the Compton frequency because the framework already fixes it there, not because it has been computed from the HVBK coefficients. The open piece is to show that the reactive mutual-friction coefficient B', acting between fluid 1 and fluid 2 at the inner rim, produces exactly \hbar\omega=\alpha_{mf}\,m_\text{eff}\,c^2=mc^2. That is now listed with a definite target in Open Problems.
Where the two programs part company is also worth one line. Danielewski and Sapa put their imaginary unit on the diagonal (i+j+k)/\sqrt3, all three twist axes weighted equally, because an isotropic elastic solid has nothing to prefer one axis. A vortex lattice does: the vortex axis \hat{\mathbf s} that the HVBK force is built around. The substrate’s i has a direction, and it is the one the mutual friction already singles out.
Deriving \hbar from Mutual Friction
In superfluid helium (He-II), the two-fluid equations include a mutual friction force between normal and superfluid components. (For the full laboratory background — where the HVBK equations come from, what is measured, and the microscopic origin of the coefficients — see The HVBK Bridge.) The standard HVBK form is:
\mathbf{F}_{ns} = \frac{B\,\rho_n\,\rho_s}{2\rho}\;\hat{\mathbf{s}} \times \bigl[\hat{\mathbf{s}} \times (\mathbf{v}_n - \mathbf{v}_s - \mathbf{v}_L)\bigr] + \frac{B'\,\rho_n\,\rho_s}{2\rho}\;\hat{\mathbf{s}} \times (\mathbf{v}_n - \mathbf{v}_s - \mathbf{v}_L)
where B, B' are dimensionless mutual friction coefficients, \rho_n and \rho_s are the normal and superfluid densities, \hat{\mathbf{s}} is the unit vector along the vortex line, \mathbf{v}_n is the normal fluid velocity, \mathbf{v}_s is the superfluid velocity, and \mathbf{v}_L is the vortex line velocity.
The substrate identification:
- \mathbf{v}_s \to velocity field of co-rotating vortices of dc1 (Fluid 1 above — the coherent “particle” flow, the superfluid-like component)
- \mathbf{v}_n \to velocity field of counter-rotating dc1 boundary eddies (Fluid 2 above — the diffusive, excitation-carrying layer, the normal-like component)
- \hat{\mathbf{s}} \to direction along the axis of each orbital system (the “vortex line”)
- \mathbf{v}_L \to drift velocity of the orbital system complexes themselves
The B' term (reactive/Hall component) does no work — it only redirects flow. The B term (dissipative component) transfers energy between the two fluids. These two channels become the SU(2)_L and U(1)_Y gauge couplings in the electroweak identification (see Weinberg Angle).
For a superfluid with quantized circulation \kappa_q = h/m_\text{eff}, the effective diffusivity of vortex-mediated transport is:
D_{sf} = \frac{\kappa_q}{4\pi \cdot \alpha_{mf}}
Setting D_{sf} = \hbar/(2m) (the quantum diffusion constant) and substituting \kappa_q = 2\pi\hbar/m_\text{eff}:
\frac{2\pi\hbar}{m_\text{eff} \cdot 4\pi \cdot \alpha_{mf}} = \frac{\hbar}{2m} \qquad\Rightarrow\qquad \frac{\hbar}{2\,m_\text{eff} \cdot \alpha_{mf}} = \frac{\hbar}{2m}
This yields the central mass relation:
\boxed{m_\text{eff} \cdot \alpha_{mf} = m}
The effective mass of the substrate quantum times the mutual friction coupling equals the particle mass. This is constraint C2 — the origin of Planck’s constant in the substrate framework. The quantum of action \hbar is not fundamental; it is 2m \cdot D, where D is the diffusion constant of the counter-rotating boundary layer.
The claim that \hbar is composite has an independent ally, arrived at by a completely different route. Volovik argues from tetrad gravity ([R156]) that \hbar is not a fundamental constant but an element of the Minkowski tetrad — in the Akama–Diakonov–Wetterich reading all diffeomorphism-invariant quantities are dimensionless, \hbar carries dimension of time, \hbar c of length, and c^2 is a ratio of two Planck constants. He works the construction out explicitly for superfluid ^4He, building the “acoustic Planck constants” of the helium vacuum from the atomic mass and density. His route is geometric (the tetrad), the substrate’s is hydrodynamic (2mD, the boundary layer’s diffusivity); both conclude that the quantum of action is a property of the medium. Neither derivation depends on the other — which is exactly what one wants of a claim this radical.
The Mass Hierarchy
Applying the mass relation to the electron and proton:
m_\text{eff} \cdot \alpha_{mf}^{(e)} = m_e = 9.109 \times 10^{-31}\;\text{kg} m_\text{eff} \cdot \alpha_{mf}^{(N)} = m_p = 1.673 \times 10^{-27}\;\text{kg}
Since m_\text{eff} is a substrate property (the same effective quantum in both regimes), the ratio gives:
\frac{\alpha_{mf}^{(N)}}{\alpha_{mf}^{(e)}} = \frac{m_p}{m_e} \approx 1836
The mutual friction coupling is ~1836× stronger in the nuclear regime than the electronic regime. This is not an arbitrary ratio — it is the proton-to-electron mass ratio, emerging from the same boundary physics operating at different scales. In He-3 (where the superfluid has internal structure analogous to the substrate’s particle vortices), \alpha_{mf} varies by orders of magnitude between temperature/pressure regimes, so this large ratio is physically natural.
When \alpha_{mf} = 1 (observed in He-II near the lambda point), the substrate quantum mass equals the particle mass — the particle is “made of” one quantum of circulation. In the electron’s regime (\alpha_{mf} = 0.3008), the effective quantum is heavier than the electron by 1/\alpha_{mf} \approx 3.3, giving m_\text{eff} = 1.70 MeV/c^2.
Kinetic Theory Cross-Check
The superfluid derivation can be cross-checked against kinetic theory. The counter-rotating dc1 particles in the boundary layer move at the inner-scale velocity v_\text{rot,inner} = c\sqrt{2\alpha_{mf}} = 0.776\,c, with a mean free path \lambda \sim 1/(n_1 \cdot \sigma) where \sigma is the dc1-dc1 collision cross section. The kinetic theory diffusivity is:
D_\text{substrate} = \frac{v_\text{rot,inner}}{3\,n_1\,\sigma}
Setting this equal to \hbar/(2m_e):
\frac{v_\text{rot,inner}}{n_1 \cdot \sigma} = \frac{3\hbar}{2\,m_e} \qquad\Rightarrow\qquad n_1 \cdot \sigma = \frac{2\,m_e\,v_\text{rot,inner}}{3\,\hbar} \approx 4.0 \times 10^6\;\text{m}^{-1}
This is constraint C2(b) — a relation linking the dc1 number density and collision cross section to the electron’s reduced Compton wavelength. The product n_1 \sigma sets the “optical depth” of the substrate per unit length: roughly 4 \times 10^6 collisions per meter, or one collision every 0.25\;\mum. With n_1 \approx 6.6 \times 10^{11} m^{-3}, this implies \sigma \sim 6 \times 10^{-6} m^2 — a macroscopically large cross section, consistent with a delocalized BEC where dc1 particles overlap across many coherence lengths. (The 0.25\;\mum figure is a momentum-exchange bookkeeping length for these overlapping, delocalized wavefunctions — not the spacing of localized particles, which sit at about one per 100\;\mum cell; the two lengths describe different things and do not conflict.)
The diffusion constant D = \hbar/(2m) thus has two equivalent substrate expressions — one from superfluid vortex dynamics (\kappa_q/(4\pi\alpha_{mf}), giving the mass relation) and one from kinetic theory (v_\text{rot,inner}/(3n_1\sigma), constraining the collision cross section). Both must hold simultaneously, providing an internal consistency check on the substrate parameters.
The quantum potential established here — the reaction force of the counter-rotating boundary layer — acts on every orbital system at every scale. At the macroscopic scale, the same boundary-crossing mechanism produces gravity: not as curvature of spacetime, but as a net dc1 current leaking through boundaries.