Molecular dynamics simulation of O₂ diffusion in polydimethylsiloxane (PDMS) and end-linked PDMS networks

Abstract

Molecular dynamics simulations are used to compute diffusion coefficients for O₂ molecules in polydimethylsiloxane (PDMS) and end-linked PDMS networks. The PDMS chains and penetrants are modelled using a hybrid interatomic potential which treats the Si and O atoms along the chain backbone explicitly while coarse-graining the methyl side groups and penetrants. In PDMS models with different molecular weights, diffusivity of the O₂ penetrants is found to modestly decrease with an increase in chain length. To match typical experimental conditions, the end-linked PDMS networks are constructed with a PDMS to crosslinking (CL) molecule mass ratio of 5:1 or 10:1, demanding that the number of CL molecules exceeds the number of PDMS chains in each model. Despite end-linking, the presence of non-bonded CL molecules promotes increased O₂ diffusivity in comparison with uncrosslinked PDMS. Temperature dependence is captured using the Williams–Landel–Ferry equation.

Keywords::

• molecular dynamics
• diffusion
• polydimethylsiloxane
• crosslinking

Acknowledgements

This work was supported by the National Science Foundation under grant CMMI#0800718. MD simulations were carried out on cyberinfrastructure resources at the University of Arkansas which are supported in part by the National Science Foundation under Grants ARI#0963249, MRI#0959124 (Razor), EPS #0918970 (CI TRAIN) and a grant from the Arkansas Science and Technology Authority.