동적 몬테카를로를 이용한 고정밀 강성 퍼콜레이션 연구

초록

삼각 격자에서 결합을 순차적으로 추가하며 강성 클러스터의 형성과 병합을 실시간으로 추적하는 동적 pebble game 알고리즘을 개발하였다. 시간적 자기유사성, 대규모 연쇄 병합 현상을 발견하고, 이벤트 기반 앙상블을 이용해 임계점 (p_c = 0.6602778(10)), 역상관길이 지수 (1/\nu = 0.850(3)), 프랙털 차원 (d_f = 1.850(2))를 기존보다 훨씬 높은 정밀도로 측정하였다.

상세 분석

본 논문은 강성 퍼콜레이션(RP)의 동적 특성을 정량적으로 규명하기 위해, 기존 정적 pebble game(PG) 알고리즘을 실시간으로 적용할 수 있는 동적 PG 알고리즘을 설계하였다. 핵심은 두 가지 혁신이다. 첫째, 각 결합 추가 시 클러스터 소속 정보를 양방향으로 유지하는 영구 데이터 구조를 도입해 결합 독립성 판단을 O(1) 시간에 수행한다. 이는 매 단계마다 전체 PG를 재실행하는 전통적 접근보다 크게 효율적이다. 둘째, “가장 큰 클러스터를 제외한다(let go of the largest)”는 최적화를 통해 병합 과정에서 가장 큰 클러스터를 탐색 대상에서 배제함으로써, 임계점 부근과 초임계 영역에서 발생하는 대규모 클러스터 탐색 비용을 현저히 감소시켰다. 결과적으로 전체 시뮬레이션 복잡도는 O(V^1.15) 수준으로, 정적 PG와 동등한 효율을 달성하였다.

동적 시뮬레이션에서는 결합 추가를 시간 변수 (t = T/E) 로 정의하고, 각 단계에서 강성 클러스터의 크기 변화(갭)와 병합된 클러스터 수를 기록한다. 분석 결과, 이러한 양들은 모두 파워‑law 형태의 분포를 보이며, 시간적 자기유사성을 나타낸다. 특히, 단일 결합이 트리거하는 대규모 연쇄 병합(cascade) 현상이 관찰되었으며, 최대 병합 클러스터 수 (M)는 시스템 크기 (L)에 대해 (M \sim L^{1/\nu}) 로 스케일링한다. 이는 역상관길이 지수와 동일한 지수값을 갖는다는 점에서, 강성 퍼콜레이션의 비국소적 제약이 강한 유한크기 보정의 근본 원인임을 시사한다.

이러한 동적 이벤트를 기반으로 “이벤트 기반 앙상블(event‑based ensemble)”을 구성하였다. 각 시뮬레이션에서 특정 이벤트(예: 최초의 시스템‑스팬닝 강성 클러스터 형성) 발생 시점 (t_L)을 pseudo‑critical point 로 정의하고, (t_L - t_c \sim L^{-1/\nu}) 와 (\sigma_{t_L} \sim L^{-1/\nu}) 를 검증하였다. 또한, 최대 클러스터 크기 (C_1) 와 그 변동성을 동일한 스케일링 형태로 분석함으로써, 기존 정적 앙상블에서 나타났던 비표준 FSS 현상을 제거하고, 순수한 임계 거동을 고정밀도로 추정할 수 있었다.

결과적으로, 본 연구는 (1) 동적 PG 알고리즘을 통한 대규모 2D RP 시뮬레이션 가능성, (2) 시간적 자기유사성과 연쇄 병합 현상의 정량적 규명, (3) 이벤트 기반 앙상블을 활용한 임계점 및 임계 지수의 고정밀 추정이라는 세 가지 주요 성과를 제공한다. 이는 강성 퍼콜레이션이 베르누이 퍼콜레이션과 다른 보편적 클래스에 속함을 재확인하고, 비국소적 제약이 유한크기 보정에 미치는 영향을 명확히 하는 데 기여한다.