Estimation of mutual diffusion coefficients in polymer/penetrant systems using nonequilibrium molecular dynamics simulations

Predicting diffusion coefficients of chemicals in and through packaging materials

Most of the physicochemical properties in polymers such as activity and partition coefficients, diffusion coefficients, and their activation with temperature are accessible to direct calculations from first principles. Such predictions are particularly relevant for food packaging as they can be used (1) to demonstrate the compliance or safety of numerous polymer materials and of their constitutive substances (e.g. additives, residues…), when they are used: as containers, coatings, sealants, gaskets, printing inks, etc. (2) or to predict the indirect contamination of food by pollutants (e.g. from recycled polymers, storage ambiance…) (3) or to assess the plasticization of materials in contact by food constituents (e.g. fat matter, aroma…). This review article summarizes the classical and last mechanistic descriptions of diffusion in polymers and discusses the reliability of semi-empirical approaches used for compliance testing both in EU and US. It is concluded that simulation of diffusion in or through polymers is not limited to worst-case assumptions but could also be applied to real cases for risk assessment, designing packaging with low leaching risk or to synthesize plastic additives with low diffusion rates.

Slippage and viscosity predictions in carbon micropores and their influence on CO2 and CH4 transport

Non-equilibrium molecular dynamics simulations of pure carbon dioxide and methane and their equimolar mixtures have been carried out with an external driving force imposed on carbon slit pores to investigate gas slippage and Klinkenberg effects. Simulations were conducted to determine the effect of pore size and exposure to an external potential on the velocity profile and slip-stick boundary conditions. The simulations indicate that molecule-wall collisions influence the velocity profile, which deviates significantly from the Navier-Stokes hydrodynamic prediction for micro- and mesopores. Also, the shape of the velocity profile is found to be independent of the applied pressure gradient in micropores. The results indicate that the velocity profile is uniform for pore sizes less than 2 nm (micropores) where the transport is mainly due to molecular streaming or Knudsen diffusion and, to a lesser extent, molecular diffusion. As pore sizes increase to 10 nm, parabolic profiles are observed due to the reduced interaction of gas molecules with the pore walls. A 3D pore network, representative of porous carbon-based materials, has been generated atomistically using the Voronoi tessellation method. Simulations have been carried out to determine the effect of the pore structure and modeled viscosity on permeability and Klinkenberg parameters. The use of the bulk-phase viscosity for estimating the permeability of CO(2) in units of Darcy in a 3D micropore network is not an appropriate assumption as it significantly underestimates the CO(2) permeability. On the other hand, since the transport properties of CH(4) are less influenced by the pore walls compared with CO(2), the use of the bulk-phase CH(4) viscosity estimates are a reasonable assumption.

Brownian dynamics simulation of tracer diffusion in a cross-linked network

Brownian dynamics simulation is employed to study the self-diffusion of tracer particles in a cross-linked gel network based on a coarse-grained bead-spring lattice model. Several effects are investigated including the network porosity, flexibility, degree of cross linking, and electrostatic interaction. For uncharged systems, the tracer long-time diffusivity is found to decrease with deceasing porosity and flexibility, but with increasing degree of cross linking. For charged systems, the diffusion is further hindered by the electrostatic interaction, regardless of whether the tracer particle and the network are oppositely or similarly charged. However, there exists a difference in the hindrance mechanism between the two cases. For the former, a substantial decrease in the diffusivity can occur at high network porosities due to electrostatic entrapment.

Molecular dynamics simulations of transport and separation of carbon dioxide–alkane mixtures in carbon nanopores

The configurational-bias Monte Carlo method, which is used for efficient generation of molecular models of n-alkane chains, is combined for the first time with the dual control-volume grand-canonical molecular-dynamics simulation, which has been developed for studying transport of molecules in pores under an external potential gradient, to investigate transport and separation of binary mixtures of n-alkanes, as well as mixtures of CO2 and n-alkanes, in carbon nanopores. The effect of various factors, such as the temperature of the system, the composition of the mixture, and the pore size, on the separation of the mixtures is investigated. We also report the preliminary results of an experimental study of transport and separation of some of the same mixtures in a carbon molecular-sieve membrane with comparable pore sizes. The results indicate that, for the mixtures considered in this paper, even in very small carbon nanopores the energetic effects still play a dominant role in the transport and separation properties of the mixtures, whereas in a real membrane they are dominated by the membrane's morphological characteristics. As a result, for the mixtures considered, a single pore may be a grossly inadequate model of a real membrane, and hence one must resort to three-dimensional molecular pore network models of the membrane.