패턴화된 확장형 나노포어를 통한 고분자 전이동 시간 스케일링 연구

초록

2차원 분자동역학 시뮬레이션으로 원통형·원뿔형 패턴 포어를 조사하였다. 고분자 사슬 길이 N과 포어 길이 Lₚ에 대해 평균 전이동 시간 ⟨τ⟩이 ⟨τ⟩ ∼ N³ F(Lₚ N¹·⁵) 로 스케일링됨을 확인했으며, γ = 3.00 ± 0.05, ϕ = 1.50 ± 0.05의 보편적 지수를 얻었다. 포어 패턴(전부 흡착, 중간 반발, 출구 반발)과 원뿔 각도 α에 관계없이 동일한 스케일링이 유지된다.

상세 분석

본 논문은 2차원 Langevin 동역학을 이용해 유연한 고분자 사슬이 확장된 패턴 포어를 통과하는 전이동 현상을 정량적으로 규명한다. 포어 형상은 원통형(α = 0)과 원뿔형(α > 0) 두 가지를 고려했으며, 포어 내부를 흡착·반발 구간으로 나누어 세 가지 패턴(A: 전구역 흡착, B: 양쪽 흡착·중간 반발, C: 입구 흡착·출구 반발)을 구현하였다. 외부 전기장에 대응하는 구동력 f_ext(x) = f₀ d/(d + 2x tanα) 를 적용해 원뿔형 포어에서는 구동력이 점진적으로 감소하도록 설계하였다.

시뮬레이션 파라미터는 σ = 1, ε = 1, K = 500, η = 1, k_BT = 1 로 설정했으며, 각 데이터 포인트는 1500–2000회의 성공적인 전이동 사건을 평균하였다. 먼저 포어 길이 Lₚ가 5σ~10σ 범위에서 α를 변화시켰을 때, 원통형 포어에서는 Lₚ가 길어질수록 ⟨τ⟩가 크게 감소한다는 점을 확인했다. 이는 포어 내부에서 구동력이 일정하게 작용해 전이동이 전반적으로 가속화되기 때문이다. 반면 원뿔형 포어에서는 짧은 Lₚ에서는 α가 증가할수록 출구가 넓어져 엔트로피적 완화가 지배, ⟨τ⟩가 감소한다. 그러나 Lₚ가 충분히 길어지면 구동력 감소가 지배하여 ⟨τ⟩가 다시 증가한다는 비단조적 거동을 보였다.

다음으로 사슬 길이 N에 대한 스케일링을 조사하였다. 동일한 Lₚ와 α = 0 조건에서 f₀ = 0.2, 0.6, 1.0 세 가지 구동력을 적용했을 때, ⟨τ⟩ ∝ N^β 형태의 지수 β가 구동력과 포어 길이에 따라 변한다는 것을 발견했다. 약한 구동력(f₀ = 0.2)에서는 β가 1.31.6 사이로 ν(=0.75)보다 작아 비평형 텐션 전파가 제한적임을 시사한다. 구동력이 강해질수록 β는 1.81.95까지 상승해, 기존 이론에서 제시된 1 + ν(≈1.75)보다 크게 되며, 이는 포어 내부에서의 강제적인 스트레칭이 사슬 전체에 전파되는 효과를 반영한다.

핵심적인 스케일링 가설 ⟨τ⟩ ∼ N^γ F(Lₚ N^ϕ) 를 검증하기 위해 ⟨τ⟩/N^γ를 Lₚ N^ϕ에 대해 플롯하였다. γ = 3.00 ± 0.05, ϕ = 1.50 ± 0.05 로 최적화했을 때, 모든 데이터(다양한 Lₚ, N, α, 포어 패턴, 구동력)가 하나의 마스터 커브에 성공적으로 수축되었다. 이는 포어 길이와 사슬 길이가 결합된 복합 스케일 변수 x = Lₚ N^ϕ 가 전이동 시간의 주요 제어인자를 제공한다는 강력한 증거이다. 또한, β와 γ, ϕ 사이의 관계 β ≈ γ − (ϕ p)/2 (p는 스케일 함수 F의 지수) 가 경험적으로 만족됨을 보고, 포어 패턴이 스케일링 지수에 미치는 영향이 미미함을 확인했다.

결론적으로, 포어의 기하학적 확장(Lₚ)과 고분자 사슬 길이(N)의 결합 효과가 전이동 시간에 지배적이며, 포어 내부의 흡착·반발 패턴이나 원뿔 각도 α에 관계없이 동일한 스케일링 법칙이 적용된다. 이는 실험적 나노포어 설계 시, 포어 길이와 사슬 길이의 비례 관계만을 조절하면 전이동 속도를 예측·제어할 수 있음을 의미한다.