Suppose that we have a tile $T$ in say $\mathbb{Z}^2$, meaning a finite subset of $\mathbb{Z}^2$. It may or may not be the case that $T$ tiles $\mathbb{Z}^2$, in the sense that $\mathbb{Z}^2$ can be partitioned into copies of $T$. But is there always some higher dimension $\mathbb{Z}^d$ that can be tiled with copies of $T$? We prove that this is the case: for any tile in $\mathbb{Z}^2$ (or in $\mathbb{Z}^n$, any $n$) there is a $d$ such that $\mathbb{Z}^d$ can be tiled with copies of it. This proves a conjecture of Chalcraft.

The first FFT-based algorithm for numerical homogenization from high-resolution images was proposed by Moulinec and Suquet in 1994 as an alternative to finite elements and twenty years later, it is still widely used in computational micromechanics of materials. The method is based on an iterative solution to an integral equation of the Lippmann-Schwinger type, whose kernel can be explicitly expressed in the Fourier domain. Only recently, it has been recognized that the algorithm has a variational structure arising from a Fourier-Galerkin method. In this talk, I will show how this insight can be used to significantly improve the performance of the original Moulinec-Suquet solver. In particular, I will focus on (i) influence of an iterative solver used to solve the system of linear algebraic equations, (ii) effects of numerical integration of the Galerkin weak form, and (iii) convergence of an a-posteriori bound on the solution during iterations.

Hardy and Littlewood's approximate functional equation for quadratic Weyl sums (theta sums) provides, by iterative application, a powerful tool for the asymptotic analysis of such sums. The classical Jacobi theta function, on the other hand, satisfies an exact functional equation, and extends to an automorphic function on the Jacobi group. In the present study we construct a related, almost everywhere non-differentiable automorphic function, which approximates quadratic Weyl sums up to an error of order one, uniformly in the summation range. This not only implies the approximate functional equation, but allows us to replace Hardy and Littlewood's renormalization approach by the dynamics of a certain homogeneous flow. The great advantage of this construction is that the approximation is global, i.e., there is no need to keep track of the error terms accumulating in an iterative procedure. Our main application is a new functional limit theorem, or invariance principle, for theta sums. The interesting observation here is that the paths of the limiting process share a number of key features with Brownian motion (scale invariance, invariance under time inversion, non-differentiability), although time increments are not independent and the value distribution at each fixed time is distinctly different from a normal distribution. Joint work with Francesco Cellarosi.