Influence of solvent polarization and non-uniform ion size on electrostatic properties between charged surfaces in an electrolyte solution

Development of a coarse-grained model for surface-functionalized gold nanoparticles: towards an accurate description of their aggregation behavior

Understanding the dispersion stability and aggregation propensity of self-assembled monolayer gold NPs at a molecular level is crucial to guide their rational design and to inform about the optimal surface functionalization for specific applications. To reach this goal, in silico modeling via coarse-grained (CG) molecular dynamics (MD) simulations is a fundamental tool to complement the information acquired from experimental studies since CG modeling allows to get a deep knowledge of the molecular interactions that take place at the nanoscale in this kind of systems. Unfortunately, current CG models of monolayer-protected AuNPs present several drawbacks that limit their accuracy in certain scenarios. We here develop a CG model that is fully compatible and extends the SPICA/SDK (Shinoda-DeVane-Klein) force field. Our model allows reproducing the behavior of AuNPs functionalized with hydrophobic as well as charged and more hydrophilic ligands. This model improves upon results obtained with previously derived CG force fields and successfully describes NPs aggregation and self-assembly in aqueous solution.

Monovalent ion-mediated charge-charge interactions drive aggregation of surface-functionalized gold nanoparticles

Monolayer-protected metal nanoparticles (NPs) are not only promising materials with a wide range of potential industrial and biological applications, but they are also a powerful tool to investigate the behaviour of matter at nanoscopic scales, including the stability of dispersions and colloidal systems. This stability is dependent on a delicate balance between attractive and repulsive interactions that occur in the solution, and it is described in quantitative terms by the classic Derjaguin-Landau-Vewey-Overbeek (DLVO) theory, that posits that aggregation between NPs is driven by van der Waals interactions and opposed by electrostatic interactions. To investigate the limits of this theory at the nanoscale, where the continuum assumptions required by the DLVO theory break down, here we investigate NP dimerization by computing the Potential of Mean Force (PMF) of this process using fully atomistic MD simulations. Serendipitously, we find that electrostatic interactions can lead to the formation of metastable NP dimers at physiological ion concentrations. These dimers are stabilized by complexes formed by negatively charged ligands belonging to distinct NPs that are bridged by positively charged monovalent ions present in solution. We validate our findings by collecting tomographic EM images of NPs in solution and by quantifying their radial distribution function, that shows a marked peak at interparticle distance comparable with that of MD simulations. Taken together, our results suggest that not only van der Waals interactions, but also electrostatic interactions mediated by monovalent ions at physiological concentrations, contribute to attraction between nano-sized charged objects at very short length scales.

Structural and electrostatic properties between pH-responsive polyelectrolyte brushes studied by augmented strong stretching theory

In this paper, we study electrostatic and structural properties between pH-responsive polyelectrolyte brushes by using a strong stretching theory accounting for excluded volume interactions, the density of polyelectrolyte chargeable sites and the Born energy difference between the inside and outside of the brush layer.In a free energy framework, we obtain self-consistent field equations to determine electrostatic properties between two pH-responsive polyelectrolyte brushes. We elucidate that in the region between two pH-responsive polyelectrolyte brushes, electrostatic potential at the centerline and osmotic pressure increase not only with excluded volume interaction, but also with density of chargeable sites on a polyelectrolyte molecule. Importantly, we clarify that when two pH-responsive polyelectrolyte brushes approach each other, the brush thickness becomes short and that a large excluded volume interaction and a large density of chargeable sites yield the enhanced contract of polyelectrolyte brushes. In addition, we also demonstrate how the influence of such quantities as pH, the number of Kuhn monomers, the density of charged sites, the lateral separation between adjacent polyelectrolyte brushes, Kuhn length on the electrostatic and structural properties between the two polyelectrolyte brushes is affected by the exclusion volume interaction. Finally, we investigate the influence of Born energy difference on the thickness of polyelectrolyte brushes and the osmotic pressure between two pH-responsive polyelectrolyte brushes.

Effect of Bjerrum pairs on the electrostatic properties of an electrolyte solution near charged surfaces: a mean-field approach

In this paper, we investigate the consequences of ion association, coupled with the considerations of finite size effects and orientational ordering of Bjerrum pairs as well as ions and water molecules, on the electric double layer near charged surfaces. Based on the lattice statistical mechanics accounting for finite sizes and dipole moments of ions, Bjerrum pairs and solvent molecules, we consider the formation of Bjerrum pairs and derive the mathematical expressions for Bjerrum pair number density as well as cation/anion number density and water molecule number density. We reveal several significant phenomena. Firstly, it is shown that our approach naturally yields the equilibrium constant for dissociation-association equilibrium between Bjerrum pairs and ions. Secondly, at low surface charge densities, an increase in the bulk concentration of Bjerrum pairs enhances the permittivity and decreases the differential capacitance. Next, for the cases where Bjerrum pairs in an alcohol electrolyte solution have a high value of dipole moment, the Bjerrum pair number density increases with decreasing distance from the charged surface, and the differential capacitance and permittivity are high compared to those for the cases with lower values of Bjerrum-pair dipole moments. Finally, we show that the difference in the concentration and dipole moment of Bjerrum pairs can lead to some variation in osmotic pressure between two similarly charged surfaces.

Multi-scale approach for modeling stability, aggregation, and network formation of nanoparticles suspended in aqueous solutions

Suspensions of nanoparticles (NPs) in aqueous solutions hold promise in many research fields, including energy applications, water desalination, and nanomedicine. The ability to tune NP interactions, and thereby to modulate the NP self-assembly process, holds the key to rationally synthesize NP suspensions. However, traditional models obtained by coupling the DLVO (Derjaguin, Landau, Verwey, and Overbeek) theory of NP interactions, or suitable modifications of it, with the kinetic theory of colloidal aggregation are inadequate to precisely model NP self-assembly because they neglect hydration forces and discrete-size effects predominant at the nanoscale. By synergistically blending molecular dynamics and stochastic dynamics simulations with continuum theories, we develop a multi-scale (MS) model, which is able to accurately predict suspension stability, timescales for NP aggregation, and macroscopic properties (e.g., the thermal conductivity) of bare and surfactant-coated NP suspensions, in good agreement with the experimental data. Our results enable the formulation of design rules for engineering NP aqueous suspensions in a wide range of applications.