2차원 사이트‑퍼콜레이션에서 격자 가스 자동자를 이용한 확산 전이 현상의 구성 의존성

초록

**
본 연구는 2차원 정사각형 사이트‑퍼콜레이션 격자 위에서 격자 가스 자동자(LGA) 모델을 활용해 무작위 보행을 시뮬레이션하였다. 점유율이 임계값을 초과할 때 초기의 서브‑디퓨전에서 장기적인 정상 확산으로 전이하는 이론적 예측을 검증했으며, 전이는 개별 격자 구성에 따라 나타났지만, 다수 구성의 평균을 취하면 전이가 사라지는 현상을 확인하였다. 이는 평균화 과정에서 짧은 거리의 불균일성이 소멸되기 때문으로 해석된다.

**

상세 요약

**
퍼콜레이션 이론에 따르면, 점유된 사이트 비율(p)이 임계점(p_c)보다 클 경우 무한 클러스터가 형성되어 장거리 전도가 가능해진다. 이때 무작위 보행자는 초기 시간 구간에서 클러스터 내부의 복잡한 구조에 갇히면서 평균제곱변위 ⟨r²(t)⟩∝t^{2ν} (ν<½) 형태의 서브‑디퓨전(비정상 확산) 현상을 보인다. 시간이 충분히 흐르면 보행자는 클러스터의 전역적인 연결성을 이용해 정상 확산 ⟨r²(t)⟩∝t 로 전이한다는 것이 기존의 스케일링 이론이다.

연구진은 격자 가스 자동자(LGA)를 선택하였다. LGA는 입자 충돌 규칙과 이동 규칙을 통해 미시적 흐름을 구현하면서도 계산 효율성이 높아 대규모 2D 퍼콜레이션 격자(10⁴10⁵ 셀)에서 수천 번의 보행을 동시에 수행할 수 있다. 점유율을 p_c≈0.5927보다 약간 높은 0.600.80 구간으로 설정하고, 각 p에 대해 500개의 서로 다른 격자 구성을 생성하였다. 각 구성마다 10⁴ 단계의 보행을 수행하고, 평균제곱변위와 확산계수 D(p)∝(p−p_c)^{μ}를 추정하였다.

결과는 두드러진 두 가지 현상을 보여준다. 첫째, D(p)의 임계 지수 μ≈1.30±0.05 로, 기존 이론값(μ≈1.30)과 일치한다. 이는 LGA가 퍼콜레이션 네트워크의 전도 특성을 정확히 포착함을 의미한다. 둘째, ⟨r²(t)⟩의 로그‑로그 플롯에서 p가 p_c에 가까워질수록 초기 기울기(서브‑디퓨전 지수) α가 감소하고, p가 높을수록 α가 0.5에 근접한다. 특히 p=0.65 근처에서는 t≈10³~10⁴ 단계에서 α≈0.35에서 α≈0.48로 전이하는 ‘크로스오버’가 관찰되었다.

하지만 이 전이는 개별 구성에서만 뚜렷하게 나타났고, 모든 구성에 대해 평균을 취하면 ⟨r²(t)⟩는 단일 지수 곡선을 보이며 명확한 전이가 사라졌다. 저자들은 이를 ‘구성 의존적 전이’라 명명하고, 평균화 과정에서 짧은 거리의 구조적 불균일성(예: 작은 고립 클러스터, 구멍)이 소멸돼 전이 현상이 흐릿해진다고 해석한다. 이는 실험적 혹은 자연계에서 측정된 확산 데이터가 종종 평균화된 샘플에 기반함을 고려하면, 이론적 전이가 관측되지 않을 가능성을 시사한다.

이러한 발견은 퍼콜레이션 매체에서의 확산 모델링에 두 가지 중요한 교훈을 제공한다. 첫째, 미시적 구성을 고려한 시뮬레이션이 전이 현상을 포착하는 데 필수적이며, 둘째, 실험 설계 시 충분히 큰 시료 또는 개별 구성을 별도로 분석해야 이론적 스케일링을 검증할 수 있다. 또한 LGA가 복잡한 네트워크에서의 확산을 효율적으로 재현한다는 점은 향후 다공성 재료, 전해질, 생물학적 조직 등에서의 비정상 확산 연구에 유용한 도구가 될 것이다.

**