Two-Periodic Aztec Diamond Dimer Model

• The two-periodic Aztec diamond is a weighted planar bipartite graph model defined by 2×2 periodic edge weights and analyzed through perfect t-embeddings and origami maps.
• Recursive constructions via gauge transformations and explicit update rules link the discrete wave equation’s basis solutions with local dimer density functions.
• In scaling limits, the model converges to space-like maximal surfaces in Minkowski space, revealing phase transitions with frozen regions and light-like cusp singularities.

The two-periodic Aztec diamond is the Aztec-diamond dimer model with $2\times 2$-periodic edge weights, studied as a weighted planar bipartite graph whose large-scale geometry is encoded by a perfect $t$-embedding and its associated origami map. In recent work, the coordinates of both the perfect $t$-embedding and the origami map were expressed as sums of density functions coming from the octahedron equation, while the broader theory of doubly periodic Aztec diamonds showed that the corresponding scaled $t$-surfaces converge to space-like maximal surfaces in Minkowski space, with frozen regions collapsing to boundary points and gas regions collapsing to interior light-like cusps (Berggren et al., 8 Aug 2025, Berggren et al., 6 Aug 2025, Chelkak et al., 2021).

1. Model and weighted dimer formulation

The Aztec diamond $A_{n+1}$ is defined as the set of faces $(j,k)\in\mathbb Z^2$ of the square grid $(\mathbb Z+\frac12)^2$ such that

$|j|+|k|\le n.$

In the two-periodic setting, the dimer weights are assigned by

$$\nu_e= \begin{cases} a, &n=1,2 \mod 4 \text{ and } j \text{ even} \quad \text{or}\quad n=3,0 \mod 4 \text{ and } j \text{ odd}, \ 1, & \text{ otherwise}, \end{cases}$$

where $a>0$, and $t$0 is adjacent to a face $t$1 with $t$2 odd. This is the weighted Aztec-diamond sequence used in the two-periodic theory (Berggren et al., 8 Aug 2025).

A key technical step is to replace the original graph by a reduced Aztec diamond $t$3, obtained by gauge transformations, contraction of degree-$t$4 vertices, and merging of parallel edges. This reduction is essential because the shuffling and spider-move dynamics become compatible with the $t$5-embedding formalism. The reduction also makes it possible to transfer local graph transformations into explicit recurrences for the embedding and the origami map (Berggren et al., 8 Aug 2025).

Within the later theory of doubly periodic Aztec diamonds, the two-periodic Aztec diamond appears as the $t$6-periodic specialization. In that setting, it serves as the genus-one example in which the same spectral curve and the same limit shape can nevertheless produce different limiting $t$7-surfaces, depending on how the periodic weights are shifted (Berggren et al., 6 Aug 2025).

2. Perfect $t$8-embeddings and origami maps

Given a weighted planar bipartite graph $t$9, a $t$0-embedding $t$1 is a proper embedding of the augmented dual graph into $t$2 such that two conditions hold: the sum of the angles at each inner vertex at corners corresponding to black faces is $t$3 and similarly for white faces, and the geometric dual-edge lengths are gauge equivalent to the original dimer weights. A perfect $t$4-embedding is a $t$5-embedding whose outer face is a tangential polygon to a circle and whose non-boundary edges emanating from boundary vertices are bisectors of the corresponding boundary angles (Berggren et al., 8 Aug 2025).

In the general theory, this “perfectness” is a boundary regularity condition added to the ordinary finite $t$6-embedding axioms. The angle condition is the discrete integrability condition behind the construction, and the geometric edge-length condition identifies the embedding with a complex Kasteleyn realization of the dimer model. The associated domain can be written in terms of the boundary data $t$7 as

$t$8

and the boundary vertices satisfy

$t$9

The terminology “perfect” therefore refers to a specific finite-volume boundary condition, not to any extremal or optimality property (Chelkak et al., 2021).

Every $t$0-embedding carries an origami map. Choosing an origami square root function $t$1, one defines

$t$2

Equivalently, the origami map is obtained by folding the plane along the edges of the $t$3-embedding. In the perfect case, the pair $t$4 is the basic discrete geometric object underlying the asymptotic theory of height fluctuations and Minkowski-space limits (Berggren et al., 8 Aug 2025).

3. Recursive construction on the reduced two-periodic Aztec diamond

For the reduced two-periodic Aztec diamond, the perfect $t$5-embedding is constructed recursively. If $t$6 denotes the position in $t$7 of the dual vertex indexed by $t$8 in the perfect $t$9-embedding of $A_{n+1}$0, then $A_{n+1}$1 is obtained from $A_{n+1}$2 by explicit update rules derived from the shuffling algorithm and preservation under elementary graph moves. The outer boundary is embedded in a rhombus with diagonals $A_{n+1}$3 and $A_{n+1}$4, with boundary conditions

$A_{n+1}$5

The origami map $A_{n+1}$6 satisfies the same update rules, but with boundary conditions

$A_{n+1}$7

This gives a local linear evolution for both $A_{n+1}$8 and $A_{n+1}$9 (Berggren et al., 8 Aug 2025).

The recurrence is reformulated on

$(j,k)\in\mathbb Z^2$0

For boundary parameters $(j,k)\in\mathbb Z^2$1, one introduces basis solutions $(j,k)\in\mathbb Z^2$2 of the universal discrete wave equation, so that

$(j,k)\in\mathbb Z^2$3

The perfect $(j,k)\in\mathbb Z^2$4-embedding and origami map are then expressed by

$(j,k)\in\mathbb Z^2$5

and

$(j,k)\in\mathbb Z^2$6

These formulas clarify the structural relation between the two objects. The $(j,k)\in\mathbb Z^2$7-embedding is the antisymmetric combination of the directional components in the east-west and north-south directions, whereas the origami map is the corresponding symmetric combination. In the two-periodic case, this decomposition replaces the more elementary single-fundamental-solution description available in the uniform model (Berggren et al., 8 Aug 2025).

4. Octahedron equation and density-function representation

The central integrable input is the octahedron equation

$(j,k)\in\mathbb Z^2$8

with initial data specified on the parity classes $(j,k)\in\mathbb Z^2$9. The associated density functions are defined by

$(\mathbb Z+\frac12)^2$0

and admit the probabilistic interpretation

$(\mathbb Z+\frac12)^2$1

where $(\mathbb Z+\frac12)^2$2 is the number of dimers around face $(\mathbb Z+\frac12)^2$3 in a random dimer configuration (Berggren et al., 8 Aug 2025).

For a special $(\mathbb Z+\frac12)^2$4-periodic choice of octahedron initial data,

$(\mathbb Z+\frac12)^2$5

under the constraints

$(\mathbb Z+\frac12)^2$6

the linear recurrence satisfied by the densities matches the recurrence for the $(\mathbb Z+\frac12)^2$7-embedding basis solutions. The particular choice

$(\mathbb Z+\frac12)^2$8

produces the two-periodic Aztec-diamond weights on the octahedron side (Berggren et al., 8 Aug 2025).

The decisive identification is that the shifted fundamental solutions of the discrete wave equation coincide with the density functions: $(\mathbb Z+\frac12)^2$9 Consequently, the directional components $|j|+|k|\le n.$0, and therefore the coordinates of both $|j|+|k|\le n.$1 and $|j|+|k|\le n.$2, are explicit sums of density functions. The main theorem states that these formulas define a perfect $|j|+|k|\le n.$3-embedding and its origami map of the two-periodic Aztec diamond. In this form, the real and imaginary coordinates of the embedding become linear combinations of expected local dimer deficiencies, giving the geometry a direct probabilistic interpretation (Berggren et al., 8 Aug 2025).

5. Scaling limits, phases, and maximal surfaces

For periodic Aztec diamonds, the discrete pair

$|j|+|k|\le n.$4

is viewed as a $|j|+|k|\le n.$5-surface. The large-scale theory proves that these $|j|+|k|\le n.$6-surfaces converge to space-like maximal surfaces in Minkowski space. In the no-gas case the limit lies in $|j|+|k|\le n.$7, whereas in the general gas case it naturally lives in $|j|+|k|\le n.$8. The global conformal structure remains robust and coincides with the Kenyon–Okounkov conformal structure, but the detailed embedding geometry depends on the precise placement of the periodic weights (Berggren et al., 6 Aug 2025).

The geometric effect of the phase structure is highly specific. All frozen regions collapse to four boundary points, regardless of how many frozen regions are present, while each gas region collapses to a distinct interior light-like cusp. The limiting parametrization satisfies the harmonicity and conformality identities appropriate to maximal surfaces, and the spacelike inequality holds in the interior with equality on the boundary and at cusp degeneration. In the gas regime, the cusps are therefore genuine geometric singularities of the maximal surface rather than removable artifacts (Berggren et al., 6 Aug 2025).

The two-periodic Aztec diamond is the $|j|+|k|\le n.$9-periodic, genus-one case of this theory. It is singled out in the doubly periodic analysis because there are four shifted versions with the same spectral curve and the same limit shape but different limiting $$\nu_e= \begin{cases} a, &n=1,2 \mod 4 \text{ and } j \text{ even} \quad \text{or}\quad n=3,0 \mod 4 \text{ and } j \text{ odd}, \ 1, & \text{ otherwise}, \end{cases}$$0-surfaces. For the four weight choices

$$\nu_e= \begin{cases} a, &n=1,2 \mod 4 \text{ and } j \text{ even} \quad \text{or}\quad n=3,0 \mod 4 \text{ and } j \text{ odd}, \ 1, & \text{ otherwise}, \end{cases}$$1

the limiting surfaces satisfy

$$\nu_e= \begin{cases} a, &n=1,2 \mod 4 \text{ and } j \text{ even} \quad \text{or}\quad n=3,0 \mod 4 \text{ and } j \text{ odd}, \ 1, & \text{ otherwise}, \end{cases}$$2

but for $$\nu_e= \begin{cases} a, &n=1,2 \mod 4 \text{ and } j \text{ even} \quad \text{or}\quad n=3,0 \mod 4 \text{ and } j \text{ odd}, \ 1, & \text{ otherwise}, \end{cases}$$3 they cannot be embedded into $$\nu_e= \begin{cases} a, &n=1,2 \mod 4 \text{ and } j \text{ even} \quad \text{or}\quad n=3,0 \mod 4 \text{ and } j \text{ odd}, \ 1, & \text{ otherwise}, \end{cases}$$4 by global shifts or rotations of the origami map, and the cusps are not locally contained in any lower-dimensional subspace of $$\nu_e= \begin{cases} a, &n=1,2 \mod 4 \text{ and } j \text{ even} \quad \text{or}\quad n=3,0 \mod 4 \text{ and } j \text{ odd}, \ 1, & \text{ otherwise}, \end{cases}$$5. A common misconception is therefore that periodicity or the spectral curve alone determines the limiting $$\nu_e= \begin{cases} a, &n=1,2 \mod 4 \text{ and } j \text{ even} \quad \text{or}\quad n=3,0 \mod 4 \text{ and } j \text{ odd}, \ 1, & \text{ otherwise}, \end{cases}$$6-surface; the published results show instead that the cusp positions and even the ambient Minkowski subspace can depend on the specific shift of the periodic weight pattern (Berggren et al., 6 Aug 2025).

6. Fluctuation theory, conformal structure, and open directions

The motivation for perfect $$\nu_e= \begin{cases} a, &n=1,2 \mod 4 \text{ and } j \text{ even} \quad \text{or}\quad n=3,0 \mod 4 \text{ and } j \text{ odd}, \ 1, & \text{ otherwise}, \end{cases}$$7-embeddings comes from the general dimer scaling theory. The foundational result states that if a sequence of perfect $$\nu_e= \begin{cases} a, &n=1,2 \mod 4 \text{ and } j \text{ even} \quad \text{or}\quad n=3,0 \mod 4 \text{ and } j \text{ odd}, \ 1, & \text{ otherwise}, \end{cases}$$8-embeddings converges, together with its origami maps, to a Lorentz-minimal surface and satisfies the technical assumptions $$\nu_e= \begin{cases} a, &n=1,2 \mod 4 \text{ and } j \text{ even} \quad \text{or}\quad n=3,0 \mod 4 \text{ and } j \text{ odd}, \ 1, & \text{ otherwise}, \end{cases}$$9 and $a>0$0, then the gradients of the dimer height correlation functions converge to those of the Dirichlet Gaussian free field in the intrinsic metric of that surface. The theorem is explicitly a gradient statement rather than full field convergence, and the authors note that their proof does not exclude convergence of the height field to “GFF plus an independent global random variable” (Chelkak et al., 2021).

Uniform Aztec diamonds provided the first verified example of a model satisfying all conditions of this general theorem: perfect $a>0$1-embeddings of uniformly weighted Aztec diamonds were used to prove convergence of gradients of height fluctuations to those of the Gaussian free field, thereby confirming a prediction from earlier work (Berggren et al., 2023). The two-periodic theory does not merely parallel that example. By expressing the perfect $a>0$2-embedding and origami map through octahedron-equation densities, it provides an explicit integrable realization of the geometry in a genuinely periodic setting and suggests a route to asymptotic analysis even when inverse-Kasteleyn formulas are unavailable (Berggren et al., 8 Aug 2025).

The doubly periodic theory adds a further conjectural layer. For models with gas regions, global height fluctuations depend on the Kenyon–Okounkov conformal structure and on a shift parameter $a>0$3; in the present setting, $a>0$4, the spectral-divisor parameter. Since the cusp locations of the limiting $a>0$5-surface depend nontrivially on $a>0$6, the conjecture is that the cusp apices encode this discrete-Gaussian shift. This suggests that the limiting $a>0$7-surface may carry not only the conformal geometry of the liquid region but also the extra global data needed in the gas regime (Berggren et al., 6 Aug 2025).

At the level of the general theory, several questions remain open: existence of perfect $a>0$8-embeddings for sufficiently nondegenerate weighted planar bipartite graphs, uniqueness up to the natural Lorentz symmetries, formulation of a discrete notion of Lorentz-minimality, and a fuller description of how frozen, liquid, and gaseous regions are represented by the geometry of perfect $a>0$9-embeddings and their limiting surfaces (Chelkak et al., 2021). In that broader program, the two-periodic Aztec diamond occupies a central position because it is simultaneously an explicit integrable model, a test case for periodic $t$00-surface geometry, and a prototype for the interaction between density observables, conformal structure, and cusp singularities.