Consider the graph given by taking a rectangular lattice with $m$ rows and $n$ columns and joining each vertex to its four nearest neighbors, where vertices on the boundary are connected periodically (for example, $(1,2)$ is connected to $(1,1),(1,3),(2,2)$ and $(m,2)$). Are the exact eigenvalues of its Laplacian matrix $L$ known?
I know how to handle the problem in the 1D case. Here except for the first and last rows, the matrix is tridiagonal, with its diagonal entries being $2$ and its superdiagonal and subdiagonal entries being $-1$. (Here I am using the positive semidefinite convention for the Laplacian, as usual in graph theory but reversed from the usual for PDE). Additionally there is a $-1$ in the top right and bottom left corners. This means that we have the circulant matrix corresponding to the vector $(2,-1,0,0,\dots,-1)$, where there are $n-3$ zeros. In this case as for any circulant matrix, the eigenvectors are the columns of the discrete Fourier transform and the eigenvalues can be read off by substitution.
This suggests that the columns of an appropriate 2D Fourier transform would be the eigenvectors of our matrix here...and in fact that is correct. To be specific, under the column major ordering of the vertices, the eigenvectors of the Laplacian matrix are the columns of $F_n \otimes F_m$ where $F_n$ and $F_m$ are the discrete Fourier matrices and $\otimes$ is the Kronecker product.
Is there an easy way to read off the actual eigenvalues from this representation of the eigenvectors? In other words, I can see that we have, for example
$$\lambda_i=\frac{\sum_{j=1}^{mn} L_{1j} (F_n \otimes F_m)_{ji}}{(F_n \otimes F_m)_{1i}}$$
but does this have some simple explicit formula?