확률적 셀룰러 오토마톤 p182‑q200의 멸종‑생존 전이와 지향 퍼콜레이션

초록

본 연구는 Wolfram 규칙 182와 200을 확률적으로 혼합한 셀룰러 오토마톤(p182‑q200)을 대상으로 평균장 이론과 대규모 몬테카를로 시뮬레이션을 수행하였다. 파라미터를 조절해 한 규칙을 끄고 다른 규칙을 켤 때, 시스템은 멸종(활성도 0)에서 지속적인 활성 상태로 전이하는 임계 현상을 보이며, 그 임계 지수는 지향 퍼콜레이션(Directed Percolation) 보편군에 속한다. 특히 규칙 200에 가까워질수록 정지 상태의 밀도 프로파일이 특이하게 형성되고, 확산적(다이프루시브) 동역학이 지배한다. 논문 말미에서는 이 모델과 연관된 확률적 격자 가스 모델에 대한 흥미로운 연결 고리도 제시한다.

상세 분석

본 논문은 두 개의 1차원 이진 셀룰러 오토마톤(CA) 규칙, 즉 Wolfram 번호 182와 200을 확률적으로 섞어 만든 새로운 확률적 셀룰러 오토마톤(p182‑q200)을 정의한다. 규칙 182는 주변 3칸의 상태에 따라 중앙칸을 1로 만드는 비교적 활성화된 전이 규칙이며, 규칙 200은 주변이 모두 0일 때만 중앙을 1로 만드는 억제적인 규칙이다. 두 규칙을 혼합하는 파라미터 p(=1‑q)를 도입함으로써, p=1일 때는 순수 182, p=0일 때는 순수 200이 된다. 저자는 먼저 1‑site, 2‑site, 3‑site 평균장 근사를 적용해 정적 밀도 ρ의 자기일관 방정식을 도출하고, 이론적으로 임계점 p_c≈0.48을 예측한다. 그러나 평균장 근사는 차원 의존적 상관을 무시하므로 정확한 임계값을 제공하지 못한다는 한계가 있다.

이를 보완하기 위해 저자는 L=10^4~10^5 크기의 1차원 격자에 대해 비동시식 업데이트 방식을 사용한 몬테카를로 시뮬레이션을 수행한다. 초기 조건은 무작위로 ρ_0=0.5인 상태를 사용하고, 충분한 열화 시간(10^5 스텝) 후에 정적 밀도와 시간 의존적 활성도 A(t)를 측정한다. 결과는 p가 약 0.475±0.005에서 급격히 변하는 멸종‑생존 전이를 보이며, 이는 전통적인 지향 퍼콜레이션(DP) 전이와 동일한 특성을 가진다. 특히, 임계점 근처에서 A(t)∼t^{−δ} (δ≈0.159)와 ρ∼(p−p_c)^{β} (β≈0.276)라는 지수값이 DP 보편군의 알려진 값과 일치한다는 점이 강조된다.

또한, p가 0에 가까워질수록(규칙 200에 가까워질수록) 시스템은 매우 느린 확산적 동역학을 보인다. 이는 활성 영역이 드물게 생성되고, 생성된 영역이 주변에 퍼지는 과정이 랜덤 워크와 유사하게 진행되기 때문이다. 저자는 이 현상을 “diffusive bottleneck”이라 부르며, 정적 밀도 프로파일이 중앙에 높은 밀도를 유지하고 양쪽으로 급격히 감소하는 형태를 띤다고 보고한다. 이러한 프로파일은 규칙 200이 자체적으로 “absorbing” 상태(모든 셀 0)로 수렴하는 경향이 강함을 반영한다.

마지막으로, 논문은 p182‑q200과 구조적으로 유사한 확률적 격자 가스 모델을 제시한다. 이 모델은 입자들이 인접한 사이트에 이동하거나 소멸·생성되는 과정을 규정하며, p와 q가 각각 입자 생성·소멸 확률에 대응한다. 저자는 두 모델이 동일한 마코프 연산자를 공유하므로, p182‑q200의 임계 현상이 격자 가스에서도 동일한 universality class를 따른다고 추론한다. 이는 확률적 CA와 격자 가스 사이의 이론적 연결 고리를 제공하며, 향후 비평형 통계물리학에서 중요한 연구 방향을 제시한다.

전반적으로, 본 연구는 간단한 규칙 혼합을 통해 복잡한 비평형 상전이를 구현하고, 평균장 이론과 수치 시뮬레이션을 조합해 정확한 임계 특성을 규명한 점이 의의이다. 특히, DP 보편군에 속한다는 결론은 다양한 물리·생물 시스템에서 관찰되는 멸종‑생존 전이를 이해하는 데 유용한 모델 사례를 제공한다.