\section{Introduction}

The foundation of the Standard Model of particle physics (SM) is closely connected with the Dirac field and its symmetries. The experimental verification of the SM and the elegance of the Dirac equation, have implicitly led us to assume that all massive spin-half fields must satisfy the Dirac equation and are of mass-dimension three-half. However, this assumption has now been called into question due to the works of Ahluwalia and Grumiller~\cite{Ahluwalia:2004ab,Ahluwalia:2004sz} and it has since been further developed in refs.~\cite{Ahluwalia:2008xi,Ahluwalia:2009rh}. These works proved by construction, the existence of self-interacting massive spin-half fields of mass-dimension one and that they satisfy the Klein-Gordon but not the Dirac equation. In this paper, we generalise these results by constructing self-interacting mass-dimension one fields for any spin.

The  massive spin-half fields of mass-dimension one are constructed using spin-half Elko as expansion coefficients. Elko is a German acronym for \textit{Eigenspinoren des Ladungskonjugationsoperators}. In English, it means eigenspinor of the charge conjugation operator. The precise mathematical definition of Elko for any spin and the proof of its existence are given in the next section.

Elko is closely related to the Majorana spinors. The term Elko was only later introduced to avoid possible confusions. We will not elaborate on the similarities and differences between the two as a detailed exposition can be found in ref.~\cite{Ahluwalia:2009rh}. Prior to the introduction of Elko, the properties of the higher-spin generalisation of Majorana spinors have been investigated by Ahluwalia~\cite{Ahluwalia:1994uy}. While many results obtained by Ahluwalia are in agreement with the results presented in this paper, due to additional insights, the kinematics and the physical interpretation of Elko and the associated massive spin-half fields now differ from the original ones given in ref.~\cite{Ahluwalia:1994uy}.

Although the spin-half Elko have the correct transformation under the $(\frac{1}{2},0)\oplus(0,\frac{1}{2})$ representation of the Lorentz group, their quantum fields are not Lorentz-covariant~\cite{Ahluwalia:2008xi,Ahluwalia:2009rh}.
This result was later generalised to Elko for any spin~\cite{Lee:2012thesis}\footnote{We are using the term \textit{spin} loosely. By a spin-$j$ field, it is meant that the spinors transform under the finite dimensional $(j,0)\oplus(0,j)$ representation.}.  This immediately raises the question on whether the quantum fields under consideration are well-defined. A partial answer was provided by Ahluwalia and Horvath~\cite{Ahluwalia:2010zn}. Their results suggest that the quantum fields satisfy the symmetry of very special relativity (VSR) proposed by Cohen and Glashow~\cite{Cohen:2006ky}.

Here we will not pursue this direction any further as our primary goal is to construct and study the general properties of Elko and their quantum fields for any spin. But as we shall see, the results obtained in this paper do suggest that our constructs are well-defined in the sense that the fields are local and have positive-definite free Hamiltonians.

We restrict ourselves to Lorentz instead of VSR transformations. This is because the properties of spin-half Elko remain unchanged under Lorentz and VSR transformations and we expect this to hold for any spin. Additionally, the general structure of Lorentz transformations is better understood and hence easier to manipulate when performing calculations.

The Lorentz violation for the massive spin-half fields of mass-dimension one can be best characterised by a non Lorentz-covariant term that appear in the Elko spin-sums. Ahluwalia~\cite{Ahluwalia:2013uxa} and Speran\c{c}a~\cite{Speranca:2013hqa} have shown that this term provides an interesting decomposition of $\gamma^{\mu}p_{\mu}$ into a product of two non Lorentz-covariant matrices such that all three are mutually commutative. In this paper, we show that their results can be generalised to any spin thus giving us a decomposition of $\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}$.

The spin-half Elko and their fields of mass-dimension one have attracted interests from various areas of research. In cosmology, it was shown by various authors that Elko has the properties to generate inflation~\cite{Boehmer:2010ma,Boehmer:2010tv,Boehmer:2009aw,Boehmer:2008ah,Boehmer:2008rz,Boehmer:2007ut,
Boehmer:2007dh,Boehmer:2006qq,Chee:2010ju,Wei:2010ad,Shankaranarayanan:2010st,Shankaranarayanan:2009sz,Gredat:2008qf}. The mathematical properties of Elko have been studied in detail by da Rocha et al.~\cite{HoffdaSilva:2009is,
daRocha:2008we,daRocha:2007pz,daRocha:2005ti,daRocha:2011yr,Bernardini:2012sc}.
Wunderle and Dick have used Elko to construct
supersymmetric Lagrangians for spin-half fields of mass dimension one~\cite{Wunderle:2010yw}. Fabbri has shown that the massive spin-half field of mass-dimension one does not violate causality ~\cite{Fabbri:2010va,Fabbri:2009aj,Fabbri:2009ka} while Ahluwalia, Lee and Schritt studied its symmetry~\cite{Ahluwalia:2008xi,Ahluwalia:2009rh}.

The paper is organised as follow. In sec.~\ref{Elko_any_spin}, we construct Elko for any spin. This includes all the essential properties of Elko that are needed to construct a quantum field. This task is carried out in sec.~\ref{field_operators}. From the derived propagators and Lagrangian densities, we conclude that the fields are of mass-dimension one independent of their spin. They are local, but the definitions of the conjugate momenta must be modified for the free Hamiltonians to be physical. The resulting free bosonic and fermionic Hamiltonians are positive-definite but now have positive and vanishing vacuum energies respectively.
Finally, gathering all the relevant results, we show that $\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}$ can be decomposed into a product of two non Lorentz-covariant matrices. This decomposition generalises the result obtained by Ahluwalia~\cite{Ahluwalia:2013uxa} and Speran\c{c}a~\cite{Speranca:2013hqa} and points toward a possible non-trivial connection between mass-dimension one fields and Lorentz-invariant fields.

\section{Elko for any spin}\label{Elko_any_spin}
We follow the procedure outlined in refs.~\cite{Ahluwalia:2004ab,Ahluwalia:2004sz,Ahluwalia:2008xi,Ahluwalia:2009rh} to construct Elko for any spin. The first task is to define a suitable charge conjugation operator $\mathcal{C}$ in the $(j,0)\oplus(0,j)$ representation of the Lorentz group. Here, we define this operator as~\footnote{See app.~\ref{A} for details on the derivation of $\mathcal{C}$.}
\begin{equation}
\mathcal{C}\equiv\left(\begin{matrix}
O & -i\Theta^{-1} \\
-i\Theta & O \end{matrix}\right)\mathcal{K}
\end{equation}
where $\mathcal{K}$ is the complex conjugation operator acting to the right. The matrix $\Theta$ is the Wigner time-reversal operator of dimension $(2j+1)\times(2j+1)$ defined by
\begin{equation}
\Theta\mbox{\boldmath$\displaystyle\mathbf{J}$}\Theta^{-1}=-\mbox{\boldmath$\displaystyle\mathbf{J}$}^{*}\label{eq:Theta}
\end{equation}
where $\mbox{\boldmath$\displaystyle\mathbf{J}$}=(J_{1},J_{2},J_{3})$ are the rotation generators of the $(j,0)$ and $(0,j)$ representations. We adopt the basis where $J_{3}$ is diagonal so that $\mbox{\boldmath$\displaystyle\mathbf{J}$}$ are given by~\cite[eqs.~(2.5.21-2.5.22)]{Weinberg:1995mt}
\begin{equation}
(J_{1}\pm iJ_{2})_{\sigma\bar{\sigma}}
=\delta_{\sigma,\bar{\sigma}\pm1}\sqrt{(j\mp\bar{\sigma})(j\pm\bar{\sigma}+1)},\nonumber
\end{equation}
\begin{equation}
(J_{3})_{\sigma\bar{\sigma}}=\sigma\delta_{\sigma\bar{\sigma}}\label{eq:r3}
\end{equation}
where $\sigma=-j\cdots j$. To prove the existence of $\Theta$, we solve for each of the component of $\mbox{\boldmath$\displaystyle\mathbf{J}$}$ and apply the complex conjugation operation. This gives us the following identity
\begin{equation}
\mbox{\boldmath$\displaystyle\mathbf{J}$}^{*}_{\bar{\sigma}\sigma}=-(-1)^{\bar{\sigma}-\sigma}\mbox{\boldmath$\displaystyle\mathbf{J}$}_{-\bar{\sigma},-\sigma}.
\end{equation}
Rewriting this equation in matrix form, we see that it is equivalent to eq.~(\ref{eq:Theta}). The general solution for $\Theta$ (up to a sign) is then\footnote{We define $\Theta$ such that when $j=\frac{1}{2}$, $\Theta=-i\sigma_{2}$ where $\sigma_{2}$ is the second component of the Pauli matrix.}
\begin{equation}
\Theta_{\bar{\sigma}\sigma}=(-1)^{-j-\bar{\sigma}}\delta_{-\bar{\sigma},\sigma}\label{eq:Theta_soln}
\end{equation}
which guarantees the existence of the Wigner time-reversal operator for any spin.

We are now ready to construct Elko for any spin. We start the construction by adopting the standard polarisation basis, introducing a set of functions $\phi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)$ of the $(0,j)$ representation as eigenfunctions of $J_{3}$~\footnote{In order to avoid confusion, we refer to objects that transform under the $(0,j)$ or $(j,0)$ representation simply as `functions'. The term spinors are strictly reserved for objects that transform under the $(\frac{1}{2},0)$ or $(0,\frac{1}{2})$ representation. Consequently, the German acronym for higher-spin Elko should become
\textit{Eigenfunktionen des Ladungskonjugationsoperators}.}
\begin{equation}
J_{3}\phi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=\sigma\phi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\label{eq:phi0}
\end{equation}
where $\sigma=-j,\cdots, j$ spanning the eigenvalues of $J_{3}$. In the component form, we take the solutions of $\phi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)$ (up to a constant) to be
\begin{equation}
\phi_{\ell}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=\sqrt{\frac{m}{2}}\delta_{\ell\sigma}
\end{equation}
with $m$ being the mass of the particle. The functions of arbitrary momentum $\phi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)$ are obtained by
\begin{eqnarray}
\phi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)&=&\exp\left[i\mathbf{K}^{(0,j)}\cdot\mbox{\boldmath$\displaystyle\mathbf{\varphi}$}\right]\nonumber\\
&=&\exp(-\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\mbox{\boldmath$\displaystyle\mathbf{\varphi}$})\phi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\label{eq:phi}
\end{eqnarray}
where $\mathbf{K}^{(0,j)}=i\mbox{\boldmath$\displaystyle\mathbf{J}$}$ in the $(0,j)$ representation. The rapidity parameter $\mbox{\boldmath$\displaystyle\mathbf{\varphi}$}$ is defined as
\begin{equation}
\cosh\varphi=\frac{\sqrt{|\mbox{\boldmath$\displaystyle\mathbf{p}$}|^{2}+m^{2}}}{m},\hspace{0.5cm}
\sinh\varphi=\frac{|\mbox{\boldmath$\displaystyle\mathbf{p}$}|}{m},\hspace{0.5cm}
\mbox{\boldmath$\displaystyle\mathbf{\varphi}$}=\varphi\hat{\mbox{\boldmath$\displaystyle\mathbf{p}$}}.
\end{equation}
It follows from the existence of $\Theta$ and eq.~(\ref{eq:Theta}) that the set of complex conjugated functions $[\vartheta\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)]$ where $\vartheta$ is a phase to be determined transform under the $(j,0)$ representation
\begin{eqnarray}
[\vartheta\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)]&=&\exp\left[i\mathbf{K}^{(j,0)}\cdot\mbox{\boldmath$\displaystyle\mathbf{\varphi}$}\right][\vartheta\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)]\nonumber\\
&=&\exp\left[\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\mbox{\boldmath$\displaystyle\mathbf{\varphi}$}\right][\vartheta\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)]\label{eq:Theta_phi}
\end{eqnarray}
where $\mathbf{K}^{(j,0)}=-i\mbox{\boldmath$\displaystyle\mathbf{J}$}$.
Acting $J_{3}$ on $[\vartheta\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)]$, we see that they have opposite spin-projection with respect to $\phi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)$
\begin{equation}
J_{3}[\vartheta\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)]=-\sigma [\vartheta\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)].\label{eq:opposite_projection}
\end{equation}
Combining the functions given by eqs.~(\ref{eq:phi}) and (\ref{eq:Theta_phi}), we obtain a set of $2\times(2j+1)$-component functions of the $(j,0)\oplus(0,j)$ representation of the form
\begin{equation}
\chi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=\left(\begin{matrix}
\vartheta\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\\
\phi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\end{matrix}\right).
\end{equation}
Their transformation under boost is given by
\begin{eqnarray}
\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)&=&\kappa(\mbox{\boldmath$\displaystyle\mathbf{p}$})\chi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\nonumber\\
&=&\left(\begin{matrix}
\exp(\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\mbox{\boldmath$\displaystyle\mathbf{\varphi}$}) & O \\
O & \exp(-\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\mbox{\boldmath$\displaystyle\mathbf{\varphi}$}) \end{matrix}\right)\chi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\label{eq:boost}.
\end{eqnarray}
These functions become eigenfunctions of the charge conjugation operator and hence Elko for any spin when $\vartheta=\pm i$
\begin{equation}
\mathcal{C}\chi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\vert_{\vartheta=\pm i}=\mp(-1)^{2j}\chi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\vert_{\vartheta=\pm i}.\label{eq:cchi}
\end{equation}
This relation holds for arbitrary momentum since $\mathcal{C}$ commutes with $\kappa(\mbox{\boldmath$\displaystyle\mathbf{p}$})$ defined in eq.~(\ref{eq:boost}). We note, when $j$ takes half-integral or integral values, the eigenvalues of  Elko are even or odd with respect to the signs of the phase $\vartheta$.
It is important to note that this is not in contradiction with the result obtained in ref.~\cite{Ahluwalia:1994uy}. The spin-one bosonic charge-conjugation matrix of~\citep[eq.~(8)]{Ahluwalia:1994uy} differs from $\mathcal{C}^{(j=1)}$ up to a multiplication of a constant matrix. Once this differences is accounted for, we see that eq.~(\ref{eq:cchi}) is equivalent to~\cite[eq.~(12)]{Ahluwalia:1994uy}.

In order to construct quantum fields, we follow the general structure that a field operator must have two sets of expansion coefficients associated with its particle and anti-particle. Therefore, we introduce the self-conjugate ($\vartheta=i$)  and anti-self-conjugate ($\vartheta=-i$) functions
\begin{subequations}
\begin{equation}
\xi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=\chi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\vert_{\vartheta=i,\phi(\mathbf{0},\sigma)\rightarrow\rho(\mathbf{0},\sigma)}\label{eq:xi}
\end{equation}
\begin{equation}
\zeta(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=\alpha(\sigma)\chi(\mbox{\boldmath$\displaystyle\mathbf{0}$},-\sigma)\vert_{\vartheta=-i,\phi(\mathbf{0},\sigma)\rightarrow\varrho(\mathbf{0},\sigma)}\label{eq:zeta}
\end{equation}
\end{subequations}
with
\begin{equation}
\rho_{\ell}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=\sqrt{\frac{m}{2}}c\delta_{\ell\sigma},\hspace{0.5cm}
\varrho_{\ell}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=\sqrt{\frac{m}{2}}d\delta_{\ell\sigma}
\end{equation}
where $c$ and $d$ are proportionality constants to be determined and $\alpha(\sigma)$ is a phase. The reason behind the introduction of $\rho(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)$ and $\varrho(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)$ and the factor of $\frac{1}{\sqrt{2}}$ will become apparent in sec.~\ref{field_operators} when we study the locality structure of the fields and compute the free Hamiltonians.

The spinor labels and the phase $\alpha(\sigma)$ are non-trivial. For the spin-half fields, their solutions are carefully chosen by the demand of locality when the particles are taken to be Majorana fermions~\cite{Ahluwalia:2008xi,Ahluwalia:2009rh}. As for quantum fields of any spin, consistency demands the general solutions for Elko to coincide with the spin-half case when $j=\frac{1}{2}$. However, it is beyond the scope of this paper to analyse the general locality structure and determine $\alpha(\sigma)$ when the particles are Majorana. Instead, we focus our attention on the quantum fields with distinct particles and anti-particles where the locality structure are insensitive to the choice of $\alpha(\sigma)$.

It is important to note that the original spin-half Elko were constructed in the helicity basis where the two component Weyl spinors of the $(0,\frac{1}{2})$ representation are eigenspinors of $\mbox{\boldmath$\displaystyle\mathbf{\sigma}$}\cdot\hat{\mbox{\boldmath$\displaystyle\mathbf{p}$}}$. It was later discovered that in the polarisation basis, in addition to the labels and phases, the rest spinors must also contain information of a preferred plane under boost for the fields to be local~\cite[eqs.~(17-18)]{Ahluwalia:2009rh}~\footnote{See app.~\ref{B} for details.}. We will have more to say about this issue in sec.~\ref{spin_sums}. But so far, for the sake of simplifying the calculation, it is not necessary to introduce the required phases in the polarisation basis. Later, we will obtain Elko in helicity basis by performing a unitary transformation.

\subsection{Dual coefficients}

In the $(j,0)\oplus(0,j)$ representation, the dual coefficients of an arbitrary function $\psi(\mbox{\boldmath$\displaystyle\mathbf{p}$})$ is usually defined as
\begin{equation}
\overline{\psi}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\equiv\psi^{\dag}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\Gamma\label{eq:dirac_dual}.
\end{equation}
By the demand of Lorentz-invariant norm for the functions, the matrix $\Gamma$ is constrained by the following conditions
\begin{equation}
[\Gamma,\mbox{\boldmath$\displaystyle\mathbf{J}$}]=\{\Gamma,\mathbf{K}\}=O
\end{equation}
where we have used the fact that $\mbox{\boldmath$\displaystyle\mathbf{J}$}$ and $\mathbf{K}$ are Hermitian and anti-Hermitian respectively. Solving the equations and setting the appropriate proportionality constants to unity, we obtain
\begin{equation}
\Gamma=\left(\begin{matrix}
O & I \\
I & O \end{matrix}\right)\label{eq:gamma}
\end{equation}
where the $I$ represents the $(2j+1)\times(2j+1)$ identity matrix. Under this dual, the inner-product between the functions as well as their quantum fields would be Lorentz-invariant thus allowing us to construct Lorentz-invariant Lagrangian densities.

However, Lorentz invariance of the dual alone does not imply that it is universally applicable. In fact, eq.~(\ref{eq:dirac_dual}) is not the correct dual for Elko. Explicit computations show that the following inner-products identically vanish
\begin{equation}
\overline{\xi}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\xi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=\overline{\zeta}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\zeta(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=0.\label{eq:vanishing_norm}
\end{equation}
For the purpose of constructing Lagrangian densities, apart from the Lorentz invariance of the inner-product, the norms of the coefficients must also be non-vanishing and orthonormal. Equation~(\ref{eq:vanishing_norm}) shows that eq.~(\ref{eq:dirac_dual}) is inadequate for Elko and must therefore be replaced. The solution to this problem was first found by Ahluwalia in ref.\cite{Ahluwalia:1994uy} and further developed by Ahluwalia and Grumiller~\cite{Ahluwalia:2004ab,Ahluwalia:2004sz}. Here, we show that their result can be trivially generalised to higher-spin. Their solution is based on the observation that not all the inner-products under eq.~(\ref{eq:dirac_dual}) are vanishing. The general inner-products are given by
\begin{subequations}
\begin{equation}
\overline{\xi}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\xi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma')=\frac{im(-1)^{-j-\sigma}}{2}\left[c^{*2}-(-1)^{-2\sigma}c^{2}\right]\delta_{\sigma,-\sigma'},\label{eq:xi_dual}
\end{equation}
\begin{equation}
\overline{\zeta}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\zeta(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma')=\frac{-im(-1)^{-j-\sigma}}{2}
\left[d^{*2}-(-1)^{-2\sigma}d^{2}\right]\delta_{\sigma,-\sigma'}. \label{eq:zeta_dual}
\end{equation}
\end{subequations}
In light of this observation, we define the dual coefficients for Elko as
\begin{subequations}
\begin{equation}
\gdual{\xi}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\equiv\beta(\sigma)\xi^{\dag}(\mbox{\boldmath$\displaystyle\mathbf{0}$},-\sigma)\Gamma
\end{equation}
\begin{equation}
\gdual{\zeta}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\equiv\beta(\sigma)\zeta^{\dag}(\mbox{\boldmath$\displaystyle\mathbf{0}$},-\sigma)\Gamma
\end{equation}
\end{subequations}
where $\beta(\sigma)$ is a phase to be determined. From eqs.~(\ref{eq:xi_dual}) and (\ref{eq:zeta_dual}), we find that
\begin{equation}
\overline{\xi}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\xi(\mbox{\boldmath$\displaystyle\mathbf{0}$},-\sigma)=-\overline{\xi}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma\pm 1)\xi(\mbox{\boldmath$\displaystyle\mathbf{0}$},-(\sigma\pm 1))\label{eq:Elko_dual}
\end{equation}
and likewise for the anti-self-conjugate functions. For reasons to become apparent in sec.~\ref{kinematics}, we demand all the self conjugate functions to have the same norm and likewise for their anti-self-conjugate counterpart. This is achieved by demanding $\beta(\sigma)$ to satisfy the recursive relation
\begin{equation}
\beta(\sigma)=-\beta(\sigma\pm1).\label{eq:recursive_beta}
\end{equation}
Taking $\sigma=j$ as a point of reference, we obtain
\begin{equation}
\beta(\sigma)=(-1)^{j-\sigma}\beta(j)\label{eq:recursive_beta2}
\end{equation}
where $\sigma=-j,\cdots,+j$.

Computing the Elko norms using eqs.~(\ref{eq:Elko_dual}) and (\ref{eq:recursive_beta2}), we find
\begin{subequations}
\begin{equation}
\gdual{\xi}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\xi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma')=
\begin{cases}
im\beta(j)(c^{*2}+c^{2})\delta_{\sigma\sigma'}/2 & j=\frac{1}{2},\frac{3}{2},\cdots \\
im\beta(j)(c^{*2}-c^{2})\delta_{\sigma\sigma'}/2 & j=1,2,\cdots \label{eq:xi_norms}
\end{cases}
\end{equation}
\begin{equation}
\gdual{\zeta}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\zeta(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma')=
\begin{cases}
-im\beta(j)(d^{*2}+d^{2})\delta_{\sigma\sigma'}/2 & j=\frac{1}{2},\frac{3}{2},\cdots \\
-im\beta(j)(d^{*2}-d^{2})\delta_{\sigma\sigma'}/2 & j=1,2,\cdots. \label{eq:zeta_norms}
\end{cases}
\end{equation}
\end{subequations}
We want the Elko norms to be real, so without the loss of generality, we restrict the square of the proportionality constants $c^{2}$ and $d^{2}$ to be either purely real or purely imaginary. For the norms to be non-vanishing, the proportionality constants must then satisfy
\begin{equation}
c^{*2}=-(-1)^{2j}c^{2},\hspace{0.5cm}
d^{*2}=-(-1)^{2j}d^{2}.\label{eq:cd_constants}
\end{equation}
Substituting eq.~(\ref{eq:cd_constants}) into eqs.~(\ref{eq:xi_norms}) and (\ref{eq:zeta_norms}), the Elko norms simplify to
\begin{subequations}
\begin{equation}
\gdual{\xi}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\xi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma')=\frac{im\beta(j)c^{*2}}{2}\delta_{\sigma\sigma'},\label{eq:xi_norm_final}
\end{equation}
\begin{equation}
\gdual{\zeta}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\zeta(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma')=\frac{-im\beta(j)d^{*2}}{2}\delta_{\sigma\sigma'}.\label{eq:zeta_norm_final}
\end{equation}
\end{subequations}
The phase $\beta(j)$ and the proportionality constants will be determined in sec.~\ref{field_operators} when we study the properties of the quantum fields.

\subsection{Spin-sums and preferred plane}\label{spin_sums}
In the rest frame, using the functions defined by eqs.(\ref{eq:xi}) and (\ref{eq:zeta}), the spin-sums of interests to us are
\begin{subequations}
\begin{equation}
N(\mbox{\boldmath$\displaystyle\mathbf{0}$})=\sum_{\sigma}\xi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\gdual{\xi}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma),\label{eq:xi_spin_sum}
\end{equation}
\begin{equation}
M(\mbox{\boldmath$\displaystyle\mathbf{0}$})=\sum_{\sigma}\zeta(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\gdual{\zeta}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma).\label{eq:zeta_spin_sum}
\end{equation}
\end{subequations}
These spin-sums are important because they determine the propagator and locality structure of the associated quantum fields. We first compute them in the polarisation and then the helicity basis at rest and at arbitrary momentum respectively.

\subsubsection*{Polarisation basis}
As we have mentioned at the end of sec.~\ref{Elko_any_spin}, the spin-half Elko in the polarisation basis must have non-trivial phases to ensure the quantum fields are local. At this stage, we do not know what these specific phases are for Elko of any spin. But fortunately, this is not necessary. By computing $N(\mbox{\boldmath$\displaystyle\mathbf{0}$})$ and $M(\mbox{\boldmath$\displaystyle\mathbf{0}$})$, we can later obtain the correct functions and spin-sums of arbitrary momentum in the helicity basis by performing a unitary transformation. The required phases for Elko in the polarisation basis can then be obtained as a limit of their counterpart in the helicity basis.

Computing the spin-sums using eqs.~(\ref{eq:xi}) and (\ref{eq:zeta}), we get
\begin{subequations}
\begin{equation}
N(\mbox{\boldmath$\displaystyle\mathbf{0}$})=\frac{1}{2}\left[im\beta(j)c^{*2}\left(\begin{matrix}
I & O \\
O & I \end{matrix}\right)
+m\beta(j)|c|^{2}(-1)^{2j}\left(\begin{matrix}
O & \Theta\\
\Theta & O \end{matrix}\right)\right],\label{eq:N0}
\end{equation}
\begin{equation}
M(\mbox{\boldmath$\displaystyle\mathbf{0}$})=\frac{1}{2}\left[-im\beta(j)d^{*2}\left(\begin{matrix}
I & O \\
O & I \end{matrix}\right)
+m\beta(j)|d|^{2}(-1)^{2j}\left(\begin{matrix}
O & \Theta\\
\Theta & O \end{matrix}\right)\right].
\end{equation}
\end{subequations}
To understand why non-trivial phases are required in the polarisation basis, it is sufficient to study the behaviour of $N(\mbox{\boldmath$\displaystyle\mathbf{p}$})$
\begin{eqnarray}
N(\mbox{\boldmath$\displaystyle\mathbf{p}$})&=&\sum_{\sigma}\xi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\gdual{\xi}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\nonumber\\
&=&\kappa(\mbox{\boldmath$\displaystyle\mathbf{p}$})N(\mbox{\boldmath$\displaystyle\mathbf{0}$})\kappa^{-1}(\mbox{\boldmath$\displaystyle\mathbf{p}$})
\end{eqnarray}
where we have used the general identity $\Gamma\kappa^{\dag}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\Gamma=\kappa^{-1}(\mbox{\boldmath$\displaystyle\mathbf{p}$})$. Expanding $N(\mbox{\boldmath$\displaystyle\mathbf{0}$})$ using eq.(\ref{eq:boost}) and eq.~(\ref{eq:N0}) we obtain
\begin{equation}
N(\mbox{\boldmath$\displaystyle\mathbf{p}$})=\frac{1}{2}\left[im\beta(j)c^{*2}\left(\begin{matrix}
I & O \\
O & I \end{matrix}\right)+m\beta(j)|c|^{2}(-1)^{2j}
\left(\begin{matrix}
O & e^{\mathbf{J}\cdot\pmb{\varphi}}\Theta e^{\mathbf{J}\cdot\pmb{\varphi}}\\
e^{-\mathbf{J}\cdot\pmb{\varphi}}\Theta e^{-\mathbf{J}\cdot\pmb{\varphi}} & O \end{matrix}\right)\right].
\end{equation}
For general momentum $\mbox{\boldmath$\displaystyle\mathbf{p}$}$, the sub-matrices $e^{\pm\mathbf{J}\cdot\pmb{\varphi}}\Theta e^{\pm\mathbf{J}\cdot\pmb{\varphi}}$ are not Lorentz covariant. In fact, due to its non-covariance, the resulting quantum fields would be non-local. These observations are in agreement with the previous results in the literature that the massive spin-half fields violate Lorentz symmetry~\cite{Ahluwalia:2008xi,Ahluwalia:2009rh}. In spite of these difficulties, the work of Ahluwalia and Horvath suggests the fields may satisfy VSR symmetry~\cite{Ahluwalia:2010zn}.
Instead of attempting to determine the underlying symmetry group, we choose to focus on higher-spin fields and their implications on particle physics and the general structure of quantum field theory.

Returning to the problem at hand, we note that due to eq.~(\ref{eq:Theta}), the spin-sum is in fact invariant and covariant if it is boosted along the 1, 3-axis and 2-axis respectively
\begin{equation}
N(p^{1})=N(p^{3})=N(\mbox{\boldmath$\displaystyle\mathbf{0}$}),
\end{equation}
\begin{equation}
N(p^{2})=\frac{1}{2}\left[im\beta(j)c^{*2}\left(\begin{matrix}
I & O \\
O & I \end{matrix}\right)
+m\beta(j)|c|^{2}
\left(\begin{matrix}
O & e^{2J_{2}\varphi}\\
e^{-2J_{2}\varphi} & O \end{matrix}\right)\right]\label{eq:1_3_axis_spin_sum}
\end{equation}
where the sub-matrices $e^{\pm 2J_{2}\varphi}$ are Lorentz-covariant.\footnote{See ref.~\cite[app.~A]{Weinberg:1964cn} for more details.}
According to eq.~(\ref{eq:1_3_axis_spin_sum}), when the spin-projection of the particle coincides with the direction of the boost (along the 3-axis), the spin-sum remains invariant.\footnote{The functions of the $(0,j)$ representation are eigenfunctions of $J_{3}$, see eq.~(\ref{eq:phi0}).} Therefore, in the helicity basis where the spin-projections are aligned along the direction of motion, one would expect the spin-sum to be invariant under a general boost.

\subsubsection*{Helicity basis}

Here we determine the spin-sums in the helicity basis by applying a similarity transformation on their counterpart in the polarisation basis. But to achieve this, we first have to determine the relationships between Elko in both basis.

In the helicity basis, the functions of the $(0,j)$ representation are defined as
\begin{equation}
\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\hat{\mbox{\boldmath$\displaystyle\mathbf{p}$}}\phi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)=\sigma\phi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)\label{eq:Jphat_eigenvectors}
\end{equation}
where
\begin{equation}
\hat{\mbox{\boldmath$\displaystyle\mathbf{p}$}}=(\cos\phi\sin\theta,\sin\phi\sin\theta,\cos\theta)
\end{equation}
and $\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$}=\mbox{\boldmath$\displaystyle\mathbf{p}$}\vert_{|\mathbf{p}|\rightarrow0}$ instead of $\mbox{\boldmath$\displaystyle\mathbf{p}$}=0$
so the spin-projections are aligned along $\hat{\mbox{\boldmath$\displaystyle\mathbf{p}$}}$. Since $\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\hat{\mbox{\boldmath$\displaystyle\mathbf{p}$}}$ is Hermitian and has the same eigenvalue spectrum as $J_{3}$, it follows that there exists a unitary matrix $S(\theta,\phi)$ (whose precise form is not of interest to us) such that
\begin{equation}
SJ_{3}S^{\dag}=\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\hat{\mbox{\boldmath$\displaystyle\mathbf{p}$}}.\label{eq:J_3_diagonalise}
\end{equation}
Therefore, instead of solving eq.~(\ref{eq:Jphat_eigenvectors}) for the eigenfunctions, we obtain $\phi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)$ using eq.~(\ref{eq:J_3_diagonalise}) by performing the following unitary transformation
\begin{equation}
\phi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)=S\phi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma).\label{eq:0_j_spinors_helicity}
\end{equation}
Following the same computation that gave us eq.~(\ref{eq:Theta_phi}), we see that
$\vartheta\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)$
transforms under the $(j,0)$ representation and has opposite helicity with respect to $\phi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)$
\begin{equation}
\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\hat{\mbox{\boldmath$\displaystyle\mathbf{p}$}}\left[\vartheta \Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)\right]=-\sigma\left[\vartheta \Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)\right].\label{eq:j_0_spinors_helicity}
\end{equation}
From eq.~(\ref{eq:J_3_diagonalise}) we also obtain
\begin{equation}
J_{3}\left[S^{-1}\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\right]=-\sigma\left[S^{-1}\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)\right]
\end{equation}
thus allowing us to make the following identification\footnote{Rigorously speaking, this identification holds up to a constant. But for $\chi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)$ to be an eigenspinor of the charge conjugation operator, we must take the constant to be unity.}
\begin{equation}
\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)=S\Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma).
\end{equation}
Combining the functions given by eqs.~(\ref{eq:0_j_spinors_helicity}) and (\ref{eq:j_0_spinors_helicity}), we obtain
\begin{equation}
\chi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)=\left(\begin{matrix}
\vartheta \Theta\phi^{*}(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)\\
\phi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)
\end{matrix}\right).
\end{equation}
Setting $\vartheta=\pm i$ and choosing the relevant phases we obtain Elko in the helicity basis. Since the definitions of Elko were already given by eqs.~(\ref{eq:xi}) and (\ref{eq:zeta}) and remain unchanged, we refrain from repeating ourselves. Elko in the helicity and polarisation basis are therefore related by the following unitary transformation
\begin{equation}
\xi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)=\mathcal{S}\xi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma),\hspace{0.5cm}
\zeta(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)=\mathcal{S}\zeta(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)
\end{equation}
where
\begin{equation}
\mathcal{S}=\left(\begin{matrix}
S & O \\
O & S \end{matrix}\right).
\end{equation}
The self-conjugate spin-sum in the helicity basis is then given by
\begin{eqnarray}
N(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$})&=&\mathcal{S}N(\mbox{\boldmath$\displaystyle\mathbf{0}$})\mathcal{S}^{\dag}\nonumber\\
&=&\frac{1}{2}\left[im\beta(j)c^{*2}\left(\begin{matrix}
I & O \\
O & I \end{matrix}\right)
+m\beta(j)|c|^{2}(-1)^{2j}
\left(\begin{matrix}
O & S\Theta S^{\dag} \\
S\Theta S^{\dag} & O \end{matrix}\right)\right].
\end{eqnarray}

Now we need to boost $N(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$})$ to arbitrary momentum and determine the precise form of the sub-matrix $S\Theta S^{\dag}$. The former task can be accomplished as follow. Noting that the boost along the 3-axis is
\begin{equation}
\kappa(p^{3})=\left(\begin{matrix}
e^{J_{3}\varphi} & O \\
O & e^{-J_{3}\varphi} \end{matrix}\right),
\end{equation}
using eq.~(\ref{eq:J_3_diagonalise}), the general boost $\kappa(\mbox{\boldmath$\displaystyle\mathbf{p}$})$ can then be obtained by a similarity transformation
\begin{equation}
\kappa(\mbox{\boldmath$\displaystyle\mathbf{p}$})=\mathcal{S}\kappa(p^{3})\mathcal{S}^{\dag}
\end{equation}
where $|\mbox{\boldmath$\displaystyle\mathbf{p}$}|=p^{3}$. But since $N(p^{3})=N(\mbox{\boldmath$\displaystyle\mathbf{0}$})$, a straightforward calculation shows that the spin-sum in the helicity basis is boost -invariant
\begin{equation}
N_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$})=N(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$}).
\end{equation}
The subscript $h$ indicates the spin-sum of arbitrary momentum is computed in the helicity basis which differs from $N(\mbox{\boldmath$\displaystyle\mathbf{p}$})$. Whenever confusions may arise, any functions in the helicity basis will be accompanied by the subscript $h$.

While there is a well-defined procedure to compute $S(\theta,\phi)$ by diagonalising $\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\hat{\mbox{\boldmath$\displaystyle\mathbf{p}$}}$,  this is not a practical task for large matrices. Instead, it is easier to determine $S\Theta S^{\dag}$ directly. Using eq.~(\ref{eq:Theta}) where $\Theta J_{3}\Theta^{-1}=-J_{3}$, we obtain the following condition
\begin{equation}
\{A,\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\hat{\mbox{\boldmath$\displaystyle\mathbf{p}$}}\}=O \label{eq:anticommute}
\end{equation}
where $A=S\Theta S^{\dag}$ must anti-commute with $\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\hat{\mbox{\boldmath$\displaystyle\mathbf{p}$}}$. However, eq.~(\ref{eq:anticommute}) can only determine $A$ up to a constant. To determine the constant, we use the $j=\frac{1}{2}$ spin-sum as an initial theoretical data and then assume these properties hold for any spin. When $j=\frac{1}{2}$, the self-conjugate spin-sum, with all the proportionality constants determined is\footnote{We will determine all the proportionality constants of Elko in sec.~\ref{field_operators} such that the resulting self-conjugate spin-sum coincides with eq.~(\ref{eq:spin_half_sums}) when $j=\frac{1}{2}$.}
\begin{equation}
N_{h}^{(j=1/2)}(\mbox{\boldmath$\displaystyle\mathbf{p}$})=\frac{m}{2}\left[I+\mathcal{A}^{(j=1/2)}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\right]\label{eq:spin_half_sums}
\end{equation}
where $\mathcal{A}^{(j=1/2)}(\phi)$ is an off-diagonal matrix of the form
\begin{equation}
\mathcal{A}^{(j=1/2)}(\mbox{\boldmath$\displaystyle\mathbf{p}$})=i\left(\begin{matrix}
0 & 0 & 0 & -e^{-i\phi} \\
0 & 0 & e^{i\phi} & 0 \\
0 & -e^{-i\phi} & 0 & 0\\
e^{i\phi} & 0 & 0 & 0 \end{matrix}\right).\label{eq:g_phi}
\end{equation}
As a result of our assumption, the matrix $A$ must be off-diagonal for any spin. An evidence in support of this result is that $\lim_{\theta,\phi\rightarrow 0}A(\theta,\phi)=\Theta$ is an off-diagonal matrix. Therefore, we take $A$ to be
\begin{equation}
A_{\ell m}=g_{\ell}(\theta,\phi)\delta_{\ell,-m}.\label{eq:component_A}
\end{equation}
Substituting eqs.~(\ref{eq:component_A}) and (\ref{eq:r3}) into eq.~(\ref{eq:anticommute}), we obtain the recursive relation
\begin{equation}
g_{-\sigma}(\theta,\phi)=-g_{-\sigma+1}(\theta,\phi)e^{2i\phi}.
\end{equation}
Equating $N_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$})$ with $N_{h}^{(j=1/2)}(\mbox{\boldmath$\displaystyle\mathbf{p}$})$ of eq.~(\ref{eq:spin_half_sums}), modulo their respective undetermined proportionality constants and the factor of $i$ in $\mathcal{A}^{(j=1/2)}(\mbox{\boldmath$\displaystyle\mathbf{p}$})$, we get $g_{-1/2}=e^{i\phi}$. Extrapolate this to all spin, we obtain
\begin{equation}
g_{-j}=e^{i\phi}.\label{eq:g_-j}
\end{equation}
Solving the recursive relation using eq.~(\ref{eq:g_-j}), the general solution of $g_{\ell}(\phi)$ is given by
\begin{equation}
g_{\ell}=(-1)^{j+\ell}e^{-2i(j+\ell-\frac{1}{2})\phi}
\end{equation}
where $\ell=-j,\cdots j$ so that
\begin{equation}
A_{\ell m}=(-1)^{j+\ell}e^{-2i(j+\ell-\frac{1}{2})\phi}\delta_{\ell,-m}.\label{eq:A_soln}
\end{equation}
For example, when $j=\frac{1}{2}$, we have
\begin{equation}
A^{(j=1/2)}=\left(\begin{matrix}
0 & -e^{-i\phi} \\
e^{i\phi} & 0 \end{matrix}\right)
\end{equation}
and for $j=1$,
\begin{equation}
A^{(j=1)}=\left(\begin{matrix}
0 & 0 & e^{-3i\phi}\\
0 & -e^{-i\phi} & 0 \\
e^{i\phi} & 0 & 0 \end{matrix}\right).
\end{equation}
Therefore, in the helicity basis, the self-conjugate and anti-self-conjugate spin-sums are
\begin{subequations}
\begin{equation}
N_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$})=\sum_{\sigma}\xi_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\gdual{\xi}_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=\frac{1}{2}\left[ im\beta(j)c^{*2}I
+m\beta(j)|c|^{2}(-1)^{2j}
\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\right],
\end{equation}
\begin{equation}
M_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$})=\sum_{\sigma}\zeta_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\gdual{\zeta}_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=\frac{1}{2}\left[-im\beta(j)d^{*2}I
+m\beta(j)|d|^{2}(-1)^{2j}
\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\right]
\end{equation}
\end{subequations}
where
\begin{equation}
\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})=\left(\begin{matrix}
O & A \\
A & O \end{matrix}\right).\label{eq:mathcalA}
\end{equation}
The momentum dependence of $\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})$ is due to the fact that $e^{\pm i\phi}=(p^{1}\pm ip^{2})/(|\mbox{\boldmath$\displaystyle\mathbf{p}$}|\sin\theta)$.

Now  we return to the problem of determining Elko and its spin-sums in the polarisation basis. This can now be solved by noting that the spin-sums in the helicity basis specify a preferred plane with only $\phi$-dependence so the limit $\theta\rightarrow0$ does not affect the spin-sum. In this limit, $\phi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)\vert_{\theta\rightarrow0}$ becomes a spinor in the polarisation basis in the sense that it is now an eigenspinor of $J_{3}$
\begin{equation}
J_{3}\phi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)\vert_{\theta\rightarrow0}=\sigma \phi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)\vert_{\theta\rightarrow0}.
\end{equation}
The remaining $\phi$-dependence for $\phi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)\vert_{\theta\rightarrow0}$ now appears as non-trivial phases. In the spin-half case, these phases ensure the locality of the fields.  Therefore, the self-conjugate and anti-self-conjugate Elko for any spin in the polarisation basis are given by
\begin{subequations}
\begin{equation}
\xi(\mbox{\boldmath$\displaystyle\mathbf{0}$}_{\phi},\sigma)=\xi(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)\vert_{\theta\rightarrow0} ,
\end{equation}
\begin{equation}
\zeta(\mbox{\boldmath$\displaystyle\mathbf{0}$}_{\phi},\sigma)=\zeta(\mbox{\boldmath$\displaystyle\mathbf{\epsilon}$},\sigma)\vert_{\theta\rightarrow0}
\end{equation}
\end{subequations}
where subscripts denote the directional-dependence of the rest functions. Similarly the corresponding Elko of arbitrary momentum is denoted as
\begin{subequations}
\begin{equation}
\xi(\mbox{\boldmath$\displaystyle\mathbf{p}$}_{\phi},\sigma)=\xi_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\vert_{\theta\rightarrow0},
\end{equation}
\begin{equation}
\zeta(\mbox{\boldmath$\displaystyle\mathbf{p}$}_{\phi},\sigma)=\zeta_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\vert_{\theta\rightarrow0}.
\end{equation}
\end{subequations}
In the helicity basis, Elko of arbitrary momentum satisfy
\begin{subequations}
\begin{equation}
\xi_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=\mathcal{S}\kappa(p^{3})\xi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=\mathcal{S}\xi(p^{3},\sigma),\label{eq:helicity_xi_p}
\end{equation}
\begin{equation}
\zeta_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=\mathcal{S}\kappa(p^{3})\zeta(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=\mathcal{S}\xi(p^{3},\sigma).\label{eq:helicity_zeta_p}
\end{equation}
\end{subequations}
On the left-hand side we have Elko of arbitrary momentum in the helicity basis while Elko on the right-hand side are defined in the polarisation basis according to eq.~(\ref{eq:phi0}). From now onwards, we drop the $h$ subscript as the results obtained for the rest of the paper hold in both helicity and polarisation basis.

To conclude this section, the introduction of direction-dependent non-trivial phases ensure the Elko spin-sums remain the same in both helicity and polarisation basis. Our results show that the Elko spin-sums for any spin are not Lorentz covariant and contain a preferred plane. Consequently, the resulting quantum fields violate Lorentz symmetry. Determining the underlying symmetry for such quantum fields is undoubtedly an important task, but it is not one we will undertake in this paper since our starting point is not suitable for this task.
A potentially more appropriate approach is to start with a symmetry group and derive the particle states and quantum fields following the formalism developed by Wigner and Weinberg~\cite{Wigner:1939cj,Weinberg:1964cn,Weinberg:1964ev}.

We should add that on this front, significant progress has been made by Ahluwalia and Horvath~\cite{Ahluwalia:2010zn}. They showed that the Elko spin-sums are covariant under VSR transformations proposed by Cohen and Glashow~\cite{Cohen:2006ky} thus suggesting Elko satisfies VSR symmetry. However, to show that the fields for any spin satisfy VSR symmetry, one must be able to derive Elko from first principle. This task has not yet been accomplished. Therefore, we leave this as an open problem for future investigation. For now, in the context of this paper, the next task is to derive the equations of motion for Elko in the momentum space.

\subsection{Higher-spin field equations}

The field equation can be derived by exploiting the following general identity in the $(j,0)\oplus(0,j)$ representation
\begin{equation}
\kappa(p)\Gamma\kappa^{-1}(\mbox{\boldmath$\displaystyle\mathbf{p}$})=\left(\begin{matrix}
O & \exp(2\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\mbox{\boldmath$\displaystyle\mathbf{\varphi}$}) \\
\exp(-2\mbox{\boldmath$\displaystyle\mathbf{J}$}\cdot\mbox{\boldmath$\displaystyle\mathbf{\varphi}$}) & O \end{matrix}\right)\equiv
\frac{(-1)^{2j}}{m^{2j}}\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}},\label{eq:gamma_2j}
\end{equation}
\begin{equation}
\gamma^{\mu_{1}\cdots\mu_{2j}}\gamma^{\nu_{1}\cdots\nu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}p_{\nu_{1}}\cdots p_{\nu_{2j}}=
(p^{\mu}p_{\mu})^{2j}=m^{4j}.\label{eq:generalised_dispersion}
\end{equation}
where $\gamma^{\mu_{1}\cdots\mu_{2j}}$ is a generalisation of the spin-half $\gamma^{\mu}$ matrices~\cite[app.~A]{Weinberg:1964cn}. These matrices are completely symmetric in all indices and is a rank-$2j$ tensor with respect to the Lorentz transformation of the $(j,0)\oplus(0,j)$ representation.
Here we shall work with $\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)$ leaving the phase $\vartheta$ unfixed for convenience. This is possible because it does not affect the final results. To see how eq.~(\ref{eq:gamma_2j}) helps us, we compute $\Gamma\chi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)$ and find
\begin{equation}
\Gamma\chi(\mbox{\boldmath$\displaystyle\mathbf{0}$},\sigma)=\pm\vartheta f(j)(-1)^{-j+\sigma} \chi(\mbox{\boldmath$\displaystyle\mathbf{0}$},-\sigma)\label{eq:rest_equation_chi}
\end{equation}
where $f(j)$ is defined as
\begin{equation}
f(j)=\begin{cases}
1 & j=\frac{1}{2},\frac{3}{2},\cdots\\
i & j=1,2,\cdots
\end{cases}
\end{equation}
To derive eq.~(\ref{eq:rest_equation_chi}), it is easier to consider the bosonic and fermionic Elko separately. The $\pm$ sign originates from taking the square root of eq.~(\ref{eq:cd_constants})
\begin{equation}
c^{*}=\pm(-1)^{j+1/2}c,\hspace{0.5cm}
d^{\,*}=\pm(-1)^{j+1/2}d.
\end{equation}
For generality, we keep the freedom of the $\pm$ sign but the choice is inconsequential since the field equation is insensitive to them. Multiply eq.~(\ref{eq:rest_equation_chi}) from the left by $\kappa(\mbox{\boldmath$\displaystyle\mathbf{p}$})$ and using eq.~(\ref{eq:gamma_2j}), we obtain
\begin{subequations}
\begin{equation}
\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=\pm f(j)(-1)^{-j+\sigma}\vartheta m^{2j}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},-\sigma) \label{eq:chi_eq1}
\end{equation}
\begin{equation}
\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},-\sigma)=\mp f^{-1}(j)(-1)^{j-\sigma}\vartheta m^{2j}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma).\label{eq:chi_eq2}
\end{equation}
\end{subequations}
The above equations show that Elko for any spin does not satisfy the usual Lorentz-invariant higher-spin equations~\footnote{See app.~\ref{A} for details.}. Instead the operator $\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}$ maps $\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)$ to $\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},-\sigma)$. Substituting eq.~(\ref{eq:chi_eq2}) into eq.~(\ref{eq:chi_eq1}) or vice versa, we obtain
\begin{equation}
\left(\gamma^{\mu_{1}\cdots\mu_{2j}}\gamma^{\nu_{1}\cdots\nu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}p_{\nu_{1}}\cdots p_{\nu_{2j}}-m^{4j}I\right)\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=0.
\end{equation}
From eq.~(\ref{eq:generalised_dispersion}), this is becomes a higher-spin generalisation of the Klein-Gordon equation
\begin{equation}
[(p^{\mu}p_{\mu})^{2j}-m^{4j}]\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=0
\end{equation}
and consequently
\begin{equation}
[(p^{\mu}p_{\mu})^{2j}-m^{4j}]\xi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=[(p^{\mu}p_{\mu})^{2j}-m^{4j}]\zeta(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=0.\label{eq:Elko_KG}
\end{equation}

While eq.~(\ref{eq:Elko_KG}) is trivial, it is important to note that unlike the expansion coefficients for a massive Lorentz-invariant quantum field of the $(j,0)\oplus(0,j)$ representation, Elko does not satisfy the generalised higher-spin Dirac equation. Instead it satisfies the higher-spin generalisation of the Klein-Gordon equation for any spin. This poses an interesting problem, how do we write down a Lagrangian for the corresponding field when $j>\frac{1}{2}$? We answer this question in sec.~\ref{kinematics} by computing the propagator of the spin-$j$ field.
As a result, it is shown that the fields for any spin have mass-dimension one and renormalisable self-interactions.
\section{Mass-dimension one fields} \label{field_operators}

In this section, we explicitly construct quantum fields using Elko in the helicity basis as expansion coefficients. The unusual properties of Elko that we have encountered are naturally inherited by the field operators and will be revealed once we derive the free propagator and Lagrangian density. During the process, without losing generality, we will also determine the remaining proportionality constants of Elko.

We take the Elko field with the appropriate normalisation to be
\begin{equation}
\Lambda(x)=(2\pi)^{-3/2}\int\frac{d^{3}p}{\sqrt{2mp^{0}}}\sum_{\sigma}
\left[e^{-ip\cdot x}\xi_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)a(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)+e^{ip\cdot x}\zeta_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)b^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\right]
\end{equation}
where the creation and annihilation operators satisfy the canonical algebra
\begin{equation}
[a(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma),a^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$}',\sigma')]_{\pm}=[b(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma),b^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$}',\sigma')]_{\pm}=\delta_{\sigma\sigma'}\delta^{3}(\mbox{\boldmath$\displaystyle\mathbf{p}$}-\mbox{\boldmath$\displaystyle\mathbf{p}$}')
\label{eq:creation_annihilation}
\end{equation}
while all other combinations vanish. The particles and anti-particles are created by acting $a^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)$ and $b^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)$ on the vacuum $|\,\,\rangle$. The operator $\ddag$ has the same operational meaning as Hermitian conjugation when applied to the particle state but its action on the coefficients is different. We define its operation on Elko as
\begin{subequations}
\begin{equation}
\xi^{\ddag}_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\equiv\gdual{\xi}_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\Gamma,
\end{equation}
\begin{equation}
\zeta^{\ddag}_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\equiv\gdual{\zeta}_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\Gamma.
\end{equation}
\end{subequations}
Therefore, the adjoint of the quantum field is given by
\begin{equation}
\dual{\Lambda}(x)\equiv(2\pi)^{-3/2}\int\frac{d^{3}p}{\sqrt{2mp^{0}}}\sum_{\sigma}
\left[e^{ip\cdot x}\gdual{\xi}_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)a^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)+e^{-ip\cdot x}\gdual{\zeta}_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)b(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\right].
\label{eq:L_Lad}
\end{equation}

The work of Ahluwalia and Horvath~\cite{Ahluwalia:2010zn} suggested that Elko satisfies VSR symmetry proposed by Cohen and Glashow~\cite{Cohen:2006ky}. Since the VSR groups are subgroups of the Poincar\'{e} group, we expect $\Lambda(x)$ and $\dual{\Lambda}(x)$ to preserve certain properties of Lorentz invariant field theories. In this spirit, we compute $[\Lambda(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\dual{\Lambda}(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{\pm}$ to examine the spin-statistics and locality structure of the quantum field. Using eq.~(\ref{eq:creation_annihilation}), we get
\begin{equation}
[\Lambda(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\dual{\Lambda}(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{\pm}=(2\pi)^{-3}
\int\frac{d^{3}p}{2mp^{0}}e^{i\mathbf{p\cdot(x-y)}}\left[N_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\pm M_{h}(-\mbox{\boldmath$\displaystyle\mathbf{p}$})\right].
\end{equation}
The term in the square bracket reads
\begin{equation}
N_{h}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\pm M_{h}(-\mbox{\boldmath$\displaystyle\mathbf{p}$})=\frac{1}{2}\left[im\beta(j)\left(c^{*2}\mp d^{*2}\right)I+m\beta(j)(-1)^{2j}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\left(|c|^{2}\mp |d|^{2}\right)\right]
\label{eq:N_pm_M}
\end{equation}
where we have used the fact that $\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})=-\mathcal{A}(-\mbox{\boldmath$\displaystyle\mathbf{p}$})$ and
\begin{equation}
-\mbox{\boldmath$\displaystyle\mathbf{p}$}=\mbox{\boldmath$\displaystyle\mathbf{p}$}\vert_{\phi\rightarrow\phi\pm\pi,\theta\rightarrow\pi-\theta}.
\end{equation}
The plus and minus sign apply when the momentum along the 2-axis is positive or negative respectively. Substituting eq.~(\ref{eq:N_pm_M}) into eq.~(\ref{eq:L_Lad}), we get
\begin{equation}
[\Lambda(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\dual{\Lambda}(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{\pm}=(2\pi)^{-3}\int\frac{d^{3}p}{4p^{0}}e^{i\mathbf{p\cdot(x-y)}}
\left[i\beta(j)\left(c^{*2}\mp d^{*2}\right)I+\beta(j)(-1)^{2j}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\left(|c|^{2}\mp |d|^{2}\right)\right].
\end{equation}
For the fermionic field, the anti-commutator identically vanishes
\begin{equation}
[\Lambda(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\dual{\Lambda}(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{+}=0,\hspace{0.5cm} j=\frac{1}{2},\frac{3}{2},\cdots
\end{equation}
provided that we choose
\begin{equation}
c^{*2}=d^{*2},\hspace{0.5cm} j=\frac{1}{2},\frac{3}{2},\cdots.
\end{equation}
On the other hand, for the bosonic field, if $c^{*2}=d^{*2}$ then the first term of the commutator would be proportional to the standard function
\begin{equation}
\Delta_{+}(\mathbf{x-y})=(2\pi)^{-3}\int\frac{d^{3}p}{2p^{0}}e^{i\mathbf{p\cdot(x-y)}}
\end{equation}
which is non-vanishing for all $\mathbf{x-y}$. Therefore, we must choose
\begin{equation}
c^{*2}=-d^{*2},\hspace{0.5cm} j=1,2,\cdots.
\end{equation}
The resulting bosonic commutator becomes
\begin{equation}
[\Lambda(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\dual{\Lambda}(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{-}=|c|^{2}\beta(j)(2\pi)^{-3}\int\frac{d^{3}p}{2p^{0}}e^{i\mathbf{p\cdot(x-y)}}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$}),\hspace{0.5cm}j=1,2,\cdots\label{eq:A_integral}
\end{equation}
Generally, such an integral
\begin{equation}
\int\frac{d^{3}p}{2p^{0}}e^{i\mathbf{p\cdot(x-y)}}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$}),\label{eq:Fourier_transform_A}
\end{equation}
for any spin, is non-vanishing unless $\mathbf{x-y}$ is aligned to the 3-axis. Therefore, eq.~(\ref{eq:A_integral}) defines a preferred axis of locality for the bosonic fields~\footnote{I am grateful to D.~V. Ahluwalia for discussions on this point. Equation~(\ref{eq:Fourier_transform_A}) is formally a Fourier transform involving $\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})$. For the rest of the paper, we will encounter various Fourier transforms of $\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})$. They are all non-vanishing unless $\mathbf{x-y}$ is aligned to the preferred 3-axis.}. We leave the evaluation of the integral for future investigation.

The existence of the preferred axis is not new, it was first identified for spin-half fermions in refs.~\cite{Ahluwalia:2008xi,Ahluwalia:2009rh}. Here we show that the preferred axis features prominently for Elko and its quantum fields of any spin. They appear in the free propagators and the equal-time locality anti-commutators and commutators between the fields and their conjugate momenta.

\subsection{Free propagators}\label{kinematics}

We are now ready to derive the free propagator. Using the quantum field and its adjoint, the propagator may be obtained by computing $\langle\,\,|T[\Lambda(x),\gdual{\Lambda}(y)]|\,\,\rangle$. Substituting the Elko spin-sums into the generic expression of the propagator given by eq.~(\ref{eq:g_prop}), we obtain the bosonic and fermionic propagator
\begin{subequations}
\begin{equation}
S(y,x)=\frac{i}{2}\int\frac{d^{4}q}{(2\pi)^{4}}e^{-iq\cdot(x-y)}\left[\frac{i\beta(j)c^{*2}I+(q^{0}/\sqrt{|\mathbf{q}|^{2}+m^{2}})\beta(j)|c|^{2}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{q}$})}{q^{\mu}q_{\mu}-m^{2}+i\epsilon}\right],\,j=1,2,\cdots \label{eq:bosonic_propagator}
\end{equation}
\begin{equation}
S(y,x)=
\frac{i}{2}\int\frac{d^{4}q}{(2\pi)^{4}}e^{-iq\cdot(x-y)}\left[\frac{i\beta(j)c^{*2}I-\beta(j)|c|^{2}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{q}$})}{q^{\mu}q_{\mu}-m^{2}+i\epsilon}\right], \, j=\frac{1}{2},\frac{3}{2},\cdots \label{eq:fermionic_propagator}
\end{equation}
\end{subequations}
Equations~(\ref{eq:bosonic_propagator}) and (\ref{eq:fermionic_propagator}) are the main important result of this paper. The integrals over $\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{q}$})$ for the bosonic and fermionic propagators are generally non-zero unless $\mathbf{x-y}$ is aligned along the preferred 3-axis. A simple dimensional analysis shows that the fields for any spin all have mass-dimension one. This result is in stark contrast with the well-known result that the massive Lorentz-invariant fields of the $(j,0)\oplus(0,j)$ representation are of mass-dimension $j+1$. But in hindsight, since Lorentz symmetry is violated, the general properties of Lorentz-invariant field theories may no longer be applicable.

Because the integral over $\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{q}$})$ is non-zero, the propagators are no longer the Green's function of the Klein-Gordon operator $\partial^{\mu}\partial_{\mu}+m^{2}$ where $\partial_{\mu}=\partial/\partial x^{\mu}$.
In spite of this, given that Elko satisfies the higher-spin generalisation of the Klein-Gordon equation and that this equation is satisfied due to the dispersion relation $p^{\mu}p_{\mu}=m^{2}$, taking into account of the term proportional to $\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{q}$})$ that appears in the propagator, we propose that the propagator is a Green's function to the following operator of the form
\begin{equation}
(\partial^{\mu}\partial_{\mu}+m^{2}+\mathcal{F})S(y,x)=i\delta^{4}(x-y)
\end{equation}
where $\mathcal{F}$ is a non-local term~\footnote{By non-local, we mean that the term has more than one space-time coordinate dependence. A term is local if it only depends on only one space-time coordinate}. The Lagrangian density for $\Lambda(x)$ and $\dual{\Lambda}(x)$ for any spin with renormalisable self-interaction is then
\begin{eqnarray}
\mathscr{L}&=&\mathscr{L}_{0}-\frac{g}{4}(\gdual{\Lambda}\Lambda)^{2}+\mathcal{F}'(\Lambda,\dual{\Lambda})\nonumber\\
&=&\partial^{\mu}\gdual{\Lambda}\partial_{\mu}\Lambda-m^{2}\gdual{\Lambda}\Lambda-\frac{g}{4}(\gdual{\Lambda}\Lambda)^{2}+\mathcal{F}'(\Lambda,\dual{\Lambda})
\label{eq:Elko_Lagrangian}
\end{eqnarray}
where $g$ is the coupling constant. We leave the task of determining the non-local operators $\mathcal{F}$ and $\mathcal{F}'$ for future investigation.

At this stage, the physical interpretations for $\mathcal{F}$ and $\mathcal{F}'$ are unclear, but we expect them to affect the full Hamiltonian. As for the free Hamiltonian and conjugate momentum, they will only depend on $\mathscr{L}_{0}$. In order to convince the reader that $\mathscr{L}_{0}$ is the correct free Lagrangian density, in the next two sections, we derive the spin-statistics relations of the fields and compute the free Hamiltonians.

\subsection{Canonical variables}\label{spin_statistics}

The particles and anti-particles we have introduced are by construction, physically distinct. As a result, the corresponding Lagrangian density has two independent variables $\Lambda(x)$ and $\dual{\Lambda}(x)$ and are accompanied by their respective conjugate momenta $\Pi(x)$ and $\dual{\Pi}(x)$.

We define the conjugate momenta to be
\begin{equation}
\Pi(x)\equiv\frac{\partial\mathscr{L}}{\partial(\partial\Lambda/\partial t)},\hspace{0.5cm}
\dual{\Pi}(x)\equiv(-1)^{2j}\frac{\partial\mathscr{L}}{\partial(\partial\gdual{\Lambda}/\partial t)}.
\label{eq:conjugate_momenta}
\end{equation}
In sec.~\ref{free_hamiltonian}, we show that the factor $(-1)^{2j}$ is required in order to obtain physical free Hamiltonians. Here this is necessary if we want the canonical variables to satisfy
\begin{equation}
[Q_{n}(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),P_{m}(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{\pm}\sim i\delta_{nm}\delta^{3}(\mbox{\boldmath$\displaystyle\mathbf{x}$}-\mbox{\boldmath$\displaystyle\mathbf{y}$}),\hspace{0.5cm}n,m=1,2\label{eq:QP_algebra}
\end{equation}
up to a real positive constant with
\begin{equation}
Q_{1}(x)=\Lambda(x),\hspace{0.5cm}
Q_{2}(x)=\dual{\Lambda}(x)
\end{equation}
\begin{equation}
P_{1}(x)=\Pi(x),\hspace{0.5cm}
P_{2}(x)=\dual{\Pi}(x).
\end{equation}
A straightforward computation gives us the following equal-time bosonic commutator
\begin{equation}
[\Lambda(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\Pi(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{-}=[\dual{\Lambda}(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\dual{\Pi}(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{-}=
-\beta(j)c^{*2}\delta^{3}(\mbox{\boldmath$\displaystyle\mathbf{x}$}-\mbox{\boldmath$\displaystyle\mathbf{y}$}),\hspace{0.5cm} j=1,2,\cdots
\end{equation}
and fermionic anti-commutator
\begin{subequations}
\begin{eqnarray}
\hspace{-1cm}[\Lambda(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\Pi(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{+}&=&-\frac{\beta(j)c^{*2}}{2}\delta^{3}(\mbox{\boldmath$\displaystyle\mathbf{x}$}-\mbox{\boldmath$\displaystyle\mathbf{y}$})I  \nonumber\\
&&-\frac{i\beta(j)|c|^{2}}{2(2\pi)^{3}}\int d^{3}p\, e^{-i\mathbf{p\cdot(x-y)}}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$}),\hspace{0.5cm}j=\frac{1}{2},\frac{3}{2},\cdots\label{eq:anti_commutator1}
\end{eqnarray}
\begin{eqnarray}
\hspace{-1cm}[\dual{\Lambda}(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\dual{\Pi}(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{+}&=&\frac{-\beta(j)c^{*2}}{2}\delta^{3}(\mbox{\boldmath$\displaystyle\mathbf{x}$}-\mbox{\boldmath$\displaystyle\mathbf{y}$})I  \nonumber\\
&&-\frac{i\beta(j)|c|^{2}}{2(2\pi)^{3}}\int d^{3}p\,e^{i\mathbf{p\cdot(x-y)}}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$}),\hspace{0.5cm}j=\frac{1}{2},\frac{3}{2},\cdots.\label{eq:anti_commutator2}
\end{eqnarray}
\end{subequations}
We see that the bosonic commutator is consistent with our expectation while the fermionic anti-commutator contains an additional integral over $\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})$. This integral is formally a Fourier transform of $\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{q}$})$ so it is generally non-vanishing unless $\mathbf{x-y}$ is aligned to the 3-axis. For $j=\frac{1}{2}$, the integral has been explicitly evaluated~\footnote{D.~V.~Ahluwalia, private communication.}.

Comparing the above expressions with eq.(\ref{eq:QP_algebra}), we see that the coefficient $-\beta(j)c^{*2}$ of the Dirac delta function must be an imaginary number. Had we not introduced the factor $(-1)^{2j}$ for $\dual{\Pi}(x)$, the first term of eqs.~(\ref{eq:anti_commutator1}) and (\ref{eq:anti_commutator2}) would differ by a minus sign and becomes inconsistent with eq.~(\ref{eq:QP_algebra}).

\subsection{Free Hamiltonian}\label{free_hamiltonian}
The signs of the proportionality constant in eq.~(\ref{eq:QP_algebra}) may not be particularly important. But the overall sign of a Hamiltonian is important since it determines whether the theory is physical or not. Below we show that using the conjugate momenta defined in eq.~(\ref{eq:conjugate_momenta}), we obtain the desired positive-definite free Hamiltonians for any spin.

We compute the free Hamiltonian using the Legendre transformation
\begin{eqnarray}
H_{0}&=&\int d^{3}x \sum_{n}P_{n}(x)\dot{Q}_{n}(x)-\mathscr{L}_{0} \nonumber \\
&=&\int d^{3}x \left[
(-1)^{2j}\frac{\partial\Lambda}{\partial t}\frac{\partial\dual{\Lambda}}{\partial t} -\partial_{i}\dual{\Lambda}\partial^{i}\Lambda+m^{2}\dual{\Lambda}\Lambda\right]
\end{eqnarray}
Substituting $\Lambda(x)$ and $\dual{\Lambda}(x)$ into $H_{0}$, we obtain
\begin{eqnarray}
H_{0}&=&\int\frac{d^{3}p}{m}p^{0}\sum_{\sigma}\Big[\dual{\xi}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\xi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)[a^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)a(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)+\delta(\mbox{\boldmath$\displaystyle\mathbf{0}$})]\nonumber\\
&&\hspace{3cm}+\dual{\zeta}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\zeta(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)[(-1)^{2j}b^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)b(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)+\delta(\mbox{\boldmath$\displaystyle\mathbf{0}$})]\Big].
\end{eqnarray}
Although we do not show the explicit calculation here, it is important to remark that if the $(-1)^{2j}$ factor is absent from $\dual{\Pi}(x)$, the cross-terms proportional to $a(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)b(-\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma')$ and $a^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)b^{\ddag}(-\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma')$ would be non-vanishing. Therefore, the introduction of $(-1)^{2j}$ to $\dual{\Pi}(x)$ is crucial in ensuring that the free fermionic Hamiltonian is diagonal. Also, we see that the norms of the self-conjugate and anti-self- conjugate Elko must be identical within their respective sectors so they contribute equally to the free Hamiltonian.

Using eq.~(\ref{eq:cd_constants}) and the Elko norms given in eqs.~(\ref{eq:xi_norm_final}) and (\ref{eq:zeta_norm_final}), we see that for the fermionic Hamiltonian, the divergent terms proportional to $\delta^{3}(\mbox{\boldmath$\displaystyle\mathbf{0}$})$ cancel giving us
\begin{equation}
H_{0}=2i\beta(j)c^{*2}\int d^{3}p\sum_{\sigma}p^{0}\left[a^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)a(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)
+b^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)b(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\right],\hspace{0.5cm}j=\frac{1}{2},\frac{3}{2},\cdots
\end{equation}
Therefore, the fermionic field for Elko has vanishing vacuum energy. Setting the proportionality constants of the fermionic field and its adjoint to be
\begin{equation}
\beta(j)=-i,\hspace{0.5cm} c=d=1,\hspace{0.5cm}j=\frac{1}{2},\frac{3}{2},\cdots
\end{equation}
we obtain the desired free fermionic Hamiltonian with no vacuum energy
\begin{equation}
H_{0}=\int d^{3}p\sum_{\sigma}p^{0}\left[a^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)a(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)
+b^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)b(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\right],\hspace{0.5cm}j=\frac{1}{2},\frac{3}{2},\cdots.
\end{equation}
Performing a similar computation for the bosonic field, in order to obtain a physical Hamiltonian, we set the proportionality constants to
\begin{equation}
\beta(j)=1,\hspace{0.5cm} c=d^{*}=e^{i\pi/4},\hspace{0.5cm} j=1,2,\cdots
\end{equation}
which gives us
\begin{equation}
H_{0}=\int d^{3}p\sum_{\sigma}p^{0}\left[a^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)a(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)
+b^{\ddag}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)b(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)+2\delta^{3}(\mbox{\boldmath$\displaystyle\mathbf{0}$})\right],\hspace{0.5cm}j=1,2,\cdots.
\end{equation}
Unlike the fermions, the vacuum energy is positive and non-vanishing, but since only the shift in vacuum energy is measurable, the divergent term can be ignored. It is helpful to remind ourselves that it is unusual for the fermionic field to have vanishing vacuum energy. In Lorentz invariant field theories, the free Hamiltonians of bosons and fermions have positive and negative divergent vacuum energies respectively.

We have now determine all the proportionality constants. Substituting them into the canonical commutators and anti-commutators we obtain
\begin{equation}
[\Lambda(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\Pi(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{-}=[\dual{\Lambda}(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\dual{\Pi}(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{-}=
\frac{i}{2}\delta^{3}(\mbox{\boldmath$\displaystyle\mathbf{x}$}-\mbox{\boldmath$\displaystyle\mathbf{y}$})I,\hspace{0.5cm} j=1,2,\cdots
\end{equation}
\begin{subequations}
\begin{eqnarray}
\hspace{-1cm}[\Lambda(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\Pi(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{+}&=&\frac{i}{2}\delta^{3}(\mbox{\boldmath$\displaystyle\mathbf{x}$}-\mbox{\boldmath$\displaystyle\mathbf{y}$})I  \nonumber\\
&&-\frac{1}{2(2\pi)^{3}}\int d^{3}p\, e^{-i\mathbf{p\cdot(x-y)}}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$}),\hspace{0.5cm}j=\frac{1}{2},\frac{3}{2},\cdots
\end{eqnarray}
\begin{eqnarray}
\hspace{-1cm}[\dual{\Lambda}(\mbox{\boldmath$\displaystyle\mathbf{x}$},t),\dual{\Pi}(\mbox{\boldmath$\displaystyle\mathbf{y}$},t)]_{+}&=&\frac{i}{2}\delta^{3}(\mbox{\boldmath$\displaystyle\mathbf{x}$}-\mbox{\boldmath$\displaystyle\mathbf{y}$})I  \nonumber\\
&&-\frac{1}{2(2\pi)^{3}}\int d^{3}p\,e^{i\mathbf{p\cdot(x-y)}}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$}),\hspace{0.5cm}j=\frac{1}{2},\frac{3}{2},\cdots.
\end{eqnarray}
\end{subequations}
We see that the demand for physical Hamiltonians forces the canonical algebra to be modified by a factor of 1/2. This factor cannot be absorbed into the fields for it would then change the free Hamiltonians by a factor of 2.  Substituting these constants into the propagators, we obtain
\begin{subequations}
\begin{equation}
S(y,x)=\frac{i}{2}\int\frac{d^{4}q}{(2\pi)^{4}}e^{-iq\cdot(x-y)}\left[\frac{I-i(q^{0}/\sqrt{|\mathbf{q}|^{2}+m^{2}})\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{q}$})}{q^{\mu}q_{\mu}-m^{2}+i\epsilon}\right],\hspace{0.5cm}j=1,2,\cdots\label{eq:final_bosonic_prop}
\end{equation}
\begin{equation}
S(y,x)=
\frac{i}{2}\int\frac{d^{4}q}{(2\pi)^{4}}e^{-iq\cdot(x-y)}\left[\frac{I+i\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{q}$})}{q^{\mu}q_{\mu}-m^{2}+i\epsilon}\right],  \hspace{0.5cm} j=\frac{1}{2},\frac{3}{2},\cdots.\label{eq:final_fermionic_prop}
\end{equation}
\end{subequations}
A complete investigation on the properties of propagators would require us to compute scattering amplitudes. But without first determining the underlying symmetry, any brute force computation of the scattering amplitudes would be premature.
In spite of this, the propagators given by eqs.~(\ref{eq:final_bosonic_prop}) and (\ref{eq:final_fermionic_prop}) already provide us important information, that the fields for any spin have mass-dimension one and renormalisable self-interactions.

\subsection{A decomposition of $\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}$}

We have now determined all the proportionality constants so it is only appropriate that the important results obtained in this paper are presented in their final forms. It is certainly possible and without any loss of generality, to choose these constants at the beginning of our construction. However, this would have left the reader wondering how they are determined and whether the choices are arbitrary. Our derivation shows that, apart from trivial global phases, there are no additional freedoms.

Using the results we have obatined so far, we are able to generalise the results obtained by Ahluwalia~\cite{Ahluwalia:2013uxa} and Speran\c{c}a~\cite{Speranca:2013hqa} by decomposing $\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}$.
Starting with the Elko spin-sums, we have

\begin{equation}
\sum_{\sigma}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\dual{\chi}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=-\frac{im}{2}\left[\vartheta I-\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\right],\, j=\frac{1}{2},\frac{3}{2},\cdots
\end{equation}
and
\begin{equation}
\sum_{\sigma}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\dual{\chi}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=\frac{m}{2}[I+\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})],\, j=1,2,\cdots
\end{equation}
For the fermions, the self-conjugate and anti-self-conjugate spin-sums can be obtained by taking $\vartheta=i$ and $\vartheta=-i$ respectively. For bosons, the self-conjugate and anti-self-conjugate spin-sums are identical so we simply write them using $\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)$. There is in fact an elegant way to relate the spin-sums to $\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}\cdots\mu_{2j}}$. For this purpose, we go back to eqs.~(\ref{eq:chi_eq1}) and~(\ref{eq:chi_eq2}) which now read
\begin{subequations}
\begin{equation}
\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=(-1)^{j+\sigma}\vartheta m^{2j}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},-\sigma),\, j=\frac{1}{2},\frac{3}{2},\cdots
\end{equation}
\begin{equation}
\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=(-1)^{j+\sigma} m^{2j}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},-\sigma),\, j=1,2,\cdots
\end{equation}
\end{subequations}
Taking the Hermitian conjugation of the above expressions and using the definition of the dual coefficients, we get
\begin{subequations}
\begin{equation}
\bar{\chi}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}=-i\vartheta^{*}(-1)^{j+\sigma} m^{2j}\dual{\chi}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma),\, j=\frac{1}{2},\frac{3}{2},\cdots
\end{equation}
\begin{equation}
\bar{\chi}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}= m^{2j}\dual{\chi}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma),\, j=1,2,\cdots
\end{equation}
\end{subequations}
Multiply both sides of the equations with the self-conjugate or anti-self-conjugate functions and sum over $\sigma$, we obtain the following identities
\begin{equation}
\sum_{\sigma}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\bar{\chi}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=\frac{1}{2m^{2j-1}}
[I-i\vartheta\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})]\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}},\, j=\frac{1}{2},\frac{3}{2},\cdots
\end{equation}
and
\begin{equation}
\sum_{\sigma}\chi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\bar{\chi}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)=\frac{1}{2m^{2j-1}}
[I+\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})]\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}},\, j=1,2,\cdots.
\end{equation}

Let us now introduce the operator~\cite{Speranca:2013hqa,Ahluwalia:2013uxa}
\begin{eqnarray}
\Xi(\mbox{\boldmath$\displaystyle\mathbf{p}$})&=&\frac{i}{m}\sum_{\sigma}\left[\xi(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\bar{\xi}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)-\zeta(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\bar{\zeta}(\mbox{\boldmath$\displaystyle\mathbf{p}$},\sigma)\right]\nonumber\\
&=&\frac{i}{m^{2j}}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}.
\end{eqnarray}
This identically vanishes for bosons but it is non zero for fermions. Therefore, while the following results apply to both bosons and fermions, they are only meaningful for the later case. We can now derive some interesting identities. Consider
\begin{equation}
\Xi^{2}(\mbox{\boldmath$\displaystyle\mathbf{p}$})=-\frac{1}{m^{4j}}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})
\gamma^{\nu_{1}\cdots\nu_{2j}}p_{\nu_{1}}\cdots p_{\nu_{2j}}
\end{equation}
Using eqs.~(\ref{eq:anticommute}) and (\ref{eq:mathcalA}), it is straightforward to show that
\begin{equation}
[\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$}),\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}]=O.\label{eq:A_gamma}
\end{equation}
It follows that
\begin{equation}
\Xi^{2}(\mbox{\boldmath$\displaystyle\mathbf{p}$})=-\mathcal{A}^{2}(\mbox{\boldmath$\displaystyle\mathbf{p}$})=(-1)^{2j+1}e^{-2i(j-1/2)\phi} I
\end{equation}
and
\begin{equation}
[\Xi(\mbox{\boldmath$\displaystyle\mathbf{p}$}),\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})]=[\Xi(\mbox{\boldmath$\displaystyle\mathbf{p}$}),\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}]=O.\label{eq:Xi_A_Gamma}
\end{equation}
Equations (\ref{eq:A_gamma}) and (\ref{eq:Xi_A_Gamma}) were initially derived by Ahluwalia for $j=\frac{1}{2}$~\cite{Ahluwalia:2013uxa}. Our result shows that these identities are valid for any spin.

An interesting way to look these identities is that $\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}$ can be decomposed into a product of two commuting matrices
\begin{equation}
\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}=im^{2j}(-1)^{2j+1}e^{2i(j-1/2)\phi}\Xi(\mbox{\boldmath$\displaystyle\mathbf{p}$})\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$}).
\end{equation}
The matrices $\Xi(\mbox{\boldmath$\displaystyle\mathbf{p}$})$ and $\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{p}$})$ are non-covariant but the violating terms cancel leaving their product Lorentz-covariant. Our observation suggest that there are potentially interesting mathematics behind the decomposition.
For now, we leave their transformation properties as a future subject of research.

\section{Conclusions}

In refs.~\cite{Ahluwalia:2004ab,Ahluwalia:2004sz,Ahluwalia:2008xi,Ahluwalia:2009rh}, it was shown that contrary to conventional wisdom, there exists fermionic fields of mass-dimension one that satisfy the Klein-Gordon but not the Dirac equation. These fields constructed using spin-half Elko as expansion coefficients are Lorentz violating and have renormalisable self-interactions. The results obtained in this paper are a direct generalisation of these works. The properties of spin-half Elko are naturally inherited by their higher-spin generalisations.

It should be noted that our work may not be of pure academic interest. In the above references, it was argued that the fermions described by the massive spin-half fields of mass-dimension one are potential dark matter candidates due to renormalisable self-interactions and limited electromagnetic interactions. Since these properties remain unchanged for the higher-spin fields, it follows that they are also likely dark matter candidates. If dark matter exists, they are unlikely to be comprised entirely of one single specie. Our construct thus offers a possible solution since it contains both bosonic and fermionic particles.

While there are stringent limits on Lorentz symmetry violation, they cannot rule out our construct as these experiments were conducted using only the SM particles. There are no direct evidence suggesting that dark matter must satisfy Lorentz symmetry~\cite{Blas:2012vn}. Towards determining the symmetries of the mass-dimension one fields, the works of Ahluwalia~\cite{Ahluwalia:2013uxa}, Speran\c{c}a~\cite{Speranca:2013hqa} and our generalisation of their results on the decomposition of $\gamma^{\mu_{1}\cdots\mu_{2j}}p_{\mu_{1}}\cdots p_{\mu_{2j}}$, have provided additional insights on Lorentz violation. We believe further investigation in this direction will not only give us a deeper understanding of Elko but also elucidate on a possible connection between mass-dimension one fields and Lorentz-invariant fields.

\section*{Acknowledgement}

This research was supported by CNPq. I would like to thank D. V. Ahluwalia for numerous discussions and for communicating the result on the Fourier transform of $\mathcal{A}(\mbox{\boldmath$\displaystyle\mathbf{q}$})$.