A Reaction-Diffusion Model for Capturing Mass Loss and Microstructure Evolution in Enzymatic Degradation of...

The microstructure of semicrystalline bioresorbable polymers is central to their biomedical performance because the crystalline content influences both the mechanical stability and the degradation behaviour. Experimental studies have shown that crystallinity evolves concurrently with mass loss during enzymatic degradation. However, most existing models represent the material as a single homogeneous structure, preventing them from capturing this microstructural evolution or the state-selective mechanisms that drive it. We present a one-dimensional partial differential equation model for the enzymatic degradation of thin films, which treats the crystalline and amorphous states as distinct reactive components. Calibrated to poly(ε-caprolactone) (PCL) degraded by Candida antarctica lipase in vitro, the model accurately reproduces both the observed weight-loss profile and the concurrent decline in crystallinity. Parameter uncertainty analysis indicates that while there are varying degrees of confidence in individual parameter values, the overall model predictive uncertainty is well constrained. Parameter sensitivity analysis shows that the amorphous catalytic rate (the rate at which the enzyme degrades the amorphous region) is the dominant driver of degradation dynamics. The identified model parameters are used to explore the role of film thickness on the rates of mass and crystallinity loss. It was found that thin films remain largely reaction-limited, whereas thicker specimens become increasingly transport-influenced, with slower degradation and delayed structural evolution in the material interior. The model provides a useful tool to explore the effect of changing PCL film thickness on degradation rate and crystallinity-related properties without extensive experimentation.