Let
$$M = \left(\begin{matrix} a & a &... & a\\a & a & ... & a\\\vdots & \vdots & \ddots & \vdots\\a & a & ... &a\end{matrix}\right)$$
For unitary diagonalization: $U^{\dagger}MU =D$. What are $U$ and $D$?
My generalization from examples tell me that
$$D = \left(\begin{matrix} Na & 0 &... & 0\\0 & 0 & ... & 0\\\vdots & \vdots & ... & \vdots\\0 & 0 & ... &0\end{matrix}\right)$$
where $N$ is the dimension of $M$. $Na$ is the only nonzero eigenvalue of $M$! Then for $U$,
$$U = \left(\begin{matrix} \frac{1}{\sqrt{N}} & 0 &... & 0\\\frac{1}{\sqrt{N}} & 0 & ... & 0\\\vdots & \vdots & ... & \vdots\\\frac{1}{\sqrt{N}} & 0 & ... &0\end{matrix}\right)$$
But how to calculate the rest of the columns to make it unitary?
The nicest matrix to use here (in my opinion) is the DFT matrix. Note that it's a unitary matrix whose first column is $\frac 1{\sqrt{N}}(1,\dots,1)$, which is all we really need here.