where the first term is for taking into account the entropy of the counterions residing outside brush, the second term for restricting the number of counterions and the last term for inclusion of the electrostatic interactions throughout the region outside brush. In this region, Eq. ( ) is coupled to the following Poisson-Boltzmann differential equation,