To address the limitations of SSA in simulating large-scale systems, the 𝜏-leaping algorithm has been developed as one of the most notable methods. 𝜏-leaping allows multiple reactions to occur within a single time step by approximating reaction counts as Poisson random variables, significantly reducing computational effort. However, due to the extremely low concentration of radicals in polymerization systems, 𝜏-leaping is not well-suited for accurately modeling polymerization processes.