Physics as an Optimization Problem over All Measurements

Preprint

Article

Physics as an Optimization Problem over All Measurements

Altmetrics

Downloads

449

Views

186

Comments

Alexandre Harvey-Tremblay^*

This version is not peer-reviewed

Submitted:

21 June 2024

Posted:

25 June 2024

Read the latest preprint version here

Alerts

Abstract

We propose a novel approach to quantum theory construction that involves solving a maximization problem on the Shannon entropy of all possible measurements of a system, relative to its initial preparation. This maximization problem is additionally constrained by a phase condition that vanishes under measurements. Specifically, enforcing a vanishing U(1)-valued phase constraint leads to standard quantum mechanics, while a vanishing Spin^c(3,1)-valued phase constraint extends the theory to relativistic quantum mechanics and to quantum gravity. The latter scenario is found to incorporate the SU(3)xSU(2)xU(1) symmetries of the Standard Model as invariants of the probability measure itself, and to construct the metric tensor as an operator via a double-copy mechanism applied to Dirac currents. Significantly, this solution is consistent exclusively with a 3+1-dimensional spacetime configuration---all other dimensional settings are shown to lead to fundamental obstructions. This framework seamlessly integrates fundamental concepts from quantum mechanics, relativistic quantum mechanics, quantum gravity, the dimensional specificity of spacetime, and particle physics symmetries as the solution to a simple entropy maximization problem.

Keywords:

Subject: Physical Sciences - Quantum Science and Technology

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 ( $\bar{E}$ ):

$\begin{matrix} \bar{E} = \sum_{i} ρ_{i} E_{i} \end{matrix}$

(1)

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

Definition 1

(Fundamental Lagrange Multiplier Equation of SM).

\begin{matrix} L = \underset{Boltzmann entropy}{\underset{︸}{- k_{B} \sum_{i} ρ_{i} ln ρ_{i}}} + \underset{N o r m a l i z a t i o n C o n s t r a i n t}{\underset{︸}{λ (1 - \sum_{i} ρ_{i})}} + \underset{A v e r a g e E n e r g y C o n s t r a i n t}{\underset{︸}{β (\bar{E} - \sum_{i} ρ_{i} E_{i})}} \end{matrix}

(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{matrix} ρ_{i} = \underset{Microcanonical Ensemble}{\underset{︸}{\frac{1}{\sum_{i} exp (- β E_{i})}}} exp (- β E_{i}) \end{matrix}

(3)

Proof.

This is an well-known result by E. T. Jaynes [3,4]. As a convenience, we replicate the proof in Annex 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{matrix} 0 = tr \sum_{i} ρ_{i} [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}] \end{matrix}$

(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{matrix} L = \underset{R e l a t i v e S h a n n o n E n t r o p y}{\underset{︸}{- \sum_{i} ρ_{i} ln \frac{ρ_{i}}{p_{i}}}} + \underset{N o r m a l i z a t i o n C o n s t r a i n t}{\underset{︸}{λ (1 - \sum_{i} ρ_{i})}} + \underset{V a n i s h i n g C o m p l e x - P h a s e}{\underset{︸}{τ (- tr \sum_{i} ρ_{i} [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}])}} \end{matrix}

(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

ρ

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{matrix} ρ_{i} = \underset{Unitarily Invariant Ensemble}{\underset{︸}{\frac{1}{\sum_{i} p_{i} ∥ exp (- i t E_{i} / ℏ) ∥}}} \underset{Born Rule}{\underset{︸}{∥ exp (- i t E_{i} / ℏ) ∥}} \underset{Initial Preparation}{\underset{︸}{p_{i}}} \end{matrix}

(6)

where we have defined

τ = t / ℏ

(analogous to

β = 1 / (k_{B} T)

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

u = a + x + f + v + b

, where a is a scalar,

x

is a vector,

f

is a bivector,

v

is a pseudo-vector and

b

is a pseudo-scalar, be a multivector of the geometric algebra

GA (3, 1)

, and let

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 ${Spin}^{c} (3, 1)$ group. The constraint is:

$\begin{matrix} 0 = tr \frac{1}{2} \sum_{i} ρ_{i} M_{i} \end{matrix}$

(7)

where $M_{i}$ is the matrix representation of the multivector $u_{i} = f_{i} + b_{i}$ of $GA (3, 1)$ . Using the real Majorana representation of the gamma matrices, the representation is as follows:

$\begin{matrix} u ≅ M = [\begin{matrix} 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{matrix}] \end{matrix}$

(8)

Similarly to the QM case, here the matrix representation engenders a ${Spin}^{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{matrix} L = \underset{R e l a t i v e S h a n n o n E n t r o p y}{\underset{︸}{- \sum_{i} ρ_{i} ln \frac{ρ_{i}}{p_{i}}}} + \underset{N o r m a l i z a t i o n C o n s t r a i n t}{\underset{︸}{λ (1 - \sum_{i} ρ_{i})}} + \underset{V a n i s h i n g S p i n^{c} (3, 1) P h a s e}{\underset{︸}{ζ (- tr \frac{1}{2} \sum_{i} ρ_{i} M_{i})}} \end{matrix}

(9)

where λ and ζ are the Lagrange multipliers.

Theorem 3.

The least biased probability measure that connects an initial preparation

p_{i}

to its final measurement

ρ_{i}

, under the constraint of the vanishing relativistic phase, is:

\begin{matrix} ρ (q) = \underset{Spin ​^{c} (3, 1) Invariant Ensemble}{\underset{︸}{\frac{1}{\sum_{i} p_{i} det exp (- ζ \frac{1}{2} M_{i})}}} \underset{Spin ​^{c} (3, 1) Born Rule}{\underset{︸}{det exp (- ζ \frac{1}{2} M_{i})}} \underset{Initial Preparation}{\underset{︸}{p_{i}}} \end{matrix}

(10)

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

In the results section, we aim to demonstrate that this solution establishes the foundation for a quantum mechanical theory of observables transforming in 3+1-dimensional spacetime. In this framework,

ζ

acts as a one-parameter generator for boosts, rotations, and phase transformations.

Incorporating additional structures, such as the Dirac equation, serves as a constraint that narrows the space of possible wavefunctions to those that satisfy it. This constitutes a specific realization of a RQM theory within this general framework.

1.3. Quantum Gravity

The maximization problem employed for RQM, as described in Section 1.2, will also be found to encompass QG. While RQM was situated within a subset of the solution space, fully leveraging the generality of the solution is necessary for QG. Specifically, the solution will be found to incorporate a double-copy mechanism applied to the Dirac current, such that two Dirac currents are multiplied together to yield a metric tensor expectation value. In this context, a single copy of the Dirac current pertains to RQM, whereas the double-copy leads to QG. The double-copy of the Dirac current facilitates the construction of the metric tensor as an observable from basis vectors, thereby establishing a connection between the probability measure and the geometrical 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

U (1)

phase), the entropy maximization technique yields a solution only in 3+1-dimensional spacetime (vanishing

{Spin}^{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

U (1)

-Phase Anti-Constraint:

\begin{matrix} 0 = tr \sum_{i} ρ_{i} [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}] \end{matrix}

(11)

where

E_{i}

are scalar-valued. The usage of the matrix generates a

U (1)

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

ρ_{i}

. For that, we deploy the following Lagrange multiplier equation:

\begin{matrix} L = \underset{Relative Shannon Entropy}{\underset{︸}{- \sum_{i} ρ_{i} ln \frac{ρ_{i}}{p_{i}}}} + \underset{Normalization Constraint}{\underset{︸}{λ (1 - \sum_{i} ρ_{i})}} + \underset{Vanishing U (1) - Phase}{\underset{︸}{τ (tr \sum_{i} ρ_{i} [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}])}} \end{matrix}

(12)

where

λ

and

τ

are the Lagrange multipliers.

We solve the maximization problem as follows:

\begin{matrix} 0 & = \frac{\partial L [ρ_{1}, \dots, ρ_{n}]}{\partial ρ_{i}} \end{matrix}

(13)

\begin{matrix} = - ln \frac{ρ_{i}}{p_{i}} - p_{i} - λ - τ tr [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}] \end{matrix}

(14)

\begin{matrix} = ln \frac{ρ_{i}}{p_{i}} + p_{i} + λ - τ tr [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}] \end{matrix}

(15)

\begin{matrix} \Rightarrow ln \frac{ρ_{i}}{p_{i}} & = - p_{i} - λ - τ tr [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}] \end{matrix}

(16)

\begin{matrix} \Rightarrow ρ_{i} & = p_{i} exp (- p_{i} - λ) exp (- τ tr [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}]) \end{matrix}

(17)

\begin{matrix} = \frac{1}{Z (τ)} p_{i} exp (- τ tr [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}]) \end{matrix}

(18)

The partition function is obtained as follows:

\begin{matrix} 1 & = \sum_{i} p_{i} exp (- p_{i} - λ) exp (- τ tr [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}]) \end{matrix}

(19)

\begin{matrix} \Rightarrow {(exp (- p_{i} - λ))}^{- 1} & = \sum_{i} p_{i} exp (- τ tr [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}]) \end{matrix}

(20)

\begin{matrix} Z (τ) & : = \sum_{i} p_{i} exp (- τ tr [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}]) \end{matrix}

(21)

Finally, the least biased probability measure that connects an initial preparation

p_{i}

to its final measurement

ρ_{i}

, under the constraint of the vanishing U(1) phase, is:

\begin{matrix} ρ_{i} = \frac{1}{\sum_{i} p_{i} exp (- τ tr [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}])} exp (- τ tr [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}]) p_{i} \end{matrix}

(22)

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{matrix} ρ_{i} = \frac{p_{i}}{\sum_{i} p_{i}} \end{matrix}

(23)

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

[\begin{matrix} a & - b \\ b & a \end{matrix}]

a + i b

. For instance, the exponential of a complex matrix is:

\begin{matrix} exp [\begin{matrix} a & - b \\ b & a \end{matrix}] & = r [\begin{matrix} cos (b) & - sin (b) \\ sin (b) & cos (b) \end{matrix}], where r = exp a \end{matrix}

(24)

Then, we associate the exponential trace to the complex norm using

exp tr M \equiv det exp M

\begin{matrix} exp tr [\begin{matrix} a & - b \\ b & a \end{matrix}] = det exp [\begin{matrix} a & - b \\ b & a \end{matrix}] & = r^{2} det [\begin{matrix} cos (b) & - sin (b) \\ sin (b) & cos (b) \end{matrix}], where r = exp a \end{matrix}

(25)

\begin{matrix} = r^{2} ({cos}^{2} (b) + {sin}^{2} (b)) \end{matrix}

(26)

\begin{matrix} = ∥ r (cos (b) + i sin (b)) ∥ \end{matrix}

(27)

\begin{matrix} = ∥ r exp (i b) ∥ \end{matrix}

(28)

Finally, substituting

τ = t / ℏ

analogously to

β = 1 / (k_{B} T)

, and applying the complex-norm representation to both the numerator and to the denominator, consolidates the Born rule, normalization, and initial prepration into :

\begin{matrix} ρ_{i} & = \underset{Unitarily Invariant Partition Function}{\underset{︸}{\frac{1}{\sum_{i} p_{i} ∥ exp (- i t E_{i} / ℏ) ∥}}} \underset{Born Rule}{\underset{︸}{∥ exp (- i t E_{i} / ℏ) ∥}} \underset{Initial Preparation}{\underset{︸}{p_{i}}} \end{matrix}

(29)

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{matrix} \sum_{i} p_{i} ∥ exp (- i t E_{i} / ℏ) ∥ = Z = 〈 ψ | ψ 〉 \end{matrix}

(30)

where

\begin{matrix} [\begin{matrix} ψ_{1} (t) \\ ⋮ \\ ψ_{n} (t) \end{matrix}] = [\begin{matrix} exp (- i t E_{1} / ℏ) \\ ⋱ \\ exp (- i t E_{n} / ℏ) \end{matrix}] [\begin{matrix} ψ_{1} (0) \\ ⋮ \\ ψ_{n} (0) \end{matrix}] \end{matrix}

(31)

We clarify that

p_{i}

represents the probability associated with the initial preparation of the wavefunction, where

p_{i} = 〈 ψ_{i} (0) | ψ_{i} (0) 〉

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 $| ψ 〉$ with $1 / Z = 1 / \sqrt{〈 ψ | ψ 〉}$ . This normalization links $| ψ 〉$ 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{matrix} \bar{O} = \sum_{i} p_{i} O_{i} ∥ exp (- i t E_{i} / ℏ) ∥ \end{matrix}$

(32)

Since $Z = 〈 ψ | ψ 〉$ , then any self-adjoint operator satisfying the condition $〈 O ψ | ϕ 〉 = 〈 ψ | O ϕ 〉$ will equate the above equation, simply because $〈 O 〉 = 〈 ψ | O | ψ 〉$ . This demonstrates QM Axiom 2 of 5.
Upon transforming Equation 31 out of its eigenbasis through unitary operations, we find that the energy, $E_{i}$ , typically transforms in the manner of a Hamiltonian operator:

$\begin{matrix} | ψ (t) 〉 = exp (- i t H / ℏ) | ψ (0) 〉 \end{matrix}$

(33)

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

$\begin{matrix} \frac{\partial}{\partial t} | ψ (t) 〉 & = \frac{\partial}{\partial t} (exp (- i t H / ℏ) | ψ (0) 〉) \end{matrix}$

(34)

$\begin{matrix} = - i H / ℏ exp (- i t H / ℏ) | ψ (0) 〉 \end{matrix}$

(35)

$\begin{matrix} = - i H / ℏ | ψ (t) 〉 \end{matrix}$

(36)

$\begin{matrix} \Rightarrow H | ψ (t) 〉 & = i ℏ \frac{\partial}{\partial t} | ψ (t) 〉 \end{matrix}$

(37)

which is the Schrödinger equation. This demonstrates QM Axiom 3 of 5.
From Equation 31 it follows that the possible microstates $E_{i}$ of the system correspond to specific eigenvalues of $H$ . An observation can thus be conceptualized as sampling from $ρ$ , 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 $H$ . Measured in the eigenbasis, the probability measure is:

$\begin{matrix} ρ_{i} (t) = \frac{1}{〈 ψ | ψ 〉} {(ψ_{i} (t))}^{†} ψ_{i} (t) . \end{matrix}$

(38)

In scenarios where the probability measure $ρ_{i} (τ)$ is expressed in a basis other than its eigenbasis, the probability $P (λ_{i})$ of obtaining the eigenvalue $λ_{i}$ is given as a projection on a eigenstate:

$\begin{matrix} P (λ_{i}) = {| 〈 λ_{i} | ψ 〉 |}^{2} \end{matrix}$

(39)

Here, $| 〈 λ_{i} | ψ 〉 |^{2}$ signifies the squared magnitude of the amplitude of the state $| ψ 〉$ when projected onto the eigenstate $| λ_{i} 〉$ . As this argument hold for any observables, this demonstrates QM Axiom 4 of 5.
Finally, since the probability measure (Equation 29) 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 11) is sufficient to entail the foundations of quantum mechanics (Axiom 1, 2, 3, 4 and 5) through the principle of entropy maximization. Equation 11 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 dimensional 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{matrix} L = \underset{Relative Shannon Entropy}{\underset{︸}{- \sum_{i} ρ_{i} ln \frac{ρ_{i}}{p_{i}}}} + \underset{Normalization Constraint}{\underset{︸}{λ (1 - \sum_{i} ρ_{i})}} + \underset{Vanishing Relativistic Phase}{\underset{︸}{θ (- tr \frac{1}{2} \sum_{i} ρ_{i} M_{i})}} \end{matrix}

(40)

where

λ

and

θ

are the Lagrange multipliers, and where

M_{i}

is the matrix representation of a multivector

b_{i}

GA (2)

, where

b_{i}

is a pseudo-scalar. In general a multivector

u = a + x + b

GA (2)

, where a is a scalar,

x

is a vector and

b

a pseudo-scalar, is represented as follows:

\begin{matrix} [\begin{matrix} a + x & y - b \\ y + b & a - x \end{matrix}] ≅ a + x σ_{x} + y σ_{y} + b σ_{x} \land σ_{y} \end{matrix}

(41)

The basis elements are defined as:

\begin{matrix} σ_{x} = [\begin{matrix} 1 & 0 \\ 0 & - 1 \end{matrix}], σ_{y} = [\begin{matrix} 0 & 1 \\ 1 & 0 \end{matrix}], σ_{x} \land σ_{y} = [\begin{matrix} 0 & - 1 \\ 1 & 0 \end{matrix}] \end{matrix}

(42)

If we take

a \to 0, x \to 0

then

M

reduces as follows:

\begin{matrix} {u = a + x + b |}_{a \to 0, x \to 0} = b \Rightarrow M = [\begin{matrix} 0 & - b \\ b & 0 \end{matrix}] \end{matrix}

(43)

The Lagrange multiplier equation can be solved as follows:

\begin{matrix} 0 & = \frac{\partial L [ρ_{1}, \dots, ρ_{n}]}{\partial ρ_{i}} \end{matrix}

(44)

\begin{matrix} = - ln \frac{ρ_{i}}{p_{i}} - p_{i} - λ - θ tr \frac{1}{2} [\begin{matrix} 0 & - b_{i} \\ b_{i} & 0 \end{matrix}] \end{matrix}

(45)

\begin{matrix} = ln \frac{ρ_{i}}{p_{i}} + p_{i} + λ + θ tr \frac{1}{2} [\begin{matrix} 0 & - b_{i} \\ b_{i} & 0 \end{matrix}] \end{matrix}

(46)

\begin{matrix} \Rightarrow ln \frac{ρ_{i}}{p_{i}} & = - p_{i} - λ - θ tr \frac{1}{2} [\begin{matrix} 0 & - b_{i} \\ b_{i} & 0 \end{matrix}] \end{matrix}

(47)

\begin{matrix} \Rightarrow ρ_{i} & = p_{i} exp (- p_{i} - λ) exp (- θ tr \frac{1}{2} [\begin{matrix} 0 & - b_{i} \\ b_{i} & 0 \end{matrix}]) \end{matrix}

(48)

\begin{matrix} = \frac{1}{Z (θ)} p_{i} exp (- θ tr \frac{1}{2} [\begin{matrix} 0 & - b_{i} \\ b_{i} & 0 \end{matrix}]) \end{matrix}

(49)

The partition function

Z (θ)

, serving as a normalization constant, is determined as follows:

\begin{matrix} 1 & = \sum_{i} p_{i} exp (- p_{i} - λ) exp (- θ tr \frac{1}{2} [\begin{matrix} 0 & - b_{i} \\ b_{i} & 0 \end{matrix}]) \end{matrix}

(50)

\begin{matrix} \Rightarrow {(exp (- p_{i} - λ))}^{- 1} & = \sum_{i} p_{i} exp (- θ tr \frac{1}{2} [\begin{matrix} 0 & - b_{i} \\ b_{i} & 0 \end{matrix}]) \end{matrix}

(51)

\begin{matrix} Z (θ) & : = \sum_{i} p_{i} exp (- θ tr \frac{1}{2} [\begin{matrix} 0 & - b_{i} \\ b_{i} & 0 \end{matrix}]) \end{matrix}

(52)

Consequently, the least biased probability measure that connects an initial preparation

p_{i}

to a final measurement

ρ_{i}

, under the constraint of the vanishing relativistic phase in 2D is:

\begin{matrix} ρ_{i} = \underset{Spin (2) Invariant Ensemble}{\underset{︸}{\frac{1}{\sum_{i} p_{i} det exp (- \frac{1}{2} θ [\begin{matrix} 0 & - b_{i} \\ b_{i} & 0 \end{matrix}])}}} \underset{Spin (2) Born Rule}{\underset{︸}{det exp (- \frac{1}{2} θ [\begin{matrix} 0 & - b_{i} \\ b_{i} & 0 \end{matrix}])}} \underset{Initial Preparation}{\underset{︸}{p_{i}}} \end{matrix}

(53)

where

det exp M = exp tr M

In 2D, the Lagrange multiplier

θ

correspond to an angle of rotation, and in 1+1D it would correspond to the rapidity

ζ

\begin{matrix} 2 D : & exp θ [\begin{matrix} 0 & - 1 \\ 1 & 0 \end{matrix}] = [\begin{matrix} cos θ & - sin θ \\ sin θ & cos θ \end{matrix}] & θ is the angle of rotation \end{matrix}

(54)

\begin{matrix} 1 + 1 D : & exp ζ [\begin{matrix} 0 & 1 \\ 1 & 0 \end{matrix}] = [\begin{matrix} cosh ζ & sinh ζ \\ sinh ζ & cosh ζ \end{matrix}] & ζ is the rapidity \end{matrix}

(55)

The 2D solution may appear equivalent to the QM case because they are related by an isomorphism

Spin (2) ≅ SO (2) ≅ U (1)

and under the replacement

θ \to τ

. 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

u = a + x + b

be a multi-vector of the geometric algebra over the reals in two dimensions

GA (2)

. The multivector conjugate is defined as:

\begin{matrix} u^{‡} = a - x - b \end{matrix}

(56)

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{matrix} u^{‡} u = det M \end{matrix}

(57)

Proof.

Let

u = a + x σ_{x} + y σ_{y} + b σ_{x} \land σ_{y}

, and let

M

be its matrix representation

[\begin{matrix} a + x & y - b \\ y + b & a - x \end{matrix}]

. Then:

\begin{matrix} 1 : & u^{‡} u \end{matrix}

(58)

\begin{matrix} = {(a + x σ_{x} + y σ_{y} + b σ_{x} \land σ_{y})}^{‡} (a + x σ_{x} + y σ_{y} + b σ_{x} \land σ_{y}) \end{matrix}

(59)

\begin{matrix} = (a - x σ_{x} - y σ_{y} - b σ_{x} \land σ_{y}) (a + x σ_{x} + y σ_{y} + b σ_{x} \land σ_{y}) \end{matrix}

(60)

\begin{matrix} = a^{2} - x^{2} - y^{2} + b^{2} \end{matrix}

(61)

\begin{matrix} 2 : & det M \end{matrix}

(62)

\begin{matrix} = det [\begin{matrix} a + x & y - b \\ y + b & a - x \end{matrix}] \end{matrix}

(63)

\begin{matrix} = (a + x) (a - x) - (y - b) (y + b) \end{matrix}

(64)

\begin{matrix} = a^{2} - x^{2} - y^{2} + b^{2} \end{matrix}

(65)

□

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

|V〉 {〉 \in (GA (2))}^{n}

\begin{matrix} |V〉 〉 = [\begin{matrix} a_{1} + x_{1} + b_{1} \\ ⋮ \\ a_{n} + x_{n} + b_{n} \end{matrix}] \end{matrix}

(66)

The multivector conjugate transpose of

|V〉 〉

is defined as first taking the transpose and then the element-wise multivector conjugate:

\begin{matrix} 〈 〈V| = [\begin{matrix} a_{1} - x_{1} - b_{1} & \dots & a_{n} - x_{n} - b_{n} \end{matrix}] \end{matrix}

(67)

Definition 6

(Bilinear Form). Let

|V〉 〉

and

|W〉 〉

be two vectors valued in

GA (2)

. We introduce the following bilinear form:

\begin{matrix} 〈 〈V | W〉 〉 = (a_{1} - x_{1} - b_{1}) (a_{1} + x_{1} + b_{1}) + \dots (a_{n} - x_{n} - b_{n}) (a_{n} + x_{n} + b_{n}) \end{matrix}

(68)

Theorem 5

(Inner Product). Restricted to the even sub-algebra of

GA (2)

, the bilinear form is an inner product.

Proof.

\begin{matrix} {〈 〈V | W〉 〉}_{x \to 0} & = (a_{1} - b_{1}) (a_{1} + b_{1}) + \dots (a_{n} - b_{n}) (a_{n} + b_{n}) \end{matrix}

(69)

This is isomorphic to the inner product of a complex Hilbert space, with the identification

i ≅ σ_{x} \land σ_{y}

. □

Definition 7.

(Spin(2)-valued Wavefunction).

\begin{matrix} |ψ〉 〉 = [\begin{matrix} e^{\frac{1}{2} (a_{1} + b_{1})} \\ ⋮ \\ e^{\frac{1}{2} (a_{n} + b_{n})} \end{matrix}] = [\begin{matrix} \sqrt{ρ_{1}} R_{1} \\ ⋮ \\ \sqrt{ρ_{2}} R_{2} \end{matrix}] \end{matrix}

(70)

where

\sqrt{ρ_{i}} = e^{\frac{1}{2} a_{i}}

representing the square root of the probability and

R_{i} = e^{\frac{1}{2} 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 = 〈 〈ψ | ψ〉 〉

Proof.

\begin{matrix} 〈 〈ψ | ψ〉 〉 = \sum_{i} ψ_{i}^{‡} ψ_{i} = \sum_{i} ρ_{i} R_{i}^{‡} R_{i} = \sum_{i} ρ_{i} = Z \end{matrix}

(71)

□

Definition 8.

(Spin(2)-valued Evolution Operator).

\begin{matrix} T = [\begin{matrix} e^{- \frac{1}{2} θ b_{1}} \\ ⋱ \\ e^{- \frac{1}{2} θ b_{n}} \end{matrix}] \end{matrix}

(72)

Theorem 7.

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

Proof.

We note that:

\begin{matrix} 〈 〈T v | T v〉 〉 = 〈 〈v | v〉 〉 = v^{‡} T^{‡} T v \Rightarrow T^{‡} T = I \end{matrix}

(73)

then, since

[\begin{matrix} e^{\frac{1}{2} θ b_{1}} \\ ⋱ \\ e^{\frac{1}{2} θ b_{n}} \end{matrix}] [\begin{matrix} e^{- \frac{1}{2} θ b_{1}} \\ ⋱ \\ e^{- \frac{1}{2} θ b_{n}} \end{matrix}] = I

, the relation

T^{‡} T = I

is satisfied. □

We note that the even sub-algebra of

GA (2)

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

θ

Definition 9.

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

\begin{matrix} \frac{d}{d θ} [\begin{matrix} ψ_{1} (θ) \\ ⋮ \\ ψ_{n} (θ) \end{matrix}] = [\begin{matrix} - \frac{1}{2} b_{1} \\ ⋱ \\ - \frac{1}{2} b_{n} \end{matrix}] [\begin{matrix} ψ_{1} (θ) \\ ⋮ \\ ψ_{n} (θ) \end{matrix}] \end{matrix}

(74)

Here,

θ

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

ψ = \sqrt{ρ} R e^{i b / 2}

, where

R = e^{f / 2}

is a Lorentz boost or rotation and where

e^{i b / 2}

is a phase. In 2D, as the algebra only admits a bivector, his formulation would reduce to

ψ = \sqrt{ρ} 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

σ_{x}

and

σ_{y}

, the Dirac current for the 2D theory is defined as:

\begin{matrix} J_{1} \equiv ψ^{‡} σ_{x} ψ = ρ \underset{SO (2)}{\underset{︸}{R^{‡} σ_{x} R}} = ρ {\tilde{σ}}_{x} \end{matrix}

(75)

\begin{matrix} J_{2} \equiv ψ^{‡} σ_{y} ψ = ρ \underset{SO (2)}{\underset{︸}{R^{‡} σ_{y} R}} = ρ {\tilde{σ}}_{y} \end{matrix}

(76)

where

{\tilde{σ}}_{x}

and

{\tilde{σ}}_{y}

are a SO(2) rotated basis vectors.

Obstructions

As stated, all dimensional 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 e_{0} \land e_{1}) (a + b e_{0} \land e_{1})$ , which yields negative probabilities: $a^{2} - b^{2} \in R$ for certain wavefunction states, in contrast to the non-negative probabilities $a^{2} + b^{2} \in R^{\geq 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 ( $σ_{x}$ and $σ_{y}$ in 2D, and $e_{0}$ and $e_{1}$ in 1+1D) are not self-adjoint. Although useable in the context 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{matrix} {(e_{μ} u)}^{‡} u \neq u^{‡} e_{μ} u \end{matrix}$

(77)

because ${(e_{μ} u)}^{‡} u = u^{‡} e_{μ}^{‡} u = u^{‡} (- e_{μ}) 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{matrix} L = \underset{Relative Shannon Entropy}{\underset{︸}{- \sum_{i} ρ_{i} ln \frac{ρ_{i}}{p_{i}}}} + \underset{Normalization Constraint}{\underset{︸}{λ (1 - \sum_{i} ρ_{i})}} + \underset{Vanishing {Spin}^{c} (3, 1) - Phase}{\underset{︸}{ζ (- tr \frac{1}{2} \sum_{i} ρ_{i} M_{i})}} \end{matrix}

(78)

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

\begin{matrix} ρ_{i} = \underset{Spin ​^{c} (3, 1) Invariant Ensemble}{\underset{︸}{\frac{1}{\sum_{i} p_{i} det exp (- ζ \frac{1}{2} M_{i})}}} \underset{{Spin}^{c} (3, 1) Born Rule}{\underset{︸}{det exp (- ζ \frac{1}{2} M_{i})}} \underset{Initial Preparation}{\underset{︸}{p_{i}}} \end{matrix}

(79)

where

ζ

is a "twisted-phase" rapidity. (If the invariance group was Spin(3,1) instead of Spin

​^{c}

(3,1), obtainable by posing

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

GA (3, 1)

Definition 12

(Multivector). Let

u

be a multivector of

GA (3, 1)

. Its general form is:

\begin{matrix} u & = a \end{matrix}

(80)

\begin{matrix} + t γ_{0} + x γ_{1} + y γ_{2} + z γ_{3} \end{matrix}

(81)

\begin{matrix} + f_{01} γ_{0} \land γ_{1} + f_{02} γ_{0} \land γ_{2} + f_{03} γ_{0} \land γ_{3} + f_{12} γ_{1} \land γ_{2} + f_{13} γ_{1} \land γ_{3} + f_{23} γ_{2} \land γ_{3} \end{matrix}

(82)

\begin{matrix} + p γ_{1} \land γ_{2} \land γ_{3} + q γ_{0} \land γ_{2} \land γ_{3} + v γ_{0} \land γ_{1} \land γ_{3} + w γ_{0} \land γ_{1} \land γ_{2} \end{matrix}

(83)

\begin{matrix} + b γ_{0} \land γ_{1} \land γ_{2} \land γ_{3} \end{matrix}

(84)

where

γ_{0}, γ_{1}, γ_{2}, γ_{3}

are the basis vectors in the real Majorana representation.

A more compact notation for

u

\begin{matrix} u = a + x + f + v + b \end{matrix}

(85)

where a is a scalar,

x

a vector,

f

a bivector,

v

is pseudo-vector and

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

u

\begin{matrix} M = [\begin{matrix} 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{matrix}] \end{matrix}

(86)

To manipulate and analyze multivectors in

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{matrix} u^{‡} = a - x - f + v + b \end{matrix}

(87)

Definition 15

(3,4 Blade Conjugate). The 3,4 blade conjugate of

u

\begin{matrix} {⌊ u ⌋}_{3, 4} = a + x + f - v - b \end{matrix}

(88)

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 (u v) = N (u) N (v)

to the domain of multivectors:

Definition 16.

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

\begin{matrix} GA (0, 1) : & φ^{†} φ \end{matrix}

(89)

\begin{matrix} GA (2, 0) : & φ^{‡} φ \end{matrix}

(90)

\begin{matrix} GA (3, 0) : & {⌊ φ^{‡} φ ⌋}_{3} φ^{‡} φ \end{matrix}

(91)

\begin{matrix} GA (3, 1) : & {⌊ φ^{‡} φ ⌋}_{3, 4} φ^{‡} φ \end{matrix}

(92)

\begin{matrix} GA (4, 1) : & {({⌊ φ^{‡} φ ⌋}_{3, 4} φ^{‡} φ)}^{†} ({⌊ φ^{‡} φ ⌋}_{3, 4} φ^{‡} φ) \end{matrix}

(93)

We can now express the determinant of the matrix representation of a multivector via the self-product

{⌊ φ^{‡} φ ⌋}_{3, 4} φ^{‡} φ

. 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

GA (3, 1)

Theorem 8

(Determinant as a Multivector Self-Product).

\begin{matrix} {⌊ u^{‡} u ⌋}_{3, 4} u^{‡} u = det M \end{matrix}

(94)

Proof.

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

Definition 17

(

GA (3, 1)

-valued Vector).

\begin{matrix} |V〉 〉 = [\begin{matrix} u_{1} \\ ⋮ \\ u_{n} \end{matrix}] = [\begin{matrix} a_{1} + x_{1} + f_{1} + v_{1} + b_{1} \\ ⋮ \\ a_{n} + x_{n} + f_{n} + v_{n} + b_{n} \end{matrix}] \end{matrix}

(95)

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{matrix} 〈 〈V | V | V | V〉 〉 & = \sum_{i} {⌊ \underset{copy}{\underset{︸}{ψ_{i}^{‡} ψ_{i}}} 1 ⌋}_{3, 4} \underset{copy}{\underset{︸}{ψ_{i}^{‡} ψ_{i}}} 2 \end{matrix}

(96)

\begin{matrix} = {⌊ \underset{copy}{\underset{︸}{[\begin{matrix} u_{1}^{‡} & \dots & u_{n} \end{matrix}] [\begin{matrix} u_{1} & \dots & 0 \\ ⋮ & ⋱ & ⋮ \\ 0 & \dots & u_{n} \end{matrix}]}} 1 ⌋}_{3, 4} \underset{copy}{\underset{︸}{[\begin{matrix} u_{1}^{‡} & \dots & 0 \\ ⋮ & ⋱ & ⋮ \\ 0 & \dots & u_{n}^{‡} \end{matrix}] [\begin{matrix} u_{1} \\ ⋮ \\ u_{n} \end{matrix}]}} 2 \end{matrix}

(97)

Theorem 9

(Partition Function).

Z = 〈 〈V | V | V | V〉 〉

Proof.

\begin{matrix} 〈 〈V | V | V | V〉 〉 \end{matrix}

(98)

\begin{matrix} = {⌊ [\begin{matrix} u_{1}^{‡} & \dots & u_{n} \end{matrix}] [\begin{matrix} u_{1} & \dots & 0 \\ ⋮ & ⋱ & ⋮ \\ 0 & \dots & u_{n} \end{matrix}] ⌋}_{3, 4} [\begin{matrix} u_{1}^{‡} & \dots & 0 \\ ⋮ & ⋱ & ⋮ \\ 0 & \dots & u_{n}^{‡} \end{matrix}] [\begin{matrix} u_{1} \\ ⋮ \\ u_{n} \end{matrix}] \end{matrix}

(99)

\begin{matrix} = {⌊ [\begin{matrix} u_{1}^{‡} u_{1} & \dots & u_{n} u_{n} \end{matrix}] ⌋}_{3, 4} [\begin{matrix} u_{1}^{‡} u_{1} \\ ⋮ \\ u_{n}^{‡} u_{n} \end{matrix}] \end{matrix}

(100)

\begin{matrix} = {⌊ u_{1}^{‡} u_{1} ⌋}_{3, 4} u_{1}^{‡} u_{1} + \dots + {⌊ u_{n}^{‡} u_{n} ⌋}_{3, 4} u_{n}^{‡} u_{n} \end{matrix}

(101)

\begin{matrix} = \sum_{i = 1}^{n} det M_{u_{i}} \end{matrix}

(102)

\begin{matrix} = Z \end{matrix}

(103)

□

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

GA (3, 1)

is always non-negative.

Proof.

Let

|V〉 〉 = [\begin{matrix} a_{1} + f_{1} + b_{1} \\ ⋮ \\ a_{n} + f_{n} + b_{n} \end{matrix}]

. Then,

\begin{matrix} 〈 〈V | V | V | V〉 〉 \end{matrix}

(104)

\begin{matrix} = {⌊ [\begin{matrix} {(a_{1} + f_{1} + b_{1})}^{‡} (a_{1} + f_{1} + b_{1}) & \dots \end{matrix}] ⌋}_{3, 4} [\begin{matrix} {(a_{1} + f_{1} + b_{1})}^{‡} (a_{1} + f_{1} + b_{1}) \\ ⋮ \end{matrix}] \end{matrix}

(105)

\begin{matrix} = {⌊ [\begin{matrix} (a_{1} - f_{1} + b_{1}) (a_{1} + f_{1} + b_{1}) & \dots \end{matrix}] ⌋}_{3, 4} [\begin{matrix} (a_{1} - f_{1} + b_{1}) (a_{1} + f_{1} + b_{1}) \\ ⋮ \end{matrix}] \end{matrix}

(106)

\begin{matrix} = {⌊ [\begin{matrix} a_{1}^{2} + a_{1} f_{1} + a_{1} b_{1} - f_{1} a_{1} - f_{1}^{2} - f_{1} b_{1} + b_{1} a_{1} + b_{1} f_{1} + b_{1}^{2} & \dots \end{matrix}] ⌋}_{3, 4} \dots \end{matrix}

(107)

\begin{matrix} = {⌊ [\begin{matrix} a_{1}^{2} - f_{1}^{2} + b_{1}^{2} & \dots \end{matrix}] ⌋}_{3, 4} \dots \end{matrix}

(108)

We note 1)

b^{2} = {(b I)}^{2} = - b^{2}

and 2)

f^{2} = - E_{1}^{2} - E_{2}^{2} - E_{3}^{2} + B_{1}^{2} + B_{2}^{2} + B_{3}^{2} + 4 e_{0} e_{1} e_{2} e_{3} (E_{1} B_{1} + E_{2} B_{2} + E_{3} B_{3})

\begin{matrix} = {⌊ [\begin{matrix} a_{1}^{2} - b_{1}^{2} + E_{1}^{2} + E_{2}^{2} + E_{3}^{2} - B_{1}^{2} - B_{2}^{2} - B_{3}^{2} - 4 e_{0} e_{1} e_{2} e_{3} (E_{1} B_{1} + E_{2} B_{2} + E_{3} B_{3}) & \dots \end{matrix}] ⌋}_{3, 4} \dots \end{matrix}

(109)

We note that the terms are now complex numbers, which we rewrite as

ℜ (z) = 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

ℑ (z) = - 4 (E_{1} B_{1} + E_{2} B_{2} + E_{3} B_{3})

\begin{matrix} = {⌊ [\begin{matrix} z_{1} & \dots & z_{2} \end{matrix}] ⌋}_{3, 4} [\begin{matrix} z_{n} \\ ⋮ \\ z_{n} \end{matrix}] \end{matrix}

(110)

\begin{matrix} = [\begin{matrix} z_{1}^{†} & \dots & z_{2}^{†} \end{matrix}] [\begin{matrix} z_{n} \\ ⋮ \\ z_{n} \end{matrix}] \end{matrix}

(111)

\begin{matrix} = z_{1}^{‡} z_{1} + \dots + z_{n}^{‡} z_{n} \end{matrix}

(112)

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

{Spin}^{c} (3, 1)

-valued wavefunction, which is valued in the even sub-algebra of

GA (3, 1)

, as follows:

Definition 19

(

{Spin}^{c} (3, 1)

-valued Wavefunction).

\begin{matrix} |ψ〉 〉 = [\begin{matrix} e^{\frac{1}{2} (a_{1} + f_{1} + b_{1})} \\ ⋮ \\ e^{\frac{1}{2} (a_{n} + f_{n} + b_{n})} \end{matrix}] = [\begin{matrix} \sqrt{ρ_{1}} R_{1} B_{1} \\ ⋮ \\ \sqrt{ρ_{n}} R_{n} B_{n} \end{matrix}] \end{matrix}

(113)

where

R_{i} = e^{\frac{1}{2} f_{i}}

is a rotor,

B_{i} = e^{\frac{1}{2} b_{i}}

is a phase, and where

\sqrt{ρ_{i}} = e^{\frac{1}{2} a_{i}} \geq 0

Any even multivectors of GA(3,1) admits a unique exponential representation, except when

ρ_{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

(

{Spin}^{c} (3, 1)

Evolution Operator).

\begin{matrix} T = [\begin{matrix} e^{- \frac{1}{2} ζ (f_{1} + b_{1})} \\ ⋱ \\ e^{- \frac{1}{2} ζ (f_{n} + b_{n})} \end{matrix}] \end{matrix}

(114)

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

ζ

Definition 21

(

{Spin}^{c} (3, 1)

-valued Schrödinger equation).

\begin{matrix} \frac{d}{d ζ} [\begin{matrix} ψ_{1} (ζ) \\ ⋮ \\ ψ_{n} (ζ) \end{matrix}] = [\begin{matrix} - \frac{1}{2} (f_{1} + b_{1}) \\ ⋱ \\ - \frac{1}{2} (f_{n} + b_{n}) \end{matrix}] [\begin{matrix} ψ_{1} (ζ) \\ ⋮ \\ ψ_{n} (ζ) \end{matrix}] \end{matrix}

(115)

In this case

ζ

represents a one-parameter evolution parameter akin to time, which is able to transform the measurement basis under action of the

{Spin}^{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.

2.3.2. RQM

Definition 22

(David Hestenes’ Formulation). The

{Spin}^{c} (3, 1)

-valued wavefunction we have recovered is formulated identically to David Hestenes’ [7] formulation of the wavefunction within GA(3,1).

\begin{matrix} ψ = e^{\frac{1}{2} (a + f + b)} = \sqrt{ρ} R e^{- i b / 2} \end{matrix}

(116)

where

e^{\frac{1}{2} a} = \sqrt{ρ}

e^{\frac{1}{2} f} = R

and

e^{\frac{1}{2} b} = e^{- i b / 2}

Before we continue the RQM investigation, let us note that the double-copy product contains two copies of a bilinear form

ψ^{‡} ψ

\begin{matrix} {⌊ \underset{copy 1}{\underset{︸}{ψ^{‡} ψ}} ⌋}_{3, 4} \underset{copy 2}{\underset{︸}{ψ^{‡} ψ}} \end{matrix}

(117)

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{matrix} J & \equiv ψ^{‡} γ_{μ} ψ = ρ R^{‡} B^{‡} γ_{μ} B R = ρ R^{‡} γ_{μ} B^{- 1} B R = ρ \underset{SO (3, 1)}{\underset{︸}{R^{‡} γ_{μ} R}} = ρ {\tilde{γ}}_{μ} \end{matrix}

(118)

where

{\tilde{γ}}_{μ}

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

U (1)

SU (2)

, and

SU (3)

symmetries and the

Spin (3, 1)

and unitary

U^{†} U = I

symmetries which play a fundamental role in the standard model of particle physics.

Theorem 11.

(U(1) Invariance). Let

e^{\frac{1}{2} b}

be a general element of U(1). Then, the equality

\begin{matrix} {⌊ ψ^{‡} γ_{0} ψ ⌋}_{3, 4} ψ^{‡} γ_{0} ψ = {⌊ \underset{copy 1}{\underset{︸}{{(e^{\frac{1}{2} b} ψ)}^{‡} γ_{0} e^{\frac{1}{2} b} ψ}} ⌋}_{3, 4} \underset{copy 2}{\underset{︸}{{(e^{\frac{1}{2} b} ψ)}^{‡} γ_{0} e^{\frac{1}{2} b} ψ}} \end{matrix}

(119)

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

Proof.

Equation 119 is invariant if this expression is satisfied:

\begin{matrix} e^{\frac{1}{2} b} γ_{0} e^{\frac{1}{2} b} = γ_{0} \end{matrix}

(120)

This is always satisfied simply because

e^{\frac{1}{2} b} γ_{0} e^{\frac{1}{2} b} = γ_{0} e^{- \frac{1}{2} b} e^{\frac{1}{2} b} = γ_{0}

□

Theorem 12.

(SU(2) Invariance). Let

e^{\frac{1}{2} f}

be a general element of Spin(3,1). Then, the equality:

\begin{matrix} {⌊ ψ^{‡} γ_{0} ψ ⌋}_{3, 4} ψ^{‡} γ_{0} ψ = {⌊ \underset{copy 1}{\underset{︸}{{(e^{\frac{1}{2} f} ψ)}^{‡} γ_{0} e^{\frac{1}{2} f} ψ}} ⌋}_{3, 4} \underset{copy 2}{\underset{︸}{{(e^{\frac{1}{2} f} ψ)}^{‡} γ_{0} e^{\frac{1}{2} f} ψ}} \end{matrix}

(121)

is satisfied for if

f = θ_{1} γ_{0} γ_{1} + θ_{2} γ_{0} γ_{2} + θ_{3} γ_{0} γ_{3}

(which generates SU(2)), yielding a SU(2) symmetry for each copied bilinear form.

Proof.

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

\begin{matrix} e^{- \frac{1}{2} f} γ_{0} e^{\frac{1}{2} f} = γ_{0} \end{matrix}

(122)

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

\begin{matrix} e^{- θ_{1} γ_{0} γ_{1} - θ_{2} γ_{0} γ_{2} - θ_{3} γ_{0} γ_{3} - B_{1} γ_{2} γ_{3} - B_{2} γ_{1} γ_{3} - B_{3} γ_{1} γ_{2}} γ_{0} e^{\frac{1}{2} f} \end{matrix}

(123)

\begin{matrix} = γ_{0} e^{- θ_{1} γ_{0} γ_{1} - θ_{2} γ_{0} γ_{2} - θ_{3} γ_{0} γ_{3} + B_{1} γ_{2} γ_{3} + B_{2} γ_{1} γ_{3} + B_{3} γ_{1} γ_{2}} e^{\frac{1}{2} f} \end{matrix}

(124)

Therefore, the product

e^{- \frac{1}{2} f} γ_{0} e^{\frac{1}{2} f}

reduces to

γ_{0}

if and only if

B_{1} = B_{2} = B_{3} = 0

, leaving

f = θ_{1} γ_{0} γ_{1} + θ_{2} γ_{0} γ_{2} + θ_{3} γ_{0} γ_{3}

Finally, we note that

e^{θ_{1} γ_{0} γ_{1} + θ_{2} γ_{0} γ_{2} + θ_{3} γ_{0} γ_{3}}

generates

SU (2)

. □

Theorem 13.

(SU(3) invariance). Let

f = E_{1} γ_{0} γ_{1} + E_{2} γ_{0} γ_{2} + E_{3} γ_{0} γ_{3} + B_{1} γ_{2} γ_{3} + B_{2} γ_{1} γ_{3} + B_{3} γ_{1} γ_{2}

be a general bivector. Then, the equality:

\begin{matrix} {⌊ ψ^{‡} γ_{0} ψ ⌋}_{3, 4} ψ^{‡} γ_{0} ψ = {⌊ \underset{copy 1}{\underset{︸}{{(f ψ)}^{‡} γ_{0} f ψ}} ⌋}_{3, 4} \underset{copy 2}{\underset{︸}{{(f ψ)}^{‡} γ_{0} f ψ}} \end{matrix}

(125)

is satisfied if

(E_{1}^{2} + B_{1}^{2}) + (E_{2}^{2} + B_{2}^{2}) + (E_{3}^{2} + B_{3}^{2}) = 1

, where

E_{1} γ_{0} γ_{1} + B_{1} γ_{0} γ_{1} I ≅ E_{1} + i B_{1}

E_{2} γ_{0} γ_{2} + B_{2} γ_{0} γ_{2} I ≅ E_{2} + i B_{2}

and

E_{3} γ_{0} γ_{3} + B_{3} γ_{0} γ_{3} I ≅ E_{3} + i B_{3}

which are the defining conditions for the SU(3) symmetry group. This yields a SU(3) symmetry for each copy of the bilinear form.

Proof.

Equation 125 is invariant if this expression is satisfied [9,10]:

\begin{matrix} - f γ_{0} f = γ_{0} \end{matrix}

(126)

which we can rewrite as follows:

\begin{matrix} - (E_{1} γ_{0} γ_{1} + E_{2} γ_{0} γ_{2} + E_{3} γ_{0} γ_{3} + B_{1} γ_{2} γ_{3} + B_{2} γ_{1} γ_{3} + B_{3} γ_{1} γ_{2}) γ_{0} f \end{matrix}

(127)

The first three terms anticommute with

γ_{0}

, while the last three commute with

γ_{0}

\begin{matrix} = γ_{0} (E_{1} γ_{0} γ_{1} + E_{2} γ_{0} γ_{2} + E_{3} γ_{0} γ_{3} - B_{1} γ_{2} γ_{3} - B_{2} (q) γ_{1} γ_{3} - B_{3} (q) γ_{1} γ_{2}) f (q) \end{matrix}

(128)

This can be written as:

\begin{matrix} γ_{0} (E - B) (E + B) \end{matrix}

(129)

\begin{matrix} = γ_{0} (E^{2} + E B - B E - B^{2}) \end{matrix}

(130)

where

E = E_{1} γ_{0} γ_{1} + E_{2} γ_{0} γ_{2} + E_{3} γ_{0} γ_{3}

and

B = B_{1} γ_{2} γ_{3} + B_{2} γ_{1} γ_{3} + B_{3} γ_{1} γ_{2}

Thus, for

- f γ_{0} f = γ_{0}

, we require: 1)

E^{2} - B^{2} = 1

and 2)

E B = B E

. The second requirement means that

E

and

B

must commute (and thus be isomorphic to three complex numbers), and the first implies:

\begin{matrix} E^{2} - B^{2} = (E_{1}^{2} + B_{1}^{2}) + (E_{2}^{2} + B_{2}^{2}) + (E_{3}^{2} + B_{3}^{2}) = 1 \end{matrix}

(131)

which are the defining conditions for the

SU (3)

symmetry group. □

Theorem 14.

(Spin(3,1) invariance). Let

e^{\frac{1}{2} f}

be a general element of Spin(3,1). Then, the equality:

\begin{matrix} {⌊ ψ^{‡} ψ ⌋}_{3, 4} ψ^{‡} ψ = {⌊ \underset{copy 1}{\underset{︸}{{(e^{\frac{1}{2} f} ψ)}^{‡} e^{\frac{1}{2} f} ψ}} ⌋}_{3, 4} \underset{copy 2}{\underset{︸}{{(e^{\frac{1}{2} f} ψ)}^{‡} e^{\frac{1}{2} f} ψ}} \end{matrix}

(132)

is always satisfied, yielding a Spin(3,1) symmetry for each copied bilinear form.

Proof.

Equation 132 is invariant if this expression is satisfied

\begin{matrix} e^{- \frac{1}{2} f} e^{\frac{1}{2} f} = 1 \end{matrix}

(133)

which is always the case. □

Theorem 15

(Unitary invariance). Let U be

n \times n

unitary matrices. Then unitary invariance:

\begin{matrix} 〈 〈ψ | γ_{μ} ψ | ψ | γ_{ν} ψ〉 〉 = 〈 〈U ψ | γ_{μ} U ψ | U ψ | γ_{ν} U ψ〉 〉 \Rightarrow U^{†} U = I \end{matrix}

(134)

is individually satisfied for each copied bilinear form.

Proof.

Equation 134 is satisfied if

U^{‡} γ_{μ} U = γ_{μ}

. Since U is valued in complex numbers, then

U^{‡} = U^{T}

, and since

γ_{μ} γ_{0} γ_{1} γ_{2} γ_{3} = - γ_{0} γ_{1} γ_{2} γ_{3} γ_{μ}

, it follows that:

\begin{matrix} γ_{μ} U^{†} U = γ_{μ} \end{matrix}

(135)

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{matrix} {(e^{- i θ (x)} ψ (x))}^{†} e^{- i θ (x)} ψ (x) = ψ {(x)}^{†} ψ (x) \end{matrix}

(136)

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. This suggests that the double-copy product framework may provide a more unified and consistent approach to deriving the fundamental symmetries of nature.

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 and the standard model gauge symmetries. The goal of this section is to extend the methodology to arbitrary basis vectors, in which the metric tensor exists as an observable. To formulate the theory, we will exploit the features of the double-copy product, which will enable us to express the metric tensor as an observable. To a degree, this formulation is reminiscent of the Bern-Carrasco-Johansson (BCJ) [11] double-copy gauge theory of perturbatively expanded quantum gravity, but our double-copy is formulated at the level of the Dirac current directly. By applying the double-copy product to the Dirac current, we aim to establish a connection between the wavefunction and the geometrical structure of spacetime, providing a novel approach to quantize gravity.

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

\begin{matrix} g_{μ ν} = \frac{1}{2} (e_{μ} e_{ν} + e_{ν} e_{μ}) \end{matrix}

(137)

Then, we note that the double-copy product acts on a pair of basis element

e_{μ}

and

e_{ν}

, as follows:

\begin{matrix} \frac{1}{2} ({⌊ \underset{copy 1}{\underset{︸}{ψ^{‡} e_{μ} ψ}} ⌋}_{3, 4} \underset{copy 1}{\underset{︸}{ψ^{‡} e_{ν} ψ}} + {⌊ \underset{copy 1}{\underset{︸}{ψ^{‡} e_{ν} φ}} ⌋}_{3, 4} \underset{copy 1}{\underset{︸}{ψ^{‡} e_{μ} ψ}}) \end{matrix}

(138)

\begin{matrix} = \frac{1}{2} (\tilde{R} \underset{Born rule copy 1}{\underset{︸}{ρ e^{i b / 2} e^{- i b / 2}}} e_{μ} R \tilde{R} \underset{Born rule copy 2}{\underset{︸}{ρ e^{- i b / 2} e^{i b / 2}}} e_{ν} R + \tilde{R} \underset{Born rule copy 2}{\underset{︸}{ρ e^{i b / 2} e^{- i b / 2}}} e_{ν} R \tilde{R} \underset{Born rule copy 1}{\underset{︸}{ρ e^{- i b / 2} e^{i b / 2}}} e_{μ} R) \end{matrix}

(139)

\begin{matrix} = \frac{1}{2} ρ^{2} (\tilde{R} e_{μ} R \tilde{R} e_{ν} R + \tilde{R} e_{ν} R \tilde{R} e_{ν} R) \end{matrix}

(140)

\begin{matrix} = ρ^{2} \underset{metric tensor}{\underset{︸}{\frac{1}{2} ({\tilde{e}}_{μ} {\tilde{e}}_{ν} + {\tilde{e}}_{ν} {\tilde{e}}_{μ})}} \end{matrix}

(141)

where

{\tilde{e}}_{μ}

and

{\tilde{e}}_{ν}

are SO(3,1) rotated basis vectors.

As one can swap

e_{μ}

and

e_{ν}

and obtain the same metric tensor, the double-copy product guarantees that

g_{μ ν}

is symmetric.

Furthermore, since

e_{μ}^{‡} = - e_{μ}

, we get:

\begin{matrix} {⌊ {(e_{μ} ψ)}^{‡} ψ ⌋}_{3, 4} {(e_{ν} ψ)}^{‡} ψ \end{matrix}

(142)

\begin{matrix} = {⌊ ψ^{‡} (- 1) e_{μ}^{‡} ψ ⌋}_{3, 4} ψ^{‡} (- 1) e_{ν}^{‡} ψ \end{matrix}

(143)

\begin{matrix} = {⌊ ψ^{‡} e_{μ} ψ ⌋}_{3, 4} ψ^{‡} e_{ν} ψ \end{matrix}

(144)

which allows us to conclude that

e_{μ}

and

e_{ν}

are self-adjoint within the double-copy product, entailing the interpretation of

g_{μ ν}

as an observable.

In the double-copy product, the metric tensor is a double copy of Dirac currents (we recall that the Dirac current is defined as

ρ {\tilde{γ}}_{μ} = ψ^{‡} γ_{μ} ψ

). As such, we intuit that it encodes the probabilistic structure of a general relativistic quantum theory in the form of a metric tensor analogously to how the Dirac current encodes the probabilistic structure of a special relativistic quantum theory in the form of a vector.

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

\begin{matrix} T = [\begin{matrix} e^{- \frac{1}{2} ζ (f_{1} + b_{1})} \\ ⋱ \\ e^{- \frac{1}{2} ζ (f_{n} + b_{n})} \end{matrix}] \end{matrix}

(145)

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

\begin{matrix} {⌊ \underset{copy 1}{\underset{︸}{ψ^{‡} T^{‡} e_{μ} T ψ}} ⌋}_{3, 4} \underset{copy 2}{\underset{︸}{ψ^{‡} T^{‡} e_{ν} T ψ}} \end{matrix}

(146)

which realizes an

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{matrix} {⌊ \underset{copy 1}{\underset{︸}{ψ^{‡} \underset{SO (3, 1) evolution}{\underset{︸}{e^{\frac{1}{2} ζ f} e_{μ} e^{- \frac{1}{2} ζ f}}} \underset{unitary evolution}{\underset{︸}{e^{\frac{1}{2} ζ b} e^{- \frac{1}{2} ζ b}}} ψ}} ⌋}_{3, 4} \underset{copy 2}{\underset{︸}{ψ^{‡} \underset{SO (3, 1) evolution}{\underset{︸}{e^{\frac{1}{2} ζ f} e_{μ} e^{- \frac{1}{2} ζ f}}} \underset{unitary evolution}{\underset{︸}{e^{\frac{1}{2} ζ b} e^{- \frac{1}{2} ζ b}}} ψ}} \end{matrix}

(147)

In summary, this initial investigation has identified a scenario in which the wavefunction measures the metric tensor by acting on the basis vectors. In this context, 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 given by a double-copy of unitary quantum theories operating on the basis vectors. This formulation respects both the SO(3,1) symmetry, which is essential for describing the structure of spacetime, and the unitary symmetry, which is a fundamental aspect of quantum mechanics. We will now investigate this formulation in the context of gravity.

2.3.5. Linearized Gravity (Sketch)

It is well known that the Einstein-Hilbert action:

\begin{matrix} S [g_{μ ν}] = \frac{c^{4}}{16 π G} \int R \sqrt{- | g |} d^{4} x \end{matrix}

(148)

under the assumption of a small perturbation

g_{μ ν} = η_{μ ν} + h_{μ ν}

, and working in de Donder gauge

\partial^{α} h_{α μ} - \frac{1}{2} \partial_{μ} h = 0

, where

h = η^{μ ν} h_{μ ν}

, can be reduced to its linearized form which is:

\begin{matrix} S_{E H}^{(1)} [h_{μ ν}] = \int d^{4} x (\frac{1}{2} \partial_{μ} h_{ρ σ} \partial^{μ} h^{ρ σ} - \frac{1}{4} \partial_{μ} h \partial^{μ} h) \end{matrix}

(149)

Furthermore, varying this action with respect to

h_{μ ν}

and applying the transverse-traceless gauge, yields the wave equation as the equation of motion:

\begin{matrix} \frac{δ S}{δ h_{μ ν}} = 0 \Rightarrow □ h_{μ ν} = 0 \end{matrix}

(150)

We now wish to express the wave equation for

h_{μ ν}

in terms of basis vectors

h_{μ}

and

h_{ν}

, such that

h_{μ ν} = \frac{1}{2} (h_{μ} h_{ν} + h_{ν} h_{μ})

. The expression of the wave equation becomes:

\begin{matrix} □ \frac{1}{2} (h_{μ} h_{ν} + h_{ν} h_{μ}) = 0 \end{matrix}

(151)

We identify the solution for

h_{μ}

and

h_{ν}

by an ansatz:

\begin{matrix} h_{μ} (\vec{x}, t) = Re \int d^{3} k (A_{+} (\vec{k}) e_{μ}^{+} (\vec{k}) + A_{\times} (\vec{k}) e_{μ}^{\times} (\vec{k})) e^{- i (\vec{k} \cdot \vec{x} - ω_{k} t)} \end{matrix}

(152)

\begin{matrix} h_{ν} (\vec{x}, t) = Re \int d^{3} k (A_{+} (\vec{k}) e_{ν}^{+} (\vec{k}) + A_{\times} (\vec{k}) e_{ν}^{\times} (\vec{k})) e^{- i (\vec{k} \cdot \vec{x} - ω_{k} t)} \end{matrix}

(153)

The × and + symbol designed two polarizations.

Finally, we promote

h_{μ} (\vec{x}, t)

and

h_{ν} (\vec{x}, t)

to operators:

\begin{matrix} {\hat{h}}_{μ} (\vec{x}, t) = Re \int d^{3} k ({\hat{A}}_{+} (\vec{k}) e_{μ}^{+} (\vec{k}) + {\hat{A}}_{\times} (\vec{k}) e_{μ}^{\times} (\vec{k})) e^{- i (\vec{k} \cdot \vec{x} - ω_{k} t)} \end{matrix}

(154)

\begin{matrix} {\hat{h}}_{ν} (\vec{x}, t) = Re \int d^{3} k ({\hat{A}}_{+} (\vec{k}) e_{ν}^{+} (\vec{k}) + {\hat{A}}_{\times} (\vec{k}) e_{ν}^{\times} (\vec{k})) e^{- i (\vec{k} \cdot \vec{x} - ω_{k} t)} \end{matrix}

(155)

In the conventional perturbative approach to quantum gravity, the metric tensor

{\hat{h}}_{μ ν}

is quantized, and its expectation value is calculated using the Born rule as follows:

〈 {\hat{h}}_{μ ν} 〉 = 〈 ψ | {\hat{h}}_{μ ν} | ψ 〉

. In contrast, our approach quantizes the two basis vectors

{\hat{h}}_{μ}

and

{\hat{h}}_{ν}

separately and utilizes a double-copy product to obtain the expectation value of the metric tensor as

〈 {\hat{h}}_{μ ν} 〉 = \frac{1}{2} (〈 〈 ψ | {\hat{h}}_{μ} ψ | ψ | {\hat{h}}_{ν} ψ 〉 〉 + 〈 〈 ψ | {\hat{h}}_{ν} ψ | ψ | {\hat{h}}_{μ} ψ 〉 〉)

. In the latter case, the operators realize a double copy of Dirac currents, each copy requiring an individual application of the Born rule. The resulting probability calculation differs from the conventional approach, as it involves the product of two separate Born rule applications. Remarkably, our probability measure appears unable to assign an expectation value to a quantized metric tensor unless it is done via a double-copy mechanism, making the double-copy mechanism fundamental to our construction.

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

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 dimensional configurations except those we have explored here (e.g.

GA (0) ≅ R

GA (0, 1) ≅ C

and

GA (3, 1)

) are obstructed:

\begin{matrix} Dimensions & Obstruction \end{matrix}

\begin{matrix} GA (0) ≅ R & Unobstructed \Rightarrow statistical mechanics \end{matrix}

(156)

\begin{matrix} GA (0, 1) ≅ C & Unobstructed \Rightarrow quantum mechanics \end{matrix}

(157)

\begin{matrix} GA (1, 0) & Negative probabilities in the RQM \end{matrix}

(158)

\begin{matrix} GA (2, 0) & No metric measurement \end{matrix}

(159)

\begin{matrix} GA (1, 1) & Negative probabilities in the RQM \end{matrix}

(160)

\begin{matrix} GA (0, 2) & Not isomorphic to a real matrix algebra \end{matrix}

(161)

\begin{matrix} GA (3, 0) & Not isomorphic to a real matrix algebra \end{matrix}

(162)

\begin{matrix} GA (2, 1) & Not isomorphic to a real matrix algebra \end{matrix}

(163)

\begin{matrix} GA (1, 2) & Not isomorphic to a real matrix algebra \end{matrix}

(164)

\begin{matrix} GA (0, 3) & Not isomorphic to a real matrix algebra \end{matrix}

(165)

\begin{matrix} GA (4, 0) & Not isomorphic to a real matrix algebra \end{matrix}

(166)

\begin{matrix} GA (3, 1) & Unobstructed \Rightarrow quantum gravity \land SU (3) \times SU (2) \times U (1) \end{matrix}

(167)

\begin{matrix} GA (2, 2) & Negative probabilities in the RQM \end{matrix}

(168)

\begin{matrix} GA (1, 3) & Not isomorphic to a real matrix algebra \end{matrix}

(169)

\begin{matrix} GA (0, 4) & Not isomorphic to a real matrix algebra \end{matrix}

(170)

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

(171)

\begin{matrix} GA (6, 0) & No probability measure as a self - product \\ ⋮ & ⋮ \end{matrix}

(172)

\begin{matrix} \infty \end{matrix}

(173)

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{matrix} GA (0, 2) ≅ H \end{matrix}

(174)

\begin{matrix} GA (3, 0) ≅ M_{2} (C) \end{matrix}

(175)

\begin{matrix} GA (2, 1) ≅ M_{2}^{2} (R) \end{matrix}

(176)

\begin{matrix} GA (1, 2) ≅ M_{2} (C) \end{matrix}

(177)

\begin{matrix} GA (0, 3) ≅ H^{2} \end{matrix}

(178)

\begin{matrix} GA (4, 0) ≅ M_{2} (H) \end{matrix}

(179)

\begin{matrix} GA (1, 3) ≅ M_{2} (H) \end{matrix}

(180)

\begin{matrix} GA (0, 4) ≅ M_{2} (H) \end{matrix}

(181)

\begin{matrix} GA (5, 0) ≅ M_{2}^{2} (H) \end{matrix}

(182)

The determinant of these objects is valued in

C

or in

H

, where

C

are the complex numbers, and where

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:

$GA (1, 0)$ :: Let $ψ = a + b e_{1}$ , then:

$\begin{matrix} {(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{matrix}$

(183)

which is valued in $R$ .
$GA (1, 1)$ :: Let $ψ = a + b e_{0} e_{1}$ , then:

$\begin{matrix} {(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{matrix}$

(184)

which is valued in $R$ .
$GA (2, 2)$ :: Let $ψ = a + b e_{0} e_{\emptyset} e_{1} e_{2}$ , where $e_{0}^{2} = - 1, e_{\emptyset}^{2} = - 1, e_{1}^{2} = 1, e_{2}^{2} = 1$ , then:

$\begin{matrix} {⌊ {(a + b)}^{‡} (a + b) ⌋}_{3, 4} {(a + b)}^{‡} (a + b) \end{matrix}$

(185)

$\begin{matrix} = {⌊ a^{2} + 2 a b + b^{2} ⌋}_{3, 4} (a^{2} + 2 a b + b^{2}) \end{matrix}$

(186)

We note that $b^{2} = b^{2} e_{0} e_{\emptyset} e_{1} e_{2} e_{0} e_{\emptyset} e_{1} e_{2} = b^{2}$ , therefore:

$\begin{matrix} = (a^{2} + b^{2} - 2 a b) (a^{2} + b^{2} + 2 a b) \end{matrix}$

(187)

$\begin{matrix} = {(a^{2} + b^{2})}^{2} - 4 a^{2} b^{2} \end{matrix}$

(188)

$\begin{matrix} = {(a^{2} + b^{2})}^{2} - 4 a^{2} b^{2} \end{matrix}$

(189)

which is valued in $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{matrix} GA (0, 1) : & φ^{†} φ \end{matrix}

(190)

\begin{matrix} GA (2, 0) : & φ^{‡} φ \end{matrix}

(191)

\begin{matrix} GA (3, 0) : & {⌊ φ^{‡} φ ⌋}_{3} φ^{‡} φ \end{matrix}

(192)

\begin{matrix} GA (3, 1) : & {⌊ φ^{‡} φ ⌋}_{3, 4} φ^{‡} φ \end{matrix}

(193)

\begin{matrix} GA (4, 1) : & {({⌊ φ^{‡} φ ⌋}_{3, 4} φ^{‡} φ)}^{†} ({⌊ φ^{‡} φ ⌋}_{3, 4} φ^{‡} φ) \end{matrix}

(194)

Theorem 18

(No Metric Measurements). This obstruction applies to

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{matrix} g_{μ ν} = \frac{1}{2} (e_{μ} e_{ν} + e_{ν} e_{μ}) \end{matrix}

(195)

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

GA (2, 0)

cannot produce a metric measurement as a quantum observable, thus its geometry is not observationally complete with respect to its geometry. □

Conjecture 1.

(No probability measures as a self-product (in 6D)). The multivector representation of the norm in 6D cannot satisfy any observables.

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{matrix} s (B) = b_{1} B f_{5} (f_{4} (B) f_{3} (f_{2} (B) f_{1} (B))) + b_{2} B g_{5} (g_{4} (B) g_{3} (g_{2} (B) g_{1} (B))) \end{matrix}

(196)

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{matrix} a_{0}^{4} - 2 a_{0}^{2} a_{47}^{2} + b_{2} a_{0}^{2} a_{47}^{2} p_{412} p_{422} + 〈 72 monomials 〉 = 0 \end{matrix}

(197)

\begin{matrix} b_{1} a_{0}^{3} a_{52} + 2 b_{2} a_{0} a_{47}^{2} a_{52} p_{412} p_{422} p_{432} p_{442} p_{452} + 〈 72 monomials 〉 = 0 \end{matrix}

(198)

\begin{matrix} 〈 74 monomials 〉 = 0 \end{matrix}

(199)

\begin{matrix} 〈 74 monomials 〉 = 0 \end{matrix}

(200)

From Equation 196, it is possible to see that no observable

O

can satisfy this equation because the linear combination does not allow one to factor it out of the equation.

\begin{matrix} b_{1} O B f_{5} (f_{4} (B) f_{3} (f_{2} (B) f_{1} (B))) + b_{2} B g_{5} (g_{4} (B) g_{3} (g_{2} (B) g_{1} (B))) = b_{1} B f_{5} (f_{4} (B) f_{3} (f_{2} (B) f_{1} (B))) + b_{2} O B g_{5} (g_{4} (B) g_{3} (g_{2} (B) g_{1} (B))) \end{matrix}

(201)

Any equality of the above type between

b_{1} O

and

b_{2} O

is frustrated by the factors

b_{1}

and

b_{2}

, forcing

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

GA (1, 3)

even if in relativity, the signature of the metric

(+, -, -, -)

versus

(-, -, -, +)

does not influence the physics. However, in geometric algebra,

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.

GA (1, 3) ≅ M_{2} (H)

and

det M_{2} (H) \in H

Consequently, 3+1D is the only dimensional configuration (other than the "non-geometric" configurations of

GA (0) ≅ R

and

GA (0, 1) ≅ 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

ρ

. 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{matrix} Constraint & Vanishing Phase & Inferred Theory & Wavefunction \\ \bar{E} = \sum_{i} ρ_{i} E_{i} & none & SM & N . A . \\ 0 = tr \sum_{i} ρ_{i} [\begin{matrix} 0 & - E_{i} \\ E_{i} & 0 \end{matrix}] & U (1) & QM & C^{n} \\ 0 = \frac{1}{2} tr \sum_{i} ρ_{i} M_{i} & {Spin}^{c} (3, 1) & RQM / QG & {(R \times {Spin}^{c} (3, 1))}^{n} \end{matrix}

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.

Interpretation

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 dimensional configurations, as demonstrated in Section 2.4.

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 dimensional 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.

Data Availability Statement

No datasets were generated or analyzed during the current study. During the preparation of this manuscript, we utilized a Large Language Model (LLM), for assistance with spelling and grammar corrections, as well as for minor improvements to the text to enhance clarity and readability. This AI tool did not contribute to the conceptual development of the work, data analysis, interpretation of results, or the decision-making process in the research. Its use was limited to language editing and minor textual enhancements to ensure the manuscript met the required linguistic standards.

Conflicts of Interest

The author declares that he has no competing financial or non-financial interests that are directly or indirectly related to the work submitted for publication.

Appendix A. SM

Here, we solve the Lagrange multiplier equation of SM.

\begin{matrix} L = \underset{Boltzmann Entropy}{\underset{︸}{- k_{B} \sum_{i} ρ_{i} ln ρ_{i}}} + \underset{Normalization Constraint}{\underset{︸}{λ (1 - \sum_{i} ρ_{i})}} + \underset{Average Energy Constraint}{\underset{︸}{β (\bar{E} - \sum_{i} ρ_{i} E_{i})}} \end{matrix}

(A1)

We solve the maximization problem as follows:

\begin{matrix} 0 & = \frac{\partial L (ρ_{i}, \dots, ρ_{n})}{\partial ρ_{i}} \end{matrix}

(A2)

\begin{matrix} = - ln ρ_{i} - 1 - λ - β E_{i} \end{matrix}

(A3)

\begin{matrix} = ln ρ_{i} + 1 + λ + β E_{i} \end{matrix}

(A4)

\begin{matrix} \Rightarrow ln ρ_{i} & = - 1 - λ - β E_{i} \end{matrix}

(A5)

\begin{matrix} \Rightarrow ρ_{i} & = exp (- 1 - λ) exp (- β E_{i}) \end{matrix}

(A6)

\begin{matrix} = \frac{1}{Z (τ)} exp (- β E_{i}) \end{matrix}

(A7)

The partition function, is obtained as follows:

\begin{matrix} 1 & = \sum_{i} exp (- 1 - λ) exp (- β E_{i}) \end{matrix}

(A8)

\begin{matrix} \Rightarrow {(exp (- 1 - λ))}^{- 1} & = \sum_{i} exp (- β E_{i}) \end{matrix}

(A9)

\begin{matrix} Z (τ) & : = \sum_{i} exp (- β E_{i}) \end{matrix}

(A10)

Finally, the probability measure is:

\begin{matrix} ρ_{i} = \frac{1}{\sum_{i} exp (- β E_{i})} exp (- β E_{i}) \end{matrix}

(A11)

Appendix B. RQM in 3+1D

\begin{matrix} L = \underset{Relative Shannon Entropy}{\underset{︸}{- \sum_{i} ρ_{i} ln \frac{ρ_{i}}{p_{i}}}} + \underset{Normalization Constraint}{\underset{︸}{λ (1 - \sum_{i} ρ_{i})}} + \underset{Vanishing Relativistic - Phase Anti - Constraint}{\underset{︸}{ζ (- tr \frac{1}{2} \sum_{i} ρ_{i} M_{i})}} \end{matrix}

(A12)

The solution is obtained using the same step-by-step process as the 2D case, and yields:

\begin{matrix} ρ_{i} = \underset{Spin ​^{c} (3, 1) Invariant Ensemble}{\underset{︸}{\frac{1}{\sum_{i} p_{i} det exp (- ζ \frac{1}{2} M_{i})}}} \underset{{Spin}^{c} (3, 1) Born Rule}{\underset{︸}{det exp (- ζ \frac{1}{2} M_{i})}} \underset{Initial Preparation}{\underset{︸}{p_{i}}} \end{matrix}

(A13)

Proof.

The Lagrange multiplier equation can be solved as follows:

\begin{matrix} 0 & = \frac{\partial L (ρ_{1}, \dots, ρ_{n})}{\partial ρ_{i}} \end{matrix}

(A14)

\begin{matrix} = - ln \frac{ρ_{i}}{p_{i}} - p_{i} - λ - ζ tr \frac{1}{2} M_{i} \end{matrix}

(A15)

\begin{matrix} = ln \frac{ρ_{i}}{p_{i}} + p_{i} + λ + ζ tr \frac{1}{2} M_{i} \end{matrix}

(A16)

\begin{matrix} \Rightarrow ln \frac{ρ_{i}}{p_{i}} & = - p_{i} - λ - ζ tr \frac{1}{2} M_{i} \end{matrix}

(A17)

\begin{matrix} \Rightarrow ρ_{i} & = p_{i} exp (- p_{i} - λ) exp (- ζ tr \frac{1}{2} M_{i}) \end{matrix}

(A18)

\begin{matrix} = \frac{1}{Z (ζ)} p_{i} exp (- ζ tr \frac{1}{2} M_{i}) \end{matrix}

(A19)

The partition function

Z (ζ)

, serving as a normalization constant, is determined as follows:

\begin{matrix} 1 & = \sum_{i} p_{i} exp (- p_{i} - λ) exp (- ζ tr \frac{1}{2} M_{i}) \end{matrix}

(A20)

\begin{matrix} \Rightarrow {(exp (- p_{i} - λ))}^{- 1} & = \sum_{i} p_{i} exp (- ζ tr \frac{1}{2} M_{i}) \end{matrix}

(A21)

\begin{matrix} Z (ζ) & : = \sum_{i} p_{i} exp (- ζ tr \frac{1}{2} M_{i}) \end{matrix}

(A22)

□

Appendix C. SageMath program showing ⌊u ‡ u⌋ 3,4 u ‡ u=detM u

[H] from sage.algebras.clifford_algebra import CliffordAlgebra

from sage.quadratic_forms.quadratic_form import QuadraticForm

from sage.symbolic.ring import SR

from sage.matrix.constructor import Matrix

# Define the quadratic form for GA(3,1) over the Symbolic Ring

Q = QuadraticForm(SR, 4, [-1, 0, 0, 0, 1, 0, 0, 1, 0, 1])

# Initialize the GA(3,1) algebra over the Symbolic Ring

algebra = CliffordAlgebra(Q)

# Define the basis vectors

e0, e1, e2, e3 = algebra.gens()

# Define the scalar variables for each basis element

a = var(’a’)

t, x, y, z = var(’t x y z’)

f01, f02, f03, f12, f23, f13 = var(’f01 f02 f03 f12 f23 f13’)

v, w, q, p = var(’v w q p’)

b = var(’b’)

# Create a general multivector

udegree0=a

udegree1=t*e0+x*e1+y*e2+z*e3

udegree2=f01*e0*e1+f02*e0*e2+f03*e0*e3+f12*e1*e2+f13*e1*e3+f23*e2*e3

udegree3=v*e0*e1*e2+w*e0*e1*e3+q*e0*e2*e3+p*e1*e2*e3

udegree4=b*e0*e1*e2*e3

u=udegree0+udegree1+udegree2+udegree3+udegree4

u2 = u.clifford_conjugate()*u

u2degree0 = sum(x for x in u2.terms() if x.degree() == 0)

u2degree1 = sum(x for x in u2.terms() if x.degree() == 1)

u2degree2 = sum(x for x in u2.terms() if x.degree() == 2)

u2degree3 = sum(x for x in u2.terms() if x.degree() == 3)

u2degree4 = sum(x for x in u2.terms() if x.degree() == 4)

u2conj34 = u2degree0+u2degree1+u2degree2-u2degree3-u2degree4

mdegree0 = a

mdegree1 = t*y0+x*y1+y*y2+z*y3

mdegree2 = f01*y0*y1+f02*y0*y2+f03*y0*y3+f12*y1*y2+f13*y1*y3+f23*y2*y3

mdegree3 = v*y0*y1*y2+w*y0*y1*y3+q*y0*y2*y3+p*y1*y2*y3

mdegree4 = b*y0*y1*y2*y3

m=mdegree0+mdegree1+mdegree2+mdegree3+mdegree4

print(u2conj34*u2 == m.det())

The program outputs

True

showing, by computer assisted symbolic manipulations, that the determinant of the real Majorana representation of a multivector u is equal to the double-copy form:

det M_{u} = {⌊ u^{‡} u ⌋}_{3, 4} u^{‡} u

References

Dirac, P.A.M. The principles of quantum mechanics; Number 27, Oxford university press, 1981.
Von Neumann, J. Mathematical foundations of quantum mechanics: New edition; Vol. 53, Princeton university press, 2018.
Jaynes, E.T. Information theory and statistical mechanics. Physical review 1957, 106, 620. [Google Scholar] [CrossRef]
Jaynes, E.T. Information theory and statistical mechanics. II. Physical review 1957, 108, 171. [Google Scholar] [CrossRef]
Kullback, S.; Leibler, R.A. On information and sufficiency. The annals of mathematical statistics 1951, 22, 79–86. [Google Scholar] [CrossRef]
Shannon, C.E. A mathematical theory of communication. Bell system technical journal 1948, 27, 379–423. [Google Scholar] [CrossRef]
Hestenes, D. Spacetime physics with geometric algebra. American Journal of Physics 2003, 71, 691–714. [Google Scholar] [CrossRef]
Lundholm, D. Geometric (Clifford) algebra and its applications. arXiv preprint math/0605280 2006.
Hestenes, D. Space-time structure of weak and electromagnetic interactions. Foundations of Physics 1982, 12, 153–168. [Google Scholar] [CrossRef]
Lasenby, A. Some recent results for SU(3) and Octonions within the Geometric Algebra approach to the fundamental forces of nature. 2022. arXiv:2202.06733 2022.
Bern, Z.; Carrasco, J.J.M.; Johansson, H. Perturbative quantum gravity as a double copy of gauge theory. Physical Review Letters 2010, 105, 061602. [Google Scholar] [CrossRef] [PubMed]
Acus, A.; Dargys, A. Inverse of multivector: Beyond p+ q= 5 threshold. 2017. arXiv:1712.05204 2017.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Physics as an Optimization Problem over All Measurements

Abstract

1. Introduction

1.1. Quantum Mechanics

1.2. Relativistic Quantum Mechanics

1.3. Quantum Gravity

1.4. Dimensional Obstructions

2. Results

2.1. Quantum Mechanics

2.2. RQM in 2D

Obstructions

2.3. RQM in 3+1D

2.3.1. Preliminaries

2.3.2. RQM

2.3.3. Standard Model Gauge Symmetries

2.3.4. Quantum Gravity

2.3.5. Linearized Gravity (Sketch)

2.4. Dimensional Obstructions

3. Discussion

Interpretation

4. Conclusion

Data Availability Statement

Conflicts of Interest

Appendix A. SM

Appendix B. RQM in 3+1D

Appendix C. SageMath program showing ⌊u ‡ u⌋ 3,4 u ‡ u=detM u

References

MDPI Initiatives

Important Links

Subscribe