[论文解读] Dielectric response of confined water films: Insights from classical DFT

Daniel Borgis, Damien Laage+2 more|arXiv (Cornell University)|Feb 7, 2023

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 $ 的非单调行为，归因于单层或双层膜中极化响应与硬球堆积的相互作用。
• 该模型成功再现了实验中低厚度下介电响应的饱和现象，并为现象学解释中使用的三电容模型提供了物理解释基础。