Phantom Chain Simulations for the Fracture of Energy-Minimized Tetra- and Tri-Branched Networks

Abstract

Effects of branch functionality on mechanical properties of polymer networks are yet to be fully elucidated, although multi-functional approaches have been mainly attempted. For instance, Fujiyabu et al. [Sci. Adv. 2022, 8, abk0010] recently reported that polymer networks made from tri-branch prepolymers exhibit superior mechanical properties to tetra-branch analogues. Although they attribute the difference to stretch-induced crystallization observed in tri-branch poly (ethylene glycol) networks, the mechanism still needs to be clarified. In this study, we performed coarse-grained molecular simulations to extract the effect of branch functionality. The prepolymers were replaced by bead-spring phantom chains, and gelation was simulated by a Brownian dynamics scheme. We subjected the resultant networks to energy minimization and uniaxial stretch by introducing breakage for elongated segments. In the stress–strain relation thus obtained, stress and strain at the break were larger for tri-branch networks than that for tetra-branch analogues, consistent with the experiment. The superiority of tri-branch networks is observed in a wide range of the conversion ratio in gelation, molecular weights of prepolymers, and polymer concentrations. The result implies that the mechanical superiority of tri-branch networks to tetra-branch ones is due to a fundamental structural difference generated during gelation.