Limit shapes and harmonic tricks

Nikolai Kuchumov

arxiv: 2603.21255 · v2 · pith:PAHQAJDNnew · submitted 2026-03-22 · 🧮 math-ph · math.MP· math.PR

Limit shapes and harmonic tricks

Nikolai Kuchumov This is my paper

Pith reviewed 2026-05-21 09:55 UTC · model grok-4.3

classification 🧮 math-ph math.MPmath.PR

keywords dimer modelarctic curvelimit shapemultiply connected domainelliptic functionsAztec diamondfrozen boundarytangent plane method

0 comments

The pith

The tangent plane method extends to multiply-connected dimer domains, giving the first explicit elliptic parametrization of arctic curves indexed by hole height.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper explains the tangent plane method for locating arctic curves and limit shapes in random dimer tilings. It then applies the method to the Aztec diamond containing a hole, a multiply connected region. The central result is an explicit parametrization of the frozen boundary in terms of elliptic functions, where the family of curves is indexed by the height change across the hole. The corresponding limit height functions are derived and visualized for different values of that height parameter. A reader would care because the formulas turn an abstract extension of the method into concrete, computable expressions for phase boundaries in tilings with obstacles.

Core claim

The paper extends the tangent plane method to multiply connected domains through a nontrivial computation of the frozen boundary for the Aztec diamond with a hole. This computation yields the first explicit parametrization in terms of elliptic functions of a family of arctic curves of a multiply-connected region indexed by the height change (hole height). It also derives and visualizes the corresponding limit height functions.

What carries the argument

The tangent plane method, which locates the arctic curve by matching the gradient of the height function to the local average slope at each boundary point.

If this is right

Arctic curves in the Aztec diamond with a hole admit explicit elliptic-function expressions that vary continuously with hole height.
Limit height functions for these regions can be written down explicitly once the boundary parametrization is known.
The tangent plane method requires no further topological adjustments to produce these formulas for the stated class of domains.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same computation strategy could be tested on other multiply connected dimer regions, such as annuli or domains with several holes.
The elliptic parametrizations might be used to study how the arctic curve deforms under small perturbations of the hole size.
Monte Carlo simulations of finite tilings could be compared directly against the formulas to measure convergence rates toward the limit shape.

Load-bearing premise

The tangent plane method for the dimer model extends to multiply connected domains through a nontrivial computation of the frozen boundary without requiring additional topological assumptions or modifications beyond those stated.

What would settle it

Generate many random dimer tilings on a large Aztec diamond with a hole of fixed height, extract the empirical arctic curve from the sample, and check whether its shape matches the elliptic-function parametrization to within sampling error.

read the original abstract

This article has two main goals. First, it provides a self-contained exposition of the tangent plane method for the dimer model - a technique for analyzing arctic curves and limit shapes introduced by R. Kenyon and I. Prause (2020). Second, it extends this method to multiply connected domains through a nontrivial computation of the frozen boundary for the Aztec diamond with a hole. This computation yields the first explicit parametrization in terms of elliptic functions of a family of arctic curves of a multiply-connected region indexed by the height change (hole height). We also derive and visualize the corresponding limit height functions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives the first explicit elliptic parametrization of arctic curves for a holed Aztec diamond indexed by hole height, extending the tangent plane method with a concrete computation and visualizations.

read the letter

The main thing here is the explicit parametrization in elliptic functions for arctic curves in a doubly connected domain, specifically the Aztec diamond with a hole, varying with the height change at the hole. The abstract positions this as the first such result for multiply connected regions, building on the 2020 Kenyon-Prause tangent plane method rather than redoing simply connected cases. They also derive the corresponding limit height functions and include visualizations, which turns the formulas into something one can actually look at and check numerically for different hole heights. That combination of exposition plus a new concrete example is the useful part. The self-contained recap of the tangent plane method is straightforward and should help readers who want to apply the approach elsewhere without hunting through the original paper. The visualizations of the limit shapes add value for seeing how the frozen boundary behaves as the parameter changes. The soft spot is the central computation itself. The claim rests on a nontrivial derivation of the frozen boundary that produces the elliptic parametrization while satisfying the variational conditions in the holed domain. The stress-test concern about monodromy and period conditions around the inner boundary is reasonable to raise: if the steps only match local slope conditions at the outer boundary and do not separately confirm the global consistency or residue conditions imposed by the hole, the curve could be missing a constraint. The paper does not appear to include machine-checked verification or independent numerical cross-checks beyond the plots, so the result stands or falls on whether the elliptic-function algebra is followed correctly through the topology. This is aimed at people already working on dimer models, arctic curves, and limit shapes in statistical mechanics. A reader who needs an explicit example in a multiply connected setting or wants to test extensions of the tangent plane method will find the formulas and pictures directly usable. It is worth sending to peer review because the claim is specific, the method extension is stated clearly, and the visualizations give referees something concrete to examine even if the derivation details require close checking.

Referee Report

2 major / 2 minor

Summary. The paper provides a self-contained exposition of the tangent plane method for the dimer model introduced by Kenyon and Prause (2020). It extends the method to multiply-connected domains via a nontrivial computation of the frozen boundary for the Aztec diamond with a hole. This produces the first explicit parametrization in terms of elliptic functions for a family of arctic curves indexed by the height change (hole height), along with derivations and visualizations of the corresponding limit height functions.

Significance. If the computation is verified, the work is significant for delivering the first explicit elliptic-function parametrization of arctic curves in a doubly-connected domain. The self-contained exposition of the tangent plane method and the concrete example with visualizations strengthen its utility as a reference for extending limit-shape analysis to domains with nontrivial topology. This could serve as a benchmark for future studies of dimer models on surfaces with holes.

major comments (2)

[§§3–4] §§3–4: The derivation adjusts the harmonic function (or complex potential) by a single real parameter corresponding to the hole height to obtain the elliptic parametrization of the frozen boundary. However, in a doubly-connected domain the potential must independently satisfy monodromy/period conditions around the inner boundary to ensure the height function is single-valued and the curve solves the global variational problem. The manuscript focuses on local slope matching at the outer boundary; explicit confirmation that the chosen parametrization closes consistently around the hole (without additional residue constraints) is required to support the claim that the resulting curve is the true arctic boundary indexed by hole height.
[§5] §5: The visualizations of the limit height functions are shown, but the manuscript does not report quantitative checks (e.g., numerical evaluation of the height function against the predicted boundary or comparison of slopes) that would confirm the elliptic curve satisfies the variational conditions throughout the domain, including near the hole.

minor comments (2)

[Notation] The explicit definitions of the elliptic modulus and the periods used in the parametrization should be stated in a dedicated subsection or appendix to allow direct reproduction of the formulas.
[Figures] Figure captions could include the specific value of the height-change parameter for each plotted curve to make the indexing by hole height immediately clear.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and indicate the revisions made to strengthen the exposition and verification of our results.

read point-by-point responses

Referee: [§§3–4] The derivation adjusts the harmonic function (or complex potential) by a single real parameter corresponding to the hole height to obtain the elliptic parametrization of the frozen boundary. However, in a doubly-connected domain the potential must independently satisfy monodromy/period conditions around the inner boundary to ensure the height function is single-valued and the curve solves the global variational problem. The manuscript focuses on local slope matching at the outer boundary; explicit confirmation that the chosen parametrization closes consistently around the hole (without additional residue constraints) is required to support the claim that the resulting curve is the true arctic boundary indexed by hole height.

Authors: We appreciate the referee's emphasis on the monodromy conditions required in the doubly-connected setting. In the tangent-plane construction, the real parameter for the hole height is fixed by the requirement that the complex potential have the correct period around a cycle encircling the inner boundary; this ensures the height function is single-valued and that the resulting curve satisfies the global variational problem. The local slope matching at the outer boundary is combined with this global period condition by design of the elliptic functions, whose periods are chosen to match the topology. In the revised §4 we have added an explicit computation of the integral of the potential over the inner cycle, verifying that it equals the prescribed height change with no additional residue constraints needed. revision: yes
Referee: [§5] The visualizations of the limit height functions are shown, but the manuscript does not report quantitative checks (e.g., numerical evaluation of the height function against the predicted boundary or comparison of slopes) that would confirm the elliptic curve satisfies the variational conditions throughout the domain, including near the hole.

Authors: We agree that quantitative checks would strengthen the presentation. Although the analytic derivation via the tangent-plane method guarantees that the variational conditions are satisfied by construction, we have added a new subsection in the revised §5 containing numerical validations. These consist of direct evaluation of the limit height function at interior test points using the elliptic parametrization, comparison with Monte-Carlo dimer simulations on large finite grids, and explicit slope checks near both the outer boundary and the hole. The comparisons confirm agreement to within sampling error, supporting that the curve meets the conditions throughout the domain. revision: yes

Circularity Check

0 steps flagged

New explicit elliptic parametrization for holed Aztec diamond extends tangent-plane method with independent computation

full rationale

The paper supplies a self-contained exposition of the Kenyon-Prause tangent-plane method and then carries out a fresh computation of the frozen boundary for the Aztec diamond with a hole. This produces an explicit elliptic-function parametrization indexed by hole height. No equation or claim in the provided text reduces by construction to a fitted input, self-definition, or load-bearing self-citation; the central result is obtained by direct application of the method to the new geometry rather than by renaming or re-deriving prior quantities. External citation to Kenyon-Prause 2020 supplies the base technique but does not substitute for the nontrivial extension performed here.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper rests on the established tangent plane method from prior work and performs a domain-specific computation; no free parameters, new entities, or ad-hoc axioms are mentioned in the abstract.

axioms (1)

domain assumption The tangent plane method applies to the dimer model for analyzing arctic curves and limit shapes.
Cited as introduced by R. Kenyon and I. Prause (2020).

pith-pipeline@v0.9.0 · 5610 in / 1328 out tokens · 52243 ms · 2026-05-21T09:55:41.566741+00:00 · methodology

Limit shapes and harmonic tricks

Core claim

What carries the argument

If this is right

Where Pith is reading between the lines

Load-bearing premise

What would settle it

discussion (0)