Surface-grafted polyelectrolytes provide a versatile way to create functionalized interfaces and nanochannels with externally controllable properties. Understanding the behavior of ions within the brush-like assemblies is crucial for the further development of these devices. We demonstrate that the ion transport through the brushes is governed by the interplay of electrostatic ion–polymer binding and steric effects, leading to a mobility that depends nonmonotonically on grafting density. However, the ion–polymer binding can be modulated by the dielectric properties of the substrate. As a result, surface polarization suppresses ion mobility near insulating interfaces and enhances it near conducting interfaces, even causing a shift from nonmonotonic to monotonic variation with grafting density.

The net charge of solvated entities, ranging from polyelectrolytes and biomolecules to charged nanoparticles and membranes, depends on the local dissociation equilibrium of individual ionizable groups. Incorporation of this phenomenon, charge regulation (CR), in theoretical and computational models requires dynamic, configuration-dependent recalculation of surface charges and is therefore typically approximated by assuming constant net charge on particles. Various computational methods exist that address this. We present an alternative, particularly efficient CR Monte Carlo method (CR-MC), which explicitly models the redistribution of individual charges and accurately samples the correct grand-canonical charge distribution. In addition, we provide an open-source implementation in the large-scale Atomic/Molecular Massively Parallel Simulator molecular dynamics (MD) simulation package, resulting in a hybrid MD/CR-MC simulation method. This implementation is designed to handle a wide range of implicit-solvent systems that model discreet ionizable groups or surface sites. The computational cost of the method scales linearly with the number of ionizable groups, thereby allowing accurate simulations of systems containing thousands of individual ionizable sites. By matter of illustration, we use the CR-MC method to quantify the effects of CR on the nature of the polyelectrolyte coil–globule transition and on the effective interaction between oppositely charged nanoparticles.

Understanding the structural response of polyelectrolyte brushes to variation in both intrinsic and external properties is highly relevant for their application as functionalized interfaces and components of nanodevices. Using coarse-grained simulations, we examine an aspect that is largely unexplored, namely the dielectric mismatch between the solvent and the substrate. We systematically study how this permittivity contrast alters the brush structure over a range of Bjerrum length, polymer charge density, counterion size and valency, salt concentration, polymer grafting density, and external electric field. In addition to the expected brush contraction near metallic substrates and expansion on low-permittivity substrates, we find various regimes where variation of the substrate properties qualitatively alters the brush response.

We have developed an accurate and efficient method for molecular dynamics simulations of charged particles confined by planar dielectric interfaces. The algorithm combines the image-charge method for near field with the harmonic surface mapping, which converts the contribution of infinite far-field charges into a finite number of charges on an auxiliary spherical surface. We approximate the electrostatic potential of far-field charges via spherical harmonic expansion and determine the coefficients by fitting the Dirichlet-to-Neumann boundary condition, which only requires the potential within the simulation cell. Instead of performing the direct evaluation of spherical harmonic series expansion, we use Green’s second identity to transform the series expansion into a spherical integral, which can be accurately represented by discrete charges on the sphere. Therefore, the fast multipole method can be readily employed to sum over all charges within and on the sphere, achieving truly linear O(N) complexity. Our algorithm can be applied to a broad range of charged complex fluids under dielectric confinement.