Physics as the Solution to an Optimization Problem on Entropy

1. Introduction

The canonical formalism of quantum mechanics (QM) is based on five principal axioms[1,2]:

QM Axiom 1 of 5

State Space: Each physical system corresponds to a complex Hilbert space, with the system’s state represented by a ray in this space.

QM Axiom 2 of 5

Observables: Physical observables correspond to Hermitian operators within the Hilbert space.

QM Axiom 3 of 5

Dynamics: The time evolution of a quantum system is dictated by the Schrödinger equation, where the Hamiltonian operator signifies the system’s total energy.

QM Axiom 4 of 5

Measurement: The act of measuring an observable results in the system’s transition to an eigenstate of the associated operator, with the measurement value being one of the eigenvalues.

QM Axiom 5 of 5

Probability Interpretation: The likelihood of a specific measurement outcome is determined by the squared magnitude of the state vector’s projection onto the relevant eigenstate.

Contrastingly, statistical mechanics (SM), the other statistical pillar of physics, derives its probability measure through entropy maximization, constrained by the following expression:

SM Constraint 1 of 1:

Average Energy Constraint: The average of energy measurements of a system at thermodynamic equilibrium converge to a specific value ($\overline{E}$):

$$\begin{array}{c}\hfill \overline{E}=\sum _i{\rho }_i{E}_i\end{array}$$

(1)

To maximize entropy while satisfying this constraint, the theory uses a Lagrange multiplier approach.

Definition 1 

(Fundamental Lagrange Multiplier Equation of SM).

$$\begin{array}{c}\hfill \mathcal{L}=\underset{\mathit{Boltzmann}\mathit{entropy}}{\underbrace{-k_B{\sum }_i{\rho }_{i}ln{\rho }_i}}+\underset{NormalizationConstraint}{\underbrace{\lambda \left(1-{\sum }_i{\rho }_i\right)}}+\underset{AverageEnergyConstraint}{\underbrace{\beta \left(\overline{E}-{\sum }_i{\rho }_i{E}_i\right)}}\end{array}$$

(2)

where λ and β are the Lagrange multipliers.

Theorem 1 

(Gibbs Measure). The solution to the Lagrange multiplier equation of SM is the Gibbs measure.

$$\begin{array}{c}\hfill {\rho }_i=\underset{\mathit{Microcanonical}\mathit{Ensemble}}{\underbrace{{\displaystyle \frac{1}{{\sum }_{i}exp(-\beta E_i)}}}}exp(-\beta E_i)\end{array}$$

(3)

Proof. 

This is an well-known result by E. T. Jaynes [3,4]. As a convenience, we replicate the proof in Appendix A.    □

As evident from E. T. Jaynes’ methodological innovation, SM relies on a single constraint related to the nature of the measurements under consideration, which allows the formulation of an optimization problem sufficient to derive the relevant probability measure. This is an exceptionally parsimonious formulation of a physical theory.

We propose a generalization of E. T. Jaynes’ approach to the realms of Quantum Mechanics (QM), relativistic Quantum Mechanics (RQM) and Quantum Gravity (QG). For each domain, we will introduce a single constraint related to measurements, formulate a corresponding entropy maximization problem, and present a main theorem that encapsulates the theory. This formulation reduces fundamental physics to its most parsimonious expression, deriving the core theories as optimal solutions to a well-defined entropy maximization problem.

1.1. Quantum Mechanics

To reformulate QM as the solution to an entropy maximization problem, we propose the following constraint:

QM Constraint 1 of 1

Vanishing Complex-Phase: Quantum measurements admit a vanishing complex phase. The constraint is:

$$\begin{array}{c}\hfill 0=\mathrm{tr}\sum _i{\rho }_i\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\end{array}$$

(4)

Here, the matrix representation engenders the complex phase, and the trace will cause it to vanish under measurement.

It associates to the follow equation:

Definition 2 

(Fundamental Lagrange Multiplier Equation of QM).

$$\begin{array}{c}\hfill \mathcal{L}=\underset{RelativeShannonEntropy}{\underbrace{-{\sum }_i{\rho }_{i}ln{\displaystyle \frac{{\rho }_i}{p_i}}}}+\underset{NormalizationConstraint}{\underbrace{\lambda \left(1-{\sum }_i{\rho }_i\right)}}+\underset{VanishingComplex\text{-}Phase}{\underbrace{\tau \left(-\mathrm{tr}{\sum }_i{\rho }_i\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\right)}}\end{array}$$

(5)

where λ and τ are the Lagrange multipliers.

The relative Shannon entropy[5,6] is utilized because we are solving for the least biased theory that connects an initial preparation p to its final measurement $\rho$.

Theorem 2. 

The least biased probability measure that connects an initial preparation p to its final measurement ρ, under the constraint of the vanishing complex-phase, is:

$$\begin{array}{c}\hfill {\rho }_i=\underset{\mathit{Unitarily}\mathit{Invariant}\mathit{Ensemble}}{\underbrace{{\displaystyle \frac{1}{{\sum }_i{p}_i\parallel exp(-it{E}_i/\hslash )\parallel }}}}\underset{\mathit{Born}\mathit{Rule}}{\underbrace{\parallel exp(-it{E}_i/\hslash )\parallel }}\underset{\mathit{Initial}\mathit{Preparation}}{\underbrace{p_i}}\end{array}$$

(6)

where we have defined $\tau =t/\hslash$ (analogous to $\beta =1/\left(k_{B}T\right)$ in SM).

The proof of this theorem will be presented in the results section. We will show that this solution entails the five axioms of QM, which are now promoted to theorems, yielding a parsimonious formulation of QM.

1.2. Relativistic Quantum Mechanics

Before we can discuss RQM, we first need to introduce a notation. Let $\mathbf{u}=a+\mathbf{x}+\mathbf{f}+\mathbf{v}+\mathbf{b}$, where a is a scalar, $\mathbf{x}$ is a vector, $\mathbf{f}$ is a bivector, $\mathbf{v}$ is a pseudo-vector and $\mathbf{b}$ is a pseudo-scalar, be a multivector of the geometric algebra $\mathrm{GA}(3,1)$, and let $\mathbf{M}$ be its matrix representation. Then, the fundamental constraint is:

QG Constraint 1 of 1

Vanishing Relativistic Phase: Our formulation of RQM is based around a vanishing phase spanning the ${\mathrm{Spin}}^{\mathrm{c}}(3,1)$ group. The constraint is:

$$\begin{array}{c}\hfill 0=\mathrm{tr}{\displaystyle \frac{1}{2}}\sum _i{\rho }_i{\mathbf{M}}_i\end{array}$$

(7)

where ${\mathbf{M}}_i$ is the matrix representation of the multivector ${\mathbf{u}}_i={\mathbf{f}}_i+{\mathbf{b}}_i$ of $\mathrm{GA}(3,1)$. Using the real Majorana representation of the gamma matrices, the representation is as follows:

$$\begin{array}{cc}\hfill \mathbf{M}& =\left[\begin{array}{cccc}{f}_{02}& b-f_{13}& -f_{01}+f_{12}& f_{03}+f_{23}\\ -b+f_{13}& f_{02}& f_{03}+f_{23}& f_{01}-f_{12}\\ -f_{01}-f_{12}& f_{03}-f_{23}& -f_{02}& -b-f_{13}\\ f_{03}-f_{23}& f_{01}+f_{12}& b+f_{13}& -f_{02}\end{array}\right]\hfill \end{array}$$

(8)

$$\begin{array}{cc}\hfill & \cong f_{01}{\gamma }_0{\gamma }_1+f_{02}{\gamma }_0{\gamma }_2+f_{03}{\gamma }_0{\gamma }_3+f_{12}{\gamma }_1{\gamma }_2+f_{13}{\gamma }_1{\gamma }_3+f_{23}{\gamma }_2{\gamma }_3+b{\gamma }_0{\gamma }_1{\gamma }_2{\gamma }_3\hfill \end{array}$$

(9)

Similarly to the QM case, here the matrix representation engenders a ${\mathrm{Spin}}^{\mathrm{c}}(3,1)$-phase and the trace will cause it to vanish under measurement.

The Lagrange multiplier equation is as follows:

Definition 3 

(Fundamental Lagrange Multiplier Equation of RQM).

$$\begin{array}{c}\hfill \mathcal{L}=\underset{RelativeShannonEntropy}{\underbrace{-{\sum }_i{\rho }_{i}ln{\displaystyle \frac{{\rho }_i}{p_i}}}}+\underset{NormalizationConstraint}{\underbrace{\lambda \left(1-{\sum }_i{\rho }_i\right)}}+\underset{VanishingSpi{n}^c(3,1)Phase}{\underbrace{\zeta \left(-\mathrm{tr}{\displaystyle \frac{1}{2}}{\sum }_i{\rho }_i{\mathbf{M}}_i\right)}}\end{array}$$

(10)

where λ and ζ are the Lagrange multipliers.

Theorem 3. 

The least biased probability measure that connects an initial preparation $p_i$ to its final measurement ${\rho }_i$, under the constraint of the vanishing relativistic phase, is:

$$\begin{array}{c}\hfill \rho \left(q\right)=\underset{{\mathit{Spin}}^c(3,1)\mathit{Invariant}\mathit{Ensemble}}{\underbrace{{\displaystyle \frac{1}{{\sum }_i{p}_{i}detexp(-\zeta \frac{1}{2}{\mathbf{M}}_i)}}}}\underset{{\mathit{Spin}}^c(3,1)\mathit{Born}\mathit{Rule}}{\underbrace{detexp(-\zeta {\displaystyle \frac{1}{2}}{\mathbf{M}}_i)}}\underset{\mathit{Initial}\mathit{Preparation}}{\underbrace{p_i}}\end{array}$$

(11)

The proof of this theorem is presented in the results section.

In the results section, we demonstrate that the solution to this optimization problem provides a foundation for a quantum mechanical theory in 3+1-dimensional spacetime. The parameter $\zeta$ emerges as a the generator of boosts, rotations, and phase transformations. This single parameter accounts for all change-of-basis transformations an observer can perform in spacetime before measurement. Furthermore, we show that the Dirac current resulting from this formulation exhibits invariance under $\mathrm{SU}\left(3\right)\times \mathrm{SU}\left(2\right)\times \mathrm{U}\left(1\right)$ gauge symmetries, aligning with the fundamental symmetries of the Standard Model.

1.3. Quantum Gravity

Our approach to RQM extends naturally to QG. While RQM utilizes only a portion of the solution space, QG leverages its full generality. The key innovation is the introduction of a double-copy mechanism applied to the Dirac current: the probability measure multiply two Dirac currents to obtain a metric tensor expectation value. This mechanism allows us to construct the metric tensor as an observable from basis vectors, thereby establishing a direct link between the probability measure, metric measurements, and the geometric structure of spacetime.

1.4. Dimensional Obstructions

We conclude the results section with a series of theorems demonstrating that, except for SM (no vanishing phase) and QM (vanishing $\mathrm{U}\left(1\right)$ phase), the entropy maximization technique yields a solution only in 3+1-dimensional spacetime (vanishing ${\mathrm{Spin}}^{\mathrm{c}}(3,1)$ phase). In all other dimensional configurations, various obstructions are encountered. These findings suggest an intriguing connection between the entropy maximization approach and the specific dimensionality of our universe, the implications of which are discussed.

2. Results

2.1. Quantum Mechanics

In statistical mechanics, the founding observation is that energy measurements of a thermally equilibrated system tend towards an average value. Comparatively, in QM, the founding observation involves the interplay between the systematic elimination of complex phases in measurement outcomes and the presence of interference effects in repeated measurement outcomes. To represent this observation, we introduce the Vanishing $\mathrm{U}\left(1\right)$-Phase Anti-Constraint:

$$\begin{array}{c}\hfill 0=\mathrm{tr}\sum _i{\rho }_i\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\end{array}$$

(12)

where $E_i$ are scalar-valued. The usage of the matrix generates a $\mathrm{U}\left(1\right)$ phase, and the trace causes it to vanish under measurements.

At first glance, this expression may seem to reduce to a tautology equating zero with zero, suggesting it imposes no restriction on energy measurements. However, this appearance is deceptive. Unlike a conventional constraint that limits the solution space, this expression serves as a formal device to expand it, allowing for the incorporation of complex phases into the probability measure. The expression’s role in broadening, rather than restricting, the solution space leads to its designation as an ”anti-constraint.”

In general, usage of anti-constraints expand classical probability measures into larger domains, such as quantum probabilities.

Its significance will become evident upon the completion of the optimization problem. For the moment, this expression can be conceptualized as an ansatz that, when incorporated as an anti-constraint within an entropy-maximization problem, resolves into the axioms of quantum mechanics.

Our next procedural step involves solving the corresponding Lagrange multiplier equation, mirroring the methodology employed in statistical mechanics by E. T. Jaynes. We utilize the relative Shannon entropy because we wish to solve for the least biased probability measure that connects an initial preparation $p_i$ to its final measurement ${\rho }_i$. For that, we deploy the following Lagrange multiplier equation:

$$\begin{array}{c}\hfill \mathcal{L}=\underset{RelativeShannonEntropy}{\underbrace{-{\sum }_i{\rho }_{i}ln{\displaystyle \frac{{\rho }_i}{p_i}}}}+\underset{NormalizationConstraint}{\underbrace{\lambda \left(1-{\sum }_i{\rho }_i\right)}}+\underset{VanishingU\left(1\right)\text{-}Phase}{\underbrace{\tau \left(tr{\sum }_i{\rho }_i\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\right)}}\end{array}$$

(13)

Where $\lambda$ and $\tau$ are the Lagrange multipliers.

We solve the maximization problem as follows:

$$\begin{array}{cc}\hfill 0& ={\displaystyle \frac{\partial \mathcal{L}[{\rho }_1,\cdots ,{\rho }_n]}{\partial {\rho }_i}}\hfill \end{array}$$

(14)

$$\begin{array}{cc}\hfill & =-ln{\displaystyle \frac{{\rho }_i}{p_i}}-p_i-\lambda -\tau tr\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\hfill \end{array}$$

(15)

$$\begin{array}{cc}\hfill & =ln{\displaystyle \frac{{\rho }_i}{p_i}}+p_i+\lambda -\tau tr\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\hfill \end{array}$$

(16)

$$\begin{array}{cc}\hfill \Rightarrow ln{\displaystyle \frac{{\rho }_i}{p_i}}& =-p_i-\lambda -\tau tr\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\hfill \end{array}$$

(17)

$$\begin{array}{cc}\hfill \Rightarrow {\rho }_i& =p_{i}exp(-p_i-\lambda )exp\left(-\tau tr\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\right)\hfill \end{array}$$

(18)

$$\begin{array}{cc}\hfill & ={\displaystyle \frac{1}{Z\left(\tau \right)}}{p}_{i}exp\left(-\tau tr\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\right)\hfill \end{array}$$

(19)

The partition function is obtained as follows:

$$\begin{array}{cc}\hfill 1& =\sum _i{p}_{i}exp(-p_i-\lambda )exp\left(-\tau tr\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\right)\hfill \end{array}$$

(20)

$$\begin{array}{cc}\hfill \Rightarrow {\left(exp(-p_i-\lambda )\right)}^{-1}& =\sum _i{p}_{i}exp\left(-\tau tr\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\right)\hfill \end{array}$$

(21)

$$\begin{array}{cc}\hfill Z\left(\tau \right)& :=\sum _i{p}_{i}exp\left(-\tau tr\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\right)\hfill \end{array}$$

(22)

Finally, the least biased probability measure that connects an initial preparation $p_i$ to its final measurement ${\rho }_i$, under the constraint of the vanishing U(1) phase, is:

$$\begin{array}{c}\hfill {\rho }_i={\displaystyle \frac{1}{{\sum }_i{p}_{i}exp\left(-\tau tr\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\right)}}exp\left(-\tau tr\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\right)p_i\end{array}$$

(23)

Though initially unfamiliar, this form effectively establishes a comprehensive formulation of quantum mechanics, as we will demonstrate.

Upon examination, we find that phase elimination is manifestly evident in the probability measure: since the trace evaluates to zero, the probability measure simplifies to classical probabilities, aligning precisely with the Born rule’s exclusion of complex phases:

$$\begin{array}{c}\hfill {\rho }_i={\displaystyle \frac{p_i}{{\sum }_i{p}_i}}\end{array}$$

(24)

However, the significance of this phase elimination extends beyond this mere simplicity. As we will soon see, the partition function Z gains unitary invariance, allowing for the emergence of interference patterns and other quantum characteristics under appropriate basis changes.

We will begin by aligning our results with the conventional quantum mechanical notation. As such, we transform the representation of complex numbers from $\left[\begin{array}{cc}a& -b\\ b& a\end{array}\right]$ to $a+ib$. For instance, the exponential of a complex matrix is:

$$\begin{array}{cc}\hfill exp\left[\begin{array}{cc}a& -b\\ b& a\end{array}\right]& =r\left[\begin{array}{cc}cos\left(b\right)& -sin\left(b\right)\\ sin\left(b\right)& cos\left(b\right)\end{array}\right],\mathrm{where}r=expa\hfill \end{array}$$

(25)

Then, we associate the exponential trace to the complex norm using $exp\mathrm{tr}\mathbf{M}\equiv detexp\mathbf{M}$:

$$\begin{array}{cc}\hfill exp\mathrm{tr}\left[\begin{array}{cc}a& -b\\ b& a\end{array}\right]=detexp\left[\begin{array}{cc}a& -b\\ b& a\end{array}\right]& =r^{2}det\left[\begin{array}{cc}cos\left(b\right)& -sin\left(b\right)\\ sin\left(b\right)& cos\left(b\right)\end{array}\right],\mathrm{where}r=expa\hfill \end{array}$$

(26)

$$\begin{array}{cc}\hfill & =r^2({cos}^2\left(b\right)+{sin}^2\left(b\right))\hfill \end{array}$$

(27)

$$\begin{array}{cc}\hfill & =\parallel r(cos(b)+isin(b\left)\right)\parallel \hfill \end{array}$$

(28)

$$\begin{array}{cc}\hfill & =\parallel rexp\left(ib\right)\parallel \hfill \end{array}$$

(29)

Finally, substituting $\tau =t/\hslash$ analogously to $\beta =1/\left(k_{B}T\right)$, and applying the complex-norm representation to both the numerator and to the denominator, consolidates the Born rule, normalization, and initial prepration into :

$$\begin{array}{cc}\hfill {\rho }_i& =\underset{\mathrm{Unitarily}\mathrm{Invariant}\mathrm{Partition}\mathrm{Function}}{\underbrace{{\displaystyle \frac{1}{{\sum }_i{p}_i\parallel exp(-it{E}_i/\hslash )\parallel }}}}\underset{\mathrm{Born}\mathrm{Rule}}{\underbrace{\parallel exp(-it{E}_i/\hslash )\parallel }}\underset{\mathrm{Initial}\mathrm{Preparation}}{\underbrace{p_i}}\hfill \end{array}$$

(30)

We are now in a position to explore the solution space.

The wavefunction is delineated by decomposing the complex norm into a complex number and its conjugate. It is then visualized as a vector within a complex n-dimensional Hilbert space. The partition function acts as the inner product. This relationship is articulated as follows:

$$\begin{array}{c}\hfill \sum _i{p}_i\parallel exp(-it{E}_i/\hslash )\parallel =Z=\langle \psi |\psi \rangle \end{array}$$

(31)

where

$$\begin{array}{c}\hfill \left[\begin{array}{c}{\psi }_1\left(t\right)\\ ⋮\\ {\psi }_n\left(t\right)\end{array}\right]=\left[\begin{array}{c}exp(-it{E}_1/\hslash )\\ & \ddots \end{array}\right]\begin{array}{cc}\hfill & exp(-it{E}_n/\hslash )\hfill \end{array}\left[\begin{array}{c}{\psi }_1\left(0\right)\\ ⋮\\ {\psi }_n\left(0\right)\end{array}\right]\end{array}$$

(32)

We clarify that $p_i$ represents the probability associated with the initial preparation of the wavefunction, where $p_i=\langle {\psi }_i\left(0\right)|{\psi }_i\left(0\right)\rangle$.

We also note that Z is invariant under unitary transformations.

Let us now investigate how the axioms of quantum mechanics are recovered from this result:

• The entropy maximization procedure inherently normalizes the vectors $|\psi \rangle$ with $1/Z=1/\sqrt{\langle \psi |\psi \rangle }$. This normalization links $|\psi \rangle$ to a unit vector in Hilbert space. Furthermore, as physical states associate to the probability measure, and the probability is defined up to a phase, we conclude that physical states map to Rays within Hilbert space. This demonstrates QM Axiom 1 of 5.
• In Z, an observable must satisfy:

  $$\begin{array}{c}\hfill \overline{O}=\sum _i{p}_i{O}_i\parallel exp(-it{E}_i/\hslash )\parallel \end{array}$$

  (33)
  Since $Z=\langle \psi |\psi \rangle$, then any self-adjoint operator satisfying the condition $\langle \mathbf{O}\psi |\varphi \rangle =\langle \psi |\mathbf{O}\varphi \rangle$ will equate the above equation, simply because $\langle \mathbf{O}\rangle =\langle \psi \left|\mathbf{O}\right|\psi \rangle$. This demonstrates QM Axiom 2 of 5.
• Upon transforming Equation 32 out of its eigenbasis through unitary operations, we find that the energy, $E_i$, typically transforms in the manner of a Hamiltonian operator:

  $$\begin{array}{c}\hfill \left|\psi \right(t)\rangle =exp(-it\mathbf{H}/\hslash \left)\right|\psi \left(0\right)\rangle \end{array}$$

  (34)
  The system’s dynamics emerge from differentiating the solution with respect to the Lagrange multiplier. This is manifested as:

  $$\begin{array}{cc}\hfill {\displaystyle \frac{\partial }{\partial t}}|\psi \left(t\right)\rangle & ={\displaystyle \frac{\partial }{\partial t}}(exp(-it\mathbf{H}/\hslash )|\psi \left(0\right)\rangle )\hfill \end{array}$$

  (35)

  $$\begin{array}{cc}\hfill & =-i\mathbf{H}/\hslash exp(-it\mathbf{H}/\hslash )\left|\psi \right(0)\rangle \hfill \end{array}$$

  (36)

  $$\begin{array}{cc}\hfill & =-i\mathbf{H}/\hslash \left|\psi \right(t)\rangle \hfill \end{array}$$

  (37)

  $$\begin{array}{cc}\hfill \Rightarrow \mathbf{H}\left|\psi \right(t)\rangle & =i\hslash {\displaystyle \frac{\partial }{\partial t}}|\psi \left(t\right)\rangle \hfill \end{array}$$

  (38)

  which is the Schrödinger equation. This demonstrates QM Axiom 3 of 5.
• From Equation 32 it follows that the possible microstates $E_i$ of the system correspond to specific eigenvalues of $\mathbf{H}$. An observation can thus be conceptualized as sampling from $\rho$, with the measured state being the occupied microstate i. Consequently, when a measurement occurs, the system invariably emerges in one of these microstates, which directly corresponds to an eigenstate of $\mathbf{H}$. Measured in the eigenbasis, the probability measure is:

  $$\begin{array}{c}\hfill {\rho }_i\left(t\right)={\displaystyle \frac{1}{\langle \psi |\psi \rangle }}{\left({\psi }_i\left(t\right)\right)}^{†}{\psi }_i\left(t\right).\end{array}$$

  (39)
  In scenarios where the probability measure ${\rho }_i\left(\tau \right)$ is expressed in a basis other than its eigenbasis, the probability $P\left({\lambda }_i\right)$ of obtaining the eigenvalue ${\lambda }_i$ is given as a projection on a eigenstate:

  $$\begin{array}{c}\hfill P\left({\lambda }_i\right)={\left|\langle {\lambda }_i|\psi \rangle \right|}^2\end{array}$$

  (40)
  Here, $|\langle {\lambda }_i|\psi \rangle {|}^2$ signifies the squared magnitude of the amplitude of the state $|\psi \rangle$ when projected onto the eigenstate $|{\lambda }_i\rangle$. As this argument hold for any observables, this demonstrates QM Axiom 4 of 5.
• Finally, since the probability measure (Equation 30) replicates the Born rule, QM Axiom 5 of 5 is also demonstrated.

Revisiting quantum mechanics with this perspective offers a coherent and unified narrative. Specifically, the vanishing U(1) phase constraint (Equation 12) is sufficient to entail the foundations of quantum mechanics (Axiom 1, 2, 3, 4 and 5) through the principle of entropy maximization. Equation 12 becomes the formulation’s new singular foundation, and Axioms 1, 2, 3, 4, and 5 are now promoted to theorems.

2.2. RQM in 2D

In this section, we investigate RQM in 2D. Although all geometric configurations except 3+1D contain obstructions, which will be discussed later in this section, the 2D case provides a valuable starting point before addressing the more complex 3+1D case. In RQM 2D, the fundamental Lagrange Multiplier Equation is:

$$\begin{array}{c}\hfill \mathcal{L}=\underset{\mathrm{Relative}\mathrm{Shannon}\mathrm{Entropy}}{\underbrace{-{\sum }_i{\rho }_{i}ln{\displaystyle \frac{{\rho }_i}{p_i}}}}+\underset{\mathrm{Normalization}\mathrm{Constraint}}{\underbrace{\lambda \left(1-{\sum }_i{\rho }_i\right)}}+\underset{\mathrm{Vanishing}\mathrm{Relativistic}\mathrm{Phase}}{\underbrace{\theta \left(-\mathrm{tr}{\displaystyle \frac{1}{2}}{\sum }_i{\rho }_i{\mathbf{M}}_i\right)}}\end{array}$$

(41)

where $\lambda$ and $\theta$ are the Lagrange multipliers, and where ${\mathbf{M}}_i$ is the matrix representation of a multivector ${\mathbf{b}}_i$ of $\mathrm{GA}\left(2\right)$, where ${\mathbf{b}}_i$ is a pseudo-scalar. In general a multivector $\mathbf{u}=a+\mathbf{x}+\mathbf{b}$ of $\mathrm{GA}\left(2\right)$, where a is a scalar, $\mathbf{x}$ is a vector and $\mathbf{b}$ a pseudo-scalar, is represented as follows:

$$\begin{array}{c}\hfill \left[\begin{array}{cc}a+x& y-b\\ y+b& a-x\end{array}\right]\cong a+x{\sigma }_x+y{\sigma }_y+b{\sigma }_x\wedge {\sigma }_y\end{array}$$

(42)

The basis elements are defined as:

$$\begin{array}{c}\hfill {\sigma }_x=\left[\begin{array}{cc}1& 0\\ 0& -1\end{array}\right],{\sigma }_y=\left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right],{\sigma }_x\wedge {\sigma }_y=\left[\begin{array}{cc}0& -1\\ 1& 0\end{array}\right]\end{array}$$

(43)

If we take $a\to 0,\mathbf{x}\to 0$ then $\mathbf{M}$ reduces as follows:

$$\begin{array}{c}\hfill {\mathbf{u}=a+\mathbf{x}+\mathbf{b}|}_{a\to 0,\mathbf{x}\to 0}=\mathbf{b}\Rightarrow \mathbf{M}=\left[\begin{array}{cc}0& -b\\ b& 0\end{array}\right]\end{array}$$

(44)

The Lagrange multiplier equation can be solved as follows:

$$\begin{array}{cc}\hfill 0& ={\displaystyle \frac{\partial \mathcal{L}[{\rho }_1,\cdots ,{\rho }_n]}{\partial {\rho }_i}}\hfill \end{array}$$

(45)

$$\begin{array}{cc}\hfill & =-ln{\displaystyle \frac{{\rho }_i}{p_i}}-p_i-\lambda -\theta tr{\displaystyle \frac{1}{2}}\left[\begin{array}{cc}0& -b_i\\ b_i& 0\end{array}\right]\hfill \end{array}$$

(46)

$$\begin{array}{cc}\hfill & =ln{\displaystyle \frac{{\rho }_i}{p_i}}+p_i+\lambda +\theta tr{\displaystyle \frac{1}{2}}\left[\begin{array}{cc}0& -b_i\\ b_i& 0\end{array}\right]\hfill \end{array}$$

(47)

$$\begin{array}{cc}\hfill \Rightarrow ln{\displaystyle \frac{{\rho }_i}{p_i}}& =-p_i-\lambda -\theta tr{\displaystyle \frac{1}{2}}\left[\begin{array}{cc}0& -b_i\\ b_i& 0\end{array}\right]\hfill \end{array}$$

(48)

$$\begin{array}{cc}\hfill \Rightarrow {\rho }_i& =p_{i}exp(-p_i-\lambda )exp\left(-\theta tr{\displaystyle \frac{1}{2}}\left[\begin{array}{cc}0& -b_i\\ b_i& 0\end{array}\right]\right)\hfill \end{array}$$

(49)

$$\begin{array}{cc}\hfill & ={\displaystyle \frac{1}{Z\left(\theta \right)}}{p}_{i}exp\left(-\theta tr{\displaystyle \frac{1}{2}}\left[\begin{array}{cc}0& -b_i\\ b_i& 0\end{array}\right]\right)\hfill \end{array}$$

(50)

The partition function $Z\left(\theta \right)$, serving as a normalization constant, is determined as follows:

$$\begin{array}{cc}\hfill 1& =\sum _i{p}_{i}exp(-p_i-\lambda )exp\left(-\theta tr{\displaystyle \frac{1}{2}}\left[\begin{array}{cc}0& -b_i\\ b_i& 0\end{array}\right]\right)\hfill \end{array}$$

(51)

$$\begin{array}{cc}\hfill \Rightarrow {\left(exp(-p_i-\lambda )\right)}^{-1}& =\sum _i{p}_{i}exp\left(-\theta tr{\displaystyle \frac{1}{2}}\left[\begin{array}{cc}0& -b_i\\ b_i& 0\end{array}\right]\right)\hfill \end{array}$$

(52)

$$\begin{array}{cc}\hfill Z\left(\theta \right)& :=\sum _i{p}_{i}exp\left(-\theta tr{\displaystyle \frac{1}{2}}\left[\begin{array}{cc}0& -b_i\\ b_i& 0\end{array}\right]\right)\hfill \end{array}$$

(53)

Consequently, the least biased probability measure that connects an initial preparation $p_i$ to a final measurement ${\rho }_i$, under the constraint of the vanishing relativistic phase in 2D is:

$$\begin{array}{c}\hfill {\rho }_i=\underset{\mathrm{Spin}\left(2\right)\mathrm{Invariant}\mathrm{Ensemble}}{\underbrace{{\displaystyle \frac{1}{{\sum }_i{p}_{i}detexp\left(-\frac{1}{2}\theta \left[\begin{array}{cc}0& -b_i\\ b_i& 0\end{array}\right]\right)}}}}\underset{\mathrm{Spin}\left(2\right)}{\underbrace{detexp\left(-{\displaystyle \frac{1}{2}}\theta \left[\begin{array}{cc}0& -b_i\\ b_i& 0\end{array}\right]\right)}}\mathrm{Born}\mathrm{Rule}\underset{\mathrm{Initial}\mathrm{Preparation}}{\underbrace{p_i}}\end{array}$$

(54)

where $detexpM=exp\mathrm{tr}M$.

In 2D, the Lagrange multiplier $\theta$ correspond to an angle of rotation, and in 1+1D it would correspond to the rapidity $\zeta$:

$$\begin{array}{cccccc}\hfill & 2\mathrm{D}:\hfill & \hfill & exp\theta \left[\begin{array}{cc}0& -1\\ 1& 0\end{array}\right]=\left[\begin{array}{cc}cos\theta & -sin\theta \\ sin\theta & cos\theta \end{array}\right]\hfill & \hfill & \theta \mathrm{is}\mathrm{the}\mathrm{angle}\mathrm{of}\mathrm{rotation}\hfill \end{array}$$

(55)

$$\begin{array}{cccccc}\hfill & 1+1\mathrm{D}:\hfill & \hfill & exp\zeta \left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right]=\left[\begin{array}{cc}cosh\zeta & sinh\zeta \\ sinh\zeta & cosh\zeta \end{array}\right]\hfill & \hfill & \zeta \mathrm{is}\mathrm{the}\mathrm{rapidity}\hfill \end{array}$$

(56)

The 2D solution may appear equivalent to the QM case because they are related by an isomorphism $\mathrm{Spin}\left(2\right)\cong \mathrm{SO}\left(2\right)\cong \mathrm{U}\left(1\right)$ and under the replacement $\theta \to \tau$. However, an isomorphism does not mean identical, and in Spin(2) we gain extra structures related to a relativistic description, which are not available in the QM case.

To investigate the solution in more detail, we introduce the multivector conjugate, also known as the Clifford conjugate, which generalizes the concept of complex conjugation to multivectors.

Definition 4 

(Multivector conjugate (a.k.a Clifford conjugate)). Let $\mathbf{u}=a+\mathbf{x}+\mathbf{b}$ be a multi-vector of the geometric algebra over the reals in two dimensions $\mathrm{GA}\left(2\right)$. The multivector conjugate is defined as:

$$\begin{array}{c}\hfill {\mathbf{u}}^{‡}=a-\mathbf{x}-\mathbf{b}\end{array}$$

(57)

The determinant of the matrix representation of a multivector can be expressed as a self-product:

Theorem 4 

(Determinant as a Multivector Self-Product).

$$\begin{array}{c}\hfill {\mathbf{u}}^{‡}\mathbf{u}=det\mathbf{M}\end{array}$$

(58)

Proof. 

Let $\mathbf{u}=a+x{\sigma }_x+y{\sigma }_y+b{\sigma }_x\wedge {\sigma }_y$, and let $\mathbf{M}$ be its matrix representation $\left[\begin{array}{cc}a+x& y-b\\ y+b& a-x\end{array}\right]$. Then:

$$\begin{array}{cc}\hfill 1:& {\mathbf{u}}^{‡}\mathbf{u}\hfill \end{array}$$

(59)

$$\begin{array}{cc}\hfill & ={(a+x{\sigma }_x+y{\sigma }_y+b{\sigma }_x\wedge {\sigma }_y)}^{‡}(a+x{\sigma }_x+y{\sigma }_y+b{\sigma }_x\wedge {\sigma }_y)\hfill \end{array}$$

(60)

$$\begin{array}{cc}\hfill & =(a-x{\sigma }_x-y{\sigma }_y-b{\sigma }_x\wedge {\sigma }_y)(a+x{\sigma }_x+y{\sigma }_y+b{\sigma }_x\wedge {\sigma }_y)\hfill \end{array}$$

(61)

$$\begin{array}{cc}\hfill & =a^2-x^2-y^2+b^2\hfill \end{array}$$

(62)

$$\begin{array}{cc}\hfill 2:& det\mathbf{M}\hfill \end{array}$$

(63)

$$\begin{array}{cc}\hfill & =det\left[\begin{array}{cc}a+x& y-b\\ y+b& a-x\end{array}\right]\hfill \end{array}$$

(64)

$$\begin{array}{cc}\hfill & =(a+x)(a-x)-(y-b)(y+b)\hfill \end{array}$$

(65)

$$\begin{array}{cc}\hfill & =a^2-x^2-y^2+b^2\hfill \end{array}$$

(66)

   □

Building upon the concept of the multivector conjugate, we introduce the multivector conjugate transpose, which serves as an extension of the Hermitian conjugate to the domain of multivectors.

Definition 5 

(Multivector Conjugate Transpose). Let $\left|V\right.〉{\rangle \in \left(\mathrm{GA}\left(2\right)\right)}^n$:

$$\begin{array}{c}\hfill \left|V\right.〉\rangle =\left[\begin{array}{c}{a}_1+{\mathbf{x}}_1+{\mathbf{b}}_1\\ ⋮\\ a_n+{\mathbf{x}}_n+{\mathbf{b}}_n\end{array}\right]\end{array}$$

(67)

The multivector conjugate transpose of $\left|V\right.〉\rangle$ is defined as first taking the transpose and then the element-wise multivector conjugate:

$$\begin{array}{c}\hfill \langle \left.〈V\right|=\left[\begin{array}{ccc}{a}_1-{\mathbf{x}}_1-{\mathbf{b}}_1& \cdots & a_n-{\mathbf{x}}_n-{\mathbf{b}}_n\end{array}\right]\end{array}$$

(68)

Definition 6 

(Bilinear Form). Let $\left|V\right.〉\rangle$ and $\left|W\right.〉\rangle$ be two vectors valued in $\mathrm{GA}\left(2\right)$. We introduce the following bilinear form:

$$\begin{array}{c}\hfill \langle 〈V|W〉\rangle =(a_1-{\mathbf{x}}_1-{\mathbf{b}}_1)(a_1+{\mathbf{x}}_1+{\mathbf{b}}_1)+\cdots (a_n-{\mathbf{x}}_n-{\mathbf{b}}_n)(a_n+{\mathbf{x}}_n+{\mathbf{b}}_n)\end{array}$$

(69)

Theorem 5 

(Inner Product). Restricted to the even sub-algebra of $\mathrm{GA}\left(2\right)$, the bilinear form is an inner product.

Proof. 

$$\begin{array}{cc}\hfill {\langle 〈V|W〉\rangle }_{\mathbf{x}\to 0}& =(a_1-{\mathbf{b}}_1)(a_1+{\mathbf{b}}_1)+\cdots (a_n-{\mathbf{b}}_n)(a_n+{\mathbf{b}}_n)\hfill \end{array}$$

(70)

This is isomorphic to the inner product of a complex Hilbert space, with the identification $i\cong {\sigma }_x\wedge {\sigma }_y$.    □

Definition 7 

(Spin(2)-valuedWavefunction).

$$\begin{array}{c}\hfill \left|\psi \right.〉\rangle =\left[\begin{array}{c}{e}^{{\displaystyle \frac{1}{2}}(a_1+{\mathbf{b}}_1)}\\ ⋮\\ e^{{\displaystyle \frac{1}{2}}(a_n+{\mathbf{b}}_n)}\end{array}\right]=\left[\begin{array}{c}\sqrt{{\rho }_1}{R}_1\\ ⋮\\ \sqrt{{\rho }_2}{R}_2\end{array}\right]\end{array}$$

(71)

where $\sqrt{{\rho }_i}=e^{{\displaystyle \frac{1}{2}}{a}_i}$ representing the square root of the probability and $R_i=e^{{\displaystyle \frac{1}{2}}{\mathbf{b}}_i}$ representing a rotor in 2D (or boost in 1+1D).

The partition function of the probability measure can be expressed using the bilinear form applied to the Spin(2)-valued Wavefunction:

Theorem 6 

(Partition Function). $Z=\langle 〈\psi |\psi 〉\rangle$

Proof. 

$$\begin{array}{c}\hfill \langle 〈\psi |\psi 〉\rangle =\sum _i{\psi }_i^{‡}{\psi }_i=\sum _i{\rho }_i{R}_i^{‡}{R}_i=\sum _i{\rho }_i=Z\end{array}$$

(72)

   □

Definition 8 

(Spin(2)-valued Evolution Operator).

$$\begin{array}{c}\hfill T=\left[\begin{array}{c}{e}^{-{\displaystyle \frac{1}{2}}\theta {\mathbf{b}}_1}\\ & \ddots & \\ & & e^{-{\displaystyle \frac{1}{2}}\theta {\mathbf{b}}_n}\end{array}\right]\end{array}$$

(73)

Theorem 7. 

The partition function is invariant with respect to the Spin(2)-valued evolution operator.

Proof. 

We note that:

$$\begin{array}{c}\hfill \langle 〈T\mathbf{v}|T\mathbf{v}〉\rangle =\langle 〈\mathbf{v}|\mathbf{v}〉\rangle ={\mathbf{v}}^{‡}{T}^{‡}T\mathbf{v}\Rightarrow T^{‡}T=I\end{array}$$

(74)

then, since $\left[\begin{array}{c}{e}^{{\displaystyle \frac{1}{2}}\theta {\mathbf{b}}_1}\\ & \ddots & \\ & & e^{{\displaystyle \frac{1}{2}}\theta {\mathbf{b}}_n}\end{array}\right]\left[\begin{array}{c}{e}^{-{\displaystyle \frac{1}{2}}\theta {\mathbf{b}}_1}\\ & \ddots & \\ & & e^{-{\displaystyle \frac{1}{2}}\theta {\mathbf{b}}_n}\end{array}\right]=I$, the relation $T^{‡}T=I$ is satisfied.    □

We note that the even sub-algebra of $\mathrm{GA}\left(2\right)$, being closed under addition and multiplication and constituting an inner product through its bilinear form, allows for the construction of a Hilbert space. In this context, the Hilbert space is Spin(2)-valued. The primary distinction between a wavefunction in a complex Hilbert space and one in a Spin(2)-valued Hilbert space lies in the subject matter of the theory. Specifically, in the latter, the construction governs the change in orientation experienced by an observer, which in turn dictates the measurement basis used in the experiment, consistently with the rotational symmetry and freedom of the system.

The dynamics of observer orientation transformations are described by the Schrödinger equation, which is derived by taking the derivative of the wavefunction with respect to the Lagrange multiplier, $\theta$:

Definition 9 

(Spin(2)-valued Schrödinger Equation).

$$\begin{array}{c}\hfill {\displaystyle \frac{\mathrm{d}}{\mathrm{d}\theta }}\left[\begin{array}{c}{\psi }_1\left(\theta \right)\\ ⋮\\ {\psi }_n\left(\theta \right)\end{array}\right]=\left[\begin{array}{c}-{\displaystyle \frac{1}{2}}{\mathbf{b}}_1\\ & \ddots & \\ & & -{\displaystyle \frac{1}{2}}{\mathbf{b}}_n\end{array}\right]\left[\begin{array}{c}{\psi }_1\left(\theta \right)\\ ⋮\\ {\psi }_n\left(\theta \right)\end{array}\right]\end{array}$$

(75)

Here, $\theta$ represents a global one-parameter evolution parameter akin to time, which is able to transform the wavefunction under the Spin(2), locally across the states of the Hilbert space. This is an extremely general equation that captures all transformations that can be done consistently with the symmetries of the wavefunction.

Definition 10 

(David Hestenes’ Formulation). In 3+1D, the David Hestenes’ formulation [7] of the wavefunction is $\psi =\sqrt{\rho }R{e}^{ib/2}$, where $R=e^{\mathbf{f}/2}$ is a Lorentz boost or rotation and where $e^{ib/2}$ is a phase. In 2D, as the algebra only admits a bivector, his formulation would reduce to $\psi =\sqrt{\rho }R$, which is the form we have recovered.

The definition of the Dirac current applicable to our wavefunction follows the formulation of David Hestenes:

Definition 11 

(Dirac Current). Given the basis ${\sigma }_x$ and ${\sigma }_y$, the Dirac current for the 2D theory is defined as:

$$\begin{array}{c}\hfill J_1\equiv {\psi }^{‡}{\sigma }_x\psi =\rho \underset{\mathrm{SO}\left(2\right)}{\underbrace{R^{‡}{\sigma }_{x}R}}=\rho {\tilde{\sigma }}_x\end{array}$$

(76)

$$\begin{array}{c}\hfill J_2\equiv {\psi }^{‡}{\sigma }_y\psi =\rho \underset{\mathrm{SO}\left(2\right)}{\underbrace{R^{‡}{\sigma }_{y}R}}=\rho {\tilde{\sigma }}_y\end{array}$$

(77)

where ${\tilde{\sigma }}_x$ and ${\tilde{\sigma }}_y$ are a SO(2) rotated basis vectors.

2.2.1. Obstructions

As stated, all geometric configurations except 3+1D contain obstructions. Specifically, in 1+1D and 2D, we identify two obstructions:

• In 1+1D: The 1+1D theory results in a split-complex quantum theory due to the bilinear form $(a-b{\mathbf{e}}_0\wedge {\mathbf{e}}_1)(a+b{\mathbf{e}}_0\wedge {\mathbf{e}}_1)$, which yields negative probabilities: $a^2-b^2\in \mathbb{R}$ for certain wavefunction states, in contrast to the non-negative probabilities $a^2+b^2\in {\mathbb{R}}^{\ge 0}$ obtained in the Euclidean 2D case. (This is why we had to use 2D instead of 1+1D in this two-dimensional introduction...)
• In 1+1D and in 2D: The basis vectors (${\sigma }_x$ and ${\sigma }_y$ in 2D, and ${\mathbf{e}}_0$ and ${\mathbf{e}}_1$ in 1+1D) are not self-adjoint. Although useable in the context of defining the Dirac current, their non-self-adjointness prevents the construction of the metric tensor as an observable. The benefits of having the basis vectors self-adjoint will become obvious in the 3+1D case, where we will be able to construct the metric tensor from basis measurements. Specifically, in 2D:

  $$\begin{array}{c}\hfill {\left({\mathbf{e}}_{\mu }\mathbf{u}\right)}^{‡}\mathbf{u}\ne {\mathbf{u}}^{‡}{\mathbf{e}}_{\mu }\mathbf{u}\end{array}$$

  (78)

  because ${\left({\mathbf{e}}_{\mu }\mathbf{u}\right)}^{‡}\mathbf{u}={\mathbf{u}}^{‡}{\mathbf{e}}_{\mu }^{‡}\mathbf{u}={\mathbf{u}}^{‡}(-{\mathbf{e}}_{\mu })\mathbf{u}$.

In the following section, we will explore the obstruction-free 3+1D case.

2.3. RQM in 3+1D

In this section, we extend the concepts and techniques developed for multivector amplitudes in 2D to the more physically relevant case of 3+1D dimensions. The Lagrange multiplier equation is as follows:

$$\begin{array}{c}\hfill \mathcal{L}=\underset{\mathrm{Relative}\mathrm{Shannon}\mathrm{Entropy}}{\underbrace{-{\sum }_i{\rho }_{i}ln{\displaystyle \frac{{\rho }_i}{p_i}}}}+\underset{\mathrm{Normalization}\mathrm{Constraint}}{\underbrace{\lambda \left(1-{\sum }_i{\rho }_i\right)}}+\underset{\mathrm{Vanishing}{\mathrm{Spin}}^c(3,1)\text{-}\mathrm{Phase}}{\underbrace{\zeta \left(-\mathrm{tr}{\displaystyle \frac{1}{2}}{\sum }_i{\rho }_i{\mathbf{M}}_i\right)}}\end{array}$$

(79)

The solution (proof in Appendix B) is obtained using the same step-by-step process as the 2D case, and yields:

$$\begin{array}{c}\hfill {\rho }_i=\underset{{\mathrm{Spin}}^c(3,1)\mathrm{Invariant}\mathrm{Ensemble}}{\underbrace{{\displaystyle \frac{1}{{\sum }_i{p}_{i}detexp(-\zeta \frac{1}{2}{\mathbf{M}}_i)}}}}\underset{{\mathrm{Spin}}^c(3,1)\mathrm{Born}\mathrm{Rule}}{\underbrace{detexp(-\zeta {\displaystyle \frac{1}{2}}{\mathbf{M}}_i)}}\underset{\mathrm{Initial}\mathrm{Preparation}}{\underbrace{p_i}}\end{array}$$

(80)

where $\zeta$ is a "twisted-phase" rapidity. (If the invariance group was Spin(3,1) instead of Spin^c(3,1), obtainable by posing $\mathbf{b}\to 0$, then it would simply be the rapidity).

2.3.1. Preliminaries

Our initial goal will be to express the partition function as a self-product of elements of the vector space. As such, we begin by defining a general multivector in the geometric algebra $\mathrm{GA}(3,1)$.

Definition 12 

(Multivector). Let $\mathbf{u}$ be a multivector of $\mathrm{GA}(3,1)$. Its general form is:

$$\begin{array}{cc}\hfill \mathbf{u}& =a\hfill \end{array}$$

(81)

$$\begin{array}{cc}\hfill & +t{\gamma }_0+x{\gamma }_1+y{\gamma }_2+z{\gamma }_3\hfill \end{array}$$

(82)

$$\begin{array}{cc}\hfill & +f_{01}{\gamma }_0\wedge {\gamma }_1+f_{02}{\gamma }_0\wedge {\gamma }_2+f_{03}{\gamma }_0\wedge {\gamma }_3+f_{12}{\gamma }_1\wedge {\gamma }_2+f_{13}{\gamma }_1\wedge {\gamma }_3+f_{23}{\gamma }_2\wedge {\gamma }_3\hfill \end{array}$$

(83)

$$\begin{array}{cc}\hfill & +p{\gamma }_1\wedge {\gamma }_2\wedge {\gamma }_3+q{\gamma }_0\wedge {\gamma }_2\wedge {\gamma }_3+v{\gamma }_0\wedge {\gamma }_1\wedge {\gamma }_3+w{\gamma }_0\wedge {\gamma }_1\wedge {\gamma }_2\hfill \end{array}$$

(84)

$$\begin{array}{cc}\hfill & +b{\gamma }_0\wedge {\gamma }_1\wedge {\gamma }_2\wedge {\gamma }_3\hfill \end{array}$$

(85)

where ${\gamma }_0,{\gamma }_1,{\gamma }_2,{\gamma }_3$ are the basis vectors in the real Majorana representation.

A more compact notation for $\mathbf{u}$ is

$$\begin{array}{c}\hfill \mathbf{u}=a+\mathbf{x}+\mathbf{f}+\mathbf{v}+\mathbf{b}\end{array}$$

(86)

where a is a scalar, $\mathbf{x}$ a vector, $\mathbf{f}$ a bivector, $\mathbf{v}$ is pseudo-vector and $\mathbf{b}$ a pseudo-scalar.

This general multivector can be represented by a $4\times 4$ real matrix using the real Majorana representation:

Definition 13 

(Matrix Representation of $\mathbf{u}$).

$$\begin{array}{c}\hfill \mathbf{M}=\left[\begin{array}{cccc}a+f_{02}-q-z& b-f_{13}+w-x& -f_{01}+f_{12}-p+v& f_{03}+f_{23}+t+y\\ -b+f_{13}+w-x& a+f_{02}+q+z& f_{03}+f_{23}-t-y& f_{01}-f_{12}-p+v\\ -f_{01}-f_{12}+p+v& f_{03}-f_{23}+t-y& a-f_{02}+q-z& -b-f_{13}-w-x\\ f_{03}-f_{23}-t+y& f_{01}+f_{12}+p+v& b+f_{13}-w-x& a-f_{02}-q+z\end{array}\right]\end{array}$$

(87)

To manipulate and analyze multivectors in $\mathrm{GA}(3,1)$, we introduce several important operations, such as the multivector conjugate, the 3,4 blade conjugate, and the multivector self-product.

Definition 14 

(Multivector Conjugate (in 4D)).

$$\begin{array}{c}\hfill {\mathbf{u}}^{‡}=a-\mathbf{x}-\mathbf{f}+\mathbf{v}+\mathbf{b}\end{array}$$

(88)

Definition 15 

(3,4 Blade Conjugate). The 3,4 blade conjugate of $\mathbf{u}$ is

$$\begin{array}{c}\hfill {\lfloor \mathbf{u}\rfloor }_{3,4}=a+\mathbf{x}+\mathbf{f}-\mathbf{v}-\mathbf{b}\end{array}$$

(89)

The results of Lundholm[8], demonstrates that the multivector norms in the following definition, are the unique forms which carries the properties of the determinants such as $N\left(\mathbf{u}\mathbf{v}\right)=N\left(\mathbf{u}\right)N\left(\mathbf{v}\right)$ to the domain of multivectors:

Definition 16. 

The self-products associated with low-dimensional geometric algebras are:

$$\begin{array}{cccc}\hfill & \mathrm{GA}(0,1):\hfill & \hfill & {\phi }^{†}\phi \hfill \end{array}$$

(90)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(2,0):\hfill & \hfill & {\phi }^{‡}\phi \hfill \end{array}$$

(91)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(3,0):\hfill & \hfill & {\lfloor {\phi }^{‡}\phi \rfloor }_3{\phi }^{‡}\phi \hfill \end{array}$$

(92)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(3,1):\hfill & \hfill & {\lfloor {\phi }^{‡}\phi \rfloor }_{3,4}{\phi }^{‡}\phi \hfill \end{array}$$

(93)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(4,1):\hfill & \hfill & {\left({\lfloor {\phi }^{‡}\phi \rfloor }_{3,4}{\phi }^{‡}\phi \right)}^{†}\left({\lfloor {\phi }^{‡}\phi \rfloor }_{3,4}{\phi }^{‡}\phi \right)\hfill \end{array}$$

(94)

We can now express the determinant of the matrix representation of a multivector via the self-product ${\lfloor {\phi }^{‡}\phi \rfloor }_{3,4}{\phi }^{‡}\phi$. Again, this choice is not arbitrary, but the unique choice with allows us to represent the determinant of the matrix representation of a multivector within $\mathrm{GA}(3,1)$:

Theorem 8 

(Determinant as a Multivector Self-Product).

$$\begin{array}{c}\hfill {\lfloor {\mathbf{u}}^{‡}\mathbf{u}\rfloor }_{3,4}{\mathbf{u}}^{‡}\mathbf{u}=det\mathbf{M}\end{array}$$

(95)

Proof. 

Please find a computer assisted proof of this equality in Appendix C.    □

Definition 17 

($\mathrm{GA}(3,1)$-valued Vector).

$$\begin{array}{c}\hfill \left|V\right.〉\rangle =\left[\begin{array}{c}{\mathbf{u}}_1\\ ⋮\\ {\mathbf{u}}_n\end{array}\right]=\left[\begin{array}{c}{a}_1+{\mathbf{x}}_1+{\mathbf{f}}_1+{\mathbf{v}}_1+{\mathbf{b}}_1\\ ⋮\\ a_n+{\mathbf{x}}_n+{\mathbf{f}}_n+{\mathbf{v}}_n+{\mathbf{b}}_n\end{array}\right]\end{array}$$

(96)

These constructions allow us to express the partition function in terms of the multivector self-product:

Definition 18 

(Double-Copy Product). Instead of an inner product, we obtain what we call a double-copy product:

$$\begin{array}{cc}\hfill \langle 〈V\left|V\right|V|V〉\rangle & =\sum _i{\lfloor \underset{\mathrm{copy}1}{\underbrace{{\psi }_i^{‡}{\psi }_i}}\rfloor }_{3,4}\underset{\mathrm{copy}2}{\underbrace{{\psi }_i^{‡}{\psi }_i}}\hfill \end{array}$$

(97)

$$\begin{array}{cc}\hfill & ={\lfloor \underset{\mathrm{copy}1}{\underbrace{\left[\begin{array}{ccc}{\mathbf{u}}_1^{‡}& \cdots & {\mathbf{u}}_n\end{array}\right]\left[\begin{array}{ccc}{\mathbf{u}}_1& \cdots & 0\\ ⋮& \ddots & ⋮\\ 0& \cdots & {\mathbf{u}}_n\end{array}\right]}}\rfloor }_{3,4}\underset{\mathrm{copy}2}{\underbrace{\left[\begin{array}{ccc}{\mathbf{u}}_1^{‡}& \cdots & 0\\ ⋮& \ddots & ⋮\\ 0& \cdots & {\mathbf{u}}_n^{‡}\end{array}\right]\left[\begin{array}{c}{\mathbf{u}}_1\\ ⋮\\ {\mathbf{u}}_n\end{array}\right]}}\hfill \end{array}$$

(98)

Theorem 9 

(Partition Function). $Z=\langle 〈V\left|V\right|V|V〉\rangle$

Proof. 

$$\begin{array}{cc}\hfill & \langle 〈V\left|V\right|V|V〉\rangle \hfill \end{array}$$

(99)

$$\begin{array}{cc}\hfill & ={\lfloor \left[\begin{array}{ccc}{\mathbf{u}}_1^{‡}& \cdots & {\mathbf{u}}_n\end{array}\right]\left[\begin{array}{ccc}{\mathbf{u}}_1& \cdots & 0\\ ⋮& \ddots & ⋮\\ 0& \cdots & {\mathbf{u}}_n\end{array}\right]\rfloor }_{3,4}\left[\begin{array}{ccc}{\mathbf{u}}_1^{‡}& \cdots & 0\\ ⋮& \ddots & ⋮\\ 0& \cdots & {\mathbf{u}}_n^{‡}\end{array}\right]\left[\begin{array}{c}{\mathbf{u}}_1\\ ⋮\\ {\mathbf{u}}_n\end{array}\right]\hfill \end{array}$$

(100)

$$\begin{array}{cc}\hfill & ={\lfloor \left[\begin{array}{ccc}{\mathbf{u}}_1^{‡}{\mathbf{u}}_1& \cdots & {\mathbf{u}}_n{\mathbf{u}}_n\end{array}\right]\rfloor }_{3,4}\left[\begin{array}{c}{\mathbf{u}}_1^{‡}{\mathbf{u}}_1\\ ⋮\\ {\mathbf{u}}_n^{‡}{\mathbf{u}}_n\end{array}\right]\hfill \end{array}$$

(101)

$$\begin{array}{cc}\hfill & ={\lfloor {\mathbf{u}}_1^{‡}{\mathbf{u}}_1\rfloor }_{3,4}{\mathbf{u}}_1^{‡}{\mathbf{u}}_1+\cdots +{\lfloor {\mathbf{u}}_n^{‡}{\mathbf{u}}_n\rfloor }_{3,4}{\mathbf{u}}_n^{‡}{\mathbf{u}}_n\hfill \end{array}$$

(102)

$$\begin{array}{cc}\hfill & =\sum _{i=1}^{n}det{\mathbf{M}}_{{\mathbf{u}}_i}\hfill \end{array}$$

(103)

$$\begin{array}{cc}\hfill & =Z\hfill \end{array}$$

(104)

   □

Desirable properties for the double-copy product are introduced by reducing multivectors to its subgroups. First, non-negativity:

Theorem 10 

(Non-negativity). The double-copy product, applied to the even sub-algebra of $\mathrm{GA}(3,1)$ is always non-negative.

Proof. 

Let $\left|V\right.〉\rangle =\left[\begin{array}{c}{a}_1+{\mathbf{f}}_1+{\mathbf{b}}_1\\ ⋮\\ a_n+{\mathbf{f}}_n+{\mathbf{b}}_n\end{array}\right]$. Then,

$$\begin{array}{cc}\hfill & \langle 〈V\left|V\right|V|V〉\rangle \hfill \end{array}$$

(105)

$$\begin{array}{cc}\hfill & ={\lfloor \left[\begin{array}{cc}{(a_1+{\mathbf{f}}_1+{\mathbf{b}}_1)}^{‡}(a_1+{\mathbf{f}}_1+{\mathbf{b}}_1)& \cdots \end{array}\right]\rfloor }_{3,4}\left[\begin{array}{c}{(a_1+{\mathbf{f}}_1+{\mathbf{b}}_1)}^{‡}(a_1+{\mathbf{f}}_1+{\mathbf{b}}_1)\\ ⋮\end{array}\right]\hfill \end{array}$$

(106)

$$\begin{array}{cc}\hfill & ={\lfloor \left[\begin{array}{cc}(a_1-{\mathbf{f}}_1+{\mathbf{b}}_1)(a_1+{\mathbf{f}}_1+{\mathbf{b}}_1)& \cdots \end{array}\right]\rfloor }_{3,4}\left[\begin{array}{c}(a_1-{\mathbf{f}}_1+{\mathbf{b}}_1)(a_1+{\mathbf{f}}_1+{\mathbf{b}}_1)\\ ⋮\end{array}\right]\hfill \end{array}$$

(107)

$$\begin{array}{cc}\hfill & ={\lfloor \left[\begin{array}{cc}{a}_1^2+a_1{\mathbf{f}}_1+a_1{\mathbf{b}}_1-{\mathbf{f}}_1{a}_1-{\mathbf{f}}_1^2-{\mathbf{f}}_1{\mathbf{b}}_1+{\mathbf{b}}_1{a}_1+{\mathbf{b}}_1{\mathbf{f}}_1+{\mathbf{b}}_1^2& \cdots \end{array}\right]\rfloor }_{3,4}\cdots \hfill \end{array}$$

(108)

$$\begin{array}{cc}\hfill & ={\lfloor \left[\begin{array}{cc}{a}_1^2-{\mathbf{f}}_1^2+{\mathbf{b}}_1^2& \cdots \end{array}\right]\rfloor }_{3,4}\cdots \hfill \end{array}$$

(109)

We note 1) ${\mathbf{b}}^2={\left(bI\right)}^2=-b^2$ and 2) ${\mathbf{f}}^2=-E_1^2-E_2^2-E_3^2+B_1^2+B_2^2+B_3^2+4e_0{e}_1{e}_2{e}_3(E_1{B}_1+E_2{B}_2+E_3{B}_3)$

$$\begin{array}{c}={\lfloor \left[\begin{array}{cc}{a}_1^2-b_1^2+E_1^2+E_2^2+E_3^2-B_1^2-B_2^2-B_3^2-4e_0{e}_1{e}_2{e}_3(E_1{B}_1+E_2{B}_2+E_3{B}_3)& \cdots \end{array}\right]\rfloor }_{3,4}\cdots \hfill \end{array}$$

(110)

We note that the terms are now complex numbers, which we rewrite as $\Re \left(z\right)=a_1^2-b_1^2+E_1^2+E_2^2+E_3^2-B_1^2-B_2^2-B_3^2$ and $\Im \left(z\right)=-4(E_1{B}_1+E_2{B}_2+E_3{B}_3)$

$$\begin{array}{cc}\hfill 1-1& ={\lfloor \left[\begin{array}{ccc}{z}_1& \cdots & z_2\end{array}\right]\rfloor }_{3,4}\left[\begin{array}{c}{z}_n\\ ⋮\\ z_n\end{array}\right]\hfill \end{array}$$

(111)

$$\begin{array}{cc}\hfill & =\left[\begin{array}{ccc}{z}_1^{†}& \cdots & z_2^{†}\end{array}\right]\left[\begin{array}{c}{z}_n\\ ⋮\\ z_n\end{array}\right]\hfill \end{array}$$

(112)

$$\begin{array}{cc}\hfill & =z_1^{‡}{z}_1+\cdots +z_n^{‡}{z}_n\hfill \end{array}$$

(113)

which is always non-negative.    □

Then, positive-definiteness of the double-copy product is obtained by creating an equivalence class between the zero vector and any non-zero vector of length zero, and taking the zero vector as the representative of the class. To realize the equivalence class, we define the ${\mathrm{Spin}}^{\mathrm{c}}(3,1)$-valued wavefunction, which is valued in the even sub-algebra of $\mathrm{GA}(3,1)$, as follows:

Definition 19 

(${\mathrm{Spin}}^{\mathrm{c}}(3,1)$-valued Wavefunction).

$$\begin{array}{c}\hfill \left|\psi \right.〉\rangle =\left[\begin{array}{c}{e}^{{\displaystyle \frac{1}{2}}(a_1+{\mathbf{f}}_1+{\mathbf{b}}_1)}\\ ⋮\\ e^{{\displaystyle \frac{1}{2}}(a_n+{\mathbf{f}}_n+{\mathbf{b}}_n)}\end{array}\right]=\left[\begin{array}{c}\sqrt{{\rho }_1}{R}_1{B}_1\\ ⋮\\ \sqrt{{\rho }_n}{R}_n{B}_n\end{array}\right]\end{array}$$

(114)

where $R_i=e^{{\displaystyle \frac{1}{2}}{\mathbf{f}}_i}$ is a rotor, $B_i=e^{{\displaystyle \frac{1}{2}}{\mathbf{b}}_i}$ is a phase, and where $\sqrt{{\rho }_i}=e^{{\displaystyle \frac{1}{2}}{a}_i}\ge 0$.

Any even multivectors of GA(3,1) admits a unique exponential representation, except when ${\rho }_i=0$ in which it is surjective. Consequently, in this representation the double-copy product yields 0 only for the zero vector, rendering the double-copy product positive-definite.

Now, let us turn our attention to the evolution operator, which leaves the partition function invariant:

Definition 20 

(${\mathrm{Spin}}^{\mathrm{c}}(3,1)$ Evolution Operator).

$$\begin{array}{c}\hfill T=\left[\begin{array}{c}{e}^{-{\displaystyle \frac{1}{2}}\zeta ({\mathbf{f}}_1+{\mathbf{b}}_1)}\\ & \ddots \\ & & e^{-{\displaystyle \frac{1}{2}}\zeta ({\mathbf{f}}_n+{\mathbf{b}}_n)}\end{array}\right]\end{array}$$

(115)

In turn, this leads to a Schrödinger equation obtained by taking the derivative of the wavefunction with respect to the Lagrange multiplier $\zeta$:

Definition 21 

(${\mathrm{Spin}}^{\mathrm{c}}(3,1)$-valued Schrödinger equation).

$$\begin{array}{c}\hfill {\displaystyle \frac{\mathrm{d}}{\mathrm{d}\zeta }}\left[\begin{array}{c}{\psi }_1\left(\zeta \right)\\ ⋮\\ {\psi }_n\left(\zeta \right)\end{array}\right]=\left[\begin{array}{c}-{\displaystyle \frac{1}{2}}({\mathbf{f}}_1+{\mathbf{b}}_1)\\ & & \ddots \\ & & & -{\displaystyle \frac{1}{2}}({\mathbf{f}}_n+{\mathbf{b}}_n)\end{array}\right]\left[\begin{array}{c}{\psi }_1\left(\zeta \right)\\ ⋮\\ {\psi }_n\left(\zeta \right)\end{array}\right]\end{array}$$

(116)

In this case $\zeta$ represents a one-parameter evolution parameter akin to time, which is able to transform the measurement basis under action of the ${\mathrm{Spin}}^{\mathrm{c}}(3,1)$ group. This is an extremely general equation that captures all transformations that can be done consistently with the symmetries of the wavefunction.

Theorem 11 

(Spin^c(3,1) invariance). Let $e^{{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}$ be a general element of Spin^c(3,1). Then, the equality:

$$\begin{array}{c}\hfill {\lfloor {\psi }^{‡}\psi \rfloor }_{3,4}{\psi }^{‡}\psi ={\lfloor {\left(e^{{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \right)}^{‡}{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \rfloor }_{3,4}{\left(e^{{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \right)}^{‡}{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \end{array}$$

(117)

is always satisfied.

Proof. 

$$\begin{array}{cc}\hfill & {\lfloor {\left(e^{{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \right)}^{‡}{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \rfloor }_{3,4}{\left(e^{{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \right)}^{‡}{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \hfill \end{array}$$

(118)

$$\begin{array}{cc}\hfill & ={\lfloor {\psi }^{‡}{e}^{-{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \rfloor }_{3,4}{\psi }^{‡}{e}^{-{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \hfill \end{array}$$

(119)

$$\begin{array}{cc}\hfill & ={\lfloor {\psi }^{‡}{e}^{\mathbf{b}}\psi \rfloor }_{3,4}{\psi }^{‡}{e}^{\mathbf{b}}\psi \hfill \end{array}$$

(120)

$$\begin{array}{cc}\hfill & ={\lfloor {\psi }^{‡}\psi \rfloor }_{3,4}{e}^{-\mathbf{b}}{e}^{\mathbf{b}}{\psi }^{‡}\psi \hfill \end{array}$$

(121)

$$\begin{array}{cc}\hfill & ={\lfloor {\psi }^{‡}\psi \rfloor }_{3,4}{\psi }^{‡}\psi \hfill \end{array}$$

(122)

   □

2.3.2. RQM

Definition 22 

(David Hestenes’ Wavefunction). The ${\mathrm{Spin}}^{\mathrm{c}}(3,1)$-valued wavefunction we have recovered is formulated identically to David Hestenes’[7] formulation of the wavefunction within GA(3,1).

$$\begin{array}{c}\hfill \psi =e^{{\displaystyle \frac{1}{2}}(a+\mathbf{f}+\mathbf{b})}=\sqrt{\rho }R{e}^{-ib/2}\end{array}$$

(123)

where $e^{{\displaystyle \frac{1}{2}}a}=\sqrt{\rho }$, $e^{{\displaystyle \frac{1}{2}}\mathbf{f}}=R$ and $e^{{\displaystyle \frac{1}{2}}\mathbf{b}}=e^{-ib/2}$.

Before we continue the RQM investigation, let us note that the double-copy product contains two copies of a bilinear form ${\psi }^{‡}\psi$:

$$\begin{array}{c}\hfill {\lfloor \underset{\mathrm{copy}1}{\underbrace{{\psi }^{‡}\psi }}\rfloor }_{3,4}\underset{\mathrm{copy}2}{\underbrace{{\psi }^{‡}\psi }}\end{array}$$

(124)

In the present section, we will investigate the properties of each copy individually, leaving the properties specific to the double-copy for the section on quantum gravity.

Taking a single copy, the Dirac current is obtained directly from the gamma matrices, as follows:

Definition 23 

(Dirac Current). The definition of the Dirac current is the same as Hestenes’:

$$\begin{array}{cc}\hfill J& \equiv {\psi }^{‡}{\gamma }_{\mu }\psi =\rho R^{‡}{B}^{‡}{\gamma }_{\mu }BR=\rho R^{‡}{\gamma }_{\mu }{B}^{-1}BR=\rho \underset{\mathrm{SO}(3,1)}{\underbrace{R^{‡}{\gamma }_{\mu }R}}=\rho {\tilde{\gamma }}_{\mu }\hfill \end{array}$$

(125)

where ${\tilde{\gamma }}_{\mu }$ is a SO(3,1) rotated basis vector.

2.3.3. Standard Model Gauge Symmetries

We will now demonstrate that the copied bilinear form is automatically invariant with respect to the $\mathrm{U}\left(1\right)$, $\mathrm{SU}\left(2\right)$, and $\mathrm{SU}\left(3\right)$ symmetries and the unitary $U^{†}U=I$ symmetry which play a fundamental role in the standard model of particle physics. These constitute the set of symmetries that stabilize the Dirac current ${\left(T\psi \right)}^{‡}{\gamma }_{0}T\psi ={\psi }^{‡}{\gamma }_0\psi$.

Theorem 12 

(U(1) Invariance). Let $e^{{\displaystyle \frac{1}{2}}\mathbf{b}}$ be a general element of U(1). Then, the equality

$$\begin{array}{c}\hfill {\lfloor {\psi }^{‡}{\gamma }_0\psi \rfloor }_{3,4}{\psi }^{‡}{\gamma }_0\psi ={\lfloor \underset{\mathrm{copy}1}{\underbrace{{\left(e^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \right)}^{‡}{\gamma }_0{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi }}\rfloor }_{3,4}\underset{\mathrm{copy}2}{\underbrace{{\left(e^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi \right)}^{‡}{\gamma }_0{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}\psi }}\end{array}$$

(126)

is satisfied, yielding a U(1) symmetry for each copied bilinear form.

Proof. 

Equation 126 is invariant if this expression is satisfied:

$$\begin{array}{c}\hfill e^{{\displaystyle \frac{1}{2}}\mathbf{b}}{\gamma }_0{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}={\gamma }_0\end{array}$$

(127)

This is always satisfied simply because $e^{{\displaystyle \frac{1}{2}}\mathbf{b}}{\gamma }_0{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}={\gamma }_0{e}^{-{\displaystyle \frac{1}{2}}\mathbf{b}}{e}^{{\displaystyle \frac{1}{2}}\mathbf{b}}={\gamma }_0$    □

Theorem 13 

(SU(2) Invariance). Let $e^{{\displaystyle \frac{1}{2}}\mathbf{f}}$ be a general element of Spin(3,1). Then, the equality:

$$\begin{array}{c}\hfill {\lfloor {\psi }^{‡}{\gamma }_0\psi \rfloor }_{3,4}{\psi }^{‡}{\gamma }_0\psi ={\lfloor \underset{\mathrm{copy}1}{\underbrace{{\left(e^{{\displaystyle \frac{1}{2}}\mathbf{f}}\psi \right)}^{‡}{\gamma }_0{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}\psi }}\rfloor }_{3,4}\underset{\mathrm{copy}2}{\underbrace{{\left(e^{{\displaystyle \frac{1}{2}}\mathbf{f}}\psi \right)}^{‡}{\gamma }_0{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}\psi }}\end{array}$$

(128)

is satisfied for if $\mathbf{f}={\theta }_1{\gamma }_2{\gamma }_3+{\theta }_2{\gamma }_1{\gamma }_3+{\theta }_3{\gamma }_1{\gamma }_2$ (which generates SU(2)), yielding a SU(2) symmetry for each copied bilinear form.

Proof. 

Equation 128 is invariant if this expression is satisfied[9]:

$$\begin{array}{c}\hfill e^{-{\displaystyle \frac{1}{2}}\mathbf{f}}{\gamma }_0{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}={\gamma }_0\end{array}$$

(129)

We now note that moving the left-most term to the right of the gamma matrix yields:

$$\begin{array}{cc}\hfill & e^{-E_1{\gamma }_0{\gamma }_1-E_2{\gamma }_0{\gamma }_2-E_3{\gamma }_0{\gamma }_3-{\theta }_1{\gamma }_2{\gamma }_3-{\theta }_2{\gamma }_1{\gamma }_3-{\theta }_3{\gamma }_1{\gamma }_2}{\gamma }_0{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}\hfill \end{array}$$

(130)

$$\begin{array}{cc}\hfill & ={\gamma }_0{e}^{E_1{\gamma }_0{\gamma }_1+E_2{\gamma }_0{\gamma }_2+E_3{\gamma }_0{\gamma }_3-{\theta }_1{\gamma }_2{\gamma }_3-{\theta }_2{\gamma }_1{\gamma }_3-{\theta }_3{\gamma }_1{\gamma }_2}{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}\hfill \end{array}$$

(131)

Therefore, the product $e^{-{\displaystyle \frac{1}{2}}\mathbf{f}}{\gamma }_0{e}^{{\displaystyle \frac{1}{2}}\mathbf{f}}$ reduces to ${\gamma }_0$ if and only if $E_1=E_2=E_3=0$, leaving $\mathbf{f}={\theta }_1{\gamma }_2{\gamma }_3+{\theta }_2{\gamma }_1{\gamma }_3+{\theta }_3{\gamma }_1{\gamma }_2$:

Finally, we note that $e^{{\theta }_1{\gamma }_2{\gamma }_3+{\theta }_2{\gamma }_1{\gamma }_3+{\theta }_3{\gamma }_1{\gamma }_2}$ generates $\mathrm{SU}\left(2\right)$.    □

Theorem 14 

(SU(3)). The generators of SU(3) in GA(3,1) are given by Anthony Lesenby in [10] and are as follows:

$$\begin{array}{cccc}\hfill {\widehat{E}}_{ij}& ={\widehat{e}}_i{\widehat{e}}_j-{\widehat{f}}_i{\widehat{f}}_j\hfill & \hfill & \mathrm{where}i<j\hfill \end{array}$$

(132)

$$\begin{array}{cccc}\hfill {\widehat{F}}_{ij}& ={\widehat{e}}_i{\widehat{f}}_j+{\widehat{e}}_j{\widehat{f}}_i\hfill & \hfill & \mathrm{where}i<j\hfill \end{array}$$

(133)

$$\begin{array}{cccc}\hfill \widehat{J}& ={\widehat{e}}_i{\widehat{f}}_i\hfill & \hfill & \mathrm{where}i=1,2,3\hfill \end{array}$$

(134)

where

$$\begin{array}{cc}\hfill {\widehat{e}}_i& =\mathrm{multiplication}\mathrm{on}\mathrm{the}\mathrm{left}\mathrm{by}{\sigma }_i,\mathrm{so}\mathrm{that}{\widehat{e}}_i\left(F\right)={\sigma }_{i}F\hfill \end{array}$$

(135)

$$\begin{array}{cc}\hfill {\widehat{f}}_i& =\mathrm{multiplication}\mathrm{on}\mathrm{the}\mathrm{right}\mathrm{by}I{\sigma }_i,\mathrm{so}\mathrm{that}{\widehat{f}}_i\left(F\right)=I{\sigma }_{i}F\hfill \end{array}$$

(136)

This defines the 9 generators of U(3).

With the additional restriction on $\widehat{J}$

$$\begin{array}{c}\hfill {\alpha }_1{\widehat{J}}_1+{\alpha }_2{\widehat{J}}_2+{\alpha }_3{\widehat{J}}_3,\mathrm{with}{\alpha }_1+{\alpha }_2+{\alpha }_3=0\end{array}$$

(137)

the number generators is reduced to 8, consistently with SU(3).

We now must show that the following equation is satisfied for all 8 generators:

$$\begin{array}{c}\hfill {\lfloor {\psi }^{‡}{\gamma }_0\psi \rfloor }_{3,4}{\psi }^{‡}{\gamma }_0\psi ={\lfloor \underset{\mathrm{copy}1}{\underbrace{{\left(e^{{\theta }^i{\lambda }_i}\psi \right)}^{‡}{\gamma }_0{e}^{{\theta }^i{\lambda }_i}\psi }}\rfloor }_{3,4}\underset{\mathrm{copy}2}{\underbrace{{\left(e^{{\theta }^i{\lambda }_i}\psi \right)}^{‡}{\gamma }_0{e}^{{\theta }^i{\lambda }_i}\psi }}\end{array}$$

(138)

Proof. 

First, we note the following action:

$$\begin{array}{c}\hfill -\mathbf{f}{\gamma }_0\mathbf{f}={\gamma }_0\end{array}$$

(139)

which we can rewrite as follows:

$$\begin{array}{c}\hfill -(E_1{\gamma }_0{\gamma }_1+E_2{\gamma }_0{\gamma }_2+E_3{\gamma }_0{\gamma }_3+B_1{\gamma }_2{\gamma }_3+B_2{\gamma }_1{\gamma }_3+B_3{\gamma }_1{\gamma }_2){\gamma }_0\mathbf{f}\end{array}$$

(140)

The first three terms anticommute with ${\gamma }_0$, while the last three commute with ${\gamma }_0$:

$$\begin{array}{c}\hfill ={\gamma }_0(E_1{\gamma }_0{\gamma }_1+E_2{\gamma }_0{\gamma }_2+E_3{\gamma }_0{\gamma }_3-B_1{\gamma }_2{\gamma }_3-B_2\left(q\right){\gamma }_1{\gamma }_3-B_3\left(q\right){\gamma }_1{\gamma }_2)\mathbf{f}\left(q\right)\end{array}$$

(141)

This can be written as:

$$\begin{array}{cc}\hfill & {\gamma }_0(\mathbf{E}-\mathbf{B})(\mathbf{E}+\mathbf{B})\hfill \end{array}$$

(142)

$$\begin{array}{cc}\hfill & ={\gamma }_0({\mathbf{E}}^2+\mathbf{E}\mathbf{B}-\mathbf{B}\mathbf{E}-{\mathbf{B}}^2)\hfill \end{array}$$

(143)

where $\mathbf{E}=E_1{\gamma }_0{\gamma }_1+E_2{\gamma }_0{\gamma }_2+E_3{\gamma }_0{\gamma }_3$ and $\mathbf{B}=B_1{\gamma }_2{\gamma }_3+B_2{\gamma }_1{\gamma }_3+B_3{\gamma }_1{\gamma }_2$.

Thus, for $-\mathbf{f}{\gamma }_0\mathbf{f}={\gamma }_0$, we require: 1) ${\mathbf{E}}^2-{\mathbf{B}}^2=1$ and 2) $\mathbf{E}\mathbf{B}=\mathbf{B}\mathbf{E}$. The first requirement expands as follows:

$$\begin{array}{c}\hfill {\mathbf{E}}^2-{\mathbf{B}}^2=(E_1^2+B_1^2)+(E_2^2+B_2^2)+(E_3^2+B_3^2)=1\end{array}$$

(144)

which is the defining conditions for the $\mathrm{SU}\left(3\right)$ symmetry group.

Finally, as the SU(3) norm is a consequence of preserving the Dirac current, it follows that the SU(3) generators provided by Lasenby, acting on $\mathbf{f}$, cannot change the SU(3) norm, hence must also preserve the Dirac current.    □

Theorem 15 

(Unitary invariance). Let U be $n\times n$ unitary matrices. Then unitary invariance:

$$\begin{array}{c}\hfill \langle 〈\psi |{\gamma }_{\mu }\psi \left|\psi \right|{\gamma }_{\nu }\psi 〉\rangle =\langle 〈U\psi |{\gamma }_{\mu }U\psi \left|U\psi \right|{\gamma }_{\nu }U\psi 〉\rangle \Rightarrow U^{†}U=I\end{array}$$

(145)

is individually satisfied for each copied bilinear form.

Proof. 

Equation 145 is satisfied if $U^{‡}{\gamma }_{\mu }U={\gamma }_{\mu }$. Since U is valued in complex numbers, then $U^{‡}=U^T$, and since ${\gamma }_{\mu }{\gamma }_0{\gamma }_1{\gamma }_2{\gamma }_3=-{\gamma }_0{\gamma }_1{\gamma }_2{\gamma }_3{\gamma }_{\mu }$, it follows that:

$$\begin{array}{c}\hfill {\gamma }_{\mu }{U}^{†}U={\gamma }_{\mu }\end{array}$$

(146)

which is satisfied when $U^{†}U=I$.    □

The invariances SU(3), SU(2) and U(1) discussed above can be promoted to local symmetries using the usual gauge symmetry construction techniques, along with the Dirac equation or field Lagrangian.

In conventional QM, the Born rule naturally leads to a U(1)-valued gauge theory due to the following symmetry:

$$\begin{array}{c}\hfill {\left(e^{-i\theta \left(x\right)}\psi \left(x\right)\right)}^{†}{e}^{-i\theta \left(x\right)}\psi \left(x\right)=\psi {\left(x\right)}^{†}\psi \left(x\right)\end{array}$$

(147)

However, the SU(3) and SU(2) symmetries do not emerge from the probability measure in the same way and must instead be introduced manually, justified by experimental considerations. This raises the question: why these specific symmetries and not others? In contrast, within the double-copy product framework, all three symmetry groups–U(1), SU(2), and SU(3)–as well as the Spin(3,1) and unitary symmetries, follow naturally from the invariance of the probability measure, in the same way that U(1) symmetry follows from the Born rule.

2.3.4. Quantum Gravity

In the previous section, we developed a quantum theory valued in Spin^c(3,1), which served as the arena for RQM. We then demonstrated how a single copy of this theory leads to the gauge symmetries of the standard model. The goal of this section is to extend this methodology to arbitrary basis vectors, in which the metric tensor emerges as an observable. To achieve this, we will utilize both copies.

Our formulation is reminiscent of the Bern-Carrasco-Johansson (BCJ) double-copy approach to perturbatively expanded quantum gravity [11]. However, our double-copy is applied directly at the level of the Dirac current, rather than to gauge theory amplitudes.

By applying the double-copy product to the Dirac current, we establish a connection between quantum theory and the geometrical structure of spacetime.

We recall the definition of the metric tensor in terms of basis vectors of geometric algebra, as follows:

$$\begin{array}{c}\hfill g_{\mu \nu }={\displaystyle \frac{1}{2}}({\mathbf{e}}_{\mu }{\mathbf{e}}_{\nu }+{\mathbf{e}}_{\nu }{\mathbf{e}}_{\mu })\end{array}$$

(148)

Then, we note that the double-copy product acts on a pair of basis element ${\mathbf{e}}_{\mu }$ and ${\mathbf{e}}_{\nu }$, as follows:

$$\begin{array}{cc}\hfill & {\displaystyle \frac{1}{2}}\left({\lfloor \underset{\mathrm{copy}1}{\underbrace{{\psi }^{‡}{\mathbf{e}}_{\mu }\psi }}\rfloor }_{3,4}\underset{\mathrm{copy}2}{\underbrace{{\psi }^{‡}{\mathbf{e}}_{\nu }\psi }}+{\lfloor \underset{\mathrm{copy}2}{\underbrace{{\psi }^{‡}{\mathbf{e}}_{\nu }\phi }}\rfloor }_{3,4}\underset{\mathrm{copy}1}{\underbrace{{\psi }^{‡}{\mathbf{e}}_{\mu }\psi }}\right)\hfill \end{array}$$

(149)

$$\begin{array}{cc}\hfill & ={\displaystyle \frac{1}{2}}\left(\tilde{R}\underset{\mathrm{Born}\mathrm{rule}\mathrm{copy}1}{\underbrace{\rho e^{ib/2}{e}^{-ib/2}}}{\mathbf{e}}_{\mu }R\tilde{R}\underset{\mathrm{Born}\mathrm{rule}\mathrm{copy}2}{\underbrace{\rho e^{-ib/2}{e}^{ib/2}}}{\mathbf{e}}_{\nu }R+\tilde{R}\underset{\mathrm{Born}\mathrm{rule}\mathrm{copy}2}{\underbrace{\rho e^{ib/2}{e}^{-ib/2}}}{\mathbf{e}}_{\nu }R\tilde{R}\underset{\mathrm{Born}\mathrm{rule}\mathrm{copy}1}{\underbrace{\rho e^{-ib/2}{e}^{ib/2}}}{\mathbf{e}}_{\mu }R\right)\hfill \end{array}$$

(150)

$$\begin{array}{cc}\hfill & ={\displaystyle \frac{1}{2}}{\rho }^2\left(\tilde{R}{\mathbf{e}}_{\mu }R\tilde{R}{\mathbf{e}}_{\nu }R+\tilde{R}{\mathbf{e}}_{\nu }R\tilde{R}{\mathbf{e}}_{\nu }R\right)\hfill \end{array}$$

(151)

$$\begin{array}{cc}\hfill & ={\rho }^2\underset{\mathrm{metric}\mathrm{tensor}}{\underbrace{{\displaystyle \frac{1}{2}}({\tilde{\mathbf{e}}}_{\mu }{\tilde{\mathbf{e}}}_{\nu }+{\tilde{\mathbf{e}}}_{\nu }{\tilde{\mathbf{e}}}_{\mu })}}\hfill \end{array}$$

(152)

where ${\tilde{\mathbf{e}}}_{\mu }$ and ${\tilde{\mathbf{e}}}_{\nu }$ are SO(3,1) rotated basis vectors.

As one can swap ${\mathbf{e}}_{\mu }$ and ${\mathbf{e}}_{\nu }$ and obtain the same metric tensor, the double-copy product guarantees that $g_{\mu \nu }$ is symmetric.

Furthermore, since ${\mathbf{e}}_{\mu }^{‡}=-{\mathbf{e}}_{\mu }$, we get:

$$\begin{array}{cc}\hfill & {\lfloor {\left({\mathbf{e}}_{\mu }\psi \right)}^{‡}\psi \rfloor }_{3,4}{\left({\mathbf{e}}_{\nu }\psi \right)}^{‡}\psi \hfill \end{array}$$

(153)

$$\begin{array}{cc}\hfill & ={\lfloor {\psi }^{‡}(-1){\mathbf{e}}_{\mu }^{‡}\psi \rfloor }_{3,4}{\psi }^{‡}(-1){\mathbf{e}}_{\nu }^{‡}\psi \hfill \end{array}$$

(154)

$$\begin{array}{cc}\hfill & ={\lfloor {\psi }^{‡}{\mathbf{e}}_{\mu }\psi \rfloor }_{3,4}{\psi }^{‡}{\mathbf{e}}_{\nu }\psi \hfill \end{array}$$

(155)

which allows us to conclude that ${\mathbf{e}}_{\mu }$ and ${\mathbf{e}}_{\nu }$ are self-adjoint within the double-copy product, entailing the interpretation of $g_{\mu \nu }$ as an observable.

In the double-copy product, the metric tensor emerges as a double copy of Dirac currents. This formulation suggests that the metric tensor encodes the probabilistic structure of a quantum theory of gravity in the form of a rank-2 tensor, analogous to how the Dirac current encodes the probabilistic structure of a special relativistic quantum theory in the form of a 4-vector.

Let us now investigate the dynamics. We recall that the evolution operator (Definition 20) is:

$$\begin{array}{c}\hfill T=\left[\begin{array}{c}{e}^{-{\displaystyle \frac{1}{2}}\zeta ({\mathbf{f}}_1+{\mathbf{b}}_1)}\\ & \ddots \\ & & e^{-{\displaystyle \frac{1}{2}}\zeta ({\mathbf{f}}_n+{\mathbf{b}}_n)}\end{array}\right]\end{array}$$

(156)

Acting on the wavefunction, the effect of this operator cascades down to the basis vectors via the double-copy product:

$$\begin{array}{c}\hfill {\lfloor \underset{\mathrm{copy}1}{\underbrace{{\psi }^{‡}{T}^{‡}{\mathbf{e}}_{\mu }T\psi }}\rfloor }_{3,4}\underset{\mathrm{copy}2}{\underbrace{{\psi }^{‡}{T}^{‡}{\mathbf{e}}_{\nu }T\psi }}\end{array}$$

(157)

which realizes an $\mathrm{SO}(3,1)$ transformation of the metric tensor via action of the exponential of a bivector, and a double-copy unitary invariant transformation via action of the exponential of a pseudo-scalar:

$$\begin{array}{c}\hfill {\lfloor \underset{\mathrm{copy}1}{\underbrace{{\psi }^{‡}\underset{\mathrm{SO}(3,1)\mathrm{evolution}}{\underbrace{e^{{\displaystyle \frac{1}{2}}\zeta \mathbf{f}}{\mathbf{e}}_{\mu }{e}^{-{\displaystyle \frac{1}{2}}\zeta \mathbf{f}}}}\underset{\mathrm{unitary}\mathrm{evolution}}{\underbrace{e^{{\displaystyle \frac{1}{2}}\zeta \mathbf{b}}{e}^{-{\displaystyle \frac{1}{2}}\zeta \mathbf{b}}}}\psi }}\rfloor }_{3,4}\underset{\mathrm{copy}2}{\underbrace{{\psi }^{‡}\underset{\mathrm{SO}(3,1)\mathrm{evolution}}{\underbrace{e^{{\displaystyle \frac{1}{2}}\zeta \mathbf{f}}{\mathbf{e}}_{\mu }{e}^{-{\displaystyle \frac{1}{2}}\zeta \mathbf{f}}}}\underset{\mathrm{unitary}\mathrm{evolution}}{\underbrace{e^{{\displaystyle \frac{1}{2}}\zeta \mathbf{b}}{e}^{-{\displaystyle \frac{1}{2}}\zeta \mathbf{b}}}}\psi }}\end{array}$$

(158)

In summary, this initial investigation has identified a scenario in which the metric tensor is measured using basis vectors. The evolution operator, governed by the Schrödinger equation, dynamically realizes SO(3,1) transformations on the metric tensor. Furthermore, the amplitudes associated with possible metric tensors are derived from a double-copy of unitary quantum theories acting on the basis vectors. This formulation simultaneously preserves the SO(3,1) symmetry, essential for describing spacetime structure, and the unitary symmetry, fundamental to quantum mechanics. It describes all changes of basis transformations that an observer in 3+1D spacetime can perform prior to measuring a quantum system.

2.3.5. Starting Point for a Quantum Theory of Gravity

The symmetries of interest can be approached through two distinct strategies:

• Particle Physics Approach: We impose the condition $e^{-{\displaystyle \frac{1}{2}}(\mathbf{f}+\mathbf{b})}{\gamma }_0{e}^{{\displaystyle \frac{1}{2}}(\mathbf{f}+\mathbf{b})}={\gamma }_0$ on the double-copy Dirac current (Equation 158). This constraint leads to the symmetries of the standard model of particle physics, as detailed in Section 2.3.3.
• Gravitational Approach: We allow the double-copy Dirac current (which is equivalent to the metric tensor) to transform freely under SO(3,1). Instead of constraining the current itself, we focus on constructing SO(3,1)-invariant quantities from it. These invariants are typically formed as specific combinations of the double-copy Dirac current and its derivatives. A key example is the Einstein tensor, which remains invariant under SO(3,1) transformations of the metric. The Einstein tensor is particularly significant because it arises from the variation of the Einstein-Hilbert action, which is simplest action leading to such an invariant.

The first strategy provides the gauge symmetries for the standard model of particle physics, while the second offers a path towards a quantum theory of gravity. In the following section, we will explore the gravitational approach in more detail.

2.3.6. Gravitons

Since the double-copy product of the Dirac current holds for any non-degenerate symmetric rank-2 tensor, it should be possible to show that gravitons can be expressed as a special case of this double-copy mechanism. As such, let us now investigate the wave equation in linearized gravity.

It is well known that the Einstein-Hilbert action:

$$\begin{array}{c}\hfill S\left[g_{\mu \nu }\right]={\displaystyle \frac{c^4}{16\pi G}}\int R\sqrt{-\left|g\right|}{\mathrm{d}}^{4}x\end{array}$$

(159)

under the assumption of a small perturbation $g_{\mu \nu }={\eta }_{\mu \nu }+h_{\mu \nu }$, and working in de Donder gauge ${\partial }^{\alpha }{h}_{\alpha \mu }-{\displaystyle \frac{1}{2}}{\partial }_{\mu }h=0$, where $h={\eta }^{\mu \nu }{h}_{\mu \nu }$, can be reduced to its linearized form which is:

$$\begin{array}{c}\hfill S_{EH}^{\left(1\right)}\left[h_{\mu \nu }\right]=\int {\mathrm{d}}^{4}x\left({\displaystyle \frac{1}{2}}{\partial }_{\mu }{h}_{\rho \sigma }{\partial }^{\mu }{h}^{\rho \sigma }-{\displaystyle \frac{1}{4}}{\partial }_{\mu }h{\partial }^{\mu }h\right)\end{array}$$

(160)

Furthermore, varying this action with respect to $h_{\mu \nu }$ and applying the transverse-traceless gauge, yields the wave equation as the equation of motion:

$$\begin{array}{c}\hfill {\displaystyle \frac{\delta S}{\delta h_{\mu \nu }}}=0\Rightarrow \square h_{\mu \nu }=0\end{array}$$

(161)

We now wish to express the wave equation for $h_{\mu \nu }$ in terms of basis vectors ${\mathbf{h}}_{\mu }$ and ${\mathbf{h}}_{\nu }$, such that $h_{\mu \nu }={\displaystyle \frac{1}{2}}({\mathbf{h}}_{\mu }{\mathbf{h}}_{\nu }+{\mathbf{h}}_{\nu }{\mathbf{h}}_{\mu })$. The expression of the wave equation becomes:

$$\begin{array}{c}\hfill \square {\displaystyle \frac{1}{2}}({\mathbf{h}}_{\mu }{\mathbf{h}}_{\nu }+{\mathbf{h}}_{\nu }{\mathbf{h}}_{\mu })=0\end{array}$$

(162)

We identify the solution for ${\mathbf{h}}_{\mu }$ and ${\mathbf{h}}_{\nu }$ by an ansatz:

$$\begin{array}{c}\hfill {\mathbf{h}}_{\mu }(\overrightarrow{x},t)=\mathrm{Re}\int {\mathrm{d}}^{3}k\left(\sqrt{A_{+}}\left(\overrightarrow{k}\right)e_{\mu }^{+}\left(\overrightarrow{k}\right)+\sqrt{A_{\times }}\left(\overrightarrow{k}\right)e_{\mu }^{\times }\left(\overrightarrow{k}\right)\right)e^{-{\displaystyle \frac{1}{2}}i(\overrightarrow{k}·\overrightarrow{x}-{\omega }_{k}t)}\end{array}$$

(163)

$$\begin{array}{c}\hfill {\mathbf{h}}_{\nu }(\overrightarrow{x},t)=\mathrm{Re}\int {\mathrm{d}}^{3}k\left(\sqrt{A_{+}}\left(\overrightarrow{k}\right)e_{\nu }^{+}\left(\overrightarrow{k}\right)+\sqrt{A_{\times }}\left(\overrightarrow{k}\right)e_{\nu }^{\times }\left(\overrightarrow{k}\right)\right)e^{-{\displaystyle \frac{1}{2}}i(\overrightarrow{k}·\overrightarrow{x}-{\omega }_{k}t)}\end{array}$$

(164)

The × and + symbol designed two polarizations.

Then, we promote ${\mathbf{h}}_{\mu }(\overrightarrow{x},t)$ and ${\mathbf{h}}_{\nu }(\overrightarrow{x},t)$ to operators:

$$\begin{array}{c}\hfill {\widehat{\mathbf{h}}}_{\mu }(\overrightarrow{x},t)=\mathrm{Re}\int {\mathrm{d}}^{3}k\left(\sqrt{{\widehat{A}}_{+}}\left(\overrightarrow{k}\right)e_{\mu }^{+}\left(\overrightarrow{k}\right)+\sqrt{{\widehat{A}}_{\times }}\left(\overrightarrow{k}\right)e_{\mu }^{\times }\left(\overrightarrow{k}\right)\right)e^{-{\displaystyle \frac{1}{2}}i(\overrightarrow{k}·\overrightarrow{x}-{\omega }_{k}t)}\end{array}$$

(165)

$$\begin{array}{c}\hfill {\widehat{\mathbf{h}}}_{\nu }(\overrightarrow{x},t)=\mathrm{Re}\int {\mathrm{d}}^{3}k\left(\sqrt{{\widehat{A}}_{+}}\left(\overrightarrow{k}\right)e_{\nu }^{+}\left(\overrightarrow{k}\right)+\sqrt{{\widehat{A}}_{\times }}\left(\overrightarrow{k}\right)e_{\nu }^{\times }\left(\overrightarrow{k}\right)\right)e^{-i{\displaystyle \frac{1}{2}}(\overrightarrow{k}·\overrightarrow{x}-{\omega }_{k}t)}\end{array}$$

(166)

Finally, the probabilities associated to these operators, corresponding to a metric tensor expectation value, are given using the double-copy product as follows:

$$\begin{array}{c}\hfill \langle {\widehat{h}}_{\mu \nu }\rangle ={\displaystyle \frac{1}{2}}(\langle \langle \psi |{\widehat{\mathbf{h}}}_{\mu }\psi \left|\psi \right|{\widehat{\mathbf{h}}}_{\nu }\psi \rangle \rangle +\langle \langle \psi |{\widehat{\mathbf{h}}}_{\nu }\psi \left|\psi \right|{\widehat{\mathbf{h}}}_{\mu }\psi \rangle \rangle )\end{array}$$

(167)

We note that each copy individually applies the Born rule to one of two operators. In contrast, in the conventional perturbative approach to quantum gravity, the metric tensor ${\widehat{h}}_{\mu \nu }$ is quantized, and its expectation value is calculated using the Born rule as follows: $\langle {\widehat{h}}_{\mu \nu }\rangle =\langle \psi |{\widehat{h}}_{\mu \nu }|\psi \rangle$. Thus, our probability calculation differs from the conventional approach, as it involves the product of two separate Born rule applications rather than just one.

The double-copy mechanism, identified by Bern, Carrasco, and Johansson (BCJ) [11], has been shown to simplify calculations of scattering amplitudes for gravitons in many cases. However, it remains an open question whether our double-copy mechanism carries similar advantages.

While we have described gravitons in terms of perturbations $h_{\mu \nu }$ to the metric, it’s important to note that these perturbations transform under SO(3,1) in a way that preserves the SO(3,1) invariance of the Einstein tensor. The Einstein tensor, constructed from $h_{\mu \nu }$ and its derivatives, remains invariant under these transformations. Thus completing (a special case of) the second strategy.

2.4. Dimensional Obstructions

In this section, we explore the dimensional obstructions that arise when attempting to resolve the entropy maximization problem for other dimensional configurations. We found that all geometric configurations except those we have explored here (e.g. $\mathrm{GA}\left(0\right)\cong \mathbb{R}$, $\mathrm{GA}(0,1)\cong \mathbb{C}$ and $\mathrm{GA}(3,1)$) are either obstructed or incomplete. By obstructed, we mean that the solution to the entropy maximization problem, $\rho$, does not satisfy the properties of a probability measure, and by incomplete we refer to the $\mathrm{GA}(2,0)$ case where the metric is not an observable.

$$\begin{array}{cccc}\hfill & \mathit{Dimensions}\hfill & \hfill & \mathit{Obstruction}\hfill \end{array}$$

$$\begin{array}{cccc}\hfill & \mathrm{GA}\left(0\right)\hfill & \hfill & \mathbf{Unobstructed}\Rightarrow \mathrm{statistical}\mathrm{mechanics}\hfill \end{array}$$

(168)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(0,1)\hfill & \hfill & \mathbf{Unobstructed}\Rightarrow \mathrm{quantum}\mathrm{mechanics}\hfill \end{array}$$

(169)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(1,0)\hfill & \hfill & \mathrm{Negative}\mathrm{probabilities}\mathrm{in}\mathrm{the}\mathrm{RQM}\hfill \end{array}$$

(170)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(2,0)\hfill & \hfill & \mathrm{Incomplete}/\mathrm{No}\mathrm{metric}\mathrm{measurement}\hfill \end{array}$$

(171)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(1,1)\hfill & \hfill & \mathrm{Negative}\mathrm{probabilities}\mathrm{in}\mathrm{the}\mathrm{RQM}\hfill \end{array}$$

(172)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(0,2)\hfill & \hfill & \mathrm{Not}\mathrm{isomorphic}\mathrm{to}\mathrm{a}\mathrm{real}\mathrm{matrix}\mathrm{algebra}\hfill \end{array}$$

(173)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(3,0)\hfill & \hfill & \mathrm{Not}\mathrm{isomorphic}\mathrm{to}\mathrm{a}\mathrm{real}\mathrm{matrix}\mathrm{algebra}\hfill \end{array}$$

(174)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(2,1)\hfill & \hfill & \mathrm{Not}\mathrm{isomorphic}\mathrm{to}\mathrm{a}\mathrm{real}\mathrm{matrix}\mathrm{algebra}\hfill \end{array}$$

(175)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(1,2)\hfill & \hfill & \mathrm{Not}\mathrm{isomorphic}\mathrm{to}\mathrm{a}\mathrm{real}\mathrm{matrix}\mathrm{algebra}\hfill \end{array}$$

(176)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(0,3)\hfill & \hfill & \mathrm{Not}\mathrm{isomorphic}\mathrm{to}\mathrm{a}\mathrm{real}\mathrm{matrix}\mathrm{algebra}\hfill \end{array}$$

(177)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(4,0)\hfill & \hfill & \mathrm{Not}\mathrm{isomorphic}\mathrm{to}\mathrm{a}\mathrm{real}\mathrm{matrix}\mathrm{algebra}\hfill \end{array}$$

(178)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(3,1)\hfill & \hfill & \mathbf{Unobstructed}\Rightarrow \mathrm{quantum}\mathrm{gravity}\wedge \mathrm{SU}\left(3\right)\times \mathrm{SU}\left(2\right)\times \mathrm{U}\left(1\right)\hfill \end{array}$$

(179)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(2,2)\hfill & \hfill & \mathrm{Negative}\mathrm{probabilities}\mathrm{in}\mathrm{the}\mathrm{RQM}\hfill \end{array}$$

(180)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(1,3)\hfill & \hfill & \mathrm{Not}\mathrm{isomorphic}\mathrm{to}\mathrm{a}\mathrm{real}\mathrm{matrix}\mathrm{algebra}\hfill \end{array}$$

(181)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(0,4)\hfill & \hfill & \mathrm{Not}\mathrm{isomorphic}\mathrm{to}\mathrm{a}\mathrm{real}\mathrm{matrix}\mathrm{algebra}\hfill \end{array}$$

(182)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(5,0)\hfill & \hfill & \mathrm{Not}\mathrm{isomorphic}\mathrm{to}\mathrm{a}\mathrm{real}\mathrm{matrix}\mathrm{algebra}\hfill \\ \hfill & ⋮\hfill & \hfill & ⋮\hfill \end{array}$$

(183)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(6,0)\hfill & \hfill & \mathrm{No}\mathrm{probability}\mathrm{measure}\mathrm{as}\mathrm{a}\mathrm{self}-\mathrm{product}\hfill \\ \hfill & ⋮\hfill & \hfill & ⋮\hfill \end{array}$$

(184)

$$\begin{array}{cc}\hfill & \infty \hfill \end{array}$$

(185)

Let us now demonstrate the obstructions mentioned above.

Theorem 16 

(Not isomorphic to a real matrix algebra). The determinant of the matrix representation of the geometric algebras in this category is either complex-valued or quaternion-valued, making them unsuitable as a probability.

Proof. 

These geometric algebras are classified as follows:

$$\begin{array}{cc}\hfill & \mathrm{GA}(0,2)\cong \mathbb{H}\hfill \end{array}$$

(186)

$$\begin{array}{cc}\hfill & \mathrm{GA}(3,0)\cong {\mathbb{M}}_2\left(\mathbb{C}\right)\hfill \end{array}$$

(187)

$$\begin{array}{cc}\hfill & \mathrm{GA}(2,1)\cong {\mathbb{M}}_2^2\left(\mathbb{R}\right)\hfill \end{array}$$

(188)

$$\begin{array}{cc}\hfill & \mathrm{GA}(1,2)\cong {\mathbb{M}}_2\left(\mathbb{C}\right)\hfill \end{array}$$

(189)

$$\begin{array}{cc}\hfill & \mathrm{GA}(0,3)\cong {\mathbb{H}}^2\hfill \end{array}$$

(190)

$$\begin{array}{cc}\hfill & \mathrm{GA}(4,0)\cong {\mathbb{M}}_2\left(\mathbb{H}\right)\hfill \end{array}$$

(191)

$$\begin{array}{cc}\hfill & \mathrm{GA}(1,3)\cong {\mathbb{M}}_2\left(\mathbb{H}\right)\hfill \end{array}$$

(192)

$$\begin{array}{cc}\hfill & \mathrm{GA}(0,4)\cong {\mathbb{M}}_2\left(\mathbb{H}\right)\hfill \end{array}$$

(193)

$$\begin{array}{cc}\hfill & \mathrm{GA}(5,0)\cong {\mathbb{M}}_2^2\left(\mathbb{H}\right)\hfill \end{array}$$

(194)

The determinant of these objects is valued in $\mathbb{C}$ or in $\mathbb{H}$, where $\mathbb{C}$ are the complex numbers, and where $\mathbb{H}$ are the quaternions.    □

Theorem 17 

(Negative Probabilities in the RQM). The even sub-algebra, which associates to the RQM part of the theory, of these dimensional configurations allows for negative probabilities, making them unsuitable as a RQM.

Proof. 

This category contains three dimensional configurations:

$\mathrm{GA}(1,0)$:

Let $\psi =a+b{e}_1$, then:

$$\begin{array}{c}\hfill {(a+b{e}_1)}^{‡}(a+b{e}_1)=(a-b{e}_1)(a+b{e}_1)=a^2-b^2{e}_1{e}_1=a^2-b^2\end{array}$$

(195)

which is valued in $\mathbb{R}$.

$\mathrm{GA}(1,1)$:

Let $\psi =a+b{e}_0{e}_1$, then:

$$\begin{array}{c}\hfill {(a+b{e}_0{e}_1)}^{‡}(a+b{e}_0{e}_1)=(a-b{e}_0{e}_1)(a+b{e}_0{e}_1)=a^2-b^2{e}_0{e}_1{e}_0{e}_1=a^2-b^2\end{array}$$

(196)

which is valued in $\mathbb{R}$.

$\mathrm{GA}(2,2)$:

Let $\psi =a+b{e}_0{e}_{\varnothing }{e}_1{e}_2$, where $e_0^2=-1,e_{\varnothing }^2=-1,e_1^2=1,e_2^2=1$, then:

$$\begin{array}{cc}\hfill & {\lfloor {(a+\mathbf{b})}^{‡}(a+\mathbf{b})\rfloor }_{3,4}{(a+\mathbf{b})}^{‡}(a+\mathbf{b})\hfill \end{array}$$

(197)

$$\begin{array}{cc}\hfill & ={\lfloor a^2+2a\mathbf{b}+{\mathbf{b}}^2\rfloor }_{3,4}(a^2+2a\mathbf{b}+{\mathbf{b}}^2)\hfill \end{array}$$

(198)

We note that ${\mathbf{b}}^2=b^2{e}_0{e}_{\varnothing }{e}_1{e}_2{e}_0{e}_{\varnothing }{e}_1{e}_2=b^2$, therefore:

$$\begin{array}{cc}\hfill 1-1& =(a^2+b^2-2a\mathbf{b})(a^2+b^2+2a\mathbf{b})\hfill \end{array}$$

(199)

$$\begin{array}{cc}\hfill & ={(a^2+b^2)}^2-4a^2{\mathbf{b}}^2\hfill \end{array}$$

(200)

$$\begin{array}{cc}\hfill & ={(a^2+b^2)}^2-4a^2{b}^2\hfill \end{array}$$

(201)

which is valued in $\mathbb{R}$.

In all of these cases the RQM probability can be negative.    □

We repeat the following self-products[8] (Definition 16), which will help us demonstrate the next theorem:

$$\begin{array}{cccc}\hfill & \mathrm{GA}(0,1):\hfill & \hfill & {\phi }^{†}\phi \hfill \end{array}$$

(202)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(2,0):\hfill & \hfill & {\phi }^{‡}\phi \hfill \end{array}$$

(203)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(3,0):\hfill & \hfill & {\lfloor {\phi }^{‡}\phi \rfloor }_3{\phi }^{‡}\phi \hfill \end{array}$$

(204)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(3,1):\hfill & \hfill & {\lfloor {\phi }^{‡}\phi \rfloor }_{3,4}{\phi }^{‡}\phi \hfill \end{array}$$

(205)

$$\begin{array}{cccc}\hfill & \mathrm{GA}(4,1):\hfill & \hfill & {\left({\lfloor {\phi }^{‡}\phi \rfloor }_{3,4}{\phi }^{‡}\phi \right)}^{†}\left({\lfloor {\phi }^{‡}\phi \rfloor }_{3,4}{\phi }^{‡}\phi \right)\hfill \end{array}$$

(206)

Theorem 18 

(No Metric Measurements). This obstruction applies to $\mathrm{GA}(2,0)$. A probability measure of at least four self-products are required for the theory to be observationally complete with respect to its geometry.

Proof. 

A metric measurement requires a probability measure of 4 self products because the metric tensor is defined using 2 self-products of the gamma matrices:

$$\begin{array}{c}\hfill g_{\mu \nu }={\displaystyle \frac{1}{2}}({\mathbf{e}}_{\mu }{\mathbf{e}}_{\nu }+{\mathbf{e}}_{\nu }{\mathbf{e}}_{\mu })\end{array}$$

(207)

Each pair of wavefunction products fixes one basis elements. Thus, two pairs of wavefunction products are required to fix the geometry from the wavefunction. As probability measures of four self-products begin to appear in 3D, then the $\mathrm{GA}(2,0)$ cannot produce a metric measurement as a quantum observable, thus its geometry is not observationally complete with respect to its geometry.    □

Theorem 18 

(No Metric Measurements). The multivector representation of the norm in 6D cannot satisfy any observables.

Proof 

(Argument). In six dimensions and above, the self-product patterns found in Definition 16 collapse. The research by Acus et al.[12] in 6D geometric algebra demonstrates that the determinant, so far defined through a self-products of the multivector, fails to extend into 6D. The crux of the difficulty is evident in the reduced case of a 6D multivector containing only scalar and grade-4 elements:

$$\begin{array}{c}\hfill s\left(B\right)=b_{1}B{f}_5\left(f_4\left(B\right)f_3\left(f_2\left(B\right)f_1\left(B\right)\right)\right)+b_{2}B{g}_5\left(g_4\left(B\right)g_3\left(g_2\left(B\right)g_1\left(B\right)\right)\right)\end{array}$$

(208)

This equation is not a multivector self-product but a linear sum of two multivector self-products[12].

The full expression is given in the form of a system of 4 equations, which is too long to list in its entirety. A small characteristic part is shown:

$$\begin{array}{cc}\hfill & a_0^4-2a_0^2{a}_{47}^2+b_2{a}_0^2{a}_{47}^2{p}_{412}{p}_{422}+\langle 72\mathrm{monomials}\rangle =0\hfill \end{array}$$

(209)

$$\begin{array}{cc}\hfill & b_1{a}_0^3{a}_{52}+2b_2{a}_0{a}_{47}^2{a}_{52}{p}_{412}{p}_{422}{p}_{432}{p}_{442}{p}_{452}+\langle 72\mathrm{monomials}\rangle =0\hfill \end{array}$$

(210)

$$\begin{array}{cc}\hfill & \langle 74\mathrm{monomials}\rangle =0\hfill \end{array}$$

(211)

$$\begin{array}{cc}\hfill & \langle 74\mathrm{monomials}\rangle =0\hfill \end{array}$$

(212)

From Equation 208, it is possible to see that no observable $\mathbf{O}$ can satisfy this equation because the linear combination does not allow one to factor it out of the equation.

$$\begin{array}{c}\hfill b_1\mathbf{O}B{f}_5\left(f_4\left(B\right)f_3\left(f_2\left(B\right)f_1\left(B\right)\right)\right)+b_{2}B{g}_5\left(g_4\left(B\right)g_3\left(g_2\left(B\right)g_1\left(B\right)\right)\right)=b_{1}B{f}_5\left(f_4\left(B\right)f_3\left(f_2\left(B\right)f_1\left(B\right)\right)\right)+b_2\mathbf{O}B{g}_5\left(g_4\left(B\right)g_3\left(g_2\left(B\right)g_1\left(B\right)\right)\right)\end{array}$$

(213)

Any equality of the above type between $b_1\mathbf{O}$ and $b_2\mathbf{O}$ is frustrated by the factors $b_1$ and $b_2$, forcing $\mathbf{O}=1$ as the only satisfying observable. Since the obstruction occurs within grade-4, which is part of the even sub-algebra it is questionable that a satisfactory theory (with non-trivial observables) be constructible in 6D, suing this method.    □

This conjecture proposes that the multivector representation of the determinant in 6D does not allow for the construction of non-trivial observables, which is a crucial requirement for a relevant quantum formalism. The linear combination of multivector self-products in the 6D expression prevents the factorization of observables, limiting their role to the identity operator.

Conjecture 2 

(No probability measures as a self-product (above 6D)). The norms beyond 6D are progressively more complex than the 6D case, which is already obstructed.

These theorems and conjectures provide additional insights into the unique role of the unobstructed 3+1D signature in our proposal.

It is also interesting that our proposal is able to rule out $\mathrm{GA}(1,3)$ even if in relativity, the signature of the metric $(+,-,-,-)$ versus $(-,-,-,+)$ does not influence the physics. However, in geometric algebra, $\mathrm{GA}(1,3)$ represents 1 space dimension and 3 time dimensions. Therefore, it is not the signature itself that is ruled out but rather the specific arrangement of 3 time and 1 space dimensions, as this configuration yields quaternion-valued "probabilities" (i.e. $\mathrm{GA}(1,3)\cong {\mathbb{M}}_2\left(\mathbb{H}\right)$ and $det{\mathbb{M}}_2\left(\mathbb{H}\right)\in \mathbb{H}$).

Consequently, 3+1D is the only dimensional configuration (other than the "non-geometric" configurations of $\mathrm{GA}\left(0\right)\cong \mathbb{R}$ and $\mathrm{GA}(0,1)\cong \mathbb{C}$) in which a ’least biased’ solution to the problem of maximizing the Shannon entropy of quantum measurements relative to an initial preparation, exists. This is an extremely constraining result regarding the possible spacetime configurations of the universe, and our ability (or inability) to construct a least biased theory to investigate it.

3. Discussion

The principle of maximum entropy[3] states that the probability measure that best represents the current state of knowledge about a system is the one with the largest entropy, constrained by prior data. In QM, an experiment typically consists of three stages: an initial preparation, followed by some transformations, and concluding with a final measurement of the system, which yields the result of the experiment. Consistent with the maximum entropy principle, our aim is to derive the ’least biased’ theory that connects the initial preparation p to its final measurement $\rho$. By formulating the theory as a solution to a maximization problem, rather than merely by axiomatic stipulation, we ensure (by mathematical proof) that the resulting framework is as unbiased as possible given the available information.

Using this methodology, fundamental physics can be formulated as the general solution to a maximization problem involving the Shannon entropy of all possible measurements of an arbitrary system relative to its initial preparation, subject to a vanishing phase constraint. The structure of the inferred theory is thus determined by the nature and generality of the employed constraint. In this paper, we have investigated three specific entropy maximization problems, each characterized by a different constraint and corresponding to a distinct level of description in physics:

$$\begin{array}{cccccccc}\hfill & \mathit{Constraint}\hfill & \hfill & \mathit{Vanishing}\mathit{Phase}\hfill & \hfill & \mathit{Inferred}\mathit{Theory}\hfill & \hfill & \mathit{Wavefunction}\hfill \\ \hfill & \overline{E}=\sum _i{\rho }_i{E}_i\hfill & \hfill & \mathrm{none}\hfill & \hfill & \mathrm{SM}\hfill & \hfill & \mathrm{N}.\mathrm{A}.\hfill \\ \hfill & 0=\mathrm{tr}\sum _i{\rho }_i\left[\begin{array}{cc}0& -E_i\\ E_i& 0\end{array}\right]\hfill & \hfill & \mathrm{U}\left(1\right)\hfill & \hfill & \mathrm{QM}\hfill & \hfill & {\mathbb{C}}^n\hfill \\ \hfill & 0={\displaystyle \frac{1}{2}}\mathrm{tr}\sum _i{\rho }_i{\mathbf{M}}_i\hfill & \hfill & {\mathrm{Spin}}^{\mathrm{c}}(3,1)\hfill & \hfill & \mathrm{RQM}/\mathrm{QG}\hfill & \hfill & {(\mathbb{R}\times {\mathrm{Spin}}^{\mathrm{c}}(3,1))}^n\hfill \end{array}$$

where n represents the size of the ensemble.

Despite the differences in constraints, the three theories here-so formulated share a common logical genesis, adhere to the same principle of maximum entropy, and qualify as the least biased theory for their given constraint.

3.1. Guarantee of Epistemic Soundness

The Born rule is the least biased probability measure operating on a complex Hilbert space, as established in Theorem 2. However, when the framework is extended to include geometric considerations, the Born rule is no longer the least biased measure. Instead, as demonstrated in Theorem 3, the double-copy product emerges as the least biased probability measure in the 3+1D setting. Interestingly, no solutions exist for other geometric configurations, as demonstrated in Section 2.4.

3.2. Guarantee of Ontological Soundness

Our approach inverts the traditional theory construction paradigm: instead of postulating abstract mathematical entities like wavefunctions and Hilbert spaces, we derive them from the measurements—the theory’s only constraint. Remarkably, the logical foundation of the theory is its ontology. This guarantees that the theory’s foundation matches what transpired in the lab and also renders it inherently resistant to falsification by the very measurements used in its construction.

4. Conclusion

This paper presents a novel approach to quantum theory construction by solving a maximization problem on the Shannon entropy of all possible measurements of a system relative to its initial preparation, under the constraint of a vanishing phase. By selecting the appropriate group of the vanishing phase, the solution resolves to quantum mechanics, relativistic quantum mechanics, or a candidate for a theory of quantum gravity. The resulting measure is invariant under a wide range of geometric transformations, including those generated by the gauge groups of the Standard Model, and leads to the metric tensor as an operator involving a double copy of Dirac currents, without additional assumptions. Again without additional assumptions, the theory automatically fails in geometric configurations other than 3+1D. This result aligns with the observed dimensionality and gauge symmetries of the universe, suggesting a possible explanation for its specific structure. This approach may offer a promising avenue for the unification of fundamental physical theories and may provide new insights into the underlying principles governing the structure of our universe.