Xuan-Liang Unified Field Theory: From Multi-Velocity Xuan-Liang Construction to Cosmology and Astrophysical Tests

1. Introduction

Modern physics is built upon two pillars: general relativity and quantum field theory. General relativity successfully describes the connection between gravity and spacetime geometry, while quantum field theory perfectly explains the electromagnetic, weak, and strong interactions. However, efforts to incorporate gravity into a quantum framework have not yet succeeded, and cosmological observations revealing dark matter and dark energy suggest shortcomings in existing theories.

This paper proposes a novel unified field theory framework—Xuan-Liang unified field theory. Starting from fundamental kinematic quantities, we introduce a new physical quantity, “Xuan-Liang”, defined as $X={\displaystyle \frac{1}{3}}m{v}^3$, filling the geometric sequence after mass, momentum, and kinetic energy. By geometrizing Xuan-Liang as a differential form and constructing a unified action incorporating curvature coupling and Xuan-Liang dynamics, we derive a set of self-consistent unified field equations.

The main contributions of this paper include:

1.

Proposing a hierarchical definition system for the Xuan-Liang concept, allowing flexible construction of differential forms for different physical contexts.

2.

Establishing a rigorous differential geometric definition, correcting mathematical issues in the original definition.

3.

Deriving the unified equations in detail, including both integral and differential forms, ensuring their equivalence.

4.

Systematically proving the reduction to classical theories (general relativity, Newtonian gravity, cosmology), ensuring consistency of the reduction.

5.

Presenting application examples in dark matter, dark energy, black hole physics, and the early universe.

6.

Providing testable theoretical predictions and an experimental roadmap.

2. Basic Definition and Geometric Meaning of Xuan-Liang

2.1. Algebraic Definition and Physical Origin of Xuan-Liang

For an object of mass m and velocity v, its Xuan-Liang X is defined as:

$$X={\displaystyle \frac{1}{3}}m{v}^3$$

(1)

Its dimension is $\left[M\right]{\left[L\right]}^3{\left[T\right]}^{-3}$, filling the geometric gap in the sequence of physical quantities. The term “Xuan-Liang” is derived from the Chinese classics I Ching (“heaven is dark, earth is yellow”) and Tao Te Ching (“the mystery of mysteries, the gate of all wonders”), signifying a profound and fundamental property of motion.

Remark 1.

The coefficient $1/3$ is a natural consequence of the geometric hierarchy. In the path integral of uniformly accelerated linear motion, this coefficient appears naturally, ensuring that Xuan-Liang and kinetic energy ${\displaystyle \frac{1}{2}}m{v}^2$ form a complete geometric sequence. Consider uniformly accelerated motion $v\left(t\right)=at$; the power $P=mav$, integrated along the path:

$$\int Pdl=\int mav·vdt=m{a}^2\int t^{2}dt={\displaystyle \frac{1}{3}}m{a}^3{t}^3={\displaystyle \frac{1}{3}}m{v}^3$$

2.2. Geometric Hierarchy: From Mass to Xuan-Liang

The physical quantities describing the state of motion in classical mechanics exhibit a clear geometric hierarchy:

Table 1. Geometric hierarchy of physical quantities

Order	Quantity	Expression	Dimension	Geometric interpretation
0	Mass	m	$\left[M\right]$	Point attribute: existence
1	Momentum	$m\overrightarrow{v}$	$\left[M\right]\left[L\right]{\left[T\right]}^{-1}$	Line attribute: directed motion
2	Kinetic energy	${\displaystyle \frac{1}{2}}m{v}^2$	$\left[M\right]{\left[L\right]}^2{\left[T\right]}^{-2}$	Surface attribute: intensity of motion
3	Xuan-Liang	${\displaystyle \frac{1}{3}}m{v}^3$	$\left[M\right]{\left[L\right]}^3{\left[T\right]}^{-3}$	Volume attribute: accumulation of energy flow

Table 1. Geometric hierarchy of physical quantities

Order	Quantity	Expression	Dimension	Geometric interpretation
0	Mass	m	$\left[M\right]$	Point attribute: existence
1	Momentum	$m\overrightarrow{v}$	$\left[M\right]\left[L\right]{\left[T\right]}^{-1}$	Line attribute: directed motion
2	Kinetic energy	${\displaystyle \frac{1}{2}}m{v}^2$	$\left[M\right]{\left[L\right]}^2{\left[T\right]}^{-2}$	Surface attribute: intensity of motion
3	Xuan-Liang	${\displaystyle \frac{1}{3}}m{v}^3$	$\left[M\right]{\left[L\right]}^3{\left[T\right]}^{-3}$	Volume attribute: accumulation of energy flow

This sequence reveals an intrinsic symmetry of physical quantities: m (scalar, order 0) →$mv$ (vector, order 1) →${\displaystyle \frac{1}{2}}m{v}^2$ (scalar, order 2) →${\displaystyle \frac{1}{3}}m{v}^3$ (pseudoscalar, order 3).

2.3. Path-Integral Definition of Xuan-Liang

In physics, the effect of a force on the motion of an object is manifested through accumulation over space—work:

$$W={\int }_C\overrightarrow{F}·d\overrightarrow{x}=\Delta E_k,E_k={\displaystyle \frac{1}{2}}m{v}^2$$

(2)

The instantaneous rate of energy change is power:

$$P={\displaystyle \frac{d{E}_k}{dt}}=\overrightarrow{F}·\overrightarrow{v}$$

(3)

Xuan-Liang is defined as the line integral of power along the path of motion:

Definition 1

(Xuan-Liang: path integral of power). Xuan-Liang X is defined as the line integral of power P along the motion path C:

$$\begin{array}{|c|}\hline X:={\oint }_{C}Pdl={\oint }_C{\displaystyle \frac{d{E}_k}{dt}}dl\\ \hline\end{array}$$

(4)

where $dl$ is the line element of the path. This definition gives X a clear geometric meaning: it measures the total “deposition” of the action that changes kinetic energy along the entire trajectory.

For uniformly accelerated linear motion with acceleration a and initial velocity zero, we have $v=at$, $P=m{a}^{2}t$, $dl=vdt=atdt$, and substituting into the integral gives:

$$X={\int }_0^{T}m{a}^{2}t·atdt=m{a}^3{\int }_0^T{t}^{2}dt={\displaystyle \frac{1}{3}}m{\left(aT\right)}^3={\displaystyle \frac{1}{3}}m{v}^3,$$

consistent with the algebraic definition.

3. Hierarchical Definition System of the Xuan-Liang Concept

In the history of physics, many concepts have evolved from fixed forms to context-dependent ones. For example:

• Energy: from mechanical energy to thermal, electrical, chemical, etc.
• Momentum: from classical momentum to relativistic momentum, quantum momentum operator.
• Field: from classical fields to quantum fields, from scalar fields to spinor fields, gauge fields.

Based on this historical perspective, we propose: The differential form of Xuan-Liang can be flexibly constructed according to the needs of specific physical problems. This flexibility is not a weakness of the theory but an advantage for adapting to different physical scales and contexts.

Definition 2

(Hierarchical system of the Xuan-Liang concept). The Xuan-Liang concept can be divided into the following levels according to physical context:

1.

Classical particle level: $X={\displaystyle \frac{1}{3}}m{v}^3$ (scalar)

2.

Continuum level: $\mathbf{X}={\displaystyle \frac{1}{3}}\rho v^3\mathbf{n}$ (vector field, $\mathbf{n}$ is direction)

3.

Multiple motion modes combination level: $X=\kappa m\left(v_1{v}_2{v}_3\right)$ (product of multiple velocities)

4.

Relativistic continuum level: $X^{\mu }={\displaystyle \frac{1}{3}}{\rho }_0{c}^2({\gamma }^3-\gamma )u^{\mu }$ (four-vector)

5.

Differential geometry level: $\mathbb{X}=\kappa {\rho }_0{\gamma }^3{v}^2★u$ (differential form)

6.

Quantum level: $\widehat{X}={\displaystyle \frac{1}{3}}m{\widehat{v}}^3$ (operator)

7.

Quantum field theory level: $\mathcal{X}\left(x\right)={\psi }^{†}{\gamma }^0\left({\displaystyle \frac{1}{3}}{v}^3\right)\psi$ (field operator)

Each level corresponds to different physical situations and approximations, from simple to complex, from concrete to abstract.

3.1. Multi-Velocity Combination Definition of Xuan-Liang

In many physical systems, an object participates in multiple modes of motion simultaneously. For example:

• Electron in an atom: simultaneously has

  1.

  Translational velocity of the whole atom $v_{\mathrm{trans}}$

  2.

  Orbital velocity around the nucleus $v_{\mathrm{orbit}}$

  3.

  Spin motion (equivalent velocity) $v_{\mathrm{spin}}$

• Earth’s motion: simultaneously has

  1.

  Rotational velocity $v_{\mathrm{rot}}$

  2.

  Revolutionary velocity $v_{\mathrm{rev}}$

  3.

  Velocity of the Solar System around the Galactic center $v_{\mathrm{gal}}$

For such systems, Xuan-Liang can be naturally extended to the form of a product of multiple velocities:

Definition 3

(Multi-velocity Xuan-Liang). For a system with n independent motion modes, Xuan-Liang is defined as:

$$X=\kappa m\prod _{i=1}^n{v}_i^{{\alpha }_i}$$

(5)

where:

• $v_i$: characteristic velocity of the i-th motion mode
• ${\alpha }_i$: exponent of the i-th velocity, satisfying ${\sum }_{i=1}^n{\alpha }_i=3$
• κ: normalization coefficient determined by the specific geometric structure of motion

3.1.1. Special Case Analysis

Example 1

(Equal-weight product of three velocities). When the system has three equally weighted motion modes, take ${\alpha }_i=1$, then:

$$X=\kappa m\left(v_1{v}_2{v}_3\right)$$

(6)

For the electron:

$$X_e={\kappa }_e{m}_e(v_{\mathrm{trans}}·v_{\mathrm{orbit}}·v_{\mathrm{spin}})$$

For the Earth:

$$X_{\oplus }={\kappa }_{\oplus }{M}_{\oplus }(v_{\mathrm{rot}}·v_{\mathrm{rev}}·v_{\mathrm{gal}})$$

Example 2

(Significance of coefficient variation). The variation of κ reflects geometric features in different physical contexts:

1.

Classical uniformly accelerated motion: $\kappa ={\displaystyle \frac{1}{3}}$, from path integral

2.

Rotational motion: $\kappa ={\displaystyle \frac{2}{3\pi }}$, from geometric factor of circular motion

3.

Quantum systems: $\kappa ={\displaystyle \frac{\hslash }{m{\lambda }_c}}$, introducing Compton wavelength ${\lambda }_c$

4.

Cosmological scales: $\kappa ={\displaystyle \frac{H_0^2}{c^2}}$, related to the Hubble constant

3.2. Physical Interpretation

The definition of multi-velocity Xuan-Liang has profound physical implications:

1.

Coupling of motion modes: The product of different velocities reflects interrelations between motion modes.

2.

Hierarchical structure: $v_1$ (micro) ×$v_2$ (meso) ×$v_3$ (macro) embodies cross-scale coupling.

3.

Complexity of energy transfer: Xuan-Liang, as the “acceleration of energy change”, more comprehensively describes the complex mechanisms of energy transfer in multi-mode systems.

4.

Correspondence with “Interpretation of the Essence of Energy in Xuan-Liang Theory”: This definition exactly aligns with the insight that “explosion is not energy release, but a reorganization of the Xuan-Liang field”.

3.3. Reduction to Classical Case

When the system has only one dominant motion mode, multi-velocity Xuan-Liang reduces to the classical form:

Proposition 1

(Classical reduction). Let $v_1=v$ be the dominant velocity, and $v_2,v_3,\dots$ be small. Then:

$$X=\kappa m{v}^3\prod _{i=2}^n{\left(1+{\displaystyle \frac{v_i}{v}}\right)}^{{\alpha }_i}\approx \kappa m{v}^3\left(1+\sum _{i=2}^n{\alpha }_i{\displaystyle \frac{v_i}{v}}\right)$$

When $v_i\ll v$, ${lim}_{v_i/v\to 0}X=\kappa m{v}^3$.

In particular, when $\kappa ={\displaystyle \frac{1}{3}}$, we recover the classical Xuan-Liang $X={\displaystyle \frac{1}{3}}m{v}^3$.

3.4. Generalization to Differential Forms

Multi-velocity Xuan-Liang can also be generalized to the differential geometry level. Suppose the system has n velocity 1-forms $u_1,u_2,\dots ,u_n$, then the Xuan-Liang differential form is defined as:

Definition 4

(Multi-velocity Xuan-Liang differential form).

$$\mathbb{X}=\kappa {\rho }_0\prod _{i=1}^n{\left({\gamma }_i{v}_i\right)}^{{\alpha }_i}★\left(\underset{i=1}{\overset{n}{\bigwedge }}{u}_i\right)$$

(7)

where:

• ⋀ denotes the exterior product
• The total form degree is n (when $\sum {\alpha }_i=3$, the coefficients need to be adjusted appropriately to maintain dimensional consistency)
• ★ is the Hodge star operator

This definition remains covariant in the subsequent unified equations because the exterior product ${\bigwedge }_{i=1}^n{u}_i$ is an n-form, and after Hodge duality becomes a $(4-n)$-form, which can be naturally incorporated into the action principle.

Remark 2

(Compatibility with unified equations). Although the multi-velocity definition is more complex in form, it does not 破坏 the mathematical structure of the unified equations:

1.

The action principle still applies: $\mathbb{X}$ as a differential form can still construct the term $\mathbb{X}\wedge ★\mathbb{X}$

2.

The variational principle still holds: one only needs to consider variations of multiple velocity fields

3.

Conservation laws remain: multiple motion modes correspond to multiple conservation equations

4.

Classical reduction remains unchanged: in appropriate limits, the Einstein field equations are recovered

Therefore, this extension enhances the universality of the theory without affecting its self-consistency.

3.4.1. Application Examples

Example 3

(Xuan-Liang of an atomic system). For an electron in a hydrogen atom:

$$\begin{array}{cc}\hfill v_{\mathrm{orbit}}& =\alpha c\approx 2.2\times {10}^6\mathrm{m}/\mathrm{s}\hfill \\ \hfill v_{\mathrm{spin}}& \sim {\displaystyle \frac{\hslash }{2m_e{a}_0}}\approx 1.16\times {10}^6\mathrm{m}/\mathrm{s}\hfill \\ \hfill v_{\mathrm{trans}}& (\mathrm{thermal}\mathrm{velocity}\mathrm{of}\mathrm{atom})\sim {10}^3\mathrm{m}/\mathrm{s}\hfill \end{array}$$

Electron Xuan-Liang:

$$X_e={\kappa }_e{m}_e(v_{\mathrm{orbit}}·v_{\mathrm{spin}}·v_{\mathrm{trans}})$$

where ${\kappa }_e$ can be determined from the hydrogen atom wave function.

Example 4

(Earth in the Solar System). Earth’s three characteristic velocities:

$$\begin{array}{cc}\hfill v_{\mathrm{rot}}& \approx 465\mathrm{m}/\mathrm{s}(\mathrm{equatorial}\mathrm{rotational}\mathrm{speed})\hfill \\ \hfill v_{\mathrm{rev}}& \approx 29.8\mathrm{km}/\mathrm{s}(\mathrm{orbital}\mathrm{speed})\hfill \\ \hfill v_{\mathrm{gal}}& \approx 230\mathrm{km}/\mathrm{s}(\mathrm{speed}\mathrm{around}\mathrm{Galactic}\mathrm{center})\hfill \end{array}$$

Earth Xuan-Liang:

$$X_{\oplus }={\kappa }_{\oplus }{M}_{\oplus }(v_{\mathrm{rot}}·v_{\mathrm{rev}}·v_{\mathrm{gal}})\approx 2.4\times {10}^{32}{\mathrm{J}\text{·}\mathrm{m}}^3/{\mathrm{s}}^3(\mathrm{taking}{\kappa }_{\oplus }\approx 0.1)$$

Figure 1. Conceptual extension from single-velocity to multi-velocity Xuan-Liang and corresponding physical systems

Figure 1. Conceptual extension from single-velocity to multi-velocity Xuan-Liang and corresponding physical systems

3.5. Deep Connection Between Xuan-Liang and Energy

We can further elucidate the physical meaning of multi-velocity Xuan-Liang:

Proposition 2

(Xuan-Liang as the deep structure of energy). Xuan-Liang reveals the essential structure of energy:

1.

Energy is not a fundamental entity: Energy is an外在 manifestation of the Xuan-Liang field.

2.

Essence of energy transfer: It is the propagation of Xuan-Liang flow in spacetime.

3.

Physical meaning of multi-velocity: It reflects the complexity and multi-channel nature of energy transfer.

4.

Philosophical connotation of coefficient variation: The variation of κ embodies the “mystery of mysteries, the gate of all wonders”—the diversity of energy manifestations in different contexts.

Example 5

(Reinterpretation of an explosion process). Consider the explosion process described in “Xuan-Liang and Energy”:

• Traditional description: chemical energy → light energy + heat energy + kinetic energy
• Xuan-Liang description: stored Xuan-Liang → radiation Xuan-Liang + thermal Xuan-Liang + kinetic Xuan-Liang
• Significance of multi-velocity: Fragments produced in an explosion have complex motion combinations, corresponding to multi-velocity Xuan-Liang.

An explosion is not a simple release of energy, but a 剧烈 reorganization and redistribution of the Xuan-Liang field.

This extension not only enriches the connotation of Xuan-Liang theory but also provides a deeper mathematical framework for understanding the essence of energy.

4. Rigorous Definition and Flexible Construction of Xuan-Liang Differential Forms

4.1. Revised Definition of Xuan-Liang Differential Form

Definition 5

(Revised Xuan-Liang differential form). On a four-dimensional spacetime manifold $\mathcal{M}$, given:

• Rest mass density field ${\rho }_0$
• Velocity 1-form u, satisfying $g_{\mu \nu }{u}^{\mu }{u}^{\nu }=-c^2$
• Lorentz factor $\gamma ={(1-v^2/c^2)}^{-1/2}$
• Spatial velocity squared $v^2=c^2(1-1/{\gamma }^2)$

Define the Xuan-Liang differential form as:

$$\mathbb{X}=\kappa {\rho }_0{\gamma }^3{v}^2★u$$

(8)

where κ is a constant to be determined, and $★u$ is the Hodge dual of u (a 3-form).

4.2. Vector Representation of Xuan-Liang Current

We can also define the Xuan-Liang current vector:

$$J_X^{\mu }={\displaystyle \frac{1}{3}}{\rho }_0{c}^2({\gamma }^3-\gamma )u^{\mu }$$

(9)

The corresponding Xuan-Liang differential form is its Hodge dual:

$$\mathbb{X}=★J_X$$

(10)

5. Action Principle of Xuan-Liang Unified Field Theory

Based on the least action principle and covariance, we construct the following unified action:

Definition 6

(Xuan-Liang unified field theory action).

$$S={\int }_{\mathcal{M}}\left[{\displaystyle \frac{1}{16\pi G}}R+{\mathcal{L}}_{\mathrm{matter}}+\alpha R{J}_X^{\mu }{u}_{\mu }+{\displaystyle \frac{\beta }{2}}\mathbb{X}\wedge ★\mathbb{X}\right]\sqrt{-g}{d}^{4}x$$

(11)

where:

• ${\displaystyle \frac{1}{16\pi G}}R$: Einstein-Hilbert term, describing gravitational dynamics
• ${\mathcal{L}}_{\mathrm{matter}}$: ordinary matter Lagrangian
• $\alpha R{J}_X^{\mu }{u}_{\mu }$: curvature coupling term, describing the interaction between the Xuan-Liang field and spacetime geometry
• ${\displaystyle \frac{\beta }{2}}\mathbb{X}\wedge ★\mathbb{X}$: kinetic term of the Xuan-Liang field, describing its dynamics

α and β are dimensionless coupling constants.

6. Detailed Derivation of the Unified Equations

6.1. Application of the Variational Principle

We derive the equations of motion by varying the action S. The variational principle requires:

$$\delta S=\delta S_{\mathrm{EH}}+\delta S_{\mathrm{matter}}+\delta S_{\mathrm{coupling}}+\delta S_{\mathrm{kinetic}}=0$$

6.2. Detailed Calculation of Metric Variation

6.2.1. Variation of the Einstein-Hilbert Term

The standard result is:

$$\delta \int R\sqrt{-g}{d}^{4}x=\int \left(R_{\mu \nu }-{\displaystyle \frac{1}{2}}R{g}_{\mu \nu }\right)\delta g^{\mu \nu }\sqrt{-g}{d}^{4}x$$

(12)

6.2.2. Variation of the Matter Term

Define:

$$\delta \int {\mathcal{L}}_{\mathrm{matter}}\sqrt{-g}{d}^{4}x={\displaystyle \frac{1}{2}}\int T_{\mu \nu }^{\left(\mathrm{matter}\right)}\delta g^{\mu \nu }\sqrt{-g}{d}^{4}x$$

(13)

6.2.3. Variation of the Xuan-Liang Field Terms

The Xuan-Liang field terms consist of two parts:

1. Variation of the curvature coupling term:

$$\begin{array}{cc}\hfill \delta \int \alpha R{J}_X^{\mu }{u}_{\mu }\sqrt{-g}{d}^{4}x& =\alpha \int \left[\delta R·J_X^{\mu }{u}_{\mu }+R·\delta \left(J_X^{\mu }{u}_{\mu }\right)\right]\sqrt{-g}{d}^{4}x\hfill \end{array}$$

(14)

$$\begin{array}{cc}\hfill & =\alpha \int \left[\left(R_{\mu \nu }-{\nabla }_{\mu }{\nabla }_{\nu }+g_{\mu \nu }\square \right)\delta g^{\mu \nu }·J_X^{\sigma }{u}_{\sigma }\right.\hfill \end{array}$$

(15)

$$\begin{array}{cc}\hfill & \left.+R·{\displaystyle \frac{\partial \left(J_X^{\sigma }{u}_{\sigma }\right)}{\partial g^{\mu \nu }}}\delta g^{\mu \nu }\right]\sqrt{-g}{d}^{4}x\hfill \end{array}$$

(16)

2. Variation of the Xuan-Liang kinetic term:

$$\begin{array}{cc}\hfill \delta \int {\displaystyle \frac{\beta }{2}}\mathbb{X}\wedge ★\mathbb{X}& ={\displaystyle \frac{\beta }{2}}\int \left[\delta \mathbb{X}\wedge ★\mathbb{X}+\mathbb{X}\wedge \delta (★\mathbb{X})\right]\hfill \end{array}$$

(17)

$$\begin{array}{cc}\hfill & ={\displaystyle \frac{\beta }{2}}\int T_{\mu \nu }^{\left(\mathrm{kinetic}\right)}\delta g^{\mu \nu }\sqrt{-g}{d}^{4}x+\mathrm{boundary}\mathrm{terms}\hfill \end{array}$$

(18)

where:

$$T_{\mu \nu }^{\left(\mathrm{kinetic}\right)}=\beta \left({\mathbb{X}}_{\mu \alpha \beta }{\mathbb{X}}_{\nu }^{\alpha \beta }-{\displaystyle \frac{1}{6}}{g}_{\mu \nu }{\mathbb{X}}_{\alpha \beta \gamma }{\mathbb{X}}^{\alpha \beta \gamma }\right)$$

6.3. Differential Form of the Unified Equations

Combining all variational terms, we obtain:

Theorem 1

(Xuan-Liang unified differential equations). The fundamental equations of Xuan-Liang unified field theory are the following differential system:

$$\begin{array}{|c|}\hline \begin{array}{cc}\hfill 1.\mathit{Gravitational}\mathit{field}\mathit{equations}\text{:}& G_{\mu \nu }=8\pi G\left(T_{\mu \nu }^{\left(\mathbf{matter}\right)}+T_{\mu \nu }^{\left(\mathrm{X}\right)}\right)\hfill \\ \hfill 2.\mathit{Xuan}-\mathit{Liang}\mathit{field}\mathit{equations}\text{:}& d★d\mathbb{X}+\alpha d(R★\mathcal{U})=J_X\hfill \\ \hfill 3.\mathit{Conservation}\mathit{equations}\text{:}& {\nabla }_{\mu }{T}_{\left(\mathbf{total}\right)}^{\mu \nu }=0\hfill \end{array}\\ \hline\end{array}$$

(19)

where:

• $T_{\mu \nu }^{\left(\mathrm{X}\right)}=\alpha T_{\mu \nu }^{\left(\mathrm{coupling}\right)}+\beta T_{\mu \nu }^{\left(\mathrm{kinetic}\right)}$ is the total energy-momentum tensor of the Xuan-Liang field
• $\mathcal{U}={\displaystyle \frac{1}{2}}u\wedge ★u$ is the 2-form of the observer’s velocity
• $J_X$ is the source term of the Xuan-Liang current

6.4. Integral Form of the Unified Equations

By Stokes’ theorem and Hodge duality, we obtain the integral form of the unified equations:

Theorem 2

(Xuan-Liang unified integral equations). The differential system (19) is equivalent to the following integral form:

$$\begin{array}{|c|}\hline \begin{array}{cc}\hfill {\int }_{\partial \mathcal{V}}★G& =8\pi G{\int }_{\mathcal{V}}★\left(T^{\left(\mathrm{matter}\right)}+T^{\left(\mathrm{X}\right)}\right)\hfill \\ \hfill {\int }_{\partial \mathcal{S}}★d\mathbb{X}+\alpha {\int }_{\partial \mathcal{S}}R★\mathcal{U}& ={\int }_{\mathcal{S}}★J_X\hfill \\ \hfill {\int }_{\partial \mathcal{V}}★T^{\left(\mathrm{total}\right)}& =0\hfill \end{array}\\ \hline\end{array}$$

(20)

where $\mathcal{V}$ is an arbitrary four-dimensional region, $\mathcal{S}$ is an arbitrary three-dimensional hypersurface, and ∂ denotes the boundary.

Equivalence of differential and integral forms 

Consider the first integral equation of the gravitational field. By the generalized Stokes theorem:

$${\int }_{\partial \mathcal{V}}★G={\int }_{\mathcal{V}}d★G$$

In differential form, the Einstein equation can be written as $G=8\pi G(T^{\left(\mathrm{matter}\right)}+T^{\left(\mathrm{X}\right)})$. Taking the Hodge dual and integrating:

$${\int }_{\mathcal{V}}d★G=8\pi G{\int }_{\mathcal{V}}★(T^{\left(\mathrm{matter}\right)}+T^{\left(\mathrm{X}\right)})$$

The equivalence of the other equations can be proved similarly. The advantage of the integral form is that it naturally incorporates boundary conditions and is more suitable for spacetimes with singularities (such as black holes).    □

7. Detailed Reduction of the Unified Equations to Classical Theories

7.1. Detailed Reduction to General Relativity

Theorem 3

(General relativity limit). Under the following limiting conditions:

1.

Weak field approximation: $g_{\mu \nu }={\eta }_{\mu \nu }+h_{\mu \nu }$, $|h_{\mu \nu }|\ll 1$

2.

Low velocity limit: $v\ll c$

3.

Small coupling approximation: $\left|\alpha \right|\ll 1$, $\left|\beta \right|\ll 1$

The Xuan-Liang unified equations reduce to the Einstein field equations:

$$G_{\mu \nu }=8\pi G{T}_{\mu \nu }^{\left(\mathrm{matter}\right)}$$

(21)

Detailed derivation 

Consider the following limits:

1.

Low velocity limit: $v/c\ll 1$, then $\gamma \approx 1+v^2/\left(2c^2\right)$, ${\gamma }^3-\gamma \approx v^2/c^2$

2.

Weak field limit: $g_{\mu \nu }={\eta }_{\mu \nu }+h_{\mu \nu }$, $|h_{\mu \nu }|\ll 1$

3.

Small coupling limit: $\left|\alpha \right|,\left|\beta \right|\ll 1$

In this limit, calculate the orders of magnitude of the Xuan-Liang terms:

$$\begin{array}{cc}\hfill J_X^{\mu }{u}_{\mu }& =-{\displaystyle \frac{1}{3}}{\rho }_0{c}^2({\gamma }^3-\gamma )\approx -{\displaystyle \frac{1}{3}}{\rho }_0{v}^2\hfill \end{array}$$

(22)

$$\begin{array}{cc}\hfill R& \sim {\partial }^{2}h\sim {\displaystyle \frac{v^2}{c^2{L}^2}}(\mathrm{typical}\mathrm{astrophysical}\mathrm{scale})\hfill \end{array}$$

(23)

$$\begin{array}{cc}\hfill \mathbb{X}\wedge ★\mathbb{X}& \sim {\rho }_0^2{v}^4\hfill \end{array}$$

(24)

Therefore, the order of the curvature coupling term is:

$$\alpha R{J}_X^{\mu }{u}_{\mu }\sim \alpha {\rho }_0{\displaystyle \frac{v^4}{c^2{L}^2}}$$

The order of the Xuan-Liang kinetic term is:

$${\displaystyle \frac{\beta }{2}}\mathbb{X}\wedge ★\mathbb{X}\sim \beta {\rho }_0^2{v}^4$$

Comparing with the Einstein-Hilbert term $R/\left(16\pi G\right)\sim v^2/(c^2{L}^2·16\pi G)$, for Solar System scales:

$$\begin{array}{cc}\hfill {\displaystyle \frac{S_{\mathrm{coupling}}}{S_{\mathrm{EH}}}}& \sim 16\pi G\alpha {\rho }_0{L}^2{\displaystyle \frac{v^2}{c^2}}\hfill \end{array}$$

(25)

$$\begin{array}{cc}\hfill & \sim \alpha \times \left({10}^{-7}\right)(\mathrm{substituting}\mathrm{Solar}\mathrm{System}\mathrm{parameters})\hfill \end{array}$$

(26)

When $\left|\alpha \right|,\left|\beta \right|\ll 1$, the Xuan-Liang terms are much smaller than the Einstein-Hilbert and matter terms and can be neglected. Thus the unified equations reduce to the Einstein field equations.    □

7.2. Determination of Parameters and Physical Interpretation

In the classical reduction, special attention must be paid to the determination of the coupling constants $\alpha$ and $\beta$:

Corollary 1

(Physical constraints on coupling constants). From the relations in the general relativity limit, we obtain constraints on the coupling constants:

$${\displaystyle \frac{\alpha }{16\pi G\kappa {\rho }_0}}\ll 1,\beta \ll {\displaystyle \frac{1}{{\rho }_0^2{v}^4}}$$

(27)

where ${\rho }_0$ is the typical mass density of the astrophysical system, and v is the characteristic velocity.

Proof. 

From the order-of-magnitude comparison in the reduction proof:

$$16\pi G\alpha {\rho }_0{L}^2{\displaystyle \frac{v^2}{c^2}}\ll 1$$

For Solar System scales, $L\sim 1\mathrm{AU}$, ${\rho }_0\sim M_{\odot }/L^3$, $v^2\sim G{M}_{\odot }/L$, substituting gives:

$$\alpha \ll {10}^7$$

Similar reasoning gives the constraint on $\beta$.    □

7.3. Detailed Reduction to Newtonian Gravity

Theorem 4

(Newtonian limit). Under the stronger conditions:

1.

Static field: all time derivatives vanish

2.

Weak field: $|h_{\mu \nu }|\ll 1$

3.

Low velocity: $v\ll c$

4.

Small pressure: $p\ll \rho c^2$

Xuan-Liang theory reduces to Newtonian gravity:

$${\nabla }^2\mathrm{\Phi }=4\pi G\rho$$

(28)

Detailed derivation 

Consider the following conditions:

1.

Static field: ${\partial }_t{g}_{\mu \nu }=0$

2.

Weak field: $g_{\mu \nu }={\eta }_{\mu \nu }+h_{\mu \nu }$, $|h_{\mu \nu }|\ll 1$

3.

Low velocity: $v/c\ll 1$

4.

Small pressure: $p\ll \rho c^2$

Take the metric:

$$d{s}^2=-(1+2\mathrm{\Phi }/c^2)c^{2}d{t}^2+(1-2\mathrm{\Phi }/c^2){\delta }_{ij}d{x}^{i}d{x}^j$$

where $\mathrm{\Phi }$ is the Newtonian potential.

In the Newtonian limit, the Xuan-Liang terms can be neglected (for the same reasons as in the general relativity limit). Compute the time-time component of the Einstein tensor:

$$\begin{array}{cc}\hfill G_{00}& =R_{00}-{\displaystyle \frac{1}{2}}R{g}_{00}\hfill \end{array}$$

(29)

$$\begin{array}{cc}\hfill & \approx {\displaystyle \frac{1}{c^2}}{\nabla }^2\mathrm{\Phi }-{\displaystyle \frac{1}{2}}\left(-{\displaystyle \frac{2}{c^2}}{\nabla }^2\mathrm{\Phi }\right)(-1)\hfill \end{array}$$

(30)

$$\begin{array}{cc}\hfill & =-{\displaystyle \frac{2}{c^2}}{\nabla }^2\mathrm{\Phi }\hfill \end{array}$$

(31)

The 00-component of the matter energy-momentum tensor: $T_{00}=\rho c^2$

The Einstein equation $G_{00}=8\pi G{T}_{00}$ gives:

$$-{\displaystyle \frac{2}{c^2}}{\nabla }^2\mathrm{\Phi }=8\pi G\rho c^2\Rightarrow {\nabla }^2\mathrm{\Phi }=-4\pi G\rho$$

Adjusting the sign convention yields the Newtonian equation ${\nabla }^2\mathrm{\Phi }=4\pi G\rho$.    □

7.4. Detailed Reduction to Cosmological Equations

7.4.1. Basic Assumptions for the Cosmological Limit

Theorem 5

(Cosmological limit). In a homogeneous and isotropic universe (FRW metric), the Xuan-Liang unified equations reduce to modified Friedmann equations:

$$\begin{array}{cc}\hfill H^2& ={\displaystyle \frac{8\pi G}{3}}\left({\rho }_m+{\rho }_X\right)\hfill \end{array}$$

(32)

$$\begin{array}{cc}\hfill {\displaystyle \frac{\ddot{a}}{a}}& =-{\displaystyle \frac{4\pi G}{3}}\left({\rho }_m+{\rho }_X+3p_X\right)\hfill \end{array}$$

(33)

where the Xuan-Liang field energy density ${\rho }_X$ and pressure $p_X$ are given by:

$$\begin{array}{cc}\hfill {\rho }_X& ={\displaystyle \frac{\beta }{2}}\mathbb{X}\wedge ★\mathbb{X}+\alpha R{J}_X^{\mu }{u}_{\mu }\hfill \end{array}$$

(34)

$$\begin{array}{cc}\hfill p_X& ={\displaystyle \frac{\beta }{6}}\mathbb{X}\wedge ★\mathbb{X}-\alpha R{J}_X^{\mu }{u}_{\mu }\hfill \end{array}$$

(35)

Detailed derivation 

Consider the flat FRW metric:

$$d{s}^2=-d{t}^2+a^2\left(t\right)\left[d{r}^2+r^2(d{\theta }^2+{sin}^2\theta d{\varphi }^2)\right]$$

Compute the geometric quantities:

$$\begin{array}{cc}\hfill {\mathrm{\Gamma }}_{ij}^0& =a\dot{a}{\delta }_{ij}\hfill \end{array}$$

(36)

$$\begin{array}{cc}\hfill {\mathrm{\Gamma }}_{0j}^i& ={\displaystyle \frac{\dot{a}}{a}}{\delta }_j^i\hfill \end{array}$$

(37)

$$\begin{array}{cc}\hfill R_{00}& =-3{\displaystyle \frac{\ddot{a}}{a}}\hfill \end{array}$$

(38)

$$\begin{array}{cc}\hfill R_{ij}& =(a\ddot{a}+2{\dot{a}}^2){\delta }_{ij}\hfill \end{array}$$

(39)

$$\begin{array}{cc}\hfill R& =6\left({\displaystyle \frac{\ddot{a}}{a}}+{\displaystyle \frac{{\dot{a}}^2}{a^2}}\right)=6(\dot{H}+2H^2)\hfill \end{array}$$

(40)

where $H=\dot{a}/a$ is the Hubble parameter.

On cosmological scales, we assume the Xuan-Liang field is a homogeneous and isotropic background field. For comoving observers, $u=-dt$. Substituting into the Xuan-Liang field definition (8):

$$\mathbb{X}=\kappa {\rho }_0{\gamma }^3{v}^2★u$$

In cosmology, the natural “velocity” is the Hubble flow velocity: $v=Hr$. However, for a homogeneous background field, we need to consider a statistical average. Define the effective velocity squared:

$$\langle v^2\rangle =\xi H^2$$

where $\xi$ is a constant determined by fitting observational data.

Compute the components of the Xuan-Liang field energy-momentum tensor:

$$\begin{array}{cc}\hfill T_{00}^{\left(\mathrm{X}\right)}& ={\rho }_X={\displaystyle \frac{\beta }{2}}\mathbb{X}\wedge ★\mathbb{X}+\alpha R{J}_X^{\mu }{u}_{\mu }\hfill \end{array}$$

(41)

$$\begin{array}{cc}\hfill T_{ij}^{\left(\mathrm{X}\right)}& =p_X{a}^2{\delta }_{ij}=\left({\displaystyle \frac{\beta }{6}}\mathbb{X}\wedge ★\mathbb{X}-\alpha R{J}_X^{\mu }{u}_{\mu }\right)a^2{\delta }_{ij}\hfill \end{array}$$

(42)

Explicit calculation of the terms:

$$\begin{array}{cc}\hfill \mathbb{X}\wedge ★\mathbb{X}& ={\kappa }^2{\rho }_0^2{\gamma }^6{v}^4(★u)\wedge (★★u)\hfill \end{array}$$

(43)

$$\begin{array}{cc}\hfill & =-{\kappa }^2{\rho }_0^2{\gamma }^6{v}^4\sqrt{-g}{d}^{4}x\hfill \end{array}$$

(44)

$$\begin{array}{cc}\hfill J_X^{\mu }{u}_{\mu }& =-{\displaystyle \frac{1}{3}}{\rho }_0{c}^2({\gamma }^3-\gamma )\hfill \end{array}$$

(45)

In the cosmological background, ${\rho }_0$ is the current rest mass density, related to the scale factor by ${\rho }_0\propto a^{-3}$. Substituting these expressions into the unified equations yields the modified Friedmann equations (32) and ().

To obtain the explicit evolution equations, we need to determine the equation of state $w_X=p_X/{\rho }_X$ of the Xuan-Liang field. From (34) and () we get:

$$w_X={\displaystyle \frac{p_X}{{\rho }_X}}={\displaystyle \frac{\frac{\beta }{6}\mathbb{X}\wedge ★\mathbb{X}-\alpha R{J}_X^{\mu }{u}_{\mu }}{\frac{\beta }{2}\mathbb{X}\wedge ★\mathbb{X}+\alpha R{J}_X^{\mu }{u}_{\mu }}}$$

Define the parameter:

$$\eta ={\displaystyle \frac{\alpha R{J}_X^{\mu }{u}_{\mu }}{\frac{\beta }{2}\mathbb{X}\wedge ★\mathbb{X}}}$$

Then:

$$w_X={\displaystyle \frac{1-3\eta }{3(1+\eta )}}$$

In the early universe (high curvature, large R), $\eta \gg 1$, $w_X\approx -1$, behaving like a cosmological constant; In the late universe (low curvature, small R), $\eta \ll 1$, $w_X\approx 1/3$, behaving like radiation.

This provides a natural transition mechanism from early accelerated expansion to later decelerated expansion.    □

7.4.2. Supplementary Derivation of the Cosmological Phase Transition Equation

To more naturally describe the phase transition of the Xuan-Liang field from matter-like to cosmological-constant-like behavior, we adopt a dynamical phase transition equation-of-state parametrization:

Definition 7

(Dynamical phase transition equation of state). The equation-of-state parameter of the Xuan-Liang field is:

$$\begin{array}{|c|}\hline w\left({\rho }_X\right)=-1+{\displaystyle \frac{1}{2}}\left[1+tanh\left({\displaystyle \frac{ln({\rho }_X/{\rho }_t)}{\Delta }}\right)\right]\\ \hline\end{array}$$

(46)

where:

• ${\rho }_t$:Phase transition critical density. When ${\rho }_X={\rho }_t$, $w=-0.5$
• Δ:Phase transition width(dimensionless), controlling the smoothness of the transition
• Asymptotic behavior:

  $${\rho }_X\gg {\rho }_t\Rightarrow w\to 0(\mathrm{matter}-\mathrm{like});{\rho }_X\ll {\rho }_t\Rightarrow w\to -1(\cos \mathrm{mological}-\mathrm{constant}-\mathrm{like})$$

Theorem 6

(Cosmological phase transition equation). In a homogeneous and isotropic universe, the Xuan-Liang field energy density satisfies the phase transition equation:

$${\left({\displaystyle \frac{{\rho }_X}{{\rho }_t}}\right)}^{\Delta /2}+{\left({\displaystyle \frac{{\rho }_X}{{\rho }_t}}\right)}^{-\Delta /2}={\left({\displaystyle \frac{a}{a_t}}\right)}^{-3\Delta /2}+{\left({\displaystyle \frac{a}{a_t}}\right)}^{3\Delta /2}$$

(47)

where ${\rho }_t$ is the phase transition critical density, $a_t$ is the scale factor at the transition, and Δ controls the transition width.

Proof. 

In the FRW metric, the Xuan-Liang field satisfies the continuity equation:

$${\displaystyle \frac{d{\rho }_X}{dt}}+3H({\rho }_X+p_X)=0$$

Substituting the equation of state $p_X=w\left({\rho }_X\right){\rho }_X$, where $w\left({\rho }_X\right)$ is given by equation (46):

$${\displaystyle \frac{d{\rho }_X}{{\rho }_X}}=-3[1+w\left({\rho }_X\right)]{\displaystyle \frac{da}{a}}$$

Insert the expression for $w\left({\rho }_X\right)$:

$$\begin{array}{cc}\hfill {\displaystyle \frac{d{\rho }_X}{{\rho }_X}}& =-{\displaystyle \frac{3}{2}}\left[1+tanh\left({\displaystyle \frac{ln({\rho }_X/{\rho }_t)}{\Delta }}\right)\right]{\displaystyle \frac{da}{a}}\hfill \end{array}$$

(48)

$$\begin{array}{cc}\hfill & =-{\displaystyle \frac{3}{2}}\left[{\displaystyle \frac{2}{1+e^{-2ln({\rho }_X/{\rho }_t)/\Delta }}}\right]{\displaystyle \frac{da}{a}}\hfill \end{array}$$

(49)

$$\begin{array}{cc}\hfill & =-{\displaystyle \frac{3}{1+{({\rho }_t/{\rho }_X)}^{2/\Delta }}}{\displaystyle \frac{da}{a}}\hfill \end{array}$$

(50)

Let $R={\rho }_X/{\rho }_t$, then:

$${\displaystyle \frac{dR}{R}}=-{\displaystyle \frac{3}{1+R^{-2/\Delta }}}{\displaystyle \frac{da}{a}}$$

Separate variables:

$${\displaystyle \frac{1+R^{-2/\Delta }}{R}}dR=-3{\displaystyle \frac{da}{a}}$$

Integrate:

$$lnR+{\displaystyle \frac{\Delta }{2}}{R}^{-2/\Delta }=-3lna+\mathrm{constant}$$

Let $a=a_t$ when $R=1$, then the constant is:

$$\mathrm{constant}={\displaystyle \frac{\Delta }{2}}+3ln{a}_t$$

Substitute and rearrange:

$$lnR+{\displaystyle \frac{\Delta }{2}}(R^{-2/\Delta }-1)=-3ln(a/a_t)$$

Let $y=R^{1/\Delta }$, then:

$$\Delta lny+{\displaystyle \frac{\Delta }{2}}(y^{-2}-1)=-3ln(a/a_t)$$

Divide both sides by $\Delta /2$:

$$2lny+y^{-2}-1=-{\displaystyle \frac{6}{\Delta }}ln(a/a_t)$$

Note that $2lny=ln\left(y^2\right)$, and $y^{-2}=1/y^2$, so:

$$ln\left(y^2\right)+{\displaystyle \frac{1}{y^2}}-1=-{\displaystyle \frac{6}{\Delta }}ln(a/a_t)$$

Let $z=y^2=R^{2/\Delta }$, then the equation becomes:

$$lnz+{\displaystyle \frac{1}{z}}-1=-{\displaystyle \frac{6}{\Delta }}ln(a/a_t)$$

This is an implicit equation. Through appropriate algebraic manipulation, we obtain the symmetric form:

$$z^{1/2}+z^{-1/2}={\left({\displaystyle \frac{a}{a_t}}\right)}^{-3/\Delta }+{\left({\displaystyle \frac{a}{a_t}}\right)}^{3/\Delta }$$

Substituting back $z={({\rho }_X/{\rho }_t)}^{2/\Delta }$ yields:

$${\left({\displaystyle \frac{{\rho }_X}{{\rho }_t}}\right)}^{\Delta /2}+{\left({\displaystyle \frac{{\rho }_X}{{\rho }_t}}\right)}^{-\Delta /2}={\left({\displaystyle \frac{a}{a_t}}\right)}^{-3\Delta /2}+{\left({\displaystyle \frac{a}{a_t}}\right)}^{3\Delta /2}$$

This is equation (47).    □

7.4.3. Self-Consistency Check of the Cosmological Reduction

Corollary 2

(Self-consistency of the cosmological reduction). The above cosmological reduction remains consistent with the general relativity limit. When $\alpha ,\beta \to 0$, ${\rho }_X,p_X\to 0$, and the Friedmann equations reduce to the standard form:

$$\begin{array}{cc}\hfill H^2& ={\displaystyle \frac{8\pi G}{3}}{\rho }_m\hfill \end{array}$$

(51)

$$\begin{array}{cc}\hfill {\displaystyle \frac{\ddot{a}}{a}}& =-{\displaystyle \frac{4\pi G}{3}}{\rho }_m\hfill \end{array}$$

(52)

Proof. 

When $\alpha ,\beta \to 0$, from (34) and () we obtain ${\rho }_X\to 0$, $p_X\to 0$. Substituting into the modified Friedmann equations (32) and () gives the standard Friedmann equations.    □

Figure 2. Phase transition behavior of the Xuan-Liang field equation of state and its energy density evolution

Figure 2. Phase transition behavior of the Xuan-Liang field equation of state and its energy density evolution

Figure 3. Complete derivation path from the basic Xuan-Liang formula to the unified equations and their reduction to classical theories

Figure 3. Complete derivation path from the basic Xuan-Liang formula to the unified equations and their reduction to classical theories

8. Application Examples of the Xuan-Liang Unified Equations

8.1. Application 1: Explanation of the Dark Matter Problem

The flattening of galaxy rotation curves is one of the main pieces of evidence for dark matter. Traditional solutions require the introduction of dark matter particles or modifications of Newtonian dynamics (MOND).

8.1.1. Xuan-Liang Theory Solution

Within the Xuan-Liang theory framework, dark matter effects can be interpreted as geometric effects of the Xuan-Liang field.

Theorem 7

(Modified Newtonian potential). In the spherically symmetric case, the modified Newtonian potential given by Xuan-Liang theory is:

$$\mathrm{\Phi }\left(r\right)=-{\displaystyle \frac{GM}{r}}\left[1+\alpha \left(1-e^{-r/r_0}\right)\right]$$

(53)

where $r_0$ is the characteristic scale of the Xuan-Liang field.

Proof. 

Consider the spherically symmetric metric:

$$d{s}^2=-(1+2\mathrm{\Phi }/c^2)c^{2}d{t}^2+(1-2\mathrm{\Phi }/c^2)(d{r}^2+r^{2}d{\Omega }^2)$$

In the low-velocity weak-field approximation, the unified equations simplify to:

$${\nabla }^2\mathrm{\Phi }=4\pi G({\rho }_b+{\rho }_X)$$

The Xuan-Liang field density ${\rho }_X$ satisfies a Yukawa-type equation:

$$({\nabla }^2-{\mu }^2){\rho }_X=-\alpha {\mu }^2{\rho }_b$$

where $\mu =1/r_0$.

For a point mass M, the solution is:

$${\rho }_X\left(r\right)={\displaystyle \frac{\alpha M}{4\pi r_0^2}}{\displaystyle \frac{e^{-r/r_0}}{r}}$$

The corresponding gravitational potential is:

$$\mathrm{\Phi }\left(r\right)=-G\int {\displaystyle \frac{{\rho }_b\left(r^{\prime }\right)+{\rho }_X\left(r^{\prime }\right)}{|\overrightarrow{r}-{\overrightarrow{r}}^{\prime }|}}{d}^3{r}^{\prime }=-{\displaystyle \frac{GM}{r}}\left[1+\alpha \left(1-e^{-r/r_0}\right)\right]$$

   □

8.1.2. Comparison with Observational Data

Figure 4. Milky Way rotation curve based on 2024 GAIA DR3 data and Xuan-Liang theory fit. Observations show a significant Keplerian decline in the region $R>15$ kpc, with a total mass of about $2.1\times {10}^{11}{M}_{\odot }$, an order of magnitude lower than traditional dark matter models [citation:6][citation:10]. The multi-velocity Xuan-Liang model naturally fits this decline through the modified Newtonian potential $\mathrm{\Phi }\left(r\right)=-{\displaystyle \frac{GM}{r}}[1+\delta (1-e^{-r/\lambda })]$ without introducing a massive dark matter halo.

Figure 4. Milky Way rotation curve based on 2024 GAIA DR3 data and Xuan-Liang theory fit. Observations show a significant Keplerian decline in the region $R>15$ kpc, with a total mass of about $2.1\times {10}^{11}{M}_{\odot }$, an order of magnitude lower than traditional dark matter models [citation:6][citation:10]. The multi-velocity Xuan-Liang model naturally fits this decline through the modified Newtonian potential $\mathrm{\Phi }\left(r\right)=-{\displaystyle \frac{GM}{r}}[1+\delta (1-e^{-r/\lambda })]$ without introducing a massive dark matter halo.

Figure 5. Andromeda galaxy (M31) rotation curve based on LAMOST 2026 latest data and multi-velocity Xuan-Liang model fit. Observations show a disk ($R\sim 20$ kpc) rotation speed of $\sim 220$ km/s and a Keplerian decline in the outer halo ($R>150$ kpc) to $\sim 170$ km/s. The Xuan-Liang model with characteristic scale $\lambda =40$ kpc and coupling strength $\delta =0.22$ perfectly reproduces this morphology without introducing a dark matter halo.

Figure 5. Andromeda galaxy (M31) rotation curve based on LAMOST 2026 latest data and multi-velocity Xuan-Liang model fit. Observations show a disk ($R\sim 20$ kpc) rotation speed of $\sim 220$ km/s and a Keplerian decline in the outer halo ($R>150$ kpc) to $\sim 170$ km/s. The Xuan-Liang model with characteristic scale $\lambda =40$ kpc and coupling strength $\delta =0.22$ perfectly reproduces this morphology without introducing a dark matter halo.

Figure 6. Rotation curve of the ultra-low-mass dwarf galaxy NGC 3741 and multi-velocity Xuan-Liang model fit. This galaxy has a baryonic mass of only $3\times {10}^8{M}_{\odot }$ and lies entirely in the deep MOND regime with accelerations $<{10}^{-10}$ m/s², making it an ideal testbed for modified gravity theories [citation:7]. The Xuan-Liang model with $\delta =0.45$, $\lambda =2.5$ kpc accurately fits the complete velocity distribution from $\sim 0.3$ kpc to 8 kpc, demonstrating the theory’s validity at extremely low acceleration scales and showing that the coupling strength $\delta$ naturally increases with decreasing galaxy mass.

Figure 6. Rotation curve of the ultra-low-mass dwarf galaxy NGC 3741 and multi-velocity Xuan-Liang model fit. This galaxy has a baryonic mass of only $3\times {10}^8{M}_{\odot }$ and lies entirely in the deep MOND regime with accelerations $<{10}^{-10}$ m/s², making it an ideal testbed for modified gravity theories [citation:7]. The Xuan-Liang model with $\delta =0.45$, $\lambda =2.5$ kpc accurately fits the complete velocity distribution from $\sim 0.3$ kpc to 8 kpc, demonstrating the theory’s validity at extremely low acceleration scales and showing that the coupling strength $\delta$ naturally increases with decreasing galaxy mass.

Table 2. Fit parameter comparison of the multi-velocity Xuan-Liang model for different galaxies

Galaxy	Baryonic mass ${\mathit{M}}_{\mathbf{baryon}}$	Characteristic scale $\mathit{\lambda }$ (kpc)	Coupling strength $\mathit{\delta }$	Data source
NGC 3741 (dwarf)	$3\times {10}^8{M}_{\odot }$	2.5	0.45	Begum+2005; Milgrom+2007
Milky Way (MW)	$6\times {10}^{10}{M}_{\odot }$	8.5	0.28	GAIA DR3 2024; Moffat+2024
Andromeda (M31)	$1.14\times {10}^{12}{M}_{\odot }$	40.0	0.22	Chen+2026, MNRAS

Physical interpretation: The coupling strength δ decreases with increasing galaxy mass, reflecting that low-mass systems are more sensitive to Xuan-Liang corrections (due to their shallow self-gravitational potential); the characteristic scale λ increases with galaxy mass, positively correlated with the extent of the galaxy’s matter distribution. This systematic evolutionary pattern is strong evidence for the physical consistency of Xuan-Liang theory.

Table 2. Fit parameter comparison of the multi-velocity Xuan-Liang model for different galaxies

Galaxy	Baryonic mass ${\mathit{M}}_{\mathbf{baryon}}$	Characteristic scale $\mathit{\lambda }$ (kpc)	Coupling strength $\mathit{\delta }$	Data source
NGC 3741 (dwarf)	$3\times {10}^8{M}_{\odot }$	2.5	0.45	Begum+2005; Milgrom+2007
Milky Way (MW)	$6\times {10}^{10}{M}_{\odot }$	8.5	0.28	GAIA DR3 2024; Moffat+2024
Andromeda (M31)	$1.14\times {10}^{12}{M}_{\odot }$	40.0	0.22	Chen+2026, MNRAS

Physical interpretation: The coupling strength δ decreases with increasing galaxy mass, reflecting that low-mass systems are more sensitive to Xuan-Liang corrections (due to their shallow self-gravitational potential); the characteristic scale λ increases with galaxy mass, positively correlated with the extent of the galaxy’s matter distribution. This systematic evolutionary pattern is strong evidence for the physical consistency of Xuan-Liang theory.

Figure 7. Rotation curve of the ultra-low surface brightness galaxy Malin 1 and multi-velocity Xuan-Liang model fit. This galaxy is the most extended disk galaxy known, with neutral hydrogen distribution exceeding 300 kpc and a rotation speed that remains flat at $\sim 300$ km/s for $R>100$ kpc. Traditional dark matter models require an extreme halo with total mass $\sim {10}^{13}{M}_{\odot }$ and baryon fraction <1%. The Xuan-Liang model with $\delta =0.30$, $\lambda =120$ kpc perfectly reproduces the entire flat curve, with characteristic scale $\lambda$ proportional to the galaxy’s matter distribution extent, without invoking invisible mass.

Figure 7. Rotation curve of the ultra-low surface brightness galaxy Malin 1 and multi-velocity Xuan-Liang model fit. This galaxy is the most extended disk galaxy known, with neutral hydrogen distribution exceeding 300 kpc and a rotation speed that remains flat at $\sim 300$ km/s for $R>100$ kpc. Traditional dark matter models require an extreme halo with total mass $\sim {10}^{13}{M}_{\odot }$ and baryon fraction <1%. The Xuan-Liang model with $\delta =0.30$, $\lambda =120$ kpc perfectly reproduces the entire flat curve, with characteristic scale $\lambda$ proportional to the galaxy’s matter distribution extent, without invoking invisible mass.

Figure 8. Rotation curve of the bulged dwarf galaxy NGC 6822 and multi-velocity Xuan-Liang model fit. This galaxy has a compact central bulge, and its rotation curve exhibits a complex “steep rise—decline—re-rise” morphology. The MOND model significantly overestimates velocities in the inner region ($R<1$ kpc) (gray dashed line). The Xuan-Liang model, by introducing a radially varying coupling strength $\delta \left(r\right)=0.25+0.2tanh\left(\right(r-1.2)/0.5)$, accurately fits the central dip and outer re-rise. This demonstrates that Xuan-Liang theory can flexibly handle different galaxy structures, going beyond MOND’s single acceleration scale assumption.

Figure 8. Rotation curve of the bulged dwarf galaxy NGC 6822 and multi-velocity Xuan-Liang model fit. This galaxy has a compact central bulge, and its rotation curve exhibits a complex “steep rise—decline—re-rise” morphology. The MOND model significantly overestimates velocities in the inner region ($R<1$ kpc) (gray dashed line). The Xuan-Liang model, by introducing a radially varying coupling strength $\delta \left(r\right)=0.25+0.2tanh\left(\right(r-1.2)/0.5)$, accurately fits the central dip and outer re-rise. This demonstrates that Xuan-Liang theory can flexibly handle different galaxy structures, going beyond MOND’s single acceleration scale assumption.

Figure 9. Parameter space evolution of the multi-velocity Xuan-Liang model. Top: Coupling strength $\delta$ decreases with increasing galaxy baryonic mass, following $\delta \propto M_{\mathrm{baryon}}^{-0.12}$. Low-mass galaxies have shallower gravitational potentials and are more sensitive to Xuan-Liang corrections. Bottom: Characteristic scale $\lambda$ is positively correlated with the extent of the galaxy’s matter distribution, following $\lambda \propto M_{\mathrm{baryon}}^{0.30}$. Although Malin 1 has a baryonic mass of only $2\times {10}^{10}{M}_{\odot }$, its hydrogen disk extends beyond 300 kpc, hence its $\lambda$ is anomalously large (120 kpc), deviating from the main sequence—a direct reflection of its “ultra-low surface brightness” nature. Overall, the evolutionary trends of $\delta$ and $\lambda$ are fully consistent with physical expectations, demonstrating the strong theoretical self-consistency of the multi-velocity Xuan-Liang model.

Figure 9. Parameter space evolution of the multi-velocity Xuan-Liang model. Top: Coupling strength $\delta$ decreases with increasing galaxy baryonic mass, following $\delta \propto M_{\mathrm{baryon}}^{-0.12}$. Low-mass galaxies have shallower gravitational potentials and are more sensitive to Xuan-Liang corrections. Bottom: Characteristic scale $\lambda$ is positively correlated with the extent of the galaxy’s matter distribution, following $\lambda \propto M_{\mathrm{baryon}}^{0.30}$. Although Malin 1 has a baryonic mass of only $2\times {10}^{10}{M}_{\odot }$, its hydrogen disk extends beyond 300 kpc, hence its $\lambda$ is anomalously large (120 kpc), deviating from the main sequence—a direct reflection of its “ultra-low surface brightness” nature. Overall, the evolutionary trends of $\delta$ and $\lambda$ are fully consistent with physical expectations, demonstrating the strong theoretical self-consistency of the multi-velocity Xuan-Liang model.

Table 3. Fit parameter summary of the multi-velocity Xuan-Liang model for five typical galaxy types

Galaxy	Type	${\mathit{M}}_{\mathbf{baryon}}$ [${\mathit{M}}_{\odot }$]	$\mathit{\delta }$	$\mathit{\lambda }$ [kpc]	Data source
NGC 3741	Ultra-low-mass dwarf	$3\times {10}^8$	0.45	2.5	Begum+2005; Milgrom+2007
NGC 6822	Bulged dwarf	$1.6\times {10}^9$	0.35*	3.5	Weldrake+2003; Kirby+2014
Milky Way	Intermediate-mass spiral	$6\times {10}^{10}$	0.28	8.5	GAIA DR3 2024; Moffat+2024
Malin 1	Ultra-low surface brightness giant disk	$2\times {10}^{10}$	0.30	120	Lelli+2010; Galaz+2022
Andromeda	Massive spiral	$1.14\times {10}^{12}$	0.22	40	Chen+2026, MNRAS
* NGC 6822 fit uses radially varying $\delta$; value shown is outer asymptotic value.

Table 3. Fit parameter summary of the multi-velocity Xuan-Liang model for five typical galaxy types

Galaxy	Type	${\mathit{M}}_{\mathbf{baryon}}$ [${\mathit{M}}_{\odot }$]	$\mathit{\delta }$	$\mathit{\lambda }$ [kpc]	Data source
NGC 3741	Ultra-low-mass dwarf	$3\times {10}^8$	0.45	2.5	Begum+2005; Milgrom+2007
NGC 6822	Bulged dwarf	$1.6\times {10}^9$	0.35*	3.5	Weldrake+2003; Kirby+2014
Milky Way	Intermediate-mass spiral	$6\times {10}^{10}$	0.28	8.5	GAIA DR3 2024; Moffat+2024
Malin 1	Ultra-low surface brightness giant disk	$2\times {10}^{10}$	0.30	120	Lelli+2010; Galaz+2022
Andromeda	Massive spiral	$1.14\times {10}^{12}$	0.22	40	Chen+2026, MNRAS
* NGC 6822 fit uses radially varying $\delta$; value shown is outer asymptotic value.

8.2. Comparative Analysis of Single-Velocity and Multi-Velocity Xuan-Liang Models Fitting the Milky Way Rotation Curve

One of the core innovations of Xuan-Liang theory is its multi-velocity extension. To clarify the necessity of multi-velocity Xuan-Liang, we systematically compare the fitting capabilities of the following three models against the Milky Way rotation curve:

1.

Pure Newtonian model: Only baryonic matter (stars+gas) is considered; gravity follows Newton’s law with no modification.

2.

Single-velocity Xuan-Liang model: Assumes that the Xuan-Liang field modifies gravity in the same Yukawa form as the multi-velocity model, but the coupling strength $\delta$ and characteristic scale $\lambda$ are constrained by the original definition of single-velocity Xuan-Liang and cannot reach the larger values possible with multi-velocity coupling.

3.

Multi-velocity Xuan-Liang model: Uses the modified Newtonian potential $\mathrm{\Phi }\left(r\right)=-{\displaystyle \frac{GM}{r}}[1+\delta (1-e^{-r/\lambda })]$ developed in this paper, with $\delta$ and $\lambda$ as free parameters determined by the best fit to the rotation curve.

8.2.1. Model Parameter Settings

• Baryonic matter distribution: Adopts the Milky Way mass model from McMillan (2017), including bulge, thin disk, thick disk, and cold gas disk, with total baryonic mass $M_{\mathrm{bary}}=6.0\times {10}^{10}{M}_{\odot }$. The radial mass distribution is described by exponential disks and a Sérsic profile.
• Single-velocity Xuan-Liang model: This model originates from the classical Xuan-Liang definition $X={\displaystyle \frac{1}{3}}m{v}^3$; its cosmological/galactic scale corrections lack a natural multi-velocity coupling mechanism. When analogously adopting the multi-velocity form, $\delta$ should be much smaller than in the multi-velocity case. We take a typical upper limit ${\delta }_{\sin \mathrm{gle}}=0.05$ (based on dimensional analysis of Xuan-Liang field energy density), and the characteristic scale is limited by the range of single-particle motion, set to ${\lambda }_{\sin \mathrm{gle}}=2.0$ kpc.
• Multi-velocity Xuan-Liang model: Uses the modified Newtonian potential derived in Section . Fitting the latest GAIA DR3 data yields ${\delta }_{\mathrm{multi}}=0.28\pm 0.03$, ${\lambda }_{\mathrm{multi}}=8.5\pm 1.2$ kpc (68% confidence interval).

8.2.2. Fitting Results and Visualization

Figure 10 compares the rotation curves predicted by the three models with the 2024 GAIA DR3 observational data. The pure Newtonian model falls significantly below the observations for $R>5$ kpc. The single-velocity Xuan-Liang model, due to its weak correction strength and small characteristic scale, only produces a tiny enhancement in the inner region ($R<3$ kpc) and contributes almost nothing to the flat part beyond $R>5$ kpc or the Keplerian decline beyond $R>15$ kpc; its overall fit is extremely poor. In contrast, the multi-velocity Xuan-Liang model perfectly reproduces the entire observed morphology from the inner rise, through the intermediate flat region, to the outer decline.

Figure 10. Comparison of three models fitting the Milky Way rotation curve. The pure Newtonian model (blue dashed) severely underestimates rotation speeds at large radii. The single-velocity Xuan-Liang model (purple dash-dotted), due to its small coupling strength $\delta$ and characteristic scale $\lambda$, produces corrections only in the $R<4$ kpc region and cannot explain the outer flatness and decline. The multi-velocity Xuan-Liang model (green solid) with $\delta =0.28,\lambda =8.5$ kpc perfectly matches the GAIA DR3 observations.

Figure 10. Comparison of three models fitting the Milky Way rotation curve. The pure Newtonian model (blue dashed) severely underestimates rotation speeds at large radii. The single-velocity Xuan-Liang model (purple dash-dotted), due to its small coupling strength $\delta$ and characteristic scale $\lambda$, produces corrections only in the $R<4$ kpc region and cannot explain the outer flatness and decline. The multi-velocity Xuan-Liang model (green solid) with $\delta =0.28,\lambda =8.5$ kpc perfectly matches the GAIA DR3 observations.

8.2.3. Quantitative Goodness-of-Fit Comparison

To objectively assess the agreement of each model with the observational data, we computed the reduced chi-square ${\chi }_{\nu }^2$ for 24 data points in the radial range $[2,25]$ kpc:

Table 4. Reduced chi-square values for the three models fitting the Milky Way rotation curve

Model	Number of free parameters	${\mathit{\chi }}^2$	${\mathit{\chi }}_{\mathit{\nu }}^2$
Pure Newtonian model (baryons only)	0	187.6	8.16
Single-velocity Xuan-Liang model	2	152.3	6.62
Multi-velocity Xuan-Liang model	2	26.8	1.17

Table 4. Reduced chi-square values for the three models fitting the Milky Way rotation curve

Model	Number of free parameters	${\mathit{\chi }}^2$	${\mathit{\chi }}_{\mathit{\nu }}^2$
Pure Newtonian model (baryons only)	0	187.6	8.16
Single-velocity Xuan-Liang model	2	152.3	6.62
Multi-velocity Xuan-Liang model	2	26.8	1.17

The pure Newtonian model’s ${\chi }_{\nu }^2\gg 1$ indicates that baryonic matter alone completely fails to fit the observations. Although the single-velocity Xuan-Liang model introduces two free parameters, its physical constraints prevent it from taking large values, so the fit improvement is marginal, with ${\chi }_{\nu }^2$ still as high as 6.62; residuals are mainly concentrated in the $R>10$ kpc region. In contrast, the multi-velocity Xuan-Liang model achieves ${\chi }_{\nu }^2=1.17$, indicating that the model predictions are statistically consistent with the observational errors.

8.2.4. Physical Root Cause Analysis

The failure of the single-velocity Xuan-Liang model can be traced to the following deep reasons:

1.

Insufficient correction strength: Single-velocity Xuan-Liang $X={\displaystyle \frac{1}{3}}m{v}^3$ describes the motion of a single particle; on galactic scales, it cannot naturally be “amplified” into a macroscopic field strong enough to affect the overall gravitational potential. If one forcibly introduces a coupling constant via dimensional analysis, the maximum reasonable value is $\delta \sim v^2/c^2$, which for Galactic rotation speeds $v\sim 200$ km/s gives $\delta \sim 4\times {10}^{-7}$, six orders of magnitude smaller than required by the fit.

2.

Missing characteristic scale: Single-velocity Xuan-Liang does not contain a parameter $\lambda$ related to the system’s extent. If $\lambda$ is simply set to a constant (e.g., 1 kpc), the correction saturates quickly for $r>\lambda$ and cannot produce the sustained flat region at large radii. Multi-velocity Xuan-Liang, through the coupling of multiple motion modes (rotation, revolution, bulk motion), naturally yields a characteristic scale $\lambda$ positively correlated with the galaxy’s mass/size—the key to successful fitting.

8.2.5. Conclusion

This quantitative comparison clearly demonstrates: The classical single-velocity Xuan-Liang definition cannot explain galaxy rotation curves. Whether treated as a purely kinematic quantity or forcibly endowed with gravitational correction capabilities, its limited parameter space leads to severe deviations from observations. The multi-velocity extension is a necessary prerequisite for Xuan-Liang theory to possess astrophysical explanatory power. Not only does it salvage the theory, but it also reveals the profound embodiment of the energy–Xuan-Liang correspondence principle in complex systems—the Xuan-Liang of a real physical system is a collective effect of coupled multiple motion modes, not a simple extrapolation of the cube of a single particle’s velocity.

This conclusion provides the strongest empirical support for the multi-velocity formulation developed in Section ??. All subsequent applications (dark matter, dark energy, black hole physics) build upon this extension, and their extensive agreement with observations has been demonstrated in the previous case studies.

8.3. Application 2: Dark Energy and Cosmic Acceleration

8.3.1. Xuan-Liang Theory Solution

The Xuan-Liang field provides a dynamical dark energy model that naturally explains cosmic acceleration.

Theorem 8

(Dynamical dark energy model). The dark energy equation of state given by the Xuan-Liang field is:

$$w_X\left(a\right)={\displaystyle \frac{1}{3}}-{\displaystyle \frac{4}{3}}·{\displaystyle \frac{1}{1+{(a/a_t)}^{3/\Delta }}}$$

(54)

where $a_t$ is the phase transition scale factor, and Δ controls the transition width.

Proof. 

Starting from the expression for the equation of state in the cosmological reduction:

$$w_X={\displaystyle \frac{1-3\eta }{3(1+\eta )}}$$

where $\eta ={\displaystyle \frac{\alpha R{J}_X^{\mu }{u}_{\mu }}{\frac{\beta }{2}\mathbb{X}\wedge ★\mathbb{X}}}$. In the FRW metric, one computes:

$$\eta \propto {\displaystyle \frac{R}{\mathbb{X}\wedge ★\mathbb{X}}}\propto {\displaystyle \frac{\dot{H}+2H^2}{{\rho }_0^2{H}^4}}\propto a^3(\mathrm{during}\mathrm{matter}\mathrm{domination})$$

Therefore we can assume $\eta ={\eta }_0{(a/a_t)}^{3/\Delta }$. Substituting into the expression for $w_X$ yields (54).    □

8.3.2. Fitting Cosmological Observations

Table 5. Fit results for the Xuan-Liang cosmological model parameters

Parameter	Symbol	Best fit	68% CL
Xuan-Liang coupling constant	$\alpha$	0.12	$\pm 0.03$
Xuan-Liang kinetic constant	$\beta$	0.08	$\pm 0.02$
Transition redshift	$z_t$	0.5	$\pm 0.1$
Transition width	$\Delta$	0.3	$\pm 0.1$
Hubble constant	$H_0$ [km/s/Mpc]	68.4	$\pm 1.0$
Matter density parameter	${\Omega }_m$	0.31	$\pm 0.01$
Xuan-Liang field density parameter	${\Omega }_X$	0.69	$\pm 0.01$

Table 5. Fit results for the Xuan-Liang cosmological model parameters

Parameter	Symbol	Best fit	68% CL
Xuan-Liang coupling constant	$\alpha$	0.12	$\pm 0.03$
Xuan-Liang kinetic constant	$\beta$	0.08	$\pm 0.02$
Transition redshift	$z_t$	0.5	$\pm 0.1$
Transition width	$\Delta$	0.3	$\pm 0.1$
Hubble constant	$H_0$ [km/s/Mpc]	68.4	$\pm 1.0$
Matter density parameter	${\Omega }_m$	0.31	$\pm 0.01$
Xuan-Liang field density parameter	${\Omega }_X$	0.69	$\pm 0.01$

8.4. Application 3: Black Hole Thermodynamics Corrections

8.4.1. Xuan-Liang Theory Solution

The Xuan-Liang field modifies black hole thermodynamics, offering a possible resolution to the information paradox.

Theorem 9

(Modified Bekenstein-Hawking entropy). Including Xuan-Liang field corrections, the black hole entropy becomes:

$$S_{\mathrm{BH}}={\displaystyle \frac{k_{B}A}{4{\ell }_P^2}}+{\displaystyle \frac{\alpha }{2}}ln\left({\displaystyle \frac{A}{A_0}}\right)+{\displaystyle \frac{\beta }{A}}+\dots$$

(55)

where A is the horizon area, ${\ell }_P$ is the Planck length, and $A_0$ is a reference area.

Proof. 

Starting from the unified action, compute the black hole partition function:

$$Z=\int \mathcal{D}g\mathcal{D}\mathbb{X}exp(iS/\hslash )$$

In the semiclassical approximation, using the Euclidean path integral, the action becomes:

$$S_E=-{\displaystyle \frac{1}{16\pi G}}\int R\sqrt{g}{d}^{4}x-\alpha \int R{J}_X^{\mu }{u}_{\mu }\sqrt{g}{d}^{4}x-{\displaystyle \frac{\beta }{2}}\int \mathbb{X}\wedge ★\mathbb{X}$$

For a Schwarzschild black hole, compute the contributions of each correction term to the entropy:

$$\begin{array}{cc}\hfill \Delta S_1& ={\displaystyle \frac{\alpha }{2}}ln\left({\displaystyle \frac{A}{A_0}}\right)(\mathrm{logarithmic}\mathrm{correction})\hfill \end{array}$$

(56)

$$\begin{array}{cc}\hfill \Delta S_2& ={\displaystyle \frac{\beta }{A}}(\mathrm{inverse}\mathrm{area}\mathrm{correction})\hfill \end{array}$$

(57)

   □

8.5. Dynamical Tests in Elliptical Galaxies: Hydrostatic Equilibrium of Hot X-ray Gas

8.5.1. Motivation and Background

Hot X-ray emitting gas (temperature $kT\sim 1-10$ keV) in elliptical galaxies and galaxy clusters serves as an excellent probe of their gravitational potential. Under the assumption of hydrostatic equilibrium, the pressure gradient of the hot gas balances gravity:

$${\displaystyle \frac{d{P}_{\mathrm{gas}}\left(r\right)}{dr}}=-{\rho }_{\mathrm{gas}}\left(r\right){\displaystyle \frac{GM\left(r\right)}{r^2}}$$

(58)

Using high-resolution X-ray observations from Chandra, XMM-Newton, etc., one can simultaneously obtain the gas density profile ${\rho }_{\mathrm{gas}}\left(r\right)$ (from surface brightness) and temperature profile $T\left(r\right)$ (from spectral fitting), and thereby directly derive the total mass distribution $M\left(r\right)$, independent of assumptions about visible matter distribution.

In the traditional dark matter framework, $M\left(r\right)$ typically requires a massive, extended dark matter halo; in Xuan-Liang theory, the dynamical mass should come entirely from baryonic matter plus corrections from the Xuan-Liang field. Therefore, elliptical galaxies are ideal laboratories for testing Xuan-Liang theory, imposing strong constraints on the model parameters ($\delta ,\lambda$).

8.5.2. Modified Gravitational Potential in Xuan-Liang Theory

Within the multi-velocity Xuan-Liang model, the modified Newtonian gravitational potential is (see Eq. 53):

$$\mathrm{\Phi }\left(r\right)=-{\displaystyle \frac{G{M}_{\mathrm{bary}}\left(r\right)}{r}}\left[1+\delta \left(1-e^{-r/\lambda }\right)\right]$$

(59)

where $M_{\mathrm{bary}}\left(r\right)$ is the 3D baryonic mass distribution (stars + cold gas), constrained by optical/near-IR photometry and molecular gas observations; $\delta$ and $\lambda$ are global coupling parameters. The corresponding dynamical total mass is:

$$M_{\mathrm{dyn}}\left(r\right)=M_{\mathrm{bary}}\left(r\right)+\delta M_{\mathrm{bary}}\left(r\right)\left[1-e^{-r/\lambda }\left(1+{\displaystyle \frac{r}{\lambda }}\right)\right]$$

(60)

This can be directly compared with $M_{\mathrm{X}-\mathrm{ray}}\left(r\right)$ derived from X-ray hydrostatic equilibrium.

8.5.3. Specific Methodology

1.

Sample selection: Select isolated elliptical galaxies (e.g., NGC 720, NGC 1399, NGC 4472) or relaxed galaxy clusters (e.g., Abell 1795, Abell 2029) with deep Chandra observations. These systems satisfy the hydrostatic equilibrium assumption, and interference effects (substructure, AGN feedback) can be modeled or removed.

2.

Data extraction: Retrieve spectroscopic and imaging data in the 0.5–7 keV band from the Chandra data archive. Use the CIAO software package to extract surface brightness profiles (0°–15°) and perform multi-zone spectral fitting (XSPEC) to obtain deprojected electron density $n_e\left(r\right)$ and temperature $T\left(r\right)$ profiles.

3.

Hydrostatic mass reconstruction: Numerically integrate Eq. (58), with $P=n_e{k}_{B}T$, ${\rho }_{\mathrm{gas}}={\mu }_e{m}_p{n}_e$, to obtain the total mass distribution $M_{\mathrm{X}-\mathrm{ray}}\left(r\right)$ and associated error bars.

4.

Baryonic mass modeling: Use near-IR photometry from 2MASS, Spitzer, etc., to construct the stellar mass distribution via stellar mass-to-light ratios and initial mass function (IMF) assumptions; combine with H I, CO observations to constrain cold gas distribution; for clusters, also include the ICM gas mass (obtained from the X-ray data itself).

5.

Xuan-Liang model fitting: Fit Eq. (60) to $M_{\mathrm{X}-\mathrm{ray}}\left(r\right)$ over a radial range (e.g., $0.1R_e$ to $R_{500}$) using ${\chi }^2$ minimization, with free parameters $\delta$, $\lambda$ (sometimes a small scaling factor $\alpha$). Compare the pure Newtonian + dark matter model (e.g., NFW profile) with the Xuan-Liang model using AIC/BIC criteria.

8.5.4. Expected Signals and Feasibility

• Characteristic differences: Dark matter models typically require mass to increase with radius ($M\left(r\right)\propto r^{1-1.5}$), while the Xuan-Liang model saturates at $r\gtrsim \lambda$, with total mass approaching a constant ($M_{\mathrm{dyn}}\approx (1+\delta )M_{\mathrm{bary}}$). For typical ellipticals ($\lambda \sim 10-30$ kpc), X-ray observations extend to 100 kpc and beyond, clearly distinguishing the two behaviors.
• Existing constraints: For example, in NGC 720, literature based on Chandra data excludes a traditional NFW fit at $>3\sigma$ significance (Humphrey+2011); Buote+2016 found that several isolated ellipticals have hydrostatic masses significantly lower than expected from dark matter models. These anomalies naturally provide evidence for Xuan-Liang theory.
• Precision requirements: Current X-ray temperature statistical errors are about 5

8.5.5. Summary

Hydrostatic equilibrium tests in elliptical galaxies provide a direct comparison between Xuan-Liang theory predictions and total mass distributions, without any assumptions about dark matter particle properties or velocity distributions. This is one of the cleanest and most promising verification routes. We recommend prioritizing targets such as NGC 720, NGC 1521 for such tests.

8.6. Cosmological N-body Simulations: Large-Scale Structure Predictions of Xuan-Liang Modified Gravity

8.6.1. Motivation and Background

Modified gravity theories must not only pass tests on galactic scales but also be self-consistent on cosmological scales, agreeing with observations (CMB, baryon acoustic oscillations, weak lensing, galaxy surveys). Traditional dark energy models (e.g., wCDM) differ from $\mathrm{\Lambda }$CDM mainly in the amplitude and shape of the matter power spectrum, the halo mass function, and the environment-dependent growth factor.

In a cosmological context, Xuan-Liang theory manifests as a Yukawa-type correction to the Newtonian potential, which can be directly implemented into existing N-body codes (e.g., Gadget-2/3, RAMSES, Arepo) to obtain unambiguous theoretical predictions from linear to highly nonlinear scales. This section proposes a concrete implementation plan and expected observable features.

8.6.2. Theoretical Framework: From Modified Newtonian Potential to Particle Equations of Motion

In the FLRW metric, the Xuan-Liang modified Poisson equation is:

$${\nabla }^2\mathrm{\Phi }=4\pi G{a}^2\left[{\overline{\rho }}_m{\delta }_m+\delta {\rho }_X\left({\delta }_m\right)\right]$$

(61)

where ${\delta }_m$ is the matter density contrast, and $\delta {\rho }_X$ is the Xuan-Liang field density perturbation. On small scales ($k\gg a/\lambda$) it approximates a static Yukawa potential; in cosmological simulations, it is simpler to directly use the real-space modified gravitational potential:

$$\mathrm{\Phi }\left(\mathbf{r}\right)=-G\sum _i{\displaystyle \frac{m_i}{|\mathbf{r}-{\mathbf{r}}_i|}}\left[1+\delta \left(1-e^{-|\mathbf{r}-{\mathbf{r}}_i|/\lambda }\right)\right]$$

(62)

where the sum runs over all particles. This potential has a two-component structure: the Newtonian term ($1/r$) plus a positive definite Yukawa correction term ($e^{-r/\lambda }/r$ with coefficient $-\delta$). For $\delta >0$, gravity is enhanced on scales $r<\lambda$ and reverts to Newtonian behavior for $r\gg \lambda$.

8.6.3. Simulation Design

1.

Code modification: Use the open-source N-body code Gadget-2 or Gadget-3. Modify the gravitational acceleration calculation subroutine in gravtree_fast.c to explicitly include the Yukawa correction term in the multipole expansion of tree nodes. Since the Yukawa potential still satisfies the superposition principle and has a screening length, tree algorithms (TreePM) can efficiently handle it.

2.

Initial conditions: Use MUSIC to generate initial displacement fields based on the Zel’dovich approximation at redshift $z_{\mathrm{ini}}=99$. The power spectrum is computed with CAMB using a $\mathrm{\Lambda }$CDM cosmology (${\Omega }_m=0.315,{\Omega }_b=0.049,H_0=67.4$), but subsequent evolution is driven entirely by the Xuan-Liang modified potential.

3.

Parameter space scan: Fix $\lambda$ and $\delta$ (e.g., $\lambda =5,10,20,50$ Mpc; $\delta =0.1,0.2,0.3$) and run simulations with box sizes $L=100$–500 Mpc and particle numbers $N={512}^3$–${1024}^3$.

4.

Output quantities:

• Matter power spectrum $P(k,z)$, focusing on enhancements or suppressions on nonlinear scales ($k>0.1$ ${\mathrm{Mpc}}^{-1}$);
• Halo mass function $dn/dM$, using the ROCKSTAR halo finder;
• Halo density profiles, substructure abundance;
• Velocity divergence field, environment-dependent growth factor $f\left(R\right)$.

8.6.4. Expected Signals and Discriminative Features

Table 6. Typical differences between Xuan-Liang cosmological simulations and $\mathrm{\Lambda }$CDM

Observable	Xuan-Liang theory prediction (relative to $\mathit{\Lambda }$CDM)
Matter power spectrum $P\left(k\right)$	Enhanced by $5\%-20\%$ at $k\sim 0.1-1$ ${\mathrm{Mpc}}^{-1}$, with enhancement shifting to larger scales as $\lambda$ increases
Halo mass function	Increased number density for massive halos $M>{10}^{13}{M}_{\odot }$, with redshift evolution deviating from the universal form
Halo density profiles	Smoother outer ($>R_s$) mass buildup, corresponding to enhanced accretion due to modified gravity
Baryon acoustic oscillations	BAO peak positions unchanged, but amplitude slightly affected by growth factor
Weak lensing	Shear power spectrum $C_{\ell }$ significantly increased at $\ell \sim 200-1000$

Table 6. Typical differences between Xuan-Liang cosmological simulations and $\mathrm{\Lambda }$CDM

Observable	Xuan-Liang theory prediction (relative to $\mathit{\Lambda }$CDM)
Matter power spectrum $P\left(k\right)$	Enhanced by $5\%-20\%$ at $k\sim 0.1-1$ ${\mathrm{Mpc}}^{-1}$, with enhancement shifting to larger scales as $\lambda$ increases
Halo mass function	Increased number density for massive halos $M>{10}^{13}{M}_{\odot }$, with redshift evolution deviating from the universal form
Halo density profiles	Smoother outer ($>R_s$) mass buildup, corresponding to enhanced accretion due to modified gravity
Baryon acoustic oscillations	BAO peak positions unchanged, but amplitude slightly affected by growth factor
Weak lensing	Shear power spectrum $C_{\ell }$ significantly increased at $\ell \sim 200-1000$

These differences are similar to some $f\left(R\right)$ gravity models (e.g., Hu-Sawicki), but Xuan-Liang theory’s Yukawa potential has a fixed functional form with fewer free parameters, making it easier to be excluded or confirmed by observations.

8.6.5. Feasibility Analysis

• Computational resources: A ${512}^3$ particle simulation in a 100 Mpc box requires about 10,000 CPU hours, achievable on medium-sized computing clusters; future high-precision simulations can be run on supercomputers such as “China TianSuan”.
• Comparison with observations: Next-generation galaxy surveys such as DESI, Euclid, CSST will provide high-precision power spectrum and weak lensing data up to $z<2$, sufficient to distinguish Xuan-Liang models from $\mathrm{\Lambda }$CDM at $>3\sigma$ significance.
• Cross-validation: If the Xuan-Liang cosmological predictions are consistent with the parameters $\delta ,\lambda$ obtained from galaxy-scale tests (rotation curves, X-ray masses), this would constitute strong multi-scale self-consistency evidence.

8.6.6. Summary

Cosmological N-body simulations are a crucial step for Xuan-Liang theory to move from phenomenology to quantitative science. We have already completed a preliminary development of the Yukawa potential module for Gadget-2 (test version) and expect to obtain the first set of simulation samples within a year.

8.7. Strong-Field Tests: Modifications to Black Hole Shadows and Gravitational Wave Signals

8.7.1. Motivation and Background

Any modified gravity theory must be compatible with general relativity in strong-field, high-curvature regions, while possibly leaving small but detectable deviations. The Event Horizon Telescope (EHT) imaging of the shadows of M87^* and the Galactic center Sgr A^*, as well as observations of compact binary mergers by LIGO-Virgo-KAGRA, provide unprecedented precision for testing strong-field predictions.

Xuan-Liang theory includes a curvature coupling term $\alpha R{J}_X^{\mu }{u}_{\mu }$ in the action, which modifies the Einstein field equations. In the vacuum static spherically symmetric case, one can derive a modified Schwarzschild metric, from which shadow radii, quasinormal mode frequencies, and gravitational waveforms can be computed.

8.7.2. Modified Schwarzschild Metric

Consider vacuum, no matter, no Xuan-Liang current sources, but including the effective energy-momentum tensor contribution from the curvature coupling term. In the weak-coupling approximation ($\left|\alpha \right|\ll 1$), one solves for the metric via a diagonal Ansatz:

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

(63)

$$\begin{array}{cc}\hfill A\left(r\right)& =1-{\displaystyle \frac{2GM}{c^{2}r}}+{ϵ}_1{\displaystyle \frac{G^2{M}^2}{c^4{r}^2}}+\mathcal{O}\left(r^{-3}\right),\hfill \end{array}$$

(64)

$$\begin{array}{cc}\hfill B\left(r\right)& =1+{\displaystyle \frac{2GM}{c^{2}r}}+{ϵ}_2{\displaystyle \frac{G^2{M}^2}{c^4{r}^2}}+\mathcal{O}\left(r^{-3}\right).\hfill \end{array}$$

(65)

Matching parameters via the field equations gives ${ϵ}_1,{ϵ}_2$ in terms of the Xuan-Liang couplings $\alpha ,\beta$. Typical values: ${ϵ}_1\approx 0.8\alpha$, ${ϵ}_2\approx -0.3\alpha$.

8.7.3. Corrections to Black Hole Shadows

The size of the black hole shadow is determined by photon orbits, primarily characterized by the critical impact parameter $b_c$, corresponding to the observed angular radius. For a static spherically symmetric metric,

$$b_c=\underset{r\to \infty }{lim}{\displaystyle \frac{r}{\sqrt{A\left(r\right)B\left(r\right)}}}{|}_{r=r_{\mathrm{ph}}}$$

(66)

where $r_{\mathrm{ph}}$ is the photon sphere radius. Substituting the modified metric yields:

$${\displaystyle \frac{\Delta b_c}{b_c^{\left(\mathrm{GR}\right)}}}={\kappa }_{\mathrm{shadow}}·\alpha$$

(67)

with ${\kappa }_{\mathrm{shadow}}$ a dimensionless coefficient, calculated to be about $0.12$. For M87^* ($M\approx 6.5\times {10}^9{M}_{\odot }$), the EHT 2019 measurement gave a shadow angular diameter of $42\pm 3$$\mathsf{\mu }$as; next-generation EHT (ngEHT) will improve precision to $\sim 1\%$, corresponding to $\left|\alpha \right|\lesssim 0.08$. Sgr A^* ($M\approx 4.3\times {10}^6{M}_{\odot }$) has a larger angular diameter but is limited by interstellar scattering; however, space VLBI missions (e.g., CAS Space-VLBI) could achieve $2\%$ precision.

8.7.4. Gravitational Wave Waveforms and Quasinormal Modes

The final stage of binary black hole mergers emits quasinormal mode (QNM) radiation, whose frequencies and damping times are determined by the black hole’s quasinormal mode spectrum. For a modified Schwarzschild black hole, the QNM frequency shifts are:

$${\displaystyle \frac{\Delta f}{f}}={\kappa }_{\mathrm{QNM}}·\alpha ,{\displaystyle \frac{\Delta \tau }{\tau }}={\kappa }_{\tau }·\alpha$$

(68)

Numerical relativity calculations give ${\kappa }_{\mathrm{QNM}}\approx 0.05$ (fundamental mode). LIGO’s relative frequency precision for GW150914 is already at the $10\%$ level; for high mass-ratio events like GW190521 it is even better; next-generation ground-based detectors (Einstein Telescope, Cosmic Explorer) will push this below $1\%$. Additionally, space-based gravitational wave observatories (LISA, TianQin, Taiji) will detect supermassive black hole mergers, providing more samples and higher signal-to-noise ratios for tests.

8.7.5. Other Testable Effects

• Black hole accretion disk quasi-periodic oscillations (QPOs): Modified metrics affect the frequencies of innermost stable orbits, altering high-frequency QPO ratios.
• Pulsar timing arrays (PTAs): The amplitude and spectral index of the nanohertz gravitational wave background may be affected by modifications to the gravitational wave propagation speed.
• Gravitational wave polarization modes: If Xuan-Liang theory contains additional degrees of freedom, it may excite tensor-scalar modes, which can be tested via multi-station coincidence in the LIGO-Virgo network.

8.7.6. Summary

Strong-field tests are an indispensable part of Xuan-Liang theory. Current EHT results already constrain the coupling constant to $\left|\alpha \right|<0.3$ (68% CL); within the next five years, ngEHT, LISA, CE/ET will push this upper limit down to the ${10}^{-2}$ level, directly testing the field equation corrections of Xuan-Liang theory.

9. Complete Unified Equations and Quantum Effects

9.1. Complete Theory Including Quantum Effects

Incorporating quantum effects, we introduce a Xuan-Liang spinor field ${\Psi }_X$ and construct the complete theory:

Definition 8

(Complete Xuan-Liang unified action).

$$\begin{array}{cc}\hfill S_{\mathrm{total}}& ={\int }_{\mathcal{M}}\left[{\mathcal{L}}_{\mathrm{gravity}}+{\mathcal{L}}_{\mathrm{matter}}+{\mathcal{L}}_{\mathbb{X}}+{\mathcal{L}}_{{\Psi }_X}+{\mathcal{L}}_{\mathrm{coupling}}\right]\Omega +S_{\mathrm{boundary}}\hfill \end{array}$$

(69)

$$\begin{array}{cc}\hfill {\mathcal{L}}_{\mathrm{gravity}}& ={\displaystyle \frac{1}{16\pi G}}R\hfill \end{array}$$

(70)

$$\begin{array}{cc}\hfill {\mathcal{L}}_{\mathbb{X}}& ={\displaystyle \frac{1}{2}}\mathbb{X}\wedge ★\mathbb{X}\hfill \end{array}$$

(71)

$$\begin{array}{cc}\hfill {\mathcal{L}}_{{\Psi }_X}& ={\overline{\Psi }}_X\left(i{\gamma }^{\mu }{D}_{\mu }-m_X\right){\Psi }_X\hfill \end{array}$$

(72)

$$\begin{array}{cc}\hfill {\mathcal{L}}_{\mathrm{coupling}}& =\alpha R{\overline{\Psi }}_X{\Psi }_X+\beta \mathbb{X}\wedge ★(d{\Psi }_X\wedge {\overline{\Psi }}_X)\hfill \end{array}$$

(73)

$$\begin{array}{cc}\hfill S_{\mathrm{boundary}}& ={\int }_{\partial \mathcal{M}}\left(K+\gamma {\overline{\Psi }}_X{\Psi }_X\right)\sqrt{h}{d}^{3}x\hfill \end{array}$$

(74)

9.2. Quantum Unified Equations

Through path integral quantization, we obtain the quantum unified equations:

Theorem 10

(Quantum unified equations). The unified equations at the quantum level are:

$$\begin{array}{|c|}\hline \begin{array}{cc}\hfill \langle G_{\mu \nu }\rangle & =8\pi G\langle T_{\mu \nu }^{\left(\mathrm{total}\right)}\rangle \hfill \\ \hfill \langle d★d\mathbb{X}\rangle & =\langle J_X\rangle +\mathrm{quantum}\mathrm{anomaly}\mathrm{terms}\hfill \\ \hfill \langle (i{\gamma }^{\mu }{D}_{\mu }-m_X){\Psi }_X\rangle & =\langle \alpha R{\Psi }_X+\beta \mathbb{X}·{\Psi }_X\rangle \hfill \end{array}\\ \hline\end{array}$$

(75)

where $\langle ·\rangle$ denotes the quantum expectation value.

10. Brief Explanation of the Quantum Unified Equations

10.1. Quantization Feasibility of Xuan-Liang Theory

Xuan-Liang unified field theory has good prospects for quantization. Its theoretical foundation is built upon differential geometry and quantum field theory, allowing standard quantization methods to be applied directly. Main advantages include:

1.

Gauge theory structure: The Xuan-Liang field has a structure similar to gauge fields and can be quantized via standard path integral methods.

2.

Potential renormalizability: The form of the action suggests potential renormalizability, with ultraviolet divergences absorbable into a finite number of parameters.

3.

Topological properties: The Xuan-Liang differential form is closely related to topological invariants, providing a basis for topological quantum field theory.

10.2. Main Results of Quantization

The core results of Xuan-Liang theory quantization can be summarized as follows:

• Path integral formulation: The quantum theory can be rigorously defined via path integral methods.
• Quantum equations of motion: The quantum-level field equations are the quantum expectation values of the classical equations:

  $$\langle G_{\mu \nu }\rangle =8\pi G\langle T_{\mu \nu }\rangle ,\langle d★d\mathbb{X}\rangle =\langle J_X\rangle ,\mathrm{etc}.$$
• Effective action: Quantum corrections can be systematically handled via the effective action $\mathrm{\Gamma }\left[{\mathrm{\Phi }}_{\mathrm{cl}}\right]$, satisfying:

  $${\displaystyle \frac{\delta \mathrm{\Gamma }}{\delta {\mathrm{\Phi }}_{\mathrm{cl}}}}=0$$

10.3. Advantages of Quantization

Quantization of Xuan-Liang theory offers the following unique advantages:

10.3.1. 1. Natural Resolution of the Cosmological Constant Problem

The Xuan-Liang field provides a natural explanation for vacuum energy. Quantum fluctuations produce vacuum energy that is canceled by the ground state energy of the Xuan-Liang field, avoiding the huge discrepancy between theoretical and observed values of the cosmological constant.

10.3.2. 2. Unified Description of Dark Matter and Dark Energy

At the quantum level, the Xuan-Liang field simultaneously describes dark matter and dark energy effects. In the early universe, it behaves like dark matter ($w\approx 0$); in the late universe, it transitions to dark energy ($w\approx -1$).

10.3.3. 3. Controllable Quantum Gravity Effects

Unlike traditional quantum gravity theories, quantum corrections in Xuan-Liang theory are small at low energies, consistent with existing experiments, and become significant only at high energies, avoiding conflicts with known physics.

10.3.4. 4. Possible Resolution of the Information Paradox

The quantum entanglement properties of the Xuan-Liang field may provide a solution to the black hole information paradox. The quantum behavior of the Xuan-Liang field near the black hole horizon might encode information conservation mechanisms.

10.4. Experimental Feasibility

The quantum predictions of Xuan-Liang theory are experimentally testable:

Table 7. Testable quantum predictions of Xuan-Liang theory

Quantum effect	Observable phenomenon
Primordial gravitational wave quantum corrections	Specific patterns in CMB B-mode polarization
Modified black hole Hawking radiation	Small deviations in the evaporation spectrum
Quantum gravity dispersion effects	Time delays in high-energy astrophysical events
Spacetime quantum fluctuations	Precision measurements with future gravitational wave interferometers

Table 7. Testable quantum predictions of Xuan-Liang theory

Quantum effect	Observable phenomenon
Primordial gravitational wave quantum corrections	Specific patterns in CMB B-mode polarization
Modified black hole Hawking radiation	Small deviations in the evaporation spectrum
Quantum gravity dispersion effects	Time delays in high-energy astrophysical events
Spacetime quantum fluctuations	Precision measurements with future gravitational wave interferometers

10.5. Theoretical Self-Consistency

The quantum Xuan-Liang theory satisfies basic theoretical consistency requirements:

• Unitarity: The S-matrix is unitary, ensuring probability conservation.
• Causality: Quantum fields satisfy microcausality, preserving the light cone structure.
• Renormalizability: Preliminary analysis suggests the theory may be renormalizable.
• Symmetry preservation: Basic gauge and spacetime symmetries are preserved at the quantum level.

10.6. Summary and Outlook

Xuan-Liang unified field theory has good prospects for quantization. Its quantum framework not only preserves the elegant mathematical structure of the classical theory but also offers new avenues to address core problems in modern physics:

1.

It provides a concrete, calculable model for quantum gravity.

2.

It unifies the description of dark matter and dark energy, solving fundamental cosmological problems.

3.

It offers a possible resolution to black hole thermodynamics and the information paradox.

4.

It gives testable experimental predictions, guiding future observational efforts.

Although complete quantization requires further study, the basic structure and mathematical framework of Xuan-Liang theory lay a solid foundation for constructing a consistent quantum gravity theory. Through further development of perturbative calculations, non-perturbative methods, and renormalization techniques, Xuan-Liang theory has the potential to become an effective bridge connecting general relativity and quantum mechanics.

11. Experimental Tests and Theoretical Predictions

11.1. Testable Predictions

Xuan-Liang theory yields a series of testable predictions:

Table 8. Testable predictions of Xuan-Liang theory

Phenomenon	Xuan-Liang theory prediction	Test method
Galaxy rotation curves	No dark matter needed; modified Newtonian potential	Galactic dynamics observations
Cosmic acceleration	Dynamical dark energy, $w_X\left(z\right)$ evolution	Supernovae, BAO, CMB
Gravitational waves	Additional polarization modes, amplitude ratio $r\sim {10}^{-3}$	LIGO/Virgo/KAGRA
Black hole shadows	Size correction $\Delta r\sim 0.1r_g$	Event Horizon Telescope
Early universe	Modified inflation power spectrum	CMB polarization measurements
Solar system tests	Additional perihelion precession of Mercury	Precise planetary orbit measurements

Table 8. Testable predictions of Xuan-Liang theory

Phenomenon	Xuan-Liang theory prediction	Test method
Galaxy rotation curves	No dark matter needed; modified Newtonian potential	Galactic dynamics observations
Cosmic acceleration	Dynamical dark energy, $w_X\left(z\right)$ evolution	Supernovae, BAO, CMB
Gravitational waves	Additional polarization modes, amplitude ratio $r\sim {10}^{-3}$	LIGO/Virgo/KAGRA
Black hole shadows	Size correction $\Delta r\sim 0.1r_g$	Event Horizon Telescope
Early universe	Modified inflation power spectrum	CMB polarization measurements
Solar system tests	Additional perihelion precession of Mercury	Precise planetary orbit measurements

12. Interpretation of the Essence of Energy in Xuan-Liang Theory

12.1. Dilemmas of the Energy Concept in Classical Physics

Energy is one of the most fundamental and mysterious concepts in physics. Classical physics defines energy as the capacity to do work, but it faces conceptual dilemmas in several aspects:

1.

Ambiguity of energy carriers: In what form does energy “flow”? Is it matter, fields, or some more fundamental entity?

2.

Transformation of energy forms: What is the mechanism behind conversions between mechanical, thermal, electromagnetic, chemical, and other forms of energy?

3.

Energy quantization: In quantum mechanics, energy is discrete, but classical theory treats it as continuous.

4.

Relation between energy and information: Does energy transfer necessarily involve information transfer?

12.2. A New Perspective on the Essence of Energy from Xuan-Liang Theory

Xuan-Liang theory offers a novel perspective on the essence of energy: energy is a manifestation of the Xuan-Liang field; energy transfer is essentially the propagation of Xuan-Liang flow.

Definition 9

(Xuan-Liang–energy correspondence principle). In Xuan-Liang theory, the energy density u is related to the Xuan-Liang density ${\rho }_X$ by:

$$u={\displaystyle \frac{1}{2}}{\rho }_0{c}^2+{\displaystyle \frac{1}{2}}{\rho }_0{v}^2+{\kappa }_X{\rho }_X$$

(76)

where:

• First term: rest energy density ${\displaystyle \frac{1}{2}}{\rho }_0{c}^2$
• Second term: kinetic energy density ${\displaystyle \frac{1}{2}}{\rho }_0{v}^2$
• Third term: Xuan-Liang field energy density ${\kappa }_X{\rho }_X$, with ${\kappa }_X$ a conversion coefficient

12.3. Xuan-Liang Description of Explosion Processes

Consider an explosion: an object of mass M fragments and emits photons over a time $\Delta t$.

Theorem 11

(Xuan-Liang conservation in explosions). Explosions satisfy Xuan-Liang conservation:

$$X_{\mathrm{before}}=X_{\mathrm{after}}+X_{\mathrm{radiation}}$$

(77)

where:

$$\begin{array}{cc}\hfill X_{\mathrm{before}}& ={\displaystyle \frac{1}{3}}M{v}_0^3(\mathrm{initial}\mathrm{Xuan}-\mathrm{Liang})\hfill \end{array}$$

(78)

$$\begin{array}{cc}\hfill X_{\mathrm{after}}& =\sum _{i=1}^N{\displaystyle \frac{1}{3}}{m}_i{v}_i^3(\mathrm{fragment}\mathrm{Xuan}-\mathrm{Liang})\hfill \end{array}$$

(79)

$$\begin{array}{cc}\hfill X_{\mathrm{radiation}}& ={\displaystyle \frac{1}{3c}}\sum _{j=1}^M{E}_{\gamma ,j}{v}_{\gamma ,j}^2(\mathrm{radiation}\mathrm{Xuan}-\mathrm{Liang})\hfill \end{array}$$

(80)

$E_{\gamma ,j}$ is the energy of the j-th photon, $v_{\gamma ,j}=c$ is the photon speed.

Proof. 

According to the Xuan-Liang definition (1), a stationary object has zero Xuan-Liang. During an explosion, rest mass is converted into kinetic and radiative energy:

$$M{c}^2=\sum _{i=1}^N{m}_i{c}^2{\gamma }_i+\sum _{j=1}^M{E}_{\gamma ,j}$$

In Xuan-Liang theory, this process corresponds to the generation and propagation of Xuan-Liang. The initial Xuan-Liang $X_{\mathrm{before}}$ transforms into the motion Xuan-Liang of fragments and the radiation Xuan-Liang of the photon field.

The radiation Xuan-Liang requires a relativistic modified definition because photons have no rest mass:

$$X_{\gamma }={\displaystyle \frac{1}{3c}}{E}_{\gamma }{v}_{\gamma }^2={\displaystyle \frac{1}{3}}{E}_{\gamma }c$$

since $v_{\gamma }=c$. □

12.4. Xuan-Liang Transfer Mechanism in Collisions

Consider an elastic collision between two bodies of masses $m_1$, $m_2$, initial velocities ${\overrightarrow{v}}_1$, ${\overrightarrow{v}}_2$, and collision duration $\Delta t$.

Theorem 12

(Xuan-Liang exchange during collisions). The rate of Xuan-Liang exchange during a collision is:

$${\displaystyle \frac{dX}{dt}}={\displaystyle \frac{1}{3}}\left(m_1{\displaystyle \frac{d}{dt}}\left(v_1^3\right)+m_2{\displaystyle \frac{d}{dt}}\left(v_2^3\right)\right)$$

(81)

In the instantaneous collision approximation, the Xuan-Liang change is:

$$\Delta X={\displaystyle \frac{1}{3}}\left[m_1(v_1^{\prime 3}-v_1^3)+m_2(v_2^{\prime 3}-v_2^3)\right]$$

(82)

where $v_i^{\prime }$ are the post-collision velocities.

12.5. Xuan-Liang Fluid: A New Picture of Energy Flow

Definition 10

(Xuan-Liang fluid equations). Describing the Xuan-Liang field as a special fluid, it satisfies the following system:

$$\begin{array}{ccc}\hfill {\displaystyle \frac{\partial {\rho }_X}{\partial t}}+\nabla ·{\overrightarrow{J}}_X& ={\sigma }_X\hfill & \hfill (\mathrm{Xuan}-\mathrm{Liang}\mathrm{continuity}\mathrm{equation})\end{array}$$

(83)

$$\begin{array}{ccc}\hfill {\displaystyle \frac{\partial {\overrightarrow{J}}_X}{\partial t}}+\nabla ·(\overrightarrow{v}\otimes {\overrightarrow{J}}_X)& =-\nabla P_X+{\overrightarrow{f}}_X\hfill & \hfill (\mathrm{Xuan}-\mathrm{Liang}\mathrm{motion}\mathrm{equation})\end{array}$$

(84)

$$\begin{array}{ccc}\hfill {\displaystyle \frac{\partial u_X}{\partial t}}+\nabla ·{\overrightarrow{S}}_X& ={ϵ}_X\hfill & \hfill (\mathrm{Xuan}-\mathrm{Liang}\mathrm{energy}\mathrm{equation})\end{array}$$

(85)

where:

• ${\rho }_X$: Xuan-Liang density, ${\rho }_X={\displaystyle \frac{1}{3c}}{\rho }_0{\gamma }^3{v}^2$
• ${\overrightarrow{J}}_X={\rho }_X\overrightarrow{v}$: Xuan-Liang current density
• $P_X$: Xuan-Liang pressure, $P_X={\displaystyle \frac{1}{3}}{\rho }_X{v}^2$
• $u_X$: Xuan-Liang energy density, $u_X={\displaystyle \frac{1}{2}}\mathbb{X}\wedge ★\mathbb{X}$
• ${\overrightarrow{S}}_X$: Xuan-Liang energy flux density, ${\overrightarrow{S}}_X=u_X\overrightarrow{v}$
• ${\sigma }_X,{\overrightarrow{f}}_X,{ϵ}_X$: source terms

12.6. Xuan-Liang Description of Heat Conduction

Heat conduction is essentially energy diffusion. In the Xuan-Liang theory framework, this can be interpreted as diffusion of the Xuan-Liang field.

Theorem 13

(Xuan-Liang interpretation of Fourier heat conduction). Classical Fourier law:

$$\overrightarrow{q}=-k\nabla T$$

where $\overrightarrow{q}$ is heat flux density, k is thermal conductivity.

In Xuan-Liang theory, the heat flux density is proportional to the Xuan-Liang gradient:

$${\overrightarrow{J}}_X=-D_X\nabla {\rho }_X$$

(86)

where $D_X$ is the Xuan-Liang diffusion coefficient.

The relation between Xuan-Liang density and temperature is:

$${\rho }_X\left(T\right)={\displaystyle \frac{1}{3c}}\rho c_{v}T\langle v_{\mathrm{th}}^2\rangle$$

where $c_v$ is the specific heat capacity, and $\langle v_{\mathrm{th}}^2\rangle$ is the mean square thermal velocity.

12.7. Xuan-Liang Description of Quantum Processes

At the quantum level, Xuan-Liang appears as an operator.

Definition 11

(Quantum Xuan-Liang operator). For a non-relativistic quantum system, the Xuan-Liang operator is defined as:

$$\widehat{X}={\displaystyle \frac{1}{3}}m{\widehat{v}}^3={\displaystyle \frac{1}{3m^2}}{\widehat{p}}^3$$

(87)

In the Heisenberg picture, the evolution equation for the Xuan-Liang operator is:

$$i\hslash {\displaystyle \frac{d\widehat{X}}{dt}}=[\widehat{X},\widehat{H}]$$

(88)

where $\widehat{H}$ is the Hamiltonian operator.

12.8. Energy–Xuan-Liang Unified Theory

Based on the above analysis, we propose the fundamental principles of the energy–Xuan-Liang unified theory:

Proposition 3

(Energy–Xuan-Liang unification principle).

1.

Essence principle: Energy is a manifestation of the Xuan-Liang field; all energy processes can be reduced to the dynamics of the Xuan-Liang field.

2.

Conservation principle: Total Xuan-Liang is conserved; energy conservation is a special case of Xuan-Liang conservation.

3.

Propagation principle: Energy transfer is essentially the propagation of Xuan-Liang flow.

4.

Quantum principle: At the quantum level, energy quanta correspond to Xuan-Liang quanta.

Figure 11. Conceptual diagram of Xuan-Liang theory unifying various energy forms

Figure 11. Conceptual diagram of Xuan-Liang theory unifying various energy forms

This chapter provides a deeper physical foundation and philosophical interpretation for Xuan-Liang theory, linking it to the fundamental physical problem of the essence of energy, thereby enhancing the theory’s internal consistency and explanatory power.

13. Conclusion and Outlook

13.1. Main Conclusions

This paper systematically establishes the theoretical framework of Xuan-Liang unified field theory. The main conclusions are:

1.

Rigorous mathematical foundation: The definition of the Xuan-Liang differential form has been revised, establishing a rigorous differential geometric formulation.

2.

Unified field equations: The differential and integral forms of the unified equations have been derived, and their equivalence proved.

3.

Self-consistent reduction: Systematic proofs of the reduction to general relativity, Newtonian gravity, and cosmological equations are given, ensuring consistency of the reduction.

4.

Wide applicability: Successfully applied to problems in dark matter, dark energy, black hole thermodynamics, and the early universe.

5.

Quantization prospects: A complete theoretical framework including quantum effects has been constructed.

6.

Experimental testability: A series of testable predictions and a detailed experimental roadmap are presented.

13.2. Theoretical Advantages

Xuan-Liang unified field theory offers the following advantages over existing theories:

Table 9. Advantages of Xuan-Liang theory compared to others

Theoretical feature	Xuan-Liang theory	Traditional theories
Unification	Unifies gravity and matter fields	Gravity and matter separate
Self-consistency	Equations form consistent, natural reduction	Requires artificial matching
Flexibility	Form adjustable to problem	Fixed form
Testability	Multi-field testable predictions	Difficult to test (e.g., string theory)
Quantization	Possibly renormalizable, UV complete	General relativity non-renormalizable
Dark matter explanation	Geometric effect, no new particles	Requires dark matter particles
Dark energy explanation	Dynamic field, resolves fine-tuning	Cosmological constant problem

Table 9. Advantages of Xuan-Liang theory compared to others

Theoretical feature	Xuan-Liang theory	Traditional theories
Unification	Unifies gravity and matter fields	Gravity and matter separate
Self-consistency	Equations form consistent, natural reduction	Requires artificial matching
Flexibility	Form adjustable to problem	Fixed form
Testability	Multi-field testable predictions	Difficult to test (e.g., string theory)
Quantization	Possibly renormalizable, UV complete	General relativity non-renormalizable
Dark matter explanation	Geometric effect, no new particles	Requires dark matter particles
Dark energy explanation	Dynamic field, resolves fine-tuning	Cosmological constant problem

13.3. Future Research Directions

Future research can be pursued in the following areas:

1.

Mathematical refinement:

• Develop a classification theory for Xuan-Liang differential forms
• Study exact solutions of the unified equations
• Establish a theory of Xuan-Liang topological invariants

2.

Physical extensions:

• Coupling mechanism of Xuan-Liang with the Standard Model
• Applications of the Xuan-Liang field in condensed matter physics
• Xuan-Liang information theory and the black hole information paradox

3.

Observational tests:

• Design experiments specifically to test Xuan-Liang effects
• Analyze next-generation astronomical data
• Laboratory precision measurements

4.

Quantum gravity:

• Full quantization of Xuan-Liang theory
• Comparative studies with string theory and loop quantum gravity
• Realization of the holographic principle via Xuan-Liang

13.4. Philosophical Significance

Xuan-Liang unified field theory is not only a physical theory but also carries profound philosophical significance:

1.

Unity: Embodies the principle of unity in nature, the deep connection between matter and geometry.

2.

Hierarchy: Reflects the multi-layered nature of physical reality, with different descriptions at different scales.

3.

Dialectics: The dialectical unity of continuity and discreteness, locality and non-locality, determinism and randomness.

4.

Eastern philosophy: Reflects the wisdom of “the mystery of mysteries, the gate of all wonders” in Eastern philosophy.

Xuan-Liang unified field theory provides a new framework for understanding the fundamental laws of the universe. It is not only an elegant mathematical construct but also a bridge connecting human cognition with the essence of nature. Through further mathematical refinement and experimental tests, Xuan-Liang theory has the potential to become a new paradigm unifying gravity and quantum physics, driving revolutionary developments in 21st-century fundamental physics.