표면 특성이 나노채널 내 오일 수송을 지배: 분자동역학 시뮬레이션 통찰

초록

본 연구는 분자동역학 시뮬레이션을 활용하여 표면 상호작용이 나노채널 내 오일 수송에 미치는 영향을 규명했습니다. 채널 폭, 오일 극성, 재료 조성, 표면 거칠기가 모두 오일 분자의 질량중심(COM) 변위와 확산에 현저한 영향을 미침을 발견했습니다. 6 nm 채널에서 COM 변위는 2 nm 채널보다 30배 이상 컸으며, 채널 중심의 확산 계수는 표면 근처보다 거의 2배 높았습니다. 극성 오일 성분은 실리카 기판과 더 강한 정전기적 상호작용을 겪어 비극성 옥타데케인에 비해 이동도가 크게 감소했습니다. 또한 금 기판은 훨씬 높은 상호작용 에너지를 부과하여 오일을 거의 고정화시킨 반면, 방해석은 중간 정도의 변위를 허용했습니다. 마지막으로 나노스케일 공동은 국소 수송을 수정합니다: 얕은 공동은 채널을 효과적으로 넓혀 변위를 약간 증가시키는 반면, 깊은 공동은 흡착 증가로 흐름을 방해합니다. 이러한 발견은 제한된 기하학에서 오일 이동을 지배하는 메커니즘을 명확히 하고 효율적인 셰일 오일 추출 전략을 안내합니다.

배경

전 세계 에너지 수요 증가와 재래식 자원 고갈로 비재래식 셰일 오일이 방대한 매장량으로 인해 집중적인 연구 관심을 끌고 있습니다[1,2]. 셰일 오일 자원이 재래식 오일의 3배 이상임에도[2], 2-100 nm 범위의 나노스코픽 기공 내 갇혀 있어 개발 가능한 비율은 제한적입니다[4,5]. 이러한 구속은 밀도, 젖음성, 확산 같은 유체 특성을 변화시키는 현저한 표면 효과를 유도하여, 거시적 시스템과 현저히 다른 수송 거동을 초래합니다[10-12]. 예를 들어, 물은 동등한 거시적 튜브보다 탄소 나노튜브에서 더 빠르게 흐릅니다[6,7], 탄소 나노튜브의 친수/소수 특성은 직경에 따라 변합니다[9]. 선행 분자동역학 연구는 전단 응력과 기공 크기가 탄소 및 실리카 나노슬릿에서 물과 오일 역학에 영향을 미침을 보였습니다[13-15]. 그러나 표면 화학과 지형이 셰일 오일 수송을 어떻게 지배하는지에 대한 체계적 연구는 드뭅니다. 이러한 효과 이해는 추출 기술 최적화에 필수적이며 막 분리, 생체분자 수송, 에너지 저장 장치에도 광범위한 함의를 가집니다[16-20].

방법

모든 시뮬레이션은 Materials Studio(Accelrys Inc.)의 Discover 모듈과 COMPASS 전원자 힘장으로 수행되었으며, 이는 응축상에서 결합, 교차항, 비결합 상호작용을 정확히 포착합니다:

$$ {E}_{\mathrm{total}}={E}_{\mathrm{valence}}+{E}_{\mathrm{cross}\text{-}term}+{E}_{\mathrm{nonbond}} \tag{1} $$

각 항의 전체 표현은 원본 논문에 제시되어 있습니다. 반데르발스 힘은 르나르-존스 퍼텐셜(절단 15.5 Å)로 모델링했고, 정전기학은 에발트 방법으로 처리했습니다. 시뮬레이션은 앤더슨 서모스탯으로 298 K에서 NVT 앙상블로 실행했습니다. 모든 방향에 주기적 경계 조건을 적용했습니다. 데이터는 5 ps마다 샘플링했습니다.

셰일의 주요 광물인 실리카 기판은 (001) 표면에서 구성하고 지질 조건을 모방하기 위해 수산기화했습니다[26]. 두 개의 대향 표면(각 1.5 × 7 × 0.85 nm³)이 2, 4, 6 nm 폭의 슬릿 채널을 형성했습니다. 옥타데케인, 피리딘, 페놀—일반적인 셰일 오일 성분—을 도입하여 0.8 g cm⁻³ 벌크 밀도를 달성했습니다(2 nm 채널 40분자; 4 nm 및 6 nm 채널에서 피리딘 407, 페놀 344, 옥타데케인 80 및 120분자). 에너지 최소화와 500 ps 평형화 후, 채널 축을 따라 원자당 3.1×10⁻¹⁴ N의 일정한 체적 힘을 가해 분자를 구동했습니다. 이는 MD 수송 연구의 표준 접근법입니다[27-29].

결과 및 고찰

채널 폭의 효과

가해진 힘 하에서 채널을 통과하는 분자 수가 꾸준히 증가했습니다(보충 정보, Fig. S1). 그림 2는 더 넓은 채널이 2 ns 후 더 큰 COM 변위를 산출함을 보여줍니다. COM 위치는 다음과 같이 정의됩니다:

$$ {z}_{\mathrm{COM}}={\displaystyle \sum_i\frac{m_i}{M}{r}_i} \tag{5} $$

2 nm 채널에서 COM은 0.85 nm만 이동한 반면, 6 nm 채널에서는 변위가 30 nm를 초과했습니다. 이 현저한 차이는 더 좁은 기공에서 더 강한 흡착에서 기인하며, 평균 상호작용 에너지로 정량화됩니다:

$$ {E}_{\mathrm{average}\ \mathrm{interaction}}=\frac{E_{\mathrm{total}}-(E_{\mathrm{oil}}+E_{\mathrm{substrate}})}{N} \tag{6} $$

그림 3은 명확한 역상관을 보여줍니다: 상호작용 에너지가 높을수록 변위가 감소하며, 에너지 3배 증가가 COM 이동을 30배 이상 감소시킵니다. 상호작용 효과는 채널이 넓어짐에 따라 감소하여 나노스케일 수송의 크기 의존적 특성을 강조합니다.

층상 밀도 프로파일(그림 2)은 벽에 인접한 5 Å 두께의 계면층을 드러냅니다. 이 층 내 옥타데케인 배향은 주로 표면과 평행(θ ≈ 80°-90°)인 반면, 채널 중심의 분자는 무작위 배향을 보입니다. 평균제곱변위(식 7)에서 유도된 확산 계수는 방사 위치에 따라 변합니다: 중심에서 가장 높고 벽 쪽으로 감소합니다(그림 5a). 6 nm 채널은 가파른 포물선 프로파일을 보이는 반면, 2 nm 채널은 거의 균일하고 낮은 확산을 보여 피스톤 같은 전면을 연상시킵니다. 이러한 경향은 벽 흡착(표면 근처 "그림자" 영역에서 운동 억제)과 분자간 운동량 전달 간의 경쟁에 기인합니다.

극성의 효과

비슷한 분자 수의 페놀, 피리딘, 옥타데케인을 4 nm 실리카 채널에 넣은 시뮬레이션은 분자 극성의 역할을 보여줍니다. 동일한 힘에도 불구하고 옥타데케인의 COM 변위는 페놀이나 피리딘의 약 16배였습니다(그림 8). 총 상호작용 에너지(반데르발스 + 정전기)는 극성 종에서 현저히 높았으며, 정전기학이 총합의 약 절반을 기여했습니다(그림 8). 제1원리 계산된 쌍극자 모멘트는 옥타데케인(0.032 Debye)이 사실상 비극성인 반면, 페놀(1.306 Debye)과 피리딘(2.245 Debye)은 유의미한 쌍극자를 가짐을 확인합니다. 결과적으로 극성 분자는 실리카 표면에 더 강하게 결합하여 이동성을 감소시킵니다.

재료 조성의 효과

실리카, 금, 방해석 채널 비교(그림 9)는 금이 옥타데케인을 강하게 고정시키는 반면, 방해석은 중간 변위를 허용함을 보여줍니다. 그림 10은 COM 변위가 평균 상호작용 에너지와 상관관계가 있음을 나타냅니다: 실리카가 가장 낮은 에너지와 가장 높은 변위, 금이 반대, 방해석이 중간입니다. 실리카와 방해석의 에너지가 유사함에도 변위가 다른 것은 표면 원자 유형과 텍스처도 수송에 영향을 미침을 시사합니다.

표면 거칠기의 효과

실리카 표면에 나노스케일 공동(3 Å 또는 6 Å 깊이)을 도입하면 국소 흐름이 변경됩니다. 공동은 갇힌 오일 분자를 수용하여 인접 유체 층의 속도를 감소시킵니다(그림 11). 3 Å와 6 Å 공동의 COM 변위는 각각 3.95 Å와 3.07 Å로, 평평한 표면 값(3.17 Å)보다 얕은 공동에서 약간 높지만 깊은 공동에서 더 낮습니다. 얕은 공동에서의 향상은 채널 효과적 확장에 기인할 가능성이 높으며, 깊은 공동은 더 강한 흡착을 부과하여 흐름을 방해합니다.

결론

이 분자동역학 연구는 나노채널 폭, 오일 극성, 기판 재료, 표면 거칠기가 오일 수송을 집단적으로 결정함을 입증했습니다. 6 nm 채널은 2 nm 채널보다 30배 이상 큰 COM 변위를 제공하며, 채널 중심 확산 계수는 벽 근처보다 거의 2배 높습니다. 극성 분자와 고상호작용 기판은 이동성을 현저히 억제하는 반면, 얕은 표면 공동은 수송을 약간 향상시킬 수 있습니다. 이러한 통찰은 나노구조 추출 시스템 설계와 제한된 기하학에서 유체 수송이 중요한 광범위한 응용에 필수적입니다.

약어

COM:

질량중심

COMPASS:

응축상 최적화 분자 퍼텐셜(원자론적 시뮬레이션 연구용)

MD:

분자동역학

vdW:

반데르발스