A quaternionic matrix is a matrix whose elements are quaternions.

The quaternions form a noncommutative ring, and therefore addition and multiplication can be defined for quaternionic matrices as for matrices over any ring.

Addition. The sum of two quaternionic matrices A and B is defined in the usual way by element-wise addition:

${\displaystyle (A+B)_{ij}=A_{ij}+B_{ij}.\,}$

Multiplication. The product of two quaternionic matrices A and B also follows the usual definition for matrix multiplication. For it to be defined, the number of columns of A must equal the number of rows of B. Then the entry in the ith row and jth column of the product is the dot product of the ith row of the first matrix with the jth column of the second matrix. Specifically:

${\displaystyle (AB)_{ij}=\sum _{s}A_{is}B_{sj}.\,}$

For example, for

${\displaystyle U={\begin{pmatrix}u_{11}&u_{12}\\u_{21}&u_{22}\\\end{pmatrix}},\quad V={\begin{pmatrix}v_{11}&v_{12}\\v_{21}&v_{22}\\\end{pmatrix}},}$

the product is

${\displaystyle UV={\begin{pmatrix}u_{11}v_{11}+u_{12}v_{21}&u_{11}v_{12}+u_{12}v_{22}\\u_{21}v_{11}+u_{22}v_{21}&u_{21}v_{12}+u_{22}v_{22}\\\end{pmatrix}}.}$

Since quaternionic multiplication is noncommutative, care must be taken to preserve the order of the factors when computing the product of matrices.

The identity for this multiplication is, as expected, the diagonal matrix I = diag(1, 1, ... , 1). Multiplication follows the usual laws of associativity and distributivity. The trace of a matrix is defined as the sum of the diagonal elements, but in general

${\displaystyle \operatorname {trace} (AB)\neq \operatorname {trace} (BA).}$

Left scalar multiplication, and right scalar multiplication are defined by

${\displaystyle (cA)_{ij}=cA_{ij},\qquad (Ac)_{ij}=A_{ij}c.\,}$

Again, since multiplication is not commutative some care must be taken in the order of the factors.^[1]

There is no natural way to define a determinant for (square) quaternionic matrices so that the values of the determinant are quaternions.^[2] Complex valued determinants can be defined however.^[3] The quaternion a + bi + cj + dk can be represented as the 2×2 complex matrix

${\displaystyle {\begin{bmatrix}~~a+bi&c+di\\-c+di&a-bi\end{bmatrix}}.}$

This defines a map Ψ_mn from the m by n quaternionic matrices to the 2m by 2n complex matrices by replacing each entry in the quaternionic matrix by its 2 by 2 complex representation. The complex valued determinant of a square quaternionic matrix A is then defined as det(Ψ(A)). Many of the usual laws for determinants hold; in particular, an n by n matrix is invertible if and only if its determinant is nonzero.

Quaternionic matrices are used in quantum mechanics^[4] and in the treatment of multibody problems.^[5]