Coarse-Grained Simulations of Polyrotaxane Hydrogels under Quiescent and Shear Conditions

Cyclodextrin-based polyrotaxanes (PR) form hydrogels in water through cyclodextrin (CD) aggregation and crystallization. These networks can break under shear flows, making them versatile platforms for extrusion-based 3D printing. To optimize the material properties of 3D-printed PR gels, a microscopic understanding of the structural evolution during 3D printing is necessary. Here, we employ coarse-grained (CG) simulations to reveal the PR assembly process at the molecular level. Our simulations reproduce the experimental crystal morphologies of PR at varying concentrations and chain lengths under quiescent conditions. Using nonequilibrium simulations, we show shear flow ruptures crystalline domain connectivity in PR gels and stacks the lamellae in the gradient direction, allowing the materials to flow during 3D printing. After the cessation of flow, the anisotropic crystal alignment and absence of available dangling PR diminish the intercrystal connectivity. The printed materials are therefore mechanically weaker than the pristine hydrogels, in agreement with experimental results. Nonetheless, by relating the microscopic structural evolution with viscoelastic properties of PR gels and solutions, we elucidate how flow conditions and sample composition affect the 3D printing performance of PR hydrogels.