거시진화 셀룰러 오토마타 모델

초록

본 논문은 단일 조상으로부터 돌연변이와 지역 상호작용을 통해 종이 생성·소멸하는 다종 생태계를 셀룰러 오토마타(CA)로 구현한다. Monte‑Carlo 시뮬레이션 결과, 종 다양성의 급격한 폭발, 주기적 대멸종, 군집 구조 형성 등 실제 진화 현상을 재현한다.

상세 분석

이 연구는 전통적인 생태·진화 모델이 고정된 종 수와 연속적인 미분 방정식에 의존하는 한계를 극복하고자, 이산적인 공간-시간 격자 위에 셀룰러 오토마타 규칙을 부여한다. 격자의 각 셀은 ‘생물체’ 혹은 ‘빈 공간’ 상태를 가질 수 있으며, 생물체는 종을 나타내는 정수 코드와 에너지(또는 적합도) 값을 보유한다. 초기 조건은 하나의 조상 종만 존재하도록 설정하고, 매 시간 단계마다 무작위로 선택된 셀에 대해 다음과 같은 로컬 규칙이 적용된다. 첫째, 인접 8셀(모르렐리 이웃) 중 살아있는 셀이 일정 수 이상이면 번식이 가능하고, 번식 시 자식 셀은 부모와 동일한 종 코드를 물려받는다. 둘째, 번식 과정에서 사전 정의된 돌연변이 확률 μ가 적용되어, 일정 확률로 종 코드가 변형(새로운 정수 코드)된다. 이는 새로운 종의 탄생을 의미한다. 셋째, 자원의 제한을 모사하기 위해 각 셀은 에너지 양을 소모하며, 인접 셀과의 경쟁에서 승리한 개체만 에너지를 유지하고, 패배자는 사라진다. 네번째, 환경 변동을 구현하기 위해 전역 파라미터(예: 자원 풍부도, 경쟁 강도)를 시간에 따라 조정한다.

Monte‑Carlo 시뮬레이션은 위 규칙을 수천에서 수만 세대에 걸쳐 반복함으로써, 종 다양성(Species Richness), 평균 수명, 군집 크기 분포 등의 통계량을 추출한다. 결과는 다음과 같은 특징적인 거시진화 현상을 보여준다. (1) 초기 단계에서 급격한 종 폭발이 일어나며, 이는 ‘적응 방사’와 유사한 패턴이다. (2) 일정 주기 혹은 외부 충격(환경 파라미터 급변) 후 대규모 멸종이 발생하고, 이후 새로운 종군이 재편성된다. (3) 종 간 경쟁과 공생 관계가 공간적으로 클러스터를 형성해, 군집 내 다양성은 높지만 군집 간 차별화가 뚜렷하다. (4) 돌연변이 확률 μ와 경쟁 강도 β의 조합에 따라 시스템은 ‘활성’(다양성 유지)와 ‘정체’(종 수 고정) 두 상전이 상태 사이를 전이한다.

이 모델의 강점은 (i) 종 수가 동적으로 변한다는 점, (ii) 로컬 상호작용만으로도 복잡한 거시적 패턴이 자발적으로 나타난다는 점, (iii) 파라미터 조정을 통해 다양한 생태·진화 시나리오를 실험할 수 있다는 점이다. 반면 한계점으로는 (a) 생물학적 현실성을 높이기 위한 유전적 구조(게놈)와 표현형 매핑이 단순화되어 있다, (b) 공간 차원(2D 격자)과 경계 조건이 실제 지구 규모의 복잡성을 충분히 반영하지 못한다, (c) 시간 스케일이 세대 단위로 고정돼 있어 장기적인 지질학적 변화를 정밀히 모델링하기 어렵다. 향후 연구에서는 다층 격자, 가변 환경 흐름, 그리고 유전 알고리즘 기반의 돌연변이 모델을 도입해 모델의 생물학적 타당성을 강화할 수 있다.