Weak and strong coupling theories for polarizable colloids and nanoparticles

Criteria of Distribution Transitions in Dispersed Multiphase Systems Based on an Extended Lattice Model

Dispersed multiphase systems are ubiquitous in biological systems, energy industries, and medical science. The distribution transition of the dispersed phase is critical to the properties and functions of the multiphase systems, among which the agglomeration, adsorption, and extraction processes are of most significance due to their impact on the colloidal stability, interface modification, and particle synthesis. To reveal fundamental correlations between the macroscopic particle distributions and the microscopic interactions, general thermodynamic models of the dispersed multiphase systems are needed. Here, based on Meyer's model, which is restricted to binary isotropic mixtures, we propose a novel extended lattice model that can be applied to multicomponent anisotropic mixtures with interfaces considered. For agglomeration, adsorption, and extraction processes, the approximate free energy differences between the initial distribution and the final distribution are obtained. Based on the minimum free energy principle, the above free energy differences are used to derive three criteria for the prediction of the preferable distribution of the system with given interparticle interaction potentials. While the quasi-uniform number density assumption is still required for all the previous lattice models, the long-range interactions neglected by previous lattice models are incorporated. The validity of our model and criteria is verified by many-body dissipative particle dynamics (mDPD) simulations. By tuning the interaction coefficients between mDPD particles, the simulated distribution transitions for all the agglomeration, adsorption, and extraction cases perfectly match the predictions from the three criteria. The good agreement between the theoretical predictions and the mDPD simulation results shows the great potential of our model for applications in various dispersed multiphase systems.

Effects of electrostatic coupling and surface polarization on polyelectrolyte brush structure

In this work, we perform molecular dynamics simulations to study a spherical polyelectrolyte brush. We explore the effects of surface polarization and electrostatic coupling on brush size and distribution of counterions. The method of image charges is considered to take into account surface polarization, considering a metallic, an unpolarizable, and a dielectric nano-core. It is observed that, for all cases, a moderate shrinking-swelling effect appears with an increase in the electrostatic coupling parameter. This effect occurs under high Manning ratios. The curves relating the average size of polyelectrolyte brush as a function of coupling show a minimum. The results show that the grafting density of polyelectrolytes on the nano-core surface plays an important role in the polarization effect. We consider a modified Poisson-Boltzmann theory to describe the counterion profiles around the brush in the case of unpolarizable nano-cores and weak electrostatic coupling.

Like-Charge Attraction at the Nanoscale: Ground-State Correlations and Water Destructuring

Like-charge attraction, driven by ionic correlations, challenges our understanding of electrostatics both in soft and hard matter. For two charged planar surfaces confining counterions and water, we prove that, even at relatively low correlation strength, the relevant physics is the ground-state one, oblivious of fluctuations. Based on this, we derive a simple and accurate interaction pressure that fulfills known exact requirements and can be used as an effective potential. We test this equation against implicit-solvent Monte Carlo simulations and against explicit-solvent simulations of cement and several types of clays. We argue that water destructuring under nanometric confinement drastically reduces dielectric screening, enhancing ionic correlations. Our equation of state at reduced permittivity therefore explains the exotic attractive regime reported for these materials, even in the absence of multivalent counterions.

Reactive Monte Carlo simulations for charge regulation of colloidal particles

We use a reactive Monte Carlo simulation method and the primitive model of electrolyte to study acid-base equilibrium that controls charge regulation in colloidal systems. The simulations are performed in a semi-grand canonical ensemble in which colloidal suspension is in contact with a reservoir of salt and strong acid. The interior of colloidal particles is modeled as a low dielectric medium, different from the surrounding water. The effective colloidal charge is calculated for different numbers of surface acidic groups, pH, salt concentrations, and types of electrolyte. In the case of potassium chloride, the titration curves are compared with the experimental measurements obtained using potentiometric titration. A good agreement is found between simulations and experiments. In the case of lithium chloride, the specific ionic adsorption is taken into account through the partial dehydration of lithium ion.

Electrolytes in regimes of strong confinement: surface charge modulations, osmotic equilibrium and electroneutrality

In the present work, we study an electrolyte solution confined between planar surfaces with nanopatterned charged domains, which has been connected to a bulk ionic reservoir. The system is investigated through an improved Monte Carlo (MC) simulation method, suitable for simulation of electrolytes in the presence of modulated surface charge distributions. We also employ a linear approach in the spirit of the classical Debye-Hückel approximation, which allows one to obtain explicit expressions for the averaged potentials, ionic profiles, effective surface interactions and the net ionic charge confined between the walls. Emphasis is placed on the limit of strongly confined electrolytes, in which case local electroneutrality in the inter-surface space might not be fulfilled. In order to access the effects of such a lack of local charge neutrality on the ion-induced interactions between surfaces with modulated charge domains, we consider two distinct model systems for the confined electrolyte: one in which a salt reservoir is explicitly taken into account via the osmotic equilibrium with an electrolyte of fixed bulk concentration, and a second one in which the equilibrium with a charge neutral ionic reservoir is implicitly considered. While in the former case the osmotic ion exchange might lead to non-vanishing net charges, in the latter model charge neutrality is enforced through the appearance of an implicit Donnan potential across the charged interfaces. A strong dependence of the ion-induced surface interactions on the employed model system is observed at all surface separations. These findings strongly suggest that due care is to be taken while choosing among different scenarios to describe the ion exchange in electrolytes confined between charged surfaces, even in cases when the monopole (non zero net charge) surface contributions are absent.

Enhancement of the ion concentration in a salt solution near a wall due to electrical image potentials and enhancement of surface tension due to the presence of salt

Effects of electrical image potentials on the salt ion concentration near a solid wall are studied using a one-loop approximation treatment of the grand canonical partition function, which is the Debye-Hückel approximation. Electrical image potentials resulting from both metallic and dielectric walls of dielectric constant larger than that of water near the wall are considered. Our treatment of this problem supports the conclusions of an earlier publication by one of the authors which shows that near a solid wall there should be a high concentration of ions, resulting from image potentials. We have also applied our treatment to the increase of the surface tension of a liquid that occurs when salt is dissolved in the liquid. Our treatment gives the -c_{s}log(c_{s}) dependence of the surface tension found by Onsager and Samarasa in the small c_{s} limit, where c_{s} is the salt concentration.

Strong-coupling theory of counterions with hard cores between symmetrically charged walls

By a combination of Monte Carlo simulations and analytical calculations, we investigate the effective interactions between highly charged planar interfaces, neutralized by mobile counterions (salt-free system). While most previous analysis have focused on pointlike counterions, we treat them as charged hard spheres. We thus work out the fate of like-charge attraction when steric effects are at work. The analytical approach partitions counterions in two subpopulations, one for each plate, and integrates out one subpopulation to derive an effective Hamiltonian for the remaining one. The effective Hamiltonian features plaquette four-particle interactions, and it is worked out by computing a Gibbs-Bogoliubov inequality for the free energy. At the root of the treatment is the fact that under strong electrostatic coupling, the system of charges forms an ordered arrangement, that can be affected by steric interactions. Fluctuations around the reference positions are accounted for. To dominant order at high coupling, it is found that steric effects do not significantly affect the interplate effective pressure, apart at small distances where hard-sphere overlap are unavoidable, and thus rule out configurations.

Effects of electrical image forces on salt dissolved in water

It will be shown that for a solution of salt dissolved in water in contact with a metallic or dielectric wall, the concentration of salt ions (both positive and negative) near the wall can be large enough to exceed the salt's solubility limit, as a result of electrical image charge forces. In addition, since the dielectric constant of water increases from 2.1 at the wall to 81 at about 1 nm from a solid wall, there will be an attractive image potential near the plane on which this increase of the dielectric constant occurs. Such large ion concentrations could possibly be used in water desalination.

Charge polarization near dielectric interfaces and the multiple-scattering formalism

Interfacial charge polarization is ubiquitous in systems with sharp dielectric contrast. Fully resolving the interfacial charges often relies on demanding numerical algorithms to solve the boundary value problem. The recent development of an analytical multiple-scattering formalism to solve the interfacial charge polarization problem for particles carrying monopolar, dipolar, and multipolar charges is reviewed. Every term produced in this formalism has a simple interpretation, and terms for spherical particles can be rapidly evaluated using an image-line construction. Several practical applications of this formalism are illustrated. A dielectric virial expansion for polarizable particles based on this formalism is also described. The origins of singular polarization charges for particles in close contact are explained and evaluated for both dielectric and conducting spheres.

An improved theory of the electric conductance of ionic solutions based on the concept of the ion-atmosphere's smaller-ion shell

Although the electric conductivity of ionic solutions is of great importance in science and technology, complete understanding of the physical factors and principles underlying ion and charge transport in solution has so far been missing. The Debye-Hückel-Onsager (DHO) theory of the equivalent ionic conductivity (λi for single ion i, Λ for a binary electrolyte) is a limiting law of only a few electrolyte valence families, and its extension to solution of finite concentration has not been entirely successful; at most, it has led to semiempirical expressions of λi or Λ vs. concentration (c), incorporating adjustable parameters of no clear physical meaning. Improving the description of the ion atmosphere and including in the conductivity model the smaller-ion shell (SiS) of this atmosphere (D. Fraenkel, Mol. Phys., 2010, 108, 1435) allow now an extension of the DHO theory to a "DHO-SiS" theory that offers analytic expressions for λi and Λ as a function of c without adjustable parameters. The fit with experiment is very good, often excellent, for many electrolytes of the various valence families. λi values can be converted to single-ion activity coefficients, γis, suggesting that the latter ionic factors do have physical validity, contrary to what many literature sources claim. Together with the recently proposed theory of the limiting conductivities, λ0i and Λ0 (D. Fraenkel, Mol. Phys., 2018, 116, 2271), the new theory offers a simple and comprehensive ion conductivity treatment that improves upon existing theoretical treatments, with λi and Λ expressed by ionic sizes instead of λ0i and Λ0.

Simulations of Coulomb systems confined by polarizable surfaces using periodic Green functions

We present an efficient approach for simulating Coulomb systems confined by planar polarizable surfaces. The method is based on the solution of the Poisson equation using periodic Green functions. It is shown that the electrostatic energy arising from the surface polarization can be decoupled from the energy due to the direct Coulomb interaction between the ions. This allows us to combine an efficient Ewald summation method, or any other fast method for summing over the replicas, with the polarization contribution calculated using Green function techniques. We apply the method to calculate density profiles of ions confined between the charged dielectric and metal surfaces.

Ionic size effects on the Poisson-Boltzmann theory

In this paper, we develop a simple theory to study the effects of ionic size on ionic distributions around a charged spherical particle. We include a correction to the regular Poisson-Boltzmann equation in order to take into account the size of ions in a mean-field regime. The results are compared with Monte Carlo simulations and a density functional theory based on the fundamental measure approach and a second-order bulk expansion which accounts for electrostatic correlations. The agreement is very good even for multivalent ions. Our results show that the theory can be applied with very good accuracy in the description of ions with highly effective ionic radii and low concentration, interacting with a colloid or a nanoparticle in an electrolyte solution.

Ionic structure around polarizable metal nanoparticles in aqueous electrolytes

Metal nanoparticles are receiving increased scientific attention owing to their unique physical and chemical properties that make them suitable for a wide range of applications in diverse fields, such as electrochemistry, biochemistry, and nanomedicine. Their high metallic polarizability is a crucial determinant that defines their electrostatic character in various electrolyte solutions. Here, we introduce a continuum-based model of a metal nanoparticle with explicit polarizability in the presence of different kinds of electrolytes. We employ several, variously sophisticated, theoretical approaches, corroborated by Monte Carlo simulations in order to elucidate the basic electrostatics principles of the model. We investigate how different kinds of asymmetries between the ions result in non-trivial phenomena, such as charge separation and a build-up of a so-called zero surface-charge double layer.

Double-layer force suppression between charged microspheres

In this paper we propose a protocol to suppress double-layer forces between two microspheres immersed in a dielectric medium, being one microsphere metallic at a controlled potential ψ_{M} and the other a charged one either metallic or dielectric. The approach is valid for a wide range of distances between them. We show that, for a given distance between the two microspheres, the double-layer force can be totally suppressed by simply tuning ψ_{M} up to values dictated by the linearized Poisson-Boltzmann equation. Our key finding is that such values can be substantially different from the ones predicted by the commonly used proximity force approximation, also known as the Derjaguin approximation, even in situations where the latter is expected to be accurate. The proposed procedure can be used to suppress the double-layer interaction in force spectroscopy experiments, thus paving the way for measurements of other surface interactions, such as Casimir dispersion forces.

On the physics of both surface overcharging and charge reversal at heterophase interfaces

A series of Monte Carlo simulations are employed to reveal the physics of both surface overcharging and charge reversal at a negatively charged dielectric interface exposed to a bulk solution containing a +2:−1 electrolyte in the absence and presence of a monovalent salt.

Mean-field beyond mean-field: the single particle view for moderately to strongly coupled charged fluids

In a counter-ion only charged fluid, Coulomb coupling is quantified by a single dimensionless parameter. Yet, the theoretical treatment of moderately to strongly coupled charged fluids is a difficult task, central to the understanding of a wealth of soft matter problems, including biological systems. We show that the corresponding coupling regime can be remarkably well described by a single particle treatment, which, at variance with previous works, takes due account of inter-ionic interactions. To this end, the prototypical problem of a planar charged dielectric interface is worked out. Testing our predictions against Monte Carlo simulation data reveals an excellent agreement.

Inhomogeneous screening near the dielectric interface

Screening is one of the most important concepts in the study of charged systems. Near a dielectric interface, the ion distribution in a salt solution can be highly nonuniform. Here, we develop a theory that self-consistently treats the inhomogeneous screening effects. At higher concentrations when the bulk Debye screening length is comparable to the Bjerrum length, the double layer structure and interfacial properties are significantly affected by the inhomogeneous screening. In particular, the depletion zone is considerably wider than that predicted by the bulk screening approximation or the WKB approximation. The characteristic length of the depletion layer in this regime scales with the Bjerrum length, resulting in a linear increase of the negative adsorption of ions with concentration, in agreement with experiments. For asymmetric salts, inhomogeneous screening leads to enhanced charge separation and surface potential.

Charge reversal at a planar boundary between two dielectrics

Despite the ubiquitous character and relevance of the electric double layer in the entire realm of interface and colloid science, very little is known of the effect that surface heterogeneity exerts on the underlying mechanisms of ion adsorption. Herein, computer simulations offer a perspective that, in sharp contrast to the homogeneously charged surface, discrete groups promote multivalent counterion binding, leading to charge reversal but possibly having not a sign change of the electrophoretic mobility. Counterintuitively, the introduction of dielectric images yields a significantly greater accumulation of counterions, which further facilitates the magnitude of charge reversal. The reported results are very sensitive to both the degree of ion hydration and the representation of surface charges. Our findings shed light on the mechanism for charge reversal over a broad range of coupling regimes operating the adsorption of counterions through surface group bridging attraction with their own images and provide opportunities for experimental studies and theoretical development.

Electrolytes between dielectric charged surfaces: Simulations and theory

We present a simulation method to study electrolyte solutions in a dielectric slab geometry using a modified 3D Ewald summation. The method is fast and easy to implement, allowing us to rapidly resum an infinite series of image charges. In the weak coupling limit, we also develop a mean-field theory which allows us to predict the ionic distribution between the dielectric charged plates. The agreement between both approaches, theoretical and simulational, is very good, validating both methods. Examples of ionic density profiles in the strong electrostatic coupling limit are also presented. Finally, we explore the confinement of charge asymmetric electrolytes between neutral surfaces.

On the theoretical description of weakly charged surfaces

It is widely accepted that the Poisson-Boltzmann (PB) theory provides a valid description for charged surfaces in the so-called weak coupling limit. Here, we show that the image charge repulsion creates a depletion boundary layer that cannot be captured by a regular perturbation approach. The correct weak-coupling theory must include the self-energy of the ion due to the image charge interaction. The image force qualitatively alters the double layer structure and properties, and gives rise to many non-PB effects, such as nonmonotonic dependence of the surface energy on concentration and charge inversion. In the presence of dielectric discontinuity, there is no limiting condition for which the PB theory is valid.

Are antagonistic salts surfactants?

It is well known that surfactants decrease both water/air and water/oil interfacial tensions whereas in contrast inorganic salts increase both. We study a new, third class of surface-active ionic solutes, which have been called antagonistic salts, consisting of an organic group with a small inorganic counterion. These show decreased interfacial tension at the oil/water interface due to a redistribution of the organic group in the oil but do not show any surface activity at the air/water interface and are consequently different from surfactants that lower both tensions. We use a simple modeling using Poisson-Boltzmann theory that accounts for the surface activity of the antagonistic salt at the water/oil interface.

A simulation assessment of the thermodynamics of dense ion-dipole mixtures with polarization

Molecular dynamics (MD) simulations are employed to ascertain the relative importance of various electrostatic interaction contributions, including induction interactions, to the thermodynamics of dense, hot ion-dipole mixtures. In the absence of polarization, we find that an MD-constrained free energy term accounting for the ion-dipole interactions, combined with well tested ionic and dipolar contributions, yields a simple, fairly accurate free energy form that may be a better option for describing the thermodynamics of such mixtures than the mean spherical approximation (MSA). Polarization contributions induced by the presence of permanent dipoles and ions are found to be additive to a good approximation, simplifying the thermodynamic modeling. We suggest simple free energy corrections that account for these two effects, based in part on standard perturbative treatments and partly on comparisons with MD simulation. Even though the proposed approximations likely need further study, they provide a first quantitative assessment of polarization contributions at high densities and temperatures and may serve as a guide for future modeling efforts.

Self-energy-modified Poisson-Nernst-Planck equations: WKB approximation and finite-difference approaches

We propose a modified Poisson-Nernst-Planck (PNP) model to investigate charge transport in electrolytes of inhomogeneous dielectric environment. The model includes the ionic polarization due to the dielectric inhomogeneity and the ion-ion correlation. This is achieved by the self energy of test ions through solving a generalized Debye-Hückel (DH) equation. We develop numerical methods for the system composed of the PNP and DH equations. Particularly, toward the numerical challenge of solving the high-dimensional DH equation, we developed an analytical WKB approximation and a numerical approach based on the selective inversion of sparse matrices. The model and numerical methods are validated by simulating the charge diffusion in electrolytes between two electrodes, for which effects of dielectrics and correlation are investigated by comparing the results with the prediction by the classical PNP theory. We find that, at the length scale of the interface separation comparable to the Bjerrum length, the results of the modified equations are significantly different from the classical PNP predictions mostly due to the dielectric effect. It is also shown that when the ion self energy is in weak or mediate strength, the WKB approximation presents a high accuracy, compared to precise finite-difference results.

Correlation potential of a test ion near a strongly charged plate

We analytically calculate the correlation potential of a test ion near a strongly charged plate inside a dilute m:-n electrolyte. We do this by calculating the electrostatic Green's function in the presence of a nonlinear background potential, the latter having been obtained using the nonlinear Poisson-Boltzmann theory. We consider the general case where the dielectric constants of the plate and the electrolyte are distinct. The following generic results emerge from our analyses: (1) If the distance to the plate Δz is much larger than a Gouy-Chapman length, the plate surface will behave effectively as an infinitely charged surface, and the dielectric constant of the plate effectively plays no role. (2) If Δz is larger than a Gouy-Chapman length but shorter than a Debye length, the correlation potential can be interpreted in terms of an image charge that is three times larger than the source charge. This behavior is independent of the valences of the ions. (3) The Green's function vanishes inside the plate if the surface charge density is infinitely large; hence the electrostatic potential is constant there. In this respect, a strongly charged plate behaves like a conductor plate. (4) If Δz is smaller than a Gouy-Chapman length, the correlation potential is dominated by the conventional image charge due to the dielectric discontinuity at the interface. (5) If Δz is larger than a Debye length, the leading order behavior of the correlation potential will depend on the valences of the ions in the electrolyte. Furthermore, inside an asymmetric electrolyte, the correlation potential is singly screened, i.e., it undergoes exponential decay with a decay width equal to the Debye length.

Large counterions boost the solubility and renormalized charge of suspended nanoparticles

Colloidal particles are ubiquitous in biology and in everyday products such as milk, cosmetics, lubricants, paints, or drugs. The stability and aggregation of colloidal suspensions are of paramount importance in nature and in diverse nanotechnological applications, including the fabrication of photonic materials and scaffolds for biological assemblies, gene therapy, diagnostics, targeted drug delivery, and molecular labeling. Electrolyte solutions have been extensively used to stabilize and direct the assembly of colloidal particles. In electrolytes, the effective electrostatic interactions among the suspended colloids can be changed over various length scales by tuning the ionic concentration. However, a major limitation is gelation or flocculation at high salt concentrations. This is explained by classical theories, which show that the electrostatic repulsion among charged colloids is significantly reduced at high electrolyte concentrations. As a result, these screened colloidal particles are expected to aggregate due to short-range attractive interactions or dispersion forces as the salt concentration increases. We discuss here a robust, tunable mechanism for colloidal stability by which large counterions prevent highly charged nanoparticles from aggregating in salt solutions with concentrations up to 1 M. Large counterions are shown to generate a thicker ionic cloud in the proximity of each charged colloid, which strengthens short-range repulsions among colloidal particles and also increases the corresponding renormalized colloidal charge perceived at larger separation distances. These effects thus provide a reliable stabilization mechanism in a broad range of biological and synthetic colloidal suspensions.

Electrostatic interactions mediated by polarizable counterions: weak and strong coupling limits

We investigate the statistical mechanics of an inhomogeneous Coulomb fluid composed of charged particles with static polarizability. We derive the weak- and the strong-coupling approximations and evaluate the partition function in a planar dielectric slab geometry with charged boundaries. We investigate the density profiles and the disjoining pressure for both approximations. Comparison to the case of non-polarizable counterions shows that polarizability brings important differences in the counterion density distribution as well as the counterion mediated electrostatic interactions between charged dielectric interfaces.

Surface tension of an electrolyte-air interface: a Monte Carlo study

We present a new method for calculating the surface tension of an electrolyte-air interface using Monte Carlo (MC) simulations with an implicit solvent in a spherical drop geometry. The boundary conditions for the electric field at the interface are accounted for using image and counter-image charges. The density profiles obtained from the simulations are used to calculate the excess surface tension of the electrolyte-air interface using the Gibbs adsorption isotherm equation. The results are found to be in good agreement with experiments and the earlier theoretical calculations.

Computational evidence of two driving mechanisms for overcharging in an electric double layer near the point of zero charge

We have adopted an ensemble Monte Carlo simulation method to systematically verify two physical driving mechanisms responsible for overcharging which refers to the adsorption of an effective charge onto a like-charged planar surface around the point of zero charge within the primitive model of mixed electrolytes with varying salt concentrations. One is electrostatic in character dominated by dielectric images and the other is purely entropic in origin by ionic size asymmetry effects, of which the former has never been reported both theoretically and experimentally and the latter could be interpreted satisfactorily in terms of available theoretical approaches. The electrostatically driven mechanism is found to critically depend on the ionic sizes while the entropically driven mechanism occurs with almost the same efficiency in a relative wide range of surface charge density. Depending on the delicate interplay between charge and steric correlations, the two distinct driving mechanisms may cooperatively give rise to a more pronounced overcharging process.

Models of membrane electrostatics

Formulas are derived for the electrostatic potential of a charge in or near a membrane modeled as one or more dielectric slabs lying between two semi-infinite dielectrics. One can use these formulas in Monte Carlo codes to compute the distribution of ions near cell membranes more accurately than by using Poisson-Boltzmann theory or its linearized version. Here I use them to discuss the electric field of a uniformly charged membrane, the image charges of an ion, the distribution of salt ions near a charged membrane, the energy of a zwitterion near a lipid slab, and the effect of including the phosphate head groups as thin layers of high electric permittivity.

Macroion solutions in the cell model studied by field theory and Monte Carlo simulations

Aqueous solutions of charged spherical macroions with variable dielectric permittivity and their associated counterions are examined within the cell model using a field theory and Monte Carlo simulations. The field theory is based on separation of fields into short- and long-wavelength terms, which are subjected to different statistical-mechanical treatments. The simulations were performed by using a new, accurate, and fast algorithm for numerical evaluation of the electrostatic polarization interaction. The field theory provides counterion distributions outside a macroion in good agreement with the simulation results over the full range from weak to strong electrostatic coupling. A low-dielectric macroion leads to a displacement of the counterions away from the macroion.

Solubility and transport of cationic and anionic patterned nanoparticles

We analyze bulk diffusion and transport through hydrophobic nanochannels of nanoparticles (NPs) with different hydrophobic-hydrophilic patterns achieved by coating a fraction of the NP sites with positive or negative charges via explicit solvent molecular dynamics simulations. Ten different charge pattern types including Janus charged-hydrophobic NPs are studied. The cationic NPs are more affected by the patterns and have higher diffusion constants and fluxes than their anionic NPs counterparts. The NP-water interaction dependence on surface pattern and field strength explains these observations. The NP-water Coulomb interaction of anionic NPs in the bulk, which are much stronger than the hydrophobic NP-water interactions, are stronger for NPs with higher localized charge, and stronger than in the cationic NPs counterparts. The diffusion and transport of anionic NPs such as proteins and protein charge ladders with the same total charge but different surface charge patterns are slowest for the highest localized charge pattern, which also adsorb strongest onto surfaces. Our model demonstrates the separation (by reverse osmosis, capillary electrophoresis, or chromatography) of cationic NPs, including proteins with equal net charge but different surface charge distributions.