Abstract:The diffusion behavior of low molecular weight solvents within 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 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.
What problem does this paper attempt to address?
This paper aims to explore the diffusion behavior of low - molecular - weight solvents in unentangled polystyrene films, especially during the transition process between the glassy and rubbery states. Specifically, the research focuses on the Case II or Class II diffusion phenomenon, which is different from the conventional Fickian diffusion, manifested as the solvent concentration front moving at a constant speed, forming a swollen rubber region and an unaffected glass region.
The main questions in the paper are:
- **How to observe and understand the Case II diffusion phenomenon through all - atom molecular dynamics simulations?**
- **What are the kinetic mechanisms of solvent diffusion in glassy polymers?**
To answer these questions, the authors carried out detailed molecular dynamics simulations to study the diffusion behavior of ethyl acetate, toluene and their mixtures in 10 - mer polystyrene films at different temperatures (especially below and above the glass transition temperature). Through these simulations, the authors hope to reveal the details of the interactions between the solvent and the polymer and how these interactions affect the diffusion process.
### Key Point Summary
1. **Differences in Diffusion Behavior**: In the glassy state (e.g., at 298 K), solvent diffusion exhibits Case II diffusion characteristics, that is, the solvent concentration front moves at a constant speed; while in the rubbery state (e.g., at 450 K), the diffusion behavior conforms to Fickian diffusion, with no obvious concentration front.
2. **Simulation Method**: All - atom molecular dynamics simulations were used to study the diffusion behavior of solvents in polystyrene films under different conditions, including different temperatures, different solvent types and their mixtures.
3. **Result Analysis**: Through the analysis of density profiles and the time - dependence of the front position, the existence of Case II diffusion was verified, and the factors affecting the diffusion behavior, such as the re - absorption of polymer chains and the size of solvent molecules, were explored.
### Formulas
- **Fickian Diffusion**: \[ C(x,t) = \frac{C_0}{2} \left[ \operatorname{erfc} \left( \frac{x - vt}{2\sqrt{Dt}} \right) + \exp \left( \frac{vx}{D} \right) \operatorname{erfc} \left( \frac{x + vt}{2\sqrt{Dt}} \right) \right] \]
- **Case II Diffusion**: \[ Z = v \cdot t^n \]
- \( Z \): Front position (unit: Å)
- \( t \): Time (unit: ns)
- \( v \): Linear propagation speed (unit: Å/ns)
- \( n \): Fitting parameter, for typical Case II diffusion, \( n \) is close to 1
Through these studies, the authors hope to provide new insights for understanding and optimizing solvent diffusion in glassy polymers, which has important application values in fields such as drug release, packaging materials and gas separation.