Local statistics for random domino tilings of the Aztec diamond

We prove an asymtotic formula for the probability that, if one chooses a domino tiling of a large Aztec diamond at random according to the uniform distribution on such tilings, the tiling will contain a domino covering a given pair of adjacent lattice squares. This formula quantifies the effect of the diamond's boundary conditions on the behavior of typical tilings; in addition, it yields a new proof of the arctic circle theorem of Jockusch, Propp, and Shor. Our approach is to use the saddle point method to estimate certain weighted sums of squares of Krawtchouk polynomials (whos relevance to domino tilings is demonstrated elsewhere), and to combine these estimates with some exponential sum bounds to deduce our final result. This approach generalizes straightforwardly to the case in which the probability distribution on the set of tilings incorporates bias favoring horizontal over vertical tiles or vice versa. We also prove a fairly general large deviation estimate for domino tilings of simply connected planar regions that implies that some of our results on Aztec diamonds apply to many other similar regions as well.