The Pauli matrices are $$e = \begin{bmatrix}1&&0\\0&&1\end{bmatrix},\ X = \begin{bmatrix}0&&1\\1&&0\end{bmatrix},\ Y = i\begin{bmatrix}0&&-1\\1&&0\end{bmatrix},\ Z = \begin{bmatrix}1&&0\\0&&-1\end{bmatrix}$$ I want to create a group that contains all these matrices. I started drawing a Cayley Table as
| * | e | X | Y | Z |
|---|---|---|---|---|
| e | e | X | Y | Z |
| X | X | e | | |
| Y | Y | | e | |
| Z | Z | | | e |
At this point I realize that XY is not equal to YX, i.e. the group is non abelian. So, the rest of the group should have distinct element i.e. 6 more elements, i.e, XY, YX, XZ, ZX, ZY and YZ. Now the total number of elements come out to be 10 now, i.e. e, X, Y, Z, XY, YX, XZ, ZX, ZY and YZ.
But the Pauli matrices is a group with 16 elements. I don't know what I am missing. I am a beginner in Group theory and I am making a colossal mistake somewhere. Please help.
Yes, the total number of elements is $16$. Consider the elements $XYZ$ and its six permutations.