마이크로채널 내 반강성 고분자의 흐름 유도 이동과 회전 동역학

초록

본 연구는 원통형 마이크로채널에서 포아송 흐름에 노출된 반강성 고분자를 다중입자 충돌 동역학(MPC)와 분자동역학(MD) 결합 시뮬레이션으로 조사한다. 고분자 길이가 지속길이의 절반인 경우, 채널 중심부에서는 U자형 구조가 나타나고, 반경에 따라 정렬 및 구부러짐이 크게 변한다. 수소역학적 상호작용(HI)이 존재할 때 강한 횡방향(크로스‑스트림) 이동이 발생하며, 이는 고분자 방향성에 따른 확산 이방성에 기인한다. 또한, 고분자는 회전(tumbling) 운동을 보이며, 회전 시간은 Peclet 수에 대해 전단 흐름과 동일한 스케일링을 따른다.

상세 분석

이 논문은 반강성 고분자(길이 L_r ≈ 13a, 지속길이 L_p ≈ 2L_r)의 미세채널 내 동역학을 정량적으로 규명한다. 시뮬레이션은 MPC 기반 유체 모델에 고분자 구성을 MD로 결합한 하이브리드 방식을 사용했으며, HI를 켜고 끄는 두 경우를 비교함으로써 수소역학적 효과를 명확히 분리했다. 채널 반경 R = 8.5a, Peclet 수(Pe = γ̇τ_R) 를 0부터 360까지 변화시키며 고분자의 구조와 이동을 관찰하였다. 주요 결과는 다음과 같다.

1. 구조적 변이: 낮은 Pe에서는 고분자가 거의 등방적으로 배치되지만, Pe ≳ 50에서 중심부 근처에 U‑shaped(곡선) 구조가 빈번히 나타난다. 이는 굽힘 에너지 U_b가 중심부에서 상승함으로써 확인되며, 고분자 끝부분이 서로 가까워지는 ‘hairpin’ 형태가 전이적으로 관찰된다.

2. 방향성 및 정렬: 정렬 지표 S = (3cos²θ−1)/2 를 이용한 분석에서, 채널 벽에 가까울수록 S → 1(완전 정렬)이며, 중심부에서는 Pe가 증가함에 따라 S가 감소하고, Pe ≳ 100에서는 음의 값(≈−0.5)까지 내려가 U‑shaped 구조와 연관된다.

3. 횡방향 이동: 모노머 밀도 P_mo(r)와 중심질량 밀도 P_cm(r) 모두 Pe가 증가하면 중심부에서 감소하고, 일정 반경에서 피크를 형성한다. HI가 존재할 때는 피크가 더 뚜렷하고, 분포 폭 h_r²가 Pe≈50 이후 감소한다(내향 이동). HI를 제거하면 피크는 존재하지만 폭이 감소하지 않으며, 오히려 외향 이동이 강화된다. 이는 HI가 벽에 의한 ‘리프트’ 힘을 제공해 고분자를 중심부로 끌어당기는 역할을 함을 의미한다.

4. 회전(tumbling) 동역학: 고분자는 흐름에 의해 주기적인 회전 운동을 보이며, 회전 시간 τ_tumble ∝ Pe^−2/3 로 전단 흐름에서 관찰된 스케일링과 일치한다. 이는 고분자 길이가 지속길이와 비교적 짧아 전단에 의해 강하게 정렬·비틀림이 일어나기 때문이다.

5. 강성 의존성: 동일한 체인 길이에서 강성을 증가시키면 외향 이동이 더욱 두드러진다. 이는 강한 굽힘 저항이 고분자를 더 긴 ‘로드’ 형태로 유지시켜, 흐름에 의한 정렬이 강화되고, 결과적으로 횡방향 확산 이방성이 커지기 때문이다.

전반적으로, 반강성 고분자는 유체와의 수소역학적 상호작용에 의해 복합적인 구조·동역학 변화를 겪으며, Peclet 수와 채널 반경에 따라 중심부와 벽 근처에서 상반된 행동을 보인다. 이러한 현상은 마이크로플루이딕 디바이스에서 고분자(예: 액틴 필라멘트, 짧은 DNA) 운반 및 정렬을 설계하는 데 중요한 물리적 통찰을 제공한다.