Spacetime geometry of acoustics and electromagnetism

Both acoustics and electromagnetism represent measurable fields in terms of dynamical potential fields. Electromagnetic force-fields form a spacetime bivector that is represented by a dynamical energy-momentum 4-vector potential field. Acoustic pressure and velocity fields form an energy-momentum density 4-vector field that is represented by a dynamical action scalar potential field. Surprisingly, standard field theory analyses of spin angular momentum based on these traditional potential representations contradict recent experiments, which motivates a careful reassessment of both theories. We analyze extensions of both theories that use the full geometric structure of spacetime to respect essential symmetries enforced by vacuum wave propagation. The resulting extensions are geometrically complete and phase-invariant (i.e., dual-symmetric) formulations that span all five grades of spacetime, with dynamical potentials and measurable fields spanning complementary grades that are related by a spacetime vector derivative (i.e., the quantum Dirac operator). These complete representations correct the equations of motion, energy-momentum tensors, forces experienced by probes, Lagrangian densities, and allowed gauge freedoms, while making manifest the deep structural connections to relativistic quantum field theories. Finally, we discuss the implications of these corrections to experimental tests.


Introduction
Recent work in locally measuring the angular momentum of acoustic and optical fields using small probes [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15] has prompted a reexamination of the theoretical treatment of intrinsic spin in relativistic fields.The corrections needed for the field theory calculations to agree with experiment have important implications for the foundations of both classical and quantum field theory descriptions of natural phenomena.In particular, dynamical potentials that represent the measurable fields in each theory must be generalized in a geometrically motivated way that respects structural symmetries.
The energy-momentum and angular momentum tensors that are predicted in classical field theory follow from Noether's celebrated theorem that identifies conserved quantities from the continuous symmetries of the Lagrangian [16], so their prediction critically depends upon the functional structure of the Lagrangian density for the theory.The canonical tensors derived from the Noether theorem, however, contain contributions from intrinsic spin angular momentum that are neglected in standard textbook treatments [16,17].In each reference frame, the total angular momentum J = L + S splits into two distinct parts: an extrinsic orbital angular momentum part L = r × P involving the conserved canonical momentum P , and an intrinsic spin angular momentum part S that is formally local to the infinitesimal volume around each point of the field.Similarly, the derived canonical momentum P = P K + P S splits into a kinetic contribution P K and a spin-dependent contribution P S = −( ∇ × S)/2 involving the same spin vector S [18,19,20,21,22].Traditionally, the local spindependent parts of the canonical momentum and angular momentum are removed using a symmetrization procedure developed by Belinfante [23], leaving behind just the kinetic momentum P K and its contribution to the orbital angular momentum r × P K .The historical justification for neglecting the spin contributions has been that localized angular momentum integrates to zero over a volume due to Stokes' theorem, so does not contribute to the integrated total angular momentum of the field [16,23].
However, in the past few decades it has become clear that this omission of the predicted spin contributions was premature.Localized probes occupying sufficiently small volumes can have boundaries with suitably broken symmetries that they retain nonzero spin contributions to the integrated momentum and angular momentum transferred from the field.The spin contributions to the net force and torque on the local probe can therefore be physically measured.Indeed, local spin angular momentum densities have recently been measured for both acoustic [1,2,3,5,6] and electromagnetic [7,24,8,9,10,11,12,13,14,15] fields, allowing a direct check for the field theory predictions.Surprisingly, the historically accepted Lagrangian formulations have failed to correctly predict these measured experimental results for the spin angular momentum densities in both acoustics and electromagnetism, despite their many other successes over a century of use.
In the case of acoustics, a vanishing spin angular momentum is predicted by the traditional approach using a dynamical action scalar potential [16,25,26,27].Even without the experimental refutation of this prediction, microscopic considerations readily show why this prediction is inadequate.The isotropic medium in acoustics is composed of individual molecules bouncing around in a chaotic way, with the acoustic field describing the mean time-averaged molecular motion.If a localized region of the medium has molecules undergoing elliptical orbits on average, one would expect the corresponding mean field to have intrinsic spin in those regions that locally describes the mean angular momentum of those underlying elliptical orbits [2,3,28,29].Creating such a situation physically is straightforward, since propagating sound waves cause longitudinal oscillations of the molecules on average, so superposing two sound waves propagating in perpendicular directions will result in elliptical mean orbits locally with the ellipticity dependent upon the relative phase of the two sound waves.Indeed, experimental measurements used precisely such an acoustic setup to measure nonzero acoustic spin with local probes [1,2,5,6], thus contradicting traditional field theory predictions.
In the case of electromagnetism, an electrically biased spin angular momentum is predicted by the traditional Lagrangian approach that uses electric scalar and vector potentials describing energy and momentum per unit electric charge [18,19,20,21,22], meaning that the predicted spin vector depends only upon the local electric field and not the magnetic field.However, optical fields far from sources show an egalitarianism between the electric and magnetic fields, known as dual (electricmagnetic) symmetry [30,7,31,32,33,18,19,21], that is related to the conserved helicity of the propagating optical field [34,21].An electrically biased prediction thus violates the expected symmetries of the propagating field.Indeed, contributions from both the electric and magnetic fields are required to correctly predict the measured spin angular momentum [7,18,19,10,11,21].
Fixing the Lagrangian approach to these theories so that they correctly predict the measured spin angular momentum densities is a subtle affair, since the many successes of the existing theories tightly constrain any proposed modification.The key to resolving the dilemma is to reassess how the measurable fields are represented in terms of dynamical potentials.Changing the potential representation alters the functional form of the Lagrangian and thus the predicted conserved quantities from Noether's theorem.However, there is substantial gauge freedom in such potential representations that leave the equations of motion or other predictions of the theory invariant.
In electromagnetism, a dual-symmetric representation involving both electric and magnetic versions of the scalar and vector potentials successfully describes the measured spin angular momentum [7,31,32,18,34,19,35,21], which causes an interesting tension with the prior neglect of magnetic potentials that has been motivated by the apparent lack of magnetic charge.In acoustics, the correct spin angular momentum was first calculated from microscopic arguments [2,27,3], and agreed with the experimentally measured results.The field theory justification came afterwards [29] and modified the representation by augmenting the scalar potential with a second bivector potential field (i.e., with a structure similar to the electromagnetic field) containing information about both the mean displacement and rotational vorticity of the underlying acoustic medium.To agree with experimental measurements of spin, the representations of the measurable fields in each theory must involve at least two dynamical potentials.
In this paper, we explore the important role played by the dynamical potentials in both acoustics and elec-tromagnetism, while highlighting their many structural analogies for clarity.We build on our previous work [21,29] to explain the emergence of the multiple dynamical potentials from fundamentally geometric considerations that both constrain and inform any needed modifications of the traditional theory.As part of this goal, we deliberately minimize reliance upon particular Lagrangian treatments and instead focus on understanding which quantities necessarily appear from considerations of the geometry alone.Both acoustics and electromagnetism describe waves with constant speed, so naturally adhere to the symmetries enforced by the geometry of relativistic spacetime, which tightly constrains the allowed structure of the possible potential representations.We find that expressing both theories in the geometric language of a Clifford bundle over spacetime significantly clarifies the structure of each theory while illuminating both their many similarities and their key geometric differences.
We show that the physical content of both acoustics and electromagnetism is tightly constrained when the entire geometry of spacetime is taken into account.The expanded potential representations and other critical modifications to each theory can be motivated without appeal to a particular Lagrangian, which in turn constrains which choices of Lagrangian are physically consistent.Gauge freedoms of each theory also become tightly constrained, with some traditional gauge choices actually modifying the equations of motion and the forces experienced by probes, making them experimentally testable and thus not true gauge freedoms.In particular, causal gauge choices (like the Lorenz-FitzGerald contraint in electromagnetism) become experimentally motivated, rather than an optional choice.
We also highlight that the geometric language used here for classical acoustic and electromagnetic field theories has thought-provoking connections to relativistic fields more broadly.Issues with potential representations, gauge fixing, and degenerate Lagrangian constructions have plagued quantum theories as much as classical theories for many decades.Our work here may shed some light on possible routes forward towards more geometrically complete treatments that resolve several long-standing issues in both both classical and quantum relativistic field theories.

Acoustic and electromagnetic theories
We first review the traditional formulations of linear acoustics and electromagnetism to highlight their similarities and geometric structure.

Acoustic equations in 3D
Recall that the measurable fields in linear acoustics are the pressure (P ) and velocity ( v) densities of an isotropic medium with equilibrium mass-density ρ and compressibility β.The acoustic equations of motion for small perturbations are then, with an extra irrotational constraint, Note that ρ v is a momentum per volume and P is an energy per volume, so the first two equations can be interpreted as the force and power density constraint equations for the medium, while the third constraint guarantees longitudinal flow.The Lagrangian density, has familiar forms of kinetic and potential energy densities that are analogous to a spring.Taking time derivatives and rearranging the equations of motion yields two wave equations, for the momentum and energy densities of the medium, which describe acoustic energy and momentum density waves with a constant wave speed c ≡ 1/ √ ρβ.

Electromagnetism in 3D
Similarly, recall that the measurable fields in electromagnetism are the electric E = D/ǫ and magnetic H = B/µ fields in a medium with permittivity ǫ and permeability µ.The equations of motion are, with divergence-free constraints imposing transverse flow, The Lagrangian density has a form similar to acoustics, Taking time derivatives and rearranging the equations of motion and constraints yields two wave equations, that describe electromagnetic waves with a constant wave speed c ≡ 1/ √ ǫµ.

Traditional potential representations
Curiously, both theories cannot derive their respective equations of motion using the Lagrangians as defined in Eqs. ( 3) and (7) in terms of the measurable fields directly.Instead, they must introduce related dynamical potential fields from which the measurable fields are derived.These potential fields must be varied in the Lagrangian in order to produce the correct equations of motion.Thus, at least from the perspective of the field Lagrangians, the measurable fields are less fundamental than the dynamical potential fields.
In acoustics, it is customary to introduce a scalar potential φ, satisfying The acoustic scalar field φ has units of action density, so these definitions intuitively match the Hamilton-Jacobi definitions of energy and momentum from an action.In electromagnetism it is customary to introduce a scalar potential φ e and vector potential A e , satisfying The electromagnetic scalar and vector potentials (φ e , A e ) have units of energy and momentum per unit charge, respectively.
The electromagnetic potentials have an additional gauge freedom such that the transformations, φ e → φ e − ∂ t χ and A e → A e + ∇χ, for some scalar field χ, will leave the measurable electromagnetic fields ( E, B) invariant.This freedom generalizes the introductory physics observation that one can shift the zero of potential energy without affecting the measurable energy differences.Despite this gauge freedom, however, one still must vary the potential fields in the Lagrangians after making these substitutions in order to derive the equations of motion.

The problem with spin angular momentum
Getting the correct equations of motion is not the only purpose of a Lagrangian formulation of the theories, however.A Lagrangian also enables the identification of the conserved quantities of the field from its continuous symmetries according to Noether's theorem.For example, the energy-momentum stress tensor arises from translational symmetry, while the angular momentum tensor arises from from rotational symmetries.These conserved quantities can be experimentally verified, providing a check on the validity of the framework.
Without dwelling on the derivation details, the traditional formulations of both theories predict the following.For electromagnetism, the derived spin angular momentum tensor is fully characterized by a spin vector of local angular momentum.Choosing the radiation gauge that sets the scalar potential (energy per charge) to zero everywhere φ e = 0 to focus solely on the vector potential (momentum per charge) yields a spin vector with a simple form [16,18,9,20,21,22,13], that asymmetrically depends on only the electric field E, not the magnetic field, and also depends on the vector potential A e itself.Despite appearances, this spin vector is in fact gauge-invariant and measurable, since the transverse part of the vector potential that contributes to the local spin is not affected by the choice of gauge in the reference frame of the probe [21].For monochromatic fields of frequency ω, in particular, the cycle-averaged spin vector can be expressed with the cycle-averaged complex representation E of the electric field alone as S e = ǫ Im(E * × E)/2ω, which avoids the awkward question of apparent gauge-dependence [18,21,7,9,10,11].In contrast, acoustics is described with a scalar potential that predicts a vanishing spin vector [16,25,26,27], Neither prediction is experimentally correct!For electromagnetism, both the electric and magnetic fields contribute symmetrically to the total spin angular momentum when carefully probed in the laboratory far from sources [7,24,8,9,10,11,12,13,14], It also involves both electric and magnetic vector potentials, or their cycle-averaged monochromatic forms, ferent choices of probe have been used to isolate and verify the distinct electric and magnetic contributions to this local spin vector [15].The traditionally predicted electric-biased spin vector is experimentally refuted.For acoustics, recent measurements of intrinsic spin show it to be nonzero [1,2,27,3,5,6], and indeed match our recently derived spin-vector expression [29], that involves a mean displacement potential x satisfying v = ∂ t x in a vorticity-free choice of gauge.The cycleaveraged and monochromatic version of this modified spin density, S a = ρ Im(v * × v)/4ω, was also microscopically derived [2,27,3,28] by observing that small circular orbits of the molecules in the medium should appear like localized spin angular momentum after taking the mean field continuum limit.The vanishing spin vector predicted by the traditional scalar potential is thus experimentally refuted.
These experimental contradictions force the reevaluation of the formalism for both theories so that they naturally support these conclusions about the intrinsic spin.From the perspective of 3D space, however, it is not obvious how to resolve this theoretical dilemma.The following sections review the 4D spacetime approach to both theories, which will make the structural gaps in the traditional theories more readily apparent.By analyzing both theories without appeal to a Lagrangian formalism a priori, it will become clear how the geometry of spacetime places strong constraints on the admissible structure, even without focusing specifically on these measured contradictions regarding the spin angular momentum.The geometric constraints of spacetime give an important clue for how to proceed rigorously to resolve these longstanding theoretical dilemmas.

Spacetime formulations
Accommodating the natural symmetries of electromagnetic waves with constant speed c was the original motivation for introducing the framework of spacetime with a Minkowski metric in special relativity.That framework naturally unifies both electromagnetic fields into a single bivector F and unifies both energy and momentum into a single 4-vector.Notably, acoustics has a very similar structure, but with a different wave speed c, and even contains a natural pair of energy and momentum densities, which similarly motivates the construction of an analogous acoustic spacetime that respects the symmetries of the wave equation in the acoustic medium [36].Introducing similar spacetime geometries in both theories considerably clarifies their common structure.
The following treatment expresses fields as part of the tangent bundle of Clifford algebras over spacetime [37,38,39,40,41,42,43,44,45,46,47,48,49,50], since this mathematical formalism helps highlight the geometric structure of the theories in a particularly efficient and illuminating way while maintaining a clear correspondence to standard 3D vector calculus and differential forms.A brief overview of this mathematical framework is provided in Appendix A for completeness, with a more expanded pedagogical treatment in Ref. [21].Tables 1 and 2 also summarize the essential features of the algebras as a quick reference.µ=0 with Minkowski metric γµ •γν = ηµν of signature (+, −, −, −).The Clifford product γµγν = γµ • γν + γµ ∧ γν unifies the metric and Grassmann wedge product as the symmetric and antisymmetric parts of the same associative and invertible product, so orthogonal vectors naturally anticommute γµγν = −γν γµ.The pseudoscalar I ≡ γ 0 γ 1 γ 2 γ 3 satisfies I 2 = −1, commutes with even grades, and anticommutes with odd grades.Multiplication by I enacts the Hodge-star operation that maps to the geometric complement.The grade-2 units γ k γ 0 = γ k ∧ γ 0 ≡ σ k are space-time planes that are perceived as 3-dimensional spatial unit vectors in a particular reference frame being dragged along the temporal direction γ 0 , so are notated with 3-vector arrows accordingly.The grade-2 units like γ 2 γ 3 = γ 2 ∧ γ 3 ≡ −I σ 1 = − σ 2 σ 3 are spatial planes orthogonal to an effective 3-vector σ 1 in the frame γ 0 .The Lorentz group is the group of rotations in spacetime planes, with each unit plane generating a boost or spatial rotation upon exponentiation.Notably, the Dirac matrices used to represent Dirac spinors in relativistic quantum theory are matrix representations of the unit vectors γµ in spacetime Clifford algebra.

Acoustic spacetime
In the case of acoustics, the pressure P is an energy density, while ρ v is a momentum density, so should naturally combine into an energy-momentum density 4-vector that keeps the equilibrium speed of sound c = 1/ √ ρβ invariant.However, this mapping should be done carefully so that role of the background medium is correctly preserved.When the medium is completely at rest and at equilibrium there is still an equilibrium pressure P 0 = ρc 2 in the medium that acts as a background energy density.Denoting γ 0 as the distinguished timelike (γ 2 0 = 1) unit vector of the rest frame, the equilibrium 4-momentum should be a constant pressure at every point in the medium, with zero velocity, Now consider a reference frame of a pointlike observer moving through the medium at a velocity v relative to that equilibrium frame.Defining the boost angle tanh α = | v|/c and the unit space-time plane v = v/| v| for the boost rotation, a Lorentz transformation to the frame of this observer is represented by a half-angle spinor exp(−αv/2) that acts as a double-sided group The even-graded subalgebra of the spacetime algebra is closed and has 2 3 graded basis elements.Choosing a particular reference frame with timelike direction γ 0 splits the 6 planar directions of spacetime into 3 space-time planes { σ k ≡ γ k γ 0 } 3 k=1 and 3 purely spatial planes {I σ k = −ǫ kij γ i γ j } 3 k=1 .The directions σ k are perceived by the evolving frame as spatial 3D unit vectors since each γ k is being dragged along the temporal axis γ 0 .The apparent unit 3D volume element is thus counter-intuitively identical to the same unit 4D spacetime volume element in every reference frame, I = σ 1 σ 2 σ 3 = γ 0 γ 1 γ 2 γ 3 , which satisfies I 2 = −1 and commutes with the entire 3D subalgebra.Each spatial plane I σ k is geometrically complementary to an orthogonal 3D spatial axis σ k and is oriented to rotate around that axis in accordance with the right-hand rule by convention, indicated algebraically by (I σ 1 )(I σ 2 )(I σ 3 ) = I 4 = 1.They are notably related to the quaternionic imaginary units, which are the spatial unit planes with left-handed orientation, I σ 1 ≡ i, I σ 2 ≡ −j, I σ 3 ≡ k, indicated algebraically by i 2 = j 2 = k 2 = ijk = −1.Each 3D Clifford subalgebra is thus equivalent to a biquaternion algebra that augments the quaternions with an extra imaginary unit I. Notably, the Pauli matrices used to represent both nonrelativistic spin and lightcone Weyl spinors are matrix representations of the unit vectors σ k in 3D Clifford algebra that conflate the pseudoscalar I with a generic scalar imaginary unit i.

Signature
+ grade-3 automorphism on any element of the spacetime algebra.Performing this hyperbolic boost rotation by angle α in the unit plane v on the 4-vector p 0 yields, e −αv/2 p 0 e αv/2 = P 0 c e −αv γ 0 , The moving observer sees an increased pressure, P ′ = γ v P 0 , with dilation factor as well as a momentum density p ′ = −γ v ρ v of the medium flowing past the observer at a velocity − v with an effectively increased mass density ρ = P ′ /c 2 .This increase in apparent pressure and mass density makes sense physically, since the motion of the observer shortens the effective distances between particles in the medium on average, as long as the inertial observer is treated as pointlike to avoid symmetry-breaking drag from particle collisions while in motion.In the limit that the observer moves through the medium near the speed of sound, the background medium will appear to have no space between its particles, so the effective density and pressure of the medium will appear to diverge for the moving observer.At the speed of sound c, the observer reaches a shock wave discontinuity where on average all particles seem to have zero spacing between them and the spacetime description fails.Thus, from an equilibrium perspective this spacetime formulation of acoustics describes the right behavior for speeds below the speed of sound c.
The pressure P (x) and velocity v(x) density fields considered in the 3D case are then deviations away from this equilibrium background of the medium while still in the stationary reference frame γ 0 of the medium, making the total pressure field P tot (x) = P 0 + P (x).Crucially, this construction allows for the deviation pressure P to become locally negative during pressure wave propagation, physically indicating a local decrease of the background particle density relative to the equilibrium state of the medium as a locally negative energy density.The total energy-momentum 4-vector should then have the form p tot ≡ p 0 + p, with deviation, expressed here as a paravector in the reference frame of the medium γ 0 .As a deviation from equilibrium, the pressure P and velocity v are independent components of the energy-momentum deviation field p with a relationship that will be determined by their sources.Thus, the deviation field need not have any particular signature a priori, with p 2 = (P/c) 2 − |ρ v| 2 allowed to be positive, negative, or zero.Given this spacetime structure, the source-free acoustic equation of motion and Lagrangian density for this deviation field then take the laconic forms, where γ 0 is the spacetime vector derivative (known as the Dirac operator in relativistic quantum mechanics).The acoustic equation is thus notably identical in form to the massless Dirac equation, but involving the 4-vector p.
The equation of motion in Eq. ( 18) expands into three independent grades in the rest frame, thus reproducing the expected equations.Here I is the pseudoscalar (unit 4-volume) of spacetime satisfying Moreover, taking another derivative immediately produces a scalar wave equation, Similarly, the acoustic scalar potential satisfies, which implies that the original equation of motion is itself a wave equation for the scalar potential, The corresponding Lagrangian density similarly expands into the expected form in the rest frame, Achieving mathematical expressions that are so compact rarely happens accidentally in physics, which gives strong support to the formal introduction of acoustic spacetime as being conceptually beneficial.

Electromagnetic spacetime formulation
In the case of electromagnetism in standard spacetime, the electric and magnetic fields naturally combine into an electromagnetic field (Faraday) bivector [21,39,48], expressed here as a complex split into its polar and axial 3-vector parts, with speed of light c = 1/ √ ǫµ.This naturally complex decomposition of the electromagnetic field is equivalent to the Riemann-Silberstein approach [51], and is precisely equal to the single photon wave function in quantum electrodynamics [52,53], with the important replacement of the generic scalar imaginary i with the geometrically meaningful pseudoscalar I to make F properly frame-invariant [39,21].
The spacetime bivector F is geometrically planar, so is one dimension greater than the linear energymomentum density p in acoustics.Nevertheless, the electromagnetic equation of motion and standard Lagrangian density have nearly identical laconic forms to the acoustic Eqs.(18), with the same form as the massless Dirac equation for the quantum electron.
The equation of motion expands into four distinct and independent contributions [21], thus reproducing the expected equations.A second derivative immediately yields the expected wave equation, The Hodge-star complement of the (Faraday) bivector F is its geometrically dual (Maxwell) bivector, such that F = ζ GI with wave impedance ζ.Notably, the roles of the magnetic and electric fields are exchanged for the geometric complement field, yet it obeys identical vacuum equations of motion ∇G = 0, which is a symmetry called (electric-magnetic) dual symmetry [30,7,31,32,33,18,34,9].This dual symmetry extends to sources, provided that the electric and magnetic charges are also swapped in tandem, and plays an important role in constructing a Lagrangian formulation of the theory that correctly predicts the experimentally confirmed conserved Noether currents for both spin density and helicity [21].The electromagnetic scalar and vector potentials combine into an energy-momentum-per-unit-charge 4-vector potential field a e = γ 0 (φ e /c − A e ), satisfying, which reproduces the expected relations, provided that the Lorenz-FitzGerald gauge condition is satisfied, which enforces causal evolution for the potential field.With this causal gauge constraint, the original equation of motion also is itself a wave equation for the potential, The corresponding Lagrangian density also expands into the expected form, but notably requires a different overall sign and an extra scalar projection to neglect the pseudoscalar part of the invariant square [21].Notably, this omitted pseudoscalar part involving E • H has been revived in recent treatments of axion contributions to the field [54,55,56], which is an interesting topic for future study.

Dynamical potential representations
Let us now revisit how to introduce and justify dynamical potential representations of the measurable fields.
Recall that the standard motivation for introducing potential fields in 3D starts from the vector identities curl(grad) = 0 and div(curl) = 0.That is, in acoustics the irrotational constraint in Eq. ( 2) is automatically satisfied if v is the gradient of a scalar field φ.Similarly, the magnetic divergence-free condition in Eq. ( 6) is automatically satisfied if H is the curl of a vector field A e , while Eq. ( 5) implies −∂ t A e must contribute to E, yielding a vanishing curl condition for the remainder that is automatically satisfied by the gradient of a scalar field.In 4D these basic motivations can be considerably generalized in a principled way by examining what is allowed by the geometry of spacetime.

Hodge decomposition into potentials
As detailed in the Appendix A.3, the spacetime vector derivative, ∇ = ∇ • + ∇ ∧ ∼ δ + d, naturally splits into a grade-raising 4-curl, ∇ ∧ , which is completely equivalent to the exterior derivative d on forms, as well as a grade-lowering 4-divergence, ∇ • , which is completely equivalent to the codifferential δ = ⋆ −1 d⋆ on forms, with multiplication by the spacetime unit volume pseudoscalar I being the equivalent of the grade-inverting Hodge star ⋆ operation on forms.Thus, these independent parts of the derivative satisfy the identities, that generalize the 3D vector identities used above and imply that the d'Alembertian (scalar wave operator), has the form of the Hodge Laplacian (δ + d) 2 = δd + dδ.Given these equivalences, the Hodge decomposition theorem from differential forms [57] also holds, which states that any multivector object in the spacetime Clifford algebra has an exact decomposition into a 4-curl of a lower-grade object A − , a 4divergence of a higher-grade object A + , and a harmonic part A 0 annihilated by the Laplacian ∇ 2 A 0 = 0.This theorem generalizes the Helmholtz theorem for 3D vector calculus and justifies the introduction of three potential fields (A − , A 0 , A + ) that are generally needed to represent any measurable field A.

Electromagnetic potentials
Starting with electromagnetism, for which these extensions have been more explored historically, the measurable field F is a grade-2 bivector field, so we expect a complete Hodge decomposition, into a grade-1 4-vector potential a e , a grade-2 harmonic bivector F 0 satisfying ∇ 2 F 0 = 0, and a grade-3 pseudo-4-vector potential a m I, allowing for undetermined constants λ − , λ 0 , λ + , that determine the relative fractions of each potential contribution to the total measurable field F .The 4-vector potential a e is precisely the standard electric vector potential representation in Eq. ( 29), with units of energy-momentum per electric charge.Similarly, the pseudo-4-vector potential, is identifiable as the magnetic potential representation, with impedance ζ as in Eq. ( 28), and units of energymomentum per magnetic charge (using the Amperemeter convention).This potential contains the magnetic scalar φ m and vector A m potentials used in magnetostatics as well as more general treatments of electromagnetism that include magnetic sources.This is precisely the quantity that is needed to reproduce the experimentally confirmed results for the electromagnetic spin density [7,18,8,9,10,11,12,13,14].
To compute what the contribution of this magnetic potential looks like, first consider its full derivative, where product associativity makes the expressions for the grade-2 and grade-4 parts manifestly similar to those from the electric potential, Thus the representation of the electromagnetic field with the pseudo-4-vector potential alone a m I is, which match the expected magnetic expressions.
Similarly to the electric potential, setting the grade-4 part of the derivative to zero, enforces a causal constraint like the Lorenz-FitzGerald constraint in Eq. (30).In this case the field representation takes the analogous laconic form F = ∇a m I and Maxwell's equation becomes a wave equation for the potential, ∇F = (∇ 2 a m )I = 0. Thus, the full Hodge decomposition of the measurable field into the three potential fields has the 3D expansion, with the electric and magnetic potentials contributing naturally dual contributions.The roles of the three potentials is clarified by including sources in the equation of motion, In principle, the sources can have both electric 4-vector charge-currents and magnetic pseudo-4-vector chargecurrents, using the Ampere-meter convention for magnetic charge.
Expanding the equation of motion into potentials yields, which has scalar Laplacian operators and thus shows that the electric potential is only affected by electric charge, while the magnetic potential is only affected by magnetic charge.
In principle the third harmonic field contribution F 0 can couple to both types of charge since ∇F 0 = ∇•F 0 + ∇ ∧ F 0 spans both grades; however, from a Lagrangian perspective, directly varying a bivector field like F 0 as a dynamical field will not reproduce the correct equation of motion and so cannot produce such a contribution coupled to sources.Thus, only two potentials couple to the distinct types of charge independently, while the third potential can only describe a homogeneous background contribution, λ 0 ∇F 0 = 0. Notably, if the homogeneous potential F 0 is dropped as problematic for the Lagrangian formalism, the remaining potentials can be bundled together as a naturally complex 4-vector potential representation of the field F [21], The relative proportionality factors are kept arbitrary here, but will generally become constrained by how the field couples to sources and probes [21,29].We will further explore this subtle issue in future work.

Acoustic potentials
In the case of acoustics, the measurable field p is a grade-1 4-vector field, one grade lower than electromagnetism, so has an asymmetric Hodge decomposition, for a grade-0 scalar potential field φ, a grade-1 harmonic 4-vector contribution p ′ 0 satisfying ∇ 2 p ′ 0 = 0, and a grade-2 bivector potential field M .Analogously to the electromagnetic case, the coefficients λ −,0,+ are arbitrary constants that determine the representation proportions.
The scalar potential contribution is precisely the standard representation in Eq. ( 21).The harmonic contribution p ′ 0 should play a role similar to F 0 in electromagnetism and provide a homogeneous contribution to the field, such as the background pressure p 0 of the medium, but is problematic for a Lagrangian treatment of the dynamics.The representation of an acoustic field with a bivector potential, however, is not commonly seen in the literature.Since p has units of energy-momentum density, φ has units of action density, while M has units of angular momentum density, which are formally the same as action but with an important conceptual difference.The factor of 1/3 arises from the observation that a derivative of the orbital angular momentum, L = p ∧ r, does recover the linear momentum, −∇ • L = 3p, but with an extra factor of 3. The bivector potential relation p = −(∇•M )/3 in Eq. ( 50) anticipates this overcounting so that M can be interpreted as an angular momentum density.
Recall more generally that a relativistic angular momentum bivector has two distinct pieces, M = L + S. The extrinsic orbital angular momentum depends on r and the linear momentum p, and splits in a particular frame into a polar mass-moment vector N L = t p − (E/c 2 ) r and axial rotational angular momentum LI = r × pI.In contrast, the intrinsic spin angular momentum S = c N S + SI is independent of the linear motion.However, it does contribute spin-massmoment N S and spin-angular momentum S contributions to the total angular momentum, The acoustic potential field M in Eq. ( 51) generally describes such a total angular momentum density for the medium.In the rest frame, its polar mass-moment density vector c N = c(ρ x) is a space-time boost rotation plane that is directly proportional to a massdensity displacement field x describing linear shifts of the density away from equilibrium.Its rotational angular momentum J = (ρc) y is a purely spatial rotation plane that can be similarly described by an axial rotational displacement field yI.The axial displacement y has units of length and is directed along the axis of rotation, with its magnitude being the effective radius at which the frame-invariant reference momentum density of the medium, ρc, would produce the angular momentum density | J|.
Expanding the bivector representation into 3D is illuminating.The derivative of a bivector has the same form as in Eq. ( 26) [29], yielding the expansions, and thus the acoustic field correspondence, These expressions confirm that the polar 3-vector x indeed acts as the usual mean displacement field of the medium away from its isotropic equilibrium state in the rest frame.This mean displacement field x has been used productively for numerical methods in acoustics [58,59,60,61] and plays a key role in microscopic derivations of the equations of motion [16], but its critical role as part of the acoustic angular momentum density potential M seems underappreciated.The time derivative of this displacement contributes one factor of the mean velocity field v of the medium, as expected.
The other two factors of the velocity come from curl of the rotational displacement field y.Since y has units of length, its curl ∇ × y is a unitless fraction of the speed c that contributes to the net linear velocity v. Indeed, note that if (ρc) yI → r × p y I were interpreted explicitly as an orbital angular momentum density with some constant linear momentum p y , then its curl, −c ∇ × y = −( ∇ × ( r × p y )/ρ = 2 p y /ρ = 2 v, would precisely yield twice the velocity.
Similarly to the electromagnetic vector potentials, there is also a gauge freedom in the representation p = ∇ • M/3, since ∇ • (∇ • (bI)) = 0 for any pseudo-4-vector field bI.Adding such a field bI = (ρc)γ 0 (b 0 /c− b)I shifts the bivector potential according to M → M + ∇ • (bI), analogously to the gauge transformations a e → a e + ∇χ e and a m I → a m I + ∇ • (χ m I), and expands to the 3D gauge transformation expressions, Moreover, imposing the condition ∇ ∧ M = 0 yields the simpler form p = −∇M/3 and the 3D constraints, which are the equivalents of the Lorenz-FitzGerald gauge conditions in Eq. ( 30) and simplify the acoustic equation to a causal wave equation for the bivector potential, ∇p = −∇ 2 M/3 = 0 [21].Now consider adding the geometrically motivated types of source to the acoustic equation of motion, as done previously for electromagnetism.The acoustic field p is an energy-momentum density 4-vector, so can have both scalar and bivector sources, Similar to electromagnetism, the scalar potential φ couples only to the scalar source ν, while the bivector potential M couples only to the bivector source N as component-wise wave equations.
Writing the measurable field equations with sources ν = ρ and N = F /c + ρ ΩI in 3D yields, These equations can be immediately interpreted to identify the physical meaning of each type of source: ρ is a mass-density scalar source that directly affects the local particle concentration of the medium, and thus the pressure as an energy density; F is a force density polar vector source, e.g., a directional speaker, that directly displaces the velocity field to produce boost angular momentum; and, Ω is a vorticity density axial vector source, e.g. a spinning propeller, that directly causes local rotations in the velocity field to produce rotational angular momentum.In summary, the full Hodge decomposition of the acoustic energy and momentum density fields in terms of the three possible potentials is, Neglecting the homogeneous background potentials P ′ 0 and v ′ 0 , the sources couple directly to each independent component of the three dynamical potential fields,

Completing the geometric descriptions
To this point the discussion has been primarily guided by the traditional representations of both acoustics and electromagnetism, including the measurable fields and their equations of motion.However, given the natural embedding of both theories into the geometric constraints imposed by the symmetries of spacetime, we can explore whether there is any additional freedom in each theory that has been traditionally overlooked but is physically meaningful.Examining the potential decompositions in the preceding section for each theory motivates their completion to fully use the available spacetime geometry.For convenience, we summarize our main results in the Tables 3 and 4.

Electromagnetic spinor field
Observe in the electromagnetic case that the potential representation in Eq. ( 36) motivates augmenting the description of F with two scalar fields that together form an auxiliary complex scalar field that keeps track of deviations away from the Lorenz-FitzGerald gauge constraints in Eqs.(30) and (42).That is, given the complex vector potential z em = λ − a e + λ + a m I in Eq. ( 49), a full vector derivative yields parts of three distinct grades, Bundled together with F , the three fields form an evengraded electromagnetic spinor field ψ em as a complete description of the measurable electromagnetic field.
Examining the form of these equations in 3D and recalling that the potential a e describes the field energymomentum per electric charge, the scalar field W e can be identified as a power per electric charge, This power vanishes as a causal continuity condition in the absence of external sources, giving a physical motivation for imposing the Lorenz-FitzGerald gauge constraint in Eq. (30).In principle, a nonzero W e would provide a mechanism to keep track of energy fluctuations in the field in the presence of sources, or when noncausal gauges are chosen.Similarly, the pseudoscalar field W m I can be identified as a power per magnetic charge, with an interpretation similar to the electric power potential.This field also vanishes as a continuity condition in absence of sources, giving a physical motivation to imposing the magnetic version of the Lorenz-FitzGerald gauge constraint in Eq. ( 42).
The complete spinor field ψ em thus describes both the electromagnetic power and force fields per electric and magnetic charge.That is, in the presence of a charged probe, these fields will determine the flow of energy and momentum between the probe and the electromagnetic field.Since this flow of energy and momentum to a probe can be monitored, the entire spinor field ψ em is a measurable physical field for electromagnetism, which implies that the power fields W e and W m must be tightly constrained by existing experimental data.
To see what these experimental constraints imply theoretically, we first show how the inclusion of non-zero power fields would modify the traditional theory, then show why existing experiments tightly constrain these fields so that they do in fact vanish in agreement with the Maxwell theory.
An important consequence of including the power fields W e and W m in the theory, even when they vanish, is that they restrict the potential gauge freedom of the electromagnetic field.That is, some potential gauge freedoms for the Faraday bivector F considered in isolation are not in fact gauge freedoms for the full spinor field ψ em .More precisely, the transformations a e → a e + ∇χ e and a m I → a m I + ∇χ m I, do not affect F , but they do affect the power fields, W e → W e + λ − c 2 ∇ 2 χ e and W m → W m + λ + c∇ 2 χ m I, unless the gauge generators additionally satisfy vacuum wave equations ∇ 2 χ e = ∇ 2 χ m = 0. Thus, even vanishing power fields imply nontrivial additional restrictions for the allowed gauge transformations of the theory compared to the traditional Maxwell theory.
Importantly, the potential equation of motion, yields the same causal wave equation as in Eq. ( 48), so this extension to the electromagnetic theory does not change the dynamics of potential fields satisfying a causal gauge.However, the inclusion of non-zero power fields does change the force field equations, making them a non-trivial extension to the Maxwell theory for potential gauges that do not satisfy the causal Lorenz-FitzGerald conditions.Expanding to a particular 3D frame yields the power-corrected Maxwell equations, That these corrections have been so successfully neglected in practice is strong experimental evidence that the Lorentz-FitzGerald causal gauge conditions, W e = W m = 0, are actually required physically and are not optional partial gauge constraints.
The standard electromagnetic Lagrangian density also gets corrections from the power fields, However, this standard Lagrangian density does not respect all symmetries of the full complex vector potential z em , so requires more fundamental correction.An alternative Lagrangian density that is properly invariant under the dual-symmetric gauge symmetry, z em → z em e Iθ , of its vacuum equation of motion ∇ 2 z em = 0, is [21], Note that this alternative is properly a pseudoscalar density that can be integrated over a volume as part of a 4-form, as discussed in Appendix A.4. Critically, this alteration involves both the original spinor field ψ em = ∇z em = λ − ∇a e + λ + ∇a m I and a dual spinor field, where ζ = µc is the wave impedance.This dual field exactly exchanges the roles of the electric and magnetic potentials and is a generalization of the Maxwell bivector, in Eq. ( 28).Importantly, this dual field ψ em,dual vanishes when the electric, ∇a e , and magnetic, ∇a m I, contributions to the total field ψ em are equal, unlike the Maxwell bivector G.The dual-symmetric Lagrangian L ems thus vanishes on shell [21], unlike the traditional Lagrangian, which is precisely what one should expect if the electromagnetic field were defined analogously to an acoustic field, i.e., as an effective mean field that averages over the microscopic dynamics of massless point particles with null energy-momenta.The full implications of this alternative Lagrangian will be explored in forthcoming work.
Returning to the Lagrangian-independent part of electromagnetic theory, the inclusion of nonzero power fields would also correct the measurable energy and momentum of the field, and thus the Lorentz force on probe charges.Specifically, the symmetrized (Belinfante) energy-momentum tensor has an intuitive bilinear form, satisfying −c T em (1) 0 = L em .Notably, this has precisely the same form as the spinor automorphism seen in Dirac's theory of the relativistic quantum electron.Picking a timelike flux direction γ 0 , this tensor yields,

Measurable fields
Force equation Table 3 Summary of electromagnetic structure, including all geometrically permitted fields and sources.The spinor field ψem augments the Faraday bivector F with a complex scalar field We + WmI that tracks measurable electromagnetic power per unit charge and add nontrivial corrections to the standard theory, while the electric and magnetic vector potentials compose a complex vector potential zem.The constants λ ± → 1/2 become symmetrized in vacuum, but are determined more generally by the coupling to sources [29].The standard Lagrangian density is shown along with a dual-symmetric alternative explored in Ref. [21], which introduces a dual spinor field ψ em,dual that generalizes the Maxwell bivector G and vanishes on shell.
resulting in frame-dependent energy and momentum, with the energy modified by the power fields.
The energy-momentum continuity condition for a probe along an arbitrary direction b is then [21], where p p = mu = mγ v (c + v)γ 0 is the probe particle momentum with 4-velocity u satisfying u 2 = c 2 , τ is its proper time, and is the Lorentz factor relative to a frame moving at velocity v.After using ∇ψ em = µ j p , this condition yields the modified Lorentz force, Here j p = qu = qγ v (c + v)γ 0 is the charge-current of the probe with a complex charge q = µ(q e − q m I/c), that includes both electric and magnetic parts in general [21].
Expanding this in a particular reference frame yields modified power and force equations for the probe, where, F em = q e E + (q e v) × (µ H) are the observed power and force on the probe in the frame γ 0 .Here the roles of W e and W m as transferable power per unit charge are particularly clear.These relations are summarized for convenience in Table 3.
To reiterate, these modified force and power equations imply that any non-zero power fields W e and W m should be experimentally measurable by weak probe charges just like the usual electromagnetic force fields E and H.The extra force terms that depend on the power fields are parallel to the probe velocity, so should be qualitatively distinguishable from the magnetic forces that are perpendicular to the velocity and the electric forces that are independent of velocity.These extra force terms would cause a probe charge to brake or accelerate along its direction of motion.Evidently, such an effect has not been necessary to include in order to describe experimental observations, implying that the power fields vanish physically, which is no surprise given their omission from the traditional Maxwell theory.
Importantly, however, the absence of these technically possible force terms has an added significance in this extended formulation, since verifying that the power fields are zero is an experimental confirmation of the Lorenz-Fitzgerald causal gauge conditions ∇ • a e = ∇ ∧ a m = 0 that set W e = W m = 0.In the traditional Maxwell theory, these conditions are optional partial gauge constraints; however, in this extended formulation the causal gauge constraints are an experimentally motivated necessity.Put another way, this geometrically complete electromagnetic theory is more tightly constrained by existing experimental data than the traditional electromagnetic theory.
As an intriguing side note, the possibility of consistently having a nonzero power-induced self-braking force term in the theory could be related to the longstanding problem of radiation back-reaction on a charged particle from its own emitted field.Including such selfinteraction in the traditional Maxwell theory yields a problematic self-acceleration term that seems to contradict experiment [62], and is notably derived by integrating the radiated power.It may be the case that including the nonzero power fields W e and W m near the charged particle, along with their corresponding force term parallel to the velocity, could help resolve this consistency issue for the traditional electromagnetic theory, which is an intriguing question for future investigation.

Acoustic complex 4-vector field
Similarly to electromagnetism, the acoustic potential representation in Eq. ( 50) motivates a geometric completion by augmenting the 4-vector field p with a pseudo-4-vector field wI.This addition corresponds to the unused bivector potential contribution ∇∧M and becomes part of a complex 4-vector field, that is analogous to the electromagnetic spinor field ψ em in Eq. (66).The Hodge decomposition of the introduced field wI then motivates adding two more potentials, w 0 I and φ w I, to yield a geometrically complete representation of z ac , with five acoustic potentials.However, much like the electromagnetic potential F 0 , the two potentials p ′ 0 and w ′ 0 must obey homogeneous equations and are not dynamical variables in a Lagrangian description, leaving only the three evengraded potentials (φ, M/3, φ w ) as dynamically fundamental, similarly to how the odd-graded potentials (a e , a m ) are fundamental in the electromagnetic case.Keeping only the dynamical potentials yields a spinor potential representation for the acoustic complex 4-vector field, Taken together, the five fields in z ac and ψ ac fully span all five geometric grades of the acoustic spacetime, just like the five fields in z em and ψ em for electromagnetism span the grades of normal spacetime.The two theories are thus grade-complements of one another.
To understand the meaning of wI, consider its relation to the potential fields, The contribution of M to this field has the form in Eq. ( 54), which immediately clarifies the meaning of the components of wI = γ 0 (P w /c − ρ w)I, Since yI is a rotational displacement field, wI evidently describes a rotational velocity field.That is, the axial vector w is a characteristic speed of rotational deformations directed along the axis of rotation.It is also notably characterized by the rotation vector ∇ × x that is used in geophysics to describe S-wave propagation from earthquakes [63].The quantity (ρc) wI can also be understood as a torque density for the medium, much like (ρc) v is a directed work density.The pseudoscalar pressure P w correspondingly describes a rotational energy density for the medium, much like the pressure P is an energy density for linear motion (work).
The pseudoscalar potential φ w I, in turn, evidently acts as a pseudoscalar action density.That is, the pseudovector (rotational) momentum density ρ wI = ∇(φ w I) is defined via a gradient while the pseudoscalar (rotational) energy density P w I = −∂ t (φ w I) is defined via a negative time derivative, in the same way that the vector momentum and scalar energy are defined by derivatives of a scalar action density in Hamilton-Jacobi theory.Incidentally, this natural appearance of both scalar and pseudoscalar action densities in acoustics gives strong motivation for considering complex scalar actions more generally.Indeed, such a complex action also naturally appears in the Hamiltonian-Jacobi approach to quantum theory [64], which is an interesting connection that will be explored in future work.
Like with the electromagnetic case, an important corollary of this geometric completion is that the gauge freedom of the theory is more tightly constrained.More specifically, given p = −∇ • M/3, there is a gauge freedom M → M + ∇ • (bI) that leaves p unaffected.However, this transformation does affect wI, which contains the contribution, unless the transformation generator bI also satisfies ∇ • (∇ ∧ b) = 0.This replacement of some gauge freedom with the explicit inclusion of previously neglected information using additional fields is entirely analogous to the electric and magnetic external power densities W e and W m that appear in the extended spinor description ψ em for electromagnetism in Eq. ( 69).The geometrically complete acoustic equation that includes all five grades of spacetime is, and accommodates a previously neglected pseudoscalar source ν w I.The primary potentials couple directly to these distinct types of sources, yielding independent component-wise wave equations.
In terms of the measurable energy-momentum fields, this source correspondence becomes, Expanding into 3D fields, using ν = ρ, N = F /c + ρ ΩI, and ν w I = ρw I yields the equations of motion, which are the acoustic equivalents of the extended 3D Maxwell Eqs.( 72)-(75) for electromagnetism.Notably, the physical meaning of these equations can be unambiguously interpreted as power and force equations.The traditional acoustic Lagrangian similarly acquires corrections from the pseudovector field wI, However, much like in the electromagnetic case, the vacuum acoustic equation ∇z ac = −∇ 2 ψ ac = 0 now supports an additional gauge symmetry, ψ ac → ψ ac e Iθ that is not respected by this scalar Lagrangian.To address this, we proposed an alternative (pseudoscalar) Lagrangian that does respect this symmetry in Ref. [29], .
As in the electromagnetic case, this involves both the vector field z ac = −∇ψ ac and a dual vector field, where the momentum density ζ ac = ρc = ρ/β is the acoustic version of the electromagnetic wave impedance.When the contributions of the complex scalar potential and the bivector potential to the total vector field z ac are equal, this dual field vanishes on shell, just as in the electromagnetic case [29].The Lagrangian correspondingly vanishes, as expected for massless phonons with null energy-momenta composing the acoustic field.The full implications of this alternative acoustic Lagrangian will be explored in forthcoming work.
Returning to Lagrangian-independent aspects of the acoustic theory, and leveraging the analogies between acoustics and electromagnetism, we can observe that the symmetrized (Belinfante) acoustic energy-momentum tensor should have a bilinear form analogous to Eq. (79), involving an automorphism by the complete measurable field z ac and satisfying −c T ac (1) 0 = L ac .Computing the energy-momentum density from this symmetrized tensor indeed yields, with the expected energy and momentum densities, that include contributions from both linear motion in p and rotational motion in wI.
The energy-momentum continuity condition for a probe particle then yields an acoustic force analogous to the Lorentz force.After using −∇z ac = ψ N = ν + N + ν w I for sources associated with the probe, is the net force felt by the probe from an acoustic field.
Expanding this relation into 3D using ν = ρp , N = F p /c + (ρ Ω) p I, and ν w I = ρw,p I yields, Thus, up to a density factor coupling the probe to the medium, the force applied to the probe is opposite the force it applies to the medium, as expected.Using more specific expressions relating these applied forces to the motion of the probe will then yield the appropriate probe equation of motion, in complete analogy to the electromagnetic Lorentz force prescription.These relations are summarized for convenience in Table 4.

Wave solutions
We now give explicit examples of wave solutions far from sources, to show the essential interplay between the various potential and measurable fields in the complete geometric arena of spacetime.These explicit solutions give strong support for the complete geometric picture being physically motivated and necessary, while also clarifying the role of the pseudoscalar phase freedom of the geometrically complete fields.

Potential fields
Force equation Summary of acoustics structure, including all geometrically permitted fields and sources.In the spinor potential ψac, the scalars φ + φwI form a complex action density, while the bivector M is an angular momentum density.The field zac includes an energy-momentum density 4-vector p dual to a rotational energy-momentum pseudovector wI.The proportionality constants λ ± = λ 4 → 1/2 symmetrize in vacuum, but are determined more generally by source coupling [29].The standard Lagrangian density is shown along with a generalization of the dual-symmetric alternative explored in Ref. [29], which involves the dual field zac that vanishes on shell.In the force equation, the probe interacts with the medium as an effective source ψ N that determines the coupling and thus the back-reaction on the probe.

Complex exponential waves
Wave solutions will rely upon an interesting property of spacetime that has been used implicitly up to this point: complementary pairs of geometric objects can be combined into complex objects that support intrinsic phases.The behavior of these phases is quite subtle for non-scalar grades and warrants a more careful discussion here.In particular, non-scalar grades admit null objects that have a global phase freedom that only becomes constrained locally by the equations of motion in the form of propagating waves.
A general spacetime multivector M can be factored, = (α + βI) + (a + bI) + F, where α, β are real scalars, a, b are 4-vectors, F is a bivector, and I is the spacetime pseudoscalar.The oddgraded part, z, can thus be understood as a complex 4-vector, while the even-graded part, ψ = ζ + F , can be understood as a spinor decomposed into a complex scalar ζ and bivector F .That these even and odd grades naturally combine in these ways has been essential for the preceding discussion of potentials and measurable fields, which underscores their physical importance.The complex scalar part behaves exactly as a standard complex scalar, so it can be written in polar form, with an intrinsic magnitude and phase, A conjugate can then be defined, ζ * ≡ |ζ| e −Iϕ ζ , by flipping the sign of the phase so that ζ * ζ = |ζ| 2 is a pure scalar with no pseudoscalar part.At least for scalars, this flip in sign of the phase is equivalent to flipping the sign of I in the Cartesian form as a sum of the two grades.This linear property of the conjugate will fail, however, for other grades.The complex vector part behaves analogously; however, I anti-commutes with vectors, so z = a + bI = a− Ib.The reverse z = a+ Ib = a− bI then seems like a conjugate (as we erroneously thought in Ref. [21]); however, it contracts with z to produce a complex scalar in two distinct ways, Each scalar is reversion invariant, ( zz) ∼ = zz, and can be written in a polar form with the same squared magnitude It follows that z should have a polar decomposition much like a complex scalar, with some canonical part z 0 satisfying, z 0 z 0 = z 0 z 0 = |z| 2 .Interestingly, however, squaring z yields, which is phase-invariant but with a bivector part, (z 2 − z 2 )/2 = 2(a ∧ b)I.The structure of z 0 is thus not generally a single-grade of vector, as one might expect.Instead, a zero phase implies that a b = 0, so the two halves of the remaining factor z 0 must not share any parallel overlap.This canonical vector can always be computed explicitly by inverting the phase of z: z 0 = z e −Iϕz/2 .A bivector F behaves analogously to a complex vector z, except I commutes with bivectors and the reversion F = −F is a simple negation.Consider a 3D expansion of a bivector F = A + BI for simplicity.It then follows that the square, is a complex scalar similar to zz, so can be written similarly in terms of a canonical part F 0 and a phase tan . Similarly to z, the canonical part can be computed directly by removing the phase, F 0 = F exp(−Iϕ F /2). Unlike for z, however, the canonical part F 0 always squares either to zero or a pure scalar by construction.
Combining the three independently complex grade sectors, a general multivector thus decomposes, into a sum of distinct polar decompositions.Common phases can be factored out as an overall global phase for the multivector.The relative phases between each sector of M then determine the canonical part of the total multivector.Critically, both complex vectors z and bivectors F can have null factors as their canonical parts, making their intrinsic phases degenerate.That is, if zz = 0 then z = z 0 e Iϕz/2 with null z 0 satisfying z 0 z 0 = 0 and arbitrary ϕ z .Similarly, if F 2 = 0 then F = F 0 e IϕF /2 with null F 0 satisfying F 2 0 = 0 and arbitrary ϕ F .It is precisely this global phase freedom of null fields that permits wave solutions for first-order vacuum equations.
Their phases are initially unconstrained until they are related locally to phases at nearby points via the equations of motion and then fixed by boundary conditions.Indeed, consider a simple spacetime-dependent phase like ϕ = k • r for some wave vector k = (ω/c + k)γ 0 and displacement r = (ct+ r)γ 0 .A complex scalar field with a constant amplitude and exponential phase, contains traveling waves that oscillate with constant overall amplitude α between complementary geometric grades as part of a rotation in an effective complex phase plane.A derivative of this phase factor produces, and cannot vanish unless k = 0, which prevents wave behavior.However, a second derivative, shows that ζ can satisfy the source-free wave equation if k 2 = 0.That is, the wave vector k must be null, and determined only by its frequency ω/c and a spatial unit vector k satisfying k2 = 1 in a particular frame γ 0 .Boosting to a different frame can scale the frequency ω/c (as a Doppler effect) and rotate the unit vector k, but not affect the nullity k 2 = 0.It then follows from, that we can identify the source-free dispersion relation ω = c| k| and wave vector k = (ω/c) k in frame γ 0 .A consequence of this necessary form for a complex scalar solution to a source-free wave equation is that its first derivative is a null complex vector, since zz = −α 2 e −I(k•r) k 2 e −I(k•r) = 0.This null factor enables the first-order source-free equation, ∇z = 0, to have an nontrivial wave solution.
Unlike the scalar imaginary i used in traditional ad hoc introductions of complex exponentials to describe wave behavior, the spacetime pseudoscalar I acts as a Hodge star operation to relate two geometrically complementary grades, making both "real" and "imaginary" parts of a complex exponential physically significant.The following sections explore the consequences of this natural complex exponential wave structure of spacetime more concretely in both electromagnetism and acoustics, which in turn clarifies why it is essential to include fields of all spacetime grades.

Electromagnetic waves
Consider the electromagnetic equation far from sources, ∇ψ em = ∇ 2 z em = 0. Since the quadratic wave operator acts on the complex vector potential, any wave solution with constant amplitude must have the form, with an intrinsic phase ϕ 0 when r = 0 and a null wave vector ±k = ±(ω/c)[1 + k]γ 0 that has two choices of sign for its characteristic frequency ±ω.Its canonical vector potential, then has zero phase by construction, meaning z 0 z 0 = z 0 z 0 , and thus a e,0 • a m,0 = 0 according to Eq. ( 117).
Since the phase factor causes rotations between the two parts of the canonical potentials, their contributions to the total field should symmetrize, which forces The electromagnetic spinor field then has the form, with constant canonical spinor, so the second derivative yields, with the prefactor vanishing as required, Notably, the scalar power fields in ζ 0 will vanish only when the vector potentials are transverse to the wave propagation, a e,0 • k = a m,0 • k = 0, making the triplet, (k, a e,0 , a m,0 ), an orthogonal set.This transversality is preserved by the phase factor, which only rotates between two vector potentials without altering their directionality.However, these solutions also admit longitudinal waves where the power fields in ζ 0 do not vanish.This possibility of longitudinal waves is neglected in traditional treatments of source-free electromagnetism that omit the possibility of scalar power fields.
To understanding this point more clearly in a particular frame, consider the terms involving the electric potential, a e,0 = (φ e,0 /c + A e,0 )γ 0 , In the transverse case with k • a e,0 = 0, it follows that φ e,0 /c = k • A e,0 , so = F e,0 = E e,0 /c + µ H e,0 I, is orthogonal to the propagation direction k.The polar electric field E e,0 is the frequency-scaled component of A e,0 orthogonal to k, while the rotational magnetic field H e,0 = µc k × E e,0 is an axial vector orthogonal to both k and A e,0 .The phase evolution of the wave then rotates these electric and magnetic fields around the propagation axis, naturally producing a circularly polarized electromagnetic wave [39,21], = with the direction of helical circulation of the fields around the propagation axis determined by the sign of the frequency ±ω in the wave vector k.That the two circularly polarized helicities of an electromagnetic wave can be expressed so naturally as complex exponentials with signed frequencies is encouraging.
In the case that a longitudinal part of the wave is admitted, the power field no longer vanishes, That is, solving for φ e,0 and simplifying Eq. ( 138) shows that such a nonzero power field affects the electric field, by restoring a longitudinal part parallel to k.
The complex structure of both the spinor field ψ em and the vector potential z em is critical for the treatment of electromagnetic waves using only invariant spacetime objects.Since the wave propagation itself exchanges the roles of the geometrically dual fields, all grades of field must be considered together on equal footing.The intrinsic phase-symmetry of the theory (i.e., electricmagnetic dual-symmetry) cannot be neglected, which in turn justifies correcting the fundamental Lagrangian density of the theory to the dual-symmetric form in Eq. (77) that correctly respects the global phase freedom of vacuum null field propagation.

Acoustic waves
The treatment of acoustic waves is equally enlightening.Consider the acoustic equation far from sources, ∇z ac = −∇ 2 ψ ac = 0.The complex vector field vanishes a single derivative, which implies that its constant canonical part z 0,ac should be annihilated by the null wave vector k = (ω/c)(1 + k)γ 0 for the wave motion, It follows that the canonical part has the form, with a pressure amplitude P for the longitudinal part of the wave aligned with the wave direction k and a pressure amplitude P⊥ for the transverse part of the wave aligned with the direction â such that k • â = 0. Noting that −kâ = −(ω/c)(1 + k)γ 0 â = (ω/c)(â + k × âI)γ 0 , the transverse wave helically rotates in the spatial plane spanned by the reference unit vector â and its orthogonal complement k × â, in close analogy to the transverse electromagnetic waves in Eq. ( 140).After factoring out a common reference null fourvector energy-momentum density for the wave, using real proportions σ ,⊥ such that, the constant wave amplitude has the laconic form, with the relative amplitudes, phases, and reference transverse direction relegated to a spinor transformation factor scaling p.This spinor can be understood as a natural extension of the Jones polarization vector for transverse electromagnetic waves, so characterizes both longitudinal and transverse parts of an acoustic wave.
The associated potential spinor is related with a derivative, z ac (r) = −∇ψ ac (r).Assuming a wave solution ψ ac (r) = ψ 0 exp(∓I(k • r)) with constant spinor amplitude ψ 0 thus yields the condition, For purely longitudinal waves with σ = 0 and σ ⊥ = 0, this constraint has two natural solutions.First, ψ 0, can be a complex scalar, Second, ψ ac, can be a pure bivector potential, after noting the eigenvalue relation that is proportional to the complex scalar action potential but explicitly aligned with the direction of wave propagation k.For purely transverse waves with c ⊥ = 1 and c = 0, on the other hand, there is no scalar potential solution but there is a pure bivector potential, that is also proportional to the form of the scalar action potential but instead aligned with a transverse direction â orthogonal to k.The scalar potential is not sufficient in isolation to include this important directional information.A general wave will be a linear combination of these longitudinal and transverse contributions.Notably, however, there is also a nontrivial internal gauge freedom in the spinor potential that doesn't contribute to the measurable energy-momentum field z 0,ac .Noting that k 2 = 0, satisfies kψ 0,g = 0 for any constant 4-vectors r n and r s with units of length.This degree of freedom is an irreducible spinor so has both a complex scalar part, p • r n + (p • r s )I and a corresponding bivector part, p ∧ r n + (p ∧ r s )I, that only together satisfy kψ 0,g = 0.This extra freedom is notable because the canonical bivector part matches the anticipated form of an angular momentum, but as an r-independent intrinsic (spin) part that does not contribute directly to the pressure and velocity fields.The complete potential representation of a constant amplitude acoustic wave thus has the form, where the λ s and λ b factors determine the representation proportions for the longitudinal part of the wave.For simplicity, let us focus on the bivector-biased representation with λ s = 0 and λ b = 1 and fix the relative phases so ∓Ie Iϕ 0, = ∓Ie Iϕ 0,⊥ = 1 in order to examine the basic structure of the bivector magnitude as a natural angular momentum density.With these simplifications, the bivector part has the form, Assuming the suggestive forms r n = (cτ n + r n )γ 0 and r s = (cτ s + r s )γ 0 , the gauge terms become, Thus, after adding the non-gauge terms, the displacement fields in the potential ψ 0 2 = (ρc)[ x 0 + y 0 I] are, These linear and rotational displacements determine the corresponding mass-moment N 0 = ρ x 0 and rotational angular momentum J 0 = (ρc) y 0 described by the bivector potential, and are precisely what one would expect as angular momentum given the momentum density ρ v 0 in the characteristic null energy-momentum density p for the wave augmented by an internal spinangular momentum density.
Finally, if we also consider non-constant amplitudes in the potential relation z ac = −∇ψ ac , then longitudinal waves z 0,ac, = pe Iϕ 0, also admit an r-dependent irreducible spinor potential, using ∇(kr) = −2k, since ∇(rk) = 4k and ∇(r ∧ k) = 3k.This r-dependent longitudinal spinor potential augments the scalar potential to include orbital angular momentum.That is, for ϕ 0, = 0, the bivector part of this potential has the form (p∧r)/3 = (c N + LI)/3 with mass-moment N = ρ v 0 t − ρ| v 0 | r and rotational angular momentum L = r × (ρ v 0 ) in the expected orbital forms.All five grades of spacetime play important roles in this analysis of acoustic waves.Similarly to the electromagnetic case, the phase evolution of the wave makes manifest the need for dual (linear-rotational) exchangesymmetry for the physical quantities in the theory.This observation in turn motivates correcting the acoustic Lagrangian to the dual-symmetric form in Eq. ( 104) that respects this phase freedom of null fields in vacuum, which will be further explored in future work.Most importantly, the potentials and measurable fields of every grade of spacetime have physically intuitive meanings.The acoustic spinor potential has a particularly rich structure that can accommodate both orbital and spin-angular momentum contributions.Which forms of the spinor potential are physically relevant will depend on its direct coupling to the sources and probes as boundary conditions.

Conclusion
In this paper we have provided a detailed overview of both acoustics and electromagnetism, using the same mathematical language of a spacetime Clifford bundle to highlight their many structural similarities and key differences.This convenient formalism naturally respects the strict geometric constraints of each theory, while also simplifying derivations and streamlining translations between frame-invariant and frame-dependent descriptions for clarity.Notably, acoustics and electromagnetism have a complementary grade structure, so taken together the two theories provide a comprehensive overview of how measurable relativistic fields can be represented by associated potential fields.
In the process of carefully exploring the complete formulations of each theory permitted by the constraints of spacetime, we highlighted a number of important generalizations and corrections to each theory.Some of these corrections are quite subtle, so have only become apparent by systematically following what is required by the geometry.These extensions also suggest several intriguing avenues of future research that may shed light on long-standing controversies and problems with each theory.We now briefly summarize a few of the notable observations and corrections.
The complete set of measurable fields for acoustics form an odd-graded complex 4-vector field, which includes a vector energy-momentum density combining the usual pressure and velocity fields, as well as a pseudovector rotational energy-momentum density.
The complete set of dynamical potential field of acoustics form an even-graded spinor field, which includes a complex scalar action potential density and an angular momentum potential density.Thus, all five grades of spacetime house meaningful field contributions that are physically relevant to the description of acoustic phenomena.Notably, the concreteness of acoustics as a mean field model for microscopic dynamics in the medium makes the interpretation of the derived fields straightforward, which thus provides a useful analogy to other relativistic field theories are traditionally less easily interpretable.
In particular, the dynamical gauge fields have clear physical interpretations at every stage of our development of acoustics, which directly refutes the commonly held belief that gauge fields are physically meaningless just because they have gauge freedoms.As a concrete example, one of the gauge fields for acoustics is a massdensity displacement vector field away from the equilibrium configuration of the medium, with its gauge freedom corresponding to the freedom of choice of origin for displacements.Just as introductory physics students are not told that spatial displacements are physically meaningless just because the choice of origin is arbitrary, we should also not tell ourselves that gauge fields are physically meaningless just because they have similarly arbitrary choices of reference.
The grade structure of electromagnetism is geometrically complementary to acoustics, making the comparison of the two theories illuminating.The complete set of dynamical potentials form an odd-graded complex vector field, which includes an electric vector potential describing energy-momentum per charge and a magnetic pseudovector potential describing rotational energy-momentum per charge.The complete set of measurable electromagnetic fields form an even-graded spinor field, which includes the expected Faraday bivector field describing force per charge, but may also include two additional scalar and pseudoscalar fields describing power per electric and magnetic charge, respectively.These added power fields directly keep track of any local deviations away from the causal Lorenz-FitzGerald gauge constraints.
The technical possibility of these extra power fields has been seemingly overlooked in standard treatments of electromagnetism, and their inclusion modifies the equations of motion as well as the Lorentz force felt by probe charges to make the electromagnetic force and power fields both experimentally testable.Importantly, the additional Lorentz force terms are parallel to the velocity, so are qualitatively different from the other forces and thus easily distinguishable as anomalous braking or self-acceleration effects, which implies that existing experimental data can already tightly constrain the values of the power fields.The clear absence of these power correction terms in existing lab data thus provides experimental support for the Lorenz-FitzGerald causal gauge conditions being satisfied in this geometrically complete treatment of electromagnetism, which contrasts with their interpretation as optional partial gauge constraints in traditional Maxwell electromagnetism.
As important examples, we analyzed in detail how wave solutions appear in each theory, without appealing to ad hoc complex exponentials using a scalar imaginary.Instead, we emphasize that wave solutions naturally involve the spacetime pseudoscalar as an effective imaginary unit, such that the waves produce phase rotations that continuously exchange information between pairs of complementary grades, in the sense of a Hodge-star grade-inversion duality.We highlighted that these intrinsic and frame-independent phase rotations are enabled by the phase-degeneracy of null bivector fields and null complex vector fields, which complement the usual complex scalar fields to give three distinct phase-rotation sectors of the spacetime geometry.Our detailed analysis of such wave propagation also shows why a geometrically complete description of each theory that includes all complementary grades is not just aesthetically pleasing, but physically necessary.
This global phase-freedom of fields in spacetime is an important gauge symmetry of each theory, which is underappreciated.In electromagnetism, this dual (phaserotation) symmetry exchanges the electric and magnetic sectors of the theory in a continuous way.Properly preserving this gauge symmetry in the Lagrangian densities requires nontrivial modifications that symmetrize the contributions of each potential to the measured fields.The predictions of physically conserved currents are determined by the symmetries of the Lagrangians, so these needed modifications have nontrivial experimental consequences.Indeed, recent measurements of local spin density in both electromagnetic and acoustic fields have been performed using small probe particles, and refute the spin angular momentum predictions of the standard dual-asymmetric Lagrangian densities as being experimentally incorrect.The common structure of the corrected dual-symmetric Lagrangians is intriguing, with each theory requiring the introduction of a dual field that vanishes on shell in vacuum.The dualsymmetric Lagrangian then also vanishes on shell, as should be expected for a mean-field average of massless (null) particles like phonons and photons.
As a final observation, the many important connections and clarifications that we presented throughout this work are a testament to the calculational efficiency, conceptual clarity, and comprehensive scope provided by the spacetime Clifford bundle formalism.Indeed, important concepts like the grade-rotating phase freedom of dual fields are difficult to see and express using standard mathematical frameworks for each theory, such as tensor component analysis or differential forms.We hope that this work will provide a useful reference for how to apply the Clifford bundle formalism to relativistic field theories more broadly, in both classical and quantum settings.
Acknowledgements The authors thank Konstantin Bliokh and Franco Nori for valuable discussions.This work was supported by NSF-BSF Grant Award No. 1915015.SA additionally thanks Chapman University for a Presidential Fellowship that helped to foster cross-disciplinary collaboration.

A The geometric algebra of spacetime
To discuss the geometric structure of both electromagnetism and acoustics, it is instructive to use a mathematical language of geometrically invariant objects that is coordinatefree.The Clifford algebra Cl 1,3 [R] of spacetime is a powerful choice of mathematical language suitable for this purpose [37,38,39,40,41,42,43,44,45,46,47,48,49,50].It is the largest associative algebra generated by unit vectors {γµ} 3 µ=0 over the real numbers R that satisfy the Minkowski spacetime metric η(γµ, γν ) = ηµν with signature (+, −, −, −).As such, it conveniently contains and unifies many other disparate mathematical frameworks while also clarifying subtle connections between them.This algebraic framework also makes manifest the deep connections between geometry and the treatment of spin in quantum mechanics, since the Dirac matrices are matrix representations of spacetime unit vectors and the Clifford product.A detailed and pedagogical exposition of this spacetime algebra as applied to electromagnetism can be found in Ref. [21].

A.1 Spacetime Clifford algebra
As a brief reminder, the Clifford product between spacetime vectors a = µ a µ γµ and b = ν b ν γν is associative, (ab)c = a(bc) and contains both the Minkowski metric and the Grassmann wedge product as its symmetric and antisymmetric parts, respectively, The full product ab is thus generally noncommutative, ab = ba, and is invertible for factors of nonzero magnitude, a −1 = a/a 2 .The metric a • b contracts the vectors to a scalar point by projecting b onto the direction of a, then projecting that product of lengths with shared direction onto the point at the origin of a. Conversely, the wedge product a ∧ b expands the vectors by dragging b along a and attaching its tail to the head of a, resulting in a plane segment with magnitude equal to the dragged parallelogram area and a sign that indicates a circulating orientation along a, then along b, then back along −a, then back along −b.
The full product ab thus spans grades and can be written as a rotor generalization of a complex number in polar form, that indicates that a rotates in the unit plane C by an angle θ to reach b, with the signature (C 2 ) of C dictating whether the rotation is hyperbolic or elliptic.The inverse of a product is (ab , where (ab) ∼ ≡ ba is the reversion operation that reverses (or transposes) the order of all Clifford products.Successive wedge products like a ∧ b ∧ c are associative (a ∧ b) ∧ c = a ∧ (b ∧ c) and create higher-dimensional volumes in spacetime with similarly defined circulating orientations.Successive Clifford products like abc thus generally have multi-grade content, but also have natural inverses (abc) −1 = (abc) ∼ /(a 2 b 2 c 2 ).For reference, Table 1 gives a compact summary of the full graded basis of 2 4 orthogonal unit elements that span the five distinct grades of spacetime, from points to 4-volumes.
The spacetime algebra formalism has the notable advantage of simplifying the translation between geometrically invariant 4D quantities and their relative 3D decompositions observed in particular reference frames.The spatial frame {γ k } 3 k=1 of each observer in spacetime is dragged along a particular temporal direction γ 0 , creating 3 space-time planes γ k γ 0 = γ k ∧ γ 0 ≡ σ k that are perceived as 3D spatial unit vectors by observers within the evolving frame.Each of the remaining 3 spatial planes γ i γ j ≡ −ǫ ijk (I σ k ) = − σ i σ j = − σ i ∧ σ j is orthogonal to a 3D axis σ k around which its orientation rotates.For reference, Table 2 gives a compact overview of the closed 3D Clifford subalgebra of a particular frame and highlights its formal equivalence to the algebra of biquaternions, as well as the complex Pauli matrix algebra for nonrelativistic spin.
This natural embedding of 3D reference frames within the algebra permits invariant 4-vectors to be easily factored and expressed as frame-dependent paravectors.For example, energy-momentum factors, into its relative scalar energy E/c and 3-vector momentum p = k p k σ k components.Similarly, the six components of invariant bivectors straightforwardly split into a complex representation involving relative polar and axial 3-vectors.For example, the electromagnetic field F = 1 2 3 µ,ν=0 F µν γµ ∧ γν is a bivector that splits, into a polar electric field E and axial magnetic field H.This formalism also handles group transformations in a natural way.For example, the Hodge duality between a spatial plane and an orthogonal axis in relative 3D space defines the Gibbs 3-vector cross product, that is closed for grade-1 3D vectors at the expense of breaking the associativity of the wedge product, but this definition is identical to the spin-(1/2) rotation group commutator relations for the Pauli spin-matrices, if the pseudoscalar I is identified with the scalar imaginary i when restricted to 3D space.More generally, the six unit planes of spacetime form a closed subalgebra under the wedge product, which is then precisely the Lie bracket for the Lorentz group of spacetime rotations generated by the unit planes directly.That is, the three space-time planes σ k ∼ −iK k with positive signature generate hyperbolic boost rotations, while the three purely spatial planes I σ k ∼ −iJ k with negative signature generate elliptical spatial rotations, exp(θ Thus, Eq. ( 165) is also equivalent to the Lorentz group Lie bracket relations [K i , K j ] = −iǫ ijk J k , with the other relations, [J i , J j ] = iǫ ijk J k and [J i , K j ] = iǫ ijk K k , obtained as variations with extra factors of the pseudoscalar I.It follows that any Lorentz transformation of an element A of the algebra has the form of a group inner automorphism with half-angle rotors generated by a unit spacetime plane B, where the rotor inverse is R(θ) −1 = R(θ) ∼ = R(−θ).In particular, rotations in purely spatial planes R(θ) = exp(I n θ/2) around a 3D spatial unit vector n = k n k σ k have precisely the same form used for Pauli spin-rotations.

A.2 Spacetime Clifford bundle
The preceding treatment holds for a particular spacetime volume containing all geometric content in the spacetime Clifford algebra.However, each point in spacetime has an infinitesimal volume around it that can contain local geometric objects at that point, forming a Clifford bundle over spacetime in which fields Φ(x) can be defined of arbitrary geometric grade.The dynamics of these spacetime fields both enact local geometric transformations at each point x and connect the geometric content of nearby points by parallel transporting elements of their tangent algebras along curves connecting those points.More precisely, given a flat spacetime manifold M, each point x ∈ M has a tangent space T M(x) spanned by a local tangent vector basis {γµ(x)} 3 µ=0 .Products of these tangent vectors construct a spacetime algebra Cl 1,3 [R](x) that is local to the point x.Each tangent space formally has a dual cotangent space of functions T * M(x) = T M(x) → R that is spanned by a basis of one-forms {ω µ (x)} 3 µ=0 defined from the tangent vectors γµ(x) to jointly satisfy the Euclidean metric ω µ (x)[γν (x)] = δ µ ν .The one-form basis is in one-toone correspondence with a reciprocal vector basis {γ µ (x) = (γµ(x)) −1 } 3 µ=0 in the tangent space, defined such that γ µ (x) • γν (x) = δ µ ν , according to the identity ω µ (x)[v(x)] ≡ γ µ (x)•v(x) for any tangent vector v(x) ∈ T M(x).This correspondence means that any one-form α(x)[v(x)] = µ αµω µ (x)[v(x)] = ( µ αµγ µ (x)) • v(x) ∈ T * M(x) can be identified with a tangent vector a(x) = µ αµγ µ (x) ∈ T M(x) with the same components αµ in the reciprocal basis.This corresponding vector also has the component expansion a(x) = µ a µ γµ(x), so γµ(x) • a(x) = aµ = a µ ηµµ and the two component expansions of the same geometric object are related by the metric as expected from tensor component analysis.In a similar way, the entire Grassmann algebra of forms in the cotangent space can be identified with geometric objects in the graded structure of the tangent space, which is a useful simplification.
These independent tangent Clifford algebras and associated cotangent spaces of forms at each point x must then be connected by linear translations between nearby points.The is the appropriate generator of translations along the spacetime manifold, potentially modified by adding an appropriate connection Dx = ∇x + A(x) that enacts the parallel transport of relevant geometric content during translation.Note the natural expansion of ∇x into reciprocal (dual) basis vectors, which will be clarified in the next section.Each ∂/∂x µ (x) is the partial derivative along a particular coordinate function x µ (x) evaluated at the point x.Also note that the point r = 3 µ=0 x µ (x)γµ(x) can itself be expressed in terms of these same coordinate functions as a 4-vector displacement r from the origin, assuming a flat manifold with parallelized tangent vector directions for simplicity, which conveniently lifts a vector representation of the manifold points into the same tangent Clifford algebra for ease of calculation.
It is customary to omit the point x for brevity in ∇x and write the partial derivatives more compactly as ∂µ, yielding the simpler form ∇ = µ γ µ ∂µ.Further noting that γ k = (γ k ) −1 = −γ k for k = 1, 2, 3 and γ 0 = γ 0 then yields convenient relative frame expansions of the (covariant) vector derivative, which differ in sign from similar expansions of (contravariant) 4-vectors like the displacement r = γ 0 (ct − x) = (ct + x)γ 0 in the tangent space.

A.3 Spacetime differential forms
Given a scalar field φ(x) on the manifold, the exterior derivative d creates a differential one-form, that takes a tangent vector v(x) and returns a linear approximation to a change in φ(x), in the form of a directional (partial) derivative along v(x).Note that the differential is entirely characterized by the total vector derivative ∇φ(x) of the field.This identification between the exterior derivative d acting in the cotangent space T * M and the vector derivative ∇ acting in the tangent space T M is convenient, and is why the partial derivative components in Eq. ( 170) are naturally expressed in the reciprocal basis.
Taking the exterior derivative d of a one-form field such as ω(x)[v(x)] = a(x) • v(x) similarly raises its grade, and can be identified with the 4-curl ∇∧a(x) of the object a(x) in the tangent space that corresponds to ω(x) in the cotangent space.This identification between the exterior derivative d and the 4-curl ∇∧ holds generally for differential forms of any grade in the cotangent bundle, including grade-0 scalar fields where ∇ ∧ φ(x) = ∇φ(x) is identified with the gradient.The exterior derivative satisfies the useful identity d 2 (•) = 0 (equivalent to a Bianchi identity), which is also satisfied by the corresponding four-curl, ∇ ∧ (∇ ∧ (•)) = 0. Given the metric and the Hodge star operation ⋆, the codifferential δ ≡ ⋆ −1 d⋆ can also be defined on differential forms, which lowers the grade by one rather than raising the grade like the exterior derivative.Using the correspondence between ⋆ in the cotangent space and I in the tangent space, this definition yields so is equivalent to the grade-lowering 4-divergence ∇ • .Moreover, the codifferential satisfies another useful (Bianchi) identity δ 2 (•) = ⋆ −1 d 2 ⋆ (•) = 0, which is also satisfied by the corresponding four-divergence, ∇ • (∇ • (•)) = 0.It then follows that the vector derivative under the Clifford product precisely unifies both the grade-raising exterior derivative and the grade-lowering codifferential of the cotangent space of forms, but is expressed in the tangent algebra, (δ + d)ω(x) ∼ ∇a(x) = ∇ • a(x) + ∇ ∧ a(x). (175) It then naturally follows that the square of the vector derivative gives the correct 4-Laplacian, or scalar wave operator, (δ + d) 2 ω(x) = (δd + dδ)ω(x) (176) By working with the vector derivative in the tangent Clifford bundle, all the convenience and clarity provided by differential forms is preserved, but without the complication of excess function arguments as is necessary with forms.This simplification enables much closer correspondence to standard notations and vector calculus methods, as used in the main text.In fact, one can argue that the unification of the exterior derivative and codifferential encourages more structurally transparent calculations than with differential forms alone.
For example, the Maxwell electromagnetic equations of motion can be written as the single Clifford algebra equation ∇F = j, where F is a bivector and j = je + (jm/c)I is a complex 4-vector containing both electric and magnetic sources.This contrasts with the pair of equations commonly seen in treatments that use differential forms, δF = µ je and δG = µ jm/c, where F is the Faraday 2-form and G = ⋆F is the dual Maxwell 2-form treated separately, with the currents je and jm both treated as one-forms despite the magnetic current jm properly having the character of a pseudo-4-vector geometrically.One can fix this standard differential forms treatment by writing instead, (δ +d)F = µ (je +⋆jm/c), which then becomes identical in the cotangent bundle to the content of the spacetime Clifford algebra formulation in the tangent bundle, ∇F = µ j, albeit using less compact notation.

A.4 Lagrangian densities
Recall that a Lagrangian density properly is a 4-volume integration measure that is a function of the dynamical fields; extremizing the action as the integrated Lagrangian density with respect to variations of the dynamical fields then produce the classical equations of motion.Thus, the Lagrangian density should geometrically correspond to a volume 4-form on spacetime that is invariant under the correct physical symmetries, including those of spacetime itself.
To make explicit how this idea connects to the tangent Clifford bundle used here, a 4-form is a local linear map d 4 S(x)[v 0 (x), v 1 (x), v 2 (x), v 3 (x)] = L(x)•(v 0 (x)∧v 1 (x)∧v 2 (x)∧ v 3 (x)) in the cotangent space at a point x that takes four 4vector arguments in the tangent space at x, constructs an invariant pseudoscalar 4-volume in the tangent space out of them, then contracts that volume with another local 4-volume L(x) to produce a scalar measured value.The 4-volume L(x) is a grade-4 object in the tangent Clifford algebra, so is proportional to the local unit pseudoscalar I(x) at the point x.A spacetime Riemann sum, k d 4 S(x)[δV k (x)], implicitly passes in small vector increments (δx) µ γµ(x) along independent vector directions to construct small invariant 4-volumes δV k (x) = ( µ (δx) µ )I(x) that are proportional to the local pseudoscalar I(x), then evaluates the 4-form measure d 4 S(x)[δV k (x)] on each such local volume, and sums those measured values over a grid of small volumes that partition the spacetime volume being integrated.Taking the limit as the small volumes δV k (x) become infinitesimal in magnitude yields a spacetime action integral S = d 4 S(x) = L(x) • d 4 V (x) as the limiting value of the Riemann sum, with resulting directed volume measure The important takeaway is that such a integration measure 4-form is fully characterized by a grade-4 pseudoscalar field L(x) in the tangent Clifford bundle, which we can identify as the geometric content of the Lagrangian density.For the action S to yield an invariant scalar value, as usually assumed, the density L(x) must be proportional to the unit 4-volume field I(x).The traditional scalar Lagrangians in Eqs.(23) and (32) leave out this volume factor as implicit.

4 -
vector derivative (or Dirac operator, familiar from the Dirac equation for the quantum electron),