 Solution

For large positive dimensionless surface potential $u_0=e{\Phi }_0/(k_{B}T)\gtrsim 1$, use the nonlinear Boltzmann distribution rather than their linear expansion:The counterions become strongly enriched near the wall and co-ions strongly depleted. A planar mean-field solution isIt approaches the linear screening form far from the wall. Nonlinearity preserves positive densities and produces a much stronger near-wall asymmetry; an extrapolated linear co-ion profile could incorrectly become negative. The total local ion density is $2n_0\cosh u$, no longer constant.

The original figure compares weak- and strong-potential profiles. At very large local concentrations, finite ion size, specific adsorption and ion correlations limit this ideal point-ion mean-field model; it does not predict unlimited physically realizable counterion packing.