[論文レビュー] 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）の実験データとの比較により、本手法の妥当性が検証され、先行するMDシミュレーションとも整合するが、それらの収束問題を回避している。

実験結果

リサーチクエスチョン

• RQ1 界面付近での微視的定義の曖昧さがあるにもかかわらず、閉じ込められた水における意味のある局所誘電率を定義できるか？
• RQ2 膜厚が1 nm未満で誘電率が約2–3に低下するのはなぜか。この挙動は水に特有のものなのか、それとも極性流体に一般的な現象なのか？
• RQ3 閉じ込められた膜において、縦方向（垂直）と横方向（平行）の電場配置で誘電応答にどのような違いが生じるか？
• RQ4 低有効誘電率を示す界面層が、薄膜のマクロな誘電応答を決定づける役割を果たすのか？
• RQ5 水素結合の破壊よりも、幾何的閉じ込めやハードスフィアパッキング効果が、観察された誘電異常を説明する上でどの程度重要か？

主な発見

• 縦方向の有効誘電率は、界面層厚さ $ h_i \approx 10\,\text{\AA} $ の2倍以上も大きなスラブ厚さに達するまでバルク値に達しない。これは、膜厚が増加するに従ってゆっくりと増加する理由を説明する。
• 膜厚が1 nm未満のとき観察される約2–3の誘電率は、水固有の性質ではなく、界面層（$ \epsilon_i \approx 2-3 $）に起因する一般の閉じ込め効果である。これはバルク水の高い極性とは無関係である。
• 横方向配置では、有効誘電率がスラブ幅約10 nmでバルク値に達する。これは縦方向の場合に比べてはるかに短い厚さである。
• 縦方向の場合、誘電率 $ 1/\epsilon $ の依存性により、低誘電率応答を示す界面領域が全静電容量を支配し、システムが界面特性に極めて敏感であることが明らかになった。
• 1 nm未満の膜厚で $ \bar{\epsilon}_\perp(h) \approx 2 $ の非単調な挙動は、単層または二層膜における極性応答とハードスフィアパッキングの相互作用に起因する。
• 本モデルは、膜厚が小さい場合の誘電応答の飽和をうまく再現し、現象論的解釈で用いられる三コンデンサモデルの物理的根拠を提供している。