Molecular dynamics simulations of the chain dynamics in monodisperse oligomer melts and of the oligomer tracer diffusion in an entangled polymer matrix

Diffusion Coefficients of Variable-Size Amphiphilic Additives in a Glass-Forming Polyethylene Matrix

Diffusion of additives in polymers is an important issue in the plastics industry since migratory-type molecules are widely used to tune the properties of polymeric composites. Predicting the diffusional behavior of new additives can minimize the need for repetitive experiments. This work presents molecular dynamics simulations at the microsecond time scale and uses the MARTINI force field to estimate self-diffusion coefficients, D, of six monounsaturated amides and their analogs carboxylic acids in polyethylene matrices (PE, MW = 5600 Da). The results are strongly influenced by the glass-forming properties of the PE matrix, which we characterize by three distinct temperatures. The metastability region (T < 325 K), the glass transition temperature (Tg = 256-260 K), and the end of the transition (T ≅ 200 K). Self-diffusion mechanisms are inferred from the results of the dependence of D on the molecular mass of the additive, observing a Rouse-like behavior at high temperatures and deviations from it within the metastability region of the matrix. Interestingly, D values are nonsensitive to the nature of the considered polar head for additives of similar size. The temperature-dependent behavior of D follows, at fixed additive size, a linear Arrhenius pattern at high temperatures and a super Arrhenius trend at lower temperatures, which is well represented with a power law equation as predicted by the Mode Coupling Theory (MCT). We offer a conceptual explanation for the observed super-Arrhenius behavior. This explanation draws on Truhlar and Kohen's interpretation of the available energies at both the initial and the transition states along the diffusion pathway. The matrix's mobility significantly affects solute self-diffusion, yielding equal activation enthalpies for the Arrhenius region or the same power law parameters for the super-Arrhenius regime. Finally, we establish a one-to-one time-equivalence of the self-diffusion processes between CG and all-atom systems for the largest additives and the PE matrix in the high-temperature regime.

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.

Modeling in food across the scales: towards a universal mass transfer simulator of small molecules in food

Multiscale modeling in food is the cutting-edge strategy to revisit food structure and food composition to meet specific targets such as bioavailability, oral perception, or to evaluate the contamination of food by chemicals. A special implementation of Langevin dynamics is proposed to describe mass transfer in structured food. The concepts of random walks over discrete times and physicochemical interactions are connected via an exact solution of the Fokker–Planck equation across interfaces. The methodology is illustrated on the calculation of effective diffusivities of small solutes in emulsions in relationship with their polydispersity, the volume fraction of dispersed phase d = [0.1, 0.4], the ratio of diffusion coefficients between the two phases, rD = [10−2, 102], and the partition coefficients between the continuous and disperse phases, K = [10−2, + ∞[. Simulated diffusion paths are detailed in 2D emulsions and the effective diffusivities compared with the core–shell model of Kalnin and Kotomin (J Phys A Math Gen 31(35):7227–7234, 1998). The same effects are finally tabulated for 3D emulsions covering the full range of food applications. The methodology is comprehensive enough to enable various extensions such as chemisorption, adsorption in the surfactant layer, local flows, flocculation/creaming.

Boson peak, elasticity, and glass transition temperature in polymer glasses: Effects of the rigidity of chain bending

The excess low-frequency vibrational spectrum, called boson peak, and non-affine elastic response are the most important particularities of glasses. Herein, the vibrational and mechanical properties of polymeric glasses are examined by using coarse-grained molecular dynamics simulations, with particular attention to the effects of the bending rigidity of the polymer chains. As the rigidity increases, the system undergoes a glass transition at a higher temperature (under a constant pressure), which decreases the density of the glass phase. The elastic moduli, which are controlled by the decrease of the density and the increase of the rigidity, show a non-monotonic dependence on the rigidity of the polymer chain that arises from the non-affine component. Moreover, a clear boson peak is observed in the vibrational density of states, which depends on the macroscopic shear modulus G. In particular, the boson peak frequency ω_BP is proportional to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{G}$$\end{document}G. These results provide a positive correlation between the boson peak, shear elasticity, and the glass transition temperature.

A blob model to parameterize polymer hole free volumes and solute diffusion

Solute diffusion in solid polymers has tremendous applications in packaging, reservoir, and biomedical technologies but remains poorly understood. Diffusion of non-entangled linear solutes with chemically identical patterns (blobs) deviates dramatically in polymers in the solid-state (α_lin > 1, Macromolecules 2013, 46, 874) from their behaviors in the molten state (α_lin = 1, Macromolecules, 2007, 40, 3970). This work uses the scale invariance of the diffusivities, D, of linear probes D(N·M_blob + M_anchor,T,T_g) = N^{-α_lin(T,T_g)}D(M_blob + M_anchor,T,T_g) comprising N identical blobs of mass M_blob and possibly one different terminal pattern (anchor of mass M_anchor) to evaluate the amounts of hole-free volume in seven polymers (aliphatic, semi-aromatic and aromatic) over a broad range of temperatures (-70 K ≤T-T_g≤ 160 K). The new parameterization of the concept of hole-free volumes opens the application of the free-volume theory (FVT) developed by Vrentas and Duda to practically any polymer, regardless of the availability of free-volume parameters. The quality of the estimations was tested with various probes including n-alkanes, 1-alcohols, n-alkyl acetates, and n-alkylbenzene. The effects of enthalpic and entropic effects of the blobs and the anchor were analyzed and quantified. Blind validation of the reformulated FVT was tested successfully by predicting from first principles the diffusivities of water and toluene in amorphous polyethylene terephthalate from 4 °C to 180 °C and in various other polymers. The new blob model would open the rational design of additives with controlled diffusivities in thermoplastics.

Structure and dynamics of stereo-regular poly(methyl-methacrylate) melts through atomistic molecular dynamics simulations

Poly(methyl-methacrylate), PMMA, is a disubstituted vinyl polymer whose properties depend significantly on its tacticity. Here we present a detailed study of the structure, conformation, and dynamics of syndiotactic, atactic, and isotactic PMMA melts at various temperatures (580, 550, 520, and 490 K) via all-atom molecular dynamics simulations. The calculated volumetric properties are close to experimental data. The T_g and chain dimensions of PMMA model systems are found to depend strongly on tacticity in agreement with experimental findings. The backbone bonds in trans (t), diads in tt, and inter-diads in t|t torsional states are the most populated for all PMMA stereo-chemistries and their fractions increase with the number of syndiotactic sequences. Also, the effective torsional barrier heights for the backbone, ester side group, and α-methyl group are larger for syndiotactic PMMA compared to the isotactic one. The structure of the PMMA chains is studied by computing the intra- and inter-chain static structure factors, S(q), and the radial pair distribution functions. In the first peak of S(q), both intra- and inter-chain components contribute, whereas the second and third peaks mainly come from inter- and intra-chain parts, respectively. For all PMMA stereo-isomers a clear tendency of ester-methyl groups to aggregate is observed. The local dynamics are studied by analyzing torsional autocorrelation functions for various dihedral angles. A wide spectrum of correlation times and different activation energies are observed for the motions of different parts of PMMA chains. The stereo-chemistry affects the backbone, ester side group, and α-methyl motions, whereas the ester-methyl rotation remains unaffected. The dynamic heterogeneity of the PMMA chains is also studied in detail for the different stereo-chemistries via the temperature dependence of the stretching exponent. Furthermore, the reorientational dynamics at the chain level and translational dynamics for monomer and chain centers-of-mass are analyzed.

Molecular thermodynamics for food science and engineering

We argue that thanks to molecular modeling approaches, many thermodynamic properties required in Food Science and Food Engineering will be calculable within a few hours from first principles in a near future. These new possibilities will enable to bridge via multiscale modeling composition, process and storage effects to reach global optimization, innovative concepts for food or its packaging. An outlook of techniques and a series of examples are given in this perspective. We emphasize solute chemical potentials in polymers, liquids and their mixtures as they cannot be understood and estimated without theory. The presented atomistic and coarse-grained methods offer a natural framework to their conceptualization in polynary systems, entangled or crosslinked homo- or heteropolymers.

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.

A Review of Multiscale Computational Methods in Polymeric Materials

Polymeric materials display distinguished characteristics which stem from the interplay of phenomena at various length and time scales. Further development of polymer systems critically relies on a comprehensive understanding of the fundamentals of their hierarchical structure and behaviors. As such, the inherent multiscale nature of polymer systems is only reflected by a multiscale analysis which accounts for all important mechanisms. Since multiscale modelling is a rapidly growing multidisciplinary field, the emerging possibilities and challenges can be of a truly diverse nature. The present review attempts to provide a rather comprehensive overview of the recent developments in the field of multiscale modelling and simulation of polymeric materials. In order to understand the characteristics of the building blocks of multiscale methods, first a brief review of some significant computational methods at individual length and time scales is provided. These methods cover quantum mechanical scale, atomistic domain (Monte Carlo and molecular dynamics), mesoscopic scale (Brownian dynamics, dissipative particle dynamics, and lattice Boltzmann method), and finally macroscopic realm (finite element and volume methods). Afterwards, different prescriptions to envelope these methods in a multiscale strategy are discussed in details. Sequential, concurrent, and adaptive resolution schemes are presented along with the latest updates and ongoing challenges in research. In sequential methods, various systematic coarse-graining and backmapping approaches are addressed. For the concurrent strategy, we aimed to introduce the fundamentals and significant methods including the handshaking concept, energy-based, and force-based coupling approaches. Although such methods are very popular in metals and carbon nanomaterials, their use in polymeric materials is still limited. We have illustrated their applications in polymer science by several examples hoping for raising attention towards the existing possibilities. The relatively new adaptive resolution schemes are then covered including their advantages and shortcomings. Finally, some novel ideas in order to extend the reaches of atomistic techniques are reviewed. We conclude the review by outlining the existing challenges and possibilities for future research.

The unique role of bond length in the glassy dynamics of colloidal polymers

Bond length is generally not considered as a controllable variable for molecular polymers. Hence, no experimental, simulation or theoretical research, to our knowledge, has examined the influence of bond length on the glassy dynamics of polymers. Recently, a new class of assembling materials called "colloidal polymers" has been synthesized. These colloidal polymers have advantages over molecular polymers in the visibility and flexibility of tuning, for example, the size and shape of the "monomers", the interaction, and the bond length. Dense suspension of colloidal polymers will become a very promising ideal model system for exploring the fundamental problems in the glass transition of chain "molecules". Here, we study the static structure and activated dynamics of hard-sphere colloidal polymers by generalizing the colloidal nonlinear Langevin equation theory to colloidal polymers. Surprisingly, we find that the bond length plays a critical and unique role in many aspects. For instance, the universal relations of the characteristic local lengths and the activated barrier versus the "degree of supercooling", and the structural relaxation versus local vibrational motion are found to be dependent on bond length and independent of chain length and rigidity. We hope that our findings inspire future experimental and simulation research studies on the glassy dynamics of colloidal polymers.

Simulated glass-forming polymer melts: dynamic scattering functions, chain length effects, and mode-coupling theory analysis

We present molecular-dynamics simulations for a fully flexible model of polymer melts with different chain length N ranging from short oligomers (N = 4) to values near the entanglement length (N = 64). For these systems we explore the structural relaxation of the supercooled melt near the critical temperature T c of mode-coupling theory (MCT). Coherent and incoherent scattering functions are analyzed in terms of the idealized MCT. For temperatures T > T c we provide evidence for the space-time factorization property of the β relaxation and for the time-temperature superposition principle (TTSP) of the α relaxation, and we also discuss deviations from these predictions for T ≈ T c. For T larger than the smallest temperature where the TTSP holds we perform a quantitative analysis of the dynamics with the asymptotic MCT predictions for the late β regime. Within MCT a key quantity, in addition to T c, is the exponent parameter λ. For the fully flexible polymer models studied we find that λ is independent of N and has a value (λ = 0.735 ) typical of simple glass-forming liquids. On the other hand, the critical temperature increases with chain length toward an asymptotic value T c (∞) . This increase can be described by T c (∞) - T c(N) ∼ 1/N and may be interpreted in terms of the N dependence of the monomer density ρ, if we assume that the MCT glass transition is ruled by a soft-sphere-like constant coupling parameter Γ c = ρ c T c (-1/4), where ρ c is the monomer density at T c. In addition, we also estimate T c from a Hansen-Verlet-like criterion and MCT calculations based on structural input from the simulation. For our polymer model both the Hansen-Verlet criterion and the MCT calculations suggest T c to decrease with increasing chain length, in contrast to the direct analysis of the simulation data.

Molecular dynamics simulation of melting and crystallization processes of polyethylene clusters confined in armchair single-walled carbon nanotubes

The confined interaction is important to understand the melting and crystallization of polymers within single-wall carbon tube (SWNT). However, it is difficult for us to observe this interaction. In the current work, the structures and behaviors of melting and crystallization for polyethylene (PE) clusters confined in armchair single-walled carbon nanotubes ((n,n)-SWNTs) are investigated and examined based on molecular dynamics (MD) simulations. The nonbonded energies, structures, Lindemman indices, radial density distributions, and diffusion coefficients are used to demonstrate the features of melting phase transition for PE clusters confined in (n,n)-SWNTs. The chain end-to-end distance (R(n)) and chain end-to-end distribution are used to examine the flexibility of the PE chain confined in SWNT. The global orientational order parameter (P2) is employed to reveal the order degree of whole PE polymer. The effect of polymerization degree on melting temperature and the influence of SWNT chirality on structure of PE cluster are examined and discussed. Results demonstrate that within the confined environment of SWNT, PE clusters adopt novel co-axial crystalline layer structure, in which parallel chains of each layer are approximately vertical to tube axis. The disordered-ordered transformation of PE chains in each layer is an important structural feature for crystallization of confined PE clusters. SWNTs have a considerable effect on the structures and stabilities of the confined PE clusters.

Some Insight into Stability of Amorphous Poly(ethylene glycol) Dimethyl Ether Polymers Based on Molecular Dynamics Simulations

Poly(ethylene glycol) dimethyl ether (PEGDME) polymers are widely used as drug solid dispersion reagents. They can cause the amorphization of drugs and improve their solubility, stability, and bioavailability. However, the mechanism about why amorphous PEGDME 2000 polymer is highly stable is unclear so far. Molecular dynamics (MD) simulation is a unique key technique to solve it. In the current work, we systematically investigate structure, aggregate state, and thermodynamic and kinetic behaviors during the phase-transition processes of the PEGDME polymers with different polymerization degree in terms of MD simulations. The melting and glass-transition temperatures of the polymers are in good agreement with experimental values. The amorphous PEGDME2000 exhibits high stability, which is consistent with the recent experiment results and can be ascribed to a combination of two factors, that is, a high thermodynamic driving force for amorphization and a relatively low molecular mobility.

Hydrodynamic and viscoelastic effects in polymer diffusion

We develop a fluctuating hydrodynamics approach to study the impact of the hydrodynamic and viscoelastic interactions on the motion of the center of mass of a polymer as well as on the relaxation of Rouse modes, either in a Θ solvent or in a melt of identical unentangled chains. We show that this method allows us to describe the effect of hydrodynamic interactions beyond the Zimm (for a single chain in a Θ solvent) or the Rouse models (for an unentangled melt). In the latter case, we recover the same important effect of the viscoelastic hydrodynamic interactions on the center-of-mass diffusion, first described in Farago et al (2011 Phys. Rev. Lett. 107 178301).

Simulated glass-forming polymer melts: glass transition temperature and elastic constants of the glassy state

By means of molecular-dynamics simulation we study a flexible and a semiflexible bead-spring model for a polymer melt on cooling through the glass transition. Results for the glass transition temperature T(g) and for the elastic properties of the glassy state are presented. We find that T(g) increases with chain length N and is for all N larger for the semiflexible model. The N dependence of T(g) is compared to experimental results from the literature. Furthermore, we characterize the polymer glass below T(g) via its elastic properties, i.e., via the Lamé coefficients λ and μ. The Lamé coefficients are determined from the fluctuation formalism which allows to split λ and μ into affine (Born term) and nonaffine (fluctuation term) contributions. We find that the fluctuation term represents a substantial correction to the Born term. Since the Born terms for λ and μ are identical, the fluctuation terms are responsible for the different temperature dependence of the Lamé coefficients. While λ decreases linearly on approaching T(g) from below, the shear modulus μ displays a much stronger decrease near T(g). From the present simulation data it is not possible to decide whether μ takes a finite value at T(g), as would be expected from mode-coupling theory, or vanishes continuously, as suggested by recent work from replica theory.