Dynamics of small molecules in dense polymers subject to thermal motion

Molecular simulation of the effect of temperature and architecture on polyethylene barrier properties

We present a multiscale approach for calculating the low-concentration solubility, diffusivity, and selectivity of small molecules through polymer matrixes. The proposed modeling scheme consists of two main stages; first, thoroughly equilibrated and representative poly(ethylene) (PE) atomistic melt configurations were obtained through the application of a Monte Carlo (MC) scheme based on advanced chain-connectivity altering moves (linear architectures) or the combination of localized MC moves followed by molecular dynamics. In the second phase, transition-state theory (TST), as proposed by Gusev and Suter [Gusev, A. A.; Suter, U. W. J. Chem. Phys. 1993, 99, 2228], was invoked in a coarser level of description to calculate the barrier properties of the studied macromolecules to small gas molecules at infinite dilution. The multiscale methodology was successfully applied on PE melts characterized by various molecular weights (MW) (from C78 up to C1000) and polydispersity indices at a wide range of temperature conditions. The effect of molecular architecture on the barrier properties was examined through the comparison between linear and short-chain branched structures bearing the same total number of carbon atoms. Simulation results were found to be in very good agreement with available experimental data. Additionally, the new scheme has been further validated by comparing the qualitative behavior of solubility, diffusivity, and selectivity with previously reported trends in the literature based on both experimental and simulation studies. The present study concludes that density plays a dominant role that determines the behavior of the polymer as a barrier material, especially in terms of diffusivity. Additionally, it is evidenced that short-chain branching has a small effect on the barrier properties of PE when the comparison is performed on purely amorphous samples. The hierarchical method presented here not only is faster when compared against conventional molecular dynamics simulations, but in some cases, like the vicinity of the glass transition temperature or for long polymer chain melts, it opens the way to the calculation of the barrier properties at realistic simulation times.

Gas diffusion in glasses via a probabilistic molecular dynamics

A probabilistic protocol which makes possible the calculation of the diffusivity of light gases in amorphous materials from limited Monte Carlo and molecular dynamics data is presented. Diffusion coefficients are calculated for helium and methane in polystyrene, and for helium, neon, and methane in three pairs of polysulfone isomers. Results include diffusion coefficients as small as 10(-9) cm2/s and are in good agreement with results obtained from traditional molecular dynamics and with available experimental data.

A generalized direct-particle-deletion scheme for the calculation of chemical potential and solubilities of small- and medium-sized molecules in amorphous polymers

The development, validation, and first applications of a generalized version of an inverse Widom method are described. It permits the calculation of solubility coefficients for molecules as large as, e.g., benzene in all polymers for which reasonable forcefield parameters exist. Predicting the solubility is a key to the knowledge-based design of materials utilized to solve permeability related problems. For long time, particle insertion methods, such as the Widom method, were the only way to predict solubilities from molecular models, but they, in most cases, only worked well for rather small penetrants (e.g., H2, O2, N2). Therefore, a few years ago, a new particle deletion algorithm "DPD" was introduced by Boulougouris, Economou, and Theodorou to overcome this problem in principle. The related computer code was, however, only applicable to special, relatively simple model systems. As application examples for the generalized version described here, solubility calculations for nitrogen, oxygen, and benzene in poly(dimethyl siloxane) are presented.

Molecular dynamics of polystyrene solutions in microwave fields

Equilibrium and nonequilibrium molecular dynamics simulation techniques were used to assess the influence of an applied microwave field on the dynamics of methylamine-methanol and methylamine-dimethylformamide (DMF) solutions bound within atactic polystyrene over a range of polymer densities from 35 to 96 wt % polymer. Atomistically detailed systems were studied, ranging from 3000 to 10 644 particles, using previously established potential models. Structural and dynamical properties were determined in the canonical (NVT) ensemble at 298 K. The simulated DMF self-diffusion coefficients in polystyrene solutions were compared with the zero-field experimental results established with pulsed-gradient spin-echo NMR spectrometry. A simulated external microwave field, with a rms electric field intensity of 0.1 VA, was applied to these systems and the simulated dynamical results over field frequencies up to 10(4) GHz were compared with the zero-field values. Simulated evidence of athermal effects on the diffusive characteristics of these mixtures is reported.

Morphogen gradient formation in a complex environment: an anomalous diffusion model

Current models of morphogen-induced patterning assume that morphogens undergo normal, or Fickian, diffusion, although the validity of this assumption has never been examined. Here we argue that the interaction of morphogens with the complex extracellular surrounding may lead to anomalous diffusion. We present a phenomenological model that captures this interaction, and derive the properties of the morphogen profile under conditions of anomalous (non-Fickian) diffusion. In this context we consider the continuous time random walk formalism and extend its application to account for degradation of morphogen particles. We show that within the anomalous diffusion model, morphogen profiles are fundamentally distinct from the corresponding Fickian profiles. Differences were found in several key aspects, including the role of degradation in determining the profile, the rate by which it spreads in time and its long-term behavior. We analyze the effect of an abrupt change in the extracellular environment on the concentration profiles. Furthermore, we discuss the robustness of the morphogen distribution to fluctuations in morphogen production rate, and describe a feedback mechanism that can buffer such fluctuations. Our study also provides rigorous criteria to distinguish experimentally between Fickian and anomalous modes of morphogen transport.

EVEBAT: a fast strategy for the examination of the empty space in polymer matrices

A very efficient strategy to evaluate the "free volume" and the "unoccupied volume" in polymer matrices is proposed. The method, which has been denoted EVEBAT (Empty Volume Evaluation Based on Atom-Types), emerges from the concept of atom type, which is the basis of a force field. Accordingly, the numerical treatment is minimized by analyzing the atom types used to describe the polymer matrix, which are always a few, rather than all the explicit, atoms involved in it. The performance of the new method has been pointed out by a detailed comparison with the classical insertion algorithms.