electromagnetism - Can I transform electromagnetic tensors by matrix multiplication?

I know that the eletromagnetic field tensor $F^{\mu\nu}$, can be transfomed to another reference frame by $$F^{\alpha\beta} = \varLambda^{\alpha}_{\mu}\varLambda^{\beta}_{\nu}F^{\mu\nu}$$

Since these tensors can be represented by matrices so I thought that I could represent the electromagnetic field tensor in another inertial reference frame by doing matrix multiplications, but I ended up with:

$$F^{\alpha\beta}= \begin{bmatrix} \gamma & -\gamma\beta & 0 & 0 \\ -\gamma\beta & \gamma & 0 & 0\\ 0 & 0 & 1 & 0\\ 0 & 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} \gamma & -\gamma\beta & 0 & 0 \\ -\gamma\beta & \gamma & 0 & 0\\ 0 & 0 & 1 & 0\\ 0 & 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} 0 & -\frac{1}{c}E_x & -\frac{1}{c}E_y & -\frac{1}{c}E_z \\ \frac{1}{c}E_x & 0 & -B_z & B_y\\ \frac{1}{c}E_y & B_z &0 & -B_x\\ \frac{1}{c}E_z & -B_y & B_x & 0 \end{bmatrix} =$$ $$=\begin{bmatrix} -2\frac{\gamma^2\beta}{c}E_x & -\frac{\gamma^2}{c}E_x(1+\beta^2)& -\frac{\gamma^2}{c}E_y(1+\beta^2)+2\gamma^2\beta^2B_z & -\frac{\gamma^2}{c}E_z(1+\beta^2)-2\gamma^2\beta B_y \\ \frac{1}{c}E_x & 2\frac{\gamma^2\beta}{c}E_x & 2\frac{\gamma^2\beta}{c}E_y-\gamma^2B_z(\beta^2+1) & 2\frac{\gamma^2\beta}{c}E_z+\gamma^2B_y(\beta^2+1)\\ \frac{1}{c}E_y & B_z &0 & -B_x\\ \frac{1}{c}E_z & -B_y & B_x & 0 \end{bmatrix}$$ But as it can easily be seen, this matix is not anti-symmetric as an eletromagnetic field tensor should be, by definition, so does this mean that I cannot use matrix multiplication on tensors or does it mean that I made a mistake somewhere in the calculations of the matrix products?

Answer

^{$$\begin{align} F^{\alpha\beta} & = \begin{bmatrix} \hphantom{-}\gamma & -\gamma\beta & \hphantom{-} 0 & \hphantom{-}0\hphantom{-} \vphantom{\dfrac12}\\ -\gamma\beta & \hphantom{-}\gamma & \hphantom{-}0 & \hphantom{-} 0 \hphantom{-} \vphantom{\dfrac12}\\ \hphantom{-}0 & \hphantom{-}0 & \hphantom{-}1 & \hphantom{-} 0 \hphantom{-} \vphantom{\dfrac12}\\ \hphantom{-} 0 & \hphantom{-} 0 & \hphantom{-} 0 & \hphantom{-} 1 \hphantom{-} \vphantom{\dfrac12} \end{bmatrix} \begin{bmatrix} 0 & -E_x & -E_y & -E_z \vphantom{\dfrac12}\\ E_x & \hphantom{-} 0 & -cB_z & \hphantom{-} cB_y \vphantom{\dfrac12} \\ E_y & \hphantom{-} cB_z & \hphantom{-} 0 & -cB_x \vphantom{\dfrac12} \\ E_z & -cB_y & \hphantom{-} cB_x & \hphantom{-} 0 \vphantom{\dfrac12} \end{bmatrix} \begin{bmatrix} \hphantom{-}\gamma & -\gamma\beta & \hphantom{-} 0 & \hphantom{-}0\hphantom{-} \vphantom{\dfrac12}\\ -\gamma\beta & \hphantom{-}\gamma & \hphantom{-}0 & \hphantom{-} 0 \hphantom{-} \vphantom{\dfrac12}\\ \hphantom{-}0 & \hphantom{-}0 & \hphantom{-}1 & \hphantom{-} 0 \hphantom{-} \vphantom{\dfrac12}\\ \hphantom{-} 0 & \hphantom{-} 0 & \hphantom{-} 0 & \hphantom{-} 1 \hphantom{-} \vphantom{\dfrac12} \end{bmatrix} \nonumber\\ & = \begin{bmatrix} \!\!-\gamma\beta E_x &\!\! -\gamma E_x & -\gamma(E_y-\beta cB_z) & -\gamma(E_z+\beta cB_y)\vphantom{\dfrac12} \hphantom{-} \\ \hphantom{-}\gamma E_x & \!\gamma\beta E_x & \hphantom{-}\gamma(\beta E_y- cB_z) & \hphantom{-}\gamma(\beta E_z+cB_y)\vphantom{\dfrac12}\hphantom{-}\\ \hphantom{\gamma\beta} E_y & \hphantom{-} cB_z &0 & -cB_x\vphantom{\dfrac12}\hphantom{-}\\ \hphantom{\gamma\beta} E_z & -cB_y & cB_x & \hphantom{-}0\vphantom{\dfrac12}\hphantom{-} \end{bmatrix} \begin{bmatrix} \hphantom{-}\gamma & -\gamma\beta & \hphantom{-} 0 & \hphantom{-}0\hphantom{-} \vphantom{\dfrac12}\\ -\gamma\beta & \hphantom{-}\gamma & \hphantom{-}0 & \hphantom{-} 0 \hphantom{-} \vphantom{\dfrac12}\\ \hphantom{-}0 & \hphantom{-}0 & \hphantom{-}1 & \hphantom{-} 0 \hphantom{-} \vphantom{\dfrac12}\\ \hphantom{-} 0 & \hphantom{-} 0 & \hphantom{-} 0 & \hphantom{-} 1 \hphantom{-} \vphantom{\dfrac12} \end{bmatrix} \nonumber\\ & =\begin{bmatrix} \hphantom{-}0 & -E_x & -\gamma(E_y-\beta cB_z) &-\gamma(E_z+\beta cB_y)\hphantom{-}\vphantom{\dfrac12} \\ \hphantom{-} E_x & 0 & \hphantom{-}\gamma(\beta E_y- cB_z) & \hphantom{-}\gamma(\beta E_z+cB_y)\hphantom{-}\vphantom{\dfrac12}\\ \hphantom{-}\gamma(E_y-\beta cB_z) & -\gamma(\beta E_y- cB_z) &\hphantom{-} 0 & -cB_x\hphantom{-}\vphantom{\dfrac12}\\ \hphantom{-}\gamma(E_z+\beta cB_y) & -\gamma(\beta E_z+cB_y) & \hphantom{-} cB_x & \hphantom{-}0\hphantom{-}\vphantom{\dfrac12} \end{bmatrix} \nonumber\\ & =\begin{bmatrix} \hphantom{-}0 & -E'_x & -E'_y & -E'_z \hphantom{-}\vphantom{\dfrac12}\\ \hphantom{-}E'_x & \hphantom{-} 0 & -cB'_z & \hphantom{-} cB'_y \hphantom{-}\vphantom{\dfrac12} \\ \hphantom{-}E'_y & \hphantom{-} cB'_z & \hphantom{-} 0 & -cB'_x \hphantom{-}\vphantom{\dfrac12} \\ \hphantom{-}E'_z & -cB'_y & \hphantom{-} cB'_x & \hphantom{-} 0 \hphantom{-}\vphantom{\dfrac12} \end{bmatrix} \tag{01} \end{align}$$} Since $\:\beta=\upsilon/c$ $$\begin{align} E'_x & = \hphantom{-}E_x \tag{02.x}\vphantom{\dfrac12}\\ E'_y & = \gamma(E_y-\upsilon B_z) \tag{02.y}\vphantom{\dfrac12}\\ E'_z & = \gamma(E_z+\upsilon B_y) \tag{02.z}\vphantom{\dfrac12}\\ B'_x & = \hphantom{-}B_x \tag{03.x}\vphantom{\dfrac12}\\ B'_y & = \gamma(B_y+\dfrac{\upsilon}{c^{2}} E_z) \tag{03.y}\\ B'_z & = \gamma(B_z-\dfrac{\upsilon}{c^{2}} E_y) \tag{03.z} \end{align}$$

---

^{For Your Information :}

^{The equations of a more general Lorentz Transformation between two systems $\:\mathrm S(\mathbf{x},t)\:$ and $\:\mathrm S'(\mathbf{x}',t')\:$, the latter translating with constant velocity $\:\mathbf{v}\!=\!\upsilon\mathbf{n}\,,\Vert\mathbf{n}\Vert=1\,, \upsilon \in (-c,+c)$, with respect to the former, are : $$\begin{align} \mathbf{x}'& \!=\!\mathbf{x}\!\boldsymbol{+}\!(\gamma\!\boldsymbol{-}\!1)(\mathbf{n}\boldsymbol{\cdot}\mathbf{x})\mathbf{n}\!\boldsymbol{-}\!\gamma\mathbf{v}t \tag{ft-01a}\\ t' & \!=\! \gamma\left(t\!\boldsymbol{-}\!\dfrac{\mathbf{v}\boldsymbol{\cdot} \mathbf{x}}{c^{2}}\right) \tag{ft-01b}\\ \gamma & \!=\!\left(1\!\boldsymbol{-}\!\dfrac{\upsilon^{2}}{c^{2}}\right)^{\boldsymbol{-}1/2} \tag{ft-01c} \end{align}$$ see Figure.(1)}

^{Under (ft-01) the vectors $\:\mathbf{E},\mathbf{B}\:$ of the electromagnetic field in empty space are transformed as follows :
$$\begin{align} \mathbf{E}'& \!=\!\gamma\mathbf{E}\!\boldsymbol{-}\!(\gamma\!\boldsymbol{-}\!1)(\mathbf{n}\boldsymbol{\cdot}\mathbf{E})\mathbf{n}\boldsymbol{+}\:\gamma\left(\mathbf{v}\boldsymbol{\times}\mathbf{B}\right) \tag{ft-02a}\\ \mathbf{B}'& \!=\!\gamma\mathbf{B}\!\boldsymbol{-}\!(\gamma\!\boldsymbol{-}\!1)(\mathbf{n}\boldsymbol{\cdot}\mathbf{B})\mathbf{n}\!\boldsymbol{-}\!\dfrac{\gamma}{c^{2}}\left(\mathbf{v}\boldsymbol{\times}\mathbf{E}\right) \tag{ft-02b} \end{align}$$ Equations (02),(03) are a special case of (ft-02) for $\:\mathbf{n}=(1,0,0)$.}

---

^{(1) See a 3D version of this Figure here : Figure 3D version}