공정 순차 모드의 일반화 커뮤니케이션 P‑시스템과 인구 프로토콜의 통합

초록

**
이 논문은 일반화 커뮤니케이션 P‑시스템(GC‑P‑시스템)에 새로운 파생 모드인 공정 순차 모드(fs‑mode)를 도입하고, 이를 인구 프로토콜(PP)과 연결한다. 공정성을 만족하는 순차적 적용 규칙을 통해 PP를 GC‑P‑시스템의 한 변형으로 표현함을 보이며, 확률적 실행 전략이 Gillespie 알고리즘과 동등함을 증명한다. 무한 인구 한계에서 시스템의 동역학을 ODE로 기술하고, 역변환 조건과 Lotka‑Volterra·대수수 근사 등 사례를 제시한다.

**

상세 분석

**
논문은 먼저 기존의 GC‑P‑시스템 정의를 재정리한다. 각 셀은 라벨 i(0≤i≤n)와 다중집합 z_i 로 표현되고, 규칙은 (a,i)(b,j)→(a,k)(b,l) 형태의 동시 이동을 허용한다. 전통적인 파생 모드는 ‘최대 병렬(maximally parallel)’이었지만, 저자는 ‘공정 순차 모드(fs‑mode)’를 제안한다. 이 모드에서는 매 단계마다 하나의 규칙만 적용되며, 실행 전체가 공정성(fairness) 조건을 만족해야 한다. 즉, 어떤 구성 C가 무한히 반복되면 C에서 도달 가능한 모든 후속 구성도 무한히 나타나야 한다는 의미다.

이 공정 순차 모드는 인구 프로토콜(PP)의 파생 방식과 구조적으로 동일하다. PP는 유한 상태 집합 Q와 전이 관계 δ⊆Q⁴ 로 정의되며, 두 에이전트가 만나 (q₁,q₂)→(q₀₁,q₀₂) 로 상태를 바꾼다. 저자는 각 상태 q∈Q 를 라벨 q 를 가진 셀에 토큰 • 하나씩 배치함으로써 PP를 단일 알파벳 GC‑P‑시스템으로 인코딩한다. 규칙 집합 R은 δ의 모든 전이를 (•,q₁)(•,q₂)→(•,q₀₁)(•,q₀₂) 로 변환한다. 이렇게 하면 초기 다중집합 w_q = •^{|C₀|_q} 로 초기 구성을 설정하고, 출력 셀은 필요 없으므로 h=1 로 잡는다. 결과적으로 PP의 모든 공정 실행은 GC‑P‑시스템의 공정 순차 실행과 일대일 대응한다.

다음으로 저자는 확률적 파생 전략을 도입한다. 각 규칙에 반응 상수 c_μ 를 부여하고, 현재 구성 S 에서 가능한 적용 횟수 h_S^μ 를 계산하면, ‘propensity’ a_S^μ = c_μ·h_S^μ 가 정의된다. 이때 전체 propensities a_S⁰ = Σ_μ a_S^μ 를 이용해 다음 반응과 발생 시간을 샘플링하면, 이는 Gillespie의 SSA와 정확히 동일한 확률 과정을 만든다. 따라서 GC‑P‑시스템을 이 방식으로 실행하면, 연속‑시간 마코프 체인 모델을 구현하게 된다.

무한 인구 한계(N→∞)에서는 각 셀의 토큰 수를 농도 x_i = |z_i|/N 로 정규화하고, 질량 작용 법칙을 가정하면 a_S^μ/N → k_μ·Π_i x_i^{α_i} 형태가 된다. 여기서 α_i 는 규칙 μ 에서 소비되는 토큰 수이다. 이렇게 하면 확률적 미분 방정식이 결정론적 ODE 시스템으로 수렴한다는 표준 결과를 그대로 적용할 수 있다. 저자는 이 과정을 역으로 이용해, 주어진 ODE 시스템이 특정 형태의 GC‑P‑시스템(특히 농도‑의존적 확률 규칙)으로 구현 가능한 충분조건을 제시한다.

마지막으로 논문은 두 가지 응용 사례를 제시한다. 첫 번째는 Lotka‑Volterra 포식자‑피식자 모델을 GC‑P‑시스템 규칙으로 구현하고, 두 번째는 대수적 수(예: √2)의 근사값을 무한히 반복되는 공정 실행을 통해 계산하는 방법이다. 두 사례 모두 확률적 실행이 Gillespie SSA와 일치함을 실험적으로 확인하고, ODE 근사와의 차이를 정량화한다.

이러한 기여는 (1) PP와 GC‑P‑시스템 사이의 구조적·동역학적 동등성을 명확히 하고, (2) GC‑P‑시스템에 확률적·연속적 해석 도구를 도입함으로써 생물·화학 시스템 모델링에 새로운 접근법을 제공한다는 점에서 의의가 크다.

**