Domino tiling

Thurston's height condition

William Thurston (1990) describes a test for determining whether a simply-connected region, formed as the union of unit squares in the plane, has a domino tiling. He forms an undirected graph that has as its vertices the points (x,y,z) in the three-dimensional integer lattice, where each such point is connected to four neighbors: if x + y is even, then (x,y,z) is connected to (x + 1,y,z + 1), (x − 1,y,z + 1), (x,y + 1,z − 1), and (x,y − 1,z − 1), while if x + y is odd, then (x,y,z) is connected to (x + 1,y,z − 1), (x − 1,y,z − 1), (x,y + 1,z + 1), and (x,y − 1,z + 1). The boundary of the region, viewed as a sequence of integer points in the (x,y) plane, lifts uniquely (once a starting height is chosen) to a path in this three-dimensional graph. A necessary condition for this region to be tileable is that this path must close up to form a simple closed curve in three dimensions, however, this condition is not sufficient. Using more careful analysis of the boundary path, Thurston gave a criterion for tileability of a region that was sufficient as well as necessary.

Counting tilings of regions

The number of ways to cover an ${\displaystyle m\times n}$ rectangle with ${\displaystyle \frac{mn}{2}}$ dominoes, calculated independently by Temperley & Fisher (1961) and Kasteleyn (1961), is given by $${\displaystyle \prod _{j=1}^{\lceil \frac{m}{2}\rceil }\prod _{k=1}^{\lceil \frac{n}{2}\rceil }\left(4{\cos }^2\frac{\pi j}{m+1}+4{\cos }^2\frac{\pi k}{n+1}\right)}$$ (sequence A099390 in the OEIS).

When both m and n are odd, the formula correctly reduces to zero possible domino tilings.

A special case occurs when tiling the ${\displaystyle 2\times n}$ rectangle with n dominoes: the sequence reduces to the Fibonacci sequence.

Another special case happens for squares with m = n = 0, 2, 4, 6, 8, 10, 12, ... is

1, 2, 36, 6728, 12988816, 258584046368, 53060477521960000, ... (sequence A004003 in the OEIS).

These numbers can be found by writing them as the Pfaffian of an ${\displaystyle mn\times mn}$ skew-symmetric matrix whose eigenvalues can be found explicitly. This technique may be applied in many mathematics-related subjects, for example, in the classical, 2-dimensional computation of the dimer-dimer correlator function in statistical mechanics.

The number of tilings of a region is very sensitive to boundary conditions, and can change dramatically with apparently insignificant changes in the shape of the region. This is illustrated by the number of tilings of an Aztec diamond of order n, where the number of tilings is 2^{(n + 1)n/2}. If this is replaced by the "augmented Aztec diamond" of order n with 3 long rows in the middle rather than 2, the number of tilings drops to the much smaller number D(n,n), a Delannoy number, which has only exponential rather than super-exponential growth in n. For the "reduced Aztec diamond" of order n with only one long middle row, there is only one tiling.

Tatami

Tatami are Japanese floor mats in the shape of a domino (1x2 rectangle). They are used to tile rooms, but with additional rules about how they may be placed. In particular, typically, junctions where three tatami meet are considered auspicious, while junctions where four meet are inauspicious, so a proper tatami tiling is one where only three tatami meet at any corner. The problem of tiling an irregular room by tatami that meet three to a corner is NP-complete.

Applications in statistical physics

There is a one-to-one correspondence between a periodic domino tiling and a ground state configuration of the fully frustrated Ising model on a two-dimensional periodic lattice. At the ground state, each plaquette of the spin model must contain exactly one frustrated interaction. Therefore, viewing from the dual lattice, each frustrated edge must be "covered" by a 1x2 rectangle, such that the rectangles span the entire lattice and do not overlap, or a domino tiling of the dual lattice.