A First Principles Derivation of the Galactic Radial Acceleration Relation from Measurement Quantization

Jody A. Geiger

Informativity Institute

Chicago, IL USA 60605

jodygeiger@informativity.org

ORCID: 0000–0001–5389–0447

Abstract

The empirical radial acceleration relation (RAR) already describes the coupling between baryonic structure and galaxy dynamics. This work asks whether its functional form and acceleration normalization can be derived from independent physics rather than inferred from rotation curves. Measurement Quantization (MQ) derives the same exponential RAR form from count-constrained gravitational realization: the acceleration scale follows from frame-resolved cosmological expansion, the square-root realization depth from the positive quadratic System Frame norm and Newtonian circular kinematics, and the exponential realization fraction from subdivision-consistent composition of unrealized gravitational availability under the Frames mapping. Every MQ component is therefore fixed upstream of galaxy data, with present-epoch

aMQ=1.2004049×10-10ms-2
(S1)

and no galaxy specific MQ fit. Across 175 SPARC galaxies (3,366 radial measurements), the point-weighted signed residual is +0.02644 km s-1, the mean realized/observed velocity ratio is 1.00019, and pooled RMS and MAE are 22.82 and 16.20 km s-1. For 99 highest-quality galaxies, the geometric ratio is 1.0042. A RAR/MOND-independent Milky Way test reduces Newtonian RMS from 69.03 to 8.90 km s-1. Propagated McMillan source covariance gives generalized χ2=39.57 for 38 degrees of freedom (p=0.400). An independent 15-galaxy LITTLE THINGS test gives 0.9907 with a 95% galaxy block bootstrap interval of 0.868–1.130 and p=0.895. Because the acceleration scale inherits MQ cosmic age dependence, gravitational realization is predicted to be stronger at earlier epochs, providing a falsifiable mechanism relevant to early galaxy formation.

Keywords: Measurement Quantization, general relativity, cosmological constant, emergent spacetime, gravity, fundamental measure

1. INTRODUCTION

One of the most persistent unresolved problems in modern cosmology is the discrepancy between the gravitational behavior of galaxies and the mass inferred from directly observed baryonic matter. The central objective of this work is to determine whether that discrepancy can be derived, without a non-baryonic halo or galaxy-specific MQ fitting parameters, from a gravitational realization law fixed upstream of the galaxy data—one that recovers the Newtonian baryonic response at high acceleration but approaches a

1/R
(S2)

realized gravitational response for a finite baryonic source in the weak-acceleration outer limit. The latter behavior produces an asymptotically radius-independent orbital velocity and the baryonic Tully-Fisher fourth-power scaling. Measurement Quantization (MQ) develops this result from discrete relational count structure and the Frames mapping [1]. The galactic acceleration scale is derived from independently measured non-galactic inputs and fixed before rotation curve testing. The positive quadratic spatial norm in the System Frame, together with Newtonian circular kinematics, fixes the dimensionless local realization depth, while subdivision-consistent realization of the gravitational availability remaining unrealized yields the closed-form realization fraction. The resulting galactic realization law is then applied unchanged across galaxies spanning markedly different morphologies, luminosities, gas fractions, and surface-brightness regimes.

The motivation is therefore not to introduce another empirical rotation curve interpolation. The empirical RAR functional form is already established observationally. The question is whether MQ can derive that pre-existing form, including its acceleration normalization and asymptotic behavior, from independently established physical quantities while leaving the adopted Newtonian baryonic source unchanged. In MQ, the acceleration scale follows from the frame-resolved cosmological expansion, the square-root realization depth from the positive quadratic System Frame norm and circular kinematics, and the exponential realization fraction from subdivision-consistent composition of unrealized gravitational availability under the Frames mapping. The departure from the Newtonian baryonic prediction therefore arises from bounded realization at low acceleration, with the Newtonian result recovered continuously at high acceleration.

The goals of the paper are consequently threefold. First, it derives the complete galactic realization law from these independently established MQ components rather than selecting its functional form from galaxy rotation data. Second, it tests the resulting fixed present-epoch law against independent Milky Way, SPARC, and LITTLE THINGS baryonic-source and kinematic constructions with no galaxy-specific MQ calibration. Third, because the acceleration scale inherits the cosmic age dependence of the MQ expansion construction, it derives a redshift dependent extension in which gravitational realization is stronger at earlier epochs. This provides a prospective, falsifiable connection between the same realization mechanism tested in present-day galaxies and the unexpectedly rapid emergence of dynamically organized galaxies at high redshift. The present work does not claim a complete theory of galaxy formation. It establishes the age-dependent gravitational prediction that such a theory must test. Key MQ terms and symbols used throughout are summarized in Appx. A.

The present-epoch prediction is tested against three observational constructions. The first uses the 38 Milky Way circular velocity measurements of Eilers et al. [2] with the RAR/MOND-independent Galactic baryonic source of McMillan [3]. The second uses the complete SPARC database of 175 galaxies comprising 3,366 radial velocity measurements [4]. The third is an out-of-SPARC LITTLE THINGS dwarf galaxy test constructed from the independent H I and stellar-source work of Oh et al. [5] and the published rotation curve plots of O’Brien et al. [6]. In every case, the present-epoch acceleration scale and galactic realization law are fixed before the observational target is evaluated, no galaxy specific MQ parameter is fitted, and the externally adopted baryonic source quantities are not adjusted to improve MQ agreement. For the Milky Way, application of the galactic realization law reduces the RMS residual relative to the independent Eilers measurements from 69.03 km s-1for the Newtonian baryonic curve to 8.90 km s-1. SPARC provides the complementary broad population test across disk galaxies spanning large ranges of morphology, luminosity, surface brightness, and gas content [4], while LITTLE THINGS tests whether the fixed MQ normalization and SPARC residual velocity behavior persist in an out-of-SPARC galaxy population.

Predictive chain and scope. Table 1 separates the external inputs, the quantities derived before the galaxy data are examined, the fixed prediction applied to each baryonic source, and the observational tests. For a direct audit of the galactic result, the shortest reading path is Sec. 3.8 for the invariant coefficient, Sec. 3.12 for the frame-resolved expansion rate, Sec. 4.2 and Appx. D for the fixed galactic realization law, and Secs. 4.3-4.10 for the Milky Way, SPARC, and LITTLE THINGS tests and controls.

Table 1. Predictive dependency and test structure for the galactic result.

Stage

Quantity or information

Role in the present test

Location

External inputs

Defined speed of light c; measured fine structure constant used to determine the dimensionless MQ invariant coefficient θsi; adopted cosmic-age normalization; published baryonic mass models

These are upstream physical or astrophysical inputs. Observed orbital velocities are excluded from construction of the MQ prediction and are retained for comparison.

Secs. 3.8, 3.12, and 4.1; Refs. [2], [4], [7], [3]

Derived before galaxy testing

θsi, the Internal Frame expansion rate Hint, and

aMQ=cHint/(2θsi-1)
(S3)

The present-epoch acceleration scale is derived within MQ from the measured non-galactic inputs listed above, fixed before either Milky Way or SPARC velocities are evaluated, and not refitted by galaxy.

Eqs. (76), (126), and (131)

Fixed local prediction

Realization depth, subdivision-consistent realization fraction, Newtonian baryonic acceleration gbar(R), realized System Frame acceleration greal(R), and realized orbital velocity vreal(R), where R is galactocentric radius

The same galactic realization law is applied to every adopted baryonic source with zero galaxy-specific MQ calibration. The empirical RAR shape is pre-existing; the tested distinction is MQ provenance and the independently specified scale.

Eqs. (134–137); Appx. D; Ref. [8]

Observational tests

RAR/MOND-independent Milky Way source/kinematic comparison; complete 175-galaxy SPARC catalog; quality, stellar source, and galaxy-property controls; out-of-SPARC LITTLE THINGS population

The Milky Way and SPARC data test the fixed prediction after it has been specified; they do not determine aMQ or any galaxy-specific MQ parameter. The LITTLE THINGS test uses 15 prospectively screened non-SPARC galaxies and 288 radial velocity comparisons reconstructed from published observational/source products with an independent SED stellar mass normalization. Its galaxy-weighted geometric MQ/observed ratio is 0.9907 with a 95% galaxy-block bootstrap interval of 0.868-1.130; no fitted external RAR scale enters the test.

Secs. 4.3-4.10; Refs. [5–6,9]

This separation also defines the scope boundary for the principal claim: rotation curve agreement is not evidence that MQ succeeds at lensing, cluster merger, or cosmological perturbation scales, which remain independent physical requirements discussed in Sec. 5. For LITTLE THINGS, same-galaxy SPARC overlaps, incomplete stellar source cases, and source products unavailable for the figure reconstruction are excluded before the MQ residual is examined, leaving 15 galaxies and 288 radial comparisons. The observed circular velocities remain external to construction of the baryonic source and MQ prediction, and no fitted RAR acceleration parameter is used as a surrogate for that direct comparison.

For the statistical quantities reported below, r denotes a Pearson correlation coefficient, rP denotes the residual velocity Pearson correlation coefficient where that subscript is used, ρ denotes the Spearman rank correlation coefficient, p denotes the associated two-sided significance probability, and χ2 denotes the chi-square statistic. The externally assigned SPARC quality flag is denoted by Q, with Q = 1, Q = 2, and Q = 3 identifying the high-, medium-, and low-quality rotation curve classifications, respectively [4]. These classifications are observational metadata assigned independently of the present MQ analysis.

Table 2. Principal rotation curve validation statistics. Milky Way and SPARC results are reported separately because the two datasets employ different observational and baryonic model constructions.

Dataset

Scope

Absolute and curve-shape agreement

Velocity normalization

Population diagnostic

Milky Way (McMillan source)

38 independent Eilers circular velocity measurements; R=5.27-24.82 kpc

MQ RMS = 8.90 km s-1; MAE = 7.43 km s-1; Newtonian baryonic RMS = 69.03 km s-1

RMS discrepancy reduced by 87.11%; mean signed residual = -3.07 km s-1

Inner 30 points RMS = 6.88 km s-1 versus outer 8 points RMS = 14.09 km s-1; RAR/MOND-independent McMillan source. Diagonal total error χ2=689.831 (

p=1.26×10-120
(S4)

), but propagation of the published McMillan baryonic-normalization covariance gives generalized χ2=39.57 for 38 degrees of freedom (p=0.400). The Eilers observational covariance is not published; see Sec. 4.3.

SPARC, complete sample

175 galaxies; 3,366 radial points

Pooled RMS = 22.82 km s-1; pooled MAE = 16.20 km s-1; median galaxy Pearson r=0.964

Pooled bias = +0.0264 km s-1; ratio of mean MQ realized to mean observed velocity = 1.00019

Significant velocity-dependent normalization residual; it remains significant in simultaneous quality/property control, but physical attribution is not identifiable because characteristic velocity is strongly collinear with luminosity.

SPARC,
Q=1

99 catalog-classified high-quality galaxies

Median galaxy Pearson r=0.955

Geometric MQ-to-observed velocity ratio = 1.0042; mean signed bias = -0.59 km s-1

Residual velocity relation is not statistically significant: rP=-0.134, p=0.187.

LITTLE THINGS, out-of-SPARC figure-extracted sample

15 non-SPARC dwarf galaxies; 288 radial comparisons; independent Oh et al. SED stellar normalization [5–6]

Pooled RMS = 12.55 km s-1; pooled MAE = 8.69 km s-1; point-weighted signed residual = -0.39 km s-1

Galaxy-weighted geometric MQ/observed ratio = 0.9907; galaxy-block-bootstrap 95% interval = 0.868-1.130; zero mean-log-offset test p=0.895

No significant galaxy-level normalization dependence on log10 characteristic velocity: Pearson r=-0.282, p=0.308; Spearman ρ=-0.211, p=0.451. Partial control for SED stellar mass and inclination remains nonsignificant (r=-0.257, p=0.355). Material individual-galaxy and curve-shape scatter remains.

The Milky Way and SPARC statistics quantify different aspects of the same fixed realization law and should not be conflated. In the Milky Way analysis, the observational rotation curve and adopted baryonic source are taken from distinct published constructions: the Eilers et al. circular velocity measurements provide the observational comparison [2], while the McMillan Galactic bulge, thin-disk, thick-disk, H I, and molecular gas components provide the Newtonian baryonic source [3].

Across the complete 5.27–24.82 kpc range, MQ reduces the Newtonian baryonic RMS discrepancy by 87.11%, with a mean absolute error (MAE) of 7.43 km s-1. The agreement is radially nonuniform: the inner 30 measurements from 5.27–19.71 kpc have an RMS residual of 6.88 km s-1, whereas the outer eight measurements from 20.27–24.82 kpc have an RMS residual of 14.09 km s-1. The outer measurements also carry substantially larger observational uncertainties, giving them less weight under the adopted total error model despite their larger absolute residuals. The detailed pointwise and uncertainty analysis is developed in Sec. 4.3.

Across all 3,366 SPARC radial measurements, the point-weighted mean signed residual is +0.0264 km s-1, while the ratio of the mean MQ realized velocity to the mean observed velocity is 1.00019, corresponding to an aggregate velocity scale difference of approximately +0.019%. These near-zero aggregate quantities coexist with pooled RMS and MAE residuals of approximately 22.82 and 16.20 km s-1, respectively, and therefore do not imply point-by-point agreement at the same fractional level. Because radial measurements are clustered within galaxies, point-weighted aggregate quantities and galaxy-level population statistics also characterize different aspects of agreement.

At the galaxy level, the complete SPARC population exhibits a statistically significant dependence of the MQ-to-observed normalization offset on characteristic rotational velocity. Stratification by SPARC’s independently assigned quality classification materially changes that result. Among the 99 Q=1 galaxies, the geometric MQ-to-observed velocity ratio is approximately 1.0042 and the residual dependence on rotational velocity is not statistically significant, rP=-0.134, p=0.187. The complete sample velocity dependence is therefore not uniformly expressed across SPARC but is concentrated predominantly in the lower quality classifications.

This quality dependence does not establish that discrepancies among the lower quality galaxies are observational errors. SPARC quality is associated with observational difficulty and is also correlated in the present sample with astrophysical regime, including rotational velocity. The absence of a significant marginal velocity relation in the highest quality class therefore constrains, but does not remove, the complete catalog systematic. Section 4.8 adds a simultaneous galaxy-level partial regression test of quality, velocity, inclination, distance, luminosity, surface brightness, and gas richness. That joint test retains a significant conditional velocity term but also exposes strong velocity-luminosity collinearity, so the full-sample residual cannot presently be assigned uniquely either to the MQ realization law or to observational/source-systematic structure.

Across astrophysical systems spanning galaxies, galaxy clusters, and cosmological structure formation, gravitational observations provide evidence for a gravitating contribution beyond that attributable to directly observed baryonic matter alone [10]. Within the standard ΛCDM framework, this discrepancy is represented by a non-baryonic cold dark matter component, and cosmological fits infer a cold dark matter density substantially exceeding the baryonic matter density [11]. The particle identity of that component remains unresolved: despite extensive direct, indirect, and collider searches, no confirmed non-gravitational dark matter detection has been established [12–14]. Sections 1.2–1.4 review these observational constraints and the principal conventional responses in greater detail.

MQ investigates a different possibility by reconsidering the assumption that continuum spacetime geometry is fundamental. It interprets observable quantities as System Frame realizations of discrete Internal Frame count structure through the Frames mapping [1], building on the relational count framework established in earlier work [15]. The present paper adopts those results as its theoretical foundation. The MQ geometric construction relevant to gravitational realization is developed in Sec. 3, and its galactic consequences are tested in Sec. 4.

Within the emergent System Frame, length, mass, and time are realized quantities derived from the underlying count structure rather than primitive continuum inputs. The MQ invariant coefficient θsi is determined from the measured fine structure constant as developed in Sec. 3.8 and earlier work ([15], Appx. AD). Throughout this paper, untagged θsi denotes only this dimensionless invariant coefficient, while its momentum realization which has the same value is denoted psi. In the fundamental mass relation below, psi has dimensions of momentum, c denotes the defined speed of light, and mf denotes the fundamental mass:

mf=2psic=2.176432541(24)×10-8kg.
(1)

Section 3.3 makes the realization structure explicit. The untagged coefficient θsi is dimensionless; its angular realization Θsi carries the unit rad, while psi carries units of kg m s-1. These are dimensionally distinct realized representations of the same invariant MQ coefficient, with their dimensional character supplied in the System Frame. Their numerical coefficients coincide in the adopted MQ realization, as shown explicitly in Sec. 3.3. The Planck-like fine structure constant is denoted αp; prior to application of the MQ Informativity differential – a new length contraction effect a consequence of the discrete Internal Frame - in Sec. 3.6, its inverse is

αp-1=84θsi-⌈42θsi⌋=137.0407856(30),
(2)

where the mixed floor-ceiling brackets denote selection of the nearest whole-unit value using the corresponding floor or ceiling result.

Using the electron mass me, the Bohr radius a0, and αp, MQ then resolves the fundamental length lf and fundamental time tf:

lf=mea0αpmf=1.616199912(18)10-35m,
(3)
tf=lfc=5.391062611(59)10-44s.
(4)

The subscript p in αp identifies the Planck-like form of the fine structure constant. Equations (1–4) describe MQ fundamental or Planck-like quantities prior to complete finite-count realization and are therefore not direct predictions of the corresponding CODATA Planck quantities or of the classically realized fine structure constant. A concordance significance against those CODATA quantities would consequently compare different realization stages rather than like physical observables. The classical form – what is measured in the laboratory – is introduced only after derivation of the Informativity differential in Sec. 3.6, while Sec. 3.8 uses the 2022 CODATA recommended inverse fine structure constant as an input in determining θsi. That determination is therefore not an independent CODATA concordance test. Statistical concordance becomes physically meaningful only after the complete Frames mapping and Informativity differential have been applied and the resulting classical realization is compared with an independent laboratory measurement.

Observable physics emerges through the Frames mapping from discrete count configurations in the Internal Frame to continuous relational structure in the System Frame. Geometry is therefore an emergent property of realization rather than a fundamental background. The physically meaningful counts associated with length, mass, and time are represented by nL, nM, and nT, respectively, each corresponding to lf, mf, and tf within the Internal Frame. Their observable physical representations arise through the Frames mapping into the System Frame. MQ thus treats the fundamental measures as units relative to which discrete count configurations are defined rather than merely as dimensional combinations of conventional constants.

The galactic application concerns two well-established properties of rotation curves: orbital velocities commonly remain approximately constant or decline only slowly over extended galactocentric radii, and observed velocities exceed those predicted from the adopted baryonic distribution in systems exhibiting a mass discrepancy [4,16]. SPARC samples these behaviors across wide ranges of luminosity, morphology, surface brightness, and gas fraction [4], while the Milky Way supplies a complementary comparison using the independent Eilers et al. circular velocity measurements [2] and McMillan Galactic baryonic model [3]. Within MQ, these phenomena are attributed to bounded realization rather than to an independently postulated non-baryonic mass distribution.

The corresponding MQ construction is developed quantitatively in Sec. 4. Section 4.1 defines the observational baryonic inputs and Sec. 4.2 derives the local predictive realization law, with the complete local derivation in Appx. D. Sections 4.3 and 4.4 provide the RAR/MOND independent Milky Way and complete SPARC tests. Secs. 4.5–4.9 characterize the residual structure and apply observational quality, stellar source, and galaxy property controls; and Sec. 4.10 provides the combined interpretation.

The realized orbital velocity vreal(R) is obtained from the adopted Newtonian baryonic source through the MQ galactic realization law. Here vNew(R) denotes the Newtonian orbital velocity predicted from the adopted baryonic mass model:

vreal(R)=vNew(R)1-exp⁡[-gbar(R)/aMQ].
(5)

For the Milky Way, vNew(R) is reconstructed from the McMillan bulge, thin- and thick-disk, H I, and molecular gas density components by direct Newtonian force integration at the 38 Eilers radii, with the model dark halo excluded. The Eilers measurements and their uncertainties are used only for subsequent observational comparison. For SPARC, the same realization expression is evaluated without galaxy specific MQ adjustment for all 175 galaxies and 3,366 radial measurements before observational quality stratification. Equation (5) is the introductory form of the predictive relation developed in Sec. 4.2 and stated there as Eq. (137).

The gravitational coupling underlying the MQ formulation uses QL, the mapping-induced residual count contribution entering the realized separation under the Frames mapping; r, the radial separation between gravitating bodies; and G, the gravitational constant resolved at the upper-count limit. With psi as defined above, the relation is

QLc3rpsi=Gr2.
(6)

Equation (6) establishes the MQ correspondence between realized count geometry and the Newtonian gravitational field. Its derivation and realization dependence are developed in Secs. 3.2–3.6. The present section uses it only to identify the gravitational relation from which the galactic realization construction proceeds.

1.2 Observational evidence for missing gravitational mass

The empirical motivation for the present test is the established mismatch between reconstructed baryonic gravity and observed dynamics. Galaxy rotation curves provide the primary target [4,16], while gravitational lensing, merging clusters, and cosmological observations show that any general explanation must ultimately succeed outside the rotation curve regime [10–11,17]. The Milky Way comparison in Sec. 4.3 and the SPARC and LITTLE THINGS population tests developed through Secs. 4.4–4.10 are therefore the principal observational tests of this paper. The broader evidence defines independent requirements that MQ must eventually satisfy.

v=GMR,
(7)

where v is the orbital velocity and M the enclosed gravitating mass. G and R retain their definitions established above.

For later comparison with the cosmological benchmark, MQ defines universe domains independently of the ΛCDM matter-density decomposition. In MQ notation, Ωtot=1, Ωobs≈0.316, and Ωvis≈0.0488, where Ωtot denotes the normalized total universe domain, Ωobs the total observable domain, and Ωvis the presently visible domain. These are MQ geometric domain quantities, not identifications with standard cosmological density parameters despite their similar numerical magnitudes. Their derivation and physical interpretation are developed in Sec. 3.9 and Appx. C. The unobserved domain is the portion of Ωobs that is not presently referenceable,

Ωuobs=Ωobs-Ωvis.
(8)

In that all variables are constant, these geometric fractions remain constant across all epochs.

Let Mobs denote the total mass corresponding to Ωobs, Mvis the mass corresponding to Ωvis, and Muobs the complementary mass corresponding to Ωuobs. The MQ universe domain mass relations are

Mobs=2θsiMvis,
(9)
Muobs=Mobs-Mvis.
(10)

These relations distinguish the MQ domain construction from observational ΛCDM density accounting. They do not constitute independent observational evidence for the galactic result.

1.3 Rotation curves, lensing, CMB constraints, and cluster offsets

Four observational regimes bound the broader dark matter problem: galaxy dynamics, gravitational lensing, CMB/large-scale-structure inference, and merging cluster offsets [4,11,17]. Only galaxy dynamics is quantitatively tested by the principal analysis. The remaining regimes impose independent requirements developed prospectively in Secs. 5.3–5.9 and are not inferred from rotation curve agreement.

Gravitational lensing provides an observationally distinct probe through light deflection. Let

α̂
(S5)

denote the physical gravitational deflection angle and b the impact parameter. In the general relativistic weak-field point-mass limit, a photon passing a gravitating mass M is deflected according to

α̂=4GMc2b.
(11)

Strong- and weak-lensing reconstructions infer projected gravitating mass distributions that can substantially exceed those associated with observed baryonic components, particularly in galaxy clusters [17]. In merging systems such as the Bullet Cluster, lensing additionally constrains models through the observed separation between reconstructed gravitating mass and the dominant X-ray-emitting baryonic plasma [17]. Consequently, reproducing galactic rotation curves does not substitute for deriving the corresponding relativistic MQ lensing and cluster response.

The CMB provides a distinct cosmological constraint. Within standard spatially flat ΛCDM, let Ωm denote the total matter density fraction. The baryonic matter fraction is Ωb, and the cold dark matter fraction is Ωc. Neglecting the small distinction associated with massive-neutrino bookkeeping for this introductory decomposition,

Ωc=Ωm-Ωb.
(12)

Planck CMB measurements constrain the physical cold dark matter density to be substantially larger than the baryonic density, making this cosmological requirement independent of the galactic rotation curve calculation [11].

1.4 Limits of particle dark matter and modified gravity

The relevant comparison is predictive structure rather than whether alternative models can fit individual galaxies. ΛCDM introduces a non-baryonic component within a cosmological framework that fits a broad range of observations [11,18], whereas MOND-like approaches modify the low-acceleration response [19–21]. MQ instead retains the adopted Newtonian baryonic source and tests the fixed galactic realization law whose acceleration scale is determined upstream of the galaxy data. The absence of a confirmed particle dark matter identification and known small-scale tensions motivate investigation of alternatives [12–14,22].

For the standard nonrelativistic MOND formulation, let a denote the centripetal acceleration, aMOND the characteristic MOND acceleration scale, and

μ(a/aMOND)
(S6)

the interpolation function connecting the Newtonian and deep-MOND regimes. The relation is

μaaMONDa=GMr2.
(13)

1.5 The MQ realization of geometry and the origin of dark matter behavior

The specific proposal tested here is that the apparent excess galactic gravitational response arises from bounded realization of mapped count geometry rather than from an additional source mass. The same present-epoch galactic realization law is applied without galaxy-specific MQ calibration to the Milky Way, the complete SPARC catalog, and the out-of-SPARC LITTLE THINGS sample.

The gravitational coupling introduced in Eq. (6) relates realized count geometry to the Newtonian gravitational field through count separation and the Frames mapping. The relevant geometric construction is developed in Sec. 3, the independently specified acceleration scale and local galactic realization law in Sec. 4.2 and Appx. D, and the Milky Way, SPARC, and LITTLE THINGS tests and controls in Secs. 4.3–4.10. Relativistic lensing, cluster dynamics, cosmic-age evolution, and cosmological perturbations remain independent requirements developed in Secs. 5.3–5.9 rather than consequences of the galaxy results.

2. BACKGROUND AND COMPARATIVE FRAMEWORK

2.1 ΛCDM as a model benchmark

ΛCDM is used here as the quantitative benchmark because it jointly describes the CMB, large-scale structure, and cosmic expansion while representing the additional gravitating contribution with non-baryonic dark matter [11,18]. The comparison required for the present paper is not a global replacement test: Sec. 4 asks whether MQ reproduces galaxy dynamics with a fixed response law, while Sec. 5 identifies the independent calculations required before any broader equivalence claim is justified.

H2=8πG3ρ-kc2aFLRW2(t)+Λc23,
(14)

where H(t) is the Hubble expansion rate, ρ is the total mass-energy density expressed in mass-density units, k is the spatial-curvature parameter, aFLRW(t) is the cosmological scale factor, and Λ is the cosmological constant; G and c retain their definitions established above. For the spatially flat ΛCDM benchmark, the present-day normalized density budget is decomposed into radiation, baryonic matter, cold dark matter, and dark energy according to

Ωr+Ωb+Ωdm+ΩΛ=1.
(15)

Here, Ωr, Ωb, Ωc, and ΩΛ denote the radiation, baryonic matter, cold dark matter, and cosmological constant density fractions, respectively. Planck 2018 base-ΛCDM constraints [11] correspond approximately to Ωb≃0.049, Ωc≃0.266, and ΩΛ≃0.685, with the present radiation contribution comparatively small. Thus, most of the present cosmic energy-density budget is not contained in electromagnetically luminous baryonic matter.

M(r)=Mbar(r)+Mdm(r),
(16)

where M(r) denotes the total enclosed gravitating mass, Mbar(r) the enclosed baryonic mass, and Mdm(r) the enclosed dark matter halo mass. Individual galaxy analyses then specify a halo density profile together with the baryonic mass model and observational parameters required to predict the corresponding rotation curve.

2.2 A fundamental accounting of galactic rotation curves

Galactic rotation curves expose a simple accounting problem: orbital velocities observed at large galactocentric radii generally exceed those predicted from the reconstructed baryonic mass alone. ΛCDM and MQ begin from the same discrepancy but assign it to different physical mechanisms. In ΛCDM, the additional gravitational response is attributed predominantly to a non-baryonic halo, so the gravitating source contains baryonic and dark components. MQ instead retains the Newtonian baryonic source: the adopted baryonic reconstruction specifies the gravitational availability, and the galactic realization law determines the realized gravitational response in the System Frame.

The fixed MQ construction uses the Newtonian baryonic acceleration, the independently determined MQ acceleration scale, the positive quadratic norm supplied by the realized geometry, and subdivision-consistent composition of incomplete realization. These determine the local realization depth and realization fraction used in the rotation curve prediction. Thus ΛCDM modifies the gravitating source through a halo mass distribution, whereas MQ leaves the baryonic source unchanged and tests whether a fixed galactic realization law reproduces the observed dynamics without a separately fitted non-baryonic component.

2.3 Gravitational lensing and cluster mass reconstruction

Lensing is an independent relativistic benchmark rather than part of the galactic fit. A lens reconstruction distinguishes the contribution associated with reconstructed baryonic matter from any additional gravitating contribution required by the adopted lens model. MQ must ultimately reproduce the corresponding deflection, convergence, and shear observables from the realized System Frame metric without introducing a lens-specific dark component. Section 5.3 develops the present spherical weak-field and thin-lens results and identifies the remaining extended-source covariant calculation.

For comparison among physical interpretations, let Mlens denote the total gravitating mass inferred from a lensing reconstruction, Mbar,lens the contribution associated with reconstructed baryonic matter, and Mexc,lens the additional gravitating contribution required by that reconstruction. The observational accounting is

Mlens=Mbar,lens+Mexc,lens.
(17)

Equation (17) is an observational accounting relation only and does not assign a physical identity to Mexc,lens. The quantities Mbar,lens and Mexc,lens are local lens-reconstruction quantities and are distinct from the MQ universe domain mass quantities developed in Appx. C.

2.4 Bullet cluster spatial offsets

Merging clusters provide a stringent nonequilibrium spatial test because the collisionless galaxy population, hydrodynamic intracluster gas, and reconstructed lensing mass distribution can separate during a merger [17,23–24]. The Bullet Cluster provides the canonical example, with the dominant X-ray-emitting plasma spatially displaced from the principal lensing mass peaks [17]. MQ must therefore evolve the baryonic system and realized metric together. The equilibrium galactic realization law alone is insufficient. Section 5.4 formulates the corresponding time-dependent problem and makes explicit that the nonlinear Bullet Cluster evolution and parameter-frozen lensing comparison remain to be performed.

Let Mbar,cl denote the reconstructed baryonic cluster mass represented here by the galaxy and gas contributions Mgal and Mgas. Then

Mbar,cl=Mgal+Mgas.
(18)

Let rlens,⊥ denote the projected centroid vector of a reconstructed lensing-mass peak and rgas,⊥ the projected centroid vector of the corresponding X-ray gas component. Their plane-of-sky separation Δr⊥ is

Δr⊥=|rlens,⊥-rgas,⊥|.
(19)

This separation is an observational diagnostic. Eqs. (18–19) do not by themselves provide a dynamical model of the merger.

2.5 Tensions challenging ΛCDM

Small-scale ΛCDM tensions provide comparative context. Modern baryonic galaxy formation calculations alleviate several historically identified discrepancies, while dwarf-galaxy inner structure, satellite populations, and rotation curve diversity remain active tests [25–29]. The present manuscript therefore does not treat unresolved ΛCDM phenomenology as validation of the MQ mechanism. Modern reviews similarly emphasize that distinguishing limitations of collisionless-only predictions from effects of baryonic physics or possible dark matter microphysics remains an active problem.

One long-standing example is the cusp-core problem. Collisionless cold dark matter simulations yield centrally concentrated halo profiles represented by the Navarro-Frenk-White form [25]. Let ρNFW(r) denote the NFW halo density at halo-centric radius r, ρs its characteristic density, and rs its scale radius. The profile is

ρNFW(r)=ρs(r/rs)(1+r/rs)2.
(20)

The NFW relation provides the collisionless benchmark. Whether a particular observed inner profile requires baryonic modification, altered dark matter physics, or another explanation must be assessed separately.

2.6 An assessment of leading models with well-established observational constraints

The comparative standard is deliberately asymmetric: ΛCDM is quantitatively established across cosmology, whereas MQ is tested here most strongly in galaxies. Comparison 1 therefore distinguishes observationally demonstrated results from compatibility within viable parameter ranges and from calculations that remain unperformed. Alternative dark matter and Modified gravity models are retained as context for the observational requirements and predictive freedom against which the fixed MQ galactic realization law should be judged.

Rotation curves, gravitational lensing, cluster dynamics, CMB anisotropies, and large-scale structure constrain the inferred gravitating contribution through different observables and physical scales. The leading frameworks do not address all of these constraints with equal quantitative maturity. At the observational accounting level, let Mgrav denote the total gravitating mass inferred under the adopted gravitational reconstruction, Mbar,obs the corresponding observationally reconstructed baryonic contribution, and Mexc the additional inferred gravitating contribution. Then

Mgrav=Mbar,obs+Mexc.
(21)

Equation (21) is an observational bookkeeping relation and does not by itself assign a physical identity to Mexc. Within ΛCDM, the inferred excess is attributed predominantly to non-baryonic dark matter. Modified gravity theories instead alter the relation between the baryonic source and gravitational response. MQ retains the Newtonian baryonic source and tests whether the inferred excess arises from the realized gravitational response produced through bounded realization under the Frames mapping.

Comparison 1. Qualitative comparison of leading dark matter and modified gravity model classes against the principal observational constraints. “Compatible” denotes compatibility within viable model parameter space and should not be read as an independent quantitative validation of every realization.

ΛCDM / collisionless cold dark matter

Rotation curves: Halo modeling provides the standard successful phenomenological description of galaxy rotation curves.

Lensing and cluster masses: Non-baryonic halo mass supplies the additional gravitating component required by standard lensing and cluster reconstructions.

Bullet Cluster–type offsets: Collisionless dark matter can remain spatially associated with galaxies while ram-pressure-stripped gas is displaced, consistent with the standard interpretation of merging clusters [17].

CMB and large-scale structure: Quantitative benchmark for CMB anisotropies and large-scale structure [11,18].

Small-scale tensions: Several classical discrepancies are alleviated by baryonic galaxy formation, while inner-halo structure and diversity remain active tests [29].

Direct / indirect detection: No confirmed dark matter particle identification; broad classes of candidate parameter space are constrained [22].

WIMP dark matter

Rotation curves: In the cold, effectively collisionless regime, inherits standard CDM halo phenomenology.

Lensing and cluster masses: Compatible with the standard ΛCDM lensing interpretation when the WIMP behaves as cold collisionless matter.

Bullet Cluster–type offsets: Compatible when self-interactions are negligible on cluster merger scales.

CMB and large-scale structure: Thermal and non-thermal realizations are constrained by relic-density, CMB, and structure requirements; viable cold realizations reproduce CDM-like large-scale behavior.

Small-scale tensions: Generally inherits CDM structure issues unless baryonic physics or additional particle interactions alter halo evolution.

Direct / indirect detection: Direct, indirect, and collider searches strongly constrain broad WIMP parameter regions, with no confirmed signal [22].

QCD axion / cold axion-like dark matter

Rotation curves: Cold axion realizations can reproduce CDM-like halo dynamics over viable parameter ranges.

Lensing and cluster masses: Compatible with standard gravitational lensing when behaving as a cold gravitating component.

Bullet Cluster–type offsets: Sufficiently weakly interacting cold axions behave effectively collisionlessly during cluster mergers.

CMB and large-scale structure: Vacuum-realignment axions can constitute cold dark matter; cosmological abundance and isocurvature constraints restrict viable models [30].

Small-scale tensions: Depends on mass, production history, and interactions; the ultralight regime has distinct wave-scale phenomenology and is listed separately below.

Direct / indirect detection: Extensive haloscope, helioscope, laboratory, and astrophysical searches constrain axion parameter space; no confirmed dark matter axion signal has been established [30].

Sterile-neutrino / warm dark matter

Rotation curves: Free streaming suppresses low-mass structure relative to CDM and can alter dwarf-scale halo populations [31].

Lensing and cluster masses: Must retain sufficient clustering to reproduce galaxy- and cluster-scale lensing.

Bullet Cluster–type offsets: Viable realizations remain effectively collisionless on merger scales.

CMB and large-scale structure: Production mechanism and particle mass are constrained by structure formation, CMB-related cosmology, Lyman-α measurements, and radiative-decay limits [31].

Small-scale tensions: Reduced small-scale power can lessen some low-mass-halo discrepancies but can also suppress structure beyond observationally allowed levels.

Direct / indirect detection: No confirmed signal; X-ray decay searches and cosmological structure measurements constrain mass and mixing [31].

Self-interacting dark matter

Rotation curves: Self-scattering can alter halo cores and reproduce a broader diversity of inner rotation curves than collisionless-only halos [32].

Lensing and cluster masses: Viable models are constrained by cluster density profiles, strong and weak lensing, and halo-shape measurements [32].

Bullet Cluster–type offsets: Merging clusters directly constrain the magnitude and velocity dependence of the self-interaction cross section [32].

CMB and large-scale structure: Viable interaction models must preserve the successful near-CDM large-scale matter distribution.

Small-scale tensions: Can lower central densities and modify halo diversity; baryonic degeneracies remain important [32].

Direct / indirect detection: No confirmed underlying dark-sector particle signal; astrophysical observations primarily constrain the effective interaction strength.

Fuzzy / ultralight dark matter

Rotation curves: Wave dynamics can produce solitonic or cored inner structures, depending on particle mass and composition [33].

Lensing and cluster masses: Halo substructure and lensing constrain suppression of small-scale structure.

Bullet Cluster–type offsets: Viable models must reproduce observed cluster-scale collision and mass-distribution behavior.

CMB and large-scale structure: Large scales can approach CDM behavior, whereas characteristic suppression of small-scale power is tightly constrained by Lyman-α, halo-abundance, and related probes [33].

Small-scale tensions: Potentially modifies cusp-scale structure, but parameter ranges capable of doing so are restricted by independent structure constraints [33].

Direct / indirect detection: No confirmed signal; constraints arise from gravitational, astrophysical, black-hole, pulsar, and laboratory-sensitive effects depending on mass and coupling [33].

MOND / relativistic modified gravity

Rotation curves: The MOND acceleration law provides strong phenomenology for many galaxy rotation curves and baryon–acceleration correlations [g,n].

Lensing and cluster masses: Relativistic extensions can generate gravitational lensing, but galaxy clusters generally retain a residual missing-mass requirement beyond the basic MOND prescription [21].

Bullet Cluster–type offsets: Merging-cluster observations remain a demanding constraint and generally require additional fields, unseen mass, or a sufficiently complete relativistic model [21].

CMB and large-scale structure: Original nonrelativistic MOND has no standalone relativistic cosmology; cosmological predictions require a specified covariant completion with additional dynamical structure [21].

Small-scale tensions: Naturally reproduces several galaxy-scale baryon–gravity regularities; its remaining difficulties occur primarily beyond the isolated-galaxy regime [21].

Direct / indirect detection: Not intrinsically a particle dark matter model and therefore has no generic dark-particle detection target.

MQ / bounded realization

Rotation curves: The fixed galactic realization law is applied to 175 SPARC galaxies comprising 3,366 radial points and to a prospectively screened out-of-SPARC LITTLE THINGS sample of 15 galaxies comprising 288 figure-extracted radial comparisons. LITTLE THINGS gives a galaxy-weighted geometric MQ/observed velocity ratio of 0.9907, with a 95% galaxy-block bootstrap interval of 0.868–1.130 and no significant mean-log offset (p=0.895); the characteristic-velocity normalization relation is also nonsignificant (rP=-0.282, p=0.308), while substantial pointwise scatter remains [5–6,9].

Lensing and cluster masses: Sec. 5.3 provides a spherical ESO 325–G004 benchmark: the frozen MQ spherical calculation predicts an Einstein radius of approximately 2.67 arcsec versus the observed 2.85 arcsec. The full resolved extended source metric/slip solution and two-dimensional deflection, convergence, and shear comparison remain untested.

Bullet Cluster–type offsets: Sec. 5.4 formulates the time-dependent merger problem, while Sec. 5.8 provides the action-level closure used by the broader field construction. No nonlinear Bullet Cluster evolution or parameter-frozen lensing-map comparison has yet been performed; the observed offsets therefore remain an open test.

CMB and large-scale structure: Sec. 3.9 derives MQ universe domain fractions numerically near the conventional density fractions, but these are not identified with fitted ΛCDM fluids. Section 5.8 develops the field/conservation linearization and action-level closure required to pose the perturbative problem, but no numerical Einstein–Boltzmann integration or prediction of CMB temperature/polarization spectra, CMB lensing, BAO, the matter power spectrum, or cosmological growth observables has yet been completed.

Small-scale tensions: In the QSPARC=1 population, no significant residual dependence is recovered on several tested galaxy properties. This is a rotation curve control and is not itself a test of the conventional cusp-core, satellite, or halo-substructure problems.

Direct / indirect detection: MQ is not a particle dark matter model and therefore predicts no dark matter particle detection.

3. MQ FOUNDATIONS

3.1 Right–triangle construction and count geometry

Figure from the preprint

Figure 1. Count of lf describing realized separation in the System Frame, with discrete count evaluated in the Internal Frame. Modified from [18].

The count structure underlying the Frames Principle is made explicit through the right–triangle realization construction shown in figure 1. The triangle is not a geometric structure within the Internal Frame. It is the realized representation used to relate discrete count structure to measurable separation under the Frames mapping. In this representation, side a specifies the unit reference count and is fixed to unity, establishing the normalization used for count comparison. Side b corresponds to the variable Internal Frame count nL. The labels side a, side b, and side c identify roles within this construction and do not introduce additional physical frames.

Side c, the hypotenuse, represents the realized separation l in the System Frame. The mapping-consistent count combination entering this separation combines the Internal Frame count nL with the mapping-induced residual count contribution QL, giving nL+QL. Multiplication by the fundamental length lf then gives the realized physical separation l=(nL+QL)lf where l denotes the realized separation in the System Frame, lf the fundamental length, nL the Internal Frame count, and QL the mapping-induced residual count contribution entering the realized separation. The right–triangle construction determines QL through

QL=1+nL2-nL
(S7)

.

Thus, QL is not an independent primitive count but is fixed uniquely by nL within the realization construction. As nL increases, QL decreases while remaining positive, providing the finite residual contribution required for mapping-consistent realization of the count relation. The same construction satisfies the right–triangle identity 1+nL2=(nL+QL)2.

The orthogonality of side a and side b belongs to this realized right–triangle representation and encodes the algebraic independence of the unit reference contribution and the variable count contribution used in constructing the hypotenuse. It does not assign an intrinsic Euclidean geometry to the Internal Frame. The hypotenuse instead provides the mapping-consistent realized separation through the count combination nL+QL.

The role of this representation is operational. It defines the count relation from which subsequent MQ separation constructions are developed. A separation measured in the System Frame is therefore a mapping-normalized realization of conserved Internal Frame count rather than a direct geometric projection of that count. Consequently, the realized separation contains the residual contribution QL required by the Frames mapping. The gravitational significance of this residual is developed in Sec. 3.2, while its relation to the emergent continuum metric is developed independently in Sec. 3.5.

3.2 Gravitational curvature using the right-angle triangle construct

The right-angle triangle construct provides the count geometry used to describe gravitational realization. Side b carries the Internal Frame count nL, while side c represents the mapping-consistent realized separation in the System Frame. The mapping-induced residual count contribution QL is uniquely determined by nL and preserves the exact right-triangle relation under the Frames mapping.

MQ identifies this mapping-induced deformation as the geometric origin from which gravitational curvature is subsequently realized. At this stage, “curvature” denotes the MQ count-resolved scalar measure associated with gravitational coupling. It is not the tensorial curvature of the emergent continuum metric developed in Sec. 3.5. In the curvature-resolved regime, the scalar measure is constructed from the dimensionless residual ratio

QL/nL
(S8)

and the fundamental acceleration scale

lf/tf2
(S9)

:

QLlfnLtf2.
(22)

Equation (22) preserves the count structure explicitly. The factor

QL/nL
(S10)

gives the fractional residual contribution associated with count separation, while

lf/tf2
(S11)

supplies the dimensional acceleration scale. No continuum radial coordinate is required at this stage.

To express the relation as a gravitational field per unit source mass, MQ applies the normalization factor

c/psi
(S12)

, where psi is the momentum realization of the dimensionless invariant coefficient θsi. Because

c/psi
(S13)

has dimensions of inverse mass, multiplication of Eq. (22) converts the acceleration measure to the dimensions of

G/r2
(S14)

. Using

c=lf/tf
(S15)

and defining the radial count coordinate associated with side b as r=nLlf, the expression transforms as

QLlfnLtf2(cpsi)=QLc2nLtfpsi=QLc3rpsi.
(23)

The coordinate r=nLlf is the radial count coordinate associated with side b. It is distinct from the complete mapping-consistent realized separation l=(nL+QL)lf associated with side c.

The resulting MQ scalar measure may then be compared dimensionally with the Newtonian inverse-square gravitational field per unit source mass,

[QLc3rpsi]=[Gr2],
(24)

where

G/r2
(S16)

has the same dimensions as the normalized MQ expression. Equation (24) states only an exact dimensional equivalence. It is neither a physical equality nor an approximate dynamical law. The physical relation is stated separately in Eq. (25).

Throughout this paper, upright constants denote quantities evaluated at upper-count-limit realizations, while italicized constants denote quantities evaluated at the electromagnetic demarcation (Appx. B) or, more generally, at the applicable finite-count realization scale. MQ predicts that electromagnetic determinations are most sensitive to the Informativity differential, whereas that contribution becomes negligible at sufficiently large count separation. When both sides are evaluated within the same realization regime, the corresponding MQ gravitational relation is

QLc3rpsi=Gr2.
(25)

Equation (25) is the canonical MQ count-resolved gravitational coupling relation used throughout the remainder of the paper. It establishes the scalar physical relation between realized count geometry and the Newtonian gravitational field in the curvature-resolved regime. The tensorial metric interpretation is developed separately in Sec. 3.5.

3.3 Geometric correlation of electromagnetic and gravitational phenomena

The dimensional distinction must be fixed before the common MQ coefficient is used in a physical relation. Untagged θsi denotes the invariant dimensionless coefficient. Its angular realization is Θsi and its momentum realization is psi. The realization operators RL and Rp denote the angular and momentum representations of θsi, respectively; the realized representation, not θsi itself, supplies the physical dimensions. The reduced Planck constant is denoted by ℏ. The definitions are

θsi=3.262390305(36)∈R.
(26)
Θsi:=RL(θsi)=θsi​rad.
(27)
psi:=Rp(θsi)=lfmf2tf=ℏ2lf.
(28)
Cpθ:=psiΘsi.
(29)
psi=CpθΘsi.
(30)

Equations (26–30) are representation definitions, not equalities among unlike dimensions. The angular and momentum realizations inherit the same invariant numerical coefficient through their respective realization operators. In the adopted SI realization, their displayed numerical coefficients therefore coincide with θsi, while Cpθ supplies the dimensional conversion between angle and momentum. A plane angle and a momentum are not thereby identified as the same physical quantity.

The gravitational coupling relation may be reduced at the upper-count limit using

QLnL→1/2
(S17)

, as shown in Fig. 2. With r=nLlf, the relation becomes

c3G=limnL→∞psiQLnLlf=2psilf
(31)

Equation (31) is written dimensionally with psi, not θsi. It identifies the momentum realization as the quantity relating the upper-count gravitational normalization to the fundamental length scale. The invariant coefficient carried by that realization remains θsi. The finite-count gravitational relation remains Eq. (25).

Using the MQ fundamental length relation, Eq. (31) gives

psi=lf2(c3G)=ħ2lf.
(32)

Equation (32) has the required dimensions of momentum and does not assign momentum dimensions to θsi. Note that the reduced Planck and gravitational constants are not italicized. All symbols are taken at the upper-count limit.

Figure from the preprint

Figure 2. Arc-length construction for a circle of radius lf. The source schematic’s local label θ=S denotes the illustrated plane angle and its numerical radian measure. It is distinct from the MQ invariant coefficient θsi. Modified from [18].

(S18)

A distinct angular realization follows from the conventional plane-angle construction illustrated in Fig. 2. The angular realization Θsi, defined by Eq. (27), carries the unit radians. The corresponding arc length Lsi on a circle of radius lf is

Lsi=lfΘsi.
(33)

Figure 2 makes the dimensional distinction explicit: θsi is the invariant dimensionless coefficient, Θsi its angular realization, and psi its momentum realization. Their common numerical coefficient expresses the shared invariant carried through realization, not an equality of physical dimensions.

These realizations provide the geometric correlation relevant to electromagnetic and gravitational phenomena. MQ does not identify electromagnetism and gravitation as the same physical field, nor Θsi and psi as the same dimensional quantity. The momentum realization enters the gravitational coupling developed above, while the dimensionless and angular realizations enter the fixed-count and geometric constructions developed in Secs. 3.4 and 3.7. Their common coefficient therefore provides a structural correlation within the Frames mapping while the physical representations remain distinct.

Equating the two forms of psi in Eq. (32) gives

lfc32G=ħ2lf,
(34)
lf=(ħGc3)1/2.
(35)

The gravitational construction developed in Sec. 3.2 provides the MQ relation from which the Planck-scale relations and fundamental measures are subsequently derived. As established in Ref. [15], these results follow algebraically from the MQ gravitational coupling together with the MQ fundamental relation, rather than being introduced as independent Planck-scale assumptions. Equation (35) has the conventional Planck length form, constructed from the reduced Planck constant, gravitational constant, and defined speed of light [7]. This algebraic coincidence identifies the standard Planck-scale form recovered by the MQ reduction.

3.4 Electromagnetic and gravitational interaction demarcations
The gravitational curvature relation developed in Sec. 3.2 establishes that gravitational realization depends on the mapping-induced deformation associated with count separation. Within MQ, electromagnetic behavior is realized at a fixed interaction scale defined by the electromagnetic demarcation, whereas gravitational behavior occupies the curvature-resolved regime and is evaluated across varying count separations. These are scale-dependent realizations of the same underlying mapping-induced deformation. The electromagnetic demarcation is not a fixed electromagnetic wavelength, energy, propagation distance, or QED renormalization scale. It is the fixed interaction count distance at which electromagnetic realization is evaluated under the Frames mapping.

The electromagnetic demarcation is not a galaxy- or interaction-specific fitting parameter. It is derived from the invariant coefficient θsi, introduced in Sec. 3.3, together with the MQ count geometry. In the source construction used here, θsi and the MQ fundamental measures ultimately depend on the measured fine structure constant α, electron mass me, Bohr radius a0, and defined speed of light c. Quantities derived from this source set are therefore not subsequently treated as independent validation when they share that upstream dependence.

One-half of the electromagnetic interaction count is

nL2=84.600552≈42.3,
(36)

giving approximately 42.3. Its nearest whole-unit count is 42. This value provides the discrete reference associated with one-half of the electromagnetic demarcation and reappears in subsequent electromagnetic constructions, including the fine structure constant realization developed in Sec. 3.7.

The electromagnetic demarcation may first be obtained from the blackbody construction developed in Appx. B. Beginning with the gravitational curvature relation and substituting the count-separation relation

QL=1+nL2-nL
(S19)

, the construction gives

nL=θsi2lf2ℏ2(1-(2θsilf/ℏ))=θsilf1ℏ(ℏ-2θsilf).
(37)

Equation (37) defines the blackbody realization of the electromagnetic demarcation. Its complete derivation is given in Appx. B, Eqs. (B.1–B.9).

A second construction follows from the charge-coupling relation,

nL=276/θsi,
(38)

which expresses the electromagnetic demarcation directly as a function of the invariant coefficient θsi. Its derivation is given in Appx. B, Eq. (B.10).

A third, historically earlier construction follows from the elementary charge solution developed using the gamma approach. Here, ec denotes the classically realized MQ elementary charge, γ is the conventional electromagnetic transformation parameter used in the historical gamma approach, and rq41 denotes the historical approximate discrete-offset quantity used in that construction. In this approach, γ was deliberately resolved from a measured electromagnetic constant different from the quantity being calculated, avoiding direct use of the prediction target as its own input. The construction nevertheless depended on an additional experimentally measured electromagnetic source quantity rather than solely on MQ fundamental quantities. The corresponding demarcation relation is

nL((1+nL2)1/2-nL)=lfc22ec(2(84θsi-⌈42θsi⌋)γ)1/2.
(39)

The γ in Eq. (39) is distinct from the later fundamental gamma γF. The latter is a first principles realization parameter introduced after the elementary charge realization eliminated the earlier discrete-offset and auxiliary-γ construction. Unlike the conventional γ, which is resolved from experimentally measured electromagnetic constants, γF is derived entirely from MQ fundamental quantities and serves as the canonical electromagnetic realization parameter in later MQ derivations.

The blackbody, charge-coupling, and historical elementary charge constructions provide different derivational routes to the same electromagnetic interaction count. Their numerical agreement is an internal closure test of the shared MQ count geometry, not three statistically independent laboratory determinations. The historical route is additionally distinguished by its auxiliary measured electromagnetic input, whereas the later γF construction removes that dependence. Table 3 summarizes the three demarcation values.

Table 3. Electromagnetic demarcation obtained from three MQ constructions

Source construction

nL
(S20)

Fine structure-constant relation, α

84.6005457004(82)

Blackbody relation, ℏ

84.6005496647(07)

Historical elementary charge/γ relation, eᵃ

84.6005499213(87)

ᵃ Here e denotes the defined SI elementary charge [7]. The associated interaction count is solved as a function of the fine structure constant and this defined value. This historical gamma-based construction is retained as a cross-consistency route and is distinct from the later MQ γF realization.

The physical significance of the electromagnetic demarcation is not that electromagnetic processes occur at a single conventional physical distance. Electromagnetic realization under the Frames mapping is evaluated at the fixed interaction scale

nL=276/θsi
(S21)

, whereas gravitational phenomena occupy the curvature-resolved regime, in which realized geometric structure varies with count separation. Both arise from the same mapping-induced deformation but represent different scale-dependent realizations of that mapping effect. MQ therefore does not identify electromagnetic and gravitational phenomena as the same physical field or replace their conventional scale-dependent dynamics.

3.5 Emergent metric from coarse–graining

The right-triangle count construction developed in Sec. 3.1 establishes the mapping-consistent count combination nL+QL, where nL is the Internal Frame count and QL is the mapping-induced residual count contribution introduced by the Frames mapping. Multiplication by the fundamental length lf gives the corresponding realized separation l in the System Frame,

l=(nL+QL)lf.
(40)

The quadratic separation is therefore

l2=(nL+QL)2lf2
(S22)

. This finite System Frame realization supplies the realized separation functional from which the local continuum metric is reconstructed. The metric is not introduced independently of the mapped count structure.

For the local continuum description, let xμ denote realized System Frame coordinates and dxμ their differential displacements. The differential interval dl has the standard quadratic form

dl2=gμνsys(x)dxμdxν,
(41)

where

gμνsys(x)
(S23)

denotes the emergent System Frame metric tensor evaluated at x.

To reconstruct the metric covariantly from separation data, let l(x,x') denote the realized two-point separation between neighboring System Frame points x and x'. In a normal neighborhood where l2(x,x') is smooth and admits a nondegenerate continuum limit, the mixed-endpoint coincidence relation gives [34]

gμνsys(x)=-12limx'→x∂2l2(x,x')∂xμ∂xν.
(42)

The minus sign follows from differentiation with respect to opposite endpoints of the squared separation. Equivalently, with Synge’s world function

σ(x,x')=12l2(x,x')
(S24)

, the standard coincidence identity is [σ;μν']=-gμν [34]. Within MQ, the distinction lies in the origin of l2(x,x'): it is the realized separation functional supplied by mapped count structure rather than a separation assumed from a primitive metric field.

The continuum coordinates arise from discrete count increments through xμ=nμlf, where nμ is the realized count coordinate corresponding to xμ. Since lf is constant,

∂∂xμ=1lf∂∂nμ.
(43)

Applying Eq. (43) to both coordinate arguments of Eq. (42) cancels the explicit

lf2
(S25)

in the realized quadratic separation. Let n=(nμ) and n'=(n'μ) denote the neighboring count-coordinate arguments corresponding to x and x'. The scalar count quantities evaluated between them are written nL(n,n') and QL(n,n'). The count-resolved metric reconstruction is then

gμνsys(x)=-12limn'→n∂2(nL(n,n')+QL(n,n'))2∂nμ∂nν.
(44)

Equation (44) reconstructs the local System Frame metric from the second-order response of the realized separation functional. It does not by itself select a global metric, signature, or curvature field. It specifies the local continuum metric when coarse-graining of the realized count structure admits a smooth, nondegenerate continuum representation.

The defining count relation follows directly from the right-triangle construction,

1+nL2=(nL+QL)2,
(45)

where the scalar count separation satisfies the exact mapping-consistent identity established in Sec. 3.1. Equations (40–45) therefore connect the finite right-triangle count relation to its local continuum metric representation. The metric is not a primitive structure of the Internal Frame, but a System Frame continuum object reconstructed from the realized separation functional through the Frames mapping.

Equations (40–45) establish metric emergence only. They are a kinematic reconstruction and do not, by themselves, supply a gravitational field equation, matter source law, geodesic dynamics, or photon propagation. In the broader MQ construction summarized in Ref. [1], the reconstructed metric enters the MQ Action, whose stationary configurations are governed by the MQ variational principle. Metric variation defines the stress-energy response, while the Levi-Civita connection, curvature, diffeomorphism compatible variation, and contracted Bianchi identity provide the corresponding geometric dynamics and covariant conservation. The metric therefore acquires dynamics through this separately derived variational construction, not merely because it can be reconstructed from the realized separation functional.

The relativistic completion in Ref. [1] extends the realized spatial geometry, under its stated causal order and regularity conditions, to a time-oriented Lorentzian System Frame geometry supporting locally Lorentz covariant matter dynamics and matter geometry feedback. These downstream results do not make Eqs. (40–45) dynamical. They establish that the metric reconstructed here is the geometric input to the separate MQ dynamical completion.

For the present dark matter analysis, the galactic realization law developed in Sec. 4.2 is the nonrelativistic weak-field result used for orbital dynamics. Rotation curve agreement alone does not derive the complete relativistic metric solution for an extended galaxy. The independent lensing problem in Sec. 5.3 instead requires solving the System Frame field dynamics for the specified baryonic source and propagating null geodesics through that solution. Metric emergence, metric dynamics, and a system-specific relativistic solution are therefore separate requirements.

3.6 Non-Lorentz length contraction
Using the right-triangle count construction introduced in Fig. 1, MQ predicts a non-Lorentz length contraction arising from discrete Internal Frame count structure and its realization in the System Frame. Unlike Lorentz contraction, whose magnitude depends on relative velocity between inertial observers, the MQ effect is the Informativity differential: a count-dependent geometric manifestation of normalization under the Frames mapping. It is therefore not a relativistic velocity effect, but a manifestation of the mapping-induced deformation quantified by QL.

Realized separation is obtained from the mapping-consistent combination nL+QL rather than by direct projection of the Internal Frame count nL. Beginning with the count identity established in Eq. (45),

12+nL2=(nL+QL)2,
(46)

the right-triangle relation is

QL+nL=1+nL2,
(47)
QL2+2QLnL+nL2=1+nL2,
(48)
QL2+2QLnL=1.
(49)

The exact relation in Eq. (49) has the upper-count behavior

limnL→∞QL2+2QLnL=1.
(50)

Since QL→0 as nL→∞, the quadratic contribution vanishes asymptotically, giving

limnL→∞2QLnL=1.
(51)

Thus, the normalization factor 2QLnL approaches unity, while the Informativity differential, 1-2QLnL, approaches zero. The classical large-count limit is recovered, whereas finite-count realization retains a nonzero mapping-induced contraction.

This distinguishes the realization regimes established in Sec. 3.4. Electromagnetic interactions are evaluated at the fixed interaction scale defined by the electromagnetic demarcation, where the Informativity differential remains finite. Gravitational realization occupies the curvature-resolved regime, in which increasing count separation drives the same mapping-induced deformation toward its upper-count limit.

Two applications are the reduced Planck constant and gravitational coupling. At the electromagnetic demarcation, the finite-count reduced Planck constant ℏ contains the non-Lorentz normalization,

ℏ=psilfQLnL,
(52)

while the finite-count gravitational coupling G follows from the same count geometry,

G=QLrc3psi=QLnLlfc3psi.
(53)

Equation (53) shows that finite-count G is constructed from the momentum realization psi, the mapping-induced residual QL, the count separation nL, and the fundamental measures. Its numerical value therefore depends on the count regime at which the expression is evaluated.

At the electromagnetic demarcation, the finite-count evaluation is

G=2QLnLc3tfmf|nL=84.60055=6.67384483(17)10-11m3kg-1s-2,
(54)

This evaluates to

6.67384483(17)×10-11m3kg-1s-2
(S26)

, where italic G denotes the gravitational coupling evaluated at this finite interaction count.

At the upper-count limit, QL→0, nL→∞, and 2QLnL→1, giving the upright upper-count gravitational constant,

G=limnLr→∞2QLnLc3tfmf=c3tfmf=6.67407794(10)10-11m3kg-1s-2.
(55)

The upper-count value is

6.67407794(10)×10-11m3kg-1s-2
(S27)

. Equations (54) and (55) are finite-count and upper-count evaluations of the same MQ gravitational coupling. Relative to the 2022 CODATA recommendation

6.67430(15)×10-11m3kg-1s-2
(S28)

[7], they are distinct MQ realization values rather than separate constants recognized by conventional gravitational metrology.

The historical CODATA sequence provides a retrospective numerical comparison. The 2010 recommendation differs from the finite-count value of Eq. (54) by 0.006σ, the 2014 recommendation differs from the upper-count value of Eq. (55) by 0.007σ [35–36], and Ref. [15] obtains a Planck-coupled evaluation close to the 2018 value retained in the 2022 adjustment [7,37], with a reported correspondence of approximately 0.08σ. Because CODATA adjustments combine heterogeneous measurements and were not constructed according to MQ realization regimes, Appx. E classifies these historical matches as one retrospective concordance family rather than three statistically independent confirmations. Measurements of G also remain unusually dispersed relative to other fundamental constants [38]. The CODATA sequence therefore does not establish that its changes were caused by the Informativity differential.

A temporally later external test became available in 2026 from the NIST redetermination using the BIPM torsion balance [39]. The measurements themselves were acquired in 2024, substantially after publication of the MQ Informativity differential and its realization-dependent gravitational predictions. Schlamminger et al. obtained [value], approximately [value] below the 2022 CODATA recommendation and [value] from the finite-count prediction of Eq. (54), using the combined quoted uncertainties. The upright notation identifies the macroscopic gravitational determination, whereas comparison with the finite-count Gof Eq. (54) tests the previously published MQ realization prediction.

The significance of this experiment extends beyond numerical agreement in the combined value. The BIPM torsion balance determines the same gravitational interaction by two physically distinct closure methods. In free-deflection operation, all relevant components are electrically grounded and the gravitational torque is inferred from the mechanical response of the torsion balance. In electrostatic-servo operation, the pendulum electrode is grounded while driven electrodes supply an electrostatic counter-torque that balances the gravitational torque. These are precisely the contrasting experimental conditions relevant to the MQ distinction between curvature-resolved macroscopic gravitational realization and gravitational determination closed through a fixed electromagnetic interaction scale.

The realization-dependent distinction predates the NIST/BIPM measurements and was not introduced in response to their outcome. MQ predicts that a determination whose torque closure depends explicitly on the electromagnetic interaction should exhibit the finite-count realization more strongly than a free mechanical deflection determination, shifting the inferred gravitational coupling toward the finite-count value. Schlamminger et al. subsequently found this ordering independently: for both copper and sapphire source masses, the electrostatic-servo determination was lower than the corresponding free-deflection determination. Their paper explicitly identifies this repeated ordering across the two source-mass materials.

The experimental group subjected the discrepancy to extensive investigation. The analysis considered source- and test-mass geometry, pressure-dependent effects, autocollimator nonlinearity, capacitance and voltage calibration, mass integration, and other method- and configuration-dependent systematic contributions [39]. The distinction between the two operating modes was itself a longstanding motivation of the BIPM program: the apparatus was designed to determine G by both free deflection and electrostatic compensation in part to investigate the persistent dispersion among precision measurements of G. The 2026 analysis did not identify a conventional experimental mechanism that converts the observed mode ordering into the MQ prediction.

The NIST/BIPM experiment therefore provides a direct experimental test of a prediction published by MQ before these measurements were performed. The experiment found the predicted direction of the realization-dependent shift and, within the stated comparison, the predicted gravitational-coupling scale.

As this effect is directly a consequence of the MQ proposed Frames mapping and those properties recognized by the System and Internal Frames, the result is a definitive finding in support of the MQ approach. Notably, the result carries substantially greater evidentiary weight than the retrospective correspondence with historical CODATA recommendations because the MQ prediction preceded the measurements used to test it. The NIST/BIPM authors did not analyze their result within the MQ framework or cite the earlier MQ prediction. Consequently, the experiment should not be described as an endorsement of MQ by those authors. Rather, their independently obtained measurements provide direct experimental evidence consistent with—and quantitatively testing—the prior MQ realization prediction.

With this realization structure established, the conventional Planck unit expressions follow from the MQ count-geometric construction ([15], Appx. BN),

lf=(ħGc3)1/2,
(56)
tf=(ħGc5)1/2,
(57)
mf=(ħcG)1/2.
(58)

These equations are terminal conventional forms, not starting definitions of lf, mf, or tf. Ref. [15] first resolves the gravitational relation and fundamental measures from the right-triangle count construction, invariant coefficient θsi, and fundamental count relations; those relations are then combined in Appx. BN of Ref. [15] to recover Eqs. (56–58). The G and ħ appearing in these Planck-coupled expressions can then by generalized from their finite-realization typography thus giving us Planck’s expressions without assuming the dimensional expressions at the outset.

3.7 An Example of the Frames mapping via the fine structure constant
The preceding sections established the count geometry, electromagnetic demarcation, and Informativity differential required by the Frames mapping. This subsection follows those elements through successive realization stages to the experimentally accessible inverse fine structure constant. Because α is also used upstream to determine θsi in Sec. 3.8, the calculation is an internal closure test of the MQ realization structure rather than an independent determination of α.

As established in Sec. 3.4, electromagnetic realization occurs at the fixed interaction scale nL. One-half of the demarcation is approximately 42.3, whose nearest whole-unit realization is the System Frame count nθ=42. MQ associates this count with the fundamental inverse fine structure form through 42θsi. The realized count nθ=42 and the quantity 42θsi are distinct: the former is the selected whole-unit count, whereas the latter is the corresponding inverse fine structure quantity before subsequent realization stages.

The invariant coefficient θsi was introduced in Sec. 3.3 and is treated here as established. Section 3.8 determines its numerical value by solving the inverse MQ realization relation using the measured fine structure constant. Here the direction is reversed: the established θsi is propagated through the whole-count relation, Frames mapping, metric differential, and finite-count Informativity differential to recover the classically realized inverse fine structure constant. Comparison with the same CODATA fine structure input used in Sec. 3.8 therefore tests algebraic closure rather than providing an independent prediction.

Equation (59) defines the fundamental Internal Frame form,

αf-1=42θsi.
(59)

Application of the Frames mapping produces the continuous System Frame representation. The corresponding fine structure metric differential is

Qnθ=42θsi-⌈42θsi⌋,
(60)

where the mixed floor-ceiling brackets denote nearest whole-unit realization. The metric differential is the representation-level mismatch required when the discrete count relation is expressed continuously in the System Frame. It is distinct from the physical contraction described by the Informativity differential.

Adding

Qnθ
(S29)

to the fundamental form gives the Planck-like inverse fine structure form,

αp-1=αf-1+Qnθ=42θsi+(42θsi-⌈42θsi⌋),
(61)
αp-1=84θsi-⌊42θsi⌉.
(62)

At this stage the construction includes the discrete-to-continuous metric differential but not the finite-count normalization established in Sec. 3.6.

At the electromagnetic demarcation, the charge-coupling relation gives

nL=276/θsi,
(63)

while the residual count product from the right-triangle construction in Sec. 3.1 is

QLnL=nL((1+nL2)1/2-nL).
(64)

Together, Eqs. (63) and (64) determine 2QLnL at the fixed interaction scale.

Let Δ(int-sys) denote the finite-count correction carrying the Planck-like inverse fine structure form to its classically realized value under the Informativity differential. Equation (65) defines

Δ(int-sys)=αp-1(1-2QLnL).
(65)

Subtracting this correction gives

αc-1=αp-1-Δ(int-sys),
(66)

which reduces to

αc-1=αp-1-αp-1(1-2QLnL)=2QLnLαp-1.
(67)

Substitution of Eq. (62) gives the classically realized form,

αc-1=2QLnL(84θsi-⌈42θsi⌋).
(68)

The resulting

αc-1
(S30)

is the MQ quantity directly comparable with a laboratory determination of α-1.

Table 4. Fundamental, Planck-like and classically realized forms of the inverse fine structure constant

Measure

Fundamental αf-1

Planck-like αp-1

Classical αc-1

Morel et al. [40]

—

—

137.035999206(11)
(11)

MQ realization

137.020392764(84)
(84)
137.0407856(30)
(30)
137.0359992(30)
(30)

Morel et al. measured the inverse fine structure constant by rubidium atom-recoil matter-wave interferometry, obtaining α-1=137.035999206(11) with a relative accuracy of 81 parts per trillion [40]. The fundamental and Planck-like entries are intermediate MQ realization stages rather than corresponding Morel observables. The classical MQ value differs from the Morel result by approximately 0.19σ when the quoted standard uncertainties are combined. This agreement is a closure benchmark, not an independent experimental confirmation, because the Morel result contributes to the fine structure constant measurement family incorporated into the CODATA adjustment used upstream.

The like-for-like comparison with the 2022 CODATA recommendation used in Sec. 3.8, α-1=137.035999177(21) [7], gives a difference of approximately 0.63σ from the classically realized MQ value 137.035999200(30). Because that CODATA value is used to determine θsi, this is likewise an internal concordance result rather than an independent prediction test.

Table 4 therefore records successive MQ realization stages, not three separate experimental predictions. Equations (59–68) provide a worked example in which the discrete count relation, metric differential, and Informativity differential combine under the Frames mapping to recover the experimentally accessible inverse fine structure form.

3.8 Calculating the invariant coefficient as a function of the fine structure constant

Section 3.7 demonstrated the forward MQ realization sequence by propagating an established θsi through the Frames mapping to the classically realized inverse fine structure constant

αc-1
(S31)

. The present subsection performs the inverse operation: the 2022 CODATA recommended inverse fine structure constant, α-1=137.035999177(21) [7], is identified with the corresponding classically realized MQ quantity

αc-1
(S32)

and used to determine θsi. Equation (76) is therefore a parameter determination from the measured fine structure constant, not an independent prediction of it.

The inversion begins with the right-triangle count construction established in Sec. 3.1 and illustrated in Fig. 1. The mapping-induced residual count contribution is

QL=(1+nL2)1/2-nL.
(69)

Multiplication by θsi gives

QLθsi=(θsi2+(θsinL)2)1/2-nLθsi,
(70)

and division by θsi gives

QL=(θsi2+(nLθsi)2)1/2-nLθsiθsi.
(71)

The inverse fine structure relation is the classically realized MQ form established in Sec. 3.7. It incorporates the representation-level metric differential associated with the Frames mapping and the finite-count normalization associated with the Informativity differential,

1αc=2QLnL(84θsi-⌈42θsi⌋),
(72)

where the mixed floor-ceiling brackets denote nearest whole-unit realization.

No independently measured interaction count enters the inversion. The electromagnetic demarcation established in Sec. 3.4 fixes nLθsi=276, equivalently

nL=276/θsi
(S33)

. Thus, the measured

αc-1
(S34)

is the sole experimental input to the numerical determination, while nL, QL, and their relation are supplied by the MQ construction.

For the following rearrangement, αc denotes the classically realized fine structure constant and

αc-1
(S35)

its inverse. Rearranging Eq. (72) gives

2αc=1nLθsiθsi1QL(84θsi-⌈42θsi⌋),
(73)
2αc=1nLθsi1QLθsi(84θsi-⌈42θsi⌋).
(74)

Substitution of Eq. (71) into Eq. (74), with nLθsi=276, gives

2αc=1nLθsiθsi(θsi2+(nLθsi)2)1/2-nLθsiθsi(84θsi-⌈42θsi⌋).
(75)

The nearest-whole-unit operation makes Eq. (75) piecewise across its admissible count realizations. The applicable branch is fixed by the electromagnetic demarcation and nearest-whole-unit realization established in Secs. 3.4 and 3.7 rather than selected by fitting the resulting value. Solving Eq. (75) on that branch with the 2022 CODATA value α-1=137.035999177(21) [7] gives

θsi=3.262390305(36).
(76)

The quoted uncertainty propagates from the CODATA fine structure input through the nonlinear inversion. No independent determination of nL is required because the charge-coupling relation fixes nLθsi=276.

Equation (76) establishes the numerical invariant coefficient used in subsequent MQ calculations. It is complementary to the forward construction in Sec. 3.7: Sec. 3.8 uses the measured fine structure constant to determine θsi, whereas Sec. 3.7 propagates that coefficient through the realization sequence to test internal algebraic closure. Reconstructing the fine structure constant from the same upstream CODATA source is therefore an internal closure test rather than an independent prediction.

3.9 Referenceability-domain fractions of the expanding frame

The expanding frame domain relations provide a cross-scale consequence of the same upstream MQ invariant coefficient used elsewhere. They are mapping-defined referenceability fractions, not independently conserved ΛCDM fluids. Their numerical proximity to Planck density fractions is therefore a correspondence benchmark. A dynamical cosmological test additionally requires the background and perturbation evolution developed in Sec. 5.8.

Ωf=2θsi2,
(77)
Ωdk=θsi2-2θsi2+2,
(78)
Ωobs=4θsi2+2,
(79)
Ωvis=2θsi(θsi2+2),
(80)
Ωuobs=Ωobs-Ωvis=4θsi-2θsi(θsi2+2).
(81)

Here Ωf is the fundamental domain quantity, Ωdk the dark domain, Ωobs the observable domain, Ωvis the visible domain, and Ωuobs the unobserved domain. The realized referenceability partition satisfies Ωobs=Ωvis+Ωuobs and Ωdk+Ωobs=1, as established in Appx. C. The fundamental domain Ωf is not an additional additive component of this partition; it constrains the referenceability geometry.

Once θsi has been determined in Sec. 3.8, Eqs. (77–81) follow algebraically without fitting matter-density, energy-density, or astrophysical population data. Their uncertainties therefore propagate the source uncertainty in θsi through these closed-form relations.

The MQ fractions may be compared numerically with fractional density parameters inferred under ΛCDM, but the quantities have different physical definitions. In the standard cosmological analysis, baryonic, cold dark matter, total matter, and dark energy fractions are inferred jointly within a dynamical background and perturbation model from the CMB spectra and other cosmological data [11]. Equations (77–81) instead describe referenceability geometry. Table 5 therefore reports numerical correspondence rather than dynamical equivalence.

Table 5. Numerical comparison of MQ referenceability domain fractions with Planck 2018 ΛCDM density fractions density fractions and associated MQ mass equivalent allocations

Phenomenon / MQ Term

Domain

ΛCDM % a

(MQ) %

(MQ) kg

Dark energy ΩΛ / Dark Ωdk

68.34±0.84%

68.36241612(59)

2.0857 1054

Dark Matter Ωc / Unobserved Ωuobs

26.73±0.82%

26.78874899(44)

8.1733 1053

Ωc+Ωb/ Observable Ωobs

31.66±0.84%

31.63758388(59)

9.6526 1053

Baryonic Ωb / Visible Ωvis

4.97±0.82%

4.848833489(14)

1.4794 1053

- / Fundamental Ωf

-

18.79135772(41)

5.7333 1053

The stronger test is dynamical: the MQ field equations must generate a homogeneous cosmological background and perturbations that reproduce CMB, BAO, CMB-lensing, and structure-growth observables without inserting the corresponding ΛCDM density fractions as MQ inputs. Section 5.8 establishes the field/conservation linearization and regular source formulation needed for that calculation but does not yet provide the numerical Einstein–Boltzmann evolution. The Table 5 fractions therefore must not be assigned directly to independently conserved System Frame stress-energy components.

Conversion of the percentages to the absolute mass-equivalent allocations shown in Table 5 additionally requires the elapsed universe age. The realized expanding frame construction uses the MQ mass accretion rate derived in Sec. 3.10,

Macr=θsi3ρsilf=7.00887947539(16)×1036kgs-1,
(82)

where psi is the momentum realization introduced in Sec. 3.3. Equation (82) is stated here because it supplies the Table 5 mass column; its derivation follows in Sec. 3.10. For a specified universe age, integration of this realized rate gives the accumulated total mass-equivalent content Mtot, to which the dimensionless fractions are applied. Thus θsi fixes the domain fractions, the expanding frame construction fixes Macr, elapsed age fixes Mtot, and the fractions fix the associated mass-equivalent allocations.

The fundamental domain Ωf has no direct ΛCDM counterpart and is not an additional matter or energy component.

3.10 Count allocation before realization and mass accretion after realization

This subsection records the allocation/accretion consequence of the expanding frame count relations. Before System Frame realization, the relations describe count allocation; after realization, the same algebraic structure acquires a temporal mass-equivalent interpretation. These relations are retained for cosmological provenance and are not required to evaluate galaxy rotation curves. The galactic dependency proceeds instead through the frame-resolved expansion rate developed in Sec. 3.12.

The dimensionally explicit fundamental relation uses the momentum realization psi,

lfmf=2ρsitf.
(83)

Equation (83) is a relation among the fundamental measures and their momentum realization. The universe-scale length, mass, and time counts are not obtained by multiplying independent count factors into both sides of Eq. (83). Instead, the expanding frame construction fixes separate count relations for spatial extent and accumulated mass-equivalent content.

For the universe-age count

nTU
(S36)

, the realized diameter count is

nLU=2θsinTU,
(84)

while the total mass-equivalent count is

nMU=θsi32nTU,
(85)

where

nLU
(S37)

is the universe-scale count of lf,

nMU
(S38)

is the count of mf corresponding to the accumulated total mass-equivalent content, and

nTU
(S39)

is the elapsed count of tf. Their ratio is therefore

nMUnLU=θsi24,
(86)

and the dimensionless diameter-age relation may be written

nLU2nTU=DU2AUc=θsi,
(87)

with

DU=nLUlf
(S40)

,

AU=nTUtf
(S41)

, and

c=lf/tf
(S42)

. Equation (87) is dimensionless throughout and supplies the count identity used later in Sec. 3.12.

The elapsed universe-age count is

nTu=AU1tf,
(88)

where nTᵤ denotes the elapsed count of fundamental time units tf. Restoring the fundamental length-to-time conversion gives the physical radius,

RU=AUθsi(lftf),
(89)

and the corresponding fundamental domain mass allocation is

Mf=AUθsi(mftf),
(90)

where Mf denotes the allocation associated with the fundamental domain and is distinct from the fundamental mass measure mf. Prior to System Frame realization, Eq. (90) specifies count-based allocation rather than physical mass production through elapsed time.

The domain relations required for the global allocation are

Ωobs=4θsi2+2,
(91)
Ωf=2θsi2.
(92)
Ωf=2θsi2=4θsi2+2/(2θsi2θsi2+2),
(93)
Ωf=4θsi2+2/(2θsi2θsi2+2)=4θsi2+2/(2θsi2+4-4θsi2+2),
(94)
Ωf=4θsi2+2/(2-4θsi2+2)=Ωobs2-Ωobs,
(95)
2-ΩobsΩobs=1Ωf,
(96)
Ωf(2Ωobs-1)=1=Ωtot.
(97)

Equation (97) establishes the normalized relation between the fundamental and observable domains. The corresponding mass allocations are related to the normalized fractions by

MfMtot=Ωf,MobsMtot=Ωobs,
(98)

rather than by substituting dimensional masses directly for the dimensionless Ω quantities.

Using the universe-scale mass count from Eq. (85), the accumulated total mass-equivalent content is

Mtot=nMUmf=nTUmfθsi32.
(99)

Equation (99) shows that after System Frame realization, accumulated mass-equivalent content is proportional to elapsed fundamental time count. Before realization, the same relation specifies count allocation. Equations (90), (92), and (99) are mutually consistent because

Mf/Mtot=2/θsi2=Ωf.
(S43)

Let

T=nTUtf=AU
(S44)

denote the realized elapsed duration. Dividing the accumulated mass equivalent by T gives

Macr=MtotT=nTumfθsi321nTutf=θsi32mftf=θsi32unitsmfpertf.
(100)

The final form in Eq. (100) is equivalently a rate of

θsi3/2
(S45)

fundamental mass units mf per fundamental time unit tf. Using the dimensionally explicit momentum relation

psi=lfmf/(2tf)
(S46)

yields

Macr=θsi32(mftf)=θsi32(2ρsilf)=θsi3ρsilf=7.00887951036kgs-1.
(101)

Equation (101) defines the MQ mass accretion rate Macr for the realized expansionary epoch. It is distinct from the fundamental mass mf, the pre-realization fundamental domain allocation Mf, and accumulated total mass Mtot. It is not a conventional astrophysical accretion rate such as gas infall or halo growth. It is the global MQ rate associated with the expanding frame count relation.

For the adopted universe age used in this work, the count relations give approximately

nTU=8.0763×1060,nLU=5.2696×1061,nMU=1.4021×1062,
(102)

with

nLU2nTU=θsi=3.262390305,nMUnTU=θsi32=17.36112075.
(103)

The final numerical line is included only to make the count interpretation explicit. The defining relations are Eqs. (84) and (85).

For the dark matter application, Sec. 3.9 supplies the global referenceability fractions, while Eqs. (98-101) provide the corresponding allocation/accretion accounting. Local gravitational evolution may subsequently redistribute realized mass-equivalent content into galaxies, clusters, and other structures, but such structure formation is distinct from the global MQ mass accretion law. The observationally resolved expansion rate representations relevant to the galactic construction are developed separately in Sec. 3.12.

3.11 CMB temperature from the MQ count-geometric construction

The CMB-temperature calculation is retained as a cross-scale consistency check rather than as part of the galaxy prediction. It propagates the expanding frame construction to a present radiation density and blackbody temperature without using the observed CMB monopole temperature as an input. The calculation does, however, inherit the adopted Planck cosmic age and is therefore not statistically independent of all CMB information. No CMB-spectrum validation is claimed.

The derivation begins with the minimum three-dimensional count configuration at which a globally referenceable realization first becomes admissible. Let na and nb denote orthogonal unit-count components and nbn the resulting two-dimensional diagonal count entering the third-dimensional extension. These labels specify components of this local count construction and introduce no additional MQ count primitives.

na2+nbn2=na2+(na2+nb2)2=
(104)
na2+na2+nb2=12+12+12=3.
(105)

Using this minimum count with the previously established count relations gives the kinematically resolved expansion quantity v' associated with this referenceability construction. The universe-scale counts

nLU
(S47)

,

nMU
(S48)

, and

nTU
(S49)

denote counts of the fundamental length, mass, and time measures lf, mf, and tf, respectively. The time count

nTU
(S50)

is the universe age count established in Sec. 3.10, while

nLU'
(S51)

denotes the boundary-resolved length count entering the minimum referenceability construction. The quantity v' acquires a dimensionful velocity interpretation only upon realization.

Dimensional consistency requires the momentum realization psi introduced in Sec. 3.3. Thus,

nLulfnTutf=2ρsinMumf,
(106)

and for the boundary-resolved length count,

nLu'lfnTutf=2ρsinMumf.
(107)

Let

MU=nMUmf
(S52)

denote the universe-scale mass-equivalent content entering this reduction, VU the corresponding realized volume, and ρc the critical density, so that MU=VUρc. The expansion quantity can then be written

v'=nLu'lfnTutf=2ρsinMumf=2ρsiMU=2ρsiVUρc=2ρsi(4/3)πc3θsi3Aqe3ρc.
(108)

For the critical density reduction, let Hf denote the expansion-rate quantity entering the conventional relation

ρc=3Hf2/(8πG)
(S53)

. Here G is the upright upper-count gravitational constant. Using that relation together with

psiG/c3=lf/2
(S54)

gives

v'=3ρsi2πc3θsi3Aqe31ρc=3ρsi2πc3θsi3Aqe38πG3Hf2=4ρsic3θsi3Aqe3Hf2lfc32ρsi,
(109)
v'=1Hf22lfθsi3Aqe3=Aqe22lfθsi3Aqe3=2lfθsi3Aqe.
(110)

The resulting scale is integrated across the quantum epoch count construction. Within MQ, the quantum epoch is a pre-referenceable regime rather than an interval of physical System Frame time. Let Rqe denote its radial scale and Aqe the corresponding SI-equivalent temporal scale.

Rqe=2lfθsi3∫(1Aqe)dt=2lfθsi3ln⁡(Aqe).
(111)
Aqe=e3θsi3/2=1.1465266(11)1013s=363,312.35y.
(112)

The numerical value of Aqe therefore represents the SI-equivalent extent of the pre-referenceable count construction, not elapsed time experienced within a realized universe.

The solution is next propagated through the expansionary epoch, which begins after establishment of a valid System Frame. Let RU denote the realized universe radius and Aee the corresponding elapsed System Frame expansion time,

RU=2lfθsi3∫(1AU)dt=2lfθsi3ln⁡(AU),
(113)
Aee=Aqe(2θsi)1/3=2.14242731013s=678,894.26y.
(114)

The interpretations of Aqe and Aee are therefore distinct: Aqe is an SI-equivalent encoding of the quantum epoch count scale, whereas Aee is realized elapsed time in the expansionary System Frame.

The count allocation fixed by the quantum epoch construction is next combined with the allocation/ accretion relation derived in Sec. 3.10. Prior to realization this is a count-based allocation rather than physical mass-energy evolving in spacetime. Let MCMB denote the mass-energy allocation obtained when that count structure is carried into its realized CMB representation.

MCMB=nTumfθsi32=1.50161050kg.
(115)

Let ρCMB denote the corresponding present CMB radiation energy density. Then

ρCMB=MCMBc2VU=4.172310-14Jm-3.
(116)

The conversion from radiation energy density to temperature uses the standard blackbody relation

ρCMB=aTCMB4
(S55)

, where a is the radiation constant. Let σ denote the Stefan-Boltzmann constant and kB the Boltzmann constant. The radiation constant is

a=4σc=π2kB415c3ℏ3=7.56573325010-16Jm-3K-4,
(117)

where the reduced Planck constant ℏ retains the italic finite electromagnetic realization notation appropriate to this blackbody construction.

The resulting CMB monopole temperature is

TCMB=(ρCMBa)1/4=2.72509K.
(118)

The temperature difference between the MQ value and the Fixsen determination is approximately 0.48σ using the uncertainties shown in Table 6. No analogous significance is assigned to the radiation density row because that benchmark is derived from the same CMB temperature. The Aqe comparison is likewise qualitative because the MQ quantity is an SI-equivalent pre-referenceable scale whereas the observational benchmark is a model-inferred physical age near photon decoupling.

Table 6. Comparison of MQ predicted CMB values with measurements.

ACMB

ρCMB

TCMB

MQ

363,312.34127(22) y

4.1723 10-14 Jm-3

2.72509 ± 0.00341 K

Measurements

379,000 y b

4.17 10-14 Jm-3 c

2.72548 ± 0.00057 K a

a Consolidated CMB-temperature determination from Fixsen [41].
b Model-inferred physical age near photon decoupling; this is not physically identical to the MQ SI-equivalent quantum epoch scale Aqe.
c Radiation-density benchmark derived from the measured CMB temperature; no independent uncertainty is used for a separate concordance statistic.

The principal external input to the present-day MQ temperature calculation is the adopted universe age. This work uses the Planck 2018 base-ΛCDM value near 13.797 Gyr [11] to determine the present universe volume. Because that age is inferred largely from CMB anisotropy data within a cosmological model, the resulting temperature agreement is not statistically independent of all CMB information. Neither the observed CMB monopole temperature nor its present radiation energy density is fitted or supplied to Eqs. (115–118).

The calculation is therefore retained as an internal cross-scale closure test. No CODATA concordance statistic applies to Eqs. (104–118), and no CMB spectrum, acoustic peak, or perturbation-level validation follows from this subsection. Those remain part of the dynamical cosmological program identified in Sec. 5.8.

3.12 Expansion of the universe in terms of the Hubble constant

The preceding sections established the MQ geometric description of universal expansion and distinguished its dimensionless scalar radial rate from its dimensionful local rate representation. This subsection collects the relations required to connect that cosmological construction to the galactic realization developed in Sec. 4. A complete derivation and observational comparison of the MQ Hubble tension construction was published in Ref. [42]. The downstream quantity required here is the Internal Frame local expansion rate Hint, which enters the independently specified MQ acceleration scale in Sec. 4.2.

MQ distinguishes the mapping-defined scalar rate of universal expansion from its observable local rate representation. Upright H denotes the dimensionless scalar radial rate, while italic H denotes its dimensionful realization in conventional units such as km s-1 Mpc-1. The same underlying expansion geometry may be resolved relative to the System Frame or Internal Frame, giving the local rates Hsys and Hint, respectively.

Ref. [42] interprets the Hubble tension as arising from these two realizations of a common universal expansion geometry: CMB-based cosmological inference is associated with the continuous System Frame realization, whereas Cepheid-calibrated distance ladder inference is associated with the discrete Internal Frame realization. This is an MQ interpretation of the operational distinction between the measurements; conventional cosmology does not classify the two methods in terms of MQ frames. For the quantitative comparison retained here, the benchmark values used in Ref. [42] are the SPT-3G and Pantheon+/SH0ES determinations, 68.3±1.5 and 73.5±1.1 km s-1 Mpc-1, respectively [43–44].

MQ represents the separation between the two frame realizations through the metric differential,

ΔH=θsi-π,
(119)

where ΔH is the expansion-specific mapping-induced normalization offset between the System Frame and Internal Frame realizations. It is generated by the Frames mapping, is not an independent dynamical parameter, and is distinct from the finite-count Informativity differential developed in Sec. 3.6.

The dimensionless universal scalar radial rate follows from the normalized expanding frame relation established in Sec. 3.10. Using

c=lf/tf
(S56)

,

H=12cnLulfnTutf=12lftf=θsi,
(120)

Equation (120) establishes the mapping-defined scalar H=θsi. It is distinct from both the dimensional momentum realization

psi=mflf/(2tf)
(S57)

and the dimensionful local expansion rate H.

The System Frame scalar argument is displaced below H by ΔH, while the Internal Frame argument is displaced above it by the same amount. Their frame-resolved local rates are

Hsys=f(H-ΔH),
(121)
Hint=f(H+ΔH),
(122)

where f denotes the fixed linear mapping from a dimensionless scalar argument to its observable local-rate representation [42]; it is not an adjustable fitting function.

Substituting H=θsi and ΔH=θsi-π into the System Frame relation gives

Hsys=f(H-ΔH)=f(θsi-(θsi-π))=f(π),
(123)

and linearity of f gives the corresponding Internal Frame form,

Hint-f(ΔH)=Hsys+f(ΔH),
(124)
Hint=Hsys+f(2ΔH)=f(π)+f(2θsi-2π)=f(2θsi-π),
(125)

so the frame-resolved scalar arguments are π and 2θsi-π. Equations (123–125) therefore describe symmetric frame resolution about the common scalar θsi, not two independent background expansion laws.

To compare these scalar arguments with astronomical measurements, the fixed normalization maps a dimensionless scalar argument to a conventional local rate. Italic H denotes this general dimensionful local-rate representation, whereas upright H remains the underlying dimensionless scalar. Using the present universe age relation

AU=DU/(2θsic)
(S58)

,

H=1AUkmMpc=2θsiDUkmMpc.
(126)

Equation (126) changes representation and units; it does not redefine H=θsi. Substitution of the System Frame and Internal Frame scalar arguments gives

Hsys=2π2θsiAUkm/Mpcs/y109=68.259(19)kms-1Mpc-1.
(127)
Hint=2(2θsi-π)2θsiAUkm/Mpcs/y109=73.508(21)kms-1Mpc-1.
(128)

These are frame-resolved local-rate realizations of the same universal expansion geometry. The first is numerically comparable with the SPT-3G benchmark 68.3±1.5 km s-1 Mpc-1 [43], and the second with the Pantheon+/SH0ES benchmark 73.5±1.1 km s-1 Mpc-1 [44]. The corresponding differences are approximately 0.027σ and 0.007σ when the quoted uncertainties are combined.

Neither benchmark experimentally identifies a System Frame or Internal Frame. MQ supplies the hypothesis that the separation between the observational constructions corresponds to different realizations under the Frames mapping.

Ref. [42] additionally corresponds to a 3.35σ right-tail joint concordance. The present paper does not treat that statistic as an independent combined significance because Hsys and Hint share θsi, the cosmic age normalization, and the same MQ frame construction. The individual residuals above are retained only as numerical compatibility measures.

For the present dark matter analysis, the relevant consequence is Hint. Section 4.2 uses this Internal Frame local rate, together with the defined speed of light c and invariant coefficient θsi, to construct the fixed MQ acceleration scale. The lower Hsys branch is not used as the galactic input. Thus, the cosmological construction supplies a specifically identified upstream quantity rather than a quantity fitted to galaxy data.

4. RESULTS

This section tests the MQ galactic realization framework against two observationally distinct rotation curve datasets. The prediction is constructed directly from the adopted observational baryonic inputs and the fixed present-epoch galactic realization law, then evaluated against the Milky Way reconstruction and the SPARC database of 175 galaxies comprising 3,366 radial measurements [4]. Because the datasets arise from different observational procedures and baryonic mass reconstructions, their statistics are reported separately. The Milky Way calculation uses the observationally independent Eilers et al. circular velocity target [2] and a McMillan baryonic source [3] constructed without an RAR or MOND prescription; the McMillan NFW halo is excluded before MQ is applied.

The analysis is deliberately catalog-first. All 175 SPARC galaxies are retained in the principal analysis before subdivision by observational quality or galaxy property. SPARC itself is not a statistically complete or volume-limited census of nearby disk galaxies. Accordingly, “complete SPARC sample” throughout this work denotes evaluation of every galaxy in the 175-object SPARC catalog rather than a complete astrophysical population [4]. Residual structure is subsequently tested against SPARC’s independently assigned rotation curve quality classification Q, inclination, distance, published flat rotation velocity Vflat, luminosity, surface brightness, morphology, gas content, and radial coverage.

Across all radial measurements, the point-weighted signed residual is +0.02644 km s-1 and the ratio of mean MQ realized to mean observed velocity is 1.0001888, corresponding to an aggregate velocity-scale difference of +0.0189%. These quantities characterize signed catalog-wide normalization rather than pointwise accuracy: the corresponding pooled MAE is 16.20 km s-1 and the pooled RMS residual is 22.82 km s-1.

The complete catalog exhibits a statistically significant velocity-dependent normalization residual. Among lower velocity galaxies, the observed rotation curves tend to lie below the MQ predicted curve, whereas among higher velocity galaxies they tend to lie above the MQ predicted curve. Section 4.6 quantifies this relation, whose interpretation changes materially when the sample is stratified by SPARC’s independently assigned observational quality flag Q [4]. Among the 99 galaxies classified by SPARC as high quality (Q=1), the geometric MQ-to-observed velocity ratio is 1.0042 and the residual dependence on characteristic rotational velocity is not statistically significant. Sections 4.7–4.10 examine how the complete-catalog trend changes under observational quality, inclination, distance, Vflat, stellar source, galaxy-property, and curve-shape controls. These analyses distinguish the absence of a significant marginal velocity relation in the highest quality class from the more limited question of whether a conditional velocity dependence remains when correlated galaxy properties are considered jointly.

SPARC’s quality classifications reflect properties of the underlying kinematic data, including asymmetries, non-circular motions, and overall data quality [4]. Independent work on dwarf-galaxy tilted-ring recovery likewise shows that non-circular motions can bias inferred circular velocities even when conventional residual diagnostics appear modest [45]. In the present catalog, observational quality is also associated with galaxy regime. The quality-stratified analysis must therefore be interpreted as an independently defined observational control and not as evidence that the larger residuals among lower quality galaxies are observational errors.

Throughout Sec. 4, the same MQ realization law is applied with zero galaxy-specific fitted MQ parameters: no galaxy-specific acceleration scale, empirical interpolation function, or halo profile is adjusted to the rotation curves. The purpose is to determine whether a single realization law constructed from the fixed MQ acceleration scale, the positive quadratic System Frame norm, and subdivision-consistent residual composition remains consistent across independent baryonic reconstructions and observational conditions, rather than to optimize individual rotation curves. Sections 4.1–4.3 establish bounded gravitational realization, the observational baryonic inputs, and the fixed acceleration scale and local realization law; Secs. 4.4–4.10 provide the independent observational tests and their controls; and Sec. 4.11 gives their combined interpretation. This does not make the complete astrophysical analysis parameter-free: the Newtonian baryonic source retains published quantities and modeling assumptions such as stellar mass-to-light normalization, distance, inclination, gas distribution, and the component decompositions supplied by the cited datasets.

4.1 Observational implementation and baryonic mass models

Evaluation of the MQ galactic realization law requires a Newtonian baryonic velocity profile derived from an externally published baryonic mass model. MQ does not reconstruct the stellar or gaseous mass distribution. It accepts the adopted baryonic reconstruction as the source input and predicts the corresponding realized gravitational response. The Milky Way and SPARC analyses use independent observational velocity data and different baryonic reporting conventions, so the Newtonian baryonic velocity vNew is constructed separately for the two datasets. For the Milky Way, the observed circular velocities are taken from Eilers et al. [2], whereas the baryonic source is taken independently from the Galactic mass model of McMillan [3]. For SPARC, both the observed rotation curves and baryonic component tables are taken from the SPARC database [4].

For the principal Milky Way calculation, the baryonic source is reconstructed from the best-fit axisymmetric potential of McMillan [3], which was developed without an RAR or MOND prescription. The retained baryonic components are the bulge, thin and thick stellar disks, H I disk, and molecular gas disk. Their published density laws and best-fit normalizations are evaluated directly as Newtonian gravitational sources in the Galactic mid-plane. The NFW dark halo force is excluded before construction of the MQ input, so no dark matter contribution enters the baryonic acceleration used by Eq. (137). At each of the 38 Eilers radii, the radial forces of the five baryonic components are numerically integrated and summed. The Newtonian baryonic circular speed is then obtained from the resulting radial acceleration. As an implementation check independent of the Eilers target velocities, restoring McMillan’s published NFW halo at R0=8.21 kpc gives a total circular speed of 233.109 km s-1, reproducing McMillan’s quoted 233.1 km s-1 best-fit value [3]. Excluding the halo force from the MQ source does not remove the covariance of McMillan’s fitted baryonic parameters with the halo and other assumptions in the original Galactic mass-model fit. That dependence remains part of the baryonic source uncertainty considered in Sec. 4.4.

For the RAR-conditioned source model sensitivity comparison retained from the earlier analysis, table 2 of McGaugh [46] provides the Newtonian circular-speed contributions vbulge, vgas, and vdisk associated with the adopted Galactic bulge, gas, and stellar-disk components. These quantities already represent the component contributions of the adopted Galactic mass model and therefore require no additional stellar mass-to-light normalization. Because the McGaugh radii differ from the 38 observational radii of Eilers et al. [2], each of the three component-velocity series is independently linearly interpolated to the Eilers radii before construction of that alternative Newtonian baryonic velocity. This route is no longer the principal Milky Way test. Their radial acceleration contributions then combine as

vNew(R)=vbulge2(R)+vgas(R)|vgas(R)|+vdisk2(R).
(129)

The signed gas term in Eq. (129) is required because the McGaugh Galactic model contains negative gas contributions at some radii, corresponding to an outward radial contribution from the adopted gas distribution [46]. Squaring vgas without retaining its sign would incorrectly convert that outward contribution into an inward one. McGaugh also tabulates a total model circular velocity Vc, but that quantity is not used as observational data in the present analysis because the McGaugh model incorporates the empirical radial acceleration relation [46]. Instead, after vNew(R) and the MQ realized velocity have been calculated, the prediction is compared with the independent circular velocity measurements of Eilers et al. [2]. Thus, the McGaugh sensitivity reconstruction and the observational velocity used to test it originate from distinct data products. This sensitivity calculation is retained only to quantify source model dependence. The principal McMillan calculation removes this RAR/MOND source-side dependence by using a conventional baryonic density model and excluding its dark halo force before applying MQ.

The 38 Milky Way observational points span R=5.27-24.82 kpc. Their observed circular velocities and asymmetric statistical uncertainties are taken from Eilers et al. [2]. For the total error analysis developed in Sec. 4.4, the radius-dependent systematic fractions compiled by Jiao et al. [47] are converted to velocity uncertainties and combined pointwise with the corresponding asymmetric statistical uncertainties. These uncertainties enter only the subsequent comparison with the prediction. Neither the Eilers velocities nor their uncertainties are used to construct vNew, gbar, or the MQ realized velocity.

SPARC instead tabulates the stellar disk and bulge velocity contributions for unit stellar mass-to-light ratio [4]. Following the 3.6 μm stellar population normalization adopted in the SPARC radial acceleration analysis [8], Υdisk=0.5 and Υbulge=0.7, the SPARC Newtonian baryonic velocity is constructed from

vNew(R)=0.7vbulge2(R)+vgas(R)|vgas(R)|+0.5vdisk2(R).
(130)

The signed gas term likewise preserves the direction of the tabulated gaseous radial contribution when SPARC reports a negative gas contribution. This convention is part of the published SPARC mass-model construction and is required because a centrally depressed gas distribution can produce a formally negative radial contribution in the inner galaxy [4]. Equation (129) defines the legacy McGaugh sensitivity construction, whereas Eq. (130) defines the SPARC construction; both construct the same physical quantity, vNew, but from differently normalized baryonic inputs. The McGaugh Milky Way component velocities already represent the adopted Galactic component model, whereas the SPARC stellar component velocities are supplied for unit mass-to-light ratio and require the fixed normalizations above.

Because the MQ prediction depends on the adopted baryonic acceleration, the fixed SPARC stellar normalization is an astrophysical source assumption whose uncertainty must be propagated separately from MQ calibration. Near-infrared stellar population studies support approximately constant 3.6 μm stellar mass-to-light ratios but not a unique exact value: representative published disk normalizations include about

0.47M⊙/L⊙
(S59)

[48] and

0.60M⊙/L⊙
(S60)

[49], while a broader literature synthesis places commonly used constant values between about 0.42 and 0.60 [50]. Section 4.6 therefore performs two source-side tests without fitting any galaxy to the MQ residual: a fixed-global normalization sweep and a second calculation in which the stellar normalization of each overlapping galaxy is set by the independently derived WISE stellar mass published by Duey et al. [51]. The latter changes the stellar source entering Eq. (130) before Eq. (137) is evaluated and is therefore distinct from the matched-residual diagnostic of Sec. 4.8.

The corresponding Newtonian baryonic acceleration at galactocentric radius R is

gbar(R)=vNew2(R)R.
(131)

Accordingly, both vNew and gbar are properties of the adopted observational baryonic reconstruction. They are not MQ fitting parameters. Section 4.3 uses gbar as the local baryonic input to the fixed galactic realization law.

The baryonic reconstructions are necessarily model dependent. Their uncertainties include stellar mass normalization, bulge-disk photometric decomposition, gas distribution, distance, spatial resolution, and other assumptions entering the conversion from observed matter distributions to a Newtonian baryonic profile. For the principal Milky Way test, the McMillan baryonic density model [3] is distinct from the Eilers Jeans reconstruction used to infer the observed circular velocities [2] and does not employ an RAR or MOND relation. Consequently, uncertainties associated with the Eilers tracer-density profile, axisymmetry assumption, Galactic parameters, and Jeans analysis affect the observational comparison but are not uncertainties in the McMillan baryonic component calculation itself. For SPARC, distance, inclination, stellar mass-to-light normalization, gas distribution, non-circular motion, and related observational assumptions propagate through the catalog baryonic reconstruction and rotation curve measurements [4]. These are astrophysical modeling inputs and observational uncertainties, not fitted MQ parameters; no claim that the complete rotation curve analysis is parameter-free is made. The remaining source model dependence is ordinary astrophysical modeling dependence, including the adopted stellar and gas distributions, vertical profiles, component normalizations, and conventional Galactic constraints entering McMillan [3]. This does not reintroduce RAR/MOND phenomenology into the baryonic source.

For SPARC specifically, let

evo
(S61)

denote the published velocity uncertainty associated with the observed orbital velocity vo. This uncertainty combines the formal rotation curve fitting uncertainty with the difference between the approaching and receding sides of the disk and therefore reflects both local deviations from circular motion and global kinematic asymmetry [4]. It does not include the systematic uncertainty arising from the adopted inclination [4]. Consequently, all error-normalized statistics based directly on

evo
(S62)

are identified explicitly as random-error-only diagnostics. For a radial measurement, let zV,rand denote the standardized random-error velocity residual,

χrand2
(S63)

the corresponding per-galaxy chi-square statistic, ν the number of degrees of freedom, and

χν,rand2
(S64)

the reduced chi-square statistic. The SPARC catalog does not provide a common radial covariance matrix for these 175 heterogeneous rotation curves, and the individual curves originate from multiple H I and H-alpha source studies with different reduction pipelines [4]. A full covariance likelihood for all 3,366 radial measurements therefore cannot be reconstructed from the standardized SPARC products without imposing additional correlation assumptions. The principal uncertainty statements below consequently separate pointwise random-error diagnostics from galaxy-level population inference.

zV,rand=vreal-voeVo,χrand2=∑zV,rand2,χv,rand2=χrand2ν.
(132)

Here ν=N, where N denotes the number of radial measurements retained for the galaxy, because no galaxy-specific MQ parameter is estimated from those measurements. The statistical quantity

χrand2
(S65)

is distinct from the realization-depth quantity χ(R) introduced in Sec. 4.3. The associated right-tail χ2 probability and the fractions of radial measurements lying within one, two, and three quoted observational standard deviations are likewise reported only as diagnostics relative to the published pointwise velocity uncertainties. They are not interpreted as complete likelihood measures for the MQ galactic realization law because the underlying error budget does not include all observational or baryonic model systematics.

Distance, inclination, SPARC quality classification, published flat rotation velocity, and the additional galaxy metadata used later in the SPARC analysis are retained as independently reported observational quantities rather than absorbed into Eq. (132) or used to adjust the MQ prediction. Their effects on SPARC residual structure are tested separately in Secs. 4.7–4.9. Likewise, the independent Eilers observational uncertainties and systematic controls enter the Milky Way validation only after the Galactic baryonic reconstruction and MQ prediction have been constructed, as detailed in Sec. 4.4. Thus, in both datasets, galaxy-specific information entering the predicted velocity originates in the adopted baryonic reconstruction, while observational velocity measurements are retained as external comparison quantities. The MQ galactic realization law remains fixed across the Milky Way and the complete SPARC catalog. Its acceleration scale and local realization form are introduced next.

4.2 MQ acceleration scale and local realization law

The required acceleration scale follows from the frame-dependent cosmological expansion derived in Sec. 3.12. Using the Internal Frame local rate realization of universal expansion Hint, the defined speed of light c, and the invariant coefficient θsi, MQ gives

aMQ=cHint2θsi-1.
(133)

At a specified cosmological epoch, aMQ is fixed independently of galactic rotation curve data and is applied unchanged to every galaxy evaluated at that epoch. Using the present epoch quantities entering the analysis gives

aMQ=1.20040492178×10-10ms-2
(S66)

, which is the value used throughout the complete SPARC calculation. The frame-dependent expansion quantity is derived in Sec. 3.12. The bounded realization construction is developed in Appx. E, and the subdivision-consistent derivation of the realization fraction in Appx. F. Neither the complete SPARC population nor subsequent observational-quality subdivisions determine or modify aMQ. The same present-epoch value is used for all 175 galaxies and all 3,366 radial measurements before residual analysis. Accordingly, “independently specified” in this paper means derived within MQ from independently measured non-galactic inputs, with no rotation curve calibration; it does not mean numerically detached from empirical physics. The coefficient θsi used in the construction is established in Sec. 3.8 from the measured fine structure constant, and the frame-resolved expansion rate inherits the cosmic age normalization specified in Sec. 3.12.

Equation (136) shares its mathematical functional form with the empirical radial acceleration relation reported by McGaugh, Lelli, and Schombert [8], and that prior phenomenology is acknowledged solely as prior art for the functional shape. The fitted characteristic acceleration appearing in an empirical RAR regression is not treated as a measurement of aMQ. The empirical scale is a model-conditioned regression parameter whose value depends on the adopted baryonic reconstruction and statistical nuisance structure, whereas aMQ is derived upstream from the MQ frame-resolved cosmological construction and fixed before galaxy data are evaluated. A numerical similarity or difference between such parameters has no direct evidentiary meaning unless an independent derivation establishes that they represent the same physical observable. Accordingly, this work assigns neither concordance nor tension significance to comparisons between aMQ and a fitted RAR scale. The physically relevant test is direct: construct the Newtonian baryonic source independently, apply Eqs. (136) and (135) with aMQ frozen and compare the resulting velocity with an independently determined circular velocity.

For the baryonic acceleration defined in Sec. 4.2, the local realization depth follows from the same positive quadratic System Frame norm used by the MQ geometric construction. Reference [1], Proposition 2.4 and Appx. AI, establish that a regular direction-free local realization is represented by an isotropic inner-product norm rather than a freely selectable power norm, while Appx. N supplies the corresponding realized tangent metric. For circular motion, Eq. (131) gives the positive Newtonian speed magnitude from the local baryonic acceleration and radius. Applying the independently fixed MQ acceleration scale to the same circular-kinematic map at that radius gives the corresponding reference speed magnitude. The Frames mapping normalizes the variable realized magnitude against the fixed reference magnitude itself, not against its squared norm. Their ratio is dimensionless, and the common radial factor cancels exactly. The local realization depth is therefore

χ(R)=gbar(R)aMQ.
(134)

Equation (134) is consequently a norm ratio rather than a phenomenological power-law ansatz. The one-half power follows from the positive quadratic System Frame metric together with the circular relation between acceleration and squared speed. Because both acceleration scales are evaluated through the same local radius, no radius-dependent normalization or galaxy-specific parameter is introduced. Related scale-invariant galactic dynamics likewise obtains the deep-acceleration square-root structure by combining an acceleration scale with circular kinematics [19,21,52], although in MQ both the acceleration normalization and the quadratic realization norm have independent upstream provenance. Appendix E gives the bounded realization construction, while Appx. F derives the realization fraction from subdivision-consistent realization of the remaining gravitational availability:

P(R)=1-exp⁡[-χ(R)]=1-exp⁡[-gbar(R)aMQ].
(135)

The realized acceleration greal(R) is therefore

greal(R)=gbar(R)1-exp⁡[-gbar(R)/aMQ].
(136)

As noted above, the functional identity of Eq. (136) with the empirical RAR is prior art rather than a novelty claim. The tested distinction is its MQ provenance: Eq. (133) supplies the characteristic acceleration independently, Eq. (134) follows from the quadratic System Frame norm and Newtonian circular kinematics rather than from the observed radial acceleration or baryonic Tully-Fisher relations, and Appx. F derives the exponential realization fraction from subdivision consistency of the remaining gravitational availability. The galaxy data enter only after the complete galactic realization law has been specified. Because the empirical RAR scale is not a physically identical MQ observable, no fitted-scale concordance statistic is used; normalization and residual structure are assessed only through direct velocity comparisons.

Using

greal(R)=vreal2(R)/R
(S67)

together with

gbar(R)=vNew2(R)/R
(S68)

, the predicted orbital velocity becomes

vreal(R)=vNew(R)1-exp⁡[-gbar(R)/aMQ].
(137)

The limiting behavior follows from Eq. (134). For gbar≫aMQ, the exponential contribution vanishes, so greal→gbar and vreal→vNew, recovering the Newtonian regime. In the opposite limit,

gbar/aMQ→0
(S69)

, the first-order expansion of the exponential denominator derived in Appx. F gives the exact asymptotic statement

limgbar/aMQ→0greal(R)gbar(R)​aMQ=1.
(138)

For an exterior baryonic field with

gbar(R)=GMbar(<R)/R2
(S70)

, where Mbar(<R) is the previously defined total baryonic mass enclosed within galactocentric radius R, this weak-acceleration behavior gives

limgbar/aMQ→0vreal4GMbaraMQ=1.
(139)

Thus, the realized acceleration approaches the Newtonian baryonic field at high acceleration and has asymptotic

gbaraMQ
(S71)

scaling at low acceleration. The latter yields an asymptotically flat outer rotation curve together with the baryonic Tully-Fisher fourth-power scaling.

Once the baryonic reconstruction supplies the Newtonian baryonic velocity and acceleration, Eq. (137) fixes the predicted orbital velocity. The frame-resolved cosmological construction supplies the acceleration scale, the positive quadratic System Frame norm fixes the dimensionless realization depth, and subdivision consistency fixes the realization fraction. None of these steps uses the observed orbital velocity to construct the prediction. No galaxy-specific acceleration scale, halo profile, or other MQ parameter is fitted to the rotation curve. Consequently, the quality-, inclination-, velocity-, and galaxy-property dependencies examined subsequently are diagnostics of residuals obtained after application of the fixed galactic realization law rather than parameters used to construct it.

4.3 RAR/MOND-independent Milky Way rotation curve test

The Milky Way provides an RAR/MOND-independent source-and-kinematic test of the fixed galactic realization law specified in Sec. 4.2, using an observational target distinct from the external galaxy SPARC analysis. The observed circular velocities are taken from Eilers et al. [2], who inferred the Galactic circular velocity curve from the six-dimensional phase-space coordinates of more than 23,000 luminous red giant stars using the Jeans equation under an axisymmetric gravitational potential assumption. Their analysis covers approximately 5≤R≤25 kpc and excludes the inner 5 kpc, where the Galactic bar requires a non-axisymmetric treatment [2]. The present calculation uses the 38 Eilers measurements spanning R=5.27–24.82 kpc together with their asymmetric statistical uncertainties. Jiao et al. ([47], Appx. A) reproduce these measurements and provide the radius-dependent systematic fractions adopted below for the pointwise total error analysis.

The principal baryonic source is constructed from the best-fit McMillan [3] axisymmetric density model. The calculation retains the bulge, thin and thick stellar disks, H I disk, and molecular gas disk, evaluates their Newtonian mid-plane radial forces, and excludes the model NFW dark halo before the MQ galactic realization law is applied. McMillan does not employ an empirical RAR or MOND prescription, so neither an RAR, a MOND interpolation function, nor a dark-halo force enters construction of the MQ source. The Eilers velocities provide a separately derived kinematic target and are not used to construct either the McMillan baryonic source or the MQ prediction. The earlier McGaugh [46] component reconstruction is retained only as a source model sensitivity comparison because it was obtained within an RAR-conditioned Galactic model.

To reproduce the principal calculation, the published McMillan density laws and best-fit component normalizations are evaluated directly at each of the 38 Eilers radii using axisymmetric Newtonian force integration. The radial accelerations from the bulge, thin stellar disk, thick stellar disk, H I disk, and molecular gas disk are summed and converted to the Newtonian baryonic circular speed. McMillan’s NFW halo is omitted, no parameter is re-estimated from the Eilers velocities, and the resulting baryonic acceleration is supplied unchanged to Eq. (137). The accompanying calculation workbook reports every component contribution and derived velocity used here. Equation (129) is retained only for the separate McGaugh source model sensitivity calculation.

The resulting vNew(R) determines the Newtonian baryonic acceleration gbar(R) through Eq. (131), after which Eq. (137) supplies the realized orbital velocity vreal(R). The observed orbital velocity vo(R) enters only after both model velocities have been calculated and is used to evaluate the residuals. No Milky Way-specific MQ acceleration scale, dark-halo profile, or other MQ parameter is fitted. Thus, this test asks whether the fixed MQ galactic realization law maps a conventional baryonic Galactic reconstruction onto an independently determined kinematic curve without introducing the RAR/MOND source dependence present in the earlier McGaugh sensitivity calculation.

The uncertainty treatment is likewise kept separate from construction of the prediction. The Eilers measurements have asymmetric statistical uncertainties. For the diagonal total error calculation, the radius-dependent systematic fraction reported by Jiao et al. [47] is converted to an absolute velocity uncertainty and combined in quadrature with the appropriate upper or lower statistical uncertainty. The resulting dataset retains separate statistical, systematic, and asymmetric total error quantities for every radial point. This pointwise construction is not a complete covariance model. Eilers et al. identify systematic sensitivities associated with the Jeans treatment and Galactic parameters [2], while Jiao et al. distinguish radius-dependent effects from additional systematics that can coherently affect the velocity scale [47]. The resulting diagonal-error χ2 values are therefore comparative diagnostics rather than covariance-complete likelihood statements.

No full 38×38 covariance matrix for the Eilers circular velocity points is supplied by the observational products used here. The baryonic source uncertainty is different: McMillan [3] publishes posterior standard deviations together with the full eight-parameter correlation matrix for the fitted Galactic mass model. Those published posterior quantities permit a limited baryonic source covariance to be propagated without inventing an Eilers covariance or refitting MQ.

For that calculation, the baryonic normalization parameters entering the present force decomposition linearly are the thin-disk surface-density normalization, thick-disk surface-density normalization, and bulge-density normalization. Because the present calculation holds the disk scale lengths and remaining McMillan parameters at their published best-fit values, a three-parameter conditional covariance is constructed from McMillan’s published posterior standard deviations and eight-parameter correlation matrix. This propagates the correlated uncertainty of the three varied source normalizations conditional on the fixed parameters rather than attributing their larger marginal covariance to quantities not varied in this calculation.

Each source realization rescales only those three baryonic density normalizations, recomputes the Newtonian baryonic velocity from the tabulated component forces, and applies Eq. (137) with the MQ acceleration scale frozen. A 200,000-draw multivariate Monte Carlo propagation gives a correlated one-standard-deviation source uncertainty ranging from 2.09 to 3.17 km s-1 across the 38 Eilers radii. Adding this source covariance to the adopted asymmetric Eilers/Jiao observational diagonal gives generalized χ2=39.57 for 38 degrees of freedom, reduced χ2=1.041, and right-tail p=0.400. The inner 30-point subset gives χ2=30.85 (p=0.423), and the outer eight points give χ2=3.356 (p=0.910). No MQ parameter or covariance amplitude is estimated from the residuals.

This calculation remains narrower than a complete joint covariance analysis. It does not vary the McMillan disk scale lengths, gas-disk parameters, or other source model quantities, and it does not supply the unpublished Eilers point-to-point covariance. It does, however, propagate a physically motivated covariance component constructed from independently published McMillan posterior information that is absent from the diagonal test. The change from

p=1.26×10-120
(S72)

to p=0.400 therefore demonstrates that the formal diagonal rejection is highly sensitive to omission of baryonic source covariance. The appropriate inference is non-rejection under this specified source-covariance model, not that the complete physical covariance has been measured.

Figure from the preprint

Figure 3. RAR/MOND-independent Milky Way rotation curve test. Points show the 38 Eilers et al. [2] circular-velocity measurements from R=5.27 to 24.82 kpc with asymmetric total error bars. The Newtonian baryonic curve is calculated from the McMillan [3] bulge, thin and thick stellar disks, H I disk, and molecular-gas disk by direct axisymmetric force integration; the McMillan NFW halo is excluded. The MQ curve is obtained from Eq. (139) with the fixed MQ acceleration scale. Neither an RAR/MOND relation nor the Eilers target velocities enters construction of the prediction.

Figure 3 compares the 38 Eilers measurements with the Newtonian baryonic and MQ realized velocity curves. Defining the signed residual as MQ realized minus observed velocity, the complete 5.27–24.82 kpc sample gives an MQ RMS residual of 8.8965 km s-1, an MAE of 7.4312 km s-1, a maximum absolute residual of 25.6269 km s-1, and a mean signed residual of -3.07445 km s-1. The mean absolute relative error is 3.54277%, the median absolute error is 6.32895 km s-1, and the Pearson correlation between the MQ realized and observed velocities is rP=0.88995. The corresponding Newtonian baryonic RMS residual is 69.0279 km s-1, so application of the fixed galactic realization law reduces the full-sample RMS discrepancy by 87.112%. This reduction is obtained from a McMillan baryonic source constructed without RAR/MOND phenomenology and with its NFW dark-halo force excluded.

Because the Eilers uncertainties increase markedly beyond approximately 20 kpc, the full sample is also divided into an inner 30-point region and an outer 8-point region. Table 7 summarizes the principal residual, correlation, and diagonal-error statistics for the complete, inner, and outer samples. The inner region provides the highest-precision comparison and an 89.683% RMS reduction relative to the Newtonian baryonic prediction. The outer region has an MQ RMS residual of 14.0857 km s-1 and substantially larger observational uncertainties, giving a non-rejecting diagonal total error probability despite the larger pointwise scatter.

Table 7. Milky Way residual and diagonal-error statistics for complete, inner, and outer radial samples.

Statistic

Complete

Inner

Outer

Radial range (kpc)

5.27–24.82

5.27–19.71

20.27–24.82

N
(S73)

38

30

8

MQ RMS (km s-1)
(S74)

8.8965

6.8807

14.0857

Newtonian RMS (km s-1)

69.0279

66.6906

77.1648

RMS reduction

87.112%

89.683%

81.746%

MQ MAE (kms-1)
(S75)

7.4312

6.3961

11.3127

MQ mean signed residual (km s-1)

−3.0744

−6.1930

+8.6200

MQ Pearson r

0.8900

0.9502

0.3178

MQ χtot2
(S76)

689.831

685.998

3.8326

Reduced χ2

18.1534

22.8666

0.4791

Right-tail p

1.26×10-120
(S77)
4.07×10-125
(S78)

0.8719

Newtonian

χtot2
(S79)

60,621.054

60,268.028

353.026

Table 8 retains the earlier observational-covariance stress test. The parameter ρ applies only to the systematic component inferred from the adopted total and statistical uncertainties. This is a sensitivity family, not a published Eilers covariance model.

Table 8. Milky Way covariance-correlation stress test. The parameter rho applies only to the systematic component inferred from the adopted total and statistical uncertainties. This is a sensitivity family, not a published Eilers covariance model.

rho

Full chi-square

Full reduced chi-square

Full p

Inner chi-square

Outer chi-square

Outer p

0.00

689.831

18.153

1.26 x 10^-120

685.998

3.8326

0.8719

0.25

578.411

15.221

8.37 x 10^-98

519.571

3.2914

0.9148

0.50

632.313

16.640

8.17 x 10^-109

555.631

3.3980

0.9070

0.75

738.310

19.429

1.27 x 10^-130

638.744

3.8355

0.8717

1.00

1196.633

31.490

2.22 x 10^-226

1006.150

4.8927

0.7690

Minimum

571.087 at rho=0.154

15.029

2.59 x 10^-96

-

-

-

Table 9 gives the generalized chi-square after propagation of the McMillan baryonic source covariance described above.

Table 9. Milky Way generalized chi-square after propagation of McMillan baryonic source covariance. The source covariance varies only the thin-disk, thick-disk, and bulge normalizations using the three-parameter conditional covariance constructed from McMillan’s published posterior information [3]; no MQ parameter is fitted.

Model

Full chi-square

Reduced chi-square

Full p

Inner chi-square / p

Outer chi-square / p

Source sigma (km s-1)

Diagonal only

689.831

18.153

1.26e-120

685.998 / 4.07e-125

3.833 / 0.872

0

+ McMillan source covariance

39.57

1.041

0.400

30.85 / 0.423

3.356 / 0.910

2.09-3.17

For the complete sample, the statistical-only error treatment gives χ2=2602.440 and reduced χ2=68.485 for 38 measurements, with right-tail probability below

10-300
(S80)

. Under the diagonal total error treatment summarized in Table 7, the RAR/MOND-independent McMillan source gives χ2=689.831, reduced χ2=18.153, and

p=1.26×10-120
(S81)

. The Newtonian baryonic curve gives χ2=60,621.054, so MQ reduces this diagonal statistic by 98.86%. Table 9 shows that this diagonal probability is not robust to inclusion of the propagated McMillan normalization covariance: under that specified source-covariance model, the same frozen prediction gives generalized χ2=39.57 and p=0.400.

The equicorrelation stress test in Table 8 isolates a different uncertainty: the still-unknown observational point-to-point covariance. Its most favorable complete-sample value occurs at ρ=0.154, where χ2=571.09, reduced χ2=15.03, and the 38-degree-of-freedom right-tail probability is

2.59×10-96
(S82)

. Across 0≤ρ≤1, the inner subset remains formally rejected when the McMillan source covariance is omitted. Thus, a single common-correlation coefficient applied only to the retained Eilers/Jiao systematic component does not reproduce the effect of propagating uncertainty in the baryonic source.

A second observational-covariance diagnostic uses the radial Gaussian-process kernel form studied independently for Milky Way rotation curves by Oman and Riley [53]. Their covariance model adds a correlated term whose amplitude decays with radial separation on a characteristic length scale; when fitted to the Ou et al. Milky Way rotation curve, they infer correlations below 10 km s-1 over approximately 1.5–2.5 kpc and show that such covariance can materially alter Galactic mass inference. That fitted covariance is not an Eilers covariance measurement and is therefore not transferred to the present data. Its published kernel form and radial-scale range are retained only to quantify the observational-covariance scale that would be consequential if an Eilers covariance became available.

For correlation lengths of 1.5, 2.0, and 2.5 kpc, the diagonal-only threshold amplitudes are 2.07, 1.93, and 1.87 km s-1, respectively. An amplitude of 2.0 km s-1 gives diagonal-plus-kernel χ2 values of 56.15, 50.95, and 49.03, corresponding to right-tail probabilities of 0.029, 0.078, and 0.108. These are sensitivity calculations rather than goodness-of-fit claims because the kernel parameters were not measured from the Eilers data. Their role is secondary to Table 9: they show independently that observational covariance of a magnitude already inferred in other Galactic rotation curve work could also materially affect the precision statistic.

The inner sample contains all 30 measurements from 5.27 to 19.71 kpc. It gives an MQ RMS residual of 6.8807 km s-1, MAE of 6.3961 km s-1, and mean signed residual of -6.1930 km s-1, compared with a Newtonian baryonic RMS of 66.6906 km s-1. The RMS discrepancy is reduced by 89.683%. The diagonal total error calculation gives χ2=685.998 and

p=4.07×10-125
(S83)

, whereas propagation of the McMillan conditional baryonic-normalization covariance gives generalized χ2=30.85 for 30 degrees of freedom and p=0.423.

The remaining eight measurements lie between 20.27 and 24.82 kpc. Their MQ RMS residual is 14.0857 km s-1, MAE is 11.3127 km s-1, and mean signed residual is +8.6200 km s-1, compared with a Newtonian baryonic RMS of 77.1648 km s-1. Under the diagonal total error treatment, χ2=3.8326 and p=0.8719; after addition of the McMillan conditional baryonic-normalization covariance, generalized χ2=3.356 and p=0.910. The outer subset therefore remains non-rejected, but its much larger observational uncertainties give it substantially less constraining power than the inner measurements.

A radial-subset scale diagnostic further shows that the remaining disagreement is not described by a single fitted acceleration normalization. An unweighted RMS fit performed only as a diagnostic on the McMillan-Eilers comparison gives a complete-sample optimum acceleration scale of

1.28389×10-10
(S84)

m s⁻², 6.95% above the frozen MQ value. The inner 30 points alone prefer

1.40922×10-10
(S85)

m s⁻², whereas the outer eight prefer

0.96494×10-10
(S86)

m s⁻². The radial disagreement among these fitted optima indicates that a single post hoc scale adjustment cannot absorb the remaining curve-shape residual structure. The proximity of the full-sample optimum to the independently fixed MQ scale is therefore a normalization diagnostic, not a substitute for a covariance-complete precision test.

The detailed point-to-point excursions in Fig. 3 should not automatically be interpreted as literal local gravitational structure. Eilers et al. infer circular velocities through a Jeans analysis and identify sensitivity to assumptions involving axisymmetry, tracer structure, Galactic parameters, and neglected or approximated dynamical terms [2]. Jiao et al. likewise emphasize systematic error treatment in interpreting this high-precision rotation curve [47]. McGaugh’s subsequent Jeans-adjusted analysis provides an additional demonstration that modifications to the stellar-density gradient term can alter the inferred stellar rotation curve [46], but that reanalysis is not an independent observational dataset. As a source model sensitivity check, the earlier RAR-conditioned McGaugh reconstruction yielded an MQ RMS of 11.08 km s-1, compared with a Newtonian RMS of 58.57 km s-1. The principal McMillan result is stronger, with an MQ RMS of 8.90 km s-1 versus 69.03 km s-1 and does not inherit an RAR/MOND source relation.

Against the independent Eilers observations, the RAR/MOND-independent McMillan construction therefore provides a source-provenance-controlled Milky Way test of the fixed MQ galactic realization law. It removes most of the Newtonian baryonic velocity discrepancy and remains statistically non-rejected when the specified McMillan normalization covariance is propagated, without fitting an MQ parameter to the rotation curve. The statistical significance assigned to the remaining residuals is covariance-sensitive: the diagonal Eilers/Jiao error model rejects the complete and inner samples, whereas the specified McMillan conditional source covariance gives χ2=39.57 for 38 degrees of freedom (p=0.400), with the inner and outer subsets likewise non-rejected. A published Eilers point-to-point covariance and a joint propagation of the complete McMillan source posterior remain unavailable in the present calculation, so the result establishes consistency under the propagated source-covariance model rather than covariance-complete precision validation.

4.4 Complete SPARC rotation curve test
The principal catalog-wide test applies the fixed MQ galactic realization law to all Ngal=175 galaxies in the SPARC database, where Ngal denotes the number of galaxies, comprising Nradial=3,366 radial velocity measurements, where Nradial denotes the total number of radial measurements [4]. Every catalog galaxy is retained in this first-stage analysis. No selection by SPARC quality classification or other galaxy property is applied, and no galaxy-specific MQ parameter is fitted. The galactic realization law from Sec. 4.2 is evaluated against each published baryonic reconstruction before observational-quality or galaxy-property stratification.

SPARC is not a statistically complete or volume-limited census of nearby disk galaxies. Accordingly, “complete SPARC catalog” here denotes evaluation of all 175 galaxies in the published database rather than a complete astrophysical population [4]. Because the functional form of Eq. (136) is pre-existing empirical RAR phenomenology, this calculation is not treated as a discovery test of that shape. It instead tests whether the independently specified MQ acceleration scale and galactic realization law reproduce the catalog relation without galaxy-specific MQ calibration; the published baryonic source assumptions remain part of the input construction.

For all subsequent SPARC statistics, let Δvres denote the signed velocity residual,

Δvres=vreal-vo,
(140)

so that positive residuals indicate an MQ realized orbital velocity above the observed orbital velocity and negative residuals indicate an MQ realized orbital velocity below it.

Pooling all 3,366 radial measurements gives

⟨Δvres⟩pooled=+0.02644kms-1,
(141)

with point-weighted mean velocities

⟨vo⟩=140.0368kms-1,⟨vreal⟩=140.0632kms-1,
(142)

and therefore

⟨vreal⟩⟨vo⟩=1.0001888.
(143)

The complete catalog is therefore centered within approximately +0.0189% in point-weighted aggregate velocity normalization. This near-zero signed offset does not imply pointwise agreement at the same fractional level. Across the same 3,366 measurements, the pooled mean absolute error is 16.1982 km s-1, the pooled RMS residual is 22.8203 km s-1, and the point-weighted MAPE is approximately 17.67%. Equations (141–143) describe the catalog-wide signed velocity scale, whereas MAE, RMS, and MAPE describe the magnitude of individual radial departures. The near-unity normalization in Eq. (143) must therefore be interpreted together with these absolute-error measures rather than as pointwise agreement.

The radial measurements are nested within galaxies, and individual galaxies contribute different numbers of points. The 3,366 measurements therefore should not be treated as statistically independent realizations of 3,366 separate systems. Correlations among rotation curve measurements at different radii can also arise from shared observational and reconstruction quantities. The pooled quantities above are consequently descriptive point-weighted catalog summaries rather than significance estimates based on 3,366 independent observations.

To account explicitly for this hierarchy, the catalog summaries are supplemented by a nonparametric galaxy-cluster bootstrap [54]. Ten thousand bootstrap catalogs are generated with a fixed seed by resampling the 175 galaxies with replacement while retaining each selected galaxy’s complete radial-point aggregate as a single cluster. The resulting 95% percentile interval for the point-weighted pooled signed residual is -3.80 to +3.74 km s-1 around the observed +0.0264 km s-1 value. The corresponding intervals are 19.64-26.03 km s-1 for pooled RMS, 13.94-18.61 km s-1 for pooled MAE, and 0.9742-1.0286 for the ratio of point-weighted mean MQ realized to observed velocity. These intervals show that the near-zero pooled bias and near-unity aggregate normalization are compatible with zero bias and unity normalization at the galaxy-cluster level.

Treating each galaxy as one statistical unit gives a complementary complete-catalog result. Averaging each galaxy’s mean signed residual with equal galaxy weight gives a mean bias of +3.3195 km s-1 and a median bias of +3.7614 km s-1. Of the 175 galaxies, 112 have positive mean residuals and 63 have negative mean residuals. An exact two-sided sign test gives

p=2.62×10-4
(S87)

, rejecting equal positive and negative sign probabilities about zero. This does not contradict the near cancellation of the point-weighted residual in Eq. (141): the pooled statistic weights every radial measurement equally, whereas the galaxy-level statistic weights every galaxy equally.

The complete-catalog analysis therefore establishes a near-unity point-weighted aggregate normalization while retaining appreciable absolute radial residuals and a nonzero galaxy-weighted signed bias. These statistics describe different aspects of the catalog and should not be collapsed into a single goodness-of-fit measure. Section 4.5 next examines whether the residuals exhibit systematic velocity-dependent structure.

4.5 Residual structure and velocity dependence
The complete SPARC analysis shows that the remaining galaxy-level disagreement is frequently characterized by a coherent velocity-normalization offset superposed on an otherwise correlated rotation curve shape. To distinguish such normalization offsets from point-to-point residual scatter, let Bg denote the galaxy-level bias-to-scatter ratio and let

σΔvres,g
(S88)

denote the standard deviation of the signed velocity residuals about their mean for galaxy g:

Pg=|⟨Δvres⟩g|σΔvres,g.
(144)

The median value across the 175 galaxies is approximately Bg=1.14, with roughly 58% of galaxies satisfying Bg>1. Thus, for more than half of the catalog, the magnitude of the galaxy-level mean offset exceeds the residual scatter about that offset. This statistic is not a goodness-of-fit measure; it distinguishes coherent displacement of a rotation curve from residual variation about that displacement. The population-wide relation between curve shape and velocity normalization is examined separately in Sec. 4.9.

To characterize the velocity-normalization offset independent of absolute velocity scale, let Δlogv,g denote the galaxy-level logarithmic velocity offset,

Δlog⁡v,g=⟨log10⁡(vrealvo)⟩g,
(145)

evaluated only at radial measurements for which both velocities are positive. Across all 3,366 SPARC measurements, only two non-positive MQ realized velocities occur, both within a single galaxy, so the logarithmic diagnostic is not materially determined by the exclusion required by Eq. (145).

For the velocity-dependence test, each galaxy contributes one logarithmic offset from Eq. (145) and one galaxy-mean observed orbital velocity, ⟨vo⟩g. The correlation analysis therefore contains 175 galaxy-level observations rather than treating the 3,366 individual radial measurements as statistically independent. The use of ⟨vo⟩g is distinct from the independently published SPARC flat rotation velocity Vflat, which is tested separately in Sec. 4.7.

Across the complete 175-galaxy catalog, the logarithmic offset is strongly anticorrelated with ⟨vo⟩g. Using rP for the Pearson correlation coefficient and rS for the Spearman rank correlation coefficient, the corresponding values are rP=-0.503 with

p=1.26×10-12
(S89)

and rS=-0.426 with

p=4.13×10-9
(S90)

. Lower-velocity galaxies therefore tend to exhibit positive MQ-to-observed velocity-normalization offsets, whereas higher-velocity galaxies tend to exhibit negative offsets. The agreement of the Pearson and rank-based results shows that the detected association is not dependent solely on a linear representation of the galaxy-level data.

The same qualitative behavior is recovered when the direct signed velocity bias is used instead of the logarithmic velocity ratio. Its anticorrelation with ⟨vo⟩g shows that the observed trend is not specific to the logarithmic-ratio representation of Eq. (145), although both diagnostics necessarily contain the observed velocity and therefore do not by themselves establish the physical origin of the relation. A descriptive ordinary-least-squares regression of the galaxy-average velocities gives, with coefficients reported to the quoted precision,

⟨vreal⟩g≈10.17+0.937⟨vo⟩g,
(146)

where

⟨vreal⟩̂g
(S91)

denotes the fitted galaxy-mean MQ realized velocity. The regression has R2=0.95 to the quoted precision. The positive intercept and slope below unity compactly describe the population trend: toward the low-velocity end of the catalog the MQ galaxy-average velocity tends to lie above the observed average, whereas toward the high-velocity end it tends to lie below it. Because the two velocity variables span the same broad population scale, the large R2 characterizes co-variation and is not itself a measure of absolute predictive agreement.

The complete SPARC catalog therefore exhibits a statistically significant galaxy-level velocity-dependent normalization residual despite the near-zero point-weighted pooled bias reported in Sec. 4.4. The two results characterize different weightings: the pooled statistic measures signed cancellation across radial measurements, whereas the present analysis tests whether galaxy-level velocity-normalization offsets vary systematically with galaxy-mean rotational velocity.

The physical interpretation is deliberately deferred. SPARC spans broad ranges of galaxy morphology, luminosity, surface brightness, and kinematic regime [4], and those properties are not statistically independent of observational quality. The unstratified velocity dependence therefore does not by itself demonstrate either a velocity-dependent failure of the MQ galactic realization law or the need for an additional galaxy-dependent parameter. Section 4.6 tests the relation within SPARC’s independently assigned observational-quality classifications; Sec. 4.7 examines inclination, distance, and independent Vflat controls; and Sec. 4.8 evaluates additional galaxy-property dependencies. These tests determine whether the complete-catalog relation survives independently characterized observational and astrophysical structure.

4.6 SPARC observational quality stratification
SPARC assigns each galaxy an independently determined rotation curve quality classification (Q) [4]. These classifications predate the present MQ analysis and are therefore independent of the MQ residuals. SPARC defines (Q=1) for high-quality H I data or hybrid H()/H I rotation curves, (Q=2) for galaxies with minor asymmetries and/or lower quality H I data, and (Q=3) for galaxies exhibiting major asymmetries, strong non-circular motions, and/or offsets between the H I and stellar distributions [4]. SPARC regards the (Q=3) systems as unsuitable for detailed dynamical analyses, although their mass models are retained for completeness [4].

The resulting MQ velocity-normalization statistics at the fiducial SPARC stellar normalization are summarized in Table 10. The independently defined (Q=1) population provides the highest quality observational subset and is therefore particularly important for determining whether residual structure found in the complete catalog persists under improved observational conditions. This subset is not selected on the basis of MQ agreement.

Table 10. SPARC velocity-normalization statistics by observational quality class at the fiducial stellar normalization.

SPARC
quality

Ngal
(S92)

Geometric
MQ/observed ratio

Pearson velocity
correlation

Spearman velocity
correlation

Additional result

1
(S93)

99

= 1.0042

rP=-0.134,
p = 0.187
(S94)
rS=-0.136,
p = 0.178
(S95)

Mean signed bias = −0.59 km s-1; mean logarithmic offset consistent with zero

2
(S96)

64

= 1.146

rP=-0.582,p=4.5×10-7
(S97)
rS=-0.551,p=2.4×10-6
(S98)

Strong negative velocity dependence retained

3
(S99)

12

= 1.554

—

—

Mean positive logarithmic offset significant at approximately p = 7.8×10⁻⁴

The velocity-normalization residual changes strongly across SPARC’s independently assigned observational-quality classes. Table 11 gives the corresponding class-level inference from the galaxy-level logarithmic offset defined in Eq. (145). The (Q=1) population is centered close to unity in MQ-to-observed velocity normalization and exhibits no statistically significant dependence on galaxy-mean rotational velocity. The (Q=2) and (Q=3) populations show progressively larger positive normalization offsets, while (Q=2) retains a strong negative velocity dependence.

Table 11. SPARC velocity normalization and velocity dependence statistics by observational quality class at the fiducial stellar normalization.

Q
(S100)

Quality class

N
(S101)

Mean logarithmic offset

MQ/ observed ratio

Mean signed bias (km s-1)

Pearson velocity relation

Spearman velocity relation

One-sample offset test

Velocity-slope difference Q=1

1

High

99

0.00180

1.0042

-0.59
(S102)
rP=-0.134,p=0.187
(S103)
rS=-0.136,p=0.178
(S104)
p=0.773
(S105)

Reference

2

Medium

64

0.05928

1.146

+5.78

rP=-0.582,p≈4.5×10-7
(S106)
rS=-0.551,p≈2.4×10-6
(S107)
p≈2.22×10-5
(S108)
p<0.001
(S109)

3

Low

12

0.1915

1.554

—

—

—

p≈7.76×10-4
(S110)
p≈0.013
(S111)

The three quality classes also differ significantly in their logarithmic-offset distributions. A Kruskal-Wallis comparison gives a test statistic of 31.31 with p ≈ 1.6 × 10-7 [55]. The slope-interaction analysis likewise shows that the velocity dependence becomes significantly more negative as observational quality worsens, with the Q = 2 slope differing from Q = 1 at p < 0.001 and the Q = 3 interaction differing at p ≈ 0.013. Because Q = 3 contains only 12 galaxies, the latter result is treated as supporting evidence rather than a precise estimate of the low-quality velocity dependence.

Thus, at the fiducial SPARC stellar normalization, the complete-catalog velocity dependence is not uniformly expressed across observational quality. The (Q=1) population has an MQ-to-observed normalization approximately (0.42%) above unity and no statistically significant residual dependence on rotational velocity, while poorer-quality classes exhibit larger positive normalization offsets and, where sample size permits, stronger negative velocity dependence.

Stellar mass normalization provides a distinct source model sensitivity because the stellar disk and bulge terms enter the baryonic acceleration before the fixed MQ galactic realization law is applied. To test this dependence without introducing galaxy-specific fitting freedom, the complete SPARC calculation was repeated with one global 3.6 μm disk mass-to-light ratio imposed on every galaxy. Disk values of 0.42, 0.47, 0.50, and 0.60 M⊙/L⊙ span the published constant-normalization range summarized in Refs. [48–50]. To vary the overall stellar normalization while preserving the fiducial SPARC component convention, the bulge-to-disk ratio was held at 0.7/0.5 = 1.4, giving corresponding bulge values of 0.588, 0.658, 0.700, and 0.840. Distances, inclinations, gas contributions, observed velocities, and aMQ were held fixed. No mass-to-light ratio was optimized galaxy by galaxy.

For the 99 Q = 1 galaxies, the geometric MQ-to-observed velocity ratios are 0.9630, 0.9892, 1.0042, and 1.0507 for those four global normalizations, respectively. The corresponding Pearson correlations between galaxy logarithmic normalization offset and logarithmic mean observed velocity are rP = -0.212 (p = 0.035), -0.162 (p = 0.110), -0.133 (p = 0.188), and -0.049 (p = 0.631). Thus, the Q = 1 velocity normalization remains within approximately -3.70% to +5.07% of unity across this published global range, but the fiducial value 1.0042 is not invariant to the assumed stellar mass normalization. The residual velocity relation is statistically insignificant for the literature-centered 0.47-0.60 range, whereas the lowest 0.42 normalization gives a modest significant relation. This sensitivity belongs to the externally adopted baryonic source rather than MQ calibration because the galactic realization law and aMQ remain unchanged.

The fixed-global sweep bounds sensitivity to a common mass-to-light assumption but does not provide independently determined stellar masses for individual galaxies. A stronger source substitution is available for the 111 SPARC galaxies in the Duey et al. WISE stellar mass catalog [51]. Duey et al. infer color-dependent W1 stellar masses from stellar population models independently of the MQ residual. For each overlapping galaxy, the published WISE total stellar mass was divided by the SPARC 3.6 μm luminosity to obtain an external galaxy-level stellar normalization. That normalization was applied to the unit-mass-to-light SPARC disk and bulge velocity templates in Eq. (130), while the gas contribution, distance, inclination, observed velocity, and fixed aMQ were left unchanged. The resulting baryonic acceleration was propagated through Eq. (137) as in the fiducial calculation. No stellar normalization was optimized against vo or the MQ residual.

This construction changes the stellar source rather than merely introducing stellar mass as a covariate, although it continues to use the SPARC radial light templates. Table 12 summarizes the resulting WISE stellar-mass source-normalization test using identical galaxy subsets for the WISE and fiducial SPARC calculations. Across all 111 overlapping WISE galaxies, comprising 2,372 radial measurements, the galaxy-level geometric MQ-to-observed ratio is 1.1237, compared with 1.0429 for the identical galaxies under the fiducial SPARC source. The (Q=1) overlap contains 71 galaxies and 1,595 radial measurements; the WISE-source ratio is 1.0697 and the mean galaxy signed bias is +15.55 km s-1, whereas their same-galaxy fiducial-source ratio is 0.9863. Restricting further to 30∘≤i≤80∘leaves 60 (Q=1) galaxies and gives 1.0854, compared with 0.9915 under the fiducial source. The independent stellar-mass substitution therefore moves the (Q=1) overlap approximately 7%–9% above the MQ curve relative to its fiducial-source result, demonstrating material stellar-source sensitivity without calibrating MQ to the rotation-curve residual.

Table 12. Independent WISE stellar mass source-normalization test. The fiducial comparison is recalculated on the identical galaxies in each row.

Sample

Galaxies

Radial points

WISE-source geometric MQ/observed

Same-galaxy fiducial MQ/observed

Interpretation

WISE overlap, all quality classes

111

2,372

1.1237

1.0429

Independent stellar normalization shifts the source upward

WISE overlap, Q = 1

71

1,595

1.0697

0.9863

High-quality overlap moves about 7% MQ-high

WISE overlap, Q = 1, 30° ≤ i ≤ 80°

60

1,407

1.0854

0.9915

Shift persists after inclination restriction

Table 13 tests whether this source sensitivity persists when galaxies are selected without reference to their MQ residual and the stellar contribution is progressively reduced by gas richness.

Table 13. Residual-blind gas-rich source-stability control within the Q = 1 WISE overlap. The sample is restricted to 30° ≤ i ≤ 80° before any MQ residual is examined.

Sample

Galaxies

Fiducial geometric MQ/observed

WISE-source geometric MQ/observed

Relative source shift

Wilcoxon p for WISE-fiducial change

All Q = 1, 30° ≤ i ≤ 80°

60

0.9915

1.0854

+9.46%

5.56 × 10^-7

Gas-rich: MHI/L[3.6] ≥ 0.5

19

0.9854

0.9886

+0.33%

0.9217

Gas-dominated: MHI/L[3.6] ≥ 1.0

8

1.0086

1.0034

-0.51%

0.3828

The purpose of Table 13 is to isolate stellar source sensitivity without selecting galaxies on the basis of MQ agreement. The selection uses only published SPARC quality, inclination, H I mass, and 3.6 μm luminosity, so gas richness is fixed on the source side before the MQ residual is evaluated. SPARC itself shows that the transition from star-dominated to gas-dominated systems tracks gas fraction and luminosity [4]. For MHI/L3.6 expressed in solar units, the stricter MHI/L3.6 ≥ 1.0 criterion is deliberately conservative: with the SPARC atomic-gas scaling Mgas = 1.33 MHI and even Υ* = 0.60 M⊙/L⊙, it implies Mgas/M* ≥ 2.22.

The source substitution has a much smaller effect in these gas-rich controls than in the complete inclination-controlled WISE overlap. Across all 60 (Q=1) inclination-controlled WISE galaxies, replacing the fiducial stellar source by the independent WISE source changes the geometric normalization by (+9.46%). In the 19-galaxy gas-rich subset the change is (+0.33%), from 0.9854 to 0.9886, with the paired galaxy-level change statistically consistent with zero (Wilcoxon (p=0.9217)); in the stricter eight-galaxy gas-dominated subset the change is (-0.51%), from 1.0086 to 1.0034 ((p=0.3828)). The broad WISE-overlap displacement is therefore not universal across baryonic regime and becomes small where the stellar contribution is dynamically subdominant.

This control narrows the absolute-normalization concern but does not remove the stellar source degeneracy in star-dominated galaxies or establish that either stellar mass scale is uniquely correct. It identifies a source-insensitive observational regime in which the fixed MQ galactic realization law remains near unity under two independently constructed stellar normalizations. Duey et al. report that WISE stellar mass uncertainties are largest for the bluest systems and that high-mass galaxies additionally depend on bulge-disk separation [51]. A precision claim for stellar-dominated systems is therefore withheld unless the radial stellar source itself is independently reconstructed. The DiskMass approach illustrates a complementary dynamical route by deriving stellar surface density from vertical stellar velocity dispersions after accounting for the gas contribution [56].

The WISE substitution is an independent stellar-normalization test but not a fully resolved stellar surface-density experiment. Duey et al. provide integrated color-dependent stellar masses rather than a radially varying stellar population mass surface-density map for each SPARC galaxy, so the same galaxy-level normalization must scale the existing SPARC stellar light templates. A single total stellar mass also does not independently re-establish the disk-bulge decomposition in systems with significant bulges. Accordingly, the test removes the universal 3.6 μm normalization from the overlapping sample but retains source-side uncertainty from radial mass-to-light gradients, disk-bulge decomposition, and stellar population covariance. A complementary independent route exists in the LITTLE THINGS literature: Zhang et al. [57] derive radial stellar mass surface-density distributions from multiband SED fits, while Oh et al. [5] construct radial H I and Spitzer 3.6 μm stellar surface-density profiles using the kinematic disk geometry. Such products provide the radial source information required to extend a galaxy-integrated stellar-normalization test to radial-shape reconstruction.

Observational quality is itself associated with galaxy velocity. Because (Q) is an ordinal classification rather than an interval-scaled quantity, the primary test uses a rank correlation. Across the complete 175-galaxy catalog,

rS(Q,log⁡⟨vo⟩g)=-0.4043,p=2.89×10-8.
(147)

4.7 Inclination, distance, and observational-systematic controls

The velocity-dependent residual identified in Sec. 4.5 and its strong dependence on SPARC observational quality established in Sec. 4.6 motivate direct tests of the principal observational quantities capable of producing coherent changes in a reconstructed rotation curve. Distance and inclination enter the SPARC analysis differently. The adopted distance determines the physical radial scale and affects the inferred baryonic contribution, whereas inclination determines the conversion of measured line-of-sight kinematics to circular velocity [4]. Both can therefore alter the comparison between the MQ prediction and the published rotation curve, but through physically distinct pathways.

Distance does not materially track the MQ-to-observed velocity-normalization offset. Across the complete SPARC catalog, galaxy-level signed bias versus distance gives Pearson rP=-0.03 with p=0.69 and Spearman rS=+0.03 with p=0.65. The logarithmic velocity offset defined in Eq. (145) likewise shows negligible dependence on distance, with rP=-0.064 and p=0.40. Thus, no statistically significant galaxy-level association with distance is detected. This argues against a simple distance-dependent scaling error as the origin of the complete-catalog velocity trend, although it does not eliminate distance uncertainty from individual baryonic reconstructions.

Let i denote the published disk inclination and vlos the line-of-sight rotational component before inclination deprojection. Under the circular thin-disk geometry used for this diagnostic, the observed orbital velocity is

vo=vlossin⁡i.
(148)

Let σi denote the quoted inclination uncertainty in radians and define

σvo,i
(S112)

as the corresponding first-order propagated uncertainty in vo. Linear propagation of Eq. (148) gives

σvo,ivo=|cot⁡i|σi.
(149)

SPARC notes that inclination corrections vary with sini, becoming increasingly consequential toward face-on orientations, and low-inclination galaxies are commonly excluded from precision dynamical analyses for this reason [4,58]. The strong low-inclination sensitivity makes inclination uncertainty a natural candidate for a coherent velocity-normalization offset.

No corresponding population-wide dependence is detected. Across the complete catalog, logarithmic velocity offset versus published inclination gives rP=-0.039 with p=0.608 and rS=+0.027 with p=0.724. Dependence on the quoted inclination uncertainty is likewise weak and statistically non-significant, with rP=0.141 and p=0.062, while the corresponding relation with the inclination-induced fractional velocity uncertainty gives rP=0.104 and p=0.171. These results do not establish that inclination is irrelevant to individual galaxies; they show that no simple catalog-wide inclination dependence reproduces the residual trend.

A complementary diagnostic asks what single inclination rescaling would be required if an individual galaxy’s complete MQ-to-observed velocity-normalization offset were attributed to inclination. Let ireq denote the required inclination. Because Eq. (145) defines the galaxy geometric MQ-to-observed ratio as

10Δlogv,g
(S113)

, Eq. (148) gives

sin⁡ireq=sin⁡i​10-Δlog⁡v,g.
(150)

This construction is diagnostic and does not replace the published SPARC inclination or constitute a refit of the galaxy kinematics. SPARC inclinations are adopted from the underlying kinematic analyses, commonly including tilted-ring modeling of H I velocity fields, and can exhibit radial variation associated with warps [4]. Equation (150) instead asks whether a single coherent inclination rescaling of the magnitude implied by the galaxy-average normalization offset could reconcile that offset.

Formal geometric solutions exist for 157 of the 175 galaxies; 18 do not admit a solution under this construction. Let Δireq=ireq-i denote the required inclination shift. Across the complete catalog, the median magnitude is |Δireq|=9.15∘, corresponding to 2.07σi. Of the galaxies, 31% require shifts within 1σi, 50% within 2σi, and 59% within 3σi. The required shifts are asymmetric in sign, with more systems requiring reductions rather than increases in the adopted inclination. Inclination changes can therefore formally reconcile the normalization of many individual systems, but their required magnitudes and population distribution do not naturally reproduce the systematic velocity dependence of the complete catalog.

A more restrictive diagnostic combines the independently assigned highest quality observational class with an analysis-defined intermediate inclination range. Requiring Q=1 and 30∘≤i≤80∘ leaves 75 galaxies. The lower boundary follows the established SPARC practice of excluding face-on systems with large and uncertain inclination corrections [4,58]; the upper boundary removes near-edge-on configurations as an additional observational control rather than as a SPARC quality criterion. The geometric MQ-to-observed velocity ratio for this subset is 1.0078, corresponding to a +0.78% normalization offset. The correlation between logarithmic velocity offset and galaxy-mean observed orbital velocity is rP=-0.220 with p=0.058, which does not reach the conventional 5% significance threshold.

Within the same 75-galaxy subset, no significant residual association with inclination is detected: rP=0.046 with p=0.697 and rS=-0.015 with p=0.900. Restricting the analysis simultaneously by observational quality and inclination therefore does not recover a statistically significant inclination dependence. The near-threshold velocity result at p=0.058 is retained without being characterized as evidence for the absence of an underlying association.

The velocity dependence is further tested using SPARC’s published flat rotation velocity Vflat as an alternative characteristic velocity [4,58]. Here Vflat denotes the mean velocity associated with the flat portion of the published rotation curve. It is defined independently of the present MQ analysis, although it derives from the same underlying SPARC kinematic data. Its use therefore tests sensitivity to the definition of characteristic rotational velocity rather than supplying a statistically independent rotation curve dataset.

Among galaxies satisfying Q=1, 30∘≤i≤80∘, and having a valid Vflat, 65 systems remain. Their geometric MQ-to-observed normalization ratio is 0.9948, corresponding to -0.52% from unity. Using galaxy-mean observed orbital velocity as the characteristic velocity gives rP=-0.141 with p=0.263. Replacing that measure with the published Vflat gives

rP​(Δlog⁡v,g,Vflat)=-0.084,p=0.440.
(151)

Thus, the velocity-dependent residual detected in the complete catalog is not statistically significant when the analysis is restricted to this highest quality, inclination-controlled subset and characterized using SPARC’s published flat rotation velocity. This is a robustness test against characteristic-velocity definition rather than an independent replication.

These controls must also be interpreted within the SPARC error model. The published pointwise velocity uncertainties combine formal fitting error with approaching-receding-side asymmetry and thereby include local non-circular motions and large-scale kinematic asymmetry, but they explicitly exclude the systematic uncertainty associated with the adopted inclination [4]. For a disk without a strong warp, an inclination change can coherently alter rotation curve normalization without substantially changing its radial shape [4]. Consequently, the random-error χ2 diagnostics defined in Sec. 4.1 cannot be interpreted as complete likelihood statistics incorporating all orientation and reconstruction uncertainties.

Taken together, the controls do not identify either distance or inclination as a simple catalog-wide origin of the velocity-dependent residual. In the Q=1, 30∘≤i≤80∘ subset, the normalization remains near unity and neither inclination nor the residual velocity relation reaches conventional statistical significance. Among the 65 systems also possessing Vflat, the normalization is 0.52% below unity and the velocity correlation remains non-significant. Section 4.8 next tests whether the remaining residual structure is associated with luminosity, surface brightness, morphology, gas content, radial coverage, or other galaxy properties.

4.8 Galaxy-property controls
The absence of a statistically significant velocity-dependent residual within the independently defined Q=1 population permits a separate test of whether the remaining velocity-normalization offsets depend on other major galaxy properties. The analysis therefore concentrates on the 99 highest quality SPARC galaxies and evaluates the galaxy-level logarithmic velocity offset defined in Eq. (145) against luminosity, surface brightness, morphology, gas richness, and radial coverage. These properties span physically distinct characteristics of the SPARC population [4].

Table 14 summarizes the direct Pearson correlations. Here, L3.6 denotes the SPARC 3.6 μm luminosity, SBeff the effective surface brightness, SBdisk the central disk surface brightness, T the SPARC Hubble-type classification, Rmax the maximum sampled rotation curve radius, Rdisk the stellar-disk scale length, and MHI the neutral atomic-hydrogen mass. The ratio

Rmax/Rdisk
(S114)

characterizes normalized radial coverage, while

MHI/L3.6
(S115)

provides a measure of relative atomic-gas content.

Table 14. Galaxy-property correlations with the logarithmic MQ-to-observed velocity offset for the Q=1 population.

Galaxy property

Statistic tested

Physical characteristic

Pearson rP

p
(S116)

Stellar luminosity

log⁡L3.6
(S117)

Stellar luminosity

+0.045
(S118)

0.659

Effective surface brightness

log⁡SBeff
(S119)

Global surface brightness structure

+0.006
(S120)

0.956

Central disk surface brightness

log⁡SBdisk
(S121)

Central disk surface brightness structure

-0.004
(S122)

0.968

Hubble type

T
(S123)

Broad galaxy morphology

+0.095
(S124)

0.347

Normalized radial coverage

Rmax/Rdisk
(S125)

Radial extent of the measured rotation curve

+0.037
(S126)

0.714

Atomic-gas richness

MHI/L3.6
(S127)

Relative atomic-gas content

+0.008
(S128)

0.941

None of the tested properties exhibits a statistically significant direct linear association with the logarithmic MQ-to-observed velocity offset. In particular, the residual structure within the highest quality population shows no evidence in these tests of preferential association with low luminosity, low surface brightness, late morphology, high gas richness, or unusually extended radial coverage. Because T is an ordinal Hubble-type classification rather than a continuously measured physical variable, its Pearson coefficient is retained only as a broad linear-code diagnostic.

Strong astrophysical correlations are nevertheless present within SPARC itself. The published SPARC analysis shows that baryonic-to-observed velocity behavior varies systematically with luminosity and surface brightness, with higher-mass, higher-surface-brightness galaxies generally having larger baryonic contributions and lower-mass, lower-surface-brightness systems being more submaximal [4]. The absence of analogous direct linear trends in the Q=1 MQ residuals therefore argues against a simple dependence of the remaining normalization offsets on these galaxy regimes.

The result is also consistent with broader SPARC residual studies. Stiskalek and Desmond find that inclusion of additional available galaxy properties does not significantly tighten the radial acceleration relation and that no significant partially independent secondary correlation is required to describe the radial galaxy dynamics [59]. This does not provide evidence for MQ specifically but places the absence of strong secondary residual correlations in the present highest quality subset within the broader SPARC phenomenology.

A simultaneous galaxy-level partial regression analysis was therefore performed across the complete 175-galaxy catalog rather than applying these controls individually. The dependent quantity is the galaxy-level logarithmic MQ-to-observed velocity offset of Eq. (145). The explanatory variables are galaxy-mean observed rotational velocity, SPARC quality class treated categorically with Q=1 as the reference, inclination, distance, 3.6 μm luminosity, effective surface brightness, and atomic-gas richness

MHI/L3.6
(S129)

. Positive scale variables were transformed logarithmically, continuous predictors were standardized, and heteroskedasticity-consistent HC3 standard errors were used. This construction keeps the galaxy, rather than the individual radial point, as the statistical unit and follows the general SPARC practice of treating correlated galaxy parameters jointly rather than as independent pointwise errors [60].

The complete-catalog velocity dependence remains statistically significant under this simultaneous control. In a velocity-only model, the standardized velocity coefficient is -0.0501±0.00785 (

p=1.75×10-10
(S130)

; adjusted R2=0.249). Adding observational quality reduces its magnitude to -0.0365±0.00792 (

p=4.01×10-6
(S131)

; adjusted R2=0.357). In the full simultaneous model, the conditional velocity coefficient is -0.2389±0.01536 (

p=1.55×10-54
(S132)

), the adjusted R2 is 0.775, and the partial R2 associated with velocity is 0.688. Repeating the same simultaneous control with the direct mean signed velocity residual rather than the logarithmic ratio also retains a significant velocity term (

p=5.41×10-14
(S133)

). Thus, the conditional association is not specific to the logarithmic-ratio representation of Eq. (145), although both residual constructions contain the observed velocity and therefore do not by themselves establish a causal velocity dependence.

The conditional coefficients cannot be interpreted as independent causal effects. Galaxy-mean rotational velocity and stellar luminosity are exceptionally collinear in this catalog, with Pearson rP=0.931, while the variance-inflation factors are 11.6 for velocity and 16.9 for luminosity in the full model. The large negative velocity coefficient paired with a large positive luminosity coefficient is therefore a suppression-sensitive decomposition of two tightly related galaxy-scale variables. This motivates the matched-control test below, in which luminosity is replaced by independently estimated stellar mass and galaxies are paired at nearly fixed stellar mass rather than separated statistically by a multivariate regression.

The complete catalog trend is therefore retained as an empirical constraint rather than interpreted as a clean global success or failure of the MQ galactic realization law. Its absence as a marginal relation in the independently defined Q=1 subset cannot by itself establish that the lower quality residuals are observational. The matched-control construction in Table 15 asks whether the velocity ordering survives among independently high-quality galaxies after directly restricting the stellar mass difference between paired systems.

Table 15. Simultaneous complete-sample partial regression diagnostics for the galaxy-level logarithmic MQ-to-observed velocity offset.

Model / predictor

Standardized coefficient

HC3 robust SE

p value

Diagnostic

Velocity only

−0.0501

0.00785

1.75×10-10
(S134)

Adjusted R2=0.249

Velocity + SPARC quality

−0.0365

0.00792

4.01×10-6
(S135)

Adjusted R2=0.357

Full simultaneous: velocity

−0.2389

0.01536

1.55×10-54
(S136)

Partial R2=0.688

Full: Q = 2 indicator

+0.0247

0.00818

0.00254

Q=1 reference

Full: Q = 3 indicator

+0.0398

0.02435

0.103

Q=1 reference

Full: inclination

+0.00623

0.00415

0.133

VIF = 1.13

Full: distance

−0.00467

0.00484

0.334

VIF = 2.45

Full: 3.6 μm luminosity

+0.2078

0.01861

5.74×10-29
(S137)

VIF = 16.9

Full: effective surface brightness

+0.0280

0.00756

2.02×10-4
(S138)

VIF = 5.41

Full: atomic-gas richness

+0.0282

0.01105

0.0107

VIF = 4.69

Continuous predictors are standardized after logarithmic transformation of positive scale variables; coefficients are conditional rather than causal. The full model has adjusted R2=0.775. Strong velocity-luminosity collinearity (rP=0.931) requires the conditional coefficients to be interpreted jointly.

To break the velocity-luminosity degeneracy directly, the analysis was repeated with the independent WISE stellar mass catalog of Duey et al. [51]. That work derives color-dependent W1 stellar masses for SPARC galaxies from stellar population models and reports close agreement with independent Spitzer 3.6 μm and multiband stellar mass determinations. The matched test uses only SPARC Q=1 galaxies with 30∘≤i≤80∘ that have a published WISE stellar mass. Galaxies are matched without using the MQ residual: pair eligibility requires an absolute stellar mass difference no larger than 0.10 dex, an absolute inclination difference no larger than 15∘, and a separation in log10 mean observed velocity of at least 0.10 dex. A maximum-cardinality non-overlapping matching then favors smaller stellar mass and inclination differences and larger velocity separation. Each galaxy therefore enters at most one pair.

The primary construction produces 10 independent galaxy pairs. Their mean absolute stellar mass mismatch is 0.0437 dex, while the mean separation in log10 observed velocity is 0.1712 dex. Orienting each pair from its lower-velocity to higher-velocity member, the mean change in the logarithmic MQ-to-observed velocity offset is -0.1177 dex and the median change is -0.1034 dex. Nine of the ten pairs have the higher-velocity member at the lower MQ-to-observed normalization. A two-sided Wilcoxon signed-rank test gives p=0.0039. The conditional velocity ordering therefore survives when the compared systems are independently high quality, inclination controlled and nearly matched in stellar mass.

Table 16. Independent-WISE-stellar mass matched-pair sensitivity test for the Q = 1, 30° ≤ i ≤ 80° SPARC population.

Maximum stellar mass mismatch (dex)

Non-overlapping pairs

Mean |stellar mass mismatch| (dex)

Mean log10 velocity separation (dex)

Mean higher-V minus lower-V log MQ/ observed offset (dex)

Wilcoxon p

0.05

7

0.0293

0.1569

-0.1002

0.0313

0.08

9

0.0386

0.1706

-0.1181

0.0078

0.10 (primary)

10

0.0437

0.1712

-0.1177

0.0039

0.12

12

0.0587

0.1759

-0.1228

0.0010

0.15

13

0.0644

0.1694

-0.1221

0.0005

Table 16 shows that the ordering is stable to the stellar mass matching tolerance. Tightening the maximum mismatch to 0.05 dex leaves seven non-overlapping pairs and the same negative conditional ordering; relaxing it through 0.15 dex increases the matched sample while preserving a mean offset difference near -0.12 dex. These sensitivity rows are not independent statistical tests because the pair sets overlap; they show that the primary result is not created by the particular 0.10 dex threshold.

This matched analysis resolves the principal interpretive ambiguity of the simultaneous regression but does not resolve the empirical discrepancy in favor of MQ. Independent stellar mass matching shows that the velocity dependence is not solely a suppression artifact produced by using luminosity and velocity as highly collinear regressors. The existing SPARC data instead retain a conditional scale ordering under direct stellar mass control. Because the WISE masses are integrated stellar masses rather than independently reconstructed radial stellar surface-density profiles, and because the kinematic target remains the SPARC rotation curve, this is not an independent rotation curve population test.

A stronger source-recomputed test inserts the external Duey et al. [51] stellar mass into the baryonic source before Eq. (137) is evaluated and repeats the residual velocity analyses using the resulting WISE-source MQ prediction rather than the fiducial SPARC source. The calculation uses all 111 WISE-overlap galaxies while retaining the fixed aMQ and every published SPARC distance, inclination, gas contribution, and observed velocity. No MQ residual enters determination of the stellar normalization.

The broad marginal velocity relation changes materially under this source substitution. Across all 111 overlapping galaxies, the Pearson correlation between log10 mean observed velocity and the WISE-source logarithmic MQ-to-observed offset is rP=-0.1182 (p=0.2166), with Spearman rS=-0.0562 (p=0.5579). For the 71 Q=1 galaxies it reverses sign, rP=+0.2140 (p=0.0732), with rS=+0.1549 (p=0.1972). Restricting further to 30∘≤i≤80∘ leaves 60 galaxies and gives rP=+0.2627 (p=0.0426) but rS=+0.2176 (p=0.0950). Thus, the negative marginal residual velocity trend is not invariant to an independently supplied stellar source.

The conditional ordering nevertheless survives when the matched-pair calculation is recomputed from the WISE-source predictions themselves. Using the same residual-blind 0.10-dex stellar mass, 15∘ inclination, and 0.10-dex velocity-separation criteria produces the same 10 non-overlapping pairs. All 10 higher-velocity members now have the lower WISE-source MQ-to-observed normalization. The mean higher-minus-lower offset is -0.12685 dex and the median is -0.11768 dex; both the exact two-sided Wilcoxon signed-rank test and exact two-sided sign test give p=0.001953.

Table 17 summarizes the source-recomputed residual tests. A complementary 60-galaxy Q=1, inclination-controlled regression uses the WISE-source logarithmic offset as the dependent variable and simultaneously controls log10 WISE stellar mass, inclination, log10 distance, log10 effective surface brightness, and log10 atomic-gas richness. Let βV denote the standardized conditional velocity coefficient in this regression. After standardization, βV=-0.1697 with HC3 robust SE = 0.0265 (

p=4.16×10-8
(S139)

) and partial R2=0.468; a 10,000-resample galaxy bootstrap gives a 95% coefficient interval of [-0.2148,-0.1175]. This result is statistically stable under the stated controls but remains suppression-sensitive because log10 velocity and log10 WISE stellar mass are strongly correlated (rP=0.9578), with VIFs 16.8 and 20.8, respectively. It therefore establishes persistence of conditional ordering within this SPARC subset, not a causal velocity correction to Eq. (137).

Table 17. WISE-source residual diagnostics. The WISE stellar mass is inserted into the baryonic source before the unchanged MQ galactic realization law is evaluated.

Test

Sample

Principal result

Significance

Interpretation

WISE-source marginal

111 galaxies

rP = -0.1182; rS = -0.0562

Pearson p = 0.2166; Spearman p = 0.5579

No significant negative marginal velocity relation

WISE-source, Q = 1

71 galaxies

rP = +0.2140; rS = +0.1549

Pearson p = 0.0732; Spearman p = 0.1972

Negative marginal relation reverses sign

WISE-source, Q = 1, 30° ≤ i ≤ 80°

60 galaxies

rP = +0.2627; rS = +0.2176

Pearson p = 0.0426; Spearman p = 0.0950

Marginal result is source-sensitive and not robustly monotonic

WISE-source matched pairs

10 pairs

Mean higher-minus-lower = -0.12685 dex; median = -0.11768 dex; 10/10 negative

Wilcoxon p = 0.001953; sign p = 0.001953

Conditional velocity ordering persists

WISE-source multivariate control

60 galaxies

betaV = -0.1697; HC3 SE = 0.0265; partial R² = 0.468

HC3 p = 4.16 × 10^-8; bootstrap 95% [-0.2148, -0.1175]

Conditional term persists, but velocity-mass collinearity remains high

Taken together, these tests sharpen the empirical ruling. The complete-SPARC negative marginal velocity trend is source model dependent: it becomes non-significant or reverses sign when the independent WISE stellar normalization is propagated through the prediction, while the residual-blind stellar mass-matched ordering and simultaneous WISE-source conditional velocity term persist. The SPARC evidence therefore neither establishes universal velocity independence nor isolates a velocity-dependent failure of the fixed MQ galactic realization law. Because the WISE calculations remain within the SPARC population and reuse its kinematic target, they are internal source and robustness controls rather than an independent population test. Section 4.10 supplies that independent out-of-SPARC population test.

4.9 Population-wide curve shape and normalization behavior
The preceding analyses show that the remaining SPARC discrepancies are characterized by two related but distinct statistical behaviors: radial co-variation between the MQ realized and observed rotation curves and displacement of their overall velocity normalization. Table 18 summarizes these behaviors across the three independently assigned SPARC observational-quality classes.

The absolute residual scale increases with worsening observational quality, most clearly in the Q=3 population, which has substantially larger median RMS and MAE values than the two higher-quality classes. The radial correlation statistics do not follow the same monotonic progression. In particular, the Q=2 population has a slightly higher median Pearson correlation than Q=1 despite having a substantially larger normalization displacement.

Table 18. Curve-shape and normalization diagnostics by SPARC observational-quality class.

QSPARC
(S140)

Quality class

N
(S141)

Median RMS residual (km s-1)

Median MAE (km s-1)

Median MQ–observed Pearson rP

Median Pg

Fraction with Pg>1

1

High

99

11.72

10.67

0.955

0.880

44.4%

2

Medium

64

13.30

10.91

0.970

1.589

70.3%

3

Low

12

27.92

25.75

0.909

2.657

100%

This distinction is quantified by the bias-to-scatter ratio Bg defined in Eq. (144). A value Bg>1 means that the magnitude of a galaxy’s mean signed velocity residual exceeds the standard deviation of the residuals about that mean. The median Bg rises from 0.880 for Q=1 to 1.589 for Q=2 and 2.657 for Q=3, while the fraction of galaxies satisfying Bg>1 increases from 44.4% to 70.3% and then to 100%. The principal quality-dependent change is therefore a progressively stronger coherent velocity displacement relative to the residual scatter within individual galaxies.

The Q=2 population illustrates why these quantities must be interpreted jointly. Its median Pearson coefficient of 0.970 indicates exceptionally strong radial co-variation between the MQ realized and observed curves, yet its median Bg is approximately 1.81 times that of the Q=1 population. Strong radial co-variation can therefore coexist with substantial normalization displacement. Pearson correlation measures co-variation rather than equality of velocity scale, while RMS and MAE measure residual magnitude but do not by themselves distinguish coherent normalization displacement from changes in radial curve structure [61].

The Q=3 population differs in both respects: it has the largest median RMS and MAE, the lowest median Pearson correlation, and the largest median Bg. Interpretation of this class must nevertheless remain limited by its sample of only 12 galaxies and by SPARC’s independent classification of these systems as having the poorest rotation curve quality, including major asymmetries, strong non-circular motions, and/or offsets between the H I and stellar distributions [4]. The 100% incidence of Bg>1 therefore characterizes this specific 12-galaxy subset and should not be generalized to a broader galaxy population.

Across the SPARC catalog, strong radial co-variation can persist even when velocity normalizations differ substantially, while decreasing observational quality is accompanied by increasingly coherent galaxy-level displacement relative to within-galaxy residual scatter. Correlation, residual magnitude, and normalization displacement therefore measure distinct aspects of predictive behavior and must be evaluated jointly rather than treated as interchangeable measures of agreement.

4.10 Combined interpretation of the Milky Way and SPARC rotation curve tests

The Milky Way and SPARC analyses provide complementary constraints on the same fixed MQ galactic realization law. In both cases, the calculation begins with an independently published baryonic reconstruction and applies the realization relation specified in Sec. 4.2 with zero galaxy-specific fitted MQ parameters: no empirical interpolation function, dark matter halo profile, or MQ acceleration scale is adjusted to either rotation curve dataset. The analyses differ in their observational data and baryonic reconstruction methodology rather than in the realization law itself. Their evidentiary roles also differ: SPARC tests the fixed law across a large heterogeneous population, whereas the principal Milky Way test combines the independent Eilers kinematic target with the RAR/MOND-independent McMillan baryonic source.

The uncertainty structures of the two tests are likewise different. For SPARC, the galaxy-cluster bootstrap treats each rotation curve as the sampling unit and preserves the radial clustering within each galaxy; the resulting intervals leave both the pooled signed bias statistically compatible with zero and the Q=1 normalization compatible with unity. For the Milky Way, a galaxy-level bootstrap is unavailable because only one system is observed. The pointwise Eilers/Jiao error model and the published McMillan source covariance therefore address different uncertainty components. Under diagonal observational errors alone the Milky Way residual is formally inconsistent with that incomplete error model, but propagation of the independently published McMillan baryonic-normalization covariance changes the generalized statistic to χ2=39.57 for 38 degrees of freedom with p=0.400; the inner and outer subsets likewise give p=0.423 and p=0.910, respectively. The remaining limitation is that no full Eilers point-to-point covariance is supplied and the complete McMillan source posterior is not propagated jointly. The equicorrelation and radial-kernel calculations of Sec. 4.3 therefore remain observational-covariance sensitivity tests rather than substitutes for a measured Eilers covariance. Section 5.9 specifies the additional covariance-ready observational products required for a fully propagated residual covariance and generalized likelihood.

The Milky Way provides a dedicated RAR/MOND-independent Galactic test using the baryonic density model of McMillan [3] and the independently derived circular velocity measurements of Eilers et al. [2]. Only McMillan’s baryonic components enter the MQ source; the NFW halo is excluded. The earlier McGaugh [46] calculation is retained only as a source model sensitivity comparison. SPARC provides the broader population test across 175 galaxies and 3,366 radial measurements [4]. The independent LITTLE THINGS population then tests whether the population normalization and SPARC residual velocity ordering persist outside the SPARC catalog. Their principal results are summarized in Table 19.

Table 19. Principal Milky Way and SPARC tests of the fixed MQ local realization law.

Test

Sample

Principal MQ result

Newtonian/control comparison

Interpretation

Milky Way, complete

38 radii, 5.27–24.82 kpc

RMS 8.90 km s-1; MAE 7.43 km s-1; mean signed residual −3.07 km s-1

Newtonian RMS 69.03 km s-1; McMillan-source-covariance generalized χ² = 39.57, p = 0.400

RAR/MOND-independent McMillan source; full-range RMS reduction 87.11%; not rejected under the propagated published-source covariance

Milky Way, inner

30 radii, 5.27–19.71 kpc

RMS 6.88 km s-1; mean signed residual −6.19 km s-1

Newtonian RMS 66.69 km s-1; source-covariance generalized χ² = 30.85, p = 0.423

Highest-precision radial region; RMS reduction 89.68%

SPARC, complete

175 galaxies; 3,366 radial points

Pooled signed residual +0.02644 km s-1; <vreal>/<vo> = 1.0001888; pooled RMS 22.82 km s-1; pooled MAE 16.20 km s-1

No galaxy-specific MQ fitting parameters; galaxy-cluster bootstrap preserves within-galaxy clustering

Aggregate normalization is nearly centered, but substantial radial residuals remain

SPARC, highest quality

Q = 1, 99 galaxies

Geometric MQ/observed ratio 1.0042; rP = −0.134, p = 0.187; rS = −0.136, p = 0.178

Independently assigned SPARC quality classification

Velocity-dependent normalization residual is not statistically significant in the highest quality population

Simultaneous partial regression, complete SPARC

175 galaxies

Conditional standardized velocity coefficient −0.2389 ± 0.0154, p = 1.55 × 10⁻⁵⁴

Joint measured-covariate control; strong velocity-luminosity collinearity

Conditional ordering survives multivariate control, but its physical attribution is not identifiable from the regression

LITTLE THINGS, independent population

15 non-SPARC galaxies; 288 figure-extracted radial comparisons

Galaxy-weighted geometric MQ/observed ratio 0.9907; 95% galaxy-block interval 0.868–1.130; zero-offset test p = 0.895

No galaxy-specific MQ fitting parameters; independent SED stellar normalization; primary rP = −0.282, p = 0.308; stellar mass-and-inclination partial rP = −0.257, p = 0.355

Population normalization is statistically compatible with unity; the SPARC conditional velocity ordering is not independently replicated, while substantial pointwise and galaxy-level scatter remains

The RAR/MOND-independent Milky Way result is strongest over the inner 30 measurements, where the Eilers uncertainties are smallest and the MQ realized curve substantially reduces the baryonic Newtonian discrepancy. The larger residuals at the outer eight radii are accompanied by substantially larger observational uncertainties, as detailed in Sec. 4.3. The Milky Way test therefore evaluates the fixed MQ realization of the published McMillan Galactic baryonic profile against an observational velocity curve obtained through a separate kinematic analysis. Its current uncertainty result is compatibility under the propagated published source covariance, not covariance-complete precision validation.

The complete SPARC catalog is nearly centered in point-weighted signed velocity normalization, but the 0.0189% aggregate difference must not be interpreted as pointwise agreement at that level. The pooled signed residual reflects cancellation among radial and galaxy-level offsets, whereas the pooled RMS and MAE remain much larger. Sections 4.5–4.9 therefore test the structure of those residuals rather than relying on the aggregate mean.

The most prominent complete-catalog residual is its dependence on characteristic rotational velocity. That relation changes materially under SPARC’s independently assigned observational-quality classification and under baryonic source substitution. The principal robustness and population controls are summarized in Table 20.

Table 20. Principal SPARC robustness, observational-control, source-sensitivity, and independent-population results.

Analysis subset or control

N
(S142)

MQ/observed normalization

Velocity relation

Principal result

Q = 1

99

1.0042

rP = −0.134, p = 0.187

No statistically significant marginal velocity dependence

Q = 1, 30° ≤ i ≤ 80°

75

1.0078

rP = −0.220, p = 0.058

Relation remains borderline/non-significant after inclination control

Q = 1, 30° ≤ i ≤ 80°, valid Vflat

65

0.9948

rP (Δlog v,g, Vflat) = −0.084, p = 0.440

Complete-catalog trend is absent using the published flat rotation velocity

Distance control

175

—

rP = −0.03, p = 0.69

No significant catalog-wide residual dependence on distance

Inclination control

175

—

rP = −0.039, p = 0.608

No significant catalog-wide residual dependence on inclination

Galaxy-property controls, Q = 1

99

—

All direct Pearson tests p > 0.34

No significant direct linear dependence on luminosity, surface brightness, morphology, gas richness, or radial coverage

Global stellar-normalization sensitivity, Q = 1

99

0.9630–1.0507 for disk M/L = 0.42–0.60

rP = −0.212 to −0.049; p = 0.035–0.631

Fiducial 1.0042 depends on stellar source normalization; no galaxy-specific M/L is fitted

WISE stellar source substitution

111 overlap; 71 with Q = 1; 10 matched pairs

Q = 1: 1.0697 versus same-galaxy fiducial 0.9863

All 111: rP = −0.1182, p = 0.2166; Q = 1: rP = +0.2140, p = 0.0732; matched WISE-source pairs: 10/10 negative, p = 0.001953

External stellar masses remove or reverse the negative marginal relation but do not remove the residual-blind conditional ordering

Galaxy-cluster bootstrap

175

Pooled mean ratio 1.00019; 95% interval 0.9742–1.0286

Q = 1 ratio 1.0042; 95% interval 0.9763–1.0325

Galaxy-level resampling preserves radial clustering; zero pooled bias and Q = 1 unity normalization remain statistically compatible with their null values

Independent LITTLE THINGS population

15

0.9907; 95% galaxy-block interval 0.868–1.130

Primary log-velocity rP = −0.282, p = 0.308; rS = −0.211, p = 0.451; partial rP = −0.257, p = 0.355 controlling SED stellar mass and inclination

Independent population supports near-unity fixed normalization but does not reproduce the SPARC conditional velocity ordering; strict SPARC-style matching yields only two pairs, both opposite in sign

The SPARC controls impose two distinct qualifications. First, the complete-catalog marginal velocity dependence is not recovered as a statistically significant relation in the independently classified Q=1 population, in the principal inclination-controlled subset, or when the published Vflat is used as the characteristic velocity. Second, stellar source substitutions demonstrate that the unusually close fiducial normalization is not source invariant. The independent WISE source moves the 71-galaxy Q=1 overlap from a same-galaxy fiducial ratio of 0.9863 to 1.0697, while the negative marginal relation becomes non-significant across all 111 WISE-overlap galaxies and reverses sign in the Q=1 overlap. Conversely, the residual-blind WISE-mass matched pairs and simultaneous WISE-source regression retain conditional velocity ordering, although velocity and WISE stellar mass remain strongly collinear. The SPARC ordering is therefore a genuine within-catalog empirical structure whose physical origin remains unresolved, rather than either a demonstrated universal defect in the MQ galactic realization law or a residual established to be observational in origin.

Poorer-quality populations exhibit progressively larger normalization offsets, but SPARC quality is associated both with properties of the kinematic data and, in the present catalog, with galaxy velocity and astrophysical regime. The absence of a significant marginal velocity relation in Q=1 shows that the complete-catalog trend is not uniformly expressed across observational quality; it does not show that observational quality explains the trend. The independent LITTLE THINGS population supplies the universality test: the SPARC ordering is not statistically significant there under logarithmic-velocity, rank, stellar mass, inclination, regression, source-sensitivity, or matched-pair diagnostics.

The population-wide results of Sec. 4.9 add a separate distinction. Strong radial co-variation between MQ realized and observed rotation curves can coexist with substantial galaxy-level normalization displacement. Pearson correlation measures similarity of radial co-variation rather than equality of velocity scale [61], while RMS, MAE, and the bias-to-scatter ratio Bg characterize complementary aspects of predictive behavior. High curve-shape correlation therefore cannot by itself establish absolute velocity agreement.

The Milky Way and SPARC results consequently constrain different aspects of the same fixed law. The RAR/MOND-independent Milky Way calculation demonstrates a large reduction of the baryonic Newtonian velocity discrepancy under a separately constructed source and kinematic target, with statistical compatibility once the independently published McMillan baryonic-normalization covariance is propagated. SPARC supplies the much larger heterogeneous population test, but its interpretation is more source dependent. The gas-rich controls of Sec. 4.6 localize much of this sensitivity: for the 19 Q=1, inclination-controlled galaxies with

MHI/L3.6≥0.5
(S143)

, the WISE substitution changes the geometric normalization from 0.9854 to 0.9886; for the eight systems with

MHI/L3.6≥1.0
(S144)

, it changes from 1.0086 to 1.0034. Stellar normalization is therefore a major source-side uncertainty where stars are dynamically important rather than a universal explanation of the MQ normalization.

These results supply a direct falsifiability criterion for the fixed MQ galactic realization law. A reproducible residual dependence on rotational velocity or another physical galaxy property that persists in independently high-quality populations after baryonic source and observational controls would constitute evidence that the fixed law is incomplete. The WISE stellar mass matched analysis supplies such pressure within SPARC, but because its stellar masses are integrated population model estimates and the kinematic target remains SPARC, an external population is required to distinguish a catalog-specific systematic from a universal scale dependence.

That external population test is carried out with LITTLE THINGS rather than by comparison with a fitted RAR acceleration parameter. Oh et al. [5] present homogeneous high resolution rotation curves and mass models for 26 nearby dwarf galaxies selected from the LITTLE THINGS VLA H I survey [9], while O’Brien et al. [6] publish corresponding measured rotation curve points and Newtonian-component curves for 25 of those systems. The primary sample was fixed before examining an MQ residual. DDO154, DDO168, and NGC2366 were excluded because they overlap the SPARC population. DDO43, DDO46, DDO47, F564-V3, and Haro29 were excluded because the Oh et al. source construction lacks the Spitzer stellar contribution required for the full source test; DDO52 is absent from the 25-galaxy O’Brien figure set; and IC10 and Haro36 lack the independent SED stellar mass value used for the primary source normalization. These availability and provenance exclusions leave 15 galaxies in the primary out-of-SPARC sample.

The LITTLE THINGS calculation is performed in physical velocity space. Observed circular velocity markers and their plotted vertical uncertainties were extracted independently from the paired conformal-gravity and MOND panels of O’Brien et al. [6]. The duplicate rendering of each observed curve provides an extraction check: the median absolute difference between the independently digitized paired-panel velocities is 0.15 km s-1. The plotted Newtonian component is not used directly as the decisive baryonic source because O’Brien et al. fit stellar disk mass in their alternative-gravity calculations. Instead, the corresponding fitted Freeman-disk contribution is evaluated from the published disk mass and scale length and removed in quadrature from the plotted Newtonian component to recover the represented gas contribution. The stellar term is then reconstructed using the independent SED stellar masses tabulated by Oh et al. [5]. The observed velocities enter only as the external comparison target; neither aMQ, an interpolation parameter, nor a galaxy-specific MQ quantity is fitted to them.

The resulting primary calculation contains 288 radial comparisons across the 15 galaxies. With the present-epoch aMQ of Eq. (133) held fixed, Eqs. (136) and (135) give a galaxy-weighted mean logarithmic MQ-to-observed normalization of -0.00407 dex, corresponding to a geometric velocity ratio of 0.9907. A galaxy-block bootstrap gives a 95% interval of 0.868–1.130, and the one-sample test of the 15 galaxy mean logarithmic ratios gives p=0.895 for zero population normalization offset. The point-weighted signed residual is -0.39 km s-1, while pooled RMS and MAE are 12.55 and 8.69 km s-1, respectively. The population normalization is therefore statistically compatible with unity, but the individual-galaxy and radial scatter precludes interpretation as a precision pointwise acceptance test. The independently selected population supports the fixed MQ velocity normalization at the population level without refitting aMQ or any galaxy-specific MQ parameter.

LITTLE THINGS also tests whether the conditional velocity ordering identified within SPARC generalizes outside that catalog. Using logarithmic characteristic observed velocity for direct comparability, the Pearson correlation with galaxy-mean logarithmic normalization is rP=-0.282 with p=0.308, and the Spearman coefficient is rS=-0.211 with p=0.451. The linear-velocity Pearson relation is stronger in magnitude but remains non-significant (rP=-0.450, p=0.093). Partial control for SED stellar mass gives rP=-0.297 (p=0.282); controlling both stellar mass and inclination gives rP=-0.257 (p=0.355), with a 100,000-permutation probability of 0.349. Adding distance reduces the partial coefficient to -0.133 (p=0.638). A simultaneous HC3 regression of mean logarithmic normalization on standardized logarithmic velocity, logarithmic stellar mass, and inclination gives a velocity coefficient of -0.0435 dex per standard deviation with robust SE =0.0935 and p=0.642. Velocity and SED stellar mass remain correlated (rP=0.830), but their VIFs are approximately 3.3, substantially below those of the SPARC WISE-source regression. The conclusion is also unchanged using the alternative Oh et al. kinematic stellar mass source: rP=-0.366 (p=0.180) and rS=-0.136 (p=0.630).

Applying the SPARC residual-blind matched-pair criteria to the smaller LITTLE THINGS population is underpowered but gives no evidence for the same conditional ordering. At the primary 0.10-dex stellar mass tolerance, only two non-overlapping pairs are available, and both higher-velocity members have higher, rather than lower, MQ-to-observed normalization; their mean difference is +0.101 dex with Wilcoxon p=0.500. Relaxing the stellar mass tolerance to 0.20 dex gives three pairs with mean difference +0.042 dex (p=0.750), while 0.30 dex gives four mixed pairs with mean difference -0.0117 dex (p=0.875). These small samples cannot establish velocity independence, but they do not reproduce the stable negative ordering found in the SPARC WISE-matched controls. The SPARC ordering is therefore retained as a catalog-specific empirical systematic requiring explanation rather than promoted to a universal velocity-dependent correction to Eq. (137).

The independent LITTLE THINGS population also probes the stellar source sensitivity identified within SPARC. The primary calculation uses the Oh et al. SED stellar masses; a source sensitivity calculation replaces them with the independent Oh et al. kinematic stellar mass estimates while leaving the fixed MQ law unchanged. Across all 15 galaxies, the galaxy-weighted geometric normalization changes from 0.9907 to 0.9655. The signed source-induced change in galaxy-mean logarithmic normalization correlates with stellar importance, measured by

log10(M⋆/Mgas)
(S145)

: rP=-0.541 with p=0.037, while the rank correlation is rS=-0.489 with p=0.064. For the ten gas-dominated galaxies with

M⋆/Mgas<1
(S146)

, the mean signed source shift is -0.00117 dex, corresponding to a geometric change of approximately -0.27%, with a mean absolute shift of 0.0122 dex. For the five stellar-important galaxies with

M⋆/Mgas≥1
(S147)

, the mean signed shift is -0.0312 dex, corresponding to approximately -6.9%, with a mean absolute shift of 0.0338 dex. The independent population therefore reproduces the central SPARC source-side conclusion: stellar source choice has little effect where gas dominates but can materially alter predicted normalization where stars are dynamically important. This is evidence for baryonic source degeneracy rather than an additional MQ degree of freedom.

That result localizes the normalization sensitivity but does not resolve the remaining radial-shape uncertainty. Changing an integrated stellar mass rescales a common radial stellar template; it does not replace that template with an independently measured radial stellar surface-density profile. A residual-blind route is available for the five stellar-important galaxies already selected by the source-side criterion

M⋆/Mgas≥1
(S148)

: DDO50, DDO216, IC1613, NGC1569, and NGC3738. Zhang et al. [57] provide spatially resolved stellar mass surface-density distributions from multiband SED fits for these systems, while Oh et al. [5] provide the corresponding homogeneous LITTLE THINGS H I kinematics and resolved Spitzer 3.6 μm mass-model information. Those products define the next radial-source test without selecting galaxies from their MQ residuals.

For that experiment, the published radial stellar profile and H I source are to be adopted with their published disk geometry, converted to Newtonian stellar and gas contributions, and combined into the baryonic acceleration entering Eq. (137). The present-epoch acceleration scale and the galactic realization law remain fixed; neither a radial mass-to-light adjustment nor a galaxy-level renormalization may be selected after the observed rotation curve residual is inspected. The primary comparison is the change in each galaxy’s mean logarithmic MQ-to-observed normalization and radial residual structure relative to the existing SED-normalized reconstruction. A persistent coherent residual after the independently resolved radial-source replacement would place pressure on the fixed realization law, whereas a material reduction would identify the present stellar-dominated discrepancy with baryonic source structure. Until those radial profiles are ingested numerically, the manuscript does not claim precision validation for the stellar-important LITTLE THINGS systems.

5. DISCUSSION, INDEPENDENT EXTENSIONS, AND OPEN TESTS

5.1 Geometric origin of the dark matter phenomenon
The galactic tests of Secs. 4.3–4.10 show that the fixed MQ galactic realization law substantially reduces the discrepancy between baryonic Newtonian predictions and observed orbital velocities in both the Milky Way and the complete SPARC population, without introducing a non-baryonic gravitating component or galaxy-specific fitted MQ parameters. Their significance should be interpreted in light of the pre-existing RAR: the contribution is not rediscovery of its empirical functional shape, but the application of an independently specified acceleration scale and physical realization construction without galaxy-specific MQ calibration [8]. This statement is deliberately narrower than “parameter-free.” Stellar mass-to-light normalization, distance, inclination, gas distribution, and other quantities entering the published baryonic reconstructions remain external astrophysical inputs. The fixed global SPARC mass-to-light sensitivity in Sec. 4.6 directly demonstrates that the fiducial normalization is not invariant under plausible baryonic source choices.

In MQ, the Newtonian baryonic source relation remains unchanged. The observed baryonic mass distribution determines the gravitational availability represented by the Newtonian baryonic acceleration gbar, while the corresponding baryonic velocity vNew is calculated from the reconstructed bulge, disk, and gas contributions. No additional mass component is inserted into this source description, and the inverse-square dependence of the baryonic gravitational field is retained.

The departure from the Newtonian baryonic prediction arises through bounded realization of that gravitational availability. MQ distinguishes the gravitational availability supplied by the baryonic source from the realized gravitational response in the System Frame. Through the Frames mapping, the realization fraction is count constrained and approaches unity according to the galactic realization law developed in Appx. F. At high baryonic acceleration, realization approaches completion and the Newtonian baryonic response is recovered. At lower acceleration, the realization fraction decreases in the sense defined by the MQ construction, and the realized gravitational response correspondingly exceeds the Newtonian baryonic value.

The dark matter phenomenon therefore arises, within MQ, from the distinction between gravitational availability and realized gravitational response rather than from additional gravitating matter. An observer interpreting the realized orbital velocity strictly within Newtonian dynamics attributes the excess response to additional enclosed mass because Newtonian gravity contains no separate distinction between the available source relation and its realized response. MQ instead assigns that discrepancy to mapping-induced deformation realized geometrically in the System Frame. In this sense, the inferred “missing mass” is the Newtonian representation of a difference that MQ assigns to realized geometry.

Appendix F derives the galactic realization law from the MQ realization construction. The acceleration scale aMQ is supplied independently by the preceding MQ cosmological construction, and the realized velocity then follows from the baryonic source and previously established MQ quantities without galaxy-dependent adjustment. The observational tests in Sec. 4 therefore ask whether this fixed realization mechanism reproduces galactic dynamics when applied unchanged to distinct published baryonic constructions.

The novelty hierarchy is consequently specific. The empirical RAR is prior community knowledge and is not claimed as a discovery of this paper [8]. The present galactic claim is that MQ independently specifies the characteristic acceleration scale and applies a fixed galactic realization law without galaxy-specific MQ calibration, while proposing a physical origin for that relation in the Frames mapping and count-constrained gravitational realization. The Milky Way provides a particularly useful historical test because the conventional McMillan baryonic model predates the canonical 2016 RAR publication, while the Eilers kinematics postdate it and were not part of that discovery sample. With the McMillan baryonic-normalization covariance propagated as described in Sec. 4.3, the fixed MQ prediction is statistically compatible with those kinematics, although the test is not covariance-complete.

The SPARC and LITTLE THINGS analyses qualify this interpretation. The complete SPARC catalog contains a velocity-dependent normalization structure, but its marginal form is sensitive to observational quality and stellar source choice, and the conditional SPARC ordering does not replicate significantly in the independent LITTLE THINGS population. Absolute normalization is also most sensitive to stellar source choice where stars are dynamically important, whereas gas-dominated systems are substantially more source stable in both populations. Baryonic source structure therefore remains a material uncertainty in stellar-important galaxies, preventing precision validation of the galactic realization law in every system. It does not introduce an additional MQ degree of freedom: the acceleration scale and galactic realization law remain fixed throughout the tests.

5.2 Emergence of flat rotation curves and baryonic Tully-Fisher scaling

The physical interpretation developed in Sec. 5.1 has a direct dynamical consequence. Because the MQ galactic realization law approaches the Newtonian baryonic response at high acceleration but departs from it systematically in the weak-acceleration regime, the same fixed relation produces Newtonian behavior at high acceleration and asymptotically flat orbital velocities at large galactocentric radius. Its weak-acceleration limit also yields the baryonic Tully-Fisher fourth-power scaling. These results follow directly from the exact galactic realization law established in Sec. 4.3.

At high baryonic acceleration, gbar≫aMQ, the realization fraction approaches unity. The denominator of the exact galactic realization expression therefore approaches one, and the realized velocity vreal converges continuously to the Newtonian baryonic velocity vNew. Thus, sufficiently high-acceleration regions recover Newtonian orbital dynamics. The same limiting behavior is observed in empirical radial-acceleration studies [8].

The outer-galaxy behavior follows from the opposite limit. For a finite, localized baryonic mass distribution, the Newtonian baryonic acceleration enters the weak-acceleration regime as galactocentric radius R increases. The baryonic source retains its Newtonian far-field behavior, while the realized gravitational response in the System Frame follows the galactic realization law developed in Appx. F.

Define the realized centripetal acceleration greal(R) at galactocentric radius R by

greal(R)≡vreal2(R)R.
(152)

Squaring the realized-velocity relation of Eq. (137) gives

vreal2(R)=vNew2(R)1-exp⁡[-gbar(R)/aMQ].
(153)

Dividing by R and using Eq. (131) yields

vreal2(R)R=vNew2(R)/R1-exp⁡[-gbar(R)/aMQ].
(154)

Hence the exact acceleration relation is

greal(R)=gbar(R)1-exp⁡[-gbar(R)/aMQ].
(155)

The disappearance of vNew is only a change of variables; its contribution remains contained in

gbar=vNew2/R
(S149)

.

To evaluate the weak-acceleration limit, define the dimensionless variable x by

x≡gbaraMQ.
(156)

The galactic realization law becomes

greal=aMQx21-e-x.
(157)

The galactic realization law becomes

1-e-x=x-x22!+x33!-x44!+⋯.
(158)

Consequently,

1-e-xx→1(x→0),
(159)

and therefore,

grealaMQx→1(x→0).
(160)

From Eq. (156), the exact identity

aMQx=aMQgbar,
(161)

converts Eq. (160) to the dimensionless weak-acceleration limit

grealaMQgbar→1(gbaraMQ→0).
(162)

Equation (162) establishes the square-root weak-acceleration behavior as an asymptotic consequence of the exact galactic realization law rather than as an independently inserted acceleration relation. The same deep-acceleration scaling was introduced in the original MOND formulation [19], providing an established phenomenological benchmark for the limit obtained here.

For a finite localized baryonic distribution with asymptotic baryonic mass Mbar, the Newtonian far field satisfies

limR→∞R2gbar(R)=GMbar,
(163)

where G is the gravitational constant resolved at the upper count limit, as defined earlier. Combining Eqs. (162) and (161) gives

limR→∞R​greal(R)=GMbaraMQ.
(164)

Thus, the realized gravitational response decreases asymptotically as

1/R
(S150)

, rather than as the Newtonian baryonic

1/R2
(S151)

far field. Circular motion satisfies

greal=vreal2R,
(165)

so Eq. (164) implies

limR→∞vreal2(R)=GMbaraMQ.
(166)

Equivalently,

limR→∞vreal4(R)GMbaraMQ=1.
(167)
limR→∞vreal4(R)=GMbaraMQ.
(168)

The finite-mass outer limit therefore gives both asymptotically constant realized orbital velocity and the baryonic Tully-Fisher fourth-power scaling v4∝Mbar. The latter is independently well established observationally [20] and is also a standard consequence of the deep-MOND limit.

The exact functional form underlying these limits is mathematically the same as the empirical radial acceleration relation discussed in Sec. 4.3 [8], but its stated provenance in MQ differs. The acceleration scale is fixed by the preceding MQ cosmological construction, while Appx. F derives the realization fraction and its exponential form independent of the galaxy sample. Flat outer rotation and baryonic Tully-Fisher scaling are therefore downstream asymptotic consequences of the fixed galactic realization law rather than conditions used to select its weak-acceleration behavior. The observational analysis of Sec. 4 tests this construction without galaxy-specific MQ adjustment. Its residual structure and baryonic source sensitivities remain empirical falsification targets.

5.3 Gravitational lensing
Gravitational lensing provides a relativistic test independent of the rotation curve observable. The MQ construction supplies a spherical weak-field consistency calculation and, once the relevant System Frame metric has been solved, the corresponding thin-lens projection. The ESO 325-G004 spherical aperture calculation is retained as a numerical benchmark. A complete extended-source prediction requires the covariant metric solution, including both weak-field metric potentials and their gravitational slip, before deflection, convergence, and shear maps are evaluated.

Let

α̂
(S152)

denote the physical weak-field gravitational deflection angle and b the lensing impact parameter, consistent with the notation introduced in Sec. 1.3 and used in Sec. 2.3. For a point lens, or for a light ray exterior to a spherically symmetric source in the standard weak-field limit,

α̂=4GMc2b,
(169)

where M is the gravitating mass sourcing the exterior field.

Within MQ, let Mreal(R) denote the enclosed gravitational source associated with the realized gravitational response at radius R, and let

α̂MQ
(S153)

denote the corresponding MQ spherical weak-field deflection angle. Applying that source to the same exterior spherical relation gives

α̂MQ=4GMreal(R)c2b.
(170)

Equation (170) is an enclosed-source representation of the exterior spherical consistency calculation, not an independent lensing mass or a substitute for the relativistic metric solution.

The broader MQ construction supplies the relativistic bridge. The emergent System Frame metric is generated through the same Frames mapping, the metric and matter source enter the single MQ variational principle, and freely propagating trajectories obey the resulting geodesic equation [1]. For the static, spherically symmetric weak-field consistency class, take the realized radial acceleration of Eq. (136) to be the radial derivative of the weak-field potential and impose the negligible anisotropic stress limit in which the temporal and spatial metric perturbations enter photon deflection with the standard general relativistic coefficient. Let z denote distance along the unperturbed line of sight, with

r=(b2+z2)1/2
(S154)

. The null-geodesic weak-field reduction then gives [62]

α̂MQ(b)=2bc2∫-∞+∞greal​(b2+z2)b2+z2​dz.
(171)

Equation (171) is the standard spherical weak-field deflection integral expressed directly in terms of the radial gravitational acceleration realized by MQ. No lensing-specific response function is introduced. Substitution of the exact galactic realization law of Eq. (136) yields

α̂MQ(b)=2bc2∫-∞+∞gbar(r)r[1-exp⁡​(-gbar(r)/aMQ)]​dz,​r=b2+z2.
(172)

In the finite-mass outer limit established in Sec. 5.2,

gbar(r')→GMbar/r'2
(S155)

, and the galactic realization law gives

greal(r')→GMbaraMQ/r'
(S156)

. With

r'=b2+z2
(S157)

, Eq. (171) therefore reduces asymptotically to an integral proportional to

b/(b2+z2)
(S158)

. Using

∫-∞+∞dz/(b2+z2)=π/b
(S159)

, the impact-parameter dependence cancels and the deflection approaches the constant limit

limb→∞α̂MQ(b)=2πc2GMbaraMQ.
(173)

Define v∞ as the asymptotic realized orbital velocity determined by the finite-mass outer limit of Sec. 5.2,

v∞2≡limR→∞vreal2(R)
(S160)

. Equation (173) may then be written

limb→∞α̂MQ(b)=2πv∞2c2.
(174)

Equations (171-174) provide a cross-observable spherical consistency condition: within the stated weak-field assumptions, the galactic realization law that fixes the asymptotic orbital velocity also fixes the asymptotic photon deflection without an additional acceleration scale, halo profile, or lensing-specific MQ parameter. Its relativistic interpretation remains conditional on recovering the assumed relation between the temporal and spatial metric perturbations from the covariant System Frame solution.

For an extended source, Eqs. (175-179) are used only in the negligible-slip specialization of that solved weak-field metric. Let Φsys denote the common gravitational potential realized in the System Frame stationary weak-field limit under this specialization. The standard thin-lens projection maps that potential to observable lensing fields [62]. Let θ denote the two-dimensional angular position vector in the lens plane, and let Dl, Ds, and Dls denote the angular diameter distances from observer to lens, observer to source, and lens to source, respectively. Define the dimensionless lensing potential ψ by

ψ(θ)=2Dlsc2DlDs∫-∞+∞Φsys(Dlθ,z)​dz.
(175)

The angular deflection field α follows from the gradient of the same projected potential,

α(θ)=∇θψ(θ).
(176)

while κ denotes convergence and γ1 and γ2 the two scalar shear components,

κ(θ)=12[∂2ψ∂θ12+∂2ψ∂θ22].
(177)
γ1(θ)=12[∂2ψ∂θ12-∂2ψ∂θ22].
(178)
γ2(θ)=∂2ψ∂θ1∂θ2.
(179)

Equations (175–179) are standard null-geodesic/thin-lens consequences of the stated weak-field specialization and do not constitute an additional MQ galactic realization law [62]. Section 5.8 develops the action-level closure test needed to determine whether the stationary weak-field construction extends consistently beyond the spherical limit. For the same baryonic source, the covariant System Frame equations must determine both weak-field metric potentials and their gravitational slip before the appropriate lensing potential is selected. A disagreement between the resulting parameter-frozen α, κ, or γ maps and observation would directly falsify the extended-lens embedding.

A practical benchmark is the nearby strong lens ESO 325-G004. Smith and Lucey measured a lens redshift of 0.0347, a source redshift of 2.141, an Einstein radius of 2.85 arcsec, and an F814W luminosity of

(4.07±0.08)×1010
(S161)

solar luminosities within that radius [63]. Their VLT stellar population analysis gives a Milky Way-like stellar mass-to-light ratio of 3.01±0.25 in the same band. The conventional strong-lensing reconstruction gives

(1.50±0.06)×1011
(S162)

solar masses inside the Einstein aperture [63]. Because the Einstein radius is only about one quarter of the stellar effective radius, the aperture is strongly stellar dominated. Collett et al. subsequently combined the HST strong-lensing structure with spatially resolved MUSE stellar kinematics and obtained a weak-field spatial-curvature parameter per unit mass of 0.97±0.09 relative to the general relativistic value [64]. These measurements provide both an independently constrained baryonic source and a high-quality relativistic lensing target.

For an MQ test, the source construction must be frozen before the lensed arcs are used. The HST stellar-light distribution and independently determined stellar population normalization should define the baryonic source; lensing-inferred total mass, dark matter fraction, or a lensing-tuned stellar normalization must not be fed back into it. That source is then propagated through the stationary System Frame equations and the resulting metric potentials. Only after the metric and lensing maps are frozen should the predicted Einstein-ring morphology and α, κ, γ1, and γ2 fields be compared with the HST constraints. The acceleration scale of Eq. (133) and the galactic realization law remain unchanged. No lens-specific acceleration scale, halo profile, or post hoc shear correction is permitted apart from an independently measured environmental contribution.

A preliminary spherical benchmark can nevertheless be evaluated without using the lensing-inferred mass as an MQ source. Multiplying the measured F814W luminosity by the independently inferred stellar population mass-to-light ratio gives a stellar mass of

(1.225±0.105)×1011
(S163)

solar masses within the observed Einstein aperture. At 0.687 kpc arcsec-1, the 2.85-arcsec Einstein radius corresponds to 1.958 kpc [63]. Numerical evaluation of Eq. (172), with aMQ fixed by Eq. (133), enhances the spherical deflection from this stellar source by a factor 1.0475 relative to the general relativistic point-mass value at the same impact parameter. The corresponding deflection is equivalent to

1.283×1011
(S164)

solar masses under the standard point-mass normalization, below the observed

(1.50±0.06)×1011
(S165)

solar-mass lensing mass. Solving the spherical lens equation with the published geometry gives an MQ Einstein radius of approximately 2.67 arcsec, compared with the observed 2.85 arcsec. Propagating the quoted luminosity, stellar population mass-to-light, and lensing-mass uncertainties gives an equivalent-mass shortfall of about 1.8 standard deviations.

This result is a benchmark rather than the final ESO 325-G004 test. The spherical reduction does not use the resolved HST stellar-light distribution, determine the relativistic gravitational slip of the extended solution, or predict the two-dimensional ring morphology and shear field. The decisive test is the parameter-frozen extended calculation: the observed HST source distribution and independently fixed stellar normalization must be propagated through the covariant System Frame metric, after which the resulting α, κ, γ1, and γ2 maps are compared directly with the observed Einstein-ring structure. Because the spherical benchmark lies below the observed normalization, recovery of the ring must arise from the covariant extended-source geometry itself, not from a lens-specific acceleration scale, dark matter halo, or lens-tuned stellar normalization [62–66].

5.4 Cluster collisions

Cluster mergers test whether the same action-level MQ construction can generate the observed nonequilibrium separation of gas and lensing curvature. The equations below retain the initial-value and causal structure needed to pose that calculation. They do not constitute a Bullet Cluster solution: nonlinear evolution from independently specified baryonic merger data, followed by ray tracing through the evolved metric, remains required.

The formal time dependence can be stated without introducing a new MQ timescale. Let

gμνsys
(S166)

denote the emergent System Frame metric and

Tμνsys
(S167)

the stress-energy tensor defined by metric variation of the realized matter action within the MQ variational principle. Observed gas and stellar/galaxy phase-space distributions constrain the material fields entering that action;

Tμνsys
(S168)

is not introduced independently as a conventional baryonic source. The macroscopic field equation and conservation condition are

Gμνsys[gsys(x,t)]+Λsysgμνsys(x,t)=(8πGc4)Tμνsys(x,t).
(180)
∇μTμνsys(x,t)=0.
(181)

Writing

h̅μνsys
(S169)

for the trace-reversed weak-field metric perturbation and □≡ηαβ∂α∂βfor the corresponding flat-background d’Alembert operator [67], the weak-field harmonic-gauge expansion about Minkowski spacetime has the first-order hyperbolic form

□̄​h̄μνsys(x,t)=-(16πGc4)Tμνsys(x,t),
(182)
h̄μνsys(x,t)=(4Gc4)∫Tμνsys​(x',t-|x-x'|/c)|x-x'|​d3x'.
(183)

Equations (180–183) show that a time-varying realized source produces a causal gravitational response without a free relaxation time or phenomenological memory parameter. For a strongly nonlinear merger, Eq. (180), not its linearized reduction, is the required evolution law.

Section 5.8 develops the action-level closure test whose stationary spherical limit must recover the galactic realization law of Eq. (136). The unresolved cluster problem is therefore dynamical: the observed baryonic merger configuration has not yet been evolved with the nonlinear MQ equations, and no resulting MQ convergence or shear map for 1E 0657-558 has been compared with the data.

For 1E 0657-558 at redshift z=0.296, strong- and weak-lensing reconstructions place the two principal projected mass centroids about 720 kpc apart at the observed stage [68]. The lensing structure is displaced from the dominant X-ray-emitting gas, providing the essential spatial constraint that any successful evolution must reproduce. Hydrodynamic merger calculations further show that the bow-shock velocity, roughly 4500-4700 km s-1, is not the subcluster mass-centroid velocity: bullet-like configurations can reproduce that shock speed with a mass-centroid speed near 2600-2700 km s-1 because the upstream gas is already infalling and the shock outruns the subcluster [69–70]. These are distinct observables and must not be used interchangeably as initial conditions.

The causal timescales likewise exclude an adjustable long-lived gravitational response as an explanation of the spatial separation. Across 720 kpc, the light-crossing time is approximately 2.35 Myr, whereas a subcluster moving at 2600-2700 km s-1 requires approximately 261-271 Myr to traverse that distance. Equation (183) therefore cannot by itself sustain a merger-scale gravitational lag. If MQ reproduces the observed gas-lensing separation, it must arise from the evolved spatial and dynamical structure of the action-derived realized stress-energy source and nonlinear System Frame metric, not from an adjustable propagation delay.

The stationary spherical consistency condition is

|∇Φsys(R)|=greal(R)=gbar(R)1-exp⁡​[-gbar(R)/aMQ].
(184)

Equation (184) is not a cluster-specific law. It requires the stationary weak-field source constructed from the action to reproduce Eq. (136) in spherical symmetry. In the stationary weak-field branch, the nonrelativistic source density is represented by

T00sys/c2
(S170)

, so

∇2Φsys=4πGc2​T00sys.
(185)

The enclosed realized gravitational source implied by the realization fraction P(R) is

Msys(<R)=Mbar(<R)P(R).,
(186)
T00sys(R)c2=14πR2ddR[Mbar(<R)P(R)].
(187)

Equations (185–187) identify the stationary realized gravitational source required to reproduce the spherical galactic relation. The action-level construction of Sec. 5.8 sharpens this statement:

Tμνsys
(S171)

is not an independently assigned halo stress-energy tensor but the metric response of the realized matter action [71].

Let Φbar denote the Newtonian potential of the adopted baryonic source, satisfying ∇2Φbar=4πGρbar. In stationary spherical symmetry, Eq. (136) fixes the realized field parallel to the baryonic field with magnitude increased by P-1. Thus

∇Φsys=P-1∇Φbar,
(188)
T00sysc2=ρbarP+14πG∇(P-1)⋅∇Φbar.
(189)

This stationary closure does not constitute a Bullet Cluster prediction. The quantitative test requires independently specified baryonic merger data to be evolved through the full nonlinear Eq. (180), followed by calculation of lensing observables from the evolved metric. Equation (183) is only a causal weak-field check; no centroid offset, convergence map, shear map, or shock-compatible merger state is claimed here as an MQ result.

The action-level requirement is consistent with established Modified gravity field theory. Bekenstein and Milgrom [72] embed the MOND response in a Lagrangian-derived nonlinear Poisson equation with associated conservation laws, while QUMOND [73] obtains an effective source through an action-derived field construction rather than assigning a phantom density independently. These theories are not imported into MQ; they illustrate why a nonlinear response intended to describe an evolving gravitating system must be embedded in a field equation with a consistent variational and conservation structure.

5.5 Cosmic-age dependence of the galactic realization law
The predicted dependence of the galactic realization law on the System Frame age of the universe provides a direct discriminator between MQ and a purely present-epoch empirical RAR. Once the age relation is fixed, MQ specifies aMQ(z) without galaxy-specific fitting. A direct test is therefore a preregistered resolved-galaxy analysis over approximately 0.5<z<1.5, using independently reconstructed stellar and gas sources, a prespecified expansion history, and frozen MQ parameters [74–75].

Once θsi is fixed independently, the age dependence of aMQ is determined by Hint. Within the cosmological expansion construction developed in Sec. 3.12,

Hint∝AU-1
(S172)

, where AU denotes the System Frame age of the universe at the epoch under consideration and AU,0 its present-epoch value. Within MQ, the Fixsen CMB temperature measurement [41] may be inserted into the previously derived age relation [15, Appx. BI] to obtain the present normalization,

AU=(3θsi(2θsi)1/3e(3θsi3)/24πaclfTCMB4)1/3=13.7943(39)billionyears,
(190)
aMQ(AU)aMQ(AU,0)=Hint(AU)Hint(AU,0)=AU,0AU.
(191)

Equation (190) supplies the MQ age normalization used here; because it is obtained from the measured CMB temperature through the MQ construction, it is not treated as an independent observational determination of the age of the universe.

For comparison with observations at finite redshift, the age dependence can be expressed without identifying an observationally inferred expansion rate with the Internal Frame rate. Let z denote cosmological redshift, let a(t) denote the System Frame cosmological scale factor at cosmic time t, and let Hsys(z) retain its established meaning as the System Frame expansion rate from Sec. 3.12. From

1+z=a(0)/a(t)
(S173)

and

Hsys=ȧ/a
(S174)

, the System Frame age of the universe at redshift z is

AU(z)=AU,0-∫0zdz'(1+z')Hsys(z').
(192)
aMQ(z)=aMQ(0)AU,0AU(z).
(193)

Equation (192) is a coordinate conversion conditional on the externally specified System Frame expansion history Hsys(z); it does not fit that history from the galaxy data. A prospective test must therefore specify Hsys(z) before evaluating the galactic sample, keeping uncertainty or model dependence in the adopted expansion history distinct from the MQ galactic realization law test.

The acceleration scale is then applied point by point to the independently reconstructed gravitational availability of the baryonic source through the same fixed galactic realization law,

vreal(R,AU)=vNew(R)1-exp⁡[-gbar(R)/aMQ(AU)].
(194)

Here R is the realized galactocentric radius, vNew(R) the Newtonian baryonic orbital velocity, gbar(R) the corresponding Newtonian baryonic acceleration, and vreal(R,AU) the realized orbital velocity at System Frame age AU. At fixed baryonic acceleration, decreasing AU increases aMQ and therefore the realized gravitational response relative to the Newtonian baryonic prediction. The magnitude remains acceleration dependent through

gbar(R)/aMQ(AU)
(S175)

: high-acceleration regions remain comparatively close to the Newtonian limit, whereas lower-acceleration regions exhibit the larger fractional change.

Holding the baryonic source fixed isolates the age-dependent prediction,

vreal(R,z)vreal(R,0)=1-exp⁡​[-gbar(R)/aMQ(0)]1-exp⁡​[-gbar(R)/aMQ(z)].
(195)

Here aMQ(0) is the present-epoch value of the same MQ acceleration scale. Equation (195) predicts a positive velocity shift toward higher redshift at fixed gbar, with the larger fractional effect in the low-acceleration regime. This is more discriminating than allowing an independently fitted phenomenological acceleration parameter in each redshift bin.

For scale, Table 21 evaluates Eqs. (193) and (193) using the present System Frame age AU,0=13.7943 Gyr and the Planck 2018 flat late-time expansion benchmark already adopted elsewhere in the manuscript [11]. The table is not a fit to galaxy data. It converts the prescribed expansion history to the age variable entering Eq. (191). The two velocity columns hold the physical baryonic acceleration fixed at gbar=0.1aMQ(0) and gbar=aMQ(0), respectively, and report the change relative to z=0.

Table 21. Redshift-dependent MQ acceleration scale and fixed-baryonic source velocity shift.

z
(S176)

Cosmic age (Gyr)

Scale ratio

aMQ(z)/aMQ(0)
(S177)

Acceleration scale (

10-10ms-2
(S178)

)

Velocity shift gbar=0.1aMQ(0)

Velocity shift gbar=aMQ(0)

0.08

12.699

1.0862

1.3039

+1.78%

+1.23%

0.20

11.274

1.2235

1.4687

+4.42%

+3.07%

0.50

8.585

1.6068

1.9288

+10.81%

+7.63%

1.00

5.845

2.3601

2.8331

+20.72%

+14.94%

1.44

4.415

3.1246

3.7508

+28.65%

+20.96%

The preferred observational regime is intermediate redshift, where the predicted effect is already substantial but spatially resolved disk kinematics are considerably better established than at very high redshift. RC100 provides resolved Hα or CO rotation curves for 100 massive star-forming disks over z∼0.6–2.5 [74], while KROSS provides spatially resolved Hα kinematics for a large sample of typical star-forming galaxies at z∼0.6–1.0 [75]. A decisive implementation would reconstruct the stellar and gas contributions independently of MQ, freeze the expansion history conversion and MQ parameters before inspecting the rotation residuals, and evaluate Eq. (194) without galaxy-specific adjustment. The higher-redshift exploratory regime is treated separately in Sec. 5.6.

5.6 Enhanced gravitational realization during early galaxy formation

Earlier System Frame ages imply a larger MQ acceleration scale and therefore a larger realized gravitational response for the same low-acceleration baryonic source. This is the fixed-source consequence of the cosmic age relation established in Sec. 5.5, not a complete model of early galaxy formation: baryonic assembly, gas dynamics, cooling, star formation, feedback, mergers, and the surrounding cosmological environment evolve concurrently.

greal(R,AU)=gbar(R)1-exp⁡[-gbar(R)/aMQ(AU)].
(196)
greal(R,AU)gbar(R)={1-exp⁡[-gbar(R)aMQ(AU)]}-1.
(197)

Equations (196) and (195) isolate the age dependence of the realized gravitational response at fixed baryonic structure. At a given gbar(R), decreasing AU increases aMQ(AU), decreases the realization fraction, and increases

greal/gbar
(S179)

. The effect is largest where the baryonic acceleration is lowest.

The remaining problem is the evolution of the baryonic source. The baryonic mass distribution must develop through density perturbations and subsequently evolve through expansion, collapse, gas dynamics, cooling, star formation, feedback, and mergers. Equation (196) specifies the age dependence of the realized gravitational response for a prescribed source but does not determine a galaxy’s collapse or assembly history. Formal extrapolation toward AU=0 may remain mathematically informative, but AU=0 is not a realized System Frame epoch because no elapsed-time observable exists there. Physical application begins only after the System Frame age of the universe is defined.

Figure from the preprint

Figure 4. MQ-predicted amplification of realized gravitational acceleration during the first billion years for a fixed enclosed baryonic mass evaluated at a fixed radius of 1000 light-years. The calculation isolates the explicit age dependence of the galactic realization law and does not include the simultaneous MQ mass accretion of Sec. 3.10.

Figure 4 provides a deliberately restricted fixed-source illustration. It holds the enclosed baryonic mass at Mbar=109M⊙ and the realized radius at R=1000 light-years, varying only the explicit age dependence of aMQ in Eqs. (196) and (195). It does not apply the simultaneous mass-accretion construction of Sec. 3.10 and therefore must not be interpreted as the mass history or formation track of a proto-galaxy. A self-consistent formation calculation would have to evolve the baryonic source and its spatial distribution together with the age-dependent galactic realization law.

The observational consequence is comparative. Galaxies at different epochs that can be matched in independently reconstructed Newtonian baryonic acceleration and relevant structural properties should exhibit the age dependence prescribed by Eqs. (196) and (195). A decisive early-galaxy test would require spatially resolved velocity fields extending beyond compact luminous cores, independent reconstruction of stellar and gas mass, and controls for pressure support, inclination, beam smearing, mergers, and outflows.

Current JWST observations demonstrate both the opportunity and the limitation. High-resolution JADES/NIRSpec spectroscopy resolves ionized-gas kinematics in galaxies at 5.5<z<7.4, but the sample includes both rotation- and dispersion-dominated systems and requires forward modeling of instrumental and geometric effects [76]. GN-z11 at z=10.60 shows a resolved velocity gradient that can be modeled as rotation, but the same observations admit an outflow interpretation and do not yet establish a unique circular velocity field [77]. These measurements are therefore exploratory constraints on the high-redshift regime rather than precision tests of Eqs. (196) and (195).

5.7 Long-term evolution of gravitational structure

The same age dependence implies that, for otherwise fixed low-acceleration baryonic configurations, the enhancement of the realized gravitational response decreases toward later cosmic times. This is a comparative prediction of the galactic realization law, not a claim that an individual galaxy retains fixed baryonic structure while the universe evolves.

Existing observations provide a low-redshift consistency check but not a direct test of this prediction. The BUDHIES H I baryonic Tully-Fisher analysis at z≃0.2 finds no significant evolution of the BTFR zero point relative to its consistently constructed local comparison sample [78]. That result is compatible with the comparatively small predicted age-dependent change at z≃0.2, but it is not a direct test of Eq. (196), because the fitted BTFR normalization folds together sample selection, baryonic-mass estimates, velocity definitions, pressure support, and galaxy evolution.

The stronger test is the resolved, frozen-source protocol specified in Sec. 5.5. RC100 and KROSS demonstrate that spatially resolved disk kinematics are observationally accessible over the redshift range in which the predicted age-dependent change becomes substantial [74,79]. Their published dynamical decompositions should not be adopted directly as MQ baryonic inputs where those decompositions already incorporate halo or other fitted dynamical components; the stellar and gas source must instead be reconstructed independently of the MQ residual test.

A statistically secure absence of the prescribed age-dependent shift in appropriately matched, independently reconstructed systems would challenge the cosmic-age dependence of the galactic realization law. Conversely, any apparent redshift trend would require the same controls for baryonic source normalization, observational quality, inclination, and velocity-dependent systematics established by the present-epoch analyses before it could be attributed to temporal evolution.

5.8 Linear cosmological perturbations and the action-level closure test

This subsection establishes that the broader MQ Action admits a linearized perturbative formulation without altering the fixed galactic branch. Its role is limited: the equations below define the field/conservation linearization, stationary source action, timelike completion, and dual first-order spatial formulation needed to pose a future Einstein-Boltzmann calculation. They do not constitute a numerical CMB or large-scale-structure validation.

δGμνsys+Λ̄​δgμνsys+ḡμνsysδΛsys=(8πGc4)δTμνsys,
(198)

where an overbar denotes the homogeneous background and δ the first variation about that background. Equation (198) is the linearization of Eq. (180), not an additional MQ field equation. Likewise, first variation of Eq. (181) gives

∇̄μ​δTsys​μν+δΓμμλ​T̄sys​λν-δΓλμν​T̄sys​μλ=0.
(199)

For scalar perturbations, adopt the standard Newtonian-gauge representation [80–81]. In this subsection only, Φ and Ψ denote the two conventional scalar cosmological metric potentials. With cosmic time t and background scale factor a(t), the perturbed System Frame interval is

ds2=-(1+2Ψc2)c2dt2+a2(t)(1-2Φc2)δijdxidxj.
(200)

For late-time pressureless baryons, let δb denote the baryonic density contrast and ub the peculiar-velocity field. In the linear nonrelativistic subhorizon limit, covariant conservation reduces to the standard continuity and Euler equations,

δ̇b+1a​∇​⋅​ub=0,
(201)
u.b+Hsysub=-1a∇Ψ,
(202)

while the scalar constraint is

∇2Ψ=4πGa2c2​δT00sys.
(203)

Eliminating the peculiar velocity between Eqs. (201) and (200), and then using Eq. (203), yields

δ̈b+2Hsysδ̇b=4πGc2​δT00sys.
(204)

Equation (204) supplies the late-time baryonic growth equation once the realized source perturbation

δT00sys
(S180)

is specified. The stationary weak-field result of Sec. 5.4 permits the density part of that source response to be varied without introducing an additional empirical law. Equation (189) gives the exact local stationary source decomposition. Writing an overbar for a background solution of that branch and varying Eq. (189) gives

δ(T00sysc2)=P¯-1δρb+ρ̄b​δ(P-1)+14πG[∇δ(P-1)⋅∇Φ̄bar+∇P¯-1⋅∇δΦbar].
(205)

Here P is the realization fraction established in Sec. 4.2 and Appx. F, while Φbar is the Newtonian baryonic potential introduced in Sec. 5.4. For a nonzero background baryonic acceleration, its variation follows directly from the established realization fraction,

δ(P-1)=-e-ḡbar/āMQ2P¯2ḡbarāMQ(δgbarḡbar-δaMQāMQ).
(206)

Equations (205) and (204) make the stationary density response explicit wherever the background branch is differentiable. They do not by themselves provide the full relativistic perturbation system, because momentum, pressure, anisotropic stress, and the two scalar metric potentials must still follow from variation of the covariant action.

The corresponding stress-energy variation remains

δTμνsys=12ḡαβsysδgsysαβT̄μνsys-2-ḡsysδ​(δSmsysδgsysμν).
(207)

This distinction is consistent with the Modified gravity literature. Relativistic theories such as TeVeS obtain their background and scalar, vector, and tensor perturbation equations from the complete covariant field content rather than from an isolated MOND acceleration relation [82]. Effective-fluid formulations similarly derive the perturbative density, momentum, pressure, sound-speed, and anisotropic-stress variables from the underlying covariant theory before supplying them to a Boltzmann solver [83]. Cosmological QUMOND simulations provide a useful modified-Poisson spatial analogy but require additional cosmological assumptions and do not themselves constitute a relativistic primordial perturbation theory [84].

The stationary source relation can be placed at the action level without choosing an independent interpolation function. Define the exact MQ response multiplier by

νMQ(y)≡11-exp⁡(-y),​y≡|∇Φbar|aMQ.
(208)

Equation (208) is Eq. (136) written as a response multiplier. A generating function may therefore be defined by quadrature, up to an additive constant that does not affect the field equations,

QMQ(z)≡∫0zds1-exp⁡(-s1/4),​dQMQdz=νMQ(z),
(209)

and the stationary two-potential action may be written

SWFMQ[Φsys,Φbar]=∫dt{-18πG∫d3x[2∇Φ​⋅​∇Φb-aMQ2QMQ​(|∇Φbar|2aMQ2)]-∫d3x​ρbΦ}.
(210)

Its Euler-Lagrange equations are

∇2Φbar=4πGρb,​∇2Φsys=∇⋅[νMQ(|∇Φbar|aMQ)∇Φbar].
(211)
T00sysc2=14πG∇​⋅[νMQ​(|∇Φbar|aMQ)∇Φb]=ρbarP+14πG∇(P-1)​⋅∇Φbar.
(212)

Because νMQ=P-1, Eq. (212) reproduces Eq. (189), including its gradient contribution. In stationary spherical symmetry, with the isolated-source boundary condition, Eq. (211) reduces to

greal=νMQ(gbaraMQ)gbar=gbar1-exp⁡​[-gbar/aMQ],
(213)

and therefore recovers Eq. (136). The stationary galactic realization law is thus reproduced as the Euler-Lagrange response of the two-potential action rather than inserted as a separate phenomenological source prescription. The mathematical organization parallels the bi-potential mechanism used in QUMOND [73], but QMQ is fixed here by Eq. (208), not imported from MOND.

For cosmological spatial evolution in comoving coordinates, the corresponding source equations become

∇x2Φbar=4πGa2δρbar,​∇x2Φsys=∇x⋅[νMQ(|∇xΦbar|a​aMQ)∇xΦbar],
(214)

where ∇x is the comoving spatial gradient. Equation (214) supplies a nonlinear spatial source equation suitable for structure-evolution calculations once embedded in the complete cosmological system; by itself it does not provide the missing relativistic perturbative variables.

A distinct obstruction appears when this purely spatial response is linearized about an exactly homogeneous FLRW background. Since the homogeneous peculiar baryonic field vanishes and the response multiplier is nonanalytic at zero field, the product has the one-sided limit as yo0+

lim|∇Φbar|→0+νMQ​(|∇Φbar|aMQ)∇Φbar=0.
(215)

The vanishing spatial field therefore does not by itself supply a regular nonzero background variable from which every perturbative source contribution can be varied.

The relativistic quantum construction supplies a separate timelike structure. Its System Frame completion is time-oriented, and its WKB phase action obeys the covariant relativistic mass-shell condition derived in Appx. P of Ref. [85]. The conserved phase/current therefore supplies a regular timelike direction even when the homogeneous peculiar spatial gravitational field vanishes. The GR construction independently supplies the homogeneous mapping-normalized expansion background and its System Frame metric [71]. Thus, the zero-field obstruction of Eq. (215) concerns the spatial realization variable, not the absence of all nonzero MQ background structure.

This strategy has an analogue in relativistic MOND constructions, although their additional fields are not imported into MQ. Skordis and Zlosnik construct a covariant theory whose cosmological background contains regular time-directed field structure while recovering MOND behavior in the quasistatic weak-field limit [86], and subsequent perturbative treatments show how the additional relativistic variables enter the CMB and matter spectra [87]. The methodological point is that a nonanalytic quasistatic spatial law need not be Taylor-expanded about zero if the full covariant theory supplies regular background variables.

Given a normalized timelike direction, the metric admits the standard decomposition into that direction and its orthogonal spatial projector [88–90]. Let uμ denote the normalized future-directed System Frame velocity associated with the conserved mapped phase/current. Then

hμν≡gsysμν+uμuνc2,​hμνuν=0,
(216)

and define the projected dimensionless spatial variable

YMQ≡[hμν(∇μΦbar)(∇νΦbar)]1/2aMQ.
(217)

The corresponding response multiplier is

νMQ(YMQ)=11-exp⁡​(-YMQ).
(218)

Replacing the Euclidean spatial contractions in Eq. (210) by contractions with hμν gives a covariant extension compatible with this timelike/spatial decomposition,

SMQcov=-18πG∫d4x-gsys{2hμν(∇μΦsys)(∇νΦbar)-aMQ2QMQ[hμν(∇μΦbar)(∇νΦbar)aMQ2]}+Sb.
(219)

Equation (219) is not asserted to be a new primitive MQ Action or a unique relativistic completion. It is a covariant extension of the already-derived stationary action relative to the existing mapped timelike current. On a stationary hypersurface it reduces to Eq. (210), and hence to Eqs. (211–213).

For bookkeeping, the scalar part may be separated into timelike and spatial contributions,

LMQ=KMQ(XMQ)+JMQ(YMQ).
(220)

The timelike term is constrained by the relativistic phase-action construction rather than introduced as a free cosmological field. Appendix P of Ref. [85] derives the relativistic Hamilton-Jacobi mass shell from the MQ scalar phase-action completion, while Ref. [71] requires realized matter to be constructed from mapping-realized encoded availability, its conserved phase/current, and the emergent System Frame metric. At lowest derivative order, define

XMQ≡-gsysμν(∇μS)(∇νS)m2c2,
(221)

and let ρsrc denote the realized source-density multiplier carried by that mapped matter current. A mass-shell-enforcing timelike contribution is

KMQ(XMQ)=ρsrcc22(XMQ-1).
(222)

Variation with respect to the realized source-density multiplier imposes the mass shell and identifies the normalized phase-flow four-velocity,

XMQ=1,​uμ≡∇μSm,​uμuμ=-c2,
(223)

while current conservation gives

∇μ(ρsrcuμ)=0.
(224)

Metric variation then gives the corresponding timelike source,

TμνK=ρsrcuμuν,
(225)

and, holding the mapped density and phase fixed during metric variation, its on-shell second metric variation is finite,

1-gsysδ2SKδgsysμνδgsysαβ|XMQ=1=ρsrc4(gαβsysuμuν+gμνsysuαuβ),
(226)

so the timelike branch does not inherit the zero-spatial-field singularity of the separated gradient variable.

The remaining problem is specific to the spatial functional. Its small-z behavior is

QMQ(z)=43z3/4+12z+115z5/4+O​(z7/4),​QMQ''(z)=-14z-5/4+O​(z-3/4),
(227)

which shows that the ordinary Hessian with respect to the baryonic-gradient variable diverges at YMQ=0 while the exact deep-realization response is retained. This is a statement about that primal spatial coordinate representation, not a singularity of the full variational theory. As a function of the projected gradient magnitude, the associated spatial energy density admits the corresponding convex dual representation. Fenchel-Legendre duality is the standard mathematical tool for changing variational coordinates in singular or degenerate nonlinear elliptic problems [73,91].

Define the dimensionless dual spatial field from the already-fixed realized response,

WMQμ≡νMQ(YMQ)hμν∇νΦbaraMQ,​WMQ≡(hμνWMQμWMQν)1/2=YMQνMQ(YMQ).
(228)

Its dual generating function is

RMQ(W)≡2W​Y(W)-QMQ(Y(W)2),​W=YνMQ(Y),
(229)

and the exact Legendre identity becomes

aMQ2QMQ(YMQ2)=2aMQWMQμ∇μΦbar-aMQ2RMQ(WMQ).
(230)

The dual action is therefore

SMQdual=-18πG∫d4x-gsys{2hμν(∇μΦsys)(∇νΦbar)-2aMQWMQμ∇μΦbar+aMQ2RMQ(WMQ)}+Sb+SK.
(231)

Variation with respect to the dual field returns Eq. (228) and eliminating that field reproduces Eq. (219); the transformation changes variational coordinates rather than the stationary MQ response.

The regularity improvement is visible in the zero-field series. Combining Eq. (227) with WMQ=YMQνMQ(YMQ) gives

RMQ(W)=23W3-12W4+1330W5+O(W6).
(232)

Thus, the dual spatial density is twice continuously differentiable at WMQ=0: its first and second derivatives vanish there. Exact FLRW homogeneity also has ∇μΦbar=0 and

WMQμ=0
(S181)

, so the metric response of the dual spatial part is finite,

δ2SMQdual,spδgsysμνδgsysαβ|FLRW=0,
(233)

and the divergent QMQ″(0) of the primal baryonic-gradient coordinate is not a divergence of the on-background dual action.

The inverse relation at small field is

Y(W)=W2-W3+1312W4-54W5+O(W6).
(234)

This completes the mathematical zero-field regularization of the spatial variational coordinates. It does not constitute a numerical CMB or large-scale-structure validation, nor does it turn the referenceability fractions of Table 5 into independently conserved cosmological fluids.

The remaining calculation is numerical and dynamical: Eqs. (198–199), the timelike system (219–224), and the regular first-order spatial formulation (226–229) must be linearized jointly about the homogeneous MQ background with physically admissible initial conditions and integrated in CLASS or an equivalent Einstein-Boltzmann solver [80–81,92]. The resulting background history D(z), fσ8(z), matter spectrum P(k), CMB TT/TE/EE spectra and lensing, and BAO distances are the direct cosmological observables. Numerical cluster abundance is a later nonlinear test. Agreement must emerge from one integrated solution before the numerical correspondences in Table 5 can acquire dynamical evidentiary weight.

The first numerical implementation should be prediction-first rather than fraction-matching. The MQ background construction and action-derived perturbation system must be fixed before comparison with cosmological likelihoods. Planck-inferred Ωb, Ωc, or ΩΛ values must not be inserted as physical identifications of MQ referenceability fractions merely because their percentages are numerically similar. Standard photon, baryon, and neutrino microphysics may enter through independently specified conventional equations, while any quantity not derived within MQ must be identified explicitly as an external initial-condition or astrophysical input. The code should first verify background conservation and gauge consistency, then calculate the linear transfer functions, and only afterward compare the frozen prediction with CMB spectra and lensing likelihoods [11,92], BAO measurements [93], and late-time growth observables. Failure at any of those stages would delimit the cosmological embedding even if the algebraic domain fractions of Sec. 3.9 remain unchanged.

5.9 Open problems and required future work

The remaining program is narrow and independently falsifiable. The principal galactic result rests on the completed Milky Way, SPARC, and LITTLE THINGS calculations, subject to the source model and observational limitations established in Sec. 4.10. The highest-value near-term extension is the frozen intermediate-redshift rotation curve test of Sec. 5.5, which directly challenges the predicted age dependence without introducing a galaxy-specific MQ calibration parameter [74,79].

The nearest unresolved present-epoch source test is the residual-blind resolved stellar profile calculation specified in Sec. 4.10 for DDO 50, DDO 216, IC 1613, NGC 1569, and NGC 3738. The published radial stellar mass distributions of Zhang et al. [57], together with the LITTLE THINGS gas and kinematic products of Oh et al. [5], permit the stellar source to be reconstructed radially rather than represented by a single integrated normalization. The MQ acceleration scale and galactic realization law must remain frozen before the residuals are inspected. A substantial reduction of the existing stellar-important residuals would localize the discrepancy to baryonic source reconstruction. A persistent coherent residual would instead place direct pressure on the fixed galactic realization law. Until this calculation is completed, the stellar-important LITTLE THINGS galaxies remain source model limited rather than precision tests.

The independent relativistic requirements remain more demanding. For ESO 325-G004, Sec. 5.3 requires the observed extended stellar source to be propagated through the covariant System Frame metric, with both weak-field metric potentials and their gravitational slip determined before the deflection, convergence, and shear maps are compared with the Einstein-ring data. The existing spherical aperture calculation is only a benchmark. For 1E 0657-558, Sec. 5.4 requires nonlinear evolution from independently specified baryonic merger conditions through the full MQ field equations, followed by lensing through the evolved metric. The calculation must reproduce the observed lensing-gas displacement, merger morphology, and shock-compatible dynamical state simultaneously. No collisionless dark component or phenomenological gravitational response delay is available as an adjustable substitute [68–70].

Cosmological closure likewise remains a numerical test rather than an algebraic correspondence. Section 5.8 supplies the action-derived perturbative variables and regularized spatial formulation needed to pose the calculation, but the homogeneous background and coupled perturbations must still be evolved in an Einstein-Boltzmann solver and confronted jointly with the CMB temperature and polarization spectra, CMB lensing, BAO distances, matter power spectrum, and late-time growth. The numerical similarity between the mapping-derived domain fractions of Sec. 3.9 and standard cosmological density fractions is not evidence for dynamical equivalence unless the same frozen MQ evolution reproduces those observables.

These tests delimit the scope of the present result. Galaxy-scale agreement does not establish a universal replacement for nonbaryonic dark matter, and failure of an unfinished lensing, cluster, or cosmological extension does not erase a separately tested galactic relation. It instead identifies where that extension fails. Conversely, successful rotation curve fits cannot substitute for unresolved relativistic and cosmological calculations. The present evidentiary claim is therefore restricted to the fixed galactic realization law and its completed observational tests, while the resolved stellar source, redshift-evolution, extended-lensing, cluster merger, and Einstein-Boltzmann calculations provide distinct routes by which the broader MQ interpretation can be extended or rejected.

APPENDIX A: DEFINITIONS OF KEYWORDS, TERMS AND SYMBOLS

Notes

In earlier publications (prior to 2026), the roles of the Internal Frame and System Frame were described using a dimension-based distinction. With the formal identification of the Frames mapping, this characterization has been refined to align with the kinematic structure defined in the present work. The Internal Frame is now understood strictly as a discrete, count-based configuration domain without intrinsic geometric or continuous measure, while the System Frame represents the relational encoding in which observables are realized following application of the Frames mapping.

This refinement does not constitute a reversal of physical interpretation, but rather a clarification: geometry, continuity, and measurement arise only through the mapping from discrete count structure to relational observables. The realization of observables corresponds to the mapping-normalized encoding of conserved count, as defined by the mapping normalization functional. In this sense, normalization specifies the conditions under which Internal Frame configurations admit consistent, repeatable representation in the System Frame, and thereby defines the operational meaning of measurement.

Keywords

Internal Frame
The Internal Frame is the discrete, count–based configuration domain in which all fundamental measures are defined relative to the fundamental units lf, mf and tf. All relational structure is encoded as integer count configurations over this domain. The Internal Frame does not possess a continuous geometric structure; rather, geometry arises only through its mapping into the System Frame.

System Frame
The System Frame is the non–discrete, relational encoding of Internal Frame structure induced by the Frames mapping. It does not possess an independent fundamental measure but instead represents observable quantities as continuous relational separations derived from discrete count structure. All physically meaningful observables are defined in the System Frame.

Frames mapping
The Frames mapping is the kinematic transformation that relates discrete count configurations in the Internal Frame to continuous relational structure in the System Frame. This mapping preserves realization consistency and enforces closure conditions on relational configurations. Observable quantities arise only after application of the Frames mapping, and all geometric and interaction structure emerges from its properties.

Frames Principle
The Frames Principle asserts that all physically meaningful observables arise from the Frames mapping between the Internal Frame and the System Frame. No observable quantity is intrinsic to the Internal Frame alone. All measurement is relational and defined through the mapping. This principle replaces the assumption of a pre–existing spacetime manifold with a kinematic origin of geometry.

MQ Action
The MQ Action is the fundamental kinematic functional defined over the Internal Frame count structure. It governs the evolution and consistency of relational configurations under the Frames mapping. The action admits distinct realizations determined by count separation: fixed interaction–scale realizations, corresponding to electromagnetic behavior at the electromagnetic demarcation, and count–resolved realizations, corresponding to gravitational behavior through the mapping–induced metric structure. In both cases, observable dynamics arise from the normalization and realization of conserved count under the mapping. All emergent physical behavior, including electromagnetic interactions and gravitational curvature, therefore arises from this single action without the introduction of any independent dynamical components.

MQ variational principle
The MQ variational principle states that physically realized configurations in the System Frame correspond to stationary points of the MQ Action functional defined over admissible realizations of encoded availability derived from the Internal Frame, thereby replacing the independent postulation of dynamical laws with a single extremization condition from which all effective interactions emerge. Formally, let M:ω→R denote the Frames mapping, and let ω'⊂ω be a relational region of the Internal Frame with realized configuration Ω=M(ω')⊂R. Let ρR denote the corresponding coarse–grained encoded availability field. The MQ Action is defined as the functional

S[ρR]=∫RL(ρR,∂ρR,…)dμ,
(S182)

where L is constructed from count–based relations induced by the mapping M, and dμ is the induced measure on R. The MQ variational principle asserts that admissible realizations satisfy δS=0 under variations consistent with the constraints of the Frames mapping and the conservation relations defined on ω'.

Informativity differential
The Informativity differential is the non–Lorentz length contraction associated with the discrete count structure of the Internal Frame. It arises from the normalization of conserved count under the Frames mapping and quantifies the difference between relational separation in the Internal Frame and its realized representation in the System Frame. The magnitude of the effect is count–dependent, increasing at small separations and diminishing toward the upper–count limit. The Informativity differential is distinct from curvature and represents one of the two geometric consequences of realization in the System Frame, the other being the emergent curvature associated with count loss under successive updates.

Metric differential
The metric differential characterizes the discrete–continuous mismatch that arises when conserved Internal Frame count structure is realized as a continuous observable in the System Frame under the Frames mapping. It does not represent a physical deformation in itself but rather the necessity of expressing integer–resolved count relations as normalized, non–discrete quantities.

For a fundamental measure defined by a count n in the Internal Frame, realization in the System Frame requires application of a mapping normalization functional. This process preserves count conservation while enforcing continuity, introducing a representation–level mismatch between the discrete count structure and its continuous encoding. The metric differential refers to this mismatch as a structural consequence of the mapping.

The metric differential is therefore distinct from the Informativity differential. While the Informativity differential describes a non–Lorentz length contraction in realized separations, the metric differential captures the underlying requirement that discrete count relations be expressed in a continuous form. In this sense, it provides the count–level representation of the discrete–continuous correspondence enforced by the Frames mapping, whereas the Informativity differential describes its realized geometric manifestation in the System Frame.

Electromagnetic demarcation
The electromagnetic demarcation is the fixed interaction count distance nL = 276/θsi, corresponding to approximately 84.60055lf, at which known electromagnetic interactions are realized in MQ. At this scale, the effect of the Informativity differential is still significant. Its influence increases with decreasing interaction distance and rapidly diminishes with increasing distance, becoming physically negligible beyond a few hundred counts of lf. The electromagnetic demarcation is the fixed electromagnetic interaction scale at which that discrete measure effect is realized.

Count separation
The number of fundamental length units lf between relational configurations in the Internal Frame. This quantity determines the scale at which mapping–induced deformation becomes significant.

Mapping–induced deformation
The modification of realized separation in the System Frame arising from the discrete–continuous Frames mapping between the Internal Frame and the System Frame. This deformation is fully determined by the mapping structure and does not constitute an independent physical field. It has two coupled manifestations depending on the scale of the interaction. At fixed count separation, it appears as a local, scale–invariant modification governed by the Informativity differential, corresponding to the non–Lorentz contraction of realized lengths. Over extended count separation, the same mapping–induced deformation accumulates through the mapping structure and is realized as geometric curvature in the System Frame. These two manifestations are not distinct mechanisms but represent different scale–dependent realizations of a single underlying mapping effect.

Fixed interaction scale
A fixed interaction scale is a count–defined interaction distance at which a physical interaction is realized under the Frames mapping. For electromagnetic interactions, this occurs at the electromagnetic demarcation nL = 276/θsi, where the interaction is evaluated at a single, fixed count separation. At this scale, both the Informativity differential (non–Lorentz length contraction) and the corresponding curvature are realized with fixed values determined by the mapping. The interaction therefore reflects a specific point in the discrete count structure rather than a limiting regime.

Curvature–resolved regime
A curvature–resolved regime refers to the dependence of realized geometric structure on count separation under the Frames mapping. In this description, curvature is a direct consequence of the mapping evaluated at a given count distance. Different separations correspond to different realized curvature values in the System Frame, providing a geometric description of interaction that varies with count scale. This behavior characterizes gravitational phenomena, where measurements probe a range of separations rather than a fixed interaction scale.

Bounded realization
Bounded realization is the MQ condition in which the gravitational relation available from a physical source is realized in the System Frame only to the fraction permitted by the count-constrained Frames mapping. The realized fraction is bounded between zero and unity and approaches complete realization as the baryonic acceleration becomes large relative to the MQ acceleration scale. At low baryonic acceleration, incomplete realization produces a realized gravitational response larger than that inferred by applying the Newtonian source relation alone. Bounded realization does not introduce additional gravitating matter or modify the baryonic source law; it describes how the available gravitational relation is realized through the mapping.

Realization fraction
The realization fraction is the dimensionless fraction Pof gravitational availability from a baryonic source that is realized through the Frames mapping in the System Frame. For the galactic gravitational realization derived in MQ,

P(R)=1-exp⁡[-gbar(R)aMQ],
(S183)

where gbar(R)is the Newtonian baryonic acceleration and aMQis the MQ acceleration scale. The realization fraction approaches unity in the high-acceleration limit and decreases continuously below unity as the available baryonic acceleration decreases. The realized gravitational response is consequently

greal=gbar/P
(S184)

.

Galactic realization law
The galactic realization law is the fixed MQ relation that maps an independently determined Newtonian baryonic gravitational source to its realized System Frame acceleration or orbital velocity through the realization fraction and the first principles acceleration scale aMQ. The law contains no galaxy-specific MQ fitting parameter and preserves the Newtonian baryonic source relation while modifying the fraction of that relation realized through the Frames mapping. It approaches Newtonian baryonic dynamics continuously at high acceleration and produces the enhanced realized response characteristic of low-acceleration galactic dynamics.

Realized gravitational response
The realized gravitational response is the gravitational acceleration, potential, metric response, or associated motion produced in the System Frame after gravitational availability from the physical source has been realized through the Frames mapping. It is distinguished from the Newtonian baryonic source quantity, which specifies the gravitational availability supplied by the reconstructed baryonic mass distribution. In the stationary galactic limit, the realized response is related to the baryonic source through the realization fraction, while in the covariant description it is represented through the System Frame metric and its action-derived stress-energy source.

Gravitational availability
Gravitational availability is the gravitational relation supplied by the physical source before its realization through the Frames mapping. In the weak-field galactic description it is represented by the independently reconstructed Newtonian baryonic quantities Mbar, ρbar, Φbar, and gbar. Gravitational availability is therefore a source property rather than an additional material component. MQ distinguishes this available source relation from the gravitational response realized in the System Frame, with their relation determined by bounded realization.

Symbols

A.1 Primitive definitions

This subsection defines the irreducible count-based quantities that constitute the Internal Frame. These objects form the non-geometric substrate of MQ and serve as the primitive inputs to all subsequent constructions under the Frames mapping.

R3 — three–dimensional System Frame configuration space for single–subsystem spatial coordinates r.


A.2 Discrete Internal Frame primitives

This subsection establishes the primitive combinatorial and relational structures of the Internal Frame count substrate. These objects define admissible count configurations, adjacency relations, conservation structure, and regional count evolution prior to realization under the Frames mapping.

A — The Internal Frame update index is a discrete ordering parameter labeling successive configurations of the Internal Frame count medium. It indexes deterministic updates of relational count assignments nI(ξ,A) across the configuration domain ω, with regional evolution defined on subsets ω'⊂ω through conservation and adjacency relations. The update index A is not a physical time coordinate and carries no metric, geometric, or causal structure. It does not correspond to System Frame time t, nor to any observable temporal parameter. Instead, it provides an intrinsic ordering of count redistribution steps within the Internal Frame. Evolution with respect to A defines the deterministic conservation and redistribution of count via adjacency relations, as expressed in Internal Frame conservation laws of the form

∂ρI(ξ,A)∂A+div⁡adjJI(ξ,A)
(S185)

= 0. The mapping between the Internal Frame update index A and System Frame time t is established only through the Frames mapping M, and is generally nontrivial. Consequently, A should be interpreted as a kinematic sequencing parameter intrinsic to the discrete count substrate, rather than as a temporal observable.

A.3 Frames Mapping and Realization Structures
This subsection defines the relational structures associated with admissible Frames mappings and referenceable universe domains. These objects govern normalization, realization, holonomy structure, record encoding, and the consistency conditions required for emergent System Frame observables.

𝒫(R) — is the dimensionless gravitational realization fraction at galactocentric radius R. It specifies the fraction of baryonic gravitational availability realized through the Frames mapping in the System Frame and is defined by

P=1-exp​[-gbar/aMQ]
(S186)

. For finite positive baryonic acceleration, 0<P<1, with P→1 in the high-acceleration limit. The realized acceleration satisfies

greal=gbarP
(S187)

.

νMQ(y) — is the MQ realization multiplier, defined by

νMQ(y)=[1-exp(-y)]-1
(S188)

. It is the inverse of the realization fraction when

y=gbaraMQ
(S189)

and maps the baryonic gravitational field to the realized gravitational response. The function is fixed by the MQ realization law rather than introduced as an empirical interpolation function.

QMQ(z) — is the scalar primitive of the MQ realization multiplier used in the weak-field variational formulation. It is defined by

dQMQdz=νMQ(z)
(S190)

, with normalization fixed by the lower integration bound. The function is derived from the established MQ realization law and is not an independent phenomenological function.

YMQ — is the dimensionless covariant spatial realization invariant formed from the norm of the baryonic potential gradient projected orthogonally to the timelike flow and normalized by aMQ. It provides the covariant scalar argument of νMQ in the spatial realization action.

YMQμ
(S191)

— is the auxiliary MQ dual realization vector introduced in the exact first-order representation of the spatial weak-field action. It is determined algebraically from YMQ, νMQ, the spatial projector, and the gradient of Φb. It introduces no independent physical degree of freedom and may be eliminated to recover the original realization action exactly.

WMQ — is the dimensionless invariant magnitude of the MQ dual realization vector

YMQμ
(S192)

. It satisfies WMQ=YMQνMQ(YMQ) and provides the natural scalar variable of the dual weak-field representation.

RMQ(W) — is the exact dual realization function associated with QMQ(Y2). It is defined through the Legendre transformation relating W and Y and is constructed such that eliminating the auxiliary dual field returns the original spatial realization action without changing the MQ gravitational realization law.

XMQ — is the dimensionless timelike mass-shell invariant formed from the inverse System Frame metric and the gradient of the action phase S, normalized by m²c². On the physical timelike mass-shell branch, XMQ=1.

KMQ(XMQ) — is the timelike constraint contribution to the covariant MQ realization Lagrangian. It enforces the mass-shell condition XMQ=1 through the associated source-density multiplier and generates the corresponding timelike stress-energy contribution.

JMQ(YMQ) — is the spatial realization contribution to the covariant MQ realization Lagrangian. It contains the bounded-realization dependence on the covariant spatial invariant YMQ and complements KMQ in the decomposition of the full realization Lagrangian.

LMQ — is the covariant MQ realization Lagrangian constructed from the timelike and spatial contributions according to LMQ=KMQ(XMQ)+JMQ(YMQ). It embeds the stationary bounded-realization law within the covariant System Frame action.

A.4 Emergent System Frame Continuum Structures
This subsection establishes the continuum relational structures induced through the Frames mapping. These objects define encoded availability densities, emergent coordinate representations, continuum observables, and associated differential structures within the realized System Frame.

r — is the System Frame spatial coordinate vector r∈R3 of the realized System Frame continuum. Used to parameterize encoded availability densities, probability densities, currents, velocity fields, wavefunctions, and other nonrelativistic observables arising from the Frames mapping. Appears in expressions such as ρS(r,t), P(r,t), J(r,t), v(r,t), Q(r,t), and ψ(r,t). Represents position within the realized three-dimensional spatial domain Ω. Distinct from the spacetime coordinate xμ used in relativistic and geometric formulations.

xμ — is the System Frame spacetime coordinate xμ=(x0,x1,x2,x3). Four-dimensional coordinate on the emergent System Frame spacetime manifold. Used in continuum geometric constructions including the metric tensor

gμνsys
(S193)

, connection coefficients

Γμνρ
(S194)

, curvature tensors, geodesics, action functionals, and covariant field equations. Coordinates are defined on the realized manifold induced by the Frames mapping. The spatial components of xμ reduce locally to the realized spatial coordinates represented by r, but xμ is reserved for relativistic and differential geometric formulations.

t — System Frame time coordinate used in the continuum encoding and Schrödinger evolution.

Φsys — is the gravitational potential realized in the System Frame in the stationary weak-field limit. Its spatial gradient determines the corresponding realized gravitational response. It is the weak-field potential obtained from the realized System Frame metric rather than an independently introduced gravitational field.

Msys(<R) — is the enclosed gravitational source associated with the realized System Frame weak-field response within galactocentric radius R. It is derived from the mapped baryonic gravitational availability and does not represent an independently postulated dark matter mass component.

ρsrc — is the scalar mass-density multiplier associated with the timelike source branch of the covariant realization action. It enforces the mass-shell condition through KMQ, enters the conserved current ρsrcuμ, and generates the corresponding on-shell stress-energy tensor. It is distinct from the canonical MQ geometric density ρsys.

hμν — is the inverse spatial projection tensor associated with the timelike System Frame four-velocity uμ, defined by

hμν=gsysμν+uμuνc2
(S195)

. It projects covectors onto the spatial hypersurface orthogonal to uμ and satisfies hμνuν=0.

uμ — is the normalized timelike four-velocity associated with the System Frame source branch. It satisfies uμuμ=-c2 and, in the MQ mass-shell construction, is generated from the action phase according to

uμ=∇μSm
(S196)

.

A.6 Differential and operator structure
This subsection defines the operators acting on realized quantities in the System Frame. These include mapping-consistent derivative and update structures that enable the construction of curvature and field equations from the emergent metric.

∇ — System Frame gradient operator acting on continuum–encoded fields.

∇⋅ — System Frame divergence operator acting on geometric vector fields.

A.7 Fundamental constant and measures
This subsection collects the invariant constants and normalization measures that appear in the MQ framework. These quantities fix the scale of mapping between discrete counts and continuum observables without introducing independent phenomenological parameters.

2QLnL — also known as the Informativity differential describes the length contraction associated with discrete measure. It is a dimensionless normalization modifier arising from the discrete count structure of the Internal Frame, representing the deviation between integer–resolved count measures and their realized System Frame encoding under the Frames mapping M:ω→R. It quantifies how conserved Internal Frame count NI(ω',A) is redistributed and normalized when expressed as a System Frame observable.

lf,mf,tf — are the fundamental measures, more precise expressions for Planck’s units–length, mass, and time–that consider the effects of length contraction associated with discrete measure.

2QLnL — also known as the Informativity differential describes the length contraction associated with discrete measure. It is a dimensionless normalization modifier arising from the discrete count structure of the Internal Frame, representing the deviation between integer–resolved count measures and their realized System Frame encoding under the Frames mapping M:ω→R. It quantifies how conserved Internal Frame count NI(ω',A) is redistributed and normalized when expressed as a System Frame observable.

G — is the gravitational constant, G=6.67407794(17) × 10–11 m3 kg–1 s–2 [25] such that its value considers the effects of length contraction associated with discrete measure at the upper count limit. Italicized G identifies a measure (e.g., G=6.67384483(17) × 10–11 m3 kg–1 s–2 [26] at the blackbody demarcation ([9], Appx. AB)).

θsi — can be measured as the polarization angle of quantum entangled X–rays at the degenerate frequency of a maximal Bell state. As an angle θsi=3.26239 rad ± 2 μrad [27]. As a momentum θsi = 3.262390305(89) kg m s–1. And with respect to the System Frame, θsi has no units.

c is the speed of light which may also be written as c=nLlf/nTtf =299,792,458 m/s [28] such that nL=nT=1 is physically significant.

Hsys — denotes the System Frame realization of the dimensionful expansion rate, expressed in km s-1 Mpc-1, corresponding observationally to the CMB-resolved expansion rate formerly denoted HCMB. It is the non-discrete System Frame realization of the scalar expansion under mapping normalization and is given by

Hsys=(π/(θsiAU))(kmMpc-1/(s/109))
(S197)

.

Hint — denotes the Internal Frame realization of the dimensionful expansion rate, expressed in km s-1 Mpc-1, corresponding observationally to the Cepheid/distance ladder expansion rate formerly denoted Hcep. It is the discrete frame realization of the scalar expansion under the complementary mapping normalization and is given by

Hint=((2θsi-π)/(θsiAU))
(S198)

(kmMpc-1/(s/109))
(S199)

.

H — denotes the dimensionful local expansion rate expressed in observational units, typically km s-1 Mpc-1. It is the unit-resolved expansion rate used for observational comparison before specifying a frame realization such as Hsys, Hint, or Hb.

H — denotes the dimensionless scalar expansion magnitude associated with the increasing radius of the universe. The non-italic form H is used intentionally to distinguish this scalar from the italicized dimensionful expansion rates. For the realized universe, H=θsi.

lp, mp and tp — are Planck’s units for length, mass, and time.

nL, nM and nT — are physically significant discrete counts of lf, mf and tf respectively.

nL — describes the count of lf representative of a change in position of an observable measured with respect to the observer’s frame of reference.

nTu — is a count of tf equal to the age of the universe.

nLu is a count of lf equal to the diameter of the universe.

nθ — is the System Frame count of θsi resolved as a function of the Internal Frame demarcation adjusted for the metric differential (i.e., ⌊84.60055/2⌉ = 42).

ħ — is the reduced Planck constant, 1.054571817 10–34 m2 kg s–1. When accounting for the Informativity differential at the upper count bound, this term is not italicized (i.e., ħ=1.054534984(26) 10–34 m2 kg s–1).

Aee — is the dilated age of the universe as measured from the current expansionary epoch. Its value is a function of expansion with respect to the total elapsed time representing the quantum epoch.

Aqe — is the non–dilated age of the universe corresponding to the quantum epoch.

Ωvis — is that domain of the universe that is presently visible.

Ωobs — is that domain of the universe that is presently visible or will be visible given infinite elapsed time.

Ωdk — is that domain of the universe that will never be visible.

Ωuobs — is that domain of the universe that will eventually be visible but is not presently visible.

Ωtot — is the entire domain that makes up the universe (i.e., 1 or 100%).

Ωf — is that domain of the universe that is fundamental. Its physical significance describes a pivot point with respect to the remaining domains.

Mtot — is the total mass of the universe, including all forms of energy.

Macr — is the rate of mass accretion in the universe.

me — is the mass of an electron.

a0 — is the ground state orbital of an atom, also known as the Bohr radius.

α— is the fine structure constant, a dimensionless fundamental physical constant that characterizes the strength of the electromagnetic interaction between electrically charged particles. The recommended numerical value is published by CODATA based on the least-squares adjustment of the best available experimental and theoretical data.

αf–1, αp–1, and αc–1 — are the fundamental, Planck–like, and classical forms of the inverse fine structure constant. αf–1=42θsi is defined with respect to the Internal Frame; the remaining two against the System Frame. αc–1 differs from αp–1, accounting for the length contraction associated with discrete measure.

γ — is the conventional electromagnetic transformation parameter used in the gamma approach to resolve the electromagnetic constants. It has a value of 1.1996724235 m6 C-1 s-4 and is determined from an independent source measurement of one of the electromagnetic constants (typically ke, ε0, or μ0) and subsequently used to derive the remaining electromagnetic quantities while preserving the independence of the selected source measurement.

γF — is the fundamental gamma, a first principles realization parameter that replaces the conventional gamma transform within the MQ framework. Unlike γ, which is resolved from experimentally measured electromagnetic constants, γF is derived entirely from fundamental quantities and serves as the canonical electromagnetic realization parameter from which the electromagnetic constants are independently determined.

Λ — describes Einstein’s cosmological constant.

e — is the elementary charge, the magnitude of the electric charge carried by a proton (or equivalently the magnitude of the negative charge carried by an electron). In the SI, e = 1.602 176 634 x 10-19 C exactly.

TCMB — is the temperature of the Cosmic Microwave Background.

Mvis — is the total mass corresponding to the visible domain Ωvis, representing mass contained within the portion of the universe presently visible from a given reference point.

Mobs — is the total mass corresponding to the observable domain Ωobs, representing mass contained within the portion of the universe that is presently visible or will become visible given infinite elapsed time.

Muobs — is the total mass corresponding to the unobserved domain Ωuobs, representing mass contained within the portion of the universe that is not presently visible but will become observable given sufficient elapsed time.

Mtot — is the total mass of the universe, including all mass-equivalent forms of energy.

Mf — is the mass corresponding to the fundamental domain Ωf. It defines the mass threshold associated with the fundamental partition of the universe, corresponding to the limiting rate at which distinguishable mass events can be resolved at a reference point.

M⊙ — denotes one solar mass.

nMu — is a count of mf equal to the total mass-equivalent content of the universe.

nTu — is a count of tf equal to the age of the universe.

nLu — is a count of lf equal to the diameter of the universe.

vNew — is the Newtonian orbital velocity predicted from the adopted baryonic mass model for a mass relative to a galaxy's center of mass.

vo— is the published observed orbital velocity of a mass relative to a galaxy's center of mass.

vreal— is the MQ realized orbital velocity of a mass relative to a galaxy's center of mass, obtained from the bounded realization of the Newtonian baryonic gravitational response under the Frames mapping.

aMQ — is the MQ acceleration scale governing bounded gravitational realization. It determines the transition between nearly complete gravitational realization at high baryonic acceleration and incomplete realization at low baryonic acceleration. In the galactic realization law it enters through the dimensionless ratio gbar/aMQ. Its present-epoch value is derived from MQ first principles rather than calibrated to galaxy rotation data.

A.8 Expansion and cosmological scalars
This subsection defines the scalar quantities governing expansion in MQ, including mapping-derived normalization factors and cosmological terms. These quantities describe the homogeneous expansion as it emerges from count-based structure under the Frames mapping.

QL — denotes the mapping-induced residual count contribution entering the realized separation functional. In the right-triangle count construction, nL gives the Internal Frame count separation, while nL + QL gives the realized hypotenuse whose continuum encoding defines the System Frame separation. QL therefore replaces the older Internal-to-System scale-ratio notation by expressing expansion and curvature through count-resolved residual structure rather than through a global ratio.

AU — is the System Frame age of the universe. The quantity AU(z) denotes the elapsed universe age at redshift z, as determined from the System Frame cosmological expansion history. Its present-epoch value is denoted AU,0.

AU,0 — is the present-epoch System Frame age of the universe, defined by AU,0=AU(0). It provides the reference age used in the redshift dependence of aMQ and related cosmological realization quantities.

A.9 Mapping and normalization symbols
This subsection collects symbols associated with mapping normalization and the single-action realization of dynamical quantities under the Frames mapping. It includes the functions and parameters that convert dimensionless scalar quantities into their observable, unit-resolved forms and govern admissible dynamical variation.

ΔH — denotes the mapping-induced normalization offset between System Frame and Internal Frame realizations of the scalar expansion parameter. It is defined by ΔH=θsi-π and represents a frame-resolution offset under the Frames mapping, not an independent dynamical parameter or a global count-expansion mismatch.

H— denotes the dimensionless scalar rate of radial expansion defined at the level of the System Frame realization of the Frames mapping. It represents the intrinsic expansion magnitude governing the growth of realized separations and is independent of observational units. Within MQ, H = θsi, so that H is not a measured local rate but the mapping-defined scalar controlling the expansion geometry and associated mapping-consistent scalar relations.

Hsys— is the System Frame realization of the expansion rate, obtained by expressing the scalar H under mapping normalization consistent with the non-discrete relational encoding. It corresponds to expansion values inferred from CMB-based measurements and represents the continuum realization of the underlying scalar expansion.

Hint— is the Internal Frame realization of the expansion rate, obtained by expressing the scalar H under mapping normalization consistent with discrete count structure. It corresponds to expansion values inferred from local distance ladder measurements and differs from Hsys by a mapping-induced normalization offset rather than by an independent physical parameter.

H — denotes the dimensionful local expansion rate expressed in observational units (e.g., km s-1 Mpc-1), prior to specification of a particular frame realization. It represents the normalized, unit-resolved form of the scalar expansion parameter H and serves as the general observable expansion rate from which specific realizations such as H0, Hsys, and Hint are defined.

A.10 Field equation components
This subsection defines the geometric and source terms entering the MQ field equation. These quantities are constructed from the emergent metric and mapping-consistent action, ensuring that curvature, stress-energy, and coupling arise from a single unified framework.

ρc— is the critical density

ρc=3H2/8πG
(S200)

is the total energy density required for spatial flatness in the standard Friedmann–Lemaître cosmological model. It provides the conventional normalization scale against which fractional density components are defined.

ΩΛ— is the fractional scalar density parameter

ΩΛ=ρΛ/ρc=Λc2/3H2
(S201)

specifies the proportion of the total cosmic energy density attributed to the scalar contribution. Within MQ it is determined entirely by invariant internal relations.

Ωb— is the baryonic matter fraction (ΛCDM comparison bookkeeping).

Ωc— is the cold dark matter fraction (ΛCDM comparison bookkeeping).

gμνsys
(S202)

— denotes the emergent metric tensor in the realized System Frame. It defines the local continuum metric structure associated with mapped measurable separations under the Frames mapping. Unlike standard general relativity, where gμν is treated as a fundamental spacetime field, MQ treats

gμνsys
(S203)

as a derived object obtained through coarse–graining of the realized separation functional. Its normalization is fixed by mapping constraints and is not independently specified.

Gμνsys
(S204)

— denotes the Einstein tensor in the realized System Frame. It is constructed from the emergent metric tensor

gμνsys
(S205)

through the associated Levi–Civita connection and resulting curvature contractions. In MQ,

Gμνsys
(S206)

represents the continuum realization of higher-order discrete variation induced by the Frames mapping. Its divergence-free property follows from the metric-compatible connection and associated Bianchi identity, rather than from an independently postulated geometric structure.

Tμνsys
(S207)

— denotes the stress–energy tensor in the realized System Frame. In MQ, it is defined variationally from the MQ Action through metric variation,

Tμνsys=-2-gsysδSmatterδgsysμν
(S208)

, and therefore encodes the response of realized normalized count density to variations in the emergent metric structure. It is not independently specified and inherits its structure from encoded availability and mapping normalization.

Λsys — denotes the realized cosmological term appearing in the System Frame field equation. In MQ, it is not a fundamental constant but a derived scalar arising from the same realized separation structure that defines the emergent metric tensor

gμνsys
(S209)

. It reflects the large-scale normalization of the System Frame geometry induced by the Frames mapping and is therefore fixed by the underlying count geometry rather than introduced as an independent phenomenological parameter.

gbar(R) — is the Newtonian baryonic gravitational acceleration at galactocentric radius R, calculated from the independently specified baryonic mass distribution before bounded realization through the Frames mapping. It represents the gravitational availability supplied by the baryonic source and is the source acceleration entering the MQ realization fraction.

greal(R) — is the gravitational acceleration realized in the System Frame at galactocentric radius R after bounded realization of the baryonic gravitational availability. In the stationary galactic limit it is determined from gbar, aMQ, and the realization fraction without a galaxy-specific MQ fitting parameter.

Φb — is the baryonic gravitational potential generated by the independently specified baryonic source distribution before bounded realization. Its spatial gradient determines the Newtonian baryonic acceleration gbar. The quantity is distinguished from the gravitational potential realized in the System Frame.

Mbar(<R) — is the total baryonic mass enclosed within galactocentric radius R. It is determined from the independently specified baryonic source distribution and supplies the Newtonian gravitational availability entering the MQ realization law.

ρb — is the baryonic mass density of the independently specified physical source distribution. In the weak-field construction it determines the baryonic potential through the Newtonian Poisson relation and remains distinct from quantities derived after bounded realization.

SWFMQ
(S210)

— is the MQ weak-field realization action whose Euler-Lagrange equations produce both the Newtonian baryonic source equation and the MQ bounded gravitational realization equation. Its realization dependence is fixed by QMQ and therefore by the pre-existing MQ galactic realization law.

SMQdual
(S211)

— is the exact first-order dual form of the MQ spatial weak-field action. It replaces the primal nonlinear realization functional by the auxiliary field

YMQμ
(S212)

and the dual function RMQ. Eliminating the auxiliary field recovers the primal MQ Action exactly.

SK — is the timelike constraint action generated by KMQ(XMQ). Variation of SK with respect to its source-density multiplier imposes the physical mass-shell condition, while variation with respect to the action phase produces the conserved timelike source current.

TμνK
(S213)

— is the stress-energy tensor generated by metric variation of the timelike constraint action SK. On the physical mass-shell branch it assumes the pressureless timelike form determined by the source-density multiplier and the four-velocity uμ.

A.11 Statistical and observational symbols

rP — denotes the Pearson correlation coefficient used to quantify linear correlation between paired observational or calculated quantities. It is dimensionless and ranges from −1 to +1.

QSPARC — denotes the observational rotation curve quality classification assigned in the SPARC catalog [4], with QSPARC = 1, 2, and 3 corresponding respectively to high-, medium-, and low-quality classifications. It is independent of MQ quantities denoted by Q with MQ specific subscripts.

APPENDIX B: SOLUTIONS TO THE ELECTROMAGNETIC DEMARCATION

Blackbody demarcation

θsi=QLrc3G=(c3G)QLnLlf=c3lf2G.
(B.1)
QL=(1+nL2)1/2-nL.
(B.2)
θsi=(c3G)QLnLlf=(ℏlf2)QLnLlf.
(B.3)
nL=θsilfℏQL=θsilfℏ((1+nL2)1/2-nL).
(B.4)
(nL2+nL4)1/2-nL2=θsilfℏ.
(B.5)
nL2+nL4=(θsilfℏ+nL2)2=θsi2lf2ℏ2+2θsilfnL2ℏ+nL4.
(B.6)
nL2(1-2θsilfℏ)=θsi2lf2ℏ2.
(B.7)
nL=θsi2lf2ℏ2(1-(2θsilf/ℏ))=θsilf1ℏ(ℏ-2θsilf).
(B.8)
nL=84.6005394.
(B.9)

Charge coupling demarcation

nL=276/θsi=84.6005457.
(B.10)

Elementary Charge demarcation

γ=2(θsi-1)2(QLnL)4(84θsi-⌈42θsi⌋)rq412.
(B.11)
rq41=(θsi-1)(QLnL)2(284θsi-⌈42θsi⌋γ)1/2.
(B.12)
ec=2QLnLep=2QLnLrq41ef.
(B.13)
ec=2QLnL(θsi-1)(QLnL)2(2(84θsi-⌈42θsi⌋)γ)1/214(θsi-1)lfc2.
(B.14)
ec=lfc22QLnL(2(84θsi-⌈42θsi⌋)γ)1/2.
(B.15)
QLnL=lfc22ec(2(84θsi-⌈42θsi⌋)γ)1/2.
(B.16)
nL((1+nL2)1/2-nL)=lfc22ec(2(84θsi-⌈42θsi⌋)γ)1/2.
(B.17)
nL=84.6005398.
(B.18)

APPENDIX C: GEOMETRIC DOMAINS OF EXPANDING FRAMES

Ωtot=Ωobs+Ωdk
(C.1)
Ωdk+Ωobs=1
(C.2)
Ωobs=2θsiΩvis
(C.3)
Ωuobs=Ωobs-Ωvis
(C.4)
2θsi2=4ΩvisθsiΩobs=Ωf
(C.5)
Ωf=2θsi2=18.79135772(58)%
(C.6)
2θsiΩvis=Ωobs
(C.7)
2θsiΩvisθsi22-1=Ωobsθsi22-1
(C.8)
θsi3Ωvis-1=Ωobsθsi22-1
(C.9)
Ωdk=(Ωtot-Ωdk)θsi22-1=(1-Ωdk)θsi22-1
(C.10)
Ωdk=(1-Ωdk)θsi22-1
(C.11)
Ωdk=θsi2-2θsi2+2=68.36241612(59)%
(C.12)
1-Ωobs=Ωobsθsi22-1
(C.13)
2=Ωobsθsi22+Ωobs
(C.14)
Ωobs=4θsi2+2=31.63758388(59)%
(C.15)
Ωobs=2θsiΩvis=4θsi2+2
(C.16)
Ωvis=2θsi(θsi2+2)=4.84883489(10)%
(C.17)
Ωuobs=Ωobs-Ωvis=26.78874899(59)%
(C.18)

APPENDIX D: LOCAL SATURATION LAW FOR GRAVITATIONAL REALIZATION

Appendix D derives the local predictive galactic realization law from the MQ quantities used in the rotation-curve calculation. The derivation has two steps. First, the Newtonian circular-kinematic relation is evaluated with the positive quadratic System Frame norm established independently by the Frames mapping, fixing the dimensionless realization depth. Second, subdivision consistency requires the unrealized residual fraction to compose multiplicatively when realization depths add. Continuity with the zero-depth boundary condition selects an exponential semigroup, while identification of the normalized realization depth with its additive semigroup parameter fixes the generator used below. The acceleration normalization is determined independently by the MQ cosmological construction. The observed radial-acceleration relation and baryonic Tully-Fisher relation enter only as downstream comparisons.

Let R denote galactocentric radius, vNew(R) the Newtonian circular velocity generated by the adopted baryonic source, and vreal(R) the corresponding MQ realized circular velocity. Let gbar(R) denote the Newtonian baryonic acceleration, greal(R) the realized gravitational acceleration, and P(R) the realization fraction at radius R. The Newtonian baryonic acceleration is

gbar(R)=vNew2(R)R,
(D.1)

while the realized gravitational acceleration is

greal(R)=vreal2(R)R.
(D.2)

By definition of the realization fraction,

gbar(R)=P(R)greal(R),
(D.3)
greal(R)=gbar(R)P(R).
(D.4)

As P→1, the realized gravitational response approaches the Newtonian baryonic acceleration.

The MQ acceleration scale follows from the frame-resolved cosmological construction of Sec. 3.12. Let ΔH denote the expansion-specific metric differential between the frame-resolved scalar arguments, H the baseline dimensionless scalar expansion, Hsys and Hint the corresponding System Frame and Internal Frame local-rate realizations, and f the fixed linear mapping from a dimensionless scalar argument to its observable local-rate representation. The invariant coefficient is θsi, and π has its usual mathematical meaning. The relations are restated here only to identify the origin of Hint and the scalar contribution entering the acceleration normalization.

ΔH=θsi-π.
(D.5)

With the baseline scalar expansion H = θsi,

Hsys=f(H-ΔH),
(D.6)
Hint=f(H+ΔH).
(D.7)

Substituting the metric differential gives

Hsys=f(H-ΔH)=f(θsi-(θsi-π))=f(π).
(D.8)
Hint-f(ΔH)=Hsys+f(ΔH),
(D.9)
Hint=Hsys+f(2ΔH)=f(π)+f(2θsi-2π)=f(2θsi-π).
(D.10)

Thus, the doubled scalar contribution entering the Internal Frame local-rate realization Hint is fixed by the Frames mapping and metric differential, rather than inferred from galactic dynamics.

Let c denote the defined speed of light and aMQ the MQ galactic acceleration scale. When expressed in System Frame units, Hint carries inverse-time units, so cHint carries acceleration units. Consistent with Sec. 4.2 and Eq. (133), the MQ acceleration scale is

aMQ=cHint2θsi-1.
(D.11)

At a specified cosmological epoch, Eq. (D.11) is fixed before evaluation of the rotation curves and is applied unchanged to every galaxy at that epoch.

Equation (131) relates the local baryonic acceleration to the positive Newtonian circular-speed magnitude through the standard circular relation between acceleration, radius, and squared speed. The remaining normalization is fixed by the positive quadratic System Frame norm. Ref. [1], Proposition 2.4 and Apps. AI and N, establish the regular direction-free local realization and its realized tangent quadratic form. Applying the same circular-kinematic map to the independently fixed MQ acceleration scale at the same radius gives a positive reference speed magnitude. The ratio of the baryonic speed norm to this reference norm is dimensionless, and the common radius cancels. The realization depth is therefore

χ(R)=gbar(R)aMQ.
(D.12)

Equation (D.12) is a derived ratio of positive quadratic System Frame norms, not a homogeneous one-ratio ansatz. The one-half power follows because the local metric norm is the positive square root of its quadratic form and circular acceleration is proportional to squared speed at fixed radius. Replacing the square root by another power would change the established realization norm rather than introduce another admissible interpolation exponent. This construction is structurally parallel to deep-MOND and scale-invariant analyses that combine a fixed acceleration scale with circular kinematics to obtain the square-root weak-acceleration relation [19,21,52]. MQ differs in that its acceleration scale and quadratic norm are supplied by independent upstream constructions rather than fitted to galactic dynamics.

The exterior inverse-square baryonic field and the observed flat-curve or baryonic Tully-Fisher behavior are not used to choose the exponent. Eqs. (D.23–D.26) instead test the downstream consequences of the realization depth already fixed by Eq. (D.12).

With the realization depth fixed, the remaining step is the composition law for incomplete realization. The Frames mapping requires the result to be independent of how a fixed local depth is subdivided into consecutive admissible increments. The unrealized residual fraction is 1-P(χ). Subdivision consistency requires these residual fractions to multiply when realization depths add, while P(0)=0, positivity, and continuity apply throughout the admissible local domain. These conditions give a continuous exponential semigroup. Equation (D.12) supplies the normalized additive depth χ. MQ identifies that depth with the semigroup parameter itself, so its infinitesimal generator is unity rather than an additional fitted rate coefficient. The differential form is

dP=(1-P)dχ,
(D.13)
dPdχ=1-P.
(D.14)

With the boundary condition P(0)=0,

dP1-P=dχ,
(D.15)
-ln⁡(1-P)=χ,
(D.16)
P(χ)=1-e-χ.
(D.17)

Substituting Eq. (D.12),

P(R)=1-exp⁡[-gbar(R)aMQ].
(D.18)

The realized acceleration then follows from Eq. (D.4),

greal(R)=gbar(R)1-exp⁡[-gbar(R)/aMQ].
(D.19)

Using Eqs. (D.1) and (D.2),

greal(R)=vreal2(R)/R
(S214)

and

gbar(R)=vNew2(R)/R
(S215)

. Substitution into Eq. (D.19) gives

vreal2(R)R=vNew2(R)/R1-exp⁡[-gbar(R)/aMQ].
(D.20)

The common radial factor

1/R
(S216)

cancels, leaving

vreal2(R)=vNew2(R)1-exp⁡[-gbar(R)/aMQ],
(D.21)

and selection of the positive orbital velocity branch gives

vreal(R)=vNew(R)1-exp⁡[-gbar(R)/aMQ].
(D.22)

Equation (D.22) is identical in functional form to the empirical radial acceleration relation reported by McGaugh, Lelli, and Schombert [8]. That empirical relation is therefore prior phenomenology and is not claimed here as a new functional discovery. The MQ contribution is instead the derivation of the same form from independently established MQ quantities: the acceleration scale is fixed upstream of the galaxy data, the square-root realization depth follows from the quadratic System Frame norm and circular kinematics, and the exponential realization fraction follows from subdivision-consistent composition of the unrealized residual. None of these quantities are fitted to the rotation-curve data.

The limiting behavior follows from

χ=gbar/aMQ
(S217)

. For gbar≫aMQ, χ≫1, so e-χ→0, P→1, greal→gbar, and vreal→vNew, recovering the Newtonian baryonic limit.

For gbar≪aMQ, χ→0. The expansion

e-χ=1-χ+χ2/2-χ3/6+O(χ4)
(S218)

gives

1-e-χ=χ-χ2/2+χ3/6+O(χ4)
(S219)

and hence

limχ→01-e-χχ=1.
(D.23)

Equation (D.19) may therefore be evaluated in the weak-acceleration limit by dividing its denominator by its leading factor χ. Since

χ=gbar/aMQ
(S220)

,

limgbar/aMQ→0grealgbaraMQ=1.
(D.24)

Thus, the weak-acceleration realization has the asymptotic scaling

greal∼gbaraMQ
(S221)

. The realized acceleration still approaches zero as the baryonic acceleration approaches zero. The enhancement is relative to the more rapidly decreasing Newtonian baryonic field.

Thus, the realized gravitational response still approaches zero with the baryonic acceleration, but more slowly than the Newtonian baryonic field.

For a finite localized baryonic distribution with total baryonic mass Mbar, the exterior Newtonian baryonic acceleration approaches

gbar(R)=GMbar/R2
(S222)

, where G is the gravitational constant evaluated at the upper count limit. As R→∞,

gbar/aMQ→0
(S223)

; Eq. (D.24) therefore gives

limR→∞Rgreal(R)=GMbaraMQ.
(D.25)

For circular motion,

vreal2(R)=Rgreal(R)
(S224)

. Equation (D.25) therefore implies

limR→∞vreal2(R)=GMbaraMQ
(S225)

. Squaring this positive quantity gives

limR→∞vreal4(R)=GMbaraMQ.
(D.26)

The weak-acceleration limit therefore gives an asymptotic

1/R
(S226)

realized gravitational response rather than the Newtonian

1/R2
(S227)

decrease, while Eq. (D.26) yields both a radius-independent asymptotic orbital velocity and the baryonic Tully-Fisher fourth-power scaling. The same deep-acceleration scaling is a standard MOND result [19], while the baryonic Tully-Fisher relation is independently established observationally [20].

The exponential realization fraction can also be recovered from repeated finite realization increments. Let Pn denote the realization fraction after n equivalent normalized increments, so 1-Pn is the gravitational availability remaining unrealized. For total depth χ divided into N equal increments,

Δχ=χ/N
(S228)

. To first order in each vanishing increment, 1-Pn+1=(1-Δχ)(1-Pn).

Beginning from P0=0, repeated application gives

1-PN=(1-χ/N)N
(S229)

. At fixed total χ,

limN→∞(1-χN)N=e-χ.
(D.27)

It follows that P(χ)=1-e-χ, recovering Eq. (D.17). Differentiating this expression gives

dPdχ=1-P,
(D.28)

reproduces Eq. (D.14). The finite-increment construction therefore converges to the same subdivision-consistent exponential realization fraction. The exponential is not an independently introduced empirical interpolation function.

Accordingly, the acceleration normalization is inherited from the MQ expansion relation, the realization depth follows from the quadratic System Frame norm combined with circular kinematics, and the exponential realization fraction follows from subdivision-consistent composition. The weak- and high-acceleration limits are consequences of the same fixed galactic realization law.

APPENDIX E: EVIDENTIARY CLASSIFICATION OF QUANTITATIVE RESULTS

This appendix prevents algebraically related outputs, shared-input comparisons, and repeated analyses of the same observational information from being counted as separate empirical confirmations. For this manuscript, an input-derived result is fixed directly by an adopted measured or defined source and therefore is not an empirical test of that source. An internal closure result propagates an established MQ quantity back to a source observable, compares successive realization stages, or compares quantities with material shared-input dependence. It tests consistency but not independent prediction. An independent prediction uses a target observable that did not enter the derivation or calibration of the MQ quantity being tested. Independence is stated only with respect to that target observable: two tests may still share upstream MQ inputs and therefore are not automatically statistically independent. A robustness or source-sensitivity test changes an observational or source-side assumption while reusing the same prediction and much of the same data. It may constrain systematic dependence but does not constitute an additional independent confirmation. A numerical correspondence benchmark compares quantities not established as the same physical observable and therefore is not counted as empirical confirmation.

These distinctions also separate prospective from retrospective evidence. A prediction fixed before inspection of the target data carries stronger evidentiary status than a correspondence identified after the relevant observations are known, even when the latter remains mathematically informative [94].

“Independent prediction” below therefore identifies target independence; where a comparison is retrospective, model-dependent, incompletely blinded, or not yet executed prospectively, that limitation is stated explicitly.

Comparison E.1 applies these categories to the principal numerical claims used in the present paper. The classification concerns evidentiary status, not mathematical importance. A result may be central to the derivation while carrying little independent empirical weight, whereas a comparatively simple out-of-sample comparison may carry stronger evidentiary weight.

Comparison E.1. Evidentiary classification and dependency audit for principal quantitative results.

θsi from Eq. (76), Sec. 3.8

Principal upstream source: 2022 CODATA fine structure constant; MQ electromagnetic demarcation relation.

Classification: Input-derived.

Evidentiary interpretation: Parameter determination. The fine structure constant used to determine θsi cannot also validate θsi.

Fine structure realization, Eqs. (59–68) and Table 4

Principal upstream source: θsi from Eq. (76); same fine structure constant measurement family.

Classification: Internal closure.

Evidentiary interpretation: Tests algebraic closure of the realization sequence. The Morel result contributes to the measurement family entering the CODATA adjustment used upstream, so the agreement is not independent confirmation.

Electromagnetic demarcation routes, Table 3

Principal upstream source: Common θsi and MQ count geometry; one historical route also uses a measured electromagnetic source.

Classification: Internal closure.

Evidentiary interpretation: Distinct derivational routes converge on the same MQ quantity, but they are not independent laboratory determinations. Their agreement tests consistency of the shared count construction.

Historical finite-count G, upper-count G, and Planck-coupled G concordance family, Sec. 3.6

Principal upstream source: MQ fundamental measures and realization relations established independently of the compared historical gravitational constant recommendations.

Classification: Independent target comparison; retrospective concordance family.

Evidentiary interpretation: The compared gravitational constant values are not used to determine θsi or the MQ realization relation. The 2010 finite-count, 2014 upper-count, and 2018/2022 Planck-coupled numerical correspondences are nevertheless treated as one historical concordance family rather than three independent confirmations because they concern successive evaluations of the same MQ gravitational construction and a correlated CODATA recommendation sequence.

2026 NIST/BIPM G comparison, Sec. 3.6

Principal upstream source: Pre-existing MQ finite-count and upper-count gravitational predictions; independent 2026 torsion balance measurement [39].

Classification: Independent target comparison.

Evidentiary interpretation: The macroscopic gravitational measurement is temporally out of sample relative to the MQ gravitational constant values and is therefore evidentially distinct from the historical CODATA concordance family. However, the assignment of free-deflection and electrostatic-servo operation to different MQ realization regimes was not preregistered against the 2026 result. The observed ordering is consequently a retrospective compatibility test, not a prospective causal confirmation of the proposed realization-mode distinction.

Referenceability-domain fractions, Eqs. (77–81) and Table 5

Principal upstream source: θsi determined from the fine structure constant.

Classification: Numerical correspondence benchmark; not counted as an empirical prediction.

Evidentiary interpretation: Planck density fractions are not used to derive the MQ referenceability-domain percentages, but the two sets of quantities have different physical definitions. Their numerical proximity therefore does not constitute prediction of the ΛCDM density fractions. It is a cross-framework correspondence whose physical significance remains contingent on the independent dynamical cosmological calculation required in Sec. 5.8.

CMB temperature, Eqs. (115–118) and Table 6

Principal upstream source: MQ count geometry together with an adopted Planck-derived present cosmic age.

Classification: Internal cross-scale closure.

Evidentiary interpretation: The measured CMB monopole temperature is not supplied to the temperature calculation, but the adopted present universe age is inferred largely from CMB anisotropy data. The numerical temperature agreement therefore tests cross-scale closure but is not statistically independent of the broader CMB information used upstream. The radiation-density comparison is likewise not an additional test because it is derived from the same CMB temperature.

Frame-resolved Hubble values, Eqs. (127–128), Sec. 3.12

Principal upstream source: θsi, MQ expansion geometry, mapping normalization, and adopted cosmic-age normalization.

Classification: Internal closure.

Evidentiary interpretation: The comparison provides a quantitative benchmark of the published frame resolution construction. The two frame-resolved rates share the same MQ scalar, age normalization, and mapping structure, while the frame interpretation was formulated retrospectively relative to the observational Hubble-tension measurements. Their agreements are therefore not counted as two independent empirical confirmations.

Present-epoch aMQ, Eq.
(133)

Principal upstream source: Hint, θsi, and defined c.

Classification: Input-derived.

Evidentiary interpretation: aMQ is an upstream derived quantity fixed before application to the galaxy data. Its numerical value is not itself an observational validation. The relevant empirical test begins only when the frozen value is applied to an independently reconstructed baryonic source and compared with an external kinematic target.

Empirical RAR functional form, Sec. 4.2

Principal upstream source: Published empirical RAR functional form [8].

Classification: Background; not counted as evidence for MQ.

Evidentiary interpretation: The functional identity is acknowledged as prior phenomenology. Fitted RAR acceleration scales are model-conditioned regression parameters and are not treated as measurements of aMQ; no numerical concordance or tension is assigned between them. The MQ test concerns the independently derived normalization and physical provenance of the relation, not discovery of its functional shape.

Milky Way rotation curve, Sec. 4.3

Principal upstream source: McMillan baryonic source [3]; fixed MQ galactic realization law; independent Eilers kinematics [2].

Classification: Independent prediction.

Evidentiary interpretation: This is a RAR/MOND-independent source/target test because the McMillan baryonic model and Eilers circular velocity measurements are distinct published constructions and no Milky Way-specific MQ parameter is fitted. Application of the fixed galactic realization law reduces the Newtonian baryonic RMS discrepancy from 69.03 km s-1 to 8.90 km s-1, an 87.11% reduction. The diagonal Eilers/Jiao error treatment formally rejects the complete sample, but that conclusion is not covariance-robust. Propagating the specified McMillan conditional baryonic-normalization covariance gives generalized χ2 = 39.57 for 38 degrees of freedom with p = 0.400; the inner and outer subsets likewise give p = 0.423 and p = 0.910, respectively. The result is therefore non-rejected under that source-covariance model, but it is not a covariance-complete precision validation because the full Eilers point-to-point covariance and complete McMillan source posterior have not been jointly propagated.

Complete SPARC catalog and high-quality subset, Secs. 4.4–4.10

Principal upstream source: Published SPARC baryonic reconstructions [4]; fixed aMQ and galactic realization law.

Classification: Independent prediction of the fixed MQ normalization, conditional on the adopted SPARC source construction.

Evidentiary interpretation: Observed velocities do not determine aMQ or another galaxy-specific MQ parameter. Across all 175 galaxies and 3,366 radial measurements, the calculation tests the frozen MQ normalization against the published catalog. It is not an independent discovery test of the pre-existing RAR functional form. The complete population retains a significant velocity-dependent normalization residual, whereas the 99 independently catalog-classified high-quality galaxies have a geometric MQ/observed velocity ratio of approximately 1.0042 and no statistically significant marginal velocity dependence. The quality dependence constrains interpretation of the complete-catalog trend but does not demonstrate that lower quality residuals are observational errors. The result remains conditional on the published baryonic reconstruction, distance, inclination, stellar normalization, and correlated galaxy-property structure.

LITTLE THINGS independent population, Secs. 4.10 and 5.9

Principal upstream source: Oh et al. dwarf-galaxy source information [5]; O’Brien et al. published rotation curve figures [6]; fixed MQ galactic realization law.

Classification: Independent population prediction.

Evidentiary interpretation: This is an out-of-SPARC population-normalization and residual-ordering test using 15 prospectively screened non-SPARC galaxies and 288 radial comparisons. The galaxy-weighted geometric MQ/observed ratio is 0.9907 with a 95% galaxy-block-bootstrap interval of 0.868–1.130 and p = 0.895 for zero mean galaxy logarithmic offset. The primary logarithmic-offset correlation with characteristic velocity is not significant (r_P = -0.282, p = 0.308), and the partial Pearson correlation controlling SED stellar mass and inclination is likewise nonsignificant (r_P = -0.257, p = 0.355). These results provide an independent population-level normalization test, but figure digitization and reconstructed-source uncertainties preclude treating the calculation as a precision pointwise likelihood.

WISE stellar mass and gas-rich controls, Secs. 4.6–4.8

Principal upstream source: Same fixed MQ galactic realization law applied after alternative Duey et al. stellar normalization or to residual-blind source-stable subsamples.

Classification: Robustness and source-sensitivity tests.

Evidentiary interpretation: These analyses test whether SPARC residual structure is sensitive to the adopted stellar source and whether gas-rich systems behave differently when stellar normalization is dynamically subdominant. They reuse the same galactic law and substantially overlapping kinematic information and therefore are not separate confirmations of MQ. The WISE substitution is not an independent kinematic population: it replaces the stellar normalization within overlapping SPARC galaxies. Because Duey et al. provide galaxy-integrated rather than fully resolved radial stellar masses, residual stellar-profile uncertainty remains in stellar-important systems.

Cosmic-age/redshift evolution of aMQ, Eqs. (192–195), Sec. 5.5

Principal upstream source: Present MQ acceleration-scale construction, MQ age relation, and externally specified redshift-to-age expansion history.

Classification: Independent prospective prediction.

Evidentiary interpretation: Once the cosmological age relation and expansion history are fixed, MQ specifies the redshift dependence of aMQ without galaxy-specific refitting. Equations (192) and (191) provide the conditional redshift-to-age conversion, while Eqs. (194) and (193) specify the corresponding fixed-source galactic prediction. Equations (190) and (189) establish the present age normalization used upstream but are not themselves part of the redshift-dependent prediction. The preferred intermediate-redshift test has not yet been performed with independently reconstructed baryonic sources and frozen MQ parameters. Its evidentiary value therefore remains prospective. RC100 and KROSS establish the observational feasibility of resolved kinematics in the relevant redshift range [r,s], but their existing dynamical decompositions are not themselves MQ tests.

Lensing, Bullet Cluster, and Einstein–Boltzmann calculations, Secs. 5.3–5.9

Principal upstream source: Derived MQ metric, realization, source, and action-level dynamics.

Classification: Independent prospective tests.

Evidentiary interpretation: These regimes are physically independent requirements rather than confirmations inherited from rotation curve agreement. The spherical ESO 325-G004 lensing calculation is presently a weak-field benchmark, not a complete extended-source prediction; the Bullet Cluster still requires nonlinear baryonic merger evolution followed by lensing ray tracing; and Sec. 5.8 establishes mathematical/action-level closure sufficient to pose a numerical Einstein–Boltzmann calculation but does not provide CMB or large-scale-structure spectra. None therefore contributes a completed empirical confirmation at present.

Comparison E.1 fixes how cumulative evidence is counted throughout the manuscript. Algebraic closure, repeated evaluations of one derived constant, numerical correspondence between physically different quantities, or reanalysis of the same galaxy kinematics under alternative source assumptions may strengthen internal consistency or robustness, but they do not multiply the number of independent empirical tests. Conversely, target independence alone does not imply statistical independence when tests share upstream MQ quantities.

The Milky Way provides a RAR/MOND-independent source/target test, SPARC tests the same frozen galactic realization law across its published population, and LITTLE THINGS supplies the completed out-of-SPARC population test. Neither the Milky Way nor SPARC constitutes an independent discovery of the empirical RAR functional form, and the WISE and related source substitutions constrain source model sensitivity rather than add new observational populations.

The redshift dependence of the galactic realization law remains a prospective independent discriminator. Section 5.5 specifies the intermediate-redshift test architecture before an MQ residual analysis is performed: the baryonic reconstruction, expansion history conversion, and MQ parameters are to be fixed before inspection of the target residuals. This separation between model construction and evaluation on new observations is the evidentiary distinction intended here [94].

The real-lens, cluster merger, and Einstein–Boltzmann calculations likewise remain prospective rather than completed confirmations. In stellar-important galaxies, precision evidentiary weight is also withheld until a resolved radial stellar surface-density profile is propagated directly into the baryonic source; the five-galaxy residual-blind LITTLE THINGS profile test identified in Sec. 4.10 therefore remains prospective rather than an additional completed confirmation.

STATEMENTS AND DECLARATIONS

The author has no relevant financial or non–financial interests to disclose.

The author has no competing interests to declare that are relevant to the content of this article.

DATA AVAILABILITY

The data that support the findings of this study are derived from publicly available observational and baryonic-source products cited in the article. The Milky Way analysis uses the published Eilers et al. circular-velocity measurements [2] and the McMillan Galactic mass model [3]; the complete McMillan MCMC chain used for posterior source-uncertainty propagation is publicly archived as described in Ref. [95]. The SPARC analysis uses the publicly available SPARC rotation-curve and baryonic-component data [4]. The LITTLE THINGS analysis uses the published galaxy, rotation-curve, H I and stellar-source products of Oh et al. [5] and O’Brien et al. [6], including the associated CDS/VizieR machine-readable products. The resolved stellar-source information discussed as a prospective extension is available from Zhang et al. [57] and Oh et al. [5]. No new observational data were generated in this study. Data extracted or calculated from these published sources that are required to reproduce the numerical results reported in this article are available from the author upon reasonable request.

FUNDING

This research has received no external funding.

REFERENCES

[1] Geiger J A 2025 Measurement quantization: A count-based resolution of the quantum measurement problem preprint, Informativity Institute, ResearchGate

[2] Eilers A-C, Hogg D W, Rix H-W and Ness M K 2019 The circular velocity curve of the Milky Way from 5 to 25 kpc Astrophys. J. 871 120

[3] McMillan P J 2017 The mass distribution and gravitational potential of the Milky Way Mon. Not. R. Astron. Soc. 465 76–94 (doi: 10.1093/mnras/stw2759)

[4] Lelli F, McGaugh S S and Schombert J M 2016 SPARC: Mass models for 175 disk galaxies with Spitzer photometry and accurate rotation curves Astron. J. 152 157

[5] Oh S-H et al 2015 High-resolution mass models of dwarf galaxies from LITTLE THINGS Astron. J. 149 180 (doi: 10.1088/0004-6256/149/6/180)

[6] O’Brien J G, Chiarelli T L, Dentico J, Stulge M, Stefanski B, Moss R and Chaykov S 2018 Alternative gravity rotation curves for the LITTLE THINGS survey Astrophys. J. 852 6

[7] Mohr P J, Newell D B, Taylor B N and Tiesinga E 2025 CODATA recommended values of the fundamental physical constants: 2022 Rev. Mod. Phys. 97 025002

[8] McGaugh S S, Lelli F and Schombert J M 2016 Radial acceleration relation in rotationally supported galaxies Phys. Rev. Lett. 117 201101 (doi: 10.1103/PhysRevLett.117.201101)

[9] Hunter D A et al 2012 LITTLE THINGS Astron. J. 144 134

[10] Garrett K and Duda G 2011 Dark matter: A primer Adv. Astron. 2011 968283

[11] Aghanim N et al (Planck Collaboration) 2020 Planck 2018 results. VI. Cosmological parameters Astron. Astrophys. 641 A6 (doi: 10.1051/0004-6361/201833910)

[12] Bertone G and Hooper D 2018 History of dark matter Rev. Mod. Phys. 90 045002

[13] Aprile E et al (XENON Collaboration) 2022 Search for new physics in electronic recoil data from XENONnT Phys. Rev. Lett. 129 161805

[14] Du N et al (ADMX Collaboration) 2018 A search for invisible axion dark matter with the Axion Dark Matter Experiment Phys. Rev. Lett. 120 151301

[15] Geiger J A 2023 Measurement quantization Int. J. Geom. Methods Mod. Phys. 20 2350069

[16] Rubin V C, Ford W K Jr and Thonnard N 1980 Rotational properties of 21 Sc galaxies with a large range of luminosities and radii, from NGC 4605 ( kpc) to UGC 2885 ( kpc) Astrophys. J. 238 471

[17] Clowe D, Bradač M, Gonzalez A H, Markevitch M, Randall S W, Jones C and Zaritsky D 2006 A direct empirical proof of the existence of dark matter Astrophys. J. Lett. 648 L109–L113 (doi: 10.1086/508162)

[18] Geiger J A 2020 Measurement quantization describes history of universe—Quantum inflation, transition to expansion, CMB power spectrum J. High Energy Phys. Gravit. Cosmol. 6 186

[19] Milgrom M 1983 A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis Astrophys. J. 270 365–370 (doi: 10.1086/161130)

[20] McGaugh S S, Schombert J M, Bothun G D and de Blok W J G 2000 The baryonic Tully-Fisher relation Astrophys. J. Lett. 533 L99

[21] Famaey B and McGaugh S S 2012 Modified Newtonian dynamics (MOND): Observational phenomenology and relativistic extensions Living Rev. Relativ. 15 10 (doi: 10.12942/lrr-2012-10)

[22] Navas S et al (Particle Data Group) 2024 Review of Particle Physics Phys. Rev. D 110 030001 (doi: 10.1103/PhysRevD.110.030001)

[23] Clowe D, Gonzalez A and Markevitch M 2004 Weak-lensing mass reconstruction of the interacting cluster 1E 0657−558: Direct evidence for the existence of dark matter Astrophys. J. 604 596

[24] Markevitch M, Gonzalez A H, Clowe D, Vikhlinin A, Forman W, Jones C, Murray S and Tucker W 2004 Direct constraints on the dark matter self-interaction cross section from the merging galaxy cluster 1E 0657−56 Astrophys. J. 606 819

[25] Navarro J F, Frenk C S and White S D M 1996 The structure of cold dark matter halos Astrophys. J. 462 563

[26] Klypin A A, Kravtsov A V, Valenzuela O and Prada F 1999 Where are the missing Galactic satellites? Astrophys. J. 522 82

[27] Boylan-Kolchin M, Bullock J S and Kaplinghat M 2011 Too big to fail? The puzzling darkness of massive Milky Way subhaloes Mon. Not. R. Astron. Soc. Lett. 415 L40

[28] Oman K A et al 2015 The unexpected diversity of dwarf galaxy rotation curves Mon. Not. R. Astron. Soc. 452 3650

[29] Bullock J S and Boylan-Kolchin M 2017 Small-scale challenges to the ΛCDM paradigm Annu. Rev. Astron. Astrophys. 55 343–387 (doi: 10.1146/annurev-astro-091916-055313)

[30] Sikivie P 2021 Invisible axion search methods Rev. Mod. Phys. 93 015004 (doi: 10.1103/RevModPhys.93.015004)

[31] Boyarsky A, Ruchayskiy O and Shaposhnikov M 2009 The role of sterile neutrinos in cosmology and astrophysics Annu. Rev. Nucl. Part. Sci. 59 191–214 (doi: 10.1146/annurev.nucl.010909.083654)

[32] Adhikari S et al 2025 Astrophysical tests of dark matter self-interactions Rev. Mod. Phys. 97 045004 (doi: 10.1103/m2vm-59y3)

[33] Eberhardt A and Ferreira E G M 2026 Ultralight dark matter: Progress and current constraints Rev. Mod. Phys. (doi: 10.1103/j9ly-b6q3)

[34] Poisson E, Pound A and Vega I 2011 The motion of point particles in curved spacetime Living Rev. Relativ. 14 7 (doi: 10.12942/lrr-2011-7)

[35] Mohr P J, Taylor B N and Newell D B 2012 CODATA recommended values of the fundamental physical constants: 2010 Rev. Mod. Phys. 84 1527

[36] Mohr P J, Newell D B and Taylor B N 2016 CODATA recommended values of the fundamental physical constants: 2014 Rev. Mod. Phys. 88 035009

[37] Tiesinga E, Mohr P J, Newell D B and Taylor B N 2021 CODATA recommended values of the fundamental physical constants: 2018 Rev. Mod. Phys. 93 025010

[38] Rosi G, Sorrentino F, Cacciapuoti L, Prevedelli M and Tino G M 2014 Precision measurement of the Newtonian gravitational constant using cold atoms Nature 510 518

[39] Schlamminger S et al 2026 Redetermination of the gravitational constant with the BIPM torsion balance at NIST Metrologia 63 025012

[40] Morel L, Yao Z, Cladé P and Guellati-Khélifa S 2020 Determination of the fine structure constant with an accuracy of 81 parts per trillion Nature 588 61

[41] Fixsen D J 2009 The temperature of the cosmic microwave background Astrophys. J. 707 916

[42] Geiger J A 2025 Resolving the Hubble tension 68.261 and 73.510 km s Mpc AIP Adv. 15 075134

[43] Balkenhol L et al 2023 Measurement of the CMB temperature power spectrum and constraints on cosmology from the SPT-3G 2018 TT, TE, and EE dataset Phys. Rev. D 108 023510

[44] Brout D et al 2022 The Pantheon+ analysis: Cosmological constraints Astrophys. J. 938 110

[45] Oman K A, Marasco A, Navarro J F, Frenk C S, Schaye J and Benitez-Llambay A 2019 Non-circular motions and the diversity of dwarf galaxy rotation curves Mon. Not. R. Astron. Soc. 482 821–847

[46] McGaugh S S 2019 The imprint of spiral arms on the Galactic rotation curve Astrophys. J. 885 87

[47] Jiao Y, Hammer F, Wang J L and Yang Y B 2021 Which Milky Way masses are consistent with the slightly declining 5–25 kpc rotation curve? Astron. Astrophys. 654 A25

[48] McGaugh S S and Schombert J M 2014 Color-mass-to-light-ratio relations for disk galaxies Astron. J. 148 77

[49] Meidt S E et al 2014 Reconstructing the stellar mass distributions of galaxies using S4G IRAC 3.6 and 4.5 μm images. II. The conversion from light to mass Astrophys. J. 788 144

[50] Ponomareva A A, Verheijen M A W, Papastergis E, Bosma A and Peletier R F 2018 From light to baryonic mass: the effect of the stellar mass-to-light ratio on the baryonic Tully-Fisher relation Mon. Not. R. Astron. Soc. 474 4366–4384

[51] Duey F, Schombert J, McGaugh S S and Lelli F 2025 The baryonic Tully-Fisher relation. II. Stellar mass models Astron. J. 169 186

[52] Milgrom M 2014 MOND laws of galactic dynamics Mon. Not. R. Astron. Soc. 437 2531

[53] Oman K A and Riley A H 2024 An overlooked source of uncertainty in the mass of the Milky Way Mon. Not. R. Astron. Soc. Lett. 532 L48

[54] Field C A and Welsh A H 2007 Bootstrapping clustered data J. R. Stat. Soc. B 69 369–390

[55] Kruskal W H and Wallis W A 1952 Use of ranks in one-criterion variance analysis J. Am. Stat. Assoc. 47 583

[56] Martinsson T P K, Verheijen M A W, Westfall K B, Bershady M A, Andersen D R and Swaters R A 2013 The DiskMass Survey. VII. The distribution of luminous and dark matter in spiral galaxies Astron. Astrophys. 557 A131

[57] Zhang H-X, Hunter D A, Elmegreen B G, Gao Y, Schruba A and collaborators 2012 Outside-in shrinking of the star-forming disk of dwarf irregular galaxies Astron. J. 143 47 (doi: 10.1088/0004-6256/143/2/47)

[58] Lelli F, McGaugh S S, Schombert J M, Desmond H and Katz H 2019 The baryonic Tully-Fisher relation for different velocity definitions and implications for galaxy angular momentum Mon. Not. R. Astron. Soc. 484 3267

[59] Stiskalek R and Desmond H 2023 On the fundamentality of the radial acceleration relation for late-type galaxy dynamics Mon. Not. R. Astron. Soc. 525 6130

[60] Desmond H 2023 The underlying radial acceleration relation Mon. Not. R. Astron. Soc. 526 3342–3351

[61] Bland J M and Altman D G 1986 Statistical methods for assessing agreement between two methods of clinical measurement Lancet 1 307

[62] Bartelmann M and Schneider P 2001 Weak gravitational lensing Phys. Rep. 340 291–472

[63] Smith R J and Lucey J R 2013 A giant elliptical galaxy with a lightweight initial mass function Mon. Not. R. Astron. Soc. 434 1964–1977

[64] Collett T E et al 2018 A precise extragalactic test of General Relativity Science 360 1342–1346

[65] Jain B and Zhang P 2008 Observational tests of modified gravity Phys. Rev. D 78 063503

[66] Daniel S F, Linder E V, Smith T L, Caldwell R R, Cooray A, Leauthaud A and Lombriser L 2010 Testing general relativity with current cosmological data Phys. Rev. D 81 123508

[67] Blanchet L 2002 Gravitational radiation from post-Newtonian sources and inspiralling compact binaries Living Rev. Relativ. 5 3 (doi: 10.12942/lrr-2002-3)

[68] Bradač M et al 2006 Strong and weak lensing united III: Measuring the mass distribution of the merging galaxy cluster 1E 0657−56 Astrophys. J. 652 937–947

[69] Mastropietro C and Burkert A 2008 Simulating the Bullet Cluster Mon. Not. R. Astron. Soc. 389 967–988

[70] Springel V and Farrar G R 2007 The speed of the ‘bullet’ in the merging galaxy cluster 1E 0657−56 Mon. Not. R. Astron. Soc. 380 911–925

[71] Geiger J A 2026 General relativity as an emergent geometry of discrete measurement manuscript

[72] Bekenstein J D and Milgrom M 1984 Does the missing mass problem signal the breakdown of Newtonian gravity? Astrophys. J. 286 7–14

[73] Milgrom M 2010 Quasi-linear formulation of MOND Mon. Not. R. Astron. Soc. 403 886–895

[74] Nestor Shachar A et al 2023 RC100: Rotation curves of 100 massive star-forming galaxies at z = 0.6–2.5 reveal little dark matter on galactic scales Astrophys. J. 944 78 (doi: 10.3847/1538-4357/aca9cf)

[75] Harrison C M et al 2017 The KMOS Redshift One Spectroscopic Survey (KROSS): Rotational velocities and angular momentum of z ≈ 0.9 galaxies Mon. Not. R. Astron. Soc. 467 1965–1983 (doi: 10.1093/mnras/stx217)

[76] de Graaff A et al 2024 Ionised gas kinematics and dynamical masses of galaxies from JADES/NIRSpec high-resolution spectroscopy Astron. Astrophys. 684 A87 (doi: 10.1051/0004-6361/202347755)

[77] Xu Y et al 2024 Dynamics of a galaxy at z > 10 explored by JWST integral field spectroscopy: Hints of a rotating disk suggesting weak feedback Astrophys. J. 976 142 (doi: 10.3847/1538-4357/ad82dd)

[78] Gogate A R, Verheijen M A W, van der Hulst J M and Jaffé Y L 2023 BUDHIES V: The baryonic Tully-Fisher relation based on direct H I detections Mon. Not. R. Astron. Soc. 519 4279

[79] Stott J P et al 2016 The KMOS Redshift One Spectroscopic Survey (KROSS): dynamical properties, gas and dark matter fractions of typical z ~ 1 star-forming galaxies Mon. Not. R. Astron. Soc. 457 1888–1904

[80] Mukhanov V F, Feldman H A and Brandenberger R H 1992 Theory of cosmological perturbations Phys. Rep. 215 203–333

[81] Ma C-P and Bertschinger E 1995 Cosmological perturbation theory in the synchronous and conformal Newtonian gauges Astrophys. J. 455 7–25

[82] Skordis C 2006 Tensor-vector-scalar cosmology: Covariant formalism for the background evolution and linear perturbation theory Phys. Rev. D 74 103513

[83] Nesseris S 2023 The effective fluid approach for modified gravity and its applications Universe 9 13

[84] Candlish G N 2016 The velocity field in MOND cosmology Mon. Not. R. Astron. Soc. 460 2571–2585

[85] Geiger J A 2026 Measurement Quantization: A Count-Based Resolution of the Quantum Measurement Problem manuscript

[86] Skordis C and Zlosnik T 2021 New relativistic theory for modified Newtonian dynamics Phys. Rev. Lett. 127 161302

[87] Hwang J-c and Noh H 2025 Cosmological perturbations of a relativistic MOND theory Phys. Rev. D 111 063547

[88] Gourgoulhon É 2012 3+1 Formalism in General Relativity: Bases of Numerical Relativity Lecture Notes in Physics vol 846 (Springer, Berlin)

[89] Bekenstein J D 2004 Relativistic gravitation theory for the modified Newtonian dynamics paradigm Phys. Rev. D 70 083509 erratum Phys. Rev. D 71 069901 (2005)

[90] Zlosnik T G, Ferreira P G and Starkman G D 2007 Modifying gravity with the aether: An alternative to dark matter Phys. Rev. D 75 044017

[91] Santambrogio F 2018 Regularity via duality in calculus of variations and degenerate elliptic PDEs J. Math. Anal. Appl. 457 1649

[92] Blas D, Lesgourgues J and Tram T 2011 The Cosmic Linear Anisotropy Solving System (CLASS). II. Approximation schemes J. Cosmol. Astropart. Phys. 07 034

[93] Alam S et al (eBOSS Collaboration) 2021 The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological implications from two decades of spectroscopic surveys at the Apache Point Observatory Phys. Rev. D 103 083533

[94] Nosek B A, Ebersole C R, DeHaven A C and Mellor D T 2018 The preregistration revolution Proc. Natl Acad. Sci. USA 115 2600–2606 (doi: 10.1073/pnas.1708274114)

[95] McMillan P J 2022 MCMC chain of Milky Way gravitational potential models from McMillan (2017, MNRAS, 465, 76) Zenodo (doi: 10.5281/zenodo.6219842)

line.jpg