반대·수직 구동 이진 입자계의 래닝·히스테리시스와 음의 드래그 현상

초록

입자 간 반발력이 있는 이진 시스템에서 한 종류는 +x, 다른 종류는 –x(반대 구동) 혹은 +y(수직 구동) 방향으로 외부 힘을 가한다. 시뮬레이션 결과, 반대 구동에서는 높은 구동력에서 상분리된 래닝 상태가 형성되고, 구동을 증가‑감소 사이클링할 때 속도‑구동 곡선과 6배 조정 결함 비율에 강한 히스테리시스가 나타난다. 수직 구동에서는 저구동에서 정체(jammed) 상태가 존재하고, 구동이 증가하면 무질서 흐름 또는 기울어진 래닝 상태로 전이되며, 이 전이 역시 히스테리시스를 보인다. 특히, 래닝이 기울어지면서 한 종류는 다른 종류와 반대 방향으로 움직이는 ‘음의 드래그’ 현상이 관찰된다. 입자 밀도와 구동 강도에 따라 전이 구간과 히스테리시스 폭이 달라지며, 고정된 x‑구동에 y‑구동을 점진적으로 가할 경우 락‑언락 전이가 여러 차례 발생한다.

상세 분석

본 연구는 2차원 주기적 경계 조건을 갖는 Coulomb 상호작용 입자 N개를 두 종류(A, B)로 나누어, A는 +x, B는 –x(반대 구동) 혹은 +y(수직 구동) 방향으로 외부 힘 F_D를 가하는 모델을 구축하였다. 과잉 감쇠(damped) 운동 방정식 η·dR_i/dt = –∑_{j≠i}∇V(R_ij)+F_D·σ_i·â+F_D·(1–σ_i)·b̂을 이용해 분자동역학 시뮬레이션을 수행했으며, 구동을 0→F_max→0 순으로 0.001 단위씩 증가·감소시키는 사이클을 적용하였다. 주요 관측 변수는 종(species)별 평균 속도 ⟨V_x⟩, ⟨V_y⟩와 Voronoi 분석을 통한 6배 조정 입자 비율 P_6이다.

반대 구동 경우, 저구동(F_D<0.365)에서는 입자들이 서로 교차하지 못하고 정체(jammed) 상태를 이루며 ⟨V⟩≈0이다. 구동이 임계값을 초과하면 입자들이 풀려 흐름(fluctuating) 상태로 전이하고, 이때 P_6는 0.5 수준으로 감소한다. 구동이 충분히 커지면 상분리된 래닝(laned) 구조가 형성되며, 같은 방향으로 이동하는 입자들이 좁은 띠를 이루어 충돌이 최소화된다. 이때 ⟨V⟩는 급격히 상승하고 P_6는 0.8에 근접한다. 중요한 점은 구동을 감소시킬 때 래닝 상태가 훨씬 낮은 구동력(F_D≈0.28)까지 유지되어 히스테리시스 루프가 형성된다는 것이다. 밀도 ρ가 증가할수록 래닝 전이 구간이 오른쪽으로 이동하고, 고밀도(ρ≈0.93)에서는 정체 구간이 넓어지며, 래닝이 사라지는 구동값도 크게 늘어난다.

수직 구동에서는 초기 저구동에서 입자들이 45° 방향으로 이동하려는 경향을 보이며 정체 상태가 유지된다. 구동이 증가하면 정체‑플로우(decoupling) 전이가 일어나 무질서 흐름이 나타나고, 이후에는 기울어진 래닝(lane tilted) 상태가 나타난다. 이때 래닝은 x‑축과 y‑축 사이의 각도가 변하면서, 한 종류는 다른 종류와 반대 방향으로 이동하는 음의 드래그 현상이 관찰된다. 특히, 고정된 x‑구동(F_A) 위에 y‑구동(F_B)을 서서히 증가시킬 때, 락‑언락 전이가 연속적으로 발생하며, 속도-구동 곡선에 급격한 점프와 하강이 교차한다. 이러한 전이들은 P_6의 급격한 변동과 일치하며, 전이 전후의 구조적 차이를 Voronoi 셀 분포로 확인하였다.

전반적으로, 이진 시스템에서 구동 방향을 바꾸는 것만으로도 다양한 비평형 위상(정체, 무질서 흐름, 래닝, 기울어진 래닝)과 그 사이의 1차 전이(히스테리시스가 뚜렷함)를 유도할 수 있음을 보여준다. 이는 스키르미온‑스키르미온늄 혼합물, 보행자 흐름, 전기‑중력 구동 콜로이드 등 실험적 시스템에 적용 가능하며, 특히 음의 드래그 현상은 마찰 제어 및 물질 운반에 새로운 설계 원리를 제공한다.