Geopolymers offer a potential alternative to ordinary Portland cement owing to their performance in mechanical and thermal properties, as well as environmental benefits stemming from a reduced carbon footprint. This paper endeavors to build upon prior atomistic computational work delving deeper into the intricate relationship between pH levels and the resulting material's properties, including pore size distribution, geopolymer nucleate cluster dimensions, total system energy, and monomer poly-condensation behavior. Coarse-grained Monte Carlo (CGMC) simulation inputs include tetrahedral geometry and binding energy parameters derived from DFT simulations for aluminate and silicate monomers. Elevated pH values may can alter reactivity and phase stability, or, in the structural concrete application, may passivate the embedded steel reinforcement. Thus, we examine the effects of pH values set at 11, 12, and 13 (based on silicate speciation chemistry), investigating their respective contributions to the nucleation of geopolymers. To simulate a larger system to obtain representative results, we propose the numerical implementation of an Octree cell. Finally, we further digitize the resulting expanded structure to ascertain pore size distribution, facilitating a comparative analysis. The novelty of this study is underscored by its expansion in both system size, more accurate monomer representation, and pH range when compared to previous CGMC simulation approaches. The results unveil a discernible correlation between the number of clusters and pores under specific pH levels. This links geopolymerization mechanisms under varying pH conditions to the resulting chemical properties and final structural state.