두 종 비셰크 모델에서 인구·속도·소음 이질성이 군집 안정성에 미치는 영향

초록

본 연구는 상호 반대 정렬(inter‑species anti‑alignment) 상호작용을 갖는 두 종 비셰크 모델에 인구 비율, 자가 추진 속도, 공간적 활동 구배, 그리고 소음 강도 차이와 같은 네 가지 이질성을 도입하여, 평행(Parallel) 및 반평행(Anti‑parallel) 군집 상태의 존재와 전이를 정량적으로 조사한다. 인구 불균형이 커지면 소수 종은 밴드 형성을 잃고 다수 종의 단일 종 군집으로 전이되며, 속도 차이는 반평행 군집을 크게 촉진한다. 그러나 속도 차이가 공간적으로 구분된 경우, 빠른 영역이 넓을수록 평행 군집이 유지된다. 소음 이질성에서는 낮은 소음(‘차가운’) 종이 더 빠른 밴드 속도를 보이며 일시적으로 평행 군집을 파괴하지만, 장기적으로는 다시 회복된다. 이러한 결과는 내재적·외재적 이질성이 다중 종 활성 물질의 집단 운동을 근본적으로 재구성할 수 있음을 보여준다.

상세 분석

논문은 기존의 reciprocal two‑species Vicsek model(TSVM)을 기반으로 네 가지 이질성을 체계적으로 탐구한다. 첫 번째는 인구 이질성으로, 두 종의 입자 수 N_A와 N_B가 다를 때를 의미한다. 인구 비율 m_0=(N_A−N_B)/N이 증가하면 소수 종(B)의 밴드 수가 급격히 감소하고, 일정 임계값(m_0≈0.6) 이상에서는 B가 전혀 밴드를 형성하지 못하고 가스 상태에 머문다. 이때 전체 시스템은 다수 종(A)의 단일 종 군집(SSF)으로 전이되며, 평행(PF)과 반평행(APF) 상태가 동시에 존재하던 확률 분포 P(v_s,v_a)도 하나의 피크로 수렴한다.

두 번째는 속도 이질성이다. 동일한 인구(N_A=N_B)에서 v_A≠v_B일 경우, 빠른 종이 형성하는 밴드가 상대 종보다 더 강하게 전파되어 반평행(APF) 상태가 우세해진다. 시뮬레이션 결과는 v_A/v_B 비율이 1.5 이상이면 APF의 전역 정렬 파라미터 ⟨v_a⟩이 ⟨v_s⟩보다 크게 되며, PF는 거의 사라진다. 이는 서로 반대 방향으로 정렬하려는 상호작용이 속도 차이에 의해 비대칭적으로 강화되기 때문이다.

세 번째는 공간적 활동 구배(스페이셜 이질성)이다. 시뮬레이션 박스를 두 영역으로 나누어 한쪽에서는 A가 빠르고 B가 느리며, 반대쪽에서는 그 순서가 뒤바뀌는 설정을 적용한다. 이때 전체 평균 속도는 동일하지만, 빠른 영역이 전체 면적의 60 % 이상을 차지하면 PF 상태가 유지된다. 즉, 빠른 종이 차지하는 공간이 넓을수록 두 종이 같은 방향으로 이동하는 평행 밴드가 안정화된다. 이는 환경적 제약이 군집 형태를 결정짓는 중요한 매개변수임을 시사한다.

네 번째는 소음 이질성이다. η_A≠η_B로 두 종에 서로 다른 각도 잡음 강도를 부여한다. 낮은 소음(‘차가운’) 종은 더 큰 정렬 파라미터 ⟨v⟩와 높은 밴드 전파 속도를 보이며, 초기에는 PF 상태가 깨지고 APF‑like 전이 현상이 나타난다. 그러나 장시간 시뮬레이션( t≈10^6)에서는 두 종이 다시 같은 방향으로 정렬해 ⟨v_s⟩>⟨v_a⟩인 PF 상태로 복귀한다. 이는 소음 차이가 일시적인 동역학적 불안정을 야기하지만, 장기적으로는 상호 반대 정렬 상호작용이 시스템을 평행 정렬로 끌어당긴다는 점을 보여준다.

전반적으로 논문은 order parameter v_s=|∑_i σ_i|/N와 v_a=|∑_i s_i σ_i|/N을 이용해 PF와 APF를 구분하고, 각 이질성에 대한 파라미터 스캔을 통해 상전이 경계와 안정성 영역을 정량화한다. 특히 인구 불균형과 속도 차이는 시스템을 APF‑dominant 혹은 SSF‑dominant 영역으로 강하게 이동시키는 반면, 공간적·소음 이질성은 환경적 조건에 따라 PF를 유지하거나 일시적으로 파괴한다는 복합적인 메커니즘을 제시한다.