Es una buena pregunta. Para responderla vamos a empezar con el álgebra de Clifford generada por \gamma matrices. \begin{equation} \gamma_{\mu}\gamma_{\nu}+ \gamma_{\mu}\gamma_{\nu}=2\eta_{\mu\nu} \end{equation} con \mu,\nu=0,1,2,\cdots N con la firma métrica \eta_{\mu\nu}=\text{diag}(+,-,-,-,\cdots,-) . Utilizando I y \gamma_{\mu} podemos construir un conjunto de matrices como las siguientes \begin{equation} I, \gamma_{\mu},\gamma_{\mu}\gamma_{\nu}\quad(\mu<\nu), \gamma_{\mu}\gamma_{\nu}\gamma_{\lambda}\quad(\mu<\nu<\lambda),\cdots,\gamma_{1}\gamma_{2}\cdots\gamma_{N}. \end{equation}
Hay \begin{equation} \sum_{p=0}^{N}\binom{N}{p} = 2^{N} \end{equation} tales matrices. Llamémoslas \Gamma_{A} , donde A corre de 0 a 2^{N}-1 . Ahora dejemos que \gamma_{\mu} son d\times d matrices irreducibles dimensionales. Nuestro objetivo es encontrar una relación entre d y N . Para ello vamos a definir una matriz \begin{equation} S = \sum_{A=0}^{2^N-1}(\Gamma_{A})^{-1}Y\Gamma_{A} \end{equation} . Donde Y es algo arbitrario d\times d matriz. Se deduce que \begin{equation} (\Gamma_{B})^{-1}S\Gamma_{B} = \sum_{A=0}^{2^N-1}(\Gamma_{A}\Gamma_{B})^{-1}Y\Gamma_{A}\Gamma_{B} =\sum_{C=0}^{2^N-1}(\Gamma_{C})^{-1}Y\Gamma_{C}=S \end{equation} Donde hemos utilizado \Gamma_{A}\Gamma_{B}=\epsilon_{AB}\Gamma_{C} con \epsilon_{AB}^{2}=1
Por lo tanto, \begin{equation}S\Gamma_{A}=\Gamma_{A}S\end{equation} Desde S conmuta con todas las matrices del conjunto, por el lema de Schur concluimos que S debe ser proporcional a la matriz identidad para que podamos escribir \begin{equation} S = \sum_{A=0}^{2^N-1}(\Gamma_{A})^{-1}Y\Gamma_{A} = \lambda I \end{equation}
Tomando la traza obtenemos \begin{eqnarray} \text{Tr} S & = & \sum_{A=0}^{2^N-1} \text{Tr} Y = \lambda d\\ \Rightarrow \lambda & = & \frac{2^{N}}{d}\text{Tr} Y \end{eqnarray} o \begin{equation} \sum_{A=0}^{2^N-1}(\Gamma_{A})^{-1}Y\Gamma_{A} = \frac{2^{N}}{d}\text{Tr} Y \end{equation}
Tomando la (j; m) El elemento de la matriz de ambos lados de la última ecuación da como resultado \begin{equation} \sum_{A=0}^{2^N-1}((\Gamma_{A})^{-1})_{jk}(\Gamma_{A})_{km} = \frac{2^{N}}{d}\delta_{jm} \delta_{kl} \end{equation} donde j; k; l; m = 1; 2;\cdots; d y hemos utilizado el hecho de que Y es un arbitrario d \times d matriz. Si establecemos j = k; l = m y sumar sobre estos dos índices, lo que da \begin{equation} \sum_{A=0}^{2^N-1} \text{Tr}[(\Gamma_{A})^{-1}] \text{Tr}[\Gamma_{A}] = 2^{N}\end{equation} Hay dos casos a considerar, a saber, N incluso y N impar. Para N = 2M (incluso), \text{Tr} \Gamma_{A} = 0 excepto en el caso de \Gamma_{0} = 1 para lo cual \text{Tr} \Gamma_{0} = d . Lo que da \begin{equation} d^2 = 2^N\qquad \text{or} \quad \boxed{d = 2^{N/2}} \end{equation} Este es el resultado principal. Para un espacio de Minkowski de cuatro dimensiones el tiempo N=4 cosequntly la dimensión de la representación irreducible es d = 2^{4/2} =4 .