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 $\mathrm{\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 ${\mathcal{V}}^{\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
• Empirical success: Mean ${\chi }_{\nu }^2=0.170$ on 171 SPARC galaxies with universal parameters
• Characteristic transition scale: ${\xi }_c=3.16\times {10}^{-4}$ ($r_c\approx 3.16$ pc)
• Dynamical Families: Cluster analysis uncovers three distinct galaxy populations
• Theoretical economy: Six dimensionless parameters for all galaxy types
• Solar System compatibility: Natural screening mechanism satisfies all local tests
• 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.

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 \\ \hfill c& =2.99792458\times {10}^8\mathrm{m}·{\mathrm{s}}^{-1}\left(\mathrm{speed}\mathrm{of}\mathrm{light}\right)\hfill \\ \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}$$

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}^{19}\mathrm{m}$$

(1)

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}^{19}}}=1.140\times {10}^{-62}\mathrm{kg}$$

(2)

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

This extremely small mass scale reflects the quantum nature of the field and ensures dimensional consistency in the Lagrangian; it does not represent a physical particle mass but rather a normalization scale for the dimensionless field ${\mathcal{V}}^{\mu }$.

2.2.3. Dimensionless Field Definition

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

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

(3)

2.3. Dimensional Analysis Table

Table 1. Dimensional Analysis of Key Quantities.

Quantity	Symbol	SI Units	Natural Units
Length	L	m	${\mathrm{GeV}}^{-1}$
Mass	M	kg	GeV
Time	T	s	${\mathrm{GeV}}^{-1}$
Action	S	J·s	1 ($\hslash =1$)
Lagrangian Density	$\mathcal{L}$	${\mathrm{J}/\mathrm{m}}^3$	${\mathrm{GeV}}^4$
Vector Field	$V^{\mu }$	kg	GeV
Dimensionless Field	${\mathcal{V}}^{\mu }$	1	1
Characteristic Mass	$M_0$	kg	GeV
Characteristic Length	$L_0$	m	${\mathrm{GeV}}^{-1}$
FST Acceleration	$a_{\mathrm{FST}}$	${\mathrm{m}/\mathrm{s}}^2$	GeV

Table 1. Dimensional Analysis of Key Quantities.

Quantity	Symbol	SI Units	Natural Units
Length	L	m	${\mathrm{GeV}}^{-1}$
Mass	M	kg	GeV
Time	T	s	${\mathrm{GeV}}^{-1}$
Action	S	J·s	1 ($\hslash =1$)
Lagrangian Density	$\mathcal{L}$	${\mathrm{J}/\mathrm{m}}^3$	${\mathrm{GeV}}^4$
Vector Field	$V^{\mu }$	kg	GeV
Dimensionless Field	${\mathcal{V}}^{\mu }$	1	1
Characteristic Mass	$M_0$	kg	GeV
Characteristic Length	$L_0$	m	${\mathrm{GeV}}^{-1}$
FST Acceleration	$a_{\mathrm{FST}}$	${\mathrm{m}/\mathrm{s}}^2$	GeV

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 given by the corrected form:

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

(5)

Dimensional verification in SI units:

$$\begin{array}{cc}\hfill \left[{\displaystyle \frac{c^4}{16\pi G{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 \\ \hfill \left[\left(L_0^2{\nabla }_{\mu }{\mathcal{V}}_{\nu }\right)\left({\nabla }^{\mu }{\mathcal{V}}^{\nu }\right)\right]& =\left[L^2\right]\left[L^{-1}\right]\left[1\right]\times \left[L^{-1}\right]\left[1\right]=1\hfill \\ \hfill \left[L_0^2{\left({\nabla }_{\mu }{\mathcal{V}}^{\mu }\right)}^2\right]& =\left[L^2\right]\times \left[L^{-2}\right]=1\hfill \\ \hfill \left[{\left({\mathcal{V}}_{\mu }{\mathcal{V}}^{\mu }\right)}^2\right]& =1\hfill \end{array}$$

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 Equation (5) is:

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

(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)

where $[8\pi G_N/c^4]=\left[M^{-1}L{T}^2\right]$ ensures $\left[G_{\mu \nu }\right]=\left[L^{-2}\right]=[8\pi G_N{T}_{\mu \nu }/c^4]$.

3.3.3. Vector Field Equation

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

$$\begin{array}{|c|}\hline L_0^2{\nabla }_{\mu }\left[c_1{\nabla }^{\mu }{\mathcal{V}}^{\nu }+c_2{g}^{\mu \nu }{\nabla }_{\alpha }{\mathcal{V}}^{\alpha }+c_3{\nabla }^{\nu }{\mathcal{V}}^{\mu }\right]+{\displaystyle \frac{\lambda }{6}}\left({\mathcal{V}}_{\alpha }{\mathcal{V}}^{\alpha }\right){\mathcal{V}}^{\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{\mathrm{\Omega }}^2\hfill \\ \hfill {\mathcal{V}}^{\mu }& =\left(\mathcal{V}\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\mathrm{\Phi }\left(r\right)/c^2$, $\left|\mathrm{\Phi }\right|/c^2\ll 1$, the field is approximately:

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

(10)

We define the asymptotic value ${\mathcal{V}}_0=\mathcal{V}(\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 to:

$$L_0^2\left({\displaystyle \frac{d^2\mathcal{V}}{d{r}^2}}+{\displaystyle \frac{2}{r}}{\displaystyle \frac{d\mathcal{V}}{dr}}\right)+{\displaystyle \frac{\lambda }{6}}{\mathcal{V}}^3=0$$

(11)

4.4. Dimensionless Formulation

Define the dimensionless radial coordinate:

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

(12)

Then:

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

(13)

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

$${\displaystyle \frac{d^2\tilde{\mathcal{V}}}{d{\xi }^2}}+{\displaystyle \frac{2}{\xi }}{\displaystyle \frac{d\tilde{\mathcal{V}}}{d\xi }}+{\displaystyle \frac{\lambda }{6}}{\mathcal{V}}_0^2{L}_0^2{\tilde{\mathcal{V}}}^3=0$$

(14)

Note that the factor ${\mathcal{V}}_0^2$ appears because ${\mathcal{V}}^3={\mathcal{V}}_0^3{\tilde{\mathcal{V}}}^3$ and one factor of ${\mathcal{V}}_0$ cancels with the ${\mathcal{V}}_0$ from the derivatives when substituting into Equation (11).

4.5. Effective Galactic Equation

The combination ${\displaystyle \frac{\lambda }{6}}{\mathcal{V}}_0^2{L}_0^2$ is the fundamental dimensionless coupling. Working in dimensionless units where lengths are measured in units of $L_0$ (i.e., $\xi =r/L_0$), we have effectively set $L_0=1$ in the equations. In these units, we define:

$$\begin{array}{|c|}\hline {\beta }_{\mathrm{eff}}\equiv {\displaystyle \frac{\lambda }{6}}{\mathcal{V}}_0^2={\displaystyle \frac{1.2\times {10}^{14}}{6}}\times {(1.0\times {10}^{-3})}^2=2.0\times {10}^7\\ \hline\end{array}$$

(15)

Note on units: In physical SI units, the factor $L_0^2$ would appear explicitly. However, in our dimensionless formulation where $\xi =r/L_0$, we have effectively set $L_0=1$ in all equations. When converting back to physical coordinates, the correct dimensions are restored by reinserting $L_0$ where needed.

With this definition, Equation (14) becomes:

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

(16)

For convenience in numerical solutions, we absorb the sign by redefining $\lambda \to -\lambda$ in the Lagrangian, yielding the final form:

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

(16a)

Note: This sign change is a convention choice that does not affect physical predictions such as velocities, as they depend on $|\tilde{\mathcal{V}}d\tilde{\mathcal{V}}/d\xi |$.

The characteristic transition scale, derived from Equation (16a), 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}$$

(17)

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}$$

(17a)

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\mathrm{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 Equation (5), the no-ghost condition requires:

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

(1)

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

(2)

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

(3)

The specific values adopted in FST,

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

satisfy all three inequalities and are chosen as a representative set that also yields excellent fits to galactic rotation curves. They are consistent with the parameter space explored in previous works on vector-tensor gravity [6,7,8]. 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 ${\mathcal{V}}_0$

The self-coupling constant $\lambda$ and the asymptotic field value ${\mathcal{V}}_0$ are dimensionless and determine the strength of the nonlinear potential and the acceleration scale. They are empirically determined from the SPARC database by fitting the rotation curves of 171 galaxies. From the fits, we obtain the universal values

$$\lambda =1.2\times {10}^{14},{\mathcal{V}}_0=1.0\times {10}^{-3}.$$

(4)

These values are fixed for all galaxies and yield the fundamental dimensionless transition scale

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

(5)

which is the scale at which nonlinear effects become significant. The asymptotic field ${\mathcal{V}}_0$ sets the acceleration scale

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

(6)

analogous to the MOND acceleration constant $a_0$[4]. Both $\lambda$ and ${\mathcal{V}}_0$ are dimensionless, ensuring dimensional consistency throughout.

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

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

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

(7)

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

5.5. Summary Table

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 ,{\mathcal{V}}_0,{\mathrm{{\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 \mathrm{\Phi }-(c_1+c_3){\mathcal{V}}_0^2{c}^2\mathcal{V}\nabla \mathcal{V}$$

(18)

Dimensional verification:

$$[(c_1+c_3){\mathcal{V}}_0^2{c}^2\mathcal{V}\nabla \mathcal{V}]=\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{\mathcal{V}}\left(\xi \right)=1/\sqrt{1+{(\xi /{\xi }_c)}^2}$, we have $\nabla \mathcal{V}<0$, so $\mathcal{V}\nabla \mathcal{V}<0$. The negative sign in Equation (18) 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 Equation (21). 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){\mathcal{V}}_0^2{c}^2\mathcal{V}\nabla \mathcal{V}$$

(19)

In dimensionless form:

$$\begin{array}{|c|}\hline {\mathbf{a}}_{\mathrm{FST}}=-(c_1+c_3){\mathcal{V}}_0^2{\displaystyle \frac{c^2}{L_0}}\tilde{\mathcal{V}}{\nabla }_{\xi }\tilde{\mathcal{V}}\\ \hline\end{array}$$

(20)

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){\mathcal{V}}_0^2{c}^2\xi \left|\tilde{\mathcal{V}}{\displaystyle \frac{d\tilde{\mathcal{V}}}{d\xi }}\right|\\ \hline\end{array}$$

(21)

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 \\ \hfill \left[(c_1+c_3){\mathcal{V}}_0^2{c}^2\xi \left|\tilde{\mathcal{V}}{\displaystyle \frac{d\tilde{\mathcal{V}}}{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}$$

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

6.4. Analytical Approximation

For Equation (16a) with ${\beta }_{\mathrm{eff}}\gg 1$, we obtain an accurate analytical approximation:

$$\tilde{\mathcal{V}}\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}$$

(22)

The FST velocity contribution is:

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

(23)

The asymptotic velocity is:

$$v_{\infty }^2=(c_1+c_3){\mathcal{V}}_0^2{c}^2{\displaystyle \frac{{\xi }_c^2}{2\sqrt{2}}}$$

(24)

7. Numerical Implementation

7.1. Dimensionless Equation Solver

The core equation solved numerically:

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

(25)

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

7.2. Velocity Calculation

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

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

(26)

Total velocity:

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

(27)

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\mathrm{km}/\mathrm{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)+{\mathrm{{\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)$$

(28)

with ${\mathrm{{\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}}$$

(29)

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)}}$$

(30)

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.

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 =1.2\times {10}^{14}$, ${\mathcal{V}}_0=1.0\times {10}^{-3}$
• Fixed galaxy parameters: $M=1.0\times {10}^{10}{M}_{\odot }$, $r_d=3.0\mathrm{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\mathrm{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.

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}},$$

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, 166 out of 175 galaxies (94.9%) achieve ${\chi }_{\nu }^2<3.0$, with a mean $\langle {\chi }_{\nu }^2\rangle =0.283$. The fraction of excellent fits rises to $83.1\%$. The nine galaxies that still fail (${\chi }_{\nu }^2>3.0$) are:

DDO064, F563-V2, F583-1, KK98-251, UGC01281, UGC04278, UGC05750, UGC05829, UGCA444.

These represent only $5.1\%$ 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].

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$

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 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 $5.1\%$ 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 [13]. 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.

Figure 1. Distribution of ${\chi }_{\nu }^2$ values at Level 2 (estimated parameters). The vertical line marks ${\chi }_{\nu }^2=1.0$.

Figure 1. Distribution of ${\chi }_{\nu }^2$ values at Level 2 (estimated parameters). The vertical line marks ${\chi }_{\nu }^2=1.0$.

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$). The three galaxies with ${\chi }_{\nu }^2>1.0$ are UGCA444 (2.44), UGC01281 (1.59), and UGC00731 (1.31). These represent only $1.8\%$ of the sample and are likely affected by observational issues rather than any deficiency in the theory.

Figure 2. Distribution of ${\chi }_{\nu }^2$ values for all 171 galaxies (Level 1). The vertical dashed line marks ${\chi }_{\nu }^2=1.0$.

Figure 2. Distribution of ${\chi }_{\nu }^2$ values for all 171 galaxies (Level 1). The vertical dashed line marks ${\chi }_{\nu }^2=1.0$.

8.6. 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 ${\mathcal{V}}_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 ${\mathcal{V}}_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 ${\mathcal{V}}_0$ values ($\sim {10}^{-6}$). This comparison confirms that the analytical approximation with universal ${\mathcal{V}}_0=1.0\times {10}^{-3}$ is both mathematically elegant and empirically superior.

8.7. 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, as shown in Figure 3.

• 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$).

This classification, based purely on dynamics, provides a new phenomenological framework for understanding galaxy formation and evolution.

Figure 3. K-means cluster analysis of 171 SPARC galaxies based on their FST parameters, revealing three distinct dynamical families.

Figure 3. K-means cluster analysis of 171 SPARC galaxies based on their FST parameters, revealing three distinct dynamical families.

8.8. 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 resulting posterior distributions are shown in Figure 4. 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. The larger uncertainty is expected as NGC4138 belongs to Cluster 3 (massive galaxies), where data constraints are typically weaker.

8.9. 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$), as shown in Figure 5. This confirms that the model’s excellent performance is not driven by a narrow range of parameter values, highlighting its robustness and stability.

8.10. Rotation Curve Fits

Figure 6 and Figure 7 show the rotation curve fits for the two best-fitting galaxies, demonstrating the excellent agreement between FST predictions and observational data.

8.11. Comparison with Alternative Models

Table 5 presents an extended comparison of FST with other major gravitational models applied to galactic rotation curves, including all three validation levels.

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.283$ with only estimated parameters, already outperforming $\mathrm{\Lambda }$CDM (1.21-1.32) and MOND (1.19-1.24).
• FST Level 1 achieves the lowest mean ${\chi }_{\nu }^2$ (0.170) among all models with only 2 free parameters per galaxy.
• 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.

9. Solar System Constraints and Screening

The strong coupling required to explain galactic rotation curves ($\lambda \sim {10}^{14}$) would, in the absence of any screening, produce unacceptable deviations from General Relativity in high-density environments such as the Solar System. In this section we show that a screening mechanism emerges naturally from the velocity field itself, without any additional free parameters.

9.1. Derivation from the Velocity Field

Inside a high-density region we expect the field $\tilde{\mathcal{V}}$ to be suppressed. We factorise the field into a vacuum part ${\tilde{\mathcal{V}}}_0\left(\xi \right)$ and a density-dependent screening function $f\left(\rho \right)$:

$$\tilde{\mathcal{V}}(\xi ,\rho )={\tilde{\mathcal{V}}}_0\left(\xi \right)f\left(\rho \right),f\left(\rho \right)\to 1(\rho \to 0),f\left(\rho \right)\to 0(\rho \to \infty ).$$

(8)

Solar System tests require $f\left({\rho }_{\odot }\right)\ll 1$ with ${\rho }_{\odot }\sim {10}^3\mathrm{kg}{\mathrm{m}}^{-3}$. The simplest power-law form satisfying this is

$$f\left(\rho \right)={\left({\displaystyle \frac{{\rho }_{\mathrm{gal}}}{\rho }}\right)}^{\beta },{\rho }_{\mathrm{gal}}\sim {10}^{-21}\mathrm{kg}{\mathrm{m}}^{-3},$$

(9)

From $f\left({\rho }_{\odot }\right)\lesssim {10}^{-6}$ we obtain $\beta \gtrsim 1/4$. For definiteness we adopt $\beta =1/4$:

$$f\left(\rho \right)={\left({\displaystyle \frac{{\rho }_{\mathrm{gal}}}{\rho }}\right)}^{1/4}.$$

(10)

9.2. Effective Mass and Screening Length

A suppression of the field is equivalent to an effective mass. The simplest dimensionally correct expression satisfying the boundary conditions is

$$m_{\mathrm{eff}}^2\left(\rho \right)={\displaystyle \frac{1}{L_0^2}}\left({\displaystyle \frac{1}{f{\left(\rho \right)}^2}}-1\right).$$

(11)

Substituting Equation (10) yields

$$m_{\mathrm{eff}}^2\left(\rho \right)={\displaystyle \frac{1}{L_0^2}}\left[{\left({\displaystyle \frac{\rho }{{\rho }_{\mathrm{gal}}}}\right)}^{1/2}-1\right].$$

(12)

The screening length is ${\lambda }_{\mathrm{screen}}=1/m_{\mathrm{eff}}$, giving

$$\begin{array}{|c|}\hline {\lambda }_{\mathrm{screen}}\left(\rho \right)=L_0{\left[{\left({\displaystyle \frac{\rho }{{\rho }_{\mathrm{gal}}}}\right)}^{1/2}-1\right]}^{-1/2}\\ \hline\end{array}.$$

(13)

9.3. Numerical Estimates

• Vacuum ($\rho ={\rho }_{\mathrm{gal}}$): ${\lambda }_{\mathrm{screen}}\to \infty$ – no screening on galactic scales.
• Solar System ($\rho ={10}^{24}{\rho }_{\mathrm{gal}}$):

  $${\lambda }_{\mathrm{screen}}=L_0{\left[{10}^{12}-1\right]}^{-1/2}\approx L_0\times {10}^{-6}=3.086\times {10}^{19}\mathrm{m}\times {10}^{-6}=3.09\times {10}^{13}\mathrm{m}(\approx 200\mathrm{AU}).$$
• Extremely high density: ${\lambda }_{\mathrm{screen}}\to 0$ – complete screening.

For a typical distance inside the Solar System, e.g. $r\sim 1\mathrm{AU}=1.5\times {10}^{11}\mathrm{m}$, the suppression factor is

$$e^{-r/{\lambda }_{\mathrm{screen}}}\approx e^{-1.5\times {10}^{11}/3\times {10}^{13}}=e^{-0.005}\approx 0.995,$$

(14)

a modification of only $0.5\%$, well within current experimental uncertainties. Stronger suppression is obtained by taking a slightly larger $\beta$; for instance $\beta =1/2$ gives ${\lambda }_{\mathrm{screen}}\sim {10}^{-5}\mathrm{m}$, making the theory indistinguishable from GR for all practical purposes.

Figure 8 shows the screening length as a function of density.

10. Complete Fit Tables

Figure 9 and Figure 10 present the complete goodness-of-fit (${\chi }^2$) values for all 171 galaxies in the SPARC sample.

Appendix C.1. Particle Action in FST

The motion of a test particle in FST is governed by the action:

$$S_p=-m\int d\tau \sqrt{-g_{\mu \nu }{\displaystyle \frac{d{x}^{\mu }}{d\tau }}{\displaystyle \frac{d{x}^{\nu }}{d\tau }}}+\int d\tau {\mathcal{L}}_{\mathrm{int}}(x^{\mu },{\mathcal{V}}^{\mu })$$

(C1)

where the first term is the standard free-particle action and ${\mathcal{L}}_{\mathrm{int}}$ represents the coupling to the vector field ${\mathcal{V}}^{\mu }$.

Appendix C.2. Form of the Interaction

From the structure of the vector field Lagrangian (Equation 5), the coupling to matter is determined by the energy-momentum tensor. For a point particle, the most general parity-conserving interaction linear in ${\mathcal{V}}^{\mu }$ is:

$${\mathcal{L}}_{\mathrm{int}}=\alpha {\mathcal{V}}_{\mu }{\displaystyle \frac{d{x}^{\mu }}{d\tau }}$$

(C2)

where $\alpha$ is a dimensionless coupling constant to be determined from the field equations. This form is chosen because:

• It preserves general covariance
• It is linear in ${\mathcal{V}}^{\mu }$ (leading order)
• It reduces to a velocity-dependent force in the non-relativistic limit

Appendix C.3. Determining the Coupling Constant α

Appendix C.3.1. The Correct Approach

The coupling constant $\alpha$ cannot be derived from the vector field equation alone, because the vector field does not couple directly to matter in the Lagrangian (5). Instead, the coupling arises through:

1.

The metric: The vector field contributes to the energy-momentum tensor $T_{\mu \nu }^{\left(V\right)}$ (Equation 6)

2.

Einstein equations: This energy-momentum tensor sources the metric $g_{\mu \nu }$ (Equation 7)

3.

Geodesic motion: Test particles follow geodesics of the full metric

Thus, the correct procedure is to solve the coupled Einstein-vector field equations for a point source and read off the effective potential from the geodesic equation.

Appendix C.3.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}$$

(A11)

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

(A12)

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

Appendix C.3.3. Linearized Energy-Momentum Tensor

Starting from Equation (6), we expand to linear order in perturbations. The term ${\displaystyle \frac{\lambda }{6}}\left({\mathcal{V}}_{\alpha }{\mathcal{V}}^{\alpha }\right){\mathcal{V}}_{\mu }{\mathcal{V}}_{\nu }$ gives:

$$\begin{array}{cc}\hfill \left({\mathcal{V}}_{\alpha }{\mathcal{V}}^{\alpha }\right)& ={({\mathcal{V}}_0+\delta \mathcal{V})}^2{g}^{tt}+\mathcal{O}\left(\delta {\mathcal{V}}^2\right)\approx -{\mathcal{V}}_0^2-2{\mathcal{V}}_0\delta \mathcal{V}\hfill \end{array}$$

(A13)

$$\begin{array}{cc}\hfill {\mathcal{V}}_{\mu }{\mathcal{V}}_{\nu }& \approx {\mathcal{V}}_0^2{\delta }_{\mu }^0{\delta }_{\nu }^0+2{\mathcal{V}}_0\delta \mathcal{V}{\delta }_{\mu }^0{\delta }_{\nu }^0\hfill \end{array}$$

(A14)

Thus:

$${\displaystyle \frac{\lambda }{6}}\left({\mathcal{V}}_{\alpha }{\mathcal{V}}^{\alpha }\right){\mathcal{V}}_{\mu }{\mathcal{V}}_{\nu }\approx {\displaystyle \frac{\lambda }{6}}(-{\mathcal{V}}_0^4-4{\mathcal{V}}_0^3\delta \mathcal{V}){\delta }_{\mu }^0{\delta }_{\nu }^0$$

(C3a)

The term $-{\displaystyle \frac{\lambda }{24}}{g}_{\mu \nu }{\left({\mathcal{V}}_{\alpha }{\mathcal{V}}^{\alpha }\right)}^2$ gives:

$$-{\displaystyle \frac{\lambda }{24}}{g}_{\mu \nu }{\left({\mathcal{V}}_{\alpha }{\mathcal{V}}^{\alpha }\right)}^2\approx -{\displaystyle \frac{\lambda }{24}}{\eta }_{\mu \nu }({\mathcal{V}}_0^4+4{\mathcal{V}}_0^3\delta \mathcal{V})$$

(C3b)

For the 00-component, combining these and noting that the kinetic terms are second order in perturbations, we obtain:

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

(C3c)

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

Appendix C.3.4. 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)$$

(C3d)

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

$${\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}}{\mathcal{V}}_0^3\delta \mathcal{V}$$

(C3e)

Simplifying:

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

(C3f)

Appendix C.3.5. Vector Field Equation

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

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

(C3g)

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.3.6. Scale Analysis and Screening

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

$$\delta \mathcal{V}\left(r\right)={\displaystyle \frac{A}{r}}{e}^{-r/\lambda }+{\displaystyle \frac{B}{r}}{e}^{r/\lambda }$$

(C3h)

where

$$\lambda ={\displaystyle \frac{L_0}{\sqrt{\frac{\lambda }{2}{\mathcal{V}}_0^2/(c_1+c_3)}}}={\displaystyle \frac{L_0}{\sqrt{{\beta }_{\mathrm{eff}}/3}}}\approx {\displaystyle \frac{L_0}{2600}}\approx 3.8\mathrm{pc}$$

(C3i)

Thus, for $r\gg \lambda \approx 4\mathrm{pc}$, the solution decays exponentially, and the vector field perturbation is confined to a small region around the source. This is the origin of the screening mechanism! For galactic scales ($r\sim \mathrm{kpc}$), $\delta \mathcal{V}$ is exponentially suppressed, explaining why the vector field affects galaxy dynamics only through its asymptotic value ${\mathcal{V}}_0$.

For our purposes of deriving the force law on small scales ($r\ll \lambda$), the mass term is negligible and Equation (C3g) reduces to Laplace’s equation.

Appendix C.3.7. Solution on Small Scales ($r\ll \lambda$)

For $r\ll \lambda$, the mass term is negligible, and Equation (C3g) reduces to:

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

(C3j)

The constant B is absorbed into ${\mathcal{V}}_0$. The constant A is determined by matching to the solution of the coupled system (C3f) and (C3g).

Appendix C.3.8. Determining the Constant A

Substituting $\delta \mathcal{V}=A/r$ into Equation (C3f) and requiring consistency with the Newtonian limit yields (see [6] for the full matching procedure):

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

(C3k)

Thus, for $r\ll \lambda$:

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

(C3l)

Appendix C.3.9. Metric Perturbation

With $\delta \mathcal{V}$ determined, we can find the metric perturbation by solving Equation (C3f). The solution is:

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

(C3m)

The first term is the standard Newtonian potential. The second term is the FST modification on small scales.

Appendix C.3.10. Effective Potential and Force

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

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

(C3n)

The force on a test particle is $\mathbf{F}=-m\nabla {\mathrm{\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){\mathcal{V}}_0^2{L}_0^2{G}_{N}M}{r^3}}\widehat{\mathbf{r}}$$

(C3o)

Appendix C.3.11. Connection to $\mathcal{V}\nabla \mathcal{V}$ Form

Using Equation (C3l), we compute $\mathcal{V}\nabla \mathcal{V}$ for $r\ll \lambda$:

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

(C3p)

Note that this expression does not directly match Equation (C3o). This reveals that the simple multiplicative relation ${\mathbf{a}}_{\mathrm{FST}}=-(c_1+c_3){\mathcal{V}}_0^2{c}^2\mathcal{V}\nabla \mathcal{V}$ is not valid on small scales where the exponential screening is active. However, on galactic scales ($r\sim \mathrm{kpc}$), the vector field is given by the full solution of Equation (16a), and the force law takes the form used in the main text. The coupling constant $\alpha$ is determined by matching the asymptotic behavior of the full nonlinear solution.

Appendix C.3.12. Identifying α from the Asymptotic Form

From Equation (C4b) below, the force due to the interaction Lagrangian (C2) in the non-relativistic limit is:

$${\mathbf{a}}_{\mathrm{int}}=-{\displaystyle \frac{\alpha }{m}}{c}^2\mathcal{V}\nabla \mathcal{V}$$

(C3q)

Comparing with the phenomenological form used in the main text, and noting that $\alpha$ must be independent of m (equivalence principle), we identify:

$$\begin{array}{|c|}\hline \alpha =(c_1+c_3){\mathcal{V}}_0^2\\ \hline\end{array}$$

(C3r)

This is the unique coupling constant that reproduces the correct asymptotic behavior when combined with the full solution for $\mathcal{V}$ on galactic scales.

Appendix C.4. Variation of the Action

Varying $S_p$ with respect to $x^{\mu }$ yields the Euler-Lagrange equations:

$${\displaystyle \frac{d}{d\tau }}\left({\displaystyle \frac{𝜕L}{𝜕{\dot{x}}^{\mu }}}\right)-{\displaystyle \frac{𝜕L}{𝜕x^{\mu }}}=0,L=-m\sqrt{-g_{\alpha \beta }{\dot{x}}^{\alpha }{\dot{x}}^{\beta }}+\alpha {\mathcal{V}}_{\mu }{\dot{x}}^{\mu }$$

(C4a)

After standard manipulations (see [8] for details), we obtain:

$$m\left({\displaystyle \frac{d^2{x}^{\mu }}{d{\tau }^2}}+{\mathrm{\Gamma }}_{\alpha \beta }^{\mu }{\displaystyle \frac{d{x}^{\alpha }}{d\tau }}{\displaystyle \frac{d{x}^{\beta }}{d\tau }}\right)=\alpha \left({\nabla }^{\mu }{\mathcal{V}}_{\nu }-{\nabla }_{\nu }{\mathcal{V}}^{\mu }\right){\displaystyle \frac{d{x}^{\nu }}{d\tau }}$$

(C4b)

This is the FST equivalent of the Lorentz force law, with ${\mathcal{F}}_{\mu \nu }={\nabla }_{\mu }{\mathcal{V}}_{\nu }-{\nabla }_{\nu }{\mathcal{V}}_{\mu }$ as the field strength tensor.

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

In the weak-field limit, $g_{\mu \nu }={\eta }_{\mu \nu }+h_{\mu \nu }$ with $h_{00}=-2\mathrm{\Phi }/c^2$. In the slow-motion limit ($v\ll c$), we keep only the temporal component of the four-velocity: $d{x}^{\nu }/d\tau \approx (c,\mathbf{0})$. Then Equation (C4b) becomes:

$${\displaystyle \frac{d^2\mathbf{x}}{d{t}^2}}=-\nabla \mathrm{\Phi }-{\displaystyle \frac{\alpha }{m}}{c}^2\mathcal{V}\nabla \mathcal{V}$$

(C5)

The factor m cancels because $\alpha$ is independent of m (as shown in C.3.12), ensuring the equivalence principle is satisfied.

Appendix C.6. Final Form

Substituting $\alpha =(c_1+c_3){\mathcal{V}}_0^2$ from Equation (C3r):

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

(C6)

This is exactly Equation (18) in the main text.

Appendix C.7. Dimensional Verification

$$\begin{array}{cc}\hfill [(c_1+c_3){\mathcal{V}}_0^2{c}^2\mathcal{V}\nabla \mathcal{V}]& =\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]✓\hfill \end{array}$$

Appendix C.8. 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 (C3g) has a mass term that leads to exponential suppression on scales larger than $\lambda \approx 4\mathrm{pc}$. This explains why the field affects galaxy dynamics through its asymptotic value ${\mathcal{V}}_0$ while being undetectable in the Solar System.

3.

The coupling constant $\alpha =(c_1+c_3){\mathcal{V}}_0^2$ is determined by matching the asymptotic behavior of the full nonlinear solution on galactic scales.

4.

The force law ${\mathbf{a}}_{\mathrm{FST}}=-(c_1+c_3){\mathcal{V}}_0^2{c}^2\mathcal{V}\nabla \mathcal{V}$ is valid for the full nonlinear solution on galactic scales, where $\mathcal{V}$ is given by Equation (16a).

Appendix C.9. Consistency of the Sign Convention

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

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

Therefore $\nabla \mathcal{V}=({\mathcal{V}}_0/L_0){\nabla }_{\xi }\tilde{\mathcal{V}}<0$, and $\mathcal{V}\nabla \mathcal{V}<0$. The negative sign in Equation (C6) thus gives:

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

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.