Understanding the Response of Poly(ethylene glycol) diacrylate (PEGDA) Hydrogel Networks: A Statistical Mechanics-Based Framework

Thanks to many promising properties, including biocompatibility and the ability to experience large deformations, poly(ethylene glycol) diacrylate (PEGDA) hydrogels are excellent candidate materials for a wide range of applications. Interestingly, the polymerization of PEGDA leads to a network microstructure that is fundamentally different from that of the "classic" polymeric gels. Specifically, PEGDA hydrogels comprise PEG chains that are interconnected by multifunctional densely grafted rod-like polyacrylates (PAs), which serve as cross-linkers. In this work, we derive a microstructurally motivated model that captures the essential features which enable deformation in PEGDA hydrogels: (1) entropic elasticity of PEG chains, (2) deformation of PA rods, and (3) PA-PA interactions. Expressions for the energy-density functions and the stress associated with each of the three contributions are derived. The model demonstrates the microstructural evolution of the network during loading and reveals the role of key microscopic quantities. To validate the model, we fabricate and compress PEGDA hydrogel discs. The model is in excellent agreement with our experimental findings for a broad range of PEGDA compositions. Interestingly, we show that the response of PEGDA hydrogels with short PEG chains and long PA rods is governed by PA-PA interactions, whereas networks with longer PEG chains are dominated by entropy. To enable design, we employ the model to investigate the influence of key microstructural quantities, such as the length of the PEG and the PA chains, on the macroscopic properties and response. The findings from this work pave the way to the efficient design of PEGDA hydrogels with tunable properties and behaviors, which will enable the optimization of their performance in various applications.