체인 빨아들이기: 흐름 구동 나노채널로의 폴리머 전이 동역학

초록

2차원 좁은 채널에 흐름에 의해 구동되는 폴리머 사슬이 채널 안으로 들어가는 과정을, 전단력 F가 채널 내 단위 모노머당 일정하다고 가정하고 Langevin 동역학 시뮬레이션과 이론적 스케일링을 통해 분석하였다. 채널 폭 R가 충분히 넓어 사슬이 채널 안에서 정상 상태 속도를 유지하면 전이 시간 τ는 τ ∝ N/F 로 스케일링한다. R이 작아 정상 상태가 형성되지 못하면 τ는 N·F에 대해 비보편적 의존성을 보인다. 또한 각 모노머가 채널 입구를 통과하는 대기시간 분포와 전이 좌표 s에 대한 최대 대기시간 s_max의 변화를 조사하였다.

상세 분석

본 연구는 전통적인 ‘단일 포어’ 전이 모델과 달리, 긴 나노채널 내부에서 흐름에 의해 발생하는 구동력이 전이된 모노머 수 s(t)에 비례한다는 새로운 물리적 상황을 설정하였다. 이를 위해 2차원에서 Lennard‑Jones 및 FENE 포텐셜로 연결된 bead‑spring 사슬을 사용하고, Langevin 방정식으로 각 bead에 보존력, 마찰력, 랜덤 힘을 포함시켜 시뮬레이션하였다. 구동력은 F · s(t) ẑ 형태로 적용되며, 이는 전기장 혹은 유체 흐름에 의해 채널 내부에서만 작용한다.

스케일링 분석에서는 먼저 사슬이 채널 폭 R보다 큰 블롭(blob) 구조를 형성한다는 가정을 도입하였다. 블롭당 포함되는 monomer 수 g ∼ (R/σ)^{1/ν_2D}이며, 전체 블롭 수 n_b = N/g ∼ N(σ/R)^{1/ν_2D}이다. 이때 자유에너지 비용은 F_conf ∼ N(σ/R)^{1/ν_2D} k_BT 로 추정된다. 블롭 picture를 이용해 사슬의 종축 길이 R_k ∼ n_b R ∼ N R^{-1/3} 로 얻어지며, 이는 약한 구동 하에서 전이 시간 τ_ss ∼ R_k/F ∼ N R^{-1/3}/F 로 전개된다. 구동력이 충분히 강하면 블롭이 파괴되고 사슬이 거의 완전히 늘어나므로 R_k ∼ N이 되고, 최종적으로 τ ∝ N/F 가 된다.

또한 전이 과정에서 발생하는 마찰력과 엔트로피적 저항을 정량화하였다. 전이쪽 마찰력 F_trans,f ≈ ξ s(t) dR_k/dt ≈ C R^{-1/3} ξ s(t) ds/dt 로 표현하고, 엔트로피 저항 F_trans,e ≈ 2.12 k_BT/R 로 근사하였다. 구동력과 마찰·엔트로피 저항의 균형식 C R^{-1/3} ξ s ds/dt + F_trans,e = s F 를 풀면 τ에 대한 식 (7)·(9)를 얻는다. 여기서 중요한 점은 s(t) 가 충분히 빨리 정상 상태 속도에 도달하면 F_trans,e 를 무시하고 τ ∝ N/F 로 단순화된다는 것이다.

시뮬레이션 결과는 이론적 예측을 뒷받침한다. 채널 폭 R = 4.5σ, 구동력 F = 0.1–0.5 ε/σ에 대해 전이 확률은 F에 거의 선형적으로 증가했으며, 긴 사슬(N ≥ 128)에서는 전이 시간이 N에 대해 선형 스케일링을 보였다. 반면 R이 작아 블롭 형성이 억제되면 전이 시간은 N·F에 대해 비선형적인 의존성을 나타냈다. 대기시간 분포를 조사한 결과, 전이 좌표 s에 대한 대기시간 w(s)는 초기 구간에서 급격히 증가했다가 s_max 에서 최대에 도달하고, 이후 감소한다. s_max 은 구동력 F와 채널 폭 R이 커질수록 좌측으로 이동하며, 이는 더 많은 모노머가 동시에 채널 안에 있을수록 전이 속도가 가속된다는 물리적 의미와 일치한다.

요약하면, 이 연구는 “채널 내부에서 구동력이 모노머 수에 비례하는” 새로운 전이 메커니즘을 제시하고, 블롭 이론과 마찰·엔트로피 균형을 통해 전이 시간의 스케일링을 정량적으로 설명하였다. 또한 대기시간 분석을 통해 전이 과정의 미시적 동역학을 상세히 파악함으로써, 나노플루이딕스 및 전기구동 나노채널 설계에 실질적인 가이드라인을 제공한다.