2차원 아벨리안 모래더미 모델의 최적 지역 개입 전략

본 연구는 자기 조직화 임계성(self‑organized criticality, SOC)의 대표 모델인 2차원 아벨리안 모래더미(Abelian Sandpile Model, ASM)를 대상으로, 외부 제어자가 특정 정점에서 모래 알갱이를 제거함으로써 눈덩이 붕괴(avalanches)의 규모를 감소시키는 전략을 이론적으로 분석한다. 1. **모델 정의 및 기본 개념** - 격자 \(V=\{1,\dots ,L\}^2\) 에 대해 경계 정점을 싱크(sink)와 연결한 ‘와이어드’ 그래프를 사용한다. 각 정점 \(v\) 는 높이 \(\eta(v)\in\{0,1,2,3\}\) 을 갖고, \(3\) 인 경우를 ‘임계(critical)’라 부른다. - 모래 알갱이가 무작위 정점에 떨어지면, 해당 정점이 \(4\) 이상이 되면 ‘톱핑(toppling)’이 발생하고, 인접 정점으로 알갱이가 전달된다. 톱핑 연산 \(\tau_i\) 와 추가 연산 \(\alpha_i\) 는 교환법칙을 만족한다(Dhar, 1990). - ‘눈덩이 붕괴 크기’ \(x(\eta,v_i)\) 는 \(\alpha_i\) 후 안정화 과정에서 발생한 전체 톱핑 횟수이며, 이는 확률적 마코프 체인의 전이 확률에 직접 대응한다. 2. **기존 파동 기반 기대값 계산 방법의 한계** - Dorso와 Dadamia(2020)는 ‘파동(wave)’ 개념을 이용해, 임계 정점 집합 \(A\) (‘생성자’)에서 시작되는 붕괴의 첫 파동 \(W_A(\eta)\) 이 모든 시작점에 대해 동일하다고 가정하고, 이후 파동도 동일하게 취급해 기대 붕괴 크기를 재귀적으로 구했다. - 그러나 실제로는 첫 파동 후 남은 임계 정점 집합이 여러 연결 컴포넌트로 분리될 수 있다(Figure 1). 이 경우 두 번째 파동의 크기가 시작점에 따라 달라지며, 기존 식(20) 등은 적용되지 않는다. 3. **파동 분할을 포함한 일반화된 방법** - **Lemma 3.1**: 임계 컴포넌트 \(A\) 내 어느 정점에서 시작하든 첫 파동 \(W_{i1}(\eta)\) 은 동일함을 증명한다. 증명은 톱핑·추가 연산의 교환성, 그리고 임계 정점들의 연결성을 이용한 구성적 논증이다. - 첫 파동 후 남은 임계 정점 집합을 다시 연결 컴포넌트들로 분할하고, 각 컴포넌트에 대해 독립적인 두 번째 파동을 정의한다. 이렇게 하면 전체 파동 수 \(n(\eta,v_i)\)와 각 파동 크기 \(w_{ij}(\eta)\) 를 정확히 구할 수 있다. - 파동 분할 과정에서 발생할 수 있는 ‘음수 높이’ 문제를 방지하기 위해 ‘가상 톱핑 순서’를 도입한다. 이는 실제 물리적 톱핑과 동등함을 보이며, 알고리즘이 모든 생성자에 대해 기대 붕괴 크기를 반환함을 정리한다. 4. **정사각형 생성자에 대한 구체적 분석** - 정사각형 \(S_a\) (한 변의 길이 \(a\)) 가 비임계 정점으로 둘러싸여 있으면, 붕괴는 외부로 전파되지 않으므로 내부 구조만을 분석하면 된다. - **파동 구조**: 첫 파동은 정사각형 전체를 포함하고 크기는 \(a^2\). 두 번째 파동부터는 정사각형 내부의 ‘코너’와 ‘가장자리’ 정점이 차례로 탈락하면서 파동 크기가 감소한다. 파동 크기의 정확한 식은 \(w_{2}=4(a-1)\), \(w_{3}=4(a-2)\) 등으로 전개된다. - **개입 연산** \(\gamma_v\) (정점 \(v\) 의 모든 모래 제거) 의 효과를 정량화한다. 기대 붕괴 크기 감소량 \(\Delta E(v)=\mathbb{E}