Accelerated ReaxFF simulations of biodegradation of polymers

Biodegradable polymers with elastic properties have gained tremendous attention in the field of tissue engineering. Polyesters are prevalent biomaterials used for fabricating scaffolds for bone tissues. Recently, the reactivity of hexanediol-co-citric acid was studied using accelerated reactive molecular dynamics simulations.

Due to the high activation energy barrier of ester and ether hydrolysis reaction we use the bond boost scheme within the ReaxFF framework (see also our tutorial). The selectivity of the bond boost parameter was established by tuning the bond boost parameters. Ester hydrolysis occurs with a lower barrier than ether hydrolysis. On reducing the bond boost parameters, the ether hydrolysis almost stops. The polymers were also mechanically tested to find its tensile modulus, showing strain-dependent behavior (see also tutorial on stress-strain). Polyesters are found to be more ductile than polyesther-ethers.

ReaxFF transition state structure of reaction between hexanediol-co-citric acid and two water molecules (top), degraded polymer bundle containing multiple chains on hydrolysis, and then applying mechanical strain (middle), single degraded polymer chain (bottom).

Let us know if you would like an example input!

Do you want try for yourself? Request a free 30-day evaluation for the whole Amsterdam Modeling Suite.

Newsletter: tips & tricks, highlights, events

Would you like to keep up to date with the latest developments in the Amsterdam Modeling Suite and the SCM team, learn more about new applications and functionality?
Subscribe to our newsletter!

Nabankur Dasgupta, Dundar E. Yilmaz, and Adri van Duin, Simulations of the Biodegradation of Citrate-Based Polymers for Artificial Scaffolds Using Accelerated Reactive Molecular Dynamics, J. Phys. Chem. B 124, 5311–5322 (2020)

Key concepts