모래더미 모델에서 손상 전파와 임계 현상

초록

본 연구는 자기조직화 임계성(SOC) 모델인 모래더미(sandpile) 시스템에서 손상 전파(damage spreading)를 조사한다. 두 개의 동일한 초기 상태를 가진 시스템을 만든 뒤, 한쪽에 작은 국소 교란을 가하고 시간에 따라 두 시스템 사이의 차이(손상된 사이트 수)를 추적한다. 손상이 최초로 전역적으로 퍼지는 순간을 ‘임계 시간’이라 정의하고, 시스템 크기 L에 대한 임계 시간 τ(L)와 임계 시점의 손상된 사이트 총수 D(L)의 규모법칙을 측정한다. 수치 실험 결과 τ(L)∝L^α, D(L)∝L^β 형태의 전형적인 멱법칙을 보이며, α와 β는 각각 약 1.5와 2.0 정도의 값을 가진다. 이는 SOC 시스템에서도 전통적인 스핀 모델과 유사한 손상 전파 특성이 나타남을 시사한다.

상세 분석

본 논문은 SOC의 대표적 모델인 BTW(Bak‑Tang‑Wiesenfeld) 모래더미에 손상 전파 개념을 도입함으로써, 임계 상태에서의 동적 민감성을 정량화하고자 한다. 손상 전파는 일반적으로 두 복제 시스템을 동시 진화시키면서 초기 미세 교란이 어떻게 확대되는지를 측정하는 방법으로, 물리학에서는 이소스핀계, 퍼콜레이션, 네트워크 전파 등 다양한 분야에 적용돼 왔다. 저자들은 먼저 L×L 격자에 임계 높이(z_c=4)를 설정하고, 무작위로 입자를 떨어뜨려 평형 SOC 상태에 도달시킨다. 그 후, 동일한 초기 구성을 복제한 두 시스템 중 하나에 단일 사이트에 입자를 추가하는 국소 교란을 가한다. 이후 두 시스템을 동일한 랜덤 시드 하에 동시 업데이트함으로써, 외부 구동은 동일하게 유지하면서 내부 상태 차이만이 진화하도록 설계하였다.

손상은 각 시간 단계에서 두 격자 사이의 높이 차이가 존재하는 사이트 수로 정의되며, 시간에 따라 손상 영역이 확산하거나 수축한다. 저자들은 손상이 처음으로 전역적으로 퍼지는 순간, 즉 손상된 사이트가 시스템 전체에 걸쳐 연속적인 클러스터를 형성하는 시점을 ‘임계 시간 τ(L)’이라 명명한다. 이 τ(L)를 다양한 시스템 크기 L(예: L=16, 32, 64, 128, 256)에서 측정한 결과, 로그‑로그 플롯에서 직선 형태를 보이며 멱법칙 τ(L)∝L^α를 만족한다는 것을 확인하였다. α는 약 1.5±0.1로, 이는 1차원 확산(α=2)보다 느리지만, 2차원 전통적인 스핀계의 손상 전파 지수와 유사한 값을 가진다.

또한, 임계 시점에서 손상된 전체 사이트 수 D(L) 역시 L에 대한 멱법칙 D(L)∝L^β를 따랐으며, β는 약 2.0±0.05로 측정되었다. 이는 손상이 시스템 전체에 거의 균등하게 퍼졌음을 의미한다. 저자들은 이러한 결과가 SOC 시스템이 내부적으로 장거리 상관성을 갖고 있어, 작은 교란이 전역적인 변화를 일으킬 수 있음을 보여준다고 주장한다.

통계적 오류를 최소화하기 위해 각 L에 대해 수천 번의 독립 시뮬레이션을 수행했으며, 초기 교란 위치와 랜덤 입자 투입 순서가 결과에 미치는 영향을 분석하였다. 결과는 초기 교란 위치에 크게 의존하지 않으며, 평균적인 스케일링 지수는 시스템 전반에 걸쳐 일관되게 나타난다.

마지막으로, 저자들은 손상 전파 지수가 SOC 모델마다 다를 수 있음을 언급하며, 차원, 임계 규칙, 경계 조건 등에 따라 α와 β가 변할 가능성을 제시한다. 이는 향후 다양한 SOC 모델(예: 연속형 오일러-라그랑주 모델, 지진 모델 등)에서 손상 전파를 비교 연구함으로써, SOC의 보편적 동역학 특성을 규명하는 데 중요한 출발점이 될 수 있다.