Atomistically Modeling the Chemical Potential of Small Molecules in Dense Polymer Microstructures. 1. Method

Modelling Sorption and Transport of Gases in Polymeric Membranes across Different Scales: A Review

Professor Giulio C. Sarti has provided outstanding contributions to the modelling of fluid sorption and transport in polymeric materials, with a special eye on industrial applications such as membrane separation, due to his Chemical Engineering background. He was the co-creator of innovative theories such as the Non-Equilibrium Theory for Glassy Polymers (NET-GP), a flexible tool to estimate the solubility of pure and mixed fluids in a wide range of polymers, and of the Standard Transport Model (STM) for estimating membrane permeability and selectivity. In this review, inspired by his rigorous and original approach to representing membrane fundamentals, we provide an overview of the most significant and up-to-date modeling tools available to estimate the main properties governing polymeric membranes in fluid separation, namely solubility and diffusivity. The paper is not meant to be comprehensive, but it focuses on those contributions that are most relevant or that show the potential to be relevant in the future. We do not restrict our view to the field of macroscopic modelling, which was the main playground of professor Sarti, but also devote our attention to Molecular and Multiscale Hierarchical Modeling. This work proposes a critical evaluation of the different approaches considered, along with their limitations and potentiality.

Molecular View into the Cyclodextrin Cavity: Structure and Hydration

We find, through atomistic molecular dynamics simulation of native cyclodextrins (CDs) in water, that although the outer surface of a CD appears like a truncated cone, the inner cavity resembles a conical hourglass because of the inward protrusion of the glycosidic oxygens. Furthermore, the conformations of the constituent α-glucose molecules are found to differ significantly from a free monomeric α-glucose molecule. This is the first computational study that maps the conformational change to the preferential hydrogen bond donating capacity of one of the secondary hydroxyl groups of CD, in consensus with an NMR experiment. We have developed a simple and novel geometry-based technique to identify water molecules occupying the nonspherical CD cavity, and the computed water occupancies are in close agreement with the experimental and density functional theory studies. Our analysis reveals that a water molecule in CD cavity loses out about two hydrogen bonds and remains energetically frustrated but possesses higher orientational degree of freedom compared to bulk water. In the context of CD-drug complexation, these imply a nonclassical, that is, enthalpically driven hydrophobic association of a drug in CD cavity.

Molecular Modeling Investigations of Sorption and Diffusion of Small Molecules in Glassy Polymers

With a wide range of applications, from energy and environmental engineering, such as in gas separations and water purification, to biomedical engineering and packaging, glassy polymeric materials remain in the core of novel membrane and state-of the art barrier technologies. This review focuses on molecular simulation methodologies implemented for the study of sorption and diffusion of small molecules in dense glassy polymeric systems. Basic concepts are introduced and systematic methods for the generation of realistic polymer configurations are briefly presented. Challenges related to the long length and time scale phenomena that govern the permeation process in the glassy polymer matrix are described and molecular simulation approaches developed to address the multiscale problem at hand are discussed.

Role of solvation dynamics and local ordering of water in inducing conformational transitions in poly(N-isopropylacrylamide) oligomers through the LCST

Conformational transitions in thermo-sensitive polymers are critical in determining their functional properties. The atomistic origin of polymer collapse at the lower critical solution temperature (LCST) remains a fundamental and challenging problem in polymer science. Here, molecular dynamics simulations are used to establish the role of solvation dynamics and local ordering of water in inducing conformational transitions in isotactic-rich poly(N-isopropylacrylamide) (PNIPAM) oligomers when the temperature is changed through the LCST. Simulated atomic trajectories are used to identify stable conformations of the water-molecule network in the vicinity of polymer segments, as a function of the polymer chain length. The dynamics of the conformational evolution of the polymer chain within its surrounding water molecules is evaluated using various structural and dynamical correlation functions. Around the polymer, water forms cage-like structures with hydrogen bonds. Such structures form at temperatures both below and above the LCST. The structures formed at temperatures above LCST, however, are significantly different from those formed below LCST. Short oligomers consisting of 3, 5, and 10 monomer units (3-, 5-, and 10-mer), are characterized by significantly higher hydration level (water per monomer ~ 16). Increasing the temperature from 278 to 310 K does not perturb the structure of water around the short oligomers. In the case of 3-, 5-, and 10-mer, a distinct coil-to-globule transition was not observed when the temperature was raised from 278 to 310 K. For a PNIPAM polymer chain consisting of 30 monomeric units (30-mer), however, there exist significantly different conformations corresponding to two distinct temperature regimes. Below LCST, the water molecules in the first hydration layer (~12) around hydrophilic groups arrange themselves in a specific ordered manner by forming a hydrogen-bonded network with the polymer, resulting in a solvated polymer acting as hydrophilic. Above LCST, this arrangement of water is no longer stable, and the hydrophobic interactions become dominant, which contributes to the collapse of the polymer. Thus, this study provides atomic-scale insights into the role of solvation dynamics in inducing coil-to-globule phase transitions through the LCST for thermo-sensitive polymers like PNIPAM.

Predicting water sorption and volume swelling in dense polymer systems via computer simulation

Atomistic model structures of amorphous polyamide 6 (PA-6) and of an adhesive system consisting of the diglycidyl ether of bisphenol A (DGEBA) as epoxy resin and isophorone diamine (IPD) as a curing agent are generated. For the adhesive, we use a new approach for the generation of the cross-linked polymer networks. It takes into account the chemical reaction kinetics of the curing reaction and, therefore, results in more realistic network structures. On the basis of the corresponding model structures, the equilibrium water content and the swelling ratio of amorphous PA-6 and of the DGEBA+IPD networks are calculated via computer simulation for different thermodynamic conditions. A hybrid method is used combining the molecular dynamics technique with an accelerated test particle insertion method. The results are in reasonable agreement with experiments and, in the case of the PA-6 system, with results obtained via other computer simulation methods.

Grand canonical ensemble molecular dynamics simulation of water solubility in polyamide-6,6

Grand canonical ensemble molecular dynamics simulation is employed to calculate the solubility of water in polyamide-6,6. It is shown that performing two separate simulations, one in the polymeric phase and one in the gaseous phase, is sufficient to find the phase coexistence point. In this method, the chemical potential of water in the polymer phase is expanded as a first-order Taylor series in terms of pressure. Knowing the chemical potential of water in the polymer phase in terms of pressure, another simulation for water in the gaseous phase, in the grand canonical ensemble, is done in which the target chemical potential is set in terms of pressure in the gas phase. The phase coexistence point can easily be calculated from the results of these two independent simulations. Our calculated sorption isotherms and solubility coefficients of water in polyamide-6,6, over a wide range of temperatures and pressures, agree with experimental data.

A comparative study of the properties of polar and nonpolar solvent/solute/polystyrene solutions in microwave fields via molecular dynamics

The influence of an applied microwave field on the dynamics of methylamine-dichloromethane (DCM) mixtures bound within atactic polystyrene (a-PS) over a range of polymer densities from 30 to 94 wt % polymer was examined using atomistic molecular dynamics simulations. This study is an extension of previous studies on methylamine transport in relatively polar polystyrene solutions of methanol and dimethylformamide [M. J. Purdue et al., J. Chem. Phys. 124, 204904 (2006)]. A direct comparison is made across the three types of polystyrene solutions. Consideration is given to both solvent and reagent transport within the polymer solutions under zero-field conditions and in an external electromagnetic field in the canonical ensemble (NVT) at 298.0 K. Various frequencies up to 10(4) GHz and a rms electric field intensity of 0.1 VA were applied. The simulation studies were validated by comparison of the simulated zero-field self-diffusion coefficients of DCM in a-PS with those obtained using pulsed-gradient spin-echo NMR spectrometry. Athermal effects of microwave fields on solute transport behavior within polymer solutions are discussed.

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.