9 Apr 2026

While studying for my matrix algebra midterm I came across the exercise of proving Fuglede's theorem for finite-dimensional, complex, inner product spaces, which I found particularly exciting, and whose solution I will share in a moment.

But first, even before stating the result, I will show a lemma that will be employed throughout the proof:

Let $A \in \mathscr M_n(\mathbf C)$ be a normal matrix. There exists a polynomial $p \in \mathbf C[x]$ such that $$ p(A) = \adj A. $$
Arguably the most important fact about normal matrices is the spectral theorem, which tells that they are unitarily diagonalizable; that is, there exist $U, \Lambda \in \mathscr M_n(\mathbf C)$ with $U \adj U = \adj U U = \mathbf 1$ and $\Lambda$ diagonal such that $$ A = U\Lambda \adj U; $$ hence, we look for a polynomial $p$ such that $$ \begin{align*} p(A) &= Up(\Lambda) \adj U = U \begin{pmatrix} p(\lambda_1) & & &\\ & p(\lambda_2) & &\\ & & \ddots &\\ & & & p(\lambda_n) \end{pmatrix} \adj U = U \begin{pmatrix} \bar \lambda_1 & & &\\ & \bar \lambda_2 & &\\ & & \ddots &\\ & & & \bar \lambda_n \end{pmatrix} \adj U = \adj A\\ \implies p(\lambda_i) &= \bar \lambda_i. \end{align*} $$ Let now, without loss of generality, $\lambda_1, \lambda_2, \ldots, \lambda_k$ be the $k$ distinct eigenvalues of $A$, and set $$ p(t) = q_0 + q_1 t + q_2 t^2 + \cdots + q_{k-1} t^{k-1}; $$ we want $p$ to be the interpolating polynomial of the points $(\lambda_i, \bar \lambda_i)$, so its coefficients are given by the equation $$ \begin{pmatrix} 1 & \lambda_1 & \lambda_1^2 & \cdots & \lambda_1^{k-1}\\ 1 & \lambda_2 & \lambda_2^2 & \cdots & \lambda_2^{k-1}\\ 1 & \lambda_3 & \lambda_3^2 & \cdots & \lambda_3^{k-1}\\ \vdots & \vdots & \vdots & \ddots & \vdots\\ 1 & \lambda_k & \lambda_k^2 & \cdots & \lambda_k^{k-1}\\ \end{pmatrix} \begin{pmatrix} q_0\\ q_1\\ q_2\\ \vdots\\ q_{k-1} \end{pmatrix} = \begin{pmatrix} \bar \lambda_1\\ \bar \lambda_2\\ \bar \lambda_3\\ \vdots\\ \bar \lambda_k \end{pmatrix}, $$ which has a solution because the system matrix is a Vandermonde matrix with distinct entries, and so it is invertible.

With this tool in mind, we can now tackle the simplified version of Fuglede's theorem:

Let $N \in \mathscr M_n(\mathbf C)$ be a normal matrix, and $A \in \mathscr M_n(\mathbf C)$ an arbitrary complex matrix. If $A$ commutes with $N$, so does with $\adj N$, i.e. $$ NA = AN \implies \adj N A = A \adj N. $$
Notice first that $A$ commutes with any power of $N$; we can prove this claim inductively, where the base case is the initial assumption, and for the inductive step $$ N^{n+1} A = N N^n A = N A N^n = A N N^n = A N^{n+1}. $$ Extending this result, because matrix product is distributive over addition ($\mathscr M_n(\mathbf C)$ is a ring) $A$ commutes with any linear combination of powers of $N$, in other words, with any polynomial acting on $N$.
Finally, invoking the previous lemma settles the proof: $$ A \adj N = A p(N) = p(N) A = \adj N A. $$

And that is about it. I think this is a nice result, so I wanted to share it here.