Electrolyte distribution around two like-charged rods: their effective attractive interaction and angular dependent charge reversal

A simple model for two like-charged parallel rods immersed in an electrolyte solution is considered. We derived the three point extension (TPE) of the hypernetted chain/mean spherical approximation (TPE-HNC/MSA) and Poisson-Boltzmann (TPE-PB) integral equations. We numerically solve these equations and compare them to our results of Monte Carlo (MC) simulations. The effective interaction force, F(T), the charge distribution profiles, rho(el)(x,y), and the angular dependent integrated charge function, P(theta), are calculated for this system. The analysis of F(T) is carried out in terms of the electrostatic and entropic (depletion) contributions, F(E) and F(C). We studied several cases of monovalent and divalent electrolytes, for which the ionic size and concentration are varied. We find good qualitative agreement between TPE-HNC/MSA and MC in all the cases studied. The rod-rod force is found to be attractive when immersed in large size, monovalent or divalent electrolytes. In general, the TPE-PB has poor agreement with the MC. For large monovalent and divalent electrolytes, we find angular dependent charge reversal charge inversion and polarizability. We discuss the intimate relationship between this angular dependent charge reversal and rod-rod attraction.