Influence of System Size in Molecular Dynamics Simulations of Gas Permeation in Glassy Polymers

Effect of solvent on the emulsion and morphology of polyfluorene films: all-atom molecular dynamics approach

The morphology of conjugated polymer thin films deposited by the resonant infrared matrix-assisted pulsed laser evaporation (RIR-MAPLE) process is related to the emulsion characteristics. However, a fundamental understanding of how and why the emulsion characteristics control the film properties and device performance is yet unclear. We performed all-atom molecular dynamics simulations of emulsions containing a mixture of polyfluorene (PFO) polymer, various primary solvents, secondary solvent, and water. The emulsion properties were then examined as a function of variable primary solvent and correlated with the morphology of deposited PFO thin films. The examination of the explicit interactions between all components of the emulsion indicated that using a primary solvent with a lower solubility-in-water and a higher non-bonded interaction energy ratio, between the solvent, polymer, and water in the emulsion recipe, produced the best result with smoother and denser films. Additionally, our simulation results are consistent with the AFM experimental results, indicating that interactions driven by trichlorobenzene (TCB) primary solvent within the emulsion are responsible for high-quality, smooth, and continuous thin film surfaces. Overall, this study can support the choice of a suitable primary solvent and provides the computational framework for predictions of new recipes for polymeric emulsion systems.

In Silico Prediction of Food Properties: A Multiscale Perspective

Several open software packages have popularized modeling and simulation strategies at the food product scale. Food processing and key digestion steps can be described in 3D using the principles of continuum mechanics. However, compared to other branches of engineering, the necessary transport, mechanical, chemical, and thermodynamic properties have been insufficiently tabulated and documented. Natural variability, accented by food evolution during processing and deconstruction, requires considering composition and structure-dependent properties. This review presents practical approaches where the premises for modeling and simulation start at a so-called “microscopic” scale where constituents or phase properties are known. The concept of microscopic or ground scale is shown to be very flexible from atoms to cellular structures. Zooming in on spatial details tends to increase the overall cost of simulations and the integration over food regions or time scales. The independence of scales facilitates the reuse of calculations and makes multiscale modeling capable of meeting food manufacturing needs. On one hand, new image-modeling strategies without equations or meshes are emerging. On the other hand, complex notions such as compositional effects, multiphase organization, and non-equilibrium thermodynamics are naturally incorporated in models without linearization or simplifications. Multiscale method’s applicability to hierarchically predict food properties is discussed with comprehensive examples relevant to food science, engineering and packaging. Entropy-driven properties such as transport and sorption are emphasized to illustrate how microscopic details bring new degrees of freedom to explore food-specific concepts such as safety, bioavailability, shelf-life and food formulation. Routes for performing spatial and temporal homogenization with and without chemical details are developed. Creating a community sharing computational codes, force fields, and generic food structures is the next step and should be encouraged. This paper provides a framework for the transfer of results from other fields and the development of methods specific to the food domain.

Polyimides Containing Phosphaphenanthrene Skeleton: Gas-Transport Properties and Molecular Dynamics Simulations

A series of new semifluorinated polyimide (PI) films with phosphaphenanthrene skeleton were prepared by thermal imidization of poly(amic acid)s derived from a diamine monomer: 1,1-bis[2'-trifluoromethyl-4'-(4″-aminophenyl)phenoxy]-1-(6-oxido-6H-dibenz⟨c,e⟩⟨1,2⟩oxaphosphorin-6-yl)ethane on reaction with four structurally different aromatic dianhydrides. The chemical structures of the polymers were established by Fourier transform infrared and 1H NMR spectroscopy techniques. The polymers showed a good combination of thermal and mechanical properties (T d10 up to 416 °C under synthetic air and tensile strength up to 91 MPa), low dielectric constant (2.10-2.55 at 1 MHz), and T g values as high as 261 °C. Gas permeabilities of these films were investigated for four different gases CO2, O2, N2, and CH4. The PI films showed high gas permeability (P CO2 up to 175 and P O2 up to 64 barrer) with high permselectivity (P CO2 /P CH4 up to 51 and P O2 /P N2 up to 7.1), and the values are better than those of many other similar polymers reported earlier. For the O2/N2 gas pair, the PIs (PI A) surpassed the present upper boundary limit drawn by Robeson. A detailed molecular dynamics (MD) simulation study has been conducted to understand better the gas-transport properties. The effect of phosphaphenanthrene skeleton, its spatial arrangement, and size distribution function of the free volume were studied using molecular dynamics (MD) simulation and the results are correlated with the experimental data.

Soluble, optically transparent polyamides with a phosphaphenanthrene skeleton: synthesis, characterization, gas permeation and molecular dynamics simulations

Soluble, optically transparent polyamides with a phosphaphenanthrene skeleton: synthesis, characterization, gas permeation and molecular dynamics simulations.

Computational insights on the role of film thickness on the physical properties of ultrathin polysulfone membranes

A pioneering work to elucidate physical properties of ultrathin membrane films from atomistic point of view in Materials Studio.

A comparison of homopolymer and block copolymer structure in 6FDA-based polyimides

Two homopolyimides and the corresponding block copolyimide, all based on the 4,4'-(hexafluoroisopropylidene)diphthalic dianhydride (6FDA), have been synthesized and fully atomistic models have been studied using molecular dynamics (MD) simulation. The respective diamines were 1,3-phenylenediamine (mPDA) and 2,3,5,6-tetramethyl-1,4-phenylenediamine (durene). These polyimides are potential candidates for gas separation applications. The synthesized polymers were processed as dense flat membranes. The effects of diamine structure were investigated at the molecular level and an attempt to compare the structural features of homo- and block copolyimides was made. Amorphous models were generated using a hybrid pivot Monte Carlo-MD sampling preparation technique. Average model densities were validated against experimental measurements on the dense films. Cohesive energies, Hildebrand solubility parameters, conformational characteristics, intermolecular interactions and available void spaces were analysed for each system. The durene diamine was found to hinder stacking and increase the available space. This is associated with the steric effect of the methyl substituents. In general, 6FDA-mPDA/durene exhibits an intermediate behaviour with respect to its base polyimides. For most of the examined properties, the differences between different size simulated systems were minor with the exception of the free volume distribution.

Sorption and diffusion of carbon dioxide and nitrogen in poly(methyl methacrylate)

Molecular dynamics simulations are performed to determine the solubility and diffusion coefficient of carbon dioxide and nitrogen in poly(methyl methacrylate) (PMMA). The solubilities of CO2 in the polymer are calculated employing our grand canonical ensemble simulation method, fixing the target excess chemical potential of CO2 in the polymer and varying the number of CO2 molecules in the polymer matrix till establishing equilibrium. It is shown that the calculated sorption isotherms of CO2 in PMMA, employing this method well agrees with experiment. Our results on the diffusion coefficients of CO2 and N2 in PMMA are shown to obey a common hopping mechanism. It is shown that the higher solubility of CO2 than that of N2 is a consequence of more attractive interactions between the carbonyl group of polymer and the sorbent. While the residence time of CO2 beside the carbonyl group of polymer is about three times higher than that of N2, the diffusion coefficient of CO2 in PMMA is higher than that of N2. The higher diffusion coefficient of CO2, compared to N2, in PMMA is shown to be due to the higher (≈3 times) swelling of polymer upon CO2 uptake.

Thermal properties of bulk polyimides: insights from computer modeling versus experiment

Due to the great importance for many industrial applications it is crucial from the point of view of theoretical description to reproduce thermal properties of thermoplastic polyimides as accurate as possible in order to establish "chemical structure-physical properties" relationships of new materials. In this paper we employ differential scanning calorimetry, dilatometry, and atomistic molecular dynamics (MD) simulations to explore whether the state-of-the-art computer modeling can serve as a precise tool for probing thermal properties of polyimides with highly polar groups. For this purpose the polyimide R-BAPS based on dianhydride 1,3-bis(3',4-dicarboxyphenoxy)benzene (dianhydride R) and diamine 4,4'-bis(4''-aminophenoxy)biphenyl sulphone) (diamine BAPS) was synthesized and extensively studied. Overall, our findings show that the widely used glass-transition temperature Tg evaluated from MD simulations should be employed with great caution for verification of the polyimide computational models against experimental data: in addition to the well-known impact of the cooling rate on the glass-transition temperature, correct definition of Tg requires cooling that starts from very high temperatures (no less than 800 K for considered polyimides) and accurate evaluation of the appropriate cooling rate, otherwise the errors in the measured values of Tg become undefined. In contrast to the glass-transition temperature, the volumetric coefficient of thermal expansion (CTE) does not depend on the cooling rate in the low-temperature domain (T < Tg) so that comparison of computational and experimental values of CTE provides a much safer way for proper validation of the theoretical model when electrostatic interactions are taken into account explicitly. Remarkably, this conclusion is most likely of generic nature: we show that it also holds for the commercial polyimide EXTEM, another polyimide with a similar chemical structure.