나노채널에서의 비정상적 모세관 충전과 습윤성 전환

초록

본 연구는 물과 모델 액체 금속을 탄소 나노튜브(CNT) 내부에 삽입해 분자동역학(MD) 시뮬레이션을 수행함으로써, 반나노미터 이하의 채널에서 모세관 충전 속도와 상승 높이가 튜브 반경에 따라 비단조적으로 변하고, 특정 반경에서는 습윤성이 친수성에서 소수성으로 전환되는 현상을 발견하였다. 이러한 현상은 전통적인 연속체 모델이 설명하지 못하는 ‘분리압(디스조인 프레셔)’에 기인한다는 점을 제시한다.

상세 분석

본 논문은 나노스케일에서의 모세관 현상이 연속체 유체역학과 표면장력 모델을 넘어서는 복잡한 현상임을 체계적으로 입증한다. 먼저, 전통적인 연속체 접근법을 정리하면서, 모세관 압력 차 Δp_men = 2γ cosθ / a와 입구 손실을 나타내는 Sampson식 Δp_ent = C η Q / a³을 결합해 속도 v_c = 4 Δγ / (π C η) 를 도출한다. 여기서 C는 원자 규모 입구 형상에 따라 반경 의존성을 갖는 계수이며, 이는 유한 요소 해석을 통해 사전에 계산된다. 이 식은 나노튜브 내부에서 플러그 흐름(plug‑flow) 가정하에 반경에 무관한 일정 속도를 예측한다는 점에서, 실제 MD 결과와의 차이를 검증하는 기준이 된다.

MD 시뮬레이션에서는 TIP4P/2005 물 모델과 AMBER96 탄소‑산소 상호작용 파라미터를 사용해 3.9 Å~24 Å 반경의 CNT를 10 nm 길이로 설정하고, 물이 한쪽 저장소에서 다른쪽 저장소로 흐르는 과정을 관찰한다. 동적 단계에서는 물 분자 수 N(t)의 선형 증가율을 통해 충전 속도 v_c를 측정하고, 정적 단계에서는 피스톤에 작용하는 힘을 이용해 실제 메니스쿠스 압력 차 Δp_men을 직접 계산한다.

결과적으로 a_c > 15 Å(≈1 nm)에서는 연속체 예측과 실험값이 일치하지만, a_c < 7.5 Å 구간에서는 속도가 비단조적 변동을 보이며, a_c ≈ 5.1 Å와 3.9 Å에서 국소 최대값을, a_c ≈ 4.3 Å와 4.7 Å에서는 음의 속도(역류)를 나타낸다. 정적 측정에서도 동일한 비정상적 압력 변동이 관찰되어, 원인이 입구 손실이 아니라 메니스쿠스 압력 자체의 변동임을 확인한다. 이는 물이 나노튜브 내부에서 층화·구조화되면서 발생하는 ‘분리압(디스조인 프레셔)’이 기존의 라플라스 압력에 추가되어, 효과적인 접촉각을 반전시키는 메커니즘으로 작용한다는 것을 의미한다.

이 현상의 일반성을 검증하기 위해, 방향성 결합이 없는 모델 액체 금속(임베디드 원자 모델 기반)으로 동일한 시뮬레이션을 수행하였다. 금속의 점도와 표면장력은 물과 크게 차이나지만, 반경이 1 nm 이하인 경우 동일한 비정상적 속도와 압력 변화를 보였으며, 특히 a_c ≈ 4.3 Å~4.7 Å 구간에서 역류가 발생한다. 따라서 물 특유의 수소결합이 아니라, 유체가 강하게 구조화되는 일반적인 현상이 원인임을 확인한다.

추가로, 슬릿(평판) 형태의 나노채널에서도 동일한 비단조적 압력 변화를 관찰했으며, 이를 기존의 ‘분리압’ 이론(디스조인 프레셔)과 정량적으로 일치시켰다. 즉, 유체가 제한된 공간에 들어갈 때 발생하는 격자 구조와 밀도 진동이 자유 에너지에 기여해, 효과적인 접촉각을 변동시키고, 때로는 친수성에서 소수성으로 전환시킨다.

이러한 발견은 나노다공성 매질, 멤브레인 선택성, 나노유체 에너지 저장 등 다양한 응용 분야에 중대한 영향을 미친다. 기존의 연속체 기반 설계법으로는 예측이 불가능한 ‘반전 습윤성’ 현상을 고려해야 하며, 설계 단계에서 분리압을 정량화하는 새로운 모델링 접근이 필요함을 시사한다.