Is every square matrix a sum of a diagonal matrix and a nilpotent matrix

The question is in the title. I am working with complex matrices. This is true for $$2\times 2$$ matrices, but becomes complicated already for $$3\times 3$$ matrices if we try to brute force it. By brute force, I mean taking an arbitrary $$3\times 3$$ matrix and imposing conditions on the diagonal entries so that the matrix is nilpotent.

A follow up question for matrices which can be expressed as a sum of diagonal and nilpotent matrices is whether there are infinitely many ways of doing this? This feels unlikely because the final answer seems to depend on solving $$n$$ independent polynomial equations in $$n$$ variables, but I might be wrong.

• If it were "diagonalisable" rather than "diagonal", then this would be a good, standard question regarding Jordan Normal Form. Are you sure it asks about diagonal matrices? I am not sure one way or another, but I am skeptical that the result is true as stated. Oct 8, 2021 at 0:59
• This isn't a question I found in a textbook, it's just something I was thinking about. I know about the "diagonalisable" result, but I noticed in examples I tried that I could in fact do diagonal, and so I got curious. Oct 8, 2021 at 1:02
• Good to know. I thought I'd check in case people went clambering up the wrong path. Oct 8, 2021 at 1:05
• This seems to be related to Hilbert's Nullstellensatz. Let $A$ be your given $n\times n$ matrix and let $D$ be the diagonal matrix with variables $X_1,\ldots,X_n$ on the diagonal. Then equation $(A-D)^n=0$ gives you $n^2$ polynomial equations. By Hilbert's theorem it suffices to show that these polynomials do not generate to whole polynomial ring as an ideal. Oct 8, 2021 at 6:06
• On the other hand, it suffices to solve $n$ polynomial equations $tr(A-D)=tr((A-D)^2)=\cdots=tr((A-D)^n)=0$ in $n$ unknowns. This seems possible over an algebraically closed field. Oct 8, 2021 at 7:27

This is always possible. Given $$A \in M_n(\mathbb{C})$$, we need to find a diagonal matrix $$D$$ such that $$A - D$$ is nilpotent. By a well-known result, this is equivalent to requiring that $$\operatorname{tr}((D - A)^k) = 0$$ for $$1 \leq k \leq n$$. Denote the diagonal entries of $$D$$ by $$x_1, \dots, x_n$$. Written explicitly, the conditions $$\operatorname{tr}(D-A) = \operatorname{tr}((D-A)^2) = \dots = \operatorname{tr}((D-A)^n) = 0$$ gives us a system of $$n$$ polynomial equations for $$x_1,\dots,x_n$$ of degrees $$1,\dots,n$$: $$\operatorname{tr}(D) - \operatorname{tr}(A) = 0, \\ \operatorname{tr}(D^2) - 2\operatorname{tr}(DA) + \operatorname{tr}(A^2) = 0, \\ \operatorname{tr}(D^3) - 3\operatorname{tr}(D^2A) + 3\operatorname{tr}(DA^2) - \operatorname{tr}(A^3) = 0, \\ \vdots \\ \operatorname{tr}(D^n) + \textrm{ lower order terms } + (-1)^n \operatorname{tr}(A^n) = 0.$$ We can homogenize the equations by adding an extra variable $$w$$ and then they will have the form: $$\operatorname{tr}(D) - w\operatorname{tr}(A) = 0, \\ \operatorname{tr}(D^2) - 2w\operatorname{tr}(DA) + w^2\operatorname{tr}(A^2) = 0, \\ \operatorname{tr}(D^3) - 3w\operatorname{tr}(D^2A) + 3w^2\operatorname{tr}(DA^2) - w^3\operatorname{tr}(A^3) = 0, \\ \vdots \\ \operatorname{tr}(D^n) + w \cdot \left( \textrm{lower order terms} \right) = 0.$$ Now, Bézout's theorem guarantees that the homogenized system has either infinitely many solutions in $$\mathbb{P}^n = \{ [x_1 : \dots : x_n : w] \, | \, (x_1, \dots, x_n,w) \in \mathbb{C}^{n+1} \setminus \{ \vec{0} \} \},$$ or, if it has finitely many solutions, then the number of solutions is $$1 \cdot 2 \cdot \dots \cdot n = n!$$, counted with multiplicities. Any solution $$[x_1 : \dots : x_n : w]$$ of the homogenized system for which $$w \neq 0$$ gives us a solution $$\left( \frac{x_1}{w}, \dots, \frac{x_n}{w} \right)$$ to the original system of equations while solutions of the form $$[x_1 \dots : x_n : 0]$$ correspond to extra solutions we have possibly added by moving to the projective space. However, if $$w = 0$$, the homogenized system of equations becomes $$\operatorname{tr}(D) = \operatorname{tr}(D^2) = \dots = \operatorname{tr}(D^n) = 0.$$ This implies that $$D$$ is both nilpotent and diagonal so $$D = 0$$ which means that the homogenized system has no solutions of the form $$[x_1 : \dots : x_n : 0]$$. Hence, we haven't added any solutions by moving to the projective space and the bottom line is that the original system has either infinitely many solutions or $$n!$$ solutions, counted with multiplicities. In any case, such a $$D$$ always exist.

1. If you only care about the existence of $$D$$, you only need to know that the zero set of $$n$$ homogeneous equations in $$\mathbb{P}^n$$ is always non-empty. This is certainly weaker than Bézout's theorem but I'm not an expert and have no idea if this result by itself has an "elementary" proof.
2. The argument depends heavily on the fact that $$\mathbb{C}$$ is algebraically closed. You can verify by a direct calculation that the result is false for $$2 \times 2$$ real matrices.
3. Generically, one expects the number of different ways to write $$A$$ as a sum of nilpotent and a diagonal matrix to be $$n!$$. It can certainly be less (this already happens for $$2 \times 2$$ matrices and corresponds to solutions with multiplicities) but I suspect that it will always be finite. Geometrically, this should correspond to the fact that the $$n$$ different hypersurfaces defined by the equations have "no common components" but I don't know how to justify this rigourously.
• Brilliant, thanks! I am quite intrigued by the $n!$, wondering if there’s a way to relate all the solutions or generate some other solutions if we already know some. Oct 11, 2021 at 8:26
• Regarding your point (3): Any positive dimensional subvariety of $\mathbb{P}^n$ intersects the hyperplane at infinity. So, by showing there are no solutions at infinity, you have shown that there are only finitely many solutions. PS: Very nice answer! Oct 11, 2021 at 17:38
• @MathManiac: I'm not sure, this depends on whether the resulting system of equations has a nice symmetry and this depends heavily on the form of $A$. For example, if $A = \begin{pmatrix} 1 & 1 & 1 \\ 1 & 1 & 1 \\ 1 & 1 & 1 \end{pmatrix}$ you can check that the resulting system of equations is equivalent to $x + y + z = 3, x^2 + y^2 + z^2 = -3, x^3 + y^3 + z^3 = -9$. Since all the polynomials involved are symmetric, if $(a,b,c)$ is a solution then any permutation of it will give you another solution so if $a,b,c$ are distinct (which is indeed the case) then this gives you all $3!$ solutions. Oct 12, 2021 at 12:16
• Low tech explanation: Let $X \subset \mathbb{P}^n$ have dimension $d>0$ and let $H \subset \mathbb{P}^n$ be a projective hyperplane. Let $\tilde{X}$ and $\tilde{H}$, in $\mathbb{A}^{n+1}$, be the cones on $X$ and $H$. Then $\tilde{X}$ has dimension $d+1$ so $\tilde{X} \cap \tilde{H}$ either has dimension $\geq d$ or else is empty. But $0 \in \tilde{X} \cap \tilde{H}$, so the intersection isn't empty. So $\tilde{X} \cap \tilde{H}$ contains a point other than $0$, and $X \cap H$ is nonempty. Oct 12, 2021 at 13:07