as the resin prepolymer and DDS (4,4'-diaminodiphenyl sulfone) as the amine hardener (Figure ). While presented on a specific example, it can be, however, easily extended to account for other amine hardeners or adapted to other epoxy prepolymers. We first develop a block chemistry approach by parametrizing all fragments that can appear in the system. This force field allows us to accurately model the dynamic evolution of electrostatic interactions without on-the-fly charge re-computation, all along respecting the overall charge neutrality of the system. We then expand on the existing bond/react procedures for epoxies and establish a three step protocol that can achieve high curing without ad-hoc iterative protocols, pre-activated reactive groups or excessive distance cut-offs. This allows us to simulate large systems with high accuracies, as determined by applying rigorous statistical procedures for extracting macroscopic properties. The obtained results are found to excellently match experimentally available data, which shows that the procedures established herein could be utilized in the broader context of quantitative modeling of thermoset polymers.