We consider the question of when the Laplace eigenfunctions on an arbitrary flat torus T Γ := R d /Γ are flexible enough to approximate, over the natural length scale of order 1/ √ λ where λ ≫ 1 is the eigenvalue, an arbitary solution of the Helmholtz equation ∆h + h = 0 on R d . This problem is motivated by the fact that, by the asymptotics for the local Weyl law, "approximate Laplace eigenfunctions" do have this approximation property on any compact Riemannian manifold. What we find is that the answer depends solely on the arithmetic properties of the spectrum. Specifically, recall that the eigenvalues of T Γ are of the form λ k = Q Γ (k), where Q Γ is a quadratic form and k ∈ Z d . Our main result is that the eigenfunctions of T Γ have the desired approximation property if and only Q Γ is a multiple of a quadratic form with integer coefficients. In particular, the set of lattices Γ for which this approximation property holds has measure zero but includes all rational lattices. A consequence of this fact is that when Q Γ is a multiple of a quadratic form with integer coefficients, Laplace eigenfunctions exhibit an extremely flexible behavior over scales of order 1/ √ λ. In particular, there are eigenfunctions of arbitrarily high energy that exhibit nodal components diffeomorphic to any compact hypersurface of diameter O(1/ √ λ).