좌절된 정렬이 만든 정지 육각 격자

초록

본 논문은 Vicsek‑Kuramoto 모델에 위상 지연 파라미터 α를 도입하여, α가 π/2 이하일 때는 전통적인 전역 동기화와 스워밍을 보이다가, α가 π/2를 초과하면 Hopf‑Turing 분기점을 통해 정지 육각 격자와 소용돌이·이중‑클러스터 격자 등 복합적인 시공간 패턴이 나타난다. 격자 형성은 결합 강도 K와 상호작용 반경 d₀에 의해 지배되며, 순수한 방향성 상호작용만으로도 대칭 격자를 만들 수 있음을 보여준다.

상세 분석

본 연구는 자기구동 입자들의 방향을 Kuramoto 위상으로 동일시한 Vicsek‑Kuramoto 프레임워크에 ‘좌절(frustration)’ 파라미터 α를 추가함으로써, 위상 지연이 집단 행동에 미치는 영향을 체계적으로 탐구한다. 모델식(1)은 입자 위치 rᵢ와 위상 θᵢ의 동역학을 각각 자기구동 속도 v와 회전 확산 Dᵣ, 그리고 이웃 입자와의 정렬 상호작용 K·F(θⱼ−θᵢ) 로 기술한다. 여기서 F(θ)=sin(θ+α)−sinα는 α=0이면 순수 정렬, α=π이면 완전 반정렬을 재현하도록 설계되었으며, −sinα 항은 θⱼ=θᵢ가 평형이 되도록 보정한다.

시뮬레이션 결과는 α에 따라 네 가지 주요 영역을 만든다. (1) α<π/2에서는 전역 위상 동기화(p_std≈1)와 공간적 균일성을 유지하는 전통적인 Vicsek 상태가 관찰된다. (2) α가 π/2에 가까워지면 스워밍 클러스터가 형성돼 위상 동기화는 유지되지만 입자 밀도가 국소적으로 증가한다. (3) α가 π/2를 초과하면 시스템은 Hopf‑Turing 분기점에 도달해 복소 고유값이 양의 실수부를 갖게 되고, 이는 주기적인 시공간 패턴, 즉 정지 육각 격자를 유도한다. 격자 내부에서는 두 종류의 호흡‑유사 운동이 나타난다. 첫 번째는 ‘소용돌이 격자’에서 입자들이 평균 반경 v/|⟨θ̇⟩|≈v/K 주변을 주기적으로 팽창·수축하고, 두 번째는 ‘이중‑클러스터 격자’에서 위상이 π 차이인 두 집단이 반대 방향으로 회전하면서 거의 고정된 원형 궤도를 그린다. 이러한 호흡 운동은 식(4)에서 유도된 평균 회전 속도 ⟨θ̇⟩=−K sinα에 의해 정량화된다.

선형 안정성 분석에서는 연속 분포 ρ(r,θ,t)의 연속 방정식(5)와 비선형 연산자 T를 도입해, 균일 해 ρ₀에 대한 작은 섭동 ε·e^{λt+i k·r}Φ(θ) 를 전개한다. 고유값 λ(k)의 실수부가 양이 되는 조건은 J₁(k d₀)·cosα<0 (식 14)이며, 이는 α=π/2가 임계값임을 직접적으로 보여준다. α<π/2에서는 k=0 모드가 가장 불안정해 전역 동기화만 발생하고, α>π/2에서는 특정 비영( k* )에서 복소 고유값이 발생해 Hopf‑Turing 패턴이 나타난다. 흥미롭게도 k*는 α에 크게 의존하지 않으며, 이는 격자 간격이 α>π/2 구간 전체에서 거의 일정함을 설명한다(그림 4b).

격자 구조의 ‘우세’와 ‘열위’는 K와 d₀에 의해 구분된다. (K, d₀) 파라미터 평면에서 격자 우세 영역은 높은 K와 중간~큰 d₀에서 나타나며, 여기서는 격자 상수 a가 d₀에 비례하고, 단위 셀 부피는 K가 커질수록 감소한다(식 15). 이는 격자 상수와 입자 이동 반경 v/K 사이의 기하학적 맞춤 조건 2π/k* = 2vK 로 정량화된다. 따라서 격자 형성은 패턴의 공간 주기와 입자 동역학 사이의 균형에 의해 결정된다.

결론적으로, 순수한 방향성 상호작용만으로도 복잡한 대칭 격자를 생성할 수 있음을 입증했으며, 좌절 파라미터 α가 새로운 설계 변수로 작용해 패턴 토폴로지를 제어한다는 점에서 활성 물질 설계에 중요한 통찰을 제공한다. 다만, 격자 선택 메커니즘, 호흡 주기의 정확한 예측, 그리고 잡음·주파수 분산이 포함된 실제 시스템에서의 안정성 등에 대한 추가 연구가 필요하다.