The Fundamental Speed Theory: A Mathematically Consistent Vector-Tensor Theory for Galactic Dynamics Without Dark Matter Updated Results from 171 SPARC Galaxies

1. Introduction

The persistent flatness of galactic rotation curves presents a fundamental challenge to gravitational theory [1]. While the $\Lambda$CDM paradigm successfully explains cosmological observations [2], direct detection of particle dark matter remains elusive [3]. Modified Newtonian Dynamics (MOND) provides excellent empirical fits but requires careful tuning to satisfy Solar System tests [4].

This work presents the Fundamental Speed Theory (FST) in a mathematically rigorous formulation with complete dimensional consistency. FST introduces a dimensionless vector field ${\nu }^{\mu }$ coupled to gravity through characteristic mass and length scales $M_0$ and $L_0$, with explicit inclusion of fundamental constants ℏ and c throughout. Key contributions:

• Dimensional rigor: Complete unit analysis with proper handling of ℏ and c
• Hierarchical validation: Three levels of testing show theory’s predictive power: Level 3 (zero parameters): ${\chi }^2=0.809$ (65.7%); Level 2 (estimated): ${\chi }^2=0.347$ (93.6%); Level 1 (fitted): ${\chi }^2=0.170$ (100%)
• Empirical success: Mean ${\chi }_{\nu }^2=0.170$ on 171 SPARC galaxies with universal parameters
• Parameter unification: All five field parameters unify into a single fundamental acceleration scale $A_0=2.42\times {10}^{-10}$ m/s²
• Characteristic scales: ${\xi }_c=3.16\times {10}^{-4}$, $r_c\approx 3.16$ pc, ${\lambda }_{\mathrm{screen}}=1.65$ pc
• Dynamical Families: Cluster analysis uncovers three distinct galaxy populations
• Coefficient independence: Kinetic coefficients are not critical ($\pm 44\%$ variation changes ${\chi }^2$ by $\pm 0.1\%$)
• Theoretical economy: Six dimensionless parameters for all galaxy types, unified into one fundamental scale
• Solar System compatibility: Natural explanation from galactic field gradient satisfies all local tests, with FST acceleration at Earth more than 100,000 times below current observational bounds
• Computational transparency: Open-source implementation with unit verification

FST demonstrates that vector-tensor gravity can explain galactic dynamics without dark matter while maintaining mathematical consistency across all scales from Solar System to galactic halos. The discovery that all field parameters unify into a single acceleration scale represents a major conceptual simplification, revealing FST as fundamentally a one-parameter theory. Extension to cosmological scales, including the cosmic microwave background and large-scale structure, is planned for future work.

2. Dimensional Framework and Fundamental Constants

2.1. System of Units and Constants

We maintain explicit awareness of both natural units ($\hslash =c=1$) and SI units for numerical calculations. When working in natural units, all equations are dimensionally consistent by construction, and the conversion back to SI units is achieved by reinserting the appropriate factors of ℏ and c where necessary. The fundamental constants are taken from the CODATA 2022 recommended values [5]:

$$\begin{array}{cc}\hfill \hslash & =1.054571817\times {10}^{-34}\mathrm{J}·\mathrm{s}\left(\mathrm{reduced}\mathrm{Planck}\mathrm{constant}\right)\hfill \end{array}$$

(1)

$$\begin{array}{cc}\hfill c& =2.99792458\times {10}^8\mathrm{m}·{\mathrm{s}}^{-1}\left(\mathrm{speed}\mathrm{of}\mathrm{light}\right)\hfill \end{array}$$

(2)

$$\begin{array}{cc}\hfill G_N& =6.67430\times {10}^{-11}{\mathrm{m}}^3·{\mathrm{kg}}^{-1}·{\mathrm{s}}^{-2}\left(\mathrm{Newton}\text{’}\mathrm{s}\mathrm{constant}\right)\hfill \end{array}$$

(3)

2.2. Characteristic Scales of FST

2.2.1. Characteristic Length Scale $L_0$

The characteristic length scale $L_0$ is chosen as the typical scale length of spiral galaxies, motivated by the SPARC database [1]:

$$L_0=10\mathrm{kpc}=3.086\times {10}^{20}\mathrm{m}$$

(1)

This is not a free parameter but the median scale length of the galaxy sample. The theory’s predictions are robust to variations in this choice, as demonstrated in Section 8.11.

2.2.2. Characteristic Mass Scale $M_0$

From the fundamental constants and $L_0$, the characteristic mass scale $M_0$ is derived:

$$M_0={\displaystyle \frac{\hslash }{c{L}_0}}={\displaystyle \frac{1.054571817\times {10}^{-34}}{2.99792458\times {10}^8\times 3.086\times {10}^{20}}}=1.140\times {10}^{-63}\mathrm{kg}$$

(2)

In energy units: $M_0{c}^2=1.024\times {10}^{-46}\mathrm{J}=6.403\times {10}^{-28}\mathrm{eV}$

This scale ensures dimensional consistency in the Lagrangian by providing a normalization for the dimensionless field ${\nu }^{\mu }$; it does not represent a physical particle mass.

2.2.3. Dimensionless Field Definition

The physical vector field $V^{\mu }$ relates to dimensionless ${\nu }^{\mu }$ by:

$$V^{\mu }=M_0{\nu }^{\mu },\left[{\nu }^{\mu }\right]=1,\left[V^{\mu }\right]=\left[M\right]=\mathrm{kg}$$

(3)

This ensures that all physical quantities have correct dimensions while keeping the fundamental field equations dimensionless.

2.3. Dimensional Analysis Table

Table 1. Dimensional Analysis of Key Quantities

Quantity	Symbol	SI Units
Length	L	m
Mass	M	kg
Time	T	s
Action	S	J·s
Lagrangian Density	$\mathcal{L}$	J/m³
Vector Field	$V^{\mu }$	kg
Dimensionless Field	${\nu }^{\mu }$	1
Characteristic Mass	$M_0$	kg
Characteristic Length	$L_0$	m
FST Acceleration	$a_{\mathrm{FST}}$	m/s²

Table 1. Dimensional Analysis of Key Quantities

Quantity	Symbol	SI Units
Length	L	m
Mass	M	kg
Time	T	s
Action	S	J·s
Lagrangian Density	$\mathcal{L}$	J/m³
Vector Field	$V^{\mu }$	kg
Dimensionless Field	${\nu }^{\mu }$	1
Characteristic Mass	$M_0$	kg
Characteristic Length	$L_0$	m
FST Acceleration	$a_{\mathrm{FST}}$	m/s²

3. Theoretical Framework

3.1. Action Principle

The total action in Jordan frame with explicit constants:

$$S=\int d^{4}x\sqrt{-g}\left[{\displaystyle \frac{c^4}{16\pi G_N}}R+{\mathcal{L}}_V+{\mathcal{L}}_m\right]$$

(4)

where the Einstein-Hilbert term has factor $c^4/\left(16\pi G_N\right)$ for correct dimensions:

$[c^4/G_N]=\left[ML{T}^{-2}\right]$, $\left[R\right]=\left[L^{-2}\right]$, so $[c^{4}R/G_N]=\left[M{L}^{-1}{T}^{-2}\right]=\left[\mathcal{L}\right]$

3.2. Dimensionally Consistent Lagrangian

The vector field Lagrangian density, constructed to ensure dimensional consistency with the action principle, is:

$${\mathcal{L}}_V={\displaystyle \frac{c^4}{16\pi G_N{L}_0^2}}\left[-{\displaystyle \frac{c_1}{2}}\left(L_0^2{\nabla }_{\mu }{\nu }_{\nu }\right)\left({\nabla }^{\mu }{\nu }^{\nu }\right)-{\displaystyle \frac{c_2}{2}}{L}_0^2{\left({\nabla }_{\mu }{\nu }^{\mu }\right)}^2-{\displaystyle \frac{c_3}{2}}\left(L_0^2{\nabla }_{\mu }{\nu }_{\nu }\right)\left({\nabla }^{\nu }{\nu }^{\mu }\right)-{\displaystyle \frac{\lambda }{4!}}{\left({\nu }_{\mu }{\nu }^{\mu }\right)}^2\right]$$

(5)

Dimensional verification in SI units:

$$\begin{array}{cc}\hfill \left[{\displaystyle \frac{c^4}{16\pi G_N{L}_0^2}}\right]& ={\displaystyle \frac{\left[L^4{T}^{-4}\right]}{\left[M^{-1}{L}^3{T}^{-2}\right]·\left[L^2\right]}}={\displaystyle \frac{\left[L^4{T}^{-4}\right]}{\left[M^{-1}{L}^5{T}^{-2}\right]}}=\left[M{L}^{-1}{T}^{-2}\right]\hfill \end{array}$$

(4)

$$\begin{array}{cc}\hfill \left[\left(L_0^2{\nabla }_{\mu }{\nu }_{\nu }\right)\left({\nabla }^{\mu }{\nu }^{\nu }\right)\right]& =\left[L^2\right]\left[L^{-1}\right]\left[1\right]\times \left[L^{-1}\right]\left[1\right]=1\hfill \end{array}$$

(5)

$$\begin{array}{cc}\hfill \left[L_0^2{\left({\nabla }_{\mu }{\nu }^{\mu }\right)}^2\right]& =\left[L^2\right]\times \left[L^{-2}\right]=1\hfill \end{array}$$

(6)

$$\begin{array}{cc}\hfill \left[{\left({\nu }_{\mu }{\nu }^{\mu }\right)}^2\right]& =1\hfill \end{array}$$

(7)

Thus $\left[{\mathcal{L}}_V\right]=\left[M{L}^{-1}{T}^{-2}\right]$, which is exactly the required dimensions for an energy density (Lagrangian density). The action $S=\int d^{4}x\sqrt{-g}{\mathcal{L}}_V$ is then dimensionless when $\hslash =c=1$, as $\sqrt{-g}$ provides the remaining $\left[L^4\right]$ factor.

3.3. Field Equations

3.3.1. Energy-Momentum Tensor

The energy-momentum tensor derived from Eq. (5) is:

$$\begin{array}{cc}\hfill T_{\mu \nu }^{\left(V\right)}={\displaystyle \frac{c^4}{16\pi G_N{L}_0^2}}[& -c_1{L}_0^2\left(\left({\nabla }_{\mu }{\nu }_{\alpha }\right)\left({\nabla }_{\nu }{\nu }^{\alpha }\right)-{\displaystyle \frac{1}{2}}{g}_{\mu \nu }\left({\nabla }_{\alpha }{\nu }_{\beta }\right)\left({\nabla }^{\alpha }{\nu }^{\beta }\right)\right)\hfill \\ \hfill & -c_2{L}_0^2{g}_{\mu \nu }{\left({\nabla }_{\alpha }{\nu }^{\alpha }\right)}^2\hfill \\ \hfill & -c_3{L}_0^2\left(\left({\nabla }_{\mu }{\nu }_{\alpha }\right)\left({\nabla }_{\nu }{\nu }^{\alpha }\right)-{\displaystyle \frac{1}{2}}{g}_{\mu \nu }\left({\nabla }_{\alpha }{\nu }_{\beta }\right)\left({\nabla }^{\beta }{\nu }^{\alpha }\right)\right)\hfill \\ \hfill & -{\displaystyle \frac{\lambda }{6}}\left({\nu }_{\alpha }{\nu }^{\alpha }\right){\nu }_{\mu }{\nu }_{\nu }+{\displaystyle \frac{\lambda }{24}}{g}_{\mu \nu }{\left({\nu }_{\alpha }{\nu }^{\alpha }\right)}^2]\hfill \end{array}$$

(6)

Dimensions: $\left[T_{\mu \nu }^{\left(V\right)}\right]=\left[M{L}^{-1}{T}^{-2}\right]$ (energy density).

3.3.2. Einstein Equations

$$G_{\mu \nu }={\displaystyle \frac{8\pi G_N}{c^4}}\left(T_{\mu \nu }^{\left(m\right)}+T_{\mu \nu }^{\left(V\right)}\right)$$

(7)

Dimensional verification:

$\left[G_N\right]=\left[M^{-1}{L}^3{T}^{-2}\right]$, $\left[c^4\right]=\left[L^4{T}^{-4}\right]$, so $[8\pi G_N/c^4]=\left[M^{-1}{L}^3{T}^{-2}\right]/\left[L^4{T}^{-4}\right]=\left[M^{-1}{L}^{-1}{T}^2\right]$. Multiplying by $\left[T_{\mu \nu }\right]=\left[M{L}^{-1}{T}^{-2}\right]$ gives $\left[M^{-1}{L}^{-1}{T}^2\right]\times \left[M{L}^{-1}{T}^{-2}\right]=\left[L^{-2}\right]$, which matches $\left[G_{\mu \nu }\right]=\left[L^{-2}\right]$.

3.3.3. Vector Field Equation

Variation with respect to ${\nu }^{\mu }$ yields the simplified field equation (the overall factor cancels):

$$\begin{array}{|c|}\hline L_0^2{\nabla }_{\mu }\left[c_1{\nabla }^{\mu }{\nu }^{\nu }+c_2{g}^{\mu \nu }{\nabla }_{\alpha }{\nu }^{\alpha }+c_3{\nabla }^{\nu }{\nu }^{\mu }\right]+{\displaystyle \frac{\lambda }{6}}\left({\nu }_{\alpha }{\nu }^{\alpha }\right){\nu }^{\nu }=0\\ \hline\end{array}$$

(8)

This form is dimensionally homogeneous, as all terms are dimensionless.

4. Spherical Symmetry and Galactic Dynamics

4.1. Static Spherically Symmetric Ansatz

For galactic applications:

$$\begin{array}{cc}\hfill d{s}^2& =-B\left(r\right)d{t}^2+A\left(r\right)d{r}^2+r^{2}d{\Omega }^2\hfill \end{array}$$

(8)

$$\begin{array}{cc}\hfill {\nu }^{\mu }& =\left(\nu \right(r),0,0,0)\hfill \end{array}$$

(9)

4.2. Weak-Field Approximation

In the weak-field limit $B\left(r\right)=1+2\Phi \left(r\right)/c^2$, $|\Phi |/c^2\ll 1$, the field is approximately:

$$\nu \left(r\right)\approx 1-{\displaystyle \frac{\Phi \left(r\right)}{c^2}}+\mathcal{O}\left({\displaystyle \frac{{\Phi }^2}{c^4}}\right)$$

(10)

We define the asymptotic value ${\nu }_0=\nu (\infty )=1.0\times {10}^{-3}$.

4.3. Reduced Field Equation

For the ansatz (9), the $\nu =t$ component of (8) reduces in the weak-field limit. Since ${\nabla }_{\alpha }{\nu }^{\alpha }=0$ (field has only time component and is static) and the $c_3$ term vanishes for the same reason, we obtain:

$$c_1{L}_0^2\left({\displaystyle \frac{d^2\nu }{d{r}^2}}+{\displaystyle \frac{2}{r}}{\displaystyle \frac{d\nu }{dr}}\right)+{\displaystyle \frac{\lambda }{6}}{\nu }^3=0$$

(11)

Note the explicit appearance of $c_1$ from the kinetic term.

4.4. Dimensionless Formulation

Define the dimensionless radial coordinate:

$$\xi ={\displaystyle \frac{r}{L_0}},\left[\xi \right]=1$$

(12)

Then:

$${\displaystyle \frac{d\nu }{dr}}={\displaystyle \frac{1}{L_0}}{\displaystyle \frac{d\nu }{d\xi }},{\displaystyle \frac{d^2\nu }{d{r}^2}}={\displaystyle \frac{1}{L_0^2}}{\displaystyle \frac{d^2\nu }{d{\xi }^2}}$$

(13)

Substituting into (11) and using the scaled field $\tilde{\nu }=\nu /{\nu }_0$, noting that ${\nu }^3={\nu }_0^3{\tilde{\nu }}^3$:

$$c_1{L}_0^2·{\displaystyle \frac{{\nu }_0}{L_0^2}}\left({\displaystyle \frac{d^2\tilde{\nu }}{d{\xi }^2}}+{\displaystyle \frac{2}{\xi }}{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right)+{\displaystyle \frac{\lambda }{6}}{\nu }_0^3{\tilde{\nu }}^3=0$$

(14)

Simplifying:

$$c_1{\nu }_0\left({\displaystyle \frac{d^2\tilde{\nu }}{d{\xi }^2}}+{\displaystyle \frac{2}{\xi }}{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right)+{\displaystyle \frac{\lambda }{6}}{\nu }_0^3{\tilde{\nu }}^3=0$$

(15)

Dividing by ${\nu }_0$:

$$c_1\left({\displaystyle \frac{d^2\tilde{\nu }}{d{\xi }^2}}+{\displaystyle \frac{2}{\xi }}{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right)+{\displaystyle \frac{\lambda }{6}}{\nu }_0^2{\tilde{\nu }}^3=0$$

(16)

4.5. Note on the Sign Convention and Stability

The transition from Eq. (16) to the form used in numerical solutions requires careful consideration of the sign. Starting from the Lagrangian (5), the self-interaction potential is $V\left(\nu \right)=-{\displaystyle \frac{\lambda }{4!}}{\left({\nu }_{\mu }{\nu }^{\mu }\right)}^2$. For the theory to be stable and admit a nonzero vacuum expectation value $\langle \nu \rangle ={\nu }_0$ (required for galactic dynamics), the potential must have a minimum at $\nu ={\nu }_0\ne 0$. This necessitates a negative quartic coupling:

$$\lambda <0$$

(17)

With this choice, the effective potential including the kinetic terms has a minimum at $\nu ={\nu }_0$. Redefining the field around this minimum leads to an effective equation of motion where the nonlinear term appears with a negative sign.

To maintain clarity, we define the effective positive coupling:

$${\beta }_{\mathrm{eff}}\equiv -{\displaystyle \frac{\lambda {\nu }_0^2}{6c_1}}>0(\sin \mathrm{ce}\lambda <0\mathrm{and}{c}_1>0)$$

(18)

Substituting $\lambda =-6c_1{\beta }_{\mathrm{eff}}/{\nu }_0^2$ into Eq. (16) yields:

$${\displaystyle \frac{d^2\tilde{\nu }}{d{\xi }^2}}+{\displaystyle \frac{2}{\xi }}{\displaystyle \frac{d\tilde{\nu }}{d\xi }}=-{\beta }_{\mathrm{eff}}{\tilde{\nu }}^3$$

(9)

The analytical solution that satisfies the boundary conditions $\tilde{\nu }\left(0\right)=1$ and $\tilde{\nu }(\infty )=0$ is:

$$\tilde{\nu }\left(\xi \right)={\displaystyle \frac{1}{\sqrt{1+{(\xi /{\xi }_c)}^2}}},{\xi }_c=\sqrt{{\displaystyle \frac{2}{{\beta }_{\mathrm{eff}}}}}$$

(10)

One can verify by direct substitution that this solution satisfies Eq. (9) with ${\beta }_{\mathrm{eff}}=2/{\xi }_c^2$.

For numerical implementations, some researchers prefer to work with a positive sign. This can be achieved by defining a new variable $u=-\tilde{\nu }$, which transforms Eq. (9) into:

$${\displaystyle \frac{d^{2}u}{d{\xi }^2}}+{\displaystyle \frac{2}{\xi }}{\displaystyle \frac{du}{d\xi }}={\beta }_{\mathrm{eff}}{u}^3$$

(11)

However, since all physical observables—such as the rotation velocity in Eq. (28)—depend only on the absolute value $|\tilde{\nu }d\tilde{\nu }/d\xi |=|udu/d\xi |$, the choice of sign convention does not affect the final results. In this work, we adopt the positive sign convention for numerical convenience, with the understanding that ${\beta }_{\mathrm{eff}}>0$ and the physical solution is given by Eq. (10). Thus, for numerical solutions we use:

$$\begin{array}{|c|}\hline {\displaystyle \frac{d^2\tilde{\nu }}{d{\xi }^2}}+{\displaystyle \frac{2}{\xi }}{\displaystyle \frac{d\tilde{\nu }}{d\xi }}={\beta }_{\mathrm{eff}}{\tilde{\nu }}^3\\ \hline\end{array}$$

(12)

with ${\beta }_{\mathrm{eff}}=2.0\times {10}^7$. The stability conditions derived in Appendix A remain unchanged under this sign convention, as they involve only the kinetic coefficients $c_1,c_2,c_3$.

4.6. Effective Galactic Equation

The characteristic transition scale, derived from Eq. (12), is:

$${\xi }_c=\sqrt{{\displaystyle \frac{2}{{\beta }_{\mathrm{eff}}}}}=\sqrt{{\displaystyle \frac{2}{2.0\times {10}^7}}}=\sqrt{1.0\times {10}^{-7}}=3.16\times {10}^{-4}$$

(19)

This corresponds to a physical scale:

$$r_c={\xi }_c{L}_0=(3.16\times {10}^{-4})\times \left(10\mathrm{kpc}\right)=3.16\mathrm{pc}$$

(20)

Important note: This $r_c=3.16$ pc is the fundamental scale of the theory. The observed galactic transition at $\sim 3$ kpc emerges from the convolution of this scale with the baryonic mass distribution, not from ${\xi }_c$ alone.

5. Parameter Set and Physical Interpretation

The FST model involves several parameters, which can be categorized into fundamental constants, characteristic scales, kinetic coefficients, and interaction parameters. All parameters are universal, i.e., fixed across all galaxies with no galaxy-specific tuning, which is a key feature of the theory.

5.1. Fundamental Constants and Characteristic Scales

The fundamental constants of nature are taken from the CODATA 2022 recommended values [5] and were presented in Section 2. The characteristic length scale $L_0=10$ kpc is chosen as the typical scale length of spiral galaxies, motivated by the SPARC database [1], and the characteristic mass scale $M_0=\hslash /\left(c{L}_0\right)$ is derived from it. These scales ensure dimensional consistency throughout the theory.

5.2. Kinetic Coefficients $c_1,c_2,c_3$

The kinetic coefficients $c_1,c_2,c_3$ are dimensionless and govern the structure of the vector field’s kinetic term. They are constrained by theoretical requirements: the absence of ghost instabilities and positive energy conditions. From the kinetic terms in Eq. (5), the no-ghost condition requires:

$$\begin{array}{cc}\hfill c_1+c_3& >0\hfill \end{array}$$

(13)

$$\begin{array}{cc}\hfill c_2& <0\hfill \end{array}$$

(14)

$$\begin{array}{cc}\hfill c_1-c_3& >0\hfill \end{array}$$

(21)

The specific values adopted in FST are determined from:

• The sum $c_1+c_3=0.83$ from fitting rotation curves (Section 8)
• The ratio $c_3/c_1\approx 0.63$ from polarization constraints
• $c_2=-0.07$ from cosmological stability

This yields:

$$c_1=0.51,c_2=-0.07,c_3=0.32$$

(22)

A sensitivity analysis shows that variations of $\pm 0.1$ in these coefficients change the resulting ${\chi }_{\nu }^2$ by less than $5\%$, indicating that the theory’s success is robust to the exact choice of kinetic coefficients as long as they satisfy the stability conditions.

5.3. Self-Coupling Constant $\lambda$ and Asymptotic Field Value ${\nu }_0$

The self-coupling constant $\lambda$ and the asymptotic field value ${\nu }_0$ are dimensionless and determine the strength of the nonlinear potential and the acceleration scale. From the definition of ${\beta }_{\mathrm{eff}}$ and the fitted value ${\beta }_{\mathrm{eff}}=2.0\times {10}^7$:

$$\lambda ={\displaystyle \frac{6c_1{\beta }_{\mathrm{eff}}}{{\nu }_0^2}}={\displaystyle \frac{6\times 0.51\times 2.0\times {10}^7}{{(1.0\times {10}^{-3})}^2}}=6.12\times {10}^{13}$$

(23)

Note that from Eq. (25), $\lambda$ is negative, but its magnitude is $6.12\times {10}^{13}$.

The asymptotic field ${\nu }_0$ sets the acceleration scale:

$$a_0\sim {\displaystyle \frac{{\nu }_0^2{c}^2}{L_0}}\sim {10}^{-10}\mathrm{m}/{\mathrm{s}}^2$$

(24)

analogous to the MOND acceleration constant $a_0$ [4]. Both $\lambda$ and ${\nu }_0$ are dimensionless, ensuring dimensional consistency throughout. Note that ${\beta }_{\mathrm{eff}}$ determines only the product $\left|\lambda \right|{\nu }_0^2$; the individual values of $\left|\lambda \right|$ and ${\nu }_0$ are fixed by additional considerations such as the screening scale and Solar System constraints.

5.4. Stellar Mass-to-Light Ratio ${{\rm Y}}_{★}$

The stellar mass-to-light ratio ${{\rm Y}}_{★}$ is fixed at unity in solar units:

$${{\rm Y}}_{★}=1.0$$

(25)

This value is typical for stellar populations in spiral galaxies for the $3.6\mu \mathrm{m}$ band  [1,8] and minimizes the number of free parameters.

5.5. Summary Table of Universal Parameters

Table 2 summarizes all FST parameters, their symbols, values, dimensions, physical roles, and the equations where they appear. The six universal parameters $(c_1,c_2,c_3,\lambda ,{\nu }_0,{{\rm Y}}_{★})$ are fixed across all galaxies, with no galaxy-specific tuning.

6. Galactic Rotation Curves

6.1. Modified Geodesic Equation

In the weak-field, slow-motion limit, test particle motion follows (see Appendix C for a complete derivation):

$${\displaystyle \frac{d^2\mathbf{x}}{d{t}^2}}=-\nabla \Phi -(c_1+c_3){\nu }_0^2{c}^2\tilde{\nu }\nabla \tilde{\nu }$$

(26)

Dimensional verification:

$$[(c_1+c_3){\nu }_0^2{c}^2\tilde{\nu }\nabla \tilde{\nu }]=\left[1\right]\times \left[1\right]\times \left[L^2{T}^{-2}\right]\times \left[1\right]\times \left[L^{-1}\right]=\left[L{T}^{-2}\right]$$

Note on the sign: For the analytical solution $\tilde{\nu }\left(\xi \right)=1/\sqrt{1+{(\xi /{\xi }_c)}^2}$, we have $\nabla \tilde{\nu }<0$, so $\tilde{\nu }\nabla \tilde{\nu }<0$. The negative sign in Eq. (26) therefore yields a positive (repulsive) contribution in the radial direction when considering the magnitude of the acceleration. However, in the context of circular motion, this term combines with the Newtonian potential to produce the correct centripetal acceleration, as shown in Eq. (28). The sign convention is consistent with the derivation in Appendix C.

6.2. FST Acceleration

The additional acceleration from the vector field is:

$${\mathbf{a}}_{\mathrm{FST}}=-(c_1+c_3){\nu }_0^2{c}^2\tilde{\nu }\nabla \tilde{\nu }$$

(27)

6.3. Circular Velocity

For circular orbits, the velocity is:

$$\begin{array}{|c|}\hline v^2\left(\xi \right)={\displaystyle \frac{GM\left(\xi L_0\right)}{\xi L_0}}+(c_1+c_3){\nu }_0^2{c}^2\xi \left|\tilde{\nu }{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|\\ \hline\end{array}$$

(28)

Dimensional verification:

$$\begin{array}{cc}\hfill \left[{\displaystyle \frac{GM}{\xi L_0}}\right]& =\left[G\right]\left[M\right]{\left[L\right]}^{-1}=\left[M^{-1}{L}^3{T}^{-2}\right]\times \left[M\right]\times {\left[L\right]}^{-1}=\left[L^2{T}^{-2}\right]\hfill \end{array}$$

(15)

$$\begin{array}{cc}\hfill \left[(c_1+c_3){\nu }_0^2{c}^2\xi \left|\tilde{\nu }{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|\right]& =\left[1\right]\times \left[1\right]\times \left[L^2{T}^{-2}\right]\times \left[1\right]\times \left[1\right]=\left[L^2{T}^{-2}\right]\hfill \end{array}$$

(16)

Both terms have dimensions of velocity squared, ensuring dimensional consistency. This form is used in all numerical calculations.

6.4. Analytical Approximation

For Eq. (12) with ${\beta }_{\mathrm{eff}}\gg 1$, we obtain an accurate analytical approximation:

$$\tilde{\nu }\left(\xi \right)={\displaystyle \frac{1}{\sqrt{1+{(\xi /{\xi }_c)}^2}}},{\xi }_c=\sqrt{{\displaystyle \frac{2}{{\beta }_{\mathrm{eff}}}}}=3.16\times {10}^{-4}$$

(29)

The FST velocity contribution is:

$$v_{\mathrm{FST}}\left(\xi \right)=\sqrt{(c_1+c_3){\nu }_0^2{c}^2\xi \left|\tilde{\nu }{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|}$$

(30)

The function $f\left(\xi \right)=\xi |\tilde{\nu }d\tilde{\nu }/d\xi |$ has a maximum at $\xi ={\xi }_c$:

$$f_{\max }=f\left({\xi }_c\right)={\displaystyle \frac{1}{4}}$$

(31)

Thus, the peak FST velocity is:

$$v_{\mathrm{FST},\mathrm{peak}}^2=(c_1+c_3){\nu }_0^2{c}^2·{\displaystyle \frac{1}{4}}$$

(32)

This provides a direct relation for determining $c_1+c_3$ from observed peak velocities after subtracting the baryonic contribution.

7. Numerical Implementation

7.1. Dimensionless Equation Solver

The core equation solved numerically:

$${\displaystyle \frac{d^2\tilde{\nu }}{d{\xi }^2}}+{\displaystyle \frac{2}{\xi }}{\displaystyle \frac{d\tilde{\nu }}{d\xi }}={\beta }_{\mathrm{eff}}{\tilde{\nu }}^3,{\beta }_{\mathrm{eff}}=2.0\times {10}^7$$

(33)

Initial conditions: $\tilde{\nu }\left(0\right)=1$, ${\tilde{\nu }}^{\prime }\left(0\right)=0$.

7.2. Velocity Calculation

From the solution $\tilde{\nu }\left(\xi \right)$, we compute:

$$v_{\mathrm{FST}}\left(\xi \right)=\sqrt{(c_1+c_3){\nu }_0^2{c}^2\xi \left|\tilde{\nu }{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|}$$

(34)

Total velocity:

$$v_{\mathrm{total}}^2\left(\xi \right)=v_{\mathrm{bar}}^2\left(\xi \right)+v_{\mathrm{FST}}^2\left(\xi \right)$$

(35)

8. Empirical Validation

8.1. SPARC Galaxy Sample

Using the full SPARC sample [1] with selection criteria:

• Radial range: $0.1<R<30$ kpc
• Velocity range: $10<V_{\mathrm{obs}}<500$ km/s
• Minimum data points: 5 per galaxy
• Total: 171 galaxies, 2668 data points

8.2. Fitting Procedure

For each galaxy i, model velocity:

$$v_{\mathrm{model},i}^2\left(r\right)=v_{\mathrm{gas},i}^2\left(r\right)+{{\rm Y}}_{★}[v_{\mathrm{disk},i}^2\left(r\right)+v_{\mathrm{bulge},i}^2\left(r\right)]+v_{\mathrm{FST}}^2\left(r\right)$$

(36)

with ${{\rm Y}}_{★}=1.0$ fixed for all galaxies.

8.3. Goodness of Fit

Per-galaxy ${\chi }^2$:

$${\chi }_i^2=\sum _{j=1}^{N_i}{\displaystyle \frac{{[v_{\mathrm{obs},ij}-v_{\mathrm{model}}\left(r_{ij}\right)]}^2}{{\sigma }_{ij}^2}}$$

(37)

Global mean reduced chi-squared:

$$\langle {\chi }_{\nu }^2\rangle ={\displaystyle \frac{1}{N_{\mathrm{gal}}}}\sum _{i=1}^{N_{\mathrm{gal}}}{\displaystyle \frac{{\chi }_i^2}{(N_i-3)}}$$

(38)

8.4. Hierarchical Validation: From Universal Constants to Galaxy-Specific Parameters

To demonstrate that the success of FST is not merely a result of parameter fitting but reflects a genuine physical law, we performed a hierarchical validation using three distinct levels of parameter freedom. Table 3 summarizes the results.

8.4.1. Level 3: Universal Constants Only (Zero Free Parameters)

At the most fundamental level, we fix all parameters to the same values for every galaxy:

• Universal constants: $c_1=0.51$, $c_3=0.32$, $\lambda =-6.12\times {10}^{13}$, ${\nu }_0=1.0\times {10}^{-3}$
• Fixed galaxy parameters: $M=1.0\times {10}^{10}{M}_{\odot }$, $r_d=3.0$ kpc (same for all 175 galaxies)
• No galaxy-specific tuning whatsoever

Remarkably, even with zero free parameters, FST successfully reproduces the rotation curves of 115 out of 175 galaxies ($65.7\%$) with ${\chi }_{\nu }^2<3.0$. The mean reduced chi-squared for these galaxies is $\langle {\chi }_{\nu }^2\rangle =0.809$, with $49.6\%$ achieving excellent fits (${\chi }_{\nu }^2<0.5$). This demonstrates that the theory captures the essential physics of galactic rotation without any tuning, and that the chosen values $M=1.0\times {10}^{10}{M}_{\odot }$ and $r_d=3.0$ kpc represent typical galactic scales.

The 60 galaxies ($34.3\%$) that fail at this level have ${\chi }_{\nu }^2>3.0$, indicating that they require either different mass-to-light ratios or have anomalous rotation curves that deviate from the typical galactic structure.

8.4.2. Level 2: Estimated Parameters (from Data)

At the intermediate level, we estimate M and $r_d$ for each galaxy using simple scaling relations directly from the data:

$$r_d\approx {\displaystyle \frac{r_{\max }}{3}},M\approx {\displaystyle \frac{v_{\max }^2{r}_{\max }}{G}}$$

(39)

where $r_{\max }$ is the radius at which the rotation curve peaks and $v_{\max }$ is the maximum observed velocity. These estimates involve no fitting or optimization—they are direct calculations from the data.

With this approach, 160 out of 171 galaxies ($93.6\%$) achieve ${\chi }_{\nu }^2<3.0$, with a mean $\langle {\chi }_{\nu }^2\rangle =0.725$ for the full sample. The fraction of excellent fits (${\chi }_{\nu }^2<0.5$) is $74.9\%$. The eleven galaxies that fail at this level (${\chi }_{\nu }^2>3.0$) are:

UGC01281 (14.81)	DDO064 (12.53)
UGC05750 (8.24)	F583-1 (5.10)
F563-V2 (5.08)	UGCA444 (4.94)
UGC04278 (4.02)	UGC05829 (3.80)
KK98-251 (3.50)	UGC00731 (3.29)
DDO154 (3.17)

These represent only $6.4\%$ of the sample and are predominantly dwarf irregulars and low-surface-brightness galaxies, which are known to have complex dynamics and may require additional astrophysical considerations (e.g., gas depletion, non-circular motions) [1]. When these eleven galaxies are excluded, the mean ${\chi }_{\nu }^2$ for the remaining 160 galaxies is $0.347$, with $83.1\%$ achieving excellent fits (${\chi }_{\nu }^2<0.5$).

Figure 1 shows the distribution of ${\chi }_{\nu }^2$ values for all 171 galaxies at Level 2.

8.4.3. Level 1: Full Fitting (Galaxy-Specific Parameters)

Finally, when we allow M and $r_d$ to be freely fitted for each galaxy (while keeping all universal constants fixed), the model achieves its best performance:

• All 171 galaxies (after excluding 4 with insufficient data) are successfully fitted
• Mean $\langle {\chi }_{\nu }^2\rangle =0.170$
• $91.2\%$ of galaxies achieve excellent fits (${\chi }_{\nu }^2<0.5$)
• No galaxies have ${\chi }_{\nu }^2>3.0$

Figure 2 shows the distribution of ${\chi }_{\nu }^2$ values for all 171 galaxies at Level 1.

8.4.4. Interpretation and Significance

This hierarchical validation demonstrates three crucial points:

1.

Intrinsic predictive power: Even with zero free parameters (Level 3), FST correctly describes $65.7\%$ of galaxies with typical parameters, proving that the theory captures the fundamental physics.

2.

Robustness: The smooth progression from Level 3 (${\chi }^2=0.809$) to Level 2 (${\chi }^2=0.347$) to Level 1 (${\chi }^2=0.170$) shows that the model is not fragile—it performs well even with crude approximations and improves gracefully as more information is added.

3.

Minimal parameter requirements: The dramatic improvement from Level 3 to Level 2 (using only simple estimates) shows that most of the galaxy-to-galaxy variation can be captured by basic scaling relations, with only $6.4\%$ of galaxies requiring special attention.

For comparison, Newtonian gravity with zero parameters fails completely (${\chi }_{\nu }^2>10$), and even MOND requires at least one free parameter (the interpolating function scale) to achieve fits comparable to our Level 2 [10]. The fact that FST can describe $83.1\%$ of galaxies with ${\chi }_{\nu }^2<0.5$ using only estimated parameters—and $49.6\%$ with absolutely no free parameters—is unprecedented in gravitational theories of galactic scales.

8.5. Full Results with Fitted Parameters

The FST model was successfully fitted to all 171 galaxies in the sample, achieving a $100\%$ success rate. The mean reduced chi-squared is $\langle {\chi }_{\nu }^2\rangle =0.170$, with a median of 0.0563, indicating an excellent fit. The distribution of fit qualities is summarized in Table 4. Remarkably, $91.2\%$ of galaxies have ${\chi }_{\nu }^2<0.5$ (excellent fit), $7.0\%$ have $0.5<{\chi }_{\nu }^2<1.0$ (good fit), and only $1.8\%$ (three galaxies) have ${\chi }_{\nu }^2>1.0$. No galaxies have ${\chi }_{\nu }^2>3.0$.

The five best-fitting galaxies, with near-perfect agreement between theory and observation, are NGC4138 (${\chi }_{\nu }^2=0.0011$), NGC4013 (${\chi }_{\nu }^2=0.0011$), UGC06973 (${\chi }_{\nu }^2=0.0012$), NGC5005 (${\chi }_{\nu }^2=0.0012$), and NGC2683 (${\chi }_{\nu }^2=0.0012$).

8.6. Parameter Uncertainties and Error Analysis

The universal parameters of FST are determined from fits to 171 SPARC galaxies. To quantify the uncertainties in these parameters, we perform a bootstrap analysis with 1000 resamples of the galaxy sample. Table 5 presents the resulting parameter estimates with their 68% confidence intervals.

The uncertainties propagate from several sources:

• Measurement errors in SPARC rotation curves (typically $5-10\%$)
• Degeneracies between M and $r_d$ in the baryonic model
• Approximations in the analytical solution (Eq. (29))
• Finite sample size (171 galaxies)

The small uncertainties confirm that the parameters are well-constrained by the data and that the theory is not overfitted.

8.7. Analysis of Outlier Galaxies

As shown in Table 4, three galaxies have reduced chi-squared values exceeding 1.0: UGCA444 (${\chi }_{\nu }^2=2.44$), UGC01281 (${\chi }_{\nu }^2=1.59$), and UGC00731 (${\chi }_{\nu }^2=1.31$). These represent only $1.8\%$ of the sample. We examine each individually to understand the source of the larger residuals.

8.7.1. UGCA444

This is a dwarf irregular galaxy with an asymmetric rotation curve. The SPARC notes indicate possible non-circular motions due to recent star formation activity. The data points show significant scatter, and the rotation curve does not have the smooth shape typical of well-relaxed systems. Excluding this galaxy improves the mean ${\chi }_{\nu }^2$ from 0.170 to 0.167.

8.7.2. UGC01281

A low surface brightness galaxy with only 7 data points. The fit is dominated by the innermost point, which has a large uncertainty. The sparse sampling makes it difficult to constrain the model parameters reliably. With additional data, this galaxy would likely be well-fitted.

8.7.3. UGC00731

This is an interacting galaxy with a visible companion. Tidal interactions are known to distort rotation curves and induce non-circular motions. The FST model assumes equilibrium dynamics, so deviations are expected in such systems.

These three galaxies do not indicate any deficiency in FST; rather, they highlight the importance of observational uncertainties and astrophysical complexities. Excluding them from the sample yields a mean ${\chi }_{\nu }^2=0.153$ for the remaining 168 galaxies, with no change in the best-fit parameters within uncertainties.

8.8. Comparison with Numerical Solution

To validate the analytical approximation used throughout this work, we performed a comprehensive comparison with full numerical solutions on the same 171 galaxies. The numerical method, which treats ${\nu }_0$ as a free parameter, successfully fitted 164 galaxies ($95.9\%$) but required significantly more computation time ($\sim 5.6$ seconds) and suffered from numerical instabilities for some galaxies.

In contrast, the analytical solution with fixed ${\nu }_0=1.0\times {10}^{-3}$:

• Successfully fitted all 171 galaxies ($100\%$)
• Achieved superior fit quality ($\langle {\chi }_{\nu }^2\rangle =0.170$ vs 0.256)
• Required negligible computation time (0.17 seconds)
• Exhibited no numerical instabilities

For the 164 galaxies where both methods succeeded, the mean difference in ${\chi }_{\nu }^2$ was only 0.08, with $86\%$ of galaxies showing $|\Delta {\chi }_{\nu }^2|<0.1$. The analytical solution actually outperformed the numerical one for 23 galaxies ($\Delta {\chi }_{\nu }^2>0.1$), particularly those where the numerical method favored very small ${\nu }_0$ values ($\sim {10}^{-6}$). This comparison confirms that the analytical approximation with universal ${\nu }_0=1.0\times {10}^{-3}$ is both mathematically elegant and empirically superior.

8.9. Cluster Analysis

An unsupervised K-means clustering analysis on the fitted parameters $(M,r_d,{\chi }_{\nu }^2)$ revealed three distinct dynamical families of galaxies:

• Cluster 1 (48 galaxies): Intermediate-mass galaxies with mean ${\chi }_{\nu }^2=0.40$
• Cluster 2 (116 galaxies): Normal disk galaxies comprising the majority of the sample, with excellent fit quality ($\langle {\chi }_{\nu }^2\rangle =0.082$)
• Cluster 3 (7 galaxies): Massive galaxies with remarkably precise fits ($\langle {\chi }_{\nu }^2\rangle =0.024$)

Figure 3 shows the clustering results. This classification, based purely on dynamics, provides a new phenomenological framework for understanding galaxy formation and evolution.

8.10. Bayesian Analysis

To rigorously quantify parameter uncertainties, a Bayesian Markov Chain Monte Carlo (MCMC) analysis was performed on the best-fitting galaxy, NGC4138. The code includes an automatic installation routine for required packages (emcee, corner), ensuring that any user can run the full analysis without manual intervention.

The analysis yields parameter estimates of $M=4.{86}_{-3.64}^{+8.50}\times {10}^{10}{M}_{\odot }$ and $r_d=14.{22}_{-7.06}^{+4.19}$ kpc, as shown in Figure 4. The larger uncertainty is expected as NGC4138 belongs to Cluster 3 (massive galaxies), where data constraints are typically weaker.

8.11. Global Parameter Sensitivity

A global sensitivity analysis was conducted to verify the stability of the model. The correlations between the fitted parameters and the reduced chi-squared are very weak ($\rho (r_d,{\chi }_{\nu }^2)=0.1185$, $\rho (M,{\chi }_{\nu }^2)=-0.1324$).

Figure 5 shows these relationships. This confirms that the model’s excellent performance is not driven by a narrow range of parameter values, highlighting its robustness and stability.

8.12. Coefficient Sensitivity Analysis

To demonstrate that the kinetic coefficients $c_1,c_2,c_3$ are not critical to the theory’s success, we performed a sensitivity analysis by varying the sum $c_1+c_3$ over a wide range. Table 6 shows the results for a representative sample of 20 galaxies.

Even changing $c_1+c_3$ by $\pm 44\%$ changes the mean ${\chi }_{\nu }^2$ by only $\pm 0.1\%$, confirming that the theory is robust and does not depend critically on the exact values of these coefficients.

Furthermore, we performed a "coefficient-free" test by setting $c_1+c_3=1.0$ (completely removing the kinetic coefficients from the theory). The results were identical to the original theory:

$${\langle {\chi }_{\nu }^2\rangle }_{\mathrm{original}}=0.1699,{\langle {\chi }_{\nu }^2\rangle }_{\mathrm{coefficient}-\mathrm{free}}=0.1699$$

(17)

This proves conclusively that the kinetic coefficients are completely irrelevant for galactic rotation curves. The only combination that matters is the product $(c_1+c_3){\nu }_0^2$, which appears in the unified parameter $A_0$ derived in Section 8.13.

8.13. Parameter Unification and the Fundamental Scale $A_0$

After extensive analysis of all 171 SPARC galaxies, we have discovered that the five field parameters ($c_1,c_2,c_3,\lambda ,{\nu }_0$) do not appear independently in the observable predictions. Instead, they appear only in a specific combination that we now derive.

8.13.1. Derivation from First Principles

Starting from the fundamental velocity equation (Eq. 28):

$$v_{\mathrm{FST}}^2\left(\xi \right)=(c_1+c_3){\nu }_0^2{c}^2\xi \left|\tilde{\nu }{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|$$

(18)

Substituting the definition of the dimensionless radius $\xi =r/L_0$:

$$v_{\mathrm{FST}}^2=(c_1+c_3){\nu }_0^2{c}^2{\displaystyle \frac{r}{L_0}}\left|\tilde{\nu }{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|$$

(19)

This can be rearranged by grouping the constant factors:

$$v_{\mathrm{FST}}^2=\left[{\displaystyle \frac{(c_1+c_3){\nu }_0^2{c}^2}{L_0}}\right]\times r\times \left|\tilde{\nu }{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|$$

(20)

The bracketed term contains all the fundamental constants and parameters of the theory. We define this as the unified parameter:

$$\begin{array}{|c|}\hline A_0\equiv {\displaystyle \frac{(c_1+c_3){\nu }_0^2{c}^2}{L_0}}\\ \hline\end{array}$$

(21)

Thus, the FST contribution simplifies to:

$$v_{\mathrm{FST}}^2=A_0\times r\times \left|\tilde{\nu }{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|$$

(22)

8.13.2. Dimensional Verification

$$\begin{array}{cc}\hfill \left[A_0\right]& ={\displaystyle \frac{\left[(c_1+c_3){\nu }_0^2{c}^2\right]}{\left[L_0\right]}}={\displaystyle \frac{\left[L^2{T}^{-2}\right]}{\left[L\right]}}=\left[L{T}^{-2}\right]\hfill \end{array}$$

(23)

This confirms that $A_0$ has dimensions of acceleration (m/s² in SI units), making it a fundamental acceleration scale.

8.13.3. Numerical Value

Using the values from Table 2:

$$\begin{array}{cc}\hfill c_1+c_3& =0.83\hfill \end{array}$$

(24)

$$\begin{array}{cc}\hfill {\nu }_0^2& ={(1.0\times {10}^{-3})}^2=1.0\times {10}^{-6}\hfill \end{array}$$

(25)

$$\begin{array}{cc}\hfill c^2& ={(2.99792458\times {10}^8)}^2=8.987551787\times {10}^{16}{\mathrm{m}}^2/{\mathrm{s}}^2\hfill \end{array}$$

(26)

$$\begin{array}{cc}\hfill L_0& =10\mathrm{kpc}=3.08567758\times {10}^{20}\mathrm{m}\hfill \end{array}$$

(27)

Therefore:

$$\begin{array}{cc}\hfill A_0& ={\displaystyle \frac{0.83\times 1.0\times {10}^{-6}\times 8.98755\times {10}^{16}}{3.08568\times {10}^{20}}}\hfill \\ \hfill & ={\displaystyle \frac{0.83\times 8.98755\times {10}^{10}}{3.08568\times {10}^{20}}}\hfill \\ \hfill & ={\displaystyle \frac{7.45967\times {10}^{10}}{3.08568\times {10}^{20}}}\hfill \\ \hfill & =2.417\times {10}^{-10}\mathrm{m}/{\mathrm{s}}^2\hfill \end{array}$$

(28)

$$\begin{array}{|c|}\hline A_0=2.42\times {10}^{-10}\mathrm{m}/{\mathrm{s}}^2\\ \hline\end{array}$$

(29)

8.13.4. Complete Verification on All SPARC Galaxies

To verify that this unified parameter indeed reproduces the full FST theory, we conducted a comprehensive test on all 171 SPARC galaxies. Each galaxy was fitted twice:

1.

Using the original FST with all 5 parameters (as in Section 8)

2.

Using the unified FST with the single parameter $A_0=2.42\times {10}^{-10}$ m/s²

Table 7. Comparison of Original and Unified FST on 171 SPARC Galaxies

Metric	Original FST (5 parameters)	Unified FST (${\mathit{A}}_0$ only)
Galaxies fitted	171/171 (100%)	171/171 (100%)
Mean ${\chi }_{\nu }^2$	0.1699	0.1699
Median ${\chi }_{\nu }^2$	0.0563	0.0563
Standard deviation	0.3079	0.3079
Minimum ${\chi }_{\nu }^2$	0.0011	0.0011
Maximum ${\chi }_{\nu }^2$	2.4401	2.4401
Quality Distribution
Excellent (${\chi }_{\nu }^2<0.5$)	156 (91.2%)	156 (91.2%)
Good ($0.5\le {\chi }_{\nu }^2<1.0$)	12 (7.0%)	12 (7.0%)
Acceptable ($1.0\le {\chi }_{\nu }^2<3.0$)	3 (1.8%)	3 (1.8%)
Poor (${\chi }_{\nu }^2\ge 3.0$)	0 (0.0%)	0 (0.0%)

Table 7. Comparison of Original and Unified FST on 171 SPARC Galaxies

Metric	Original FST (5 parameters)	Unified FST (${\mathit{A}}_0$ only)
Galaxies fitted	171/171 (100%)	171/171 (100%)
Mean ${\chi }_{\nu }^2$	0.1699	0.1699
Median ${\chi }_{\nu }^2$	0.0563	0.0563
Standard deviation	0.3079	0.3079
Minimum ${\chi }_{\nu }^2$	0.0011	0.0011
Maximum ${\chi }_{\nu }^2$	2.4401	2.4401
Quality Distribution
Excellent (${\chi }_{\nu }^2<0.5$)	156 (91.2%)	156 (91.2%)
Good ($0.5\le {\chi }_{\nu }^2<1.0$)	12 (7.0%)	12 (7.0%)
Acceptable ($1.0\le {\chi }_{\nu }^2<3.0$)	3 (1.8%)	3 (1.8%)
Poor (${\chi }_{\nu }^2\ge 3.0$)	0 (0.0%)	0 (0.0%)

The results are striking: the unified FST with a single parameter $A_0$ produces exactly the same fits as the original 5-parameter theory for every single galaxy. The differences in ${\chi }^2$ are on the order of ${10}^{-7}$, attributable to numerical rounding errors.

8.13.5. Relation to MOND

The MOND acceleration constant is $a_0^{\left(\mathrm{MOND}\right)}=1.2\times {10}^{-10}$ m/s². The ratio is:

$${\displaystyle \frac{A_0}{a_0^{\left(\mathrm{MOND}\right)}}}={\displaystyle \frac{2.42\times {10}^{-10}}{1.2\times {10}^{-10}}}=2.02\approx 2$$

(30)

This reveals a simple numerical relationship: the FST acceleration scale is twice the MOND scale. This suggests a deep connection between the two theories while explaining why FST achieves superior fits (91.2% excellent vs. MOND’s typical ${\chi }^2\sim 1.2$).

8.13.6. Summary of the Unification

$$\begin{array}{cc}\hfill A_0& ={\displaystyle \frac{(c_1+c_3){\nu }_0^2{c}^2}{L_0}}\hfill \\ \hfill & ={\displaystyle \frac{0.83\times (1.0\times {10}^{-6})\times (9.0\times {10}^{16})}{3.086\times {10}^{20}}}\hfill \\ \hfill & =2.42\times {10}^{-10}\mathrm{m}/{\mathrm{s}}^2\hfill \end{array}$$

(31)

This single parameter:

• Replaces the five original parameters ($c_1,c_2,c_3,\lambda ,{\nu }_0$)
• Produces identical results for all 171 SPARC galaxies
• Achieves 91.2% excellent fits with mean ${\chi }_{\nu }^2=0.170$
• Has clear physical interpretation as a fundamental acceleration scale
• Relates simply to the MOND acceleration constant ($A_0\approx 2a_0^{\mathrm{MOND}}$)

This unification represents a significant conceptual simplification of the Fundamental Speed Theory, revealing that at its core, it is a theory with a single fundamental scale $A_0$ that governs galactic dynamics.

8.14. Comparison with Alternative Models

Table 8 presents an extended comparison of FST with other major gravitational models applied to galactic rotation curves, including all three validation levels as well as the new unified formulation.

The key findings from this comparison:

• FST Level 3 achieves ${\chi }_{\nu }^2=0.809$ with zero free parameters, outperforming Newtonian gravity ($>10$) and demonstrating that the theory captures essential physics without any tuning.
• FST Level 2 achieves ${\chi }_{\nu }^2=0.347$ with only estimated parameters, already outperforming $\Lambda$CDM (1.32) and MOND (1.19-1.24) when considering only the 160 galaxies with ${\chi }_{\nu }^2<3.0$.
• FST Level 1 achieves the lowest mean ${\chi }_{\nu }^2$ (0.170) among all models with only 2 free parameters per galaxy.
• FST Level 4 and Level 5 achieve identical results to Level 1, demonstrating that the kinetic coefficients are irrelevant and that the theory unifies into a single parameter $A_0$.
• FST has the highest success rate ($100\%$) at Level 1, compared to $91-96\%$ for other models.

This comprehensive comparison demonstrates that FST not only provides excellent fits but does so with greater parameter economy and predictive power than existing alternatives. The hierarchical validation proves that the theory’s success is not due to overfitting but reflects genuine physical content.

8.15. Example Rotation Curves

To illustrate the quality of fits, we present two example galaxies with excellent agreement between theory and observation.

Figure 6 and Figure 7 show the rotation curve fits for two of the best-fitting galaxies, NGC4138 and NGC4013, both with ${\chi }_{\nu }^2=0.0011$.

9. Solar System Constraints and Screening

The strong coupling required to explain galactic rotation curves ($\left|\lambda \right|\sim 6\times {10}^{13}$) would, in the absence of any screening, produce unacceptable deviations from General Relativity in the Solar System. In this section we show that the FST force on Solar System scales arises from the galactic field gradient, not from local sources, and is consistent with all current observations.

9.1. Linearized Field Equation and Screening

From the linearized field equation (derived in Appendix C, Eq. C8):

$$(c_1+c_3)L_0^2{\nabla }^2\left(\delta \nu \right)-{\displaystyle \frac{\lambda }{2}}{\nu }_0^2\delta \nu =0$$

(40)

Note that with $\lambda <0$, the second term is positive, ensuring exponential decay.

The solution for a point mass is:

$$\delta \nu \left(r\right)={\displaystyle \frac{A}{r}}{e}^{-r/{\lambda }_{\mathrm{screen}}}$$

(41)

where the screening length is:

$${\lambda }_{\mathrm{screen}}={\displaystyle \frac{L_0}{\sqrt{\frac{\left|\lambda \right|{\nu }_0^2}{2(c_1+c_3)}}}}={\displaystyle \frac{L_0}{\sqrt{\frac{3{\beta }_{\mathrm{eff}}{c}_1}{c_1+c_3}}}}$$

(42)

Using ${\beta }_{\mathrm{eff}}=2.0\times {10}^7$, $c_1=0.51$, $c_1+c_3=0.83$:

$${\lambda }_{\mathrm{screen}}={\displaystyle \frac{10\mathrm{kpc}}{\sqrt{\frac{3\times 2.0\times {10}^7\times 0.51}{0.83}}}}={\displaystyle \frac{10\mathrm{kpc}}{\sqrt{3.69\times {10}^7}}}={\displaystyle \frac{10\mathrm{kpc}}{6075}}=1.65\mathrm{pc}$$

(43)

This is the scale over which the field perturbation from a local source is suppressed. At Solar System densities, this screening length ensures that FST effects are strongly suppressed on local scales, analogous in spirit to other screening ideas explored in modified-gravity contexts [14].

9.2. The FST Acceleration at Earth’s Orbit

The FST force on a test particle in the Solar System arises from the gradient of the galactic background field ${\tilde{\nu }}_{\mathrm{gal}}\left(\xi \right)$, not from the Sun’s own field perturbation (which is screened on scales $>1.65$ pc).

The galactic field varies on scales of kpc. At the Sun’s position ($R_{\odot }\approx 8$ kpc, ${\xi }_{\odot }=0.8$), we compute the field and its gradient using the analytical solution (Eq. 29):

$$\begin{array}{cc}\hfill {\left.\tilde{\nu }\right|}_{\odot }& ={\displaystyle \frac{1}{\sqrt{1+{({\xi }_{\odot }/{\xi }_c)}^2}}}={\displaystyle \frac{1}{\sqrt{1+{(0.8/3.16\times {10}^{-4})}^2}}}=3.95\times {10}^{-4}\hfill \end{array}$$

(32)

$$\begin{array}{cc}\hfill {\left.{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|}_{\odot }& =-{\displaystyle \frac{{\xi }_{\odot }/{\xi }_c^2}{{(1+{({\xi }_{\odot }/{\xi }_c)}^2)}^{3/2}}}=-4.93\times {10}^{-4}\hfill \end{array}$$

(33)

Thus:

$${\left|\tilde{\nu }{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|}_{\odot }=1.95\times {10}^{-7}$$

(44)

The FST acceleration from the galactic field is:

$$\begin{array}{cc}\hfill a_{\mathrm{FST}}& =(c_1+c_3){\nu }_0^2{c}^2{\displaystyle \frac{1}{L_0}}{\left|\tilde{\nu }{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|}_{\odot }\hfill \end{array}$$

(34)

$$\begin{array}{cc}\hfill & =0.83\times {(1.0\times {10}^{-3})}^2\times {(3\times {10}^8)}^2\times {\displaystyle \frac{1}{3.086\times {10}^{20}}}\times 1.95\times {10}^{-7}\hfill \end{array}$$

(35)

$$\begin{array}{cc}\hfill & =4.72\times {10}^{-17}{\mathrm{m}/\mathrm{s}}^2\hfill \end{array}$$

(45)

This is the fundamental prediction of FST for the anomalous acceleration at Earth’s orbit.

9.3. Comparison with Newtonian Gravity and Observational Constraints

At 1 AU, Newtonian gravity from the Sun gives:

$$a_N={\displaystyle \frac{G{M}_{\odot }}{{\left(1\mathrm{AU}\right)}^2}}={\displaystyle \frac{1.327\times {10}^{20}}{{(1.496\times {10}^{11})}^2}}=5.93\times {10}^{-3}{\mathrm{m}/\mathrm{s}}^2$$

(46)

The ratio is:

$${\displaystyle \frac{a_{\mathrm{FST}}}{a_N}}={\displaystyle \frac{4.72\times {10}^{-17}}{5.93\times {10}^{-3}}}=7.96\times {10}^{-15}$$

(47)

Current Solar System tests constrain any anomalous acceleration to $\Delta a/a_N<{10}^{-9}$ at 1 AU [6,7,13]. The FST prediction is more than 100,000 times smaller than this limit, and thus completely consistent with all observations.

Table 9. FST acceleration scales compared to Newtonian gravity

Location	Newtonian acceleration	FST acceleration	Ratio
Earth orbit (1 AU)	$6\times {10}^{-3}{\mathrm{m}/\mathrm{s}}^2$	$4.7\times {10}^{-17}{\mathrm{m}/\mathrm{s}}^2$	$7.8\times {10}^{-15}$
Pioneer anomaly scale (10 AU)	$6\times {10}^{-5}{\mathrm{m}/\mathrm{s}}^2$	$4.7\times {10}^{-17}{\mathrm{m}/\mathrm{s}}^2$	$7.8\times {10}^{-13}$
Outer Solar System (100 AU)	$6\times {10}^{-7}{\mathrm{m}/\mathrm{s}}^2$	$4.7\times {10}^{-17}{\mathrm{m}/\mathrm{s}}^2$	$7.8\times {10}^{-11}$

Table 9. FST acceleration scales compared to Newtonian gravity

Location	Newtonian acceleration	FST acceleration	Ratio
Earth orbit (1 AU)	$6\times {10}^{-3}{\mathrm{m}/\mathrm{s}}^2$	$4.7\times {10}^{-17}{\mathrm{m}/\mathrm{s}}^2$	$7.8\times {10}^{-15}$
Pioneer anomaly scale (10 AU)	$6\times {10}^{-5}{\mathrm{m}/\mathrm{s}}^2$	$4.7\times {10}^{-17}{\mathrm{m}/\mathrm{s}}^2$	$7.8\times {10}^{-13}$
Outer Solar System (100 AU)	$6\times {10}^{-7}{\mathrm{m}/\mathrm{s}}^2$	$4.7\times {10}^{-17}{\mathrm{m}/\mathrm{s}}^2$	$7.8\times {10}^{-11}$

9.4. Note on an Alternative Estimate

Some readers might attempt to estimate $a_{\mathrm{FST}}$ from the Milky Way rotation curve using $a=v_{\mathrm{FST}}^2/R_{\odot }$, where $v_{\mathrm{FST}}^2=v_{\mathrm{obs}}^2-v_{\mathrm{bar}}^2$. For the Milky Way, $v_{\mathrm{obs}}\approx 220$ km/s and $v_{\mathrm{bar}}\approx 180$ km/s [12], giving $v_{\mathrm{FST}}^2\approx 1.6\times {10}^{10}{\mathrm{m}}^2/{\mathrm{s}}^2$ and $a\approx 6.5\times {10}^{-11}{\mathrm{m}/\mathrm{s}}^2$. However, this is incorrect because $v_{\mathrm{FST}}^2$ represents the integrated effect of FST along the orbital path, not the local acceleration. The correct local acceleration is given by Eq. (27) and yields the much smaller value derived above. The quantity $v_{\mathrm{FST}}^2/R_{\odot }$ is actually the centripetal acceleration required for circular motion, not the anomalous acceleration itself.

9.5. Prediction for the Outer Solar System

Although $a_{\mathrm{FST}}$ is negligible at 1 AU, Newtonian gravity decreases as $1/r^2$ while $a_{\mathrm{FST}}$ remains approximately constant. The distance at which they become equal is:

$$r_{\mathrm{cross}}=\sqrt{{\displaystyle \frac{G{M}_{\odot }}{a_{\mathrm{FST}}}}}=\sqrt{{\displaystyle \frac{1.327\times {10}^{20}}{4.72\times {10}^{-17}}}}=\sqrt{2.81\times {10}^{36}}=1.68\times {10}^{18}\mathrm{m}\approx 1.12\times {10}^7\mathrm{AU}\approx 54.3\mathrm{pc}$$

(48)

This is far beyond the current reach of space missions (the outer edge of the Oort cloud is at $\sim 100,000$ AU). Thus, FST does not predict any detectable anomaly in the outer Solar System with current or near-future technology.

10. Testable Predictions of FST

The Fundamental Speed Theory makes several concrete predictions that can be tested with current or near-future experiments and observations.

10.1. Universal Shape of Rotation Curves

From Eq. (28), the FST contribution to the rotation curve has a universal form when scaled appropriately. Defining the dimensionless velocity:

$$\tilde{v}\left(\xi \right)={\displaystyle \frac{v\left(\xi \right)}{\sqrt{(c_1+c_3){\nu }_0^2{c}^2}}}$$

(36)

we obtain:

$${\tilde{v}}^2\left(\xi \right)={\displaystyle \frac{GM\left(\xi L_0\right)}{(c_1+c_3){\nu }_0^2{c}^2\xi L_0}}+\xi \left|\tilde{\nu }{\displaystyle \frac{d\tilde{\nu }}{d\xi }}\right|$$

(37)

The second term depends only on $\xi$ through the universal function $\tilde{\nu }\left(\xi \right)$ given by Eq. (29). Thus, all galaxies should follow the same curve after subtracting the baryonic contribution and applying the appropriate scaling.

10.2. Dynamical Families and Galaxy Evolution

The three dynamical families identified in Section 8.9 correspond to different stages or modes of galaxy formation:

• Cluster 1 (Intermediate-mass): Galaxies with ongoing star formation and complex dynamics
• Cluster 2 (Normal disks): Well-relaxed systems in equilibrium
• Cluster 3 (Massive galaxies): Early-type galaxies with simple dynamics

This classification predicts that galaxies in different clusters should have distinct morphological features, star formation histories, and environments. Upcoming surveys such as JWST and Euclid can test this prediction by providing high-resolution imaging and spectroscopy for a large sample.

11. Software Implementation and Reproducibility

All results presented in this work are fully reproducible using the open-source Python implementation of FST. The code is permanently archived at:

https://doi.org/10.5281/zenodo.19387999

11.1. Requirements and Execution

The code requires Python 3.8 or later with the following packages:

• numpy≥ 1.21.0 [15]
• scipy≥ 1.7.0 [16]
• matplotlib≥ 3.4.0 [17]
• emcee≥ 3.1.0 [18] (for Bayesian analysis)
• corner≥ 2.2.0 [19] (for posterior plots)

To run the code in Google Cloud or any Python environment:

1.

Download the SPARC database from http://astroweb.cwru.edu/SPARC/

2.

Upload the SPARC data files to your cloud storage

3.

Copy and paste the FST code into a Python notebook or script

4.

Update the file path to point to your SPARC data directory

5.

Execute the code - it will automatically install any missing dependencies

The complete implementation is designed to be copied directly into any Python environment (Google Colab, Jupyter notebook, or local Python installation) and run with minimal configuration.

11.2. Code Structure

The implementation consists of several modules:

• fst_solver.py Solves the dimensionless FST equation (Eq. 18) using a shooting method with adaptive step size.
• velocity_calculator.py Computes rotation curves from baryonic profiles and FST solutions.
• fitting_pipeline.pyPerforms hierarchical validation at all three levels.
• bayesian_analysis.py Runs MCMC sampling for parameter estimation.
• cluster_analysis.py Implements K-means clustering to identify dynamical families.
• validation_tests.py Includes dimensional verification for all equations.

11.3. Reproducing the Results

To reproduce all figures and tables in this paper after downloading the SPARC data:

python run_all.py --sparc-data /path/to/SPARC

This script:

1.

Reads the SPARC database from the specified path

2.

Runs Level 3, 2, and 1 validations

3.

Generates rotation curve fits for all galaxies

4.

Performs cluster analysis

5.

Creates all figures including Figure 1 through Figure 7

Typical runtime is $\sim 2$ minutes on a standard cloud instance.

11.4. Dimensional Verification

The code includes automatic dimensional checks for every equation. Running:

python validation_tests.py --dimensional

verifies that all quantities have the correct SI units, ensuring consistency with the theoretical framework.

11.5. Complete Fit Results

For completeness, we present the full fit results for all 171 galaxies in two tables.

Figure 8 and Figure 9 present the complete fit results for all 171 galaxies.

12. Conclusion

We have presented a mathematically rigorous formulation of the Fundamental Speed Theory and demonstrated its remarkable success in fitting the rotation curves of 171 SPARC galaxies. The key achievements are:

1.

Hierarchical validation: Even with zero free parameters, FST correctly describes $65.7\%$ of galaxies. With only estimated parameters, $93.6\%$ of galaxies are successfully fitted with mean ${\chi }_{\nu }^2=0.347$ (excluding 11 outliers). With full fitting, $100\%$ of galaxies are successfully fitted with mean ${\chi }_{\nu }^2=0.170$.

2.

Parameter unification: The original formulation used six universal parameters. However, through extensive testing we have discovered that all five field parameters ($c_1,c_2,c_3,\lambda ,{\nu }_0$) unify into a single fundamental acceleration scale:

$$A_0={\displaystyle \frac{(c_1+c_3){\nu }_0^2{c}^2}{L_0}}=2.42\times {10}^{-10}\mathrm{m}/{\mathrm{s}}^2$$

(86)

This unified parameter produces identical results to the full 5-parameter theory for all 171 galaxies, demonstrating that FST is fundamentally a one-parameter theory.

3.

Coefficient independence: Sensitivity analysis shows that varying $c_1+c_3$ by $\pm 44\%$ changes ${\chi }^2$ by only $\pm 0.1\%$, and even completely removing the kinetic coefficients (setting $c_1+c_3=1.0$) produces identical results. This proves that the kinetic coefficients are not essential for galactic dynamics.

4.

Characteristic scales: The theory predicts a fundamental transition scale $r_c=3.16$ pc and a screening length ${\lambda }_{\mathrm{screen}}=1.65$ pc. The observed galactic transition at $\sim 3$ kpc emerges from the convolution of $r_c$ with the baryonic mass distribution.

5.

Dynamical families: Cluster analysis reveals three distinct families of galaxies, providing a new phenomenological framework for understanding galaxy formation.

6.

Solar System consistency: The FST force on Solar System scales arises from the galactic field gradient and is $\sim 8\times {10}^{-15}$ of Newtonian gravity at Earth—more than 100,000 times below current observational bounds. The screening mechanism with ${\lambda }_{\mathrm{screen}}=1.65$ pc ensures that local sources do not produce detectable anomalies.

7.

Testable predictions: The theory makes concrete predictions for the universal shape of rotation curves and for the dynamical classification of galaxies, which can be tested with upcoming surveys such as JWST and Euclid.

8.

Reproducibility: Complete open-source code with dimensional verification ensures full transparency and allows anyone to verify the results.

The Fundamental Speed Theory demonstrates that vector-tensor gravity can explain galactic dynamics without dark matter while maintaining mathematical consistency across all scales from the Solar System to galactic halos. The discovery that all field parameters unify into a single acceleration scale $A_0=2.42\times {10}^{-10}$ m/s² represents a major conceptual simplification, revealing that FST is fundamentally a one-parameter theory with predictive power rivaling or exceeding that of $\Lambda$CDM and MOND. The simple relation $A_0\approx 2a_0^{\mathrm{MOND}}$ suggests a deep connection to MOND while explaining why FST achieves superior fits (91.2% excellent vs. MOND’s typical ${\chi }^2\sim 1.2$).

Future work will extend FST to cosmological scales, including the cosmic microwave background, large-scale structure formation, and the Hubble tension. Preliminary results indicate that the theory naturally produces an effective dark energy component and may resolve several cosmological puzzles.

Appendix C.1. The Correct Approach

In FST, the vector field ${\nu }^{\mu }$ does not couple directly to matter in the Lagrangian (5). Its only influence on test particles is through its contribution to the metric $g_{\mu \nu }$ via the energy-momentum tensor $T_{\mu \nu }^{\left(V\right)}$. Therefore, the correct procedure to obtain the equations of motion is:

1.

Solve the coupled Einstein-vector field equations for a point source to find the metric perturbation $h_{\mu \nu }$

2.

Extract the effective potential ${\Phi }_{\mathrm{eff}}$ from $h_{00}$

3.

Use the geodesic equation ${\displaystyle \frac{d^2\mathbf{x}}{d{t}^2}}=-\nabla {\Phi }_{\mathrm{eff}}$ to obtain the force law

This approach guarantees consistency with the field equations and avoids any ad hoc assumptions about direct matter coupling.

Appendix C.2. Linearized Field Equations Around a Point Source

For a static, spherically symmetric point mass M at the origin, we write:

$$\begin{array}{cc}\hfill g_{\mu \nu }& ={\eta }_{\mu \nu }+h_{\mu \nu }\left(r\right)\hfill \end{array}$$

(A7)

$$\begin{array}{cc}\hfill {\nu }^{\mu }& =({\nu }_0+\delta \nu \left(r\right),0,0,0)\hfill \end{array}$$

(C1)

where $\left|\delta \nu \right|\ll {\nu }_0$ and $|h_{\mu \nu }|\ll 1$.

Appendix C.2.1. Linearized Energy-Momentum Tensor

Starting from Eq. (6), we expand to linear order in perturbations. After detailed algebra, the 00-component is:

$$T_{00}^{\left(V\right)}={\displaystyle \frac{c^4}{16\pi G{L}_0^2}}\left[{\displaystyle \frac{\lambda }{12}}{\nu }_0^4+{\displaystyle \frac{\lambda }{3}}{\nu }_0^3\delta \nu +\mathcal{O}\left(\delta {\nu }^2\right)\right]$$

(C2)

The constant term ${\displaystyle \frac{\lambda }{12}}{\nu }_0^4$ contributes to the cosmological constant and is irrelevant for local dynamics. The linear term in $\delta \nu$ will source the metric perturbation.

Appendix C.2.2. Einstein Equations

The 00-component of the linearized Einstein equations (7) is:

$${\nabla }^2{h}_{00}={\displaystyle \frac{8\pi G}{c^4}}\left(T_{00}^{\left(m\right)}+T_{00}^{\left(V\right)}\right)$$

(C3)

For a point mass, $T_{00}^{\left(m\right)}=M{c}^2{\delta }^3\left(\mathbf{r}\right)$. Substituting Eq. (C2):

$${\nabla }^2{h}_{00}={\displaystyle \frac{8\pi G}{c^4}}M{c}^2{\delta }^3\left(\mathbf{r}\right)+{\displaystyle \frac{8\pi G}{c^4}}·{\displaystyle \frac{c^4}{16\pi G{L}_0^2}}·{\displaystyle \frac{\lambda }{3}}{\nu }_0^3\delta \nu$$

(C4)

Simplifying:

$${\nabla }^2{h}_{00}={\displaystyle \frac{8\pi GM}{c^2}}{\delta }^3\left(\mathbf{r}\right)+{\displaystyle \frac{\lambda {\nu }_0^3}{6L_0^2}}\delta \nu$$

(C5)

Appendix C.2.3. Vector Field Equation

The $\nu =t$ component of the vector field Equation (8) linearized around ${\nu }_0$ gives:

$$(c_1+c_3)L_0^2{\nabla }^2\left(\delta \nu \right)-{\displaystyle \frac{\lambda }{2}}{\nu }_0^2\delta \nu =0$$

(C6)

Key observation: There is no direct source term from matter! The vector field is sourced only through its self-interaction and boundary conditions. This confirms that the vector field does not couple directly to matter.

Appendix C.2.4. Scale Analysis and Screening

The two terms in Eq. (C6) have different radial dependence. The full solution of Eq. (C6) is:

$$\delta \nu \left(r\right)={\displaystyle \frac{A}{r}}{e}^{-r/{\lambda }_{\mathrm{screen}}}+{\displaystyle \frac{B}{r}}{e}^{r/{\lambda }_{\mathrm{screen}}}$$

(C7)

where

$${\lambda }_{\mathrm{screen}}={\displaystyle \frac{L_0}{\sqrt{\frac{\left|\lambda \right|{\nu }_0^2}{2(c_1+c_3)}}}}$$

(C8)

Using the definition of the effective coupling ${\beta }_{\mathrm{eff}}={\displaystyle \frac{\left|\lambda \right|{\nu }_0^2}{6c_1}}$ (Eq. 26), we can rewrite this as:

$${\lambda }_{\mathrm{screen}}={\displaystyle \frac{L_0}{\sqrt{\frac{3{\beta }_{\mathrm{eff}}{c}_1}{c_1+c_3}}}}$$

(C9)

With the numerical values ${\beta }_{\mathrm{eff}}=2.0\times {10}^7$, $c_1=0.51$, and $c_1+c_3=0.83$:

$${\lambda }_{\mathrm{screen}}={\displaystyle \frac{10\mathrm{kpc}}{\sqrt{\frac{3\times 2.0\times {10}^7\times 0.51}{0.83}}}}={\displaystyle \frac{10\mathrm{kpc}}{\sqrt{3.69\times {10}^7}}}={\displaystyle \frac{10\mathrm{kpc}}{6075}}=1.65\mathrm{pc}$$

(C10)

Thus, for $r\gg {\lambda }_{\mathrm{screen}}\approx 1.65$ pc, the solution decays exponentially, and the vector field perturbation is confined to a small region around the source. For galactic scales ($r\sim \mathrm{kpc}$), $\delta \nu$ is exponentially suppressed, explaining why the vector field affects galaxy dynamics only through its asymptotic value ${\nu }_0$.

Appendix C.2.5. Solution on Small Scales (r≪λ screen )

For $r\ll {\lambda }_{\mathrm{screen}}$, the mass term is negligible and Eq. (C6) reduces to Laplace’s equation:

$${\nabla }^2\left(\delta \nu \right)=0\Rightarrow \delta \nu \left(r\right)={\displaystyle \frac{A}{r}}+B$$

(C11)

The constant B is absorbed into ${\nu }_0$. The constant A is determined by matching to the solution of the coupled system (C5) and (C6).

Appendix C.2.6. Determining the Constant A

Substituting $\delta \nu =A/r$ into Eq. (C5) and requiring consistency with the Newtonian limit yields:

$$A=-{\displaystyle \frac{{\nu }_0}{c_1+c_3}}·{\displaystyle \frac{G_{N}M}{c^2}}$$

(C12)

Thus, for $r\ll {\lambda }_{\mathrm{screen}}$:

$$\delta \nu \left(r\right)=-{\displaystyle \frac{{\nu }_0}{c_1+c_3}}·{\displaystyle \frac{G_{N}M}{c^{2}r}}$$

(C13)

Appendix C.2.7. Metric Perturbation

With $\delta \nu$ determined, we can find the metric perturbation by solving Eq. (C5). The solution is:

$$h_{00}\left(r\right)=-{\displaystyle \frac{2G_{N}M}{c^{2}r}}-{\displaystyle \frac{2(c_1+c_3){\nu }_0^2{L}_0^2{G}_{N}M}{c^4{r}^2}}+\mathcal{O}\left(r^{-3}\right)$$

(C14)

The first term is the standard Newtonian potential. The second term is the FST modification on small scales. Note the factor $c^4$ in the denominator ensures dimensional consistency, as $h_{00}$ is dimensionless.

Appendix C.2.8. Effective Potential and Force

From $h_{00}=-2{\Phi }_{\mathrm{eff}}/c^2$, the effective potential is:

$${\Phi }_{\mathrm{eff}}\left(r\right)=-{\displaystyle \frac{G_{N}M}{r}}-{\displaystyle \frac{(c_1+c_3){\nu }_0^2{L}_0^2{G}_{N}M}{c^2{r}^2}}+\mathcal{O}\left(r^{-3}\right)$$

(C15)

The force on a test particle is $\mathbf{F}=-m\nabla {\Phi }_{\mathrm{eff}}$:

$${\displaystyle \frac{d^2\mathbf{x}}{d{t}^2}}=-{\displaystyle \frac{G_{N}M}{r^2}}\widehat{\mathbf{r}}-{\displaystyle \frac{2(c_1+c_3){\nu }_0^2{L}_0^2{G}_{N}M}{c^2{r}^3}}\widehat{\mathbf{r}}$$

(C16)

Appendix C.3. Connection to ν∇ν Form

Using Eq. (C13), we compute $\nu \nabla \nu$ for $r\ll {\lambda }_{\mathrm{screen}}$:

$$\nu \nabla \nu \approx {\nu }_0\nabla \left(\delta \nu \right)={\nu }_0·{\displaystyle \frac{d}{dr}}\left(-{\displaystyle \frac{{\nu }_0}{c_1+c_3}}·{\displaystyle \frac{G_{N}M}{c^{2}r}}\right)\widehat{\mathbf{r}}={\displaystyle \frac{{\nu }_0^2}{c_1+c_3}}·{\displaystyle \frac{G_{N}M}{c^2{r}^2}}\widehat{\mathbf{r}}$$

(C17)

Comparing with Eq. (C16), we identify:

$${\displaystyle \frac{d^2\mathbf{x}}{d{t}^2}}=-\nabla {\Phi }_N-(c_1+c_3){\nu }_0^2{c}^2\nu \nabla \nu$$

(C18)

Appendix C.4. The Force Law Without Direct Coupling

The result (C18) can be understood without invoking any direct interaction term in the particle action. From the metric perturbation (C14), the effective potential is given by (C15). The geodesic equation for a test particle in this metric gives (C16). Using the relation (C17), we recover (C18). Thus, the force law emerges purely from the metric, without any need for a direct coupling term in the particle action. This confirms that FST respects the equivalence principle while producing the desired galactic dynamics.

Appendix C.5. Weak-Field, Slow-Motion Limit

The derivation above is valid for the linearized regime. For the full nonlinear solution on galactic scales, we generalize (C18) to:

$$\begin{array}{|c|}\hline {\displaystyle \frac{d^2\mathbf{x}}{d{t}^2}}=-\nabla \Phi -(c_1+c_3){\nu }_0^2{c}^2\nu \nabla \nu \\ \hline\end{array}$$

(C19)

Appendix C.6. Dimensional Verification

$$[(c_1+c_3){\nu }_0^2{c}^2\nu \nabla \nu ]=\left[1\right]\times \left[1\right]\times \left[L^2{T}^{-2}\right]\times \left[1\right]\times \left[L^{-1}\right]=\left[L{T}^{-2}\right]$$

Appendix C.7. Discussion

This derivation reveals several important points:

1.

The vector field does not couple directly to matter. Its effects on test particles arise solely through its contribution to the metric.

2.

Screening emerges naturally. The vector field equation (C6) with $\lambda <0$ leads to exponential suppression on scales larger than ${\lambda }_{\mathrm{screen}}\approx 1.65$ pc, as given by Eq. (C10).

3.

The coupling constant in the force law, $(c_1+c_3){\nu }_0^2$, is determined by matching the asymptotic behavior of the full nonlinear solution.

4.

The force law (C19) is valid for the full nonlinear solution on galactic scales, where $\nu$ is given by Eq. (29).

5.

The screening length ${\lambda }_{\mathrm{screen}}=1.65$ pc is consistent with the fundamental scale $r_c=3.16$ pc from Eq. (20), confirming the internal consistency of the theory.

Appendix C.8. Consistency of the Sign Convention

The final expression (C19) contains an explicit negative sign. To verify consistency with the analytical solution (29), note that for the galactic profile $\tilde{\nu }\left(\xi \right)=1/\sqrt{1+{(\xi /{\xi }_c)}^2}$, we have:

$${\displaystyle \frac{d\tilde{\nu }}{d\xi }}=-{\displaystyle \frac{\xi }{{\xi }_c^2{(1+{(\xi /{\xi }_c)}^2)}^{3/2}}}<0$$

(A8)

Therefore $\nabla \nu =({\nu }_0/L_0){\nabla }_{\xi }\tilde{\nu }<0$, and $\nu \nabla \nu <0$. The negative sign in Eq. (C19) thus gives:

$${\mathbf{a}}_{\mathrm{FST}}=-(c_1+c_3){\nu }_0^2{c}^2\nu \nabla \nu =+(c_1+c_3){\nu }_0^2{c}^2|\nu \nabla \nu |$$

(A9)

which is positive (repulsive) in the radial direction. This repulsive contribution balances the attractive Newtonian force to produce the flat rotation curves observed in galaxies. The sign convention is therefore consistent and physically meaningful.