The Rarita–Schwinger Equation: Spin 3/2 in Biquaternionic Form

Introduction

The Proca equation carries spin one, and the Dirac equation carries spin one-half, and both are first-order equations for a single field. Spin three-halves is the first case that requires a field with two indices: the Rarita–Schwinger field $\psi_\mu$ is a four-vector whose components are Dirac spinors, sixteen complex components in all, and its equation removes twelve of them, leaving the four states of a spin-$\tfrac32$ particle. The equation was written down by Rarita and Schwinger in 1941, and it is the field equation of the gravitino in supergravity; it is also the place where the standard field theory of higher spin first becomes delicate, since its minimal coupling to electromagnetism propagates superluminally.

This article writes the Rarita–Schwinger system in biquaternionic form and reduces it in the rest frame, where the spin-$\tfrac32$ content is transparent. The framework treats the field as a biquaternion-valued four-vector, that is, as a family $\tilde{\Psi}_\mu$, $\mu=0,1,2,3$, of elements of $\mathbb{B}$, and the equation as two statements: a Dirac equation on each component, and a trace constraint that contracts the vector index with the Clifford generators. In the rest frame the two together force the time component to vanish, force the spatial components to their upper (positive-frequency) parts, and impose the single vector condition $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=0$ on the spatial amplitudes. The surviving space is exactly four-dimensional, and it carries the spin-$\tfrac32$ quartet. The construction also exhibits the spin-$\tfrac32$ projector: the transverse projector on the vector index, $P_{ij} = \delta_{ij} - \tfrac13\sigma_i\sigma_j$, which is idempotent and selects precisely the constrained spatial amplitudes.

Two features of the standard theory are recorded and read in the algebra. First, in the massless case the system has a gauge freedom $\psi_\mu\to\psi_\mu+\partial_\mu\epsilon$, with $\epsilon$ an arbitrary Dirac spinor; the freedom removes the spin-$\tfrac12$ branches that the representation theory places in the field, leaving the two helicities $\pm\tfrac32$. Second, in the massive case the constraint structure is more rigid and the minimal coupling is known to be inconsistent: the field propagates acausally in an external electromagnetic background, the Velo–Zwanziger obstruction, and a consistent interacting theory requires either a gravitational background or a gauge principle. The framework represents these facts; it does not remove them.

The conventions are those of the companion articles. The algebra is $\mathbb{B}=\mathbb{C}\otimes_\mathbb{R}\mathbb{H}$ with basis $e_0=1,e_1,e_2,e_3$, $e_k^2=-e_0$, central scalar imaginary $i$, gradient $\tilde{\nabla}=e_0\partial_{ict}+\nabla$ with conjugate $\tilde{\nabla}^{\natural}$ and d'Alembertian $\Box=\tilde{\nabla}\tilde{\nabla}^{\natural}=\partial_{ict}^2+\Delta$. The Clifford generators are used as a translation tool with the mostly-minus metric $g=\mathrm{diag}(+1,-1,-1,-1)$, so that $(\gamma^0)^2=+I_4$, $(\gamma^k)^2=-I_4$, and $\{\gamma^\mu,\gamma^\nu\}=2g^{\mu\nu}I_4$, and the Dirac representation $\gamma^0=\mathrm{diag}(I_2,-I_2)$, $\gamma^k=\big(\begin{smallmatrix}0&\sigma^k\\-\sigma^k&0\end{smallmatrix}\big)$ is used for the explicit rest-frame computation. The chirality operator is $\gamma_5 = i_{\mathrm{Cl}}\gamma^0\gamma^1\gamma^2\gamma^3$, with $i_{\mathrm{Cl}}$ the scalar imaginary of the complexified Clifford algebra rather than the biquaternion imaginary, and the total antisymmetrizations are $\gamma^{\mu\nu}=\tfrac12[\gamma^\mu,\gamma^\nu]$ and $\gamma^{\mu\nu\rho}=\tfrac1{3!}\sum_{\mathrm{perms}}\pm\gamma^\mu\gamma^\nu\gamma^\rho$.

The Standard Rarita–Schwinger Field

The field and its equation

The Rarita–Schwinger field is a vector-spinor: a four-vector index and a Dirac spinor index, written $\psi_\mu{}^\alpha$ with $\mu=0,1,2,3$ and $\alpha=1,\dots,4$. In the massless case the equation is

$$ \epsilon^{\mu\nu\rho\sigma}\gamma_5\gamma_\nu\partial_\rho\psi_\sigma = 0, $$

equivalently $\gamma^{\mu\nu\rho}\partial_\rho\psi_\nu=0$, the constant between the two standard writings being the one carried by the $\epsilon$–$\gamma_5$ duality of the companion dictionary and being immaterial when $m=0$. In the massive case the equation is

$$ \epsilon^{\mu\nu\rho\sigma}\gamma_5\gamma_\nu\partial_\rho\psi_\sigma + m\,\gamma^{\mu\nu}\psi_\nu = 0, $$

equivalently, in the trivector form, $(i\gamma^{\mu\nu\rho}\partial_\rho - m\gamma^{\mu\nu})\psi_\nu=0$; the two writings are related by the duality identity

$$ \epsilon^{\mu\nu\rho\sigma}\gamma_5\gamma_\nu\partial_\rho\psi_\sigma = -i\,\gamma^{\mu\sigma\rho}\partial_\rho\psi_\sigma, $$

which the companion dictionary fixes once the chirality operator is $\gamma_5=i_{\mathrm{Cl}}\gamma^0\gamma^1\gamma^2\gamma^3$ and the antisymmetrized products are taken with unit weight. That identity also fixes the relative sign of the mass terms: the pair displayed above is the one in which the trivector form carries the same minus as the Dirac operator of the companion article, $(i\not\partial-m)\psi=0$, and it is the pair whose positive-frequency solutions have the rest-frame frequency $\omega=m$, as the reduction below shows.

The equation is the Euler–Lagrange equation of the Rarita–Schwinger Lagrangian, whose kinetic term is the same antisymmetrized combination and whose mass term is the antisymmetrized pair $\bar{\psi}_\mu\gamma^{\mu\nu}\psi_\nu$,

$$ \mathcal{L} = -\tfrac12\,\bar{\psi}_\mu\,\epsilon^{\mu\nu\rho\sigma}\gamma_5\gamma_\nu\partial_\rho\psi_\sigma - \tfrac12\,m\,\bar{\psi}_\mu\gamma^{\mu\nu}\psi_\nu , $$

up to the overall normalisation, which is a convention; varying with respect to $\bar{\psi}_\mu$ returns the equation displayed above. The index structure is fixed by the requirement that the equation be first order and that its free solutions carry spin $\tfrac32$. No first-order equation for a single lower-spin field can do that: the two-index structure is necessary, and the antisymmetrizations are what make the equation carry an irreducible spin-$\tfrac32$ rather than a reducible sum. The Lagrangian also shows where the mass term sits: it is the antisymmetrized pair $\gamma^{\mu\nu}$, not the Clifford metric $g^{\mu\nu}$, and this is what distinguishes the spin-$\tfrac32$ mass term from the Proca mass term $\tfrac12m^2A_\mu A^\mu$.

The constraint and the reduction

The massive equation implies the trace in two contractions, and the same two contractions deliver the reduced dynamics. Write $S := \gamma^\sigma\psi_\sigma$ for the trace spinor. Taking the divergence of the equation kills the kinetic term, because $\gamma^{\mu\sigma\rho}$ is antisymmetric in its first and third indices, and leaves, using $\gamma^{\mu\sigma} = \gamma^\mu\gamma^\sigma - g^{\mu\sigma}$,

$$ 0 = \gamma^{\mu\sigma}\partial_\mu\psi_\sigma = \not\partial\,S - \partial^\sigma\psi_\sigma . $$

Contracting the equation with $\gamma_\mu$ instead, and using the two traces $\gamma_\mu\gamma^{\mu\sigma\rho} = 2\gamma^{\sigma\rho}$ and $\gamma_\mu\gamma^{\mu\sigma} = 3\gamma^\sigma$, leaves

$$ 2i\,\gamma^{\sigma\rho}\partial_\sigma\psi_\rho = 3m\,S . $$

The left-hand side vanishes as well, again by $\gamma^{\sigma\rho}=\gamma^\sigma\gamma^\rho-g^{\sigma\rho}$ together with the relation just obtained, so for $m\neq0$ the trace vanishes,

$$ \gamma^\mu\psi_\mu = 0 , $$

and with it the divergence $\partial^\mu\psi_\mu = 0$ follows from the first relation. The trace is therefore a consequence of the massive equation and not an independent condition.

The trace in turn linearises the double product. Multiplying the identity $\gamma^\mu\gamma^\sigma = \gamma^{\mu\sigma}+g^{\mu\sigma}$ by $\psi_\sigma$ and summing over $\sigma$ gives $\gamma^{\mu\sigma}\psi_\sigma = -g^{\mu\mu}\psi_\mu$ when the trace vanishes, with no sum on $\mu$. Substituting this, and the vanishing of $\gamma^{\sigma\rho}\partial_\rho\psi_\sigma$ and of $\gamma^\sigma\partial^\mu\psi_\sigma = \partial^\mu S$, into the equation written through $\gamma^{\mu\sigma\rho} = \gamma^\mu\gamma^{\sigma\rho}-g^{\mu\sigma}\gamma^\rho+g^{\mu\rho}\gamma^\sigma$, leaves

$$ E_\mu = \bigl(i\not\partial - m\bigr)\psi_\mu \qquad\text{(no sum on }\mu\text{)} , $$

which vanishes exactly when the component Dirac equation holds,

$$ (i\not\partial - m)\psi_\mu = 0 \quad\text{for each }\mu . $$

The antisymmetrized equation is thus equivalent to the Dirac equation on each component together with the trace, and, through them, the divergence. The mass term regenerates the Dirac operator on each component, so that the reduced dynamics is the one the companion Dirac article supplies. The trace constraint removes one Dirac spinor's worth of components from the sixteen, and it is the part of the reduction that the algebra can state as a single contraction.

The massless case separates the constraints, and the separation is what the gauge freedom is for. For $m=0$ the same two contractions no longer force the trace: $\gamma^\mu\psi_\mu=0$ is not contained in the equation, and adjoining it lowers the solution space from six dimensions to four. The divergence is a different matter, because its four components are linear combinations of the components of the equation: $\partial^\mu\psi_\mu=0$ remains a consequence of the equation even when the mass vanishes. It is in the massless case that the gauge freedom

$$ \psi_\mu \longrightarrow \psi_\mu + \partial_\mu\epsilon, $$

appears, with $\epsilon$ an arbitrary Dirac spinor: the antisymmetrized kinetic structure is built from the field strength $\partial_\rho\psi_\sigma$, which changes by a symmetric derivative that the antisymmetrization discards. The gauge freedom is the spin-$\tfrac32$ analogue of the electromagnetic gauge freedom, it is what removes the spin-$\tfrac12$ branch that the vector index introduces, and it is also what makes the trace usable: the trace values that occur on the massless solution space are exactly the ones a gauge transformation can shift, so the gauge $\gamma^\mu\psi_\mu=0$ is attainable and the trace is a complete gauge choice rather than a derived constraint.

Representation content

The Lorentz transformation character of the field is read off from the two indices. The vector index carries $(\tfrac12,\tfrac12)$ and the spinor index carries $(\tfrac12,0)\oplus(0,\tfrac12)$, so the field transforms in the tensor product

$$ \left(\tfrac12,\tfrac12\right)\otimes\left[\left(\tfrac12,0\right)\oplus\left(0,\tfrac12\right)\right] = \left(1,\tfrac12\right)\oplus\left(0,\tfrac12\right)\oplus\left(\tfrac12,1\right)\oplus\left(\tfrac12,0\right), $$

of complex dimension $6+2+6+2=16$, as the component count requires. On restriction to the rotation subgroup, with $j\otimes j'$ evaluated at the diagonal,

$$ \left(1,\tfrac12\right)\big|_{SU(2)} = \tfrac32\oplus\tfrac12, \qquad \left(0,\tfrac12\right)\big|_{SU(2)} = \tfrac12, $$

and conjugate statements for the other two summands. The field therefore contains

$$ 2\times\tfrac32 \;\oplus\; 4\times\tfrac12 $$

under $SU(2)$, of dimensions $2\cdot4 + 4\cdot2 = 16$. Two spin-$\tfrac32$ multiplets — one for the particle and one for the antiparticle — and four spin-$\tfrac12$ multiplets. The whole content of the Rarita–Schwinger equation and of its gauge structure is the removal of those four spin-$\tfrac12$ multiplets, and the rest-frame computation below exhibits the removal explicitly.

Counting the physical states

Two independent conditions remove the spin-$\tfrac12$ content, and because they act on different structures their effects add. The algebraic trace $\gamma^\mu\psi_\mu=0$ removes one Dirac spinor's worth of components, which is two spin-$\tfrac12$ multiplets. The gauge freedom in the massless case removes another Dirac spinor's worth, another two spin-$\tfrac12$ multiplets. What is left is the two spin-$\tfrac32$ multiplets, of total dimension eight, and for each frequency sign that is the four-dimensional quartet of a spin-$\tfrac32$ particle. In the massive case the gauge freedom is absent, and the second pair of spin-$\tfrac12$ multiplets is removed instead by the subsidiary conditions that follow from the equation; the physical count is the same four states per charge.

The counting is therefore representation-theoretic and does not depend on the detailed form of the antisymmetrizations: any first-order equation whose solutions are pure spin-$\tfrac32$ must remove exactly the four spin-$\tfrac12$ multiplets from the tensor product, and the trace constraint and the gauge freedom are the two mechanisms that do it. This is the reason the field is the minimal carrier of spin $\tfrac32$: the next tensor product that could carry it, the three-index object, would carry still more lower-spin admixture, and the two-index object is the smallest that contains $\tfrac32$ at all.

The known difficulties

In the massive case the constraint structure is rigid: the components of $\gamma^\mu\psi_\mu$ and of $\gamma_\mu(\text{equation})$ form a chain of algebraic conditions with no gauge freedom, and the resulting theory has the right free-field content but is known to be inconsistent upon minimal coupling to an external electromagnetic field. The propagation becomes acausal for sufficiently strong backgrounds — the Velo–Zwanziger obstruction — and a consistent interacting spin-$\tfrac32$ theory requires a gravitational or gauge principle rather than a fixed background. This is the standard status of the field, and it is recorded here because the biquaternion formulation does not alter it: the algebra rewrites the free equation faithfully, and the difficulties belong to the coupling.

The Vector-Spinor Field in the Biquaternion Framework

The carrier

The biquaternion carrier of a spin-$\tfrac32$ field is a biquaternion-valued four-vector,

$$ \tilde{\Psi}_\mu \in \mathbb{B}, \qquad \mu=0,1,2,3, $$

that is, a family of four elements of the algebra, one for each spacetime index. In components, $\tilde{\Psi}_\mu = \Psi_{\mu,0}e_0 + \Psi_{\mu,1}e_1 + \Psi_{\mu,2}e_2 + \Psi_{\mu,3}e_3$ with complex coefficients, so that the carrier has $4\times4 = 16$ complex components, matching the standard field. The vector index is carried by the algebra's own four-vector structure, as it is for the Proca potential, and the spinor index by left multiplication on the algebra, as it is for the Dirac field, whose carrier is a minimal left ideal of $\mathbb{B}$. The Rarita–Schwinger field is therefore a mixed object: neither an element of a module over $\mathbb{B}$ nor a pure element of $\mathbb{B}$ transforming by conjugation, but a tensor product of the two. This is the structural reason it lies outside the algebra's native spin content, a point the final article of the subcategory develops.

The components of the equation

The standard equation is transcribed component by component, using the correspondence between the Clifford generators and the algebra's units. The correspondence is the dictionary's frame identity, in which the four generators are written in terms of the fixed odd frame $\gamma^0$ and four representatives carrying the directions,

$$ \gamma^\mu = \gamma^0\,\Phi(e^\mu), \qquad e^0 = e_0, \qquad e^k = i e_k , $$

with $\Phi$ the dictionary's map onto the even Clifford algebra; equivalently, a four-vector contracted with the generators is the fixed frame times the image of the corresponding element of the informational sector $\mathbb{M}_+$, whose basis is $\{e_0, ie_1, ie_2, ie_3\}$. Products of generators are governed by the dictionary's multiplicativity, $\Phi(\tilde{P}\tilde{Q}) = \Phi(\tilde{P})\Phi(\tilde{Q})$, together with one sign that the frame introduces: the images of $e_0$ and of the three $e_k$ commute with $\gamma^0$, while the images of the three $ie_k$ anticommute with it. Hence

$$ \gamma^\mu\gamma^\nu = \epsilon_\mu\,\Phi(e^\mu e^\nu), \qquad \epsilon_0 = +1, \quad \epsilon_k = -1 , $$

a relation verified for all sixteen index pairs in the Dirac representation. The mass term of the equation, the antisymmetrized pair, comes out of the same rule as $\gamma^{0k} = \Phi(ie_k)$ for the timelike pairs and $\gamma^{jk} = \Phi(\pm e_l)$, $l$ the remaining spatial index, for the spacelike pairs; the ambiguous sign is the one fixed by the orientation $e_1e_2 = e_3$. This is the dictionary's statement that the timelike bivectors are the images of the informational sector and the spacelike ones the images of the material sector. The trivectors, which carry the kinetic structure, are $\gamma^{\mu\nu\rho} = \gamma^0\,\Phi(w)$ with $w \in \mathbb{M}_-$, the four of them being the images of $e_3$, $-e_2$, $e_1$ and $-i$, again as the dictionary records. Chirality is the one place where the correspondence needs a word. The dictionary assigns the biquaternion imaginary the Clifford image of minus the pseudoscalar, $\Phi(i)=-\omega$, so the chirality operator is the product of the two commuting square roots of $-1$,

$$ \gamma_5 = i_{\mathrm{Cl}}\,\omega = -\Phi_{\mathbb{C}}\!\left(i_{\mathrm{Cl}}\,i\right), $$

the external complex unit $i_{\mathrm{Cl}}$ of the complexified Clifford algebra and the biquaternion imaginary $i$. It is not an element of $\mathbb{B}$: it needs both units. It anticommutes with every generator and commutes with every even element, hence with the whole biquaternion image, and its eigenspaces are the chiral halves of the spinor space. The transcriptions are the same ones used for the Dirac equation and are the content of the companion dictionary between the Dirac algebra and the biquaternions.

The equivalent constraint form

For the free field the antisymmetrized equation is equivalent to the pair of equations

$$ (i\not\partial-m)\psi_\mu = 0 \quad\text{for each }\mu, \qquad \gamma^\mu\psi_\mu = 0 , $$

in the following sense. If each component satisfies the Dirac equation, if the trace constraint holds and if the divergence constraint $\partial^\mu\psi_\mu=0$ holds, then the antisymmetrized equation follows; conversely, the antisymmetrized equation implies the component Dirac equations and the trace, and in the massive case the divergence as well. The reduction is standard for the free theory, and it is the form adapted to the algebra, because the Dirac equation on a biquaternion component is exactly the biquaternionic Dirac equation of the companion article, and the trace constraint is an algebraic condition on the index contraction. In the rest frame the two descriptions have the same four-dimensional solution space, and the explicit computation below verifies this.

The biquaternion trace constraint is the vanishing of the contraction of the vector index with the basis units,

$$ \tilde{\Psi}^{\text{tr}} := \sum_\mu e_\mu\,\tilde{\Psi}_\mu = 0, $$

up to the index-position conventions that translate $\gamma^\mu$ into the basis, and its independence is what makes the constraint algebraic. In the explicit rest-frame computation below the contraction is carried out with the Dirac matrices directly, so that no convention is left implicit.

The Rest-Frame Reduction

Setting up

Take the field at rest and of positive frequency, so that each component is proportional to $e^{-imt/\hbar}$ in natural units, and the derivative is $\partial_0\to -im$. The Dirac operator becomes

$$ i\not\partial - m \;\longrightarrow\; i\gamma^0(-im) - m = m\left(\gamma^0 - I_4\right), $$

acting on each component. Since $\gamma^0 = \mathrm{diag}(I_2,-I_2)$, the operator $m(\gamma^0-I_4)$ annihilates exactly the upper components and annihilates nothing else; the Dirac equation on the spatial components therefore forces them to their upper two-component parts.

The trace constraint at rest

Write $\psi_0$ and $\psi_k$ for the components, and use the trace constraint $\gamma^\mu\psi_\mu = \gamma^0\psi_0 + \gamma^k\psi_k = 0$. The Dirac equation forces $\psi_0$ upper and $\psi_k$ upper. Writing $\psi_k = (\chi_k, 0)$ with $\chi_k$ a two-component spinor and using $\gamma^k = \big(\begin{smallmatrix}0&\sigma^k\\-\sigma^k&0\end{smallmatrix}\big)$,

$$ \gamma^k\psi_k = \begin{pmatrix} \sigma^k\chi_k \\ -\sigma^k\chi_k \end{pmatrix}, $$

so the trace constraint reads, in its two chiral halves,

$$ \psi_0 + \boldsymbol{\sigma}\cdot\boldsymbol{\chi} = 0 \quad\text{(upper)}, \qquad -\boldsymbol{\sigma}\cdot\boldsymbol{\chi} = 0 \quad\text{(lower)} . $$

The lower equation gives $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=0$; the upper then gives $\psi_0 = -\boldsymbol{\sigma}\cdot\boldsymbol{\chi} = 0$. The trace constraint therefore forces the time component to vanish and imposes the single vector condition $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=0$ on the three spatial amplitudes. The combined system of the Dirac equation on all four components and the trace constraint has nullity four, the constraint $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=0$ has nullity four, and the two solution spaces coincide.

The spin-one-half content and its removal

The three spatial amplitudes $\chi_k$ carry $3\times2 = 6$ complex components before the condition. Under the rotation group the spatial index is a spin-one object and the two-component spinor a spin-$\tfrac12$ object, so

$$ 1\otimes\tfrac12 = \tfrac32 \oplus \tfrac12, \qquad 3\times2 = 4 + 2 , $$

and the condition $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=0$ removes exactly the spin-$\tfrac12$ piece of dimension two. What remains is a four-dimensional space, the spin-$\tfrac32$ quartet, with states of spin projection $\pm\tfrac32$ and $\pm\tfrac12$ along any axis. The representation-theory count of the field had two spin-$\tfrac32$ multiplets; the one found here is the positive-frequency or particle multiplet, and the negative-frequency or antiparticle multiplet is obtained by the conjugate choice of frequency and has the same four states. The four spin-$\tfrac12$ multiplets of the field are removed, two by the trace constraint and two by the Dirac equation, exactly as the general consistency analysis requires.

Verification

The reduction was checked by explicit complex-matrix computation in the Dirac representation, with $\gamma^0=\mathrm{diag}(I_2,-I_2)$, $\gamma^k=\big(\begin{smallmatrix}0&\sigma^k\\-\sigma^k&0\end{smallmatrix}\big)$, and $m$ a fixed nonzero parameter. At rest the Dirac operator on the sixteen amplitudes is $m(\gamma^0-I_4)$ on each component block; the trace constraint was assembled as the $4\times16$ matrix of $\gamma^\mu$ on the four component blocks; and the null space of the combined $20\times16$ system was computed by direct elimination. Its dimension is four. The null space of the single condition $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=0$ on the six spatial upper components is also four, and the two spaces agree element by element. The antisymmetrized form of the equation gives the same four-dimensional null space, confirming the equivalence of the trace form and the antisymmetrized form at the level of the free solutions. The check was performed on a basis of the solution space, not on a single wavefunction, and the representation is the one declared above.

Three further checks were run at a generic on-shell momentum, where $\mathbf{k}$ has three distinct nonzero components, so that no rest-frame degeneracy can hide a mismatch. The duality identity was verified as a matrix identity in the free vector index at four momenta, generic and lightlike, with residual zero, and the two writings of the massive equation were verified to have the same null space. The antisymmetrized equation has null space of dimension $4$, and it coincides with the null space of the component-Dirac-plus-trace system and with that of the same system augmented by the contracted divergence condition $\partial^\mu\psi_\mu=0$, so that in the massive case the trace and the divergence are consequences of the equation rather than additional conditions; a single-component amplitude was checked to fail the equation by an amount of order $m$, so that the solution space is genuinely a constrained subspace and not the whole space. For the massless equation at a generic lightlike momentum the same computation gives a null space of dimension $6$, of which the four-dimensional space of pure-gauge modes $\psi_\mu\propto\partial_\mu\epsilon$ is a subspace; the trace condition is not implied, the divergence condition is, and the two-dimensional quotient is the two helicities. The trace values of the massless solutions span exactly the two-dimensional space of trace shifts that the gauge transformation can produce, which is the precise sense in which the trace condition is a gauge choice and the counting is unchanged.

The Spin-Three-Halves Projector

The constraint $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=0$ defines a subspace of the six-dimensional space of spatial amplitudes — three vector directions, two spinor components each — and the projector onto it is elementary. Define, on the vector index,

$$ P_{ij} = \delta_{ij} - \tfrac13\,\sigma_i\sigma_j . $$

Then $P$ is idempotent,

$$ P_{ik}P_{kj} = \delta_{ij} - \tfrac23\sigma_i\sigma_j + \tfrac19\sigma_i\sigma_k\sigma_k\sigma_j = \delta_{ij} - \tfrac23\sigma_i\sigma_j + \tfrac13\sigma_i\sigma_j = \delta_{ij} - \tfrac13\sigma_i\sigma_j = P_{ij}, $$

using $\sigma_k\sigma_k = 3I_2$, and it annihilates exactly the longitudinal direction: for $\chi_k$ proportional to $\sigma_k$ acting on any spinor, $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=\sigma_k\sigma_k\xi = 3\xi \neq 0$, so the longitudinal combination is removed, while for $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=0$ one has $P_{ij}\chi_j = \chi_i$. Its trace in the vector index is $3-1=2$, so it projects the three spatial directions onto two, and combined with the two spinor components it leaves the four states of spin $\tfrac32$.

The projector is the rest-frame form of the covariant spin-$\tfrac32$ projector of the standard theory, which is obtained by replacing the rest-frame $\gamma^0$ by $\not p/m$ and the spatial delta by the transverse projector with respect to $p^\mu$. In the biquaternion framework the boost is effected by rotor conjugation, so the covariant projector is the rotor conjugate of the rest-frame combination above; the algebra supplies the boost without a separate construction.

The Massless Case and the Gauge Freedom

The gauge transformation in the algebra

In the massless case the antisymmetrized equation and the trace constraint leave a gauge freedom, and in the biquaternion form the transformation is

$$ \tilde{\Psi}_\mu \;\longrightarrow\; \tilde{\Psi}_\mu + \tilde{\nabla}_\mu\,\tilde{\epsilon}, $$

where $\tilde{\epsilon}$ is a biquaternion playing the role of the Dirac-spinor parameter and $\tilde{\nabla}_\mu$ is the $\mu$-th component of the gradient. The transformation shifts the field by a gradient, and the antisymmetrized kinetic structure of the equation is insensitive to it, just as the electromagnetic field strength is insensitive to a gradient of the potential. The number of gauge parameters is the number of components of $\tilde{\epsilon}$, namely four, which is one Dirac spinor; and the gauge freedom removes two of the four spin-$\tfrac12$ multiplets of the representation content. The remaining two are removed by the trace condition and the equations of motion — with the difference, recorded above, that in the massless case the trace is a gauge choice and the divergence condition is the consequence of the equation, so that the two mechanisms do not simply add but interlock. What is left is a single spin-$\tfrac32$ multiplet per frequency sign, that is, the two helicities $\pm\tfrac32$.

Helicities and the two physical states

The physical content of the massless field is therefore two states per wavevector, the helicities $\pm\tfrac32$. There is no longitudinal or scalar helicity, and the count is the spin-$\tfrac32$ analogue of the two transverse helicities of the photon: in both cases the gauge freedom and the constraint together halve-and-halve the naive index content. The parallel is exact at the level of the representation theory. The photon's two-form is the antisymmetric square of the four-vector, $(1,0)\oplus(0,1)$; the massless Rarita–Schwinger field's two physical states are the two helicities of the $(1,\tfrac12)$ and $(\tfrac12,1)$ summands after the spin-$\tfrac12$ pieces are gauged away. In both cases the framework's self-dual split is the algebraic trace of the helicity structure, and in both cases the massless limit is the case in which the third and lower states are removed by gauge rather than retained by a mass.

The Massive Case

For $m\neq0$ the mass term

$$ - m\,\gamma^{\mu\nu}\tilde{\Psi}_\nu $$

is the antisymmetrized pair of Clifford generators acting on the field, and its biquaternion image is the same antisymmetrized pair of representatives, contracted with the corresponding component of the field; in the trace form of the equation it is simply the Dirac mass on each component. The mass term plays the same structural role as in the Proca and Dirac cases: it is linear in the field, it is what remains when the gauge freedom is absent, and it fixes the dispersion relation on the physical shell through the Dirac operator. In the rest frame the massive reduction above is the whole story: time component zero, spatial components upper, $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=0$, four states.

The massive case is also where the known difficulty lives. The constraint analysis that produces $\gamma^\mu\psi_\mu=0$ and $\partial^\mu\psi_\mu=0$ is rigid when $m\neq0$, and its consequences for the interacting theory are the Velo–Zwanziger obstruction: minimal coupling to a background electromagnetic field makes some modes propagate faster than light, and the theory is not causal as a fixed-background field theory. The biquaternion formulation reproduces the free constraint structure exactly and says nothing different about the coupling, because the coupling is an additional structure that no rewriting of the free equation can supply. This is the analogue, for spin $\tfrac32$, of the statement recorded for the Proca field: the algebra represents the free field faithfully, and what the free algebra cannot decide is the interaction.

The dispersion relation

Because each component satisfies the Dirac equation, each component satisfies the Klein–Gordon equation on squaring, exactly as in the spin-$\tfrac12$ case. For a biquaternion plane wave of the component,

$$ \tilde{\Psi}_\mu = \tilde{\Psi}_{\mu,0}\,e^{\,i\,\mathrm{Sc}(\tilde{K}\tilde{Q}^{\natural})}, $$

the Dirac equation on the component requires $\tilde{K}\tilde{K}^{\natural} = -m^2c^2/\hbar^2$ in the normalization of the companion Dirac article, that is,

$$ \frac{\omega^2}{c^2} - \mathbf{k}^2 = \frac{m^2c^2}{\hbar^2}, \qquad\text{i.e.}\qquad \omega^2 = c^2\mathbf{k}^2 + \frac{m^2c^4}{\hbar^2}, $$

the same physical massive shell as the Proca field. The trace constraint and the subsidiary conditions then reduce the sixteen components to the four states of the quartet, with the three-momentum on the shell. The dispersion relation is therefore common to the massive spin-one and spin-$\tfrac32$ fields, as it must be: the mass shell is a property of the representation's Casimir and not of the index structure.

Comparison with Spin One and Spin One-Half

The three field equations of this subcategory and its lower-spin neighbours form a sequence, and the sequence is the clearest statement of what the biquaternion framework does and does not carry natively.

Field	Standard carrier	Components	Constraints	Physical components
Dirac, spin $\tfrac12$	Dirac spinor	$4$	the equation alone, removing $2$	$2$
Proca, spin $1$	four-vector	$4$	Lorenz, removing $1$	$3$
Rarita–Schwinger, spin $\tfrac32$	vector-spinor	$16$	trace and subsidiary, removing $12$	$4$

In each row the physical component count is $2s+1$ for the particle, and the equation together with its constraints removes the remainder. The role that the field's index structure plays is to supply a carrier large enough to contain the spin: the Dirac spinor is the defining module of the algebra and needs no index; the Proca field adds a four-vector index, which the algebra carries by conjugation on $\mathbb{M}_-$; the Rarita–Schwinger field adds a vector index to the spinor, which the algebra carries only as a tensor product of the two. Each step multiplies the carrier by the four-vector and each step therefore adds lower-spin admixture that a constraint has to remove.

The biquaternion description follows the same sequence. The Dirac field is an element of the spinor module or of the algebra's minimal left ideal; the Proca field is an element of the material sector $\mathbb{M}_-$ with the field strength in the vector part; the Rarita–Schwinger field is a family of elements of $\mathbb{B}$ indexed by the spacetime direction. The pattern is that the algebra carries the four-vector index by conjugation and the spinor index by left multiplication, and that no carrier beyond these two is available. This is the structural observation that the final article of the subcategory turns into a statement about the algebra's module category: the spin content accessible to the algebra's own modules is bounded, and it is $\{0,\tfrac12\}$.

What Is Standard and What the Algebra Adds

The standard part. The vector-spinor field, the antisymmetrized equation, the trace constraint, the gauge freedom, the reduction to the Dirac equation, the spin-$\tfrac32$ content, and the Velo–Zwanziger difficulty are all standard and are transcribed here. The representation-theory count $2\times\tfrac32\oplus4\times\tfrac12$ and the reduction $1\otimes\tfrac12=\tfrac32\oplus\tfrac12$ are standard.

The algebra's part. Three statements are the framework's. First, the natural carrier of the field is the family of biquaternion-valued components $\tilde{\Psi}_\mu$, and the two structural conditions of the free theory — the Dirac equation on each component and the trace constraint $\sum_\mu e_\mu\tilde{\Psi}_\mu=0$ — are each single algebraic statements rather than lists of component equations. Second, chirality is carried by the two commuting square roots of $-1$ that the Dirac dictionary identifies, the biquaternion central imaginary and the external complex unit of the complexified Clifford algebra, $\gamma_5=-\Phi_{\mathbb{C}}(i_{\mathrm{Cl}}i)$; it commutes with the whole biquaternion image, so the vector-spinor's chirality structure needs no independent operator beyond the two units the dictionary already supplies. Third, the rest-frame projector $P_{ij}=\delta_{ij}-\tfrac13\sigma_i\sigma_j$ is obtained from the constraint by elementary algebra, and its boosted form is the rotor conjugate of the rest-frame object, so the covariant spin-$\tfrac32$ projector is available without a separate construction.

What is imported. The identification of the field with a spin-$\tfrac32$ particle, the value of the mass, the existence of a consistent interaction, and the resolution of the Velo–Zwanziger problem are inputs beyond the free algebra. The canonical quantisation of the field and its propagator belong to the field-theoretic companion subcategory.

The structural point that this article makes for the subcategory is that the field is a tensor product of the algebra's vector structure with its spinor structure, so its spin-$\tfrac32$ content is not native to the algebra but arises from the product. The final article of the subcategory states this in general and proves that the algebra's own modules carry only spin zero and spin one-half.

The conventions are those of three companion articles:

Companion article The Dirac Equation in Biquaternionic Form, for the biquaternionic Dirac operator and the component Dirac equation.
Companion article The Dirac Algebra and Biquaternions — A Dictionary, for the translation of the Clifford generators and their antisymmetrized products into the algebra's units.
Companion article The Self-Dual and Anti-Self-Dual Split: Spin 1 from the Biquaternion Material Sector, for the vector part of the algebra and the four-vector representation from which the vector-spinor is built.

Summary

The Rarita–Schwinger field is a vector-spinor of sixteen complex components, with the massless equation $\epsilon^{\mu\nu\rho\sigma}\gamma_5\gamma_\nu\partial_\rho\psi_\sigma=0$ and the massive equation the same with the mass term $+m\gamma^{\mu\nu}\psi_\nu$, equivalently, in the trivector form, $(i\gamma^{\mu\nu\rho}\partial_\rho-m\gamma^{\mu\nu})\psi_\nu=0$. The two writings are related by the duality identity $\epsilon^{\mu\nu\rho\sigma}\gamma_5\gamma_\nu\partial_\rho\psi_\sigma=-i\gamma^{\mu\sigma\rho}\partial_\rho\psi_\sigma$, which fixes the relative sign of their mass terms; the overall sign is the free convention, because $\psi_\mu\to\gamma_5\psi_\mu$ exchanges the two possibilities. In the massive case the equation implies the algebraic trace $\gamma^\mu\psi_\mu=0$ and the divergence condition, and it is equivalent to the component Dirac equation together with the trace; the massless theory instead has the gauge freedom $\psi_\mu\to\psi_\mu+\partial_\mu\epsilon$ with a single Dirac-spinor parameter, under which the trace is a complete gauge choice while the divergence remains a consequence of the equation.

In the biquaternion framework the carrier is a family of biquaternion-valued components $\tilde{\Psi}_\mu\in\mathbb{B}$, one for each spacetime index, and the free system is the pair consisting of the biquaternionic Dirac equation on each component and the algebraic trace constraint $\sum_\mu e_\mu\tilde{\Psi}_\mu=0$. Chirality is the product of the two square roots of $-1$ of the dictionary, $\gamma_5=-\Phi_{\mathbb{C}}(i_{\mathrm{Cl}}i)$, and the Clifford generators are the fixed frame times the images of the basis units of the informational sector, $\gamma^\mu=\gamma^0\Phi(e^\mu)$ with $e^0=e_0$ and $e^k=ie_k$.

In the rest frame the system forces the time component to vanish and the spatial amplitudes to satisfy $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=0$. The surviving space is four-dimensional and carries the spin-$\tfrac32$ quartet; the projector onto it is $P_{ij}=\delta_{ij}-\tfrac13\sigma_i\sigma_j$, idempotent and of vector-index rank two. The representation content $2\times\tfrac32\oplus4\times\tfrac12$ loses its four spin-$\tfrac12$ multiplets to the constraint and the gauge freedom, and the massless field retains the two helicities $\pm\tfrac32$.

The field is a tensor product of the algebra's vector structure and its spinor structure, so spin $\tfrac32$ is not native to the algebra. The free equation is represented faithfully; the known difficulty of the interacting theory, the Velo–Zwanziger acausal propagation under minimal coupling, is a property of the coupling and is not removed by the rewriting.

Summary of Notation

Symbol	Meaning
$\mathbb{B}=\mathbb{C}\otimes_\mathbb{R}\mathbb{H}$	Biquaternion algebra
$e_0=1,e_1,e_2,e_3$	Quaternion basis, $e_k^2=-e_0$
$i$	Central scalar imaginary of $\mathbb{B}$, $i^2=-1$
$\mathbb{M}_-$, $\mathbb{M}_+$	Material (anti-Hermitian) and informational (Hermitian) sectors
$\Phi$, $\gamma^\mu=\gamma^0\Phi(e^\mu)$	Dictionary map onto the even Clifford algebra; $e^0=e_0$, $e^k=ie_k$
$i_{\mathrm{Cl}}$, $\gamma_5=-\Phi_{\mathbb{C}}(i_{\mathrm{Cl}}i)$	External imaginary of the complexified Clifford algebra and the chirality operator
$\tilde{\Psi}_\mu\in\mathbb{B}$	Biquaternion-valued components of the vector-spinor
$\gamma^\mu$, $\{\gamma^\mu,\gamma^\nu\}=2g^{\mu\nu}I_4$	Clifford generators, $g=\mathrm{diag}(+1,-1,-1,-1)$
$\gamma^{\mu\nu}=\tfrac12[\gamma^\mu,\gamma^\nu]$	Antisymmetrized pair, the mass-term structure
$\gamma^{\mu\nu\rho}$	Antisymmetrized triple, kinetic structure of the equation
$\epsilon^{\mu\nu\rho\sigma}\gamma_5\gamma_\nu\partial_\rho\psi_\sigma=-i\,\gamma^{\mu\sigma\rho}\partial_\rho\psi_\sigma$	Duality identity relating the two writings of the equation
$\tilde{\Psi}^{\text{tr}}=\sum_\mu e_\mu\tilde{\Psi}_\mu$	Trace constraint, image of $\gamma^\mu\psi_\mu=0$
$\psi_\mu=(\chi_\mu,\eta_\mu)$	Upper and lower two-component parts in the Dirac representation
$\boldsymbol{\chi}$	Spatial amplitudes; constraint $\boldsymbol{\sigma}\cdot\boldsymbol{\chi}=0$
$P_{ij}=\delta_{ij}-\tfrac13\sigma_i\sigma_j$	Spin-$\tfrac32$ projector on the vector index
$\tilde{\nabla}_\mu$	$\mu$-th component of the biquaternionic gradient
$2\times\tfrac32\oplus4\times\tfrac12$	$SU(2)$ content of the vector-spinor