Suppose we are given a family of $k$ vectors $(\vec{v}_1,\dots,\vec{v}_k)$ in $\mathbb{R}^k$ such that $\vec{v}_i\cdot\vec{v}_i=1$ for all $1\leq i\leq k$ and $||\vec{v}_i-\vec{v}_j||=1$ for all $i\neq j$. How do we compute the eigen values and their multiplicity of the Gram matrix
$$GM(\vec{v}_1,\dots,\vec{v}_k)= \begin{pmatrix} \vec{v}_1\cdot\vec{v}_1 & \vec{v}_1\cdot\vec{v}_2 &\dots &\vec{v}_1\cdot\vec{v}_k\\ \vec{v}_2\cdot\vec{v}_1 & \vec{v}_2\cdot\vec{v}_2 & \dots &\vec{v}_2\cdot\vec{v}_k\\ \vdots & \vdots &\ddots &\vdots\\ \vec{v}_k\cdot\vec{v}_1 & \vec{v}_k\cdot\vec{v}_2&\dots&\vec{v}_k\cdot\vec{v}_k \end{pmatrix} $$
I understand that, by hypothesis, we have $\vec{v}_i\cdot\vec{v}_i$ for all $1\leq i\leq k$ which in turn implies $$\begin{align*} ||\vec{v}_i-\vec{v}_j||^2 =\vec{v}_i\cdot\vec{v}_i - 2(\vec{v}_i\cdot\vec{v}_j)+\vec{v}_j\cdot\vec{v}_j =2-2(\vec{v}_i\cdot\vec{v}_j) \end{align*} $$
so that $$\vec{v}_i\cdot\vec{v}_j=\frac{1}{2}$$ for all $i\neq j$.
This means $$GM(\vec{v}_1,\dots,\vec{v}_k)= \begin{pmatrix} 1 & \frac{1}{2} &\dots & \frac{1}{2}\\ \frac{1}{2} & 1 & \dots &\frac{1}{2}\\ \vdots & \vdots &\ddots &\vdots\\ \frac{1}{2} & \frac{1}{2}&\dots&1 \end{pmatrix} $$
but I am unsure how to continue from here. More generally, is there a method to computing the eigenvales/characteristic polynomials of large symmetric matrices?
Any help would be appreciated. Thank you. (For context, see this paper on page four).