미생물 성장과 사멸이 만든 진동: 자원 적응과 시간 스케일의 교차점

초록

본 논문은 단일 미생물 집단이 지속적으로 공급되는 기질을 소비하면서 성장·사멸하는 최소 소비자‑자원 모델을 제시한다. 분석을 통해 미생물이 자원 변동에 적응하는 시간 스케일(τ₍food₎)과 사멸·회복 시간 스케일(τ₍death₎)의 비율이 β = τ₍food₎/τ₍death₎ ≥ 1일 때만 감쇠 진동이 발생함을 밝혀냈다. 영양 재활용, 다중 자원 이용 등 생물물리적 요소가 τ₍food₎를 단축하면 진동 영역이 축소되고 감쇠 속도가 빨라진다. 반대로 여러 필수 자원을 동시에 요구하면 진동 가능성이 확대된다.

상세 분석

이 연구는 미생물-기질 시스템을 두 변수 X(t) (미생물 밀도)와 S(t) (기질 농도)로 기술하고, 성장률 r·µ(S)와 사멸률 λ을 포함한 연속 방정식(1‑2)을 사용한다. 비차원화 과정에서 τ = λt, x = X·K_S/γ, s = S/K_S 로 변환하면 시스템은 두 비차원 파라미터 r/λ와 n/λ에만 의존하게 된다. 고정점 (x*, s*)는 µ⁻¹(λ/r) 형태로 정의되며, 여기서 µ(s)는 단조 증가 함수(예: Monod)이다. Jacobian 행렬을 전개하면 고유값 λ̃± = −β ± √(β²−4β)/2 로 나타나며, β = (r/λ)·µ′(s*)·x* 로 정의된다. β가 4 이하이면 복소 고유값을 갖고 감쇠 진동이 가능하고, β ≤ 2이면 진동 주파수가 감쇠율보다 커서 실제 눈에 띄는 진동을 만든다. β > 4이면 고유값이 모두 실수이며 과감쇠(over‑damped) 상태가 된다. 따라서 β는 “자원 적응 시간(τ₍food₎) 대비 사멸·복구 시간(τ₍death₎)”의 비율로 해석될 수 있다.

다양한 모델 확장을 통해 이 시간 스케일 비율이 어떻게 변하는지를 조사하였다. (i) necromass 재활용을 도입하면 추가적인 기질 공급 항이 생겨 µ′(s*)가 증가하고, 결과적으로 τ₍food₎가 감소한다. 이는 β를 감소시켜 진동 영역을 축소하고 감쇠 속도를 높인다. (ii) 보완적 다중 자원(서로 대체 가능한) 사용 역시 총 자원 공급 효율을 높여 τ₍food₎를 단축한다. 반면, 여러 자원이 동시에 필수적(필수적 동시 사용)인 경우, 각 자원의 제한이 누적되어 τ₍food₎가 늘어나 β가 커지고, 진동 가능 영역이 확대된다. (iii) 외부 흐름에 의한 기질 제거(chemostat‑like dilution)를 포함하면 자원 감소 속도가 빨라져 τ₍food₎가 다시 짧아지며, 진동이 억제된다.

또한, 결정론적 모델에서 β ≤ 4인 영역에서는 외부 잡음(예: 기질 유입 변동)으로 인해 지속적인 진동이 관찰될 수 있음을 시뮬레이션으로 확인했다. 이는 실제 실험실 연속배양(CSTR)이나 자연계 토양·수중 환경에서 흔히 보고되는 미세한 진동 현상과 일맥상통한다.

핵심적인 통찰은 복잡한 생태·생리학적 메커니즘을 단일 비차원 파라미터 β로 압축할 수 있다는 점이다. 이는 실험 설계 시 성장률 r, 사멸률 λ, 기질 공급 n을 조절함으로써 원하는 동적 특성(예: 빠른 안정화 vs. 진동 유지)을 예측하고 제어할 수 있음을 의미한다.