국소적으로 초균일한 상태에서의 거대한 밀도 변동

초록

활성 네마틱 입자들이 서로 접촉했을 때만 움직이는 규칙을 도입하면, 큰 스케일에서는 전형적인 거대한 밀도 변동(GNF)이 나타나면서도 중간 스케일에서는 밀도 변동이 억제되는 초균일성(hyperuniformity)이 동시에 발생한다. 저자들은 이 현상을 수치 시뮬레이션과 연속체 이론으로 분석하고, 임계점 근처에서 두 현상이 교차하는 길이 ξ가 발산한다는 점을 밝혀냈다.

상세 분석

본 연구는 두 가지 상반된 비평형 현상, 즉 방향성 질서가 있는 활성 물질에서 흔히 관찰되는 거대한 밀도 변동(GNF)과 흡수 전이(absorbing phase transition) 근처에서 나타나는 초균일성(hyperuniformity)을 하나의 모델에 결합한다. 이를 위해 저자들은 ‘네마틱 랜덤 오가니제이션 모델(NROM)’을 제안했는데, 핵심은 입자들이 서로 겹쳐 있을 때만 일정한 스텝 δ₀만큼 자신의 축 방향으로 이동하도록 하는 것이다. 입자들의 방향은 Vicsek‑type 규칙에 따라 잡음 σ를 포함해 업데이트되며, σ=0.1, δ₀=0.3, 상호작용 반경 Rₐₗ=1, 입자 직경 D=0.46 등 파라미터가 선택된다.

시뮬레이션 결과는 전통적인 ROM과 유사하게 입자 밀도 φ가 임계값 φ_c≈0.30984를 초과하면 활동 상태가 지속되고, φ<φ_c에서는 시스템이 흡수 상태에 고정된다. 활동 정도 A는 φ−φ_c에 대해 A∝(φ−φ_c)^β 형태로 증가하며, β≈0.63으로 2차원 보존된 유향 퍼콜레이션(CDP) 클래스와 일치한다. 흥미롭게도 φ>φ_c 구간에서도 네마틱 정렬이 강하게 나타나 스칼라 순서 파라미터 S>0을 유지한다.

밀도 변동을 박스 크기 ℓ에 대한 분산 ⟨Δn²⟩와 평균 입자 수 ⟨n⟩의 관계로 분석하면, 중간 스케일(ℓ≈ξ)에서는 ⟨Δn²⟩∝⟨n⟩^{α_hu} (α_hu≈0.77) 로 초균일성을 보이고, 큰 스케일에서는 ⟨Δn²⟩∝⟨n⟩^{α_gnf} (α_gnf≈1.65) 로 전형적인 GNF가 나타난다. 교차 길이 ξ는 φ→φ_c⁺에서 ξ∝(φ−φ_c)^{-μ} (μ≈0.625) 로 발산한다. 구조인자 S(q) 역시 q→0에서 λ_gnf<0 (GNF)와 중간 q에서 λ_hu>0 (초균일성) 두 개의 지수적 스케일을 보이며, 동일한 스케일링 관계 S(q)∝q^{λ}와 ξ의 발산을 확인한다.

연속체 이론적 접근에서는 밀도 ρ, 네마틱 텐서 Q, 활동 밀도 A 세 개의 필드를 도입하고, 보존된 밀도 흐름 J = D_ρ∇A + χ₁ Q·∇A + χ₂ A∇·Q 로 전개한다. ρ와 A는 각각 보존된 흐름과 비보존된 반응항(κρ−a)A−λA²을 갖는 CDP 형태의 동역학을 따른다. Q는 활성 입자의 확산에 의해 운반되며, χ₃∇∇A 항이 추가된다. 선형화 후 푸리에 변환을 수행하면 구조인자 S(q)의 스케일링식 (7)이 도출되며, 여기서 ξ∝Δφ^{-3/4} 로 예측된다. 이론적 지수 μ≈0.75와 ζ’≈0.5는 수치 결과와 비교적 근접해, 선형 이론이 핵심 메커니즘을 포착함을 시사한다.

핵심 메커니즘은 다음과 같다. 활동 입자만이 밀도 흐름을 생성하므로, 활동이 억제되는 중간 스케일에서는 흐름이 감소해 밀도 변동이 억제되고 초균일성이 나타난다. 반면, 네마틱 질서가 장거리 상관을 제공해 흐름이 비정상적으로 증폭되면 큰 스케일에서 GNF가 재현된다. 따라서 두 현상은 활동-밀도 결합과 네마틱 텐서의 곡률‑구동 전류가 경쟁하면서 발생한다.

이 연구는 (i) 환경‑의존적 활동 규칙이 흡수 전이와 초균일성을 유도하고, (ii) 네마틱 질서가 존재할 때 GNF와 초균일성이 공존할 수 있음을 실증적으로 보여준다. 이는 마이크로 로봇 군집, 사이클링 전단 입자, 퀸크 롤러 등 실험적 시스템에 적용 가능하며, 향후 흡수 전이와 방향성 질서가 동시에 존재하는 새로운 보편성 클래스를 탐구하는 길을 열어준다.