Density functional theory for charged fluids

GCMe: Efficient Implementation of the Gaussian Core Model with Smeared Electrostatic Interactions for Molecular Dynamics Simulations of Soft Matter Systems

In recent years, molecular dynamics (MD) simulations have emerged as an essential tool for understanding the structure, dynamics, and phase behavior of charged soft matter systems. To explore phenomena across greater length and time scales in MD simulations, molecules are often coarse-grained for better computational performance. However, commonly used force fields represent particles as hard-core interaction centers with point charges, which often overemphasizes the packing effect and short-range electrostatics, especially in systems with bulky deformable organic molecules and systems with strong coarse-graining. This underscores the need for an efficient soft-core model to physically capture the effective interactions between coarse-grained particles. To this end, we implement a soft-core model uniting the Gaussian core model with smeared electrostatic interactions that is phenomenologically equivalent to recent theoretical models. We first parametrize it generically using water as the model solvent. Then, we benchmark its performance in the OpenMM toolkit for different boundary conditions to highlight a computational speedup of up to 34 × compared to commonly used force fields and existing implementations. Finally, we demonstrate its utility by investigating how boundary polarizability affects the adsorption behavior of a polyelectrolyte solution on perfectly conducting and nonmetal boundaries.

Boosting Salinity Energy Extraction Efficiency in Capacitive Mixing by Polyelectrolyte Surface Coating

The extraction of salinity gradient energy in the capacitive mixing (CapMix) technique can be enhanced by using polyelectrolyte-coated electrodes. The micromechanism of polyelectrolyte (PE) coating enhancing the salinity energy extraction is studied by using a statistical thermodynamic theory. When PE takes same charge sign as the coated electrodes, the extraction efficiency can be boosted owing to the enhanced response of electrical double layer (EDL) to external cell voltage (V0). For the optimal case studied, the extraction efficiency was boosted from 0.25 to 1.25% by PE coatings. Owing to counterion adsorption and the enhanced response of EDL, the extraction energy density presented a local maximum at V0 = 0, which is higher than another local maximum value when V0 ≠ 0. This provides important guidance on the two approaches of CapMix in terms of capacitive Donnan potential (CDP, V0 = 0) and capacitive double-layer expansion (CDLE, V0 ≠ 0). Under the effects of PE coating, the extraction efficiency by CDLE can be improved to about 11% by CDP for the optimal studied case. The synergistic effect of grafting conditions can significantly elevate the energy density and extraction efficiency of the CDP process.

Statistical field theory of polarizable polymer chains with nonlocal dipolar interactions

The electromechanical response of polymeric soft matter to applied electric fields is of fundamental scientific interest as well as relevant to technologies for sensing and actuation. Several existing theoretical and numerical approaches for polarizable polymers subject to a combined applied electric field and stretch are based on discrete monomer models. In these models, accounting for the interactions between the induced dipoles on monomers is challenging due to the nonlocality of these interactions. On the other hand, the framework of statistical field theory provides a continuous description of polymer chains that potentially enables a tractable way to account for these interactions. However, prior formulations using this framework have been restricted to the case of weak anisotropy of the monomer polarizability. This paper formulates a general approach based in the framework of statistical field theory to account for the nonlocal nature of the dipolar interactions without any restrictions on the anisotropy or nonlinearity of the polarizability of the monomer. The approach is based on three key elements: (1) the statistical field theory framework, in which the discrete monomers are regularized to a continuous dipole distribution, (2) a replacement of the nonlocal dipole-dipole interactions by the local electrostatics partial differential equation with the continuous dipole distribution as the forcing, and (3) the use of a completely general relation between the polarization and the local electric field. Rather than treat the dipole-dipole interactions directly, the continuous description in the field theory enables the computationally tractable nonlocal-to-local transformation. Further, it enables the use of a realistic statistical-mechanical ensemble wherein the average far-field applied electric field is prescribed, rather than prescribing the applied field at every point in the polymer domain. The model is applied, using the finite element method, to study the electromechanical response of a polymer chain in the ensemble with fixed far-field applied electric field and fixed chain stretch. The nonlocal dipolar interactions are found to increase, over the case where dipole-dipole interactions are neglected, the magnitudes of the polarization and electric field by orders of magnitude as well as significantly change their spatial distributions. Next, the effect of the relative orientation between the applied field and the chain on the local electric field and polarization is studied. The model predicts that the elastic response of the polymer chain is linear, consistent with the Gaussian approximation, and largely unchanged by the orientation of the applied electric field, though the polarization and local electric field distributions are significantly impacted.

Unraveling Flow Effect on Capacitive Energy Extraction from Salinity Gradients

The harvesting of salinity gradient energy through a capacitive double-layer expansion (CDLE) technique is directly associated with ion adsorption and desorption in electrodes. Herein, we show that energy extraction can be modulated by regulating ion adsorption/desorption through water flow. The flow effects on the output energy, capacitance, and energy density under practical conditions are systematically investigated from a theoretical perspective, upon which the optimal operating condition is identified for energy extraction. We demonstrate that the net charge accumulation displays a negative correlation with the water flow velocity and so does the surface charge density, and this causes a nontrivial variation in the magnitude of output energy when water flows are introduced. When high water flows are introduced in both the charging and discharging processes, the energy extraction can be significantly reduced by 47.69-49.32%. However, when a high flow is solely exerted in the discharging process, the energy extraction can be enhanced by 12.94-14.49% even at low operation voltages. This study not only offers a comprehensive understanding of the microscopic mechanisms of surface-engineered energy extraction with water flows but also provides a novel direction for energy extraction enhancement.

Spiers Memorial Lecture: Towards understanding of iontronic systems: electroosmotic flow of monovalent and divalent electrolyte through charged cylindrical nanopores

In many practical applications, ions are the primary charge carrier and must move through either semipermeable membranes or through pores, which mimic ion channels in biological systems. In analogy to electronic devices, the "iontronic" ones use electric fields to induce the charge motion. However, unlike the electrons that move through a conductor, motion of ions is usually associated with simultaneous solvent flow. A study of electroosmotic flow through narrow pores is an outstanding challenge that lies at the interface of non-equilibrium statistical mechanics and fluid dynamics. In this paper, we will review recent works that use dissipative particle dynamics simulations to tackle this difficult problem. We will also present a classical density functional theory (DFT) based on the hypernetted-chain approximation (HNC), which allows us to calculate the velocity of electroosmotic flows inside nanopores containing 1 : 1 or 2 : 1 electrolyte solution. The theoretical results will be compared with simulations. In simulations, the electrostatic interactions are treated using the recently introduced pseudo-1D Ewald summation method. The zeta potentials calculated from the location of the shear plane of a pure solvent are found to agree reasonably well with the Smoluchowski equation. However, the quantitative structure of the fluid velocity profiles deviates significantly from the predictions of the Smoluchowski equation in the case of charged pores with 2 : 1 electrolyte. For low to moderate surface charge densities, the DFT allows us to accurately calculate the electrostatic potential profiles and the zeta potentials inside the nanopores. For pores with 1 : 1 electrolyte, the agreement between theory and simulation is particularly good for large ions, for which steric effects dominate over the ionic electrostatic correlations. The electroosmotic flow is found to depend very strongly on the ionic radii. In the case of pores containing 2 : 1 electrolyte, we observe a reentrant transition in which the electroosmotic flow first reverses and then returns to normal as the surface change density of the pore is increased.

Enabling High-Capacitance Supercapacitors by Polyelectrolyte Brushes

Polyelectrolyte brushes (PEBs) hold excellent potential for designing high-capacitance electrical double-layer capacitors (EDLCs), a crucial component of supercapacitors. Both experiments and computational simulations have shown their energy-storage advantage. However, the effect of PEBs on the energy storage of EDLCs is not yet fully understood. Herein, we systematically study the energy-storage effects of polyanionic (PA) and polycationic (PC) brushes using polymer density functional theory (DFT). First, the application of polymer DFT in polyelectrolyte-grafted EDLCs is successfully validated using molecular dynamics simulations. With the help of polymer DFT, an interfacial adhesion microstructure of the PA/PC brushes is observed. Most importantly, the results show that polyelectrolyte-grafted EDLCs achieve a significant increase in capacitance at low salt concentrations and surface voltages, offering an excellent energy-storage advantage over traditional EDLCs. However, this advantage is considerably diminished at high salt concentrations or surface voltages, showing unusual salt- and voltage-dependent behaviors of energy-storage capacity. Nonetheless, the PC-grafted EDLCs maintain their outstanding energy-storage performance, even at relatively high salt concentrations and surface voltages. These findings deepen our comprehension of PEBs at the molecular level and provide insights for the molecular design of high-capacitance supercapacitors.

Phase Behavior of Ion-Containing Polymers in Polar Solvents: Predictions from a Liquid-State Theory with Local Short-Range Interactions

The thermodynamic phase behavior of charged polymers is a crucial property underlying their role in biology and various industrial applications. A complete understanding of the phase behaviors of such polymer solutions remains challenging due to the multi-component nature of the system and the delicate interplay among various factors, including the translational entropy of each component, excluded volume interactions, chain connectivity, electrostatic interactions, and other specific interactions. In this work, the phase behavior of partially charged ion-containing polymers in polar solvents is studied by further developing a liquid-state (LS) theory with local shortrange interactions. This work is based on the LS theory developed for fully-charged polyelectrolyte solutions. Specific interactions between charged groups of the polymer and counterions, between neutral segments of the polymer, and between charged segments of the polymer are incorporated into the LS theory by an extra Helmholtz free energy from the perturbed-chain statistical associating fluid theory (PC-SAFT). The influence of the sequence structure of the partially charged polymer is modeled by the number of connections between bonded segments. The effects of chain length, charge fraction, counterion valency, and specific short-range interactions are explored. A computational App for salt-free polymer solutions is developed and presented, which allows easy computation of the binodal curve and critical point by specifying values for the relevant model parameters.

Structural features and electrostatic energy storage of electric double layers in confined polyelectrolyte solutions under low-salt conditions

Self-consistent field theory is used to systematically study the characteristics and electrostatic energy storage of electric double layers in confined polyelectrolyte solutions for salt-free and low salt concentration systems.

Ising density functional theory for weak polyelectrolytes with strong coupling of ionization and intrachain correlations

We report a theoretical framework for weak polyelectrolytes by combining the polymer density functional theory with the Ising model for charge regulation. The so-called Ising density functional theory provides an accurate description of the effects of polymer conformation on the ionization of individual segments and is able to account for both the intra- and interchain correlations due to the excluded-volume effects, chain connectivity, and electrostatic interactions. Theoretical predictions of the titration behavior and microscopic structure of ionizable polymers are found to be in excellent agreement with the experiment.

Nonelectrostatic Adsorption of Polyelectrolytes and Mediated Interactions between Solid Surfaces

Polymer-mediated interaction between two solid surfaces is directly connected to the properties of the adsorbed polymer layers. Nonelectrostatic interactions with a surface can significantly impact the adsorption of polyelectrolytes to charged surfaces. We use a classical density functional theory to study the effect of various polyelectrolyte solution properties on the adsorption and interaction between two like-charged surfaces. Our results show that nonelectrostatic interactions not only enhance polyelectrolyte adsorption but can also result in qualitatively different salt effects with respect to the adsorbed amount. In particular, we observe decreasing, increasing, and a previously unreported nonmonotonic behavior in the adsorbed amount of polymer with added salt under the conditions studied, although the nonmonotonic regime only occurs for a narrow range in the parameter space. With sufficient nonelectrostatic adsorption, the adsorbed polymer layers produce a long-range repulsive barrier that is strong enough to overcome dispersive interactions that cause surfaces to attract. Concurrently, a short-range bridging attraction is observed when the two polyelectrolyte layers span both the surfaces. Both the repulsive barrier and bridging attraction depend on the charge density of the polymer backbone and the bulk salt concentration but not on the chain length in the semidilute regime studied.

Maxwell Equations without a Polarization Field, Using a Paradigm from Biophysics

When forces are applied to matter, the distribution of mass changes. Similarly, when an electric field is applied to matter with charge, the distribution of charge changes. The change in the distribution of charge (when a local electric field is applied) might in general be called the induced charge. When the change in charge is simply related to the applied local electric field, the polarization field P is widely used to describe the induced charge. This approach does not allow electrical measurements (in themselves) to determine the structure of the polarization fields. Many polarization fields will produce the same electrical forces because only the divergence of polarization enters Maxwell's first equation, relating charge and electric forces and field. The curl of any function can be added to a polarization field P without changing the electric field at all. The divergence of the curl is always zero. Additional information is needed to specify the curl and thus the structure of the P field. When the structure of charge changes substantially with the local electric field, the induced charge is a nonlinear and time dependent function of the field and P is not a useful framework to describe either the electrical or structural basis-induced charge. In the nonlinear, time dependent case, models must describe the charge distribution and how it varies as the field changes. One class of models has been used widely in biophysics to describe field dependent charge, i.e., the phenomenon of nonlinear time dependent induced charge, called 'gating current' in the biophysical literature. The operational definition of gating current has worked well in biophysics for fifty years, where it has been found to makes neurons respond sensitively to voltage. Theoretical estimates of polarization computed with this definition fit experimental data. I propose that the operational definition of gating current be used to define voltage and time dependent induced charge, although other definitions may be needed as well, for example if the induced charge is fundamentally current dependent. Gating currents involve substantial changes in structure and so need to be computed from a combination of electrodynamics and mechanics because everything charged interacts with everything charged as well as most things mechanical. It may be useful to separate the classical polarization field as a component of the total induced charge, as it is in biophysics. When nothing is known about polarization, it is necessary to use an approximate representation of polarization with a dielectric constant that is a single real positive number. This approximation allows important results in some cases, e.g., design of integrated circuits in silicon semiconductors, but can be seriously misleading in other cases, e.g., ionic solutions.

Size and charge correlations in spherical electric double layers: a case study with fully asymmetric mixed electrolytes within the solvent primitive model

Size and charge correlations in spherical electric double layers are investigated through Monte Carlo simulations and density functional theory, through a solvent primitive model representation. A fully asymmetric mixed electrolyte is used for the small ions, whereas the solvent, apart from being a continuum dielectric, is also treated as an individual component. A partially perturbative density functional theory is adopted here, and for comparison, a standard canonical ensemble Monte Carlo simulation is used. The hard-sphere free energy is treated within a weighted density approach and the residual ionic contribution is estimated through perturbation around the uniform density. The results from both methods corroborate each other quantitatively over a wide range of physical parameters. The importance of structural correlations is envisaged through the size and charge asymmetry of the supporting electrolytes that includes the solvent as a component.

Size and charge correlations in spherical electric double layers are investigated through Monte Carlo simulations and density functional theory, through a solvent primitive model representation.

Charge regulation mechanism in end-tethered weak polyampholytes

Weak polyampholytes, containing oppositely charged dissociable groups, are expected to be responsive to changes in ionic conditions. Here, we determine structural and thermodynamic properties, including the charged groups' degrees of dissociation, of end-tethered weak polyampholyte layers as a function of salt concentration, pH, and the solvent quality. For diblock weak polyampholytes grafted by their acidic blocks, we find that the acidic monomers increase their charge while the basic monomers decrease their charge with decreasing salt concentration for pH values less than the pKa value of both monomers and vice versa when the pH > pKa. This complex charge regulation occurs because the electrostatic attraction between oppositely charged blocks is stronger than the repulsion between monomers with the same charge in both good and poor solvents when the screening by salt ions is weak. This is evidenced by the retraction of the top block into the bottom layer. In the case of poor solvent conditions to the basic block (the top block), we find lateral segregation of basic monomers into micelles, forming a two-dimensional hexagonal pattern on the surface at intermediate and high pH values for monovalent salt concentrations from 0.01 to 0.1 M. When the solvent is poor to both blocks, we find lateral segregation of the grafted acidic block into lamellae with longitudinal undulations of low and high acidic monomer density. By exploiting weak block polyampholytes, our work expands the parameter space for creating responsive surfaces stable over a wide range of pH and salt concentration.

Theoretical Modeling of Chemical Equilibrium in Weak Polyelectrolyte Layers on Curved Nanosystems

Surface functionalization with end-tethered weak polyelectrolytes (PE) is a versatile way to modify and control surface properties, given their ability to alter their degree of charge depending on external cues like pH and salt concentration. Weak PEs find usage in a wide range of applications, from colloidal stabilization, lubrication, adhesion, wetting to biomedical applications such as drug delivery and theranostics applications. They are also ubiquitous in many biological systems. Here, we present an overview of some of the main theoretical methods that we consider key in the field of weak PE at interfaces. Several applications involving engineered nanoparticles, synthetic and biological nanopores, as well as biological macromolecules are discussed to illustrate the salient features of systems involving weak PE near an interface or under (nano)confinement. The key feature is that by confining weak PEs near an interface the degree of charge is different from what would be expected in solution. This is the result of the strong coupling between structural organization of weak PE and its chemical state. The responsiveness of engineered and biological nanomaterials comprising weak PE combined with an adequate level of modeling can provide the keys to a rational design of smart nanosystems.

On the origin of oscillatory interactions between surfaces mediated by polyelectrolyte solution

We use a numerical implementation of polymer classical density functional theory with an incompressibility condition to study the system consisting of nonadsorbing polyelectrolytes confined by two planar surfaces and quantify the effective interaction between the two planar surfaces as a function of the polyelectrolyte and salt concentrations. Our results indicate that for the uncharged surfaces (and weakly charged surfaces), the effective interaction primarily consists of a short-range attraction due to the depletion followed by repulsion due to the electric double layer overlapping and electrostatic correlations. For salt-free and low salt concentration systems, the magnitude of the repulsion barrier is determined by the overlap between the electric double layers, while at relatively high salt concentrations, the magnitude of the repulsion barrier is determined by the electrostatic correlations. Due to the competition between the electric double layer and the electrostatic correlations, the magnitude of the repulsion barrier varies nonmonotonically. In contrast, a mean-field Poisson-Boltzmann treatment of the electrostatics predicts a monotonically decreasing repulsion barrier with increasing salt concentration. At moderate salt concentrations, our theory predicts oscillatory interaction profiles. A comparison with the mean-field Poisson-Boltzmann treatment of electrostatics suggests that the oscillations are due primarily to electrostatic correlations.

Complex Coacervation in Asymmetric Solutions of Polycation and Polyanion

We study coacervation upon mixing two aqueous solutions of polyelectrolytes (PEs) with opposite charge, by considering asymmetric effects of PE composition and charge valency. The phase behavior, interfacial structure, and coacervate composition are investigated by a classical density-functional theory. We find two types of coacervation that are different in their density. Supernatant phase in low-density coacervation (LDCA) fully consists of small ions, while in high-density coacervation (HDCA) it contains a considerable amount of PE chains. Asymmetric PE composition could generate an electric double layer at the interface of coacervate. For HDCA, ion density changes monotonically, while for LDCA it shows a global minimum at the double layer, giving a low tension value. Charged species of high charge valency enhance the existence of double layer. Our results explained the coacervate structure of low interfacial tension, which is important for experiments and industrial applications.

Dynamic Adsorption of Ions into Like-Charged Nanospace: A Dynamic Density Functional Theory Study

The adsorption processes of ions into charged nanospace are associated with many practical applications. Whereas a large number of microporous materials have been prepared toward efficient adsorption of ions from solutions, theoretical models that allow for capturing the characteristics of ion dynamic adsorption into like-charged nanopores are still few. The difficulty originates from the overlapping of electric potentials inside the pores. Herein, a theoretical model is proposed by incorporating dynamic density functional theory with modified Poisson equation for investigating the dynamic adsorption of ions into like-charged nanoslits. This model is rationalized by comparing the theoretical predictions with corresponding simulation results. Afterward, by analyzing the adsorption dynamics, we show that the overlapping effect is associated with the pore size, ion bulk concentration, and surface charge density, and it plays a dominant role in the coupling between the total adsorption amount of ions and total adsorption time. Specifically, with weak overlapping effect, the total adsorption amount is intuitively proportional to the total adsorption time; however, when the overlapping effect is strong, the total adsorption amount may be inversely proportional to the total adsorption time, indicating that both high adsorption amount and short adsorption time can be achieved simultaneously. This work provides a meaningful insight toward the rational design and optimization of microporous materials for efficient ion adsorption.

Density functional approach to the description of the structure of dimer nanoparticles at liquid–liquid interfaces

A density functional approach to the description of the structure of dimer nanoparticles at liquid–liquid interfaces.