All-atom molecular dynamics simulation of solvent diffusion in an unentangled polystyrene film

The diffusion behavior of low molecular weight solvents within an unentangled polystyrene film below and above its glass transition temperature is investigated. The diffusion behavior in the glassy state exhibits a distinct behavior known as case II or class II diffusion, noticeably diverging from conventional Fickian diffusion observed above the glass transition temperature of the polymer film. In the context of case II diffusion, the primary experimental observation entails the emergence of a well-defined concentration front moving at a constant speed, delineating a swollen, rubbery region from a glassy region within the polymer system. Despite the prevalence of this phenomenon in experimental settings, simulating case II diffusion has posed a significant challenge, primarily due to the computationally intensive nature of the diffusion process. To address this, we have developed an all-atom molecular dynamics simulation approach for the observation of case II diffusion in glassy polymers. This method aims to unravel the intricacies of the diffusion process, providing valuable insights into the dynamic interactions between solvents and the polymer matrix.