In accordance with the above introduction, we get a clear understanding of end-grafted PBs, i.e., the interfacial properties including the microstructures and excess free energies of grafted PBs are significantly influenced by the physical constraints of the grafted monomer. In the theoretical numerical implementations of grafted PBs, the density profiles of monomers and counterions are set to be step-function density distributions in the first iterative step (i.e., t = 1) according to the physical constraints of grafted PBs. Then, the new density profiles at the next iterative step are calculated with the help of two key self-consistent determined parameters (i.e., P and C ). The numerical details are introduced as follows.