If $Q$ is the sum of squares quadratic form $\sum_1^n x_i^2$ over some lattice, then $r_Q(m)$, the number of representations of an integer $m$ by $Q$ (order/sign matter) is sometimes given in a nice formula, as in the case with Jacobi's formula in the case of $n=4$. Can we say something moderate about how $r_Q$ grows? It seems like it shouldn't grow faster than polynomially, but I am struggling to see why in a rigorous way. Since $\sum_1^n x_i^2=m$ is the equation of a sphere (or I suppose an ellipsoid if we transform our lattice to $Z^n$ under a change of variable?), we should be able to bound such solutions by some function of the surface area since there should be some small upper bound for the proportion of lattice points over the surface area and the surface area is a polynomial function of the radius.
I would like to make the argument a bit more rigorous, any advice?
Edit: To be clear, I'm interested in how $r_Q(m)$ grows with $m$.
For each solution to $\sum_{i=1}^n x_i^2 = m$, you can place an open cube of side-length $1$ centered on that solution. Doing so, you obtain a family of disjoint open cubes covering the solutions, with each open cube having volume $1$. So your next move should be to bound the volume that you can possibly cover with those cubes, polynomially in $m$.
Since every cube centered on $x$ is included in the open ball of radius $\sqrt n/2$ centered on $x$ (using the euclidean metric), and $d(0,x) = \sqrt m$ whenever $x$ is a representation of $m$, you obtain with the triangular inequality that all those cubes are included in $\{y \in \Bbb R^n, \sqrt m - \sqrt n/2 < d(0,y) < \sqrt m + \sqrt n/2 \} = B(0,\sqrt m + \sqrt n/2) - B(0,\sqrt m - \sqrt n/2)$
Since the volume of a ball of radius $R$ is $k_nR^n$ for some constant $k_n$, you get that the volume of this domain is $k_n((\sqrt m + \sqrt n/2)^n - (\sqrt m - \sqrt n/2)^n) \sim k_n(n^{3/2}m^{(n-1)/2}) $ (you can also argue intuitively that this should behave like $\sqrt n$ times the area of the sphere of radius $\sqrt m$, and again find a volume bounded by a $K_n m^{(n-1)/2}$).
Since the number of representations of $m$ is bounded by the volume of this domain, and since this volume is a $O(m^{(n-1)/2})$, that number $r_n(m)$ grows polynomially at worst.
You can also use this reasoning with volumes to show that the degree $(m-1)/2$ is the smallest possible for this kind of result : if it wasn't, then putting cubes around every solution of $\sum x_i < R^2$ would tell you that the volume of the sphere of radius $R$ would grow smaller than $R^n$, which is absurd.