772a

Theodora Spyriouni^{1}, Georgios C. Boulougouris^{2}, and **Doros N. Theodorou**^{2}. (1) Institute of Physical Chemistry, NCSR "Demokritos", 15310 Aghia Paraskevi Attikis, Athens, Greece, (2) School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Zografou Campus, Athens, 157 80, Greece

A molecular simulation strategy allowing the computation of sorption isotherms, up to high penetrant activities, in glassy polymer matrices will be presented. The new scheme permits an equilibrium repartition of the penetrant molecules in the polymer matrix. Glassy matrices of polystyrene (PS) were obtained from reverse mapping of equilibrated coarse-grained configurations (Spyriouni et al., *Macromolecules* **2007**, *40*, 3876). The sorption isotherms of CO_{2} in atactic PS and the induced polymer swelling have been calculated in the temperature range from 308 K to 405 K (below and above the glass transition temperature of the pure polymer) and for pressures up to 300 bar. The results compare favorably with available experimental data. Derivative properties, such as the partial molar volume and the partial molar enthalpy of CO_{2}, have been calculated in good agreement with experimental values. The segmental dynamics of PS have been analyzed and the pressure induced glass transition has been calculated for various temperatures below the *T*_{g} of pure PS.

The glassy matrices are loaded with CO_{2} at various concentrations by insertions of CO_{2} molecules, after mapping the accessible free volume. The loaded matrices are subjected to molecular dynamics (MD) simulations at constant temperature, pressure and number of molecules (*N*_{1}*N*_{2}*PT* ensemble). The penetrant is redistributed among the accessible cavities in order to achieve an equilibrium repartitioning in the thermally fluctuating polymer matrix. This is accomplished by selecting periodically one configuration and attempting a certain number of CO_{2} displacements between the accessible cavities. The latter are determined by Delaunay tessellation followed by volume and connectivity analysis, as implemented by Greenfield and Theodorou (*Macromolecules* **1993**, *26*, 5461). Following acceptance or rejection of the attempted repartitioning moves, we proceed with the MD simulation of the polymer-penetrant system.

The excess chemical potential of CO_{2} is calculated using the Direct Particle Deletion (DPD) method, a generalization of Staged Particle Deletion (Boulougouris et al. *Mol. Phys.* **1999**, *96*, 905; *J. Chem. Phys.*, **2001**, *115*, 8231). This method is very powerful for dense systems and/or large solute molecules in comparison to traditional Widom insertion. Moreover, the DPD version of the method permits the calculation of the excess chemical potential from a single simulation run.

The fugacity of CO_{2} in the polymer matrix is calculated from the excess chemical potential via the relationship: *f*_{CO2} = *N*_{2}*k*_{B}*T* exp(bm^{ex})/<V>, where *N*_{2} is the number of CO_{2} molecules in the polymer matrix and <V> is the average volume under *NPT* conditions. A new pressure is calculated via an equation of state for pure gaseous CO_{2}, such as SAFT, by setting the gas-phase fugacity equal to *f*_{CO2}. The calculations are repeated at the new pressure, until the CO_{2} fugacity in the polymer matrix and in the gas are equal. The scheme converges rapidly, since the CO_{2} fugacity in the polymer is not very sensitive to the pressure.

Web Page: comse.chemeng.ntua.gr

See more of #772 - Thermodynamics of Polymers II (08A04)

See more of Materials Engineering and Sciences Division

See more of The 2008 Annual Meeting

See more of Materials Engineering and Sciences Division

See more of The 2008 Annual Meeting