We finally determined the Young's modulus of our system by performing a tensile stress simulation. To make sure that the system was not trapped in a metastable equilibrium, which is known to affect the mechanical properties of MD polymeric networks, we subjected each cured sample to a post-curing annealing procedure consisting of 10 ns at a temperature of 1000 K, after which the temperature was quickly lowered to 300 K. Each system was then allowed to equilibrate in the NPT ensemble at 298 K using isotropic pressure coupling until its density reached steady state, with equilibration times reported in Table . To rule out the possibility of insufficient annealing time, we also extended the post-curing annealing period up to 20 𝑛𝑠, but without observing significant changes in the mechanical response of the system. After completing equilibration at room temperature, each cross-linked sample underwent a uniaxial tensile test along the x, y, and z axes, at an engineering strain rate of 2 10 -8 s -1 .