CaltechAUTHORS
  A Caltech Library Service

Gas Sorption and Barrier Properties of Polymeric Membranes from Molecular Dynamics and Monte Carlo Simulations

Cozmuta, Ioana and Blanco, Mario and Goddard, William A., III (2007) Gas Sorption and Barrier Properties of Polymeric Membranes from Molecular Dynamics and Monte Carlo Simulations. Journal of Physical Chemistry B, 111 (12). pp. 3151-3166. ISSN 1520-6106. doi:10.1021/jp062942h. https://resolver.caltech.edu/CaltechAUTHORS:20170609-071851389

[img] PDF - Supplemental Material
See Usage Policy.

680kB

Use this Persistent URL to link to this item: https://resolver.caltech.edu/CaltechAUTHORS:20170609-071851389

Abstract

It is important for many industrial processes to design new materials with improved selective permeability properties. Besides diffusion, the molecule's solubility contributes largely to the overall permeation process. This study presents a method to calculate solubility coefficients of gases such as O_2, H_2O (vapor), N_2, and CO_2 in polymeric matrices from simulation methods (Molecular Dynamics and Monte Carlo) using first principle predictions. The generation and equilibration (annealing) of five polymer models (polypropylene, polyvinyl alcohol, polyvinyl dichloride, polyvinyl chloride-trifluoroethylene, and polyethylene terephtalate) are extensively described. For each polymer, the average density and Hansen solubilities over a set of ten samples compare well with experimental data. For polyethylene terephtalate, the average properties between a small (n = 10) and a large (n = 100) set are compared. Boltzmann averages and probability density distributions of binding and strain energies indicate that the smaller set is biased in sampling configurations with higher energies. However, the sample with the lowest cohesive energy density from the smaller set is representative of the average of the larger set. Density-wise, low molecular weight polymers tend to have on average lower densities. Infinite molecular weight samples do however provide a very good representation of the experimental density. Solubility constants calculated with two ensembles (grand canonical and Henry's constant) are equivalent within 20%. For each polymer sample, the solubility constant is then calculated using the faster (10×) Henry's constant ensemble (HCE) from 150 ps of NPT dynamics of the polymer matrix. The influence of various factors (bad contact fraction, number of iterations) on the accuracy of Henry's constant is discussed. To validate the calculations against experimental results, the solubilities of nitrogen and carbon dioxide in polypropylene are examined over a range of temperatures between 250 and 650 K. The magnitudes of the calculated solubilities agree well with experimental results, and the trends with temperature are predicted correctly. The HCE method is used to predict the solubility constants at 298 K of water vapor and oxygen. The water vapor solubilities follow more closely the experimental trend of permeabilities, both ranging over 4 orders of magnitude. For oxygen, the calculated values do not follow entirely the experimental trend of permeabilities, most probably because at this temperature some of the polymers are in the glassy regime and thus are diffusion dominated. Our study also concludes large confidence limits are associated with the calculated Henry's constants. By investigating several factors (terminal ends of the polymer chains, void distribution, etc.), we conclude that the large confidence limits are intimately related to the polymer's conformational changes caused by thermal fluctuations and have to be regardedat least at microscaleas a characteristic of each polymer and the nature of its interaction with the solute. Reducing the mobility of the polymer matrix as well as controlling the distribution of the free (occupiable) volume would act as mechanisms toward lowering both the gas solubility and the diffusion coefficients.


Item Type:Article
Related URLs:
URLURL TypeDescription
http://dx.doi.org/10.1021/jp062942hDOIArticle
http://pubs.acs.org/doi/abs/10.1021/jp062942hPublisherArticle
http://pubs.acs.org/doi/suppl/10.1021/jp062942hPublisherSupporting Information
ORCID:
AuthorORCID
Goddard, William A., III0000-0003-0097-5716
Additional Information:© 2007 American Chemical Society. Received 14 May 2006. Published online 8 March 2007. Published in print 1 March 2007. We would like to thank Alejandro Strachan for his contribution and valuable suggestions in the initial stages of this project. We thank R. Niemer and Y. Zhou with help in building some of the polymer samples used in this study. This work was supported in part by a grant from the Avery Dennison Corporation. The facilities of the MSC are partly funded by NSF, DOE-ASCI, Chevron, NIH, ONR, Seiko-Epson, Kellogg's, General Motors, Beckman Institute, Asahi Chemical, and Nippon Steel.
Funders:
Funding AgencyGrant Number
Avery-DennisonUNSPECIFIED
NSFUNSPECIFIED
Department of Energy (DOE)UNSPECIFIED
ChevronUNSPECIFIED
NIHUNSPECIFIED
Office of Naval Research (ONR)UNSPECIFIED
Seiko-EpsonUNSPECIFIED
Kellogg'sUNSPECIFIED
General MotorsUNSPECIFIED
Caltech Beckman InstituteUNSPECIFIED
Asahi ChemicalUNSPECIFIED
Nippon SteelUNSPECIFIED
Issue or Number:12
DOI:10.1021/jp062942h
Record Number:CaltechAUTHORS:20170609-071851389
Persistent URL:https://resolver.caltech.edu/CaltechAUTHORS:20170609-071851389
Official Citation:Gas Sorption and Barrier Properties of Polymeric Membranes from Molecular Dynamics and Monte Carlo Simulations Ioana Cozmuta, Mario Blanco, and William A. Goddard III The Journal of Physical Chemistry B 2007 111 (12), 3151-3166 DOI: 10.1021/jp062942h
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:78040
Collection:CaltechAUTHORS
Deposited By: Ruth Sustaita
Deposited On:09 Jun 2017 16:10
Last Modified:15 Nov 2021 17:36

Repository Staff Only: item control page