Molecular dynamics modeling and simulations are employed to study the immobilization of ligands on the surface of the pores of a base porous polymeric matrix. The results show the significant effects that the counter-ions have on the spatial distribution of the density of immobilized ligands as well as on the pore size and pore connectivity distributions of the porous adsorbent medium being constructed. The results for the systems studied in this work indicate that by using doubly charged counter-ions whose numbers during ligand immobilization are half to those of singly charged counter-ions, the ligand immobilization process proceeds faster and the magnitude of local nonelectroneutrality becomes smaller. More importantly, the pore structures of the adsorbent media resulting from the system using doubly charged counter-ions have porous structures that are characterized by more mid-sized pores and higher pore connectivity than the porous adsorbent structures generated by the system employing singly charged counter-ions and, furthermore, the density distribution of the immobilized ligands in the porous structures where doubly charged counter-ions are employed tends to be more uniform laterally and the ligands are surrounded by fewer counter-ions. These characteristics affected by the use of doubly charged counter-ions could provide important advantages with respect to the transport and adsorption of adsorbate biomolecules of interest. Furthermore, the results of this work indicate that the type of counter-ions being used in the ligand immobilization process could represent an additional control variable for affecting the ligand density distribution as well as the pore size and pore connectivity distributions of the porous structure of the adsorbent medium being constructed.