Spontaneous charging and crystallization of water droplets in oil

Charge-Dipole Attraction as a Surface Interaction between Water Droplets Immersed in Organic Phases

The dynamic behavior of emulsion droplets during their interactions with one another or with solid surfaces plays a paramount role in their ultimate stability in various applications. While the interaction of oil droplets through a surrounding aqueous phase is well understood, recent studies on the interaction of water droplets through a surrounding pure organic phase showed the presence of an unexplained attraction between water droplets at relatively long ranges. In this research study, we propose fixed-surface-charge-bulk-dipole attraction as a new interaction force between water-in-oil droplets and then derive an equation for its disjoining pressure. The behavior of water droplets in the presence and absence of this charge-dipole interaction was numerically quantified using the Stokes-Reynolds-Young-Laplace model and compared to the experimental data. Numerically calculated net force curves are in excellent agreement with experimental data from the literature when charge-dipole attraction is included, while they deviate in its absence. In addition, the water droplet and thin oil film profiles in the presence and absence of charge-dipole attraction were calculated and compared. This research indicates that charge-dipole attraction can adequately explain the mysterious force observed in some studies, which demonstrates its unexplored potential to capture the physical properties and dynamic behavior of water droplets in organic phases with useful implications to unravel unidentified interactions between emulsion droplets in different industries.

Aqueous Nanoclusters Govern Ion Partitioning in Dense Polymer Membranes

The uptake and sorption of charged molecules by responsive polymer membranes and hydrogels in aqueous solutions is of key importance for the development of soft functional materials. Here, we investigate the partitioning of simple monatomic (Na⁺, K⁺, Cs⁺, Cl^-, I^-) and one molecular ion (4-nitrophenolate; NP^-) within a dense, electroneutral poly(N-isopropylacrylamide) membrane using explicit-water molecular dynamics simulations. Inside the predominantly hydrophobic environment, water distributes in a network of polydisperse water nanoclusters. The average cluster size determines the mean electrostatic self-energy of the simple ions, which preferably reside deeply inside them; we therefore find substantially larger partition ratios K ≃10^-1 than expected from a simple Born picture using a uniform dielectric constant. Despite their irregular shapes, we observe that the water clusters possess a universal negative electrostatic potential with respect to their surroundings, as is known for aqueous liquid-vapor interfaces. This potential, which we find concealed in cases of symmetric monatomic salts, can dramatically impact the transfer free energies of larger charged molecules because of their weak hydration and increased affinity to interfaces. Consequently, and in stark contrast to the simple ions, the molecular ion NP^- can have a partition ratio much larger than unity, K ≃10-30 (depending on the cation type) or even 10³ in excess of monovalent salt, which explains recent observations of enhanced reaction kinetics of NP^- reduction catalyzed within dense polymer networks. These results also suggest that ionizing a molecule can even enhance the partitioning in a collapsed, rather hydrophobic gel, which strongly challenges the traditional simplistic reasoning.

Simulation of charged systems in heterogeneous dielectric media via a true energy functional

For charged systems in heterogeneous dielectric media, a key obstacle for molecular dynamics (MD) simulations is the need to solve the Poisson equation in the media. This obstacle can be bypassed using MD methods that treat the local polarization charge density as a dynamic variable, but such approaches require access to a true free energy functional, one that evaluates to the equilibrium electrostatic energy at its minimum. In this Letter, we derive the needed functional. As an application, we develop a Car-Parrinello MD method for the simulation of free charges present near a spherical emulsion droplet separating two immiscible liquids with different dielectric constants. Our results show the presence of nonmonotonic ionic profiles in the dielectric with a lower dielectric constant.

Precipitation in aqueous mixtures with addition of a strongly hydrophilic or hydrophobic solute

We examine phase separation in aqueous mixtures due to preferential solvation with a low-density solute (hydrophilic ions or hydrophobic particles). For hydrophilic ions, preferential solvation can stabilize water domains enriched with ions. This precipitation occurs above a critical solute density n(p) in wide ranges of the temperature and the average composition, where the mixture solvent would be in a one-phase state without solute. The volume fraction of precipitated domains tends to zero as the average solute density n is decreased to np or as the interaction parameter χ is decreased to a critical value χ(p). If we start with one-phase states with n>n(p) or χ>χ(p), precipitation proceeds via homogeneous nucleation or via heterogeneous nucleation, for example, around suspended colloids. In the latter case, colloid particles are wrapped by thick wetting layers. We also predict a first-order prewetting transition for n or χ slightly below np or χ(p) for neutral colloids.

Curvature dependence of the electrolytic liquid-liquid interfacial tension

The interfacial tension of a liquid droplet surrounded by another liquid in the presence of microscopic ions is studied as a function of the droplet radius. An analytical expression for the interfacial tension is obtained within a linear Poisson-Boltzmann theory and compared with numerical results from nonlinear Poisson-Boltzmann theory. The excess liquid-liquid interfacial tension with respect to the pure salt-free liquid-liquid interfacial tension is found to decompose into a curvature-independent part due to short-ranged interfacial effects and a curvature-dependent electrostatic contribution. Several curvature-dependent regimes of different scalings of the electrostatic excess interfacial tension are identified. Symmetry relations of the interfacial tension upon swapping droplet and bulk liquid are found to hold in the low-curvature limit, which, e.g., lead to a sign change of the excess Tolman length. For some systems a low-curvature expansion up to the second order turns out to be applicable if and only if the droplet size exceeds the Debye screening length in the droplet, independent of the Debye length in the bulk.