Interfacial solvation can explain attraction between like-charged objects in aqueous solution

A charge-dependent long-ranged force drives tailored assembly of matter in solution

The interaction between charged objects in solution is generally expected to recapitulate two central principles of electromagnetics: (1) like-charged objects repel, and (2) they do so regardless of the sign of their electrical charge. Here we demonstrate experimentally that the solvent plays a hitherto unforeseen but crucial role in interparticle interactions, and importantly, that interactions in the fluid phase can break charge-reversal symmetry. We show that in aqueous solution, negatively charged particles can attract at long range while positively charged particles repel. In solvents that exhibit an inversion of the net molecular dipole at an interface, such as alcohols, we find that the converse can be true: positively charged particles may attract whereas negatives repel. The observations hold across a wide variety of surface chemistries: from inorganic silica and polymeric particles to polyelectrolyte- and polypeptide-coated surfaces in aqueous solution. A theory of interparticle interactions that invokes solvent structuring at an interface captures the observations. Our study establishes a nanoscopic interfacial mechanism by which solvent molecules may give rise to a strong and long-ranged force in solution, with immediate ramifications for a range of particulate and molecular processes across length scales such as self-assembly, gelation and crystallization, biomolecular condensation, coacervation, and phase segregation.

Electrostatic interactions between charge regulated spherical macroions

We study the interaction between two charge regulating spherical macroions with dielectric interior and dissociable surface groups immersed in a monovalent electrolyte solution. The charge dissociation is modelled via the Frumkin-Fowler-Guggenheim isotherm, which allows for multiple adsorption equilibrium states. The interactions are derived from the solutions of the mean-field Poisson-Boltzmann type theory with charge regulation boundary conditions. For a range of conditions we find symmetry breaking transitions from symmetric to asymmetric charge distribution exhibiting annealed charge patchiness, which results in like-charge attraction even in a univalent electrolyte-thus fundamentally modifying the nature of electrostatic interactions in charge-stabilized colloidal suspensions.

Electrokinetic, electrochemical, and electrostatic surface potentials of the pristine water liquid-vapor interface

Although conceptually simple, the air-water interface displays rich behavior and is subject to intense experimental and theoretical investigations. Different definitions of the electrostatic surface potential as well as different calculation methods, each relevant for distinct experimental scenarios, lead to widely varying potential magnitudes and sometimes even different signs. Based on quantum-chemical density-functional-theory molecular dynamics (DFT-MD) simulations, different surface potentials are evaluated and compared to force-field (FF) MD simulations. As well explained in the literature, the laterally averaged electrostatic surface potential, accessible to electron holography, is dominated by the trace of the water molecular quadrupole moment, and using DFT-MD amounts to +4.35 V inside the water phase, very different from results obtained with FF water models which yield negative values of the order of -0.4 to -0.6 V. Thus, when predicting potentials within water molecules, as relevant for photoelectron spectroscopy and non-linear interface-specific spectroscopy, DFT simulations should be used. The electrochemical surface potential, relevant for ion transfer reactions and ion surface adsorption, is much smaller, less than 200 mV in magnitude, and depends specifically on the ion radius. Charge transfer between interfacial water molecules leads to a sizable surface potential as well. However, when probing electrokinetics by explicitly applying a lateral electric field in DFT-MD simulations, the electrokinetic ζ-potential turns out to be negligible, in agreement with predictions using continuous hydrodynamic models. Thus, interfacial polarization charges from intermolecular charge transfer do not lead to significant electrokinetic mobility at the pristine vapor-liquid water interface, even assuming these transfer charges are mobile in an external electric field.

Single-Molecule Trapping and Measurement in a Nanostructured Lipid Bilayer System

The repulsive electrostatic force between a biomolecule and a like-charged surface can be geometrically tailored to create spatial traps for charged molecules in solution. Using a parallel-plate system composed of silicon dioxide surfaces, we recently demonstrated single-molecule trapping and high precision molecular charge measurements in a nanostructured free energy landscape. Here we show that surfaces coated with charged lipid bilayers provide a system with tunable surface properties for molecular electrometry experiments. Working with molecular species whose effective charge and geometry are well-defined, we demonstrate the ability to quantitatively probe the electrical charge density of a supported lipid bilayer. Our findings indicate that the fraction of charged lipids in nanoslit lipid bilayers can be significantly different from that in the precursor lipid mixtures used to generate them. We also explore the temporal stability of bilayer properties in nanofluidic systems. Beyond their relevance in molecular measurement, such experimental systems offer the opportunity to examine lipid bilayer formation and wetting dynamics on nanostructured surfaces.

Far-Field Electrostatic Signatures of Macromolecular 3D Conformation

In solution as in vacuum, the electrostatic field distribution in the vicinity of a charged object carries information on its three-dimensional geometry. We report on an experimental study exploring the effect of molecular shape on long-range electrostatic interactions in solution. Working with DNA nanostructures carrying approximately equal amounts of total charge but each in a different three-dimensional conformation, we demonstrate that the geometry of the distribution of charge in a molecule has substantial impact on its electrical interactions. For instance, a tetrahedral structure, which is the most compact distribution of charge we tested, can create a far-field effect that is effectively identical to that of a rod-shaped molecule carrying half the amount of total structural charge. Our experiments demonstrate that escape-time electrometry (ETe) furnishes a rapid and facile method to screen and identify 3D conformations of charged biomolecules or molecular complexes in solution.

The Role of Surface Chemistry in the Orientational Behavior of Water at an Interface

Molecular dynamics studies have demonstrated that molecular water at an interface, with either a gas or a solid, displays anisotropic orientational behavior in contrast to its bulk counterpart. This effect has been recently implicated in the like-charge attraction problem for colloidal particles in solution. Here, negatively charged particles in solution display a long-ranged attraction where continuum electrostatic theory predicts monotonically repulsive interactions, particularly in solutions with monovalent salt ions at low ionic strength. Anisotropic orientational behavior of solvent molecules at an interface gives rise to an excess interfacial electrical potential which we suggest generates an additional solvation contribution to the total free energy that is traditionally overlooked in continuum descriptions of interparticle interactions in solution. In the present investigation we perform molecular dynamics simulation based calculations of the interfacial potential using realistic surface models representing various chemistries as well as different solvents. Similar to previous work that focused on simple model surfaces constructed by using oxygen atoms, we find that solvents at more realistic model surfaces exhibit substantial anisotropic orientational behavior. We explore the dependence of the interfacial solvation potential on surface properties such as surface group chemistry and group density at silica and carboxylated polystyrene interfaces. For water, we note surprisingly good agreement between results obtained for a simple O-atom wall and more complex surface models, suggesting a general qualitative consistency of the interfacial solvation effect for surfaces in contact with water. In contrast, for an aprotic solvent such as DMSO, surface chemistry appears to exert a stronger influence on the sign and magnitude of the interfacial solvation potential. The study carries broad implications for molecular-scale interactions and may find relevance in explaining a range of phenomena in soft-matter physics and cell biology.

Nonmonotonic Pair Potentials in the Interaction of Like-Charged Objects in Solution

We consider the long-standing like-charge attraction problem, wherein under certain conditions, similarly charged spheres suspended in aqueous electrolyte have been observed to display a minimum in their interaction potential, contrary to the intuitively expected monotonically varying repulsion. Recently, we described an interfacial mechanism invoking the molecular nature of the solvent that explains this anomalous experimental observation. In our model for the interaction of negatively charged particles in water, the minimum in the pair potential results from the superposition of competing contributions to the total free energy. One of these contributions is the canonical repulsive electrostatic term, whereas the other is a solvation-induced attractive contribution. We find that whereas both contributions grow approximately exponentially with decreasing interparticle separation, the occurrence of a stable, long-ranged minimum in the pair potential arises from differences in the precise interparticle separation dependence of the two terms. Specifically, the interfacial solvation term exhibits a more gradual decay with distance than the electrostatic repulsion, permitting the attractive contribution to dominate the interaction at large distances. Importantly, these disparities become evident in quantities calculated from exact numerical solutions to the governing nonlinear Poisson-Boltzmann (PB) equation for the spatial electrical potential distribution in the system. In marked contrast, we find that the linearized PB equation, applicable in the regime of low surface electrical potentials, does not support nonmonotonic trends in the total interaction free energy within the present model. Our results point to the importance of exact descriptions of electrostatic interactions in real systems that most often do not subscribe to particular mathematical limits where analytical approximations may provide a sufficiently accurate description of the problem.

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.

Systematic optimization of a fragment-based force field against experimental pure-liquid properties considering large compound families: application to oxygen and nitrogen compounds

The CombiFF approach is a workflow for the automated refinement of force-field parameters against experimental condensed-phase data, considering entire classes of organic molecules constructed using a fragment library via combinatorial isomer enumeration. One peculiarity of this approach is that it relies on an electronegativity-equalization scheme to account for induction effects within molecules, with values of the atomic hardness and electronegativity as electrostatic parameters, rather than the partial charges themselves. In a previous article [M. P. Oliveira, M. Andrey, S. R. Rieder, L. Kern, D. F. Hahn, S. Riniker, B. A. C. Horta and P. H. Hünenberger, J. Chem. Theory. Comput. 2020, 16, 7525], CombiFF was introduced and applied to calibrate a GROMOS-compatible united-atom force field for the saturated acyclic (halo-)alkane family. Here, this scheme is employed for the construction of a corresponding force field for saturated acyclic compounds encompassing eight common chemical functional groups involving oxygen and/or nitrogen atoms, namely: ether, aldehyde, ketone, ester, alcohol, carboxylic acid, amine, and amide. Monofunctional as well as homo-polyfunctional compounds are considered. A total of 1712 experimental liquid densities ρliq and vaporization enthalpies ΔHvap concerning 1175 molecules are used for the calibration (339 molecules) and validation (836 molecules) of the 102 non-bonded interaction parameters of the force field. Using initial parameter values based on the GROMOS 2016H66 parameter set, convergence is reached after five iterations. Given access to one processor per simulated system, this operation only requires a few days of wall-clock computing time. After optimization, the root-mean-square deviations from experiment are 29.9 (22.4) kg m-3 for ρliq and 4.1 (5.5) kJ mol-1 for ΔHvap for the calibration (validation) set. Thus, a very good level of agreement with experiment is achieved in terms of these two properties, although the errors are inhomogeneously distributed across the different chemical functional groups.

Local structure of DNA toroids reveals curvature-dependent intermolecular forces

In viruses and cells, DNA is closely packed and tightly curved thanks to polyvalent cations inducing an effective attraction between its negatively charged filaments. Our understanding of this effective attraction remains very incomplete, partly because experimental data is limited to bulk measurements on large samples of mostly uncurved DNA helices. Here we use cryo electron microscopy to shed light on the interaction between highly curved helices. We find that the spacing between DNA helices in spermine-induced DNA toroidal condensates depends on their location within the torus, consistent with a mathematical model based on the competition between electrostatic interactions and the bending rigidity of DNA. We use our model to infer the characteristics of the interaction potential, and find that its equilibrium spacing strongly depends on the curvature of the filaments. In addition, the interaction is much softer than previously reported in bulk samples using different salt conditions. Beyond viruses and cells, our characterization of the interactions governing DNA-based dense structures could help develop robust designs in DNA nanotechnologies.

Charge regulation of colloidal particles in aqueous solutions

We study the charge regulation of colloidal particles inside aqueous electrolyte solutions. To stabilize a colloidal suspension against precipitation, colloidal particles are synthesized with either acidic or basic groups on their surface. On contact with water, these surface groups undergo proton transfer reactions, resulting in colloidal surface charge. The charge is determined by the condition of local chemical equilibrium between hydronium ions inside the solution and at the colloidal surface. We use a model of Baxter sticky spheres to explicitly calculate the equilibrium dissociation constants and to construct a theory which is able to quantitatively predict the effective charge of colloidal particles with either acidic or basic surface groups. The predictions of the theory for the model are found to be in excellent agreement with the results of Monte Carlo simulations. This theory is further extended to treat colloidal particles with a mixture of both acidic and basic surface groups.