Development of a radical polymerization algorithm for molecular dynamics simulations of antifreezing hydrogels with double-network structures

The development and understanding of antifreezing hydrogels are crucial both in principle and practice for the design and delivery of new materials. The current antifreezing mechanisms in hydrogels are almost exclusively derived from their incorporation of antifreezing additives, rather than from the inherent properties of the polymers themselves. Moreover, developing a computational model for the independent yet interconnected double-network (DN) structures in hydrogels has proven to be an exceptionally difficult task. Here, we develop a multiscale simulation platform, integrating ‘random walk reactive polymerization’ (RWRP) with molecular dynamics (MD) simulations, to computationally construct a physically-chemically linked PVA/PHEAA DN hydrogels from monomers that mimic a radical polymerization and to investigate water structures, dynamics, and interactions confined in PVA/PHEAA hydrogels with various water contents and temperatures, aiming to uncover antifreezing mechanism at atomic levels. Collective simulation results indicate that the antifreezing property of PVA/PHEAA hydrogels arises from a combination of intrinsic, strong water-binding networks and crosslinkers and tightly crosslinked and interpenetrating double-network structures, both of which enhance polymer-water interactions for competitively inhibiting ice nucleation and growth. These computational findings provide atomic-level insights into the interplay between polymers and water molecules in hydrogels, which may determine their resistance to freezing.