Continuum model of the simple dielectric fluid: consistency between density based and continuum mechanics methods

Force density functional theory in- and out-of-equilibrium

When a fluid is subject to an external field, as is the case near an interface or under spatial confinement, then the density becomes spatially inhomogeneous. Although the one-body density provides much useful information, a higher level of resolution is provided by the two-body correlations. These give a statistical description of the internal microstructure of the fluid and enable calculation of the average interparticle force, which plays an essential role in determining both the equilibrium and dynamic properties of interacting fluids. We present a theoretical framework for the description of inhomogeneous (classical) many-body systems, based explicitly on the two-body correlation functions. By consideration of local Noether-invariance against spatial distortion of the system we demonstrate the fundamental status of the Yvon-Born-Green (YBG) equation as a local force-balance within the fluid. Using the inhomogeneous Ornstein-Zernike equation we show that the two-body correlations are density functionals and, thus, that the average interparticle force entering the YBG equation is also a functional of the one-body density. The force-based theory we develop provides an alternative to standard density functional theory for the study of inhomogeneous systems both in- and out-of-equilibrium. We compare force-based density profiles to the results of the standard potential-based (dynamical) density functional theory. In-equilibrium, we confirm both analytically and numerically that the standard approach yields profiles that are consistent with the compressibility pressure, whereas the force-density functional gives profiles consistent with the virial pressure. For both approaches we explicitly prove the hard-wall contact theorem that connects the value of the density profile at the hard-wall with the bulk pressure. The structure of the theory offers deep insights into the nature of correlation in dense and inhomogeneous systems.

Chemomechanical equilibrium at the interface between a simple elastic solid and its liquid phase

Applying diffusion coupled deformation theory, we investigate how the elastic properties of a solid body are modified when forced to keep its chemical potential aligned with that of its melt. The theory is implemented at the classical level of continuum mechanics, treating materials as simple continua defined by uniform constitutive relations. A phase boundary is a sharp dividing surface separating two continua in mechanical and chemical equilibrium. We closely follow the continuum theory of the swelling of elastomers (gels) but now applied to a simple two phase one-component system. The liquid is modeled by a local free energy density defining a chemical potential and hydrostatic pressure as usual. The model is extended to a solid by adding a non-linear shear elastic energy term with an effective modulus depending on density. Imposing chemomechanical equilibrium with the liquid reservoir reduces the bulk modulus of the solid to zero. The shear modulus remains finite. The stability of the hyper-compressible solid is investigated in a thought experiment. A mechanical load is applied to a rectangular bar under the constraint of fixed lateral dimensions. The linear elastic modulus for axial loading is evaluated and found to be larger than zero, implying that the bar, despite the zero bulk modulus, can support a weight placed on its upper surface. The weight is stabilized by the induced shear stress. The density dependence of the shear modulus is found to be a second order effect reducing the density of the stressed solid (chemostriction).

Electric-field-based Poisson-Boltzmann theory: Treating mobile charge as polarization

Mobile charge in an electrolytic solution can in principle be represented as the divergence of ionic polarization. After adding explicit solvent polarization a finite volume of an electrolyte can then be treated as a composite nonuniform dielectric body. Writing the electrostatic interactions as an integral over electric-field energy density we show that the Poisson-Boltzmann functional in this formulation is convex and can be used to derive the equilibrium equations for electric potential and ion concentration by a variational procedure developed by Ericksen for dielectric continua [J. L. Ericksen, Arch. Rational Mech. Anal. 183, 299 (2007)AVRMAW0003-952710.1007/s00205-006-0042-4]. The Maxwell field equations are enforced by extending the set of variational parameters by a vector potential representing the dielectric displacement which is fully transverse in a dielectric system without embedded external charge. The electric-field energy density in this representation is a function of the vector potential and the sum of ionic and solvent polarization making the mutual screening explicit. Transverse polarization is accounted for by construction, lifting the restriction to longitudinal polarization inherent in the electrostatic potential based formulation of Poisson-Boltzmann mean field theory.