[논문 리뷰] Dielectric response of confined water films: Insights from classical DFT

Daniel Borgis, Damien Laage외 2명|arXiv (Cornell University)|2023. 02. 07.

Electrostatics and Colloid InteractionsChemistry인용 수 3

한 줄 요약

이 연구는 수소 결합의 붕괴가 아닌 일반적인 협착 효과에 기인한 것으로 나타나는, 나노스케일 협착 상태에서의 물의 유전 응답을 분석하기 위해 수소 결합 밀도와 극성 밀도를 연계한 고전적 밀도함수이론(DFT)을 사용한다. 얇은 필름(≤1 nm)에서 관찰되는 명백한 낮은 유전율(약 2–3)은 물의 고유한 성질 때문이 아니라 일반적인 협착 효과에 기인하며, 수직 방향 필드에서는 표면계층이 응답을 지배한다. 반면, 평행 필드 방향에서는 유전율이 훨씬 더 작은 두께에서 균일한 상태에 도달한다.

ABSTRACT

We re-examine the problem of the dielectric response of highly polar liquids such as water in confinement between two walls using a simple two-variable density functional theory involving number and polarisation densities. In the longitudinal polarisation case where a perturbing field is applied perpendicularly to the walls, we show that the notion of local dielectric constant, although ill-defined at a microscopic level, makes sense when a coarse-graining over the typical size of a particle is introduced. The approach makes it possible to study the effective dielectric response of thin liquid films of various thicknesses in connection to the recent experiments of [Fumagalli et al. , Science, 2018, 360, 1339-1342], and to discuss the notion interfacial dielectric constant. We argue that the observed properties as function of slab dimension, in particular the very low dielectric constants of the order of 2-3 measured for thin slabs of 1 nm thickness do not highlight any special property of water but can be recovered for a generic polar solvent having similar particle size and the same high dielectric constant. Regarding the transverse polarisation case where the perturbing field is parallel to the walls, the associated effective dielectric constant as a function of the slab dimension reaches bulk-like values at much shorter widths than in the longitudinal case.

연구 동기 및 목표

• 수소 결합 밀도와 극성 밀도를 기반으로 한 단순화된 이변수 고전적 DFT 프레임워크를 사용하여 협착된 물의 유전 응답을 재구성한다.
• 협착된 비균일한 유체에서 국소 유전율의 물리적 의미와 타당성을 명확히 한다.
• 박막 물 필름에서 관측된 낮은 유전율과 균일한 물의 높은 극성에 기대하는 기대치 사이의 명백한 모순을 해결한다.
• 유전 응답에서 고유한 용매 거동과 기하학적/협착 효과를 구분한다.

제안 방법

• 수소 결합 밀도와 극성 밀도를 연계하여 협착 조건에서 비균일한 극성 유체를 기술하는 이변수 고전적 DFT 형식을 적용한다.
• 용매 입자 크기의 평균화를 통해 군집화된 유전율을 유도함으로써, 벽 근처에서 미세 구조적으로 불명확한 국소 유전율 응답에도 의미 있는 정의를 가능하게 한다.
• 종방향(수직) 필드의 경우, 세 개의 커패시터를 직렬로 연결한 등가 모델을 사용하며, 표면계층 두께 $ h_i \approx 10\,\text{\AA} $ 와 효과적 유전율 $ \epsilon_i \approx 2-3 $ 를 고려하여 균일한 값으로 수렴하는 데 느린 경향을 설명한다.
• 횡방향(수평) 필드의 경우, 세 개의 커패시터를 병렬로 연결한 모델을 분석하여, 약 10 nm의 필름 두께에서 균일한 유전율 응답으로의 신속한 수렴을 보여준다.
• 통계적 노이즈 없이 해석적 해와 수치적 해를 유도함으로써 다양한 필름 두께에 대한 체계적 분석이 가능하다.
• 이전의 분자 동역학 시뮬레이션과 일치하며 수렴 문제를 피하는 바, Fumagalli 등(2018)의 실험 데이터와의 비교를 통해 모델의 타당성을 검증하였다.

실험 결과

연구 질문

• RQ1 벽 근처에서 미세 구조적으로 불명확한 국소 유전율이지만, 협착 상태의 물에 대해 의미 있는 국소 유전율을 정의할 수 있는가?
• RQ2 두께가 1 nm 이하일 때 물 필름의 유전율이 약 2–3으로 떨어지는 이유는 무엇이며, 이는 물에 특유한 현상인가, 아니면 극성 유체에 일반적인 현상인가?
• RQ3 협착 필름에서 종방향(수직)과 횡방향(수평) 필드 방향 간의 유전 응답은 어떻게 다를까?
• RQ4 낮은 효과적 유전율를 가지는 표면계층이 얇은 필름의 거시적 유전 응답을 결정하는 데 어떤 역할을 하는가?
• RQ5 수소 결합의 붕괴가 아니라 기하학적 협착과 하드 스피어 패킹 효과가 관측된 유전 이방성 현상을 설명하는 데 얼마나 기여하는가?

주요 결과

• 종방향에서의 효과적 유전율는 표면계층 두께 $ h_i \approx 10\,\text{\AA} $ 의 두 배보다 훨씬 큰 필름 두께에서야 균일한 값에 도달하며, 이는 두께 증가에 따른 느린 증가 경향을 설명한다.
• 두께가 1 nm 이하인 필름에서 관측된 약 2–3의 유전율는 물의 고유한 성질이 아니라 일반적인 협착 효과의 결과이며, $ \epsilon_i \approx 2-3 $ 인 표면계층에 기인하며, 균일한 물의 높은 극성과는 무관하다.
• 횡방향 구성에서는 약 10 nm의 필름 두께에서 효과적 유전율가 균일한 값에 도달하며, 종방향의 경우보다 훨씬 짧은 두께에서 이루어진다.
• 종방향의 경우, $ 1/\epsilon $ 의 의존성으로 인해 낮은 유전율 응답을 보이는 표면계층이 총 용량을 지배하므로, 시스템은 표면계층의 특성에 매우 민감하다.
• 1 nm 이하 두께에서의 $ \bar{\epsilon}_\perp(h) \approx 2 $ 인 비단조화적 행동은 단일 또는 이중층 필름에서 극성 응답과 하드 스피어 패킹의 상호작용에 기인한다.
• 모델은 낮은 두께에서의 유전 응답 포화를 성공적으로 재현하며, 현상학적 해석에서 사용된 세 커패시터 모델의 물리적 기반을 제공한다.