속도 정렬에 의한 집단 운동의 평균장 이론

초록

본 논문은 2차원에서 속도 정렬 상호작용을 갖는 자가구동 브라운 입자(Active Brownian Particles)를 모델링하고, 미시적 라그랑지안 방정식으로부터 비선형 포커-플랑크 방정식을 유도한 뒤, 확률분포의 모멘트 전개를 통해 평균장 이론을 구축한다. 두 종류의 추진(마찰) 함수—Rayleigh‑Helmholtz형과 선형 Schienbein‑Gruler형—에 대해 정렬 전이의 연속·불연속 특성을 분석하고, 평균 속도와 별도로 정의한 ‘유효 온도’(속도 변동성) 필드가 정렬 정도에 따라 어떻게 감소하는지를 보여준다.

상세 분석

이 연구는 개체 기반 모델과 연속체(평균장) 모델 사이의 연결 고리를 체계적으로 제시한다는 점에서 이론 생태학·물리학 분야에 큰 의미가 있다. 먼저, 입자 i의 위치 r_i와 속도 v_i에 대해 질량 m=1로 정규화한 Langevin 방정식(1)·(2)를 설정한다. 여기서 핵심은 속도‑의존성 마찰 γ(v)·v이며, 두 형태를 고려한다. (i) Rayleigh‑Helmholtz 마찰 −γ(v)v=(α−βv²)v (α,β≥0) 로, 저속에서는 음의 마찰(가속) → 자가구동, 고속에서는 양의 마찰(감속) → 속도 제한을 만든다. (ii) 선형 Schienbein‑Gruler 마찰은 γ(v)=γ₀−γ₁v 형태로, 비슷한 자가구동 특성을 갖지만 비선형 항이 없으므로 속도 분포가 Gaussian에 가깝다.

속도 정렬 상호작용은 ‘극성(polar)’ 형태로, 각 입자는 반경 ε 내 평균 속도 u_ε,i와의 차이에 비례하는 힘 µ(u_ε,i−v_i) 를 받는다. µ=1/τ_a 로 정의되며, τ_a는 정렬 시간 스케일이다. 이때 u_ε,i는 입자들의 속도 평균을 적분으로 표현한 비선형 연산자이며, 평균장 근사(N_ε≫1) 하에서 u_ε,i≈u(r,t) 로 치환한다.

포커-플랑크 방정식(8)은 확률밀도 P(r,v,t)의 시간 진화를 기술하고, 비선형 항 µ(u−v)·∇_v P 로 인해 비선형 Fokker‑Planck 형태가 된다. 여기서 모멘트 정의(4)와 연산을 통해 무한 계층의 모멘트 방정식(12)을 얻는다. n번째 모멘트는 (n+1)번째 모멘트에 의존하므로 폐쇄(close) 가 필요하다. 저자는 평균 속도 u와 속도 편차 δv=v−u 로 분해하고, x와 y 방향의 편차가 독립이며 교차 평균이 0임을 가정한다(13). 이 가정은 대칭성(무질서 상태와 정렬 상태)과 수치 검증에 의해 정당화된다.

그 결과, 공분산 행렬은 대각 형태가 되며, 대각 원소를 ‘온도’ T_x, T_y 로 정의한다. 이는 전통적인 열역학 온도와는 달리 비평형 정지 상태에서 서로 다를 수 있다. 모멘트와 u, T 사이의 관계식(15a‑d)은 가우시안 근사와 유사하지만, 4차 모멘트에 θ_k 라는 고차 변동성을 추가한다.

Rayleigh‑Helmholtz 마찰을 적용한 경우, 평균장 방정식은 다음과 같이 정리된다.

• 질량 보존: ∂_t ρ + ∇·(ρu)=0
• 운동량 방정식: ∂_t(ρu) + ∇·(ρu⊗u) = ρ