[Paper Review] 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 InteractionsChemistry3 citations

TL;DR

This study uses classical density functional theory (DFT) with coupled number and polarization densities to analyze the dielectric response of water in nanoscale confinement. It shows that the apparent low dielectric constant (~2–3) in thin films (≤1 nm) arises not from unique water properties but from general confinement effects, with interfacial layers dominating the response in perpendicular fields, while transverse responses reach bulk-like behavior at much smaller thicknesses due to parallel field alignment.

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.

Motivation & Objective

• To re-express the dielectric response of confined water using a simplified two-variable classical DFT framework based on number and polarization densities.
• To clarify the physical meaning and validity of the local dielectric constant in confined, inhomogeneous fluids.
• To resolve the apparent contradiction between experimental observations of low dielectric constants in thin water films and the expectation of bulk water’s high polarity.
• To distinguish between intrinsic solvent behavior and geometric/confinement effects in dielectric response.

Proposed method

• A two-variable classical DFT formalism is employed, coupling number density and polarization density to describe inhomogeneous polar fluids under confinement.
• The theory derives a coarse-grained dielectric constant by smoothing over a solvent particle size, enabling a meaningful local dielectric response despite microscopic ill-definition near walls.
• The longitudinal (perpendicular) field case is modeled using a three-capacitor-in-series equivalent, with interfacial thickness $ h_i \approx 10\,\text{\AA} $ and effective dielectric constant $ \epsilon_i \approx 2-3 $, explaining slow convergence to bulk values.
• The transverse (parallel) field case is analyzed via a three-capacitor-in-parallel model, showing rapid convergence to bulk dielectric response at slab widths of ~10 nm.
• Analytical solutions and numerical solutions are derived without statistical noise, enabling systematic study across various slab thicknesses.
• The approach is validated by comparison with experimental data from Fumagalli et al. (2018) and consistent with prior MD simulations, while avoiding their convergence issues.

Experimental results

Research questions

• RQ1 Can a meaningful local dielectric constant be defined for water in confinement, despite its microscopic ill-definition near interfaces?
• RQ2 Why does the dielectric constant of water films drop to ~2–3 at thicknesses below 1 nm, and is this behavior specific to water or general to polar fluids?
• RQ3 How does the dielectric response differ between longitudinal (perpendicular) and transverse (parallel) field orientations in confined films?
• RQ4 What is the role of interfacial layers with low effective dielectric constants in determining the macroscopic dielectric response of thin films?
• RQ5 To what extent do geometric confinement and hard-sphere packing effects, rather than hydrogen bonding disruption, explain the observed dielectric anomalies?

Key findings

• The effective dielectric constant in the longitudinal direction reaches bulk values only for slab thicknesses much larger than twice the interfacial thickness $ h_i \approx 10\,\text{\AA} $, explaining the slow increase with thickness.
• The observed dielectric constant of ~2–3 for films ≤1 nm is not a unique property of water but a general consequence of confinement, arising from interfacial layers with $ \epsilon_i \approx 2-3 $, independent of bulk water’s high polarity.
• In the transverse configuration, the effective dielectric constant reaches bulk values at slab widths of approximately 10 nm, significantly shorter than in the longitudinal case.
• The interfacial region with low dielectric response dominates the total capacitance in the longitudinal case due to the $ 1/\epsilon $ dependence, making the system highly sensitive to interfacial properties.
• The non-monotonic behavior of $ \bar{\epsilon}_\perp(h) \approx 2 $ at sub-1 nm thicknesses is attributed to the interplay between polarization response and hard-sphere packing in single- or double-layer films.
• The model successfully reproduces the experimental saturation of dielectric response at low thicknesses and provides a physical basis for the three-capacitor models used in phenomenological interpretations.