Sorption of light gases in glassy poly(ethyl methacrylate)

Although gas sorption in glassy polymers is a well‐studied phenomenon, no general microscopical model is developed which is able to describe the gas sorption in a wide temperature range using only characteristics of polymer and gas molecule. In this work, sorption isotherms and desorption kinetics of O₂, Ar, and N₂ for glassy poly(ethyl methacrylate) have been measured in the temperature range from 160 to 308 K. To describe both the phenomena, the model is developed which postulates that, in the frozen structure of glassy polymer, any cavities between macromolecules are the sorption sites for small molecules. The cavities of small size can expand elastically to accommodate a gas molecule. The sorption sites are considered to be the potential wells and their depths are distributed according to Gaussian law. The concentration of sorption sites, their mean depth and depths dispersion, and the frequency of molecules oscillations in the sorption sites are the only parameters which determine both the gas transport and sorption. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2018 , 56, 288–296     

Sorption and diffusion of low molecular weight gases and vapors in glassy polymers at high concentration of sorbate

A theoretical approach has been developed to describe the sorption and diffusion processes of low weight molecular gases and vapors in polymers at wide ranges of sorbate concentration. The equation of an S‐shaped gas sorption isotherm in glassy polymer matrix has been derived. The concentration dependence of the sorbate molecule diffusion coefficient has been established. For an S‐shaped sorption isotherm, this dependence is nonmonotonous. The conditions of cluster formation of sorbate molecules have been analyzed within the proposed approach, in which it is possible to determine a correlation between these conditions and parameters of sorption isotherm. The comparison of the experimental and theoretical data provides an assessment of the microscopic characteristics of investigated polymer–vapor systems, such as the distances between vapor molecules in a matrix corresponding to intermolecular repulsion and attraction. © 1999 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 37: 2314–2323, 1999     

Evaluation of gas sorption parameters and prediction of sorption isotherms in glassy polymers

A model of continuous‐site distribution for gas sorption in glassy polymers is examined with sorption data of CO₂ and Ar in polycarbonate. A procedure is presented for determining from a measured isotherm the number of sorption sites in a polymer, an important parameter that previously had to be assumed. With this parameter value and solubility data obtained at zero pressure, the model can reasonably predict sorption isotherms of CO₂ in glassy polycarbonate for a wide temperature range. The number of sorption sites and the average site volume evaluated from CO₂ sorption isotherms are employed for the prediction of Ar sorption isotherms with zero‐pressure solubility data and the independently measured partial molar volume of Ar. A reasonable fit to the measured isotherms of Ar is achieved. With the proposed procedure, the continuous‐site model shows several advantages over the conventional dual‐mode sorption model. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 883–888, 2000     

Molecular motions in a polycarbonate composite as studied by thermally stimulated recovery and dynamic mechanical analysis

The viscoelastic properties of a short carbon-fibre reinforced polycarbonate are studied in the glassy state and in the glass transition region by means of Thermally Stimulated Recovery (TSR) and Dynamic Mechanical Analysis (DMA) techniques in the flexural mode. The temperature calibration of the apparatus is discussed. Each thermal sampling experiment is analysed according to a simple Voigt-Kelvin model in which the temperature dependence of the retardation time follows the Arrhenius equation with two adjustable thermokinetic parameters. Non-linear fittings of the experimental data are performed with the corresponding constitutive equation. Both the activation energy and the pre-exponential factor follow the usual trends with temperature and present the well-known compensation phenomena. A divergence between the TSR results of the two materials is found above ± 145°C and may be a result of the effect of the interface. The observed movements in the glassy state which are still found to have a cooperative character are analysed according to the Adam-Gibbs theory. Dynamic-mechanical relaxation properties of the composite are reported and analysed in terms of Takayanagi's block model.     

The Cubic-Two-State Equation of State: Cross-associating mixtures and Monte Carlo study of self-associating prototypes

A new equation of state for associating fluids has recently been presented by Medeiros and Tellez-Arredondo, the Cubic-Two-State Equation of State (CTS EoS) [Ind. Eng. Chem. Res. 47 (2008) 5723]. This equation arises from the coupling of the Soave–Redlich–Kwong EoS (SRK) with an association term from a two-state association model. The CTS EoS is polynomial in volume and it is able to describe vapor pressures and molar volume of associating fluids such as water, alcohol and phenol, among others. The equation is also able to describe the liquid–vapor equilibria of their mixtures with alkanes. In this paper, the physical and thermodynamic foundations of the CTS EoS are further investigated. In order to verify its applicability for cross-associating systems, the equation was employed in the prediction of phase equilibria behavior of binary alcohol–alcohol and water–alcohol mixtures. Very good agreement between predictions and experimental phase equilibria data was obtained with very simple combining rules and only one adjustable binary parameter. No additional parameters were necessary to describe ternary systems. With the purpose of checking the model's hypothesis and limitations, the two-state association term was coupled with the hard sphere Carnahan–Starling EoS, forming the CS-TS equation and the association characteristic parameters were determined theoretically for prototype association fluids. Monte Carlo NPT simulations of such fluids were performed and the results were compared with the equation's predictions. The CS-TS was able to describe qualitatively the pvT$pvT$ behavior of the prototype; nevertheless, it is not as accurate as those predictions obtained from the combination CS with Wertheim's association term. It seems that, when adjusting parameters of the CTS EoS to real substances, the discrepancies between the predicted and the real association contribution are dissipated among other adjustable parameters, specially on the dispersive term of the SRK equation. Finally, it is shown that CTS EoS isotherms can only have one or three real bigger roots than the co-volume for positive pressures, similar to cubic equations of state, and then it has the desirable form to describe vapor–liquid phase equilibria of associating compounds mixtures.     

Solving the adsorption integral equation with Langmuir-kernel and the influence of temperature on the stability of the solution

Adsorption energy distribution functions can be calculated from measured adsorption isotherms by solving the adsorption integral equation. In this context, it is common practice to use general regularization methods, which are independent of the kernel of the adsorption integral equation, but do not permit error estimation. In order to overcome this disadvantage, we present in this paper a solution theory which is tailor-made for the Langmuir kernel of the adsorption integral equation. The presented theory by means of differentiation and Fourier series is the basis for a regularization method with explicit terms for error amplification. By means of simple and complicated adsorption energy distribution functions we show for ideal gas adsorption isotherms without measurement error that reliable distribution functions can be obtained from the isotherms. Furthermore we show how the stability of the solution depends on temperature.     

On the sorption of gas mixtures by polymer glasses

It was realized recently that various phenomena, related to the sorption of small molecules in polymer glasses could be described within the framework of a site distribution (SD) model. According to the SD model, non-equilibrium structure of glassy polymer leads to the distribution of sorption energies within the interchain holes. The parameters of the distribution for the given polymer–gas system could be expressed through the polymer–gas characteristics or evaluated from the experimental pressure–concentration isotherms. In this work we show how these parameters could be used to predict the sorption isotherms for gas mixtures. The suppression of solubility of each component by the other components, which is the main feature of mixed sorption by polymer glasses can be described within the SD model through the competitive occupancy of low-energy sorption sites. The clear physical meaning of the energy distribution parameters allows to analyze the role of different factors on the competitive sorption from gas mixtures. The comparison of SD model with the other theoretical approaches are given and new experiments, which could check the validity of our approach are proposed.     

The effect of lattice compressibility on the thermodynamics of gas sorption in polymers

A compressible lattice model with holes, the glassy polymer lattice sorption model (GPLSM), was used to model the sorption of carbon dioxide, methane, and ethylene in glassy polycarbonate and carbon dioxide in glassy tetramethyl polycarbonate. For glassy polymers, an incompressible lattice model, such as the Flory–Huggins theory, requires concentration-dependent and physically unrealistic values for the lattice site volumes in order to satisfy lattice incompressibility. Rather than forcing lattice incompressibility, GPLSM was used and reasonable parameter values were obtained. The effect of conditioning on gas sorption in glassy polymers was analyzed quantitatively with GPLSM. The Henry's law constant decreases significantly upon gas conditioning, reflecting changes in the polymer matrix at infinite dilution. Treating the Henry's law constant as a hypothetical vapor pressure at infinite dilution, gas molecules in the conditioned polymer are less “volatile” than those in the unconditioned polymer. Flory–Huggins theory was used to model the sorption of carbon dioxide, methane, and ethylene in silicone rubber. Above the glass transition temperature, the criterion of lattice incompressibility for Flory-Huggins theory was satisfied with physically realistic and constant values for the lattice site volumes. © 1992 John Wiley & Sons, Inc.     

Sorption of methanol/MTBE and diffusion of methanol in 6FDA–ODA polyimide

This article discusses the diffusion and solubility behavior of methanol/methyl tert‐butyl ether (MTBE) in glassy 6FDA–ODA polyimide prepared from hexafluoroisopropylidene 2,2‐bis(phthalic anhydride) (6FDA) and oxydianiline (ODA). The diffusion coefficients and sorption isotherm of methanol vapor in 6FDA–ODA polyimide at various pressures and film thicknesses were obtained with a McBain‐type vapor sorption apparatus. Methanol/MTBE mixed‐liquid sorption isotherms were obtained by head‐space chromatography and compared with a pure methanol sorption isotherm obtained with a quartz spring balance. Methanol sorption isotherms obtained with the two methods were almost identical. Both methanol sorption isotherms obeyed the dual‐mode model at a lower activity, which is typical for glassy polymer behavior. The MTBE was readily sorbed into the polymer in the presence of methanol, but the MTBE sorption isotherm exhibited a highly nonideal behavior. The MTBE sorption levels were a strong function of the methanol sorption level. Methanol diffusion in the polymer was analyzed in terms of the partial immobilization model with model parameters obtained from average diffusion coefficients and the dual‐mode sorption parameters. Simple average diffusion coefficients were obtained from sorption kinetics experiments, whereas the dual‐mode sorption parameters were obtained from equilibrium methanol sorption experiments. An analysis of the mobility and solubility data for methanol indicated that methanol tends to form clusters at higher sorption levels. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 2254–2267, 2000     

Gas solubility of carbon dioxide in poly(lactic acid) at high pressures: Thermal treatment effect

The sorption of carbon dioxide in glassy Poly(lactic acid) (PLA) films was studied by quartz crystal microbalance (QCM) at high pressures. Two thermal treatments, melted and quenched, were performed in PLA with two different L:D contents, 80:20 and 98:2, films and compared with a third thermal protocol, annealed, and used in a previous work. The results obtained show that for pressures higher than 2 MPa, the carbon dioxide solubility is larger in PLA 80:20 than in PLA 98:2, indicating that the L:D plays a dominant role on this property. The thermal treatments only affect the gas solubility in PLA 98:2. Sorption isotherms at temperatures 303, 313, and 323 K, below the glass transition temperature of the polymer, and pressures up to 5 MPa were measured and analyzed with three different models, the dual‐mode sorption model, the Flory–Huggins equation, and a modified dual‐mode sorption model where the Henry's law term was substituted by the Flory–Huggins equation. This last model performs especially well for CO₂ in PLA 80:20, due to the convex upward curvature of the solubility isotherms for that system. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 616–625, 2007     

Permeation of CO2 and N2 through glassy poly(dimethyl phenylene) oxide under steady- and presteady-state conditions

Glassy polymers are often used for gas separations because of their high selectivity. Although the dual-mode permeation model correctly fits their sorption and permeation isotherms, its physical interpretation is disputed, and it does not describe permeation far from steady state, a condition expected when separations involve intermittent renewable energy sources. To develop a more comprehensive permeation model, we combine experiment, molecular dynamics, and multiscale reaction–diffusion modeling to characterize the time-dependent permeation of N₂ and CO₂ through a glassy poly(dimethyl phenylene oxide) membrane, a model system. Simulations of experimental time-dependent permeation data for both gases in the presteady-state and steady-state regimes show that both single- and dual-mode reaction–diffusion models reproduce the experimental observations, and that sorbed gas concentrations lag the external pressure rise. The results point to environment-sensitive diffusion coefficients as a vital characteristic of transport in glassy polymers.     

On the crystallinity effect on the gas sorption in semicrystalline linear low density polyethylene (LLDPE)

This article describes the solubility of carbon dioxide, ethylene and propane in 1‐octene based polyethylene of 0.94, 0.92, 0.904, and 0.87 densities. The isotherms obtained in the gas sorption experimental device display a sorption behavior similar to that of glassy polymers. We apply the dual model to semicrystalline polymers assuming that Henry's sites are related to the amorphous phase, which decreases when the crystallinity percentage increases, whereas the surface of the crystalline phase acts as a Langmuir site with higher gas‐polymer affinity than glassy polymers. The good concordance of the calculated k_D values, using the Flory‐Huggins theory of polymer diluent mixtures, with the experimental results suggest that Henry's gas sorption fulfills this theory and, therefore, it may be a suitable way to estimate polymer‐gas enthalpic interactions. Particularly, the variation of k_D with the crystallinity fraction is exponential and the proportionality of the total sorption with the amorphous content seems only apparent. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1798–1807, 2007     

Adsorption of vapors and gases by carbon sorbents with homogeneous microporous structures in the initial region of the adsorption isotherm

Examination of the adsorption isotherm equation derived from the theory of volume filling of micropores (TVFM) showed that the region of low pressures of experimental isotherms is not always described by this equation. A method for describing the initial region of adsorption isotherms using the parameters of TVFM was proposed.     

Application of water-activated carbon isotherm models to water adsorption isotherms of single-walled carbon nanotubes

The objective of this study is to understand the interactions of water with novel nanocarbons by implementing semiempirical models that were developed to interpret adsorption isotherms of water in common carbonaceous adsorbents. Water adsorption isotherms were gravimetrically determined on several single-walled carbon nanotube (SWNT) and activated carbon samples. Each isotherm was fitted to the Dubinin-Serpinsky (DS) equation, the Dubinin-Astakov equation, the cooperative multimolecular sorption theory, and the Do and Do equations. The applicability of these models was evaluated by high correlation coefficients and the significance of fitting parameters, especially those that delineate the concentration of hydrophilic functional groups, micropore volume, and the size of water clusters. Samples were also characterized by spectroscopic and adsorption techniques, and properties complementary to those quantified by the fitting parameters were extracted from the data collected. The comparison of fitting parameters with sample characterization results was used as the methodology for selecting the most informative and the best-fitting model. We conclude that the Do equation, as modified by Marban et al., is the most suitable semiempirical equation for predicting from experimental isotherms alone the size of molecular clusters that facilitate adsorption in SWNTs, deconvoluting the experimental isotherms into two subisotherms: adsorption onto hydrophilic groups and filling of micropores, and quantifying the concentration of hydrophilic functional groups, as well as determining the micropore volume explored by water. With the exception of the DS equation, the application of other water isotherm models to SWNTs is not computationally tractable. The findings from this research should aid studies of water adsorption in SWNTs by molecular simulation, which remains the most popular tool for understanding the microscopic behavior of water in nanocarbons.     

A density functional theory for Lennard-Jones fluids in cylindrical pores and its applications to adsorption of nitrogen on MCM-41 materials

A density functional theory (DFT) constructed from the modified fundamental-measure theory and the modified Benedict-Webb-Rubin equation of state is presented. The Helmholtz free energy functional due to attractive interaction is expressed as a functional of attractive weighted-density in which the weight function is a mean-field-like type. An obvious advantage of the present theory is that it reproduces accurate bulk properties such as chemical potential, bulk pressure, vapor-liquid interfacial tension, and so forth when compared with molecular simulations and experiments with the same set of molecular parameters. Capabilities of the present DFT are demonstrated by its applicability to adsorption of argon and nitrogen on, respectively, a model cylindrical pore and mesoporous MCM-41 materials. Comparison of the theoretical results of argon in the model cylindrical pore with those from the newly published molecular simulations indicates that the present DFT predicts accurate average densities in the pore, slightly overestimates the pore pressure, and correctly describes the effect of the fluid-pore wall interaction on average densities and pressures in the pore. Application to adsorption of nitrogen on MCM-41 at 77.4 K shows that the present DFT predicts density profiles and adsorption isotherms in good agreement with those from molecular simulations and experiments. In contrast, the hysteresis loop of adsorption calculated from the mean-field theory shifts toward the low pressure region because a low bulk saturated pressure is produced from the mean-field equation of state. The present DFT offers a good way to describe the adsorption isotherms of porous materials as a function of temperature and pressure.     

A simple adsorption model for ionic surfactants

In the past, few theoretical attempts have been made to describe quantitatively the adsorption of ionic surfactants at liquid interfaces. Well-known adsorption isotherms due to Frumkin or Hill–de Boer cannot respond to the specific electrostatic and geometric properties of the surfactant molecules. Our approach is based on a combination of the Gouy–Chapman theory with a modified Frumkin isotherm. The modification implies that the system is free to choose an optimal head group area and an optimal arrangement of the surfactant molecules in the interface as a function of bulk concentration. Interaction energies between neighbouring adsorbed surfactant molecules and between surfactant and water molecules are taken into consideration. The minimum of the Gibbs free energy of the system is equivalent to a minimal interfacial tension. Thus, the thermodynamically stable isotherm can be obtained as the lower envelope of the family of σ versus ln c isotherms resulting from different choices of the model parameters, including the area per molecule. According to the Gibbs equation, the Γ versus ln c adsorption isotherm is obtained as the derivative of this envelope. By variation of the model parameters, the envelope of the calculated adsorption isotherms can be fitted to experimental data of the interfacial tension versus bulk concentration. A computer program is used to calculate the σ versus c and the Γ versus ln c curves as well as to fit the parameters. Received: 28 October 1999/Accepted: 8 February 2000