공간 비대칭을 고려한 포식자‑피식자 확률 셀룰러 오토마톤

초록

우리는 포식자‑피식자 시스템의 확률적 동역학을 분석하기 위해 확률 셀룰러 오토마톤을 도입한다. 각 종의 개체는 격자상의 사이트에 존재하며, 로트카‑볼테라 모델에 기반한 마코프식 로컬 규칙에 따라 상호작용한다. 이때 이웃 구조는 비대칭적(비등방성)이며, 즉 한 방향으로만 영향을 미치는 특성을 가진다. 우리는 이러한 공간적 비대칭성이 종들의 개체밀도 진동 특성에 미치는 영향을 조사한다. 연구는 단순 평균장(mean‑field)과 쌍 평균장(pair mean‑field) 근사를 이용하여 수행되었으며, 분석 과정에서 비대칭성을 명시적으로 고려하였다.

상세 요약

본 논문은 전통적인 포식자‑피식자 모델을 격자 기반의 확률 셀룰러 오토마톤(Probabilistic Cellular Automaton, PCA) 형태로 구현함으로써, 미시적 상호작용이 거시적 집단 동역학에 미치는 영향을 정량적으로 탐구한다. 가장 큰 특징은 이웃 정의에 비대칭성을 도입한 점이다. 일반적인 동등방성(isotropic) 모델에서는 각 사이트가 네 방향(또는 여덟 방향)으로 동일하게 영향을 주고받지만, 여기서는 예를 들어 오른쪽과 위쪽 이웃만을 고려하는 ‘편향된’ 이웃을 설정한다. 이러한 설계는 실제 생태계에서 서식지의 흐름, 풍향, 물류 경로 등 방향성 요인이 존재할 때 발생할 수 있는 비대칭적 상호작용을 모사한다는 점에서 의미가 크다.

모델의 동역학은 로트카‑볼테라 방정식의 확률적 버전을 격자에 적용한 것으로, 세 가지 기본 전이 규칙(피식자 번식, 포식자 사냥·번식, 포식자 사망)이 각각 일정 확률로 수행된다. 전이 확률은 해당 사이트와 선택된 이웃 사이트의 상태에만 의존하므로 마코프성(Markovian) 특성을 유지한다.

분석 방법으로는 먼저 ‘단순 평균장(Mean‑Field)’ 근사를 적용해 전체 시스템의 평균 밀도 방정식을 도출한다. 이 단계에서는 공간 상관을 완전히 무시하고, 각 사이트가 전체 평균 밀도에만 의존한다는 가정을 둔다. 결과적으로 비대칭성이 평균장 수준에서는 사라지는 경향이 나타나며, 전통적인 로트카‑볼테라의 주기적 진동과 유사한 해를 얻는다. 그러나 실제 시뮬레이션에서는 이러한 평균장 해가 과도하게 진동 주기를 과소평가하거나, 진동 폭을 과대평가하는 등 차이를 보인다.

이를 보완하기 위해 ‘쌍 평균장(Pair Mean‑Field)’ 근사를 도입한다. 여기서는 인접한 두 사이트 사이의 상관 함수 (P_{AB}) (A와 B는 각각 피식자, 포식자, 빈칸 등)의 동역학을 명시적으로 기술한다. 비대칭 이웃 구조가 쌍 평균장 방정식에 직접적으로 나타나며, 예를 들어 오른쪽 이웃과 위쪽 이웃에 대한 전이 확률이 서로 다르게 적용된다. 이때 얻어지는 고유값 분석 결과는 비대칭성에 의해 진동 주기가 방향에 따라 달라지는 ‘전파형 진동(Traveling Wave)’ 현상을 예측한다. 즉, 포식자와 피식자 밀도의 파동이 격자 상에서 한 방향으로 이동하면서 지속적인 주기를 유지한다는 것이다.

시뮬레이션 결과는 쌍 평균장 근사의 예측과 정량적으로 일치한다. 비대칭성이 강할수록 파동 속도가 증가하고, 특정 파라미터 구간에서는 파동이 소멸하고 고정점(steady state)으로 수렴한다. 이는 비대칭성이 시스템의 안정성 경계(stability boundary)를 이동시켜, 기존 동등방성 모델에서는 관찰되지 않던 새로운 동역학적 위상을 만든다.

이러한 발견은 이론적 생태학뿐 아니라, 전염병 확산, 화학 반응 확산, 사회적 의견 전파 등 방향성 상호작용이 중요한 복합 시스템에 대한 모델링에 시사점을 제공한다. 특히, 공간 비대칭성을 명시적으로 포함한 평균장 접근법은 계산 비용을 크게 늘리지 않으면서도, 실제 현상에 가까운 예측을 가능하게 한다는 점에서 실용적이다.