희소 상호작용이 만든 비가우시안 풍경: 일반화 로트카‑볼테라 모델의 위상과 다중 평형

초록

본 논문은 무작위·희소·대칭인 상호작용 행렬을 갖는 일반화 로트카‑볼테라(gLV) 모델의 평형 상태를 분석한다. 베일리프 전파와 캐비티 방법을 이용해 베타-정규화된 희소 네트워크(무작위 정규 그래프)에서 정확한 종 풍부도 분포를 구하고, 완전 연결 모델과는 달리 강한 비가우시안(감마형) 분포와 새로운 ‘위상적 유리상태(topological glass)’를 발견한다. 이 위상은 다중 평형이 존재하지만, 그 원천이 네트워크 토폴로지와 멸종 현상에 크게 의존한다는 점에서 기존 연구와 차별화된다.

상세 분석

논문은 먼저 전통적인 gLV 모델이 완전 연결(fully‑connected) 가정하에 해석될 때는 평균 상호작용 μ, 변동성 σ, 온도 T 세 파라미터에 의해 다중 평형(multiple‑equilibria) 단계가 나타난다는 점을 상기한다. 그러나 실제 생태계는 상호작용이 희소하고 대칭이 아닌 경우가 많아 이러한 가정이 현실을 반영하지 못한다는 문제 제기를 하고, 대칭이면서도 희소한 상호작용을 갖는 무작위 정규 그래프(k‑regular graph, 여기서는 k=3)를 모델링 대상으로 선택한다.

핵심 이론적 도구는 캐비티 방법과 베일리프 전파(BP) 방정식이다. 희소 그래프는 루프 길이가 log N에 비례해 무한대 N에서 사실상 트리 구조를 이루므로, 각 변수를 제외한 ‘캐비티’ 그래프에서의 마진을 재귀적으로 계산할 수 있다. 식(7)‑(9)는 이 마진 ηᵢ→ⱼ(nᵢ)와 최종 마진 ηᵢ(nᵢ)를 정의하고, 이를 통해 전체 Gibbs‑Boltzmann 분포 P(n)∝e^{−H_eff/T}의 정확한 1‑점 분포를 얻는다. 여기서 H_eff는 식(3)‑(6)에서 유도된 효과 해밀토니안이며, 상호작용 계수 α_{ij}는 평균 μ/k와 분산 σ²/k를 갖는 정규분포에서 추출된다.

수치적으로는 정수형 종 개체수 n_i∈ℕ을 사용해 계산 효율성을 높였으며, ε=10^{-4}와 같은 작은 정규화 파라미터를 도입해 n_i=0에서의 로그 항을 회피한다. 이산화가 결과에 미치는 영향은 부록 B에서 검증하였다.

주요 결과는 두 가지로 요약된다. 첫째, σ를 증가시킬수록 단일 평형 단계에서도 종 풍부도 분포 η(n)이 가우시안에서 감마형(또는 로그‑정규형)으로 급격히 변한다. 이는 희소 토폴로지가 평균‑제곱 변동성 외에 추가적인 ‘구조적 잡음’을 제공해 비가우시안 꼬리를 생성한다는 점을 시사한다. 둘째, σ=0(즉, 상호작용이 모두 동일)인 경우에도 μ와 T의 조합에 따라 ‘위상적 다중‑어트랙터’ 단계가 나타난다. 이 단계는 기존 완전 연결 모델에서 관찰되는 다중 평형과는 달리, 서로 다른 살아남은 종들의 서브그래프(연결성) 구조가 서로 다른 고정점으로 이어지는 ‘위상적 유리상태’를 의미한다. 특히 멸종(λ≈0) 조건이 존재해야 이 현상이 뚜렷해지며, 온도 T가 유한할 때도 지속된다.

이러한 발견은 (i) 희소성 자체가 비가우시안 풍부도와 새로운 위상 전이를 유도한다는 점, (ii) 기존의 ‘σ‑driven’ 다중 평형 메커니즘이 반드시 필요하지 않으며, (iii) 실제 생태계 데이터에서 관찰되는 로그‑정규·감마형 풍부도 분포를 설명하는 데 희소 네트워크가 핵심 역할을 할 수 있음을 뒷받침한다. 또한, 베일리프 전파와 캐비티 방법이 희소·대칭 gLV 모델의 정확한 정적 해를 제공한다는 점에서 이론적·계산적 방법론의 확장성을 보여준다.