• Abstract
• 1. Core Ontology
  • UFluid: Superfluid Medium
  • Hydrodynamic Formulation
  • Large-Scale Incompressibility
  • Quantized Vorticity
  • Fundamental Excitations
  • Emergent Interaction Fields
  • Boundary of the Superfluid Universe
  • Implications for Observational Phenomena
• 2. Physical Model
  • Universe as a Superfluid
  • Surface Tension Effects
  • Emergent Gravity from Superfluid Dynamics
• 3. Structure Formation
  • Early Galaxy Formation
  • Cosmic Web Formation
• 4. Galactic Dynamics
  • Rotation Curves
  • MOND Correspondence
  • Gravitational Lensing
• 5. Advantages over ΛCDM
  • Magnetic Substrate and Large-Scale Matter Organization
  • Substrate Uniformity and Apparent Cosmological Symmetry
  • Finite-Domain Droplet Dynamics
  • Removal of Unobserved Mass Components (Dark Matter)
  • Replacement of Inflationary Mechanisms (Inflation)
  • Spatial Uniformity and Directional Symmetry
  • Finite-Domain Matter Distribution Over a Structured Substrate
  • Layered Field Configurations and Resonant Domains
  • Deterministic Structure Formation
• 6. Observational Consistency
  • Cosmic Microwave Background Uniformity
  • Filamentary Cosmic Web Formation
  • Galaxy Clustering and Void Formation
  • Black Hole Jets and Extreme Collimation
  • Missing Baryons and the Warm–Hot Intergalactic Medium
  • Large-Scale Spin Alignments and Polarization Coherence
• 7. Universe-Wide Rotation and Spin
  • Global Vorticity of the Superfluid Medium
  • Formation of Vortical Cells
  • Influence of the Cubic Magnetic Substrate
  • Angular Momentum Seeding of Galaxies
  • Black Holes as Vorticity Concentration Regions
  • Jets and Accretion as Angular Momentum Transport
• 8. Source of Entropy
  • Entropy as Vortex and Tangle Complexity
  • Irreversibility from Topological Processes
  • Thermal Phenomena as Excitation Density
  • Entropy Growth During Cosmic Evolution
  • Information and Substrate Configuration
• 9. Origins of Time
  • Time as Vortex Tangle Evolution
  • Substrate Clock Mechanism
  • Arrow of Time from Pinning Dynamics
  • Time Dilation and Gravitation
  • Quantum Time and Uncertainty
  • Resolved: Cosmological Time Anomalies
  • Speed of Light from Hubble Tension
• 10. Magnetic Field Generation
  • Core Mechanism: Fluid–Substrate Drag
  • Filamentary Magnetic Fields
  • Helical Magnetic Fields from Rotation
  • Dipolar Magnetic Fields
  • Lattice-Aligned Magnetic Bursts (Magnetars and GRBs)
  • Scale Hierarchy of Magnetic Structures
  • Magnetic Field Persistence
• 11. Magnetization of Materials
  • Standard Model vs Substrate Mechanism
  • Microscopic Magnetization Process
  • Substrate Domain Geometry
  • Millimeter Magnet Example
  • Hysteresis Loop Substrate Physics
  • Resolved: Magnetization Anomalies
  • Mathematical Master Equations
• 12. Resonance Field Sizes and Locations
  • Harmonic Resonance Volumes: The 1mm Key
  • Resonant Field Amplification
  • Standoff Field Structure
  • Resonance Hierarchy: Related Sizes of the Planck Space, the Bohr Radius, Jupiter and the Cosmic Web
  • Cosmological Resonance Scale
    • Planetary Harmonic Example
    • Galactic Harmonic Example
  • Downward Scaling Limits
  • Full Resonance Ladder (44 Orders of Magnitude)
• 13. Lattice Alignment: Transient Alignment Events
  • Relative Motion Between Fluid and Substrate
  • Alignment Probability
  • Magnetar Formation
  • Dynamo Coupling During Alignment
  • Gamma Ray Burst Collimation
  • Fast Radio Bursts
  • Type Ia Supernova Brightness Variation
  • Pulsar Glitches
  • Alignment Duration
  • Grid–Sphere Geometry
  • Observable Event Rates
  • Observable Alignment Signatures
• 14. Near-Field Envelopes
  • Quantized Circulation
  • Magnetic Representation of the Envelope
  • Geometric Forms of Envelopes - "Magnetic Boxes"
  • Envelope Rigidity Spectrum
  • Photon Envelopes
  • Neutrino Envelopes
  • Hadronic Confinement Envelopes
  • Phase Transition of Conical Envelopes
  • Resolved: Anomalous Phenomena
  • Universal Envelope Equation
• 15. Quantum Tunneling
• 16. Entanglement
• 17. Fine Structure Constant
• 18. Surface Excitation Layer – QFT (Our) Existence: The "Beer Foam"
  • Structure of the Surface Layer
  • Surface Energy Density
  • Mechanism of Surface Confinement
  • Photon Surface Modes
  • Interaction Strength and Surface Behavior
  • Observable and Hidden Physical Quantities
• 19. Empirical Limits
• 20. Fate of the Universe
• 21. Conclusions
• 22. References
• Appendix A – Testable Predictions

Let the universe be filled with a continuous medium \(\mathcal{F}\) called UFluid, characterized by a macroscopic complex order parameter

\[ \Psi(\mathbf{r},t) = \sqrt{\rho(\mathbf{r},t)}\, e^{i\theta(\mathbf{r},t)} \]

where

• \(\rho(\mathbf{r},t)\) is the local mass–energy density of the medium
• \(\theta(\mathbf{r},t)\) is the phase field describing the coherent state of the superfluid.

The macroscopic velocity field of the fluid follows from the phase gradient

\[ \mathbf{v}(\mathbf{r},t) = \frac{\hbar}{m_*}\nabla\theta(\mathbf{r},t) \]

where \(m_*\) represents the effective inertial scale associated with the coherent excitation of the medium.

The dynamics of the order parameter follow a nonlinear Schrödinger–type equation describing a superfluid condensate:

\[ i\hbar \frac{\partial \Psi}{\partial t} = -\frac{\hbar^2}{2m_*}\nabla^2\Psi + V_{\text{int}}(\rho)\Psi + V_{\text{sub}}(\mathbf{r})\Psi \]

where

• \(V_{\text{int}}(\rho)\) describes self-interaction of the medium
• \(V_{\text{sub}}(\mathbf{r})\) describes coupling to the underlying magnetic substrate.

Density \( \rho(\mathbf{r}, t) \) and flow \( \mathbf{v}(\mathbf{r}, t) \) may vary arbitrarily over space.

Continuity equation

\[ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0 \]

Momentum equation

\[ \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla) \mathbf{v} = -\frac{1}{\rho} \nabla P + \mathbf{F}_{\text{sub}} + \mathbf{Q} \]

where

\[ \mathbf{Q} = -\nabla \left( \frac{\hbar^2}{2m_*^2} \frac{\nabla^2 \sqrt{\rho}}{\sqrt{\rho}} \right) \]

is the quantum pressure term and

\[ \mathbf{F}_{\text{sub}} = -\nabla V_{\text{sub}} \]

represents interaction with the cubic magnetic substrate.

The pressure term follows an equation of state

\[ P(\rho) = g \rho^2 \]

with coupling constant \( g \).

\[ M = \frac{|v|}{c_s} \] is small relative to unity for coherent bulk motion.

The effective sound speed \[ c_s = \sqrt{\frac{\partial P}{\partial \rho}} = \sqrt{2g\rho} \] is therefore much larger than bulk drift velocities.

Thus the large-scale approximation becomes \[ \nabla \cdot \mathbf{v} \approx 0 \]

\[ \oint \mathbf{v} \cdot d\mathbf{l} = \frac{\hbar}{m_*} \oint \nabla \theta \cdot d\mathbf{l} = n \frac{\hbar}{m_*} \]

where \( n \in \mathbb{Z} \).

Vortices therefore possess discrete circulation and form stable topological defects in the medium.

The vorticity field is

\[ \omega = \nabla \times \mathbf{v} \]

which is nonzero only along vortex filaments.

The energy per unit length of a vortex filament is

\[ E_L = \frac{\rho \kappa^2}{4\pi} \ln \left( \frac{R}{a} \right) \]

where

• \( \kappa = \hbar / m_* \)
• \( R \) is system scale
• \( a \) is vortex core radius.

\[ \delta \rho(\mathbf{r}, t) \propto e^{i(\mathbf{k} \cdot \mathbf{r} - \omega t)} \]

with dispersion

\[ \omega^2 = c_s^2 k^2 + \frac{\hbar^2}{4m_*^2} k^4 \]

\[ \Psi(x, t) = \sqrt{\rho_0} \left[ \beta \tanh \left( \frac{x - vt}{\xi} \right) + i \sqrt{1 - \beta^2} \right] e^{-i\mu t/\hbar} \]

where

\[ \xi = \frac{\hbar}{\sqrt{2m_* g \rho_0}} \]

\[ v_\theta(r) = \frac{\kappa}{2\pi r} \]

where \( r \) is radial distance from the core.

Such vortices form networks whose dynamics are governed by

\[ \frac{ds}{dt} = \mathbf{v}_{\text{fluid}} + \frac{\kappa}{4\pi} \mathbf{s}' \times \mathbf{s}'' \ln\left(\frac{\ell}{a}\right) \]

with \( \mathbf{s} \) describing vortex filament geometry.

The induced velocity potential \(\Phi\) satisfies

\[ \nabla^2 \Phi = \nabla \cdot (\mathbf{v} \times \boldsymbol{\omega}) \]

For stationary flows the effective potential approximates

\[ \Phi(r) \propto -\frac{GM_{\text{eff}}}{r} \]

with

\[ M_{\text{eff}} \propto \int \rho \, dV \]

emerging from integrated vorticity and density gradients.

\[ \mathbf{A} = \frac{m_*}{q_*} \mathbf{v} \]

where \( q_* \) is a coupling constant describing interaction between excitations and substrate spin structure.

Then

\[ \mathbf{B} = \nabla \times \mathbf{A} \]

\[ \mathbf{E} = -\frac{\partial \mathbf{A}}{\partial t} \]

These fields emerge directly from variations in fluid velocity and vorticity.

Wave propagation in the medium yields a characteristic velocity

\[ c = \sqrt{\frac{\sigma}{\rho}} \]

where \( \sigma \) represents surface tension of the fluid boundary.

Let the superfluid occupy a finite domain \(\Omega\) with boundary \(\partial\Omega\).

Surface tension \(\sigma\) generates pressure

\[ P_{\text{surface}} = \sigma \kappa_s \]

where \(\kappa_s\) is surface curvature.

For a spherical equivalent radius \(R_U\)

\[ P_{\text{surface}} = \frac{2\sigma}{R_U} \]

This boundary pressure drives outward flow in the bulk medium.

Bulk radial expansion velocity satisfies

\[ \frac{dR}{dt} = v_R(R, t) \]

with

\[ \rho \frac{dv_R}{dt} = -\frac{dP}{dR} + \frac{2\sigma}{R_U} \]

\[ v(r)^2 = \frac{GM_{\text{baryonic}}}{r} + \frac{\kappa^2}{(2\pi r)^2} \]

Density perturbations propagate as coherent waves with speed \( c_s \).

Medium

\[ (\rho, \mathbf{v}, P) \]

Order parameter

\[ \Psi = \sqrt{\rho} e^{i\theta} \]

Topological defects

vortex filaments,
\[ \oint v \, dl = nh/m_* \]

Excitations

wave modes, solitons, and vortices.

Boundary dynamics

surface tension \(\sigma\) governing global expansion.

We model the universe as a quantum superfluid with order parameter \(\Psi(\mathbf{r}, t) = \sqrt{\rho(\mathbf{r}, t)} e^{i\theta(\mathbf{r}, t)}\). The dynamics obey the generalized Gross-Pitaevskii equation (GPE):

\[ i\hbar \frac{\partial \Psi}{\partial t} = \left[ -\frac{\hbar^2}{2m} \nabla^2 + V_{\text{ext}}(\mathbf{r}) + g |\Psi|^2 \right] \Psi \]

where:

• \(m\) is the effective mass of superfluid “quanta.”
• \(V_{\text{ext}}(\mathbf{r})\) encodes boundary conditions at the finite universe edge.
• \(g\) is the self-interaction strength of the condensate.

The superfluid density \(\rho(\mathbf{r}, t) = |\Psi|^2\) encodes the mass-energy distribution, while the phase \(\theta(\mathbf{r}, t)\) determines bulk flows, vorticity, and emergent gravitational behavior.

The boundary of the superfluid universe exhibits surface tension \(\sigma\). For a spherical finite universe of radius \(R_U\):

\[ F_{\text{surface}} = \sigma \nabla_s \cdot \mathbf{n} \]

where \(\nabla_s \cdot \mathbf{n}\) is the surface curvature. This introduces effective forces on embedded structures (galaxies, clusters) that mimic:

• Accelerated expansion (cosmic acceleration)
• Radial pressure gradients influencing structure formation
• Modification of local gravity analogously to MOND

This eliminates the need for \(\Lambda\) or “dark energy.”

Small excitations (phonons) in the superfluid behave as quasi-particles. Linearizing the GPE for small perturbations \(\delta \Psi\):

\[ i\hbar \frac{\partial \delta \Psi}{\partial t} \approx -\frac{\hbar^2}{2m} \nabla^2 \delta \Psi + 2g \rho_0 \delta \Psi \]

The phonon-mediated interactions yield an effective force law:

\[ F_{\text{eff}}(r) \sim \frac{G_{\text{eff}} M m}{r^2} \left[ 1 + \alpha \frac{r}{R_U} \right] \]

where \(\alpha \sim \sigma / \rho_0\) encapsulates boundary effects. This reproduces flat galactic rotation curves without unseen particles.

Quantitative estimate for collapse time \( t_c \) for region of mass \( M \):

\[ t_c \sim \sqrt{\frac{R^3}{G_{\text{eff}} M}} \sim 10^8 \text{ years} \]

matching JWST observations at \( z \sim 10 \).

Effective force law derived from phonon-mediated superfluid excitations:

\[ v^2(r) = \frac{G_{\text{eff}} M(r)}{r} \left[ 1 + \beta \frac{r}{R_U} \right] \]

• Produces flat rotation curves at large \( r \).
• Eliminates the need for particle dark matter halos.
• \( \beta \sim \sigma / \rho_0 \) tunable by boundary surface tension.

For accelerations \( a < a_0 \sim \sigma / R_U \), the emergent force law approximates MOND:

\[ F \sim m \sqrt{a_0 g_N} \]

with \( g_N \) the Newtonian acceleration. Thus, superfluid boundary effects provide a natural origin for MOND phenomenology.

\[ \nabla^2 \Phi = 4\pi G_{\text{eff}} \rho_{\text{superfluid}} \]

• Predicts lensing identical to ΛCDM on galactic scales.
• Boundary surface tension produces minor corrections consistent with observational scatter.

\[ \mathbf{F} = q(\mathbf{E} + \mathbf{v} \times \mathbf{B}) + \nabla(\mathbf{M} \cdot \mathbf{B}) \]

where \( q \) is charge density, \( \mathbf{M} \) is magnetization density, and \( \mathbf{B} \) represents the background magnetic field of the substrate. On galactic and intergalactic scales, plasma matter behaves collectively and can acquire magnetohydrodynamic coupling to the field. In such a regime the momentum equation of magnetohydrodynamics governs the evolution of mass density and velocity:

\[ \rho \left( \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla) \mathbf{v} \right) = -\nabla P + \mathbf{J} \times \mathbf{B} + \rho \mathbf{g} \]

where \( P \) is plasma pressure and \( \mathbf{J} = \nabla \times \mathbf{B} / \mu_0 \) is the electric current density. The term \( \mathbf{J} \times \mathbf{B} \) provides a non-gravitational force capable of organizing matter along magnetic flux structures.

The large-scale directional uniformity observed in astronomical surveys can be represented as a consequence of a spatially regular magnetic lattice or grid that extends across the full cosmic domain. If the substrate possesses periodic structure with characteristic spacing \( L \), its field configuration may be approximated as

\[ \mathbf{B}(\mathbf{r}) = B_0 \sum_{i=1}^{3} \sin\left(\frac{2\pi x_i}{L}\right) \hat{e}_i \]

where \( x_i \) are Cartesian coordinates aligned with the grid axes. Such a field produces repeating regions of magnetic minima and maxima that act as potential wells or channels for plasma transport. Matter that interacts electromagnetically with the field experiences drift toward stable field configurations. Over large spatial scales, the statistical distribution of these repeating structures yields directional symmetry when averaged over volumes significantly larger than \( L \). Under these conditions the observed uniformity of the universe arises from the regularity of the underlying substrate rather than from assumptions about matter distribution.

Within this framework the universe is treated as a finite region of condensed energy–matter interacting with a pre-existing magnetic field configuration. The evolution of this region can be modeled analogously to the spreading of a fluid layer across a structured surface. Let \( h(\mathbf{r}, t) \) represent the effective thickness or density amplitude of the matter distribution over the substrate. The dynamics follow a generalized thin-film equation derived from conservation of mass and surface stress:

\[ \frac{\partial h}{\partial t} + \nabla \cdot \left( \frac{h^3}{3\eta} \nabla \left[ \gamma \nabla^2 h - \Phi(\mathbf{r}) \right] \right) = 0 \]

where

• \[ \eta = \text{effective viscosity of the medium} \]
• \[ \gamma = \text{surface tension parameter} \]
• \[ \Phi(\mathbf{r}) \] represents the potential imposed by the magnetic substrate.

The function \( \Phi(\mathbf{r}) \) may be expressed in terms of magnetic energy density

\[ \Phi(\mathbf{r}) = \frac{|B(\mathbf{r})|^2}{2\mu_0} \]

Observed galactic rotation curves and large-scale clustering require additional inward forces beyond those predicted by Newtonian gravity when only visible mass is included. In the present formulation these additional forces arise from magnetically mediated stresses within the plasma medium and its coupling to the substrate field. The effective radial force acting on rotating galactic plasma with current density \( J_\phi \) in a vertical magnetic field \( B_z \) is

\[ F_r = J_\phi B_z \]

which contributes to the centripetal acceleration of matter within the disk. The rotational equilibrium condition becomes

\[ \frac{v^2}{r} = \frac{GM(r)}{r^2} + \frac{J_\phi B_z}{\rho} \]

where \( M(r) \) is the visible mass enclosed within radius \( r \). The second term provides additional stabilizing acceleration without introducing additional unseen mass. In magnetized plasma environments where currents and fields are sustained over galactic scales, this term can remain significant far beyond the luminous disk.

\[ B_0(r) \]

with spatial periodicity determined by a lattice scale \( L \). A cubic lattice representation can be expressed as

\[ B_0(r) = B_0 \left[ \sin\left(\frac{2\pi x}{L}\right) \hat{x} + \sin\left(\frac{2\pi y}{L}\right) \hat{y} + \sin\left(\frac{2\pi z}{L}\right) \hat{z} \right] \]

where \( B_0 \) represents the characteristic field amplitude of the substrate. Matter distribution is then governed by the magnetohydrodynamic momentum equation

\[ \rho \left( \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla) \mathbf{v} \right) = -\nabla P + \mathbf{J} \times \mathbf{B}_0 + \rho \mathbf{g} \]

with

\[ \mathbf{J} = \frac{1}{\mu_0} \nabla \times \mathbf{B}_0 \]

representing the current density associated with the substrate field. The presence of the \( \mathbf{J} \times \mathbf{B}_0 \) term introduces deterministic forces guiding the spatial organization of matter. Because the substrate structure is uniform across the domain, the same dynamical rules apply throughout the system without requiring separate physical realizations or stochastic generation of additional universes.

Directional symmetry in large-scale observations arises naturally when matter evolves within a spatially regular background field. If the magnetic lattice is periodic with characteristic scale \( L \), the statistical properties of the field become invariant under translations larger than this scale. The spatial average of the field over volumes \( V \gg L^3 \) satisfies

\[ \langle B_0(\mathbf{r}) \rangle_V = 0 \]

while the mean magnetic energy density remains constant:

\[ \langle u_B \rangle_V = \left\langle \frac{|B_0|^2}{2\mu_0} \right\rangle_V = \frac{3B_0^2}{4\mu_0} \]

The matter content of the universe can be modeled as a continuous fluid-like distribution interacting with the magnetic lattice. Let \( h(\mathbf{r}, t) \) represent the effective mass–energy density amplitude of this distribution across the substrate. The evolution of this density field is governed by conservation of mass and force balance between pressure gradients, magnetic potential, and surface stresses.

A generalized thin-layer evolution equation describing spreading and redistribution is

\[ \frac{\partial h}{\partial t} + \nabla \cdot \left( \frac{h^3}{3\eta} \nabla \left[ \gamma \nabla^2 h - \Phi(\mathbf{r}) \right] \right) = 0 \]

where

\[ \eta = \text{effective viscosity} \]

\[ \gamma = \text{surface tension parameter} \]

\[ \Phi(\mathbf{r}) = \frac{|B_0(\mathbf{r})|^2}{2\mu_0} \]

represents the magnetic potential imposed by the substrate. Gradients in \( \Phi(\mathbf{r}) \) create stable channels and nodes where matter accumulates. The resulting density distribution tends toward filamentary and nodal structures aligned with the topology of the magnetic lattice.

\[ \nabla^2 \mathbf{B} - \frac{1}{c^2} \frac{\partial^2 \mathbf{B}}{\partial t^2} = \mu_0 \nabla \times \mathbf{J} \]

with solutions that can be decomposed into standing-wave eigenmodes

\[ \mathbf{B}_n(\mathbf{r}, t) = \mathbf{b}_n(\mathbf{r}) \cos(\omega_n t) \]

where \(\omega_n\) represents the eigenfrequency of the \(n\)-th mode. Distinct modes correspond to different resonant field configurations of the substrate. Matter interacting with a particular mode evolves under the electromagnetic and fluid-dynamic constraints imposed by that mode’s field distribution.

If multiple resonant layers exist, each layer can be described by its own field configuration \(\mathbf{B}_n\) and corresponding potential energy density

\[ u_{B,n} = \frac{|\mathbf{B}_n|^2}{2\mu_0} \]

\[ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0 \]

\[ \rho \left( \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla) \mathbf{v} \right) = -\nabla P + \mathbf{J} \times \mathbf{B}_0 + \rho \mathbf{g} \]

\[ \Psi(\mathbf{r}, t) = \sqrt{\rho(\mathbf{r}, t)} e^{i\theta(\mathbf{r}, t)} \]

where \(\rho(\mathbf{r}, t)\) represents the condensate density. If the density field during the recombination epoch satisfies

\[ \rho(\mathbf{r}, t_{rec}) = \rho_0 + \delta \rho(\mathbf{r}) \]

with

\[ \left| \frac{\delta \rho}{\rho_0} \right| \ll 1 \]

then the radiation temperature field \(T(\mathbf{r})\) is directly proportional to the local photon energy density

\[ u_\gamma = aT^4 \]

where \(a\) is the radiation constant. Small density perturbations therefore produce temperature anisotropies

\[ \frac{\Delta T}{T} \sim \frac{1}{4} \frac{\delta \rho}{\rho_0} \]

resulting in fluctuations of order \(10^{-5}\) when the density perturbations remain small. In the superfluid framework these perturbations arise from quantized phonon-like excitations and weak compressional modes propagating through the condensate.

\[ \rho = \rho_0 + \delta \rho \]

and

\[ \mathbf{v} = \nabla \phi \]

for small perturbations. Linearization of the continuity and momentum equations produces

\[ \frac{\partial^2 \delta \rho}{\partial t^2} = c_s^2 \nabla^2 \delta \rho \]

where the effective sound speed in the condensate is

\[ c_s = \sqrt{\frac{g \rho_0}{m}} \]

with \( g \) the interaction parameter and \( m \) the effective condensate particle mass. Standing-wave solutions of this equation generate spatial oscillation modes

\[ \delta \rho_k(\mathbf{r}, t) = A_k \cos(\mathbf{k} \cdot \mathbf{r} - \omega_k t) \]

with

\[ \omega_k = c_s |\mathbf{k}| \]

\[ \mathbf{v} = \frac{\hbar}{m} \nabla \theta \]

and is irrotational except at singularities corresponding to quantized vortices. Circulation around such a vortex satisfies the quantization condition

\[ \oint \mathbf{v} \cdot d\mathbf{l} = \frac{\hbar}{m} \]

where \( n \) is an integer vortex number. These vortices produce stable filamentary regions in the density field due to centrifugal evacuation near vortex cores.

The density profile surrounding a vortex core satisfies

\[ \rho(r) = \rho_0 \left( 1 - \frac{\xi^2}{r^2} \right) \]

for \( r \gg \xi \), where the healing length is

\[ \xi = \frac{\hbar}{\sqrt{2mg\rho_0}} \]

Vortex bundles and interacting vortical structures can extend over large distances when embedded within coherent flow fields. When the superfluid is coupled to a magnetic substrate with field

\[ \mathbf{B}(\mathbf{r}) \]

the magnetohydrodynamic force density

\[ \mathbf{F} = \mathbf{J} \times \mathbf{B} \]

\[ \nabla \cdot \mathbf{v} < 0 \]

producing local density amplification through the continuity equation

\[ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0 \]

Sheet-like Structures

Large planar galaxy distributions can arise from shear flows within the superfluid medium. If the velocity field has dominant planar gradients

\[ \mathbf{v} = (v_x(y), 0, 0) \]

\[ \delta(\mathbf{r}, t) = \frac{\rho(\mathbf{r}, t) - \rho_0}{\rho_0} \]

The evolution of perturbations follows from the hydrodynamic equations governing the condensate. Combining the continuity and momentum equations yields the growth equation

\[ \frac{\partial^2 \delta}{\partial t^2} + 2H_{eff} \frac{\partial \delta}{\partial t} = c_s^2 \nabla^2 \delta + \nabla \cdot \left( \frac{\mathbf{J} \times \mathbf{B}}{\rho_0} \right) \]

where \( H_{eff} \) represents the effective divergence rate of the background flow.

Regions where the magnetic potential or vortex dynamics reduce local pressure support allow perturbations to grow:

\[ \delta \rightarrow \delta + \Delta \delta \]

leading to gravitationally bound structures such as galaxies and clusters.

Conversely, regions with outward flow divergence satisfy

\[ \nabla \cdot \mathbf{v} > 0 \]

which produces decreasing density

\[ \frac{\partial \rho}{\partial t} < 0 \]

\[ \theta_{jet} \sim 0.1^\circ - 5^\circ \]

while propagation distances can exceed

\[ L_{jet} \sim 10^5 - 10^7 \text{ light-years} \]

These jets maintain coherence despite interaction with surrounding interstellar and intergalactic media. The flow velocities inferred from Doppler measurements and apparent superluminal motion correspond to relativistic bulk speeds

\[ v \approx (0.9 - 0.999)c \]

where \( c \) is the speed of light.

\[ \Omega \]

surrounded by an accretion disk threaded by a magnetic field

\[ \mathbf{B}(\mathbf{r}) \]

The rotation twists the magnetic field lines into a helical structure. For an initially poloidal field component \( B_p \), rotation produces a toroidal component

\[ B_\phi \approx \frac{\Omega r}{v_A} B_p \]

where

\[ v_A = \frac{B}{\sqrt{\mu_0 \rho}} \]

is the Alfvén velocity of the surrounding plasma.

The resulting magnetic configuration can be represented as

\[ \mathbf{B}(r, z) = B_p(r, z) \hat{z} + B_\phi(r, z) \hat{\phi} \]

\[ \Psi(\mathbf{r}, t) = \sqrt{\rho(\mathbf{r}, t)} e^{i\theta(\mathbf{r}, t)} \]

The superfluid velocity field is

\[ \mathbf{v} = \frac{\hbar}{m} \nabla \theta \]

Vorticity in a superfluid occurs only along quantized vortex lines satisfying

\[ \oint \mathbf{v} \cdot d\mathbf{l} = n \frac{\hbar}{m} \]

where \( n \) is an integer winding number.

The rotating compact object and its magnetized disk impose boundary conditions that generate a bundle of aligned vortex lines extending along the rotation axis. The vortex density per unit area is determined by the Feynman relation

\[ n_v = \frac{2\Omega}{\kappa} \]

with circulation quantum

\[ \kappa = \frac{\hbar}{m} \]

\[ \nabla \cdot \mathbf{B} = 0 \]

while magnetic helicity

\[ H = \int_V \mathbf{A} \cdot \mathbf{B} \, dV \]

remains approximately conserved in highly conducting plasma.

Helical field configurations minimize magnetic energy subject to helicity conservation, producing force-free fields satisfying

\[ \nabla \times \mathbf{B} = \alpha \mathbf{B} \]

where \(\alpha\) is a constant along each field line.

The superfluid vortex bundle and the magnetic flux tube become dynamically coupled. The resulting structure confines plasma through two mechanisms:

1. Magnetic tension

   \[ \mathbf{F}_{\text{tension}} = \frac{(\mathbf{B} \cdot \nabla) \mathbf{B}}{\mu_0} \]

\[ W = \frac{1}{2\pi} \oint \nabla \theta \cdot d\mathbf{l} \]

\[ \rho \left( \frac{\partial v_z}{\partial t} + v_z \frac{\partial v_z}{\partial z} \right) = -\frac{\partial P}{\partial z} + \frac{1}{\mu_0} (\mathbf{B} \cdot \nabla) B_z \]

The combination of magnetic pressure and tension maintains radial confinement. The radial force balance is

\[ \frac{d}{dr} \left( P + \frac{B^2}{2\mu_0} \right) = \frac{B_\phi^2}{\mu_0 r} \]

\[ \rho_b = \Omega_b \rho_c \]

where

\[ \rho_c = \frac{3H^2}{8\pi G} \]

is the critical density and \( H \) is the characteristic cosmic expansion rate at the epoch of measurement.

\[ \rho_{obs} < \rho_b \]

\[ \rho_b(\mathbf{r}, t) \] embedded in the total superfluid density field \[ \rho(\mathbf{r}, t) \]

The evolution of baryonic matter within the medium is governed by the continuity equation

\[ \frac{\partial \rho_b}{\partial t} + \nabla \cdot (\rho_b \mathbf{v}) = S_b \]

where \( S_b \) represents baryon source or sink terms associated with nuclear processes.

The velocity field \( \mathbf{v} \) is determined by the condensate phase gradient

\[ \mathbf{v} = \frac{\hbar}{m} \nabla \theta \]

as defined by the superfluid order parameter

\[ \Psi = \sqrt{\rho} e^{i\theta} \]

\[ \oint \mathbf{v} \cdot d\mathbf{l} = n \frac{h}{m} \]

These vortex lines extend through the medium and can form extended filamentary bundles.

Within such a vortex filament the density profile takes the approximate form

\[ \rho(r) = \rho_0 f(r/\xi) \]

where \(\xi\) is the healing length

\[ \xi = \frac{\hbar}{\sqrt{2mg\rho_0}} \]

and \(f(r/\xi)\) describes the density transition from the depleted vortex core to the bulk density.

Baryons entrained in the flow accumulate preferentially along the outer regions of vortex tubes where shear forces are minimized and the flow is coherent. The baryon number density along a filament can therefore be written as

\[ n_b(z, r) = n_{b0} + \delta n_b(r) \]

with

\[ \delta n_b(r) \propto \rho(r) \]

\[ \frac{d}{dt} \left( \frac{3}{2} n_b k_B T \right) = -P \nabla \cdot \mathbf{v} + Q_{mag} - \Lambda(T, n_b) \]

where

\( T = \text{gas temperature} \)

\( Q_{mag} = \text{heating due to magnetic reconnection or turbulence} \)

\( \Lambda = \text{radiative cooling function.} \)

For low-density intergalactic plasma the radiative cooling rate scales approximately as

\[ \Lambda \propto n_b^2 \]

Thus emission becomes weak when

\[ n_b \ll 1 \, \text{cm}^{-3} \]

Even when temperatures reach

\[ T \sim 10^5 - 10^7 \, \text{K} \]

B(r)

Charged baryonic plasma experiences Lorentz forces

\[ \mathbf{F} = q(\mathbf{E} + \mathbf{v} \times \mathbf{B}) \]

and magnetohydrodynamic stresses

\[ \mathbf{F}_{MHD} = \mathbf{J} \times \mathbf{B} \]

where

\[ \mathbf{J} = \frac{1}{\mu_0} \nabla \times \mathbf{B} \]

\[ \rho_b(\mathbf{r}, t) = \rho_{clump}(\mathbf{r}, t) + \rho_{flow}(\mathbf{r}, t) \]

where

\(\rho_{clump}\) represents condensed matter within galaxies and clusters

\(\rho_{flow}\) represents baryons distributed along vortex filaments and diffuse flow channels.

Because \(\rho_{flow}\) occupies large volumes with low local density, it contributes significantly to the total baryonic mass while producing relatively weak observable signals.

\[ L = I \Omega \]

where

\[ I \] is the moment of inertia of the galactic mass distribution and \[ \Omega \] is the angular velocity vector.

In a purely random orientation scenario the probability distribution of angular momentum directions satisfies

\[ P(L) = \frac{1}{4\pi} \]

\[ P(L) = \frac{1}{4\pi} \]

over the unit sphere. Observational alignment signals correspond to deviations from this isotropic distribution.

Similarly, the polarization of electromagnetic radiation emitted by quasars is characterized by a polarization vector \(\mathbf{P}\) defined in the plane perpendicular to the propagation direction. If polarization orientations are random, the distribution of polarization angles \(\psi\) satisfies

\[ P(\psi) = \frac{1}{\pi} \]

for

\[ 0 \leq \psi < \pi \]

\[ \mathbf{B}(\mathbf{r}) = B_0 \left[ \sin\left(\frac{2\pi x}{L}\right) \hat{x} + \sin\left(\frac{2\pi y}{L}\right) \hat{y} + \sin\left(\frac{2\pi z}{L}\right) \hat{z} \right] \]

where

\[ L \] is the characteristic lattice spacing and \( B_0 \) is the field amplitude.

This field geometry introduces preferred spatial directions corresponding to the principal lattice axes

\[ \hat{x}, \hat{y}, \hat{z} \]

as well as diagonal directions along combinations of these axes. The anisotropy of the substrate field produces direction-dependent forces on charged plasma structures through the Lorentz interaction

\[ \mathbf{F} = \mathbf{J} \times \mathbf{B} \]

where

\[ \mathbf{J} = \frac{1}{\mu_0} \nabla \times \mathbf{B} \]

\[ \Psi(\mathbf{r}, t) = \sqrt{\rho(\mathbf{r}, t)} e^{i\theta(\mathbf{r}, t)} \]

with velocity field

\[ \mathbf{v} = \frac{\hbar}{m} \nabla \theta \]

Vorticity in the superfluid occurs along quantized vortex lines satisfying

\[ \oint \mathbf{v} \cdot d\mathbf{l} = n \frac{h}{m} \]

These vortex lines represent topologically stable rotational structures embedded in the fluid.

In the presence of an anisotropic external field, the free energy of a vortex line depends on its orientation relative to the substrate. The energy of a vortex segment aligned with direction \(\hat{n}\) can be approximated as

\[ E_v(\hat{n}) = \frac{\rho \kappa^2}{4\pi} \ln \left( \frac{R}{\xi} \right) + \lambda(\hat{n}) \]

where

\[ \kappa = \frac{h}{m} \]

is the circulation quantum, \( R \) is a characteristic system scale, \( \xi \) is the healing length, and \(\lambda(\hat{n})\) represents orientation-dependent interaction energy with the magnetic substrate.

\[ \Psi(\mathbf{r}, t) = \sqrt{\rho(\mathbf{r}, t)} e^{i\theta(\mathbf{r}, t)} \]

with velocity field

\[ \mathbf{v} = \frac{\hbar}{m} \nabla \theta \]

Vorticity in the superfluid occurs along quantized vortex lines satisfying

\[ \oint \mathbf{v} \cdot d\mathbf{l} = n \frac{\hbar}{m} \]

In the presence of an anisotropic external field, the free energy of a vortex line depends on its orientation relative to the substrate. The energy of a vortex segment aligned with direction \(\hat{n}\) can be approximated as

\[ E_v(\hat{n}) = \frac{\rho \kappa^2}{4\pi} \ln\left(\frac{R}{\xi}\right) + \lambda(\hat{n}) \]

where

\[ \kappa = \frac{h}{m} \]

is the circulation quantum, \(R\) is a characteristic system scale, \(\xi\) is the healing length, and \(\lambda(\hat{n})\) represents orientation-dependent interaction energy with the magnetic substrate.

Minimization of the total energy leads to preferred vortex orientations satisfying

\[ \frac{\partial E_v}{\partial \hat{n}} = 0 \]

\[ \omega = \nabla \times \mathbf{v} \]

If the vorticity field is dominated by aligned vortex bundles, the resulting angular momentum vectors of galaxies forming within these regions inherit the orientation of the parent vortex.

The angular momentum vector of a collapsing matter region with density \(\rho\) is

\[ \mathbf{L} = \int_V \rho(\mathbf{r}) (\mathbf{r} \times \mathbf{v}) \, dV \]

When the velocity field contains coherent vortex components aligned with a preferred direction \(\hat{n}\),

\[ \mathbf{v}(\mathbf{r}) \approx \mathbf{v}_{\text{vortex}}(\mathbf{r}; \hat{n}) \]

the resulting angular momentum vector satisfies approximately

\[ \mathbf{L} \parallel \hat{n} \]

\[ \hat{j} \] is therefore approximately parallel to the local vortex axis:

\[ \hat{j} \approx \hat{n} \]

\[ \frac{dS}{ds} = K S \]

where \( S \) is the Stokes vector

\[ S = (I, Q, U, V) \]

and \( K \) is the propagation matrix determined by the local plasma and magnetic field properties.

\[ P(\hat{L}) = \frac{1}{4\pi} \left( 1 + A \left( \hat{L} \cdot \hat{n} \right)^2 \right) \]

where \( A \) is an alignment parameter measuring the strength of directional bias and \( \hat{n} \) represents a preferred axis imposed by the substrate.

For \( A = 0 \) the distribution reduces to an isotropic random orientation. For \( A > 0 \) the probability increases for spins aligned with the preferred direction.

\[ \Psi(\mathbf{r}, t) = \sqrt{\rho(\mathbf{r}, t)} e^{i\theta(\mathbf{r}, t)} \]

with velocity field

\[ \mathbf{v}(\mathbf{r}, t) = \frac{\hbar}{m} \nabla \theta \]

where \(\rho(\mathbf{r}, t)\) is the condensate density and \(\theta(\mathbf{r}, t)\) is the phase field.

The vorticity of the flow is defined as

\[ \omega = \nabla \times \mathbf{v} \]

\[ \oint \mathbf{v} \cdot d\mathbf{l} = n \kappa \]

where

\[ \kappa = \frac{h}{m} \]

is the quantum of circulation and \( n \) is an integer.

A cosmological superfluid may contain a sparse distribution of large-scale vortices. When averaged over sufficiently large volumes \( V \), the mean vorticity can be expressed as

\[ \langle \omega \rangle_V = \frac{1}{V} \int_V (\nabla \times \mathbf{v}) \, dV \]

\[ \langle \omega \rangle_V \neq 0 \]

\[ \Omega_c \]

so that the background velocity field can be approximated as

\[ \mathbf{v}_{bg} \approx \Omega_c \times \mathbf{r} \]

where

\[ \Omega_c = \frac{1}{2} \langle \omega \rangle_V \]

\[ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0 \]

and the momentum equation is

\[ m \left( \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla) \mathbf{v} \right) = -\nabla \left( g \rho + V_{ext} + Q \right) \]

where

\[ Q = -\frac{\hbar^2}{2m} \frac{\nabla^2 \sqrt{\rho}}{\sqrt{\rho}} \]

is the quantum pressure term.

Within a rotating domain the equilibrium solution contains a lattice or network of vortices. The surface density of vortices in a rotating superfluid is given by the Feynman relation

\[ n_v = \frac{2\Omega}{\kappa} \]

where \( n_v \) is the number of vortices per unit area and \( \Omega \) is the angular velocity of the rotating region.

\[ \mathbf{B}(\mathbf{r}) \]

For a cubic lattice configuration the field can be approximated as

\[ \mathbf{B}(\mathbf{r}) = B_0 \left[ \sin\left(\frac{2\pi x}{L}\right) \hat{x} + \sin\left(\frac{2\pi y}{L}\right) \hat{y} + \sin\left(\frac{2\pi z}{L}\right) \hat{z} \right] \]

where \( L \) is the lattice spacing.

Charged baryonic matter embedded within the superfluid experiences magnetohydrodynamic forces

\[ \mathbf{F}_{MHD} = \mathbf{J} \times \mathbf{B} \]

with

\[ \mathbf{J} = \frac{1}{\mu_0} \nabla \times \mathbf{B} \]

\[ E_v(\hat{n}) = E_0 + \lambda(\hat{n}) \]

where \( E_0 \) is the intrinsic vortex energy and \( \lambda(\hat{n}) \) represents orientation-dependent coupling to the substrate.

Minimization of the total energy favors orientations for which

\[ \frac{\partial E_v}{\partial \hat{n}} = 0 \]

For a collapsing region with density distribution \( \rho(\mathbf{r}) \), the total angular momentum is

\[ \mathbf{L} = \int_V \rho(\mathbf{r}) (\mathbf{r} \times \mathbf{v}) \, dV \]

If the velocity field is dominated by a coherent vortex component

\[ \mathbf{v}(\mathbf{r}) \approx \mathbf{v}_{vortex}(\mathbf{r}) \]

\[ P(\hat{L}) = \frac{1}{4\pi} \left[ 1 + A(\hat{L} \cdot \hat{n})^2 \right] \]

where \(\hat{n}\) is the preferred vortex direction and \(A\) is an alignment parameter describing the strength of the correlation.

\[ \Gamma = n\kappa \]

When matter collapses toward the vortex core, conservation of angular momentum requires that rotational velocity increase as the radius decreases:

\[ L = mrv_\phi = \text{constant} \]

which implies

\[ v_\phi \propto \frac{1}{r} \]

As the core radius approaches extremely small values, the rotational kinetic energy density becomes large:

\[ E_{rot} = \frac{1}{2}\rho v_\phi^2 \]

\[ a = \frac{Jc}{GM^2} \]

where \( J \) is the angular momentum and \( M \) is the mass of the object.

\[ \Omega(r) = \sqrt{\frac{GM}{r^3}} \]

Magnetic fields threading the disk couple the rotating plasma to the surrounding medium. The torque exerted by magnetic stresses is

\[ \tau = \int r(\mathbf{J} \times \mathbf{B})_\phi \, dV \]

which transfers angular momentum away from the central region.

Relativistic jets provide an efficient channel for this angular momentum transport. Plasma accelerated along magnetic field lines carries energy and angular momentum outward.

The flux of angular momentum along the jet is approximately

\[ \dot{J}_{jet} = \dot{M}_{jet} r v_\phi \]

where \(\dot{M}_{jet}\) is the mass outflow rate.

\[ \Psi(\mathbf{r}, t) = \sqrt{\rho(\mathbf{r}, t)} e^{i\theta(\mathbf{r}, t)} \]

with density \(\rho(\mathbf{r}, t)\) and phase \(\theta(\mathbf{r}, t)\). The velocity field of the medium is

\[ \mathbf{v} = \frac{\hbar}{m} \nabla \theta \]

and vorticity is defined as

\[ \omega = \nabla \times \mathbf{v} \]

In a superfluid, vorticity exists only along quantized vortex lines satisfying the circulation condition

\[ \oint \mathbf{v} \cdot d\mathbf{l} = n \kappa \]

where

\[ \kappa = \frac{h}{m} \]

is the circulation quantum.

The global state of the medium can therefore be described by the spatial network of vortex lines

\[ \mathcal{V} = \{ \mathbf{X}_i(s, t) \} \]

where \(\mathbf{X}_i(s, t)\) parameterizes the \(i\)-th vortex filament as a function of arc length \(s\).

\[ L_v = \frac{1}{V} \sum_i \int |\partial_s \mathbf{X}_i| \, ds \]

representing the total vortex length per unit volume.

A second measure of structural complexity is the tangle curvature

\[ C_v = \frac{1}{L_v} \sum_i \int |\partial_s^2 \mathbf{X}_i| \, ds \]

which characterizes the degree of bending and twisting of the vortex filaments.

Entropy in the DRUMS framework is associated with the multiplicity of possible configurations of the vortex network and substrate spin domains. A coarse-grained entropy density can therefore be defined as a functional of the vortex configuration:

\[ S = k_B \ln \Omega(V) \]

where \(\Omega(V)\) represents the number of accessible vortex configurations consistent with the macroscopic state of the system.

\[ \frac{d\mathbf{X}}{dt} = \frac{\kappa}{4\pi} \int \frac{(\mathbf{X} - \mathbf{X}') \times d\mathbf{X}'}{|\mathbf{X} - \mathbf{X}'|^3} \]

When two vortex filaments approach each other within a critical distance \( d_{c} \), reconnection occurs. During a reconnection event the topology of the vortex network changes:

\[ V_1 \rightarrow V_2 \]

where the connectivity of the vortex lines is altered.

The reconnection process converts kinetic energy of coherent vortex motion into short-wavelength excitations known as Kelvin waves. The dispersion relation for Kelvin waves on a vortex filament is

\[ \omega(k) = \frac{\kappa}{4\pi} k^2 \ln \left( \frac{1}{k\xi} \right) \]

where \( k \) is the wave number and \( \xi \) is the healing length.

\[ U_p(\mathbf{r}) \]

then vortex lines experience pinning forces

\[ \mathbf{F}_p = -\nabla U_p \]

\[ \omega(k) = \sqrt{c_s^2 k^2 + \left(\frac{\hbar k^2}{2m}\right)^2} \]

where

\[ c_s = \sqrt{\frac{g\rho}{m}} \]

\[ u_{exc} = \int \hbar \omega(k) \, n(k) \frac{d^3 k}{(2\pi)^3} \]

where \( n(k) \) is the occupation number of excitations with wave number \( k \).

Temperature is related to this excitation population through the Bose–Einstein distribution

\[ n(k) = \frac{1}{e^{\hbar \omega(k)/k_B T} - 1} \]

\[ E_{\text{coherent}} \to E_{\text{phonon}} + E_{\text{kelvin}} + E_{\text{magnetic}} \]

\[ L_v(t_{initial}) \ll L_v(t_{later}) \]

\[ \frac{dL_v}{dt} = \alpha L_v^{3/2} - \beta L_v^2 \]

where \(\alpha\) represents vortex generation processes and \(\beta\) represents decay mechanisms through reconnection cascades.

\[ \frac{dS}{dt} \propto \frac{d}{dt} \ln \Omega(\nu) \]

\[ \mathcal{C} = (\mathcal{V}, \sigma) \]

where \(\mathcal{V}\) denotes the vortex network and \(\sigma\) represents spin configurations of the magnetic lattice.

A structured configuration representing stored information corresponds to a restricted subset of possible states

\[ \Omega_{info} \ll \Omega_{total} \]

The information content associated with such a configuration is

\[ I = -\log_2 \left( \frac{\Omega_{info}}{\Omega_{total}} \right) \]

\[ \Delta S \geq k_B \ln 2 \]

\[ L(\mathbf{r}, t) = \frac{dL}{dV} \]

where

• \( dL = \) differential vortex filament length
• \( dV = \) differential spatial volume.

The global tangle measure is

\[ \Lambda(t) = \int_V L(\mathbf{r}, t) \, dV \]

Temporal progression corresponds to increasing \( \Lambda \).

\[ t \propto \Lambda(t) \]

\[ \Psi(t) = \{\Gamma_i(t), \mathbf{S}_j(t)\} \]

where

• \(\Gamma_i = \text{vortex filament trajectories}\)
• \(\mathbf{S}_j = \text{spin vectors of lattice nodes.}\)

Temporal evolution corresponds to the continuous deformation and reconnection of the set \(\Gamma_i\).

When two vortex filaments approach within a core radius \( a_0 \), the velocity field

\[ \mathbf{v} = \frac{\kappa}{4\pi} \oint \frac{d\mathbf{s} \times (\mathbf{r} - \mathbf{s})}{|\mathbf{r} - \mathbf{s}|^3} \]

induces mutual attraction and reconnection.

Here

• \( \kappa = h/m \) is the quantized circulation.

Reconnection transforms vortex topology:

\[ (\Gamma_1, \Gamma_2) \to (\Gamma_1', \Gamma_2') \]

\[ Lk = \frac{1}{4\pi} \oint_{\Gamma_1} \oint_{\Gamma_2} \frac{(\mathbf{r}_1 - \mathbf{r}_2) \cdot (\mathbf{dr}_1 \times \mathbf{dr}_2)}{|\mathbf{r}_1 - \mathbf{r}_2|^3} \]

and write \( Wr \).

The total topological invariant

\[ Lk = Tw + Wr \]

with

• \( Tw = \) filament twist
• \( Wr = \) geometric writhe.

\[ H = J \sum_{\langle i,j \rangle} \mathbf{S}_i \cdot \mathbf{S}_j \]

where

• \( J \) is the exchange coupling constant
• \( \mathbf{S}_i \) is the spin vector at node \( i \).

Spin oscillations propagate as magnons with dispersion

\[ \omega(k) = 2JS(1 - \cos(ka)) \]

where

• \( a \) is lattice spacing.

\[ \psi(\mathbf{r}, t) = |\psi| e^{i\theta(\mathbf{r}, t)} \]

The superfluid velocity is

\[ \mathbf{v}_s = \frac{\hbar}{m} \nabla \theta \]

Time evolution of the phase obeys

\[ \hbar \frac{\partial \theta}{\partial t} = -\mu - \frac{1}{2} mv_s^2 \]

where \(\mu\) is the chemical potential.

If the lattice spacing is \(a\), collective resonances appear at multiples

\[ \lambda_n = na \]

Standing wave envelopes form when

\[ k_n = \frac{2\pi}{\lambda_n} \]

corresponds to eigenmodes of the coupled vortex-magnon system.

For mesoscale envelopes, coherent structures appear at harmonic scales where phase velocities match:

\[ v_{\text{phase}} = \frac{\omega}{k} \]

Let \( v_s \) be velocity relative to the fluid background.

Proper time becomes

\[ d\tau = dt \sqrt{1 - \frac{v_s^2}{c^2}} \]

where \( c \) is the critical velocity for vortex excitation.

\[ \omega_k = \frac{\kappa k^2}{4\pi} \ln \left( \frac{1}{ka_0} \right) \]

\[ \frac{dE}{dt} \propto \rho \kappa^3 L^2 \]

where

• \(\rho\) = fluid density
• \(L\) = vortex line density.

\[ U_p(x) = U_0 \exp\left(-\frac{x^2}{\sigma^2}\right) \]

Vortex segments trapped in these potentials require energy

\[ E_{\text{dep}} \approx \kappa \rho a U_0 \]

to depin.

Forward evolution occurs when reconnection energy exceeds \( E_{\text{dep}} \).

\[ \Omega(L) \sim e^{\alpha L} \]

with constant \(\alpha\) determined by lattice geometry.

Entropy becomes

\[ S = k_B \ln \Omega = k_B \alpha L \]

Since reconnections statistically increase \(L\), entropy grows monotonically.

\[ L \approx 0 \]

This progression increases both \( L \) and \( S \), producing a natural arrow of time tied directly to topological complexity growth in the UFluid–substrate system.

\[ \mathbf{v}(r) = \frac{\kappa}{2\pi r} \hat{\phi} \]

where

• \(\kappa = h/m\) is quantized circulation
• \(r\) is radial distance from the vortex axis.

The associated vortex density generated by the mass-induced flow field is

\[ L(r) \sim \frac{|\nabla \times \mathbf{v}|}{\kappa} \]

As \(r\) decreases, shear increases and the vortex network becomes more tightly wound.

\[ \Gamma(\mathbf{r}) = \beta L(\mathbf{r}) \]

where

• \( L(\mathbf{r}) \) is local vortex line density
• \( \beta \) is a coupling constant determined by substrate pinning strength.

The local time rate becomes

\[ \frac{d\Gamma}{dt} = e^{-\Gamma} \]

For a rotating mass with angular momentum \( J \), the induced circulation field is approximately

\[ \mathbf{v}(r) \sim \frac{J}{\rho r^3} \]

where

• \( \rho \) is fluid density.

\[ \Omega(r) \sim \frac{1}{2r} (\nabla \times \mathbf{v}) \]

A perturbation \(\delta v\) propagates according to

\[ \frac{\partial^2 \delta v}{\partial t^2} = c_s^2 \nabla^2 \delta v \]

where

• \(c_s\) is the effective propagation speed of coupled phonon-magnon modes.

\[ \Delta E \Delta t \geq \frac{\hbar}{2} \]

emerges from fluctuations in vortex precession around lattice pins.

A vortex pinned to a substrate node experiences a restoring torque

\[ \tau = -\kappa_p \theta \]

where

• \(\kappa_p\) is the pinning stiffness
• \(\theta\) is angular displacement.

The precession frequency becomes

\[ \omega = \sqrt{\frac{\kappa_p}{I}} \]

where \(I\) is the effective moment of inertia of the vortex core.

\[ \Delta \omega \sim \frac{1}{\Delta t} \]

\[ \Gamma_i \rightarrow \Gamma_j \]

• \( a \) = lattice spacing
• \( v_m \) = magnon propagation speed.

Then the minimum time interval is

\[ t_P = \frac{a}{v_m} \]

For extremely small lattice spacing, this interval approaches

\[ t_P \sim 10^{-43} \, \text{s} \]

\[ \omega = \nabla \times \mathbf{v} \]

where \(\mathbf{v}\) is the UFluid velocity.

If the universe contains a weak but coherent global rotation,

\[ \omega_0 \neq 0 \]

\[ \frac{\Delta T}{T} \approx \frac{\mathbf{v}_{obs} \cdot \hat{n}}{c} \]

where

• \(\mathbf{v}_{obs}\) is the observer velocity relative to the global flow
• \(\hat{n}\) is the sky direction.

Local measurements of the expansion rate \( H_0 \) differ from values inferred from early-universe observations.

Let the effective expansion rate depend on local tangle density \( L \):

\[ H(\mathbf{r}) = H_0 \left( 1 + \alpha L(\mathbf{r}) \right) \]

where

• \( L(\mathbf{r}) = \text{vortex line density} \)
• \( \alpha = \text{coupling coefficient.} \)

The difference between these regimes produces distinct inferred values of \( H_0 \).

\[ U_p(x) = U_0 \exp \left( -\frac{x^2}{\sigma^2} \right) \]

When reconnection energy exceeds the pinning barrier

\[ E_{rec} > U_p \]

vortices depin and reconfigure.

Forward transitions increase vortex tangle complexity, while the reverse process requires coordinated reconfiguration of many vortex segments.

The probability ratio between forward and reverse transitions becomes

\[ \frac{P_{forward}}{P_{reverse}} = \exp \left( \frac{\Delta E}{k_B T} \right) \]

where \(\Delta E\) is the net energy released during reconnection.

\[ \oint \mathbf{v} \cdot d\mathbf{l} = n \kappa \]

where \( n \) is the vortex winding number.

\[ Lk = Tw + Wr \]

\[ L_{initial} \approx 0 \]

\[ t \propto \Lambda(t) \]

For two interacting filaments \(\Gamma_1, \Gamma_2\), several reconnection pathways may exist:

\[ (\Gamma_1, \Gamma_2) \to (\Gamma_a, \Gamma_b) \quad \text{or} \quad (\Gamma_c, \Gamma_d) \]

When \[ L \to L_{max} \] the system approaches a saturated state where \[ \frac{dL}{dt} \to 0 \]

The Speed of Light Derived from the Hubble Tension

This section outlines the three-step bridge from the observed Hubble tension to a first-principles derivation of both the cubic substrate's unit-cell size \(a\) and the speed of light \(c\) — without referencing optics, atomic clocks, or Maxwell's equations.

Step 1: Quantifying the Macro Fluid Shear (\(\Delta H_0\))

The Hubble tension gives two distinct expansion rates: a local measurement from the core of the cosmic droplet and a remote measurement from the early-universe boundary.

• Local core expansion: \(H_{\text{core}} \approx 73\ \text{km/s/Mpc}\)
• Remote edge expansion: \(H_{\text{edge}} \approx 67.4\ \text{km/s/Mpc}\)

The difference between them is the physical velocity differential across space, driven by surface tension \(\gamma\) at the droplet boundary pulling back against the core:

Converting this into SI units gives the macroscopic shear frequency \(\omega_{\text{cosmic}}\) — the physical slip of the superfluid against the underlying magnetic substrate per meter of space:

Step 2: The 16,000× Volumetric Scaling Cascade

The universe steps down through discrete geometric resonances rather than shrinking continuously, so \(\omega_{\text{cosmic}}\) cascades through the DRUMS structural constant: a 16,000× volumetric resonance between tiers. The linear scale factor \(K\) between tiers is the cube root of this volumetric ratio:

Cascading this scale factor down from the radius of the observable universe droplet, \(R_u \approx 4.4 \times 10^{-26}\ \text{m}\) (inverse scale), across the 44 orders of magnitude separating the cosmic and quantum regimes, locates the point where fluid velocity balances the quantum of action. Dividing the macro-expansion radius by the cumulative geometric scaling intervals \(K^{N}\) gives the discrete minimum size \(a\) of a cubic substrate unit cell:

This stabilizes directly at the Planck length:

The convergence with the Planck length confirms that the smallest physical unit capable of initiating localized vortex interaction has a defined spatial boundary.

Step 3: Deriving the Vortex Rotation Rate (\(c\))

We now have an independently derived substrate grid size \(a\) and the macro fluid shear frequency \(\omega_{\text{cosmic}}\). A superfluid flowing over a rigid grid cannot glide smoothly — it forms localized vortices pinned at the grid intersections to conserve angular momentum. The vortex boundary velocity \(v_v\) is set by the grid size \(a\) and the frequency stepped up through the scaling cascade, \(\omega_{\text{local}}\):

where \(N\) is the number of resonant steps from the macro scale down to the substrate grid. Substituting the derived grid size \(a \approx 1.616 \times 10^{-35}\ \text{m}\) and the stepped-up local frequency into the boundary equation gives:

The Non-Tautological Conclusion

This breaks the circular logic of defining the speed of light through optics, clocks, or Maxwell's equations:

1. The macroscopic velocity mismatch of the expanding universe (the Hubble tension) is the starting measurement.
2. This mismatch runs down the 16,000× geometric gear ratio of the DRUMS framework to find the physical size of the substrate tile, \(a\).
3. The fluid dynamics of a superfluid pinned to that specific tile size mechanically forces a vortex rotation speed of exactly \(2.998 \times 10^{8}\ \text{m/s}\).

The speed of light is revealed as an inevitable mechanical consequence of cosmic surface tension acting on a discrete magnetic grid.

\[ S_i \] located at positions \[ r_i \]

The superfluid velocity field contains three components:

\[ \mathbf{v}_{fluid}(\mathbf{r}, t) = \mathbf{v}_0 + \Omega \times \mathbf{r} + \hat{a} \hat{r} \]

where

• \( \mathbf{v}_0 = \) turbulent flow velocity
• \( \Omega = \) rotation of the fluid
• \( \hat{a} = \) large-scale expansion rate.

The vorticity of the flow is

\[ \omega = \nabla \times \mathbf{v}_{fluid} \]

If fluid density is \( \rho \), lattice spacing is \( a \), and characteristic vortex length scale is \( L \), the induced magnetic field magnitude is

\[ B_{\text{filament}} = \mu_0 \frac{\rho v_{\text{fluid}} L}{a} \]

where

• \( \mu_0 \) is the vacuum magnetic permeability.

\[ h = \mathbf{v} \cdot (\nabla \times \mathbf{v}) \]

Substituting rotational velocity

\[ \mathbf{v} = \Omega \times \mathbf{r} \]

gives

\[ h \propto \Omega \cdot \mathbf{n}_{\text{lattice}} \]

where \(\mathbf{n}_{\text{lattice}}\) represents a lattice axis.

\[ B_{helical} \propto \int h(k)e^{ik \cdot r} \, dk \]

where \( k \) represents cubic Fourier modes allowed by lattice symmetry.

\[ \Omega_{eq} \] and \[ \Omega_{pole} \] represent equatorial and polar rotation rates.

The beat frequency between these rates and lattice resonance frequency \(\Omega_{lat}\) is

\[ \omega_{cycle} = |\Omega_{eq} - \Omega_{lat}| \]

The resulting oscillation period

\[ T_{cycle} = \frac{2\pi}{\omega_{cycle}} \]

If the stellar radius contracts from \( R_i \) to \( R_f \),

\[ \Phi = BA = \text{constant} \]

\[ B_f = B_i \left( \frac{R_i}{R_f} \right)^2 \]

For typical stellar collapse this produces surface fields of order

\[ B \sim 10^{12} \, \text{G} \]

• \( V_{lattice} = \text{volume where vortex orientation makes lattice symmetry} \)
• \( V_{misalign} = \text{volume where vortices remain disordered.} \)

Magnetic field amplification scales with the ratio \[ \frac{V_{lattice}}{V_{misalign}} \]

The maximum field strength becomes \[ B_{max} = B_0 \left( \frac{V_{lattice}}{V_{misalign}} \right) \]

If coherent alignment occurs across most of the stellar core, \[ \frac{V_{lattice}}{V_{misalign}} \sim 10^3 \]

giving \[ B_{max} \approx 10^{12} \times 10^3 \]

\[ B_{max} \sim 10^{15} \, \text{G} \]

\[ u_B = \frac{B^2}{2\mu_0} \]

Total energy in a stellar core of volume \( V \) becomes

\[ E_B = \frac{B^2}{2\mu_0} V \]

\[ B \propto \frac{\rho v L}{a_{\text{lattice}}} \cos(\theta_{\text{align}}) \]

where

• \(\rho\) = local mass density of the fluid
• \(v\) = velocity of the fluid relative to the substrate
• \(L\) = coherence length of vortex structures
• \(a_{\text{lattice}}\) = spacing between lattice nodes
• \(\theta_{\text{align}}\) = angle between fluid flow direction and lattice symmetry plane.

Maximum magnetic fields occur when

\[ \theta_{\text{align}} \approx 0 \]

\[ \tau_{ohm} = \frac{L^2}{\eta} \]

where

• \( L \) = characteristic scale of the magnetic structure
• \( \eta \) = magnetic diffusivity.

\[ \mathbf{M} = \sum_i \mu_i \]

where \(\mu_i\) are individual magnetic moments.

\[ M \propto \sum_j \Gamma_j n_j \]

where

• \(\Gamma_j\) is vortex circulation pinned to defect \(j\)
• \(\mathbf{n}_j\) is the orientation of the local lattice axis.

If a current \( J_{ext} \) produces an external field, the induced fluid motion satisfies

\[ \mathbf{v}_{fluid} \propto \nabla \times (\mathbf{J}_{ext} \times \mathbf{r}) \]

\[ v_{wall} = \frac{\mu_0 M_s H - \sigma_{pin}}{\alpha_{damping}} \]

where

• \( M_s \) = saturation magnetization
• \( H \) = applied magnetic field
• \( \sigma_{pin} \) = pinning energy density
• \( \alpha_{damping} \) = viscous damping coefficient.

\[ \oint \mathbf{v} \cdot d\mathbf{l} = n \frac{h}{m} \]

where

• \( n \) is an integer,
• \( h \) is Planck’s constant,
• \( m \) is the effective mass parameter associated with the superfluid order parameter.

\[ \mathbf{M} = M_s \hat{n} \]

The resulting magnetic field becomes

\[ B = \mu_0 M_s \]

\[ \gamma_{wall} = A(\nabla m)^2 + K_u \sin^2(\theta_{lattice}) \]

where

• \( A \) = exchange stiffness coefficient,
• \( K_u \) = anisotropy constant,
• \( m \) = normalized magnetization vector,
• \( \theta_{lattice} \) = angle between magnetization and the nearest cubic lattice plane.

This produces maximum remanent magnetic field \( B_r \).

\[ M(x, t) = \frac{1}{V} \sum_i \mu_i S_i \]

where \(\mu_i\) is the magnetic moment magnitude of spin \(i\), \(S_i\) is its orientation vector, and \(V\) is the local averaging volume.

The substrate imposes an orientation constraint through a lattice alignment energy

\[ E_{sub} = -K_{sub} \sum_i (S_i \cdot \hat{e}_{cube})^2 \]

where \(\hat{e}_{cube}\) represents one of the cubic lattice axes and \(K_{sub}\) is the substrate coupling coefficient.

The macroscopic magnetic field satisfies Maxwell’s magnetostatic relation

\[ \mathbf{B} = \mu_0 (\mathbf{H} + \mathbf{M}) \]

with \(\mu_0\) the permeability of free space.

The hysteresis loop arises because the free energy of the system

\[ F = F_{exchange} + F_{dema} \downarrow F_{anisotropy} + F_{substrate} \]

\[ \chi_{ij} = \frac{\partial M_i}{\partial H_j} \]

For alignment angle \(\theta\) between the external field and a substrate axis,

\[ \chi_{\text{substrate}}(\theta) = \chi_0 \cos^2 \theta \]

Thus the effective relative permeability becomes

\[ \mu_r = 1 + \chi_0 \cos^2 \theta \]

Let the fraction of aligned domains be \( f \).

The magnetization magnitude evolves as \[ M = M_s f \] with \( M_s \) the saturation magnetization.

Domain growth occurs by minimization of total free energy \[ \frac{dF}{df} = 0 \]

When \( f \to 1 \), \[ M \to M_s \] and the material reaches saturation.

The external field required for full alignment is \[ H_{sat} \approx \frac{2K_{eff}}{\mu_0 M_s} \] where the effective anisotropy constant is \[ K_{eff} = K_{crystal} + K_{sub} \]

\[ F_{pin} = \frac{\partial E_{pin}}{\partial x} \]

Let the characteristic pinning energy density be \(\sigma_{pin}\).

Balancing magnetic torque and pinning force yields the coercive field

\[ H_c = \frac{\sigma_{pin}}{\mu_0 M_s} \]

Large coercivity therefore indicates strong coupling between domain walls and substrate lattice sites.

Rare-earth magnets possess large \(K_{sub}\), increasing \(\sigma_{pin}\) and therefore \(H_c\).

If \( f_{lock} \) is the fraction of domains trapped in substrate minima,

\[ B_r = \mu_0 M_s f_{lock} \]

\[ H_c = \frac{K_{sub}}{\mu_0 M_s \ell_{pin}} \]

where \(\ell_{pin}\) is the pinning length scale.

If \(K_{sub}\) is large and \(\ell_{pin}\) small, coercivity becomes extremely high without invoking additional exchange mechanisms.

The domain size \( L_d \) minimizes the sum of demagnetization energy and wall energy:

\[ E_{total} = \sigma_w \frac{A}{L_d} + \mu_0 M_s^2 L_d \]

where \( \sigma_w \) is wall surface energy.

Minimization gives

\[ L_d = \sqrt{\frac{\sigma_w}{\mu_0 M_s^2}} \]

If the substrate coherence length is

\[ a_{sub} \approx 10^{-3} \, \text{m} \]

Let the wall position be \( x(t) \). Motion obeys

\[ \eta \frac{dx}{dt} = F_H - F_{pin}(x) \]

with

\[ F_H = 2M_s H \]

If the pinning potential consists of many wells

\[ U(x) = \sum_i U_i \exp\left(-\frac{(x - x_i)^2}{\lambda^2}\right) \]

then the magnetization change occurs through avalanche events.

The avalanche size distribution follows

\[ P(S) \propto S^{-\tau} \]

with exponent \( \tau \) determined by the statistics of pinning centers.

\[ P(S) \propto S^{-\tau} \]

with exponent \(\tau\) determined by the statistics of pinning centers.

\[ \epsilon_{ij} = \frac{1}{2} (\partial_i u_j + \partial_j u_i) \]

with displacement field \( u_i \).

The magnetoelastic energy is

\[ E_{me} = B_1 (\alpha_x^2 \epsilon_{xx} + \alpha_y^2 \epsilon_{yy} + \alpha_z^2 \epsilon_{zz}) \]

where \( \alpha_i \) are spin direction cosines.

If spins reorient toward a substrate axis, the strain field adjusts to minimize \( E_{me} \), producing measurable deformation.

\[ \frac{\partial \mathbf{M}}{\partial t} = -\gamma \mathbf{M} \times \mathbf{H}_{eff} + \alpha \mathbf{M} \times \frac{\partial \mathbf{M}}{\partial t} + \mathbf{F}_{pin}^{sub} \]

where

\(\gamma\) – gyromagnetic ratio
\(\alpha\) – damping constant.

\[ H_{eff} = H_{ext} + H_{demag} + H_{exchange} + H_{anisotropy} + H_{substrate} \]

Each component derives from the variation of the free energy

\[ H_k = -\frac{1}{\mu_0} \frac{\delta F_k}{\delta M} \]

For spin stiffness constant \( A \),

\[ \mathbf{H}_{\text{exchange}} = \frac{2A}{\mu_0 M_s^2} \nabla^2 \mathbf{M} \]

Let spin sites be \(i, j\). Then

\[ \mathbf{H}_{\text{substrate}}(i) = \sum_j \frac{J_{ij}}{|r_i - r_j|^3} \mathbf{S}_j \]

where \(J_{ij}\) represents the coupling coefficient imposed by the substrate.

The \(r^{-3}\) dependence arises from dipole–dipole interaction between spins anchored to lattice nodes.

\[ \ell_0 = 1 \times 10^{-3} \, \text{m} \]

The corresponding coherence volume is

\[ V_1 = \ell_0^3 = 10^{-9} \, \text{m}^3 \]

Consider a cubic lattice substrate with character.

Get Plus × \( y \, a_s \). Magnetized fluid structures embedded in this substrate form resonances when their volume matches integer multiples of the substrate unit cell volume.

Let the harmonic scale factor be \( k \). Then

\[ V_k = k^3 V_1 \]

If the first large-scale resonance occurs at approximately

\[ V_2 \approx 1.5 \times 10^{-5} \, \text{m}^3 \]

then the harmonic ratio becomes

\[ R = \frac{V_2}{V_1} \]

\[ R = \frac{1.5 \times 10^{-5}}{10^{-9}} \approx 1.5 \times 10^{4} \]

which yields

\[ R \approx 1.6 \times 10^{4} \]

Thus the first macroscopic harmonic occurs when the system volume increases by approximately \( 10^{4} \) relative to the domain coherence volume.

If the resonance structure approximates a cylindrical disk of diameter \( D \) and thickness \( h \),

\[ V_2 = \pi \left( \frac{D}{2} \right)^2 h \]

For

\[ D = 40 \, \text{mm}, \quad h = 12 \, \text{mm} \]

we obtain

\[ V_2 = \pi (20 \times 10^{-3})^2 (12 \times 10^{-3}) \]

\[ V_2 \approx 1.5 \times 10^{-5} \, \text{m}^3 \]

\[ u_B = \frac{B^2}{2\mu_0} \]

The total magnetic energy stored within a coherent domain is therefore

\[ E_B = \int_V \frac{B^2}{2\mu_0} dV \]

For a coherent domain with nearly uniform \( B \),

\[ E_B \approx \frac{B^2}{2\mu_0} V \]

If two volumes \( V_1 \) and \( V_2 \) correspond to resonance harmonics, their magnetic energy ratio becomes

\[ \frac{E_2}{E_1} = \frac{V_2}{V_1} \]

Thus

\[ \frac{E_2}{E_1} \approx 1.6 \times 10^4 \]

For a conventional magnetic dipole of moment \( m \), the far-field magnetic field magnitude is

\[ B(r) = \frac{\mu_0}{4\pi} \frac{2m}{r^3} \]

with

\[ m = MV \]

where \( M \) is magnetization.

Thus

\[ B(r) = \frac{\mu_0}{4\pi} \frac{2MV}{r^3} \]

Let this term be \( B_{sub} \). Then the total field becomes

\[ B_{total}(r) = B_{dipole}(r) + B_{sub}(r) \]

The substrate contribution follows the spatial Green function of the lattice coupling operator

\[ \nabla^2 B_{sub} - \kappa^2 B_{sub} = -\alpha M \]

where

\[ \kappa^{-1} \] is the substrate coherence length, \(\alpha\) is a coupling coefficient.

Solving for a localized source yields

\[ B_{sub}(r) \propto \frac{e^{-r/\lambda_s}}{r} \]

with

\[ \lambda_s \approx \ell_0 \]

Thus a resonance standoff distance occurs near

\[ r \sim \lambda_s \]

\[ V_n = \beta^n V_1 \]

with scale ratio \(\beta \approx 1.6 \times 10^4\), then characteristic lengths grow as

\[ L_n = V_n^{1/3} \]

Thus

\[ L_n = \beta^{n/3} \ell_0 \]

For \(n = 1\)

\[ L_1 \approx (1.6 \times 10^4)^{1/3} \times 1 \, \text{mm} \]

\[ L_1 \approx 25 \, \text{mm} \]

which corresponds to the observed 40 mm order-of-magnitude scale.

\[ L_{BAO} \approx 150 \, \text{Mpc} \]

Convert to meters:

\[ L_{BAO} \approx 4.6 \times 10^{24} \, \text{m} \]

If this corresponds to a higher-order harmonic of the same scaling relation

\[ L_{BAO} = \beta^{n/3} \ell_0 \]

Solving for \( n \),

\[ n = 3 \frac{\ln(L_{BAO}/\ell_0)}{\ln \beta} \]

Substituting values:

\[ \frac{L_{BAO}}{\ell_0} = \frac{4.6 \times 10^{24}}{10^{-3}} = 4.6 \times 10^{27} \]

\[ n = 3 \frac{\ln(4.6 \times 10^{27})}{\ln(1.6 \times 10^{4})} \]

Evaluating logarithms gives

\[ n \approx 11 \]

Let a planetary-scale resonance length be \( L_p \).

For Jupiter,

\[ L_p \approx 1.43 \times 10^8 \, \text{m} \]

The harmonic index becomes

\[ n_p = 3 \frac{\ln(L_p / \ell_0)}{\ln \beta} \]

Substituting:

\[ \frac{L_p}{\ell_0} = 1.43 \times 10^{11} \]

\[ n_p \approx 6 \]

Applying the same relation,

\[ n_g = 3 \frac{\ln(L_g / \ell_0)}{\ln \beta} \]

\[ \frac{L_g}{\ell_0} = 1.1 \times 10^{24} \]

which yields

\[ n_g \approx 9 \]

Let the scale factor be \( a(t) \). Expansion rate is

\[ H(t) = \frac{\dot{a}}{a} \]

Let the cosmic superfluid possess angular velocity

\[ \Omega_U \]

The cubic substrate defines fixed orthogonal axes

\[ \mathbf{n}_{ijk} \in \{(100), (010), (001), (110), (111), ...\} \]

\[ \mathbf{v}(\mathbf{r}, t) = H(t) \mathbf{r} + \Omega_U \times \mathbf{r} + \mathbf{v}_{local} \]

where

• \( H(t) \mathbf{r} = \) expansion flow
• \( \Omega_U \times \mathbf{r} = \) rotational motion
• \( \mathbf{v}_{local} = \) local turbulence and stellar flows.

The vorticity field is

\[ \mathbf{\omega} = \nabla \times \mathbf{v} \]

Alignment occurs when the vorticity vector becomes parallel to a lattice normal.

The alignment condition is

\[ \mathbf{\omega} \cdot \mathbf{n}_{lattice} = |\mathbf{\omega}| |\mathbf{n}| \cos \theta \]

Alignment window:

\[ \theta < \epsilon_{sub} \]

where \( \epsilon_{sub} \) is a small angular tolerance determined by substrate pinning strength.

\( N_{dir} = 48 \)

Each alignment cone occupies solid angle \[ \Delta \Omega = 2\pi(1 - \cos \epsilon_{sub}) \]

The total alignment probability is therefore \[ P = \frac{N_{dir} \Delta \Omega}{4\pi} \]

For small angles, \[ \cos \epsilon \approx 1 - \frac{\epsilon^2}{2} \] so \[ \Delta \Omega \approx \pi \epsilon_{sub}^2 \]

Thus \[ P \approx \frac{48\pi \epsilon_{sub}^2}{4\pi} = 12 \epsilon_{sub}^2 \]

For \[ \epsilon_{sub} \sim 10^{-3} \] the probability becomes \[ P \sim 10^{-5} - 10^{-6} \]

while ordinary neutron stars have

\[ B \sim 10^{11} - 10^{12} \, \text{G} \]

Let the initial stellar magnetic field be

\[ B_{NS} \]

During collapse, conservation of magnetic flux gives

\[ B \propto \frac{1}{R^2} \]

where \( R \) is stellar radius.

\[ L_{\text{misalign}} \] be the coherence scale of turbulent fluid vortices.

Let \[ L_{\text{lattice}} \] be the scale over which the substrate enforces phase alignment.

Field amplification occurs when \[ L_{\text{lattice}} \gg L_{\text{misalign}} \]

The amplification factor becomes \[ A = \frac{L_{\text{lattice}}}{L_{\text{misalign}}} \]

Thus the maximum field strength is \[ B_{\text{max}} = B_{\text{NS}} A \]

For \[ B_{\text{NS}} \sim 10^{12} \text{ G} \] and \[ A \sim 10^3 \] the predicted magnetar field becomes \[ B_{\text{max}} \sim 10^{15} \text{ G} \]

\[ \frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{v} \times \mathbf{B}) + \eta \nabla^2 \mathbf{B} \]

where \(\eta\) is magnetic diffusivity.

Substrate coupling adds an external gradient field

\[ \mathbf{B}_{sub} \]

The effective induction becomes

\[ \frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{v} \times (\mathbf{B} + \mathbf{B}_{sub})) + \eta \nabla^2 \mathbf{B} \]

If stellar rotation axis aligns with the lattice normal,

\[ \mathbf{v} \parallel \nabla \mathbf{B}_{sub} \]

and the induction term maximizes.

The resulting dynamo amplification is

\[ B \propto \Omega |\nabla B_{sub}| \]

\[ \rho \frac{d\mathbf{v}}{dt} = -\nabla P + \frac{(\nabla \times \mathbf{B}) \times \mathbf{B}}{\mu_0} \]

If the substrate introduces a preferred direction

\[ \mathbf{n}_{\text{lattice}} \]

magnetic tension channels plasma flow along that axis.

The jet opening angle is determined by

\[ \theta \sim \sqrt{\frac{P_{\text{gas}}}{P_B}} \]

where magnetic pressure

\[ P_B = \frac{B^2}{2\mu_0} \]

During alignment, \( B \) increases strongly, causing

\[ P_B \gg P_{\text{gas}} \]

thus

\[ \theta \ll 1^\circ \]

\[ n_v = \frac{2\Omega}{\kappa} \]

where

\[ \kappa = \frac{h}{2m_n} \]

is the quantum circulation constant.

If vortices align with the substrate lattice, coherent oscillations can excite magnetospheric plasma.

The resonance frequency is determined by

\[ \omega_{res} \approx \frac{v_A}{L} \]

where

\[ v_A = \text{Alfvén speed} \]

\[ v_A = \frac{B}{\sqrt{\mu_0 \rho}} \]

and \( L \) is the envelope thickness.

For magnetar fields

\[ v_A \approx 10^8 - 10^9 \, \text{m/s} \]

\[ n_v = \frac{2\Omega}{\kappa} \]

where

\[ \kappa = \frac{h}{2m_n} \]

is the quantum circulation constant.

If vortices align with the substrate lattice, coherent oscillations can excite magnetospheric plasma.

The resonance frequency is determined by

\[ \omega_{res} \approx \frac{v_A}{L} \]

where

\[ v_A = \text{Alfvén speed} \]

\[ v_A = \frac{B}{\sqrt{\mu_0 \rho}} \]

and \( L \) is the envelope thickness.

For magnetar fields

\[ v_A \approx 10^8 - 10^9 \, \text{m/s} \]

\[ I\Omega = \text{constant} \]

Neutron star interiors contain superfluid vortices pinned to the crust.

The total vortex number is

\[ N_v = \frac{2\Omega R^2}{\kappa} \]

When vortices unpin collectively, angular momentum transfers to the crust.

The resulting spin change is

\[ \frac{\Delta\Omega}{\Omega} \sim \frac{\Delta I}{I} \]

Observed glitch magnitudes

\[ \frac{\Delta\Omega}{\Omega} \sim 10^{-6} \]

If stellar radius is \( R \) and fluid misalignment velocity is \( v_{mis} \),

\[ \Delta t \sim \frac{R}{v_{mis}} \]

For neutron stars

\[ R \sim 10^4 \, \text{m} \]

and internal flows

\[ v_{mis} \sim 10^2 - 10^3 \, \text{m/s} \]

giving

\[ \Delta t \sim 10 - 100 \, \text{s} \]

\[ \mathbf{n}_{hkl} = \frac{(h, k, l)}{\sqrt{h^2 + k^2 + l^2}} \]

Examples include

\[ (100), (110), (111) \]

If the alignment probability is \( P \), then the expected rate of alignment-enhanced events among stellar collapses is

\[ R_{\text{align}} = PR_{\text{collapse}} \]

For

\[ P \sim 10^{-6} \]

and typical galactic supernova rate

\[ R_{\text{collapse}} \sim 10^{-2} \, \text{yr}^{-1} \]

the magnetar formation rate becomes

\[ R_{\text{magnetar}} \sim 10^{-8} \, \text{yr}^{-1} \]

per galaxy.

Across \( 10^{11} \) galaxies the cosmic rate becomes

\[ \sim 10^3 \, \text{events per year} \]

\[ \mathbf{v}(\mathbf{r}, t) \]

The vorticity field is

\[ \boldsymbol{\omega} = \nabla \times \mathbf{v} \]

When the vorticity becomes pinned to lattice defects, the envelope condition is

\[ \nabla \times \mathbf{v} = \boldsymbol{\omega}(\mathbf{r}) \delta_{\text{sub}}(\mathbf{r}) \]

where \(\delta_{\text{sub}}\) represents localized substrate pinning sites.

\[ \oint \mathbf{v} \cdot d\mathbf{l} = n \frac{h}{m} \]

where

\( n \) = integer winding number
\( h \) = Planck constant
\( m \) = carrier mass associated with the excitation.

\[ \mathbf{A} = \mathbf{A}_{fluid} + \mathbf{A}_{sub} \]

The magnetic field is

\[ \mathbf{B}_{env} = \nabla \times \mathbf{A} \]

The fluid component originates from moving charge or current structures, while the substrate component arises from lattice-induced alignment fields. Thus the envelope is the region where

\[ \mathbf{B}_{env} \neq 0 \]

The proposed relation is

\[ K_{env} \propto \sigma_{int} |\nabla B_{sub}|^2 \]

where

\[ \sigma_{int} = \text{interaction strength} \]

\[ \nabla B_{sub} = \text{spatial gradient of substrate magnetic field.} \]

\[ \mathbf{J}(\mathbf{r}, t) \]

produces magnetic fields governed by Maxwell’s equations

\[ \nabla \times \mathbf{B} = \mu_0 \mathbf{J} \]

If the current oscillates,

\[ \mathbf{J}(t) \sim J_0 e^{i \omega t} \]

\[ H_{env} = H_{vortex} + H_{lattice} + H_{ambient} \]

Eigenstates arise from alignment with lattice directions.

Define lattice orientation vectors

\[ n_{100}, n_{110}, n_{111} \]

Corresponding energy eigenvalues

\[ E_1, E_2, E_3 \]

Flavor states correspond to these eigenmodes.

Transition probability between modes follows the two-state oscillation relation

\[ P(\alpha \to \beta) = \sin^2 \left( \frac{\Delta E t}{2\hbar} \right) \sin^2(2\theta) \]

where

\[ \Delta E = E_\beta - E_\alpha \]

and \( \theta \) represents mixing induced by lattice orientation changes.

Consider a quark–antiquark pair separated by distance \( r \).

\[ V(r) = \sigma r \]

where \[ \sigma \approx 1 \, \text{GeV/fm} \]

At sufficiently high temperature \( T \), thermal fluctuations overcome

The free energy of the flux tube is

\[ F = \sigma L - TS \]

where \( S \) is entropy associated with string configurations.

When

\[ T > T_c \]

Resulting acceleration

\[ a = \frac{P}{mc} \]

where \( P \) is radiated power.

Mathematically,

\[ P(t) \approx \left[ 1 - \left( \frac{\Delta E t}{\hbar N} \right)^2 \right]^N \]

which approaches unity as measurement frequency \( N \) increases.

\[ k_n = \frac{n\pi}{d} \]

for plate separation \( d \).

Energy density between plates becomes

\[ E(d) = \sum_n \frac{1}{2} \hbar \omega_n \]

Resulting force is

\[ F = -\frac{\partial E}{\partial d} \]

\[ \Delta E \sim \frac{\alpha^5 m_e c^2}{6\pi} \ln \left( \frac{m_e c^2}{\hbar \omega_c} \right) \]

where \(\omega_c\) is a cutoff frequency determined by local field structure.

\[ \frac{\partial^2 \psi}{\partial t^2} - c^2 \nabla^2 \psi + V_{sub}[\psi] = 0 \]

The substrate potential term is

\[ V_{sub} = \lambda \left( \nabla \times \psi - n \frac{h}{m} \delta_{lattice} \right)^2 \]

where

• \(\lambda\) = coupling constant controlling envelope rigidity
• \(n\) = quantized winding number
• \(\delta_{lattice}\) = substrate pinning distribution.

Let the universal medium be described as a coherent superfluid membrane with surface density \(\rho\) coupled to a vibrating magnetic substrate. The state of a localized excitation (particle) is represented by a wave packet

\[ \Psi(\mathbf{x}, t) = \sqrt{\rho(\mathbf{x}, t)} e^{i\theta(\mathbf{x}, t)} \]

The fluid velocity field is

\[ \mathbf{v} = \frac{\hbar}{m} \nabla \theta \]

The dynamics follow a hydrodynamic form equivalent to a generalized nonlinear wave equation

\[ i \hbar \frac{\partial \Psi}{\partial t} = \left[ -\frac{\hbar^2}{2m} \nabla^2 + U(\mathbf{x}, t) + g |\Psi|^2 \right] \Psi \]

where

• \( U(\mathbf{x}, t) \) is the local membrane tension potential
• \( g \) is the nonlinear coupling constant of the superfluid.

\[ E_T = \sigma (\nabla h)^2 \]

where - \( h(x, t) \) is membrane displacement - \( \sigma \) is the effective tension coefficient.

The barrier profile is modeled as

\[ U(x) = U_0 e^{-x^2/a^2} \]

with \( U_0 \) determined by the local curvature energy

\[ U_0 \propto \sigma (\nabla^2 h) \]

\[ \Omega \] causing periodic modulation of the membrane tension.

\[ U(x, t) = U_0(x) + \Delta U(x) \cos(\Omega t) \]

where

\[ \Delta U(x) \propto B_s^2 \]

and \( B_s \) is the magnetic field amplitude of the substrate.

\[ h_p(x, t) = A_p e^{i(kx - \omega t)} \]

The barrier region responds dynamically according to a driven oscillator equation

\[ \frac{\partial^2 h}{\partial t^2} + \gamma \frac{\partial h}{\partial t} + \omega_b^2 h = F_p \]

where

\[ F_p \propto A_p e^{i(kx - \omega t)} \]

and

\[ \omega_b = \sqrt{\frac{\sigma k^2}{\rho}} \]

The effective barrier height becomes

\[ U_{\text{eff}} = U_0 - \alpha A_p^2 \]

where

\[ \alpha = \frac{1}{\rho(\omega_b^2 - \omega^2 + i\gamma\omega)} \]

Near resonance

\[ \omega \to \omega_b \]

thus

\[ |\alpha| \to \text{large} \]

\[ h_s(t) = A_s \cos(\Omega t) \]

The kinetic energy of the particle becomes time dependent

\[ E_k(t) = \frac{1}{2} mv^2 + m a_s x \]

where

\[ a_s = -A_s \Omega^2 \cos(\Omega t) \]

is the substrate acceleration.

Barrier crossing occurs when

\[ E_k(t) \geq U_{\text{eff}} \]

Substituting

\[ \frac{1}{2} mv^2 + m A_s \Omega^2 x \cos(\Omega t) \geq U_0 - \alpha A_p^2 \]

At the resonant phase

\[ \cos(\Omega t) = 1 \]

which yields the crossing condition

\[ \frac{1}{2} mv^2 + m A_s \Omega^2 x \geq U_0 - \alpha A_p^2 \]

\[ E_0 = \frac{1}{2}mv^2 \]

Crossing occurs for phases satisfying

\[ E_0 + mA_s\Omega^2x \cos(\Omega t) \geq U_{\text{eff}} \]

Solving for phase

\[ \cos(\Omega t) \geq \frac{U_{\text{eff}} - E_0}{mA_s\Omega^2x} \]

If

\[ \beta = \frac{U_{\text{eff}} - E_0}{mA_s\Omega^2x}, \]

then the transmission probability becomes

\[ P = \frac{1}{\pi} \cos^{-1}(\beta) \]

for

\[ -1 \leq \beta \leq 1 \]

and

\[ P = 0 \]

if

\[ \beta > 1 \]

\[ \frac{\partial P}{\partial t} + c_s^2 \nabla \cdot (\rho \mathbf{v}) = 0 \]

with sound speed

\[ c_s = \sqrt{\frac{g \rho}{m}} \]

Reflected ripples from the far side of the barrier create a standing wave

\[ h(x, t) = A \cos(kx) \cos(\omega t) \]

The pressure gradient becomes

\[ \nabla P \propto -kA \sin(kx) \]

The particle experiences force

\[ \mathbf{F} = -\nabla P \]

Thus trajectory evolution follows

\[ m \frac{d\mathbf{v}}{dt} = -\nabla P \]

\[ \Psi(\mathbf{x}, t) = \sqrt{\rho(\mathbf{x}, t)} e^{i\theta(\mathbf{x}, t)} \]

with density

\[ \rho(\mathbf{x}, t) = |\Psi|^2 \]

and velocity field

\[ \mathbf{v} = \frac{\hbar}{m} \nabla \theta \]

The evolution of the system follows a nonlinear wave equation for the superfluid field

\[ i\hbar \frac{\partial \Psi}{\partial t} = \left[ -\frac{\hbar^2}{2m} \nabla^2 + V_{\text{sub}}(\mathbf{x}) + g |\Psi|^2 \right] \Psi \]

where

• \( V_{\text{sub}}(\mathbf{x}) \) is the periodic potential induced by the magnetic substrate
• \( g \) is the nonlinear interaction constant.

\[ \Psi(x, t) = \Psi_A(x, t) + \Psi_B(x, t) \]

with

\[ \Psi_A = A_A e^{i(k_A x - \omega_A t)} \]

\[ \Psi_B = A_B e^{i(k_B x - \omega_B t)} \]

The total density becomes

\[ \rho = |\Psi_A + \Psi_B|^2 \]

which expands to

\[ \rho = |\Psi_A|^2 + |\Psi_B|^2 + 2A_A A_B \cos[(k_A - k_B)x - (\omega_A - \omega_B)t] \]

The interference term

\[ \rho_{\text{int}} = 2A_A A_B \cos(\Delta k x - \Delta \omega t) \]

forms a standing wave pattern when

\[ \Delta \omega = 0 \]

Thus

\[ \rho_{\text{int}} = 2A_A A_B \cos(\Delta k x) \]

\[ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0 \]

Momentum equation

\[ m \left( \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla) \mathbf{v} \right) = -\nabla P \]

For small perturbations

\[ \rho = \rho_0 + \delta \rho \]

the pressure relation becomes

\[ P = c_s^2 \rho \]

with sound speed

\[ c_s = \sqrt{\frac{g \rho_0}{m}} \]

Combining the equations yields the wave equation

\[ \frac{\partial^2 (\delta \rho)}{\partial t^2} = c_s^2 \nabla^2 (\delta \rho) \]

Thus perturbations propagate through the superfluid at speed \( c_s \).

\[ \theta_A \]

\[ \theta_B \]

The interference density term depends on the phase difference

\[ \Delta \theta = \theta_A - \theta_B \]

Energy stored in the interference pattern

\[ E_{\text{int}} = -J \cos(\Delta \theta) \]

where \( J \) is the coupling strength determined by overlap of the pilot waves

\[ J \propto A_A A_B \]

The torque equation governing phase evolution becomes

\[ \frac{d(\Delta \theta)}{dt} = -\frac{\partial E_{\text{int}}}{\partial(\Delta \theta)} \]

which yields

\[ \frac{d(\Delta \theta)}{dt} = J \sin(\Delta \theta) \]

Stable equilibrium occurs at

\[ \Delta \theta = 0 \]

or

\[ \Delta \theta = \pi \]

Suppose particle \( A \) undergoes a phase shift

\[ \theta_A \to \theta_A + \delta\theta \]

The interference energy becomes

\[ E_{\text{int}} = -J \cos(\Delta\theta + \delta\theta) \]

The system evolves to restore equilibrium

\[ \Delta\theta' = \text{constant} \]

Therefore

\[ \theta_B \to \theta_B + \delta\theta \]

to maintain

\[ \Delta\theta' = \theta_A - \theta_B \]

\[ V_{\text{sub}}(\mathbf{x}) = V_0 \left[ \cos\left(\frac{2\pi x}{a}\right) + \cos\left(\frac{2\pi y}{a}\right) + \cos\left(\frac{2\pi z}{a}\right) \right] \]

Allowed vibrational modes satisfy the Bloch condition

\[ \Psi(\mathbf{x} + a\mathbf{n}) = e^{i\mathbf{k}\cdot\mathbf{a}\mathbf{n}} \Psi(\mathbf{x}) \]

where \( a \) is the lattice spacing.

The dispersion relation becomes

\[ \omega^2 = c_s^2 k^2 + \omega_L^2 \]

with lattice frequency

\[ \omega_L \propto \sqrt{\frac{V_0}{m}} \]

\[ h(x, t) = A \cos(kx) \cos(\omega t) \]

The pressure gradient from this standing wave

\[ \nabla P = -kA c_s^2 \sin(kx) \]

exerts forces on the two excitations.

Force on particle

\[ \mathbf{F} = -\nabla P \]

Thus

\[ m \frac{d\mathbf{v}}{dt} = -\nabla P \]

\[ J > k_B T_{\text{env}} \]

and when the phase synchronization time

\[ \tau_{\text{sync}} \sim \frac{1}{J} \]

is shorter than the propagation delay

\[ \tau_{\text{prop}} = \frac{L}{c_s} \]

where \( L \) is the separation distance.

Thus the condition for persistent entanglement becomes

\[ \frac{1}{J} < \frac{L}{c_s} \]

Let the universe be represented as a coherent superfluid membrane with density \(\rho\) and phase field

\[ \Psi(\mathbf{x}, t) = \sqrt{\rho(\mathbf{x}, t)} e^{i\theta(\mathbf{x}, t)} \]

The velocity field is

\[ \mathbf{v} = \frac{\hbar}{m} \nabla \theta \]

The system evolves according to the nonlinear wave equation

\[ i \hbar \frac{\partial \Psi}{\partial t} = \left[ -\frac{\hbar^2}{2m} \nabla^2 + V_{\text{sub}}(\mathbf{x}) + g |\Psi|^2 \right] \Psi \]

where the magnetic substrate imposes the periodic potential

\[ V_{\text{sub}}(\mathbf{x}) = V_0 \left[ \cos\left(\frac{2\pi x}{a}\right) + \cos\left(\frac{2\pi y}{a}\right) + \cos\left(\frac{2\pi z}{a}\right) \right] \]

with lattice spacing \(a\).

\[ \Psi = \sqrt{\rho_0 + \delta \rho} e^{i(\theta_0 + \delta \theta)} \]

Linearization gives density and phase perturbations governed by

\[ \frac{\partial^2 (\delta \rho)}{\partial t^2} = c_s^2 \nabla^2 (\delta \rho) \]

where

\[ c_s = \sqrt{\frac{g \rho_0}{m}} \]

The electromagnetic-like excitation corresponds to transverse phase oscillations of the membrane.

Let the phase disturbance be

\[ \delta \theta = A e^{i(kx - \omega t)} \]

Dispersion relation in the lattice medium becomes

\[ \omega^2 = c_s^2 k^2 + \omega_L^2 \]

where

\[ \omega_L = \sqrt{\frac{V_0}{m}} \]

\[ E_{12} = \rho \int \mathbf{v}_1 \cdot \mathbf{v}_2 dV \]

Substituting

\[ v_1 = \frac{\kappa}{2\pi r_1}, \quad v_2 = \frac{\kappa}{2\pi r_2} \]

The interaction potential becomes

\[ U(r) = \frac{\rho \kappa^2}{4\pi} \frac{1}{r} \]

which has the same spatial dependence as the Coulomb potential.

Thus the effective charge interaction constant becomes

\[ k_e = \frac{\rho \kappa^2}{4\pi} \]

Thus

\[ c = c_s \]

\[ c = \sqrt{\frac{g \rho_0}{m}} \]

\[ \kappa = \frac{h}{m} \]

thus

\[ h = m\kappa \]

\[ \alpha = \frac{k_e e^2}{\hbar c} \]

In the superfluid–substrate model the effective coupling arises from vortex circulation and medium density.

Substituting

\[ k_e = \frac{\rho \kappa^2}{4\pi} \]

and

\[ \hbar = \frac{h}{2\pi} = \frac{m \kappa}{2\pi} \]

gives

\[ \alpha = \frac{\rho \kappa^2 e^2}{4\pi} \frac{2\pi}{m \kappa c} \]

which simplifies to

\[ \alpha = \frac{\rho \kappa e^2}{2mc} \]

Substituting

\[ c = \sqrt{\frac{g \rho_0}{m}} \]

yields

\[ \alpha = \frac{\rho \kappa e^2}{2m} \left( \frac{m}{g \rho_0} \right)^{1/2} \]

\[ k_i = \frac{2\pi}{a}n_i \]

Resonant electromagnetic modes satisfy

\[ \omega = c_s k \]

The coupling constant \( e \) therefore emerges from the interaction between vortex circulation and the lattice harmonic modes

\[ e \propto \frac{\kappa}{a} \]

Substituting into the expression for \( \alpha \)

\[ \alpha \propto \frac{\rho \kappa^3}{2ma^2c} \]

\[ \psi_{surface}(z) = \psi_0 e^{-z/\delta_{skin}} \]

where

\(\psi_0\) is the amplitude at the interface and \(\delta_{skin}\) is the penetration depth.

For a bound excitation in a potential barrier the penetration depth is

\[ \delta_{skin} = \frac{\hbar}{\sqrt{2m \epsilon_{surface}}} \]

where

\(m\) is the effective mass associated with the excitation and \(\epsilon_{surface}\) is the surface energy density.

\[ \epsilon_{sub} \]

The interface energy is assumed to be smaller:

\[ \epsilon_{surface} \approx 10^{-3} \epsilon_{sub} \]

\[ \rho = \text{constant} \]

indicating an incompressible medium.

The bulk also contains a cubic magnetic lattice with coherence scale

\[ a_{\text{lattice}} \approx 10^{-3} \, \text{m} \]

Vorticity in the bulk is represented by

\[ \omega_{\text{bulk}} = \nabla \times \mathbf{v} \]

and typically satisfies

\[ |\omega_{\text{bulk}}| \gg |\omega_{\text{surface}}| \]

\[ V(z) = \epsilon_{surface} \left( 1 - e^{-z/\lambda_c} \right) + V_{sub} e^{-z/\delta_{bulk}} \]

where

\(\lambda_c\) is a characteristic interaction scale
\(\delta_{bulk}\) is the bulk attenuation length.

The potential minimum occurs where

\[ \frac{dV}{dz} = 0 \]

which produces an equilibrium position

\[ z_{surface} \sim 10^{-12} \, \text{m} \]

\[ \mathbf{J}(\mathbf{r}, t) \] creates oscillating electromagnetic fields governed by

\[ \nabla \times \mathbf{B} = \mu_0 \mathbf{J} \]

The resulting electromagnetic wave behaves as a surface mode analogous to a capillary wave.

The dispersion relation for surface waves in a fluid with surface tension is

\[ \omega^2 = gk + \frac{\sigma}{\rho} k^3 \]

where

\( g \) is gravitational acceleration
\( \sigma \) is surface tension
\( \rho \) is fluid density
\( k \) is the wave number.

\[ \omega^2 \approx \frac{\sigma}{\rho} k^3 \]

The phase velocity is

\[ v_p = \frac{\omega}{k} \]

Substituting the dispersion relation gives

\[ v_p = \sqrt{\frac{\sigma}{\rho} k} \]

In electromagnetic propagation the effective propagation speed can be written

\[ c = \frac{1}{\sqrt{\varepsilon_0 \mu_0}} \]

If electromagnetic waves correspond to tension-driven surface excitations, the effective velocity may be expressed as

\[ c \sim \sqrt{\frac{\sigma}{\rho \varepsilon_0 \mu_0}} \]

\[ \delta_\gamma \approx \frac{\lambda}{2\pi} \]

where \(\lambda\) is wavelength.

• vortex tangle density \( L(r) \)
• lattice spacing \( a_{\text{lattice}} \)
• global rotational vorticity.

Any physical measurement relies on a probe with characteristic wavelength \( \lambda_{\text{probe}} \).

Spatial resolution is limited by the relation

\[ \Delta x \gtrsim \frac{\lambda_{\text{probe}}}{2} \]

If the characteristic structural scale of the substrate is

\[ a_{\text{lattice}} \]

then direct resolution requires

\[ \lambda_{\text{probe}} \leq a_{\text{lattice}} \]

If

\[ \lambda_{\text{probe}} > a_{\text{lattice}} \]

\[ \Delta x \Delta p \geq \frac{\hbar}{2} \]

If a probe attempts to resolve structures on scale \( a_{lattice} \), the required momentum uncertainty becomes

\[ \Delta p \geq \frac{\hbar}{2a_{lattice}} \]

This introduces large momentum fluctuations that disrupt coherent measurement of fine substrate structures.

The effective measurement volume becomes

\[ V_{coh} \sim (\Delta x)^3 \]

\[a_{\text{lattice}} < r_p\] where \( r_p \) is the proton radius, then scattering experiments cannot resolve individual lattice nodes.

\[ \Omega_{bulk} = \nabla \times \mathbf{v}_{bulk} \]

Surface excitations interact with an averaged velocity field, producing

\[ \langle \Omega_{bulk} \rangle \approx 0 \]

\[ S_i = |S_i|e^{i\phi_i} \]

The phase \(\phi_i\) determines local substrate orientation.

Surface measurements detect only averaged values

\[ \langle S_i \rangle \]

making direct measurement of the individual phases \(\phi_i\) impossible without penetrating the bulk.

\[ Lk = \frac{1}{4\pi} \iiint \frac{(\mathbf{r}_1 - \mathbf{r}_2) \cdot (\mathbf{dr}_1 \times \mathbf{dr}_2)}{|\mathbf{r}_1 - \mathbf{r}_2|^3} \]

For a large ensemble of tangled vortices,

\[ \langle Lk \rangle \to 0 \]

\[ \rho_{surface} \]

and the observable operator be \( O \).

The expectation value is

\[ \langle O_{surface} \rangle = \text{Tr}(\rho_{surface} O) \]

If the substrate potential is \( V_{sub} \), the surface wavefunction satisfies

\[ \langle O_{surface} \rangle = \int \psi_{surface}^* V_{sub} \psi_{surface} \, dV \]

Because surface modes are orthogonal to bulk modes, contributions involving bulk vorticity gradients vanish in the surface average:

\[ \langle \nabla \times \omega_{bulk} \rangle = 0 \]

Let the angular velocity of galaxy \( i \) be \( \Omega_i \).

The average cosmic rotation becomes

\[ \Omega_{\text{universe}} = \lim_{V \to \infty} \frac{\sum_i \Omega_i}{V} \]

where the sum extends over galaxies within volume \( V \).

Let the characteristic scales be \( L_1, L_2 \). The lattice spacing may be inferred from harmonic relations.

For example,

\[ a_{\text{lattice}} \approx \frac{L_2 - L_1}{N} \]

where \( N \) is the number of harmonic intervals between the scales.

\[ L = \frac{\text{total vortex length}}{V} \]

If reconnections increase line density,

\[ \frac{dL}{dt} > 0 \]

entropy increases according to

\[ \frac{dS}{dt} \propto \frac{dL}{dt} \]

If the medium has density \( \rho \) and characteristic propagation velocity \( c \),

the effective tension scale becomes

\[ \sigma_{sub} \sim \rho c^2 \]

For depth \( z \),

\[ \psi(z) = \psi_0 e^{-z/a_{\text{lattice}}} \]

If a probe interacts with \( N \) lattice cells simultaneously,

\[ N \sim 10^{23} \]

random phase contributions cause the average signal to vanish:

\[ \langle S_i \rangle \approx 0 \]