The Ehrhart series of magic squares of orders seven and eight
Let $\mathrm{IMS}_n(m)$ count the $n\times n$ nonnegative integer matrices whose row sums, column sums, main-diagonal sum, and antidiagonal sum are all $m$. We determine the Ehrhart series $F_n(q)=\sum_{m\geq0}\mathrm{IMS}_n(m)q^m$ as reduced rational functions for $n=7$ and $n=8$. Their numerator--denominator degrees are respectively $(366,373)$ and $(540,548)$. Both numerators have positive integer coefficients, are palindromic and strictly unimodal. The proofs share one finite architecture: a signed SimpCone decomposition is evaluated by quotient characters over finite fields, a certified common denominator and Ehrhart reciprocity reduce the rational identity to finitely many coefficients, an explicit counting bound lifts modular congruences to integer equalities, and exact gcd computations prove reducedness. For order eight, a face-index pole certificate gives a degree-$598$ common denominator without enumerating the full face lattice, leaving $296$ independent coefficients in degrees $0$ through $295$. Once the candidate rational function is known, the first six production primes certify this finite prefix by the same bounded-coefficient argument used for order seven.