A discrete version of plane wave solutions of the Dirac equation in the Joyce form

TL;DR

This paper develops a discrete geometric scheme for plane wave solutions of the Dirac-Kähler equation in Joyce form, analyzing conditions for their validity in a discretized setting.

Contribution

It introduces a novel discrete discretization scheme combining forward and backward differences for the Dirac-Kähler equation in Joyce form.

Findings

A geometric discretization scheme for the Dirac-Kähler equation is proposed.

Conditions for discrete plane wave solutions to satisfy the discrete Joyce equation are identified.

The scheme bridges continuous and discrete formulations of Dirac equations.

Abstract

We construct a discrete version of the plane wave solution to a discrete Dirac-K\"{a}hler equation in the Joyce form. A geometric discretisation scheme based on both forward and backward difference operators is used. The conditions under which a discrete plane wave solution satisfies a discrete Joyce equation are discussed.

Equations

$$(\omega ,\phi )={\int }_M\omega \wedge \ast \overline{\phi },$$

$$\delta =\ast d\ast .$$

$$\Lambda (M)={\Lambda }^0(M)\oplus {\Lambda }^1(M)\oplus {\Lambda }^2(M)\oplus {\Lambda }^3(M)\oplus {\Lambda }^4(M).$$

$$\Omega =\sum _{r=0}^4\stackrel{r}{\omega },$$

$$i(d+\delta )\Omega =m\Omega ,$$

$$i\delta \stackrel{1}{\omega }=m\stackrel{0}{\omega },$$

$$i(d\stackrel{0}{\omega }+\delta \stackrel{2}{\omega })=m\stackrel{1}{\omega },$$

$$i(d\stackrel{1}{\omega }+\delta \stackrel{3}{\omega })=m\stackrel{2}{\omega },$$

$$i(d\stackrel{2}{\omega }+\delta \stackrel{4}{\omega })=m\stackrel{3}{\omega },$$

$$id\stackrel{3}{\omega }=m\stackrel{4}{\omega }.$$

$$\nabla =\sum _{\mu =0}^3{\gamma }_{\mu }{\partial }^{\mu },\mu =0,1,2,3,$$

$$i\nabla \Omega =m\Omega ,\Omega \in \text{C}\ell (1,3).$$

$$i\nabla {\Omega }^{ev}=m{\Omega }^{ev}{\gamma }_0,{\Omega }^{ev}\in \text{C}{\ell }^{ev}(1,3).$$

$$\Phi =A{e}^{\pm ip\cdot x},$$

$$i(d+\delta ){\Omega }^{ev}=m{\Omega }^{ev}{e}^0,{\Omega }^{ev}\in {\Lambda }^{ev}(M).$$

$$\partial x_{\kappa }=0,\partial e_{\kappa }=x_{\kappa +1}-x_{\kappa }.$$

$$s_k=s_{k_0}\otimes s_{k_1}\otimes s_{k_2}\otimes s_{k_3},$$

$$e_k^{01}=e_{k_0}\otimes e_{k_1}\otimes x_{k_2}\otimes x_{k_3},e_k^{12}=x_{k_0}\otimes e_{k_1}\otimes e_{k_2}\otimes x_{k_3},$$

$$e_k^{02}=e_{k_0}\otimes x_{k_1}\otimes e_{k_2}\otimes x_{k_3},e_k^{13}=x_{k_0}\otimes e_{k_1}\otimes x_{k_2}\otimes e_{k_3},$$

$$e_k^{03}=e_{k_0}\otimes x_{k_1}\otimes x_{k_2}\otimes e_{k_3},e_k^{23}=x_{k_0}\otimes x_{k_1}\otimes e_{k_2}\otimes e_{k_3}.$$

$$\partial (s_{k_{\mu }}\otimes s_k)=\partial s_{k_{\mu }}\otimes s_k+Q(k_{\mu })s_{k_{\mu }}\otimes \partial s_k,$$

$$K(4)=K^0(4)\oplus K^1(4)\oplus K^2(4)\oplus K^3(4)\oplus K^4(4).$$

$$\stackrel{r}{\omega }=\sum _k\sum _{(r)}{\omega }_k^{(r)}{s}_{(r)}^k,$$

$$\langle\varepsilon_{k},\ s^{k}\rangle=\left\{\begin{array}[]{l}0,\quad\varepsilon_{k}\neq s_{k}\\ 1,\quad\varepsilon_{k}=s_{k}.\end{array}\right.$$

$$⟨\partial a,\stackrel{r}{\omega }⟩=⟨a,d^c\stackrel{r}{\omega }⟩,$$

$$d^c\stackrel{4}{\omega }=0\text{mbox}and{d}^c{d}^c\stackrel{r}{\omega }=0\text{mbox}foranyr.$$

$${\tau }_{\mu }k=(k_0,...k_{\mu }+1,...k_3),{\sigma }_{\mu }k=(k_0,...k_{\mu }-1,...k_3),\mu =0,1,2,3.$$

$${\Delta }_{\mu }^{+}{\omega }_k^{(r)}={\omega }_{{\tau }_{\mu }k}^{(r)}-{\omega }_k^{(r)},$$

$${\Delta }_{\mu }^{-}{\omega }_k^{(r)}={\omega }_k^{(r)}-{\omega }_{{\sigma }_{\mu }k}^{(r)},$$

$${\Delta }_{\mu }^{-}{\omega }_k^{(r)}={\Delta }_{\mu }^{+}{\omega }_{{\sigma }_{\mu }k}^{(r)},{\Delta }_{\mu }^{+}{\omega }_k^{(r)}={\Delta }_{\mu }^{-}{\omega }_{{\tau }_{\mu }k}^{(r)}.$$

$$\stackrel{0}{\omega }=\sum _k{\stackrel{0}{\omega }}_k{x}^k,\stackrel{2}{\omega }=\sum _k\sum _{\mu <\nu }{\omega }_k^{\mu \nu }{e}_{\mu \nu }^k,\stackrel{4}{\omega }=\sum _k{\stackrel{4}{\omega }}_k{e}^k,$$

$$\stackrel{1}{\omega }=\sum _k\sum _{\mu =0}^3{\omega }_k^{\mu }{e}_{\mu }^k,\stackrel{3}{\omega }=\sum _k\sum _{\iota <\mu <\nu }{\omega }_k^{\iota \mu \nu }{e}_{\iota \mu \nu }^k,$$

$$x^k=x^{k_0}\otimes x^{k_1}\otimes x^{k_2}\otimes x^{k_3},e^k=e^{k_0}\otimes e^{k_1}\otimes e^{k_2}\otimes e^{k_3}$$

Full text

**Key words and phrases: **geometric discretisation, Dirac-Kähler equation, Joyce equation, difference equations, discrete models, plane wave solution

**Math. Subj. Clas.: **39A12, 39A70, 81Q05

1 Introduction

Discrete models od Dirac type equations based on the Dirac-Kähler formulation have been of interest recently from both the applied and the theoretical point of view. In the Dirac-Kähler approach, a discretisation scheme is geometric in nature and rests upon the use of the differential forms calculus. This means that the geometric properties and the algebraic relationships between the differential, the exterior product and the Hodge star operator are expected to be captured in the case of their discrete counterparts. This work is a continuation of that studied in the papers [15, 16, 17, 18]. In this paper, we are mainly interested in a discrete model in which both forward and backward differences operators are used. This model relies on the use a geometric discretisation scheme proposed in [15], and a discrete Clifford calculus framework on discrete forms described in [18]. For a review of discrete Clifford calculus frameworks on lattices, we refer the reader to [3, 5, 6, 7, 13, 19]. Our purpose here is to construct a discrete version of the plane wave solution to a discrete Dirac-Kähler equation in the Joyce form.

We first briefly review some definitions and basic facts on the Dirac-Kähler equation [12, 14] and the Dirac equation in the spacetime algebra [8, 9]. Let $M={\mathbb{R}}^{1,3}$ be Minkowski space. Denote by $\Lambda^{r}(M)$ the vector space of smooth differential $r$-forms, $r=0,1,2,3,4$. We consider $\Lambda^{r}(M)$ over $\mathbb{C}$. Let $\omega,\ \varphi\in\Lambda^{r}(M)$. The inner product is defined by

[TABLE]

where $\wedge$ is the exterior product, $\overline{\varphi}$ denotes the complex conjugate of the form $\varphi$ and $\ast$ is the Hodge star operator $\ast:\Lambda^{r}(M)\rightarrow\Lambda^{4-r}(M)$ with respect to the Lorentz metric. Let $d:\Lambda^{r}(M)\rightarrow\Lambda^{r+1}(M)$ be the exterior differential and let $\delta:\Lambda^{r}(M)\rightarrow\Lambda^{r-1}(M)$ be the formal adjoint of $d$ with respect to (1.1). We have

[TABLE]

Denote by $\Lambda(M)$ the set of all differential forms on $M$. We have

[TABLE]

Let $\Omega\in\Lambda(M)$ be an inhomogeneous differential form, i.e.

[TABLE]

where $\overset{r}{\omega}\in\Lambda^{r}(M)$. The Dirac-Kähler equation for a free electron is given by

[TABLE]

where $i$ is the usual complex unit and $m$ is a mass parameter. It is easy to show that Eq. (1.3) is equivalent to the set of equations

[TABLE]

The operator $d+\delta$ is the analogue of the gradient operator in Minkowski spacetime

[TABLE]

where $\gamma_{\mu}$ is the Dirac gamma matrix. If one think of $\{\gamma_{0},\gamma_{1},\gamma_{2},\gamma_{3}\}$ as a vector basis in spacetime, then the gamma matrices $\gamma_{\mu}$ can be considered as generators of the Clifford algebra of spacetime $\mathcal{C}\ell(1,3)$ [1, 2]. Hestenes [9] calls this algebra the spacetime algebra. It is known that an inhomogeneous form $\Omega$ can be represented as element of $\mathcal{C}\ell(1,3)$. Then the Dirac-Kähler equation can be written as the algebraic equation

[TABLE]

Equation (1.5) is equivalent to the four Dirac equations (traditional column-spinor equations) for a free electron. Let $\emph{C}\ell^{ev}(1,3)$ be the even subalgebra of the algebra $\emph{C}\ell(1,3)$. Consider the equation

[TABLE]

In [10], this equation is called the ”generalized bivector Dirac equation”. Following Baylis [2] we call Eq. (1.6) the Joyce equation. This equation is equivalent to two copies of the usual Dirac equation. For a deeper discussion of equivalence of Dirac formulations we refer the reader to [11]. Equation (1.6) admits the plane wave solution of the form

[TABLE]

where $A\in\mathcal{C}\ell^{ev}(1,3)$ is a constant element, $p=\{p_{0},p_{1},p_{2},p_{3}\}$ is a four-momentum and $p\cdot x=p_{\mu}x^{\mu}$ (see [2] for more details).

It should be noted that the graded algebra $\Lambda(M)$ endowed with the Clifford multiplication is an example of a Clifford algebra. In this case the basis covectors $e^{\mu}=dx^{\mu}$ are considered as generators of the Clifford algebra. Let $\Lambda^{ev}(M)=\Lambda^{0}(M)\oplus\Lambda^{2}(M)\oplus\Lambda^{4}(M)$. Denote by $\Omega^{ev}$ the even part of the form (1.2). Then Eq. (1.6) can be rewritten in terms of inhomogeneous forms as

[TABLE]

We call this equation the Dirac-Kähler equation in the Joyce form.

In this paper, we focus on the construction of a discrete version of the plane wave solution (1.7) for a discrete counterpart of Eq. (1.8). In [17], the same problem was considered by using a discretisation scheme based on a combinatorial double complex construction. However, this construction generates difference operators of the forward type only and it is not enough to obtain an exact geometric counterpart of the Clifford algebra $\emph{C}\ell(1,3)$. For this purpose, both forward and backward differences are needed. This fact is well-known [14]. As a new result, we announce here the construction of a discrete plane wave solution based on both forward and backward difference operators. This seems to suggest that such a discrete model will be more suitable for applied purposes.

2 Combinatorial model and difference operators

A combinatorial model of Minkowski space and a discretisation scheme are adopted from [16]. For the convenience of the reader we briefly repeat the relevant material from [16] without proofs, thus making our presentation self-contained. Following [4], let the tensor product $C(4)=C\otimes C\otimes C\otimes C$ of a 1-dimensional complex be a combinatorial model of Euclidean space ${\mathbb{R}}^{4}$. The 1-dimensional complex $C$ is defined in the following way. Introduce the sets $\{x_{\kappa}\}$ and $\{e_{\kappa}\}$, $\kappa\in{\mathbb{Z}}$. Let $C^{0}$ and $C^{1}$ be the free abelian groups of 0-dimensional and 1-dimensional chains generated by $\{x_{\kappa}\}$ and $\{e_{\kappa}\}$. The free abelian group is understood as the direct sum of infinity cyclic groups generated by $\{x_{\kappa}\}$, $\{e_{\kappa}\}$. The boundary operator $\partial:C^{0}\rightarrow 0,\quad\partial:C^{1}\rightarrow C^{0}$ is given by

[TABLE]

The definition is extended to arbitrary chains by linearity. The direct sum $C=C^{0}\oplus C^{1}$ with the boundary operator $\partial$ defines the 1-dimensional complex. It is known that a free abelian group is an abelian group with basis. One can regard the sets $\{x_{\kappa}\}$, $\{e_{\kappa}\}$ as sets of basis elements of the groups $C^{0}$ and $C^{1}$. Geometrically we can interpret the 0-dimensional basis elements $x_{\kappa}$ as points of the real line and the 1-dimensional basis elements $e_{\kappa}$ as open intervals between points, i.e. $e_{\kappa}=(x_{\kappa},x_{\kappa+1})$. We call the complex $C$ a combinatorial real line.

Multiplying the basis elements $x_{\kappa}$, $e_{\kappa}$ in various way we obtain basis elements of $C(4)$. Let $s_{k}$ be an arbitrary basis element of $C(4)$. Then we have

[TABLE]

where $s_{k_{\mu}}$ is either $x_{k_{\mu}}$ or $e_{k_{\mu}}$ and $k=(k_{0},k_{1},k_{2},k_{3})$, $k_{\mu}\in\mathbb{Z}$. Let us denote by $s_{k}^{(r)}$ the $r$-dimensional basis element of $C(4)$. The dimension $r$ of $s_{k}^{(r)}$ is given by the number of factors $e_{k_{\mu}}$ that appear in it. The notation $s_{k}^{(r)}$ means that the product (2.2) contains exactly $r$ $1$-dimensional elements $e_{k_{\mu}}$ and $4-r$ [math]-dimensional elements $x_{k_{\mu}}$, and the superscript $(r)$ indicates also a position of $e_{k_{\mu}}$ in $s_{k}^{(r)}$. For example, the 2-dimensional basis elements of $C(4)$ can be written as

[TABLE]

Let $C(p)$ is the tensor product of $p$ factors of $C=C(1)$, $p=1,2,3$. The definition (2.1) of $\partial$ is extended to an arbitrary basis element of $C(4)$ by induction on $p$. Suppose that the boundary operator has been defined for any basis element $s_{k}\in C(p)$. Then we introduce it for the basis element $s_{k_{\mu}}\otimes s_{k}\in C(p+1)$ by the rule

[TABLE]

where $s_{k_{\mu}}\in C$ and $Q(k_{\mu})$ is equal to $+1$ if $s_{k_{\mu}}=x_{k_{\mu}}$ and to $-1$ if $s_{k_{\mu}}=e_{k_{\mu}}$. The operation (2.3) is linearly extended to arbitrary chains. It is easy to check that $\partial\partial a=0$ for any chain $a\in C(4)$.

Let $C(4)$ be a combinatorial model of Minkowski space. In what follows we assume that the index $k_{0}$ in (2.2) corresponds to the time coordinate of $M$. Hence, the basis elements of the indicated 1-dimensional complex $C$ will be written as $x_{k_{0}}$ and $e_{k_{0}}$.

Let us now consider a dual complex to $C(4)$. We define it as the complex of cochains $K(4)$ with complex coefficients. The complex $K(4)$ has the same structure as $C(4)$, namely ${K(4)=K\otimes K\otimes K\otimes K}$, where $K$ is a dual complex to the 1-dimensional complex $C$. We will use a superscript to indicate a basis element of $K$. Let $x^{\kappa}$ and $e^{\kappa}$, $\kappa\in{\mathbb{Z}}$, be the 0- and 1-dimensional basis elements of $K$. Then an arbitrary $r$-dimensional basis element of $K(4)$ can be written as $s_{(r)}^{k}=s^{k_{0}}\otimes s^{k_{1}}\otimes s^{k_{2}}\otimes s^{k_{3}}$, where $s^{k_{\mu}}$ is either $x^{k_{\mu}}$ or $e^{k_{\mu}}$ and $k=(k_{0},k_{1},k_{2},k_{3})$. We will call cochains forms, emphasizing their relationship with differential forms. Denote by $K^{r}(4)$ the set of all $r$-forms. Then $K(4)$ can be expressed by

[TABLE]

The complex $K(4)$ is a discrete analogue of $\Lambda(M)$. Let $\overset{r}{\omega}\in K^{r}(4)$, then we have

[TABLE]

where $\omega_{k}^{(r)}\in\mathbb{C}$.

As in [4], we define the pairing (chain-cochain) operation for any basis elements $\varepsilon_{k}\in C(4)$, $s^{k}\in K(4)$ by the rule

[TABLE]

The operation (2.5) is linearly extended to arbitrary chains and cochains.

The coboundary operator $d^{c}:K^{r}(4)\rightarrow K^{r+1}(4)$ is defined by

[TABLE]

where $a\in C(4)$ is an $r+1$ dimensional chain. The operator $d^{c}$ is an analog of the exterior differential. From the above it follows that

[TABLE]

It is convenient to introduce the shift operators $\tau_{\mu}$ and $\sigma_{\mu}$ in the set of indices by

[TABLE]

Let the difference operators $\Delta^{+}_{\mu}$ and $\Delta^{-}_{\mu}$ be defined by

[TABLE]

[TABLE]

where $\omega_{k}^{(r)}\in\mathbb{C}$ is a component of $\overset{r}{\omega}\in K^{r}(4)$. It is clear, that

[TABLE]

Note that it is enough to use one of the two above defined difference operators to describe discrete analogs of $d$ and $\delta$, as it is shown in [15], but in the proposed approach both forward and backward differences are needed to construct a discrete version of the plane wave solution. It is helpful here to write (2.4) more explicitly as

[TABLE]

[TABLE]

where

[TABLE]

are the 0-, 4-dimensional basis elements of $K(4)$, and $e_{\mu}^{k}$, $e_{\mu\nu}^{k}$ and $e_{\iota\mu\nu}^{k}$ are the 1-, 2- and 3-dimensional basis elements of $K(4)$.

Using (2.3), (2.6) and (2.8) we can calculate

[TABLE]

[TABLE]

[TABLE]

Let us now introduce a $\cup$-multiplication of discrete forms which is an analog of the exterior multiplication for differential forms. Denote by $K(p)$ the tensor product of $p$ factors of the 1-dimensional complex $K=K(1)$. For the basis elements of $K$ the $\cup$-multiplication is defined as follows

[TABLE]

supposing the product to be zero in all other case. To arbitrary basis elements of $K(p)$, $p=2,3,4$, the definition (2.14) is extended by induction on $p$ (see [4] and [15] for more details). Again, to arbitrary discrete forms the $\cup$-multiplication is extended linearly.

It is important to note that the definition above is suitable to deal with a discrete version of the Leibniz rule. The following result was proven in [4].

Proposition 2.1**.**

Let $\varphi\in K^{r}(4)$ and $\psi\in K^{q}(4)$ be arbitrary discrete forms. Then

[TABLE]

Consider a discrete analog of the Hodge star operator $\ast$. For simplicity of notation, we continue to write $\ast$ in the discrete case. Define the operation $\ast:K^{r}(4)\rightarrow K^{4-r}(4)$ for an arbitrary basis element $s^{k}=s^{k_{0}}\otimes s^{k_{1}}\otimes s^{k_{2}}\otimes s^{k_{3}}$ by the rule

[TABLE]

where $Q(k_{0})$ is equal to $+1$ if $s^{k_{0}}=x^{k_{0}}$ and to $-1$ if $s^{k_{0}}=e^{k_{0}}$, and $e^{k}$ is given by (2.10). For example, for the 1-dimensional basis elements $e_{\mu}^{k}$ we have $e_{0}^{k}\cup\ast e_{0}^{k}=-e^{k}$ and $e_{\mu}^{k}\cup\ast e_{\mu}^{k}=e^{k}$  for $\mu=1,2,3$. It is clear that the definition (2.16) preserves the Lorentz signature of metric in our discrete model. From (2.16) a more detailed calculation leads to

[TABLE]

[TABLE]

[TABLE]

[TABLE]

Here $\tau_{\mu\nu}k=\tau_{\mu}\tau_{\nu}k$ and $\tau_{\mu\nu\iota}k=\tau_{\mu}\tau_{\nu}\tau_{\iota}k$, where $\tau_{\mu}$ is defined by (2.7), and

[TABLE]

The operation $\ast$ is linearly extended to arbitrary forms. It is easy to check that

[TABLE]

where $s^{k}_{(r)}$ is an $r$-dimensional basic element of $K(4)$. Then if we perform the operation $\ast$ twice on any $r$-form $\overset{r}{\omega}\in K(4)$, we obtain

[TABLE]

where $\tau k$ is given by (2.21) and

[TABLE]

Hence the operation $(\ast)^{2}$ is equivalent to a shift with corresponding sign. This is slightly different from the continuum case, where applying the Hodge star operator twice leaves a differential form unchanged up to sign, i.e. $(\ast)^{2}=\pm 1$.

Let us consider the $4$-dimensional finite chain $e_{n}\subset C(4)$ of the form:

[TABLE]

where $n_{\mu}\in\mathbb{N}$ is a fixed number for each $\mu=0,1,2,3$ and $e_{k}$ is given by (2.10). The finite sum (2.23) is the set defined by a finite number of $4$-dimensional basis elements of $C(4)$. This set imitates a domain of $M$.

Suppose that the $r$-form (2.4) is vanished on $C(4)\setminus e_{n}$, i.e. if $k_{\mu}<1$ or $k_{\mu}>n_{\mu}$ then $\omega_{k}^{(r)}=0$ for any $r$ and $k$. Then for forms $\varphi,\ \omega\in K^{r}(4)$ of the same degree $r$ the inner product over the set $e_{n}$ (2.23) is defined by the rule

[TABLE]

where $\overline{\omega}$ denotes the complex conjugate of the form $\omega$. For forms of different degrees the product (2.24) is set equal to zero. See also [15]. The definition (2.24) imitates correctly the continuum case (1.1) and the Lorentz metric structure is still captured here. Using (2.5), (2.14) and (2.17)–(2.20) we obtain

[TABLE]

The following result is taken from [15].

Proposition 2.2**.**

Let $\overset{r}{\varphi}\in K^{r}(4)$ and $\overset{r+1}{\omega}\in K^{r+1}(4)$, $r=0,1,2,3$. Then we have

[TABLE]

where

[TABLE]

is the operator formally adjoint of $d^{c}$.

Here $\ast^{-1}$ is the inverse of $\ast$, i.e. $\ast\ast^{-1}=1$.

The operator $\delta^{c}:K^{r+1}(4)\rightarrow K^{r}(4)$ is a discrete analog of the codifferential $\delta$. For the 0-form $\overset{0}{\omega}\in K^{0}(4)$ we have $\delta^{c}\overset{0}{\omega}=0$. It is obvious from (2.25) that $\delta^{c}\delta^{c}\overset{r}{\omega}=0$  for any $r=1,2,3,4.$ Using (2.11)–(2.13) and (2.25) we can calculate

[TABLE]

[TABLE]

[TABLE]

[TABLE]

Note that formulas (2.27)–(2.30) are essentially the same as the corresponding formulas from [15] up to notation.

As in the continuum case, the discrete codifferencial does not depend on the signature of the inner product, i.e. the operator $\delta^{c}$ is the same in both $(+---)$ and $(-+++)$ cases. Recall that in the continuum case we have $\delta=\ast d\ast.$ However in the discrete case the operator $\delta^{c}$ is not equal to $\ast d^{c}\ast$. Indeed, let

[TABLE]

But an easy computation shows that

[TABLE]

where $\sigma k$ is given by (2.22). Therefore in comparison with $\delta^{c}\overset{r+1}{\omega}$, the form $\ast d^{c}\ast\overset{r+1}{\omega}$ has the same components which are shifting to the left by one by all indexes. For example,

[TABLE]

while $\delta^{c}\overset{1}{\omega}$ has the view (2.27).

3 Discrete Dirac-Kähler and Joyce equations

The linear map

[TABLE]

is called a discrete analogue of the Laplacian. It is clear that

[TABLE]

Then a square root of the discrete Laplacian can be written either as $d^{c}-\delta^{c}$ or as $i(d^{c}+\delta^{c})$, where $i$ is the usual complex unit. Here we consider a discrete version of the Dirac-Kähler equation with $i(d^{c}+\delta^{c})$. Note that in [15], a discrete model of the Dirac-Kähler equation has been constructed by using the operator $d^{c}-\delta^{c}$. Let $\Omega$ be a discrete inhomogeneous form, that is

[TABLE]

where $\overset{r}{\omega}$ is given by (2.4). An alternative notation includes $(\Omega)_{r}=\overset{r}{\omega}$ for the $r$-form part of an inhomogeneous form.

Introduce a discrete analog of the Dirac-Kähler equation (1.3) by the rule

[TABLE]

where $m$ is a positive number (mass parameter). We can write this equation more explicitly by separating its homogeneous components as

[TABLE]

Using (2.11)–(2.13) and (2.27)–(2.30) Eqs. (3) can be written as the set of following difference equations

[TABLE]

Note that the equations above contain the difference operators both the forward and backward type. This is in contrast with the situation described in [15], where difference operators of the forward type only are used for the discrete construction.

Let us define the Clifford multiplication in $K(4)$ by the following rules:

[TABLE]

supposing the product to be zero in all other cases.

The operation is linearly extended to arbitrary discrete forms.

Consider the following unit forms

[TABLE]

Note that the unit 0-form $x$ plays a role of the unit element in $K(4)$ with respect to the Clifford multiplication, i.e. for any $r$-form $\overset{r}{\omega}$ we have

[TABLE]

The following is straightforward.

Proposition 3.1**.**

[TABLE]

Proposition 3.2**.**

For any inhomogeneous form $\Omega\in K(4)$ we have

[TABLE]

where $\Delta^{+}_{\mu}$ and $\Delta^{-}_{\mu}$ are the difference operators which act on each component of $\overset{r}{\omega}$ by the rules (2.8) and (2.9).

Proof.

By (2.11) it is clear that

[TABLE]

In the case of the 2-form $\overset{2}{\omega}$ we compute

[TABLE]

From this using (2) and (2) we have

[TABLE]

and

[TABLE]

For $\overset{4}{\omega}$ we obtain

[TABLE]

By (2.30), we have

[TABLE]

Similar calculations apply to the case of odd forms. Hence

[TABLE]

and

[TABLE]

Combining (3.7) and (3.8) yields (3.6). ∎

Clearly, the discrete Dirac-Kähler equation (3.2) can be rewritten in the form

[TABLE]

Let $K^{ev}(4)=K^{0}(4)\oplus K^{2}(4)\oplus K^{4}(4)$ and let $\Omega^{ev}\in K^{ev}(4)$ be a complex-valued even inhomogeneous form, i.e. $\Omega^{ev}=\overset{0}{\omega}+\overset{2}{\omega}+\overset{4}{\omega}$. A discrete analogue of the Joyce equation (1.8) is defined by

[TABLE]

where $e_{0}$ is given by (3.4). From (3.6) it follows that Eq. (3.9) is equivalent to

[TABLE]

A more detailed calculation leads to the following system of 8 difference equations

[TABLE]

for each $k=(k_{0},k_{1},k_{2},k_{3})$.

4 Plane Wave Solutions

In this section, we consider solutions of the discrete Joyce equation which imitate the plane wave solutions for the continuum counterpart. We will mainly follow the strategy of [18], but there are extra difficulties here. These difficulties arise in the construction of eigenvectors of the operator $i(d^{c}+\delta^{c})$ since this operator consists of difference operators of the two types.

Let us consider the complex-valued 0-forms

[TABLE]

where

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

and $p_{\mu}\in\mathbb{R}$. It is easy to check that

[TABLE]

[TABLE]

As a consequence we obtain

[TABLE]

[TABLE]

[TABLE]

Let $A\in K^{ev}(4)$ is a constant complex-valued form. Hence $A$ can be expanded as

[TABLE]

where $\alpha^{0},\alpha^{\mu\nu},\alpha^{4}\in\mathbb{C}$ and $x$, $e_{\mu\nu}$, $e$ are the unit forms given by (3.4). Consider the form

[TABLE]

where

[TABLE]

and $\psi^{0}$, $\psi^{\mu\nu}$ and $\psi^{4}$ are given by (4.1).

Now we apply the operators $d^{c}$ and $\delta^{c}$ to the forms $\overset{0}{\varphi}$, $\overset{2}{\varphi}$ and $\overset{4}{\varphi}$. Using (3.7), (3.8) and (4.2) we obtain

[TABLE]

and

[TABLE]

By (4.3) and (4.4), as in the proof Proposition 3.2 we calculate

[TABLE]

and

[TABLE]

This yields

[TABLE]

Therefore, we have the following result.

Proposition 4.1**.**

For the form (4.6) we have that

[TABLE]

Substituting (4.7) into Eq. (3.9) we obtain

[TABLE]

It is not difficult to show that Eq. (4.8) is equivalent to the following system of equations

[TABLE]

From (3.5) it follows that $e_{0}e_{0}=x$. Hence Eq. (4.8) can be written as

[TABLE]

A direct computation shows that

[TABLE]

Then multiplying both sides of (4.17) by the same factor: -\big{(}p_{0}x-\sum_{\mu=1}^{3}p_{\mu}e_{0}e_{\mu}\big{)} gives

[TABLE]

By (3.5) this gives the equation

[TABLE]

Applying (4.8) to the right-hand side we obtain

[TABLE]

or equivalently,

[TABLE]

Thus exactly as in [18], we have the following assertion.

Proposition 4.2**.**

The form (4.6) is a non-trivial solution of Eq. (3.9) if and only if

[TABLE]

Therefore, if we define $p=\{p_{0},p_{1},p_{2},p_{3}\}$ to be the energy-momentum vector of a particle with (proper) mass $m$, then the relation (4.18) is the energy-momentum relation. In view of this, it makes sense to state that the form (4.6) is a discrete version of the plane wave solution. Recall that the Joyce equation (1.6) admits the plane wave solutions of the form (1.7). Thus the 0-form $\psi=\psi^{0}+\psi^{\mu\nu}+\psi^{4}$ given by (4.1) plays a role of the function $e^{ip\cdot x}$ in the discrete case.

It should be noted that the same proposition for a discrete analog of the Hestenes equation in the case of a discrete model based on the double complex construction is proven in [17].

Let us represent the even complex-valued form (4.5) as

[TABLE]

where

[TABLE]

[TABLE]

It is easy to check that $A_{+}$ commutes with $e_{0}$ and $A_{-}$ anticommutes with it, i.e.

[TABLE]

Then the form (4.6) can be represent also as

[TABLE]

where

[TABLE]

[TABLE]

Since the 0-forms (4.1) commute with $e_{0}$ then Eq. (4.19) is true also for $\Phi_{\pm}$, i.e.

[TABLE]

As is shown in [18], the proof of the following lemma follows from a direct computation.

Lemma 4.3**.**

The form $e_{0\mu}\Phi_{-}$ commutes with $e_{0}$ and has the view (4.20), while $e_{0\mu}\Phi_{+}$ anticommutes with $e_{0}$ and has the view (4.21) for any $\mu=1,2,3$.

By Lemma 4.3, we have the following result whose proof may be found in [18].

Theorem 4.4**.**

The form $\Phi$ given by (4.6) is a non-trivial solution of the discrete Joyce equation if and only if the condition

[TABLE]

holds, or equivalently,

[TABLE]

It is not difficult to show that the condition (4.22) is equivalent to the system of Eqs. (4.13)–(4.16). Similarly, the condition (4.23) gives the system of Eqs. (4.9)–(4.12). Moreover, using Eqs. (4.13)–(4.16) the form (4.22) can be written as

[TABLE]

Hence

[TABLE]

where

[TABLE]

It turns out that the discrete general plane wave solution of Eq. (3.9) consists of four linearly independent solutions for given $p_{\mu}$, $\mu=1,2,3$. To be more precise, there are four linearly independent solutions for each positive and negative $p_{0}=\pm\sqrt{m^{2}+p_{1}^{2}+p_{2}^{2}+p_{3}^{2}}$. Thus the result is the same as for the continuum counterpart. The presence of eight such solutions in the continuum case is interpreted in details in [10].

Similarly, again, if we take the condition (4.23) and use Eqs. (4.9)–(4.12) then we obtain

[TABLE]

This leads to the general solution

[TABLE]

where

[TABLE]

It is clear that the equivalence of (4.22) and (4.23) implies the equivalence of the solutions (4) and (4).