Polyelectrolyte brush bilayers in weak interpenetration regime: Scaling theory and molecular dynamics simulations.

We employ molecular dynamics (MD) simulations and develop scaling theories to quantify the equilibrium behavior of polyelectrolyte (PE) brush bilayers (BBLs) in the weakly interpenetrated regime, which is characterized by d_{0}<d_{g}<2d_{0}, where d_{g} is the gap between the opposing plates where the PE brushes are grafted and d_{0} is the unperturbed height of a PE brush grafted at a single plate. Scaling predictions establish that, for the weakly interpenetrated osmotic PE BBLs δ∼N^{1/2}(2-d_{g}/d_{0})^{1/2} (where δ is the interpenetration length and N is the number of Kuhn segments in PE brush). MD simulations excellently recover this dependence of δ on N and the extent of interpenetration (quantified by d_{g}/d_{0}). These predictions, unlike the existing studies, establish a finite interpenetration for all values of d_{g}/d_{0} as long as d_{g}<2d_{0}. Finally, we quantify the monomer and counterion concentration distributions and point out that these two distributions may quantitatively deviate from each other at locations very close to the channel centerline, where the interpenetration-induced monomer concentration can be significantly low.