활성·수동 입자 쌍의 결합 유도 운동 모델

초록

활성 입자와 수동 입자를 선형 스프링으로 연결하고, 활성 입자는 현재 속도 방향으로 자체 추진력을, 수동 입자는 활성 입자로부터 반발력을 받는 2차원 모델을 제안한다. 수치 시뮬레이션과 선형 안정성 분석을 통해 직선 이동, 원형 이동, 슬라럼(지그재그) 이동이라는 세 가지 대표적인 궤적을 확인하고, 자체 추진력 크기에 따라 직선‑원형 전이가 피치포크 분기 형태로 일어남을 보여준다.

상세 분석

본 논문은 두 입자(활성·수동) 사이의 상호작용을 최소화한 수학적 모델을 구축함으로써, 복잡한 화학·유체역학적 메커니즘을 추상화한다. 모델의 핵심 가정은 (1) 두 입자는 자연 길이 R과 강성 k를 갖는 선형 스프링으로 연결, (2) 활성 입자는 자신의 속도 벡터 방향으로 크기 f₂의 자기 추진력을 지속적으로 받으며, (3) 수동 입자는 활성 입자로부터 크기 f₁의 비상호적 반발력을 받는다는 점이다. 질량 m과 점성 저항 η는 양 입자에 동일하게 적용한다. 비차원화 과정을 거쳐 파라미터 η, f₁, f₂만을 남겨 두었으며, 이는 실험적 시스템(예: 캄퍼 디스크와 금속 와셔)에서 직접 조절 가능한 물리량에 대응한다.

운동 방정식은 중심질량(CM) 좌표 r, 상대 위치 ℓ, CM 속도 v, 상대 속도 w 로 재구성되어 6차원 동역학을 형성한다. 특히 ℓ=0(입자 겹침)과 v=w/2(활성 입자 정지)에서 특이점이 발생하지만, 실제 시뮬레이션에서는 이러한 상황을 회피한다. 수치 해석에서는 4차 Runge‑Kutta 방법을 사용해 Δt=10⁻⁴, t=10⁴까지 적분했으며, 초기 조건은 약간의 비대칭을 주어 다양한 궤적을 탐색했다.

시뮬레이션 결과는 f₂/f₁ 비율에 따라 네 가지 동역학적 패턴을 보여준다. (i) PPS(Passive‑Particle‑Preceding‑Straight)에서는 수동 입자가 앞서며 두 입자는 일직선으로 이동한다. (ii) PPC(Passive‑Particle‑Preceding‑Circular)와 (iii) APC(Active‑Particle‑Preceding‑Circular)에서는 각각 수동·활성 입자가 앞서며 원형 궤적을 그린다. (iv) SL(Slalom)에서는 활성 입자가 앞서며 파동형(지그재그) 궤적을 만든다. 각 패턴은 COM 속도와 ℓ과 v 사이 각도 φ, COM 속도 방향 ξ의 시간 변화를 통해 정량적으로 구분된다. φ가 0에 수렴하면 PPS, 0<φ<π/2이면 PPC, φ>π/2이면 APC, φ가 π를 중심으로 진동하면 SL으로 정의한다.

선형 안정성 분석을 통해 직선‑원형 전이가 피치포크 분기로 설명된다. f₂가 임계값을 초과하면 φ=0 해가 불안정해지고, 대칭이 깨진 두 개의 비대칭 해(φ≠0)가 나타난다. 이때 φ의 평균값이 π/2를 기준으로 앞서거나 뒤처지는 입자에 따라 PPC와 APC가 구분된다. 또한, SL과 APC 사이에 다중안정성(bistability) 영역이 존재함을 확인했으며, 경계 근처에서는 모호한(ambiguous) 동역학이 나타난다. 이러한 결과는 실험적 시스템에서 관찰된 전이 현상을 이론적으로 뒷받침한다.

논문은 모델이 매우 단순함에도 불구하고, 비대칭 추진력과 비상호적 반발력이라는 두 기본 물리량만으로 복잡한 집단 운동을 재현할 수 있음을 보여준다. 이는 화학적 구배, 표면 장력 구배, 혹은 유체 흐름에 의해 구동되는 다양한 활성 물질 시스템에 일반화될 수 있다. 향후 연구에서는 비선형 스프링, 시간 의존성 추진력, 다입자 상호작용 등을 포함해 더 풍부한 패턴(예: 클러스터링, 스웜 형성)과 연계할 수 있을 것이다.