Using molecular dynamic simulation, we developed an algorithm to protonate the surface of amorphous porous silica grain particle model and study its effect. In this work, can be used cosmic dust coagulation. The cluster was protonated by placing H atoms on oxygen having only a single bond, namely, non-bridging oxygens. are placed opposite Si–O bond with distance around 1 Å form silanol (Si–O–H)...