아지텍 다이아몬드 구멍의 아틱 곡선과 극한 형태 분석

본 논문은 두 가지 주요 목표를 가지고 있다. 첫 번째는 도미노 모델의 아틱 곡선과 극한 형태를 분석하기 위해 케니언·프라우즈(2020)가 제시한 ‘접평면 방법(tangent plane method)’을 완전하고 자가 포함적인 형태로 정리·설명하는 것이다. 두 번째는 이 방법을 다중 연결 영역, 특히 중앙에 구멍이 있는 아지텍 다이아몬드(‘Aztec diamond with a hole’)에 적용하여, 구멍의 크기(또는 ‘hole height’)를 매개변수로 하는 아틱 곡선 군을 엘립틱 함수로 명시적으로 파라미터화하고, 대응하는 극한 높이 함수(limit height function)를 시각화한다. 1. **배경 및 변분 원리** 도미노 타일링은 이중 격자 위의 완전 매칭으로 모델링되며, 각 타일링은 높이 함수 H : V(Γ)→ℤ 으로 코딩된다. 대규모 한계(N→∞)에서 재스케일된 높이 H/N 은 확률적으로 거의 확정적인 연속 함수 h 로 수렴한다. 이 h 는 ‘표면 장력’ σ(∇h) 의 적분을 최소화하는 변분 문제의 해이며, σ는 뉴턴 다각형 N 위에서 볼록하고, 그 헤시안은 양정정칙 메트릭을 만든다. 이 변분 원리는 ‘동결 영역(frozen region)’과 ‘액체 영역(liquid region)’을 구분하는 아틱 곡선(‘frozen boundary’)을 정의한다. 2. **복소 기울기와 스펙트럴 커브** 케니언·오쿠노프(Kenyon‑Okounkov)와 코언·케니언·프라우즈(Cohn‑Kenyon‑Propp) 등은 복소 기울기 z 을 도입해 복소 버거스 방정식과 스펙트럴 커브 P(z,w)=0 을 연결하였다. 여기서 z,w 는 스펙트럴 커브 위의 좌표이며, Harnack 커브일 경우 실축 위에 두 개의 복소 공액점이 존재한다. 복소 기울기와 실제 기울기 (s,t) 는 s= (1/π)·arg w, t=−(1/π)·arg z 로 연결된다. 3. **접평면 방법의 핵심 아이디어** 접평면 방법은 h 의 그래프를 각 점 (x₀,y₀) 에서의 접평면 P_{x₀,y₀}: s(x₀,y₀) x + t(x₀,y₀) y + c(x₀,y₀)=z 으로 표현한다. 여기서 s=∂h/∂x, t=∂h/∂y, c=h−xs−yt. 중요한 사실은 s,t,c 가 복소 매개변수 u (‘intrinsic coordinate’) 위에서 조화함수라는 점이다. 이 조화성은 스펙트럴 커브가 Harnack 곡선일 때, Ronkin 함수의 라그랑주 방정식과 Harnack‑정리의 결합을 통해 증명된다. 4. **내재 좌표와 매핑** 액체 영역 L 을 복소 좌표 u (보통 단일 연결이면 반평면 H, 다중 연결이면 원환 Σ⁺)에 균일화한다. 이때 z(u), w(u) 는 u 에 대한 엘립틱(또는 초월) 함수이며, s(u), t(u), c(u) 는 각각 Re 또는 Im 연산을 통해 얻어진 조화함수이다. 매핑 ∇h: L→N 을 복소 좌표 z 와 연결하면, ∇h = (1/π)(arg w, −arg z) 가 된다. 5. **아틱 곡선의 파라미터화** 경계 ∂L 위에서 복소 접선 방정식     s(u) x + t(u) y + c(u)=0 을 실수화하면 두 실수 자유도 x(u), y(u) 가 결정된다. 따라서 u∈∂Σ 에 대해 (x(u),y(u)) 가 바로 아틱 곡선 ∂L 의 파라미터화가 된다. 내부 u∈Σ 에 대해서는 동일한 식을 이용해 h(u)=s(u)x(u)+t(u)y(u)+c(u) 를 계산하면, 극한 높이 함수가 얻어진다. 6. **다중 연결 영역: 구멍이 있는 아지텍 다이아몬드** 본 논문의 핵심 응용은 중앙에 작은 아지텍 다이아몬드(크기 ⌊κ(τ)N⌋)가 삽입된 큰 아지텍 다이아몬드(차수 N)이다. 이 경우 액체 영역은 두 개의 연결 성분(외부와 내부)으로 나뉘며, 복소 매개변수 u 는 모듈러 파라미터 τ 에 의해 정의된 원환 Σ⁺ ≃