Malek O. Khan

The influence of discrete surface charges on the force between charged surfaces.

The force between two parallel charged flat surfaces, with discrete surface charges, has been calculated with Monte Carlo simulations for different values of the electrostatic coupling. For low electrostatic coupling (small counterion valence, small surface charge, high dielectric constant, and high temperature) the total force is dominated by the entropic contribution and can be described by mean field theory, independent of the character of the surface charges. For moderate electrostatic coupling, counterion correlation effects lead to a smaller repulsion than predicted by mean field theory. This correlation effect is strengthened by discrete surface charges and the repulsive force is further reduced. For large electrostatic coupling the total force for smeared out surface charges is known to be attractive due to counterion correlations. If discrete surface charges are considered the attractive force is weakened and can even be turned into a repulsive force. This is due to the counterions being strongly correlated to the discrete surface charges forming effective, oppositely directed, dipoles on the two walls.

Interaction between charged surfaces mediated by rodlike counterions: the influence of discrete charge distribution in the solution and on the surfaces.

The interaction between two charged surfaces, with discrete or uniform charge distributions, embedded in a solution of rodlike counterions has been studied. Monte Carlo simulations and density functional theory have been applied to study the concentration profiles of counterions and the force between the charged surfaces. We show that for low surface charge densities the repulsive force between like-charged surfaces is observed regardless of the rod length. Where high surface charge densities are present, attractive forces at surface separations related to the rod length are observed.

A scalable parallel Monte Carlo method for free energy simulations of molecular systems

We present a method of parallelizing flat histogram Monte Carlo simulations, which give the free energy of a molecular system as an output. In the serial version, a constant probability distribution, as a function of any system parameter, is calculated by updating an external potential that is added to the system Hamiltonian. This external potential is related to the free energy. In the parallel implementation, the simulation is distributed on to different processors. With regular intervals the modifying potential is summed over all processors and distributed back to every processor, thus spreading the information of which parts of parameter space have been explored. This implementation is shown to decrease the execution time linearly with added number of processors.

Decreased Osmotic Pressure via Interfacial Charge Clustering

Behaviors of a model interfacial system featuring the clustering of discrete, mobile wall charges and a counterion solution are investigated. The results demonstrate that even very small localized charge clusters produce significant effects for the osmotic pressure, effects that are not adequately represented in common colloidal models. We observe a pronounced decrease in osmotic pressure where a certain level of clustering is attained, with potentially significant implications for theories of colloidal stability and biochemical processes at microscopic length scales. The stable level of wall charge clustering, and thus the osmotic pressure, is found to be variable in the addition of an attractive potential, as measured via the excess free energy of clustering.