The porous morphology drives the application of nano-and mesoporous materials in a variety of areas. However, limited information is available on the porosity development during their synthesis at various times due to challenges in experimental and simulation techniques. In this work, we have probed the porosity development in silica particles using Monte Carlo simulation techniques. We have developed an algorithm to measure the porosity of small, irregular shaped, finite-size particles formed during the polymerization process. We observed that smaller and denser clusters are formed initially, which later aggregate to form porous clusters in the range of 500−1000 h for the system having a concentration of silica precursor = 104 mg/cm 3 in a reactive solvent. In the case of nonreactive solvents, the porosity development is significantly different and occurs from initial stages of polymerization and lasts until the aggregation stage, in the range of 2−2000 h for similar concentrations. This significant change is due to faster kinetics of polymerization. Further polymerization leads to the formation of denser clusters due to aging in both cases. The mechanism of porosity formation in reactive solvent systems is due to random aggregation of denser clusters, whereas in a nonreactive solvent, it occurs by merging of smaller porous clusters. We also prepared a phase diagram of porosity evolution at various concentrations and observed that it has a tremendous effect on porosity development. We find three concentration ranges where silica cluster transforms from small denser to large denser particles in multiple stages. For a nonreactive solvent, the porosity evolution phase diagram shifts and stretches toward lower concentration ranges. We believe that this detailed understanding of porosity development will be useful to control the porous morphology of nano-and mesoporous materials during their synthesis.