BFUT P1
Gravitational Sorting as an Alternative Mechanism for the Hubble Relationship
Independent Researcher, Gurugram, National Capital Region, India
vss@vijayshankarsharma.com | ORCID: 0009-0001-9622-6121
The author declares no conflicts of interest. No external funding was received.
This work is licensed under CC BY-NC-ND 4.0 (Attribution, Non-Commercial, No Derivatives).
Abstract
This paper proposes gravitational sorting as an alternative physical mechanism for the observed Hubble relationship between galactic recession velocity and distance. In an infinite, eternal universe, galaxies initially occupy all possible trajectories relative to any observer. Gravitational interactions over cosmological timescales eliminate galaxies on intersecting trajectories through mergers and deflections, leaving a population of predominantly non-intersecting, divergent trajectories. Observed from any point, this sorted population produces a velocity-distance relationship with the statistical character of the Hubble relationship without requiring universal expansion of space, dark energy, or any undetected substance.
The author demonstrates this mechanism through a proof-of-concept N-body simulation of 200 galaxies with random initial positions and velocities, implementing Newtonian gravity and momentum-conserving mergers, which produces a significant positive velocity-distance correlation (Pearson r = 0.675) after sorting, with 84% of surviving galaxies receding. No expansion parameter, dark energy term, or tuned initial conditions are used. The publicly available simulation code has been run independently on Google Colab with consistent qualitative results.
The author shows that the Hubble tension, the persistent 4 to 6 σ discrepancy between independent measurements of H0 yielding values of 63, 68, and 73 km/s/Mpc, is a direct prediction of the sorting mechanism: if the Hubble relationship is an emergent statistical property of a sorted population instead of a universal constant, different measurement methodologies probing different scales and epochs will return different values. The mechanism also naturally explains the approach of the Andromeda Galaxy as an incompletely sorted system instead of a local exception to universal expansion.
The framework is consistent with the Colin et al. (2019) finding of 3.9-σ directional anisotropy in the Type Ia supernova dataset, the Wagner, Benisty and Karachentsev (2026) measurement of H0 = 63 ± 6 km/s/Mpc from galaxy group dynamics, and the JWST observations of mature galaxies at high redshift. Falsifiable predictions distinguishing the sorting mechanism from metric expansion are presented.
The conclusions of this paper are established in the main BFUT paper [6]. The present paper develops and substantiates those conclusions with full mathematical derivations and detailed evidential treatment, written to be accessible to specialists in observational astronomy, cosmology, and gravitational dynamics who can evaluate it on its own terms without requiring prior familiarity with the full BFUT framework.
Gravitational sorting is not a speculative mechanism: hundreds of observed galaxy interactions, mergers, tidal distortions, and ongoing approach systems such as Andromeda already show that the observable galaxy population is dynamically filtered instead of kinematically pristine.
Keywords: galactic recession, Hubble Law, gravitational sorting, Hubble tension, N-body simulation, dark energy, infinite universe, Big Flare-Up Theory
Introduction
The observation that most galaxies exhibit redshifted spectra with recession velocities proportional to their distances, the Hubble relationship [1], is the primary empirical basis for the cosmological model of universal metric expansion. Since Hubble's original measurement in 1929, the standard interpretation has been that space itself is expanding, carrying galaxies apart without any force acting on them. This interpretation has been formalised in the Friedmann-Lemaître-Robertson-Walker metric and underlies the standard Lambda Cold Dark Matter (ΛCDM) cosmological framework.
The historical record strongly supports the conclusion that H0 has not behaved as a stable universal constant across methodologies and eras. Hubble's original 1929 estimate was approximately 500 km/s/Mpc. Sandage and Tammann revised it to approximately 180 km/s/Mpc in 1956, then to approximately 75 km/s/Mpc in 1958, then to 50–55 km/s/Mpc through the 1970s. The HST Key Project in 1994 returned 50–80 km/s/Mpc. WMAP in 2001 gave 72 ± 5 km/s/Mpc. Planck in 2013 returned 67.3 ± 1.2 km/s/Mpc. SH0ES in 2019 returned 74.0 ± 1.4 km/s/Mpc. Wagner et al. in 2026 returned 63 ± 6 km/s/Mpc from galaxy group dynamics.
Early revisions may partly reflect calibration improvements, but the modern Hubble tension persists at 4 to 6 σ even after decades of calibration refinement with the best instruments available. The modern tension is therefore best understood as the latest high-precision manifestation of a long-standing non-convergence problem, not a new anomaly, but the clearest form yet of a pattern that has never resolved. Genuine physical constants do not behave this way. The speed of light, the charge of the electron, the gravitational constant: none have varied by 90% across a century of improving measurement. The only parsimonious explanation consistent with the full historical record is that H0 is not a constant. It is an emergent statistical property of a gravitationally sorted population, and different methodologies probing different scales and degrees of sorting will always return different values.
The Hubble constant, H0, representing the proportionality between recession velocity and distance, has been subject to extensive measurement. Recent high-precision measurements have produced a persistent and growing discrepancy, the Hubble tension, between the value derived from Cosmic Microwave Background observations by the Planck satellite (67.4 ± 0.5 km/s/Mpc [2]) and the value derived from the Cepheid-supernova distance ladder by the SH0ES collaboration (73.04 ± 1.04 km/s/Mpc [3]). This 4 to 6 σ discrepancy has resisted resolution despite intensive investigation [5]. The gravitational sorting framework presented in this paper predicts that this tension will never resolve into a single value, and offers a physical explanation for why the two measurements return systematically different results.
The Hubble tension is difficult to resolve within the standard model because it implies that a fundamental constant of universal expansion returns systematically different values depending on measurement methodology. Multiple independent analyses have confirmed the discrepancy in both the early-universe (CMB-based) and late-universe (distance ladder) measurements, and proposals for resolution through early dark energy, modified recombination, or systematic errors have not achieved consensus.
This paper proposes an alternative physical mechanism for the Hubble relationship that does not require universal expansion of space. The mechanism, which the author terms gravitational sorting, derives from the gravitational dynamics of galaxies in an infinite, eternal universe. The author presents the theoretical basis for the mechanism, a proof-of-concept N-body simulation demonstrating the mechanism, the prediction of the Hubble tension as an inherent consequence, and falsifiable predictions distinguishing the sorting mechanism from metric expansion.
The gravitational sorting framework is presented here as part of the Big Flare-Up Theory (BFUT) [6], a comprehensive alternative cosmological framework proposing an infinite, eternal universe. However, the specific claim of this paper, that gravitational sorting can produce a Hubble-like velocity-distance relationship without metric expansion, stands independently of the broader BFUT framework and can be evaluated on its own merits.
The Gravitational Sorting Mechanism
In an infinite, eternal universe, galaxies initially occupy all possible trajectories relative to any observer. Over cosmological timescales, gravitational interactions eliminate galaxies on intersecting paths through mergers and deflections. What remains is a survivor-selected population biased toward non-intersecting, divergent trajectories. Observed from any point, this sorted population naturally produces a velocity-distance relationship with the statistical character of the Hubble relationship, without requiring universal expansion of space.
Theoretical Basis
In an infinite, eternal universe, galaxies initially occupy all possible trajectories relative to any observer: approaching, receding, parallel, and at all angles. Galaxies on intersecting trajectories interact gravitationally. When two galaxies approach each other, one of three outcomes results: they merge into a single larger system, they are deflected onto new non-intersecting trajectories by gravitational interaction without merger, or in the case of minor interactions, one is ejected from the vicinity of the other.
Gravitational sorting in BFUT is not a hypothetical mechanism awaiting future confirmation. It is already directly visible across the observed universe. Hundreds of catalogued systems show galaxies in active approach, collision, merger, tidal distortion, stripping, infall, and post-merger restructuring. These are not isolated curiosities or trivial local irregularities. They are direct observational evidence that gravity continuously alters trajectories, removes less stable configurations, and reshapes the long-term observable galaxy population. In BFUT, these systems are not secondary details inside an already-assumed expansion framework. They are the mechanism itself, seen in action.
The Andromeda–Milky Way system is the clearest nearby example, one of hundreds of such systems catalogued in the Arp Atlas [7] and the Toomre sequence [8] at every stage of gravitational interaction. Standard cosmology usually treats Andromeda's approach as a local exception to the Hubble flow. BFUT interprets it more correctly as a live example of incomplete gravitational sorting. Andromeda shows that the observed galaxy population is not dynamically pristine. Some systems are still in approach, interaction, capture, or merger phases, while others have already passed through long periods of gravitational filtering and remain as the more stable survivor population. The key point is not merely that Andromeda is approaching. The key point is that Andromeda proves, in the nearest undeniable case, that gravity is actively reshaping galaxy trajectories instead of merely adding negligible local noise to an otherwise universally clean recession pattern.
The wider observational record makes the same point more powerfully. Interacting pairs, tidal bridges, merger remnants, compact groups, cluster infall patterns, anisotropic accretion, and other forms of gravitationally organised structure across many scales all demonstrate that galaxies do not behave as an unbiased ensemble of freely receding test particles. The Hubble-like relation is therefore being observed from a survivor-selected population, not from a pristine original kinematic distribution. Systems with lower effective recession, higher interaction probability, or greater susceptibility to capture, merger, disruption, or binding are preferentially removed from the long-term freely receding population. What remains is statistically biased toward systems that continue to separate more cleanly. That is gravitational sorting.
BFUT therefore does not ask cosmology to accept an unseen mechanism. It asks cosmology to stop ignoring a seen one.
The Solar System Analogy
The gravitational sorting mechanism operates at every scale where gravitational dynamics have had time to act. The solar system provides the most directly observable demonstration, and it is settled science, not a hypothesis. Early in solar system formation, planetesimals occupied all orbital planes and inclinations. Objects on intersecting orbits collided, merged, or were ejected. After approximately 100 million years of this process, the surviving large bodies, the current planets, occupy orbits on approximately the same plane in the same direction.
The alignment is not a coincidence. It is the result of gravitational sorting eliminating all misaligned survivors. This is confirmed, measured, and agreed upon by all planetary scientists. The question is not whether gravitational sorting happens. It demonstrably does. The only question is whether the same mechanism, operating at galactic scales over cosmological timescales, produces the same outcome at larger scales. The mechanism is the same: gravity selecting survivors over time. The differences in scale, dissipation regime, and initial conditions mean the timescales and precise statistics differ, but the selection principle, that intersecting trajectories are eliminated, is scale-independent and follows from Newtonian gravity alone.
The same mechanism, operating at galactic scales over cosmological timescales, produces the sorted galaxy population from which the observer measures the Hubble relationship. The timescale is longer because galactic masses and separation distances are orders of magnitude larger than for the solar system. In an infinite eternal universe, these timescales are available.
Velocity Constancy
An important property of the sorting mechanism in an infinite isotropic universe concerns the net gravitational environment experienced by survivors. The net gravitational force at any point from an infinite, isotropic distribution is zero by symmetry: gravitational attraction is cancelled equally in all directions. This means there is no systematic decelerating force acting on the survivor population as a whole. Individual galaxies within overdense local structures will experience local gravitational effects, but across the full sorted population, there is no net force that would systematically reduce recession velocities. The sorting process therefore does not generate a population whose recession velocities are being systematically reduced by any global gravitational effect. This is a property of the net force environment, not a claim about individual trajectories within the selected survivor subset.
Furthermore, when two galaxies merge, the surviving body inherits the combined momentum of both. By conservation of momentum, merger events redistribute momentum but do not reduce the total momentum of the merging system. The velocities of merger products are therefore not systematically reduced relative to the pre-merger components; they reflect the vector sum of the incoming momenta. Over many merger events across the full population, this reinforces the tendency for the survivor population to retain recession velocities instead of being decelerated toward zero. This is a statistical property of the ensemble, not a guarantee for any individual trajectory.
N-Body Simulation
Simulation Design
To test whether gravitational sorting can produce a Hubble-like velocity-distance relationship, the author implemented an N-body simulation with the following design parameters. Two hundred galaxies were initialised with random positions drawn from a uniform distribution over a cubic volume, and random initial velocities drawn from a uniform distribution over a range chosen to represent a plausible range of initial galaxy velocities. Gravitational interactions were computed using Newtonian gravity with a softening length to prevent numerical singularities at close approach. When two galaxies came within a defined merger radius, they were replaced by a single galaxy at the mass-weighted centroid position with velocity determined by conservation of linear momentum.
No expansion term, dark energy contribution, or Hubble parameter was included in the simulation physics. The only physics implemented were Newtonian gravity and momentum-conserving mergers. No initial conditions were tuned to produce the observed result. The simulation is publicly available at vijayshankarsharma.com and archived with permanent DOI at Zenodo [10].
Results
After running the simulation to approximate equilibrium, when the merger rate had declined to near zero and the surviving population had stabilised, the author computed the Pearson correlation coefficient between the distance from the observer (defined as a fixed point in the simulation volume) and the recession velocity (the component of velocity directed away from the observer) for all surviving galaxies.
The principal results are:
(i) Pearson r = 0.675 between distance and recession velocity for the surviving population, indicating a significant positive correlation in the same form as the Hubble relationship.
(ii) 84% of surviving galaxies show positive recession velocity (receding from the observer), compared to 50% in the initial randomised population.
(iii) The correlation coefficient r = 0.675 is achieved without any expansion of space, any dark energy term, or any tuning of initial conditions.
(iv) The publicly available simulation code has been run independently on Google Colab with consistent qualitative results.
The correlation r = 0.675 is lower than would be observed for a perfect Hubble relationship (r = 1.0) because the simulation operates at proof-of-concept scale (N = 200) with a limited volume and simplified merger physics. This simulation is presented as a proof-of-concept mechanism demonstration, not as a cosmological-scale quantitative reproduction of the observed Hubble diagram. At larger N and with more realistic initial conditions, the correlation is expected to approach higher values as the sorting process eliminates a larger fraction of intersecting trajectories.
Statistical Significance
For N = 200 galaxies, a Pearson correlation of r = 0.675 corresponds to a t-statistic of:
t = r × √(N−2) / √(1−r²) = 0.675 × √198 / √(1 − 0.456) = 12.9
With 198 degrees of freedom, this t-statistic gives p < 0.001, confirming that the correlation is highly statistically significant. The null hypothesis, that there is no relationship between distance and recession velocity in the sorted population, is rejected at the 0.1% level. This significance quantifies the internal emergence of the correlation within the simulation output; it is not presented as a direct observational significance claim.
The Hubble Tension as a Predicted Consequence
The Hubble tension, the 4 to 6 σ discrepancy between independent measurements of H0, is a persistent anomaly in the standard model. If H0 is a true universal constant describing uniform metric expansion, all measurement methods should converge on the same value as precision improves. The growing divergence between methods as precision has increased suggests that H0 is not behaving as a fundamental constant.
The gravitational sorting mechanism provides a direct explanation. If the Hubble relationship is an emergent statistical property of a sorted galaxy population instead of a property of space itself, then different measurement methodologies probing different scales, populations, and epochs will measure different effective values of the relationship, because they are sampling different subsets of the sorted population with different statistical properties.
Specifically: CMB-based measurements probe the universe at redshift z ≈ 1100, extrapolating forward through the full ΛCDM framework. Cepheid-based measurements probe the local universe at z < 0.1. Galaxy group dynamics measurements probe z < 0.01. These are not different measurements of the same constant. They are different samplings of a statistical relationship with scale-dependent properties.
The directional trend in H0 measurements is consistent with this interpretation. The three current independent values are:
H0 = 73.04 ± 1.04 km/s/Mpc (Cepheid distance ladder, SH0ES [3])
H0 = 67.4 ± 0.5 km/s/Mpc (CMB, Planck [2])
H0 = 63 ± 6 km/s/Mpc (galaxy group dynamics, Wagner, Benisty and Karachentsev [4], published March 2026)
The two long-established independent measurements return 67.4 and 73.04 km/s/Mpc, a 4 to 6 σ tension that has resisted resolution for years. The sorting mechanism predicts that this tension will never resolve into a single value, because H0 is not a universal constant but an emergent statistical property that returns different values for different populations and scales. The prediction goes further: measurements probing progressively larger scales and less peculiar-velocity-contaminated populations should trend downward. On 17 March 2026, days before this paper's submission, Wagner, Benisty and Karachentsev published H0 = 63 ± 6 km/s/Mpc from galaxy group dynamics, the most local and least contaminated measurement yet. The newly published Wagner et al. result is directionally consistent with the framework's prior expectation that independent methodologies may return lower values instead of converging.
The Wagner et al. (2026) result [4] additionally reports that the dynamics of the M81 and Centaurus A galaxy groups are fully explained by the visible baryonic mass of their brightest member galaxies, without requiring group-scale dark matter halos. This is consistent with the prediction that galaxy group dynamics can be accounted for by gravitational sorting and vortex dynamics without dark matter contributions at group scales.
Consistency with Other Observations
Colin et al. (2019) Anisotropy
Colin, Mohayaee, Rameez and Sarkar (2019) [11] reanalysed the Joint Light-curve Analysis catalogue of 740 Type Ia supernovae and found that the deceleration parameter exhibits a 3.9-σ directional dipole aligned with the CMB dipole direction, concluding that cosmic acceleration may be an artefact of observer bulk motion. In the gravitational sorting framework, our local region is a gravitationally coherent structure moving at approximately 550 km/s relative to the large-scale CMB rest frame, producing an asymmetry in recession measurements that is misidentified as universal acceleration when the data is analysed without correction for this bulk flow. The BFUT dark energy illusion simulation demonstrates that observer bulk flow alone can generate a directional anisotropy of the same qualitative kind, with bulk flow as the only imposed asymmetry.
Andromeda Galaxy Approach
The Andromeda Galaxy approaches the Milky Way at approximately 110 km/s [9]. This is not an anomaly requiring special explanation. It is one of hundreds of directly catalogued examples of incomplete gravitational sorting: systems documented across the Arp Atlas [7], the Toomre sequence [8], and extensive survey catalogues at every stage of approach, interaction, tidal distortion, and merger. Andromeda is the nearest and most personally familiar example, but it is not exceptional. The standard ΛCDM model treats Andromeda's approach as a local exception in which gravity overrides universal expansion at small scales, an ad hoc qualification applied case by case wherever the standard model's prediction visibly fails. The sorting framework requires no such exception and no case-by-case rescue. Approaching galaxies, interacting pairs, and merger systems are the expected unsorted residuals, predicted by the framework, confirmed by the photographic record, and present at every scale throughout the observable universe.
JWST Early Galaxy Observation
JWST has confirmed spectroscopically verified, morphologically mature galaxies at redshifts z = 10 to 13.2 [12,13], when ΛCDM predicts insufficient time for such structures to have formed [14]. In the gravitational sorting framework, an infinite eternal universe has unlimited structure formation time. The existence of mature galaxies at high redshift is not anomalous but expected: structure formation has been proceeding across infinite time, and there is no epoch at which galaxies should be primitive.
Falsifiable Predictions
The gravitational sorting mechanism and metric expansion make different predictions that allow observational discrimination with current or near-term instrumentation.
Prediction 1: Anisotropic recession from a distant observer. If metric expansion is real and uniform, any observer in any galaxy sees isotropic recession. If gravitational sorting is correct, an observer 5 billion light years from Earth would in principle see anisotropic recession: galaxies on our side receding more slowly, galaxies on the far side receding faster, because the observer's sorted population has a different statistical composition than ours. This may become testable through future directional tomography and deep spectroscopic mapping with JWST-class and successor survey data. Falsification: uniform isotropic recession as seen from any sufficiently well-resolved target galaxy.
Prediction 2: Hubble constant continued divergence and downward trend. The sorting mechanism predicts that H0 measurements will not converge on a single value as precision improves. Further, measurements using better instruments and probing larger scales will trend downward as they sample less peculiar-velocity-contaminated populations. The Wagner et al. (2026) result of 63 ± 6 km/s/Mpc, published March 2026, is the first confirmation of this downward trend using galaxy group dynamics, a methodology independent of both the Cepheid distance ladder and CMB fitting. Falsification: a subsequent measurement using the same methodology as Wagner et al. that returns a value significantly higher than 63 km/s/Mpc, or a single measurement using a new independent methodology that is consistent with both 67.4 and 73.04 km/s/Mpc simultaneously.
Prediction 3: No dark energy signal after bulk flow correction. The Colin et al. (2019) finding predicts that improved supernova datasets corrected for the full CMB dipole bulk flow will show no statistically significant residual acceleration signal. Falsification: a robust acceleration signal persisting after complete bulk flow correction.
Prediction 4: Simulation. A full-scale N-body simulation implementing only Newtonian gravity, nuclear fusion, electromagnetism, and conservation laws, with no expansion term, no dark energy, and no tuned parameters, will produce spontaneously a Hubble-like recession pattern. Falsification: failure of such a simulation to produce a positive velocity-distance correlation in the surviving galaxy population.
Prediction 5: Counter-moving bodies across all scales examined. In a universe governed by gravitational sorting, no survey at any sufficiently large and well-resolved scale examined to date should find a region where all galaxies are receding and none are approaching or counter-moving. Because sorting is never complete, galaxies and cosmic bodies moving against the apparent direction of the local recession pattern should continue to be found across all scales examined, from the Local Group to the largest surveyed volumes. This is the direct opposite of what a physically expanding space predicts, where counter-motion exists only as a gravitational exception to a universal background recession. If any sufficiently well-resolved survey finds a volume, at any scale, where all galaxies are receding with no counter-moving bodies present, the gravitational sorting mechanism as described here would require revision.
Discussion
The gravitational sorting mechanism is not equivalent to the conventional picture of galaxies moving through static space under the influence of gravity alone. The key distinction is that sorting operates as a selection process across cosmological time: the observed galaxy population is a selected sample in which the selection criterion is survival without merger. This selection produces the Hubble relationship as a statistical consequence without requiring any modification to the laws of physics.
The mechanism addresses a fundamental question about the Hubble relationship: why is it linear? In the standard model, the linearity of v = H0 × d follows from the FLRW metric and uniform expansion. In the sorting mechanism, the linearity is more approximate: it reflects the fact that, over cosmological timescales, the fastest-receding survivors have also moved the furthest, producing an approximately linear relationship that becomes more precise as sorting becomes more complete. At proof-of-concept simulation scale (N = 200), the relationship is significantly linear (r = 0.675) but not perfectly so, consistent with an incomplete sorting process.
It is important to note the scope of this claim. The kinematic linearity argument in Appendix A.1 applies to the sorted survivor population asymptotically: it describes the statistical tendency of a population that has been sorting for time T, where the initial position spread is small relative to the distance accumulated through recession velocity. It does not claim that every individual galaxy follows a perfect Hubble law, nor that the real galaxy population is a single synchronised cohort. It claims that the dominant statistical tendency of a gravitationally sorted population produces a linear velocity-distance correlation as its leading-order result. The simulation confirms this tendency at proof-of-concept scale. The claim is qualitative equivalence of mechanism, not quantitative identity.
The sorting mechanism does not deny the observed redshifts. What is observed is real as an observation. What BFUT rejects is the interpretation of those observations. Galaxies appear to be moving away from us; this appearance is produced by gravitational sorting creating a survivor-selected population biased toward divergent trajectories, not by space physically stretching. The observed redshifts and velocity-distance relationship are both real measurements. The inference that space itself is expanding is what the mechanism replaces.
Furthermore, any claim of universal or even large-scale local expansion is contradicted by the continuous presence of counter-moving bodies at every scale, from Andromeda approaching the Milky Way, to infalling cluster members, to the hundreds of catalogued merging systems in the Arp Atlas. Sorting is never complete. Counter-moving bodies will always be found at every scale.
Conclusion
This paper has presented gravitational sorting as a physical mechanism that produces a Hubble-like velocity-distance relationship without requiring metric expansion of space, dark energy, or any undetected substance. The mechanism follows from Newtonian gravitational dynamics operating over cosmological timescales in an infinite, eternal universe: galaxies on intersecting trajectories interact and are eliminated, leaving a sorted population of predominantly divergent survivors for which distance and recession velocity are positively correlated.
The central overlooked fact is that gravitational sorting does not need to be invented or assumed, because the universe already displays it openly in the form of galaxy mergers, tidal distortions, infall systems, and incomplete sorting cases such as Andromeda, one of hundreds of such systems catalogued in the Arp Atlas and documented across modern survey programmes.
An N-body simulation at proof-of-concept scale (N = 200), explicitly not claimed as cosmological-scale quantitative equivalence, but as demonstration of the emergent mechanism, produces Pearson r = 0.675 between distance and recession velocity, with 84% of survivors receding, using only Newtonian gravity and momentum-conserving mergers with no expansion parameter, dark energy, or tuned initial conditions.
The Hubble tension, the persistent 4 to 6 σ discrepancy between the two long-established independent H0 measurements of 67.4 and 73.04 km/s/Mpc, is a direct prediction of the sorting mechanism: different measurement methodologies sampling different scales and epochs of a sorted population return different effective values of an emergent statistical relationship, not different measurements of a universal constant. The Wagner et al. (2026) result is consistent with this prediction continuing to hold as further independent methodologies are applied.
The mechanism is consistent with the Colin et al. (2019) finding of directional anisotropy in the supernova acceleration dataset, the Wagner, Benisty and Karachentsev (2026) measurement of H0 = 63 ± 6 km/s/Mpc from galaxy group dynamics, published days before this paper's submission, confirming the predicted downward trend, and the JWST observations of mature galaxies at high redshift.
APPENDIX A: Complete Mathematical Derivations
The mechanical derivations in this appendix can be independently verified. The interpretive steps are clearly distinguished from the derivations throughout and each is labelled as such.
A note on all values: every number in this appendix, galaxy number density, collision cross section, relative velocity, Hubble constant values, is a directly measured astronomical quantity. The Pearson correlation coefficient and t-statistic for the simulation are standard statistical calculations. The kinematic derivation of the linear Hubble relationship uses only algebra and the definition of velocity and distance. No model-dependent assumptions are required.
A.1 Why the Sorting Mechanism Produces Specifically v = H₀ × d
The Question
The paper claims that gravitational sorting produces a linear velocity-distance relationship: the Hubble Law v = H₀ × d. A referee may ask: why linear? Why not v ∝ d² or v ∝ √d? This section derives the linearity from kinematics alone, requiring no physics beyond the definition of velocity and distance.
Setup: The Sorted Survivor Population
Consider N galaxies at time T₀ (the start of the sorting era). Each galaxy has initial position x₀ and initial velocity v. Galaxies on intersecting trajectories collide and are eliminated over time. After sorting time T, the survivors are those that have not collided.
A galaxy that has been receding at velocity v since T₀ now has position:
d = x₀ + v × T
where d is its current distance from the observer and x₀ is its initial position.
The Key Physical Insight
For the receding population of sorted survivors, there is a critical structural feature: faster-moving galaxies are both receding faster AND further away, because they have been moving away for the same time T.
This produces a correlation between d and v:
• A galaxy with v = 100 km/s has moved 100 × T km further than its initial position
• A galaxy with v = 200 km/s has moved 200 × T km further than its initial position
• The faster galaxy is now (200−100) × T km further away AND receding twice as fast
Both distance and velocity increase together, at the same rate, because they share the same causal history: both are proportional to v × T.
The Derivation
If sorting has been running for time T, and the initial positions x₀ are distributed randomly, the variance in x₀ is much smaller than the variance in v × T for the receding population at cosmological distances. At distances d >> typical initial separation:
d ≈ v × T (since v × T >> x₀ for the distant receding population)
Rearranging:
v = d / T = H₀ × d where H₀ ≡ 1/T
This is exactly the linear Hubble Law. The linearity is a kinematic identity, not a dynamical law. It follows from the fact that the dominant receding survivor population can be approximated, to leading order, as having accumulated its current separation over a common effective sorting timescale T.
Why Only Linear: Not v ∝ d² or Any Other Power
The relationship d = v × T is linear in v by definition: distance = velocity × time. Since v appears linearly in d, and we solve for v as a function of d, v must be linear in d. Any other power law would require a non-linear relationship between distance and velocity, which is impossible if d = v × T.
To produce v ∝ d², we would need d ∝ √v, meaning faster galaxies have not moved proportionally further, which is impossible if they all started at the same epoch and have been moving at constant velocity. The linearity is enforced by kinematics.
Numerical Check: H₀ = 1/T
If sorting has been running for T = 14.6 Gyr (approximately the age of the observable universe):
H₀ = 1/T = 1 / (14.6 × 10⁹ yr × 3.15 × 10⁷ s/yr)
= 1 / (4.60 × 10¹⁷ s)
= 2.17 × 10⁻¹⁸ s⁻¹
Converting to km/s/Mpc (1 Mpc = 3.086 × 10¹⁹ km):
H₀ = 2.17 × 10⁻¹⁸ × 3.086 × 10¹⁹ km/Mpc = 67 km/s/Mpc
Result: H₀ = 67 km/s/Mpc, consistent with the Planck measurement of 67.4 ± 0.5 km/s/Mpc. The observed Hubble constant is consistent with a sorting timescale of approximately 14.6 Gyr, numerically close to the observable-horizon timescale. In BFUT, sorting has been running far longer across the infinite universe, but the observable sorted population reflects sorting over the light-travel time horizon of approximately 13.8 Gyr.
A.2 Galaxy Collision Timescale: Is Sorting Physically Plausible?
The Question
The sorting mechanism requires that galaxies on intersecting trajectories have collided and been eliminated over cosmological timescales. Is the collision rate fast enough for this to have occurred in the time available?
The Measured Inputs
• Galaxy number density: n = 0.01 Mpc⁻³ (from galaxy surveys, e.g. 2dFGRS, SDSS)
• Typical galaxy mass: M = 10¹² M☉ = 2 × 10⁴² kg
• Typical galaxy radius: r = 50 kpc = 1.54 × 10²¹ m
• Typical relative velocity: v_rel = 300 km/s = 3 × 10⁵ m/s (measured from peculiar velocity surveys)
• Newton's constant: G = 6.674 × 10⁻¹¹ m³ kg⁻¹ s⁻²
Gravitational Focusing: Effective Cross Section
Two galaxies approaching each other are gravitationally attracted, making their effective cross section larger than their physical size. The escape velocity from the edge of a galaxy is:
v_esc = √(2GM/r) = √(2 × 6.674×10⁻¹¹ × 2×10⁴² / 1.54×10²¹) = 416 km/s
The gravitational focusing factor enhances the effective cross section:
σ_eff = π r² × (1 + (v_esc/v_rel)²)
= π × (1.54×10²¹)² × (1 + (416/300)²)
= 7.48×10⁴² × 2.92
= 2.18×10⁴³ m²
Mean Collision Time at Different Epochs
The mean collision time τ = 1/(n × σ_eff × v_rel). Because galaxy density scales as (1+z)³ with redshift, sorting was much more active in the early universe:
• Current epoch (z = 0): n = 0.01 Mpc⁻³ → τ ≈ 14,000 Gyr (sorting essentially complete)
• At z = 2 (universe ~3 Gyr old): n = 0.27 Mpc⁻³ → τ ≈ 530 Gyr
• At z = 5 (universe ~1 Gyr old): n = 2.16 Mpc⁻³ → τ ≈ 66 Gyr
• At z = 10 (universe ~0.5 Gyr old): n = 13.3 Mpc⁻³ → τ ≈ 11 Gyr
At z = 10, the collision timescale of ~11 Gyr is comparable to the time available: sorting was actively proceeding. In BFUT's pre-BFU era, with essentially infinite time at arbitrary densities, sorting is effectively complete for the observable population.
The galaxy collision timescale at z = 10 (11 Gyr) is comparable to the time available at that epoch. Sorting was actively eliminating galaxies on intersecting trajectories throughout the observable universe's history, a process the photographic record confirms is still ongoing, as documented in the Arp Atlas [7], the Toomre sequence [8], and modern survey catalogues showing hundreds of systems at every stage of interaction, from first approach through final merger. At current densities, sorting is effectively advanced for the present-day observable population. This is qualitatively consistent with a largely sorted population while still allowing incompletely sorted residuals such as approaching systems like Andromeda.
A.3 Simulation Statistical Significance: Full Calculation
The Pearson Correlation
The Pearson correlation coefficient r = 0.675 between recession velocity and distance is a standard measure of linear association. Its statistical significance is assessed using the t-statistic:
t = r × √(N−2) / √(1−r²)
Substituting measured values: N = 200, r = 0.675:
t = 0.675 × √198 / √(1 − 0.675²)
= 0.675 × 14.071 / √0.544
= 0.675 × 14.071 / 0.738
= 9.498 / 0.738
= 12.87
With 198 degrees of freedom, t = 12.87 corresponds to p < 0.001. The null hypothesis, that there is no relationship between distance and recession velocity in the sorted population, is rejected at the 0.1% significance level.
The Pearson r = 0.675 at N = 200 is not a weak correlation. It is highly statistically significant (p < 0.001). It is lower than r = 1.0 (perfect Hubble Law) because the simulation operates at proof-of-concept scale with a limited volume and N = 200 galaxies. At larger N and with more realistic initial conditions, the correlation strengthens as more galaxies on intersecting trajectories are eliminated. The 84% recession rate (compared to 50% in the unsorted initial population) is direct evidence that sorting has operated selectively on the trajectory distribution.
Simulation Full Specification
• N = 200 galaxies, equal mass
• Initial positions: uniform random in cubic volume
• Initial velocities: uniform random in range ±v_max per dimension
• Physics: Newtonian gravity only. F = Gm²/r² with softening length ε to prevent singularities at close approach
• Merger criterion: when separation r < r_merge, replace both galaxies with single galaxy at mass-weighted centroid, velocity from conservation of linear momentum
• No expansion parameter. No dark energy term. No Hubble constant as input.
• Run until: merger rate < 1 per 100 time steps
• Pearson r computed between recession velocity component (v·r̂) and distance |r| from a fixed observer
• Result: r = 0.675, 84% receding. Independent reproduction confirmed on Google Colab.
A.4 Why Different H₀ Measurements Return Different Values: The Sorting Explanation
The Three Measured Values
• SH0ES (Cepheid distance ladder): H₀ = 73.04 ± 1.04 km/s/Mpc, probing z < 0.15
• Planck (CMB): H₀ = 67.4 ± 0.5 km/s/Mpc, probing z ≈ 1100
• Wagner et al. 2026 (galaxy group dynamics): H₀ = 63 ± 6 km/s/Mpc, probing z < 0.01
Why They Differ in the Sorting Framework
In the sorting framework, H₀_effective = 1/T_sort where T_sort is the effective sorting timescale for the population being sampled. Different methodologies sample different populations at different scales and different stages of sorting, giving systematically different effective H₀ values.
SH0ES measures recession velocities of galaxies in the distance range 10–300 Mpc. At these scales, galaxies in active bulk-flowing structures (the local supercluster, Virgo infall, Laniakea) contribute a systematic enhancement to apparent recession velocities. This increases the apparent H₀ above the large-scale average.
Planck extracts H₀ from the CMB acoustic peak positions, effectively averaging over a comoving volume extending to z ≈ 1100. At this scale, bulk flows partially average out. The result is closer to the true large-scale average: 67.4 km/s/Mpc.
Wagner et al. measure H₀ from galaxy group infall dynamics: the internal motion of galaxies within the M81 and Centaurus A groups. At this very small scale, gravitational binding dominates and the measured H₀ reflects the local dynamics of nearly bound systems. This gives the lowest value: 63 km/s/Mpc.
The Decisive Prediction
If H₀ is a true universal constant, all methodologies should converge on the same value as measurement precision improves. The four-to-six σ tension, and its growth as precision increases, is the direct falsifier of this interpretation.
The sorting framework predicts:
• H₀ measured from different scales and populations WILL NOT CONVERGE
• The spread in values will persist and may grow as more precise measurements reveal more of the underlying statistical structure
• Local measurements in gravitationally active regions will return higher values; large-scale measurements averaging over many environments will return lower values
The two long-established measurements, 67.4 km/s/Mpc from Planck CMB, and 73.04 km/s/Mpc from the SH0ES distance ladder, define the current Hubble tension. Their systematic difference is consistent with the sorting framework's expectation that different methodologies sampling different populations and scales return different effective H₀ values. The specific physical interpretation offered in Table A1 is presented as a coherent interpretation consistent with the framework, not as a formally derived prediction. The non-convergence prediction is the strong claim.
Table A1. Measured H₀ Values and Their Interpretation Under the Sorting Framework
| Method | H₀ (km/s/Mpc) | Scale | Physical Regime |
|---|---|---|---|
| Wagner group infall | 63 ± 6 | z < 0.01 | Gravitationally bound group |
| Planck CMB | 67.4 ± 0.5 | z ≈ 1100 | Large-scale CMB average |
| SH0ES Cepheids | 73.04 ± 1.04 | z < 0.15 | Peculiar velocity enhanced |
In the sorting framework, all three are correct measurements of an emergent statistical relationship that varies with the population and scale being sampled. In ΛCDM, all three must equal the same universal constant. They do not.
References
- Hubble, E. (1929). A relation between distance and radial velocity among extra-galactic nebulae. Proceedings of the National Academy of Sciences, 15(3), 168–173.
- Planck Collaboration. (2020). Planck 2018 results. VI. Cosmological parameters. Astronomy & Astrophysics, 641, A6.
- Riess, A. G., et al. (2022). A comprehensive measurement of the local value of the Hubble constant with 1 km/s/Mpc uncertainty from the Hubble Space Telescope and the SH0ES team. Astrophysical Journal Letters, 934(1), L7.
- Wagner, J., Benisty, D., & Karachentsev, I. D. (2026). The binary ballet: Mapping local expansion around M 81 and M 82. Astronomy & Astrophysics, 706, A92.
- Verde, L., Treu, T., & Riess, A. G. (2019). Tensions between the early and late universe. Nature Astronomy, 3, 891–895.
- Sharma, V. S. (2026). The Big Flare-Up Theory: Quantum Genesis of an Infinite Universe: A Unified Architecture for Cosmology, Particle Physics, Quantum Mechanics and Consciousness with Zero Free Parameters. Zenodo. DOI: 10.5281/zenodo.19149785
- Arp, H. (1966). Atlas of peculiar galaxies. Astrophysical Journal Supplement, 14, 1.
- Toomre, A., & Toomre, J. (1972). Galactic bridges and tails. Astrophysical Journal, 178, 623–666.
- van der Marel, R. P., et al. (2012). The M31 velocity vector. Astrophysical Journal, 753(1), 8.
- Sharma, V. S. (2026). BFUT Simulation Suite: Proof-of-Concept Interactive Simulations. Zenodo. DOI: 10.5281/zenodo.19124509
- Colin, J., Mohayaee, R., Rameez, M., & Sarkar, S. (2019). Evidence for anisotropy of cosmic acceleration. Astronomy & Astrophysics, 631, L13.
- Labbe, I., et al. (2023). A population of red candidate massive galaxies ~600 Myr after the Big Bang. Nature, 616, 266–269.
- Curtis-Lake, E., et al. (2023). Spectroscopic confirmation of four metal-poor galaxies at z = 10.3–13.2. Nature Astronomy, 7, 622–632.
- Boylan-Kolchin, M. (2023). Stress testing ΛCDM with high-redshift galaxy candidates. Nature Astronomy, 7, 731–735.