Mallick and Mishra [‘Links between the solvability of matrix and scalar Diophantine equations’, Bull. Aust. Math. Soc., 10.1017/S0004972726101087] observed that solutions of the matrix Diophantine equation $X^n + Y^n = 2Z^n$
, where $X, Y, Z$
are $2 \times 2$
matrices in a certain matrix class $G_2(d, l)$
, correspond to solutions of the scalar equation $\alpha ^n + \beta ^n = 2\gamma ^n$
over $\mathbb {Z}[\sqrt {d}]$
and assert that there are only trivial solutions. We show that this fails when $d = -1$
with $4 \mid n$
and $d = -3$
with $3 \mid n$
, and determine the values of $l$
that admit solutions. We extend the analysis to $k \times k$
matrices, and show that for every totally real field $K$
and every even $n$
, all solutions are trivial. We obtain some results for odd $n$
, but the problem for odd $n$
in totally real fields of degree $k \geqslant 3$
remains open.