Canonical Quantization of the Biquaternion Dirac Field

Introduction

The companion article Quantum Mechanics in Biquaternionic Form raises, as the second of its list of open questions, a precise one:

"The formalism is first-quantized. To describe creation and annihilation of particles, one needs a Fock space and operator-valued fields. How does the biquaternion framework extend to this setting?"

The companion article The Dirac Equation in Biquaternionic Form asks the same question from the field side, and the article The Dirac Equation in Biquaternionic Form — Solutions and the Non-Relativistic Limit supplies the classical material a mode expansion needs: the positive- and negative-frequency plane-wave solutions, their covariant normalizations, and the spin sums. The article The KMS Condition and the Biquaternion Framework supplies the thermal motivation and identifies exactly what a fermionic extension must add: a $\mathbb{Z}/2$ grading under which fermionic operators anticommute, which the biquaternion algebra does not by itself provide.

This article is a first attempt at that extension. Its subject is the canonical quantization of the biquaternion Dirac field: the promotion of the classical field to an operator-valued field, the equal-time anticommutation relations, the expansion in the plane-wave solutions established by the parent article, the algebra of the mode operators, the Hamiltonian and its normal-ordering constant, the reason the modes anticommute rather than commute, and the relation of the resulting operator algebra to the thermal-state framework of the KMS article.

The article has to be read with a clear division in mind, and we state it at the outset because it governs everything that follows.

The standard part. The spinor-module form of the biquaternion Dirac equation is the ordinary Dirac equation, $(i\gamma^\mu\partial_\mu-m)\psi=0$ in the notation of the companion articles. Its canonical quantization — the Lagrangian, the conjugate momentum, the equal-time anticommutators, the mode expansion, the mode algebra, the Hamiltonian, and the normal-ordering constant — is textbook quantum field theory, transcribed into the notation of this series. None of it is new, and none of it depends on the biquaternion structure beyond the kinematical conventions already fixed by the companion articles.
The open part. What is not standard, and what this article cannot settle, is the intrinsic formulation of the quantization in the algebra $\mathbb{B}$ itself: whether the operator-valued field should be taken to be $\mathbb{B}$-valued or valued in the spinor module on which $\mathbb{B}\cong M_2(\mathbb{C})$ acts; what the correct $\mathbb{B}$-intrinsic Lagrangian and conjugate momentum are; what the biquaternion Fock space is; and how the $\mathbb{Z}/2$ grading required by fermionic statistics is to be embedded in $\mathbb{B}$.

We therefore make no claim that the extension is canonical in the biquaternion framework. The claim is narrower: the extension can be carried out on the spinor module with the parent article's solutions, and it supplies exactly the operator structure — anticommuting modes and a fermion-parity grading — that the KMS article identifies as missing. Whether that structure lifts to a genuinely biquaternionic one is the open question this article is meant to expose.

The article is organized as follows. The next section fixes the classical field and its conjugate momentum, and identifies the constraint that forces the canonical structure to be imposed rather than assumed. The following section states the equal-time anticommutation relations and the reasons they, rather than commutators, are the consistent choice. The next section carries out the expansion in the parent article's plane-wave solutions. The section after that derives the anticommutators of the mode operators and builds the Fock space. The next section gives the Hamiltonian, the normal-ordering constant, and the four-momentum and charge. A section on spin and statistics records why anticommutation is forced here, without developing the full spin–statistics theorem, which is the subject of a planned companion article. A section relates the result to the KMS thermal-state framework. The article closes with an explicit accounting of what is standard and what is open.

Conventions. We use those of the companion articles throughout. The biquaternion algebra is $\mathbb{B}=\mathbb{C}\otimes_\mathbb{R}\mathbb{H}$, the quaternion basis is $e_0=1,e_1,e_2,e_3$ with $e_k^2=-e_0$, and $i$ is the scalar imaginary with $i^2=-1$. The biquaternionic gradient is $\tilde{\nabla}=e_0\partial_{ict}+e_1\partial_x+e_2\partial_y+e_3\partial_z$, with $\Box=\tilde{\nabla}\tilde{\nabla}^{\natural}=\tilde{\nabla}^{\natural}\tilde{\nabla}$. The Clifford metric of the gamma matrices is $g=\mathrm{diag}(+1,-1,-1,-1)$, through $\{\gamma^\mu,\gamma^\nu\}=2g^{\mu\nu}I_4$, and the spacetime metric of the $ict$ gradient is $\eta=\mathrm{diag}(-1,+1,+1,+1)=-g$; the relative sign is the convention the companion article flags. The Dirac adjoint is $\bar\psi=\psi^\dagger\gamma^0$. In the solution and quantization sections we work in natural units $\hbar=c=1$, as the parent solution article does, restoring $\hbar$ and $c$ only where they carry meaning, namely in the mass-shell relation and in the thermal frequency.

The Classical Field and Its Conjugate Momentum

The biquaternion Dirac field of the companion article is the map $\tilde{\Psi}:\mathbb{R}^{1,3}\to\mathbb{B}$ satisfying the linear, chirality-off-diagonal pair

$$ \tilde{\nabla}\tilde{\Psi}_R=m\tilde{\Psi}_L,\qquad \tilde{\nabla}^{\natural}\tilde{\Psi}_L=m\tilde{\Psi}_R, $$

with massless limit $\tilde{\nabla}\tilde{\Psi}=0$; here $\tilde{\Psi}=\tilde{\Psi}_L+\tilde{\Psi}_R$ carries one component per chirality.

The companion solutions article works with its spinor-module representative, the four-component complex Dirac spinor $\psi$, on whose underlying two-component spinor module the algebra $\mathbb{B}\cong M_2(\mathbb{C})$ acts. There the equation is the standard Dirac equation

$$ (i\gamma^\mu\partial_\mu-m)\psi=0, $$

whose solutions, normalizations, and spin sums are established in the parent article and are used below without rederivation.

The Lagrangian

The classical field theory is defined by the Lagrangian density

$$ \mathcal{L}=\bar{\psi}\,(i\gamma^\mu\partial_\mu-m)\,\psi, $$

a real Lorentz scalar. Variation with respect to $\bar{\psi}$ returns the Dirac equation $(i\gamma^\mu\partial_\mu-m)\psi=0$, and variation with respect to $\psi$ returns its adjoint; the two are equivalent under $\bar\psi=\psi^\dagger\gamma^0$. This is the standard Lagrangian of a free Dirac field, written in the gamma-matrix convention of the companion articles.

The Conjugate Momentum and Its Constraint

The canonical momentum conjugate to $\psi$ is

$$ \pi=\frac{\partial\mathcal{L}}{\partial(\partial_0\psi)}=i\,\psi^\dagger . $$

The step is worth displaying, because it is the origin of the entire fermionic structure. The only term in $\mathcal{L}$ containing the time derivative is $\bar{\psi}\,i\gamma^0\partial_0\psi$. Since $\bar{\psi}\gamma^0=\psi^\dagger$ and $(\gamma^0)^2=I_4$, that term is $i\,\psi^\dagger\partial_0\psi$. Differentiating with respect to $\partial_0\psi_a$ therefore returns $i\psi^\dagger_a$, i.e. $\pi=i\psi^\dagger$.

The result is not $\dot{\psi}$, and this is the essential difference from a scalar field. The Dirac equation is first order in time, so the momentum is algebraically tied to the field by

$$ \chi=\pi-i\psi^\dagger\approx 0 . $$

Equation $\chi\approx0$ is a second-class constraint: the Legendre transform is singular, the field and its momentum are not independent canonical variables, and the naive Poisson bracket must be replaced by the Dirac bracket. The standard resolution is not to quantize the brackets of the constrained system by hand but to impose the canonical (anti)commutation relations directly and verify that they are consistent with the field equation. The constraint is the technical reason a Dirac field is quantized by imposing relations on $\psi$ and $\pi$ rather than by the ordinary rule $[\phi,\pi]=i\delta$; it is also the reason the choice of bracket — commutator versus anticommutator — is not a free convention, as the next section explains. We do not reproduce the Dirac-bracket construction; it is textbook, and the companion articles do not need it.

The Biquaternion Reading and the Open Point

Written on the spinor module, the Lagrangian and the canonical momentum are the standard ones. Written in the algebra, they are not yet fixed. The parent article's mass term is the linear, chirality-off-diagonal coupling, whose spinor-module representative is the ordinary mass term $m\psi$; the anti-Hermitian conjugation $\flat=-{}^{*}$ is the algebra's real structure, not the mass. The $\mathbb{B}$-intrinsic Lagrangian — the real scalar built from $\tilde{\nabla}\tilde{\Psi}$, $\tilde{\Psi}$, their conjugates, and the trace or norm pairing of the algebra — is not determined by the algebra alone, and neither is the intrinsic conjugate momentum. Two inequivalent transcriptions are available: one may treat the operator field as $\mathbb{B}$-valued and pair it with the real part of the trace, or one may treat it as valued in the spinor module and quantize there. We adopt the second, which is the parent article's working setting, and we return to the first in the closing section on what is open. The reader should keep in mind that the mode structure below is the standard one; what is being tested here is whether the biquaternion framework naturally reproduces it.

Equal-Time Anticommutation Relations

Quantization promotes the classical field and its momentum to operators $\hat{\psi}$, $\hat{\pi}$ and imposes equal-time brackets. For the Dirac field the brackets are anticommutators:

$$ \{\hat{\psi}_a(\mathbf{x},t),\hat{\psi}_b^\dagger(\mathbf{y},t)\}=\delta_{ab}\,\delta^{(3)}(\mathbf{x}-\mathbf{y}), $$

$$ \{\hat{\psi}_a(\mathbf{x},t),\hat{\psi}_b(\mathbf{y},t)\}=0,\qquad \{\hat{\psi}_a^\dagger(\mathbf{x},t),\hat{\psi}_b^\dagger(\mathbf{y},t)\}=0, $$

where $a,b$ are spinor indices. In terms of the conjugate momentum $\hat{\pi}=i\hat{\psi}^\dagger$, the first relation reads

$$ \{\hat{\psi}_a(\mathbf{x},t),\hat{\pi}_b(\mathbf{y},t)\}=i\,\delta_{ab}\,\delta^{(3)}(\mathbf{x}-\mathbf{y}), $$

which is the fermionic analogue of the canonical relation $[\hat{\phi},\hat{\pi}]=i\delta^{(3)}$ of a scalar field. These are the standard canonical relations, and they are imposed rather than derived.

The choice of anticommutator over commutator is forced, and three reasons can be given at the level of this article. They are stated here because they are operative in the construction; the full spin–statistics theorem is the subject of a planned companion article and is not developed here.

1. Boundedness of the energy. As shown in the section on the Hamiltonian, the mode expansion of the energy contains a negative-frequency contribution proportional to $\hat{b}\,\hat{b}^\dagger$. If the modes commuted, $\hat{b}\hat{b}^\dagger=\hat{b}^\dagger\hat{b}+e_0$, and the energy would be $H=\sum_{\mathbf p,r}E_{\mathbf p}(\hat{a}^\dagger\hat{a}-\hat{b}^\dagger\hat{b})+\text{const}$, unbounded below: adding antiparticles would lower the energy without limit, and no stable vacuum would exist. If the modes anticommute, $\hat{b}\hat{b}^\dagger=e_0-\hat{b}^\dagger\hat{b}$, and the same expression becomes a positive sum plus a c-number. The anticommutator is what makes the negative-frequency modes into positive-energy antiparticles. This is the operator content of the Dirac sea.

2. Microcausality. For a field of spin $\tfrac12$, the commutator $[\hat{\psi}(x),\bar{\hat{\psi}}(y)]$ does not vanish for spacelike separated $x,y$, so quantizing with commutators would correlate measurements in causally disconnected regions. The anticommutator $\{\hat{\psi}(x),\bar{\hat{\psi}}(y)\}$ does vanish there, and it is the quantity that must be local for a fermion field. This is the content of Pauli's 1940 argument: the requirement that the fields commute or anticommute outside the light cone, together with the transformation properties of the spinor representation, selects anticommutation for half-integer spin.

3. Positivity of the inner product. Quantizing a half-integer-spin field with commutators produces states of negative norm; the anticommutator is what makes the Fock inner product positive-definite. This is the Hilbert-space form of the same obstruction.

The three reasons are aspects of one theorem. We note only that they all point the same way for a field whose one-particle states transform in the spinor representation — which is the case here, since the biquaternion algebra $\mathbb{B}$ acts on a two-dimensional complex spinor module.

Field Expansion in the Plane-Wave Solutions

The classical equation is linear, so its general solution is a superposition of the plane waves constructed in the parent article. Those solutions come in two branches: the positive-frequency branch spanned by the spinors $u^{(r)}(\mathbf{p})e^{-ip\cdot x}$ and the negative-frequency branch spanned by $v^{(r)}(\mathbf{p})e^{+ip\cdot x}$, with $r=1,2$,

$$ E_{\mathbf p}=+\sqrt{\mathbf p^2+m^2},\qquad p\cdot x=E_{\mathbf p}t-\mathbf p\cdot\mathbf x, $$

and with the covariant normalizations and spin sums established there:

$$ \bar u^{(r)}u^{(s)}=2m\,\delta^{rs},\quad \bar v^{(r)}v^{(s)}=-2m\,\delta^{rs},\quad u^{(r){}^{*}}u^{(s)}=v^{(r){}^{*}}v^{(s)}=2E_{\mathbf p}\,\delta^{rs}, $$

$$ \sum_{r=1}^{2}u^{(r)}(\mathbf p)\bar u^{(r)}(\mathbf p)=\not p+m,\qquad \sum_{r=1}^{2}v^{(r)}(\mathbf p)\bar v^{(r)}(\mathbf p)=\not p-m . $$

Quantization consists in replacing the constant coefficients of the classical expansion by operators:

$$ \hat{\psi}(x)=\int\!\frac{d^3p}{(2\pi)^3}\,\frac{1}{\sqrt{2E_{\mathbf p}}} \sum_{r=1}^{2}\Big[\hat a_r(\mathbf p)\,u^{(r)}(\mathbf p)\,e^{-ip\cdot x} +\hat b_r^\dagger(\mathbf p)\,v^{(r)}(\mathbf p)\,e^{+ip\cdot x}\Big]. $$

The field and its adjoint,

$$ \hat{\psi}^\dagger(y)=\int\!\frac{d^3q}{(2\pi)^3}\,\frac{1}{\sqrt{2E_{\mathbf q}}} \sum_{s=1}^{2}\Big[\hat a_s^\dagger(\mathbf q)\,u^{(s){}^{*}}(\mathbf q)\,e^{+iq\cdot y} +\hat b_s(\mathbf q)\,v^{(s){}^{*}}(\mathbf q)\,e^{-iq\cdot y}\Big], $$

together define the operator-valued field. The annihilation operators $\hat a_r(\mathbf p)$ multiply the positive-frequency solutions, and the creation operators $\hat b_r^\dagger(\mathbf p)$ multiply the negative-frequency solutions. The attribution of a creation operator to the negative-frequency branch is the operator form of the reinterpretation of the Dirac sea: the modes of "negative frequency" are the antiparticles of positive energy.

Two structural remarks belong here. First, $\hat{\psi}$ is not Hermitian, and it should not be: a complex Dirac field carries a conserved charge, and particle and antiparticle are distinct. Second, the expansion is written on the spinor module, following the parent article. In the algebra, each mode function is the module representative of a biquaternion plane wave $\tilde{\Psi}_0\exp\!\left(i\,\mathrm{Sc}(\tilde{k}\tilde{Q}^{\natural})\right)$; the two branches are the two roots of the single biquaternion mass-shell condition

$$ \tilde k\tilde k^{\natural}=-\frac{m^2c^2}{\hbar^2}, $$

and the polarization equation is the algebraic form of the momentum-space equation. Promoting the coefficients to operators is thus a transcription of the standard mode expansion onto the parent article's biquaternion plane waves; it is not, by itself, a derivation of the expansion from the algebra.

Anticommutators of the Mode Operators

For the equal-time relations of Section 3 to hold, the mode operators must satisfy

$$ \{\hat a_r(\mathbf p),\hat a_s^\dagger(\mathbf q)\}=(2\pi)^3\,\delta_{rs}\,\delta^{(3)}(\mathbf p-\mathbf q), $$

$$ \{\hat b_r(\mathbf p),\hat b_s^\dagger(\mathbf q)\}=(2\pi)^3\,\delta_{rs}\,\delta^{(3)}(\mathbf p-\mathbf q), $$

with all other anticommutators vanishing:

$$ \{\hat a_r(\mathbf p),\hat a_s(\mathbf q)\}=0,\qquad \{\hat b_r(\mathbf p),\hat b_s(\mathbf q)\}=0,\qquad \{\hat a_r(\mathbf p),\hat b_s(\mathbf q)\}=0,\qquad \{\hat a_r(\mathbf p),\hat b_s^\dagger(\mathbf q)\}=0 . $$

These are the standard anticommutation relations of the free Dirac field, in the continuum normalization inherited from the parent article's covariant spinors.

Consistency with the Equal-Time Relations

It is worth verifying that these relations reproduce $\{\hat\psi_a(\mathbf x,t),\hat\psi_b^\dagger(\mathbf y,t)\}=\delta_{ab}\delta^{(3)}(\mathbf x-\mathbf y)$, because the verification exposes a point of care about the relative phases of the two branches. Inserting the expansions and using the mode algebra gives, at equal times,

$$ \{\hat{\psi}_a(\mathbf x,t),\hat{\psi}_b^\dagger(\mathbf y,t)\} =\int\!\frac{d^3p}{(2\pi)^3}\frac{1}{2E_{\mathbf p}} \Big[\sum_r u^{(r)}_a u^{(r){}^{*}}_b\,e^{+i\mathbf p\cdot\boldsymbol\Delta} +\sum_r v^{(r)}_a v^{(r){}^{*}}_b\,e^{-i\mathbf p\cdot\boldsymbol\Delta}\Big], $$

where $\boldsymbol\Delta=\mathbf x-\mathbf y$. The two phases are opposite: the positive-frequency term carries $e^{+i\mathbf p\cdot\boldsymbol\Delta}$ and the negative-frequency term $e^{-i\mathbf p\cdot\boldsymbol\Delta}$. From the spin sums one obtains the two useful combinations

$$ \sum_r u^{(r)}u^{(r){}^{*}}=E_{\mathbf p}I_4+m\gamma^0+\gamma^0\boldsymbol\gamma\cdot\mathbf p,\qquad \sum_r v^{(r)}v^{(r){}^{*}}=E_{\mathbf p}I_4-m\gamma^0+\gamma^0\boldsymbol\gamma\cdot\mathbf p, $$

so that, with $\theta=\mathbf p\cdot\boldsymbol\Delta$,

$$ \sum_r\Big[u^{(r)}u^{(r){}^{*}}e^{+i\theta}+v^{(r)}v^{(r){}^{*}}e^{-i\theta}\Big] =E_{\mathbf p}\cos\theta\,I_4\cdot 2 +2\cos\theta\,\gamma^0\boldsymbol\gamma\cdot\mathbf p +2im\sin\theta\,\gamma^0 . $$

The first term is even in $\mathbf p$ and gives $2E_{\mathbf p}\delta^{(3)}(\boldsymbol\Delta)$ after dividing by $2E_{\mathbf p}$; the second and third terms are odd in $\mathbf p$ and vanish under symmetric integration. Hence the equal-time relation holds. The identity was checked numerically for a representative momentum ($m=1.3$, $|\mathbf p|=0.9$, $\mathbf p$ along $e_3$): the spin sums above and their two combinations reproduce the right-hand sides to a maximum absolute error of $4.4\times10^{-16}$. The cancellation of the odd terms is exactly what the opposite phases of the two branches accomplish; had both phases been equal, the cross terms would have survived and the canonical relation would have failed.

The Fock Space and Pauli Exclusion

The mode algebra generates the Fock space. Define the vacuum $|0\rangle$ by

$$ \hat a_r(\mathbf p)|0\rangle=0=\hat b_r(\mathbf p)|0\rangle \quad\text{for all } r,\ \mathbf p , $$

and build states by acting with the creation operators. The one-particle states are $\hat a_r^\dagger(\mathbf p)|0\rangle$ and the one-antiparticle states $\hat b_r^\dagger(\mathbf p)|0\rangle$. At fixed $\mathbf p$, the positive-frequency solutions span a two-dimensional complex space, and the full one-particle space is the $L^2$ completion of that solution space; the many-particle space is its antisymmetric tensor algebra (exterior algebra), because the mode operators anticommute. In particular,

$$ \hat a_r^\dagger(\mathbf p)\hat a_s^\dagger(\mathbf q)|0\rangle =-\,\hat a_s^\dagger(\mathbf q)\hat a_r^\dagger(\mathbf p)|0\rangle, $$

so the two-particle state is antisymmetric, and

$$ \big(\hat a_r^\dagger(\mathbf p)\big)^2=0,\qquad \big(\hat b_r^\dagger(\mathbf p)\big)^2=0, $$

which is the Pauli exclusion principle as an algebraic identity: the square of a creation operator vanishes because $\{\hat a_r(\mathbf p),\hat a_s(\mathbf q)\}=0$. This is the sharpest sense in which the statistics is carried by the mode algebra. The construction of the biquaternion Fock space as an object in its own right — rather than as the exterior algebra of the spinor-module solution space — is the subject of a planned companion article; here we need only the mode algebra.

The Hamiltonian and the Normal-Ordering Constant

The conserved energy of the free Dirac field is

$$ \hat H=\int d^3x\;\hat{\psi}^\dagger\big(-i\gamma^0\boldsymbol\gamma\cdot\nabla+m\gamma^0\big)\hat{\psi} =\int d^3x\;\bar{\hat\psi}\big(-i\boldsymbol\gamma\cdot\nabla+m\big)\hat{\psi}, $$

the integrand being the $00$ component of the canonical energy–momentum tensor, equivalently the Legendre transform of the Lagrangian of the preceding section. Substituting the mode expansion and using the parent article's normalizations and orthogonality gives the standard result

$$ \hat H=\sum_{r}\int\!\frac{d^3p}{(2\pi)^3}\,E_{\mathbf p} \Big(\hat a_r^\dagger(\mathbf p)\hat a_r(\mathbf p) -\hat b_r(\mathbf p)\hat b_r^\dagger(\mathbf p)\Big). $$

The negative sign on the second term is the operator trace of the Dirac sea. Using the anticommutator $\hat b_r(\mathbf p)\hat b_r^\dagger(\mathbf p)=e_0-\hat b_r^\dagger(\mathbf p)\hat b_r(\mathbf p)$ (more precisely, its continuum form) this becomes

$$ \hat H=\;:\!\hat H\!:\;+\;E_0, $$

where the normal-ordered Hamiltonian is

$$ :\!\hat H\!:\;=\sum_{r}\int\!\frac{d^3p}{(2\pi)^3}\,E_{\mathbf p} \Big(\hat a_r^\dagger(\mathbf p)\hat a_r(\mathbf p) +\hat b_r^\dagger(\mathbf p)\hat b_r(\mathbf p)\Big), $$

manifestly non-negative, and the normal-ordering constant is the c-number

$$ E_0=-\sum_{r}\int\!\frac{d^3p}{(2\pi)^3}\,(2\pi)^3\delta^{(3)}(\mathbf 0)\,E_{\mathbf p} =-\,2V\!\int\!\frac{d^3p}{(2\pi)^3}\,E_{\mathbf p}. $$

Here $V$ is the volume in which the theory is placed, and $(2\pi)^3\delta^{(3)}(\mathbf 0)=V$ is the standard transcription of the continuum normalization. The constant is negative and quartically divergent, in contrast with the positive half-integer zero-point energy of a scalar field, and it is a c-number: it commutes with everything and does not affect the operator equations of motion. Normal ordering is the convention that removes it. For fermions the normal-ordering prescription moves every creation operator to the left of every annihilation operator and inserts a minus sign for each transposition of two fermion operators; for example $:\!\hat b_r(\mathbf p)\hat b_r^\dagger(\mathbf q)\!:=-\hat b_r^\dagger(\mathbf q)\hat b_r(\mathbf p)$. With this definition the constant term in $\hat H$ is removed by construction.

The same mode algebra gives the other conserved charges. The normal-ordered four-momentum is

$$ :\!\hat P^\mu\!:\;=\sum_{r}\int\!\frac{d^3p}{(2\pi)^3}\,p^\mu \Big(\hat a_r^\dagger\hat a_r+\hat b_r^\dagger\hat b_r\Big),\qquad p^\mu=(E_{\mathbf p},\mathbf p), $$

and the normal-ordered electric charge, from $j^\mu=\bar\psi\gamma^\mu\psi$ with $j^0=\psi^\dagger\psi$, is

$$ :\!\hat Q\!:\;=\sum_{r}\int\!\frac{d^3p}{(2\pi)^3} \Big(\hat a_r^\dagger\hat a_r-\hat b_r^\dagger\hat b_r\Big), $$

so that particles and antiparticles carry opposite charge, as they must. All three expressions are the standard ones; their derivation uses only the parent article's spinor relations.

The status of $E_0$ deserves one cautionary sentence. Its divergence is an ultraviolet feature of the free field, removed by the conventional normal-ordering prescription, and nothing in the biquaternion structure as used here changes that. Whether the biquaternion framework offers a distinguished or geometrically meaningful regularization of the vacuum energy is an open question, and we do not claim one.

Spin and Statistics

The anticommutation relations are not a convention; they are selected. The three reasons given in Section 3 — boundedness of the energy, microcausality, and positivity of the norm — are all consequences of one structural fact: the one-particle states of the field transform in the spinor representation of the Lorentz group, and the biquaternion algebra realizes that representation on its module. For a field of half-integer spin the representation is double-valued: a rotation by $2\pi$ in the vector representation corresponds to a rotation by $2\pi$ in the spinor representation, and only a rotation by $4\pi$ returns the spinor to itself. It is this two-valuedness, together with the requirement that the local observables be built from bilinears, that forces the anticommutator: the commutator of two spinor fields fails to vanish at spacelike separation, while the anticommutator does.

Equivalently, the operator algebra generated by the modes is $\mathbb{Z}/2$-graded by fermion parity. Define

$$ (-1)^F\,\hat a_r(\mathbf p)\,(-1)^F=-\hat a_r(\mathbf p),\qquad (-1)^F\,\hat b_r(\mathbf p)\,(-1)^F=-\hat b_r(\mathbf p),\qquad (-1)^F|0\rangle=|0\rangle , $$

so that $(-1)^F$ is $+1$ on states with an even number of particles and antiparticles and $-1$ on states with an odd number. The operators split into an even (bosonic) part commuting with $(-1)^F$ and an odd (fermionic) part anticommuting with it, and the rule that fermionic operators anticommute is precisely the statement that the algebra is graded-commutative. This grading is the structure the KMS article identifies as absent from the biquaternion algebra; the present construction supplies it on the operator algebra generated by the mode operators, which is a partial answer.

This article deliberately stops here. The full spin–statistics theorem — that integer-spin fields must be quantized with commutators and half-integer-spin fields with anticommutators, and the precise relation of that dichotomy to the topology of the Lorentz group — is the subject of a dedicated companion article. We have recorded only the reasons that are operative in the canonical construction, and we have not attempted to reproduce the theorem. In particular, the question of whether the biquaternion framework yields a derivation of spin–statistics rather than a transcription of the standard argument is not settled here.

Relation to the KMS Thermal Framework

The KMS article characterizes thermal equilibrium by an analytic property of correlation functions. For operators $\hat A,\hat B$ in the algebra and a state $\omega_\beta$ at inverse temperature $\beta=\hbar/(k_BT)$, the thermal correlation function is

$$ F_{\hat A\hat B}(t)=\omega_\beta\big(\hat A\,\alpha_t(\hat B)\big), $$

with $\alpha_t$ the Heisenberg evolution; the KMS condition states that $F_{\hat A\hat B}(t)$ extends analytically to the strip $0<\mathrm{Im}(t)<\beta$ and satisfies

$$ F_{\hat A\hat B}(t+i\beta)=F_{\hat B\hat A}(-t) $$

on its boundary. The KMS article observes that this boundary relation is the same for bosons and fermions, that the statistics enters instead through the algebra of the operators, and that a fermionic reading requires a $\mathbb{Z}/2$ grading together with the rule that fermionic operators anticommute — a rule the biquaternion algebra does not automatically supply.

The quantized Dirac field supplies exactly that rule. The mode algebra of Section 5 is graded by $(-1)^F$, and the fermionic operators anticommute. The KMS condition can therefore be stated for the thermal state

$$ \omega_\beta(\,\cdot\,)=\frac{1}{Z}\,\mathrm{Tr}\big(e^{-\beta(\hat H-\mu\hat Q)}\,\cdot\,\big) $$

on the Fock space, with $\hat H$ and $\hat Q$ the normal-ordered operators above, $Z$ the partition function, and $\mu$ the chemical potential. The point of contact with the KMS article is sharp and elementary. Take a single mode of energy $E$ and normalize it so that $\{\hat a,\hat a^\dagger\}=e_0$. In the Heisenberg picture $\alpha_t(\hat a)=e^{-iEt}\hat a$ and $\alpha_t(\hat a^\dagger)=e^{+iEt}\hat a^\dagger$, so

$$ F_{\hat a\hat a^\dagger}(t)=e^{+iEt}\,\omega_\beta(\hat a\hat a^\dagger),\qquad F_{\hat a^\dagger\hat a}(-t)=e^{+iEt}\,\omega_\beta(\hat a^\dagger\hat a). $$

The KMS condition $F_{\hat a\hat a^\dagger}(t+i\beta)=F_{\hat a^\dagger\hat a}(-t)$ then gives $e^{-\beta E}\,\omega_\beta(\hat a\hat a^\dagger)=\omega_\beta(\hat a^\dagger\hat a)$. Writing $n=\omega_\beta(\hat a^\dagger\hat a)$ and using the anticommutator in the form $\hat a\hat a^\dagger=e_0-\hat a^\dagger\hat a$ yields $1-n=e^{\beta E}n$, i.e.

$$ n=\frac{1}{e^{\beta E}+1}, $$

the Fermi–Dirac distribution. The same computation for a bosonic mode, where the relevant relation is the commutator $\hat a\hat a^\dagger=e_0+\hat a^\dagger\hat a$, gives $n=1/(e^{\beta E}-1)$, the Bose–Einstein distribution. This is the cleanest sense in which the statistics is carried by the (anti)commutator and not by the KMS boundary relation, exactly as the KMS article states. In the continuum, the free-field thermal two-point function is

$$ \omega_\beta\big(\hat a_r^\dagger(\mathbf p)\hat a_s(\mathbf q)\big) =(2\pi)^3\,\delta_{rs}\,\delta^{(3)}(\mathbf p-\mathbf q)\,\frac{1}{e^{\beta E_{\mathbf p}}+1}, $$

and the time-ordered fermionic correlation functions built from it are antiperiodic in imaginary time with period $\beta$, hence periodic with period $2\beta$; the Fourier modes of the compact imaginary-time direction are consequently the fermionic Matsubara frequencies

$$ \omega_n=\frac{(2n+1)\pi}{\beta}, $$

in the KMS article's notation. The antiperiodicity is the fermionic twist that the KMS article ties to the spin structure of the field and to the missing grading; here it follows from the same anticommutator that produced the Fermi–Dirac distribution.

The structural facts the KMS article records remain intact. The imaginary time direction is intrinsic to the material sector $\mathbb{M}_-$; the complexified time direction in which the KMS continuation takes place is the sum of the imaginary time direction of $\mathbb{M}_-$ and the real time direction of $\mathbb{M}_+$; and the modular Hamiltonian $K=-\log\rho$ is a Hermitian element of the algebra, hence lies in the informational sector $\mathbb{M}_+$. What the present article adds is the operator algebra on which the thermal state is defined, together with the grading the KMS article flagged as absent. What it does not add is a genuinely biquaternionic construction of the thermal state: the state $\omega_\beta$ above is a complex-valued functional on the Fock space built from the spinor-module field, and whether it extends to a functional intrinsically defined on $\mathbb{B}$ — with the grading embedded in $\mathbb{B}$ rather than in the mode algebra, and with a spinor-module form of the fermionic twist — remains open. The relation established here is a consistency of structures, not a derivation of the fermionic KMS condition from the biquaternion algebra.

What Is Standard and What Is Open

Because the question addressed here is itself an open question of the framework, we separate the two kinds of statement explicitly.

Standard quantum field theory, transcribed.

The Lagrangian $\mathcal{L}=\bar\psi(i\gamma^\mu\partial_\mu-m)\psi$ and the conjugate momentum $\pi=i\psi^\dagger$, together with the second-class constraint $\pi-i\psi^\dagger\approx0$.
The equal-time anticommutation relations, imposed and verified, and the three reasons anticommutation rather than commutation is forced.
The mode expansion in the parent article's plane-wave solutions, with the covariant normalization $\bar u u=2m$, $u^\dagger u=2E$, and the parent's spin sums.
The mode-operator anticommutators, the Fock space, and Pauli exclusion as the algebraic identity $(\hat a^\dagger)^2=0$.
The Hamiltonian $\hat H=\sum_r\int\frac{d^3p}{(2\pi)^3}E_{\mathbf p}(\hat a^\dagger\hat a+\hat b^\dagger\hat b)$ after normal ordering, the normal-ordering constant $E_0=-2V\int\frac{d^3p}{(2\pi)^3}E_{\mathbf p}$, and the four-momentum and charge.
The Fermi–Dirac distribution as the consequence of the KMS condition together with the anticommutator, and the fermionic Matsubara frequencies.

Open in the biquaternion framework.

The intrinsic operator field. Whether the quantized field should be $\mathbb{B}$-valued, $M_2(\mathbb{C})$-valued, or spinor-module-valued, and to what extent these give inequivalent quantizations.
The intrinsic Lagrangian and momentum. The $\mathbb{B}$-intrinsic Lagrangian, its real-scalar pairing, and the conjugate momentum are not fixed by the algebra alone; the spinor-module transcription is the one adopted above.
The biquaternion Fock space. Whether it is the exterior algebra of the spinor-module solution space, or a construction native to $\mathbb{B}$; the planned companion article on the Fock space is the natural place to settle this.
The grading. The $\mathbb{Z}/2$ grading is supplied here on the mode algebra; embedding it in $\mathbb{B}$ itself, as the KMS article requires, is not done.
The spin–statistics theorem. Whether the framework derives spin–statistics or only transcribes it is left to the dedicated article.
The vacuum energy. Whether the biquaternion structure singles out a regularization of $E_0$, or assigns it a geometric meaning, is open.
Empirical content. As with the rest of the framework, whether any of this yields a prediction distinguishing it from standard quantum field theory is the unanswered question.

The extension attempted here is therefore not canonical in the framework: the part that works is the standard canonical quantization of the Dirac field, expressed in the parent article's plane-wave solutions and the material- and informational-sector notation; the part that is genuinely biquaternionic is not yet fixed. This article is the first attempt at the companion article's second-quantization question, and it should be read as delimiting the problem rather than closing it.

Summary

The biquaternion Dirac field is quantized by promoting its spinor-module representative to an operator-valued field. The classical theory has Lagrangian $\mathcal{L}=\bar\psi(i\gamma^\mu\partial_\mu-m)\psi$ and conjugate momentum $\pi=i\psi^\dagger$; the momentum is algebraically tied to the field by a second-class constraint, so the canonical structure is imposed rather than assumed. The equal-time relations are anticommutators,

$$ \{\hat\psi_a(\mathbf x,t),\hat\psi_b^\dagger(\mathbf y,t)\}=\delta_{ab}\delta^{(3)}(\mathbf x-\mathbf y), $$

with all other equal-time anticommutators vanishing; anticommutation rather than commutation is forced by the boundedness of the energy, by microcausality, and by positivity of the norm, all three being consequences of the spinor character of the field.

The field expands in the plane-wave solutions established by the parent article,

$$ \hat{\psi}(x)=\int\!\frac{d^3p}{(2\pi)^3}\frac{1}{\sqrt{2E_{\mathbf p}}}\sum_r \Big[\hat a_r(\mathbf p)u^{(r)}(\mathbf p)e^{-ip\cdot x}+\hat b_r^\dagger(\mathbf p)v^{(r)}(\mathbf p)e^{+ip\cdot x}\Big], $$

with the mode anticommutators $\{\hat a_r(\mathbf p),\hat a_s^\dagger(\mathbf q)\}=\{\hat b_r(\mathbf p),\hat b_s^\dagger(\mathbf q)\}=(2\pi)^3\delta_{rs}\delta^{(3)}(\mathbf p-\mathbf q)$ and all others zero. These reproduce the equal-time relation: the positive- and negative-frequency branches enter with opposite phases, the even part of the spin sums integrates to the delta function, and the odd parts cancel. The resulting Fock space carries Pauli exclusion as the algebraic identity $(\hat a_r^\dagger)^2=0$.

The Hamiltonian is $\hat H=\sum_r\int\frac{d^3p}{(2\pi)^3}E_{\mathbf p}(\hat a_r^\dagger\hat a_r+\hat b_r^\dagger\hat b_r)$ after normal ordering, with the negative quartically divergent normal-ordering constant $E_0=-2V\int\frac{d^3p}{(2\pi)^3}E_{\mathbf p}$ removed. The normal-ordered charge has particles and antiparticles of opposite sign. The mode algebra is $\mathbb{Z}/2$-graded by fermion parity, which is the structure the KMS article identifies as missing; with it, the KMS condition plus the anticommutator yields the Fermi–Dirac distribution and the antiperiodic imaginary-time correlators whose Matsubara frequencies are $(2n+1)\pi/\beta$.

The extension is a transcription of standard canonical quantization onto the parent article's biquaternion plane waves, not a derivation of the quantized structure from the algebra $\mathbb{B}$. The intrinsic $\mathbb{B}$-valued Lagrangian and momentum, the biquaternion Fock space, and the embedding of the grading in $\mathbb{B}$ remain open.

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$	Scalar imaginary, $i^2=-1$
$\mathbb{M}_+,\mathbb{M}_-$	Informational (Hermitian) and material (anti-Hermitian) sectors
$\tilde{\nabla}=e_0\partial_{ict}+\sum_k e_k\partial_k$	Biquaternionic gradient (Dirac operator)
$\Box=\tilde{\nabla}\tilde{\nabla}^{\natural}$	d'Alembertian
$\tilde{\Psi}$	Classical biquaternion Dirac field
$\tilde{\Psi}^\flat=-\tilde{\Psi}^{*}$	Anti-Hermitian conjugate (the algebra's real structure; not the mass term)

| $\psi=(\psi_L,\psi_R)$, $\bar\psi=\psi^\dagger\gamma^0$ | Dirac spinor (spinor module) and its adjoint | | $\gamma^\mu$, $\gamma_5=i\gamma^0\gamma^1\gamma^2\gamma^3$ | Gamma matrices and chirality operator | | $g=\mathrm{diag}(+1,-1,-1,-1)$ | Clifford metric of the $\gamma^\mu$; $g=-\eta$ | | $\eta=\mathrm{diag}(-1,+1,+1,+1)$ | Spacetime metric of the $ict$ gradient | | $\mathcal{L}=\bar\psi(i\gamma^\mu\partial_\mu-m)\psi$ | Dirac Lagrangian density | | $\pi=\partial\mathcal{L}/\partial(\partial_0\psi)=i\psi^\dagger$ | Conjugate momentum; $\pi-i\psi^\dagger\approx0$ is second class | | $u^{(r)}(\mathbf p),v^{(r)}(\mathbf p)$ | Parent article's positive- and negative-frequency spinors | | $\bar u u=2m$, $u^\dagger u=2E_{\mathbf p}$ | Covariant and Hermitian normalizations | | $p^\mu=(E_{\mathbf p},\mathbf p)$, $\not p=\gamma^0E_{\mathbf p}-\boldsymbol\gamma\cdot\mathbf p$ | Four-momentum and Feynman slash | | $\tilde k\tilde k^{\natural}=-m^2c^2/\hbar^2$ | Biquaternionic mass-shell condition | | $\hat\psi$, $\hat a_r(\mathbf p)$, $\hat b_r(\mathbf p)$ | Quantized field; particle and antiparticle annihilation operators | | $\{\hat a_r(\mathbf p),\hat a_s^\dagger(\mathbf q)\}=(2\pi)^3\delta_{rs}\delta^{(3)}(\mathbf p-\mathbf q)$ | Mode anticommutator | | $\hat H$, $:\!\hat H\!:$ | Hamiltonian and its normal-ordered form | | $E_0=-2V\int\frac{d^3p}{(2\pi)^3}E_{\mathbf p}$ | Normal-ordering constant (vacuum energy) | | $(-1)^F$ | Fermion-parity grading | | $\beta=\hbar/(k_BT)$, $\omega_\beta$, $\alpha_t$ | Inverse temperature, thermal state, Heisenberg evolution | | $F_{\hat A\hat B}(t+i\beta)=F_{\hat B\hat A}(-t)$ | KMS condition | | $K=-\log\rho$ | Modular Hamiltonian (Hermitian, in $\mathbb{M}_+$) | | $\omega_n=(2n+1)\pi/\beta$ | Fermionic Matsubara frequencies |