Si $b = 0$ o $c = 0$ entonces la matriz es triangular y los valores propios son sólo las entradas diagonales $a^2, d^2$ (que en particular son racionales). Por lo tanto, cualquiera de los $12$ opciones con $\{a, b, c, d\} = \{0, 1, 2, 3\}$ y $b = 0$ o $c = 0$ da una solución, y estos ejemplos evidentemente minimizan $\max\{a, b, c, d\}$ entre las soluciones no negativas.
Si excluimos $0 \in \{a, b, c, d\}$ para evitar estas soluciones triviales, las soluciones mínimas y sus valores propios son \begin{array}{cc}\hline (a, b, c, d) & \lambda \\ \hline(1, 2, 5, 4) & -4, 21 \\ (4, 2, 3, 5) & 13, 28 \\ \hline \end{array} y los seis ejemplos obtenidos a partir de ellos utilizando las simetrías evidentes $a \leftrightarrow d$ y $b \leftrightarrow c$ .
Para el $3 \times 3$ caso ya la matriz $$\pmatrix{0^2 & 1^2 & 2^2 \\ 3^2 & 4^2 & 5^2 \\ 6^2 & 7^2 & 8^2}$$ tiene valores propios irracionales: Su polinomio característico, $c(t) = t^3 - 80 t^2 - 354 t + 216$ tiene discriminante positivo, por lo que tiene tres raíces reales. Por otro lado, $c(t) \equiv t^3 + t + 1 \pmod 5$ pero este último polinomio no tiene raíces módulo $5$ Por lo tanto $c(t)$ es irreducible sobre $\Bbb Q$ es decir, sus tres raíces reales son irracionales.
En un sentido que pueda precisarse, la mayoría de los polinomios racionales no tienen todas las raíces racionales, y veo pocas razones para esperar que los polinomios característicos de matrices con entradas cuadradas (distintas) sean especiales en este sentido, por lo que quizá sea más interesante preguntar por un ejemplo cuyos valores propios son racional (y, por tanto, integral). Un rápido script de Maple (transcrito a continuación) encuentra muchos ejemplos con entradas $0^2, \ldots, 8^2$ . El primero de ellos lexicográficamente es $$\pmatrix{0^2 & 2^2 & 3^2 \\ 5^2 & 8^2 & 1^2 \\ 7^2 & 6^2 & 4^2} , \quad \textrm{which has eigenvalues} \quad {-13}, 24, 69 .$$ Es poco frecuente incluso entre las matrices con entradas $0^2, \ldots, 8^2$ para tener todos los valores propios racionales: Sólo ocurre para $252$ de la $9! = 362880$ casos, $180$ de los cuales no tienen valores propios nulos.
restart;
with(combinat): with(LinearAlgebra):
m := 3;
N := 9;
for numberSet in choose(N, m^2) do
shifted := map(U -> U - 1, numberSet);
print([shifted]);
for ordering in permute(shifted) do
map(U -> U^2, ordering);
A := convert([seq(%[((i - 1) * m + 1)..(i * m)], i=1..m)], Matrix);
c := CharacteristicPolynomial(A, t);
if (convert(map(degree, map(U -> U[1], factors(c)[2]), t), set) = {1}) then
print(ordering, A, solve(c));
fi:
od:
od:
Aquí la constante $m$ es el tamaño de la matriz, $[0, \ldots, N - 1]$ es el intervalo del que se eligen los números al cuadrado.